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

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.

Strong mitonuclear discordance in the phylogeny of Neodermata and evolutionary rates of Polyopisthocotylea

The genomic evolution of Polyopisthocotylea remains poorly understood in comparison to the remaining three classes of Neodermata: Monopisthocotylea, Cestoda, and Trematoda. Moreover, the evolutionary sequence of major events in the phylogeny of Neodermata remains unresolved. Herein we sequenced the mitogenome and transcriptome of the polyopisthocotylean Diplorchis sp., and conducted comparative evolutionary analyses using nuclear (nDNA) and mitochondrial (mtDNA) genomic datasets of Neodermata. We found strong mitonuclear discordance in the phylogeny of Neodermata. Polyopisthocotylea exhibited striking mitonuclear discordance in relative evolutionary rates: the fastest-evolving mtDNA in Neodermata and a comparatively slowly-evolving nDNA genome. This was largely attributable to its very long stem branch in mtDNA topologies, not exhibited by the nDNA data. We found indications that the fast evolution of mitochondrial genomes of Polyopisthocotylea may be driven both by relaxed purifying selection pressures and elevated levels of directional selection. We identified mitochondria-associated genes encoded in the nuclear genome: they exhibited unique evolutionary rates, but not correlated with the evolutionary rate of mtDNA, and there is no evidence for compensatory evolution (they evolved slower than the rest of the genome). Finally, there appears to exist an exceptionally large (≈6.3 kb) nuclear mitochondrial DNA segment (numt) in the nuclear genome of newly sequenced Diplorchis sp. A 3'-end segment of the 16S rRNA gene encoded by the numt was expressed, suggesting that this gene acquired novel, regulatory functions after the transposition to the nuclear genome. In conclusion, Polyopisthocotylea appears to be the lineage with the fastest-evolving mtDNA sequences among all of Bilateria, but most of the substitutions were accumulated deep in the evolutionary history of this lineage. As the nuclear genome does not exhibit a similar pattern, the circumstances underpinning this evolutionary phenomenon remain a mystery.

Complete Mitogenome sequencing of the fish louse Argulus japonicus (Crustacea: Branchiura): Comparative analyses and phylogenetic implications

The fish louse Argulus japonicus, a branchiuran crustacean of the Argulidae family, is attracting increasing attention because of its parasitic tendencies and significant health threats to global fish farming. The mitogenomes can yield a foundation for studying epidemiology, genetic diversity, and molecular ecology and therefore may be used to assist in the surveillance and control of A. japonicus. In this study, we sequenced and assembled the complete mitogenome of A. japonicus to shed light on its genetic and evolutionary blueprint. Our investigation indicated that the 15,045-bp circular genome of A. japonicus encodes 13 protein-coding genes (PCGs), 22 transfer RNAs (tRNAs), and 2 ribosomal RNAs (rRNAs) with significant AT and GC skews. Comparative genomics provided an evolutionary scenario for the genetic diversity of 13 PCGs: all were under purifying selection, with cox1 and nad6 having the lowest and highest evolutionary rates, respectively. Genome-wide phylogenetic trees established a close relationship between species of the families Argulidae (Arguloida) and Armilliferidae (Porocephalida) within Crustacea, and further, A. japonicus and Argulus americanus were determined to be more closely related to each other than to others within the family Argulidae. Single PCG-based phylogenies supported nad1 and nad6 as the best genetic markers for evolutionary and phylogenetic studies for branchiuran crustaceans due to their similar phylogenetic topologies with those of genome-based phylogenetic analyses. To sum up, these comprehensive mitogenomic data of A. japonicus and related species refine valuable marker resources and should contribute to molecular diagnostic methods, epidemiological investigations, and ecological studies of the fish ectoparasites in Crustacea.

Comprehensive molecular characterisation of the complete mitogenome of Ergasilus tumidus and phylogenetic relationships of Copepoda inferred from mitogenomes

Although parasitic copepods of the genus Ergasilus von Nordmann, 1832 are globally distributed parasites of fish, their phylogenetic relationships with other Copepoda are not clear, and the characteristics of their mitochondrial genomes (mitogenomes) are not thoroughly understood. The objective of this study was to address these knowledge gaps by sequencing the complete mitogenome of Ergasilus tumidus Markevich, 1940. The complete mitogenome (GenBank Acc. No. OQ596537) was 14,431 bp long and it comprised 13 protein-coding genes (PCGs), 22 tRNAs, two tRNAs, and two control regions (CRs). Phylogenetic analyses, conducted using concatenated nucleotide and amino acid sequences of 13 protein-coding genes, produced two partially incongruent topologies. While the order Calanoida was consistently resolved as the sister lineage to the other three orders, topological instability was observed in the relationships of the orders Cyclopoida, Siphonostomatoida and Harpacticoida. Siphonostomatoida clustered with Cyclopoida in the nucleotide-based phylogeny, but with Harpacticoida in the amino acid-based phylogeny. The latter topology conforms to the widely accepted relationships, but we speculate that the former topology is more likely to be the correct one. Our study provides a complete mitogenome sequence of E. tumidus, which helps us better understand the molecular evolution of the genus Ergasilus. Additionally, we suggest a different perspective on the controversial phylogenetic relationships among Siphonostomatoida, Cyclopoida and Harpacticoida, diverging from previously accepted views.

Mitochondrial Genome Evolution in Annelida-A Systematic Study on Conservative and Variable Gene Orders and the Factors Influencing its Evolution

The mitochondrial genomes of Bilateria are relatively conserved in their protein-coding, rRNA, and tRNA gene complement, but the order of these genes can range from very conserved to very variable depending on the taxon. The supposedly conserved gene order of Annelida has been used to support the placement of some taxa within Annelida. Recently, authors have cast doubts on the conserved nature of the annelid gene order. Various factors may influence gene order variability including, among others, increased substitution rates, base composition differences, structure of noncoding regions, parasitism, living in extreme habitats, short generation times, and biomineralization. However, these analyses were neither done systematically nor based on well-established reference trees. Several focused on only a few of these factors and biological factors were usually explored ad-hoc without rigorous testing or correlation analyses. Herein, we investigated the variability and evolution of the annelid gene order and the factors that potentially influenced its evolution, using a comprehensive and systematic approach. The analyses were based on 170 genomes, including 33 previously unrepresented species. Our analyses included 706 different molecular properties, 20 life-history and ecological traits, and a reference tree corresponding to recent improvements concerning the annelid tree. The results showed that the gene order with and without tRNAs is generally conserved. However, individual taxa exhibit higher degrees of variability. None of the analyzed life-history and ecological traits explained the observed variability across mitochondrial gene orders. In contrast, the combination and interaction of the best-predicting factors for substitution rate and base composition explained up to 30% of the observed variability. Accordingly, correlation analyses of different molecular properties of the mitochondrial genomes showed an intricate network of direct and indirect correlations between the different molecular factors. Hence, gene order evolution seems to be driven by molecular evolutionary aspects rather than by life history or ecology. On the other hand, variability of the gene order does not predict if a taxon is difficult to place in molecular phylogenetic reconstructions using sequence data or not. We also discuss the molecular properties of annelid mitochondrial genomes considering canonical views on gene evolution and potential reasons why the canonical views do not always fit to the observed patterns without making some adjustments. [Annelida; compositional biases; ecology; gene order; life history; macroevolution; mitochondrial genomes; substitution rates.].

Evolutionary Rates, Divergence Rates, and Performance of Individual Mitochondrial Genes Based on Phylogenetic Analysis of Copepoda

Copepoda is a large and diverse group of crustaceans, which is widely distributed worldwide. It encompasses roughly 9 orders, whose phylogeny remains unresolved. We sequenced the complete mitochondrial genome (mitogenome) of Sinergasilus major (Markevich, 1940) and used it to explore the phylogeny and mitogenomic evolution of Copepoda. The mitogenome of S. major (14,588 bp) encodes the standard 37 genes as well as a putative control region, and molecular features are highly conserved compared to other Copepoda mitogenomes. Comparative analyses indicated that the nad2 gene has relatively high nucleotide diversity and evolutionary rate, as well as the largest amount of phylogenetic information. These results indicate that nad2 may be a better marker to investigate phylogenetic relationships among closely related species in Copepoda than the commonly used cox1 gene. The sister-group relationship of Siphonostomatoida and Cyclopoida was recovered with strong support in our study. The only topological ambiguity was found within Cyclopoida, which might be caused by the rapid evolution and sparse taxon sampling of this lineage. More taxa and genes should be used to reconstruct the Copepoda phylogeny in the future.

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.

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.

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.

Sequencing and phylogenomics of the complete mitochondrial genome of Allodiplogaster sp. (Rhabditida: Diplogasteridae): A new gene order and its phylogenetic implications

Gene order has been utilized as a phylogenetic signal for many taxa. However, its phylogenetic performance has not been evaluated in Nematoda. As there is only one nematode mitogenome available to date, in the Diplogasteridae family, we sequenced the mitogenome of Allodiplogaster sp. and constructed a phylogeny for Nematoda using this updated mitogenome dataset. We then compared this phylogeny to one constructed using gene order. The complete mitochondrial genome of Allodiplogaster sp. was 13,953 bp in size and included 22 tRNAs, two rRNAs, and 12 protein-coding genes. To assess how Allodiplogaster sp. is related to other nematode species, we used Bayesian inference and maximum likelihood algorithms to construct phylogenetic trees of the Nematoda. We found that: 1) The target species Allodiplogaster sp. is closely related to Allodiplogaster sudhausi. The topology of the mitogenome based phylogeny was nearly identical to previous phylogenies created using 18S rRNA data, except for the placement of the Strongyloididae family. 2) The maximum likelihood tree constructed using gene order was roughly consistent with the mitogenome-based tree at the family level, but not at the species level. 3) Protein-coding genes were ordered differently in Allodiplogaster sp. versus Allodiplogaster sudhausi; this represents the first report of such a reordering in the class Chromadorea in our study. Our study confirms that gene order represents useful phylogenetic information for the Nematoda: the maximum likelihood tree based on gene order provided additional support for the nematode phylogeny constructed using molecular data.

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.

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.