The complete mitochondrial genome of Eriocheir japonica sinensis (Decapoda: Varunidae) and its phylogenetic analysis

The complete mitochondrial genome (mitogenome) can provide novel insights into understanding the mechanisms underlying mitogenome evolution. In this study, the complete mitogenome of Eriocheir japonica sinensis (Decapoda: Varunidae) was determined to be 16,378 bp, including 13 protein-coding genes (PCGs), two rRNA genes, 22 tRNA genes and a D-loop region. The AT skew of the E. j. sinensis mitogenome was slightly negative (−0.016), indicating a higher number of T compared with A nucleotides. The nucleotide composition of the mitogenome was also biased toward A + T nucleotides (71.6%). All PCGs were initiated by ATN codons. Eight of the 13 PCGs harbored the incomplete termination codon by T, or TA. All other tRNA genes displayed a typical clover-leaf structure of mitochondrial tRNA. The D-loop region of the E. j. sinensis mitogenome was 918 bp in length. Based on 13 PCGs, phylogenetic analysis confirmed the placement of E. j. sinensis within the Varunidae.     

Characterization of the complete mitochondrial genome of Drawida gisti (Metagynophora,Moniligastridae) and comparison with other Metagynophora species

Here, the complete mitochondrial genome (mitogenome) of Drawida gisti was sequenced and compared with the mitogenomes of other Metagynophora species. The circular mitogenome was 14,648 bp in length and contained two ribosomal RNA genes (rRNAs), 13 protein-coding genes (PCGs), and 22 transfer RNA genes (tRNAs). The types of constitutive genes and the direction of the coding strand that appeared in Drawida mitogenome were identical to those observed in other Metagynophora species, except for a missing lengthy non-coding region. The conservative relationships between Drawida species were supported by the overall analyses of 13 PCGs, two rRNAs, and 22 tRNAs. A comparison of the Metagynophora mitogenomes revealed that the ATP8 gene possessed the highest polymorphism among the 13 PCGs and two rRNAs. Phylogenetic analysis suggested that the Moniligastridae contained Drawida, which is a primitive Metagynophora group. Our study provides a step forward toward elucidating the evolutionary linkages within Drawida and even Metagynophora.     

Characterization of the complete mitochondrial genome of Hydrotaea spinigera (Diptera: Muscidae) with phylogenetic implications

Hydrotaea spinigera Stein is a muscid species with high forensic importance. To further supplement the genome-level features of related species, the complete mitochondrial genome of H. spinigera is amplified, sequenced, annotated, analyzed, and compared with 13 other species of the family Muscidae. This mitogenome is 15,517?bp in length, with a standard set of 13 protein-coding genes (PCGs), 22 transfer RNA (tRNA) genes, two ribosomal RNA (rRNA) genes, and a typical non-coding control region. The gene organization, base composition, and secondary structures of tRNA genes of this mitogenome are similar to other dipteran species. Phylogenetic analyses based on all the mitochondrial PCGs and rRNA genes for 14 muscid and one outgroup taxa are performed using Bayesian inference (BI) and maximum likelihood (ML) analysis. The inferred trees indicate that the Muscidae, the subfamily Muscinae, and the tribes Azeliini and Reinwardtiini are monophyletic, whereas the monophyly of the subfamily Azeliinae and the tribe Muscini is not supported.     

The complete mitochondrial genome of the sycamore lace bug Corythucha ciliata (Hemiptera: Tingidae)

The complete mitochondrial genome of the sycamore lace bug, Corythucha ciliata, was sequenced in this study. It represents the first sequenced mitogenome of family Tingidae in Heteroptera. The mitogenome of C. ciliata is 15,257 bp and contains 37 genes including 13 protein-coding genes (PCGs), 22 tRNA genes, two rRNA genes and a large non-coding region. Gene arrangement, nucleotide content, codon usage, and amino acid composition and asymmetry indicate a high degree of conservation with six other species of Cimicomorpha. The 13 PCGs initiated with ATN as the start codon and terminated with TAA, TA or T as stop codon. The evolutionary rate of each PCG was different, among which ATP8 showed the highest rate while ATP6 indicated the lowest rate. The 22 tRNAs genes apparently fold into a typical cloverleaf structure; however, the anticodon (TTC) of trnSer (AGN) differs from other Heteropteran insects. Secondary structure modeling of rRNA genes revealed similarity to other insects, except for two incomplete helices (H1648 and H2735) in lrRNA. The predicted secondary structure of lrRNA indicates 45 helices in six domains, whereas srRNA has 27 helices in three domains. Three potential stem–loops and two tandem repeats (–TCTAAT–) were identified in the A+T-rich region. Phylogenetic analysis indicated that C. ciliata is a sister group to other Heteroptera species based on analysis of the 13 PCGs.     

Characterization of the complete mitochondrial genome of the giant silkworm moth, Eriogyna pyretorum (Lepidoptera: Saturniidae)

The complete mitochondrial genome (mitogenome) of Eriogyna pyretorum (Lepidoptera: Saturniidae) was determined as being composed of 15,327 base pairs (bp), including 13 protein-coding genes (PCGs), 2 rRNA genes, 22 tRNA genes, and a control region. The arrangement of the PCGs is the same as that found in the other sequenced lepidopteran. The AT skewness for the E. pyretorum mitogenome is slightly negative (-0.031), indicating the occurrence of more Ts than As. The nucleotide composition of the E. pyretorum mitogenome is also biased toward A + T nucleotides (80.82%). All PCGs are initiated by ATN codons, except for cytochrome c oxidase subunit 1 and 2 (cox1 and cox2). Two of the 13 PCGs harbor the incomplete termination codon by T. All tRNA genes have a typical clover-leaf structure of mitochondrial tRNA, with the exception of trnS1(AGN) and trnS2(UCN). Phylogenetic analysis among the available lepidopteran species supports the current morphology-based hypothesis that Bombycoidea, Geometroidea, Notodontidea, Papilionoidea and Pyraloidea are monophyletic. As has been previously suggested, Bombycidae (Bombyx mori and Bombyx mandarina), Sphingoidae (Manduca sexta) and Saturniidae (Antheraea pernyi, Antheraea yamamai, E. pyretorum and Caligula boisduvalii) formed a group.     

The Complete Mitochondrial Genome of the Beet Webworm,Spoladea recurvalis (Lepidoptera: Crambidae) and Its Phylogenetic Implications

The complete mitochondrial genome (mitogenome) of the beet webworm, Spoladea recurvalis has been sequenced. The circular genome is 15,273 bp in size, encoding 13 protein-coding genes (PCGs), two rRNA genes, and 22 tRNA genes and containing a control region with gene order and orientation identical to that of other ditrysian lepidopteran mitogenomes. The nucleotide composition of the mitogenome shows a high A+T content of 80.9%, and the AT skewness is slightly negative (-0.023). All PCGs start with the typical ATN codons, except for COX1, which may start with the CGA codon. Nine of 13 PCGs have the common stop codon TAA; however, COX1, COX2 and ND5 utilize the T nucleotide and ND4 utilizes TA nucleotides as incomplete termination codons. All tRNAs genes are folded into the typical cloverleaf structure of mitochondrial tRNAs, except for the tRNA^Ser(AGY) gene, in which the DHU arm fails to form a stable stem-loop structure. A total of 157 bp intergenic spacers are scattered in 17 regions. The overlapping sequences are 42 bp in total and found in eight different locations. The 329 bp AT-rich region is comprised of non-repetitive sequences, including the motif ATAG, which is followed by a 14 bp poly-T stretch, a (AT₁₁ microsatellite-like repeat, which is adjacent to the motif ATTTA, and a 9 bp poly-A, which is immediately upstream from the tRNA^Met gene. Phylogenetic analysis, based on 13 PCGs and 13 PCGs+2 rRNAs using Bayesian inference and Maximum likelihood methods, show that the classification position of Pyraloidea is inconsistent with the traditional classification. Hesperioidea is placed within the Papilionoidea rather than as a sister group to it. The Pyraloidea is placed within the Macrolepidoptera with other superfamilies instead of the Papilionoidea.     

The complete mitochondrial genome of <Emphasis Type="Italic">Damora sagana</Emphasis> and phylogenetic analyses of the family Nymphalidae

The monotypic genus Damora (Nymphalidae, Heliconiinae) contains a single species, Damora sagana, which is widely distributed across southern China. Herein, its complete mitogenome was sequenced to further understand lepidopteran mitogenome characteristics, reconstruct the nymphalid family phylogeny, and infer the subdivision of Heliconiinae species. The circular mitogenome was 15,151 bp long, abundant in A and T, and comprised of 13 protein-coding genes (PCGs), 22 transfer RNA (tRNA) genes, 2 ribosomal RNA (rRNA) genes, and one control region with a gene arrangement typical of lepidopteran mitogenomes. ATN codons initiated all PCGs, except cytochrome c oxidase subunit 1 (COX1), which was initiated by a CGA sequence as has been observed in other lepidopterans. Three PCGs (COX1, COX2 and ND4) employed a single T termination signal, whereas others had the typical complete termination codon (TAA). All tRNA genes were folded into the typical cloverleaf structure except for tRNA-Ser (AGN). The A+T-rich region included the conserved motif ‘ATAGA’ followed by a 17 bp poly-T stretch, which was also observed in tribe Argynnini mitogenomes. A phylogenetic tree was constructed via multiple methods using the 13 PCGs data of D. sagana and other available mitogenomes of nymphalid species. All three phylogenetic trees yielded the same topology. These results were consistent with those from previous studies of most major nymphalid groups, except those regarding tribe subdivision in certain subfamilies such as Argynnini + (Acraeini + Heliconiini) for Heliconiine. Furthermore, our analyses identified that the genus Cethosia was grouped with the genus Acraea composing the tribe Acraeini with strong support.     

The first complete mitogenome of Picumnus innominatus (Aves,Piciformes, Picidae) and phylogenetic inference within the Picidae

The avian family Picidae, which is nearly global in distribution, contains the piculets (Picumninae and Nesoctitinae), the woodpeckers (Picinae), and the wrynecks (Jynginae). However, the phylogenetic relationships within the Picidae remain obscure for most genera. In the present study, the complete mitochondrial genome of Picumnus innominatus was determined and described, which was the first complete mitogenome reported in the Picumnus. The circular mitogenome of P. innominatus was 17,180 bp in size and consisted of 13 protein-coding genes (PCGs), 22 tRNAs, 2 rRNA genes, a control region (CR), and a noncoding region. The gene order and arrangement of the P. innominatus mitogenome were identical to other mitogenomes of the Picidae. Moreover, strikingly large tandem repeats were found in the noncoding region of the P. innominatus mitogenome, which have not yet been covered in other picid species to date. At the family level (Picidae), the highest dN/dS ratio was detected for the ND1 gene (1.38726) among 13 PCGs, indicating that positive selection was powerful for this gene. Bayesian and Maximum Likelihood phylogenetic analyses based on the combination of 12S rRNA and CYTB gene supported strongly that the Picumninae is monophyletic.     

The complete mitochondrial genome of cricket Sclerogryllus punctatus (Orthoptera: Gryllidae) and phylogenetic analysis

The crickets of genus Sclerogryllus Gorochov, 1985 belongs to subfamily Sclerogryllinae of family Gryllidae. In this study, we report the first complete mitogenome sequences of the genus Sclerogryllus, and analyze the features of mitogenomes of S. punctatus. The mitogenome of S. punctatus was 15,438 bp and consisted of 37 genes, coding for 13 proteins, 2 ribosomal RNA (rRNA) and 22 transfer RNA (tRNA), and a control region. S. punctatus shares the arrangement of trnE-trnS-trnN with most mitogenomes of Grylloidea. Besides, the tRNAs possess the typical cloverleaf secondary structure except for the trnS1 (AGN) gene. The phylogenetic analysis using 13 protein-coding genes and 2 rRNA represents that genus Sclerogryllus is included in subfamily Gryllinae. Our results uncover the phylogenetic position of genus Sclerogryllus by mitogenome data within the family Gryllidae.     

Characterization of the complete mitochondrial genome of Diaphania pyloalis (Lepidoptera: Pyralididae)

The complete mitochondrial genome (mitogenome) of Diaphania pyloalis (Lepidoptera: Pyralididae) was determined to be 15,298 bp and has the typical gene organization of mitogenomes from lepidopteran insects. It consists of 13 protein-coding genes (PCGs), two rRNA genes, 22 tRNA genes and an A + T-rich region. The A + T content of this mitogenome is 80.83% and the AT skew is slightly positive. All PCGs are initiated by ATN codons, except for cytochrome c oxidase subunit 1 (cox1) gene which is initiated by CGA. Only the cox2 gene has an incomplete stop codon consisting of just a T. All the tRNA genes display a typical clover-leaf structure of mitochondrial tRNA. The A + T-rich region of the mitogenome is 332 bp in length, including several common features found in lepidopteran mitogenomes. Phylogenetic analysis showed that the D. pyloalis is close to Pyralididae.     

The complete mitogenome of Apocheima cinerarius (Lepidoptera: Geometridae: Ennominae) and comparison with that of other lepidopteran insects

The complete mitochondrial genome (mitogenome) of a female flightless geometrid moth Apocheima cinerarius was found to be 15,722 bp in length, containing 13 protein-coding genes (PCGs), 22 transfer RNA (tRNA) genes, 2 ribosomal RNA (rRNA) genes, and a control region. The A + T content of the complete mitogenome is 80.83%. The AT skew value ([A − T] / [A + T]) is 0.027. The 13 PCGs of the mitogenome start with typical ATN codons, except for cox1 with the start codon CGA. All the tRNA genes have typical cloverleaf secondary structures, except for trnSer(AGN). The secondary structures of rrnL and rrnS were predicted. Six structural domains including conserved regions (IV, V) and variable regions (I, II, III, VI) were identified in the secondary structure of rrnL. The secondary structure of rrnS consists of 3 structural domains. The control region of A. cinerarius begins with conserved motifs of “ATAGA” + 19-bp poly T. It also contains a microsatellite-like (TA)₂₆, a stem-and-loop structure, and a poly-A stretch. Phylogenetic analysis showed that Geometroidea is more closely related to Bombycoidea than to Noctuoidea. A. cinerarius is more closely related to Biston panterinaria than to Phthonandria atrilineata, which is in accordance with the conventional morphology-based classification.     

The complete mitochondrial genome of fall armyworm <Emphasis Type="Italic">Spodoptera frugiperda</Emphasis> (Lepidoptera:Noctuidae)

The mitochondrial genome (mitogenome) can provide important information for understanding molecular evolution and phylogenetic analyses. The complete mitogenome of Spodoptera frugiperda (Lepidoptera:Noctuidae) was determined to be 15,365 bp in length and has the typical gene order found in Noctuidae mitogenomes, it includes 13 protein-coding genes (PCGs), two rRNA genes, 22 tRNA genes and a A+T-rich region. The nucleotide composition was biased toward A+T nucleotides (81.09 %) and the AT skew of this mitogenome was slightly positive (0.004). All PCGs were initiated by ATN codons, except for cytochrome c oxidase subunit 1 (cox1) gene which was initiated by CGA. Eight of the 13 PCGs have the incomplete termination codon, T or TA. All the tRNA genes displayed the typical clover-leaf structure of mitochondrial tRNAs, with the exception of trnS1 (AGN). The A+T-rich region was 328 bp in length and consisted of several features common to the Noctuidae insects. Phylogenetic analysis showed that the S. frugiperda was within the Noctuidae.     

Complete mitochondrial genome of Pycnonotus xanthorrhous (Passeriformes,Pycnonotidae) and phylogenetic consideration

The complete mitochondrial genome (mitogenome) of Pycnonotus xanthorrhous was sequenced via next generation sequencing. The full length of the circular genome is 16,952 bp. It consists of 37 typical animal mitochondrial genes including 13 protein-coding genes (PCGs), 22 transfer RNA (tRNA) and 2 ribosomal RNA (rRNA) genes. P. xanthorrhous also contains one control region (CR) and one pseudo-control region, and shares the identical gene arrangements with sequenced Pycnonotus spp. which differs from the typical vertebrates gene order. Phylogenetic analyses indicates that Passerida sensu stricto contains three major clades and the core Sylvioidea form a monophyletic group. Furthermore, we investigated the evolution of control region within this lineage and revealed the multiple independent origins of duplicate control region.     

Characterization of the complete mitochondrial genome of Bombyx mori strain H9 (Lepidoptera: Bombycidae)

The complete mitochondrial genome (mitogenome) of Bombyx mori strain H9 (Lepidoptera: Bombycidae) is 15,670 base pairs (bp) in length, encoding 13 protein-coding genes (PCGs), two rRNA genes, 22 tRNA genes and a control region. The nucleotide composition of the genome is highly A + T biased, accounting for 81.31%, with a slightly positive AT skewness (0.059). The arrangement of 13 PCGs is similar to that of other sequenced lepidopterans. All the PCGs are initiated by ATN codons, except for the cytochrome c oxidase subunit 1 (cox1) gene, which is proposed by the TTAG sequence as observed in other lepidopterans. Unlike the other PCGs, the cox1 and cytochrome c oxidase subunit 2 (cox2) genes have incomplete stop codons consisting of just a T. All tRNAs have typical structures of insect mitochondrial tRNAs, which is different from other sequenced lepidopterans. The structure of A + T-rich region is similar to that of other sequenced lepidopterans, including non-repetitive sequences, the ATAGA binding domain, a 18 bp poly-T stretch and a poly-A element upstream of transfer RNA M (trnM) gene. Phylogenetic analysis shows that the domesticated silkmoth B. mori originated from the Chinese Bombyx mandarina.     

The complete mitochondrial genome of Parasesarma pictum (Brachyura: Grapsoidea: Sesarmidae) and comparison with other Brachyuran crabs

Mitochondrial DNA (mtDNA) is an extrachromosomal genome which can provide important information for evolution and phylogenetic analysis. In this study, we assembled a complete mitogenome of a crab Parasesarma pictum (Brachyura: Grapsoidea: Sesarmidae) from next generation sequencing reads at the first time. P. pictum is a mudflat crab, belonging to the Sesarmidae family (subfamily Sesarminae), which is perched on East Asia. The 15,716 bp mitogenome covers 13 protein-coding genes (PCGs), 22 transfer RNA genes (tRNAs), 2 ribosomal RNA genes (rRNAs), and one control region (CR). The control region spanns 420 bp. The genome composition was highly A+T biased 75.60% and showed negative AT-skew (−0.03) and negative GC-skew (−0.47). Compared with the ancestor of Brachyura, the gene order of Sesarmidae has several differences and the gene order of P. pictum is typical for mitogenomes of Sesarmidae. Phylogenetic tree based on nucleotide sequences of mitochondrial 13 PCGs using BI and ML determined that P. pictum has a sister group relationship with Parasesarma tripectinis and belongs to Sesarmidae.     

The complete nucleotide sequence of the mitochondrial genome of <Emphasis Type="Italic">Phthonandria atrilineata</Emphasis> (Lepidoptera: Geometridae)

Using long-polymerase chain reaction (Long-PCR) method, we determined the complete nucleotide sequence of the mitochondrial genome (mitogenome) of Phthonandria atrilineata. The complete mtDNA from P. atrilineata was 15,499 base pairs in length and contained 13 protein-coding genes (PCGs), 2 rRNA genes, 22 tRNA genes, and a control region. The P. atrilineata genes were in the same order and orientation as the completely sequenced mitogenomes of other lepidopteran species. The nucleotide composition of P. atrilineata mitogenome was biased toward A + T nucleotides (81.02%), and the 13 PCGs show different A + T contents that range from 73.25% (cox1) to 92.12% (atp8). Phthonandria had the canonical set of 22 tRNA genes, that fold in the typical cloverleaf structure described for metazoan mt tRNAs, with the unique exception of trnS(AGN). The phylogenetic relationships were reconstructed with the concatenated sequences of the 13 PCGs of the mitochondrial genome, which confirmed that P. atrilineata is most closely related to the superfamily Bombycoidea.     

The complete mitochondrial genome of Pardosa pusiola (Araneae,Lycosidae) and its phylogenetic implications

The mitochondrial genome (mitogenome) is useful for identification and phylogenetic analyses among arthropods, but there are no sufficient mitogenome data for wolf spiders. To enrich the mitogenome database of wolf spiders, the complete mitogenome of Pardosa pusiola was sequenced by high-throughput sequencing. It is 14,284 bp, comprising 13 protein-coding genes (PCGs), 22 transfer RNA genes (tRNAs), two ribosomal RNA genes (rRNAs), and a control region (CR). It represents a high bias toward A and T nucleotides with an A + T content of 76.49%. The mitogenome exhibited a negative AT skew (−0.13) and a positive GC skew (0.32). Most PCGs started with ATN codons and ended with TAA, TAG, or an incomplete T. In addition, most tRNAs had aberrant secondary structures with the absence of DHU arm or TΨC arm. Analysis performed with CREx software demonstrated that large-scale rearrangements of tRNAs were observed in the mitogenome of P. pusiola as compared with the putative ancestral mitogenome. The Bayesian inference (BI) and maximum likelihood (ML) phylogenetic trees based on the 13 PCGs of 25 spiders had the same topology, which could be presented as (Araneidae + (Agelenidae + (Dictynidae + Desidae)) + (Salticidae + (Thomisidae + (Oxyopidae + (Pisauridae + Lycosidae))))). This study offers a useful genetic resource for the taxonomy and phylogeny of spiders.     

Complete mitochondrial genome of Dacus conopsoides (Insecta: Tephritidae) with tRNA gene duplication and molecular phylogeny of Dacini tribe

To date there is only a single report on the complete mitochondrial genome of the Dacus fruit flies. We report here the whole mitogenome of Dacus conopsoides with first report of tRNA gene duplication in tephritid fruit flies determined using next-generation sequencing and discuss the molecular phylogeny of Dacini tribe. It had a total length of 15,852 bp, comprising 13 protein coding genes, 2 rRNA genes, 23 tRNA genes, and a non-coding region (A + T-rich control region). The 65-bp trnF gene was duplicated, and the 68-bp trnE gene was partially duplicated resulting in a 31-bp pseudogene. The cloverleaf structure for trnN, trnH, and trnF lacked the TΨC-loop, while trnS lacked the D-stem. The start codons for the protein coding genes included 6 ATG, 3 ATC, 2 ATA, and 1 each of ATT and TCG. Seven PCGs had TAA stop codon, two had TAG and four had incomplete T stop codon. Molecular phylogeny based on 15 mt-genes (13 PCGs +2 rRNA genes) and 30 taxa of Tephritidae indicated D. conopsoides forming a monophyletic sister group with D. longicornis supported by high bootstrap value. The lineage containing also the monophyletic genus Zeugodacus. The Dacini and Ceratitidini tribes of the subfamily Dacinae were monophyletic but the subfamilies Dacinae and Trypetinae were paraphyletic. A broader taxa sampling of the Tephritidae is needed to better elucidate the phylogenetics and systematics of the tribes and subfamilies of tephritid fruit flies.     

Complete mitochondrial genome of the atlas moth, Attacus atlas (Lepidoptera: Saturniidae) and the phylogenetic relationship of Saturniidae species

Mitochondrial genome (mitogenome) can provide information for genomic structure as well as for phylogenetic analysis and evolutionary biology. In this study, we present the complete mitogenome of the atlas moth, Attacus atlas (Lepidoptera: Saturniidae), a well-known silk-producing and ornamental insect with the largest wing surface area of all moths. The mitogenome of A. atlas is a circular molecule of 15,282 bp long, and its nucleotide composition shows heavily biased towards As and Ts, accounting for 79.30%. This genome comprises 13 protein-coding genes (PCGs), two ribosomal RNA genes (rRNAs), 22 transfer RNA genes (tRNAs), and an A + T-rich region. It is of note that this genome exhibits a slightly positive AT skew, which is different from the other known Saturniidae species. All PCGs are initiated by ATN codons, except for COI with CGA instead. Only six PCGs use a common stop codon of TAA or TAG, whereas the remaining seven use an incomplete termination codon T or TA. All tRNAs have the typical clover-leaf structure, with an exception for tRNA^Ser(AGN). The A. atlas A + T-rich region contains non-repetitive sequences, but harbors several features common to the Bombycoidea insects. The phylogenetic relationships based on Maximum Likelihood method provide a well-supported outline of Saturniidae, which is in accordance with the traditional morphological classification and recent molecular works.     

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,359 bp 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.