Dataset Information

Cancer progression models and fitness landscapes: a many-to-many relationship.

ABSTRACT: Motivation:The identification of constraints, due to gene interactions, in the order of accumulation of mutations during cancer progression can allow us to single out therapeutic targets. Cancer progression models (CPMs) use genotype frequency data from cross-sectional samples to identify these constraints, and return Directed Acyclic Graphs (DAGs) of restrictions where arrows indicate dependencies or constraints. On the other hand, fitness landscapes, which map genotypes to fitness, contain all possible paths of tumor progression. Thus, we expect a correspondence between DAGs from CPMs and the fitness landscapes where evolution happened. But many fitness landscapes-e.g. those with reciprocal sign epistasis-cannot be represented by CPMs. Results:Using simulated data under 500 fitness landscapes, I show that CPMs' performance (prediction of genotypes that can exist) degrades with reciprocal sign epistasis. There is large variability in the DAGs inferred from each landscape, which is also affected by mutation rate, detection regime and fitness landscape features, in ways that depend on CPM method. Using three cancer datasets, I show that these problems strongly affect the analysis of empirical data: fitness landscapes that are widely different from each other produce data similar to the empirically observed ones and lead to DAGs that infer very different restrictions. Because reciprocal sign epistasis can be common in cancer, these results question the use and interpretation of CPMs. Availability and implementation:Code available from Supplementary Material. Contact:ramon.diaz@iib.uam.es. Supplementary information:Supplementary data are available at Bioinformatics online.

SUBMITTER: Diaz-Uriarte R

PROVIDER: S-EPMC6031050 | biostudies-literature | 2018 Mar

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

Cancer progression models and fitness landscapes: a many-to-many relationship.

Diaz-Uriarte Ramon R

Bioinformatics (Oxford, England) 20180301 5

<h4>Motivation</h4>The identification of constraints, due to gene interactions, in the order of accumulation of mutations during cancer progression can allow us to single out therapeutic targets. Cancer progression models (CPMs) use genotype frequency data from cross-sectional samples to identify these constraints, and return Directed Acyclic Graphs (DAGs) of restrictions where arrows indicate dependencies or constraints. On the other hand, fitness landscapes, which map genotypes to fitness, con ...[more]

PMID: 29048486

Similar Datasets

Project description:Recombination is a fundamental biological process with profound evolutionary implications. Theory predicts that recombination increases the effectiveness of selection in natural populations. Yet, direct tests of this prediction have been restricted to qualitative trends due to the lack of detailed characterization of recombination rate variation across genomes and within species. The use of imprecise recombination rates can also skew population genetic analyses designed to assess the presence and mode of selection across genomes. Here we report the first integrated high-resolution description of genomic and population variation in recombination, which also distinguishes between the two outcomes of meiotic recombination: crossing over (CO) and gene conversion (GC). We characterized the products of 5,860 female meioses in Drosophila melanogaster by genotyping a total of 139 million informative SNPs and mapped 106,964 recombination events at a resolution down to 2 kilobases. This approach allowed us to generate whole-genome CO and GC maps as well as a detailed description of variation in recombination among individuals of this species. We describe many levels of variation in recombination rates. At a large-scale (100 kb), CO rates exhibit extreme and highly punctuated variation along chromosomes, with hot and coldspots. We also show extensive intra-specific variation in CO landscapes that is associated with hotspots at low frequency in our sample. GC rates are more uniformly distributed across the genome than CO rates and detectable in regions with reduced or absent CO. At a local scale, recombination events are associated with numerous sequence motifs and tend to occur within transcript regions, thus suggesting that chromatin accessibility favors double-strand breaks. All these non-independent layers of variation in recombination across genomes and among individuals need to be taken into account in order to obtain relevant estimates of recombination rates, and should be included in a new generation of population genetic models of the interaction between selection and linkage.

Project description:BACKGROUND: Tumorigenesis requires multiple genetic changes. Mutator mutations are mutations that increase genomic instability, and according to the mutator hypothesis, accelerate tumorigenesis by facilitating oncogenic mutations. Alternatively, repeated lineage selection and expansion without increased mutation frequency may explain observed cancer incidence. Mutator lineages also risk increased deleterious mutations, leading to extinction, thus providing another counterargument to the mutator hypothesis. Both selection and extinction involve changes in lineage fitness, which may be represented as "trajectories" through a "fitness landscape" defined by genetics and environment. METHODOLOGY/PRINCIPAL FINDINGS: Here I systematically analyze the relative efficiency of tumorigenesis with and without mutator mutations by evaluating archetypal fitness trajectories using deterministic and stochastic mathematical models. I hypothesize that tumorigenic mechanisms occur clinically in proportion to their relative efficiency. This work quantifies the relative importance of mutator pathways as a function of experimentally measurable parameters, demonstrating that mutator pathways generally enhance efficiency of tumorigenesis. An optimal mutation rate for tumor evolution is derived, and shown to differ from that for species evolution. CONCLUSIONS/SIGNIFICANCE: The models address the major counterarguments to the mutator hypothesis, confirming that mutator mechanisms are generally more efficient routes to tumorigenesis than non-mutator mechanisms. Mutator mutations are more likely to occur early, and to occur when more oncogenic mutations are required to create a tumor. Mutator mutations likely occur in a minority of premalignant lesions, but these mutator premalignant lesions are disproportionately likely to develop into malignant tumors. Tumor heterogeneity due to mutator mutations may contribute to therapeutic resistance, and the degree of heterogeneity of tumors may need to be considered when therapeutic strategies are devised. The model explains and predicts important biological observations in bacterial and mouse systems, as well as clinical observations.

Project description:Loss of habitat, specifically deforestation, is a major driver of biodiversity loss. Species-area relationship (SAR) models traditionally have been used for estimating species richness, species loss as a function of habitat loss, and extrapolation of richness for given areas. Sampling-species relationships (SSRs) are interrelated yet separate drivers for species richness estimates. Traditionally, however, SAR and SSR models have been used independently and not incorporated into a single approach. We developed and compared predictive models that incorporate sampling effort species-area relationships (SESARS) along the entire Atlantic Forest of South America, and then applied the best-fit model to estimate richness in forest remnants of Interior Atlantic Forest of eastern Paraguay. This framework was applied to non-volant small mammal assemblages that reflect different tolerances to forest loss and fragmentation. In order to account for differences in functionality we estimated small mammal richness of 1) the entire non-volant small mammal assemblage, including introduced species; 2) the native species forest assemblage; and 3) the forest-specialist assemblage, with the latter two assemblages being subsets of the entire assemblage. Finally, we geospatially modeled species richness for each of the three assemblages throughout eastern Paraguay to identify remnants with high species richness. We found that multiple regression power-law interaction-term models that only included area and the interactions of area and sampling as predictors, worked best for predicting species richness for the entire assemblage and the native species forest assemblage, while several traditional SAR models (logistic, power, exponential, and ratio) best described forest-specialist richness. Species richness was significantly different between assemblages. We identified obvious remnants with high species richness in eastern Paraguay, and these remnants often were geographically isolated. We also found relatively high predicted species richness (in relation to the entire range of predicted richness values) in several geographically-isolated, medium-size forest remnants that likely have not been considered as possible priority areas for conservation. These findings highlight the importance of using an empirical dataset, created using sources representing diverse sampling efforts, to develop robust predictive models. This approach is particularly important in geographic locations where field sampling is limited yet the geographic area is experiencing rapid and dramatic land cover changes. When combined, area and sampling are powerful modeling predictors for questions of biogeography, ecology, and conservation, especially when addressing habitat loss and fragmentation.

Dataset Information

Cancer progression models and fitness landscapes: a many-to-many relationship.

Publications

Cancer progression models and fitness landscapes: a many-to-many relationship.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets