- Open Access
Identification of miRNA-mRNA-TF regulatory networks in peripheral blood mononuclear cells of type 1 diabetes
BMC Endocrine Disorders volume 22, Article number: 119 (2022)
Type 1 diabetes (T1D) is a T lymphocyte-mediated and B lymphocyte-assisted autoimmune disease. We aimed to identify abnormally expressed genes in peripheral blood mononuclear cells (PBMCs) of T1D and explore their possible molecular regulatory network.
Expression datasets were downloaded from the Gene Expression Omnibus (GEO) database. Then, the differentially expressed genes (DEGs) and differentially expressed miRNAs (DEmiRNAs) were identified, and functional enrichment and immune cell infiltration analysis were performed. The starBase, miRTarBase, TarBase, JASPAR, ENCODE, and TRRUST databases constructed the miRNA-mRNA-TF regulatory network. The ROC curves were plotted to evaluate the sensitivity and specificity of miRNAs and mRNAs.
A total of 216 DEGs directly or indirectly related to type I diabetes mellitus, natural killer cell-mediated cytotoxicity, Th1, and Th2 cell differentiation, and the IL-17 and TNF signaling pathways were obtained. The miRNA-mRNA-TF network indicates that miR-320a and SOX5 are the only miRNAs and TFs that both target ADM and RRAGD. The ROC curves showed that ADM (0.9375), RRAGD (0.8958), and hsa-mir-320a (0.9417) had high accuracy in T1D diagnosis.
The constructed regulatory networks, including miR-320a/ADM/SOX5 and miR-320a/RRAGD/SOX5, may provide new insight into the mechanisms of development and progression in T1D.
Diabetes is a chronic disease. Compared with type 2 diabetes (T2D), type 1 diabetes (T1D) is an autoimmune-mediated disease in which selective destruction of pancreatic islet β cells is caused by multiple complex crosstalks of genetic and environmental determinants . Emerging evidence indicates that miRNA, as a key regulator of β-cell physiology, is essential for fine-tuning gene expression regulation in the differentiation of insulin-producing cells and contributing to the acquisition and management of unique phenotypes in diabetes . Namely, dysregulation of miRNA expression contributes to β-cell dysfunction and facilitates the development of multiple forms of diabetes mellitus.
As one of the critical biological players, transcription factors (TFs) are crucial for the pathophysiological process of multiple diseases, including cancer, autoimmunity, and diabetes . Given that the therapeutic regulation of TFs by drugs may affect significant changes in gene expression patterns and could improve disease conditions, identifying the key TFs in the pathophysiology of T1D may promise better antidiabetic drugs. In addition, early research has shown that disorders of components in peripheral blood mononuclear cells (PBMCs), including T cells and B cells, and any factors that can affect their functional balance may be involved in the occurrence and development of T1D .
Although the importance of miRNAs in diabetes has been slowly recognized in recent years [2, 5, 6], little is known about the complex molecular regulation mechanisms of miRNAs, mRNAs, and TFs in the development of T1D. Therefore, it is crucial to explore the molecular regulatory mechanisms of miRNAs, mRNAs, and TFs in the PBMCs of patients with T1D, which may be valuable for diagnosis and treatment.
Materials and methods
The expression datasets of GSE55098 , GSE55099, and GSE33440  were downloaded from the Gene Expression Omnibus (GEO) database (website: http://www.ncbi.nlm.nih.gov/geo/) . The three datasets contain a total of 40 newly diagnosed T1D patients and 26 normal controls. The details of GSE55098, GSE55099, and GSE33440 are listed in Table 1. In this study, the GSE55098 dataset was utilized for DEG analysis, functional enrichment analysis, ssGSEA, and receiver operating characteristic (ROC) curve analysis; the GSE55099 dataset was mainly used for differential microRNA expression analysis. The DEGs of GSE33440 were utilized only to take the intersection with the DEGs of GSE55099 to obtain the hub genes.
Based on the annotation file of the platform, the corresponding genes with the probes were further matched, and probes that did not map to genes were removed. The expression of the same gene in multiple probes is expressed as the mean value. The differentially expressed genes (DEGs) and differentially expressed miRNAs (DEmiRNAs) between the peripheral blood mononuclear cells of T1D and non-T1D patients under the criteria of |log fold change|> 0.5 and P < 0.05 were obtained by using the “limma”  package.
Functional enrichment and immune cell infiltration analysis
To further investigate the major biological functions of DEGs, the Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways  were explored in GSE55098 based on the "clusterProfiler"  package with the criteria of q-value < 0.05. Moreover, the single-sample gene set enrichment analysis (ssGSEA) algorithm in the “GSVA” package  was utilized to quantify the relative abundance of 24 immune cell types  infiltration within the entire PBMC subset of profile GSE55098.
Hub gene selected and evaluation
In this study, intersecting genes between GSE33440 and GSE55098 were defined as hub genes. Then, we utilized Pearson correlation and Wilcoxon rank-sum tests to assess the association of immune cell infiltration with the expression of these hub genes. Moreover, receiver operating characteristic (ROC) curves and related area under the curve (AUC) values were performed to evaluate the sensitivity and specificity of hub genes and miRNAs for T1D diagnosis. The STATS package calculated the correlation between seven hub genes.
Construction of the miRNA-mRNA-TF target regulatory network
Based on the current miRNA-mRNA interaction regulatory network prediction database, starBase (Version 2.0), miRTarBase (Version 8.0), and TarBase (Version 8.0) were used to predict the potential upstream regulatory miRNAs of hub genes [14,15,16]. Moreover, the JASPAR, ENCODE, and TRRUST [17,18,19] databases were utilized to predict the possible transcription factors of seven hub genes. The intersecting TFs that appeared in two transcription factor prediction databases were selected as the targeted TFs, and the targeted miRNAs were obtained in the same way. Based on starBase, miRTarBase, TarBase, JASPAR, ENCODE, TRRUST, and DEmiRNAs, we further explored the effect of the miRNAs and transcription factors (TFs) that targeted the hub genes related to T1D. Visualization of the miRNA-mRNA-TF target regulatory network was performed by Cytoscape software .
Identification of DEGs
Based on the criteria of |log FC|> 0.5 and P < 0.05 of the data preprocessing specified in the Materials and Methods section, we identified 21 DEmiRNAs (10 up- and 11 downregulated miRNAs), 216 DEGs (99 up- and 117 downregulated genes) and 545 DEGs (296 up- and 249 downregulated genes) in GSE55098, GSE55099 and GSE33440, respectively (Supplementary Table 1). Moreover, we visualized the top 10 up-and downregulated DEmiRNAs and DEGs in related datasets by volcano mapping and heat mapping (Fig. 1).
Functional enrichment and immune cells infiltration
To further explore the functional implications, we performed GO and KEGG functional enrichment analysis of DEGs in the GSE55098 dataset (Supplementary Table 2). The GO results indicated that these DEGs were mainly associated with neutrophil activation involved in the immune response, natural killer cell-mediated immunity regulation, and lymphocyte-mediated immunity in the biological process analysis (Fig. 2B). KEGG results showed that these DEGs were directly or indirectly related to type I diabetes mellitus, natural killer cell-mediated cytotoxicity, Th1, and Th2 cell differentiation, cytokine-cytokine receptor interaction, Th17 cell differentiation, the IL-17 signaling pathway, and the TNF signaling pathway (Fig. 2C). Moreover, we noticed that there was a significant difference in the expression of multiple immune cell types between the T1D and no-T1D groups (Fig. 2A), including activated CD8 T cells, immature dendritic cells, mast cells, regulatory T cells, type 1 T helper cells, and type 17 T helper cells (P < 0.001).
Hub gene evaluation
By taking the intersecting genes of the DEGs in GSE33440 and GSE55098 as the target genes of this study, we finally obtained seven upregulated DEGs as hub genes, including ADM, FFAR2, HLA-DQA1, ID3, OGFRL1, RRAGD, and TNF (Fig. 3H). ssGSEA is a deconvolution algorithm that assesses the level of immune cell infiltration in samples based on the expression level of immune cell-specific marker genes . Based on this algorithm, we found that these hub genes are mainly related to the expression of multiple immune cell types, including neutrophils, type 2 T helper cells, macrophages, regulatory T cells, gamma delta T cells, activated CD8 T cells, mast cells, immature dendritic cells, eosinophils and type 17 T helper cells (Fig. 3A-G). ROC curves and AUC values were plotted to evaluate the sensitivity and specificity of hub genes for T1D diagnosis. The results indicated that these hub genes were highly accurate in T1D diagnosis (Fig. 4). The AUCs were 0.9375 (ADM), 0.9792 (TNF), 0.9688 (OGFRL1), 0.9271 (HLA-DQA1), 0.8958 (RRAGD), 0.8250 (ID3), and 0.7604 (FFAR2). Similarly, miR-320a also has high sensitivity and specificity for T1D diagnosis, and the AUC was 0.9417 (Fig. 4H).
miRNA—mRNA-TF target regulatory network
To further explore the regulatory relationship in hub genes, TFs, and miRNAs, we constructed a miRNA-mRNA-TF target regulatory network in this study and then visualized the network with Cytoscape software. Among them, 25 miRNAs targeted ADM, 33 miRNAs targeted RRAGD, 26 miRNAs targeted OGFRL1, 15 miRNAs targeted ID3, 2 miRNAs targeted FFAR2, and 6 miRNAs targeted HLA-DQA2, and 12 miRNAs targeted TNF. miR-320a, as a DEmiRNA, appeared in the miRNA prediction database (Fig. 5A; Supplement Table 3). A total of 27 TFs targeted seven hub genes in the TRRUST database, 66 TFs in the JASPAR database, and 139 TFs in the ENCODE database; among them, 19 TFs appeared in more than two databases (Fig. 5B; Supplement Table 4). Remarkably, miR-320a and SOX5 are the only miRNAs and TFs targeting ADM and RRAGD. Moreover, there was a high correlation between the ADM and RRAGD (Fig. 5C).
This study analyzed a PBMC microarray dataset involving T1D patients and normal controls and obtained 216 DEGs. The functional enrichment analysis results indicate that these DEGs are closely related to various immunity regulation pathways and biological processes. Since T1D is a T lymphocyte-mediated and B lymphocyte-assisted autoimmune disease , it is not surprising that these DEGs are closely associated with immune processes. Then, we obtained seven hub genes: ADM, FFAR2, HLA-DQA1, ID3, OGFRL1, RRAGD, and TNF. These hub genes are not only associated with the expression of multiple immune cell types but also have a high accuracy in T1D diagnosis, with AUCs of 0.9375 (ADM), 0.9792 (TNF), 0.9688 (OGFRL1), 0.9271 (HLA-DQA1), 0.8958 (RRAGD), 0.8250 (ID3), and 0.7604 (FFAR2). Finally, based on the miRNA prediction databases and transcription factor prediction databases, we constructed a miRNA-mRNA-TF target regulatory network that included 49 miRNAs, 7 mRNAs, and 19 TFs. Remarkably, miR-320a and SOX5 are the only miRNAs and TFs targeting ADM and RRAGD.
Type 1 diabetes (T1D) is a complex autoimmune disease in which the development of pathology requires the crosstalk between various immune cells . As one of the essential innate immune cells, Neutrophils have been required for initiating diabetogenic T cell responses in non-obese diabetic (NOD) mice . Neutrophils infiltrated the pancreas and could form neutrophil extracellular traps that induce dendritic cell activation leading to Th1 polarization of T1D patients  and escalating the activating macrophages-related inflammation through modulating inflammatory cytokines directly or indirectly . Interestingly, Macrophages can stimulate β-cells by secreting IL-1β to produce CXCR2 ligands, which are contributed to neutrophil recruitment in pancreatic islets . However, some questions remain: Are there some trinsic mechanisms that can different immune subsets recruiting (such as Neutrophils and Macrophages), affect their pathogenic properties, prolong and enhance their inflammatory activity, and perhaps promote β-cell destruction?
This study found that ADM and RRAGD in the regulatory network are closely associated with Macrophages and/or Neutrophils (P < 0.01). Early studies reported that ADM is elevated in the plasma of patients with pancreatic dysfunctions, including T1D . Zudaire et al.  further indicated that elevated ADM increases circulating glucose levels, while a blocking monoclonal antibody of ADM has the opposite effect and ameliorates postprandial recovery. Given the localization of ADM in the pancreas and experimental evidence, we hypothesize that ADM may play an important role in maintaining insulin homeostasis and normoglycemia and altering pancreatic physiologies. Although no evidence directly indicates that ADM mediates pancreatic β-cell damage by macrophages, ADM production is related to increased expression of monocyte/macrophage differentiation . The study by Xu et al. shows that ADM can enhance the migration and invasion of myelomonocytic cells by activating multiple signaling pathways, including MAPK, PI3K/Akt, and eNOS in pancreatic cancer, as well as promote myelomonocytic cells trans-endothelial migration through increased expression of VCAM-1 and ICAM-1. Conversely, knockout or antagonist of ADM significantly inhibited monocyte recruitment . As for neutrophils, ADM can exert pro-inflammatory effects in the cutaneous microvasculature to mediate vasodilation and promote edema formation related to acute inflammatory processes, thereby enhancing neutrophil accumulation in the skin microvasculature . Another hub gene, RRAGD, is a monomeric guanine nucleotide-binding protein that can act as a molecular switch to interact with the C terminal region of leucyl tRNA synthase, which acts as a direct sensor of the amino acid leucine and is contributed to the activation of mTORC1 . Studies have shown that mTORC1 is closely related to glucose metabolism, lipid metabolism, and islet resistance . In the state of excess nutrients, the activation of mTORC1 signaling may promote oxidative stress, ER stress, and inflammation, thereby damaging pancreatic β cells . The activation of mTORC1 signaling may also impair autophagic flux, consequently triggering the death of pancreatic β-cells  and leading to T2D pancreatic islet failure or T1D. Namely, RRAGD may participate in the occurrence and development of diabetes by activating mTOR1. Moreover, activated mTORC1 regulates inflammatory immune responses in various innate immune cells ( such as monocytes, macrophages, and dendritic cells) in the . Manipulation of genes in its pathway has been shown to alter macrophage polarization and the production of inflammatory and immunoregulatory cytokines in mice . Specifically, mTORC1 inhibits apoptosis and induces macrophage proliferation by inducing the expression and metabolic reprogramming of the cell cycle kinase CDK4 .
In further investigation of regulatory mechanisms of hub genes in T1D, we found that miR-320a and SOX5 are the only miRNA and TF that both target ADM and RRAGD. Knowledgeably, abnormal miRNA expression can cause some pathological conditions. Previous research indicated that miR-320 is lowly expressed in the serum of people with diabetes , especially those with diabetic retinopathy . Similarly, decreased expression of miR-320 was observed in the kidneys of glucose-induced diabetic rats and human umbilical vein endothelial cells . Notably, miR-320a is significantly downregulated in the human islets of T2D subjects, which is accompanied by an increase in glucagon, and the increased glucagon as a marker of diabetes will worsen hyperglycemia . Moreover, the downregulation of miR-320a participates in regulating multiple signaling pathways, including PI3K/AKT and MAPK/ERK, by targeting IGF-1R, which is associated with obesity and T2D . However, there are inconsistencies in the results of some studies both on miR-320a and diabetes. Ling et al. revealed that miR-320 is upregulated in insulin-resistant adipocytes, and anti-miR-320 oligonucleotides can reverse insulin resistance and increase adipocyte insulin sensitivity by insulin PI3K signaling pathways . Given the inconsistent role of miR-320a in diabetes, a future longitudinal study of miR-320 and diabetes-related progression is needed to understand the specific expression and molecular regulatory mechanisms.
SOX5 is a member of the SOX transcription factor family that participates in embryonic development and cell fate regulation and is expressed in both the nucleus and cytoplasm of human α- cells and β- cells . Cytosolic expression of SOX5 was especially evident in α-cells; conversely, nuclear SOX5 was reduced by 67% in T2D β-cells compared with nondiabetic β-cells. Moreover, Axelsson et al. further indicated that knockdown of SOX5 attenuates glucose-stimulated insulin secretion and reduces the expression of L-type Ca2 + channels. In contrast, overexpressed SOX5 can restore impaired glucose-stimulated insulin secretion and increase insulin secretion by 18% in T2D islets .
Although understanding the regulatory network of miRNAs, TFs, and mRNA is crucial for revealing the underlying mechanisms of T1D and identifying and constructing a molecular regulatory network, some limitations remain in this study. This study is only based on the GEO database to conduct preliminary exploration at the level of bioinformatics, and further comprehensive experiments in vitro and in vivo are needed to verify these constructed regulatory mechanisms. Therefore, obtaining insight into the gene interactions and molecular mechanisms of diabetes still has a long way to go.
Generally, the constructed regulatory mechanism networks, including mir-320a/ADM/SOX5 and mir-320a/RRAGD/SOX5, are robust and promising and may provide new insight into the mechanisms of development and progression of T1D.
Availability of data and materials
The entire sequencing profile data in this study come from the Gene Expression Omnibus (GEO) database (website: http://www.ncbi.nlm.nih.gov/geo/).
Atkinson MA, Eisenbarth GS, Michels AW. Type 1 diabetes. Lancet. 2014;383(9911):69–82.
Brozzi F. Role of microRNA in pancreatic beta cell function. Int Rev Cell Mol Biol. 2021;359:257–86.
Henley MJ, Koehler AN. Advances in targeting “undruggable” transcription factors with small molecules. Nat Rev Drug Discov. 2021;20(9):669–88.
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.
Januszewski AS, Cho YH, Joglekar MV, Farr RJ, Scott ES, Wong WKM, Carroll LM, Loh YW, Benitez-Aguirre PZ, Keech AC, et al. Insulin micro-secretion in Type 1 diabetes and related microRNA profiles. Sci Rep. 2021;11(1):11727.
Tamara A, Coulson DJ, Latief JS, Bakhashab S, Weaver JU. Upregulated anti-angiogenic miR-424-5p in type 1 diabetes (model of subclinical cardiovascular disease) correlates with endothelial progenitor cells, CXCR1/2 and other parameters of vascular health. Stem Cell Res Ther. 2021;12(1):249.
Irvine KM, Gallego P, An X, Best SE, Thomas G, Wells C, Harris M, Cotterill A, Thomas R. Peripheral blood monocyte gene expression profile clinically stratifies patients with recent-onset type 1 diabetes. Diabetes. 2012;61(5):1281–90.
Barrett T, Wilhite SE, Ledoux P, Evangelista C, Kim IF, Tomashevsky M, Marshall KA, Phillippy KH, Sherman PM, Holko M, et al. NCBI GEO: archive for functional genomics data sets–update. Nucleic Acids Res. 2013;41(Database issue):D991-995.
Ritchie ME, Phipson B, Wu D, Hu Y, Law CW, Shi W, Smyth GK. limma powers differential expression analyses for RNA-sequencing and microarray studies. Nucleic Acids Res. 2015;43(7):e47.
Subramanian A, Tamayo P, Mootha VK, Mukherjee S, Ebert BL, Gillette MA, Paulovich A, Pomeroy SL, Golub TR, Lander ES, et al. Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles. Proc Natl Acad Sci U S A. 2005;102(43):15545–50.
Yu G, Wang LG, Han Y, He QY. clusterProfiler: an R package for comparing biological themes among gene clusters. OMICS. 2012;16(5):284–7.
Hänzelmann S, Castelo R, Guinney J. GSVA: gene set variation analysis for microarray and RNA-seq data. BMC Bioinformatics. 2013;14:7.
Bindea G, Mlecnik B, Tosolini M, Kirilovsky A, Waldner M, Obenauf AC, Angell H, Fredriksen T, Lafontaine L, Berger A, et al. Spatiotemporal dynamics of intratumoral immune cells reveal the immune landscape in human cancer. Immunity. 2013;39(4):782–95.
Li JH, Liu S, Zhou H, Qu LH, Yang JH. starBase v2.0: decoding miRNA-ceRNA, miRNA-ncRNA and protein-RNA interaction networks from large-scale CLIP-Seq data. Nucleic Acids Res. 2014;42(Database issue):D92-97.
Huang HY, Lin YC, Li J, Huang KY, Shrestha S, Hong HC, Tang Y, Chen YG, Jin CN, Yu Y, et al. miRTarBase 2020: updates to the experimentally validated microRNA-target interaction database. Nucleic Acids Res. 2020;48(D1):D148-d154.
Karagkouni D, Paraskevopoulou MD, Chatzopoulos S, Vlachos IS, Tastsoglou S, Kanellos I, Papadimitriou D, Kavakiotis I, Maniou S, Skoufos G, et al. DIANA-TarBase v8: a decade-long collection of experimentally supported miRNA-gene interactions. Nucleic Acids Res. 2018;46(D1):D239-d245.
Fornes O, Castro-Mondragon JA, Khan A, van der Lee R, Zhang X, Richmond PA, Modi BP, Correard S, Gheorghe M, Baranašić D, et al. JASPAR 2020: update of the open-access database of transcription factor binding profiles. Nucleic Acids Res. 2020;48(D1):D87-d92.
Moore JE, Purcaro MJ, Pratt HE, Epstein CB, Shoresh N, Adrian J, Kawli T, Davis CA, Dobin A, Kaul R, et al. Expanded encyclopaedias of DNA elements in the human and mouse genomes. Nature. 2020;583(7818):699–710.
Han H, Cho JW, Lee S, Yun A, Kim H, Bae D, Yang S, Kim CY, Lee M, Kim E, et al. TRRUST v2: an expanded reference database of human and mouse transcriptional regulatory interactions. Nucleic Acids Res. 2018;46(D1):D380-d386.
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(11):2498–504.
Kitamura K, Kangawa K, Kawamoto M, Ichiki Y, Nakamura S, Matsuo H, Eto T. Adrenomedullin: a novel hypotensive peptide isolated from human pheochromocytoma. Biochem Biophys Res Commun. 1993;192(2):553–60.
Lehuen A, Diana J, Zaccone P, Cooke A. Immune cell crosstalk in type 1 diabetes. Nat Rev Immunol. 2010;10(7):501–13.
Diana J, Simoni Y, Furio L, Beaudoin L, Agerberth B, Barrat F, Lehuen A. Crosstalk between neutrophils, B-1a cells and plasmacytoid dendritic cells initiates autoimmune diabetes. Nat Med. 2013;19(1):65–73.
Parackova Z, Zentsova I, Vrabcova P, Klocperk A, Sumnik Z, Pruhova S, Petruzelkova L, Hasler R, Sediva A. Neutrophil extracellular trap induced dendritic cell activation leads to Th1 polarization in type 1 diabetes. Front Immunol. 2020;11:661.
Warnatsch A, Ioannou M, Wang Q, Papayannopoulos V, Inflammation. Neutrophil extracellular traps license macrophages for cytokine production in atherosclerosis. Science. 2015;349(6245):316–20.
Diana J, Lehuen A. Macrophages and β-cells are responsible for CXCR2-mediated neutrophil infiltration of the pancreas during autoimmune diabetes. EMBO Mol Med. 2014;6(8):1090–104.
Hayashi M, Shimosawa T, Isaka M, Yamada S, Fujita R, Fujita T. Plasma adrenomedullin in diabetes. Lancet. 1997;350(9089):1449–50.
Zudaire E, Cuttitta F, Martínez A. Regulation of pancreatic physiology by adrenomedullin and its binding protein. Regul Pept. 2003;112(1–3):121–30.
Kubo A, Minamino N, Isumi Y, Kangawa K, Dohi K, Matsuo H. Adrenomedullin production is correlated with differentiation in human leukemia cell lines and peripheral blood monocytes. FEBS Lett. 1998;426(2):233–7.
Xu M, Qi F, Zhang S, Ma X, Wang S, Wang C, Fu Y, Luo Y. Adrenomedullin promotes the growth of pancreatic ductal adenocarcinoma through recruitment of myelomonocytic cells. Oncotarget. 2016;7(34):55043–56.
Chu DQ, Choy M, Foster P, Cao T, Brain SD. A comparative study of the ability of calcitonin gene-related peptide and adrenomedullin(13–52) to modulate microvascular but not thermal hyperalgesia responses. Br J Pharmacol. 2000;130(7):1589–96.
Sancak Y, Bar-Peled L, Zoncu R, Markhard AL, Nada S, Sabatini DM. Ragulator-Rag complex targets mTORC1 to the lysosomal surface and is necessary for its activation by amino acids. Cell. 2010;141(2):290–303.
Szwed A, Kim E, Jacinto E. Regulation and metabolic functions of mTORC1 and mTORC2. Physiol Rev. 2021;101(3):1371–426.
Ardestani A, Maedler K. Deathly triangle for pancreatic β-cells: Hippo pathway-MTORC1-autophagy. Autophagy. 2021;17(12):4494–6.
Weichhart T, Hengstschläger M, Linke M. Regulation of innate immune cell function by mTOR. Nat Rev Immunol. 2015;15(10):599–614.
Byles V, Covarrubias AJ, Ben-Sahra I, Lamming DW, Sabatini DM, Manning BD, Horng T. The TSC-mTOR pathway regulates macrophage polarization. Nat Commun. 2013;4:2834.
Linke M, Pham HT, Katholnig K, Schnöller T, Miller A, Demel F, Schütz B, Rosner M, Kovacic B, Sukhbaatar N, et al. Chronic signaling via the metabolic checkpoint kinase mTORC1 induces macrophage granuloma formation and marks sarcoidosis progression. Nat Immunol. 2017;18(3):293–302.
Zampetaki A, Kiechl S, Drozdov I, Willeit P, Mayr U, Prokopi M, Mayr A, Weger S, Oberhollenzer F, Bonora E, et al. Plasma microRNA profiling reveals loss of endothelial miR-126 and other microRNAs in type 2 diabetes. Circ Res. 2010;107(6):810–7.
Prado MSG, de Jesus ML, de Goes TC, Mendonça LSO, Kaneto CM. Downregulation of circulating miR-320a and target gene prediction in patients with diabetic retinopathy. BMC Res Notes. 2020;13(1):155.
Feng B, Chakrabarti S. miR-320 regulates glucose-induced gene expression in diabetes. ISRN Endocrinol. 2012;2012:549875.
Jo S, Xu G, Jing G, Chen J, Shalev A. Human glucagon expression is under the control of miR-320a. Endocrinology. 2021;162(3):bqaa238.
Guo T, Feng Y, Liu Q, Yang X, Jiang T, Chen Y, Zhang Q. MicroRNA-320a suppresses in GBM patients and modulates glioma cell functions by targeting IGF-1R. Tumour Biol. 2014;35(11):11269–75.
Ling HY, Ou HS, Feng SD, Zhang XY, Tuo QH, Chen LX, Zhu BY, Gao ZP, Tang CK, Yin WD, et al. CHANGES IN microRNA (miR) profile and effects of miR-320 in insulin-resistant 3T3-L1 adipocytes. Clin Exp Pharmacol Physiol. 2009;36(9):e32-39.
Axelsson AS, Mahdi T, Nenonen HA, Singh T, Hänzelmann S, Wendt A, Bagge A, Reinbothe TM, Millstein J, Yang X, et al. Sox5 regulates beta-cell phenotype and is reduced in type 2 diabetes. Nat Commun. 2017;8:15652.
This research was supported by a grant from The Youth Fund Project of the First Affiliated Hospital of Shihezi University School of Medicine (QN202104).
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Wang, W., Wang, H., Liu, Y. et al. Identification of miRNA-mRNA-TF regulatory networks in peripheral blood mononuclear cells of type 1 diabetes. BMC Endocr Disord 22, 119 (2022). https://doi.org/10.1186/s12902-022-01038-y
- Type 1 diabetes (T1D)
- Differentially expressed genes (DEGs)
- Transcription factors (TFs)
- Regulatory networks