Moritz Gerstung

The temporal order of genetic and pathway alterations in tumorigenesis.

Cancer evolves through the accumulation of mutations, but the order in which mutations occur is poorly understood. Inference of a temporal ordering on the level of genes is challenging because clinically and histologically identical tumors often have few mutated genes in common. This heterogeneity may at least in part be due to mutations in different genes having similar phenotypic effects by acting in the same functional pathway. We estimate the constraints on the order in which alterations accumulate during cancer progression from cross-sectional mutation data using a probabilistic graphical model termed Hidden Conjunctive Bayesian Network (H-CBN). The possible orders are analyzed on the level of genes and, after mapping genes to functional pathways, also on the pathway level. We find stronger evidence for pathway order constraints than for gene order constraints, indicating that temporal ordering results from selective pressure acting at the pathway level. The accumulation of changes in core pathways differs among cancer types, yet a common feature is that progression appears to begin with mutations in genes that regulate apoptosis pathways and to conclude with mutations in genes involved in invasion pathways. H-CBN models provide a quantitative and intuitive model of tumorigenesis showing that the genetic events can be linked to the phenotypic progression on the level of pathways.

Quantifying cancer progression with conjunctive Bayesian networks

Motivation: Cancer is an evolutionary process characterized by accumulating mutations. However, the precise timing and the order of genetic alterations that drive tumor progression remain enigmatic. Results: We present a specific probabilistic graphical model for the accumulation of mutations and their interdependencies. The Bayesian network models cancer progression by an explicit unobservable accumulation process in time that is separated from the observable but error-prone detection of mutations. Model parameters are estimated by an Expectation-Maximization algorithm and the underlying interaction graph is obtained by a simulated annealing procedure. Applying this method to cytogenetic data for different cancer types, we find multiple complex oncogenetic pathways deviating substantially from simplified models, such as linear pathways or trees. We further demonstrate how the inferred progression dynamics can be used to improve genetics-based survival predictions which could support diagnostics and prognosis. Availability: The software package ct-cbn is available under a GPL license on the web site cbg.ethz.ch/software/ct-cbn Contact: moritz.gerstung@bsse.ethz.ch

A Quantitative and Dynamic Model for Plant Stem Cell Regulation

Plants maintain pools of totipotent stem cells throughout their entire life. These stem cells are embedded within specialized tissues called meristems, which form the growing points of the organism. The shoot apical meristem of the reference plant Arabidopsis thaliana is subdivided into several distinct domains, which execute diverse biological functions, such as tissue organization, cell-proliferation and differentiation. The number of cells required for growth and organ formation changes over the course of a plants life, while the structure of the meristem remains remarkably constant. Thus, regulatory systems must be in place, which allow for an adaptation of cell proliferation within the shoot apical meristem, while maintaining the organization at the tissue level. To advance our understanding of this dynamic tissue behavior, we measured domain sizes as well as cell division rates of the shoot apical meristem under various environmental conditions, which cause adaptations in meristem size. Based on our results we developed a mathematical model to explain the observed changes by a cell pool size dependent regulation of cell proliferation and differentiation, which is able to correctly predict CLV3 and WUS over-expression phenotypes. While the model shows stem cell homeostasis under constant growth conditions, it predicts a variation in stem cell number under changing conditions. Consistent with our experimental data this behavior is correlated with variations in cell proliferation. Therefore, we investigate different signaling mechanisms, which could stabilize stem cell number despite variations in cell proliferation. Our results shed light onto the dynamic constraints of stem cell pool maintenance in the shoot apical meristem of Arabidopsis in different environmental conditions and developmental states.

Genomic Profiling of Advanced-Stage Oral Cancers Reveals Chromosome 11q Alterations as Markers of Poor Clinical Outcome

Identifying oral cancer lesions associated with high risk of relapse and predicting clinical outcome remain challenging questions in clinical practice. Genomic alterations may add prognostic information and indicate biological aggressiveness thereby emphasizing the need for genome-wide profiling of oral cancers. High-resolution array comparative genomic hybridization was performed to delineate the genomic alterations in clinically annotated primary gingivo-buccal complex and tongue cancers (nu200a=u200a60). The specific genomic alterations so identified were evaluated for their potential clinical relevance. Copy-number changes were observed on chromosomal arms with most frequent gains on 3q (60%), 5p (50%), 7p (50%), 8q (73%), 11q13 (47%), 14q11.2 (47%), and 19p13.3 (58%) and losses on 3p14.2 (55%) and 8p (83%). Univariate statistical analysis with correction for multiple testing revealed chromosomal gain of region 11q22.1–q22.2 and losses of 17p13.3 and 11q23–q25 to be associated with loco-regional recurrence (Pu200a=u200a0.004, Pu200a=u200a0.003, and Pu200a=u200a0.0003) and shorter survival (Pu200a=u200a0.009, Pu200a=u200a0.003, and P 0.0001) respectively. The gain of 11q22 and loss of 11q23-q25 were validated by interphase fluorescent in situ hybridization (I-FISH). This study identifies a tractable number of genomic alterations with few underlying genes that may potentially be utilized as biological markers for prognosis and treatment decisions in oral cancers.

Kinetic, Energetic, and Mechanical Differences between Dark-State Rhodopsin and Opsin

Rhodopsin, the photoreceptor pigment of the retina, initiates vision upon photon capture by its covalently linked chromophore 11-cis-retinal. In the absence of light, the chromophore serves as an inverse agonist locking the receptor in the inactive dark state. In the absence of chromophore, the apoprotein opsin shows low-level constitutive activity. Toward revealing insight into receptor properties controlled by the chromophore, we applied dynamic single-molecule force spectroscopy to quantify the kinetic, energetic, and mechanical differences between dark-state rhodopsin and opsin in native membranes from the retina of mice. Both rhodopsin and opsin are stabilized by ten structural segments. Compared to dark-state rhodopsin, the structural segments stabilizing opsin showed higher interaction strengths and mechanical rigidities and lower conformational variabilities, lifetimes, and free energies. These changes outline a common mechanism toward activating G-protein-coupled receptors. Additionally, we detected that opsin was more pliable and frequently stabilized alternate structural intermediates.

A High-Density Map for Navigating the Human Polycomb Complexome

Polycomb group (PcG) proteins are major determinants of gene silencing and epigenetic memory in higher eukaryotes. Here, we systematically mapped the human PcG complexome using a robust affinity purification mass spectrometry approach. Our high-density protein interaction network uncovered a diverse range of PcG complexes. Moreover, our analysis identified PcG interactors linking them to the PcG system, thus providing insight into the molecular function of PcG complexes and mechanisms of recruitment to target genes. We identified two human PRC2 complexes and two PR-DUB deubiquitination complexes, which contain the O-linked N-acetylglucosamine transferase OGT1 and several transcription factors. Finally, genome-wide profiling of PR-DUB components indicated that the human PR-DUB and PRC1 complexes bind distinct sets of target genes, suggesting differential impact on cellular processes in mammals.

Genome-wide expression and copy number analysis identifies driver genes in gingivobuccal cancers.

The molecular mechanisms contributing to the development and progression of gingivobuccal complex (GBC) cancers—a sub‐site of oral cancer, comprising the buccal mucosa, the gingivobuccal sulcus, the lower gingival region, and the retromolar trigone–remain poorly understood. Identifying the GBC cancer‐related gene expression signature and the driver genes residing on the altered chromosomal regions is critical for understanding the molecular basis of its pathogenesis. Genome‐wide expression profiling of 27 GBC cancers with known chromosomal alterations was performed to reveal differentially expressed genes. Putative driver genes were identified by integrating copy number and gene expression data. A total of 315 genes were found differentially expressed (P ≤ 0.05, logFC > 2.0) of which 11 genes were validated by real‐time quantitative reverse transcriptase‐PCR (qRT‐PCR) in tumors (n = 57) and normal GBC tissues (n = 18). Overexpression of LY6K, in chromosome band 8q24.3, was validated by immunohistochemical (IHC) analysis. We found that 78.5% (2,417/3,079) of the genes located in regions of recurrent chromosomal alterations show copy number dependent expression indicating that copy number alteration has a direct effect on global gene expression. The integrative analysis revealed BIRC3 in 11q22.2 as a candidate driver gene associated with poor clinical outcome. Our study identified previously unreported differentially expressed genes in a homogeneous subtype of oral cancer and the candidate driver genes that may contribute to the development and progression of the disease.

Waiting Time Models of Cancer Progression

Cancer progression is an evolutionary process driven by mutation and selection in a population of tumor cells. In multistage models of cancer progression, each stage is associated with the occurrence of genetic alterations and their fixation in the population. The accumulation of mutations is described using conjunctive Bayesian networks, an exponential family of waiting time models in which the occurrence of mutations is constrained by a partial temporal order. Two opposing limit cases arise if mutations either follow a linear order or occur independently. Exact analytical expressions for the waiting time until a specific number of mutations have accumulated are derived in these limit cases as well as for the general conjunctive Bayesian network. In a stochastic population genetics model that accounts for mutation and selection, waves of clonal expansions sweep through the population at equidistant intervals. An approximate analytical expression for the waiting time is compared to the results obtained with conjunctive Bayesian networks.

Evolutionary Games with Affine Fitness Functions: Applications to Cancer

We analyze the dynamics of evolutionary games in which fitness is defined as an affine function of the expected payoff and a constant contribution. The resulting inhomogeneous replicator equation has an homogeneous equivalent with modified payoffs. The affine terms also influence the stochastic dynamics of a two-strategy Moran model of a finite population. We then apply the affine fitness function in a model for tumor–normal cell interactions to determine which are the most successful tumor strategies. In order to analyze the dynamics of concurrent strategies within a tumor population, we extend the model to a three-strategy game involving distinct tumor cell types as well as normal cells. In this model, interaction with normal cells, in combination with an increased constant fitness, is the most effective way of establishing a population of tumor cells in normal tissue.

The BET protein FSH functionally interacts with ASH1 to orchestrate global gene activity in Drosophila

BackgroundThe question of how cells re-establish gene expression states after cell division is still poorly understood. Genetic and molecular analyses have indicated that Trithorax group (TrxG) proteins are critical for the long-term maintenance of active gene expression states in many organisms. A generally accepted model suggests that TrxG proteins contribute to maintenance of transcription by protecting genes from inappropriate Polycomb group (PcG)-mediated silencing, instead of directly promoting transcription.Results and discussionHere we report a physical and functional interaction in Drosophila between two members of the TrxG, the histone methyltransferase ASH1 and the bromodomain and extraterminal family protein FSH. We investigated this interface at the genome level, uncovering a widespread co-localization of both proteins at promoters and PcG-bound intergenic elements. Our integrative analysis of chromatin maps and gene expression profiles revealed that the observed ASH1-FSH binding pattern at promoters is a hallmark of active genes. Inhibition of FSH-binding to chromatin resulted in global down-regulation of transcription. In addition, we found that genes displaying marks of robust PcG-mediated repression also have ASH1 and FSH bound to their promoters.ConclusionsOur data strongly favor a global coactivator function of ASH1 and FSH during transcription, as opposed to the notion that TrxG proteins impede inappropriate PcG-mediated silencing, but are dispensable elsewhere. Instead, our results suggest that PcG repression needs to overcome the transcription-promoting function of ASH1 and FSH in order to silence genes.