The complete mitochondrial genomes of two common shrimps (Litopenaeus vannamei and Fenneropenaeus chinensis) and their phylogenomic considerations

Characterization of the Complete Mitochondrial Genome of the Spotted Catfish Arius maculatus (Thunberg, 1792) and Its Phylogenetic Implications

The spotted catfish, Arius maculatus (Siluriformes), is an important economical aquaculture species inhabiting the Indian Ocean, as well as the western Pacific Ocean. The bioinformatics data in previous studies about the phylogenetic reconstruction of Siluriformes were insufficient and incomplete. In the present study, we presented a newly sequenced A. maculatus mitochondrial genome (mtDNA). The A. maculatus mtDNA was 16,710 bp in length and contained two ribosomal RNA (rRNA) genes, thirteen protein-coding genes (PCGs), twenty-two transfer RNA (tRNA) genes, and one D-loop region. The composition and order of these above genes were similar to those found in most other vertebrates. The relative synonymous codon usage (RSCU) of the 13 PCGs in A. maculatus mtDNA was consistent with that of PCGs in other published Siluriformes mtDNA. Furthermore, the average non-synonymous/synonymous mutation ratio (Ka/Ks) analysis, based on the 13 PCGs of the four Ariidae species, showed a strong purifying selection. Additionally, phylogenetic analysis, according to 13 concatenated PCG nucleotide and amino acid datasets, showed that A. maculatus and Netuma thalassina (Netuma), Occidentarius platypogon (Occidentarius), and Bagre panamensis (Bagre) were clustered as sister clade. The complete mtDNA of A. maculatus provides a helpful dataset for research on the population structure and genetic diversity of Ariidae species.

The complete mitogenome of Lysmata vittata (Crustacea: Decapoda: Hippolytidae) with implication of phylogenomics and population genetics

In this study, the complete mitogenome of Lysmata vittata (Crustacea: Decapoda: Hippolytidae) has been determined. The genome sequence was 22003 base pairs (bp) and it included thirteen protein-coding genes (PCGs), twenty-two transfer RNA genes (tRNAs), two ribosomal RNA genes (rRNAs) and three putative control regions (CRs). The nucleotide composition of AT was 71.50%, with a slightly negative AT skewness (-0.04). Usually the standard start codon of the PCGs was ATN, while cox1, nad4L and cox3 began with TTG, TTG and GTG. The canonical termination codon was TAA, while nad5 and nad4 ended with incomplete stop codon T, and cox1 ended with TAG. The mitochondrial gene arrangement of eight species of the Hippolytidae were compared with the order of genes of Decapoda ancestors, finding that the gene arrangement order of the Lebbeus groenlandicus had not changed, but the gene arrangement order of other species changed to varying degrees. The positions of the two tRNAs genes (trnA and trnR) of the L. vittata had translocations, which also showed that the Hippolytidae species were relatively unconserved in evolution. Phylogenetic analysis of 50 shrimp showed that L. vittata formed a monophyletic clade with Lysmata/Exhippolysmata species. This study should be helpful to better understand the evolutionary status, and population genetic diversity of L. vittata and related species.

Characterization of Penaeus vannamei mitogenome focusing on genetic diversity

The diversity of the Penaeus vannamei mitochondrial genome has still been poorly characterized, there are no validated mitochondrial markers available for populational studies, and the heteroplasmy has not yet been investigated in this species. In this study, metagenomic reads extracted from the muscle of a single individual were used to assemble the mitochondrial genome (mtDNA). These data associated with mitochondrial genomes previously described allowed to evaluate the inter-individual variability and heteroplasmy. Comparison among 45 mtDNA control regions led to the detection of conserved and variable segments and the characterization of two hypervariable regions. The analysis of diversity revealed mostly low frequency polymorphisms, and heteroplasmy was found in practically all mitochondrial genes, with a high occurrence of indels. These results indicate that the design of mitochondrial markers for P. vannamei must be done with caution. The mapping of conserved and variable regions and the characterization of heteroplasmy presented here will contribute to increasing the efficiency of mitochondrial markers for population or individual studies.

Analysis of transcriptome difference between rapid-growing and slow-growing in Penaeus vannamei

Penaeus vannamei is the principle cultured shrimp species in China. However, with the increase of culture density, the growth difference between individuals is also expanding. Here, we make use of RNA-seq to study the growth mechanisms of P. vannamei. After 120 days, we examined the transcriptomes of rapid-growing individuals (RG) and slow-growing individuals (SG). A total of 2116 and 176 differentially expressed genes (DEGs) were found in SG and RG, respectively. Moreover, the main DEGs are opsin, heat shock protein (HSP), actin, myosin, superoxide dismutase (SOD), cuticle protein, and chitinase. GO analysis further revealed that the DEGs were enriched in biological processes significantly, such as "sensory perception," "sensory perception of light stimulus," "response to stimulus," and "response to stress." Additionally, KEGG enrichment analysis showed that the DEGs were mainly enriched in "pentose and glucuronate interconversions," "amino sugar and nucleotide sugar metabolism," "glycophospholipid biosynthesis," and "glutathione metabolism." Interestingly, the upstream genes in the ecdysone signaling pathway, including molting inhibition hormone (MIH) and crustacean hyperglycemic hormone (CHH), did not differ significantly between RG and SG, which suggests that the cause for the inconsistent growth performance is due to the stress levels rather than the ecdysone signal pathway. In summary, this work provides data that will be useful for future studies on shrimp growth and development.

A chromosome-level genome assembly of the yellowfin seabream (Acanthopagrus latus; Hottuyn, 1782) provides insights into its osmoregulation and sex reversal

The yellowfin seabream Acanthopagrus latus is the economically most important Sparidae fish in the northern South China Sea. As euryhaline fish, they are perfect model for investigating osmoregulatory mechanisms in teleosts. Moreover, the reproductive biology of hermaphrodites has long been intriguing; however, little information is known about the molecular pathways underlying their sex change. Here, we report a chromosome level reference genome of A. latus generated by employing the PacBio single molecule sequencing technique (SMRT) and high-throughput chromosome conformation capture (Hi-C) technologies. The draft genome of yellowfin seabream was 806 Mb, with 732 Mb scaffolds anchored on 24 chromosomes. The contig N50 and scaffold N50 were 2.6 Mb and 30.17 Mb, respectively. The assembly is of high integrity and includes 92.23% universal single-copy orthologues based on benchmarking universal single-copy orthologs (BUSCO) analysis. A total of 19,631 protein-coding genes were functionally annotated in the reference genome. Moreover, ARRDC3 and GSTA gene families which related to osmoregulation underwent an extensive expansion in two euryhaline sparids fish genomes compared to other teleost genomes. Moreover, integrating sex-specific transcriptome analyses, several genes related to the transforming growth factor beta (TGF-β) signalling pathway involved in sex differentiation and development. This genomic resource will not only be valuable for studying the osmoregulatory mechanisms in estuarine fish and sex determination in hermaphrodite vertebrate species, but also provide useful genomic tools for facilitating breeding of the yellowfin seabream.

Characterization of the complete mitochondrial genomes of two species of Penaeidae (Decapoda: Dendrobranchiata) and the phylogenetic implications for Penaeoidea

In the present study, mitogenomes of the species Trachypenaeus curvirostris and Parapenaeus fissuroides (Decapoda: Dendrobranchiata: Penaeidae) were sequenced. The total lengths of the two species were 15,956 bp and 15,937 bp in length with A + T biases of 67.08% and 67.69%, respectively. Both two species showed positive AT skews (0.016, 0.058) and negative GC skews (-0.254, -0.310). Both mitogenomes contained 13 protein-coding genes, 22 transfer RNA genes, and 2 ribosomal RNA genes. Results of phylogenetic analyses support close relationships among Aristeidae, Benthesicymidae and Solenoceridae. The family Sicyoniidae was observed to be deeply nested within Penaeidae. Within Penaeidae, T. curvirostris and P. fissuroides were most closely related to the genus Parapenaeopsis and Metapenaeopsis, respectively, indicated that these two species belong to Penaeidae. These results will help to better understand the evolutionary position of Penaeidae and provide reference for further phylogenetic research on Penaeoidea species.

Genetic analysis of whole mitochondrial genome of Lateolabrax maculatus (Perciformes: Moronidae) indicates the presence of two populations along the Chinese coast

The whole mitochondrial genome ofLateolabrax maculatus(Cuvier, 1828) was used to investigate the reasons for the observed patterns of genetic differentiation among 12 populations in northern and southern China. The haplotype diversity and nucleotide diversity ofL. maculatuswere 0.998 and 0.00169, respectively. Pairwise FSTvalues between populations ranged from 0.001 to 0.429, correlating positively with geographic distance. Genetic structure analysis and haplotype network analysis indicated that these populations were split into two groups, in agreement with geographic segregation and environment. Tajima’s D values, Fu’s Fs tests and Bayesian skyline plot (BSP) indicated that a demographic expansion event may have occurred in the history ofL. maculatus. Through selection pressure analysis, we found evidence of significant negative selection at the ATP6, ND3, Cytb, COX3, COX2 and COX1 genes. In our hypotheses, this study implied that demographic events and selection of local environmental conditions, including temperature, are responsible for population divergence. These findings are a step forward toward the understanding of the genetic basis of differentiation and adaptation, as well as conservation ofL. maculatus.

Comparative Mitogenomic Analysis of Forensically Important Sarcophagid Flies (Diptera: Sarcophagidae) and Implications of Species Identification

The flesh flies (Diptera: Sarcophagidae) are significant in forensic investigations. The mitochondrial genome (mitogeome) has been widely used as genetic markers for phylogenetic analysis and species identification. To further understand the mitogenome-level features in Sarcophagidae, the complete mitogenome of Sarcophaga formosensis (Kirneret Lopes, 1961) (Diptera: Sarcophagidae) and Sarcophaga misera (Walker, 1849) (Diptera: Sarcophagidae) was firstly sequenced, annotated, and compared with other 13 Sarcophagidae species. The result indicated that the gene arrangement, gene content, base composition, and codon usage were conserved in the ancestral arthropod. Evolutionary rate of the mitogenome fragments revealed that the nonsynonymous and synonymous substitution rates (Ka and Ks) ratio was less than 1.00, indicating these variable sites under strong purifying selection. Almost all transfer RNA genes (tRNAs) have typical clover-leaf structures within these sarcophagid mitogenomes, except tRNA-Ser (AGN) is lack of the dihydrouridine arm. This comparative mitogenomic analysis sheds light on the architecture and evolution of mitogenomes in the Sarcophagidae. Phylogenetic analyses containing the interspecific distances from different regions in these species provided us new insights into the application of these effective genetic markers for species identification of flesh flies.

The complete mitochondrial genome of Melon thrips, Thrips palmi (Thripinae): Comparative analysis

The melon thrips, Thrips palmi is a serious pest and vector for plant viruses on a wide range of economically important crops. DNA barcoding evidenced the presence of cryptic diversity in T. palmi and that warrants exhaustive molecular studies. Our present study is on decoding the first complete mitochondrial genome of T. palmi (15,333 bp) through next-generation sequencing (NGS). The T. palmi mt genome contains 37 genes, including 13 Protein coding genes (PCGs), two ribosomal RNA (rRNAs), 22 transfer RNA (tRNAs), and two control regions (CRs). The majority strand of T. palmi revealed 78.29% A+T content, and 21.72% G+C content with positive AT skew (0.09) and negative GC skew (-0.06). The ATN initiation codons were observed in 12 PCGs except for cox1 which have unique start codon (TTG). The relative synonymous codon usage (RSCU) analysis revealed Phe, Leu, Ile, Tyr, Asn, Lys and Met were the most frequently used amino acids in all PCGs. The codon (CGG) which is assigned to Arginine in most insects but absent in T. palmi. The Ka/Ks ratio ranges from 0.078 in cox1 to 0.913 in atp8. We observed the typical cloverleaf secondary structure in most of the tRNA genes with a few exceptions; absence of DHU stem and loop in trnV and trnS, absence of DHU loop in trnE, lack of T-arm and loop in trnN. The T. palmi gene order (GO) was compared with ancestral GO and observed an extensive gene arrangement in PCGs, tRNAs and rRNAs. The cox2 gene was separated from the gene block 'cox2-trnL2' in T. palmi as compared with the other thrips mt genomes, including ancestor GO. Further, the nad1, trnQ, trnC, trnL1, trnV, trnF, rrnS, and rrnL were inversely transpositioned in T. palmi GO. The gene blocks 'trnQ-trnS2-trnD' and 'trnN-trnE-trnS1-trnL1' seems to be genus specific. The T. palmi mt genome contained 24 intergenic spacer regions and 12 overlapping regions. The 62 bp of CR2 shows the similarity with CR1 indicating a possible duplication. The occurrence of multiple CRs in thrips mt genomes seems to be a derived trait which needs further investigation. Although, the study depicted extensive gene rearrangements in T. palmi mt genome, but the negative GC skew reflects only strand asymmetry. Both the ML and BI phylogenetic trees revealed the close relationships of Thrips with Scirtothrips as compared to Frankliniella. Thus, more mt genomes of the diverse thrips species are required to understand the in-depth phylogenetic and evolutionary relationships.

The complete mitochondrial genome of Somanniathelphusa boyangensis and phylogenetic analysis of Genus Somanniathelphusa (Crustacea: Decapoda: Parathelphusidae)

In this study, the authors first obtained the mitochondrial genome of Somanniathelphusa boyangensis. The results showed that the mitochondrial genome is 17,032bp in length, included 13 protein-coding genes, 2 rRNAs genes, 22 tRNAs genes and 1 putative control region, and it has the characteristics of the metazoan mitochondrial genome A+T bias. All tRNA genes display the typical clover-leaf secondary structure except tRNASer(AGN), which has lost the dihydroxyuridine arm. The GenBank database contains the mitochondrial genomes of representatives of approximately 22 families of Brachyura, comprising 56 species, including 4 species of freshwater crab. The authors established the phylogenetic relationships using the maximum likelihood and Bayesian inference methods. The phylogenetic relationship indicated that the molecular taxonomy of S. boyangensis is consistent with current morphological classification, and Parathelphusidae and Potamidae are derived within the freshwater clade or as part of it. In addition, the authors used the COX1 sequence of Somanniathelphusa in GenBank and the COX1 sequence of S. boyangensis to estimated the divergence time of this genus. The result displayed that the divergence time of Somanniathelphusa qiongshanensis is consistent with the separation of Hainan Island from mainland China in the Beibu Gulf, and the divergence time for Somanniathelphusa taiwanensis and Somanniathelphusa amoyensis is consistent with the separation of Taiwan Province from Mainland China at Fujian Province. These data indicate that geologic events influenced speciation of the genus Somanniathelphusa.

Sequencing and characterization of the complete mitochondrial genome of Japanese Swellshark (Cephalloscyllium umbratile)

To further comprehend the genome features of Cephalloscyllium umbratile (Carcharhiniformes), an endangered species, the complete mitochondrial DNA (mtDNA) was firstly sequenced and annotated. The full-length mtDNA of C. umbratile was 16,697 bp and contained ribosomal RNA (rRNA) genes, 13 protein-coding genes (PCGs), 23 transfer RNA (tRNA) genes, and a major non-coding control region. Each PCG was initiated by an authoritative ATN codon, except for COX1 initiated by a GTG codon. Seven of 13 PCGs had a typical TAA termination codon, while others terminated with a single T or TA. Moreover, the relative synonymous codon usage of the 13 PCGs was consistent with that of other published Carcharhiniformes. All tRNA genes had typical clover-leaf secondary structures, except for tRNA-Ser (GCT), which lacked the dihydrouridine 'DHU' arm. Furthermore, the analysis of the average Ka/Ks in the 13 PCGs of three Carcharhiniformes species indicated a strong purifying selection within this group. In addition, phylogenetic analysis revealed that C. umbratile was closely related to Glyphis glyphis and Glyphis garricki. Our data supply a useful resource for further studies on genetic diversity and population structure of C. umbratile.

Complete mitochondrial genome of Clistocoeloma sinensis (Brachyura: Grapsoidea): Gene rearrangements and higher-level phylogeny of the Brachyura

Deciphering the animal mitochondrial genome (mitogenome) is very important to understand their molecular evolution and phylogenetic relationships. In this study, the complete mitogenome of Clistocoeloma sinensis was determined. The mitogenome of C. sinensis was 15,706 bp long, and its A+T content was 75.7%. The A+T skew of the mitogenome of C. sinensis was slightly negative (−0.020). All the transfer RNA genes had the typical cloverleaf structure, except for the trnS1 gene, which lacked a dihydroxyuridine arm. The two ribosomal RNA genes had 80.2% A+T content. The A+T-rich region spanned 684 bp. The gene order within the complete mitogenome of C. sinensis was identical to the pancrustacean ground pattern except for the translocation of trnH. Additionally, the gene order of trnI-trnQ-trnM in the pancrustacean ground pattern becomes trnQ-trnI-trnM in C. sinensis. Our phylogenetic analysis showed that C. sinensis and Sesarmops sinensis cluster together with high nodal support values, indicating that C. sinensis and S. sinensis have a sister group relationship. The results support that C. sinensis belongs to Grapsoidea, Sesarmidae. Our findings also indicate that Varunidae and Sesarmidae species share close relationships. Thus, mitogenomes are likely to be valuable tools for systematics in other groups of Crustacea.

The complete mitochondrial genome of Sesarmops sinensis reveals gene rearrangements and phylogenetic relationships in Brachyura

Mitochondrial genome (mitogenome) is very important to understand molecular evolution and phylogenetics. Herein, in this study, the complete mitogenome of Sesarmops sinensis was reported. The mitogenome was 15,905 bp in size, and contained 13 protein-coding genes (PCGs), two ribosomal RNA (rRNA) genes, 22 transfer RNA (tRNA) genes, and a control region (CR). The AT skew and the GC skew are both negative in the mitogenomes of S. sinensis. The nucleotide composition of the S. sinensis mitogenome was also biased toward A + T nucleotides (75.7%). All tRNA genes displayed a typical mitochondrial tRNA cloverleaf structure, except for the trnS1 gene, which lacked a dihydroxyuridine arm. S. sinensis exhibits a novel rearrangement compared with the Pancrustacean ground pattern and other Brachyura species. Based on the 13 PCGs, the phylogenetic analysis showed that S. sinensis and Sesarma neglectum were clustered on one branch with high nodal support values, indicating that S. sinensis and S. neglectum have a sister group relationship. The group (S. sinensis + S. neglectum) was sister to (Parasesarmops tripectinis + Metopaulias depressus), suggesting that S. sinensis belongs to Grapsoidea, Sesarmidae. Phylogenetic trees based on amino acid sequences and nucleotide sequences of mitochondrial 13 PCGs using BI and ML respectively indicate that section Eubrachyura consists of four groups clearly. The resulting phylogeny supports the establishment of a separate subsection Potamoida. These four groups correspond to four subsections of Raninoida, Heterotremata, Potamoida, and Thoracotremata.

Digging deeper: new gene order rearrangements and distinct patterns of codons usage in mitochondrial genomes among shrimps from the Axiidea, Gebiidea and Caridea (Crustacea: Decapoda)

BACKGROUND

Whole mitochondrial DNA is being increasingly utilized for comparative genomic and phylogenetic studies at deep and shallow evolutionary levels for a range of taxonomic groups. Although mitogenome sequences are deposited at an increasing rate into public databases, their taxonomic representation is unequal across major taxonomic groups. In the case of decapod crustaceans, several infraorders, including Axiidea (ghost shrimps, sponge shrimps, and mud lobsters) and Caridea (true shrimps) are still under-represented, limiting comprehensive phylogenetic studies that utilize mitogenomic information.

METHODS

Sequence reads from partial genome scans were generated using the Illumina MiSeq platform and mitogenome sequences were assembled from these low coverage reads. In addition to examining phylogenetic relationships within the three infraorders, Axiidea, Gebiidea, and Caridea, we also investigated the diversity and frequency of codon usage bias and mitogenome gene order rearrangements.

RESULTS

We present new mitogenome sequences for five shrimp species from Australia that includes two ghost shrimps, Callianassa ceramica and Trypaea australiensis, along with three caridean shrimps, Macrobrachium bullatum, Alpheus lobidens, and Caridina cf. nilotica. Strong differences in codon usage were discovered among the three infraorders and significant gene order rearrangements were observed. While the gene order rearrangements are congruent with the inferred phylogenetic relationships and consistent with taxonomic classification, they are unevenly distributed within and among the three infraorders.

DISCUSSION

Our findings suggest potential for mitogenome rearrangements to be useful phylogenetic markers for decapod crustaceans and at the same time raise important questions concerning the drivers of mitogenome evolution in different decapod crustacean lineages.

Characterization of complete mitochondrial genome of fives tripe wrasse (Thalassoma quinquevittatum, Lay & Bennett, 1839) and phylogenetic analysis

To further supplement the genome-level features in related species, T. quinquevittatum complete mtDNA was firstly sequenced and de novo assembled by next-generation sequencing. The full-length mtDNA of T. quinquevittatum was a 16,896bp fragment, which was atypical of Labridae, with 2 ribosomal RNA (rRNA) genes, 13 protein-coding genes (PCGs), 23 transfer RNA (tRNA) genes, and a major non-coding control region (D-loop region). Additionally, the mtDNA of T. quinquevittatum exhibited characteristics of A (27.1%), T (29.3%), G (17.8%), and C (25.8%) with a high A+T content (56.4%). Furthermore, the analysis of the average Ka/Ks in the 13 PCGs of three Labridae species indicated a strong purifying selection within this group. Additionally, the phylogenetic analysis based on 13 concatenated PCGs nucleotide and amino acid datasets, showed high value support for the following sister clade among the four genera (T. quinquevittatum, Halichoeres trimaculatus, Halichoeres melanurus, Parajulis poecilepterus). The complete mtDNA of the T. quinquevittatum provided important information for the study in population genetics and evolutionary theory.

Complete mitochondrial genome of Rhynchocinetes durbanensis (Rhynchocinetidae: Rhynchocinetes)

The complete mitochondrial genome of Rhynchocinetes durbanensis (Rhynchocinetidae: Rhynchocinetes) was sequenced in this study. The genome sequence was 17,695 bp in size, with the base composition of 35.14% A, 32.98% T, 20.34% G and 11.55% C of the light strand. The gene order and genes were the same as that found in other previously reported shrimps, including 13 protein-coding genes, 24 transfer RNA genes and two ribosomal RNA genes. Except for ND5, ND4, ND4L, ND1 genes and eight tRNA genes and two ribosomal RNA genes, all other mitochondrial genes were encoded on the heavy strand. The start codon of COX1 gene was not determined. These complete mitogenome data provide the basis for taxonomic and conservation research of Rhynchocinetes durbanensis, and closely related species.

Identification of the complete mitochondrial genome of the pacific white shrimp Litopenaeus vannamei (Decapoda: Penaeidae)

The complete mitochondrial genome (mitogenome) of Litopenaeus vannamei was determined to be 15 989 bp in size, it contains 22 tRNA genes, 13 protein-coding genes (PCGs), two rRNA genes and a D-loop region. The nucleotide composition was biased toward A + T nucleotides (67.8%) and the AT skew of this mitogenome was slightly negative (-0.025). All tRNA genes displayed a typical clover-leaf structure, except for trnS1 (AGN). Thirteen PCGs were initiated by ATN codons, except for cox1 gene which was initiated by ACG. Four of the 13 PCGs had the incomplete termination codon by T. The rrnL was 1371 bp and the rrnS was 853 bp, and the AT content of the combined rRNA genes was 70.9%. The D-loop region of the mitogenome was 985 bp in length and the A + T content was 82.7%, and no tandem repeat was found in this region. Phylogenetic analysis showed that the placement of L. vannamei was within the Penaeidae.

Shotgun assembly of the mitochondrial genome from Fenneropenaeus penicillatus with phylogenetic consideration

The complete mitochondrial genome is of great importance for better understanding of the genome-level characteristics and phylogenetic relationships among related species. In this study, Fenneropenaeus penicillatus mitochondrial genome sequence was determined by next-generation sequencing. The complete genome DNA was 16,040 bp in length and consisted of a typical set of 13 protein-coding genes, 22 transfer RNA (tRNA) genes, two ribosomal RNA (rRNA) genes and a putative control region (CR). The gene arrangement is identical to the pancrustacean pattern. The overall base composition of its mitochondrial genome is estimated to be 34.1% for A, 34.1% for T, 12.5% for G and 19.3% for C with a high A+T content (68.2%). The analysis of the average Ka/Ks in the 13 mitochondrial protein-coding genes of penaeid shrimps indicated a strong purifying selection within this group. The phylogenetic analysis based on mitochondrial sequences and 13 concatenated protein-coding genes showed strong statistic support for the following relationship among the five genera ((Penaeus s.s+Fenneropenaeus)+(Litopenaeus+Farfantepenaeus))+Marsupenaeus. The sequence data of F. penicillatus can provide useful information for the studies on molecular systematics, population structure, stock evaluation and conservation genetics.

First mitochondrial genome for the red crab (Charybdis feriata) with implication of phylogenomics and population genetics

In this study, we first described the complete mitochondrial genome for the red crab (Charybdis feriata), elucidated its phylogenetic relationship among 20 species within Decapoda, and estimated the population genetic diversity. The mitochondrial genome was 15,660 bp in size and encoded 13 protein-coding genes, 22 transfer RNA (tRNA) genes, and two ribosomal RNA genes. The gene arrangement of the mitochondrial genome was the same as that of its sister species, C. japonica. Phylogenomic analysis suggested that genus Charybdis should be classified into subfamily Portuninae but not into subfamily Thalamitinae. Moreover, a total of 33 haplotypes of complete cytochrome c oxidase subunit I gene were defined in 70 individuals of C. feriata derived from three localities. Haplotype diversity and nucleotide diversity values among three localities indicated a high level of genetic diversity in C. feriata. AMOVA analysis suggested a low level of genetic differentiation among the three localities (FST = 0.0023, P > 0.05). Neutrality tests and mismatch analysis revealed that C. feriata might have undergone a population expansion event that possibly occurred in the last 61,498 to 43,814 years. This study should be helpful to better understand the evolutionary status, and population genetic diversity of C. feriata and related species.

Mitochondrial genome of the intertidal acorn barnacle Tetraclita serrata Darwin, 1854 (Crustacea: Sessilia): Gene order comparison and phylogenetic consideration within Sessilia

The complete mitochondrial genome of the intertidal barnacle Tetraclita serrata Darwin, 1854 (Crustacea: Maxillopoda: Sessilia) is presented. The genome is a circular molecule of 15,200 bp, which encodes 13 PCGs, 2 ribosomal RNA genes, and 22 transfer RNA genes. All non-coding regions are 591 bp in length, with the longest one speculated as the control region (389 bp), which is located between srRNA and trnK. The overall A+T content of the mitochondrial genome of T. serrata is 65.4%, which is lowest among all the eight mitochondrial genomes reported from sessile barnacles. There are variations of initiation and stop codons in the reported sessile barnacle mitochondrial genomes. Large-scale gene rearrangements are found in these genomes as compared to the pancrustacean ground pattern. ML and Bayesian analyses of all 15 complete mitochondrial genomes available from Maxillopoda lead to identical phylogenies. The phylogenetic tree based on mitochondrial PCGs shows that Argulus americanus (Branchiura) cluster with Armillifer armillatus (Pentastomida), distinct from all ten species from Cirripedia. Within the order Sessilia, Amphibalanus amphitrite (Balanidae) clusters with Striatobalanus amaryllis (Archaeobalanidae), and Nobia grandis (Pyrgomatidae). However, the two Megabalanus (Balanidae) are separated from the above grouping, resulting in non-monophyly of the family Balanidae. Moreover, the two Megabalanus have large-scale rearrangements as compared to the gene order shared by former three species. Therefore, both phylogenetic analysis using PCG sequences and gene order comparison suggest that Balanidae is not a monophyletic group. Given the limited taxa and moderate support values of the internal branches, the non-monophyly of the family Balanidae requires further verification.

The mitochondrial genome of Nobia grandis Sowerby, 1839 (Cirripedia: Sessilia): the first report from the coral-inhabiting barnacles family Pyrgomatidae

This work presents the coral-inhabiting barnacle Nobia grandis Sowerby, 1839 complete mitochondrial genome, which is the first report from the family Pyrgomatidae (Cirripedia: Sessilia). The N. grandis mitochondrial genome is 15,032 bp in length, containing a total of 469 bp of non-coding nucleotides spreading in 11 intergenic regions (with the largest region of 376 bp). Compared with the pancrustacean ground pattern, there are not less than seven tRNAs rearranged in the N. grandis mitochondrial genome. Gene overlaps are founded in eight places. Nine PCGs (COX1-3, ATP6, ATP8, CYTB, ND2, ND3 and ND6) are encoded on the heavy strand while the remaining 4 PCGs and the two rRNAs are located on the light strand. As the first representative from the family Pyrgomatidae, the N. grandis mitochondrial genome will help us to explore the evolutionary history and molecular evolution of coral barnacles and Sessilia in future studies.

Complete mitochondrial genome of the acorn barnacle Striatobalanus amaryllis (Crustacea: Maxillopoda): the first representative from Archaeobalanidae

The mitochondrial genome of the barnacle Striatobalanus amaryllis (Sessilia: family Archaeobalanidae) is 15,063 bp in length. All the 13 protein-coding genes (PCGs) initiate with ATD codon (ATG, ATA or ATT). Four PCGs (COX3, ND3, ND4 and ND4L) end with incomplete stop codon (T- -). Four PCGs (ND1, ND4, ND4L and ND5) are encoded on the light strand (underlined below). Refer to the pancrustacean ground pattern, there are not less than seven tRNAs rearranged in the S. amaryllis mitochondrial genome, including tRNA(Ala), tRNA(Glu)/tRNA(Ser)((AGY)), tRNA(Pro)/tRNA(Thr), tRNA(Pro)/tRNA(Thr), tRNA(Tyr), tRNA(Lys), tRNA(Gln) and tRNA(Cys). Three tRNAs (tRNA(Lys), tRNA(Gln) and tRNA(Cys)) are rearranged between S. amaryllis and Tetraclita japonica (Sessilia: Tetraclitidae), meanwhile one tRNA (tRNA(Cys)) inverted from one strand to another. Compared with Megabalanus volcano (Sessilia: Balanidae), an inversion of one large gene block is identified (including three PCGs and three tRNAs) in S. amaryllis mitochondrial genome: tRNA(Phe)-ND5-tRNA(His)-ND4-ND4L-tRNA(Pro).

RNA-Seq analysis reveals genes associated with resistance to Taura syndrome virus (TSV) in the Pacific white shrimp Litopenaeus vannamei

Outbreak of Taura syndrome virus (TSV) is one of the major pathogens of the Pacific white shrimp Litopenaeus vannamei. Although selective breeding for improvement of TSV resistance in L. vannamei has been successfully developed and has led to a great benefit to the shrimp farming industry worldwide. The molecular mechanisms underlying the viral resistance in shrimp remain largely unknown. In the present study, we conducted the first transcriptomic profiling of host responses in hemolymph and hemocytes in order to identify the differentially expressed genes associated with resistance to TSV in L. vannamei. High-throughput RNA-Seq was employed, obtaining 193.6 and 171.2 million high-quality Illumina reads from TSV-resistant and susceptible L. vannamei lines respectively. A total of 61,937 contigs were generated with an average length of 546.26 bp. BLASTX-based gene annotation (E-value < 10(-5)) allowed the identification of 12,398 unique proteins against the NCBI non-redundant NR database. In addition, comparison of digital gene expression between resistant and susceptible strains revealed 1374 significantly differentially expressed contigs (representing 697 unigenes). Gene pathway analysis of the differentially expressed gene set highlighted several putative genes involved in the immune response activity including (1) pathogen/antigen recognition including immune regulator, adhesive protein and signal transducer; (2) coagulation; (3) proPO pathway cascade; (4) antioxidation; and (5) protease. The expression patterns of 22 differentially expressed genes involving immune response were validated by quantitative real-time RT-PCR (average correlation coefficients 0.94, p-value < 0.001). Our results provide valuable information on gene functions associated with resistance to TSV in L. vannamei.

The first mitochondrial genome from Mysida (Crustacea: Malacostraca) reveals an unusual gene arrangement

This is the first report to present the Neomysis orientalis mitochondrial genome as a representative from the order Mysida. While mitochondrial protein-coding genes (PCGs) commonly use several alternatives to ATN as start codons, all 13 PCGs in N. orientalis mitochondrial genome initiate with ATG or ATA. Five PCGs (atp6. atp8. cob. nad4 and nad4L) start with ATG, while the other genes (cox1-3. nad1-3. nad5 and nad6) start with ATA. Only two PCGs (cox2 and nad2) in the N. orientalis mitochondrial genome end with incomplete stop codons (T- or TA-), and all the remaining ones have TAA or TAG stop codon. Only one PCG (nad4L) is encoded on the light strand and all other 12 PCGs are located at the heavy strand. Both rRNAs (srRNA and lrRNA) are encoded on the light strand. In common with 15 of the other 18 mitochondrial genomes from Peracarida, the major gene arrangement in the N. orientalis mitochondrial genome is different from the pancrustacean ground pattern. The largest conserved gene block in N. orientalis only contains two genes but those in the other 18 peracarid mitochondrial genomes contain more than five genes. Thus, the N. orientalis mitochondrial genome, as the first mitochondrial genome from the order Mysida, reveals an unusual gene arrangement that is unique compared with the other malacostracan mitochondrial genomes.

Three nucleus-encoded subunits of mitochondrial cytochrome c oxidase of the whiteleg shrimp Litopenaeus vannamei: cDNA characterization, phylogeny and mRNA expression during hypoxia and reoxygenation

The mitochondrial cytochrome c oxidase (COX) catalyzes the reduction of oxygen to water playing a key role in the respiratory chain and ATP synthesis. The nucleus-encoded COX subunits do not participate in catalysis, but some are known to play a role in the expression, assembly and activity of the enzyme. Since hypoxia continuously affects the shrimp environment, it is important to study COX to understand their ability to deal with low oxygen levels. The goal of this research was to characterize the complementary DNA (cDNA) sequences of three nucleus-encoded subunits -coxIV, coxVa, and coxVb- and to evaluate the shrimp COX response to hypoxia by measuring their gene expression. The cDNA sequence of coxIV consisted of 532bp, which encodes a 17.47kDa protein, while coxVa cDNA consisted of 460bp and coded a protein of 17.11kDa, and the coxVb coding sequence consisted of 364bp encoding a 13.74kDa protein. Shrimp subunits do not have isoforms, and they are not differentially expressed during hypoxia, as observed in mammals. Coordinated changes were detected in the mRNA amounts of nuclear and mitochondrial subnits; these changes, at the transcriptional level, are suggested to be controlled through transcriptional factors Sp1 and NRF2.

The complete mitochondrial genome sequence and gene organization of the mud crab (Scylla paramamosain) with phylogenetic consideration

The complete mitochondrial genome is of great importance for better understanding the genome-level characteristics and phylogenetic relationships among related species. In the present study, we determined the complete mitochondrial genome DNA sequence of the mud crab (Scylla paramamosain) by 454 deep sequencing and Sanger sequencing approaches. The complete genome DNA was 15,824 bp in length and contained a typical set of 13 protein-coding genes, 22 transfer RNA (tRNA) genes, two ribosomal RNA (rRNA) genes and a putative control region (CR). Of 37 genes, twenty-three were encoded by the heavy strand (H-strand), while the other ones were encoded by light strand (L-strand). The gene order in the mitochondrial genome was largely identical to those obtained in most arthropods, although the relative position of gene tRNA(His) differed from other arthropods. Among 13 protein-coding genes, three (ATPase subunit 6 (ATP6), NADH dehydrogenase subunits 1 (ND1) and ND3) started with a rare start codon ATT, whereas, one gene cytochrome c oxidase subunit I (COI) ended with the incomplete stop codon TA. All 22 tRNAs could fold into a typical clover-leaf secondary structure, with the gene sizes ranging from 63 to 73 bp. The phylogenetic analysis based on 12 concatenated protein-coding genes showed that the molecular genetic relationship of 19 species of 11 genera was identical to the traditional taxonomy.

Analysis of Litopenaeus vannamei transcriptome using the next-generation DNA sequencing technique

BACKGROUND

Pacific white shrimp (Litopenaeus vannamei), the major species of farmed shrimps in the world, has been attracting extensive studies, which require more and more genome background knowledge. The now available transcriptome data of L. vannamei are insufficient for research requirements, and have not been adequately assembled and annotated.

METHODOLOGY/PRINCIPAL FINDINGS

This is the first study that used a next-generation high-throughput DNA sequencing technique, the Solexa/Illumina GA II method, to analyze the transcriptome from whole bodies of L. vannamei larvae. More than 2.4 Gb of raw data were generated, and 109,169 unigenes with a mean length of 396 bp were assembled using the SOAP denovo software. 73,505 unigenes (>200 bp) with good quality sequences were selected and subjected to annotation analysis, among which 37.80% can be matched in NCBI Nr database, 37.3% matched in Swissprot, and 44.1% matched in TrEMBL. Using BLAST and BLAST2Go softwares, 11,153 unigenes were classified into 25 Clusters of Orthologous Groups of proteins (COG) categories, 8171 unigenes were assigned into 51 Gene ontology (GO) functional groups, and 18,154 unigenes were divided into 220 Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways. To primarily verify part of the results of assembly and annotations, 12 assembled unigenes that are homologous to many embryo development-related genes were chosen and subjected to RT-PCR for electrophoresis and Sanger sequencing analyses, and to real-time PCR for expression profile analyses during embryo development.

CONCLUSIONS/SIGNIFICANCE

The L. vannamei transcriptome analyzed using the next-generation sequencing technique enriches the information of L. vannamei genes, which will facilitate our understanding of the genome background of crustaceans, and promote the studies on L. vannamei.

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 mitochondrial genome of the Japanese snapping shrimp Alpheus japonicus (Crustacea: Decapoda: Caridea): gene rearrangement and phylogeny within Caridea

The complete sequence of the mitochondrial genome of the Japanese snapping shrimp Alpheus japonicus Miers (Crustacea: Decapoda: Caridea) is presented here. A comparative analysis based on the currently available mitochondrial genomic data revealed many previously unknown characteristics of the mitochondrial genomes of caridean shrimps. The A. japonicus mitochondrial genome is 16487 bp long and contains the typical set of 37 metazoan genes. The gene arrangements in the mitochondrial genomes of four previously studied carideans (Macrobrachium rosenbergii, M. nipponense, M. lanchesteri and Halocaridina rubra) were found to be identical to the pancrustacean ground pattern; thus, it was considered that gene rearrangements probably did not occur in the suborder Caridea. In the present study, a translocation of the trnE gene involving inversion was found in Alpheus mitochondrial genomes. This phenomenon has not been reported in any other crustacean mitochondrial genome that has been studied so far; however, the translocation of one transfer RNA gene (trnP or trnT) was reported in the mitochondrial genome of Exopalaemon carinicauda. When the ratios of the nonsynonymous and synonymous substitutions rates (Ka/Ks) for the 13 protein coding genes from two Alpheus species (A. japonicus and A. distinguendus) and three Macrobrachium species (M. rosenbergii, M. nipponense, M. lanchesteri) were calculated, the Ka/Ks values for all the protein coding genes in Alpheus and Macrobrachium mitochondrial genomes were found to be less than 1 (between 0.0048 and 0.2057), indicating that a strong purification selection had occurred. The phylogenetic tree that was constructed based on the mitochondrial protein coding genes in the genomes of nine related species indicated that Palaemonidae and Alpheidae formed a monophyly and shared a statistically significant relationship, (Palaemonidae+Alpheidae)+Atyidae, at the family level.

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.

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 mantid shrimp Oratosquilla oratoria (Crustacea: Malacostraca: Stomatopoda): Novel non-coding regions features and phylogenetic implications of the Stomatopoda

The complete mitochondrial (mt) genome sequence of Oratosquilla oratoria (Crustacea: Malacostraca: Stomatopoda) was determined; a circular molecule of 15,783 bp in length. The gene content and arrangement are consistent with the pancrustacean ground pattern. The mt control region of O. oratoria is characterized by no GA-block near the 3' end and different position of [TA(A)]n-blocks compared with other reported Stomatopoda species. The sequence of the second hairpin structure is relative conserved which suggests this region may be a synapomorphic character for the Stomatopoda. In addition, a relative large intergenic spacer (101 bp) with higher A+T content than that in control region was identified between the tRNA(Glu) and tRNA(Phe) genes. Phylogenetic analyses based on the current dataset of complete mt genomes strongly support the Stomatopoda is closely related to Euphausiacea. They in turn cluster with Penaeoidea and Caridea clades while other decapods form a separate group, which rejects the monophyly of Decapoda. This challenges the suitability of Stomatopoda as an outgroup of Decapoda in phylogenetic analyses. The basal position of Stomatopoda within Eumalacostraca according to the morphological characters is also questioned.

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

The genetics and molecular biology of the commercially important Chinese spiny lobster, Panulirus stimpsoni are little known. Here, we present the complete mitochondrial genome sequence of P. stimpsoni, determined by the long polymerase chain reaction and primer walking sequencing method. The entire genome is 15,677 bp in length, encoding the standard set of 13 protein-coding genes, two ribosomal RNA genes, and 22 transfer RNA genes. The overall A+T content of the genome is 65.6%, lower than most malacostracan species. The gene order is consistent with the pancrustacean ground pattern. Several conserved elements were identified from P. stimpsoni control region, viz. one [TA(A)]n-block, two GA-blocks and three hairpin structures. However, the position of [TA(A)]n-block and number of hairpin structure are different from those in the congeneric P. japonicus and other decapods. Phylogenetic analyses using the concatenated nucleotide and amino acid sequences of 13 protein-coding genes do not support the monophyly of suborder Pleocyemata, which is in contrast to most morphological and molecular results. However, the position of Palinura and Astacidea is unstable, as represented by the basal or sister branches to other Reptantia species. P. stimpsoni, as the second species of Palinura with complete mitochondrial genome available, will provide important information on both genomics and conservation biology of the group.

Complete mitochondrial genome of the Asian paddle crab Charybdis japonica (Crustacea: Decapoda: Portunidae): gene rearrangement of the marine brachyurans and phylogenetic considerations of the decapods

Given the commercial and ecological importance of the Asian paddle crab, Charybdis japonica, there is a clearly need for genetic and molecular research on this species. Here, we present the complete mitochondrial genome sequence of C. japonica, determined by the long-polymerase chain reaction and primer walking sequencing method. The entire genome is 15,738 bp in length, encoding a standard set of 13 protein-coding genes, two ribosomal RNA genes, and 22 transfer RNA genes, plus the putative control region, which is typical for metazoans. The total A+T content of the genome is 69.2%, lower than the other brachyuran crabs except for Callinectes sapidus. The gene order is identical to the published marine brachyurans and differs from the ancestral pancrustacean order by only the position of the tRNA ( His ) gene. Phylogenetic analyses using the concatenated nucleotide and amino acid sequences of 13 protein-coding genes strongly support the monophyly of Dendrobranchiata and Pleocyemata, which is consistent with the previous taxonomic classification. However, the systematic status of Charybdis within subfamily Thalamitinae of family Portunidae is not supported. C. japonica, as the first species of Charybdis with complete mitochondrial genome available, will provide important information on both genomics and molecular ecology of the group.

The complete mitogenome of the hydrothermal vent crab Xenograpsus testudinatus (Decapoda, Brachyura) and comparison with brachyuran crabs

In this study, we analyzed the complete mitochondrial (mt) genome of a hydrothermal vent crab Xenograpsus testudinatus (Decapoda: Brachyura) obtained from the hydrothermal vents off Kueishantao Island, Taiwan, which extend from the deep sea Okinawa Trench. The mitogenome of X. testudinatus was 15,796 bp in length and contained the same 37 genes (e.g. 2 rRNAs, 22 tRNAs, and 13 PCGs) found in other metazoan mitogenomes. Analysis of the structural mt gene order in X. testudinatus revealed that the 13 PCGs, excluding a translocation of ND6-Cyt b cluster, were similarly ordered when compared to the pancrustacean ground pattern; however the tRNAs were severely rearranged. Phylogenetic analysis of decapod mitogenomes showed that the molecular taxonomy of the vent crab was in accordance with its morphological systematics. Together, these findings suggest that the vent crab studied here has little mitochondrial genetic variation when compared with morphologically defined conspecifics from other marine habitats.

A close phylogenetic relationship between Sipuncula and Annelida evidenced from the complete mitochondrial genome sequence of Phascolosoma esculenta

Background

There are many advantages to the application of complete mitochondrial (mt) genomes in the accurate reconstruction of phylogenetic relationships in Metazoa. Although over one thousand metazoan genomes have been sequenced, the taxonomic sampling is highly biased, left with many phyla without a single representative of complete mitochondrial genome. Sipuncula (peanut worms or star worms) is a small taxon of worm-like marine organisms with an uncertain phylogenetic position. In this report, we present the mitochondrial genome sequence of Phascolosoma esculenta, the first complete mitochondrial genome of the phylum.

Results

The mitochondrial genome of P.esculenta is 15,494 bp in length. The coding strand consists of 32.1% A, 21.5% C, 13.0% G, and 33.4% T bases (AT = 65.5%; AT skew = -0.019; GC skew = -0.248). It contains thirteen protein-coding genes (PCGs) with 3,709 codons in total, twenty-two transfer RNA genes, two ribosomal RNA genes and a non-coding AT-rich region (AT = 74.2%). All of the 37 identified genes are transcribed from the same DNA strand. Compared with the typical set of metazoan mt genomes, sipunculid lacks trnR but has an additional trnM. Maximum Likelihood and Bayesian analyses of the protein sequences show that Myzostomida, Sipuncula and Annelida (including echiurans and pogonophorans) form a monophyletic group, which supports a closer relationship between Sipuncula and Annelida than with Mollusca, Brachiopoda, and some other lophotrochozoan groups.

Conclusion

This is the first report of a complete mitochondrial genome as a representative within the phylum Sipuncula. It shares many more similar features with the four known annelid and one echiuran mtDNAs. Firstly, sipunculans and annelids share quite similar gene order in the mitochondrial genome, with all 37 genes located on the same strand; secondly, phylogenetic analyses based on the concatenated protein sequences also strongly support the sipunculan + annelid clade (including echiurans and pogonophorans). Hence annelid "key-characters" including segmentation may be more labile than previously assumed.

The complete mitochondrial genome of the ridgetail white prawn Exopalaemon carinicauda Holthuis, 1950 (Crustacean: Decapoda: Palaemonidae) revealed a novel rearrangement of tRNA genes

The complete mitochondrial (mt) DNA sequence was determined for a ridgetail white prawn, Exopalaemon carinicauda Holthuis, 1950 (Crustacea: Decopoda: Palaemonidae). The mt genome is 15,730 bp in length, encoding a standard set of 13 protein-coding genes, 2 ribosomal RNA genes, and 22 transfer RNA genes, which is typical for metazoans. The majority-strand consists of 33.6% A, 23.0% C, 13.4% G, and 30.0% T bases (AT skew=0.057; GC skew=-0.264). A total of 1045 bp of non-coding nucleotides were observed in 16 intergenic regions, including a major A+T rich (79.7%) noncoding region (886 bp). A novel translocation of tRNA(Pro) and tRNA(Thr) was found when comparing this genome with the pancrustacean ground pattern indicating that gene order is not conserved among caridean mitochondria. Furthermore, the rate of Ka/Ks in 13 protein-coding genes between three caridean species is much less than 1, which indicates a strong purifying selection within this group. To investigate the phylogenetic relationship within Malacostraca, phylogenetic trees based on currently available malacostracan complete mitochondrial sequences were built with the maximum likelihood and Bayesian models. All analyses based on nucleotide and amino acid data strongly support the monophyly of Decapoda. The Penaeidae, Reptantia, Caridea, and Meiura clades were also recovered as monophyletic groups with strong statistical support. However, the phylogenetic relationships within Pleocyemata are unstable, as represented by the inclusion or exclusion of Caridea.

The complete mitochondrial genomes of the yellowleg shrimp Farfantepenaeus californiensis and the blue shrimp Litopenaeus stylirostris (Crustacea: Decapoda)

Mitochondria play key roles in many cellular processes. Description of penaeid shrimp genes, including mitochondrial genomes are fairly recent and some are focusing on commercially important shrimp as the Pacific shrimp Litopenaeus vannamei that is being used for aquaculture not only in America, but also in Asia. Much less is known about other Pacific shrimp such as the yellowleg shrimp Farfantepenaeus californiensis and the blue shrimp Litopenaeus stylirostris. We report the complete mitogenomes from these last two Pacific shrimp species. Long DNA fragments were obtained by PCR and then used to get internal fragments for sequencing. The complete F. californiensis and L. stylirostris mtDNAs are 15,975 and 15,988 bp long, containing the 37 common sequences and a control region of 990 and 999 bp, respectively. The gene order is identical to that of the tiger shrimp Penaeus monodon. Secondary structures for the 22 tRNAs are proposed and phylogenetic relationships for selected complete crustacean mitogenomes are included. Phylogenomic relationships among five shrimp show strong statistical support for the monophyly of the genus across the analysis. Litopenaeus species define a clade, with close relationship to Farfantepenaeus, and both clade with the sister group of Penaeus and Fenneropenaeus.

Nuclear and mitochondrial subunits from the white shrimp Litopenaeus vannamei F(0)F(1) ATP-synthase complex: cDNA sequence, molecular modeling, and mRNA quantification of atp9 and atp6

We studied for the first time the ATP-synthase complex from shrimp as a model to understand the basis of crustacean bioenergetics since they are exposed to endogenous processes as molting that demand high amount of energy. We analyzed the cDNA sequence of two subunits of the Fo sector from mitochondrial ATP-synthase in the white shrimp Litopenaeus vannamei. The nucleus encoded atp9 subunit presents a 773 bp sequence, containing a signal peptide sequence only observed in crustaceans, and the mitochondrial encoded atp6 subunit presents a sequence of 675 bp, and exhibits high identity with homologous sequences from invertebrate species. ATP9 and ATP6 protein structural models interaction suggest specific functional characteristics from both proteins in the mitochondrial enzyme. Differences in the steady-state mRNA levels of atp9 and atp6 from five different tissues correlate with tissue function. Moreover, significant changes in the mRNA levels of both subunits at different molt stages were detected. We discussed some insights about the enzyme structure and the regulation mechanisms from both ATP-synthase subunits related to the energy requirements of shrimp.

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

Background

Metazoan mitochondrial genomes usually consist of the same 37 genes. Such genes contain useful information for phylogenetic analyses and evolution modelling. Although complete mitochondrial genomes have been determined for over 1,000 animals to date, hydrothermal vent species have, thus far, remained excluded due to the scarcity of collected specimens.

Results

The mitochondrial genome of the hydrothermal vent galatheid crab Shinkaia crosnieri is 15,182 bp in length, and is composed of 13 protein-coding genes, two ribosomal RNA genes and only 18 transfer RNA genes. The total AT content of the genome, as is typical for decapods, is 72.9%. We identified a non-coding control region of 327 bp according to its location and AT-richness. This is the smallest control region discovered in crustaceans so far. A mechanism of cytoplasmic tRNA import was addressed to compensate for the four missing tRNAs. The S. crosnieri mitogenome exhibits a novel arrangement of mitochondrial genes. We investigated the mitochondrial gene orders and found that at least six rearrangements from the ancestral pancrustacean (crustacean + hexapod) pattern have happened successively. The codon usage, nucleotide composition and bias show no substantial difference with other decapods. Phylogenetic analyses using the concatenated nucleotide and amino acid sequences of the 13 protein-coding genes prove consistent with the previous classification based upon their morphology.

Conclusion

The present study will supply considerable data of use for both genomic and evolutionary research on hydrothermal vent ecosystems. The mitochondrial genetic characteristics of decapods are sustained in this case of S. crosnieri despite the absence of several tRNAs and a number of dramatic rearrangements. Our results may provide evidence for the immigrating hypothesis about how vent species originate.