An Effective Method to Identify Shared Pathways and Common Factors among Neurodegenerative Diseases

Groups of distinct but related diseases often share common symptoms, which suggest likely overlaps in underlying pathogenic mechanisms. Identifying the shared pathways and common factors among those disorders can be expected to deepen our understanding for them and help designing new treatment strategies effected on those diseases. Neurodegeneration diseases, including Alzheimer''s disease (AD), Parkinson''s disease (PD) and Huntington''s disease (HD), were taken as a case study in this research. Reported susceptibility genes for AD, PD and HD were collected and human protein-protein interaction network (hPPIN) was used to identify biological pathways related to neurodegeneration. 81 KEGG pathways were found to be correlated with neurodegenerative disorders. 36 out of the 81 are human disease pathways, and the remaining ones are involved in miscellaneous human functional pathways. Cancers and infectious diseases are two major subclasses within the disease group. Apoptosis is one of the most significant functional pathways. Most of those pathways found here are actually consistent with prior knowledge of neurodegenerative diseases except two cell communication pathways: adherens and tight junctions. Gene expression analysis showed a high probability that the two pathways were related to neurodegenerative diseases. A combination of common susceptibility genes and hPPIN is an effective method to study shared pathways involved in a group of closely related disorders. Common modules, which might play a bridging role in linking neurodegenerative disorders and the enriched pathways, were identified by clustering analysis. The identified shared pathways and common modules can be expected to yield clues for effective target discovery efforts on neurodegeneration.     

Network Analysis Identifies SOD2 mRNA as a Potential Biomarker for Parkinson's Disease

Increasing evidence indicates that Parkinson''s disease (PD) and type 2 diabetes (T2DM) share dysregulated molecular networks. We identified 84 genes shared between PD and T2DM from curated disease-gene databases. Nitric oxide biosynthesis, lipid and carbohydrate metabolism, insulin secretion and inflammation were identified as common dysregulated pathways. A network prioritization approach was implemented to rank genes according to their distance to seed genes and their involvement in common biological pathways. Quantitative polymerase chain reaction assays revealed that a highly ranked gene, superoxide dismutase 2 (SOD2), is upregulated in PD patients compared to healthy controls in 192 whole blood samples from two independent clinical trials, the Harvard Biomarker Study (HBS) and the Diagnostic and Prognostic Biomarkers in Parkinson''s disease (PROBE). The results from this study reinforce the idea that shared molecular networks between PD and T2DM provides an additional source of biologically meaningful biomarkers. Evaluation of this biomarker in de novo PD patients and in a larger prospective longitudinal study is warranted.     

Identification of Risk Pathways and Functional Modules for Coronary Artery Disease Based on Genome-wide SNP Data

Coronary artery disease(CAD) is a complex human disease, involving multiple genes and their nonlinear interactions, which often act in a modular fashion. Genome-wide single nucleotide polymorphism(SNP) profiling provides an effective technique to unravel these underlying genetic interplays or their functional involvements for CAD. This study aimed to identify the susceptible pathways and modules for CAD based on SNP omics. First, the Wellcome Trust Case Control Consortium(WTCCC) SNP datasets of CAD and control samples were used to assess the jointeffect of multiple genetic variants at the pathway level, using logistic kernel machine regression model. Then, an expanded genetic network was constructed by integrating statistical gene–gene interactions involved in these susceptible pathways with their protein–protein interaction(PPI)knowledge. Finally, risk functional modules were identified by decomposition of the network. Of 276 KEGG pathways analyzed, 6 pathways were found to have a significant effect on CAD. Other than glycerolipid metabolism, glycosaminoglycan biosynthesis, and cardiac muscle contraction pathways, three pathways related to other diseases were also revealed, including Alzheimer's disease, non-alcoholic fatty liver disease, and Huntington's disease. A genetic epistatic network of 95 genes was further constructed using the abovementioned integrative approach. Of 10 functional modules derived from the network, 6 have been annotated to phospholipase C activity and cell adhesion molecule binding, which also have known functional involvement in Alzheimer's disease.These findings indicate an overlap of the underlying molecular mechanisms between CAD and Alzheimer's disease, thus providing new insights into the molecular basis for CAD and its molecular relationships with other diseases.     

Disease gene interaction pathways: a potential framework for how disease genes associate by disease-risk modules

Background

Disease genes that interact cooperatively play crucial roles in the process of complex diseases, yet how to analyze and represent their associations is still an open problem. Traditional methods have failed to represent direct biological evidences that disease genes associate with each other in the pathogenesis of complex diseases. Molecular networks, assumed as ‘a form of biological systems’, consist of a set of interacting biological modules (functional modules or pathways) and this notion could provide a promising insight into deciphering this topic.

Methodology/Principal Findings

In this paper, we hypothesized that disease genes might associate by virtue of the associations between biological modules in molecular networks. Then we introduced a novel disease gene interaction pathway representation and analysis paradigm, and managed to identify the disease gene interaction pathway for 61 known disease genes of coronary artery disease (CAD), which contained 46 disease-risk modules and 182 interaction relationships. As demonstrated, disease genes associate through prescribed communication protocols of common biological functions and pathways.

Conclusions/Significance

Our analysis was proved to be coincident with our primary hypothesis that disease genes of complex diseases interact with their neighbors in a cooperative manner, associate with each other through shared biological functions and pathways of disease-risk modules, and finally cause dysfunctions of a series of biological processes in molecular networks. We hope our paradigm could be a promising method to identify disease gene interaction pathways for other types of complex diseases, affording additional clues in the pathogenesis of complex diseases.     

Complex Disease Interventions from a Network Model for Type 2 Diabetes

There is accumulating evidence that the proteins encoded by the genes associated with a common disorder interact with each other, participate in similar pathways and share GO terms. It has been anticipated that the functional modules in a disease related functional linkage network are informative to reveal significant metabolic processes and disease’s associations with other complex disorders. In the current study, Type 2 diabetes associated functional linkage network (T2DFN) containing 2770 proteins and 15041 linkages was constructed. The functional modules in this network were scored and evaluated in terms of shared pathways, co-localization, co-expression and associations with similar diseases. The assembly of top scoring overlapping members in the functional modules revealed that, along with the well known biological pathways, circadian rhythm, diverse actions of nuclear receptors in steroid and retinoic acid metabolisms have significant occurrence in the pathophysiology of the disease. The disease’s association with other metabolic and neuromuscular disorders was established through shared proteins. Nuclear receptor NRIP1 has a pivotal role in lipid and carbohydrate metabolism, indicating the need to investigate subsequent effects of NRIP1 on Type 2 diabetes. Our study also revealed that CREB binding protein (CREBBP) and cardiotrophin-1 (CTF1) have suggestive roles in linking Type 2 diabetes and neuromuscular diseases.     

Proteomic and bioinformatic analysis of membrane proteome in type 2 diabetic mouse liver

Type 2 diabetes mellitus (T2DM) is the most prevalent and serious metabolic disease affecting people worldwide. T2DM results from insulin resistance of the liver, muscle, and adipose tissue. In this study, we used proteomic and bioinformatic methodologies to identify novel hepatic membrane proteins that are related to the development of hepatic insulin resistance, steatosis, and T2DM. Using FT‐ICR MS, we identified 95 significantly differentially expressed proteins in the membrane fraction of normal and T2DM db/db mouse liver. These proteins are primarily involved in energy metabolism pathways, molecular transport, and cellular signaling, and many of them have not previously been reported in diabetic studies. Bioinformatic analysis revealed that 16 proteins may be related to the regulation of insulin signaling in the liver. In addition, six proteins are associated with energy stress‐induced, nine proteins with inflammatory stress‐induced, and 14 proteins with endoplasmic reticulum stress‐induced hepatic insulin resistance. Moreover, we identified 19 proteins that may regulate hepatic insulin resistance in a c‐Jun amino‐terminal kinase‐dependent manner. In addition, three proteins, 14–3‐3 protein beta (YWHAB), Slc2a4 (GLUT4), and Dlg4 (PSD‐95), are discovered by comprehensive bioinformatic analysis, which have correlations with several proteins identified by proteomics approach. The newly identified proteins in T2DM should provide additional insight into the development and pathophysiology of hepatic steatosis and insulin resistance, and they may serve as useful diagnostic markers and/or therapeutic targets for these diseases.     

Molecular differences in Alzheimer's disease between male and female patients determined by integrative network analysis

Alzheimer's disease (AD) is a complex neurodegenerative disease and the most common cause of dementia among the elderly. There has been increasing recognition of sex differences in AD prevalence, clinical manifestation, disease course and prognosis. However, there have been few studies on the molecular mechanism underlying these differences. To address this issue, we carried out global gene expression and integrative network analyses based on expression profiles (GSE84422) across 17 cortical regions of 125 individuals with AD. There were few genes that were differentially expressed across the 17 regions between the two sexes, with only four (encoding glutamate metabotropic receptor 2, oestrogen‐related receptor beta, kinesin family member 26B, and aspartoacylase) that were differentially expressed in three regions. A pan‐cortical brain region co‐expression network analysis identified pathways and genes (eg, glycogen synthase kinase 3β) that were significantly associated with clinical characteristics of AD (such as neurofibrillary score) in males only. Similarity analyses between region‐specific networks indicated that male patients exhibited greater variability, especially in the superior parietal lobule, dorsolateral prefrontal cortex and occipital visual cortex. A network module analysis revealed an association between clinical traits and crosstalk of sex‐specific modules. An examination of temporal and spatial patterns of sex differences in AD showed that molecular networks were more conserved in females than in males in different cortical regions and at different AD stages. These findings provide insight into critical molecular pathways governing sex differences in AD pathology.     

Identification of gene expression profiles and key genes in subchondral bone of osteoarthritis using weighted gene coexpression network analysis

Osteoarthritis (OA) significantly influences the quality life of people around the world. It is urgent to find an effective way to understand the genetic etiology of OA. We used weighted gene coexpression network analysis (WGCNA) to explore the key genes involved in the subchondral bone pathological process of OA. Fifty gene expression profiles of GSE51588 were downloaded from the Gene Expression Omnibus database. The OA‐associated genes and gene ontologies were acquired from JuniorDoc. Weighted gene coexpression network analysis was used to find disease‐related networks based on 21756 gene expression correlation coefficients, hub‐genes with the highest connectivity in each module were selected, and the correlation between module eigengene and clinical traits was calculated. The genes in the traits‐related gene coexpression modules were subject to functional annotation and pathway enrichment analysis using ClusterProfiler. A total of 73 gene modules were identified, of which, 12 modules were found with high connectivity with clinical traits. Five modules were found with enriched OA‐associated genes. Moreover, 310 OA‐associated genes were found, and 34 of them were among hub‐genes in each module. Consequently, enrichment results indicated some key metabolic pathways, such as extracellular matrix (ECM)‐receptor interaction (hsa04512), focal adhesion (hsa04510), the phosphatidylinositol 3'‐kinase (PI3K)‐Akt signaling pathway (PI3K‐AKT) (hsa04151), transforming growth factor beta pathway, and Wnt pathway. We intended to identify some core genes, collagen (COL)6A3, COL6A1, ITGA11, BAMBI, and HCK, which could influence downstream signaling pathways once they were activated. In this study, we identified important genes within key coexpression modules, which associate with a pathological process of subchondral bone in OA. Functional analysis results could provide important information to understand the mechanism of OA.     

Identification of prognostic biomarkers for breast cancer brain metastases based on the bioinformatics analysis

PurposeThe prognosis of breast cancer (BC) patients who develop into brain metastases (BMs) is very poor. Thus, it is of great significance to explore the etiology of BMs in BC and identify the key genes involved in this process to improve the survival of BC patients with BMs.Patients and methodsThe gene expression data and the clinical information of BC patients were downloaded from TCGA and GEO database. Differentially expressed genes (DEGs) in TCGA-BRCA and GSE12276 were overlapped to find differentially expressed metastatic genes (DEMGs). The protein-protein interaction (PPI) network of DEMGs was constructed via STRING database. ClusterProfiler R package was applied to perform the gene ontology (GO) enrichment analysis of DEMGs. The univariate Cox regression analysis and the Kaplan-Meier (K-M) curves were plotted to screen DEMGs associated with the overall survival and the metastatic recurrence survival, which were identified as the key genes associated with the BMs in BC. The immune infiltration and the expressions of immune checkpoints for BC patients with brain relapses and BC patients with other relapses were analyzed respectively. The correlations among the expressions of key genes and the differently infiltrated immune cells or the differentially expressed immune checkpoints were calculated. The gene set enrichment analysis (GSEA) of each key gene was conducted to investigate the potential mechanisms of key genes involved in BC patients with BMs. Moreover, CTD database was used to predict the drug-gene interaction network of key genes.ResultsA total of 154 DEGs were identified in BC patients at M0 and M1 in TCGA database. A total of 667 DEGs were identified in BC patients with brain relapses and with other relapses. By overlapping these DEGs, 17 DEMGs were identified, which were enriched in the cell proliferation related biological processes and the immune related molecular functions. The univariate Cox regression analysis and the Kaplan-Meier curves revealed that CXCL9 and GPR171 were closely associated with the overall survival and the metastatic recurrence survival and were identified as key genes associated with BMs in BC. The analyses of immune infiltration and immune checkpoint expressions showed that there was a significant difference of the immune microenvironment between brain relapses and other relapses in BC. GSEA indicated that CXCL9 and GPR171 may regulate BMs in BC via the immune-related pathways.ConclusionOur study identified the key genes associated with BMs in BC patients and explore the underlying mechanisms involved in the etiology of BMs in BC. These findings may provide a promising approach for the treatments of BC patients with BMs.     

Associations between HIV and human pathways revealed by protein-protein interactions and correlated gene expression profiles

Background

AIDS is one of the most devastating diseases in human history. Decades of studies have revealed host factors required for HIV infection, indicating that HIV exploits host processes for its own purposes. HIV infection leads to AIDS as well as various comorbidities. The associations between HIV and human pathways and diseases may reveal non-obvious relationships between HIV and non-HIV-defining diseases.

Principal Findings

Human biological pathways were evaluated and statistically compared against the presence of HIV host factor related genes. All of the obtained scores comparing HIV targeted genes and biological pathways were ranked. Different rank results based on overlapping genes, recovered virus-host interactions, co-expressed genes, and common interactions in human protein-protein interaction networks were obtained. Correlations between rankings suggested that these measures yielded diverse rankings. Rank combination of these ranks led to a final ranking of HIV-associated pathways, which revealed that HIV is associated with immune cell-related pathways and several cancer-related pathways. The proposed method is also applicable to the evaluation of associations between other pathogens and human pathways and diseases.

Conclusions

Our results suggest that HIV infection shares common molecular mechanisms with certain signaling pathways and cancers. Interference in apoptosis pathways and the long-term suppression of immune system functions by HIV infection might contribute to tumorigenesis. Relationships between HIV infection and human pathways of disease may aid in the identification of common drug targets for viral infections and other diseases.     

Carbohydrate metabolic pathway genes associated with quantitative trait loci (QTL) for obesity and type 2 diabetes: Identification by data mining

Increasing consumption of refined carbohydrates is now being recognized as a primary contributor to the development of nutritionally related chronic diseases such as obesity and type 2 diabetes mellitus (T2DM). A data mining approach was used to evaluate the role of carbohydrate metabolic pathway genes in the development of obesity and T2DM. Data from public databases were used to map the position of the carbohydrate metabolic pathway genes to known quantitative trait loci (QTL) for obesity and T2DM and for examining the pathway genes for the presence of sequence and structural genetic variants such as single nucleotide polymorphisms (SNPs) and copy number variants (CNS), respectively. The results demonstrated that a majority of the genes of the carbohydrate metabolic pathways are associated with QTL for obesity and many for T2DM. In addition, some key genes of the pathways also encode non-synonymous SNPs that exhibit significant differences in population frequencies. This study emphasizes the significance of the metabolic pathways genes in the development of disease phenotypes, its differential occurrence across populations and between individuals, and a strategy for interpreting an individuals' risk for disease.     

A shared comparison of diabetes mellitus and neurodegenerative disorders

Diabetes mellitus (DM), one of the most prevalent metabolic diseases in the world population, is associated with a number of comorbid conditions including obesity, pancreatic endocrine changes, and renal and cardio-cerebrovascular alterations, coupled with peripheral neuropathy and neurodegenerative disease, some of these disorders are bundled into metabolic syndrome. Type 1 DM (T1DM) is an autoimmune disease that destroys the insulin-secreting islet cells. Type 2 DM (T2DM) is diabetes that is associated with an imbalance in the glucagon/insulin homeostasis that leads to the formation of amyloid deposits in the brain, pancreatic islet cells, and possibly in the kidney glomerulus. There are several layers of molecular pathologic alterations that contribute to the DM metabolic pathophysiology and its associated neuropathic manifestations. In this review, we describe the general signature metabolic features of DM and the cross-talk with neurodegeneration. We will assess the underlying molecular key players associated with DM-induced neuropathic disorders that are associated with both T1DM and T2DM. In this context, we will highlight the role of tau and amyloid protein deposits in the brain as well in the pancreatic islet cells, and possibly in the kidney glomerulus. Furthermore, we will discuss the central role of mitochondria, oxidative stress, and the unfolded protein response in mediating the DM-associated neuropathic degeneration. This study will elucidate the relationship between DM and neurodegeneration which may account for the evolution of other neurodegenerative diseases, particularly Alzheimer's disease and Parkinson's disease as discussed later.     

Shared molecular and functional frameworks among five complex human disorders: a comparative study on interactomes linked to susceptibility genes

Background

Genome-wide association studies (gwas) are invaluable in revealing the common variants predisposing to complex human diseases. Yet, until now, the large volumes of data generated from such analyses have not been explored extensively enough to identify the molecular and functional framework hosting the susceptibility genes.

Methodology/Principal Findings

We investigated the relationships among five neurodegenerative and/or autoimmune complex human diseases (Parkinson''s disease-Park, Alzheimer''s disease-Alz, multiple sclerosis-MS, rheumatoid arthritis-RA and Type 1 diabetes-T1D) by characterising the interactomes linked to their gwas-genes. An initial study on the MS interactome indicated that several genes predisposing to the other autoimmune or neurodegenerative disorders may come into contact with it, suggesting that susceptibility to distinct diseases may converge towards common molecular and biological networks. In order to test this hypothesis, we performed pathway enrichment analyses on each disease interactome independently. Several issues related to immune function and growth factor signalling pathways appeared in all autoimmune diseases, and, surprisingly, in Alzheimer''s disease. Furthermore, the paired analyses of disease interactomes revealed significant molecular and functional relatedness among autoimmune diseases, and, unexpectedly, between T1D and Alz.

Conclusions/Significance

The systems biology approach highlighted several known pathogenic processes, indicating that changes in these functions might be driven or sustained by the framework linked to genetic susceptibility. Moreover, the comparative analyses among the five genetic interactomes revealed unexpected genetic relationships, which await further biological validation. Overall, this study outlines the potential of systems biology to uncover links between genetics and pathogenesis of complex human disorders.     

Detecting disease associated modules and prioritizing active genes based on high throughput data

Background  

The accumulation of high-throughput data greatly promotes computational investigation of gene function in the context of complex biological systems. However, a biological function is not simply controlled by an individual gene since genes function in a cooperative manner to achieve biological processes. In the study of human diseases, rather than to discover disease related genes, identifying disease associated pathways and modules becomes an essential problem in the field of systems biology.     

Genetic evidence for common pathways in human age‐related diseases

Aging is the single largest risk factor for chronic disease. Studies in model organisms have identified conserved pathways that modulate aging rate and the onset and progression of multiple age‐related diseases, suggesting that common pathways of aging may influence age‐related diseases in humans as well. To determine whether there is genetic evidence supporting the notion of common pathways underlying age‐related diseases, we analyzed the genes and pathways found to be associated with five major categories of age‐related disease using a total of 410 genomewide association studies (GWAS). While only a small number of genes are shared among all five disease categories, those found in at least three of the five major age‐related disease categories are highly enriched for apoliprotein metabolism genes. We found that a more substantial number of gene ontology (GO) terms are shared among the 5 age‐related disease categories and shared GO terms include canonical aging pathways identified in model organisms, such as nutrient‐sensing signaling, translation, proteostasis, stress responses, and genome maintenance. Taking advantage of the vast amount of genetic data from the GWAS, our findings provide the first direct evidence that conserved pathways of aging simultaneously influence multiple age‐related diseases in humans as has been demonstrated in model organisms.     

A comprehensive approach to the molecular determinants of lifespan using a Boolean model of geroconversion

Altered molecular responses to insulin and growth factors (GF) are responsible for late‐life shortening diseases such as type‐2 diabetes mellitus (T2DM) and cancers. We have built a network of the signaling pathways that control S‐phase entry and a specific type of senescence called geroconversion. We have translated this network into a Boolean model to study possible cell phenotype outcomes under diverse molecular signaling conditions. In the context of insulin resistance, the model was able to reproduce the variations of the senescence level observed in tissues related to T2DM's main morbidity and mortality. Furthermore, by calibrating the pharmacodynamics of mTOR inhibitors, we have been able to reproduce the dose‐dependent effect of rapamycin on liver degeneration and lifespan expansion in wild‐type and HER2–neu mice. Using the model, we have finally performed an in silico prospective screen of the risk–benefit ratio of rapamycin dosage for healthy lifespan expansion strategies. We present here a comprehensive prognostic and predictive systems biology tool for human aging.     

Network medicine: a network-based approach to human disease

Given the functional interdependencies between the molecular components in a human cell, a disease is rarely a consequence of an abnormality in a single gene, but reflects the perturbations of the complex intracellular and intercellular network that links tissue and organ systems. The emerging tools of network medicine offer a platform to explore systematically not only the molecular complexity of a particular disease, leading to the identification of disease modules and pathways, but also the molecular relationships among apparently distinct (patho)phenotypes. Advances in this direction are essential for identifying new disease genes, for uncovering the biological significance of disease-associated mutations identified by genome-wide association studies and full-genome sequencing, and for identifying drug targets and biomarkers for complex diseases.