Deciphering Signaling Pathway Networks to Understand the Molecular Mechanisms of Metformin Action

A drug exerts its effects typically through a signal transduction cascade, which is non-linear and involves intertwined networks of multiple signaling pathways. Construction of such a signaling pathway network (SPNetwork) can enable identification of novel drug targets and deep understanding of drug action. However, it is challenging to synopsize critical components of these interwoven pathways into one network. To tackle this issue, we developed a novel computational framework, the Drug-specific Signaling Pathway Network (DSPathNet). The DSPathNet amalgamates the prior drug knowledge and drug-induced gene expression via random walk algorithms. Using the drug metformin, we illustrated this framework and obtained one metformin-specific SPNetwork containing 477 nodes and 1,366 edges. To evaluate this network, we performed the gene set enrichment analysis using the disease genes of type 2 diabetes (T2D) and cancer, one T2D genome-wide association study (GWAS) dataset, three cancer GWAS datasets, and one GWAS dataset of cancer patients with T2D on metformin. The results showed that the metformin network was significantly enriched with disease genes for both T2D and cancer, and that the network also included genes that may be associated with metformin-associated cancer survival. Furthermore, from the metformin SPNetwork and common genes to T2D and cancer, we generated a subnetwork to highlight the molecule crosstalk between T2D and cancer. The follow-up network analyses and literature mining revealed that seven genes (CDKN1A, ESR1, MAX, MYC, PPARGC1A, SP1, and STK11) and one novel MYC-centered pathway with CDKN1A, SP1, and STK11 might play important roles in metformin’s antidiabetic and anticancer effects. Some results are supported by previous studies. In summary, our study 1) develops a novel framework to construct drug-specific signal transduction networks; 2) provides insights into the molecular mode of metformin; 3) serves a model for exploring signaling pathways to facilitate understanding of drug action, disease pathogenesis, and identification of drug targets.     

Association of Single Nucleotide Polymorphisms in Wnt Signaling Pathway Genes with Breast Cancer in Saudi Patients

Breast cancer is a complex heterogeneous disease involving genetic and epigenetic alterations in genes encoding proteins that are components of various signaling pathways. Candidate gene approach have identified association of genetic variants in the Wnt signaling pathway genes and increased susceptibility to several diseases including breast cancer. Due to the rarity of somatic mutations in key genes of Wnt pathway, we investigated the association of genetic variants in these genes with predisposition to breast cancers. We performed a case-control study to identify risk variants by examining 15 SNPs located in 8 genes associated with Wnt signaling. Genotypic analysis of individual locus showed statistically significant association of five SNPs located in β-catenin, AXIN2, DKK3, SFRP3 and TCF7L2 with breast cancers. Increased risk was observed only with the SNP in β-catenin while the other four SNPs conferred protection against breast cancers. Majority of these associations persisted after stratification of the cases based on estrogen receptor status and age of on-set of breast cancer. The rs7775 SNP in exon 6 of SFRP3 gene that codes for either arginine or glycine exhibited very strong association with breast cancer, even after Bonferroni''s correction. Apart from these five variants, rs3923086 in AXIN2 and rs3763511 in DKK4 that did not show any association in the overall population were significantly associated with early on-set and estrogen receptor negative breast cancers, respectively. This is the first study to utilize pathway based approach to identify association of risk variants in the Wnt signaling pathway genes with breast cancers. Confirmation of our findings in larger populations of different ethnicities would provide evidence for the role of Wnt pathway as well as screening markers for early detection of breast carcinomas.     

A systems genetics approach provides a bridge from discovered genetic variants to biological pathways in rheumatoid arthritis

Genome-wide association studies (GWAS) have yielded novel genetic loci underlying common diseases. We propose a systems genetics approach to utilize these discoveries for better understanding of the genetic architecture of rheumatoid arthritis (RA). Current evidence of genetic associations with RA was sought through PubMed and the NHGRI GWAS catalog. The associations of 15 single nucleotide polymorphisms and HLA-DRB1 alleles were confirmed in 1,287 cases and 1,500 controls of Japanese subjects. Among these, HLA-DRB1 alleles and eight SNPs showed significant associations and all but one of the variants had the same direction of effect as identified in the previous studies, indicating that the genetic risk factors underlying RA are shared across populations. By receiver operating characteristic curve analysis, the area under the curve (AUC) for the genetic risk score based on the selected variants was 68.4%. For seropositive RA patients only, the AUC improved to 70.9%, indicating good but suboptimal predictive ability. A simulation study shows that more than 200 additional loci with similar effect size as recent GWAS findings or 20 rare variants with intermediate effects are needed to achieve AUC = 80.0%. We performed the random walk with restart (RWR) algorithm to prioritize genes for future mapping studies. The performance of the algorithm was confirmed by leave-one-out cross-validation. The RWR algorithm pointed to ZAP70 in the first rank, in which mutation causes RA-like autoimmune arthritis in mice. By applying the hierarchical clustering method to a subnetwork comprising RA-associated genes and top-ranked genes by the RWR, we found three functional modules relevant to RA etiology: “leukocyte activation and differentiation”, “pattern-recognition receptor signaling pathway”, and “chemokines and their receptors”.These results suggest that the systems genetics approach is useful to find directions of future mapping strategies to illuminate biological pathways.     

Pooled Sequencing of 531 Genes in Inflammatory Bowel Disease Identifies an Associated Rare Variant in BTNL2 and Implicates Other Immune Related Genes

The contribution of rare coding sequence variants to genetic susceptibility in complex disorders is an important but unresolved question. Most studies thus far have investigated a limited number of genes from regions which contain common disease associated variants. Here we investigate this in inflammatory bowel disease by sequencing the exons and proximal promoters of 531 genes selected from both genome-wide association studies and pathway analysis in pooled DNA panels from 474 cases of Crohn’s disease and 480 controls. 80 variants with evidence of association in the sequencing experiment or with potential functional significance were selected for follow up genotyping in 6,507 IBD cases and 3,064 population controls. The top 5 disease associated variants were genotyped in an extension panel of 3,662 IBD cases and 3,639 controls, and tested for association in a combined analysis of 10,147 IBD cases and 7,008 controls. A rare coding variant p.G454C in the BTNL2 gene within the major histocompatibility complex was significantly associated with increased risk for IBD (p = 9.65x10⁻¹⁰, OR = 2.3[95% CI = 1.75–3.04]), but was independent of the known common associated CD and UC variants at this locus. Rare (<1%) and low frequency (1–5%) variants in 3 additional genes showed suggestive association (p<0.005) with either an increased risk (ARIH2 c.338-6C>T) or decreased risk (IL12B p.V298F, and NICN p.H191R) of IBD. These results provide additional insights into the involvement of the inhibition of T cell activation in the development of both sub-phenotypes of inflammatory bowel disease. We suggest that although rare coding variants may make a modest overall contribution to complex disease susceptibility, they can inform our understanding of the molecular pathways that contribute to pathogenesis.     

Pathway Network Analyses for Autism Reveal Multisystem Involvement,Major Overlaps with Other Diseases and Convergence upon MAPK and Calcium Signaling

We used established databases in standard ways to systematically characterize gene ontologies, pathways and functional linkages in the large set of genes now associated with autism spectrum disorders (ASDs). These conditions are particularly challenging—they lack clear pathognomonic biological markers, they involve great heterogeneity across multiple levels (genes, systemic biological and brain characteristics, and nuances of behavioral manifestations)—and yet everyone with this diagnosis meets the same defining behavioral criteria. Using the human gene list from Simons Foundation Autism Research Initiative (SFARI) we performed gene set enrichment analysis with the Kyoto Encyclopedia of Genes and Genomes (KEGG) Pathway Database, and then derived a pathway network from pathway-pathway functional interactions again in reference to KEGG. Through identifying the GO (Gene Ontology) groups in which SFARI genes were enriched, mapping the coherence between pathways and GO groups, and ranking the relative strengths of representation of pathway network components, we 1) identified 10 disease-associated and 30 function-associated pathways 2) revealed calcium signaling pathway and neuroactive ligand-receptor interaction as the most enriched, statistically significant pathways from the enrichment analysis, 3) showed calcium signaling pathways and MAPK signaling pathway to be interactive hubs with other pathways and also to be involved with pervasively present biological processes, 4) found convergent indications that the process “calcium-PRC (protein kinase C)-Ras-Raf-MAPK/ERK” is likely a major contributor to ASD pathophysiology, and 5) noted that perturbations associated with KEGG’s category of environmental information processing were common. These findings support the idea that ASD-associated genes may contribute not only to core features of ASD themselves but also to vulnerability to other chronic and systemic problems potentially including cancer, metabolic conditions and heart diseases. ASDs may thus arise, or emerge, from underlying vulnerabilities related to pleiotropic genes associated with pervasively important molecular mechanisms, vulnerability to environmental input and multiple systemic co-morbidities.     

Whole Exome Sequencing Identifies Novel Recurrently Mutated Genes in Patients with Splenic Marginal Zone Lymphoma

The pathogenesis of splenic marginal zone lymphoma (SMZL) remains largely unknown. Recent high-throughput sequencing studies have identified recurrent mutations in key pathways, most notably NOTCH2 mutations in >25% of patients. These studies are based on small, heterogeneous discovery cohorts, and therefore only captured a fraction of the lesions present in the SMZL genome. To identify further novel pathogenic mutations within related biochemical pathways, we applied whole exome sequencing (WES) and copy number (CN) analysis to a biologically and clinically homogeneous cohort of seven SMZL patients with 7q abnormalities and IGHV1-2*04 gene usage. We identified 173 somatic non-silent variants, affecting 160 distinct genes. In additional to providing independent validation of the presence of mutation in several previously reported genes (NOTCH2, TNFAIP3, MAP3K14, MLL2 and SPEN), our study defined eight additional recurrently mutated genes in SMZL; these genes are CREBBP, CBFA2T3, AMOTL1, FAT4, FBXO11, PLA2G4D, TRRAP and USH2A. By integrating our WES and CN data we identified three mutated putative candidate genes targeted by 7q deletions (CUL1, EZH2 and FLNC), with FLNC positioned within the well-characterized 7q minimally deleted region. Taken together, this work expands the reported directory of recurrently mutated cancer genes in this disease, thereby expanding our understanding of SMZL pathogenesis. Ultimately, this work will help to establish a stratified approach to care including the possibility of targeted therapy.     

Combinations of Susceptibility Genes Are Associated with Higher Risk for Multiple Sclerosis and Imply Disease Course Specificity

Multiple sclerosis (MS) is a chronic autoimmune disease of the central nervous system that predominantly affects young adults. The genetic contributions to this multifactorial disease were underscored by a genome wide association study (GWAS) conducted by the International Multiple Sclerosis Genetic Consortium in a multinational cohort prompting the discovery of 57 non-MHC MS-associated common genetic variants. Hitherto, few of these newly reported variants have been replicated in larger independent patient cohorts. We genotyped a cohort of 1033 MS patients and 644 healthy controls with a consistent genetic background for the 57 non-MHC variants reported to be associated with MS by the first large GWAS as well as the HLA DRB1*1501 tagging SNP rs3135388. We robustly replicated three of the 57 non-MHC reported MS-associated single nucleotide polymorphisms (SNPs). In addition, our study revealed several genotype-genotype combinations with an evidently higher degree of disease association than the genotypes of the single SNPs. We further correlated well-defined clinical phenotypes, i.e. ataxia, visual impairment due to optic neuritis and paresis with single SNPs and genotype combinations, and identified several associations. The results may open new avenues for clinical implications of the MS associated genetic variants reported from large GWAS.     

Detecting shared pathways linked to rheumatoid arthritis with other autoimmune diseases in a in silico analysis

Pathway-based analysis approach has exploded in use during the last several years. It is successful in recognizing additional biological insight of disease and finding groupings of risk genes that represent disease developing processes. Therefore, shared pathways, with pleiotropic effects, are important for understanding similar pathogenesis and indicating the common genetic origin of certain diseases. Here, we present a pathway analysis to reveal the potential disease associations between RA and three potential RA-related autoimmune diseases: psoriasis, diabetes mellitus, type 1 (T1D) and systemic lupus erythematosus (SLE). First, a comprehensive knowledge mining of public databases is performed to discover risk genes associated with RA, T1D, SLE and psoriasis; then by enrichment test of these genes, disease-related risk pathways are detected to recognize the pathways common for RA and three other diseases. Finally, the underlying disease associations are evaluated with the association rules mining method. In total, we identify multiple RA risk pathways with significant pleiotropic effects, the most unsurprising of which are the immunology related pathways. Meanwhile for the first time we highlight the involvement of the viral myocarditis pathway related to cardiovascular disease (CVD) in autoimmune diseases such as RA, psoriasis, T1D and SLE. Further Association rule mining results validate the strong association between RA and T1D and RA and SLE. It is clear that pleiotropy is a common property of pathways associated with disease traits. We provide novel pathway associations among RA and three autoimmune diseases. These results ascertain that there are shared genetic risk profiles that predispose individuals to autoimmune diseases.     

Identification of Novel Type 2 Diabetes Candidate Genes Involved in the Crosstalk between the Mitochondrial and the Insulin Signaling Systems

Type 2 Diabetes (T2D) is a highly prevalent chronic metabolic disease with strong co-morbidity with obesity and cardiovascular diseases. There is growing evidence supporting the notion that a crosstalk between mitochondria and the insulin signaling cascade could be involved in the etiology of T2D and insulin resistance. In this study we investigated the molecular basis of this crosstalk by using systems biology approaches. We combined, filtered, and interrogated different types of functional interaction data, such as direct protein–protein interactions, co-expression analyses, and metabolic and signaling dependencies. As a result, we constructed the mitochondria-insulin (MITIN) network, which highlights 286 genes as candidate functional linkers between these two systems. The results of internal gene expression analysis of three independent experimental models of mitochondria and insulin signaling perturbations further support the connecting roles of these genes. In addition, we further assessed whether these genes are involved in the etiology of T2D using the genome-wide association study meta-analysis from the DIAGRAM consortium, involving 8,130 T2D cases and 38,987 controls. We found modest enrichment of genes associated with T2D amongst our linker genes (p = 0.0549), including three already validated T2D SNPs and 15 additional SNPs, which, when combined, were collectively associated to increased fasting glucose levels according to MAGIC genome wide meta-analysis (p = 8.12×10⁻⁵). This study highlights the potential of combining systems biology, experimental, and genome-wide association data mining for identifying novel genes and related variants that increase vulnerability to complex diseases.     

Maternal coding variants in complement receptor 1 and spontaneous idiopathic preterm birth

Preterm birth (PTB) is a major global public health concern. However, little is known about the pathophysiology of spontaneous idiopathic PTB. We tested the hypothesis that rare variants in families would target specific genes and pathways that contribute to PTB risk in the general population. Whole-exome sequencing was performed on 10 PTB mothers from densely affected families including two mother–daughter pairs. We identified novel variants shared between the two mother–daughter pairs when compared to a 1000 Genomes Project background exome file and investigated these genes for pathway aggregation using the Kyoto Encyclopedia of Genes and Genomes (KEGG). Genes in enriched pathways were then surveyed in the other six PTB exomes and tested for association in a larger number of nuclear families. The KEGG complement and coagulation cascade was one of the most enriched pathways in our two mother–daughter pairs. When the six genes found in this pathway (CFH, CR1, F13B, F5, CR2, and C4BPA) were examined for novel missense variants, half of all the exomes harbored at least one. Association analysis of variants in these six gene regions in nuclear families from Finland (237 cases and 328 controls) found statistically significant associations after multiple test corrections in three CR1 SNPs; the strongest in an exonic missense SNP, rs6691117, p value = 6.91e?5, OR = 1.71. Our results demonstrate the importance of the complement and coagulation cascades in the pathophysiology of PTB, and suggest potential screening and intervention approaches to prevent prematurity that target this pathway.     

Gene-Based Association Analysis Identified Novel Genes Associated with Bone Mineral Density

Genetic factors contribute to the variation of bone mineral density (BMD), which is a major risk factor of osteoporosis. The aim of this study was to identify more “novel” genes for BMD. Based on the publicly available SNP-based P values, we performed an initial gene-based analysis in a total of 32,961 individuals. Furthermore, we performed differential expression, pathway and protein-protein interaction analyses to find supplementary evidence to support the significance of the identified genes. About 21,695 genes for femoral neck (FN)-BMD and 21,683 genes for lumbar spine (LS)-BMD were analyzed using gene-based association analysis. A total of 35 FN-BMD associated genes and 53 LS-BMD associated genes were identified (P < 2.3×10^-6) after Bonferroni correction. Among them, 64 genes have not been reported in previous SNP-based genome-wide association studies. Differential expression analysis further supported the significant associations of 14 genes with FN-BMD and 19 genes with LS-BMD. Especially, WNT3 and WNT9B in the Wnt signaling pathway for FN-BMD were further supported by pathway analysis and protein-protein interaction analysis. The present study took the advantage of gene-based association method to perform a supplementary analysis of the GWAS dataset and found some BMD-associated genes. The evidence taken together supported the importance of Wnt signaling pathway genes in determining osteoporosis. Our findings provided more insights into the genetic basis of osteoporosis.     

Common Inherited Variation in Mitochondrial Genes Is Not Enriched for Associations with Type 2 Diabetes or Related Glycemic Traits

Mitochondrial dysfunction has been observed in skeletal muscle of people with diabetes and insulin-resistant individuals. Furthermore, inherited mutations in mitochondrial DNA can cause a rare form of diabetes. However, it is unclear whether mitochondrial dysfunction is a primary cause of the common form of diabetes. To date, common genetic variants robustly associated with type 2 diabetes (T2D) are not known to affect mitochondrial function. One possibility is that multiple mitochondrial genes contain modest genetic effects that collectively influence T2D risk. To test this hypothesis we developed a method named Meta-Analysis Gene-set Enrichment of variaNT Associations (MAGENTA; http://www.broadinstitute.org/mpg/magenta). MAGENTA, in analogy to Gene Set Enrichment Analysis, tests whether sets of functionally related genes are enriched for associations with a polygenic disease or trait. MAGENTA was specifically designed to exploit the statistical power of large genome-wide association (GWA) study meta-analyses whose individual genotypes are not available. This is achieved by combining variant association p-values into gene scores and then correcting for confounders, such as gene size, variant number, and linkage disequilibrium properties. Using simulations, we determined the range of parameters for which MAGENTA can detect associations likely missed by single-marker analysis. We verified MAGENTA''s performance on empirical data by identifying known relevant pathways in lipid and lipoprotein GWA meta-analyses. We then tested our mitochondrial hypothesis by applying MAGENTA to three gene sets: nuclear regulators of mitochondrial genes, oxidative phosphorylation genes, and ∼1,000 nuclear-encoded mitochondrial genes. The analysis was performed using the most recent T2D GWA meta-analysis of 47,117 people and meta-analyses of seven diabetes-related glycemic traits (up to 46,186 non-diabetic individuals). This well-powered analysis found no significant enrichment of associations to T2D or any of the glycemic traits in any of the gene sets tested. These results suggest that common variants affecting nuclear-encoded mitochondrial genes have at most a small genetic contribution to T2D susceptibility.     

Genetic Variants in Epidermal Growth Factor Receptor Pathway Genes and Risk of Esophageal Squamous Cell Carcinoma and Gastric Cancer in a Chinese Population

The epidermal growth factor receptor (EGFR) signaling pathway regulates cell proliferation, differentiation, and survival, and is frequently dysregulated in esophageal and gastric cancers. Few studies have comprehensively examined the association between germline genetic variants in the EGFR pathway and risk of esophageal and gastric cancers. Based on a genome-wide association study in a Han Chinese population, we examined 3443 SNPs in 127 genes in the EGFR pathway for 1942 esophageal squamous cell carcinomas (ESCCs), 1758 gastric cancers (GCs), and 2111 controls. SNP-level analyses were conducted using logistic regression models. We applied the resampling-based adaptive rank truncated product approach to determine the gene- and pathway-level associations. The EGFR pathway was significantly associated with GC risk (P = 2.16×10⁻³). Gene-level analyses found 10 genes to be associated with GC, including FYN, MAPK8, MAP2K4, GNAI3, MAP2K1, TLN1, PRLR, PLCG2, RPS6KB2, and PIK3R3 (P<0.05). For ESCC, we did not observe a significant pathway-level association (P = 0.72), but gene-level analyses suggested associations between GNAI3, CHRNE, PAK4, WASL, and ITCH, and ESCC (P<0.05). Our data suggest an association between specific genes in the EGFR signaling pathway and risk of GC and ESCC. Further studies are warranted to validate these associations and to investigate underlying mechanisms.     

Genetic signature of human longevity in PKC and NF‐κB signaling

Gene variants associated with longevity are also associated with protection against cognitive decline, dementia and Alzheimer''s disease, suggesting that common physiologic pathways act at the interface of longevity and cognitive function. To test the hypothesis that variants in genes implicated in cognitive function may promote exceptional longevity, we performed a comprehensive 3‐stage study to identify functional longevity‐associated variants in ~700 candidate genes in up to 450 centenarians and 500 controls by target capture sequencing analysis. We found an enrichment of longevity‐associated genes in the nPKC and NF‐κB signaling pathways by gene‐based association analyses. Functional analysis of the top three gene variants (NFKBIA, CLU, PRKCH) suggests that non‐coding variants modulate the expression of cognate genes, thereby reducing signaling through the nPKC and NF‐κB. This matches genetic studies in multiple model organisms, suggesting that the evolutionary conservation of reduced PKC and NF‐κB signaling pathways in exceptional longevity may include humans.