Evolutionary history of inversions in the direction of architecture-driven mutational pressures in crustacean mitochondrial genomes

Inversions of the origin of replication (ORI) of mitochondrial genomes produce asymmetrical mutational pressures that can cause artefactual clustering in phylogenetic analyses. It is therefore an absolute prerequisite for all molecular evolution studies that use mitochondrial data to account for ORI events in the evolutionary history of their dataset. The number of ORI events in crustaceans remains unknown; several studies reported ORI events in some crustacean lineages on the basis of fully inversed (e.g. negative vs. positive) GC skew patterns, but studies of isolated lineages could have easily overlooked ORI events that produced merely a reduction in the skew magnitude. In this study, we used a comprehensive taxonomic approach to systematically study the evolutionary history of ORI events in crustaceans using all available mitogenomes and combining signals from lineage-specific skew magnitude and direction (+ or -), cumulative skew diagrams, and gene rearrangements. We inferred 24 putative ORI events (14 of which have not been proposed before): 17 with relative confidence, and 7 speculative. Most of these were located at lower taxonomic levels, but there are indications of ORIs that occurred at or above the order-level: Copepoda, Isopoda, and putatively in Branchiopoda and Poecilostomatida+Cyclopoida. Several putative ORI events did not result in fully inversed skews. In many lineages skew plots were not informative for the prediction of replication origin and direction of mutational pressures, but inversions of the mitogenome fragment comprising the ancestral CR (rrnS-CR-trnI) were rather good predictors of skew inversions. We also found that skew plots can be a useful tool to indirectly infer the relative strengths of mutational/purifying pressures in some crustacean lineages: when purifying pressures outweigh mutational, GC skew plots are strongly affected by the strand distribution of genes, and when mutational > purifying, GC skew plots can be even completely (apparently) unaffected by the strand distribution of genes. This observation has very important repercussions for phylogenetic and evolutionary studies, as it implies that not only the relatively rare ORI events, but also much more common gene strand switches and same-strand rearrangements can produce mutational bursts, which in turn affect phylogenetic and evolutionary analyses. We argue that such compositional biases may produce misleading signals not only in phylogenetic but also in other types of evolutionary analyses (dN/dS ratios, codon usage bias, base composition, branch length comparison, etc.), and discuss several such examples. Therefore, all studies aiming to study the evolution of mtDNA sequences should pay close attention to architectural rearrangements.