Abstract:Background Type 1 diabetes mellitus (T1DM) is caused by immune cell-mediated β-cell dysfunction. In recent decades, N6-methyladenosine (m6A) has attracted widespread attention in the scientific research field because it plays vital roles in the pathogenesis of immunity-related diseases, including autoimmune diseases. However, neither the m6A modification profile nor the potential role it plays in T1DM pathogenesis has been investigated to date. Materials and Methods An m6A mRNA epitranscriptomic microarray analysis was performed to analyze m6A regulator expression patterns and m6A methylation patterns in immune cells of T1DM patients (n=6) and healthy individuals (n=6). A bioinformatics analysis was subsequently performed to explore the potential biological functions and signaling pathways underlying T1DM pathogenesis. Furthermore, mRNA expression and m6A methylation levels were subsequently verified by qRT–PCR and methylated RNA immunoprecipitation–qPCR (MeRIP–qPCR), respectively, in the T1DM and healthy groups (n=6 per group). Results Among the multiple m6A regulators, METTL3 and IGF2BP2 had significantly downregulated expression, and YTHDC1 and HNRNPA2B1 had significantly upregulated expression in the T1DM group relative to the healthy group. The microarray analysis revealed 4247 differentially methylated transcripts, including 932 hypermethylated and 3315 hypomethylated transcripts, and 4264 differentially expressed transcripts, including 1818 upregulated transcripts and 2446 downregulated transcripts in the T1DM group relative to the healthy group. An association analysis between methylation and gene expression demonstrated that the expression of 590 hypermethylated transcripts was upregulated, and that of 1890 hypomethylated transcripts was downregulated. Pearson correlation analysis showed significant correlations between the expression levels of differentially expressed m6A regulators and the methylation levels of differentially methylated transcripts and significant correlations between the expression levels of differentially expressed m6A regulators and that of differentially expressed transcripts. Moreover, Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses demonstrated that differentially methylated transcripts were involved in pathways related to immunity, including some closely associated with T1DM. Conclusions Our study presents m6A regulator expression patterns and m6A methylation patterns of immune cells in T1DM, showing that the m6A mark and m6A regulators are promising targets for T1DM diagnosis and treatment.

Genome-Wide Identification of N6-Methyladenosine Associated SNPs As Potential Functional Variants for Type 1 Diabetes

Integrative Analysis Predicts M 6A-Snps in Type 2 Diabetes Pathogenesis

Pgpm_a_334346 1369..1378

The m6A methylation profiles of immune cells in type 1 diabetes mellitus

SNPS in the Exons of Toll-Like Receptors Are Associated with Susceptibility to Type 1 Diabetes in Chinese Population

A regulatory variant rs9379874 in T1D risk region 6p22.2 affects BTN3A1 expression regulating T cell function

Polymorphisms in NLRP1 Gene Are Associated with Type 1 Diabetes

Single-Nucleotide Polymorphisms Inside Microrna Target Sites Influence The Susceptibility To Type 2 Diabetes

Micro-ribonucleic Acid-Binding Site Variants of Type 2 Diabetes Candidate Loci Predispose to Gestational Diabetes Mellitus in Chinese Han Women.

Risk Variants Disrupting Enhancers of T(H)1 and T-REG Cells in Type 1 Diabetes

Identification of genetic variants in m6A modification genes associated with pancreatic cancer risk in the Chinese population

Novel Association of SNP Rs2297828 in PRDM16 Gene with Predisposition to Type 2 Diabetes.

Identification of genetic variants in m 6 A modification genes associated with pancreatic cancer risk in the Chinese population

A functional variant rs912304 for late-onset T1D risk contributes to islet dysfunction by regulating proinflammatory cytokine-responsive gene STXBP6 expression

Effects of Genetic Variants of Nuclear Receptor Y on the Risk of Type 2 Diabetes Mellitus.

Genetic Variants in Peroxisome Proliferator-Activated Receptor-γ and Retinoid X Receptor-α Gene and Type 2 Diabetes Risk: A Case-Control Study of a Chinese Han Population

Genetic Variants in AC092159.2 and Risk of Gestational Diabetes Mellitus in a Chinese Population.

Integrative Genomic Analysis Predicts Regulatory Role of N6-Methyladenosine-Associated SNPs for Adiposity

U‐shaped Association Between Serum IGF2BP3 and T2DM: A Cross‐sectional Study in Chinese Population

In silico genome-wide identification of m6A-associated SNPs as potential functional variants for periodontitis

Identification and functional study of GATA4 gene regulatory variants in type 2 diabetes mellitus