Selfish Mitochondrial DNA Proliferates and Diversifies in Small, but not Large, Experimental Populations of Caenorhabditis briggsae

Dissecting the sequential evolution of a selfish mitochondrial genome in Caenorhabditis elegans

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.

Comparative mitochondrial genomics in Nematoda reveal astonishing variation in compositional biases and substitution rates indicative of multi-level selection

BACKGROUND

Nematodes are the most abundant and diverse metazoans on Earth, and are known to significantly affect ecosystem functioning. A better understanding of their biology and ecology, including potential adaptations to diverse habitats and lifestyles, is key to understanding their response to global change scenarios. Mitochondrial genomes offer high species level characterization, low cost of sequencing, and an ease of data handling that can provide insights into nematode evolutionary pressures.

RESULTS

Generally, nematode mitochondrial genomes exhibited similar structural characteristics (e.g., gene size and GC content), but displayed remarkable variability around these general patterns. Compositional strand biases showed strong codon position specific G skews and relationships with nematode life traits (especially parasitic feeding habits) equal to or greater than with predicted phylogeny. On average, nematode mitochondrial genomes showed low non-synonymous substitution rates, but also high clade specific deviations from these means. Despite the presence of significant mutational saturation, non-synonymous (dN) and synonymous (dS) substitution rates could still be significantly explained by feeding habit and/or habitat. Low ratios of dN:dS rates, particularly associated with the parasitic lifestyles, suggested the presence of strong purifying selection.

CONCLUSIONS

Nematode mitochondrial genomes demonstrated a capacity to accumulate diversity in composition, structure, and content while still maintaining functional genes. Moreover, they demonstrated a capacity for rapid evolutionary change pointing to a potential interaction between multi-level selection pressures and rapid evolution. In conclusion, this study helps establish a background for our understanding of the potential evolutionary pressures shaping nematode mitochondrial genomes, while outlining likely routes of future inquiry.

Evolutionary codependency: insights into the mitonuclear interaction landscape from experimental and wild Caenorhabditis nematodes

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.

Mitochondrial Inheritance Following Nuclear Transfer: From Cloned Animals to Patients with Mitochondrial Disease

Mitochondria are indispensable power plants of eukaryotic cells that also act as a major biochemical hub. As such, mitochondrial dysfunction, which can originate from mutations in the mitochondrial genome (mtDNA), may impair organism fitness and lead to severe diseases in humans. MtDNA is a multi-copy, highly polymorphic genome that is uniparentally transmitted through the maternal line. Several mechanisms act in the germline to counteract heteroplasmy (i.e., coexistence of two or more mtDNA variants) and prevent expansion of mtDNA mutations. However, reproductive biotechnologies such as cloning by nuclear transfer can disrupt mtDNA inheritance, resulting in new genetic combinations that may be unstable and have physiological consequences. Here, we review the current understanding of mitochondrial inheritance, with emphasis on its pattern in animals and human embryos generated by nuclear transfer.

Stochastic survival of the densest and mitochondrial DNA clonal expansion in aging

The expansion of mitochondrial DNA molecules with deletions has been associated with aging, particularly in skeletal muscle fibers; its mechanism has remained unclear for three decades. Previous accounts have assigned a replicative advantage (RA) to mitochondrial DNA containing deletion mutations, but there is also evidence that cells can selectively remove defective mitochondrial DNA. Here we present a spatial model that, without an RA, but instead through a combination of enhanced density for mutants and noise, produces a wave of expanding mutations with speeds consistent with experimental data. A standard model based on RA yields waves that are too fast. We provide a formula that predicts that wave speed drops with copy number, consonant with experimental data. Crucially, our model yields traveling waves of mutants even if mutants are preferentially eliminated. Additionally, we predict that mutant loads observed in single-cell experiments can be produced by de novo mutation rates that are drastically lower than previously thought for neutral models. Given this exemplar of how spatial structure (multiple linked mtDNA populations), noise, and density affect muscle cell aging, we introduce the mechanism of stochastic survival of the densest (SSD), an alternative to RA, that may underpin other evolutionary phenomena.

Modulating mitochondrial DNA mutations: factors shaping heteroplasmy in the germ line and somatic cells

Until recently it was thought that most humans only harbor one type of mitochondrial DNA (mtDNA), however, deep sequencing and single-cell analysis has shown the converse - that mixed populations of mtDNA (heteroplasmy) are the norm. This is important because heteroplasmy levels can change dramatically during transmission in the female germ line, leading to high levels causing severe mitochondrial diseases. There is also emerging evidence that low level mtDNA mutations contribute to common late onset diseases such as neurodegenerative disorders and cardiometabolic diseases because the inherited mutation levels can change within developing organs and non-dividing cells over time. Initial predictions suggested that the segregation of mtDNA heteroplasmy was largely stochastic, with an equal tendency for levels to increase or decrease. However, transgenic animal work and single-cell analysis have shown this not to be the case during germ-line transmission and in somatic tissues during life. Mutation levels in specific mtDNA regions can increase or decrease in different contexts and the underlying molecular mechanisms are starting to be unraveled. In this review we provide a synthesis of recent literature on the mechanisms of selection for and against mtDNA variants. We identify the most pertinent gaps in our understanding and suggest ways these could be addressed using state of the art techniques.

A comparison of mtDNA deletion mutant proliferation mechanisms

In this paper we use simulation methods to investigate the proliferation of deletion mutations of mitochondrial DNA in neurons. We simulate three mtDNA proliferation mechanisms, namely, random drift, replicative advantage and vicious cycle. For each mechanism, we investigated the effect mutation rates have on neuron loss within a human host. We also compare heteroplasmy of each mechanism at mutation rates that yield the levels neuron loss that would be associated with dementia. Both random drift and vicious cycle predicted high levels of heteroplasmy, while replicative advantage showed a small number of dominant clones with a low background of heteroplasmy.

Mitochondrial Short-Term Plastic Responses and Long-Term Evolutionary Dynamics in Animal Species

How do species respond or adapt to environmental changes? The answer to this depends partly on mitochondrial epigenetics and genetics, new players in promoting adaptation to both short- and long-term environmental changes. In this review, we explore how mitochondrial epigenetics and genetics mechanisms, such as mtDNA methylation, mtDNA-derived noncoding RNAs, micropeptides, mtDNA mutations, and adaptations, can contribute to animal plasticity and adaptation. We also briefly discuss the challenges in assessing mtDNA adaptive evolution. In sum, this review covers new advances in the field of mitochondrial genomics, many of which are still controversial, and discusses processes still somewhat obscure, and some of which are still quite speculative and require further robust experimentation.

Extreme heterogeneity of human mitochondrial DNA from organelles to populations

Contrary to the long-held view that most humans harbour only identical mitochondrial genomes, deep resequencing has uncovered unanticipated extreme genetic variation within mitochondrial DNA (mtDNA). Most, if not all, humans contain multiple mtDNA genotypes (heteroplasmy); specific patterns of variants accumulate in different tissues, including cancers, over time; and some variants are preferentially passed down or suppressed in the maternal germ line. These findings cast light on the origin and spread of mtDNA mutations at multiple scales, from the organelle to the human population, and challenge the conventional view that high percentages of a mutation are required before a new variant has functional consequences.

Disentangling the intertwined roles of mutation, selection and drift in the mitochondrial genome

Understanding and quantifying the rates of change in the mitochondrial genome is a major component of many areas of biological inquiry, from phylogenetics to human health. A critical parameter in understanding rates of change is estimating the mitochondrial mutation rate (mtDNA MR). Although the first direct estimates of mtDNA MRs were reported almost 20 years ago, the number of estimates has not grown markedly since that time. This is largely owing to the challenges associated with time- and labour-intensive mutation accumulation (MA) experiments. But even MA experiments do not solve a major problem with estimating mtDNA MRs-the challenge of disentangling the role of mutation from other evolutionary forces acting within the cell. Now that it is widely understood that any newly generated mutant allele in the mitochondria will initially be at very low frequency (1/N, where N is the number of mtDNA molecules in the cell), the importance of understanding the effective population size (Ne) of the mtDNA and the size of genetic bottlenecks during gametogenesis and development has come into the spotlight. In addition to these factors regulating the role of genetic drift, advances in our understanding of mitochondrial replication and turnover allow us to more easily envision how natural selection within the cell might favour or purge mutations in multi-copy organellar genomes. Here, we review the unique features of the mitochondrial genome that pose a challenge for accurate MR estimation and discuss ways to overcome those challenges. Estimates of mtDNA MRs remain one of the most widely used parameters in biology, thus accurate quantification and a deeper understanding of how and why they may vary within and between individuals, populations and species is an important goal. This article is part of the theme issue 'Linking the mitochondrial genotype to phenotype: a complex endeavour'.

Mitochondrial DNA Variation and Selfish Propagation Following Experimental Bottlenecking in Two Distantly Related Caenorhabditis briggsae Isolates

Understanding mitochondrial DNA (mtDNA) evolution and inheritance has broad implications for animal speciation and human disease models. However, few natural models exist that can simultaneously represent mtDNA transmission bias, mutation, and copy number variation. Certain isolates of the nematode Caenorhabditis briggsae harbor large, naturally-occurring mtDNA deletions of several hundred basepairs affecting the NADH dehydrogenase subunit 5 (nduo-5) gene that can be functionally detrimental. These deletion variants can behave as selfish DNA elements under genetic drift conditions, but whether all of these large deletion variants are transmitted in the same preferential manner remains unclear. In addition, the degree to which transgenerational mtDNA evolution profiles are shared between isolates that differ in their propensity to accumulate the nduo-5 deletion is also unclear. We address these knowledge gaps by experimentally bottlenecking two isolates of C. briggsae with different nduo-5 deletion frequencies for up to 50 generations and performing total DNA sequencing to identify mtDNA variation. We observed multiple mutation profile differences and similarities between C. briggsae isolates, a potentially species-specific pattern of copy number dysregulation, and some evidence for genetic hitchhiking in the deletion-bearing isolate. Our results further support C. briggsae as a practical model for characterizing naturally-occurring mtgenome variation and contribute to the understanding of how mtgenome variation persists in animal populations and how it presents in mitochondrial disease states.

Intracellular quality control of mitochondrial DNA: evidence and limitations

Eukaryotic cells can harbour mitochondria with markedly different transmembrane potentials. Intracellular mitochondrial quality-control mechanisms (e.g. mitophagy) rely on this intracellular variation to distinguish functional and damaged (depolarized) mitochondria. Given that intracellular mitochondrial DNA (mtDNA) genetic variation can induce mitochondrial heterogeneity, mitophagy could remove deleterious mtDNA variants in cells. However, the reliance of mitophagy on the mitochondrial transmembrane potential suggests that mtDNAs with deleterious mutations in ATP synthase can evade the control. This evasion is possible because inhibition of ATP synthase can increase the mitochondrial transmembrane potential. Moreover, the linkage of the mtDNA genotype to individual mitochondrial performance is expected to be weak owing to intracellular mitochondrial intercomplementation. Nonetheless, I reason that intracellular mtDNA quality control is possible and crucial at the zygote stage of the life cycle. Indeed, species with biparental mtDNA inheritance or frequent 'leakage' of paternal mtDNA can be vulnerable to invasion of selfish mtDNAs at the stage of gamete fusion. Here, I critically review recent findings on intracellular mtDNA quality control by mitophagy and discuss other mechanisms by which the nuclear genome can affect the competition of mtDNA variants in the cell. This article is part of the theme issue 'Linking the mitochondrial genotype to phenotype: a complex endeavour'.

The conflict within: origin, proliferation and persistence of a spontaneously arising selfish mitochondrial genome

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'.

A Genome-wide Screen Reveals that Reducing Mitochondrial DNA Polymerase Can Promote Elimination of Deleterious Mitochondrial Mutations

A mutant mitochondrial genome arising amid the pool of mitochondrial genomes within a cell must compete with existing genomes to survive to the next generation. Even weak selective forces can bias transmission of one genome over another to affect the inheritance of mitochondrial diseases and guide the evolution of mitochondrial DNA (mtDNA). Studies in several systems suggested that purifying selection in the female germline reduces transmission of detrimental mitochondrial mutations [1, 2, 3, 4, 5, 6, 7]. In contrast, some selfish genomes can take over despite a cost to host fitness [8, 9, 10, 11, 12, 13]. Within individuals, the outcome of competition is therefore influenced by multiple selective forces. The nuclear genome, which encodes most proteins within mitochondria, and all external regulators of mitochondrial biogenesis and dynamics can influence the competition between mitochondrial genomes [14, 15, 16, 17, 18], yet little is known about how this works. Previously, we established a Drosophila line transmitting two mitochondrial genomes in a stable ratio enforced by purifying selection benefiting one genome and a selfish advantage favoring the other [8]. Here, to find nuclear genes that impact mtDNA competition, we screened heterozygous deletions tiling ∼70% of the euchromatic regions and examined their influence on this ratio. This genome-wide screen detected many nuclear modifiers of this ratio and identified one as the catalytic subunit of mtDNA polymerase gene (POLG), tam. A reduced dose of tam drove elimination of defective mitochondrial genomes. This study suggests that our approach will uncover targets for interventions that would block propagation of pathogenic mitochondrial mutations.

•

Multiple nuclear factors affect selective transmission of mitochondrial genomes

•

Reducing mtDNA polymerase restricts the transmission of detrimental mtDNA mutants

Complex Transmission Patterns and Age-Related Dynamics of a Selfish mtDNA Deletion

Despite wide-ranging implications of selfish mitochondrial DNA (mtDNA) elements for human disease and topics in evolutionary biology (e.g., speciation), the forces controlling their formation, age-related accumulation, and offspring transmission remain largely unknown. Selfish mtDNA poses a significant challenge to genome integrity, mitochondrial function, and organismal fitness. For instance, numerous human diseases are associated with mtDNA mutations; however, few genetic systems can simultaneously represent pathogenic mitochondrial genome evolution and inheritance. The nematode Caenorhabditis briggsae is one such system. Natural C. briggsae isolates harbor varying levels of a large-scale deletion affecting the mitochondrial nduo-5 gene, termed nad5Δ. A subset of these isolates contains putative compensatory mutations that may reduce the risk of deletion formation. We studied the dynamics of nad5Δ heteroplasmy levels during animal development and transmission from mothers to offspring in genetically diverse C. briggsae natural isolates. Results support previous work demonstrating that nad5Δ is a selfish element and that heteroplasmy levels of this deletion can be quite plastic, exhibiting high degrees of inter-family variability and divergence between generations. The latter is consistent with a mitochondrial bottleneck effect, and contrasts with previous findings from a laboratory-derived model uaDf5 mtDNA deletion in C. elegans. However, we also found evidence for among-isolate differences in the ability to limit nad5Δ accumulation, the pattern of which suggested that forces other than the compensatory mutations are important in protecting individuals and populations from rampant mtDNA deletion expansion over short time scales.

Selfish Mitonuclear Conflict

Mitochondria, a nearly ubiquitous feature of eukaryotes, are derived from an ancient symbiosis. Despite billions of years of cooperative coevolution - in what is arguably the most important mutualism in the history of life - the persistence of mitochondrial genomes also creates conditions for genetic conflict with the nucleus. Because mitochondrial genomes are present in numerous copies per cell, they are subject to both within- and among-organism levels of selection. Accordingly, 'selfish' genotypes that increase their own proliferation can rise to high frequencies even if they decrease organismal fitness. It has been argued that uniparental (often maternal) inheritance of cytoplasmic genomes evolved to curtail such selfish replication by minimizing within-individual variation and, hence, within-individual selection. However, uniparental inheritance creates conditions for cytonuclear conflict over sex determination and sex ratio, as well as conditions for sexual antagonism when mitochondrial variants increase transmission by enhancing maternal fitness but have the side-effect of being harmful to males (i.e., 'mother's curse'). Here, we review recent advances in understanding selfish replication and sexual antagonism in the evolution of mitochondrial genomes and the mechanisms that suppress selfish interactions, drawing parallels and contrasts with other organelles (plastids) and bacterial endosymbionts that arose more recently. Although cytonuclear conflict is widespread across eukaryotes, it can be cryptic due to nuclear suppression, highly variable, and lineage-specific, reflecting the diverse biology of eukaryotes and the varying architectures of their cytoplasmic genomes.

Sex and Mitonuclear Adaptation in Experimental Caenorhabditis elegans Populations

To reveal phenotypic and functional genomic patterns of mitonuclear adaptation, a laboratory adaptation study with Caenorhabditis elegans nematodes was conducted in which independently evolving lines were initiated from a low-fitness mitochondrial electron transport chain (ETC) mutant, gas-1 Following 60 generations of evolution in large population sizes with competition for food resources, two distinct classes of lines representing different degrees of adaptive response emerged: a low-fitness class that exhibited minimal or no improvement compared to the gas-1 mutant ancestor, and a high-fitness class containing lines that exhibited partial recovery of wild-type fitness. Many lines that achieved higher reproductive and competitive fitness levels were also noted to evolve high frequencies of males during the experiment, consistent with adaptation in these lines having been facilitated by outcrossing. Whole-genome sequencing and analysis revealed an enrichment of mutations in loci that occur in a gas-1-centric region of the C. elegans interactome and could be classified into a small number of functional genomic categories. A highly nonrandom pattern of mitochondrial DNA mutation was observed within high-fitness gas-1 lines, with parallel fixations of nonsynonymous base substitutions within genes encoding NADH dehydrogenase subunits I and VI. These mitochondrial gene products reside within ETC complex I alongside the nuclear-encoded GAS-1 protein, suggesting that rapid adaptation of select gas-1 recovery lines was driven by fixation of compensatory mitochondrial mutations.

Comparative mitochondrial genomics reveals a possible role of a recent duplication of NADH dehydrogenase subunit 5 in gene regulation

Mitochondrial genome (mtDNA) carries not only well-conserved protein coding, tRNA and rRNA genes, but also highly variable non-coding regions (NCRs). However, the NCRs show poor conservation across species, making their function and evolution elusive. Identification and functional characterization of NCRs across species would be critical for addressing these questions. To this end, we devised a computational pipeline and performed de novo assembly and annotation of mtDNA from 19 Caenorhabditis species using next-generation sequencing (NGS) data. The mtDNAs for 14 out of the 19 species are reported for the first time. Comparison of the 19 genomes reveals species-specific sampling of partial displacement-loop (D-loop) sequence as a novel NCR inserted into a unique tRNA cluster, suggesting an important role of the D-loop and the tRNA cluster in shaping NCR evolution. Intriguingly, RNA-Seq analysis suggests that a novel NCR resulting from a recent duplication of NADH dehydrogenase subunit 5 (ND5) could be utilized as a 3' UTR for up-regulation of its upstream gene. The expression analysis shows a species- and sex-specific expression of mitochondrial genes encoded by mtDNA and nucleus, respectively. Our analyses provide important insights into the function and evolution of mitochondrial NCRs and pave the way for further studying the function and evolution of mitochondrial genome.

Experimental Evolution with Caenorhabditis Nematodes

The hermaphroditic nematode Caenorhabditis elegans has been one of the primary model systems in biology since the 1970s, but only within the last two decades has this nematode also become a useful model for experimental evolution. Here, we outline the goals and major foci of experimental evolution with C. elegans and related species, such as C. briggsae and C. remanei, by discussing the principles of experimental design, and highlighting the strengths and limitations of Caenorhabditis as model systems. We then review three exemplars of Caenorhabditis experimental evolution studies, underlining representative evolution experiments that have addressed the: (1) maintenance of genetic variation; (2) role of natural selection during transitions from outcrossing to selfing, as well as the maintenance of mixed breeding modes during evolution; and (3) evolution of phenotypic plasticity and its role in adaptation to variable environments, including host-pathogen coevolution. We conclude by suggesting some future directions for which experimental evolution with Caenorhabditis would be particularly informative.

The on-again, off-again relationship between mitochondrial genomes and species boundaries

The study of reproductive isolation and species barriers frequently focuses on mitochondrial genomes and has produced two alternative and almost diametrically opposed narratives. On one hand, mtDNA may be at the forefront of speciation events, with co-evolved mitonuclear interactions responsible for some of the earliest genetic incompatibilities arising among isolated populations. On the other hand, there are numerous cases of introgression of mtDNA across species boundaries even when nuclear gene flow is restricted. We argue that these seemingly contradictory patterns can result from a single underlying cause. Specifically, the accumulation of deleterious mutations in mtDNA creates a problem with two alternative evolutionary solutions. In some cases, compensatory or epistatic changes in the nuclear genome may ameliorate the effects of mitochondrial mutations, thereby establishing coadapted mitonuclear genotypes within populations and forming the basis of reproductive incompatibilities between populations. Alternatively, populations with high mitochondrial mutation loads may be rescued by replacement with a more fit, foreign mitochondrial haplotype. Coupled with many nonadaptive mechanisms of introgression that can preferentially affect cytoplasmic genomes, this form of adaptive introgression may contribute to the widespread discordance between mitochondrial and nuclear genealogies. Here, we review recent advances related to mitochondrial introgression and mitonuclear incompatibilities, including the potential for cointrogression of mtDNA and interacting nuclear genes. We also address an emerging controversy over the classic assumption that selection on mitochondrial genomes is inefficient and discuss the mechanisms that lead lineages down alternative evolutionary paths in response to mitochondrial mutation accumulation.

Mitochondrial depolarization in yeast zygotes inhibits clonal expansion of selfish mtDNA

Non-identical copies of mitochondrial DNA (mtDNA) compete with each other within a cell and the ultimate variant of mtDNA present depends on their relative replication rates. Using yeast Saccharomyces cerevisiae cells as a model, we studied the effects of mitochondrial inhibitors on the competition between wild-type mtDNA and mutant selfish mtDNA in heteroplasmic zygotes. We found that decreasing mitochondrial transmembrane potential by adding uncouplers or valinomycin changes the competition outcomes in favor of the wild-type mtDNA. This effect was significantly lower in cells with disrupted mitochondria fission or repression of the autophagy-related genes ATG8, ATG32 or ATG33, implying that heteroplasmic zygotes activate mitochondrial degradation in response to the depolarization. Moreover, the rate of mitochondrially targeted GFP turnover was higher in zygotes treated with uncoupler than in haploid cells or untreated zygotes. Finally, we showed that vacuoles of zygotes with uncoupler-activated autophagy contained DNA. Taken together, our data demonstrate that mitochondrial depolarization inhibits clonal expansion of selfish mtDNA and this effect depends on mitochondrial fission and autophagy. These observations suggest an activation of mitochondria quality control mechanisms in heteroplasmic yeast zygotes.

Mitochondrial Genomes of Kinorhyncha: trnM Duplication and New Gene Orders within Animals

Many features of mitochondrial genomes of animals, such as patterns of gene arrangement, nucleotide content and substitution rate variation are extensively used in evolutionary and phylogenetic studies. Nearly 6,000 mitochondrial genomes of animals have already been sequenced, covering the majority of animal phyla. One of the groups that escaped mitogenome sequencing is phylum Kinorhyncha-an isolated taxon of microscopic worm-like ecdysozoans. The kinorhynchs are thought to be one of the early-branching lineages of Ecdysozoa, and their mitochondrial genomes may be important for resolving evolutionary relations between major animal taxa. Here we present the results of sequencing and analysis of mitochondrial genomes from two members of Kinorhyncha, Echinoderes svetlanae (Cyclorhagida) and Pycnophyes kielensis (Allomalorhagida). Their mitochondrial genomes are circular molecules approximately 15 Kbp in size. The kinorhynch mitochondrial gene sequences are highly divergent, which precludes accurate phylogenetic inference. The mitogenomes of both species encode a typical metazoan complement of 37 genes, which are all positioned on the major strand, but the gene order is distinct and unique among Ecdysozoa or animals as a whole. We predict four types of start codons for protein-coding genes in E. svetlanae and five in P. kielensis with a consensus DTD in single letter code. The mitochondrial genomes of E. svetlanae and P. kielensis encode duplicated methionine tRNA genes that display compensatory nucleotide substitutions. Two distant species of Kinorhyncha demonstrate similar patterns of gene arrangements in their mitogenomes. Both genomes have duplicated methionine tRNA genes; the duplication predates the divergence of two species. The kinorhynchs share a few features pertaining to gene order that align them with Priapulida. Gene order analysis reveals that gene arrangement specific of Priapulida may be ancestral for Scalidophora, Ecdysozoa, and even Protostomia.

Positive Selection in Rapidly Evolving Plastid-Nuclear Enzyme Complexes

Rates of sequence evolution in plastid genomes are generally low, but numerous angiosperm lineages exhibit accelerated evolutionary rates in similar subsets of plastid genes. These genes include clpP1 and accD, which encode components of the caseinolytic protease (CLP) and acetyl-coA carboxylase (ACCase) complexes, respectively. Whether these extreme and repeated accelerations in rates of plastid genome evolution result from adaptive change in proteins (i.e., positive selection) or simply a loss of functional constraint (i.e., relaxed purifying selection) is a source of ongoing controversy. To address this, we have taken advantage of the multiple independent accelerations that have occurred within the genus Silene (Caryophyllaceae) by examining phylogenetic and population genetic variation in the nuclear genes that encode subunits of the CLP and ACCase complexes. We found that, in species with accelerated plastid genome evolution, the nuclear-encoded subunits in the CLP and ACCase complexes are also evolving rapidly, especially those involved in direct physical interactions with plastid-encoded proteins. A massive excess of nonsynonymous substitutions between species relative to levels of intraspecific polymorphism indicated a history of strong positive selection (particularly in CLP genes). Interestingly, however, some species are likely undergoing loss of the native (heteromeric) plastid ACCase and putative functional replacement by a duplicated cytosolic (homomeric) ACCase. Overall, the patterns of molecular evolution in these plastid-nuclear complexes are unusual for anciently conserved enzymes. They instead resemble cases of antagonistic coevolution between pathogens and host immune genes. We discuss a possible role of plastid-nuclear conflict as a novel cause of accelerated evolution.

The mitochondrial genome of Globodera ellingtonae is composed of two circles with segregated gene content and differential copy numbers

Background

The evolution of animal mitochondrial (mt) genomes has resulted in a highly conserved structure: a single compact circular chromosome approximately 14 to 20 kb long. Within the last two decades exceptions to this conserved structure, such as the division of the genome into multiple chromosomes, have been reported in a diverse set of metazoans. We report on the two circle multipartite mt genome of a newly described cyst nematode, Globodera ellingtonae.

Results

The G. ellingtonae mt genome was found to be comprised of two circles, each larger than any other multipartite circular mt chromosome yet reported, and both were larger than the single mt circle of the model nematode Caenorhabditis elegans. The genetic content of the genome was disproportionately divided between the two circles, although they shared a ~6.5 kb non-coding region. The 17.8 kb circle (mtDNA-I) contained ten protein-coding genes and two tRNA genes, whereas the 14.4 kb circle (mtDNA-II) contained two protein-coding genes, 20 tRNA genes and both rRNA genes. Perhaps correlated with this division of genetic content, the copy number of mtDNA-II was more than four-fold that of mtDNA-I in individual nematodes. The difference in copy number increased between second-stage and fourth-stage juveniles.

Conclusions

The segregation of gene types to different mt circles in G. ellingtonae could provide benefit by localizing gene functional types to independent transcriptional units. This is the first report of both two-circle and several-circle mt genomes within a single genus. The differential copy number associated with this multipartite mt organization could provide a model system for deconstructing mechanisms regulating mtDNA copy number both in somatic cells and during germline development.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-3047-x) contains supplementary material, which is available to authorized users.

Homeostatic Responses Regulate Selfish Mitochondrial Genome Dynamics in C. elegans

Mutant mitochondrial genomes (mtDNA) can be viewed as selfish genetic elements that persist in a state of heteroplasmy despite having potentially deleterious metabolic consequences. We sought to study regulation of selfish mtDNA dynamics. We establish that the large 3.1-kb deletion-bearing mtDNA variant uaDf5 is a selfish genome in Caenorhabditis elegans. Next, we show that uaDf5 mutant mtDNA replicates in addition to, not at the expense of, wild-type mtDNA. These data suggest the existence of a homeostatic copy-number control that is exploited by uaDf5 to "hitchhike" to high frequency. We also observe activation of the mitochondrial unfolded protein response (UPR(mt)) in uaDf5 animals. Loss of UPR(mt) causes a decrease in uaDf5 frequency, whereas its constitutive activation increases uaDf5 levels. UPR(mt) activation protects uaDf5 from mitophagy. Taken together, we propose that mtDNA copy-number control and UPR(mt) represent two homeostatic response mechanisms that play important roles in regulating selfish mitochondrial genome dynamics.

Selfish drive can trump function when animal mitochondrial genomes compete

Mitochondrial genomes compete for transmission from mother to progeny. We explored this competition by introducing a second genome into Drosophila melanogaster to follow transmission. Competitions between closely related genomes favored those functional in electron transport, resulting in a host-beneficial purifying selection. In contrast, matchups between distantly related genomes often favored those with negligible, negative or lethal consequences, indicating selfish selection. Exhibiting powerful selfish selection, a genome carrying a detrimental mutation displaced a complementing genome, leading to population death after several generations. In a different pairing, opposing selfish and purifying selection counterbalanced to give stable transmission of two genomes. Sequencing of recombinant mitochondrial genomes showed that the noncoding region, containing origins of replication, governs selfish transmission. Uniparental inheritance prevents encounters between distantly related genomes. Nonetheless, in each maternal lineage, constant competition among sibling genomes selects for super-replicators. We suggest that this relentless competition drives positive selection, promoting change in the sequences influencing transmission.

Paths of Heritable Mitochondrial DNA Mutation and Heteroplasmy in Reference and gas-1 Strains of Caenorhabditis elegans

Heteroplasmy—the presence of more than one mitochondrial DNA (mtDNA) sequence type in a cell, tissue, or individual—impacts human mitochondrial disease and numerous aging-related syndromes. Understanding the trans-generational dynamics of mtDNA is critical to understanding the underlying mechanisms of mitochondrial disease and evolution. We investigated mtDNA mutation and heteroplasmy using a set of wild-type (N2 strain) and mitochondrial electron transport chain (ETC) mutant (gas-1) mutant Caenorhabditis elegans mutation-accumulation (MA) lines. The N2 MA lines, derived from a previous experiment, were bottlenecked for 250 generations. The gas-1 MA lines were created for this study, and bottlenecked in the laboratory for up to 50 generations. We applied Illumina-MiSeq DNA sequencing to L1 larvae from five gas-1 MA lines and five N2 MA lines to detect and characterize mtDNA mutation and heteroplasmic inheritance patterns evolving under extreme drift. mtDNA copy number increased in both sets of MA lines: three-fold on average among the gas-1 MA lines and five-fold on average among N2 MA lines. Eight heteroplasmic single base substitution polymorphisms were detected in the gas-1 MA lines; only one was observed in the N2 MA lines. Heteroplasmy frequencies ranged broadly in the gas-1 MA lines, from as low as 2.3% to complete fixation (homoplasmy). An initially low-frequency (<5%) heteroplasmy discovered in the gas-1 progenitor was observed to fix in one gas-1 MA line, achieve higher frequency (37.4%) in another, and be lost in the other three lines. A similar low-frequency heteroplasmy was detected in the N2 progenitor, but was lost in all five N2 MA lines. We identified three insertion-deletion (indel) heteroplasmies in gas-1 MA lines and six indel variants in the N2 MA lines, most occurring at homopolymeric nucleotide runs. The observed bias toward accumulation of single nucleotide polymorphisms in gas-1 MA lines is consistent with the idea that impaired mitochondrial activity renders mtDNA more vulnerable to this type of mutation. The consistent increases in mtDNA copy number implies that extreme genetic drift provides a permissive environment for elevated organelle genome copy number in C. elegans reference and gas-1 strains. This study broadens our understanding of the heteroplasmic mitochondrial mutation process in a multicellular model organism.

Mitogenomics reveals high synteny and long evolutionary histories of sympatric cryptic nematode species

Species with seemingly identical morphology but with distinct genetic differences are abundant in the marine environment and frequently co-occur in the same habitat. Such cryptic species are typically delineated using a limited number of mitochondrial and/or nuclear marker genes, which do not yield information on gene order and gene content of the genomes under consideration. We used next-generation sequencing to study the composition of the mitochondrial genomes of four sympatrically distributed cryptic species of the Litoditis marina species complex (PmI, PmII, PmIII, and PmIV). The ecology, biology, and natural occurrence of these four species are well known, but the evolutionary processes behind this cryptic speciation remain largely unknown. The gene order of the mitochondrial genomes of the four species was conserved, but differences in genome length, gene length, and codon usage were observed. The atp8 gene was lacking in all four species. Phylogenetic analyses confirm that PmI and PmIV are sister species and that PmIII diverged earliest. The most recent common ancestor of the four cryptic species was estimated to have diverged 16 MYA. Synonymous mutations outnumbered nonsynonymous changes in all protein-encoding genes, with the Complex IV genes (coxI-III) experiencing the strongest purifying selection. Our mitogenomic results show that morphologically similar species can have long evolutionary histories and that PmIII has several differences in genetic makeup compared to the three other species, which may explain why it is better adapted to higher temperatures than the other species.