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

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.

The Mitochondrial Genome of a Freshwater Pelagic Amphipod Macrohectopus branickii Is among the Longest in Metazoa

There are more than 350 species of amphipods (Crustacea) in Lake Baikal, which have emerged predominantly through the course of endemic radiation. This group represents a remarkable model for studying various aspects of evolution, one of which is the evolution of mitochondrial (mt) genome architectures. We sequenced and assembled the mt genome of a pelagic Baikalian amphipod species Macrohectopus branickii. The mt genome is revealed to have an extraordinary length (42,256 bp), deviating significantly from the genomes of other amphipod species and the majority of animals. The mt genome of M. branickii has a unique gene order within amphipods, duplications of the four tRNA genes and Cox2, and a long non-coding region, that makes up about two thirds of the genome’s size. The extension of the mt genome was most likely caused by multiple duplications and inversions of regions harboring ribosomal RNA genes. In this study, we analyzed the patterns of mt genome length changes in amphipods and other animal phyla. Through a statistical analysis, we demonstrated that the variability in the mt genome length may be a characteristic of certain phyla and is primarily conferred by expansions of non-coding regions.

Selfing is the safest sex for Caenorhabditis tropicalis

Mating systems have profound effects on genetic diversity and compatibility. The convergent evolution of self-fertilization in three Caenorhabditis species provides a powerful lens to examine causes and consequences of mating system transitions. Among the selfers, Caenorhabditis tropicalis is the least genetically diverse and most afflicted by outbreeding depression. We generated a chromosomal-scale genome for C. tropicalis and surveyed global diversity. Population structure is very strong, and islands of extreme divergence punctuate a genomic background that is highly homogeneous around the globe. Outbreeding depression in the laboratory is caused largely by multiple Medea-like elements, genetically consistent with maternal toxin/zygotic antidote systems. Loci with Medea activity harbor novel and duplicated genes, and their activity is modified by mito-nuclear background. Segregating Medea elements dramatically reduce fitness, and simulations show that selfing limits their spread. Frequent selfing in C. tropicalis may therefore be a strategy to avoid Medea-mediated outbreeding depression.

Mitochondrial dysfunction induces RNA interference in C. elegans through a pathway homologous to the mammalian RIG-I antiviral response

RNA interference (RNAi) is an antiviral pathway common to many eukaryotes that detects and cleaves foreign nucleic acids. In mammals, mitochondrially localized proteins such as mitochondrial antiviral signaling (MAVS), retinoic acid-inducible gene I (RIG-I), and melanoma differentiation-associated protein 5 (MDA5) mediate antiviral responses. Here, we report that mitochondrial dysfunction in Caenorhabditis elegans activates RNAi-directed silencing via induction of a pathway homologous to the mammalian RIG-I helicase viral response pathway. The induction of RNAi also requires the conserved RNA decapping enzyme EOL-1/DXO. The transcriptional induction of eol-1 requires DRH-1 as well as the mitochondrial unfolded protein response (UPR^mt). Upon mitochondrial dysfunction, EOL-1 is concentrated into foci that depend on the transcription of mitochondrial RNAs that may form double-stranded RNA (dsRNA), as has been observed in mammalian antiviral responses. Enhanced RNAi triggered by mitochondrial dysfunction is necessary for the increase in longevity that is induced by mitochondrial dysfunction.

Surveillance of mitochondrial dysfunction in the nematode Caenorhabditis elegans triggers the activation of an RNA interference pathway to mediate antiviral defense, in a manner homologous to the mammalian RIG-I helicase viral response pathway.

Direct full-length RNA sequencing reveals unexpected transcriptome complexity during Caenorhabditis elegans development

Massively parallel sequencing of the polyadenylated RNAs has played a key role in delineating transcriptome complexity, including alternative use of an exon, promoter, 5′ or 3′ splice site or polyadenylation site, and RNA modification. However, reads derived from the current RNA-seq technologies are usually short and deprived of information on modification, compromising their potential in defining transcriptome complexity. Here, we applied a direct RNA sequencing method with ultralong reads using Oxford Nanopore Technologies to study the transcriptome complexity in Caenorhabditis elegans. We generated approximately six million reads using native poly(A)-tailed mRNAs from three developmental stages, with average read lengths ranging from 900 to 1100 nt. Around half of the reads represent full-length transcripts. To utilize the full-length transcripts in defining transcriptome complexity, we devised a method to classify the long reads as the same as existing transcripts or as a novel transcript using sequence mapping tracks rather than existing intron/exon structures, which allowed us to identify roughly 57,000 novel isoforms and recover at least 26,000 out of the 33,500 existing isoforms. The sets of genes with differential expression versus differential isoform usage over development are largely different, implying a fine-tuned regulation at isoform level. We also observed an unexpected increase in putative RNA modification in all bases in the coding region relative to the UTR, suggesting their possible roles in translation. The RNA reads and the method for read classification are expected to deliver new insights into RNA processing and modification and their underlying biology in the future.