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

The complete mitochondrial genome of Choroterpes (Euthralus) yixingensis (Ephemeroptera: Leptophlebiidae) and its mitochondrial protein-coding gene expression under imidacloprid stress.

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

The first complete mitochondrial genome of Euroleon coreanus (Okamoto, 1926) (Neuroptera: Myrmeleontidae) and its phylogeny.

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.

Novel tRNA gene rearrangements in the mitochondrial genomes of praying mantises (Mantodea: Mantidae): Translocation, duplication and pseudogenization.

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.

The first complete mitochondrial genome of Zoodes fulguratus (Gahan 1906) (Coleoptera: Cerambycidae: Cerambycinae) and its phylogeny.

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.

Mitogenome Analysis of Four Lamiinae Species (Coleoptera: Cerambycidae) and Gene Expression Responses by Monochamus alternatus When Infected with the Parasitic Nematode, Bursaphelenchus mucronatus.

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.

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

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