Proper mitochondrial function is reliant on favourable mitonuclear epistatic interactions between two genomes, the mitochondrial (mtDNA) and the nuclear with contrasting rates of mutation, copy-number and modes of replication and transmission. The compactness and limited coding capacity of animal mtDNA along with their high mutation rates, uniparental inheritance and lack of recombination has contributed to the view that mtDNA is prone to accumulating deleterious mutations, with mitonuclear adaptation primarily attributed to the nuclear genome. Employing a prospective evolution framework in Caenorhabditis elegans , we investigated the adaptive capacity of the mtDNA in ETC-deficient lines bearing four focal ETC mutants across three different breeding systems following 60 generations of adaptation. The mtDNA genome can adapt rapidly to mitochondrial dysfunction caused by deleterious mutations in nuclear genes. Nonsynonymous changes in mtDNA protein-coding genes had significantly higher frequencies and were positively correlated with fitness, indicative of positive selection. Mitochondrial gene nd-1 emerged as a key driver of compensatory mitonuclear adaptation. Notably, compensatory mutations reach higher frequency when they affected ETC subunits in the same complex as the original deleterious mutation. Our results argue for a disproportionately outsized role for mitochondria in mitonuclear adaptation despite their diminutive genome size and protein-coding capacity.

Spontaneous mutations display biases in their relative frequencies with important consequences for genome structure and composition. While laboratory studies have provided important insights into the spontaneous mutation spectrum, laboratory environments for optimal growth may engender biases that are not representative of natural populations. We analyzed the mitochondrial genomes of 1,524 Caenorhabditis elegans natural isolates comprising 531 unique isotypes to investigate mtDNA polymorphism in the wild. Ancestral reconstruction was used to polarize 2,464 variants (88 indels, 2,376 SNPs) and the results were compared with mutations identified in experimental lines under relaxed selection. MtDNA variant distribution in natural isolates is strongly dependent on site-degeneracy in a manner consistent with purifying selection. There is significant variation in the proportion of synonymous and nonsynonymous polymorphism between genes. Specifically, electron transport chain complex I genes are enriched for nonsynonymous polymorphism. The probability of synonymous mutation is higher at sites with neighboring G/C nucleotides and the per gene synonymous polymorphism is negatively correlated with A+T-content at the first and second codon positions. Furthermore, the 5' and 3' ends of genes have both higher A+T-content and less synonymous polymorphism than central regions. There is evidence of natural selection for preferred codons. We identify the first cases of large heteroplasmic mtDNA structural variants in C. elegans natural isolates, comprising deletions and duplications. Although some patterns of mtDNA mutational bias are similar between laboratory and natural populations, there exist significant differences. In particular, G/C → T/A transversions typically associated with oxidative damage and 8-oxoguanine are strikingly rarer at four-fold degenerate sites in natural populations relative to laboratory populations suggesting that the latter are more prone to mtDNA oxidative damage. Hence, mutational spectra observed in laboratory strains may differ in important aspects from those of natural populations.

Metabolic functioning in nearly all eukaryotes relies on molecular machinery dual-encoded by mitochondrial DNA (mtDNA) and nuclear DNA (nDNA) genomes. The two genomes have sustained an extraordinary degree of cooperation across evolutionary time, preserving the capacity for indispensable processes including oxidative phosphorylation and ATP production, which in turn influence many fitness-related traits. How this cooperation is maintained when one member of the pair is debilitated by deleterious mutation is poorly understood, as is the influence of mutation location (mtDNA or nDNA), mating system, or the potentially compensatory effects of mtDNA copy number changes on the process. We asked whether and to what extent populations experiencing mitonuclear mismatch can recover ancestral levels of fitness by allowing C. elegans nematodes containing either mitochondrial or nuclear mutations of electron transport chain (ETC) genes to evolve under three mating systems–facultatively outcrossing (wildtype), obligately selfing, and obligately outcrossing–for 60 generations. In alignment with evolutionary theory, we observed an inverse relationship between the magnitude of fitness recovery and the ancestral fitness level of strains with the latter outweighing any effect of mating system. We interpret these findings in light of previously reported male frequency evolution in the same mutant lines. The relationship between the amount of fitness evolution and change in mtDNA copy number was influenced by strains’ ETC mutant background and its interaction with mating system. To our knowledge, this work provides the first direct test of the effects of reproductive mode and evolution under mitonuclear mismatch on the population dynamics of mtDNA genomes.

Selfish mitochondrial DNA (mtDNA) mutations are variants that can proliferate within cells and enjoy a replication or transmission bias without fitness benefits for the host. mtDNA deletions in Caenorhabditis elegans can reach high heteroplasmic frequencies despite significantly reducing fitness, illustrating how new mtDNA variants can give rise to genetic conflict between different levels of selection and between the nuclear and mitochondrial genomes. During a mutation accumulation experiment in C. elegans, a 1,034-bp deletion originated spontaneously and reached an 81.7% frequency within an experimental evolution line. This heteroplasmic mtDNA deletion, designated as meuDf1, eliminated portions of 2 protein-coding genes (coxIII and nd4) and tRNA-thr in entirety. mtDNA copy number in meuDf1 heteroplasmic individuals was 35% higher than in individuals with wild-type mitochondria. After backcrossing into a common genetic background, the meuDf1 mitotype was associated with reduction in several fitness traits and independent competition experiments found a 40% reduction in composite fitness. Experiments that relaxed individual selection by single individual bottlenecks demonstrated that the deletion-bearing mtDNA possessed a strong transmission bias, thereby qualifying it as a novel selfish mitotype.

Mitochondrial genomes exist in a nested hierarchy of populations where mitochondrial variants are subject to genetic drift and selection at each level of organization, sometimes engendering conflict between different levels of selection, and between the nuclear and mitochondrial genomes. Deletion mutants in the Caenorhabditis elegans mitochondrial genome can reach high intracellular frequencies despite strongly detrimental effects on fitness. During a mutation accumulation (MA) experiment in C. elegans, a 499 bp deletion in ctb-1 rose to 90% frequency within cells while significantly reducing fitness. During the experiment, the deletion-bearing mtDNA acquired three additional mutations in nd5, namely two single insertion frameshift mutations in a homopolymeric run, and a base substitution. Despite an additional fitness cost of these secondary mutations, all deletion-bearing molecules contained the nd5 mutations at the termination of the MA experiment. The presence of mutant mtDNA was associated with increased mtDNA copy-number. Variation in mtDNA copy-number was greater in the MA lines than in a wildtype nuclear background, including a severe reduction in copy-number at one generational timepoint. Evolutionary replay experiments using different generations of the MA experiment as starting points suggests that two of the secondary mutations contribute to the proliferation of the original ctb-1 deletion by unknown mechanisms.

Aided by new technologies, the upsurgence of research into mitochondrial genome biology during the past 15 years suggests that we have misunderstood, and perhaps dramatically underestimated, the ongoing biological and evolutionary significance of our long-time symbiotic partner. While we have begun to scratch the surface of several topics, many questions regarding the nature of mutation and selection in the mitochondrial genome, and the nature of its relationship to the nuclear genome, remain unanswered. Although best known for their contributions to studies of developmental and aging biology, Caenorhabditis nematodes are increasingly recognized as excellent model systems to advance understanding in these areas. We review recent discoveries with relevance to mitonuclear coevolution and conflict and offer several fertile areas for future work.

We provide a partial test of the mitonuclear sex hypothesis with the first controlled study of how male frequencies and rates of outcrossing evolve in response to mitonuclear mismatch by allowing replicate lineages of C. elegans nematodes containing either mitochondrial or nuclear mutations of electron transport chain (ETC) genes to evolve under three sexual systems: facultatively outcrossing (wildtype), obligately selfing, and obligately outcrossing. Among facultatively outcrossing lines, we found evolution of increased male frequency in at least one replicate line of all four ETC mutant backgrounds tested-nuclear isp-1, mitochondrial cox-1 and ctb-1, and an isp-1 IV; ctb-1M mitonuclear double mutant-and confirmed for a single line set (cox-1) that increased male frequency also resulted in successful outcrossing. We previously found the same result for lines evolved from another nuclear ETC mutant, gas-1. For several lines in the current experiment, however, male frequency declined to wildtype levels (near 0%) in later generations. Male frequency did not change in lines evolved from a wildtype control strain. Additional phenotypic assays of lines evolved from the mitochondrial cox-1 mutant indicated that evolution of high male frequency was accompanied by evolution of increased male sperm size and mating success with tester females, but that it did not translate into increased mating success with coevolved hermaphrodites. Rather, hermaphrodites' self-crossed reproductive fitness increased, consistent with sexually antagonistic coevolution. In accordance with evolutionary theory, males and sexual outcrossing may be most beneficial to populations evolving from a state of low ancestral fitness (gas-1, as previously reported) and less beneficial or deleterious to those evolving from a state of higher ancestral fitness (cox-1). In support of this idea, the obligately outcrossing fog-2 V; cox-1 M lines exhibited no fitness evolution compared to their ancestor, while facultatively outcrossing lines showed slight upward evolution of fitness, and all but one of the obligately selfing xol-1 X; cox-1 M lines evolved substantially increased fitness-even beyond wildtype levels. This work provides a foundation to directly test the effect of reproductive mode on the evolutionary dynamics of mitonuclear genomes, as well as whether compensatory mutations (nuclear or mitochondrial) can rescue populations from mitochondrial dysfunction.

The DNA transposon Tc1 was the first transposable element to be characterized in Caenorhabditis elegans and to date, remains the best-studied transposable element in Caenorhabditis worms. While Tc1 copy-number is regulated at approximately 30 copies in the laboratory Bristol N2 and the vast majority of C. elegans strains, the Bergerac strain and its derivatives have experienced a marked Tc1 proliferation. Given the historical importance of the Bergerac strain in the development of the C. elegans model, we implemented a modern genomic analysis of three Bergerac strains (CB4851, RW6999, and RW7000) in conjunction with multiple phenotypic assays to better elucidate the (1) genomic distribution of Tc1 and (2) phenotypic consequences of transposable element deregulation for the host organism. The median estimates of Tc1 copy-number in the Bergerac strains ranged from 451 to 748, which is both (1) greater than previously estimated and (2) likely to be an underestimate of the actual copy-numbers since coverage-based estimates and digital droplet polymerase chain reaction results both suggest higher Tc1 numbers. All three Bergerac strains had significantly reduced trait means compared with the N2 control for each of four fitness-related traits, with specific traits displaying significant differences between Bergerac strains. Tc1 proliferation was genome-wide, specific to Tc1, and particularly high on chromosomes V and X. There were fewer Tc1 insertions in highly expressed chromatin environments than expected by chance. Furthermore, Tc1 integration motifs were also less frequent in exon than noncoding sequences. The source of the proliferation of Tc1 in the Bergerac strains is specific to Tc1 and independent of other transposable elements. The Bergerac strains contain none of the alleles that have previously been found to derepress transposable element activity in C. elegans. However, the Bergerac strains had several Tc1 insertions near or within highly germline-transcribed genes which could account for the recent germline proliferation.

DNA mismatch repair (MMR), an evolutionarily conserved repair pathway shared by prokaryotic and eukaryotic species alike, influences molecular evolution by detecting and correcting mismatches, thereby protecting genetic fidelity, reducing the mutational load, and preventing lethality. Herein we conduct the first genome-wide evaluation of the alterations to the mutation rate and spectrum under impaired activity of the MutS alpha homolog, msh-2, in Caenorhabditis elegans male-female fog-2(lf) lines. We performed mutation accumulation (MA) under RNAi-induced knockdown of msh-2 for up to 50 generations, followed by next-generation sequencing of 19 MA lines and the ancestral control. msh-2 impairment in the male-female background substantially increased the frequency of nuclear base substitutions (similar to 23x) and small indels (similar to 328x) relative to wildtype hermaphrodites. However, we observed no increase in the mutation rates of mtDNA, and copy-number changes of single-copy genes. There was a marked increase in copy-number variation of rDNA genes under MMR impairment. In C. elegans, msh-2 repairs transitions more efficiently than transversions and increases the AT mutational bias relative to wildtype. The local sequence context, including sequence complexity, G + C-content, and flanking bases influenced the mutation rate. The X chromosome exhibited lower substitution and higher indel rates than autosomes, which can either result from sex-specific mutation rates or a nonrandom distribution of mutable sites between chromosomes. Provided the observed difference in mutational pattern is mostly due to MMR impairment, our results indicate that the specificity of MMR varies between taxa, and is more efficient in detecting and repairing small indels in eukaryotes relative to prokaryotes.

ABSTRACTTo ascertain the effect of relaxation of selection on global gene expression,Caenorhabditis elegansmutation accumulation (MA) lines were propagated under varying degrees of efficiency of selection determined by their different population sizes (N= 1, 10 and 100). Both the mutational variance (Vm), and the residual variance (Vr) were greatest in MA lines with the lowest efficiency of natural selection. The results suggest that gene expression is under strong balancing selection. Furthermore, mutations resulting in increased transcriptional noise or sensitivity to microenvironmental variation accumulate most under extreme genetic drift. In contrast, theVm/Vrratio was lowest in theN=1 lines. Chromatin domains associated with broad gene silencing and active transcription exhibited the greatest and the smallest increase in transcriptional variation, respectively. Furthermore, the preponderance of overexpressed genes was especially pronounced in mitochondrial respiration, stress response, and immune system pathways, especially in low fitnessN= 1 lines.

BackgroundTransposable elements (TEs) are an almost universal constituent of eukaryotic genomes. In animals, Piwi-interacting small RNAs (piRNAs) and repressive chromatin often play crucial roles in preventing TE transcription and thus restricting TE activity. Nevertheless, TE content varies widely across eukaryotes and the dynamics of TE activity and TE silencing across evolutionary time is poorly understood.ResultsHere, we used experimentally evolved populations of C. elegans to study the dynamics of TE expression over 409 generations. The experimental populations were evolved at population sizes of 1, 10 and 100 individuals to manipulate the efficiency of natural selection versus genetic drift. We demonstrate increased TE expression relative to the ancestral population, with the largest increases occurring in the smallest populations. We show that the transcriptional activation of TEs within active regions of the genome is associated with failure of piRNA-mediated silencing, whilst desilenced TEs in repressed chromatin domains retain small RNAs. Additionally, we find that the sequence context of the surrounding region influences the propensity of TEs to lose silencing through failure of small RNA-mediated silencing.ConclusionsOur results show that natural selection in C. elegans is responsible for maintaining low levels of TE expression, and provide new insights into the epigenomic features responsible.

Mitochondrial genomes can sustain mutations that are simultaneously detrimental to individual fitness and yet, can proliferate within individuals owing to a replicative advantage. We analysed the fitness effects and population dynamics of a mitochondrial genome containing a novel 499 bp deletion in the cytochrome b(1) ( ctb-1 ) gene (Δ ctb-1 ) encoding the cytochrome b of complex III in Caenorhabditis elegans. Δ ctb-1 reached a high heteroplasmic frequency of 96% in one experimental line during a mutation accumulation experiment and was linked to additional spontaneous mutations in nd5 and tRNA-Asn . The Δ ctb-1 mutant mitotype imposed a significant fitness cost including a 65% and 52% reduction in productivity and competitive fitness, respectively, relative to individuals bearing wild-type (WT) mitochondria. Deletion-bearing worms were rapidly purged within a few generations when competed against WT mitochondrial DNA (mtDNA) bearing worms in experimental populations. By contrast, the Δ ctb-1 mitotype was able to persist in large populations comprising heteroplasmic individuals only, although the average intracellular frequency of Δ ctb-1 exhibited a slow decline owing to competition among individuals bearing different frequencies of the heteroplasmy. Within experimental lines subjected to severe population bottlenecks ( n = 1), the relative intracellular frequency of Δ ctb-1 increased, which is a hallmark of selfish drive. A positive correlation between Δ ctb-1 and WT mtDNA copy-number suggests a mechanism that increases total mtDNA per se , and does not discern the Δ ctb-1 mitotype from the WT mtDNA. This study demonstrates the selfish nature of the Δ ctb-1 mitotype, given its transmission advantage and substantial fitness load for the host, and highlights the importance of population size for the population dynamics of selfish mtDNA. This article is part of the theme issue ‘Linking the mitochondrial genotype to phenotype: a complex endeavour’.

Experimental investigations into the rates and fitness effects of spontaneous mutations are fundamental to our understanding of the evolutionary process. To gain insights into the molecular and fitness consequences of spontaneous mutations, we conducted a mutation accumulation (MA) experiment at varying population sizes in the nematode Caenorhabditis elegans, evolving 35 lines in parallel for 409 generations at three population sizes (N = 1, 10, and 100 individuals). Here, we focus on nuclear SNPs and small insertion/deletions (indels) under minimal influence of selection, as well as their accrual rates in larger populations under greater selection efficacy. The spontaneous rates of base substitutions and small indels are 1.84 (95% C.I. ± 0.14) × 10-9 substitutions and 6.84 (95% C.I. ± 0.97) × 10-10 changes/site/generation, respectively. Small indels exhibit a deletion bias with deletions exceeding insertions by threefold. Notably, there was no correlation between the frequency of base substitutions, nonsynonymous substitutions, or small indels with population size. These results contrast with our previous analysis of mitochondrial DNA mutations and nuclear copy-number changes in these MA lines, and suggest that nuclear base substitutions and small indels are under less stringent purifying selection compared to the former mutational classes. A transition bias was observed in exons as was a near universal base substitution bias toward A/T. Strongly context-dependent base substitutions, where 5'-Ts and 3'-As increase the frequency of A/T → T/A transversions, especially at the boundaries of A or T homopolymeric runs, manifest as higher mutation rates in (i) introns and intergenic regions relative to exons, (ii) chromosomal cores vs. arms and tips, and (iii) germline-expressed genes.

Mutations spawn genetic variation which, in turn, fuels evolution. Hence, experimental investigations into the rate and fitness effects of spontaneous mutations are central to the study of evolution. Mutation accumulation (MA) experiments have served as a cornerstone for furthering our understanding of spontaneous mutations for four decades. In the pregenomic era, phenotypic measurements of fitness-related traits in MA lines were used to indirectly estimate key mutational parameters, such as the genomic mutation rate, new mutational variance per generation, and the average fitness effect of mutations. Rapidly emerging next-generating sequencing technology has supplanted this phenotype-dependent approach, enabling direct empirical estimates of the mutation rate and a more nuanced understanding of the relative contributions of different classes of mutations to the standing genetic variation. Whole-genome sequencing of MA lines bears immense potential to provide a unified account of the evolutionary process at multiple levels—the genetic basis of variation, and the evolutionary dynamics of mutations under the forces of selection and drift. In this review, we have attempted to synthesize key insights into the spontaneous mutation process that are rapidly emerging from the partnering of classical MA experiments with high-throughput sequencing, with particular emphasis on the spontaneous rates and molecular properties of different mutational classes in nuclear and mitochondrial genomes of diverse taxa, the contribution of mutations to the evolution of gene expression, and the rate and stability of transgenerational epigenetic modifications. Future advances in sequencing technologies will enable greater species representation to further refine our understanding of mutational parameters and their functional consequences.

Significance Copy-number variants are ubiquitous in nature, yet their immediate functional consequences are obscure. We conducted a spontaneous mutation accumulation experiment at varying sizes in Caenorhabditis elegans , thereby enabling the simultaneous investigation of the mutational input and strength of selection on the evolution of copy-number changes. Whole-genome sequencing reveals the highest genome-wide rate of gene duplication for any species thus far. Our transcriptome analysis further demonstrates that gene duplication frequently results in a greater than two-fold change in transcription. Despite the adaptive role of duplication as the primary source of novel genes, we find duplications and deletions of highly transcribed genes to be more detrimental to fitness and evidence for selection against increase in transcript abundance.

Mitochondrial genomes of metazoans, given their elevated rates of evolution, have served as pivotal markers for phylogeographic studies and recent phylogenetic events. In order to determine the dynamics of spontaneous mitochondrial mutations in small populations in the absence and presence of selection, we evolved mutation accumulation (MA) lines of Caenorhabditis elegans in parallel over 409 consecutive generations at three varying population sizes of N = 1, 10, and 100 hermaphrodites. The N =1 populations should have a minimal influence of natural selection to provide the spontaneous mutation rate and the expected rate of neutral evolution, whereas larger population sizes should experience increasing intensity of selection. New mutations were identified by Illumina paired-end sequencing of 86 mtDNA genomes across 35 experimental lines and compared with published genomes of natural isolates. The spontaneous mitochondrial mutation rate was estimated at 1.05 × 10-7/site/generation. A strong G/C→A/T mutational bias was observed in both the MA lines and the natural isolates. This suggests that the low G + C content at synonymous sites is the product of mutation bias rather than selection as previously proposed. The mitochondrial effective population size per worm generation was estimated to be 62. Although it was previously concluded that heteroplasmy was rare in C. elegans, the vast majority of mutations in this study were heteroplasmic despite an experimental regime exceeding 400 generations. The frequencies of frameshift and nonsynonymous mutations were negatively correlated with population size, which suggests their deleterious effects on fitness and a potent role for selection in their eradication.

AbstractHigh‐throughput genomics have revealed widespread copy‐number differences within and among populations of species belonging to diverse taxonomic groups. Experimental evolution in model organisms with minimal natural selection demonstrates that genome‐wide empirical estimates of the spontaneous rates of gene duplication and deletion are extremely high and contribute to the abundance of copy‐number variation (CNV). CNVs are, on average, deleterious with respect to fitness, but their high spontaneous rates of origin can also facilitate rapid adaptation to novel environmental challenges. Duplications and deletions constitute opposing forces that shape genome complexity and size. Gene duplications are the ultimate source of new genes that confer novel phenotypes, whereas deletions remove superfluous genetic material. Furthermore, CNVs can contribute to the evolution of genetic incompatibility and speciation. Future challenges for understanding the evolutionary potential of CNVs include elucidating the relative roles of genetic drift and natural selection for the maintenance of CNVs in populations.Key ConceptsGene duplications are naturally occurring mutations within genomes wherein genic material is duplicated, resulting in additional copies of the duplicated gene.Gene deletions are mutations where genic material has been deleted, thereby reducing the number of copies of the deleted gene.Gene duplications and deletions have resulted in extensive gene copy‐number variation (CNV) in populations across all domains of life.The evolution of gene content and genome size is the net result of duplications adding new genetic material and contributing to the evolution of new genes, and deletions, which are continuously removing genetic information.The study of copy‐number variants at single loci has a long history in population genetics, but a systematic analysis of CNVs across whole genomes became possible only after technical breakthroughs in DNA microarray technology and whole‐genome sequencing.The rates of spontaneous gene duplication and deletion per gene are extraordinarily high and usually much higher than the nucleotide substitution rates.Most CNVs in natural populations are deleterious.Both duplications and deletions can contribute to adaptive genetic variation in natural and experimental populations.The long‐term maintenance of duplicated genes, which contribute to the evolution of novel genes, can be achieved by both positive selection on gene copy‐number and subfunctionalisation resulting from the loss of partial functions of complementary gene copies.Gene duplication and subfunctionalisation can contribute to reproductive isolation and speciation.A comprehensive understanding of how CNVs contribute to the evolution of genomes will require a combination of high‐throughput genomics with analysis of the functional and fitness consequences of CNVs in experimental and natural populations.

Background Gene copy-number variation (CNVs), which provides the raw material for the evolution of novel genes, is widespread in natural populations. We investigated whether CNVs constitute a common mechanism of genetic change during adaptation in experimental Caenorhabditis elegans populations. Outcrossing C. elegans populations with low fitness were evolved for >200 generations. The frequencies of CNVs in these populations were analyzed by oligonucleotide array comparative genome hybridization, quantitative PCR, PCR, DNA sequencing across breakpoints, and single-worm PCR. Results Multiple duplications and deletions rose to intermediate or high frequencies in independent populations. Several lines of evidence suggest that these changes were adaptive: (i) copy-number changes reached high frequency or were fixed in a short time, (ii) many independent populations harbored CNVs spanning the same genes, and (iii) larger average size of CNVs in adapting populations relative to spontaneous CNVs. The latter is expected if larger CNVs are more likely to encompass genes under selection for a change in gene dosage. Several convergent CNVs originated in populations descended from different low fitness ancestors as well as high fitness controls. Conclusions We show that gene copy-number changes are a common class of adaptive genetic change. Due to the high rates of origin of spontaneous duplications and deletions, copy-number changes containing the same genes arose readily in independent populations. Duplications that reached high frequencies in these adapting populations were significantly larger in span. Many convergent CNVs may be general adaptations to laboratory conditions. These results demonstrate the great potential borne by CNVs for evolutionary adaptation.