Mitochondrial Protein-Coding Gene Expression in the Lizard Sphenomorphus incognitus (Squamata:Scincidae) Responding to Different Temperature Stresses.

In the context of global warming, the frequency of severe weather occurrences, such as unexpected cold spells and heat waves, will grow, as well as the intensity of these natural disasters. Lizards, as a large group of reptiles, are ectothermic. Their body temperatures are predominantly regulated by their environment and temperature variations directly impact their behavior and physiological activities. Frequent cold periods and heat waves can affect their biochemistry and physiology, and often their ability to maintain their body temperature. Mitochondria, as the center of energy metabolism, are crucial for maintaining body temperature, regulating metabolic rate, and preventing cellular oxidative damage. Here, we used RT-qPCR technology to investigate the expression patterns and their differences for the 13 mitochondrial PCGs in Sphenomorphus incognitus (Squamata:Scincidae), also known as the brown forest skink, under extreme temperature stress at 4 °C, 8 °C, 34 °C, and 38 °C for 24 h, compared to the control group at 25 °C. In southern China, for lizards, 4 °C is close to lethal, and 8 °C induces hibernation, while 34/38 °C is considered hot and environmentally realistic. Results showed that at a low temperature of 4 °C for 24 h, transcript levels of ATP8, ND1, ND4, COI, and ND4L significantly decreased, to values of 0.52 ± 0.08, 0.65 ± 0.04, 0.68 ± 0.10, 0.28 ± 0.02, and 0.35 ± 0.02, respectively, compared with controls. By contrast, transcript levels of COIII exhibited a significant increase, with a mean value of 1.86 ± 0.21. However, exposure to 8 °C for 24 h did not lead to an increase in transcript levels. Indeed, transcript levels of ATP6, ATP8, ND1, ND3, and ND4 were significantly downregulated, to 0.48 ± 0.11, 0.68 ± 0.07, 0.41 ± 0.08, 0.54 ± 0.10, and 0.52 ± 0.07, respectively, as compared with controls. Exposure to a hot environment of 34 °C for 24 h led to an increase in transcript levels of COI, COII, COIII, ND3, ND5, CYTB, and ATP6, with values that were 3.3 ± 0.24, 2.0 ± 0.2, 2.70 ± 1.06, 1.57 ± 0,08, 1.47 ± 0.13, 1.39 ± 0.56, and 1.86 ± 0.12, respectively, over controls. By contrast, ND4L exhibited a significant decrease (to 0.31 ± 0.01) compared with controls. When exposed to 38 °C, the transcript levels of the 13 PCGs significantly increased, ranging from a 2.04 ± 0.23 increase in ND1 to a 6.30 ± 0.96 rise in ND6. Under two different levels of cold and heat stress, the expression patterns of mitochondrial genes in S. incognitus vary, possibly associated with different strategies employed by this species in response to low and high temperatures, allowing for rapid compensatory adjustments in mitochondrial electron transport chain proteins in response to temperature changes. Furthermore, this underscores once again the significant role of mitochondrial function in determining thermal plasticity in reptiles.

Differential Mitochondrial Genome Expression of Four Hylid Frog Species under Low-Temperature Stress and Its Relationship with Amphibian Temperature Adaptation.

Extreme weather poses huge challenges for animals that must adapt to wide variations in environmental temperature and, in many cases, it can lead to the local extirpation of populations or even the extinction of an entire species. Previous studies have found that one element of amphibian adaptation to environmental stress involves changes in mitochondrial gene expression at low temperatures. However, to date, comparative studies of gene expression in organisms living at extreme temperatures have focused mainly on nuclear genes. This study sequenced the complete mitochondrial genomes of five Asian hylid frog species: Dryophytes japonicus, D. immaculata, Hyla annectans, H. chinensis and H. zhaopingensis. It compared the phylogenetic relationships within the Hylidae family and explored the association between mitochondrial gene expression and evolutionary adaptations to cold stress. The present results showed that in D. immaculata, transcript levels of 12 out of 13 mitochondria genes were significantly reduced under cold exposure (p < 0.05); hence, we put forward the conjecture that D. immaculata adapts by entering a hibernation state at low temperature. In H. annectans, the transcripts of 10 genes (ND1, ND2, ND3, ND4, ND4L, ND5, ND6, COX1, COX2 and ATP8) were significantly reduced in response to cold exposure, and five mitochondrial genes in H. chinensis (ND1, ND2, ND3, ND4L and ATP6) also showed significantly reduced expression and transcript levels under cold conditions. By contrast, transcript levels of ND2 and ATP6 in H. zhaopingensis were significantly increased at low temperatures, possibly related to the narrow distribution of this species primarily at low latitudes. Indeed, H. zhaopingensis has little ability to adapt to low temperature (4 °C), or maybe to enter into hibernation, and it shows metabolic disorder in the cold. The present study demonstrates that the regulatory trend of mitochondrial gene expression in amphibians is correlated with their ability to adapt to variable climates in extreme environments. These results can predict which species are more likely to undergo extirpation or extinction with climate change and, thereby, provide new ideas for the study of species extinction in highly variable winter climates.

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.

Differential Mitochondrial Genome Expression of Three Sympatric Lizards in Response to Low-Temperature Stress.

Ecological factors related to climate extremes have a significant influence on the adaptability of organisms, especially for ectotherms such as reptiles that are sensitive to temperature change. Climate extremes can seriously affect the survival and internal physiology of lizards, sometimes even resulting in the loss of local populations or even complete extinction. Indeed, studies have shown that the expression levels of the nuclear genes and mitochondrial genomes of reptiles change under low-temperature stress. At present, the temperature adaptability of reptiles has rarely been studied at the mitochondrial genome level. In the present study, the mitochondrial genomes of three species of lizards, Calotes versicolor, Ateuchosaurus chinensis, and Hemidactylus bowringii, which live in regions of sympatry, were sequenced. We used RT-qPCR to explore the level of mitochondrial gene expression under low-temperature stress, as compared to a control temperature. Among the 13 protein-coding genes (PCGs), the steady-state transcript levels of ND4L, ND1, ATP6, and COII were reduced to levels of 0.61 ± 0.06, 0.50 ± 0.08, 0.44 ± 0.16, and 0.41 ± 0.09 in C. versicolor, respectively, compared with controls. The transcript levels of the ND3 and ND6 genes fell to levels of just 0.72 ± 0.05 and 0.67 ± 0.05 in H. bowringii, compared with controls. However, the transcript levels of ND3, ND5, ND6, ATP6, ATP8, Cytb, and COIII genes increased to 1.97 ± 0.15, 2.94 ± 0.43, 1.66 ± 0.07, 1.59 ± 0.17, 1.46 ± 0.04, 1.70 ± 0.16, and 1.83 ± 0.07 in A. chinensis. Therefore, the differences in mitochondrial gene expression may be internally related to the adaptative strategy of the three species under low-temperature stress, indicating that low-temperature environments have a greater impact on A. chinensis, with a small distribution area. In extreme environments, the regulatory trend of mitochondrial gene expression in reptiles is associated with their ability to adapt to extreme climates, which means differential mitochondrial genome expression can be used to explore the response of different lizards in the same region to low temperatures. Our experiment aims to provide one new research method to evaluate the potential extinction of reptile species in warm winter climates.

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.

Nine Mitochondrial Genomes of Phasmatodea with Two Novel Mitochondrial Gene Rearrangements and Phylogeny.

The classification of stick and leaf insects (Order Phasmatodea) is flawed at various taxonomic ranks due to a lack of robust phylogenetic relationships and convergent morphological characteristics. In this study, we sequenced nine new mitogenomes that ranged from 15,011 bp to 17,761 bp in length. In the mitogenome of Carausis sp., we found a translocation of trnR and trnA, which can be explained by the tandem duplication/random loss (TDRL) model. In the Stheneboea repudiosa Brunner von Wattenwyl, 1907, a novel mitochondrial structure of 12S rRNA-CR1-trnI-CR2-trnQ-trnM was found for the first time in Phasmatodea. Due to the low homology of CR1 and CR2, we hypothesized that trnI was inverted through recombination and then translocated into the middle of the control region. Control region repeats were frequently detected in the newly sequenced mitogenomes. To explore phylogenetic relationships in Phasmatodea, mtPCGs from 56 Phasmatodean species (composed of 9 stick insects from this study, 31 GenBank data, and 16 data derived from transcriptome splicing) were used for Bayesian inference (BI), and maximum likelihood (ML) analyses. Both analyses supported the monophyly of Lonchodinae and Necrosciinae, but Lonchodidae was polyphyletic. Phasmatidae was monophyletic, and Clitumninae was paraphyletic. Phyllidae was located at the base of Neophasmatodea and formed a sister group with the remaining Neophasmatodea. Bacillidae and Pseudophasmatidae were recovered as a sister group. Heteroptergidae was monophyletic, and the Heteropteryginae sister to the clade (Obriminae + Dataminae) was supported by BI analysis and ML analysis.

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