The complete mitochondrial genome ofMacrotermes barneyiLight (Isoptera: Termitidae)

Novel Gene Rearrangements in Mitochondrial Genomes of four families of Praying Mantises (Insecta, Mantodea) and Phylogenetic Relationships of Mantodea

The mitochondrial genome (mitogenome) plays an important role in phylogenetic studies of many species. The mitogenomes of many praying mantis groups have been well-studied, but mitogenomes of special mimic praying mantises, especially Acanthopoidea and Galinthiadoidea species, are still sorely lacking in the NCBI database. The present study analyzes five mitogenomes from four species of Acanthopoidea (Angela sp., Callibia diana, Coptopteryx sp., Raptrix fusca) and one of Galinthiadoidea (Galinthias amoena) that were sequenced by the primer-walking method. Among Angela sp. and Coptopteryx sp., three gene rearrangements were detected in ND3-A-R-N-S-E-F and COX1-L2-COX2 gene regions, two of which were novel. In addition, individual tandem repeats were found in control regions of four mitogenomes (Angela sp., C. diana, Coptopteryx sp., G. amoena). For those, plausible explanations were derived from the tandem duplication-random loss (TDRL) model and the slipped-strand mispairing model. One potential motif was found in Acanthopidae that was seen as a synapomorphy. Several conserved block sequences (CBSs) were detected within Acanthopoidea that paved the way for the design of specific primers. Via BI and ML analysis, based on four datasets (PCG12, PCG12R, PCG123, PCG123R), the merged phylogenetic tree within Mantodea was reconstructed. This showed that the monophyly of Acanthopoidea was supported and that the PCG12R dataset was the most suitable for reconstructing the phylogenetic tree within Mantodea.

Characterization of the complete mitochondrial genome of Cryptotermes domesticus (Blattodea: Kalotermitidae): Genome description and phylogenetic implications

The complete mitochondrial genome of Cryptotermes domesticus (Haviland) was sequenced and annotated to study its characteristics and the phylogenetic relationship of C. domesticus to other termite species. The mitogenome of C. domesticus is a circular, close, and double-stranded molecule with a length of 15,655 bp. The sequenced mitogenome contains 37 typical genes, which are highly conserved in gene size, organization, and codon usage. Transfer RNA genes (tRNAs) also have typical secondary structures. All of the 13 protein-coding genes (PCGs) start with an ATN codon, except for nad4, which starts with GTG and terminates with the terminal codon TAA and TAG or the incomplete form T-- (cox2 and nad5). Most tRNAs have a typical cloverleaf structure, except for trnS1, in which this form is replaced by a simple loop and lacks the dihydrouridine (DHU) arm. The nucleotide diversity (Pi) and nonsynonymous (Ka)/synonymous (Ks) mutation rate ratios indicate that nad1, cox1, and cox3 are the most conserved genes, and that cox1 has the lowest rate of evolution. In addition, an 89 bp repeated sequence was found in the A + T-rich region. Phylogenetic analysis was performed using Bayesian inference (BI) and maximum likelihood (ML) methods based on 13 PCGs, and the monophyly of Kalotermitidae was supported.

Novel Mitochondrial Gene Rearrangement and Intergenic Regions Exist in the Mitochondrial Genomes from Four Newly Established Families of Praying Mantises (Insecta: Mantodea)

Simple Summary

Mantodea is regarded as an excellent material to study the gene rearrangements and large non-coding regions (LNCRs) in mitochondrial genomes. Meanwhile, as a result of the convergent evolution and parallelism, the gene rearrangements and LNCRs are specific to some taxonomic groups within Mantodea, which play an important role in phylogenetic relationship research. Nine mitochondrial genomes (mitogenomes) from four newly established families of praying mantises are obtained and annotated. Eight types of gene rearrangements, including four novel types of gene rearrangements in Mantodea, are detected, which can be explained by the tandem replication-random loss (TDRL) model. Moreover, one conserved motif between trnI-trnQ is detected in Toxoderidae. This study shed light on the formation mechanisms of these gene rearrangements and LNCRs in four newly established families of praying mantises.

Abstract

Long non-coding regions (NCRs) and gene rearrangements are commonly seen in mitochondrial genomes of Mantodea and are primarily focused on three regions: CR-I-Q-M-ND2, COX2-K-D-ATP8, and ND3-A-R-N-S-E-F-ND5. In this study, eight complete and one nearly complete mitochondrial genomes of praying mantises were acquired for the purpose of discussing mitochondrial gene rearrangements and phylogenetic relationships within Mantodea, primarily in the newly established families Haaniidae and Gonypetidae. Except for Heterochaeta sp. JZ-2017, novel mitochondrial gene arrangements were detected in Cheddikulama straminea, Sinomiopteryx graham, Pseudovates chlorophaea, Spilomantis occipitalis. Of note is the fact that one type of novel arrangement was detected for the first time in the Cyt b-S2-ND1 region. This could be reliably explained by the tandem replication-random loss (TDRL) model. The long NCR between trnT and trnP was generally found in Iridopteryginae and was similar to the ND4L or ND6 gene. Combined with gene rearrangements and intergenic regions, the monophyly of Haaniidae was supported, whereas the paraphyly of Gonypetidae was recovered. Furthermore, several synapomorphies unique to some clades were detected that conserved block sequences between trnI and trnQ and gaps between trnT and trnP in Toxoderidae and Iridopteryginae, respectively.

The complete mitochondrial genome of Leptomantella tonkinae (Hebard, 1920) (Mantodea: Leptomantellidae) and its phylogeny

The complete mitochondrial (mt) genome of Leptomantella tonkinae (Hebard, 1920) was 15,527 bp in length and contained 13 protein-coding genes, 22 transfer RNAs, two ribosomal RNAs, and one control region. The gene arrangement of mt genome of L. tonkinae was identical to the primitive mantis. The overall AT content of the mt genome was 74%. In ML and BI phylogenetic analyses, the monophyly of Leptomantellidae was robustly supported and the clade of Leptomantellidae is a sister clade to the group of (Gonypetidae+(Leptomantellidae+(Amorphoscelidae+Nanomantidae))).

The complete mitochondrial genome of Xanthomantis bimaculata (Mantodea: Iridopterygidae) and its phylogeny

The mitochondrial genome sequence of Xanthomantis bimaculata (Mantodea: Iridopterygidae) from Yunnan, China is a circular molecule with the typical insect mitochondrial gene arrangement, which is 15,941 bp in length and contains 22 tRNAs, two rRNAs, 13 protein-coding genes, and one control region. The overall AT content of the mitogenome is 73.12% (A = 37.58%, T = 35.54%, C = 16.54%, G = 10.34%). In BI and ML phylogenetic analyses, X. bimaculata was a sister clade to Sceptuchus simplex. The monophyly of the families Iridopterygidae, Thespidae and Liturgusidae were supported.

Data for praying mantis mitochondrial genomes and phylogenetic constructions within Mantodea

In this data article, we provide five datasets of mantis mitochondrial genomes: (1) PCG123: nucleotide sequences of 13 protein-coding genes including all codon positions; (2) PCG123R: nucleotide sequences of two rRNAs and 13 protein-coding genes including all codon positions; (3) PCG12: nucleotide sequences of 13 protein-coding genes without third codon positions; (4) PCG12R: nucleotide sequences of two rRNAs and 13 protein-coding genes without third codon positions, and (5) PCGAA: amino acid sequences of 13 protein-coding genes. These were used to construct phylogenetic relationships within Mantodea and the phylogenetic trees inferred from Bayesian analysis using two data sets (PCG12R, PCGAA) and Maximum Likelihood analysis using four data sets (PCG123, PCG12, PCG12R and PCGAA). We also provide initiation codon, termination codon, amino acid length and nucleotide diversity (Pi) of protein-coding genes among 27 mantises. The whole mitochondrial genomes of 27 praying mantises were submitted to GenBank with the accession numbers KY689112-KY689138.

Gene characteristics of the complete mitochondrial genomes of Paratoxodera polyacantha and Toxodera hauseri (Mantodea: Toxoderidae)

The family Toxoderidae (Mantodea) contains an ecologically diverse group of praying mantis species that have in common greatly elongated bodies. In this study, we sequenced and compared the complete mitochondrial genomes of two Toxoderidae species, Paratoxodera polyacantha and Toxodera hauseri, and compared their mitochondrial genome characteristics with another member of the Toxoderidae, Stenotoxodera porioni (KY689118). The lengths of the mitogenomes of T. hauseri and P. polyacantha were 15,616 bp and 15,999 bp, respectively, which is similar to that of S. porioni (15,846 bp). The size of each gene as well as the A+T-rich region and the A+T content of the whole genome were also very similar among the three species as were the protein-coding genes, the A+T content and the codon usages. The mitogenome of T. hauseri had the typical 22 tRNAs, whereas that of P. polyacantha had 26 tRNAs including an extra two copies of trnA-trnR. Intergenic regions of 67 bp and 76 bp were found in T. hauseri and P. polyacantha, respectively, between COX2 and trnK; these can be explained as residues of a tandem duplication/random loss of trnK and trnD. This non-coding region may be synapomorphic for Toxoderidae. In BI and ML analyses, the monophyly of Toxoderidae was supported and P. polyacantha was the sister clade to T. hauseri and S. porioni.

The Complete Mitochondrial Genome of Coptotermes 'suzhouensis' (syn. Coptotermes formosanus) (Isoptera: Rhinotermitidae) and Molecular Phylogeny Analysis

Coptotermes suzhouensis (Isoptera: Rhinotermitidae) is a significant subterranean termite pest of wooden structures and is widely distributed in southeastern China. The complete mitochondrial DNA sequence of C. suzhouensis was analyzed in this study. The mitogenome was a circular molecule of 15,764 bp in length, which contained 13 protein-coding genes (PCGs), 22 transfer RNA genes, two ribosomal RNA genes, and an A+T-rich region with a gene arrangement typical of Isoptera mitogenomes. All PCGs were initiated by ATN codons and terminated by complete termination codons (TAA), except COX2, ND5, and Cytb, which ended with an incomplete termination codon T. All tRNAs displayed a typical clover-leaf structure, except for tRNASer(AGN), which did not contain the stem-loop structure in the DHU arm. The A+T content (69.23%) of the A+T-rich region (949 bp) was higher than that of the entire mitogenome (65.60%), and two different sets of repeat units (A+B) were distributed in this region. Comparison of complete mitogenome sequences with those of Coptotermes formosanus indicated that the two taxa have very high genetic similarity. Forty-one representative termite species were used to construct phylogenetic trees by maximum likelihood, maximum parsimony, and Bayesian inference methods. The phylogenetic analyses also strongly supported (BPP, MLBP, and MPBP = 100%) that all C. suzhouensis and C. formosanus samples gathered into one clade with genetic distances between 0.000 and 0.002. This study provides molecular evidence for a more robust phylogenetic position of C. suzhouensis and inferrs that C. suzhouensis was the synonymy of C. formosanus.

Complete mitochondrial genome of Cacopsylla coccinae (Hemiptera: Psyllidae)

In this study, the first complete mitochondrial genome (mitogenome) sequence of Cacopsylla coccinae was determined by long PCR and primer walking methods. The complete mitochondrial genome is 14,832 bp in length and contains 13 protein-coding genes, 2 ribosomal RNA genes, 22 transfer RNA genes as well as a control region. The overall base composition of the genome is A (38.16%), T (33.88%), C (17.95%) and G (10.01%). Stop codon was incomplete for coxII gene and ND1 gene. The gene overlaps were suggested between 13 pairs of the contiguous genes in C. coccinae. The mitogenome would contribute to resolving phylogenetic position and interrelationships of Cacopsylla.

The complete mitochondrial genomes of three bristletails (Insecta: Archaeognatha): the paraphyly of Machilidae and insights into archaeognathan phylogeny

The order Archaeognatha was an ancient group of Hexapoda and was considered as the most primitive of living insects. Two extant families (Meinertellidae and Machilidae) consisted of approximately 500 species. This study determined 3 complete mitochondrial genomes and 2 nearly complete mitochondrial genome sequences of the bristletail. The size of the 5 mitochondrial genome sequences of bristletail were relatively modest, containing 13 protein-coding genes (PCGs), 2 ribosomal RNA (rRNA) genes, 22 transfer RNA (tRNA) genes and one control region. The gene orders were identical to that of Drosophila yakuba and most bristletail species suggesting a conserved genome evolution within the Archaeognatha. In order to estimate archaeognathan evolutionary relationships, phylogenetic analyses were conducted using concatenated nucleotide sequences of 13 protein-coding genes, with four different computational algorithms (NJ, MP, ML and BI). Based on the results, the monophyly of the family Machilidae was challenged by both datasets (W12 and G12 datasets). The relationships among archaeognathan subfamilies seemed to be tangled and the subfamily Machilinae was also believed to be a paraphyletic group in our study.

The evolutionary history of termites as inferred from 66 mitochondrial genomes

Termites have colonized many habitats and are among the most abundant animals in tropical ecosystems, which they modify considerably through their actions. The timing of their rise in abundance and of the dispersal events that gave rise to modern termite lineages is not well understood. To shed light on termite origins and diversification, we sequenced the mitochondrial genome of 48 termite species and combined them with 18 previously sequenced termite mitochondrial genomes for phylogenetic and molecular clock analyses using multiple fossil calibrations. The 66 genomes represent most major clades of termites. Unlike previous phylogenetic studies based on fewer molecular data, our phylogenetic tree is fully resolved for the lower termites. The phylogenetic positions of Macrotermitinae and Apicotermitinae are also resolved as the basal groups in the higher termites, but in the crown termitid groups, including Termitinae + Syntermitinae + Nasutitermitinae + Cubitermitinae, the position of some nodes remains uncertain. Our molecular clock tree indicates that the lineages leading to termites and Cryptocercus roaches diverged 170 Ma (153-196 Ma 95% confidence interval [CI]), that modern Termitidae arose 54 Ma (46-66 Ma 95% CI), and that the crown termitid group arose 40 Ma (35-49 Ma 95% CI). This indicates that the distribution of basal termite clades was influenced by the final stages of the breakup of Pangaea. Our inference of ancestral geographic ranges shows that the Termitidae, which includes more than 75% of extant termite species, most likely originated in Africa or Asia, and acquired their pantropical distribution after a series of dispersal and subsequent diversification events.

The complete mitochondrial genome of fungus-growing termite, Macrotermes natalensis (Isoptera: Macrotermitinae)

In this study, the complete mitochondrial genome of a fungus-growing termite, Macrotermes natalensis, with the total length of 16,325 bp, is reported for the first time. This mtgenome harbors 13 protein-coding genes, 22 transfer RNA genes, 2 ribosomal RNA genes, and 1 control region (D-loop). The total base composition is 44.1% of A, 21.5% of T, 11.5% of G, and 22.9% of C, so the percentage of A and T (65.6%) is much higher than that of G and C. Most of the genes are distributed on H-strand, except for four subunit genes (ND1, ND4, ND4L, ND5) and eight tRNA genes. All protein-coding genes start with an ATN codon, and terminate with the canonical stop codon (TAA/TAG) or a single T (T- -). The complete mitochondrial genome sequence reported here adds a new genetic resource for the genus Macrotermes and might be useful for phylogenetic and systematic analyses within the genus Macrotermes.

The complete mitogenome of the dampwood termite Zootermopsis nevadensis (Insecta: Isoptera: Termopsidae)

The complete mitogenome of the dampwood termite Zootermopsis nevadensis was reconstructed from whole-genome Illumina sequencing data with an average coverage of 7052×. The circular genome is 15,444 bp in length, and consists of 22 transfer RNAs (tRNAs), 13 protein-coding genes (PCGs), two ribosomal RNAs (rRNAs) and one D-loop region. All PCGs are initiated with ATN codons, except for the ND1 and ND5 genes with the start codon GTG. Some PCGs harbor TAG (ND1) or incomplete stop codon T (COII, ND5 & Cytb), while the others use TAA as their stop codons. The nucleotide composition is asymmetric (45.8% A, 19.8% C, 10.9% G, 23.5% T) with an overall GC content of 30.7%. These data would contribute to the design of novel molecular markers for population and evolutionary studies of this and related termite species.