Cytosine-5 RNA methylation links protein synthesis to cell metabolism

Abstract

Posttranscriptional modifications in transfer RNA (tRNA) are often critical for normal development because they adapt protein synthesis rates to a dynamically changing microenvironment. However, the precise cellular mechanisms linking the extrinsic stimulus to the intrinsic RNA modification pathways remain largely unclear. Here, we identified the cytosine-5 RNA methyltransferase NSUN2 as a sensor for external stress stimuli. Exposure to oxidative stress efficiently repressed NSUN2, causing a reduction of methylation at specific tRNA sites. Using metabolic profiling, we showed that loss of tRNA methylation captured cells in a distinct catabolic state. Mechanistically, loss of NSUN2 altered the biogenesis of tRNA-derived noncoding fragments (tRFs) in response to stress, leading to impaired regulation of protein synthesis. The intracellular accumulation of a specific subset of tRFs correlated with the dynamic repression of global protein synthesis. Finally, NSUN2-driven RNA methylation was functionally required to adapt cell cycle progression to the early stress response. In summary, we revealed that changes in tRNA methylation profiles were sufficient to specify cellular metabolic states and efficiently adapt protein synthesis rates to cell stress.

About the Authors:

Nikoletta A. Gkatza

Roles Conceptualization, Formal analysis, Investigation, Methodology, Validation, Visualization, Writing – original draft, Writing – review & editing

Affiliation: Department of Genetics, University of Cambridge, Cambridge, United Kingdom

ORCID logo http://orcid.org/0000-0001-6360-7940

Cecilia Castro

Roles Formal analysis, Investigation, Methodology, Validation

Affiliation: Department of Biochemistry, University of Cambridge, Cambridge, United Kingdom

Robert F. Harvey

Roles Formal analysis, Validation

Affiliation: Medical Research Council Toxicology Unit, University of Cambridge, Cambridge, United Kingdom

Matthias Heiß

Roles Formal analysis, Validation

Affiliation: Department of Chemistry, Ludwig-Maximilians-University Munich, Munich, Germany

ORCID logo http://orcid.org/0000-0002-6770-9511

Martyna C. Popis

Roles Formal analysis, Investigation, Validation

Affiliation: Department of Genetics, University of Cambridge, Cambridge, United Kingdom

Sandra Blanco

Roles Conceptualization, Formal analysis, Investigation

Affiliations Cancer Cell Signaling and Metabolism Lab, Proteomics Unit CIC bioGUNE, Derio, Spain, Molecular Mechanisms Program, Centro de Investigación del Cáncer and Instituto de Biología Molecular y Celular del Cáncer, Consejo Superior de Investigaciones Científicas (CSIC)-University of Salamanca, Salamanca, Spain

ORCID logo http://orcid.org/0000-0001-7352-9392

Susanne Bornelöv

Roles Data curation, Formal analysis, Methodology, Software

Affiliation: Wellcome–Medical Research Council Cambridge Stem Cell Institute, University of Cambridge, Cambridge, United Kingdom

ORCID logo http://orcid.org/0000-0001-9276-9981

Abdulrahim A. Sajini

Roles Investigation, Validation

Affiliation: Department of Biomedical Engineering, Khalifa University of Science and Technology, Abu Dhabi, United Arab Emirates

ORCID logo http://orcid.org/0000-0002-8600-3899

Joseph G. Gleeson

Roles Methodology, Resources

Affiliation: Department of Neurosciences, San Diego School of Medicine, University of California, La Jolla, California, United States of America

Julian L. Griffin

Roles Methodology, Resources

Affiliation: Department of Biochemistry, University of Cambridge, Cambridge, United Kingdom

James A. West

Roles Formal analysis, Methodology

Affiliation: Department of Biochemistry, University of Cambridge, Cambridge, United Kingdom

ORCID logo http://orcid.org/0000-0002-1535-7737

Stefanie Kellner

Roles Formal analysis, Methodology

Affiliation: Department of Chemistry, Ludwig-Maximilians-University Munich, Munich, Germany

ORCID logo http://orcid.org/0000-0003-3224-7502

Anne E. Willis

Roles Conceptualization, Methodology, Supervision, Validation

Affiliation: Medical Research Council Toxicology Unit, University of Cambridge, Cambridge, United Kingdom

Sabine Dietmann

Roles Data curation, Formal analysis, Methodology, Validation

Affiliation: Wellcome–Medical Research Council Cambridge Stem Cell Institute, University of Cambridge, Cambridge, United Kingdom

Michaela Frye

Roles Conceptualization, Formal analysis, Funding acquisition, Investigation, Resources, Supervision, Validation, Writing – original draft

* E-mail: [email protected]

Affiliations Department of Genetics, University of Cambridge, Cambridge, United Kingdom, German Cancer Center (Deutsches Krebsforschungszntrum), Heidelberg, Germany

ORCID logo http://orcid.org/0000-0002-5636-6840

Introduction

Clinical and genetic heterogeneity in diseases such as metabolic disorders and cancer remains a major challenge for targeted therapies. Phenotypic disease variation can be caused by cell type–specific modulation of gene products via both the transcription and translation machineries [1,2]. Recently, the formation of a variety of chemical modifications in RNA emerged as an additional regulatory layer of gene expression programmes [3].

Of the over 170 known RNA modifications, methylation is the most common [4]. In RNA, 5-methylcytosine (m5C) is often required for normal development, and its formation is mediated by at least eight highly conserved enzymes (NSUN1–7 and DNMT2) [5]. For instance, loss-of-function mutations in the human NSUN2 and NSUN3 genes cause neurodevelopmental and mitochondrial disorders, respectively [6–9]. NSUN2 methylates transfer RNAs (tRNAs) site-specifically at either the anticodon or variable loop (VL), thereby protecting tRNAs from endonucleolytic cleavage [10,11]. This protection is important to prevent the accumulation of tRNA-derived fragments (tRFs), which would otherwise inhibit protein synthesis [11–13].

Phenotypically, the loss of NSUN2 thus leads to an expansion of low translating stem and progenitor populations in skin and brain [11,12,14]. Low protein synthesis in slowly cycling or quiescent stem cells saves energy and prevents premature exhaustion [12,15–17] and is required to maintain a fully functional stem cell state by enhancing resilience towards differentiation cues [12].

The precise mechanisms linking external stresses to RNA modifications and protein synthesis have remained largely unknown in mammals. Here, we reveal that one such mechanism is the endogenous and dynamic methylation of RNA by NSUN2. Our data show that expression of NSUN2 is required to metabolically support high protein synthesis rates. By performing a time course that mapped endogenous, stress-induced tRNA methylation, we observed for the first time, to our knowledge, that loss of methylation occurs in a highly site-specific manner as early as 2 hours following an insult. Stress-induced loss of tRNA methylation altered both the fragmentation pattern among tRNAs and protein synthesis rates. Thus, our data reveal how RNA methylation is mechanistically integrated into metabolic homeostasis.

Results

NSUN2 regulation of stem cell differentiation is not controlled via transcription

Disruption of the Nsun2 gene in mice causes global hypomethylation of tRNAs and a developmental growth retardation [5]. The abnormal development of tissues including brain and skin is the result of impaired stem cell differentiation [11,14,18] (S1A and S1B Fig). The expression of NSUN2 is highly dynamic within tissues. For instance, NSUN2 is absent in quiescent stem cells in hair follicle bulges (BGs), steadily increases in progenitor cells in the hair germ (HG), and is highest in the growing (anagen) hair bulb (HB) (Fig 1A and 1B).

[Figure omitted. See PDF.]

Fig 1. Loss of NSUN2 triggers a shift of the metabolic state towards catabolism.

(A, B) Detection of Nsun2 RNA in Nsun2+/+ and Nsun2−/− mouse skin in early (A) and late (B) anagen. Scale bar: 50 μm. (C-E) Transcriptional changes in skin of wild-type (Nsun2+/+) and Nsun2 knockout (Nsun2−/−) mice. Highlighted in red are significant FC expression differences (FDR < 0.05) in hair follicle stem cells (CD34+/ITGA6high) (C), progenitor cells (PCADhigh/ITGA6low) (D), and anagen skin (E). (n = 3–4 mice per genotype and condition). (F-H) Multivariate analyses of data obtained from MS (F) or NMR spectroscopy–based metabolic profiling (G) using mouse back skin (n = 3–5 mice) or human dermal fibroblasts (n = 5 samples per genotype) (H). Model parameters: R2X = 94.5%, R2Y = 99.9%, and Q2 = 95.8% using partial least square discriminant analysis (F), R2X = 70% and Q2 = 30% (G), and R2X = 85.9% and Q2 = 78.6% (H), using principal component analysis. (I-K) Metabolic differences between NSUN2+/− and NSUN2−/− normalised to NSUN2+/+ human dermal fibroblasts relating to the methionine cycle (I), free amino acids (J), and the TCA cycle (K). The underlying data for this figure can be found in S1–S3 Data and S1 File. BG, hair follicle bulge; FC, fold-change; FDR, false discovery rate; HB, hair bulb; IFE, interfollicular epidermis; ITGA6, integrin alpha-6; MS, mass spectrometry; NMR, nuclear magnetic resonance; PC1, Principal Component 1; PC2, Principal Component 2; PCAD, P-cadherin; TCA, tricarboxylic acid.

https://doi.org/10.1371/journal.pbio.3000297.g001

To dissect the underlying molecular pathways that decrease stem cell sensitivity towards differentiation stimuli in the absence of NSUN2, we transcriptionally profiled mouse stem and progenitor cell populations in the hair follicles, as well as hair follicles in the resting (telogen) or growing (anagen) state. Stem and progenitor cells were isolated by flow sorting using the cell surface marker integrin alpha-6 (ITGA6) and then further separated into quiescent stem cells (ITGA6high/CD34+) and activated, cycling progenitors (ITGA6low/P-cadherin [PCAD]+) (S1B–S1D Fig) [19–21]. In Nsun2−/− mouse skin, quiescent stem cells increase at the expense of activated progenitors (S1B Fig) [12,18]. Transcriptional comparison of stem and progenitor cells isolated from wild-type and Nsun2−/− mice revealed that this increase does not appear to be transcriptionally driven (Fig 1C and 1D). Transcriptionally, these cell populations were highly similar, regardless of their expression of NSUN2 (S1 Data).

Similarly, loss of NSUN2 in actively growing hair follicles (anagen) that highly express NSUN2 in wild-type cells (Fig 1B) resulted in the differential expression of just over 100 genes (red) (Fig 1E). Metabolic genes were highly enriched in these few differentially expressed genes, suggesting that changes in metabolism might account for the insensitivity of stem cells to differentiation cues (S1E and S1F Fig).

NSUN2 regulates metabolism by promoting an anabolic cell state

NSUN2 methylates most tRNAs and a smaller number of noncoding and coding RNAs using S-adenosyl-methionine (SAM) as the methyl donor (S1G Fig) [5,22–24]. SAM and its downstream metabolite, S-adenosyl-homocysteine (SAH), are integral to the one-carbon metabolism encompassing the folate and methionine cycles (S1H Fig). The one-carbon metabolism supports multiple physiological processes, including nucleotide biosynthesis (purines and thymidine), amino acid homeostasis, and the redox defence system [25]. Similar to NSUN2 deletion, dysregulation of the one-carbon metabolism impairs foetal growth and has been linked to neurodevelopmental disorders and cancer [25–27].

We therefore asked whether NSUN2 controlled the one-carbon metabolic cycles that generate the methyl donor required for NSUN2 to modify RNA. We used mass spectrometry (MS) and nuclear magnetic resonance (NMR) spectroscopy to measure how metabolites were affected by loss of NSUN2. We compared wild-type and Nsun2−/− mouse skin in anagen, as well as human dermal fibroblasts expressing (NSUN2+/+ and NSUN2+/−) or lacking NSUN2 (NSUN2−/−) (S2 Data) [7]. Multivariate analysis of either the MS or NMR spectroscopic data clearly separated the genotypes (Fig 1F–1H). Therefore, loss of NSUN2 established a distinct cellular metabolic state.

Three major metabolic pathways were affected by deletion of NSUN2: (1) the methionine cycle, (2) amino acid synthetic pathways, and (3) the tricarboxylic acid (TCA) cycle (Fig 1I–1K; S1I–S1N Fig; S2 Data; S3 Data). First, higher levels of metabolites of the methionine cycle strongly indicated that protein degradation was enhanced in the absence of NSUN2. Cellular metabolism up-regulated methionine and SAM, but not SAH, in response to loss of NSUN2 (Fig 1I). Such an alteration in the stochastic ratio of SAM to SAH can reshape the landscape of protein methylation [28]. In particular, increased levels of free symmetric dimethylarginine (SDMA) and asymmetric dimethylarginine (ADMA) can only be caused by enhanced protein degradation, because they are generated solely upon proteolysis of methylated proteins [29]. Second, we found substantial enhancement of free amino acid levels in the absence of NSUN2 (Fig 1J; S1J and S1M Fig), indicating a reduced rate of protein synthesis. Third, the changes we observed in the metabolites of the TCA cycle indicated that upon loss of NSUN2, metabolism rebalanced from oxidative phosphorylation towards glycolysis (Fig 1K).

In conclusion, the metabolic changes we observed included increased protein degradation, suppressed protein synthesis, and enhanced glycolysis and together showed that cells lacking NSUN2 are maintained in a catabolic state (S1O Fig).

Methylation-dependent and -independent NSUN2 functions regulate metabolism

NSUN2-expressing and -lacking mouse epidermal cell populations were metabolically different, but their transcriptomes were highly similar. In contrast, RNA sequencing (RNA-seq) of human NSUN2-expressing and -lacking dermal fibroblasts identified 2,867 differentially expressed genes (S2A and S2B Fig; S4A Data). However, the transcriptome of NSUN2−/− cells remained largely unaltered when we reexpressed the NSUN2 protein (S2C–S2E Fig; S4B Data), indicating that the differences in gene expression were cell line–specific rather than driven by the presence or absence of NSUN2. To confirm that we indeed rescued tRNA methylation, we reexpressed the wild type (NSUN2), an enzymatic dead version of NSUN2 (K190M), or the empty vector (e.vector) as a control in NSUN2−/− cells [30]. Quantitative MS and RNA bisulfite sequencing (BS-seq) confirmed remethylation of NSUN2-specific sites in tRNAs (Fig 2A–2E; S5 Data).

[Figure omitted. See PDF.]

Fig 2. Methylation-dependent and -independent functions of NSUN2.

(A) Schematic representation of NSUN2-methylated tRNA sites in the anticodon loop (C34) and the VL (C46, C47). (B, C) Number of m5C per tRNA in all tRNAs (B) or tRNA leucine (C) quantified by mass spectrometry in NSUN2−/ cells reexpressing NSUN2, the enzymatic dead version of NSUN2 (K190M), or the empty (‘e.’) vector control. *padj < 0.05; ****padj < 0.0001 (ordinary one-way ANOVA, multiple comparisons). (D) Quantification of m5C levels in all rescued tRNAs (padj < 0.05) using RNA bisulfite sequencing. (E) Heatmaps of example tRNAs showing the rescued m5C sites in five replicates of NSUN2+/+ cells or NSUN2−/− cells reexpressing NSUN2, K190M, or the empty vector (‘e.v.’). (F) PCAs of tRFs differentially abundant in NSUN2-overexpressing NSUN2−/− cells. (G, H) Log2 coverage of tRFs smaller than 46 nucleotides (G) or larger than 46 nucleotides (H). (I) Polysome profile of NSUN2-expressing (NSUN2+/+) and -lacking (NSUN2−/−) cells. Shown is one out three replicates. (J, K) Protein synthesis levels measured by flow cytometry using OP-puro in the indicated cells. CHX served as a control. Data represent mean, and error bars are ±SD. Student’s t test. *p < 0.05, **p < 0.01, ****p < 0.0001. (L-N) Metabolic differences between NSUN2−/− cells overexpressing the NSUN2 or K190 protein normalised (‘norm.’) to NSUN2−/− cell infected with the empty vector control (‘e.V.’) relating to the methionine cycle (L), free amino acids (M), and the TCA cycle (N). The underlying data for this figure can be found in S5–S7 Data and S1 File. ADMA, asymmetric dimethylarginine; CHX, cycloheximide; FC, fold-change; m5C, 5-methylcytosine; OP-puro, O-propargyl-puromycin; PCA, principle component analysis; SAH, S-adenosyl-homocysteine; SAM, S-adenosyl-methionine; SDMA, symmetric dimethylarginine; TCA, tricarboxylic acid; tRF, tRNA-derived fragment; tRNA, transfer RNA; VL, variable loop.

https://doi.org/10.1371/journal.pbio.3000297.g002

One functional role of m5C in tRNAs is to protect from endonucleolytic cleavage by angiogenin [10,11]. To measure the production of tRFs, we performed small RNA-seq in the rescued NSUN2−/− cells. Next, we identified all significantly different tRFs between K190M- and NSUN2-expressing cells (padj < 0.05) that were unaltered when K190M was compared to the empty vector control (padj > 0.75) (Fig 2F; S6 Data). Overexpression of the wild-type but not the enzymatic dead NSUN2 protein rescued the formation of tRFs smaller than 46 nucleotides (Fig 2G and 2H), demonstrating that the biogenesis of distinct tRFs was driven by NSUN2-specific methylation.

Angiogenin-mediated tRNA cleavage inhibits global protein synthesis [31]. Therefore, we next asked how mRNA translation was affected by removal of NSUN2. Polysome profiling confirmed differences in the global abundance of polysomes following NSUN2-removal (Fig 2I; S2F Fig). To quantify de novo protein synthesis, we measured the incorporation of the reporter molecule O-propargyl-puromycin (OP-puro) into nascent polypeptides (S2G Fig) [32]. OP-puro forms covalent conjugates that can be imaged by microscopy and quantified by flow cytometry (S2H Fig). Treatment with cycloheximide (CHX), a potent inhibitor of mRNA translation, served as a positive control [33]. Deletion of NSUN2 repressed protein synthesis (Fig 2J), whereas inhibition of angiogenin was sufficient to up-regulate protein synthesis (S2I Fig). Inhibition of mammalian target of rapamycin complex 1 (mTORC1) using rapamycin down-regulated NSUN2 protein expression but reduced protein synthesis similarly in NSUN2-expressing and -lacking cells (S2J–S2L Fig). Thus, loss of NSUN2 led to decreased protein synthesis rate. Accordingly, reexpression of NSUN2 enhanced de novo protein synthesis, but this effect was independent of its methylation activity (Fig 2K). One explanation for why expression of K190M also led to the up-regulation of protein synthesis might be that its binding protected the unmethylated tRNAs from processing.

Finally, we asked whether restoring NSUN2-specific methylation sites also affected the one-carbon metabolism. Although reexpression of NSUN2 or K190M significantly altered the metabolic profile of NSUN2−/− cells (S2M and S2N Fig), NSUN2 expression was not sufficient to reverse the catabolic cell state (Fig 2L–2N; S7 Data).

Together, NSUN2-dependent methylation protected tRNAs from processing into tRFs, and reexpression of NSUN2 in NSUN2−/− cells enhanced protein synthesis in a methylation-independent manner. However, expression of NSUN2 was not sufficient to switch from a catabolic to anabolic cell state.

NSUN2 functions in the dynamic adaptation of protein synthesis in response to stress

One explanation for why up-regulation of NSUN2 was not sufficient to induce an anabolic cell state was that both tRNA methylation and metabolism normally function in response to external cues. Therefore, we next sought to determine the importance of NSUN2-mediated RNA methylation in response to a changing microenvironment. As a stimulus, we chose stress because NSUN2 is required for the proper cellular response upon stress signals in brain and skin in vivo [11,12].

The dynamic regulation of protein synthesis is an integral part of the cellular stress response [34]. Global protein synthesis is repressed in response to oxidative stress, and translation of mRNAs that encode specific stress-related proteins is enhanced (Fig 3A). Therefore, we analysed how global protein synthesis changed over time in response to stress by measuring OP-puro incorporation (Fig 3B). In NSUN2+/+ cells, protein synthesis rates changed dynamically upon sodium arsenite (NaAsO2) treatment and recovered within 4 hours (Fig 3C; S3A Fig). NSUN2-depleted cells showed attenuated changes to protein synthesis rates (Fig 3C and 3D; S3A Fig), which we confirmed using fluorescence imaging (S3B and S3C Fig). Protein synthesis rates of rescued NSUN2−/− cells (NSUN2) were comparable to NSUN2+/+ cells and slightly, but not significantly, reduced when the enzymatic dead version of NSUN2 (K190M) was expressed (Fig 3E).

[Figure omitted. See PDF.]

Fig 3. NSUN2 functions in the cell cycle to adapt dynamic protein synthesis in response to stress.

(A) Schematic representation how oxidative stress modulates global and gene-specific translation. (B) Treatment regime using arsenite to induce stress and OP-puro to measure protein synthesis. (C) Log2 FC of protein synthesis in NSUN2+/+, NSUN2+/− and NSUN2−/− cells in response to stress compared to the untreated controls (‘Ctr’). CHX served as a control. (n = 2–3 samples per time point). (D) Relative protein synthesis levels in response to stress in NSUN2+/+ and −/− cells measured as FC compared to CHX control. (E) Log2 FC of protein synthesis in NSUN2−/− cells rescued with NSUN2 or the enzymatic dead version K190M after exposure to stress at the indicated time points. (F-I) Polysome profile of NSUN2+/+ (F) and NSUN2−/− cells rescued with wt (I) or mutated NSUN2 (K190M) (H). The empty vector (‘e.V.’)-infected cells served as control (G). Shown is one out of two replicates. (J) Gating used for cell cycle analyses using DAPI incorporation. (K-M) Percentage of NSUN2+/+ (grey) and NSUN2−/− (red) cells in G0/G1- (K), S- (L), and G2/M- (M) phases of the cell cycle after treatment with sodium arsenite for the indicated time. (n = 3 samples per time point). Data presented as mean, error bars ± SD. p-Value: two-way ANOVA calculating row (grey; treatment) and column (black; genotype) factor variation. The underlying data for this figure can be found in S1 File. CHX, cycloheximide; eIF2, eukaryotic initiation factor 2; FC, fold-change; Norm. abs., normalised absorbance; OP-puro, O-propargyl-puromycin; wt, wild-type.

https://doi.org/10.1371/journal.pbio.3000297.g003

Polysome profiling confirmed a recovery of mRNA translation between 2 and 4 hours of stress (Fig 3F). We also measured a difference of the heavy polysomes after 4 hours in NSUN2+/+ cells, which appeared to be less evident in NSUN2−/− control or K190M-overexpressing cells (Fig 3F–3H; S3D Fig). NSUN2−/− cells rescued with the wild-type NSUN2 protein responded similarly to NSUN2+/+ cells and recovered the heavy polysomes after 4 hours of stress (Fig 3F and 3I; S3D Fig).

In summary, cell stress caused a strong but temporary reduction of protein synthesis, which was attenuated by loss of NSUN2.

NSUN2 functions in the oxidative stress response by altering the cell cycle phases

To test for the functional relevance of NSUN2-regulated protein synthesis rates, we asked what impact NSUN2 deletion had on cell survival, as oxidative modifications to RNA can result in cell death [35]. In both NSUN2-expressing and -lacking cells, cell survival began to decrease between 2 and 4 hours after stress exposure, but at these early time points there were no differences in cell death (S3E–S3G Fig). We then asked whether cell division was affected by the absence of NSUN2, since the lengths of the cell cycle phases can play important roles in cellular adaptation and response to external stress stimuli [36]. For example, cells exposed to sodium arsenite often slow cell cycle progression to facilitate repair of oxidative lesions [37]. To test whether RNA methylation was required for adapting cell cycle to stress, we measured the cell cycle progression of NSUN2-expressing and -lacking cells upon exposure to arsenite (Fig 3J–3M; S3H Fig). In response to stress, the percentage of NSUN2+/+ cells decreased in the G1/G0-phase but increased in the S-phase and G2/M-phase of the cell cycle (Fig 3K–3M; grey). In contrast, the cell cycle progression of NSUN2−/− cells remained stable (Fig 3K–3M; red), indicating that NSUN2−/− cells failed to adapt the cell cycle phases to the stress stimulus. A tight regulation of global protein synthesis might be needed to avoid accumulation of proteins during cell cycle arrest and repair.

Cytosine-5 RNA methylation is a metabolic sensor of external stress

Since the dynamic regulation of protein synthesis was reduced when cells lacked NSUN2 and the down-regulation of global translation is integral to the cellular stress response, we next asked whether NSUN2 acted as a sensor of the external stress stimulus. We measured how endogenous NSUN2 expression was affected by oxidative stress induced by exposure to sodium arsenite (Fig 4A). NSUN2 expression decreased sharply on RNA and protein levels between 1 and 2 hours of sodium arsenite treatment (Fig 4B and 4C), which coincided with the formation of stress granules in the cytoplasm (S4A Fig). Similarly, one stimulus of UVB exposure was sufficient to rapidly reduce Nsun2 RNA levels in human epidermal and dermal cells (S4B Fig), whereas vehicle control treatments did not affect NSUN2 expression (S4C Fig).

[Figure omitted. See PDF.]

Fig 4. Levels of m5C changes site-specifically and dynamically in response to oxidative stress.

(A) Time course of sodium arsenite treatment. (B) Log2 FC of Nsun2 RNA expression in NSUN2+/+ and NSUN2+/− cells relative to GAPDH and normalised to the untreated control (‘Ctr’). Shown are 3 replicates. (C) Western blot analysis of the indicated proteins using whole cell lysates from NSUN2+/+ and NSUN2−/− cells. Hsp90 served as a loading control. (D,E) Detection of m5C in sodium arsenite–treated and untreated (‘ctr’) NSUN2+/+ and −/− cells using mass spectrometry. (n = 3 samples per time point). (F) Quantification of tRNA methylation percentage using RNA bisulfite sequencing of NSUN2+/+ and NSUN2−/− cells (n = 4 samples per time point). (G) Heatmap of methylation status of individual tRNA molecules shown in (F). (H,I) Quantification (H) and heatmap (I) of methylation changes in the tRNAs LeuCAA and AspGTC in NSUN2+/+ and NSUN2−/− cells. (J) Quantification of methylation in non-tRNA targets. Data represent median in F, H, and J. Error bars are ±SD. p-Values: Student’s t test, *p < 0.05 and **p < 0.01. ***p < 0.001. The underlying data for this figure can be found in S8 and S9 Data and S1 File. FC, fold-change; GAPDH, glyceraldehyde 3-phosphate dehydrogenase; HSP90, heat shock protein 90; m5C, 5-methylcytosine; NPMI, nucleophosmin; tRNA, transfer RNA; VL, variable loop.

https://doi.org/10.1371/journal.pbio.3000297.g004

The nucleolus, where NSUN2 resides, can act as a stress sensor [38,39]. Nucleophosmin (NPMI) is a marker for nucleolar stress, and we observed a rapid, strong down-regulation of both NPMI and NSUN2 upon arsenite treatment (Fig 4C). Additional NSUN family members residing in the mitochondria (NSUN3, NSUN4) and cytoplasm (NSUN6) were similarly repressed in response to arsenite stress (Fig 4C). In conclusion, expression of NSUN2 and at least three of its m5C RNA methyltransferase family members were repressed in response to oxidative stress.

Stress induces a site-specific and dynamic loss of m5C

We then asked how stress-related loss of NSUN2 altered tRNA methylation (Fig 4D). Indeed, MS quantifying total m5C in tRNAs confirmed that cells chronically lacking NSUN2 (NSUN2−/−) show low levels of m5C (Fig 4E) [11]. However, acute depletion of NSUN2 by oxidative stress exposure resulted in no significant differences in total m5C levels in tRNAs (Fig 4E). Because tRNAs are highly abundant and stable and commonly contain m5C, we considered the possibility that early site-specific changes might not be detectable by MS. Therefore, we performed RNA BS-seq to quantify m5C at single-nucleotide resolution in tRNAs (S4D Fig) [11,40].

RNA BS-seq revealed only a modest overall reduction of cytosine-5 tRNA methylation at NSUN2-dependent sites in two independent experiments (Fig 4F; S4E Fig; S8A and S8B Data). NSUN2-independent methylated sites were unaffected (S4F Fig). Furthermore, our data allowed a high-resolution inspection of specific sites that require NSUN2 activity, which showed dynamic and reproducible changes in response to stress (Fig 4G–4I; S4G–S4I Fig). For example, methylation at C34 in LeuCAA decreased gradually with arsenite treatment (Fig 4H and 4I; S4H Fig). Methylation of AspGTC at position C47 only significantly changed at 2 hours of arsenite treatment, whereas nearby sites (C38, C48) remained unchanged (Fig 4H and 4I; S4I Fig). We observed a similar change in methylation levels when we pooled all other potentially non–tRNA targeted sites by NSUN2, whereas other m5C sites were unaffected by stress (Fig 4J; S9A and S9B Data). Together, our data revealed that NSUN2-mediated deposition of m5C at distinct sites in tRNAs changed dynamically in response to oxidative stress.

Dynamic changes of site-specific m5C levels require NSUN2

Next, we asked whether NSUN2 was solely responsible for causing the site-specific methylation changes in tRNAs in response to oxidative stress. We rescued NSUN2−/− cells with NSUN2, K190M, or the empty vector as a control and performed RNA BS-seq. We selected for all sites with a minimum coverage of 100 and more than 5% methylation in the pooled replicates (n = 5) of untreated NSUN2-overexpressing cells (approximately 7,400 sites) (S10 Data). Reexpression of NSUN2 significantly restored 525 methylation sites when compared to the empty vector control (Fig 5A) and 431 sites when compared to K190M-overexpressing cells (Fig 5B). Reexpression of NSUN2 restored the level of tRNA methylation similarly to endogenous NSUN2 (S5A Fig).

[Figure omitted. See PDF.]

Fig 5. NSUN2-mediated tRNA methylation is dynamic and site-specific.

(A, B) Volcano plot depicting the significant methylation changes when NSUN2 was reexpressed in NSUN2−/− cells compared to empty (‘E.’) vector (A) or K190M controls (B). (C) Global methylation levels of all m5C sites identified in (B) after treatment with arsenite for 0, 2, or 4 hours. Shown are all sites >20% methylation in NSUN2-rescued cells. (D, E) Mass spectrometry analyses to quantify the number of methylated sites (m5C) in all tRNA (D) or only tRNA LeuCAA (E) in response to stress. (F) Heatmap showing all significantly different m5C sites (p < 0.05) changing upon stress in NSUN2-overexpressing cells. (G-L) Examples of m5C sites identified in (F). *padj < 0.05; **padj < 0.01; ****padj < 0.0001 (ordinary one-way ANOVA, multiple comparisons). (M) Mass spectrometry analyses to quantify m5C in tRNA leucine (upper panel) and all tRNAs (lower panel) in the presence of an angiogenin inhibitor (‘Angi’) and arsenite. Data presented as mean (n = 3), error bars ± SD. p-Value: padj ANOVA. (N-P) Methylation levels (pooled from 5 replicates) of cytosines along tRNA 74-Leu CAA (upper panels) and 145-Leu CAA (lower panels) detecting all m5C sites within the tRNA molecule with different dynamic changes in response to stress. The underlying data for this figure can be found in S10 Data and S1 File. m5C, 5-methylcytosine; tRNA, transfer RNA; VL, variable loop.

https://doi.org/10.1371/journal.pbio.3000297.g005

We confirmed that the overexpressed proteins behaved like endogenous NSUN2 and were down-regulated upon exposure to arsenite (S5B Fig), and then we assessed how global methylation levels changed. We selected all significant sites with more than 20% methylation in the NSUN2-rescued cells when compared to K190M-infected cells and found that the overall median methylation levels remained largely unaltered (Fig 5C). MS for m5C using all tRNAs or tRNA LeuCAA confirmed similar methylation levels after exposure stress (Fig 5D and 5E). Thus, global methylation levels in NSUN2−/− rescued cells remained largely stable in response to stress.

Then, we identified all individual m5C sites showing significantly different methylation levels in NSUN2-rescued cells after 2 or 4 hours of stress (Fig 5F). A majority of these specific sites required NSUN2 (Fig 5F), and their changes in methylation were highly comparable to cells expressing endogenous NSUN2 (NSUN2+/+; Fig 4G–4I). For instance, C34 in the anticodon loop of tRNA LeuCAA exhibited up to 4-fold reduction of methylation (Fig 5G and 5H). Levels of m5C located to other positions showed a more modest reduction or were unaltered (Fig 5I–5L). Thus, we identified specific m5C sites that recapitulated the same changes in methylation when exposed to stress as described for NSUN2+/+ cells. Our rescue experiment therefore demonstrates that distinct m5C sites in tRNAs changed dynamically in response to stress and that these changes in methylation levels directly depended on NSUN2.

Methylation levels within the same tRNA molecule are independent from each other

To explain why the site-specific methylation levels were not detected by MS, we first considered the possibility that unmethylated tRNAs escaped the MS analysis because of the long half-life of methylated tRNAs and the enhanced biogenesis of tRFs in response to stress. To test this hypothesis, we measured m5C in stress-exposed tRNAs in the presence of an angiogenin inhibitor (N65828) [41]. Indeed, we measured a reduction of m5C in all tRNAs and more than 2-fold reduction in tRNA LeuCAA (Fig 5M), confirming that angiogenin was a major endonuclease cleaving unmethylated tRNAs (S5C Fig). However, the number of m5C per tRNA LeuCAA decreased only slightly after 4 hours of stress (Fig 5M; upper panel). tRNA LeuCAA contained two methylation sites, one in the anticodon C34 and one in the VL (Fig 5N–5P), and these sites were affected differently by oxidative stress. Only C34 decreased upon stress (Fig 5N–5P). We observed a similar dynamic change of methylation patterns within tRNA HisGTG (S5E Fig). Moreover, other leucine isotype tRNAs such as LeuCAG were not methylated at C34 (S5D Fig). Thus, our data revealed that site-specific methylation of tRNAs changes dynamically during the cellular stress response even within the same tRNA molecule.

Site-specific tRNA methylation determines tRFs biogenesis in response to oxidative stress

Next, we asked how altered expression of NSUN2 during the cellular stress response affected the biogenesis of tRFs. We performed small RNA-seq to identify all tRNA-derived sequences that significantly (padj < 0.01) changed in NSUN2−/− cells after exposure to sodium arsenite for 2 hours (S11A and S11B Data). As expected, most tRNA-derived sequences started from the first tRNA nucleotide. Furthermore, we observed an enrichment of fragments starting at around positions 20 and 36, corresponding to the D-loop and the anticodon loop, respectively [42,43] (Fig 6A). Because some RNA modifications can stall reverse transcription (RT) and potentially cause sequencing biases, we only analysed tRNA-derived sequences with consistent coverage in all four replicates and omitted all fragments with missing read values in any of the conditions from the analyses (S11C Data).

[Figure omitted. See PDF.]

Fig 6. Site-specific tRNA methylation determines biogenesis of distinct tRFs.

(A) Density plot of tRNA-derived sequences beginning at the indicated positions. (B) Clustering of tRNA-derived fragments < 40 nucleotides in NSUN2+/+ and NSUN2−/− cells untreated (‘Ctr’) or treated with sodium arsenite for 2 or 4 hours. Shown are 3 out of 4 replicates per time point. (C, D) Heatmap (C) and log2 FC of tRFs shown in (C). (E, F) PCAs (E) and heatmap (F) of significantly different tRFs in NSUN2−/− cells rescued with a wild-type (NSUN2) or point mutated (K190M) NSUN2 construct after 4 hours of exposure to stress compared to the untreated control (‘0h’). (G, H) Violin plots showing the read distribution of the tRFs shown in (E, F). (I) Log2 FC of the up-regulated tRFs when NSUN2-overexpressing cells are exposed to stress for 4 hours. tRNA glutamic acid–derived tRFs are highlighted in red. Line indicates the mean. The underlying data for this figure can be found in S11 and S12 Data and S1 File. FC, fold-change; PCA, principle component analysis; tRF, tRNA-derived fragment; tRNA, transfer RNA.

https://doi.org/10.1371/journal.pbio.3000297.g006

Next, we focused specifically on tRFs that were smaller than 40 nucleotides, which capture the products of endonucleolytic cleavage (S11D Data). Principal component analysis (PCA) confirmed a clear separation by both genotype and treatment (Fig 6B), and tRFs from NSUN2−/− cells clustered closer to NSUN2+/+ cells exposed to stress for 4 hours (Fig 6B). Furthermore, tRFs clearly separated the 2- and 4-hours-treated NSUN2−/− cells from their untreated control. In contrast, NSUN2+/+ cells only showed distinct changes in tRFs when exposed for 4 hours to arsenite (Fig 6B). When subjected to stress, we discovered three distinct patterns of tRF production (Fig 6C; Cluster 1–3). Each cluster was dominated by distinct tRNA isoacceptors (Fig 6D), indicating that tRNAs were subjected to processing based on their sequence and potentially other corresponding internal modifications. Our data support a model in which all cells produce tRFs but the tRNA fragmentation pattern depends on the identity of the isoacceptor and the presence of modifications at distinct sites within the tRNA molecule.

In line with the reduction of NSUN2 protein expression in response to stress, tRFs derived from NSUN2-methylated tRNAs significantly increased after 4 hours of stress (S6A–S6C Fig; red and S5C Fig). The same tRFs remained unchanged or decreased in NSUN2-depleted cells (S6A–S6C Fig; blue and S5C Fig). In conclusion, the production of tRFs in response to stress differed substantially in the absence of NSUN2, confirming that tRF biogenesis was influenced by tRNA methylation, isoacceptor identity, and cellular stress.

Biogenesis of tRF subsets is directly determined by tRNA methylation

To identify tRFs whose biogenesis directly required the enzymatic activity of NSUN2, we performed small RNA-seq analyses in the rescued NSUN2−/− cells that were either untreated or exposed to sodium arsenite. We identified all tRFs with differential abundance in NSUN2-rescued cells after 4 hours of stress (Fig 6E; S12 Data). Untreated NSUN2-rescued cells clustered close to K190M-overexpressing cells (Fig 6E), which was due to a distinct set of up-regulated and down-regulated tRFs (Fig 6F; ‘NSUN2 4h’). The production of these tRFs required NSUN2-mediated methylation because they remained unchanged in K190M-overexpressing cells (Fig 6G and 6H). As described for cells expressing endogenous levels of NSUN2 (NSUN2+/+), the up- and down-regulated groups of tRFs were dominated by specific tRNA isoacceptors (Fig 6I; S6D Fig). We concluded that NSUN2-mediated methylation at distinct sites within the tRNA molecule protected from processing into tRFs and thereby promoted efficient mRNA translation.

In particular, the formation of tRNA GluCTC-derived tRFs depended on NSUN2 (Fig 6I). tRFs inhibit protein synthesis via several mechanisms including through direct inhibition of the ribosome or displacement of RNA-binding proteins [13,44,45]. To demonstrate that tRFs directly link stress sensing by NSUN2 to repression of protein synthesis, we independently added synthetic 5′ and 3′ tRNA GluCTC fragments to NSUN2+/+ and −/− cells and measured protein synthesis (S6E–S6H Fig). A scrambled noncoding RNA served as a negative control. Within 2 hours, 5′ and 3′ derived tRFs reduced protein synthesis in NSUN2+/+ cells (S6F Fig). Similar to the stress response, the protein synthesis rate was only temporarily repressed (S6F Fig). In contrast, protein synthesis in NSUN2−/− cells was unaffected by the introduction of additional tRFs, since NSUN2−/− cells are already saturated (S6G Fig; S6C Fig; S5C Fig). We confirmed that the transfection efficiency was comparable in NSUN2+/+ and −/− cells (S6H Fig). Thus, our data demonstrate that the intrinsic formation of tRFs is sufficient to dynamically regulate protein synthesis.

m5C is required to balance anabolic and catabolic pathways during the stress response

To better understand how the failure to adapt protein synthesis rates to the external stress stimulus affected metabolic pathways, we identified all differentially expressed genes in NSUN2−/− cells after exposure to sodium arsenite (Fig 7A; S13A–S13C Data). A total of 2,799 genes were differentially expressed in NSUN2−/− cells when compared with NSUN2+/+ cells in all conditions and thus not influenced by the stress signal. We identified 884 and 1,584 misregulated genes in NSUN2−/− cells after 2 and 4 hours of stress, respectively (Fig 7A).

[Figure omitted. See PDF.]

Fig 7. Loss of NSUN2 alters mitochondrial function and catabolic pathways in response to stress.

(A) Venn diagram showing all significantly expressed genes in NSUN2−/− cells compared to NSUN2+/+ when untreated (‘Ctr’) or treated for 2 and 4 hours with sodium arsenite. (B, C) Gene enrichment analysis for biological processes (GOrilla) using the 884 uniquely changed genes in NSUN2−/− cells after 2 hours (B) or the 1,584 uniquely changed genes after 4 hours (C) of stress exposure. Colour code indicates p-value, and size reflects enrichment. The underlying data for this figure can be found in S13 and S14 Data.

https://doi.org/10.1371/journal.pbio.3000297.g007

Two hours after stress exposure, regulators controlling mitochondrial function were significantly enriched in the differentially expressed genes (Fig 7B; S14A Data). In particular, genes encoding for complex I subunit assembly significantly changed in the absence of NSUN2. Complex I is the major entry point for electrons into the respiratory chain and likely to act as the rate-limiting step in respiration [46]. In line with these results, we confirmed that NSUN2−/− cells contained less-active mitochondria in particular during the early stress response (S7A and S7B Fig).

Four hours after stress exposure, the differentially expressed genes were significantly enriched in regulators of RNA and protein catabolic pathways (Fig 7C; S14B Data). These results closely parallel those of our metabolic assay that showed how loss of m5C captures cells in a catabolic state and affected mitochondrial function (Fig 1).

Finally, we asked whether factors promoting the stable catabolic state in NSUN2−/− cells were translationally regulated and could therefore be identified using ribosome profiling (Ribo seq). We used our previous published dataset (phs000645.v5.p1), in which we profiled NSUN2−/− cells that were rescued with NSUN2 or two different enzymatic dead versions of the protein (C271A, C321A) and compared them to the empty vector control cells (S15 Data) [12]. We confirmed that genes specifically enriched in the NSUN2-rescued cells were regulators of mRNA translation and RNA catabolic processes (S7C–S7E Fig). We now have further demonstrated that the formation of m5C balances protein synthesis with the metabolic requirements of stress responses.

In summary, our data revealed a highly dynamic regulation of protein synthesis rates in response to stress that is tightly coordinated by tRNA methylation and cleavage. We further discovered that loss of a major RNA methyltransferase shaped the biogenesis of tRNA fragments and thereby induced a catabolic cell state.

Discussion

mRNA translation is a critical step for all gene expression programmes and also represents the most energy-consuming processes within cells [47]. Therefore, global and transcript-specific translation are continuously adapted to environmental cues including nutrients, growth factors, and stress stimuli. How precisely the external cues are sensed by the mRNA translation machinery, integrated into metabolic pathways, and then coordinated to trigger the appropriate cell response is less well understood. Here, we identify NSUN2 as an important sensor for oxidative stress. NSUN2 further links the environmental cue to the protein synthesis machinery via tRNA methylation. Loss of NSUN2 repressed global protein synthesis and thereby induced a catabolic cell state without affecting gene transcription.

External stress stimuli repress global protein synthesis, causing a switch of translation to stress- and cell type–specific regulatory proteins [48]. tRNAs play multiple regulatory roles in the adaptation of protein synthesis to the cellular stress response. For instance, the content of internal tRNA modifications and also the level of charged tRNAs rapidly changes in response to stress [49]. Loss of tRNA methylation at the VL causes global reduction of protein synthesis [11,12]. Changing the modifications at the wobble anticodon position alters transcript-specific translation [50–52].

Mature mammalian tRNAs are extremely stable, with an estimated half-life of 2–4 days [49,53], making the detection of dynamically changing m5C sites within hours of a stress response challenging. Nevertheless, we detected site-specific changes by RNA BS-seq. How m5C was site-specifically removed in distinct tRNAs is unknown, as no eraser protein has been identified so far. Alternatively, unmodified tRNAs might have a higher turnover and are therefore more difficult to capture than fully modified tRNAs.

Nucleotide modifications outside the anticodon loop are often linked to differential processing and cleavage of tRNAs [10,11,54,55]. Accordingly, NSUN2-expressing and -lacking cell populations differed in their production of tRFs under normal conditions and when exposed to cell stress. Depending on NSUN2-expression and stress, we observed distinct tRNA processing patterns. For instance, tRFs derived from tRNA GluCTC were specifically enriched in NSUN2-depleted cells and decreased in abundance with stress. Similarly, synthetic tRNA GluCTC fragments efficiently repressed protein synthesis, yet the effect was only short-lived. How the half-life of tRFs is regulated is currently unknown, but our data strongly indicate a dynamic turnover of the tRFs in response to stress. A high turnover of distinct tRFs can also explain the only temporal repression of global protein synthesis after stress exposure.

In summary, NSUN2-mediated formation of m5C is critical for integrating the cellular metabolic state with global protein synthesis and thereby triggers the appropriate cellular responses to external stress stimuli.

Methods

Ethics statement

The research including mice has been regulated under the Animals (Scientific Procedures) Act 1986 Amendment Regulations 2012 following ethical review and approval by the University of Cambridge Animal Welfare and Ethical Review Body (AWERB) under the terms of the United Kingdom Home Office licences PPL80/2231, PPL80/2619 and PPL_P36B3A804.

Transgenic mice

The NSUN2 knockout mice (homozygous Nsun2Gt[D014D11]Wrst) were generated and genotyped as previously described [18].

Cell sorting and analysis

Mouse keratinocytes from skin in telogen at postnatal day (P)49 of wild-type or NSUN2−/− male mice were isolated as follows: whole mouse back skin was sterilised with 10% Betadine and 70% ethanol and washed in phosphate-buffered saline (PBS). The dermal side was floated on 0.25% trypsin without EDTA (Thermo Fisher Scientific) for 2 hours at 37°C. The epidermis was subsequently scraped from the dermis and disaggregated by gentle pipetting in low-calcium medium with 10% FCS and filtered through a 70 μm cell strainer. The cells were pelleted and resuspended in the following antibodies to cell surface markers in 2% bovine serum albumin (BSA): PE-conjugated anti-ITGA6 (1:500, clone GoH3, eBiosciences), Alexa Fluor 647-conjugated anti-CD34 (1:50, RAM34, eBiosciences), and goat anti-PCAD (1:50, R&D Systems). After incubation for 30 minutes at 4°C, cells were washed twice in PBS. For detection of PCAD, cells were incubated for 10 minutes at 4°C with anti-goat Alexa Fluor 488-congugated secondary antibody (1:500, Thermo Fisher Scientific). Cells were gated using forward versus side scatter to eliminate debris. Doublet discrimination was carried out using pulse width. The viable cells were then gated by their exclusion of sytox. Cells were sorted with a MoFlo high-speed sorter (Beckman Coulter) as follows: ITGA6high/CD34+ as hair follicle BG stem cells and ITGA6low/P-CADhigh as HG progenitor cells.

Amplification of RNA from flow-sorted cells

Whole back skin of wild-type or NSUN2−/− male mice at P49 was collected as skin in telogen and at P32–34 as skin in anagen. Hair follicle stages in telogen and anagen were verified by HE staining for each skin biopsy. RNA from flow-sorted cells was purified using Pure-Link RNA Micro Isolation Kit (Thermo Fisher Scientific). RNA was purified from 105 to 3 × 105 ITGA6high/CD34+ per biological replicate and 104 to 5 × 104 ITGA6low/PCADhigh per biological replicate. For each biological replicate, flow-sorted cells from 1 to 3 mice were pooled. Total RNA (75 pg; with an RIN of 8 or above) were amplified using an adapted version of the Kurimoto protocol [56]. The main changes to their published protocol included the use of an increased amount of Superscript III, primers, and dNTPs as well as a longer RT reaction. Samples were amplified for a total of 29 cycles, and each sample was split into 4 during this amplification process and combined at the end to avoid errors introduced by PCR being overrepresented in the final library. Low-DNA/RNA bind tips and tubes were used throughout the experiment.

Gene expression arrays and analyses

RNA samples were amplified using the Genechip WT Plus kit (Thermo Fisher Scientific). Briefly, the RNA was converted into cDNA and amplified using a mix of primers to target polyA and non-polyA mRNAs generating biotin-labelled cRNA. The cRNA was then hybridised to the Gene ST array, stained, and scanned using the GeneTitan Instrument (Thermo Fisher Scientific). After purification, quality control, and quantity normalisation, the cRNAs of four samples per genotype and condition were hybridised to the Affymetrix MouseWG-6 v2.0 Expression BeadChip (Illumina). Hybridisation, washing, staining, and scanning were performed according to standard Illlumina protocols (Illumina Whole-Genome Gene Expression DirectHyb Assay). Microarray hybridisation, washing, and scanning were performed at the Cambridge Genomic Services at the Department of Pathology (University of Cambridge, Cambridge, UK). Differential expression analysis was done with the R limma package using default settings.

In situ hybridisation via RNAscope

Nsun2 RNA in mouse skin was labelled using the RNAscope in situ hybridisation technology [57]. Freshly cut paraffin-embedded tissue sections were first heated at 56°C for 1 hour, dewaxed, and dehydrated. Endogenous peroxidases were blocked, followed by target retrieval steps (antigen retrieval and protease digestion). The Nsun2 probe (Cat.426721) was hybridised for 2 hours at 40°C, followed by the six-step amplification protocol. The signal was developed with DAB, and slides were counterstained with 6% Mayer’s Haematoxylin, further dehydrated, cleared in Xylene, and mounted in DPX.

Cell lines and culture conditions

NSUN2+/+ human dermal fibroblasts were purchased from Invitrogen (C0135C; Thermo Fisher). NSUN2+/− and two lines of NSUN2−/− human dermal fibroblasts were established as previously described [7]. All cells were cultured in Minimum Essential Medium (MEM, 31095–052; Thermo Fisher) supplemented with 20% HyClone FBS (SV30180.03; GE Healthcare) and 1% Penicillin-Streptomycin (P0781; Sigma-Aldrich). For the rescue experiments, NSUN2−/− cells were stably infected with the wild-type or an enzymatic dead version of NSUN2 carrying a single point mutation (K190M) [30]. Empty vector (pBABE)–infected cells served as a control. All NSUN2−/− rescue cells were cultured in MEM (31095–052; Thermo Fisher) supplemented with 20% HyClone FBS (SV30180.03; GE Healthcare) and no antibiotics.

Cells were passaged using Trypsin-EDTA (0.25%, 25200072; Thermo Fisher) diluted in PBS to a 1:1 ratio and split as required. Cells were maintained at 37°C in a humidified incubator with a 5% CO2 atmosphere. Cells were grown on plastic dishes or flasks of tissue culture grade, depending on the experiment (Falcon; Corning and Nunc; Thermo Fisher).

To induce oxidative stress, cells at 80% confluency were incubated with fresh prewarmed media containing 200 μM sodium arsenite (stock solution 200 mM in PBS; NaAsO2, Sodium [meta]arsenite, S7400; Sigma-Aldrich) for the indicated time. To induce UV radiation stress, cells at 80%–90% confluency without medium were exposed to 100 J/m2 of UV light in a CL-1000 Ultraviolet Crosslinker (UVP). Fresh media were added directly after the exposure. To label the active mitochondria, the MitoTracker Red CMXRos reagent was used (M7512; Invitrogen) following the manual’s recommendations. Cells were incubated with 200 mM MitoTracker reagent (in DMSO) dissolved in cell culture medium for 30 minutes. To inhibit the mTOR pathway, rapamycin was used at a 500 nm or 1 μM concentration diluted in cell culture medium (stock solution: 2.74 mM in DMSO, R8781; Sigma-Aldrich). To inhibit angiogenin, cells were exposed to 43 μM of the small-molecule inhibitor N65828 (8-amino-5-[40-hydroxybiphenyl-4-yl azo] naphthalene-2-sulphonate) obtained from the National Cancer Institute (http://dtp.cancer.gov) (stock aliquots of 1 mg/mL: 1 mg powder in 100 μl DMSO and 900 μl PBS). Cells incubated with 0.5% (v/v) PBS or 0.5% (v/v) DMSO in cell culture medium served as controls.

For transfection with tRNA fragments, cells were grown in 6-well plates and additionally transfected with synthetic 5′ and 3′ tRNA fragments. For this, 10 μM of tRNA fragments were mixed with antibiotic-free medium and the DharmaFECT1 Transfection Reagent, by following manufacturer's recommendations (T-2001, Dharmacon; GE Healthcare). The tRNA fragments used were 5′-Glu-CTC tRNA-F synthetic sense (5′ UCC CUG GUG GUC UAG UGG UUA GGA UUC GGC GCU CUC) and 3′-Glu-CTC-tRNA-F synthetic sense (5′ CCG CCG CGG CCC GGG UUC GAU UCC CGG UCA GGG AA) (Thermo Fisher). Mock control samples with just the transfection reagent and 10 μM of fluorescently labelled control siRNA (Qiagen) served as negative controls. Samples were collected in a time course for up to 24 hours posttransfection.

Metabolic analysis

For LC-MS/MS, NSUN2+/+, NSUN2+/−, and two biological replicates of NSUN2−/− cells were cultured in 150 mm round dishes and in five technical replicates each. In total, 8 to 10 million cells were collected per sample, pelleted, and flash-frozen in liquid nitrogen until further analysis. In addition, back skin from NSUN2−/− was collected from P26–P27 male mice in anagen. A total of 5 wild-type and 3 NSUN2−/− mice were used. Fat and connective tissue was scraped off, and the skin samples were snap frozen in liquid nitrogen. Each sample (70 mg) was homogenised on ice in methanol/chloroform using a tissue homogeniser (Polytron PT2500, Kinematica). To assess the metabolic content of these samples, LC-MS/MS was performed and optimised with internal standards, for the targeted analysis of aqueous metabolites including methionine and TCA cycle intermediates, nucleotides, and amino acids [58,59]. For NMR spectroscopy, skin samples from 4 wild-type and 3 NSUN2−/− mice were prepared as described before and processed as previously described [60].

RNA isolation, RT, and RT-qPCR

Total RNA was prepared using Trizol reagent (Thermo Fisher Scientific) according to manufacturer's instructions. For further purification, total RNA was subjected to 30 minutes TURBO DNase treatment (AM2239; Invitrogen) at 37°C following manufacturer’s instructions. To inactivate the reaction, the samples were Phenol:Chloroform (77617; Sigma-Aldrich) extracted and precipitated using sodium acetate (pH 5.5) and 1–2 μl Glycoblue (AM9516; Ambion) overnight at −80°C. The RNA was washed with ethanol and resuspended. The concentration of each sample was assessed using a Qubit Fluorometer and the Qubit RNA HS Assay Kit (Q32855; Invitrogen) following the kit's instructions.

Double-stranded cDNA was synthesised from 1 μg of RNA and the Superscript III reverse transcriptase (18080085; Life Technologies) following manufacturer’s instructions using Random Primers (C1181; Promega). Each RT quantitative PCR (RT-qPCR) reaction was set up in MicroAmp Optical 96-well plates (N8010560) using 1 μl cDNA, 5 μl TaqMan Fast Universal PCR Master Mix (4366073), 3.5 μl RNAse-free water, and 0.5 μl of a predesigned Nsun2 probe (Hs00214829_m1; Applied Biosystems). A human GAPDH probe was used for normalisation using the ΔCt method (4333764T; Applied Biosystems). RT-qPCR and data acquisition were conducted using the StepOnePlus Real-Time PCR System (Applied Biosystems).

Protein extraction and western blotting

Cells were first rinsed with PBS and lysed in ice-cold RIPA buffer (50 mM Tris-HCL [pH 7.4], 1% NP-40, 150 mM NaCl, 0.1% SDS, 0.5% Sodium deoxycholate– 1 mL RIPA per T75 or 100 mm culture dish). RIPA was supplemented with cOmplete Mini EDTA-free Protease Inhibitor Cocktail tablets (11836170001, Roche), and cells were collected using a cell scraper. The lysates were centrifuged for 15 minutes at maximum speed in a precooled centrifuge at 4°C and their supernatant collected and kept on ice. The concentration of each protein sample was assessed using the Pierce BCA Protein Assay kit (23225; Thermo Fisher) according to the manufacturer's instructions and measured using a spectrophotometer and the Softmax software for protein quantification.

Cell protein lysates were mixed with NuPAGE LDS Sample Buffer (4X) (NP0007; Invitrogen) and run on 7.5% or 12.5% separating polyacrylamide gels. PageRuler Prestained Protein Ladder (26616; Thermo Fisher) was loaded as a size marker. All gels were run in 1x SDS Running Buffer (30 g Tris Base, 144 g Glycine, 10 g SDS, and H2O to 1 L; for 10x buffer) at 160V for approximately 1 hour, after which the proteins were transferred to a nitrocellulose or PVDF membrane (GE Healthcare) in 1x Transfer Buffer (30.3 g Tris Base, 144 g Glycine, and H2O to 1 L; for 10x buffer) containing 15% (v/v) methanol at 90 V, for 1.5 hours on ice. Membranes were blocked for a minimum of 1 hour at room temperature in 5% (w/v) nonfat milk or 5% (w/v) BSA (A4503-50G; Sigma-Aldrich) in 1x TBS and 0.1% Tween-20 (TBS-T) (48.4 g Tris Base, 160 g NaCl, and H2O to 1 L; for 20x TBS buffer at adjusted pH 7.6) and incubated with primary antibody in blocking solution overnight at 4°C. Each membrane was washed three times for 10 minutes in TBS-T prior to incubation with the appropriate horseradish peroxidase (HRP)-labelled secondary antibody (1:10,000) in TBS-T at room temperature for 1 hour (anti-mouse NA931 and anti-rabbit IgG HRP NA934; GE Healthcare–anti-goat IgG HRP sc-2020; Santa Cruz). Finally, the membranes were washed as before and the antibodies detected by using the Amersham ECL Prime Western Blotting Detection Reagent (RPN2232; GE Healthcare). The primary antibodies used were NSUN2 (1:1,000, 20854-1-AP, Proteintech), NSUN3 (1:500, GTX46175, GenTex), NSUN4 (1:500, ab101625, Abcam), NSUN6 (1:500, 17240-1-AP, Proteintech), NPMI (1:500, B0556, Sigma-Aldrich), and HSP90 (1:500, sc-13119, Santa Cruz).

Immunofluorescence

Cells were grown directly on autoclaved glass coverslips placed in 12-well culture dishes. Cell culture and sodium arsenite or recombinant angiogenin treatments were carried out as described before. To fix the cells, the glass coverslips were first briefly rinsed with wash buffer comprising 0.1% Tween-20 in 1x PBS and incubated for 10 minutes at room temperature with paraformaldehyde (4% PFA in PBS; Santa Cruz). Cells were then washed three times with ice-cold wash buffer for a couple of minutes. To examine intracellular proteins, cells were permeabilised for 10 minutes with PBS containing 0.25% Triton X-100 at room temperature and washed again three times for 5 minutes in wash buffer. This step was omitted for the study of membrane proteins. To block nonspecific antibody binding, cells were incubated with blocking buffer comprising 1% BSA and 22.52 mg/mL glycine in PBS with 0.1% Tween-20 (PBST) for 1 hour. To detect specific proteins of interest, cells were then incubated with primary antibodies diluted in 1% BSA in PBST at 4°C overnight. The cells were then washed three times in wash buffer for 5 minutes each. To label the detected proteins, cells were incubated with the respective Alexa Fluor secondary antibodies diluted in 1% BSA in PBST for 1 hour at room temperature on a shaker and protected from light (1:1,000; Life Technologies). Cells were washed as before and their nuclei counterstained with DAPI (1:3,000 in PBS; Sigma) for 5 minutes. Finally, cells were rinsed with wash buffer and the glass coverslips mounted on slides with PBS:glycerol solution (1:1). The primary antibodies used were Angiogenin (1:500, PC317L, Calbiochem), eIF4A1 (1:200, sc-14211, Santa Cruz), and p-eIF2a (1:1,000, 9721S, Cell Signaling).

Detection of m5C by MS

Human dermal fibroblasts were grown and treated with sodium arsenite as described above. NSUN2+/+ and NSUN2−/− cells were cultured in 100 mm dishes. Three replicates per condition and genotype were used for the experiments. Cell pellets were collected and lysed directly in Trizol for RNA extraction. To measure the m5C modifications, optimised MS analysis was performed [61]. Total RNA was loaded on a size-exclusion column (Agilent Bio SEC-3, 3 μm, 300 Å, 7.8 × 300 mm, Agilent, Waldbronn, Germany) and RNA fractions eluted with 100 mM ammonium acetate at pH 7 as the mobile phase. tRNA was separated from rRNA and small RNAs, vacuum concentrated, and reconstituted in water. tRNA concentrations were determined by UV spectroscopy at 260 nm. rRNA and other large RNA species were excluded from the fractions.

Purified tRNA (approximately 400 ng) was digested using a mixture of benzonase (2 U), bacterial alkaline phosphatase (2 U), and phosphodiesterase I (0.2 U) in a final reaction volume of 20 μL. The reaction mixture was supplemented with MgCl2 to a final concentration of 1 mM and Tris-HCl (pH 8.0) to a final concentration of 50 mM. Nucleobase deaminase inhibitor coformycin and tetrahydrouridine were added at a final concentration of 10 μg/mL and 50 μg/ml, respectively, and butylated hydroxytoluene (an antioxidant) was added at a final concentration of 0.5 mM (for further detail, see [62]). The digestion was allowed to proceed for 2 hours at 37°C and was stopped by filtering through a 10 kDa MWCO filter (AcroPrep Advance, 350 μl, Omega 10K MWCO, Pall, Dreieich, Germany) at 3,000g for 30 minutes. After addition of 10 μl pure water for salt dilution purposes, 18 μL of filtrate was mixed with 2 μL of internal standard (produced as recently described [63]). Each sample (10 μL) was injected for LC-MS/MS analysis (corresponding to around 150 ng tRNA digest). Calibration solutions for absolute quantification were prepared as recently described [63].

For quantification, an Agilent 1290 Infinity II equipped with a DAD combined with an Agilent Technologies G6470A Triple Quad system and electro-spray ionisation (ESI-MS, Agilent Jetstream) was used. Operating settings were as follows: positive ion mode, skimmer voltage 15 V, cell accelerator voltage 5 V, N2 gas temperature 230°C and N2 gas flow 6 L/minute, sheath gas (N2) temperature 400°C with a flow of 12 L/minute, capillary voltage of 2,500 V, nozzle voltage of 0 V, and the nebuliser at 40 psi. The instrument was operated in dynamic MRM mode, and the individual MS parameters for the nucleosides are summarised in S16 Data. The mobile phases were A as 5 mM NH4OAc (≥99%, HiPerSolv CHROMANORM, VWR) aqueous buffer, brought to pH 5.6 with glacial acetic acid (≥99%, HiPerSolv CHROMANORM, VWR), and B as pure acetonitrile (Roth, LC-MS grade, purity ≥ 99.95). A Synergi Fusion-RP column (Phenomenex, Torrance, CA, United States; Synergi 2.5 μm Fusion-RP 100 Å, 150 × 2.0 mm) at 35°C and a flow rate of 0.35 ml/minute were used. The gradient began with 100% A for 1 minute and increased to 10% B by 5 minutes and to 40% B by 7 minutes. The column was flushed with 40% B for 1 minute and returned to starting conditions to 100% A by 8.5 minutes followed by re-equilibration at 100% A for 2.5 additional minutes.

For the isolation of single tRNA isoacceptors, a modified approach was used to the one previously described in [64]. Briefly, tRNA LeuCAA was hybridised with a reverse complementary, biotinylated DNA-oligonucleotide (AGT CTG GCG CCT TAG ACC ACT CGG CCA TCC TGA CAA A [Biotin]) followed by immobilisation on streptavidin-coated magnetic beads (Dynabeads MyOne Streptavidin T1, Life Technologies). The hybridisation step was performed in 5× SSC buffer (20x: 3 M NaCl, 300 mM trisodium citrate [pH 7.0]) using 1 μl biotinylated oligonucleotide (100 μM) and 1 μg of SEC purified total tRNA per 25 μl beads. Samples were denatured at 90°C for 3 minutes, hybridised at 65°C for 10 minutes, and cooled to room temperature. The magnetic Dynabeads were washed as described in [64]. Immobilisation of the hybrid was performed at 25°C under shaking for 30 minutes. The supernatant containing nontarget tRNAs was removed, and the beads were washed once in 1× SSC buffer and three times in 0.1× SSC buffer. Finally, the beads were resuspended in 20 μl MilliQ water and heated to 80°C for 3 minutes to elute the target tRNA. No further DNase treatment was carried out. The isolated tRNA was digested with 10 μl digestion mix (final volume: 30 μL). Finally, 9 μl of the digested and isolated tRNA was coinjected with 1 μl stable isotope-labelled internal standards for measurement by MS [63].

Sucrose density centrifugation

Sucrose gradients were used to separate subpolysomal and polysomal ribosomes. Sucrose gradients (10%–50% [w/v]) were prepared in gradient buffer (300 mM NaCl2, 15 mM MgCl2, 15 mM Tris-HCL [pH 7.5], 1 mM DTT, 0.1 mg/ml CHX).

Human dermal fibroblasts were grown and treated with sodium arsenite as described above. At 70%–80% confluency, cells were washed twice in PBS-CHX (100 μg/ml) and scraped into 500 μl lysis buffer (300 mM NaCl2, 15 mM MgCl2, 15 mM tris-HCL [pH 7.5], 1 mM DTT, 0.2 M sucrose, 0.1 mg/ml CHX, 0.5% IGEPAL, 100 U RNasin Plus RNase Inhibitor [N2615; Promega] per 0.5 ml). Cell lysates were incubated on ice for 3 minutes prior to pelleting the cells at 1,300g for 5 minutes. The supernatant was then layered on top of the gradient and centrifuged at 38,000 rpm (acceleration 9, deceleration 6) for 2 hours at 4°C using an ultracentrifuge (Beckman Coulter). Gradients were fractioned using a gradient fractionation system (Presearch), and fractions were collected at 1-minute intervals using the FOXY Jr collection system (Presearch) at a flow rate of 1 ml/minute. Absorbance was measured constantly at 254 nm using a UA-6 UV-VIS detector (Presearch).

Sequencing library preparation

For all sequencing experiments, total RNA was isolated from human dermal fibroblasts, as described. The quality of the RNA and prepared libraries was assessed using the appropriate high-sensitivity RNA and DNA microfluidics chip on the Agilent Bioanalyzer platform (Agilent Technologies). All cDNA libraries were subjected to a cleanup step with the use of Agencourt AMPure XP beads (Beckman Coulter) and further multiplexed for sequencing on a HiSeq platform (Illumina).

RNA BS-seq

RNA BS-seq libraries were generated from untreated and sodium arsenite–treated human dermal fibroblasts for 2 or 4 hours. Two independent experiments were performed; the first consisted of 4 and the second of 5 replicates per genotype and time point. Thus, a total of 9 replicates per condition and genotype were processed and analysed. To generate the samples, RNA was prepared and BS-treated as previously described [11]. Briefly, 15 μg of total RNA was depleted from the ribosomal RNA population by using the Ribo-zero kit (MRZH11124, Illumina). Approximately 1–2 μg of rRNA-depleted samples were then BS-treated by mixing them with 42.5 μl 40% sodium BS solution (pH 5.0) and 17.5 μl DNA protection buffer supplied with the EpiTect Bisulfite Kit (59104; Qiagen). The reaction mixture was then incubated for three to four cycles of 5 minutes at 70°C, followed by 1 hour at 60°C on a thermal cycler. To desalt the reaction, all samples were passaged through Micro Bio-Spin 6 chromatography columns, following the manufacturer's instructions (732–6221; Bio-Rad), and next desulfonated by adding an equal volume of 1 M Tris (pH 9.0) to the reaction mixture and incubating for 1 hour at 37°C. BS-treated RNA samples were then precipitated overnight with 2.5 volumes of 100% ethanol, 0.1 volumes of 3 M sodium acetate (pH 5.5), and 1–2 μl Glycoblue (AM9516; Ambion). To repair the 2′,3′-cyclic phosphate and 5'-hydroxyl termini produced during the BS/desulfonation reaction, T4 Polynucleotide Kinase (PNK) was used by mixing the RNA with 10x T4 PNK reaction buffer, 10 mM ATP, 10 U of T4 PNK enzyme, and RNase-free H2O to a final volume of 50 μl (M0201S; New England Biolabs). This was incubated at 37°C for 30 minutes, followed by a heat inactivation of the enzyme at 65°C for 20 minutes and an overnight precipitation of the RNA, as before. Approximately 400 ng of BS-converted RNA was finally used to generate BS-seq libraries. Taking into account that the above BS treatment and desulfonation protocol cleaves long RNAs into 100-nucleotide fragments, the TruSeq Small RNA preparation kit was used (Illumina) to generate platform-compatible libraries. In brief, kit-provided RNA-seq adapters were ligated to the BS-converted RNAs, reverse transcribed at 50°C for 1 hour with SuperScript III and 2 mM of each dNTP (18080085; Thermo Fisher Scientific), and followed by a PCR amplification programme, following the user manual's recommendations, to finalise the process.

Standard and small RNA-seq

Standard RNA-seq libraries were generated from human dermal fibroblasts previously treated with sodium arsenite for 2 or 4 hours. All cell samples were processed in four replicates per condition. Total RNA was extracted as described and depleted from ribosomal RNA using the Ribo-zero rRNA Removal Kit (MRZH11124, Illumina). The rRNA-depleted RNA was then used to generate RNA-seq libraries using the NEXTflex Rapid Directional RNA-seq Kit (5138–08; Illumina) or RNAHyper with Riboerase Kit (NEB).

Small RNA libraries were generated from human dermal fibroblasts treated with sodium arsenite for 2 or 4 hours, as described. At least four replicates per condition and genotype were processed and analysed. To generate the small RNA libraries, 10 μg of total RNA was subjected to a deaminoacylation reaction by mixing the purified RNA with 0.1 M Tris-HCl (pH 9.0), 1 mM EDTA, and RNAse-free H2O to a final volume of 300 μl. This reaction was incubated for 30 minutes at 37°C to deaminoacylate mature tRNAs. To further size-select the RNA, samples were loaded on Novex TBE Urea 6% gels (EC68652BOX; Thermo Fisher) and the desired RNA fractions purified. TBE Urea gels were pre-run in 1x TBE buffer (for 5x stock: 54 g Tris base, 27.5 g boric acid, 20 mL of 0.5 M EDTA [pH 8.0]) at 200 V for 30 minutes, while each sample was mixed with 2x Denaturing Sample Buffer (LC6876; Thermo Fisher) and denatured at 80°C for 5 minutes. To avoid renaturation, all samples were kept on ice. The polyacrylamide gels were then loaded with the RNA samples and the RNA Marker Low Easy ladder (R0002; Abnova) and run at 180 V for 1 hour. Each gel was next stained with ethidium bromide in 1x TBE buffer on a shaker for a few minutes, prior to being visualised over a UV transluminator. Using a scalpel and the RNA ladder as a guide, 20–150 bp fragments were excised. The gel piece was placed in a 0.5 ml tube pierced with a needle to create perforations, fitted in a 2 ml RNAse-free tube. To shatter the gel, the tubes were centrifuged for 5 minutes at maximum speed in a precooled 4°C benchtop centrifuge, and the gel pieces were eluted in 600 μl NaCl and 1 μl RNasin Plus RNase Inhibitor (N2615; Promega) overnight while rotating in a cold room. To remove any gel remnants and salts, each gel eluate was passed through a sterile Costar Spin-X centrifuge tube filter (CLS8160; Corning). The RNA flow-through was precipitated overnight, with ethanol and sodium acetate in −80°C. To collect the size-selected RNA, all samples were centrifuged for 30 minutes at maximum speed in a precooled 4°C benchtop centrifuge, the pellets washed with 70% ethanol, air dried, and finally resuspended in 12 μl RNAse-free H2O. To prepare the ends of the fragments for adapter ligation and library preparation, the samples were first heated at 70°C for 10 minutes and then mixed with 10x T4 PNK Reaction Buffer, 20 U of T4 PNK enzyme, and RNase-free H2O to a final volume of 50 μl (M0201S; New England Biolabs). This was incubated at 37°C for 1 hour, followed by a heat inactivation of the enzyme at 65°C for 20 minutes and a final phenol:chloroform cleanup step and overnight ethanol precipitation. The final RNA pellet was resuspended in 7 μl RNAse-free H2O. The libraries were generated using the TruSeq Small RNA Preparation Kit (Illumina), for which 3′ adenylated and 5′ phosphorylated adapters, suitable for Illumina RNA-seq, were ligated to an average of 400 ng small purified RNA fractions. Samples were then reverse transcribed at 50°C for 1 hour with the SuperScript III cDNA synthesis kit (18080085; Thermo Fisher Scientific), followed by a PCR amplification programme with Phusion DNA polymerase (F530S; Thermo Fisher Scientific).

Sequencing data analyses

For all sequencing datasets, automatically generated FastQC reports (http://www.bioinformatics.babraham.ac.uk/projects/fastqc) were used for the initial assessment of the quality and basic processing of the reads. Sequencing adapters were trimmed from the 5′ and 3′ ends of the reads using cutadapt (v.1.8.1; https://pypi.python.org/pypi/cutadapt/1.8.1).

For standard RNA-seq data, adapters were removed and paired-end RNA-seq reads were aligned to the human reference genome (GRCh38/g38) using Tophat2 (v.2.0.9; options:–read-mismatch 2 –max-multihits 1 –GTF) guided by ENSEMBL gene models (release 82), while allowing for two mismatches per read and unique alignments only. In order to determine the mRNA abundance, RNA-seq read counts for the full transcript were obtained using featureCounts (v1.5.0-p1; http://subread.sourceforge.net; options: -O—minOverlap 10 -p -C -M -B -t exon -g gene_id -T 6 -a). The datasets were then normalised, and the statistical significance of differential expression was evaluated by using the R/Bioconductor DESeq2 package (https://bioconductor.org/packages/release/bioc/html/DESeq2.html).

To determine RNA methylation levels using RNA BS-seq data, a tailored approach to the one previously published was followed [11]. BS-seq reads were aligned to the human reference genome (GRCh38/hg38) by using Bismark (http://www.bioinformatics.babraham.ac.uk/projects/bismark; v.0.14.4; options: ‘–directional–n 1 –l 100’). To complete the mapping, tRNA gene predictions were obtained from GtRNAdb (http://lowelab.ucsc.edu/GtRNAdb). Heatmaps displaying either C or T in the aligned reads at each cytosine position were generated using custom PERL scripts and matrix2png (http://www.chibi.ubc.ca/matrix2png/) for visualisation. Cytosine positions on the heatmaps were reported relative to the annotated transcriptional start sites of the documented beginning of each tRNA.

For small RNA-seq data and to determine the abundance of tRNA fragments, a tailored approach was followed [11] in which first adapter-trimmed paired-end reads were mapped to the human reference genome (GRCh38/hg38) using bowtie (v.2.1.0; http://bowtie-bio.sourceforge.net/bowtie2; options:—no-mixed—no-discordant—end-to-end), while considering only reads that mapped uniquely to the genome. To annotate the tRNA-seq fragments, tRNA genes were downloaded from GtRNAdb (http://lowelab.ucsc.edu/GtRNAdb). Counts per tRNA fragments were normalised, and their respective abundances were statistically evaluated using the R/Bioconductor DESeq2 package (https://bioconductor.org/packages/release/bioc/html/DESeq2.html) and represented as log2(DESeq2-normalised counts).

For tRNA fragment analyses, four replicates per condition and genotype were used from the second biological replicate. To identify tRNA-derived sequences, all significant fragments (p < 0.01) in NSUN2−/− cells versus NSUN2+/+ cells after 4 hours of sodium arsenite treatment were selected and filtered for fragments smaller than 40 nucleotides. All rows with missing values (0) in any condition were excluded. The PCAs and heatmaps show all four or the most similar three replicates. PCA and heatmaps were generated using https://biit.cs.ut.ee/clustvis/.

Gene set enrichment analyses were done using Enrichr (http://amp.pharm.mssm.edu/Enrichr/) [65,66]. To analyse the affected biological processes and cellular pathways in response to stress, we used the RNA-seq data obtained from human dermal fibroblasts and performed gene enrichment analyses using GOrilla (http://cbl-gorilla.cs.technion.ac.il/) [67]. As running mode, we used ‘two ranked lists of genes’ using the 2,799 commonly regulated as a background list.

Cell cycle and apoptosis

For cell cycle analyses, human dermal fibroblasts untreated and treated with sodium arsenite were briefly rinsed with PBS and collected with Trypsin-EDTA (1:1 in PBS). Following two consecutive washes with cold PBS, the cells were fixed and resuspended in ice-cold 70% ethanol. All samples were kept at 4°C until further processing. The samples were centrifuged at 12,000g for 5 minutes and rinsed twice with PBS. All pellets were resuspended in 3 ml PBS containing DAPI dye (1:3,000) and incubated at room temperature for 30–45 minutes while protected from light. The fluorescence of each sample was measured on a flow cytometer at 450/50 405 nm. To analyse the data, DAPI fluorescence measurements were visualised by histogram and the cell cycle phases determined by the curves and x-axis values.

For cell viability analyses, human dermal fibroblasts untreated and treated with sodium arsenite were assayed with the FITC Annexin V kit following the manufacturer's recommendations (556419; BD Biosciences). Briefly, cells were rinsed with PBS and collected with Trypsin-EDTA (1:1 in PBS). The fluorescence of each sample was measured at 530/30 488 nm and 610/20 561 nm. To analyse the data, fluorescence measurements were visualised by scatterplot, while compensation and quadrants were set up based on control samples of unstained cells, cells only stained with FITC Annexin V (no PI), and cells stained with PI (no FITC Annexin V).

Global mRNA translation analysis

To investigate global protein synthesis, human dermal fibroblasts untreated and treated with sodium arsenite, recombinant angiogenin, or tRNA fragments were complemented with OP-puro and further labelled as previously described [16]. Reconstituted OP-puro (50 μM; 10 mM reconstituted stock [pH 6.4]; Medchem Source) was added to each sample in culture medium precisely 1 hour prior to collection. In each assay, a sample not treated with OP-puro served as a negative control, and as a positive control, samples treated with 50 μg/mL CHX (100 mg/ml in DMSO, C4859; Sigma-Aldrich) for 15 minutes were used. Cells were rinsed with PBS and then collected with Trypsin-EDTA (1:1 in PBS). After pelleting the collected cells and washing them once with PBS, these were resuspended in 0.5 ml PFA (1% w/v in PBS; Santa Cruz) and kept for 15 minutes on ice in the dark. Following fixation, samples were washed in PBS and permeabilised in PBS supplemented with 3% FBS and 0.1% Saponin (47036; Sigma-Aldrich) for 5 minutes at room temperature. To conjugate OP-puro to a fluorochrome, an azide-alkyne cycloaddition was performed for 30 minutes at room temperature in the dark. For this, the Click-iT Cell Reaction Buffer Kit (C10269; Life Technologies) and 5 μM of Alexa Fluor 555-Azide (A20012; Life Technologies) were used. To remove excess reagents and reduce the background signal, the cells were washed twice in PBS supplemented with 3% FBS and 0.1% Saponin. Finally, all samples were resuspended in PBS containing DAPI (1:3,000) and kept at 4°C until further analysis. The fluorescence of each sample was measured at 450/50 405 nm for DAPI and 585/15 561 nm for OP-puro. To analyse the data, fluorescence measurements were visualised by histogram, and the raw fluorescence values were extracted.

To visualise polypeptide synthesis in cultured cells, human dermal fibroblasts were grown directly on autoclaved glass coverslips placed in 12-well culture plates. Cell culture, sodium arsenite treatments, and the incorporation of OP-puro were carried out as described. To fix the cells, the glass coverslips were first briefly rinsed with 1x PBS and incubated for 15 minutes at room temperature with PFA (4% PFA in PBS; Santa Cruz). To permeabilise the cells and additionally block unspecific antibody binding, all coverslips were first incubated with 0.4% Triton X-100 in PBS for 15 minutes and next with 3% BSA in PBS for 1 hour at room temperature. To conjugate OP-puro to a fluorochrome, an azide-alkyne cycloaddition was performed for 30 minutes at room temperature and in the dark. For this, the Click-iT Cell Reaction Buffer Kit (C10269; Life Technologies) and 5 μM of Alexa Fluor 555-Azide (A20012; Life Technologies) were used. To remove excess reagents and reduce the background signal, the cells were washed once with 3% BSA in PBS. Op-puro-labelled cells were counterstained for their nuclei by using Hoechst (working dilution 0.5 μg/mL) for 10 minutes and mounted on glass slides with a PBS:glycerol solution (1:1). Glass slides were stored in the dark at 4°C until imaging.

Flow cytometry data acquisition and analysis

Flow cytometry analysis of all cell assays described was performed with the LSRFortessa Flow Cytometer (BD Biosciences). Data were analysed using the FlowJo software. All samples were gated using forward versus side scatter to eliminate debris.

Image acquisition and analysis

All bright-field and immunofluorescence images were acquired using an Axio Imager 2 Upright Research Microscope (Zeiss) with ORCA-Flash4.0 digital camera (Hamamatsu) or AxioCam MRc colour digital camera (Zeiss). All images were initially processed using the ZEN Microscope and Imaging Software and later the Adobe Photoshop CS6 package. The CellProfiler software packages were used for image quantifications and analysis. Finally, all images and figures were arranged using Adobe Illustrator CS7.

Data analysis and statistical tests

Data analysis and plotting were performed with the RStudio package and the GraphPad Prism 7 software. Statistical tests were performed with GraphPad Prism 7. If not stated otherwise, data are presented as mean values, and error bars represent standard deviation (SD).

Supporting information

S1 Fig. Metabolic profiles of NSUN2-expressing and -lacking cells.

(A) Schematic representation of stem cell differentiation in the absence or presence of NSUN2. (B) Marker expression and cellular differences in hair follicles expressing or lacking NSUN2. (C) Gating used for flow cytometry sorting of BG stem cells (ITGA6high/CD34+) in telogen (P49) wild-type and NSUN2−/− mice. (D) Gating used for flow cytometry sorting of HG progenitor cells (ITGA6low/PCAD+) in telogen (P49) wild-type and NSun2−/− mice. (E,F) Log2 FC (E) and Gene Ontology categories (F) of combined differential expressed genes (FDR < 0.05) in anagen and progenitor (ITGA6low/PCAD+) populations from skin of NSUN2+/+ and −/− mice. (G) Schematic representation of NSUN-dependent methylation at cytosine-5. (H) Overview of the one-carbon metabolism network. (I-N) Metabolic differences between NSUN2+/+ and NSUN2−/− mice relating to the methionine cycle (I,L), free amino acids (J,M), and free nucleotides (K,N) measured by NMR-based (I-K) or MS-based (L-N) metabolic profiling (n = 3–5 mice). (O) Model of how protein homeostasis changes the balance between protein synthesis and degradation in NSUN+/+ (upper panel) and NSUN2−/− (lower panel) cells. The underlying data for this figure can be found in S2 Data and S1 File. BG, bulge; DP, dermal papilla; FC, fold-change; FDR, false discovery rate; HG, hair germ; IFE, interfollicular epidermis; ITGA6, integrin alpha-6; MS, mass spectrometry; NMR, nuclear magnetic resonance; PCAD, P-cadherin; SAH, S-adenosyl-homocysteine; SAM, S-adenosyl-methionine; SG, sebaceous gland.

https://doi.org/10.1371/journal.pbio.3000297.s001

(TIF)

S2 Fig. Rescue for loss of NSUN2 by reexpressing the wild-type or enzymatic dead protein.

(A, B) Differentially expressed genes in NSUN2−/− compared to NSUN2+/+ cells (A) and Nsun2 RNA levels in NSUN2+/+, +/−, and −/− cells (B) measured by RNA sequencing. (C, D) The transcriptional profile of NSUN2−/− cells overexpressing the NSUN2 protein is largely unaltered (C) although Nsun2 is highly expressed (D). Expression of the empty (‘e.’) vector served as a control. (E) Venn diagram of differentially expressed genes (padj < 0.01) in NSUN2−/− versus +/+ compared to NSUN2-rescued cells. (F) Two out of three replicates of polysome profiles using NSUN2+/+ and −/− cells. (G) Schematic representation of OP-puro incorporation in actively translating ribosomes. OP-puro mimics an amino-acyl-loaded tRNA molecule. (H) Example raw data outputs from OP-puro fluorescence analysis using a flow cytometer. CHX served as a control. (I) Protein synthesis measured by OP-puro incorporation in NSUN2+/+ and −/− cells after incubation with an angiogenin inhibitor (ANGi). (J) Western blot for NSUN2 and tubulin after incubation with 500 or 1,000 nm RAPA for 12 or 24 hours (h). (K) Quantification of protein expression shown in (J). (L) De novo protein synthesis in NSUN2+/+ and −/− after incubation with RAPA or CHX. DMSO served as a vehicle control (J-L). (M, N) Metabolic differences of NSUN2−/− cells rescued with the empty vector (‘e.v.’), K190M, or the NSUN2 protein shown as a PCA plot (M) or as Log2 FC differences of the significant different (p < 0.01 NSUN2 versus e.v.) metabolites (N). The underlying data for this figure can be found in S4 and S7 Data and S1 File. CHX, cycloheximide; OP-puro, O-propargyl-puromycin; PCA, principle component analysis; RAPA, rapamycin; tRNA, transfer RNA.

https://doi.org/10.1371/journal.pbio.3000297.s002

(TIF)

S3 Fig. NSUN2 regulates cell cycle phases and global protein synthesis during the cellular stress response.

(A) Example raw data outputs from OP-puro fluorescence analysis using a flow cytometer for human dermal fibroblasts treated with sodium arsenite. Dotted line represents the mean level of OP-puro positive control. (B) Immunofluorescence detection of OP-puro incorporation in human dermal fibroblasts. DAPI: nuclear counterstain. Scale bar: 20 μm. (C) Measurement of OP-puro fluorescence intensity in cells using microscope-acquired images. Each dot represents one cell. Data are represented as median. (D) Second replicate of polysome profiling of NSUN2+/+ and NSUN2−/− cells rescued with wt or mutated NSUN2 (K190M). The empty vector (‘e.V.’)-infected cells served as control (see Fig 3F–3I). (E) Example of raw data output from AnV and PI analysis to measure cell death. (F, G) Percentage of cells that are viable, apoptotic, or necrotic in NSUN2+/+ and NSUN2−/− cells exposed to sodium arsenite for the indicated hours (hr) (n = 3 samples per time point). (H) Summary of cell cycle distribution shown in Fig 3A–3D. Data represented as mean in (K-H). Error bars are ±SD. The underlying data for this figure can be found in S1 File. AnV, AnnexinV; OP-puro, O-propargyl-puromycin; PI, propidium iodide; wt, wild-type.

https://doi.org/10.1371/journal.pbio.3000297.s003

(TIF)

S4 Fig. RNA methylation levels change dynamically in response to oxidative stress.

(A) Immunofluorescence detection of the stress granules markers eIF4A1 (upper panels) and p-eIF2A (lower panels) in untreated (control) or sodium arsenite–treated NSUN2+/+ and NSUN2−/− cells. DAPI: nuclear counterstain. Scale, 20 μm. (B) Nsun2 RNA levels in response to UVB exposure in primary human keratinocytes and dermal fibroblasts. (C) Western blot for NSUN2 in NSUN2+/+ and −/− cells incubated with vehicle control (DMSO, PBS). (D) Experimental outline of sample collection and RNA BS sequencing. (E,F) Quantification of tRNA methylation percentage of NSUN2-dependent (E) and -independent (F) sites in a second independent experiment (n = 5 samples per time point). (G) Second independent RNA BS-seq data shown as heatmap of methylation status of individual tRNA molecules in NSUN2+/+ and NSUN2−/− cells. (H, I) Quantification of methylation changes in the tRNAs LeuCAA and AspGTC in NSUN2+/+ and NSUN2−/− cells shown in (E). Data represented as median in (E, F, H, I). One-way ANOVA adjusted p-value (H,I), **p < 0.005, ***p < 0.0005, ****p < 0.0001. The underlying data for this figure can be found in S8 Data and S1 File. BS, bisulfite; BS-seq, BS sequencing; eIF2, eukaryotic Initiation Factor 2; tRNA, transfer RNA.

https://doi.org/10.1371/journal.pbio.3000297.s004

(TIF)

S5 Fig. Dynamic levels of NSUN2, tRFs, and methylation in response to stress.

(A) Reexpression of NSUN2 but not K190M in NSUN2−/− cells restores methylation to similar levels of endogenous NSUN2 (NSUN2+/+). (B) Western blot for NSUN2 in NSUN2−/− cells infected with an empty (‘E.’) vector control, the enzymatic dead K190M, or the wild-type NSUN2 construct in cells untreated (0) or treated for 2 and 4 hours (h) with sodium arsenite. HSP90 served as a loading control. (C) Raw data (reads) for the indicated tRNAs obtained from small RNA sequencing in NSUN2-expressing (+/+; black) or -lacking (−/−; red) untreated (0h) or (h) treated 2 and 4 hours with arsenite. (D, E) Methylation levels (pooled from 5 replicates) of cytosines along tRNA 71-Leu CAG (C) and 1-His GTG (F) detecting m5C sites, in the variable loop (D) and C70 (E). The underlying data for this figure can be found in S10 Data and S1 File. HSP90, heat shock protein 90; m5C, 5-methylcytosine; tRF, tRNA-derived fragment; tRNA, transfer RNA.

https://doi.org/10.1371/journal.pbio.3000297.s005

(TIF)

S6 Fig. Site-specific methylation shapes tRF formation to regulate protein synthesis.

(A-C) Comparison of site-specific tRNA methylation and fragmentation in NSUN2+/+ (red) and NSUN2−/− (blue) cells. (n = 4 samples per time point). Data represent median and range. Adjusted p-value: one-way ANOVA. **p < 0.005, ***p = 0.0005. (D) Log2 FC of the down-regulated tRFs when NSUN2-overexpressing cells are exposed to stress for 4 hours. tRNA lysine-derived tRFs are highlighted in red; tRNA histidine-derived tRFs are highlighted in blue. Line indicates the mean. (E) Treatment regime to measure global protein synthesis of NSUN2+/+ and −/− cells transfected with tRNA GluCTC-derived 5′ and 3′ tRFs after exposure to sodium arsenite. (F, G) Log2 FC of protein synthesis in NSUN2+/+ (F) and NSUN2−/− (G) in response to synthetic 5′ or 3′ tRFs. (n = 3 samples per time point). (H) A fluorescence siRNA was used as a control for transfection efficiency. The underlying data for this figure can be found in S11 and S12 Data and the S1 File. FC, fold-change; siRNA, small interfering RNA; tRF, tRNA-derived fragment; tRNA, transfer RNA.

https://doi.org/10.1371/journal.pbio.3000297.s006

(TIF)

S7 Fig. Mitochondrial activity is reduced and catabolic pathways enhanced in the absence of NSUN2.

(A, B) Mitochondrial activity (‘mito’) and protein synthesis (‘OP-puro’) after exposure to arsenite for the indicated time (hours) or CHX in NSUN2+/+ (A) and NSUN2−/− (B) cells. (C) GO analyses using Ribo seq data in NSUN2−/− cells rescued with NSUN2 or the enzymatic dead versions of NSUN2 C321A (left panel) and C271A (right panel). (D, E) PCA plot (D) and heatmap (E) of genes belonging to the GO: nuclear-transcribed mRNA catabolic process, nonsense-mediated decay. The underlying data for this figure can be found in S15 Data and S1 File. CHX, cycloheximide; GO, Gene Ontology; OP-puro, O-propargyl-puromycin; PCA, principle component analysis; seq, sequencing.

https://doi.org/10.1371/journal.pbio.3000297.s007

(TIF)

S1 Data. Mouse skin microarray complete dataset.

https://doi.org/10.1371/journal.pbio.3000297.s008

(XLSX)

S2 Data. Mass spectrometry raw data.

S2A: Mass spectrometry data for human dermal fibroblasts; S2B: Mass spectrometry data for mouse skin; S2C: NMR data for mouse skin. NMR, nuclear magnetic resonance.

https://doi.org/10.1371/journal.pbio.3000297.s009

(XLSX)

S3 Data. Multiple t test for metabolites shown in Fig 1I–1K.

https://doi.org/10.1371/journal.pbio.3000297.s010

(XLSX)

S4 Data. RNA sequencing data.

S4A: RNA sequencing data of NSUN2+/+, NSUN2+/−, and two lines of NSUN2−/− human dermal fibroblasts. S4B: RNA sequencing data NSUN2−/− cells reexpressing the NSUN2 protein or the empty vector as a control.

https://doi.org/10.1371/journal.pbio.3000297.s011

(XLSX)

S5 Data. RNA bisulfite sequencing of NSUN2−/− cells rescued with the wild-type NSUN2 construct or the enzymatic dead NSUN2 construct K190M, and the empty vector as a control.

https://doi.org/10.1371/journal.pbio.3000297.s012

(XLSX)

S6 Data. Small RNA sequencing of rescued NSUN2−/− cells using the wild-type or enzymatic dead version of NSUN2 or the empty vector as control.

https://doi.org/10.1371/journal.pbio.3000297.s013

(XLSX)

S7 Data. Mass spectrometry data for rescued NSUN2−/− cells using the wild-type or enzymatic dead version of NSUN2 or the empty vector as control.

https://doi.org/10.1371/journal.pbio.3000297.s014

(XLSX)

S8 Data. RNA BS sequencing data.

S8A: BS-seq tRNA methylated sites (n = 4 conversion assays) from first experimental replicate. S8B: BS-seq tRNA methylated sites (n = 5 conversion assays) from second independent experimental replicate. BS, bisulfite; BS-seq, BS sequencing; tRNA, transfer RNA.

https://doi.org/10.1371/journal.pbio.3000297.s015

(XLSX)

S9 Data. RNA BS sequencing data identifying other potentially non–tRNA targeted sites by NSUN2.

S9A: BS-seq to detect other potentially non–tRNA targeted sites by NSUN2 from 1 replicate (n = 4 conversion assays). S9B: BS-seq to detect other all sites not targeted by NSUN2 from 1 replicate (n = 4 conversion assays). BS, bisulfite; tRNA, transfer RNA.

https://doi.org/10.1371/journal.pbio.3000297.s016

(XLSX)

S10 Data. RNA BS sequencing data of unstressed and stressed NSUN2−/− cells infected with the empty vector control, the wild-type NSUN2, or the enzymatic dead construct K190M.

BS, bisulfite.

https://doi.org/10.1371/journal.pbio.3000297.s017

(XLSX)

S11 Data. Small RNA sequencing in unstressed and stressed cells.

S11A: tRNA fragments found in the small RNA-seq dataset. S11B: tRNA fragment sequencing and analysis (p < 0.01 NSUN2−/− at 2 hours of stress). S11C: tRNA fragment sequencing and analysis (samples with missing values removed). S11D: tRNA fragment sequencing and analysis (fragments smaller than 40 nucleotides). RNA-seq, RNA sequencing; tRNA, transfer RNA.

https://doi.org/10.1371/journal.pbio.3000297.s018

(XLSX)

S12 Data. tRNA fragments in NSUN2, K190M, or empty vector infected cells after 0, 2, and 4 hours of treatment with sodium arsenite.

https://doi.org/10.1371/journal.pbio.3000297.s019

(XLSX)

S13 Data. RNA sequencing data.

S13A: RNA-seq data from human dermal fibroblasts untreated (‘ctr’) or treated with sodium arsenite for 2 or 4 hours. Shown are transcriptional differences between NSUN2−/− and NSUN2 +/+ cells in the control. S13B: RNA-seq data from human dermal fibroblasts untreated (‘ctr’) or treated with sodium arsenite for 2 or 4 hours. Shown are transcriptional differences between NSUN2−/− and NSUN2 +/+ cells after 2 hours of stress. S13C: RNA-seq data from human dermal fibroblasts untreated (‘ctr’) or treated with sodium arsenite for 2 or 4 hours. Shown are transcriptional differences between NSUN2−/− and NSUN2 +/+ cells after 4 hours of stress. RNA-seq, RNA sequencing.

https://doi.org/10.1371/journal.pbio.3000297.s020

(XLSX)

S14 Data. Gene ontology analyses.

S14A: Gene enrichment (Cellular processes_Gorilla) for differentially expressed genes in NSUN2+/+ versus NSUN2−/− cells after 2 hours of exposure to sodium arsenite. S14B: Gene enrichment (Cellular processes_Gorilla) for differentially expressed genes in NSUN2+/+ versus NSUN2−/− cells after 4 hours of exposure to sodium arsenite. S14C: Gene enrichment (Cellular processes_Gorilla) for differentially expressed genes in NSUN2+/+ versus NSUN2−/− cells in untreated condition revealed no significant enrichment.

https://doi.org/10.1371/journal.pbio.3000297.s021

(XLSX)

S15 Data. Ribo seq data of NSUN2−/− cells rescued with NSUN2 (wt) or enzymatic dead versions of the NSUN2 protein (C321C, C271A) or the empty vector (‘empty’) as control.

wt, wild type.

https://doi.org/10.1371/journal.pbio.3000297.s022

(XLSX)

S16 Data. Individual mass spectrometric parameters for the nucleosides.

https://doi.org/10.1371/journal.pbio.3000297.s023

(XLSX)

S1 File. Underlying data for Figs 1C–1E, 1I–1K, 2B–2D, 2D–2N, 3C–3E, 3K–3M, 4E, 4F, 4H, 4J, 5A–5E, 5G–5P, 6D and 6G–6I and S1E, S1F, S1I–S1N, S2B, S2D, S2I–S2L, S2N, S3C, S3F–S3H, S4B, S4E, S4F, S4H, S4I, S5A, S5D, S5E, S6A–S6D, S6F–S6H, S7A–S7C and S7E Figs.

https://doi.org/10.1371/journal.pbio.3000297.s024

(XLSX)

Acknowledgments

We gratefully acknowledge the support of the Wellcome-MRC Cambridge Stem Cell Institute facilities and, in particular for this work, Maike Paramor, Irina Pshenichnaya, Gemini Chu, Peter Humphreys, and Andy Riddell. We also thank the Cancer Research UK Cambridge Institute Genomics Core facility for their support. S. K. thanks Thomas Carell and his group for instrument time (LC-MS/MS).

Citation: Gkatza NA, Castro C, Harvey RF, Heiß M, Popis MC, Blanco S, et al. (2019) Cytosine-5 RNA methylation links protein synthesis to cell metabolism. PLoS Biol 17(6): e3000297. https://doi.org/10.1371/journal.pbio.3000297

References

1. Van Heyningen V, Yeyati PL. Mechanisms of non-Mendelian inheritance in genetic disease. Hum Mol Genet. 2004;13 Spec No 2:R225–33. pmid:15358729.

2. Proudfoot NJ, Furger A, Dye MJ. Integrating mRNA processing with transcription. Cell. 2002;108(4):501–12. pmid:11909521.

3. Frye M, Harada BT, Behm M, He C. RNA modifications modulate gene expression during development. Science. 2018;361(6409):1346–9. Epub 2018/09/29. pmid:30262497.

4. Boccaletto P, Machnicka MA, Purta E, Piatkowski P, Baginski B, Wirecki TK, et al. MODOMICS: a database of RNA modification pathways. 2017 update. Nucleic Acids Res. 2018;46(D1):D303–D7. pmid:29106616; PubMed Central PMCID: PMC5753262.

5. Frye M, Blanco S. Post-transcriptional modifications in development and stem cells. Development. 2016;143(21):3871–81. pmid:27803056.

6. Khan MA, Rafiq MA, Noor A, Hussain S, Flores JV, Rupp V, et al. Mutation in NSUN2, which Encodes an RNA Methyltransferase, Causes Autosomal-Recessive Intellectual Disability. Am J Hum Genet. 2012;90(5):856–63. Epub 2012/05/01. S0002-9297(12)00203-0 [pii] pmid:22541562.

7. Martinez FJ, Lee JH, Lee JE, Blanco S, Nickerson E, Gabriel S, et al. Whole exome sequencing identifies a splicing mutation in NSUN2 as a cause of a Dubowitz-like syndrome. J Med Genet. 2012;49(6):380–5. pmid:22577224.

8. Van Haute L, Dietmann S, Kremer L, Hussain S, Pearce SF, Powell CA, et al. Deficient methylation and formylation of mt-tRNA(Met) wobble cytosine in a patient carrying mutations in NSUN3. Nat Commun. 2016;7:12039. pmid:27356879.

9. Nakano S, Suzuki T, Kawarada L, Iwata H, Asano K, Suzuki T. NSUN3 methylase initiates 5-formylcytidine biogenesis in human mitochondrial tRNA(Met). Nature chemical biology. 2016;12(7):546–51. pmid:27214402.

10. Schaefer M, Pollex T, Hanna K, Tuorto F, Meusburger M, Helm M, et al. RNA methylation by Dnmt2 protects transfer RNAs against stress-induced cleavage. Genes Dev. 2010;24(15):1590–5. Epub 2010/08/04. 24/15/1590 [pii] pmid:20679393.

11. Blanco S, Dietmann S, Flores JV, Hussain S, Kutter C, Humphreys P, et al. Aberrant methylation of tRNAs links cellular stress to neuro-developmental disorders. The EMBO journal. 2014;33(18):2020–39. pmid:25063673; PubMed Central PMCID: PMC4195770.

12. Blanco S, Bandiera R, Popis M, Hussain S, Lombard P, Aleksic J, et al. Stem cell function and stress response are controlled by protein synthesis. Nature. 2016;534(7607):335–40. pmid:27306184.

13. Ivanov P, Emara MM, Villen J, Gygi SP, Anderson P. Angiogenin-induced tRNA fragments inhibit translation initiation. Mol Cell. 2011;43(4):613–23. Epub 2011/08/23. S1097-2765(11)00524-7 [pii] pmid:21855800; PubMed Central PMCID: PMC3160621.

14. Flores JV, Cordero-Espinoza L, Oeztuerk-Winder F, Andersson-Rolf A, Selmi T, Blanco S, et al. Cytosine-5 RNA Methylation Regulates Neural Stem Cell Differentiation and Motility. Stem Cell Reports. 2017;8(1):112–24. pmid:28041877; PubMed Central PMCID: PMC5233436.

15. Llorens-Bobadilla E, Zhao S, Baser A, Saiz-Castro G, Zwadlo K, Martin-Villalba A. Single-Cell Transcriptomics Reveals a Population of Dormant Neural Stem Cells that Become Activated upon Brain Injury. Cell Stem Cell. 2015;17(3):329–40. pmid:26235341.

16. Signer RA, Magee JA, Salic A, Morrison SJ. Haematopoietic stem cells require a highly regulated protein synthesis rate. Nature. 2014;509(7498):49–54. pmid:24670665; PubMed Central PMCID: PMC4015626.

17. Baser A, Skabkin M, Kleber S, Dang Y, Gulculer Balta GS, Kalamakis G, et al. Onset of differentiation is post-transcriptionally controlled in adult neural stem cells. Nature. 2019;566(7742):100–4. Epub 2019/02/01. pmid:30700908.

18. Blanco S, Kurowski A, Nichols J, Watt FM, Benitah SA, Frye M. The RNA-methyltransferase Misu (NSun2) poises epidermal stem cells to differentiate. PLoS Genet. 2011;7(12):e1002403. Epub 2011/12/07. PGENETICS-D-11-00619 [pii]. pmid:22144916; PubMed Central PMCID: PMC3228827.

19. Trempus CS, Morris RJ, Bortner CD, Cotsarelis G, Faircloth RS, Reece JM, et al. Enrichment for living murine keratinocytes from the hair follicle bulge with the cell surface marker CD34. J Invest Dermatol. 2003;120(4):501–11. Epub 2003/03/22. 12088 [pii] pmid:12648211.

20. Blanpain C, Lowry WE, Geoghegan A, Polak L, Fuchs E. Self-renewal, multipotency, and the existence of two cell populations within an epithelial stem cell niche. Cell. 2004;118(5):635–48. pmid:15339667.

21. Greco V, Chen T, Rendl M, Schober M, Pasolli HA, Stokes N, et al. A two-step mechanism for stem cell activation during hair regeneration. Cell Stem Cell. 2009;4(2):155–69. pmid:19200804.

22. Bujnicki JM, Feder M, Ayres CL, Redman KL. Sequence-structure-function studies of tRNA:m5C methyltransferase Trm4p and its relationship to DNA:m5C and RNA:m5U methyltransferases. Nucleic Acids Res. 2004;32(8):2453–63. pmid:15121902.

23. Hussain S, Sajini AA, Blanco S, Dietmann S, Lombard P, Sugimoto Y, et al. NSun2-mediated cytosine-5 methylation of vault noncoding RNA determines its processing into regulatory small RNAs. Cell reports. 2013;4(2):255–61. pmid:23871666; PubMed Central PMCID: PMC3730056.

24. Yang X, Yang Y, Sun BF, Chen YS, Xu JW, Lai WY, et al. 5-methylcytosine promotes mRNA export—NSUN2 as the methyltransferase and ALYREF as an m(5)C reader. Cell Res. 2017;27(5):606–25. pmid:28418038; PubMed Central PMCID: PMC5594206.

25. Ducker GS, Rabinowitz JD. One-Carbon Metabolism in Health and Disease. Cell metabolism. 2017;25(1):27–42. pmid:27641100; PubMed Central PMCID: PMC5353360.

26. Kalhan SC. One carbon metabolism in pregnancy: Impact on maternal, fetal and neonatal health. Molecular and cellular endocrinology. 2016;435:48–60. pmid:27267668; PubMed Central PMCID: PMC5014566.

27. Yang M, Vousden KH. Serine and one-carbon metabolism in cancer. Nat Rev Cancer. 2016;16(10):650–62. pmid:27634448.

28. Mentch SJ, Mehrmohamadi M, Huang L, Liu X, Gupta D, Mattocks D, et al. Histone Methylation Dynamics and Gene Regulation Occur through the Sensing of One-Carbon Metabolism. Cell metabolism. 2015;22(5):861–73. pmid:26411344; PubMed Central PMCID: PMC4635069.

29. Teerlink T, Luo Z, Palm F, Wilcox CS. Cellular ADMA: regulation and action. Pharmacol Res. 2009;60(6):448–60. pmid:19682580; PubMed Central PMCID: PMC2767414.

30. Hussain S, Benavente SB, Nascimento E, Dragoni I, Kurowski A, Gillich A, et al. The nucleolar RNA methyltransferase Misu (NSun2) is required for mitotic spindle stability. J Cell Biol. 2009;186(1):27–40. pmid:19596847.

31. Lyons SM, Fay MM, Akiyama Y, Anderson PJ, Ivanov P. RNA biology of angiogenin: Current state and perspectives. RNA Biol. 2017;14(2):171–8. Epub 2016/12/25. pmid:28010172; PubMed Central PMCID: PMC5324756.

32. Liu J, Xu Y, Stoleru D, Salic A. Imaging protein synthesis in cells and tissues with an alkyne analog of puromycin. Proc Natl Acad Sci U S A. 2012;109(2):413–8. pmid:22160674; PubMed Central PMCID: PMC3258597.

33. Schneider-Poetsch T, Ju J, Eyler DE, Dang Y, Bhat S, Merrick WC, et al. Inhibition of eukaryotic translation elongation by cycloheximide and lactimidomycin. Nature chemical biology. 2010;6(3):209–17. pmid:20118940; PubMed Central PMCID: PMC2831214.

34. Orioli A. tRNA biology in the omics era: Stress signalling dynamics and cancer progression. Bioessays. 2017;39(3). pmid:28026052.

35. Kong Q, Lin CL. Oxidative damage to RNA: mechanisms, consequences, and diseases. Cell Mol Life Sci. 2010;67(11):1817–29. pmid:20148281; PubMed Central PMCID: PMC3010397.

36. Fulda S, Gorman AM, Hori O, Samali A. Cellular stress responses: cell survival and cell death. Int J Cell Biol. 2010;2010:214074. pmid:20182529; PubMed Central PMCID: PMC2825543.

37. McCollum G, Keng PC, States JC, McCabe MJ Jr., Arsenite delays progression through each cell cycle phase and induces apoptosis following G2/M arrest in U937 myeloid leukemia cells. J Pharmacol Exp Ther. 2005;313(2):877–87. pmid:15722406.

38. Yang K, Wang M, Zhao Y, Sun X, Yang Y, Li X, et al. A redox mechanism underlying nucleolar stress sensing by nucleophosmin. Nat Commun. 2016;7:13599. pmid:27886181; PubMed Central PMCID: PMC5133708.

39. Frye M, Watt FM. The RNA methyltransferase Misu (NSun2) mediates Myc-induced proliferation and is upregulated in tumors. Curr Biol. 2006;16(10):971–81. pmid:16713953.

40. Pollex T, Hanna K, Schaefer M. Detection of cytosine methylation in RNA using bisulfite sequencing. Cold Spring Harbor Protocols. 2010;5(10). pmid:20889702

41. Kao RY, Jenkins JL, Olson KA, Key ME, Fett JW, Shapiro R. A small-molecule inhibitor of the ribonucleolytic activity of human angiogenin that possesses antitumor activity. Proc Natl Acad Sci U S A. 2002;99(15):10066–71. Epub 2002/07/16. [pii]. pmid:12118120; PubMed Central PMCID: PMC126625.

42. Cole C, Sobala A, Lu C, Thatcher SR, Bowman A, Brown JW, et al. Filtering of deep sequencing data reveals the existence of abundant Dicer-dependent small RNAs derived from tRNAs. RNA. 2009;15(12):2147–60. pmid:19850906; PubMed Central PMCID: PMC2779667.

43. Fu H, Feng J, Liu Q, Sun F, Tie Y, Zhu J, et al. Stress induces tRNA cleavage by angiogenin in mammalian cells. FEBS Lett. 2009;583(2):437–42. Epub 2008/12/31. S0014-5793(08)01035-1 [pii] pmid:19114040.

44. Goodarzi H, Liu X, Nguyen HC, Zhang S, Fish L, Tavazoie SF. Endogenous tRNA-Derived Fragments Suppress Breast Cancer Progression via YBX1 Displacement. Cell. 2015;161(4):790–802. pmid:25957686; PubMed Central PMCID: PMC4457382.

45. Gebetsberger J, Wyss L, Mleczko AM, Reuther J, Polacek N. A tRNA-derived fragment competes with mRNA for ribosome binding and regulates translation during stress. RNA Biol. 2017;14(10):1364–73. pmid:27892771; PubMed Central PMCID: PMC5711459.

46. Sharma LK, Lu J, Bai Y. Mitochondrial respiratory complex I: structure, function and implication in human diseases. Curr Med Chem. 2009;16(10):1266–77. pmid:19355884; PubMed Central PMCID: PMC4706149.

47. Lindqvist LM, Tandoc K, Topisirovic I, Furic L. Cross-talk between protein synthesis, energy metabolism and autophagy in cancer. Curr Opin Genet Dev. 2018;48:104–11. pmid:29179096; PubMed Central PMCID: PMC5869074.

48. Buszczak M, Signer RA, Morrison SJ. Cellular differences in protein synthesis regulate tissue homeostasis. Cell. 2014;159(2):242–51. pmid:25303523; PubMed Central PMCID: PMC4222182.

49. Pan T. Modifications and functional genomics of human transfer RNA. Cell Res. 2018;28(4):395–404. pmid:29463900.

50. Chan CT, Pang YL, Deng W, Babu IR, Dyavaiah M, Begley TJ, et al. Reprogramming of tRNA modifications controls the oxidative stress response by codon-biased translation of proteins. Nat Commun. 2012;3:937. Epub 2012/07/05. ncomms1938 [pii] pmid:22760636; PubMed Central PMCID: PMC3535174.

51. Endres L, Dedon PC, Begley TJ. Codon-biased translation can be regulated by wobble-base tRNA modification systems during cellular stress responses. RNA Biol. 2015;12(6):603–14. pmid:25892531.

52. Ranjan N, Rodnina MV. tRNA wobble modifications and protein homeostasis. Translation (Austin). 2016;4(1):e1143076. pmid:27335723; PubMed Central PMCID: PMC4909423.

53. Phizicky EM, Hopper AK. tRNA biology charges to the front. Genes Dev. 2010;24(17):1832–60. Epub 2010/09/03. 24/17/1832 [pii] pmid:20810645.

54. Guzzi N, Ciesla M, Ngoc PCT, Lang S, Arora S, Dimitriou M, et al. Pseudouridylation of tRNA-Derived Fragments Steers Translational Control in Stem Cells. Cell. 2018. pmid:29628141.

55. Tuorto F, Liebers R, Musch T, Schaefer M, Hofmann S, Kellner S, et al. RNA cytosine methylation by Dnmt2 and NSun2 promotes tRNA stability and protein synthesis. Nat Struct Mol Biol. 2012;19(9):900–5. Epub 2012/08/14. nsmb.2357 [pii] pmid:22885326.

56. Kurimoto K, Yabuta Y, Ohinata Y, Saitou M. Global single-cell cDNA amplification to provide a template for representative high-density oligonucleotide microarray analysis. Nat Protoc. 2007;2(3):739–52. pmid:17406636

57. Wang F, Flanagan J, Su N, Wang LC, Bui S, Nielson A, et al. RNAscope: A novel in situ RNA analysis platform for formalin-fixed, paraffin-embedded tissues. Journal of Molecular Diagnostics. 2012;14(1):22–9. pmid:22166544

58. Ament Z, West JA, Stanley E, Ashmore T, Roberts LD, Wright J, et al. PPAR-pan activation induces hepatic oxidative stress and lipidomic remodelling. Free Radical Biology and Medicine. 2016;95:357–68. pmid:26654758

59. Le Belle JE, Harris NG, Williams SR, Bhakoo KK. A comparison of cell and tissue extraction techniques using high-resolution 1H-NMR spectroscopy. NMR in Biomedicine. 2002;15(1):37–44. pmid:11840551

60. Mazzon M, Castro C, Thaa B, Liu L, Mutso M, Liu X, et al. Alphavirus-induced hyperactivation of PI3K/AKT directs pro-viral metabolic changes. PLoS Pathog. 2018;14(1):e1006835. Epub 2018/01/30. pmid:29377936; PubMed Central PMCID: PMC5805360.

61. Heiss M, Reichle VF, Kellner S. Observing the fate of tRNA and its modifications by nucleic acid isotope labeling mass spectrometry: NAIL-MS. RNA Biol. 2017;14(9):1260–8. Epub 2017/05/11. pmid:28488916; PubMed Central PMCID: PMC5699550.

62. Cai WM, Chionh YH, Hia F, Gu C, Kellner S, McBee ME, et al. A Platform for Discovery and Quantification of Modified Ribonucleosides in RNA: Application to Stress-Induced Reprogramming of tRNA Modifications. Methods Enzymol. 2015;560:29–71. Epub 2015/08/09. pmid:26253965; PubMed Central PMCID: PMC4774897.

63. Borland K, Diesend J, Ito-Kureha T, Heissmeyer V, Hammann C, Buck AH, et al. Production and Application of Stable Isotope-Labeled Internal Standards for RNA Modification Analysis. Genes (Basel). 2019;10(1). Epub 2019/01/10. pmid:30621251; PubMed Central PMCID: PMC6356711.

64. Hauenschild R, Tserovski L, Schmid K, Thuring K, Winz ML, Sharma S, et al. The reverse transcription signature of N-1-methyladenosine in RNA-Seq is sequence dependent. Nucleic Acids Res. 2015;43(20):9950–64. pmid:26365242.

65. Kuleshov MV, Jones MR, Rouillard AD, Fernandez NF, Duan Q, Wang Z, et al. Enrichr: a comprehensive gene set enrichment analysis web server 2016 update. Nucleic Acids Res. 2016;44(W1):W90–7. pmid:27141961; PubMed Central PMCID: PMC4987924.

66. Chen EY, Tan CM, Kou Y, Duan Q, Wang Z, Meirelles GV, et al. Enrichr: interactive and collaborative HTML5 gene list enrichment analysis tool. BMC Bioinformatics. 2013;14:128. pmid:23586463; PubMed Central PMCID: PMC3637064.

67. Eden E, Navon R, Steinfeld I, Lipson D, Yakhini Z. GOrilla: a tool for discovery and visualization of enriched GO terms in ranked gene lists. BMC Bioinformatics. 2009;10:48. pmid:19192299; PubMed Central PMCID: PMC2644678.