Improved systematic tRNA gene annotation allows new insights into the evolution of mitochondrial tRNA structures and into the mechanisms of mitochondrial genome rearrangements

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.

Complete mitochondrial genomes of three vulnerable cave bat species and their phylogenetic relationships within the order Chiroptera

The IUCN Red List of Threatened Species contains 175 Brazilian bat species that are threatened by extinction in some degree. From this perspective, it is essential to expand the knowledge about the genetic diversity of vulnerable bats. Genomic sequencing can be useful to generate robust and informative genetic references, increasing resolution when analyzing relationships among populations, species, or higher taxonomic levels. In this study, we sequenced and characterized in detail the first complete mitochondrial genomes of Furipterus horrens, Lonchorhina aurita, and Natalus macrourus, and investigated their phylogenetic position based on amino acid sequences of protein-coding genes (PCGs). The mitogenomes of these species are 16,516, 16,697, and 16,668 bp in length, respectively, and each comprises 13 PCGs, 22 tRNA genes, two rRNA genes, and a putative control region (CR). In the three species, genes were arranged similarly to all other previously described bat mitogenomes, and nucleotide composition was also consistent with the reported range. The length and arrangement of rrnS and rrnL were also consistent with those of other bat species, showing a positive AT-skew and a negative GC-skew. Except for trnS1, for which we did not observe the DHU arm, all other tRNAs showed the cloverleaf secondary structure in the three species. In addition, the mitogenomes showed minor differences in start and stop codons, and in all PCGs, codons ending in adenine were more common compared to those ending in guanine. We found that PCGs of the three species use multiple codons to encode each amino acid, following the previously documented pattern. Furthermore, all PCGs are under purifying selection, with atp8 experiencing the most relaxed purifying selection. Considering the phylogenetic reconstruction, F. horrens was recovered as sister to Noctilio leporinus, L. aurita and Tonatia bidens shared a node within Phyllostomidae, and N. macrourus appeared as sister to Molossidae and Vespertilionidae.

Comparative mitochondrial genomics of endemic Mexican vesper yellow bats genus Rhogeessa (Chiroptera: Vespertilionidae) and insights into internal relationships in the family Vespertilionidae

In the species-rich family Vespertilionidae, vesper yellow bats in the genus Rhogeessa include eleven species, three of them endemic to Mexico. These insectivorous bats provide important ecosystem services, including pest control. Even though some aspects of their biology are well- known, only a few genomic resources are available for these species, which limits our understanding of their biology. In this study, we assembled and annotated the mitochondrial genome of four species: R. aenea, R. genowaysi, R. mira, and R. parvula. We generated a phylomitogenomic hypothesis based on translated protein-coding genes for a total of 52 species in the family Vespertilionidae and examined the phylogenetic position of the genus Rhogeessa and species within the family. The AT-rich mitogenomes of R. aenea, R. genowaysi, R. mira, and R. parvula are 16,763, 16,781, 16,807, and 16,794 pb in length, respectively. Each studied mitogenome encodes 13 Protein Coding Genes (PCGs), 22 transfer RNA genes, and 2 rRNA genes, and contains a putative control region (CR). All tRNAs exhibit a 'cloverleaf' secondary structure, except tRNA-Serine-1 that lacked the DHU arm in all studied mitogenomes. Selective pressure analyses indicated that all protein-coding genes are exposed to purifying selection. The phylomitogenomic analysis supported the monophyletic status of the family Vespertilionidae, confirmed the placement of Rhogeessa within the tribe Antrozoini, and clarified phylogenetic relationships within and among subfamilies and tribes in this family. Our results indicate that phylomitogenomics are useful to explore the evolutionary history of vesper bats. The assembly and comprehensive analysis of mitochondrial genomes offer the potential to generate molecular references and resources beneficial for genetic analyses aimed at understanding the ecology and evolution of these remarkable bats.

Comparative mitogenomics of marine angelfishes (F: Pomacanthidae)

The targeted capture of ultraconserved elements (UCEs) has substantially increased the amount of genetic data available for phylogenomic reconstructions. These capture datasets frequently contain mitochondrial DNA as a by-product, often in the form of complete mitogenomes. These can be efficiently harvested to expand existing datasets without additional costs. Here, we present new mitochondrial genomes for six marine angelfish species (F: Pomacanthidae), assembled and annotated from off-target UCE reads. We provide the first comparative analysis of all mitochondrial genomes available for the Pomacanthidae. Results showed that the average length of pomacanthid mitogenomes is 16.8 kbp. Total GC and AT content varied between 44.5% and 46.3%, and 53.7% and 55.5%, respectively. The architecture of angelfish mitogenomes was comparable to that seen in other fish species with 13 protein-coding genes (PCGs), 22 transfer RNA genes, two ribosomal RNA genes and the control region. All 13 PCGs evolved under purifying selection, highlighting a high level of selection pressure and gene expression to preserve genetic integrity. The ND6 and ATP8 genes had the highest ratio of non-synonymous (dN) to synonymous (dS) substitutions, indicating a relaxation of purifying selection constraints. Finally, these newly assembled mitogenomes will allow further investigations of the population genetics, systematics and evolutionary biology of one of the most prominent reef fish family in the aquarium trade.

Detailed characterization of the complete mitochondrial genome of the oceanic whitetip shark Carcharhinus longimanus (Poey, 1861)

BACKGROUND

The oceanic whitetip shark Carcharhinus longimanus (family Carcharhinidae) is one of the largest sharks inhabiting all tropical and subtropical oceanic regions. Due to their life history traits and mortality attributed to pelagic longline fishing practices, this species is experiencing substantial population decline. Currently, C. longimanus is considered by the IUCN Red List of Threatened Species as "vulnerable" throughout its range and "critically endangered" in the western north Atlantic. This study sequences and describes the complete mitochondrial genome of C. longimanus in detail.

METHODS AND RESULTS

The mitochondrial genome of C. longimanus was assembled through next-generation sequencing and then analyzed using specialized bioinformatics tools. The circular, double-stranded AT-rich mitogenome of C. longimanus is 16,704 bp long and contains 22 tRNA genes, 2 rRNA genes, 13 protein coding genes and a 1,065 bp long control region (CR). Out of the 22 tRNA genes, only one (tRNA-Ser1) lacked a typical 'cloverleaf' secondary structure. The prevalence of TTA (Leu), ATT (Ile) and CTA (Leu) codons in the PCGs likely contributes to the AT-rich nature of this mitogenome. In the CR, ten microsatellites were detected but no tandem repeats were found. Stem-and-loop secondary structures were common along the entire length of the CR. Ka/Ks values estimated for all PCGs were < 1, indicating that all the PCGs experience purifying selection. A phylomitogenomic analysis based on translated PCGs confirms the sister relationship between C. longimanus and C. obscurus. The analysis did not support the monophyly of the genus Carcharhinus.

CONCLUSIONS

The assembled mitochondrial genome of this pelagic shark can provide insight into the phylogenetic relationships in the genus Carcharhinus and aid conservation and management efforts in the Central Pacific Ocean.

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.

Complete mitochondrial genome of the Minois paupera Alphéraky, 1888 (nymphalidae: satyrinae) and its phylogenetic analysis

The present study firstly reported a complete mitochondrial genome of Minois paupera (Alphéraky, 1888), a Satyrinae species endemic to China. This mitogenome is circular, 15,213 bp in length, and consists of 37 typical mitochondrial genes, including 13 protein-coding genes (PCGs), 22 tRNAs, and two rRNAs. The phylogenetic position was inferred using 31 previously published complete mitogenomes, and the results reveal that M. paupera is the most closely related to Minois dryas. The complete mitogenome of M. paupera provides useful genetic information for further research on the phylogeography and phylogeny of the genus Minois.

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.

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.

Characterization of the complete mitochondrial genome and phylogenetic analyses of Haemaphysalis tibetensis Hoogstraal, 1965 (Acari: Ixodidae)

Ticks are specialized ectoparasites that feed on blood, causing physical harm to the host and facilitating pathogen transmission. The genus Haemaphysalis contains vectors for numerous infectious agents. These agents cause various diseases in humans and animals. Mitochondrial genome sequences serve as reliable molecular markers, forming a crucial basis for evolutionary analyses, studying species origins, and exploring molecular phylogeny. We extracted mitochondrial genome from the enriched mitochondria of Haemaphysalis tibetensis and obtained a 14,714-bp sequence. The mitochondrial genome consists of 13 protein-coding genes (PCGs), two ribosomal RNA, 22 transfer RNAs (tRNAs), and two control regions. The nucleotide composition of H. tibetensis mitochondrial genome was 38.38 % for A, 9.61 % for G, 39.32 % for T, and 12.69 % for C. The A + T content of H. tibetensis mitochondrial genome was 77.7 %, significantly higher than the G + C content. The repeat units of H. tibetensis exhibited two identical repeat units of 33 bp in length, positioned downstream of nad1 and rrnL genes. Furthermore, phylogenetic analyses based on the 13 PCGs indicated that Haemaphysalis tibetensis (subgenus Allophysalis) formed a monophyletic clade with Haemaphysalis nepalensis (subgenus Herpetobia) and Haemaphysalis danieli (subgenus Allophysalis). Although the species Haemaphysalis inermis, Haemaphysalis kitaokai, Haemaphysalis kolonini, and Haemaphysalis colasbelcouri belong to the subgenus Alloceraea, which were morphologically primitive hemaphysalines just like H. tibetensis, these four tick species cannot form a single clade with H. tibetensis. In this study, the whole mitochondrial genome sequence of H. tibetensis from Tibet was obtained, which enriched the mitochondrial genome data of ticks and provided genetic markers to study the population heredity and molecular evolution of the genus Haemaphysalis.

Insights into the nuclear and mitochondrial genome of the Lemon shark Negaprion brevirostris using low-coverage sequencing: Genome size, repetitive elements, mitochondrial genome, and phylogenetic placement

The Lemon shark Negaprion brevirostris is an important species experiencing conservation issues that is in need of genomic resources. Herein, we conducted a genome survey sequencing in N. brevirostris and determined genome size, explored repetitive elements, assembled and annotated the 45S rRNA DNA operon, and assembled and described in detail the mitochondrial genome. Lastly, the phylogenetic position of N. brevirostris in the family Carcharhinidae was examined using translated protein coding genes. The estimated haploid genome size ranged between 2.29 and 2.58 Gbp using a k-mer analysis, which is slightly below the genome size estimated for other sharks belonging to the family Carcharhinidae. Using a k-mer analysis, approx. 64-71 % of the genome of N. brevirostris was composed of repetitive elements. A relatively large proportion of the 'repeatome' could not be annotated. Taking into account only annotated repetitive elements, Class I - Long Interspersed Nuclear Element (LINE) were the most abundant repetitive elements followed by Class I - Penelope and Satellite DNA. The nuclear ribosomal operon was fully assembled. The AT-rich complete mitochondrial genome was 16,703 bp long and encoded 13 protein coding genes, 2 ribosomal RNA genes, and 22 transfer RNA genes. Negaprion brevirostris is closely related to the genera Carcharhinus, Glyphis and Lamiopsis in the family Carcharinidae. This new genomic resources will aid with the development of conservation plans for this large coastal shark.

Comparative analyses of mitogenomes in the social bees with insights into evolution of long inverted repeats in the Meliponini

The insect mitogenome is typically a compact circular molecule with highly conserved gene contents. Nonetheless, mitogenome structural variations have been reported in specific taxa, and gene rearrangements, usually the tRNAs, occur in different lineages. Because synapomorphies of mitogenome organizations can provide information for phylogenetic inferences, comparative analyses of mitogenomes have been given increasing attention. However, most studies use a very few species to represent the whole genus, tribe, family, or even order, overlooking potential variations at lower taxonomic levels, which might lead to some incorrect inferences. To provide new insights into mitogenome organizations and their implications for phylogenetic inference, this study conducted comparative analyses for mitogenomes of three social bee tribes (Meliponini, Bombini, and Apini) based on the phylogenetic framework with denser taxonomic sampling at the species and population levels. Comparative analyses revealed that mitogenomes of Apini and Bombini are the typical type, while those of Meliponini show diverse variations in mitogenome sizes and organizations. Large inverted repeats (IRs) cause significant gene rearrangements of protein coding genes (PCGs) and rRNAs in Indo-Malay/Australian stingless bee species. Molecular evolution analyses showed that the lineage with IRs have lower d N/ d S ratios for PCGs than lineages without IRs, indicating potential effects of IRs on the evolution of mitochondrial genes. The finding of IRs and different patterns of gene rearrangements suggested that Meliponini is a hotspot in mitogenome evolution. Unlike conserved PCGs and rRNAs whose rearrangements were found only in the mentioned lineages within Meliponini, tRNA rearrangements are common across all three tribes of social bees, and are significant even at the species level, indicating that comprehensive sampling is needed to fully understand the patterns of tRNA rearrangements, and their implications for phylogenetic inference.

Novel insight into lepidopteran phylogenetics from the mitochondrial genome of the apple fruit moth of the family Argyresthiidae

BACKGROUND

The order Lepidoptera has an abundance of species, including both agriculturally beneficial and detrimental insects. Molecular data has been used to investigate the phylogenetic relationships of major subdivisions in Lepidoptera, which has enhanced our understanding of the evolutionary relationships at the family and superfamily levels. However, the phylogenetic placement of many superfamilies and/or families in this order is still unknown. In this study, we determine the systematic status of the family Argyresthiidae within Lepidoptera and explore its phylogenetic affinities and implications for the evolution of the order. We describe the first mitochondrial (mt) genome from a member of Argyresthiidae, the apple fruit moth Argyresthia conjugella. The insect is an important pest on apples in Fennoscandia, as it switches hosts when the main host fails to produce crops.

RESULTS

The mt genome of A. conjugella contains 16,044 bp and encodes all 37 genes commonly found in insect mt genomes, including 13 protein-coding genes (PCGs), two ribosomal RNAs, 22 transfer RNAs, and a large control region (1101 bp). The nucleotide composition was extremely AT-rich (82%). All detected PCGs (13) began with an ATN codon and terminated with a TAA stop codon, except the start codon in cox1 is ATT. All 22 tRNAs had cloverleaf secondary structures, except trnS1, where one of the dihydrouridine (DHU) arms is missing, reflecting potential differences in gene expression. When compared to the mt genomes of 507 other Lepidoptera representing 18 superfamilies and 42 families, phylogenomic analyses found that A. conjugella had the closest relationship with the Plutellidae family (Yponomeutoidea-super family). We also detected a sister relationship between Yponomeutoidea and the superfamily Tineidae.

CONCLUSIONS

Our results underline the potential importance of mt genomes in comparative genomic analyses of Lepidoptera species and provide valuable evolutionary insight across the tree of Lepidoptera species.

A comprehensive phylogenomic study unveils evolutionary patterns and challenges in the mitochondrial genomes of Carcharhiniformes: A focus on Triakidae

The complex evolutionary patterns in the mitochondrial genome (mitogenome) of the most species-rich shark order, the Carcharhiniformes (ground sharks) has led to challenges in the phylogenomic reconstruction of the families and genera belonging to the order, particularly the family Triakidae (houndsharks). The current state of Triakidae phylogeny remains controversial, with arguments for both monophyly and paraphyly within the family. We hypothesize that this variability is triggered by the selection of different a priori partitioning schemes to account for site and gene heterogeneity within the mitogenome. Here we used an extensive statistical framework to select the a priori partitioning scheme for inference of the mitochondrial phylogenomic relationships within Carcharhiniformes, tested site heterogeneous CAT + GTR + G4 models and incorporated the multi-species coalescent model (MSCM) into our analyses to account for the influence of gene tree discordance on species tree inference. We included five newly assembled houndshark mitogenomes to increase resolution of Triakidae. During the assembly procedure, we uncovered a 714 bp-duplication in the mitogenome of Galeorhinus galeus. Phylogenetic reconstruction confirmed monophyly within Triakidae and the existence of two distinct clades of the expanded Mustelus genus. The latter alludes to potential evolutionary reversal of reproductive mode from placental to aplacental, suggesting that reproductive mode has played a role in the trajectory of adaptive divergence. These new sequences have the potential to contribute to population genomic investigations, species phylogeography delineation, environmental DNA metabarcoding databases and, ultimately, improved conservation strategies for these ecologically and economically important species.

The Phylogenetic Relationship of Lamiinae (Coleoptera: Cerambycidae) Using Mitochondrial Genomes

Lamiinae is the largest subfamily of the Cerambycidae (longhorn beetles), with approximately 21,863 described species. Previous phylogenetic studies of Lamiinae showed that this subfamily was monophyletic, but the relationship between the tribes of Lamiinae is still controversial. Partial molecular data and species morphological characteristics are not sufficient to resolve species phylogenetic studies perfectly. At the same time, the full mitochondrial genome contains more comprehensive genetic data. Benefiting from the development of next-generation sequencing (NGS), mitochondrial genomes can be easily acquired and used as reliable molecular markers to investigate phylogenetic relationships within Cerambycidae. Using NGS technology, we obtained 11 mitochondrial genome sequences of Lamiinae species. Based on this newly generated mitochondrial genome dataset matrix, we reconstructed the phylogeny of Lamiinae. The Bayesian Inference and Maximum Likelihood analyses strongly support the monophyly of four tribes (Lamiini, Batocerini, Mesosini, and Saperdini), whereas the tribe Acanthocinini was identified as paraphyletic. Other mitochondrial structural features were also observed: the start codon in the nad1 gene of all 11 mitochondrial genomes is TTG; 17-22 bp intergenic spacers (IGS) with a 'TACTA' motif were found between trnS2 and nad1. Moreover, two long IGS were found in Mesosa myops and Batocera sp. Tandem repeats were found in the IGS of Batocera sp.

The complete mitochondrial genome of Echinolaelaps fukienensis provide insights into phylogeny and rearrangement in the superfamily Dermanyssoidea

BACKGROUND

Echinolaelaps fukienensis is the dominant mite species parasitic on the body surface of the genus Niviventer. The mitochondrial genome (mitogenome) has its own independent genetic material and genetic system, and is now widely used in population genetics, genealogical biogeography, phylogeny and molecular evolution studies. Species diversity of the superfamily Dermanyssoidea is very rich, but its mitogenomes AT content is high, and it is difficult to amplify the complete mitogenome by routine PCR. To date, we have only obtained the mitogenomes of 6 species, scarcity on sequence data has greatly impeded the studies in the superfamily Dermanyssoidea.

METHODS

Echinolaelaps fukienensis were collected in 2019 from the body surface of Niviventer confucianus (Rodentia, Muridae) in Yunnan Province. The E. fukienensis mitogenome was determined and analyzed for the first time using the Illumina Novoseq 6000 platform. Phylogenetic analyses of the superfamily Dermanyssoidea were conducted based on the entire mitogenome sequences.

RESULTS

The E. fukienensis mitogenome was 14,402 bp, which is known the smallest genome of the superfamily Dermanyssoidea, encoding a total of 37 genes, including 13 PCGs, 22 tRNAs, 2 rRNAs and 1 control region. Most protein-coding genes use ATN as the start codon and TAN as the stop codon. AT and GC skew of atp8 genes in E. fukienensis were both 0. The average length of 22 tRNA genes of E. fukienensis was 64 bp, and secondary structures of tRNAs showed base mismatches and missing D-arms in many places. Compared with gene arrangement pattern of the hypothetical ancestor of arthropods, the E. fukienensis mitogenome shows a novel arrangement pattern. Phylogenetic tree supported the monophyly of the superfamily Dermanyssoidea. Echinolaelaps fukienensis being the least genetic distant (0.2762) and most closely related to Varroa destructor.

CONCLUSIONS

This study analyzed comprehensive the structure and evolution of the E. fukienensis mitogenome for the first time, enriches molecular data of the genus Echinolaelaps, which will contribute to further understand phylogeny and rearrangement patterns of the superfamily Dermanyssoidea.

Whole genome assembly and annotation of the endangered Caribbean coral Acropora cervicornis

Coral species in the genus Acropora are key ecological components of coral reefs worldwide and represent the most diverse genus of scleractinian corals. While key species of Indo-Pacific Acropora have annotated genomes, no annotated genome has been published for either of the two species of Caribbean Acropora. Here we present the first fully annotated genome of the endangered Caribbean staghorn coral, Acropora cervicornis. We assembled and annotated this genome using high-fidelity nanopore long-read sequencing with gene annotations validated with mRNA sequencing. The assembled genome size is 318 Mb, with 28,059 validated genes. Comparative genomic analyses with other Acropora revealed unique features in A. cervicornis, including contractions in immune pathways and expansions in signaling pathways. Phylogenetic analysis confirms previous findings showing that A. cervicornis diverged from Indo-Pacific relatives around 41 million years ago, with the closure of the western Tethys Sea, prior to the primary radiation of Indo-Pacific Acropora. This new A. cervicornis genome enriches our understanding of the speciose Acropora and addresses evolutionary inquiries concerning speciation and hybridization in this diverse clade.

Exploring the mitochondrial genome of Caridina pseudogracilirostris: a comparative analysis within the Atyidae Family

BACKGROUND

Caridina pseudogracilirostris is a highly adaptive estuarine species found in brackish waters and marshes along the southwestern and southern coastal regions of India.

METHODS AND RESULTS

The whole mitochondrial genome of C. pseudogracilirostris is 15,451 bp in length with 59.3% AT content and encodes 37 genes, including 22 tRNAs, 13 protein-coding genes, and two rRNAs, which are arranged in a distinctive pattern similar to most crustaceans. ML and BI methods were used for phylogenetic analysis of C. pseudogracilirostris clustered with other Caridina species, supporting the monophyly of the Caridina genus within the Atyidae family. The fully annotated mitochondrial genome of C. pseudogracilirostris was submitted to GenBank under accession number OQ534868.1.

CONCLUSIONS

We are the first to report on the C. pseudogracilirostris whole mitochondrial genome, which provides a valuable resource for future research on genetics, evolution, phylogenetics, etc., among Caridina species and other species. The phylogenetic investigation supports the monophyly of the Caridina genus within the Atyidae family and emphasizes the value of mitochondrial genome data in determining the evolutionary relationships among crustaceans.

The mitochondrial genome of the mountain wooly tapir, Tapirus pinchaque and a formal test of the effect of altitude on the adaptive evolution of mitochondrial protein coding genes in odd-toed ungulates

BACKGROUND

The harsh conditions of high-altitude environments are known to drive the evolution of physiological and morphological traits in endothermic animals. These conditions are expected to result in the adaptive evolution of protein coding genes encoded in mitochondrial genomes that are vital for the oxidative phosphorylation pathway. In this study, we formally tested for signatures of adaptive evolution on mitochondrial protein coding genes in Tapirus pinchaque and other odd-toed ungulates inhabiting high-elevation environments.

RESULTS

The AT-rich mitochondrial genome of T. pinchaque is 16,750 bp long. A phylomitogenomic analysis supports the monophyly of the genus Tapirus and families in the Perissodactyla. The ratio of non-synonymous to synonymous substitutions demonstrated that all mitochondrial genes undergo purifying selection in T. pinchaque and other odd ungulates living at high elevations. Over this negative background selection, Branch Models suggested that cox3 and nad6 might be undergoing stronger purifying selection than other mitochondrial protein coding genes. Furthermore, Site Models suggested that one and four sites in nad2 and nad5, respectively, could be experiencing positive selection. However, these results were supported by Likelihood Ratio Tests but not Bayesian Empirical Bayes posterior probabilities. Additional analyses (in DataMonkey) indicated a relaxation of selection strength in nad6, evidence of episodic diversifying selection in cob, and revealed episodic positive/diversifying selection signatures for two sites in nad1, and one site each in nad2 and nad4.

CONCLUSION

The mitochondrial genome of T. pinchaque is an important genomic resource for conservation of this species and this study contributes to the understanding of adaptive evolution of mitochondrial protein coding genes in odd-toed ungulates inhabiting high-altitude environments.

MitoHiFi: a python pipeline for mitochondrial genome assembly from PacBio high fidelity reads

BACKGROUND

 PacBio high fidelity (HiFi) sequencing reads are both long (15-20 kb) and highly accurate (> Q20). Because of these properties, they have revolutionised genome assembly leading to more accurate and contiguous genomes. In eukaryotes the mitochondrial genome is sequenced alongside the nuclear genome often at very high coverage. A dedicated tool for mitochondrial genome assembly using HiFi reads is still missing.

RESULTS

 MitoHiFi was developed within the Darwin Tree of Life Project to assemble mitochondrial genomes from the HiFi reads generated for target species. The input for MitoHiFi is either the raw reads or the assembled contigs, and the tool outputs a mitochondrial genome sequence fasta file along with annotation of protein and RNA genes. Variants arising from heteroplasmy are assembled independently, and nuclear insertions of mitochondrial sequences are identified and not used in organellar genome assembly. MitoHiFi has been used to assemble 374 mitochondrial genomes (368 Metazoa and 6 Fungi species) for the Darwin Tree of Life Project, the Vertebrate Genomes Project and the Aquatic Symbiosis Genome Project. Inspection of 60 mitochondrial genomes assembled with MitoHiFi for species that already have reference sequences in public databases showed the widespread presence of previously unreported repeats.

CONCLUSIONS

 MitoHiFi is able to assemble mitochondrial genomes from a wide phylogenetic range of taxa from Pacbio HiFi data. MitoHiFi is written in python and is freely available on GitHub ( https://github.com/marcelauliano/MitoHiFi ). MitoHiFi is available with its dependencies as a Docker container on GitHub (ghcr.io/marcelauliano/mitohifi:master).

Comparative mitochondrial genomics of American nectar-feeding long-nosed bats Leptonycteris spp. with insights into the phylogeny of the family Phyllostomidae

Among leaf-nosed bats (family Phyllostomidae), the genus Leptonycteris (subfamily Glossophaginae), contains three migratory and obligate nectar-dwelling species of great ecological and economic importance; the greater long-nosed bat L. nivalis, the lesser long-nosed bat L. yerbabuenae, and the southern long-nosed bat L. curasoae. According to the IUCN, the three species are categorized as 'vulnerable', 'endangered', and 'near threatened', respectively. In this study, we assembled and characterized in detail the mitochondrial genome of Leptonycteris spp. and examined the phylogenetic position of this genus in the family Phyllostomidae based on protein coding genes (PCGs). The mitogenomes of L. nivalis, L. curasoae, and L. yerbabuenae are 16,708, 16,758, and 16,729 bp in length and each encode 13 PCGs, 22 transfer RNA genes, 2 rRNA genes, and a putative control region (CR). Mitochondrial gene order is identical to that reported before for the family Phyllostomidae. All tRNAs exhibit a 'cloverleaf' secondary structure, except tRNA-Serine-1 that is missing the DHU arm in the three species. All PCGs are exposed to purifying selection with atp8 experiencing the most relaxed purifying selection as the ω ratio was higher for this gene compared to the other PCGs in each species. The CR of each species contains three functional domains: extended termination associated sequence (ETAS), Central, and conserved sequence block domain (CSB). A phylomitogenomic analysis revealed that Leptonycteris is monophyletic and most closely related to the genus Glossophaga. The analysis also supported the monophyly of the family Glossophaginae in the speciose family Phyllostomidae. The mitochondria characterization of these species provides relevant information to develop molecular markers for conservation purposes.

Detecting gene breakpoints in noisy genome sequences using position-annotated colored de-Bruijn graphs

BACKGROUND

Identifying the locations of gene breakpoints between species of different taxonomic groups can provide useful insights into the underlying evolutionary processes. Given the exact locations of their genes, the breakpoints can be computed without much effort. However, often, existing gene annotations are erroneous, or only nucleotide sequences are available. Especially in mitochondrial genomes, high variations in gene orders are usually accompanied by a high degree of sequence inconsistencies. This makes accurately locating breakpoints in mitogenomic nucleotide sequences a challenging task.

RESULTS

This contribution presents a novel method for detecting gene breakpoints in the nucleotide sequences of complete mitochondrial genomes, taking into account possible high substitution rates. The method is implemented in the software package DeBBI. DeBBI allows to analyze transposition- and inversion-based breakpoints independently and uses a parallel program design, allowing to make use of modern multi-processor systems. Extensive tests on synthetic data sets, covering a broad range of sequence dissimilarities and different numbers of introduced breakpoints, demonstrate DeBBI 's ability to produce accurate results. Case studies using species of various taxonomic groups further show DeBBI 's applicability to real-life data. While (some) multiple sequence alignment tools can also be used for the task at hand, we demonstrate that especially gene breaks between short, poorly conserved tRNA genes can be detected more frequently with the proposed approach.

CONCLUSION

The proposed method constructs a position-annotated de-Bruijn graph of the input sequences. Using a heuristic algorithm, this graph is searched for particular structures, called bulges, which may be associated with the breakpoint locations. Despite the large size of these structures, the algorithm only requires a small number of graph traversal steps.

The mitochondrial genomes of big-eared bats, Macrotus waterhousii and Macrotus californicus (Chiroptera: Phyllostomidae: Macrotinae)

In the species-rich family Phyllostomidae, the genus Macrotus ('big eared' bats) contains only two species; Macrotus waterhousii, distributed in western, central, and southern Mexico, Guatemala and some Caribbean Islands, and Macrotus californicus, distributed in the southwestern USA, and in the Baja California peninsula and the state of Sonora in Mexico. In this study, we sequenced and assembled the mitochondrial genome of Macrotus waterhousii and characterized in detail this genome and that of the congeneric M. californicus. Then, we examined the phylogenetic position of Macrotus in the family Phyllostomidae based on protein coding genes (PCGs). The AT-rich mitochondrial genomes of M. waterhousii and M. californicus are 16,792 and 16,691 bp long, respectively, and each encode 13 PCGs, 22 tRNA genes, 2 rRNA genes, and a putative non-coding control region 1,336 and 1,232 bp long, respectively. Mitochondrial synteny in Macrotus is identical to that reported before for all other cofamilial species. In the two studied species, all tRNAs exhibit a 'typical' cloverleaf secondary structure with the exception of trnS1, which lacks the D arm. A selective pressure analysis demonstrated that all PCGs are under purifying selection. The CR of the two species feature three domains previously reported in other mammals, including bats: extended terminal associated sequences (ETAS), central (CD), and conserved sequence block (CSB). A phylogenetic analysis based on the 13 mitochondrial PCGs demonstrated that Macrotus is monophyletic and the subfamily Macrotinae is a sister group of all remaining phyllostomids in our analysis, except Micronycterinae. The assembly and detailed analysis of these mitochondrial genomes represents a step further to continue improving the understanding of phylogenetic relationships within the species-rich family Phyllostomidae.

The Complete Mitochondrial Genome of Mytilisepta virgata (Mollusca: Bivalvia), Novel Gene Rearrangements, and the Phylogenetic Relationships of Mytilidae

The circular mitochondrial genome of Mytilisepta virgata spans 14,713 bp, which contains 13 protein-coding genes (PCGs), 2 ribosomal RNA genes, and 22 transfer RNA genes. Analysis of the 13 PCGs reveals that the mitochondrial gene arrangement of Mytilisepta is relatively conserved at the genus level. The location of the atp8 gene in Mytilisepta keenae differs from that of other species. However, compared with the putative molluscan ancestral gene order, M. virgata exhibits a high level of rearrangement. We constructed phylogenetic trees based on concatenated 12 PCGs from Mytilidae. As a result, we found that M. virgata is in the same clade as other Mytilisepta spp. The result of estimated divergence times revealed that M. virgata and M. keenae diverged around the early Paleogene period, although the oldest Mytilisepta fossil was from the late or upper Eocene period. Our results provide robust statistical evidence for a sister-group relationship within Mytilida. The findings not only confirm previous results, but also provide valuable insights into the evolutionary history of Mytilidae.

Complete Mitochondrial Genome of Piophila casei (Diptera: Piophilidae): Genome Description and Phylogenetic Implications

Piophila casei is a flesh-feeding Diptera insect that adversely affects foodstuffs, such as dry-cured ham and cheese, and decaying human and animal carcasses. However, the unknown mitochondrial genome of P. casei can provide information on its genetic structure and phylogenetic position, which is of great significance to the research on its prevention and control. Therefore, we sequenced, annotated, and analyzed the previously unknown complete mitochondrial genome of P. casei. The complete mt genome of P. casei is a typical circular DNA, 15,785 bp in length, with a high A + T content of 76.6%. It contains 13 protein-coding genes (PCG), 2 ribosomal RNA (rRNA) genes, 22 transfer RNA (tRNA) genes, and 1 control region. Phylogenetic analysis of 25 Diptera species was conducted using Bayesian and maximum likelihood methods, and their divergence times were inferred. The comparison of the mt genomes from two morphologically similar insects P. casei and Piophila megastigmata indicates a divergence time of 7.28 MYA between these species. The study provides a reference for understanding the forensic medicine, taxonomy, and genetics of P. casei.

Comparative analysis of mitochondrial genomes reveals family-specific architectures and molecular features in scorpions (Arthropoda: Arachnida: Scorpiones)

Scorpions are a group of arachnids with great evolutionary success that comprise more than 2,000 described species. Mitochondrial genomes have been little studied in this clade. We describe and compare different scorpion mitochondrial genomes and analyze their architecture and molecular characteristics. We assembled eight new scorpion mitochondrial genomes from transcriptomic datasets, annotated them, predicted the secondary structures of tRNAs, and compared the nucleotide composition, codon usage, and relative synonymous codon usage of 16 complete scorpion mitochondrial genomes. Lastly, we provided a phylogeny based on all mitochondrial protein coding genes. We characterized the mitogenomes in detail and reported particularities such as dissimilar synteny in the family Buthidae compared to other scorpions, unusual tRNA secondary structures, and unconventional start and stop codons in all scorpions. Our comparative analysis revealed that scorpion mitochondrial genomes exhibit different architectures and features depending on taxonomic identity. We highlight the parvorder Buthida, particularly the family Buthidae, as it invariably exhibited different mitogenome features such as synteny, codon usage, and AT-skew compared to the parvorder Iurida that included the rest of the scorpion families we analyzed in this study. Our results provide a better understanding of the evolution of mitogenome features and phylogenetic relationships in scorpions.

The Complete Mitochondrial Genome of the Chinese White Wax Scale Insect, Ericerus pela Chavannes (Hemiptera: Coccidae), with Novel Gene Arrangement and Truncated tRNA Genes

The Chinese white wax scale insect, Ericerus pela Chavannes (Hemiptera: Coccidae), is one of the scale insects with great economic value and has been dispersed and reared in China for over one thousand years. Its mitochondrial genome provides essential information for the molecular identification and genetic study of this species. We assembled the complete mitochondrial genome of E. pela based on PacBio sequencing and analyzed its genomic features. The genome was 17,766 bp in length with 13 protein-coding genes, 22 tRNAs, and two rRNA genes. The analysis results showed E. pela had significant gene rearrangements involving tRNAs compared with other Coccoidea species. Furthermore, E. pela's nine tRNAs were identified to have obvious truncated structures. The phylogenetic tree compiled of the species showed a long branch of the Coccoidea lineage, which indicated the high evolutionary rate in this group. Our study revealed the mitochondrial characteristics of E. pela and enriched the mitochondrial genetic information on Coccoidea species. It also determined the occurrence of gene rearrangement for the species in this superfamily.

The Complete Mitochondrial Genome of Box Tree Moth Cydalima perspectalis and Insights into Phylogenetics in Pyraloidea

To resolve and reconstruct phylogenetic relationships within Pyraloidea based on molecular data, the mitochondrial genome (mitogenome) was widely applied to understand phylogenetic relations at different taxonomic levels. In this research, a complete mitogenome of Cydalima perspectalis was recorded, and the phylogenetic position of C. perspectalis was inferred based on the sequence in combination with other available sequence data. According to the research, the circular mitochondrial genome is 15,180 bp in length. It contains 22 transfer RNA genes (tRNAs), two ribosomal RNA genes (rRNAs), 13 typical protein-coding genes (PCGs), and a non-coding control region. The arrangement of a gene of the C. perspectalis mitogenome is not the same as the putative ancestral arthropod mitogenome. All of the PCGs are initiated by ATN codons, except for the cytochrome c oxidase subunit 1 (cox1) gene, which is undertaken by CGA. Five genes have incomplete stop codons that contain only ‘T’. All tRNA genes display a typical clover–leaf structure of mitochondrial tRNA, except for trnS1 (AGN). The control region contained an ‘ATAGG(A)’-like motif followed by a poly-T stretch. Based on the mitochondrial data, phylogenetic analysis within Pyraloidea was carried out using Bayesian inference (BI) and maximum likelihood (ML) analyses. Phylogenetic analysis showed that C. perspectalis is more closely related to Pygospila tyres within Spilomelinae than those of Crambidae and Pyraloidea.

MitoFish, MitoAnnotator, and MiFish Pipeline: Updates in 10 Years

MitoFish, MitoAnnotator, and MiFish Pipeline are comprehensive databases of fish mitochondrial genomes (mitogenomes), accurate annotation software of fish mitogenomes, and a web platform for metabarcoding analysis of fish mitochondrial environmental DNA (eDNA), respectively. The MitoFish Suite currently receives over 48,000 visits worldwide every year; however, the performance and usefulness of the online platforms can still be improved. Here, we present essential updates on these platforms, including an enrichment of the reference data sets, an enhanced searching function, substantially faster genome annotation and eDNA analysis with the denoising of sequencing errors, and a multisample comparative analysis function. These updates have made our platform more intuitive, effective, and reliable. These updated platforms are freely available at http://mitofish.aori.u-tokyo.ac.jp/.

The mitochondrial genome of Heterosentis pseudobagri (Wang & Zhang, 1987) Pichelin & Cribb, 1999 reveals novel aspects of tRNA genes evolution in Acanthocephala

BACKGROUND

Acanthocephala is a clade of obligate endoparasites whose mitochondrial genomes (mitogenomes) and evolution remain relatively poorly understood. Previous studies reported that atp8 is lacking from acanthocephalan mitogenomes, and that tRNA genes often have nonstandard structures. Heterosentis pseudobagri (Arhythmacanthidae) is an acanthocephalan fish endoparasite for which no molecular data are currently available, and biological information is unavailable in the English language. Furthermore, there are currently no mitogenomes available for Arhythmacanthidae.

METHODS

We sequenced its mitogenome and transcriptome, and conducted comparative mitogenomic analyses with almost all available acanthocephalan mitogenomes.

RESULTS

The mitogenome had all genes encoded on the same strand and unique gene order in the dataset. Among the 12 protein-coding genes, several genes were highly divergent and annotated with difficulty. Moreover, several tRNA genes could not be identified automatically, so we had to identify them manually via a detailed comparison with orthologues. As common in acanthocephalans, some tRNAs lacked either the TWC arm or the DHU arm, but in several cases, we annotated tRNA genes only on the basis of the conserved narrow central segment comprising the anticodon, while the flanking 5' and 3' ends did not exhibit any resemblance to orthologues and they could not be folded into a tRNA secondary structure. We corroborated that these are not sequencing artefacts by assembling the mitogenome from transcriptomic data. Although this phenomenon was not observed in previous studies, our comparative analyses revealed the existence of highly divergent tRNAs in multiple acanthocephalan lineages.

CONCLUSIONS

These findings indicate either that multiple tRNA genes are non-functional or that (some) tRNA genes in (some) acanthocephalans might undergo extensive posttranscriptional tRNA processing which restores them to more conventional structures. It is necessary to sequence mitogenomes from yet unrepresented lineages and further explore the unusual patterns of tRNA evolution in Acanthocephala.

Pathways to polar adaptation in fishes revealed by long-read sequencing

Long-read sequencing is driving a new reality for genome science in which highly contiguous assemblies can be produced efficiently with modest resources. Genome assemblies from long-read sequences are particularly exciting for understanding the evolution of complex genomic regions that are often difficult to assemble. In this study, we utilized long-read sequencing data to generate a high-quality genome assembly for an Antarctic eelpout, Ophthalmolycus amberensis, the first for the globally distributed family Zoarcidae. We used this assembly to understand how O. amberensis has adapted to the harsh Southern Ocean and compared it to another group of Antarctic fishes: the notothenioids. We showed that selection has largely acted on different targets in eelpouts relative to notothenioids. However, we did find some overlap; in both groups, genes involved in membrane structure, thermal tolerance and vision have evidence of positive selection. We found evidence for historical shifts of transposable element activity in O. amberensis and other polar fishes, perhaps reflecting a response to environmental change. We were specifically interested in the evolution of two complex genomic loci known to underlie key adaptations to polar seas: haemoglobin and antifreeze proteins (AFPs). We observed unique evolution of the haemoglobin MN cluster in eelpouts and related fishes in the suborder Zoarcoidei relative to other Perciformes. For AFPs, we identified the first species in the suborder with no evidence of afpIII sequences (Cebidichthys violaceus) in the genomic region where they are found in all other Zoarcoidei, potentially reflecting a lineage-specific loss of this cluster. Beyond polar fishes, our results highlight the power of long-read sequencing to understand genome evolution.

A novel, conserved and possibly functional motif "WHWGHTW" in mitochondrial transcription across Bilateria

The animal mitogenomes which undergone a reductive evolution has an obvious loss of coding capacity compared to their known closest relatives, but it has not yet been fully investigated why and how the intergenic regions do not encode protein and have no known functions, are stably maintained, replicated, and transmitted by the genome. These relatively small intergenic regions may not be under neutral evolution and they may have functional and/or regulatory roles that have yet to be identified. Here, the distribution pattern, sequence content and location of a novel sequence motif of 'WWWGHTW' were bioinformatically investigated and characterised by constructing a sampling mitogenome dataset of 1889 species from 14 phyla representing the clade of Bilateria. This motif is reverse complementary of the previously described DmTTF binding sequence and found in the nd4L- (X) -trnT gene cluster. This cluster commonly exhibits a strand displacement region and an intergenic region among the bilaterian superphylums, particularly in Ecdysozoa. This motif may be accepted as a substrate providing binding sites for the specific interaction with transcription factors because of (i) its reverse complementarity of previously described DmTTF binding sequence, and (ii) the possession of G and T nucleotides in the fourth and sixth positions, (iii) the bias on T and G nucleotides instead of C and A in the degenerated positions. This suggestion is also supported by the presence of a strand displacement region in the nd4L- (X) -trnT gene cluster, particularly in Ecdysozoa consisting of the most rearranged mitogenomes among the bilaterian superphylums.

On distinguishing between canonical tRNA genes and tRNA gene fragments in prokaryotes

Automated genome annotation is essential for extracting biological information from sequence data. The identification and annotation of tRNA genes is frequently performed by the software package tRNAscan-SE, the output of which is listed for selected genomes in the Genomic tRNA database (GtRNAdb). Here, we highlight a pervasive error in prokaryotic tRNA gene sets on GtRNAdb: the mis-categorization of partial, non-canonical tRNA genes as standard, canonical tRNA genes. Firstly, we demonstrate the issue using the tRNA gene sets of 20 organisms from the archaeal taxon Thermococcaceae. According to GtRNAdb, these organisms collectively deviate from the expected set of tRNA genes in 15 instances, including the listing of eleven putative canonical tRNA genes. However, after detailed manual annotation, only one of these eleven remains; the others are either partial, non-canonical tRNA genes resulting from the integration of genetic elements or CRISPR-Cas activity (seven instances), or attributable to ambiguities in input sequences (three instances). Secondly, we show that similar examples of the mis-categorization of predicted tRNA sequences occur throughout the prokaryotic sections of GtRNAdb. While both canonical and non-canonical prokaryotic tRNA gene sequences identified by tRNAscan-SE are biologically interesting, the challenge of reliably distinguishing between them remains. We recommend employing a combination of (i) screening input sequences for the genetic elements typically associated with non-canonical tRNA genes, and ambiguities, (ii) activating the tRNAscan-SE automated pseudogene detection function, and (iii) scrutinizing predicted tRNA genes with low isotype scores. These measures greatly reduce manual annotation efforts, and lead to improved prokaryotic tRNA gene set predictions.

Complete Mitogenome of Oreolalax Omeimontis Reveals Phylogenetic Status and Novel Gene Arrangement of Archaeobatrachia

Species of the genus Oreolalax displayed crucial morphological characteristics of vertebrates transitioning from aquatic to terrestrial habitats; thus, they can be regarded as a representative vertebrate genus for this landing phenomenon. But the present phylogenetic status of Oreolalax omeimontis has been controversial with morphological and molecular approaches, and specific gene rearrangements were discovered in all six published Oreolalax mitogenomes, which are rarely observed in Archaeobatrachia. Therefore, this study determined the complete mitogenome of O. omeimontis with the aim of identifying its precise phylogenetic position and novel gene arrangement in Archaeobatrachia. Phylogenetic analysis with Bayesian inference and maximum likelihood indicates O. omeimontis is a sister group to O. lichuanensis, which is consistent with previous phylogenetic analysis based on morphological characteristics, but contrasts with other studies using multiple gene fragments. Moreover, although the duplication of trnM occurred in all seven Oreolalax species, the translocation of trnQ and trnM occurred differently in O. omeimontis to the other six, and this unique rearrangement would happen after the speciation of O. omeimontis. In general, this study sheds new light on the phylogenetic relationships and gene rearrangements of Archaeobatrachia.

The complete mitochondrial genome of the roosterfish Nematistius pectoralis Gill 1862: purifying selection in protein coding genes, organization of the control region, and insights into family-level phylogenomic relationships in the recently erected order Carangiformes

The roosterfish Nematistius pectoralis is considered as one of the most magnificent sportfishes worldwide. This study developed the first genomic resource for this trophy-fish that is heavily targeted by the fly-fishing industry. The complete mitochondrial genome of N. pectoralis was assembled using short read sequences and analyzed in detail. The mitochondrial genome of N. pectoralis is 16,537 bp in length and comprises 13 protein-coding genes (PCGs), 2 ribosomal RNA genes (12S and 16S), and 22 transfer RNA genes. A long intergenic space 770 bp in length was assumed to be the D-loop or Control Region (CR). Most of the PCGs and tRNA genes are encoded in the L-strand. All PCGs are under purifying selection and atp8 and nad6 experienced the least selective pressure. All tRNAs exhibit a cloverleaf secondary structure except tRNA-Serine 1 that lacked the D-arm loop. The D-loop of N. pectoralis exhibits three domains commonly described in other fishes; extended terminal associated sequences (ETAS), central, and conserved sequence block (CSB). A ML phylogenetic reconstruction of the newly recognized order Carangiformes based on all 13 mitochondrial PCGs did not support the monophyly of this clade but recognized several families as monophyletic, including Bothidae, Carangidae, Istiophoridae, Latidae, Paralichthyidae, Polynemidae, and Rhombosoleidae. Nematistius pectoralis was sister to a clade composed of Toxotes chatareus (fam. Toxotidae) + Lactarius lactarius (fam. Lactariidae). This genomic resource developed for N. pectoralis will aid in improving our understanding of the population genomics of and strengthen conservation and management strategies in this remarkable trophy-fish.

The genome sequence of the scarce swallowtail, Iphiclides podalirius

The scarce swallowtail, Iphiclides podalirius (Linnaeus, 1758), is a species of butterfly in the family Papilionidae. Here, we present a chromosome-level genome assembly for Iphiclides podalirius as well as gene and transposable element annotations. We investigate how the density of genomic features differs between the 30 Iphiclides podalirius chromosomes. We find that shorter chromosomes have higher heterozygosity at four-fold-degenerate sites and a greater density of transposable elements. While the first result is an expected consequence of differences in recombination rate, the second suggests a counter-intuitive relationship between recombination and transposable element evolution. This high-quality genome assembly, the first for any species in the tribe Leptocircini, will be a valuable resource for population genomics in the genus Iphiclides and comparative genomics more generally.

High-quality, chromosome-scale genome assemblies: comparisons of three Diaphorina citri (Asian citrus psyllid) geographic populations

The Asian citrus psyllid, Diaphorina citri, is the insect vector of the causal agent of huanglongbing (HLB), a devastating bacterial disease of commercial citrus. Presently, few genomic resources exist for D. citri. In this study, we utilized PacBio HiFi and chromatin confirmation contact (Hi-C) sequencing to sequence, assemble, and compare three high-quality, chromosome-scale genome assemblies of D. citri collected from California, Taiwan, and Uruguay. Our assemblies had final sizes of 282.67 Mb (California), 282.89 Mb (Taiwan), and 266.67 Mb (Uruguay) assembled into 13 pseudomolecules—a reduction in assembly size of 41–45% compared with previous assemblies which we validated using flow cytometry. We identified the X chromosome in D. citri and annotated each assembly for repetitive elements, protein-coding genes, transfer RNAs, ribosomal RNAs, piwi-interacting RNA clusters, and endogenous viral elements. Between 19,083 and 20,357 protein-coding genes were predicted. Repetitive DNA accounts for 36.87–38.26% of each assembly. Comparative analyses and mitochondrial haplotype networks suggest that Taiwan and Uruguay D. citri are more closely related, while California D. citri are closely related to Florida D. citri. These high-quality, chromosome-scale assemblies provide new genomic resources to researchers to further D. citri and HLB research.

205 newly assembled mitogenomes provide mixed evidence for rivers as drivers of speciation for Amazonian primates

Mitochondrial DNA remains a cornerstone for molecular ecology, especially for study species from which high‐quality tissue samples cannot be easily obtained. Methods using mitochondrial markers are usually reliant on reference databases, but these are often incomplete. Furthermore, available mitochondrial genomes often lack crucial metadata, such as sampling location, limiting their utility for many analyses. Here, we assembled 205 new mitochondrial genomes for platyrrhine primates, most from the Amazon and with known sampling locations. We present a dated mitogenomic phylogeny based on these samples along with additional published platyrrhine mitogenomes, and use this to assess support for the long‐standing riverine barrier hypothesis (RBH), which proposes that river formation was a major driver of speciation in Amazonian primates. Along the Amazon, Negro, and Madeira rivers, we found mixed support for the RBH. While we identified divergences that coincide with a river barrier, only some occur synchronously and also overlap with the proposed dates of river formation. The most compelling evidence is for the Amazon river potentially driving speciation within bearded saki monkeys (Chiropotes spp.) and within the smallest extant platyrrhines, the marmosets and tamarins. However, we also found that even large rivers do not appear to be barriers for some primates, including howler monkeys (Alouatta spp.), uakaris (Cacajao spp.), sakis (Pithecia spp.), and robust capuchins (Sapajus spp.). Our results support a more nuanced, clade‐specific effect of riverine barriers and suggest that other evolutionary mechanisms, besides the RBH and allopatric speciation, may have played an important role in the diversification of platyrrhines.

The genome sequence of the lesser marbled fritillary, Brenthis ino, and evidence for a segregating neo-Z chromosome

The lesser marbled fritillary, Brenthis ino (Rottemburg, 1775), is a species of Palearctic butterfly. Male Brenthis ino individuals have been reported to have between 12 and 14 pairs of chromosomes, a much-reduced chromosome number than is typical in butterflies. Here, we present a chromosome-level genome assembly for Brenthis ino, as well as gene and transposable element annotations. The assembly is 411.8 Mb in length with a contig N50 of 9.6 Mb and a scaffold N50 of 29.5 Mb. We also show evidence that the male individual from which we generated HiC data was heterozygous for a neo-Z chromosome, consistent with inheriting 14 chromosomes from one parent and 13 from the other. This genome assembly will be a valuable resource for studying chromosome evolution in Lepidoptera, as well as for comparative and population genomics more generally.

PltRNAdb: Plant transfer RNA database

Transfer RNAs (tRNAs) are intermediate-sized non-coding RNAs found in all organisms that help translate messenger RNA into protein. Recently, the number of sequenced plant genomes has increased dramatically. The availability of this extensive data greatly accelerates the study of tRNAs on a large scale. Here, 8,768,261 scaffolds/chromosomes containing 229,093 giga-base pairs representing whole-genome sequences of 256 plant species were analyzed to identify tRNA genes. As a result, 331,242 nuclear, 3,216 chloroplast, and 1,467 mitochondrial tRNA genes were identified. The nuclear tRNA genes include 275,134 tRNAs decoding 20 standard amino acids, 1,325 suppressor tRNAs, 6,273 tRNAs with unknown isotypes, 48,475 predicted pseudogenes, and 37,873 tRNAs with introns. Efforts also extended to the creation of PltRNAdb (https://bioinformatics.um6p.ma/PltRNAdb/index.php), a data source for tRNA genes from 256 plant species. PltRNAdb website allows researchers to search, browse, visualize, BLAST, and download predicted tRNA genes. PltRNAdb will help improve our understanding of plant tRNAs and open the door to discovering the unknown regulatory roles of tRNAs in plant genomes.

Describing novel mitochondrial genomes of Antarctic amphipods

To date, only one mitogenome from an Antarctic amphipod has been published. Here, novel complete mitochondrial genomes (mitogenomes) of two morphospecies are assembled, namely, Charcotia amundseni and Eusirus giganteus. For the latter species, we have assembled two mitogenomes from different genetic clades of this species. The lengths of Eusirus and Charcotia mitogenomes range from 15,534 to 15,619 base pairs and their mitogenomes are composed of 13 protein coding genes, 22 transfer RNAs, 2 ribosomal RNAs, and 1 putative control region CR. Some tRNAs display aberrant structures suggesting that minimalization is also ongoing in amphipod mitogenomes. The novel mitogenomes of the two Antarctic species have features distinguishing them from other amphipod mitogenomes such as a lower AT-richness in the whole mitogenomes and a negative GC- skew in both strands of protein coding genes. The genetically most variable mitochondrial regions of amphipods are nad6 and atp8, while cox1 shows low nucleotide diversity among closely and more distantly related species. In comparison to the pancrustacean mitochondrial ground pattern, E. giganteus shows a translocation of the nad1 gene, while cytb and nad6 genes are translocated in C. amundseni. Phylogenetic analysis based on mitogenomes illustrates that Eusirus and Charcotia cluster together with other species belonging to the same amphipod superfamilies. In the absence of reference nuclear genomes, mitogenomes can be useful to develop markers for studying population genetics or evolutionary relationships at higher taxonomic levels.

Chromosome-level genome assembly of the bioluminescent cardinalfish Siphamia tubifer, an emerging model for symbiosis research

The bioluminescent symbiosis involving the sea urchin cardinalfish Siphamia tubifer and the luminous bacterium Photobacterium mandapamensis is an emerging vertebrate model for the study of microbial symbiosis. However, little genetic data are available for the host, limiting the scope of research that can be implemented with this association. We present a chromosome-level genome assembly for S. tubifer using a combination of PacBio HiFi sequencing and Hi–C technologies. The final assembly was 1.2 Gb distributed on 23 chromosomes and contained 32,365 protein coding genes with a BUSCO score of 99%. A comparison of the S. tubifer genome to that of another nonluminous species of cardinalfish revealed a high degree of synteny, whereas a comparison to a more distant relative in the sister order Gobiiformes revealed the fusion of two chromosomes in the cardinalfish genomes. The complete mitogenome of S. tubifer was also assembled, and an inversion in the vertebrate WANCY tRNA genes as well as heteroplasmy in the length of the control region were discovered. A phylogenetic analysis based on whole the mitochondrial genome indicated that S. tubifer is divergent from the rest of the cardinalfish family, highlighting the potential role of the bioluminescent symbiosis in the initial divergence of Siphamia. This high-quality reference genome will provide novel opportunities for the bioluminescent S. tubifer–P. mandapamensis association to be used as a model for symbiosis research.

First mitochondrial genomes of Chrysopetalidae (Annelida) from shallow-water and deep-sea chemosynthetic environments

Among Annelida, Chrysopetalidae is an ecologically and morphologically diverse group, which includes shallow-water, deep-sea, free-living, and symbiotic species. Here, the four first mitochondrial genomes of this group are presented and described. One of the free-living shallow-water species Chrysopetalum debile (Chrysopetalinae), one of the yet undescribed free-living deep-sea species Boudemos sp., and those of the two deep-sea bivalve endosymbionts Craseoschema thyasiricola and Iheyomytilidicola lauensis (Calamyzinae). An updated phylogeny of Chrysopetalidae is performed, which supports previous phylogenetic hypotheses within Chrysopetalinae and indicates a complex ecological evolution within Calamyzinae. Additionally, analyses of natural selection pressure in the four mitochondrial genomes and additional genes from the two shallow-water species Bhawania goodei and Arichlidon gathofi were performed. Relaxed selection pressure in the mitochondrion of deep-sea and symbiotic species was found, with many sites under selection identified in the COX3 gene of deep-sea species.

Characterization of The Complete Mitochondrial Genome of the Bromeliad Crab Metopaulias depressus (Rathbun, 1896) (Crustacea: Decapoda: Brachyura: Sesarmidae)

Metopaulias depressus is a non-marine crab endemic to Jamaica that dwells in rainforest bromeliads and exhibits elaborate active parental care behavior. Current genomic resources on M. depressus are rare, limiting the understanding of its adaptation to terrestrial life in species that evolved from marine ancestors. This study reports the complete mitochondrial genome of M. depressus assembled using Sanger sequencing. The AT-rich mitochondrial genome of M. depressus is 15,765 bp in length and comprises 13 protein-coding genes (PCGs), 2 ribosomal RNA genes, and 22 transfer RNA genes. A single 691 bp-long intergenic space is assumed to be the control region (CR) or D-loop. A set of selective pressure analyses indicate that the entirety of the PCGs experience purifying selection. Cox1, cox2, nad5, cox3, and atp6 experience strong purifying selection, and atp8 experiences weak purifying selection compared to the rest of the PCGs. The secondary structures of most tRNA genes exhibit a standard ‘cloverleaf’ structure, with the exception of trnS1, which lacks the dihydroxyuridine (DHU) arm but not the loop, the trnH gene, which lacks the thymine pseudouracil cytosine (T) loop but not the arm, and trnM, which exhibits an overly developed T loop. A maximum likelihood phylogenetic analysis based on all PCGs indicated that M. depressus is more closely related to the genera Clistocoeloma, Nanosesarma, and Parasesarma than to Chiromantes, Geosesarma, and Orisarma. This study contributes to deciphering the phylogenetic relationships within the family Sesarmidae and represents a new genomic resource for this iconic crab species.

Complete Mitogenome of Cruznema Tripartitum Confirms Highly Conserved Gene Arrangement within Family Rhabditidae

Mitochondrial genomes have widely been used as molecular markers in understanding the patterns and processes of nematode evolution. The species in genus Cruznema are free-living bacterivores as well as parasites of crickets and mollusks. The complete mitochondrial genome of C. tripartitum was determined through high-throughput sequencing as the first sequenced representative of the genus Cruznema. The genome is comprised of 14,067 bp nucleotides, and includes 12 protein-coding, two rRNA, and 22 tRNA genes. Phylogenetic analyses based on amino acid data support C. tripartitum as a sister to the clade containing Caenorhabditis elegans and Oscheius chongmingensis. The analysis of gene arrangement suggested that C. tripartitum shares the same gene order with O. chongmingensis, Litoditis marina, Diplocapter coronatus, genus Caenorhabditis, and Pristionchus pacificus. Thus, the mitochondrial gene arrangement is highly conserved in the family Rhabditidae as well as some species in Diplogasteridae.

Novel gene rearrangement in the mitochondrial genome of Anastatus fulloi (Hymenoptera Chalcidoidea) and phylogenetic implications for Chalcidoidea

The genus Anastatus comprises a large group of parasitoids, including several biological control agents in agricultural and forest systems. The taxonomy and phylogeny of these species remain controversial. In this study, the mitogenome of A. fulloi Sheng and Wang was sequenced and characterized. The nearly full-length mitogenome of A. fulloi was 15,692 bp, compromising 13 protein-coding genes (PCGs), 2 rRNA genes, 22 tRNA genes and a control region (CR). The total A + T contents were 83.83%, 82.18%, 87.58%, 87.27%, and 82.13% in the whole mitogenome, 13 PCGs, 22 tRNA genes, 2 rRNA genes, and CR, respectively. The mitogenome presented negative AT skews and positive GC skews, except for the CR. Most PCGs were encoded on the heavy strand, started with ATN codons, and ended with TAA codons. Among the 3736 amino acid-encoding codons, TTA (Leu1), CGA (Arg), TCA (Ser2), and TCT (Ser2) were predominant. Most tRNAs had cloverleaf secondary structures, except trnS1, with the absence of a dihydrouridine (DHU) arm. Compared with mitogenomes of the ancestral insect and another parasitoid within Eupelmidae, large-scale rearrangements were found in the mitogenome of A. fulloi, especially inversions and inverse transpositions of tRNA genes. The gene arrangements of parasitoid mitogenomes within Chalcidoidea were variable. A novel gene arrangement was presented in the mitogenome of A. fulloi. Phylogenetic analyses based on the 13 protein-coding genes of 20 parasitoids indicated that the phylogenetic relationship of 6 superfamilies could be presented as Mymaridae + (Eupelmidae + (Encyrtidae + (Trichogrammatidae + (Pteromalidae + Eulophidae)))). This study presents the first mitogenome of the Anastatus genus and offers insights into the identification, taxonomy, and phylogeny of these parasitoids.

The complete mitochondrial genome of an endemic cichlid Etroplus canarensis from Western Ghats, India (Perciformes: Cichlidae) and molecular phylogenetic analysis

BACKGROUND

The Indian endemic cichlid Etroplus canarensis (Canara pearl spot) is an endangered fish and is one among the three Indian cichlids (Etroplinae) that had a restricted distribution in the South Canara region of Karnataka, India. Despite considerable investigations, the phylogeny of Indian Cichlids and its biogeographical origin is still ambiguous and remains a question under discussion which is scrutinized based on whole mitogenomes in the present study.

METHODS AND RESULTS

We report the 16,339 bp complete mitochondrial genome of E. canarensis for the first time using the next-generation sequencing methods. Comparison of gene arrangement and genome characterization was found to commensurate with the previous reports on two Indian cichlid fishes, E. suratensis and E. maculatus. ND6 has been identified as a gene with the highest evolutionary rate and COI and COII is the most conserved gene based on p-genetic distance calculation. Substitution rate (ka/ks) was found to be very low indicating a reduced rate of evolution among subfamily Etroplinae accounting for its subsided species divergence of Indian cichlids. Phylogenetic analysis of Indian cichlids based on a combined dataset of 12 protein-coding genes representing cichlids generated high posterior probability values pillaring paraphyletic nature of Indian Malagasy lineage and monophyletic Indian genus Etroplus.

CONCLUSION

The mitogenome sequence of E. canarensis may provide fundamental molecular data useful for further researches on genetic diversity, endemicity and the conservation of this endangered freshwater fish.