tRNA punctuation model of RNA processing in human mitochondria

Mitogenomic recognition of incognito lineages in the mud spiny lobster Panulirus polyphagus (Herbst, 1793): A tale of unique genetic structuring and diversification

This study provides the first documentation of three deep conspecific lineages within Panulirus polyphagus in the Indian Ocean, bridging the gap in genetic research. Comparative mitogenomics between lineages (L) at both species and family levels, evolutionary relationships and heterogeneity of sequence divergence within Decapoda, and divergence time estimation were performed. The characterized mitogenomes ranged from 15,685-15,705 bp in size and exhibited a typical pancrustacean pattern. Among the three lineages, L1 predominated the Bay of Bengal, L2 the Arabian Sea, and L2.a, a less common lineage genetically closer to L2, was restricted to the latter region. A minor lineage L1.a, was observed in the Coral Triangle area. All PCGs displayed evidence of purifying selection across species and family levels. The largest genetic distance (K2P) between lineages was 9 %, notably between L1.a and L2.a. The phylogenetic tree subdivided the Achelates into Palinuridae and Scyllaridae, and the topology demonstrated a distinct pattern of lineage diversification within P. polyphagus. AliGROOVE analysis revealed no discernible divergence in Decapoda. The diversification of P. polyphagus appears to have occurred during Miocene, with further diversification in Pliocene. Furthermore, genetic stocks and population connectivity recognized here will provide valuable insight for spatial management planning of this dwindling resource.

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

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

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

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

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

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

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

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

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

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

RNA 5-methylcytosine marks mitochondrial double-stranded RNAs for degradation and cytosolic release

Mitochondria are essential regulators of innate immunity. They generate long mitochondrial double-stranded RNAs (mt-dsRNAs) and release them into the cytosol to trigger an immune response under pathological stress conditions. Yet the regulation of these self-immunogenic RNAs remains largely unknown. Here, we employ CRISPR screening on mitochondrial RNA (mtRNA)-binding proteins and identify NOP2/Sun RNA methyltransferase 4 (NSUN4) as a key regulator of mt-dsRNA expression in human cells. We find that NSUN4 induces 5-methylcytosine (m5C) modification on mtRNAs, especially on the termini of light-strand long noncoding RNAs. These m5C-modified RNAs are recognized by complement C1q-binding protein (C1QBP), which recruits polyribonucleotide nucleotidyltransferase to facilitate RNA turnover. Suppression of NSUN4 or C1QBP results in increased mt-dsRNA expression, while C1QBP deficiency also leads to increased cytosolic mt-dsRNAs and subsequent immune activation. Collectively, our study unveils the mechanism underlying the selective degradation of light-strand mtRNAs and establishes a molecular mark for mtRNA decay and cytosolic release.

Context-dependent structure formation of hairpin motifs in bacteriophage MS2 genomic RNA

Many functions of ribonucleic acid (RNA) rely on its ability to assume specific sequence-structure motifs. Packaging signals found in certain RNA viruses are one such prominent example of functional RNA motifs. These signals are short hairpin loops that interact with coat proteins and drive viral self-assembly. As they are found in different positions along the much longer genomic RNA, the formation of their correct structure occurs as a part of a larger context. Any changes to this context can consequently lead to changes in the structure of the motifs themselves. In fact, previous studies have shown that structure and function of RNA motifs can be highly context sensitive to the flanking sequence surrounding them. However, in what ways different flanking sequences influence the structure of an RNA motif they surround has yet to be studied in detail. We focus on a hairpin-rich region of the RNA genome of bacteriophage MS2-a well-studied RNA virus with a wide potential for use in biotechnology-and systematically examine context-dependent structural stability of 14 previously identified hairpin motifs, which include putative and confirmed packaging signals. Combining secondary and tertiary RNA structure prediction of the hairpin motifs placed in different contexts, ranging from the native genomic sequence to random RNA sequences and unstructured poly-U sequences, we determine different measures of motif structural stability. In this way, we show that while some motif structures can be stable in any context, others require specific context provided by the genome. Our results demonstrate the importance of context in RNA structure formation and how changes in the flanking sequence of an RNA motif sometimes lead to drastic changes in its structure. Structural stability of a motif in different contexts could provide additional insights into its functionality as well as assist in determining whether it remains functional when intentionally placed in other contexts.

Implications of p53 in mitochondrial dysfunction and Parkinson's disease

Purpose: To study the underlying molecular mechanisms of p53 in the mitochondrial dysfunction and the pathogenesis of Parkinson's disease (PD), and provide a potential therapeutic target for PD treatment.

Methods: We review the contributions of p53 to mitochondrial changes leading to apoptosis and the subsequent degeneration of dopaminergic neurons in PD.

Results: P53 is a multifunctional protein implicated in the regulation of diverse cellular processes via transcription-dependent and transcription-independent mechanisms. Mitochondria are vital subcellular organelles for that maintain cellular function, and mitochondrial defect and impairment are primary causes of dopaminergic neuron degeneration in PD. Increasing evidence has revealed that mitochondrial dysfunction-associated dopaminergic neuron degeneration is tightly regulated by p53 in PD pathogenesis. Neurodegenerative stress triggers p53 activation, which induces mitochondrial changes, including transmembrane permeability, reactive oxygen species production, Ca2+ overload, electron transport chain defects and other dynamic alterations, and these changes contribute to neurodegeneration and are linked closely with PD occurrence and development. P53 inhibition has been shown to attenuate mitochondrial dysfunction and protect dopaminergic neurons from degeneration under conditions of neurodegenerative stress.

Conclusions: p53 appears to be a potential target for neuroprotective therapy of PD.

The Mitogenome of the Subarctic Octocoral Alcyonium digitatum Reveals a Putative tRNAPro Gene Nested within MutS

We sequenced and analyzed the complete mitogenome of a Norwegian isolate of the octocoral Alcyonium digitatum using the Ion Torrent sequencing technology. The 18,790 bp circular mitochondrial genome was found to harbor the same set of 17 genes, which encode 14 protein subunits, two structural ribosomal RNAs and one tRNA, as reported in other octocorals. In addition, we detected a new tRNAPro-like gene sequence nested within the MutS protein coding region. This putative tRNA gene feature appears to be conserved among the octocorals but has not been reported previously. The A. digitatum mitogenome was also shown to harbor an optional gene (ORFA) that encodes a putative protein of 191 amino acids with unknown function. A mitogenome-based phylogenetic analysis, presented as a maximum likelihood tree, showed that A. digitatum clustered with high statistical confidence with two other Alcyonium species endemic to the Mediterranean Sea and the Southeast Pacific Ocean.

The human mitochondrial mRNA structurome reveals mechanisms of gene expression

The human mitochondrial genome encodes crucial oxidative phosphorylation system proteins, pivotal for aerobic energy transduction. They are translated from nine monocistronic and two bicistronic transcripts whose native structures remain unexplored, posing a gap in understanding mitochondrial gene expression. In this work, we devised the mitochondrial dimethyl sulfate mutational profiling with sequencing (mitoDMS-MaPseq) method and applied detection of RNA folding ensembles using expectation-maximization (DREEM) clustering to unravel the native mitochondrial messenger RNA (mt-mRNA) structurome in wild-type (WT) and leucine-rich pentatricopeptide repeat-containing protein (LRPPRC)-deficient cells. Our findings elucidate LRPPRC's role as a holdase contributing to maintaining mt-mRNA folding and efficient translation. mt-mRNA structural insights in WT mitochondria, coupled with metabolic labeling, unveil potential mRNA-programmed translational pausing and a distinct programmed ribosomal frameshifting mechanism. Our data define a critical layer of mitochondrial gene expression regulation. These mt-mRNA folding maps provide a reference for studying mt-mRNA structures in diverse physiological and pathological contexts.

Complete mitochondrial genome of the introduced Indian walking stick Carausius morosus (Lonchodidae, Insecta) from California

We present the complete mitochondrial genome of Carausius morosus from Salinas, CA. The mitochondrial genome of C. morosus is circular, AT rich (78.1%), and 16,671 bp in length. It consists of 13 protein-coding, 22 transfer RNA, and 2 ribosomal RNA genes and is identical in gene content to Carausius sp.

A novel deep-benthic sea cucumber species of Benthodytes (Holothuroidea, Elasipodida, Psychropotidae) and its comprehensive mitochondrial genome sequencing and evolutionary analysis

BACKGROUND

The holothurians, commonly known as sea cucumbers, are marine organisms that possess significant dietary, nutritional, and medicinal value. However, the National Center for Biotechnology Information (NCBI) currently possesses only approximately 70 complete mitochondrial genome datasets of Holothurioidea, which poses limitations on conducting comprehensive research on their genetic resources and evolutionary patterns. In this study, a novel species of sea cucumber belonging to the genus Benthodytes, was discovered in the western Pacific Ocean. The genomic DNA of the novel sea cucumber was extracted, sequenced, assembled and subjected to thorough analysis.

RESULTS

The mtDNA of Benthodytes sp. Gxx-2023 (GenBank No. OR992091) exhibits a circular structure spanning 17,386 bp, comprising of 13 protein-coding genes (PCGs), 24 non-coding RNAs (2 rRNA genes and 22 tRNA genes), along with two putative control regions measuring 882 bp and 1153 bp, respectively. It exhibits a high AT% content and negative AT-skew, which distinguishing it from the majority of sea cucumbers in terms of environmental adaptability evolution. The mitochondrial gene homology between Gxx-2023 and other sea cucumbers is significantly low, with less than 91% similarity to Benthodytes marianensis, which exhibits the highest level of homology. Additionally, its homology with other sea cucumbers is below 80%. The mitogenome of this species exhibits a unique pattern in terms of start and stop codons, featuring only two types of start codons (ATG and ATT) and three types of stop codons including the incomplete T. Notably, the abundance of AT in the Second position of the codons surpasses that of the First and Third position. The gene arrangement of PCGs exhibits a relatively conserved pattern, while there exists substantial variability in tRNA. Evolutionary analysis revealed that it formed a distinct cluster with B. marianensis and exhibited relatively distant phylogenetic relationships with other sea cucumbers.

CONCLUSIONS

These findings contribute to the taxonomic diversity of sea cucumbers in the Elasipodida order, thereby holding significant implications for the conservation of biological genetic resources, evolutionary advancements, and the exploration of novel sea cucumber resources.

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

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

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

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

Long-read RNA sequencing can probe organelle genome pervasive transcription

40 years ago, organelle genomes were assumed to be streamlined and, perhaps, unexciting remnants of their prokaryotic past. However, the field of organelle genomics has exposed an unparallel diversity in genome architecture (i.e. genome size, structure, and content). The transcription of these eccentric genomes can be just as elaborate - organelle genomes are pervasively transcribed into a plethora of RNA types. However, while organelle protein-coding genes are known to produce polycistronic transcripts that undergo heavy posttranscriptional processing, the nature of organelle noncoding transcriptomes is still poorly resolved. Here, we review how wet-lab experiments and second-generation sequencing data (i.e. short reads) have been useful to determine certain types of organelle RNAs, particularly noncoding RNAs. We then explain how third-generation (long-read) RNA-Seq data represent the new frontier in organelle transcriptomics. We show that public repositories (e.g. NCBI SRA) already contain enough data for inter-phyla comparative studies and argue that organelle biologists can benefit from such data. We discuss the prospects of using publicly available sequencing data for organelle-focused studies and examine the challenges of such an approach. We highlight that the lack of a comprehensive database dedicated to organelle genomics/transcriptomics is a major impediment to the development of a field with implications in basic and applied science.

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

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

Isoform and pathway-specific regulation of post-transcriptional RNA processing in human cells

Steady-state levels of RNA transcripts are controlled by their rates of synthesis and degradation. Here we used nascent RNA Bru-seq and BruChase-seq to profile RNA dynamics across 16 human cell lines as part of ENCODE4 Deeply Profiled Cell Lines collection. We show that RNA turnover dynamics differ widely between transcripts of different genes and between different classes of RNA. Gene set enrichment analysis (GSEA) revealed that transcripts encoding proteins belonging to the same pathway often show similar turnover dynamics. Furthermore, transcript isoforms show distinct dynamics suggesting that RNA turnover is important in regulating mRNA isoform choice. Finally, splicing across newly made transcripts appears to be cooperative with either all or none type splicing. These data sets generated as part of ENCODE4 illustrate the intricate and coordinated regulation of RNA dynamics in controlling gene expression to allow for the precise coordination of cellular functions.

SLIRP promotes autoimmune diseases by amplifying antiviral signaling via positive feedback regulation

The abnormal innate immune response is a prominent feature underlying autoimmune diseases. One emerging factor that can trigger dysregulated immune activation is cytosolic mitochondrial double-stranded RNAs (mt-dsRNAs). However, the mechanism by which mt-dsRNAs stimulate immune responses remains poorly understood. Here, we discover SRA stem-loop interacting RNA binding protein (SLIRP) as a key amplifier of mt-dsRNA-triggered antiviral signals. In autoimmune diseases, SLIRP is commonly upregulated, and targeted knockdown of SLIRP dampens the interferon response. We find that the activation of melanoma differentiation-associated gene 5 (MDA5) by exogenous dsRNAs upregulates SLIRP, which then stabilizes mt-dsRNAs and promotes their cytosolic release to activate MDA5 further, augmenting the interferon response. Furthermore, the downregulation of SLIRP partially rescues the abnormal interferon-stimulated gene expression in autoimmune patients' primary cells and makes cells vulnerable to certain viral infections. Our study unveils SLIRP as a pivotal mediator of interferon response through positive feedback amplification of antiviral signaling.

Mechanisms and pathologies of human mitochondrial DNA replication and deletion formation

Human mitochondria possess a multi-copy circular genome, mitochondrial DNA (mtDNA), that is essential for cellular energy metabolism. The number of copies of mtDNA per cell, and their integrity, are maintained by nuclear-encoded mtDNA replication and repair machineries. Aberrant mtDNA replication and mtDNA breakage are believed to cause deletions within mtDNA. The genomic location and breakpoint sequences of these deletions show similar patterns across various inherited and acquired diseases, and are also observed during normal ageing, suggesting a common mechanism of deletion formation. However, an ongoing debate over the mechanism by which mtDNA replicates has made it difficult to develop clear and testable models for how mtDNA rearrangements arise and propagate at a molecular and cellular level. These deletions may impair energy metabolism if present in a high proportion of the mtDNA copies within the cell, and can be seen in primary mitochondrial diseases, either in sporadic cases or caused by autosomal variants in nuclear-encoded mtDNA maintenance genes. These mitochondrial diseases have diverse genetic causes and multiple modes of inheritance, and show notoriously broad clinical heterogeneity with complex tissue specificities, which further makes establishing genotype-phenotype relationships challenging. In this review, we aim to cover our current understanding of how the human mitochondrial genome is replicated, the mechanisms by which mtDNA replication and repair can lead to mtDNA instability in the form of large-scale rearrangements, how rearranged mtDNAs subsequently accumulate within cells, and the pathological consequences when this occurs.

Comparative Mitogenomics and Phylogenetic Implications for Nine Species of the Subfamily Meconematinae (Orthoptera: Tettigoniidae)

Currently, the subfamily Meconematinae encompasses 1029 species, but whole-mitochondrial-genome assemblies have only been made available for 13. In this study, the whole mitochondrial genomes (mitogenomes) of nine additional species in the subfamily Meconematinae were sequenced. The size ranged from 15,627 bp to 17,461 bp, indicating double-stranded circular structures. The length of the control region was the main cause of the difference in mitochondrial genome length among the nine species. All the mitogenomes including 13 protein-coding genes (PCGs), 22 transfer RNA genes (tRNAs), 2 ribosomal RNA genes (rRNAs) and a control region (CR). The majority strand encoded 23 genes, and the minority strand encoded 14 genes. A phylogenetic analysis reaffirmed the monophyletic status of each subfamily, but the monophysitism of Xizicus, Xiphidiopsis and Phlugiolopsis was not supported.

Structural basis for human mitochondrial tRNA maturation

The human mitochondrial genome is transcribed into two RNAs, containing mRNAs, rRNAs and tRNAs, all dedicated to produce essential proteins of the respiratory chain. The precise excision of tRNAs by the mitochondrial endoribonucleases (mt-RNase), P and Z, releases all RNA species from the two RNA transcripts. The tRNAs then undergo 3'-CCA addition. In metazoan mitochondria, RNase P is a multi-enzyme assembly that comprises the endoribonuclease PRORP and a tRNA methyltransferase subcomplex. The requirement for this tRNA methyltransferase subcomplex for mt-RNase P cleavage activity, as well as the mechanisms of pre-tRNA 3'-cleavage and 3'-CCA addition, are still poorly understood. Here, we report cryo-EM structures that visualise four steps of mitochondrial tRNA maturation: 5' and 3' tRNA-end processing, methylation and 3'-CCA addition, and explain the defined sequential order of the tRNA processing steps. The methyltransferase subcomplex recognises the pre-tRNA in a distinct mode that can support tRNA-end processing and 3'-CCA addition, likely resulting from an evolutionary adaptation of mitochondrial tRNA maturation complexes to the structurally-fragile mitochondrial tRNAs. This subcomplex can also ensure a tRNA-folding quality-control checkpoint before the sequential docking of the maturation enzymes. Altogether, our study provides detailed molecular insight into RNA-transcript processing and tRNA maturation in human mitochondria.

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

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

Malaria parasites require a divergent heme oxygenase for apicoplast gene expression and biogenesis

Malaria parasites have evolved unusual metabolic adaptations that specialize them for growth within heme-rich human erythrocytes. During blood-stage infection, Plasmodium falciparum parasites internalize and digest abundant host hemoglobin within the digestive vacuole. This massive catabolic process generates copious free heme, most of which is biomineralized into inert hemozoin. Parasites also express a divergent heme oxygenase (HO)-like protein (PfHO) that lacks key active-site residues and has lost canonical HO activity. The cellular role of this unusual protein that underpins its retention by parasites has been unknown. To unravel PfHO function, we first determined a 2.8 Å-resolution X-ray structure that revealed a highly α-helical fold indicative of distant HO homology. Localization studies unveiled PfHO targeting to the apicoplast organelle, where it is imported and undergoes N-terminal processing but retains most of the electropositive transit peptide. We observed that conditional knockdown of PfHO was lethal to parasites, which died from defective apicoplast biogenesis and impaired isoprenoid-precursor synthesis. Complementation and molecular-interaction studies revealed an essential role for the electropositive N-terminus of PfHO, which selectively associates with the apicoplast genome and enzymes involved in nucleic acid metabolism and gene expression. PfHO knockdown resulted in a specific deficiency in levels of apicoplast-encoded RNA but not DNA. These studies reveal an essential function for PfHO in apicoplast maintenance and suggest that Plasmodium repurposed the conserved HO scaffold from its canonical heme-degrading function in the ancestral chloroplast to fulfill a critical adaptive role in organelle gene expression.

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

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

Molecular pathways in mitochondrial disorders due to a defective mitochondrial protein synthesis

Mitochondria play a central role in cellular metabolism producing the necessary ATP through oxidative phosphorylation. As a remnant of their prokaryotic past, mitochondria contain their own genome, which encodes 13 subunits of the oxidative phosphorylation system, as well as the tRNAs and rRNAs necessary for their translation in the organelle. Mitochondrial protein synthesis depends on the import of a vast array of nuclear-encoded proteins including the mitochondrial ribosome protein components, translation factors, aminoacyl-tRNA synthetases or assembly factors among others. Cryo-EM studies have improved our understanding of the composition of the mitochondrial ribosome and the factors required for mitochondrial protein synthesis and the advances in next-generation sequencing techniques have allowed for the identification of a growing number of genes involved in mitochondrial pathologies with a defective translation. These disorders are often multisystemic, affecting those tissues with a higher energy demand, and often present with neurodegenerative phenotypes. In this article, we review the known proteins required for mitochondrial translation, the disorders that derive from a defective mitochondrial protein synthesis and the animal models that have been established for their study.

Mitochondrial protein synthesis quality control

Human mitochondrial DNA is one of the most simplified cellular genomes and facilitates compartmentalized gene expression. Within the organelle, there is no physical barrier to separate transcription and translation, nor is there evidence that quality control surveillance pathways are active to prevent translation on faulty mRNA transcripts. Mitochondrial ribosomes synthesize 13 hydrophobic proteins that require co-translational insertion into the inner membrane of the organelle. To maintain the integrity of the inner membrane, which is essential for organelle function, requires responsive quality control mechanisms to recognize aberrations in protein synthesis. In this review, we explore how defects in mitochondrial protein synthesis can arise due to the culmination of inherent mistakes that occur throughout the steps of gene expression. In turn, we examine the stepwise series of quality control processes that are needed to eliminate any mistakes that would perturb organelle homeostasis. We aim to provide an integrated view on the quality control mechanisms of mitochondrial protein synthesis and to identify promising avenues for future research.

Illuminating mitochondrial translation through mouse models

Mitochondria are hubs of metabolic activity with a major role in ATP conversion by oxidative phosphorylation (OXPHOS). The mammalian mitochondrial genome encodes 11 mRNAs encoding 13 OXPHOS proteins along with 2 rRNAs and 22 tRNAs, that facilitate their translation on mitoribosomes. Maintaining the internal production of core OXPHOS subunits requires modulation of the mitochondrial capacity to match the cellular requirements and correct insertion of particularly hydrophobic proteins into the inner mitochondrial membrane. The mitochondrial translation system is essential for energy production and defects result in severe, phenotypically diverse diseases, including mitochondrial diseases that typically affect postmitotic tissues with high metabolic demands. Understanding the complex mechanisms that underlie the pathologies of diseases involving impaired mitochondrial translation is key to tailoring specific treatments and effectively targeting the affected organs. Disease mutations have provided a fundamental, yet limited, understanding of mitochondrial protein synthesis, since effective modification of the mitochondrial genome has proven challenging. However, advances in next generation sequencing, cryoelectron microscopy, and multi-omic technologies have revealed unexpected and unusual features of the mitochondrial protein synthesis machinery in the last decade. Genome editing tools have generated unique models that have accelerated our mechanistic understanding of mitochondrial translation and its physiological importance. Here we review the most recent mouse models of disease pathogenesis caused by defects in mitochondrial protein synthesis and discuss their value for preclinical research and therapeutic development.

RNA degradation in human mitochondria: the journey is not finished

Mitochondria are vital organelles present in almost all eukaryotic cells. Although most of the mitochondrial proteins are nuclear-encoded, mitochondria contain their own genome, whose proper expression is necessary for mitochondrial function. Transcription of the human mitochondrial genome results in the synthesis of long polycistronic transcripts that are subsequently processed by endonucleases to release individual RNA molecules, including precursors of sense protein-encoding mRNA (mt-mRNA) and a vast amount of antisense noncoding RNAs. Because of mitochondrial DNA (mtDNA) organization, the regulation of individual gene expression at the transcriptional level is limited. Although transcription of most protein-coding mitochondrial genes occurs with the same frequency, steady-state levels of mature transcripts are different. Therefore, post-transcriptional processes are important for regulating mt-mRNA levels. The mitochondrial degradosome is a complex composed of the RNA helicase SUV3 (also known as SUPV3L1) and polynucleotide phosphorylase (PNPase, PNPT1). It is the best-characterized RNA-degrading machinery in human mitochondria, which is primarily responsible for the decay of mitochondrial antisense RNA. The mechanism of mitochondrial sense RNA decay is less understood. This review aims to provide a general picture of mitochondrial genome expression, with a particular focus on mitochondrial RNA (mtRNA) degradation.

The molecular machinery for maturation of primary mtDNA transcripts

Human mitochondria harbour a circular, polyploid genome (mtDNA) encoding 11 messenger RNAs (mRNAs), two ribosomal RNAs (rRNAs) and 22 transfer RNAs (tRNAs). Mitochondrial transcription produces long, polycistronic transcripts that span almost the entire length of the genome, and hence contain all three types of RNAs. The primary transcripts then undergo a number of processing and maturation steps, which constitute key regulatory points of mitochondrial gene expression. The first step of mitochondrial RNA processing consists of the separation of primary transcripts into individual, functional RNA molecules and can occur by two distinct pathways. Both are carried out by dedicated molecular machineries that substantially differ from RNA processing enzymes found elsewhere. As a result, the underlying molecular mechanisms remain poorly understood. Over the last years, genetic, biochemical and structural studies have identified key players involved in both RNA processing pathways and provided the first insights into the underlying mechanisms. Here, we review our current understanding of RNA processing in mammalian mitochondria and provide an outlook on open questions in the field.

Mitochondrial RelA empowers mtDNA G-quadruplex formation for hypoxia adaptation in cancer cells

Mitochondrial DNA (mtDNA) G-quadruplexes (G4s) have important regulatory roles in energy metabolism, yet their specific functions and underlying regulatory mechanisms have not been delineated. Using a chemical-genetic screening strategy, we demonstrated that the JAK/STAT3 pathway is the primary regulatory mechanism governing mtDNA G4 dynamics in hypoxic cancer cells. Further proteomic analysis showed that activation of the JAK/STAT3 pathway facilitates the translocation of RelA, a member of the NF-κB family, to the mitochondria, where RelA binds to mtDNA G4s and promotes their folding, resulting in increased mtDNA instability, inhibited mtDNA transcription, and subsequent mitochondrial dysfunction. This binding event disrupts the equilibrium of energy metabolism, catalyzing a metabolic shift favoring glycolysis. Collectively, the results provide insights into a strategy employed by cancer cells to adapt to hypoxia through metabolic reprogramming.

The first mitogenomic phylogenetic framework of Dorcus sensu lato (Coleoptera: Lucanidae), with an emphasis on generic taxonomy in Eastern Asia

BACKGROUND

Dorcus stag beetles in broad sense are one of the most diverse group in Lucanidae and important saproxylic insects playing a crucial role in nutrient recycling and forest biomonitoring. However, the dazzling morphological differentiations have caused numerous systematic confusion within the big genus, especially the puzzlingly generic taxonomy. So far, there is lack of molecular phylogenetic study to address the chaotic situation. In this study, we undertook mitochondrial genome sequencing of 42 representative species including 18 newly-sequenced ones from Eastern Asia and reconstructed the phylogenetic framework of stag beetles in Dorcus sensu lato for the first time.

RESULTS

The mitogenome datasets of Dorcus species have indicated the variable mitogenomic lengths ranged from 15,785 to 19,813 bp. Each mitogenome contained 13 PCGs, 2 rRNAs, 22 tRNAs, and a control region, and all PCGs were under strong purifying selection (Ka/Ks < 1). Notably, we have identified the presence of a substantial intergenic spacer (IGS) between the trnAser (UCN) and NAD1 genes, with varying lengths ranging from 129 bp (in D. hansi) to 158 bp (in D. tityus). The mitogenomic phylogenetic analysis of 42 species showed that Eastern Asia Dorcus was monophyletic, and divided into eight clades with significant genetic distance. Four of them, Clade VIII, VII, VI and I are clustered by the representative species of Serrognathus Motschulsky, Kirchnerius Schenk, Falcicornis Séguy and Dorcus s.s. respectively, which supported their fully generic positions as the previous morphological study presented. The topology also showed the remaining clades were distinctly separated from the species of Dorcus sensu lato, which implied that each of them might demonstrate independent generic status. The Linnaeus nomenclatures were suggested as Eurydorcus Didier stat. res., Eurytrachellelus Didier stat. res., Hemisodorcus Thomson stat. res. and Velutinodorcus Maes stat. res. For Clade V, IV, III and II respectively.

CONCLUSION

This study recognized the monophyly of Dorcus stag beetles and provided a framework for the molecular phylogeny of this group for the first time. The newly generated mitogenomic data serves as a valuable resource for future investigations on lucanid beetles. The generic relationship would facilitate the systematics of Dorcus stag beetles and thus be useful for exploring their evolutionary, ecological, and conservation aspects.

Molecular identification and phylogenetic analysis of the mitogenome in endangered giant nuthatch Sitta magna (Passeriformes, Sittidae)

The Giant Nuthatch Sitta magna (family Sittidae) is a passerine bird, the quantification of the number of habitats and species on a global scale remains low. Most species are restricted to low elevations in southwest China, eastern Myanmar, and northern Thailand. To characterize the mitochondrial genome sequence of S. magna and its phylogenetic relationships with other members within the genus Sitta, the mitochondrial genome of S. magna was sequenced using the whole genome shotgun method. The sequencing results showed that the mitochondrial genome was 16,829 bp long and consisted of 13 protein-coding genes (PCGs), 22 transfer RNA genes (tRNAs), 2 ribosomal RNA genes (rRNAs), and one control region (D-loop). All tRNAs were predicted to form a typical clover secondary structure. Among the 13 PCGs, only the start codon in COI was ATC, the start codon by the remaining 12 PCGs was ATG, and the stop codons were TAG, TAA, AGG, AGA, and TA. Bayesian inference and maximum likelihood phylogenetic analysis of the sequences of 17 species generated consistent well-supported phylogenies. The family Polioptilidae and the family Troglodytidae were closely related, and the family Sittidae was confined to a single branch. The genus Sitta in the family Sittidae was mainly clustered into three branches. Our findings provide new mitochondrial genomic data that could be used for phylogenetic and taxonomic studies; our results also certificate into the phylogenetic relationships within the genus Sitta ((S. himalayensi+(S. nagaensis + S. europaea))+(S. villosa + S. yunnanensis))+(S. carolinensis + S. magna).

Characterization of the Mitogenome of the Genus Dendrocerus Ratzeburg (Hymenoptera: Megaspilidae) with the Specific Designed Primers

In Hymenoptera, the monophyly of Evaniomorpha has been the focus of debate among different scholars. In this study, we sequenced two mitochondrial genomes of Dendrocerus (Hymenoptera: Megaspilidae) to analyze the mitochondrial genomic features of Dendrocerus and provide new molecular data for phylogenetic studies of Evaniomorpha. The mitogenome sizes of D. bellus and D. anisodontus were 15,445 bp and 15,373 bp, respectively, with the trnG of D. bellus missing. The nucleotide composition was significantly biased toward adenine and thymine, with A + T contents of 81.2% (D. bellus) and 82.4% (D. anisodontus). Using Ceraphron sp. (Ceraphronidae) as reference, the Ka/Ks values of NAD4L and NAD6 in D. anisodontus were both greater than one, indicating that non-synonymous mutations are favored by Darwinian selection, which is rare in other hymenopteran species. Compared with Ceraphon sp. gene order, nine operations were identified in D. anisodontus, including four reversals, four TDRLs (tandem duplication random losses) and one transposition, or four reversals and five TDRLs. Phylogenetic analysis of 40 mitochondrial genomes showed that Evaniomorpha was not a monophyletic group, which was also supported by the PBD values. Ceraphronoidea is a monophyletic group and is a sister to Aulacidae + Gasteruptiidae. Based on the conserved region of the newly sequenced mitochondrial genomes, a pair of specific primers MegaF/MegaR was designed for sequencing the COX1 genes in Megaspilidae and a 60% rate of success was achieved in the genus Dendrocerus.

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

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

Role of Rmd9p in 3'-end processing of mitochondrial 15S rRNA in Saccharomyces cerevisiae

Ribosome biogenesis, involving processing/assembly of rRNAs and r-proteins is a vital process. In Saccharomyces cerevisiae mitochondria, ribosomal small subunit comprises 15S rRNA (15S). While the 15S 5'-end processing uses Ccm1p and Pet127p, the mechanisms of the 3'-end processing remain unclear. We reveal involvement of Rmd9p in safeguarding/processing 15S 3'-end. Rmd9p deficiency results in a cleavage at a position 183 nucleotides upstream of 15S 3'-end, and in the loss of the 3'-minor domain. Rmd9p binds to the sequences in the 3'-end region of 15S, and a genetic interaction between rmd9 and dss1 indicates that Rmd9p regulates/limits mtEXO activity during the 3'-end spacer processing.

Phylogenetic placement and comparative analysis of the mitochondrial genomes of Idiostoloidea (Hemiptera: Heteroptera)

The classification system and the higher level phylogenetic relationships of Pentatomomorpha, the second largest infraorder of Heteroptera (Insecta: Hemiptera), have been debated and remain controversial over decades. In particular, the placement and phylogenetic relationship of Idiostoloidea are not well resolved, which hampers a better understanding of the evolutionary history of Pentatomomorpha. In this study, for the first time, we reported the complete mitochondrial genome for two narrowly distributed families of Idiostoloidea (including Idiostolidae and Henicocoridae), respectively. The length of the mitochondrial genome of Monteithocoris hirsutus and Henicocoris sp. is 16,632 and 16,013 bp, respectively. The content of AT is ranging from 75.15% to 80.48%. The mitogenomic structure of Idiostoloidea is highly conservative and there are no gene arrangements. By using the Bayesian inference, maximum likelihood, and Bayesian site-heterogeneous mixture model, we inferred the phylogenetic relationships within Pentatomomorpha and estimated their divergence times based on concatenated mitogenomes and nuclear ribosomal genes. Our results support the classification system of six superfamilies within Pentatomomorpha and confirm the monophyletic groups of each superfamily, with the following phylogenetic relationships: (Aradoidea + (Pentatomoidea + (Idiostoloidea + (Coreoidea + (Pyrrhocoroidea + Lygaeoidea))))). Furthermore, estimated divergence times revealed that most pentatomomorphan superfamilies and families diverged during the Late Jurassic to Early Cretaceous, which coincides with the explosive radiation of angiosperms.

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

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

A kinetic dichotomy between mitochondrial and nuclear gene expression processes

Oxidative phosphorylation (OXPHOS) complexes, encoded by both mitochondrial and nuclear DNA, are essential producers of cellular ATP, but how nuclear and mitochondrial gene expression steps are coordinated to achieve balanced OXPHOS subunit biogenesis remains unresolved. Here, we present a parallel quantitative analysis of the human nuclear and mitochondrial messenger RNA (mt-mRNA) life cycles, including transcript production, processing, ribosome association, and degradation. The kinetic rates of nearly every stage of gene expression differed starkly across compartments. Compared with nuclear mRNAs, mt-mRNAs were produced 1,100-fold more, degraded 7-fold faster, and accumulated to 160-fold higher levels. Quantitative modeling and depletion of mitochondrial factors LRPPRC and FASTKD5 identified critical points of mitochondrial regulatory control, revealing that the mitonuclear expression disparities intrinsically arise from the highly polycistronic nature of human mitochondrial pre-mRNA. We propose that resolving these differences requires a 100-fold slower mitochondrial translation rate, illuminating the mitoribosome as a nexus of mitonuclear co-regulation.

Comparative Analysis of the Complete Mitochondrial Genomes of Three Sisoridae (Osteichthyes, Siluriformes) and the Phylogenetic Relationships of Sisoridae

The family Sisoridae is one of the largest and most diverse Asiatic catfish families, with most species occurring in the water systems of the Qinhai-Tibetan Plateau and East Himalayas. At present, the phylogenetic relationship of the Sisoridae is relatively chaotic. In this study, the mitochondrial genomes (mitogenomes) of three species Creteuchiloglanis kamengensis, Glaridoglanis andersonii, and Exostoma sp. were systematically investigated, the phylogenetic relationships of the family were reconstructed and to determine the phylogenetic position of Exostoma sp. within Sisoridae. The lengths of the mitogenomes' sequences of C. kamengensis, G. andersonii, and Exostoma sp. were 16,589 bp, 16,531 bp, and 16,529 bp, respectively. They all contained one identical control region (D-loop), two ribosomal RNAs (rRNAs), 13 protein-coding genes (PCGs) and 22 transfer RNA (tRNA) genes. We applied two approaches, Bayesian Inference (BI) and Maximum Likelihood (ML), to construct phylogenetic trees. Our findings revealed that the topological structure of both ML and BI trees exhibited significant congruence. Specifically, the phylogenetic tree strongly supports the monophyly of Sisorinae and Glyptosternoids and provides new molecular biological data to support the reconstruction of phylogenetic relationships with Sisoridae. This study is of great scientific value for phylogenetic and genetic variation studies of the Sisoridae.

Long-read direct RNA sequencing of the mitochondrial transcriptome of Saccharomyces cerevisiae reveals condition-dependent intron abundance

Mitochondria fulfil many essential roles and have their own genome, which is expressed as polycistronic transcripts that undergo co- or posttranscriptional processing and splicing. Due to the inherent complexity and limited technical accessibility of the mitochondrial transcriptome, fundamental questions regarding mitochondrial gene expression and splicing remain unresolved, even in the model eukaryote Saccharomyces cerevisiae. Long-read sequencing could address these fundamental questions. Therefore, a method for the enrichment of mitochondrial RNA and sequencing using Nanopore technology was developed, enabling the resolution of splicing of polycistronic genes and the quantification of spliced RNA. This method successfully captured the full mitochondrial transcriptome and resolved RNA splicing patterns with single-base resolution and was applied to explore the transcriptome of S. cerevisiae grown with glucose or ethanol as the sole carbon source, revealing the impact of growth conditions on mitochondrial RNA expression and splicing. This study uncovered a remarkable difference in the turnover of Group II introns between yeast grown in either mostly fermentative or fully respiratory conditions. Whether this accumulation of introns in glucose medium has an impact on mitochondrial functions remains to be explored. Combined with the high tractability of the model yeast S. cerevisiae, the developed method enables to monitor mitochondrial transcriptome responses in a broad range of relevant contexts, including oxidative stress, apoptosis and mitochondrial diseases.

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.

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.

Complete Mitochondrial Genomes of Nedyopus patrioticus: New Insights into the Color Polymorphism of Millipedes

There has been debate about whether individuals with different color phenotypes should have different taxonomic status. In order to determine whether the different color phenotypes of Nedyopus patrioticus require separate taxonomic status or are simply synonyms, here, the complete mitochondrial genomes (mitogenomes) of two different colored N. patrioticus, i.e., red N. patrioticus and white N. patrioticus, are presented. The two mitogenomes were 15,781 bp and 15,798 bp in length, respectively. Each mitogenome contained 13 PCGs, 19 tRNAs, 2 rRNAs, and 1 CR, with a lack of trnI, trnL2, and trnV compared to other Polydesmida species. All genes were located on a single strand in two mitogenomes. Mitochondrial DNA analyses revealed that red N. patrioticus and white N. patrioticus did not show clear evolutionary differences. Furthermore, no significant divergence was discovered by means of base composition analysis. As a result, we suggest that white N. patrioticus might be regarded as a synonym for red N. patrioticus. The current findings confirmed the existence of color polymorphism in N. patrioticus, which provides exciting possibilities for future research. It is necessary to apply a combination of molecular and morphological methods in the taxonomy of millipedes.

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.

The genome and transcriptome of the snail Biomphalaria sudanica s.l.: immune gene diversification and highly polymorphic genomic regions in an important African vector of Schistosoma mansoni

BACKGROUND

Control and elimination of schistosomiasis is an arduous task, with current strategies proving inadequate to break transmission. Exploration of genetic approaches to interrupt Schistosoma mansoni transmission, the causative agent for human intestinal schistosomiasis in sub-Saharan Africa and South America, has led to genomic research of the snail vector hosts of the genus Biomphalaria. Few complete genomic resources exist, with African Biomphalaria species being particularly underrepresented despite this being where the majority of S. mansoni infections occur. Here we generate and annotate the first genome assembly of Biomphalaria sudanica sensu lato, a species responsible for S. mansoni transmission in lake and marsh habitats of the African Rift Valley. Supported by whole-genome diversity data among five inbred lines, we describe orthologs of immune-relevant gene regions in the South American vector B. glabrata and present a bioinformatic pipeline to identify candidate novel pathogen recognition receptors (PRRs).

RESULTS

De novo genome and transcriptome assembly of inbred B. sudanica originating from the shoreline of Lake Victoria (Kisumu, Kenya) resulted in a haploid genome size of ~ 944.2 Mb (6,728 fragments, N50 = 1.067 Mb), comprising 23,598 genes (BUSCO = 93.6% complete). The B. sudanica genome contains orthologues to all described immune genes/regions tied to protection against S. mansoni in B. glabrata, including the polymorphic transmembrane clusters (PTC1 and PTC2), RADres, and other loci. The B. sudanica PTC2 candidate immune genomic region contained many PRR-like genes across a much wider genomic region than has been shown in B. glabrata, as well as a large inversion between species. High levels of intra-species nucleotide diversity were seen in PTC2, as well as in regions linked to PTC1 and RADres orthologues. Immune related and putative PRR gene families were significantly over-represented in the sub-set of B. sudanica genes determined as hyperdiverse, including high extracellular diversity in transmembrane genes, which could be under pathogen-mediated balancing selection. However, no overall expansion in immunity related genes was seen in African compared to South American lineages.

CONCLUSIONS

The B. sudanica genome and analyses presented here will facilitate future research in vector immune defense mechanisms against pathogens. This genomic/transcriptomic resource provides necessary data for the future development of molecular snail vector control/surveillance tools, facilitating schistosome transmission interruption mechanisms in Africa.

The First Complete Mitochondrial Genome of the Genus Litostrophus: Insights into the Rearrangement and Evolution of Mitochondrial Genomes in Diplopoda

This study presents the complete mitochondrial genome (mitogenome) of Litostrophus scaber, which is the first mitogenome of the genus Litostrophus. The mitogenome is a circular molecule with a length of 15,081 bp. The proportion of adenine and thymine (A + T) was 69.25%. The gene ND4L used TGA as the initiation codon, while the other PCGs utilized ATN (A, T, G, C) as the initiation codons. More than half of the PCGs used T as an incomplete termination codon. The transcription direction of the L. scaber mitogenome matched Spirobolus bungii, in contrast to most millipedes. Novel rearrangements were found in the L. scaber mitogenome: trnQ -trnC and trnL1- trnP underwent short-distance translocations and the gene block rrnS-rrnL-ND1 moved to a position between ND4 and ND5, resulting in the formation of a novel gene order. The phylogenetic analysis showed that L. scaber is most closely related to S. bungii, followed by Narceus magnum. These findings enhance our understanding of the rearrangement and evolution of Diplopoda mitogenomes.

Mitogenomes Provide Insights into the Species Boundaries and Phylogenetic Relationships among Three Dolycoris Sloe Bugs (Hemiptera: Pentatomidae) from China

(1) Background: The three sloe bugs, Dolycoris baccarum, Dolycoris indicus, and Dolycoris penicillatus, are found in the Chinese mainland and are morphologically similar. The species boundaries and phylogenetic relationships of the three species remain uncertain; (2) Methods: In this study, we generated multiple mitochondrial genomes (mitogenomes) for each of the three species and conducted comparative mitogenomic analysis, species delimitation, and phylogenetic analysis based on these data; (3) Results: Mitogenomes of the three Dolycoris species are conserved in nucleotide composition, gene arrangement, and codon usage. All protein-coding genes (PCGs) were found to be under purifying selection, and the ND4 evolved at the fastest rate. Most species delimitation analyses based on the COI gene and the concatenated 13 PCGs retrieved three operational taxonomic units (OTUs), which corresponded well with the three Dolycoris species identified based on morphological characters. A clear-cut barcode gap was discovered between the interspecific and intraspecific genetic distances of the three Dolycoris species. Phylogenetic analyses strongly supported the monophyly of Dolycoris, with interspecific relationship inferred as (D. indicus + (D. baccarum + D. penicillatus)); (4) Conclusions: Our study provides the first insight into the species boundaries and phylogenetic relationships of the three Dolycoris species distributed across the Chinese mainland.

Comparative Mitogenome Analyses of Fifteen Ramshorn Snails and Insights into the Phylogeny of Planorbidae (Gastropoda: Hygrophila)

Ramshorn snails from the family Planorbidae are important freshwater snails due to their low trophic level, and some of them act as intermediate hosts for zoonotic trematodes. There are about 250 species from 40 genera of Planorbidae, but only 14 species from 5 genera (Anisus, Biomphalaria, Bulinus, Gyraulus, and Planorbella) have sequenced complete mitochondrial genomes (mitogenomes). In this study, we sequenced and assembled a high-quality mitogenome of a ramshorn snail, Polypylis sp. TS-2018, which represented the first mitogenome of the genus. The mitogenome of Polypylis sp. TS-2018 is 13,749 bp in length, which is shorter than that of most gastropods. It contains 13 protein-coding genes (PCGs), 22 transfer RNA (tRNA) genes, and 2 ribosomal RNA (rRNA). We compared mitogenome characteristics, selection pressure, and gene rearrangement among all of the available mitogenomes of ramshorn snails. We found that the nonsynonymous and synonymous substitution rates (Ka/Ks) of most PCGs indicated purifying and negative selection, except for atp8 of Anisus, Biomphalaria, and Gyraulus, which indicated positive selection. We observed that transpositions and reverse transpositions occurred on 10 tRNAs and rrnS, which resulted in six gene arrangement types. We reconstructed the phylogenetic trees using the sequences of PCGs and rRNAs and strongly supported the monophyly of each genus, as well as three tribes in Planorbidae. Both the gene rearrangement and phylogenetic results suggested that Polypylis had a close relationship with Anisus and Gyraulus, while Bulinus was the sister group to all of the other genera. Our results provide useful data for further investigation of species identification, population genetics, and phylogenetics among ramshorn snails.