Mitochondrial Architecture Rearrangements Produce Asymmetrical Nonadaptive Mutational Pressures That Subvert the Phylogenetic Reconstruction in Isopoda

Complete mitochondrial genomes of Dactylogyrus crucifer and Dactylogyrus zandti reveal distinct patterns of codon usage within Dactylogyrus

Monogeneans of the genus Dactylogyrus Diesing, 1850, the largest genus in the family Dactylogyridae, mostly parasitize the gills of cyprinoid hosts; however, only 3 Dactylogyrus' mitochondrial genomes (mitogenomes) are studied so far. The aim of this research is to extend our understanding of the mitogenomes of Dactylogyrus. We sequenced the mitogenomes of D. crucifer and D. zandti isolated from Rutilus rutilus and Abramis brama orientalis in northwest China, and then we compared these mitogenomes with other monogeneans. We used Illumina NovaSeq to sequence the entire mitochondrial genomes of D. crucifer and D. zandti and characterized the mitogenomes to understand the gene structure, gene identity, the secondary structures of the 22 tRNA genes, and relative synonymous codon usage. We used the analytic Bayesian Information and Maximum Likelihood methods to determine their associated phylogenetic trees. The mitogenomes of D. crucifer and D. zandti were 14,403 and 18,584 bp, respectively. Organization and positioning of these genes were in accordance with Dactylogyrus lamellatus and Dactylogyrus tuba. The nucleotide composition of Dactylogyridae was different from other families of Monogenea, and the A+T count of genus Dactylogyrus (54 - 58.4 %) was lower than other genus species of the family Dactylogyridea (63.9 - 78.4 %) in protein-coding genes. Dactylogyrus members displayed a codon usage bias. The relative synonymous codon used by Dactylogyrus was not conserved and was lower than other monogeneans. The codon use patterns of closely-related species isolated from closely-related hosts were identical. Phylogenetic analyses using mitogenomic dataset produced Dactylogyrus isolated from host subfamily Leuciscinae formed a sister-group. Our results contributed significantly to an increased database of mitogenomes, more than 50 %, for Dactylogyrus that may help future studies of mitochondrial genes and codon uses for the analysis of monogenean phylogenetics.

Complete mitochondrial genome and phylogenetic analysis of Dollfustrema vaneyi (Trematoda: Bucephalidae)

BACKGROUND

The Bucephalidae is a large family of digenean trematodes but most previous analyses of its phylogenetic position have relied on a single mitochondrial gene or morphological features. Mitochondrial genomes (mitogenomes) remain unavailable for the entire family. To address this, we sequenced the complete mitogenome of Dollfustrema vaneyi and analyzed the phylogenetic relationships with other trematodes.

RESULTS

The circular genome of Dollfustrema vaneyi spanned 14,959 bp and contained 12 protein-coding genes, 22 transfer RNA genes, 2 ribosomal RNA genes, and a major non-coding region. We used concatenated amino acid and nucleotide sequences of all 36 genes for phylogenetic analyses, conducted using MrBayes, IQ-TREE and PhyloBayes. We identified pronounced topological instability across different analyses. The addition of recently sequenced two mitogenomes for the Aspidogastrea subclass along with the use of a site-heterogeneous model stabilized the topology, particularly the positions of Azygiidae and Bucephalidae. The stabilized results indicated that Azygiidae was the closest lineage to Bucephalidae in the available dataset, and together, they clustered at the base of the Plagiorchiida.

CONCLUSIONS

Our study provides the first comprehensive description and annotation of the mitochondrial genome for the Bucephalidae family. The results indicate a close phylogenetic relationship between Azygiidae and Bucephalidae, and reveal their basal placement within the order Plagiorchiida. Furthermore, the inclusion of Aspidogastrea mitogenomes and the site-heterogeneous model significantly improved the topological stability. These data will provide key molecular resources for future taxonomic and phylogenetic studies of the family Bucephalidae.

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 mitogenomic analysis of subterranean and surface amphipods (Crustacea, Amphipoda) with special reference to the family Crangonyctidae

Mitochondrial genomes play important roles in studying genome evolution, phylogenetic analyses, and species identification. Amphipods (Class Malacostraca, Order Amphipoda) are one of the most ecologically diverse crustacean groups occurring in a diverse array of aquatic and terrestrial environments globally, from freshwater streams and lakes to groundwater aquifers and the deep sea, but we have a limited understanding of how habitat influences the molecular evolution of mitochondrial energy metabolism. Subterranean amphipods likely experience different evolutionary pressures on energy management compared to surface-dwelling taxa that generally encounter higher levels of predation and energy resources and live in more variable environments. In this study, we compared the mitogenomes, including the 13 protein-coding genes involved in the oxidative phosphorylation (OXPHOS) pathway, of surface and subterranean amphipods to uncover potentially different molecular signals of energy metabolism between surface and subterranean environments in this diverse crustacean group. We compared base composition, codon usage, gene order rearrangement, conducted comparative mitogenomic and phylogenomic analyses, and examined evolutionary signals of 35 amphipod mitogenomes representing 13 families, with an emphasis on Crangonyctidae. Mitogenome size, AT content, GC-skew, gene order, uncommon start codons, location of putative control region (CR), length of rrnL and intergenic spacers differed between surface and subterranean amphipods. Among crangonyctid amphipods, the spring-dwelling Crangonyx forbesi exhibited a unique gene order, a long nad5 locus, longer rrnL and rrnS loci, and unconventional start codons. Evidence of directional selection was detected in several protein-encoding genes of the OXPHOS pathway in the mitogenomes of surface amphipods, while a signal of purifying selection was more prominent in subterranean species, which is consistent with the hypothesis that the mitogenome of surface-adapted species has evolved in response to a more energy demanding environment compared to subterranean amphipods. Overall, gene order, locations of non-coding regions, and base-substitution rates points to habitat as an important factor influencing the evolution of amphipod mitogenomes.

Mitogenomic evolutionary rates in bilateria are influenced by parasitic lifestyle and locomotory capacity

The evidence that parasitic animals exhibit elevated mitogenomic evolutionary rates is inconsistent and limited to Arthropoda. Similarly, the evidence that mitogenomic evolution is faster in species with low locomotory capacity is limited to a handful of animal lineages. We hypothesised that these two variables are associated and that locomotory capacity is a major underlying factor driving the elevated rates in parasites. Here, we study the evolutionary rates of mitogenomes of 10,906 bilaterian species classified according to their locomotory capacity and parasitic/free-living life history. In Bilateria, evolutionary rates were by far the highest in endoparasites, much lower in ectoparasites with reduced locomotory capacity and free-living lineages with low locomotory capacity, followed by parasitoids, ectoparasites with high locomotory capacity, and finally micropredatory and free-living lineages. The life history categorisation (parasitism) explained ≈45%, locomotory capacity categorisation explained ≈39%, and together they explained ≈56% of the total variability in evolutionary rates of mitochondrial protein-coding genes in Bilateria. Our findings suggest that these two variables play major roles in calibrating the mitogenomic molecular clock in bilaterian animals.

The complete mitochondrial genomes of the Leuciscus baicalensis and Rutilus rutilus: a detailed genomic comparison among closely related species of the Leuciscinae subfamily

Cyprinidae is the largest family in the order of freshwater fish Cypriniformes. Increased subfamily members of Cyprinidae have been suggested to be re-classified for decades. In this study, we sequenced the mitochondrial genomes (mitogenomes) of Leuciscus baicalensis and Rutilus rutilus collected from northwest China and compared with other closely related species to determine their associated family or subfamily. We used Illumina NovaSeq to sequence the entire mitochondrial genomes of Leuciscus baicalensis and Rutilus rutilus and characterized the mitogenomes by the gene structure, gene order, and the secondary structures of the 22 tRNA genes. We compared mitogenome features of Leuciscinae with other subfamilies in Cyprinidae. We used the analytic Bayesian Information and Maximum Likelihood methods to determine phylogenetic trees of 13 PCGs. The mitogenomes of Leuciscus baicalensis and Rutilus rutilus were 16,607 bp and 16,606 bp, respectively. Organization and location of these genes were consistent with already studied Leuciscinae fishes. Synonymous codon usage was conservative in Leuciscinae as compared with other subfamilies in Cyprinidae. Phylogenetic analysis indicated that Leuciscinae was a monophyletic group, and genus Leuciscus was a paraphyletic group. Our approach, for the first time, of studying comparative mitochondrial genomics and phylogenetics together provided a supportive platform to the analysis of population genetics and phylogeny for Leuciscinae. Our results indicated a promising potential of comparative mitochondrial genomics in the manifestation of phylogenetic relationships between fishes, leading us to a suggestion that mitogenomes should be routinely considered in clarifying phylogenetics of family and subfamily members of fish.

Comparative Mitogenome Analyses Uncover Mitogenome Features and Phylogenetic Implications of the Parrotfishes (Perciformes: Scaridae)

In order to investigate the molecular evolution of mitogenomes among the family Scaridae, the complete mitogenome sequences of twelve parrotfish species were determined and compared with those of seven other parrotfish species. The comparative analysis revealed that the general features and organization of the mitogenome were similar among the 19 parrotfish species. The base composition was similar among the parrotfishes, with the exception of the genus Calotomus, which exhibited an unusual negative AT skew in the whole mitogenome. The PCGs showed similar codon usage, and all of them underwent a strong purifying selection. The gene rearrangement typical of the parrotfishes was detected, with the tRNAMet inserted between the tRNAIle and tRNAGln, and the tRNAGln was followed by a putative tRNAMet pseudogene. The parrotfish mitogenomes displayed conserved gene overlaps and secondary structure in most tRNA genes, while the non-coding intergenic spacers varied among species. Phylogenetic analysis based on the thirteen PCGs and two rRNAs strongly supported the hypothesis that the parrotfishes could be subdivided into two clades with distinct ecological adaptations. The early divergence of the sea grass and coral reef clades occurred in the late Oligocene, probably related to the expansion of sea grass habitat. Later diversification within the coral reef clade could be dated back to the Miocene, likely associated with the geomorphology alternation since the closing of the Tethys Ocean. This work provided fundamental molecular data that will be useful for species identification, conservation, and further studies on the evolution of parrotfishes.

The mitochondrial genome of Heterosentis pseudobagri (Wang & Zhang, 1987) Pichelin & Cribb, 1999 reveals novel aspects of tRNA genes evolution in Acanthocephala

BACKGROUND

Acanthocephala is a clade of obligate endoparasites whose mitochondrial genomes (mitogenomes) and evolution remain relatively poorly understood. Previous studies reported that atp8 is lacking from acanthocephalan mitogenomes, and that tRNA genes often have nonstandard structures. Heterosentis pseudobagri (Arhythmacanthidae) is an acanthocephalan fish endoparasite for which no molecular data are currently available, and biological information is unavailable in the English language. Furthermore, there are currently no mitogenomes available for Arhythmacanthidae.

METHODS

We sequenced its mitogenome and transcriptome, and conducted comparative mitogenomic analyses with almost all available acanthocephalan mitogenomes.

RESULTS

The mitogenome had all genes encoded on the same strand and unique gene order in the dataset. Among the 12 protein-coding genes, several genes were highly divergent and annotated with difficulty. Moreover, several tRNA genes could not be identified automatically, so we had to identify them manually via a detailed comparison with orthologues. As common in acanthocephalans, some tRNAs lacked either the TWC arm or the DHU arm, but in several cases, we annotated tRNA genes only on the basis of the conserved narrow central segment comprising the anticodon, while the flanking 5' and 3' ends did not exhibit any resemblance to orthologues and they could not be folded into a tRNA secondary structure. We corroborated that these are not sequencing artefacts by assembling the mitogenome from transcriptomic data. Although this phenomenon was not observed in previous studies, our comparative analyses revealed the existence of highly divergent tRNAs in multiple acanthocephalan lineages.

CONCLUSIONS

These findings indicate either that multiple tRNA genes are non-functional or that (some) tRNA genes in (some) acanthocephalans might undergo extensive posttranscriptional tRNA processing which restores them to more conventional structures. It is necessary to sequence mitogenomes from yet unrepresented lineages and further explore the unusual patterns of tRNA evolution in Acanthocephala.

Mitochondrial phylogenomics provides conclusive evidence that the family Ancyrocephalidae is deeply paraphyletic

BACKGROUND

Unresolved taxonomic classification and paraphyly pervade the flatworm class Monogenea: the class itself may be paraphyletic and split into Polyopisthocotylea and Monopisthocotylea; there are some indications that the monopisthocotylean order Dactylogyridea may also be paraphyletic; single-gene markers and some morphological traits indicate that the family Ancyrocephalidae is paraphyletic and intertwined with the family Dactylogyridae.

METHODS

To attempt to study the relationships of Ancyrocephalidae and Monopisthocotylea using a phylogenetic marker with high resolution, we sequenced mitochondrial genomes of two fish ectoparasites from the family Dactylogyridae: Dactylogyrus simplex and Dactylogyrus tuba. We conducted phylogenetic analyses using three datasets and three methods. Datasets were ITS1 (nuclear) and nucleotide and amino acid sequences of almost complete mitogenomes of almost all available Monopisthocotylea mitogenomes. Methods were maximum likelihood (IQ-TREE), Bayesian inference (MrBayes) and CAT-GTR (PhyloBayes).

RESULTS

Both mitogenomes exhibited the ancestral gene order for Neodermata, and both were compact, with few and small intergenic regions and many and large overlaps. Gene sequences were remarkably divergent for nominally congeneric species, with only trnI exhibiting an identity value > 80%. Both mitogenomes had exceptionally low A + T base content and AT skews. We found evidence of pervasive compositional heterogeneity in the dataset and indications that base composition biases cause phylogenetic artefacts. All six mitogenomic analyses produced unique topologies, but all nine analyses produced topologies that rendered Ancyrocephalidae deeply paraphyletic. Mitogenomic data consistently resolved the order Capsalidea as nested within the Dactylogyridea.

CONCLUSIONS

The analyses indicate that taxonomic revisions are needed for multiple Polyopisthocotylea lineages, from genera to orders. In combination with previous findings, these results offer conclusive evidence that Ancyrocephalidae is a paraphyletic taxon. The most parsimonious solution to resolve this is to create a catch-all Dactylogyridae sensu lato clade comprising the current Ancyrocephalidae, Ancylodiscoididae, Pseudodactylogyridae and Dactylogyridae families, but the revision needs to be confirmed by another marker with a sufficient resolution.

Integrative taxonomy of the subfamily Orbioninae (Crustacea: Isopoda) based on mitochondrial and nuclear data with evidence that supports Epicaridea as a suborder

Epicaridea is a group of isopods with high morphological diversity, reduction and loss of characters, and strong sexual dimorphism due to their parasitic lifestyles but their systematics is not well understood. Despite the use of nuclear and mitochondrial genes to test the phylogeny of many invertebrate groups, few molecular data from epicarideans are known, especially from the subfamily Orbioninae. Species in this group are obligate penaeoid shrimp parasites and the lack molecular data has hampered studies on the phylogeny of Orbioninae. To rectify this, mitochondrial and nuclear genes of 9 orbionine species are sequenced here. Compared to the isopod ground pattern, the sequences of orbionines seem to be more plastic near the control region and major translocations are located between rrns and cob. A phylogenetic analysis based on three data sets showed strong support for a monophyletic Orbioninae and that Epicaridea should be accepted at the rank of a suborder within Isopoda. The monophyly of Parapenaeon and Orbione is in doubt based on morphological and molecular data. The genus Parapenaeon is revised and a new genus Aparapenaeon is erected for Parapenaeon japonica and three closely related species.

Mitochondrial phylogenomics of Acanthocephala: nucleotide alignments produce long-branch attraction artefacts

BACKGROUND

Classification of the Acanthocephala, a clade of obligate endoparasites, remains unresolved because of insufficiently strong resolution of morphological characters and scarcity of molecular data with a sufficient resolution. Mitochondrial genomes may be a suitable candidate, but they are available for a small number of species and their suitability for the task has not been tested thoroughly.

METHODS

Herein, we sequenced the first mitogenome for the large family Rhadinorhynchidae: Micracanthorhynchina dakusuiensis. These are also the first molecular data generated for this entire genus. We conducted a series of phylogenetic analyses using concatenated nucleotides (NUC) and amino acids (AAs) of all 12 protein-coding genes, three different algorithms, and the entire available acanthocephalan mitogenomic dataset.

RESULTS

We found evidence for strong compositional heterogeneity in the dataset, and Micracanthorhynchina dakusuiensis exhibited a disproportionately long branch in all analyses. This caused a long-branch attraction artefact (LBA) of M. dakusuiensis resolved at the base of the Echinorhynchida clade when the NUC dataset was used in combination with standard phylogenetic algorithms, maximum likelihood (ML) and Bayesian inference (BI). Both the use of the AA dataset (BI-AAs and ML-AAs) and the CAT-GTR model designed for suppression of LBA (CAT-GTR-AAs and CAT-GTR-NUC) at least partially attenuated this LBA artefact. The results support Illiosentidae as the basal radiation of Echinorhynchida and Rhadinorhynchidae forming a clade with Echinorhynchidae and Pomporhynchidae. The questions of the monophyly of Rhadinorhynchidae and its sister lineage remain unresolved. The order Echinorhynchida was paraphyletic in all of our analyses.

CONCLUSIONS

Future studies should take care to attenuate compositional heterogeneity-driven LBA artefacts when applying mitogenomic data to resolve the phylogeny of Acanthocephala.

Integrative taxonomy on the rare sky-island Ligidium species from southwest China (Isopoda, Oniscidea, Ligiidae)

Background

The sky-island Ligidium species fauna in southwest China is poorly known. Before this study, six of the seven sky-island species of the genus were known to be endemic to southwest China. In morphology, Ligidium species are often difficult to identify, and an appraisal of integrative taxonomy is needed.

Results

We integrated morphology and molecular analyses to delimit Ligidium species. Molecular species delimitation based on distance- and evolutionary models recovered seven-candidate lineages from five gene markers (COI, 12S rRNA, 18S rRNA, 28S rRNA and NAK). We also estimated that the species divergences of sky-island Ligidium in southwest China started in late Eocene (40.97 Mya) to middle Miocene (15.19 Mya).

Four new species (L. duospinatum Li, sp. nov., L. acuminatum Li, sp. nov., L. rotundum Li, sp. nov. and L. tridentatum Li, sp. nov.) are described. Morphological confusion among L. denticulatum Shen, 1949, L. inerme Nunomura & Xie, 2000 and L. sichuanense Nunomura, 2002 is clarified by integrative taxonomy.

Conclusion

This work confirms that an integrative approach to Ligidium taxonomy is fundamental for objective classification, and deduced the uplift of Qinghai-Tibetan Plateau in the late Eocene and middle Miocene as one of the principal reasons for the species divergences of sky-island Ligidium in southwest China. We also inferred that sky-island mountains have a huge reserve of higher Ligidium species diversity.

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.

The complete mitochondrial genomes of Paradiplozoon yarkandense and Paradiplozoon homoion confirm that Diplozoidae evolve at an elevated rate

BACKGROUND

Diplozoidae are monogenean (Monogenea: Polyopisthocotylea) fish parasites characterised by a unique life history: two larvae permanently fuse into an X-shaped "Siamese" organism. Taxonomy and phylogeny of Diplozoidae and Polyopisthocotylea remain unresolved due to the unavailability of molecular markers with sufficiently high resolution. Mitogenomes may be a suitable candidate, but there are currently only 12 available for the Polyopisthocotylea (three for Diplozoidae). The only available study of diplozoid mitogenomes found unique base composition patterns and elevated evolution rates in comparison with other Monogenean mitogenomes.

METHODS

To further explore their evolution and generate molecular data for evolutionary studies, we sequenced the complete mitogenomes of two Diplozoidae species, Paradiplozoon homoion and Paradiplozoon yarkandense, and conducted a number of comparative mitogenomic analyses with other polyopisthocotyleans.

RESULTS

We found further evidence that mitogenomes of Diplozoidae evolve at a unique, elevated rate, which was reflected in their exceptionally long branches, large sizes, unique base composition, skews, and very low gene sequence similarity levels between the two newly sequenced species. They also exhibited remarkably large overlaps between some genes. Phylogenetic analysis of Polyopisthocotylea resolved all major taxa as monophyletic, and Mazocraeidea was split into two major clades: (Diplozoidae) + (all four remaining families: Diclidophoridae, Chauhaneidae, Mazocraeidae and Microcotylidae). It also provided further confirmation that the genus Paradiplozoon is paraphyletic and requires a taxonomic revision, so the two species may have to be renamed Indodiplozoon homoion and Diplozoon yarkandense comb. nov.

CONCLUSIONS

Although our findings indicate that mitogenomes may be a promising tool for resolving the phylogeny of Polyopisthocotylea, elevated evolutionary rates of Diplozoidae may cause phylogenetic artefacts, so future studies should pay caution to this problem. Furthermore, as the reason for their elevated evolution remains unknown, Diplozoidae are a remarkably interesting lineage for other types of evolutionary mitogenomic studies.

Sequencing of the Complete Mitochondrial Genome of Pingus sinensis (Spirurina: Quimperiidae): Gene Arrangements and Phylogenetic Implications

Despite several decades of intensive research on spirurine nematodes, molecular data on some of the main lineages are still absent, which makes taxonomic classification insufficiently resolved. In the present study, we sequenced the first complete mitogenome for the family Quimperiidae, belonging to P. sinensis (Spirurina: Quimperiidae), a parasite living in the intestines of snakehead (Ophiocephalus argus). The circular mitogenome is 13,874 bp long, and it contains the standard nematode gene set: 22 transfer RNAs, 2 ribosomal RNAs and 12 protein-coding genes. There are also two long non-coding regions (NCR), in addition to only 8 other intergenic regions, ranging in size from 1 to 58 bp. To investigate its phylogenetic position and study the relationships among other available Spirurina, we performed the phylogenetic analysis using Bayesian inference and maximum likelihood approaches by concatenating the nucleotide sequences of all 36 genes on a dataset containing all available mitogenomes of the suborder Spirurina from NCBI and compared with gene order phylogenies using the MLGO program. Both supported the closer relationship of Ascaridoidea to Seuratoidea than to Spiruroidea. Pingus formed a sister-group with the Cucullanus genus. The results provide a new insights into the relationships within Spirurina.

Building on 150 Years of Knowledge: The Freshwater Isopod Asellus aquaticus as an Integrative Eco-Evolutionary Model System

Interactions between organisms and their environments are central to how biological diversity arises and how natural populations and ecosystems respond to environmental change. These interactions involve processes by which phenotypes are affected by or respond to external conditions (e.g., via phenotypic plasticity or natural selection) as well as processes by which organisms reciprocally interact with the environment (e.g., via eco-evolutionary feedbacks). Organism-environment interactions can be highly dynamic and operate on different hierarchical levels, from genes and phenotypes to populations, communities, and ecosystems. Therefore, the study of organism-environment interactions requires integrative approaches and model systems that are suitable for studies across different hierarchical levels. Here, we introduce the freshwater isopod Asellus aquaticus, a keystone species and an emerging invertebrate model system, as a prime candidate to address fundamental questions in ecology and evolution, and the interfaces therein. We review relevant fields of research that have used A. aquaticus and draft a set of specific scientific questions that can be answered using this species. Specifically, we propose that studies on A. aquaticus can help understanding (i) the influence of host-microbiome interactions on organismal and ecosystem function, (ii) the relevance of biotic interactions in ecosystem processes, and (iii) how ecological conditions and evolutionary forces facilitate phenotypic diversification.

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.

Evolutionary History of Inversions in Directional Mutational Pressures in Crustacean Mitochondrial Genomes: Implications for Evolutionary Studies

Inversions of the origin of replication (ORI) in mitochondrial genomes produce asymmetrical mutational pressures that can cause strong base composition skews. Due to skews often being overlooked, the total number of crustacean lineages that underwent ORI events remains unknown. We analysed skews, cumulative skew plots, conserved sequence motifs, and mitochondrial architecture of all 965 available crustacean mitogenomes (699 unique species). We found indications of an ORI in 159 (22.7%) species, and mapped these to 23 ORI events: 16 identified with confidence and 7 putative (13 newly proposed, and for 5 we improved the resolution). Two ORIs occurred at or above the order level: Isopoda and Copepoda. Shifts in skew plots are not a precise tool for identifying the replication mechanism. We discuss how ORIs can produce mutational bursts in mitogenomes and show how these can interfere with various types of evolutionary studies. Phylogenetic analyses were plagued by artefactual clustering, and ORI lineages exhibited longer branches, a higher number of synonymous substitutions, higher mutational saturation, and higher compositional heterogeneity. ORI events also affected codon usage and protein properties. We discuss how this may have caused erroneous interpretation of data in previous studies that did not account for skew patterns.

Mitogenomics of Cladocera (Branchiopoda): Marked gene order rearrangements and independent predation roots

Cladocera (Crustacea: Branchiopoda) is a key group of invertebrates. Despite a long history of phylogenetic research, relationships within this group remain disputed. We here provide new insights based on 15 new mitochondrial genomes obtained from high-throughput sequencing (HTS) and 40 mitogenomes extracted from published HTS datasets. Together with 25 mitogenomes from GenBank, we generated a matrix of 80 mitogenomes, 44 of them belonging to Cladocera. We also obtained a matrix with 168 nuclear orthologous genes to further assess the phylogenetic result from mitogenomes based on published data and one new HTS data ofLeptodora. Maximum likelihood and Bayesian phylogenetic analyses recovered all Branchiopoda orders as monophyletic and supported a sister-group relationship between Anomopoda and Onychopoda, making the taxon Gymnomera paraphyletic and supporting an independent origin of predatory Haplopoda and Onychopoda. The nuclear phylogeny and topological tests also support Gymnomera as paraphyletic, and the nuclear phylogeny strongly supports a sister-group relationship between Ctenopoda and Haplopoda. We provide a fossil-calibrated time tree, congruent with a Carboniferous origin for Cladocera and a subsequent diversification of the crown group of Anomopoda, Onychopoda, and Ctenopoda, at least in the Triassic. Despite their long evolutionary history, non-Cladoceran Branchiopoda exhibited high mitogenome structural stability. On the other hand, 21 out of 24 gene rearrangements occurred within the relatively younger Cladocera. We found the differential base compositional skewness patterns between Daphnia s.s. and Ctenodaphnia, which might be related to the divergence between these taxa. We also provide evidence to support the recent finding that Spinicaudata possesses mitogenomes with inversed compositional skewness without gene rearrangement. Such a pattern has only been reported in Spinicaudata.

Mitochondrial genomes of two Polydora (Spionidae) species provide further evidence that mitochondrial architecture in the Sedentaria (Annelida) is not conserved

Contrary to the early evidence, which indicated that the mitochondrial architecture in one of the two major annelida clades, Sedentaria, is relatively conserved, a handful of relatively recent studies found evidence that some species exhibit elevated rates of mitochondrial architecture evolution. We sequenced complete mitogenomes belonging to two congeneric shell-boring Spionidae species that cause considerable economic losses in the commercial marine mollusk aquaculture: Polydora brevipalpa and Polydora websteri. The two mitogenomes exhibited very similar architecture. In comparison to other sedentarians, they exhibited some standard features, including all genes encoded on the same strand, uncommon but not unique duplicated trnM gene, as well as a number of unique features. Their comparatively large size (17,673 bp) can be attributed to four non-coding regions larger than 500 bp. We identified an unusually large (putative) overlap of 14 bases between nad2 and cox1 genes in both species. Importantly, the two species exhibited completely rearranged gene orders in comparison to all other available mitogenomes. Along with Serpulidae and Sabellidae, Polydora is the third identified sedentarian lineage that exhibits disproportionally elevated rates of mitogenomic architecture rearrangements. Selection analyses indicate that these three lineages also exhibited relaxed purifying selection pressures.

Disrupted architecture and fast evolution of the mitochondrial genome of Argeia pugettensis (Isopoda): implications for speciation and fitness

BACKGROUND

Argeia pugettensis is an isopod species that parasitizes other crustaceans. Its huge native geographic range spans the Pacific from China to California, but molecular data are available only for a handful of specimens from North-American populations. We sequenced and characterised the complete mitogenome of a specimen collected in the Yellow Sea.

RESULTS

It exhibited a barcode (cox1) similarity level of only 87-89% with North-American populations, which is unusually low for conspecifics. Its mitogenome is among the largest in isopods (≈16.5 Kbp), mostly due to a large duplicated palindromic genomic segment (2 Kbp) comprising three genes. However, it lost a segment comprising three genes, nad4L-trnP-nad6, and many genes exhibited highly divergent sequences in comparison to isopod orthologues, including numerous mutations, deletions and insertions. Phylogenetic and selection analyses corroborated that this is one of the handful of most rapidly evolving available isopod mitogenomes, and that it evolves under highly relaxed selection constraints (as opposed to positive selection). However, its nuclear 18S gene is highly conserved, which suggests that rapid evolution is limited to its mitochondrial genome. The cox1 sequence analysis indicates that elevated mitogenomic evolutionary rates are not shared by North-American conspecifics, which suggests a breakdown of cox1 barcoding in this species.

CONCLUSIONS

A highly architecturally disrupted mitogenome and decoupling of mitochondrial and nuclear rates would normally be expected to have strong negative impacts on the fitness of the organism, so the existence of this lineage is a puzzling evolutionary question. Additional studies are needed to assess the phylogenetic breadth of this disrupted mitochondrial architecture and its impact on fitness.

Comparison of the complete mitochondrial genome of Phyllophorus liuwutiensis (Echinodermata: Holothuroidea: Phyllophoridae) to that of other sea cucumbers

Sea cucumber species are abundant (>1400 species) and widely distributed globally. mtDNA sequencing is frequently used to identify the phylogenetic and evolutionary relationships among species. However, there are no reports on the mitochondrial genome of Phyllophorus liuwutiensis. Here, we performed mtDNA sequencing of P. liuwutiensis to examine its phylogenetic relationships with other echinoderms. Its mitochondrial genome (15 969 bp) contains 37 coding genes, including 13 protein‐coding genes, 22 tRNA genes and 2 rRNA genes. Except for one protein‐coding gene (nad6) and five tRNA genes encoded on the negative strand, all other genes were encoded on the positive strand. The mitochondrial bases of P. liuwutiensis were composed of 29.55% T, 22.16% C, 35.64% A and 12.64% G. The putative control region was 703 bp in length. Seven overlapping regions (1–10 bp) were found. The noncoding region between the genes ranged from 1 to 130 bp in length. One putative control region has been found in the P. liuwutiensis mitogenome. All of the tRNA genes were predicted to fold into a cloverleaf structure. In addition, we compared the gene arrangements of six echinoderms, revealing that the gene order of P. liuwutiensis was a new arrangement.

Architectural instability, inverted skews and mitochondrial phylogenomics of Isopoda: outgroup choice affects the long-branch attraction artefacts

The majority strand of mitochondrial genomes of crustaceans usually exhibits negative GC skews. Most isopods exhibit an inversed strand asymmetry, believed to be a consequence of an inversion of the replication origin (ROI). Recently, we proposed that an additional ROI event in the common ancestor of Cymothoidae and Corallanidae families resulted in a double-inverted skew (negative GC), and that taxa with homoplastic skews cluster together in phylogenetic analyses (long-branch attraction, LBA). Herein, we further explore these hypotheses, for which we sequenced the mitogenome of Asotana magnifica (Cymothoidae), and tested whether our conclusions were biased by poor taxon sampling and inclusion of outgroups. (1) The new mitogenome also exhibits a double-inverted skew, which supports the hypothesis of an additional ROI event in the common ancestor of Cymothoidae and Corallanidae families. (2) It exhibits a unique gene order, which corroborates that isopods possess exceptionally destabilized mitogenomic architecture. (3) Improved taxonomic sampling failed to resolve skew-driven phylogenetic artefacts. (4) The use of a single outgroup exacerbated the LBA, whereas both the use of a large number of outgroups and complete exclusion of outgroups ameliorated it.

Mitochondrial genomes and 28S rDNA contradict the proposed obsoletion of the order Tetraonchidea (Platyhelminthes: Monogenea)

Due to the incongruence of morphology-based hypotheses and scarcity of molecular data, validity of the order Tetraonchidea remains contentious. The only complete mitogenome currently available for the entire order is that of Paratetraonchoides inermis (Tetraonchoididae). To study the phylogeny of Tetraonchidea from mitogenomic perspective, we sequenced the first mitogenome for the family Tetraonchidae: Tetraonchus monenteron (Tetraonchidea). To get a nuclear-data perspective, we also sequenced nuclear 28S rDNA gene of both species. The mitogenome of T. monenteron does not have high A + T content, nor tRNA pseudo-genes, both of which were unique features reported in P. inermis. However, T. monenteron exhibits a unique gene order, with a large number of tRNA rearrangements in comparison to P. inermis and other monogeneans. Phylogenetic analyses conducted using Bayesian inference and maximum likelihood methods, complemented with partitioning, consistently support the sister-group relationship of T. monenteron (Tetraonchidae) and P. inermis (Tetraonchoididae). This is also partially supported by the 28S rDNA data and two morphologic apomorphies. This close relationship of Tetraonchidae and Tetraonchoididae challenges the latest major morphology-based classification, which proposed obsoletion of the Tetraonchidea order, and grouped Tetraonchoididae into the Gyrodactylidea clade. The validity of this order shall have to be further confirmed with more data.

Disentangling the interplay of positive and negative selection forces that shaped mitochondrial genomes of Gammarus pisinnus and Gammarus lacustris

We hypothesized that the mitogenome of Gammarus lacustris (GL), native to the Qinghai-Tibet Plateau, might exhibit genetic adaptations to the extreme environmental conditions associated with high altitudes (greater than 3000 m). To test this, we also sequenced the mitogenome of Gammarus pisinnus (GP), whose native range is close to the Tibetan plateau, but at a much lower altitude (200-1500 m). The two mitogenomes exhibited conserved mitochondrial architecture, but low identity between genes (55% atp8 to 76.1% cox1). Standard (homogeneous) phylogenetic models resolved Gammaridae as paraphyletic, but 'heterogeneous' CAT-GTR model as monophyletic. In indirect support of our working hypothesis, GL, GP and Gammarus fossarum exhibit evidence of episodic diversifying selection within the studied Gammaroidea dataset. The mitogenome of GL generally evolves under a strong purifying selection, whereas GP evolves under directional (especially pronounced in atp8) and/or relaxed selection. This is surprising, as GP does not inhabit a unique ecological niche compared to other gammarids. We propose that this rapid evolution of the GP mitogenome may be a reflection of its relatively recent speciation and heightened non-adaptive (putatively metabolic rate-driven) mutational pressures. To test these hypotheses, we urge sequencing mitogenomes of remaining Gammarus species populating the same geographical range as GP.

Complete mitochondrial genomes from transcriptomes: assessing pros and cons of data mining for assembling new mitogenomes

Thousands of eukaryotes transcriptomes have been generated, mainly to investigate nuclear genes expression, and the amount of available data is constantly increasing. A neglected but promising use of this large amount of data is to assemble organelle genomes. To assess the reliability of this approach, we attempted to reconstruct complete mitochondrial genomes from RNA-Seq experiments of Reticulitermes termite species, for which transcriptomes and conspecific mitogenomes are available. We successfully assembled complete molecules, although a few gaps corresponding to tRNAs had to be filled manually. We also reconstructed, for the first time, the mitogenome of Reticulitermes banyulensis. The accuracy and completeness of mitogenomes reconstruction appeared independent from transcriptome size, read length and sequencing design (single/paired end), and using reference genomes from congeneric or intra-familial taxa did not significantly affect the assembly. Transcriptome-derived mitogenomes were found highly similar to the conspecific ones obtained from genome sequencing (nucleotide divergence ranging from 0% to 3.5%) and yielded a congruent phylogenetic tree. Reads from contaminants and nuclear transcripts, although slowing down the process, did not result in chimeric sequence reconstruction. We suggest that the described approach has the potential to increase the number of available mitogenomes by exploiting the rapidly increasing number of transcriptomes.

Mitochondrial Genomes of Two Thaparocleidus Species (Platyhelminthes: Monogenea) Reveal the First rRNA Gene Rearrangement among the Neodermata

Phylogenetic framework for the closely related Ancylodiscoidinae and Ancyrocephalinae subfamilies remains contentious. As this issue was never studied using a large molecular marker, we sequenced the first two Ancylodiscoidinae mitogenomes: Thaparocleidus asoti and Thaparocleidus varicus. Both mitogenomes had two non-coding regions (NCRs) that contained a number of repetitive hairpin-forming elements (RHE). Due to these, the mitogenome of T. asoti (16,074 bp) is the longest among the Monogenea; especially large is its major NCR, with 3500 bp, approximately 1500 bp of which could not be sequenced (thus, the total mitogenome size is ≈ 17,600 bp). Although RHEs have been identified in other monopisthocotyleans, they appear to be independently derived in different taxa. The presence of RHEs may have contributed to the high gene order rearrangement rate observed in the two mitogenomes, including the first report of a transposition of rRNA genes within the Neodermata. Phylogenetic analyses using mitogenomic dataset produced Dactylogyrinae embedded within the Ancyrocephalinae (paraphyly), whereas Ancylodiscoidinae formed a sister-group with them. This was also supported by the gene order analysis. 28S rDNA dataset produced polyphyletic Dactylogyridae and Ancyrocephalinae. The phylogeny of the two subfamilies shall have to be further evaluated with more data.