Animal mitochondrial genomes

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.

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.

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.

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

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.

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.

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.

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.

Sequence comparison of the mitochondrial genomes of Plesionika species (Caridea: Pandalidae), gene rearrangement and phylogenetic relationships of Caridea

Background

Despite the Caridean shrimps' vast species richness and ecological diversity, controversies persist in their molecular classification. Within Caridea, the Pandalidae family exemplifies significant taxonomic diversity. As of June 25, 2023, GenBank hosts only nine complete mitochondrial genomes (mitogenomes) for this family. The Plesionika genus within Pandalidae is recognized as polyphyletic. To improve our understanding of the mitogenome evolution and phylogenetic relationships of Caridea, this study introduces three novel mitogenome sequences from the Plesionika genus: P.  ortmanni, P. izumiae and P. lophotes.

Methods

The complete mitochondrial genomes of three Plesionika species were sequenced utilizing Illumina's next-generation sequencing (NGS) technology. After assembling and annotating the mitogenomes, we conducted structural analyses to examine circular maps, sequence structure characteristics, base composition, amino acid content, and synonymous codon usage frequency. Additionally, phylogenetic analysis was performed by integrating existing mitogenome sequences of true shrimp available in GenBank.

Results

The complete mitogenomes of the three Plesionika species encompass 37 canonical genes, comprising 13 protein-coding genes (PCGs), 22 transfer RNAs (tRNAs), two ribosomal RNAs (rRNAs), and one control region (CR). The lengths of these mitogenomes are as follows: 15,908 bp for P. ortmanni, 16,074 bp for P. izumiae and 15,933 bp for P. lophotes. Our analyses extended to their genomic features and structural functions, detailing base composition, gene arrangement, and codon usage. Additionally, we performed selection pressure analysis on the PCGs of all Pandalidae species available in Genbank, indicating evolutionary purification selection acted on the PCGs across Pandalidae species. Compared with the ancestral Caridea, translocation of two tRNA genes, i.e., trnP or trnT, were found in the two newly sequenced Plesionika species-P. izumiae and P. lophotes. We constructed a phylogenetic tree of Caridea using the sequences of 13 PCGs in mitogenomes. The results revealed that family Pandalidae exhibited robust monophyly, while genus Plesionika appeared to be a polyphyletic group.

Conclusions

Gene rearrangements within the Pandalidae family were observed for the first time. Furthermore, a significant correlation was discovered between phylogenetics of the Caridea clade and arrangement of mitochondrial genes. Our findings offer a detailed exploration of Plesionika mitogenomes, laying a crucial groundwork for subsequent investigations into genetic diversity, phylogenetic evolution, and selective breeding within this genus.

Complete Mitochondrial Genomes and Phylogenetic Analysis of Genus Henricia (Asteroidea: Spinulosida: Echinasteridae)

The genus Henricia is known to have intraspecific morphological variations, making species identification difficult. Therefore, molecular phylogeny analysis based on genetic characteristics is valuable for species identification. We present complete mitochondrial genomic sequences of Henricia longispina aleutica, H. reniossa, and H. sanguinolenta for the first time in this study. This study will make a significant contribution to our understanding of Henricia species and its relationships within the class Asteroidea. Lengths of mitochondrial genomes of the three species are 16,217, 16,223, and 16,194 bp, respectively, with a circular form. These genomes contained 13 protein-coding genes, two ribosomal RNA genes, 22 transfer RNA genes, and a D-loop. The gene order and direction aligned with other asteroid species. Phylogenetic relationship analysis showed that our Henricia species were in a monophyletic clade with other Henricia species and in a large clade with species (Echinaster brasiliensis) from the same family. These findings provide valuable insight into understanding the phylogenetic relationships of species in the genus Henricia.

Comparative Mitogenome of Phylogenetic Relationships and Divergence Time Analysis within Potamanthidae (Insecta: Ephemeroptera)

Potamanthidae belongs to the superfamily Ephemeroidea but has no complete mt genome released in the NCBI (except for two unchecked and one partial mt genome). Since the sister clade to Potamanthidae has always been controversial, we sequenced seven mt genomes of Potamanthidae (two species from Rhoenanthus and five species from Potamanthus) in order to rebuild the phylogenetic relationships of Potamanthidae in this study. The divergence time of Potamanthidae was also investigated by utilizing five fossil calibration points because of the indeterminate origin time. In addition, because Rhoenanthus coreanus and Potamanthus luteus are always in low-temperature environments, we aimed to explore whether these two species were under positive selection at the mt genome level. Amongst the 13 PCGs, CGA was used as the start codon in COX1, whereas other genes conformed to initiating with an ATN start codon. From this analysis, UUA (L), AUU (I), and UUU (F) had the highest usage. Furthermore, the DHU arm was absent in the secondary structure of S1 in all species. By combining the 13 PCGs and 2 rRNAs, we reconstructed the phylogenetic relationship of Potamanthidae within Ephemeroptera. The monophyly of Potamanthidae and the monophyly of Rhoenanthus and Potamanthus were supported in the results. The phylogenetic relationship of Potamanthidae + (Ephemeridae + Polymitarcyidae) was also recovered with a high prior probability. The divergence times of Potamanthidae were traced to be 90.44 Mya (95% HPD, 62.80-121.74 Mya), and the divergence times of Rhoenanthus and Potamanthus originated at approximately 64.77 Mya (95% HPD, 43.82-88.68 Mya), thus belonging to the late Pliocene Epoch or early Miocene Epoch. In addition, the data indicated that R. coreanus was under negative selection and that ATP8 and ND2 in Potamanthidae had a high evolutionary rate.

Characterization of the complete mitochondrial genome of an endemic species in China, Aulocera merlina (Lepidoptera: Nymphalidae: Satyrinae) and phylogenetic analysis within Satyrinae

The mitochondrial genome (mitogenome) has been extensively used as molecular markers in determining the insect phylogenetic relationships. In order to resolve the relationships among tribes and subtribes of Satyrinae at the mitochondrial genomic level, we obtained the complete mitogenome of Aulocera merlina (Oberthür, 1890) (Lepidoptera: Nymphalidae: Satyrinae) with a size of 15,259 bp. The mitogenome consisted of 37 typical genes, including 13 protein-coding genes (PCGs), 2 ribosomal RNA genes (rRNAs), 22 transfer RNA genes (tRNAs), and an A + T-rich region. The gene organization and arrangement were similar to those of all other known Satyrinae mitogenomes. All PCGs were initiated with the canonical codon pattern ATN, except for the cox1 gene, which used an atypical CGA codon. Nine PCGs used the complete stop codon TAA, while the remaining PCGs (cox1, cox2, nad4, and nad5) were terminated with a single T nucleotide. The canonical cloverleaf secondary structures were found in all tRNAs, except for trnS1 which lacked a dihydrouridine arm. The 448 bp A + T-rich region was located between rrnS and trnM, and it included the motif ATAGA followed by a 19-bp poly-T stretch and a microsatellite-like (TA)6 element preceded by the ATTTA motif. The phylogenetic tree, inferred using Bayesian inference and maximum likelihood methods, generated similar tree topologies, revealing well-supported monophyletic groups at the tribe level and recovering the relationship ((Satyrini + Melanitini) + ((Amathusiini + Elymniini) + Zetherini)). The close relationship between Satyrina and Melanargiina within the Satyrini was widely accepted. Additionally, Lethina, Parargina, and Mycalesina were closely related and collectively formed a sister group to Coenonymphina. Moreover, A. merlina was closely related to Oeneis buddha within the Satyrina. These findings will provide valuable information for future studies aiming to elucidate the phylogenetic relationships of Satyrinae.

First report of the complete mitochondrial genome of Carpomya pardalina (Bigot) (Diptera: Tephritidae) and phylogenetic relationships with other Tephritidae

Carpomya pardalina is known for its potential invasiveness, which poses a significant and alarming threat to Cucurbitaceae crops. It is considered a highly perilous pest species that requires immediate attention for quarantine and prevention. Due to the challenges in distinguishing pests of the Tephritidae family based on morphological characteristics, it is imperative to elucidate the mitochondrial genomic information of C. pardalina. In this study, the mitochondrial genome sequence of C. pardalina was determined and analyzed using next-generation sequencing. The results revealed that the mitogenome sequence had a total length of 16,257 bp, representing a typical circular molecule. It consisted of 13 PCGs, two rRNA genes, 22 tRNA genes and a non-coding region. The structure and organization of the mitochondrial genome of C. pardalina were found to be typical and similar to the published homologous sequences of other fruit flies in the Tephritidae family. Phylogenetic analysis confirmed that C. pardalina belongs to the Carpomya genus, which is consistent with traditional morphological taxonomy. Additionally, Carpomya and Rhagoletis were identified as sister groups. This study presents the first report of the complete mitochondrial genome of C. pardalina, which can serve as a valuable resource for future investigations in species diagnosis, evolutionary biology, prevention and control measures.

The characterization of the mitochondrial genome of Graptemys ouachitensis

Graptemys ouachitensis (CAGLE, 1953) belongs to the Graptemys genus, the Emydidae family, and the Testudines order. This study involved sequencing the complete mitochondrial genome (mitogenome) of G. ouachitensis using next-generation sequencing, and analyzing the essential characteristics, and phylogenetic relationship. The results revealed that the G. ouachitensis mitogenome was 16,674 bp in length (A: 34.1%, C: 26.0%, G: 13.0%, T: 26.9%) and included 22 tRNAs, 13 protein-coding genes, two ribosomal RNA genes, and a non-coding control region (GenBank accession: NC071766). The genome composition of G. ouachitensis presented a slight A + T bias (61.0%) and exhibited a positive AT skew (0.118) and a negative GC skew (-0.333). A phylogenetic analysis based on the complete mitogenome indicated that the G. ouachitensis was more closely associated with Malaclemys terrapin than the other eight known Emydidae species. Thus, our findings present a novel mitogenome at the species level. This study introduces the first complete mitogenome of G. ouachitensis, providing valuable molecular information for phylogenetic and conservation genetics analyses of G. ouachitensis.

Complete Mitochondrial Genome of Four Peristediidae Fish Species: Genome Characterization and Phylogenetic Analysis

The systematic revision of the family Peristediidae remains an unresolved issue due to their diverse and unique morphology. Despite the popularity of using mitochondrial genome research to comprehensively understand phylogenetic relationships in fish, genetic data for peristediid fish need to be included. Therefore, this study aims to investigate the mitochondrial genomic characteristics and intra-family phylogenetic relationships of Peristediidae by utilizing mitochondrial genome analysis. Therefore, this study aims to investigate the phylogenetic relationship of Peristediidae by utilizing mitochondrial genome analysis. The mitochondrial genome of four species of Peristediidae (Peristedion liorhynchus, Satyrichthys welchi, Satyrichthys rieffeli, and Scalicus amiscus) collected in the East China Sea was studied. The mitochondrial gene sequence lengths of four fish species were 16,533 bp, 16,526 bp, 16,527 bp, and 16,526 bp, respectively. They had the same mitochondrial structure and were all composed of 37 genes and one control region. Most PCGs used ATG as the start codon, and a few used GTG as the start codon. An incomplete stop codon (TA/T) occurred. The AT-skew and GC-skew values of 13 PCGs from four species were negative, and the GC-skew amplitude was greater than that of AT-skew. All cases of D-arm were found in tRNA-Ser (GCT). The Ka/Ks ratio analysis indicated that 13 PCGs were suffering purifying selection. Based on 12 PCGs (excluding ND6) sequences, a phylogenetic tree was constructed using Bayesian inference (BI) and maximum likelihood (ML) methods, providing a further supplement to the scientific classification of Peristediidae fish. According to the results of divergence time, the four species of fish had apparent divergence in the Early Cenozoic, which indicates that the geological events at that time caused the climax of species divergence and evolution.

Complete Mitochondrial Genome for Lucilia cuprina dorsalis (Diptera: Calliphoridae) from the Northern Territory, Australia

The Australian sheep blowfly, Lucilia cuprina dorsalis, is a major sheep ectoparasite causing subcutaneous myiasis (flystrike), which can lead to reduced livestock productivity and, in severe instances, death of the affected animals. It is also a primary colonizer of carrion, an efficient pollinator, and used in maggot debridement therapy and forensic investigations. In this study, we report the complete mitochondrial (mt) genome of L. c. dorsalis from the Northern Territory (NT), Australia, where sheep are prohibited animals, unlike the rest of Australia. The mt genome is 15,943 bp in length, comprising 13 protein-coding genes (PCGs), two ribosomal RNAs (rRNAs), 22 transfer RNAs (tRNAs), and a non-coding control region. The gene order of the current mt genome is consistent with the previously published L. cuprina mt genomes. Nucleotide composition revealed an AT bias, accounting for 77.5% of total mt genome nucleotides. Phylogenetic analyses of 56 species/taxa of dipterans indicated that L. c. dorsalis and L. sericata are the closest among all sibling species of the genus Lucilia, which helps to explain species evolution within the family Luciliinae. This study provides the first complete mt genome sequence for L. c. dorsalis derived from the NT, Australia to facilitate species identification and the examination of the evolutionary history of these blowflies.

Phylomitogenomics of two Neotropical species of long-legged crickets Endecous Saussure, 1878 (Orthoptera: Phalangopsidae)

Mitochondrial genomes have provided significant insights into the evolution of several insects. A typical mitogenome contains 37 genes, and variations in gene order can indicate evolutionary relationships between species. In this study, we have assembled the first complete mitogenomes of Endecous chape and E. onthophagus and analyzed the phylogenetic implications for the Gryllidea infraorder. We performed DNA extractions and genome sequencing for both Endecous species. Subsequently, we searched for raw data in the Sequence Read Archive (SRA) in NCBI. Using the SRA data, we assembled the partial mitogenome of Dianemobius nigrofasciatus and annotated the protein-coding genes (PCGs) for nine species. Phylogenomic relationships were reconstructed using Maximum Likelihood (ML) and Bayesian Inference (BI), utilizing the PCGs from 49 Gryllidea species. The mitogenome lengths of E. chape and E. onthophagus are 16,266 bp and 16,023 bp, respectively, while D. nigrofasciatus has a length of 15,359 bp. Our results indicate that species within the infraorder exhibit four types of gene order arrangements that align with their phylogenetic relationships. Both phylogenomic trees displayed strong support, and the ML corroborated with the literature. Gryllidea species have significantly contributed to various fields, and studying their mitogenomes can provide valuable insights into this infraorder evolution.

Evaluating the mitochondrial genomic diversity, global distribution and niche overlap of two invasive Phthorimaea species

This study sought to evaluate the genetic diversity of two invasive Phthorimaea species (Phthorimaea operculella Zeller and Phthorimaea absoluta Meyrick), and identify potential niche overlap of both species. The complete mitogenomes of P. operculella and P. absoluta were sequenced and compared. Furthermore, the diversity within the family Gelechiidae was assessed. Subsequently, two species distribution models (MaxEnt and BIOCLIM) were used to identify niche overlaps of both species globally. The complete mitogenomes of both species were similar in size and structure, with a pairwise identity of 92.3%. The models showed a niche overlap of both species and revealed areas of marginal to high suitability for both pests in countries where they have not been reported. Therefore, these results present a compelling case for a deeper genetic and ecological investigation of the Gelechiidae family for taxonomic harmonization, an early warning for surveillance, stricter phytosanitary considerations and preventive management against the spread of the pests.

Codon usage bias analysis in the mitochondrial genomes of five Rhingia Scopoli (Diptera, Syrphidae, Eristalinae) species

This study presents the sequencing and annotation of mitochondrial genomes from five Rhingia species of the family Syrphidae, focusing on codon bias. Each species possessed 22 tRNAs genes, 13 protein-coding genes, 2 rRNAs genes, and a control region, without any observed gene rearrangements. Nucleotide composition analysis revealed a higher AT content compared with GC content, indicating AT enrichment. Neutrality plot, Parity rule 2 bias, and effective number of codons plot analyses collectively indicated that natural selection primarily influences the codon usage bias in the five Rhingia species. Relative synonymous codon usage analysis identified the optimal codons for Rhingia binotata, R. fromosana, R. campestris, R. louguanensis, and R. xanthopoda as 10, 14, 10, 11, and 12, respectively, all ending with A/U and exhibiting AT preference. Phylogenetic analysis, based on maximum likelihood and Bayesian inference methods applied to three datasets, confirmed the monophyly of Rhingia. In conclusion, this research establishes a foundation for understanding the phylogenetic evolution and codon usage patterns in Rhingia, offering valuable for future studies.

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

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

Mitogenomic Analysis and Phylogenetic Implications for the Deltocephaline Tribe Chiasmini (Hemiptera: Cicadellidae: Deltocephalinae)

The grassland leafhopper tribe Chiasmini (Cicadellidae: Deltocephalinae) presently comprises 324 described species worldwide, with the highest species diversity occurring in the Nearctic region but a greater diversity of genera occurring in the Old World. In China, this tribe comprises 39 described species in 11 genera, but the fauna remains understudied. The complete mitogenomes of three species of this tribe have been sequenced previously. In order to better understand the phylogenetic position of Chiasmini within the subfamily Deltocephalinae and to investigate relationships among Chiasmini genera and species, we sequenced and analyzed the complete mitogenomes of 13 species belonging to seven genera from China. Comparison of the newly sequenced mitogenomes reveals a closed circular double-stranded structure containing 37 genes with a total length of 14,805 to 16,269 bp and a variable number of non-coding A + T-rich regions. The gene size, gene order, gene arrangement, base composition, codon usage, and secondary structure of tRNAs of the newly sequenced mitogenomes of these 13 species are highly conserved in Chiasmini. The ATN codon is commonly used as the start codon in protein-coding genes (PCGs), except for ND5 in Doratura sp. and ATP6 in Nephotettix nigropictus, which use the rare GTG start codon. Most protein-coding genes have TAA or TAG as the stop codon, but some genes have an incomplete T stop codon. Except for the tRNA for serine (trnS1(AGN)), the secondary structure of the other 21 tRNAs is a typical cloverleaf structure. In addition to the primary type of G-U mismatch, five other types of tRNA mismatches were observed: A-A, A-C, A-G, U-C, and U-U. Chiasmini mitochondrial genomes exhibit gene overlaps with three relatively stable regions: the overlapping sequence between trnW and trnC is AAGTCTTA, the overlapping sequence between ATP8 and ATP6 is generally ATGATTA, and the overlapping sequence between ND4 and ND4L is generally TTATCAT. The largest non-coding region is the control region, which exhibits significant length and compositional variation among species. Some Chiasmini have tandem repeat structures within their control regions. Unlike some other deltocephaline leafhoppers, the sequenced Chiasmini lack mitochondrial gene rearrangements. Phylogenetic analyses of different combinations of protein-coding and ribosomal genes using maximum likelihood and Bayesian methods under different models, using either amino acid or nucleotide sequences, are generally consistent and also agree with results of prior analyses of nuclear and partial mitochondrial gene sequence data, indicating that complete mitochondrial genomes are phylogenetically informative at different levels of divergence within Chiasmini and among leafhoppers in general. Apart from Athysanini and Opsiini, most of the deltocephaline tribes are recovered as monophyletic. The results of ML and BI analyses show that Chiasmini is a monophyletic group with seven monophyletic genera arranged as follows: ((Zahniserius + (Gurawa + (Doratura + Aconurella))) + (Leofa + (Exitianus + Nephotettix))).

Characteristics of the complete mitochondrial genome of Gerres limbatus (Cuvier, 1830) (Perciformes: Gerreidae)

The saddleback silver-biddy Gerres limbatus (Cuvier 1830) is distributed in Indo-West Pacific Oceans and associated with shallow coastal marine waters and estuaries. In this study, the complete mitochondrial genome of G. limbatus was firstly documented, which is 16,730 bp in length, including 13 protein-coding genes, 22 transfer RNA genes and 2 ribosomal RNA genes. The overall base composition of the mitochondrial genome is 26.42% A, 28.68% C, 27.32% T, and 17.58% G. The Maximum Likelihood phylogenetic tree was constructed based on COI gene of the 31 species from the family Gerreidae, with Heteroclinus puellarum and Hypopterus macropterus as outgroups. It revealed that G. erythrourus was placed as the sister group to G. limbatus.

Mitochondrial genome comparison and phylogenetic position of Fannia pusio among the Calyptratae flies

Fannia pusio, the chicken dung fly species, remains unexplored despite its forensic, sanitary, and veterinary importance in the Nearctic and Neotropical regions. In this study, we obtained the complete mitochondrial genome of Fannia pusio for the first time using next-generation sequencing. We compared it with previously published mitogenomes of the genus from the Palearctic region, and its phylogenetic position was studied based on the concatenated protein-coding genes (PCGs) dataset of Calyptratae flies. The circular mitochondrial genome of F. pusio is 16,176 bp in length, with a high A + T content (78.3%), whose gene synteny, codon usage analysis, and amino acid frequency are similar to previously reported Fannia mitogenomes. All PCGs underwent purifying selection except the nad2 gene. Interspecific K2P distances of PCGs of Fannia yielded an average of 12.4% (8.1%-21.1%). The Fannia genus is monophyletic and closely related to Muscidae based on molecular data. Further taxonomic sampling is required to deep into the phylogenetic relationships of the originally proposed species-groups and subgroups within the genus. These results provide a valuable dataset for studying the mitochondrial genome evolution and a resource for the taxonomy and systematics of Fannia.

Describe the morphology and mitochondrial genome of Mecidea indica Dallas, 1851 (Hemiptera, Pentatomidae), with its phylogenetic position

We here describe the external morphology and complete mitochondrial genome characteristics of Mecidea indica Dallas, 1851, and clarify the evolutionary rate and divergence time. The M. indica mitochondrial genome length is 15,670 bp, and it exhibits a typical high A+T-skew (76.31%). The sequence shows strong synteny with the original gene arrangement of Drosophila yakuba Burla, 1954 without rearrangement. The M. indica mitochondrial genome characteristics were analyzed, and phylogenetic trees of Pentatomidae were reconstructed using Bayesian methods based on different datasets of the mitochondrial genome datasets. Phylogenetic analysis shows that M. indica belongs to Pentaotominae and form a sister-group with Anaxilaus musgravei Gross, 1976, and Asopinae is highly supported as monophyletic. Molecular clock analysis estimates a divergence time of Pentatomidae of 122.75 Mya (95% HPD: 98.76-145.43 Mya), within the Mesozoic Cretaceous; the divergence time of M. indica and A. musgravii was no later than 50.50 Mya (95% HPD: 37.20-64.80 Mya). In addition, the divergence time of Asopinae was 62.32 Mya (95% HPD: 47.08-78.23 Mya), which was in the Paleogene of the Cenozoic era. This study is of great significance for reconstructing the phylogeny of Pentatomidae and providing insights into its evolutionary history.

The complete mitochondrial genome of Leucoptera coffeella (Lepidoptera: Lyonetiidae) and phylogenetic relationships within the Yponomeutoidea superfamily

The coffee leaf miner (Leucoptera coffeella) is one of the major pests of coffee crops in the neotropical regions, and causes major economic losses. Few molecular data are available to identify this pest and advances in the knowledge of the genome of L. coffeella will contribute to improving pest identification and also clarify taxonomy of this microlepidoptera. L. coffeella DNA was extracted and sequenced using PacBio HiFi technology. Here we report the complete L. coffeella circular mitochondrial genome (16,407 bp) assembled using Aladin software. We found a total of 37 genes, including 13 protein-coding genes (PCGs), 22 transfer RNA genes (tRNAs), 2 ribosomal RNA genes (rRNAs) and an A + T rich-region and a D-loop. The L. coffeella mitochondrial gene organization is highly conserved with similarities to lepidopteran mitochondrial gene rearrangements (trnM-trnI-trnQ). We concatenated the 13 PCG to construct a phylogenetic tree and inferred the relationship between L. coffeella and other lepidopteran species. L. coffeella is found in the Lyonetiidae clade together with L. malifoliella and Lyonetia clerkella, both leaf miners. Interestingly, this clade is assigned in the Yponomeutoidea superfamily together with Gracillariidae, and both superfamilies displayed species with leaf-mining feeding habits.

The complete mitochondrial genome of Chibiraga houshuaii (Lepidoptera, Limacodidae) and its phylogenetic implications

Chibiraga is a mall East Asian genus in the family Limacodidae (slug-moths). The latter includes many agricultural pests. Mitochondrial genome analysis is an important tool for studying insect molecular identification and phylogenetics. However, there are very few mitogenome sequences available for Limacodidae species, and none for the genus Chibiraga at all. To explore the mitogenome features of Chibiraga and verify its phylogenetic position, the complete mitogenome of Chibiraga houshuaii was sequenced and annotated. The complete 15,487 bp genome encoded 37 mitochondrial genes, including 13 protein-coding genes (PCGs), 22 transfer RNA (tRNA) genes, two ribosomal RNA (rRNA) genes, and a control region (CR). Most of the PCGs had typical ATN start codons and terminated with TAA or a single T residue. UUA (Leu2), AUU (Ile), UUU (Phe), AUA (Met) and AAU (Asn) were the five most frequently used codons. All tRNAs were folded into cloverleaf secondary structure, except for trnS1, which lacked the DHU arm. Phylogenetic analyses within the superfamily Zygaenoidea were performed based on multiple datasets from mitochondrial genes. The results showed that the families Phaudidae, Limacodidae and Zygaenidae were respectively recovered as monophyly; C. houshuaii was clustered in a clade with nettle type larvae in Limacodidae.

Phylogenomics resolves the higher-level phylogeny of herbivorous eriophyoid mites (Acariformes: Eriophyoidea)

BACKGROUND

Eriophyoid mites (Eriophyoidea) are among the largest groups in the Acariformes; they are strictly phytophagous. The higher-level phylogeny of eriophyoid mites, however, remains unresolved due to the limited number of available morphological characters-some of them are homoplastic. Nevertheless, the eriophyoid mites sequenced to date showed highly variable mitochondrial (mt) gene orders, which could potentially be useful for resolving the higher-level phylogenetic relationships.

RESULTS

Here, we sequenced and compared the complete mt genomes of 153 eriophyoid mite species, which showed 54 patterns of rearranged mt gene orders relative to that of the hypothetical ancestor of arthropods. The shared derived mt gene clusters support the monophyly of eriophyoid mites (Eriophyoidea) as a whole and the monophylies of six clades within Eriophyoidea. These monophyletic groups and their relationships were largely supported in the phylogenetic trees inferred from mt genome sequences as well. Our molecular dating results showed that Eriophyoidea originated in the Triassic and diversified in the Cretaceous, coinciding with the diversification of angiosperms.

CONCLUSIONS

This study reveals multiple molecular synapomorphies (i.e. shared derived mt gene clusters) at different levels (i.e. family, subfamily or tribe level) from the complete mt genomes of 153 eriophyoid mite species. We demonstrated the use of derived mt gene clusters in unveiling the higher-level phylogeny of eriophyoid mites, and underlines the origin of these mites and their co-diversification with angiosperms.

PPR596 Is Required for nad2 Intron Splicing and Complex I Biogenesis in Arabidopsis

Mitochondria are essential organelles that generate energy via oxidative phosphorylation. Plant mitochondrial genome encodes some of the respiratory complex subunits, and these transcripts require accurate processing, including C-to-U RNA editing and intron splicing. Pentatricopeptide repeats (PPR) proteins are involved in various organellar RNA processing events. PPR596, a P-type PPR protein, was previously identified to function in the C-to-U editing of mitochondrial rps3 transcripts in Arabidopsis. Here, we demonstrate that PPR596 functions in the cis-splicing of nad2 intron 3 in mitochondria. Loss of the PPR596 function affects the editing at rps3eU1344SS, impairs nad2 intron 3 splicing and reduces the mitochondrial complex I's assembly and activity, while inducing alternative oxidase (AOX) gene expression. This defect in nad2 intron splicing provides a plausible explanation for the slow growth of the ppr595 mutants. Although a few P-type PPR proteins are involved in RNA C-to-U editing, our results suggest that the primary function of PPR596 is intron splicing.

Comparative mitogenomic analysis of subterranean and surface amphipods (Crustacea, Amphipoda) with special reference to the family Crangonyctidae

Mitochondrial genomes play important roles in studying genome evolution, phylogenetic analyses, and species identification. Amphipods (Class Malacostraca, Order Amphipoda) are one of the most ecologically diverse crustacean groups occurring in a diverse array of aquatic and terrestrial environments globally, from freshwater streams and lakes to groundwater aquifers and the deep sea, but we have a limited understanding of how habitat influences the molecular evolution of mitochondrial energy metabolism. Subterranean amphipods likely experience different evolutionary pressures on energy management compared to surface-dwelling taxa that generally encounter higher levels of predation and energy resources and live in more variable environments. In this study, we compared the mitogenomes, including the 13 protein-coding genes involved in the oxidative phosphorylation (OXPHOS) pathway, of surface and subterranean amphipods to uncover potentially different molecular signals of energy metabolism between surface and subterranean environments in this diverse crustacean group. We compared base composition, codon usage, gene order rearrangement, conducted comparative mitogenomic and phylogenomic analyses, and examined evolutionary signals of 35 amphipod mitogenomes representing 13 families, with an emphasis on Crangonyctidae. Mitogenome size, AT content, GC-skew, gene order, uncommon start codons, location of putative control region (CR), length of rrnL and intergenic spacers differed between surface and subterranean amphipods. Among crangonyctid amphipods, the spring-dwelling Crangonyx forbesi exhibited a unique gene order, a long nad5 locus, longer rrnL and rrnS loci, and unconventional start codons. Evidence of directional selection was detected in several protein-encoding genes of the OXPHOS pathway in the mitogenomes of surface amphipods, while a signal of purifying selection was more prominent in subterranean species, which is consistent with the hypothesis that the mitogenome of surface-adapted species has evolved in response to a more energy demanding environment compared to subterranean amphipods. Overall, gene order, locations of non-coding regions, and base-substitution rates points to habitat as an important factor influencing the evolution of amphipod mitogenomes.

Circulating mitochondrial long non-coding 7S RNA in primary health care patients with depression/anxiety

BACKGROUND

The significant role of long non-coding 7S RNA in controlling mitochondrial transcription highlights its importance in mitochondrial function. Considering the suggested connection between mitochondrial dysfunction and the onset of mental disorders, this study aimed to explore the potential involvement of 7S RNA in the context of depression/anxiety.

RESULTS

A total of 181 patients in primary health care (age 20-64 years) with depression/anxiety and 59 healthy controls were included in the study. 7S RNA was measured using quantitative real-time PCR in plasma samples collected before (baseline) and after 8 weeks of treatment (mindfulness or cognitive-based behavioral therapy). Upon adjustment for age and sex, the baseline plasma levels of 7S RNA were significantly higher in patients than in healthy controls (p < 0.001). Notably, post-treatment, there was a significant reduction in 7S RNA levels (p = 0.03). These changes in 7S RNA were related to the treatment response, as indicated by HADS-D (Hospital Anxiety and Depression Scale) scores (ß = -0.04, p = 0.04), even after accounting for baseline scores and other cofounders.

CONCLUSION

The findings of this study indicate an association between plasma 7S RNA levels and depression/anxiety, as well as treatment response. While further confirmatory analyses are necessary, plasma 7S RNA holds promise as a potential predictive biomarker for both depression/anxiety and the treatment response within these disorders.

First mitogenomic characterization of Macromotettixoides (Orthoptera, Tetrigidae), with the descriptions of two new species

Classification of species is commonly based on morphological, molecular, and distribution depending on the taxa. Macromotettixoides Zheng, Wei & Jiang, 2005 and Formosatettix Tinkham, 1937 are both wingless types of Tetrigidae with extremely similar morphological characteristics, and in the current taxonomic system they are placed in two different subfamilies, Metrodorinae and Tetriginae, respectively. It is difficult to clearly identify the species of these two genera by morphological characteristics, and molecular data is often needed to assist identification. Here, the complete mitogenomes of two new species were sequenced and assembled, with that of Macromotettixoidesorthomargina. Molecular data of species of Formosatettix were used to test the monophyly of Macromotettixoides and to re-assess the generic characters, and also to test whether Macromotettixoides belongs to the Asian Metrodorinae or Tetriginae. Furthermore, mitochondrial characteristics were analyzed and the phylogeny of the Tetrigidae reconstructed based on mitochondrial protein-coding genes (PCGs). The results indicated that the two new species were clustered with Macromotettixoides rather than Formosatettix, and the anterior margin of the fastigium and pronotum of the two new species usually had the humeral angle different from that of Formosatettix. Therefore, after integrating morphological and molecular data, the two new species were placed in the genus Macromotettixoides, M.maoershanensissp. nov. and M.brachycornasp. nov. Finally, a phylogenetic reconstruction supported Macromotettixoides being assigned to Tetriginae rather than Metrodorinae, in contrast to the previous classification of this genus.

Complete Mitochondrial Genome and Phylogenetic Position of Nurudea zhengii Ren (Insecta, Hemiptera, Aphididae)

Nurudea zhengii Ren was identified by aphid morphological characteristics as well as the gall shape and host plant species, and placed in the tribe Fordini (Hemiptera, Aphididae, Eriosomatinae). Here, its whole genome was firstly sequenced by a genome-skimming method and its mitochondrial genome (mitogenome) was assembled to examine its genetic variation and phylogenetic position. The complete mitogenome of Nurudea zhengii is 15,392 bp in length, and consists of 13 protein-coding genes, 22 tRNAs, two rRNAs and one D-loop region. The gene order follows the mitogenomes of the other Rhus gall aphids, and similarly has an AT bias with the content of 83.9%. The majority strand is A-skewed and C-skewed, and shows opposite skewness for G-skewed in the minority strands. The ratios of nonsynonymous to synonymous substitution rates of protein-coding genes are lower than one except for ATP8, which indicated that ATP8 was undergoing positive selection. Phylogenetic analysis among the Rhus gall aphids based on 13 protein-coding genes and two rRNA genes showed that N. zhengii was sister to N. shiraii, and then clustered with N. yanoniella as a group with high support value. The two species, N. shiraii and N. yanoniella, share the same host plant Rhus chinensis, while the host of N. zhengii is R. hypoleuca. However, the phylogenetic relationship indicated that the taxa sharing the same host plant were not absolutely clustered as the closest taxa at least at species level.

Comparative Analysis of the Systematics and Evolution of the Pampus Genus of Fish (Perciformes: Stromateidae) Based on Osteology, Population Genetics and Complete Mitogenomes

Pampus is a widespread species of fish in the western Pacific and Indian Oceans that has significant commercial worth. Its evolutionary history and phylogenetics are still poorly understood, and details on its intraspecific taxonomy are debatable, despite some morphological and molecular research. Here, we analyzed this species using skeletal structure data as well as nuclear (S7 gene) and mitochondrial genetic information (COI, D-loop and mitogenomes). We found that the genetic distance between P. argenteus and P. echinogaster was much smaller than that between other Pampus species, and both maximum likelihood and Bayesian phylogenetic trees yielded almost identical tree topologies. An additional and adjacent M repeat was found in the downstream region of the IQM gene cluster of P. argenteus and P. echinogaster, and the trnL2 gene of P. minor was translocated. The genus Pampus experienced early rapid radiation during the Palaeocene with major lineages diversifying within a relatively narrow timescale. Additionally, three different methods were conducted to distinguish the genus Pampus species, proving that P. argenteus and P. echinogaster are the same species, and P. liuorum is speculated to be a valid species. Overall, our study provides new insights not only into the evolutionary history of Pampus but its intraspecific taxonomy as well.

Comparative Analysis and Phylogenetic Study of Dawkinsia filamentosa and Pethia nigrofasciata Mitochondrial Genomes

Smiliogastrinae are recognized for their high nutritional and ornamental value. In this study, we employed high-throughput sequencing technology to acquire the complete mitochondrial genome sequences of Dawkinsia filamentosa and Pethia nigrofasciata. The gene composition and arrangement order in these species were similar to those of typical vertebrates, comprising 13 protein-coding genes, 22 tRNA genes, 2 rRNA genes, and 1 non-coding region. The mitochondrial genomes of D. filamentosa and P. nigrofasciata measure 16,598 and 16,948 bp, respectively. Both D. filamentosa and P. nigrofasciata exhibit a significant preference for AT bases and an anti-G bias. Notably, the AT and GC skew values of the ND6 gene fluctuated markedly, suggesting that the selection and mutation pressures on this gene may differ from those affecting other genes. Phylogenetic analysis, based on the complete mitochondrial genomes of 23 Cyprinidae fishes, revealed that D. filamentosa is closely related to the sister group comprising Dawkinsia denisonii and Sahyadria chalakkudiensis. Similarly, P. nigrofasciata forms a sister group with Pethia ticto and Pethia stoliczkana.

Climate-driven mitochondrial selection in lacertid lizards

The mitochondrion, which is an intracellular organelle responsible for most of the energy-producing pathways, can have its genome targeted for climate-driven selection. However, climate-driven mitochondrial selection remains a sparsely studied area in reptiles. Here, we reported the complete mitochondrial genome sequence of a lacertid lizard (Takydromus intermedius) and used mitogenomes from 54 species of lacertid lizards to study their phylogenetic relationships and to identify the mitochondrial genes under positive selection by climate. The length of the complete mitochondrial genome sequence of T. intermedius was 17,713 bp, which was within the range of lengths (17,224-18,943) ever reported for Takydromus species. The arrangement of mitochondrial genes in T. intermedius was the same as in other congeneric species. The 54 lacertid species could be divided into three geographically and climatically different clades. We identified three mitochondrial genes (ATP6, ATP8, and ND3) under positive selection by climate, and found that isothermality, temperature seasonality, precipitation of wettest month, and precipitation seasonality were the most important climatic variables contributing to the gene selection.

Mitochondrially mediated RNA interference, a retrograde signaling system affecting nuclear gene expression

Several functional classes of short noncoding RNAs are involved in manifold regulatory processes in eukaryotes, including, among the best characterized, miRNAs. One of the most intriguing regulatory networks in the eukaryotic cell is the mito-nuclear crosstalk: recently, miRNA-like elements of mitochondrial origin, called smithRNAs, were detected in a bivalve species, Ruditapes philippinarum. These RNA molecules originate in the organelle but were shown in vivo to regulate nuclear genes. Since miRNA genes evolve easily de novo with respect to protein-coding genes, in the present work we estimate the probability with which a newly arisen smithRNA finds a suitable target in the nuclear transcriptome. Simulations with transcriptomes of 12 bivalve species suggest that this probability is high and not species specific: one in a hundred million (1 × 10-8) if five mismatches between the smithRNA and the 3' mRNA are allowed, yet many more are allowed in animals. We propose that novel smithRNAs may easily evolve as exaptation of the pre-existing mitochondrial RNAs. In turn, the ability of evolving novel smithRNAs may have played a pivotal role in mito-nuclear interactions during animal evolution, including the intriguing possibility of acting as speciation trigger.

Facilitating taxonomy and phylogenetics: An informative and cost-effective protocol integrating long amplicon PCRs and third-generation sequencing

Phylogenetic inference has become a standard technique in integrative taxonomy and systematics, as well as in biogeography and ecology. DNA barcodes are often used for phylogenetic inference, despite being strongly limited due to their low number of informative sites. Also, because current DNA barcodes are based on a fraction of a single, fast-evolving gene, they are highly unsuitable for resolving deeper phylogenetic relationships due to saturation. In recent years, methods that analyse hundreds and thousands of loci at once have improved the resolution of the Tree of Life, but these methods require resources, experience and molecular laboratories that most taxonomists do not have. This paper introduces a PCR-based protocol that produces long amplicons of both slow- and fast-evolving unlinked mitochondrial and nuclear gene regions, which can be sequenced by the affordable and portable ONT MinION platform with low infrastructure or funding requirements. As a proof of concept, we inferred a phylogeny of a sample of 63 spider species from 20 families using our proposed protocol. The results were overall consistent with the results from approaches based on hundreds and thousands of loci, while requiring just a fraction of the cost and labour of such approaches, making our protocol accessible to taxonomists worldwide.