Experimental Evolution of Enhanced Growth by Bacillus subtilis at Low Atmospheric Pressure: Genomic Changes Revealed by Whole-Genome Sequencing

Knowledge of how microorganisms respond and adapt to low-pressure (LP) environments is limited. Previously, Bacillus subtilis strain WN624 was grown at the near-inhibitory LP of 5 kPa for 1,000 generations and strain WN1106, which exhibited increased relative fitness at 5 kPa, was isolated. Genomic sequence differences between ancestral strain WN624 and LP-evolved strain WN1106 were identified using whole-genome sequencing. LP-evolved strain WN1106 carried amino acid-altering mutations in the coding sequences of only seven genes (fliI, parC, ytoI, bacD, resD, walK, and yvlD) and a single 9-nucleotide in-frame deletion in the rnjB gene that encodes RNase J2, a component of the RNA degradosome. By using a collection of frozen stocks of the LP-evolved culture taken at 50-generation intervals, it was determined that (i) the fitness increase at LP occurred rapidly, while (ii) mutation acquisition exhibited complex kinetics. A knockout mutant of rnjB was shown to increase the competitive fitness of B. subtilis at both LP and standard atmospheric pressure.     

Properties of selected mutations and genotypic landscapes under Fisher's geometric model

The fitness landscape—the mapping between genotypes and fitness—determines properties of the process of adaptation. Several small genotypic fitness landscapes have recently been built by selecting a handful of beneficial mutations and measuring fitness of all combinations of these mutations. Here, we generate several testable predictions for the properties of these small genotypic landscapes under Fisher's geometric model of adaptation. When the ancestral strain is far from the fitness optimum, we analytically compute the fitness effect of selected mutations and their epistatic interactions. Epistasis may be negative or positive on average depending on the distance of the ancestral genotype to the optimum and whether mutations were independently selected, or coselected in an adaptive walk. Simulations show that genotypic landscapes built from Fisher's model are very close to an additive landscape when the ancestral strain is far from the optimum. However, when it is close to the optimum, a large diversity of landscape with substantial roughness and sign epistasis emerged. Strikingly, small genotypic landscapes built from several replicate adaptive walks on the same underlying landscape were highly variable, suggesting that several realizations of small genotypic landscapes are needed to gain information about the underlying architecture of the fitness landscape.     

Fitness change in relation to mutation number in spontaneous mutation accumulation lines of Chlamydomonas reinhardtii

Although all genetic variation ultimately stems from mutations, their properties are difficult to study directly. Here, we used multiple mutation accumulation (MA) lines derived from five genetic backgrounds of the green algae Chlamydomonas reinhardtii that have been previously subjected to whole genome sequencing to investigate the relationship between the number of spontaneous mutations and change in fitness from a nonevolved ancestor. MA lines were on average less fit than their ancestors and we detected a significantly negative correlation between the change in fitness and the total number of accumulated mutations in the genome. Likewise, the number of mutations located within coding regions significantly and negatively impacted MA line fitness. We used the fitness data to parameterize a maximum likelihood model to estimate discrete categories of mutational effects, and found that models containing one to two mutational effect categories (one neutral and one deleterious category) fitted the data best. However, the best‐fitting mutational effects models were highly dependent on the genetic background of the ancestral strain.     

ChimeriVax-West Nile virus live-attenuated vaccine: preclinical evaluation of safety, immunogenicity, and efficacy

The availability of ChimeriVax vaccine technology for delivery of flavivirus protective antigens at the time West Nile (WN) virus was first detected in North America in 1999 contributed to the rapid development of the vaccine candidate against WN virus described here. ChimeriVax-Japanese encephalitis (JE), the first live- attenuated vaccine developed with this technology has successfully undergone phase I and II clinical trials. The ChimeriVax technology utilizes yellow fever virus (YF) 17D vaccine strain capsid and nonstructural genes to deliver the envelope gene of other flaviviruses as live-attenuated chimeric viruses. Amino acid sequence homology between the envelope protein (E) of JE and WN viruses facilitated targeting attenuating mutation sites to develop the WN vaccine. Here we discuss preclinical studies with the ChimeriVax-WN virus in mice and macaques. ChimeriVax-WN virus vaccine is less neurovirulent than the commercial YF 17D vaccine in mice and nonhuman primates. Attenuation of the virus is determined by the chimeric nature of the construct containing attenuating mutations in the YF 17D virus backbone and three point mutations introduced to alter residues 107, 316, and 440 in the WN virus E protein gene. The safety, immunogenicity, and efficacy of the ChimeriVax-WN(02) vaccine in the macaque model indicate the vaccine candidate is expected to be safe and immunogenic for humans.     

AbrB is a regulator of the sigma(W) regulon in Bacillus subtilis

The Bacillus subtilis global regulator AbrB was found to negatively control expression of sigW and genes of the sigma(W) regulon. AbrB bound to DNA regions in the autoregulatory sigW promoter and to some, but not all, of the other sigma(W)-dependent promoters in B. subtilis. Defects in antibiotic resistance properties caused by spo0A mutations are at least partially correlated with AbrB repression of the sigma(W) regulon.     

Mutation-selection balance: ancestry,load, and maximum principle

We analyze the equilibrium behavior of deterministic haploid mutation-selection models. To this end, both the forward and the time-reversed evolution processes are considered. The stationary state of the latter is called the ancestral distribution, which turns out as a key for the study of mutation-selection balance. We find that the ancestral genotype frequencies determine the sensitivity of the equilibrium mean fitness to changes in the corresponding fitness values and discuss implications for the evolution of mutational robustness. We further show that the difference between the ancestral and the population mean fitness, termed mutational loss, provides a measure for the sensitivity of the equilibrium mean fitness to changes in the mutation rate. The interrelation of the loss and the mutation load is discussed. For a class of models in which the number of mutations in an individual is taken as the trait value, and fitness is a function of the trait, we use the ancestor formulation to derive a simple maximum principle, from which the mean and variance of fitness and the trait may be derived; the results are exact for a number of limiting cases, and otherwise yield approximations which are accurate for a wide range of parameters. These results are applied to threshold phenomena caused by the interplay of selection and mutation (known as error thresholds). They lead to a clarification of concepts, as well as criteria for the existence of error thresholds.     

Epistatic interactions between ancestral genotype and beneficial mutations shape evolvability in Pseudomonas aeruginosa

The idea that interactions between mutations influence adaptation by driving populations to low and high fitness peaks on adaptive landscapes is deeply ingrained in evolutionary theory. Here, we investigate the impact of epistasis on evolvability by challenging populations of two Pseudomonas aeruginosa clones bearing different initial mutations (in rpoB conferring rifampicin resistance, and the type IV pili gene network) to adaptation to a medium containing l ‐serine as the sole carbon source. Despite being initially indistinguishable in fitness, populations founded by the two ancestral genotypes reached different fitness following 300 generations of evolution. Genome sequencing revealed that the difference could not be explained by acquiring mutations in different targets of selection; the majority of clones from both ancestors converged on one of the following two strategies: (1) acquiring mutations in either PA2449 (gcsR, an l ‐serine‐metabolism RpoN enhancer binding protein) or (2) protease genes. Additionally, populations from both ancestors converged on loss‐of‐function mutations in the type IV pili gene network, either due to ancestral or acquired mutations. No compensatory or reversion mutations were observed in RNA polymerase (RNAP) genes, in spite of the large fitness costs typically associated with mutations in rpoB. Although current theory points to sign epistasis as the dominant constraint on evolvability, these results suggest that the role of magnitude epistasis in constraining evolvability may be underappreciated. The contribution of magnitude epistasis is likely to be greatest under the biologically relevant mutation supply rates that make back mutations probabilistically unlikely.     

Mutation accumulation in populations of varying size: the distribution of mutational effects for fitness correlates in Caenorhabditis elegans

The consequences of mutation for population-genetic and evolutionary processes depend on the rate and, especially, the frequency distribution of mutational effects on fitness. We sought to approximate the form of the distribution of mutational effects by conducting divergence experiments in which lines of a DNA repair-deficient strain of Caenorhabditis elegans, msh-2, were maintained at a range of population sizes. Assays of these lines conducted in parallel with the ancestral control suggest that the mutational variance is dominated by contributions from highly detrimental mutations. This was evidenced by the ability of all but the smallest population-size treatments to maintain relatively high levels of mean fitness even under the 100-fold increase in mutational pressure caused by knocking out the msh-2 gene. However, we show that the mean fitness decline experienced by larger populations is actually greater than expected on the basis of our estimates of mutational parameters, which could be consistent with the existence of a common class of mutations with small individual effects. Further, comparison of the total mutation rate estimated from direct sequencing of DNA to that detected from phenotypic analyses implies the existence of a large class of evolutionarily relevant mutations with no measurable effect on laboratory fitness.     

Reducing the cost of resistance; experimental evolution in the filamentous fungus Aspergillus nidulans

We have studied compensatory evolution in a fludioxonil resistant mutant of the filamentous fungus Aspergillus nidulans. In an evolution experiment lasting for 27 weeks (about 3000 cell cycles) 35 parallel strains of this mutant evolved in three different environmental conditions. Our results show a severe cost of resistance (56%) in the absence of fludioxonil and in all conditions the mutant strain was able to restore fitness without loss of the resistance. In several cases, the evolved strain reached a higher fitness than the original sensitive ancestor. Fitness compensation occurred in one, two or three discrete steps. Genetic analysis of crosses between different evolved strains and between evolved and ancestral strains revealed interaction between compensatory mutations and provided information on the number of loci involved in fitness compensation. In addition, we discuss the opportunities for the experimental study of evolutionary processes provided by the filamentous fungus A. nidulans.     

Estimating numbers of EMS-induced mutations affecting life history traits in Caenorhabditis elegans in crosses between inbred sublines

Inbred lines of the nematode Caenorhabditis elegans containing independent EMS-induced mutations were crossed to the ancestral wild-type strain (N2). Replicated inbred sublines were generated from the F1 offspring under conditions of minimal selection and, along with the N2 and mutant progenitor lines, were assayed for several fitness correlates including relative fitness (w). A modification of the Castle-Wright estimator and a maximum-likelihood (ML) method were used to estimate the numbers and effects of detectable mutations affecting these characters. The ML method allows for variation in mutational effects by fitting either one or two classes of mutational effect, and uses a Box-Cox power transformation of residual values to account for a skewed distribution of residuals. Both the Castle-Wright and the ML analyses suggest that most of the variation among sublines was due to a few (approximately 1.5-2.5 on average) large-effect mutations. Under ML, a model with two classes of mutational effects, including a class with small effects, fitted better than a single mutation class model, although not significantly better. Nonetheless, given that we expect there to be many mutations induced per line, our results support the hypothesis that mutations vary widely in their effects.     

Measuring selection coefficients below 10(-3): method, questions, and prospects

Measuring fitness with precision is a key issue in evolutionary biology, particularly in studying mutations of small effects. It is usually thought that sampling error and drift prevent precise measurement of very small fitness effects. We circumvented these limits by using a new combined approach to measuring and analyzing fitness. We estimated the mutational fitness effect (MFE) of three independent mini-Tn10 transposon insertion mutations by conducting competition experiments in large populations of Escherichia coli under controlled laboratory conditions. Using flow cytometry to assess genotype frequencies from very large samples alleviated the problem of sampling error, while the effect of drift was controlled by using large populations and massive replication of fitness measures. Furthermore, with a set of four competition experiments between ancestral and mutant genotypes, we were able to decompose fitness measures into four estimated parameters that account for fitness effects of our fluorescent marker (α), the mutation (β), epistasis between the mutation and the marker (γ), and departure from transitivity (τ). Our method allowed us to estimate mean selection coefficients to a precision of 2 × 10(-4). We also found small, but significant, epistatic interactions between the allelic effects of mutations and markers and confirmed that fitness effects were transitive in most cases. Unexpectedly, we also detected variation in measures of s that were significantly bigger than expected due to drift alone, indicating the existence of cryptic variation, even in fully controlled experiments. Overall our results indicate that selection coefficients are best understood as being distributed, representing a limit on the precision with which selection can be measured, even under controlled laboratory conditions.     

Dominance and overdominance of mildly deleterious induced mutations for fitness traits in Caenorhabditis elegans

We estimated the average dominance coefficient of mildly deleterious mutations (h, the proportion by which mutations in the heterozygous state reduce fitness components relative to those in the homozygous state) in the nematode Caenorhabditis elegans. From 56 worm lines that carry mutations induced by the point mutagen ethyl methanesulfonate (EMS), we selected 19 lines that are relatively high in fitness and estimated the viabilities, productivities, and relative fitnesses of heterozygotes and homozygotes compared to the ancestral wild type. There was very little effect of homozygous or heterozygous mutations on egg-to-adult viability. For productivity and relative fitness, we found that the average dominance coefficient, h, was approximately 0.1, suggesting that mildly deleterious mutations are on average partially recessive. These estimates were not significantly different from zero (complete recessivity) but were significantly different from 0.5 (additivity). In addition, there was a significant amount of variation in h among lines, and analysis of average dominance coefficients of individual lines suggested that several lines showed overdominance for fitness. Further investigation of two of these lines partially confirmed this finding.     

Fitness decline under osmotic stress in Caenorhabditis elegans populations subjected to spontaneous mutation accumulation at varying population sizes

The consequences of mutations for population fitness depends on their individual selection coefficients and the effective population size. An earlier study of Caenorhabditis elegans spontaneous mutation accumulation lines evolved for 409 generations at three population sizes found that N_e  = 1 populations declined significantly in fitness whereas the fitness of larger populations (N_e  = 5, 50) was indistinguishable from the ancestral control under benign conditions. To test if larger MA populations harbor a load of cryptic deleterious mutations that are obscured under benign laboratory conditions, we measured fitness under osmotic stress via exposure to hypersaline conditions. The fitness of N_e  = 1 lines exhibited a further decline under osmotic stress compared to benign conditions. However, the fitness of larger populations remained indistinguishable from that of the ancestral control. The average effects of deleterious mutations in N_e  = 1 lines were estimated to be 22% for productivity and 14% for survivorship, exceeding values previously detected under benign conditions. Our results suggest that fitness decline is due to large effect mutations that are rapidly removed via selection even in small populations, with implications for conservation practices. Genetic stochasticity may not be as potent and immediate a threat to the persistence of small populations as other demographic and environmental stochastic factors.     

The experimental evolution of herbicide resistance in Chlamydomonas reinhardtii results in a positive correlation between fitness in the presence and absence of herbicides

Pleiotropic fitness trade-offs will be key determinants of the evolutionary dynamics of selection for pesticide resistance. However, for herbicide resistance, empirical support for a fitness cost of resistance is mixed, and it is therefore also questionable what further ecological trade-offs can be assumed to apply to herbicide resistance. Here, we test the existence of trade-offs by experimentally evolving herbicide resistance in Chlamydomonas reinhardtii. Although fitness costs are detected for all herbicides, we find that, counterintuitively, the most resistant populations also have the lowest fitness costs as measured by growth rate in the ancestral environment. Furthermore, after controlling for differences in the evolutionary dynamics of resistance to different herbicides, we also detect significant positive correlations between resistance, fitness in the ancestral environment and cross-resistance to other herbicides. We attribute this to the highest levels of nontarget-site resistance being achieved by fixing mutations that more broadly affect cellular physiology, which results in both more cross-resistance and less overall antagonistic pleiotropy on maximum growth rate. Consequently, the lack of classical ecological trade-offs could present a major challenge for herbicide resistance management.     

High-resolution mutation mapping reveals parallel experimental evolution in yeast

Understanding the genetic basis of evolutionary adaptation is limited by our ability to efficiently identify the genomic locations of adaptive mutations. Here we describe a method that can quickly and precisely map the genetic basis of naturally and experimentally evolved complex traits using linkage analysis. A yeast strain that expresses the evolved trait is crossed to a distinct strain background and DNA from a large pool of progeny that express the trait of interest is hybridized to oligonucleotide microarrays that detect thousands of polymorphisms between the two strains. Adaptive mutations are detected by linkage to the polymorphisms from the evolved parent. We successfully tested our method by mapping five known genes to a precision of 0.2–24 kb (0.1–10 cM), and developed computer simulations to test the effect of different factors on mapping precision. We then applied this method to four yeast strains that had independently adapted to a fluctuating glucose–galactose environment. All four strains had acquired one or more missense mutations in GAL80, the repressor of the galactose utilization pathway. When transferred into the ancestral strain, the gal80 mutations conferred the fitness advantage that the evolved strains show in the transition from glucose to galactose. Our results show an example of parallel adaptation caused by mutations in the same gene.