Mutator Mutations Enhance Tumorigenic Efficiency across Fitness Landscapes

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.     

The strength of genetic interactions scales weakly with mutational effects

Background

Genetic interactions pervade every aspect of biology, from evolutionary theory, where they determine the accessibility of evolutionary paths, to medicine, where they can contribute to complex genetic diseases. Until very recently, studies on epistatic interactions have been based on a handful of mutations, providing at best anecdotal evidence about the frequency and the typical strength of genetic interactions. In this study, we analyze a publicly available dataset that contains the growth rates of over five million double knockout mutants of the yeast Saccharomyces cerevisiae.

Results

We discuss a geometric definition of epistasis that reveals a simple and surprisingly weak scaling law for the characteristic strength of genetic interactions as a function of the effects of the mutations being combined. We then utilized this scaling to quantify the roughness of naturally occurring fitness landscapes. Finally, we show how the observed roughness differs from what is predicted by Fisher''s geometric model of epistasis, and discuss the consequences for evolutionary dynamics.

Conclusions

Although epistatic interactions between specific genes remain largely unpredictable, the statistical properties of an ensemble of interactions can display conspicuous regularities and be described by simple mathematical laws. By exploiting the amount of data produced by modern high-throughput techniques, it is now possible to thoroughly test the predictions of theoretical models of genetic interactions and to build informed computational models of evolution on realistic fitness landscapes.     

A Many-Body Field Theory Approach to Stochastic Models in Population Biology

Background

Many models used in theoretical ecology, or mathematical epidemiology are stochastic, and may also be spatially-explicit. Techniques from quantum field theory have been used before in reaction-diffusion systems, principally to investigate their critical behavior. Here we argue that they make many calculations easier and are a possible starting point for new approximations.

Methodology

We review the many-body field formalism for Markov processes and illustrate how to apply it to a ‘Brownian bug’ population model, and to an epidemic model. We show how the master equation and the moment hierarchy can both be written in particularly compact forms. The introduction of functional methods allows the systematic computation of the effective action, which gives the dynamics of mean quantities. We obtain the 1-loop approximation to the effective action for general (space-) translation invariant systems, and thus approximations to the non-equilibrium dynamics of the mean fields.

Conclusions

The master equations for spatial stochastic systems normally take a neater form in the many-body field formalism. One can write down the dynamics for generating functional of physically-relevant moments, equivalent to the whole moment hierarchy. The 1-loop dynamics of the mean fields are the same as those of a particular moment-closure.     

In Vitro HIV-1 Evolution in Response to Triple Reverse Transcriptase Inhibitors & In Silico Phenotypic Analysis

Background

Effectiveness of ART regimens strongly depends upon complex interactions between the selective pressure of drugs and the evolution of mutations that allow or restrict drug resistance.

Methods

Four clinical isolates from NRTI-exposed, NNRTI-naive subjects were passaged in increasing concentrations of NVP in combination with 1 µM 3 TC and 2 µM ADV to assess selective pressures of multi-drug treatment. A novel parameter inference procedure, based on a stochastic viral growth model, was used to estimate phenotypic resistance and fitness from in vitro combination passage experiments.

Results

Newly developed mathematical methods estimated key phenotypic parameters of mutations arising through selective pressure exerted by 3 TC and NVP. Concentrations of 1 µM 3 TC maintained the M184V mutation, which was associated with intrinsic fitness deficits. Increasing NVP concentrations selected major NNRTI resistance mutations. The evolutionary pathway of NVP resistance was highly dependent on the viral genetic background, epistasis as well as stochasticity. Parameter estimation indicated that the previously unrecognized mutation L228Q was associated with NVP resistance in some isolates.

Conclusion

Serial passage of viruses in the presence of multiple drugs may resemble the selection of mutations observed among treated individuals and populations in vivo and indicate evolutionary preferences and restrictions. Phenotypic resistance estimated here “in silico” from in vitro passage experiments agreed well with previous knowledge, suggesting that the unique combination of “wet-” and “dry-lab” experimentation may improve our understanding of HIV-1 resistance evolution in the future.     

Divergent selection in trailing- versus leading-edge populations of Biscutella laevigata

Background and Aims

Knowledge on how climate-induced range shifts might affect natural selection is crucial to understand the evolution of species ranges.

Methods

Using historical demographic perspectives gathered from regional-scale phylogeography on the alpine herb Biscutella laevigata, indirect inferences on gene flow and signature of selection based on AFLP genotyping were compared between local populations persisting at the trailing edge and expanding at the leading edge.

Key Results

Spatial autocorrelation revealed that gene flow was two times more restricted at the trailing edge and genome scans indicated divergent selection in this persisting population. In contrast, no pattern of selection emerged in the expanding population at the leading edge.

Conclusions

Historical effects may determine different architecture of genetic variation and selective patterns within local populations, what is arguably important to understand evolutionary processes acting across the species ranges.     

Quantitative Trait Evolution and Environmental Change

Background

Given the recent changes in climate, there is an urgent need to understand the evolutionary ability of populations to respond to these changes.

Methodology/Principal Findings

We performed individual-based simulations with different shapes of the fitness curve, different heritabilities, different levels of density compensation, and different autocorrelation of environmental noise imposed on an environmental trend to study the ability of a population to adapt to changing conditions. The main finding is that when there is a positive autocorrelation of environmental noise, the outcome of the evolutionary process is much more unpredictable compared to when the noise has no autocorrelation. In addition, we found that strong selection resulted in a higher load, and more extinctions, and that this was most pronounced when heritability was low. The level of density-compensation was important in determining the variance in load when there was strong selection, and when genetic variance was lower when the level of density-compensation was low.

Conclusions

The strong effect of the details of the environmental fluctuations makes predictions concerning the evolutionary future of populations very hard to make. In addition, to be able to make good predictions we need information on heritability, fitness functions and levels of density compensation. The results strongly suggest that patterns of environmental noise must be incorporated in future models of environmental change, such as global warming.     

Viral and epidemiological determinants of the invasion dynamics of novel dengue genotypes

Background

Dengue has become a major concern for international public health. Frequent epidemic outbreaks are believed to be driven by a complex interplay of immunological interactions between its four co-circulating serotypes and large fluctuations in mosquito densities. Viral lineage replacement events, caused for example by different levels of cross-protection or differences in viral fitness, have also been linked to a temporary change in dengue epidemiology. A major replacement event was recently described for South-East Asia where the Asian-1 genotype of dengue serotype 2 replaced the resident Asian/American type. Although this was proposed to be due to increased viral fitness in terms of enhanced human-to-mosquito transmission, no major change in dengue epidemiology could be observed.

Methods/Results

Here we investigate the invasion dynamics of a novel, advantageous dengue genotype within a model system and determine the factors influencing the success and rate of fixation as well as their epidemiological consequences. We find that while viral fitness overall correlates with invasion success and competitive exclusion of the resident genotype, the epidemiological landscape plays a more significant role for successful emergence. Novel genotypes can thus face high risks of stochastic extinction despite their fitness advantage if they get introduced during episodes of high dengue prevalence, especially with respect to that particular serotype.

Conclusion

The rarity of markers for positive selection has often been explained by strong purifying selection whereby the constraints imposed by dengue''s two-host cycle are expected to result in a high rate of deleterious mutations. Our results demonstrate that even highly beneficial mutants are under severe threat of extinction, which would suggest that apart from purifying selection, stochastic effects and genetic drift beyond seasonal bottlenecks are equally important in shaping dengue''s viral ecology and evolution.     

Evidence for Repeated Independent Evolution of Migration in the Largest Family of Bats

Background

How migration evolved represents one of the most poignant questions in evolutionary biology. While studies on the evolution of migration in birds are well represented in the literature, migration in bats has received relatively little attention. Yet, more than 30 species of bats are known to migrate annually from breeding to non-breeding locations. Our study is the first to test hypotheses on the evolutionary history of migration in bats using a phylogenetic framework.

Methods and Principal Findings

In addition to providing a review of bat migration in relation to existing hypotheses on the evolution of migration in birds, we use a previously published supertree to formulate and test hypotheses on the evolutionary history of migration in bats. Our results suggest that migration in bats has evolved independently in several lineages potentially as the need arises to track resources (food, roosting site) but not through a series of steps from short- to long-distance migrants, as has been suggested for birds. Moreover, our analyses do not indicate that migration is an ancestral state but has relatively recently evolved in bats. Our results also show that migration is significantly less likely to evolve in cave roosting bats than in tree roosting species.

Conclusions and Significance

This is the first study to provide evidence that migration has evolved independently in bat lineages that are not closely related. If migration evolved as a need to track seasonal resources or seek adequate roosting sites, climate change may have a pivotal impact on bat migratory habits. Our study provides a strong framework for future research on the evolution of migration in chiropterans.     

Simulation of the evolution of root water foraging strategies in dry and shallow soils

Background and Aims

The dynamic structural development of plants can be seen as a strategy for exploiting the limited resources available within their environment, and we would expect that evolution would lead to efficient strategies that reduce costs while maximizing resource acquisition. In particular, perennial species endemic to habitats with shallow soils in seasonally dry environments have been shown to have a specialized root system morphology that may enhance access to water resources in the underlying rock. This study aimed to explore these hypotheses by applying evolutionary algorithms to a functional–structural root growth model.

Methods

A simulation model of a plant''s root system was developed, which represents the dynamics of water uptake and structural growth. The model is simple enough for evolutionary optimization to be computationally feasible, yet flexible enough to allow a range of structural development strategies to be explored. The model was combined with an evolutionary algorithm in order to investigate a case study habitat with a highly heterogeneous distribution of resources, both spatially and temporally – the situation of perennial plants occurring on shallow soils in seasonally dry environments. Evolution was simulated under two contrasting fitness criteria: (1) the ability to find wet cracks in underlying rock, and (2) maximizing above-ground biomass.

Key Results

The novel approach successfully resulted in the evolution of more efficient structural development strategies for both fitness criteria. Different rooting strategies evolved when different criteria were applied, and each evolved strategy made ecological sense in terms of the corresponding fitness criterion. Evolution selected for root system morphologies which matched those of real species from corresponding habitats.

Conclusions

Specialized root morphology with deeper rather than shallower lateral branching enhances access to water resources in underlying rock. More generally, the approach provides insights into both evolutionary processes and ecological costs and benefits of different plant growth strategies.     

The potential of plant viruses to promote genotypic diversity via genotype x environment interactions

Background and Aims

Genotype by environment (G × E) interactions are important for the long-term persistence of plant species in heterogeneous environments. It has often been suggested that disease is a key factor for the maintenance of genotypic diversity in plant populations. However, empirical evidence for this contention is scarce. Here virus infection is proposed as a possible candidate for maintaining genotypic diversity in their host plants.

Methods

The effects of White clover mosaic virus (WClMV) on the performance and development of different Trifolium repens genotypes were analysed and the G × E interactions were examined with respect to genotype-specific plant responses to WClMV infection. Thus, the environment is defined as the presence or absence of the virus.

Key Results

WClMV had a negative effect on plant performance as shown by a decrease in biomass and number of ramets. These effects of virus infection differ greatly among host genotypes, representing a strong G × E interaction. Moreover, the relative fitness and associated ranking of genotypes changed significantly between control and virus treatments. This shift in relative fitness among genotypes suggests the potential for WClMV to provoke differential selection on T. repens genotypes, which may lead to negative frequency-dependent selection in host populations.

Conclusions

The apparent G × E interaction and evident repercussions for relative fitness reported in this study stress the importance of viruses for ecological and evolutionary processes and suggest an important role for viruses in shaping population dynamics and micro-evolutionary processes.     

Genome size reduction can trigger rapid phenotypic evolution in invasive plants

Background and Aims

The study of rapid evolution in invasive species has highlighted the fundamental role played by founder events, emergence of genetic novelties through recombination and rapid response to new selective pressures. However, whether rapid adaptation of introduced species can be driven by punctual changes in genome organization has received little attention. In plants, variation in genome size, i.e. variation in the amount of DNA per monoploid set of chromosomes through loss or gain of repeated DNA sequences, is known to influence a number of physiological, phenological and life-history features. The present study investigated whether change in genome size has contributed to the evolution of greater potential of vegetative growth in invasive populations of an introduced grass.

Methods

The study was based on the recent demonstration that invasive genotypes of reed canarygrass (Phalaris arundinacea) occurring in North America have emerged from recombination between introduced European strains. The genome sizes of more than 200 invasive and native genotypes were measured and their genome size was related to their phenotypic traits measured in a common glasshouse environment. Population genetics data were used to infer phylogeographical relationships between study populations, and the evolutionary history of genome size within the study species was inferred.

Key Results

Invasive genotypes had a smaller genome than European native genotypes from which they are derived. This smaller genome size had phenotypic effects that increased the species'' invasive potential, including a higher early growth rate, due to a negative relationship between genome size and rate of stem elongation. Based on inferred phylogeographical relationships of invasive and native populations, evolutionary models were consistent with a scenario of genome reduction by natural selection during the invasion process, rather than a scenario of stochastic change.

Conclusions

Punctual reduction in genome size could cause rapid changes in key phenotypic traits that enhance invasive ability. Although the generality of genome size variation leading to phenotypic evolution and the specific genomic mechanisms involved are not known, change in genome size may constitute an important but previously under-appreciated mechanism of rapid evolutionary change that may promote evolutionary novelties over short time scales.Key words: Biological invasion, evolutionary models, genome size, Phalaris arundinacea, quantile regression, relative growth rate, rapid evolution     

Adaptive evolution and functional redesign of core metabolic proteins in snakes

Background

Adaptive evolutionary episodes in core metabolic proteins are uncommon, and are even more rarely linked to major macroevolutionary shifts.

Methodology/Principal Findings

We conducted extensive molecular evolutionary analyses on snake mitochondrial proteins and discovered multiple lines of evidence suggesting that the proteins at the core of aerobic metabolism in snakes have undergone remarkably large episodic bursts of adaptive change. We show that snake mitochondrial proteins experienced unprecedented levels of positive selection, coevolution, convergence, and reversion at functionally critical residues. We examined Cytochrome C oxidase subunit I (COI) in detail, and show that it experienced extensive modification of normally conserved residues involved in proton transport and delivery of electrons and oxygen. Thus, adaptive changes likely altered the flow of protons and other aspects of function in CO, thereby influencing fundamental characteristics of aerobic metabolism. We refer to these processes as “evolutionary redesign” because of the magnitude of the episodic bursts and the degree to which they affected core functional residues.

Conclusions/Significance

The evolutionary redesign of snake COI coincided with adaptive bursts in other mitochondrial proteins and substantial changes in mitochondrial genome structure. It also generally coincided with or preceded major shifts in ecological niche and the evolution of extensive physiological adaptations related to lung reduction, large prey consumption, and venom evolution. The parallel timing of these major evolutionary events suggests that evolutionary redesign of metabolic and mitochondrial function may be related to, or underlie, the extreme changes in physiological and metabolic efficiency, flexibility, and innovation observed in snake evolution.     

A Minimal Model of Signaling Network Elucidates Cell-to-Cell Stochastic Variability in Apoptosis

Background

Signaling networks are designed to sense an environmental stimulus and adapt to it. We propose and study a minimal model of signaling network that can sense and respond to external stimuli of varying strength in an adaptive manner. The structure of this minimal network is derived based on some simple assumptions on its differential response to external stimuli.

Methodology

We employ stochastic differential equations and probability distributions obtained from stochastic simulations to characterize differential signaling response in our minimal network model. Gillespie''s stochastic simulation algorithm (SSA) is used in this study.

Conclusions/Significance

We show that the proposed minimal signaling network displays two distinct types of response as the strength of the stimulus is decreased. The signaling network has a deterministic part that undergoes rapid activation by a strong stimulus in which case cell-to-cell fluctuations can be ignored. As the strength of the stimulus decreases, the stochastic part of the network begins dominating the signaling response where slow activation is observed with characteristic large cell-to-cell stochastic variability. Interestingly, this proposed stochastic signaling network can capture some of the essential signaling behaviors of a complex apoptotic cell death signaling network that has been studied through experiments and large-scale computer simulations. Thus we claim that the proposed signaling network is an appropriate minimal model of apoptosis signaling. Elucidating the fundamental design principles of complex cellular signaling pathways such as apoptosis signaling remains a challenging task. We demonstrate how our proposed minimal model can help elucidate the effect of a specific apoptotic inhibitor Bcl-2 on apoptotic signaling in a cell-type independent manner. We also discuss the implications of our study in elucidating the adaptive strategy of cell death signaling pathways.     

Functional evolution of leptin of Ochotona curzoniae in adaptive thermogenesis driven by cold environmental stress

Background

Environmental stress can accelerate the directional selection and evolutionary rate of specific stress-response proteins to bring about new or altered functions, enhancing an organism''s fitness to challenging environments. Plateau pika (Ochotona curzoniae), an endemic and keystone species on Qinghai-Tibetan Plateau, is a high hypoxia and low temperature tolerant mammal with high resting metabolic rate and non-shivering thermogenesis to cope in this harsh plateau environment. Leptin is a key hormone related to how these animals regulate energy homeostasis. Previous molecular evolutionary analysis helped to generate the hypothesis that adaptive evolution of plateau pika leptin may be driven by cold stress.

Methodology/Principal Findings

To test the hypothesis, recombinant pika leptin was first purified. The thermogenic characteristics of C57BL/6J mice injected with pika leptin under warm (23±1°C) and cold (5±1°C) acclimation is investigated. Expression levels of genes regulating adaptive thermogenesis in brown adipose tissue and the hypothalamus are compared between pika leptin and human leptin treatment, suggesting that pika leptin has adaptively and functionally evolved. Our results show that pika leptin regulates energy homeostasis via reduced food intake and increased energy expenditure under both warm and cold conditions. Compared with human leptin, pika leptin demonstrates a superior induced capacity for adaptive thermogenesis, which is reflected in a more enhanced β-oxidation, mitochondrial biogenesis and heat production. Moreover, leptin treatment combined with cold stimulation has a significant synergistic effect on adaptive thermogenesis, more so than is observed with a single cold exposure or single leptin treatment.

Conclusions/Significance

These findings support the hypothesis that cold stress has driven the functional evolution of plateau pika leptin as an ecological adaptation to the Qinghai-Tibetan Plateau.     

Molecular Evolution of the Primate Antiviral Restriction Factor Tetherin

Background

Tetherin is a recently identified antiviral restriction factor that restricts HIV-1 particle release in the absence of the HIV-1 viral protein U (Vpu). It is reminiscent of APOBEC3G and TRIM5a that also antagonize HIV. APOBEC3G and TRIM5a have been demonstrated to evolve under pervasive positive selection throughout primate evolution, supporting the red-queen hypothesis. Therefore, one naturally presumes that Tetherin also evolves under pervasive positive selection throughout primate evolution and supports the red-queen hypothesis. Here, we performed a detailed evolutionary analysis to address this presumption.

Methodology/Principal Findings

Results of non-synonymous and synonymous substitution rates reveal that Tetherin as a whole experiences neutral evolution rather than pervasive positive selection throughout primate evolution, as well as in non-primate mammal evolution. Sliding-window analyses show that the regions of the primate Tetherin that interact with viral proteins are under positive selection or relaxed purifying selection. In particular, the sites identified under positive selection generally focus on these regions, indicating that the main selective pressure acting on the primate Tetherin comes from virus infection. The branch-site model detected positive selection acting on the ancestral branch of the New World Monkey lineage, suggesting an episodic adaptive evolution. The positive selection was also found in duplicated Tetherins in ruminants. Moreover, there is no bias in the alterations of amino acids in the evolution of the primate Tetherin, implying that the primate Tetherin may retain broad spectrum of antiviral activity by maintaining structure stability.

Conclusions/Significance

These results conclude that the molecular evolution of Tetherin may be attributed to the host–virus arms race, supporting the Red Queen hypothesis, and Tetherin may be in an intermediate stage in transition from neutral to pervasive adaptive evolution.     

Corolla morphology influences diversification rates in bifid toadflaxes (Linaria sect. Versicolores)

Background and Aims

The role of flower specialization in plant speciation and evolution remains controversial. In this study the evolution of flower traits restricting access to pollinators was analysed in the bifid toadflaxes (Linaria sect. Versicolores), a monophyletic group of ∼30 species and subspecies with highly specialized corollas.

Methods

A time-calibrated phylogeny based on both nuclear and plastid DNA sequences was obtained using a coalescent-based method, and flower morphology was characterized by means of morphometric analyses. Directional trends in flower shape evolution and trait-dependent diversification rates were jointly analysed using recently developed methods, and morphological shifts were reconstructed along the phylogeny. Pollinator surveys were conducted for a representative sample of species.

Key Results

A restrictive character state (narrow corolla tube) was reconstructed in the most recent common ancestor of Linaria sect. Versicolores. After its early loss in the most species-rich clade, this character state has been convergently reacquired in multiple lineages of this clade in recent times, yet it seems to have exerted a negative influence on diversification rates. Comparative analyses and pollinator surveys suggest that the narrow- and broad-tubed flowers are evolutionary optima representing divergent strategies of pollen placement on nectar-feeding insects.

Conclusions

The results confirm that different forms of floral specialization can lead to dissimilar evolutionary success in terms of diversification. It is additionally suggested that opposing individual-level and species-level selection pressures may have driven the evolution of pollinator-restrictive traits in bifid toadflaxes.     

Floral polymorphism and the fitness implications of attracting pollinating and florivorous insects

Background and Aims

Floral polymorphism is frequently attributed to pollinator-mediated selection. Multiple studies, however, have revealed the importance of non-pollinating visitors in floral evolution. Using the polymorphic annual daisy Ursinia calenduliflora, this study investigated the importance of different insect visitors, and their effects on fitness, in the maintenance of floral polymorphism.

Methods

The spatial structure of a discrete floral polymorphism was characterized based on the presence/absence of anthocyanin floret spots in U. calenduliflora. A 3-year observational study was then conducted in polymorphic populations to investigate differences in visitation rates of dominant visitors to floral morphs. Experiments were performed to explore the floral preference of male and female Megapalpus capensis (the dominant insect visitor) and their effectiveness as pollinators. Next, floral damage by antagonistic florivores and the reproductive success of the two floral morphs were surveyed in multiple populations and years.

Key Results

Floral polymorphism in U. calenduliflora was structured spatially, as were insect visitation patterns. Megapalpus capensis males were the dominant visitors and exhibited strong preference for the spotted morph in natural and experimental observations. While this may indicate potential fitness benefits for the spotted morph, female fitness did not differ between floral morphs. However, as M. capensis males are very efficient at exporting U. calenduliflora pollen, their preference may likely increase the reproductive fitness of the spotted morph through male fitness components. The spotted morph, however, also suffered significantly greater costs due to ovule predation by florivores than the spotless morph.

Conclusions

The results suggest that pollinators and florivores may potentially exert antagonistic selection that could contribute to the maintenance of floral polymorphism across the range of U. calenduliflora. The relative strength of selection imposed by each agent is potentially determined by insect community composition and abundance at each site, highlighting the importance of community context in the evolution of floral phenotypes.