The mutational structure of metabolism in Caenorhabditis elegans

A properly functioning organism must maintain metabolic homeostasis. Deleterious mutations degrade organismal function, presumably at least in part via effects on metabolic function. Here we present an initial investigation into the mutational structure of the Caenorhabditis elegans metabolome by means of a mutation accumulation experiment. We find that pool sizes of 29 metabolites vary greatly in their vulnerability to mutation, both in terms of the rate of accumulation of genetic variance (the mutational variance, VM) and the rate of change of the trait mean (the mutational bias, ΔM). Strikingly, some metabolites are much more vulnerable to mutation than any other trait previously studied in the same way. Although we cannot statistically assess the strength of mutational correlations between individual metabolites, principal component analysis provides strong evidence that some metabolite pools are genetically correlated, but also that there is substantial scope for independent evolution of different groups of metabolites. Averaged over mutation accumulation lines, PC3 is positively correlated with relative fitness, but a model in which metabolites are uncorrelated with fitness is nearly as good by Akaike's Information Criterion.     

RAPID DECLINE IN FITNESS OF MUTATION ACCUMULATION LINES OF GONOCHORISTIC (OUTCROSSING) CAENORHABDITIS NEMATODES

Evolutionary theory predicts that the strength of natural selection to reduce the mutation rate should be stronger in self‐fertilizing than in outcrossing taxa. However, the relative efficacy of selection on mutation rate relative to the many other factors influencing the evolution of any species is poorly understood. To address this question, we allowed mutations to accumulate for ∼100 generations in several sets of “mutation accumulation” (MA) lines in three species of gonochoristic (dieocious) Caenorhabditis (C. remanei, C. brenneri, C. sp. 5) as well as in a dioecious strain of the historically self‐fertile hermaprohodite C. elegans. In every case, the rate of mutational decay is substantially greater in the gonochoristic taxa than in C. elegans (∼4× greater on average). Residual heterozygosity in the ancestral controls of these MA lines introduces some complications in interpreting the results, but circumstantial evidence suggests the results are not primarily due to inbreeding depression resulting from residual segregating variation. The results suggest that natural selection operates to optimize the mutation rate in Caenorhabditis and that the strength (or efficiency) of selection differs consistently on the basis of mating system, as predicted by theory. However, context‐dependent environmental and/or synergistic epistasis could also explain the results.     

SPONTANEOUS MUTATION ACCUMULATION IN MULTIPLE STRAINS OF THE GREEN ALGA,CHLAMYDOMONAS REINHARDTII

Estimates of mutational parameters, such as the average fitness effect of a new mutation and the rate at which new genetic variation for fitness is created by mutation, are important for the understanding of many biological processes. However, the causes of interspecific variation in mutational parameters and the extent to which they vary within species remain largely unknown. We maintained multiple strains of the unicellular eukaryote Chlamydomonas reinhardtii, for approximately 1000 generations under relaxed selection by transferring a single cell every ～10 generations. Mean fitness of the lines tended to decline with generations of mutation accumulation whereas mutational variance increased. We did not find any evidence for differences among strains in any of the mutational parameters estimated. The overall change in mean fitness per cell division and rate of input of mutational variance per cell division were more similar to values observed in multicellular organisms than to those in other single‐celled microbes. However, after taking into account differences in genome size among species, estimates from multicellular organisms and microbes, including our new estimates from C. reinhardtii, become substantially more similar. Thus, we suggest that variation in genome size is an important determinant of interspecific variation in mutational parameters.     

Short-term heritable variation overwhelms 200 generations of mutational variance for metabolic traits in Caenorhabditis elegans

Metabolic disorders have a large heritable component, and have increased markedly in human populations over the past few generations. Genome-wide association studies of metabolic traits typically find a substantial unexplained fraction of total heritability, suggesting an important role of spontaneous mutation. An alternative explanation is that epigenetic effects contribute significantly to the heritable variation. Here, we report a study designed to quantify the cumulative effects of spontaneous mutation on adenosine metabolism in the nematode Caenorhabditis elegans, including both the activity and concentration of two metabolic enzymes and the standing pools of their associated metabolites. The only prior studies on the effects of mutation on metabolic enzyme activity, in Drosophila melanogaster, found that total enzyme activity presents a mutational target similar to that of morphological and life-history traits. However, those studies were not designed to account for short-term heritable effects. We find that the short-term heritable variance for most traits is of similar magnitude as the variance among MA lines. This result suggests that the potential heritable effects of epigenetic variation in metabolic disease warrant additional scrutiny.     

Quantifying natural seasonal variation in mutation parameters with mutation accumulation lines

Mutations create novel genetic variants, but their contribution to variation in fitness and other phenotypes may depend on environmental conditions. Furthermore, natural environments may be highly heterogeneous. We assessed phenotypes associated with survival and reproductive success in over 30,000 plants representing 100 mutation accumulation lines of Arabidopsis thaliana across four temporal environments at a single field site. In each of the four assays, environmental variance was substantially larger than mutational variance. For some traits, whether mutational variance was significantly varied between seasons. The founder genotype had mean trait values near the mean of the distribution of the mutation accumulation lines in all field experiments. New mutations also contributed more phenotypic variation than would be predicted, given phenotypic and sequence‐level divergence among natural populations of A. thaliana. The combination of large environmental variance with a mean effect of mutation near zero suggests that mutations could contribute substantially to standing genetic variation.     

MUTATION,SELECTION, AND THE MAINTENANCE OF LIFE-HISTORY VARIATION IN A NATURAL POPULATION

In an effort to provide insight into the role of mutation in the maintenance of genetic variance for life-history traits, we accumulated spontaneous mutations in 10 sets of clonal replicates of Daphnia pulex for approximately 30 generations and compared the variance generated by mutation with the standing level of variation in the wild population. Mutations for quantitative traits appear to arise at a fairly high rate in this species, on the order of at least 0.6 per character per generation, but have relatively small heterozygous effects, changing the phenotype by less than 2.5% of the mean. The mean persistence time of a new mutation affecting life-history/body-size traits is approximately 40 generations in the natural population, which requires an average selection coefficient against new mutations of approximately 3% in the heterozygous state. These data are consistent with the idea that the vast majority of standing genetic variance for life-history characters may be largely a consequence of the recurrent introduction of transient cohorts of mutations that are at least conditionally deleterious and raise issues about the meaning of conventional measures of standing levels of variation for fitness-related traits.     

Effect of migration and environmental heterogeneity on the maintenance of quantitative genetic variation: a simulation study

The paradox of high genetic variation observed in traits under stabilizing selection is a long‐standing problem in evolutionary theory, as mutation rates appear too low to explain observed levels of standing genetic variation under classic models of mutation–selection balance. Spatially or temporally heterogeneous environments can maintain more standing genetic variation within populations than homogeneous environments, but it is unclear whether such conditions can resolve the above discrepancy between theory and observation. Here, we use individual‐based simulations to explore the effect of various types of environmental heterogeneity on the maintenance of genetic variation (V_A) for a quantitative trait under stabilizing selection. We find that V_A is maximized at intermediate migration rates in spatially heterogeneous environments and that the observed patterns are robust to changes in population size. Spatial environmental heterogeneity increased variation by as much as 10‐fold over mutation–selection balance alone, whereas pure temporal environmental heterogeneity increased variance by only 45% at max. Our results show that some combinations of spatial heterogeneity and migration can maintain considerably more variation than mutation–selection balance, potentially reconciling the discrepancy between theoretical predictions and empirical observations. However, given the narrow regions of parameter space required for this effect, this is unlikely to provide a general explanation for the maintenance of variation. Nonetheless, our results suggest that habitat fragmentation may affect the maintenance of V_A and thereby reduce the adaptive capacity of populations.     

PERSPECTIVE: SPONTANEOUS DELETERIOUS MUTATION

Mildly deleterious mutation has been invoked as a leading explanation for a diverse array of observations in evolutionary genetics and molecular evolution and is thought to be a significant risk of extinction for small populations. However, much of the empirical evidence for the deleterious-mutation process derives from studies of Drosophila melanogaster, some of which have been called into question. We review a broad array of data that collectively support the hypothesis that deleterious mutations arise in flies at rate of about one per individual per generation, with the average mutation decreasing fitness by about only 2% in the heterozygous state. Empirical evidence from microbes, plants, and several other animal species provide further support for the idea that most mutations have only mildly deleterious effects on fitness, and several other species appear to have genomic mutation rates that are of the order of magnitude observed in Drosophila. However, there is mounting evidence that some organisms have genomic deleterious mutation rates that are substantially lower than one per individual per generation. These lower rates may be at least partially reconciled with the Drosophila data by taking into consideration the number of germline cell divisions per generation. To fully resolve the existing controversy over the properties of spontaneous mutations, a number of issues need to be clarified. These include the form of the distribution of mutational effects and the extent to which this is modified by the environmental and genetic background and the contribution of basic biological features such as generation length and genome size to interspecific differences in the genomic mutation rate. Once such information is available, it should be possible to make a refined statement about the long-term impact of mutation on the genetic integrity of human populations subject to relaxed selection resulting from modern medical procedures.     

Pseudomonas aeruginosa in a hydropathic facility: diversity, susceptibility and imipenem resistance mutation

Aims: To detect Pseudomonas aeruginosa in water and treatment equipment biofilms of a thermae hydropathic facility and to study antibiotic susceptibility and genetic diversity. Methods and Results: One hundred and fifty‐four planktonic isolates were obtained from 2220 water samples during 4 years. Seventy‐two biofilm isolates were obtained from 23 samples of inner parts of three inhalation equipments. Antibiotic susceptibility was determined by disc diffusion. All isolates were susceptible to tested antimicrobials, except two biofilm isolates and one planktonic isolate. Twenty‐one resistant mutants were observed (nine from biofilms), mostly with imipenem (IP) resistance (81%), by diminished expression of OprD porin, as it was observed by sodium dodecyl sulfate–polyacrylamide gel electrophoresis (SDS–PAGE). Random amplification polymorphic DNA showed a genetically heterogeneous population that is spread through the entire system and persistent in time. IP resistance mutation ability was spread through the population. Conclusions: The permanent assessment of Ps. aeruginosa is necessary not only in water, as expressed in official programmes, but also in equipments where biofilms are evident. Ps. aeruginosa was more prevalent in biofilm populations and presented higher ability to adapt to antibiotic pressure. Significance and Impact of the Study: Twenty‐one million people use thermae in Europe. Official microbiological quality control programmes only consider water surveillance. Present study proves the need of a review on current official programmes.     

Enterohaemorrhagic Escherichia coli O157:H7 Shiga‐like toxin 1 is required for full pathogenicity and activation of the p38 mitogen‐activated protein kinase pathway in Caenorhabditis elegans

Enterohaemorrhagic Escherichia coli (EHEC) causes life‐threatening infections in humans as a consequence of the production of Shiga‐like toxins. Lack of a good animal model system currently hinders in vivo study of EHEC virulence by systematic genetic methods. Here we applied the genetically tractable animal, Caenorhabditis elegans, as a surrogate host to study the virulence of EHEC as well as the host immunity to this human pathogen. Our results show that E. coli O157:H7, a serotype of EHEC, infects and kills C. elegans. Bacterial colonization and induction of the characteristic attaching and effacing (A/E) lesions in the intact intestinal epithelium of C. elegans by E. coli O157:H7 were concomitantly demonstrated in vivo. Genetic analysis indicated that the Shiga‐like toxin 1 (Stx1) of E. coli O157:H7 is a virulence factor in C. elegans and is required for full toxicity. Moreover, the C. elegans p38 mitogen‐activated protein kinase (MAPK) pathway, anevolutionarily conserved innate immune and stress response signalling pathway, is activated in the regulation of host susceptibility to EHEC infection in a Stx1‐dependent manner. Our results validate the EHEC–C. elegans interaction as suitable for future comprehensive genetic screens for both novel bacterial and host factors involved in the pathogenesis of EHEC infection.     

Adaptive introgression in animals: examples and comparison to new mutation and standing variation as sources of adaptive variation

Adaptive genetic variation has been thought to originate primarily from either new mutation or standing variation. Another potential source of adaptive variation is adaptive variants from other (donor) species that are introgressed into the (recipient) species, termed adaptive introgression. Here, the various attributes of these three potential sources of adaptive variation are compared. For example, the rate of adaptive change is generally thought to be faster from standing variation, slower from mutation and potentially intermediate from adaptive introgression. Additionally, the higher initial frequency of adaptive variation from standing variation and lower initial frequency from mutation might result in a higher probability of fixation of the adaptive variants for standing variation. Adaptive variation from introgression might have higher initial frequency than new adaptive mutations but lower than that from standing variation, again making the impact of adaptive introgression variation potentially intermediate. Adaptive introgressive variants might have multiple changes within a gene and affect multiple loci, an advantage also potentially found for adaptive standing variation but not for new adaptive mutants. The processes that might produce a common variant in two taxa, convergence, trans‐species polymorphism from incomplete lineage sorting or from balancing selection and adaptive introgression, are also compared. Finally, potential examples of adaptive introgression in animals, including balancing selection for multiple alleles for major histocompatibility complex (MHC), S and csd genes, pesticide resistance in mice, black colour in wolves and white colour in coyotes, Neanderthal or Denisovan ancestry in humans, mimicry genes in Heliconius butterflies, beak traits in Darwin's finches, yellow skin in chickens and non‐native ancestry in an endangered native salamander, are examined.     

GENETIC ADAPTATION AND ACCLIMATION OF PHYTOPLANKTON ALONG A STRESS GRADIENT IN THE EXTREME WATERS OF THE AGRIO RIVER–CAVIAHUE LAKE (ARGENTINA)1

We tested if different adaptation strategies were linked to a stress gradient in phytoplankton cells. For this purpose, we studied the adaptation and acclimation of Dictyosphaerium chlorelloides (Naumann) Komárek et Perman (Chlorophyta) and Microcystis aeruginosa (Kütz.) Kütz. (Cyanobacteria) to different water samples (from extremely acid, metal‐rich water to moderate stressful conditions) of the Agrio River–Caviahue Lake system (Neuquén, Argentina). Both experimental strains were isolated from pristine, slightly alkaline waters. To distinguish between physiological acclimation and genetic adaptation (an adaptive evolution event), a modified Luria‐Delbrück fluctuation analysis was carried out with both species by using as selective agent sample waters from different points along the stress gradient. M. aeruginosa did not acclimate to any of the waters tested from different points along the stress gradient nor did D. chlorelloides to the two most acidic and metal‐rich waters. However, D. chlorelloides proliferated by rapid genetic adaptation, as the consequence of a single mutation (5.4 × 10^?7 resistant mutants per cell per division) at one locus, in less extreme water and also by acclimation in the least extreme water. It is hypothesized that the stress gradient resulted in different strategies of adaptation in phytoplankton cells from nonextreme waters. Thus, very extreme conditions were lethal for both organisms, but as stressful conditions decreased, adaptation of D. chlorelloides cells was possible by the selection of resistant mutants, and in less extreme conditions, by acclimation.     

Inbreeding depression maintained by recessive lethal mutations interacting with stabilizing selection on quantitative characters in a partially self‐fertilizing population

The bimodal distribution of fitness effects of new mutations and standing genetic variation, due to early‐acting strongly deleterious recessive mutations and late‐acting mildly deleterious mutations, is analyzed using the Kondrashov model for lethals (K), with either the infinitesimal model for selfing (IMS) or the Gaussian allele model (GAM) for quantitative genetic variance under stabilizing selection. In the combined models (KIMS and KGAM) high genomic mutation rates to lethals and weak stabilizing selection on many characters create strong interactions between early and late inbreeding depression, by changing the distribution of lineages selfed consecutively for different numbers of generations. Alternative stable equilibria can exist at intermediate selfing rates for a given set of parameters. Evolution of quantitative genetic variance under multivariate stabilizing selection can strongly influence the purging of nearly recessive lethals, and sometimes vice versa. If the selfing rate at the purging threshold for quantitative genetic variance in IMS or GAM alone exceeds that for nearly recessive lethals in K alone, then in KIMS and KGAM stabilizing selection causes selective interference with purging of lethals, increasing the mean number of lethals compared to K; otherwise, stabilizing selection causes selective facilitation in purging of lethals, decreasing the mean number of lethals.     

Adaptation of the toxic freshwater cyanobacterium Microcystis aeruginosa to salinity is achieved by the selection of spontaneous mutants

The cyanobacterium Microcystis aeruginosa causes most of the harmful toxic blooms in freshwater ecosystems. Some strains of M. aeruginosa tolerate low‐medium levels of salinity, and because salinization of freshwater aquatic systems is increasing worldwide it is relevant to know what adaptive mechanisms allow tolerance to salinity. The mechanisms involved in the adaptation of M. aeruginosa to salinity (acclimation vs. genetic adaptation) were tested by a fluctuation analysis design, and then the maximum capacity of adaptation to salinity was studied by a ratchet protocol experiment. Whereas a dose of 10 g NaCl L^?1 completely inhibited the growth of M. aeruginosa, salinity‐resistant genetic variants, capable of tolerating up to 14 g NaCl L^?1, were isolated in the fluctuation analysis experiment. The salinity‐resistant cells arose by spontaneous mutations at a rate of 7.3 × 10^?7 mutants per cell division. We observed with the ratchet protocol that three independent culture populations of M. aeruginosa were able to adapt to up to 15.1 g L^?1 of NaCl, suggesting that successive mutation‐selection processes can enhance the highest salinity level to which M. aeruginosa cells can initially adapt. We propose that increasing salinity in water reservoirs could lead to the selection of salinity‐resistant mutants of M. aeruginosa.     

Slow Recovery from Inbreeding Depression Generated by the Complex Genetic Architecture of Segregating Deleterious Mutations

The deleterious effects of inbreeding have been of extreme importance to evolutionary biology, but it has been difficult to characterize the complex interactions between genetic constraints and selection that lead to fitness loss and recovery after inbreeding. Haploid organisms and selfing organisms like the nematode Caenorhabditis elegans are capable of rapid recovery from the fixation of novel deleterious mutation; however, the potential for recovery and genomic consequences of inbreeding in diploid, outcrossing organisms are not well understood. We sought to answer two questions: 1) Can a diploid, outcrossing population recover from inbreeding via standing genetic variation and new mutation? and 2) How does allelic diversity change during recovery? We inbred C. remanei, an outcrossing relative of C. elegans, through brother-sister mating for 30 generations followed by recovery at large population size. Inbreeding reduced fitness but, surprisingly, recovery from inbreeding at large populations sizes generated only very moderate fitness recovery after 300 generations. We found that 65% of ancestral single nucleotide polymorphisms (SNPs) were fixed in the inbred population, far fewer than the theoretical expectation of ∼99%. Under recovery, 36 SNPs across 30 genes involved in alimentary, muscular, nervous, and reproductive systems changed reproducibly across replicates, indicating that strong selection for fitness recovery does exist. Our results indicate that recovery from inbreeding depression via standing genetic variation and mutation is likely to be constrained by the large number of segregating deleterious variants present in natural populations, limiting the capacity for recovery of small populations.     

The effects of spontaneous mutation on competitive fitness in Drosophila melanogaster

Abstract We have analysed the effect of 288 generations of mutation accumulation (MA) on chromosome II competitive fitness in 21 full‐sib lines of Drosophila melanogaster and in a large control population, all derived from the same isogenic base. The rate of mean log‐fitness decline and that of increase of the between‐line variance were consistent with a low rate (λ ≈ 0.03 per gamete and generation), and moderate average fitness effect [E(s) ≈ 0.1] of deleterious mutation. Subsequently, crosses were made between pairs of MA lines, and these were maintained with effective size on the order of a few tens. In these crosses, MA recombinant chromosomes quickly recovered to about the average fitness level of control chromosomes. Thus, deleterious mutations responsible for the fitness decline were efficiently selected against in relatively small populations, confirming that their effects were larger than a few percent.     

Sexual selection on spontaneous mutations strengthens the between‐sex genetic correlation for fitness

A proposed benefit to sexual selection is that it promotes purging of deleterious mutations from populations. For this benefit to be realized, sexual selection, which is usually stronger on males, must purge mutations deleterious to both sexes. Here, we experimentally test the hypothesis that sexual selection on males purges deleterious mutations that affect both male and female fitness. We measured male and female fitness in two panels of spontaneous mutation‐accumulation lines of the fly, Drosophila serrata, each established from a common ancestor. One panel of mutation accumulation lines limited both natural and sexual selection (LS lines), whereas the other panel limited natural selection, but allowed sexual selection to operate (SS lines). Although mutation accumulation caused a significant reduction in male and female fitness in both the LS and SS lines, sexual selection had no detectable effect on the extent of the fitness reduction. Similarly, despite evidence of mutational variance for fitness in males and females of both treatments, sexual selection had no significant impact on the amount of mutational genetic variance for fitness. However, sexual selection did reshape the between‐sex correlation for fitness: significantly strengthening it in the SS lines. After 25 generations, the between‐sex correlation for fitness was positive but considerably less than one in the LS lines, suggesting that, although most mutations had sexually concordant fitness effects, sex‐limited, and/or sex‐biased mutations contributed substantially to the mutational variance. In the SS lines this correlation was strong and could not be distinguished from unity. Individual‐based simulations that mimick the experimental setup reveal two conditions that may drive our results: (1) a modest‐to‐large fraction of mutations have sex‐limited (or highly sex‐biased) fitness effects, and (2) the average fitness effect of sex‐limited mutations is larger than the average fitness effect of mutations that affect both sexes similarly.     

Constraints on the evolution of function‐valued traits: a study of growth in Tribolium castaneum

Growth trajectories often impact individual fitness. They are continuous by nature and so are amenable to analysis using a function‐valued (FV) trait framework to reveal their underlying genetic architecture. Previous studies have found high levels of standing additive genetic (co)variance for growth trajectories despite the expectation that growth should be responding to frequent strong directional selection. In this study, the FV framework is used to estimate the additive genetic covariance function for growth trajectories in larval Tribolium castaneum to address questions about standing additive genetic (co)variance and possible evolutionary constraints on growth and to predict responses to four plausible selection regimes. Results show that additive genetic (co)variance is high at the early ages, but decreases towards later ages in the larval period. A selection gradient function of the same size and in the same direction of the first eigenfunction of the G‐function should give the maximal response. However, evolutionary constraints may be acting to keep this maximal response from being realized, through either conflicting effects on survivability and fecundity of larger body size, few evolutionary directions having sufficient additive variance for a response, genetic trade‐offs with other traits or physiological regulatory mechanisms. More light may be shed on these constraints through the development of more sophisticated statistical approaches and implementation of additional empirical studies to explicitly test for specific types of constraints.     

THE CONTRIBUTION OF SPONTANEOUS MUTATIONS TO THERMAL SENSITIVITY CURVE VARIATION IN DROSOPHILA SERRATA

Many traits studied in ecology and evolutionary biology change their expression in response to a continuously varying environmental factor. One well‐studied example are thermal performance curves (TPCs); continuous reaction norms that describe the relationship between organismal performance and temperature and are useful for understanding the trade‐offs involved in thermal adaptation. We characterized curves describing the thermal sensitivity of voluntary locomotor activity in a set of 66 spontaneous mutation accumulation lines in the fly Drosophila serrata. Factor‐analytic modeling of the mutational variance–covariance matrix, M , revealed support for three axes of mutational variation in males and two in females. These independent axes of mutational variance corresponded well to the major axes of TPC variation required for different types of thermal adaptation; “faster‐slower” representing changes in performance largely independent of temperature, and the “hotter‐colder” and “generalist‐specialist” axes, representing trade‐offs. In contrast to its near‐absence from standing variance in this species, a “faster‐slower” axis, accounted for most mutational variance (75% in males and 66% in females) suggesting selection may easily fix or remove these types of mutations in outbred populations. Axes resembling the “hotter‐colder” and “generalist‐specialist” modes of variation contributed less mutational variance but nonetheless point to an appreciable input of new mutations that may contribute to thermal adaptation.