Comparative Mitogenome of Phylogenetic Relationships and Divergence Time Analysis within Potamanthidae (Insecta: Ephemeroptera).

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.

How Does Mitochondrial Protein-Coding Gene Expression in Fejervarya kawamurai (Anura: Dicroglossidae) Respond to Extreme Temperatures?

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.

Two-Fold ND5 Genes, Three-Fold Control Regions, lncRNA, and the "Missing" ATP8 Found in the Mitogenomes of Polypedates megacephalus (Rhacophridae: Polypedates).

In prior research on the mitochondrial genome (mitogenome) of Polypedates megacephalus, the one copy of ND5 gene was translocated to the control region (CR) and the ATP8 gene was not found. Gene loss is uncommon among vertebrates. However, in this study, we resequenced the mitogenomes of P. megacephalus from different regions using a "primer bridging" approach with Sanger sequencing technologies, which revealed the "missing" ATP8 gene in P. megacephalus as well as three other previously published Polypedates. The mitogenome of this species was found to contain two copies of the ND5 genes and three copies of the control regions. Furthermore, multiple tandem repeats were identified in the control regions. Notably, we observed that there was no correlation between genetic divergence and geographic distance. However, using the mitogenome, gene expression analysis was performed via RT-qPCR of liver samples and it was thus determined that COIII, ND2, ND4, and ND6 were reduced to 0.64 ± 0.24, 0.55 ± 0.34, 0.44 ± 0.21 and 0.65 ± 0.17, respectively, under low-temperature stress (8 °C) as compared with controls (p < 0.05). Remarkably, the transcript of long non-coding RNA (lncRNA) between positions 8029 and 8612 decreased significantly with exposure to low-temperature stress (8 °C). Antisense ND6 gene expression showed a downward trend, but this was not significant. These results reveal that modulations of protein-coding mitochondrial genes and lncRNAs of P. megacephalus play a crucial role in the molecular response to cold stress.

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 Two Mitogenomes of Hyla sanchiangensis (Anura: Hylidae), with Phylogenetic Relationships and Selection Pressure Analyses of Hylidae.

Hyla sanchiangensis (Anura: Hylidae) is endemic to China and is distributed across Anhui, Zhejiang, Fujian, Guangdong, Guangxi, Hunan, and Guizhou provinces. The mitogenomes of H. sanchiangensis from two different sites (Jinxiu, Guangxi, and Wencheng, Zhejiang) were sequenced. Phylogenetic analyses were conducted, including 38 mitogenomes of Hylidae from the NCBI database, and assessed the phylogenetic relationship of H. sanchiangensis within the analyzed dataset. Two mitogenomes of H. sanchiangensis showed the typical mitochondrial gene arrangement with 13 protein-coding genes (PCGs), two ribosomal RNA genes (12S rRNA and 16S rRNA), 22 transfer RNA (tRNA) genes, and one non-coding control region (D-loop). The lengths of the 12S rRNA and 16S rRNA genes from both samples (Jinxiu and Wencheng) were 933 bp and 1604 bp, respectively. The genetic distance (p-distance transformed into percent) on the basis of the mitogenomes (excluding the control region) of the two samples was calculated as 4.4%. Hyla sanchiangensis showed a close phylogenetic relationship with the clade of (H. annectans + H. tsinlingensis), which was supported by ML and BI analyses. In the branch-site model, five positive selection sites were found in the clade of Hyla and Dryophytes: Cytb protein (at position 316), ND3 protein (at position 85), and ND5 protein (at position 400) have one site, respectively, and two sites in ND4 protein (at positions 47 and 200). Based on the results, we hypothesized that the positive selection of Hyla and Dryophytes was due to their experience of cold stress in historical events, but more evidence is needed to support this conclusion.

The Phylogenetic Relationship of Lamiinae (Coleoptera: Cerambycidae) Using Mitochondrial Genomes.

Lamiinae is the largest subfamily of the Cerambycidae (longhorn beetles), with approximately 21,863 described species. Previous phylogenetic studies of Lamiinae showed that this subfamily was monophyletic, but the relationship between the tribes of Lamiinae is still controversial. Partial molecular data and species morphological characteristics are not sufficient to resolve species phylogenetic studies perfectly. At the same time, the full mitochondrial genome contains more comprehensive genetic data. Benefiting from the development of next-generation sequencing (NGS), mitochondrial genomes can be easily acquired and used as reliable molecular markers to investigate phylogenetic relationships within Cerambycidae. Using NGS technology, we obtained 11 mitochondrial genome sequences of Lamiinae species. Based on this newly generated mitochondrial genome dataset matrix, we reconstructed the phylogeny of Lamiinae. The Bayesian Inference and Maximum Likelihood analyses strongly support the monophyly of four tribes (Lamiini, Batocerini, Mesosini, and Saperdini), whereas the tribe Acanthocinini was identified as paraphyletic. Other mitochondrial structural features were also observed: the start codon in the nad1 gene of all 11 mitochondrial genomes is TTG; 17-22 bp intergenic spacers (IGS) with a 'TACTA' motif were found between trnS2 and nad1. Moreover, two long IGS were found in Mesosa myops and Batocera sp. Tandem repeats were found in the IGS of Batocera sp.

The Genetic Diversity and the Divergence Time in Extant Primitive Mayfly, Siphluriscus chinensis Ulmer, 1920 Using the Mitochondrial Genome.

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.

The Complete Mitochondrial Genomes of Three Sphenomorphinae Species (Squamata: Scincidae) and the Selective Pressure Analysis on Mitochondrial Genomes of Limbless Isopachys gyldenstolpei.

In order to adapt to diverse habitats, organisms often evolve corresponding adaptive mechanisms to cope with their survival needs. The species-rich family of Scincidae contains both limbed and limbless species, which differ fundamentally in their locomotor demands, such as relying on the movement of limbs or only body swing to move. Locomotion requires energy, and different types of locomotion have their own energy requirements. Mitochondria are the energy factories of living things, which provide a lot of energy for various physiological activities of organisms. Therefore, mitochondrial genomes could be tools to explore whether the limb loss of skinks are selected by adaptive evolution. Isopachys gyldenstolpei is a typical limbless skink. Here, we report the complete mitochondrial genomes of I. gyldenstolpei, Sphenomorphus indicus, and Tropidophorus hainanus. The latter two species were included as limbed comparator species to the limbless I. gyldenstolpei. The results showed that the full lengths of the mitochondrial genomes of I. gyldenstolpei, S. indicus, and T. hainanus were 17,210, 16,944, and 17,001 bp, respectively. Three mitochondrial genomes have typical circular double-stranded structures similar to other reptiles, including 13 protein-coding genes, 22 transfer RNAs, 2 ribosomal RNAs, and the control region. Three mitochondrial genomes obtained in this study were combined with fifteen mitochondrially complete genomes of Scincidae in the NCBI database; the phylogenetic relationship between limbless I. gyldenstolpei and limbed skinks (S. indicus and T. hainanus) is discussed. Through BI and ML trees, Sphenomorphinae and Mabuyinae were monophyletic, while the paraphyly of Scincinae was also recovered. The limbless skink I. gyldenstolpei is closer to the species of Tropidophorus, which has formed a sister group with (T. hainanus + T. hangman). In the mitochondrial genome adaptations between limbless I. gyldenstolpei and limbed skinks, one positively selected site was found in the branch-site model analysis, which was located in ND2 (at position 28, BEB value = 0.907). Through analyzing the protein structure and function of the selected site, we found it was distributed in mitochondrial protein complex I. Positive selection of some mitochondrial genes in limbless skinks may be related to the requirement of energy to fit in their locomotion. Further research is still needed to confirm this conclusion though.

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

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 Genetic Diversity in Thereuonema tuberculata (Wood, 1862) (Scutigeromorpha: Scutigeridae) and the Phylogenetic Relationship of Scutigeromorpha Using the Mitochondrial Genome.

Based on morphological characteristics to make species identification, the cryptic species of the Scutigeromorpha can be greatly underestimated. The mitochondrial genome provides a desirable tool for the biological identifications and the discovery of the cryptic species. The capacity to acquire mitochondrial genome sequences has substantially improved in recent years using next-generation sequencing (NGS) technology. On the basis of the next-generation sequencing, we obtained four complete mitochondrial genomes of Thereuonema tuberculata (Wood, 1862) from Nanyang, Henan Province (NY), Nanchang, Jiangxi Province (NC), Jinan, Shandong Province (JN), and Dali, Yunnan Province (DL) in China with GenBank numbers OK513221, OL449685, ON058988 and ON058989, respectively. The lengths of the four mitochondrial genomes ranged from 14,903 to 14,909 bp. The composition and order of genes of the four mitochondrial genomes were identical to the published mitochondrial genome of Scutigera coleoptrata (Linnaeus, 1758) (Scutigeromorpha: Scutigerdae). It was the first time that the tandem repeats in the control region were detected in Scutigeromorpha. We also calculated the corrected pairwise genetic distance of four complete mitochondrial genomes of T. tuberculata, ranging from 7.7 to 15.2%. The results showed that the T.tuberculata NC belonged to the typical sample of T. tuberculata, and T. tuberculata DL was hypothesized as a cryptic species of T. tuberculata. Meanwhile, T. tuberculata NY and T. tuberculata JN were hypothesized as potential cryptic species of T. tuberculata in this study. In both BI and ML trees, the monophyly of Scutigeromorpha, Scolopendromorpha, Geophilomorpha, and Lithobiomorpha was forcefully advocated. Moreover, Scutigeromorpha was recovered as the sister clade of (Scolopendromorpha + (Lithobiomorpha + Geophilomorpha)). Four specimens of T. tuberculata were clustered into one clade, which was the sister to the clade of S. coleoptrata.

The first complete mitochondrial genome of Hexagenia rigida Mc Dunnough, 1924 (Ephemeroptera: Ephemeridae) and its phylogeny.

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.

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

Cryptic Species Exist in Vietnamella sinensis Hsu, 1936 (Insecta: Ephemeroptera) from Studies of Complete Mitochondrial Genomes.

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.

Mitochondrial gene expression in different organs of Hoplobatrachus rugulosus from China and Thailand under low-temperature stress.

BACKGROUND: Hoplobatrachus rugulosus (Anura: Dicroglossidae) is distributed in China and Thailand and the former can survive substantially lower temperatures than the latter. The mitochondrial genomes of the two subspecies also differ: Chinese tiger frogs (CT frogs) display two identical ND5 genes whereas Thai tiger frogs (TT frogs) have two different ND5 genes. Metabolism of ectotherms is very sensitive to temperature change and different organs have different demands on energy metabolism at low temperatures. Therefore, we conducted studies to understand: (1) the differences in mitochondrial gene expression of tiger frogs from China (CT frogs) versus Thailand (TT frogs); (2) the differences in mitochondrial gene expression of tiger frogs (CT and TT frogs) under short term 24 h hypothermia exposure at 25 °C and 8 °C; (3) the differences in mitochondrial gene expression in three organs (brain, liver and kidney) of CT and TT frogs. RESULTS: Utilizing RT-qPCR and comparing control groups at 25 °C with low temperature groups at 8 °C, we came to the following results. (1) At the same temperature, mitochondrial gene expression was significantly different in two subspecies. The transcript levels of two identical ND5 of CT frogs were observed to decrease significantly at low temperatures (P < 0.05) whereas the two different copies of ND5 in TT frogs were not. (2) Under low temperature stress, most of the genes in the brain, liver and kidney were down-regulated (except for COI and ATP6 measured in brain and COI measured in liver of CT frogs). (3) For both CT and TT frogs, the changes in overall pattern of mitochondrial gene expression in different organs under low temperature and normal temperature was brain > liver > kidney. CONCLUSIONS: We mainly drew the following conclusions: (1) The differences in the structure and expression of the ND5 gene between CT and TT frogs could result in the different tolerances to low temperature stress. (2) At low temperatures, the transcript levels of most of mitochondrial protein-encoding genes were down-regulated, which could have a significant effect in reducing metabolic rate and supporting long term survival at low temperatures. (3) The expression pattern of mitochondrial genes in different organs was related to mitochondrial activity and mtDNA replication in different organs.

The complete mitochondrial genome of Rana johnsi (Smith, 2009) (Anura: Ranidae) and its phylogeny.

Rana johnsi (Smith 2009) firstly considered as the member of genus Pseudorana, has been moved into the genus Rana. In this study, we sequenced the complete mitochondrial (mt) genome of R. johnsi using the Sanger method. The circular mt genome was 17,873 bp in length and contains 13 protein-coding genes (PCGs), 22 transfer RNA (tRNA) genes, two ribosome RNA genes, and one control region. The overall nucleotide composition in majority-strand was 28% A, 29% T, 29% C, and 14% G. We discussed the phylogenetic relationship of R. johnsi in genus Rana using ML and BI analyses based on 13 PCGs. Excluding the clade of subgenus Lithobates, Rana draytonii was the basal clade to all other Rana species, which included R. johnsi as the basal clade. The monophyly of genus Rana was supported, whereas Pseudorana was failed to support.

The Mitochondrial Genomes of 18 New Pleurosticti (Coleoptera: Scarabaeidae) Exhibit a Novel trnQ-NCR-trnI-trnM Gene Rearrangement and Clarify Phylogenetic Relationships of Subfamilies within Scarabaeidae.

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 mitochondrial genome of Eurycantha calcarata Lucas, 1869 (Phasmatodea: Lonchodinae) and its phylogeny.

The Lonchodinae (Phasmatodea: Phasmatidae) is rich in insect species with more than 330 species of 40 genera. The phylogenetic relationships within Lonchodinae have been under debate. We successfully sequenced the complete mitogenome of Eurycantha calcarata Lucas, 1869 (Phasmatodea: Lonchodinae) with a length of 16,280 bp, which had the same genes and gene arrangements as those of various published papers on stick insects. The whole mitogenome and control region of E. calcarata had a high AT content of 78.2 and 85.9%, respectively. All PCGs used ATN as the start codon, and most PCGs used TAA/TAG as the stop codons excluding COX2 (T), COX3 (TA), and ND5 (TA). To discuss the phylogeny of Lonchodinae, we reconstructed the phylogenetic relationships of 27 species of Phasmatodea including E. calcarata and two species of Embioptera used as outgroups. In BI and ML trees, the monophyly of Lonchodinae and Necrosciinae was well supported, whereas the monophyly of Clitumninae was not recovered. These results indicated that Lonchodinae was a sister clade to Phylliinae and E. calcarata was a sister clade to Phraortes genus.

Three Complete Mitochondrial Genomes of Orestes guangxiensis, Peruphasma schultei, and Phryganistria guangxiensis (Insecta: Phasmatodea) and Their Phylogeny.

Insects of the order Phasmatodea are mainly distributed in the tropics and subtropics and are best known for their remarkable camouflage as plants. In this study, we sequenced three complete mitochondrial genomes from three different families: Orestes guangxiensis, Peruphasma schultei, and Phryganistria guangxiensis. The lengths of the three mitochondrial genomes were 15,896 bp, 16,869 bp, and 17,005 bp, respectively, and the gene composition and structure of the three stick insects were identical to those of the most recent common ancestor of insects. The phylogenetic relationships among stick insects have been chaotic for a long time. In order to discuss the intra- and inter-ordinal relationship of Phasmatodea, we used the 13 protein-coding genes (PCGs) of 85 species for maximum likelihood (ML) and Bayesian inference (BI) analyses. Results showed that the internal topological structure of Phasmatodea had a few differences in both ML and BI trees and long-branch attraction (LBA) appeared between Embioptera and Zoraptera, which led to a non-monophyletic Phasmatodea. Consequently, after removal of the Embioptera and Zoraptera species, we re-performed ML and BI analyses with the remaining 81 species, which showed identical topology except for the position of Tectarchus ovobessus (Phasmatodea). We recovered the monophyly of Phasmatodea and the sister-group relationship between Phasmatodea and Mantophasmatodea. Our analyses also recovered the monophyly of Heteropterygidae and the paraphyly of Diapheromeridae, Phasmatidae, Lonchodidae, Lonchodinae, and Clitumninae. In this study, Peruphasma schultei (Pseudophasmatidae), Phraortes sp. YW-2014 (Lonchodidae), and species of Diapheromeridae clustered into the clade of Phasmatidae. Within Heteropterygidae, O. guangxiensis was the sister clade to O. mouhotii belonging to Dataminae, and the relationship of (Heteropteryginae + (Dataminae + Obriminae)) was recovered.

Comparative Mitogenomes of Two Coreamachilis Species (Microcoryphia: Machilidae) along with Phylogenetic Analyses of Microcoryphia.

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.

Insight into the Phylogenetic Relationships among Three Subfamilies within Heptageniidae (Insecta: Ephemeroptera) along with Low-Temperature Selection Pressure Analyses Using Mitogenomes.

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.