Conserved and highly expressed tRNA derived fragments in zebrafish
© Soares et al. 2015
Received: 27 May 2015
Accepted: 9 December 2015
Published: 22 December 2015
Small non-coding RNAs (sncRNAs) are a class of transcripts implicated in several eukaryotic regulatory mechanisms, namely gene silencing and chromatin regulation. Despite significant progress in their identification by next generation sequencing (NGS) we are still far from understanding their full diversity and functional repertoire.
Here we report the identification of tRNA derived fragments (tRFs) by NGS of the sncRNA fraction of zebrafish. The tRFs identified are 18–30 nt long, are derived from specific 5′ and 3′ processing of mature tRNAs and are differentially expressed during development and in differentiated tissues, suggesting that they are likely produced by specific processing rather than random degradation of tRNAs. We further show that a highly expressed tRF (5′tRF-ProCGG) is cleaved in vitro by Dicer and has silencing ability, indicating that it can enter the RNAi pathway. A computational analysis of zebrafish tRFs shows that they are conserved among vertebrates and mining of publicly available datasets reveals that some 5′tRFs are differentially expressed in disease conditions, namely during infection and colorectal cancer.
tRFs constitute a class of conserved regulatory RNAs in vertebrates and may be involved in mechanisms of genome regulation and in some diseases.
KeywordstRNA derived fragments Zebrafish Small non coding RNAs tRNAs
Small non-coding regulatory RNAs (sncRNAs) play fundamental roles in many aspects of biology and their classification is based on their size, structure, biogenesis and function . MicroRNAs (miRNAs) constitute the most extensively studied class of sncRNAs and are known to regulate the expression of target genes at the translational level . The advent of next generation sequencing (NGS) allowed for the identification and production of miRNA profiles [3–6], identification of piRNAs , 21-U RNAs  and rasiRNAs . More recently, sncRNAs up to 30 nucleotides derived from tRNAs—the tRNA derived fragments (tRFs)—have also been identified [10, 11]. Fragments of tRNA have been retrieved by many computational analysis of NGS datasets, but they were initially considered random degradation products of tRNAs and were discarded from further analysis . However, recent experimental data showed that they are a stable and functional class of sncRNAs produced by specific cleavage of certain tRNAs [10, 11, 13–15], and can be classified according to their origin, namely 5′tRFs, also known as tRF-5 series, which derive from 5′ processing of the mature tRNAs; 3′tRFs, also known as tRF-3 series, which derive from 3′ processing of the mature tRNA and 3′U-tRFs, or tRF-1 series, which derive from pre-tRNA processing by RNAse Z cleavage [11, 16]. Different studies have also demonstrated that some 5′- and 3′tRFs are generated by Dicer cleavage, similarly to miRNAs [10, 17], and that some of them are incorporated into Argonaute (Ago) proteins [13, 18] or are associated with piwi proteins . Indeed, deep sequencing of HeLa cells identified 5′tRFs generated by Dicer cleavage at the D-loop, both in vitro and in vivo . Dicer-dependent generation of 3′tRFs also occurs in a human embryonic kidney 293 cell line and luciferase reporter assays demonstrated that they silence target genes, suggesting that they are involved in gene regulation . Moreover, a Dicer dependent 3′tRF is down-regulated in B cell lymphoma and represses the RPA1 gene, which is involved in DNA repair, similarly to miRNAs . This functional similarity between tRFs and miRNAs may happen because some tRNAs have the potential to form a long hairpin as an alternative to the typical tRNA cloverleaf secondary structure, thus functioning as a Dicer substrate, as is the case of the tRNAIle in mouse embryonic stem cells .
Besides tRFs, longer tRNA fragments (30–38 nt) corresponding to tRNA halves are also produced by cleavage at the 5′ or 3′ sides of the tRNA anticodon by specific endonucleases, rather than Dicer, in response to nutritional deprivation and other stress conditions [20, 21], as is the case in Tetrahymena thermophila during amino acid starvation . The 5′ derived tRNA halves induced by stress (tiRNAs) in human cells have the ability to inhibit translation initiation  and trigger the formation of stress granules , due to a terminal oligoguanine motif , indicating that these molecules are involved in gene expression regulation.
We have previously applied NGS to the discovery of zebrafish sncRNAs and identified novel miRNAs in this organism . A new analysis of the sequencing datasets described herein identified 10 new tRFs that originate from specific cleavage of tRNAs. Expression analysis by northern blot shows that these tRFs are differentially expressed at different developmental stages and in certain tissues and their abundance is higher than the corresponding mature tRNA. Our data show that a 5′tRF, namely 5′tRF-ProCGG can be generated by Dicer and has trans-silencing ability, indicating that it may enter the RNAi pathway, controlling gene expression of specific targets. Northern blot and computational analysis also demonstrate that those tRFs are conserved in vertebrates and are differentially expressed in some disease states, namely during infection and cancer, suggesting their involvement in the mechanisms underlying diseases and their potential use as disease biomarkers.
Identification of tRFs in zebrafish
Zebrafish tRFs identified by deep sequencing
Nr of reads
Development, adult tissues
Development, adult tissues
Development, adult tissues
Development, adult tissues
In zebrafish 12292 genes code for the different tRNAs, according to the genomic tRNA database. Those with highest copy numbers are: Lysine (1478 genes), Glycine (1162 genes), Leucine (1132 genes) and Serine (1065 genes). Although the 10 tRFs identified derived from different classes of tRNAs, there seems to be no correlation between the number of tRNA genes and the generation of tRFs, as most of these abundant tRNAs did not generate tRFs, with exception of Lysine (Fig. 1b). There was a slight enrichment in tRFs generated from Lysine, Glutamine and Proline tRNAs (Fig. 1b), suggesting that the generation of tRFs was not dependent on tRNA expression levels. 5′tRF-LysTTT and 3′tRF-LysCTT derived from different tRNALys isoacceptors, 5′tRF-GluCTC and 3′tRF-GluCTC derived from tRNAGlu, whereas 5′tRF-ProCGG and 3′tRF-ProAGG derived from tRNAPro isoacceptors.
In order to confirm the deep sequencing data, the expression patterns of four tRFs, namely 5′tRF-LysTTT, 5′tRF-GluCTC, 5′tRF-ProCGG and 3′tRF-ProAGG, were studied in different developmental stages [24 h post fertilization (hpf), 48, 72 and 96 hpf] and in different zebrafish adult tissues (brain, fins, bone/muscle, skin) using northern blot analysis, as described in the “Methods”.
The other tRFs tested were not detected by northern blotting during development (Fig. 2a, d), but were detected in adult tissues. The abundance of the 5′tRF-LysTTT was similar in fins, muscle and skin (~1.5 tRF/U6 relative expression), whereas the abundance of the corresponding mature tRNA was higher in muscle (7 tRNA/U6 relative expression) than in other tissues (Fig. 2a). 3′tRF-ProAGG was also detected in adult tissues only, at low levels—maximum 0.3 tRF/U6 relative expression (Fig. 2d).
The lack of correlation between the mature tRNAs and the respective tRFs abundance in different tissues suggested that tRF generation is a regulated process, rather than a random degradation process. Moreover, almost no bands of intermediate molecular weight were detected in the northern blots, indicating that tRNA cleavage occurs at specific cleavage sites. In some tissues, namely brain, fins and skin a band of intermediate size ~30 nt was detected for 5′tRF-ProCGG (Fig. 2c), however the smaller band, corresponding to 5′tRF-ProCGG was always the most prominent one, indicating preferential accumulation of this tRF.
These tRFs are also differentially produced in a variety of tissues (Fig. 3b). The expression of both tRFs in bone is approximately twofold higher than in muscle. In fact, the levels of the 5′tRF-GluCTC were highest in bone and for comparative purposes its relative expression was considered 100 % in this tissue. Its relative expression was also high in skin (60 %) and gut (58 %) and lower in eyes (18 %), brain (15 %) and gills (20 %). This data confirmed the initial profiling data which showed that 5′tRF-GluCTC generation was low in brain and high in skin (Fig. 2). The 5′tRF-ProCGG showed the maximal abundance in gills and for this reason we considered its relative expression as 100 % in this sample. The relative abundance of this tRF was also variable between tissues: eyes (75 %), bone (98 %), skin (70 %) and gut (98 %) and lower in brain and fins (<30 %), confirming the initial profiling data (Figs. 2c, 3b). Therefore, different tRFs are differentially produced and accumulated in the various tissues of zebrafish, in particular both 5′tRF-GluCTC and 5′tRF-ProCGG are generated at lower levels in the brain than in any other tissues tested.
Biogenesis of tRFs
Angiogenin, an RNAseA family member has also been implicated in the generation of some tRNA derived fragments, namely tRNA halves, in response to stress [20, 23, 24]. Another study also shows that angiogenin is involved in the production of 20–25 nt tRFs in vitro in HEK293 cells . Our experiments confirmed that angiogenin is not involved in the generation of the tRFs tested, (Fig. 4b), reinforcing the role of Dicer in the biogenesis of 5′tRF-GluCTC or 5′tRF-ProCGG and indicating that 3′tRF-AlaAGC generation is probably dependent on the action of alternative ribonucleases. Besides we did not detect any increase in 5′tRF-GluCTC or 5′tRF-ProCGG production upon exposure to stress (data not shown), but the tRNAs were cleaved into short non-specific tRNA fragments after 10 min of incubation with Angiogenin (Fig. 4b). Longer incubation times increased the degradation of tRNA and after 2 h the complete fraction of mature tRNA was degraded. Moreover, angiogenin degraded completely the fragments of the tRNA(Glu), but was not able to degrade the 5′tRF-ProCGG, indicating that these tRFs are not angiogenin targets (Fig. 4b).
Recent computational analysis also suggested that Dicer cleavage of tRNAs may occur when the primary transcripts form long hairpins, as an alternative fold to the standard tRNA cloverleaf secondary structure . To test this hypothesis, we determined the probability of the tRNA transcripts forming pre-miRNA like folds using RNAshapes . The Shape Probability option was used to calculate the shape probabilities based on the partition function where the probability of a shape is the sum of the probabilities of all structures that fall into this shape. 3′tRF-AlaAGC can only derive from one tRNA(Ala) locus, while 5′tRF-GluCTC and 5′tRF-ProCGG can be assigned to different tRNA(Glu) and tRNA(Pro) locus in zebrafish (>100 tRNA locus), as found in the genomic tRNA database . We have tested the folding of all of them using RNAShapes and verified that 5 % of tRNA(Glu) and 50 % of tRNA(Pro) transcripts are more prone to acquire a hairpin like folding than any other folding. In most of these cases, the probability of hairpin-like folding is higher than 70 %. Figure 4c shows some of the examples of alternative folds of tRNA(Glu) and tRNA(Pro) obtained, supporting the hypothesis that 5′tRF-GluCTC and 5′tRF-ProCGG can be produced by Dicer. On the other hand the mature tRNA(Ala)AGC showed 98 % probability of forming a typical cloverleaf tRNA-like structure (Fig. 4c).
5′tRF-ProCGG has the ability to silence gene expression
There was a slight decrease in RFP intensity levels after microinjection of DFRS-5′tRF-GluCTC relative to GFP, but the RFP signal was not abolished (Fig. 5m–o). Since 5′tRF-GluCTC levels are low at 24 and 48 hpf we checked for RFP expression after 72 hpf, but no considerable alterations were observed (Additional file 2: Figure S2G–I), indicating that the endogenous 5′tRF-GluCTC does not silence gene expression as efficiently as 5′tRF-ProCGG. Approximately 50 embryos were microinjected per condition and biological replicate. Three biological replicates were performed for each condition tested.
tRFs are conserved between zebrafish and humans
5′tRF levels are affected during infection and in colorectal cancer
We have also analysed the presence of tRFs in human sequencing datasets deposited in the GEO database. From the different datasets analysed, we were able to identify 5′tRF-LysTTT, 5′tRF-ValCAC, 5′tRF-GluCTC, 5′tRF-ProCGG, 3′tRF-GluCTC and 3′tRF-ProAGG in 5 of them, namely GSM1293576 , GSE33584 , GSE29173 , GSE22918  and GSE43550. Low number of reads of 5′tRF-ValCAC and 5′tRF-GluCTC (maximum 185 reads) were identified in serum of old and young individuals (GSM1293576), with no particular age preference (Additional file 3: Table S1). Both tRFs identified were 4 nt longer than the original zebrafish tRFs, probably indicating variation in processing of these molecules in human serum.
3′tRF-ProAGG was the only zebrafish tRF identified in breast tissue (normal and tumour tissue)—GSE29173—, but no significant change in read numbers was found between samples (Additional file 3: Table S1). Only 5′tRF-GluCTC reads were found in the nucleus and in the cytoplasm of 5-8F cells—a nasopharyngeal carcinoma cell line (GSE22918). Reads corresponding to other tRFs, namely 5′tRF-LysTTT, 5′tRF-ProCGG, 3′tRF-GluCTC and 3′tRF-ProAGG were only found in the cytoplasm of 5-8F cells, indicating that most tRFs are present in the cytoplasm rather than in the nucleus (Additional file 3: Table S1). This is not surprising as mature cytoplasmic tRNAs are the precursors of most tRFs, namely 3′ and 5′tRFs . The fact that 5′tRF-GluCTC was found in the nucleus indicates that tRFs migrate to this organelle, similarly to what has already been found for other sncRNAs, namely miRNAs [35, 36].
5′tRF-LysTTT, 5′tRF-ValCAC and 5′tRF-GluCTC were identified in 24 and 72 h non-infected and cytomegalovirus infected primary human fibroblasts (GSE33584). The three identified tRFs were always 2 nt longer when compared to the original zebrafish tRF. Read numbers were relatively low for 5′tRF-LysTTT and 5′tRF-ValCAC, but still 5′tRF-LysTTT reads were fourfold higher in infected cells 72 h post-infection while 5′tRF-ValCAC reads were fourfold lower in 72 h post-infected cells. The number of 5′tRF-GluCTC reads was ~threefold higher 24 h post-infection and ~tenfold higher 72 h post-infection, indicating a role for this tRF during infection (Additional file 3: Table S1; Fig. 7c, d). This particular tRF and 5′tRF-LysTTT, among others, were up-regulated upon respiratory syncytial virus infection in airway epithelial cells (A549 cells) and were implicated in RSV replication . In this case, 5′tRF-GluCTC was only 1 nt longer than the zebrafish 5′tRF-GluCTC and, consequently, 1 nt shorter than the 5′tRF-GluCTC detected in human fibroblasts. Since 5′tRF-GluCTC variability was found in different cell lines it is likely that the 3′end variability is tissue specific and reflect tissue specific variation in processing. 5′tRF-LysTTT had no variation in A549 cells relative to the zebrafish tRF and 5′tRF-GluCTC was abundant after rickettsia infection in mice and exhibited no sequence variation when compared to the equivalent zebrafish tRF .
tRFs are highly expressed in adult zebrafish
In this study we have identified 10 zebrafish tRFs—4 belonging to the tRF-5 series and 6 belonging to the tRF-3 series—distributed throughout development and adult tissues. Since our deep sequencing experiment was designed to identify miRNAs with a NGS read length cutoff between 15 and 30 nt, it is likely that longer tRFs that are generally induced by nutritional deprivation and stress conditions were missed [20, 22, 39]. Studies are needed to clarify this question.
Several tRFs were previously identified in zebrafish early developmental stages, up to 24 hpf . Sequencing of sncRNAs in different developmental stages, namely 256-cell, sphere, shield and 24 hpf showed that tRFs are expressed at low level during early development, but there is a slight enrichment at the 256 cell stage and at 24 hpf . Our analysis started at 24 hpf, when most tRFs are not expressed and progressed to zebrafish tissues where tRFs are abundantly expressed. Most tRFs detected in previous studies  were not identified by us, probably due to the higher sequencing coverage used [26, 37]. Alternatively, those tRFs are present during early development and were missed because we did not analyze the early developmental stages. It is also possible that those tRFs are expressed at exceptionally low levels and were missed by our NGS experiment, as we used a chip that retrieves a maximum of 100,000 reads. Nevertheless, 3′tRF-ProAGG was present in 256-cell, sphere stages and 24 hpf. 5′tRF-ProCGG was also present at 24 hpf . Furthermore, analysis of zebrafish sequencing datasets available in GEO showed that these two tRFs are present in endothelial cells of 24 hpf embryos .
We have found that the cellular concentration of tRFs is not correlated with the abundance of the corresponding tRNAs, which is an indication that these molecules are not products of random tRNA degradation, as previously shown . However, in some cases, namely 5′tRF-ProCGG, an alternative band of approximately 30 nt was detected in adult ZF tissues, indicating that tRF processing may differ across different tissues or developmental stages and thus, different regulation mechanisms of tRF generation may exist and should be further explored in the future.
Northern blot analysis of the tRFs identified by NGS showed that they are highly expressed in mature tissues, such as bone and skin. One of the most striking results was obtained for 5′tRF-ProCGG whose higher abundance in the skin and bone/muscle relative to the corresponding mature tRNA may suggest a role in these particular tissues that needs to be further explored. Its in vitro cleavage by Dicer supports this hypothesis as it allows for incorporation into the RNAi pathway [10, 17, 18], but in vivo cleavage should be tested in the future with zebrafish dicer mutants . How Dicer recognizes these tRNAs still needs to be clarified. One possibility is that the pre-tRNA transcripts may have alternative folds with long hairpins that are recognized and cleaved by Dicer [13, 19]. Coincidently, part of the mature tRNAs that originate 5′tRF-GluCTC and 5′tRF-ProCGG can form long hairpins recognizable by Dicer. It will be interesting to experimentally validate these folding predictions to clarify the role of tRNA folding in the biogenesis of tRFs. Moreover, the folding of tRNAs into alternative secondary and tertiary structures may also explain the existence of tRNA gene duplications in the genomes of vertebrates and the lack of correlation between the abundance of mature tRNAs and their respective tRFs.
Zebrafish tRFs are conserved in vertebrates
We have shown by northern blot analysis that 5′tRF-GluCTC and 5′tRF-ProCGG are conserved in vertebrates, namely in mouse and in human, which is supported by published studies. Indeed, the 5′tRF-GluCTC is described in a human fetus hepatic tissue , a human monocytic cell line (THP-1)  and in adenocarcinomic human alveolar basal epithelial cells (A549) . The tRNAGlu fragment cloned from human fetus hepatic tissue is 3 nt longer than the zebrafish 5′tRF-GluCTC, whereas the 5′tRF-GluCTC identified in THP-1 cells and A549 is 2 nucleotides and 1 nucleotide shorter, respectively, suggesting that the 3′-end of this fragment is heterogeneous, similarly to miRNAs. This may be due to mismatches between the retrieved sequences and the genomic loci, but the cleavage position may also vary with the cell type or tissue [4, 37] due to differential processing . Besides, our analysis of publically available NGS datasets identified several of the zebrafish tRFs described in our study, namely 5′tRF-LysTTT, 5′tRF-ValCAC, 5′tRF-GluCTC, 5′tRF-ProCGG, 3′tRF-GluCTC and 3′tRF-ProAGG. Taken together, our data confirms the conservation of tRFs in vertebrates, suggests conserved tRF processing pathways and tissue specific expression of the tRFs.
5′tRF-ProCGG has silencing ability
Dicer-dependent tRFs share some features of the RNAi pathway and may represent a class of sncRNA molecules with specific functions, including gene silencing. Silencing has been demonstrated with standard reporter assays for a Dicer-dependent tRF  and a recent study showed that a Dicer dependent 3′tRF can repress expression of a set of endogenous genes, including RPA1, by binding to target sites on the mRNA 3′UTR, similarly to miRNAs . Our study shows that 5′tRF-ProCGG functions in the silencing pathway, as the endogenous 5′tRF-ProCGG is able to induce silencing of a RFP reporter by binding to complementary target sites present on its 3′UTR. This miRNA-like feature was further confirmed by loss of silencing of a target gene mutated on the seed region of the tRF (between nucleotides 2 and 8). Since this tRF is highly expressed in bone and skin it will be interesting to validate putative gene targets to understand its functional role. We have predicted computationally putative targets for this tRF, as described before , by maintaining the seed match between nucleotides 2 and 7, and allowing up to 6 mismatches in the remaining sequence. Some of the predicted target genes are involved in embryonic patterning, cartilage and skeletal development, namely sec23b and myst3, which is consistent with the expression pattern of 5′tRF-ProCGG (Additional file 4: Table S2).
Endogenous 5′tRF-GluCTC did not silence efficiently the reporter used in this study, despite its apparent Dicer dependent production in vitro. This could be explained by low levels of endogenous 5′tRF-GluCTC at 24 hpf, however even at 72 hpf there was no obvious RFP silencing, meaning that this fragment does not have miRNA-like trans-silencing ability. Recently 5′tRF-GluCTC was identified in A549 cells after respiratory syncytial virus infection and the authors demonstrated that silencing of a reporter occurred through a mechanism distinct from the miRNA/siRNA pathway . The silencing mechanism was not investigated, but other studies also showed that 5′tRFs can inhibit translation by alternative mechanisms. For example, Sobala and colleagues demonstrated that translational inhibition by 5′tRFs is dependent on a conserved dinucleotide (GG) motif present at position 19 and does not need complementary target sites in the mRNA . These authors also found that endogenous 5′tRFs were not able to silence reporters transfected into cells and that silencing was only observed when a duplex mimic of those tRFs was co-transfected, which is not surprising as the synthetic duplex functions as a siRNA and induces miRNA-like-silencing. General inhibition of translation was achieved when different synthetic 5′tRFs bearing the GG dinucleotide, not in the duplex form, were transfected into cells . A synthetic tRF similar to 5′tRF-GluCTC was the only 5′tRF that did not exhibit silencing ability, despite the GG motif .
tRFs are differentially expressed during infection and in cancer
5′tRF-GluCTC is up-regulated after respiratory syncytial virus  and cytomegalovirus infection  and is down regulated in colon rectal cancer, as shown in the results section. Other tRFs are down regulated in cancer. For example, a 3′tRF (also called CU1276) is down-regulated in B cell lymphoma, similarly to 5′tRF-ValCAC, 5′tRF-LysTTT and 5′tRF-GluCTC in colorectal cancer . These data show a trend for down regulation of tRFs in cancer and up regulation during infection, but further validation experiments are needed to clarify whether these molecules can be used as disease biomarkers. For example, it would be interesting to analyze the tRF profile in different cancers and under different types of infections and also in different populations and experimentally validate the sequencing data by northern blot and qPCR techniques.
Our data show that tRFs are conserved in vertebrates and confirms that tRFs are expressed in a cell and tissue specific manner. For example, in zebrafish the identified tRFs are mostly expressed in muscle, bone and skin and less expressed in the brain and during development, which is probably related to its biological function. Moreover, our data shows that 5′tRF-ProCGG is more expressed in these tissues than its corresponding mature tRNA and that it can play a role in gene expression regulation as it exhibits silencing ability. It is now important to experimentally validate the biological targets of these molecules and determine if those targets are correlated with its high expression. Our data also shows that besides being conserved in vertebrates, tRF expression is affected in specific disease conditions, namely infection and cancer. The differential expression of tRFs in these conditions indicates that tRFs have the potential to be used as disease biomarkers. It would be interesting to analyze the available NGS datasets of human diseases to identify all the tRFs present in specific samples and identify situations where these molecules are deregulated. Their similarity with miRNAs may allow them to recruit the miRNA machinery, but they may have their own machinery or cooperate with other sncRNAs to control important biological processes.
Wild type AB zebrafish strain was obtained from the fish facility at the Department of Biology, University of Aveiro and maintained at 28 °C on a 14 h-light/10 h-dark cycle. Zebrafish maintenance followed the Portuguese law for animal experimentation (Regulatory Guideline no 113/2013, 7th August, 2013) and the experiments were approved by the National Food and Veterinary Authority (DGAV) in Portugal and by the committee for animal experimentation and well-being of the Biology Department, University of Aveiro.
Computational analysis of sequencing reads
Base calling and quality trimming of sequence reads was carried out as described before . For the identification of tRFs, reads that did not match miRNAs were aligned against a small RNA database extracted from Biomart/Ensembl. Up to two mismatches were allowed in alignments to ensure that sequences with sequencing errors or post-transcriptional modifications, which could produce reverse transcription errors during cDNA library construction, were not discarded. tRFs identified were then aligned against their mature tRNA to verify if the sequences were not randomly distributed through the mature sequence.
Total RNA was extracted from zebrafish embryos, adult tissues and cultured cells with TRIZOL® (Invitrogen). Twenty micrograms of total RNA was fractionated on 10 % denaturing polyacrylamide (PAA) gels. RNA was then transferred during 30 min to Hybond-N membranes (GE Healthcare) using a semidry transfer system and was UV-crosslinked. Antisense oligonucleotides complementary to predicted tRF candidates were radio labeled with [32P]-ATP and T4 polynucleotide kinase (Takara) and were used as hybridization probes. Membranes were pre-hybridized for 4 h in pre-hybridization/hybridization buffer containing 5× Denhardt’s solution, 1 % SDS and 6.6× SSPE at 64 °C (5′tRF-LysTTT and 5′tRF-GluCTC), 40 °C (5′tRF-ProCGG), 57 °C (3′tRF-ProAGG and 3′tRF-AlaAGC probe) and 56 °C (U6). Membranes were stripped twice and probed with a maximum of three probes. For the stripping, membranes were incubated with a solution of 50 % formamide and 2xSSPE at 65 °C for 1 h, followed by 15 min incubation with TE buffer to neutralize the membrane. [32P]-ATP labeled probes were added to the hybridization chamber and incubated with the membranes overnight at the mentioned temperatures. Membranes were then washed twice with washing solution, containing 2xSSPE and 0.1 % SDS at room temperature and twice with washing solution at 57 °C for 3 min, and were exposed to a phosphor screen (Biorad®) overnight and scanned using a Molecular Imager® FX (Biorad), equipped with Quantity One FX software.
Angiogenin and dicer in vitro assays
For the Angiogenin cleavage assay, 20 µg of total RNA extracted from 72 h post fertilization (hpf) zebrafish embryos were incubated with 1 µM of the enzyme in PBS + 0.1 % BSA for 10 min, 30 min, 1 and 2 h at 37 °C. The cleaved products were recovered using phenol/chlorophorm extraction followed by ethanol precipitation. Approximately 20 µg of each sample were fractionated on 10 % PAA gels and transferred onto Hybond-N membranes for northern blot analysis as described previously.
For the Dicer assay, recombinant human Dicer enzyme kit from Genlantis was used according to the manufacturer’s instructions with minor changes. Briefly, 20 µg of total RNA extracted from 72 hpf zebrafish embryos were incubated with 2 U of recombinant Dicer in 20 µL reactions, for 30 min, 1, 2, 4 and 6 h. Reactions were stopped with Dicer Stop Solution and were electrophoresed on 10 % PAA gels and transferred onto Hybond-N membranes for northern blot analysis, as described previously.
Dual reporter assay
To test tRFs silencing ability, a dual fluorescence reporter system (DFRS)—pDSV2-eGFP-mRFP was used. This reporter bears two fluorescent proteins, GFP and RFP, controlled by the same promoter. GFP identifies the tissues expressing the plasmid and the RFP, which contains a 3′UTR cassette complementary to the tRF of interest, functions as a silencing sensor. DFRS plasmids were obtained from DR. Wieland Huttner  and slightly modified. Briefly, primers containing multiple cloning sites (XhoI-SaII_NheI), as described in  were annealed and inserted into ECORI-NotI sites of the DFRS plasmid. Next, two different DFRS plasmids containing complementary sites for 5′tRF-GluCTC and 5′tRF-ProCGG and a control DFRS were generated. DFRS-5′tRF-GluCTC, DFRS-5′tRF-ProCGG and DFRS-control were each cloned into XhoI-NheI site by annealing nucleotides A–B, C–D and E–F, respectively.
To study the 5′tRF-ProCGG functional domains DFRS-5′tRF-ProCGG bearing mutations on the 5′ and on the 3′ end were generated. 5′ end mutations of DFRS-5′tRF-ProCGG were obtained after annealing primers G–H. 3′ end mutations were obtained after annealing primers I–J. Mutated sequences were cloned into XhoI-NheI site.
Primers (written in direction 5′–3′):
Zebrafish embryos microinjection
One cell zebrafish eggs were microinjected with ~1000 pL of a solution containing 20 ng/µL of the reporter plasmid, phenol red and 0.9 M KCl. Embryos were kept at 28 °C and analyzed at 24 and 72 hpf under an epifluorescence microscope Imager.Z1 (Zeiss), a GFP and mRFP filter, AxioCam HRm camera (Zeiss) and AxioVision software (Zeiss). Approximately 50 embryos were used per replica and per condition. Three biological replicates were performed.
The 3′UTR sequences of zebrafish mRNAs were extracted from Biomart (http://www.biomart.org) and blasted against the antisense 5′tRF-ProCGG sequence. Sequences with perfect seed match between nucleotides 2 and 7, and no more than six mismatches in the remaining sequence were retained for further analysis. Targets were considered positive whenever RNAhybrid confirmed them thermodynamically. Targets were discarded when RNAhybrid did not retrieve the targets obtained in the first approach.
Computational analysis of publicly available sequencing datasets
Publicly available NGS sequencing datasets deposited in GEO were analyzed to identify the tRFs uncovered in zebrafish. Briefly, datasets were downloaded from GEO and a blast search for each tRF was performed. Lists of identified tRFs and respective number of reads were generated for each sample. Only sequences with more than 50 reads were kept for further analysis. All data used corresponded to the processed sequence along with the cloning frequency that was available in all datasets used. Normalized data against the total number of reads corresponding to sncRNAs for each sample was used to compare the relative levels of tRFs in different samples in the same experiment (sample vs control).
Availability of supporting data
The datasets supporting the results of this article have been deposited in NCBI’s Gene Expression Omnibus  and are accessible through GEO Series accession number GSE66718 (http://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?&acc=GSE66718).
ARS and MASS conceived and design the study. ARS, NF and MR performed the experiments. ARS, HRA, JLO and GM performed the computational analysis of NGS datasets. ARS wrote the manuscript. All authors read and approved the final manuscript.
The authors are most grateful to the Portuguese Foundation for Science and Technology (FCT) for funding our work through project COMPETE/FEDER FCT-ANR/IMI-MIC/0041/2012 and Grant SFRH/BPD/77528/2011 to A.R.S. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
The authors declare that they have no competing interests.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Kim VN, Han J, Siomi MC. Biogenesis of small RNAs in animals. Nat Rev Mol Cell Biol. 2009;10(2):126–39.PubMedView ArticleGoogle Scholar
- Bartel DP. MicroRNAs: genomics, biogenesis, mechanism, and function. Cell. 2004;116(2):281–97.PubMedView ArticleGoogle Scholar
- Burnside J, Ouyang M, Anderson A, Bernberg E, Lu C, Meyers BC, Green PJ, Markis M, Isaacs G, Huang E et al. Deep sequencing of chicken microRNAs. Bmc Genom. 2008;9:185.View ArticleGoogle Scholar
- Landgraf P, Rusu M, Sheridan R, Sewer A, Iovino N, Aravin A, Pfeffer S, Rice A, Kamphorst AO, Landthaler M, et al. A mammalian microRNA expression atlas based on small RNA library sequencing. Cell. 2007;129(7):1401–14.PubMedPubMed CentralView ArticleGoogle Scholar
- Morin RD, O’Connor MD, Griffith M, Kuchenbauer F, Delaney A, Prabhu AL, Zhao Y, McDonald H, Zeng T, Hirst M, et al. Application of massively parallel sequencing to microRNA profiling and discovery in human embryonic stem cells. Genome Res. 2008;18(4):610–21.PubMedPubMed CentralView ArticleGoogle Scholar
- Sunkar R, Zhou XF, Zheng Y, Zhang WX, Zhu JK. Identification of novel and candidate miRNAs in rice by high throughput sequencing. Bmc Plant Biol. 2008;8:25.PubMedPubMed CentralView ArticleGoogle Scholar
- Houwing S, Kamminga LM, Berezikov E, Cronembold D, Girard A, van den Elst H, Filippov DV, Blaser H, Raz E, Moens CB, et al. A role for Piwi and piRNAs in germ cell maintenance and transposon silencing in Zebrafish. Cell. 2007;129(1):69–82.PubMedView ArticleGoogle Scholar
- Ruby JG, Jan C, Player C, Axtell MJ, Lee W, Nusbaum C, Ge H, Bartel DP. Large-scale sequencing reveals 21U-RNAs and additional microRNAs and endogenous siRNAs in C. elegans. Cell. 2006;127(6):1193–207.PubMedView ArticleGoogle Scholar
- Aravin AA, Lagos-Quintana M, Yalcin A, Zavolan M, Marks D, Snyder B, Gaasterland T, Meyer J, Tuschl T. The small RNA profile during Drosophila melanogaster development. Dev Cell. 2003;5(2):337–50.PubMedView ArticleGoogle Scholar
- Cole C, Sobala A, Lu C, Thatcher SR, Bowman A, Brown JW, Green PJ, Barton GJ, Hutvagner G. Filtering of deep sequencing data reveals the existence of abundant Dicer-dependent small RNAs derived from tRNAs. RNA. 2009;15(12):2147–60.PubMedPubMed CentralView ArticleGoogle Scholar
- Lee YS, Shibata Y, Malhotra A, Dutta A. A novel class of small RNAs: tRNA-derived RNA fragments (tRFs). Genes Dev. 2009;23(22):2639–49.PubMedPubMed CentralView ArticleGoogle Scholar
- Kawaji H, Nakamura M, Takahashi Y, Sandelin A, Katayama S, Fukuda S, Daub CO, Kai C, Kawai J, Yasuda J, et al. Hidden layers of human small RNAs. BMC Genom. 2008;9:157.View ArticleGoogle Scholar
- Burroughs AM, Ando Y, Hoon ML, Tomaru Y, Suzuki H, Hayashizaki Y, Daub CO. 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(1):158–77.PubMedPubMed CentralView ArticleGoogle Scholar
- Wang Q, Lee I, Ren J, Ajay SS, Lee YS, Bao X. Identification and functional characterization of tRNA-derived RNA fragments (tRFs) in respiratory syncytial virus infection. Mol Ther J Am Soc Gene Ther. 2013;21(2):368–79.View ArticleGoogle Scholar
- Keam SP, Young PE, McCorkindale AL, Dang TH, Clancy JL, Humphreys DT, Preiss T, Hutvagner G, Martin DI, Cropley JE, et al. The human Piwi protein Hiwi2 associates with tRNA-derived piRNAs in somatic cells. Nucleic Acids Res. 2014;42(14):8984–95.PubMedPubMed CentralView ArticleGoogle Scholar
- Sobala A, Hutvagner G. Small RNAs derived from the 5′ end of tRNA can inhibit protein translation in human cells. RNA Biol. 2013;10(4):553–63.PubMedPubMed CentralView ArticleGoogle Scholar
- Maute RL, Schneider C, Sumazin P, Holmes A, Califano A, Basso K, Dalla-Favera R. tRNA-derived microRNA modulates proliferation and the DNA damage response and is down-regulated in B cell lymphoma. Proc Natl Acad Sci USA. 2013;110(4):1404–9.PubMedPubMed CentralView ArticleGoogle Scholar
- Haussecker D, Huang Y, Lau A, Parameswaran P, Fire AZ, Kay MA. Human tRNA-derived small RNAs in the global regulation of RNA silencing. RNA. 2010;16(4):673–95.PubMedPubMed CentralView ArticleGoogle Scholar
- Babiarz JE, Ruby JG, Wang Y, Bartel DP, Blelloch R. Mouse ES cells express endogenous shRNAs, siRNAs, and other Microprocessor-independent, Dicer-dependent small RNAs. Genes Dev. 2008;22(20):2773–85.PubMedPubMed CentralView ArticleGoogle Scholar
- Fu H, Feng J, Liu Q, Sun F, Tie Y, Zhu J, Xing R, Sun Z, Zheng X. Stress induces tRNA cleavage by angiogenin in mammalian cells. FEBS Lett. 2009;583(2):437–42.PubMedView ArticleGoogle Scholar
- Thompson DM, Parker R. The RNase Rny1p cleaves tRNAs and promotes cell death during oxidative stress in Saccharomyces cerevisiae. J Cell Biol. 2009;185(1):43–50.PubMedPubMed CentralView ArticleGoogle Scholar
- Lee SR, Collins K. Starvation-induced cleavage of the tRNA anticodon loop in Tetrahymena thermophila. J Biol Chem. 2005;280(52):42744–9.PubMedView ArticleGoogle Scholar
- Yamasaki S, Ivanov P, Hu GF, Anderson P. Angiogenin cleaves tRNA and promotes stress-induced translational repression. J Cell Biol. 2009;185(1):35–42.PubMedPubMed CentralView ArticleGoogle Scholar
- Emara MM, Ivanov P, Hickman T, Dawra N, Tisdale S, Kedersha N, Hu GF, Anderson P. Angiogenin-induced tRNA-derived stress-induced RNAs promote stress-induced stress granule assembly. J Biol Chem. 2010;285(14):10959–68.PubMedPubMed CentralView ArticleGoogle Scholar
- Ivanov P, Emara MM, Villen J, Gygi SP, Anderson P. Angiogenin-induced tRNA fragments inhibit translation initiation. Mol Cell. 2011;43(4):613–23.PubMedPubMed CentralView ArticleGoogle Scholar
- Soares AR, Pereira PM, Santos B, Egas C, Gomes AC, Arrais J, Oliveira JL, Moura GR, Santos MAS. Parallel DNA pyrosequencing unveils new zebrafish microRNAs. Bmc Genom. 2009;10:195.View ArticleGoogle Scholar
- Li Z, Ender C, Meister G, Moore PS, Chang Y, John B. Extensive terminal and asymmetric processing of small RNAs from rRNAs, snoRNAs, snRNAs, and tRNAs. Nucleic Acids Res. 2012;40(14):6787–99.PubMedPubMed CentralView ArticleGoogle Scholar
- Steffen P, Voss B, Rehmsmeier M, Reeder J, Giegerich R. RNAshapes: an integrated RNA analysis package based on abstract shapes. Bioinformatics. 2006;22(4):500–3.PubMedView ArticleGoogle Scholar
- Chan PP, Lowe TM. GtRNAdb: a database of transfer RNA genes detected in genomic sequence. Nucleic Acids Res. 2009;37(Database issue):D93–97.
- Noren Hooten N, Fitzpatrick M, Wood WH, 3rd, De S, Ejiogu N, Zhang Y, Mattison JA, Becker KG, Zonderman AB, Evans MK. Age-related changes in microRNA levels in serum. Aging. 2013;5(10):725–40.
- Stark TJ, Arnold JD, Spector DH, Yeo GW. High-resolution profiling and analysis of viral and host small RNAs during human cytomegalovirus infection. J Virol. 2012;86(1):226–35.PubMedPubMed CentralView ArticleGoogle Scholar
- Farazi TA, Horlings HM, Ten Hoeve JJ, Mihailovic A, Halfwerk H, Morozov P, Brown M, Hafner M, Reyal F, van Kouwenhove M, et al. MicroRNA sequence and expression analysis in breast tumors by deep sequencing. Cancer Res. 2011;71(13):4443–53.PubMedPubMed CentralView ArticleGoogle Scholar
- Liao JY, Ma LM, Guo YH, Zhang YC, Zhou H, Shao P, Chen YQ, Qu LH. Deep sequencing of human nuclear and cytoplasmic small RNAs reveals an unexpectedly complex subcellular distribution of miRNAs and tRNA 3′ trailers. PLoS One. 2010;5(5):e10563.PubMedPubMed CentralView ArticleGoogle Scholar
- Phizicky EM, Hopper AK. tRNA biology charges to the front. Genes Dev. 2010;24(17):1832–60.PubMedPubMed CentralView ArticleGoogle Scholar
- Jeffries CD, Fried HM, Perkins DO. Nuclear and cytoplasmic localization of neural stem cell microRNAs. RNA. 2011;17(4):675–86.PubMedPubMed CentralView ArticleGoogle Scholar
- Huang V, Li LC. miRNA goes nuclear. RNA Biol. 2012;9(3):269–73.PubMedPubMed CentralView ArticleGoogle Scholar
- Wei C, Salichos L, Wittgrove CM, Rokas A, Patton JG. Transcriptome-wide analysis of small RNA expression in early zebrafish development. RNA. 2012;18(5):915–29.PubMedPubMed CentralView ArticleGoogle Scholar
- Gong B, Lee YS, Lee I, Shelite TR, Kunkeaw N, Xu G, Lee K, Jeon SH, Johnson BH, Chang Q, et al. Compartmentalized, functional role of angiogenin during spotted fever group rickettsia-induced endothelial barrier dysfunction: evidence of possible mediation by host tRNA-derived small noncoding RNAs. BMC Infect Dis. 2013;13:285.PubMedPubMed CentralView ArticleGoogle Scholar
- Thompson DM, Lu C, Green PJ, Parker R. tRNA cleavage is a conserved response to oxidative stress in eukaryotes. RNA. 2008;14(10):2095–103.PubMedPubMed CentralView ArticleGoogle Scholar
- Nicoli S, Knyphausen CP, Zhu LJ, Lakshmanan A, Lawson ND. miR-221 is required for endothelial tip cell behaviors during vascular development. Dev Cell. 2012;22(2):418–29.PubMedPubMed CentralView ArticleGoogle Scholar
- Giraldez AJ, Cinalli RM, Glasner ME, Enright AJ, Thomson JM, Baskerville S, Hammond SM, Bartel DP, Schier AF. MicroRNAs regulate brain morphogenesis in zebrafish. Science. 2005;308(5723):833–8.PubMedView ArticleGoogle Scholar
- Soares AR, Pereira PM, Santos B, Egas C, Gomes AC, Arrais J, Oliveira JL, Moura GR, Santos MA. Parallel DNA pyrosequencing unveils new zebrafish microRNAs. BMC Genom. 2009;10:195.View ArticleGoogle Scholar
- De Pietri Tonelli D, Calegari F, Fei JF, Nomura T, Osumi N, Heisenberg CP, Huttner WB. Single-cell detection of microRNAs in developing vertebrate embryos after acute administration of a dual-fluorescence reporter/sensor plasmid. Biotechniques. 2006;41(6):727–32.
- Mishima Y, Abreu-Goodger C, Staton AA, Stahlhut C, Shou C, Cheng C, Gerstein M, Enright AJ, Giraldez AJ. Zebrafish miR-1 and miR-133 shape muscle gene expression and regulate sarcomeric actin organization. Genes Dev. 2009;23(5):619–32.PubMedPubMed CentralView ArticleGoogle Scholar
- Edgar R, Domrachev M, Lash AE. Gene expression omnibus: NCBI gene expression and hybridization array data repository. Nucleic Acids Res. 2002;30(1):207–10.PubMedPubMed CentralView ArticleGoogle Scholar