Decreased Serum microRNA-21, microRNA-25, microRNA-146a, and microRNA-181a in Autoimmune Diabetes: Potential Biomarkers for Diagnosis and Possible Involvement in Pathogenesis

Objective Previous studies have revealed dysregulated circulating microRNAs (miRNAs) in patients with type 1 diabetes (T1D). Here, we explored the serum levels of miR-21, miR-25, miR-146a, and miR-181a in patients with autoimmune diabetes (T1D and latent autoimmune diabetes of adults (LADA)) compared with type 2 diabetes (T2D) and nondiabetic individuals. Design, patients, and measurements The serum levels of miR-21, miR-25, miR-146a, and miR-181a in patients with T1D (n = 29), LADA (n = 16), and T2D (n = 31) and in nondiabetic individuals (n = 19) were determined by quantitative real-time polymerase chain reaction, and receiver-operating characteristic (ROC) curves were evaluated to determine the discriminatory performances of these four miRNAs. Furthermore, target genes and pathways potentially modulated by these four miRNAs were predicted by bioinformatics analysis to investigate the possible functions of these miRNAs in autoimmune diabetes. Subsequently, multiple logistic regression analysis was performed to identify independent predictors for autoimmune diabetes, and a nomogram was established. Results miR-21, miR-25, miR-146a, and miR-181a were significantly downregulated in the serum of patients with autoimmune diabetes compared with those in T2D patients and nondiabetic individuals (p < 0.001). The areas under the ROC curves of these four miRNAs were greater than 0.80 (p < 0.001). Bioinformatics analysis suggested that miR-21, miR-25, miR-146a, and miR-181a regulated multiple genes in pathways associated with immunity, inflammatory responses, hyperglycemia, and metabolism, which are involved in the pathogenesis of autoimmune diabetes. Multiple logistic regression analysis identified miR-25 (odds ratio (OR): 0.001, p < 0.05), miR-146a (OR: 0.136, p < 0.05), and fasting C-peptide levels (OR: 0.064, p < 0.05) as independent predictors of autoimmune diabetes. Conclusions miR-25 and miR-146a may serve as potential circulating biomarkers and provide insights into the pathogenesis of autoimmune diabetes.


Introduction
Type 1 diabetes (T1D) is a chronic progressive autoimmune disease characterized by T-cell-mediated pancreatic β-cell destruction, ultimately leading to absolute insulin deficiency and exogenous insulin-dependent hyperglycemia [1]. Latent autoimmune diabetes of adults (LADA), another specific subtype of autoimmune diabetes with less intensive pancreatic β-cell destruction compared with T1D at its initial stage, shares several common immunological and genetic features with T1D, such as the existence of islet autoantibodies, the important role of cellular immunity in disease pathogenesis, and shared risk loci [2]. Despite much effort, autoimmune diabetes still faces enormous challenges, including delayed diagnosis owing to the overlap of clinical characteristics with type 2 diabetes (T2D), insufficient reliable biomarkers, lack of appropriate targets, and lack of strategies to reverse or delay islet β-cell failure, given the unclear mechanism. erefore, biomarkers with the potential to facilitate diagnosis and provide insights into the pathogenesis of autoimmune diabetes are urgently needed.
MicroRNAs (miRNAs) are endogenous, evolutionarily conserved, small, double-stranded, noncoding RNAs measuring approximately 19-24 nucleotides in length; these molecules act as negative posttranscriptional modulators of mRNA expression [3]. Accumulating evidence indicates that multiple circulating miRNAs are dysregulated in T1D and participate in the pathogenesis of this disease, thus serving as potential noninvasive biomarkers and therapeutic targets for T1D. For example, miR-21, miR-25, miR-146a, and miR-181a [4] are associated with the regulation of immune responses, β-cell apoptosis and proliferation, and insulin biosynthesis and secretion, which are all involved in the pathogenesis of autoimmune diabetes. However, the roles of these miRNAs remain to be elucidated. miR-146a has essential regulatory roles in T-cell biology; however, most studies have used peripheral blood mononuclear cells (PBMCs) or T cells as samples, necessitating extended storage and processing before measurements. Moreover, few studies have evaluated serum miR-146a expression patterns in samples from patients with T1D. miR-21, which was previously reported to have critical roles in T-cell activation, is upregulated in serum from patients with T1D and speculated to have dual effects on pancreatic β cells [5,6]; therefore, miR-21 may dynamically change during different stages of T1D. miR-25 has previously been reported to have prognostic value with regard to the functions of residual β cells and glycemic control several months later in patients with T1D [7], necessitating additional studies to explore its association with residual β-cell function and glycemic control. Moreover, miR-181a has been reported to play a prominent role in T-cell activation, which is important in the pathogenesis of T1D [8]. Furthermore, few studies have investigated alterations in these four miRNAs in the circulation of patients with LADA, another important subtype of autoimmune diabetes.
Accordingly, in this study, we examined alterations in the levels of miR-21, miR-25, miR-146a, and miR-181a in the serum of patients with T1D and LADA to identify potential circulating biomarkers and gain insights into the pathogenesis of autoimmune diabetes.

Study Populations.
Using protocols and consent procedures approved by the ethics committee of the Peking Union Medical College Hospital, 95 individuals attending the Peking Union Medical College Hospital from January 2014 to May 2016 were recruited to the current study, including patients with T1D (n � 29), LADA (n � 16), and T2D (n � 31) and nondiabetic individuals (n � 19). All participants provided written informed consent for the collection of blood samples. Blood samples were collected and centrifuged for 10 minutes at 3000 rpm to separate serum, which was subsequently refrigerated at − 80°C before further analysis. Diabetes was diagnosed based on the World Health Organization criteria [9]. e classification of T1D was made according to the characteristics as follows: dependence on prompt insulin treatment at a rapid onset of diabetes, fasting C-peptide or 2-h postprandial C-peptide lower than 0.8 ng/mL, and positivity for islet autoantibodies. e diagnostic criteria of LADA were as follows: 30-70 years old at diagnosis; positivity of one of the islet autoantibodies, including glutamic acid decarboxylase antibody (GADA), islet cell antibody (ICA), and protein tyrosine phosphatase antibody (IA-2A); and independence from insulin therapy for at least 6 months after diagnosis [10]. Exclusion criteria included suspected other types of diabetes and possibility of organ dysfunction, infection, and pregnancy. All patients with T1D or LADA were treated with insulin, and some patients also received oral antidiabetic agents, including metformin (T1D: 7/29, LADA: 5/16) and α-glycosidase inhibitors (T1D: 5/29, LADA: 6/16). T2D patients also received oral antidiabetic agents, including metformin (22/31), α-glycosidase inhibitors (24/31), dipeptidyl peptidase inhibitors (10/31), and sulfonylureas (13/31); few of these patients were treated with insulin (7/31). Demographic data were collected, including sex, age, age at diagnosis, diabetes duration, weight, height, and body mass index (BMI).

miRNA Target Gene Prediction and Pathway Analysis.
Target gene prediction for these four miRNAs was performed using four web-based prediction tools, including MiRWALK [12], miRTarBase [13], miRDB [14], and Tar-getScan [15]. To control the false-positive rate, target genes were selected on the basis of at least three adopted prediction tools. Subsequently, functional enrichment analysis of target genes for these four miRNAs was performed with pathway annotations from the Kyoto Encyclopedia of Genes and Genomes (KEGG) database [16]. Significantly targeted pathways enriched for target genes were identified based on Fisher exact tests (p < 0.01).

Statistical Analysis.
All analyses were implemented using SPSS Statistics software (Version 25.0; SPSS, Chicago, IL, USA), R (Version 3.5.0), and GraphPad Prism 6.0 (http://www. graphpad.com). Two-sided tests were used, and statistical significance was established at a P value of 0.05. Continuous data were presented as means ± standard deviations. Normal distributions were evaluated by Kolmogorow-Smirnow and Shapiro-Wilk tests. Differences between groups were tested by nonparametric Mann-Whitney or Kruskal-Wallis tests. Receiver-operating characteristic (ROC) curves were established, and the areas under the ROC curves (AUC-ROCs) were calculated to evaluate the discriminatory power of the four miRNAs to distinguish T1D and LADA from T2D and nondiabetic individuals. Spearman correlations were conducted to determine the associations among these four miRNAs and clinical variables. Multiple logistic regression analysis was performed to identify independent predictors for autoimmune diabetes, and a nomogram was constructed based on the regression coefficients (β) from the multivariate logistic model.

Clinical Characteristics of Study Populations.
In total, 95 participants were studied. As indicated in Table 1, patients with T1D were younger and had a younger age of diabetes onset compared with patients with LADA and T2D. No significant differences were found among the three diabetic groups with regard to diabetes duration and BMI. Patients with autoimmune diabetes (T1D and LADA) suffered from less favorable glycemic control and residual β-cell function, with significantly elevated islet autoantibody titers compared with patients with T2D.

Serum miR-21, miR-25, miR-146a, and miR-181a
Were Downregulated in Autoimmune Diabetes. Differences in the serum levels of these four miRNAs were not significant between T2D patients and nondiabetic individuals. In contrast, the levels of these four miRNAs were downregulated in patients with T1D and LADA compared with those in patients with T2D and nondiabetic individuals, although the difference for patients with LADA was not significant. Additionally, decreases in the levels of these four miRNAs in individuals with LADA were weaker than those in patients with T1D, consistent with the notion that LADA is an intermediate form between T1D and T2D ( Figure 1). Next, ROC curves were performed to evaluate the discriminatory performance of the four miRNAs in distinguishing T1D and LADA from T2D and nondiabetic individuals. As shown in Figure 2, AUCs of the four miRNAs were all greater than 0.8 (miR-21: 0.834, 95% confidence interval (CI): 0.751-0.918; miR-25: 0.872, 95% CI: 0.799-0.945; miR-146a: 0.909, 95% CI: 0.850-0.968; miR-181a: 0.844, 95% CI: 0.764-0.924), indicating the favorable performance of these miRNAs for distinguishing T1D and LADA from T2D and nondiabetic individuals and their potential to serve as circulating biomarkers for autoimmune diabetes.

Serum Levels of miR-21, miR-25, miR-146a, and miR-181a Did Not Correlate with Glycemic Control, Residual β-Cell Function, and Islet Autoantibodies in Autoimmune Diabetes.
Spearman correlation analysis was conducted to explore the correlations of serum levels of miR-21, miR-25, miR-146a, and miR-181a with clinical parameters in autoimmune diabetes (Table 2). Notably, these miRNAs showed no correlations with HbA1c, FCP, or 2hCP, indicating no associations between serum levels of these miRNAs and glycemic control or residual β-cell function. Additionally, we failed to identify associations between serum levels of these miRNAs and titers or numbers of positive islet autoantibodies. Nonetheless, serum levels of miR-25 were positively correlated with age (correlation coefficient: 0.576, p < 0.01), BMI (correlation coefficient: 0.296, p < 0.05), and age of diabetes onset (correlation coefficient: 0.483, p < 0.01). -21, miR-25, miR-146a, and miR-181a. To investigate the possible functions of these four miRNAs in autoimmune diabetes, bioinformatics analysis was implemented to predict target genes and pathways potentially affected by these miRNAs. In total, 1,481 genes were predicted by at least three prediction tools for these four miRNAs. Subsequently, a multitude of genes were found to be targeted by multiple miRNAs; after excluding duplicates, 1,342 unique genes were predicted to be targeted by these four miRNAs (Figure 3(a)). Forty-five, 26, 25, and 47 pathways were identified as enriched with genes targeted by miR-21, miR-25, miR-146a, and miR-181a, respectively, after excluding malignancy-associated pathways. Some pathways overlapped, yielding 51 unique pathways.

Discussion
Reliable biomarkers are required to achieve more accurate diagnosis and provide insights into the regulation of the pathogenesis of autoimmune diabetes. In this study, we demonstrated significant downregulation of serum miR-21, miR-25, miR-146a, and miR-181a in T1D and LADA patients compared with T2D patients and nondiabetic individuals. Bioinformatics analysis revealed that these four miRNAs targeted 1,342 genes, which could impact 51 KEGG pathways, most of which were related to immunity, inflammation responses, diabetes, and metabolism and were therefore involved in autoimmune diabetes. miR-25, miR-146a, and FCP were independent predictors of autoimmune diabetes, suggesting their potential to serve as circulating biomarkers of autoimmune diabetes, which may be of immense value in clinical practice.
Several previous studies reported increased miR-21 in the plasma/serum of patients with T1D compared with that in nondiabetic individuals, whereas one study found no differences in miR-21 expression in the plasma between patients with T1D and patients with T2D [7,[17][18][19]. In contrast, we found a significant decrease in miR-21 levels in the serum of patients with T1D compared with that in nondiabetic individuals and patients with T2D in this study. is difference could be related to differences in blood samples or clinical features of the participants included in the different studies (e.g., sex, age, and disease duration). miR-21 has previously  been identified as a critical regulator of immune responses, participating in the regulation of anti-inflammatory responses [20], T lymphocyte activation [21], and regulatory T-cell (Treg) development [22]. Additionally, a growing number of studies have indicated the dual regulatory effects of miR-21 in pancreatic β-cell apoptosis, including antiapoptotic effects via suppression of tumor-suppressive programmed cell death [4], which promotes cell death [5] and pro-apoptosis effect via inhibition of the antiapoptotic BCL2 [6]. Additionally, miR-21 is associated with impaired glucose-stimulated insulin secretion [23].
Changes in circulating miR-25 levels reported in previous studies have not been conclusive. One publication showed downregulation [18], whereas others showed upregulation of miR-25 in the serum/plasma of patients with T1D compared with that in nondiabetic individuals [7,11,24]. Notably, serum miR-25 measured 1 month after T1D onset was positively correlated with stimulated C-peptide levels and negatively correlated with HbA1c measured 3 months after T1D onset, suggesting its potential prognostic role in predicting β-cell functions and glycemic control in the progression of T1D [7]. However, we failed to find an association between serum miR-25 and glycemic control or residual β-cell function, probably resulting from the different disease course between the study populations different analytical platforms applied and the absence of follow-up data after diabetes onset in the current study. Several previous publications have shown that miR-25 has multiple effects on the regulation of apoptosis via multiple pathways in various types of cancer [25,26]. However, previous studies have rarely investigated the associations between miR-25 and pancreatic β-cell apoptosis in T1D.
Similarly, alterations in miR-146a in different studies have also been contradictory. miR-146a was previously reported to be reduced in PBMCs of patients with T1D [11,27], whereas miR-146a was increased in T lymphocytes of patients with T1D in another publication [28]. miR-146a is essential for the immunosuppressive function of Tregs and is thus recognized as a canonical regulator of immune responses [29]. Loss of miR-146a could generate spontaneous autoimmunity. Interestingly, a previous study showed that decreased miR-146a in PBMCs was correlated with high GADA titers and GADA positivity, supporting that miR-146a may be associated with the ongoing islet autoimmunity in T1D [11,27]; however, we did not replicate this observation. Interestingly, despite their important roles as biomarkers of autoimmunity, islet autoantibodies do not have major roles in pancreatic β-cell destruction.
In previous studies, miR-181a was found to be upregulated in the serum of patients with T1D [7,30]. However, in this study, we found the opposite results. Moreover, loss of miR-181a increases the reactivity of peripheral T cells against self-antigens, indicating its essential role in modulation of T-cell activity [8]. miR-181a also has anti-inflammatory effects given its ability to inhibit increases in TNF, interleukin (IL)-1β, and IL-6 in macrophages treated with lipopolysaccharides [31].
In the current study, we found markedly decreased levels of serum miR-21, miR-25, miR-146a, and miR-181a in autoimmune diabetes, inconsistent with some previous studies. is discrepancy could be attributed to differences in blood samples collected, detection methods of miRNAs, and clinical features of the participants (e.g., sex, age, and disease duration) as well as a lack of replicates in previous investigations. Given the diverse regulatory roles of miRNAs in multiple pathways, it is possible that the alterations and regulatory functions of miRNAs may vary depending on the disease stage. erefore, changes in the expression levels of miRNAs may not necessarily be persistent during the course of disease progression.
Bioinformatics analysis suggested that miR-21, miR-25, miR-146a, and miR-181a regulate multiple genes in pathways associated with immunity and inflammatory responses related to the pathogenesis of autoimmune diabetes, such as TNF, Toll-like receptor, MAPK, Fc epsilon RI, thyroid hormone, osteoclast differentiation, T-cell receptor, NODlike receptor, and AGE-RAGE signaling pathways. Additionally, multiple signaling pathways associated with metabolism and diabetes/hyperglycemia were enriched with target genes of these four miRNAs, such as sphingolipid, cGMP-PKG, FoxO, AMPK, Rap1, and insulin signaling pathways. Furthermore, multiple signaling pathways have been reported to participate in the development of diabetic complications. e AGE-RAGE signaling pathway was recently reported to promote diabetes-mediated vascular calcification via increasing oxidative stress [32]. e axon guidance signaling pathway was previously shown to be involved in the development of diabetic nephropathy and retinopathy [33,34]. Additionally, the Rap1 signaling pathway is inhibited by hyperglycemia, resulting in tubular injury in individuals with diabetic nephropathy [35].
Multiple logistic regression analysis identified miR-25, miR-146a, and FCP as independent predictors for autoimmune diabetes, and a nomogram constructed based on these data exhibited favorable performance in estimation of autoimmune diabetes, suggesting the potential applications of these targets as circulating biomarkers for autoimmune diabetes. However, because most patients with autoimmune diabetes included in this study had a disease duration of several years, we did not directly demonstrate the predictive roles of miR-25 and miR-146a for autoimmune diabetes. erefore, further studies comparing individuals at risk for autoimmune diabetes with healthy controls are necessary to validate the predictive value of miR-25 and miR-146a.
ere were several limitations to this study. First, the results of the current study differed from those of previous publications, and we speculate that differences in study populations and methods used to measure miRNAs may explain these inconsistencies, which could lead to the failure of direct comparison with the results of the previous study. Second, the restricted generalizability of the results and decreased statistical power need to be considered, given the relatively small sample size. ird, to the best of our knowledge, for the first time, we found similar decreases in miR-146a and miR-181a levels in patients with LADA, suggesting the possibility that LADA and T1D share similar pathogenic mechanisms. However, further studies International Journal of Endocrinology are required to validate these results. Lastly, our study was a cross-sectional observational study; thus, we were able to identify only associations, not causal relationships of the four miRNAs with autoimmune diabetes.

Conclusions
In summary, our findings suggested that serum miR-21, miR-25, miR-146a, and miR-181a were significantly downregulated in autoimmune diabetes and may be associated with regulation of the pathogenesis of autoimmune diabetes from various aspects, including autoimmunity, inflammation responses, pancreatic β-cell apoptosis, and impaired insulin secretion. Furthermore, serum miR-25 and miR-146a may serve as independent predictors for estimation of autoimmune diabetes. ese findings could have relevance for clinical practice and suggest that these miRNAs could serve as novel circulating biomarkers of the pathogenesis of autoimmune diabetes. Further investigations are necessary to illuminate their functional roles in the pathogenesis of autoimmune diabetes.

Data Availability
e SPSS Statistics data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest
e authors declare that there are no conflicts of interest regarding the publication of this paper.