Complete mitochondrial genome of the Chinese spiny lobster Panulirus stimpsoni (Crustacea: Decapoda): genome characterization and phylogenetic considerations

Mitogenomic recognition of incognito lineages in the mud spiny lobster Panulirus polyphagus (Herbst, 1793): A tale of unique genetic structuring and diversification

This study provides the first documentation of three deep conspecific lineages within Panulirus polyphagus in the Indian Ocean, bridging the gap in genetic research. Comparative mitogenomics between lineages (L) at both species and family levels, evolutionary relationships and heterogeneity of sequence divergence within Decapoda, and divergence time estimation were performed. The characterized mitogenomes ranged from 15,685-15,705 bp in size and exhibited a typical pancrustacean pattern. Among the three lineages, L1 predominated the Bay of Bengal, L2 the Arabian Sea, and L2.a, a less common lineage genetically closer to L2, was restricted to the latter region. A minor lineage L1.a, was observed in the Coral Triangle area. All PCGs displayed evidence of purifying selection across species and family levels. The largest genetic distance (K2P) between lineages was 9 %, notably between L1.a and L2.a. The phylogenetic tree subdivided the Achelates into Palinuridae and Scyllaridae, and the topology demonstrated a distinct pattern of lineage diversification within P. polyphagus. AliGROOVE analysis revealed no discernible divergence in Decapoda. The diversification of P. polyphagus appears to have occurred during Miocene, with further diversification in Pliocene. Furthermore, genetic stocks and population connectivity recognized here will provide valuable insight for spatial management planning of this dwindling resource.

DNA barcoding and morphological identification of spiny lobsters in South Korean waters: a new record of Panulirus longipes and Panulirus homarus homarus

To date, 19 species of spiny lobsters from the genus Panulirus have been discovered, of which only P. japonicus, P. penicilatus, P. stimpsoni, and P. versicolor have been documented in South Korean waters. In this study, we aimed to identify and update the current list of spiny lobster species that inhabit South Korean waters based on the morphological features and the phylogenetic profile of cytochrome oxidase I (COI) of mitochondrial DNA (mtDNA). Spiny lobsters were collected from the southern and eastern coasts of Jeju Island, South Korea. Phylogenetic analyses were performed using neighbor-joining (NJ), maximum likelihood (ML), and Bayesian inference (BI) methods. The ML tree was used to determine the spiny lobster lineages, thereby clustering the 17 specimens collected in this study into clades A, B, C, and D, which were reciprocally monophyletic with P. japonicus, P. homarus homarus, P. longipes, and P. stimpsoni, respectively. These clades were also supported by morphological examinations. Interestingly, morphological variations, including the connected pleural and transverse groove at the third abdominal somite, were observed in four specimens that were genetically confirmed as P. japonicus. This finding is novel within the P. japonicus taxonomical reports. Additionally, this study updates the documentation of spiny lobsters inhabiting South Korean waters as P. longipes and P. homarus homarus were recorded for the first time in this region.

The complete mitochondrial genome of the red-jointed brackish-water fiddler crab Minuca minax (LeConte 1855) (Brachyura: Ocypodidae): New family gene order, and purifying selection and phylogenetic informativeness of protein coding genes

Minuca minax is a semi-terrestrial crustacean that commonly lives in low salinity, riverine habitats along the shores of the eastern United States. This study reports, for the first time, the complete mitochondrial genome of M. minax. The AT-rich mitochondrial genome of M. minax is 15,937 bp in length and comprised of 13 protein-coding genes (PCGs), 2 ribosomal RNA genes, and 22 transfer RNA genes. A single 737 bp long intergenic space is assumed to be the D-loop. Most of the PCGs and tRNA genes are encoded in the L-strand. The gene order observed in the mitochondrial genome of M. minax is new although almost identical to that reported in confamiliar species. In all other confamiliar species to which M. minax is compared, the positions of the trnQ gene and the trnI gene are switched. K_A/K_S ratios calculated for all mitochondrial PCGs show values of <1, indicating that these PCGs are evolving under purifying selection. A maximum likelihood phylogenetic analysis (concatenated PCGs [n = 13], 15 species) supports the monophyly of the subfamilies Ocypodinae and Gelaminidae. Mitochondrial PCGs have enough phylogenetic information to reveal relationships supporting higher taxonomic levels within this family. The knowledge of a complete mitochondrial genome from the red-jointed brackish-water fiddler crab M. minax contributes to the better understanding of meta-population connectivity and the mechanisms involved in the adaptation of marine organisms to near-limnic conditions.

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.

Mitochondrial genome of Chinese grass shrimp, Palaemonetes sinensis and comparison with other Palaemoninae species

The mitogenome of Chinese grass shrimp, Palaemonetes sinensis, was determined through Illumina sequencing, and the basic characteristics and gene arrangement were analyzed. The mitogenome of P. sinensis was 15955 bp in length, consisting of 13 protein-coding genes (PCGs), 22 tRNA genes, 2 rRNA genes and one control region, with tightly packed. 33 of these genes were encoded on the heavy strand, and the remainders encoded on the light strand. The composition of P. sinensis mitogenome presented a strong A + T bias, which account for 66.7%. All PCGs were initiated by a canonical ATN codon, except nad5, which was initiated by GTG. The termination codons of the PCGs were TAA, TAG and T-. The secondary structures of 22 tRNAs of P. sinensis had the typical clover structure, except of trnS1 owing to the lack of dihydroxyuridine (DHU) arm. Gene order comparison of P. sinensis and previously-sequenced Palaemoninae revealed a unique translocation between trnT and trnP in Macrobrachium. The phylogenetic analyses showed that three Exopalaemon species formed a monophyletic group and then clustered with two Palaemon species and P. sinensis successively whereas Macrobrachium clustered with Palaemon capensis in the other clade.

Characterization of the mitochondrial genome of an ancient amphipod Halice sp. MT-2017 (Pardaliscidae) from 10,908 m in the Mariana Trench

Small amphipods (Halice sp. MT-2017) with body length <1 cm were collected from the Challenger Deep (~10,920 m below sea level). The divergence time of their lineage was approximately 109 Mya, making this group ancient compared to others under study. The mitochondrial genome of Halice sp. shared the usual gene components of metazoans, comprising 13 protein coding genes (PCGs), 22 transfer RNAs (tRNAs), and 2 ribosomal RNAs (rRNAs). The arrangement of these genes, however, differed greatly from that of other amphipods. Of the 15 genes that were rearranged with respect to the pancrustacean gene pattern, 12 genes (2 PCGs, 2 rRNAs, and 8 tRNAs) were both translocated and strand-reversed. In contrast, the mitochondrial genomes in other amphipods never show so many reordered genes, and in most instances, only tRNAs were involved in strand-reversion-coupled translocation. Other characteristics, including reversed strand nucleotide composition bias, relatively higher composition of non-polar amino acids, and lower evolutionary rate, were also identified. Interestingly, the latter two features were shared with another hadal amphipod, Hirondellea gigas, suggesting their possible associations with the adaptation to deep-sea extreme habitats. Overall, our data provided a useful resource for future studies on the evolutionary and adaptive mechanisms of hadal faunas.

The complete mitochondrial genome of the Caribbean spiny lobster Panulirus argus

Panulirus argus is a keystone species and target of the most lucrative fishery in the Caribbean region. This study reports, for the first time, the complete mitochondrial genome of Panulirus argus (average coverage depth nucleotide-1 = 70×) assembled from short Illumina 150 bp PE reads. The AT-rich mitochondrial genome of Panulirus argus was 15 739 bp in length and comprised 13 protein-coding genes (PCGs), 2 ribosomal RNA genes, and 22 transfer RNA genes. A single 801 bp long intergenic space was assumed to be the D-loop. Most of the PCGs were encoded on the H-strand. The gene order observed in the mitochondrial genome of Panulirus argus corresponds to the presumed Pancrustacean ground pattern. KA/KS ratios calculated for all mitochondrial PCGs showed values < 1, indicating that all these PCGs are evolving under purifying selection. A maximum likelihood phylogenetic analysis (concatenated PCGs [n = 13], 154 arthropods) supported the monophyly of the Achelata and other infraorders within the Decapoda. Mitochondrial PCGs have enough phylogenetic informativeness to explore high-level genealogical relationships in the Pancrustacea. The complete mitochondrial genome of the Caribbean spiny lobster Panulirus argus will contribute to the better understanding of meta-population connectivity in this keystone overexploited species.

Analysis of the complete mitochondrial genome sequence of Palinura homarus

The complete mitochondrial genome sequence of Palinura homarus was obtained using PCR amplification and walking sequencing (GenBank accession no. JN_542716). The complete mitochondrial genome of P. homarus was 15,665 bp long and showed significant AT bias (67% AT content, 33% GC content). The A + T-rich region included copy-related control information and a poly (dT) structure that related to replication and transcription. In this study, the gene arrangement was consistent with other Palinura mitochondrial genomes and the sequence was strikingly similar to Panulirus ornatus, which would be useful in species identification and natural resources conservation.

Evolution of mitochondrial genomes in Baikalian amphipods

BACKGROUND

Amphipods (Crustacea) of Lake Baikal are a very numerous and diverse group of invertebrates generally believed to have originated by adaptive radiation. The evolutionary history and phylogenetic relationships in Baikalian amphipods still remain poorly understood. Sequencing of mitochondrial genomes is a relatively feasible way for obtaining a set of gene sequences suitable for robust phylogenetic inferences. The architecture of mitochondrial genomes also may provide additional information on the mechanisms of evolution of amphipods in Lake Baikal.

RESULTS

Three complete and four nearly complete mitochondrial genomes of Baikalian amphipods were obtained by high-throughput sequencing using the Illumina platform. A phylogenetic inference based on the nucleotide sequences of all mitochondrial protein coding genes revealed the Baikalian species to be a monophyletic group relative to the nearest non-Baikalian species with a completely sequenced mitochondrial genome - Gammarus duebeni. The phylogeny of Baikalian amphipods also suggests that the shallow-water species Eulimnogammarus has likely evolved from a deep-water ancestor, however many other species have to be added to the analysis to test this hypothesis. The gene order in all mitochondrial genomes of studied Baikalian amphipods differs from the pancrustacean ground pattern. Mitochondrial genomes of four species possess 23 tRNA genes, and in three genomes the extra tRNA gene copies have likely undergone remolding. Widely varying lengths of putative control regions and other intergenic spacers are typical for the mitochondrial genomes of Baikalian amphipods.

CONCLUSIONS

The mitochondrial genomes of Baikalian amphipods display varying organization suggesting an intense rearrangement process during their evolution. Comparison of complete mitochondrial genomes is a potent approach for studying the amphipod evolution in Lake Baikal.

Complete mitochondrial genome of Australian spiny lobster, Panulirus cygnus (George, 1962) (Crustacea: Decapoda: Palinuridae) from coast of Australia

We characterized the total mitochondrial genome of Australian spiny lobster, Panulirus cygnus (George, 1962), which is found along the western coast of Australia. Total mitochondrial genome length of P. cygnus was 15 724 bp, in which 13 proteins, 2 ribosomal RNAs, 22 transfer RNAs, and a putative control region were encoded. Nine and four protein-coding genes are encoded on the H-strand and on the L-strand, respectively. According to the phylogenetic analysis, P. cygnus was most closely related to Panulirus japonicus among the compared six species belonging to Palinuridae. Although overall gene organization was the same, the putative control region (between SrRNA gene and tRNA^Ilel) is least similar to one another among mitochondrial genomes from the compared six species belonging to Palinuridae.

The complete mitochondrial genome of the Japanese fan lobster Ibacus ciliatus (Crustacea, Achelata, Scyllaridae)

Slipper lobsters are a unique group of decapod crustaceans; they have cylindrical or flattened bodies and belong to the family Scyllaridae. The genus Ibacus (Leach, 1815) (Achelata, Scyllaridae, Ibacinae) consists of eight recognized species to date, all of which occur in the Indo-West Pacific Ocean and are commercially important seafood species. The complete mitochondrial genome of the Japanese fan lobster Ibacus ciliatus (Von Siebold, 1824) is 15,696 bp in size and consists of 13 protein-coding genes, 22 transfer RNAs, 2 ribosomal RNAs, and a control region of 783 bp. The base composition of I. ciliatus is 35.8% A, 34.7% T, 19.0% C, and 10.5% G, with an overall AT content of 70.5%. The mitogenome of I. ciliatus was found to have gene arrangement and transcriptional polarity identical to that of the American lobster Homarus americanus, showing the pancrustacean ground pattern. Here, we present the complete mitogenome sequence of I. ciliatus; it is the first mitogenome information from the subfamily Ibacinae, and represents the second for the family Scyllaridae.

The complete mitogenome of Cherax monticola (Crustacea: Decapoda: Parastacidae), a large highland crayfish from New Guinea

The complete mitochondrial genome of a highland freshwater crayfish, Cherax monticola, was recovered by shotgun sequencing. The mitogenome consists of 15,917 base pairs containing 13 protein-coding genes, 2 ribosomal subunit genes, 22 transfer RNAs and a non-coding AT-rich region. The base composition of C. monticola is 33.46% for T, 21.48% for C, 33.71% for A and 11.35% for G, with an AT bias of 67.17%.

Integrated shotgun sequencing and bioinformatics pipeline allows ultra-fast mitogenome recovery and confirms substantial gene rearrangements in Australian freshwater crayfishes

BACKGROUND

Although it is possible to recover the complete mitogenome directly from shotgun sequencing data, currently reported methods and pipelines are still relatively time consuming and costly. Using a sample of the Australian freshwater crayfish Engaeus lengana, we demonstrate that it is possible to achieve three-day turnaround time (four hours hands-on time) from tissue sample to NCBI-ready submission file through the integration of MiSeq sequencing platform, Nextera sample preparation protocol, MITObim assembly algorithm and MITOS annotation pipeline.

RESULTS

The complete mitochondrial genome of the parastacid freshwater crayfish, Engaeus lengana, was recovered by modest shotgun sequencing (1.2 giga bases) using the Illumina MiSeq benchtop sequencing platform. Genome assembly using the MITObim mitogenome assembler recovered the mitochondrial genome as a single contig with a 97-fold mean coverage (min. = 17; max. = 138). The mitogenome consists of 15,934 base pairs and contains the typical 37 mitochondrial genes and a non-coding AT-rich region. The genome arrangement is similar to the only other published parastacid mitogenome from the Australian genus Cherax.

CONCLUSIONS

We infer that the gene order arrangement found in Cherax destructor is common to Australian crayfish and may be a derived feature of the southern hemisphere family Parastacidae. Further, we report to our knowledge, the simplest and fastest protocol for the recovery and assembly of complete mitochondrial genomes using the MiSeq benchtop sequencer.

Molecular characterization and evolution of haemocyanin from the two freshwater shrimps Caridina multidentata (Stimpson, 1860) and Atyopsis moluccensis (De Haan, 1849)

Haemocyanin (Hc) is a copper-containing respiratory protein, floating freely dissolved in the hemolymph of many arthropod species. A typical haemocyanin is a hexamer or oligohexamer of six identical or similar subunits, with a molecular mass around 75 kDa each. In the crustaceans, the haemocyanins appear to be restricted to the remipedes and the malacostracans. We have investigated the haemocyanins of two freshwater shrimps, the Amano shrimp Caridina multidentata and the bamboo shrimp Atyopsis moluccensis. We obtained three full-length and one partial cDNA sequences of haemocyanin subunits from the Amano shrimp, which were assigned to the α- and γ-types of decapod haemocyanin subunits. Three complete and two partial haemocyanin cDNA sequences were obtained from the bamboo shrimp, which represent subunit types α, β and γ. This is the first time that sequences of all three subunit types of the decapod haemocyanins were obtained from a single species. However, mass spectrometry analyses identified only α- and γ-type subunits, suggesting that a β-subunit is not a major component of the native haemocyanin of the bamboo shrimp. Phylogenetic and molecular clock analyses showed that malacostracan haemocyanins commenced to diversify into distinct subunit types already ~515 million years ago. β-subunits diverged first, followed by α- and γ-type subunits ~396 million years ago. The haemocyanins of phyllocarids and peracarids form distinct clades within the α/γ-cluster. Within the Caridea, an early divergence of distinct α-type subunits occurred ~200 MYA. The tree of the γ-subunits suggests a common clade of the Caridea (shrimps) and Penaeidae (prawns).

Evolution and phylogeny of the mud shrimps (Crustacea: Decapoda) revealed from complete mitochondrial genomes

BACKGROUND

The evolutionary history and relationships of the mud shrimps (Crustacea: Decapoda: Gebiidea and Axiidea) are contentious, with previous attempts revealing mixed results. The mud shrimps were once classified in the infraorder Thalassinidea. Recent molecular phylogenetic analyses, however, suggest separation of the group into two individual infraorders, Gebiidea and Axiidea. Mitochondrial (mt) genome sequence and structure can be especially powerful in resolving higher systematic relationships that may offer new insights into the phylogeny of the mud shrimps and the other decapod infraorders, and test the hypothesis of dividing the mud shrimps into two infraorders.

RESULTS

We present the complete mitochondrial genome sequences of five mud shrimps, Austinogebia edulis, Upogebia major, Thalassina kelanang (Gebiidea), Nihonotrypaea thermophilus and Neaxius glyptocercus (Axiidea). All five genomes encode a standard set of 13 protein-coding genes, two ribosomal RNA genes, 22 transfer RNA genes and a putative control region. Except for T. kelanang, mud shrimp mitochondrial genomes exhibited rearrangements and novel patterns compared to the pancrustacean ground pattern. Each of the two Gebiidea species (A. edulis and U. major) and two Axiidea species (N. glyptocercus and N. thermophiles) share unique gene order specific to their infraorders and analyses further suggest these two derived gene orders have evolved independently. Phylogenetic analyses based on the concatenated nucleotide and amino acid sequences of 13 protein-coding genes indicate the possible polyphyly of mud shrimps, supporting the division of the group into two infraorders. However, the infraordinal relationships among the Gebiidea and Axiidea, and other reptants are poorly resolved. The inclusion of mt genome from more taxa, in particular the reptant infraorders Polychelida and Glypheidea is required in further analysis.

CONCLUSIONS

Phylogenetic analyses on the mt genome sequences and the distinct gene orders provide further evidences for the divergence between the two mud shrimp infraorders, Gebiidea and Axiidea, corroborating previous molecular phylogeny and justifying their infraordinal status. Mitochondrial genome sequences appear to be promising markers for resolving phylogenetic issues concerning decapod crustaceans that warrant further investigations and our present study has also provided further information concerning the mt genome evolution of the Decapoda.

The mitochondrial genomes of Euphausia pacifica and Thysanoessa raschii sequenced using 454 next-generation sequencing, with a phylogenetic analysis of their position in the Malacostracan family tree

Euphausiid krill play a critical role in coastal and oceanic food webs, linking primary producers to upper trophic levels. In addition, some species support commercial fisheries worldwide. Despite their ecological importance, the genetics of these important species remain poorly described. To improve our understanding of the genetics of these ecological links, we sequenced the mitochondrial genomes of two species of North Pacific krill, Euphausia pacifica and Thysanoessa raschii, using long-range PCR and 454 GS Junior next-generation sequencing technology. The E. pacifica mitogenome (14,692 + base pairs (bp)) encodes 13 protein-coding genes (PCGs), two ribosomal RNA (rRNA) genes, and at least 22 transfer RNA (tRNA) genes. The T. raschii mitogenome (14,240 + bp) encodes 13 PCGs, two rRNA genes, and at least 19 tRNA genes. The gene order in both species is similar to that of E. superba. Comparisons between Bering Sea and Yellow Sea E. pacifica revealed a total of 644 variable sites. The most variable protein-coding gene were atp8 (7.55 %, 12 of 159 sites variable), nad4 (6.35 %, 85 variable sites) and nad6 (6.32 %, 33 variable sites). Phylogenetic analyses to assess the phylogenetic position of the Euphausiacea, using the concatenated nucleic acid sequences of E. pacifica and T. raschii along with 46 previously published malacostracan mitogenomes, support the monophyly of the order Decapoda and indicate that the Euphausiacea share a common ancestor with the Decapoda. Future research should utilize this sequence data to explore the population genetics and molecular ecology of these species.

The complete mitochondrial genome of the Japanese mud shrimp Upogebia major (Crustacea, Decapoda)

We determined a full-length sequence of mitochondrial (mt) genome from Upogebia major. This is the first complete mt genome report for infraorder Thalassinidea in Decapoda, Crustacea. Our result showed that U. major generally followed a typical pancrustacean gene order but some tRNA genes showed a very unique gene arrangement such as duplication or translocation. Since none of the complete mt genome sequences in the infraorder Thalassinidea are available yet, this report will provide additional information in relation to mt genome diversity and evolution of the decapods.

Determination of the complete mitochondrial DNA sequence of Octopus minor

In this study, we have determined the complete nucleotide sequence of the mitochondrial genome of Octopus minor. It is 15,974 nucleotide pairs and encodes 13 proteins, two ribosomal RNAs and 22 tRNAs of the mitochondrion's own protein synthesizing system. Seven of thirteen proteins are encoded by the H-strand, while the other six proteins, as well as the two ribosomal RNAs are encoded by the L-strand. The nucleotide composition of the proteins showed a nucleotide bias against G encoded by the H-strand, while they showed a nucleotide bias against A and C encoded by the L-strand. Two of the 13 protein coding genes of O. minor began with the unorthodox translation initiation codon ATA and all others use the standard ATG. In addition, six of thirteen mt proteins of O. minor have unambiguous termination codons. There are four cases where tRNA genes appear to overlap. The long noncoding region (LNCR) of O. minor was 930 nucleotides and no repeated sequences were found in this LNCR. The gene arrangements of O. minor showed remarkable similarity to that of O. ocellatus and O. vulgaris. Phylogenetic analysis demonstrated that O. minor appears as sister taxan to the monophyletic group combined by O. ocellatus and O. vulgaris, suggesting a relative distant genetic relationship between O. minor and the other two octopus species.