Mitochondrial genome of Helice tientsinensis (Brachyura: Grapsoidea: Varunidae): Gene rearrangements and higher-level phylogeny of the Brachyura

The analyses of the complete mitochondrial genomes of three crabs revealed novel gene rearrangements and phylogenetic relationships of Brachyura

BACKGROUND

Brachyura crab is the largest branch of Decapoda crustacean. Phylogenetic relationships within Brachyura remain controversial to be investigated. The mitochondrial genome (mitogenome) is an important molecular marker for studying the phylogenetic relationships of Brachyura.

METHODS AND RESULTS

To understand the phylogeny of Brachyura, the three complete mitogenomes from Charybdis annulata, Leptodius exaratus, and Spider crab were sequenced and annotated. Their full length was 15,747, 15,716, and 16,608 bp long, respectively. The first two crabs both contained 13 protein-coding genes (PCGs), two rRNA genes, 22 tRNA genes and a control region. However, Spider crab contained 13 PCGs, two rRNA genes, 25 tRNA genes and a control region. The mitogenomes of each of the three crabs exhibited high AT content (67.8%, 69.1%, and 70.8%), with negative AT skews (-0.014, - 0.028, and - 0.017) and GC skews (-0.269, - 0.286, and - 0.341). The gene order of C. annulata was identical to the ancestor of Brachyura. Compared with the ancestor of Brachyura, L. exaratus exhibited the gene rearrangements of Val (V)-rrnS-control region, and Spider crab had the four copies of Lys (K). Phylogenetic analyses indicated that C. annulata belonged to Portunidae family, Portunoidea superfamilies, L. exaratus belonged to Xanthidae family, Xanthoidea superfamilies, and Spider crab belonged to Mithracidae family, Majoidea superfamilies. Phylogenetic analyses showed that the two species (Somanniathelphusa boyangensis and Huananpotamon lichuanense) belonging to the Potamoidea were sister groups to the Thoracotremata, thus supporting the conclusion that Heterotremata is polyphyletic.

CONCLUSION

The results of this study enriched the crab mitogenome database and enabled us to better understand the phylogenetic relationships of Brachyura.

The complete mitochondrial genome of Helicana japonica (Crustacea, Decapoda, Varunidae) from South Korea

Helicana japonica mainly inhabits burrowed holes in the mudflats and intertidal zones. Specimens from the Republic of Korea were collected and whole genomic DNA from the cheliped muscle tissue was extracted. We determined the complete mitochondrial genome using Illumina HiSeq X Ten. The mitogenome is 16,535 bp in length and consists of 13 protein-coding genes, 2 rRNA genes, and 22 tRNA genes. A phylogenetic tree was reconstructed using the maximum-likelihood of phylogeny methods. H. japonica formed a sister clade with Helicana wuana, which is another Helicana species.

The characterization, expression and activity analysis of three superoxide dismutases in Eriocheir hepuensis under azadirachtin stress

Superoxide dismutase (SOD) can effectively eliminate of excess ROS, which causes oxidative damage to lipids, proteins, and DNA. In this study, we cloned the CuZn-SOD, cMn-SOD1, and cMn-SOD2 genes in Eriocheir hepuensis, and found that the coding sequence (CDS) lengths were 627 bp, 861 bp and 1062 bp, which encoded 208, 286, and 353 amino acids, respectively. Phylogenetic analysis indicated that all SOD genes were evolutionarily conserved, while cMn-SOD2 had an extra gap (67 amino acids) in the conserved domain compared with cMn-SOD1 without huge changes in the tertiary structure of the conserved domain, suggesting that cMn-SOD2 may be a duplicate of cMn-SOD1. qRT-PCR showed that the three SOD genes were widely expressed in all the tested tissues, CuZn-SOD and cMn-SOD1 were mostly expressed in the hepatopancreas, while cMn-SOD2 was mostly expressed in thoracic ganglia. Under azadirachtin stress, the oxidation index of surviving individuals, including the T-AOC, SOD activity, and MDA contents increased in the early stage and then remained steady except for a decrease in MDA contents in the later stage. qRT-PCR showed that the three SOD genes displayed the same trends as SOD activity in surviving individuals, and the highest expressions of CuZn-SOD in the hepatopancreas, heart, and gill were 14.16, 1.41, and 30.87 times that of the corresponding control group, respectively. The changes were 1.35, 5.77 and 3.33 fold for cMn-SOD1 and 1.62, 1.71 and 1.79 fold for cMn-SOD2, respectively. However, the activity and expression of SOD genes in dead individuals were lower than that observed in surviving individuals. These results reveal that SOD plays a significant role in the defence against azadirachtin-induced oxidative stress.

Comparative Mitochondrial Genome Analyses of Sesarmid and Other Brachyuran Crabs Reveal Gene Rearrangements and Phylogeny

Mitochondrial genomes (mitogenomes) are important for understanding molecular evolution and phylogenetic relationships. The complete mitogenome of Perisesarma bidens was determined, which is 15,641 bp in length. The A + T content of P. bidens mitogenome was 74.81%. The AT skew was slightly negative (-0.021). The 22 tRNAs ranged from 65 to 73 bp and were highly A + T biased. All tRNA genes had typical cloverleaf structures, except for the trnS1 gene, which lacked a dihydrouridine (DHU) arm. The gene order within the P. bidens mitogenome was identical to the pancrustacean ground pattern, except for the translocation of the trnH. Additionally, the gene order of trnI-trnQ-trnM in pancrustacean ground pattern became trnQ-trnI-trnM in P. bidens. Phylogenetic analyses supported the inclusion of P. bidens in Sesarmidae and the promotion of Sesarminae to Sesarmidae. The results will help us to better understand the status and evolutionary history of Grapsoidea crabs.

Phylogenetic relationships of Grapsoidea and insights into the higher phylogeny of Brachyuran

Decapoda is one of the most diverse crustacean orders, and has become an important research subject. However, the phylogenetic relationships among the main lineages of Decapoda remain uncertain, especially in the order Brachyura. Herein, we sequenced the whole mitochondrial genome of V. litterata and constructed a phylogenetic tree to understand its phylogenetic relationships with other species. The results showed that the mitochondrial genome of V. litterata was generally similar to mitogenomes of Metazoa reported in the literature, with a size of 16,247 bp, 37 genes, and a control region. Both AT-skew and GC-skew were negative, indicating more abundant Cs and Ts than Gs and As. The gene arrangement of V. litterata is identical to those of Eriocheir hepuensis, Cyclograpsus granulosus, Hemigrapsus sanguineus, Helicana wuana, and Helice tientsinensis but differs from the pancrustacean ground pattern and typical arrangement of Brachyuran crabs. Phylogenetic reconstruction showed that V. litterata belongs to the Varunidae.

The complete mitochondrial genome of Calappa bilineata: The first representative from the family Calappidae and its phylogenetic position within Brachyura

In this study, we determined the complete mitogenome sequence of Calappa bilineata, which is the first mitogenome of Calappidae up to now. The total length is 15,606 bp and includes 13 protein-coding genes, 22 transfer RNAs, two ribosomal RNAs and one control region. The genome composition is highly A + T biased (68.7%), and exhibits a negative AT-skew (-0.010) and GC-skew (-0.267). As with other invertebrate mitogenomes, the PCGs start with the standard ATN and stop with the standard TAN codons or incomplete T. Phylogenetic analysis showed that C. bilineata was most closely related to Matuta planipes (Matutidae), and these two species formed a sister clade, constituting a Calappoidea group and forming a sister clade with part of Eriphioidea. The existence of the polyphyletic families raised doubts over the traditional classification system. These results will help to better understand the features of the C. bilineata mitogenome and lay foundation for further evolutionary relationships within Brachyura.

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.

Transcriptomic analysis of immune-related genes in the lipopolysaccharide-stimulated hepatopancreas of the mudflat crab Helice tientsinensis

The mudflat crab Helice tientsinensis is one of the most economically important aquaculture species in China. Nevertheless, it is susceptible to various diseases caused by viruses, bacteria and rickettsia-like organisms. A better understanding of the immune system and genes related to the responses to bacterial and viral infection is required. Herein, the hepatopancreas transcriptome of H. tientsinensis was analyzed by comparing control and lipopolysaccharide (LPS)-stimulated RNA-Seq data, yielding 91,885,038 bp and 13.78 Gb of clean reads. Following assembly and annotation, 93,207 unigenes with an average length of 883 bp were identified, of which 31,674 and 13,700 were annotated in Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) databases, respectively. Following LPS, 4845 differentially expressed genes (DEGs) were identified, of which 2491 and 2354 were up- and down-regulated, respectively. To further investigate immune-related DEGs, KEGG enrichment analysis identified immune response pathways, most notably the peroxisome and Toll-like receptor signaling pathways. Quantitative real time-PCR (qRT-PCR) confirmed the up-regulation of a random selection of DEGs. This systematic transcriptomic analysis of the innate immune pathway in H. tientsinensis expands our understanding of the immune system in crabs.

The complete mitochondrial genome of Varuna yui (Decapoda: Brachyura: Varunidae) and its phylogeny

The complete mitochondrial genome plays an important role in the research on phylogenetic relationship. Here, we reported the first complete mitochondrial genome sequence of Varuna yui Hwang & Takeda, 1986 (Varunidae). The complete mtDNA (15,915 bp in length) consisted of 13 protein-coding genes, 22 tRNAs, two rRNA genes, and a control region. The gene arrangement was identical to those observed in the Varunidae species. The phylogenetic analysis suggested that V. yui had close relationship with other Varunidae species (Helicetient sinensis, Eriocher sinesis, etc.). The newly described genome may facilitate further comparative mitogenomic analysis within Varunidae species.

The complete mitochondrial genome of Huananpotamon lichuanense (Decapoda: Brachyura) with phylogenetic implications for freshwater crabs

In the present study, we determined the complete mitochondrial genome of Huananpotamon lichuanense (Decapoda: Brachyura) for the first time. The genome is 15,380bp in length and typically consists of 37 genes. When the gene order was compared to the ancestral crustacean type, two tRNA genes (tRNA^His and tRNA^Gln) were rearranged in H. lichuanense, and the translocation of tRNA^Gln appeared only in Potamoidea crabs, such as Geothelphusa dehaani and Sinopotamon xiushuiense, supporting the monophyly of the Potamoidea superfamily. Thirteen protein-coding genes and 2 rRNA genes were divided into five complexes to perform the phylogenetic analysis, and the results showed that the trees constructed by complex I (ND1-ND6 and ND4L), complex IV (COX1-COX3) and rRNA genes better accord with the morphological classification system, suggesting that molecular markers of higher-level phylogeny can be developed in these three complexes in the future. The estimated divergence time for freshwater crabs is approximately 133.58Ma, and G. dehaani from Japan diverged from the freshwater crabs of mainland China approximately 60.66Ma. A selective pressure analysis based on current data revealed obviously increasing dN/dS ratios (except for ATP6 and ND4L) of freshwater crabs, and the accumulation of nonsynonymous mutations suggests that terrestrial habitats provide a relatively relaxed selective pressure environment for this group.

Characterisation of the complete mitochondrial genome of Helice wuana (Grapsoidea: Varunidae) and comparison with other Brachyuran crabs

The mitochondrial genome (mitogenome) provides important information for phylogenetic analysis and understanding evolutionary origins. Herein, we sequenced, annotated, and characterised the mitogenome of the crab Helice wuana to better understand its molecular evolution and phylogeny. The 16,359bp mitogenome includes 13 protein-coding genes (PCGs), 22 transfer RNA (tRNA) genes, two ribosomal RNA (rRNA) genes and one control region. The genome composition is highly A+T biased 68.42%, and exhibits a negative AT-skew (-0.036) and GC-skew (-0.269) among Brachyura crabs. Gene rearrangements were detected, as was tandem duplication followed by random loss, which explains the translocation of mitochondrial genes. Phylogenetic analysis showed that H. wuana and H. tientsinensis clustered on one branch with high nodal support values. These results confirm that the placement of H. wuana within the Varunidae family of Thoracotrematan crabs. This study will provided a better understanding for gene rearrangements and crab evolution in the further.