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.

Biochemical-free enrichment or depletion of RNA classes in real-time during direct RNA sequencing with RISER.

The heterogeneous composition of cellular transcriptomes poses a major challenge for detecting weakly expressed RNA classes, as they can be obscured by abundant RNAs. Although biochemical protocols can enrich or deplete specified RNAs, they are time-consuming, expensive and can compromise RNA integrity. Here we introduce RISER, a biochemical-free technology for the real-time enrichment or depletion of RNA classes. RISER performs selective rejection of molecules during direct RNA sequencing by identifying RNA classes directly from nanopore signals with deep learning and communicating with the sequencing hardware in real time. By targeting the dominant messenger and mitochondrial RNA classes for depletion, RISER reduces their respective read counts by more than 85%, resulting in an increase in sequencing depth of 47% on average for long non-coding RNAs. We also apply RISER for the depletion of globin mRNA in whole blood, achieving a decrease in globin reads by more than 90% as well as an increase in non-globin reads by 16% on average. Furthermore, using a GPU or a CPU, RISER is faster than GPU-accelerated basecalling and mapping. RISER's modular and retrainable software and intuitive command-line interface allow easy adaptation to other RNA classes. RISER is available at https://github.com/comprna/riser .

The mitochondrial tRNA-derived fragment, mt-tRF-LeuTAA, couples mitochondrial metabolism to insulin secretion.

OBJECTIVE: The contribution of the mitochondrial electron transfer system to insulin secretion involves more than just energy provision. We identified a small RNA fragment (mt-tRF-LeuTAA) derived from the cleavage of a mitochondrially-encoded tRNA that is conserved between mice and humans. The role of mitochondrially-encoded tRNA-derived fragments remains unknown. This study aimed to characterize the impact of mt-tRF-LeuTAA, on mitochondrial metabolism and pancreatic islet functions. METHODS: We used antisense oligonucleotides to reduce mt-tRF-LeuTAA levels in primary rat and human islet cells, as well as in insulin-secreting cell lines. We performed a joint transcriptome and proteome analysis upon mt-tRF-LeuTAA inhibition. Additionally, we employed pull-down assays followed by mass spectrometry to identify direct interactors of the fragment. Finally, we characterized the impact of mt-tRF-LeuTAA silencing on the coupling between mitochondrial metabolism and insulin secretion using high-resolution respirometry and insulin secretion assays. RESULTS: Our study unveils a modulation of mt-tRF-LeuTAA levels in pancreatic islets in different Type 2 diabetes models and in response to changes in nutritional status. The level of the fragment is finely tuned by the mechanistic target of rapamycin complex 1. Located within mitochondria, mt-tRF-LeuTAA interacts with core subunits and assembly factors of respiratory complexes of the electron transfer system. Silencing of mt-tRF-LeuTAA in islet cells limits the inner mitochondrial membrane potential and impairs mitochondrial oxidative phosphorylation, predominantly by affecting the Succinate (via Complex II)-linked electron transfer pathway. Lowering mt-tRF-LeuTAA impairs insulin secretion of rat and human pancreatic ß-cells. CONCLUSIONS: Our findings indicate that mt-tRF-LeuTAA interacts with electron transfer system complexes and is a pivotal regulator of mitochondrial oxidative phosphorylation and its coupling to insulin secretion.

Ancestral allele of DNA polymerase gamma modifies antiviral tolerance.

Mitochondria are critical modulators of antiviral tolerance through the release of mitochondrial RNA and DNA (mtDNA and mtRNA) fragments into the cytoplasm after infection, activating virus sensors and type-I interferon (IFN-I) response1-4. The relevance of these mechanisms for mitochondrial diseases remains understudied. Here we investigated mitochondrial recessive ataxia syndrome (MIRAS), which is caused by a common European founder mutation in DNA polymerase gamma (POLG1)5. Patients homozygous for the MIRAS variant p.W748S show exceptionally variable ages of onset and symptoms5, indicating that unknown modifying factors contribute to disease manifestation. We report that the mtDNA replicase POLG1 has a role in antiviral defence mechanisms to double-stranded DNA and positive-strand RNA virus infections (HSV-1, TBEV and SARS-CoV-2), and its p.W748S variant dampens innate immune responses. Our patient and knock-in mouse data show that p.W748S compromises mtDNA replisome stability, causing mtDNA depletion, aggravated by virus infection. Low mtDNA and mtRNA release into the cytoplasm and a slow IFN response in MIRAS offer viruses an early replicative advantage, leading to an augmented pro-inflammatory response, a subacute loss of GABAergic neurons and liver inflammation and necrosis. A population databank of around 300,000 Finnish individuals6 demonstrates enrichment of immunodeficient traits in carriers of the POLG1 p.W748S mutation. Our evidence suggests that POLG1 defects compromise antiviral tolerance, triggering epilepsy and liver disease. The finding has important implications for the mitochondrial disease spectrum, including epilepsy, ataxia and parkinsonism.

Ribosome Profiling and Mass Spectrometry Reveal Widespread Mitochondrial Translation Defects in a Striatal Cell Model of Huntington Disease.

Huntington disease (HD) is caused by an expanded polyglutamine mutation in huntingtin (mHTT) that promotes prominent atrophy in the striatum and subsequent psychiatric, cognitive deficits, and choreiform movements. Multiple lines of evidence point to an association between HD and aberrant striatal mitochondrial functions; however, the present knowledge about whether (or how) mitochondrial mRNA translation is differentially regulated in HD remains unclear. We found that protein synthesis is diminished in HD mitochondria compared to healthy control striatal cell models. We utilized ribosome profiling (Ribo-Seq) to analyze detailed snapshots of ribosome occupancy of the mitochondrial mRNA transcripts in control and HD striatal cell models. The Ribo-Seq data revealed almost unaltered ribosome occupancy on the nuclear-encoded mitochondrial transcripts involved in oxidative phosphorylation (SDHA, Ndufv1, Timm23, Tomm5, Mrps22) in HD cells. By contrast, ribosome occupancy was dramatically increased for mitochondrially encoded oxidative phosphorylation mRNAs (mt-Nd1, mt-Nd2, mt-Nd4, mt-Nd4l, mt-Nd5, mt-Nd6, mt-Co1, mt-Cytb, and mt-ATP8). We also applied tandem mass tag-based mass spectrometry identification of mitochondrial proteins to derive correlations between ribosome occupancy and actual mature mitochondrial protein products. We found many mitochondrial transcripts with comparable or higher ribosome occupancy, but diminished mitochondrial protein products, in HD. Thus, our study provides the first evidence of a widespread dichotomous effect on ribosome occupancy and protein abundance of mitochondria-related genes in HD.

Localization of RNAs to the mitochondria-mechanisms and functions.

The mammalian mitochondrial proteome comprises over 1000 proteins, with the majority translated from nuclear-encoded messenger RNAs (mRNAs). Mounting evidence suggests many of these mRNAs are localized to the outer mitochondrial membrane (OMM) in a pre- or cotranslational state. Upon reaching the mitochondrial surface, these mRNAs are locally translated to produce proteins that are cotranslationally imported into mitochondria. Here, we summarize various mechanisms cells use to localize RNAs, including transfer RNAs (tRNAs), to the OMM and recent technological advancements in the field to study these processes. While most early studies in the field were carried out in yeast, recent studies reveal RNA localization to the OMM and their regulation in higher organisms. Various factors regulate this localization process, including RNA sequence elements, RNA-binding proteins (RBPs), cytoskeletal motors, and translation machinery. In this review, we also highlight the role of RNA structures and modifications in mitochondrial RNA localization and discuss how these features can alter the binding properties of RNAs. Finally, in addition to RNAs related to mitochondrial function, RNAs involved in other cellular processes can also localize to the OMM, including those implicated in the innate immune response and piRNA biogenesis. As impairment of messenger RNA (mRNA) localization and regulation compromise mitochondrial function, future studies will undoubtedly expand our understanding of how RNAs localize to the OMM and investigate the consequences of their mislocalization in disorders, particularly neurodegenerative diseases, muscular dystrophies, and cancers.

MT-TN mutations lead to progressive mitochondrial encephalopathy and promotes mitophagy.

Mitochondrial encephalopathy is a neurological disorder caused by impaired mitochondrial function and energy production. One of the genetic causes of this condition is the mutation of MT-TN, a gene that encodes the mitochondrial transfer RNA (tRNA) for asparagine. MT-TN mutations affect the stability and structure of the tRNA, resulting in reduced protein synthesis and complex enzymatic deficiency of the mitochondrial respiratory chain. Our patient cohort manifests with epileptic encephalopathy, ataxia, hypotonia, and bilateral basal ganglia calcification, which differs from previously reported cases. MT-TN mutation deficiency leads to decreased basal and maximal oxygen consumption rates, disrupted spare respiratory capacity, declined mitochondrial membrane potential, and impaired ATP production. Moreover, MT-TN mutations promote mitophagy, a process of selective degradation of damaged mitochondria by autophagy. Excessive mitophagy further leads to mitochondrial biogensis as a compensatory mechanism. In this study, we provided evidence of pathogenicity for two MT-TN mutations, m.5688 T > C and m.G5691A, explored the molecular mechanisms, and summarized the clinical manifestations of MT-TN mutations. Our study expanded the genotype and phenotypic spectrum and provided new insight into mt-tRNA (Asn)-associated mitochondrial encephalopathy.

Regulation of nuclear transcription by mitochondrial RNA in endothelial cells.

Chromatin-associated RNAs (caRNAs) form a relatively poorly recognized layer of the epigenome. The caRNAs reported to date are transcribed from the nuclear genome. Here, leveraging a recently developed assay for detection of caRNAs and their genomic association, we report that mitochondrial RNAs (mtRNAs) are attached to the nuclear genome and constitute a subset of caRNA, thus termed mt-caRNA. In four human cell types analyzed, mt-caRNAs preferentially attach to promoter regions. In human endothelial cells (ECs), the level of mt-caRNA-promoter attachment changes in response to environmental stress that mimics diabetes. Suppression of a non-coding mt-caRNA in ECs attenuates stress-induced nascent RNA transcription from the nuclear genome, including that of critical genes regulating cell adhesion, and abolishes stress-induced monocyte adhesion, a hallmark of dysfunctional ECs. Finally, we report increased nuclear localization of multiple mtRNAs in the ECs of human diabetic donors, suggesting many mtRNA translocate to the nucleus in a cell stress and disease-dependent manner. These data nominate mt-caRNAs as messenger molecules responsible for mitochondrial-nuclear communication and connect the immediate product of mitochondrial transcription with the transcriptional regulation of the nuclear genome.

DEFECTIVE KERNEL 56 functions in mitochondrial RNA editing and maize seed development.

Proper seed development is essential for achieving grain production, successful seed germination, and seedling establishment in maize (Zea mays). In the past few decades, pentatricopeptide repeat (PPR) proteins have been proven to play an essential role in regulating the development of maize kernels through posttranscriptional RNA modification of mitochondrial genes. However, the underlying mechanisms remain largely unknown. Here, we characterized a mutant of DEFECTIVE KERNEL 56 (DEK56) with defective kernels that exhibited arrested development of both the embryo and endosperm. Accordingly, we isolated DEK56 through a map-based cloning strategy and found that it encoded an E subgroup PPR protein located in the mitochondria. Dysfunction of DEK56 resulted in altered cytidine (C)-to-uridine (U) editing efficiency at 48 editing sites across 21 mitochondrial transcripts. Notably, the editing efficiency of the maturase-related (matR)-1124 site was substantially reduced or abolished in the dek56 mutant. Furthermore, we found that the splicing efficiency of NADH dehydrogenase subunit 4 (nad4) Introns 1 and 3 was substantially reduced in dek56 kernels, which might be a consequence of the defective MatR function. Through a protein-protein interaction test, we hypothesized that DEK56 carries out its function by recruiting the PPR-DYW protein PPR motif, coiled-coil, and DYW domain-containing protein 1 (PCW1). This interaction is facilitated by Multiple Organellar RNA Editing Factors (ZmMORFs) and Glutamine-Rich Protein 23 (ZmGRP23). Based on these findings, we developed a working model of PPR-mediated mitochondrial processing that plays an essential role in the development of maize kernels. The present study will further broaden our understanding of PPR-mediated seed development and provide a theoretical basis for maize improvement.

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.

MINRbase: a comprehensive database of nuclear- and mitochondrial-ribosomal-RNA-derived fragments (rRFs).

We describe the Mitochondrial and Nuclear rRNA fragment database (MINRbase), a knowledge repository aimed at facilitating the study of ribosomal RNA-derived fragments (rRFs). MINRbase provides interactive access to the profiles of 130 238 expressed rRFs arising from the four human nuclear rRNAs (18S, 5.8S, 28S, 5S), two mitochondrial rRNAs (12S, 16S) or four spacers of 45S pre-rRNA. We compiled these profiles by analyzing 11 632 datasets, including the GEUVADIS and The Cancer Genome Atlas (TCGA) repositories. MINRbase offers a user-friendly interface that lets researchers issue complex queries based on one or more criteria, such as parental rRNA identity, nucleotide sequence, rRF minimum abundance and metadata keywords (e.g. tissue type, disease). A 'summary' page for each rRF provides a granular breakdown of its expression by tissue type, disease, sex, ancestry and other variables; it also allows users to create publication-ready plots at the click of a button. MINRbase has already allowed us to generate support for three novel observations: the internal spacers of 45S are prolific producers of abundant rRFs; many abundant rRFs straddle the known boundaries of rRNAs; rRF production is regimented and depends on 'personal attributes' (sex, ancestry) and 'context' (tissue type, tissue state, disease). MINRbase is available at https://cm.jefferson.edu/MINRbase/.

Pathological mutations promote proteolysis of mitochondrial tRNA-specific 2-thiouridylase 1 (MTU1) via mitochondrial caseinolytic peptidase (CLPP).

MTU1 controls intramitochondrial protein synthesis by catalyzing the 2-thiouridine modification of mitochondrial transfer RNAs (mt-tRNAs). Missense mutations in the MTU1 gene are associated with life-threatening reversible infantile hepatic failure. However, the molecular pathogenesis is not well understood. Here, we investigated 17 mutations associated with this disease, and our results showed that most disease-related mutations are partial loss-of-function mutations, with three mutations being particularly severe. Mutant MTU1 is rapidly degraded by mitochondrial caseinolytic peptidase (CLPP) through a direct interaction with its chaperone protein CLPX. Notably, knockdown of CLPP significantly increased mutant MTU1 protein expression and mt-tRNA 2-thiolation, suggesting that accelerated proteolysis of mutant MTU1 plays a role in disease pathogenesis. In addition, molecular dynamics simulations demonstrated that disease-associated mutations may lead to abnormal intermolecular interactions, thereby impairing MTU1 enzyme activity. Finally, clinical data analysis underscores a significant correlation between patient prognosis and residual 2-thiolation levels, which is partially consistent with the AlphaMissense predictions. These findings provide a comprehensive understanding of MTU1-related diseases, offering prospects for modification-based diagnostics and novel therapeutic strategies centered on targeting CLPP.

Base-resolution quantitative DAMM-seq for mapping RNA methylations in tRNA and mitochondrial polycistronic RNA.

The human AlkB family proteins, such as FTO and ALKBH5, are known to mediate RNA m6A demethylation. However, although ALKBH7 localizes in mitochondria and affects metabolism, the detailed biological function and mechanism have remained unknown for years. We developed Demethylation-Assisted Multiple Methylation sequencing (DAMM-seq) to simultaneously detect N1-methyladenosine (m1A), N3-methylcytidine (m3C), N1-methylguanosine (m1G) and N2,N2-dimethylguanosine (m22G) methylations in both steady-state RNA and nascent RNA, and discovered that human ALKBH7 demethylates m22G and m1A within mt-Ile and mt-Leu1 pre-tRNA regions, respectively, in mitochondrial polycistronic RNA. DAMM-seq quantitatively and sensitively monitors the methylation stoichiometry change at pre-tRNA junctions within nascent mt-RNA, revealing the target region where ALKBH7 regulates RNA processing and local structural switch of polycistronic mt-RNAs. A new RNA demethylase in human cells was characterized through the base-resolution quantification of multiple RNA methylations in nascent mt-RNA, resolving the long-standing question about the functional substrate of ALKBH7.

NIR Light-Activated Mitochondrial RNA Cross-Linking Strategy for H2S Monitoring and Prolonged Colorectal Tumor Imaging.

Molecular diffusion and leakage impede the long-term retention of probes/drugs and may cause potential adverse effects in theranostic fields. Spatiotemporally manipulating the organelle-immobilization behavior of probes/drugs for prolonged tumor retention is indispensable to achieving effective cancer diagnosis and therapy. Herein, we propose a rational strategy that could realize near-infrared light-activated ribonucleic acids (RNAs) cross-linking for prolonged tumor retention and simultaneously endogenous hydrogen sulfide (H2S) monitoring in colorectal tumors. Profiting from efficient singlet oxygen (1O2) generation from Cy796 under 808 nm light irradiation, the 1O2-animated furan moiety in Cy796 could covalently cross-link with cytoplasmic RNAs via a cycloaddition reaction and realize organelle immobilization. Subsequently, specific thiolysis of Cy796 assisted with H2S resulted in homologous product Cy644 with reduced 1O2 generation yields and enhanced absolute fluorescence quantum yields (from 7.42 to 27.70%) with blue-shifted absorption and emission, which avoided the molecular oxidation fluorescence quenching effect mediated by 1O2 and validated fluorescence imaging. Furthermore, studies have demonstrated that our proposed strategy possessed adequate capacity for fluorescence imaging and endogenous H2S detection in HCT116 cells, particularly accumulated at the tumor sites, and retained long-term imaging with excellent biocompatibility. The turn-on fluorescence mode and turn-off 1O2 generation efficiency in our strategy successfully realized a diminished fluorescence cross-talk and oxidation quenching effect. It is adequately envisioned that our proposed strategy for monitoring biomarkers and prolonged tumor retention will contribute tremendous dedication in the clinical, diagnostic, and therapeutic fields.

Assessing the Role of Post-Translational Modifications of Mitochondrial RNA Polymerase.

The mitochondrial proteome is subject to abundant post-translational modifications, including lysine acetylation and phosphorylation of serine, threonine, and tyrosine. The biological function of the majority of these protein modifications is unknown. Proteins required for the transcription and translation of mitochondrial DNA (mtDNA) are subject to modification. This suggests that reversible post-translational modifications may serve as a regulatory mechanism for mitochondrial gene transcription, akin to mechanisms controlling nuclear gene expression. We set out to determine whether acetylation or phosphorylation controls the function of mitochondrial RNA polymerase (POLRMT). Mass spectrometry was used to identify post-translational modifications on POLRMT. We analyzed three POLRMT modification sites (lysine 402, threonine 315, threonine 993) found in distinct structural regions. Amino acid point mutants that mimic the modified and unmodified forms of POLRMT were employed to measure the effect of acetylation or phosphorylation on the promoter binding ability of POLRMT in vitro. We found a slight decrease in binding affinity for the phosphomimic at threonine 315. We did not identify large changes in viability, mtDNA content, or mitochondrial transcript level upon overexpression of POLRMT modification mimics in HeLa cells. Our results suggest minimal biological impact of the POLRMT post-translational modifications studied in our system.

Platelet and mitochondrial RNA is decreased in plasma-derived extracellular vesicles in women with preeclampsia-an exploratory study.

BACKGROUND: Circulating extracellular vesicles (EVs) are increased in preeclampsia (PE) and are associated with severity and progression. We examined in this exploratory cohort study if the mRNAs and long noncoding RNAs (lncRNAs) in plasma-derived EVs were dysregulated in PE compared to normal pregnancy and display different temporal patterns during gestation. METHODS: We isolated EVs from plasma at weeks 22-24 and 36-38 in women with and without PE (n=7 in each group) and performed RNA-seq, focusing on mRNAs and lncRNAs. We validated highly expressed mitochondrial and platelet-derived RNAs discovered from central pathways in 60 women with/without PE. We examined further one of the regulated RNAs, noncoding mitochondrially encoded tRNA alanine (MT-TA), in leukocytes and plasma to investigate its biomarker potential and association with clinical markers of PE. RESULTS: We found abundant levels of platelet-derived and mitochondrial RNAs in EVs. Expression of these RNAs were decreased and lncRNAs increased in EVs from PE compared to without PE. These findings were further validated by qPCR for mitochondrial RNAs MT-TA, MT-ND2, MT-CYB and platelet-derived RNAs PPBP, PF4, CLU in EVs. Decreased expression of mitochondrial tRNA MT-TA in leukocytes at 22-24 weeks was strongly associated with the subsequent development of PE. CONCLUSIONS: Platelet-derived and mitochondrial RNA were highly expressed in plasma EVs and were decreased in EVs isolated from women with PE compared to without PE. LncRNAs were mostly increased in PE. The MT-TA in leukocytes may be a useful biomarker for prediction and/or early detection of PE.

A High-Efficiency Capture-Based NGS Approach for Comprehensive Analysis of Mitochondrial Transcriptome.

The transcription of the mitochondrial genome is pivotal for maintenance of mitochondrial functions, and the deregulated mitochondrial transcriptome contributes to various pathological changes. Despite substantial progress having been achieved in uncovering the transcriptional complexity of the nuclear transcriptome, many unknowns and controversies remain for the mitochondrial transcriptome, partially owing to the lack of a highly efficient mitochondrial RNA (mtRNA) sequencing and analysis approach. Here, we first comprehensively evaluated the influence of essential experimental protocols, including strand-specific library construction, two RNA enrichment strategies, and optimal rRNA depletion, on accurately profiling mitochondrial transcriptome in whole-transcriptome sequencing (WTS) data. Based on these insights, we developed a highly efficient approach specifically suitable for targeted sequencing of whole mitochondrial transcriptome, termed capture-based mtRNA seq (CAP), in which strand-specific library construction and optimal rRNA depletion were applied. Compared with WTS, CAP has a great decrease of required data volume without affecting the sensitivity and accuracy of detection. In addition, CAP also characterized the unannotated mt-tRNA transcripts whose expression levels are below the detection limits of conventional WTS. As a proof-of-concept characterization of mtRNAs, the transcription initiation sites and mtRNA cleavage ratio were accurately identified in CAP data. Moreover, CAP had very reliable performance in plasma and single-cell samples, highlighting its wide application. Altogether, the present study has established a highly efficient pipeline for targeted sequencing of mtRNAs, which may pave the way toward functional annotation of mtRNAs and mtRNA-based diagnostic and therapeutic strategies in various diseases.