You are where you live: parasitic nematode mitochondrial genome size is associated with the thermal environment generated by hosts

Morphology, Molecular Characterization, and Phylogeny of Travassosius rufus Khalil, 1922 (Strongylidea: Trichostrongylidae), a Parasite from Endangered Sino-Mongolian Beaver (Castor fiber birulai) in Xinjiang, China

The genus Travassosius Khalil, 1922, the smallest genus in the subfamily Trichostrongylinae (family Trichostrongylidae), primarily infects the only two extant beaver species worldwide and can be lethal in severe infections. However, the mitochondrial genome evolution of Travassosius remains poorly understood, and its phylogenetic placement within Trichostrongylinae is still unresolved. In this study, we applied both morphological techniques (differential interference contrast microscopy) and molecular tools (nuclear ITS2 and mitochondrial genome) to examine T. rufus Khalil, 1922. Specimens were collected from the Sino-Mongolian beaver, a subspecies of the Eurasian beaver native to the Ulungur River Basin in northern Xinjiang, China. This work presents the first complete mitochondrial genome sequence and annotation of T. rufus, and it is also the first mitochondrial genome reported for the genus Travassosius. The mitochondrial genome of T. rufus measures 13,646 bp and contains 36 genes, including 12 protein-coding genes (PCGs) (excluding atp8), 22 transfer RNA genes, and 2 ribosomal RNA genes. Phylogenetic analysis based on amino acid sequences of 12 mitochondrial PCGs strongly supports the distinctiveness of the genus Travassosius. Additionally, T. rufus appears to be closely related to Nematodirus within Trichostrongylinae. This study also addresses the possible consequences of parasitic infection for the Sino-Mongolian beaver and offers a scientific foundation for conserving this endangered subspecies and managing parasitic diseases in its population.

Drivers of interlineage variability in mitogenomic evolutionary rates in Platyhelminthes

Studies of forces driving interlineage variability in the evolutionary rates (both sequence and architecture) of mitochondrial genomes often produce contradictory results. Flatworms (Platyhelminthes) exhibit the fastest-evolving mitogenomic sequences among all bilaterian phyla. To test the effects of multiple factors previously associated with different aspects of mitogenomic evolution, we used mitogenomes of 223 flatworm species, phylogenetic multilevel regression models, and causal inference. Thermic host environment (endothermic vs. ectothermic) had nonsignificant impacts on both sequence evolution and mitogenomic size. Mitogenomic gene order rearrangements (GORR) were mostly positively correlated with mitogenomic size (R2 ≈ 20-30%). Longevity was not (negatively) correlated with sequence evolution in flatworms. The predominantly free-living "turbellaria" exhibited much shorter branches and faster-evolving mitogenomic architecture than parasitic Neodermata. As a result, "parasitism" had a strong explanatory power on the branch length variability (>90%), and there was a negative correlation between GORR and branch length. However, the stem branch of Neodermata comprised 63.6% of the total average branch length. This evolutionary period was also marked by a high rate of gene order rearrangements in the ancestral Neodermata. We discuss how this period of rapid evolution deep in the evolutionary history may have decoupled sequence evolution rates from longevity and GORR, and overestimated the explanatory power of "parasitism". This study shows that impacts of variables often vary across lineages, and stresses the importance accounting for the episodic nature of evolutionary patterns in studies of mitogenomic evolution.

The mitochondrial genome of Bottapotamon fukienense (Brachiura: Potamidae) is fragmented into two chromosomes

BACKGROUND

China is the hotspot of global freshwater crab diversity, but their wild populations are facing severe pressures associated with anthropogenic factors, necessitating the need to map their taxonomic and genetic diversity and design conservation policies.

RESULTS

Herein, we sequenced the mitochondrial genome of a Chinese freshwater crab species Bottapotamon fukienense, and found that it is fragmented into two chromosomes. We confirmed that fragmentation was not limited to a single specimen or population. Chromosome 1 comprised 15,111 base pairs (bp) and there were 26 genes and one pseudogene (pseudo-nad1) encoded on it. Chromosome 2 comprised 8,173 bp and there were 12 genes and two pseudogenes (pseudo-trnL2 and pseudo-rrnL) encoded on it. Combined, they comprise the largest mitogenome (23,284 bp) among the Potamidae. Bottapotamon was the only genus in the Potamidae dataset exhibiting rearrangements of protein-coding genes. Bottapotamon fukienense exhibited average rates of sequence evolution in the dataset and did not differ in selection pressures from the remaining Potamidae.

CONCLUSIONS

This is the first experimentally confirmed fragmentation of a mitogenome in crustaceans. While the mitogenome of B. fukienense exhibited multiple signs of elevated mitogenomic architecture evolution rates, including the exceptionally large size, duplicated genes, pseudogenisation, rearrangements of protein-coding genes, and fragmentation, there is no evidence that this is matched by elevated sequence evolutionary rates or changes in selection pressures.

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.

Diagnosis and Characterization of Ditylenchus destructor Isolated from Mazus japonicus in China

The potato rot nematode (Ditylenchus destructor) is one of the most destructive pests in the production of tuber crops, resulting in severely decreased yields and inferior product quality. In 2021, a great number of nematodes were detected in the roots of Mazus japonicus, a weed that is harmful to crop growth, in Qingdao, Shandong Province, China. The present study was undertaken to characterize and identify the nematodes isolated from M. japonicus through morphological identification and molecular approaches. Their morphological characteristics were highly consistent with the descriptions of D. destructor Thorne, 1945. The nematodes collected from M. japonicus were identified as D. destructor haplotype B using D1/D2 and sequence characterized amplified region (SCAR) primers. PCR-ITS-RFLP analysis was conducted to monitor intraspecific variations. In addition, the phylogenetic analysis of the internal transcribed spacer (ITS) demonstrated that this D. destructor population was clustered in haplotype B, supported by a 100% bootstrap value. Another assay, in which M. japonicus was inoculated with a mixture of the life stages of D. destructor, was performed. This assay showed that M. japonicus exhibited a high susceptibility to D. destructor in pots. This is the first record of D. destructor parasitizing M. japonicus in China, and it is of great importance because M. japonicus could be a potential reservoir for D. destructor in the field.

Inverted base composition skews and discontinuous mitochondrial genome architecture evolution in the Enoplea (Nematoda)

Background

Within the class Enoplea, the earliest-branching lineages in the phylum Nematoda, the relatively highly conserved ancestral mitochondrial architecture of Trichinellida is in stark contrast to the rapidly evolving architecture of Dorylaimida and Mermithida. To better understand the evolution of mitogenomic architecture in this lineage, we sequenced the mitogenome of a fish parasite Pseudocapillaria tomentosa (Trichinellida: Capillariidae) and compared it to all available enoplean mitogenomes.

Results

P. tomentosa exhibited highly reduced noncoding regions (the largest was 98 bp), and a unique base composition among the Enoplea. We attributed the latter to the inverted GC skew (0.08) in comparison to the ancestral skew in Trichinellidae (-0.43 to -0.37). Capillariidae, Trichuridae and Longidoridae (Dorylaimida) generally exhibited low negative or low positive skews (-0.1 to 0.1), whereas Mermithidae exhibited fully inverted low skews (0 to 0.05). This is indicative of inversions in the strand replication order or otherwise disrupted replication mechanism in the lineages with reduced/inverted skews. Among the Trichinellida, Trichinellidae and Trichuridae have almost perfectly conserved architecture, whereas Capillariidae exhibit multiple rearrangements of tRNA genes. In contrast, Mermithidae (Mermithida) and Longidoridae (Dorylaimida) exhibit almost no similarity to the ancestral architecture.

Conclusions

Longidoridae exhibited more rearranged mitogenomic architecture than the hypervariable Mermithidae. Similar to the Chromadorea, the evolution of mitochondrial architecture in enoplean nematodes exhibits a strong discontinuity: lineages possessing a mostly conserved architecture over tens of millions of years are interspersed with lineages exhibiting architectural hypervariability. As Longidoridae also have some of the smallest metazoan mitochondrial genomes, they contradict the prediction that compact mitogenomes should be structurally stable. Lineages exhibiting inverted skews appear to represent the intermediate phase between the Trichinellidae (ancestral) and fully derived skews in Chromadorean mitogenomes (GC skews = 0.18 to 0.64). Multiple lines of evidence (CAT-GTR analysis in our study, a majority of previous mitogenomic results, and skew disruption scenarios) support the Dorylaimia split into two sister-clades: Dorylaimida + Mermithida and Trichinellida. However, skew inversions produce strong base composition biases, which can hamper phylogenetic and other evolutionary studies, so enoplean mitogenomes have to be used with utmost care in evolutionary studies.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08607-4.

First Detection of Ditylenchus destructor Parasitizing Maize in Northeast China

Maize is one of the most important crops in the world. Heilongjiang province has the largest maize area in China. Plant-parasitic nematodes are important agricultural pests, which cause huge economic losses every year and have attracted global attention. Potato rot nematode Ditylenchus destructor is a plant-parasitic nematode with a wide range of hosts and strong survival ability in different environments, which brings risks to agricultural production. In 2020, D. destructor was detected in seven maize fields in Heilongjiang province. Morphological identification and molecular approach were used to characterize the isolated D. destructor. The observed morphological and morphometric characteristics were highly similar and consistent with the existing description. The DNA sequencing on the D2/D3 region of the ribosomal DNA 28S and the phylogenetic analysis showed that D. destructor population obtained from maize and other isolates infesting carrot, sweet potato, and potato were in subclade I supported by a 96% bootstrap value. Additionally, the phylogenetic analysis of the ITS rRNA gene sequence further indicated that this D. destructor population from maize clustered in a clade I group and belonged to ITS rRNA haplotype C. An inoculation experiment revealed that D. destructor was pathogenic on the maize seedlings in pots and caused the disease symptoms in the stem base of maize seedlings. This is the first report of D. destructor causing stem rot of maize in Heilongjiang province, and contributes additional information on disease control and safe production of maize in the region.

Slow crabs - fast genomes: Locomotory capacity predicts skew magnitude in crustacean mitogenomes

Base composition skews (G-C/G+C) of mitochondrial genomes are believed to be primarily driven by mutational pressure, which is positively correlated with metabolic rate. In marine animals, metabolic rate is also positively correlated with locomotory capacity. Given the central role of mitochondria in energy metabolism, we hypothesised that selection for locomotory capacity should be positively correlated with the strength of purifying selection (dN/dS), and thus be negatively correlated with the skew magnitude. Therefore, these two models assume diametrically opposite associations between the metabolic rate and skew magnitude: positive correlation in the prevailing paradigm, and negative in our working hypothesis. We examined correlations between the skew magnitude, metabolic rate, locomotory capacity, and several other variables previously associated with mitochondrial evolution on 287 crustacean mitogenomes. Weakly locomotory taxa had higher skew magnitude and ω (dN/dS) values, but not the gene order rearrangement rate. Skew and ω magnitudes were correlated. Multilevel regression analyses indicated that three competing variables, body size, gene order rearrangement rate, and effective population size, had negligible impacts on the skew magnitude. In most crustacean lineages selection for locomotory capacity appears to be the primary factor determining the skew magnitude. Contrary to the prevailing paradigm, this implies that adaptive selection outweighs nonadaptive selection (mutation pressure) in crustaceans. However, we found indications that effective population size (nonadaptive factor) may outweigh the impact of locomotory capacity in sessile crustaceans (Thecostraca). In conclusion, skew magnitude is a product of the interplay between adaptive and nonadaptive factors, the balance of which varies among lineages.

Mitochondrial genomes of the human pathogens Coccidioides immitis and Coccidioides posadasii

Fungal mitochondrial genomes encode genes involved in crucial cellular processes, such as oxidative phosphorylation and mitochondrial translation, and the molecule has been used as a molecular marker for population genetics studies. Coccidioides immitis and C. posadasii are endemic fungal pathogens that cause coccidioidomycosis in arid regions across both American continents. To date, approximately 150 Coccidioides isolates have been sequenced to infer patterns of variation in nuclear genomes. However, less attention has been given to the mitochondrial genomes of Coccidioides. In this report, we describe the assembly and annotation of mitochondrial reference genomes for two representative strains of C. posadasii and C. immitis, as well as assess population variation among 77 selected genomes. The sizes of the circular-mapping molecules are 68.2 Kb in C. immitis and 75.1 Kb in C. posadasii. We identify 14 mitochondrial protein-coding genes common to most fungal mitochondria, which are largely syntenic across different populations and species of Coccidioides. Both Coccidioides species are characterized by a large number of group I and II introns, harboring twice the number of elements as compared to closely related Onygenales. The introns contain complete or truncated ORFs with high similarity to homing endonucleases of the LAGLIDADG and GIY-YIG families. Phylogenetic comparisons of mitochondrial and nuclear genomes show extensive phylogenetic discordance suggesting that the evolution of the two types of genetic material is not identical. This work represents the first assessment of mitochondrial genomes among isolates of both species of Coccidioides, and provides a foundation for future functional work.

Comparative mitogenomics of the Decapoda reveals evolutionary heterogeneity in architecture and composition

The emergence of cost-effective and rapid sequencing approaches has resulted in an exponential rise in the number of mitogenomes on public databases in recent years, providing greater opportunity for undertaking large-scale comparative genomic and systematic research. Nonetheless, current datasets predominately come from small and disconnected studies on a limited number of related species, introducing sampling biases and impeding research of broad taxonomic relevance. This study contributes 21 crustacean mitogenomes from several under-represented decapod infraorders including Polychelida and Stenopodidea, which are used in combination with 225 mitogenomes available on NCBI to investigate decapod mitogenome diversity and phylogeny. An overview of mitochondrial gene orders (MGOs) reveals a high level of genomic variability within the Decapoda, with a large number of MGOs deviating from the ancestral arthropod ground pattern and unevenly distributed among infraorders. Despite the substantial morphological and ecological variation among decapods, there was limited evidence for correlations between gene rearrangement events and species ecology or lineage specific nucleotide substitution rates. Within a phylogenetic context, predicted scenarios of rearrangements show some MGOs to be informative synapomorphies for some taxonomic groups providing strong independent support for phylogenetic relationships. Additional comparisons for a range of mitogenomic features including nucleotide composition, strand asymmetry, unassigned regions and codon usage indicate several clade-specific trends that are of evolutionary and ecological interest.

Mitochondrial Recombination Reveals Mito-Mito Epistasis in Yeast

Genetic variation in mitochondrial DNA (mtDNA) provides adaptive potential although the underlying genetic architecture of fitness components within mtDNAs is not known. To dissect functional variation within mtDNAs, we first identified naturally occurring mtDNAs that conferred high or low fitness in Saccharomyces cerevisiae by comparing growth in strains containing identical nuclear genotypes but different mtDNAs. During respiratory growth under temperature and oxidative stress conditions, mitotype effects were largely independent of nuclear genotypes even in the presence of mito-nuclear interactions. Recombinant mtDNAs were generated to determine fitness components within high- and low-fitness mtDNAs. Based on phenotypic distributions of isogenic strains containing recombinant mtDNAs, we found that multiple loci contributed to mitotype fitness differences. These mitochondrial loci interacted in epistatic, nonadditive ways in certain environmental conditions. Mito-mito epistasis (i.e., nonadditive interactions between mitochondrial loci) influenced fitness in progeny from four different crosses, suggesting that mito-mito epistasis is a widespread phenomenon in yeast and other systems with recombining mtDNAs. Furthermore, we found that interruption of coadapted mito-mito interactions produced recombinant mtDNAs with lower fitness. Our results demonstrate that mito-mito epistasis results in functional variation through mitochondrial recombination in fungi, providing modes for adaptive evolution and the generation of mito-mito incompatibilities.

The complete mitochondrial genome of parasitic nematode Camallanus cotti: extreme discontinuity in the rate of mitogenomic architecture evolution within the Chromadorea class

BACKGROUND

Complete mitochondrial genomes are much better suited for the taxonomic identification and phylogenetic studies of nematodes than morphology or traditionally-used molecular markers, but they remain unavailable for the entire Camallanidae family (Chromadorea). As the only published mitogenome in the Camallanina suborder (Dracunculoidea superfamily) exhibited a unique gene order, the other objective of this research was to study the evolution of mitochondrial architecture in the Spirurida order. Thus, we sequenced the complete mitogenome of the Camallanus cotti fish parasite and conducted structural and phylogenomic comparative analyses with all available Spirurida mitogenomes.

RESULTS

The mitogenome is exceptionally large (17,901 bp) among the Chromadorea and, with 46 (pseudo-) genes, exhibits a unique architecture among nematodes. Six protein-coding genes (PCGs) and six tRNAs are duplicated. An additional (seventh) tRNA (Trp) was probably duplicated by the remolding of tRNA-Ser2 (missing). Two pairs of these duplicated PCGs might be functional; three were incomplete and one contained stop codons. Apart from Ala and Asp, all other duplicated tRNAs are conserved and probably functional. Only 19 unique tRNAs were found. Phylogenomic analysis included Gnathostomatidae (Spirurina) in the Camallanina suborder.

CONCLUSIONS

Within the Nematoda, comparable PCG duplications were observed only in the enoplean Mermithidae family, but those result from mitochondrial recombination, whereas characteristics of the studied mitogenome suggest that likely rearrangement mechanisms are either a series of duplications, transpositions and random loss events, or duplication, fragmentation and subsequent reassembly of the mitogenome. We put forward a hypothesis that the evolution of mitogenomic architecture is extremely discontinuous, and that once a long period of stasis in gene order and content has been punctuated by a rearrangement event, such a destabilised mitogenome is much more likely to undergo subsequent rearrangement events, resulting in an exponentially accelerated evolutionary rate of mitogenomic rearrangements. Implications of this model are particularly important for the application of gene order similarity as an additive source of phylogenetic information. Chromadorean nematodes, and particularly Camallanina clade (with C. cotti as an example of extremely accelerated rate of rearrangements), might be a good model to further study this discontinuity in the dynamics of mitogenomic evolution.

Molecular evolution of a widely-adopted taxonomic marker (COI) across the animal tree of life

DNA barcodes are widely used for identification and discovery of species. While such use draws on information at the DNA level, the current amassment of ca. 4.7 million COI barcodes also offers a unique resource for exploring functional constraints on DNA evolution. Here, we explore amino acid variation in a crosscut of the entire animal kingdom. Patterns of DNA variation were linked to functional constraints at the level of the amino acid sequence in functionally important parts of the enzyme. Six amino acid sites show variation with possible effects on enzyme function. Overall, patterns of amino acid variation suggest convergent or parallel evolution at the protein level connected to the transition into a parasitic life style. Denser sampling of two diverse insect taxa revealed that the beetles (Coleoptera) show more amino acid variation than the butterflies and moths (Lepidoptera), indicating fundamental difference in patterns of molecular evolution in COI. Several amino acid sites were found to be under notably strong purifying selection in Lepidoptera as compared to Coleoptera. Overall, these findings demonstrate the utility of the global DNA barcode library to extend far beyond identification and taxonomy, and will hopefully be followed by a multitude of work.

Mitochondrial Genome Supports Sibling Species of Angiostrongylus costaricensis (Nematoda: Angiostrongylidae)

Angiostrongylus costaricensis is a zoonotic parasitic nematode that causes abdominal or intestinal angiostrongyliasis in humans. It is endemic to the Americas. Although the mitochondrial genome of the Brazil taxon has been published, there is no available mitochondrial genome data on the Costa Rica taxon. We report here the complete mitochondrial genome of the Costa Rica taxon and its genetic differentiation from the Brazil taxon. The whole mitochondrial genome was obtained from next-generation sequencing of genomic DNA. It had a total length of 13,652 bp, comprising 36 genes (12 protein-coding genes—PCGs, 2 rRNA and 22 tRNA genes) and a control region (A + T rich non-coding region). It is longer than that of the Brazil taxon (13,585 bp). The larger mitogenome size of the Costa Rica taxon is due to the size of the control region as the Brazil taxon has a shorter length (265 bp) than the Costa Rica taxon (318 bp). The size of 6 PCGs and the start codon for ATP6, CYTB and NAD5 genes are different between the Costa Rica and Brazil taxa. Additionally, the two taxa differ in the stop codon of 6 PCGs. Molecular phylogeny based on 12 PCGs was concordant with two rRNA, 22 tRNA and 36 mitochondrial genes. The two taxa have a genetic distance of p = 16.2% based on 12 PCGs, p = 15.3% based on 36 mitochondrial genes, p = 13.1% based on 2 rRNA genes and p = 10.7% based on 22 tRNA genes, indicating status of sibling species. The Costa Rica and Brazil taxa of A. costaricensis are proposed to be accorded specific status as members of a species complex.

The complete mitogenome of Apocheima cinerarius (Lepidoptera: Geometridae: Ennominae) and comparison with that of other lepidopteran insects

The complete mitochondrial genome (mitogenome) of a female flightless geometrid moth Apocheima cinerarius was found to be 15,722 bp in length, containing 13 protein-coding genes (PCGs), 22 transfer RNA (tRNA) genes, 2 ribosomal RNA (rRNA) genes, and a control region. The A+T content of the complete mitogenome is 80.83%. The AT skew value ([A-T]/[A+T]) is 0.027. The 13 PCGs of the mitogenome start with typical ATN codons, except for cox1 with the start codon CGA. All the tRNA genes have typical cloverleaf secondary structures, except for trnSer(AGN). The secondary structures of rrnL and rrnS were predicted. Six structural domains including conserved regions (IV, V) and variable regions (I, II, III, VI) were identified in the secondary structure of rrnL. The secondary structure of rrnS consists of 3 structural domains. The control region of A. cinerarius begins with conserved motifs of "ATAGA"+19-bp poly T. It also contains a microsatellite-like (TA)26, a stem-and-loop structure, and a poly-A stretch. Phylogenetic analysis showed that Geometroidea is more closely related to Bombycoidea than to Noctuoidea. A. cinerarius is more closely related to Biston panterinaria than to Phthonandria atrilineata, which is in accordance with the conventional morphology-based classification.

Evolution of parasitism along convergent lines: from ecology to genomics

SUMMARY From hundreds of independent transitions from a free-living existence to a parasitic mode of life, separate parasite lineages have converged over evolutionary time to share traits and exploit their hosts in similar ways. Here, we first summarize the evidence that, at a phenotypic level, eukaryotic parasite lineages have all converged toward only six general parasitic strategies: parasitoid, parasitic castrator, directly transmitted parasite, trophically transmitted parasite, vector-transmitted parasite or micropredator. We argue that these strategies represent adaptive peaks, with the similarities among unrelated taxa within any strategy extending to all basic aspects of host exploitation and transmission among hosts and transcending phylogenetic boundaries. Then, we extend our examination of convergent patterns by looking at the evolution of parasite genomes. Despite the limited taxonomic coverage of sequenced parasite genomes currently available, we find some evidence of parallel evolution among unrelated parasite taxa with respect to genome reduction or compaction, and gene losses or gains. Matching such changes in parasite genomes with the broad phenotypic traits that define the convergence of parasites toward only six strategies of host exploitation is not possible at present. Nevertheless, as more parasite genomes become available, we may be able to detect clear trends in the evolution of parasitic genome architectures representing true convergent adaptive peaks, the genomic equivalents of the phenotypic strategies used by all parasites.