Mitochondrial Gene Expression and Beyond-Novel Aspects of Cellular Physiology.

Mitochondria are peculiar organelles whose proper function depends on the crosstalk between two genomes, mitochondrial and nuclear. The human mitochondrial genome (mtDNA) encodes only 13 proteins; nevertheless, its proper expression is essential for cellular homeostasis, as mtDNA-encoded proteins are constituents of mitochondrial respiratory complexes. In addition, mtDNA expression results in the production of RNA molecules, which influence cell physiology once released from the mitochondria into the cytoplasm. As a result, dysfunctions of mtDNA expression may lead to pathologies in humans. Here, we review the mechanisms of mitochondrial gene expression with a focus on recent findings in the field. We summarize the complex turnover of mitochondrial transcripts and present an increasing body of evidence indicating new functions of mitochondrial transcripts. We discuss mitochondrial gene regulation in different cellular contexts, focusing on stress conditions. Finally, we highlight the importance of emerging aspects of mitochondrial gene regulation in human health and disease.

1. Introduction: The Role of Mitochondria within the Cell

Mitochondria are present in the majority of eukaryotic cells, where they play a central role in many processes. They are hubs of energy production as sites of oxidative phosphorylation [1]. They are also engaged in maintaining an appropriate redox state and recycling oxidized electron carriers that are important for cell proliferation [2,3]. Moreover, mitochondria play a key role in cellular signaling, buffering calcium ions and regulating apoptosis processes [4,5]. Mitochondria are also the source of most of the cellular reactive oxygen species (ROS), which may impact on various cellular process [6,7]. Mitochondria take part in response to external stimuli, e.g., viral infection [4], and they are also the place where many basic processes related to innate immune intersect [8,9]. Mitochondrial malfunction is related to numerous pathological states in humans, such as cancer and neurodegeneration [10,11,12]. Furthermore, many human hereditary diseases are caused by mitochondrial dysfunction. Mutations in the mitochondrial genome or nuclear genes that encode mitochondrial proteins lead to primary and secondary mitochondrial diseases connected with improper mitochondrial function [13,14]. In this review, we will focus on the emerging aspects of mitochondrial biology and its implications in human health and we guide interested readers to detailed reviews on mitochondrial disorders [15,16,17].

Mitochondrial biogenesis is far more complex in comparison to other cellular organelles since they possess their own genome, which requires dedicated gene expression machinery [18,19,20]. As a consequence, mitochondrial protein production must be coordinated with the nucleocytoplasmic compartment for proper organelle homeostasis. Many nucleus-encoded proteins are needed for mitochondrial transcription, RNA processing, degradation and translation; all to produce just over a dozen proteins that are parts of the oxidative phosphorylation complexes encoded by the mitochondrial genome [18,19]. Recently, new findings deciphering the link between mitochondrial gene expression and cellular homeostasis were reported. Here, we review the basis of mitochondrial gene regulation with a focus on recent findings, highlighting new RNA-driven mechanisms by which mitochondria contribute to the regulation of cell biology in human health and disease.

2. Mitochondrial Genome—A Simple Molecule with Unorthodox Organization

Mitochondria are semiautonomous organelles possessing their own genome. Human mitochondrial DNA (mtDNA) is a ~16 kb circular molecule composed of double-stranded DNA [21]. Human mtDNA contains only a few genes, but all of them are essential for proper cell function, and mutations in mtDNA have dire consequences [13]. Although human mitochondria contain more than 1100 proteins [22], only 13 of them are encoded by mtDNA. The remaining 24 of the 37 genes present in human mtDNA encode RNAs required for the mitochondrial translational apparatus (22 tRNAs and 2 rRNAs) [19] (Figure 1). The rest of the proteins necessary for mitochondrial function are encoded in the nuclear genome and are imported into mitochondria upon synthesis in the cytoplasm [23]. These proteins are key factors in mtDNA replication, transcription and translation processes [23,24]. The presence of a distinct mitochondrial genome is linked to the endosymbiotic origin of this organelle [25,26]. Recent studies of the mitochondrial genome indicated that mitochondria evolved from a proteobacterial lineage [27]. The human mitochondrial genome was first described in the early 1960s and was one of the first genomes to be fully sequenced [28,29]. Since then, the molecular mechanisms of mtDNA expression have been intensively investigated; nevertheless, our knowledge about these processes is still not complete.

Figure 1

Schematic of human mtDNA and basic steps of mtRNA metabolism. Human mtDNA is a circular double-stranded molecule. Marked are template heavy (H-strand) and light (L-strand) strands of mtDNA. Marked are genes encoding subunits of Complex I (blue), III (orange), IV (violet) and subunits of ATP Synthase (pink). ND4/ND4L and ATP6/ATP8 open reading frames are overlapping and are included in bicistronic mRNAs. Genes encoding rRNAs are colored purple. Genes encoding tRNAs are represented as yellow dots with single letter code depicting their aminoacids. Mitochondrial transcription is initiated from L- and H- strand promoters (ITL and ITH, respectively) located within non-coding regulatory region (NCR). Transcription is driven by DNA-dependent RNA polymerase (POLRMT) with help of its co-factors: TFAM, TFB2M, TEFM and MTRES1. Mitochondrial transcription leads to the formation of three polycistronic transcripts that undergo further processing. Most of mitochondrial mRNAs and rRNAs are punctuated by tRNAs which are excised from the nascent RNA precursors by RNAse P and Z. Arisen transcripts may be subjected to various subsequent processes as maturation, stabilization, degradation, and modification. Selected factors engaged in these processes are presented.

The human mitochondrial genome is tightly packed in nucleoprotein complexes called nucleoids [30]. The application of super-resolution microscopy and mass spectrometry-based techniques significantly contributed to the determination of the nucleoid characteristics [31,32,33,34]. It was shown that each nucleoid contains one or several copies of mtDNA [33,34] and is spatially organized by DNA-binding proteins [31]. More than 20 proteins were identified to copurify with mitochondrial nucleoids [32], among which the most abundant protein mitochondrial transcription factor A (TFAM) that may serve as a mtDNA copy number regulator [35]. Discrete nucleoid complexes are dispersed among the mitochondrial network in the cell [33]. An estimated number of mtDNA molecules per mitochondrion may vary greatly between tissues and range from less than a hundred to thousands [36,37,38]. Each mtDNA molecule can be independently segregated in daughter cells, as nucleoids are exchanged within the mitochondrial network upon fusion/fission processes [39,40]. More than one mtDNA sequence can be present within an individual, leading to so-called mitochondrial heteroplasmy. This can be caused by somatic mutations, heteroplasmy of the oocyte or paternal leakage of mtDNA [41,42]. For a long time, mtDNA was considered to be strictly maternally inherited [43]; nevertheless, few reports suggest paternal contribution in the inheritance of the mitochondrial genome [44,45].

There are several models describing the mechanism of mtDNA replication. Notably, the proposed models are not mutually exclusive. Instead, they were proposed to operate in a complementary mode depending on the tissue, cellular state or energy demand [46,47,48]. For a detailed review of mtDNA replication, please see Holt and Reyes [49] and Gustafsson et al. [50].

Notably, mutations in nuclear-encoded mitochondrial proteins that participate in mtDNA maintenance are implicated in mitochondrial disorders [51]. To date, nearly 300 pathogenic point mutations have been reported in the DNA polymerase gamma (POLG) gene encoding mtDNA replicase, and POLG mutations are one of the main causes of inherited mitochondrial disorders [51,52].

3. Mitochondrial Transcription

MtDNA is composed of heavy (H-strand) and light (L-strand) strands (Figure 1) that can be distinguished by different sedimentation attributes in buoyant density ultracentrifugation due to the uneven distribution of guanines between DNA strands [21]. In humans, the G-rich H-strand serves as a template for the transcription of most mitochondrially encoded genes, while the transcription of the complementary L-strand results in the formation of mostly non-coding RNA (ncRNA) [53]. Only one protein-coding gene and 8 tRNAs are transcribed from the L-strand [53] (Figure 1). The exceptional ~1 kb, non-coding regulatory region (NCR) plays an important role both in mtDNA replication and transcription. The NCR plays an important role as a site of H-strand synthesis initiation [49,54]; moreover, transcription initiation start sites for both mtDNA strands (ITL, light-strand transcription initiation site, ITH, heavy-strand transcription initiation site) are also located within the NCR. The synthesis of RNA starts in both directions within the NCR and leads to the production of long, polycistronic transcripts (Figure 1) [55].

The human transcription apparatus appears to be simple and composed of a monomeric RNA polymerase (POLRMT) that is homologous to bacteriophage polymerases and only a few co-factors [19]. In addition to POLRMT, basic factors participating in the mitochondrial transcription process include TFAM, mitochondrial transcription factor B2 (TFB2M) and mitochondrial transcription elongation factor (TEFM) (Figure 1) [19,56]. Initial in vitro experiments showed that transcription can proceed with the presence of only two proteins: POLRMT and TFB2M. Nevertheless, the efficiency of this transcription apparatus is rather low [57]. Thanks to recent structural studies, the detailed step-by-step initiation of mitochondrial transcription was revealed [56]. It was shown that mitochondrial transcription is initiated by the binding of TFAM to the promoter region and the following recruitment of POLRMT to TFAM-bound mtDNA near the transcription start site [58,59]. Next, TFB2M is recruited to induce and stabilize the open conformation of mtDNA. After initiation of RNA synthesis, TFAM and TFB2M are subsequently released, and the elongation factor TEFM is recruited to enable the transition into transcription elongation [56]. TEFM was proposed to increase the processivity of POLRMT and to enable nearly whole-genome transcription [60,61,62]. Additionally, TEFM was shown to participate in the regulation of a replication/transcription switch [62,63,64]. Recently, TEFM was also proposed to play a role in mitochondrial RNA (mtRNA) processing [65].

Termination of mtDNA transcription is far less understood than its initiation. Mitochondrial transcription termination factor 1 (MTERF1) plays an important role in the termination of mtDNA transcription initiated from the ITL; however, it is not clear which factors take part in the termination of mtDNA transcription initiated from the ITH [56]. MTERF1 binds specific sequences within the tRNALeu gene, causing DNA unwinding and base flipping (i.e., rotation of the nucleotide base outside the DNA double helix), leading to transcription termination [66,67]. MTERF1 was also suggested to prevent interference of the transcription complexes operating in opposite directions [68] and to prevent collision of transcription and replication machineries [69]. Several other MTERF proteins are conserved in vertebrates and plants [70]. Despite playing a role in mitochondrial gene expression, none of the MTERF2-4 proteins were shown to act as a transcription terminator in mammals [71,72,73]. This raises the question whether there are other termination factors yet to be discovered in the human mitochondria.

Among mitochondrial transcription machinery constituents, only TFAM mutations are well established to contribute to human disease [17]. TFAM binds DNA in both sequence-specific and nonspecific manners. The former enables the initiation of mitochondrial transcription [58], and the latter enables the compaction of the mitochondrial genome [74,75]. As mtDNA transcription by POLRMT may serve as a source of RNA primers for mtDNA replication [76], it seems that both sequence-specific and nonspecific manners of TFAM DNA-binding contribute to the maintenance of the mtDNA copy number (reviewed in [77]). TFAM knockout in mice causes severe depletion of mtDNA and is embryonically lethal [78]. In humans, mutations in the TFAM gene were shown to cause a decreased mtDNA copy number and impaired cellular respiration underlying progressive liver failure with neonatal onset [79]. Altered TFAM levels and related changes in the mtDNA copy number were also proposed to be associated with neurodegeneration [80].

Recently, some novel factors participating in mitochondrial transcription regulation were reported. Among them, an interesting example pertains to mitochondrial transcription rescue factor 1 (MTRES1), which was shown to interact with POLRMT and TFAM and to prevent stress-induced loss of mtRNAs by acting at the mitochondrial transcription initiation level [81]. Another example refers to mitochondrial ribosomal protein L7/L12 (MRPL12), whose mutations may lead to respiratory chain deficiency that manifest as growth retardation and neurological deterioration [82]. In addition to being a constituent of the mitochondrial ribosome, MRPL12 also exists in a “free”, ribosome-unrelated matrix pool [83]. MRPL12 was shown to interact with POLRMT to regulate transcription [84]. It was proposed that MRPL12 may serve to coordinate transcription, ribosome biogenesis and/or protein synthesis processes [83,84]. Different roles of MRPL12 in mitochondrial gene regulation may be connected with the presence of two forms, the short and long forms, generated by proteolytic cleavage upon being imported into mitochondria, which may have distinct properties [85]. Several other factors, such as hormones, nuclear transcription factors and chromatin remodeling enzymes, were proposed to regulate mitochondrial transcription either by direct binding to mtDNA or by indirect regulation (reviewed in [86]).

4. Post-Transcriptional Regulation of mtRNAs

Mitochondrial transcription spans almost the entire mitochondrial genome, leading to the formation of three polycistronic transcripts (Figure 1) [19,55]. Two of them, resulting from transcription of either the H- or L-strand, encompass almost the entire genome and carry sequences corresponding to mRNAs, tRNAs, rRNAs and ncRNAs. The third transcript covers genes for two tRNAs (Phe and Val) and both rRNAs [19]. Further steps of cleavage and processing are required to obtain mature, functional RNAs [87].

In nascent precursor transcripts, most mitochondrial mRNAs and rRNAs are punctuated by tRNAs. The first stage of primary RNA processing is the excision of tRNA molecules from the polycistronic transcript, which leads to the formation of immature mRNAs and rRNAs [19,20]. This is an endonucleolytic cleavage mediated by RNAse P and elaC ribonuclease Z 2 (ELAC2), which act at the 5′ and 3′ ends of tRNAs, respectively, to release individual RNAs from polycistronic precursors [88,89]). Mitochondrial RNAse P, unlike the canonical RNase P present in the nucleus, is a complex of three proteins tRNA methyltransferase 10C (TRMT10C), hydroxysteroid 17-beta dehydrogenase 10 (HSD17B10) and protein-only RNase P catalytic subunit (PRORP) and does not contain an RNA component [88]. Released, immature transcripts undergo further processing steps or, as in the case of the majority of ncRNAs, are rapidly removed [20].

Due to the structural and transcriptional organization of the mitochondrial genome, most genes on the same strand are transcribed with equal efficiency, which leads to the formation of equal amounts of their precursors. Nevertheless, the levels of mature RNAs can differ significantly [90]. Although transcription initiated from the L-strand is more frequent than that of the H-strand [55], emerging non-coding RNAs are barely detectable [91]. These two examples show that post-transcriptional processes, especially mtRNA decay, play an important role in controlling mitochondrial gene expression to regulate steady-state levels of specific transcripts [90,92,93].

4.1. Degradation of mtRNAs

The machinery responsible for RNA decay in human mitochondria has remained unknown for many years. Studies performed within the last several years established that the mitochondrial degradosome, a complex of ATP-dependent RNA helicase SUPV3L1 (SUV3) and polynucleotide phosphorylase (PNPase also known as PNPT1), is a key player in mtRNA degradation [91,94]. The importance of mtRNA decay in maintaining mitochondrial homeostasis is underscored by the fact that the disruption of the SUV3 or PNPase gene is embryonically lethal in mice [95,96].

While SUV3 is a helicase that catalyzes the unwinding of RNA duplexes, an activity dependent on ATP hydrolysis by SUV3 [97], PNPase is a phosphorolytic 3′-5′ exoribonuclease, which catalyzes degradation of phosphodiester bonds in RNA [98]. In vitro experiments showed that PNPase is unable to degrade dsRNA substrates unless it forms a complex with SUV3 that unwinds the substrate for degradation [99]. The interaction between SUV3 and PNPase is a prerequisite for mtRNA degradation in vivo and occurs locally in D-foci [94]. The discovery of these structures shows that the RNA decay process in human mitochondria is spatially organized. Notably, components of the mitochondrial degradosome differ markedly in the submitochondrial localization. While SUV3 localizes only to the mitochondrial matrix [97], most of the PNPase is found in the mitochondrial intermembrane space [94,96]. Thus, only a fraction of the PNPase localizes to the mitochondrial matrix and cooperates with SUV3 in mtRNA degradation [94], whereas the rest of the PNPase functions in an SUV3-independent manner. Crystal structures of both proteins revealed the presence of some peculiarities [98,100]. Human PNPase, such as PNPases from other organisms, forms a trimer but has an untypical arrangement of RNA-binding domains [98]. Human SUV3 also has some distinctive features in terms of substrate binding, and it was even suggested that Suv3-like proteins may constitute a separate subfamily of helicases [100].

The main role of the degradosome appears to be clearing of non-coding mtRNA species, which arise mostly from transcription of the L-strand. Radiolabeling studies showed that under normal conditions, these RNAs are swiftly degraded [101]. Consequently, their steady-state levels are very low, and it is only when the degradosome function is impaired that they become readily detectable [91,94]. The degradosome complex was also found to be important in mt-mRNA turnover [94,102], 16S rRNA decay [103] and the exonucleolytic processing of the ND6 mRNA precursor [104]. The final products of the degradosome-mediated RNA decay are several nucleotides in length. These short RNA degradation intermediates are probably removed by RNA exonuclease 2 (REXO2), which is a postulated mitochondrial oligoribonuclease [105]. The activity of the mitochondrial degradosome is modulated by mtRNA binding proteins. While the complex of leucine rich pentatricopeptide repeat containing and SRA stem-loop interacting RNA binding proteins (LRPPRC-SLIRP complex) was suggested to suppress the degradosome-mediated decay of mitochondrial protein-coding RNAs [102]. G-rich RNA sequence binding factor 1 (GRSF1) was recently found to enhance the degradosome activity towards mtRNAs containing a G-quadruplex (G4), which are mostly non-coding mtRNAs [106].

4.2. Mitochondrial RNA-Binding Proteins (mtRBPs)

LRPPRC and SLIRP are among the best characterized noncatalytic mitochondrial RNA-binding proteins containing known RNA-interacting domains, pentatricopeptide repeats (PPR) and RRM domain, respectively. LRPPRC and SLIRP form a complex involved in the regulation of mtRNA stability, and the levels of both proteins are mutually dependent; silencing of LRPPRC results in the depletion of SLIRP and vice versa [107,108]. A recent study using RNA UV crosslinking and RNase footprinting procedures revealed that the LRPPRC-SLIRP complex modulates the secondary structures of mitochondrial transcripts, suggesting that this complex may serve as a chaperone for mtRNAs [109]. The presence of LRPPRC was shown to be important for the existence of a nontranslated, mitoribosome-unbound pool of mt-mRNAs [110]. In addition, LRPPRC was shown to be required for efficient polyadenylation of mt-mRNAs [108,110]. In mice, LRPPRC knockout is embryonically lethal [110], and in humans, mutations in the LRPPRC gene underlie Leigh syndrome, French Canadian type [111]. Recent studies suggest that LRPPRC may play a role in other various pathological states in humans, such as tumors or neurodegeneration (reviewed in detail by Cui et al. [112], emphasizing the very important role of this protein. SLIRP was shown to regulate the translation process by mediating the association of mtRNAs with mitoribosomes [113]. Surprisingly, although SLIRP knockout in mice results in extensive loss of mtRNAs, it is manifested only as a minor weight loss of the animals without any other observable phenotypes [113].

Another important RNA-binding protein, GRSF1, a member of the quasi-RRM (qRRM) family of RNA-binding proteins, was originally identified as a cytoplasmic poly(A)+ mRNA binding protein interacting with G-rich sequences [114]. Later, it was found that GRSF1 is targeted to mitochondria where it localizes to RNA-containing granules [115,116]. It was postulated that GRSF1 participates in the initial stages of polycistronic mtRNA precursor processing [116] and in the translation of some mt-mRNAs [115]. Recent findings, however, reported that GRSF1 takes part in the RNA surveillance pathway, showing that GRSF1 cooperates with the mitochondrial degradosome to regulate mtRNAs that contain G4s [106]. Vertebrates’ mitochondrial genomes have exceptional GC skews, i.e., high guanine content on one strand. As a result, transcripts that are produced by transcription of the G-poor template (i.e., L-strand) are G-rich RNAs; thus, they are prone to form G4s structures. Since G4s are stable, their presence in RNA can hinder its degradation; nevertheless, steady-state levels of mt-ncRNAs that can form G4s are extremely low. GRSF1 was found to positively regulate the degradosome-dependent decay of G4-containing mitochondrial non-coding transcripts by binding and melting G4 structures, which in turn facilitates their degradation [106,117].

Recently, two studies reported novel mtRBP engaged in the regulation of mitochondrial gene expression [81,118]. The level of C6orf203/MTRES1 protein was found to be elevated in cells under stress, and this up-regulation of MTRES1 (mitochondrial transcription rescue factor 1) was shown to prevent mitochondrial transcript loss under perturbed mitochondrial gene expression [81]. MTRES1 associates with the mitochondrial transcription machinery and acts by increasing the mitochondrial transcription without influencing the stability of mitochondrial transcripts [81]. The protective function of MTRES1 depends on its RNA-binding ability since the mutated version incapable of RNA-binding does not prevent a decrease in the mitochondrial RNA [81]. MTRES1 was also shown to associate with a large subunit of the mitochondrial ribosome and to influence mitochondrial translation [118]. Interestingly, silencing of MTRES1 causes the down-regulation of transcripts originating only from NCR without influencing other transcripts [81], which cannot explain the decrease in mitochondrial translation observed in MTRES1 knockout [118]. In contrast, it was reported that the depletion of MTRES1 leads to alterations of the mt-mRNAs’ association with the mitoribosome without influencing mitoribosome stability [118]. MTRES1 is an exciting example of mitochondrial RBP that can play a role in the regulation of mitochondrial gene expression at multiple levels. It is tempting to speculate that MTRES1 could serve as a regulatory factor coupling mitochondrial transcription and translation processes. It is possible that MTRES1 could act by interacting with mitochondrial transcription machinery and to facilitate loading of nascent mt-mRNAs on the mitoribosome. It cannot be excluded that MTRES1 could be a key regulator of mitochondrial transcription/translation coupling, especially under stress conditions. Importantly, another mtRBP, MRPL12, was shown to perform double functions in mitochondrial transcription and translation [83,84], highlighting the possible roles of MTRES1.

Another example of important players in the mitochondrial gene expression regulation concerns members of the FASTK family. Fas activated serine/threonine kinase (FASTK) and its homologs FASTKD1-5 are mitochondrially targeted RNA-binding proteins that play various roles in mtRNA metabolism as processing, translation and mitoribosome assembly proteins [104,112,119,120,121]. The FASTK family was reviewed in detail by Jourdain et al. [122]; therefore, we will not focus on these proteins here.

4.3. mtRNA Modifications

Mitochondrial RNAs undergo diverse modifications. One of the most common is adenylation of the 3’ end. Human mt-mRNAs, with the exception of mt-ND6, are adenylated at the 3’ ends, and this reaction is catalyzed by a noncanonical mitochondrial poly(A) polymerase (MTPAP) [123]. The role of mt-mRNAs’ adenylation is not fully understood. For some mt-mRNAs, the addition of adenine residues at the 3’ end is necessary to create a complete termination codon, as it is not encoded in the genome [21]. Initial studies reported some contradictory findings about the role of mt-mRNA polyadenylation, showing that changes in poly(A) tails had diverse effects on mt-mRNAs. Some of the transcripts were up-regulated, some were unaffected, and others were down-regulated [123,124,125,126]. This draws a speculation that individual transcripts may be differentially controlled, as it seems that adenylation may stabilize some transcripts while it may also direct other transcripts for degradation [20]. Notably, the polyadenylation pattern may vary within a cell type, and the same mt-mRNA can be adenylated to different extents in various cell types [127]. This individualized regulation is most likely achieved by transcript-specific protein-RNA interactions. In yeast, it was shown that each mt-mRNA has specific translation coactivators; moreover, it was proposed that they can function as a part of feedback control loops regulating translation efficiency according to current cell demands [128]. Similarly, the presence of the specific translation activator of cytochrome c oxidase I (TACO1) [129] was also detected in human cells. In addition to 3’-end adenylation, it was recently shown that human mt-mRNAs may also undergo other modifications, such as methylation [130,131] and pseudouridylation [132,133], unraveling an additional layer of mitochondrial gene expression regulation.

Mitochondrial tRNAs seem to be the most extensively modified among mtRNAs. Precursor tRNAs are modified at the 3’ end by tRNA nucleotidyl transferase 1 (TRNT1), which adds the CCA sequence not encoded in the genome [134]. They also undergo several other modifications, which are essential for their stability and proper function [135]. Similarly, mt-rRNAs require several chemical modifications for proper folding, stability and correct mitoribosome assembly, underscoring the important role of RNA modifications in orchestrating various mitochondrial processes [24].

The mitochondrial RNA-binding proteome has been intensely studied, and novel members of this group have been continuously discovered. Among them are mtRNA modifying enzymes, mtRNA processing and mitoribosome assembly factors. Here, we delineated only selected aspects of mtRBP-associated regulation of mitochondrial gene expression, mostly related to mtRNA degradation. For further reading on other aspects of mitochondrial RBP function and detailed insight into mt-tRNA and mt-rRNA editing as well as the mitochondrial epitranscriptome, we guide interested readers to other reviews [136,137,138].

5. Mitochondrial Non-Coding Transcripts and Their Implications on Human Health

5.1. Mitochondrial Double-Stranded RNA (mt-dsRNA) and the Innate Immune Response

Convergent transcription, synthesis of antisense RNA or expression of hairpin-containing RNAs can result in the formation of double-stranded RNA (dsRNA) [139]. RNA:RNA molecules seem to have important signaling or regulatory roles [140,141,142], and the function of a particular dsRNA depends on its source and form and the proteins that interacts with it. It is known that the presence of dsRNA can induce the antiviral response and that properly processed dsRNA can regulate gene expression [140,141,142].

The expression of the human mitochondrial genome is especially prone to produce dsRNA molecules, as this genome is an extraordinary example of convergent transcription. This results in the synthesis of complementary RNAs that can subsequently hybridize and form dsRNA. Under normal conditions, most L-strand transcripts are quickly removed [91,94], which prevents the formation of intermolecular dsRNA. Nevertheless, even in physiological conditions, mt-dsRNAs can be detected [143]. We and others have found that mtDNA transcription is a significant source of dsRNA in humans and have shown that mt-dsRNA can play a signaling role and trigger the interferon response [143,144].

The interferon response pathway is a complex process that is part of the innate immunity response, which together with the adaptive system, provides protection from pathogens. This pathway is far from being fully understood; nevertheless, it can be divided into three major steps: (1) detection of a pathogen, (2) induction of interferon (IFN) expression, and (3) up-regulation of IFN-stimulated genes. Consequently, the cell starts producing anti-pathogen agents and reshapes already active processes to combat a pathogen, or at least hampers its spread [145].

The detection step involves recognition of “pathogen-associated molecular patterns” (PAMPs), which are viral or bacterial nuclei acids, or other molecules specific for microorganisms [146]. The recognition is performed by several families of host sensors collectively called pattern recognition receptors (PRRs), such as interferon induced helicase C domain-containing protein 1 (MDA5) and RNA helicase RIG-I. The binding of immunogenic dsRNA by MDA5 and RIG-I initiates an orchestrated cascade of events that results in the up-regulation of interferon- and interferon-stimulated genes [147]. Normally, this pathway is kept silent due to the mechanisms that help to discriminate between host and foreign nucleic acids, remove potentially immunogenic host-produced nucleic acids or keep them away from PAMP receptors. For example, dsRNA, which can activate MDA5 or RIG-I, is localized in the nucleus and mitochondria, precluding its interaction with MDA5 and RIG-I, which are localized in the cytoplasm.

Recently it was found that mt-dsRNA can be released from mitochondria and, once localized to the cytoplasm, can activate interferon-dependent cellular pathways. Notably degradosome constituents SUV3 and PNPase were established as major actors in the regulation of mt-dsRNA [144]. A crucial role was assigned to PNPase in preventing the induction of the interferon response by mt-dsRNA. PNPase controls mt-dsRNA at two different levels that are spatially separated. First, as a component of the mitochondrial degradosome, it prevents the accumulation of dsRNA by degrading mitochondrial antisense transcripts. The mechanism that determines specificity of PNPase towards antisense transcripts has not been fully uncovered. This may result from targeting of PNPase-SUV3 to G4 mtRNAs by GRSF1. Alternatively, but not mutually exclusive, sense transcripts are protected by mtRBPs, while unprotected antisense mtRNAs are exposed to PNPase-SUV3. This function of PNPase requires cooperation with SUV3 and takes place in the mitochondrial matrix. A second function of PNPase in controlling mt-dsRNAs takes place in the intermembrane space, where a fraction of PNPase prevents the release of mt-dsRNA into the cytosol. This function is SUV3-independent because the helicase is absent from the intermembrane space [97]. The dysfunction of SUV3 thus results in the accumulation of mt-dsRNA exclusively in the matrix but does not lead to the release of mt-dsRNA into the cytosol due to the protective activity of PNPase in the intermembrane space. Therefore, silencing of SUV3 does not stimulate interferon expression, while an increase in interferon levels is observed in PNPase-depleted cells. Since inactivation of the degradosome components causes a vast increase in non-coding L-strand transcripts [91,94,106], it is most likely that the mt-dsRNA that accumulates upon degradosome dysfunction results from intermolecular RNA-RNA interactions. The mechanism by which mt-dsRNA exits mitochondria is unknown, and its identification will be an important step in deciphering a role of mt-dsRNA in cell biology.

The protective role of PNPase is of great importance, as revealed by the fact that mutations in PNPT1 (PNPase encoding gene) cause pathogenesis associated with the up-regulation of mt-dsRNA and interferon response [144,148]. A comparison of the type I interferon response across species, which involved 10 mammals and one bird, revealed PNPase as one of 62 core vertebrate interferon-stimulated genes [149]. Moreover, the PNPase protein level was observed to be higher in human cells infected by human respiratory syncytial virus (HRSV), a cause of serious infections in infants [150]. Taken together with the results on the protective role of PNPase in mt-dsRNA release, these data suggest that PNPase can be involved in the suppression of the IFN response. In fact, IFN-desensitization is an important process that enables cells to recover from IFN signaling [145]. Otherwise, prolonged activation of IFN-dependent pathways can lead to permanent dysregulation of cellular homeostasis. At the organismal level, this is manifested by the development of pathologies called type I interferonopathies. An example of such a disease is Aicardi–Goutières syndrome [151].

Interestingly, the contribution of mt-dsRNA to innate immunity was also observed in other studies. It was shown that the immune response can be modulated by mt-dsRNA via protein kinase RNA-activated (PKR). The canonical PKR-induced response comprises sensing of viral dsRNA by PKR following phosphorylation of the enzyme and altering cellular signaling pathways to promote an antiviral defense. Recent findings revealed that PKR can also be induced by endogenous RNAs, primarily by mt-dsRNAs [143]. PKR was proposed to associate with mitochondria, and mt-dsRNA sensing was suggested to take part in the mitochondrial matrix [143]. However, the proposed mechanism assumes that PKR needs to be exported to the cytoplasm to induce downstream signaling pathways. How PKR translocation to and from mitochondria is maintained remains unclear. It cannot be excluded that PKR induction may take place during mt-dsRNA escape into the cytoplasm, for example, upon aberrant PNPase function.

Notably, the role of mt-dsRNA in the modulation of innate immunity is not restricted to humans, as similar responses were reported in other species, suggesting that mt-dsRNAs are important regulators of the immune response across species. A study in a fly model demonstrated that disrupted mtRNA turnover leads to mt-dsRNA accumulation and altered immune response once mt-dsRNAs are released into the cytoplasm [152]. A study conducted in a murine model showed that the lack of p53 protein results in the activation of innate immunity by dsRNA of mitochondrial origin [153]. Consistent with previous findings on human cells [144], the immune response was shown to be dependent on MDA5 and RIG-I activation [153]. Altogether, these studies established mt-dsRNAs as important regulators of the immune response and revealed a new RNA-dependent mechanism by which mitochondria regulate cell fate.

5.2. Mitochondrial Long Non-Coding RNAs (mt-lncRNAs)

While most non-coding transcripts are swiftly degraded after synthesis, a pool of antisense mtRNAs exhibits some stability, resulting in their detection under normal conditions. These transcripts, called mirror RNAs, such as mirror ND2 [91,94], or long non-coding RNAs (ncND5, ncND6, and ncCytB) [154], are complementary to sense transcripts. Whether these mt-lncRNAs play any role remains to be elucidated; nevertheless, differential tissue-specific expression of these RNA species [154] suggests their potential regulatory function. One can imagine that by hybridizing to their coding counterparts, mt-lncRNAs could regulate the stability or translation efficiency of corresponding mtRNAs.

Other examples of lncRNAs possibly originating from mitochondrial transcription were proposed [155,156]. An example is a 2374 nt-long transcript containing an antisense 16S mt-rRNA sequence (inverted repeat) linked to the 5′ end of the 16S mt-rRNA forming a long double-stranded structure called SncmtRNA (sense non-coding mitochondrial RNA). Interestingly, this transcript was reported to be overexpressed in proliferating but not resting tumor cells [155]. Further studies from the same group reported the presence of 2 other lncRNAs containing inverted repeats linked to the 5′ end of the antisense 16S mt-rRNA, called ASncmtRNAs (antisense non-coding mitochondrial RNAs). These transcripts were found in normal proliferating cells and were reported to be down-regulated in tumor cells [156]. Detected transcripts were postulated to be hallmarks of proliferating tumor cells vs nontumor cells, which may play a role in regulating tumor progression [156,157].

A group of lncRNAs mapped to mtDNA was detected in patients with heart disease [158], and these include circulating mitochondria-derived lncRNA LIPCAR (long intergenic non-coding RNA predicting cardiac remodeling) [159]. The 5′ end of this 781nt long chimeric lncRNA maps to the antisense of the mt-Cytb gene, while its 3′ end maps to the antisense of the mt-COX2 gene. The LIPCAR transcript was shown to be up-regulated in the late stages after myocardial infarction and to be elevated in patients with chronic heart failure, serving as a potential biomarker of cardiac remodeling and cardiovascular mortality [159,160].

As presented here, ncRNAs encompass a group of potentially important molecules that can act in various ways to influence mitochondrial function and cell physiology (Figure 2). Their putative functions involve the regulation of mitochondrial translation and mtRNA stability. They may also function as sponges occupying mtRBPs. Moreover, their ability to form mt-dsRNA together with sense transcripts predisposes them to act as signaling molecules modulating the immune response. Altogether, mitochondrial non-coding RNAs extend the repertoire of potential mechanisms by which mitochondrial gene expression is regulated and influences cell physiology. Thus far, antisense transcripts were mostly studied or linked with non-physiological conditions. It remains to be seen whether this class of RNAs has any roles under normal, nonperturbed conditions.

Figure 2

Cellular functions of mitochondria-derived transcripts. MtDNA may be a source and a target for non-coding RNA. Question mark indicates that the process or mechanism is currently undefined. Nucleus-encoded miRNAs (mitomiRs) are suggested to enter mitochondria and regulate mtDNA expression, nevertheless their exact import pathway is currently unexplained. Transcription of mtDNA may result in formation of small non-coding RNAs (mitosRNAs) which may interplay with mt-mRNAs stability. Mechanism of action of mitochondrial small non-coding RNAs is currently undefined. Long non-coding RNAs (mt-lncRNAs) may arise from mtDNA transcription and may regulate stability and transcription of mt-mRNAs or may serve as baits for mitochondrial RNA-binding proteins (mtRBPs). Convergent transcription of mtDNA may result in formation of double-stranded RNA (mt-dsRNA). Under normal conditions, L-strand derived transcripts are swiftly degraded which prevents formation of mt-dsRNAs. Dysfunction of PNPase causes accumulation of mt-dsRNA and its release into cytoplasm. Once released into the cytoplasm, mt-dsRNAs can induce type I interferon response via activation of MDA5 receptors.

6. Emerging Aspects of mtDNA Maintenance and Expression

6.1. mtDNA Polymerases and mtDNA Repair

POLG, responsible for replication of mtDNA, was long believed to be the only DNA polymerase present in mitochondria [161]. In contrast, recent reports suggest that there may be more than one DNA polymerase operating in mammalian mitochondria. Recent studies proposed several additional DNA polymerases to be implicated in mtDNA metabolism; nevertheless, some of those findings require further validation [162]. The most well-documented novel mtDNA polymerase, primase and DNA directed polymerase (PrimPol), takes part both in nuclear and mitochondrial DNA maintenance [163]. PrimPol may act as a primase forming DNA and RNA primers or as a DNA polymerase to extend the DNA primers and it is also able to perform translesion synthesis [164]. It seems that in human mitochondria, the main role of PrimPol is not to prime mtDNA replication but rather to rescue stalled replication forks by bypassing lesions or repriming DNA synthesis at the site of DNA damage. Further polymerases suggested to take part in mtDNA metabolism comprise polymerases Beta (POLB), Zeta (POLZ) and Theta (POLQ) [162]. These polymerases were proposed to participate in mtDNA repair (POLB) [165] (reviewed by Kaufman and Van Houten [166], mtDNA stability (POLZ) [167] and mtDNA maintenance (POLQ) [168].

Importantly, the discovery of novel mtDNA polymerases broadens the horizon for prospective mtDNA repair pathways. Mitochondria are thought to possess a limited DNA repair system comprising mainly base excision repair (BER). In mitochondrial BER, gap filling is driven by POLG assisted by several other proteins engaged in recognition and cleavage of a damaged DNA base, elimination of the abasic site, DNA 5′-end processing to generate substrate for POLG and, finally, ligation [169,170,171,172,173,174]. Other mtDNA repair pathways are being debated; nucleotide excision repair (NER) is thought to be absent from mitochondria, and only a few studies have reported mismatch repair to be present in mitochondria [52]. It is thus possible that novel enzymes could be engaged in these processes, expanding the repertoire of possible mtDNA repair pathways. This is especially important as mtDNA damage may have significant implications for human health [175].

6.2. mtDNA Editing

As mitochondria are organelles separated with double membranes that harbor strict import machinery, for a long time, it has seemed impossible to stably edit mtDNA in vertebrates. While proteins can efficiently be targeted and delivered into mitochondria, nucleic acid uptake remains controversial. As such, import of the nucleic acid component for mtDNA transformation and editing appeared to be a hurdle [176]. This difficulty can be overcome by the use of protein-only nucleases that can efficiently be engineered to target the mitochondrial matrix and induce sequence-specific double-strand breaks in mtDNA molecules [177]. This strategy takes advantage of mtDNA heteroplasmy and is based on the cleavage of mtDNA molecules harboring certain sequences (for example, pathological mutations). As a double-strand break repair pathway is not present in mitochondria, cleaved mtDNA molecules are rapidly degraded [178,179,180]. Remaining, noncleaved mtDNA molecules with normal sequences can then be replicated to reconstitute mtDNA content [181,182,183,184]. Two DNA editing platforms were introduced into mitochondria: transcription activator-like effector nucleases (TALEN) [185] and zing-finger nucleases (ZFN) [183]; these platforms enabled efficient heteroplasmy shift. In 2018, two breakthrough studies published back-to-back by the Minczuk and Moraes labs showed successful mtDNA editing in vivo [186,187]. In a mouse model of heteroplasmic mitochondrial disease, mitochondrially targeted ZFN [186] and TALENs [187] delivered by the adeno-associated virus led to the successful elimination of mutated mtDNA and subsequent reversion of the pathological phenotype [186,187]. These studies may open a new era for the treatment of heteroplasmic mitochondrial diseases.

6.3. RNA Import into Mitochondria

As mentioned earlier, nucleic acid import into mammalian mitochondria arouses controversy. First, it is not fully clear which RNA molecules can be transported into mitochondria; second, the machinery responsible for this process is debated. PNPase, a component of the mitochondrial degradosome complex, was proposed to play a role in RNA import into mitochondria [96,188]. Indeed, a pool of PNPase is found in the mitochondrial intermembrane space [94,96,189], where it could take part in importing RNA. Nevertheless, the exact mechanism of this process is not known, and it is not certain whether PNPase plays this role alone or cooperates with other proteins. While degradation of mtRNAs is a widely accepted function of PNPase, it is unclear how the enzyme mediates the transport of RNA instead of its destruction. Another suggested pathway of RNA import concerns mitochondrial protein translocase complexes TOM and TIM [190]; however, this has not been confirmed in mammalian cells. Another controversy refers to the RNA species that could potentially be imported into mitochondria. A few studies have indicated the mitochondrial import of recombinant heterologous tRNAs comprising the introduction of motifs derived from yeast tRNA [191,192,193]. Of note, tRNAs were shown to be imported into mitochondria of trypanosomes and plants [194,195]. Nevertheless, import of endogenous RNAs H1 RNA (as an RNA component of the canonical RNaseP), 7-2 RNA (as an RNA component of the RNase MRP) and 5S rRNA (as a constituent of the mitoribosome) was recently challenged [176], as those RNAs were suggested to be dispensable in human mitochondria. Indeed, it was shown that mitochondrial RNaseP does not possess an RNA moiety [88], RNase MRP localizes predominantly to the nucleolus [196,197,198] and 5S rRNA is replaced by tRNA in the mitochondrial ribosome [199,200,201]. Establishing the molecular mechanism of mitochondrial RNA import would be of great value, as it could enable the import of Clustered Regularly Interspaced Short Palindromic Repeats (CRISPR) system constituents for mtDNA manipulation. Although obviously revolutionary and beneficial, introduction of this system in mitochondria seems to be challenging, if not impossible [176].

6.4. Small RNAs in Mitochondrial Biology

Transcription of mtDNA may be a source of distinctive mtRNA species. Recent studies reported that a various range of small ncRNAs may emerge form mitochondrial transcription. A study engaging high-throughput sequencing of small RNA (sRNA) revealed that around 3% of whole-cell sRNAs maps uniquely to the mitochondrial genome in 143B human cells [90]. Interestingly, most of reported unique small RNA sequences were derived from tRNA genes [90]. The functional role of these RNA species is not known and whether they can be a part of RNA interference pathway as it was shown for processed nuclear tRNAs is an open question [90,202,203]. Other reported mitochondrial sRNAs comprise small non-coding transcripts called mitosRNAs (mitochondrial genome-encoded small RNAs) [204] and mitochondrial microRNAs (also called mitomiRs) [205]. Notably, whether the novel sRNA species are just a consequence of mtRNA processing and whether they possess any biological function is currently debated. Small non-coding RNAs appear as a novel class of transcripts that may have regulatory functions in mitochondria, nevertheless their exact mode of action and some aspects of their biogenesis are far from well established and evoke controversies. A canonical function of cytoplasmic microRNA (miRNA) is sequence-specific gene silencing in RNA interference pathway [206]. Several reports suggested presence of microRNAs in mitochondria [207,208,209,210] proposed to be either products of mtDNA transcription or an effect of import of cytoplasmic transcripts inside mitochondria [205]. Nevertheless, miRNA biogenesis encompasses multiple steps of processing and require specialized machinery and it is not known which proteins could be engaged in these processes in mitochondria. Moreover, it is not clear how miRNAs could enter mitochondria. These questions need to be answered in order to establish functional potential of mitochondrial small ncRNAs.

Future studies will show whether identified short RNAs have functional importance or they are merely by-products of mtRNA processing and decay. For example, a recently described short RNA called tRNA-like accumulates strongly upon inhibition of mtRNA degradation by the mitochondrial degradosome [106]. This may imply that most, if not all, tRNA-like molecules are just by-products of mtRNA processing which are removed by the degradosome. Or this transcript has functional importance, but its levels need to be strictly controlled.

6.5. 37 or More?

Early studies of the human mitochondrial genome led to the mapping of mtDNA-encoded genes [21]. For a long time, human mtDNA was considered to encode 37 genes; however, a few recent reports suggested that there may be additional coding sequences within human mtDNA. Human mitochondrial genome may harbor nested small open reading frames (sORF) yielding short peptides derived from polycistronic mt-mRNAs. For example, human 16S rRNA gene was suggested to contain sORF encoding 24 amino acid peptide called humanin, which may have cytoprotective properties [211,212]. It is not known, however, whether humanin is translated inside mitochondria or in the cytoplasm and in the latter case, how humanin transcript could be exported from mitochondria and recognized by cytoplasmic translation apparatus [213]. Other mitochondria-derived peptides (MDPs) were proposed to arise from mitochondrial sORFs and were suggested to have regulatory functions within the cell. Among them are mitochondrial open reading frame of the 12S rRNA-c (MOTS-c) and small humanin-like peptides 1–6 (SHLP1-6) located within 12S rRNA and 16S rRNA genes, respectively [214,215,216]. However, synthesis of MOTS-c would require export of corresponding encoding RNA to the cytoplasm, a process of unknown mechanism. Mitochondrial sORF and sRNAs are intriguing in concept, nonetheless, initial reports of them require further substantiation. It cannot be excluded that mitochondrial sRNAs and MDPs are by-products of mitochondrial gene expression. Therefore, proposed mitochondrial sORF and sRNAs require further studies to be fully embraced by the field.