MitoZ: a toolkit for animal mitochondrial genome assembly, annotation and visualization

Characterization of the complete mitochondrial genomes of four tarantulas (Arachnida: Theraphosidae) with phylogenetic analysis

To better understand the evolution of mitochondrial genomes (mitogenomes) within the family Theraphosidae, we characterized mitogenomes of four tarantulas (Grammostola pulchripes, Phormictopus atrichomatus, Pterinochilus murinus and Pterinopelma sazimai) for the first time. The mitogenomes were all classical circular structures, with lengths ranging from 13,822 bp to 14,011 bp. The constitutive genes and the orientation of the coding strand observed in the four mitogenomes were consistent with those found in other species belonging to the Theraphosidae family. The four mitogenomes were compacted and exhibited a preference for A and T, with the rRNA sequences showing a higher A+T content. Ka/Ks and p-distances analyses showed the ND6 gene had highest evolutionary rate, while the COⅠ gene displayed relatively slower evolution. In contrast to previous phylogenetic studies, our phylogenetic analysis based on mitogenomes provides new phylogenetic relationships among subfamilies. Subfamily Theraphosinae is most closely related to Ornithoctoninae, slightly distant from Harpactirinae, and farthest from Selenocosmiinae. The new data we acquired regarding these mitogenomes will aid in understanding the complex interrelationships among species within the Theraphosidae family.

Identification and characterization of a novel Pavlovskyella (Acari: Argasidae) from Chile, parasite of the culpeo fox (Lycalopex culpaeus)

Soft ticks (Argasidae) of the subgenus Pavlovskyella Pospelova-Shtrom are worldwide distributed parasites of medical importance. However, the systematics of the subgenus are currently under debate because genetic data shows that the group is paraphyletic. Meanwhile, species of Pavlovskyella continue to be discovered. In this study a novel species of the subgenus is described from specimens collected on a fox in central Chile. The larva of this new species differentiates from other Pavlovskyella spp. by having the following combination of characters: subpyriform dorsal plate; 15 pairs of setae, 7 anterolateral, 3 central, and 5 posterolateral, and hypostome with denticles in the distal third. Nymphs and adults of the species lack cheeks, eyes or bulging structures on the flank, but exhibit dorsoventral grooves, and humps on tarsi I, II, and III. Moreover, a patch of glabrous integument appears on the distal portion of coxal folds. A phylogenetic analysis using the mitogenome indicates a monophyletic group composed by Ornithodoros (Pavlovskyella) brasiliensis Aragão, Ornithodoros (Pavlovskyella) furcosus Neumann, Ornithodoros (Pavlovskyella) improvisus Muñoz-Leal & Venzal, and Ornithodoros (Pavlovskyella) rostratus Aragão. Interestingly, a phylogeny using 18S-28S rDNA sequences shows that South American Pavlovskyella spp. are paraphyletic, as depicted in previous studies. Adding species of the subgenus from the Neotropical region to phylogenetic analyses could aid to solve this paraphyly. Furthermore, this is the fifth species of Pavlovskyella described in South America, and the second in Chile.

The first complete mitochondrial genome of Mylabris sibirica (Coleoptera: Meloidae) and its phylogenetic analysis

The complete mitochondrial genome of Mylabris sibirica Fischer von Waldheim, 1823 was sequenced and characterized. The mitogenome is 15,794 bp long with 37 annotated genes, comprising 13 protein-coding genes, 22 transfer RNA genes, and two ribosomal RNA genes. The phylogenetic analysis based on the mitochondrial genome sequences revealed that M. sibirica clustered with M. quadripunctata. This study presents the complete mitochondrial genome of M. sibirica for the first time, which could be beneficial for systematic studies on Mylabrini.

Comparative Genomics of the World's Smallest Mammals Reveals Links to Echolocation, Metabolism, and Body Size Plasticity

Originating 30 million years ago, shrews (Soricidae) have diversified into around 400 species worldwide. Shrews display a wide array of adaptations, with some species having developed distinctive traits such as echolocation, underwater diving, and venomous saliva. Accordingly, these tiny insectivores are ideal to study the genomic mechanisms of evolution and adaptation. We conducted a comparative genomic analysis of four shrew species and 16 other mammals to identify genomic variations unique to shrews. Using two existing shrew genomes and two de novo assemblies for the maritime (Sorex maritimensis) and smoky (Sorex fumeus) shrews, we identified mutations in conserved regions of the genomes, also known as accelerated regions, gene families that underwent significant expansion, and positively selected genes. Our analyses unveiled shrew-specific genomic variants in genes associated with the nervous, metabolic, and auditory systems, which can be linked to unique traits in shrews. Notably, genes suggested to be under convergent evolution in echolocating mammals exhibited accelerated regions in shrews, and pathways linked to putative body size plasticity were detected. These findings provide insight into the evolutionary mechanisms shaping shrew species, shedding light on their adaptation and divergence over time.

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.

The multifaceted role of mitochondria in cardiac function: insights and approaches

Cardiovascular disease (CVD) remains a global economic burden even in the 21st century with 85% of deaths resulting from heart attacks. Despite efforts in reducing the risk factors, and enhancing pharmacotherapeutic strategies, challenges persist in early identification of disease progression and functional recovery of damaged hearts. Targeting mitochondrial dysfunction, a key player in the pathogenesis of CVD has been less successful due to its role in other coexisting diseases. Additionally, it is the only organelle with an agathokakological function that is a remedy and a poison for the cell. In this review, we describe the origins of cardiac mitochondria and the role of heteroplasmy and mitochondrial subpopulations namely the interfibrillar, subsarcolemmal, perinuclear, and intranuclear mitochondria in maintaining cardiac function and in disease-associated remodeling. The cumulative evidence of mitochondrial retrograde communication with the nucleus is addressed, highlighting the need to study the genotype-phenotype relationships of specific organelle functions with CVD by using approaches like genome-wide association study (GWAS). Finally, we discuss the practicality of computational methods combined with single-cell sequencing technologies to address the challenges of genetic screening in the identification of heteroplasmy and contributory genes towards CVD.

A Snakemake Toolkit for the Batch Assembly, Annotation and Phylogenetic Analysis of Mitochondrial Genomes and Ribosomal Genes From Genome Skims of Museum Collections

Low coverage 'genome-skims' are often used to assemble organelle genomes and ribosomal gene sequences for cost-effective phylogenetic and barcoding studies. Natural history collections hold invaluable biological information, yet poor preservation resulting in degraded DNA often hinders polymerase chain reaction-based analyses. However, it is possible to generate libraries and sequence the short fragments typical of degraded DNA to generate genome-skims from museum collections. Here we introduce a snakemake toolkit comprised of three pipelines skim2mito, skim2rrna and gene2phylo, designed to unlock the genomic potential of historical museum specimens using genome skimming. Specifically, skim2mito and skim2rrna perform the batch assembly, annotation and phylogenetic analysis of mitochondrial genomes and nuclear ribosomal genes, respectively, from low-coverage genome skims. The third pipeline gene2phylo takes a set of gene alignments and performs phylogenetic analysis of individual genes, partitioned analysis of concatenated alignments and a phylogenetic analysis based on gene trees. We benchmark our pipelines with simulated data, followed by testing with a novel genome skimming dataset from both recent and historical solariellid gastropod samples. We show that the toolkit can recover mitochondrial and ribosomal genes from poorly preserved museum specimens of the gastropod family Solariellidae, and the phylogenetic analysis is consistent with our current understanding of taxonomic relationships. The generation of bioinformatic pipelines that facilitate processing large quantities of sequence data from the vast repository of specimens held in natural history museum collections will greatly aid species discovery and exploration of biodiversity over time, ultimately aiding conservation efforts in the face of a changing planet.

The complete mitochondrial genome of Parotis chlorochroalis (Hampson, 1912) (Lepidoptera: Crambidae)

The complete mitochondrial genome of the Parotis chlorochroalis was sequenced, revaeling a length of 15239 bp with 37 genes and an A + T-rich region. All c13 PCGs begin with typical ATN codons, except COI gene, which starts with CGA. Eleven genes terminate with TAA, two with T-. All 22 tRNA genes exhibit typical cloverleaf structure except for trnS1 P. chlorochroalis has two relatively conserved intergenic regions and two relatively conserved overlapping regions. Phylogenetic analysis support P. chlorochroalis belongs to subfamily Spilomelinae, the topologies of Crambidae are highly congruent with previous studies. This newly sequences mitochondrial genome provides valuable resources for taxonomic inference and evolutionary studies of genus Parotis.

Genomic and transcriptomic perspectives on the origin and evolution of NUMTs in Orthoptera

Nuclear mitochondrial pseudogenes (NUMTs) result from the transfer of mitochondrial DNA (mtDNA) to the nuclear genome. NUMTs, as "frozen" snapshots of mitochondria, can provide insights into diversification patterns. In this study, we analyzed the origins and insertion frequency of NUMTs using genome assembly data from ten species in Orthoptera. We found divergences between NUMTs and contemporary mtDNA in Orthoptera ranging from 0 % to 23.78 %. The results showed that the number of NUMT insertions was significantly positively correlated with the content of transposable elements in the genome. We found that 39.09 %-68.65 % of the NUMTs flanking regions (2,000 bp) contained retrotransposons, and more NUMTs originated from mitochondrial rDNA regions. Based on the analysis of the mitochondrial transcriptome, we found a potential mechanism of NUMT integration: mitochondrial transcripts are reverse transcribed into double-stranded DNA and then integrated into the genome. The probability of this mechanism occurring accounts for 0.30 %-1.02 % of total mitochondrial nuclear transfer events. Finally, based on the phylogenetic tree constructed using NUMTs and contemporary mtDNA, we provide insights into ancient evolutionary events such as species-specific "autaponumts" and "synaponumts" shared among different species, as well as post-integration duplication events.

The complete mitochondrial genome and phylogenetic analysis for Rhabdophis chiwen (squamata: colubridae)

Rhabdophis chiwen is currently found so far in Sichuan Province, China, where it predominantly feeds on earthworms and firefly larvae. In this study, we sequenced and analyzed the mitochondrial genome of R. chiwen, which measured 17,646 bp in length and encompassed 37 genes along with two control regions. The base composition revealed percentages of 33.20% A, 25.94% T, 13.27% G, and 27.59% C. Phylogenetic analyses indicate that R. chiwen belongs to the family Colubridae and forms a sister branch with R. tigrinus. This study successfully obtained the first complete mitochondrial genome of R. chiwen, offering crucial genetic data for its evolutionary history and conservation.

Complementing aculiferan mitogenomics: comparative characterization of mitochondrial genomes of Solenogastres (Mollusca, Aplacophora)

BACKGROUND

With the advances in high-throughput sequencing and bioinformatic pipelines, mitochondrial genomes have become increasingly popular for phylogenetic analyses across different clades of invertebrates. Despite the vast rise in available mitogenomic datasets of molluscs, one class of aplacophoran molluscs - Solenogastres (or Neomeniomorpha) - is still neglected.

RESULTS

Here, we present six new mitochondrial genomes from five families of Solenogastres (Amphimeniidae, Gymnomeniidae, Proneomeniidae, Pruvotinidae, Simrothiellidae), including the first complete mitogenomes, thereby now representing three of the four traditional orders. Solenogaster mitogenomes are variable in size (ranging from approximately 15,000 bp to over 17,000 bp). The gene order of the 13 protein coding genes and two rRNA genes is conserved in three blocks, but considerable variation occurs in the order of the 22 tRNA genes. Based on phylogenetic analyses and reconstruction of ancestral mitochondrial genomes of Aculifera, the position of (1) trnD gene between atp8 and atp6, (2) trnT and P genes between atp6 and nad5, and (3) trnL1 gene between G and E, resulting in a 'MCYWQGL1E'-block of tRNA genes, are all three considered synapomorphies for Solenogastres. The tRNA gene block 'KARNI' present in Polyplacophora and several conchiferan taxa is dissolved in Solenogastres.

CONCLUSION

Our study shows that mitogenomes are suitable to resolve the phylogenetic relationships among Aculifera and within Solenogastres, thus presenting a cost and time efficient compromise to approach evolutionary history in these clades.

Novel gene arrangement in the mitochondrial genome of Aspersentis megarhynchus (Acanthocephala, Echinorhynchida, Heteracanthocephalidae), and its phylogenetic implications

The Heteracanthocephalidae Petrochenko, 1956 is a rare family of acanthocephalans mainly parasitic in fishes. The pattern of mitogenomic evolution of the Heteracanthocephalidae is still unknown, and the phylogenetic relationships of the Heteracanthocephalidae with the other 14 families within the order Echinorhynchida remain unclear. In the present study, the complete mitochondrial genome of Aspersentis megarhynchus (von Linstow, 1892) Golvan, 1960 was sequenced and annotated for the first time, which represents the first mitogenomic data for the genus Aspersentis and also for the family Heteracanthocephalidae. The mitogenome of A. megarhynchus has 14,661 bp and includes 36 genes, containing 12 protein-coding genes (PCGs) (missing atp8), 22 tRNA genes, and 2 ribosomal RNAs (rrnS and rrnL), plus two non-coding regions. Comparative mitochondrial genomic analysis revealed that the presence of translocations of several tRNA genes (trnV, trnE, and trnT) and the gene arrangement in the mitogenome of A. megarhynchus represents a new type in Acanthocephala. Moreover, the mitogenomic phylogenetic results based on concatenated amino acid sequences of 12 protein-coding genes strongly supported the validity of the Heteracanthocephalidae and suggested close affinity between the Heteracanthocephalidae and Echinorhynchidae in the order Echinorhynchida.

Complete mitochondrial genome of Saldoida armata Horváth, 1911 (Heteroptera: Saldidae) and phylogenetic analysis

The complete mitochondrial genome of Saldoida armata (Heteroptera: Saldidae) is 16,049 bp in length, comprising 13 protein-coding genes (PCGs), 22 transfer RNA genes (tRNAs), 2 ribosomal RNA genes (rRNAs) and a control region. All the PCGs are initially encoded by ATN, TTG or GTG, and terminated coding with TAA or a single T. With the exception of trnS(AGN), all tRNAs exhibit a typical cloverleaf secondary structure. Phylogenetic analysis reveals the sister relationship of S. armata with other Saldidae members. The complete mitogenome of S. armata will provide useful genetic information for species identification, phylogenetic analysis and conservation of this species.

De Novo Assembly and Annotation of the Siganus fuscescens (Houttuyn, 1782) Genome: Marking a Pioneering Advance for the Siganidae Family

This study presents the first draft genome of Siganus fuscescens, and thereby establishes the first whole-genome sequence for a species in the Siganidae family. Leveraging both long and short read sequencing technologies, i.e., Oxford Nanopore and Illumina sequencing, we successfully assembled a mitogenome spanning 16.494 Kb and a first haploid genome encompassing 498 Mb. The assembled genome accounted for a 99.6% of the estimated genome size and was organized into 164 contigs with an N50 of 7.2 Mb. This genome assembly showed a GC content of 42.9% and a high Benchmarking Universal Single-Copy Orthologue (BUSCO) completeness score of 99.5% using actinopterygii_odb10 lineage, thereby meeting stringent quality standards. In addition to its structural aspects, our study also examined the functional genomics of this species, including the intricate capacity to biosynthesize long-chain polyunsaturated fatty acids (LC-PUFAs) and secrete venom. Notably, our analyses revealed various repeats elements, which collectively constituted 17.43% of the genome. Moreover, annotation of 28,351 genes uncovered both shared genetic signatures and those that are unique to S. fuscescens. Our assembled genome also displayed a moderate prevalence of gene duplication compared to other fish species, which suggests that this species has a distinctive evolutionary trajectory and potentially unique functional constraints. Taken altogether, this genomic resource establishes a robust foundation for future research on the biology, evolution, and the aquaculture potential of S. fuscescens.

Revealing the complete mtDNA genome sequence of Cemani chicken (Gallus gallus) by using Nanopore sequencing analysis

OBJECTIVE

This study aimed to identify, discover and explore the characteristics of the mtDNA genomes of Cemani chicken (Gallus gallus).

METHODS

This study used gDNA of Cemani chicken isolated from liver tissue. mtDNA sequencing was performed using WGS mtDNA analysis with nanopore technology by Oxford Nanopore Technologies GridION. Bioinformatics and data analysis were then performed.

RESULTS

This study showed that the length of the mtDNA genome is 16,789 bp, consisting of two ribosomal RNA (12S rRNA, 16S rRNA), 22 transfer RNA genes (trnR, trnG, trnK, trnD, trnS, trnY, trnC, trnN, trnA, trnW, trnM, trnQ, trnl, trnL, trnV, trnF, trnP, trnT, trnE, trnL, trnS, trnH), 13 protein-coding genes (PCGs) (ND4l, ND3, COX3, ATP6, ATP8, COX2, COX1, ND2, ND1, CYTB, ND6, ND5, ND4), and a noncoding control region (Dloop). Furthermore, analysis showed there were polymorphic sites and amino acid alterations when mtDNA Cemani chicken was aligned with references from GenBank.

CONCLUSION

Site (988T>*) in Dloop genes and (328A>G) in ND3 genes which alter glycine to stop codon, were specific markers found only in Cemani chicken.

The Complete Mitochondrial Genomes of Five Nivanini Species (Hemiptera: Cicadellidae: Evacanthinae) With Phylogenetic Analysis

To delve deeper into the phylogenetic relationships within Cicadellidae and the taxonomic arrangement of Evacanthinae, our study focuses on the mitochondrial genome sequencing of five Nivanini species: Extensus latus, Concavepiana hamulusa, Sophonia nigrilineata, Sophonia microstaina, and Sophonia fuscomarginata. The results showed that the length of the five mitochondrial genomes ranged from 15,610 to 16,032 bp and included 37 typical genes. The A + T content of Nivanini ranged from 72.5% to 78.7%, which is consistent with that of other sequenced Evacanthinae species. All transfer RNAs (tRNAs) exhibit the typical cloverleaf secondary structure, except for trnS (AGY), which lacks the dihydrouracil (DHU) arm. With regard to protein-coding genes, all started with ATN codons, except for atp8, and most of them use TAA or TAG as termination codons. Using the Bayesian inference and maximum likelihood methods, a phylogenetic tree based on all 37 genes was constructed, with a total of 57 Cicadellidae species and two outgroups included as research objects. The analyses confirmed the monophyletic nature of each subfamily, highlighting Deltocephalinae as the oldest, distinctively parallel to the others.

Complete mitochondrial genome and phylogenetic analysis of Chloris chloris (Passeriformes: Fringillidae)

In this study, we employed high-throughput metagenomic data to assemble the mitochondrial genome (mitogenome) of the European greenfinch (Chloris chloris; Linnaeus 1758). The circular mitogenome was 16,813 base pairs (bp) in length, containing 13 protein-coding genes (PCGs), 22 tRNA genes, 2 rRNA genes, and 1 control region. The base composition of the mitogenome is 30.6% A, 30.7% C, 14.2% G, and 24.5% T, resulting in a GC content of 44.9%. Phylogenetic analysis, utilizing the concatenation of the 13 mitochondrial PCGs from 32 related species of the order Passeriformes, indicated a closer relationship between C. chloris and C. sinica. Moreover, the genus Chloris was closely related to the genera of Serinus, Crithagra, Carduelis, and Acanthis. This mitogenomic data of C. chloris not only be helpful for species identification but also facilitates our understanding of the evolutionary relationship among different species in genus Chloris, which experienced rapid radiation evolution.

Chromosome-level genome assembly and annotation of Clanis bilineata tsingtauica Mell (Lepidoptera: Sphingidae)

The soybean hawkmoth Clanis bilineata tsingtauica Mell (Lepidoptera, Sphingidae; CBT), as one of the main leaf-chewing pests of soybeans, has gained popularity as an edible insect in China recently due to its high nutritional value. However, high-quality genome of CBT remains unclear, which greatly limits further research. In the present study, we assembled a high-quality chromosome-level genome of CBT using PacBio HiFi reads and Hi-C technologies for the first time. The size of the assembled genome is 477.45 Mb with a contig N50 length of 17.43 Mb. After Hi-C scaffolding, the contigs were anchored to 29 chromosomes with a mapping rate of 99.61%. Benchmarking Universal Single-Copy Orthologues (BUSCO) completeness value is 99.49%. The genome contains 252.16 Mb of repeat elements and 14,214 protein-coding genes. In addition, chromosomal synteny analysis showed that the genome of CBT has a strong synteny with that of Manduca sexta. In conclusion, this high-quality genome provides an important resource for future studies of CBT and contributes to the development of integrated pest management strategies.

Positive Selection of Mitochondrial cytochrome b Gene in the Marine Bivalve Keenocardium buelowi (Bivalvia, Cardiidae)

The mitochondrial genome provides valuable data for phylogenetic analysis and evolutionary research. In this study, we sequenced, assembled, and annotated the mitochondrial genome of Keenocardium buelowi using the Illumina platform. The genome spanned 16,967 bp and included 13 protein-coding genes (PCGs), two ribosomal RNAs, and 22 transfer RNAs. All PCGs utilized standard ATN start codons and TAN stop codons. The phylogenetic tree based on maximum likelihood and Bayesian inference analyses revealed Clinocardiinae as the sister group to Trachycardiinae, with the estimated divergence time being 44.5 million years ago (MYA) between K. buelowi and Vasticardium flavum. Notably, the cytochrome b gene (cob) exhibited a positive selection signal. Our findings provide valuable insights into the evolutionary history and molecular phylogeny of K. buelowi.

Four complete mitochondrial genomes of the subgenus Pterelachisus (Diptera, Tipulidae, Tipula) and implications for the higher phylogeny of the family Tipulidae

The complete mitochondrial genomes of Tipula (Pterelachisus) cinereocincta mesacantha Alexander, 1934, T. (P.) legalis Alexander, 1933, T. (P.) varipennis Meigen, 1818, and T. (P.) yasumatsuana Alexander, 1954 are reported, three of them being sequenced for the first time. The mitochondrial genome lengths of the four species are 15,907 bp, 15,625 bp, 15,772 bp, and 15,735 bp, respectively. All genomes exhibit a high AT base composition, with A + T content of 76.7%, 75.0%, 77.8%, and 75.4%, respectively. The newly reported mitogenomes herein show a general similarity in overall structure, gene order, base composition, and nucleotide content to those of the previously studied species within the family Tipulidae. Phylogenetic analyses were conducted to investigate the relationships within Tipulidae, using both Maximum Likelihood and Bayesian Inference approaches. The results show that the four target species of the subgenus T. (Pterelachisus) basically form a monophyletic group within Tipulidae, clustering with species of the Tipula subgenera T. (Lunatipula), T. (Vestiplex), and T. (Formotipula); however, the genus Tipula is not monophyletic. Moreover, neither the tipulid subfamily Tipulinae nor the family Limoniidae is supported to be a monophyletic group. The monophyly of the family Tipulidae, and the sister relationship between Tipulidae and Cylindrotomidae are reconfirmed. These research findings could contribute to deep insights into the systematic and evolutionary patterns of crane flies.

Molecular Phylogenetics and Mitochondrial Genomic Evolution in the Endemic Genus Pielomastax (Orthoptera: Eumastacoidea) in China

BACKGROUND/OBJECTIVES

The genus Pielomastax Chang (Orthoptera: Eumastacoidea, 1937) is endemic to China, which is mainly distributed in low- and medium-altitude areas in central and eastern China. However, there are relatively few molecular data studies on the genus Pielomastax.

METHODS

In this study, three species of the genus Pielomastax were collected from Hubei and Henan, China, namely Pielomastax sp., Pielomastax shennongjiaensis Wang (1995) and Pielomastax tenuicerca Hsia and Liu (1989). Both Pielomastax sp. and Pielomastax shennongjiaensis were collected from the Shennongjia area of Hubei, but they exhibit some differences in morphological characteristics.

RESULTS

We obtained the mitochondrial genome structures of the three species, which were similar to those of the published mitochondrial genome structures of species within Eumastacoidea with 37 typical mitochondrial genes, including 13 PCGs, 22 tRNAs, and 2 ribosomal RNAs. The results of the maximum likelihood (ML) tree and the Bayesian inference (BI) tree showed that the families Eumastacidae, Chorotypidae and Episactinae in Eumastacoidea are a monophyletic group, and Thericleinae and Episactinae are sister clades. The time-calibrated phylogeny results indicated that the divergence time between Thericleinae and Episactinae was 95.58 Ma (56.71-128.02 Ma).

CONCLUSIONS

These phylogenetic tree results indicate that Pielomastax sp. and Pielomastax shennongjiaensis are the same species. And the time-calibrated phylogeny tree and the species distribution map of the genus Pielomastax indicate that the species of the genus Pielomastax spread from eastern to central China and diversified. These studies fill the gap in molecular data for the genus Pielomastax and the taxonomic status of Episactidae.

Mitogenome-Based Phylogeny with Divergence Time Estimates Revealed the Presence of Cryptic Species within Heptageniidae (Insecta, Ephemeroptera)

Heptageniidae are known for their flat heads and bodies and are divided into three subfamilies. Despite the extensive diversity within this group and considerable efforts made to understand their evolutionary history, the internal classifications and origin time of Heptageniidae remains controversial. In this study, we newly sequenced 17 complete mitogenomes of Heptageniidae to reconstruct their phylogenetic positions within this family. Because of the ambiguous time of origin, our study also estimated the divergence time within Heptageniidae based on five fossil calibration points. The results of BI and ML trees all highly supported the monophyly of Heptageniidae and three subfamilies. The phylogenetic relationship of Rhithrogeninae + (Ecdyonurinae + Heptageniinae) was also recovered. The divergence time showed that Heptageniidae originated from 164.38 Mya (95% HPD, 150.23-181.53 Mya) in the mid-Jurassic, and Rhithrogeninae originated from 95.54 Mya (95% HPD, 73.86-120.19 Mya) in the mid-Cretaceous. Ecdyonurinae and Heptageniinae began to diverge at 90.08 Mya (95% HPD, 68.81-113.16 Mya) in the middle Cretaceous. After morphological identification, analysis of the mitogenome's composition, genetic distance calculation, phylogenetic analysis, and divergence time calculation, we suggest that two different populations of Epeorus montanus collected from Aksu, Xinjiang Uygur Autonomous Region (40°16' N, 80°26' E) and Xinyuan, Xinjiang Uygur Autonomous Region (43°20' N, 83°43' E) in China are cryptic species of E. montanus, but further detailed information on their morphological characteristics is needed to fully identify them.

Mitochondrial genome analysis across different populations reveals the intraspecific variation and phylogeography of the Caucasian soft tick relapsing fever vector, Ornithodoros (Pavlovskyella) verrucosus (Ixodida: Argasidae)

Territories in southern parts of Eastern Europe and in the Caucasus are endemic for tick-borne relapsing fever (TBRF), caused by Borrelia caucasica. This spirochete is transmitted exclusively by the bites of Ornithodoros verrucosus; however, the distribution and genetic diversity of the tick vector have not been explored. To address this, we performed a phylogeographic study of O. verrucosus specimens collected across a large geographic distribution. We sequenced and analyzed complete mitochondrial genomes of 54 individual O. verrucosus ticks representing 23 geographically diverse populations from Ukraine, Georgia, and Azerbaijan. We detected 47 unique haplotypes, with every collection site exhibiting distinct polymorphisms. This, along with other population genetic indices, suggests little evidence of gene flow between populations. The Bayesian coalescent analysis revealed the presence of four lineages that diverged in the Middle Pleistocene (770-126 kya). Two lineages were widespread and present in all study regions, while the other two were restricted to the southern foothills of the Lesser Caucasus mountain range. The sympatry of these ancient lineages suggests that isolation by environment, in addition to geographic distance, may play a role in the intraspecific divergence of tick populations. Using a phylogeographic approach, we provide a snapshot of genetic diversity in O. verrucosus and discuss the evolutionary history of the tick vector.

Identification of non-model mammal species using the MinION DNA sequencer from Oxford Nanopore

Background

The Neotropics harbors the largest species richness of the planet; however, even in well-studied groups, there are potentially hundreds of species that lack a formal description, and likewise, many already described taxa are difficult to identify using morphology. Specifically in small mammals, complex morphological diagnoses have been facilitated by the use of molecular data, particularly from mitochondrial sequences, to obtain accurate species identifications. Obtaining mitochondrial markers implies the use of PCR and specific primers, which are largely absent for non-model organisms. Oxford Nanopore Technologies (ONT) is a new alternative for sequencing the entire mitochondrial genome without the need for specific primers. Only a limited number of studies have employed exclusively ONT long-reads to assemble mitochondrial genomes, and few studies have yet evaluated the usefulness of such reads in multiple non-model organisms.

Methods

We implemented fieldwork to collect small mammals, including rodents, bats, and marsupials, in five localities in the northern extreme of the Cordillera Central of Colombia. DNA samples were sequenced using the MinION device and Flongle flow cells. Shotgun-sequenced data was used to reconstruct the mitochondrial genome of all the samples. In parallel, using a customized computational pipeline, species-level identifications were obtained based on sequencing raw reads (Whole Genome Sequencing). ONT-based identifications were corroborated using traditional morphological characters and phylogenetic analyses.

Results

A total of 24 individuals from 18 species were collected, morphologically identified, and deposited in the biological collection of Universidad EAFIT. Our different computational pipelines were able to reconstruct mitochondrial genomes from exclusively ONT reads. We obtained three new mitochondrial genomes and eight new molecular mitochondrial sequences for six species. Our species identification pipeline was able to obtain accurate species identifications for up to 75% of the individuals in as little as 5 s. Finally, our phylogenetic analyses corroborated the identifications from our automated species identification pipeline and revealed important contributions to the knowledge of the diversity of Neotropical small mammals.

Discussion

This study was able to evaluate different pipelines to reconstruct mitochondrial genomes from non-model organisms, using exclusively ONT reads, benchmarking these protocols on a multi-species dataset. The proposed methodology can be applied by non-expert taxonomists and has the potential to be implemented in real-time, without the need to euthanize the organisms and under field conditions. Therefore, it stands as a relevant tool to help increase the available data for non-model organisms, and the rate at which researchers can characterize life specially in highly biodiverse places as the Neotropics.

Characterization and phylogenetic analysis of the complete mitochondrial genome of Cotylorhiza tuberculata assembled using next-generation sequencing

In this study, the complete mitochondrial genome (mitogenome) of Cotylorhiza tuberculata (Scyphozoa; Rhizostomeae; Cepheidae) was assembled by the next-generation sequencing data. The complete mitogenome spanned 16,590 bp and contained 14 protein-coding genes, two transfer RNA genes, and two ribosomal RNA genes. Total AT% content was 67.7%, comprising A 30.22%, C 16.16%, G 17.05%, and T 36.56%. The gene arrangement exhibited consistency with the known mitogenomes of other jellyfish species. Furthermore, the phylogenetic relationship of C. tuberculata was investigated based on analysis of the 13 common protein-coding genes. Results indicated a close relationship between C. tuberculata and both Cassiopea xamachana and Cassiopea andromeda. These findings provide a valuable reference for advancing understanding of the phylogenetic relationships, taxonomic classification, and phylogeography of jellyfish 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.

Genome variations in sea cucumbers: Insights from genome survey sequencing and comparative analysis of mitochondrial genomes

Sea cucumbers, marine benthic invertebrates, play crucial roles in maintaining the stability of marine ecosystems and hold key evolutionary positions. However, information regarding their genomes remains limited. Here, we conducted genome survey analyses on seven species from four orders. Results indicated that Colochirus anceps, Colochirus quadrangularis, and Pseudocolochirus violaceus within the order Dendrochirotida have significantly larger genomes (2238-3754 Mbp) compared to conventional sea cucumber genomes, accompanied by a very high proportion of repeat sequences (69.39-72.52 %). While Holothuria edulis and Holothuria atra exhibited similar genome sizes comparable to those of other species within the order Holothuriida, heterozygosity and repeat content varied among all the six species in this order. The representative species Apostichopus californicus of the order Synallactida possesses the smallest genome size (573.45Mbp) within its order, but its heterozygosity (2.24 %) is significantly higher than that of other species. The representative species Synapta maculata of the order Apodida exhibited a normal genome size (900.97 Mbp), lower proportion of repeat sequences (42.19 %), and lower heterozygosity (0.84 %), making it the species with the least challenges for genome sequencing and assembly in the future among all surveyed species. Subsequently, we compiled genomic information from a total of 19 sea cucumber genomes, both newly sequenced and previously reported, revealing a significant linear relationship (P = 0.0001) between genome size and the proportion of repeat sequences in sea cucumbers. Additionally, phylogenetic and comparative analysis of mitochondrial genomes among them indicated extensive rearrangements within the order Apodida, leading to significant discrepancies between mitochondrial and nuclear genome phylogenies.

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.

Population genomics reveals mito-nuclear discordance and admixed populations in southern giraffe

Studying wildlife taxonomic diversity and identifying distinct populations has traditionally been largely based on morphology and geographic origin. More recently, this method has been supplemented by genetic data from the mitochondrial genome. However, this is limited as only maternally inherited and may not reflect the true nature of a population's genetics. Within the giraffe (Giraffa spp.), subspecies and unique populations were successfully characterized using both mitochondrial and genomic DNA studies, which led to new insights and, in some cases, unexpected results that required further verification. Here, we sequenced the genomes of 85 southern giraffe (G. giraffa) individuals from ten populations across southern Africa for a detailed investigation into the genetic diversity and history of its two subspecies, the Angolan (G. g. angolensis) and the South African (G. g. giraffa) giraffe. While the overall genotypes show low levels of runs of homozygosity compared to other mammals, the degree of heterozygosity is limited despite the large population size of South African giraffe. The nuclear genotype is largely congruent with the mitochondrial genotype. However, we have identified that the distribution of the Angolan giraffe is not as far east as indicated in an earlier mitochondrial DNA study. Botswana's Central Kalahari Game Reserve giraffe are unique, with a clear admixture of Angolan and South African giraffe populations. However, the enigmatic desert-dwelling giraffe of northwest Namibia is locally distinct from other Angolan giraffe yet exhibits intra-subspecies signs of admixture resulting from a recent introduction of individuals from Namibia's Etosha National Park. Whole genome sequencing is an invaluable and nearly indispensable tool for wildlife management to uncover genetic diversity that is undetectable through mitogenomic, geographical, and morphological means.

Mitogenomic profiling and gut microbial analysis of the newly identified polystyrene-consuming lesser mealworm in Kenya

Plastic waste has recently become a major global environmental concern and one of the biggest challenges has been seeking for alternative management options. Several studies have revealed the potential of several coleopteran species to degrade plastics, and this is the first research paper on plastic-degradation potential by lesser mealworms from Africa. This study evaluated the whole mitogenomic profile of the lesser mealworm to further identify the insect. The ability of the mealworm to consume Polystyrene (PS) was also evaluated alongside its associated gut microbiota diversity. Our results showed a complete circular mitochondrial genome which clustered closely to the Alphitobius genus but also suggested that our insect might be a new subspecies which require further identification. During the PS feeding trials, overall survival rates of the larvae decreased when fed a sole PS diet while PS intake was observed to increase over a 30-day period. The predominant bacteria observed in larvae fed PS diets were Kluyvera, Lactococcus, Klebsiella, Enterobacter, and Enterococcus, while Stenotrophomonas dominated the control diet. These findings demonstrated that the newly identified lesser mealworm can survive on a PS diet and has a consortium of important bacteria strongly associated with PS degradation. This work provides a better understanding of bioremediation applications, paving the way for further research into the metabolic pathways of plastic-degrading microbes and bringing hope to solving plastic waste pollution while providing high-value insect protein towards a circular economy.

The complete mitochondrial genome of Chlaenius bimaculatus Dejean, 1826 (Coleoptera: Carabidea) and its phylogenetic analyses

Chlaenius bimaculatus Dejean, 1826 (Coleoptera: Carabidea) is a predator of several lepidopteran pests, including Spodoptera frugiperda, S. litura and Helcystogramma triannulella. However, there has been little research into using C. bimaculatus to control crop pests. In this study, we sequenced the complete mitochondrial genome of C. bimaculatus. The results showed that the entire mitochondrial genome was 16,419 bp and contained 24% GC. 13 protein-coding, 22 transfer RNA, and two ribosomal RNA genes were identified. C. bimaculatus shares the same genetic arrangement and composition as other Coleoptera insects. In addition, phylogenetic analysis revealed that C. bimaculatus is closely related to Diplocheila zeelandica.

The complete mitochondrial genome of aglaeactis castelnaudii (bourcier & mulsant, 1848) (apodiformes: trochilidae: aglaeactis) and phylogenetic analysis

In this study, we employed high-throughput sequencing data to assemble the mitochondrial genome (mitogenome) of the White-tufted Sunbeam (Aglaeactis castelnaudii). The total length of the mitogenome was found to be 16,872 base pairs (bp), containing 13 protein-coding genes (PCGs), 22 transfer RNA genes, 2 ribosomal RNA genes, and 1 control region. The nucleotide composition was as follows: A 30.6%, T 24.0%, C 31.2%, and G 14.2%, resulting in a GC content of 45.4%. Phylogenetic analysis, utilizing the concatenation of the 13 mitochondrial PCGs, indicated a closer evolutionary relationship between the genus Aglaeactis and the genus Coeligena compared to other genera within the family Trochilidae investigated in this study. The mitogenome of A. castelnaudii not only contributes to species identification but also provides valuable insights for phylogenetic and conservation genetic analyses of A. castelnaudii.

The complete mitochondrial genome of Triplophysa grahami Regan 1906 (Cypriniformes: Nemacheilidae) and phylogenetic analysis

Triplophysa grahami Regan 1906 is a member of the family Nemacheilidae, Cypriniformes, and native loach in Yunnan. In this study, the complete mitochondrial genome (mitogenome) of T. grahami Regan 1906 was firstly reported and analyzed. The mitogenome of T. grahami Regan 1906 is 16,566 bp in length, including 13 protein-coding genes (PCGs), 22 transfer RNAs (tRNAs), two ribosomal RNAs (rRNAs), and one control region (D-loop). The arrangement and orientation of protein coding genes and RNAs in T. grahami Regan 1906 are identical to other species of Nemacheilidae. The base composition of T. grahami Regan 1906 mitogenome was 29.25% A, 28.55% T, 25.03% C, and 17.17% G. The phylogenetic analysis based on the mitogenome showed that T. grahami Regan 1906 belongs to the clade of genus Triplophysa and the monophyly of Triplophysa is identified. This study contributed valuable genetic data for T. grahami Regan 1906 and explored the phylogenetic relationships in Nemacheilidae.

Repetitive element expansions contribute to genome size gigantism in Pamphagidae: A comparative study (Orthoptera, Acridoidea)

Pamphagidae is a family of Acridoidea that inhabits the desert steppes of Eurasia and Africa. This study employed flow cytometry to estimate the genome size of eight species in the Pamphagidae. The results indicate that the genome size of the eight species ranged from 13.88 pg to 14.66 pg, with an average of 14.26 pg. This is the largest average genome size recorded for the Orthoptera families, as well as for the entire Insecta. Furthermore, the study explored the role of repetitive sequences in the genome, including their evolutionary dynamics and activity, using low-coverage next-generation sequencing data. The genome is composed of 14 different types of repetitive sequences, which collectively make up between 59.9% and 68.17% of the total genome. The Pamphagidae family displays high levels of transposable element (TE) activity, with the number of TEs increasing and accumulating since the family's emergence. The study found that the types of repetitive sequences contributing to the TE outburst events are similar across species. Additionally, the study identified unique repetitive elements for each species. The differences in repetitive sequences among the eight Pamphagidae species correspond to their phylogenetic relationships. The study sheds new light on genome gigantism in the Pamphagidae and provides insight into the correlation between genome size and repetitive sequences within the family.

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.

Integrative taxonomy, phylogenetics and historical biogeography of subgenus Aeschyntelus Stål, 1872 (Hemiptera: Heteroptera: Rhopalidae)

The subgenus Aeschyntelus includes six species that show variations in body color and shape, thus making it difficult to identify them based on morphological identification alone. To date, no genetic study has evaluated species within this genus. Herein, we collected 171 individuals from 90 localities of Rhopalus and employed an integrative taxonomic approach that incorporated morphological data, mitochondrial genomic data (COI, whole mitochondrial data) and nuclear genomic data (18S + 28S rRNAs, nuclear genome-wide SNPs) to delineate species boundaries. Our analyses confirmed the status of nine described species of Rhopalus and proposed the recognition of one new species known as Rhopalus qinlinganus sp. nov., which is classified within the subgenus Aeschyntelus. Discrepancies arising from nuclear and mitochondrial data suggest the presence of mito-nuclear discordance. Specifically, mitochondrial data indicated admixture within Clade A, comprising R. kerzhneri and R. latus, whereas genome-wide SNPs unambiguously identified two separate species, aligning with morphological classification. Conversely, mitochondrial data clearly distinguished Clade B- consisting of R. sapporensis into two lineages, whereas genome-wide SNPs unequivocally identified a single species. Our study also provides insights into the evolutionary history of Aeschyntelus, thus indicating that it likely originated in East Asia during the middle Miocene. The development of Aeschyntelus biodiversity in the southwestern mountains of China occurred via an uplift-driven diversification process. Our findings highlight the necessity of integrating both morphological and multiple molecular datasets for precise species identification, particularly when delineating closely related species. Additionally, it reveals the important role of mountain orogenesis on speciation within the southwestern mountains of China.

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.

Two complete mitochondrial genomes of Boulenophrys (Anura: Megophryidae: Megophryinae): characteristics and phylogenetic implications

The Chinese horned toads, Boulenophrys boettgeri (Boulenger, 1899) and Boulenophrys kuatunensis (Pope, 1929), are two captivating species within the family Megophryidae, which inhabit the mountainous streams in the Eastern of China. In this study, two new complete mitochondrial genomes of B. boettgeri and B. kuatunensis were sequenced, assembled, and annotated using next-generation sequencing. The length of mitochondrial genomes of B. boettgeri and B. kuatunensis was 16,597 and 17,921 bp, respectively, with both containing 13 protein coding genes, 22 tRNA genes, two rRNA genes, and one putative control region. Phylogenetic relationships based on protein-coding mitochondrial genes showed that the two Boulenophrys species formed a cluster with other Boulenophrys species. The two new sequences provide valuable insights into the mitochondrial genomes of these two species, offering important data for understanding the phylogenetic relationships of the genus Boulenophrys.

Morphology, complete mitochondrial genome, and molecular phylogeny of Rhabdias macrocephalum n. sp. (Nematoda: Rhabdiasidae) from Diploderma splendidum (Reptilia: Agamidae)

Species of the genus Rhabdias Stiles & Hassall, 1905 are common parasitic nematodes occurring in the lungs of amphibians and reptiles worldwide. In the present study, Rhabdias macrocephalum n. sp. is described using integrated morphological methods (light and scanning electron microscopy) and molecular approaches (sequencing of the nuclear 28S and ITS regions, and mitochondrial cox1, cox2, and 12S genes) based on specimens collected from the green striped tree dragon Diploderma splendidum (Barbour & Dunn) (Reptilia: Agamidae) in China. The complete mitochondrial genome of R. macrocephalum n. sp. was sequenced and annotated: it is 14,819 bp in length, including 12 protein coding genes (missing atp8), 22 tRNA genes, 2 rRNA genes and three non-coding regions. The gene arrangement of R. macrocephalum n. sp. is different from all of the currently available mitogenomes of nematodes and represents a novel type of mitochondrial gene arrangement reported in Nematoda. Molecular phylogenetic results based on the ITS + 28S data support the monophyly of Entomelas, Pneumonema, Serpentirhabdias, and Rhabdias, and showed R. macrocephalum n. sp. forming a most basal lineage in Rhabdias.

The mitochondrial genome of Huananpotamon koatenense (Rathbun, 1904) (Brachyura, Potamidae) and phylogenetic analysis

Freshwater crabs play essential roles in the well-functioning of the inland aquatic ecosystems. However, due to the lack of sufficient molecular resources, the study of freshwater crabs has been greatly hindered. In this study, the mitochondrial genome of Huananpotamon koatenense, a freshwater crab endemic to China, was sequenced for the first time. This mitogenome sequence is 15,528 bp long, and contains 13 protein-coding genes, 2 rRNA genes and 22 tRNA genes. Phylogenetic analyses based on 25 mitogenomes showed that H. koatenense was clustered with the known congeneric species of H. lichuanense.

Taking advantage of reference-guided assembly in a slowly-evolving lineage: Application to Testudo graeca

BACKGROUND

Obtaining de novo chromosome-level genome assemblies greatly enhances conservation and evolutionary biology studies. For many research teams, long-read sequencing technologies (that produce highly contiguous assemblies) remain unaffordable or unpractical. For the groups that display high synteny conservation, these limitations can be overcome by a reference-guided assembly using a close relative genome. Among chelonians, tortoises (Testudinidae) are considered one of the most endangered taxa, which calls for more genomic resources. Here we make the most of high synteny conservation in chelonians to produce the first chromosome-level genome assembly of the genus Testudo with one of the most iconic tortoise species in the Mediterranean basin: Testudo graeca.

RESULTS

We used high-quality, paired-end Illumina sequences to build a reference-guided assembly with the chromosome-level reference of Gopherus evgoodei. We reconstructed a 2.29 Gb haploid genome with a scaffold N50 of 107.598 Mb and 5.37% gaps. We sequenced 25,998 protein-coding genes, and identified 41.2% of the assembly as repeats. Demographic history reconstruction based on the genome revealed two events (population decline and recovery) that were consistent with previously suggested phylogeographic patterns for the species. This outlines the value of such reference-guided assemblies for phylogeographic studies.

CONCLUSIONS

Our results highlight the value of using close relatives to produce de novo draft assemblies in species where such resources are unavailable. Our annotated genome of T. graeca paves the way to delve deeper into the species' evolutionary history and provides a valuable resource to enhance direct conservation efforts on their threatened populations.

The first mitogenome of the subfamily Stenoponiinae (Siphonaptera: Ctenophthalmidae) and implications for its phylogenetic position

Fleas are the most important insect vectors that parasitize warm-blooded animals and are known vectors of zoonotic pathogens. A recent study showed that Stenoponia polyspina parasitizing Eospalax baileyi in Zoige County have carried Bartonella spp. and Spotted fever group Rickettsia (SFGR). Accurate identification and differentiation of fleas are essential for prevention and control of zoonotic pathogens. To understand phylogenetic relationship of the subfamily Stenoponiinae, we described morphological characteristics of S. polyspina and sequenced its mitogenome with 14,933 bp in size and high A + T content (~ 79%). The S. polyspina mitogenome retained the ancestral pattern of mitochondrial gene arrangement of arthropods without rearrangement. The start codons of 13 protein-coding genes (PCGs) are traditional ATN and the stop codons are TAA or TAG. Anticodon loops of all tRNA genes were 7 bp except for trnL2 and trnD had anticodon loops with 9 bp and the abnormal anticodon loops may be associated with frameshifting mutation. Genetic distance and Ka/Ks ratios indicated that all 13 PCGs of S. polyspina were subjected to purifying selection, with cox1 at the slowest rate and atp8 at the fastest rate. The mitogenomes of 24 species representing 7 families in the order Siphonaptera were selected to reconstruct phylogenetic tree based on concatenated nucleotide sequences of two datasets (PCGRNA matrix and PCG12RNA matrix) using Bayesian inference (BI) and Maximum likelihood (ML) methods. Phylogenetic tree supported that the superfamilies Ceratophylloidea, Vermipsylloidea, Pulicoidea were monophyletic and the superfamily Hystrichopsylloidea was paraphyletic. The family Ctenophthalmidae was monophyletic in PCGRNA-ML (codon partition) and paraphyletic in the remain trees. S. polyspina belongs to the subfamily Stenoponiinae was closely more related to the subfamily Rhadinopsyllinae. This paper explored phylogenetic position of diverse clades within the order Siphonaptera based on morphological and mitogenome data of S. polyspina. Our research enriched NCBI database of the order Siphonaptera.

Chromosome-level genome assembly of marmalade hoverfly Episyrphus balteatus (Diptera: Syrphidae)

Episyrphus balteatus can provide dual ecosystem services including pest control and pollination, which the larvae are excellent predators of aphid pest whereas adults are efficient pollinator. In this study, we assembled a high-quality genome of E. balteatus from northern China geographical population at the chromosome level by using Illumina, PacBio long reads, and Hi-C technologies. The 467.42 Mb genome was obtained from 723 contigs, with a contig N50 of 9.16 Mb and Scaffold N50 of 118.85 Mb, and 90.25% (431.75 Mb) of the assembly was anchored to 4 pseudo-autosomes and one pseudo-heterosome. In total, 14,848 protein-coding genes were annotated, and 95.14% of genes were fully represented in NR, GO, KEGG databases. Besides, we also obtained the mitochondrial genome of E. balteatus of 16, 837 bp in length with 37 typical mitochondrial genes. Overall, this high-quality genome is valuable for evolutionary and genetic studies of E. balteatus and other Syrphidae hoverfly species.

Identification of a longevity gene through evolutionary rate covariation of insect mito-nuclear genomes

Oxidative phosphorylation, essential for energy metabolism and linked to the regulation of longevity, involves mitochondrial and nuclear genes. The functions of these genes and their evolutionary rate covariation (ERC) have been extensively studied, but little is known about whether other nuclear genes not targeted to mitochondria evolutionarily and functionally interact with mitochondrial genes. Here we systematically examined the ERC of mitochondrial and nuclear benchmarking universal single-copy ortholog (BUSCO) genes from 472 insects, identifying 75 non-mitochondria-targeted nuclear genes. We found that the uncharacterized gene CG11837-a putative ortholog of human DIMT1-regulates insect lifespan, as its knockdown reduces median lifespan in five diverse insect species and Caenorhabditis elegans, whereas its overexpression extends median lifespans in fruit flies and C. elegans and enhances oxidative phosphorylation gene activity. Additionally, DIMT1 overexpression protects human cells from cellular senescence. Together, these data provide insights into the ERC of mito-nuclear genes and suggest that CG11837 may regulate longevity across animals.

Genome sequencing of Coryphaenoides yaquinae reveals convergent and lineage-specific molecular evolution in deep-sea adaptation

Abyssal (3501-6500 m) and hadal (>6500 m) fauna evolve under harsh abiotic stresses, characterized by high hydrostatic pressure, darkness and food shortage, providing unique opportunities to investigate mechanisms underlying environmental adaptation. Genomes of several hadal species have recently been reported. However, the genetic adaptation of deep sea species across a broad spectrum of ocean depths has yet to be thoroughly investigated, due to the challenges imposed by collecting the deep sea species. To elucidate the correlation between genetic innovation and vertical distribution, we generated a chromosome-level genome assembly of the macrourids Coryphaenoides yaquinae, which is widely distributed in the abyssal/hadal zone ranging from 3655 to 7259 m in depth. Genomic comparisons among shallow, abyssal and hadal-living species identified idiosyncratic and convergent genetic alterations underlying the extraordinary adaptations of deep-sea species including light perception, circadian regulation, hydrostatic pressure and hunger tolerance. The deep-sea fishes (Coryphaenoides Sp. and Pseudoliparis swirei) venturing into various ocean depths independently have undergone convergent amino acid substitutions in multiple proteins such as rhodopsin 1, pancreatic and duodenal homeobox 1 and melanocortin 4 receptor which are known or verified in zebrafish to be related with vision adaptation and energy expenditure. Convergent evolution events were also identified in heat shock protein 90 beta family member 1 and valosin-containing protein genes known to be related to hydrostatic pressure adaptation specifically in fishes found around the hadal range. The uncovering of the molecular convergence among the deep-sea species shed new light on the common genetic innovations required for deep-sea adaptation by the fishes.

The complete mitochondrial genome of Exoristobia philippinensis (Hymenoptera: Chalcidoidea: Encyrtidae) and phylogenetic analysis

Exoristobia philippinensis (Hymenoptera: Encyrtidae) is a worldwide parasitic wasp. This work presents the mitochondrial genome (mitogenome) of E. philippinensis for the first time. The complete mitochondrial genome of E. philippinensis was sequenced and annotated, which was 15,751 bp in length, and encoded 13 protein-coding genes (PCGs), 22 transfer RNA genes (tRNAs), and two ribosomal RNA genes (rRNAs). All 13 PCGs were initiated by the ATN (ATG, ATT, and ATA) codon, terminated with the stop codon TAA except for ND1 which ends with TAG. Phylogenetic analysis showed that E. philippinensis has a sister relationship with the genus Lamennaisia.

Phylogenetic and Comparative Genomics Study of Papilionidae Based on Mitochondrial Genomes

Most species of Papilionidae are large and beautiful ornamental butterflies. They are recognized as model organisms in ecology, evolutionary biology, genetics, and conservation biology but present numerous unresolved phylogenetic problems. Complete mitochondrial genomes (mitogenomes) have been widely used in phylogenetic studies of butterflies, but mitogenome knowledge within the family Papilionidae is limited, and its phylogeny is far from resolved. In this study, we first report the mitogenome of Byasa confusa from the subfamily Papilioninae of Papilionidae. The mitogenome of B. confusa is 15,135 bp in length and contains 13 protein-coding genes, 22 transfer RNA genes, 2 ribosomal RNA genes, and an AT-rich control region (CR), closely mirroring the genomic structure observed in related butterfly species. Comparative analysis of 77 Papilionidae mitogenomes shows gene composition and order to be identical to that of an ancestral insect, and the AT bias, Ka/Ks, and relative synonymous codon usage (RSCU) are all consistent with that of other reported butterfly mitogenomes. We conducted phylogenetic analyses using maximum-likelihood (ML) and Bayesian-inference (BI) methods, with 77 Papilionidae species as ingroups and two species of Nymphalidae and Lycaenidae as outgroups. The phylogenetic analysis indicated that B. confusa were clustered within Byasa. The phylogenetic trees show the monophyly of the subfamily Papilioninae and the tribes Leptocircini, Papilionini, and Troidini. The data supported the following relationships in tribe level on Papilioninae: (((Troidini + Papilionini) + Teinopalpini) + Leptocircini). The divergence time analysis suggests that Papilionidae originated in the late Creataceous. Overall, utilizing the largest number of Papilionidae mitogenomes sequenced to date, with the current first exploration in a phylogenetic analysis on Papilionidae (including four subfamilies), this study comprehensively reveals the mitogenome characteristics and mitogenome-based phylogeny, providing information for further studies on the mitogenome, phylogeny, evolution, and taxonomic revision of the Papilionidae family.

Mitogenomic phylogeny, biogeography, and cryptic divergence of the genus Silurus (Siluriformes: Siluridae)

The genus Silurus, an important group of catfish, exhibits heterogeneous distribution in Eurasian freshwater systems. This group includes economically important and endangered species, thereby attracting considerable scientific interest. Despite this interest, the lack of a comprehensive phylogenetic framework impedes our understanding of the mechanisms underlying the extensive diversity found within this genus. Herein, we analyzed 89 newly sequenced and 20 previously published mitochondrial genomes (mitogenomes) from 13 morphological species to reconstruct the phylogenetic relationships, biogeographic history, and species diversity of Silurus. Our phylogenetic reconstructions identified eight clades, supported by both maximum-likelihood and Bayesian inference. Sequence-based species delimitation analyses yielded multiple molecular operational taxonomic units (MOTUs) in several taxa, including the Silurus asotus complex (four MOTUs) and Silurus microdorsalis (two MOTUs), suggesting that species diversity is underestimated in the genus. A reconstructed time-calibrated tree of Silurus species provided an age estimate of the most recent common ancestor of approximately 37.61 million years ago (Ma), with divergences among clades within the genus occurring between 11.56 Ma and 29.44 Ma, and divergences among MOTUs within species occurring between 3.71 Ma and 11.56 Ma. Biogeographic reconstructions suggested that the ancestral area for the genus likely encompassed China and the Korean Peninsula, with multiple inferred dispersal events to Europe and Central and Western Asia between 21.78 Ma and 26.67 Ma and to Japan between 2.51 Ma and 18.42 Ma. Key factors such as the Eocene-Oligocene extinction event, onset and intensification of the monsoon system, and glacial cycles associated with sea-level fluctuations have likely played significant roles in shaping the evolutionary history of the genus Silurus.

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

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

Draft genome assemblies of the ponerine ant Odontoponera transversa and the carpenter ant Camponotus friedae (Hymenoptera: Formicidae)

OBJECTIVES

Ants are ecologically dominant insects in most terrestrial ecosystems, with more than 14,000 extant species in about 340 genera recorded to date. However, genomic resources are still scarce for most species, especially for species endemic in East or Southeast Asia, limiting the study of phylogeny, speciation and adaptation of this evolutionarily successful animal lineage. Here, we assemble and annotate the genomes of Odontoponera transversa and Camponotus friedae, two ant species with a natural distribution in China, to facilitate future study of ant evolution.

DATA DESCRIPTION

We obtained a total of 16 Gb and 51 Gb PacBio HiFi data for O. transversa and C. friedae, respectively, which were assembled into the draft genomes of 339 Mb for O. transversa and 233 Mb for C. friedae. Genome assessments by multiple metrics showed good completeness and high accuracy of the two assemblies. Gene annotations assisted by RNA-seq data yielded a comparable number of protein-coding genes in the two genomes (10,892 for O. transversa and 11,296 for C. friedae), while repeat annotations revealed a remarkable difference of repeat content between these two ant species (149.4 Mb for O. transversa versus 49.7 Mb for C. friedae). Besides, complete mitochondrial genomes for the two species were assembled and annotated.