Population structure, levels of selection, and the evolution of intracellular symbionts

Many obligately intracellular symbionts exhibit a characteristic set of genetic changes that include an increase in substitution rates, loss of many genes, and apparent destabilization of many proteins and structural RNAs. Authors have suggested that these changes are due to increased mutation rates, or, more commonly, decreased effective population size due to population bottlenecks at the symbiont or, perhaps, host level. I propose that the increase in substitution rates and accumulation of deleterious mutations is a consequence of the population structure imposed on the endosymbionts by strict host association, loss of horizontal transmission and potentially conflicting levels of selection. I analyze a population genetic model of endosymbiont evolution, and demonstrate that substitution rates will increase, and the effect of those substitutions on endosymbiont fitness will become more deleterious as horizontal transmission among hosts decreases. Additionally, I find that there is a critical level of horizontal transmission below which natural selection cannot effectively purge deleterious mutations, leading to an expected loss of fitness over time. This critical level varies across loci with the degree of correlation between host and endosymbiont fitness, and may help explain differential retention and loss of certain genes.     

Fitness decline in spontaneous mutation accumulation lines of Caenorhabditis elegans with varying effective population sizes

The rate and fitness effects of new mutations have been investigated by mutation accumulation (MA) experiments in which organisms are maintained at a constant minimal population size to facilitate the accumulation of mutations with minimal efficacy of selection. We evolved 35 MA lines of Caenorhabditis elegans in parallel for 409 generations at three population sizes (N = 1, 10, and 100), representing the first spontaneous long-term MA experiment at varying population sizes with corresponding differences in the efficacy of selection. Productivity and survivorship in the N = 1 lines declined by 44% and 12%, respectively. The average effects of deleterious mutations in N = 1 lines are estimated to be 16.4% for productivity and 11.8% for survivorship. Larger populations (N = 10 and 100) did not suffer a significant decline in fitness traits despite a lengthy and sustained regime of consecutive bottlenecks exceeding 400 generations. Together, these results suggest that fitness decline in very small populations is dominated by mutations with large deleterious effects. It is possible that the MA lines at larger population sizes contain a load of cryptic deleterious mutations of small to moderate effects that would be revealed in more challenging environments.     

Rapid fitness recovery in mutationally degraded lines of Caenorhabditis elegans

Abstract Deleterious mutation accumulation has been implicated in many biological phenomena and as a potentially significant threat to human health and the persistence of small populations. The vast majority of mutations with effects on fitness are known to be deleterious in a given environment, and their accumulation results in mean population fitness decline. However, whether populations are capable of recovering from negative effects of prolonged genetic bottlenecks via beneficial or compensatory mutation accumulation has not previously been tested. To address this question, long-term mutation-accumulation lines of the nematode Caenorhabditis elegans , previously propagated as single individuals each generation, were maintained in large population sizes under competitive conditions. Fitness assays of these lines and comparison to parallel mutation-accumulation lines and the ancestral control show that, while the process of fitness restoration was incomplete for some lines, full recovery of mean fitness was achieved in fewer than 80 generations. Several lines of evidence indicate that this fitness restoration was at least partially driven by compensatory mutation accumulation rather than a result of a generic form of laboratory adaptation. This surprising result has broad implications for the influence of the mutational process on many issues in evolutionary and conservation biology.     

Conservation of mRNA secondary structures may filter out mutations in Escherichia coli evolution

Recent reports indicate that mutations in viral genomes tend to preserve RNA secondary structure, and those mutations that disrupt secondary structural elements may reduce gene expression levels, thereby serving as a functional knockout. In this article, we explore the conservation of secondary structures of mRNA coding regions, a previously unknown factor in bacterial evolution, by comparing the structural consequences of mutations in essential and nonessential Escherichia coli genes accumulated over 40 000 generations in the course of the ‘long-term evolution experiment’. We monitored the extent to which mutations influence minimum free energy (MFE) values, assuming that a substantial change in MFE is indicative of structural perturbation. Our principal finding is that purifying selection tends to eliminate those mutations in essential genes that lead to greater changes of MFE values and, therefore, may be more disruptive for the corresponding mRNA secondary structures. This effect implies that synonymous mutations disrupting mRNA secondary structures may directly affect the fitness of the organism. These results demonstrate that the need to maintain intact mRNA structures imposes additional evolutionary constraints on bacterial genomes, which go beyond preservation of structure and function of the encoded proteins.     

Experimental genomic evolution: extensive compensation for loss of DNA ligase activity in a virus

Deletion of the viral ligase gene drastically reduced the fitness of bacteriophage T7 on a ligase-deficient host. Viral evolution recovered much of this fitness during long-term passage, but the final fitness remained below that of the intact virus. Compensatory changes occurred chiefly in genes involved in DNA metabolism: the viral endonuclease, helicase, and DNA polymerase. Two other compensatory changes of unknown function also occurred. Using a method to distinguish compensatory mutations from other beneficial mutations, five additional substitutions from the recovery were shown to enhance adaptation to culture conditions and were not compensatory for the deletion. In contrast to the few previous studies of viral recovery from deletions, the compensatory changes in T7 did not restore the deletion or duplicate major regions of the genome. The ability of this deleted genome to recover much of the lost fitness via mutations in its remaining genes reveals a considerable evolutionary potential to modify the interactions of its elements in maintaining an essential set of functions.     

Adaptation and ecological resistance

The notion fitness, widely used in genetics usually serves to measure a relative rate of organism reproduction. Another important character of an organism is its ecological resistance which is basically the product of macroevolution. It can be determined as a probability of an organism survival and participation in reproduction of the species. Ecological resistance determines the level of the accidental death of organisms that are genetically valuable. For the comparison of ecological resistance in different organisms and species the negative meanings of the Malthusian parameter can be used. Ecological resistance depends on the presence in genomes of essential genes and fairly complete sets of nonessential, or adaptive, genes which can reside in genomes both as "plus" and "minus" alleles. The recovery of complete sets of adaptive genes lost as a result of mutations and, thus, of a high level of ecological resistance in organisms is provided by genetic exchange between them. With respect to mutations leading to the increase in fitness the effect of genetic exchange is negative since it leads to the formation of recombination load, i.e. a decrease in fitness of the offspring. In microevolutionary processes, the elevation in ecological resistance level does not take place since it requires a long time for the formation of new genes and new elements of organization in the process of positive selection. At the same time, a constant recovery of a high level of ecological resistance of the species decreased as a result of mutations takes place in some individuals due to genetic exchange. Mutations affecting ecological resistance of an organism, as a rule, cause a decrease in its viability and they are usually excluded from populations as a result of negative (stabilizing) selection.     

Population bottlenecks constrain host microbiome diversity and genetic variation impeding fitness

It is becoming increasingly clear that microbial symbionts influence key aspects of their host’s fitness, and vice versa. This may fundamentally change our thinking about how microbes and hosts interact in influencing fitness and adaptation to changing environments. Here we explore how reductions in population size commonly experienced by threatened species influence microbiome diversity. Consequences of such reductions are normally interpreted in terms of a loss of genetic variation, increased inbreeding and associated inbreeding depression. However, fitness effects of population bottlenecks might also be mediated through microbiome diversity, such as through loss of functionally important microbes. Here we utilise 50 Drosophila melanogaster lines with different histories of population bottlenecks to explore these questions. The lines were phenotyped for egg-to-adult viability and their genomes sequenced to estimate genetic variation. The bacterial 16S rRNA gene was amplified in these lines to investigate microbial diversity. We found that 1) host population bottlenecks constrained microbiome richness and diversity, 2) core microbiomes of hosts with low genetic variation were constituted from subsets of microbiomes found in flies with higher genetic variation, 3) both microbiome diversity and host genetic variation contributed to host population fitness, 4) connectivity and robustness of bacterial networks was low in the inbred lines regardless of host genetic variation, 5) reduced microbial diversity was associated with weaker evolutionary responses of hosts in stressful environments, and 6) these effects were unrelated to Wolbachia density. These findings suggest that population bottlenecks reduce hologenomic variation (combined host and microbial genetic variation). Thus, while the current biodiversity crisis focuses on population sizes and genetic variation of eukaryotes, an additional focal point should be the microbial diversity carried by the eukaryotes, which in turn may influence host fitness and adaptability with consequences for the persistence of populations.     

Purging of inbreeding depression and fitness decline in bottlenecked populations of Drosophila melanogaster

Drastic reductions in population size, or bottlenecks, are thought to significantly erode genetic variability and reduce fitness. However, it has been suggested that a population can be purged of the genetic load responsible for reduced fitness when subjected to bottlenecks. To investigate this phenomenon, we put a number of Drosophila melanogaster isofemale lines known to differ in inbreeding depression through four ‘founder‐flush’ bottleneck cycles with flush sizes of 5 or 100 pairs and assayed for relative fitness (single‐pair productivity) after each cycle. Following the founder‐flush phase, the isofemale lines, with a large flush size and a history of inbreeding depression, recovered most of the fitness lost from early inbreeding, consistent with purging. The same isofemale lines, with a small flush size, did not regain fitness, consistent with the greater effect of genetic drift in these isofemale lines. On the other hand, the isofemale lines that did not show initial inbreeding depression declined in fitness after repeated bottlenecks, independent of the flush size. These results suggest that the nature of genetic variation in fitness may greatly influence the way in which populations respond to bottlenecks and that stochastic processes play an important role. Consequently, an attempt intentionally to purge a population of detrimental variation through inbreeding appears to be a risky strategy, particularly in the genetic management of endangered species.     

Evaluating the within-host fitness effects of mutations fixed during virus adaptation to different ecotypes of a new host

The existence of genetic variation for resistance in host populations is assumed to be essential to the spread of an emerging virus. Models predict that the rate of spread slows down with the increasing frequency and higher diversity of resistance alleles in the host population. We have been using the experimental pathosystem Arabidopsis thaliana—tobacco etch potyvirus (TEV) to explore the interplay between genetic variation in host''s susceptibility and virus diversity. We have recently shown that TEV populations evolving in A. thaliana ecotypes that differ in susceptibility to infection gained within-host fitness, virulence and infectivity in a manner compatible with a gene-for-gene model of host–parasite interactions: hard-to-infect ecotypes were infected by generalist viruses, whereas easy-to-infect ecotypes were infected by every virus. We characterized the genomes of the evolved viruses and found cases of host-driven convergent mutations. To gain further insights in the mechanistic basis of this gene-for-gene model, we have generated all viral mutations individually as well as in specific combinations and tested their within-host fitness effects across ecotypes. Most of these mutations were deleterious or neutral in their local ecotype and only a very reduced number had a host-specific beneficial effect. We conclude that most of the mutations fixed during the evolution experiment were so by drift or by selective sweeps along with the selected driver mutation. In addition, we evaluated the ruggedness of the underlying adaptive fitness landscape and found that mutational effects were mostly multiplicative, with few cases of significant epistasis.     

Mutation accumulation in transfer RNAs: molecular evidence for Muller's ratchet in mitochondrial genomes

The accumulation of deleterious mutations is thought to be a major factor preventing the long-term persistence of obligately asexual lineages relative to their sexual ancestors. This phenomenon is also of potential relevance to sexual species that harbor asexually propagating organelle genomes. A comparative study of the transfer RNA genes in animal mitochondrial and nuclear genomes demonstrates that the former accumulate nucleotide substitutions much more rapidly than do the latter, and several lines of evidence are consistent with the idea that the excess substitutions are mildly deleterious. First, the average binding stability between complementary strands in the stems of mitochondrial tRNAs is less than half that in nuclear tRNAs. Second, most loop sizes in the mitochondrial tRNAs have experienced a net reduction in size over evolutionary time, and they are nearly 50 times more variable in the mitochondrial than in the nuclear genome. Third, although nearly 20% of the nucleotides in nuclear tRNA genes (particularly those involved in tertiary interactions) are invariant across all animal taxa and all tRNA species, there are no invariant sites in the mitochondrial tRNAs. These observations, as well as results from recent laboratory experiments, are consistent with the hypothesis that nonrecombining organelle genomes are subject to gradual loss of fitness due to the cumulative chance fixation of mildly deleterious mutations.      

Fitness Is Strongly Influenced by Rare Mutations of Large Effect in a Microbial Mutation Accumulation Experiment

Our understanding of the evolutionary consequences of mutation relies heavily on estimates of the rate and fitness effect of spontaneous mutations generated by mutation accumulation (MA) experiments. We performed a classic MA experiment in which frequent sampling of MA lines was combined with whole genome resequencing to develop a high-resolution picture of the effect of spontaneous mutations in a hypermutator (ΔmutS) strain of the bacterium Pseudomonas aeruginosa. After ∼644 generations of mutation accumulation, MA lines had accumulated an average of 118 mutations, and we found that average fitness across all lines decayed linearly over time. Detailed analyses of the dynamics of fitness change in individual lines revealed that a large fraction of the total decay in fitness (42.3%) was attributable to the fixation of rare, highly deleterious mutations (comprising only 0.5% of fixed mutations). Furthermore, we found that at least 0.64% of mutations were beneficial and probably fixed due to positive selection. The majority of mutations that fixed (82.4%) were base substitutions and we failed to find any signatures of selection on nonsynonymous or intergenic mutations. Short indels made up a much smaller fraction of the mutations that were fixed (17.4%), but we found evidence of strong selection against indels that caused frameshift mutations in coding regions. These results help to quantify the amount of natural selection present in microbial MA experiments and demonstrate that changes in fitness are strongly influenced by rare mutations of large effect.     

Continuous In Vitro Evolution of Ribozymes That Operate Under Conditions of Extreme pH

Continuous in vitro evolution methods were used to study the behavior of an evolving population of RNA ligase ribozymes in response to selection pressures involving conditions of extreme pH. The starting population consisted of randomized variants of a ribozyme that had been optimized for activity at pH 8.5. The ribozymes were subjected to repeated rounds of selective amplification under progressively more acidic or more alkaline conditions. The two final evolved populations of ribozymes were able to operate at either pH 5.8 or pH 9.8, respectively. Representative individuals from the two final populations were isolated and characterized. The low-pH ribozyme exhibited a 10-fold increase in catalytic rate at pH 5.8 compared to the starting molecule. The high-pH ribozyme retained its structural integrity and activity at pH 9.8, whereas the starting molecule was denatured under this condition. These findings demonstrate that a population of functional macromolecules can adapt to stringent environmental conditions through the acquisition of relatively few mutations. The results establish continuous in vitro evolution as a useful model system for exploring the evolution of enzymatic function in extreme environments. Present address (Henriette Kühne): Cardinal Health, 2950 Trade Place, San Diego, CA 92126, USA     

Impact of epistasis and pleiotropy on evolutionary adaptation

Evolutionary adaptation is often likened to climbing a hill or peak. While this process is simple for fitness landscapes where mutations are independent, the interaction between mutations (epistasis) as well as mutations at loci that affect more than one trait (pleiotropy) are crucial in complex and realistic fitness landscapes. We investigate the impact of epistasis and pleiotropy on adaptive evolution by studying the evolution of a population of asexual haploid organisms (haplotypes) in a model of N interacting loci, where each locus interacts with K other loci. We use a quantitative measure of the magnitude of epistatic interactions between substitutions, and find that it is an increasing function of K. When haplotypes adapt at high mutation rates, more epistatic pairs of substitutions are observed on the line of descent than expected. The highest fitness is attained in landscapes with an intermediate amount of ruggedness that balance the higher fitness potential of interacting genes with their concomitant decreased evolvability. Our findings imply that the synergism between loci that interact epistatically is crucial for evolving genetic modules with high fitness, while too much ruggedness stalls the adaptive process.     

The effect of population bottlenecks on mutation rate evolution in asexual populations

In the absence of recombination, a mutator allele can spread through a population by hitchhiking with beneficial mutations that appear in its genetic background. Theoretical studies over the past decade have shown that the survival and fixation probability of beneficial mutations can be severely reduced by population size bottlenecks. Here, we use computational modelling and evolution experiments with the yeast S. cerevisiae to examine whether population bottlenecks can affect mutator dynamics in adapting asexual populations. In simulation, we show that population bottlenecks can inhibit mutator hitchhiking with beneficial mutations and are most effective at lower beneficial mutation supply rates. We then subjected experimental populations of yeast propagated at the same effective population size to three different bottleneck regimes and observed that the speed of mutator hitchhiking was significantly slower at smaller bottlenecks, consistent with our theoretical expectations. Our results, thus, suggest that bottlenecks can be an important factor in mutation rate evolution and can in certain circumstances act to stabilize or, at least, delay the progressive elevation of mutation rates in asexual populations. Additionally, our findings provide the first experimental support for the theoretically postulated effect of population bottlenecks on beneficial mutations and demonstrate the usefulness of studying mutator frequency dynamics for understanding the underlying dynamics of fitness‐affecting mutations.     

Survival Probability of Beneficial Mutations in Bacterial Batch Culture

The survival of rare beneficial mutations can be extremely sensitive to the organism’s life history and the trait affected by the mutation. Given the tremendous impact of bacteria in batch culture as a model system for the study of adaptation, it is important to understand the survival probability of beneficial mutations in these populations. Here we develop a life-history model for bacterial populations in batch culture and predict the survival of mutations that increase fitness through their effects on specific traits: lag time, fission time, viability, and the timing of stationary phase. We find that if beneficial mutations are present in the founding population at the beginning of culture growth, mutations that reduce the mortality of daughter cells are the most likely to survive drift. In contrast, of mutations that occur de novo during growth, those that delay the onset of stationary phase are the most likely to survive. Our model predicts that approximately fivefold population growth between bottlenecks will optimize the occurrence and survival of beneficial mutations of all four types. This prediction is relatively insensitive to other model parameters, such as the lag time, fission time, or mortality rate of the population. We further estimate that bottlenecks that are more severe than this optimal prediction substantially reduce the occurrence and survival of adaptive mutations.     

THE DISTRIBUTION OF PHENOTYPIC VARIANCE WITH INBREEDING

Fifty-two inbred populations of Drosophila melanogaster, each founded from a single pair, and a large number of control, outbred flies were measured for fitness and a set of six traits. A survey of the literature on the effects of inbreeding and population bottlenecks demonstrates that the commonly observed pattern of an apparent variance among characters and among species in changes of phenotypic variance may in fact be largely the result of sampling error, given the pattern of change that we demonstrate within a species for the same character. In our study, population bottlenecks on average decrease the amount of phenotypic variance for a suite of wing characteristics and size, but there is large and significant variation among lines in the amount of phenotypic variance. As a result, several lines increased in variance in spite of the average decrease. Interestingly, the changes in phenotypic variance for fitness are in sharp contrast to those seen for phenotypic variance for morphological traits. The amount of phenotypic variance for fitness varies highly significantly among lines but, on average, is increased by bottlenecks. The changes in phenotypic variance as a result of population bottlenecks are large enough to significantly affect the probability of peak shifts by the variance-induced peak shift model.     

The impact of high‐order epistasis in the within‐host fitness of a positive‐sense plant RNA virus

RNA viruses are the main source of emerging infectious diseases because of the evolutionary potential bestowed by their fast replication, large population sizes and high mutation and recombination rates. However, an equally important property, which is usually neglected, is the topography of the fitness landscape. How many fitness maxima exist and how well they are connected is especially interesting, as this determines the number of accessible evolutionary pathways. To address this question, we have reconstructed a region of the fitness landscape of tobacco etch potyvirus constituted by mutations observed during the experimental adaptation of the virus to the novel host Arabidopsis thaliana. Fitness was measured for many genotypes and showed the existence of multiple peaks and holes in the landscape. We found prevailing epistatic effects between mutations, with cases of reciprocal sign epistasis being common among pairs of mutations. We also found that high‐order epistasis was as important as pairwise epistasis in their contribution to fitness. Therefore, results suggest that the landscape was rugged due to the existence of holes caused by lethal genotypes, that a very limited number of potential neutral paths exist and that it contained a single adaptive peak.     

Mutational neighbourhood and mutation supply rate constrain adaptation in Pseudomonas aeruginosa

Understanding adaptation by natural selection requires understanding the genetic factors that determine which beneficial mutations are available for selection. Here, using experimental evolution of rifampicin-resistant Pseudomonas aeruginosa, we show that different genotypes vary in their capacity for adaptation to the cost of antibiotic resistance. We then use sequence data to show that the beneficial mutations associated with fitness recovery were specific to particular genetic backgrounds, suggesting that genotypes had access to different sets of beneficial mutations. When we manipulated the supply rate of beneficial mutations, by altering effective population size during evolution, we found that it constrained adaptation in some selection lines by restricting access to rare beneficial mutations, but that the effect varied among the genotypes in our experiment. These results suggest that mutational neighbourhood varies even among genotypes that differ by a single amino acid change, and this determines their capacity for adaptation as well as the influence of population biology processes that alter mutation supply rate.     

Molecular basis of fitness loss and fitness recovery in vesicular stomatitis virus

Viral populations subjected to repeated genetic bottleneck accumulate deleterious mutations in a process known as Muller's ratchet. Asexual viruses, such as vesicular stomatitis virus (VSV) can recover from Muller's ratchet by replication with large effective population sizes. However, mutants with a history of bottleneck transmissions often show decreased adaptability when compared to non-bottlenecked populations. We have generated a collection of bottlenecked mutants and allowed them to recover by large population passages. We have characterized fitness changes and the complete genomes of these strains. Mutations accumulated during the operation of Muller's ratchet led to the identification of two potential mutational hot spots in the VSV genome. As in other viral systems, transitions were more common than transversions. Both back mutation and compensatory mutations contributed to recovery, although a significant level of fitness increase was observed in nine of the 13 bottlenecked strains with no obvious changes in the consensus sequence. Additional replication of three strains resulted in the fixation of single point mutations. Only two mutations previously found in non-bottlenecked, high-fitness populations that had been adapting to the same environment were identified in the recovered strains.