Christoph Wierling

A Systems Biology Approach to Deciphering the Etiology of Steatosis Employing Patient-Derived Dermal Fibroblasts and iPS Cells

Non-alcoholic fatty liver disease comprises a broad spectrum of disease states ranging from simple steatosis to non-alcoholic steatohepatitis. As a result of increases in the prevalences of obesity, insulin resistance, and hyperlipidemia, the number of people with hepatic steatosis continues to increase. Differences in susceptibility to steatohepatitis and its progression to cirrhosis have been attributed to a complex interplay of genetic and external factors all addressing the intracellular network. Increase in sugar or refined carbohydrate consumption results in an increase of insulin and insulin resistance that can lead to the accumulation of fat in the liver. Here we demonstrate how a multidisciplinary approach encompassing cellular reprogramming, transcriptomics, proteomics, metabolomics, modeling, network reconstruction, and data management can be employed to unveil the mechanisms underlying the progression of steatosis. Proteomics revealed reduced AKT/mTOR signaling in fibroblasts derived from steatosis patients and further establishes that the insulin-resistant phenotype is present not only in insulin-metabolizing central organs, e.g., the liver, but is also manifested in skin fibroblasts. Transcriptome data enabled the generation of a regulatory network based on the transcription factor SREBF1, linked to a metabolic network of glycerolipid, and fatty acid biosynthesis including the downstream transcriptional targets of SREBF1 which include LIPIN1 (LPIN) and low density lipoprotein receptor. Glutathione metabolism was among the pathways enriched in steatosis patients in comparison to healthy controls. By using a model of the glutathione pathway we predict a significant increase in the flux through glutathione synthesis as both gamma-glutamylcysteine synthetase and glutathione synthetase have an increased flux. We anticipate that a larger cohort of patients and matched controls will confirm our preliminary findings presented here.

ConsensusPathDB: toward a more complete picture of cell biology

ConsensusPathDB is a meta-database that integrates different types of functional interactions from heterogeneous interaction data resources. Physical protein interactions, metabolic and signaling reactions and gene regulatory interactions are integrated in a seamless functional association network that simultaneously describes multiple functional aspects of genes, proteins, complexes, metabolites, etc. With 155 432 human, 194 480 yeast and 13 648 mouse complex functional interactions (originating from 18 databases on human and eight databases on yeast and mouse interactions each), ConsensusPathDB currently constitutes the most comprehensive publicly available interaction repository for these species. The Web interface at http://cpdb.molgen.mpg.de offers different ways of utilizing these integrated interaction data, in particular with tools for visualization, analysis and interpretation of high-throughput expression data in the light of functional interactions and biological pathways.

ConsensusPathDB—a database for integrating human functional interaction networks

ConsensusPathDB is a database system for the integration of human functional interactions. Current knowledge of these interactions is dispersed in more than 200 databases, each having a specific focus and data format. ConsensusPathDB currently integrates the content of 12 different interaction databases with heterogeneous foci comprising a total of 26 133 distinct physical entities and 74 289 distinct functional interactions (protein–protein interactions, biochemical reactions, gene regulatory interactions), and covering 1738 pathways. We describe the database schema and the methods used for data integration. Furthermore, we describe the functionality of the ConsensusPathDB web interface, where users can search and visualize interaction networks, upload, modify and expand networks in BioPAX, SBML or PSI-MI format, or carry out over-representation analysis with uploaded identifier lists with respect to substructures derived from the integrated interaction network. The ConsensusPathDB database is available at: http://cpdb.molgen.mpg.de

Primary Differentiation in the Human Blastocyst: Comparative Molecular Portraits of Inner Cell Mass and Trophectoderm Cells

The primary differentiation event during mammalian development occurs at the blastocyst stage and leads to the delineation of the inner cell mass (ICM) and the trophectoderm (TE). We provide the first global mRNA expression data from immunosurgically dissected ICM cells, TE cells, and intact human blastocysts. Using a cDNA microarray composed of 15,529 cDNAs from known and novel genes, we identify marker transcripts specific to the ICM (e.g., OCT4/POU5F1, NANOG, HMGB1, and DPPA5) and TE (e.g., CDX2, ATP1B3, SFN, and IPL), in addition to novel ICM‐ and TE‐specific expressed sequence tags. The expression patterns suggest that the emergence of pluripotent ICM and TE cell lineages from the morula is controlled by metabolic and signaling pathways, which include inter alia, WNT, mitogen‐activated protein kinase, transforming growth factor‐beta, NOTCH, integrin‐mediated cell adhesion, phosphatidylinositol 3‐kinase, and apoptosis. These data enhance our understanding of the first step in human cellular differentiation and, hence, the derivation of both embryonic stem cells and trophoblastic stem cells from these lineages.

Hedgehog‐EGFR cooperation response genes determine the oncogenic phenotype of basal cell carcinoma and tumour‐initiating pancreatic cancer cells

Inhibition of Hedgehog (HH)/GLI signalling in cancer is a promising therapeutic approach. Interactions between HH/GLI and other oncogenic pathways affect the strength and tumourigenicity of HH/GLI. Cooperation of HH/GLI with epidermal growth factor receptor (EGFR) signalling promotes transformation and cancer cell proliferation in vitro. However, the in vivo relevance of HH‐EGFR signal integration and the critical downstream mediators are largely undefined. In this report we show that genetic and pharmacologic inhibition of EGFR signalling reduces tumour growth in mouse models of HH/GLI driven basal cell carcinoma (BCC). We describe HH‐EGFR cooperation response genes including SOX2, SOX9, JUN, CXCR4 and FGF19 that are synergistically activated by HH‐EGFR signal integration and required for in vivo growth of BCC cells and tumour‐initiating pancreatic cancer cells. The data validate EGFR signalling as drug target in HH/GLI driven cancers and shed light on the molecular processes controlled by HH‐EGFR signal cooperation, providing new therapeutic strategies based on combined targeting of HH‐EGFR signalling and selected downstream target genes.

Molecular dissection of colorectal cancer in pre-clinical models identifies biomarkers predicting sensitivity to EGFR inhibitors

Colorectal carcinoma represents a heterogeneous entity, with only a fraction of the tumours responding to available therapies, requiring a better molecular understanding of the disease in precision oncology. To address this challenge, the OncoTrack consortium recruited 106 CRC patients (stages I–IV) and developed a pre-clinical platform generating a compendium of drug sensitivity data totalling >4,000 assays testing 16 clinical drugs on patient-derived in vivo and in vitro models. This large biobank of 106 tumours, 35 organoids and 59 xenografts, with extensive omics data comparing donor tumours and derived models provides a resource for advancing our understanding of CRC. Models recapitulate many of the genetic and transcriptomic features of the donors, but defined less complex molecular sub-groups because of the loss of human stroma. Linking molecular profiles with drug sensitivity patterns identifies novel biomarkers, including a signature outperforming RAS/RAF mutations in predicting sensitivity to the EGFR inhibitor cetuximab.

Synergism between Hedgehog-GLI and EGFR Signaling in Hedgehog-Responsive Human Medulloblastoma Cells Induces Downregulation of Canonical Hedgehog-Target Genes and Stabilized Expression of GLI1

Aberrant activation of Hedgehog (HH) signaling has been identified as a key etiologic factor in many human malignancies. Signal strength, target gene specificity, and oncogenic activity of HH signaling depend profoundly on interactions with other pathways, such as epidermal growth factor receptor-mediated signaling, which has been shown to cooperate with HH/GLI in basal cell carcinoma and pancreatic cancer. Our experimental data demonstrated that the Daoy human medulloblastoma cell line possesses a fully inducible endogenous HH pathway. Treatment of Daoy cells with Sonic HH or Smoothened agonist induced expression of GLI1 protein and simultaneously prevented the processing of GLI3 to its repressor form. To study interactions between HH- and EGF-induced signaling in greater detail, time-resolved measurements were carried out and analyzed at the transcriptomic and proteomic levels. The Daoy cells responded to the HH/EGF co-treatment by downregulating GLI1, PTCH, and HHIP at the transcript level; this was also observed when Amphiregulin (AREG) was used instead of EGF. We identified a novel crosstalk mechanism whereby EGFR signaling silences proteins acting as negative regulators of HH signaling, as AKT- and ERK-signaling independent process. EGFR/HH signaling maintained high GLI1 protein levels which contrasted the GLI1 downregulation on the transcript level. Conversely, a high-level synergism was also observed, due to a strong and significant upregulation of numerous canonical EGF-targets with putative tumor-promoting properties such as MMP7, VEGFA, and IL-8. In conclusion, synergistic effects between EGFR and HH signaling can selectively induce a switch from a canonical HH/GLI profile to a modulated specific target gene profile. This suggests that there are more wide-spread, yet context-dependent interactions, between HH/GLI and growth factor receptor signaling in human malignancies.

High-Throughput miRNA and mRNA Sequencing of Paired Colorectal Normal, Tumor and Metastasis Tissues and Bioinformatic Modeling of miRNA-1 Therapeutic Applications

MiRNAs are discussed as diagnostic and therapeutic molecules. However, effective miRNA drug treatments with miRNAs are, so far, hampered by the complexity of the miRNA networks. To identify potential miRNA drugs in colorectal cancer, we profiled miRNA and mRNA expression in matching normal, tumor and metastasis tissues of eight patients by Illumina sequencing. We validated six miRNAs in a large tissue screen containing 16 additional tumor entities and identified miRNA-1, miRNA-129, miRNA-497 and miRNA-215 as constantly de-regulated within the majority of cancers. Of these, we investigated miRNA-1 as representative in a systems-biology simulation of cellular cancer models implemented in PyBioS and assessed the effects of depletion as well as overexpression in terms of miRNA-1 as a potential treatment option. In this system, miRNA-1 treatment reverted the disease phenotype with different effectiveness among the patients. Scoring the gene expression changes obtained through mRNA-Seq from the same patients we show that the combination of deep sequencing and systems biological modeling can help to identify patient-specific responses to miRNA treatments. We present this data as guideline for future pre-clinical assessments of new and personalized therapeutic options.

The structural impact of cancer-associated missense mutations in oncogenes and tumor suppressors

BackgroundCurrent large-scale cancer sequencing projects have identified large numbers of somatic mutations covering an increasing number of different cancer tissues and patients. However, the characterization of these mutations at the structural and functional level remains a challenge.ResultsWe present results from an analysis of the structural impact of frequent missense cancer mutations using an automated method. We find that inactivation of tumor suppressors in cancer correlates frequently with destabilizing mutations preferably in the core of the protein, while enhanced activity of oncogenes is often linked to specific mutations at functional sites. Furthermore, our results show that this alteration of oncogenic activity is often associated with mutations at ATP or GTP binding sites.ConclusionsWith our findings we can confirm and statistically validate the hypotheses for the gain-of-function and loss-of-function mechanisms of oncogenes and tumor suppressors, respectively. We show that the distinct mutational patterns can potentially be used to pre-classify newly identified cancer-associated genes with yet unknown function.

GeNGe: Systematic Generation of Gene Regulatory Networks

Summary: The analysis of gene regulatory networks (GRNs) is a central goal of bioinformatics highly accelerated by the advent of new experimental techniques, such as RNA interference. A battery of reverse engineering methods has been developed in recent years to reconstruct the underlying GRNs from these and other experimental data. However, the performance of the individual methods is poorly understood and validation of algorithmic performances is still missing to a large extent. To enable such systematic validation, we have developed the web application GeNGe (GEne Network GEnerator), a controlled framework for the automatic generation of GRNs. The theoretical model for a GRN is a non-linear differential equation system. Networks can be user-defined or constructed in a modular way with the option to introduce global and local network perturbations. Resulting data can be used, e.g. as benchmark data for evaluating GRN reconstruction methods or for predicting effects of perturbations as theoretical counterparts of biological experiments. Availability: Available online at http://genge.molgen.mpg.de Contact: hache@molgen.mpg.de Supplementary information: Supplementary data are available at Bioinformatics online.