Search

Gene rearrangements have been found in several mitochondrial genomes of Mantodea, located in the gene blocks CR-I-Q-M-ND2, COX1-K-D-ATP8 and ND3-A-R-N-S-E-F-ND5. We have sequenced one mitogenome of Amelidae (Yersinia mexicana) and six mitogenomes of Mantidae to discuss the mitochondrial gene rearrangement and the phylogenetic relationship within Mantidae. These mitogenomes showed rearrangements of tRNA genes except for Asiadodis yunnanensis and Hierodula zhangi. These novel gene rearrangements of Mantidae were primarily concentrated in the region of CR-I-Q-M-ND2, including gene translocation, duplication and pseudogenization. For the occurrences of these rearrangements, the tandem duplication-random loss (TDRL) model and slipped-strand mispairing model were suitable to explain. Large non-coding regions (LNCRs) located in the region of CR-I-Q-M-ND2 were detected in most Mantidae species, whereas some LNCRs had high similarity to the control region (CR). Both BI and ML phylogenetic analyses supported the monophyly of Mantidae and the paraphyly of Mantinae. The phylogenetic results with the gene order and the location of NCRs acted as forceful evidence that specific gene rearrangements and special LNCRs may be synapomorphies for several groups of mantises., (Copyright © 2021 Elsevier B.V. All rights reserved.)

Simple Summary: Lamiinae is the largest subfamily among the Cerambycidae (longhorn beetles), and its members are distributed worldwide. The monophyly of Lamiinae is generally recognized, but there are still diverse ideas as to whether the tribes belonging to Lamiinae are monophylic. Ambiguous classification boundaries and the existence of synonyms are major issues leading to controversies over Lamiinae classification. It is not enough to conduct research solely on the morphological characteristics and simple molecular loci of longhorn beetles. Mitochondrial genomes have proven to be reliable markers and can shed more light on phylogenetic relationships among Lamiinae. The present study resolved infra-subfamilial relationships among Lamiinae and provides more mitochondrial data for further phylogenetic research on longhorn beetles. 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. [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]

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: Skinks are the most species-rich group of lizards and are widely distributed around the world. The family Scincidae (Reptilia: Lacertiformes) includes limbed and limbless representatives occupying diverse habitats and showing a range of morphologies. Both limbed and limbless skinks have unique locomotion patterns in their habitats. Locomotion is the process of energy consumption, of which different modes may have different energy demands. As the center of energy metabolism in organisms, mitochondria provide most of the energy for physiological and biochemical activities via oxidative phosphorylation. Here, we employed mitochondrial genomes to investigate potential selective pressures among limbless skinks. Isopachys gyldenstolpei, as a typical limbless skink, has a different locomotion pattern compared to a limbed skink. Thus, I. gyldenstolpei can be used to study whether limb loss has a positive selection on mitochondrial genes. Two typical limbed skinks, Sphenomorphus indicus and Tropidophorus hainanus, were included in this study to compare the selective pressure analysis on mitochondrial genomes. In addition, the phylogenetic relationships within Scincidae are also discussed. 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. [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]

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

The first complete mitochondrial genome of Euroleon coreanus (Okamoto, 1926) was 15,797 bp in length, and contained 13 protein-coding genes, 22 transfer RNAs, two ribosomal RNAs, and the control region. Compared to the classic insect mitochondrial genome, E. coreanus showed a gene rearrangement of ND2-C-W-Y-COX1. The overall AT content of the mitochondrial genome was 75.5%. The monophyly of Ascalaphidae, Myrmeleontidae, Nemopteridae, Nymphidae, and Psychopsidae was supported in both BI and ML trees. And E. coreanus was a sister clade to the clade of genus Myrmeleon. [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]

Simple Summary: The family Scarabaeidae is one of the largest families in the insect order Coleoptera and is comprised of two quasi-systematics groups, Pleurosticti and Laparosticti. Pleurosticti is an economically important scarab group comprising about 20,000 species, the majority of which are phytophagous. Despite the innumerable studies based on ecological, molecular, and morphological characteristics, their taxonomy is still unclear and subject to many scientific hypotheses. The mitochondrial (mt) genome can provide tangible information to resolve the phylogenetic relationships within the family Scarabaeidae. However, the available mt genomes of Scarabaeidae in GenBank are underrepresented. Thus, we sequenced and analyzed 18 new phytophagous Scarabaeidae mitochondrial genomes from two subfamilies, Cetoniinae and Dynastinae, to conduct phylogenetic analyses within Scarabaeidae. This study contributes to increasing our knowledge about phytophagous 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). [ABSTRACT FROM AUTHOR]

Simple Summary: Twenty-seven complete mitochondrial genomes of Phasmatodea have been published in the NCBI. To shed light on the intra-ordinal and inter-ordinal relationships among Phasmatodea, more mitochondrial genomes of stick insects are used to explore mitogenome structures and clarify the disputes regarding the phylogenetic relationships among Phasmatodea. We sequence and annotate the first acquired complete mitochondrial genome from the family Pseudophasmatidae (Peruphasma schultei), the first reported mitochondrial genome from the genus Phryganistria of Phasmatidae (P. guangxiensis), and the complete mitochondrial genome of Orestes guangxiensis belonging to the family Heteropterygidae. We analyze the gene composition and the structure of the three mitochondrial genomes. We recover the monophyly of Phasmatodea and show the sister-group relationship between Phasmatodea and Mantophasmatodea after removal of the Embioptera and Zoraptera species. We recover the monophyly of Heteropterygidae and the paraphyly of Diapheromeridae, Phasmatidae, Lonchodidae, Lonchodinae, and Clitumninae. 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. [ABSTRACT FROM AUTHOR]

This study determined the complete mitochondrial genome of the bristletail Songmachilis xinxiangensis (Machilidae: Machilinae). It is a circular molecule of 15,473 bp long, and contains 37 genes typically found in other insects, with 13 protein-coding genes, 2 ribosomal RNA subunits genes, and 22 tRNA genes. The mitochondrial genome arrangement of S. xinxiangensis was similar to Pedetontus silvestrii (Machilidae: Petrobiinae) and Nesomachilis australica (Meinertellidae), but different from Trigoniophthalmus alternatus (Machilidae: Machilinae) in position of tRNAAla and Petrobius brevistylis (Machilidae: Petrobiinae) in the position of tRNATyr. The A+T content of total nucleotide and control region is 73.6% and 86%, respectively. [ABSTRACT FROM AUTHOR]

• The complete mitochondrial genome of C. yixingensis was sequenced for the first time. • The steady state transcript levels of mitochondrial COX3 , ND4 and ND4L genes were reduced under 40 μg/L imidacloprid for 24 h. • Steady state transcript levels of ND4 and ND4L were also significantly reduced under different concentrations. As one of the most common benthic invertebrates in freshwater, mayflies are very sensitive to changes in water quality and have high requirements for the water environment to allow their nymphs to successfully live and grow. Neonicotinoids, such as imidacloprid, can enter fresh water and pollute the aquatic environment. The present study had two goals: (1) investigate imidacloprid effects on mayfly larvae Choroterpes (Euthralus) yixingensis , and (2) contribute to the phylogenetic status of Ephemeroptera that has always been controversial. Nymphs were collected from Jinhua, China and exposed to different concentrations imidacloprid (5, 10, 20, and 40 μg/L) in the laboratory. Survival of C. yixingensis nymphs decreased as a function of time and imidacloprid concentration with only ~ 55% survival after 72 h exposure to 40 μg/L imidacloprid. After culture under 40 μg/L imidacloprid for 24 h, the steady state transcript levels of mitochondrial COX3 , ND4 and ND4L genes were reduced to just 0.07 ± 0.11, 0.30 ± 0.16, and 0.28 ± 0.13 as compared with respective control values (P < 0.01). Steady state transcript levels of ND4 and ND4L were also significantly reduced in a dose-dependent manner (P < 0.05), suggesting that the steady state transcript pattern of these genes in mayfly nymphs can change in response to different levels of environmental contamination. Hence, the mitochondrial protein-coding genes of mayflies could potentially be developed as biomarkers for water ecotoxicity monitoring in the future. In addition, we used the mitochondrial genome sequence of C. yixingensis for an assessment of the phylogenetic tree of Ephemeroptera. The monophyly of Leptophlebiidae was supported and showed that Leptophlebiidae was a sister group to the clade (Baetidae + Caenidae). [ABSTRACT FROM AUTHOR]