ORIGINAL RESEARCH article
In-Silico Study of Immune System Associated Genes in Case of Type-2 Diabetes With Insulin Action and Resistance, and/or Obesity
- 1Department of Clinical Biochemistry, Faculty of Medicine, King Abdulaziz University, Jeddah, Saudi Arabia
- 2Stem Cell Unit, Department of Medical Laboratory Technology, Faculty of Applied Medical Sciences, King Fahd Medical Research Center, King Abdulaziz University, Jeddah, Saudi Arabia
- 3Cancer and Mutagenesis Unit, Department of Biochemistry, Cancer Metabolism and Epigenetic Unit, Faculty of Science, King Fahd Medical Research Center, King Abdulaziz University, Jeddah, Saudi Arabia
- 4SciLifeLab, Department of Oncology and Pathology, Karolinska Institutet, Stockholm, Sweden
Type-2 diabetes and obesity are among the leading human diseases and highly complex in terms of diagnostic and therapeutic approaches and are among the most frequent and highly complex and heterogeneous in nature. Based on epidemiological evidence, it is known that the patients suffering from obesity are considered to be at a significantly higher risk of type-2 diabetes. There are several pieces of evidence that support the hypothesis that these diseases interlinked and obesity may aggravate the risk(s) of type-2 diabetes. Multi-level unwanted alterations such as (epi-) genetic alterations, changes at the transcriptional level, and altered signaling pathways (receptor, cytoplasmic, and nuclear level) are the major sources that promote several complex diseases, and such a heterogeneous level of complexity is considered as a major barrier in the development of therapeutics. With so many known challenges, it is critical to understand the relationships and the shared causes between type-2 diabetes and obesity, and these are difficult to unravel and understand. For this purpose, we have selected publicly available datasets of gene expression for obesity and type-2 diabetes, have unraveled the genes and the pathways associated with the immune system, and have also focused on the T-cell signaling pathway and its components. We have applied a simplified computational approach to understanding differential gene expression and patterns and the enriched pathways for obesity and type-2 diabetes. Furthermore, we have also analyzed genes by using network-level understanding. In the analysis, we observe that there are fewer genes that are commonly differentially expressed while a comparatively higher number of pathways are shared between them. There are only 4 pathways that are associated with the immune system in case of obesity and 10 immune-associated pathways in case of type-2 diabetes, and, among them, only 2 pathways are commonly altered. Furthermore, we have presented SPNS1, PTPN6, CD247, FOS, and PIK3R5 as the overexpressed genes, which are the direct components of TCR signaling.
Type 2 diabetes (T2D) is a global epidemic that is strongly correlated with obesity (1, 2). Adipose tissue inflammation associated with obesity is known as a major cause for decreased insulin sensitivity in the case of T2D (1, 3, 4). There are a number of works where the crosstalk between the immune system and metabolism have been presented alongside evidence that supports the hypothesis that these diseases are interlinked and that obesity may aggravate the risk(s) of T2D (5–8). In such cases, it is critical to understand multi-level unwanted alterations such as (epi-)genetic alterations, changes at the transcriptional level, and altered signaling pathways (receptor, cytoplasmic, and nuclear level). These are the major sources that promote a number of complex diseases and such heterogeneous complexities are considered the major barrier in the development of therapeutic approaches. With such diverse challenges, it is crucial to understand and unravel the relationships and the common causes between type-2 diabetes and obesity (9–13).
So far, there are several studies where a number of factors have been investigated to understand the biological mechanism of obesity and T2D and the relationship between the two (1, 2, 14–16). T2D is mainly characterized by a progressive status of chronic disease, and inflammation of low-grade disease and inflammatory responses are triggered by many factors. The well-known factors contributing to T2D are age, metabolic syndrome, systematic low-grade inflammation, insulin resistance, islet cell autoantibodies, beta-cell dysfunction, and C-peptides (2, 17–23). Some of the factors that promote inflammation are altered functions of specific T lymphocyte cells, B lymphocytes, Th1, and Th17 (24–28). As mentioned above, T2D is also linked with obesity (29), and these diseases are also the case of metabolic disorder; there thus exists strong evidence for an immune-metabolic connection (10, 30, 31).
For this purpose, we have selected publicly available datasets of gene expression for obesity and T2D and have unraveled the genes and the pathways associated with the immune system; we have also focused on the T-cell signaling pathway and its components. We have applied a simplified computational approach to understanding differential gene expression and patterns and the enriched pathways for obesity and T2D. Furthermore, we have also analyzed the genes by using network-level understanding. In the analysis, we observe that there are fewer genes that are often differentially expressed while a comparatively higher number of pathways are shared between them. There are only 4 pathways that are associated with the immune system in the case of obesity and 10 immune-associated pathways in the case of T2D; among these, only two pathways are commonly altered. Here, we have studied the differentially expressed genes (DEGs), enriched pathways, and the connection between immune signaling pathways via their components between obesity and T2D and with the major focus being on the roles of TCR signaling components.
As mentioned above, we have compared immune signaling genes by using gene expression profiling of T2D and obesity for which the datasets have been collected from publicly available data. Here, the results have been presented from work from generalized analyses (i.e., DEGs and enriched pathways) to ensure precise analysis, i.e., TCR signaling.
Gene Expression Profiling for Type-2 Diabetes and Obesity
For diabetes and obesity, the data have been collected from the public database, and we have performed a comparative analysis in terms of altered gene expression patterns and their respective functions. Here, we have presented the DEGs for both the cases T2D and obesity (Figure 1A), and we observe that there are only 12 genes that are differentially expressed in both the case; 456 are T2D-specific genes, and 580 genes are obesity-specific genes. In terms of the altered functions, they share 14 pathways, and the T2D- and obesity-specific pathways have 17 and 18 (Figure 1B, Table 1), respectively. Based on this, it appears that there are only a few genes that are differentially expressed in both cases, while in terms of biological functions, a large number of pathways are affected, which means that T2D and obesity share more biological functions in terms of alterations and their source of alterations, i.e., DEGs are not shared at a large scale. The green nodes are the DEGs that either belong to T2D or obesity cases. For this reason, we have mapped out the network for common DEGs and highlighted common DEGs in red (Figure 1C). Among these common DEGs, FN1, SFRP1, and SNX5 appear to be the major sources of alteration. FN1 appears to be associated with a large number of DEGs, leading to the conclusion that either FN1 is overexpressed by those genes or it promotes the overexpression of the directly associated genes. To understand the alteration in biological function, a heatmap is shown in Figure 1D for those pathways that are commonly enriched. For this plot, the p-values in both conditions have been used. Cell adhesion molecules, PI3K—AKT, and MAPK pathways appear highly enriched in both conditions; there is regulation of the actin cytoskeleton, cAMP, and phospholipase D; Focal adhesion is exclusively and highly enriched in the case of T2D; and hematopoietic cell lineage and phagosomes are potentially enriched in cases of obesity. There are additional pathways that are enriched in both cases but have comparatively higher p-values, and these pathways are the Rap1, tight junction, leukocyte transendothelial migration, ECM-receptor interaction, and Ras signaling pathways. The complete details of the fold changes of the genes are presented in Supplementary Data S1 and S2.
Figure 1 Gene expression profiling. (A) Venn diagram for differentially expressed genes. (B) Venn diagram for enriched pathways. (C) Network representing the DEGs common between T2D and obesity. (D) Commonly enriched pathways: the lighter color or lower p-value refer to higher enrichment, and the darker blue color signifies a lower p-value (close to 0.05) and thus lesser enrichment, though all these pathways are significant in terms of p-values. The p-values are mentioned in the respective color columns and rows.
Table 1 Pathways common between T2D and obesity and specific to T2D and obesity (Venn diagram outcome).
Genes Associated With the Immune System
After analyzing the overall DEGs and all the enriched pathways, we have performed an immune-system-specific analysis for the genes and the pathways. Here, we observed that there are 10 pathways associated with the immune system in the case of T2D and four pathways in the case of obesity (Figure 2A). Hematopoietic cell lineage and leukocyte transendothelial migration pathways are commonly altered in both the cases T2D and obesity, and these two pathways are known to play roles in blood cell development and immune surveillance and inflammation leading to the conclusion that in both T2D and obesity, the blood cell development, immune surveillance, and the inflammation process are affected. Obesity-specific immune signaling pathways are oxytocin and antigen processing and presentation (APP). The oxytocin signaling pathway is known to exert a wide range of effects (central or peripheral), and this pathway is itself known to crosstalk with a number of critical pathways such as RhoA, GPCR, MAPK, ANP-cGMP, and NO-cGMP. APP is one of the most important parts of the entire immune system and is the process by which the protein antigen is ingested by an antigen-presenting cell (APC), partially digested into peptide fragments, and then displayed on the surface of the APC associated with an antigen-presenting molecule (MHC class I or MHC class II) for recognition by specific lymphocytes (T cells). There are eight immune-system-specific pathways in the case of T2D: NK-cell-mediated cytotoxicity, insulin signaling, TCR signaling, ubiquitin-mediated proteolysis, platelet activation BCR signaling, neurotrophin, and oxytocin signaling pathways. Almost all these pathways are major parts of the immune system. Based on this, we could conclude that in the case of T2D, there are more alterations in the gene expression pattern of those genes that are known to control major parts of the immune system, while in the case of obesity there less of an effect than T2D, but these few pathways are also known to be associated with controlling critical roles within the immune system. Furthermore, immune signaling pathways and their components have been presented in the form of a network, where hematopoietic cell lineage appears to be directly controlled by 12 highly significant DEGs; leukocyte transendothelial migration is directly controlled by 10 DEGs; ubiquitin-mediated proteolysis is controlled by 5 DEGs; the insulin signaling pathway is controlled by 6 DEGs; 3 DEGs control platelet activation; 9 DEGs control NK-cell-mediated cytotoxicity; 5 DEGs control TCR and BCR signaling; and 5 DEGs control neurotrophic signaling (Figure 2B).
Figure 2 Genes associated with immune systems. (A) Green-colored genes belong to the immune system when the patient is suffering from T2D and cyan-colored genes belong to the immune system in case of obesity and (B) network of immune signaling pathways and their components (DEGs).
DEGs of T2D and TCR Signaling Network
As mentioned above, there are a number of pathways that are affected as a result of T2D and obesity; in comparison to T2D, there are only four pathways that are, due to obesity, mainly affected, and the TCR signaling pathway is not affected here. The TCR signaling pathway is affected in the case of T2D. Therefore, we have now mapped out the individual genes that are associated with the TCR signaling pathway and presented the networks for all the five genes (SPNS1, PTPN6, CD247, FOS, and PIK3R5) that are components of the TCR signaling pathway and their association with the other TCR genes (Figures 3A–E). These five genes have comparatively higher fold changes, i.e., the expressions and fold changes for SPNS1, PTPN6, CD247, FOS, and PIK3R5 are -5.68, -1.65, 2.72, 1.58, and 2.61, respectively, in the case of T2D. In addition, we have also presented the TCR network (Figure 3F) for T2D and shown the overall connectivities for all the interactors, including these five genes where PTPRC has the highest connectivity, connecting with these five genes and the majority of the interactors. PTPN6 is the second most connected, and there are more genes with much higher connectivity than SPNS1, CDC247, FOS, and PIK3R5: VAV1, LCK, GRB2, FYN, CBL, ZAP70, PAK2, KRAS, and more. The majority of these top-ranked genes in terms of connectivity are considered the core components of TCR activation and signal communication, leading to the conclusion that these core genes (SPNS1, PTPN6, CD247, FOS, and PIK3R5), which are differentially expressed, may lead to the alteration of the TCR core components and may finally potentially alter the TCR signaling process.
Figure 3 DEGs of T2D and TCR signaling network. (A) SPNS1 network, (B) PTPN6 network, (C) CD247 network, (D) FOS network, (E) PIK3R5 network, and (F) the overall association between these T2D DEGs and the TCR genes and the intra-associations.
The SPNS1 gene is the direct component of TCR, NK-cell-mediated cytotoxicity, and FC-epsilon RI signaling pathways; PTPN6 belongs to TCR, BCR, NK-cell-mediated cytotoxicity, JAK-STAT, and the adherens junction; CD247 belongs to TCR and NK-cell-mediated cytotoxicity; and FOS belongs to TCR, BCR, TLR, MAPK, TNF, cAMP, estrogen, oxytocin, prolactin, circadian entrainment, and osteoclast differentiation. PIK3R5 alone belongs to a large number of critical pathways, such as TCR, BCR, NK cell-mediated cytotoxicity, JAK-STAT, ErbB, chemokine, phosphatidylinositol, mTOR, apoptosis, VEGF, focal adhesion, the regulation actin cytoskeleton, etc., and the complete details of the list of pathways for all these five genes are presented in Table 2. These genes presented in the networks (Figure 3) are known to play critical roles in a number of complex human diseases, including T2D and obesity, and the particulars have been presented in more detail in the Discussion section.
Diabetes mellitus (DM) and obesity are considered metabolic-disorder-related diseases, and DM is characterized by hyperglycemia resulting from defects in insulin secretion, insulin action, or both (1, 18, 32). Diabetic chronic hyperglycemia is associated with long-term damage, dysfunction, and failure of various organs such as the kidneys, eyes, nerves, blood vessels, and heart (1, 17, 19, 24, 32–34). T2D is most often found to be frequent in the case of obesity, and it is characterized by abnormal insulin secretion and/or a decreased sensitivity to insulin, which is also known as insulin resistance, and it results in increased blood glucose levels; T2D is known to be strongly associated with obesity (29). There are a number of biological processes and mechanisms that are considered to be the source of association between these conditions via mediating inflammation in adipose tissue and systemic insulin resistance, such as oxidative stress, endoplasmic reticulum stress, hypoxia, amyloid and lipid deposition, lipotoxicity, and glucotoxicity.
In this study, the obese samples totaled 36 samples where 18 were lean needle biopsy and 18 obese biopsy samples; these samples were age-matched to lean (n = 19) and obese (n = 18) female subjects, and the abdominal subcutaneous fat specimens have been collected. The average age (years) of the patients was 35.7 +/- 8.9 (obese-needle), 37.2 +/- 9.3 (obese-surgery), 37.2 +/- 9.9 (lean-needle), and 36.8 +/- 10.1 (lean-surgery), and the average body mass index (BMI) was 44.0 +/- 6.9 (obese-needle), 46.1 +/- 4.4 (obese-surgery), 21.9 +/- 1.9 (lean-needle), and 21.9 +/- 2.4 (lean-surgery) (35). While in the case of T2D, GSE121 (36) was used for gene expression profiling, which contains insulin resistance (IR) and insulin sensitivity (IS) for type II diabetes, while 11 were control samples, and the age of patients was not mentioned. The main goals of the authors (in the case of T2D) were to identify DEGs in skeletal muscle in insulin resistance, a major risk factor for T2D (non-insulin-dependent) diabetes, and to compare the gene expression patterns of skeletal muscle tissues from 18 insulin-sensitive versus 17 insulin-resistant, equally obese, and non-diabetic Pima Indians for which global transcript profiling has been performed (36). Different from the isolated study of obesity or T2D, our goal differed much and accomplished the following: presented a comparative study of both these two diseases; greatly explored the role of immune systems main TCR signaling pathway; and uniquely presented the network of immune signaling pathways and their components (Figure 2B).
With such levels of complexities, there exist a number of challenges for understanding the mechanism; we have therefore thoroughly analyzed the DEGs and have enriched pathways with the main focus on immune signaling pathways and the networks of DEGs belonging to the immune system. Among the top enriched pathways are the following: cell adhesion molecules, PI3K—AKT, MAPK, regulation of actin cytoskeleton, cAMP, phospholipase D, Focal adhesion, hematopoietic cell lineage, phagosome, Rap1, tight junction, leukocyte transendothelial migration, ECM-receptor interaction, and Ras signaling pathways.
For immune-system-related analysis, we observed that there are 10 pathways associated with the immune system in the case of T2D and 4 pathways in the case of obesity. Hematopoietic cell lineage and leukocyte transendothelial migration pathways are commonly altered in both the cases T2D and obesity and are known to play roles in blood cell development and immune surveillance and inflammation. Oxytocin and APP are obesity-specific immune signaling pathways and NK-cell-mediated cytotoxicity, insulin signaling, TCR signaling, ubiquitin-mediated proteolysis, platelet activation BCR signaling, neurotrophin, and oxytocin signaling pathways are T2D immune-system-specific pathways. Almost all these pathways (obese and T2D immune signaling pathways) are major parts of the immune system and are known to control major parts of the immune system. Furthermore, it is also known that altered TCR signaling is associated with T2D. Therefore, we have now mapped out the individual genes that are associated with the TCR signaling pathway and presented the networks for all the five genes (SPNS1, PTPN6, CD247, FOS, and PIK3R5), which are the components of the TCR signaling pathway, and their association with the other TCR genes (Figure 3). These five genes, SPNS1, PTPN6, CD247, FOS, and PIK3R5, are the potential components of TCR and also a number of additional immune-system-associated pathways. Among those pathways are NK-cell-mediated cytotoxicity, FC-epsilon RI signaling pathways, BCR, JAK-STAT, adherens junction, TLR, MAPK, TNF, cAMP, estrogen, oxytocin, prolactin, circadian entrainment, and osteoclast differentiation. PIK3R5 alone belongs to a large number of critical pathways such as TCR, BCR, NK cell-mediated cytotoxicity, JAK-STAT, ErbB, chemokine, phosphatidylinositol, mTOR, apoptosis, VEGF, focal adhesion, the regulation actin cytoskeleton, etc. (37). The SPNS1 gene is known to be associated with a number of human diseases, including T2D and obesity (38–42), which is mainly associated with lysosome storage disorder (LSD). PTPN6 is also known to be directly associated with controlling the inflammatory system, and altered PTPN6 (expression/mutation) may lead to potential human diseases (43–45). CD247 is known to regulate blood pressure, which is a critically affected factor in the case of T2D and obesity, by altering T-cell infiltration in the kidney; it also has an important impact on other human diseases (46–48). FOS is considered to play a role in inflammatory bone and skin disease (49, 50) and also in multiple other human diseases. PIK3R5 contributes to cell survival and is a critical pathway component to control many diseases (51–53).
As we can see in the results, PIK3R5 and FOS, which are the components of the TCR signaling pathway, control a large set of biological functions, most of which are a part of the immune system, and CD247, PTPN6, and SPNS1 comparatively control fewer pathways; in terms of connectivity, however, PTPN6 has the highest connectivity followed by SPNS1, CD247, FOS, and PIK3R5. It means that PIK3R5 simultaneously controls a large number of biological functions but does not impact more genes; in other words, the higher connectivity of the gene will increase the likelihood of affecting the expression patterns of the gene or vice versa. In terms of future perspectives, this study will be helpful for the development of therapeutic approaches and drug targeting.
Materials and Methods
For obese people, we used the GSE12050 (35) dataset; 18 of 36 samples were lean needle biopsies, and these were compared with 18 obese biopsy samples. The samples were age-matched lean (n = 19) and obese (n = 18) female subjects. For sample details, it has been mentioned that abdominal subcutaneous fat specimens have been collected by using a needle biopsy approach from 9 obese patient samples and 10 lean samples. The average ages of the patients were 35.7 (obese-needle), 37.2 (obese-surgery), 37.2 (lean-needle), and 36.8 (lean-surgery). A total of nine obese patient samples were collected from abdominal subcutaneous fat biopsies and nine lean samples from abdominal fat biopsies (35). For T2D, GSE121 (GPL101 platform) (36) has been used for T2D gene expression profiling, which contains 5 IR and 5 IS samples with type II diabetes while 11 are control samples, and the age of patients have not been mentioned. The samples were non-diabetic Pima Indians or Tohono O’Odham Indians who were classified as IS or IR (36). For this dataset, the details and the goal submitted by the submitters were as follows: to identify DEGs in skeletal muscle in insulin resistance a major risk factor for T2D (non-insulin-dependent) diabetes and to compare gene the expression pattern of skeletal muscle tissues from 18 insulin-sensitive versus 17 insulin-resistant, equally obese, non-diabetic Pima Indians for which global transcript profiling has been performed (36).
These datasets were collected from GEO and have been processed for expression value calculation, fold change, and p-value calculation by the in-built code at GEO, i.e., GEO2R (54, 55). The list of DEGs has been processed for pathway enrichment analysis and understanding (56, 57). The KEGG (37) database has been used for calculating enriched pathways, and our own code has been used for pathway and network analysis (58–61). To retrieve the DEGs association as a network, the FunCoup2.0 network database (62) was utilized throughout the work, and the network was visualized with the help of cytoscapes (63). For most of our coding and calculations, MATLAB was used.
Now, we briefly summarize all the steps used are raw file processing, intensity calculation, and normalization. As the high-throughput data of the mainly microarray-based gene expression values are calculated as intensity ranging from the lowest to highest value, it is a necessity to have optimized or scaled-down values so we normally convert them to log2 values. For normalization (64–66), GCRMA (67–71), RMA, and EB are the most commonly used approaches. Here, we have used EB for raw intensity normalization. After normalization, we proceed towards our goal, which is to understand the gene expression patterns (56, 72) and their inferred functions (56, 57). For differential gene expression prediction and statistical analysis, MATLAB functions (e.g., mattest) have been used and volcano plot. From mattest, we get the p-values, and we then go for p-value corrections, assign the p-value threshold, and use the volcano plot for fold change calculation and the selection of genes with p-values less than 0.05 and fold changes within the range >= +2.0 and <= -2.0. This range of fold change along with this p-value cutoff are commonly used, while, in our study, we selected the top-ranked genes to have the strict condition of selection based on fold change so the fold changes are either much higher or lower than the cutoff fold change. For pathway analysis, we used the KEGG (37) database and have our own code designed for the pathway and network analysis (73–75). We have described the fold change, and, by considering genes as DEGs, we have applied thresholds for FC >= +2.0 (upregulated genes) and FC <= -2.0 (downregulated genes); these DEGs have been considered for pathway enrichment analysis. But for heatmap and network presentations, the top 100 DEGs (50 up and 50 down) genes have been presented (76, 77).
Here, we observe that there are fewer genes that are commonly differentially expressed, while there are comparatively more pathways that are shared between them. There are only 4 pathways that are associated with the immune system in the case of obesity and 10 immune-associated pathways in the case of T2D, and, among them, only 2 pathways are commonly altered. Here, we have studied the differentially expressed genes (DEGs), enriched pathways, and the connection between immune signaling pathways via their components between obesity and T2D with a focus on the roles of the TCR signaling components. We conclude that in the case of T2D, there are more alterations in the gene expression pattern of those genes that are known to control the major parts of the immune system than in the case of obesity, where there is less of an effect than with T2D, but these few pathways are also known to be associated with controlling critical roles of the immune system.
Data Availability Statement
The original contributions presented in the study are included in the article/Supplementary Material. Further inquiries can be directed to the corresponding authors.
Conceptualization, BE, HS, HC, and MM. Methodology, BE, HS, HC, and MM. Software, MM. Validation, BE, HS, HC, and ZZ. Formal analysis, BE, HS, HC, and MM. Investigation: BE and MM. Resources, BE and MM. Data curation, BE and MM. Writing—original draft preparation, BE, HS, HC, and MM. Writing—review and editing, BE, HS, HC, and MM. Visualization, BE, HS, HC, and MM. Supervision, HC and MM. Project administration, BE and HC. Funding acquisition, BE. All authors contributed to the article and approved the submitted version.
This research was funded by the Deanship of Scientific Research (DSR) at King Abdulaziz University, Jeddah, grant number RG-5-140-39 and Swedish Cancer Society.
Conflict of Interest
The authors declare that the research was conducted in the absence of anycommercial or financial relationships that could be construed as a potentialconflict of interest.
The Deanship of Scientific Research (DSR) at King Abdulaziz University, Jeddah, and Swedish Cancer Society funded this project, under grant number RG-5-140-39. The work has been supported by funding by the Swedish cancer society. We are also thankful to Karolinska Institute, Solna, Sweden and King Abdulaziz University for offering us the environment to carry out the entire research work and all other supports.
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fendo.2021.641888/full#supplementary-material
DEGs, Differentially expressed genes; TCR, T cell receptor; ISPs, Immune signaling pathways; T2D, Type-2 diabetes; APP, Antigen processing and presentation; APC, Antigen-presenting cell.
10. Richard C, Wadowski M, Goruk S, Cameron L, Sharma AM, Field CJ. Individuals with obesity and type 2 diabetes have additional immune dysfunction compared with obese individuals who are metabolically healthy. BMJ Open Diabetes Res Care (2017) 5:e000379. doi: 10.1111/j.1365-2249.2010.04292.x
15. Draznin B. Molecular Mechanisms of Insulin Resistance: Serine Phosphorylation of Insulin Receptor Substrate-1 and Increased Expression of p85 : The Two Sides of a Coin. Diabetes (2006) 55:2392–7. doi: 10.2337/db06-0391
18. Rivera JF, Costes S, Gurlo T, Glabe CG, Butler PC. Autophagy defends pancreatic β cells from human islet amyloid polypeptide-induced toxicity. J Clin Invest (2014) 124:3489–500. doi: 10.1172/JCI71981DS1
24. de Candia P, Prattichizzo F, Garavelli S, De Rosa V, Galgani M, Di Rella F, et al. Type 2 Diabetes: How Much of an Autoimmune Disease? Front Endocrinol (2019) 10:348. doi: 10.1016/S2213-8587(18)30051-2
26. Sakaguchi N, Takahashi T, Hata H, Nomura T, Tagami T, Yamazaki S, et al. Altered thymic T-cell selection due to a mutation of the ZAP-70 gene causes autoimmune arthritis in mice. Nat Cell Biol (2003) 426:454–60. doi: 10.1038/nature02119
35. Mutch DM, Tordjman J, Pelloux V, Hanczar B, Henegar C, Poitou C, et al. Needle and surgical biopsy techniques differentially affect adipose tissue gene expression profiles. Am J Clin Nutr (2008) 89:51–7. doi: 10.3945/ajcn.2008.26802
36. Yang X, Pratley R, Tokraks S, Bogardus C, Permana P. Microarray profiling of skeletal muscle tissues from equally obese, non-diabetic insulin-sensitive and insulin-resistant Pima Indians. Diabetologia (2002) 45:1584–93. doi: 10.1007/s00125-002-0905-7
38. Sasaki T, Lian S, Qi J, Bayliss PE, Carr CE, Johnson JL, et al. Aberrant Autolysosomal Regulation Is Linked to The Induction of Embryonic Senescence: Differential Roles of Beclin 1 and p53 in Vertebrate Spns1 Deficiency. PloS Genet (2014) 10:e1004409. doi: 10.1371/journal.pgen.1004409.s023
39. Sasaki T, Lian S, Khan A, Llop JR, Samuelson AV, Chen W, et al. Autolysosome biogenesis and developmental senescence are regulated by both Spns1 and v-ATPase. Autophagy (2017) 13:386–403. doi: 10.1080/15548627.2016.1256934
40. Yanagisawa H, Ishii T, Endo K, Kawakami E, Nagao K, Miyashita T, et al. L-leucine and SPNS1 coordinately ameliorate dysfunction of autophagy in mouse and human Niemann-Pick type C disease. Sci Rep (2017) 7:1–9. doi: 10.1038/s41598-017-15305-9
42. Rong Y, McPhee CK, Deng S, Huang L, Chen L, Liu M, et al. Spinster is required for autophagic lysosome reformation and mTOR reactivation following starvation. Proc Natl Acad Sci (2011) 108:7826–31. doi: 10.1038/nature09406
43. Chen Z, Shojaee S, Buchner M, Geng H, Lee JW, Klemm L, et al. Signalling thresholds and negative B-cell selection in acute lymphoblastic leukaemia. Nature (2015) 521:1–16. doi: 10.1038/nature14231
44. Demosthenous C, Han JJ, Hu G, Stenson M, Gupta M. Loss of function mutations in PTPN6 promote STAT3 deregulation via JAK3 kinase in diffuse large B-cell lymphoma. Oncotarget (2015) 6:44703–13. doi: 10.18632/oncotarget.6300
45. Kanwal Z, Zakrzewska A, Hertog den J, Spaink HP, Schaaf MJM, Meijer AH. Deficiency in Hematopoietic Phosphatase Ptpn6/Shp1 Hyperactivates the Innate Immune System and Impairs Control of Bacterial Infections in Zebrafish Embryos. J Immunol (2013) 190:1631–45. doi: 10.1016/j.cell.2011.12.023
46. Rudemiller N, Lund H, Jacob HJ, Geurts AM, Mattson DL. CD247 Modulates Blood Pressure by Altering T-Lymphocyte Infiltration in the Kidney. Hypertension (2014) 63:559–64. doi: 10.1002/j.1460-2075.1993.tb06120.x
47. Ye W, Zhou Y, Xu B, Zhu D, Rui X, Xu M, et al. CD247 expression is associated with differentiation and classification in ovarian cancer. Medicine (2019) 98:e18407. doi: 10.1097/MD.0000000000018407
48. Martins M, Williams AH, Comeau M, Marion M, Ziegler JT, Freedman BI, et al. Genetic association of CD247 (CD3ζ) with SLE in a large-scale multiethnic study. Genes Immun (2019) 16:1–9. doi: 10.1038/gene.2014.73
49. Zenz R, Eferl R, Scheinecker C, Redlich K, Smolen J, Schonthaler HB, et al. Activator protein 1 (Fos/Jun) functions in inflammatory bone and skin disease. Arthritis Res Ther (2007) 10:201. doi: 10.1186/ar2338
50. Davani-Davari D, Negahdaripour M, Karimzadeh I, Seifan M, Mohkam M, Masoumi S, et al. Prebiotics: Definition, Types, Sources, Mechanisms, and Clinical Applications. Foods (2019) 8:92. doi: 10.1186/s12929-018-0493-6
52. Dornan GL, Burke JE. Disease Mediated by Oncogenic and Primary immunodeficiency Mutations in Class iA Phosphoinositide 3-Kinases. Front Immunol (2018) 9:575.s001:575.s001. doi: 10.3389/fimmu.2018.00575.s001
53. Tassan Al N, Khalil D, Shinwari J, Sharif Al L, Bavi P, Abduljaleel Z, et al. A missense mutation in PIK3R5gene in a family with ataxia and oculomotor apraxia. Hum Mutat (2011) 33:351–4. doi: 10.1186/1471-2105-9-40
56. Subramanian A, Tamayo P, Mootha VK, Mukherjee S, Ebert BL, Gillette MA, et al. Gene set enrichment analysis: A knowledge-based approach for interpreting genome-wide expression profiles. Proc Natl Acad Sci USA (2005) 102:15545. doi: 10.1073/pnas.0506580102
57. Mi H, Poudel S, Muruganujan A, Casagrande JT, Thomas PD. PANTHER version 10: expanded protein families and functions, and analysis tools. Nucleic Acids Res (2016) 44:D336–42. doi: 10.1093/nar/gkv1194
58. Warsi MK, Kamal MA, Baeshen MN, Izhari MA, Mobashir AFAM. Comparative Study of Gene Expression Profiling Unravels Functions associated with Pathogenesis of Dengue Infection. Curr Pharm Design (2020) 26(41):5293–9. doi: 10.2174/1381612826666201106093148
61. Mobashir M, Madhusudhan T, Isermann B, Beyer T, Schraven B. Negative Interactions and Feedback Regulations Are Required for Transient Cellular Response. Sci Rep (2014) 4:3718. doi: 10.1038/srep03718
66. Ideker T, Thorsson V, Siegel AF, Hood LE. Testing for differentially-expressed genes by maximum-likelihood analysis of microarray data. J Comput Biol (2000) 7:805–17. doi: 10.1089/10665270050514945
68. Chen K-H, Wang K-J, Tsai M-L, Wang K-M, Adrian AM, Cheng W-C, et al. Gene selection for cancer identification: a decision tree model empowered by particle swarm optimization algorithm. BMC Bioinf (2014) 15:1–10. doi: 10.1186/1471-2105-15-49
69. Bild AH, Parker JS, Gustafson AM, Acharya CR, Hoadley KA, Anders C, et al. An integration of complementary strategies for gene-expression analysis to reveal novel therapeutic opportunities for breast cancer. Breast Cancer Res (2009) 11:R55. doi: 10.1186/bcr2344
72. Lapointe J, Li C, Higgins JP, van de Rijn M, Bair E, Montgomery K, et al. Gene expression profiling identifies clinically relevant subtypes of prostate cancer. Proc Natl Acad Sci USA (2004) 101:811–6. doi: 10.1073/pnas.0304146101
74. Asgharzadeh MR, Pourseif MM, Barar J, Eskandani M, Niya MJ, Mashayekhi MR, et al. Functional expression and impact of testis-specific gene antigen 10 in breast cancer: a combined in vitro and in silico analysis. Tabriz Univ Med Sci (2019) 9:145–59. doi: 10.15171/bi.2019.19
75. Saadat YR, Pourseif MM, Vahed SZ, Barzegari A, Omidi Y, Barar J. Modulatory Role of Vaginal-Isolated Lactococcus lactison the Expression of miR-21, miR-200b, and TLR-4 in CAOV-4 Cells and In Silico Revalidation. Probiotics Antimicrob Proteins (2020) 12(3):1083–96. doi: 10.1007/s12602-019-09596-9
76. Kamal MA, Warsi MK, Alnajeebi A, Ali HA, Helmi N, Izhari MA, et al. Gene expression profiling and clinical relevance unravel the role hypoxia and immune signaling genes and pathways in breast cancer: Role of hypoxia and immune signaling genes in breast cancer. jimsa (2020) 1:1–9. doi: 10.36013/jimsa 2.v1i1.3
Keywords: type-2 diabetes and obesity, insulin action and resistance, immune system, enriched pathways, differentially expressed genes, TCR network, network-level understanding
Citation: Eldakhakhny BM, Al Sadoun H, Choudhry H and Mobashir M (2021) In-Silico Study of Immune System Associated Genes in Case of Type-2 Diabetes With Insulin Action and Resistance, and/or Obesity. Front. Endocrinol. 12:641888. doi: 10.3389/fendo.2021.641888
Received: 15 December 2020; Accepted: 08 March 2021;
Published: 13 April 2021.
Edited by:Undurti Narasimha Das, UND Life Sciences LLC, United States
Reviewed by:Xiuli Wu, Jinan University, China
Marcia Hiriart, Universidad Nacional Autonoma de Mexico, Mexico
Copyright © 2021 Eldakhakhny, Al Sadoun, Choudhry and Mobashir. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.