Oncogenic Signaling Pathways in The Cancer Genome Atlas.

Genetic alterations in signaling pathways that control cell-cycle progression, apoptosis, and cell growth are common hallmarks of cancer, but the extent, mechanisms, and co-occurrence of alterations in these pathways differ between individual tumors and tumor types. Using mutations, copy-number changes, mRNA expression, gene fusions and DNA methylation in 9,125 tumors profiled by The Cancer Genome Atlas (TCGA), we analyzed the mechanisms and patterns of somatic alterations in ten canonical pathways: cell cycle, Hippo, Myc, Notch, Nrf2, PI-3-Kinase/Akt, RTK-RAS, TGFß signaling, p53 and ß-catenin/Wnt. We charted the detailed landscape of pathway alterations in 33 cancer types, stratified into 64 subtypes, and identified patterns of co-occurrence and mutual exclusivity. Eighty-nine percent of tumors had at least one driver alteration in these pathways, and 57% percent of tumors had at least one alteration potentially targetable by currently available drugs. Thirty percent of tumors had multiple targetable alterations, indicating opportunities for combination therapy.

Pathogenic Germline Variants in 10,389 Adult Cancers.

We conducted the largest investigation of predisposition variants in cancer to date, discovering 853 pathogenic or likely pathogenic variants in 8% of 10,389 cases from 33 cancer types. Twenty-one genes showed single or cross-cancer associations, including novel associations of SDHA in melanoma and PALB2 in stomach adenocarcinoma. The 659 predisposition variants and 18 additional large deletions in tumor suppressors, including ATM, BRCA1, and NF1, showed low gene expression and frequent (43%) loss of heterozygosity or biallelic two-hit events. We also discovered 33 such variants in oncogenes, including missenses in MET, RET, and PTPN11 associated with high gene expression. We nominated 47 additional predisposition variants from prioritized VUSs supported by multiple evidences involving case-control frequency, loss of heterozygosity, expression effect, and co-localization with mutations and modified residues. Our integrative approach links rare predisposition variants to functional consequences, informing future guidelines of variant classification and germline genetic testing in cancer.

The somatic genomic landscape of glioblastoma.

We describe the landscape of somatic genomic alterations based on multidimensional and comprehensive characterization of more than 500 glioblastoma tumors (GBMs). We identify several novel mutated genes as well as complex rearrangements of signature receptors, including EGFR and PDGFRA. TERT promoter mutations are shown to correlate with elevated mRNA expression, supporting a role in telomerase reactivation. Correlative analyses confirm that the survival advantage of the proneural subtype is conferred by the G-CIMP phenotype, and MGMT DNA methylation may be a predictive biomarker for treatment response only in classical subtype GBM. Integrative analysis of genomic and proteomic profiles challenges the notion of therapeutic inhibition of a pathway as an alternative to inhibition of the target itself. These data will facilitate the discovery of therapeutic and diagnostic target candidates, the validation of research and clinical observations and the generation of unanticipated hypotheses that can advance our molecular understanding of this lethal cancer.

Characterizing the Impact of Communication on Cellular and Collective Behavior Using a Three-Dimensional Multiscale Cellular Model.

Communication between cells enables the coordination that drives structural and functional complexity in biological systems. Both single and multicellular organisms have evolved diverse communication systems for a range of purposes, including synchronization of behavior, division of labor, and spatial organization. Synthetic systems are also increasingly being engineered to utilize cell-cell communication. While research has elucidated the form and function of cell-cell communication in many biological systems, our knowledge is still limited by the confounding effects of other biological phenomena at play and the bias of the evolutionary background. In this work, our goal is to push forward the context-free understanding of what impact cell-cell communication can have on cellular and population behavior to more fully understand the extent to which cell-cell communication systems can be utilized, modified, and engineered. We use an in silico model of 3D multiscale cellular populations, with dynamic intracellular networks interacting via diffusible signals. We focus on two key communication parameters: the effective interaction distance at which cells are able to interact and the receptor activation threshold. We found that cell-cell communication can be divided into six different forms along the parameter axes, three asocial and three social. We also show that cellular behavior, tissue composition, and tissue diversity are all highly sensitive to both the general form and specific parameters of communication even when the cellular network has not been biased towards that behavior.

Genomic and molecular characterization of preterm birth.

Preterm birth (PTB) complications are the leading cause of long-term morbidity and mortality in children. By using whole blood samples, we integrated whole-genome sequencing (WGS), RNA sequencing (RNA-seq), and DNA methylation data for 270 PTB and 521 control families. We analyzed this combined dataset to identify genomic variants associated with PTB and secondary analyses to identify variants associated with very early PTB (VEPTB) as well as other subcategories of disease that may contribute to PTB. We identified differentially expressed genes (DEGs) and methylated genomic loci and performed expression and methylation quantitative trait loci analyses to link genomic variants to these expression and methylation changes. We performed enrichment tests to identify overlaps between new and known PTB candidate gene systems. We identified 160 significant genomic variants associated with PTB-related phenotypes. The most significant variants, DEGs, and differentially methylated loci were associated with VEPTB. Integration of all data types identified a set of 72 candidate biomarker genes for VEPTB, encompassing genes and those previously associated with PTB. Notably, PTB-associated genes RAB31 and RBPJ were identified by all three data types (WGS, RNA-seq, and methylation). Pathways associated with VEPTB include EGFR and prolactin signaling pathways, inflammation- and immunity-related pathways, chemokine signaling, IFN-γ signaling, and Notch1 signaling. Progress in identifying molecular components of a complex disease is aided by integrated analyses of multiple molecular data types and clinical data. With these data, and by stratifying PTB by subphenotype, we have identified associations between VEPTB and the underlying biology.

Function of C/EBPdelta in a regulatory circuit that discriminates between transient and persistent TLR4-induced signals.

The innate immune system is like a double-edged sword: it is absolutely required for host defense against infection, but when uncontrolled, it can trigger a plethora of inflammatory diseases. Here we use systems-biology approaches to predict and confirm the existence of a gene-regulatory network involving dynamic interaction among the transcription factors NF-kappaB, C/EBPdelta and ATF3 that controls inflammatory responses. We mathematically modeled transcriptional regulation of the genes encoding interleukin 6 and C/EBPdelta and experimentally confirmed the prediction that the combination of an initiator (NF-kappaB), an amplifier (C/EBPdelta) and an attenuator (ATF3) forms a regulatory circuit that discriminates between transient and persistent Toll-like receptor 4-induced signals. Our results suggest a mechanism that enables the innate immune system to detect the duration of infection and to respond appropriately.

Literature-based discovery of new candidates for drug repurposing.

Drug development is an expensive and time-consuming process; these could be reduced if the existing resources could be used to identify candidates for drug repurposing. This study sought to do this by text mining a large-scale literature repository to curate repurposed drug lists for different cancers. We devised a pattern-based relationship extraction method to extract disease-gene and gene-drug direct relationships from the literature. These direct relationships are used to infer indirect relationships using the ABC model. A gene-shared ranking method based on drug target similarity was then proposed to prioritize the indirect relationships. Our method of assessing drug target similarity correlated to existing anatomical therapeutic chemical code-based methods with a Pearson correlation coefficient of 0.9311. The indirect relationships ranking method achieved a significant mean average precision score of top 100 most common diseases. We also confirmed the suitability of candidates identified for repurposing as anticancer drugs by conducting a manual review of the literature and the clinical trials. Eventually, for visualization and enrichment of huge amount of repurposed drug information, a chord diagram was demonstrated to rapidly identify two novel indications for further biological evaluations.

Solving the influence maximization problem reveals regulatory organization of the yeast cell cycle.

The Influence Maximization Problem (IMP) aims to discover the set of nodes with the greatest influence on network dynamics. The problem has previously been applied in epidemiology and social network analysis. Here, we demonstrate the application to cell cycle regulatory network analysis for Saccharomyces cerevisiae. Fundamentally, gene regulation is linked to the flow of information. Therefore, our implementation of the IMP was framed as an information theoretic problem using network diffusion. Utilizing more than 26,000 regulatory edges from YeastMine, gene expression dynamics were encoded as edge weights using time lagged transfer entropy, a method for quantifying information transfer between variables. By picking a set of source nodes, a diffusion process covers a portion of the network. The size of the network cover relates to the influence of the source nodes. The set of nodes that maximizes influence is the solution to the IMP. By solving the IMP over different numbers of source nodes, an influence ranking on genes was produced. The influence ranking was compared to other metrics of network centrality. Although the top genes from each centrality ranking contained well-known cell cycle regulators, there was little agreement and no clear winner. However, it was found that influential genes tend to directly regulate or sit upstream of genes ranked by other centrality measures. The influential nodes act as critical sources of information flow, potentially having a large impact on the state of the network. Biological events that affect influential nodes and thereby affect information flow could have a strong effect on network dynamics, potentially leading to disease. Code and data can be found at: https://github.com/gibbsdavidl/miergolf.

Multiscale mutation clustering algorithm identifies pan-cancer mutational clusters associated with pathway-level changes in gene expression.

Cancer researchers have long recognized that somatic mutations are not uniformly distributed within genes. However, most approaches for identifying cancer mutations focus on either the entire-gene or single amino-acid level. We have bridged these two methodologies with a multiscale mutation clustering algorithm that identifies variable length mutation clusters in cancer genes. We ran our algorithm on 539 genes using the combined mutation data in 23 cancer types from The Cancer Genome Atlas (TCGA) and identified 1295 mutation clusters. The resulting mutation clusters cover a wide range of scales and often overlap with many kinds of protein features including structured domains, phosphorylation sites, and known single nucleotide variants. We statistically associated these multiscale clusters with gene expression and drug response data to illuminate the functional and clinical consequences of mutations in our clusters. Interestingly, we find multiple clusters within individual genes that have differential functional associations: these include PTEN, FUBP1, and CDH1. This methodology has potential implications in identifying protein regions for drug targets, understanding the biological underpinnings of cancer, and personalizing cancer treatments. Toward this end, we have made the mutation clusters and the clustering algorithm available to the public. Clusters and pathway associations can be interactively browsed at m2c.systemsbiology.net. The multiscale mutation clustering algorithm is available at https://github.com/IlyaLab/M2C.

MIIP haploinsufficiency induces chromosomal instability and promotes tumour progression in colorectal cancer.

The gene encoding migration and invasion inhibitory protein (MIIP), located on 1p36.22, is a potential tumour suppressor gene in glioma. In this study, we aimed to explore the role and mechanism of action of MIIP in colorectal cancer (CRC). MIIP protein expression gradually decreased along the colorectal adenoma-carcinoma sequence and was negatively correlated with lymph node and distant metastasis in 526 colorectal tissue samples (p < 0.05 for all). Analysis of The Cancer Genome Atlas (TCGA) data showed that decreased MIIP expression was significantly associated with MIIP hemizygous deletion (p = 0.0005), which was detected in 27.7% (52/188) of CRC cases, and associated with lymph node and distant metastasis (p < 0.05 for both). We deleted one copy of the MIIP gene in HCT116 CRC cells using zinc finger nuclease technology and demonstrated that MIIP haploinsufficiency resulted in increased colony formation and cell migration and invasion, which was consistent with the results from siRNA-mediated MIIP knockdown in two CRC cell lines (p < 0.05 for all). Moreover, MIIP haploinsufficiency promoted CRC progression in vivo (p < 0.05). Genomic instability and spectral karyotyping assays demonstrated that MIIP haploinsufficiency induced chromosomal instability (CIN). Besides modulating the downstream proteins of APC/CCdc20 , securin and cyclin B1, MIIP haploinsufficiency inhibited topoisomerase II (Topo II) activity and induced chromosomal missegregation. Therefore, we report that MIIP is a novel potential tumour suppressor gene in CRC. Moreover, we characterized the MIIP gene as a novel CIN suppressor gene, through altering the stability of mitotic checkpoint proteins and disturbing Topo II activity. Copyright © 2016 Pathological Society of Great Britain and Ireland. Published by John Wiley & Sons, Ltd.

Genomically amplified Akt3 activates DNA repair pathway and promotes glioma progression.

Akt is a robust oncogene that plays key roles in the development and progression of many cancers, including glioma. We evaluated the differential propensities of the Akt isoforms toward progression in the well-characterized RCAS/Ntv-a mouse model of PDGFB-driven low grade glioma. A constitutively active myristoylated form of Akt1 did not induce high-grade glioma (HGG). In stark contrast, Akt2 and Akt3 showed strong progression potential with 78% and 97% of tumors diagnosed as HGG, respectively. We further revealed that significant variations in polarity and hydropathy values among the Akt isoforms in both the pleckstrin homology domain (P domain) and regulatory domain (R domain) were critical in mediating glioma progression. Gene expression profiles from representative Akt-derived tumors indicated dominant and distinct roles for Akt3, consisting primarily of DNA repair pathways. TCGA data from human GBM closely reflected the DNA repair function, as Akt3 was significantly correlated with a 76-gene signature DNA repair panel. Consistently, compared with Akt1 and Akt2 overexpression models, Akt3-expressing human GBM cells had enhanced activation of DNA repair proteins, leading to increased DNA repair and subsequent resistance to radiation and temozolomide. Given the wide range of Akt3-amplified cancers, Akt3 may represent a key resistance factor.

Flexibility of Boolean Network Reservoir Computers in Approximating Arbitrary Recursive and Non-Recursive Binary Filters.

Reservoir computers (RCs) are biology-inspired computational frameworks for signal processing that are typically implemented using recurrent neural networks. Recent work has shown that Boolean networks (BN) can also be used as reservoirs. We analyze the performance of BN RCs, measuring their flexibility and identifying the factors that determine the effective approximation of Boolean functions applied in a sliding-window fashion over a binary signal, both non-recursively and recursively. We train and test BN RCs of different sizes, signal connectivity, and in-degree to approximate three-bit, five-bit, and three-bit recursive binary functions, respectively. We analyze how BN RC parameters and function average sensitivity, which is a measure of function smoothness, affect approximation accuracy as well as the spread of accuracies for a single reservoir. We found that approximation accuracy and reservoir flexibility are highly dependent on RC parameters. Overall, our results indicate that not all reservoirs are equally flexible, and RC instantiation and training can be more efficient if this is taken into account. The optimum range of RC parameters opens up an angle of exploration for understanding how biological systems might be tuned to balance system restraints with processing capacity.

Multiscale representation of genomic signals.

Genomic information is encoded on a wide range of distance scales, ranging from tens of bases to megabases. We developed a multiscale framework to analyze and visualize the information content of genomic signals. Different types of signals, such as G+C content or DNA methylation, are characterized by distinct patterns of signal enrichment or depletion across scales spanning several orders of magnitude. These patterns are associated with a variety of genomic annotations. By integrating the information across all scales, we demonstrated improved prediction of gene expression from polymerase II chromatin immunoprecipitation sequencing (ChIP-seq) measurements, and we observed that gene expression differences in colorectal cancer are related to methylation patterns that extend beyond the single-gene scale. Our software is available at https://github.com/tknijnen/msr/.

Combining dependent P-values with an empirical adaptation of Brown's method.

MOTIVATION: Combining P-values from multiple statistical tests is a common exercise in bioinformatics. However, this procedure is non-trivial for dependent P-values. Here, we discuss an empirical adaptation of Brown's method (an extension of Fisher's method) for combining dependent P-values which is appropriate for the large and correlated datasets found in high-throughput biology. RESULTS: We show that the Empirical Brown's method (EBM) outperforms Fisher's method as well as alternative approaches for combining dependent P-values using both noisy simulated data and gene expression data from The Cancer Genome Atlas. AVAILABILITY AND IMPLEMENTATION: The Empirical Brown's method is available in Python, R, and MATLAB and can be obtained from https://github.com/IlyaLab/CombiningDependentPvalues UsingEBM The R code is also available as a Bioconductor package from https://www.bioconductor.org/packages/devel/bioc/html/EmpiricalBrownsMethod.html CONTACT: Theo.Knijnenburg@systemsbiology.org SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Gene-pair expression signatures reveal lineage control.

The distinct cell types of multicellular organisms arise owing to constraints imposed by gene regulatory networks on the collective change of gene expression across the genome, creating self-stabilizing expression states, or attractors. We curated human expression data comprising 166 cell types and 2,602 transcription-regulating genes and developed a data-driven method for identifying putative determinants of cell fate built around the concept of expression reversal of gene pairs, such as those participating in toggle-switch circuits. This approach allows us to organize the cell types into their ontogenic lineage relationships. Our method identifies genes in regulatory circuits that control neuronal fate, pluripotency and blood cell differentiation, and it may be useful for prioritizing candidate factors for direct conversion of cell fate.

Structure-based predictions broadly link transcription factor mutations to gene expression changes in cancers.

Thousands of unique mutations in transcription factors (TFs) arise in cancers, and the functional and biological roles of relatively few of these have been characterized. Here, we used structure-based methods developed specifically for DNA-binding proteins to systematically predict the consequences of mutations in several TFs that are frequently mutated in cancers. The explicit consideration of protein-DNA interactions was crucial to explain the roles and prevalence of mutations in TP53 and RUNX1 in cancers, and resulted in a higher specificity of detection for known p53-regulated genes among genetic associations between TP53 genotypes and genome-wide expression in The Cancer Genome Atlas, compared to existing methods of mutation assessment. Biophysical predictions also indicated that the relative prevalence of TP53 missense mutations in cancer is proportional to their thermodynamic impacts on protein stability and DNA binding, which is consistent with the selection for the loss of p53 transcriptional function in cancers. Structure and thermodynamics-based predictions of the impacts of missense mutations that focus on specific molecular functions may be increasingly useful for the precise and large-scale inference of aberrant molecular phenotypes in cancer and other complex diseases.

Systematic measurement of transcription factor-DNA interactions by targeted mass spectrometry identifies candidate gene regulatory proteins.

Regulation of gene expression involves the orchestrated interaction of a large number of proteins with transcriptional regulatory elements in the context of chromatin. Our understanding of gene regulation is limited by the lack of a protein measurement technology that can systematically detect and quantify the ensemble of proteins associated with the transcriptional regulatory elements of specific genes. Here, we introduce a set of selected reaction monitoring (SRM) assays for the systematic measurement of 464 proteins with known or suspected roles in transcriptional regulation at RNA polymerase II transcribed promoters in Saccharomyces cerevisiae. Measurement of these proteins in nuclear extracts by SRM permitted the reproducible quantification of 42% of the proteins over a wide range of abundances. By deploying the assay to systematically identify DNA binding transcriptional regulators that interact with the environmentally regulated FLO11 promoter in cell extracts, we identified 15 regulators that bound specifically to distinct regions along ∼600 bp of the regulatory sequence. Importantly, the dataset includes a number of regulators that have been shown to either control FLO11 expression or localize to these regulatory regions in vivo. We further validated the utility of the approach by demonstrating that two of the SRM-identified factors, Mot3 and Azf1, are required for proper FLO11 expression. These results demonstrate the utility of SRM-based targeted proteomics to guide the identification of gene-specific transcriptional regulators.

Biocellion: accelerating computer simulation of multicellular biological system models.

MOTIVATION: Biological system behaviors are often the outcome of complex interactions among a large number of cells and their biotic and abiotic environment. Computational biologists attempt to understand, predict and manipulate biological system behavior through mathematical modeling and computer simulation. Discrete agent-based modeling (in combination with high-resolution grids to model the extracellular environment) is a popular approach for building biological system models. However, the computational complexity of this approach forces computational biologists to resort to coarser resolution approaches to simulate large biological systems. High-performance parallel computers have the potential to address the computing challenge, but writing efficient software for parallel computers is difficult and time-consuming. RESULTS: We have developed Biocellion, a high-performance software framework, to solve this computing challenge using parallel computers. To support a wide range of multicellular biological system models, Biocellion asks users to provide their model specifics by filling the function body of pre-defined model routines. Using Biocellion, modelers without parallel computing expertise can efficiently exploit parallel computers with less effort than writing sequential programs from scratch. We simulate cell sorting, microbial patterning and a bacterial system in soil aggregate as case studies. AVAILABILITY AND IMPLEMENTATION: Biocellion runs on x86 compatible systems with the 64 bit Linux operating system and is freely available for academic use. Visit http://biocellion.com for additional information.

MiR-506 suppresses proliferation and induces senescence by directly targeting the CDK4/6-FOXM1 axis in ovarian cancer.

Ovarian carcinoma is the most lethal gynaecological malignancy. Better understanding of the molecular pathogenesis of this disease and effective targeted therapies are needed to improve patient outcomes. MicroRNAs play important roles in cancer progression and have the potential for use as either therapeutic agents or targets. Studies in other cancers have suggested that miR-506 has anti-tumour activity, but its function has yet to be elucidated. We found that deregulation of miR-506 in ovarian carcinoma promotes an aggressive phenotype. Ectopic over-expression of miR-506 in ovarian cancer cells was sufficient to inhibit proliferation and to promote senescence. We also demonstrated that CDK4 and CDK6 are direct targets of miR-506, and that miR-506 can inhibit CDK4/6-FOXM1 signalling, which is activated in the majority of serous ovarian carcinomas. This newly recognized miR-506-CDK4/6-FOXM1 axis provides further insight into the pathogenesis of ovarian carcinoma and identifies a potential novel therapeutic agent.

A multilevel pan-cancer map links gene mutations to cancer hallmarks.

BACKGROUND: A central challenge in cancer research is to create models that bridge the gap between the molecular level on which interventions can be designed and the cellular and tissue levels on which the disease phenotypes are manifested. This study was undertaken to construct such a model from functional annotations and explore its use when integrated with large-scale cancer genomics data. METHODS: We created a map that connects genes to cancer hallmarks via signaling pathways. We projected gene mutation and focal copy number data from various cancer types onto this map. We performed statistical analyses to uncover mutually exclusive and co-occurring oncogenic aberrations within this topology. RESULTS: Our analysis showed that although the genetic fingerprint of tumor types could be very different, there were less variations at the level of hallmarks, consistent with the idea that different genetic alterations have similar functional outcomes. Additionally, we showed how the multilevel map could help to clarify the role of infrequently mutated genes, and we demonstrated that mutually exclusive gene mutations were more prevalent in pathways, whereas many co-occurring gene mutations were associated with hallmark characteristics. CONCLUSIONS: Overlaying this map with gene mutation and focal copy number data from various cancer types makes it possible to investigate the similarities and differences between tumor samples systematically at the levels of not only genes but also pathways and hallmarks.