Search

Potamanthidae belongs to the superfamily Ephemeroidea but has no complete mt genome released in the NCBI (except for two unchecked and one partial mt genome). Since the sister clade to Potamanthidae has always been controversial, we sequenced seven mt genomes of Potamanthidae (two species from Rhoenanthus and five species from Potamanthus) in order to rebuild the phylogenetic relationships of Potamanthidae in this study. The divergence time of Potamanthidae was also investigated by utilizing five fossil calibration points because of the indeterminate origin time. In addition, because Rhoenanthus coreanus and Potamanthus luteus are always in low-temperature environments, we aimed to explore whether these two species were under positive selection at the mt genome level. Amongst the 13 PCGs, CGA was used as the start codon in COX1, whereas other genes conformed to initiating with an ATN start codon. From this analysis, UUA (L), AUU (I), and UUU (F) had the highest usage. Furthermore, the DHU arm was absent in the secondary structure of S1 in all species. By combining the 13 PCGs and 2 rRNAs, we reconstructed the phylogenetic relationship of Potamanthidae within Ephemeroptera. The monophyly of Potamanthidae and the monophyly of Rhoenanthus and Potamanthus were supported in the results. The phylogenetic relationship of Potamanthidae + (Ephemeridae + Polymitarcyidae) was also recovered with a high prior probability. The divergence times of Potamanthidae were traced to be 90.44 Mya (95% HPD, 62.80–121.74 Mya), and the divergence times of Rhoenanthus and Potamanthus originated at approximately 64.77 Mya (95% HPD, 43.82–88.68 Mya), thus belonging to the late Pliocene Epoch or early Miocene Epoch. In addition, the data indicated that R. coreanus was under negative selection and that ATP8 and ND2 in Potamanthidae had a high evolutionary rate.

Unusual climates can lead to extreme temperatures. Fejervarya kawamurai, one of the most prevalent anurans in the paddy fields of tropical and subtropical regions in Asia, is sensitive to climate change. The present study focuses primarily on a single question: how do the 13 mitochondrial protein-coding genes (PCGs) respond to extreme temperature change compared with 25 °C controls? Thirty-eight genes including an extra tRNA-Met gene were identified and sequenced from the mitochondrial genome of F. kawamurai. Evolutionary relationships were assessed within the Dicroglossidae and showed that Dicroglossinae is monophyletic and F. kawamurai is a sister group to the clade of (F. multistriata + F. limnocharis). Transcript levels of mitochondrial genes in liver were also evaluated to assess responses to 24 h exposure to low (2 °C and 4 °C) or high (40 °C) temperatures. Under 2 °C, seven genes showed significant changes in liver transcript levels, among which transcript levels of ATP8, ND1, ND2, ND3, ND4, and Cytb increased, respectively, and ND5 decreased. However, exposure to 4 °C for 24 h was very different in that the expressions of ten mitochondrial protein-coding genes, except ND1, ND3, and Cytb, were significantly downregulated. Among them, the transcript level of ND5 was most significantly downregulated, decreasing by 0.28-fold. Exposure to a hot environment at 40 °C for 24 h resulted in a marked difference in transcript responses with strong upregulation of eight genes, ranging from a 1.52-fold increase in ND4L to a 2.18-fold rise in Cytb transcript levels, although COI and ND5 were reduced to 0.56 and 0.67, respectively, compared with the controls. Overall, these results suggest that at 4 °C, F. kawamurai appears to have entered a hypometabolic state of hibernation, whereas its mitochondrial oxidative phosphorylation was affected at both 2 °C and 40 °C. The majority of mitochondrial PCGs exhibited substantial changes at all three temperatures, indicating that frogs such as F. kawamurai that inhabit tropical or subtropical regions are susceptible to ambient temperature changes and can quickly employ compensating adjustments to proteins involved in the mitochondrial electron transport chain.

The phylogenetic relationship of Ephemeridae (Insect: Ephemeroptera) remains hotly debated using mitochondrial (mt) genomes. All previously reported mt genomes of Ephemeridae belong to the genus Ephemera. This study provides the first complete mt genome sequence from the genus Hexagenia with an analysis of the mitogenome of Hexagenia rigida Mc Dunnough, 1924 (Ephemeroptera: Ephemeridae) and providing new information to discuss the phylogenetic relationships within Ephemeroptera. The complete mt genome of H. rigida was a circular molecule of 16,159 bp in length, containing 37 genes (2 rRNA genes, 13 protein-coding genes, 22 tRNA genes), which showed the typical mt gene arrangement of insects. The AT content of the whole genome was 70.0% and the length of the control region was 1091 bp. All protein-coding genes used ATN as the start codon, and most PCGs used TAA/TAG as the stop codons excluding COI, COII, ND5 and Cyt b that used T as the stop codon. BI and ML phylogenetic trees constructed from 27 species of 13 families showed that Ephemeridae is a sister clade to the clade Polymitarcyidae.

Zoodes fulguratus is a common species of Cerambycidae, reported from Vietnam, Nepal, Laos, Burma, and China. To date, no mitochondrial genomes of the genus Zoodes have been reported. In this study, we sequenced and analyzed the first mitochondrial genome of Z. fulguratus to discuss its phylogenetic relationship within the subfamily Cerambycinae. This mitochondrial genome showed the typical insect gene arrangement: a circular molecule of 15,885 bp long with 13 protein-coding genes, two ribosomal RNA genes (rRNAs), and 22 transfer RNA genes (tRNAs). The AT content of the whole mitogenome was 74.2% with a high asymmetric nucleotide presentation revealed by a positive AT-skew (0.067) and a negative GC-skew (–0.178), whereas the AT content of the A + T rich region was 80%. The Maximum likelihood (ML) and Bayesian inference (BI) phylogenetic analyses showed that Z. fulguratus is a sister clade of Gnatholea eburifera.

Background The phylogenetic relationships of Odonata (dragonflies and damselflies) and Ephemeroptera (mayflies) remain unresolved. Different researchers have supported one of three hypotheses (Palaeoptera, Chiastomyaria or Metapterygota) based on data from different morphological characters and molecular markers, sometimes even re-assessing the same transcriptomes or mitochondrial genomes. The appropriate choice of outgroups and more taxon sampling is thought to eliminate artificial phylogenetic relationships and obtain an accurate phylogeny. Hence, in the current study, we sequenced 28 mt genomes from Ephemeroptera, Odonata and Plecoptera to further investigate phylogenetic relationships, the probability of each of the three hypotheses, and to examine mt gene arrangements in these species. We selected three different combinations of outgroups to analyze how outgroup choice affected the phylogenetic relationships of Odonata and Ephemeroptera. Methods Mitochondrial genomes from 28 species of mayflies, dragonflies, damselflies and stoneflies were sequenced. We used Bayesian inference (BI) and Maximum likelihood (ML) analyses for each dataset to reconstruct an accurate phylogeny of these winged insect orders. The effect of outgroup choice was assessed by separate analyses using three outgroups combinations: (a) four bristletails and three silverfish as outgroups, (b) five bristletails and three silverfish as outgroups, or (c) five diplurans as outgroups. Results Among these sequenced mitogenomes we found the gene arrangement IMQM in Heptageniidae (Ephemeroptera), and an inverted and translocated tRNA-Ile between the 12S RNA gene and the control region in Ephemerellidae (Ephemeroptera). The IMQM gene arrangement in Heptageniidae (Ephemeroptera) can be explained via the tandem-duplication and random loss model, and the transposition and inversion of tRNA-Ile genes in Ephemerellidae can be explained through the recombination and tandem duplication-random loss (TDRL) model. Our phylogenetic analysis strongly supported the Chiastomyaria hypothesis in three different outgroup combinations in BI analyses. The results also show that suitable outgroups are very important to determining phylogenetic relationships in the rapid evolution of insects especially among Ephemeroptera and Odonata. The mt genome is a suitable marker to investigate the phylogeny of inter-order and inter-family relationships of insects but outgroup choice is very important for deriving these relationships among winged insects. Hence, we must carefully choose the correct outgroup in order to discuss the relationships of Ephemeroptera and Odonata.

Simple Summary: As one of the burrowing mayfly groups with large mandibular tusks, the phylogenetic relationships within Potamanthidae are always controversial. There are at least two opposite hypotheses for mayfly grouping: Potamanthidae + (Ephemeridae + Polymitarcyidae) and (Potamanthidae + Ephemeridae) + Polymitarcyidae. Because of the indeterminate origin time of this group, the present study aimed to reconstruct the phylogenetic relationship and explore the origin time of Potamanthidae based on mitochondrial (mt) genomes. In addition, the protein-coding genes (PCGs) of these mt genomes may undergo positive selection when these species live in low-temperature environments. Potamanthidae belongs to the superfamily Ephemeroidea but has no complete mt genome released in the NCBI (except for two unchecked and one partial mt genome). Since the sister clade to Potamanthidae has always been controversial, we sequenced seven mt genomes of Potamanthidae (two species from Rhoenanthus and five species from Potamanthus) in order to rebuild the phylogenetic relationships of Potamanthidae in this study. The divergence time of Potamanthidae was also investigated by utilizing five fossil calibration points because of the indeterminate origin time. In addition, because Rhoenanthus coreanus and Potamanthus luteus are always in low-temperature environments, we aimed to explore whether these two species were under positive selection at the mt genome level. Amongst the 13 PCGs, CGA was used as the start codon in COX1, whereas other genes conformed to initiating with an ATN start codon. From this analysis, UUA (L), AUU (I), and UUU (F) had the highest usage. Furthermore, the DHU arm was absent in the secondary structure of S1 in all species. By combining the 13 PCGs and 2 rRNAs, we reconstructed the phylogenetic relationship of Potamanthidae within Ephemeroptera. The monophyly of Potamanthidae and the monophyly of Rhoenanthus and Potamanthus were supported in the results. The phylogenetic relationship of Potamanthidae + (Ephemeridae + Polymitarcyidae) was also recovered with a high prior probability. The divergence times of Potamanthidae were traced to be 90.44 Mya (95% HPD, 62.80–121.74 Mya), and the divergence times of Rhoenanthus and Potamanthus originated at approximately 64.77 Mya (95% HPD, 43.82–88.68 Mya), thus belonging to the late Pliocene Epoch or early Miocene Epoch. In addition, the data indicated that R. coreanus was under negative selection and that ATP8 and ND2 in Potamanthidae had a high evolutionary rate. [ABSTRACT FROM AUTHOR]

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.

Ephemeroptera (Insecta: Pterygota) are widely distributed all over the world with more than 3500 species. During the last decade, the phylogenetic relationships within Ephemeroptera have been a hot topic of research, especially regarding the phylogenetic relationships among Vietnamellidae. In this study, three mitochondrial genomes from three populations of Vienamella sinensis collected from Tonglu (V. sinensis TL), Chun’an (V. sinensis CN), and Qingyuan (V. sinensis QY) in Zhejiang Province, China were compared to discuss the potential existence of cryptic species. We also established their phylogenetic relationship by combining the mt genomes of 69 Ephemeroptera downloaded from NCBI. The mt genomes of V. sinensis TL, V. sinensis CN, and V. sinensis QY showed the same gene arrangement with lengths of 15,674 bp, 15,674 bp, and 15,610 bp, respectively. Comprehensive analyses of these three mt genomes revealed significant differences in mt genome organization, genetic distance, and divergence time. Our results showed that the specimens collected from Chun’an and Tonglu in Zhejiang Province, China belonged to V. sinensis, and the specimens collected from Qingyuan, Zhejiang Province, China were a cryptic species of V. sinensis. In maximum likelihood (ML) and Bayesian inference (BI) phylogenetic trees, the monophyly of the family Vietnamellidae was supported and Vietnamellidae has a close relationship with Ephemerellidae.

As a small order of Pterygota (Insecta), Ephemeroptera has almost 3,500 species around the world. Ephemerellidae is a widely distributed common group of Ephemeroptera. However, the relationship among Ephemerellidae, Vietnamellidae and Teloganellidae is still in dispute. In this study, we sequenced six complete mitogenomes of three genera from Ephemerellidae (Insecta: Ephemeroptera): Ephemerella sp. Yunnan-2018, Serratella zapekinae, Serratella sp. Yunnan-2018, Serratella sp. Liaoning-2019, Torleya grandipennis and T. tumiforceps. These mitogenomes were employed to reveal controversial phylogenetic relationships among the Ephemeroptera, with emphasis on the phylogenetic relationships among Ephemerellidae. The lengths of the six mayfly mitogenomes ranged from 15,134 bp to 15,703 bp. Four mitogenomes of Ephemerella sp. Yunnan-2018, Serratella zapekinae, Serratella sp. Yunnan-2018 and Serratella sp. Liaoning-2019 had 22 tRNAs including an inversion and translocation of trnI. By contrast, the mitogenomes of T. tumiforceps and T. grandipennis had 24 tRNAs due to an extra two copies of inversion and translocation of trnI. Within the family Ephemerellidae, disparate gene rearrangement occurred in the mitogenomes of different genera: one copy of inversion and translocation trnI in the genera Ephemerella and Serratella, and three repeat copies of inversion and translocation of trnI in the genus Torleya. A large non-coding region (≥200 bp) between trnS1 (AGN) and trnE was detected in T. grandipennis and T. tumiforceps. Among the phylogenetic relationship of the Ephemeroptera, the monophyly of almost all families except Siphlonuridae was supported by BI and ML analyses. The phylogenetic results indicated that Ephemerellidae was the sister clade to Vietnamellidae whereas Teloganellidae was not a sister clade of Ephemerellidae and Vietnamellidae.

Abstract Background Complete mitochondrial (mt) genomes have been used extensively to test hypotheses about microevolution and to study population structure, phylogeography, and phylogenetic relationships of Anura at various taxonomic levels. Large-scale mt genomic reorganizations have been observed among many fork-tongued frogs (family Dicroglossidae). The relationships among Dicroglossidae and validation of the genus Feirana are still problematic. Hence, we sequenced the complete mt genomes of Nanorana taihangnica (=F. taihangnica) and N. yunnanensis as well as partial mt genomes of six Quasipaa species (dicroglossid taxa), two Odorrana and two Amolops species (Ranidae), and one Rhacophorus species (Rhacophoridae) in order to identify unknown mt gene rearrangements, to investigate the validity of the genus Feirana, and to test the phylogenetic relationship of Dicroglossidae. Results In the mt genome of N. taihangnica two trnM genes, two trnP genes and two control regions were found. In addition, the trnA, trnN, trnC, and trnQ genes were translocated from their typical positions. In the mt genome of N. yunnanensis, three control regions were found and eight genes (ND6, trnP, trnQ, trnA, trnN, trnC, trnY and trnS genes) in the L-stand were translocated from their typical position and grouped together. We also found intraspecific rearrangement of the mitochondrial genomes in N. taihangnica and Quasipaa boulengeri. In phylogenetic trees, the genus Feirana nested deeply within the clade of genus Nanorana, indicating that the genus Feirana may be a synonym to Nanorana. Ranidae as a sister clade to Dicroglossidae and the clade of (Ranidae + Dicroglossidae) as a sister clade to (Mantellidae + Rhacophoridae) were well supported in BI analysis but low bootstrap in ML analysis. Conclusions We found that the gene arrangements of N. taihangnica and N. yunnanensis differed from other published dicroglossid mt genomes. The gene arrangements in N. taihangnica and N. yunnanensis could be explained by the Tandem Duplication and Random Loss (TDRL) and the Dimer-Mitogenome and Non-Random Loss (DMNR) models, respectively. The invalidation of the genus Feirana is supported in this study.

The relationship of Matrona and Atrocalopteryx (Odonata: Calopterygidae) is still unclear. To better understand the phylogenetic relationship of Matrona and Atrocalopteryx, we sequenced and annotated two complete mitochondrial genomes of Matrona basilaris sampled from two different locations. The length of the two complete mitochondrial genomes of M. basilaris is 16,149 bp and 15,893 bp for the specimens collected in Jinhua, Zhejiang Province and Tianmushan, Zhejiang Province, China, respectively. The two mitochondrial genomes include the typical invertebrate set of 37 genes: 13 protein-coding genes (PGCs), 22 tRNA genes, and 2 rRNA genes. The nucleotide composition of the mitogenome is similar to other odonates with high content of A + T (68.9%) and all PCGs use ATN as the start codon. Tandem repeats were detected in the control regions of the two M. basilaris samples that accounted for the different sequence lengths of the mitochondrial genomes from the two locations. Finally, BI and ML phylogenetic analysis based on the concatenated nucleotide sequences of the 13 PCGs supported the conclusion that M. basilaris is a sister clade to Atrocalopteryx melli.

The complete mitochondrial genome of Xystrocera globosa is 15,706 bp in length, containing 13 protein-coding genes, 22 transfer RNAs, two ribosomal RNAs and the A + T-rich region. The overall base composition is 72.7% AT and 27.3% GC, and the AT content of the control region is 79.3%. In ML and BI phylogenetic trees, X. globosa was a sister clade to X. grayii. The monophyly of Lamiinae and Prioninae were supported in ML analyses, but nevertheless, the monophyly of Cerambycinae was not recovered.

The availability of next-generation sequencing (NGS) in recent years has facilitated a revolution in the availability of mitochondrial (mt) genome sequences. The mt genome is a powerful tool for comparative studies and resolving the phylogenetic relationships among insect lineages. The mt genomes of phytophagous scarabs of the subfamilies Cetoniinae and Dynastinae were under-represented in GenBank. Previous research found that the subfamily Rutelinae was recovered as a paraphyletic group because the few representatives of the subfamily Dynastinae clustered into Rutelinae, but the subfamily position of Dynastinae was still unclear. In the present study, we sequenced 18 mt genomes from Dynastinae and Cetoniinae using next-generation sequencing (NGS) to re-assess the phylogenetic relationships within Scarabaeidae. All sequenced mt genomes contained 37 sets of genes (13 protein-coding genes, 22 tRNAs, and two ribosomal RNAs), with one long control region, but the gene order was not the same between Cetoniinae and Dynastinae species. All mt genomes of Dynastinae species showed the same gene rearrangement of trnQ-NCR-trnI-trnM, whereas all mt genomes of Cetoniinae species showed the ancestral insect gene order of trnI-trnQ-trnM. Phylogenetic analyses (IQ-tree and MrBayes) were conducted using 13 protein-coding genes based on nucleotide and amino acid datasets. In the ML and BI trees, we recovered the monophyly of Rutelinae, Cetoniinae, Dynastinae, and Sericinae, and the non-monophyly of Melolonthinae. Cetoniinae was shown to be a sister clade to (Dynastinae + Rutelinae).

The order Microcoryphia, commonly known as bristletails, is considered as the most primitive one among living insects. Within this order, two species, Coreamachilis coreanus and C. songi (Machilidae: Machilinae), display the following contrasting reproductive strategies: parthenogenesis occurs in C. coreanus, whereas sexual reproduction is found in C. songi. In the present study, the complete mitogenomes of C. coreanus and C. songi were sequenced to compare their mitogenome structure, analyze relationships within the Microcoryphia, and assess adaptive evolution. The length of the mitogenomes of C. coreanus and C. songi were 15,578 bp and 15,570 bp, respectively, and the gene orders were those of typical insects. A long hairpin structure was found between the ND1 and 16S rRNA genes of both species that seem to be characteristic of Machilinae and Petrobiinae species. Phylogenetic assessment of Coreamachilis was conducted using BI and ML analyses with concatenated nucleotide sequences of the 13 protein-coding genes. The results showed that the monophyly of Machilidae, Machilinae, and Petrobiinae was not supported. The genus Coreamachilis (C. coreanus and C. songi) was a sister clade to Allopsontus helanensis, and then the clade of ((C. coreanus + C. songi) + A. helanensis) was a sister clade to A. baii, which suggests that the monophyly of Allopsontus was not supported. Positive selection analysis of the 13 protein-coding genes failed to reveal any positive selection in C. coreanus or C. songi. The long hairpin structures found in Machilinae and Petrobiinae were highly consistent with the phylogenetic results and could potentially be used as an additional molecular characteristic to further discuss relationships within the Microcoryphia.

We determined 15 complete and two nearly complete mitogenomes of Heptageniidae belonging to three subfamilies (Heptageniinae, Rhithrogeninae, and Ecdyonurinae) and six genera (Afronurus, Epeorus, Leucrocuta, Maccaffertium, Stenacron, and Stenonema). Species of Rhithrogeninae and Ecdyonurinae had the same gene rearrangement of CR-I-M-Q-M-ND2, whereas a novel gene rearrangement of CR-I-M-Q-NCR-ND2 was found in Heptageniinae. Non-coding regions (NCRs) of 25–47 bp located between trnA and trnR were observed in all mayflies of Heptageniidae, which may be a synapomorphy for Heptageniidae. Both the BI and ML phylogenetic analyses supported the monophyly of Heptageniidae and its subfamilies (Heptageniinae, Rhithrogeninae, and Ecdyonurinae). The phylogenetic results combined with gene rearrangements and NCR locations confirmed the relationship of the subfamilies as (Heptageniinae + (Rhithrogeninae + Ecdyonurinae)). To assess the effects of low-temperature stress on Heptageniidae species from Ottawa, Canada, we found 27 positive selection sites in eight protein-coding genes (PCGs) using the branch-site model. The selection pressure analyses suggested that mitochondrial PCGs underwent positive selection to meet the energy requirements under low-temperature stress.

The complete mitochondrial genome of Choroterpides apiculata (Ephemeroptera: Leptophlebiidae) is typically a circular molecule of 15,199 bp in length, containing 37 genes (13 protein-coding genes, 22 tRNAs, and two rRNAs) and one control region. The overall A + T content of the whole genome is 74% and the A + T content of the control region is 79.7%. In Bayesian inference and maximum-likelihood phylogenetic trees using 24 species from 13 families of Ephemeroptera, the monophyly of the families Isonychiidae, Heptageniidae, Vietnamellidae, Leptophlebiidae, Caenidae, and Baetidae were highly supported and C. apiculata was a sister group to Habrophlebiodes zijinensis.

The phylogenetic relationship of Caenidae remains hotly debated within the Ephemeroptera. We sequenced the complete mitochondrial genome of Caenis sp. (Ephemeroptera: Caenidae) to discuss the phylogenetic relationships among the Caenidae. The mitochondrial genome of Caenis sp. collected from Jian’ou, Fujian province, China is a circular molecule of 15,392 bp in length containing 37 genes (13 protein-coding genes, 22 tRNAs, and two rRNAs), which showed the typical insect mitochondrial gene arrangement. In BI and ML phylogenetic trees using 23 species from 13 families, the monophyly of the families Caenidae, Heptageniidae, Isonychiidae, and Vietnamellidae was strongly supported. The clade of Caenidae is a sister clade to the clade of Teloganodidae and Baetidae.

The mitochondrial genome of Fejervarya kawamurai is a circular molecule of 17,650 bp in length, containing 13 protein-coding genes, two rRNA genes, 23 tRNA genes (including an extra tRNA-Met), and the control region. The AT content of the whole genome is 56.9%. In Bayesian inference (BI) and Maximum likelihood (ML) analyses, we found that F. kawamurai is a sister clade to F. multistriata and F. limnocharis. The monophyly of Fejervarya, Quasipaa, Nanorana was well supported (1.00 in BI and 100% in ML).

The phylogenetic relationship between Ephemeroptera (mayflies) and Odonata (dragonflies and damselflies) remains hotly debated in the insect evolution community. We sequenced the complete mitochondrial genome of Caenis sp. (Ephemeroptera: Caenidae) to discuss the phylogenetic relationship of Palaeoptera. The mitochondrial genome of Caenis sp. is a circular molecule of 15,254 bp in length containing 37 genes (13 protein-coding genes, 22 tRNAs, and 2 rRNAs), which showed the typical insect mitochondrial gene arrangement. In BI and ML phylogenetic trees using 71 species of 12 orders, our results support the Ephemeroptera as the basal group of winged insects.

Cerambycidae is one of the most diversified groups within Coleoptera and includes nearly 35,000 known species. The relationships at the subfamily level within Cerambycidae have not been convincingly demonstrated and the gene rearrangement of mitochondrial genomes in Cerambycidae remains unclear due to the low numbers of sequenced mitogenomes. In the present study, we determined five complete mitogenomes of Cerambycidae and investigated the phylogenetic relationship among the subfamilies of Cerambycidae based on mitogenomes. The mitogenomic arrangement of all five species was identical to the ancestral Cerambycidae type without gene rearrangement. Remarkably, however, two large intergenic spacers were detected in the mitogenome of Pterolophia sp. ZJY-2019. The origins of these intergenic spacers could be explained by the slipped-strand mispairing and duplication/random loss models. A conserved motif was found between trnS2 and nad1 gene, which was proposed to be a binding site of a transcription termination peptide. Also, tandem repeat units were identified in the A + T-rich region of all five mitogenomes. The monophyly of Lamiinae and Prioninae was strongly supported by both MrBayes and RAxML analyses based on nucleotide datasets, whereas the Cerambycinae and Lepturinae were recovered as non-monophyletic.

We determined the mitochondrial gene sequence of Monochamus alternatus and three other mitogenomes of Lamiinae (Insect: Coleoptera: Cerambycidae) belonging to three genera (Aulaconotus, Apriona and Paraglenea) to enrich the mitochondrial genome database of Lamiinae and further explore the phylogenetic relationships within the subfamily. Phylogenetic trees of the Lamiinae were built using the Bayesian inference (BI) and maximum likelihood (ML) methods and the monophyly of Monochamus, Anoplophora, and Batocera genera was supported. Anoplophora chinensis, An. glabripennis and Aristobia reticulator were closely related, suggesting they may also be potential vectors for the transmission of the pine wood pathogenic nematode (Bursaphelenchus xylophilus) in addition to M. alternatus, a well-known vector of pine wilt disease. There is a special symbiotic relationship between M. alternatus and Bursaphelenchus xylophilus. As the native sympatric sibling species of B. xylophilus, B. mucronatus also has a specific relationship that is often overlooked. The analysis of mitochondrial gene expression aimed to explore the effect of B. mucronatus on the energy metabolism of the respiratory chain of M. alternatus adults. Using RT-qPCR, we determined and analyzed the expression of eight mitochondrial protein-coding genes (COI, COII, COIII, ND1, ND4, ND5, ATP6, and Cty b) between M. alternatus infected by B. mucronatus and M. alternatus without the nematode. Expression of all the eight mitochondrial genes were up-regulated, particularly the ND4 and ND5 gene, which were up-regulated by 4–5-fold (p < 0.01). Since longicorn beetles have immune responses to nematodes, we believe that their relationship should not be viewed as symbiotic, but classed as parasitic.

Mantodea is a predatory insect group, its members occupying a diverse array of widely distributed habitats. Praying mantis species utilize hunting strategies including remarkable mimicry and unique camouflage for hiding from natural enemies while catching their prey. The emergence of a "cyclopean ear" in mantises is thought to be a morphological innovation of the group, and an "arms race" with echolocating bats is one of the hypotheses put forward to account for the emergence of the mantis ear from a coevolutionary perspective. However, this hypothesis has not been rigorously tested because of a lack of robust higher‐level phylogeny and a detailed chronogram of Mantodea. Previous phylogenetic studies found an incongruence between traditional classification and molecular phylogenetics due to the convergent evolution of various ecomorphic strategies of the lineage. Here, we performed a comprehensive phylogenetic analysis of Mantodea based on data from 61 mitogenomes. Our analyses showed that the monophyly of Acanthopidae, Haaniidae, Nanomantidae, Miomantidae and Mantidea was supported. The newly updated Gonypetidae were paraphyletic, whereas Eremiaphilidae, Deroplatyidae and Toxoderidae were polyphyletic. Our molecular dating analyses inferred that Spinomantodea originated at ca. 149 Ma (Late Jurassic), whereas the origin of hearing mantises (Cernomantodea) was inferred as Early Cretaceous (119 Ma, 95% CI: 110–129 Ma). The molecular dating results indicated that the hearing organ in mantises did not arise in response to bat predation. Our study provides a robust framework for further evolutionary comparative studies of mantises. [ABSTRACT FROM AUTHOR]

Simple Summary: The tropical and subtropical paddy frog, Fejervarya kawamurai (Anura: Dicroglossidae), is a common and widespread species in Asia. Amphibians can adapt to small changes in climate, but extreme temperature changes can lead to metabolic abnormalities affecting numerous cell functions. Damage or disruption of mitochondrial respiratory chain complexes can be lethal. The present study characterizes the mitochondrial genome of F. kawamurai to evaluate its evolutionary relationship within Dicroglossidae and to analyze the effects of extreme temperature change on mitochondrial gene expression of proteins involved in oxidative phosphorylation. Unusual climates can lead to extreme temperatures. Fejervarya kawamurai, one of the most prevalent anurans in the paddy fields of tropical and subtropical regions in Asia, is sensitive to climate change. The present study focuses primarily on a single question: how do the 13 mitochondrial protein-coding genes (PCGs) respond to extreme temperature change compared with 25 °C controls? Thirty-eight genes including an extra tRNA-Met gene were identified and sequenced from the mitochondrial genome of F. kawamurai. Evolutionary relationships were assessed within the Dicroglossidae and showed that Dicroglossinae is monophyletic and F. kawamurai is a sister group to the clade of (F. multistriata + F. limnocharis). Transcript levels of mitochondrial genes in liver were also evaluated to assess responses to 24 h exposure to low (2 °C and 4 °C) or high (40 °C) temperatures. Under 2 °C, seven genes showed significant changes in liver transcript levels, among which transcript levels of ATP8, ND1, ND2, ND3, ND4, and Cytb increased, respectively, and ND5 decreased. However, exposure to 4 °C for 24 h was very different in that the expressions of ten mitochondrial protein-coding genes, except ND1, ND3, and Cytb, were significantly downregulated. Among them, the transcript level of ND5 was most significantly downregulated, decreasing by 0.28-fold. Exposure to a hot environment at 40 °C for 24 h resulted in a marked difference in transcript responses with strong upregulation of eight genes, ranging from a 1.52-fold increase in ND4L to a 2.18-fold rise in Cytb transcript levels, although COI and ND5 were reduced to 0.56 and 0.67, respectively, compared with the controls. Overall, these results suggest that at 4 °C, F. kawamurai appears to have entered a hypometabolic state of hibernation, whereas its mitochondrial oxidative phosphorylation was affected at both 2 °C and 40 °C. The majority of mitochondrial PCGs exhibited substantial changes at all three temperatures, indicating that frogs such as F. kawamurai that inhabit tropical or subtropical regions are susceptible to ambient temperature changes and can quickly employ compensating adjustments to proteins involved in the mitochondrial electron transport chain. [ABSTRACT FROM AUTHOR]

The owlflies (Family Ascalaphidae) belong to the Neuroptera but are often mistaken as dragonflies because of morphological characters. To date, only three mitochondrial genomes of Ascalaphidae, namely Libelloides macaronius; Ascaloptynx appendiculatus; Ascalohybris subjacens, are published in GenBank, meaning that they are greatly under-represented in comparison with the 430 described species reported in this family. In this study, we sequenced and described the complete mitochondrial genome of Suhpalacsa longialata (Neuroptera, Ascalaphidae). The total length of the S. longialata mitogenome was 15,911 bp, which is the longest known to date among the available family members of Ascalaphidae. However, the size of each gene was similar to the other three Ascalaphidae species. The S. longialata mitogenome included a transposition of tRNACys and tRNATrp genes and formed an unusual gene arrangement tRNACys-tRNATrp-tRNATyr (CWY). It is likely that the transposition occurred by a duplication of both genes followed by random loss of partial duplicated genes. The nucleotide composition of the S. longialata mitogenome was as follows: A = 41.0%, T = 33.8%, C = 15.5%, G = 9.7%. Both Bayesian inference and ML analyses strongly supported S. longialata as a sister clade to (Ascalohybris subjacens + L. macaronius), and indicated that Ascalaphidae is not monophyletic.

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 Isonychia kiangsinensis is a circular molecule of 15,456 bp in length, containing 2 rRNA genes, 13 protein-coding genes, 22 tRNA genes, and a control region. The AT content of the overall base composition is 62.9%. The length of the control region for I. kiangsinensis is 745 bp with 68.6% AT content. In BI and ML phylogenetic trees, Isonychia kiangsinensis was a sister clade to I. ignota and Isonychiidae was shown to be the basal clade of Ephemeroptera excluding Siphluriscidae. The monophyly of the families Isonychiidae, Heptageniidae, Viemamellidae, and Baetidae and the genus Isonychia were well supported.

The complete mitochondrial genome sequence of the hybrid of Hoplobatrachus chinensis (♀) × H. rugulosus (♂) was obtained in this study. The circular mitochondrial genome was 20,282 bp in length (including extra ND5 genes). Compared with the complete mitogenome of the parents, the results indicated that the mitochondria of the hybrid tiger frog was consistent with a maternal inheritance. Phylogenetic analyses using concatenated nucleotide sequences of the 11 protein-coding genes with two different methods (maximum likelihood and MrBayes analysis) both highly supported a close relationship of the hybrid frogs with the Chinese tiger frog (=H. chinensis).

The mitochondrial genome of Epeorus herklotsi (Ephemeroptera: Heptageniidae) is a circular molecule of 15,801 bp in length with a base composition of 32.7% A, 32.9% T, 21.5% C, 13.0% G, including extra tRNAMet gene. The IMQM tRNA cluster is found in E. herklotsi as well as Parafornuru youi and two species of Epeorus (KM244708, KJ493406), while the typical IQM tRNA cluster is found in Paegniodes cupulatus. In BI and ML phylogenetic trees, the monophyly of the families Heptageniidae, Baetidae, and Ephemerellidae are highly supported. E. herklotsi is a sister clade to Epeorus sp2. (KJ493406).

In this study, the mitochondrial (mt) genomes of Siphluriscus chinensis (Ephemeroptera: Siphluriscidae) were evaluated in specimens collected from two sites in China: Niutou Mountain, Zhejiang Province (S. chinensis NTS) and Leigong Mountain, Guizhou Province (S. chinensis LGS) and were successfully sequenced. The lengths of the mt genomes of S. chinensis NTS and S. chinensis LGS were 15,904 bp (ON729390) and 15,212 bp (ON729391), respectively. However, an in-depth comparison of the two mt genomes showed significant differences between the specimens collected from the two sites. A detailed analysis of the genetic distance between S. chinensis NTS and S. chinensis LGS was undertaken to further achieve an accurate delimitation of S. chinensis. The genetic distance between S. chinensis NTS and the other three species within Siphluriscidae was a high value, above 12.2%. The two mt genomes were used to reconstruct phylogenetic relationships and estimate divergence time. The results demonstrated robust differences between S. chinensis NTS and S. chinensis LGS, which revealed that a kind of cryptic species existed. Maximum likelihood (ML) and Bayesian inference (BI) analyses produced well-supported phylogenetic trees that showed evolutionary relationships between Siphluriscidae (((S. chinensis HQ875717 + S. chinensis MF352165) + S. chinensis LGS) + S. chinensis NTS). The most recent common ancestor (MRCA) of four species within Siphluriscidae began to diversify during the Neogene [11.80 million years ago (Mya); 95% highest posterior densities (HPD) = 6.17–19.28 Mya], and S. chinensis NTS was first to diverge from the branches of S. chinensis LGS. In short, based on mitochondrial genomes, our results showed that the specimens collected from Leigong Mountain, Guizhou Province (S. chinensis LGS) belonged to S. chinensis, and the specimens collected from Niutou Mountain, Zhejiang Province (S. chinensis NTS) were a cryptic species of S. chinensis. [ABSTRACT FROM AUTHOR]

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. 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. [ABSTRACT FROM AUTHOR]

Simple Summary: The family Vietnamellidae (Ephemeroptera) is one of the oldest insect families in the world. However, there are still controversies about the phylogenetic relationships among Vietnamellidae, Ephemerellidae, and Teloganodidae. The mitochondrial (mt) genome can be used to discuss phylogenetic relationships and cryptic species. We sequenced and annotated three complete mt genomes of Vietnamella sinensis from three different populations, identifying a cryptic species of V. sinensis and discuss the phylogenetic relationships of Vietnamellidae. Based on the genetic distance of the whole mt genomes, the phylogenetic relationship of three populations was uncovered and the divergence time of V. sinensis QY indicated that it was a cryptic species of V. sinensis. Ephemeroptera (Insecta: Pterygota) are widely distributed all over the world with more than 3500 species. During the last decade, the phylogenetic relationships within Ephemeroptera have been a hot topic of research, especially regarding the phylogenetic relationships among Vietnamellidae. In this study, three mitochondrial genomes from three populations of Vienamella sinensis collected from Tonglu (V. sinensis TL), Chun'an (V. sinensis CN), and Qingyuan (V. sinensis QY) in Zhejiang Province, China were compared to discuss the potential existence of cryptic species. We also established their phylogenetic relationship by combining the mt genomes of 69 Ephemeroptera downloaded from NCBI. The mt genomes of V. sinensis TL, V. sinensis CN, and V. sinensis QY showed the same gene arrangement with lengths of 15,674 bp, 15,674 bp, and 15,610 bp, respectively. Comprehensive analyses of these three mt genomes revealed significant differences in mt genome organization, genetic distance, and divergence time. Our results showed that the specimens collected from Chun'an and Tonglu in Zhejiang Province, China belonged to V. sinensis, and the specimens collected from Qingyuan, Zhejiang Province, China were a cryptic species of V. sinensis. In maximum likelihood (ML) and Bayesian inference (BI) phylogenetic trees, the monophyly of the family Vietnamellidae was supported and Vietnamellidae has a close relationship with Ephemerellidae. [ABSTRACT FROM AUTHOR]

The phylogenetic relationship of Ephemeridae (Insect: Ephemeroptera) remains hotly debated using mitochondrial (mt) genomes. All previously reported mt genomes of Ephemeridae belong to the genus Ephemera. This study provides the first complete mt genome sequence from the genus Hexagenia with an analysis of the mitogenome of Hexagenia rigida Mc Dunnough, 1924 (Ephemeroptera: Ephemeridae) and providing new information to discuss the phylogenetic relationships within Ephemeroptera. The complete mt genome of H. rigida was a circular molecule of 16,159 bp in length, containing 37 genes (2 rRNA genes, 13 protein-coding genes, 22 tRNA genes), which showed the typical mt gene arrangement of insects. The AT content of the whole genome was 70.0% and the length of the control region was 1091 bp. All protein-coding genes used ATN as the start codon, and most PCGs used TAA/TAG as the stop codons excluding COI, COII, ND5 and Cyt b that used T as the stop codon. BI and ML phylogenetic trees constructed from 27 species of 13 families showed that Ephemeridae is a sister clade to the clade Polymitarcyidae. [ABSTRACT FROM AUTHOR]

Zoodes fulguratus is a common species of Cerambycidae, reported from Vietnam, Nepal, Laos, Burma, and China. To date, no mitochondrial genomes of the genus Zoodes have been reported. In this study, we sequenced and analyzed the first mitochondrial genome of Z. fulguratus to discuss its phylogenetic relationship within the subfamily Cerambycinae. This mitochondrial genome showed the typical insect gene arrangement: a circular molecule of 15,885 bp long with 13 protein-coding genes, two ribosomal RNA genes (rRNAs), and 22 transfer RNA genes (tRNAs). The AT content of the whole mitogenome was 74.2% with a high asymmetric nucleotide presentation revealed by a positive AT-skew (0.067) and a negative GC-skew (–0.178), whereas the AT content of the A + T rich region was 80%. The Maximum likelihood (ML) and Bayesian inference (BI) phylogenetic analyses showed that Z. fulguratus is a sister clade of Gnatholea eburifera. [ABSTRACT FROM AUTHOR]

The relationship of Matrona and Atrocalopteryx (Odonata: Calopterygidae) is still unclear. To better understand the phylogenetic relationship of Matrona and Atrocalopteryx, we sequenced and annotated two complete mitochondrial genomes of Matrona basilaris sampled from two different locations. The length of the two complete mitochondrial genomes of M. basilaris is 16,149 bp and 15,893 bp for the specimens collected in Jinhua, Zhejiang Province and Tianmushan, Zhejiang Province, China, respectively. The two mitochondrial genomes include the typical invertebrate set of 37 genes: 13 protein-coding genes (PGCs), 22 tRNA genes, and 2 rRNA genes. The nucleotide composition of the mitogenome is similar to other odonates with high content of A + T (68.9%) and all PCGs use ATN as the start codon. Tandem repeats were detected in the control regions of the two M. basilaris samples that accounted for the different sequence lengths of the mitochondrial genomes from the two locations. Finally, BI and ML phylogenetic analysis based on the concatenated nucleotide sequences of the 13 PCGs supported the conclusion that M. basilaris is a sister clade to Atrocalopteryx melli. [ABSTRACT FROM AUTHOR]