- Research article
- Open Access
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)
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.
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.
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.
TRG-AS1—miR-23b—PPM1L and GAS5—miR-320a—SERPINA1 regulatory axises might impact the pathogenesis of T1D.
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 .An annual 3%—4% increase in the incidence of T1D in childhood was estimated in several developed countries . Besides, T1D decreases approximately 11–13 years of life expectancy in developed countries, and even more time in developing countries . 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 . 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 . 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 . 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 . 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.
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. . 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.
Using the Release 26 (GRCh38.p10) reference genome in GENCODE database (https://www.gencodegenes.org/human/) , 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)  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 . 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/) , and the corresponding results were visualized by the R package GOplot (http://cran.r-project.org/web/packages/GOplot) . Using gene set enrichment analysis (GSEA) , 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 , the Pearson correlation coefficients  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) .
Using starbase database (version 3.0, http://starbase.sysu.edu.cn/) , 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/) . Finally, a lncRNA-miRNA-mRNA interaction network was visualized by Cytoscape software  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 .
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.
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.
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).
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).
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).
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.
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).
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).
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 , hepatocellular carcinoma  and tongue squamous cell carcinoma . It might act as a ceRNA to regulate miRNA-543-Yes-associated protein 1 in tongue squamous cell carcinoma , miR-4500-BACH1 in hepatocellular carcinoma  and miR-877-5p-SUZ12 in glioblastoma . It also identified to be involved in repeated implantation failure . 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 . 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 . 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 . 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 . The macrophage-enriched network has a causal correlation with metabolic disease traits, which involved three obesity genes (including PPM1L) . 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 . 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 . Shi et al. demonstrated that GAS5 regulates insulin signaling in adipocytes, and suggested it might be a potential target for T2DM . 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 . 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 . Wei et al. found miR-320 mimic correlated with impaired gluconeogenesis and lipid metabolism by regulating adipoR1 . 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 . Lower levels of SERPINA1 selectively impaired the ATP-binding cassette transporter A1 cholesterol efflux capacity in T2D . 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 . 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 . Thus, GAS5—miR-320a/miR-23b/miR-423—SERPINA1 regulatory axis might also function in the pathogenesis of T1D.
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/).
Type 1 diabetes
Long non-coding RNA
Differentially expressed mRNAs
Type 2 diabetes
Regulatory T cells
Nuclear factor of activated T cells 5
Metastasis Associated Lung Adenocarcinoma Transcript 1
Competing endogenous RNAs
Gene Expression Omnibus
Peripheral blood mononuclear cell
Gene set enrichment analysis
Kyoto Encyclopedia of Genes and Genomes
Quantitatively real-time PCR
Normalized enrichment score
T cell receptor gamma locus antisense RNA 1
Protein Phosphatase: Mg2 + /Mn2 + Dependent 1 l
Growth arrest specific 5
Serpin family A member 1
Adiponectin Receptor 1
Puchulu FM. Definition, diagnosis and classification of diabetes mellitus. Diabetes Care. 2018;33(Suppl 1):S62–9.
El-Sappagh S, Elmogy M, Ali F, Kwak KS. A case-base fuzzification process: diabetes diagnosis case study. Soft Comput. 2018;9:5815–34.
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.
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.
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.
Rossi G, Association AD. Diagnosis and classification of diabetes mellitus. Diabetes Care. 2018;33(Suppl 1):S62–9.
Classification of diabetes mellitus. https://www.who.int/publications/i/item/classification-of-diabetes-mellitus. Accessed 2 Dec 2020.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
Kartha RV, Subramanian S. Competing endogenous RNAs (ceRNAs): new entrants to the intricacies of gene regulation. Front Genet. 2014;5(06):8.
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.
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.
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.
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.
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.
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.
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.
Li L, Wang X, Xiao G, Gazdar A. Integrative gene set enrichment analysis utilizing isoform-specific expression. Genet Epidemiol. 2017;41(6):498–510.
Jason M. Psych issues. JEMS. 2013;38(3):14.
Schober P, Boer C, Schwarte LA. Correlation coefficients: appropriate use and interpretation. Anesth Analg. 2018;126(5):1763–8.
Kohl M, Wiese S, Warscheid B. Cytoscape: software for visualization and analysis of biological networks. Methods Mol Biol. 2011;696:291–303.
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.
Dweep H, Gretz N. miRWalk2.0: a comprehensive atlas of microRNA-target interactions. Nat Methods. 2015;12(8):697.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
Feng B, Chakrabarti S. miR-320 regulates glucose-induced gene expression in diabetes. ISRN Endocrinol. 2012;2012:549875.
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.
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.
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.
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.
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
All authors declare no conflicts of interest in this paper.
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.
About this article
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
- Type 1 diabetes
- Long non-coding RNA
- Differential expression analysis
- Interaction network