Epigenetic inheritance of diet-induced and sperm-borne mitochondrial RNAs.

Spermatozoa harbour a complex and environment-sensitive pool of small non-coding RNAs (sncRNAs)1, which influences offspring development and adult phenotypes1-7. Whether spermatozoa in the epididymis are directly susceptible to environmental cues is not fully understood8. Here we used two distinct paradigms of preconception acute high-fat diet to dissect epididymal versus testicular contributions to the sperm sncRNA pool and offspring health. We show that epididymal spermatozoa, but not developing germ cells, are sensitive to the environment and identify mitochondrial tRNAs (mt-tRNAs) and their fragments (mt-tsRNAs) as sperm-borne factors. In humans, mt-tsRNAs in spermatozoa correlate with body mass index, and paternal overweight at conception doubles offspring obesity risk and compromises metabolic health. Sperm sncRNA sequencing of mice mutant for genes involved in mitochondrial function, and metabolic phenotyping of their wild-type offspring, suggest that the upregulation of mt-tsRNAs is downstream of mitochondrial dysfunction. Single-embryo transcriptomics of genetically hybrid two-cell embryos demonstrated sperm-to-oocyte transfer of mt-tRNAs at fertilization and suggested their involvement in the control of early-embryo transcription. Our study supports the importance of paternal health at conception for offspring metabolism, shows that mt-tRNAs are diet-induced and sperm-borne and demonstrates, in a physiological setting, father-to-offspring transfer of sperm mitochondrial RNAs at fertilization.

Cryo-EM Structures Reveal Transcription Initiation Steps by Yeast Mitochondrial RNA Polymerase.

Mitochondrial RNA polymerase (mtRNAP) is crucial in cellular energy production, yet understanding of mitochondrial DNA transcription initiation lags that of bacterial and nuclear DNA transcription. We report structures of two transcription initiation intermediate states of yeast mtRNAP that explain promoter melting, template alignment, DNA scrunching, abortive synthesis, and transition into elongation. In the partially melted initiation complex (PmIC), transcription factor MTF1 makes base-specific interactions with flipped non-template (NT) nucleotides "AAGT" at -4 to -1 positions of the DNA promoter. In the initiation complex (IC), the template in the expanded 7-mer bubble positions the RNA and NTP analog UTPαS, while NT scrunches into an NT loop. The scrunched NT loop is stabilized by the centrally positioned MTF1 C-tail. The IC and PmIC states coexist in solution, revealing a dynamic equilibrium between two functional states. Frequent scrunching/unscruching transitions and the imminent steric clashes of the inflating NT loop and growing RNA:DNA with the C-tail explain abortive synthesis and transition into elongation.

Mitochondrial RNA modifications shape metabolic plasticity in metastasis.

Aggressive and metastatic cancers show enhanced metabolic plasticity1, but the precise underlying mechanisms of this remain unclear. Here we show how two NOP2/Sun RNA methyltransferase 3 (NSUN3)-dependent RNA modifications-5-methylcytosine (m5C) and its derivative 5-formylcytosine (f5C) (refs.2-4)-drive the translation of mitochondrial mRNA to power metastasis. Translation of mitochondrially encoded subunits of the oxidative phosphorylation complex depends on the formation of m5C at position 34 in mitochondrial tRNAMet. m5C-deficient human oral cancer cells exhibit increased levels of glycolysis and changes in their mitochondrial function that do not affect cell viability or primary tumour growth in vivo; however, metabolic plasticity is severely impaired as mitochondrial m5C-deficient tumours do not metastasize efficiently. We discovered that CD36-dependent non-dividing, metastasis-initiating tumour cells require mitochondrial m5C to activate invasion and dissemination. Moreover, a mitochondria-driven gene signature in patients with head and neck cancer is predictive for metastasis and disease progression. Finally, we confirm that this metabolic switch that allows the metastasis of tumour cells can be pharmacologically targeted through the inhibition of mitochondrial mRNA translation in vivo. Together, our results reveal that site-specific mitochondrial RNA modifications could be therapeutic targets to combat metastasis.

Dinucleotide Degradation by REXO2 Maintains Promoter Specificity in Mammalian Mitochondria.

Oligoribonucleases are conserved enzymes that degrade short RNA molecules of up to 5 nt in length and are assumed to constitute the final stage of RNA turnover. Here we demonstrate that REXO2 is a specialized dinucleotide-degrading enzyme that shows no preference between RNA and DNA dinucleotide substrates. A heart- and skeletal-muscle-specific knockout mouse displays elevated dinucleotide levels and alterations in gene expression patterns indicative of aberrant dinucleotide-primed transcription initiation. We find that dinucleotides act as potent stimulators of mitochondrial transcription initiation in vitro. Our data demonstrate that increased levels of dinucleotides can be used to initiate transcription, leading to an increase in transcription levels from both mitochondrial promoters and other, nonspecific sequence elements in mitochondrial DNA. Efficient RNA turnover by REXO2 is thus required to maintain promoter specificity and proper regulation of transcription in mammalian mitochondria.

MISTERMINATE Mechanistically Links Mitochondrial Dysfunction with Proteostasis Failure.

Mitochondrial dysfunction and proteostasis failure frequently coexist as hallmarks of neurodegenerative disease. How these pathologies are related is not well understood. Here, we describe a phenomenon termed MISTERMINATE (mitochondrial-stress-induced translational termination impairment and protein carboxyl terminal extension), which mechanistically links mitochondrial dysfunction with proteostasis failure. We show that mitochondrial dysfunction impairs translational termination of nuclear-encoded mitochondrial mRNAs, including complex-I 30kD subunit (C-I30) mRNA, occurring on the mitochondrial surface in Drosophila and mammalian cells. Ribosomes stalled at the normal stop codon continue to add to the C terminus of C-I30 certain amino acids non-coded by mRNA template. C-terminally extended C-I30 is toxic when assembled into C-I and forms aggregates in the cytosol. Enhancing co-translational quality control prevents C-I30 C-terminal extension and rescues mitochondrial and neuromuscular degeneration in a Parkinson's disease model. These findings emphasize the importance of efficient translation termination and reveal unexpected link between mitochondrial health and proteome homeostasis mediated by MISTERMINATE.

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.

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.

Cytosolic RNA binding of the mitochondrial TCA cycle enzyme malate dehydrogenase.

Several enzymes of intermediary metabolism have been identified to bind RNA in cells, with potential consequences for the bound RNAs and/or the enzyme. In this study, we investigate the RNA-binding activity of the mitochondrial enzyme malate dehydrogenase 2 (MDH2), which functions in the tricarboxylic acid (TCA) cycle and the malate-aspartate shuttle. We confirmed in cellulo RNA binding of MDH2 using orthogonal biochemical assays and performed enhanced cross-linking and immunoprecipitation (eCLIP) to identify the cellular RNAs associated with endogenous MDH2. Surprisingly, MDH2 preferentially binds cytosolic over mitochondrial RNAs, although the latter are abundant in the milieu of the mature protein. Subcellular fractionation followed by RNA-binding assays revealed that MDH2-RNA interactions occur predominantly outside of mitochondria. We also found that a cytosolically retained N-terminal deletion mutant of MDH2 is competent to bind RNA, indicating that mitochondrial targeting is dispensable for MDH2-RNA interactions. MDH2 RNA binding increased when cellular NAD+ levels (MDH2's cofactor) were pharmacologically diminished, suggesting that the metabolic state of cells affects RNA binding. Taken together, our data implicate an as yet unidentified function of MDH2-binding RNA in the cytosol.

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.

PKR Senses Nuclear and Mitochondrial Signals by Interacting with Endogenous Double-Stranded RNAs.

Protein kinase RNA-activated (PKR) induces immune response by sensing viral double-stranded RNAs (dsRNAs). However, growing evidence suggests that PKR can also be activated by endogenously expressed dsRNAs. Here, we capture these dsRNAs by formaldehyde-mediated crosslinking and immunoprecipitation sequencing and find that various noncoding RNAs interact with PKR. Surprisingly, the majority of the PKR-interacting RNA repertoire is occupied by mitochondrial RNAs (mtRNAs). MtRNAs can form intermolecular dsRNAs owing to bidirectional transcription of the mitochondrial genome and regulate PKR and eIF2α phosphorylation to control cell signaling and translation. Moreover, PKR activation by mtRNAs is counteracted by PKR phosphatases, disruption of which causes apoptosis from PKR overactivation even in uninfected cells. Our work unveils dynamic regulation of PKR even without infection and establishes PKR as a sensor for nuclear and mitochondrial signaling cues in regulating cellular metabolism.

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.

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/.

A conserved isoleucine in the binding pocket of RIG-I controls immune tolerance to mitochondrial RNA.

RIG-I is a cytosolic receptor of viral RNA essential for the immune response to numerous RNA viruses. Accordingly, RIG-I must sensitively detect viral RNA yet tolerate abundant self-RNA species. The basic binding cleft and an aromatic amino acid of the RIG-I C-terminal domain(CTD) mediate high-affinity recognition of 5'triphosphorylated and 5'base-paired RNA(dsRNA). Here, we found that, while 5'unmodified hydroxyl(OH)-dsRNA demonstrated residual activation potential, 5'-monophosphate(5'p)-termini, present on most cellular RNAs, prevented RIG-I activation. Determination of CTD/dsRNA co-crystal structures and mutant activation studies revealed that the evolutionarily conserved I875 within the CTD sterically inhibits 5'p-dsRNA binding. RIG-I(I875A) was activated by both synthetic 5'p-dsRNA and endogenous long dsRNA within the polyA-rich fraction of total cellular RNA. RIG-I(I875A) specifically interacted with long, polyA-bearing, mitochondrial(mt) RNA, and depletion of mtRNA from total RNA abolished its activation. Altogether, our study demonstrates that avoidance of 5'p-RNA recognition is crucial to prevent mtRNA-triggered RIG-I-mediated autoinflammation.

Interplay of endonucleolytic and exonucleolytic processing in the 3'-end formation of a mitochondrial nad2 RNA precursor in Arabidopsis.

Initiation and termination of plant mitochondrial transcription are poorly controlled steps. Precursor transcripts are thus often longer than necessary, and 3'-end processing as well as control of RNA stability are essential to produce mature mRNAs in plant mitochondria. Plant mitochondrial 3' ends are determined by 3'-to-5' exonucleolytic trimming until the progression of mitochondrial exonucleases along transcripts is stopped by stable RNA structures or RNA binding proteins. In this analysis, we investigated the function of the endonucleolytic mitochondrial stability factor 1 (EMS1) pentatricopeptide repeat (PPR) protein and showed that it is essential for the production and the stabilization of the mature form of the nad2 exons 1-2 precursor transcript, whose 3' end corresponds to the 5' half of the nad2 trans-intron 2. The accumulation of an extended rather than a truncated form of this transcript in ems1 mutant plants suggests that the role of EMS1 in 3' end formation is not strictly limited to blocking the passage of 3'-5' exonucleolytic activity, but that 3' end formation of the nad2 exons 1-2 transcript involves an EMS1-dependent endonucleolytic cleavage. This study demonstrates that the formation of the 3' end of mitochondrial transcripts may involve an interplay of endonucleolytic and exonucleolytic processing mediated by PPR proteins.

LINEAGE: Label-free identification of endogenous informative single-cell mitochondrial RNA mutation for lineage analysis.

Single-cell RNA-sequencing (scRNA-seq) has become a powerful tool for biomedical research by providing a variety of valuable information with the advancement of computational tools. Lineage analysis based on scRNA-seq provides key insights into the fate of individual cells in various systems. However, such analysis is limited by several technical challenges. On top of the considerable computational expertise and resources, these analyses also require specific types of matching data such as exogenous barcode information or bulk assay for transposase-accessible chromatin with high throughput sequencing (ATAC-seq) data. To overcome these technical challenges, we developed a user-friendly computational algorithm called "LINEAGE" (label-free identification of endogenous informative single-cell mitochondrial RNA mutation for lineage analysis). Aiming to screen out endogenous markers of lineage located on mitochondrial reads from label-free scRNA-seq data to conduct lineage inference, LINEAGE integrates a marker selection strategy by feature subspace separation and de novo "low cross-entropy subspaces" identification. In this process, the mutation type and subspace-subspace "cross-entropy" of features were both taken into consideration. LINEAGE outperformed three other methods, which were designed for similar tasks as testified with two standard datasets in terms of biological accuracy and computational efficiency. Applied on a label-free scRNA-seq dataset of BRAF-mutated cancer cells, LINEAGE also revealed genes that contribute to BRAF inhibitor resistance. LINEAGE removes most of the technical hurdles of lineage analysis, which will remarkably accelerate the discovery of the important genes or cell-lineage clusters from scRNA-seq data.

Mitochondrial double-stranded RNA triggers induction of the antiviral DNA deaminase APOBEC3A and nuclear DNA damage.

APOBEC3A is an antiviral DNA deaminase often induced by virus infection. APOBEC3A is also a source of cancer mutation in viral and nonviral tumor types. It is therefore critical to identify factors responsible for APOBEC3A upregulation. Here, we test the hypothesis that leaked mitochondrial (mt) double-stranded (ds)RNA is recognized as foreign nucleic acid, which triggers innate immune signaling, APOBEC3A upregulation, and DNA damage. Knockdown of an enzyme responsible for degrading mtdsRNA, the exoribonuclease polynucleotide phosphorylase, results in mtdsRNA leakage into the cytosol and induction of APOBEC3A expression. APOBEC3A upregulation by cytoplasmic mtdsRNA requires RIG-I, MAVS, and STAT2 and is likely part of a broader type I interferon response. Importantly, although mtdsRNA-induced APOBEC3A appears cytoplasmic by subcellular fractionation experiments, its induction triggers an overt DNA damage response characterized by elevated nuclear γ-H2AX staining. Thus, mtdsRNA dysregulation may induce APOBEC3A and contribute to observed genomic instability and mutation signatures in cancer.

The MITOCHONDRIAL TRANSCRIPT STABILITY FACTOR 4 (MTSF4) is essential for the accumulation of dicistronic rpl5-cob mRNAs in Arabidopsis thaliana.

The Arabidopsis thaliana genome harbors more than 450 nuclear genes encoding pentatricopeptide repeat (PPR) proteins that operate in the RNA metabolism of mitochondria and/or plastids. To date, the molecular function of many PPR proteins is still unknown. Here we analyzed the nucleus-encoded gene At4g19440 coding for a P-type PPR protein. Knockout of this gene interferes with normal embryo development and seed maturation. Two experimental approaches were applied to overcome lethality and to investigate the outcome of At4g19440 knockout in adult plants. These studies revealed changes in the abundance of several mitochondria-encoded transcripts. In particular, steady-state levels of dicistronic rpl5-cob RNAs were markedly reduced, whereas levels of mature ccmC and rpl2-mttB transcripts were clearly increased. Predictions according to the one repeat to one nucleotide code for PPR proteins indicate binding of the At4g19440 protein to a previously detected small RNA at the 3' termini of the dicistronic rpl5-cob transcripts. This potential interaction indicates a function of this protein in 3' end formation and stabilization of these RNA species, whereas the increase in the levels of the ccmC mRNA along with other mitochondria-encoded RNAs seems to be a secondary effect of At4g19440 knockout. Since the inactivation of At4g19440 influences the stability of several mitochondrial RNAs we call this gene MITOCHONDRIAL TRANSCRIPT STABILITY FACTOR 4 (MTSF4). This factor will be an interesting subject to study opposing effects of a single nucleus-encoded protein on mitochondrial transcript levels.

Structural insights into mammalian mitochondrial translation elongation catalyzed by mtEFG1.

Mitochondria are eukaryotic organelles of bacterial origin where respiration takes place to produce cellular chemical energy. These reactions are catalyzed by the respiratory chain complexes located in the inner mitochondrial membrane. Notably, key components of the respiratory chain complexes are encoded on the mitochondrial chromosome and their expression relies on a dedicated mitochondrial translation machinery. Defects in the mitochondrial gene expression machinery lead to a variety of diseases in humans mostly affecting tissues with high energy demand such as the nervous system, the heart, or the muscles. The mitochondrial translation system has substantially diverged from its bacterial ancestor, including alterations in the mitoribosomal architecture, multiple changes to the set of translation factors and striking reductions in otherwise conserved tRNA elements. Although a number of structures of mitochondrial ribosomes from different species have been determined, our mechanistic understanding of the mitochondrial translation cycle remains largely unexplored. Here, we present two cryo-EM reconstructions of human mitochondrial elongation factor G1 bound to the mammalian mitochondrial ribosome at two different steps of the tRNA translocation reaction during translation elongation. Our structures explain the mechanism of tRNA and mRNA translocation on the mitoribosome, the regulation of mtEFG1 activity by the ribosomal GTPase-associated center, and the basis of decreased susceptibility of mtEFG1 to the commonly used antibiotic fusidic acid.

Bi-allelic variants in the mitochondrial RNase P subunit PRORP cause mitochondrial tRNA processing defects and pleiotropic multisystem presentations.

Human mitochondrial RNase P (mt-RNase P) is responsible for 5' end processing of mitochondrial precursor tRNAs, a vital step in mitochondrial RNA maturation, and is comprised of three protein subunits: TRMT10C, SDR5C1 (HSD10), and PRORP. Pathogenic variants in TRMT10C and SDR5C1 are associated with distinct recessive or x-linked infantile onset disorders, resulting from defects in mitochondrial RNA processing. We report four unrelated families with multisystem disease associated with bi-allelic variants in PRORP, the metallonuclease subunit of mt-RNase P. Affected individuals presented with variable phenotypes comprising sensorineural hearing loss, primary ovarian insufficiency, developmental delay, and brain white matter changes. Fibroblasts from affected individuals in two families demonstrated decreased steady state levels of PRORP, an accumulation of unprocessed mitochondrial transcripts, and decreased steady state levels of mitochondrial-encoded proteins, which were rescued by introduction of the wild-type PRORP cDNA. In mt-tRNA processing assays performed with recombinant mt-RNase P proteins, the disease-associated variants resulted in diminished mitochondrial tRNA processing. Identification of disease-causing variants in PRORP indicates that pathogenic variants in all three subunits of mt-RNase P can cause mitochondrial dysfunction, each with distinct pleiotropic clinical presentations.

Trypanosoma cruzi strain and starvation-driven mitochondrial RNA editing and transcriptome variability.

Trypanosoma cruzi is a unicellular protistan parasitic species that is comprised of strains and isolates exhibiting high levels of genetic and metabolic variability. In the insect vector, it is known to be highly responsive to starvation, a signal for progression to a life stage in which it can infect mammalian cells. Most mRNAs encoded in its mitochondrion require the targeted insertion and deletion of uridines to become translatable transcripts. This study defined differences in uridine-insertion/deletion RNA editing among three strains and established the mechanism whereby abundances of edited (and, thus, translatable) mitochondrial gene products increase during starvation. Our approach utilized our custom T-Aligner toolkit to describe transcriptome-wide editing events and reconstruct editing products from high-throughput sequencing data. We found that the relative abundance of mitochondrial transcripts and the proportion of mRNAs that are edited varies greatly between analyzed strains, a characteristic that could potentially impact metabolic capacity. Starvation typically led to an increase in overall editing activity rather than affecting a specific step in the process. We also determined that transcripts CR3, CR4, and ND3 produce multiple open reading frames that, if translated, would generate different proteins. Finally, we quantitated the inherent flexibility of editing in T. cruzi and found it to be higher relative to that in a related trypanosomatid lineage. Over time, new editing domains or patterns could prove advantageous to the organism and become more widespread within individual transcriptomes or among strains.