Inferring the temporal order of cancer gene mutations in individual tumor samples

Chromosome 7 to the rescue: overcoming chromosome 10 loss in gliomas

The co-occurrence of chromosome 10 loss and chromosome 7 gain in gliomas is the most frequent loss-gain co-aneuploidy pair in human cancers, a phenomenon that has been investigated without resolution since the late 1980s. Expanding beyond previous gene-centric studies, we investigate the co-occurrence in a genome-wide manner taking an evolutionary perspective. First, by mining large tumor aneuploidy data, we predict that the more likely order is 10 loss followed by 7 gain. Second, by analyzing extensive genomic and transcriptomic data from both patients and cell lines, we find that this co-occurrence can be explained by functional rescue interactions that are highly enriched on 7, which can possibly compensate for any detrimental consequences arising from the loss of 10. Finally, by analyzing transcriptomic data from normal, non-cancerous, human brain tissues, we provide a plausible reason why this co-occurrence happens preferentially in cancers originating in certain regions of the brain.

A phylogenetic approach to inferring the order in which mutations arise during cancer progression

Although the role of evolutionary process in cancer progression is widely accepted, increasing attention is being given to the evolutionary mechanisms that can lead to differences in clinical outcome. Recent studies suggest that the temporal order in which somatic mutations accumulate during cancer progression is important. Single-cell sequencing (SCS) provides a unique opportunity to examine the effect that the mutation order has on cancer progression and treatment effect. However, the error rates associated with single-cell sequencing are known to be high, which greatly complicates the task. We propose a novel method for inferring the order in which somatic mutations arise within an individual tumor using noisy data from single-cell sequencing. Our method incorporates models at two levels in that the evolutionary process of somatic mutation within the tumor is modeled along with the technical errors that arise from the single-cell sequencing data collection process. Through analyses of simulations across a wide range of realistic scenarios, we show that our method substantially outperforms existing approaches for identifying mutation order. Most importantly, our method provides a unique means to capture and quantify the uncertainty in the inferred mutation order along a given phylogeny. We illustrate our method by analyzing data from colorectal and prostate cancer patients, in which our method strengthens previously reported mutation orders. Our work is an important step towards producing meaningful prediction of mutation order with high accuracy and measuring the uncertainty of predicted mutation order in cancer patients, with the potential to lead to new insights about the evolutionary trajectories of cancer.

Parsimonious Clone Tree Integration in cancer

BACKGROUND

Every tumor is composed of heterogeneous clones, each corresponding to a distinct subpopulation of cells that accumulated different types of somatic mutations, ranging from single-nucleotide variants (SNVs) to copy-number aberrations (CNAs). As the analysis of this intra-tumor heterogeneity has important clinical applications, several computational methods have been introduced to identify clones from DNA sequencing data. However, due to technological and methodological limitations, current analyses are restricted to identifying tumor clones only based on either SNVs or CNAs, preventing a comprehensive characterization of a tumor's clonal composition.

RESULTS

To overcome these challenges, we formulate the identification of clones in terms of both SNVs and CNAs as a integration problem while accounting for uncertainty in the input SNV and CNA proportions. We thus characterize the computational complexity of this problem and we introduce PACTION (PArsimonious Clone Tree integratION), an algorithm that solves the problem using a mixed integer linear programming formulation. On simulated data, we show that tumor clones can be identified reliably, especially when further taking into account the ancestral relationships that can be inferred from the input SNVs and CNAs. On 49 tumor samples from 10 prostate cancer patients, our integration approach provides a higher resolution view of tumor evolution than previous studies.

CONCLUSION

PACTION is an accurate and fast method that reconstructs clonal architecture of cancer tumors by integrating SNV and CNA clones inferred using existing methods.

Eutrophication as a driver of microbial community structure in lake sediments

Lake sediments are globally important carbon sinks. Although the fate of organic carbon in lake sediments depends significantly on microorganisms, only few studies have investigated controls on lake sedimentary microbial communities. Here we investigate the impact of anthropogenic eutrophication, which affects redox chemistry and organic matter (OM) sources in sediments, on microbial communities across five lakes in central Switzerland. Lipid biomarkers and distributions of microbial respiration reactions indicate strong increases in aquatic OM contributions and microbial activity with increasing trophic state. Across all lakes, 16S rRNA genes analyses indicate similar depth-dependent zonations at the phylum- and class-level that follow vertical distributions of OM sources and respiration reactions. Yet, there are notable differences, such as higher abundances of nitrifying Bacteria and Archaea in an oligotrophic lake. Furthermore, analyses at the order-level and below suggest that changes in OM sources due to eutrophication cause permanent changes in bacterial community structure. By contrast, archaeal communities are differentiated according to trophic state in recently deposited layers, but converge in older sediments deposited under different trophic regimes. Our study indicates an important role for trophic state in driving lacustrine sediment microbial communities and reveals fundamental differences in the temporal responses of sediment Bacteria and Archaea to eutrophication.

The role of HSP90 molecular chaperones in hepatocellular carcinoma

Misfolded proteins have enhanced formation of toxic oligomers and nonfunctional protein copies lead to recruiting wild-type protein types. Heat shock protein 90 (HSP90) is a molecular chaperone generated by cells that are involved in many cellular functions through regulation of folding and/or localization of large multi-protein complexes as well as client proteins. HSP90 can regulate a number of different cellular processes including cell proliferation, motility, angiogenesis, signal transduction, and adaptation to stress. HSP90 makes the mutated oncoproteins able to avoid misfolding and degradation and permits the malignant transformation. As a result, HSP90 is an important factor in several signaling pathways associated with tumorigenicity, therapy resistance, and inhibiting apoptosis. Clinically, the upregulation of HSP90 expression in hepatocellular carcinoma (HCC) is linked with advanced stages and inappropriate survival in cases suffering from this kind of cancer. The present review comprehensively assesses HSP90 functions and its possible usefulness as a potential diagnostic biomarker and therapeutic option for HCC.

Effects of ordered mutations on dynamics in signaling networks

BACKGROUND

Many previous clinical studies have found that accumulated sequential mutations are statistically related to tumorigenesis. However, they are limited in fully elucidating the significance of the ordered-mutation because they did not focus on the network dynamics. Therefore, there is a pressing need to investigate the dynamics characteristics induced by ordered-mutations.

METHODS

To quantify the ordered-mutation-inducing dynamics, we defined the mutation-sensitivity and the order-specificity that represent if the network is sensitive against a double knockout mutation and if mutation-sensitivity is specific to the mutation order, respectively, using a Boolean network model.

RESULTS

Through intensive investigations, we found that a signaling network is more sensitive when a double-mutation occurs in the direction order inducing a longer path and a smaller number of paths than in the reverse order. In addition, feedback loops involving a gene pair decreased both the mutation-sensitivity and the order-specificity. Next, we investigated relationships of functionally important genes with ordered-mutation-inducing dynamics. The network is more sensitive to mutations subject to drug-targets, whereas it is less specific to the mutation order. Both the sensitivity and specificity are increased when different-drug-targeted genes are mutated. Further, we found that tumor suppressors can efficiently suppress the amplification of oncogenes when the former are mutated earlier than the latter.

CONCLUSION

Taken together, our results help to understand the importance of the order of mutations with respect to the dynamical effects in complex biological systems.

Identification of neutral tumor evolution across cancer types

Despite extraordinary efforts to profile cancer genomes, interpreting the vast amount of genomic data in the light of cancer evolution remains challenging. Here we demonstrate that neutral tumor evolution results in a power-law distribution of the mutant allele frequencies reported by next-generation sequencing of tumor bulk samples. We find that the neutral power law fits with high precision 323 of 904 cancers from 14 types and from different cohorts. In malignancies identified as evolving neutrally, all clonal selection seemingly occurred before the onset of cancer growth and not in later-arising subclones, resulting in numerous passenger mutations that are responsible for intratumoral heterogeneity. Reanalyzing cancer sequencing data within the neutral framework allowed the measurement, in each patient, of both the in vivo mutation rate and the order and timing of mutations. This result provides a new way to interpret existing cancer genomic data and to discriminate between functional and non-functional intratumoral heterogeneity.

A computational model incorporating neural stem cell dynamics reproduces glioma incidence across the lifespan in the human population

Glioma is the most common form of primary brain tumor. Demographically, the risk of occurrence increases until old age. Here we present a novel computational model to reproduce the probability of glioma incidence across the lifespan. Previous mathematical models explaining glioma incidence are framed in a rather abstract way, and do not directly relate to empirical findings. To decrease this gap between theory and experimental observations, we incorporate recent data on cellular and molecular factors underlying gliomagenesis. Since evidence implicates the adult neural stem cell as the likely cell-of-origin of glioma, we have incorporated empirically-determined estimates of neural stem cell number, cell division rate, mutation rate and oncogenic potential into our model. We demonstrate that our model yields results which match actual demographic data in the human population. In particular, this model accounts for the observed peak incidence of glioma at approximately 80 years of age, without the need to assert differential susceptibility throughout the population. Overall, our model supports the hypothesis that glioma is caused by randomly-occurring oncogenic mutations within the neural stem cell population. Based on this model, we assess the influence of the (experimentally indicated) decrease in the number of neural stem cells and increase of cell division rate during aging. Our model provides multiple testable predictions, and suggests that different temporal sequences of oncogenic mutations can lead to tumorigenesis. Finally, we conclude that four or five oncogenic mutations are sufficient for the formation of glioma.