Losing the stem-loop structure from metazoan mitochondrial tRNAs and co-evolution of interacting factors

Evolving Escherichia coli to use a tRNA with a non-canonical fold as an adaptor of the genetic code

All known bacterial tRNAs adopt the canonical cloverleaf 2D and L-shaped 3D structures. We aimed to explore whether alternative tRNA structures could be introduced in bacterial translation. To this end, we crafted a vitamin-based genetic system to evolve Escherichia coli toward activity of structurally non-canonical tRNAs. The system reliably couples (escape frequency <10-12) growth with the activities of a novel orthogonal histidine suppressor tRNA (HisTUAC) and of the cognate ARS (HisS) via suppression of a GTA valine codon in the mRNA of an enzyme in thiamine biosynthesis (ThiN). Suppression results in the introduction of an essential histidine and thereby confers thiamine prototrophy. We then replaced HisTUAC in the system with non-canonical suppressor tRNAs and selected for growth. A strain evolved to utilize mini HisT, a tRNA lacking the D-arm, and we identified the responsible mutation in an RNase gene (pnp) involved in tRNA degradation. This indicated that HisS, the ribosome, and EF-Tu accept mini HisT ab initio, which we confirmed genetically and through in vitro translation experiments. Our results reveal a previously unknown flexibility of the bacterial translation machinery for the accepted fold of the adaptor of the genetic code and demonstrate the power of the vitamin-based suppression system.

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.

Characterization of the mitochondrial genomes of the Mexican endemic bats Corynorhinus mexicanus and Corynorhinus leonpaniaguae (Chiroptera: Vespertilionidae)

BACKGROUND

The genus Corynorhinus is composed of four recognized species: C. rafinesquii, C. townsendii, C. mexicanus, and C. leonpaniaguae, the latter two being endemic to Mexico. According to the IUCN, C. mexicanus is considered "Near Threatened", as its populations are dwindling and habitats are affected by anthropogenic disturbance. Corynorhinus leonpaniaguae has not been assigned to an IUCN Red List risk category due to its recent description.

METHODS AND RESULTS

In this study, the mitochondrial genomes of C. mexicanus and C. leonpaniaguae were assembled and characterized in detail. The mitochondrial genomes (mtDNA) of C. mexicanus and C. leonpaniaguae have lengths of 16,470 and 16,581 bp respectively, with a predominant nucleotide usage of adenine (31.670% and 31.729%, respectively) and thymine (26.15% and 26.18%, respectively). The mtDNA of C. mexicanus and C. leonpaniaguae is composed of 37 coding and non-coding elements: 22 transfer RNAs (tRNA), 13 protein-coding genes (PCGs), two ribosomal RNAs and a non-coding region, the control region, which has a length of 933 bp and 1,149 bp, respectively. All tRNAs exhibited a cloverleaf secondary structure, with the exception of trn-Ser1 which showed a deletion of the dihydrouridine arm in the two species. All PCGs are subjected to purifying selection, with atp8 being the gene showing the highest Ka/Ks value.

CONCLUSIONS

These are the first whole mitogenomic resources developed for C. mexicanus and C. leonpaniaguae and enhance our knowledge of the ecology of these species and aid in their conservation.

Molecular Identification of Acrossocheilus jishouensis (Teleostei: Cyprinidae) and Its Complete Mitochondrial Genome

Molecular identification, such as DNA barcoding, is a useful tool that is widely applied in distinguishing species. To identify the cyprinid Acrossocheilus jishouensis, which was previously known to be restricted to only its type locality, we conducted molecular identification of this species based on 23 samples in five localities. Molecular identification based on the mitochondrial COI gene sequence showed that the morphologically similar samples from the five populations were all A. jishouensis, as the mean genetic distances between populations were very small (0.1-1.6%); thus, the distribution of this species was substantially expanded. The whole mitochondrial genome of one sample was also assembled, which was 16,594 bp in length and consisted of 13 protein-coding genes (PCGs), two rRNA genes, 22 tRNA genes, and one control region. All PCGs began with ATG except the COI gene, which started with GTG; seven PCGs used the complete stop codon TAA, while four terminated in T(AA) and two ended with TAG. The overall base composition reflected a higher proportion of A+T than G+C and a positive AT-skew and negative GC-skew pattern except for the opposite in ND6. Phylogenetic relationships inferred using BI and ML methods revealed that both Acrossocheilus and Onychostoma were nonmonophyletic, which indicated that the traditional diagnoses between these two genera need to be assessed further. The results of this study not only expanded the known distribution ranges of A. jishouensis, but also provided a valuable data resource for future molecular and evolutionary studies of Acrossocheilus and other cyprinids in Barbinae.

Complete Mitochondrial Genome of the Eggplant Fruit and Shoot Borer, Leucinodes orbonalis Guenée (Lepidoptera: Crambidae), and Comparison with Other Pyraloid Moths

The eggplant fruit and shoot borer (EFSB) (Leucinodes orbonalis Guenée) is a devastating lepidopteran pest of eggplant (Solanum melongena L.) in the Philippines. Management of an insect pest like the EFSB requires an understanding of its biology, evolution, and adaptations. Genomic resources provide a starting point for understanding EFSB biology, as the resources can be used for phylogenetics and population structure studies. To date, genomic resources are scarce for EFSB; thus, this study generated its complete mitochondrial genome (mitogenome). The circular mitogenome is 15,244 bp-long. It contains 37 genes, namely 13 protein-coding, 22 tRNA, and 2 rRNA genes, and has conserved noncoding regions, motifs, and gene syntenies characteristic of lepidopteran mitogenomes. Some protein-coding genes start and end with non-canonical codons. The tRNA genes exhibit a conserved cloverleaf structure, with the exception in trnS1. Partitioned phylogenetic analysis using 72 pyraloids generated highly supported maximum likelihood and Bayesian inference trees revealing expected basal splits between Crambidae and Pyralidae, and Spilomelinae and Pyraustinae. Spilomelinae was recovered to be paraphyletic, with the EFSB robustly placed before the split of Spilomelinae and Pyraustinae. Overall, the EFSB mitogenome resource will be useful for delineations within Spilomelinae and population structure analysis.

Discovery, structure, mechanisms, and evolution of protein-only RNase P enzymes

The endoribonuclease RNase P is responsible for tRNA 5' maturation in all domains of life. A unique feature of RNase P is the variety of enzyme architectures, ranging from dual- to multi-subunit ribonucleoprotein forms with catalytic RNA subunits to protein-only enzymes, the latter occurring as single- or multi-subunit forms or homo-oligomeric assemblies. The protein-only enzymes evolved twice: a eukaryal protein-only RNase P termed PRORP and a bacterial/archaeal variant termed homolog of Aquifex RNase P (HARP); the latter replaced the RNA-based enzyme in a small group of thermophilic bacteria but otherwise coexists with the ribonucleoprotein enzyme in a few other bacteria as well as in those archaea that also encode a HARP. Here we summarize the history of the discovery of protein-only RNase P enzymes and review the state of knowledge on structure and function of bacterial HARPs and eukaryal PRORPs, including human mitochondrial RNase P as a paradigm of multi-subunit PRORPs. We also describe the phylogenetic distribution and evolution of PRORPs, as well as possible reasons for the spread of PRORPs in the eukaryal tree and for the recruitment of two additional protein subunits to metazoan mitochondrial PRORP. We outline potential applications of PRORPs in plant biotechnology and address diseases associated with mutations in human mitochondrial RNase P genes. Finally, we consider possible causes underlying the displacement of the ancient RNA enzyme by a protein-only enzyme in a small group of bacteria.

The protein-only RNase Ps, endonucleases that cleave pre-tRNA: Biological relevance, molecular architectures, substrate recognition and specificity, and protein interactomes

Protein-only RNase P (PRORP) is an essential enzyme responsible for the 5' maturation of precursor tRNAs (pre-tRNAs). PRORPs are classified into three categories with unique molecular architectures, although all three classes of PRORPs share a mechanism and have similar active sites. Single subunit PRORPs, like those found in plants, have multiple isoforms with different localizations, substrate specificities, and temperature sensitivities. Most recently, Arabidopsis thaliana PRORP2 was shown to interact with TRM1A and B, highlighting a new potential role between these enzymes. Work with At PRORPs led to the development of a ribonuclease that is being used to protect against plant viruses. The mitochondrial RNase P complex, found in metazoans, consists of PRORP, TRMT10C, and SDR5C1, and has also been shown to have substrate specificity, although the cause is unknown. Mutations in mitochondrial tRNA and mitochondrial RNase P have been linked to human disease, highlighting the need to continue understanding this complex. The last class of PRORPs, homologs of Aquifex RNase P (HARPs), is found in thermophilic archaea and bacteria. This most recently discovered type of PRORP forms a large homo-oligomer complex. Although numerous structures of HARPs have been published, it is still unclear how HARPs bind pre-tRNAs and in what ratio. There is also little investigation into the substrate specificity and ideal conditions for HARPs. Moving forward, further work is required to fully characterize each of the three classes of PRORP, the pre-tRNA binding recognition mechanism, the rules of substrate specificity, and how these three distinct classes of PRORP evolved. This article is categorized under: RNA Structure and Dynamics > RNA Structure, Dynamics and Chemistry RNA Structure and Dynamics > Influence of RNA Structure in Biological Systems.

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.

Complete mitochondrial genome of Guigarracailaoensis Wang, Chen & Zheng, 2022 (Cypriniformes, Cyprinidae) and its phylogenetic implications

Guigarracailaoensis is a member of family Cyprinidae, subfamily Labeoninae (Cypriniformes) which was recently discovered in southwestern China. Following its initial description, additional information on this species has remained notably scarce. In the current study, we assemble the complete mitochondrial genome (mitogenome) of G.cailaoensis using the Illumina sequencing platform. The mitogenome is identified as a circular, double-stranded DNA sequence of 16,593 base pairs, encompassing 13 protein-coding genes (PCGs), 22 transfer RNA genes, two ribosomal RNA genes, and a putative control region. Maximum-likelihood and Bayesian-inference approaches were used to construct phylogenetic trees for three datasets: (i) PCG sequences of the complete mitogenome (dataset 1); (ii) PCG sequences of the complete mitogenome combined with nuclear DNA (ncDNA) (Rag1) sequence (dataset 2); and (iii) ncDNA (Rag1) sequences (dataset 3). Phylogenetic analyses position G.cailaoensis as a sister taxon to the lineage consisting of Paraqianlabeolineatus Zhao, Sullivan, Zhang & Peng, 2014 and Pseudogyrinocheilusprochilus Fang, 1933 in dataset 1, and to Pseudogyrinocheilusprochilus in dataset 2, species lacking an oral disc on the lower lip. However, G.cailaoensis showed a close relationship to the lineage consisting of Discogobio and Discocheilus in dataset 3, species possessing an oral disc on the lower lip. Nonetheless, a variety of species with an oral disc on the lower lip are clustered into different lineages across the three datasets that may indicate that the development of the oral disc is homoplastic within the subfamily Labeoninae. The outcomes of this study have the potential to support conservation efforts for this species and to enrich our understanding of genetic resources in the area.

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.

Complete mitochondrial genome of Metathelazia capsulata (Pneumospiruridae) and comparison with other Spiruromorpha species

Metathelazia capsulata (family Pneumospiruridae) is a lungworm parasitizing the bronchi and bronchioles, described in four species of wild carnivores. Very little molecular data are available on this nematode and none on other species of the Pneumospiruridae family. In this work, we describe for the first time the complete mitogenome (mitochondrial genome) of M. capsulata, being the first described of the family Pneumospiruridae. The mitogenome of M. capsulata has 13,659 bp in length, an A + T content of 79.2%. The mitogenome included 12 protein-coding genes (PCGs) (lacking the atp8 gene), 22 tRNA genes, 2 rRNA genes (all the genes are coded by the heavy strand), and an AT-rich region. The PCGs varied in size (232 bp-1645 bp). Only the tRNA-Trp has the standard cloverleaf secondary structure, while the other 21 do not. The AT-rich region, with a 90.5% A + T content and a length of 389 bp, is located between the cox3 and tRNA-Ala genes. Comparison with the mitogenomes of 29 species of Spiruromorpha infraorder, belonging to different families, demonstrates that M. capsulata mitogenome shared the common characteristics of most of them. The phylogeny constructions yielded phylogenies that were in agreement with the obtained previously by using sequences and gene order data of mitogenomes.

Analysis of the Complete Mitochondrial Genome of Pteronura brasiliensis and Lontra canadensis

P. brasiliensis and L. canadensis are two otter species, which successfully occupied semi-aquatic habitats and diverged from other Mustelidae. Herein, the full-length mitochondrial genome sequences were constructed for these two otter species for the first time. Comparative mitochondrial genome, selection pressure, and phylogenetic independent contrasts (PICs) analyses were conducted to determine the structure and evolutionary characteristics of their mitochondrial genomes. Phylogenetic analyses were also conducted to confirm these two otter species' phylogenetic position. The results demonstrated that the mitochondrial genome structure of P. brasiliensis and L. canadensis were consistent across Mustelidae. However, selection pressure analyses demonstrated that the evolutionary rates of mitochondrial genome protein-coding genes (PCGs) ND1, ND4, and ND4L were higher in otters than in terrestrial Mustelidae, whereas the evolutionary rates of ND2, ND6, and COX1 were lower in otters. Additionally, PIC analysis demonstrated that the evolutionary rates of ND2, ND4, and ND4L markedly correlated with a niche type. Phylogenetic analysis showed that P. brasiliensis is situated at the base of the evolutionary tree of otters, and then L. canadensis diverged from it. This study suggests a divergent evolutionary pattern of Mustelidae mitochondrial genome PCGs, prompting the otters' adaptation to semi-aquatic habitats.

Characterization and Phylogenetic Analysis of the Complete Mitochondrial Genome of Aythya marila

Aythya marila is a large diving duck belonging to the family Anatidae. However, the phylogenetic relationship among these Aythya species remains unclear due to the presence of extensive interspecific hybridization events within the Aythya genus. Here, we sequenced and annotated the complete mitochondrial genome of A. marila, which contained 22 tRNAs, 13 protein-coding genes (PCGs), 2 ribosomal RNAs, and 1 D-loop, with a length of 16,617 bp. The sizes of the PCGs ranged from 297 to 1824 bp and were all, except for ND6, located on the heavy chain (H). ATG and TAA were the most common start and termination codons of the 13 PCGs, respectively. The fastest- and slowest-evolving genes were ATP8 and COI, respectively. Codon usage analysis indicated that CUA, AUC, GCC, UUC, CUC, and ACC were the six most frequent codons. The nucleotide diversity values indicated a high level of genetic diversity in A. marila. F_ST analysis suggested a widespread gene exchange between A. baeri and A. nyroca. Moreover, phylogenetic reconstructions using the mitochondrial genomes of all available Anatidae species showed that, in addition to A. marila, four major clades among the Anatidae (Dendrocygninae, Oxyurinae, Anserinae, and Anatinae) were closely related to A. fuligula. Overall, this study provides valuable information on the evolution of A. marila and new insights into the phylogeny of Anatidae.

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.

Analysis of the complete mitogenome of Daphnia galeata from the Han River, South Korea: structure comparison and control region evolution

Daphnia galeata is an important plankton in aquatic ecosystems. As a widely distributed species, D. galeata has been found throughout the Holarctic region. Understanding the genetic diversity and evolution of D. galeata requires the accumulation of genetic information from different locations. Even though the mitochondrial genome (mitogenome) sequence of D. galeata has already been reported, little is known about the evolution of its mitochondrial control region. In this study, D. galeata samples were collected from the Han River on the Korean Peninsula and its partial nd2 gene was sequenced for haplotype network analysis. This analysis showed that four clades of D. galeata were present in the Holarctic region. Moreover, the D. galeata examined in this study belonged to clade D and was specific to South Korea. The mitogenome of D. galeata from the Han River showed similar gene content and structure compared to sequences reported from Japan. Furthermore, the structure of control region of the Han River was similar to those of Japanese clones and differed substantially from European clone. Finally, a phylogenetic analysis based on the amino acid sequences of 13 protein-coding genes (PCGs) indicated that D. galeata from the Han River formed a cluster with clones collected from Lakes Kasumigaura, Shirakaba, and Kizaki in Japan. The differences in control region structure and stem and loop structure reflect the different evolutionary directions of the mitogenomes from Asian and European clones. These findings improve our understanding of the mitogenome structure and genetic diversity of D. galeata.

Characterization and phylogenetic analysis of the complete mitochondrial genome of Mango tilapia (Sarotherodon galilaeus: Cichlidae)

BACKGROUND

Sarotherodon galilaeus (Linné, 1758) is a member of the family Cichlidae, which is considered the most important aquaculture freshwater species endemic to Africa and the Middle East. The genetics and molecular biology of this species are rare. This requires more comprehensive mitochondrial genomes-based phylogenetics to enhance understanding of the relationship and delineate this species.

METHODS AND RESULTS

Here, we assembled the complete mitogenome of S. galilaeus using Illumina high-throughput sequencing technology. The mango tilapia mitogenome was 16,631 bp in length with an AT composition of 53.4% and 46.4% GC content. It encodes 37 genes comprising two ribosomal RNA genes (rRNAs), 22 transfer RNA genes (tRNAs), and 13 protein-coding genes (PCGs) as well as the D-loop known as the control region. The phylogenetic tree was conducted to provide a relationship within the haplotilapiine lineage based on the maximum likelihood method, and the newly sequenced S. galilaeus was clustered with other Sarotherodon species.

CONCLUSION

Our results provide a new perception of the genetic basis of S. galilaeus species for further research on systematics, evolution, population genetics, and molecular ecology.

Anticodon table of the chloroplast genome and identification of putative quadruplet anticodons in chloroplast tRNAs

The chloroplast genome of 5959 species was analyzed to construct the anticodon table of the chloroplast genome. Analysis of the chloroplast transfer ribonucleic acid (tRNA) revealed the presence of a putative quadruplet anticodon containing tRNAs in the chloroplast genome. The tRNAs with putative quadruplet anticodons were UAUG, UGGG, AUAA, GCUA, and GUUA, where the GUUA anticodon putatively encoded tRNA^Asn. The study also revealed the complete absence of tRNA genes containing ACU, CUG, GCG, CUC, CCC, and CGG anticodons in the chloroplast genome from the species studied so far. The chloroplast genome was also found to encode tRNAs encoding N-formylmethionine (fMet), Ile2, selenocysteine, and pyrrolysine. The chloroplast genomes of mycoparasitic and heterotrophic plants have had heavy losses of tRNA genes. Furthermore, the chloroplast genome was also found to encode putative spacer tRNA, tRNA fragments (tRFs), tRNA-derived, stress-induced RNA (tiRNAs), and the group I introns. An evolutionary analysis revealed that chloroplast tRNAs had evolved via multiple common ancestors and the GC% had more influence toward encoding the tRNA number in the chloroplast genome than the genome size.

Chromosome-level genome and the identification of sex chromosomes in Uloborus diversus

The orb web is a remarkable example of animal architecture that is observed in families of spiders that diverged over 200 million years ago. While several genomes exist for araneid orb-weavers, none exist for other orb-weaving families, hampering efforts to investigate the genetic basis of this complex behavior. Here we present a chromosome-level genome assembly for the cribellate orb-weaving spider Uloborus diversus. The assembly reinforces evidence of an ancient arachnid genome duplication and identifies complete open reading frames for every class of spidroin gene, which encode the proteins that are the key structural components of spider silks. We identified the 2 X chromosomes for U. diversus and identify candidate sex-determining loci. This chromosome-level assembly will be a valuable resource for evolutionary research into the origins of orb-weaving, spidroin evolution, chromosomal rearrangement, and chromosomal sex determination in spiders.

Mitochondrial genomes assembled from non-invasive eDNA metagenomic scat samples in the endangered Amur tiger Panthera tigris altaica

The Amur or Siberian tiger Panthera tigris altaica (Temminck, 1844) is currently restricted to a small region of its original geographical range in northwestern Asia and is considered 'endangered' by the IUCN Red List of Threatened Species. This solitary, territorial, and large top predator is in major need of genomic resources to inform conservation management strategies. This study formally tested if complete mitochondrial genomes of P. tigris altaica can be assembled from non-enriched metagenomic libraries generated from scat eDNA samples using the Illumina sequencing platform and open-access bioinformatics pipelines. The mitogenome of P. tigris altaica was assembled and circularized using the pipeline GetOrganelle with a coverage ranging from 322.7x to 17.6x in four different scat eDNA samples. A nearly complete mitochondrial genome (101x) was retrieved from a fifth scat eDNA sample. The complete or nearly complete mitochondrial genomes of P. tigris altaica were AT-rich and composed of 13 protein coding genes (PCGs), 22 transfer RNA genes, two ribosomal RNA genes, and a putative control region. Synteny observed in all assembled mitogenomes was identical to that reported before for P. tigris altaica and other felids. A phylogenomic analysis based on all PCGs demonstrated that the mitochondrial genomes assembled from scat eDNA reliably identify the sequenced samples as belonging to P. tigris and distinguished the same samples from closely and distantly related congeneric species. This study demonstrates that it is viable to retrieve accurate whole and nearly complete mitochondrial genomes of P. tigris altaica (and probably other felids) from scat eDNA samples without library enrichment protocols and using open-access bioinformatics workflows. This new genomic resource represents a new tool to support conservation strategies (bio-prospecting and bio-monitoring) in this iconic cat.

Comparative mitogenomes provide new insights into phylogeny and taxonomy of the subfamily Xenocyprinae (Cypriniformes: Cyprinidae)

Xenocyprinae is a cyprinid subfamily that not only has a discrete geographic distribution but also has a long history dating to the Early Miocene. However, it is controversial whether systematic classification and some species validity of Xenocyprinae exist, as well as its phylogenetic relationships and evolutionary history. In the present study, we first reviewed the description and taxonomic history of Xenocyprinae, and then the complete mitochondrial genome of Distoechodon compressus, an endemic and locally distributed species belonging to Xenocyprinae, was sequenced and annotated. Finally, all the mitogenomes of Xenocyprinae were compared to reconstruct the phylogenetic relationship and estimate the divergence time. The results showed that the mitogenomes are similar in organization and structure with 16618–16630 bp length from 12 mitogenomes of eight species. Phylogenetic analysis confirmed the monology of Xenocyprinae and illustrated three clades within the Xenocyprinae to consist of ambiguous generic classification. Plagiognathops is a valid genus located at the base of the phylogenetic tree. The genus Xenocypris was originally monophyletic, but X. fangi was excluded. Divergence time estimation revealed that the earliest divergence within Xenocyprinae occurred approximately 12.1 Mya when Plagiognathops separated from the primitive Xenocypris. The main two clades (Xenocypris and (Distoechodon + Pseudobrama + X. fangi)) diverged 10.0 Mya. The major divergence of Xenocyprinae species possibly occurred in the Middle to Late Miocene and Late Pliocene, suggesting that speciation and diversifications could be attributed to the Asian monsoon climate. This study clarifies some controversial issues of systematics and provides essential information on the taxonomy and phylogeny of the subfamily Xenocyprinae.

The Evolutionary Dynamics of the Mitochondrial tRNA in the Cichlid Fish Family

Simple Summary

Cichlids are a unique example of fish diversity and species richness which have been explained by sympatric speciation at different freshwater sources in Africa. The mitochondria contribute to cell vitality by providing energy. It contains a circular genome with an established translation system that is spatially independent of the cytosolic counterpart. The current study aimed to investigate the evolutionary dynamics of the mitochondrial tRNA and its role in cichlids’ diversity. The available cichlid mitogenomes in the public database were filtered, in addition to newly sequenced accessions from a specific cichlid group known as the haplotilapiine lineage that is widely distributed in the Egyptian sector of the Nile River. Based on the comparative analysis of mitogenomic data, we identified 22 tRNA genes, in which a single gene was D-armless, while the cloverleaf secondary structure subdivided into stem-loop formations was predicted and used to define the levels of genetic divergence for the remained tRNAs. Peculiarly, in cichlids, the formation known as “T-arm” showed the lowest polymorphism levels among other structures in contrast to other organisms (e.g., scorpions). Comparing the whole family to the specific haplotilapiine lineage showed that the tryptophan tRNA was the most conserved tRNA, with signatures of possible purifying selection.

Abstract

The mitochondrial transfer RNA genes (tRNAs) attract more attention due to their highly dynamic and rapidly evolving nature. The current study aimed to detect and evaluate the dynamics, characteristic patterns, and variations of mitochondrial tRNAs. The study was conducted in two main parts: first, the published mitogenomic sequences of cichlids mt tRNAs have been filtered. Second, the filtered mitochondrial tRNA and additional new mitogenomes representing the most prevalent Egyptian tilapiine were compared and analyzed. Our results revealed that all 22 tRNAs of cichlids folded into a classical cloverleaf secondary structure with four domains, except for trnS^GCU, missing the D domain in all cichlids. When consensus tRNAs were compared, most of the mutations were observed in the trnP at nucleotide levels (substitutions and indels), in contrast to trnL^UAA. From a structural perspective, the anticodon loop and T-loop formations were the most conserved structures among all parts of the tRNA in contrast to the A-stem and D-loop formations. The trnW was the lowest polymorphic unneutral tRNA among all cichlids (both the family and the haplotilapiine lineage), in contrast with the neutral trnD that was extremely polymorphic among and within the haplotilapiine lineage species compared to other cichlids species. From a phylogenetic perspective, the trnC was extremely hypervariable and neutral tRNA in both haplotilapiine lineage and cichlids but was unable to report correct phylogenetic signal for the cichlids. In contrast to trnI and trnY, less variable neutral tRNAs that were able to cluster the haplotilapiine lineage and cichlids species as previously reported. By observing the DNA polymorphism in the coding DNA sequences (CDS), the highest affected amino acid by non-synonymous mutations was isoleucine and was equally mutated to valine and vice versa; no correlation between mutations in CDS and tRNAs was statistically found. The current study provides an insight into the mitochondrial tRNA evolution and its effect on the cichlid diversity and speciation model at the maternal level.

Complete mitochondrial genomes of Boigakraepelini and Hebiuscraspedogaster (Reptilia, Squamata, Colubridae) and their phylogenetic implications

The complete sequence of the mitochondrial genome is a powerful tool for studying phylogenetic relationships and molecular evolution in various species. In this work, the mitogenomes of Boigakraepelini and Hebiuscraspedogaster were sequenced and characterized for the first time. The lengths of the B.kraepelini and H.craspedogaster mitogenomes were 17,124 bp and 17,120 bp, respectively, and both included 13 protein-coding genes, 22 tRNAs, two rRNAs and two control regions. The arrangements of these mitochondrial genes were the same in B.kraepelini and H.craspedogaster. In addition, both genome compositions showed A+T bias (59.03%, 60.93%) and had positive AT skews (0.179, 0.117) and negative GC skews (-0.397, -0.348). The phylogenetic results illustrated a close relationship between B.kraepelini and the genus Lycodon. Moreover, H.craspedogaster was clustered with other Hebius snakes and closely related to other Natricinae species. These results will provide references for further research on the phylogeny of Colubridae.

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.

Analysis of Complete Mitochondrial Genome of Bohadschia argus (Jaeger, 1833) (Aspidochirotida, Holothuriidae)

Bohadschia argu is a kind of sea cucumber with high economic value; it is the only undisputed species in the genus Bohadschia. In this study, the complete mitochondrial genome (mitogenome) of B. argus was acquired through high-throughput sequencing. The mitochondrial genome of B. argus was 15,656 bp in total length and contained a putative control region (CR) and 37 typical genes of animal mitochondrial genomes, including 13 protein-coding genes (PCGs), 2 ribosomal RNA genes (rrnS and rrnL) and 22 transfer RNA genes (tRNA). The sizes of the PCGs ranged from 168 bp to 1833 bp, and all PCGs except nad6 were encoded on the heavy chain (H). Both rrnS and rrnL were also encoded on the H chain. Twenty-two tRNA genes had positive AT skew and GC skew. All tRNAs had a typical cloverleaf secondary structure except for trnI, in which an arm of dihydrouridine was missing. B. argus shared the same gene arrangement order (the echinoderm ground pattern) as other species in Aspidochirotida. Phylogenetic analysis clearly revealed that B. argus belongs as a member of the Holothuriidae, and it is closely related to members of Actinopyga and Holothuria.

First genomic resource for an endangered neotropical mega-herbivore: the complete mitochondrial genome of the forest-dweller (Baird's) tapir (Tapirus bairdii)

Baird's tapir, or the Central American Tapir Tapirus bairdii (family Tapiridae), is one of the largest mammals native to the forests and wetlands of southern North America and Central America, and is categorized as 'endangered' on the 2014 IUCN Red List of Threatened Species. This study reports, for the first time, the complete mitochondrial genome of T. bairdii and examines the phylogenetic position of T. bairdii amongst closely related species in the same family and order to which it belongs using mitochondrial protein-coding genes (PCG's). The circular, double-stranded, A-T rich mitochondrial genome of T. bairdii is 16,697 bp in length consisting of 13 protein-coding genes (PCG's), two ribosomal RNA genes (rrnS (12s ribosomal RNA and rrnL (16s ribosomal RNA)), and 22 transfer RNA (tRNA) genes. A 33 bp long region was identified to be the origin of replication for the light strand (O_L), and a 1,247 bp long control region (CR) contains the origin of replication for the heavy strand (O_H). A majority of the PCG's and tRNA genes are encoded on the positive, or heavy, strand. The gene order in T. baiirdi is identical to that of T. indicus and T. terrestris, the only two other species of extant tapirs with assembled mitochondrial genomes. An analysis of Ka/Ks ratios for all the PCG's show values <1, suggesting that all these PCGs experience strong purifying selection. A maximum-likelihood phylogenetic analysis supports the monophyly of the genus Tapirus and the order Perissodactyla. The complete annotation and analysis of the mitochondrial genome of T. bairdii will contribute to a better understanding of the population genomic diversity and structure of this species, and it will assist in the conservation and protection of its dwindling populations.

The Complete Mitochondrial Genome of Stichopus naso (Aspidochirotida: Stichopodidae: Stichopus) and Its Phylogenetic Position

The mitochondrial genome is widely used to study the molecular evolution of and perform phylogenetic analyses on animals. In this study, the complete mitochondrial genome (mitogenome) of Stichopus naso was sequenced. The mitogenome was 16,239 bp in length and contained 13 protein-coding genes (PCGs), 22 transfer RNA genes (tRNAs), and 2 ribosomal RNA genes (rRNAs). The genome composition showed positive AT-skew (0.023) and negative GC-skew (−0.158). The order of the mitochondrial genes was consistent with those from the Stichopus and Isostichopus species, whereas it was different from those of other species of Aspidochirotida. The phylogenetic analysis, based on the nucleotide sequences of 13 PCGs through the methods of Bayesian inference (BI) and maximum likelihood (ML), indicated that S. naso has close relationships with S. horrens and S. monotuberculatus, and belongs to a member of Stichopodidae. Our study provides a reference mitogenome for further molecular evolution studies and phylogenetic research on sea cucumbers.

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.

The complete mitochondrial genome of the Mexican-endemic cavefish Ophisternon infernale (Synbranchiformes, Synbranchidae): insights on patterns of selection and implications for synbranchiform phylogenetics

Ophisternoninfernale is one of the 200+ troglobitic fish species worldwide, and one of the two cave-dwelling fishes endemic to the karstic aquifer of the Yucatán Peninsula, Mexico. Because of its elusive nature and the relative inaccessibility of its habitat, there is virtually no genetic information on this enigmatic fish. Herein we report the complete mitochondrial genome of O.infernale, which overall exhibits a configuration comparable to that of other synbranchiforms as well as of more distantly related teleosts. The K_A/K_S ratio indicates that most mtDNA PCGs in synbranchiforms have evolved under strong purifying selection, preventing major structural and functional protein changes. The few instances of PCGs under positive selection might be related to adaptation to decreased oxygen availability. Phylogenetic analysis of mtDNA comparative data from synbranchiforms and closely related taxa (including the indostomid Indostomusparadoxus) corroborate the notion that indostomids are more closely related to synbranchiforms than to gasterosteoids, but without rendering the former paraphyletic. Our phylogenetic results also suggest that New World species of Ophisternon might be more closely related to Synbranchus than to the remaining Ophisternon species. This novel phylogenetic hypothesis, however, should be further tested in the context of a comprehensive systematic study of the group.

Methyltransferase METTL8 is required for 3-methylcytosine modification in human mitochondrial tRNAs

A subset of eukaryotic tRNAs is methylated in the anticodon loop, forming 3-methylcytosine (m³C) modifications. In mammals, the number of tRNAs containing m³C modifications has been expanded to include mitochondrial (mt) tRNA-Ser-UGA and mt-tRNA-Thr-UGU. However, whereas the enzymes catalyzing m³C formation in nuclear-encoded tRNAs have been identified, the proteins responsible for m³C modification in mt-tRNAs are unknown. Here, we show that m³C formation in human mt-tRNAs is dependent upon the Methyltransferase-Like 8 (METTL8) enzyme. We find that METTL8 is a mitochondria-associated protein that interacts with mitochondrial seryl-tRNA synthetase, as well as with mt-tRNAs containing m³C. We demonstrate that human cells deficient in METTL8 exhibit loss of m³C modification in mt-tRNAs, but not nuclear-encoded tRNAs. Consistent with the mitochondrial import of METTL8, the formation of m³C in METTL8-deficient cells could be rescued by re-expression of wildtype METTL8, but not by a METTL8 variant lacking the N-terminal mitochondrial localization signal. Notably, we found METTL8-deficiency in human cells causes alterations in the native migration pattern of mt-tRNA-Ser-UGA, suggesting a role for m³C in tRNA folding. Altogether, these findings demonstrate that METTL8 is required for m³C formation in mitochondrial tRNAs and uncover a potential function for m³C modification in mitochondrial tRNA structure.

Mitogenome of Knodus borki (Cypriniformes: Characidae): genomic characterization and phylogenetic analysis

BACKGROUND

The taxonomic status of Knodu in the family Characidae is not yet clear. This study aimed to address this by sequencing and annotating Knodu borki Zarske, 2008.

MATERIALS AND RESULTS

K. borki Zarske, 2008 was sequenced using a Hiseq platform and the complete mitogenome was assembled in SPAdes v3.15.2 and SOAPdenovo2 v.2.01. The mitogenome of K. borki from Guangzhou, the first sequenced species of the genus Knodu, is 16,837 bp in length and contains 13 protein-coding genes (PCGs), two ribosomal (r) RNAs, 22 transfer (t) RNAs, and one D-loop. Among these 37 genes, 28 are encoded by the heavy strand, while nine are encoded by the light strand. Twenty-one of the tRNAs can form typical cloverleaf secondary structures, except tRNA-Ser1, which lacks dihydrouridine arms. All PCGs have the same start codon (ATG), with the exception of COI (GTG). Four PCGs (ND1, ATP8, ND4L, and ND5) have TAA as the stop codon, ND6 has TAG as the stop codon, COI has AGG as the stop codon, and the remaining seven genes have incomplete stop codons of TA-/T-(ND2, COII, COIII, ND3, ND4, and Cyt b as T-, ATP6 as TA-). Phylogenetic analysis showed that K. borki belongs to the family Characidae.

CONCLUSIONS

Our findings demonstrate that K. borki belongs to the family Characidae, due to consistency with the morphological identification. This study provides molecular information for further research on the phylogeny of the genus Knodus and for analyses of the taxonomic status of Characidae.

Complete mitochondrial genome of Manispentadactylapentadactyla (Mammalia: Pholidota), an endemic subspecies of Chinese pangolin: mitogenome characterisation and phylogenetic implications

The Chinese pangolin Manispentadactyla is critically endangered because of over-exploitation and illegal trafficking and includes three subspecies. However, the taxonomic status of the three subspecies of the Chinese pangolin has not been well resolved, which impedes regional conservation and illegal trade traces. In this study, the complete mitogenome sequence of M.p.pentadactyla, an endemic subspecies of the Chinese pangolin in Taiwan, was determined. The complete mitogenome of M.p.pentadactyla is 16,570 base pairs (bp) in length with 13 protein-coding genes (PCG), 23 transfer RNAs (tRNAs), two ribosomal RNAs and a 1164 bp control region. The overall base composition of the genome showed a slight A + T bias (59.9%), positive AT skew (0.1515) and negative GC skew (-0.3406), which is similar to that of other pangolins. All PCGs started with a typical ATN codon and all tRNAs were typical cloverleaf-shaped secondary structures, except for tRNA-Ser^(GCU). Phylogenetic analysis indicated a monophyletic relationship for M.p.pentadactyla and M.p.aurita and was monophyletic for M.p.pentadactyla, but paraphyletic for M.p.aurita. The paraphyly of M.p.aurita resulted from an incomplete lineage sorting. This study enriched the mitogenome database of the Chinese pangolin and the molecular information obtained should be very useful for future research on mitogenome evolution and genetic diversification in M.pentadactyla.

A genomic portrait of Sturnira parvidens: mitochondrial chromosome, repetitive elements, and microsatellite discovery

The yellow-shouldered bat, Sturnira parvidens (Phyllostomidae), is an abundant and widespread species in southern North America and Mesoamerica. Despite its important ecological role, no genomic resources exist for this species. Using low-coverage short Illumina 150 bp pair-end reads sequencing, this study reports the mitochondrial chromosome and nuclear repetitive elements, including microsatellites, in S. parvidens. The mitochondrial genome of S. parvidens is 16,612 bp in length and is comprised of 13 protein-coding genes, 2 ribosomal RNA genes, and 22 transfer RNA genes. Repetitive elements constituted ~67% of the nuclear genome while ~33% of the genome represented single- or low-copy sequences. A moderate proportion of repetitive sequences (31% putative families) could not be assigned to known repeat element families. Considering only annotated repetitive elements, the most ubiquitous repetitive elements belonged to Class I-LINE and Satellite DNA, which were considerably more abundant than Class I-LTR elements and Class II-DNA transposons (TcMar-Mariner and hAT-Charlie). A total of 193 microsatellites were identified.

The complete mitochondrial genome of the 'Zacatuche' Volcano rabbit (Romerolagus diazi), an endemic and endangered species from the Volcanic Belt of Central Mexico

BACKGROUND

The 'Zacatuche', 'Teporingo', or Volcano rabbit (Romerolagus diazi) belongs to the family Leporidae, is an endemic species restricted to the Central part of the Trans-Mexican Volcanic Belt, and is considered 'endangered' by the IUCN Red List of Threatened Species.

METHODS AND RESULTS

This study reports, for the first time, the complete mitochondrial genome of R. diazi and examined the phylogenetic position of R. diazi among other closely related co-familiar species using mitochondrial protein-coding genes (PCGs). The mitogenome of R. diazi was assembled from short Illumina 150 bp pair-end reads with a coverage of 189x. The AT-rich mitochondrial genome of R. diazi is 17,400 bp in length and is comprised of 13 PCGs, two ribosomal RNA genes, and 22 transfer RNA genes. The gene order observed in the mitochondrial genome of R. diazi is identical to that reported for other leporids. Phylogenetic analyses based on PCGs support the basal position of Romerolagus within the Leporidae, at least when compared to the genera Oryctolagus and Lepus. Nonetheless, additional mitochondrial genomes from species belonging to the genera Bunolagus, Sylvilagus, and Pronolagus, among others, are needed before a more robust conclusion about the derived vs basal placement of Romerolagus within the family Leporidae can be reached based on mitochondrial PCGs.

CONCLUSIONS

This is the first genomic resource developed for R. diazi and it represents a tool to improve our understanding about the ecology and evolutionary biology of this iconic and endangered species.

Comparative mitogenomics of Drosophilidae and the evolution of the Zygothrica genus group (Diptera, Drosophilidae)

The Zygothrica genus group of Drosophilidae encompasses more than 437 species and five genera. Although knowledge regarding its diversity has increased, uncertainties about its monophyly and position within Drosophilidae remain. Genomic approaches have been widely used to address different phylogenetic questions and analyses involving the mitogenome have revealed a cost-efficient tool to these studies. Thus, this work aims to characterize mitogenomes of three species of the Zygothrica genus group (from the Hirtodrosophila, Paraliodrosophila and Zygothrica genera), while comparing them with orthologous sequences from other 23 Drosophilidae species and addressing their phylogenetic position. General content concerning gene order and overlap, nucleotide composition, start and stop codon, codon usage and tRNA structures were compared, and phylogenetic trees were constructed under different datasets. The complete mitogenomes characterized for H. subflavohalterata affinis H002 and P. antennta present the PanCrustacea gene order with 22 transfer RNA (tRNA) genes, two ribosomal RNA (rRNA) genes, 13 protein coding genes and an A+T rich region with two T-stretched elements. Some peculiarities such as the almost complete overlap of genes tRNAH/ND4, tRNAF/ND5 and tRNAS2/ND1 are reported for different Drosophilidae species. Non-canonical secondary structures were encountered for tRNAS1 and tRNAY, revealing patterns that apply at different phylogenetic scales. According to the best depiction of the mitogenomes evolutionary history, the three Neotropical species of the Zygothrica genus group encompass a monophyletic lineage sister to Zaprionus, composing with this genus a clade that is sister to the Drosophila subgenus.

The mitochondrial genome of Faughnia haani (Stomatopoda): novel organization of the control region and phylogenetic position of the superfamily Parasquilloidea

BACKGROUND

Stomatopod crustaceans are aggressive marine predators featuring complex compound eyes and powerful raptorial appendages used for "smashing" or "spearing" prey and/or competitors. Among them, parasquilloids (superfamily Parasquilloidea) possess eyes with 2-3 midband rows of hexagonal ommatidia and spearing appendages. Here, we assembled and analyzed the complete mitochondrial genome of the parasquilloid Faughnia haani and explored family- and superfamily-level phylogenetic relationships within the Stomatopoda based on mitochondrial protein coding genes (PCGs).

RESULTS

The mitochondrial genome of F. haani is 16,089 bp in length and encodes 13 protein coding genes (PCGs), 22 transfer RNA genes, 2 ribosomal RNA genes, and a control region that is relatively well organized, containing 2 GA-blocks, 4 poly-T stretches, various [TA(A)]n-blocks, and 2 hairpin structures. This organized control region is likely a synapomorphic characteristic in the Stomatopoda. Comparison of the control region among superfamilies shows that parasquilloid species are more similar to gonodactyloids than to squilloids and lysiosquilloids given the presence of various  poly-T stretches between the hairpin structures and [TA(A)]n-blocks. Synteny is identical to that reported for other stomatopods and corresponds to the Pancrustacea ground pattern. A maximum-likelihood phylogenetic tree based on PCGs revealed that Parasquilloidea is sister to Lysiosquilloidea and Gonodactyloidea and not to Squilloidea, contradicting previous phylogenetic studies.

CONCLUSIONS

The novel phylogenetic position of Parasquilloidea revealed by our study indicates that 'spearing' raptorial appendages are plesiomorphic and that the 'smashing' type is either derived (as reported in previous studies) or apomorphic. Our results raise the possibility that the spearing raptorial claw may have independently evolved twice. The superfamily Parasquilloidea exhibits a closer relationship with other stomatopod superfamilies with a different raptorial claw type and with dissimilar numbers of midband rows of hexagonal ommatidia. Additional studies focusing on the assembly of mitochondrial genomes from species belonging to different genera, families, and superfamilies within the order Stomatopoda are warranted to reach a robust conclusion regarding the evolutionary history of this iconic clade based on mitochondrial PCGs.

tRNAscan-SE 2.0: improved detection and functional classification of transfer RNA genes

tRNAscan-SE has been widely used for transfer RNA (tRNA) gene prediction for over twenty years, developed just as the first genomes were decoded. With the massive increase in quantity and phylogenetic diversity of genomes, the accurate detection and functional prediction of tRNAs has become more challenging. Utilizing a vastly larger training set, we created nearly one hundred specialized isotype- and clade-specific models, greatly improving tRNAscan-SE’s ability to identify and classify both typical and atypical tRNAs. We employ a new comparative multi-model strategy where predicted tRNAs are scored against a full set of isotype-specific covariance models, allowing functional prediction based on both the anticodon and the highest-scoring isotype model. Comparative model scoring has also enhanced the program's ability to detect tRNA-derived SINEs and other likely pseudogenes. For the first time, tRNAscan-SE also includes fast and highly accurate detection of mitochondrial tRNAs using newly developed models. Overall, tRNA detection sensitivity and specificity is improved for all isotypes, particularly those utilizing specialized models for selenocysteine and the three subtypes of tRNA genes encoding a CAU anticodon. These enhancements will provide researchers with more accurate and detailed tRNA annotation for a wider variety of tRNAs, and may direct attention to tRNAs with novel traits.

Rapid shifts in mitochondrial tRNA import in a plant lineage with extensive mitochondrial tRNA gene loss

In most eukaryotes, transfer RNAs (tRNAs) are one of the very few classes of genes remaining in the mitochondrial genome, but some mitochondria have lost these vestiges of their prokaryotic ancestry. Sequencing of mitogenomes from the flowering plant genus Silene previously revealed a large range in tRNA gene content, suggesting rapid and ongoing gene loss/replacement. Here, we use this system to test longstanding hypotheses about how mitochondrial tRNA genes are replaced by importing nuclear-encoded tRNAs. We traced the evolutionary history of these gene loss events by sequencing mitochondrial genomes from key outgroups (Agrostemma githago and Silene [=Lychnis] chalcedonica). We then performed the first global sequencing of purified plant mitochondrial tRNA populations to characterize the expression of mitochondrial-encoded tRNAs and the identity of imported nuclear-encoded tRNAs. We also confirmed the utility of high-throughput sequencing methods for the detection of tRNA import by sequencing mitochondrial tRNA populations in a species (Solanum tuberosum) with known tRNA trafficking patterns. Mitochondrial tRNA sequencing in Silene revealed substantial shifts in the abundance of some nuclear-encoded tRNAs in conjunction with their recent history of mt-tRNA gene loss and surprising cases where tRNAs with anticodons still encoded in the mitochondrial genome also appeared to be imported. These data suggest that nuclear-encoded counterparts are likely replacing mitochondrial tRNAs even in systems with recent mitochondrial tRNA gene loss, and the redundant import of a nuclear-encoded tRNA may provide a mechanism for functional replacement between translation systems separated by billions of years of evolutionary divergence.

The complete mitochondrial genome of the longneck croaker, pseudotolithus typus Bleeker, 1863 from Sierra Leone

The complete mitochondrial DNA information of Pseudotolithus typus Bleeker, 1863, collected from Sierra Leone was determined using next-generation sequencing (NGS) and bioinfromatic analysis. Its mitogenome (16,504 bp) encoded the typical 13 protein-coding genes (PCGs), 2 ribosomal RNAs (12S & 16S), and 22 tRNAs. All 13 PCGs showed a standard start codon (ATG) but an unusual stop codon (AGA) was identified in COX1 gene. Except for ND6, all 12 PCGs were encoded on the light strand. Except for tRNA^Ser-GCT, 21 tRNAs formed the typical clover-leaf structures. Phylogenetic analysis showed three mitochondrial genomes in the genus Pseudotolithus formed a clade distinct from the other species in the same family. The mitogenome of P. typus identified in this study exhibited 96.27% and 88.86% identity to T. typus in the Guinean water and P. elongatus, respectively. Additional mitogenome sequences of Pseudotolithus species will provide useful information for their scientific management in western African countries.

Comparative mitochondrial genomics of sponge-dwelling snapping shrimps in the genus Synalpheus: Exploring differences between eusocial and non-eusocial species and insights into phylogenetic relationships in caridean shrimps

The genus Synalpheus is a cosmopolitan clade of marine shrimps found in most tropical regions. Species in this genus exhibit a range of social organizations, including pair-forming, communal breeding, and eusociality, the latter only known to have evolved within this genus in the marine realm. This study examines the complete mitochondrial genomes of seven species of Synalpheus and explores differences between eusocial and non-eusocial species considering that eusociality has been shown before to affect the strength of purifying selection in mitochondrial protein coding genes. The AT-rich mitochondrial genomes of Synalpheus range from 15,421 bp to 15,782 bp in length and comprise, invariably, 13 protein-coding genes (PCGs), two ribosomal RNA genes, and 22 transfer RNA genes. A 648 bp to 994 bp long intergenic space is assumed to be the D-loop. Mitochondrial gene synteny is identical among the studied shrimps. No major differences occur between eusocial and non-eusocial species in nucleotide composition and codon usage profiles of PCGs and in the secondary structure of tRNA genes. Maximum likelihood phylogenetic analysis of the complete concatenated PCG complement of 90 species supports the monophyly of the genus Synalpheus and its family Alpheidae. Moreover, the monophyletic status of the caridean families Alvinocaridae, Atyidae, Thoridae, Lysmatidae, Palaemonidae, and Pandalidae within caridean shrimps are fully or highly supported by the analysis. We therefore conclude that mitochondrial genomes contain sufficient phylogenetic information to resolve relationships at high taxonomic levels within the Caridea. Our analysis of mitochondrial genomes in the genus Synalpheus contributes to the understanding of the coevolution between genomic architecture and sociality in caridean shrimps and other marine organisms.

Hopeful monsters: unintended sequencing of famously malformed mite mitochondrial tRNAs reveals widespread expression and processing of sense-antisense pairs

Although tRNA structure is one of the most conserved and recognizable shapes in molecular biology, aberrant tRNAs are frequently found in the mitochondrial genomes of metazoans. The extremely degenerate structures of several mitochondrial tRNAs (mt-tRNAs) have led to doubts about their expression and function. Mites from the arachnid superorder Acariformes are predicted to have some of the shortest mt-tRNAs, with a complete loss of cloverleaf-like shape. While performing mitochondrial isolations and recently developed tRNA-seq methods in plant tissue, we inadvertently sequenced the mt-tRNAs from a common plant pest, the acariform mite Tetranychus urticae, to a high enough coverage to detect all previously annotated T. urticae tRNA regions. The results not only confirm expression, CCA-tailing and post-transcriptional base modification of these highly divergent tRNAs, but also revealed paired sense and antisense expression of multiple T. urticae mt-tRNAs. Mirrored expression of mt-tRNA genes has been hypothesized but not previously demonstrated to be common in any system. We discuss the functional roles that these divergent tRNAs could have as both decoding molecules in translation and processing signals in transcript maturation pathways, as well as how sense–antisense pairs add another dimension to the bizarre tRNA biology of mitochondrial genomes.

Post-Transcriptional Modifications of Conserved Nucleotides in the T-Loop of tRNA: A Tale of Functional Convergent Evolution

The high conservation of nucleotides of the T-loop, including their chemical identity, are hallmarks of tRNAs from organisms belonging to the three Domains of Life. These structural characteristics allow the T-loop to adopt a peculiar intraloop conformation able to interact specifically with other conserved residues of the D-loop, which ultimately folds the mature tRNA in a unique functional canonical L-shaped architecture. Paradoxically, despite the high conservation of modified nucleotides in the T-loop, enzymes catalyzing their formation depend mostly on the considered organism, attesting for an independent but convergent evolution of the post-transcriptional modification processes. The driving force behind this is the preservation of a native conformation of the tRNA elbow that underlies the various interactions of tRNA molecules with different cellular components.

Building a tRNA thermometer to estimate microbial adaptation to temperature

Because ambient temperature affects biochemical reactions, organisms living in extreme temperature conditions adapt protein composition and structure to maintain biochemical functions. While it is not feasible to experimentally determine optimal growth temperature (OGT) for every known microbial species, organisms adapted to different temperatures have measurable differences in DNA, RNA and protein composition that allow OGT prediction from genome sequence alone. In this study, we built a 'tRNA thermometer' model using tRNA sequence to predict OGT. We used sequences from 100 archaea and 683 bacteria species as input to train two Convolutional Neural Network models. The first pairs individual tRNA sequences from different species to predict which comes from a more thermophilic organism, with accuracy ranging from 0.538 to 0.992. The second uses the complete set of tRNAs in a species to predict optimal growth temperature, achieving a maximum ${r^2}$ of 0.86; comparable with other prediction accuracies in the literature despite a significant reduction in the quantity of input data. This model improves on previous OGT prediction models by providing a model with minimum input data requirements, removing laborious feature extraction and data preprocessing steps and widening the scope of valid downstream analyses.

Adaptation of the Romanomermis culicivorax CCA-Adding Enzyme to Miniaturized Armless tRNA Substrates

The mitochondrial genome of the nematode Romanomermis culicivorax encodes for miniaturized hairpin-like tRNA molecules that lack D- as well as T-arms, strongly deviating from the consensus cloverleaf. The single tRNA nucleotidyltransferase of this organism is fully active on armless tRNAs, while the human counterpart is not able to add a complete CCA-end. Transplanting single regions of the Romanomermis enzyme into the human counterpart, we identified a beta-turn element of the catalytic core that-when inserted into the human enzyme-confers full CCA-adding activity on armless tRNAs. This region, originally identified to position the 3'-end of the tRNA primer in the catalytic core, dramatically increases the enzyme's substrate affinity. While conventional tRNA substrates bind to the enzyme by interactions with the T-arm, this is not possible in the case of armless tRNAs, and the strong contribution of the beta-turn compensates for an otherwise too weak interaction required for the addition of a complete CCA-terminus. This compensation demonstrates the remarkable evolutionary plasticity of the catalytic core elements of this enzyme to adapt to unconventional tRNA substrates.

The complete mitochondrial genome of the red-jointed brackish-water fiddler crab Minuca minax (LeConte 1855) (Brachyura: Ocypodidae): New family gene order, and purifying selection and phylogenetic informativeness of protein coding genes

Minuca minax is a semi-terrestrial crustacean that commonly lives in low salinity, riverine habitats along the shores of the eastern United States. This study reports, for the first time, the complete mitochondrial genome of M. minax. The AT-rich mitochondrial genome of M. minax is 15,937 bp in length and comprised of 13 protein-coding genes (PCGs), 2 ribosomal RNA genes, and 22 transfer RNA genes. A single 737 bp long intergenic space is assumed to be the D-loop. Most of the PCGs and tRNA genes are encoded in the L-strand. The gene order observed in the mitochondrial genome of M. minax is new although almost identical to that reported in confamiliar species. In all other confamiliar species to which M. minax is compared, the positions of the trnQ gene and the trnI gene are switched. K_A/K_S ratios calculated for all mitochondrial PCGs show values of <1, indicating that these PCGs are evolving under purifying selection. A maximum likelihood phylogenetic analysis (concatenated PCGs [n = 13], 15 species) supports the monophyly of the subfamilies Ocypodinae and Gelaminidae. Mitochondrial PCGs have enough phylogenetic information to reveal relationships supporting higher taxonomic levels within this family. The knowledge of a complete mitochondrial genome from the red-jointed brackish-water fiddler crab M. minax contributes to the better understanding of meta-population connectivity and the mechanisms involved in the adaptation of marine organisms to near-limnic conditions.

Can Mitogenomes of the Northern Wheatear (Oenanthe oenanthe) Reconstruct Its Phylogeography and Reveal the Origin of Migrant Birds?

The Northern Wheatear (Oenanthe oenanthe, including the nominate and the two subspecies O. o. leucorhoa and O. o. libanotica) and the Seebohm’s Wheatear (Oenanthe seebohmi) are today regarded as two distinct species. Before, all four taxa were regarded as four subspecies of the Northern Wheatear. Their classification has exclusively been based on ecological and morphological traits, while their molecular characterization is still missing. With this study, we used next-generation sequencing to assemble 117 complete mitochondrial genomes covering O. o. oenanthe, O. o. leucorhoa and O. seebohmi. We compared the resolution power of each individual mitochondrial marker and concatenated marker sets to reconstruct the phylogeny and estimate speciation times of three taxa. Moreover, we tried to identify the origin of migratory wheatears caught on Helgoland (Germany) and on Crete (Greece). Mitogenome analysis revealed two different ancient lineages that separated around 400,000 years ago. Both lineages consisted of a mix of subspecies and species. The phylogenetic trees, as well as haplotype networks are incongruent with the present morphology-based classification. Mitogenome could not distinguish these presumed species. The genetic panmixia among present populations and taxa might be the consequence of mitochondrial introgression between ancient wheatear populations.