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

Characterization of the complete mitochondrial genome of Desmaulus extinctorium (Littorinimorpha, Calyptraeoidea, Calyptraeidae) and molecular phylogeny of Littorinimorpha

For the purpose of determining the placement of Calyptraeidae within the Littorinimorpha, we hereby furnish a thorough analysis of the mitochondrial genome (mitogenome) sequence of Desmaulus extinctorium. This mitogenome spans 16,605 base pairs and encompasses the entire set of 37 genes, including 13 PCGs, 22 tRNAs and two rRNAs, with an evident AT bias. Notably, tRNASer1 and tRNASer2 lack dihydrouracil (DHU) arms, resulting in an inability to form a secondary structure. Similarly, tRNAAla lacks a TΨC arm, rendering it incapable of forming a secondary structure. In contrast, the remaining tRNAs demonstrate a characteristic secondary structure reminiscent of a cloverleaf. A comparison with ancestral gastropods reveals distinct differences in three gene clusters (or genes), encompassing 15 tRNAs and eight PCGs. Notably, inversions and translocations represent the major types of rearrangements observed in D. extinctorium. Phylogenetic analysis demonstrates robust support for a monophyletic grouping of all Littorinimorpha species, with D. extinctorium representing a distinct Calyptraeoidea clade. In summary, this investigation provides the first complete mitochondrial dataset for a species of the Calyptraeidae, thus providing novel insights into the phylogenetic relationships within the Littorinimorpha.

Phylogenetic relationships and adaptation in deep-sea carideans revealed by mitogenomes

The deep-sea environment is characterized by extreme and inhospitable conditions, including oxygen depletion, low temperatures, high pressure, absence of light, and limited food availability. Mitochondria and mitogenomes play a crudial role in aerobic respiration to generate energy for eukaryotes. Here, using the Illumina Hiseq 4000 platform, we performed mitogenome sequencing for five deep-sea caridean species: Lebbeus shinkaiae, Lebbeus Formosus, Glyphocrangon regalis, Heterocarpus dorsalis, and Heterocarpus laevigatus, and five deep-sea caridean mitogenomes were assembled and identified. Each of the five mitogenomes contained 13 protein-coding genes, 2 rRNAs and 22 tRNAs. Specific elements, such as tandem repeats and AT-rich sequences, were observed in the control regions of Lebbeus formosus and Lebbeus shinkaiae, potentially take a role in regulating mitochondrial genome replication and transcription. The gene order of all obtained mitogenomes follows caridean ancestral type organization. Phylogenetic analysis shows a robustly supported phylogenetic tree for the infraorder Caridea. The monophyly of the families included in this study was strongly supported. This study supports the monophyly of Oplophoroidea, but rejects the monophyletic status of Nematocarcinoidea, Crangonoidea, and Alpheoidea. At the genus level, Plesionika is polyphyletic and Rimicaris is paraphyletic in our analysis. Furthermore, Paralebbeus may be considered invalid and synonymous with Lebbeus. Positive selection analysis reveals evidence for adaptive changes in the mitogenome of different deep-sea caridean lineages. Nine residues located in cox1, cox3, atp6, nad1, nad2, nad4, nad5, nad6 and cytb were determined to have undergone positive selection. Mitogenome of different deep-sea lineages experienced different positive selection, and the lineage represented by Alvinocarididae living in deep-sea hydrothermal vents experienced the strongest positive selection. This study provides valuable insights into the adaptive evolution of deep-sea shrimps at the mitochondrial, highlighting the mitogenomic strategy that contribute to their unique adaptations in the deep-sea environment.

Phylomitogenomics elucidates the evolution of symbiosis in Thoracotremata (Decapoda: Cryptochiridae, Pinnotheridae, Varunidae)

Background

Thoracotremata belong to the large group of "true" crabs (infraorder Brachyura), and they exhibit a wide range of physiological and morphological adaptations to living in terrestrial, freshwater and marine habitats. Moreover, the clade comprises various symbiotic taxa (Aphanodactylidae, Cryptochiridae, Pinnotheridae, some Varunidae) that are specialised in living with invertebrate hosts, but the evolutionary history of these symbiotic crabs is still partially unresolved.

Methods

Here we assembled and characterised the complete mitochondrial genomes (hereafter mitogenomes) of three gall crab species (Cryptochiridae): Kroppcarcinus siderastreicola, Opecarcinus hypostegus and Troglocarcinus corallicola. A phylogenetic tree of the Thoracotremata was reconstructed using 13 protein-coding genes and two ribosomal RNA genes retrieved from three new gall crab mitogenomes and a further 72 available thoracotreme mitogenomes. Furthermore, we applied a comparative analysis to characterise mitochondrial gene order arrangement, and performed a selection analysis to test for selective pressure of the protein-coding genes in symbiotic Cryptochiridae, Pinnotheridae, and Varunidae (Asthenognathus inaequipes and Tritodynamia horvathi).

Results

The results of the phylogenetic reconstruction confirm the monophyly of Cryptochiridae, which clustered separately from the Pinnotheridae. The latter clustered at the base of the tree with robust branch values. The symbiotic varunids A. inaequipes and T. horvathi clustered together in a clade with free-living Varunidae species, highlighting that symbiosis in the Thoracotremata evolved independently on multiple occasions. Different gene orders were detected in symbionts and free-living species when compared with the ancestral brachyuran gene order. Lastly, the selective pressure analysis detected two positively selected sites in the nad6 gene of Cryptochiridae, but the evidence for positive selection in Pinnotheridae and A. inaequipes and T. horvathi was weak. Adaptive evolution of mitochondrial protein-coding genes is perhaps related to the presumably higher energetic demands of a symbiotic lifestyle.

The mitochondrial genome of the pentastome parasite Raillietiella orientalis Hett, 1915 (Raillietiellida; Raillietiellidae) with notes on its phylogenetic position

In this study we sequenced and annotated the complete mitochondrial genome of the invasive reptile parasite Raillietiella orientalis using Illumina DNA sequencing. The length of the mitogenome was 15,320 bp and had a GC content of 33.1%. The mitogenome contained 13 protein-coding genes, two ribosomal RNA genes, and 22 tRNA genes, the order of which was diagnostic of Pancrustacean mitogenomes. A phylogenetic tree constructed from the 13 protein-coding genes of R. orientalis and 26 other Pancrustacean mitogenomes supported the placement of R. orientalis as part of the monophyletic subclass Pentastomida within the Maxillopoda and sister to the subclass Branchiura.

Molecular Phylogeny of Thoracotreme Crabs Including Nine Newly Determined Mitochondrial Genomes

Mitochondrial genomes are used widely for the molecular phylogenetic analysis of animals. Although phylogenetic analyses based on the mitogenomes of brachyurans often yield well-resolved phylogenies, most interfamilial phylogenetic relationships in Thoracotremata remain unclear. We determined nine new mitogenomes of Thoracotremata, including mitogenomes of Camptandriidae (Deiratonotus japonicus), Dotillidae (Ilyoplax integra, Ilyoplax pusilla, and Tmethypocoelis choreutes), Macrophthalmidae (Ilyograpsus nodulosus), Pinnotheridae (Arcotheres sp. and Indopinnixa haematosticta), Plagusiidae (Guinusia dentipes), and Percnidae (Percnon planissimum). Interestingly, Percnon planissimum (Percnidae) was found to possess ≥ 19 repeated sequences in the control region. The gene orders of Il. nodulosus, Arcotheres sp., and In. haematosticta were revealed to be unique among thoracotreme crabs. Although the results of Bayesian and maximum likelihood (ML) phylogenetic analyses of three datasets were incongruent, highly supported clades (PP ≥ 0.99 or BS ≥ 99%) were not contradictory among the analyses. All analyses suggested the paraphyly of Grapsoidea and Ocypodoidea, corroborating the findings of previous studies based on molecular phylogenies of thoracotreme crabs. The phylogenetic positions of symbiotic thoracotreme crabs, Pinnotheridae and Cryptochiridae, were highly supported (Pinnotheridae + Ocypodidae and Cryptochiridae + Grapsidae, respectively) for the Bayesian analyses but not for the ML analyses. Analyses of more thoracotreme species' mitogenome sequences in additional studies will further strengthen the framework for thoracotreme evolution.

Novel gene re-arrangement in the mitochondrial genome of Pisidiaserratifrons (Anomura, Galatheoidea, Porcellanidae) and phylogenetic associations in Anomura

To improve the taxonomy and systematics of Porcellanidae within the evolution of Anomura, we describe the complete mitochondrial genomes (mitogenomes) sequence of Pisidiaserratifrons, which is 15,344 bp in size, contains the entire set of 37 genes and has an AT-rich region. Compared with the pancrustacean ground pattern, at least five gene clusters (or genes) are significantly different with the typical genes, involving eleven tRNA genes and four PCGs and the tandem duplication/random loss and recombination models were used to explain the observed large-scale gene re-arrangements. The phylogenetic results showed that all Porcellanidae species clustered together as a group with well nodal support. Most Anomura superfamilies were found to be monophyletic, except Paguroidea. Divergence time estimation implies that the age of Anomura is over 225 MYA, dating back to at least the late Triassic. Most of the extant superfamilies and families arose during the late Cretaceous to early Tertiary. In general, the results obtained in this study will contribute to a better understanding of gene re-arrangements in Porcellanidae mitogenomes and provide new insights into the phylogeny of Anomura.

The Complete Mitochondrial Genome of the Hermit Crab Diogenes edwardsii (Anomura: Diogenidae) and Phylogenetic Relationships within Infraorder Anomura

A complete mitochondrial genome (mitogenome) can provide important information for gene rearrangement, molecular evolution and phylogenetic analysis. Currently, only a few mitogenomes of hermit crabs (superfamily Paguridae) in the infraorder Anomura have been reported. This study reports the first complete mitogenome of the hermit crab Diogenes edwardsii assembled using high-throughput sequencing. The mitogenome of Diogenes edwardsii is 19,858 bp in length and comprises 13 protein-coding genes, 2 ribosomal RNA genes, and 22 transfer RNA genes. There are 28 and six genes observed on the heavy and light strands, respectively. The genome composition was highly A + T biased (72.16%), and exhibited a negative AT-skew (-0.110) and positive GC-skew (0.233). Phylogenetic analyses based on the nucleotide dataset of 16 Anomura species indicated that D. edwardsii was closest related to Clibanarius infraspinatus in the same family, Diogenidae. Positive selection analysis showed that two residues located in cox1 and cox2 were identified as positively selected sites with high BEB value (>95%), indicating that these two genes are under positive selection pressure. This is the first complete mitogenome of the genus Diogenes, and this finding helps us to represent a new genomic resource for hermit crab species and provide data for further evolutionary status of Diogenidae in Anomura.

Checklist of decapod crustaceans from the coast of the São Paulo state (Brazil) supported by integrative molecular and morphological data: V. Dendrobranchiata and Pleocyemata [Achelata, Astacidea, Axiidea, Caridea (Alpheoidea and Processoidea excluded), Gebiidea, Stenopodidea]

This checklist is the fifth and last compilation on the decapod crustaceans reported to São Paulo (Brazil) coastal area, resulting from long-term multidisciplinary projects, which combined morphological analyses and molecular techniques. The current research includes 75 decapod species, herein referred as shrimps/lobsters-like (shrimps, ghost-shrimps, lobsters, and related groups), reported to São Paulo coastal area. These species occur in marine, estuarine, and amphidromous habitats and are classified into 21 families as follow: Aristeidae (2 spp., 2 genera), Atyidae (4 spp., 2 genera), Axianassidae (1 sp., 1 genus), Callianassidae (1 sp., 1 genus), Callichiridae (6 spp., 4 genera), Crangonidae (1 sp., 1 genus), Glyphocrangonidae (1 sp., 1 genus), Luciferidae (2 spp., 2 genera), Nephropidae (4 spp., 2 genera), Palaemonidae (15 spp., 9 genera), Palinuridae (2 spp., 1 genus), Pandalidae (1 sp., 1 genus), Pasiphaeidae (1 sp., 1 genus), Penaeidae (10 spp., 6 genera), Sergestidae (3 spp., 3 genera), Sicyoniidae (4 spp., 1 genus), Scyllaridae (5 spp., 3 genera), Solenoceridae (4 spp., 3 genera), Stenopodidae (2 spp., 1 genus), and Upogebiidae (6 spp., 1 genus). We generated new sequences of cytochrome oxidase subunit I (barcode region) and 16S genes (51 and 54, respectively) of 54 species. Our examination concluded that 75 shrimps/lobsters-like species are reported to the São Paulo coast. We excluded Leander tenuicornis (Palaemonidae), Penaeus setiferus (Penaeidae), Philocheras gorei (Crangonidae), and Rhynchocinetes typus (Rhynchocinetidae) from this list.  

Exploring mitogenome evolution in Branchiopoda (Crustacea) lineages reveals gene order rearrangements in Cladocera

The class Branchiopoda, whose origin dates back to Cambrian, includes ~ 1200 species which mainly occupy freshwater habitats. The phylogeny and systematics of the class have been debated for long time, until recent phylogenomic analyses allowed to better clarify the relationships among major clades. Based on these data, the clade Anostraca (fairy and brine shrimps) is sister to all other branchiopods, and the Notostraca (tadpole shrimps) results as sister group to Diplostraca, which includes Laevicaudata + Spinicaudata (clam shrimps) and Cladoceromorpha (water fleas + Cyclestherida). In the present analysis, thanks to an increased taxon sampling, a complex picture emerges. Most of the analyzed mitogenomes show the Pancrustacea gene order while in several other taxa they are found rearranged. These rearrangements, though, occur unevenly among taxa, most of them being found in Cladocera, and their taxonomic distribution does not agree with the phylogeny. Our data also seems to suggest the possibility of potentially homoplastic, alternative gene order within Daphniidae.

Mitogenome of a stink worm (Annelida: Travisiidae) includes degenerate group II intron that is also found in five congeneric species

Mitogenomes are useful for inferring phylogenetic relationships between organisms. Although the mitogenomes of Annelida, one of the most morphologically and ecologically diverse metazoan groups have been well sequenced, those of several families remain unexamined. This study determined the first mitogenome from the family Travisiidae (Travisia sanrikuensis), analyzed its mitogenomic features, and reconstructed a phylogeny of Sedentaria. The monophyly of the Terebellida + Arenicolida + Travisiidae clade is supported by molecular phylogenetic analysis. The placement of Travisiidae is unclear because of the lack of mitogenomes from closely related lineages. An unexpected intron appeared within the cox1 gene of T. sanrikuensis and in the same positions of five undescribed Travisia spp. Although the introns are shorter (790–1386 bp) than other group II introns, they can be considered degenerate group II introns due to type II intron maturase open reading frames, found in two of the examined species, and motifs characteristic of group II introns. This is likely the first known case in metazoans where mitochondrial group II introns obtained by a common ancestor are conserved in several descendants. Insufficient evolutionary time for intron loss in Travisiidae, or undetermined mechanisms may have helped maintain the degenerate introns.

The mitochondrial genome of Grapsus albolineatus (Decapoda: Brachyura: Grapsidae) and phylogenetic associations in Brachyura

Complete mitochondrial genomes (mitogenomes) can provide useful information for phylogenetic relationships, gene rearrangement, and evolutionary traits. In this study, we determined the complete mitochondrial DNA sequence of the herbivorous crab Grapsus albolineatus. It is a typical metazoan mitochondrial genome. The total size is 15,583 bp, contains the entire set of 37 genes, and has an AT-rich region. Then, 23 of the 37 genes were encoded by the heavy (+) strand while 14 are encoded by the light (−) strand. Compared with the pan-crustacean ground pattern, two tRNA genes (tRNA-His and tRNA-Gln) were rearranged and the tandem duplication/random loss model was used to explain the observed gene rearrangements. The phylogenetic results showed that all Grapsidae crabs clustered together as a group. Furthermore, the monophyly of each family was well supported, with the exception of Menippidae. In general, the results obtained in this study will contribute to the better understanding of gene rearrangements in Grapsidae crab mitogenomes and provide new insights into the phylogeny of Brachyura.

Morphological and histological description of the midgut caeca in true crabs (Malacostraca: Decapoda: Brachyura): origin, development and potential role

Background

The decapods are a major group of crustaceans that includes shrimps, prawns, crayfishes, lobsters, and crabs. Several studies focused on the study of the digestive system of the decapods, constituted by the oesophagus, stomach, midgut tract, midgut gland, and hindgut. Nevertheless, in the midgut tract there are associated a set of organs called “midgut caeca”, which are among the most controversial and less studied digestive organs of this group. This work used the common spider crab Maja brachydactyla Balss, 1922 as a model to resolve the origin, development, and potential role of the midgut caeca. Such organs were studied in the larvae (zoea I, zoea II, megalopa), first juveniles, and adult phases, being employed traditional and modern techniques: dissection, micro-computed tomography (Micro-CT), and light and electron microscopical analyses (TEM and SEM).

Results

The common spider crab has a pair of anterior midgut caeca and a single posterior caecum that originate from the endoderm germ layer: they develop from the midgut tract, and their epithelium is composed by secretory cells while lacking a cuticle lining. The midgut caeca are small buds in the newly hatched larvae, enlarge linearly during the larval development, and then continue growing until became elongated and coiled blind-tubules in adults. The adult midgut caeca are internally folded to increase their inner surface. The electron microscopy observations showed that the midgut caeca are highly active organs with important macroapocrine and microapocrine secretory activity. Our results suggest that the role of the caeca might be related to the digestive enzyme secretion. The secretory activity should increase as the animal grows in size.

Conclusion

The present study resolves the embryonic origin of the midgut caeca (endoderm derived organs), development (general lengthening starting from small buds), and role (active secretory organs). The secretory activity of the midgut caeca should be incorporated in the current models of the digestive physiology in different decapod taxa.

Comparative Characterization of Mitogenomes From Five Orders of Cestodes (Eucestoda: Tapeworms)

The recognized potential of using mitogenomics in phylogenetics and the more accessible use of high-throughput sequencing (HTS) offer an opportunity to investigate groups of neglected organisms. Here, we leveraged HTS to execute the most comprehensive documentation of mitogenomes for cestodes based on the number of terminals sequenced. We adopted modern approaches to obtain the complete mitogenome sequences of 86 specimens representing five orders of cestodes (three reported for the first time: Phyllobothriidea, “Tetraphyllidea” and Trypanorhyncha). These complete mitogenomes represent an increase of 41% of the mitogenomes available for cestodes (61–147) and an addition of 33% in the representativeness of the cestode orders. The complete mitochondrial genomes are conserved, circular, encoded in the same strand, and transcribed in the same direction, following the pattern observed previously for tapeworms. Their length varies from 13,369 to 13,795 bp, containing 36 genes in total. Except for the Trypanorhyncha specimen, the gene order of the other four cestode orders sequenced here suggests that it could be a synapomorphy for the acetabulate group (with a reversion for taenids). Our results also suggest that no single gene can tell all the evolutionary history contained in the mitogenome. Therefore, cestodes phylogenies based on a single mitochondrial marker may fail to capture their evolutionary history. We predict that such phylogenies would be improved if conducted under a total evidence framework. The characterization of the new mitochondrial genomes is the first step to provide a valuable resource for future studies on the evolutionary relationships of these groups of parasites.

Mitochondrial Genome of Nonmodel Marine Metazoans by Next-Generation Sequencing (NGS)

Mitochondrial genomes (mtgenome) represent an important source of information for addressing fundamental evolutionary, phylogeographic, systematic, and ecological questions in marine organisms. In the last two decades the advent of high-throughput next-generation sequencing (NGS) has provided an unprecedented possibility to access large amount of genomic data and, as such, there has been a rapid growth in mtgenome resources and studies. In particular, NGS strategies represent a great advantage for investigating nonmodel marine organisms for which no or limited genomic resources are available. Here, we describe a routinely used standardized protocol to obtain mtgenome of nonmodel marine organisms by NGS. The protocol is composed of five main steps, including DNA extraction, DNA fragmentation, library preparation, high-throughput sequencing, and bioinformatic analyses. Each of the first three steps is followed by size/quality and concentration validations. The advantages of the described protocol rely on the assumption that no a priori information on mtgenome of the studied organism is needed and on its versatility as researchers may choose several kits for DNA extraction and library preparation and adopt different methods for DNA fragmentation depending on their needs, experience, and suppliers.

High-Throughput Sequencing of faeces provides evidence for dispersal of parasites and pathogens by migratory waterbirds

Examination of faecal material has demonstrated how a broad range of organisms are distributed by bird movements. Such research has largely focused on dispersal of plant seeds by frugivores and of freshwater organisms by waterbirds. However, with few exceptions (e.g. avian influenza, Ebola virus), there is a dearth of evidence for transport of parasites and pathogens. High-throughput sequencing methods now provide a powerful means of addressing this knowledge gap by elucidating faecal contents in unprecedented detail. We collected faeces excreted by a range of migratory waterbirds in south-west Spain and pooled faecal DNA to create libraries reflective of feeding behavior. We created sets of libraries using high-throughput metagenomic and amplicon sequencing. For the latter we employed two sets of primers to broadly target the V4 region of the 18S rRNA gene (one set amplifying the region across all eukaryotes, the other excluding amplification of metazoans). Libraries revealed a wide diversity of eukaryotes, including parasites of the faecal producers themselves, parasites of food items, or those incidentally ingested. We also detected novel microbial eukaryotic taxa and found that parasite assemblage profiles were relatively distinct. Comparing the performance of the methods used supports their joint use for future studies of diversity and abundance. Because viable stages of many parasites are likely to be present in faeces, our results suggest significant levels of bird-mediated dispersal of parasites (both from avian and other hosts). Our methods revealed much hidden biodiversity, and allowed identification of the individuals who produced the faecal samples to species level, facilitating the study of interaction networks.

The mitochondrial genome of Faughnia haani (Stomatopoda): novel organization of the control region and phylogenetic position of the superfamily Parasquilloidea

BACKGROUND

Stomatopod crustaceans are aggressive marine predators featuring complex compound eyes and powerful raptorial appendages used for "smashing" or "spearing" prey and/or competitors. Among them, parasquilloids (superfamily Parasquilloidea) possess eyes with 2-3 midband rows of hexagonal ommatidia and spearing appendages. Here, we assembled and analyzed the complete mitochondrial genome of the parasquilloid Faughnia haani and explored family- and superfamily-level phylogenetic relationships within the Stomatopoda based on mitochondrial protein coding genes (PCGs).

RESULTS

The mitochondrial genome of F. haani is 16,089 bp in length and encodes 13 protein coding genes (PCGs), 22 transfer RNA genes, 2 ribosomal RNA genes, and a control region that is relatively well organized, containing 2 GA-blocks, 4 poly-T stretches, various [TA(A)]n-blocks, and 2 hairpin structures. This organized control region is likely a synapomorphic characteristic in the Stomatopoda. Comparison of the control region among superfamilies shows that parasquilloid species are more similar to gonodactyloids than to squilloids and lysiosquilloids given the presence of various  poly-T stretches between the hairpin structures and [TA(A)]n-blocks. Synteny is identical to that reported for other stomatopods and corresponds to the Pancrustacea ground pattern. A maximum-likelihood phylogenetic tree based on PCGs revealed that Parasquilloidea is sister to Lysiosquilloidea and Gonodactyloidea and not to Squilloidea, contradicting previous phylogenetic studies.

CONCLUSIONS

The novel phylogenetic position of Parasquilloidea revealed by our study indicates that 'spearing' raptorial appendages are plesiomorphic and that the 'smashing' type is either derived (as reported in previous studies) or apomorphic. Our results raise the possibility that the spearing raptorial claw may have independently evolved twice. The superfamily Parasquilloidea exhibits a closer relationship with other stomatopod superfamilies with a different raptorial claw type and with dissimilar numbers of midband rows of hexagonal ommatidia. Additional studies focusing on the assembly of mitochondrial genomes from species belonging to different genera, families, and superfamilies within the order Stomatopoda are warranted to reach a robust conclusion regarding the evolutionary history of this iconic clade based on mitochondrial PCGs.

Genome sequencing guide: An introductory toolbox to whole-genome analysis methods

To fully appreciate genetics, one must understand the link between genotype (DNA sequence) and phenotype (observable characteristics). Advances in high-throughput genomic sequencing technologies and applications, so-called "-omics," have made genetic sequencing readily available across fields in biology from applications in non-traditional study organisms to precision medicine. Thus, understanding these tools is critical for any biologist, especially those early in their career. This comprehensive review discusses the chronological development of different sequencing methods, the bioinformatics steps to analyzing this data, and social and ethical issues raised by these techniques that must be discussed and evaluated, including anticipatory guides and discussion questions for active engagement in the classroom. Additionally, the Supporting Information includes a case study to apply technical and ethical concepts from the text.

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.

The first mitochondrial genome of the genus Exhippolysmata (Decapoda: Caridea: Lysmatidae), with gene rearrangements and phylogenetic associations in Caridea

The complete mitochondrial genome (mitogenome) of animals can provide useful information for evolutionary and phylogenetic analyses. The mitogenome of the genus Exhippolysmata (i.e., Exhippolysmata ensirostris) was sequenced and annotated for the first time, its phylogenetic relationship with selected members from the infraorder Caridea was investigated. The 16,350 bp mitogenome contains the entire set of 37 common genes. The mitogenome composition was highly A + T biased at 64.43% with positive AT skew (0.009) and negative GC skew (- 0.199). All tRNA genes in the E. ensirostris mitogenome had a typical cloverleaf secondary structure, except for trnS1 (AGN), which appeared to lack the dihydrouridine arm. The gene order in the E. ensirostris mitogenome was rearranged compared with those of ancestral decapod taxa, the gene order of trnL2-cox2 changed to cox2-trnL2. The tandem duplication-random loss model is the most likely mechanism for the observed gene rearrangement of E. ensirostris. The ML and BI phylogenetic analyses place all Caridea species into one group with strong bootstrap support. The family Lysmatidae is most closely related to Alpheidae and Palaemonidae. These results will help to better understand the gene rearrangements and evolutionary position of E. ensirostris and lay a foundation for further phylogenetic studies of Caridea.

The first two complete mitogenomes of the order Apodida from deep-sea chemoautotrophic environments: New insights into the gene rearrangement, origin and evolution of the deep-sea sea cucumbers

The deep-sea ecosystem is considered as the largest and most remote biome of the world. It is meaningful and important to elucidate the life origins by exploring the origin and adaptive genetic mechanisms of the large deep-sea organisms. Sea cucumbers (Holothuroidea) are abundant and economically important group of echinoderms, living from the shallow-waters to deep-sea. In this study, we present the mitochondrial genomes of the sea cucumber Chiridota heheva and Chiridota sp. collected from the deep-sea cold seep and hydrothermal vent, respectively. This is the first reported mitochondrial genomes from the order Apodida. The mitochondrial genomes of C. heheva (17,200 bp) and Chiridota sp. (17,199 bp) display novel gene arrangements with the first protein-coding gene rearrangements in the class Holothuroidea. Bases composition analysis showed that the A + T content of deep-sea holothurians were significantly higher than that of the shallow-water groups. We compared the arrangement of genes from the 24 available holothurian mitogenomes and found that the transposition, reverse transposition and tandem-duplication-random-losses (TDRL) may be involved in the evolution of mitochondrial gene arrangements in Holothuroidea. Phylogenetic analysis revealed that the Apodida clustered with Elasipodida, forming two basal deep-sea holothurian clades. The divergence between the deep-sea and shallow-water holothurians was located at 386.93 Mya, during the Late Devonian. Mitochondrial protein-coding genes of deep-sea holothurians underwent relaxed purifying selection. There are 57 positive selected amino acids sites for some mitochondrial genes of the three deep-sea clades, implying they may involve in the adaption of deep-sea sea cucumbers.

Comparative mitochondrial genomics of sponge-dwelling snapping shrimps in the genus Synalpheus: Exploring differences between eusocial and non-eusocial species and insights into phylogenetic relationships in caridean shrimps

The genus Synalpheus is a cosmopolitan clade of marine shrimps found in most tropical regions. Species in this genus exhibit a range of social organizations, including pair-forming, communal breeding, and eusociality, the latter only known to have evolved within this genus in the marine realm. This study examines the complete mitochondrial genomes of seven species of Synalpheus and explores differences between eusocial and non-eusocial species considering that eusociality has been shown before to affect the strength of purifying selection in mitochondrial protein coding genes. The AT-rich mitochondrial genomes of Synalpheus range from 15,421 bp to 15,782 bp in length and comprise, invariably, 13 protein-coding genes (PCGs), two ribosomal RNA genes, and 22 transfer RNA genes. A 648 bp to 994 bp long intergenic space is assumed to be the D-loop. Mitochondrial gene synteny is identical among the studied shrimps. No major differences occur between eusocial and non-eusocial species in nucleotide composition and codon usage profiles of PCGs and in the secondary structure of tRNA genes. Maximum likelihood phylogenetic analysis of the complete concatenated PCG complement of 90 species supports the monophyly of the genus Synalpheus and its family Alpheidae. Moreover, the monophyletic status of the caridean families Alvinocaridae, Atyidae, Thoridae, Lysmatidae, Palaemonidae, and Pandalidae within caridean shrimps are fully or highly supported by the analysis. We therefore conclude that mitochondrial genomes contain sufficient phylogenetic information to resolve relationships at high taxonomic levels within the Caridea. Our analysis of mitochondrial genomes in the genus Synalpheus contributes to the understanding of the coevolution between genomic architecture and sociality in caridean shrimps and other marine organisms.

How to become a crab: Phenotypic constraints on a recurring body plan

A fundamental question in biology is whether phenotypes can be predicted by ecological or genomic rules. At least five cases of convergent evolution of the crab-like body plan (with a wide and flattened shape, and a bent abdomen) are known in decapod crustaceans, and have, for over 140 years, been known as "carcinization." The repeated loss of this body plan has been identified as "decarcinization." In reviewing the field, we offer phylogenetic strategies to include poorly known groups, and direct evidence from fossils, that will resolve the history of crab evolution and the degree of phenotypic variation within crabs. Proposed ecological advantages of the crab body are summarized into a hypothesis of phenotypic integration suggesting correlated evolution of the carapace shape and abdomen. Our premise provides fertile ground for future studies of the genomic and developmental basis, and the predictability, of the crab-like body form.

Phylogenetic relations and mitogenome-wide similarity metrics reveal monophyly of Penaeus sensu lato

Splitting of the genus Penaeus sensu lato into six new genera based on morphological features alone has been controversial in penaeid shrimp taxonomy. Several studies focused on building phylogenetic relations among the genera of Penaeus sensu lato. However, they lack in utilizing full mitochondrial DNA genome of shrimp representing all the six controversial genera. For the first time, the present study targeted the testing of all the six genera of Penaeus sensu lato for phylogenetic relations utilizing complete mitochondrial genome sequence. In addition, the study reports for the first time about the complete mitochondrial DNA genome sequence of Fenneropenaeus indicus, an important candidate species in aquaculture and fisheries, and utilized it for phylogenomics. The maximum likelihood and Bayesian approaches were deployed to generate and comprehend the phylogenetic relationship among the shrimp in the suborder, Dendrobranchiata. The phylogenetic relations established with limited taxon sampling considered in the study pointed to the monophyly of Penaeus sensu lato and suggested collapsing of the new genera to a single genus. Further, trends in mitogenome-wide estimates of average amino acid identity in the order Decapoda and the genus Penaeus sensu lato supported restoration of the old genus, Penaeus, rather promoting the creation of new genera.

The Pontastacus leptodactylus (Astacidae) Repeatome Provides Insight Into Genome Evolution and Reveals Remarkable Diversity of Satellite DNA

Pontastacus leptodactylus is a native European crayfish species found in both freshwater and brackish environments. It has commercial importance for fisheries and aquaculture industries. Up till now, most studies concerning P. leptodactylus have focused onto gaining knowledge about its phylogeny and population genetics. However, little is known about the chromosomal evolution and genome organization of this species. Therefore, we performed clustering analysis of a low coverage genomic dataset to identify and characterize repetitive DNA in the P. leptodactylus genome. In addition, the karyogram of P. leptodactylus (2n = 180) is presented here for the first time consisting of 75 metacentric, 14 submetacentric, and a submetacentric/metacentric heteromorphic chromosome pair. We determined the genome size to be at ~18.7 gigabase pairs. Repetitive DNA represents about 54.85% of the genome. Satellite DNA repeats are the most abundant type of repetitive DNA, making up to ~28% of the total amount of repetitive elements, followed by the Ty3/Gypsy retroelements (~15%). Our study established a surprisingly high diversity of satellite repeats in P. leptodactylus. The genome of P. leptodactylus is by far the most satellite-rich genome discovered to date with 258 satellite families described. Of the five mapped satellite DNA families on chromosomes, PlSAT3-411 co-localizes with the AT-rich DAPI positive probable (peri)centromeric heterochromatin on all chromosomes, while PlSAT14-79 co-localizes with the AT-rich DAPI positive (peri)centromeric heterochromatin on one chromosome and is also located subterminally and intercalary on some chromosomes. PlSAT1-21 is located intercalary in the vicinity of the (peri)centromeric heterochromatin on some chromosomes, while PlSAT6-70 and PlSAT7-134 are located intercalary on some P. leptodactylus chromosomes. The FISH results reveal amplification of interstitial telomeric repeats (ITRs) in P. leptodactylus. The prevalence of repetitive elements, especially the satellite DNA repeats, may have provided a driving force for the evolution of the P. leptodactylus genome.

Reanalysis and Revision of the Complete Mitochondrial Genome of Artemia urmiana Günther, 1899 (Crustacea: Anostraca)

In the previously published mitochondrial genome sequence of Artemia urmiana (NC_021382 [JQ975176]), the taxonomic status of the examined Artemia had not been determined, due to parthenogenetic populations coexisting with A. urmiana in Urmia Lake. Additionally, NC_021382 [JQ975176] has been obtained with pooled cysts of Artemia (0.25 g cysts consists of 20,000–25,000 cysts), not a single specimen. With regard to coexisting populations in Urmia Lake, and intra- and inter-specific variations in the pooled samples, NC_021382 [JQ975176] cannot be recommended as a valid sequence and any attempt to attribute it to A. urmiana or a parthenogenetic population is unreasonable. With the aid of next-generation sequencing methods, we characterized and assembled a complete mitochondrial genome of A. urmiana with defined taxonomic status. Our results reveal that in the previously published mitogenome (NC_021382 [JQ975176]), tRNA-Phe has been erroneously attributed to the heavy strand but it is encoded in the light strand. There was a major problem in the position of the ND5. It was extended over the tRNA-Phe, which is biologically incorrect. We have also identified a partial nucleotide sequence of 311 bp that was probably erroneously duplicated in the assembly of the control region of NC_021382 [JQ975176], which enlarges the control region length by 16%. This partial sequence could not be recognized in our assembled mitogenome as well as in 48 further examined specimens of A. urmiana. Although, only COX1 and 16S genes have been widely used for phylogenetic studies in Artemia, our findings reveal substantial differences in the nucleotide composition of some other genes (including ATP8, ATP6, ND3, ND6, ND1 and COX3) among Artemia species. It is suggested that these markers should be included in future phylogenetic studies.

Eocarcinus praecursor Withers, 1932 (Malacostraca, Decapoda, Meiura) is a stem group brachyuran

Beginning with the description by Withers in 1932, Eocarcinus praecursor from the Jurassic has long been considered the oldest representative of the Brachyura. In 2010 Feldmann and Schweitzer re-investigated the specimens of E. praecursor and expressed doubts about the brachyuran nature of this species. Among other characters, the suspected existence of small chelae in the 2nd or 3rd pereopods led them to the conclusion that E. praecursor must be removed from the Brachyura and rather be seen as a representative of the Anomala. However, Anomala also do not possess chelae on the 2nd and 3rd pereopods. This contradiction and other aspects initiated a new investigation of E. praecursor. It can be shown that neither the 2nd nor the 3rd pereopods of E. praecursor are chelate. Furthermore, there are no other derived characters shared with anomalans. By contrast, there are a number of apomorphies shared with Brachyura such as the shape and articulation of the large chelae and the attachment points of the last two pereopods. However, not all apomorphies of the crown group are present yet. Thus, E. praecursor is a stem group representative, which allows statements about individual steps in the evolution of the set of characters of the crown group Brachyura.

The complete mitochondrial genome of the eusocial sponge-dwelling snapping shrimp Synalpheus microneptunus

In the marine realm, eusociality is only known to have evolved within a clade of sponge-dwelling snapping shrimps in the genus Synalpheus. Deciphering the genomic underpinnings of eusociality in these marine shrimps has been limited by the sparse genomic resources in this genus. Here, we report, for a eusocial shrimp Synalpheus microneptunus, a complete mitochondrial genome (22X coverage) assembled from short Illumina 150 bp paired-end reads. The 15,603 bp long mitochondrial genome of S. microneptunus is AT-rich and includes 13 protein-coding genes (PCGs), 2 ribosomal RNA genes, 22 transfer RNA genes and an 834 bp intergenic region assumed to be the D-loop. The gene order is identical to that reported for most caridean shrimps and corresponds to the presumed Pancrustacean ground pattern. All PCGs showed signs of purifying selection, with K_A/K_S <<1 across the whole PCGs and most sliding windows within PCGs. Maximum-likelihood and Bayesian inference phylogenetic analyses of 13 PCGs and 68 terminals supports the monophyly of the Caridea and the family Alpheidae. The complete mitochondrial genome of the eusocial shrimp Synalpheus microneptunus will contribute to a better understanding of the selective pressures and rates of molecular evolution in marine eusocial animals.