Pathway-based analysis tools for complex diseases: a review

Beyond genome-wide significance: integrative approaches to the interpretation and extension of GWAS findings for alcohol use disorder

Alcohol use disorder (AUD) is a heritable complex behavior. Due to the highly polygenic nature of AUD, identifying genetic variants that comprise this heritable variation has proved to be challenging. With the exception of functional variants in alcohol metabolizing genes (e.g. ADH1B and ALDH2), few other candidate loci have been confidently linked to AUD. Genome-wide association studies (GWAS) of AUD and other alcohol-related phenotypes have either produced few hits with genome-wide significance or have failed to replicate on further study. These issues reinforce the complex nature of the genetic underpinnings for AUD and suggest that both GWAS studies with larger samples and additional analysis approaches that better harness the nominally significant loci in existing GWAS are needed. Here, we review approaches of interest in the post-GWAS era, including in silico functional analyses; functional partitioning of single nucleotide polymorphism heritability; aggregation of signal into genes and gene networks; and validation of identified loci, genes and gene networks in postmortem brain tissue and across species. These integrative approaches hold promise to illuminate our understanding of the biological basis of AUD; however, we recognize that the main challenge continues to be the extremely polygenic nature of AUD, which necessitates large samples to identify multiple loci associated with AUD liability.

Pathway-based classification of genetic diseases

Background

In medical genetics, diseases are classified according to the nature (hypothetical nature) of the underlying genetic defect. The classification is “gene-centric” and “factor-centric”; a disease may be, thereby, designated as monogenic, oligogenic or polygenic/multifactorial. Chromosomal diseases/syndromes and abnormalities are generally considered apart from these designations due to distinctly different formation mechanisms and simultaneous encompassing from several to several hundreds of co-localized genes. These definitions are ubiquitously used and are perfectly suitable for human genetics issues in historical and academic perspective. However, recent achievements in systems biology have offered a possibility to explore the consequences of a genetic defect from genomic variations to molecular/cellular pathway alterations unique to a disease. Since pathogenetic mechanisms (pathways) are more influential on our understating of disease presentation and progression than genetic defects per se, a need for a disease classification reflecting both genetic causes and molecular/cellular mechanisms appears to exist. Here, we propose an extension to the common disease classification based on the underlying genetic defects, which focuses on disease-specific molecular pathways.

Conclusion

The basic idea of our classification is to propose pathways as parameters for designating a genetic disease. To proceed, we have followed the tradition of using ancient Greek words and prefixes to create the terms for the pathway-based classification of genetic diseases. We have chosen the word “griphos” (γρῖφος), which simultaneously means “net” and “puzzle”, accurately symbolizing the term “pathway” currently used in molecular biology and medicine. Thus, diseases may be classified as monogryphic (single pathway is altered to result in a phenotype), digryphic (two pathways are altered to result in a phenotype), etc.; additionally, diseases may be designated as oligogryphic (several pathways are altered to result in a phenotype), polygryphic (numerous pathways or cascades of pathways are altered to result in a phenotype) and homeogryphic in cases of comorbid diseases resulted from shared pathway alterations. We suppose that classifying illness this way using both “gene-centric” and “pathway-centric” concepts is able to revolutionize current views on genetic diseases.

Computational Methods to Investigate the Impact of miRNAs on Pathways

Pathway analysis is a wide class of methods allowing to determine the alteration of functional processes in complex diseases. However, biological pathways are still partial, and knowledge coming from posttranscriptional regulators has started to be considered in a systematic way only recently. Here we will give a global and updated view of the main pathway and subpathway analysis methodologies, focusing on the improvements obtained through the recent introduction of microRNAs as regulatory elements in these frameworks.

PASNet: pathway-associated sparse deep neural network for prognosis prediction from high-throughput data

Background

Predicting prognosis in patients from large-scale genomic data is a fundamentally challenging problem in genomic medicine. However, the prognosis still remains poor in many diseases. The poor prognosis may be caused by high complexity of biological systems, where multiple biological components and their hierarchical relationships are involved. Moreover, it is challenging to develop robust computational solutions with high-dimension, low-sample size data.

Results

In this study, we propose a Pathway-Associated Sparse Deep Neural Network (PASNet) that not only predicts patients’ prognoses but also describes complex biological processes regarding biological pathways for prognosis. PASNet models a multilayered, hierarchical biological system of genes and pathways to predict clinical outcomes by leveraging deep learning. The sparse solution of PASNet provides the capability of model interpretability that most conventional fully-connected neural networks lack. We applied PASNet for long-term survival prediction in Glioblastoma multiforme (GBM), which is a primary brain cancer that shows poor prognostic performance. The predictive performance of PASNet was evaluated with multiple cross-validation experiments. PASNet showed a higher Area Under the Curve (AUC) and F1-score than previous long-term survival prediction classifiers, and the significance of PASNet’s performance was assessed by Wilcoxon signed-rank test. Furthermore, the biological pathways, found in PASNet, were referred to as significant pathways in GBM in previous biology and medicine research.

Conclusions

PASNet can describe the different biological systems of clinical outcomes for prognostic prediction as well as predicting prognosis more accurately than the current state-of-the-art methods. PASNet is the first pathway-based deep neural network that represents hierarchical representations of genes and pathways and their nonlinear effects, to the best of our knowledge. Additionally, PASNet would be promising due to its flexible model representation and interpretability, embodying the strengths of deep learning. The open-source code of PASNet is available at https://github.com/DataX-JieHao/PASNet.

Pathway networks generated from human disease phenome

Background

Understanding the effect of human genetic variations on disease can provide insight into phenotype-genotype relationships, and has great potential for improving the effectiveness of personalized medicine. While some genetic markers linked to disease susceptibility have been identified, a large number are still unknown. In this paper, we propose a pathway-based approach to extend disease-variant associations and find new molecular connections between genetic mutations and diseases.

Methods

We used a compilation of over 80,000 human genetic variants with known disease associations from databases including the Online Mendelian Inheritance in Man (OMIM), Clinical Variance database (ClinVar), Universal Protein Resource (UniProt), and Human Gene Mutation Database (HGMD). Furthermore, we used the Unified Medical Language System (UMLS) to normalize variant phenotype terminologies, mapping 87% of unique genetic variants to phenotypic disorder concepts. Lastly, variants were grouped by UMLS Medical Subject Heading (MeSH) identifiers to determine pathway enrichment in Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways.

Results

By linking KEGG pathways through underlying variant associations, we elucidated connections between the human genetic variant-based disease phenome and metabolic pathways, finding novel disease connections not otherwise detected through gene-level analysis. When looking at broader disease categories, our network analysis showed that large complex diseases, such as cancers, are highly linked by their common pathways. In addition, we found Cardiovascular Diseases and Skin and Connective Tissue Diseases to have the highest number of common pathways, among 35 significant main disease category (MeSH) pairings.

Conclusions

This study constitutes an important contribution to extending disease-variant connections and new molecular links between diseases. Novel disease connections were made by disease-pathway associations not otherwise detected through single-gene analysis. For instance, we found that mutations in different genes associated to Noonan Syndrome and Essential Hypertension share a common pathway. This analysis also provides the foundation to build novel disease-drug networks through their underlying common metabolic pathways, thus enabling new diagnostic and therapeutic interventions.

Electronic supplementary material

The online version of this article (10.1186/s12920-018-0386-2) contains supplementary material, which is available to authorized users.

Shared Biological Pathways Between Alzheimer's Disease and Ischemic Stroke

Alzheimer's disease (AD) and ischemic stroke (IS) are an immense socioeconomic burden worldwide. There is a possibility that shared genetic factors lead to their links at epidemiological and pathophysiological levels. Although recent genome-wide association studies (GWAS) have provided profound insights into the genetics of AD and IS, no shared genetic variants have been identified to date. This prompted us to initiate this study, which sought to identify shared pathways linking AD and IS. We took advantage of large-scale GWAS summary data of AD (17,008 AD cases and 37,154 controls) and IS (10,307 cases and 19,326 controls) to conduct pathway analyses using genetic pathways from multiple well-studied databases, including GO, KEGG, PANTHER, Reactome, and Wikipathways. Collectively, we discovered that AD and IS shared 179 GO categories (56 biological processes, 95 cellular components, and 28 molecular functions); and the following pathways: six KEGG pathways; two PANTHER pathways; four Reactome pathways; and one in Wikipathways pathway. The more fine-grained GO terms were mainly summarized into different functional categories: transcriptional and post-transcriptional regulation, synapse, endocytic membrane traffic through the endosomal system, signaling transduction, immune process, multi-organism process, protein catabolic metabolism, and cell adhesion. The shared pathways were roughly classified into three categories: immune system; cancer (NSCLC and glioma); and signal transduction pathways involving the cadherin signaling pathway, Wnt signaling pathway, G-protein signaling and downstream signaling mediated by phosphoinositides (PIPs). The majority of these common pathways linked to both AD and IS were supported by convincing evidence from the literature. In conclusion, our findings contribute to a better understanding of common biological mechanisms underlying AD and IS and serve as a guide to direct future research.

Complex Haplotypes of GSTM1 Gene Deletions Harbor Signatures of a Selective Sweep in East Asian Populations

The deletion of the metabolizing Glutathione S-transferase Mu 1 (GSTM1) gene has been associated with multiple cancers, metabolic and autoimmune disorders, as well as drug response. It is unusually common, with allele frequency reaching up to 75% in some human populations. Such high allele frequency of a derived allele with apparent impact on an otherwise conserved gene is a rare phenomenon. To investigate the evolutionary history of this locus, we analyzed 310 genomes using population genetics tools. Our analysis revealed a surprising lack of linkage disequilibrium between the deletion and the flanking single nucleotide variants in this locus. Tests that measure extended homozygosity and rapid change in allele frequency revealed signatures of an incomplete sweep in the locus. Using empirical approaches, we identified the Tanuki haplogroup, which carries the GSTM1 deletion and is found in approximately 70% of East Asian chromosomes. This haplogroup has rapidly increased in frequency in East Asian populations, contributing to a high population differentiation among continental human groups. We showed that extended homozygosity and population differentiation for this haplogroup is incompatible with simulated neutral expectations in East Asian populations. In parallel, we found that the Tanuki haplogroup is significantly associated with the expression levels of other GSTM genes. Collectively, our results suggest that standing variation in this locus has likely undergone an incomplete sweep in East Asia with regulatory impact on multiple GSTM genes. Our study provides the necessary framework for further studies to elucidate the evolutionary reasons that maintain disease-susceptibility variants in the GSTM1 locus.

Exposure to a firefighting overhaul environment without respiratory protection increases immune dysregulation and lung disease risk

Firefighting activities appear to increase the risk of acute and chronic lung disease, including malignancy. While self-contained breathing apparatuses (SCBA) mitigate exposures to inhalable asphyxiates and carcinogens, firefighters frequently remove SCBA during overhaul when the firegrounds appear clear of visible smoke. Using a mouse model of overhaul without airway protection, the impact of fireground environment exposure on lung gene expression was assessed to identify transcripts potentially critical to firefighter-related chronic pulmonary illnesses. Lung tissue was collected 2 hrs post-overhaul and evaluated via whole genome transcriptomics by RNA-seq. Although gas metering showed that the fireground overhaul levels of carbon monoxide (CO), carbon dioxide (CO2), hydrogen cyanine (HCN), hydrogen sulfide (H2S) and oxygen (O2) were within NIOSH ceiling recommendations, 3852 lung genes were differentially expressed when mice exposed to overhaul were compared to mice on the fireground but outside the overhaul environment. Importantly, overhaul exposure was associated with an up/down-regulation of 86 genes with a fold change of 1.5 or greater (p<0.5) including the immunomodulatory-linked genes S100a8 and Tnfsf9 (downregulation) and the cancer-linked genes, Capn11 and Rorc (upregulation). Taken together these findings indicate that, without respiratory protection, exposure to the fireground overhaul environment is associated with transcriptional changes impacting proteins potentially related to inflammation-associated lung disease and cancer.

Bidirectional Transcriptome Analysis of Rat Bone Marrow-Derived Mesenchymal Stem Cells and Activated Microglia in an In Vitro Coculture System

Microglia contribute to the regulation of neuroinflammation and play an important role in the pathogenesis of brain diseases. Thus, regulation of neuroinflammation triggered by activated microglia in brain diseases has become a promising curative strategy. Bone marrow-derived mesenchymal stem cells (BM-MSCs) have been shown to have therapeutic effects, resulting from the regulation of inflammatory conditions in the brain. In this study, we investigated differential gene expression in rat BM-MSCs (rBM-MSCs) that were cocultured with lipopolysaccharide- (LPS-) stimulated primary rat microglia using microarray analysis and evaluated the functional relationships through Ingenuity Pathway Analysis (IPA). We also evaluated the effects of rBM-MSC on LPS-stimulated microglia using a reverse coculture system and the same conditions of the transcriptomic analysis. In the transcriptome of rBM-MSCs, 67 genes were differentially expressed, which were highly related with migration of cells, compared to control. The prediction of the gene network using IPA and experimental validation showed that LPS-stimulated primary rat microglia increase the migration of rBM-MSCs. Reversely, expression patterns of the transcriptome in LPS-stimulated primary rat microglia were changed when cocultured with rBM-MSCs. Our results showed that 65 genes were changed, which were highly related with inflammatory response, compared to absence of rBM-MSCs. In the same way with the aforementioned, the prediction of the gene network and experimental validation showed that rBM-MSCs decrease the inflammatory response of LPS-stimulated primary rat microglia. Our data indicate that LPS-stimulated microglia increase the migration of rBM-MSCs and that rBM-MSCs reduce the inflammatory activity in LPS-stimulated microglia. The results of this study show complex mechanisms underlying the interaction between rBM-MSCs and activated microglia and may be helpful for the development of stem cell-based strategies for brain diseases.

Overexpression of matrix metalloproteinase-9 in adolescents with primary spontaneous pneumothorax for surgical intervention

OBJECTIVE

To determine gene expression profiles associated with bullae formation in patients with primary spontaneous pneumothorax and to identify candidate genes associated with surgical intervention.

METHODS

Twenty-four adolescents with primary spontaneous pneumothorax were enrolled prospectively. A global gene expression analysis of 9 paired lung biopsies (lesion and normal adjacent sites) was performed to identify differentially expressed genes associated with spontaneous pneumothorax. Pathway and network analysis was performed using the Database for Annotation, Visualization and Integrated Discovery web tool. Candidate genes and encoding proteins were assessed in blood samples and compared between patients with pneumothorax and healthy control patients.

RESULTS

A total of 1519 differentially expressed transcripts corresponding to known genes were identified comparing the lesion lung with paired adjacent normal lung. The altered genes were mainly associated with focal adhesion and extracellular matrix-receptor interaction pathways. Genes involved in proteolysis and peptidase activity were up-regulated predominantly, especially matrix metalloproteinase-1 and -9 genes. Compared with the recovery stage, blood levels of matrix metalloproteinase-9/tissue inhibitor of metalloproteinase-1 were increased at the acute stage in patients with pneumothorax and, when compared between patients treated operatively with those treated nonoperatively, were also significantly greater. In addition, ratios of their serum levels were significantly greater in patients with pneumothorax compared with healthy control patients. Furthermore, matrix metalloproteinase-9 was predominantly overexpressed in neutrophils, alveolar macrophages, and mesothelial cells of lung biopsies.

CONCLUSIONS

An imbalance of cell-extracellular matrix interactions appears to be associated with primary spontaneous pneumothorax. Matrix metalloproteinase-9 overexpression may particularly play a role in contributing to pleural porosity for surgical intervention.

Commonalities in Biological Pathways, Genetics, and Cellular Mechanism between Alzheimer Disease and Other Neurodegenerative Diseases: An In Silico-Updated Overview

BACKGROUND

Alzheimer's disease (AD) is the most common and well-studied neurodegenerative disease (ND). Biological pathways, pathophysiology and genetics of AD show commonalities with other NDs viz. Parkinson's disease (PD), Amyotrophic lateral sclerosis (ALS), Huntington's disease (HD), Prion disease and Dentatorubral-pallidoluysian atrophy (DRPLA). Many of the NDs, sharing the common features and molecular mechanisms suggest that pathology may be directly comparable and be implicated in disease prevention and development of highly effective therapies.

METHOD

In this review, a brief description of pathophysiology, clinical symptoms and available treatment of various NDs have been explored with special emphasis on AD. Commonalities in these fatal NDs provide support for therapeutic advancements and enhance the understanding of disease manifestation.

CONCLUSION

The studies concentrating on the commonalities in biological pathways, cellular mechanisms and genetics may provide the scope to researchers to identify few novel common target(s) for disease prevention and development of effective common drugs for multi-neurodegenerative diseases.

A computational framework for complex disease stratification from multiple large-scale datasets

BACKGROUND

Multilevel data integration is becoming a major area of research in systems biology. Within this area, multi-'omics datasets on complex diseases are becoming more readily available and there is a need to set standards and good practices for integrated analysis of biological, clinical and environmental data. We present a framework to plan and generate single and multi-'omics signatures of disease states.

METHODS

The framework is divided into four major steps: dataset subsetting, feature filtering, 'omics-based clustering and biomarker identification.

RESULTS

We illustrate the usefulness of this framework by identifying potential patient clusters based on integrated multi-'omics signatures in a publicly available ovarian cystadenocarcinoma dataset. The analysis generated a higher number of stable and clinically relevant clusters than previously reported, and enabled the generation of predictive models of patient outcomes.

CONCLUSIONS

This framework will help health researchers plan and perform multi-'omics big data analyses to generate hypotheses and make sense of their rich, diverse and ever growing datasets, to enable implementation of translational P4 medicine.

REDESIGN: RDF-based Differential Signaling Framework for Precision Medicine Analytics

Pathway-based analysis holds promise to be instrumental in precision and personalized medicine analytics. However, the majority of pathway-based analysis methods utilize "fixed" or "rigid" data sets that limit their ability to account for complex biological inter-dependencies. Here, we present REDESIGN: RDF-based Differential Signaling Pathway informatics framework. The distinctive feature of the REDESIGN is that it is designed to run on "flexible" ontology-enabled data sets of curated signal transduction pathway maps to uncover high explanatory differential pathway mechanisms on gene-to-gene level. The experiments on two morphoproteomic cases demonstrated REDESIGN's capability to generate actionable hypotheses in precision/personalized medicine analytics.

Molecular hydrogen upregulates heat shock response and collagen biosynthesis, and downregulates cell cycles: meta-analyses of gene expression profiles

Molecular hydrogen exerts its effect on multiple pathologies, including oxidative stress, inflammation, and apoptosis. However, its molecular mechanisms have not been fully elucidated. In order to explore the effects of molecular hydrogen, we meta-analysed gene expression profiles modulated by molecular hydrogen. We performed microarray analysis of the mouse liver with or without drinking hydrogen water. We also integrated two previously reported microarray datasets of the rat liver into meta-analyses. We used two categories of meta-analysis methods: the cross-platform method and the conventional meta-analysis method (Fisher's method). For each method, hydrogen-modulated pathways were analysed by (i) the hypergeometric test (HGT) in the class of over-representation analysis (ORA), (ii) the gene set enrichment analysis (GSEA) in the class of functional class scoring (FCS), and (iii) the signalling pathway impact analysis (SPIA), pathway regulation score (PRS), and others in the class of pathway topology-based approach (PTA). Pathways in the collagen biosynthesis and the heat-shock response were up-regulated according to (a) HGT with the cross-platform method, (b) GSEA with the cross-platform method, and (c) PRS with the cross-platform method. Pathways in cell cycles were down-regulated according to (a) HGT with the cross-platform method, (b) GSEA with the cross-platform method, and (d) GSEA with the conventional meta-analysis method. Because the heat-shock response leads to up-regulation of collagen biosynthesis and a transient arrest of cell cycles, induction of the heat-shock response is likely to be a primary event induced by molecular hydrogen in the liver of wild-type rodents.

Nuclear receptor and VEGF pathways for gene-blood lead interactions, on bone mineral density, in Korean smokers

Osteoporosis has a complex etiology and is considered a multifactorial polygenic disease, in which genetic determinants are modulated by hormonal, lifestyle, environmental, and nutritional factors. Therefore, investigating these multiple factors, and the interactions between them, might lead to a better understanding of osteoporosis pathogenesis, and possible therapeutic interventions. The objective of this study was to identify the relationship between three blood metals (Pb, Cd, and Al), in smoking and nonsmoking patients' sera, and prevalence of osteoporosis. In particular, we focused on gene-environment interactions of metal exposure, including a dataset obtained through genome-wide association study (GWAS). Subsequently, we conducted a pathway-based analysis, using a GWAS dataset, to elucidate how metal exposure influences susceptibility to osteoporosis. In this study, we evaluated blood metal exposures for estimating the prevalence of osteoporosis in 443 participants (aged 53.24 ± 8.29), from the Republic of Korea. Those analyses revealed a negative association between lead blood levels and bone mineral density in current smokers (p trend <0.01). By further using GWAS-based pathway analysis, we found nuclear receptor (FDR<0.05) and VEGF pathways (FDR<0.05) to be significantly upregulated by blood lead burden, with regard to the prevalence of osteoporosis, in current smokers. These findings suggest that the intracellular pathways of angiogenesis and nuclear hormonal signaling can modulate interactions between lead exposure and genetic variation, with regard to susceptibility to diminished bone mineral density. Our findings may provide new leads for understanding the mechanisms underlying the development of osteoporosis, including possible interventions.

Nonmonotonic Pathway Gene Expression Analysis Reveals Oncogenic Role of p27/Kip1 at Intermediate Dose

The mechanistic basis by which the level of p27^Kip1 expression influences tumor aggressiveness and patient mortality remains unclear. To elucidate the competing tumor-suppressing and oncogenic effects of p27^Kip1 on gene expression in tumors, we analyzed the transcriptomes of squamous cell papilloma derived from Cdkn1b nullizygous, heterozygous, and wild-type mice. We developed a novel functional pathway analysis method capable of testing directional and nonmonotonic dose response. This analysis can reveal potential causal relationships that might have been missed by other nondirectional pathway analysis methods. Applying this method to capture dose-response curves in papilloma gene expression data, we show that several known cancer pathways are dominated by low-high-low gene expression responses to increasing p27 gene doses. The oncogene cyclin D1, whose expression is elevated at an intermediate p27 dose, is the most responsive gene shared by these cancer pathways. Therefore, intermediate levels of p27 may promote cellular processes favoring tumorigenesis-strikingly consistent with the dominance of heterozygous mutations in CDKN1B seen in human cancers. Our findings shed new light on regulatory mechanisms for both pro- and anti-tumorigenic roles of p27^Kip1. Functional pathway dose-response analysis provides a unique opportunity to uncover nonmonotonic patterns in biological systems.

Germline variation in inflammation-related pathways and risk of Barrett's oesophagus and oesophageal adenocarcinoma

OBJECTIVE

Oesophageal adenocarcinoma (OA) incidence has risen sharply in Western countries over recent decades. Local and systemic inflammation is considered an important contributor to OA pathogenesis. Established risk factors for OA and its precursor, Barrett's oesophagus (BE), include symptomatic reflux, obesity and smoking. The role of inherited genetic susceptibility remains an area of active investigation. Here, we explore whether germline variation related to inflammatory processes influences susceptibility to BE/OA.

DESIGN

We used data from a genomewide association study of 2515 OA cases, 3295 BE cases and 3207 controls. Our analysis included 7863 single-nucleotide polymorphisms (SNPs) in 449 genes assigned to five pathways: cyclooxygenase (COX), cytokine signalling, oxidative stress, human leucocyte antigen and nuclear factor-κB. A principal components-based analytic framework was employed to evaluate pathway-level and gene-level associations with disease risk.

RESULTS

We identified a significant signal for the COX pathway in relation to BE risk (p=0.0059, false discovery rate q=0.03), and in gene-level analyses found an association with microsomal glutathione-S-transferase 1 (MGST1); (p=0.0005, q=0.005). Assessment of 36 MGST1 SNPs identified 14 variants associated with elevated BE risk (q<0.05). Four of these were subsequently confirmed (p<5.5×10-5) in a meta-analysis encompassing an independent set of 1851 BE cases and 3496 controls, and are known strong expression quantitative trait loci for MGST1. Three such variants were associated with similar elevations in OA risk.

CONCLUSIONS

This study provides the most comprehensive evaluation of inflammation-related germline variation in relation to risk of BE/OA and suggests that variants in MGST1 influence disease susceptibility.

Transcriptome profiling from adipose tissue during a low-calorie diet reveals predictors of weight and glycemic outcomes in obese, nondiabetic subjects

Background: A low-calorie diet (LCD) reduces fat mass excess, improves insulin sensitivity, and alters adipose tissue (AT) gene expression, yet the relation with clinical outcomes remains unclear.Objective: We evaluated AT transcriptome alterations during an LCD and the association with weight and glycemic outcomes both at LCD termination and 6 mo after the LCD.Design: Using RNA sequencing (RNAseq), we analyzed transcriptome changes in AT from 191 obese, nondiabetic patients within a multicenter, controlled dietary intervention. Expression changes were associated with outcomes after an 8-wk LCD (800-1000 kcal/d) and 6 mo after the LCD. Results were validated by using quantitative reverse transcriptase-polymerase chain reaction in 350 subjects from the same cohort. Statistical models were constructed to classify weight maintainers or glycemic improvers.Results: With RNAseq analyses, we identified 1173 genes that were differentially expressed after the LCD, of which 350 and 33 were associated with changes in body mass index (BMI; in kg/m²) and Matsuda index values, respectively, whereas 29 genes were associated with both endpoints. Pathway analyses highlighted enrichment in lipid and glucose metabolism. Classification models were constructed to identify weight maintainers. A model based on clinical baseline variables could not achieve any classification (validation AUC: 0.50; 95% CI: 0.36, 0.64). However, clinical changes during the LCD yielded better performance of the model (AUC: 0.73; 95% CI: 0.60, 0.87]). Adding baseline expression to this model improved the performance significantly (AUC: 0.87; 95% CI: 0.77, 0.96; Delong's P = 0.012). Similar analyses were performed to classify subjects with good glycemic improvements. Baseline- and LCD-based clinical models yielded similar performance (best AUC: 0.73; 95% CI: 0.60, 0.86). The addition of expression changes during the LCD improved the performance substantially (AUC: 0.80; 95% CI: 0.69, 0.92; P = 0.058).Conclusions: This study investigated AT transcriptome alterations after an LCD in a large cohort of obese, nondiabetic patients. Gene expression combined with clinical variables enabled us to distinguish weight and glycemic responders from nonresponders. These potential biomarkers may help clinicians understand intersubject variability and better predict the success of dietary interventions. This trial was registered at clinicaltrials.gov as NCT00390637.

MIDAS: Mining differentially activated subpaths of KEGG pathways from multi-class RNA-seq data

Pathway based analysis of high throughput transcriptome data is a widely used approach to investigate biological mechanisms. Since a pathway consists of multiple functions, the recent approach is to determine condition specific sub-pathways or subpaths. However, there are several challenges. First, few existing methods utilize explicit gene expression information from RNA-seq. More importantly, subpath activity is usually an average of statistical scores, e.g., correlations, of edges in a candidate subpath, which fails to reflect gene expression quantity information. In addition, none of existing methods can handle multiple phenotypes. To address these technical problems, we designed and implemented an algorithm, MIDAS, that determines condition specific subpaths, each of which has different activities across multiple phenotypes. MIDAS utilizes gene expression quantity information fully and the network centrality information to determine condition specific subpaths. To test performance of our tool, we used TCGA breast cancer RNA-seq gene expression profiles with five molecular subtypes. 36 differentially activate subpaths were determined. The utility of our method, MIDAS, was demonstrated in four ways. All 36 subpaths are well supported by the literature information. Subsequently, we showed that these subpaths had a good discriminant power for five cancer subtype classification and also had a prognostic power in terms of survival analysis. Finally, in a performance comparison of MIDAS to a recent subpath prediction method, PATHOME, our method identified more subpaths and much more genes that are well supported by the literature information.

AVAILABILITY

http://biohealth.snu.ac.kr/software/MIDAS/.

In silico analysis of pathways activation landscape in oral squamous cell carcinoma and oral leukoplakia

A subset of patients with oral squamous cell carcinoma (OSCC), the most common subtype of head and neck squamous cell carcinoma (HNSCC), harbor dysplastic lesions (often visually identified as leukoplakia) prior to cancer diagnosis. Although evidence suggest that leukoplakia represents an initial step in the progression to cancer, signaling networks driving this progression are poorly understood. Here, we applied in silico Pathway Activation Network Decomposition Analysis (iPANDA), a new bioinformatics software suite for qualitative analysis of intracellular signaling pathway activation using transcriptomic data, to assess a network of molecular signaling in OSCC and pre-neoplastic oral lesions. In tumor samples, our analysis detected major conserved mitogenic and survival signaling pathways strongly associated with HNSCC, suggesting that some of the pathways identified by our algorithm, but not yet validated as HNSCC related, may be attractive targets for future research. While pathways activation landscape in the majority of leukoplakias was different from that seen in OSCC, a subset of pre-neoplastic lesions has demonstrated some degree of similarity to the signaling profile seen in tumors, including dysregulation of the cancer-driving pathways related to survival and apoptosis. These results suggest that dysregulation of these signaling networks may be the driving force behind the early stages of OSCC tumorigenesis. While future studies with larger leukoplakia data sets are warranted to further estimate the values of this approach for capturing signaling features that characterize relevant lesions that actually progress to cancers, our platform proposes a promising new approach for detecting cancer-promoting pathways and tailoring the right therapy to prevent tumorigenesis.

Identifying biologically relevant putative mechanisms in a given phenotype comparison

A major challenge in life science research is understanding the mechanism involved in a given phenotype. The ability to identify the correct mechanisms is needed in order to understand fundamental and very important phenomena such as mechanisms of disease, immune systems responses to various challenges, and mechanisms of drug action. The current data analysis methods focus on the identification of the differentially expressed (DE) genes using their fold change and/or p-values. Major shortcomings of this approach are that: i) it does not consider the interactions between genes; ii) its results are sensitive to the selection of the threshold(s) used, and iii) the set of genes produced by this approach is not always conducive to formulating mechanistic hypotheses. Here we present a method that can construct networks of genes that can be considered putative mechanisms. The putative mechanisms constructed by this approach are not limited to the set of DE genes, but also considers all known and relevant gene-gene interactions. We analyzed three real datasets for which both the causes of the phenotype, as well as the true mechanisms were known. We show that the method identified the correct mechanisms when applied on microarray datasets from mouse. We compared the results of our method with the results of the classical approach, showing that our method produces more meaningful biological insights.

Pathways to smoking behaviours: biological insights from the Tobacco and Genetics Consortium meta-analysis

By running gene and pathway analyses for several smoking behaviours in the Tobacco and Genetics Consortium (TAG) sample of 74 053 individuals, 21 genes and several chains of biological pathways were implicated. Analyses were carried out using the HYbrid Set-based Test (HYST) as implemented in the Knowledge-based mining system for Genome-wide Genetic studies software. Fifteen genes are novel and were not detected with the single nucleotide polymorphism-based approach in the original TAG analysis. For quantity smoked, 14 genes passed the false discovery rate of 0.05 (corrected for multiple testing), with the top association signal located at the IREB2 gene (P=1.57E-37). Three genomic loci were significantly associated with ever smoked. The top signal is located at the noncoding antisense RNA transcript BDNF-AS (P=6.25E-07) on 11p14. The SLC25A21 gene (P=2.09E-08) yielded the top association signal in the analysis of smoking cessation. The 19q13 noncoding RNA locus exceeded the genome-wide significance in the analysis of age at initiation (P=1.33E-06). Pathways belonging to the Neuronal system pathways, harbouring the nicotinic acetylcholine receptor genes expressing the α (CHRNA 1-9), β (CHRNB 1-4), γ, δ and ɛ subunits, yielded the smallest P-values in the pathway analysis of the quantity smoked (lowest P=4.90E-42). Additionally, pathways belonging to 'a subway map of cancer pathways' regulating the cell cycle, mitotic DNA replication, axon growth and synaptic plasticity were found significantly enriched for genetic variants in ever smokers relative to never smokers (lowest P=1.61E-07). In addition, these pathways were also significantly associated with the quantity smoked (lowest P=4.28E-17). Our results shed light on one of the world's leading causes of preventable death and open a path to potential therapeutic targets. These results are informative in decoding the biological bases of other disease traits, such as depression and cancers, with which smoking shares genetic vulnerabilities.

Functional Analysis of OMICs Data and Small Molecule Compounds in an Integrated "Knowledge-Based" Platform

Analysis of NGS and other sequencing data, gene variants, gene expression, proteomics, and other high-throughput (OMICs) data is challenging because of its biological complexity and high level of technical and biological noise. One way to deal with both problems is to perform analysis with a high fidelity annotated knowledgebase of protein interactions, pathways, and functional ontologies. This knowledgebase has to be structured in a computer-readable format and must include software tools for managing experimental data, analysis, and reporting. Here, we present MetaCore™ and Key Pathway Advisor (KPA), an integrated platform for functional data analysis. On the content side, MetaCore and KPA encompass a comprehensive database of molecular interactions of different types, pathways, network models, and ten functional ontologies covering human, mouse, and rat genes. The analytical toolkit includes tools for gene/protein list enrichment analysis, statistical "interactome" tool for the identification of over- and under-connected proteins in the dataset, and a biological network analysis module made up of network generation algorithms and filters. The suite also features Advanced Search, an application for combinatorial search of the database content, as well as a Java-based tool called Pathway Map Creator for drawing and editing custom pathway maps. Applications of MetaCore and KPA include molecular mode of action of disease research, identification of potential biomarkers and drug targets, pathway hypothesis generation, analysis of biological effects for novel small molecule compounds and clinical applications (analysis of large cohorts of patients, and translational and personalized medicine).

Prioritizing biological pathways by recognizing context in time-series gene expression data

Background

The primary goal of pathway analysis using transcriptome data is to find significantly perturbed pathways. However, pathway analysis is not always successful in identifying pathways that are truly relevant to the context under study. A major reason for this difficulty is that a single gene is involved in multiple pathways. In the KEGG pathway database, there are 146 genes, each of which is involved in more than 20 pathways. Thus activation of even a single gene will result in activation of many pathways. This complex relationship often makes the pathway analysis very difficult. While we need much more powerful pathway analysis methods, a readily available alternative way is to incorporate the literature information.

Results

In this study, we propose a novel approach for prioritizing pathways by combining results from both pathway analysis tools and literature information. The basic idea is as follows. Whenever there are enough articles that provide evidence on which pathways are relevant to the context, we can be assured that the pathways are indeed related to the context, which is termed as relevance in this paper. However, if there are few or no articles reported, then we should rely on the results from the pathway analysis tools, which is termed as significance in this paper. We realized this concept as an algorithm by introducing Context Score and Impact Score and then combining the two into a single score. Our method ranked truly relevant pathways significantly higher than existing pathway analysis tools in experiments with two data sets.

Conclusions

Our novel framework was implemented as ContextTRAP by utilizing two existing tools, TRAP and BEST. ContextTRAP will be a useful tool for the pathway based analysis of gene expression data since the user can specify the context of the biological experiment in a set of keywords. The web version of ContextTRAP is available at http://biohealth.snu.ac.kr/software/contextTRAP.

Electronic supplementary material

The online version of this article (doi:10.1186/s12859-016-1335-8) contains supplementary material, which is available to authorized users.

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 joint effect 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.

MinePath: Mining for Phenotype Differential Sub-paths in Molecular Pathways

Pathway analysis methodologies couple traditional gene expression analysis with knowledge encoded in established molecular pathway networks, offering a promising approach towards the biological interpretation of phenotype differentiating genes. Early pathway analysis methodologies, named as gene set analysis (GSA), view pathways just as plain lists of genes without taking into account either the underlying pathway network topology or the involved gene regulatory relations. These approaches, even if they achieve computational efficiency and simplicity, consider pathways that involve the same genes as equivalent in terms of their gene enrichment characteristics. Most recent pathway analysis approaches take into account the underlying gene regulatory relations by examining their consistency with gene expression profiles and computing a score for each profile. Even with this approach, assessing and scoring single-relations limits the ability to reveal key gene regulation mechanisms hidden in longer pathway sub-paths. We introduce MinePath, a pathway analysis methodology that addresses and overcomes the aforementioned problems. MinePath facilitates the decomposition of pathways into their constituent sub-paths. Decomposition leads to the transformation of single-relations to complex regulation sub-paths. Regulation sub-paths are then matched with gene expression sample profiles in order to evaluate their functional status and to assess phenotype differential power. Assessment of differential power supports the identification of the most discriminant profiles. In addition, MinePath assess the significance of the pathways as a whole, ranking them by their p-values. Comparison results with state-of-the-art pathway analysis systems are indicative for the soundness and reliability of the MinePath approach. In contrast with many pathway analysis tools, MinePath is a web-based system (www.minepath.org) offering dynamic and rich pathway visualization functionality, with the unique characteristic to color regulatory relations between genes and reveal their phenotype inclination. This unique characteristic makes MinePath a valuable tool for in silico molecular biology experimentation as it serves the biomedical researchers' exploratory needs to reveal and interpret the regulatory mechanisms that underlie and putatively govern the expression of target phenotypes.

Pathway Analysis Incorporating Protein-Protein Interaction Networks Identified Candidate Pathways for the Seven Common Diseases

Pathway analysis has become popular as a secondary analysis strategy for genome-wide association studies (GWAS). Most of the current pathway analysis methods aggregate signals from the main effects of single nucleotide polymorphisms (SNPs) in genes within a pathway without considering the effects of gene-gene interactions. However, gene-gene interactions can also have critical effects on complex diseases. Protein-protein interaction (PPI) networks have been used to define gene pairs for the gene-gene interaction tests. Incorporating the PPI information to define gene pairs for interaction tests within pathways can increase the power for pathway-based association tests. We propose a pathway association test, which aggregates the interaction signals in PPI networks within a pathway, for GWAS with case-control samples. Gene size is properly considered in the test so that genes do not contribute more to the test statistic simply due to their size. Simulation studies were performed to verify that the method is a valid test and can have more power than other pathway association tests in the presence of gene-gene interactions within a pathway under different scenarios. We applied the test to the Wellcome Trust Case Control Consortium GWAS datasets for seven common diseases. The most significant pathway is the chaperones modulate interferon signaling pathway for Crohn’s disease (p-value = 0.0003). The pathway modulates interferon gamma, which induces the JAK/STAT pathway that is involved in Crohn’s disease. Several other pathways that have functional implications for the seven diseases were also identified. The proposed test based on gene-gene interaction signals in PPI networks can be used as a complementary tool to the current existing pathway analysis methods focusing on main effects of genes. An efficient software implementing the method is freely available at http://puppi.sourceforge.net.

Proteomic Screening and Lasso Regression Reveal Differential Signaling in Insulin and Insulin-like Growth Factor I (IGF1) Pathways

Insulin and insulin-like growth factor I (IGF1) influence cancer risk and progression through poorly understood mechanisms. To better understand the roles of insulin and IGF1 signaling in breast cancer, we combined proteomic screening with computational network inference to uncover differences in IGF1 and insulin induced signaling. Using reverse phase protein array, we measured the levels of 134 proteins in 21 breast cancer cell lines stimulated with IGF1 or insulin for up to 48 h. We then constructed directed protein expression networks using three separate methods: (i) lasso regression, (ii) conventional matrix inversion, and (iii) entropy maximization. These networks, named here as the time translation models, were analyzed and the inferred interactions were ranked by differential magnitude to identify pathway differences. The two top candidates, chosen for experimental validation, were shown to regulate IGF1/insulin induced phosphorylation events. First, acetyl-CoA carboxylase (ACC) knock-down was shown to increase the level of mitogen-activated protein kinase (MAPK) phosphorylation. Second, stable knock-down of E-Cadherin increased the phospho-Akt protein levels. Both of the knock-down perturbations incurred phosphorylation responses stronger in IGF1 stimulated cells compared with insulin. Overall, the time-translation modeling coupled to wet-lab experiments has proven to be powerful in inferring differential interactions downstream of IGF1 and insulin signaling, in vitro.

A systems biology-led insight into the role of the proteome in neurodegenerative diseases

INTRODUCTION

Multifactorial disorders are the result of nonlinear interactions of several factors; therefore, a reductionist approach does not appear to be appropriate. Proteomics is a global approach that can be efficiently used to investigate pathogenetic mechanisms of neurodegenerative diseases.

AREAS COVERED

Here, we report a general introduction about the systems biology approach and mechanistic insights recently obtained by over-representation analysis of proteomics data of cellular and animal models of Alzheimer's disease, Parkinson's disease and other neurodegenerative disorders, as well as of affected human tissues. Expert commentary: As an inductive method, proteomics is based on unbiased observations that further require validation of generated hypotheses. Pathway databases and over-representation analysis tools allow researchers to assign an expectation value to pathogenetic mechanisms linked to neurodegenerative diseases. The systems biology approach based on omics data may be the key to unravel the complex mechanisms underlying neurodegeneration.

Protein function in precision medicine: deep understanding with machine learning

Precision medicine and personalized health efforts propose leveraging complex molecular, medical and family history, along with other types of personal data toward better life. We argue that this ambitious objective will require advanced and specialized machine learning solutions. Simply skimming some low-hanging results off the data wealth might have limited potential. Instead, we need to better understand all parts of the system to define medically relevant causes and effects: how do particular sequence variants affect particular proteins and pathways? How do these effects, in turn, cause the health or disease-related phenotype? Toward this end, deeper understanding will not simply diffuse from deeper machine learning, but from more explicit focus on understanding protein function, context-specific protein interaction networks, and impact of variation on both.

Common pathway signature in lung and liver fibrosis

Fibrosis, a progressive accumulation of extracellular matrix components, encompasses a wide spectrum of distinct organs, and accounts for an increasing burden of morbidity and mortality worldwide. Despite the tremendous clinical impact, the mechanisms governing the fibrotic process are not yet understood, and to date, no clinically reliable therapies for fibrosis have been discovered. Here we applied Regeneration Intelligence, a new bioinformatics software suite for qualitative analysis of intracellular signaling pathway activation using transcriptomic data, to assess a network of molecular signaling in lung and liver fibrosis. In both tissues, our analysis detected major conserved signaling pathways strongly associated with fibrosis, suggesting that some of the pathways identified by our algorithm but not yet wet-lab validated as fibrogenesis related, may be attractive targets for future research. While the majority of significantly disrupted pathways were specific to histologically distinct organs, several pathways have been concurrently activated or downregulated among the hepatic and pulmonary fibrosis samples, providing new evidence of evolutionary conserved pathways that may be relevant as possible therapeutic targets. While future confirmatory studies are warranted to validate these observations, our platform proposes a promising new approach for detecting fibrosis-promoting pathways and tailoring the right therapy to prevent fibrogenesis.

Genome-wide association study of antidepressant response: involvement of the inorganic cation transmembrane transporter activity pathway

BACKGROUND

Genome-wide association studies (GWAS) represent the current frontier in pharmacogenomics. Thousands of subjects of Caucasian ancestry have been included in previous GWAS investigating antidepressant response. GWAS focused on this phenotype are lacking in Asian populations.

METHODS

A sample of 109 major depressive disorder (MDD) patients of Korean origin in antidepressant treatment was collected. Phenotypes were response and remission according to the Hamilton Rating Scale for Depression (HRSD). Genome-wide genotyping was performed using the Illumina Human Omni2.5-8 platform. The same phenotypes were used in the STAR*D level 1 (n = 1677) for independent replication. In order to corroborate findings and increase the comparability between the two datasets, three levels of analysis (SNPs, genes and pathways) were carried out. Bonferroni correction, permutations, and replication across samples were used to reduce the risk of false positives.

RESULTS

Among the genes replicated across the two samples (permutated p < 0.05 in both of them), CTNNA3 appeared promising. The inorganic cation transmembrane transporter activity pathway (GO:0022890) was associated with antidepressant response in both samples (p = 2.9e-5 and p = 0.001 in the Korean and STAR*D samples, respectively) and this pathway included CACNA1A, CACNA1C, and CACNB2 genes.

CONCLUSIONS

The present study supported the involvement of genes coding for subunits of L-type voltage-gated calcium channel in antidepressant efficacy across different ethnicities but replication of findings is required before any definitive statement.

Genetic variants involved in oxidative stress, base excision repair, DNA methylation, and folate metabolism pathways influence myeloid neoplasias susceptibility and prognosis

Myelodysplastic syndrome (MDS) and acute myeloid leukemia (AML) share common features: elevated oxidative stress, DNA repair deficiency, and aberrant DNA methylation. We performed a hospital-based case-control study to evaluate the association in variants of genes involved in oxidative stress, folate metabolism, DNA repair, and DNA methylation with susceptibility and prognosis of these malignancies. To that end, 16 SNPs (one per gene: CAT, CYBA, DNMT1, DNMT3A, DNMT3B, GPX1, KEAP1, MPO, MTRR, NEIL1, NFE2F2, OGG1, SLC19A1, SOD1, SOD2, and XRCC1) were genotyped in 191 patients (101 MDS and 90 AML) and 261 controls. We also measured oxidative stress (reactive oxygen species/total antioxidant status ratio), DNA damage (8-hydroxy-2'-deoxyguanosine), and DNA methylation (5-methylcytosine) in 50 subjects (40 MDS and 10 controls). Results showed that five genes (GPX1, NEIL1, NFE2L2, OGG1, and SOD2) were associated with MDS, two (DNMT3B and SLC19A1) with AML, and two (CYBA and DNMT1) with both diseases. We observed a correlation of CYBA TT, GPX1 TT, and SOD2 CC genotypes with increased oxidative stress levels, as well as NEIL1 TT and OGG1 GG genotypes with higher DNA damage. The 5-methylcytosine levels were negatively associated with DNMT1 CC, DNMT3A CC, and MTRR AA genotypes, and positively with DNMT3B CC genotype. Furthermore, DNMT3A, MTRR, NEIL1, and OGG1 variants modulated AML transformation in MDS patients. Additionally, DNMT3A, OGG1, GPX1, and KEAP1 variants influenced survival of MDS and AML patients. Altogether, data suggest that genetic variability influence predisposition and prognosis of MDS and AML patients, as well AML transformation rate in MDS patients. © 2016 Wiley Periodicals, Inc.

Dynamics in Transcriptomics: Advancements in RNA-seq Time Course and Downstream Analysis

Analysis of gene expression has contributed to a plethora of biological and medical research studies. Microarrays have been intensively used for the profiling of gene expression during diverse developmental processes, treatments and diseases. New massively parallel sequencing methods, often named as RNA-sequencing (RNA-seq) are extensively improving our understanding of gene regulation and signaling networks. Computational methods developed originally for microarrays analysis can now be optimized and applied to genome-wide studies in order to have access to a better comprehension of the whole transcriptome. This review addresses current challenges on RNA-seq analysis and specifically focuses on new bioinformatics tools developed for time series experiments. Furthermore, possible improvements in analysis, data integration as well as future applications of differential expression analysis are discussed.

Pharmacogenetics of major depressive disorder: top genes and pathways toward clinical applications

The pharmacogenetics of antidepressants has been not only a challenging but also frustrating research field since its birth in the 1990s. Indeed, great expectations followed the first evidence of familiar aggregation of antidepressant response. Despite the progress from candidate gene studies to genome-wide association studies (GWAS), results fell out the expectations and they were often inconsistent. Anyway, the cumulative evidence supports the involvement of some genes and molecular pathways in antidepressant efficacy. The best single genes are SLC6A4, HTR2A, BDNF, GNB3, FKBP5, ABCB1, and cytochrome P450 genes (CYP2D6 and CYP2C19). Molecular pathways involved in inflammation and neuroplasticity show the greatest support. The first studies evaluating benefits of genotype-guided antidepressant treatments provided encouraging results and confirmed the relevance of SLC6A4, HTR2A, ABCB1, and cytochrome P450 genes. Further progress in genotyping and data analysis would allow to move forward and complete the understanding of antidepressant pharmacogenetics and its translation into clinical applications.

Systems biology of neurodegenerative diseases

Neurodegenerative diseases (NDs) collectively afflict more than 40 million people worldwide. The majority of these diseases lack therapies to slow or stop progression due in large part to the challenge of disentangling the simultaneous presentation of broad, multifaceted pathophysiologic changes. Present technologies and computational capabilities suggest an optimistic future for deconvolving these changes to identify novel mechanisms driving ND onset and progression. In particular, integration of highly multi-dimensional omic analytical techniques (e.g., microarray, mass spectrometry) with computational systems biology approaches provides a systematic methodology to elucidate new mechanisms driving NDs. In this review, we begin by summarizing the complex pathophysiology of NDs associated with protein aggregation, emphasizing the shared complex dysregulation found in all of these diseases, and discuss available experimental ND models. Next, we provide an overview of technological and computational techniques used in systems biology that are applicable to studying NDs. We conclude by reviewing prior studies that have applied these approaches to NDs and comment on the necessity of combining analysis from both human tissues and model systems to identify driving mechanisms. We envision that the integration of computational approaches with multiple omic analyses of human tissues, and mouse and in vitro models, will enable the discovery of new therapeutic strategies for these devastating diseases.