Animal mitochondrial genomes

The evolution of contemporary livestock species: Insights from mitochondrial genome

The domestication of animals marks a pivotal moment in human history, profoundly influencing our demographic and cultural progress. This process has led to significant genetic, behavioral, and physical changes in livestock species compared to their wild ancestors. Understanding the evolutionary history and genetic diversity of livestock species is crucial, and mitochondrial DNA (mtDNA) has emerged as a robust marker for investigating molecular diversity in animals. Its highly conserved gene content across animal species, minimal duplications, absence of introns, and short intergenic regions make mtDNA analysis ideal for such studies. Mitochondrial DNA analysis has uncovered distinct cattle domestication events dating back to 8000 years BC in Southwestern Asia. The sequencing of water buffalo mtDNA in 2004 provided important insights into their domestication history. Caprine mtDNA analysis identified three haplogroups, indicating varied maternal origins. Sheep, domesticated 12,000 years ago, exhibit diverse mtDNA lineages, suggesting multiple domestication events. Ovine mtDNA studies revealed clades A, B, C, and a fourth lineage, group D. The origins of domestic pigs were traced to separate European and Asian events followed by interbreeding. In camels, mtDNA elucidated the phylogeographic structure and genetic differentiation between wild and domesticated species. Horses, domesticated around 3500 BC, show significant mtDNA variability, highlighting their diverse origins. Yaks exhibit unique adaptations for high-altitude environments, with mtDNA analysis providing insights into their adaptation. Chicken mtDNA studies supported a monophyletic origin from Southeast Asia's red jungle fowl, with evidence of multiple origins. This review explores livestock evolution and diversity through mtDNA studies, focusing on cattle, water buffalo, goat, sheep, pig, camel, horse, yak and chicken. It highlights mtDNA's significance in unraveling maternal lineages, genetic diversity, and domestication histories, concluding with insights into its potential application in improving livestock production and reproduction dynamics.

Species-specific variation in mitochondrial genome tandem repeat polymorphisms in hares (Lepus spp., Lagomorpha, Leporidae) provides insight into their evolution

The non-coding regions of the mitochondrial DNAs (mtDNAs) of hares, rabbits, and pikas (Lagomorpha) contain short (∼20 bp) and long (130-160 bp) tandem repeats, absent in related mammalian orders. In the presented study, we provide in-depth analysis for mountain hare (Lepus timidus) and brown hare (L. europaeus) mtDNA non-coding regions, together with a species- and population-level analysis of tandem repeat variation. Mountain hare short tandem repeats (SRs) as well as other analyzed hare species consist of two conserved 10 bp motifs, with only brown hares exhibiting a single, more variable motif. Long tandem repeats (LRs) also differ in sequence and copy number between species. Mountain hares have four to seven LRs, median value five, while brown hares exhibit five to nine LRs, median value six. Interestingly, introgressed mountain hare mtDNA in brown hares obtained an intermediate LR length distribution, with median copy number being the same as with conspecific brown hare mtDNA. In contrast, transfer of brown hare mtDNA into cultured mtDNA-less mountain hare cells maintained the original LR number, whereas the reciprocal transfer caused copy number instability, suggesting that cellular environment rather than the nuclear genomic background plays a role in the LR maintenance. Due to their dynamic nature and separation from other known conserved sequence elements on the non-coding region of hare mitochondrial genomes, the tandem repeat elements likely to represent signatures of ancient genetic rearrangements. clarifying the nature and dynamics of these rearrangements may shed light on the possible role of NCR repeated elements in mitochondria and in species evolution.

A chromosome-level genome assembly of Cape hare (Lepus capensis)

The Cape hare (Lepus capensis) is among the most widely distributed hare species globally, inhabiting extensive regions across Africa, the Middle East, and Central Asia. However, evolutionary and genetic research on L. capensis was seriously impeded by the absence of a reference genome. Here, we assembled and constructed a chromosome-level genome of L. capensis (with scaffolds anchored to 25 chromosomes and a total assembled length of 2.9 Gb, achieving a contig N50 length of 124.44 Mb) using PacBio HiFi sequencing and Hi-C assembly technology. Evaluation using BUSCO indicated the genome assembly to be 98.2% complete. The de novo prediction revealed that repetitive sequences constitute 46.13% of the entire genome, and long interspersed nuclear elements (LINEs) constituted the largest portion. We annotated a total of 13, 868 protein-coding genes using transcriptomes from two tissues (muscle and skin). This high-quality reference genome serves as a valuable genomic resource for advancing genetic studies in this species.

The Complete Mitochondrial Genome of the Siberian Scoter Melanitta stejnegeri and Its Phylogenetic Relationship in Anseriformes

The Siberian Scoter (Melanitta stejnegeri) is a medium sea duck distinct from M. deglandi due to the absence of hybridization and differences in morphological characteristics. However, knowledge of its phylogenetic relationships within Anseriformes is limited due to a lack of molecular data. In this study, the complete mitogenome of M. stejnegeri was firstly sequenced, then annotated and used to reconstruct the phylogenetic relationships of 76 Anseriformes species. The complete mitogenome of M. stejnegeri is 16,631 bp and encodes 37 typical genes: 13 protein-coding genes, 2 ribosomal RNAs, 22 transfer RNAs, and 1 non-coding control region. Its mitogenome organization is similar to that of other Anseriformes species. The phylogenetic relationships within the genus Melanitta are initially clarified, with M. americana at the base. M. stejnegeri and M. deglandi are sister groups, clustering with M. fusca and M. perspicillata in order. Phylogenetic analysis suggests that Mareca falcata and M. strepera are sister groups, differing from previous studies. Results firstly indicate that Clangula hyemalis and Somateria mollissima are sister groups, suggesting a potentially skewed phylogenetic relationship may have been overlooked in earlier analyses relying solely on mitochondrial genomes. Our results provide new mitogenome data to support further phylogenetic and taxonomic studies of Anseriformes.

The first report of complete mitogenomes of two endangered species of genus Propomacrus (Coleoptera: Scarabaeidae: Euchirinae) and phylogenetic implications

To understand the mitochondrial genome structure of two endangered and long-armed scarab beetles, Propomacrus davidi and Propomacrus bimucronatus, their complete mitogenomes were sequenced for the first time in this study. The complete mitogenomes of P. davidi and P. bimucronatus were 18, 042 bp and 18, 104 bp in length, respectively. The gene orders of their mitogenomes were highly consistent with other Coleopteran species, and the typical ATN was used as the start codon in most protein coding genes. The incomplete stop codon T was used in cox1, cox2, and nad5, and TAN was used as a complete stop codon in most protein coding genes. All predicted tRNAs could form a typical cloverleaf secondary structure, except that trnS1 lacked the dihydrouridine arm. Based on the maximum likelihood and the Bayesian inference methods, phylogenetic trees of 50 species were reconstructed. The results showed that P. davidi, P. bimucronatus, Cheirotonus jansoni and Cheirotonus gestroi clustered in the same branch, and were the most closely related. The results supported that subfamily Euchirinae is a monophyletic group of Scarabaeidae, which was consistent with the morphological classification. These molecular data enriched the complete mitogenome database of Euchirinae, and improved our understanding of the phylogenetic relationship and evolutionary characteristics of these two endangered species.

Full-Length Transcriptome Profiling of the Complete Mitochondrial Genome of Sericothrips houjii (Thysanoptera: Thripidae: Sericothripinae) Featuring Extensive Gene Rearrangement and Duplicated Control Regions

The mitochondrial genome (mitogenome) of Thysanoptera has extensive gene rearrangement, and some species have repeatable control regions. To investigate the characteristics of the gene expression, transcription and post-transcriptional processes in such extensively gene-rearranged mitogenomes, we sequenced the mitogenome and mitochondrial transcriptome of Sericothrips houjii to analyze. The mitogenome was 14,965 bp in length and included two CRs contains 140 bp repeats between COIII-trnN (CR1) and trnT-trnP (CR2). Unlike the putative ancestral arrangement of insects, S. houjii exhibited only six conserved gene blocks encompassing 14 genes (trnL2-COII, trnD-trnK, ND2-trnW, ATP8-ATP6, ND5-trnH-ND4-ND4L and trnV-lrRNA). A quantitative transcription map showed the gene with the highest relative expression in the mitogenome was ND4-ND4L. Based on analyses of polycistronic transcripts, non-coding RNAs (ncRNAs) and antisense transcripts, we proposed a transcriptional model of this mitogenome. Both CRs contained the transcription initiation sites (TISs) and transcription termination sites (TTSs) of both strands, and an additional TIS for the majority strand (J-strand) was found within antisense lrRNA. The post-transcriptional cleavage processes followed the "tRNA punctuation" model. After the cleavage of transfer RNAs (tRNAs), COI and ND3 matured as bicistronic mRNA COI/ND3 due to the translocation of intervening tRNAs, and the 3' untranslated region (UTR) remained in the mRNAs for COII, COIII, CYTB and ND5. Additionally, isoform RNAs of ND2, srRNA and lrRNA were identified. In summary, the relative mitochondrial gene expression levels, transcriptional model and post-transcriptional cleavage process of S. houjii are notably different from those insects with typical mitochondrial gene arrangements. In addition, the phylogenetic tree of Thripidae including S. houjii was reconstructed. Our study provides insights into the phylogenetic status of Sericothripinae and the transcriptional and post-transcriptional regulation processes of extensively gene-rearranged insect mitogenomes.

The Complete Mitochondrial Genome and Phylogenetic Analysis of Rhagastis binoculata (Matsumura, 1909) (Lepidoptera: Sphingidae)

The mitochondrial genome (mitogenome) Rhagastis binoculata (Matsumura, 1909), an endemic moth species in Taiwan, was sequenced and analyzed. The complete circular mitogenome of R. binoculata is 15,303 bp and contains 13 protein-coding genes, 22 transfer RNA genes, 2 ribosomal RNA genes, and an AT-rich control region. The mitogenome has an overall nucleotide composition of 41.2% A, 11.9% C, 7.5% G, and 39.4% T, with an AT content of 80.6%. Of the protein-coding genes (PCGs), 12 start with ATG, ATT, and ATC, and COX1 starts with a "CGA" codon. All of the stop codons are "TAA, TAG, or T". Our phylogenetic analysis of 21 species of Sphingidae insects suggests that R. binoculata is clustered with Rhagastis mongoliana, which belongs to the subfamily Macroglossinae.

The redescription and complete mitogenomes of two Oxycarenus species (Hemiptera, Oxycarenidae) and phylogenetic implications

In this study, the two Oxycarenidae species, O.gossypii Horváth, 1926 and Oxycarenusbicolorheraldus Distant, 1904, are redescribed, and their complete mitogenomes are sequenced and analyzed. The phylogeny of Lygaeoidea is examined using 45 complete mitogenomes of lygaeoid species and four outgroup species. The gene orientation and arrangement of the two mitogenomes are found to be consistent with typical Lygaeoidea mitochondrial features, comprising 37 genes, including 13 PCGs, 22 tRNAs, 2 rRNAs, and a control region. Nucleotide composition of the species was biased towards A and T, with the gene order identical to the putative ancestral arrangement of insects. Start codons, stop codons, RNAs, relative synonymous codon usage (RSCU), and nucleotide diversity (Pi) of Oxycarenidae exhibited characteristics similar to other families in Lygaeoidea. Bayesian-inference (BI) and maximum-likelihood (ML) methods were employed to investigate phylogenetic relationships using PCG datasets from selected species. Phylogenetic analyses reveal slightly different topologies between BI and ML methods, with variation primarily concentrated in Colobathristidae and Rhyparochromidae. Our study confirms that the two sequenced Oxycarenidae species formed a single clade, and the position of Oxycarenidae remains stable in both ML and BI phylogenetic trees. These findings expand the mitochondrial genome databases of Lygaeoidea and provide valuable insights into the phylogenetic relationships within Lygaeoidea or Pentatomomorpha.

A chromosome-level genome assembly of the disco clam, Ctenoides ales

The bivalve subclass Pteriomorphia, which includes the economically important scallops, oysters, mussels, and ark clams, exhibits extreme ecological, morphological, and behavioral diversity. Among this diversity are five morphologically distinct eye types, making Pteriomorphia an excellent setting to explore the molecular basis for the evolution of novel traits. Of pteriomorphian bivalves, Limida is the only order lacking genomic resources, greatly limiting the potential phylogenomic analyses related to eyes and phototransduction. Here, we present a limid genome assembly, the disco clam, Ctenoides ales (C. ales), which is characterized by invaginated eyes, exceptionally long tentacles, and a flashing light display. This genome assembly was constructed with PacBio long reads and Dovetail Omni-CTM proximity-ligation sequencing. The final assembly is ∼2.3Gb and over 99% of the total length is contained in 18 pseudomolecule scaffolds. We annotated 41,064 protein coding genes and reported a BUSCO completeness of 91.9% for metazoa_obd10. Additionally, we report a complete and annotated mitochondrial genome, which also had been lacking from Limida. The ∼20Kb mitogenome has 12 protein coding genes, 22 tRNAs, 2 rRNA genes, and a 1,589 bp duplicated sequence containing the origin of replication. The C. ales nuclear genome size is substantially larger than other pteriomorphian genomes, mainly accounted for by transposable element sequences. We inventoried the genome for opsins, the signaling proteins that initiate phototransduction, and found that, unlike its closest eyed-relatives, the scallops, C. ales lacks duplication of the rhabdomeric Gq-protein-coupled opsin that is typically used for invertebrate vision. In fact, C. ales has uncharacteristically few opsins relative to the other pteriomorphian families, all of which have unique expansions of xenopsins, a recently discovered opsin subfamily. This chromosome-level assembly, along with the mitogenome, is a valuable resource for comparative genomics and phylogenetics in bivalves and particularly for the understudied but charismatic limids.

Evolutionary Insights from the Mitochondrial Genome of Oikopleura dioica: Sequencing Challenges, RNA Editing, Gene Transfers to the Nucleus, and tRNA Loss

Sequencing the mitochondrial genome of the tunicate Oikopleura dioica is a challenging task due to the presence of long poly-A/T homopolymer stretches, which impair sequencing and assembly. Here, we report on the sequencing and annotation of the majority of the mitochondrial genome of O. dioica by means of combining several DNA and amplicon reads obtained by Illumina and MinIon Oxford Nanopore Technologies with public RNA sequences. We document extensive RNA editing, since all homopolymer stretches present in the mitochondrial DNA correspond to 6U-regions in the mitochondrial RNA. Out of the 13 canonical protein-coding genes, we were able to detect eight, plus an unassigned open reading frame that lacked sequence similarity to canonical mitochondrial protein-coding genes. We show that the nad3 gene has been transferred to the nucleus and acquired a mitochondria-targeting signal. In addition to two very short rRNAs, we could only identify a single tRNA (tRNA-Met), suggesting multiple losses of tRNA genes, supported by a corresponding loss of mitochondrial aminoacyl-tRNA synthetases in the nuclear genome. Based on the eight canonical protein-coding genes identified, we reconstructed maximum likelihood and Bayesian phylogenetic trees and inferred an extreme evolutionary rate of this mitochondrial genome. The phylogenetic position of appendicularians among tunicates, however, could not be accurately determined.

SmithHunter: a workflow for the identification of candidate smithRNAs and their targets

BACKGROUND

SmithRNAs (Small MITochondrial Highly-transcribed RNAs) are a novel class of small RNA molecules that are encoded in the mitochondrial genome and regulate the expression of nuclear transcripts. Initial evidence for their existence came from the Manila clam Ruditapes philippinarum, where they have been described and whose activity has been biologically validated through RNA injection experiments. Current evidence on the existence of these RNAs in other species is based only on small RNA sequencing. As a preliminary step to characterize smithRNAs across different metazoan lineages, a dedicated, unified, analytical workflow is needed.

RESULTS

We propose a novel workflow specifically designed for smithRNAs. Sequence data (from small RNA sequencing) uniquely mapping to the mitochondrial genome are clustered into putative smithRNAs and prefiltered based on their abundance, presence in replicate libraries and 5' and 3' transcription boundary conservation. The surviving sequences are subsequently compared to the untranslated regions of nuclear transcripts based on seed pairing, overall match and thermodynamic stability to identify possible targets. Ample collateral information and graphics are produced to help characterize these molecules in the species of choice and guide the operator through the analysis. The workflow was tested on the original Manila clam data. Under basic settings, the results of the original study are largely replicated. The effect of additional parameter customization (clustering threshold, stringency, minimum number of replicates, seed matching) was further evaluated.

CONCLUSIONS

The study of smithRNAs is still in its infancy and no dedicated analytical workflow is currently available. At its core, the SmithHunter workflow builds over the bioinformatic procedure originally applied to identify candidate smithRNAs in the Manila clam. In fact, this is currently the only evidence for smithRNAs that has been biologically validated and, therefore, the elective starting point for characterizing smithRNAs in other species. The original analysis was readapted using current software implementations and some minor issues were solved. Moreover, the workflow was improved by allowing the customization of different analytical parameters, mostly focusing on stringency and the possibility of accounting for a minimal level of genetic differentiation among samples.

Characterization of the complete mitochondrial genome of a newly discovered torrent catfish, Liobagrus geumgangensis, and their phylogenetic relationships

BACKGROUND

A new Liobagrus fish was reported from the Korean Peninsula, but research on this taxon is lacking. Moreover, existing research on the mitogenome of the genus Liobagrus in Korea is very limited, and no studies have been conducted on structural characteristics of transfer RNA (tRNA) or gene order comparisons between taxa; instead, research has been restricted to basic phylogeny.

OBJECTIVE

The complete mitochondrial genome of Liobagrus geumgangensis was analyzed for the first time. We then aimed to reconstruct the phylogenetic relationships of the genus Liobagrus and estimate the divergence time of speciation events.

METHODS

We used a dissected fin clip from an adult of Liobagrus geumgangensis. Genomic DNA was extracted and analyzed with whole genome sequencing (WGS) and assembled by the NOVOPlasty method. The mitogenome sequence was annotated, and a genome map, tRNA structure, and phylogenetic tree were constructed using maximum likelihood analysis. In addition, divergence time was estimated.

RESULTS

The mitochondrial genome was 16,522 bp in length and comprised 37 genes. The overall base composition was 30.5% A, 25.5% T, 28.4% C, and 15.7% G. Most tRNAs exhibited the typical clover leaf shape, except trnS1. Phylogenetic analysis revealed that Liobagrus geumgangensis clustered within a clade with four other Liobagrus species exclusive to the southern region of the Korean Peninsula. Its divergence was estimated to have occurred during the late Miocene.

CONCLUSION

Characteristics of Liobagrus geumgangensis mitogenome were consistent with those of other torrent catfish species. Time scale estimation revealed distinct groupings, with some distributed across mainland Asia and others in the southern region of the Korean Peninsula. Notably, the Korean Peninsula group was identified as its own lineage, comprising entirely endemic species.

Morphological observation, molecular identification and evolutionary analysis of Hydatigera kamiyai found in Neodon fuscus from the Qinghai-Tibetan plateau

Hydatigera kamiyai (H. kamiyai) is a new species within Hydatigera that has recently been resurrected. Voles and cats are hosts of H. kamiyai and have a certain impact on its health and economy. Moreover, the Qinghai-Tibetan plateau (QTP) is a research hotspot representing Earth's biodiversity, as its unique geographical environment and climatic conditions support the growth of a variety of mammals and provide favorable conditions for various parasites to complete their life history. The aim of this study was to reveal the phylogenetic relationships and divergence times of H. kamiyai strains isolated from Neodon fuscus on the QTP using morphological and molecular methods. In this study, we morphologically observed H. kamiyai and sequenced the whole mitochondrial genome. Then, we constructed phylogenetic trees with the maximum likelihood (ML) and Bayesian inference (BI) methods. The GTR alternative model was selected for divergence time analysis. These data demonstrated that the results were consistent with the general morphological characteristics of Hydatigera. The whole genome of H. kamiyai was 13,822 bp in size, and the A + T content (73%) was greater than the G + C content (27%). The Ka/Ks values were all <1, indicating that all 13 protein-coding genes (13 PCGs) underwent purifying selection during the process of evolution. The phylogenetic tree generated based on the 13 PCGs, cytochrom oxidase subunit I (COI), 18S rRNA and 28S rRNA revealed close phylogenetic relationships between H. kamiyai and Hydatigera, with high node support for the relationship. The divergence time based on 13 PCGs indicated that H. kamiyai diverged approximately 11.3 million years ago (Mya) in the Miocene. Interestingly, it diverged later than the period of rapid uplift in the QTP. We also speculated that H. kamiyai differentiation was caused by host differentiation due to the favorable living conditions brought about by the uplift of the QTP. As there have been relatively few investigations on the mitochondrial genome of H. kamiyai, our study could provide factual support for further studies of H. kamiyai on the QTP. We also emphasized the importance of further studies of its hosts, Neodon fuscus and cats, which will be important for further understanding the life cycle of H. kamiyai.

Genomic insights into adaptation to bipedal saltation and desert-like habitats of jerboas

Jerboas is a lineage of small rodents displaying atypical mouse-like morphology with elongated strong hindlimbs and short forelimbs. They have evolved obligate bipedal saltation and acute senses, and been well-adapted to vast desert-like habitats. Using a newly sequenced chromosome-scale genome of the Mongolian five-toed jerboa (Orientallactaga sibirica), our comparative genomic analyses and in vitro functional assays showed that the genetic innovations in both protein-coding and non-coding regions played an important role in jerboa morphological and physiological adaptation. Jerboa-specific amino acid substitutions, and segment insertions/deletions (indels) in conserved non-coding elements (CNEs) were found in components of proteoglycan biosynthesis pathway (XYLT1 and CHSY1), which plays an important role in limb development. Meanwhile, we found specific evolutionary changes functionally associated with energy or water metabolism (e.g., specific amino acid substitutions in ND5 and indels in CNEs physically near ROR2) and senses (e.g., expansion of vomeronasal receptors and the FAM136A gene family) in jerboas. Further dual-luciferase reporter assay verified that some of the CNEs with jerboa-specific segment indels exerted a significantly different influence on luciferase activity, suggesting changes in their regulatory function in jerboas. Our results revealed the potential molecular mechanisms underlying jerboa adaptation since the divergence from the Eocene-Oligocene transition, and provided more resources and new insights to enhance our understanding of the molecular basis underlying the phenotypic diversity and the environmental adaptation of mammals.

Complete mitochondrial genome and phylogenetic analysis of Mancinella alouina

BACKGROUND

The Muricidae family in the Class Gastropoda comprises numerous species with a vast range of morphological features and a worldwide presence. The phylogeny of the Muricidae has been analyzed in previous studies; however, the evolutionary relationships among the main branches of the Muricidae remain unknown.

METHODS AND RESULTS

In the present study, the mitochondrial genome of Mancinella alouina was sequenced. The mitochondrial genome was found to be 16,671 bp in length and made up of 37 genes (13 protein-coding genes, 22 transfer RNA and 2 ribosomal RNA genes). The genome has an A-T-rich region (66.5% A + T content) and all of the PCGs use the ATN start codon and the TAG or TAA stop codons. The mitochondrial gene arrangement of Mancinella alouina is similar to that of other Muricidae, except for Ocinebrellus inornatus and Ceratostoma burnetti. On the basis of a flexible molecular clock model, time-calibrated phylogenetic results indicate that the genus Mancinella diverged roughly 18.09 Mya, and that the family Muricidae emerged in the Late Cretaceous.

CONCLUSIONS

This study reveals the structural and sequence information features of the mitochondrial genome of Mancinella alouina. This study provides evidence for the relationships within the family Muricidae at the molecular level, and infer the divergence time. The results of phylogenetic analyses strongly support the current classification.

Characteristics and Comparative Analysis of Six Mitogenomes of Genus Kiefferulus Goetghebuer, 1922 (Diptera: Chironomidae)

Chironomidae is a cosmopolitan and species-rich family of insects, with many species serving as useful indicators of aquatic ecosystem health. In this study, we newly sequenced six species of Kiefferulus Goetghebuer, 1922 (Chironomidae: Chironominae) by high-throughput sequencing technology. We analyzed characters of the mitochondrial genome, including the sequence length, nucleotide composition, and evolutionary rates of this genus. The size of the newly obtained sequences ranged from 15,588 to 15,767 bp, and all of them included 22 tRNAs, 13 PCGs, 2 rRNAs, and 1 CR. The CR showed the highest AT content relative to the PCGs, rRNAs, and tRNAs. Relative synonymous codon usage analysis showed that UUA, UUU, and AUU are the preferred codons. The ratio of nonsynonymous (Ka) to synonymous (Ks) substitution rates showed that all Ka/Ks of PCGs were lower than 1, with ATP8 having the highest evolution rate, while COX1 exhibited the lowest evolution rate. We reconstructed the phylogenetic relationship of the genus Kiefferulus based on eight species (six ingroups and two outgroups), using five matrices and employing Maximum likelihood and Bayesian inference approaches. Phylogenetic analysis of the Kiefferulus showed that six species within this genus were classified into a monophyletic clade.

Phylogenetic and Comparative Analysis of Cryptochironomus, Demicryptochironomus and Harnischia Inferred from Mitogenomes (Diptera: Chironomidae)

(1) Background: Mitochondrial genomes have been extensively employed as a crucial marker in numerous dipteran families for understanding phylogenetics and systematics relations, thereby playing a pivotal role in molecular biology studies. The phylogenetic relationship of the Harnischia generic complex remains contentious due to the paucity of taxonomic and molecular data. Specifically, the evolutionary relationships among Cryptochironomus, Demicryptochironomus, and Harnischia are still unclear. (2) Methods: In this study, Polypedilum and Endochironomus were used as outgroups to analyze phylogenetic relationships among Cryptochironomus, Demicryptochironomus, and Harnischia, mitogenomes of four Cryptochironomus, two Demicryptochironomus, two Harnischia, and two Cladopelma were newly sequenced. Subsequently, we conducted a thorough analysis of the nucleotide composition, sequence length, and evolutionary rate. (3) Results: All mitogenomes exhibited structural conservation, with all genes consistently arranged in the identical order as that of the ancestral mitogenome. Nucleotide composition varied significantly among different genes, and the control region displayed the highest A + T content. All protein-coding genes undergo rigorous purification selection, with the ATP8 gene exhibiting the most rapid evolutionary rate among them. Utilizing Bayesian Inference (BI) and Maximum Likelihood (ML) methods across various databases, we reconstructed the phylogenetic relationships among the genera within the Harnischia generic complex, drawing insights from an analysis of 14 mitochondrial genomes. (4) Conclusions: Our results showed that the monophyly of the genera Harnischia was well supported in all topologies; Cryptochironomus is sister to Demicryptochironomus.

Complete mitochondrial genome of Penicillidia dufourii (Diptera: Hippoboscoidea: Nycteribiidae) and phylogenetic relationship

Penicillidia dufourii (Westwood 1834) is a specialized parasite categorized under family Nycteribiidae that prefers to parasitize the body surface of various bats under the genus Myotis. Many species of the family Nycteribiidae are carriers of various pathogens; however, research on P. dufourii remains scarce, and studies on its molecular identification and population genetic structure are still lacking. In this study, the complete mitochondrial genome of P. dufourii was elucidated for the first time using Illumina sequencing. The mitochondrial genome is 15,354 bp in size and encodes approximately 37 genes, including 13 protein-coding genes, 22 tRNA genes, 2 rRNA genes, and 1 control region. Analysis of 13 protein-coding genes revealed that UUA, UCA, CGA, and GGA were the most common codons, while nad4L had the fastest evolutionary rate and cox1 the slowest. Phylogenetic analysis based on the mitochondrial genome indicated that P. dufourii is clustered with other species of the family Nycteribiidae and is most closely related to Nycteribia parvula and Phthiridium szechuanum.

Characterization and Phylogenetic Analysis of the Complete Mitochondrial Genome of Triplophysa microphthalma

The complete mitochondrial genome has been extensively utilized in studies related to phylogenetics, offering valuable perspectives on evolutionary relationships. The mitochondrial genome of the fine-eyed plateau loach, Triplophysa microphthalma, has not attracted much attention, although this species is endemic to China. In this study, we characterized the mitochondrial genome of T. microphthalma and reassessed the classification status of its genus. The complete mitochondrial genome of T. microphthalma was 16,591 bp and contained thirty-seven genes, including thirteen protein-coding genes (PCGs), two ribosomal RNA genes (rRNAs), and twenty-two transfer RNA genes (tRNAs). All but one of the thirteen PCGs had the regular start codon ATG; the gene cox1 started with GTG. Six PCGs had incomplete stop codons (T--). These thirteen PCGs are thought to have evolved under purifying selection, and the mitogenome shared a high degree of similarity with the genomes of species within the genus Leptobotia. All tRNA genes exhibited the standard clover-shaped structure, with the exception of the trnS1 gene, which lacked a DHU stem. A phylogenetic analysis indicated that T. microphthalma was more closely related to species within the genus Triplophysa than to those in Barbatula. The present study contributes valuable genomic information for T. microphthalma, and offers new perspectives on the phylogenetic relationships among species of Triplophysa and Barbatula. The findings also provide essential data that can inform the management and conservation strategies for T. microphthalma and other species of Triplophysa and Barbatula.

The Mitochondrial Genome of Ylistrum japonicum (Bivalvia, Pectinidae) and Its Phylogenetic Analysis

The Ylistrum japonicum is a commercially valuable scallop known for its long-distance swimming abilities. Despite its economic importance, genetic and genomic research on this species is limited. This study presents the first complete mitochondrial genome of Y. japonicum. The mitochondrial genome is 19,475 bp long and encompasses 13 protein-coding genes, three ribosomal RNA genes, and 23 transfer RNA genes. Two distinct phylogenetic analyses were used to explore the phylogenetic position of the Y. japonicum within the family Pectinidae. Based on one mitochondrial phylogenetic analysis by selecting 15 Pectinidae species and additional outgroup taxa and one single gene phylogenetic analysis by 16S rRNA, two phylogenetic trees were constructed to provide clearer insights into the evolutionary placement of Y. japonicum within the family Pectinidae. Our analysis reveals that Ylistrum is a basal lineage to the Pectininae clade, distinct from its previously assigned tribe, Amusiini. This study offers critical insights into the genetic makeup and evolutionary history of Y. japonicum, enhancing our knowledge of this economically vital species.

The complete mitochondrial genome analysis of Locastra muscosalis (Walker, 1866) (Lepidoptera: Pyralidae) and phylogenetic implications

The complete mitochondrial genome of Locastra muscosalis (Walker, 1866) was sequenced and characterized in this study, which was the first reported complete mitogenome of the genus Locastra. The mitogenome of L. muscosalis has a total length of 15,177 bp, encompassing 13 protein-coding genes (PCGs), 22 transfer RNA genes (tRNAs), 2 ribosomal RNA genes (rRNAs), and an A-T rich region. Phylogenetic analysis revealed that L. muscosalis is closely associated with Orthaga euadrusalis. These data will serve as a valuable foundation for future investigations into the Epipaschiinae and Pyralidae evolutionary history.

Accumulation of Large Lineage-Specific Repeats Coincides with Sequence Acceleration and Structural Rearrangement in Plantago Plastomes

Repeats can mediate rearrangements and recombination in plant mitochondrial genomes and plastid genomes. While repeat accumulations are linked to heightened evolutionary rates and complex structures in specific lineages, debates persist regarding the extent of their influence on sequence and structural evolution. In this study, 75 Plantago plastomes were analyzed to investigate the relationships between repeats, nucleotide substitution rates, and structural variations. Extensive repeat accumulations were associated with significant rearrangements and inversions in the large inverted repeats (IRs), suggesting that repeats contribute to rearrangement hotspots. Repeats caused infrequent recombination that potentially led to substoichiometric shifting, supported by long-read sequencing. Repeats were implicated in elevating evolutionary rates by facilitating localized hypermutation, likely through DNA damage and repair processes. This study also observed a decrease in nucleotide substitution rates for loci translocating into IRs, supporting the role of biased gene conversion in maintaining lower substitution rates. Combined with known parallel changes in mitogenomes, it is proposed that potential dysfunction in nuclear-encoded genes associated with DNA replication, recombination, and repair may drive the evolution of Plantago organellar genomes. These findings contribute to understanding how repeats impact organellar evolution and stability, particularly in rapidly evolving plant lineages.

The mitochondrial genome of Bottapotamon fukienense (Brachiura: Potamidae) is fragmented into two chromosomes

BACKGROUND

China is the hotspot of global freshwater crab diversity, but their wild populations are facing severe pressures associated with anthropogenic factors, necessitating the need to map their taxonomic and genetic diversity and design conservation policies.

RESULTS

Herein, we sequenced the mitochondrial genome of a Chinese freshwater crab species Bottapotamon fukienense, and found that it is fragmented into two chromosomes. We confirmed that fragmentation was not limited to a single specimen or population. Chromosome 1 comprised 15,111 base pairs (bp) and there were 26 genes and one pseudogene (pseudo-nad1) encoded on it. Chromosome 2 comprised 8,173 bp and there were 12 genes and two pseudogenes (pseudo-trnL2 and pseudo-rrnL) encoded on it. Combined, they comprise the largest mitogenome (23,284 bp) among the Potamidae. Bottapotamon was the only genus in the Potamidae dataset exhibiting rearrangements of protein-coding genes. Bottapotamon fukienense exhibited average rates of sequence evolution in the dataset and did not differ in selection pressures from the remaining Potamidae.

CONCLUSIONS

This is the first experimentally confirmed fragmentation of a mitogenome in crustaceans. While the mitogenome of B. fukienense exhibited multiple signs of elevated mitogenomic architecture evolution rates, including the exceptionally large size, duplicated genes, pseudogenisation, rearrangements of protein-coding genes, and fragmentation, there is no evidence that this is matched by elevated sequence evolutionary rates or changes in selection pressures.

The Phylogenetic Relationships of Major Lizard Families Using Mitochondrial Genomes and Selection Pressure Analyses in Anguimorpha

Anguimorpha, within the order Squamata, represents a group with distinct morphological and behavioral characteristics in different ecological niches among lizards. Within Anguimorpha, there is a group characterized by limb loss, occupying lower ecological niches, concentrated within the subfamily Anguinae. Lizards with limbs and those without exhibit distinct locomotor abilities when adapting to their habitats, which in turn necessitate varying degrees of energy expenditure. Mitochondria, known as the metabolic powerhouses of cells, play a crucial role in providing approximately 95% of an organism's energy. Functionally, mitogenomes (mitochondrial genomes) can serve as a valuable tool for investigating potential adaptive evolutionary selection behind limb loss in reptiles. Due to the variation of mitogenome structures among each species, as well as its simple genetic structure, maternal inheritance, and high evolutionary rate, the mitogenome is increasingly utilized to reconstruct phylogenetic relationships of squamate animals. In this study, we sequenced the mitogenomes of two species within Anguimorpha as well as the mitogenomes of two species in Gekkota and four species in Scincoidea. We compared these data with the mitogenome content and evolutionary history of related species. Within Anguimorpha, between the mitogenomes of limbless and limbed lizards, a branch-site model analysis supported the presence of 10 positively selected sites: Cytb protein (at sites 183 and 187), ND2 protein (at sites 90, 155, and 198), ND3 protein (at site 21), ND5 protein (at sites 12 and 267), and ND6 protein (at sites 72 and 119). These findings suggested that positive selection of mitogenome in limbless lizards may be associated with the energy requirements for their locomotion. Additionally, we acquired data from 205 mitogenomes from the NCBI database. Bayesian inference (BI) and Maximum Likelihood (ML) trees were constructed using the 13 mitochondrial protein-coding genes (PCGs) and two rRNAs (12S rRNA and 16S rRNA) from 213 mitogenomes. Our phylogenetic tree and the divergence time estimates for Squamata based on mitogenome data are consistent with results from previous studies. Gekkota was placed at the root of Squamata in both BI and ML trees. However, within the Toxicofera clade, due to long-branch attraction, Anguimorpha and (Pleurodonta + (Serpentes + Acrodonta)) were closely related groupings, which might indicate errors and also demonstrate that mitogenome-based phylogenetic trees may not effectively resolve long-branch attraction issues. Additionally, we reviewed the origin and diversification of Squamata throughout the Mesozoic era, suggesting that Squamata originated in the Late Triassic (206.05 Mya), with the diversification of various superfamilies occurring during the Cretaceous period. Future improvements in constructing squamate phylogenetic relationships using mitogenomes will rely on identifying snake and acrodont species with slower evolutionary rates, ensuring comprehensive taxonomic coverage of squamate diversity, and increasing the number of genes analyzed.

Daphnia sp. (Branchiopoda: Cladocera) Mitochondrial Genome Gene Rearrangement and Phylogenetic Position Within Branchiopoda

In high-altitude (4500 m) freshwater lakes, Daphnia is the apex species and the dominant zooplankton. It frequently dwells in the same lake as the Gammarid. Branchiopoda, a class of Arthropoda, Crustacea, is a relatively primitive group in the subphylum Crustacea, which originated in the Cambrian period of the Paleozoic. The complete mitogenome sequence of Daphnia sp. (Branchiopoda: Cladocera) was sequenced and annotated in this study and deposited in GenBank. The sequence structure of this species was studied by comparing the original sequences with BLAST. In addition, we have also researched the mechanisms of their mitochondrial gene rearrangement by establishing a model. We have used the Bayesian inference [BI] and maximum likelihood [ML] methods to proceed with phylogenetic analysis inference, which generates identical phylogenetic topology that reveals the phylogenetic state of Daphnia. The complete mitogenome of Daphnia sp. shows that it was 15,254 bp in length and included two control regions (CRs) and 37 genes (13 protein-coding genes, 22 tRNAs and two ribosomal RNAs [16S and 12S]). In addition to tRNA-Ser (GCT), other tRNAs have a typical cloverleaf secondary structure. Meanwhile, the mitogenome of Daphnia sp. was clearly rearranged when compared to the mitogenome of typical Daphnia. In a word, we report a newly sequenced mitogenome of Daphnia sp. with a unique rearrangement phenomenon. These results will be helpful for further phylogenetic research and provide a foundation for future studies on the characteristics of the mitochondrial gene arrangement process in Daphnia.

Single-mitochondrion sequencing uncovers distinct mutational patterns and heteroplasmy landscape in mouse astrocytes and neurons

BACKGROUND

Mitochondrial (mt) heteroplasmy can cause adverse biological consequences when deleterious mtDNA mutations accumulate disrupting "normal" mt-driven processes and cellular functions. To investigate the heteroplasmy of such mtDNA changes, we developed a moderate throughput mt isolation procedure to quantify the mt single-nucleotide variant (SNV) landscape in individual mouse neurons and astrocytes. In this study, we amplified mt-genomes from 1645 single mitochondria isolated from mouse single astrocytes and neurons to (1) determine the distribution and proportion of mt-SNVs as well as mutation pattern in specific target regions across the mt-genome, (2) assess differences in mtDNA SNVs between neurons and astrocytes, and (3) study co-segregation of variants in the mouse mtDNA.

RESULTS

(1) The data show that specific sites of the mt-genome are permissive to SNV presentation while others appear to be under stringent purifying selection. Nested hierarchical analysis at the levels of mitochondrion, cell, and mouse reveals distinct patterns of inter- and intra-cellular variation for mt-SNVs at different sites. (2) Further, differences in the SNV incidence were observed between mouse neurons and astrocytes for two mt-SNV 9027:G > A and 9419:C > T showing variation in the mutational propensity between these cell types. Purifying selection was observed in neurons as shown by the Ka/Ks statistic, suggesting that neurons are under stronger evolutionary constraint as compared to astrocytes. (3) Intriguingly, these data show strong linkage between the SNV sites at nucleotide positions 9027 and 9461.

CONCLUSIONS

This study suggests that segregation as well as clonal expansion of mt-SNVs is specific to individual genomic loci, which is important foundational data in understanding of heteroplasmy and disease thresholds for mutation of pathogenic variants.

The Mitogenomic Landscape of Hexacorallia Corals: Insight into Their Slow Evolution

The utility of the mitochondrial genomes (mitogenomes) in analyzing the evolutionary history of animals has been proven. Five deep-sea corals (Bathypathes sp.1, Bathypathes sp.2, Schizopathidae 1, Trissopathes sp., and Leiopathes sp.) were collected in the South China Sea (SCS). Initially, the structures and collinearity of the five deep-sea coral mitogenomes were analyzed. The gene arrangements in the five deep-sea coral mitogenomes were similar to those in the order Antipatharia, which evidenced their conservation throughout evolutionary history. Additionally, to elucidate the slow evolutionary rates in Hexacorallia mitogenomes, we conducted comprehensive analyses, including examining phylogenetic relationships, performing average nucleotide identity (ANI) analysis, and assessing GC-skew dissimilarity combining five deep-sea coral mitogenomes and 522 reference Hexacorallia mitogenomes. Phylogenetic analysis using 13 conserved proteins revealed that species clustered together at the order level, and they exhibited interspersed distributions at the family level. The ANI results revealed that species had significant similarities (identity > 85%) within the same order, while species from different orders showed notable differences (identity < 80%). The investigation of the Hexacorallia mitogenomes also highlighted that the GC-skew dissimilarity was highly significant at the order level, but not as pronounced at the family level. These results might be attributed to the slow evolution rate of Hexacorallia mitogenomes and provide evidence of mitogenomic diversity. Furthermore, divergence time analysis revealed older divergence times assessed via mitogenomes compared with nuclear data, shedding light on significant evolutionary events shaping distinct orders within Hexacorallia corals. Those findings provide new insights into understanding the slow evolutionary rates of deep-sea corals in all lineages of Hexacorallia using their mitogenomes.

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

The first two complete mitochondrial genomes for the genus Anagyrus (Hymenoptera, Encyrtidae) and their phylogenetic implications

Anagyrus, a genus of Encyrtidae (Hymenoptera, Chalcidoidea), represents a successful group of parasitoid insects that attack various mealybug pests of agricultural and forestry plants. Until now, only 20 complete mitochondrial genomes have been sequenced, including those in this study. To enrich the diversity of mitochondrial genomes in Encyrtidae and to gain insights into their phylogenetic relationships, the mitochondrial genomes of two species of Anagyrus were sequenced, and the mitochondrial genomes of these species were compared and analyzed. Encyrtid mitochondrial genomes exhibit similarities in nucleotide composition, gene organization, and control region patterns. Comparative analysis of protein-coding genes revealed varying molecular evolutionary rates among different genes, with six genes (ATP8, ND2, ND4L, ND6, ND4 and ND5) showing higher rates than others. A phylogenetic analysis based on mitochondrial genome sequences supports the monophyly of Encyrtidae; however, the two subfamilies, Encyrtinae and Tetracneminae, are non-monophyletic. This study provides valuable insights into the phylogenetic relationships within the Encyrtidae and underscores the utility of mitochondrial genomes in the systematics of this family.

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.

JW, Vieira DBR, Brandão RCF, Medeiros DBDA, Nunes Neto JP. Sequencing and Description of the Mitochondrial Genome of Orthopodomyia fascipes (Diptera: Culicidae)

The genus Orthopodomyia Theobald, 1904 (Diptera: Culicidae) comprises 36 wild mosquito species, with distribution largely restricted to tropical and temperate areas, most of which are not recognized as vectors of epidemiological importance due to the lack of information related to their bionomy and involvement in the cycle transmission of infectious agents. Furthermore, their evolutionary relationships are not completely understood, reflecting the scarcity of genetic information about the genus. Therefore, in this study, we report the first complete description of the mitochondrial genome of a Neotropical species representing the genus, Orthopodomyia fascipes Coquillet, 1906, collected in the Brazilian Amazon region. Using High Throughput Sequencing, we obtained a mitochondrial sequence of 15,598 bp, with an average coverage of 418.5×, comprising 37 functional subunits and a final portion rich in A + T, corresponding to the control region. The phylogenetic analysis, using Maximum Likelihood and Bayesian Inference based on the 13 protein-coding genes, corroborated the monophyly of Culicidae and its two subfamilies, supporting the proximity between the tribes Orthopodomyiini and Mansoniini, partially disagreeing with previous studies based on the use of molecular and morphological markers. The information generated in this study contributes to a better understanding of the taxonomy and evolutionary history of the genus and other groups of Culicidae.

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.

Unraveling the complex evolutionary features of the Cinnamomum camphora mitochondrial genome

KEY MESSAGE

We reported the mitochondrial genome of Cinnamomum camphora for the first time, revealing frequent rearrangement events in the non-coding regions of Magnoliids mitochondrial genomes. As one of the representative species in the Lauraceae family of Magnoliids, Cinnamomum camphora holds significant economic and ecological value. In this study, the mitochondrial genome (mitogenome) of C. camphora was complete assembled and annotated using PacBio HiFi sequencing. The C. camphora mitogenome is characterized by a branch structure, spans 900,894 bp, and contains 43 protein-coding genes (PCGs), 24 tRNAs, and 3 rRNAs. Most of these PCGs are under purifying selection, with only two (ccmFc and rps7) exhibiting signs of positive selection. The C. camphora mitogenome contains numerous repetitive sequences and intracellular gene transfers, with a total of 36 mitochondrial plastid DNAs, amounting to a combined length of 23,816 bp. Comparative analysis revealed that the non-coding regions of Magnoliids mitogenomes have undergone frequent rearrangements during evolution, but the coding sequences remain highly conserved (more than 98% similarity for protein-coding sequences). Furthermore, a maximum-likelihood phylogenetic tree was reconstructed based on 25 PCGs from 23 plant mitogenomes. The analysis supports the closest relationship between C. camphora and C. chekiangense, consistent with the APG IV classification system. This study elucidates the unique evolutionary features of the C. camphora mitogenome, which will provide valuable insights into the study of genetics and evolution of the family Lauraceae.

The First Three Mitochondrial Genomes for the Characterization of the Genus Egeirotrioza (Hemiptera: Triozidae) and Phylogenetic Implications

(1) Background: Mitochondrial genomes are important markers for the study of phylogenetics and systematics. Triozidae includes some primary pests of Populus euphratica. The phylogenetic relationships of this group remain controversial due to the lack of molecular data. (2) Methods: Mitochondria of Egeirotrioza Boselli were sequenced and assembled. We analyzed the sequence length, nucleotide composition, and evolutionary rate of Triozidae, combined with the 13 published mitochondrial genomes. (3) Results: The evolutionary rate of protein-coding genes was as follows: ATP8 > ND6 > ND5 > ND2 > ND4 > ND4L > ND1 > ND3 > APT6 > CYTB > COX3 > COX2 > COX1. We reconstructed the phylogenetic relationships of Triozidae based on 16 triozid mitochondrial genomes (thirteen ingroups and three outgroups) using the maximum likelihood (ML) and Bayesian inference (BI) approaches. The phylogenetic analysis of the 16 Triozidae mitochondrial genomes showed that Egeirotrioza was closely related to Leptynoptera. (4) Conclusions: We have identified 13 PCGs, 22 tRNAs, 2 rRNAs, and 1 control region (CR) of all newly sequenced mitochondrial genomes, which were the mitochondrial gene type in animals. The results of this study provide valuable genomic information for the study of psyllid species.

The initiation of mitochondrial DNA replication

Mitochondrial DNA replication is initiated by the transcription of mitochondrial RNA polymerase (mtRNAP), as mitochondria lack a dedicated primase. However, the mechanism determining the switch between continuous transcription and premature termination to generate RNA primers for mitochondrial DNA (mtDNA) replication remains unclear. The pentatricopeptide repeat domain of mtRNAP exhibits exoribonuclease activity, which is required for the initiation of mtDNA replication in Drosophila. In this review, we explain how this exonuclease activity contributes to primer synthesis in strand-coupled mtDNA replication, and discuss how its regulation might co-ordinate mtDNA replication and transcription in both Drosophila and mammals.

Single-Mitochondrion Sequencing Uncovers Distinct Mutational Patterns and Heteroplasmy Landscape in Mouse Astrocytes and Neurons

Background

Mitochondrial (mt) heteroplasmy can cause adverse biological consequences when deleterious mtDNA mutations accumulate disrupting 'normal' mt-driven processes and cellular functions. To investigate the heteroplasmy of such mtDNA changes we developed a moderate throughput mt isolation procedure to quantify the mt single-nucleotide variant (SNV) landscape in individual mouse neurons and astrocytes In this study we amplified mt-genomes from 1,645 single mitochondria (mts) isolated from mouse single astrocytes and neurons to 1. determine the distribution and proportion of mt-SNVs as well as mutation pattern in specific target regions across the mt-genome, 2. assess differences in mtDNA SNVs between neurons and astrocytes, and 3. Study cosegregation of variants in the mouse mtDNA.

Results

1. The data show that specific sites of the mt-genome are permissive to SNV presentation while others appear to be under stringent purifying selection. Nested hierarchical analysis at the levels of mitochondrion, cell, and mouse reveals distinct patterns of inter- and intra-cellular variation for mt-SNVs at different sites. 2. Further, differences in the SNV incidence were observed between mouse neurons and astrocytes for two mt-SNV 9027:G>A and 9419:C>T showing variation in the mutational propensity between these cell types. Purifying selection was observed in neurons as shown by the Ka/Ks statistic, suggesting that neurons are under stronger evolutionary constraint as compared to astrocytes. 3. Intriguingly, these data show strong linkage between the SNV sites at nucleotide positions 9027 and 9461.

Conclusion

This study suggests that segregation as well as clonal expansion of mt-SNVs is specific to individual genomic loci, which is important foundational data in understanding of heteroplasmy and disease thresholds for mutation of pathogenic variants.

Structure and Phylogenetic Relationships of Scolopacidae Mitogenomes (Charadriiformes: Scolopacidae)

The family Scolopacidae presents a valuable subject for evolutionary research; however, molecular studies of Scolopacidae are still relatively understudied, and the phylogenetic relationships of certain species remain unclear. In this study, we sequenced and obtained complete mitochondrial DNA (mtDNA) from Actitis hypoleucos and partial mtDNA from Numenius arquata, Limosa limosa, and Limnodromus semipalmatus. The complete mtDNA contained 13 protein-coding genes (PCGs), two ribosomal RNA genes, 22 tRNA genes, and a control region. Scolopacidae contained three types of start codons and five types of stop codons (including one incomplete stop codon, T--). In 13 protein-coding genes, average uncorrected pairwise distances (Aupd) revealed that ATP8 was the least conserved while COX3 had the lowest evolutionary rate. The ratio of Ka/Ks suggested that all PCGs were under purifying selection. Using two methods (maximum likelihood and Bayesian inference) to analyze the phylogenetic relationships of the family Scolopacidae, it was found that the genera Xenus and Actitis were clustered into another sister group, while the genus Phalaropus is more closely related to the genus Tringa. The genera Limnodromus, Gallinago, and Scolopax form a monophyletic group. This study improves our understanding of the evolutionary patterns and phylogenetic relationships of the family Scolopacidae.

Comparative mitochondrial genomics in Nematoda reveal astonishing variation in compositional biases and substitution rates indicative of multi-level selection

BACKGROUND

Nematodes are the most abundant and diverse metazoans on Earth, and are known to significantly affect ecosystem functioning. A better understanding of their biology and ecology, including potential adaptations to diverse habitats and lifestyles, is key to understanding their response to global change scenarios. Mitochondrial genomes offer high species level characterization, low cost of sequencing, and an ease of data handling that can provide insights into nematode evolutionary pressures.

RESULTS

Generally, nematode mitochondrial genomes exhibited similar structural characteristics (e.g., gene size and GC content), but displayed remarkable variability around these general patterns. Compositional strand biases showed strong codon position specific G skews and relationships with nematode life traits (especially parasitic feeding habits) equal to or greater than with predicted phylogeny. On average, nematode mitochondrial genomes showed low non-synonymous substitution rates, but also high clade specific deviations from these means. Despite the presence of significant mutational saturation, non-synonymous (dN) and synonymous (dS) substitution rates could still be significantly explained by feeding habit and/or habitat. Low ratios of dN:dS rates, particularly associated with the parasitic lifestyles, suggested the presence of strong purifying selection.

CONCLUSIONS

Nematode mitochondrial genomes demonstrated a capacity to accumulate diversity in composition, structure, and content while still maintaining functional genes. Moreover, they demonstrated a capacity for rapid evolutionary change pointing to a potential interaction between multi-level selection pressures and rapid evolution. In conclusion, this study helps establish a background for our understanding of the potential evolutionary pressures shaping nematode mitochondrial genomes, while outlining likely routes of future inquiry.

Complete mitogenome and intra-family comparative mitogenomics showed distinct position of Pama Croaker Otolithoides pama

The Pama Croaker, Otolithoides pama, is an economically important fish species in Bangladesh. Intra-family similarities in morphology and typical barcode sequences of cox1 create ambiguities in its identification. Therefore, morphology and the complete mitochondrial genome of O. pama, and comparative mitogenomics within the family Sciaenidae have been studied. Extracted genomic DNA was subjected to Illumina-based short read sequencing for De-Novo mitogenome assembly. The complete mitogenome of O. pama (Accession: OQ784575.1) was 16,513 bp, with strong AC biasness and strand asymmetry. Relative synonymous codon usage (RSCU) among 13 protein-coding genes (PCGs) of O. pama was also analyzed. The studied mitogenomes including O. pama exhibited consistent sizes and gene orders, except for the genus Johnius which possessed notably longer mitogenomes with unique gene rearrangements. Different genetic distance metrics across 30 species of Sciaenidae family demonstrated 12S rRNA and the control region (CR) as the most conserved and variable regions, respectively, while most of the PCGs undergone a purifying selection. Different phylogenetic trees were congruent with one another, where O. pama was distinctly placed. This study would contribute to distinguishing closely related fish species of Sciaenidae family and can be instrumental in conserving the genetic diversity of O. pama.

Evolutionary and phylogenetic insights from the mitochondrial genomic analysis of Diceraeus melacanthus and D. furcatus (Hemiptera: Pentatomidae)

The mitochondrial genomes of D. melacanthus and D. furcatus were sequenced and used to investigate the phylogenetic relationships with 54 species of Pentatomidae. Their mitogenomes are 17,197 and 15,444 bp-long, respectively, including 13 protein-coding genes (PCGs), 2 ribosomal RNA genes, and 22/21 transfer RNA genes, with conserved gene arrangement. Leu, Lys, and Ser were the most common amino acids in their PCGs. PCGs evolutionary analysis indicated their mitogenomes are under purifying selection, and the most conserved genes are from the cytochrome complex, reinforcing their suitability as markers for molecular taxonomy. We identified 490 mtSSRs in 56 Pentatomidae species, with large variation and a positive correlation between mtSSR number and genome size. Three mtSSRs were identified in each Diceraeus species. Only the mtSSR in the nad6 (D. melacanthus) and nad4 (D. furcatus) appear to have application as molecular markers for species characterization. Phylogenetic analysis confirmed the monophyly of Pentatomidae. However, our analysis challenged the monophyly of Pentatominae and Podopinae. We also detected unexpected relationships among some tribes and genera, highlighting the complexity of the internal taxonomic structure of Pentatomidae. Both Diceraeus species were grouped in the same clade with the remaining Carpocorini analyzed.

Disruption of recombination machinery alters the mutational landscape in plant organellar genomes

Land plant organellar genomes have extremely low rates of point mutation yet also experience high rates of recombination and genome instability. Characterizing the molecular machinery responsible for these patterns is critical for understanding the evolution of these genomes. While much progress has been made towards understanding recombination activity in land plant organellar genomes, the relationship between recombination pathways and point mutation rates remains uncertain. The organellar targeted mutS homolog MSH1 has previously been shown to suppress point mutations as well as non-allelic recombination between short repeats in Arabidopsis thaliana. We therefore implemented high-fidelity Duplex Sequencing to test if other genes that function in recombination and maintenance of genome stability also affect point mutation rates. We found small to moderate increases in the frequency of single nucleotide variants (SNVs) and indels in mitochondrial and/or plastid genomes of A. thaliana mutant lines lacking radA, recA1, or recA3. In contrast, osb2 and why2 mutants did not exhibit an increase in point mutations compared to wild type (WT) controls. In addition, we analyzed the distribution of SNVs in previously generated Duplex Sequencing data from A. thaliana organellar genomes and found unexpected strand asymmetries and large effects of flanking nucleotides on mutation rates in WT plants and msh1 mutants. Finally, using long-read Oxford Nanopore sequencing, we characterized structural variants in organellar genomes of the mutant lines and show that different short repeat sequences become recombinationally active in different mutant backgrounds. Together, these complementary sequencing approaches shed light on how recombination may impact the extraordinarily low point mutation rates in plant organellar genomes.

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.

Integrative phylogenetic, phylogeographic and morphological characterisation of the Unio crassus species complex reveals cryptic diversity with important conservation implications

The global decline of freshwater mussels and their crucial ecological services highlight the need to understand their phylogeny, phylogeography and patterns of genetic diversity to guide conservation efforts. Such knowledge is urgently needed for Unio crassus, a highly imperilled species originally widespread throughout Europe and southwest Asia. Recent studies have resurrected several species from synonymy based on mitochondrial data, revealing U. crassus to be a complex of cryptic species. To address long-standing taxonomic uncertainties hindering effective conservation, we integrate morphometric, phylogenetic, and phylogeographic analyses to examine species diversity within the U. crassus complex across its entire range. Phylogenetic analyses were performed using cytochrome c oxidase subunit I (815 specimens from 182 populations) and, for selected specimens, whole mitogenome sequences and Anchored Hybrid Enrichment (AHE) data on ∼ 600 nuclear loci. Mito-nuclear discordance was detected, consistent with mitochondrial DNA gene flow between some species during the Pliocene and Pleistocene. Fossil-calibrated phylogenies based on AHE data support a Mediterranean origin for the U. crassus complex in the Early Miocene. The results of our integrative approach support 12 species in the group: the previously recognised Unio bruguierianus, Unio carneus, Unio crassus, Unio damascensis, Unio ionicus, Unio sesirmensis, and Unio tumidiformis, and the reinstatement of five nominal taxa: Unio desectusstat. rev., Unio gontieriistat. rev., Unio mardinensisstat. rev., Unio nanusstat. rev., and Unio vicariusstat. rev. Morphometric analyses of shell contours reveal important morphospace overlaps among these species, highlighting cryptic, but geographically structured, diversity. The distribution, taxonomy, phylogeography, and conservation of each species are succinctly described.

Strand-preferred base editing of organellar and nuclear genomes using CyDENT

Transcription-activator-like effector (TALE)-based tools for base editing of nuclear and organellar DNA rely on double-stranded DNA deaminases, which edit substrate bases on both strands of DNA, reducing editing precision. Here, we present CyDENT base editing, a CRISPR-free, strand-selective, modular base editor. CyDENT comprises a pair of TALEs fused with a FokI nickase, a single-strand-specific cytidine deaminase and an exonuclease to generate a single-stranded DNA substrate for deamination. We demonstrate effective base editing in nuclear, mitochondrial and chloroplast genomes. At certain mitochondrial sites, we show editing efficiencies of 14% and strand specificity of 95%. Furthermore, by exchanging the CyDENT deaminase with one that prefers editing GC motifs, we demonstrate up to 20% mitochondrial base editing at sites that are otherwise inaccessible to editing by other methods. The modular nature of CyDENT enables a suite of bespoke base editors for various applications.

Method for depletion of mitochondria DNA in human bronchial epithelial cells

Mitochondria are increasingly recognized to play a role in the airway inflammation of asthma. Model systems to study the role of mitochondrial gene expression in bronchial epithelium are lacking. Here, we create custom bronchial epithelial cell lines that are depleted of mitochondrial DNA. One week of ethidium bromide (EtBr) treatment led to ∼95 % reduction of mtDNA copy number (mtDNA-CN) in cells, which was further reduced by addition of 25 µM 2',3'-dideoxycytidin (ddC). Treatment for up to three weeks with EtBr and ddC led to near complete loss of mtDNA. The basal oxygen consumption rate (OCR) of mtDNA-depleted BET-1A and BEAS-2B cells dropped to near zero. Glycolysis measured by extracellular acidification rate (ECAR) increased ∼two-fold in cells when mtDNA was eliminated. BET-1A ρ0 and BEAS-2B ρ0 cells were cultured for two months, frozen and thawed, cultured for two more months, and maintained near zero mtDNA-CN. Mitochondrial DNA-depleted BET-1A ρ0 and BEAS-2B ρ0 cell lines are viable, lack the capacity for aerobic respiration, and increase glycolysis.•BET-1A and BEAS-2B cells were treated with ethidium bromide (EtBr) with or without 2',3'-dideoxycytidine (ddC) to create cells lacking mitochondrial DNA (mtDNA).•Cells' mtDNA copy number relative to nuclear DNA (nDNA) were verified by quantitative polymerase chain reaction (qPCR).•Cells were also assessed for oxidative phosphorylation by measures of oxygen consumption using the Seahorse analyzer.

Characterization and Molecular Phylogenetic Analysis of Subfamily Erebinae (Lepidoptera: Noctuoidea: Erebidae) Using Five Complete Mitochondrial Genomes

In this study, the complete mitogenomes of Sympis rufibasis, Lacera noctilio, Oxyodes scrobiculata, Mocis undata, and Artena dotata were newly sequenced to bring up-to-date the database using the next-generation sequencing methods. The gene order of all sequenced mitogenomes was identical consisting of 13 protein-coding genes (PCGs), 22 transfer RNA genes (tRNAs), two ribosomal RNA genes (rRNAs), and a non-coding A+T-rich region, which were common to other Lepidopteran insects. All protein-coding genes (PCGs) initiated with a canonical ATN codon and ended with TAN or an incomplete stop codon, single T. The A+T-rich region of S. rufibasis, L. noctilio, O. scrobiculata, M. undata, and A. dotata are 406 bp, 462 bp, 372 bp, 410 bp, and 406 bp long, respectively, containing number of characteristics that are distinctive to Noctuoidea moths. We analyzed concatenated amino acid sequences of protein-coding genes not including rRNAs, using Maximum Likelihood and Bayesian Inference methods. The phylogenetic analyses indicated that the tribe relationships within Erebinae were reconstructed as (Sypnini+((Erebini 1+Poaphilini 1)+((Euclidiini+Catocalini+(Hypopyrini+Erebini 2))+((Hulodini+(Poaphilini 2+Ophiusini))))). Phylogenetic analyses supported and confirmed the monophyly of the subfamilies' relationships as follows: (Hypeninae+Lymantriinae)+((Scoliopterginae+((Calpinae+Erebinae)+((Herminiinae+Aganainae)+Arctiinae)))) within Erebidae.

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.

Complete Mitochondrial Genomes of Nannostomus Pencilfish: Genome Characterization and Phylogenetic Analysis

Although the pencilfish is a globally popular economic fish in the aquarium market, its taxonomic classification could be further refined. In order to understand the taxonomy of species of the genus Nannostomus (Characiformes, Lebiasinidae) and their phylogenetic position within the order Characiformes, in this study, we characterized mitochondrial genomes (mitogenomes) from four Nannostomus species for the first time. The four mitogenomes exhibited the typical circular structure, with overall sizes varying from 16,661 bp to 16,690 bp. They contained 13 protein-coding genes (PCGs), 2 ribosomal RNA genes (rRNAs), 22 transfer RNA genes (tRNAs), and 1 control region (CR). Nucleotide composition analysis suggested that the mitochondrial sequences were biased toward A and T. Bayesian inference and maximum likelihood analyses based on PCGs support the family Lebiasinidae classification, described using four Nannostomus species, clustering together with Lebiasina multimaculata from the same family. The results of this study support the current taxonomic classification of the family Lebiasinidae. Phylogenetic analysis also suggested that gene rearrangement would not significantly impact the phylogenetic relationships within the order Characiformes. These results might provide new data regarding the phylogeny and classification of the order Characiformes, thus providing a theoretical basis for the economic development of aquarium fish markets.

Characteristics and Comparative Analysis of Mitochondrial Genomes of the Aphid Genus Hyalopterus Koch (Hemiptera: Aphididae: Aphidinae)

Using Illumina sequencing technology, we generated complete mitochondrial genomes (mitogenomes) of three constituent species of the aphid genus Hyalopterus Koch, Hyalopterus amygdali (Blanchard), Hyalopterus arundiniformis Ghulamullah, and Hyalopterus pruni (Geoffroy). The sizes of the Hyalopterus mitogenomes range from 15,306 to 15,410 bp, primarily due to variations in the length of non-coding regions. The Hyalopterus mitogenomes consist of 37 coding genes arranged in the order of the ancestral insect mitogenome, a control region, and a repeat region between trnE and trnF. According to the COI-based analysis, one previously reported mitogenome of H. pruni should be assigned to H. arundiniformis. The gene order, nucleotide composition, and codon usage in the Hyalopterus mitogenomes are highly conserved and similar to those of other species of Aphidinae. The tandem repeat units differ in nucleotide composition, length, and copy number across three Hyalopterus species. Within the widespread Eurasian species H. arundiniformis, variation in repeat units among different geographic populations is observed, indicating that the repeat region may provide valuable insights for studying the intraspecific diversification of aphids. Phylogenetic analyses based on 28 complete mitogenomes of Aphidinae supported the monophyly of Aphidinae, Aphidini, Macrosiphini, and two subtribes of Aphidini. Hyalopterus was monophyletic. H. amygdali and H. pruni formed a sister group, while H. arundiniformis was placed basally. Characterization of the mitogenomes of Hyalopterus provides valuable resources for further comparative studies and for advancing our understanding of the aphid mitogenome architecture.