Mutability of mononucleotide repeats, not oxidative stress, explains the discrepancy between laboratory-accumulated mutations and the natural allele-frequency spectrum in C. elegans

Multigenerational effects and mutagenicity of three flame retardants on germ cells in Caenorhabditis elegans

Flame retardants (FRs) have raised public concerns because of their environmental persistence and negative impacts on human health. Recent evidence has revealed that many FRs exhibit reproductive toxicities and transgenerational impacts, whereas the toxic effects of FRs on germ cells remain barely explored. Here we investigated the multigenerational effects of three flame retardants (TBBPA, TCEP and TCPP) on germ cell development in Caenorhabditis elegans, and examined the germ cell mutagenicity of these FRs by using whole genome sequencing. Parental exposure to three FRs markedly increased germ cell apoptosis, and impeded oogenesis in F1-F6 offspring. In addition, the double-increased mutation frequencies observed in progeny genomes uncover the mutagenic actions of FRs on germ cells. Analysis of mutation spectra revealed that these FRs predominantly induced point mutations at A:T base pairs, whereas both small and large indels were almost unaffected. These results revealed the long-term effects of FRs on development and genomic stability of germ cells, which may pose risks to environmental organisms and human reproductive health. Taken together, our findings suggest that germ cell mutagenicity should be carefully examined for the environmental risk assessment of FRs and other emerging pollutants.

Genomic estimates of mutation and substitution rates contradict the evolutionary speed hypothesis of the latitudinal diversity gradient

The latitudinal diversity gradient (LDG) refers to a decrease in biodiversity from the equator to the poles. The evolutionary speed hypothesis, backed by the metabolic theory of ecology, asserts that nucleotide mutation and substitution rates per site per year are higher and thereby speciation rates are higher at higher temperatures, generating the LDG. However, prior empirical investigations of the relationship between the temperature and mutation or substitution rate were based on a few genes and the results were mixed. We here revisit this relationship using genomic data. No significant correlation between the temperature and mutation rate is found in 13 prokaryotes or in 107 eukaryotes. An analysis of 234 diverse trios of bacterial taxa indicates that the synonymous substitution rate is not significantly associated with the growth temperature. The same data, however, reveal a significant negative association between the nonsynonymous substitution rate and temperature, which is explainable by a larger fraction of detrimental nonsynonymous mutations at higher temperatures due to a stronger demand for protein stability. We conclude that the evolutionary speed hypothesis of the LDG is unsupported by genomic data and advise that future mechanistic studies of the LDG should focus on other hypotheses.

Genomic diversity landscapes in outcrossing and selfing Caenorhabditis nematodes

Caenorhabditis nematodes form an excellent model for studying how the mode of reproduction affects genetic diversity, as some species reproduce via outcrossing whereas others can self-fertilize. Currently, chromosome-level patterns of diversity and recombination are only available for self-reproducing Caenorhabditis, making the generality of genomic patterns across the genus unclear given the profound potential influence of reproductive mode. Here we present a whole-genome diversity landscape, coupled with a new genetic map, for the outcrossing nematode C. remanei. We demonstrate that the genomic distribution of recombination in C. remanei, like the model nematode C. elegans, shows high recombination rates on chromosome arms and low rates toward the central regions. Patterns of genetic variation across the genome are also similar between these species, but differ dramatically in scale, being tenfold greater for C. remanei. Historical reconstructions of variation in effective population size over the past million generations echo this difference in polymorphism. Evolutionary simulations demonstrate how selection, recombination, mutation, and selfing shape variation along the genome, and that multiple drivers can produce patterns similar to those observed in natural populations. The results illustrate how genome organization and selection play a crucial role in shaping the genomic pattern of diversity whereas demographic processes scale the level of diversity across the genome as a whole.

Experimental estimates of germline mutation rate in eukaryotes: a phylogenetic meta-analysis

Mutation is the ultimate source of all genetic variation, and over the last 10 years the ready availability of whole-genome sequencing has permitted direct estimation of mutation rate for many non-model species across the tree of life. In this meta-analysis, we make a comprehensive search of the literature for mutation rate estimates in eukaryotes, identifying 140 mutation accumulation (MA) and parent-offspring (PO) sequencing studies covering 134 species. Based on these data, we revisit differences in the single-nucleotide mutation (SNM) rate between different phylogenetic lineages and update the known relationships between mutation rate and generation time, genome size, and nucleotide diversity-while accounting for phylogenetic nonindependence. We do not find a significant difference between MA and PO in estimated mutation rates, but we confirm that mammal and plant lineages have higher mutation rates than arthropods and that unicellular eukaryotes have the lowest mutation rates. We find that mutation rates are higher in species with longer generation times and larger genome sizes, even when accounting for phylogenetic relationships. Moreover, although nucleotide diversity is positively correlated with mutation rate, the gradient of the relationship is significantly less than one (on a logarithmic scale), consistent with higher mutation rates in populations with smaller effective size. For the 29 species for which data are available, we find that indel mutation rates are positively correlated with nucleotide mutation rates and that short deletions are generally more common than short insertions. Nevertheless, despite recent progress, no estimates of either SNM or indel mutation rates are available for the majority of deeply branching eukaryotic lineages-or even for most animal phyla. Even among charismatic megafauna, experimental mutation rate estimates remain unknown for amphibia and scarce for reptiles and fish.

Evaluation of toxicity and mutagenicity of oxaliplatin on germ cells in an alternative in vivo model Caenorhabditis elegans

The platinum compound oxaliplatin is a widely used chemotherapeutic drug that shows a broad spectrum of activity in various human tumors. While the treatment-related side effects of oxaliplatin on directly treated individuals have been well-documented, little is known about the influence of oxaliplatin on germ cells and non-exposed progenies. Here we investigated the reproductive toxicity of oxaliplatin in a 3R-compliant in vivo model Caenorhabditis elegans, and evaluated the germ cell mutagenicity of oxaliplatin by using whole genome sequencing. Our results indicated that oxaliplatin treatment significantly disrupts development of spermatids and oocytes. By treating parental worms with oxaliplatin for three successive generations, sequencing data unveiled the clear mutagenic effects of oxaliplatin on germ cells. Analysis of genome-wide mutation spectra showed the preferentially induction of indels by oxaliplatin. In addition, we uncovered the involvement of translesion synthesis polymerase ζ in modulating mutagenic effects of oxaliplatin. These findings suggest that germ cell mutagenicity is worthy of consideration for the health risk assessment of chemotherapeutic drugs, while the combined use of alternative in vivo models and next generation sequencing technology appears to be a promising way for the preliminary safety assessment of various drugs.

Reproductive toxicity of enrofloxacin in Caenorhabditis elegans involves oxidative stress-induced cell apoptosis

Fluoroquinolone antibiotics (FQs) that persist and bioaccumulate in the environment have aroused people's great concern. Here, we studied the adverse effects of FQs in soil animals of Caenorhabditis elegans via food-chronically exposure. The result shows C. elegans exposed to FQs exhibited reproductive toxicity with small-brood size and low-egg hatchability. To study the underlying mechanism, we conduct a deep investigation of enrofloxacin (ENR), one of the most frequently detected FQs, on nematodes which is one of commonly used animal indicator of soil sustainability. The concentration-effect curves simulated by the Hill model showed that the half effect concentrations (EC₅₀) of ENR were (494.3 ± 272.9) µmol/kg and (107.4 ± 30.9) µmol/kg for the brood size and the hatchability, respectively. Differential gene expression between the control and the ENR-exposure group enriched with the oxidative stress and cell apoptosis pathways. The results together with the enzyme activity in oxidative stress and the cell corpses suggested that ENR-induced reproductive toxicity was related to germ cell apoptosis under oxidative stress. The risk quotients of some soil and livestock samples were calculated based on the threshold value of EC₁₀ for the egg hatchability (2.65 µmol/kg). The results indicated that there was possible reproductive toxicity on the nematodes in certain agricultural soils for the FQs. This study suggested that chronic exposure to FQs at certain levels in environment would induce reproductive toxicity to the nematodes and might reduce the soil sustainability, alarming the environment risks of antibiotics abuse.

Variation in mutational (co)variances

Because of pleiotropy, mutations affect the expression and inheritance of multiple traits and, together with selection, are expected to shape standing genetic covariances between traits and eventual phenotypic divergence between populations. It is therefore important to find if the M matrix, describing mutational variances of each trait and covariances between traits, varies between genotypes. We here estimate the M matrix for six locomotion behavior traits in lines of two genotypes of the nematode Caenorhabditis elegans that accumulated mutations in a nearly neutral manner for 250 generations. We find significant mutational variance along at least one phenotypic dimension of the M matrices, but neither their size nor their orientation had detectable differences between genotypes. The number of generations of mutation accumulation, or the number of MA lines measured, was likely insufficient to sample enough mutations and detect potentially small differences between the two M matrices. We then tested if the M matrices were similar to one G matrix describing the standing genetic (co)variances of a population derived by the hybridization of several genotypes, including the two measured for M, and domesticated to a lab-defined environment for 140 generations. We found that the M and G were different because the genetic covariances caused by mutational pleiotropy in the two genotypes are smaller than those caused by linkage disequilibrium in the lab population. We further show that M matrices differed in their alignment with the lab population G matrix. If generalized to other founder genotypes of the lab population, these observations indicate that selection does not shape the evolution of the M matrix for locomotion behavior in the short-term of a few tens to hundreds of generations and suggests that the hybridization of C. elegans genotypes allows selection on new phenotypic dimensions of locomotion behavior.

Mutation and selection processes regulating short tandem repeats give rise to genetic and phenotypic diversity across species

Short tandem repeats (STRs) are units of 1-6 bp that repeat in a tandem fashion in DNA. Along with single nucleotide polymorphisms and large structural variations, they are among the major genomic variants underlying genetic, and likely phenotypic, divergence. STRs experience mutation rates that are orders of magnitude higher than other well-studied genotypic variants. Frequent copy number changes result in a wide range of alleles, and provide unique opportunities for modulating complex phenotypes through variation in repeat length. While classical studies have identified key roles of individual STR loci, the advent of improved sequencing technology, high-quality genome assemblies for diverse species, and bioinformatics methods for genome-wide STR analysis now enable more systematic study of STR variation across wide evolutionary ranges. In this review, we explore mutation and selection processes that affect STR copy number evolution, and how these processes give rise to varying STR patterns both within and across species. Finally, we review recent examples of functional and adaptive changes linked to STRs.

Whole-Genome Sequencing Reveals Germ Cell Mutagenicity of α-Endosulfan in Caenorhabditis elegans

Endosulfan is an extensively used organochlorine pesticide around the world, which was classified as a persistent organic pollutant (POP) in 2009. Although previous studies have documented the reproductive toxicity of endosulfan in a variety of organisms, little is known about the influence of endosulfan on the genome stability of germ cells and nonexposed progeny. Here we applied whole-genome sequencing to explore the germ cell mutagenicity of α-endosulfan in Caenorhabditis elegans (C. elegans). We found that, although low doses of α-endosulfan exhibited a minor effect on the reproductive capacity of C. elegans, chronic exposure to 1 μM α-endosulfan significantly increased the mutation frequencies of nonexposed progeny. Further analysis of genome-wide mutation spectra demonstrated that α-endosulfan preferentially elicited A:T → G:C substitutions and clustered mutations. By using worms deficient in DNA damage response genes, our results suggest the involvement of translesion synthesis polymerase η in modulating α-endosulfan-induced mutations in germ cells. Together, these observations reveal the germ cell mutagenicity of α-endosulfan in C. elegans and the possible underlying mechanism. In addition, our findings implicate that germ cell mutagenicity might be a necessary consideration for the health risk assessment of environmental chemicals such as POPs.

Salt stress alters the spectrum of de novo mutation available to selection during experimental adaptation of Chlamydomonas reinhardtii

The genetic basis of adaptation is driven by both selection and the spectrum of available mutations. Given that the rate of mutation is not uniformly distributed across the genome and varies depending on the environment, understanding the signatures of selection across the genome is aided by first establishing what the expectations of genetic change are from mutation. To determine the interaction between salt stress, selection, and mutation across the genome, we compared mutations observed in a selection experiment for salt tolerance in Chlamydomonas reinhardtii to those observed in mutation accumulation (MA) experiments with and without salt exposure. MA lines evolved under salt stress had a single-nucleotide mutation rate of 1.1 × 10 - 9 $1.1 \times 10^{-9}$ , similar to that of MA lines under standard conditions ( 9.6 × 10 - 10 $9.6 \times 10^{-10}$ ). However, we found that salt stress led to an increased rate of indel mutations, but that many of these mutations were removed under selection. Finally, lines adapted to salt also showed excess clustering of mutations in the genome and the co-expression network, suggesting a role for positive selection in retaining mutations in particular compartments of the genome during the evolution of salt tolerance. Our study shows that characterizing mutation rates and spectra expected under stress helps disentangle the effects of environment and selection during adaptation.

How Can We Resolve Lewontin's Paradox?

We discuss the genetic, demographic, and selective forces that are likely to be at play in restricting observed levels of DNA sequence variation in natural populations to a much smaller range of values than would be expected from the distribution of census population sizes alone-Lewontin's Paradox. While several processes that have previously been strongly emphasized must be involved, including the effects of direct selection and genetic hitchhiking, it seems unlikely that they are sufficient to explain this observation without contributions from other factors. We highlight a potentially important role for the less-appreciated contribution of population size change; specifically, the likelihood that many species and populations may be quite far from reaching the relatively high equilibrium diversity values that would be expected given their current census sizes.

Mutation, selection, and the prevalence of the Caenorhabditis elegans heat-sensitive mortal germline phenotype

Caenorhabditis elegans strains with the heat-sensitive mortal germline phenotype become progressively sterile over the course of a few tens of generations when maintained at temperatures near the upper range of C. elegans' tolerance. Mortal germline is transgenerationally heritable, and proximately under epigenetic control. Previous studies have suggested that mortal germline presents a relatively large mutational target and that mortal germline is not uncommon in natural populations of C. elegans. The mortal germline phenotype is not monolithic. Some strains exhibit a strong mortal germline phenotype, in which individuals invariably become sterile over a few generations, whereas other strains show a weaker (less penetrant) phenotype in which the onset of sterility is slower and more stochastic. We present results in which we (1) quantify the rate of mutation to the mortal germline phenotype and (2) quantify the frequency of mortal germline in a collection of 95 wild isolates. Over the course of ∼16,000 meioses, we detected one mutation to a strong mortal germline phenotype, resulting in a point estimate of the mutation rate UMrt≈ 6×10-5/genome/generation. We detected no mutations to a weak mortal germline phenotype. Six out of 95 wild isolates have a strong mortal germline phenotype, and although quantification of the weak mortal germline phenotype is inexact, the weak mortal germline phenotype is not rare in nature. We estimate a strength of selection against mutations conferring the strong mortal germline phenotype s¯≈0.1%, similar to selection against mutations affecting competitive fitness. The appreciable frequency of weak mortal germline variants in nature combined with the low mutation rate suggests that mortal germline may be maintained by balancing selection.

How sequence context-dependent mutability drives mutation rate variation in the genome

The rate of mutations varies >100-fold across the genome, altering the rate of evolution, and susceptibility to genetic diseases. The strongest predictor of mutation rate is the sequence itself, varying 75-fold between trinucleotides. The fact that DNA sequence drives its own mutation rate raises a simple but important prediction; highly mutable sequences will mutate more frequently and eliminate themselves in favor of sequences with lower mutability, leading to a lower equilibrium mutation rate. However, purifying selection constrains changes in mutable sequences, causing higher rates of mutation. We conduct a simulation using real human mutation data to test if 1) DNA evolves to a low equilibrium mutation rate and 2) purifying selection causes a higher equilibrium mutation rate in the genome’s most important regions. We explore how this simple process affects sequence evolution in the genome, and discuss the implications for modeling evolution and susceptibility to DNA damage.

Mitochondrial mutations in Caenorhabditis elegans show signatures of oxidative damage and an AT-bias

Rapid mutation rates are typical of mitochondrial genomes (mtDNAs) in animals, but it is not clear why. The difficulty of obtaining measurements of mtDNA mutation that are not biased by natural selection has stymied efforts to distinguish between competing hypotheses about the causes of high mtDNA mutation rates. Several studies which have measured mtDNA mutations in nematodes have yielded small datasets with conflicting conclusions about the relative abundance of different substitution classes (i.e., the mutation spectrum). We therefore leveraged Duplex Sequencing, a high-fidelity DNA sequencing technique, to characterize de novo mtDNA mutations in Caenorhabditis elegans. This approach detected nearly an order of magnitude more mtDNA mutations than documented in any previous nematode mutation study. Despite an existing extreme AT bias in the C. elegans mtDNA (75.6% AT), we found that a significant majority of mutations increase genomic AT content. Compared to some prior studies in nematodes and other animals, the mutation spectrum reported here contains an abundance of CG→AT transversions, supporting the hypothesis that oxidative damage may be a driver of mtDNA mutations in nematodes. Furthermore, we found an excess of G→T and C→T changes on the coding DNA strand relative to the template strand, consistent with increased exposure to oxidative damage. Analysis of the distribution of mutations across the mtDNA revealed significant variation among protein-coding genes and as well as among neighboring nucleotides. This high-resolution view of mitochondrial mutations in C. elegans highlights the value of this system for understanding relationships among oxidative damage, replication error, and mtDNA mutation.