Clonal erosion and genetic drift in cyclical parthenogens – the interplay between neutral and selective processes

The occurrence of alternating phases of clonal and sexual reproduction may strongly impact the interplay between neutral and selective genetic variation in populations. Using a physiologically structured model of the life history of Daphnia, we investigated to what extent clonal erosion associated with selection during the clonal phase affects the genetic structure as observed by neutral markers. Incorporating conservative levels of quantitative genetic variation at 11 physiological and life history traits induces strong clonal erosion, reducing clonal diversity (CD) near the end of the simulations (1000 days) to a level between 1 and 5, even in habitats with high initial CD (10⁸ clones). This strong clonal erosion caused by selection can result in reduced genetic diversity, significant excess of heterozygotes and significant genetic differentiation between populations as observed by neutral markers. Our results indicate that, especially in relatively small habitats, clonal selection may strongly impact the genetic structure and may contribute to the often observed high level of neutral genetic differentiation among natural populations of cyclical parthenogens.     

The nearly neutral and selection theories of molecular evolution under the fisher geometrical framework: substitution rate, population size, and complexity

The general theories of molecular evolution depend on relatively arbitrary assumptions about the relative distribution and rate of advantageous, deleterious, neutral, and nearly neutral mutations. The Fisher geometrical model (FGM) has been used to make distributions of mutations biologically interpretable. We explored an FGM-based molecular model to represent molecular evolutionary processes typically studied by nearly neutral and selection models, but in which distributions and relative rates of mutations with different selection coefficients are a consequence of biologically interpretable parameters, such as the average size of the phenotypic effect of mutations and the number of traits (complexity) of organisms. A variant of the FGM-based model that we called the static regime (SR) represents evolution as a nearly neutral process in which substitution rates are determined by a dynamic substitution process in which the population's phenotype remains around a suboptimum equilibrium fitness produced by a balance between slightly deleterious and slightly advantageous compensatory substitutions. As in previous nearly neutral models, the SR predicts a negative relationship between molecular evolutionary rate and population size; however, SR does not have the unrealistic properties of previous nearly neutral models such as the narrow window of selection strengths in which they work. In addition, the SR suggests that compensatory mutations cannot explain the high rate of fixations driven by positive selection currently found in DNA sequences, contrary to what has been previously suggested. We also developed a generalization of SR in which the optimum phenotype can change stochastically due to environmental or physiological shifts, which we called the variable regime (VR). VR models evolution as an interplay between adaptive processes and nearly neutral steady-state processes. When strong environmental fluctuations are incorporated, the process becomes a selection model in which evolutionary rate does not depend on population size, but is critically dependent on the complexity of organisms and mutation size. For SR as well as VR we found that key parameters of molecular evolution are linked by biological factors, and we showed that they cannot be fixed independently by arbitrary criteria, as has usually been assumed in previous molecular evolutionary models.     

Mitochondrial DNA mutations and breast tumorigenesis

Breast cancer is a heterogeneous disease and genetic factors play an important role in its genesis. Although mutations in tumor suppressors and oncogenes encoded by the nuclear genome are known to play a critical role in breast tumorigenesis, the contribution of the mitochondrial genome to this process is unclear. Like the nuclear genome, the mitochondrial genome also encodes proteins critical for mitochondrion functions such as oxidative phosphorylation (OXPHOS), which is known to be defective in cancer including breast cancer. Mitochondrial DNA (mtDNA) is more susceptible to mutations due to limited repair mechanisms compared to nuclear DNA (nDNA). Thus changes in mitochondrial genes could also contribute to the development of breast cancer. In this review we discuss mtDNA mutations that affect OXPHOS. Continuous acquisition of mtDNA mutations and selection of advantageous mutations ultimately leads to generation of cells that propagate uncontrollably to form tumors. Since irreversible damage to OXPHOS leads to a shift in energy metabolism towards enhanced aerobic glycolysis in most cancers, mutations in mtDNA represent an early event during breast tumorigenesis, and thus may serve as potential biomarkers for early detection and prognosis of breast cancer. Because mtDNA mutations lead to defective OXPHOS, development of agents that target OXPHOS will provide specificity for preventative and therapeutic agents against breast cancer with minimal toxicity.     

The role of weak selection and high mutation rates in nearly neutral evolution

Neutral dynamics occur in evolution if all types are ‘effectively equal’ in their reproductive success, where the definition of ‘effectively equal’ depends on the population size and the details of mutations. Empirically observed neutral genetic evolution in extremely large clonal populations can only be explained under current models if selection is completely absent. Such models typically consider the case where population dynamics occurs on a different timescale to evolution. However, this assumption is invalid when mutations are not rare in a whole population. We show that this has important consequences for the occurrence of neutral evolution in clonal populations. In highly connected type spaces, neutral dynamics can occur for all population sizes despite significant selective differences, via the forming of effectively neutral networks connecting rare neutral types. Biological implications include an explanation for the high diversity of rare types that survive in large clonal populations, and a theoretical justification for the use of neutral null models.     

Natural and chemotherapy-induced clonal evolution of tumors

Evolution and natural selection of tumoral clones in the process of transformation and the following carcinogenesis can be called natural clonal evolution. Its main driving factors are internal: genetic instability initiated by driver mutations and microenvironment, which enables selective pressure while forming the environment for cell transformation and their survival. We present our overview of contemporary research dealing with mechanisms of carcinogenesis in different localizations from precancerous pathologies to metastasis and relapse. It shows that natural clonal evolution establishes intratumoral heterogeneity and enables tumor progression. Tumors of monoclonal origin are of low-level intratumoral heterogeneity in the initial stages, and this increases with the size of the tumor. Tumors of polyclonal origin are of extremely high-level intratumoral heterogeneity in the initial stages and become more homogeneous when larger due to clonal expansion. In cases of chemotherapy-induced clonal evolution of a tumor, chemotherapy becomes the leading factor in treatment. The latest research shows that the impact of chemotherapy can radically increase the speed of clonal evolution and lead to new malignant and resistant clones that cause tumor metastasis. Another option of chemotherapy-induced clonal evolution is formation of a new dominant clone from a clone that was minor in the initial tumor and obtained free space due to elimination of sensitive clones by chemotherapy. As a result, in ~20% of cases, chemotherapy can stimulate metastasis and relapse of tumors due to clonal evolution. The conclusion of the overview formulates approaches to tumor treatment based on clonal evolution: in particular, precision therapy, prediction of metastasis stimulation in patients treated with chemotherapy, methods of genetic evaluation of chemotherapy efficiency and clonal-oriented treatment, and approaches to manipulating the clonal evolution of tumors are presented.     

Distinguishing driver and passenger mutations in an evolutionary history categorized by interference

In many biological scenarios, from the development of drug resistance in pathogens to the progression of healthy cells toward cancer, quantifying the selection acting on observed mutations is a central question. One difficulty in answering this question is the complexity of the background upon which mutations can arise, with multiple potential interactions between genetic loci. We here present a method for discerning selection from a population history that accounts for interference between mutations. Given sequences sampled from multiple time points in the history of a population, we infer selection at each locus by maximizing a likelihood function derived from a multilocus evolution model. We apply the method to the question of distinguishing between loci where new mutations are under positive selection (drivers) and loci that emit neutral mutations (passengers) in a Wright-Fisher model of evolution. Relative to an otherwise equivalent method in which the genetic background of mutations was ignored, our method inferred selection coefficients more accurately for both driver mutations evolving under clonal interference and passenger mutations reaching fixation in the population through genetic drift or hitchhiking. In a population history recorded by 750 sets of sequences of 100 individuals taken at intervals of 100 generations, a set of 50 loci were divided into drivers and passengers with a mean accuracy of >0.95 across a range of numbers of driver loci. The potential application of our model, either in full or in part, to a range of biological systems, is discussed.     

An analysis of selection on a colour polymorphism in the northern leopard frog

In this study, we investigated the role of selection in the maintenance of a dorsal colour polymorphism in natural populations of the northern leopard frog, Rana pipiens. We determined genetic structure both spatially and temporally from a suite of putatively neutral molecular markers and tested whether or not the colour locus exhibited patterns of genetic variation that differed from those of the neutral loci. Spatial genetic structure at the colour locus was indistinguishable from structure at neutral loci [95% confidence intervals of F(ST) (neutral) = (0.07, 0.35), F(ST) (colour locus) = 0.114]. In the temporal analysis, we found that the variance among populations in the change in allele frequency at the colour locus (equal to 0.004) lies within the 95% confidence intervals for the variance among populations in changes in allele frequencies at neutral loci. In light of our inability to show evidence for the selective maintenance of the colour polymorphism, we used computer simulations to infer the power of our spatial and temporal techniques to detect selection. The computer simulations showed that although the strength of selection (s) would need to be relatively strong to have been detected by the temporal approach (s = 0.1-0.4, depending on the model tested), the spatial analysis would have detected all but weak selection (s = 0.01-0.04, depending on the model tested). This study illustrates the importance of using a locus comparison approach to detect evidence for selective maintenance before conducting studies to measure the selective mechanisms maintaining a polymorphism.     

Simulation Study of a Multigene Family, with Special Reference to the Evolution of Compensatory Advantageous Mutations

We investigate the evolution of a multigene family incorporating the forces of drift, mutation, gene conversion, unequal crossing over and selection. The use of simulation studies is required due to the complexity of the model. Selection is modeled in two modes: positive selection as a function of the number of different beneficial alleles and negative selection against deleterious alleles. We assume that gene conversion is unbiased, and that all mutations are initially deleterious. Compensation between mutants creates beneficial and neutral alleles, and allowances are made for compensatory mutations either within or between the members of a multigene family. We find that gene conversion can enhance the rate of acquisition of compensatory advantageous mutations when genes are redundant.     

Hitchhiking in space: Ancestry in adapting,spatially extended populations*

Selective sweeps reduce neutral genetic diversity. In sexual populations, this “hitchhiking” effect is thought to be limited to the local genomic region of the sweeping allele. While this is true in panmictic populations, we find that in spatially extended populations the combined effects of many unlinked sweeps can affect patterns of ancestry (and therefore neutral genetic diversity) across the whole genome. Even low rates of sweeps can be enough to skew the spatial locations of ancestors such that neutral mutations that occur in an individual living outside a small region in the center of the range have virtually no chance of fixing in the population. The fact that nearly all ancestry rapidly traces back to a small spatial region also means that relatedness between individuals falls off very slowly as a function of the spatial distance between them.     

On the evolution of clonal plant life histories

Clonal plant life histories are special in at least four respects: (1) Clonal plants can also reproduce vegetatively, (2) vegetative reproduction can be realised with short or long spacers, (3) and it may allow to plastically place vegetative offspring in benign patches. (4) Moreover, ramets of clonal plants may remain physically and physiologically integrated. Because of the apparent utility of such traits and because ecological patterns of distribution of clonal and non-clonal plants differ, adaptation is a tempting explanation of observed clonal life-history variation. However, adaptive evolution requires (1) heritable genetic variation and (2) a trait effect on fitness, and (3) it may be constrained if other evolutionary forces are overriding selection or by constraints, costs and trade-offs. (1) The few studies undertaken so far reported broad-sense heritability for clonal traits. Variation in selectively neutral genetic markers appears as pronounced in populations of clonal as non-clonal plants. However, neutral markers may not reflect heritable variation of life-history traits. Moreover, clonal plants may have been sampled at larger spatial scales. Empirical information on the contribution of somatic mutations to heritable variation is lacking. (2) Clonal life-history traits were found to affect fitness. However, much of this evidence stems from artificial rather than natural environments. (3) The relative importance of gene flow, inbreeding, and genetic drift, compared with selection, in the evolution of clonal life histories is hardly explored. Benefits of clonal life-history traits were frequently studied and found. However, there is also evidence for constraints, trade-offs, and costs. In conclusion, though it is very likely, that clonal life-history traits are adaptive, it is neither clear to which degree this is the case, nor which clonal life-history traits constitute adaptations to which environmental factors. Moreover, evolutionary interactions among clonal life-history traits and between clonal and non-clonal ones, such as the mating system, are not well explored. There remains much interesting work to be done in this field – which will be particularly interesting if it is done in the field.     

Genealogy of Neutral Genes and Spreading of Selected Mutations in a Geographically Structured Population

In a geographically structured population, the interplay among gene migration, genetic drift and natural selection raises intriguing evolutionary problems, but the rigorous mathematical treatment is often very difficult. Therefore several approximate formulas were developed concerning the coalescence process of neutral genes and the fixation process of selected mutations in an island model, and their accuracy was examined by computer simulation. When migration is limited, the coalescence (or divergence) time for sampled neutral genes can be described by the convolution of exponential functions, as in a panmictic population, but it is determined mainly by migration rate and the number of demes from which the sample is taken. This time can be much longer than that in a panmictic population with the same number of breeding individuals. For a selected mutation, the spreading over the entire population was formulated as a birth and death process, in which the fixation probability within a deme plays a key role. With limited amounts of migration, even advantageous mutations take a large number of generations to spread. Furthermore, it is likely that these mutations which are temporarily fixed in some demes may be swamped out again by non-mutant immigrants from other demes unless selection is strong enough. These results are potentially useful for testing quantitatively various hypotheses that have been proposed for the origin of modern human populations.     

Interfering waves of adaptation promote spatial mixing

A fundamental problem of asexual adaptation is that beneficial substitutions are not efficiently accumulated in large populations: Beneficial mutations often go extinct because they compete with one another in going to fixation. It has been argued that such clonal interference may have led to the evolution of sex and recombination in well-mixed populations. Here, we study clonal interference, and mechanisms of its mitigation, in an evolutionary model of spatially structured populations with uniform selection pressure. Clonal interference is much more prevalent with spatial structure than without, due to the slow wave-like spread of beneficial mutations through space. We find that the adaptation speed of asexuals saturates when the linear habitat size exceeds a characteristic interference length, which becomes shorter with smaller migration and larger mutation rate. The limiting speed is proportional to μ(1/2) and μ(1/3) in linear and planar habitats, respectively, where the mutational supply μ is the product of mutation rate and local population density. This scaling and the existence of a speed limit should be amenable to experimental tests as they fall far below predicted adaptation speeds for well-mixed populations (that scale as the logarithm of population size). Finally, we show that not only recombination, but also long-range migration is a highly efficient mechanism of relaxing clonal competition in structured populations. Our conservative estimates of the interference length predict prevalent clonal interference in microbial colonies and biofilms, so clonal competition should be a strong driver of both genetic and spatial mixing in those contexts.     

Regions of lower crossing over harbor more rare variants in African populations of Drosophila melanogaster

A correlation between diversity levels and rates of recombination is predicted both by models of positive selection, such as hitchhiking associated with the rapid fixation of advantageous mutations, and by models of purifying selection against strongly deleterious mutations (commonly referred to as "background selection"). With parameter values appropriate for Drosophila populations, only the first class of models predicts a marked skew in the frequency spectrum of linked neutral variants, relative to a neutral model. Here, we consider 29 loci scattered throughout the Drosophila melanogaster genome. We show that, in African populations, a summary of the frequency spectrum of polymorphic mutations is positively correlated with the meiotic rate of crossing over. This pattern is demonstrated to be unlikely under a model of background selection. Models of weakly deleterious selection are not expected to produce both the observed correlation and the extent to which nucleotide diversity is reduced in regions of low (but nonzero) recombination. Thus, of existing models, hitchhiking due to the recurrent fixation of advantageous variants is the most plausible explanation for the data.     

Estimating the distribution of selection coefficients from phylogenetic data using sitewise mutation-selection models

Estimation of the distribution of selection coefficients of mutations is a long-standing issue in molecular evolution. In addition to population-based methods, the distribution can be estimated from DNA sequence data by phylogenetic-based models. Previous models have generally found unimodal distributions where the probability mass is concentrated between mildly deleterious and nearly neutral mutations. Here we use a sitewise mutation-selection phylogenetic model to estimate the distribution of selection coefficients among novel and fixed mutations (substitutions) in a data set of 244 mammalian mitochondrial genomes and a set of 401 PB2 proteins from influenza. We find a bimodal distribution of selection coefficients for novel mutations in both the mitochondrial data set and for the influenza protein evolving in its natural reservoir, birds. Most of the mutations are strongly deleterious with the rest of the probability mass concentrated around mildly deleterious to neutral mutations. The distribution of the coefficients among substitutions is unimodal and symmetrical around nearly neutral substitutions for both data sets at adaptive equilibrium. About 0.5% of the nonsynonymous mutations and 14% of the nonsynonymous substitutions in the mitochondrial proteins are advantageous, with 0.5% and 24% observed for the influenza protein. Following a host shift of influenza from birds to humans, however, we find among novel mutations in PB2 a trimodal distribution with a small mode of advantageous mutations.     

Natural selection in a postglacial range expansion: the case of the colour cline in the European barn owl

Gradients of variation—or clines—have always intrigued biologists. Classically, they have been interpreted as the outcomes of antagonistic interactions between selection and gene flow. Alternatively, clines may also establish neutrally with isolation by distance (IBD) or secondary contact between previously isolated populations. The relative importance of natural selection and these two neutral processes in the establishment of clinal variation can be tested by comparing genetic differentiation at neutral genetic markers and at the studied trait. A third neutral process, surfing of a newly arisen mutation during the colonization of a new habitat, is more difficult to test. Here, we designed a spatially explicit approximate Bayesian computation (ABC) simulation framework to evaluate whether the strong cline in the genetically based reddish coloration observed in the European barn owl (Tyto alba) arose as a by‐product of a range expansion or whether selection has to be invoked to explain this colour cline, for which we have previously ruled out the actions of IBD or secondary contact. Using ABC simulations and genetic data on 390 individuals from 20 locations genotyped at 22 microsatellites loci, we first determined how barn owls colonized Europe after the last glaciation. Using these results in new simulations on the evolution of the colour phenotype, and assuming various genetic architectures for the colour trait, we demonstrate that the observed colour cline cannot be due to the surfing of a neutral mutation. Taking advantage of spatially explicit ABC, which proved to be a powerful method to disentangle the respective roles of selection and drift in range expansions, we conclude that the formation of the colour cline observed in the barn owl must be due to natural selection.     

Noise driven evolutionary waves

Adaptation in spatially extended populations entails the propagation of evolutionary novelties across habitat ranges. Driven by natural selection, beneficial mutations sweep through the population in a "wave of advance". The standard model for these traveling waves, due to R. Fisher and A. Kolmogorov, plays an important role in many scientific areas besides evolution, such as ecology, epidemiology, chemical kinetics, and recently even in particle physics. Here, we extend the Fisher-Kolmogorov model to account for mutations that confer an increase in the density of the population, for instance as a result of an improved metabolic efficiency. We show that these mutations invade by the action of random genetic drift, even if the mutations are slightly deleterious. The ensuing class of noise-driven waves are characterized by a wave speed that decreases with increasing population sizes, contrary to conventional Fisher-Kolmogorov waves. When a trade-off exists between density and growth rate, an evolutionary optimal population density can be predicted. Our simulations and analytical results show that genetic drift in conjunction with spatial structure promotes the economical use of limited resources. The simplicity of our model, which lacks any complex interactions between individuals, suggests that noise-induced pattern formation may arise in many complex biological systems including evolution.     

Natural Selection Reduced Diversity on Human Y Chromosomes

The human Y chromosome exhibits surprisingly low levels of genetic diversity. This could result from neutral processes if the effective population size of males is reduced relative to females due to a higher variance in the number of offspring from males than from females. Alternatively, selection acting on new mutations, and affecting linked neutral sites, could reduce variability on the Y chromosome. Here, using genome-wide analyses of X, Y, autosomal and mitochondrial DNA, in combination with extensive population genetic simulations, we show that low observed Y chromosome variability is not consistent with a purely neutral model. Instead, we show that models of purifying selection are consistent with observed Y diversity. Further, the number of sites estimated to be under purifying selection greatly exceeds the number of Y-linked coding sites, suggesting the importance of the highly repetitive ampliconic regions. While we show that purifying selection removing deleterious mutations can explain the low diversity on the Y chromosome, we cannot exclude the possibility that positive selection acting on beneficial mutations could have also reduced diversity in linked neutral regions, and may have contributed to lowering human Y chromosome diversity. Because the functional significance of the ampliconic regions is poorly understood, our findings should motivate future research in this area.     

Somatic mutations in organisms with complex life histories.

A theoretical model is developed of the fate of mutations for organisms with such life-history characteristics as indeterminate growth and clonal reproduction. It focuses on how the fate of a particular mutant depends on whether it arises during mitotic cell division (somatic mutation) or during meiotic cell division (meiotic mutation). At gamete production, individuals carrying somatic mutations will produce some proportion of gametes reflecting the original, zygotic genotype and some proportion reflecting genotypes carrying the somatic mutation. Focusing on allele frequencies at gamete production allows the effects of growth and clonal reproduction to be summarized. The relative strengths of somatic and meiotic mutation can be determined, as well as the conditions under which the change in allele frequency due to one is greater than that due to the other. Examples from a published demographic study of clonal corals are used to compare somatic and meiotic mutation. When there is no selection acting on either type of mutation, only a few cell divisions per time unit on average are needed for the change in allele frequency due to somatic mutation to be greater, given empirically based mutation rates. When somatic selection is added, the most dramatic effect is seen with fairly strong negative selection acting against the somatic mutation within individuals. In this case, selection within organisms can effectively counteract the effects of somatic mutation, and the change in allele frequency due to somatic mutations will not be greater than that due to meiotic mutations for reasonable numbers of within-generation cell divisions. The majority of the mutation load, which would have been due to somatic mutation, is purged by selection within the individual organism.     

Nonclinality of molecular variation implicates selection in maintaining a morphological cline of Drosophila melanogaster

One general approach for assessing whether phenotypic variation is due to selection is to test its correlation with presumably neutral molecular variation. Neutral variation is determined by population history, the most likely alternative explanation of spatial genetic structure, whereas phenotypic variation may be influenced by the spatial pattern of selection pressure. Several methods for comparing the spatial apportionment of molecular and morphological variation have been used. Here, we present an analysis of variance framework that compares the magnitudes of latitudinal effects for molecular and morphological variation along a body size cline in Australian Drosophila populations. Explicit incorporation of the relevant environmental gradient can result in a simple and powerful test of selection. For the Australian cline, our analysis provides strong internal evidence that the cline is due to selection.     

Selection against Heteroplasmy Explains the Evolution of Uniparental Inheritance of Mitochondria

Why are mitochondria almost always inherited from one parent during sexual reproduction? Current explanations for this evolutionary mystery include conflict avoidance between the nuclear and mitochondrial genomes, clearing of deleterious mutations, and optimization of mitochondrial-nuclear coadaptation. Mathematical models, however, fail to show that uniparental inheritance can replace biparental inheritance under any existing hypothesis. Recent empirical evidence indicates that mixing two different but normal mitochondrial haplotypes within a cell (heteroplasmy) can cause cell and organism dysfunction. Using a mathematical model, we test if selection against heteroplasmy can lead to the evolution of uniparental inheritance. When we assume selection against heteroplasmy and mutations are neither advantageous nor deleterious (neutral mutations), uniparental inheritance replaces biparental inheritance for all tested parameter values. When heteroplasmy involves mutations that are advantageous or deleterious (non-neutral mutations), uniparental inheritance can still replace biparental inheritance. We show that uniparental inheritance can evolve with or without pre-existing mating types. Finally, we show that selection against heteroplasmy can explain why some organisms deviate from strict uniparental inheritance. Thus, we suggest that selection against heteroplasmy explains the evolution of uniparental inheritance.