The first complete mitochondrial genome of the genus Laelaps with novel gene arrangement reveals extensive rearrangement and phylogenetics in the superfamily Dermanyssoidea.

We collected 56 specimens of Laelaps chini from the endemic Hengduan Mountain rat species (Eothenomys miletus) and obtained the first complete mitochondrial genome of L. chini by next-generation sequencing (NGS). The L. chini mitogenome is 16,507 bp in size and contains 37 genes and a control region of 2380 bp in length. The L. chini mitogenome has a high AT content and a compact arrangement with four overlapping regions ranging from 1 to 2 bp and 16 spacer regions ranging from 1 to 48 bp. We analyzed 13 protein-coding genes of L. chini mitogenome  and found that protein-coding genes in the L. chini mitogenome preferred codons ending in A/U and codon usage pattern was mainly influenced by natural selection. Cox1 has the slowest evolution rate and cox3 has the fastest evolution rate. We combined the mitochondrial genome of eight species of gamasid mites in the superfamily Dermanyssoidea from Genbank and the L. chini mitochondrial genome to analyze its rearrangement patterns and breakpoint numbers. We found that the L. chini mitogenome showed a novel arrangement pattern and nine species of gamasid mites in the superfamily Dermanyssoidea, which have been sequenced complete mitochondrial genomes to date, all showed different degrees of rearrangement. Laelaps chini, Echinolaelaps echidninus and Echinolaelaps fukinenensis were closely related species based on genetic distance and phylogenetic analyses. Notably they are clustered with Varroa destructor of the family Varroidae, suggesting that the family Varroidae is more closely related to the family Laelapidae, but more data are needed to test whether Varroa can be classified under the family Laelapidae. The L. chini mitogenome is the first complete mitochondrial genome for the genus Laelaps, and contributes to further exploration of the mitochondrial gene rearrangements and phylogeny for the superfamily Dermanyssoidea.

Highly active repeat-mediated recombination in the mitogenome of the aquatic grass Hygroryza aristata.

BACKGROUND: Floating bamboo (Hygroryza aristata) is an endangered species with a narrow native distribution and is renowned for its unique aesthetic qualities, which holds significant ecological and ornamental value. However, the lack of genetic information research, with only one complete plastome available, significantly hampers conservation efforts and further research for this species. RESULTS: In this research, we sequenced and assembled the organelle genomes of floating bamboo, including the mitogenome (587,847 bp) and plastome (135,675 bp). The mitogenome can recombine into various configurations, which are mediated by 25 repeat pairs (13 SRs, 6 MRs, 1 LR, and 5 CRs). LR1 and SR5 are particularly notable as they have the ability to combine with other contigs, forming complex repeat units that facilitate further homologous recombination. The rate of homologous recombination varies significantly among species, yet there is still a pronounced positive correlation observed between the length of these repeat pairs and the rate of recombination they mediate. The mitogenome integrates seven intact protein-coding genes from the chloroplast. The codon usage patterns in both organelles are similar, with a noticeable bias towards C and T on the third codon. The gene map of Poales shows the entire loss of rpl6, succinate dehydrogenase subunits (sdh3 and sdh4). Additionally, the BOP clade retained more variable genes compared to the PACMAD clade. CONCLUSIONS: We provided a high-quality and well-annotated mitogenome for floating bamboo and demonstrated the presence of diverse configurations. Our study has revealed the correlation between repeat length and their corresponding recombination rate despite variations among species. Although the mitogenome can potentially exist in the form of a unicircular in vivo, this occurrence is rare and may not be stable.

Insights into phylogenetic positions and distribution patterns: Complete mitogenomes of two sympatric Asian horned toads in Boulenophrys (Anura: Megophryidae).

Boulenophrys sangzhiensis and Boulenophrys tuberogranulata, two narrow-distributed toad species within the Megophryidae family in southern China, are experiencing population declines due to habitat loss and degradation. Despite their critical conservation status, the two species remain largely overlooked in public and scientific spheres. This study presented the first sequencing, assembly, and annotation of the complete mitogenomes of both species using next-generation sequencing. The mitogenome of B. sangzhiensis was 16,950 bp, while that of B. tuberogranulata was 16,841 bp, each comprising 13 protein-coding genes (PCGs), 22 transfer RNA genes (tRNAs), two ribosomal RNA genes (rRNAs), and a noncoding control region (D-loop). The gene content, nucleotide composition, and evolutionary rates of each mitogenome were analyzed. Both mitogenomes exhibited negative AT skew and GC skew with high A + T content. ATP8 exhibited the highest evolutionary rate, while COI had the lowest. A phylogenetic analysis based on 28 mitogenomes revealed two major clades of Megophryidae, supporting the classification of two subfamilies, Megophryinae and Leptobrachiinae. Within the subfamily Megophryinae, the genus Boulenophrys was divided into two species groups. Intriguingly, despite coexisting in Zhangjiajie City, B. sangzhiensis and B. tuberogranulata exhibited distinct origins from the two different species groups, underscoring the unique role of the coexisting area Zhangjiajie in driving their speciation and preserving their current populations. A parallel pattern was also identified in the Leptobrachiinae genus Leptobrachium within the same region. This study provided valuable data references and enhanced our understanding of the molecular characteristics of these threatened amphibian species.

New geographical records and molecular investigation of the ciliate Spirostomum.

Spirostomum is a genus of large ciliates, and its species are distributed worldwide. However, there has been limited research conducted on their geographical distribution and genomics. We obtained nine samples of ciliates from eight regions in Liaoning Province, China, and conducted a study on their geographical distribution and characteristics. Morphological and second-generation high-throughput sequencing methods were applied to identify the species, and a phylogenetic tree was established to gain a deeper understanding of the geographical distribution and evolutionary relationships of Spirostomum in Northeast China. The results identified Spirostomum yagiui and Spirostomum subtilis as a newly recorded species in Northeast China region. There are now five species of Spirostomum that have been recorded in China, and new details on the genomic characteristics of Spirostomum yagiui were provided. In addition, this study also identified the main branches of Spirostomum teres and Spirostomum minus in northern China, and provided a theoretical basis for the existence of hidden species. Spirostomum yagiui is the first species in the family Spirostomidae to have undergone mitochondrial genome sequencing.

Graph-based mitochondrial genomes of three foundation species in the Saccharum genus.

KEY MESSAGE: We reported the graph-based mitochondrial genomes of three foundation species (Saccharum spontaneum, S. robustum and S. officinarum) for the first time. The results revealed pan-structural variation and evolutionary processes in the mitochondrial genomes within Saccharum. Saccharum belongs to the Andropogoneae, and cultivars species in Saccharum contribute nearly 80% of sugar production in the world. To explore the genomic studies in Saccharum, we assembled 15 complete mitochondrial genomes (mitogenome) of three foundation species (Saccharum spontaneum, S. robustum and S. officinarum) using Illumina and Oxford Nanopore Technologies sequencing data. The mitogenomes of the three species were divided into a total of eight types based on contig numbers and linkages. All mitogenomes in the three species encoded 51 unique genes, including 32 protein-coding, 3 ribosomal RNA (rRNA) and 16 transfer RNA (tRNA) genes. The existence of long and short-repeat-mediated recombinations in the mitogenome of S. officinarum and S. robustum was revealed and confirmed through PCR validation. Furthermore, employing comparative genomics and phylogenetic analyses of the organelle genomes, we unveiled the evolutionary relationships and history of the major interspecific lineages in Saccharum genus. Phylogenetic analyses of homologous fragments between S. officinarum and S. robustum showed that S. officinarum and S. robustum are phylogenetically distinct and that they were likely parallel rather than domesticated. The variations between ancient (S. sinense and S. barberi) and modern cultivated species (S. hybrid) possibly resulted from hybridization involving different S. officinarum accessions. Lastly, this project reported the first graph-based mitogenomes of three Saccharum species, and a systematic comparison of the structural organization, evolutionary processes, and pan-structural variation of the Saccharum mitogenomes revealed the differential features of the Saccharum mitogenomes.

Assemble and comparative analysis of the mitochondrial genome of Rhododendron delavayi: Insights into phylogenetic relationships and genomic variations.

Rhododendron delavayi, a notable ornamental plant primarily found in regions of China like Yunnan and Guizhou provinces, holds substantial horticultural value. To elucidate the systematic phylogenetic relationships and organelle genomic differences within R. delavayi and related Rhododendron species, we conducted sequencing and assembly of the complete mitochondrial genome of R. delavayi. The full-length mitochondrial genome of it was a singular circular molecule spanning 1,009,263 bp, comprising 53 protein-coding genes, including 18 transfer RNA (tRNA) genes, 3 ribosomal RNA (rRNA) genes, and 32 protein-coding genes. A total of 1,182 simple sequence repeats (SSRs) loci were identified in the R. delavayi mitochondrial genome, primarily consisting of single nucleotide, dinucleotide, and trinucleotide repeats. Nucleotide diversity analysis highlighted five genes (atp6, atp9, cox2, nad1, and rpl10) with the highest diversity within the mitochondrial genomes of Rhododendron genus. Comparative analysis of the mitochondrial genome of R. delavayi with those of four other Rhododendron species indicated complex rearrangements in 21 genes, including rps4, nad6, rps3, atp6, cob, atp9, nad7, among others. The mitochondrial phylogenetic tree revealed a close relationship between R. delavayi and R. decorum, forming a sister clade to R. × pulchrum and R. simsii. Furthermore, 126 plastid-to-mitochondrial gene transfers in R. delavayi were identified, ranging from 30 bp to 19,385 bp. These fragments collectively constituted 47.54 % and 9.52 % of the chloroplast and mitochondrial genomes (202,169 bp), respectively. Complex mitochondrial-to-mitochondrial transfers were also observed, with 843 identified fragments totaling 312,036 bp (30.92 % of the mitochondrial genome). Segments exceeding 10 kb may mediate homologous recombination within the mitochondrial molecules. Remarkably, our study underscores that the mitochondrial genome of R. delavayi was the largest reported within the Rhododendron genus to date. The intricate rearrangements observed in the mitochondrial genomes of Rhododendron species, alone with the identification of five potential molecular marker sites, provided valuable insights for species classification and parentage identification within the Rhododendron genus.

Case report: Echinococcus multilocularis infection in a dog showing gastrointestinal signs in Hokkaido, Japan.

Echinococcus multilocularis is a cestode that causes human alveolar echinococcosis, a lethal zoonotic disease distributed in the northern hemisphere. The life cycle of this parasite is maintained in nature by voles as intermediate hosts and foxes as definitive hosts in Hokkaido, Japan. Although dogs are also susceptible to the parasite, the infection has been considered typically asymptomatic. We report the detection of E. multilocularis eggs in the diarrheal feces of a dog with chronic gastrointestinal signs, which disappeared after anthelmintic treatment. The mitochondrial genome sequence constructed by sequencing of the overlapping PCRs using DNA from the eggs was identical to the most predominant haplotype previously reported in red foxes in Hokkaido. This case highlights that Echinococcus infection should be considered as a differential diagnosis for diarrheal dogs in the disease endemic areas. Further efforts are needed to accumulate parasite genotypes in domestic dogs as well as humans to assess the risk of human infection from dogs.

Mitochondrial genome of the fluke pond snail, Austropeplea cf. brazieri (Gastropoda: Lymnaeidae).

BACKGROUND: Lymnaeid snails of the genus Austropeplea are an important vector of the liver fluke (Fasciola hepatica), contributing to livestock production losses in Australia and New Zealand. However, the species status within Austropeplea is ambiguous due to heavy reliance on morphological analysis and a relative lack of genetic data. This study aimed to characterise the mitochondrial genome of A. cf. brazieri, an intermediate host of liver fluke in eastern Victoria. METHODS: The mitochondrial genome was assembled and annotated from a combination of second- and third-generation sequencing data. For comparative purposes, we performed phylogenetic analyses of the concatenated nucleotide sequences of the mitochondrial protein-coding genes, cytochrome c oxidase subunit 1 and 16S genes. RESULTS: The assembled mt genome was 13,757 base pairs and comprised 37 genes, including 13 protein-coding genes, 22 transfer RNA genes and 2 ribosomal RNA genes. The mt genome length, gene order and nucleotide compositions were similar to related species of lymnaeids. Phylogenetic analyses of the mt nucleotide sequences placed A. cf. brazieri within the same clade as Orientogalba ollula with strong statistical supports. Phylogenies of the cox1 and 16S mt sequences were constructed due to the wide availability of these sequences representing the lymnaeid taxa. As expected in both these phylogenies, A. cf. brazieri clustered with other Austropeplea sequences, but the nodal supports were low. CONCLUSIONS: The representative mt genome of A. cf. brazieri should provide a useful resource for future molecular, epidemiology and parasitological studies of this socio-economically important lymnaeid species.

Comparison and phylogenetic analysis of the mitochondrial genomes of Synodontis eupterus and Synodontis polli.

We aimed to distinguish Synodontis eupterus and Synodontis polli. We performed sequencing and bioinformatic analysis of their mitochondrial genomes and constructed a phylogenetic tree of Mochokidae fish using maximum likelihood and Bayesian methods based on protein-coding gene (PCG) sequences of 14 Mochokidae species. The total length of the S. eupterus mitochondrial genome was 16,579 bp, including 13 (PCGs), 22 tRNA genes, two rRNA genes, and one D-loop, with an AT-biased nucleotide composition (56.0%). The total length of the S. polli mitochondrial genome was 16,544 bp, including 13 PCGs, 22 tRNA genes, two rRNA genes, and one D-loop, with an AT-biased nucleotide composition (55.0%). In both species, except for COI, PCGs use ATG as the starting codon, the vast majority use TAG or TAA as the ending codon, and a few use incomplete codons (T - or TA -) as the ending codon. Phylogenetic analysis showed that S. eupterus and Synodontis clarias converged into one branch, S. polli and Synodontis petricola converged into one branch, Mochokiella paynei, Mochokus brevis, and nine species of the genus Synodontis converged into one branch, and M. paynei clustered with the genus Synodontis. This study lays a foundation for rebuilding a clearer Mochokidae fish classification system.

A reproducible workflow for assembling the mitochondrial genome of Acheta domesticus (Orthoptera: Gryllidae).

In this study, we report the assembly and annotation of the mitochondrial genome (mitogenome) of Acheta domesticus from breeding facility, a species commonly known as the house cricket. This species is considered to be an important edible cricket. The mitogenome was assembled using a reproducible protocol implemented on the Galaxy Europe Server, which involved uploading paired-end fastq reads for bioinformatic analysis. The resulting mitogenome is 15,784 base pairs in length and has a GC content of 29.05%. The nucleotide composition of this mitogenome is similar to that of other insect mitogenomes, with A, T, C, and G nucleotides comprising 39.2%, 31.7%, 19.6%, and 9.5% of the mitogenome, respectively. The gene organization of the A. domesticus mitogenome is identical to that of other cricket species. The mitogenome consists of 37 genes, including 13 protein-coding genes, 22 tRNA genes, and two rRNA genes. The congruence between PCA and Bayesian evolutionary tree analysis in clustering the divergent A. domesticus sequences highlights these genomes as candidates for further study to elucidate their distinct features and evolutionary history.

The complete mitogenome of the Atlantic longnose chimaera Rhinochimaera atlantica (Holt & Byrne, 1909).

Holocephali is a subclass of chondrichthyans with ample geographic distribution in marine ecosystems. Holocephalan species are organized into three families: Callorhinchidae, Chimaeridae, and Rhinochimaeridae. Despite the critical ecological and evolutionary importance, genomic information from holocephalans is still scarce, particularly from rhinochimaerids. The present study provides the first complete mitogenome of the Atlantic longnose chimaera Rhinochimaera atlantica (Holt & Byrne, 1909). The whole mitogenome was sequenced from an R. atlantica specimen, collected on the Porcupine Bank (NE Atlantic), by Illumina high-throughput sequencing. The R. atlantica mitogenome has 17,852 nucleotides with 13 protein-coding genes, 22 transfer RNA, and two ribosomal RNA genes. Nine of these genes are in the complementary strand. This mitogenome has a GC content of 41.5% and an AT content of 58.5%. The phylogenetic reconstruction provided here, using all the available complete and partial Holocephali mitogenomes, places R. atlantica in the Rhinochimaeridae family, as expected. This genomic resource will be useful in the genomic characterization of this species.

Accumulation of DNA G-quadruplex in mitochondrial genome hallmarks mesenchymal senescence.

Searching for biomarkers of senescence remains necessary and challenging. Reliable and detectable biomarkers can indicate the senescence condition of individuals, the need for intervention in a population, and the effectiveness of that intervention in controlling or delaying senescence progression and senescence-associated diseases. Therefore, it is of great importance to fulfill the unmet requisites of senescence biomarkers especially when faced with the growing global senescence nowadays. Here, we established that DNA G-quadruplex (G4) in mitochondrial genome was a reliable hallmark for mesenchymal senescence. Via developing a versatile and efficient mitochondrial G4 (mtG4) probe we revealed that in multiple types of senescence, including chronologically healthy senescence, progeria, and replicative senescence, mtG4 hallmarked aged mesenchymal stem cells. Furthermore, we revealed the underlying mechanisms by which accumulated mtG4, specifically within respiratory chain complex (RCC) I and IV loci, repressed mitochondrial genome transcription, finally impairing mitochondrial respiration and causing mitochondrial dysfunction. Our findings endowed researchers with the visible senescence biomarker based on mitochondrial genome and furthermore revealed the role of mtG4 in inhibiting RCC genes transcription to induce senescence-associated mitochondrial dysfunction. These findings depicted the crucial roles of mtG4 in predicting and controlling mesenchymal senescence.

Complete Mitochondrial Genome of Clupeonella abrau (Maliatsky, 1930) (Clupeiformes), an Endemic Freshwater Fish from the Lake Abrau (Russia).

The Abrau sprat (tyulka or sardelka) Clupeonella abrau (Maliatsky, 1930) is an endemic fish of the Lake Abrau (Krasnodar Krai, Russia). The full C. abrau mitochondrial genome (16 650 bp) showed a gene arrangement conserved in Clupeidae and 98.8% similarity with the mitochondrial genome of the related species Black and Caspian Sea sprat C. cultriventris from the Black Sea. The COX1 gene sequence was additionally studied in a museum specimen collected in the Lake Abrau in 1938. Variability in modern Abrau sprat COX1 gene locus was estimated at approximately 0.15%, the difference between C. abrau and C. cultriventris was 1.2%, and the difference between the museum and modern C. abrau specimens from the Lake Abrau was 0.92%. The study confirmed that the Abrau sprat is present in the fish community and is capable of reproducing in the lake. Various scenarios were proposed to explain colonization of the Lake Abrau by C. abrau.

Far away from home? Ancient DNA shows the presence of bicolored shrew (Crocidura leucodon) in Bronze Age Denmark.

An excavation of an Early Iron Age village near Aalborg in Denmark uncovered the jaws and skull fragments from a small mammal that were morphologically identified to the genus Crocidura (white-toothed shrews). Three Crocidura species are known from prehistoric continental Europe but none of them are distributed in Scandinavia, which is why this surprising finding warranted further analyses. The bone was radiocarbon-dated to 2840-2750 calibrated years before present (cal. BP), corresponding to the Late Bronze Age and hence earlier than the Iron Age archeological context in which it was found. Using highly optimized ancient DNA protocols, we extracted DNA from one tooth and shotgun-sequenced the sample to reconstruct a near-complete mitochondrial reference genome (17,317 bp, 32.6× coverage). Phylogenetic analyses determined this specimen as a bicolored shrew (Crocidura leucodon) but with a phylogenetic position basal to the clade of known sequences from this species. The confirmation of Crocidura presence in Denmark by the Late Bronze Age sheds new light on the prehistoric natural history of Scandinavia. We discuss the implications of this finding from both zoo-archeological and ecological perspectives. Furthermore, the mitochondrial genome reconstructed in this study offers a valuable resource for future research exploring the genetic makeup and evolutionary history of Eurasian shrew populations.

Comparative mitochondrial genome analysis provides new insights into the classification of Modiolinae.

BACKGROUND: Mitochondrial genomes have become a powerful tool for studying molecular genetics and phylogeny of mollusks. Currently, the position of Modiolinae within Mytilidae and the taxonomic and phylogenetic relationships within Modiolinae were still controversial. This study focuses on the complete mitochondrial genomes of two species: Modiolus modulaides (Röding, 1798) and Modiolus auriculatus Krauss, 1848, which have not been sequenced before. METHODS AND RESULTS: We assembled and characterized the mitochondrial genomes of M. modulaides and M. auriculatus and then analyzed the phylogenetic relationships. The mitochondrial genomes of M. modulaides and M. auriculatus were 15,422 bp and 16,027 bp, respectively. Both of them were composed of 36 functional genes, including 12 protein-coding genes, 22 transfer RNAs, and 2 ribosomal RNAs. All protein-coding genes showed A + T bias, positive GC skews, and negative AT skews in nucleotide composition. Phylogenetic analysis based on the mitochondrial genomes showed that Modiolinae and Bathymodiolinae clustered together to form a sister relationship. Seven Modiolinae species were divided into two clades: L1 (M. modulaides, M. auriculatus and Modiolus philippinarum Hanley, 1843) and L2 [Modiolus modiolus (Linnaeus, 1758), Modiolus kurilensis Bernard, 1983, Modiolus nipponicus (Oyama, 1950), and Modiolus comptus (Sowerby III, 1915)]. The divergence time of the two clades was approximately 105.75 Ma. Furthermore, the transfer RNA gene rearrangement, longer genetic distance, and greater genetic differentiation were confirmed between the L1 and L2 clades, as well as differences in the external characteristics of the shells of the two clades. CONCLUSIONS: Based on the molecular data, it was speculated that species from the L1 clade might belong to other genera or new genera. This study provides molecular information for further taxonomic and phylogenetic studies of Mytilidae.

The complete mitochondrial genome of a ground beetle Synuchus nitidus (Carabidae: Harpalinae: Sphodrini) from South Korea.

The ground beetle Synuchus nitidus (Motschulsky, 1861) (Carabidae: Harpalinae: Sphodrini) is one of the most common species in the forests of South Korea, which has the potential to be utilized as an environmental indicator. Here, we characterized the complete mitochondrial genome (mitogenome) of S. nitidus, which is the first in the harpaline tribe Sphodrini. Its genome is 16,392 bp in length and composed of 13 protein-coding genes (PCGs), 22 tRNA genes, two rRNA genes, and an A + T rich region. In addition, we reconstructed a maximum likelihood tree to elucidate the phylogenetic position of Sphodrini among the seven harpaline tribes using nucleotide sequences of the 13 PCGs. The ML tree supported a monophyletic clade of the subfamily Harpalinae and showed a close relationship between Sphodrini and Lebinii with a low bootstrap value. The complete mitogenome of S. nitidus could be helpful for molecular species identification and exploring phylogenetic relationships among carabids.

The first mitochondrial genome of Meroplius fukuharai (Diptera: Sepsidae).

Meroplius fukuharai is an important sanitary and ecological resource insect. We sequenced and annotated the mitogenome of Meroplius fukuharai which is the first representative of the genus Meroplius with nearly complete mitochondrial data. This mitogenome is 14,803 bp long, which consists of 22 transfer RNA genes, 13 protein-coding genes (PCGs), and two ribosomal RNA genes. All genes have a conservational arrangement with other published species of Sepsidae. Our results also supported the monophyly of Sepsidae, and the genus Meroplius is more closely related to genus Sepsis, Microsepsis, and Archisepsis.

Molecular phylogeny reveals Varroa mites are not a separate family but a subfamily of Laelapidae.

Varroa mites, notorious for parasitizing honeybees, are generally classified as Varroidae. Their extremely modified morphologies and behaviors have led to debates regarding their phylogenetic position and classification as an independent family. In this study, two different datasets were employed to reconstruct the phylogenies of Varroa mites and related Laelapidae species: (1) 9257 bp from the whole 13 mitochondrial protein-coding genes of 24 taxa, (2) 3158 bp from 113 taxa using Sanger sequencing of four nuclear loci. Both mitochondrial and nuclear analyses consistently place Varroa mites within the Laelapidae. Here we propose to place Varroa mites in the subfamily Varroinae stat. nov., which represents a highly morphologically adapted group within the Laelapidae. Ancestral state reconstructions reveal that bee-associated lifestyles evolved independently at least three times within Laelapidae, with most phoretic traits originating from free-living ancestors. Our revised classification and evolutionary analyses will provide new insight into understanding the Varroa mites.

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.

De Novo Assembly and Comparative Analysis of Mitochondrial Genomes of Two Pueraria montana Varieties.

Pueraria montana is a species with important medicinal value and a complex genetic background. In this study, we sequenced and assembled the mitochondrial (mt) genomes of two varieties of P. montana. The mt genome lengths of P. montana var. thomsonii and P. montana var. montana were 457,390 bp and 456,731 bp, respectively. Both P. montana mitogenomes showed a multi-branched structure consisting of two circular molecules, with 56 genes annotated, comprising 33 protein-coding genes, 18 tRNA genes (trnC-GCA and trnM-CAU are multi-copy genes), and 3 rRNA genes. Then, 207 pairs of long repeats and 96 simple sequence repeats (SSRs) were detected in the mt genomes of P. montana, and 484 potential RNA-editing sites were found across the 33 mitochondrial protein-coding genes of each variety. Additionally, a syntenic sequence analysis showed a high collinearity between the two mt genomes. This work is the first to analyze the mt genomes of P. montana. It can provide information that can be used to analyze the structure of mt genomes of higher plants and provide a foundation for future comparative genomic studies and evolutionary biology research in related species.