
==== Front
JMIR Bioinform Biotech
JMIR Bioinform Biotech
JBB
JMIR Bioinformatics and Biotechnology
2563-3570
JMIR Publications Toronto, Canada

v3i1e32437
10.2196/32437
Original Paper
Original Paper
Novel Molecular Networks and Regulatory MicroRNAs in Type 2 Diabetes Mellitus: Multiomics Integration and Interactomics Study
Mavragani Amaryllis
Chrisman Brianna Sierra
Giri Mamta
Hetti Arachchilage Madara
Khokhar Manoj MSc 1https://orcid.org/0000-0001-7106-8545

Roy Dipayan MD 1https://orcid.org/0000-0002-3429-1470

Tomo Sojit MD 1https://orcid.org/0000-0001-5381-4348

Gadwal Ashita MSc 1https://orcid.org/0000-0002-6049-8635

Sharma Praveen PhD 1https://orcid.org/0000-0002-8324-737X

Purohit Purvi PhD https://orcid.org/0000-0001-8559-2911
1Department of Biochemistry All India Institute of Medical Sciences CRL-2, first floor, Basni Industrial Area MIA 2nd Phase, Basni Jodhpur, 342005 India 91 9928388223 dr.purvipurohit@gmail.com

1 Department of Biochemistry All India Institute of Medical Sciences Jodhpur India
Corresponding Author: Purvi Purohit dr.purvipurohit@gmail.com
Jan-Dec 2022
23 2 2022
3 1 e3243728 7 2021
23 9 2021
18 11 2021
27 12 2021
©Manoj Khokhar, Dipayan Roy, Sojit Tomo, Ashita Gadwal, Praveen Sharma, Purvi Purohit. Originally published in JMIR Bioinformatics and Biotechnology (https://bioinform.jmir.org), 23.02.2022.
2022
https://creativecommons.org/licenses/by/4.0/ This is an open-access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work, first published in JMIR Bioinformatics and Biotechnology, is properly cited. The complete bibliographic information, a link to the original publication on https://bioinform.jmir.org/, as well as this copyright and license information must be included.

Background

Type 2 diabetes mellitus (T2DM) is a metabolic disorder with severe comorbidities. A multiomics approach can facilitate the identification of novel therapeutic targets and biomarkers with proper validation of potential microRNA (miRNA) interactions.

Objective

The aim of this study was to identify significant differentially expressed common target genes in various tissues and their regulating miRNAs from publicly available Gene Expression Omnibus (GEO) data sets of patients with T2DM using in silico analysis.

Methods

Using differentially expressed genes (DEGs) identified from 5 publicly available T2DM data sets, we performed functional enrichment, coexpression, and network analyses to identify pathways, protein-protein interactions, and miRNA-mRNA interactions involved in T2DM.

Results

We extracted 2852, 8631, 5501, 3662, and 3753 DEGs from the expression profiles of GEO data sets GSE38642, GSE25724, GSE20966, GSE26887, and GSE23343, respectively. DEG analysis showed that 16 common genes were enriched in insulin secretion, endocrine resistance, and other T2DM-related pathways. Four DEGs, MAML3, EEF1D, NRG1, and CDK5RAP2, were important in the cluster network regulated by commonly targeted miRNAs (hsa-let-7b-5p, hsa-mir-155-5p, hsa-mir-124-3p, hsa-mir-1-3p), which are involved in the advanced glycation end products (AGE)-receptor for advanced glycation end products (RAGE) signaling pathway, culminating in diabetic complications and endocrine resistance.

Conclusions

This study identified tissue-specific DEGs in T2DM, especially pertaining to the heart, liver, and pancreas. We identified a total of 16 common DEGs and the top four common targeting miRNAs (hsa-let-7b-5p, hsa-miR-124-3p, hsa-miR-1-3p, and has-miR-155-5p). The miRNAs identified are involved in regulating various pathways, including the phosphatidylinositol-3-kinase-protein kinase B, endocrine resistance, and AGE-RAGE signaling pathways.

type 2 diabetes mellitus
interactomics
integrative genomics
protein-protein interaction
microRNAs
miRNA
bioinformatics
multiomics
genomics
gene expression
==== Body
pmcIntroduction

Interactions among DNA, RNA, and proteins regulate their functions and have an immense effect on the underlying mechanistic processes in the pathophysiology of many diseases. Owing to the advent of newer technologies such as microarray and genome sequencing, it is now possible to investigate and analyze an enormous amount of genomic and proteomic data to predict disease pathology, outcome, and possible therapeutic targets [1]. Diabetes is a metabolic disorder characterized by hyperglycemia and glycosuria, which, if left untreated, leads to an array of complications and associated comorbidities [2]. These can include obesity, cardiomyopathy, nephropathy, retinopathy, neuropathy, and peripheral vascular disease, which have a lasting adverse effect on the quality of the patient’s life. To date, diabetes has affected almost half a billion individuals worldwide [3]. The absence of effective treatment strategies for this disease makes it a challenge to manage. The obligatory lifestyle changes and multiple treatment modalities, along with lifelong disease monitoring, depict an urgent and unmet need to develop newer and specific preventive and treatment strategies. Mortality rates in patients with type 2 diabetes mellitus (T2DM) are higher than those of individuals without diabetes and are linked to increased cardiovascular, renovascular, and neuropathic risks [4,5]. Thus, to reduce the morbidity and mortality associated with T2DM, it is important to gain a better understanding of its pathogenic pathways and regulation mechanisms. There is accumulating evidence that microRNAs (miRNAs) play an essential role in diabetes by reducing the expression of their various target genes [6,7]. It is also crucial to select the right target for disease treatment strategies in the early discovery phases, thus maximizing the drug’s success rates in the latter phases [8].

Currently, there is a vast amount of genomic data on diabetes and its complications. However, from its detection to the management of its late-stage complications, many areas still need to be explored and lacunae need to be filled. The role of molecular integration networks regulating the pathogenesis of T2DM in specific tissues is unknown. In this study, we have undertaken an in silico approach with existing tissue-specific microarray data of patients with diabetes to address this particular area by detecting novel diabetes-associated genes, their regulatory miRNAs, and their interactions to predict the pivotal pathways in tissues that are associated with disease onset and progression.

We selected five data sets from the Gene Expression Omnibus (GEO) database comprising the expression profiles of patients with diabetes and corresponding controls, and identified 16 differentially expressed genes (DEGs) overlapping the three preassigned groups. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses were further used to classify the DEGs into cellular component (CC), biological process (BP), and molecular function (MF) classes. We selected four clusters from the protein-protein interaction (PPI) network and identified the seed genes. We further investigated the miRNA and hub gene network. Finally, we explored the 16 hub genes for biological pathway enrichment and their targeting miRNAs.

Methods

Data Collection

We searched several keywords, including “type 2 diabetes mellitus,” “tissue,” “pancreas,” “liver,” “heart,” “expression profiling by array,” and “Homo sapiens,” in the GEO data sets, among which five were selected for this study: GSE38642 [9-11], GSE25724 [12], GSE20966 [13], GSE26887 [14], and GSE23343 [15] (Table 1).

Table 1 Description of Gene Expression Omnibus data sets for three groups of organs.

Sample organ	T2DMa	Control	Platform	Country	Year	

	Samples, n	Sex (M/F)	Mean age (years)	Samples, n	Sex (M/F)	Mean age (years)	
	
	
	
Pancreas	

	GSE25724 [12]	6	3/3	58.1	7	4/3	70.5	Affymetrix Human Genome U133A Array	Italy	2010	

	GSE20966 [13]	10	7/3	60.3	10	6/4	67.3	Affymetrix Human X3P Array	United States	2010	

	GSE38642 [9-11]	9	5/4	57.0	54	31/23	56.6	Affymetrix Human Gene 1.0 ST Array	Sweden	2012	
Heart (GSE26887) [14]	7	6/1	65.1	5	2/3	48.4	Affymetrix Human Gene 1.0 ST Array	Italy	2012	
Liver (GSE23343) [15]	10	4/6	—b	7	4/3	—	Affymetrix Human Genome U133 Plus 2.0 Array	Japan	2010	
aT2DM: Type 2 diabetes mellitus.

bInformation not available.

Inclusion and Exclusion Criteria

Data were restricted to human (Homo sapiens) samples, data set as the data type, expression profiling by array, tissue samples, and T2DM compared to controls (without diabetes). Thus, data from other organisms (Mus musculus, Rattus norvegicus, Xenopus laevis); series data; expression profiling by other methods (eg, massively parallel signature sequencing, reverse transcription-polymerase chain reaction, serial analysis of gene expression, genome variation or occupancy profiling by single-nucleotide polymorphism array, genome tiling array); nontissue samples (eg, blood, serum, semen, saliva, urine, body fluid); and data from patients with type 1 diabetes, gestational diabetes, or prediabetes were excluded.

The data collection process is summarized in Figure 1.

Figure 1 Flow diagram illustrating the process of data collection and number of data sets considered for inclusion.

Identification and Assortment of Differentially Expressed mRNAs

The DEGs were obtained from the five data sets using the online interactive tool GEO2R [16]. The cutoff for the selection was kept at the default of P<.05. The relaxed P-value cutoff was fixed for the initial selection because (1) we were subjecting the selected genes for a repeated analysis using ImaGEO software with a cutoff adjusted P<.05, and (2) the application of a stringent P-value cutoff in the initial selection did not enable obtaining an adequate number of genes from each data set for undertaking a meta-analysis. The overlapping DEGs among the three data sets of pancreatic tissues from patients with T2DM and controls (GSE38642, GSE25724, and GSE20966) were identified using the Venn diagram tool [17,18]. Subsequently, the common DEGs of these three data sets (GSE38642, GSE25724, and GSE20966) with those of the other two data sets for heart (GSE26887) and liver (GSE23343) samples were identified separately. The fold change expression distribution was visualized by heat maps and violin plots using the R limma (linear models for microarray data) package and Orange Data Mining software [19,20].

To check the quality of the data, quality control plots were assessed in the form of volcano plots, mean difference plots, and mean-variance trends. A volcano plot visualizes the DEGs by plotting the statistical significance against the magnitude of change, whereas the mean difference plot displays the log2 fold change against the average log2 expression level. The mean-variance trend, generated using the R packages plotSA and vooma, assesses the variance of the data. The workflow for the data processing and analysis is portrayed in Figure 2.

Figure 2 Flowchart of data processing and analysis. GEO: Gene Expression Omnibus; KEGG: Kyoto Encyclopedia of Genes and Genomes; T2DM: type 2 diabetes mellitus.

Functional Enrichment and KEGG Pathway Analysis

The DEGs were divided into three groups according to the tissue (Figure 2). The functional enrichment of each group related to T2DM was analyzed with the Database for Annotation, Visualization and Integrated Discovery (DAVID) tool for significant MF, CC, and BP GO terms. KEGG pathway analysis was performed with piNET, a versatile tool that integrates protein signatures with transcriptomic and proteomic signatures [21-23]. DAVID includes an analysis of KEGG pathways and enrichment significance of GO terms from the three categories (MF, CC, BP). We defined P<.05 as significantly enriched. The nonsignificant findings were manually removed.

PPI Network Construction and Identification of Hub Genes

The DEGs in the three groups were used to construct the PPI network using Search Tool for the Retrieval of Interacting Genes/Protein (STRING) [24]. We established the PPI network using only the overlapping DEGs with greater than 0.4 confidence score cutoffs. The “combined scores” were computed by integrating the probabilities from the various different types of evidence (by evidence channels), while correcting for the probability of randomly observing an interaction [25]. The number of interactions (by confidence level) were divided into four groups: (1) highest confidence (score≥0.90), (2) high confidence or better (score≥0.70), (3) medium confidence or better (score≥0.40), (4) low confidence links (score≥0.15). We chose medium confidence as the default setting given in STRING.

The interaction networks for each group were constructed by Cytoscape [26,27]. The Molecular Complex Detection (MCODE) [28] plugin of Cytoscape was employed to visualize significant genes in all three groups with a degree cutoff=2, node score cutoff=0.2, k-score=2, and maximum depth=100. The criteria for selecting the top 3 clusters were set as MCODE score≥3 and number of nodes≥3.

Integrative Gene Expression Meta-analysis

ImaGEO is a web tool for gene expression meta-analysis that was used to perform a comprehensive meta-analysis from all five data sets. For the retrieval and preprocessing of the data, the GEOquery package in R was used, followed by quality control, gene filtering expression, meta-analysis, and functional analysis. The meta-analysis was based on the functional modules with the MetaDE R package. For this study, we used the “effect size” parameter estimation with a fixed-effects model and an adjusted P value threshold of .05. The allowable missing values was kept at the default of 10%.

Target Prediction, Validation, and miRNA–Hub Gene Interaction

The top 10 targeting miRNAs of the hub genes were predicted by the well-established miRNA target prediction database miRNet 2.0 [29] with H. sapiens (human) as the selected organism. Default values for the degree of interaction and betweenness were retained. Common miRNAs and targeted mRNAs of all groups were sorted by the Venn diagram tool [30]. The network of all targeting miRNAs and the coexpressed mRNAs was created with FunRich and Cytoscape software. To validate the targeting miRNAs, we further sorted miRNA data sets in T2DM for comparison of differentially expressed miRNAs.

Functional Enrichment and KEGG Pathway Analysis for MiRNAs

All common miRNAs were enriched by MicroRNA Enrichment Turned Network (MIENTURNET) and KEGG pathway analysis [31]. MIENTURNET is a web tool based on the shiny package in R studio for both statistical and network-based analyses of miRNA-target enrichment. Functional enrichment was retrieved for the input list of genes, with the minimum interaction threshold set at 2 and an adjusted P value of .05. The input list infers possible experimental or computational evidence of miRNA-based regulation.

Results

Identification of DEGs in all Combined Groups

The five mRNA expression profiles of the GSE38642, GSE25724, GSE20966, GSE26887, and GSE23343 data sets, including 125 samples of the pancreas, heart, and liver tissues of patients with T2DM and controls without diabetes, were included in this study. We extracted 2852, 8631, 5501, 4210, and 3754 DEGs, respectively. The following sections describe the analysis of the DEGs derived from the datasets, and shown in Figures 3-14. In the pancreas data sets (GSE38642, GSE25724, GSE20966), a total of 321 common mRNAs were identified, 69 of which were upregulated and 95 were downregulated (Figure 3A-H and Supplementary Tables S1-S3 in Multimedia Appendix 1). The quality control plots for the DEGs are shown in Figure S1 and Figure S2 of Multimedia Appendix 1.

These Group 1 (pancreas) DEGs were then overlapped with the heart expression profile data set GSE26887, revealing a total of 70 common differentially expressed mRNAs, 5 of which were downregulated and 5 were upregulated. A total of 28 mRNAs with regulation in the opposite direction were identified (Group 2) (Figure 5A-K, Tables S4-S7 in Multimedia Appendix 1). Further, the Group 1 DEGs were overlapped with the liver data set GSE23343, and a total of 82 common differentially expressed mRNAs were obtained, out of which 8 were upregulated, 1 was downregulated, and 27 were regulated in opposite directions (Figure 7A-I, Tables S8-S11 in Multimedia Appendix 1).

DEGs for all three groups were used to establish the PPI networks (Figure 6E, Figure 8E, Figure 9A).

Figure 3 Differential mRNA expression of all three data sets (GSE38642, GSE25724, GSE20966) for Group 1 (pancreas tissues) in type 2 diabetes mellitus. (A-C) Heat maps of all, downregulated and upregulated differentially expressed genes (DEGs). Fold change expression (FCE) levels are displayed in ascending order from blue to yellow. (D-F) Venn diagrams of the total downregulated and upregulated DEGs of the three data sets. (G, H) Violin plots showing the entire FCE distribution of all three data sets for upregulated and downregulated common DEGs.

Figure 4 Differential mRNA expression of all three data sets (GSE38642, GSE25724, GSE20966) for Group 1 (pancreas tissues) in type 2 diabetes mellitus. (A-D) Enrichment analysis of common DEGs. (A) Kyoto Enclyclopedia of Genes and Genomes (KEGG) pathway enrichment: the connections are shown using red nodes (pathways) or brown nodes (DEGs) through the brown edges in a circle. The larger the size of the grey node, the more connected it is within the network. The density of red color indicates the number of connecting DEGs. (B) Gene Ontology cellular component terms. (C) Gene Ontology biological process terms. (D) Gene Ontology molecular function terms. Significant pathways represent adjusted P<.05 (false discovery rate).

Figure 5 (A-D) Heat maps of mRNA expression for the three data sets (GSE38642, GSE25724, GSE20966) of Group 1 (pancreas) and the GSE26887 (heart) data set showing all differentially expressed genes (DEGs), DEGs regulated in opposite directions, upregulated DEGs, and downregulated DEGs. (E-G) Venn diagrams of complete, upregulated, and downregulated common DEGs. The upper part of the heat map shows fold change expression (FCE) values represented by varying color densities. (H-K) Violin plots showing the entire FCE distribution of all four data sets of Group 1 (pancreas) and the heart data set.

Figure 6 (A) Enriched Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways. The connections are shown using red nodes (pathways) or brown nodes (DEGs) through the brown edges in a circle. The larger the size of the grey node, the more connected it is within the network. The density of red color indicates the number of connecting DEGs. (B-D) Gene Ontology enrichment for (B) cellular component, (C) biological process, and (D) molecular function terms. Significant pathways represent adjusted P<.05 (false discovery rate). (E) Top hub genes in the network. (F, G) Clusters using MCODE.

Figure 7 mRNA expression of three data sets (GSE38642, GSE25724, and GSE20966) of Group 1 (pancreas) and the GSE23343 data set (liver). (A-C) Venn diagrams of complete, upregulated, and downregulated common differentially expressed genes (DEGs). (D-F) Heat maps of common, oppositely regulated, and common upregulated DEGs. The upper part of the heat map shows the fold change in expression values reflected by respective color densities. (G-I) Violin plots showing the entire fold change expression (FCE) distribution of all four data sets for complete common DEGs, oppositely regulated DEGs, and common upregulated DEGs.

Figure 8 (A-D) Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway and Gene Ontology functional enrichment analysis of common DEGs. The connections are shown using red nodes (pathways) or brown nodes (DEGs) through the brown edges in a circle. The larger the size of the grey node, the more connected it is within the network. The density of red color indicates the number of connecting DEGs. (E) Protein-protein interaction networks of 82 overlapping DEGs of GSE23343 and co-expressed genes of Group 1 (321 DEGs) composed of 82 nodes and 56 edges. (F, G) Clusters from the network. Significant pathways represent adjusted P<.05 (false discovery rate).

Figure 9 (A) Common and the top hub genes (green) in the protein-protein interaction network. (B-D) Clusters of the network.

Figure 10 (A) Heat map of 16 common seed genes from the five data sets (pancreas, heart, and liver). The fold change expression levels are displayed in ascending order from blue to yellow. (B) Violin plot showing the entire fold change expression (FCE) distribution of all 16 common seed genes. (C) Venn diagram of common differentially expressed genes (DEGs). (D) Disease-gene interaction network. (E) Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways.

Figure 11 (A) Fold change expression levels of 16 common DEGs. (B) Top hub genes in the network (green) according to the criterion. (C, D, E, F) Clusters determined using MCODE.

Figure 12 Protein-microRNAs interactions (top 10 ranked) for (A) pancreas data sets, (B) heart data set, (C) liver data set, and (D) 16 common differentially expressed genes (DEGs) of all five data sets.

Figure 13 (A-E) Venn diagrams for (A) all microRNAs, (B) hsa-miR-1-3p, (C) hsa-miR-1-5p, (D) hsa-miR-155-5p, and (E) hsa-let-7b-5p with 16 common differentially expressed genes of the five data sets. (F) MicroRNA enrichment analysis of top four common targeted microRNAs.

Figure 14 Flow diagram for the process of microRNA data collection with the number of data sets considered for inclusion.

Functional Enrichment and KEGG Pathway Analysis

The enrichments for the three GO classes (BP, CC, and MF) of the 321 DEGs of Group 1 are shown in Figure 4B-D (also see Tables S12-S14 of Multimedia Appendix 1). KEGG pathway analysis showed that these genes were enriched in maturity-onset diabetes of the young, malaria, lysosome, insulin secretion, adrenergic signaling in cardiomyocytes, cell adhesion molecules, and T2DM pathways (Figure 4A, Table S15 of Multimedia Appendix 1).

The enrichments for the three GO classes of the 70 DEGs of Group 2 are shown in Figure 6B-D (also see Tables S16-S18 in Multimedia Appendix 1). The genes were mainly enriched in gap junction, melanoma, calcium signaling pathway, and GnRH signaling pathway (Figure 6A and Table S19 of Multimedia Appendix 1).

The enrichment terms for the three GO classes for the 82 DEGs in Group 3 are shown in Figure 8B-D (also see Tables S20-S22 in Multimedia Appendix 1). These genes were enriched in axon guidance (Figure 8A and Table S23 in Multimedia Appendix 1).

PPI Network and Hub Gene Identification

Group 1

The 321 overlapping DEGs of the GSE38642, GSE25724, GSE20966 pancreas data sets were used to establish the PPI network, which constituted 321 nodes, 737 edges, and a PPI enrichment P value <.001 at medium confidence (0.4) (Figure 9A). The top three significant clusters within the PPI were selected.

Cluster 1 (MCODE Score=9.556, 10 nodes, 43 edges) included the genes POLR1E, DDX10, URB1, HEATR1, DDX18, PDCD11, RSL1D1, RRP12, MAK16, and RRP15, which are mainly associated with insulin pathway, transforming growth factor (TGF)-β receptor signaling, and the mammalian target of rapamycin (mTOR) signaling pathway (Figure 9B).

Cluster 2 (MCODE score=8.000, 8 nodes, 28 edges) included the genes TRIM37, BTBD1, RNF19B, ASB4, KLHL22, SMURF1, FBXO7, and RNF6, which are mainly associated with insulin pathway, insulin-like growth factor 1 (IGF1) pathway, class I phosphatidylinositol-3-kinase (PI3K) signaling events mediated by protein kinase B (AKT), TGF-β receptor signaling, mTOR signaling pathway, platelet-derived growth factor receptor-beta signaling pathway, and epidermal growth factor (EGF) receptor (ERBB1) signaling pathway (Figure 9C).

Cluster 3 (MCODE score=4.000, 14 nodes, 26 edges) included the genes CD2, CD48, EDN3, GNAS, ITGAX, PDPN, GNRHR, RAC2, MADCAM1, WASL, GNAQ, PLEK, LAPTM5, and DVL2, which are associated with platelet activation, signaling, and aggregation; hemostasis, cell surface interactions at the vascular wall; integrin family cell surface; and IGF1 pathway (Figure 9D).

Group 2

The 70 overlapping DEGs of GSE26887 and coexpressed genes with Group 1 (321 genes) were used to establish the PPI network composed of 70 nodes, 32 edges, and a PPI enrichment P value of .05 at medium confidence (0.4). The top two significant clusters within the PPI were selected using the MCODE plugin of Cytoscape software (Figure 6E).

Cluster 1 (MCODE score=3.333, 4 nodes, 5 edges) included the genes RRP12, PDCD11, RRP15, and MRPL3 (Figure 6F). Cluster 2 (MCODE score=3.000, 3 nodes, 3 edges) included the genes PLEK, GNAQ, and IQSEC1, which are mainly associated with platelet activation, signaling, and aggregation; hemostasis; class I PI3K signaling events mediated by AKT; insulin pathway; mTOR signaling pathway; IGF1 pathway; and EGF receptor (ERBB1) signaling pathway (Figure 6G).

Group 3

The 82 overlapping DEGs of GSE23343 and the coexpressed genes of Group 1 (321 DEGs) were used to establish the PPI network composed of 82 nodes, 56 edges, and a PPI enrichment P value of .02 at medium confidence (0.4). The top two significant clusters are shown in Figure 8E.

Cluster 1 (MCODE score=03, 3 nodes, 3 edges) included the genes BTBD1, ASB4, and KLHL22, which were mainly associated with PI3K/AKT signaling in cancer (Figure 8F).

Cluster 2 (MCODE score=03, 3 nodes, 3 edges) included the genes DDX18, MAK16, and RRP15, which were mainly associated with insulin pathway, mTOR signaling pathway, IGF1 pathway, and EGF receptor (ERBB1) signaling pathway (Figure 8G).

Common Genes Among All Groups

A total of 16 overlapping DEGs were identified in all three groups. The hub genes of all data sets were ARHGEF9, CBFA2T2, CDK5RAP2, EEF1D, FGD6, FRMD4B, GNAQ, MAML3, NPC1, NRG1, OXR1, PCDH7, PFKFB2, PSIP1, SLC1A2, and TRAK1 (Table S24 in Multimedia Appendix 1). All 16 hub genes belonging to the five data sets were analyzed with the help of an expression heat map, violin plot, and Venn diagram, and their fold change expression levels were compared by bar plots and analyzed by the disease-gene interaction network and KEGG pathway (Figure 10A-E and Figure 10F; Table S25 in Multimedia Appendix 1).

The PPI network of the 16 hub genes and their related genes was established by protein STRING analysis. We selected 4 clusters from the PPI network using MCODE (Figure 11B). Cluster 1 (MCODE score=10, 10 nodes, 45 edges) included the genes MAML1, HEY2, NOTCH3, MAML3, NOTCH2, MAML2, NOTCH1, HEY1, RBPJ, and NOTCH4. The analysis also showed that cluster 1 contains MAML3 as a seed gene (Figure 11D). Cluster 2 (MCODE score=6.667, 7 nodes, 20 edges) included the genes EEF1A1, EEF1B2, EEF1G, RPLP0, RPS2, CARS, and EEF1D, with EEF1D as a seed gene (Figure 11E). Cluster 3 (MCODE score=5.714, 8 nodes, 20 edges) included the genes AKT1, NUMB, EGFR, ERBB3, GRB2, CDC42, PSEN1, and NRG1, with NRG1 as a seed gene (Figure 11C). Cluster 4 (MCODE score=4, 4 nodes, 6 edges) included the genes CDK5RAP2, CEP152, CENPJ, and PCNT, with CDK5RAP2 as the seed gene (Figure 11F).

Integrative Gene Expression and Meta-analysis

The number of genes with an adjusted P value <.05 for each data set revealed 4, 0, 3533, 171, and 1 significant genes from the meta-analysis, including ARHGEF9, SAMSN1, SLC1A2, RABGAP1L, OXR1, GNAQ, CBFA2T2, and RRP15. The 16 hub genes obtained from the gene expression meta-analysis are shown in Table S26 of Multimedia Appendix 1.

MicroRNA and Hub Gene Network

To investigate the regulatory relationship of the identified hub genes, their targeting miRNAs, and coexpressed network, the top 10 ranked DEG-targeting miRNAs were selected based on degree and betweenness values. The top 10 targeting miRNAs for the three groups were hsa-let-7b-5p, hsa-mir-107, hsa-mir-124-3p, hsa-mir-129-2-3p, hsa-mir-1-3p, hsa-mir-155-5p, hsa-mir-16-5p, hsa-mir-23b-3p, hsa-mir-27a-3p, and hsa-mir-34a-5p in Group 1 (pancreas); hsa-mir-16-5p, hsa-mir-124-3p, hsa-mir-1-3p, hsa-mir-27a-3p, hsa-let-7b-5p, hsa-mir-155-5p, hsa-mir-20a-5p, hsa-mir-26b-5p, hsa-mir-27b-3p, and hsa-mir-147a in Group 2 (heart); and hsa-mir-1-3p, hsa-mir-155-5p, hsa-mir-124-3p, hsa-let-7b-5p, hsa-mir-34a-5p, hsa-mir-101-3p, hsa-mir-15a-5p, hsa-mir-26a-5p, hsa-mir-181a-5p, and hsa-let-7a-5p in Group 3 (liver). The common hub genes were targeted by hsa-mir-16-5p, hsa-mir-27a-3p, hsa-let-7a-5p, hsa-let-7b-5p, hsa-mir-101-3p, hsa-mir-1-3p, hsa-mir-124-3p, hsa-mir-103a-3p, hsa-mir-122-5p, and hsa-mir-155-5p.

Four common miRNAs (hsa-let-7b-5p, hsa-mir-155-5p, hsa-mir-124-3p, hsa-mir-1-3p) were found in all three groups, targeting the 16 hub DEGs. The miRNAs and PPI networks representing multiple targeted nodes (DEGs) of particular miRNAs for all groups are shown in Figure 12A-D.

The common DEGs found in all three groups are targeted by hsa-miR-1-3p (GNAQ, PCDH7, CDK5RAP2, NPC1), hsa-let-7b-5p (OXR1), hsa-mir-155-5p (TRAK1, PSIP1), and hsa-mir-124-3p (NRG1). The common targeting important miRNAs (hsa-let-7b-5p, hsa-mir-155-5p, hsa-mir-124-3p, hsa-mir-1-3p) were mainly involved in the advanced glycation end products (AGE)-receptor for advanced glycation end products (RAGE) signaling pathway in diabetic complication and endocrine resistance (Figure 13A-F, Tables S27-S40 in Multimedia Appendix 1).

Target MiRNA Validation from Available Data Sets

To validate our miRNA prediction, we searched the database again and performed a thorough review of available miRNA data sets for T2DM. Our search yielded two miRNA data sets from renal tissue (GSE51674) and serum (GSE26168) samples. The flow diagram for the miRNA data set search is shown in Figure 14. However, we were not able to find any miRNA data set pertaining to the heart, pancreas, or liver tissue. Interestingly, on analysis of the data sets obtained from the renal tissue and serum, we observed a significant alteration for our predicted miRNAs in the renal tissue, which was conspicuously absent in the serum (Table 2). We assessed the expression of our predicted miRNAs in the renal tissue and serum by comparing the adjusted P values for both sample types. This analysis revealed that although the expression of miRNAs was significantly altered in renal tissues from patients with T2DM, the same was not observed in serum when compared with healthy controls. Our analysis highlights a paradoxical difference in the alteration of miRNAs in tissue and serum in T2DM.

Table 2 Validation of the fold change in expression levels of common microRNAs in the GSE51674 (kidney) and GSE26168 (serum) data sets.

MicroRNA	Adjusted P value	P value	t	B	FCa	logFC	
GSE51674 (kidney)	

	hsa-miR-124*b	<.001	<.001	–4.83	–0.85	0.50	–1.00	

	hsa-miR-1	<.001	<.001	6.06	0.92	19.42	4.28	

	hsa-miR-155	<.001	<.001	20.96	12.12	74.56	6.22	

	hsa-let-7b	<.001	<.001	4.60	–1.20	2.55	1.35	
GSE26168 (serum)	

	hsa-miR-124*	.46	.23	1.23	–6.22	1.01	0.02	

	hsa-miR-1	.86	.82	0.23	–6.95	1.00	0.00	

	hsa-miR-155	.46	.08	1.86	–5.33	1.04	0.05	

	hsa-let-7b	.46	.16	1.45	–5.94	237.61	7.89	
aFC: fold change.

b*indicates the star strand for miR-124.

Functional Enrichment of MiRNAs

The functional enrichment and pathway analysis by MIENTURNET revealed the top significant pathways for hsa-let-7b-5p, hsa-miR-124-3p, hsa-miR-1-3p, and hsa-miR-155-5p, including the PI3K-AKT signaling pathway (hsa-let-7b-5p, hsa-miR-124-3p, hsa-miR-1-3p, hsa-miR-155-5p), endocrine resistance (hsa-let-7b-5p, hsa-miR-124-3p, hsa-miR-1-3p, hsa-miR-155-5p), AGE-RAGE signaling pathway in diabetic complications (hsa-let-7b-5p, hsa-miR-124-3p, hsa-miR-1-3p, hsa-miR-155-5p), lipid and atherosclerosis (hsa-let-7b-5p, hsa-miR-1-3p, hsa-miR-155-5p), insulin signaling pathway (hsa-let-7b-5p, hsa-miR-124-3p, hsa-miR-1-3p), mitogen-activated protein kinase (MAPK) signaling pathway (hsa-let-7b-5p, hsa-miR-124-3p, hsa-miR-155-5p), fluid sheer stress and atherosclerosis (hsa-miR-124-3p, hsa-miR-1-3p), adipocytokine signaling pathway (hsa-miR-124-3p), diabetic cardiomyopathy (hsa-miR-124-3p, hsa-miR-1-3p), insulin resistance (hsa-miR-124-3p), carbohydrate digestion and absorption (hsa-miR-124-3p), regulation of lipolysis in adipocytes (hsa-miR-124-3p), glucagon signaling pathway (hsa-miR-124-3p), and TGF-β signaling pathway (hsa-miR-155-5p) (see Figure 13 and Table S41 of Multimedia Appendix 1).

Discussion

Principal Findings

Diabetes develops because of dysregulated β-cell and adipose-tissue responses to chronic fuel excess, which result in so-called nutrient spillover, insulin resistance, and metabolic stress. The latter causes multiple organ damage. However, insulin resistance, while forcing β-cells to work harder, may also have an important defensive role against nutrient-related toxic effects in tissues such as the heart [32]. The liver, which primarily regulates glucose homeostasis in the body, has a strong association with diabetes. Liver disease in diabetes can further be classified into liver disease related to diabetes, hepatogenous diabetes, and liver disease occurring coincidentally with diabetes mellitus [33]. Recently, knowledge on the pathogenesis and management of diabetes mellitus has been expanding; however, the disease is far from being effectively managed in a large proportion of patients. In silico analysis of disease pathways and exploration of various disease-related genes and their regulatory molecules have revealed unforeseen vistas. In this study, we analyzed tissue-specific microarray gene expression data sets from publicly available repositories employing a network-based bioinformatics pipeline. We identified DEGs common to different tissues of patients with T2DM and constructed diseasome networks to provide insights into the interactions of the genes. These DEGs enabled the identification of associated dysregulated molecular pathways in tissues and related GO terms. A large number of pathways and GO categories were reduced by manual curation after filtering using a P value threshold of .05.

Our analysis supports that diabetes is a multifactorial disease caused by multiple complex systems, with an abundant crossover between signaling pathways. For each data set included in the study, comprehensive analysis focusing on biological function and interaction of T2DM-related genes provided valuable information to understand the pathogenic effect of DEGs in various organs, including the heart, liver, and pancreas, of patients with diabetes. In this study, five mRNA expression profile data sets (GSE38642, GSE25724, GSE20966, GSE26887, and GSE23343), including 125 samples of the pancreas, heart, and liver tissues of patients with T2DM and controls without diabetes, were analyzed. A total of 16 seed genes were obtained after the final analysis. Some of these genes have been reported to play significant roles in T2DM and its related comorbidities. In a similar study that included DEG screening from a genome-wide association study (GWAS) catalog, Gupta and Vadde [34] identified four hub gene candidates, related signaling pathways, target miRNAs, and transcription factors. However, their selection criteria of the data sets chosen for analysis were different than those adopted in this study, which possibly accounts for the difference in results.

Neuregulin 1 (NRG1) and ERBB receptors are involved in glucose homeostasis. NRG1-ERBB pathway activation affects glucose metabolism in the liver. Mice with chronic NRG1 treatment showed increased p38 phosphorylation in the liver and improved glucose tolerance [35]. Myocardial NRG1/ERBB is altered during postmyocardial infarction heart failure associated with diabetes. NRG1 can improve the antioxidative function of the mitochondria, and thereby increase the proliferation and decrease the apoptosis of cardiomyocytes via ERBB/AKT signaling. This can explain the upregulated expression of NRG1 found in the cardiac tissue of patients with T2DM in our study. Moreover, the dysregulated insulin signaling pathway modifies titin-based cardiomyocyte tension, modulates diastolic function, impairs cyclic guanosine monophosphate (cGMP)–cGMP-dependent protein kinase signaling, and elevates protein kinase C-α activity, thereby causing titin-based cardiomyocyte stiffening in diabetic hearts. Chronic NRG1 application has shown promising results in the modulation of titin properties in T2DM-associated heart failure with a preserved ejection fraction [36]. Further, there are reports showing that hyperglycemia impairs NRG1/ERBB2 signaling by disrupting the balance between NRG1 isoforms, decreasing the expression of erbin, and correspondingly activating the MAPK pathway, ultimately aiding in the development of diabetic peripheral neuropathy [37]. Again, the downregulation of NRG1 expression in the liver found in this study points toward dysregulated glucose homeostasis.

PFKFB2 encodes 6-phosphofructo-2-kinase/fructose 2,6-bisphosphatase (PFK2/FBPase-2) isoform 2, a bifunctional enzyme involved in the synthesis and degradation of fructose 2,6-bisphosphate. Enhanced hepatic glycolysis in mice achieved by overexpressing PFK2/FBPase-2 in the liver resulted in reduced body weight and visceral fat content. PFK2/FBPase-2 is also a binding partner for glucokinase, which plays a pivotal role in the rate-limiting step of glucose-stimulated insulin secretion in pancreatic β-cells, and regulates obesity, insulin secretory dysfunction, and T2DM [38,39]. The loss of PFK2 content as a result of reduced insulin signaling impairs its regulatory function of glycolysis and elevates the levels of early glycolytic intermediates. Although this may be beneficial in the fasting state to conserve systemic glucose, it represents a pathological impairment in diabetes mellitus [40]. Interestingly, PFKFB2, among a few other genes, showed opposing expression changes in the pancreas (downregulation) and heart (upregulation). This is likely due to the impaired insulin secretion pathway in pancreatic β-cells, in which PFKFB2 plays an important role [39]. Moreover, PFKFB2 is known to alleviate myocardial injury; hence, the increased expression level in the heart is possibly a protective mechanism [41].

CDK5 regulatory subunit associated protein (CDK5RAP) 1, 2, and 3 were all found to be differentially upregulated in four data sets, except GSE23343 in which these genes were downregulated. These genes have been associated with neuronal development and spindle checkpoint function [42]. FRMD4B plays a vital role in cardiac activity regulation. However, the effect varies in different populations due to polymorphisms. FRMD4B has shown to be associated with ischemic heart failure in a European population but not in other populations [43]. The G-protein Gq, encoded by GNAQ, is a crucial key regulator of the insulin secretion pathway that is involved in glucose metabolism, and a functional GNAQ promoter haplotype was associated with altered Gq expression and with insulin resistance and obesity in women with polycystic ovary syndrome [44]. The Niemann-Pick type C1 (NPC1) protein regulates the transport of cholesterol and fatty acids from late endosomes/lysosomes and has a central role in maintaining lipid homeostasis. In humans, GWAS and post-GWAS highlighted the implication of common variants in NPC1 in adult-onset obesity, body fat mass, and T2DM. Heterozygous human carriers of rare loss-of-function coding variants in NPC1 display an increased risk of morbid adult obesity [45]. Another significant DEG pair was orexin A and B, which regulate a variety of physiological functions. The biological effects of these neuropeptides occur through OXR1, a G-protein coupled receptor. There is growing evidence that orexins regulate body weight, glucose homeostasis, and insulin sensitivity, and promote energy expenditure, thus protecting against obesity by interacting with brown adipocytes. Further, orexins control brown and white adipocytes as well as pancreatic α- and β-cell functions [46,47]. Single-cell RNA sequencing from samples of patients with gestational diabetes mellitus revealed SLC1A2 as a novel marker for syncytiotrophoblasts [48]. Such cell-type-specific marker genes in particular disease states can open new avenues of tissue-targeted therapeutic intervention. Among the other DEGs, EEF1D regulates lipid synthesis via the PI3K/AKT, PPAR, and AMPK pathways [49]. CBFA2T2 is a key regulator of adipogenic differentiation through CEBPA [50]. Further, these seed genes were analyzed as possible miRNA targets in silico, which revealed the top 10 miRNAs for each of the pancreas, liver, and heart tissues, as well as for the 16 seed genes. The role of miRNAs in the regulation of the underlying pathogenic mechanisms of diabetes and diabetic complications is well established [7,51] Several of the target miRNAs for the seed genes have already been explored in T2DM, and our in silico analysis further confirms their candidature as potential biomarkers as well as therapeutic targets. In fact, miR-124-3p was interconnected to 7 of the 16 seed genes. Pan et al [52] studied mouse primary hepatocytes and observed that regulation of miR-124-3p plays an important role in turning the hepatocytes into insulin-producing cells. A recent analysis of weighted genes in diabetic retinopathy concluded miR-124-3p to be a pivotal regulatory molecule in the underlying pathogenesis [6]. Furthermore, in isolated myocardial cells, NRG1 expression was observed to be downregulated while miR-124-3p expression was upregulated in ischemia/reperfusion injury [53], which also supports our finding of this miRNA-mRNA target interaction. The miRNA hsa-miR-124-3p affects the immune status of patients with T2DM through its interaction with the obesity-related immune cytokines [54].

Three other miRNAs, namely miR-155-5p, miR-1-3p, and let-7b-5p, were also commonly identified in all three groups. Likewise, the role of miR-155-5p in diabetes has been widely studied, especially as a marker in diabetic kidney disease (DKD) [55-57]. The expression of miR-155-5p is positively associated with urinary microalbumin and has good diagnostic and prognostic value in patients with DKD [56]. Further, dihydromyricetin attenuates renal interstitial fibrosis by regulating PTEN signaling, a critical element in the pathogenesis of DKD, through miR-155-5p [58,59]. Recently, Zhou and colleagues [60] showed that metformin can relieve inflammation and fibrosis in patients with DKD by acting through an inflammation axis involving miR-155-5p. Some recent studies have also shown that miR-155-5p interferes with immune dysregulation in COVID-19 patients with diabetes or other comorbidities [61,62]. Further, all four miRNAs were found to be involved in regulating the endocrine resistance and AGE-RAGE pathways, which is in line with recent findings [63].

The differing trend in miRNA expression observed in our comparison of miRNA data sets from serum and renal tissue in T2DM highlights the necessity to further explore the tissue-specific alterations in T2DM to better comprehend its role in various tissues.

Limitations

The main limitation of this study is that it was based on an in silico analysis; therefore, further validation of the identified novel hub genes and miRNAs is still required based on laboratory experiments with human T2DM samples. The data sets were compiled using different arrays on the Affymetrix platform, and the patient populations belong to multiple ethnic groups, which may account for some of the variability in the results. Furthermore, the predicted miRNAs in this study could not be validated within the same tissue data sets. However, the functional enrichment for the miRNAs highlighted some significant pathways related to T2DM, its complications, and its pathogenic mechanisms.

Conclusion

The aim of this study was to identify the tissue-specific differential expression of genes, especially pertaining to the heart, liver, and pancreas, in T2DM. From Group 1 (pancreas: 374 DEGs), Group 2 (heart: 86 DEGs), and Group 3 (liver: 97 DEGs), we identified a total of 16 common DEGS (ARHGEF9, CBFA2T2, CDK5RAP2, EEF1D, FGD6, FRMD4B, GNAQ, MAML3, NPC1, NRG1, OXR1, PCDH7, PFKFB2, PSIP1, SLC1A2, and TRAK1) in the selected data sets. Further, we identified the top four common miRNAs (hsa-let-7b-5p, hsa-miR-124-3p, hsa-miR-1-3p, has-miR-155-5p) targeting the 16 common DEGs. Although we were not able to find any miRNA data set pertaining to the heart, pancreas, or liver tissue, we observed significant alterations of our predicted miRNAs in renal tissue. Interestingly, this significant alteration was conspicuously absent in the serum. The miRNAs identified in this study are involved in regulating various pathways, including the PI3K-AKT signaling pathway, endocrine resistance, and the AGE-RAGE signaling pathway. Moreover, the differing trend in miRNA expression observed in our comparison of miRNA data sets from the serum and renal tissue in T2DM highlights the necessity to further explore the tissue-specific alteration in T2DM to better comprehend its role in various tissues.

The authors are grateful to All India Institute of Medical Sciences Jodhpur for providing the research facility to perform this in silico experiment. MK is supported by a senior research fellowship of The University Grants Commission of India (NOV2017-361200).

Multimedia Appendix 1 Supplementary data: Figures S1-S2, Tables S1-S41.

Abbreviations

AGE advanced glycation end products

AKT protein kinase B

BP biological process

CC cellular component

CDK5RAP CDK5 regulatory subunit associated protein

cGMP cyclic guanosine monophosphate

DAVID Database for Annotation, Visualization and Integrated Discovery

DEG differentially expressed gene

DKD diabetic kidney disease

EGF epidermal growth factor

ERBB1 epidermal growth factor receptor

GEO Gene Expression Omnibus

GO Gene Ontology

GWAS genome-wide association study

IGF1 insulin-like growth factor 1

KEGG Kyoto Encyclopedia of Genes and Genomes

limma linear models for microarray data

MAPK mitogen-activated protein kinase

MCODE Molecular Complex Detection

MF molecular function

MIENTURNET MicroRNA Enrichment Turned Network

miRNA microRNA

mTOR mammalian target of rapamycin

NPC1 Niemann-Pick type C1

NRG1 neuregulin 1

PFK2 6-phosphofructo-2-kinase/fructose 2,6-bisphosphatase isoform 2

PI3K phosphoinositide 3-kinase

PPI protein-protein interaction

RAGE receptor of advanced glycation end products

STRING Search Tool for the Retrieval of Interacting Genes/Proteins

T2DM type 2 diabetes mellitus

TGF transforming growth factor

Authors' Contributions: Concept and design: MK, DR, PP; data acquisition, analysis, and interpretation: MK, ST, AG, PP; manuscript drafting: MK, DR, ST, AG, PS, PP; manuscript revision: MK, DR, ST, PS, PP; project supervision: PP.

Conflicts of Interest: None declared.
==== Refs
1 Behera PM Behera DK Panda A Dixit A Padhi P In silico expressed sequence tag analysis in identification of probable diabetic genes as virtual therapeutic targets Biomed Res Int 2013 2013 704818 10.1155/2013/704818 23509765 23509765
2 Hruby A Hu FB The epidemiology of obesity: a big picture Pharmacoeconomics 2015 07 4 33 7 673 689 10.1007/s40273-014-0243-x 25471927 25471927
3 DF Diabetes Atlas Ninth edition 2019 International Diabetes Federation 2021-07-09 https://diabetesatlas.org/upload/resources/material/20200302_133351_IDFATLAS9e-final-web.pdf
4 Vaidya V Gangan N Sheehan J Impact of cardiovascular complications among patients with Type 2 diabetes mellitus: a systematic review Expert Rev Pharmacoecon Outcomes Res 2015 06 31 15 3 487 497 10.1586/14737167.2015.1024661 25824591 25824591
5 Pordzik J Jakubik D Jarosz-Popek J Wicik Z Eyileten C De Rosa S Indolfi C Siller-Matula JM Czajka P Postula M Significance of circulating microRNAs in diabetes mellitus type 2 and platelet reactivity: bioinformatic analysis and review Cardiovasc Diabetol 2019 08 30 18 1 113 10.1186/s12933-019-0918-x 31470851 10.1186/s12933-019-0918-x 31470851
6 You Z Zhang Y Li B Zhu X Shi K Bioinformatics analysis of weighted genes in diabetic retinopathy Invest Ophthalmol Vis Sci 2018 11 01 59 13 5558 5563 10.1167/iovs.18-25515 30480744 2716394 30480744
7 Roy D Modi A Khokhar M Sankanagoudar S Yadav D Sharma S Purohit P Sharma P MicroRNA 21 emerging role in diabetic complications: a critical update Curr Diabetes Rev 2021 02 17 2 122 135 10.2174/1573399816666200503035035 32359340 CDR-EPUB-106300 32359340
8 Plenge RM Disciplined approach to drug discovery and early development Sci Transl Med 2016 07 27 8 349 349ps15 10.1126/scitranslmed.aaf2608 27464747 8/349/349ps15
9 Taneera J Lang S Sharma A Fadista J Zhou Y Ahlqvist E Jonsson A Lyssenko V Vikman P Hansson O Parikh H Korsgren O Soni A Krus U Zhang E Jing X Esguerra JLS Wollheim CB Salehi A Rosengren A Renström Erik Groop L A systems genetics approach identifies genes and pathways for type 2 diabetes in human islets Cell Metab 2012 07 03 16 1 122 134 10.1016/j.cmet.2012.06.006 22768844 S1550-4131(12)00243-4 22768844
10 Taneera J Fadista J Ahlqvist E Zhang M Wierup N Renström E Groop L Expression profiling of cell cycle genes in human pancreatic islets with and without type 2 diabetes Mol Cell Endocrinol 2013 08 15 375 1-2 35 42 10.1016/j.mce.2013.05.003 23707792 S0303-7207(13)00206-2 23707792
11 Kanatsuna N Taneera J Vaziri-Sani F Wierup N Larsson HE Delli A Skärstrand H Balhuizen A Bennet H Steiner DF Törn C Fex M Lernmark Autoimmunity against INS-IGF2 protein expressed in human pancreatic islets J Biol Chem 2013 10 04 288 40 29013 29023 10.1074/jbc.M113.478222 23935095 S0021-9258(20)48874-X 23935095
12 Dominguez V Raimondi C Somanath S Bugliani M Loder MK Edling CE Divecha N da Silva-Xavier G Marselli L Persaud SJ Turner MD Rutter GA Marchetti P Falasca M Maffucci T Class II phosphoinositide 3-kinase regulates exocytosis of insulin granules in pancreatic beta cells J Biol Chem 2011 02 11 286 6 4216 4225 10.1074/jbc.M110.200295 21127054 S0021-9258(20)56124-3 21127054
13 Marselli L Thorne J Dahiya S Sgroi DC Sharma A Bonner-Weir S Marchetti P Weir GC Gene expression profiles of beta-cell enriched tissue obtained by laser capture microdissection from subjects with type 2 diabetes PLoS One 2010 07 13 5 7 e11499 10.1371/journal.pone.0011499 20644627 20644627
14 Greco S Fasanaro P Castelvecchio S D'Alessandra Y Arcelli D Di Donato M Malavazos A Capogrossi MC Menicanti L Martelli F MicroRNA dysregulation in diabetic ischemic heart failure patients Diabetes 2012 06 61 6 1633 1641 10.2337/db11-0952 22427379 db11-0952 22427379
15 Misu H Takamura T Takayama H Hayashi H Matsuzawa-Nagata N Kurita S Ishikura K Ando H Takeshita Y Ota T Sakurai M Yamashita T Mizukoshi E Yamashita T Honda M Miyamoto K Kubota T Kubota N Kadowaki T Kim H Lee I Minokoshi Y Saito Y Takahashi K Yamada Y Takakura N Kaneko S A liver-derived secretory protein, selenoprotein P, causes insulin resistance Cell Metab 2010 11 03 12 5 483 495 10.1016/j.cmet.2010.09.015 21035759 S1550-4131(10)00351-7 21035759
16 GEO2R National Center for Biotechnology Information (NCBI) 2022-02-04 https://www.ncbi.nlm.nih.gov/geo/geo2r/
17 Pathan M Keerthikumar S Chisanga D Alessandro R Ang C Askenase P Batagov AO Benito-Martin A Camussi G Clayton A Collino F Di Vizio D Falcon-Perez JM Fonseca P Fonseka P Fontana S Gho YS Hendrix A Hoen EN Iraci N Kastaniegaard K Kislinger T Kowal J Kurochkin IV Leonardi T Liang Y Llorente A Lunavat TR Maji S Monteleone F Øverbye A Panaretakis T Patel T Peinado H Pluchino S Principe S Ronquist G Royo F Sahoo S Spinelli C Stensballe A Théry C van Herwijnen MJC Wauben M Welton JL Zhao K Mathivanan S A novel community driven software for functional enrichment analysis of extracellular vesicles data J Extracell Vesicles 2017 12 6 1 1321455 10.1080/20013078.2017.1321455 28717418 1321455 28717418
18 Pathan M Keerthikumar S Ang C Gangoda L Quek CY Williamson NA Mouradov D Sieber OM Simpson RJ Salim A Bacic A Hill AF Stroud DA Ryan MT Agbinya JI Mariadason JM Burgess AW Mathivanan S FunRich: An open access standalone functional enrichment and interaction network analysis tool Proteomics 2015 08 17 15 15 2597 2601 10.1002/pmic.201400515 25921073 25921073
19 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 04 20 43 7 e47 10.1093/nar/gkv007 25605792 gkv007 25605792
20 Demsar J Curk T Gorup C Hocevar T Milutinovic M Mozina M Polajnar M Toplak M Staric A Stajdohar M Umek L Zagar L Zbontar J Zitnik M Zupan B Orange: data mining toolbox in Python J Machine Learn Res 2013 08 14 2349 2353
21 Huang DW Sherman BT Lempicki RA Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources Nat Protoc 2009 12 18 4 1 44 57 10.1038/nprot.2008.211 19131956 nprot.2008.211 19131956
22 Huang D Sherman B Lempicki R Bioinformatics enrichment tools: paths toward the comprehensive functional analysis of large gene lists Nucleic Acids Res 2009 01 37 1 1 13 10.1093/nar/gkn923 19033363 gkn923 19033363
23 Shamsaei B Chojnacki S Pilarczyk M Najafabadi M Niu W Chen C Ross K Matlock A Muhlich J Chutipongtanate S Zheng J Turner J Vidović D Jaffe J MacCoss M Wu C Pillai A Ma'ayan A Schürer S Kouril M Medvedovic M Meller J piNET: a versatile web platform for downstream analysis and visualization of proteomics data Nucleic Acids Res 2020 07 02 48 W1 W85 W93 10.1093/nar/gkaa436 32469073 5848491 32469073
24 Szklarczyk D Gable AL Lyon D Junge A Wyder S Huerta-Cepas J Simonovic M Doncheva NT Morris JH Bork P Jensen LJ von Mering C STRING v11: protein-protein association networks with increased coverage, supporting functional discovery in genome-wide experimental datasets Nucleic Acids Res 2019 01 08 47 D1 D607 D613 10.1093/nar/gky1131 30476243 5198476 30476243
25 von Mering C Jensen LJ Snel B Hooper SD Krupp M Foglierini M Jouffre N Huynen MA Bork P STRING: known and predicted protein-protein associations, integrated and transferred across organisms Nucleic Acids Res 2005 01 01 33 Database issue D433 D437 10.1093/nar/gki005 15608232 33/suppl_1/D433 15608232
26 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 11 13 11 2498 2504 10.1101/gr.1239303 14597658 13/11/2498 14597658
27 Otasek D Morris JH Bouças J Pico AR Demchak B Cytoscape Automation: empowering workflow-based network analysis Genome Biol 2019 09 02 20 1 185 10.1186/s13059-019-1758-4 31477170 10.1186/s13059-019-1758-4 31477170
28 Bader GD Hogue CW An automated method for finding molecular complexes in large protein interaction networks BMC Bioinformatics 2003 01 13 4 1 2 10.1186/1471-2105-4-2 12525261 12525261
29 Chang L Zhou G Soufan O Xia J miRNet 2.0: network-based visual analytics for miRNA functional analysis and systems biology Nucleic Acids Res 2020 07 02 48 W1 W244 W251 10.1093/nar/gkaa467 32484539 5850315 32484539
30 Calculate and draw custom Venn diagrams Bioinformatics & Evolutionary Genomics 2022-02-04 http://bioinformatics.psb.ugent.be/webtools/Venn/
31 Licursi V Conte F Fiscon G Paci P MIENTURNET: an interactive web tool for microRNA-target enrichment and network-based analysis BMC Bioinformatics 2019 11 04 20 1 545 10.1186/s12859-019-3105-x 31684860 10.1186/s12859-019-3105-x 31684860
32 Nolan CJ Damm P Prentki M Type 2 diabetes across generations: from pathophysiology to prevention and management Lancet 2011 07 09 378 9786 169 181 10.1016/S0140-6736(11)60614-4 21705072 S0140-6736(11)60614-4 21705072
33 Hamed AE Elsahar M Elwan NM El-Nakeep S Naguib M Soliman HH Ahmed Aboubakr A AbdelMaqsod A Sedrak H Assaad SN Elwakil R Esmat G Salh S Mostafa T Mogawer S Sadek SE Saber MM Ezelarab H Mahmoud AA Sultan S El Kassas M Kamal E ElSayed NM Moussa S Managing diabetes and liver disease association Arab J Gastroenterol 2018 12 19 4 166 179 10.1016/j.ajg.2018.08.003 30420265 S1687-1979(18)30065-0 30420265
34 Gupta MK Vadde R Identification and characterization of differentially expressed genes in type 2 diabetes using in silico approach Comput Biol Chem 2019 04 79 24 35 10.1016/j.compbiolchem.2019.01.010 30708140 S1476-9271(18)30223-8 30708140
35 Ennequin G Caillaud K Chavanelle V Teixeira A Etienne M Li X Boisseau N Sirvent P Neuregulin 1 treatment improves glucose tolerance in diabetic db/db mice, but not in healthy mice Arch Physiol Biochem 2020 10 17 126 4 320 325 10.1080/13813455.2018.1534243 30449185 30449185
36 Hopf A Andresen C Kötter S Isić M Ulrich K Sahin S Bongardt S Röll W Drove F Scheerer N Vandekerckhove L De Keulenaer GW Hamdani N Linke WA Krüger M Diabetes-induced cardiomyocyte passive stiffening is caused by impaired insulin-dependent titin modification and can be modulated by neuregulin-1 Circ Res 2018 07 20 123 3 342 355 10.1161/circresaha.117.312166 29760016
37 Pan P Dobrowsky RT Differential expression of neuregulin-1 isoforms and downregulation of erbin are associated with Erb B2 receptor activation in diabetic peripheral neuropathy Acta Neuropathol Commun 2013 07 17 1 1 39 10.1186/2051-5960-1-39 24252174 2051-5960-1-39 24252174
38 Muller Y Piaggi P Hanson R Kobes S Bhutta S Abdussamad M Leak-Johnson T Kretzler M Huang K Weil E Nelson R Knowler W Bogardus C Baier L A cis-eQTL in PFKFB2 is associated with diabetic nephropathy, adiposity and insulin secretion in American Indians Hum Mol Genet 2015 05 15 24 10 2985 2996 10.1093/hmg/ddv040 25662186 ddv040 25662186
39 Arden C Hampson L Huang G Shaw J Aldibbiat A Holliman G Manas D Khan S Lange A Agius L A role for PFK-2/FBPase-2, as distinct from fructose 2,6-bisphosphate, in regulation of insulin secretion in pancreatic beta-cells Biochem J 2008 04 01 411 1 41 51 10.1042/BJ20070962 18039179 BJ20070962 18039179
40 Bockus LB Matsuzaki S Vadvalkar SS Young ZT Giorgione JR Newhardt MF Kinter M Humphries KM Cardiac insulin signaling regulates glycolysis through phosphofructokinase 2 content and activity J Am Heart Assoc 2017 12 04 6 12 3846 10.1161/JAHA.117.007159 29203581 JAHA.117.007159
41 Gao J Feng W Lv W Liu W Fu C HIF-1/AKT signaling-activated PFKFB2 alleviates cardiac dysfunction and cardiomyocyte apoptosis in response to hypoxia Int Heart J 2021 03 30 62 2 350 358 10.1536/ihj.20-315 33678793 33678793
42 Marselli L Thorne J Dahiya S Sgroi DC Sharma A Bonner-Weir S Marchetti P Weir GC Gene expression profiles of Beta-cell enriched tissue obtained by laser capture microdissection from subjects with type 2 diabetes PLoS One 2010 07 13 5 7 e11499 10.1371/journal.pone.0011499 20644627 20644627
43 Lee HJ Yun JH Kim H Jang HB Park SI Lee H 244-LB: Glutamate is associated with type 2 diabetes through PLG regulation Diabetes 2019 06 04 68 Supplement 1 244-LB 10.2337/db19-244-lb
44 Saini C Petrenko V Pulimeno P Giovannoni L Berney T Hebrok M Howald C Dermitzakis ET Dibner C A functional circadian clock is required for proper insulin secretion by human pancreatic islet cells Diabetes Obes Metab 2016 04 22 18 4 355 365 10.1111/dom.12616 26662378 26662378
45 Lamri A Pigeyre M Garver W Meyre D The extending spectrum of NPC1-related human disorders: from Niemann-Pick C1 disease to obesity Endocr Rev 2018 04 01 39 2 192 220 10.1210/er.2017-00176 29325023 4794852 29325023
46 Skrzypski M Billert M Nowak KW Strowski MZ The role of orexin in controlling the activity of the adipo-pancreatic axis J Endocrinol 2018 08 238 2 R95 R108 10.1530/JOE-18-0122 29848609 JOE-18-0122 29848609
47 Sellayah D Sikder D Orexin receptor-1 mediates brown fat developmental differentiation Adipocyte 2012 01 01 1 1 58 63 10.4161/adip.18965 23700511 2011ADIPOCYTE020R 23700511
48 Yang Y Guo F Peng Y Chen R Zhou W Wang H OuYang J Yu B Xu Z Transcriptomic profiling of human placenta in gestational diabetes mellitus at the single-cell level Front Endocrinol 2021 5 7 12 679582 10.3389/fendo.2021.679582 34025588
49 Hou Y Xie Y Yang S Han B Shi L Bai X Liang R Dong T Zhang S Zhang Q Sun D EEF1D facilitates milk lipid synthesis by regulation of PI3K-Akt signaling in mammals FASEB J 2021 05 35 5 e21455 10.1096/fj.202000682RR 33913197 33913197
50 Luo J Dou L Yang Z Zhou Z Huang H CBFA2T2 promotes adipogenic differentiation of mesenchymal stem cells by regulating CEBPA Biochem Biophys Res Commun 2020 08 20 529 2 133 139 10.1016/j.bbrc.2020.05.120 32703401 S0006-291X(20)31042-1 32703401
51 Khokhar M Roy D Bajpai NK Bohra GK Yadav D Sharma P Purohit P Metformin mediates MicroRNA-21 regulated circulating matrix metalloproteinase-9 in diabetic nephropathy: an in-silico and clinical study Arch Physiol Biochem 2021 06 04 1 11 10.1080/13813455.2021.1922457 34087084
52 Pan G Liu Q Xin H Liu J The key regulation of miR-124-3p during reprogramming of primary mouse hepatocytes into insulin-producing cells Biochem Biophys Res Commun 2020 02 05 522 2 315 321 10.1016/j.bbrc.2019.11.058 31761319 S0006-291X(19)32178-3 31761319
53 Liu Y Ke X Guo W Wang X Peng C Liao Z Liu Q Zhou Y Circ-RHOJ.1 regulated myocardial cell proliferationand apoptosis via targeting the miR-124-3p/NRG-1 axisin myocardial ischemia/reperfusion injury Arch Med Sci 2019 8 7 1 14 10.5114/aoms.2019.87205
54 Duan J Liu H Chen J Li X Li P Zhang R Changes in gene expression of adipose tissue CD14 cells in patients with type 2 diabetes mellitus and their relationship with environmental factors Zhong Nan Da Xue Xue Bao Yi Xue Ban 2021 01 28 46 1 1 10 10.11817/j.issn.1672-7347.2021.190558 33678630 1672-7347(2021)01-0001-10 33678630
55 Wang Y Zheng Z Jia Y Yang Y Xue Y Role of p53/miR-155-5p/sirt1 loop in renal tubular injury of diabetic kidney disease J Transl Med 2018 05 30 16 1 146 10.1186/s12967-018-1486-7 29848325 10.1186/s12967-018-1486-7 29848325
56 Bai X Luo Q Tan K Guo L Diagnostic value of VDBP and miR-155-5p in diabetic nephropathy and the correlation with urinary microalbumin Exp Ther Med 2020 11 11 20 5 86 10.3892/etm.2020.9214 32968443 ETM-0-0-09214 32968443
57 Wang G Wu B Zhang B Wang K Wang H LncRNA CTBP1-AS2 alleviates high glucose-induced oxidative stress, ECM accumulation, and inflammation in diabetic nephropathy via miR-155-5p/FOXO1 axis Biochem Biophys Res Commun 2020 11 05 532 2 308 314 10.1016/j.bbrc.2020.08.073 32868076 S0006-291X(20)31670-3 32868076
58 Guo L Tan K Luo Q Bai X Dihydromyricetin promotes autophagy and attenuates renal interstitial fibrosis by regulating miR-155-5p/PTEN signaling in diabetic nephropathy Bosn J Basic Med Sci 2020 08 03 20 3 372 380 10.17305/bjbms.2019.4410 31668144 31668144
59 Khokhar M Roy D Modi A Agarwal R Yadav D Purohit P Sharma P Perspectives on the role of PTEN in diabetic nephropathy: an update Crit Rev Clin Lab Sci 2020 11 20 57 7 470 483 10.1080/10408363.2020.1746735 32306805 32306805
60 Zhou Y Ma X Han J Yang M Lv C Shao Y Wang Y Kang J Wang Q Metformin regulates inflammation and fibrosis in diabetic kidney disease through TNC/TLR4/NF-κB/miR-155-5p inflammatory loop World J Diabetes 2021 01 15 12 1 19 46 10.4239/wjd.v12.i1.19 33520106 33520106
61 Khokhar M Purohit P Roy D Tomo S Gadwal A Modi A Banerjee M Sharma P Acute kidney injury in COVID 19 - an update on pathophysiology and management modalities Arch Physiol Biochem 2020 12 15 1 14 10.1080/13813455.2020.1856141 33320717
62 Khokhar M Tomo S Purohit P MicroRNAs based regulation of cytokine regulating immune expressed genes and their transcription factors in COVID-19 Meta Gene 2022 02 31 100990 10.1016/j.mgene.2021.100990 34722158 S2214-5400(21)00141-9 34722158
63 Agarwal RG Khokhar M Purohit P Modi A Bajpai NK Bohra GK Garg MK Sharma P A clinical and in-silico study of MicroRNA-21 and growth differentiation factor-15 expression in pre-diabetes, type 2 diabetes and diabetic nephropathy Minerva Endocrinol 2022 02 01 online ahead of print 10.23736/S2724-6507.22.03646-6 35103454 S2724-6507.22.03646-6
