Systems analysis of eleven rodent disease models reveals an inflammatome signature and key drivers

Context-enriched interactome powered by proteomics helps the identification of novel regulators of macrophage activation

The role of pro-inflammatory macrophage activation in cardiovascular disease (CVD) is a complex one amenable to network approaches. While an indispensible tool for elucidating the molecular underpinnings of complex diseases including CVD, the interactome is limited in its utility as it is not specific to any cell type, experimental condition or disease state. We introduced context-specificity to the interactome by combining it with co-abundance networks derived from unbiased proteomics measurements from activated macrophage-like cells. Each macrophage phenotype contributed to certain regions of the interactome. Using a network proximity-based prioritization method on the combined network, we predicted potential regulators of macrophage activation. Prediction performance significantly increased with the addition of co-abundance edges, and the prioritized candidates captured inflammation, immunity and CVD signatures. Integrating the novel network topology with transcriptomics and proteomics revealed top candidate drivers of inflammation. In vitro loss-of-function experiments demonstrated the regulatory role of these proteins in pro-inflammatory signaling.

When human cells or tissues are injured, the body triggers a response known as inflammation to repair the damage and protect itself from further harm. However, if the same issue keeps recurring, the tissues become inflamed for longer periods of time, which may ultimately lead to health problems. This is what could be happening in cardiovascular diseases, where long-term inflammation could damage the heart and blood vessels.

Many different proteins interact with each other to control inflammation; gaining an insight into the nature of these interactions could help to pinpoint the role of each molecular actor. Researchers have used a combination of unbiased, large-scale experimental and computational approaches to develop the interactome, a map of the known interactions between all proteins in humans. However, interactions between proteins can change between cell types, or during disease. Here, Halu et al. aimed to refine the human interactome and identify new proteins involved in inflammation, especially in the context of cardiovascular disease.

Cells called macrophages produce signals that trigger inflammation whey they detect damage in other cells or tissues. The experiments used a technique called proteomics to measure the amounts of all the proteins in human macrophages. Combining these data with the human interactome made it possible to predict new links between proteins known to have a role in inflammation and other proteins in the interactome. Further analysis using other sets of data from macrophages helped identify two new candidate proteins – GBP1 and WARS – that may promote inflammation. Halu et al. then used a genetic approach to deactivate the genes and decrease the levels of these two proteins in macrophages, which caused the signals that encourage inflammation to drop.

These findings suggest that GBP1 and WARS regulate the activity of macrophages to promote inflammation. The two proteins could therefore be used as drug targets to treat cardiovascular diseases and other disorders linked to inflammation, but further studies will be needed to precisely dissect how GBP1 and WARS work in humans.

Integration of Molecular Interactome and Targeted Interaction Analysis to Identify a COPD Disease Network Module

The polygenic nature of complex diseases offers potential opportunities to utilize network-based approaches that leverage the comprehensive set of protein-protein interactions (the human interactome) to identify new genes of interest and relevant biological pathways. However, the incompleteness of the current human interactome prevents it from reaching its full potential to extract network-based knowledge from gene discovery efforts, such as genome-wide association studies, for complex diseases like chronic obstructive pulmonary disease (COPD). Here, we provide a framework that integrates the existing human interactome information with experimental protein-protein interaction data for FAM13A, one of the most highly associated genetic loci to COPD, to find a more comprehensive disease network module. We identified an initial disease network neighborhood by applying a random-walk method. Next, we developed a network-based closeness approach (C_AB) that revealed 9 out of 96 FAM13A interacting partners identified by affinity purification assays were significantly close to the initial network neighborhood. Moreover, compared to a similar method (local radiality), the C_AB approach predicts low-degree genes as potential candidates. The candidates identified by the network-based closeness approach were combined with the initial network neighborhood to build a comprehensive disease network module (163 genes) that was enriched with genes differentially expressed between controls and COPD subjects in alveolar macrophages, lung tissue, sputum, blood, and bronchial brushing datasets. Overall, we demonstrate an approach to find disease-related network components using new laboratory data to overcome incompleteness of the current interactome.

GM-CSF driven myeloid cells in adipose tissue link weight gain and insulin resistance via formation of 2-aminoadipate

In a GM-CSF driven myeloid cell deficient mouse model (Csf2^−/−) that has preserved insulin sensitivity despite increased adiposity, we used unbiased three-dimensional integration of proteome profiles, metabolic profiles, and gene regulatory networks to understand adipose tissue proteome-wide changes and their metabolic implications. Multi-dimensional liquid chromatography mass spectrometry and extended multiplex mass labeling was used to analyze proteomes of epididymal adipose tissues isolated from Csf2^+/+ and Csf2^−/− mice that were fed low fat, high fat, or high fat plus cholesterol diets for 8 weeks. The metabolic health (as measured by body weight, adiposity, plasma fasting glucose, insulin, triglycerides, phospholipids, total cholesterol levels, and glucose and insulin tolerance tests) deteriorated with diet for both genotypes, while mice lacking Csf2 were protected from insulin resistance. Regardless of diet, 30 mostly mitochondrial, branch chain amino acids (BCAA), and lysine metabolism proteins were altered between Csf2^−/− and Csf2^+/+ mice (FDR < 0.05). Lack of GM-CSF driven myeloid cells lead to reduced adipose tissue 2-oxoglutarate dehydrogenase complex (DHTKD1) levels and subsequent increase in plasma 2-aminoadipate (2-AA) levels, both of which are reported to correlate with insulin resistance. Tissue DHTKD1 levels were >4-fold upregulated and plasma 2-AA levels were >2 fold reduced in Csf2^−/− mice (p < 0.05). GM-CSF driven myeloid cells link peripheral insulin sensitivity to adiposity via lysine metabolism involving DHTKD1/2-AA axis in a diet independent manner.

Using Gene Expression to Annotate Cardiovascular GWAS Loci

Genetic variants at hundreds of loci associated with cardiovascular phenotypes have been identified by genome wide association studies. Most of these variants are located in intronic or intergenic regions rendering the functional and mechanistic follow up difficult. These non-protein-coding regions harbor regulatory sequences. Thus the study of genetic variants associated with transcription—so called expression quantitative trait loci—has emerged as a promising approach to identify regulatory sequence variants. The genes and pathways they control constitute candidate causal drivers at cardiovascular risk loci. This review provides an overview of the expression quantitative trait loci resources available for cardiovascular genetics research and the most commonly used approaches for candidate gene identification.

Integrative gene network analysis identifies key signatures, intrinsic networks and host factors for influenza virus A infections

Influenza A virus, with the limited coding capacity of 10–14 proteins, requires the host cellular machinery for many aspects of its life cycle. Knowledge of these host cell requirements not only reveals molecular pathways exploited by the virus or triggered by the immune system, but also provides further targets for antiviral drug development. To uncover novel pathways and key targets of influenza infection, we assembled a large amount of data from 12 cell-based gene-expression studies of influenza infection for an integrative network analysis. We systematically identified differentially expressed genes and gene co-expression networks induced by influenza infection. We revealed the dedicator of cytokinesis 5 (DOCK5) played potentially an important role for influenza virus replication. CRISPR/Cas9 knockout of DOCK5 reduced influenza virus replication, indicating that DOCK5 is a key regulator for the viral life cycle. DOCK5’s targets determined by the DOCK5 knockout experiments strongly validated the predicted gene signatures and networks. This study systematically uncovered and validated fundamental patterns of molecular responses, intrinsic structures of gene co-regulation, and novel key targets in influenza virus infection.

Molecular response to influenza infection involves a large number of interacting pathways in the form of complex molecular networks. Most studies on influenza infection have largely focused on testing specific molecules and hypotheses with limited data. Therefore, a global picture of molecular interactions in influenza infection is missing. In this study, we performed an integrative network analysis on a large amount of data from 12 cell-based gene expression studies of influenza infections. By combining differential expression, co-expression networks, and gene knockout experiments, we uncovered and validated fundamental patterns of molecular responses, intrinsic structures of gene co-regulation, and novel key targets in influenza infection. Our findings pave the way for other functional investigations into identifying novel therapeutic targets against influenza infection.

Gene network analysis: from heart development to cardiac therapy

Networks offer a flexible framework to represent and analyse the complex interactions between components of cellular systems. In particular gene networks inferred from expression data can support the identification of novel hypotheses on regulatory processes. In this review we focus on the use of gene network analysis in the study of heart development. Understanding heart development will promote the elucidation of the aetiology of congenital heart disease and thus possibly improve diagnostics. Moreover, it will help to establish cardiac therapies. For example, understanding cardiac differentiation during development will help to guide stem cell differentiation required for cardiac tissue engineering or to enhance endogenous repair mechanisms. We introduce different methodological frameworks to infer networks from expression data such as Boolean and Bayesian networks. Then we present currently available temporal expression data in heart development and discuss the use of network-based approaches in published studies. Collectively, our literature-based analysis indicates that gene network analysis constitutes a promising opportunity to infer therapy-relevant regulatory processes in heart development. However, the use of network-based approaches has so far been limited by the small amount of samples in available datasets. Thus, we propose to acquire high-resolution temporal expression data to improve the mathematical descriptions of regulatory processes obtained with gene network inference methodologies. Especially probabilistic methods that accommodate the intrinsic variability of biological systems have the potential to contribute to a deeper understanding of heart development.

Multiscale network modeling of oligodendrocytes reveals molecular components of myelin dysregulation in Alzheimer's disease

BACKGROUND

Oligodendrocytes (OLs) and myelin are critical for normal brain function and have been implicated in neurodegeneration. Several lines of evidence including neuroimaging and neuropathological data suggest that Alzheimer's disease (AD) may be associated with dysmyelination and a breakdown of OL-axon communication.

METHODS

In order to understand this phenomenon on a molecular level, we systematically interrogated OL-enriched gene networks constructed from large-scale genomic, transcriptomic and proteomic data obtained from human AD postmortem brain samples. We then validated these networks using gene expression datasets generated from mice with ablation of major gene expression nodes identified in our AD-dysregulated networks.

RESULTS

The robust OL gene coexpression networks that we identified were highly enriched for genes associated with AD risk variants, such as BIN1 and demonstrated strong dysregulation in AD. We further corroborated the structure of the corresponding gene causal networks using datasets generated from the brain of mice with ablation of key network drivers, such as UGT8, CNP and PLP1, which were identified from human AD brain data. Further, we found that mice with genetic ablations of Cnp mimicked aspects of myelin and mitochondrial gene expression dysregulation seen in brain samples from patients with AD, including decreased protein expression of BIN1 and GOT2.

CONCLUSIONS

This study provides a molecular blueprint of the dysregulation of gene expression networks of OL in AD and identifies key OL- and myelination-related genes and networks that are highly associated with AD.

A functional genomics predictive network model identifies regulators of inflammatory bowel disease

A major challenge in inflammatory bowel disease (IBD) is the integration of diverse IBD data sets to construct predictive models of IBD. We present a predictive model of the immune component of IBD that informs causal relationships among loci previously linked to IBD through genome-wide association studies (GWAS) using functional and regulatory annotations that relate to the cells, tissues, and pathophysiology of IBD. Our model consists of individual networks constructed using molecular data generated from intestinal samples isolated from three populations of patients with IBD at different stages of disease. We performed key driver analysis to identify genes predicted to modulate network regulatory states associated with IBD, prioritizing and prospectively validating 12 of the top key drivers experimentally. This validated key driver set not only introduces new regulators of processes central to IBD but also provides the integrated circuits of genetic, molecular, and clinical traits that can be directly queried to interrogate and refine the regulatory framework defining IBD.

Shared genetic regulatory networks for cardiovascular disease and type 2 diabetes in multiple populations of diverse ethnicities in the United States

Cardiovascular diseases (CVD) and type 2 diabetes (T2D) are closely interrelated complex diseases likely sharing overlapping pathogenesis driven by aberrant activities in gene networks. However, the molecular circuitries underlying the pathogenic commonalities remain poorly understood. We sought to identify the shared gene networks and their key intervening drivers for both CVD and T2D by conducting a comprehensive integrative analysis driven by five multi-ethnic genome-wide association studies (GWAS) for CVD and T2D, expression quantitative trait loci (eQTLs), ENCODE, and tissue-specific gene network models (both co-expression and graphical models) from CVD and T2D relevant tissues. We identified pathways regulating the metabolism of lipids, glucose, and branched-chain amino acids, along with those governing oxidation, extracellular matrix, immune response, and neuronal system as shared pathogenic processes for both diseases. Further, we uncovered 15 key drivers including HMGCR, CAV1, IGF1 and PCOLCE, whose network neighbors collectively account for approximately 35% of known GWAS hits for CVD and 22% for T2D. Finally, we cross-validated the regulatory role of the top key drivers using in vitro siRNA knockdown, in vivo gene knockout, and two Hybrid Mouse Diversity Panels each comprised of >100 strains. Findings from this in-depth assessment of genetic and functional data from multiple human cohorts provide strong support that common sets of tissue-specific molecular networks drive the pathogenesis of both CVD and T2D across ethnicities and help prioritize new therapeutic avenues for both CVD and T2D.

Molecular and genetic inflammation networks in major human diseases

It has been well-recognized that inflammation alongside tissue repair and damage maintaining tissue homeostasis determines the initiation and progression of complex diseases. Albeit with the accomplishment of having captured the most critical inflammation-involved molecules, genetic susceptibilities, epigenetic factors, and environmental factors, our schemata on the role of inflammation in complex diseases remain largely patchy, in part due to the success of reductionism in terms of research methodology per se. Omics data alongside the advances in data integration technologies have enabled reconstruction of molecular and genetic inflammation networks which shed light on the underlying pathophysiology of complex diseases or clinical conditions. Given the proven beneficial role of anti-inflammation in coronary heart disease as well as other complex diseases and immunotherapy as a revolutionary transition in oncology, it becomes timely to review our current understanding of the molecular and genetic inflammation networks underlying major human diseases. In this review, we first briefly discuss the complexity of infectious diseases and then highlight recently uncovered molecular and genetic inflammation networks in other major human diseases including obesity, type II diabetes, coronary heart disease, late onset Alzheimer's disease, Parkinson's disease, and sporadic cancer. The commonality and specificity of these molecular networks are addressed in the context of genetics based on genome-wide association study (GWAS). The double-sword role of inflammation, such as how the aberrant type 1 and/or type 2 immunity leads to chronic and severe clinical conditions, remains open in terms of the inflammasome and the core inflammatome network features. Increasingly available large Omics and clinical data in tandem with systems biology approaches have offered an exciting yet challenging opportunity toward reconstruction of more comprehensive and dynamic molecular and genetic inflammation networks, which hold great promise in transiting network snapshots to video-style multi-scale interplays of disease mechanisms, in turn leading to effective clinical intervention.

Cross-species gene modules emerge from a systems biology approach to osteoarthritis

Complexities in degenerative disorders, such as osteoarthritis, arise from multiscale biological, environmental, and temporal perturbations. Animal models serve to provide controlled representations of the natural history of degenerative disorders, but in themselves represent an additional layer of complexity. Comparing transcriptomic networks arising from gene co-expression data across species can facilitate an understanding of the preservation of functional gene modules and establish associations with disease phenotypes. This study demonstrates the preservation of osteoarthritis-associated gene modules, described by immune system and system development processes, across human and rat studies. Class prediction analysis establishes a minimal gene signature, including the expression of the Rho GDP dissociation inhibitor ARHGDIB, which consistently defined healthy human cartilage from osteoarthritic cartilage in an independent data set. The age of human clinical samples remains a strong confounder in defining the underlying gene regulatory mechanisms in osteoarthritis; however, defining preserved gene models across species may facilitate standardization of animal models of osteoarthritis to better represent human disease and control for ageing phenomena.

Cortical Spreading Depression Causes Unique Dysregulation of Inflammatory Pathways in a Transgenic Mouse Model of Migraine

Familial hemiplegic migraine type 1 (FHM1) is a rare monogenic subtype of migraine with aura caused by mutations in CACNA1A that encodes the α1A subunit of voltage-gated CaV2.1 calcium channels. Transgenic knock-in mice that carry the human FHM1 R192Q missense mutation ('FHM1 R192Q mice') exhibit an increased susceptibility to cortical spreading depression (CSD), the mechanism underlying migraine aura. Here, we analysed gene expression profiles from isolated cortical tissue of FHM1 R192Q mice 24 h after experimentally induced CSD in order to identify molecular pathways affected by CSD. Gene expression profiles were generated using deep serial analysis of gene expression sequencing. Our data reveal a signature of inflammatory signalling upon CSD in the cortex of both mutant and wild-type mice. However, only in the brains of FHM1 R192Q mice specific genes are up-regulated in response to CSD that are implicated in interferon-related inflammatory signalling. Our findings show that CSD modulates inflammatory processes in both wild-type and mutant brains, but that an additional unique inflammatory signature becomes expressed after CSD in a relevant mouse model of migraine.

Extended Multiplexing of Tandem Mass Tags (TMT) Labeling Reveals Age and High Fat Diet Specific Proteome Changes in Mouse Epididymal Adipose Tissue

The lack of high-throughput methods to analyze the adipose tissue protein composition limits our understanding of the protein networks responsible for age and diet related metabolic response. We have developed an approach using multiple-dimension liquid chromatography tandem mass spectrometry and extended multiplexing (24 biological samples) with tandem mass tags (TMT) labeling to analyze proteomes of epididymal adipose tissues isolated from mice fed either low or high fat diet for a short or a long-term, and from mice that aged on low versus high fat diets. The peripheral metabolic health (as measured by body weight, adiposity, plasma fasting glucose, insulin, triglycerides, total cholesterol levels, and glucose and insulin tolerance tests) deteriorated with diet and advancing age, with long-term high fat diet exposure being the worst. In response to short-term high fat diet, 43 proteins representing lipid metabolism (e.g. AACS, ACOX1, ACLY) and red-ox pathways (e.g. CPD2, CYP2E, SOD3) were significantly altered (FDR < 10%). Long-term high fat diet significantly altered 55 proteins associated with immune response (e.g. IGTB2, IFIT3, LGALS1) and rennin angiotensin system (e.g. ENPEP, CMA1, CPA3, ANPEP). Age-related changes on low fat diet significantly altered only 18 proteins representing mainly urea cycle (e.g. OTC, ARG1, CPS1), and amino acid biosynthesis (e.g. GMT, AKR1C6). Surprisingly, high fat diet driven age-related changes culminated with alterations in 155 proteins involving primarily the urea cycle (e.g. ARG1, CPS1), immune response/complement activation (e.g. C3, C4b, C8, C9, CFB, CFH, FGA), extracellular remodeling (e.g. EFEMP1, FBN1, FBN2, LTBP4, FERMT2, ECM1, EMILIN2, ITIH3) and apoptosis (e.g. YAP1, HIP1, NDRG1, PRKCD, MUL1) pathways. Using our adipose tissue tailored approach we have identified both age-related and high fat diet specific proteomic signatures highlighting a pronounced involvement of arginine metabolism in response to advancing age, and branched chain amino acid metabolism in early response to high fat feeding. Data are available via ProteomeXchange with identifier PXD005953.

Analyses of caspase-1-regulated transcriptomes in various tissues lead to identification of novel IL-1β-, IL-18- and sirtuin-1-independent pathways

BACKGROUND

It is well established that caspase-1 exerts its biological activities through its downstream targets such as IL-1β, IL-18, and Sirt-1. The microarray datasets derived from various caspase-1 knockout tissues indicated that caspase-1 can significantly impact the transcriptome. However, it is not known whether all the effects exerted by caspase-1 on transcriptome are mediated only by its well-known substrates. Therefore, we hypothesized that the effects of caspase-1 on transcriptome may be partially independent from IL-1β, IL-18, and Sirt-1.

METHODS

To determine new global and tissue-specific gene regulatory effects of caspase-1, we took novel microarray data analysis approaches including Venn analysis, cooperation analysis, and meta-analysis methods. We used these statistical methods to integrate different microarray datasets conducted on different caspase-1 knockout tissues and datasets where caspase-1 downstream targets were manipulated.

RESULTS

We made the following important findings: (1) Caspase-1 exerts its regulatory effects on the majority of genes in a tissue-specific manner; (2) Caspase-1 regulatory genes partially cooperates with genes regulated by sirtuin-1 during organ injury and inflammation in adipose tissue but not in the liver; (3) Caspase-1 cooperates with IL-1β in regulating less than half of the genes involved in cardiovascular disease, organismal injury, and cancer in mouse liver; (4) The meta-analysis identifies 40 caspase-1 globally regulated genes across tissues, suggesting that caspase-1 globally regulates many novel pathways; and (5) The meta-analysis identified new cooperatively and non-cooperatively regulated genes in caspase-1, IL-1β, IL-18, and Sirt-1 pathways.

CONCLUSIONS

Our findings suggest that caspase-1 regulates many new signaling pathways potentially via its known substrates and also via transcription factors and other proteins that are yet to be identified.

Traumatic Brain Injury Induces Genome-Wide Transcriptomic, Methylomic, and Network Perturbations in Brain and Blood Predicting Neurological Disorders

The complexity of the traumatic brain injury (TBI) pathology, particularly concussive injury, is a serious obstacle for diagnosis, treatment, and long-term prognosis. Here we utilize modern systems biology in a rodent model of concussive injury to gain a thorough view of the impact of TBI on fundamental aspects of gene regulation, which have the potential to drive or alter the course of the TBI pathology. TBI perturbed epigenomic programming, transcriptional activities (expression level and alternative splicing), and the organization of genes in networks centered around genes such as Anax2, Ogn, and Fmod. Transcriptomic signatures in the hippocampus are involved in neuronal signaling, metabolism, inflammation, and blood function, and they overlap with those in leukocytes from peripheral blood. The homology between genomic signatures from blood and brain elicited by TBI provides proof of concept information for development of biomarkers of TBI based on composite genomic patterns. By intersecting with human genome-wide association studies, many TBI signature genes and network regulators identified in our rodent model were causally associated with brain disorders with relevant link to TBI. The overall results show that concussive brain injury reprograms genes which could lead to predisposition to neurological and psychiatric disorders, and that genomic information from peripheral leukocytes has the potential to predict TBI pathogenesis in the brain.

EXPLORING THE REPRODUCIBILITY OF PROBABILISTIC CAUSAL MOLECULAR NETWORK MODELS

Network reconstruction algorithms are increasingly being employed in biomedical and life sciences research to integrate large-scale, high-dimensional data informing on living systems. One particular class of probabilistic causal networks being applied to model the complexity and causal structure of biological data is Bayesian networks (BNs). BNs provide an elegant mathematical framework for not only inferring causal relationships among many different molecular and higher order phenotypes, but also for incorporating highly diverse priors that provide an efficient path for incorporating existing knowledge. While significant methodological developments have broadly enabled the application of BNs to generate and validate meaningful biological hypotheses, the reproducibility of BNs in this context has not been systematically explored. In this study, we aim to determine the criteria for generating reproducible BNs in the context of transcription-based regulatory networks. We utilize two unique tissues from independent datasets, whole blood from the GTEx Consortium and liver from the Stockholm-Tartu Atherosclerosis Reverse Network Engineering Team (STARNET) study. We evaluated the reproducibility of the BNs by creating networks on data subsampled at different levels from each cohort and comparing these networks to the BNs constructed using the complete data. To help validate our results, we used simulated networks at varying sample sizes. Our study indicates that reproducibility of BNs in biological research is an issue worthy of further consideration, especially in light of the many publications that now employ findings from such constructs without appropriate attention paid to reproducibility. We find that while edge-to-edge reproducibility is strongly dependent on sample size, identification of more highly connected key driver nodes in BNs can be carried out with high confidence across a range of sample sizes.

Mergeomics: multidimensional data integration to identify pathogenic perturbations to biological systems

BACKGROUND

Complex diseases are characterized by multiple subtle perturbations to biological processes. New omics platforms can detect these perturbations, but translating the diverse molecular and statistical information into testable mechanistic hypotheses is challenging. Therefore, we set out to create a public tool that integrates these data across multiple datasets, platforms, study designs and species in order to detect the most promising targets for further mechanistic studies.

RESULTS

We developed Mergeomics, a computational pipeline consisting of independent modules that 1) leverage multi-omics association data to identify biological processes that are perturbed in disease, and 2) overlay the disease-associated processes onto molecular interaction networks to pinpoint hubs as potential key regulators. Unlike existing tools that are mostly dedicated to specific data type or settings, the Mergeomics pipeline accepts and integrates datasets across platforms, data types and species. We optimized and evaluated the performance of Mergeomics using simulation and multiple independent datasets, and benchmarked the results against alternative methods. We also demonstrate the versatility of Mergeomics in two case studies that include genome-wide, epigenome-wide and transcriptome-wide datasets from human and mouse studies of total cholesterol and fasting glucose. In both cases, the Mergeomics pipeline provided statistical and contextual evidence to prioritize further investigations in the wet lab. The software implementation of Mergeomics is freely available as a Bioconductor R package.

CONCLUSION

Mergeomics is a flexible and robust computational pipeline for multidimensional data integration. It outperforms existing tools, and is easily applicable to datasets from different studies, species and omics data types for the study of complex traits.

The relationships among host transcriptional responses reveal distinct signatures underlying viral infection-disease associations

Genome-scale DNA microarrays and computational biology facilitate new understanding of viral infections at the system level. Recent years have witnessed a major shift from microorganism-centric toward host-oriented characterization and categorization of viral infections and infection related diseases. We established host transcriptional response (HTR) relationships among 23 different types of human viral pathogens based on calculating HTR similarities using computational integration of 587 public available gene expression profiles. We further identified five virus clusters that show consensus internal HTRs and defined cluster signatures using common dysregulated genes. Individual cluster signature genes were distinguished from one another, and functional analysis revealed common and specific host cellular bioprocesses and signaling pathways involved in confronting viral infections. Through literature investigation and support from epidemiological studies, these were confirmed to be important gene factors associating viral infections with cluster-common and cluster-specific non-infectious human disease(s). Our analyses were the first to feature differential HTRs to viral infections as clusters, and they present a new perspective for understanding infection-disease associations and the underlying pathogeneses.

Genetics and Genomics of Coronary Artery Disease

Coronary artery disease (or coronary heart disease), is the leading cause of mortality in many of the developing as well as the developed countries of the world. Cholesterol-enriched plaques in the heart's blood vessels combined with inflammation lead to the lesion expansion, narrowing of blood vessels, reduced blood flow, and may subsequently cause lesion rupture and a heart attack. Even though several environmental risk factors have been established, such as high LDL-cholesterol, diabetes, and high blood pressure, the underlying genetic composition may substantially modify the disease risk; hence, genome composition and gene-environment interactions may be critical for disease progression. Ongoing scientific efforts have seen substantial advancements related to the fields of genetics and genomics, with the major breakthroughs yet to come. As genomics is the most rapidly advancing field in the life sciences, it is important to present a comprehensive overview of current efforts. Here, we present a summary of various genetic and genomics assays and approaches applied to coronary artery disease research.

Cardiometabolic risk loci share downstream cis- and trans-gene regulation across tissues and diseases

Genome-wide association studies (GWAS) have identified hundreds of cardiometabolic disease (CMD) risk loci. However, they contribute little to genetic variance, and most downstream gene-regulatory mechanisms are unknown. We genotyped and RNA-sequenced vascular and metabolic tissues from 600 coronary artery disease patients in the Stockholm-Tartu Atherosclerosis Reverse Networks Engineering Task study (STARNET). Gene expression traits associated with CMD risk single-nucleotide polymorphism (SNPs) identified by GWAS were more extensively found in STARNET than in tissue- and disease-unspecific gene-tissue expression studies, indicating sharing of downstream cis-/trans-gene regulation across tissues and CMDs. In contrast, the regulatory effects of other GWAS risk SNPs were tissue-specific; abdominal fat emerged as an important gene-regulatory site for blood lipids, such as for the low-density lipoprotein cholesterol and coronary artery disease risk gene PCSK9 STARNET provides insights into gene-regulatory mechanisms for CMD risk loci, facilitating their translation into opportunities for diagnosis, therapy, and prevention.

Reconstructing directed gene regulatory network by only gene expression data

Background

Accurately identifying gene regulatory network is an important task in understanding in vivo biological activities. The inference of such networks is often accomplished through the use of gene expression data. Many methods have been developed to evaluate gene expression dependencies between transcription factor and its target genes, and some methods also eliminate transitive interactions. The regulatory (or edge) direction is undetermined if the target gene is also a transcription factor. Some methods predict the regulatory directions in the gene regulatory networks by locating the eQTL single nucleotide polymorphism, or by observing the gene expression changes when knocking out/down the candidate transcript factors; regrettably, these additional data are usually unavailable, especially for the samples deriving from human tissues.

Results

In this study, we propose the Context Based Dependency Network (CBDN), a method that is able to infer gene regulatory networks with the regulatory directions from gene expression data only. To determine the regulatory direction, CBDN computes the influence of source to target by evaluating the magnitude changes of expression dependencies between the target gene and the others with conditioning on the source gene. CBDN extends the data processing inequality by involving the dependency direction to distinguish between direct and transitive relationship between genes. We also define two types of important regulators which can influence a majority of the genes in the network directly or indirectly. CBDN can detect both of these two types of important regulators by averaging the influence functions of candidate regulator to the other genes. In our experiments with simulated and real data, even with the regulatory direction taken into account, CBDN outperforms the state-of-the-art approaches for inferring gene regulatory network. CBDN identifies the important regulators in the predicted network: 1. TYROBP influences a batch of genes that are related to Alzheimer’s disease; 2. ZNF329 and RB1 significantly regulate those ‘mesenchymal’ gene expression signature genes for brain tumors.

Conclusion

By merely leveraging gene expression data, CBDN can efficiently infer the existence of gene-gene interactions as well as their regulatory directions. The constructed networks are helpful in the identification of important regulators for complex diseases.

Systems Genetic Analyses Highlight a TGFβ-FOXO3 Dependent Striatal Astrocyte Network Conserved across Species and Associated with Stress, Sleep, and Huntington's Disease

Recent systems-based analyses have demonstrated that sleep and stress traits emerge from shared genetic and transcriptional networks, and clinical work has elucidated the emergence of sleep dysfunction and stress susceptibility as early symptoms of Huntington's disease. Understanding the biological bases of these early non-motor symptoms may reveal therapeutic targets that prevent disease onset or slow disease progression, but the molecular mechanisms underlying this complex clinical presentation remain largely unknown. In the present work, we specifically examine the relationship between these psychiatric traits and Huntington's disease (HD) by identifying striatal transcriptional networks shared by HD, stress, and sleep phenotypes. First, we utilize a systems-based approach to examine a large publicly available human transcriptomic dataset for HD (GSE3790 from GEO) in a novel way. We use weighted gene coexpression network analysis and differential connectivity analyses to identify transcriptional networks dysregulated in HD, and we use an unbiased ranking scheme that leverages both gene- and network-level information to identify a novel astrocyte-specific network as most relevant to HD caudate. We validate this result in an independent HD cohort. Next, we computationally predict FOXO3 as a regulator of this network, and use multiple publicly available in vitro and in vivo experimental datasets to validate that this astrocyte HD network is downstream of a signaling pathway important in adult neurogenesis (TGFβ-FOXO3). We also map this HD-relevant caudate subnetwork to striatal transcriptional networks in a large (n = 100) chronically stressed (B6xA/J)F2 mouse population that has been extensively phenotyped (328 stress- and sleep-related measurements), and we show that this striatal astrocyte network is correlated to sleep and stress traits, many of which are known to be altered in HD cohorts. We identify causal regulators of this network through Bayesian network analysis, and we highlight their relevance to motor, mood, and sleep traits through multiple in silico approaches, including an examination of their protein binding partners. Finally, we show that these causal regulators may be therapeutically viable for HD because their downstream network was partially modulated by deep brain stimulation of the subthalamic nucleus, a medical intervention thought to confer some therapeutic benefit to HD patients. In conclusion, we show that an astrocyte transcriptional network is primarily associated to HD in the caudate and provide evidence for its relationship to molecular mechanisms of neural stem cell homeostasis. Furthermore, we present a unified systems-based framework for identifying gene networks that are associated with complex non-motor traits that manifest in the earliest phases of HD. By analyzing and integrating multiple independent datasets, we identify a point of molecular convergence between sleep, stress, and HD that reflects their phenotypic comorbidity and reveals a molecular pathway involved in HD progression.

Network-Based Identification and Prioritization of Key Regulators of Coronary Artery Disease Loci

OBJECTIVE

Recent genome-wide association studies of coronary artery disease (CAD) have revealed 58 genome-wide significant and 148 suggestive genetic loci. However, the molecular mechanisms through which they contribute to CAD and the clinical implications of these findings remain largely unknown. We aim to retrieve gene subnetworks of the 206 CAD loci and identify and prioritize candidate regulators to better understand the biological mechanisms underlying the genetic associations.

APPROACH AND RESULTS

We devised a new integrative genomics approach that incorporated (1) candidate genes from the top CAD loci, (2) the complete genetic association results from the 1000 genomes-based CAD genome-wide association studies from the Coronary Artery Disease Genome Wide Replication and Meta-Analysis Plus the Coronary Artery Disease consortium, (3) tissue-specific gene regulatory networks that depict the potential relationship and interactions between genes, and (4) tissue-specific gene expression patterns between CAD patients and controls. The networks and top-ranked regulators according to these data-driven criteria were further queried against literature, experimental evidence, and drug information to evaluate their disease relevance and potential as drug targets. Our analysis uncovered several potential novel regulators of CAD such as LUM and STAT3, which possess properties suitable as drug targets. We also revealed molecular relations and potential mechanisms through which the top CAD loci operate. Furthermore, we found that multiple CAD-relevant biological processes such as extracellular matrix, inflammatory and immune pathways, complement and coagulation cascades, and lipid metabolism interact in the CAD networks.

CONCLUSIONS

Our data-driven integrative genomics framework unraveled tissue-specific relations among the candidate genes of the CAD genome-wide association studies loci and prioritized novel network regulatory genes orchestrating biological processes relevant to CAD.

Variants in TRIM22 That Affect NOD2 Signaling Are Associated With Very-Early-Onset Inflammatory Bowel Disease

BACKGROUND & AIMS

Severe forms of inflammatory bowel disease (IBD) that develop in very young children can be caused by variants in a single gene. We performed whole-exome sequence (WES) analysis to identify genetic factors that might cause granulomatous colitis and severe perianal disease, with recurrent bacterial and viral infections, in an infant of consanguineous parents.

METHODS

We performed targeted WES analysis of DNA collected from the patient and her parents. We validated our findings by a similar analysis of DNA from 150 patients with very-early-onset IBD not associated with known genetic factors analyzed in Toronto, Oxford, and Munich. We compared gene expression signatures in inflamed vs noninflamed intestinal and rectal tissues collected from patients with treatment-resistant Crohn's disease who participated in a trial of ustekinumab. We performed functional studies of identified variants in primary cells from patients and cell culture.

RESULTS

We identified a homozygous variant in the tripartite motif containing 22 gene (TRIM22) of the patient, as well as in 2 patients with a disease similar phenotype. Functional studies showed that the variant disrupted the ability of TRIM22 to regulate nucleotide binding oligomerization domain containing 2 (NOD2)-dependent activation of interferon-beta signaling and nuclear factor-κB. Computational studies demonstrated a correlation between the TRIM22-NOD2 network and signaling pathways and genetic factors associated very early onset and adult-onset IBD. TRIM22 is also associated with antiviral and mycobacterial effectors and markers of inflammation, such as fecal calprotectin, C-reactive protein, and Crohn's disease activity index scores.

CONCLUSIONS

In WES and targeted exome sequence analyses of an infant with severe IBD characterized by granulomatous colitis and severe perianal disease, we identified a homozygous variant of TRIM22 that affects the ability of its product to regulate NOD2. Combined computational and functional studies showed that the TRIM22-NOD2 network regulates antiviral and antibacterial signaling pathways that contribute to inflammation. Further study of this network could lead to new disease markers and therapeutic targets for patients with very early and adult-onset IBD.

Assessing Concordance of Drug-Induced Transcriptional Response in Rodent Liver and Cultured Hepatocytes

The effect of drugs, disease and other perturbations on mRNA levels are studied using gene expression microarrays or RNA-seq, with the goal of understanding molecular effects arising from the perturbation. Previous comparisons of reproducibility across laboratories have been limited in scale and focused on a single model. The use of model systems, such as cultured primary cells or cancer cell lines, assumes that mechanistic insights derived from the models would have been observed via in vivo studies. We examined the concordance of compound-induced transcriptional changes using data from several sources: rat liver and rat primary hepatocytes (RPH) from Drug Matrix (DM) and open TG-GATEs (TG), human primary hepatocytes (HPH) from TG, and mouse liver / HepG2 results from the Gene Expression Omnibus (GEO) repository. Gene expression changes for treatments were normalized to controls and analyzed with three methods: 1) gene level for 9071 high expression genes in rat liver, 2) gene set analysis (GSA) using canonical pathways and gene ontology sets, 3) weighted gene co-expression network analysis (WGCNA). Co-expression networks performed better than genes or GSA when comparing treatment effects within rat liver and rat vs. mouse liver. Genes and modules performed similarly at Connectivity Map-style analyses, where success at identifying similar treatments among a collection of reference profiles is the goal. Comparisons between rat liver and RPH, and those between RPH, HPH and HepG2 cells reveal lower concordance for all methods. We observe that the baseline state of untreated cultured cells relative to untreated rat liver shows striking similarity with toxicant-exposed cells in vivo, indicating that gross systems level perturbation in the underlying networks in culture may contribute to the low concordance.

Gene expression studies in model systems are widely used for understanding the mechanism of drugs and other perturbations in biological systems. Other researchers have examined the reproducibility of microarray studies between laboratories, or comparing microarrays and/or RNA sequencing. However, no large scale studies have compared results from protocols which differ in minor details, or results generated in vivo vs. in vitro culture systems thought to serve as useful models. The rat liver is by far the most extensively studied model evaluating effects of drugs and other perturbations, and existing data allowed us to assess the level of concordance between rat liver and rat primary hepatocytes cultured in collagen-coated plates (i.e. “flat” culture) for hundreds of drugs. We found that the mouse liver serves as a better model of the rat liver than do rat primary hepatocytes, even after allowing for differences due to pharmacokinetics. The low concordance observed between rat liver and rat hepatocytes suggests that validating the utility of ‘omics data generated on emerging cell culture approaches (e.g. “organ-on-a-chip”, 3D-printed tissues) using rat cells and comparison to the rat liver may be necessary in order to gain confidence these approaches substantially improve on traditional culture models of human cells.

A Next Generation Multiscale View of Inborn Errors of Metabolism

Inborn errors of metabolism (IEM) are not unlike common diseases. They often present as a spectrum of disease phenotypes that correlates poorly with the severity of the disease-causing mutations. This greatly impacts patient care and reveals fundamental gaps in our knowledge of disease modifying biology. Systems biology approaches that integrate multi-omics data into molecular networks have significantly improved our understanding of complex diseases. Similar approaches to study IEM are rare despite their complex nature. We highlight that existing common disease-derived datasets and networks can be repurposed to generate novel mechanistic insight in IEM and potentially identify candidate modifiers. While understanding disease pathophysiology will advance the IEM field, the ultimate goal should be to understand per individual how their phenotype emerges given their primary mutation on the background of their whole genome, not unlike personalized medicine. We foresee that panomics and network strategies combined with recent experimental innovations will facilitate this.

Characterization of Genetic Networks Associated with Alzheimer's Disease

At the molecular level, the genetics of complex disease such as Alzheimer's disease (AD) manifests itself as series of alterations in the molecular interactions in pathways and networks that define biological processes underlying the pathophysiological states of disease. While large-scale genome-wide association (GWA) studies of late-onset alzheimer's disease (LOAD) have uncovered prominent genomic regions linked to the disease, the cause for the vast majority of LOAD cases still remains unknown. Increasingly available large-scale genomic and genetic data related to LOAD has made it possible to comprehensively uncover the mechanisms causally lined to LOAD in a completely data-driven manner. Here we review the various aspects of systems/network biology approaches and methodology in constructing genetic networks associated with AD from large sampling of postmortem brain tissues. We describe in detail a multiscale network modeling approach (MNMA) that integrates interaction and causal gene networks to analyze large-scale DNA, gene expression and pathophysiological data from multiple post-mortem brain regions of LOAD patients as well non-demented normal controls. MNMA first employs weighted gene co-expression network analysis (WGCNA) to construct multi-tissue networks that simultaneously capture intra-tissue and inter-tissue gene-gene interactions and then quantifies the change in connectivity among highly co-expressed genes in LOAD with respect to the normal state. Co-expressed gene modules are then rank ordered by relevance to pathophysiological traits and enrichment of genes differentially expressed in LOAD. Causal regulatory relationships among the genes in each module are then determined by a Bayesian network inference framework that is used to formally integrate genetic and gene expression information. MNMA has uncovered a massive remodeling of network structures in LOAD and identified novel subnetworks and key regulators that are causally linked to LOAD. In the end, we will outline the challenges in systems/network approaches to LOAD.

Efficient Test and Visualization of Multi-Set Intersections

Identification of sets of objects with shared features is a common operation in all disciplines. Analysis of intersections among multiple sets is fundamental for in-depth understanding of their complex relationships. However, so far no method has been developed to assess statistical significance of intersections among three or more sets. Moreover, the state-of-the-art approaches for visualization of multi-set intersections are not scalable. Here, we first developed a theoretical framework for computing the statistical distributions of multi-set intersections based upon combinatorial theory, and then accordingly designed a procedure to efficiently calculate the exact probabilities of multi-set intersections. We further developed multiple efficient and scalable techniques to visualize multi-set intersections and the corresponding intersection statistics. We implemented both the theoretical framework and the visualization techniques in a unified R software package, SuperExactTest. We demonstrated the utility of SuperExactTest through an intensive simulation study and a comprehensive analysis of seven independently curated cancer gene sets as well as six disease or trait associated gene sets identified by genome-wide association studies. We expect SuperExactTest developed by this study will have a broad range of applications in scientific data analysis in many disciplines.

Increased abundance of translation machinery in stem cell-derived neural progenitor cells from four schizophrenia patients

The genetic and epigenetic factors contributing to risk for schizophrenia (SZ) remain unresolved. Here we demonstrate, for the first time, perturbed global protein translation in human-induced pluripotent stem cell (hiPSC)-derived forebrain neural progenitor cells (NPCs) from four SZ patients relative to six unaffected controls. We report increased total protein levels and protein synthesis, together with two independent sets of quantitative mass spectrometry evidence indicating markedly increased levels of ribosomal and translation initiation and elongation factor proteins, in SZ hiPSC NPCs. We posit that perturbed levels of global protein synthesis in SZ hiPSC NPCs represent a novel post-transcriptional mechanism that might contribute to disease progression.

Identification of autophagy signaling network that contributes to stroke in the ischemic rodent brain via gene expression

Autophagy plays a vital role in cerebral ischemia and may be a potential target for developing novel therapy for stroke. In this study, we constructed an autophagy-related pathway network by analyzing the genes related to autophagy and ischemic stroke, and the risk genes were screened. Two autophagy-related modules were significantly up-regulated and clustered to influence cerebral ischemia. Besides, three key modular genes (NFKB1, RELA, and STAT3) were revealed. With 5-fold cross validation, the ROC curves of NFKB1, RELA, and STAT3 were 0.8256, 0.8462, and 0.8923. They formed a complex module and competitively mediated the activation of autophagy in cerebral ischemia. In conclusion, a module containing NFKB1, RELA, and STAT3 mediates autophagy, serving as a potential biomarker for the diagnosis and therapy of ischemic stroke.

Integrative network modeling approaches to personalized cancer medicine

The ability to collect millions of molecular measurements from patients is a now a reality for clinical medicine. This reality has created the challenge of how to integrate these vast amounts of data into models that accurately predict complex pathophysiology and can translate this complexity into clinically actionable outputs. Integrative informatics and data-driven approaches provide a framework for analyzing large-scale datasets and combining them into multiscale models that can be used to determine the key drivers of disease and identify optimal therapies for treating tumors. In this perspective we discuss how an integrative modeling approach is being used to inform individual treatment decisions, highlighting a recent case report that illustrates the challenges and opportunities for personalized oncology.

Bariatric Surgery Induces Disruption in Inflammatory Signaling Pathways Mediated by Immune Cells in Adipose Tissue: A RNA-Seq Study

Background

Bariatric surgery is associated to improvements in obesity-associated comorbidities thought to be mediated by a decrease of adipose inflammation. However, the molecular mechanisms behind these beneficial effects are poorly understood.

Methodology/Principal Findings

We analyzed RNA-seq expression profiles in adipose tissue from 22 obese women before and 3 months after surgery. Of 15,972 detected genes, 1214 were differentially expressed after surgery at a 5% false discovery rate. Upregulated genes were mostly involved in the basal cellular machinery. Downregulated genes were enriched in metabolic functions of adipose tissue. At baseline, 26 modules of coexpressed genes were identified. The four most stable modules reflected the innate and adaptive immune responses of adipose tissue. A first module reflecting a non-specific signature of innate immune cells, mainly macrophages, was highly conserved after surgery with the exception of DUSP2 and CD300C. A second module reflected the adaptive immune response elicited by T lymphocytes; after surgery, a disconnection was observed between genes involved in T-cell signaling and mediators of the signal transduction such as CXCR1, CXCR2, GPR97, CCR7 and IL7R. A third module reflected neutrophil-mediated inflammation; after surgery, several genes were dissociated from the module, including S100A8, S100A12, CD300E, VNN2, TUBB1 and FAM65B. We also identified a dense network of 19 genes involved in the interferon-signaling pathway which was strongly preserved after surgery, with the exception of DDX60, an antiviral factor involved in RIG-I-mediated interferon signaling. A similar loss of connection was observed in lean mice compared to their obese counterparts.

Conclusions/Significance

These results suggest that improvements of the inflammatory state following surgery might be explained by a disruption of immuno-inflammatory cascades involving a few crucial molecules which could serve as potential therapeutic targets.

Genome-wide significant loci: how important are they? Systems genetics to understand heritability of coronary artery disease and other common complex disorders

Genome-wide association studies (GWAS) have been extensively used to study common complex diseases such as coronary artery disease (CAD), revealing 153 suggestive CAD loci, of which at least 46 have been validated as having genome-wide significance. However, these loci collectively explain <10% of the genetic variance in CAD. Thus, we must address the key question of what factors constitute the remaining 90% of CAD heritability. We review possible limitations of GWAS, and contextually consider some candidate CAD loci identified by this method. Looking ahead, we propose systems genetics as a complementary approach to unlocking the CAD heritability and etiology. Systems genetics builds network models of relevant molecular processes by combining genetic and genomic datasets to ultimately identify key "drivers" of disease. By leveraging systems-based genetic approaches, we can help reveal the full genetic basis of common complex disorders, enabling novel diagnostic and therapeutic opportunities.

A systems approach identifies networks and genes linking sleep and stress: implications for neuropsychiatric disorders

Sleep dysfunction and stress susceptibility are comorbid complex traits that often precede and predispose patients to a variety of neuropsychiatric diseases. Here, we demonstrate multilevel organizations of genetic landscape, candidate genes, and molecular networks associated with 328 stress and sleep traits in a chronically stressed population of 338 (C57BL/6J × A/J) F2 mice. We constructed striatal gene co-expression networks, revealing functionally and cell-type-specific gene co-regulations important for stress and sleep. Using a composite ranking system, we identified network modules most relevant for 15 independent phenotypic categories, highlighting a mitochondria/synaptic module that links sleep and stress. The key network regulators of this module are overrepresented with genes implicated in neuropsychiatric diseases. Our work suggests that the interplay among sleep, stress, and neuropathology emerges from genetic influences on gene expression and their collective organization through complex molecular networks, providing a framework for interrogating the mechanisms underlying sleep, stress susceptibility, and related neuropsychiatric disorders.

Integrative network analysis reveals molecular mechanisms of blood pressure regulation

Genome‐wide association studies (GWAS) have identified numerous loci associated with blood pressure (BP). The molecular mechanisms underlying BP regulation, however, remain unclear. We investigated BP‐associated molecular mechanisms by integrating BP GWAS with whole blood mRNA expression profiles in 3,679 individuals, using network approaches. BP transcriptomic signatures at the single‐gene and the coexpression network module levels were identified. Four coexpression modules were identified as potentially causal based on genetic inference because expression‐related SNPs for their corresponding genes demonstrated enrichment for BP GWAS signals. Genes from the four modules were further projected onto predefined molecular interaction networks, revealing key drivers. Gene subnetworks entailing molecular interactions between key drivers and BP‐related genes were uncovered. As proof‐of‐concept, we validated SH2B3, one of the top key drivers, using Sh2b3^−/− mice. We found that a significant number of genes predicted to be regulated by SH2B3 in gene networks are perturbed in Sh2b3^−/− mice, which demonstrate an exaggerated pressor response to angiotensin II infusion. Our findings may help to identify novel targets for the prevention or treatment of hypertension.

Using SRM-MS to quantify nuclear protein abundance differences between adipose tissue depots of insulin-resistant mice

Insulin resistance (IR) underlies metabolic disease. Visceral, but not subcutaneous, white adipose tissue (WAT) has been linked to the development of IR, potentially due to differences in regulatory protein abundance. Here we investigate how protein levels are changed in IR in different WAT depots by developing a targeted proteomics approach to quantitatively compare the abundance of 42 nuclear proteins in subcutaneous and visceral WAT from a commonly used insulin-resistant mouse model, Lepr(db/db), and from C57BL/6J control mice. The most differentially expressed proteins were important in adipogenesis, as confirmed by siRNA-mediated depletion experiments, suggesting a defect in adipogenesis in visceral, but not subcutaneous, insulin-resistant WAT. Furthermore, differentiation of visceral, but not subcutaneous, insulin-resistant stromal vascular cells (SVCs) was impaired. In an in vitro approach to understand the cause of this impaired differentiation, we compared insulin-resistant visceral SVCs to preadipocyte cell culture models made insulin resistant by different stimuli. The insulin-resistant visceral SVC protein abundance profile correlated most with preadipocyte cell culture cells treated with both palmitate and TNFα. Together, our study introduces a method to simultaneously measure and quantitatively compare nuclear protein expression patterns in primary adipose tissue and adipocyte cell cultures, which we show can reveal relationships between differentiation and disease states of different adipocyte tissue types.

Systems Medicine 2.0: potential benefits of combining electronic health care records with systems science models

Background

The global burden of disease is increasingly dominated by non-communicable diseases.These diseases are less amenable to curative and preventative interventions than communicable disease. This presents a challenge to medical practice and medical research, both of which are experiencing diminishing returns from increasing investment.

Objective

Our aim was to (1) review how medical knowledge is generated, and its limitations, (2) assess the potential for emerging technologies and ideas to improve medical research, and (3) suggest solutions and recommendations to increase medical research efficiency on non-communicable diseases.

Methods

We undertook an unsystematic review of peer-reviewed literature and technology websites.

Results

Our review generated the following conclusions and recommendations. (1) Medical knowledge continues to be generated in a reductionist paradigm. This oversimplifies our models of disease, rendering them ineffective to sufficiently understand the complex nature of non-communicable diseases. (2) Some of these failings may be overcome by adopting a “Systems Medicine” paradigm, where the human body is modeled as a complex adaptive system. That is, a system with multiple components and levels interacting in complex ways, wherein disease emerges from slow changes to the system set-up. Pursuing systems medicine research will require larger datasets. (3) Increased data sharing between researchers, patients, and clinicians could provide this unmet need for data. The recent emergence of electronic health care records (EHR) could potentially facilitate this in real-time and at a global level. (4) Efforts should continue to aggregate anonymous EHR data into large interoperable data silos and release this to researchers. However, international collaboration, data linkage, and obtaining additional information from patients will remain challenging. (5) Efforts should also continue towards “Medicine 2.0”. Patients should be given access to their personal EHR data. Subsequently, online communities can give researchers the opportunity to ask patients for direct access to the patient’s EHR data and request additional study-specific information. However, selection bias towards patients who use Web 2.0 technology may be difficult to overcome.

Conclusions

Systems medicine, when combined with large-scale data sharing, has the potential to raise our understanding of non-communicable diseases, foster personalized medicine, and make substantial progress towards halting, curing, and preventing non-communicable diseases. Large-scale data amalgamation remains a core challenge and needs to be supported. A synthesis of “Medicine 2.0” and “Systems Science” concepts into “Systems Medicine 2.0” could take decades to materialize but holds much promise.

Phenotypic differences in hiPSC NPCs derived from patients with schizophrenia

Consistent with recent reports indicating that neurons differentiated in vitro from human-induced pluripotent stem cells (hiPSCs) are immature relative to those in the human brain, gene expression comparisons of our hiPSC-derived neurons to the Allen BrainSpan Atlas indicate that they most resemble fetal brain tissue. This finding suggests that, rather than modeling the late features of schizophrenia (SZ), hiPSC-based models may be better suited for the study of disease predisposition. We now report that a significant fraction of the gene signature of SZ hiPSC-derived neurons is conserved in SZ hiPSC neural progenitor cells (NPCs). We used two independent discovery-based approaches-microarray gene expression and stable isotope labeling by amino acids in cell culture (SILAC) quantitative proteomic mass spectrometry analyses-to identify cellular phenotypes in SZ hiPSC NPCs from four SZ patients. From our findings that SZ hiPSC NPCs show abnormal gene expression and protein levels related to cytoskeletal remodeling and oxidative stress, we predicted, and subsequently observed, aberrant migration and increased oxidative stress in SZ hiPSC NPCs. These reproducible NPC phenotypes were identified through scalable assays that can be applied to expanded cohorts of SZ patients, making them a potentially valuable tool with which to study the developmental mechanisms contributing to SZ.

Integrative network analysis reveals different pathophysiological mechanisms of insulin resistance among Caucasians and African Americans

Background

African Americans (AA) have more pronounced insulin resistance and higher insulin secretion than European Americans (Caucasians or CA) when matched for age, gender, and body mass index (BMI). We hypothesize that physiological differences (including insulin sensitivity [S_I]) between CAs and AAs can be explained by co-regulated gene networks in tissues involved in glucose homeostasis.

Methods

We performed integrative gene network analyses of transcriptomic data in subcutaneous adipose tissue of 99 CA and 37 AA subjects metabolically characterized as non-diabetic, with a range of S_I and BMI values.

Results

Transcripts negatively correlated with S_I in only the CA or AA subjects were enriched for inflammatory response genes and integrin-signaling genes, respectively. A sub-network (module) with TYROBP as a hub enriched for genes involved in inflammatory response (corrected p = 1.7E-26) was negatively correlated with S_I (r = −0.426, p = 4.95E-04) in CA subjects. S_I was positively correlated with transcript modules enriched for mitochondrial metabolism in both groups. Several S_I-associated co-expressed modules were enriched for genes differentially expressed between groups. Two modules involved in immune response to viral infections and function of adherens junction, are significantly correlated with S_I only in CAs. Five modules involved in drug/intracellular transport and oxidoreductase activity, among other activities, are correlated with S_I only in AAs. Furthermore, we identified driver genes of these race-specific S_I-associated modules.

Conclusions

S_I-associated transcriptional networks that were deranged predominantly in one ethnic group may explain the distinctive physiological features of glucose homeostasis among AA subjects.

Electronic supplementary material

The online version of this article (doi:10.1186/s12920-015-0078-0) contains supplementary material, which is available to authorized users.

A disease module in the interactome explains disease heterogeneity, drug response and captures novel pathways and genes in asthma

Recent advances in genetics have spurred rapid progress towards the systematic identification of genes involved in complex diseases. Still, the detailed understanding of the molecular and physiological mechanisms through which these genes affect disease phenotypes remains a major challenge. Here, we identify the asthma disease module, i.e. the local neighborhood of the interactome whose perturbation is associated with asthma, and validate it for functional and pathophysiological relevance, using both computational and experimental approaches. We find that the asthma disease module is enriched with modest GWAS P-values against the background of random variation, and with differentially expressed genes from normal and asthmatic fibroblast cells treated with an asthma-specific drug. The asthma module also contains immune response mechanisms that are shared with other immune-related disease modules. Further, using diverse omics (genomics, gene-expression, drug response) data, we identify the GAB1 signaling pathway as an important novel modulator in asthma. The wiring diagram of the uncovered asthma module suggests a relatively close link between GAB1 and glucocorticoids (GCs), which we experimentally validate, observing an increase in the level of GAB1 after GC treatment in BEAS-2B bronchial epithelial cells. The siRNA knockdown of GAB1 in the BEAS-2B cell line resulted in a decrease in the NFkB level, suggesting a novel regulatory path of the pro-inflammatory factor NFkB by GAB1 in asthma.

Considerations for Using Genetic and Epigenetic Information in Occupational Health Risk Assessment and Standard Setting

Risk assessment forms the basis for both occupational health decision-making and the development of occupational exposure limits (OELs). Although genetic and epigenetic data have not been widely used in risk assessment and ultimately, standard setting, it is possible to envision such uses. A growing body of literature demonstrates that genetic and epigenetic factors condition biological responses to occupational and environmental hazards or serve as targets of them. This presentation addresses the considerations for using genetic and epigenetic information in risk assessments, provides guidance on using this information within the classic risk assessment paradigm, and describes a framework to organize thinking about such uses. The framework is a 4 × 4 matrix involving the risk assessment functions (hazard identification, dose-response modeling, exposure assessment, and risk characterization) on one axis and inherited and acquired genetic and epigenetic data on the other axis. The cells in the matrix identify how genetic and epigenetic data can be used for each risk assessment function. Generally, genetic and epigenetic data might be used as endpoints in hazard identification, as indicators of exposure, as effect modifiers in exposure assessment and dose-response modeling, as descriptors of mode of action, and to characterize toxicity pathways. Vast amounts of genetic and epigenetic data may be generated by high-throughput technologies. These data can be useful for assessing variability and reducing uncertainty in extrapolations, and they may serve as the foundation upon which identification of biological perturbations would lead to a new paradigm of toxicity pathway-based risk assessments.

White adipose tissue reference network: a knowledge resource for exploring health-relevant relations

Optimal health is maintained by interaction of multiple intrinsic and environmental factors at different levels of complexity-from molecular, to physiological, to social. Understanding and quantification of these interactions will aid design of successful health interventions. We introduce the reference network concept as a platform for multi-level exploration of biological relations relevant for metabolic health, by integration and mining of biological interactions derived from public resources and context-specific experimental data. A White Adipose Tissue Health Reference Network (WATRefNet) was constructed as a resource for discovery and prioritization of mechanism-based biomarkers for white adipose tissue (WAT) health status and the effect of food and drug compounds on WAT health status. The WATRefNet (6,797 nodes and 32,171 edges) is based on (1) experimental data obtained from 10 studies addressing different adiposity states, (2) seven public knowledge bases of molecular interactions, (3) expert's definitions of five physiologically relevant processes key to WAT health, namely WAT expandability, Oxidative capacity, Metabolic state, Oxidative stress and Tissue inflammation, and (4) a collection of relevant biomarkers of these processes identified by BIOCLAIMS ( http://bioclaims.uib.es ). The WATRefNet comprehends multiple layers of biological complexity as it contains various types of nodes and edges that represent different biological levels and interactions. We have validated the reference network by showing overrepresentation with anti-obesity drug targets, pathology-associated genes and differentially expressed genes from an external disease model dataset. The resulting network has been used to extract subnetworks specific to the above-mentioned expert-defined physiological processes. Each of these process-specific signatures represents a mechanistically supported composite biomarker for assessing and quantifying the effect of interventions on a physiological aspect that determines WAT health status. Following this principle, five anti-diabetic drug interventions and one diet intervention were scored for the match of their expression signature to the five biomarker signatures derived from the WATRefNet. This confirmed previous observations of successful intervention by dietary lifestyle and revealed WAT-specific effects of drug interventions. The WATRefNet represents a sustainable knowledge resource for extraction of relevant relationships such as mechanisms of action, nutrient intervention targets and biomarkers and for assessment of health effects for support of health claims made on food products.

Considerations for Using Genetic and Epigenetic Information in Occupational Health Risk Assessment and Standard Setting

Risk assessment forms the basis for both occupational health decision-making and the development of occupational exposure limits (OELs). Although genetic and epigenetic data have not been widely used in risk assessment and ultimately, standard setting, it is possible to envision such uses. A growing body of literature demonstrates that genetic and epigenetic factors condition biological responses to occupational and environmental hazards or serve as targets of them. This presentation addresses the considerations for using genetic and epigenetic information in risk assessments, provides guidance on using this information within the classic risk assessment paradigm, and describes a framework to organize thinking about such uses. The framework is a 4 × 4 matrix involving the risk assessment functions (hazard identification, dose-response modeling, exposure assessment, and risk characterization) on one axis and inherited and acquired genetic and epigenetic data on the other axis. The cells in the matrix identify how genetic and epigenetic data can be used for each risk assessment function. Generally, genetic and epigenetic data might be used as endpoints in hazard identification, as indicators of exposure, as effect modifiers in exposure assessment and dose-response modeling, as descriptors of mode of action, and to characterize toxicity pathways. Vast amounts of genetic and epigenetic data may be generated by high-throughput technologies. These data can be useful for assessing variability and reducing uncertainty in extrapolations, and they may serve as the foundation upon which identification of biological perturbations would lead to a new paradigm of toxicity pathway-based risk assessments.

Evolving toward a human-cell based and multiscale approach to drug discovery for CNS disorders

A disruptive approach to therapeutic discovery and development is required in order to significantly improve the success rate of drug discovery for central nervous system (CNS) disorders. In this review, we first assess the key factors contributing to the frequent clinical failures for novel drugs. Second, we discuss cancer translational research paradigms that addressed key issues in drug discovery and development and have resulted in delivering drugs with significantly improved outcomes for patients. Finally, we discuss two emerging technologies that could improve the success rate of CNS therapies: human induced pluripotent stem cell (hiPSC)-based studies and multiscale biology models. Coincident with advances in cellular technologies that enable the generation of hiPSCs directly from patient blood or skin cells, together with methods to differentiate these hiPSC lines into specific neural cell types relevant to neurological disease, it is also now possible to combine data from large-scale forward genetics and post-mortem global epigenetic and expression studies in order to generate novel predictive models. The application of systems biology approaches to account for the multiscale nature of different data types, from genetic to molecular and cellular to clinical, can lead to new insights into human diseases that are emergent properties of biological networks, not the result of changes to single genes. Such studies have demonstrated the heterogeneity in etiological pathways and the need for studies on model systems that are patient-derived and thereby recapitulate neurological disease pathways with higher fidelity. In the context of two common and presumably representative neurological diseases, the neurodegenerative disease Alzheimer's Disease, and the psychiatric disorder schizophrenia, we propose the need for, and exemplify the impact of, a multiscale biology approach that can integrate panomic, clinical, imaging, and literature data in order to construct predictive disease network models that can (i) elucidate subtypes of syndromic diseases, (ii) provide insights into disease networks and targets and (iii) facilitate a novel drug screening strategy using patient-derived hiPSCs to discover novel therapeutics for CNS disorders.

Shared molecular pathways and gene networks for cardiovascular disease and type 2 diabetes mellitus in women across diverse ethnicities

BACKGROUND

Although cardiovascular disease (CVD) and type 2 diabetes mellitus (T2D) share many common risk factors, potential molecular mechanisms that may also be shared for these 2 disorders remain unknown.

METHODS AND RESULTS

Using an integrative pathway and network analysis, we performed genome-wide association studies in 8155 blacks, 3494 Hispanic American, and 3697 Caucasian American women who participated in the national Women's Health Initiative single-nucleotide polymorphism (SNP) Health Association Resource and the Genomics and Randomized Trials Network. Eight top pathways and gene networks related to cardiomyopathy, calcium signaling, axon guidance, cell adhesion, and extracellular matrix seemed to be commonly shared between CVD and T2D across all 3 ethnic groups. We also identified ethnicity-specific pathways, such as cell cycle (specific for Hispanic American and Caucasian American) and tight junction (CVD and combined CVD and T2D in Hispanic American). In network analysis of gene-gene or protein-protein interactions, we identified key drivers that included COL1A1, COL3A1, and ELN in the shared pathways for both CVD and T2D. These key driver genes were cross-validated in multiple mouse models of diabetes mellitus and atherosclerosis.

CONCLUSIONS

Our integrative analysis of American women of 3 ethnicities identified multiple shared biological pathways and key regulatory genes for the development of CVD and T2D. These prospective findings also support the notion that ethnicity-specific susceptibility genes and process are involved in the pathogenesis of CVD and T2D.

Periostin is temporally expressed as an extracellular matrix component in skeletal muscle regeneration and differentiation

The transcriptional events and pathways responsible for the acquisition of the myogenic phenotype during regeneration and myogenesis have been studied extensively. The modulators that shape the extracellular matrix in health and disease, however, are less understood. Understanding the components and pathways of this remodeling will aid the restoration of the architecture and prevent deterioration under pathological conditions such as fibrosis. Periostin, a matricellular protein associated with remodeling of the extracellular matrix and connective tissue architecture, is emerging in pathological conditions associated with fibrosis in adult life. Periostin also complicates fibrosis in degenerative skeletal muscle conditions such as dystrophies. This study primarily addresses the spatial and temporal involvement of periostin along skeletal muscle regeneration. In the acute skeletal muscle injury model that shows recovery without fibrosis, we show that periostin is rapidly disrupted along with the extensive necrosis and periostin mRNA is transiently upregulated during the myotube maturation. This expression is stringently initiated from the newly regenerating fibers. However, this observation is contrasting to a model that displays extensive fibrosis where upregulation of periostin expression is stable and confined to the fibrotic compartments of endomysial and perimysial space. In vitro myoblast differentiation further supports the claim that upregulation of periostin expression is a function of extracellular matrix remodeling during myofiber differentiation and maturation. We further seek to identify the expression kinetics of various periostin isoforms during the differentiation of rat and mouse myoblasts. Results depict that a singular periostin isoform dominated the rat muscle, contrasting to multiple isoforms in C2C12 myoblast cells. This study shows that periostin, a mediator with deleterious impact on conditions exhibiting fibrosis, is also produced and secreted by myoblasts and regenerating myofibers during architectural remodeling in the course of development and regeneration.

Network signatures link hepatic effects of anti-diabetic interventions with systemic disease parameters

Background

Multifactorial diseases such as type 2 diabetes mellitus (T2DM), are driven by a complex network of interconnected mechanisms that translate to a diverse range of complications at the physiological level. To optimally treat T2DM, pharmacological interventions should, ideally, target key nodes in this network that act as determinants of disease progression.

Results

We set out to discover key nodes in molecular networks based on the hepatic transcriptome dataset from a preclinical study in obese LDLR-/- mice recently published by Radonjic et al. Here, we focus on comparing efficacy of anti-diabetic dietary (DLI) and two drug treatments, namely PPARA agonist fenofibrate and LXR agonist T0901317. By combining knowledge-based and data-driven networks with a random walks based algorithm, we extracted network signatures that link the DLI and two drug interventions to dyslipidemia-related disease parameters.

Conclusions

This study identified specific and prioritized sets of key nodes in hepatic molecular networks underlying T2DM, uncovering pathways that are to be modulated by targeted T2DM drug interventions in order to modulate the complex disease phenotype.

A framework for the next generation of risk science

OBJECTIVES

In 2011, the U.S. Environmental Protection Agency initiated the NexGen project to develop a new paradigm for the next generation of risk science.

METHODS

The NexGen framework was built on three cornerstones: the availability of new data on toxicity pathways made possible by fundamental advances in basic biology and toxicological science, the incorporation of a population health perspective that recognizes that most adverse health outcomes involve multiple determinants, and a renewed focus on new risk assessment methodologies designed to better inform risk management decision making.

RESULTS

The NexGen framework has three phases. Phase I (objectives) focuses on problem formulation and scoping, taking into account the risk context and the range of available risk management decision-making options. Phase II (risk assessment) seeks to identify critical toxicity pathway perturbations using new toxicity testing tools and technologies, and to better characterize risks and uncertainties using advanced risk assessment methodologies. Phase III (risk management) involves the development of evidence-based population health risk management strategies of a regulatory, economic, advisory, community-based, or technological nature, using sound principles of risk management decision making.

CONCLUSIONS

Analysis of a series of case study prototypes indicated that many aspects of the NexGen framework are already beginning to be adopted in practice.

Integrative genomics reveals novel molecular pathways and gene networks for coronary artery disease

The majority of the heritability of coronary artery disease (CAD) remains unexplained, despite recent successes of genome-wide association studies (GWAS) in identifying novel susceptibility loci. Integrating functional genomic data from a variety of sources with a large-scale meta-analysis of CAD GWAS may facilitate the identification of novel biological processes and genes involved in CAD, as well as clarify the causal relationships of established processes. Towards this end, we integrated 14 GWAS from the CARDIoGRAM Consortium and two additional GWAS from the Ottawa Heart Institute (25,491 cases and 66,819 controls) with 1) genetics of gene expression studies of CAD-relevant tissues in humans, 2) metabolic and signaling pathways from public databases, and 3) data-driven, tissue-specific gene networks from a multitude of human and mouse experiments. We not only detected CAD-associated gene networks of lipid metabolism, coagulation, immunity, and additional networks with no clear functional annotation, but also revealed key driver genes for each CAD network based on the topology of the gene regulatory networks. In particular, we found a gene network involved in antigen processing to be strongly associated with CAD. The key driver genes of this network included glyoxalase I (GLO1) and peptidylprolyl isomerase I (PPIL1), which we verified as regulatory by siRNA experiments in human aortic endothelial cells. Our results suggest genetic influences on a diverse set of both known and novel biological processes that contribute to CAD risk. The key driver genes for these networks highlight potential novel targets for further mechanistic studies and therapeutic interventions.

Sudden death due to heart attack ranks among the top causes of death in the world, and family studies have shown that genetics has a substantial effect on heart disease risk. Recent studies suggest that multiple genetic factors each with modest effects are necessary for the development of CAD, but the genes and molecular processes involved remain poorly understood. We conducted an integrative genomics study where we used the information of gene-gene interactions to capture groups of genes that are most likely to increase heart disease risk. We not only confirmed the importance of several known CAD risk processes such as the metabolism and transport of cholesterol, immune response, and blood coagulation, but also revealed many novel processes such as neuroprotection, cell cycle, and proteolysis that were not previously implicated in CAD. In particular, we highlight several genes such as GLO1 with key regulatory roles within these processes not detected by the first wave of genetic analyses. These results highlight the value of integrating population genetic data with diverse resources that functionally annotate the human genome. Such integration facilitates the identification of novel molecular processes involved in the pathogenesis of CAD as well as potential novel targets for the development of efficacious therapeutic interventions.

Unifying immunology with informatics and multiscale biology

The immune system is a highly complex and dynamic system. Historically, the most common scientific and clinical practice has been to evaluate its individual components. This kind of approach cannot always expose the interconnecting pathways that control immune-system responses and does not reveal how the immune system works across multiple biological systems and scales. High-throughput technologies can be used to measure thousands of parameters of the immune system at a genome-wide scale. These system-wide surveys yield massive amounts of quantitative data that provide a means to monitor and probe immune-system function. New integrative analyses can help synthesize and transform these data into valuable biological insight. Here we review some of the computational analysis tools for high-dimensional data and how they can be applied to immunology.