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

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.

Exploration of CCA-added RNAs revealed the expression of mitochondrial non-coding RNAs regulated by CCA-adding enzyme.

Post-transcriptional non-template additions of nucleotides to 3'-ends of RNAs play important roles in the stability and function of RNA molecules. Although tRNA nucleotidyltransferase (CCA-adding enzyme) is known to add CCA trinucleotides to 3'-ends of tRNAs, whether other RNA species can be endogenous substrates of CCA-adding enzyme has not been widely explored yet. Herein, we used YAMAT-seq to identify non-tRNA substrates of CCA-adding enzyme. YAMAT-seq captures RNA species that form secondary structures with 4-nt protruding 3'-ends of the sequence 5'-NCCA-3', which is the hallmark structure of RNAs that are generated by CCA-adding enzyme. By executing YAMAT-seq for human breast cancer cells and mining the sequence data, we identified novel candidate substrates of CCA-adding enzyme. These included fourteen 'CCA-RNAs' that only contain CCA as non-genomic sequences, and eleven 'NCCA-RNAs' that contain CCA and other nucleotides as non-genomic sequences. All newly-identified (N)CCA-RNAs were derived from the mitochondrial genome and were localized in mitochondria. Knockdown of CCA-adding enzyme severely reduced the expression levels of (N)CCA-RNAs, suggesting that the CCA-adding enzyme-catalyzed CCA additions stabilize the expression of (N)CCA-RNAs. Furthermore, expression levels of (N)CCA-RNAs were severely reduced by various cellular treatments, including UV irradiation, amino acid starvation, inhibition of mitochondrial respiratory complexes, and inhibition of the cell cycle. These results revealed a novel CCA-mediated regulatory pathway for the expression of mitochondrial non-coding RNAs.

YAMAT-seq: an efficient method for high-throughput sequencing of mature transfer RNAs.

Besides translation, transfer RNAs (tRNAs) play many non-canonical roles in various biological pathways and exhibit highly variable expression profiles. To unravel the emerging complexities of tRNA biology and molecular mechanisms underlying them, an efficient tRNA sequencing method is required. However, the rigid structure of tRNA has been presenting a challenge to the development of such methods. We report the development of Y-shaped Adapter-ligated MAture TRNA sequencing (YAMAT-seq), an efficient and convenient method for high-throughput sequencing of mature tRNAs. YAMAT-seq circumvents the issue of inefficient adapter ligation, a characteristic of conventional RNA sequencing methods for mature tRNAs, by employing the efficient and specific ligation of Y-shaped adapter to mature tRNAs using T4 RNA Ligase 2. Subsequent cDNA amplification and next-generation sequencing successfully yield numerous mature tRNA sequences. YAMAT-seq has high specificity for mature tRNAs and high sensitivity to detect most isoacceptors from minute amount of total RNA. Moreover, YAMAT-seq shows quantitative capability to estimate expression levels of mature tRNAs, and has high reproducibility and broad applicability for various cell lines. YAMAT-seq thus provides high-throughput technique for identifying tRNA profiles and their regulations in various transcriptomes, which could play important regulatory roles in translation and other biological processes.

SHOT-RNAs: A novel class of tRNA-derived functional RNAs expressed in hormone-dependent cancers.

Sex hormones and their receptors play critical roles in the genesis and progression of breast and prostate cancers. We recently discovered that sex hormone signaling pathways promote the expression of specific tRNA halves termed Sex HOrmone-dependent TRNA-derived RNAs (SHOT-RNAs). Functional involvement of SHOT-RNAs in cell proliferation suggests a novel tRNA-engaged pathway in tumorigenesis.

Selective amplification and sequencing of cyclic phosphate-containing RNAs by the cP-RNA-seq method.

RNA digestions catalyzed by many ribonucleases generate RNA fragments that contain a 2',3'-cyclic phosphate (cP) at their 3' termini. However, standard RNA-seq methods are unable to accurately capture cP-containing RNAs because the cP inhibits the adapter ligation reaction. We recently developed a method named cP-RNA-seq that is able to selectively amplify and sequence cP-containing RNAs. Here we describe the cP-RNA-seq protocol in which the 3' termini of all RNAs, except those containing a cP, are cleaved through a periodate treatment after phosphatase treatment; hence, subsequent adapter ligation and cDNA amplification steps are exclusively applied to cP-containing RNAs. cP-RNA-seq takes ∼6 d, excluding the time required for sequencing and bioinformatics analyses, which are not covered in detail in this protocol. Biochemical validation of the existence of cP in the identified RNAs takes ∼3 d. Even though the cP-RNA-seq method was developed to identify angiogenin-generating 5'-tRNA halves as a proof of principle, the method should be applicable to global identification of cP-containing RNA repertoires in various transcriptomes.

Dissecting tRNA-derived fragment complexities using personalized transcriptomes reveals novel fragment classes and unexpected dependencies.

We analyzed transcriptomic data from 452 healthy men and women representing five different human populations and two races, and, 311 breast cancer samples from The Cancer Genome Atlas. Our studies revealed numerous constitutive, distinct fragments with overlapping sequences and quantized lengths that persist across dozens of individuals and arise from the genomic loci of all nuclear and mitochondrial human transfer RNAs (tRNAs). Surprisingly, we discovered that the tRNA fragments' length, starting and ending points, and relative abundance depend on gender, population, race and also on amino acid identity, anticodon, genomic locus, tissue, disease, and disease subtype. Moreover, the length distribution of mitochondrially-encoded tRNAs differs from that of nuclearly-encoded tRNAs, and the specifics of these distributions depend on tissue. Notably, tRNA fragments from the same anticodon do not have correlated abundances. We also report on a novel category of tRNA fragments that significantly contribute to the differences we observe across tissues, genders, populations, and races: these fragments, referred to as i-tRFs, are abundant in human tissues, wholly internal to the respective mature tRNA, and can straddle the anticodon. HITS-CLIP data analysis revealed that tRNA fragments are loaded on Argonaute in a cell-dependent manner, suggesting cell-dependent functional roles through the RNA interference pathway. We validated experimentally two i-tRF molecules: the first was found in 21 of 22 tested breast tumor and adjacent normal samples and was differentially abundant between health and disease whereas the second was found in all eight tested breast cancer cell lines.

Sex hormone-dependent tRNA halves enhance cell proliferation in breast and prostate cancers.

Sex hormones and their receptors play critical roles in the development and progression of the breast and prostate cancers. Here we report that a novel type of transfer RNA (tRNA)-derived small RNA, termed Sex HOrmone-dependent TRNA-derived RNAs (SHOT-RNAs), are specifically and abundantly expressed in estrogen receptor (ER)-positive breast cancer and androgen receptor (AR)-positive prostate cancer cell lines. SHOT-RNAs are not abundantly present in ER(-) breast cancer, AR(-) prostate cancer, or other examined cancer cell lines from other tissues. ER-dependent accumulation of SHOT-RNAs is not limited to a cell culture system, but it also occurs in luminal-type breast cancer patient tissues. SHOT-RNAs are produced from aminoacylated mature tRNAs by angiogenin-mediated anticodon cleavage, which is promoted by sex hormones and their receptors. Resultant 5'- and 3'-SHOT-RNAs, corresponding to 5'- and 3'-tRNA halves, bear a cyclic phosphate (cP) and an amino acid at the 3'-end, respectively. By devising a "cP-RNA-seq" method that is able to exclusively amplify and sequence cP-containing RNAs, we identified the complete repertoire of 5'-SHOT-RNAs. Furthermore, 5'-SHOT-RNA, but not 3'-SHOT-RNA, has significant functional involvement in cell proliferation. These results have unveiled a novel tRNA-engaged pathway in tumorigenesis of hormone-dependent cancers and implicate SHOT-RNAs as potential candidates for biomarkers and therapeutic targets.

Four-leaf clover qRT-PCR: A convenient method for selective quantification of mature tRNA.

Transfer RNAs (tRNAs) play a central role in translation and also recently appear to have a variety of other functions in biological processes beyond translation. Here we report the development of Four-Leaf clover qRT-PCR (FL-PCR), a convenient PCR-based method, which can specifically quantify individual mature tRNA species. In FL-PCR, T4 RNA ligase 2 specifically ligates a stem-loop adapter to mature tRNAs but not to precursor tRNAs or tRNA fragments. Subsequent TaqMan qRT-PCR amplifies only unmodified regions of the tRNA-adapter ligation products; therefore, FL-PCR quantification is not influenced by tRNA post-transcriptional modifications. FL-PCR has broad applicability for the quantification of various tRNAs in different cell types, and thus provides a much-needed simple method for analyzing tRNA abundance and heterogeneity.

Dumbbell-PCR: a method to quantify specific small RNA variants with a single nucleotide resolution at terminal sequences.

Recent advances in next-generation sequencing technologies have revealed that cellular functional RNAs are not always expressed as single entities with fixed terminal sequences but as multiple isoforms bearing complex heterogeneity in both length and terminal sequences, such as isomiRs, the isoforms of microRNAs. Unraveling the biogenesis and biological significance of heterogenetic RNA expression requires distinctive analysis of each RNA variant. Here, we report the development of dumbbell PCR (Db-PCR), an efficient and convenient method to distinctively quantify a specific individual small RNA variant. In Db-PCR, 5'- and 3'-stem-loop adapters are specifically hybridized and ligated to the 5'- and 3'-ends of target RNAs, respectively, by T4 RNA ligase 2 (Rnl2). The resultant ligation products with 'dumbbell-like' structures are subsequently quantified by TaqMan RT-PCR. We confirmed that high specificity of Rnl2 ligation and TaqMan RT-PCR toward target RNAs assured both 5'- and 3'-terminal sequences of target RNAs with single nucleotide resolution so that Db-PCR specifically detected target RNAs but not their corresponding terminal variants. Db-PCR had broad applicability for the quantification of various small RNAs in different cell types, and the results were consistent with those from other quantification method. Therefore, Db-PCR provides a much-needed simple method for analyzing RNA terminal heterogeneity.

PIWI Expression and Function in Cancer.

PIWI proteins, a subclade of the Argonaute family proteins, are expressed predominantly in the germline and bind to PIWI-interacting RNAs (piRNAs), which are 25-31 nucleotides in length. The PIWI/piRNA pathway plays critical roles in germline development by regulating transposons and other targets to maintain genome integrity. While the functions of PIWI in the germline have been extensively investigated, recent studies have accumulated evidence that the human PIWI proteins, HIWI and HILI, are aberrantly expressed in a variety of cancers. This review summarizes our knowledge of PIWI expression in cancer and discusses its possible role in tumorigenesis.