#PAGE_PARAMS# #ADS_HEAD_SCRIPTS# #MICRODATA#

Genome-wide identification of short 2′,3′-cyclic phosphate-containing RNAs and their regulation in aging


Authors: Megumi Shigematsu aff001;  Keisuke Morichika aff001;  Takuya Kawamura aff001;  Shozo Honda aff001;  Yohei Kirino aff001
Authors place of work: Computational Medicine Center, Sidney Kimmel Medical College, Thomas Jefferson University, Philadelphia, Pennsylvania, United States of America aff001
Published in the journal: Genome-wide identification of short 2′,3′-cyclic phosphate-containing RNAs and their regulation in aging. PLoS Genet 15(11): e32767. doi:10.1371/journal.pgen.1008469
Category: Research Article
doi: https://doi.org/10.1371/journal.pgen.1008469

Summary

RNA molecules generated by ribonuclease cleavage sometimes harbor a 2′,3′-cyclic phosphate (cP) at their 3′-ends. Those cP-containing RNAs (cP-RNAs) form a hidden layer of transcriptome because standard RNA-seq cannot capture them as a result of cP’s prevention of an adapter ligation reaction. Here we provide genome-wide analyses of short cP-RNA transcriptome across multiple mouse tissues. Using cP-RNA-seq that can exclusively sequence cP-RNAs, we identified numerous novel cP-RNA species which are mainly derived from cytoplasmic tRNAs, mRNAs, and rRNAs. Determination of the processing sites of substrate RNAs for cP-RNA generation revealed highly-specific RNA cleavage events between cytidine and adenosine in cP-RNA biogenesis. cP-RNAs were not evenly derived from the overall region of substrate RNAs but rather from specific sites, implying that cP-RNAs are not from random degradation but are produced through a regulated biogenesis pathway. The identified cP-RNAs were abundantly accumulated in mouse tissues, and the expression levels of cP-RNAs showed age-dependent reduction. These analyses of cP-RNA transcriptome unravel a novel, abundant class of non-coding RNAs whose expression could have physiological roles.

Keywords:

Sequence assembly tools – Non-coding RNA – Mammalian genomics – RNA sequencing – Ribosomal RNA – Non-coding RNA sequences – Transfer RNA – Small nucleolar RNA

Introduction

Cells express immensely diverse species of RNA molecules that play essential roles in numerous biological processes. With the advent of next-generation sequencing (NGS) technology, efforts to identify the expressed RNA molecules have greatly advanced our understanding of RNA biology [1]. Although NGS-based RNA-sequencing (RNA-seq) has become a ubiquitous tool in biological and medical research [2,3], the current standard RNA-seq methods, particularly those targeting short non-coding RNAs (ncRNAs), do not fully capture all of the RNAs expressed but allow for some “escapers” to slip by undetected. The RNAs containing a 2′,3′-cyclic phosphate (cP) at their 3′-end are one such escaper that are not ligated to a 3′-adapter during cDNA amplification procedure and thus uncaptured by standard RNA-seq. A cP end of RNA molecules, consisting of the phosphate linkage between the 2′- and 3′-positions of ribose, is mainly produced from RNA cleavage catalyzed by many endoribonucleases [4]. cP-containing RNAs (cP-RNAs) are not just accumulated as non-functional degradation products but have physiological roles in various biological processes. For example, Angiogenin (ANG), an RNase A-superfamily endoribonuclease, cleaves the anticodon-loop of tRNAs upon stress stimuli to produce functional tRNA halves known as tRNA-derived stress-induced RNAs (tiRNAs) [5]. The 5′-tiRNAs, corresponding to the 5′-tRNA halves, contain a cP at their 3′-end and promote the formation of stress granules, regulate translation, and trigger the cellular stress responses and apoptosis in neurodevelopmental disorders [6,7,8,9]. In breast and prostate cancer cells, sex hormone signaling pathways promote ANG-catalyzed tRNA cleavages, generating a distinct class of tRNA halves, the sex hormone-dependent tRNA derived RNAs (SHOT-RNAs) [10]. cP-containing 5′-SHOT-RNAs (5′-tRNA halves) have a functional significance in cell proliferation [10,11].

Given that cP-RNAs are expressed as functional molecules, capturing the entire repertoire of cP-RNAs would broaden the catalog of functional ncRNAs and reveal significant biological events that have been eluding standard RNA-seq. We recently developed “cP-RNA-seq,” which is able to specifically sequence cP-RNAs [10,12], and applied the method to a short RNA fraction of BT-474 breast cancer cells to identify 5′-SHOT-RNAs [10]. Because the previous work narrowly focused on 5′-SHOT-RNAs in only one cell line, cP-RNAs still largely remained a hidden component in the transcriptome. In this report, we expanded populations of RNA subjects and identified a comprehensive repertoire of short cP-RNAs in four different mouse tissues. A genome-wide in-depth view of cP-RNA transcriptome revealed that not only tRNAs but also other RNA species are widely, but specifically, cleaved between cytidine and adenosine for constitutive expression of diverse species of cP-RNAs. We further discovered age-dependent changes in the expression levels of cP-RNAs, implying the involvement of cP-RNA production in the aging process.

Results and discussion

cP-containing tRNA halves are abundantly and constitutively expressed in mouse tissues

Because ANG-generated 5′-tRNA halves contain a cP, exploration of the accumulation of cP-RNAs in mouse tissues began with analyses of tRNA halves. Northern blots for cytoplasmic (cyto) tRNALysCUU and tRNAAspGUC, both of which abundantly produce tRNA halves in human breast cancer cells [10], revealed constitutive expression of both the 5′- and 3′-halves in various mouse tissues (Figs 1A and S1A). Although the levels of mature tRNAs were comparable across all the examined tissues, tRNA halves showed differential expression patterns among the tissues, implying a regulated biogenesis of tRNA halves. The abundant accumulation of 5′-tRNA halves found in the spleen and less abundant accumulations in other tissues, such as the brain, heart, liver, and testis, are consistent with previous reports [13,14]. To confirm the presence of a cP in the 5′-tRNA halves, total RNA from mouse lungs was treated with T4 polynucleotide kinase (T4 PNK; removes 3′-P/cP), calf intestinal phosphatase (CIP; removes 5′-P/3′-P, but not a cP), or acid followed by CIP treatment (removes 5′-P/3′-P/cP), and then the mobility of the bands of the 5′-tRNA half was examined by northern blots, as previously described [10,12]. As shown in Fig 1B, the band of the 5′-tRNA half shifted up by CIP treatment, suggesting the presence of a P in the 5′-tRNA half. T4 PNK treatment also shifted the band up, and further up-shift was observed by the acid plus CIP treatments, indicating the additional presence of a cP in the 5′-tRNA half. These results validated that, as observed for the 5′-tRNA halves in human cell lines [10], the 5′-tRNA halves in mouse tissues contain a P at their 5′-end and a cP at their 3′-end.

Sequencing of cP-RNAs in mouse tissues.
Fig. 1. Sequencing of cP-RNAs in mouse tissues.
(A) Total RNAs extracted from mouse tissues were subjected to northern blots for the 5′-halves of cyto tRNALysCUU and tRNAAspGUC. (B) Terminal structures of the 5′-tRNA half were analyzed enzymatically. Total RNA from the mouse lung was treated with CIP, T4 PNK, or acid followed by CIP treatment (HCl + CIP). NT designates the non-treated sample used as a negative control. The treated total RNA was subjected to northern blots targeting the 5′-tRNAAspGUC half and microRNA-16 (miR-16). miR-16 was investigated as a control RNA containing 5′-P and 3′-OH ends. (C) Gel-purified 20–45-nt RNAs were subjected to cP-RNA-seq, which amplified 140–160-bp cDNA products (5′-adapter, 55 bp; 3′-adapter, 63 bp; and thereby estimated inserted sequences, 22–42 bp). The cDNAs in the region designated with a red line were purified and subjected to Illumina sequencing. (D) Proportion of cP-RNAs annotated to the indicated RNAs.

Selective amplification and sequencing of short cP-RNAs expressed in mouse tissues

To capture an entire picture of short cP-RNAs, we selected four tissues (thymus, lung, spleen, and skeletal muscle) that exhibited an abundant expression of 5′-tRNA halves (Fig 1A). The four tissues commonly showed accumulations of short RNAs that were expected to include 5′-tRNA halves (S1B Fig), whereas we did not observe such short RNA accumulation in HEK293T and HeLa cells (S1B Fig), implying that the constitutive accumulation of short RNAs might be tissue specific. We chose to focus on the 20 to 45-nucleotide (nt) RNAs which were gel-purified and subjected to cP-RNA-seq [10,12]. The method successfully amplified ~140 to 160-bp cDNA products (inserted sequences without adapters were estimated to be 22–42 bp) whose dependency on T4 PNK treatment indicates that, as expected, the amplified bands were derived from cP-RNAs (Fig 1C). Illumina sequencing of the amplified bands from the respective four tissues yielded approximately 38–47 million raw reads, of which >90%–97% were extracted as the cP-RNA reads with a length of 20–45 nt (S1 Table). Read lengths of the extracted cP-RNAs similarly showed broad distribution patterns in different tissues (S2 Fig). Mapping the cP-RNA reads to the mouse ncRNAs and genome unexpectedly revealed that not only tRNAs but also the other RNA species serve as abundant substrates for cP-RNA production (Fig 1D). Although tRNA-derived reads were the most abundant cP-RNA species and occupied a substantial proportion of cP-RNAs, ~55%–80% reads are mapped to other regions of the genome, of which messenger RNAs (mRNAs) and ribosomal RNAs (rRNAs) were particularly rich sources of cP-RNAs (Fig 1D and S1 Table).

tRNA-derived cP-RNAs are generated from focal cleavage of anticodon-loop and 3′-terminal CCA sequence of tRNAs

tRNA-derived ncRNAs can be classified into 5′-tRNA half, 3′-tRNA half, 5′-tRF, 3′-tRF, and internal-tRF (i-tRF) [15,16,17,18] (S3A Fig). tRNA-derived cP-RNAs were mostly 30–40 nt long (S2 Fig) and the majority were 5′-halves (Fig 2A) as expected. Considering that the mouse genome encodes 55 cyto tRNA isoacceptors with different anticodon sequences [19], the identified 5′-halves were derived from a rather focused subset of tRNAs, such as cyto tRNAValCAC, tRNAGlyGCC, tRNALysCUU, and tRNAGluCUC, which are in aggregates the sources of 75%–84% of the identified 5′-halves (Fig 2B). The 3′-terminal nucleotides of the 5′-halves had a strongly biased composition in being pyrimidines (C > U), which was a stable property in all four examined tissues (Fig 2C). The anticodon-loop nucleotides that are one position downstream of the 3′-end showed biases of A and U, allowing us to predict that the anticodon-loop cleavage mainly occurs between C and A/U to produce cP-containing 5′-halves (Figs 2C and 2D and S3B).

Analyses of tRNA-derived cP-RNAs.
Fig. 2. Analyses of tRNA-derived cP-RNAs.
(A) Proportion of tRNA-derived cP-RNAs classified into the indicated subgroups of tRNA-derived ncRNAs shown in S3A Fig. (B) Proportion of the 5′-half-reads derived from respective cyto tRNA species. (C) Compositions of the 3′-terminal nucleotide of the 5′-halves and its 3′-adjacent nucleotide (potential 5′-terminal nucleotide of the 5′-cleavage product) in the anticodon-loop of tRNAs. (D) Cleavage sites in the anticodon-loops in the thymus were predicted based on the 3′-terminal positions of the 5′-halves (blue arrow heads) and the 5′-terminal positions of the 3′-halves (red arrow heads). The results for the other three tissues are shown in S3B Fig. (E) Proportion of the 3′-halves classified into the indicated subgroups shown in S3A Fig. The numerical data is shown in S2 Table.

Among all species of tRNAs, tRNAHisGUG is unique in that it contains an additional nucleotide at nucleotide position (np; according to the nucleotide numbering system of tRNAs [20]) –1 of its 5′-end. Our recent TaqMan RT-qPCR analyses revealed that, in BT-474 cells, the majority (~60%) of cyto tRNAHisGUG contains G–1, but a significant proportion contains U–1 or lacks the –1 nucleotide, resulting in G1 as a 5′-terminal nucleotide [21]. As shown in S3C Fig, a majority of the identified 5′-halves of cyto tRNAHisGUG contained G–1, but the 5′-halves containing U–1 or lacking –1 nucleotide (G1) were also substantially detected, suggesting that the presence of U–1 as well as G–1 is a conserved feature in mouse cyto tRNAHisGUG, which is retained in the identified 5′-halves.

It was unexpected that a substantial proportion (9%–31%) of tRNA-derived cP-RNAs were occupied by 3′-halves. All three tRFs (5′-, 3′-, and i-tRFs) were identified as minor cP-containing species compared to the tRNA halves. When subclassified into the four subgroups (3′-CCA-/3′-CC-/3′-C-/3′-Di-halves; S3A Fig), nearly all (89%–98%) of the identified 3′-halves were 3′-CC-halves (Fig 2E and S2 Table), suggesting that the 3′-terminal cP formation resulted from the cleavage between C75 and A76 of tRNAs. The in vitro ANG-catalyzed cleavage between C75 and A76 was previously reported [22], and our results might suggest the occurrence of cleavage in vivo. The 3′-CC- and 3′-C-halves were identified as cP-RNAs also in human brain [23], suggesting the consistency between mouse and human tissues. The anticodon-loop cleavage sites predicted by the 3′-terminal position of the 5′-halves and the 5′-terminal position of the 3′-halves were sometimes consistent (Figs 2D and S3B). Therefore, specific tRNAs are likely to be selected for cP-RNA production and at least some of their 5′- and 3′-half pairs could be concurrently produced by endonucleolytic cleavage at the anticodon-loops.

mRNA-derived cP-RNAs are produced from a specific processing between adenine and cytosine

Although not anticipated, not only tRNAs but also other RNA species were identified as rich sources of cP-RNAs. mRNA-derived cP-RNAs were enriched in cP-RNA libraries from all four examined tissues, accounting for 15%–20% of all cP-RNA species (Fig 1D). The lengths of the mRNA-derived cP-RNAs were broadly distributed between 20–45 nt (S2 Fig). Both coding and untranslated regions of mRNAs became substrates for cP-RNAs (S4A Fig). Astonishing nucleotide composition biases were observed at both the 5′- and 3′-ends of mRNA-derived cP-RNAs in all four tissues (Figs 3A and S4B). The 5′-ends were mostly A (78%–89%), whereas the 3′-ends were mostly pyrimidines, and C (70%–81%) was especially enriched compared to U (17%–29%). None of the biases were observed in the nucleotides from the second position to the position preceding the 3′-ends. The nucleotides that are one position upstream of the 5′-ends and one position downstream of the 3′-ends showed clear A and C biases (Fig 3B), strongly suggesting that the mRNA-derived cP-RNAs are produced by an endoribonuclease that specifically recognizes and cleaves mRNAs between C and A residues.

Nucleotide composition analyses of mRNA- and rRNA-derived cP-RNAs.
Fig. 3. Nucleotide composition analyses of mRNA- and rRNA-derived cP-RNAs.
(A) Nucleotide compositions of the first 30 nucleotides from the 5′-end (upper) or 3′-end (lower) of the mRNA-derived cP-RNAs in the thymus and lung. The results for the other two tissues are shown in S4B Fig. (B) Nucleotide compositions around the 5′- and 3′-ends of mRNA-derived cP-RNAs. A dashed line separates upstream (-) and downstream (+) positions for the 5′-and 3′-ends, representing the cleavage site that generates mRNA-derived cP-RNAs (the regions outside of cP-RNA-generating regions are colored in grey). (C) Nucleotide compositions around the 5′- and 3′-ends of rRNA-derived cP-RNAs.

Major mRNA species producing cP-RNAs are not very consistent among different tissues (S3 Table), and the cP-RNAs were derived from many different regions throughout the mRNAs (S4C Fig). In the spleen, many abundant cP-RNAs were derived from hemoglobin mRNAs, which might be correlated to the role of the spleen in red blood cell turnover that requires the degradation of hemoglobin [24]. It should be noted that some mRNAs, such as Hist1h2br and B2m in the thymus (S4C Fig), showed very focal production of specific cP-RNAs, suggesting that the cP-RNA production and accumulation are not caused by random mRNA degradation. The most prominent example is Ncbp3 mRNA which commonly produces the most abundant mRNA-derived cP-RNA in all four tissues (S3 Table). Although there are 406 “CA” sequences (potential cleavage sites) throughout Ncbp3 mRNAs, almost all cP-RNAs are produced from only one site (position 5780–5809; the generated cP-RNA is termed cPR-Ncbp3) in the 3′-UTR (Fig 4A), where single-stranded loop regions are expected to be formed (Fig 4B). The common and specific production of cP-RNAs might suggest that these RNAs are produced as functional molecules.

TaqMan RT-qPCR quantification of cP-RNAs.
Fig. 4. TaqMan RT-qPCR quantification of cP-RNAs.
(A) The alignment patterns of cP-RNAs for Ncbp3 mRNA and 28S rRNA. The positions of representative cP-RNAs, cPR-Ncbp3 and cPR-28S, are indicated. (B) The regions from which 5′-tRNAGlyGCC half (5′-GlyGCC), cPR-Ncbp3, and cPR-28S were derived are shown in red in the secondary structure of respective substrate RNAs. Secondary structure of Ncbp3 mRNA was predicted by ViennaRNA Package 2.0 [34]. (C) The total RNA from 24-week old mouse skeletal muscle, treated with CIP, wild-type (WT) T4 PNK, or mutant (Mut) T4 PNK lacking 3′-dephosphorylation activity, was subjected to TaqMan RT-qPCR. NT designates a non-treated sample used as a negative control. The amounts from WT T4 PNK-treated RNA were set as 1, and relative amounts are indicated. Averages of three experiments with SD values are shown. (D) The expression of cP-RNAs in the lung and skeletal muscle of 24-week old mice were quantified using TaqMan RT-qPCR. The cP-RNA amounts were estimated based on the standard curves shown in S9 Fig. Averages of three independent experiments with SD values are shown.

cP-RNAs produced from other RNA species are also produced from a specific processing between adenine and cytosine

rRNA-derived cP-RNAs, whose read lengths were broadly distributed between 20–45 nt (S2 Fig), were enriched at similar levels (15%–25% of all cP-RNA species; Fig 1D) with mRNA-derived cP-RNAs. All four rRNA species (28S, 18S, 5.8S, and 5S) produce cP-RNAs (S5A Fig), and the most abundant cP-RNA species were derived from 28S rRNA (S4 Table). Nucleotide biases observed in mRNA-derived cP-RNAs were retained mostly in rRNA-derived cP-RNAs; the 3′-end and a nucleotide upstream of the 5′-end were strongly biased to C and U (C > U) (Fig 3C). The 5′-end and a nucleotide downstream of the 3′-end were not only biased to A, but also to G, and sometimes to U. Alignment visualization showed production of specific cP-RNA species from rRNAs (Figs 4A and S5B), suggesting their controlled biogenesis mechanisms. In 28S rRNA, the positions 186–217 was the most abundant cP-RNA-producing region (Fig 4A), and the representative cP-RNA from the region, termed cPR-28S, is produced from the cleavages between C185 and G186, and C217 and A218 in domain 1 (Fig 4B). Cleavage positions reside in the loop regions of 28S rRNA, possibly because the regions are highly accessible to endoribonucleases.

cP-RNAs from small nuclear RNAs (snRNAs), small nucleolar RNAs (snoRNAs), microRNAs (miRNAs), and long ncRNAs (lncRNAs) showed terminal nucleotide bias to A and C for their 5′- and 3′-end, respectively, in all four tissues (S6 and S7 Figs). Moreover, not only the cP-RNAs derived from the nuclear genome but also the cP-RNAs derived from the mitochondrial genome showed similar nucleotide biases (S8 Fig). Although those cP-RNAs were not abundant, the consistent nucleotide biases suggest that cP-RNAs are produced from a wide variety of RNA species through identical RNA processing events mainly between cytidine and adenosine.

The identified cP-RNAs are abundantly accumulated in mouse tissues

To biochemically validate that the identified sequences are expressed as cP-RNAs, we focused on the three representative cP-RNAs, 5′-half of tRNAGlyGCC (5′-GlyGCC), cPR-Ncbp3, and cPR-28S (Fig 4B), which are one of the most abundantly sequenced species among tRNA-, mRNA-, and rRNA-derived cP-RNAs, respectively. To analyze 3′-terminal structures of the three species, total RNA from mouse skeletal muscle was treated with CIP or T4 PNK. A mutant T4 PNK, which lacks 3′-dephosphorylation activity [25], was also used as a control. Subsequently, efficiency of the ligations between each target RNA and a 3′-adapter was examined by TaqMan RT-qPCR as described in [26]. As shown in Fig 4C, all three target RNAs exhibited the highest amplification signals upon T4 PNK treatment, indicating that the major 3′-terminal form of the identified sequences are a cP. CIP treatment showed certain signal levels, suggesting co-expression of 3′-P-containing RNAs as a minor form or occurrence of hydrolysis reaction converting a cP to a P during experimental procedures.

Using their specific quantification by TaqMan RT-qPCR, we explored absolute amounts of the three cP-RNAs expressed in mouse lung and skeletal muscle. The calculation of the amounts was based on the standard curve from synthetic RNAs which showed excellent linearity between input amounts and amplification signals (S9 Fig). The determined abundances of the three cP-RNAs per 1 μg of total RNA were summarized in Fig 4D. Whereas the amounts of 5′-GlyGCC and cPR-Ncbp3 were comparable between lung and muscle, the amount of cPR-28S in lung is ~14-fold higher than that in muscle. Considering previous reports estimating that 1 μg of total RNA from mouse liver contains 0.005–8.8 fmol of miRNA [27], the accumulation of cP-RNAs in mouse tissues are estimated to be much more abundant than that of miRNAs.

Age-dependent decrease in the expression levels of cP-RNAs

To explore physiological relevance of cP-RNA expression, TaqMan RT-qPCR quantification of the three cP-RNAs were performed for skeletal muscle of mice with four different ages, 2, 12, 24, and 84 weeks. While a control U6 snRNA showed consistent expression levels in all the examined mice, the expression levels of all three cP-RNAs were reduced along with aging, and the differences between 2- and 84-week old mice were especially apparent (Fig 5A). The cP-RNA levels in 12- and 24-week old mice, both of which are considered as adult mice, did not show a significant change between them and were in-between those of 2- and 84-week old mice. To ask whether the phenomenon is conserved in different tissues, we further investigated lung of mice with different ages, showing the reduction of the expression levels of all three cP-RNAs through aging (Fig 5B). Total RNA staining pattern further suggested global reduction of the expression of 30–45-nt short ncRNAs, which are expected to contain cP-RNAs, through aging (Fig 5C).

Expression of cP-RNAs in tissues from mice with different ages.
Fig. 5. Expression of cP-RNAs in tissues from mice with different ages.
(A) Total RNAs from skeletal muscle of 2–84-week old mice were subjected to TaqMan RT-qPCR for quantification of the indicated cP-RNAs. U6 snRNA was quantified as the control whose expression is unchanged through aging. **P < 0.01, ***P < 0.001 by t-test. (B) Total RNAs from lung of 2–84-week old mice were subjected to cP-RNA quantification. (C) Total RNAs from the biological duplicates of lung from 2- and 84-week old mice were subjected to denaturing PAGE and stained by SYBR Gold. The RNAs in the region designated with a red line were purified and subjected to cP-RNA-seq. (D) cDNAs amplified by cP-RNA-seq. (E) A scatter plot representing the fold changes of the read numbers for each mRNA in 84-week samples, with respect to those in 2-week samples, from cP-RNA-seq (y-axis) and mRNA-seq (x-axis). The plots of the three mRNAs, whose alignment patterns are shown in Fig 5F, are highlighted. (F) The alignment patterns of cP-RNAs for Ncbp3, Cd93, and Tsc22d1 mRNAs. The viewer height in each panel was adjusted based on the normalization by the reads of a spike-in cP-RNA.

To investigate potential relationship between cP-RNAs and their substrate RNAs in expressional regulation through aging, we focused on mRNA-derived cP-RNAs. Total RNAs from lung of 2- and 84-week old mice were subjected to cP-RNA-seq and mRNA-seq for identification of mRNA-derived cP-RNAs and their substrate mRNAs, respectively. Consistent with the downregulated expression of cP-RNAs in older mice (Fig 5A–5C), cP-RNA-seq amplified less abundant cDNAs from the 84-week old mice compared to the 2-week old mice (Fig 5D). Illumina sequencing of the cDNAs yielded ~23–29 million raw reads, of which ~3–6.5 million reads were extracted as the reads of mRNA-derived cP-RNAs (S5 Table). After normalization based on the read number information of a spike-in cP-RNA, fold changes of cP-RNA read numbers of 84 week samples, with respect to those of 2 week samples, were compared with fold changes of mRNA-seq read numbers for each mRNA in a scatter plot (Fig 5E). Strong aggregations of the spots appeared at the point where the levels of cP-RNAs were 4-fold reduced in 84-week samples compared to 2-week samples without changes in the levels of their substrate mRNAs. Consistent with the results of TaqMan RT-qPCR (Fig 5A and 5B), reduction of cPR-Ncbp3 was confirmed in the scatter plot (Fig 5E). As representative results, normalized alignment patterns of cP-RNAs in the three selected mRNAs, Ncbp3, Cd93, and Tsc22d1, are shown in Fig 5F. Specific productions of cP-RNAs through cleavage between C and A residues and their reduction in 84-week samples were observed. Together with the decreased C and A nucleotide biases in 84-week samples (S10 Fig), these results suggest that the expression of cP-RNAs is decreased through aging, which is independent of the levels of substrate RNAs but is caused by reduced C-A cleavage activity for cP-RNA productions or by reduced stability of cP-RNAs.

Conclusion

In this study, we unraveled a previously hidden layer of transcriptome by comprehensively identifying short cP-RNA species expressed in mouse tissues. Each cP-RNA library, obtained by cP-RNA-seq, commonly comprised previously uninterrogated cP-RNA species that are expected to be generated by specific RNA cleavage between cytidine and adenosine (S11 Fig). The endoribonuclease generating the identified cP-RNAs is unknown but should possess highly specific recognition and cleavage activity for the 5′-CA-3′ sequence. ANG might be the candidate enzyme because ANG produces cP-RNAs [10,28] and shows the highest activity toward the 5′-CA-3′ among four examined dinucleotides in a previous study [29]. Among other cP-generating ribonucleases, GCN4 shows high cleavage activity between C and A [30]. Because cP-RNAs are not evenly produced from overall sequences of substrate RNAs but rather derived from specific sites, the RNA cleavage activity to generate cP-RNAs might be exquisitely regulated. Further investigations are required to clarify the enzyme and its activity for the generation of the cP-RNAs identified in this study. Considering the already-proven functional significance of tRNA-derived cP-RNAs [6,7,8,9,10], it is not surprising that the identified cP-RNAs themselves have functional significance. Their abundance and expression from only specific sites of substrate RNAs might also imply their expression as functional molecules. The reduction of cP-RNA generation through aging implied its physiological relevance. It will be intriguing to examine further the cellular factors that influence, or are influenced by, the expression of cP-RNAs and their underlying molecular mechanisms, which will define cP-RNAs as an abundant class of functional non-coding RNAs.

Materials and methods

Ethics statement

All mouse experiments were conducted in compliance with the standards and guidelines of the National Institutes of Health (NIH) and were approved by the Institutional Animal Care and Use Committee at Thomas Jefferson University. Mice were euthanized by CO2 inhalation followed by tissue dissections.

RNA isolation and enzymatic treatment

For cP-RNA-seq, tissues of thymus, lung, spleen, and skeletal muscle from hindlimbs were dissected from male C57BL/6J (B6) mice at approximately 3-months ages. The dissected tissues were homogenized in TRIsure (Bioline) using a glass homogenizer on ice, after which total RNA was extracted according to the manufacturer’s instruction. For TaqMan RT-qPCR, frozen tissues were crushed using a frozen mortar and a pestle, followed by total RNA extraction by TRIsure. To investigate the terminal structures of the 5′-tRNA half, total RNA was treated with CIP (New England Biolabs), T4 PNK (New England Biolabs), or acid (incubation in 10 mM HCl at 4°C for 3 h) followed by CIP treatment as previously described [10,12].

Northern blot

Northern blot was performed as previously described [31] with the following antisense probes: 5′-tRNALysCUU half, 5′- GTCTCATGCTCTACCGACTG-3′; 3′-tRNALysCUU half, 5′-GCTCGAACCCACGACCCTGAGA-3′; 5′-tRNAAspGUC half, 5′-GGGATACTCACCACTATACTAACGAGGA-3′; 3′-tRNAAspGUC half, 5′- GAATCGAACCCCGGTCTCC-3; and miR-16, 5′-GCCAATATTTACGTGCTGCTA-3′.

cP-RNA-seq

For four tissue samples (thymus, lung, spleen, and skeletal muscle), 20 to 45-nt RNAs were gel-purified from total RNA and subjected to cP-RNA-seq procedure as previously described [10,12]. The amplified cDNAs were gel-purified and sequenced using the Illumina HiSeq 2500 system at the Next-Generation Sequencing Core of the University of Pennsylvania.

For comparison between 2- and 84-week old mice, lung 20–45-nt RNAs were gel-purified together with a synthetic spike-in cP-RNA (5′-CAGUGGUGGGCCAGAUGUAAACAUUAGAUUGUUCUUG-cP-3′, synthesized by ChemGenes). cDNAs amplified by cP-RNA-seq procedure were sequenced using Illumina NextSeq 500 at the MetaOmics Core of Sidney Kimmel Cancer Center in Thomas Jefferson University.

Bioinformatics analyses

The cP-RNA-seq libraries contain ~22–47 million raw reads (S1 and S5 Tables) and are publically available at NCBI’s Sequence Read Archive (accession No. SRR8164560, SRR8164561, SRR8164562, SRR8164563, SRR9651811, SRR9651812, SRR9651813, and SRR9651814). Before mapping, we used the cutadapt tool (DOI: http://dx.doi.org/10.14806/ej.17.1.200) to remove the 3′-adapter sequence. After selecting 20–45-nt reads, we used Rbowtie (1.15.1) [32] for the sequential mappings with one mismatch allowed. The reads were first mapped to 471 mature cyto tRNAs obtained from GtRNAdb [19], and then to mature rRNAs, to mitochondrial genome (GenBank DQ106412 sequence plus 22 mitochondrial tRNA sequences), to mRNAs of RefSeq with NM-staring accession numbers, and to whole genome. All the sequences and annotations, except for the mitochondrial genome, were derived from the Mus musculus GRCm38 (mm10) genome assembly. Nucleotide compositions were analyzed by running FastQC (https://www.bioinformatics.babraham.ac.uk/projects/fastqc/). For the nucleotide composition analyses of the sequences which are upstream or downstream of cP-RNA producing regions on the genome, 4-nt upstream/downstream sequences were retrieved by referring fasta files and subjected to FastQC analyses. The Integrative Genomics Viewer (IGV) [33] was used for visualizing the aligned reads.

For comparison between 2- and 84-week samples, the extracted 20–45-nt cP-RNA reads were first mapped to the sequence of the spike-in cP-RNA, followed by the sequential mappings as described above. The read numbers of the spike-in cP-RNA were used to calculate the normalized RPKM of each cP-RNA-producing mRNA. mRNA-seq was performed by Novogene and RPKM of each mRNA was also obtained. These RPKM values were used for calculations of fold-changes of 84-week samples, with respect to 2-week samples, and the values derived from cP-RNA-seq and mRNA-seq were compared in the scatter plot. The experiments of cP-RNA-seq and mRNA-seq were performed using biological duplicates and triplicates, respectively.

Quantification of cP-RNAs by TaqMan RT-qPCR

TaqMan RT-qPCR for specific quantification of cP-RNAs was performed as previously described [10]. Briefly, to remove a cP from cP-RNAs, total RNA was treated with T4 PNK, followed by ligation to a 3′-RNA adapter by T4 RNA ligase. Ligated RNA was then subjected to TaqMan RT-qPCR using the One Step PrimeScript RT-PCR Kit (TaKaRa), 200 nM of a TaqMan probe targeting the boundary of the target RNA and 3′-adapter, and specific forward and reverse primers. The quantified cP-RNA levels were normalized to 5S rRNA levels. The sequences of the TaqMan probes and primers are shown in S6 Table. For the quantification of absolute amounts of cP-RNAs, each synthetic RNA was mixed with 1 μg of E. coli total RNA and the mixtures were used as the substrates to obtain standard curves shown in S9 Fig. In Fig 5A and 5B, Student’s t-tests were used to evaluate the differences between 2- and 84-week samples.

Supporting information

S1 Fig [eps]
RNAs of mouse tissues.

S2 Fig [eps]
Read-length distributions of cP-RNA sequences obtained by cP-RNA-seq.

S3 Fig [eps]
Analyses of tRNA-derived cP-RNAs.

S4 Fig [eps]
Analyses of mRNA-derived cP-RNAs.

S5 Fig [eps]
Analyses of rRNA-derived cP-RNAs.

S6 Fig [eps]
5′- and 3′-terminal nucleotide compositions of the cP-RNAs derived from snRNAs, snoRNAs, miRNAs, and lncRNAs.

S7 Fig [eps]
Nucleotide compositions of the first 30 nucleotides from the 5′-end or the 3′-end of the lncRNA-derived cP-RNAs.

S8 Fig [eps]
Analyses of mitochondrial cP-RNAs.

S9 Fig [eps]
Standard curves for the quantification of cP-RNAs.

S10 Fig [eps]
Nucleotide composition analyses of mRNA-derived cP-RNAs.

S11 Fig [eps]
Schematic representation of the generation of cP-RNAs and its age-dependent changes.

S1 Table [xlsx]
Read numbers of cP-RNA libraries obtained from cP-RNA-seq for four mouse tissues.

S2 Table [xlsx]
Proportion of the 3′-halves (%).

S3 Table [xlsx]
Top 50 reads of mRNA-derived cP-RNAs.

S4 Table [xlsx]
Top 50 reads of rRNA-derived cP-RNAs.

S5 Table [xlsx]
Read numbers of cP-RNA libraries obtained from cP-RNA-seq for lung of 2 and 84 week-old mice.

S6 Table [xlsx]
Sequences of TaqMan probes and primers for RT-qPCR.


Zdroje

1. Goodwin S, McPherson JD, McCombie WR (2016) Coming of age: ten years of next-generation sequencing technologies. Nat Rev Genet 17: 333–351. doi: 10.1038/nrg.2016.49 27184599

2. Marguerat S, Bahler J (2010) RNA-seq: from technology to biology. Cell Mol Life Sci 67: 569–579. doi: 10.1007/s00018-009-0180-6 19859660

3. Costa V, Aprile M, Esposito R, Ciccodicola A (2013) RNA-Seq and human complex diseases: recent accomplishments and future perspectives. Eur J Hum Genet 21: 134–142. doi: 10.1038/ejhg.2012.129 22739340

4. Shigematsu M, Kawamura T, Kirino Y (2018) Generation of 2',3'-Cyclic Phosphate-Containing RNAs as a Hidden Layer of the Transcriptome. Front Genet 9: 562. doi: 10.3389/fgene.2018.00562 30538719

5. Yamasaki S, Ivanov P, Hu GF, Anderson P (2009) Angiogenin cleaves tRNA and promotes stress-induced translational repression. J Cell Biol 185: 35–42. doi: 10.1083/jcb.200811106 19332886

6. Emara MM, Ivanov P, Hickman T, Dawra N, Tisdale S, et al. (2010) Angiogenin-induced tRNA-derived stress-induced RNAs promote stress-induced stress granule assembly. J Biol Chem 285: 10959–10968. doi: 10.1074/jbc.M109.077560 20129916

7. Ivanov P, Emara MM, Villén J, Gygi SP, Anderson P (2011) Angiogenin-Induced tRNA Fragments Inhibit Translation Initiation. Molecular Cell 43: 613–623. doi: 10.1016/j.molcel.2011.06.022 21855800

8. Ivanov P, O'Day E, Emara MM, Wagner G, Lieberman J, et al. (2014) G-quadruplex structures contribute to the neuroprotective effects of angiogenin-induced tRNA fragments. Proc Natl Acad Sci U S A 111: 18201–18206. doi: 10.1073/pnas.1407361111 25404306

9. Blanco S, Dietmann S, Flores JV, Hussain S, Kutter C, et al. (2014) Aberrant methylation of tRNAs links cellular stress to neuro-developmental disorders. EMBO J 33: 2020–2039. doi: 10.15252/embj.201489282 25063673

10. Honda S, Loher P, Shigematsu M, Palazzo JP, Suzuki R, et al. (2015) Sex hormone-dependent tRNA halves enhance cell proliferation in breast and prostate cancers. Proc Natl Acad Sci U S A 112: E3816–3825. doi: 10.1073/pnas.1510077112 26124144

11. Honda S, Kirino Y. (2015) SHOT-RNAs: a novel class of tRNA-derived functional RNAs expressed in hormone-dependent cancers. Mol Cell Onco 3: e1079672. doi: 10.1080/23723556.2015.1079672 27308603

12. Honda S, Morichika K, Kirino Y (2016) Selective amplification and sequencing of cyclic phosphate-containing RNAs by the cP-RNA-seq method. Nat Protoc 11: 476–489. doi: 10.1038/nprot.2016.025 26866791

13. Dhahbi JM, Spindler SR, Atamna H, Yamakawa A, Boffelli D, et al. (2013) 5' tRNA halves are present as abundant complexes in serum, concentrated in blood cells, and modulated by aging and calorie restriction. BMC Genomics 14: 298. doi: 10.1186/1471-2164-14-298 23638709

14. Sharma U, Conine CC, Shea JM, Boskovic A, Derr AG, et al. (2016) Biogenesis and function of tRNA fragments during sperm maturation and fertilization in mammals. Science 351: 391–396. doi: 10.1126/science.aad6780 26721685

15. Shigematsu M, Honda S, Kirino Y (2014) Transfer RNA as a source of small functional RNA. J Mol Biol Mol Imaging 1. 26389128

16. Loher P, Telonis AG, Rigoutsos I (2017) MINTmap: fast and exhaustive profiling of nuclear and mitochondrial tRNA fragments from short RNA-seq data. Sci Rep 7: 41184. doi: 10.1038/srep41184 28220888

17. Pliatsika V, Loher P, Telonis AG, Rigoutsos I (2016) MINTbase: a framework for the interactive exploration of mitochondrial and nuclear tRNA fragments. Bioinformatics 32: 2481–2489. doi: 10.1093/bioinformatics/btw194 27153631

18. Telonis AG, Loher P, Honda S, Jing Y, Palazzo J, et al. (2015) Dissecting tRNA-derived fragment complexities using personalized transcriptomes reveals novel fragment classes and unexpected dependencies. Oncotarget 6: 24797–24822. doi: 10.18632/oncotarget.4695 26325506

19. Chan PP, Lowe TM (2009) GtRNAdb: a database of transfer RNA genes detected in genomic sequence. Nucleic Acids Res 37: D93–97. doi: 10.1093/nar/gkn787 18984615

20. Sprinzl M, Horn C, Brown M, Ioudovitch A, Steinberg S (1998) Compilation of tRNA sequences and sequences of tRNA genes. Nucleic Acids Res 26: 148–153. doi: 10.1093/nar/26.1.148 9399820

21. Shigematsu M, Kirino Y. (2017) 5'-Terminal nucleotide variations in human cytoplasmic tRNAHisGUG and its 5'-halves. RNA 23: 161–168. doi: 10.1261/rna.058024.116 27879434

22. Czech A, Wende S, Mörl M, Pan T, Ignatova Z (2013) Reversible and rapid transfer-RNA deactivation as a mechanism of translational repression in stress. PLoS Genet 9: e1003767. doi: 10.1371/journal.pgen.1003767 24009533

23. Schutz K, Hesselberth JR, Fields S (2010) Capture and sequence analysis of RNAs with terminal 2',3'-cyclic phosphates. RNA 16: 621–631. doi: 10.1261/rna.1934910 20075163

24. Mebius RE, Kraal G (2005) Structure and function of the spleen. Nat Rev Immunol 5: 606–616. doi: 10.1038/nri1669 16056254

25. Wang LK, Shuman S (2002) Mutational analysis defines the 5'-kinase and 3'-phosphatase active sites of T4 polynucleotide kinase. Nucleic Acids Res 30: 1073–1080. doi: 10.1093/nar/30.4.1073 11842120

26. Honda S, Kawamura T, Loher P, Morichika K, Rigoutsos I, et al. (2017) The biogenesis pathway of tRNA-derived piRNAs in Bombyx germ cells. Nucleic Acids Res 45: 9108–9120. doi: 10.1093/nar/gkx537 28645172

27. Bissels U, Wild S, Tomiuk S, Holste A, Hafner M, et al. (2009) Absolute quantification of microRNAs by using a universal reference. RNA 15: 2375–2384. doi: 10.1261/rna.1754109 19861428

28. Shapiro R, Riordan JF, Vallee BL (1986) Characteristic ribonucleolytic activity of human angiogenin. Biochemistry 25: 3527–3532. doi: 10.1021/bi00360a008 2424496

29. Harper JW, Vallee BL (1989) A covalent angiogenin/ribonuclease hybrid with a fourth disulfide bond generated by regional mutagenesis. Biochemistry 28: 1875–1884. doi: 10.1021/bi00430a067 2719939

30. Nikolaev Y, Deillon C, Hoffmann SR, Bigler L, Friess S, et al. (2010) The leucine zipper domains of the transcription factors GCN4 and c-Jun have ribonuclease activity. PLoS One 5: e10765. doi: 10.1371/journal.pone.0010765 20505831

31. Honda S, Kirino Y, Maragkakis M, Alexiou P, Ohtaki A, et al. (2013) Mitochondrial protein BmPAPI modulates the length of mature piRNAs. RNA 19: 1405–1418. doi: 10.1261/rna.040428.113 23970546

32. Langmead B, Trapnell C, Pop M, Salzberg SL (2009) Ultrafast and memory-efficient alignment of short DNA sequences to the human genome. Genome Biol 10: R25. doi: 10.1186/gb-2009-10-3-r25 19261174

33. Robinson JT, Thorvaldsdottir H, Winckler W, Guttman M, Lander ES, et al. (2011) Integrative genomics viewer. Nat Biotechnol 29: 24–26. doi: 10.1038/nbt.1754 21221095

34. Lorenz R, Bernhart SH, Honer Zu Siederdissen C, Tafer H, Flamm C, et al. (2011) ViennaRNA Package 2.0. Algorithms Mol Biol 6: 26. doi: 10.1186/1748-7188-6-26 22115189

Štítky
Genetika Reprodukční medicína

Článek vyšel v časopise

PLOS Genetics


2019 Číslo 11

Nejčtenější v tomto čísle
Kurzy

Zvyšte si kvalifikaci online z pohodlí domova

Hypertenze a hypercholesterolémie – synergický efekt léčby
nový kurz
Autoři: prof. MUDr. Hana Rosolová, DrSc.

Multidisciplinární zkušenosti u pacientů s diabetem
Autoři: Prof. MUDr. Martin Haluzík, DrSc., prof. MUDr. Vojtěch Melenovský, CSc., prof. MUDr. Vladimír Tesař, DrSc.

Úloha kombinovaných preparátů v léčbě arteriální hypertenze
Autoři: prof. MUDr. Martin Haluzík, DrSc.

Halitóza
Autoři: MUDr. Ladislav Korábek, CSc., MBA

Terapie roztroušené sklerózy v kostce
Autoři: MUDr. Dominika Šťastná, Ph.D.

Všechny kurzy
Přihlášení
Zapomenuté heslo

Zadejte e-mailovou adresu, se kterou jste vytvářel(a) účet, budou Vám na ni zaslány informace k nastavení nového hesla.

Přihlášení

Nemáte účet?  Registrujte se

#ADS_BOTTOM_SCRIPTS#