Next Article in Journal
Revealing Natural Intracellular Peptides in Gills of Seahorse Hippocampus reidi
Previous Article in Journal
Tutorial for Stopped-Flow Water Flux Measurements: Why a Report about “Ultrafast Water Permeation through Nanochannels with a Densely Fluorous Interior Surface” Is Flawed
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Machine Learning Model Based on Insulin Resistance Metagenes Underpins Genetic Basis of Type 2 Diabetes

1
Department of Computer Engineering & Applications, Institute of Engineering & Technology, GLA University, Mathura 281406, India
2
Department of Medicine, Sawai Man Singh Medical College and Hospital, Jaipur 302004, India
3
Department of Endocrinology, Sawai Man Singh Medical College and Hospital, Jaipur 302004, India
*
Author to whom correspondence should be addressed.
Biomolecules 2023, 13(3), 432; https://doi.org/10.3390/biom13030432
Submission received: 12 January 2023 / Revised: 12 February 2023 / Accepted: 14 February 2023 / Published: 24 February 2023
(This article belongs to the Section Biomacromolecules: Carbohydrates)

Abstract

:
Insulin resistance (IR) is considered the precursor and the key pathophysiological mechanism of type 2 diabetes (T2D) and metabolic syndrome (MetS). However, the pathways that IR shares with T2D are not clearly understood. Meta-analysis of multiple DNA microarray datasets could provide a robust set of metagenes identified across multiple studies. These metagenes would likely include a subset of genes (key metagenes) shared by both IR and T2D, and possibly responsible for the transition between them. In this study, we attempted to find these key metagenes using a feature selection method, LASSO, and then used the expression profiles of these genes to train five machine learning models: LASSO, SVM, XGBoost, Random Forest, and ANN. Among them, ANN performed well, with an area under the curve (AUC) > 95%. It also demonstrated fairly good performance in differentiating diabetics from normal glucose tolerant (NGT) persons in the test dataset, with 73% accuracy across 64 human adipose tissue samples. Furthermore, these core metagenes were also enriched in diabetes-associated terms and were found in previous genome-wide association studies of T2D and its associated glycemic traits HOMA-IR and HOMA-B. Therefore, this metagenome deserves further investigation with regard to the cardinal molecular pathological defects/pathways underlying both IR and T2D.

Graphical Abstract

1. Introduction

The global prevalence of insulin resistance (IR) has been estimated to range from 15.5 to 46.5% among adults [1]. It has been identified as a central pathophysiological factor of several endocrine–metabolic disorders, such as type 2 diabetes (T2D), high blood pressure, dyslipidemias, polycystic ovary syndrome, metabolic syndrome (MetS), cardiovascular diseases, etc. More significantly, it precedes and could lead to T2D, which now affects over 422 million people globally and accounts for 1.5 million deaths annually [2].
T2D is the most frequent subtype of diabetes and it is characterized by alterations of blood glucose levels due to varying combinations of IR and a relative deficiency of insulin secretion by the pancreatic β cells. If untreated, T2D may escalate to various microvascular and macrovascular complications such as retinopathy, diabetic kidney disease, peripheral neuropathy, atherosclerotic vascular disease, etc. Its underlying pathogenesis is complex and includes an approximately 40–70% contribution of genetic factors [3,4,5]. The majority of the genetic risk loci, however, have been found to decrease insulin secretion rather than its action [6].
IR is considered to be the major underlying pathophysiological defect in both obesity and early-stage T2D. It is also the major mechanistic link between them. This pathway is supposed to be driven by “adiposopathy” or “sick fat”, characterized by adipocyte hypertrophy and infiltration of pro-inflammatory cells (such as M1 macrophages and dendritic cells, etc.) into the expanding adipose tissue. The resulting heightened free-fatty-acid release, pro-inflammatory cytokine flux within the pathological adipose tissue, and ectopic fat deposition set the stage for heightened IR and lipotoxicity-mediated impairment in pancreatic β cells, together leading to overt T2D.
However, healthy obese individuals without complications of IR or T2D are frequently observed, indicating that the pathophysiology of IR and its transition to T2D is much more complicated than was earlier thought. Furthermore, the monogenetic disorder lipodystrophy, which is characterized by high IR status despite the substantial absence of adipose tissue, makes the pathophysiology of IR more elusive. At the molecular level, there might be some core set(s) of genes and gene-driven pathways that are shared by both IR and diabetes. These genes and pathways might also be responsible for the transition of IR to T2D. Therefore, identification of these genes would enable us not only to devise further clinical strategies to mitigate these diseases, but also to predict the transition from IR to T2D.
High-throughput gene expression methods such as DNA microarray and RNA-seq have proven useful in deciphering the pathophysiological processes occurring at the molecular level in various disease states. A huge number of gene expression datasets are now available in public databases like the Gene Expression Omnibus (GEO) and Array Express, which are routinely mined for the extraction of biological insights to further fundamental and clinical understandings. Problems associated with these studies include their low sample sizes due to ethical and other constraints, and the simultaneous measurements of tens of thousands of genes in these experiments. This high-dimensional nature makes the findings of these studies difficult to generalize to other related studies, which is a prerequisite for their clinical utility. One solution to this problem is the statistical meta-analysis of multiple gene expression datasets across related studies, which provides a robust set of metagenes that could be more confidently exploited for diagnostic and therapeutic purposes [7].
Various microarray meta-analysis studies have been conducted to decipher the molecular pathophysiology of IR and T2D. For example, Jung et al. (2018) conducted a meta-analysis of seven IR microarray studies and obtained drug signatures for two antidiabetic medicines, metformin and thiazolidinediones [8]. They also validated their signatures through cross-species analysis. Saxena et al. (2017–18) conducted combined system-level meta-analysis of IR and T2D datasets and highlighted adiposopathy—the inflammation of adipose tissue—as a central mechanism that leads to T2D [9,10]. However, one problem with large scale meta-analyses is that they yield a long list of metagenes which are difficult to interpret biologically, and it is practically difficult to exploit them further for clinical intervention due to the pleiotropic nature of the majority of genes in a genome.
In recent years, due to the availability of high volumes of clinical and molecular data, machine learning (ML) methods (also including deep learning methods) have begun to gain wide acceptance in the development of data-driven solutions to biomedical problems. These methods provide various mathematical and statistical models, which are provided with labeled datasets for supervised learning to enable them to predict the label or class of new data instances. Additionally, they can also be utilized for feature selection by identifying features or parameters in training datasets, which can significantly contribute to the models’ predictions. Machine learning has been used to develop predictive models for T2D using clinical features [11,12,13,14]; however, its use to predict T2D from gene-level features is relatively limited. Furthermore, the development of machine learning models based on IR gene expression could allow prognostication of the onset of T2D, and could be considered a novel approach with prospective clinical utility in T2D treatment.
It this study, we first attempted to obtain a robust IR metagene signature through meta-analysis of multiple microarray gene expression studies comprising samples of NGT and IR individuals.
Gene expression datasets are high-dimensional in nature and contain expression measurements of tens of thousands of genes across a limited number of samples. ML models based on these datasets may suffer problems with overfitting, which limits their predictive capability. To alleviate this problem, we first reduced the number of metagenes using an ML method—Least Absolute Shrinkage and Selection Operator (LASSO)—and obtained a core set of metagenes, hereafter termed key metagenes. These key metagenes were subsequently used to develop five baseline machine learning models based on LASSO, Support Vector Machine (SVM), eXtream Gradient Boosting (XGBoost), Random Forest (RF), and Artificial Neural Network (ANN).
We explored the potential of metagenes and key metagenes to explain the causation of T2D using GSEA-based enrichment analysis and disease-domain-specific tools, respectively. We also checked whether these key metagenes might represent the genetic basis of T2D, based on their presence in T2D and other glycemic traits identified in genome-wide association studies (GWASs). In addition to this direct evidence supporting the robustness of our ML models, we also checked for the presence of upstream transcription factors of these key metagenes in the GWASs as indirect evidence for the involvement of these genes in the genetic basis of T2D.
Our best-performing ML model demonstrated high prediction power in both cross-validation and test datasets. To sum up, the findings of present study endorse the use of our ML model in predicting the occurrence of T2D in IR individuals on the basis of the transcription signature of key metagenes. This ML model deserves further investigation of potential applications of these key metagenes in the prediction of transition from NGT to T2D in IR individuals, and of potential drug targets for the prevention and treatment of T2D.

2. Materials and Methods

2.1. Selection of Microarray Datasets and Meta-Analysis

The GEO database at the National Centre for Biotechnology Information was searched for gene expression studies including NGT and IR human tissue samples. As various insulin-responsive organs and tissues are involved in IR development, including pancreas, liver, skeletal muscle, kidneys, brain, small intestine, adipose tissue (subcutaneous and visceral), and peripheral blood mononuclear cells (PBMCs), we selected a total of nine gene expression datasets which profiled gene expression levels in these tissues (Table 1).
The GEOquery package [15] from R Bioconductor was used to download series matrix files for each selected study. Meta-analysis was carried out using the web-based tool Network Analyst [16]. To adjust for the batch effect among datasets, the ComBat function in the SVA R package was used. To derive meta-signatures, Fisher’s method was used on gene-level log-transformed p-values after adjustment of batch size. Metagenes were selected based on p-values < 0.05.

2.2. Enrichment Analysis of Metagenes

Computational validation of IR metagenes was conducted using gene set enrichment analysis [17] methods to assess whether these genes were enriched in phenotype-relevant biological processes. The GSEAPreranked method in GSEA allows the submission of metagenes ranked by their Z-scores. Enrichment analysis of metagenes was conducted using the Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) collections of the Molecular Signatures Database (MSigDB) [18] to biologically validate the obtained metasignature; the results of the enrichment analyses were later visualized using EnrichmentMap [19] in the Cytoscape [20] environment.

2.3. Selection of Key Metagenes and Their Biological Validation for T2D

The batch-effect-adjusted combined expression matrix was then filtered against the metagenes to obtain the IR meta-expression matrix. This matrix was subsequently used as the input for the ML method LASSO, which also performs feature selection and has been widely used in various biomedical studies to reduce numbers of genes [21,22,23,24]. The reduced set of metagenes obtained from LASSO was termed the “key metagenes”. As GO terms are not disease-specific, in order to explore the implications of these key metagenes in the disease domain, we used two tools—GLAD4U (Gene List Automatically Derived For You) [25] and DisGeNET [26]—to compile up-to-date gene–disease association information using text mining.

2.4. GWAS Evidence for Involvement of Key Metagenes in T2D

In order to further validate the roles of the 72 identified genes in diabetes, we checked the Type 2 Diabetes Knowledge Portal (https://t2d.hugeamp.org/; accesed on 23 November 2022) which houses genes that have been found to be associated with type 2 diabetes and other glycemic traits such as HOMA-IR and HOMA-B across multiple GWASs. The p-values for associations with phenotypes were calculated using the MAGMA (Multi-marker Analysis of GenoMic Annotation) method. In addition to finding GWAS-based direct associations of key metagenes, we also obtained a list of upstream regulators (transcription factors) of these genes using the Expression-to-Kinase (X2K) [27] webserver, which infers upstream regulatory networks from the signatures of differentially expressed genes through a combination of transcription factor enrichment analysis, protein–protein interaction network expansion, and kinase enrichment analysis.

2.5. Construction of Machine Learning Models and Their Performance Evaluation

The IR meta-expression matrix was further filtered for the key metagenes and the resulting matrix was then used to train five baseline ML models. SVM has been used to classify gene expression data for many years [28,29]; however, its use as a predictive model in the bioinformatics community has also begun to gain momentum in recent times [30]. XGBoost and RF are ensemble-based learning methods that assimilate multiple tree models to build a robust learner model. These methods have been used for the prediction and classification of gene expression data [31,32]. Deep learning models, particularly ANN, have also been used extensively in bioinformatics [33,34,35,36,37,38,39]. Details of the ML methods and their tuning parameters are presented in the Supplementary File.
The GEO dataset GSE64577, comprising 64 human adipose tissue samples, was used as a test dataset to check the performance of the ML models. These samples were drawn from a Mexican-American population under the Veterans Administration Genetic Epidemiology Study (VAGES) [40] and included 38 T2D and 26 NGT individuals as per their fasting plasma glucose level (T2D > 100 mg/dL).
We used five-fold cross-validation and the performance of the ML models was evaluated using the following parameters: accuracy (A), precession (p), recall (r), sensitivity (Sen), specificity (Spe), F1 score, and Matthews correlation coefficient (MCC) from the confusion matrix for both the validation and the test dataset. These values were calculated using the indicators TP (true positive), FP (false positive), FN (false negative), and TN (true negative).
A = T P + T N T P + T N + F P + F N
p = T P T P + F P
r = T P T P + F N
F 1 = 2 × p × r p + r
M C C = T P T N F P F N T P + F P T P + T N T N + F P T N + F N
where TP is the number of IR/T2D individuals predicted to be IR/T2D by the model; FN is the number of IR/T2D individuals incorrectly predicted to be NGT by the model; TN is the number of NGT individuals predicted to be NGT by the model; and FP is the number of NGT individuals incorrectly predicted to be IR/T2D by the model.

3. Results

NetworkAnalyst identified an IR metasignature of 2574 genes (p < 0.01), containing 1488 upregulated (Z-score > 0) and 1086 downregulated (Z-score < 0) genes. The derived IR meta-expression matrix for these metagenes had 367 expression datasets consisting of 182 NGT and 185 IR samples.
GSEA-based functional analysis of the IR datasets showed enrichment in various GO terms such as “RNA metabolism”, “chromosome”, “transmembrane transport”, and “ubiquitin protein activity” (Figure 1a). The majority of these processes were found to be downregulated, likely due to weak insulin signaling in peripheral tissues in the setting of insulin resistance. Enriched KEGG pathways were “oxidative phosphorylation”, “type 2 diabetes”, and “ribosome” (Figure 1b). Enrichment of these pathways therefore provides biological validation of the IR signature.
A total of 72 key metagenes were obtained from the meta-signature through LASSO-based feature selection. GLAD4U analysis showed robust enrichment of various diabetic terms (Figure 2) such as “diabetes mellitus”, “diabetes mellitus type 2”, “diabetes mellitus type 2 and obesity”, “diabetic ketoacidosis”, etc. (marked red). Various cardiovascular complication terms were also found, such as “ventricular premature complexes”, “shortened QT interval”, “ventricular fibrillation”, etc. (marked blue). Traditionally, diabetes mellitus is considered to be a disorder of glucose metabolism, whereas the roles of the disorders of lipid metabolism and adipose tissue dysfunction in the pathogenesis of atherosclerotic vascular disease are well known. However, this classical picture is changing rapidly due to appreciation of the fact that IR driven by adipose tissue dysfunction induces beta-cell dysfunction both directly via its exhaustion and indirectly via lipotoxicity. Therefore, adipose tissue also plays a major role in the pathogenesis of T2D.
Additionally, cardiovascular diseases and phenotypes and terms related to renal and hepatic complications such as kidney diseases, fatty liver, carcinoma hepatocellular, urogenital neoplasm, and urination disorder were also found to be enriched. Six of the key metagenes—INSR, MAP3K5, NDUFB8, NDUFS1, SDHB, and UQCRC2—are involved in nonalcoholic fatty liver disease [41] (hsa04932), which is caused by a defect in insulin suppression of free fatty acid (FAA) disposal due to the induction of insulin resistance. Furthermore, a reduction in nitric oxide production due to IR has been reported to inhibit bladder smooth muscle cell growth and consequently lead to urination disorder [42].
DisGeNET analysis highlighted that 27 of the IR key metagenes (~37%)—ALDH6A1, APOB, ARID5B, ATP2B2, ATXN1, B2M, CAT, CFB, CHI3L1, COL8A1, EIF2AK2, FHL2, GATM, HDAC4, HIPK3, HOMER1, IGFBP5, INSR, MAP3K5, MMP9, NDUFS1, RASSF7, REV3L, SLC19A2, UROD, WAS, and ZEB1—were found to be associated with T2D or its related variants such as alloxan, autoimmune, monogenic, neonatal, ketosis-prone, sudden-onset, insipidus, fibrocalculous pancreatic, phosphate, gestational, post-transplant, lipoatrophic, streptozotocin, brittle, and maturity-onset diabetes, among others (Figure 3).
To assess the pathway-based functional assignment for these genes, we selected the human collection of WikiPathways and found enrichment of the following 15 pathways (p < 0.05), including “insulin signaling”, which has a conspicuous role in the IR phenotype (Table 2).
The pathway “FTO obesity variant mechanism” is relevant to the T2D phenotype as SNPs in the Fat mass and Obesity-associated (FTO) gene have been found to be associated with adiposity and risk of obesity in multiple populations [43]. Due to involvement of FTO in energy homeostasis, it could link IR to T2D. Disturbance in the “folate metabolism” pathway has been reported to induce glucose and lipid metabolism disorders in animal studies [44], suggesting its involvement in IR. The pathway “selenium micronutrient network” also influences adipocyte physiology and modifies the risk of developing T2D. In one study, high Se and hs-CRP concentrations were found to be associated with high HbA1c levels in various BMI groups [45]. Another pathway, “vitamin B12 metabolism”, is also relevant as high prevalence of low B12 levels have been shown in European (27%) and South Indian (32%) patients with type 2 diabetes (T2D) [46]. Angiopoietin like 8 (ANGPTL8) is involved in the regulation of lipid metabolic processes and triglyceride homeostasis, so the term “Angiopoietin Like Protein 8 regulatory pathway” is also relevant to IR pathophysiology. AGE/RAGE signaling increases oxidative stress to promote diabetes-mediated vascular calcification through activation of Nox-1 and decreased expression of SOD-1. Different pharmacological interventions are underway to regulate the AGE/RAGE pathway to decrease the severity of this diabetic complication [47]. Enriched pathways therefore showed the biological validity of the key metagenes.
Key metagenes were further checked against T2D GWASs to ascertain their implications in the genetic architecture of the disease. Seventeen genes were mapped to genetic loci associated with T2D in GWASs. These genes were CFB, PRSS3, ARID5B, REV3L, CROCC, HOMER1, NDUFB8, DDX17, DCAF8, EIF2AK2, ZEB1, ATF6B, UROD, MAP3K5, MMP9, STK24, and WNT4; three genes with HOMA-IR-related loci, MATK, FHL2, and PHKB; and two with HOMA-B-related loci, HS3ST1 and ALDH6A1. The gene INSR has been found to show association with all the three traits.
The human genome consists of 24,500 protein-coding genes, of which 4049 have been identified to have T2D GWAS signals (p < 0.05), so the probability of picking a T2D gene by chance would be 0.17 (i.e., 4049/24,500).
To check the statistical robustness of our key metagenes, we estimated the probability of finding 17 T2D genes by chance across 72 key metagenes using a binomial test (n = 72, x = 17, p = 0.17) and obtained a low probability (p = 0.03), and thus concluded that our key metagene set was nonrandomly enriched with T2D genes.
We carried out further analysis by mapping of the regulatory elements of the key metagenome using the X2K webserver. The webserver analyzed the key metagenome of IR and reported 83 transcription factors (including 2 key metagenes). Interestingly, 42 of these transcription factors were found to be associated with T2D and its related traits in GWASs. In other words, we could provide in silico evidence of the mapping of key metagenome genes to T2D and related loci identified in GWASs and shed light on their potential role in the genome-to-phenome trajectory of T2D. Therefore, our ML model based on expression profiles of these 72 key metagenes could be considered highly robust for predicting diabetes.
All five machine learning models demonstrated good predictive capabilities upon five-fold cross-validation. Detailed information of the various evaluation indicators is shown in Table 3.
The classification accuracy was visualized using receiver operator characteristic (ROC) plots for all the models and all showed an area under the curve (AUC) > 75%, with ANN returning the maximum AUC of 95% (Figure 4).
The best-performing model, ANN, was subsequently evaluated for its performance on the test dataset. Out of 64 samples (38 T2D: 26 NGT), it predicted 28 samples as true positive, 21 samples as true negative, 10 samples as false negative, and 5 samples as false positive, and showed 73% sensitivity, 80% specificity, and an overall F1 score of 0.78 (overall accuracy 76%). Despite the high biological noise in the gene expression data, differences among microarray platforms, and also marked differences in the molecular pathologies of IR and T2D due to chronic hyperglycemia in the latter condition, our IR-based ML model achieved fairly good accuracy in differentiating individuals with T2D from non-T2D people.

4. Discussion

In this study, we developed a machine learning model that could differentiate an individual with T2D from a NGT individual on the basis of expression of a key metagenome of IR in insulin-responsive tissues, particularly adipose tissue, with almost 70 to 75 percent accuracy in a cross-sectional study. This key metagenome was created via the meta-analysis of publicly available databases of transcription profiles of various insulin-responsive tissues obtained from IR non-T2D individuals with diverse ethnic backgrounds. The genes in this key metagenome enrich several known pathways of not only diabetes mellitus, but also cardiovascular disease, specifically cardiac arrhythmias, kidney diseases, fatty liver, carcinoma hepatocellular, urogenital neoplasm, urination disorder, etc. The genetic loci regulating the expressions of most of these genes are mapped to genomic regions identified in GWASs as associated with diabetes and related traits. Therefore, these findings suggest that the key metagenome of insulin resistance (IR) also plays a functional role in the genetic susceptibility to T2D via the enrichment of several known pathophysiological molecular pathways of T2D in metabolically active tissues. As the IR key metagenome can be mapped to GWAS-identified genetic loci of diabetes, this metagenome is unlikely to be merely an acquired pathophysiological adaptation to impaired insulin signal transduction; rather, it is a primary genetic defect.
There are several other implications of these findings. For example, T2D and IR are complex traits and their genetic susceptibility is generally believed to be omnigenic in nature. However, in this study, out of almost 30,000 genes in the human genome, only 72 genes were found to be of functional significance at the transcriptome level to the manifestation of genetic susceptibility to T2D via the known molecular pathways of IR. However, a limitation of the present study is that the source of this cross-sectional meta-analysis was data from publicly available depositories, and the mapping and functional roles were imputed in silico. Though we have estimated that our metagenes are associated with T2D GWAS signals, due to the nature of our study being cross-sectional, these gene expression results might not be free from reverse causation. Therefore, there is a need for validation of the role of this metagenome through deciphering of the genome-to-phenome trajectory of diabetes using replication studies in different insulin-responsive tissues.
If this metagenome proves to be a gold standard transcriptomic signature of IR in diabetes, then it would also deserve further inquiry regarding clinical applications like diagnostic molecular markers and drug target discovery. However, the most important and challenging question is as follows: How should these metagenes obtained from insulin-responsive tissues in research settings be measured/assessed in clinical practice? Therefore, from the clinical point of view, there is a need to identify not only the clinical, biochemical, and radiological parameters that show associations with the key metagenome, but also the gene sequence polymorphisms and the circulatory biomarkers associated with this metagenome.
IR and beta-cell dysfunction are two major pathophysiological defects in T2D. IR has traditionally been considered to be caused by environmental factors, whereas the genetic susceptibility to develop T2D possibly plays a major role in beta-cell failure. However, identification of the key metagenome of IR, its imputed mapping to T2D GWAS loci, and functional enrichment of known T2D pathways and validation in adipose tissues of T2D individuals suggests a role of genetic factors in driving IR. Additionally, this metagenome is derived from meta-analysis of several insulin-responsive tissue transcription profiles. Therefore, it points towards a common molecular thread shared by most metabolically active insulin-responsive tissues and it could serve as “the cardinal molecular pathology” of IR, T2D, and MetS. However, this concept deserves further investigation.
Another important finding of this study is that the genes in the key metagenome of IR enrich both diabetes- and cardiac-arrhythmia-related pathways. Therefore, they shed light on a common functional genomic defect underlying these two diseases. The relationship between diabetes and cardiac arrhythmias is complex and it is not yet fully understood. However, relationships between diabetes and arrhythmias have been reported [48]. There are several potential mechanisms of arrhythmia in diabetes, such as increased blood glucose levels, glucose fluctuation, hypoglycemia, autonomic dysfunction, alterations in the architecture of the heart including fibrosis, fat deposition, hypertrophy, etc. In addition, the finding of a key metagenome linking IR with cardiac arrhythmias in the present study points towards another novel molecular mechanism shared by these two diseases.

5. Conclusions

In conclusion, a machine learning model trained with the key metagenome of IR can differentiate individuals with T2D from NGT with moderate accuracy on the basis of transcription profiles of adipose and other insulin-responsive tissues. The mapping of this metagenome to GWASs identified loci of diabetes, and the enrichment of known molecular pathways of diabetes suggest a primary role of this metagenome in the pathogenesis of diabetes rather than merely a pathophysiological response to impaired insulin signaling. Therefore, this metagenome deserves further investigation in terms of the cardinal molecular pathological defects/pathways underlying both IR and diabetes.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/biom13030432/s1.

Author Contributions

S.K.M. conceptualized the experimental design, A.S. and P.P. conducted the bioinformatics and machine learning part of study. S.K.M., P.T., N.M. and A.S. analyzed the results. All authors reviewed the manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

Correspondence and requests for data and materials should be addressed to S.K.M.

Acknowledgments

This research received no specific grant from any funding agency in the public, commercial, or not-for-profit sectors. The authors express gratitude towards the editor and anonymous reviewers for constructive suggestions and valuable comments. AS and PP are grateful to GLA University, Mathura, India, and NM and SM to SMS Medical College & Hospital, Jaipur, India.

Conflicts of Interest

The authors declare no competing interests.

References

  1. Bonora, E.; Kiechl, S.; Willeit, J.; Oberhollenzer, F.; Egger, G.; Targher, G.; Alberiche, M.; Bonadonna, R.C.; Muggeo, M. Prevalence of insulin resistance in metabolic disorders: The Bruneck Study. Diabetes 1998, 47, 1643–1649. [Google Scholar] [CrossRef] [PubMed]
  2. World Health Organization. Global Report on Diabetes. 2016. Available online: https://www.who.int/publications/i/item/9789241565257 (accessed on 1 October 2019).
  3. Meigs, J.B.; Cupples, L.A.; Wilson, P.W. Parental transmission of type 2 diabetes: The Framingham Offspring Study. Diabetes 2000, 49, 2201–2207. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Poulsen, P.; Kyvik, K.O.; Vaag, A.; Beck-Nielsen, H. Heritability of type II (non-insulin-dependent) diabetes mellitus and ab-normal glucose tolerance—A population-based twin study. Diabetologia. 1999, 42, 139–145. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Fuchsberger, C.; Flannick, J.; Teslovich, T.M.; Mahajan, A.; Agarwala, V.; Gaulton, K.J.; Ma, C.; Fontanillas, P.; Moutsianas, L.; McCarthy, D.J.; et al. The genetic architecture of type 2 diabetes. Nature 2016, 536, 41–47. [Google Scholar] [CrossRef] [Green Version]
  6. Cantley, J.; Ashcroft, F.M. Q&A: Insulin secretion and type 2 diabetes: Why do β-cells fail? BMC Biol. 2015, 13, 1–7. [Google Scholar] [CrossRef] [Green Version]
  7. Saxena, A. Bioinformatics of Meta-analyses of Genomic Data. In Bioinformatics and Human Genomics Research; CRC Press: Boca Raton, FL, USA, 2021; pp. 331–334. [Google Scholar]
  8. Jung, J.; Kim, G.W.; Lee, W.; Mok, C.; Chung, S.H.; Jang, W. Meta- and cross-species analyses of insulin resistance based on gene expression datasets in human white adipose tissues. Sci. Rep. 2018, 8, 1–13. [Google Scholar] [CrossRef] [Green Version]
  9. Saxena, A.; Sachin, K.; Bhatia, A.K. System Level Meta-analysis of Microarray Datasets for Elucidation of Diabetes Mellitus Pathobiology. Curr. Genom. 2017, 18, 298–304. [Google Scholar] [CrossRef] [Green Version]
  10. Saxena, A.; Sachin, K. A Network Biology Approach for Assessing the Role of Pathologic Adipose Tissues in Insulin Re-sistance Using Meta-analysis of Microarray Datasets. Curr. Genom. 2018, 19, 630–636. [Google Scholar] [CrossRef] [PubMed]
  11. Kavakiotis, I.; Tsave, O.; Salifoglou, A.; Maglaveras, N.; Vlahavas, I.; Chouvarda, I. Machine Learning and Data Mining Methods in Diabetes Research. Comput. Struct. Biotechnol. J. 2017, 15, 104–116. [Google Scholar] [CrossRef]
  12. Zou, Q.; Qu, K.; Luo, Y.; Yin, D.; Ju, Y.; Tang, H. Predicting Diabetes Mellitus with Machine Learning Techniques. Front. Genet. 2018, 9, 515. [Google Scholar] [CrossRef]
  13. Tigga, N.P.; Garg, S. Prediction of Type 2 Diabetes using Machine Learning Classification Methods. Procedia Comput. Sci. 2020, 167, 706–716. [Google Scholar] [CrossRef]
  14. Kaur, H.; Kumari, V. Predictive modelling and analytics for diabetes using a machine learning approach. Appl. Comput. Inform. 2020. ahead of print. [Google Scholar] [CrossRef]
  15. Davis, S.; Meltzer, P.S. GEOquery: A bridge between the Gene Expression Omnibus (GEO) and BioConductor. Bioinformatics 2007, 23, 1846–1847. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  16. Zhou, G.; Soufan, O.; Ewald, J.; Hancock, R.E.W.; Basu, N.; Xia, J. NetworkAnalyst 3.0: A visual analytics platform for comprehensive gene expression profiling and meta-analysis. Nucleic Acids Res. 2019, 47, W234–W241. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  17. Subramanian, A.; Tamayo, P.; Mootha, V.K.; Mukherjee, S.; Ebert, B.L.; Gillette, M.A.; Paulovich, A.; Pomeroy, S.L.; Golub, T.R.; Lander, E.S.; et al. Gene set enrichment analysis: A knowledge-based approach for interpreting genome-wide expression pro-files. Proc. Natl. Acad. Sci. USA 2005, 102, 15545–15550. [Google Scholar] [CrossRef] [Green Version]
  18. Liberzon, A.; Subramanian, A.; Pinchback, R.; Thorvaldsdóttir, H.; Tamayo, P.; Mesirov, J.P. Molecular signatures database (MSigDB) 3.0. Bioinformatics 2011, 27, 1739–1740. [Google Scholar] [CrossRef]
  19. Merico, D.; Isserlin, R.; Stueker, O.; Emili, A.; Bader, G.D. Enrichment Map: A Network-Based Method for Gene-Set Enrichment Visualization and Interpretation. PLOS ONE 2010, 5, e13984. [Google Scholar] [CrossRef] [Green Version]
  20. Shannon, P.; Markiel, A.; Ozier, O.; Baliga, N.S.; Wang, J.T.; Ramage, D.; Amin, N.; Schwikowski, B.; Ideker, T. Cytoscape: A software environment for integrated models of Biomolecular Interaction Networks. Genome Res. 2003, 13, 2498–2504. [Google Scholar] [CrossRef]
  21. Xiong, Y.; Ling, Q.-H.; Han, F.; Liu, Q.-H. An efficient gene selection method for microarray data based on LASSO and BPSO. BMC Bioinform. 2019, 20, 1–13. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  22. Kaneko, S.; Hirakawa, A.; Hamada, C. Enhancing the Lasso Approach for Developing a Survival Prediction Model Based on Gene Expression Data. Comput. Math. Methods Med. 2015, 2015, 1–7. [Google Scholar] [CrossRef]
  23. Fontanarosa, J.B.; Dai, Y. Using LASSO regression to detect predictive aggregate effects in genetic studies. BMC Proc. 2011, 5, S69. [Google Scholar] [CrossRef] [Green Version]
  24. Xiao, J.; Wang, R.; Cai, X.; Ye, Z. Coupling of co-expression network analysis and machine learning validation unearthed po-tential key genes involved in rheumatoid arthritis. Front. Genet. 2021, 12, 9. [Google Scholar] [CrossRef]
  25. Jourquin, J.; Duncan, D.; Shi, Z.; Zhang, B. GLAD4U: Deriving and prioritizing gene lists from PubMed literature. BMC Genom. 2012, 13 (Suppl. 8), S20. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  26. Piñero, J.; Bravo, À.; Queralt-Rosinach, N.; Gutiérrez-Sacristán, A.; Deu-Pons, J.; Centeno, E.; García-García, J.; Sanz, F.; Furlong, L.I. DisGeNET: A comprehensive platform integrating information on human disease-associated genes and variants. Nucleic Acids Res. 2017, 45, D833–D839. [Google Scholar] [CrossRef]
  27. Clarke, D.; Kuleshov, M.V.; Schilder, B.M.; Torre, D.; E Duffy, M.; Keenan, A.B.; Lachmann, A.; Feldmann, A.S.; Gundersen, G.W.; Silverstein, M.C.; et al. eXpression2Kinases (X2K) Web: Linking expression signatures to upstream cell signaling networks. Nucleic Acids Res. 2018, 46, W171–W179. [Google Scholar] [CrossRef] [Green Version]
  28. Brown, M.P.; Grundy, W.N.; Lin, D.; Cristianini, N.; Sugnet, C.W.; Furey, T.S.; Ares, M.; Haussler, D. Knowledge-based analysis of microarray gene expression data by using support vector machines. Proc. Natl. Acad. Sci. USA 2000, 97, 262–267. [Google Scholar] [CrossRef] [Green Version]
  29. George, G.; Raj, V.C. Review on feature selection techniques and the impact of SVM for cancer classification using gene ex-pression profile. arXiv 2011, arXiv:1109.1062. [Google Scholar]
  30. Chen, Z.; Li, J.; Wei, L.; Xu, W.; Shi, Y. Multiple-kernel SVM based multiple-task oriented data mining system for gene expression data analysis. Expert Syst. Appl. 2011, 38, 12151–12159. [Google Scholar] [CrossRef]
  31. Li, W.; Yin, Y.; Quan, X.; Zhang, H. Gene Expression Value Prediction Based on XGBoost Algorithm. Front. Genet. 2019, 10, 1077. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  32. Okun, O.; Priisalu, H. Random Forest for Gene Expression Based Cancer Classification: Overlooked Issues. In Proceedings of the Iberian Conference on Pattern Recognition and Image Analysis, Girona, Spain, 6–8 June 2007; pp. 483–490. [Google Scholar]
  33. Chen, Y.-C.; Ke, W.-C.; Chiu, H.-W. Risk classification of cancer survival using ANN with gene expression data from multiple laboratories. Comput. Biol. Med. 2014, 48, 1–7. [Google Scholar] [CrossRef]
  34. Khan, J.; Wei, J.S.; Ringner, M.; Saal, L.H.; Ladanyi, M.; Westermann, F.; Berthold, F.; Schwab, M.; Antonescu, C.R.; Peterson, C.; et al. Classification and diagnostic prediction of cancers using gene expression profiling and artificial neural networks. Nat. Med. 2001, 7, 673–679. [Google Scholar] [CrossRef] [PubMed]
  35. Vohradsky, J. Neural network model of gene expression. FASEB J. 2001, 15, 846–854. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  36. Sharma, A. Skin Burn Detection using Feature Extraction. Ann. Rom. Soc. Cell Biol. 2021, 25, 1656–1662. [Google Scholar]
  37. Yadav, D.P.; Sharma, A.; Singh, M.; Goyal, A. Feature Extraction Based Machine Learning for Human Burn Diagnosis from Burn Images. IEEE J. Transl. Eng. Heal. Med. 2019, 7, 1–7. [Google Scholar] [CrossRef]
  38. Yadav, D.P.; Rathor, S. Bone fracture detection and classification using deep learning approach. In Proceedings of the 2020 International Con-ference on Power Electronics & IoT Applications in Renewable Energy and its Control (PARC), Mathura, India, 28–29 February 2020; pp. 282–285. [Google Scholar]
  39. Yadav, D.P.; Saini, P.; Mittal, P. Feature Optimization Based Heart Disease Prediction using Machine Learning. In Proceedings of the 2021 5th International Conference on Information Systems and Computer Networks (ISCON), Mathura, India, 22–23 October 2021; pp. 1–5. [Google Scholar]
  40. Winnier, D.A.; Fourcaudot, M.; Norton, L.; Abdul-Ghani, M.A.; Hu, S.L.; Farook, V.S.; Coletta, D.K.; Kumar, S.; Puppala, S.; Chittoor, G.; et al. Transcriptomic Identification of ADH1B as a Novel Candidate Gene for Obesity and Insulin Resistance in Human Adipose Tissue in Mexican Americans from the Veterans Administration Genetic Epidemiology Study (VAGES). PLoS ONE 2015, 10, e0119941. [Google Scholar] [CrossRef]
  41. Ogata, H.; Goto, S.; Sato, K.; Fujibuchi, W.; Bono, H.; Kanehisa, M. KEGG: Kyoto Encyclopedia of Genes and Genomes. Nucleic Acids Res. 2000, 28, 29–34. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  42. Uzun, H.; Yılmaz, A.; Kemik, A.; Zorba, O.U.; Kalkan, M. Association of Insulin Resistance with Overactive Bladder in Female Patients. Int. Neurourol. J. 2012, 16, 181–186. [Google Scholar] [CrossRef] [Green Version]
  43. Fawcett, K.A.; Barroso, I. The genetics of obesity: FTO leads the way. Trends Genet. 2010, 26, 266–274. [Google Scholar] [CrossRef] [Green Version]
  44. Zhao, M.; Yuan, M.M.; Yuan, L.; Huang, L.L.; Liao, J.H.; Yu, X.L.; Su, C.; Chen, Y.H.; Yang, Y.Y.; Yu, H.; et al. Chronic folate deficiency induces glucose and lipid metabolism disorders and subsequent cognitive dysfunction in mice. PLoS ONE 2018, 13, e0202910. [Google Scholar] [CrossRef]
  45. Wongdokmai, R.; Shantavasinkul, P.C.; Chanprasertyothin, S.; Panpunuan, P.; Matchariyakul, D.; Sritara, P.; Sirivarasai, J. The Involvement of Selenium in Type 2 Diabetes Development Related to Obesity and Low Grade Inflammation. Diabetes Metab. Syndr. Obes. Targets Ther. 2021, 14, 1669–1680. [Google Scholar] [CrossRef]
  46. Boachie, J.; Adaikalakoteswari, A.; Samavat, J.; Saravanan, P. Low Vitamin B12 and Lipid Metabolism: Evidence from Pre-Clinical and Clinical Studies. Nutrients 2020, 12, 1925. [Google Scholar] [CrossRef] [PubMed]
  47. Kay, A.M.; Simpson, C.L.; Stewart, J.A. The Role of AGE/RAGE Signaling in Diabetes-Mediated Vascular Calcification. J. Diabetes Res. 2016, 2016, 1–8. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  48. Grisanti, L.A. Diabetes and Arrhythmias: Pathophysiology, Mechanisms and Therapeutic Outcomes. Front. Physiol. 2018, 9, 1669. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Enrichment analysis of (a) IR metagenes using Z-scores of the meta-analysis and Gene Ontology (GO) gene sets. (b) GSEA analysis of IR metagenes using KEGG in the Molecular Signatures Database (MSigDB). NES, normalized enrichment score; FDR, false discovery rate.
Figure 1. Enrichment analysis of (a) IR metagenes using Z-scores of the meta-analysis and Gene Ontology (GO) gene sets. (b) GSEA analysis of IR metagenes using KEGG in the Molecular Signatures Database (MSigDB). NES, normalized enrichment score; FDR, false discovery rate.
Biomolecules 13 00432 g001
Figure 2. Enrichment of disease terms for IR key metagenes using GLAD4U with false discovery rate <0.05 and number of genes in a category >4.
Figure 2. Enrichment of disease terms for IR key metagenes using GLAD4U with false discovery rate <0.05 and number of genes in a category >4.
Biomolecules 13 00432 g002
Figure 3. Disease–gene network for IR key metagenes for type 2 diabetes and related disorders, created using DisGeNET.
Figure 3. Disease–gene network for IR key metagenes for type 2 diabetes and related disorders, created using DisGeNET.
Biomolecules 13 00432 g003
Figure 4. Receiver operator characteristic curves for (a) Least Absolute Shrinkage and Selection Operator (LASSO), (b) Random Forest (RF), (c) Support Vector Machine (SVM), (d) eXtream Gradient Boosting (XGBoost), and (e) Artificial Neural Network (ANN), and (f) confusion matrix for best-performing ANN model.
Figure 4. Receiver operator characteristic curves for (a) Least Absolute Shrinkage and Selection Operator (LASSO), (b) Random Forest (RF), (c) Support Vector Machine (SVM), (d) eXtream Gradient Boosting (XGBoost), and (e) Artificial Neural Network (ANN), and (f) confusion matrix for best-performing ANN model.
Biomolecules 13 00432 g004
Table 1. Microarray samples used in meta-analysis for IR.
Table 1. Microarray samples used in meta-analysis for IR.
S. NoGEO SeriesTissuePlace of StudyDisease Phenotype
Insulin Resistance (IR)Normal Glucose Tolerance (NGT)
1GSE6798Skeletal muscleDepartment of Hematology in Roskilde Hospital, Roskilde, Denmark.IR = 1613
2GSE15773Subcutaneous adipose tissue (SAT) and visceral adipose tissue (VAT)Department of Molecular Medicine at University of Massachusetts, Worcester, MA, USAIR = 4 (SAT)
IR = 5 (VAT)
5 (SAT)
5 (VAT)
3GSE20950Subcutaneous adipose tissue (SAT) and visceral adipose tissue (VAT)Department of Molecular Medicine at University of Massachusetts, Worcester, MA, USA.IR = 9 (SAT)
IR = 10 (VAT)
10 (SAT)
10 (VAT)
4GSE22309Skeletal muscleDepartment of Biostatistics at University of Alabama, Birmingham, AL, USAIR = 2020
6GSE26637Subcutaneous adipose tissueDepartment of Institute for Molecular Medicine Finland (FIMM), University of Helsinki, Helsinki, Finland.IR = 5 (Fasting)
IR = 5 (Hyperinsulinemia)
5 (Fasting)
5 (Hyperinsulinemia)
7GSE34526Granulosa cellsDepartment of Zoology at University of Delhi, Delhi, India.IR = 16 (PCOS)12 (PCOS)
8GSE36297Vastus lateralis muscleDepartment of Hematology, Roskilde Hospital, Roskilde, Denmark.IR = 610
9GSE64567Fasted subcutaneous abdominal adipose tissue (FAT)Department of Medicine at University of Texas, Health Sciences Center at San Antonio, San Antonio, TX, USA.IR = 3826
Table 2. List of pathways enriched by the IR key metagenes against the human collection of WikiPathways with p < 0.05.
Table 2. List of pathways enriched by the IR key metagenes against the human collection of WikiPathways with p < 0.05.
S. No.Enrichment FDRGenes in ListTotal GenesFunctional Category
14.4E-0328FTO Obesity Variant Mechanism
24.4E-034103Electron Transport Chain
31.1E-02366Folate Metabolism
41.7E-02385Selenium Micronutrient Network
52.8E-023110DNA Damage Response (only ATM dependent)
62.9E-02238Amyotrophic lateral sclerosis (ALS)
73.0E-02250Vitamin B12 Metabolism
83.0E-023132Angiopoietin Like Protein 8 Regulatory Pathway
93.0E-02252Translation Factors
103.0E-023155Myometrial Relaxation and Contraction Pathways
113.0E-023160Insulin Signaling
123.0E-02245ATM Signaling Network in Development and Disease
133.0E-023159Epithelial to mesenchymal transition in colorectal cancer
143.0E-02260Oxidative phosphorylation
153.1E-02266AGE/RAGE pathway
Table 3. The evaluation indicators of the five machine learning models.
Table 3. The evaluation indicators of the five machine learning models.
Evaluation IndicatorsLASSOSVMXGBoostRandom ForestANN
False Positive Rate | Type I error0.220.200.080.300.04
False Negative Rate | Type II error0.270.290.210.190.05
True Negative Rate | Specificity0.770.790.910.690.95
Negative Predictive Value0.690.680.810.730.95
False Discovery Rate0.190.180.090.230.051
True Positive Rate | Recall | Sensitivity0.720.700.780.800.94
Positive Predictive Value | Precision0.800.810.900.760.94
Accuracy0.740.740.850.750.95
F1 Score0.760.750.840.780.94
Matthews Correlation Coefficient MCC0.490.500.710.500.9
ROC AUC score0.750.750.850.750.95
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Saxena, A.; Mathur, N.; Pathak, P.; Tiwari, P.; Mathur, S.K. Machine Learning Model Based on Insulin Resistance Metagenes Underpins Genetic Basis of Type 2 Diabetes. Biomolecules 2023, 13, 432. https://doi.org/10.3390/biom13030432

AMA Style

Saxena A, Mathur N, Pathak P, Tiwari P, Mathur SK. Machine Learning Model Based on Insulin Resistance Metagenes Underpins Genetic Basis of Type 2 Diabetes. Biomolecules. 2023; 13(3):432. https://doi.org/10.3390/biom13030432

Chicago/Turabian Style

Saxena, Aditya, Nitish Mathur, Pooja Pathak, Pradeep Tiwari, and Sandeep Kumar Mathur. 2023. "Machine Learning Model Based on Insulin Resistance Metagenes Underpins Genetic Basis of Type 2 Diabetes" Biomolecules 13, no. 3: 432. https://doi.org/10.3390/biom13030432

APA Style

Saxena, A., Mathur, N., Pathak, P., Tiwari, P., & Mathur, S. K. (2023). Machine Learning Model Based on Insulin Resistance Metagenes Underpins Genetic Basis of Type 2 Diabetes. Biomolecules, 13(3), 432. https://doi.org/10.3390/biom13030432

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop