PMCCPMCCPMCC

Search tips
Search criteria 

Advanced

 
Logo of narLink to Publisher's site
 
Nucleic Acids Res. 2017 May 19; 45(9): 5142–5152.
Published online 2017 February 21. doi:  10.1093/nar/gkx103
PMCID: PMC5605234

tRNA-derived small RNAs target transposable element transcripts

Abstract

tRNA-derived RNA fragments (tRFs) are 18–26 nucleotide small RNAs that are not random degradation products, but are rather specifically cleaved from mature tRNA transcripts. Abundant in stressed or viral-infected cells, the function and potential targets of tRFs are not known. We identified that in the unstressed wild-type male gamete containing pollen of flowering plants, and analogous reproductive structure in non-flowering plant species, tRFs accumulate to high levels. In the reference plant Arabidopsis thaliana, tRFs are processed by Dicer-like 1 and incorporated into Argonaute1 (AGO1), akin to a microRNA. We utilized the fact that many plant small RNAs direct cleavage of their target transcripts to demonstrate that the tRF–AGO1 complex acts to specifically target and cleave endogenous transposable element (TE) mRNAs produced from transcriptionally active TEs. The data presented here demonstrate that tRFs are bona-fide regulatory microRNA-like small RNAs involved in the regulation of genome stability through the targeting of TE transcripts.

INTRODUCTION

The development and refinement of high-throughput small RNA (sRNA) sequencing technologies have allowed the field to progressively uncover layers of non-coding RNA (ncRNA) populations involved in the regulation of diverse processes (1,2). The best-studied sRNAs are microRNAs, which are key molecules governing development and stress responses (3). In addition to microRNAs, some sRNAs derived from ncRNAs such as ribosomal RNAs (rRNAs), small nuclear RNAs (snRNAs), small nucleolar RNAs (snoRNAs) and transfer RNAs (tRNAs) are not random degradation products, but rather have specific processing and functions post-processing (46). sRNAs derived from tRNAs have gained particular attention and interest (7). tRNAs, apart from their primordial role in translation, have been suggested to exert a double regulatory function from being part of a distinct control layer through the production of different sRNAs derived from their transcripts (7). Two different classes of sRNAs derived from tRNAs have been reported, tRNA halves of 30–36 nucleotides (nt) and 18–20 nt tRNA-derived RNA fragments (termed tRFs) (8). tRNA halves inhibit protein translation in a wide range of organisms (912), while tRFs have only recently been characterized in depth and their function remains unknown (7). One striking similarity is the accumulation of both tRNA halves and tRFs in stressed cells (7,1316).

In mammals, 18–22 nt tRFs may derive from either the 5΄, 3΄ and tail regions of the precursor or mature tRNA transcripts and have been proposed to be both Dicer-dependent (1719) or -independent (5,6). tRF processing results in the precise generation 18–22 nt sRNAs that are loaded into various Argonaute (AGO) proteins (15,20). These tRF–AGO complexes are biologically active: in human cells they reduce the accumulation of a luciferase reporter construct engineered with a tRF target site (17). Interestingly, tRFs and tRNA halves accumulate to high levels during sperm maturation and affect the expression of genes transgenerationally (12). Because of their Dicer- and AGO-dependent biogenesis pathway, it has been proposed that tRFs could function akin to a microRNA to inhibit the translation or cleave partially-complementary target sites, however the identification of their endogenous targets has yet to be determined (7).

Transposable element (TE) mobilization is detrimental to host fitness as it can cause mutations via insertion, promoter modifications, abnormal genome rearrangements and chromosome breakage (21). In order to restrict their potentially detrimental effects, eukaryotic organisms have evolved multiple mechanisms based on RNA silencing and chromatin modification to control TE expression and maintain TEs in a quiescent silenced state (22). Defeat of host epigenetic silencing and TE transcriptional reactivation occurs upon certain developmental stages, stresses and mutants of various chromatin-modifying proteins (23). One such Arabidopsis mutant is ddm1 (DECREASE IN DNA METHYLATION 1), which lacks the DDM1 swi/snf family chromatin remodeler protein (24). TE transcriptional reactivation also takes place naturally in the Arabidopsis pollen vegetative nucleus, which undergoes a process termed developmental relaxation of TE silencing (DRTS) (23) due to the lack of centromeric and heterochromatin condensation in this cell type (25), which may be linked to the lack of DDM1 protein accumulation (26). In both ddm1 mutants and the pollen grain, transcriptional activation of TEs results in their degradation into small interfering RNAs (siRNAs) as a second-line post-transcriptional defence inhibiting TE activity. RNAi is triggered by sRNAs to cleave the TE mRNAs and initiate RNA-dependent RNA Polymerase (RDR) activity to create double-stranded RNA, which is further Diced into abundant secondary siRNAs (26). However, the identification of the primary (non-RDR-dependent) sRNAs that initiate TE cleavage and RNAi is incomplete.

We utilized the fact that plant sRNAs often cleave their target mRNA transcripts to determine the biogenesis pathway and endogenous targets of tRFs. We found that 19 nt tRFs derived from the 5΄ region of mature tRNA sequences (tRF-5s) accumulate to high levels specifically in pollen. tRF-5s biogenesis involves the core components of the plant microRNA pathway: Dicer-like 1 (DCL1) and AGO1, and we identify that the primary targets of tRFs are TE mRNAs. In addition, we find that the pollen/reproductive tissue-specific accumulation of 19 nt tRF-5s is a characteristic that has been deeply conserved among the plant lineage.

MATERIALS AND METHODS

Plant material

Arabidopsis plants of the ecotype Col-0 were grown under standard long day conditions at 22°C. Inflorescence tissue was used in each experiment unless otherwise noted. Mutant alleles used are described in Supplemental Table S1.

Transgenic constructs

The KRP6-tRF/sRNA target-H2B PCR fragment was cloned as in (27).

Microscopy analysis

Pollen grains of T1 plants were mounted on slides containing 50% glycerol and analyzed under a Nikon-C2 fluorescence microscope. Between four to eight independent transgenic events and 300 and 500 pollen grains per individual transformation event were analysed under the same values of exposure time and fluorescence intensity. GFP intensity in the vegetative nucleus of individual pollen grains was measured through the Nikon-Elements application.

Small RNA deep sequencing and analysis

Inflorescence sRNAs for library preparation were enriched with the mirVana miRNA Isolation Kit (Life Technologies). The sRNAs for wt Col pollen biological replicate GSM1558871 were not size selected before library production, and thus contain 19 nt tRFs. The sRNAs for wt Col pollen biological replicate GSM1558877 and ago1 pollen GSM1558878 were size selected as in (27). sRNA library production and pre-analysis was done as in (27). Mature tRNA sequences were obtained from the tRNA genomic database (http://gtrnadb.ucsc.edu/). tRNA loci and surrounding regions were extracted from the genome version indicated in the tRNA genomic database. sRNA alignments were performed using PatMaN through the UEA sRNA workbench (28). Predictions of tRF target genes were made using the psRNA-Target (29) using a mRNA/tRF pair scoring cutoff of E ≤ 5. PARE-seq libraries were analyzed with PAREsnip through the UEA sRNA workbench (28). The default high stringency settings were used, removing the t/rRNA filtering step. Protein coding genes and TEs identified through this method were further analyzed by aligning PARE-seq reads to their sequences in order to confirm the predicted target sites. The raw sequencing and genome-matched sRNAs created for this study are available from NCBI GEO repository number GSE63879.

Total RNA, sRNA Northern blot, RT-PCR and 5΄RLM RACE-PCR analysis

Total RNA was isolated using TRIzol reagent (Life Technologies). For tRF or microRNA Northern blot detection, 30 or 10 μg of sRNA-enriched RNA were loaded in each lane for inflorescence or pollen Northern blots respectively. Immunoprecipitated sRNA Northern analysis was performed as in (30). sRNA gel electrophoresis, blotting, and cross-linking were performed as described in Pall et al. (31). For mature tRNA and total RNA, 8% polyacrylamide gels were used. cDNA was generated using an oligo-dT or 3΄ region reverse primer and Tetro Reverse Transcriptase (Bioline). Cleavage RACE-PCR was performed as in (30). PCR primers and oligonucleotides used for probe synthesis are listed in Supplemental Table S2.

RESULTS

Specific accumulation of 19 nt 5΄-derived tRFs in the Arabidopsis pollen grain

To characterize ncRNA-derived sRNAs in Arabidopsis, we analyzed the accumulation of ncRNA fragments ranging from 18 to 25 nt in the pollen grain and inflorescence (flower bud) (Figure (Figure1A1A and B). This analysis revealed that both rRNA- and tRNA-derived sRNAs accumulate to high levels in pollen, while, microRNAs decrease in accumulation as previously described (32) (Figure (Figure1A;1A; see Supplemental Figure S1A for analysis of individual bioreplicates). Both sRNAs derived from rRNAs and tRNAs undergo an increase in pollen (rRNA 4.14-fold and tRNA 5.2-fold relative to inflorescence; Figure Figure1B),1B), although the rRNA-derived sRNAs increase is contributed by all sizes of sRNAs, suggesting potential increased degradation and/or transcription of rRNAs (Figure (Figure1A).1A). In contrast, tRF increase is contributed by a single size class of 19 nt sRNAs, suggesting specific processing (Figure (Figure1A1A and Supplemental Figure S1A). The tRFs size profile ranging from 18 to 25 nt in sRNA libraries from leaf, inflorescence and pollen (from additional three bioreps to confirm our initial observation) revealed that 19 nt tRFs specifically accumulate to a higher level in the pollen grain (29-fold higher, Figure Figure1C,1C, Supplemental Figure S1B). We further validated this result by Northern blot, confirming the over accumulation of a 19 nt tRF (from the Ala-AGC tRNA) (Figure (Figure1D1D).

Figure 1.
Characterization of tRFs in Arabidopsis pollen. (A) Analysis of ncRNA-derived sRNAs in inflorescences and pollen. Values are averages of two independent bioreplicates. (B) Analysis of the relative accumulation of ncRNA-derived sRNAs in pollen relative ...

tRFs can derive from both the 5΄ and the 3΄ ends of the mature tRNAs (7). In order to determine the processing features of tRF over-accumulation in pollen, we plotted the 5΄ nucleotide of 19 nt tRFs along tRNA genes and 100 nt on either side of each tRNA coding region. We found that tRFs accumulate from the 5΄ end of the mature tRNA (Supplementary Figure S1C and D). In inflorescence tissue, where few tRFs accumulate (Figure (Figure1C),1C), there is a preference for 5΄ processing of the mature tRNA transcript into tRFs (49%, Supplementary Figure S1C), while in pollen this 5΄ processing preference is increased (88%, Supplementary Figure S1D) indicating that in pollen 5΄-derived tRF biogenesis is more specific and/or efficient. We further checked if tRF specific processing from the 5΄ region is masked by the use of tRNA genomic sequences lacking the CCA sequence that is added to the 3΄ end of tRNAs post-transcriptionally (33), and this analysis revealed that tRFs are only processed from the 5΄ end of mature tRNA sequences (termed tRF-5s, Supplemental Figure S1E). From the 49 tRNAs grouped by codon type in the Arabidopsis genome, tRFs accumulate for 47 (95.92%, Figure Figure1F).1F). Intron presence (SerGCT, MetCAT and TyrGTA), or the number of genes coding for each specific tRNA codon type in the genome, does not influence tRF-5s accumulation (Supplementary Figure S1G). Together, our data suggests that tRF formation is a general phenomenon for nearly all tRNAs in Arabidopsis pollen.

To determine if tRF formation in pollen is found in other plants (34), we investigated published data in rice and maize. We found that the specific 5΄ processing of 19 nt tRF-5s in pollen was conserved in rice and maize (Figure (Figure2A,2A, ,B,B, ,D,D, ,EE and Supplemental Figure 2). Interestingly maize pollen produces a second peak of 24 nt tRFs pointing to the presence of some species-specific processing features (Figure (Figure2B).2B). In order to analyse if the 19 nt tRF accumulation phenomena was a characteristic potentially conserved also in non-flowering land plant precursors, we analysed the tRF distribution in the reproductive tissue of the bryophyte Physcomitrella patens. The Physcomitrella life cycle comprises an alteration of two generations: a filamentous structure termed protonema (a haploid gametophyte) and the gametophore (derived from spores that produce plant-like structures from their buds). These gametophores give rise to the gametes which after fusion will produce a diploid sporophyte (35). The analysis of the sRNA libraries available from gametophore-sporophyte compared to protonemata or protonemata-young gametophore revealed that non-vascular plants accumulate 5΄-derived 19 nt tRFs from their reproductive structures (Figure 2C,D,E and Supplementary Figure 2). Together, these results reveal that throughout the plant kingdom, tRFs undergo specific processing from the 5΄ region of the mature tRNA transcript (tRF-5s) in pollen or analogous male gametophyte structure.

Figure 2.
Accumulation of 19 nt 5΄-derived tRFs is conserved in land plants. (A–C) tRF accumulation profile in pollen and leaf from rice (A) and maize (B), and from gametophore–sporophyte, protonemata and young gametophore–protonometa ...

ddm1 mutants display tRF accumulation

The loss of heterochromatin and activation of TEs in wild-type pollen is mimicked in ddm1 mutant plants (26). This led us to explore whether ddm1 mutants have increased tRF accumulation in the somatic plant body. We found 19 nt tRFs also accumulate to higher levels in inflorescence tissue of ddm1 mutants (4.04 mean fold increase, Figure Figure3A),3A), although this reactivation is lower and more variable between bioreplicates (Supplemental Figure S3A) compared to wild-type pollen. This accumulation is not related with earlier development of the anthers in ddm1 mutants, as 19 nt tRFs also over-accumulate in ddm1 seedlings (Supplemental Figure S3B). The 5΄ preferential processing of tRFs in pollen was also maintained in ddm1 inflorescences (Figure (Figure3B3B and Supplemental Figure S3C). We further confirmed by Northern blot the higher level of accumulation for 19 nt tRF-5s (Figure (Figure3C),3C), indicating the validity of our sRNA high-throughput sequencing analysis. Although the distribution of 19 nt tRF over-accumulation in ddm1 was similar to the one experienced in pollen, fewer tRNAs increase tRF accumulation in ddm1 compared to wt Col inflorescences (Figure (Figure3D3D).

Figure 3.
DDM1 influences the accumulation of tRFs. (A) tRF accumulation profile in ddm1 compared to wt Col inflorescence tissue. The average value of three bioreplicates is shown. (B) Distribution of 5΄ end nucleotide of 19 nt tRFs on tRNA loci in ddm1 ...

Interestingly, we found that increased levels of the corresponding mature tRNA could explain some tRF accumulation in ddm1, while for other tRNAs this correlation did not exist (Figure (Figure3E3E and F). This suggests that some mature tRNA accumulation levels are under the control of, or at least influenced by, DDM1. Other mutants with decreased DNA methylation (met1 and the triple mutant drm1/drm2/cmt3, termed ddc) also have slightly higher accumulation levels of tRFs (Supplemental Figure S3D). Therefore, both mature tRNA levels and tRF accumulation are influenced by epigenetic regulation, and the known reprogramming of heterochromatin and DNA methylation in the pollen grain may contribute to tRF accumulation in this tissue.

tRFs are processed by components of the microRNA pathway

Several studies have shown that human tRF biogenesis depends on the miRNA pathway (18,36), although other authors suggest that DICER could be dispensable for tRF biogenesis (6,34). We find that when tRF-5s accumulate to high levels (in ddm1 and wt pollen) they display similar processing precision values compared to microRNAs, suggesting that tRF-5s were cleaved with a high precision (Supplemental Figure S4A). The observed increased precise processing of 19 nt tRF-5s in ddm1 and pollen led us to examine tRF biogenesis in those tissues. We examined Dicer-like1 (DCL1) due to its known role in microRNA biogenesis (37). We observe reduced tRF accumulation in a ddm1/dcl1 double mutant (Figure (Figure4A4A and Supplemental Figure S4B), demonstrating that the biogenesis of tRFs is similar to microRNAs (38). We next confirmed the accumulation of tRFs by Northern blot in pollen grains from dcl1 and ago1 single mutants (Figure (Figure4B).4B). We found that dcl1 reduces the accumulation of both miRNAs and tRFs (Figure (Figure4B),4B), while defective ago1 only reduces miRNA accumulation, suggesting that AGO1 functions downstream of tRF biogenesis. Furthermore, dcl1 alters to a higher extent the accumulation of tRF-5s in the pollen grain compared to inflorescence (Supplementary Figure S4F–H). This could be indicative of a pollen-exclusive channeling of excessive or reactivated tRNAs to the miRNA pathway.

Figure 4.
5΄-derived tRFs have a microRNA-like biogenesis pathway. (A) tRF accumulation profile in ddm1 and ddm1/dcl1. (B) Northern blot detection of 19 nt AlaAGC tRF-5 in wt, dcl1 and ago1 pollen grains. (C) tRF enrichment profile in wt AGO1 immunoprecipitated ...

Previous bioinformatic analysis of public sRNA libraries reported that 19 nt tRFs are loaded preferentially into AGO1 in wt Col-0 Arabidopsis (34,39). We confirmed that tRFs are incorporated into AGO1 by analysing sRNA libraries constructed from AGO1 immunoprecipitates derived from wt Col and ddm1 inflorescences (Figure (Figure4C4C and Supplemental Figures S4C and D). In our AGO1-immunoprecipitated sRNA libraries, 19 nt tRFs are enriched 9-fold relative to wt inflorescences (Figure (Figure4C),4C), higher than some miRNAs (Figure (Figure4D).4D). In fact, tRFs represented a substantial portion of AGO1-immunoprecipitated sRNAs: 11.24% (wt) and a 17.14% (ddm1 background) of the ncRNA sRNAs (Supplemental Figure S4D). In addition, we verified the presence of tRFs in AGO1 through Northern blot detection of AGO1-immunoprecipitated sRNAs (Figure (Figure4E).4E). The lack of accumulation of this tRF in AGO1-immunoprecipitated sRNAs from dcl1 mutant plants indicates the specific AGO1 loading of DCL1 processed tRFs (Figure (Figure4E),4E), while a mutation in dcl2 showed no effect on tRF accumulation.

Next, we analyzed the cellular localization of tRFs by Northern blot in RNA extracted separately from nuclear and cytoplasmic fractions. 19 nt tRF-5s have a cytoplasmic accumulation pattern (Figure (Figure4F),4F), similar to that of AGO1-loaded microRNAs. In addition, we further confirmed the lack of polyadenylated tRNA transcripts, which could potentially feed into the miRNA biogenesis pathway, through oligo dT-primed RT-PCR in ddm1, ddm1/dcl1 and ddm1/ago1 backgrounds (Supplemental Figure S4E). This data demonstrates that the precursor molecules to tRF-5s are not Pol II-derived transcripts, but rather are mature tRNAs processed by the miRNA pathway (Supplemental Figure S4E).

tRFs target transposable elements

The loading of tRF-5s into AGO1 led us to question if these sRNAs are active components of the RNA silencing network. First, we performed a bioinformatic prediction of potential tRF target-mRNAs, which showed that the majority of these were protein coding genes (78.82%) and TEs (190, 17.72%) (Figure (Figure5A).5A). To check if any of these predicted targets were biologically real, we analyzed the correlation between the predicted 19 nt tRF-5s and mRNA cleavage data from publically available PARE libraries from wt Col-0 inflorescences in a 200 nt window (Supplemental Table S3). Plant microRNAs often cleave their target transcripts between the 10/11 nt of the guide microRNA, providing a mRNA cleavage footprint detectable by PARE-sequencing (40). While no obvious preferential accumulation of PARE reads was detected at the predicted target sites for protein coding genes (Figure (Figure5B,5B, 0.98% of PARE reads at the predicted tRF target site), when the same analysis was carried out for TE sequences, an 8.95-fold increase of TE PARE reads was detected at the tRF target sites (Figure (Figure5C).5C). Interestingly, in ddm1 mutants where both TE mRNAs and tRFs accumulate to higher levels, targeting of TEs by tRFs was amplified 22.9-fold (10 nt surrounding the predicted target site, Figure Figure5D).5D). This increased targeting in ddm1 mutants is exclusive to TEs, as protein-coding gene targeting remains low in ddm1 mutants (Supplemental Figure S5A). As a control we created in silico random 19 nt tRFs not derived from the tRNA 5΄ extreme and found that they do not show the same specific cleavage at the predicted target site as real tRF-5s sequences (Supplemental Figure S5B). In order to check if this targeting was caused by wrong annotation of tRNA genes or presence of tRNA genes within TEs, we analyzed the relative distance of tRNA genes to TEs. Surprisingly, 15 tRNA genes are located within TE sequences in the Arabidopsis genome (Supplemental Figure S5C). From those, 11 correspond to the annotation of a LysCTT tRNA gene within an ATLINEIII element (Supplemental Figure S5D). We found that the TE-located tRNA genes did not induce target cleavage of TEs (Supplemental Table S5), indicating that these intra-TE tRNA genes are not tRF targets.

Figure 5.
tRFs target transposable elements. (A) Categorization of predicted tRF mRNA-targets. (B and C) wt Col inflorescence PARE read alignment along a 100 nt window 5΄ and 3΄ to the predicted target sites for protein coding genes (B) or TEs (C). ...

A majority of tRFs in the ddm1 background targeted TEs belonging to the Gypsy family (83.45%, Figure Figure5E)5E) and have a single tRF target site (81.74%) (Figure (Figure5E).5E). Therefore, our data suggests that when TEs are transcriptionally activated and produce mRNAs, tRFs are able to direct the cleavage of these TE mRNAs. Interestingly, TE sites targeted by tRF-5s are relatively well conserved between Arabidopsis, rice and maize. Both maize and rice share more than 40% of the target sites (up to three mismatches) with Arabidopsis (Supplemental Figure S5E). TE target sites for tRF-5s and miRNAs do not overlap extensively, as only 3% of miRNA targeted TEs are also targeted by tRF-5s (Supplemental Figure S5F). Moreover, the targeted position between miRNAs and tRFs does not overlap in any of the TEs (Supplemental Table S6). This supports a lack of relationship between microRNA and tRF target sites.

Validation of tRF-mediated cleavage of TE mRNAs

To confirm the cleavage ability of tRFs, we selected a 19 nt tRF-5s derived from MetCAT, which accumulates to intermediate levels (Figure (Figure1F)1F) in wt Col-0 pollen and was predicted to target a member of the Gypsy TE family, Athila6A. We analyzed the specific Athila6A cleavage through 5΄ RLM RACE PCR and detected cleavage only in the ddm1 background with a high specificity for the predicted tRF cleavage site (9 out of 10 clones sequenced, Figure Figure6A).6A). In accordance with tRF-5s biogenesis, TE cleavage by the tRF was dependent on DCL1 and partially on AGO1 (Figure (Figure6A).6A). We next generated STTM lines (41) in the ddm1 background (therefore with transcriptionally active TEs) in order to specifically sequester and inhibit the action of the MetCAT 19 nt tRF-5s using as control another transgenic line sequestering a 19 nt tRF-5 derived from AlaAGC (which does not target Athila6A). Analysis of the 5΄RACE product (Figure (Figure6C)6C) and Athila6A transcript accumulation level (Figure (Figure6D)6D) indicated that the targeting of the TE by the tRF was abolished in those transgenic lines specifically sequestering 19 nt MetCAT tRF-5.

Figure 6.
tRF targeting validation. (A) tRF alignment with Athila6A. Lines indicate perfect complementarity, and dotted interactions specify imperfect non-canonical base pairing sites. Cleavage site and frequency of observed/number of 5΄RLM RACE PCR product ...

We next determined if tRF-5s are functional in the wild-type pollen grain, where they accumulate to high levels (Figure (Figure1C).1C). We created a pollen GFP reporter line with a 5΄UTR containing a 19 nt stretch complementary to the MetCAT tRF-5 or a 21 nt ‘mock’ site not targeted by any endogenous siRNA, microRNA or tRF (Figure (Figure6E),6E), and obtained transgenic lines in both wt and dcl1 backgrounds. Analysis of GFP fluorescence indicated that the addition of the 19 nt MetCAT tRF-5 complementary sequence resulted in a statistically significant reduction in pollen fluorescence compared to the control target sequence, while when transformed into the dcl1 mutant background, this difference was abolished (Figure (Figure6E6E and F and Supplementary Figure 6A and B). Together, our data confirms the targeting of TE transcripts by tRNA-derived tRFs (Supplementary Figure S6C).

DISCUSSION

New regulatory functions continue to be discovered for tRNAs, highlighting our lack of complete knowledge of regulatory RNAs and new surprising roles for one of the most well-studied RNA types (42). Here, we have demonstrated that tRNAs are specifically cleaved in plants into 19 nt tRF-5s and that these fragments have miRNA-like biogenesis, localization and function (see model in Supplemental Figure S6C). Although tRF association with RNA silencing has been demonstrated (6), here we present evidence that tRFs are part of a genome protection mechanism that functions through the targeting and cleavage of TE mRNAs, although how this system evolved remains unknown.

Interestingly, several recent papers have linked TEs and tRNAs. Mutations in a subunit of RNaseP (involved in the processing of the 5΄ trailer of tRNA precursors) in Drosophila leads to sterility due to activation of DNA damage checkpoint proteins and TE reactivation coupled with piRNA source collapse (43). Moreover, two miRNAs (miR1511 and miR845, known to target LTR retrotransposons) (44) have precursors that have an evolutionary origin in the tRNA MetCAT (45). Interestingly, here we found that the tRNA MetCAT by itself produces a 19 nt tRF-5 that directly targets the LTR retrotransposon Athila6A. tRF production and TE activity seem to have parallel regulation, as both tRFs and TEs are known be expressed and accumulate during stress (13,14,46). In addition, both tRFs and TE mRNAs accumulate in situations with de-regulated heterochromatin such as ddm1 mutants and in pollen. TE expression is a common phenomenon in wt pollen across flowering plant species (26,47), and we find that tRF accumulation parallels this accumulation in male gametophytes. Interestingly, tRF and tRNA accumulation may be linked to coordinate translation rates. This occurs in human cells (48) and pollen grain germination is dependent on the delayed translation of transcripts accumulated at maturation (49). An excess of mature tRNAs awaiting translation to restart could lead to funnelling of tRNAs in to the miRNA pathway (19). Together, our data suggests a conserved natural tRF-5 biogenesis in the male gametophyte of plants and a potential role in gamete protection against intrinsic gametophyte-associated TE reactivation.

ACCESSION NUMBER

GEO GSE63879.

Supplementary Material

Supplementary Data

ACKNOWLEDGEMENTS

The authors thank Jay Hollick, Jane Jackman and Andrea McCue for helpful comments, and Andrea McCue & Chris DeFraia for their data contributions.

SUPPLEMENTARY DATA

Supplementary Data are available at NAR Online.

FUNDING

National Science Foundation [MCB-1252370 to R.K.S.]; European Commission (Marie Curie IOF grant) [PIOF-GA-2012-330069 to G.M.]. Funding for open access charge: Swedish University of Agricultural Sciences.

Conflict of interest statement. None declared.

REFERENCES

1. Martens-Uzunova E.S., Olvedy M., Jenster G. Beyond microRNA–novel RNAs derived from small non-coding RNA and their implication in cancer. Cancer Lett. 2013; 340:201–211. [PubMed]
2. Kawaji H., Nakamura M., Takahashi Y., Sandelin A., Katayama S., Fukuda S., Daub C.O., Kai C., Kawai J., Yasuda J. et al. Hidden layers of human small RNAs. BMC Genomics. 2008; 9:157. [PMC free article] [PubMed]
3. Jones-Rhoades M.W., Bartel D.P., Bartel B. MicroRNAS and their regulatory roles in plants. Annu. Rev. Plant Biol. 2006; 57:19–53. [PubMed]
4. Ender C., Krek A., Friedlander M.R., Beitzinger M., Weinmann L., Chen W., Pfeffer S., Rajewsky N., Meister G. A human snoRNA with microRNA-like functions. Mol. Cell. 2008; 32:519–528. [PubMed]
5. Li Z., Ender C., Meister G., Moore P.S., Chang Y., John B. Extensive terminal and asymmetric processing of small RNAs from rRNAs, snoRNAs, snRNAs, and tRNAs. Nucleic Acids Res. 2012; 40:6787–6799. [PMC free article] [PubMed]
6. Kumar P., Anaya J., Mudunuri S.B., Dutta A. Meta-analysis of tRNA derived RNA fragments reveals that they are evolutionarily conserved and associate with AGO proteins to recognize specific RNA targets. BMC Biol. 2014; 12:78. [PMC free article] [PubMed]
7. Pederson T. Regulatory RNAs derived from transfer RNA. RNA. 2010; 16:1865–1869. [PubMed]
8. Dhahbi J.M. 5΄ tRNA halves: the next generation of immune signaling molecules. Front. Immunol. 2015; 6:74. [PMC free article] [PubMed]
9. Zhang S., Sun L., Kragler F. The phloem-delivered RNA pool contains small noncoding RNAs and interferes with translation. Plant Physiol. 2009; 150:378–387. [PubMed]
10. Ivanov P., Emara M.M., Villen J., Gygi S.P., Anderson P. Angiogenin-induced tRNA fragments inhibit translation initiation. Mol. Cell. 2011; 43:613–623. [PMC free article] [PubMed]
11. Gebetsberger J., Zywicki M., Kunzi A., Polacek N. tRNA-derived fragments target the ribosome and function as regulatory non-coding RNA in Haloferax volcanii. Archaea. 2012; 2012:260909. [PMC free article] [PubMed]
12. Sharma U., Conine C.C., Shea J.M., Boskovic A., Derr A.G., Bing X.Y., Belleannee C., Kucukural A., Serra R.W., Sun F. Biogenesis and function of tRNA fragments during sperm maturation and fertilization in mammals. Science. 2015; 351:391–396. [PMC free article] [PubMed]
13. Hackenberg M., Huang P.J., Huang C.Y., Shi B.J., Gustafson P., Langridge P. A comprehensive expression profile of microRNAs and other classes of non-coding small RNAs in barley under phosphorous-deficient and -sufficient conditions. DNA Res. 2013; 20:109–125. [PMC free article] [PubMed]
14. Hsieh L.C., Lin S.I., Shih A.C., Chen J.W., Lin W.Y., Tseng C.Y., Li W.H., Chiou T.J. Uncovering small RNA-mediated responses to phosphate deficiency in Arabidopsis by deep sequencing. Plant Physiol. 2009; 151:2120–2132. [PubMed]
15. Wang Q., Lee I., Ren J., Ajay S.S., Lee Y.S., Bao X. Identification and functional characterization of tRNA-derived RNA fragments (tRFs) in respiratory syncytial virus infection. Mol. Ther. 2013; 21:368–379. [PubMed]
16. Garcia-Silva M.R., Cabrera-Cabrera F., Guida M.C., Cayota A. Hints of tRNA-derived small RNAs role in RNA silencing mechanisms. Genes (Basel). 2012; 3:603–614. [PMC free article] [PubMed]
17. Yeung M.L., Bennasser Y., Watashi K., Le S.Y., Houzet L., Jeang K.T. Pyrosequencing of small non-coding RNAs in HIV-1 infected cells: evidence for the processing of a viral-cellular double-stranded RNA hybrid. Nucleic Acids Res. 2009; 37:6575–6586. [PMC free article] [PubMed]
18. Cole C., Sobala A., Lu C., Thatcher S.R., Bowman A., Brown J.W., Green P.J., Barton G.J., Hutvagner G. Filtering of deep sequencing data reveals the existence of abundant Dicer-dependent small RNAs derived from tRNAs. RNA. 2009; 15:2147–2160. [PubMed]
19. Hasler D., Lehmann G., Murakawa Y., Klironomos F., Jakob L., Grasser F.A., Rajewsky N., Landthaler M., Meister G. The lupus autoantigen La prevents mis-channeling of tRNA fragments into the human microRNA pathway. Mol. Cell. 2016; 63:110–124. [PubMed]
20. Burroughs A.M., Ando Y., de Hoon M.J., Tomaru Y., Suzuki H., Hayashizaki Y., Daub C.O. Deep-sequencing of human Argonaute-associated small RNAs provides insight into miRNA sorting and reveals Argonaute association with RNA fragments of diverse origin. RNA Biol. 2011; 8:158–177. [PMC free article] [PubMed]
21. Lisch D., Bennetzen J.L. Transposable element origins of epigenetic gene regulation. Curr. Opin. Plant Biol. 2011; 14:156–161. [PubMed]
22. Slotkin R.K., Martienssen R. Transposable elements and the epigenetic regulation of the genome. Nat. Rev. Genet. 2007; 8:272–285. [PubMed]
23. Martinez G., Slotkin R.K. Developmental relaxation of transposable element silencing in plants: functional or byproduct. Curr. Opin. Plant Biol. 2012; 15:496–502. [PubMed]
24. Zemach A., Kim M.Y., Hsieh P.H., Coleman-Derr D., Eshed-Williams L., Thao K., Harmer S.L., Zilberman D. The Arabidopsis nucleosome remodeler DDM1 allows DNA methyltransferases to access H1-containing heterochromatin. Cell. 2013; 153:193–205. [PMC free article] [PubMed]
25. Schoft V.K., Chumak N., Mosiolek M., Slusarz L., Komnenovic V., Brownfield L., Twell D., Kakutani T., Tamaru H. Induction of RNA-directed DNA methylation upon decondensation of constitutive heterochromatin. EMBO Rep. 2009; 10:1015–1021. [PubMed]
26. Slotkin R.K., Vaughn M., Borges F., Tanurdzic M., Becker J.D., Feijo J.A., Martienssen R.A. Epigenetic reprogramming and small RNA silencing of transposable elements in pollen. Cell. 2009; 136:461–472. [PMC free article] [PubMed]
27. Martinez G., Panda K., Kohler C., Slotkin R.K. Silencing in sperm cells is directed by RNA movement from the surrounding nurse cell. Nat. Plants. 2016; 2:16030. [PubMed]
28. Stocks M.B., Moxon S., Mapleson D., Woolfenden H.C., Mohorianu I., Folkes L., Schwach F., Dalmay T., Moulton V. The UEA sRNA workbench: a suite of tools for analysing and visualizing next generation sequencing microRNA and small RNA datasets. Bioinformatics. 2012; 28:2059–2061. [PMC free article] [PubMed]
29. Dai X., Zhao P.X. psRNATarget: a plant small RNA target analysis server. Nucleic Acids Res. 2011; 39:W155–W159. [PMC free article] [PubMed]
30. McCue A.D., Nuthikattu S., Slotkin R.K. Genome-wide identification of genes regulated in trans by transposable element small interfering RNAs. RNA Biol. 2013; 10:1379–1395. [PMC free article] [PubMed]
31. Pall G.S., Hamilton A.J. Improved northern blot method for enhanced detection of small RNA. Nat. Protoc. 2008; 3:1077–1084. [PubMed]
32. Borges F., Pereira P.A., Slotkin R.K., Martienssen R.A., Becker J.D. MicroRNA activity in the Arabidopsis male germline. J. Exp. Bot. 2011; 62:1611–1620. [PMC free article] [PubMed]
33. Hopper A.K. Transfer RNA post-transcriptional processing, turnover, and subcellular dynamics in the yeast Saccharomyces cerevisiae. Genetics. 2013; 194:43–67. [PubMed]
34. Alves C.S., Vicentini R., Duarte G.T., Pinoti V.F., Vincentz M., Nogueira F.T. Genome-wide identification and characterization of tRNA-derived RNA fragments in land plants. Plant Mol. Biol. 2016; 93:35–48. [PubMed]
35. Schaefer D.G., Zryd J.P. The moss Physcomitrella patens, now and then. Plant Physiol. 2001; 127:1430–1438. [PubMed]
36. Haussecker D., Huang Y., Lau A., Parameswaran P., Fire A.Z., Kay M.A. Human tRNA-derived small RNAs in the global regulation of RNA silencing. RNA. 2010; 16:673–695. [PubMed]
37. Baulcombe D. RNA silencing in plants. Nature. 2004; 431:356–363. [PubMed]
38. Vaucheret H., Vazquez F., Crete P., Bartel D.P. The action of ARGONAUTE1 in the miRNA pathway and its regulation by the miRNA pathway are crucial for plant development. Genes Dev. 2004; 18:1187–1197. [PubMed]
39. Loss-Morais G., Waterhouse P.M., Margis R. Description of plant tRNA-derived RNA fragments (tRFs) associated with argonaute and identification of their putative targets. Biol. Direct. 2013; 8:6. [PMC free article] [PubMed]
40. German M.A., Pillay M., Jeong D.H., Hetawal A., Luo S., Janardhanan P., Kannan V., Rymarquis L.A., Nobuta K., German R. et al. Global identification of microRNA-target RNA pairs by parallel analysis of RNA ends. Nat. Biotechnol. 2008; 26:941–946. [PubMed]
41. Yan J., Gu Y., Jia X., Kang W., Pan S., Tang X., Chen X., Tang G. Effective small RNA destruction by the expression of a short tandem target mimic in Arabidopsis. Plant Cell. 2012; 24:415–427. [PubMed]
42. Stoecklin G., Diederichs S. tRNAs: new tricks from old dogs. EMBO J. 2014; 33:1981–1983. [PubMed]
43. Molla-Herman A., Valles A.M., Ganem-Elbaz C., Antoniewski C., Huynh J.R. tRNA processing defects induce replication stress and Chk2-dependent disruption of piRNA transcription. EMBO J. 2015; 34:3009–3027. [PubMed]
44. Creasey K.M., Zhai J., Borges F., Van Ex F., Regulski M., Meyers B.C., Martienssen R.A. miRNAs trigger widespread epigenetically activated siRNAs from transposons in Arabidopsis. Nature. 2014; 508:411–415. [PMC free article] [PubMed]
45. Surbanovski N., Brilli M., Moser M., Si-Ammour A. A highly specific microRNA-mediated mechanism silences LTR retrotransposons of strawberry. Plant J. 2016; 85:70–82. [PubMed]
46. McCue A.D., Nuthikattu S., Reeder S.H., Slotkin R.K. Gene expression and stress response mediated by the epigenetic regulation of a transposable element small RNA. PLoS Genet. 2012; 8:e1002474. [PMC free article] [PubMed]
47. Nobuta K., Venu R.C., Lu C., Belo A., Vemaraju K., Kulkarni K., Wang W., Pillay M., Green P.J., Wang G.L. et al. An expression atlas of rice mRNAs and small RNAs. Nat. Biotechnol. 2007; 25:473–477. [PubMed]
48. Gingold H., Tehler D., Christoffersen N.R., Nielsen M.M., Asmar F., Kooistra S.M., Christophersen N.S., Christensen L.L., Borre M., Sorensen K.D. et al. A dual program for translation regulation in cellular proliferation and differentiation. Cell. 2014; 158:1281–1292. [PubMed]
49. Honys D., Combe P.J., Twell D., Capková V. The translationally repressed pollen-specific ntp303 mRNA is stored in non-polysomal mRNPs during pollen maturation. Sexual Plant Reprod. 2000; 13:135–144.

Articles from Nucleic Acids Research are provided here courtesy of Oxford University Press