- Research article
- Open Access
- Published:
Identifying common and specific microRNAs expressed in peripheral blood mononuclear cell of type 1, type 2, and gestational diabetes mellitus patients
BMC Research Notes volume 6, Article number: 491 (2013)
Abstract
Background
Regardless the regulatory function of microRNAs (miRNA), their differential expression pattern has been used to define miRNA signatures and to disclose disease biomarkers. To address the question of whether patients presenting the different types of diabetes mellitus could be distinguished on the basis of their miRNA and mRNA expression profiling, we obtained peripheral blood mononuclear cell (PBMC) RNAs from 7 type 1 (T1D), 7 type 2 (T2D), and 6 gestational diabetes (GDM) patients, which were hybridized to Agilent miRNA and mRNA microarrays. Data quantification and quality control were obtained using the Feature Extraction software, and data distribution was normalized using quantile function implemented in the Aroma light package. Differentially expressed miRNAs/mRNAs were identified using Rank products, comparing T1DxGDM, T2DxGDM and T1DxT2D. Hierarchical clustering was performed using the average linkage criterion with Pearson uncentered distance as metrics.
Results
The use of the same microarrays platform permitted the identification of sets of shared or specific miRNAs/mRNA interaction for each type of diabetes. Nine miRNAs (hsa-miR-126, hsa-miR-1307, hsa-miR-142-3p, hsa-miR-142-5p, hsa-miR-144, hsa-miR-199a-5p, hsa-miR-27a, hsa-miR-29b, and hsa-miR-342-3p) were shared among T1D, T2D and GDM, and additional specific miRNAs were identified for T1D (20 miRNAs), T2D (14) and GDM (19) patients. ROC curves allowed the identification of specific and relevant (greater AUC values) miRNAs for each type of diabetes, including: i) hsa-miR-1274a, hsa-miR-1274b and hsa-let-7f for T1D; ii) hsa-miR-222, hsa-miR-30e and hsa-miR-140-3p for T2D, and iii) hsa-miR-181a and hsa-miR-1268 for GDM. Many of these miRNAs targeted mRNAs associated with diabetes pathogenesis.
Conclusions
These results indicate that PBMC can be used as reporter cells to characterize the miRNA expression profiling disclosed by the different diabetes mellitus manifestations. Shared miRNAs may characterize diabetes as a metabolic and inflammatory disorder, whereas specific miRNAs may represent biological markers for each type of diabetes, deserving further attention.
Background
MicroRNAs (miRNAs) are small noncoding single-stranded RNAs of approximately 22 nucleotides, well-conserved among different species. The precursors of miRNAs are transcribed into the nucleus and several steps are required until mature miRNAs can be exported to the cytoplasm and incorporated into the RNA-induced silencing complex (RISC). This process can affect gene expression at posttranslational level, leading to mRNA degradation, translational repression or both. As a consequence, the same mRNA expression can be modulated by one or several miRNAs, while one miRNA can regulate several genes. Much attention has been devoted to the role of miRNAs acting on the posttranscriptional gene control and it has been suggested that they may be used as biomarkers of disorders of varied etiologies [1].
MiRNAs present a tissue-specific expression pattern [2] and are involved in the regulation of several important cellular functions, such as cell cycle regulation, apoptosis, differentiation and maintenance of the immune system cell repertoire [3]. Although many factors are implicated on gene regulation, at least one third of the transcription control may be attributed to the action of miRNAs [4]. The probability of interaction of miRNAs with mRNAs may be evaluated by several algorithms in several databases, including TargetScan, Pictar, miRanda and others [5–7]. The advent of microarray platforms for miRNAs and mRNAs has provided a powerful tool to unveil new mRNA/miRNA interactions. In this context, considering the availability of miRNAs in tissues, plasma or in peripheral blood, aberrant expression of miRNAs has been used as biomarkers for neoplastic and non-neoplastic disorders [8–10].
Diabetes is one of the most prevalent chronic diseases, affecting 6.4% of the world’s adult population [11], and the disease may be classified into three principal types: type 1 (T1D), type 2 (T2D), and gestational diabetes mellitus (GDM). T1D is an autoimmune disease characterized by an inflammatory response, which leads to a progressive destruction of pancreatic beta cells, resulting in deregulation of glucose metabolism and insulin deficiency. T2D is usually a non-insulin-dependent diabetes, generally associated with obesity, caused by deficient insulin secretion by pancreatic beta-cell islets or deficient insulin action in peripheral tissues [12]. GDM has been defined as any degree of glucose intolerance with onset or first recognition during pregnancy, and, in most of cases, after birth delivery occur the stabilization of glucose metabolism and of the oral glucose tolerance test.
Several studies have been conducted to evaluate specifically expressed miRNAs in the diverse types of diabetes, particularly associated with the regulation of insulin production and secretion [13–18], differentiation of human pre-adipocytes [19], and association with T1D pathogenesis [20]. Other studies investigated the involvement of differentially expressed miRNAs in T1D [3] and T2D [21, 22]; however, there is no information regarding the comparisons of the major types of diabetes using the same microarray platforms. In the present study, we evaluated the mRNA and miRNA profiles of patients presenting T1D, T2D and GDM, using the same microarray platforms, intending to unveil shared and privately expressed miRNAs.
Results
The differential mRNA/miRNA expression analysis revealed similarities and dissimilarities among the different types of diabetes mellitus. Both mRNA and miRNA analyses were performed comparing two groups; i.e., T1D versus GDM, T2G versus GDM and T1D versus T2D. Statistical analysis of mRNAs by rank products comparing groups of patients yielded 523 differentially expressed transcripts when comparing T1D versus GDM, 328 transcripts for T2G versus GDM, and 477 for T1D versus T2D (P ≤ 0.001). Rank products analysis for miRNAs yielded 54 (T1D versus GDM), 28 (T2G versus GDM) and 31 (T1D versus T2D) differentially expressed miRNAs. As seen in Figure 1, the transcript profiles of mRNA and miRNA of patients clearly separated them into distinct clusters.
Hierarchical clustering of mRNA (upper dendrograms) and microRNAs (lower dendrograms). Clustering analyses refer to the comparisons of the transcript profiles (mRNA and miRNA) between T1D versus GDM (1A and 1D), between T2D versus GDM (1B and 1E), and between T1D versus T2D (1C and 1F). As observed, the mRNA and miRNA profiles were distinct for each type of diabetes.
Overall, the interaction of mRNAs with miRNAs disclosed several predicted interactions, which have been previously described in databanks [23–26]; however, most of these interactions have not been reported in association with diabetes. The comparison between T1D and GDM (523 mRNA and 54 miRNA differentially expressed) revealed 31 predicted interactions, including 21 distinct mRNAs and 13 distinct miRNAs. Among these 21 mRNAs, 8 were downregulated (ABCA13, HTRA3, SPTB, SLC6A8, ANKRD22, OXTR, GYPA, and OR51E2) and 13 were upregulated (ACTL6B, C1orf87, FEZ1, IQGAP3, DDX3Y, EIF1AY, ZFY, COL13A1, HRH4, UTY, DLG5, MS4A2, and SLC11A2). The analysis of miRNAs revealed that 5 out of 13 were down regulated (hsa-miR-636, hsa-miR-939, hsa-miR-720, hsa-miR-595, and hsa-miR-623) and 8 were up regulated (hsa-miR-338-3p, hsa-miR-342-3p, hsa-miR-30b, has-miR-30c, hsa-miR-27a, hsa-miR-27b, hsa-miR-374a, and hsa-miR-92a). These results are shown in Table 1.
The comparison T2D versus GDM (328 mRNAs and 28 miRNAs differentially expressed) yielded 42 predicted interactions, encompassing 23 transcripts and 17 miRNAs (Table 2). Among the 23 differentially expressed mRNAs, 16 were downregulated (MMP1, RASGEF1B, SLC6A8, IL1A, EGR3, OXTR, GYPA, ORS1E2, MMP8, SEZ6L2, HTRA3, C8orf4, NR4A1, EREG, LALBA, and CTTN) and 7 upregulated (EIF1AY, ZFY, LPAL2, DDX3Y, SETD5, BAGALNT3, and MLYCD). Among the 17 differentially expressed miRNAs, 6 out of 17 were downregulated (hsa-miR-451, hsa-miR-199a-3p, hsa-miR-595, hsa-miR-1268, hsa-miR-181d, and hsa-miR-486-5p) and 11 were upregulated (hsa-miR-342-3p, hsa-miR-30b, hsa-miR-144, hsa-miR-140-3p, hsa-miR-30e, hsa-miR-142-5p, hsa-miR-378, hsa-miR-181a, hsa-miR-101, hsa-miR-142-3p, and hsa-miR-324-5p) (Table 2).
Finally, the comparison T1D versus T2D (477 mRNA and 31 miRNA differentially expressed) produced 80 predicted interactions, encompassing 42 mRNAs and 23 miRNAs. Among the differentially expressed mRNAs, 12 were downregulated (BATF2, SPG20, ZSWIM3, ANKRD22, IFNG, CA1, ADAMTSL4, SETD5, TRIM36, C1orf25, C9orf40, and PARS2) and 30 were upregulated (C1orf87, PANX2, PDE4D, PITPNM2, RASGEF1B, FBLN2, GPR125, SNF1LK, ABCD3, CCL20, MPHOSPH6, KRT2, STK35, IL1A, CXCL3, SLC11A2, CRABP1, DDX3Y, SLC16A3, IRS2, SNAI1, ZNF507, ID1, SOD2, TBKBP1, TSPYL5, C1orf96, DFNB31, PER1, and EREG). Regarding miRNAs, 10 were downregulated (hsa-miR-342-3p, hsa-miR-720, hsa-miR-144, hsa-miR-140-3p, hsa-miR-30e, hsa-miR-1260, hsa-miR-1308, hsa-miR-142-5p, hsa-miR-29b, and hsa-miR-142-3p) and 13 were upregulated (hsa-miR-27a, hsa-miR-21, hsa-miR-130a, hsa-miR-150, hsa-miR-223, hsa-miR-451, hsa-let-7a, hsa-let-7e, hsa-let-7f, hsa-let-7 g, hsa-miR-20b, hsa-miR-199a-3p, and hsa-miR-103), yielding 80 interactions shown in Table 3.
As shown in Tables 1, 2 and 3, several types of mRNA/miRNA expected interactions were observed; i. e., induction of a miRNA and repression of an mRNA or vice versa, and also simultaneous induction or repression of both miRNAs and mRNAs. Several predicted interactions were observed, many of them encompassing genes involved in the pathogenesis of diabetes, for instance miR-30e targeting IL1A and IRS2, miR-181a targeting IL1A, miR-223 targeting SLC11A2, miR-29b targeting ID1, miR-21 targeting CCL20, and many others. Figure 2 shows mRNA/miRNA gene networks generated from data observed in Tables 1, 2 and 3 (upper networks show all predicted interactions, and lower networks show only the negative correlations; i. e., increased miRNA versus decreased mRNA or vice-versa).
Networks between miRNAs and mRNAs. The relationship between miRNAs and mRNAs were evaluated by constructing networks using the Cytoscape software. (A) Upper networks show all interactions described in Tables 1, 2 and 3, and (B) lower networks show only the negative correlations, i.e., increased miRNA versus decreased mRNA or vice-versa. Red circles represent miRNAs and the grey ones represent mRNA.
To identify shared and specific miRNAs involved in each type of diabetes, several Venn diagrams were constructed taking into account only the differentially expressed miRNAs (Figure 3). After performing the multiple comparisons between T1D versus T2D, T1D versus GDM and T2D versus GDM, we observed that 9 miRNAs were shared by T1D, T2D and GDM, i.e., hsa-miR-126, hsa-miR-1307, hsa-miR-142-3p, hsa-miR-142-5p, hsa-miR-144, hsa-miR-199a-5p, hsa-miR-27a, hsa-miR-29b, and hsa-miR-342-3p. These miRNAs appear in all intersections of the four Venn diagrams shown in Figure 3. Besides the shared miRNAs, specific miRNAs were identified for each type of diabetes, including: i) T1D: hsa-let-7f, hsa-let-7 g, hsa-miR-103, hsa-miR-1260, hsa-miR-1274a, hsa-miR-1274b, hsa-miR-130a, hsa-miR-150, hsa-miR-20b, hsa-miR-21 and hsa-miR-720, ii) T2D: hsa-miR-140-3p, hsa-miR-199a-3p, hsa-miR-222, hsa-miR-30e and hsa-miR-451, and iii) GDM: hsa-miR-101, hsa-miR-1180, hsa-miR-1268, hsa-miR-181a, hsa-miR-181d, hsa-miR-26a, hsa-miR-29a, hsa-miR-29c, hsa-miR-30b and hsa-miR-595.
Venn diagrams showing common and specific microRNAs for the three types of diabetes. The central intersection of the upper diagram shows the nine shared miRNAs among T1D, T2D and GDM, the upper right intersection shows the 5 miRNAs specific for T2D, the upper left intersection shows the 11 miRNAs specific for T1D, and the middle lower intersection shows the 10 miRNAs specific for GDM patients. Lower Venn diagrams identify specific miRNAs for each type of diabetes (bold letters), as well as the shared ones.
The receiver operating characteristic curve (ROC curve) is a diagram to assess sensitivity in function of the specificity. It is usually used in medicine to determine the accuracy of a test to discriminate disease cases from normal cases [27], and to establish a cut-off value for some diagnostic tests. The ROC curve can also be used to compare more than one diagnostic performance from different laboratories [28]. In the present study, the ROC analysis was used to search for potential biomarkers for each type of diabetes. Several ROC curves were constructed for specific miRNAs, considered to be of biological relevance those exhibiting values of area under curve (AUC) greater than 0.8, i.e., high sensibility and high specificity. Accordingly, the hsa-miR-1274a, hsa-miR-1274b and hsa-let-7f presented better results for T1D; hsa-miR-222, hsa-miR-30e and hsa-miR-140-3p for T2D; and hsa-miR-181a and hsa-miR-1268 for GDM patients (Figure 4).
Identification of most relevant specific miRNAs for each diabetes type. The values of the area under the curve (AUC) were estimated for all specific miRNAs obtained after the multiple comparisons among the three types of diabetes as shown in Figure 3. MiRNAs exhibiting high AUC values are highlighted within blue rectangles.
Discussion
The analysis of miRNA expression profile displayed by the three types of diabetes showed distinct hybridization signatures according to the type of diabetes, in accordance with the results reported by us in a previous study evaluating the mRNA profiling of these patients [29]. The major focuses on the analysis of mRNA/miRNA signatures encompassing T1D, T2D and GDM patients included: i) pinpoint miRNAs shared among the three types of diabetes; ii) select miRNAs specific for each type of diabetes; iii) disclose putative mRNAs targeted by the differentially expressed miRNAs; and iv) identify non-described miRNAs associated with each type of diabetes. All selected miRNAs were detected in at least 3 out of 6 algorithm analyses used in the present study.
MiRNAs play regulatory roles in many biological processes associated with diabetes, including adipocyte differentiation, metabolic integration, insulin resistance and appetite regulation [30]. The role of miRNAs in diabetes has been associated with several pathogenic features, including: i) miR-410, miR-200a and miR-130a regulate secretion of insulin in response to stimulatory levels of glucose, and overexpression of miR-410 enhances the levels of glucose-stimulated insulin secretion [18]; ii) miRNA-143 is upregulated during differentiation of human pre-adipocytes [19]; iii) miR-30d is upregulated in pancreatic beta-cells and collaborates for increasing insulin gene expression [17]; iv) miR-9 acts in the fine-tuning of glucose metabolism [14]; v) miR-375 is involved on the control of insulin gene expression and secretion [13, 15]; vi) overexpression of miR-29 represses insulin-stimulated glucose uptake and may result in insulin resistance [31].
Considering that 9 miRNAs were shared among the three types of diabetes, including hsa-miR-126, hsa-miR-144, hsa-miR-27a, hsa-miR-29b, hsa-miR-1307, hsa-miR-142-3p, hsa-miR-142-5p, hsa-miR-199a-5p, and hsa-miR-342-3p, we may suppose that these miRNAs are associated with diabetes per se.
MiR-126 has been associated with immune response of adipose tissue and macrophage infiltration into the adipose tissue [32]. It is expressed in endothelial cells and contributes to the endothelial homeostasis and vascular integrity [33]. In diabetic patients, alteration in miR-126 expression may inhibit endothelial progenitor cell proliferation, migration, and may induce apoptosis [34]. High glucose concentration is associated with reduced miR-126 content in endothelial apoptotic bodies [21], and a positive association between miR-126 expression and myocardial infarction has been observed in diabetic patients, probably due to the role of miR-126 on platelet function [35].
MiR-144, miR-27a and miR-29b present a more relevant role on diabetes. MiR-144 targets IRS1 (insulin receptor substrate-1), a gene highly involved in insulin signaling pathway, and upregulation of this miRNA exhibits a linear relationship with the glycemic status in T2D patients [22]. Thus, the control of miR-144 expression may be a potential therapeutic target for T2D patients, deserving further studies.
MiR-27a together with miR-150, miR-192, miR-320a, and miR-375 regulate several biological events related to the pathogenesis of diabetes [22, 36–38]. MiR-27a has been associated with hyperglycemia and metabolic syndrome in T2D patients: i) upregulation of this miRNA has been observed in hyperglycemic rats exhibiting T2D [39]; ii) its expression is associated with the fasting glucose level, suggesting its potential role in the early-phase hyperglycemia [40]; iii) considering that miR-27a has a potential angiogenic function, its downregulation in diabetes patients should reduce the angiogenic potential of endothelial progenitor cells in diabetes [34]. In the present study, miR-27a was highly expressed in T1D followed by T2D and GDM. Although there are no studies evaluating the role of miR-27a in T1D, this miRNA may be involved in shared mechanisms for hyperglycemia control in the major types of diabetes.
Several studies report that the miR-29 family, particularly, miR-29b has a role in diabetes: i) in the T2D rat model (Goto-Kakizaki), overexpression of miR-29 family represses insulin-stimulate glucose uptake, facilitating insulin resistance [31]; ii) miR-29 family members (miR-29a, miR-29b, and miR-29c) are expressed in mouse pancreatic beta-cells, and their expression increases with the age of prediabetic NOD mice [41], contributing to insulin resistance in animal model of diabetes T2D [31]; iii) miR-29 overexpression downregulates glucose-induced insulin secretion in human islet cells [41]; iv) miR-29b is highly expressed in neurons of the retinal ganglion cells of diabetic rats, suggesting a role of miR-29b in the pathogenesis of diabetic retinopathy [42]; v) besides neurons, miR-29b has been associated with progression of renal fibrosis, irrespective of the etiology, including diabetic nephropathy [43]. In concert, these findings indicate that the miR-29b and other members of this family may contribute to diabetes pathogenesis and diabetes microvascular complications.
No association with diabetes has been described for miR-1307, miR-142-3p, miR-142-5p, miR-199a-5p and miR-342-3p; however, miR-142-5p and miR-142-3p have been reported as negative regulators of CD4 T cells, and increased expression of miR-142-5p and miR-142-3p is reported in systemic lupus erythematosus, contributing to reverse T cell hyperactivity, inhibiting antibody production, reducing autoimmune activity and IgG production [44]. Deregulation of some of these miRNAs has been associated with carcinogenesis: i) miR-142-3p with acute myeloid leukemia [45], T-cell leukemogenesis [46], esophageal squamous cell carcionoma [47], and hepatocellular carcinoma [48]; ii) miR-142-5p with T cell chronic lymphocytic leukemia [49] and in pediatric brain tumors [50], and iii) miR-199a-5p with hepatocellular carcinoma [51, 52], ovarian cancer [53], and oral squamous cell carcinoma [54]. In addition, miR-342-3p has been associated with cellular proliferation [55], and has been suggested as a potential biomarker for multiple myeloma [56] and prion disease [57]. Interestingly, miR-342-3p was associated with obesity in mice [58], and the possible link between this miRNA and obesity is glucagon, a potential target for miR-342-3p. Furthermore, miR-342-3p is associated with immune response, and it is upregulated in activated B cells [59]. Therefore, it is possible that the role of miR-342-3p in obesity and its association with immune response may link this miRNA with diabetes.
To discuss specific miRNAs for each type of diabetes, we selected only those presenting high values of area under curve (AUC = 1.0) of the ROC curves (Figure 4). Considering that: i) the number of differentially expressed miRNAs is very high, ii) microRNA/mRNA interactions are similarly high, and iii) the number of patients evaluated in each group of diabetes of the present study is low, miRNAs exhibiting high AUC may represent those with greater potential to act as biomarker for larger studies. For T1D, let-7f, miR-1274a and miR-1274b fit this assumption. Overall, let-7 family members act as tumor suppressors [60, 61], but the Lin28/let-7 pathway regulates glucose metabolism in different organs, and the treatment with anti-let7 has been suggested as a potential therapy for T2D [62]. The miR-1274a and miR-1274b have never been studied in diabetes, and some authors believe that they represent a fragment of tRNA, without known function [63].
For T2D patients, hsa-miR-140-3p, hsa-miR-222 and hsa-miR-30e were induced and exhibited AUC =1.0. The miR-222 overexpression has also been associated with the development of GDM [64], and is overexpressed in T2D animal models [39] exhibiting obesity [65], and in pancreatic cancer [66]. Mir-140-3p targets the SIRT1 (sirtuin 1) gene, and reduction of miR-140-3p levels has been associated with calorie restriction [67]. SIRT1 is increased in response to a long-term calorie restriction in a tissue-specific manner [68, 69]. Although no association has been reported with diabetes, miR-30e is involved with adipose tissue dysfunction and obesity, as observed for T2D [70]. In addition, miR-30 family members are responsible for expression of mesenchymal proteins during pancreatic fibrosis [71], and miR-30d is upregulated by glucose and increased insulin gene expression, suggesting that miR-30d may be a negative regulator for insulin gene expression [17].
Regarding GDM, two miRNA exhibited AUC =1.0: i) hsa-miR-1268 that was induced and ii) hsa-miR-181a that was repressed in relation to T1D and T2D. The few studies evaluating miR-1268 are inconclusive about its function. On the other hand, miR-181a has been considered to be a link between adipose tissue dysfunction and the development of obesity-associated disorders including T2D [72]. Considering that the overexpression of miR-181a induces hepatic insulin resistance in T2D patients, this miRNA has been suggested as potential diagnostic marker for this type of diabetes [72]; however, miR-181a is also overexpressed in children with new onset T1D [73]. Since the induction of miR-181a has been previously associated with T1D and T2D, and since this miRNA was repressed in GDM when compared to T1D and T2D, miR-181a deserves further studies in GDM in terms of pathophysiology and disease marker.
Conclusions
The present study revealed shared miRNAs among the major types of diabetes, which were also observed in the MiRDisease database [74], including hsa-miR-29b, miR-142-3p and hsa-miR-142-5p. Many of these shared miRNAs have been associated with metabolic pathways, immunological processes and tumorigenesis. Some of these miRNAs were also shared by other immune-mediated or chronic inflammatory diseases, including psoriasis (hsa-miR-142-3p) and lupus nephritis (hsa-miR-150 and hsa-miR-142-5p).
This is the first miRNA study involving the three major groups of diabetes using the same microarray platform, showing that: i) the miRNA signature of each diabetes subset is distinct, ii) shared miRNAs in the major types of diabetes are associated with biological functions associated with disease pathogenesis, iii) privately expressed miRNAs may be used for further and larger studies as disease biomarkers.
Methods
Patients
Peripheral blood mononuclear cells (PBMCs) were obtained from 20 adult patients, i.e., 7 presenting T1D (aged 18–27 years), 7 with T2D (aged 41–61 years) and 6 with GDM (aged 29 to 39 years), followed-up at the Outpatient Clinics of the Division of Endocrinology, Faculty of Medicine of Ribeirao Preto, University of Sao Paulo, Brazil. The demographic, laboratory and treatment features of T1D, T2D, and GDM patients are shown in Table 4. The local Ethics Committee approved the protocol of the study (# 9153/2008) and informed written consents were obtained for all patients. Exclusion criteria included recent episodes of ketoacidosis, active nephropathy, proliferative retinopathy, diabetic foot, high HDL (high-density lipoprotein) levels and diagnosed cardiovascular diseases.
Blood collection and RNA samples
Twenty mL of peripheral blood cells were collected and used for isolation of PBMCs by gradient density (Ficoll-Hypaque- Sigma, St. Louis, MO). Total RNA was obtained using the Trizol reagent (Invitrogen, Carlsbad, CA), according to the manufacturer’s instructions. RNA concentrations and ratios were checked with the NanoDrop ND-1000 spectrophotometer (NanoDrop Products, Wilmington, DE) and RNA quality was assessed using the 2100 Bioanalyzer (Agilent Technologies, Santa Clara, CA). All RNA samples exhibited high integrity numbers (RIN) ≥ 9.0.
Oligo microarrays
Total mRNA was hybridized to one-color Agilent whole human genome 4x44K oligo microarray platforms (Agilent Technologies, Santa Clara, CA), containing ~41.000 transcripts, according to manufacturer’s instructions. MiRNA microarray was performed using Agilent human miRNA platform (v3) 8x15K (Agilent Technologies). After hybridization procedures, the slides were scanned using a Microarray Scanner with Surescan High-Resolution Technology (Agilent Technologies). A complete file providing the hybridization profiles of all samples as well as the quantitative data and experimental conditions are available in the ArrayExpress database [75] through the following account numbers: T1D (E-MEXP-3348), T2D (E-MEXP-3287) and GDM (E-MEXP-3349) for mRNAs, and T1D (E-MEXP-3409), T2D (E-MEXP-3373) and GDM (E-MEXP-3382) for miRNAs. Although these files were recorded into the ArrayExpress database at different dates, the hybridizations were performed at the same period but not in the same date, using the same batch of mRNA and miRNA microarray platforms.
Data quantification and normalization
Data quantification and quality control were performed using the Feature Extraction (FE) software version 11.0 (Agilent Technologies). Expression data were loaded into R environment [76]. Background adjustment was done by subtracting median background values from the median expression values obtained by FE, and data were subsequently log-transformed. Finally, all the data distribution was normalized by quantile function using the Aroma light package [77, 78].
Rank products statistical analysis
Differentially expressed mRNAs and miRNAs were identified using the Rank products non-parametric method by the R package RankProd [79]. Genes were considered to be of interest if P-values and percentage of false positive predictions (pfp) were smaller than 0.001. This test was used to perform paired analysis between T1D versus GDM, T2D versus GDM and T1D versus T2D for mRNA and miRNA microarrays. Hierarchical clustering of differentially expressed genes was performed using average linkage criterion and Pearson uncentered distance as a metric.
MicroRNA target prediction
We searched for miRNA targets using publicly available predictions. Variables selected for best targets included those present in at least three algorithms, using miRWalk, miRanda, RNAhybrid, and TargetScan [23–26]. The algorithm GenMir++ (Generative model for miRNA regulation) [80] was used to evaluate interactions of the differentially expressed mRNAs and miRNAs, according to their expression profiles, to identify the best candidate interactions using Bayesian inferences. The mRNA/miRNA networks were performed using Cytoscape software [81]. Intending to find association of differentially expressed miRNAs with complex disease, miR2Disease databases [74] were also used.
ROC curve
Differentially expressed miRNAs were further ranked according to sensitivity and specificity to find the best candidates between T1D vs. T2D, T1D vs. GDM and T2D vs. GDM groups. The criteria for biomarker selection were sensitivity and specificity ≥ 80%, as determined using general logistic models (glm) in the R environment followed by ROC curve analyses using ROCR package [82].
Abbreviations
- miRNA:
-
microRNA
- T1D:
-
Type 1 diabetes
- T2D:
-
Type 2 diabetes
- GDM:
-
Gestational diabetes mellitus
- AUC:
-
Area under curve
- PBMCs:
-
Peripheral mononuclear cells.
References
Artmann S, Jung K, Bleckmann A, Beissbarth T: Detection of simultaneous group effects in microRNA expression and related target gene sets. PLoS One. 2012, 7: e38365-
Liang Y, Ridzon D, Wong L, Chen C: Characterization of microRNA expression profiles in normal human tissues. BMC Genomics. 2007, 8: 166-
Hezova R, Slaby O, Faltejskova P, Mikulkova Z, Buresova I, Raja KR, Hodek J, Ovesna J, Michalek J: MicroRNA-342, microRNA-191 and microRNA-510 are differentially expressed in T regulatory cells of type 1 diabetic patients. Cell Immunol. 2010, 260: 70-74.
Dai R, Zhang Y, Khan D, Heid B, Caudell D, Crasta O, Ahmed SA: Identification of a common lupus disease-associated microRNA expression pattern in three different murine models of lupus. PLoS One. 2010, 5: e14302-
Grimson A, Farh KK, Johnston WK, Garrett-Engele P, Lim LP, Bartel DP: MicroRNA targeting specificity in mammals: determinants beyond seed pairing. Mol Cell. 2007, 27: 91-105.
Krek A, Grün D, Poy MN, Wolf R, Rosenberg L, Epstein EJ, MacMenamin P, da Piedade I, Gunsalus KC, Stoffel M, Rajewsky N: Combinatorial microRNA target predictions. Nat Genet. 2005, 37: 495-500.
Griffiths-Jones S, Grocock RJ, van Dongen S, Bateman A, Enright AJ: miRBase: microRNA sequences, targets and gene nomenclature. Nucleic Acids Res. 2006, 34: D140-D144.
Zhang B, Farwell MA: microRNAs: a new emerging class of players for disease diagnostics and gene therapy. J Cell Mol Med. 2008, 12: 3-21.
Pauley KM, Cha S, Chan EK: MicroRNA in autoimmunity and autoimmune diseases. J Autoimmun. 2009, 32: 189-194.
Ciesla M, Skrzypek K, Kozakowska M, Loboda A, Jozkowicz A, Dulak J: MicroRNAs as biomarkers of disease onset. Anal Bioanal Chem. 2011, 401: 2051-2061.
World Health Organization.http://www.idf.org/diabetesatlas,
Prentki M, Vischer S, Glennon MC, Regazzi R, Deeney JT, Corkey BE: Malonyl-CoA and long chain acyl-CoA esters as metabolic coupling factors in nutrient-induced insulin secretion. J Biol Chem. 1992, 267: 5802-5810.
Poy MN, Eliasson L, Krutzfeldt J, Kuwajima S, Ma X, Macdonald PE, Pfeffer S, Tuschl T, Rajewsky N, Rorsman P, Stoffel M: A pancreatic islet-specific microRNA regulates insulin secretion. Nature. 2004, 432: 226-230.
Plaisance V, Abderrahmani A, Perret-Menoud V, Jacquemin P, Lemaigre F, Regazzi R: MicroRNA-9 controls the expression of granuphilin/Slp4 and the secretory response of insulin-producing cells. J Biol Chem. 2006, 281: 26932-26942.
El Ouaamari A, Baroukh N, Martens GA, Lebrun P, Pipeleers D, van Obberghen E: MiR-375 targets 3′-phosphoinositide-dependent protein kinase-1 and regulates glucose-induced biological responses in pancreatic beta-cells. Diabetes. 2008, 57: 2708-2717.
Tang X, Tang G, Ozcan S: Role of microRNAs in diabetes. Biochim Biophys Acta. 2008, 1779: 697-701.
Tang X, Muniappan L, Tang G, Ozcan S: Identification of glucose-regulated miRNAs from pancreatic beta cells reveals a role for miR-30d in insulin transcription. RNA. 2009, 15: 287-293.
Hennessy E, Clynes M, Jeppesen PB, O’Driscoll L: Identification of microRNAs with a role in glucose stimulated insulin secretion by expression profiling of MIN6 cells. Biochem Biophys Res Commun. 2010, 396: 457-462.
Lynn FC: Meta-regulation: microRNA regulation of glucose and lipid metabolism. Trends Endocrinol Metab. 2009, 20: 452-459.
Bang-Berthelsen CH, Pedersen L, Fløyel T, Hagedorn PH, Gylvin T, Pociot F: Independent component and pathway-based analysis of miRNA-regulated gene expression in a model of type 1 diabetes. BMC Genomics. 2011, 12: 97-
Zampetaki A, Kiechl S, Drozdov I, Willeit P, Mayr U, Prokopi M, Mayr A, Weger S, Oberhollenzer F, Bonora E, Shah A, Willeit J, Mayr M: Plasma microRNA profiling reveals loss of endothelial miR-126 and other microRNAs in type 2 diabetes. Circ Res. 2010, 107: 810-817.
Karolina DS, Armugam A, Tavintharan S, Wong MT, Lim SC, Sum CF, Jeyaseelan K: MicroRNA 144 impairs insulin signaling by inhibiting the expression of insulin receptor substrate 1 in type 2 diabetes mellitus. PLoS One. 2011, 6: e22839-
RNAhybrid.http://bibiserv.techfak.uni-bielefeld.de/rnahybrid,
TargetScan.http://targetscan.org,
Metz CE: Basic principles of ROC analysis. Semin Nucl Med. 1978, 8: 283-298.
Griner PF, Mayewski RJ, Mushlin AI, Greenland P: Selection and interpretation of diagnostic tests and procedures. Principles and applications. Ann Intern Med. 1981, 94: 557-592.
Collares CV, Evangelista AF, Xavier DJ, Takahashi P, Almeida R, Macedo C, Manoel-Caetano F, Foss MC, Foss-Freitas MC, Rassi DM, Sakamoto-Hojo ET, Passos GA, Donadi EA: Transcriptome meta-analysis of peripheral lymphomononuclear cells indicates that gestational diabetes is closer to type 1 diabetes than to type 2 diabetes mellitus. Mol Biol Rep. 2013,  -in press
Krützfeldt J, Stoffel M: MicroRNAs: a new class of regulatory genes affecting metabolism. Cell Metab. 2006, 4: 9-12.
He A, Zhu L, Gupta N, Chang Y, Fang F: Overexpression of micro ribonucleic acid 29, highly up-regulated in diabetic rats, leads to insulin resistance in 3T3-L1 adipocytes. Mol Endocrinol. 2007, 21: 2785-2794.
Arner E, Mejhert N, Kulyté A, Balwierz PJ, Pachkov M, Cormont M, Lorente-Cebrián S, Ehrlund A, Laurencikiene J, Hedén P, Dahlman-Wright K, Tanti JF, Hayashizaki Y, Rydén M, Dahlman I, van Nimwegen E, Daub CO, Arner P: Adipose Tissue MicroRNAs as regulators of CCL2 production in human obesity. Diabetes. 2012, 61: 1986-1993.
Wang S, Aurora AB, Johnson BA, Qi X, McAnally J, Hill JA, Richardson JA, Bassel-Duby R, Olson EN: The endothelial- specific microRNA miR-126 governs vascular integrity and angiogenesis. Dev Cell. 2008, 15: 261-271.
Meng S, Cao JT, Zhang B, Zhou Q, Shen CX, Wang CQ: Downregulation of microRNA-126 in endothelial progenitor cells from diabetes patients, impairs their functional properties, via target gene Spred-1. J Mol Cell Cardiol. 2012, 53: 64-72.
Zampetaki A, Willeit P, Tilling L, Drozdov I, Prokopi M, Renard JM, Mayr A, Weger S, Schett G, Shah A, Boulanger CM, Willeit J, Chowienczyk PJ, Kiechl S, Mayr M: Prospective study on circulating MicroRNAs and risk of myocardial infarction. J Am Coll Cardiol. 2012, 60: 290-299.
Kato M, Zhang J, Wang M, Lanting L, Yuan H, Rossi JJ, Natarajan R: MicroRNA-192 in diabetic kidney glomeruli and its function in TGF-beta-induced collagen expression via inhibition of E-box repressors. Proc Natl Acad Sci U S A. 2007, 104: 3432-3437.
Ling HY, Ou HS, Feng SD, Zhang XY, Tuo QH, Chen LX, Zhu BY, Gao ZP, Tang CK, Yin WD, Zhang L, Liao DF: Changes in microRNA profile and effects of miR-320 in insulin-resistant 3T3–L1 adipocytes. Clin Exp Pharmacol Physiol. 2009, 36: e32-e39.
Zhao H, Guan J, Lee HM, Sui Y, He L, Siu JJ, Tse PP, Tong PC, Lai FM, Chan JC: Up-regulated pancreatic tissue microRNA-375 associates with human type 2 diabetes through beta-cell deficit and islet amyloid deposition. Pancreas. 2010, 39: 843-846.
Herrera BM, Lockstone HE, Taylor JM, Ria M, Barrett A, Collins S, Kaisaki P, Argoud K, Fernandez C, Travers ME, Grew JP, Randall JC, Gloyn AL, Gauguier D, McCarthy MI, Lindgren CM: Global microRNA expression profiles in insulin target tissues in a spontaneous rat model of type 2 diabetes. Diabetologia. 2010, 53: 1099-1109.
Karolina DS, Tavintharan S, Armugam A, Sepramaniam S, Pek SL, Wong MT, Lim SC, Sum CF, Jeyaseelan K: Circulating miRNA profiles in patients with metabolic syndrome. J Clin Endocrinol Metab. 2012, 97: E2271-E2276.
Roggli E, Gattesco S, Caille D, Briet C, Boitard C, Meda P, Regazzi R: Changes in microRNA expression contribute to pancreatic β-cell dysfunction in prediabetic NOD mice. Diabetes. 2012, 61: 1742-1751.
Silva VA, Polesskaya A, Sousa TA, Corrêa VM, André ND, Reis RI, Kettelhut IC, Harel-Bellan A, De Lucca FL: Expression and cellular localization of microRNA-29b and RAX, an activator of the RNA-dependent protein kinase (PKR), in the retina of streptozotocin-induced diabetic rats. Mol Vis. 2011, 17: 2228-2240.
Wang G, Kwan BC, Lai FM, Chow KM, Li PK, Szeto CC: Urinary miR-21, miR-29, and miR-93: novel biomarkers of fibrosis. Am J Nephrol. 2012, 36: 412-418.
Ding S, Liang Y, Zhao M, Liang G, Long H, Zhao S, Wang Y, Yin H, Zhang P, Zhang Q, Lu Q: Decreased microRNA-142-3p/5p expression causes CD4+ T cell activation and B cell hyperstimulation in systemic lupus erythematosus. Arthritis Rheum. 2012, 64: 2953-2963.
Wang XS, Gong JN, Yu J, Wang F, Zhang XH, Yin XL, Tan ZQ, Luo ZM, Yang GH, Shen C, Zhang JW: MicroRNA-29a and microRNA-142-3p are regulators of myeloid differentiation and acute myeloid leukemia. Blood. 2012, 119: 4992-5004.
Lv M, Zhang X, Jia H, Li D, Zhang B, Zhang H, Hong M, Jiang T, Jiang Q, Lu J, Huang X, Huang B: An oncogenic role of miR-142-3p in human T-cell acute lymphoblastic leukemia (T-ALL) by targeting glucocorticoid receptor-α and cAMP/PKA pathways. Leukemia. 2012, 26: 769-777.
Lin RJ, Xiao DW, Liao LD, Chen T, Xie ZF, Huang WZ, Wang WS, Jiang TF, Wu BL, Li EM, Xu LY: MiR-142-3p as a potential prognostic biomarker for esophageal squamous cell carcinoma. J Surg Oncol. 2012, 105: 175-182.
Wu L, Cai C, Wang X, Liu M, Li X, Tang H: MicroRNA-142-3p, a new regulator of RAC1, suppresses the migration and invasion of hepatocellular carcinoma cells. FEBS Lett. 2011, 585: 1322-1330.
Li S, Moffett HF, Lu J, Werner L, Zhang H, Ritz J, Neuberg D, Wucherpfennig KW, Brown JR, Novina CD: MicroRNA expression profiling identifies activated B cell status in chronic lymphocytic leukemia cells. PLoS One. 2011, 6: e16956-
Birks DK, Barton VN, Donson AM, Handler MH, Vibhakar R, Foreman NK: Survey of MicroRNA expression in pediatric brain tumors. Pediatr Blood Cancer. 2011, 56: 211-216.
Murakami Y, Yasuda T, Saigo K, Urashima T, Toyoda H, Okanoue T, Shimotohno K: Comprehensive analysis of microRNA expression patterns in hepatocellular carcinoma and non-tumorous tissues. Oncogene. 2006, 25: 2537-2545.
Jiang J, Gusev Y, Aderca I, Mettler TA, Nagorney DM, Brackett DJ, Roberts LR, Schmittgen TD: Association of MicroRNA expression in hepatocellular carcinomas with hepatitis infection, cirrhosis, and patient survival. Clin Cancer Res. 2008, 14: 419-427.
Iorio MV, Visone R, Di Leva G, Donati V, Petrocca F, Casalini P, Taccioli C, Volinia S, Liu CG, Alder H, Calin GA, Ménard S, Croce CM: MicroRNA signatures in human ovarian cancer. Cancer Res. 2007, 67: 8699-8707.
Yu T, Wang XY, Gong RG, Li A, Yang S, Cao YT, Wen YM, Wang CM, Yi XZ: The expression profile of microRNAs in a model of 7,12-dimethyl-benz[a]anthrance-induced oral carcinogenesis in Syrian hamster. J Exp Clin Cancer Res. 2009, 28: 64-
Van der Auwera I, Limame R, van Dam P, Vermeulen PB, Dirix LY, Van Laere SJ: Integrated miRNA and mRNA expression profiling of the inflammatory breast cancer subtype. Br J Cancer. 2010, 103: 532-541.
Ronchetti D, Lionetti M, Mosca L, Agnelli L, Andronache A, Fabris S, Deliliers GL, Neri A: An integrative genomic approach reveals coordinated expression of intronic miR-335, miR-342, and miR-561 with deregulated host genes in multiple myeloma. BMC Med Genomics. 2008, 1: 37-
Saba R, Goodman CD, Huzarewich RL, Robertson C, Booth SA: A miRNA signature of prion induced neurodegeneration. PLoS One. 2008, 3: e3652-
Chartoumpekis DV, Zaravinos A, Ziros PG, Iskrenova RP, Psyrogiannis AI, Kyriazopoulou VE, Habeos IG: Differential expression of microRNAs in adipose tissue after long-term high-fat diet-induced obesity in mice. PLoS One. 2012, 7: e34872-
Al-Quraishy S, Dkhil MA, Delic D, Abdel-Baki AA, Wunderlich F: Organ-specific testosterone-insensitive response of miRNA expression of C57BL/6 mice to plasmodium chabaudi malaria. Parasitol Res. 2012, 111: 1093-1101.
Johnson SM, Grosshans H, Shingara J, Byrom M, Jarvis R, Cheng A, Labourier E, Reinert KL, Brown D, Slack FJ: RAS is regulated by the let-7 microRNA family. Cell. 2005, 120: 635-647.
Lee Y, Dutta A: The tumor suppressor microRNA let-7 represses the HMGA2 oncogene. Genes Dev. 2007, 21: 1025-1030.
Frost RJ, Olson EN: Control of glucose homeostasis and insulin sensitivity by the Let-7 family of microRNAs. Proc Natl Acad Sci U S A. 2011, 108: 21075-21080.
Schopman NC, Heynen S, Haasnoot J, Berkhout B: A miRNA-tRNA mix-up: tRNA origin of proposed miRNA. RNA Biol. 2010, 7: 573-576.
Zhao C, Dong J, Jiang T, Shi Z, Yu B, Zhu Y, Chen D, Xu J, Huo R, Dai J, Xia Y, Pan S, Hu Z, Sha J: Early second-trimester serum miRNA profiling predicts gestational diabetes mellitus. PLoS One. 2011, 6: e23925-
Xie H, Lim B, Lodish HF: MicroRNAs induced during adipogenesis that accelerate fat cell development are downregulated in obesity. Diabetes. 2009, 58: 1050-1057.
Bloomston M, Frankel WL, Petrocca F, Volinia S, Alder H, Hagan JP, Liu CG, Bhatt D, Taccioli C, Croce CM: MicroRNA expression patterns to differentiate pancreatic adenocarcinoma from normal pancreas and chronic pancreatitis. JAMA. 2007, 297: 1901-1908.
Pando R, Even-Zohar N, Shtaif B, Edry L, Shomron N, Phillip M, Gat-Yablonski G: MicroRNAs in the growth plate are responsive to nutritional cues: association between miR-140 and SIRT1. J Nutr Biochem. 2012, 23: 1474-1481.
Cohen HY, Miller C, Bitterman KJ, Wall NR, Hekking B, Kessler B, Howitz KT, Gorospe M, de Cabo R, Sinclair DA: Calorie restriction promotes mammalian cell survival by inducing the SIRT1 deacetylase. Science. 2004, 305: 390-392.
Chen D, Bruno J, Easlon E, Lin SJ, Cheng HL, Alt FW, Guarente L: Tissue-specific regulation of SIRT1 by calorie restriction. Genes Dev. 2008, 22: 1753-1757.
Klöting N, Berthold S, Kovacs P, Schön MR, Fasshauer M, Ruschke K, Stumvoll M, Blüher M: MicroRNA expression in human omental and subcutaneous adipose tissue. PLoS One. 2009, 4: e4699-
Joglekar MV, Patil D, Joglekar VM, Rao GV, Reddy DN, Mitnala S, Shouche Y, Hardikar AA: The miR-30 family microRNAs confer epithelial phenotype to human pancreatic cells. Islets. 2009, 1: 137-147.
Zhou B, Li C, Qi W, Zhang Y, Zhang F, Wu JX, Hu YN, Wu DM, Liu Y, Yan TT, Jing Q, Liu MF, Zhai QW: Downregulation of miR-181a upregulates sirtuin-1 (SIRT1) and improves hepatic insulin sensitivity. Diabetologia. 2012, 55: 2032-2043.
Nielsen LB, Wang C, Sørensen K, Bang-Berthelsen CH, Hansen L, Andersen ML, Hougaard P, Juul A, Zhang CY, Pociot F, Mortensen HB: Circulating levels of microRNA from children with newly diagnosed type 1 diabetes and healthy controls: evidence that miR-25 associates to residual beta-cell function and glycaemic control during disease progression. Exp Diabetes Res. 2012, 2012: 896362-
Mir2disease.http://www.mir2disease.org,
ArrayExpress.http://www.ebi.ac.uk/arrayexpress,
R environment.http://www.r-project.org,
Bolstad B, Irizarry R, Astrand M, Speed T: A comparison of normalization methods for high density oligonucleotide array data based on variance and bias. Bioinformatics. 2003, 19: 185-193.
Gautier L, Cope L, Bolstad BM, Irizarry RA: affy–analysis of Affymetrix GeneChip data at the probe level. Bioinformatics. 2004, 20: 307-315.
Hong F, Breitling R, McEntee CW, Wittner BS, Nemhauser JL, Chory J: RankProd: a bioconductor package for detecting differentially expressed genes in meta-analysis. Bioinformatics. 2006, 22: 2825-2827.
Huang JC, Morris QD, Frey BJ: Bayesian inference of MicroRNA targets from sequence and expression data. J Comput Biol. 2007, 14: 550-563.
Shannon P, Markiel A, Ozier O, Baliga NS, Wang JT, Ramage D, Amin N, Schwikowski B, Ideker T: Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res. 2003, 13: 2498-2504.
ROCR package.http://cran.r-project.org/web/packages/ROCR/index.html,
Acknowledgements
This study was funded by the Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP - (FAPESP #2008/56594-8, FAPESP #2010/05622-1, FAPESP #210/00932-2, FAPESP #2010/12069-7), Conselho Nacional de Desenvolvimento CientÃfico e Tecnológico (CNPq # 563731/2010-9), and NAP-DIN (Núcleo de Apoio à Pesquisa em Doenças Inflamatórias).
Author information
Authors and Affiliations
Corresponding author
Additional information
Competing interests
The authors declare that they have no competing interests.
Authors’ contributions
CVAC and EAD wrote the paper. CVAC, AFE, GAP and EAD conceived the study and participated in its design and coordination. ETSK, GAP and EAD provided material and analytical tools. AFE, DJX and DMR performed patient selection, sample collection and cell separation. MCFF and MFF were responsible for diabetic patient’s treatment and clinical information data. AFE and DJX performed RNA extraction and sample quality control. CVAC and DJX were responsible for hybridization and acquisition of microarray data. CVAC, AFE, TA and DP performed data analysis. All authors contributed to and approved the final manuscript.
Cristhianna VA Collares, Adriane F Evangelista contributed equally to this work.
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.
Rights and permissions
Open Access This article is published under license to BioMed Central Ltd. This is an Open Access article is distributed under the terms of the Creative Commons Attribution License ( https://creativecommons.org/licenses/by/2.0 ), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
About this article
Cite this article
Collares, C.V., Evangelista, A.F., Xavier, D.J. et al. Identifying common and specific microRNAs expressed in peripheral blood mononuclear cell of type 1, type 2, and gestational diabetes mellitus patients. BMC Res Notes 6, 491 (2013). https://doi.org/10.1186/1756-0500-6-491
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/1756-0500-6-491
Keywords
- Type 1 diabetes mellitus
- Type 2 diabetes mellitus
- Gestational diabetes mellitus
- microRNA
- mRNA
- Microarray