Skip to main content
Download PDF
  • Research article
  • Open access
  • Published:

Competing endogenous RNA network analysis explores the key lncRNAs, miRNAs, and mRNAs in type 1 diabetes

BMC Medical Genomics volume 14, Article number: 35 (2021) Cite this article

Abstract

Background

Type 1 diabetes (T1D, named insulin-dependent diabetes) has a relatively rapid onset and significantly decreases life expectancy. This study is conducted to reveal the long non-coding RNA (lncRNA)-microRNA (miRNA)-mRNA regulatory axises implicated in T1D.

Methods

The gene expression profile under GSE55100 (GPL570 and GPL8786 datasets; including 12 T1D samples and 10 normal samples for each dataset) was extracted from Gene Expression Omnibus database. Using limma package, the differentially expressed mRNAs (DE-mRNAs), miRNAs (DE-miRNAs), and lncRNAs (DE-lncRNAs) between T1D and normal samples were analyzed. For the DE-mRNAs, the functional terms were enriched by DAVID tool, and the significant pathways were enriched using gene set enrichment analysis. The interactions among DE-lncRNAs, DE-miRNAs and DE-mRNAs were predicted using mirwalk and starbase. The lncRNA-miRNA-mRNA interaction network analysis was visualized by Cytoscape. The key genes in the interaction network were verified by quantitatively real-time PCR.

Results

In comparison to normal samples, 236 DE-mRNAs, 184 DE-lncRNAs, and 45 DE-miRNAs in T1D samples were identified. For the 236 DE-mRNAs, 16 Gene Ontology (GO)_biological process (BP) terms, four GO_cellular component (CC) terms, and 57 significant pathways were enriched. A network involving 36 DE-mRNAs, 8 DE- lncRNAs, and 15 DE-miRNAs was built, such as TRG-AS1—miR-23b/miR-423—PPM1L and GAS5—miR-320a/miR-23b/miR-423—SERPINA1 regulatory axises. Quantitatively real-time PCR successfully validated the expression levels of TRG-AS1- miR-23b -PPM1L and GAS5-miR-320a- SERPINA1.

Conclusion

TRG-AS1—miR-23b—PPM1L and GAS5—miR-320a—SERPINA1 regulatory axises might impact the pathogenesis of T1D.

Peer Review reports

Background

As a metabolic disorder characterized by hyperglycemia, diabetes is induced by secretion deficiency or impaired biological function of insulin [1, 2]. Diabetes can result in impairment and dysfunction of multiple tissues, including heart, blood vessels, eyes, nerves, and kidneys [3,4,5]. According to the WHO classification of diabetes mellitus in 2019, diabetes are classified as type 1 diabetes (T1D), type 2 diabetes (T2D), hybrid forms of diabetes, other specific types, such as monogenic diabetes, unclassified diabetes and hyperglyacemia first detected during pregnancy [6, 7]. The symptoms of T1D mainly include thirst, polydipsia, polyuria, polyphagia, fatigue, and rapid weight loss [8].An annual 3%—4% increase in the incidence of T1D in childhood was estimated in several developed countries [9]. Besides, T1D decreases approximately 11–13 years of life expectancy in developed countries, and even more time in developing countries [10]. Quality of life of T1D patients was significantly decreased than general population. Therefore, investigation in the pathogenesis of T1D is warranted to improve the prognosis of T1D patients.

Considerable efforts have been made on the molecular mechanism of T1D during the past decades. The roles of non-coding RNAs, including microRNA (miRNA), long non-coding RNAs (lncRNAs) in the regulation of T1D have only recently been recognized [11]. Through reducing C–C motif chemokine receptor 2 expression, miR-125a-5p limits the migration and function of regulatory T cells (Tregs) in the pancreas of T1D patients [12]. Activity enhancement of miR-181a promotes the expression of nuclear factor of activated T cells 5 (NFAT5) while suppresses the induction of forkhead box P3 + Tregs; therefore, miR-181a/NFAT5 axis may provide targets for limiting islet autoimmunity in T1D patients [13, 14]. The sera level of miR-375 is decreased in T1D patients, indicating that the changed circulating level of miR-375 may be a valuable marker of inflammation and metabolic alterations in T1D [15]. The crosstalk between p38 mitogen-activated protein kinase signaling pathway and the lncRNA metastasis associated lung adenocarcinoma transcript 1 (MALAT1) is correlated with the endothelial cell function in diabetic rats, and MALAT1 inhibition may be a promising approach for anti-angiogenic treatment of diabetic microvascular complications [16]. Besides, lncRNAs MEG3, TUG1 and PVT1 were also being reported to contribute in the pathophysiology of T1D and T1D-associated complications [17,18,19]. However, only a small part of non-coding RNAs in T1D have been revealed, and more comprehensive understanding of T1D should be developed to facilitate design of preventive therapeutic modalities.

MiiRNAs can silence gene expression by binding mRNAs, while competing endogenous RNAs (ceRNAs), such as lncRNAs can bind miRNAs via miRNA response elements competitively to regulate gene expression [20, 21]. In this study, the microarray dataset of T1D was searched from Gene Expression Omnibus (GEO) database. Afterwards, differential expression analysis, enrichment analysis, and lncRNA-miRNA-mRNA network (ceRNA network) analysis were executed to investigate the important regulatory axises involved in T1D.

Methods

Data source and data preprocessing

From GEO database (http://www.ncbi.nlm.nih.gov/geo/), the gene expression profile of T1D (accession number: GSE55100, including GPL570 and GPL8786 datasets) was downloaded. This dataset was deposited by Gu et al. [22]. There were 12 peripheral blood mononuclear cell (PBMC) samples from T1D patients and 10 PBMC samples from normal controls in both of GPL570 (including mRNA and lncRNA expression data) and GPL8786 (including miRNA expression data) datasets. The clinical characteristics of the 22 samples were listed in Table 1.

Table 1 Clinical characteristics of type 1 diabetes patients
Full size table

Using the Release 26 (GRCh38.p10) reference genome in GENCODE database (https://www.gencodegenes.org/human/) [23], sequence alignment was conducted. Only the unique alignment sequences were remained. Based on corresponding GTF annotation files, mRNAs (with annotation information of “protein coding”) and lncRNAs (with annotation information of “antisense”, “sense_intronic”, “lincRNA”, “sense_overlapping”, “processed_transcript”, “3prime_overlapping_ncRNA”, and “non_coding”) were identified, respectively. Subsequently, the probes without relevant gene symbols were filtered out. For multiple probes that mapped to one gene symbol, their average value was used as the final expression value of this gene.

Differential expression analysis

For GPL570 and GPL8786 datasets, differential expression analysis between T1D and normal samples were performed using limma package (http://www.bioconductor.org/packages/2.9/bioc/html/limma.html) [24] in R. The corresponding p values of all genes were obtained through statistical test and were conducted with multiple test correction using Benjamini & Hochberg method [25]. The |log fold change (FC)|≥ 0.5 and adjusted p value < 0.05 were considered as the thresholds for obtaining differentially expressed mRNAs (DE-mRNAs). Meanwhile, the p value ≤ 0.05 was regarded as the cut-off criterion for selecting differentially expressed miRNAs (DE-miRNAs) and differentially expressed lncRNAs (DE-lncRNAs).

Enrichment analysis for the DE-mRNAs

The Gene Ontology (GO)_biological process (BP), GO_molecular function (MF), and GO_cellular component (CC) functional terms involving the DE-mRNAs were analyzed using DAVID tool (version 6.7, https://david-d.ncifcrf.gov/) [26], and the corresponding results were visualized by the R package GOplot (http://cran.r-project.org/web/packages/GOplot) [27]. Using gene set enrichment analysis (GSEA) [28], Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment for the DE-mRNAs was carried out. The adjusted p value < 0.05 was used for screening the significant results.

Construction of lncRNA-miRNA-mRNA interaction network

Using corr.test (parameters: ci = F, adjust = "BH") in R package psych [29], the Pearson correlation coefficients [30] of the expression values of the DE-lncRNAs and the DE-mRNAs were calculated. The lncRNA-mRNA pairs with |r|≥ 0.7 and adjusted p value < 0.05 were selected and a lncRNA-mRNA co-expression network was developed By using Cytoscape software (http://www.cytoscape.org) [31].

Using starbase database (version 3.0, http://starbase.sysu.edu.cn/) [32], lncRNA-miRNA interactions were predicted for the DE-lncRNAs in the lncRNA-mRNA co-expression network. Simultaneously, miRNA-mRNA interactions were predicted for the DE-mRNAs in the lncRNA-mRNA co-expression network by using with mirwalk database (version 3.0, http://mirwalk.umm.uni-heidelberg.de/) [33]. Finally, a lncRNA-miRNA-mRNA interaction network was visualized by Cytoscape software [31] by integrating the DE-lncRNAs and DE-mRNA regulated by the same DE-miRNA. Furthermore, GO and KEGG enrichment analysis for the mRNAs involved in the interaction network were conducted using DAVID tool [26].

Validation of key genes using quantitatively real-time PCR (qRT-PCR)

T1D patients were recruited from the department of Endocrinology, China-Japan Union Hospital of Jilin University. Age and sex matched healthy controls were recruited from volunteers of physical examination. The experiments were approved by Ethics Committee of China-Japan Union Hospital of Jilin University [No.(2020)linshen(20,201,127)]. Written informed consent was received from all participants. Total RNA were extracted from the PBMCs of T1D patients and normal control samples using RNAiso Plus (TaKaRa, Shiga, Japan). The total RNA was reversed into cDNA with primeScript RT Master Mix (TaKaRa) and amplified on ABI ViiA7 (ThermoFisher, USA). GAPDH was used as the internal reference of DE-lncRNAs and DE-mRNAs and U6 was regarded as the internal reference of DE-miRNAs. The primers were listed in Table 2.

Table 2 The primers used for qRT-PCR
Full size table

Statistical analysis

The data of qRT-PCR was processed by GraphPad Prism 5.0 (San Diego, CA, USA). Data were presented as mean ± standard deviation. Differences between two groups were determined by t-test. P < 0.05 was regarded as statistical significance level.

Results

Differential expression analysis

A total of 1671 mRNAs, 1990 lncRNAs, and 533 miRNAs were identified from the gene expression profile of GSE55100. Compared with normal samples, there were 236 DE-mRNAs (121 up-regulated and 115 down-regulated), 184 DE-lncRNAs (106 up-regulated and 78 down-regulated), and 45 DE-miRNAs (30 up-regulated and 15 down-regulated) in T1D samples (Fig. 1).

Fig. 1
figure 1

The results of differential expression analysis. a The volcano plot of the differentially expressed mRNAs (DE-mRNAs); b The expression heatmap of the top 10 up-regulated mRNAs and down-regulated mRNAs; c The volcano plot of the differentially expressed lncRNAs (DE-lncRNAs); d The expression heatmap of the top 10 up-regulated lncRNAs and down-regulated lncRNAs; e The volcano plot of the differentially expressed miRNAs (DE-miRNAs); f The expression heatmap of the top 10 up-regulated miRNAs and down-regulated miRNAs. FC, fold change. In volcano plots, red and green separately represent up-regulation and down-regulation. In expression heatmaps, red and blue separately represent type 1 diabetes (T1D) and normal samples

Full size image

Enrichment analysis for the DE-mRNAs

For the 236 DE-mRNAs, 16 GO_BP terms (such as “defense response”, adjusted p value = 5.01E-08; “immune response”, adjusted p value = 1.11E-04; and “response to wounding”, adjusted p value = 3.46E-04) and four GO_CC terms (such as “plasma membrane part”, adjusted p value = 2.72E-02; “intrinsic to plasma membrane”, adjusted p value = 2.90E-02; and “integral to plasma membrane”, adjusted p value = 3.46E-02) were enriched (Table 3).

Table 3 The Gene Ontology (GO) functional terms enriched for the differentially expressed mRNAs
Full size table

Meanwhile, 57 significant pathways were enriched for the DE-mRNAs, including 23 activated pathways (normalized enrichment score (NES) > 0; such as “Th1 and th2 cell differentiation”, adjusted p value = 6.56E-03; “T cell receptor signaling pathway”, adjusted p value = 6.56E-03; and “Biosynthesis of unsaturated fatty acids”, adjusted p value = 6.56E-03) and 34 suppressed pathways (NES < 0; such as “Leishmaniasis”, adjusted p value = 6.56E-03; “Legionellosis”, adjusted p value = 6.56E-03; and “Phagosome”, adjusted p value = 6.56E-03) (Table 4).

Table 4 The activated pathways and suppressed pathways enriched for the differentially expressed mRNAs (top 10 listed)
Full size table

Construction of lncRNA-miRNA-mRNA interaction network

In the lncRNA-mRNA co-expression network, there were 127 nodes (37 lncRNAs and 90 mRNAs) and 174 edges (Fig. 2). After integrating the results of lncRNA-miRNA interaction prediction and miRNA-mRNA interaction prediction, the lncRNA-miRNA-mRNA interaction network (including 36 mRNAs, 8 lncRNAs, and 15 miRNAs) was constructed (Fig. 3). Importantly, T cell receptor gamma locus antisense RNA 1 (TRG-AS1)—miR-23b/miR-423—protein phosphatase, Mg2 + /Mn2 + dependent 1L (PPM1L) and growth arrest specific 5 (GAS5) —miR-320a/miR-23b/miR-423—serpin family A member 1 (SERPINA1) regulatory axises were implicated in the interaction network.

Fig. 2
figure 2

The long non-coding RNA (lncRNA)-mRNA co-expression network. Red and yellow circles separately represent mRNAs and lncRNAs

Full size image
Fig. 3
figure 3

The long non-coding RNA (lncRNA)-miRNA-mRNA interaction network. Red, yellow, and blue circles represent mRNAs, lncRNAs, and miRNAs, respectively

Full size image

For the mRNAs involved in the interaction network, two GO_BP terms (“Adaptive immune response”, p value = 3.41E-02; “Regulation of immune response”, p value = 4.76E-02), three GO_CC terms (“Endoplasmic reticulum”, p value = 1.78E-02; “Integral component of membrane”, p value = 3.56E-02; “ER to Golgi transport vesicle”, p value = 4.74E-02), one GO_MF term (“Protein homodimerization activity”, p value = 4.78E-02), and one KEGG pathway (“Natural killer cell mediated cytotoxicity”, p value = 1.83E-02) were enriched (Table 5).

Table 5 The Gene Ontology (GO) functional terms and pathways enriched for the mRNAs involved in the long non-coding RNA (lncRNA)-miRNA-mRNA interaction network
Full size table

Validation of key genes using qRT-PCR

In order to verify the results of bioinformation analysis, the expression of TRG-AS1-miR-23b/miR-423-PPM1L and GAS5-miR-320a /miR-23b/ miR-423- SERPINA1 were verified in T1D patients and normal controls. According to the differential expression analysis, TRG-AS1, miR-320a and PPM1L were upregulated in T1D patients, while GAS1, miR-423, miR-23b and SERPINA1 were downregulated in T1D patients. As shown in Fig. 4, the expression levels of TRG-AS1- miR-23b -PPM1L and GAS5-miR-320a- SERPINA1 were consistent with bioinformation analysis (p < 0.05). No significant difference was detected in miR-423 between control group and T1D group (p > 0.05).

Fig. 4
figure 4

Validation of key genes using quantitatively real-time PCR. *p < 0.05, **p < 0.01 compared with the control group

Full size image

Discussion

In recent years, the critical roles of non-coding RNAs including lncRNAs and miRNAs in diabetes were recognized. However, little is known about the molecular mechanism of them in regulating the development of T1D. In this study, 236 DE-mRNAs (121 up-regulated and 115 down-regulated), 184 DE-lncRNAs (106 up-regulated and 78 down-regulated), and 45 DE-miRNAs (30 up-regulated and 15 down-regulated) between T1D samples and normal samples were selected. For the 236 DE-mRNAs, 16 GO_BP terms, four GO_CC terms, and 57 significant pathways were enriched. Moreover, by constructing an lncRNA-miRNA-mRNA co-expression network, two interactions, TRG-AS1—miR-23b/ miR-423—PPM1L and GAS5—miR-320a/miR-23b/miR-423—SERPINA1 regulatory axises were identified.

TRG-AS1 has been reported to be oncogenic in glioblastoma [34], hepatocellular carcinoma [35] and tongue squamous cell carcinoma [36]. It might act as a ceRNA to regulate miRNA-543-Yes-associated protein 1 in tongue squamous cell carcinoma [36], miR-4500-BACH1 in hepatocellular carcinoma [35] and miR-877-5p-SUZ12 in glioblastoma [34]. It also identified to be involved in repeated implantation failure [37]. However, its role in diabetes was not characterized previously. In our study, we predicted the miRNAs and mRNAs that might be regulated by TRG-AS1. The lncRNA-miRNA-mRNA network displayed that it could directly regulate 5 miRNAs including miR-423 and miR-23b, and multiple mRNAs, including PPM1L indirectly. MiR-23b/27b expression is decreased in the muscle stem cells of T2D patients, which exerts a pro-myogenic function via the p53 pathway [38]. Nuclear factor, erythroid 2/miR-423-5p axis can induce gluconeogenesis and hyperglycemia through inhibiting the family with sequence similarity 3 member A -adenosine triphosphate (ATP)-serine/threonine kinase Akt pathway, which is involved in the progression of T2D and nonalcoholic fatty liver disease [39]. Increased miR-320 impairs lipid metabolism and gluconeogenesis by targeting adiponectin receptor 1 (AdipoR1), and thus miR-320 may be taken as a possible target for T2D therapy [40]. The elevated levels of the T lymphocytes expressing gamma-delta T cell receptor are implicated in the islet autoimmune process in individuals at high risk of T1D, which may serve as a promising indicator for the development of T1D [41, 42]. PPM1L mediates inositol-requiring protein-1 phosphorylation and endoplasmic reticulum stress signaling, which is considered as a causal gene for metabolic abnormalities [43]. The macrophage-enriched network has a causal correlation with metabolic disease traits, which involved three obesity genes (including PPM1L) [44]. Therefore, we hypothesized that TRG-AS1 and PPM1L might also play roles in the mechanisms of T1D through the TRG-AS1—miR-23b/miR-423—PPM1L regulatory axis.

GAS5 is a member of 5′ terminal oligopyrimidine class which could regulate cell growth, proliferation, and survival [45]. Reduced serum levels of GAS5 are related to diabetes; therefore, serum GAS5 levels combined with other parameters may be used for identifying people at high risk of diabetes more accurately [46]. Shi et al. demonstrated that GAS5 regulates insulin signaling in adipocytes, and suggested it might be a potential target for T2DM [47]. GAS5 knockdown causes cell cycle arrest and impairs insulin synthesis and secretion in Min6 pancreatic β-cells, and thus GAS5 may function in maintaining the identity and function of β cells [48]. In this study, GAS5 was found decreased in T1D patients in both RNA-seq results and the qRT-PCR results, which is consistent with previous studies. We found GAS5 could regulate 7 miRNAs, including miR-320a, miR-23b and miR-423, and tens of mRNAs, including SERPINA1 in T1D. MiR-320 negatively mediates the expression of fibronectin, endothelin 1, and vascular endothelial growth factor via extracellular signal-regulated kinases 1 and 2 in high glucose-treated human umbilical vein endothelial cells, which may provide a novel approach for treating diabetic complications [49]. Wei et al. found miR-320 mimic correlated with impaired gluconeogenesis and lipid metabolism by regulating adipoR1 [40]. SERPINA1 is a serine protease inhibitor that could target elastase, plasmin, thrombin, trypsin, chymotrypsin, and plasminogen activator. It was shown to be decreased in serum of obese mice and human subjects and the imbalance between SERPINA1 and neutrophil elastase contributed to insulin resistance [50]. Lower levels of SERPINA1 selectively impaired the ATP-binding cassette transporter A1 cholesterol efflux capacity in T2D [51]. Other members in this family were also reported to participate in development of diabetes. Anti-SERPINB13 antibody contributes to Reg gene expression and beta cell proliferation, and the immunological response may hinder the progression of T1D [52]. The serum concentrations of SERPINA12 (vaspin) is increased in T2D patients, which may be a candidate marker for evaluating the risk of severe macrovascular complications and the status of old T2D patients [53]. Thus, GAS5—miR-320a/miR-23b/miR-423—SERPINA1 regulatory axis might also function in the pathogenesis of T1D.

Conclusion

In conclusion, 236 DE-mRNAs, 184 DE-lncRNAs, and 45 DE-miRNAs between T1D and normal samples were identified. Besides, TRG-AS1—miR-23b/miR-423—PPM1L and GAS5—miR-320a/miR-23b/miR-423—SERPINA1 regulatory axises might be related to the pathogenesis of T1D. Though we have validated the expression of these genes in T1D patients, further mechanically mechanism of these regulatory axises should be investigated by subsequent studies.

Availability of data and materials

The dataset analysed during the current study are available in the GEO repository with accession number of GSE55110 (https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE55100), GPL570 (https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GPL570), GPL8786 (https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GPL8786). The reference genome of Release 26 (GRCh38.p10) is available in GENCODE database (https://www.gencodegenes.org/human/). The interactions between lncRNA and miRNAs were obtained from starbase (http://starbase.sysu.edu.cn/). The miRNA-mRNA interactions were obtained from mirwalk database (http://mirwalk.umm.uni-heidelberg.de/).

Abbreviations

T1D:

Type 1 diabetes

lncRNA:

Long non-coding RNA

miRNA:

MicroRNA

DEmRNA:

Differentially expressed mRNAs

GO:

Gene ontology

BP:

Biological process

CC:

Cellular component

T2D:

Type 2 diabetes

Tregs:

Regulatory T cells

NFAT5:

Nuclear factor of activated T cells 5

MALAT1:

Metastasis Associated Lung Adenocarcinoma Transcript 1

ceRNAs:

Competing endogenous RNAs

GEO:

Gene Expression Omnibus

PBMC:

Peripheral blood mononuclear cell

FC:

Fold change

GSEA:

Gene set enrichment analysis

KEGG:

Kyoto Encyclopedia of Genes and Genomes

qRT-PCR:

Quantitatively real-time PCR

NES:

Normalized enrichment score

TRG-AS1:

T cell receptor gamma locus antisense RNA 1

PPM1L:

Protein Phosphatase: Mg2 + /Mn2 + Dependent 1 l

GAS5:

Growth arrest specific 5

SERPINA1:

Serpin family A member 1

ATP:

Adenosine triphosphate

AdipoR1:

Adiponectin Receptor 1

References

  1. Puchulu FM. Definition, diagnosis and classification of diabetes mellitus. Diabetes Care. 2018;33(Suppl 1):S62–9.

    Google Scholar 

  2. El-Sappagh S, Elmogy M, Ali F, Kwak KS. A case-base fuzzification process: diabetes diagnosis case study. Soft Comput. 2018;9:5815–34.

    Google Scholar 

  3. Dong Y, Fernandes C, Liu Y, Wu Y, Wu H, Brophy ML, Deng L, Song K, Wen A, Wong S. Role of endoplasmic reticulum stress signalling in diabetic endothelial dysfunction and atherosclerosis. Diabetes Vasc Dis Res. 2017;14(1):14–23.

  4. Feinkohl I, Winterer G, Pischon T. Diabetes is associated with risk of postoperative cognitive dysfunction: A meta-analysis. Diabetes/Metab Res Rev. 2017;33(5):e2884.

    Article  CAS  Google Scholar 

  5. Mihardja L, Delima D, Massie RGA, Karyana M, Nugroho P, Yunir E. Prevalence of kidney dysfunction in diabetes mellitus and associated risk factors among productive age Indonesian. J Diabetes Metab Disord. 2018;17(4):53–61.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  6. Rossi G, Association AD. Diagnosis and classification of diabetes mellitus. Diabetes Care. 2018;33(Suppl 1):S62–9.

    Google Scholar 

  7. Classification of diabetes mellitus. https://www.who.int/publications/i/item/classification-of-diabetes-mellitus. Accessed 2 Dec 2020.

  8. Lundgren M, Jonsdottir B, Larsson HE, Group FS. Effect of screening for type 1 diabetes on early metabolic control: the DiPiS study. Diabetologia. 2018;62(1):53–7.

    Google Scholar 

  9. Patterson CC, Dahlquist GG, Gyürüs E, Green A, Soltész G. Incidence trends for childhood type 1 diabetes in Europe during 1989–2003 and predicted new cases 2005–20: a multicentre prospective registration study. Lancet. 2009;373(9680):2027–33.

    Article  PubMed  Google Scholar 

  10. Livingstone SJ, Levin D, Looker HC, Lindsay RS, Wild SH, Joss N, Leese G, Leslie P, McCrimmon RJ, Metcalfe W, et al. Estimated life expectancy in a scottish cohort with type 1 diabetes, 2008–2010. JAMA. 2015;313(1):37–44.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  11. Taheri M, Eghtedarian R, Dinger ME, Ghafouri-Fard S. Emerging roles of non-coding RNAs in the pathogenesis of type 1 diabetes mellitus. Biomed Pharmacother. 2020;129:110509.

    Article  CAS  PubMed  Google Scholar 

  12. Sebastiani G, Ventriglia G, Stabilini A, Socci C, Morsiani C, Laurenzi A, Nigi L, Formichi C, Mfarrej B, Petrelli A. Regulatory T-cells from pancreatic lymphnodes of patients with type-1 diabetes express increased levels of microRNA miR-125a-5p that limits CCR2 expression. Sci Rep. 2017;7(1):6897.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  13. Serr I, Scherm MG, Zahm AM, Schug J, Flynn VK, Hippich M, Kälin S, Becker M, Achenbach P, Nikolaev A. A miRNA181a/NFAT5 axis links impaired T cell tolerance induction with autoimmune type 1 diabetes. Sci Transl Med. 2018;10(422):eaag1782.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  14. Nabih ES, Andrawes NG. The association between circulating levels of miRNA-181a and pancreatic beta cells dysfunction via SMAD7 in type 1 diabetic children and adolescents. J Clin Lab Anal. 2016;30(5):727–31.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  15. Marchand L, Jalabert A, Meugnier E, Hende KVD, Fabien N, Nicolino M, Madec AM, Thivolet C, Rome S. miRNA-375 a sensor of glucotoxicity is altered in the serum of children with newly diagnosed type 1 diabetes. J Diabetes Res. 2016;2016:1869082.

  16. Liu JY, Yao J, Li XM, Song YC, Wang XQ, Li YJ, Yan B, Jiang Q. Pathogenic role of lncRNA-MALAT1 in endothelial cell dysfunction in diabetes mellitus. Cell Death Dis. 2014;5(10):e1506.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  17. Wang N, Zhu Y, Xie M, Wang L, Jin F, Li Y, Yuan Q, De W. Long Noncoding RNA Meg3 regulates mafa expression in mouse beta cells by inactivating Rad21, Smc3 or Sin3α. Cell Physiol Biochem. 2018;45(5):2031–43.

    Article  CAS  PubMed  Google Scholar 

  18. Krishnan P, Syed F, Jiyun Kang N, Mirmira RG, Evans-Molina C. Profiling of RNAs from human islet-derived exosomes in a model of type 1 diabetes. Int J Mol Sci. 2019;20(23):5903.

    Article  CAS  PubMed Central  Google Scholar 

  19. Li Z, Hao S, Yin H, Gao J, Yang Z. Autophagy ameliorates cognitive impairment through activation of PVT1 and apoptosis in diabetes mice. Behav Brain Res. 2016;305:265–77.

    Article  PubMed  Google Scholar 

  20. Kartha RV, Subramanian S. Competing endogenous RNAs (ceRNAs): new entrants to the intricacies of gene regulation. Front Genet. 2014;5(06):8.

    PubMed  PubMed Central  Google Scholar 

  21. Qu J, Li M, Zhong W, Hu C. Competing endogenous RNA in cancer: a new pattern of gene expression regulation. Int J Clin Exp Med. 2015;8(10):17110–6.

    CAS  PubMed  PubMed Central  Google Scholar 

  22. Yang M, Ye L, Wang B, Gao J, Liu R, Hong J, Wang W, Gu W, Ning G. Decreased miR-146 expression in peripheral blood mononuclear cells is correlated with ongoing islet autoimmunity in type 1 diabetes patients 1miR-146. J Diabetes. 2015;7(2):158–65.

    Article  CAS  PubMed  Google Scholar 

  23. Harrow J, Frankish A, Gonzalez JM, Tapanari E, Diekhans M, Kokocinski F, Aken BL, Barrell D, Zadissa A, Searle S. GENCODE: The reference human genome annotation for The ENCODE Project. Genome Res. 2012;22(9):1760–74.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  24. Ritchie ME, Belinda P, Di W, Yifang H, Law CW, Wei S, Smyth GK. limma powers differential expression analyses for RNA-sequencing and microarray studies. Nucleic Acids Res. 2015;43(7):e47.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  25. Solari A, Goeman JJ. Minimally adaptive BH: A tiny but uniform improvement of the procedure of Benjamini and Hochberg. Biometrical J Biometrische Zeitschrift. 2017;59(4):776–80.

    Article  Google Scholar 

  26. da Huang W, Sherman BT, Lempicki RA. Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat Protoc. 2009;4(1):44–57.

    Article  CAS  Google Scholar 

  27. Wencke W, Fátima SC, Mercedes R. GOplot: an R package for visually combining expression data with functional analysis. Bioinformatics. 2015;31(17):2912–4.

    Article  CAS  Google Scholar 

  28. Li L, Wang X, Xiao G, Gazdar A. Integrative gene set enrichment analysis utilizing isoform-specific expression. Genet Epidemiol. 2017;41(6):498–510.

    Article  PubMed  PubMed Central  Google Scholar 

  29. Jason M. Psych issues. JEMS. 2013;38(3):14.

    CAS  PubMed  Google Scholar 

  30. Schober P, Boer C, Schwarte LA. Correlation coefficients: appropriate use and interpretation. Anesth Analg. 2018;126(5):1763–8.

    Article  PubMed  Google Scholar 

  31. Kohl M, Wiese S, Warscheid B. Cytoscape: software for visualization and analysis of biological networks. Methods Mol Biol. 2011;696:291–303.

    Article  CAS  PubMed  Google Scholar 

  32. Yang JH, Li JH, Shao P, Zhou H, Chen YQ, Qu LH. starBase: a database for exploring microRNA-mRNA interaction maps from Argonaute CLIP-Seq and Degradome-Seq data. Nucleic Acids Res. 2011;39(Database issue):D202–9.

  33. Dweep H, Gretz N. miRWalk2.0: a comprehensive atlas of microRNA-target interactions. Nat Methods. 2015;12(8):697.

    Article  CAS  PubMed  Google Scholar 

  34. Xie H, Shi S, Chen Q, Chen Z. LncRNA TRG-AS1 promotes glioblastoma cell proliferation by competitively binding with miR-877-5p to regulate SUZ12 expression. Pathol Res Pract. 2019;215(8):152476.

    Article  CAS  PubMed  Google Scholar 

  35. Sun X, Qian Y, Wang X, Cao R, Zhang J, Chen W, Fang M. LncRNA TRG-AS1 stimulates hepatocellular carcinoma progression by sponging miR-4500 to modulate BACH1. Cancer Cell Int. 2020;20:367.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  36. He S, Wang X, Zhang J, Zhou F, Li L, Han X. TRG-AS1 is a potent driver of oncogenicity of tongue squamous cell carcinoma through microRNA-543/Yes-associated protein 1 axis regulation. Cell Cycle. 2020;19(15):1969–82.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  37. Xu H, Zhou M, Cao Y, Zhang D, Han M, Gao X, Xu B, Zhang A. Genome-wide analysis of long noncoding RNAs, microRNAs, and mRNAs forming a competing endogenous RNA network in repeated implantation failure. Gene. 2019;720:144056.

    Article  CAS  PubMed  Google Scholar 

  38. Henriksen TI, Davidsen PK, Pedersen M, Schultz HS, Hansen NS, Larsen TJ, Vaag A, Pedersen BK, Nielsen S, Scheele C. Dysregulation of a novel miR-23b/27b-p53 axis impairs muscle stem cell differentiation of humans with type 2 diabetes. Mol Metab. 2017;6(7):770–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  39. Yang W, Wang J, Chen Z, Chen J, Meng Y, Chen L, Chang Y, Geng B, Sun L, Dou L. NFE2 Induces miR-423-5p to promote gluconeogenesis and hyperglycemia by repressing hepatic FAM3A-ATP-Akt pathway. Diabetes. 2017;66(7):1819–32.

    Article  PubMed  CAS  Google Scholar 

  40. Wei G, Yi S, Yong D, Shaozhuang L, Guangyong Z, Sanyuan H. miR-320 mediates diabetes amelioration after duodenal-jejunal bypass via targeting adipoR1. Surg Obes Relat Dis. 2018;14(7):960–71.

    Article  PubMed  Google Scholar 

  41. Lapolla A, Sanzari M, Betterle C, Dalfrà MG, Masin M. Evaluation of T-cell receptor CD3+ γδ in gestational diabetes mellitus. Acta Diabetol. 2000;37(4):207–11.

    Article  CAS  PubMed  Google Scholar 

  42. Gomeztourino I, Kamra Y, Baptista R, Lorenc A, Peakman M. T cell receptor β-chains display abnormal shortening and repertoire sharing in type 1 diabetes. Nat Commun. 2017;8(1):1792.

    Article  CAS  Google Scholar 

  43. Lu G, Ota A, Ren S, Franklin S, Rau CD, Ping P, Lane TF, Zhou ZH, Reue K, Lusis AJ, et al. PPM1l encodes an inositol requiring-protein 1 (IRE1) specific phosphatase that regulates the functional outcome of the ER stress response. Mol Metab. 2013;2(4):405–16.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  44. Chen Y, Zhu J, Lum PY, Yang X, Pinto S, MacNeil DJ, Zhang C, Lamb J, Edwards S, Sieberts SK, et al. Variations in DNA elucidate molecular networks that cause disease. Nature. 2008;452(7186):429–35.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  45. Coccia EM, Cicala C, Charlesworth A, Ciccarelli C, Rossi GB, Philipson L, Sorrentino V. Regulation and expression of a growth arrest-specific gene (gas5) during growth, differentiation, and development. Mol Cell Biol. 1992;12(8):3514–21.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  46. Carter G, Miladinovic B, Patel AA, Deland L, Mastorides S, Patel NA. Circulating long noncoding RNA GAS5 levels are correlated to prevalence of type 2 diabetes mellitus. BBA Clin. 2015;4(C):102–7.

    Article  PubMed  PubMed Central  Google Scholar 

  47. Shi Y, Parag S, Patel R, Lui A, Murr M, Cai J, Patel NA. Stabilization of lncRNA GAS5 by a small molecule and its implications in diabetic adipocytes. Cell Chem Biol. 2019;26(3):319-330.e316.

    Article  CAS  PubMed  Google Scholar 

  48. Jin F, Wang N, Zhu Y, You L, Wang L, De W, Tang W. Downregulation of long noncoding RNA Gas5 affects cell cycle and insulin secretion in mouse pancreatic β cells. Cell Physiol Biochem. 2017;43(5):2062–73.

    Article  CAS  PubMed  Google Scholar 

  49. Feng B, Chakrabarti S. miR-320 regulates glucose-induced gene expression in diabetes. ISRN Endocrinol. 2012;2012:549875.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  50. Mansuy-Aubert V, Zhou QL, Xie X, Gong Z, Huang JY, Khan AR, Aubert G, Candelaria K, Thomas S, Shin DJ, et al. Imbalance between neutrophil elastase and its inhibitor α1-antitrypsin in obesity alters insulin sensitivity, inflammation, and energy expenditure. Cell Metab. 2013;17(4):534–48.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  51. He Y, Ronsein GE, Tang C, Jarvik GP, Davidson WS, Kothari V, Song HD, Segrest JP, Bornfeldt KE, Heinecke JW. Diabetes impairs cellular cholesterol efflux from ABCA1 to small HDL particles. Circ Res. 2020;127(9):1198–210.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  52. Kryvalap Y, Lo CW, Manuylova E, Baldzizhar R, Jospe N, Czyzyk J. Antibody response to serpin B13 induces adaptive changes in mouse pancreatic islets and slows down the decline in the residual beta cell function in children with recent onset in type 1 diabetes mellitus. J Biol Chem. 2016;291(1):266–78.

    Article  CAS  PubMed  Google Scholar 

  53. Yang W, Li Y, Tian T, Wang L, Lee P, Hua Q. Serum vaspin concentration in elderly patients with type 2 diabetes mellitus and macrovascular complications. BMC Endocr Disord. 2017;17(1):67.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

Download references

Acknowledgements

None.

Funding

None.

Author information

Author notes

  1. Co-first author: Chang Li and Bo Wei

Authors and Affiliations

  1. Departments of VIP Unit, China-Japan Union Hospital of Jilin University, Changchun, 130033, Jilin, China

    Chang Li

  2. Departments of Neurosurgery, China-Japan Union Hospital of Jilin University, Changchun, 130033, Jilin, China

    Bo Wei

  3. Department of Endocrinology, China-Japan Union Hospital of Jilin University, 126 Xiantai Street, Changchun, 130033, Jilin, China

    Jianyu Zhao

Authors
  1. Chang Li
    View author publications

    You can also search for this author in PubMed Google Scholar

  2. Bo Wei
    View author publications

    You can also search for this author in PubMed Google Scholar

  3. Jianyu Zhao
    View author publications

    You can also search for this author in PubMed Google Scholar

Contributions

CL and BW analyzed and interpreted the RNA-sequencing data. JYZ designed the study, collected the patients and performed the experiments. CL was a major contributor in writing the manuscript. BW and JYZ revised the manuscript. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Jianyu Zhao.

Ethics declarations

Ethics approval and consent to participate

The experiments were approved by Ethics Committee of China-Japan Union Hospital of Jilin University [No.(2020)linshen(20201127)]. Written informed consent was received from all participants. No participant under the age of 16 was involved in the validation set.

Consent for publication

Not applicable.

Competing interests

All authors declare no conflicts of interest in this paper.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Li, C., Wei, B. & Zhao, J. Competing endogenous RNA network analysis explores the key lncRNAs, miRNAs, and mRNAs in type 1 diabetes. BMC Med Genomics 14, 35 (2021). https://doi.org/10.1186/s12920-021-00877-3

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12920-021-00877-3

Keywords

BMC Medical Genomics

ISSN: 1755-8794

Contact us