The complete mitochondrial genome sequence of the hydrothermal vent galatheid crab Shinkaia crosnieri (Crustacea: Decapoda: Anomura): a novel arrangement and incomplete tRNA suite

A comprehensive study on the mitochondrial genome of Volva habei and exploring phylogenetic relationships in Littorinimorpha

In order to enrich our taxonomic and systematic comprehension of Ovulidae within the evolutionary framework of Littorinimorpha. we present a comprehensive analysis of the mitochondrial genome (mitogenome) sequence of Volva habei using next-generation sequencing technology (GenBank accession number OR492307). The mitogenome spans a total length of 16,519 bp, encompassing a complete set of 37 genes, including 13 protein-coding genes (PCGs), 22 tRNAs and two rRNAs, demonstrating a distinct AT bias. Notably, trnS2 lacks a dihydrouracil (DHU) arm, thus preventing the formation of a typical secondary structure. In contrast, the remaining tRNAs exhibit a characteristic cloverleaf-like secondary structure. Comparative analysis with ancestral gastropods reveals substantial differences in three gene clusters (or genes), incorporating fifteen tRNAs and eight PCGs. Of particular significance are the observed inversions and translocations, representing the predominant types of rearrangements in V. habei. Phylogenetic analysis strongly supports the monophyletic grouping of all Littorinimorpha species, with V. habei forming a distinct Ovulidae clade. It is noteworthy that V. habei forms a sister group with Cypraeidae, collectively belonging to the Cypraeoidea. In summary, this study not only advances our comprehension of the entire mitochondrial dataset for Calyptraeoidea but also provides novel insights into the phylogenetic relationships within Littorinimorpha.

The complete mitochondrial genome of Allogalathea elegans (Adams & White, 1848) (Decapoda: Galatheidae)

The genus Allogalathea belongs to the subfamily Galatheoidea of the family Galatheidae. Here, we report a mitogenome of Allogalathea elegans (Adams & White, 1848). In this study, we obtained the complete mitochondrial genome of Allogalathea elegans by sequencing, which was 16,263 bp in length. The mitogenome contained 37 genes, including the typical set of 13 protein-coding genes (PCGs), 22 transfer RNA (tRNA) genes, and 2 Ribosomal RNA (rRNA) genes. The nucleotides A, C, G, and T distribution was 36.40%, 19.44%, 9.09%, and 35.07%, respectively. The length of the total protein-coding genes was 11,172 bp, which accounts for 68.69% of the whole mitochondrial genome. The phylogenetic result generated by IQ-Tree based on 13 PGCs showed that the infraorder Anomura is monophyletic, and the infraorder Anomura is a sister group of the infraorder Glypheidea. The discovery of the complete mitochondrial genome of A. elegans would help to conduct in-depth research on the infraorder Anomura.

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.

The Complete Mitogenome of Xeruca formosensis (Rathbun, 1921) (Crustacea: Brachyura: Ocypodidae), a Fiddler Crab Endemic to Taiwan, with its Phylogenetic Position in the Family

Xeruca formosensis is a species and genus of fiddler crab endemic to Taiwan, with limited distribution in western Taiwan and the offshore Penghu Islands. This study reports the complete mitochondrial genome (mitogenome) of this species using next-generation sequencing. The mitogenome contains 15,684 bp, comprising 13 protein-coding genes, 22 tRNAs, 2 ribosomal RNAs and a 750-bp intergenic space (control region). The nucleotide composition is biased toward A+T (69.4%). A phylogenetic analysis based on the concatenated protein-coding genes showed that the genera Xeruca Shih, 2015 and Tubuca Bott, 1973 are sister to each other. In addition, the phylogeny of the 16 available mitogenomes in the family Ocypodidae also supports the current systematics of this family based on one nuclear and two mitochondrial markers. As this species inhabits high intertidal mudflats with high temperature and high salinity, mitogenome analyses may help us understand the mechanisms of adaptation to extreme environments, as well as the connectivity of metapopulations based on mitogenomes from different populations.

The mitochondrial genomes of five spring and groundwater amphipods of the family Crangonyctidae (Crustacea: Amphipoda) from eastern North America

We sequenced the mitochondrial genomes of one spring-dwelling (Crangonyx forbesi) and four groundwater amphipods (Bactrurus brachycaudus, Stygobromus allegheniensis, S. pizzinii, and S. t. potomacus) from eastern North America using a shotgun sequencing approach on an Illumina HiSeq 4000 (Illumina, San Diego, CA). All five mitochondrial genomes encoded 13 protein-coding genes, 22 transfer RNAs (tRNAs), and two ribosomal RNAs (rRNAs) representative of subphylum Crustacea. Although the four groundwater species exhibited gene orders nearly identical to the ancestral pancrustacean gene order, the spring-dwelling species, C. forbesi, possessed a transposition of the trnH–nad4–nad4l loci downstream after nad6–cytb–trnS2. Moreover, a long nad5 locus, longer rrnL, and rrnS loci, and unconventional start codons distinguished C. forbesi from the four groundwater amphipods. Overall, our five amphipod mitogenomes add to the increasing publicly available mitogenome resources for amphipods that are not only valuable for studying the evolutionary relationships of this diverse group of crustaceans but for exploring the evolution of mitochondrial genomes in general.

Comparative Mitogenomics in Hyalella (Amphipoda: Crustacea)

We present the sequencing and comparative analysis of 17 mitochondrial genomes of Nearctic and Neotropical amphipods of the genus Hyalella, most from the Andean Altiplano. The mitogenomes obtained comprised the usual 37 gene-set of the metazoan mitochondrial genome showing a gene rearrangement (a reverse transposition and a reversal) between the North and South American Hyalella mitogenomes. Hyalella mitochondrial genomes show the typical AT-richness and strong nucleotide bias among codon sites and strands of pancrustaceans. Protein-coding sequences are biased towards AT-rich codons, with a preference for leucine and serine amino acids. Numerous base changes (539) were found in tRNA stems, with 103 classified as fully compensatory, 253 hemi-compensatory and the remaining base mismatches and indels. Most compensatory Watson–Crick switches were AU -> GC linked in the same haplotype, whereas most hemi-compensatory changes resulted in wobble GU and a few AC pairs. These results suggest a pairing fitness increase in tRNAs after crossing low fitness valleys. Branch-site level models detected positive selection for several amino acid positions in up to eight mitochondrial genes, with atp6 and nad5 as the genes displaying more sites under selection.

The complete mitochondrial genome of a cold seep gastropod Phymorhynchus buccinoides (Neogastropoda: Conoidea: Raphitomidae)

Phymorhynchus is a genus of deep-sea snails that are most distributed in hydrothermal vent or cold seep environments. In this study, we presented the complete mitochondrial genome of P. buccinoides, a cold seep snail from the South China Sea. It is the first mitochondrial genome of a cold seep member of the superfamily Conoidea. The mitochondrial genome is 15,764 bp in length, and contains 13 protein-coding genes (PCGs), 2 rRNA genes, and 22 tRNA genes. These genes are encoded on the positive strand, except for 8 tRNA genes that are encoded on the negative strand. The start codon ATG and 3 types of stop codons, TAA, TAG and the truncated termination codon T, are used in the 13 PCGs. All 13 PCGs in the 26 species of Conoidea share the same gene order, while several tRNA genes have been translocated. Phylogenetic analysis revealed that P. buccinoides clustered with Typhlosyrinx sp., Eubela sp., and Phymorhynchus sp., forming the Raphitomidae clade, with high support values. Positive selection analysis showed that a residue located in atp6 (18 S) was identified as the positively selected site with high posterior probabilities, suggesting potential adaption to the cold seep environment. Overall, our data will provide a useful resource on the evolutionary adaptation of cold seep snails for future studies.

The complete mitochondrial genomes of two vent squat lobsters, Munidopsis lauensis and M. verrilli: Novel gene arrangements and phylogenetic implications

Hydrothermal vents are considered as one of the most extremely harsh environments on the Earth. In this study, the complete mitogenomes of hydrothermal vent squat lobsters, Munidopsis lauensis and M. verrilli, were determined through Illumina sequencing and compared with other available mitogenomes of anomurans. The mitogenomes of M. lauensis (17,483 bp) and M. verrilli (17,636 bp) are the largest among all Anomura mitogenomes, while the A+T contents of M. lauensis (62.40%) and M. verrilli (63.99%) are the lowest. The mitogenomes of M. lauensis and M. verrilli display novel gene arrangements, which might be the result of three tandem duplication-random loss (tdrl) events from the ancestral pancrustacean pattern. The mitochondrial gene orders of M. lauensis and M. verrilli shared the most similarities with S. crosnieri. The phylogenetic analyses based on both gene order data and nucleotide sequences (PCGs and rRNAs) revealed that the two species were closely related to Shinkaia crosnieri. Positive selection analysis revealed that eighteen residues in seven genes (atp8, Cytb, nad3, nad4, nad4l, nad5, and nad6) of the hydrothermal vent anomurans were positively selected sites.

Novel gene rearrangement in the mitochondrial genome of Coenobita brevimanus (Anomura: Coenobitidae) and phylogenetic implications for Anomura

The complete mitochondrial genomes (mitogenomes) can indicate phylogenetic relationships among organisms, as well as useful information about the process of molecular evolution and gene rearrangement mechanisms. However, knowledge on the complete mitogenome of Coenobitidae (Decapoda: Anomura) is quite scarce. Here, we describe in detail the complete mitogenome of Coenobita brevimanus, which is 16,393 bp in length, and contains 13 protein-coding genes, two ribosomal RNA, 22 transfer RNA genes, as well as a putative control region. The genome composition shows a moderate A + T bias (65.0%), and exhibited a negative AT-skew (-0.148) and a positive GC-skew (0.183). Five gene clusters (or genes) involving eleven tRNAs and two PCGs were found to have rearranged with respect to the pancrustacean ground pattern gene order. Duplication-random loss and recombination models were determined as most likely to explain the observed large-scale gene rearrangements. Phylogenetic analysis placed all Coenobitidae species into one clade. The polyphyly of Paguroidea was well supported, whereas the non-monophyly of Galatheoidea was inconsistence with previous findings on Anomura. Taken together, our results help to better understand gene rearrangement process and the evolutionary status of C. brevimanus and lay a foundation for further phylogenetic studies of Anomura.

Large-scale mitochondrial gene rearrangements in the hermit crab Pagurus nigrofascia and phylogenetic analysis of the Anomura

Complete mitochondrial genome (mitogenome) provides important information for better understanding of gene rearrangement, molecular evolution and phylogenetic analysis. Currently, only a few Paguridae mitogenomes have been reported. Herein, we described the complete mitogenome of hermit crab Pagurus nigrofascia. The total length was 15,423 bp, containing 13 protein-coding genes (PCGs), two ribosomal RNA, 22 transfer RNA genes, as well as an AT-rich region. The genome composition was highly A + T biased (71.4%), and exhibited a negative AT-skew (-0.006) and GC-skew (-0.138). Eight tRNA genes, two PCGs and an AT-rich region found to be rearranged with respect to the pancrustacean ground pattern gene order. Duplication-random loss and recombination model were adopted to explain the large-scale gene rearrangement events. Two phylogenetic trees of Anomura involving 12 families were constructed. The results showed that all Paguridae species were clustered into one clade except Pagurus longicarpus, which for the first time imposed raises doubt about the morphological taxonomy of this species. Furthermore, the present study found that higher- level phylogenetic relationships within Anomura were controversial, compared with the previous studies. Our results help to better understand gene rearrangements and the evolutionary status of P. nigrofascia and lay foundation for further phylogenetic study of Anomura.

Basal position of two new complete mitochondrial genomes of parasitic Cymothoida (Crustacea: Isopoda) challenges the monophyly of the suborder and phylogeny of the entire order

BACKGROUND

Isopoda is a highly diverse order of crustaceans with more than 10,300 species, many of which are parasitic. Taxonomy and phylogeny within the order, especially those of the suborder Cymothoida Wägele, 1989, are still debated. Mitochondrial (mt) genomes are a useful tool for phylogenetic studies, but their availability for isopods is very limited. To explore these phylogenetic controversies on the mt genomic level and study the mt genome evolution in Isopoda, we sequenced mt genomes of two parasitic isopods, Tachaea chinensis Thielemann, 1910 and Ichthyoxenos japonensis Richardson, 1913, belonging to the suborder Cymothoida, and conducted comparative and phylogenetic mt genomic analyses across Isopoda.

RESULTS

The complete mt genomes of T. chinensis and I. japonensis were 14,616 bp and 15,440 bp in size, respectively, with the A+T content higher than in other isopods (72.7 and 72.8%, respectively). Both genomes code for 13 protein-coding genes, 21 transfer RNA genes (tRNAs), 2 ribosomal RNA genes (rRNAs), and possess a control region (CR). Both are missing a gene from the complete tRNA set: T. chinensis lacks trnS1 and I. japonensis lacks trnI. Both possess unique gene orders among isopods. Within the CR of I. japonensis (284 bp), we identified a repetitive region with four tandem repeats. Phylogenetic analysis based on concatenated nucleotide sequences of 13 protein-coding genes showed that the two parasitic cymothoids clustered together and formed a basal clade within Isopoda. However, another parasitic cymothoid, Gyge ovalis Shiino, 1939, formed a sister group with the suborder Limnoriidea Brandt & Poore in Poore, 2002, whereas two free-living cymothoid species were located in the derived part of the phylogram: Bathynomus sp. formed a sister group with the suborder Sphaeromatidea Wägele, 1989, and Eurydice pulchra Leach, 1815 with a clade including Bathynomus sp., Sphaeromatidea and Valvifera G. O. Sars, 1883.

CONCLUSIONS

Our results did not recover the suborders Cymothoida and Oniscidea Latreille, 1802 as monophyletic, with parasitic and free-living cymothoidans forming separate clades. Furthermore, two parasitic cymothoidans formed the sister-clade to all other isopods, separated from Epicaridea Latreille, 1825, which challenges currently prevalent isopod phylogeny. Additional mt genomes of parasitic and free-living isopods might confer a sufficient phylogenetic resolution to enable us to resolve their relationships, and ultimately allow us to better understand the evolutionary history of the entire isopod order.

The first complete mitochondrial genome of the Mariana Trench Freyastera benthophila (Asteroidea: Brisingida: Brisingidae) allows insights into the deep-sea adaptive evolution of Brisingida

Starfish (phylum Echinodermata) are ecologically important and diverse members of marine ecosystems in all of the world's oceans, from the shallow water to the hadal zone. The deep sea is recognized as an extremely harsh environment on earth. In this study, we present the mitochondrial genome sequence of Mariana Trench starfish Freyastera benthophila, and this study is the first to explore in detail the mitochondrial genome of a deep-sea member of the order Brisingida. Similar to other starfish, it contained 13 protein-coding genes, two ribosomal RNA genes, and 22 transfer RNA genes (duplication of two tRNAs: trnL and trnS). Twenty-two of these genes are encoded on the positive strand, while the other 15 are encoded on the negative strand. The gene arrangement was identical to those of sequenced starfish. Phylogenetic analysis showed the deep-sea Brisingida as a sister taxon to the traditional members of the Asteriidae. Positive selection analysis indicated that five residues (8 N and 16 I in atp8, 47 D and 196 V in nad2, 599 N in nad5) were positively selected sites with high posterior probabilities. Compared these features with shallow sea starfish, we predict that variation specifically in atp8, nad2, and nad5 may play an important role in F. benthophila's adaptation to deep-sea environment.

The complete mitochondrial genome of the alvinocaridid shrimp Shinkaicaris leurokolos (Decapoda, Caridea): Insight into the mitochondrial genetic basis of deep-sea hydrothermal vent adaptation in the shrimp

Deep-sea hydrothermal vent is one of the most extreme environments on Earth with low oxygen and high levels of toxins. Decapod species from the family Alvinocarididae have colonized and successfully adapted to this extremely harsh environment. Mitochondria plays a vital role in oxygen usage and energy metabolism, thus it may be under selection in the adaptive evolution of the hydrothermal vent shrimps. In this study, the mitochondrial genome (mitogenome) of alvinocaridid shrimp Shinkaicaris leurokolos (Kikuchi & Hashimoto, 2000) was determined through Illumina sequencing. The mitogenome of S. leurokolos was 15,903bp in length, containing 13 protein-coding genes, 2 rRNAs, and 22 tRNAs. The gene order and orientation were identical to those of sequenced alvinocaridids. It has the longest concatenated sequences of protein-coding genes, tRNAs and shortest pooled rRNAs among the alvinocaridids. The control regions (CRs) of alvinocaridid were significantly longer (p<0.01) than those of the other caridaen. The alignment of the alvinocaridid CRs revealed two conserved sequence blocks (CSBs), and each of the CSBs included a noncanonical open reading frame (ORF), which may be involved in adjusting mitochondrial energy metabolism to adapt to the hydrothermal environment. Phylogenetic analysis supported that the deep-sea hydrothermal vent shrimps may have originated from those living in shallow area. Positive selection analysis reveals the evidence of adaptive change in the mitogenome of Alvinocarididae. Thirty potentially important adaptive residues were identified, which were located in atp6, cox1, cox3, cytb and nad1-5. This study explores the mitochondrial genetic basis of hydrothermal vent adaptation in alvinocaridid for the first time, and provides valuable clues regarding the adaptation.

Complete mitochondrial genomes from two species of Chinese freshwater crabs of the genus Sinopotamon recovered using next-generation sequencing reveal a novel gene order (Brachyura, Potamidae)

Recent morphological and molecular evidence has challenged classical interpretations of eubrachyuran phylogeny and evolution. Complete mitochondrial genomes of two species of potamid freshwater crabs, Sinopotamon yaanense and Sinopotamon yangtsekiense were obtained using next-generation sequencing. The results revealed a novel gene order with translocations of a five-gene block and a tRNA gene in comparison to available brachyuran mitochondrial genomes. DNA sequence comparisons position the Potamidae, a primary freshwater crab family, outside of the clade for the traditional heterotreme families, and closer to the clade that includes the thoracotreme families of grapsoid and ocypodoid crabs. Mitogenomic comparisons using rapid next-generation sequencing and a much wider taxonomic sample are required for a high-resolution examination of the phylogenetic relationships within the Eubrachyura.

Exploring the molecular basis of adaptive evolution in hydrothermal vent crab Austinograea alayseae by transcriptome analysis

Elucidating the genetic mechanisms of adaptation to the hydrothermal vent in organisms at genomic level is significant for understanding the adaptive evolution process in the extreme environment. We performed RNA-seq on four different tissues of a vent crab species, Austinograea alayseae, producing 725,461 unigenes and 134,489 annotated genes. Genes related to sensory, circadian rhythm, hormone, hypoxia stress, metal detoxification and immunity were identified. It was noted that in the degenerated eyestalk, transcription of phototransduction related genes which are important for retinal function was greatly reduced; three crucial neuropeptide hormones, one molt-inhibiting and two crustacean hyperglycemic hormone precursors were characterized with conserved domains; hypoxia-inducible factor 1 and two novel isoforms of metallothioneins in the vent crabs were discovered. An analysis of 6,932 orthologs among three crabs A. alayseae, Portunus trituberculutus and Eriocheir sinensis revealed 19 positive selected genes (PSGs). Most of the PSGs were involved in immune responses, such as crustins and anti-lipopolysaccharide factor, suggesting their function in the adaptation to environment. The characterization of the first vent crab transcriptome provides abundant resources for genetic and evolutionary studies of this species, and paves the way for further investigation of vent adaptation process in crabs.

The first complete mitogenome of the South China deep-sea giant isopod Bathynomus sp. (Crustacea: Isopoda: Cirolanidae) allows insights into the early mitogenomic evolution of isopods

In this study, the complete mitochondrial (mt) genome sequence of the South China deep‐sea giant isopod Bathynomus sp. was determined, and this study is the first to explore in detail the mt genome of a deep‐sea member of the order Isopoda. This species belongs to the genus Bathynomus, the members of which are saprophagous residents of the deep‐sea benthic environment; based on their large size, Bathynomus is included in the “supergiant group” of isopods. The mt genome of Bathynomus sp. is 14,965 bp in length and consists of 13 protein‐coding genes, two ribosomal RNA genes, only 18 transfer RNA genes, and a noncoding control region 362 bp in length, which is the smallest control region discovered in Isopoda to date. Although the overall genome organization is typical for metazoans, the mt genome of Bathynomus sp. shows a number of derived characters, such as an inversion of 10 genes when compared to the pancrustacean ground pattern. Rearrangements in some genes (e.g., cob, trnT, nad5, and trnF) are shared by nearly all isopod mt genomes analyzed thus far, and when compared to the putative isopod ground pattern, five rearrangements were found in Bathynomus sp. Two tRNAs exhibit modified secondary structures: The TΨC arm is absent from trnQ, and trnC lacks the DHU. Within the class Malacostraca, trnC arm loss is only found in other isopods. Phylogenetic analysis revealed that Bathynomus sp. (Cymothoida) and Sphaeroma serratum (Sphaeromatidea) form a single clade, although it is unclear whether Cymothoida is monophyletic or paraphyletic. Moreover, the evolutionary rate of Bathynomus sp. (dN/dS [nonsynonymous mutational rate/synonymous mutational rate] = 0.0705) is the slowest measured to date among Cymothoida, which may be associated with its relatively constant deep‐sea environment. Overall, our results may provide useful information for understanding the evolution of deep‐sea Isopoda species.

The mitochondrial genome sequence of a deep-sea, hydrothermal vent limpet, Lepetodrilus nux, presents a novel vetigastropod gene arrangement

While mitochondrial (mt) genomes are used extensively for comparative and evolutionary genomics, few mt genomes of deep-sea species, including hydrothermal vent species, have been determined. The Genus Lepetodrilus is a major deep-sea gastropod taxon that occurs in various deep-sea ecosystems. Using next-generation sequencing, we determined nearly the complete mitochondrial genome sequence of Lepetodrilus nux, which inhabits hydrothermal vents in the Okinawa Trough. The total length of the mitochondrial genome is 16,353bp, excluding the repeat region. It contains 13 protein-coding genes, 22 tRNA genes, two rRNA genes, and a control region, typical of most metazoan genomes. Compared with other vetigastropod mt genome sequences, L. nux employs a novel mt gene arrangement. Other novel arrangements have been identified in the vetigastropod, Fissurella volcano, and in Chrysomallon squamiferum, a neomphaline gastropod; however, all three gene arrangements are different, and Bayesian inference suggests that each lineage diverged independently. Our findings suggest that vetigastropod mt gene arrangements are more diverse than previously realized.

The complete mitochondrial genome of squat lobster, Munida gregaria (Anomura, Galatheoidea, Munididae)

We determined the mitogenome sequence of Munida gregaria (Fabricius 1793) (Anomura, Galatheoidea, Munididae), which is the first complete mitogenome sequence in the family Munididae Ahyong et al., 2010. The mitogenome of M. gregaria is 16 326 bp in length and contains 13 protein-coding genes (PCGs), 22 transfer RNAs (tRNAs), two ribosomal RNAs (rRNAs) and two control regions (CRs). Mitogenome analysis of M. gregaria showed an extra copy of the CR and rearrangements of two PCGs (nad2 and nad3) compared to the arthropod ground pattern. Additionally, it contains a tRNA (trnY) inversion and rearrangements of two PCGs (nad1 and nad3) when compared with that of Neopetrolisthes maculatus and Shinkaia crosnieri, respectively. The phylogenetic tree confirmed that M. gregaria belongs to the superfamily Galatheoidea within Anomura. Our results will be useful for the detailed study of mitogenome evolution and the phylogenetic relationships among the superfamily in the infraorder Anomura.

Mitochondrial genome rearrangements at low taxonomic levels: three distinct mitogenome gene orders in the genus Pseudoniphargus (Crustacea: Amphipoda)

A comparison of mitochondrial genomes of three species of the amphipod Pseudoniphargus revealed the occurrence of a surprisingly high level of gene rearrangement involving protein-coding genes that is a rare phenomenon at low taxonomic levels. The three Pseudoniphargus mitogenomes also display a unique gene arrangement with respect to either the presumed Pancrustacean order or those known for other amphipods. Relative long non-coding sequences appear adjacent to the putative breakage points involved in gene rearrangements of protein coding genes. Other details of the newly obtained mitochondrial genomes - e.g., gene content, nucleotide composition and codon usage - are similar to those found in the mitogenomes of other amphipod species studied. They all contain the typical mitochondrial genome set consisting of 13 protein-coding genes, 22 tRNAs, and two rRNAS, as well as a large control region. The secondary structures and characteristics of tRNA and ribosomal mitochondrial genes of these three species are also discussed.

Adaptations to Hydrothermal Vent Life in Kiwa tyleri, a New Species of Yeti Crab from the East Scotia Ridge, Antarctica

Hydrothermal vents in the Southern Ocean are the physiologically most isolated chemosynthetic environments known. Here, we describe Kiwa tyleri sp. nov., the first species of yeti crab known from the Southern Ocean. Kiwa tyleri belongs to the family Kiwaidae and is the visually dominant macrofauna of two known vent sites situated on the northern and southern segments of the East Scotia Ridge (ESR). The species is known to depend on primary productivity by chemosynthetic bacteria and resides at the warm-eurythermal vent environment for most of its life; its short-range distribution away from vents (few metres) is physiologically constrained by the stable, cold waters of the surrounding Southern Ocean. Kiwa tylerihas been shown to present differential life history adaptations in response to this contrasting thermal environment. Morphological adaptations specific to life in warm-eurythermal waters, as found on – or in close proximity of – vent chimneys, are discussed in comparison with adaptations seen in the other two known members of the family (K. hirsuta, K. puravida), which show a preference for low temperature chemosynthetic environments.

The complete mitogenome of the Norway lobster Nephrops norvegicus (Linnaeus, 1758) (Crustacea: Decapoda: Nephropidae)

The clawed lobster Nephrops norvegicus is an important commercial species in European waters. We have sequenced the complete mitochondrial genome of the species from a partial genome scan using Next-Gen sequencing. The N. norvegicus has a mitogenome of 16,132 base pairs (71.22% A+ T content) comprising 13 protein-coding genes, 2 ribosomal subunit genes, 21 transfer RNAs, and a putative 1259 bp non-coding AT-rich region. This mitogenome is the second fully characterized for the family Nephropidae and the first for the genus Nephrops. The mitogenome gene order is identical to the Maine lobster, Homarus americanus with the exception of the possible loss of the trnI gene.

Keeping mtDNA in shape between generations

Since the unexpected discovery that mitochondria contain their own distinct DNA molecules, studies of the mitochondrial DNA (mtDNA) have yielded many surprises. In animals, transmission of the mtDNA genome is explicitly non-Mendelian, with a very high number of genome copies being inherited from the mother after a drastic bottleneck. Recent work has begun to uncover the molecular details of this unusual mode of transmission. Many surprising variations in animal mitochondrial biology are known; however, a series of recent studies have identified a core of evolutionarily conserved mechanisms relating to mtDNA inheritance, e.g., mtDNA bottlenecks during germ cell development, selection against specific mtDNA mutation types during maternal transmission, and targeted destruction of sperm mitochondria. In this review, we outline recent literature on the transmission of mtDNA in animals and highlight the implications for human health and ageing.

Complete mitochondrial genome of Caridina nilotica gracilipes

In this study, we sequenced the complete mitochondrial genome of Caridina nilotica gracilipes. This mitochondrial genome, consisting of 15,550 base pairs, encoded 13 protein-coding genes, 2 ribosomal RNAs, 22 transfer RNAs, and a non-coding control region as those found in other Decapoda, with the gene synteny identical to that of typical invertebrates. Control region (D-Loop), of 673 bp in length, is located between 12S rRNA and tRNA(Ile). The overall base composition of the heavy strand shows T 30.4%, C 22.0%, A 33.0% and G 14.6%, with an AT bias of 63.4%.

Next-generation sequencing, phylogenetic signal and comparative mitogenomic analyses in Metacrangonyctidae (Amphipoda: Crustacea)

Background

Comparative mitochondrial genomic analyses are rare among crustaceans below the family or genus level. The obliged subterranean crustacean amphipods of the family Metacrangonyctidae, found from the Hispaniola (Antilles) to the Middle East, including the Canary Islands and the peri-Mediterranean region, have an evolutionary history and peculiar biogeography that can respond to Tethyan vicariance. Indeed, recent phylogenetic analysis using all protein-coding mitochondrial sequences and one nuclear ribosomal gene have lent support to this hypothesis (Bauzà-Ribot et al. 2012).

Results

We present the analyses of mitochondrial genome sequences of 21 metacrangonyctids in the genera Metacrangonyx and Longipodacrangonyx, covering the entire geographical range of the family. Most mitogenomes were attained by next-generation sequencing techniques using long-PCR fragments sequenced by Roche FLX/454 or GS Junior pyro-sequencing, obtaining a coverage depth per nucleotide of up to 281×. All mitogenomes were AT-rich and included the usual 37 genes of the metazoan mitochondrial genome, but showed a unique derived gene order not matched in any other amphipod mitogenome. We compare and discuss features such as strand bias, phylogenetic informativeness, non-synonymous/synonymous substitution rates and other mitogenomic characteristics, including ribosomal and transfer RNAs annotation and structure.

Conclusions

Next-generation sequencing of pooled long-PCR amplicons can help to rapidly generate mitogenomic information of a high number of related species to be used in phylogenetic and genomic evolutionary studies. The mitogenomes of the Metacrangonyctidae have the usual characteristics of the metazoan mitogenomes (circular molecules of 15,000-16,000 bp, coding for 13 protein genes, 22 tRNAs and two ribosomal genes) and show a conserved gene order with several rearrangements with respect to the presumed Pancrustacean ground pattern. Strand nucleotide bias appears to be reversed with respect to the condition displayed in the majority of crustacean mitogenomes since metacrangonyctids show a GC-skew at the (+) and (-) strands; this feature has been reported also in the few mitogenomes of Isopoda (Peracarida) known thus far. The features of the rRNAs, tRNAs and sequence motifs of the control region of the Metacrangonyctidae are similar to those of the few crustaceans studied at present.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-566) contains supplementary material, which is available to authorized users.

The complete mitochondrial genomes of two octopods Cistopus chinensis and Cistopus taiwanicus: revealing the phylogenetic position of the genus Cistopus within the order Octopoda

In the present study, we determined the complete mitochondrial DNA (mtDNA) sequences of two species of Cistopus, namely C. chinensis and C. taiwanicus, and conducted a comparative mt genome analysis across the class Cephalopoda. The mtDNA length of C. chinensis and C. taiwanicus are 15706 and 15793 nucleotides with an AT content of 76.21% and 76.5%, respectively. The sequence identity of mtDNA between C. chinensis and C. taiwanicus was 88%, suggesting a close relationship. Compared with C. taiwanicus and other octopods, C. chinensis encoded two additional tRNA genes, showing a novel gene arrangement. In addition, an unusual 23 poly (A) signal structure is found in the ATP8 coding region of C. chinensis. The entire genome and each protein coding gene of the two Cistopus species displayed notable levels of AT and GC skews. Based on sliding window analysis among Octopodiformes, ND1 and DN5 were considered to be more reliable molecular beacons. Phylogenetic analyses based on the 13 protein-coding genes revealed that C. chinensis and C. taiwanicus form a monophyletic group with high statistical support, consistent with previous studies based on morphological characteristics. Our results also indicated that the phylogenetic position of the genus Cistopus is closer to Octopus than to Amphioctopus and Callistoctopus. The complete mtDNA sequence of C. chinensis and C. taiwanicus represent the first whole mt genomes in the genus Cistopus. These novel mtDNA data will be important in refining the phylogenetic relationships within Octopodiformes and enriching the resource of markers for systematic, population genetic and evolutionary biological studies of Cephalopoda.

A surprising arrangement pattern and phylogenetic consideration: the complete mitochondrial genome of Belanger's croaker Johnius belangerii (Percoidei: Sciaenidae)

The complete mitochondrial genome of Johnius belangerii has been determined for the first time in this article. It was 19,154 base pairs in length, and is composed of 37 genes (13 protein-coding genes, 22 tRNA genes and 2 ribosomal RNA genes). Totally, 5 notable non-coding regions were observed, and a non-coding of 1091 bp was identified as control region based on its location and AT richness. An 800 bp tandem repeat sequence was identified in the fifth non-coding region. We investigated the mitochondrial gene arrangement pattern and found that that the tRNA(Val), 12SrRNA, 16SrRNA and tRNA(Phe) genes of J. belangerii mitogenome were orderly placed at the beginning of heavy strand. This order is different from other croakers. Combine with the phylogenetic reconstruction and genes arrangement pattern of J. belangerii mitochondrial genome, we consider J. belangerii is the most ancient genus within family Sciaenidae.

Comparative analysis of mitochondrial genomes of five aphid species (Hemiptera: Aphididae) and phylogenetic implications

Insect mitochondrial genomes (mitogenomes) are of great interest in exploring molecular evolution, phylogenetics and population genetics. Only two mitogenomes have been previously released in the insect group Aphididae, which consists of about 5,000 known species including some agricultural, forestry and horticultural pests. Here we report the complete 16,317 bp mitogenome of Cavariella salicicola and two nearly complete mitogenomes of Aphis glycines and Pterocomma pilosum. We also present a first comparative analysis of mitochondrial genomes of aphids. Results showed that aphid mitogenomes share conserved genomic organization, nucleotide and amino acid composition, and codon usage features. All 37 genes usually present in animal mitogenomes were sequenced and annotated. The analysis of gene evolutionary rate revealed the lowest and highest rates for COI and ATP8, respectively. A unique repeat region exclusively in aphid mitogenomes, which included variable numbers of tandem repeats in a lineage-specific manner, was highlighted for the first time. This region may have a function as another origin of replication. Phylogenetic reconstructions based on protein-coding genes and the stem-loop structures of control regions confirmed a sister relationship between Cavariella and pterocommatines. Current evidence suggest that pterocommatines could be formally transferred into Macrosiphini. Our paper also offers methodological instructions for obtaining other Aphididae mitochondrial genomes.

The mitochondrial genome of the scallop Mimachlamys senatoria (Bivalvia, Pectinidae)

The mitochondrial (mt) genome of the scallop Mimachlamys senatoria (17,383 bp), an economically and ecologically important bivalve, was newly sequenced and annotated. Comparative analyses between M. senatoria and its congeneric sister species M. nobilis revealed three new findings: (1) M. senatoria is more prone to use G-rich start/stop codon, and variation in start/stop codon usage is species-correlated rather than gene-correlated, and in some extent, bears useful phylogenetic information; (2) The A + T content is unexpectedly low (54.1%) in MNR and that is unexpectedly high (65.4%) in atp8 in both congeneric scallops, which may represent a novel evolutionary pattern of mt genomic nucleotide composition; and (3) The tRNA gene cluster "NGV" locating upstream of the nad1 in M. senatoria is replaced by "NTGV" in M. nobilis, and a parsimonious explanation for the existence of trnT in M. nobilis is that this gene was derived from a recently duplicated trnG gene via an alloacceptor tRNA gene recruitment process.

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 complete mitochondrial genome of the black mud crab, Scylla serrata (Crustacea: Brachyura: Portunidae) and its phylogenetic position among (pan)crustaceans

The black mud crab, Scylla serrata (Forskål 1775), is the most economically important edible crab in South-East Asia. In the present study, the complete mitochondrial genome of black mud crab, S. serrata, was determined with the sequential polymerase chain reaction and primer walking sequencing. The complete mitochondrial genome was 15,721 bp in length with an A+T content of 69.2 % and contained 37 mitochondrial genes (13 protein coding genes (PCGs), 2 ribosomal RNA genes and 22 transfer RNA genes) and a control region (CR). The analysis of the CR sequence shows that it contains a multitude of repetitive fragments which can fold into hairpin-like or secondary structures and conserved elements as in other arthropods. The gene order of S. serrata mainly retains as the pancrustacean ground pattern, except for a single translocation of trnH. The gene arrangement of S. serrata appears to be a typical feature of portunid crabs. Phylogenetic analyses with concatenated amino acid sequences of 12 PCGs establishes that S. serrata in a well-supported monophyletic Portunidae and is consistent with previous morphological classification. Moreover, the phylogenomic results strongly support monophyletic Pancrustacea (Hexapoda plus "Crustaceans"). Within Pancrustacea, this study identifies Malacostraca + Entomostraca and Branchiopoda as the sister group to Hexapoda, which confirms that "Crustacea" is not monophyletic. Cirripedia + Remipedia appear to be a basal lineage of Pancrustacea. The present study also provides considerable data for the application of both population and phylogenetic studies of other crab species.

Complete mitogenome of the deep-sea hydrothermal vent shrimp Alvinocaris chelys Komai and Chan, 2010 (Decapoda: Caridea: Alvinocarididae)

We reported the complete sequence of the mitochondrial (mt) genome of the deep-sea vent shrimp Alvinocaris chelys. This is the first complete mt genome obtained for the hydrothermal vent shrimp. The gene arrangement of the A. chelys mt genome is identical to the pancrustacean ground pattern as in most of the other caridean shrimp mt genomes available to date. However, there is an exceptionally long intergene spacer (86 bp in length) existed between the ND1 and tRNA(Leu)-CUN genes that is not previously reported. Our results provide further evidence that the mt gene order is highly conserved among caridean shrimp, in contrast to the other decapod infraorder such as Brachyura or Anomura which are of comparable or lower diversity.

The complete mitochondrial genome of Galba pervia (Gastropoda: Mollusca), an intermediate host snail of Fasciola spp

Complete mitochondrial (mt) genomes and the gene rearrangements are increasingly used as molecular markers for investigating phylogenetic relationships. Contributing to the complete mt genomes of Gastropoda, especially Pulmonata, we determined the mt genome of the freshwater snail Galba pervia, which is an important intermediate host for Fasciola spp. in China. The complete mt genome of G. pervia is 13,768 bp in length. Its genome is circular, and consists of 37 genes, including 13 genes for proteins, 2 genes for rRNA, 22 genes for tRNA. The mt gene order of G. pervia showed novel arrangement (tRNA-His, tRNA-Gly and tRNA-Tyr change positions and directions) when compared with mt genomes of Pulmonata species sequenced to date, indicating divergence among different species within the Pulmonata. A total of 3655 amino acids were deduced to encode 13 protein genes. The most frequently used amino acid is Leu (15.05%), followed by Phe (11.24%), Ser (10.76%) and IIe (8.346%). Phylogenetic analyses using the concatenated amino acid sequences of the 13 protein-coding genes, with three different computational algorithms (maximum parsimony, maximum likelihood and Bayesian analysis), all revealed that the families Lymnaeidae and Planorbidae are closely related two snail families, consistent with previous classifications based on morphological and molecular studies. The complete mt genome sequence of G. pervia showed a novel gene arrangement and it represents the first sequenced high quality mt genome of the family Lymnaeidae. These novel mtDNA data provide additional genetic markers for studying the epidemiology, population genetics and phylogeographics of freshwater snails, as well as for understanding interplay between the intermediate snail hosts and the intra-mollusca stages of Fasciola spp..

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.

Complete mitochondrial DNA sequence of Stenopus hispidus (Crustacea: Decapoda: Stenopodidea) and a novel tRNA gene cluster

As a phylogenetically valuable decapoda, a complete mitochondrial genome from Stenopodidea has not been reported to date. Here, we determined the complete mitochondrial DNA sequence of Stenopus hispidus (Olivier, 1811). The 15,528 bp genome is a circular molecule and consists of 13 protein-coding genes (PCGs) and two ribosomal RNA (rRNA) genes plus the putative control region (CR). This finding is similar to other metazoan animals but with the exception of 23 transfer RNA (tRNA) genes, which contain an additional tRNA-Gln compared with other crustaceans. With respect to the pancrustacean ground pattern mitochondria gene order, 5 tRNAs appear to be rearranged (tRNAs-Leu (CUN), Arg, Glu, Gln, and Met), one of which has also undergone inversion (tRNA-Leu (CUN)). Phylogenetic analyses reveal Stenopodidea and Reptantia form a clade sister to Caridea, which agrees with Abele and Felgenhauer's (1986) hypothesis. This topology contrasts with previous results based on morphological and some molecular data.

Dancing for food in the deep sea: bacterial farming by a new species of Yeti crab

Vent and seep animals harness chemosynthetic energy to thrive far from the sun's energy. While symbiont-derived energy fuels many taxa, vent crustaceans have remained an enigma; these shrimps, crabs, and barnacles possess a phylogenetically distinct group of chemosynthetic bacterial epibionts, yet the role of these bacteria has remained unclear. We test whether a new species of Yeti crab, which we describe as Kiwa puravida n. sp, farms the epibiotic bacteria that it grows on its chelipeds (claws), chelipeds that the crab waves in fluid escaping from a deep-sea methane seep. Lipid and isotope analyses provide evidence that epibiotic bacteria are the crab's main food source and K. puravida n. sp. has highly-modified setae (hairs) on its 3(rd) maxilliped (a mouth appendage) which it uses to harvest these bacteria. The ε- and γ- proteobacteria that this methane-seep species farms are closely related to hydrothermal-vent decapod epibionts. We hypothesize that this species waves its arm in reducing fluid to increase the productivity of its epibionts by removing boundary layers which may otherwise limit carbon fixation. The discovery of this new species, only the second within a family described in 2005, stresses how much remains undiscovered on our continental margins.

The complete mitochondrial genome sequence of Euphausia pacifica (Malacostraca: Euphausiacea) reveals a novel gene order and unusual tandem repeats

Euphausiid krill are dominant organisms in the zooplankton population and play a central role in marine ecosystems. Euphausia pacifica (Malacostraca: Euphausiacea) is one of the most important and dominant crustaceans in the North Pacific Ocean. In this paper, we described the gene content, organization, and codon usage of the E. pacifica mitochondrial genome. The mitochondrial genome of E. pacifica is 16 898 bp in length and contains a standard set of 13 protein-coding genes, 2 ribosomal RNA genes, and 22 transfer RNA genes. Translocation of three transfer RNAs (trnL1, trnL2, and trnW) was found in the E. pacifica mitochondrial genome when comparing with the pancrustacean ground pattern. The rate of Ka/Ks in 13 protein-coding genes among three krill is much less than 1, which indicates a strong purifying selection within this group. The largest noncoding region in the E. pacifica mitochondrial genome contains one section with tandem repeats (4.7 × 154 bp), which are the largest tandem repeats found in malacostracan mitochondrial genomes so far. All analyses based on nucleotide and amino acid data strongly support the monophyly of Stomatopoda, Penaeidae, Caridea, Brachyura, and Euphausiacea. The Bayesian analysis of nucleotide and amino acid datasets strongly supports the close relationship between Euphausiacea and Decapoda, which confirms traditional findings. The maximum likelihood analysis based on amino acid data strongly supports the close relationship between Euphausiacea and Penaeidae, which destroys the monophyly of Decapoda.

The complete mitochondrial genome of the Asiatic cavity-nesting honeybee Apis cerana (Hymenoptera: Apidae)

In the present study, we determined the complete mitochondrial DNA (mtDNA) sequence of Apis cerana, the Asiatic cavity-nesting honeybee. We present here an analysis of features of its gene content and genome organization in comparison with Apis mellifera to assess the variation within the genus Apis and among main groups of Hymenoptera. The size of the entire mt genome of A. cerana is 15,895 bp, containing 2 ribosomal RNA genes, 13 protein-coding genes, 22 transfer RNA (tRNA) genes and one control region. These genes are transcribed from both strands and have a nucleotide composition high in A and T. The contents of A+T of the complete genomes are 83.96% for A. cerana. The AT bias had a significant effect on both the codon usage pattern and amino acid composition of proteins. There are a total of 3672 codons in all 13 protein-coding genes, excluding termination codons. The most frequently used amino acid is Leu (15.52%), followed by Ile (12.85%), Phe (10.10%), Ser (9.15%) and Met (8.96%). Intergenic regions in the mt genome of A. cerana are 705 bp in total. The order and orientation of the gene arrangement pattern is identical to that of A. mellifera, except for the position of the tRNA-Ser(AGN) gene. Phylogenetic analyses using concatenated amino acid sequences of 13 protein-coding genes, with three different computational algorithms (NJ, MP and ML), all revealed two distinct groups with high statistical support, indicating that A. cerana and A. mellifera are two separate species, consistent with results of previous morphological and molecular studies. The complete mtDNA sequence of A. cerana provides additional genetic markers for studying population genetics, systematics and phylogeographics of honeybees.

The complete mitochondrial genome of the Asiatic cavity-nesting honeybee Apis cerana (Hymenoptera: Apidae)

In the present study, we determined the complete mitochondrial DNA (mtDNA) sequence of Apis cerana, the Asiatic cavity-nesting honeybee. We present here an analysis of features of its gene content and genome organization in comparison with Apis mellifera to assess the variation within the genus Apis and among main groups of Hymenoptera. The size of the entire mt genome of A. cerana is 15,895 bp, containing 2 ribosomal RNA genes, 13 protein-coding genes, 22 transfer RNA (tRNA) genes and one control region. These genes are transcribed from both strands and have a nucleotide composition high in A and T. The contents of A+T of the complete genomes are 83.96% for A. cerana. The AT bias had a significant effect on both the codon usage pattern and amino acid composition of proteins. There are a total of 3672 codons in all 13 protein-coding genes, excluding termination codons. The most frequently used amino acid is Leu (15.52%), followed by Ile (12.85%), Phe (10.10%), Ser (9.15%) and Met (8.96%). Intergenic regions in the mt genome of A. cerana are 705 bp in total. The order and orientation of the gene arrangement pattern is identical to that of A. mellifera, except for the position of the tRNA-Ser(AGN) gene. Phylogenetic analyses using concatenated amino acid sequences of 13 protein-coding genes, with three different computational algorithms (NJ, MP and ML), all revealed two distinct groups with high statistical support, indicating that A. cerana and A. mellifera are two separate species, consistent with results of previous morphological and molecular studies. The complete mtDNA sequence of A. cerana provides additional genetic markers for studying population genetics, systematics and phylogeographics of honeybees.