Complete mitochondrial genome of Neochauliodes parasparsus (Megaloptera: Corydalidae) with phylogenetic consideration

We sequenced the complete mitochondrial genome (mitogenome) of Neochauliodes parasparsus. The 15,995-bp mitogenome contained the standard set of 13 protein-coding genes, 22 transfer RNA genes (tRNAs), 2 ribosomal RNA genes (rRNAs), and a putative control region, with a gene arrangement that was identical to that reported for most other megalopteran species. We also predicted the secondary structure of all the RNA genes and analysed the preferred codon usage of the protein-coding genes. The putative 1265-bp control region contained two tandem repeated regions and several microsatellite-like elements. The phylogenetic analysis of available neuropteridan mitogenomes, based on the 13 protein-coding genes, appeared to support the current view of the neuropteridan phylogeny, and among the Neochauliodes spp., N. parasparsus was the most closely related to N. punctatolosus.     

The first mitochondrial genome for Brahmin moths (Brahmaeidae) and implications for the higher phylogeny of Bombycoidea

Bombycoidea comprises 10 families and 4723 species, and the phylogenetic relationships among families are still in debate. In this study, we have determined the complete mitochondrial genome (mitogenome) of Brahmaea porphyria. The 15,429-bp mitogenome contains a common set of 37 mitochondrial genes including 13 protein-coding genes, 22 transfer RNA genes (tRNAs), two ribosomal RNA genes (rRNAs) and an inferred control region, and shares the conserved gene rearrangement (trnM-trnI-trnQ) in most ditrysian mitogenomes. Moreover, we analysed the secondary structure for all the tRNA genes of B. porphyria and the preference of codon usage in the PCGs of B. porphyria. The putative 373-bp control region (CR) possesses three types of conserved elements, including ATAGA, Ploy-T stretch, and microsatellite-like elements. A phylogenetic analysis among available Bombycoidea mitogenomes using the concatenated 37 mitochondrial genes appears to support the hypothesis of (Sphingidae+Bombycidae)+Saturniidae and the relatively basal phylogenetic position of Brahmaeidae within Bombycoidea.     

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.     

Comparative mitochondrial genomics and phylogenetics among species of the Oriental dobsonfly genus Neoneuromus van der Weele, 1909 (Megaloptera: Corydalidae)

The dobsonfly genus Neoneuromus van der Weele, 1909 represents one of the megalopteran lineages with large body-size, and comprises 13 species all endemic to the Oriental region. In the present study, the mitochondrial genomes (mitogenomes) of 12 species of this genus were determined and analyzed for the first time. The mitogenome of the genus-type, Neoneuromus fenestralis (McLachlan, 1869), as a representative of these congeneric species, is herein described in detail. All of the mitogenomes of Neoneuromus are composed of 37 encoded genes and a control region. The evolutionary rates of the protein coding genes (PCGs) of the 13 species of Neoneuromus as well as different genera of Corydalinae are estimated. All the transfer RNA genes (tRNAs) have the typical clover-leaf secondary structure except trnS1 (AGN). Interspecific relationships within Neoneuromus were reconstructed based on different datasets generated from mitogenomic sequences. Our results indicate that tRNA and ribosomal RNA genes (rRNAs) of Neoneuromus species contribute phylogenetic signal when being concatenated with the PCGs, thus should be kept during phylogenetic analysis. The results sheds light on understanding the evolution of these aquatic and predatory insects.     

The first complete mitochondrial genome of the Indian Tent Turtle,Pangshura tentoria (Testudines: Geoemydidae): Characterization and comparative analysis

The characterization of a complete mitogenome is widely used in genomics studies for systematics and evolutionary research. However, the sequences and structural motifs contained within the mitogenome of Testudines taxa have rarely been examined. The present study decodes the first complete mitochondrial genome of the Indian Tent Turtle, Pangshura tentoria (16,657 bp) by using next‐generation sequencing. This denovo assembly encodes 37 genes: 13 protein‐coding genes (PCGs), 22 transfer RNA (tRNAs), two ribosomal RNA, and one control region (CR). Most of the genes were encoded on majority strand, except for one PCG (NADH dehydrogenase subunit 6) and eight tRNAs. Most of the PCGs were started with an ATG initiation codon, except for Cytochrome oxidase subunit 1 with “GTG” and NADH dehydrogenase subunit 5 with “ATA.” The termination codons, “TAA” and “AGA” were observed in two subunits of NADH dehydrogenase gene. The relative synonymous codon usage analysis revealed the maximum abundance of alanine, isoleucine, leucine, and threonine. The nonsynonymous/synonymous ratios were <1 in all PCGs, which indicates strong negative selection among all Geoemydid species. The study also found the typical cloverleaf secondary structure in most of the tRNA genes, except for serine with the lack of the conventional DHU arm. The comparative study of Geoemydid mitogenomes revealed the occurrence of tandem repeats was frequent in the 3′ end of CR. Further, two copies of a unique tandem repeat “TTCTCTTT” were identified in P. tentoria. The Bayesian and maximum‐likelihood phylogenetic trees using concatenation of 13 PCGs revealed the close relationships of P. tentoria with Batagur trivittata in the studied dataset. All the Geoemydid species showed distinct clustering with high bootstrap support congruent with previous evolutionary hypotheses. We suggest that the generations of more mitogenomes of Geoemydid species are required, to improve our understanding of their in‐depth phylogenetic and evolutionary relationships.     

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.     

Fourteen complete mitochondrial genomes of butterflies from the genus Lethe (Lepidoptera,Nymphalidae, Satyrinae) with mitogenome-based phylogenetic analysis

The mitochondrial genome (mitogenome) can help us understand the phylogenetic relationships within the genus Lethe and the subfamily Satyrinae. In this study, we sequenced the complete mitogenomes of 14 Lethe species, which range in size from 15,225 to 15,271 bp, with both 37 genes (13 PCGs, 22 tRNAs, 2 rRNAs) and a noncoding A + T-rich region. The gene arrangement and orientation is similar to typical mitogenomes of Lepidoptera. The Ka/Ks ratio shows that cox1 has the slowest evolutionary rate. The secondary structure of trnN lacks the Pseudouracil loop (TψC loop) in most Lethe species. The inferred phylogenetic analyses show that Lethe is a well-supported monophyletic group, and reveal 2 major clades within the genus Lethe, which is consistent with previous morphological classifications.     

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 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.     

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.     

Complete mitochondrial genome of the black‐tailed hornet,Vespa ducalis (Hymenoptera: Vespidae): Genomic comparisons in Vespoidea

We sequenced the complete mitochondrial genome (mitogenome) of the black‐tailed hornet, Vespa ducalis (Hymenoptera: Vespidae). The genome was 15,779‐bp long and contained typical sets of genes [13 protein‐coding genes (PCGs), 22 tRNAs, and 2 rRNAs]. The V. ducalis A + T‐rich region was 166‐bp long and was the shortest of all sequenced Vespoidea genomes, including Vespa. The genome was highly biased toward A/T nucleotides—80.1 % in the whole genome, 77.8 % in PCGs, 83.4–85.6 % in RNAs, and 92.8 % in the A + T‐rich region. These values are well within the typical range for genes and regions of Vespoidea mitogenomes. Start and stop codons in several Vespa species—including V. ducalis—were diversified, despite these species belonging to the same genus. In comparison with the ancestral mitogenomes, Vespa mitogenomes—including that of V. ducalis—showed substantial gene rearrangement; however, we detected no gene rearrangement among Vespa species. We conducted phylogenetic reconstruction based on concatenated sequences of 13 PCGs and two rRNAs (12,755 bp ) in available species of Vespoidea—21 species in six subfamilies in two families (Vespidae and Formicidae). The Bayesian inference and maximum likelihood (ML) methods revealed that each family formed strong monophyletic groups [Bayesian posterior probability (BPP) = 1; ML, 100 %]. Moreover, V. ducalis and V. mandarinia formed a strong sister group (BPP = 1; ML, 94 %).     

The complete mitochondrial genome of Diaphorina citri (Hemiptera: Psyllidae) and phylogenetic analysis

The Asian citrus psyllid, Diaphorina citri Kuwayama (Hemiptera: Psyllidae) is the most serious pest of citrus as the vector of Huanglongbing (HLB), the citrus greening disease. In this study, the complete mitochondrial genome (mitogenome) of D. citri has been sequenced and annotated, and a comparative analysis is provided with known Psylloidea species. The mitogenome of D. citri is a typical circular molecular of 15,038 bp in length with an A + T content of 74.56%, contains the typical 37 genes and the gene order is identical to the other Psylloidea mitogenomes. The nucleotide composition and codon usage of D. citri are similar to the four Psylloidea species. All protein-coding genes (PCGs) use standard initiation codons (TAN), stop with TAA and TAG except ND2 and ND5 which stop with incomplete termination codon T. All tRNAs have the typical clover-leaf structure, with the exception of trnS1 lacking the dihydrouridine (DHU) arm. The control region is located between rrnL and the trnI gene with the highest A + T content among the five Psylloidea species. Phylogenetic analysis is conducted based on the 13 PCGs or/and 2 rRNAs of 23 Sternorrhycha mitogenomes. Both maximum likelihood (ML) and Bayesian inference (BI) analysis suggest a clear relationship of Psylloidea, Aleyrodoidea and Aphidoidea within Sternorrhycha, which support the traditional morphological classification.     

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.     

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.     

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.     

Complete mitochondrial genome sequence of the Arctic gammarid,Onisimus nanseni (Crustacea; Amphipoda): Novel gene structures and unusual control region features

To analyze the mitogenome of the amphipod Onisimus nanseni, we amplified the complete mitogenome of O. nanseni using long-PCR and genome walking techniques. The mitogenome of O. nanseni is circular and contains all the typical mt genes (2 rRNAs, 22 tRNAs, and 13 protein-coding genes). It has two peculiar non-coding regions of 148 bp and 194 bp. The latter can be involved in replication and termination processes. The total length of the pooled protein-coding, rRNA, and tRNA genes is shorter than those of other crustaceans. In addition, the intergenic spacers of the O. nanseni mitogenome are considerably shorter in length than those of other crustaceans. Fourteen adjacent genes overlap, resulting in a compact mitogenomic structure. In the O. nanseni mitogenome, the AT composition is elevated, particularly in the control regions (78.9% AT), as has been demonstrated for two other amphipods. The tRNA order is highly rearranged compared to other arthropod mitogenomes, but the order of protein-coding genes and rRNAs is largely conserved. The gene cluster between the CO1 and CO3 genes is completely conserved among all amphipods compared. This provides insights into the evolution and gene structures of crustacean mitochondrial genomes, particularly in amphipods.     

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.     

The complete mitochondrial genomes of three grasshoppers, Asiotmethis zacharjini, Filchnerella helanshanensis and Pseudotmethis rubimarginis (Orthoptera: Pamphagidae)

The complete mitogenomes of Asiotmethis zacharjini, Filchnerella helanshanensis and Pseudotmethis rubimarginis are 15,660 bp, 15,657 bp and 15,661 bp in size, respectively. All three mitogenomes contain a standard set of 13 protein - coding genes, 22 transfer RNA genes (tRNAs), 2 ribosomal RNA genes (rRNAs) and an A + T-rich region in the same order as those of the other analysed caeliferan species, including the rearrangement of trnAsp and trnLys. The putative initiation codon for the cox1 gene in the three species is CCG. The long polythymine stretch (T-stretch) in the A + T-rich region of the three species is not adjacent to the trnIle but inside the stem–loop sequence in the majority strand. The mitogenomes of F. helanshanensis and P. rubimarginis have higher overall similarities. The characterization of the three mitogenomes will enrich our knowledge on the Pamphagidae mitogenome. The phylogenetic analyses indicated that within the Caelifera, Pyrgomorphoidea is a sister group to Acridoidea. The species from the Pamphagidae form a monophyletic group, as is the case for Acrididae. Furthermore, the two families cluster as sister groups, supporting the monophyly of Acridoidea. The relationships among eight acridid subfamilies were (Cyrtacanthacridinae + (Calliptaminae + (Catantopinae + (Oxyinae + (Melanopline + (Acridinae + (Oedipodinae + Gomphocerinae))))))).     

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.