scieee AI-readable full text Open interactive document viewer

Polyadenylation and degradation of structurally abnormal mitochondrial tRNAs in human cells

Toompuu, Marina,Tuomela, Tea,Laine, Pia,Paulin, Lars,Dufour, Eric,Jacobs, Howard

Full text

Published online 6 March 2018 Nucleic Acids Research, 2018, Vol. 46, No. 10 5209–5226 doi: 10.1093/nar/gky159 Polyadenylation and degradation of structurally abnormal mitochondrial tRNAs in human cells Marina Toompuu1, Tea Tuomela1, Pia Laine2, Lars Paulin2, Eric Dufour1and Howard T. Jacobs1,2,* 1Faculty of Medicine and Life Sciences, BioMediTech Institute and Tampere University Hospital, FI-33014 University of Tampere, Finland and 2Institute of Biotechnology, FI-00014 University of Helsinki, Finland Received January 17, 2018; Revised February 16, 2018; Editorial Decision February 19, 2018; Accepted February 20, 2018 ABSTRACT RNA 3polyadenylation is known to serve diverse purposes in biology, in particular, regulating mRNA stability and translation. Here we determined that, upon exposure to high levels of the intercalating agent ethidium bromide (EtBr), greater than those required to suppress mitochondrial transcription, mitochondrial tRNAs in human cells became polyadenylated. Relaxation of the inducing stress led to rapid turnover of the polyadenylated tRNAs. The extent, kinetics and duration of tRNA polyadenylation were EtBr dose-dependent, with mitochondrial tRNAs differentially sensitive to the stress. RNA interference and inhibitor studies indicated that ongoing mitochondrial ATP synthesis, plus the mitochondrial poly(A) polymerase and SUV3 helicase were required for tRNA polyadenylation, while polynucleotide phosphorylase counteracted the process and was needed, along with SUV3, for degradation of the polyadenylated tRNAs. Doxycycline treatment inhibited both tRNA polyadenylation and turnover, suggesting a possible involvement of the mitoribosome, although other translational inhibitors had only minor effects. The dysfunctional tRNALeu(UUR) bearing the pathological A3243G mutation was constitutively polyadenylated at a low level, but this was markedly enhanced after doxycycline treatment. We propose that polyadenylation of structurally and functionally abnormal mitochondrial tRNAs entrains their PNPase/SUV3-mediated destruction, and that this pathway could play an important role in mitochondrial diseases associated with tRNA mutations. INTRODUCTION RNA 3polyadenylation has been reported to serve many roles in biology, with a distinction usually drawn between eukaryotes, where poly(A) is considered to play a positive role, facilitating nuclear export, stability and translation, and prokaryotes, where poly(A) is typically used as a tag to mark RNAs for degradation (reviewed in (1–3)). In eukaryotic organelles, notably chloroplasts, the bacterial principle of poly(A)-dependent RNA degradation prevails, although poly(A) plays a more ambiguous and often taxonor even gene-specific role in mitochondria (3–6). In metazoan mitochondria, polyadenylation has been variously inferred to promote either mRNA turnover or stabilization (7–11), translation (9,10), tRNA maturation and/or repair (11,12) and to play a role, most likely an indirect one, in the sensitivity of nuclear DNA to double-strand breaks induced by ionizing radiation (13). The dichotomous effects of 3poly(A) on metazoan mitochondrial mRNAs (stabilization versus destabilization) appear to be transcript-specific, with some stabilized but others destabilized by the inhibition of polyadenylation (8,10,14,15). The addition of A residues is also formally necessary for the creation of some UAA stop codons (16). The enzyme responsible for the synthesis of the poly(A) tails of mitochondrial mRNAs is the mitochondrial poly(A) polymerase (mtPAP, product of the PAPD1 gene). mtPAP may also be involved in the oligouridylation of histone mRNAs in the cytoplasm at the termination of S-phase (17,18). The degradation of human mitochondrial RNAs tagged with poly(A) is catalyzed by the components of the mitochondrial ‘degradosome’, namely SUV3 helicase (SUV3L1 gene product) and polynucleotide phosphorylase (PNPase, PNPT1 gene product, (15,19)). Degradation is believed to initiate when this complex interacts with unfolded poly(A) tails: thus, in the absence of SUV3 function, abnormal polyadenylated RNAs accumulate (20). Oligoadenylated tRNAs are detected when the processing or surveillance enzyme PDE12 is knocked down (12,21). Low level adenylation of misprocessed or truncated RNAs, including tRNAs, has been reported even in control cells (7). Polyadenylation of tRNAs has been documented in bacteria, under the abnormal conditions of deficiency of tRNA processing enzymes (22,23) or deregulation of the poly(A) *To whom correspondence should be addressed. Tel: +358 50 341 2894; Email: how[email protected] C The Author(s) 2018. Published by Oxford University Press on behalf of Nucleic Acids Research. This is an Open Access article distributed under the terms of the Creative Commons Attribution Non-Commercial License (http://creativecommons.org/licenses/by-nc/4.0/), which permits non-commercial re-use, distribution, and reproduction in any medium, provided the original work is properly cited. For commercial re-use, please contact [email protected] Downloaded from https://academic.oup.com/nar/article-abstract/46/10/5209/4922456 by Tampere University and University Hospital user on 14 September 2018 5210 Nucleic Acids Research, 2018, Vol. 46, No. 10 polymerase PAP I (23), leading to tRNA destruction and severely impaired protein synthesis. Defective tRNAs are tagged for degradation by polyadenylation in Escherichia coli (24) and eukaryotic nuclei (25). Polyadenylated tRNAs have also been reported in chloroplasts (26), although their physiological meaning is unknown. In this study, we investigated the effects of ethidium bromide (EtBr), a DNA-intercalating agent that suppresses mitochondrial transcription. EtBr also intercalates into RNA, including tRNA (27), but this intercalation is restricted by the adoption of tertiary structures. In high-salt conditions tRNAs retain their native L-form structure, and are typically able to bind about tenfold less EtBr than when the tertiary structure is disturbed (28,29), with a preferred binding site at the base of the acceptor stem (28,30). In previous studies, we (31) and others (32)haveused EtBr at doses of up to 250 ng/ml to suppress mitochondrial transcription, resulting in the rapid disappearance of mitochondrial mRNAs and the gradual decay of the more stable mitochondrial rRNAs and tRNAs. Noting that higher doses of EtBr should intercalate more effectively into tRNAs, thus distorting their structure, we investigated the effects on mitochondrial tRNA metabolism of subjecting cells to tenfold greater concentrations of the drug than those used previously. Under these conditions, we revealed an unexpected propensity of mitochondrial tRNAs to acquire long poly(A) tails. After withdrawal of the EtBr, polyadenylated tRNAs were rapidly degraded, following a lag phase. We propose that this represents a mitochondrial surveillance system for abnormal tRNAs, similar to the tRNA quality-control process previously documented in bacteria. MATERIALS AND METHODS Cell lines and culture Previously described cell lines were as follows: 143B osteosarcoma cybrid cell lines homoplasmic for wild-type mitochondrial DNA (mtDNA; clone 43) or the 7472insC mutation (clone 47) (33,34), A549 lung carcinoma cells and the GTcybrid cell line derived from the A549 background, containing both np 3243 and np 12300 mutant mtDNA (35). Cells were routinely cultured in DMEM with or without uridine and pyruvate as described previously (34). Oligonucleotides and reagents Custom-designed DNA oligonucleotides were purchased from DNA Technology (Aarhus, Denmark). Their sequences and other relevant information are shown in Supplementary Table S1. Inhibitors used in polyadenylation assays at the final concentrations shown in figures and legends were cordycepin (3deoxyadenosine), CCCP, doxycycline, oligomycin (mixture A, B, C), or thiamphenicol (all from Sigma), or puromycin (InvivoGen). Antibodies were from Abcam: anti-mtPAP (ab154555, used at 1:2000); antiPDE12 (ab87738, 1:500); anti-PNPT1 (PNPase, ab96176, 1:2000) and anti-SUV3L1 (SUV3, ab176854, 1:2000). AntiGAPDH antibody was from Cell Signaling Technology (14C10, #2118, 1:3000), anti-␤-actin from Santa Cruz (C4, sc-47778, 1:200) and HRP-conjugated goat antirabbit (111-035-144) and goat anti-mouse (115-035-146) secondary antibodies (1:20 000) were from Jackson ImmunoResearch Inc., the latter being used for detection of the anti-actin signal. RNA analysis Total RNA was extracted from cells using the Trizol reagent (Life Technologies) under manufacturer’s recommended conditions. RNA was fractionated on 12% polyacrylamide (PAGE)-7 M urea/TBE gels, electroblotted to Zeta-Probe or Zeta-probe GT membranes (BioRad), and hybridized to gene-specific oligoprobes as described previously (34,35)or using Rapid-hyb buffer (GE Healthcare). Signals were visualized by sensitive X-ray film or, where indicated in legends, by phosphorimaging (Typhoon™imaging system, GE Healthcare). For aminoacylation analyses, total RNA was dissolved on ice in 0.1 M sodium acetate, pH 5.2 and fractionated on acidic 6.5% PAGE-7 M urea gels essentially as described previously (36,37). Samples were deacylated by heating for 10 min at 75◦C followed by 30 min at 37◦Cin 1.5 volumes of 0.5 M Tris–HCl, pH 9.0. To analyze the primary structure of tRNA products in EtBr-treated cells, total RNA was isolated, tRNAs deacylated and circularized by T4 RNA ligase (MBI Fermentas) as described (38), reverse-transcribed with oligonucleotides cser1 or cleu1 and PCR-amplified using oligonucleotides cser1 and cser2 for analysis of tRNASer(UCN) or cleu1 and cleu2 for analysis of tRNALeu(UUR), essentially as described previously (31). The high-molecular weight tRNALeu(UUR) products induced by doxycycline treatment were gel-extracted, circularized and amplified similarly. PCR products were cloned (TOPO TA cloning kit, Invitrogen) and inserts Sanger-sequenced using M13 forward primer and standard dye-terminator technology as described previously (39). Analysis of tRNA response to EtBr Cells were seeded at equal densities on 6 cm plates, giving 70–80% confluence, 14–16 h before the experiments. They were then incubated in fresh medium containing different concentrations of EtBr for the times indicated in figures and legends. Inhibitors of the respiratory chain or translation were added to fresh medium (concentration indicated in figures and legends) 1 h before treatment with EtBr. CCCP treatment was also performed in the reverse order, following 1 h of pre-incubation of cells with 5 ␮g/ml EtBr. For siRNA depletion experiments, cells were seeded on 10 cm plates and transfected with 100 nM siRNAs and Lipofectamine RNAiMAX (Invitrogen) under manufacturer’s recommended conditions. Two days after transfection (or on day 3, see legends), cells were seeded at equal densities on 6 cm plates and EtBr was added one day later. Verification of siRNA-mediated knockdown by Western blotting The efficiency of mRNA depletion by siRNA-based RNA interference (see Supplementary Table S1 for siRNA sequences) was analyzed at the protein level by Western blotting. Cells were lysed in PBS containing 1% N-dodecyl-␤- D-maltoside (or 1% SDS where indicated), 1 mM PMSF Downloaded from https://academic.oup.com/nar/article-abstract/46/10/5209/4922456 by Tampere University and University Hospital user on 14 September 2018 Nucleic Acids Research, 2018, Vol. 46, No. 10 5211 and Pierce protease inhibitor (Thermo Scientific), incubated on ice for 30 min and centrifuged for 20 min at 20 000 gmax (40). Protein concentrations were determined by the DC protein assay (BioRad), fractionated on 10% SDS-PAGE (30 ␮g protein per lane) and electroblotted to Whatman®Protran®nitrocellulose membrane (Perkin Elmer). Membranes were blocked in 5% non-fat dried milk/0.1% Tween/TBS followed by incubation with primary and HRP-conjugated secondary antibodies at the dilutions indicated above. Signals were analyzed using the ChemiDoc imaging system (BioRad). Cell growth and viability Cultures of 106cells were seeded in 10 cm cell plates. After 24 h, 2.5 ␮g/ml of EtBr was added and cells were grown for a further 24 h, then harvested and counted using EVE™automated cell counter (NanoEnTek), with cell viability determined by trypan-blue exclusion. Respirometry Plates of 3 ×106cells were treated for 24 h with 2.5 ␮g/ml EtBr. Intact cell respiration in fresh, EtBr-free medium at 37◦C, and in permeabilized cells supplied with cIcIIIor cIV-linked substrate mixtures, was determined using a Hansatech oxygraph as described previously (41), with all reagents sourced from Sigma-Aldrich. To assess coupling of respiration, permeabilized cells supplied with the cI-linked substrate mixture were treated successively with ADP, followed by oligomycin at concentrations ranging from 25 to 75 nM, and FCCP at concentrations ranging from 0.5 to 2 ␮M, in both cases selecting for further analysis the respective concentration leading to maximal inhibition or activation of oxygen consumption. mtDNA copy number analysis Quantification of total mtDNA was performed essentially as described (42), by comparing the efficiencies of amplification of mitochondrial ND1 or ND4 (43) with the nuclear multi-copy gene encoding 18S rRNA (44,45). DNA was extracted with QIAmp DNA mini kit (QIAgen) from cultured cells according to manufacturer’s instructions. 10 ng of DNA was used per qPCR reaction, using primers (300 nM) and fluorogenic probes (50 nM), both from Metabion International, as listed in Supplementary Table S1, and TaqMan®Universal PCR Master Mix (Applied Biosystems). Each reaction was performed in triplicate and in two independent runs, using the following profile: one cycle at 95◦C for 20 s and then 40 cycles at 95◦C for 3 s and 60◦C for 30 s. Threshold cycle numbers (Ct) were calculated with StepOnePlus v3.1 program (Applied Biosystems) and the results from the two runs were averaged. The ND4/18S and ND1/18S ratios were calculated from cycle threshold values, and normalized to the values from untreated cells. Measurement of whole-cell ROS levels by DCF fluorescence Cells were seeded in 12-well plates at 105cells per well (1 ml). After 24 h, 2.5 ␮g/ml EtBr was added and medium diluted to 1.5 ml. After a further 24 h, cells were washed, then stained in 0.5 ml of 2 ␮M CM-H2DCF-DA (ThermoFisher Scientific) in PBS for 60 min at 37◦C in the dark. After washing with 1 ml warmed PBS, cells were detached in 80 ␮l trypsin–EDTA solution (0.25%, Sigma-Aldrich), for 3 min at 37◦C, after which 0.5 ml of fresh medium was added and cells, maintained in the dark throughout, were analyzed with an Accuri®C6 flow cytometer (BD Biosciences) using 488 laser excitation, fast flow and limits set at 30 000 live cells (ROI selection based on forward & side scatter) or 2 min. Corrections for channel crosstalk were estimated from single dye measurements: DCF signal =FL1(533/30) – 2.8% of FL3; EtBr signal =FL3(670/LP) – 1.1% of FL1. Analysis of mitochondrial translation products Mitochondrial translation products were pulse-labeled with 35S-methionine in the presence of emetine, and analyzed as described previously (34) by SDS-14% PAGE. Library preparation for targeted next-generation sequencing Illumina TruSeq truncated reverse primer 5Phos-NNA GATCGGAAGAGCACACGTCTGAACTCCAGTCA C-Amino 3was first adenylated using a 5DNA Adenylation Kit (New England Biolabs). The adenylated oligo (5 pmol) was ligated to total RNA (500 ng) using truncated T4 RNA Ligase 2 (200 U, New England Biolabs) in a 10 ␮l reaction. cDNA was generated using a primer complementary to the ligated reverse primer, with Superscipt II RNA reverse transcriptase (400 U, ThermoFisher Scientific) in a 22 ␮l reaction. Primers were designed for human mitochondrial tRNAs (GenBank entry NC 012920) and the truncated Illumina TruSeq forward adapter sequence 5ACACTCTTTCCCTACACGACGCTCTTCCGATCT 3was added to the 5end of each of the specific primers (see Supplementary Table S1). Amplification of the selected cDNAs was performed using a 2-step PCR approach. In the first step the designed specific primers with overhangs and Illumina Index primers, 10 pmol each, were used for amplification, with 5 ␮l of cDNA as template, using Phusion HotStart II DNA Polymerase (ThermoFisher Scientific) with an activation step of 98◦C for 10 s and 18 cycles of 98◦C for 10 s, 55◦C for 30 s and 72◦C for 10 s. The obtained PCR products were then treated with Exonuclease I and FastaAP Phosphatase (both from ThermoFisher Scientific) to remove unused primers and nucleotides. An aliquot of the first PCR (5 ␮l) was used as template for the second PCR in order to incorporate the full-length Illumina adapters for the final library, using the same PCR cycle as above. The obtained products were pooled and purified using AMPure XP (Beckman Coulter). Size selection was done using BluePippin (Sage Sciences) with a cut-off of 220–550 bp. The final pool was checked on Fragment Analyzer (Advanced Analytical) and concentrations were measured using Qubit (Life Technologies). DNA sequence analysis The Illumina MiSeq System was used to sequence samples in paired end (R1 and R2) manner. Obtained reads Downloaded from https://academic.oup.com/nar/article-abstract/46/10/5209/4922456 by Tampere University and University Hospital user on 14 September 2018 5212 Nucleic Acids Research, 2018, Vol. 46, No. 10 of 326 bp (R1) and 286 bp (R2) were trimmed using cutadapt (v1.7.1; 46) with default parameters, except that minimum read length (m) was set to 50 bp and the quality minimum (q) to 20. In addition, Illumina’s Truseq adapter sequences were removed. Reads were trimmed in paired-end manner to keep both read pairs following the criteria defined above. Trimmed paired-end reads were overlapped using FLASH (fast length adjustment of short reads) tool (47). The BLASTN tool (48) was used to divide sequences into different tRNA sequence pools. Within each tRNA sequence pool, sequences were sorted and unique sequences were counted, applying exclusion criteria as indicated in Tables. Image processing Blot images were optimized for brightness and contrast and cropped, rotated and/or framed for clarity, but no other manipulations were applied. RESULTS High doses of EtBr induce polyadenylation of human mitochondrial tRNAs Treatment with 0.25 ␮g/ml EtBr is sufficient to inhibit mitochondrial RNA synthesis in 143B osteosarcoma cell cybrids, as we showed previously (31). In preliminary experiments, we observed that treatment of cells with a 10-fold higher dose of the drug resulted in altered mobility of mitochondrial tRNAs on denaturing polyacrylamide gels. We investigated this phenomenon in the same 143B cybrids as in our earlier study, which were homoplasmic, respectively, for wild-type mtDNA and for 7472insC mutant mtDNA (33). Exposure of cells to a concentration of 2.5 ␮g/ml EtBr for 8 h or more resulted in a shift of most of the pool of tRNASer(UCN) to species of higher apparent molecular weight (Figure 1A). The process had already started after 4 h of treatment. By 24 h hardly any of the original tRNA band remained on Northern blots, but much of the modified species had also been degraded (Figure 1A–C). The process had similar kinetics in wild-type and mutant cybrids (Figure 1B, C, Supplementary Figure S1A), despite the lower steady-state abundance of the mutant tRNA ((31), Figure 1A). In contrast, cytosolic tRNASer(UCH) was unaffected by EtBr treatment (Figure 1A). The electrophoretic mobility of a typical mitochondrial mRNA, such as for ND3, was also unmodified, but it was progressively degraded (Supplementary Figure S1B), as expected given the known effects of EtBr on mitochondrial transcription and the relatively short half-lives of mitochondrial mRNAs. The effects of such high levels of EtBr on mitochondrial functions and cell physiology have not previously been tested. We therefore confirmed that treatment with 2.5 ␮g/ml EtBr over 24 h did not impair cell viability (Supplementary Figure S1C), although cell proliferation was curtailed (Supplementary Figure S1C). Under these conditions, the copy number of mtDNA (Supplementary Figure S1D) and the rate of respiration of both intact and permeabilized cells (Supplementary Figure S1E) declined by ∼50%. Although measurements of membrane potential and mitochondrial ROS production by standard fluorimetric methods were not possible, due to the interference of EtBr fluorescence with that of TMRM or MitoSox, we were able to assess these parameters indirectly. The respiration of permeabilized EtBr-treated cells was stimulated by ADP or by the uncoupler FCCP, and the proportionate inhibition by oligomycin was the same as in control cells (Supplementary Figure S1F), indicating that the treatment did not result in uncoupling or major changes in membrane potential. Estimates of whole-cell ROS based on DCF fluorescence, which were not subject to interference from EtBr (Supplementary Figure S1G), indicated that EtBr treatment did result in substantially increased ROS (Supplementary Figure S1G). Finally, in accord with its hypothesized effects on RNA, EtBr treatment abrogated mitochondrial protein synthesis even more effectively than doxycycline (Supplementary Figure S1H). Next, we tested other mitochondrial tRNAs (Figure 2, Supplementary Figure S2A–E). These appeared to vary in their susceptibility to EtBr-induced modification and turnover, and the kinetics thereof. Of those tested, tRNAPhe was the most sensitive to EtBr treatment, with the mature tRNA having almost disappeared after 8 h (Figure 2B, Supplementary Figure S2A), but most of it already modified by 4 h of treatment. At the other end of the spectrum was tRNASer(AGY), approximately 10% of which remained as the mature tRNA even after 24 h, and where the relative amount of the modified tRNA never reached 50% or above. Of the other mitochondrial tRNAs tested, each appeared to respond to treatment with high levels of EtBr with specific characteristics. tRNAGln showed a similar susceptibility to EtBr-induced modification as tRNASer(UCN), while tRNAHis behaved like tRNASer(AGY),andtRNA Lys and tRNATrp were intermediate (Supplementary Figure S2B). Mature tRNALeu(UUR) turned over slowly like tRNASer(AGY) but a larger fraction became modified (Figure 2B, Supplementary Figure S2A). The set represents both clustered and individual tRNAs, tRNAs encoded on each strand of mtDNA, frequent pathological targets, and both the most ‘canonical’ and the most structurally aberrant human mitochondrial tRNAs (tRNALeu(UUR) and tRNASer(AGY), respectively). None of these features correlated in a straightforward manner with susceptibility to EtBr. The extent of modification was EtBr dose-dependent, with 5 or 10 ␮g/ml EtBr giving a more pronounced response than 2.5 ␮g/ml (e.g. Figure 3A, C, Supplementary Figure S2D). Apart from tRNASer(UCN), which differed only in abundance, other tRNAs behaved similarly in wild-type and 7472insC mutant cybrids (Supplementary Figure S2C). A similar response was observed in A549 lung carcinoma cells, following treatment with 5 ␮g/ml EtBr (Supplementary Figure S2E). To identify the nature of the novel tRNA products in EtBr-treated cells, total RNA from wild-type cybrid cells was isolated, tRNAs were deacylated, and sequences analyzed initially by Sanger sequencing of cloned, circularized RNAs amplified with tRNA-specific oligonucleotides. For each of tRNASer(UCN) and tRNALeu(UUR) we obtained sequences from >20 clones of 3-tailed molecules. The extensions consisted almost exclusively of poly(A), with tails of up to 61 nt for tRNASer(UCN) andupto38ntfor Downloaded from https://academic.oup.com/nar/article-abstract/46/10/5209/4922456 by Tampere University and University Hospital user on 14 September 2018 Nucleic Acids Research, 2018, Vol. 46, No. 10 5213 AB mt Ser(UCN) EtBr 2.5 μg/mlEtBr 2.5 μg/ml 0 4 6 8 24 0 4 6 8 24 h wt cybrids 7472insC cybrids wt cybrids 7472insC cybrids C wt cybrids 7472insC cybrids total tRNASer(UCN) % 0 4 8 12 16 20 24 28 h % 0 4 8 12 16 20 24 28 h 125 100 75 50 25 0 125 100 75 50 25 0 mature tRNASer(UCN), % of initial amount % modified tRNASer(UCN) total tRNASer(UCN) % of initial amount cyt Ser(UCH) 0 4 8 12 16 20 24 28 h 125 100 75 50 25 0 Figure 1. Mitochondrial tRNASer(UCN) is modified in cells exposed to high levels of EtBr. (A) Northern blots probed for mitochondrial tRNASer(UCN) and cytosolic tRNASer(UCH) from wild-type and 7472insC mutant 143B cybrid cells treated with 2.5 ␮g/ml EtBr for the indicated times. For similar blots, where the positions of size markers were determined, see Supplementary Figure S2F. (B,C) Quantitation by phosphorimaging of Northern blots probed for mitochondrial tRNASer(UCN) during 24 h of EtBr treatment. Data were first adjusted for background, then normalized against 5S rRNA, and finally against the corresponding signals at t=0: (B) levels of mature tRNASer(UCN) and proportion of signal corresponding to modified tRNA migrating at higher apparent molecular weight. For data from a second, independent experiment, as well as data averaged from the two experiments, see Supplementary Figure S1A. (C) Levels of total tRNASer(UCN), in this case averages for the two experiments, showing very similar profiles for wild-type and mutant cybrid cells. tRNALeu(UUR), some of which were truncated at or within the 3-terminal CCA. Polyadenylated tRNAs were detected only in EtBr-treated cells: for comparable analyses from untreated control cells see (39). To obtain a more reliable estimate of the composition of the tails, we employed targeted next-generation sequencing (tNGS) for these and several other mitochondrial tRNAs (Table 1), generating thousands of independent sequence reads (Supplementary Table S2). Non-A residues comprised less than 1% of the total for each of the analyzed tRNAs, but with some variability between them, with tRNASer(UCN) showing the lowest proportion amongst those analyzed, as well as the lowest proportion of truncation at or within the CCA. tNGS revealed 3tails of up to 51 nt (Table 1, Supplementary Table S2), although this is likely an underestimate, due to the inherent bias of the method towards shorter tailed molecules. Circularization and cloning introduce a similar bias. Based on the electrophoretic mobility of RNA size markers, tailed tRNASer(UCN) species ranged in apparent strandlength from the size of the mature tRNA up to approximately 150 nt, after 6 h of treatment with 2.5 ␮g/ml EtBr. Their size increased further by 11 h of such treatment (Supplementary Figure S2F). Although the markers do not permit an accurate extrapolation of the sizes of molecules with unstructured tails, these observations suggest poly(A) lengths of ≥100 nt. Polyadenylated mitochondrial tRNAs are rapidly turned over To investigate the fate of the polyadenylated tRNAs, cells were incubated in fresh medium for various times, following EtBr exposure, and sampled periodically for Northern blotting. After the removal of the drug, and following a short, EtBr-dose-dependent delay of a few hours, polyadenylated tRNAs were rapidly degraded, with a half-life of ∼1– 1.5 h (Figure 3A, B). Note that in this, and all subsequent experiments, tRNASer(UCN) was used as an example of an extensively polyadenylated tRNA, and, where appropriate, tRNALeu(UUR) as an example of a less extensively Downloaded from https://academic.oup.com/nar/article-abstract/46/10/5209/4922456 by Tampere University and University Hospital user on 14 September 2018 5214 Nucleic Acids Research, 2018, Vol. 46, No. 10 EtBr tRNAPhe 0 4 6 8 24 h A tRNALeu(UUR) tRNASer(AGY) tRNASer(AGY) tRNAPhe tRNALeu(UUR) B % 0 4 8 12 16 20 24 28 h 125 100 75 50 25 0 % 0 4 8 12 16 20 24 28 h 125 100 75 50 25 0 modified mature % 0 4 8 12 16 20 24 28 h 125 100 75 50 25 0 total Figure 2. Mitochondrial tRNAs are differentially modified in cells exposed to high levels of EtBr. (A) Northern blots probed for the indicated mitochondrial tRNAs in wild-type 143B cybrid cells treated with EtBr for the indicated times (tRNAPhe –2.5␮g/ml EtBr, tRNALeu(UUR) and tRNASer(AGY) –3.5␮g/ml EtBr. (B) Quantitation by phosphorimaging of Northern blots probed for the indicated mitochondrial tRNAs in wild-type 143B cybrid cells treated with 2.5 ␮g/ml EtBr for the indicated times. Data for mature and modified tRNAs were first adjusted for background, then normalized against 5S rRNA, and finally against the corresponding signals at t=0. In each case the amounts of mature tRNA, the proportion of modified tRNA and the total amount are shown. For data from a second, independent experiment, see Supplementary Figure S2A. Table 1. Composition of 3tails of mitochondrial tRNAs from cells treated with EtBr tRNA Total reads of tailed tRNA % truncated at CCA Maximum tail length % non-A Leu(UUR) 5883, 12073 63, 63 31, 41 0.66, 0.61 Ser(UCN) 1085, 605 6, 0 35, 32 0.30, 0.21 Leu(CUN) 4451 88 44 0.98 Ser(AGY) 12703, 16202 57, 47 48, 51 0.97, 0.78 Lys 250, 11398 84, 77 29, 41 0.60, 0.38 Notes: 1. Based on tNGS of RNA from 143B wild-type cybrid cells, treated with 5 ␮g/ml EtBr for 8 h. 2. Excluding sequences truncated internally beyond CCA (which have been shown previously to contain adenylate tails even in control cells), as well as unprocessed or partially processed transcripts, dimers or other obvious priming artefacts, e.g. containing internal or terminal primer-derived sequences. Allowing single erroneous base-calls from within the body of each tRNA, reflecting the limited accuracy of NGS. 3. Where shown, the two numbers separated by a comma represent analyses of two different cDNA samples. The first counted all tails of 10 nt or more (for primary sequence data see Supplementary Table S2, sheet 1). The second sampling produced a much greater total read number, so in this case only molecules with tails of ≥20 nt were analyzed (primary sequence data in Supplementary Table S2, sheet 2). Downloaded from https://academic.oup.com/nar/article-abstract/46/10/5209/4922456 by Tampere University and University Hospital user on 14 September 2018 Nucleic Acids Research, 2018, Vol. 46, No. 10 5215 B 0 4 8 1 2 3.5 5 h 2.5 μg/ml EtB r EtBr recovery tRNALeu(UUR) 0 4 8 1 2 3.5 5 6.5 8 A >345 121 72 5S 5S 2.5 μg/ml EtBr recovery 5 μg/ml EtBr recovery 0 4 8 16 24 40 64 4 8 16 24 40 64 h 0 4 8 1 2 3.5 5 6.5 8 h tRNASer(UCN) tRNASer(UCN) tRNALeu(UUR) EtBr recovery EtBr recovery EtBr recovery EtBr recovery EtBr recovery C EtBr recovery EtBr recovery 2.5 μg/ml 5 μg/ml Figure 3. Polyadenylated mitochondrial tRNAs are rapidly degraded after EtBr removal. Northern blots probed for mitochondrial tRNASer(UCN) or tRNALeu(UUR) as indicated, in wild-type 143B cybrid cells treated with 2.5 or 5 ␮g/ml EtBr for the indicated times, then allowed to recover in the absence of the drug as shown. Where shown, 5S rRNA was used as a loading control. RNA sizes (in nt) were extrapolated from the migration of mature tRNASer(UCN) (72 nt), 5S rRNA (121 nt) and ND3 mRNA (>345 nt). See also Supplementary Figure S3A. polyadenylated tRNA. Although some of the polyadenylated tRNASer(UCN) species may have been trimmed and repaired, their disappearance was not accompanied by rapid restoration of the level of mature tRNA (e.g. see Figure 3A), implying that most of the material was degraded. New transcription is probably sufficient to account for the gradual recovery in the level of mature tRNA seen from 16 h following EtBr removal. Detectable amounts of ND3 mRNA, for example, which is replenished by new transcription, started to reappear after approximately the same recovery time (Supplementary Figure S3B). For tRNALeu(UUR), which was less affected by EtBr, the polyadenylated species were also turned over rapidly once the drug was removed, with the mature tRNA being restored to its starting level only after 64 h of recovery (Figure 3C, Supplementary Figure S3A). The slow rate of recovery implies that few, if any, of the polyadenylated tRNAs were repaired, once EtBr was removed. Enzymes required for mitochondrial tRNA polyadenylation and turnover Our observations concerning the polyadenylation of mitochondrial tRNAs following exposure to EtBr, and their rapid turnover upon removal of the drug, are unprecedented. Many questions are raised regarding the molecular machinery that signals and executes these processes, and their physiological significance. mtPAP is the only RNA polymerase inside human mitochondria that has been previously characterized as being capable of synthesizing nontemplated poly(A) tails (on mitochondrial mRNAs, (14)). The enzyme is also responsible for the addition of the discriminator A during maturation of human mitochondrial tRNATyr (12), and is involved in the maturation of mitochondrial tRNACys in Drosophila (11). The low proportion of non-A residues (≤1%, Table 1) in the EtBr-induced tRNA poly(A) tails, similar to that seen in the 3tails of mitochondrial mRNAs (14), also suggests that mtPAP is involved, rather than a less faithful enzyme such as PNPase, which has been implicated in polyadenylation in chloroplasts (49–51). Unusually, mtPAP exhibits relative insensitivity to the chain terminator cordycepin (3deoxyadenosine), enabling its activity to be discriminated from other cellular RNA polymerases (52). We therefore tested the cordycepin sensitivity of mitochondrial tRNA polyadenylation produced by EtBr treatment. At 20 ␮g/ml, cordycepin had no effect on the polyadenylation of tRNASer(UCN) over 8 h of EtBr exposure (Figure 4A), despite the fact that 20 ␮g/ml cordycepin had a profound effect on the steady-state level of ND3 mRNA (Supplementary Figure S4A), synthesis of which Downloaded from https://academic.oup.com/nar/article-abstract/46/10/5209/4922456 by Tampere University and University Hospital user on 14 September 2018 5216 Nucleic Acids Research, 2018, Vol. 46, No. 10 EtBr 2.5 μg/ml – + + + + – + + + + – + + + + Recovery – 0 5 7 9 – 0 5 7 9 – 0 5 7 9 h A B + – + – + – + – 0 4 6 8 h EtBr 5 μg/ml – + + – + – + EtBr 5 μg/ml mtPAP RNAi con PNP mock RNAi con PNP SUV3 RNAi PDE12 mtPAP untreated EtBr 2.5 μg/ml – + + + + – + + + + – + + + + Recovery – 0 5 7 9 – 0 5 7 9 – 0 5 7 9 h C CDY Figure 4. Enzymes required for mitochondrial tRNA polyadenylation and turnover. Northern blots probed for mitochondrial tRNASer(UCN), following EtBr treatment in 143B wild-type cybrid cells also treated with cordycepin (CDY) or siRNAs against different genes as indicated. (A) Cells pre-treated for 1 h and then continuously throughout the experiment with or without 20 ␮g/ml cordycepin, as shown. (B) Cells pre-treated for 96 h with siRNAs as indicated (see Supplementary Table S1), or mock-transfected, before 6 h of treatment with or without EtBr, as shown. (C) Cells pre-treated for 72 h with siRNAs as indicated (see Supplementary Table S1), or untreated, before 6 h of treatment with or without EtBr, followed by recovery for the indicated times after removal of the drug. In each siRNA experiment shown, knockdown at the protein level was verified by western blots, representative examples of which are shown in Supplementary Figure S4B and D. 72 h or 96 h of knockdown gave essentially identical results, both on Western and Northern blots. Signals detected by phosphorimaging in parts (B, C). Note that we confirmed the effects of knockdown of mtPAP, SUV3 and PNPase using alternative siRNAs (Supplementary Figure S4E, F, G). depends on the cordycepin-sensitive mitochondrial RNA polymerase responsible for global transcription (53,54). In contrast, cordycepin treatment for 8 h had little effect on the steady-state level of tRNASer(UCN) in the absence of EtBr (Supplementary Figure S4A). To test further whether mitochondrial tRNA polyadenylation is due to mtPAP or PNPase, we used siRNA-based RNA interference combined with EtBr treatment. After verifying knockdown at the protein level (Supplementary Figure S4B, S4G), we used Northern blots to examine the effects on the length of poly(A) tails added to tRNASer(UCN) (Figure 4B, Supplementary Figure S4C). RNAi directed against mtPAP consistently resulted in drastic shortening of the poly(A) tails, but not in their complete abolition (Figure 4B, Supplementary Figure S4C, E), whereas RNAi directed against PNPase resulted in tail lengthening (Figure 4B, Supplementary Figure S4C), although this was less evident when EtBr was used at a lower dose, with RNAi for less time (Figure 4C, upper panel, Supplementary Figure S4E). When the two genes were knocked down simultaneously, tails were slightly decreased in length compared with untreated cells, or with cells treated with a control siRNA (Supplementary Figure S4C), consistent with the two enzymes acting antagonistically. We took a similar approach to identify the gene products responsible for turnover of polyadenylated tRNAs following removal of EtBr. mtPAP knockdown had no effect on the turnover of the shortened tRNASer(UCN) tails (Figure 4C, Supplementary Figure S4F), whereas PNPase knockdown delayed and largely prevented the degradation of polyadenylated tRNASer(UCN),(Figure4C, upper panel, Supplementary Figure S4F), with 3 or 4 days of RNAi almost completely inhibiting the process. PNPase has been previously implicated in general RNA turnover in human mitochondria as a component of the ‘degradosome’ (19). We therefore tested the other identified component of this machinery, the RNA helicase SUV3, again first verifying the effect of knockdown at the protein level (Supplementary Figure S4D, G). SUV3 knockdown resulted in shortened poly(A) tails in the presence of EtBr (Figure 4C, upper panel, Supplementary Figure S4F), but also inhibited their subsequent degradation (Figure 4C, upper panel, Supplementary Figure S4F). To exclude off-target effects, we conDownloaded from https://academic.oup.com/nar/article-abstract/46/10/5209/4922456 by Tampere University and University Hospital user on 14 September 2018 Nucleic Acids Research, 2018, Vol. 46, No. 10 5217 firmed these findings in each case, using a second siRNA (Supplementary Table S1, Figure S4E, F). The 2´ phosphodiesterase encoded by PDE12,previously implicated in the turnover of the poly(A) tails of human mitochondrial mRNA, in tRNA repair (12,55) and in the turnover of specific, oligoadenylated tRNAs (21), was tested similarly. Although PDE12 knockdown was effective at the protein level (Supplementary Figure S4D), this had no effect on the turnover of polyadenylated tRNASer(UCN) (Figure 4C, lower panel). However, PDE12 knockdown had reproducible effects on the EtBr-induced adenylation of two other tRNAs (Supplementary Figure S4H, S4I), that were previously shown to be susceptible to oligoadenylation in its absence (21). InthecaseoftRNA Lys, PDE12 knockdown revealed a more slowly migrating form of the tRNA that was seen even prior to exposure to EtBr (blue arrows in Supplementary Figure S4H), which we assume to correspond with the previously reported oligoadenylated species (21). Following EtBr exposure, while the accumulation of polyadenylated species was very similar to that seen in cells treated with a control siRNA, the mature form of the tRNA was massively depleted after EtBr treatment, implying that PDE12 is somehow required for its stabilization, rather than turnover. The polyadenylated tRNALys species were typically degraded more slowly during recovery than were those derived from tRNASer(UCN), but PDE12 knockdown may actually have accelerated this process. The effects on PNP, mtPAP and SUV3 knockdown on tRNALys adenylation were qualitatively similar to those seen for tRNASer(UCN) (Supplementary Figure S4H). In the case of tRNAHis, short-tailed molecules were below the detection limit in PDE12 knockdown cells, prior to EtBr exposure. However, they were induced by EtBr, with PDE12 knockdown facilitating the process, such that almost all of the tRNA was converted to a discrete species migrating more slowly than in cells treated with a control siRNA (Supplementary Figure S4I). PDE12 knockdown did not prevent the turnover of these species during recovery from EtBr. Once again, knockdown of PNP, mtPAP and SUV3 produced qualitatively similar effects on EtBr-induced adenylation and turnover of tRNAHis as seen for tRNASer(UCN) (Supplementary Figure S4I). The overall conclusion is that each of the tRNAs tested has a specific behaviour in response to PDE12 knockdown, EtBr treatment and recovery. Furthermore, the processes of tRNA polyadenylation and oligoadenylation, plus the corresponding deadenylation, appear to operate independently. Production and turnover of polyadenylated tRNAs depends on mitochondrial metabolism Macromolecular synthesis is an energy-requiring process, and the main substrate for polyadenylation is ATP. We confirmed that inhibition of mitochondrial ATP synthesis via treatment either with an inhibitor of ATP synthase (oligomycin, Figure 5A) or an uncoupler (CCCP, Figure 5B, Supplementary Figure S5A) resulted in a complete block of the mitochondrial tRNA polyadenylation seen upon EtBr exposure. CCCP had this effect regardless of whether it was added prior to or following EtBr exposure (Figure 5B, Supplementary Figure S5A). Next we tested the effect of inhibitors of mitochondrial protein synthesis. Doxycycline prevented mitochondrial tRNA polyadenylation in response to EtBr, at doses between 200 and 1000 ␮g/ml (Figure 5C, Supplementary Figure S5B). Paradoxically, at the highest doses used, it promoted a tiny accumulation of modified species comigrating with polyadenylated tRNAs, even in the absence of EtBr (Figure 5C, Supplementary Figure S5B, asterisked species). Thiamphenicol had no effect on mitochondrial tRNASer(UCN) polyadenylation (Figure 5D), while puromycin had a mild inhibitory effect (Figure 5E), but only at the highest dose used (10 ␮g/ml). Both drugs had an inhibitory effect on tRNALeu(UUR) polyadenylation at somewhat lower doses, i.e. 0.5–1 mg/ml thiamphenicol or 2–5 ␮g/ml puromycin (Supplementary Figure S5C, D). Doxycycline added during the recovery phase also blocked the turnover of tRNAs that had been polyadenylated in response to EtBr (Figure 5F, Supplementary Figure S5E). This resulted in the accumulation of longer products, both for tRNASer(UCN) (Figure 5F) and tRNALeu(UUR) (Supplementary Figure S5E). To ascertain whether tRNA polyadenylation under these various conditions was determined only by the extent of tRNA deacylation, we used acidic polyacrylamide gel-blots, probed for tRNAs whose aminoacylation status can be easily visualized by their mobility on such gels, tRNALys and tRNALeu(UUR) (Figure 6). These tRNAs are also less rapidly modified than, for example, tRNAPhe or tRNASer(UCN) (Figures 1and 2, Supplementary Figures S1A, S2), allowing us to investigate the acylation status of these tRNAs during a time when most of the tRNA remained unmodified. Both tRNAs were extensively deacylated within 3-4 h of exposure to EtBr (Figure 6A), corresponding roughly with the period when their polyadenylation started to become substantial (see Figure 2A, B, Supplementary Figure S2B). Doxycycline, despite its inhibitory effect on both polyadenylation (Figure 5C, Supplementary Figure S5B), and turnover (Figure 5F, Supplementary Figure S5E), did not abolish deacylation (Figure 6B), merely delaying it. Thiamphenicol at 200 ␮g/ml, which had no effect on tRNALeu(UUR) polyadenylation (Supplementary Figure S5C), also had no effect on deacylation (Figure 6C). Conversely, puromycin, which decreased the extent of tRNALeu(UUR) polyadenylation (Supplementary Figure S5D), also inhibited or delayed deacylation (Figure 6D), although this was pronounced only at the higher dose of the drug. Overall, we conclude that deacylated mitochondrial tRNAs can be polyadenylated, following their accumulation. A pathological mutant mitochondrial tRNA is polyadenylated The foregoing data led us to hypothesize that structural distortion of mitochondrial tRNAs by high levels of EtBr induced their exclusion from the pool of translationally competent tRNAs, entraining their polyadenylation and subsequent degradation. If supported, such a mechanism may represent a natural quality-control pathway relevant to Downloaded from https://academic.oup.com/nar/article-abstract/46/10/5209/4922456 by Tampere University and University Hospital user on 14 September 2018 5224 Nucleic Acids Research, 2018, Vol. 46, No. 10 7. Slomovic,S., Laufer,D., Geiger,D. and Schuster,G. (2005) Polyadenylation and degradation of human mitochondrial RNA: the prokaryotic past leaves its mark. Mol. Cell. Biol.,25, 6427–6435. 8. Nagaike,T., Suzuki,T., Katoh,T. and Ueda,T. (2005) Human mitochondrial mRNAs are stabilized with polyadenylation regulated by mitochondria-specific poly(A) polymerase and polynucleotide phosphorylase. J. Biol. Chem.,280, 19721–19727. 9. Wydro,M., Bobrowicz,A., Temperley,R.J., Lightowlers,R.N. and Chrzanowska-Lightowlers,Z.M. (2010) Targeting of the cytosolic poly(A) binding protein PABPC1 to mitochondria causes mitochondrial translation inhibition. Nucleic Acids Res.,38, 3732–3742. 10. Wilson,W.C., Hornig-Do,H.T., Bruni,F., Chang,J.H., Jourdain,A.A., Martinou,J.C., Falkenberg,M., Spahr,H., Larsson,N.G., Lewis,R.J. et al. (2014) A human mitochondrial poly(A) polymerase mutation reveals the complexities of post-transcriptional mitochondrial gene expression. Hum. Mol. Genet.,23, 6345–6355. 11. Bratic,A., Clemente,P., Calvo-Garrido,J., Maffezzini,C., Felser,A., Wibom,R., Wedell,A., Freyer,C. and Wredenberg,A. (2016) Mitochondrial polyadenylation is a one-step process required for mRNA integrity and tRNA maturation. PLoS Genet.,12, e1006028. 12. Fiedler,M., Rossmanith,W., Wahle,E. and Rammelt,C. (2015) Mitochondrial poly(A) polymerase is involved in tRNA repair. Nucleic Acids Res.,43, 9937–9949. 13. Martin,N.T., Nakamura,K., Paila,U., Woo,J., Brown,C., Wright,J.A., Teraoka,S.N., Haghayegh,S., McCurdy,D., Schneider,M. et al. (2014) Homozygous mutation of MTPAP causes cellular radiosensitivity and persistent DNA double-strand breaks. Cell Death Dis.,5, e1130. 14. Tomecki,R., Dmochowska,A., Gewartowski,K., Dziembowski,A. and Stepien,P.P. (2004) Identification of a novel human nuclear-encoded mitochondrial poly(A) polymerase. Nucleic Acids Res.,32, 6001–6014. 15. Slomovic,S. and Schuster,G. (2008) Stable PNPase RNAi silencing: its effect on the processing and adenylation of human mitochondrial RNA. RNA,14, 310–323. 16. Ojala,D., Montoya,J. and Attardi,G. (1981) tRNA punctuation model of RNA processing in human mitochondria. Nature,290, 470–474. 17. Mullen,T.E. and Marzluff,W.F. (2008) Degradation of histone mRNA requires oligouridylation followed by decapping and simultaneous degradation of the mRNA both 5to 3and 3to 5. Genes Dev.,22, 50–65. 18. Schmidt,M.J., West,S. and Norbury,C.J. (2011) The human cytoplasmic RNA terminal U-transferase ZCCHC11 targets histone mRNAs for degradation. RNA,17, 39–44. 19. Borowski,L.S., Dziembowski,A., Hejnowicz,M.S., Stepien,P.P. and Szczesny,R.J. (2013) Human mitochondrial RNA decay mediated by PNPase-hSuv3 complex takes place in distinct foci. Nucleic Acids Res.,41, 1223–1240. 20. Szczesny,R.J., Borowski,L.S., Brzezniak,L.K., Dmochowska,A., Gewartowski,K., Bartnik,E. and Stepien,P.P. (2010) Human mitochondrial RNA turnover caught in flagranti: involvement of hSuv3p helicase in RNA surveillance. Nucleic Acids Res.,38, 279–298. 21. Pearce,S.F., Rorbach,J., Van Haute,L., D’Souza,A.R., Rebelo-Guiomar,P., Powell,C.A, Brierley,I., Firth,A.E. and Minczuk,M. (2017) Maturation of selected human mitochondrial tRNAs requires deadenylation. eLife,6, e27596. 22. Maes,A., Gracia,C., Hajnsdorf,E. and Regnier,P. (2012) Search for poly(A) polymerase targets in E. coli reveals its implication in surveillance of Glu tRNA processing and degradation of stable RNAs. Mol. Microbiol.,83, 436–451. 23. Mohanty,B.K., Maples,V.F. and Kushner,S.R. (2012) Polyadenylation helps regulate functional tRNA levels in Escherichia coli. Nucleic Acids Res.,40, 4589–4603. 24. Li,Z., Reimers,S., Pandit,S. and Deutscher,M.P. (2002) RNA quality control: degradation of defective transfer RNA. EMBO J.,21, 1132–1138. 25. Anderson,J.T. (2005) RNA turnover: unexpected consequences of being tailed. Curr. Biol.,15, R635–R638. 26. Komine,Y., Kwong,L., Anguera,M.C., Schuster,G. and Stern,D.B. (2000) Polyadenylation of three classes of chloroplast RNA in Chlamydomonas reinhadtii. RNA,6, 598–607. 27. Bittman,R. (1969) Studies of the binding of ethidium bromide to transfer ribonucleic acid: absorption, fluorescence, ultracentrifugation and kinetic investigations. J. Mol. Biol.,46, 251–268. 28. Wells,B.D. and Cantor,C.R. (1977) A strong ethidium binding site in the acceptor stem of most or all transfer RNAs. Nucleic Acids Res.,4, 1667–1680. 29. Ghribi,S., Maurel,M.C., Rougee,M. and Favre,A. (1988) Evidence for tertiary structure in natural single stranded RNAs in solution. Nucleic Acids Res.,16, 1095–1112. 30. Chu,W.C., Liu,J.C. and Horowitz,J. (1997) Localization of the major ethidium bromide binding site on tRNA. Nucleic Acids Res.,25, 3944–3949. 31. Toompuu,M., Yasukawa,T., Suzuki,T., Hakkinen,T., Spelbrink,J.N., Watanabe,K. and Jacobs,H.T. (2002) The 7472insC mitochondrial DNA mutation impairs the synthesis and extent of aminoacylation of tRNASer(UCN) but not its structure or rate of turnover. J. Biol. Chem.,277, 22240–22250. 32. Yasukawa,T., Suzuki,T., Ueda,T., Ohta,S. and Watanabe,K. (2000) Modification defect at anticodon wobble nucleotide of mitochondrial tRNAs(Leu)(UUR) with pathogenic mutations of mitochondrial myopathy, encephalopathy, lactic acidosis, and stroke-like episodes. J. Biol. Chem.,275, 4251–4257. 33. Tiranti,V., Chariot,P., Carella,F., Toscano,A., Soliveri,P., Girlanda,P., Carrara,F., Fratta,G.M., Reid,F.M., Mariotti,C. et al. (1995) Maternally inherited hearing loss, ataxia and myoclonus associated with a novel point mutation in mitochondrial tRNASer(UCN) gene. Hum. Mol. Genet.,4, 1421–1427. 34. Toompuu,M., Tiranti,V., Zeviani,M. and Jacobs,H.T. (1999) Molecular phenotype of the np 7472 deafness-associated mitochondrial mutation in osteosarcoma cell cybrids. Hum. Mol. Genet.,8, 2275–2283. 35. El Meziane,A., Lehtinen,S.K., Hance,N., Nijtmans,L.G., Dunbar,D., Holt,I.J. and Jacobs,H.T. (1998) A tRNA suppressor mutation in human mitochondria. Nat. Genet.,18, 350–353. 36. Enriquez,J.A. and Attardi,G. (1996) Analysis of aminoacylation of human mitochondrial tRNAs. Methods Enzymol.,264, 183–196. 37. El Meziane,A., Lehtinen,S.K., Holt,I.J. and Jacobs,H.T. (1998) Mitochondrial tRNALeu isoforms in lung carcinoma cybrid cells containing the np 3243 mtDNA mutation. Hum. Mol. Genet.,7, 2141–2147. 38. Borner,G.V., Zeviani,M., Tiranti,V., Carrara,F., Hoffmann,S., Gerbitz,K.D., Lochmuller,H., Pongratz,D., Klopstock,T., Melberg,A. et al. (2000) Decreased aminoacylation of mutant tRNAs in MELAS but not in MERRF patients. Hum. Mol. Genet.,9, 467–475. 39. Toompuu,M., Levinger,L.L., Nadal,A., Gomez,J. and Jacobs,H.T. (2004) The 7472insC mtDNA mutation impairs 5and 3processing of tRNA(Ser(UCN)). Biochem. Biophys. Res. Commun.,322, 803–813. 40. Richter,U., Lahtinen,T., Marttinen,P., Suomi,F. and Battersby,B.J. (2015) Quality control of mitochondrial protein synthesis is required for membrane integrity and cell fitness. J. Cell Biol.,211, 373–389. 41. Cannino,G., El-Khoury,R., Pirinen,M., Hutz,B., Rustin,P., Jacobs,H.T. and Dufour,E. (2012) Glucose modulates respiratory complex I activity in response to acute mitochondrial dysfunction. J. Biol. Chem.,287, 38729–38740. 42. Tzoulis,C., Papingji,M., Fiskestrand,T., Roste,L.S. and Bindoff,L.A. (2009) Mitochondrial DNA depletion in progressive external ophthalmoplegia caused by POLG1 mutations. Acta Neurol. Scand., 120, 38–41. 43. He,L., Chinnery,P.F., Durham,S.E., Blakely,E.L., Wardell,T.M., Borthwick,G.M., Taylor,R.W. and Turnbull,D.M. (2002) Detection and quantification of mitochondrial DNA deletions in individual cells by real-time PCR. Nucleic Acids Res.,30, e68. 44. Fu,Y., Lian,Y., Kim,K.S., Zhang,L., Hindle,A.K., Brody,F., Siegel,R.S., McCaffrey,T.A. and Fu,S.W. (2010). BP1 homeoprotein enhances metastatic potential in ER-negative breast cancer. J. Cancer,1, 54–62. 45. Takahashi,M., Haraguchi,A., Tahara,Y., Aoki,N., Fukazawa,M., Tanisawa,K., Ito,T., Nakaoka,T., Higuchi,M. and Shibata,S (2017) Positive association between physical activity and PER3 expression in older adults. Sci. Rep.,7, 39771. 46. Marcel,M. (2011) Cutadapt removes adapter sequences from high-throughput sequencing reads. EMBnet J.,17, 10–12. 47. Magoˇ c,T. and Salzberg,S.L. (2011) FLASH: fast length adjustment of short reads to improve genome assemblies. Bioinformatics,27, 2957–2963. Downloaded from https://academic.oup.com/nar/article-abstract/46/10/5209/4922456 by Tampere University and University Hospital user on 14 September 2018 Nucleic Acids Research, 2018, Vol. 46, No. 10 5225 48. Altschul,S.F., Gish,W., Miller,W., Myers,E.W. and Lipman,D.J. (1990) Basic local alignment search tool. J. Mol. Biol.,215, 403–410. 49. Li,Q.S., Gupta,J.D. and Hunt,A.G. (1998) Polynucleotide phosphorylase is a component of a novel plant poly(A) polymerase. J. Biol. Chem.,273, 17539–17543. 50. Yehudai-Resheff,S., Hirsh,M. and Schuster,G. (2001) Polynucleotide phosphorylase functions as both an exonuclease and a poly(A) polymerase in spinach chloroplasts. Mol. Cell. Biol.,21, 5408–5416. 51. Zimmer,S.L., Schein,A., Zipor,G., Stern,D.B. and Schuster,G. (2009) Polyadenylation in Arabidopsis and Chlamydomonas organelles: the input of nucleotidyltransferases, poly(A) polymerases and polynucleotide phosphorylase. Plant J.,59, 88–99. 52. Hirsch,M. and Penman,S. (1974) Post-transcriptional addition of polyadenylic acid to mitochondrial RNA by a cordycepin-insensitive process. J. Mol. Biol.,83, 131–142. 53. Zylber,E.A., Perlman,S. and Penman,S. (1971) Mitochondrial RNA turnover in the presence of cordycepin. Biochim. Biophys. Acta,240, 588–593. 54. Gelfand,R. and Attardi,G. (1981) Synthesis and turnover of mitochondrial ribonucleic acid in HeLa cells: the mature ribosomal and messenger ribonucleic acid species are metabolically unstable. Mol. Cell. Biol.,1, 497–511. 55. Rorbach,J., Nicholls,T.J. and Minczuk,M. (2011) PDE12 removes mitochondrial RNA poly(A) tails and controls translation in human mitochondria. Nucleic Acids Res.,39, 7750–7763. 56. Kirino,Y., Yasukawa,T., Marjavaara,S.K., Jacobs,H.T., Holt,I.J., Watanabe,K. and Suzuki,T. (2006) Acquisition of the wobble modification in mitochondrial tRNALeu(CUN) bearing the G12300A mutation suppresses the MELAS molecular defect. Hum. Mol. Genet.,15, 897–904. 57. Chomyn,A., Enriquez,J.A., Micol,V., Fernandez-Silva,P. and Attardi,G. (2000) The mitochondrial myopathy, encephalopathy, lactic acidosis, and stroke-like episode syndrome-associated human mitochondrial tRNALeu(UUR) mutation causes aminoacylation deficiency and concomitant reduced association of mRNA with ribosomes. J. Biol. Chem.,275, 19198–19209. 58. Yasukawa,T., Kirino,Y., Ishii,N., Holt,I.J., Jacobs,H.T., Makifuchi,T., Fukuhara,N., Ohta,S., Suzuki,T. and Watanabe,K. (2005) Wobble modification deficiency in mutant tRNAs in patients with mitochondrial diseases. FEBS Lett.,579, 2948–2952. 59. Chopra,S. and Reader,J. (2015) tRNAs as antibiotic targets. Ibba, M., ed., Int. J. Mol. Sci.,16, 321–349. 60. Abbott,J.A., Francklyn,C.S. and Robey-Bond,S.M. (2014). Transfer RNA and human disease. Front. Genet.,5, 158. 61. Li,W., Zhang,Y., Zhang,C., Pei,X., Wang,Z. and Jia,S. (2014) Presence of poly(A) and poly(A)-rich tails in a positive-strand RNA virus known to lack 3 poly(A) tails. Virology,454–455, 1–10. 62. Schneider,A. (2011) Mitochondrial tRNA import and its consequences for mitochondrial translation. Annu. Rev. Biochem.,80, 1033–1053. 63. Bai,Y., Srivastava,S.K., Chang,J.H., Manley,J.L. and Tong,L. (2011) Structural basis for dimerization and activity of human PAPD1, a noncanonical poly(A) polymerase. Mol. Cell,41, 311–320. 64. Munoz-Tello,P., Rajappa,L., Coquille,S. and Thore,S. (2015) Polyuridylation in Eukaryotes: a 3-End Modification Regulating RNA Life. Biomed. Res. Int.,2015, 968127. 65. Borowski,L.S., Szczesny,R.J., Brzezniak,L.K. and Stepien,P.P. (2010) RNA turnover in human mitochondria: more questions than answers? Biochim. Biophys. Acta,1797, 1066–1070. 66. Szczesny,R.J., Borowski,L.S., Malecki,M., Wojcik,M.A., Stepien,P.P. and Golik,P. (2012) RNA degradation in yeast and human mitochondria. Biochim. Biophys. Acta,1819, 1027–1034. 67. Piechota,J., Tomecki,R., Gewartowski,K., Szczesny,R., Dmochowska,A., Kudla,M., Dybczynska,L., Stepien,P.P. and Bartnik,E. (2006) Differential stability of mitochondrial mRNA in HeLa cells. Acta Biochim. Pol.,53, 157–168. 68. Lizano,E., Scheibe,M., Rammelt,C., Betat,H. and Morl,M. (2008) A comparative analysis of CCA-adding enzymes from human and E. coli: differences in CCA addition and tRNA 3-end repair. Biochimie, 90, 762–772. 69. Yue,D., Maizels,N. and Weiner,A.M. (1996) CCA-adding enzymes and poly(A) polymerases are all members of the same nucleotidyltransferase superfamily: characterization of the CCA-adding enzyme from the archaeal hyperthermophile Sulfolobus shibatae. RNA,2, 895–908. 70. Hoffmann,B., Nickel,J., Speer,F. and Schafer,B. (2008) The 3ends of mature transcripts are generated by a processosome complex in fission yeast mitochondria. J. Mol. Biol.,377, 1024–1037. 71. Clemente,P., Pajak,A., Laine,I., Wibom,R., Wedell,A., Freyer,C. and Wredenberg,A. (2015) SUV3 helicase is required for correct processing of mitochondrial transcripts. Nucleic Acids Res.,43, 7398–7413. 72. Wang,D.D., Guo,X.E., Modrek,A.S., Chen,C.F., Chen,P.L. and Lee,W.H. (2014) Helicase SUV3, polynucleotide phosphorylase, and mitochondrial polyadenylation polymerase form a transient complex to modulate mitochondrial mRNA polyadenylated tail lengths in response to energetic changes. J. Biol. Chem.,289, 16727–16735. 73. Dziembowski,A., Piwowarski,J., Hoser,R., Minczuk,M., Dmochowska,A., Siep,M., van der Spek,H., Grivell,L. and Stepien,P.P. (2003) The yeast mitochondrial degradosome. Its composition, interplay between RNA helicase and RNase activities and the role in mitochondrial RNA metabolism. J. Biol. Chem.,278, 1603–1611. 74. Khidr,L., Wu,G., Davila,A., Procaccio,V., Wallace,D. and Lee,W.H. (2008) Role of SUV3 helicase in maintaining mitochondrial homeostasis in human cells. J. Biol. Chem.,283, 27064–27073. 75. Chen,H.W., Rainey,R.N., Balatoni,C.E., Dawson,D.W., Troke,J.J., Wasiak,S., Hong,J.S., McBride,H.M., Koehler,C.M., Teitell,M.A. et al. (2006) Mammalian polynucleotide phosphorylase is an intermembrane space RNase that maintains mitochondrial homeostasis. Mol. Cell. Biol.,26, 8475–8487. 76. Chujo,T., Ohira,T., Sakaguchi,Y., Goshima,N., Nomura,N., Nagao,A. and Suzuki,T. (2012) LRPPRC/SLIRP suppresses PNPase-mediated mRNA decay and promotes polyadenylation in human mitochondria. Nucleic Acids Res.,40, 8033–8047. 77. Chang,J.H. and Tong,L. (2012) Mitochondrial poly(A) polymerase and polyadenylation. Biochim. Biophys. Acta,1819, 992–997. 78. Brodersen,D.E., Clemons,W.M. Jr, Carter,A.P., Morgan-Warren,R.J., Wimberly,B.T. and Ramakrishnan,V. (2000) The structural basis for the action of the antibiotics tetracycline, pactamycin, and hygromycin B on the 30S ribosomal subunit. Cell,103, 1143–1154. 79. Pioletti,M., Schlunzen,F., Harms,J., Zarivach,R., Gluhmann,M., Avila,H., Bashan,A., Bartels,H., Auerbach,T., Jacobi,C. et al. (2001) Crystal structures of complexes of the small ribosomal subunit with tetracycline, edeine and IF3. EMBO J.,20, 1829–1839. 80. Sarkar,S. and Thach,R.E. (1968) Inhibition of formylmethionyl-transfer RNA binding to ribosomes by tetracycline. Proc. Natl. Acad. Sci. U.S.A.,60, 1479–1486. 81. Chopra,I., Hawkey,P.M. and Hinton,M. (1992) Tetracyclines, molecular and clinical aspects. J. Antimicrob. Chemother.,29, 245–277. 82. Wilson,D.N. (2009) The A-Z of bacterial translation inhibitors. Crit. Rev. Biochem. Mol. Biol.,44, 393–433. 83. Drainas,D., Kalpaxis,D.L. and Coutsogeorgopoulos,C. (1987) Inhibition of ribosomal peptidyltransferase by chloramphenicol. Kinetic studies. Eur. J. Biochem.,164, 53–58. 84. Marks,J., Kannan,K., Roncase,E.J., Klepacki,D., Kefi,A., Orelle,C., Vazquez-Laslop,N. and Mankin,A.S. (2016) Context-specific inhibition of translation by ribosomal antibiotics targeting the peptidyl transferase center. Proc. Natl. Acad. Sci. U.S.A.,113, 12150–12155. 85. Nathans,D. (1964) Puromycin inhibition of protein synthesis: incorporation of puromycin into peptide chains. Proc. Natl. Acad. Sci. U.S.A.,51, 585–592. 86. Jourdain,A.A., Koppen,M., Wydro,M., Rodley,C.D., Lightowlers,R.N., Chrzanowska-Lightowlers,Z.M. and Martinou,J.C. (2013) GRSF1 regulates RNA processing in mitochondrial RNA granules. Cell Metab.,17, 399–410. 87. Antonicka,H. and Shoubridge,E.A. (2015) Mitochondrial RNA granules are centers for posttranscriptional RNA processing and ribosome biogenesis. Cell Rep.,10, 920–932. 88. Iborra,F.J., Kimura,H. and Cook,P.R. (2004) The functional organization of mitochondrial genomes in human cells. BMC Biol.,2, 9. 89. Rorbach,J., Richter,R., Wessels,H.J., Wydro,M., Pekalski,M., Farhoud,M., Kuhl,I., Gaisne,M., Bonnefoy,N., Smeitink,J.A. et al. Downloaded from https://academic.oup.com/nar/article-abstract/46/10/5209/4922456 by Tampere University and University Hospital user on 14 September 2018 5226 Nucleic Acids Research, 2018, Vol. 46, No. 10 (2008) The human mitochondrial ribosome recycling factor is essential for cell viability. Nucleic Acids Res,36, 5787–5799. 90. Bogenhagen,D.F., Rousseau,D. and Burke,S. (2008) The layered structure of human mitochondrial DNA nucleoids. J. Biol. Chem., 283, 3665–3675. 91. Simms,C.L., Thomas,E.N. and Zaher,H.S. (2017) Ribosome-based quality control of mRNA and nascent peptides. Wiley Interdiscip. Rev. RNA,8, e1366. 92. Temperley,R.J., Seneca,S.H., Tonska,K., Bartnik,E., Bindoff,L.A., Lightowlers,R.N. and Chrzanowska-Lightowlers,Z.M. (2003) Investigation of a pathogenic mtDNA microdeletion reveals a translation-dependent deadenylation decay pathway in human mitochondria. Hum. Mol. Genet.,12, 2341–2348. 93. Rorbach,J., Yusoff,A.A., Tuppen,H., Abg-Kamaludin,D.P., Chrzanowska-Lightowlers,Z.M., Taylor,R.W., Turnbull,D.M., McFarland,R. and Lightowlers,R.N. (2008) Overexpression of human mitochondrial valyl tRNA synthetase can partially restore levels of cognate mt-tRNAVal carrying the pathogenic C25U mutation. Nucleic Acids Res.,36, 3065–3074. 94. McFarland,R., Clark,K.M., Morris,A.A., Taylor,R.W., Macphail,S., Lightowlers,R.N. and Turnbull,D.M. (2002) Multiple neonatal deaths due to a homoplasmic mitochondrial DNA mutation. Nat. Genet,30, 145–146. 95. Stewart,J.B., Freyer,C., Elson,J.L. and Larsson,N.G. (2008) Purifying selection of mtDNA and its implications for understanding evolution and mitochondrial disease. Nat. Rev. Genet.,9, 657–662. 96. Chrzanowska-Lightowlers,Z.M., Pajak,A. and Lightowlers,R.N. (2011) Termination of protein synthesis in mammalian mitochondria. J. Biol. Chem.,286, 34479–34485. 97. Huynen,M.A., Duarte,I., Chrzanowska-Lightowlers,Z.M. and Nabuurs,S.B. (2012) Structure based hypothesis of a mitochondrial ribosome rescue mechanism. Biol. Direct,7, 14. 98. Kogure,H., Hikawa,Y., Hagihara,M., Tochio,N., Koshiba,S., Inoue,Y., Guntert,P., Kigawa,T., Yokoyama,S. and Nameki,N. (2012) Solution structure and siRNA-mediated knockdown analysis of the mitochondrial disease-related protein C12orf65. Proteins,80, 2629–2642. 99. Feaga,H.A., Quickel,M.D., Hankey-Giblin,P.A. and Keiler,K.C. (2016) Human cells require non-stop ribosome rescue activity in mitochondria. PLoS Genet.,12, e1005964. Downloaded from https://academic.oup.com/nar/article-abstract/46/10/5209/4922456 by Tampere University and University Hospital user on 14 September 2018