Identification and Functional Analysis of the Pre-piRNA 3' Trimmer in Silkworms.

PIWI-interacting RNAs (piRNAs) play a crucial role in transposon silencing in animal germ cells. In piRNA biogenesis, single-stranded piRNA intermediates are loaded into PIWI-clade proteins and cleaved by Zucchini/MitoPLD, yielding precursor piRNAs (pre-piRNAs). Pre-piRNAs that are longer than the mature piRNA length are then trimmed at their 3' ends. Although recent studies implicated the Tudor domain protein Papi/Tdrkh in pre-piRNA trimming, the identity of Trimmer and its relationship with Papi/Tdrkh remain unknown. Here, we identified PNLDC1, an uncharacterized 3'-5' exonuclease, as Trimmer in silkworms. Trimmer is enriched in the mitochondrial fraction and binds to Papi/Tdrkh. Depletion of Trimmer and Papi/Tdrkh additively inhibits trimming, causing accumulation of ∼35-40-nt pre-piRNAs that are impaired for target cleavage and prone to degradation. Our results highlight the cooperative action of Trimmer and Papi/Tdrkh in piRNA maturation.

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.

The biogenesis pathway of tRNA-derived piRNAs in Bombyx germ cells.

Transfer RNAs (tRNAs) function in translational machinery and further serves as a source of short non-coding RNAs (ncRNAs). tRNA-derived ncRNAs show differential expression profiles and play roles in many biological processes beyond translation. Molecular mechanisms that shape and regulate their expression profiles are largely unknown. Here, we report the mechanism of biogenesis for tRNA-derived Piwi-interacting RNAs (td-piRNAs) expressed in Bombyx BmN4 cells. In the cells, two cytoplasmic tRNA species, tRNAAspGUC and tRNAHisGUG, served as major sources for td-piRNAs, which were derived from the 5'-part of the respective tRNAs. cP-RNA-seq identified the two tRNAs as major substrates for the 5'-tRNA halves as well, suggesting a previously uncharacterized link between 5'-tRNA halves and td-piRNAs. An increase in levels of the 5'-tRNA halves, induced by BmNSun2 knockdown, enhanced the td-piRNA expression levels without quantitative change in mature tRNAs, indicating that 5'-tRNA halves, not mature tRNAs, are the direct precursors for td-piRNAs. For the generation of tRNAHisGUG-derived piRNAs, BmThg1l-mediated nucleotide addition to -1 position of tRNAHisGUG was required, revealing an important function of BmThg1l in piRNA biogenesis. Our study advances the understanding of biogenesis mechanisms and the genesis of specific expression profiles for tRNA-derived ncRNAs.

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.

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.

Analysis of 13 cell types reveals evidence for the expression of numerous novel primate- and tissue-specific microRNAs.

Two decades after the discovery of the first animal microRNA (miRNA), the number of miRNAs in animal genomes remains a vexing question. Here, we report findings from analyzing 1,323 short RNA sequencing samples (RNA-seq) from 13 different human tissue types. Using stringent thresholding criteria, we identified 3,707 statistically significant novel mature miRNAs at a false discovery rate of ≤ 0.05 arising from 3,494 novel precursors; 91.5% of these novel miRNAs were identified independently in 10 or more of the processed samples. Analysis of these novel miRNAs revealed tissue-specific dependencies and a commensurate low Jaccard similarity index in intertissue comparisons. Of these novel miRNAs, 1,657 (45%) were identified in 43 datasets that were generated by cross-linking followed by Argonaute immunoprecipitation and sequencing (Ago CLIP-seq) and represented 3 of the 13 tissues, indicating that these miRNAs are active in the RNA interference pathway. Moreover, experimental investigation through stem-loop PCR of a random collection of newly discovered miRNAs in 12 cell lines representing 5 tissues confirmed their presence and tissue dependence. Among the newly identified miRNAs are many novel miRNA clusters, new members of known miRNA clusters, previously unreported products from uncharacterized arms of miRNA precursors, and previously unrecognized paralogues of functionally important miRNA families (e.g., miR-15/107). Examination of the sequence conservation across vertebrate and invertebrate organisms showed 56.7% of the newly discovered miRNAs to be human-specific whereas the majority (94.4%) are primate lineage-specific. Our findings suggest that the repertoire of human miRNAs is far more extensive than currently represented by public repositories and that there is a significant number of lineage- and/or tissue-specific miRNAs that are uncharacterized.

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.

Mitochondrial protein BmPAPI modulates the length of mature piRNAs.

PIWI proteins and their associated PIWI-interacting RNAs (piRNAs) protect genome integrity by silencing transposons in animal germlines. The molecular mechanisms and components responsible for piRNA biogenesis remain elusive. PIWI proteins contain conserved symmetrical dimethylarginines (sDMAs) that are specifically targeted by TUDOR domain-containing proteins. Here we report that the sDMAs of PIWI proteins play crucial roles in PIWI localization and piRNA biogenesis in Bombyx mori-derived BmN4 cells, which harbor fully functional piRNA biogenesis machinery. Moreover, RNAi screenings for Bombyx genes encoding TUDOR domain-containing proteins identified BmPAPI, a Bombyx homolog of Drosophila PAPI, as a factor modulating the length of mature piRNAs. BmPAPI specifically recognized sDMAs and interacted with PIWI proteins at the surface of the mitochondrial outer membrane. BmPAPI depletion resulted in 3'-terminal extensions of mature piRNAs without affecting the piRNA quantity. These results reveal the BmPAPI-involved piRNA precursor processing mechanism on mitochondrial outer membrane scaffolds.

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.

Concomitant microduplications of MECP2 and ATRX in male patients with severe mental retardation.

Investigations of chromosomal rearrangements in patients with mental retardation (MR) are particularly informative in the search for genes involved in MR. Here we report a family with concomitant duplications of methyl CpG binding protein 2 (MECP2) at Xq28 and ATRX (the causative gene for X-linked alpha thalassemia/mental retardation) at Xq21.1 detected by array-comparative genomic hybridization. The alterations were observed in a 25-year-old man who inherited them from his mother, who showed a normal phenotype and completely skewed X-chromosome inactivation, and also in his cousin, a 32-year-old man. The proband and his cousin showed severe MR, muscular hypotonia, recurrent respiratory infections and various other features characteristic of MECP2 duplication syndrome. However, the proband also had cerebellar atrophy never reported before in MECP2 duplication syndrome, suggesting that his phenotypes were modified through the ATRX duplication in an additive or epistatic manner.

The incidence of hypoplasia of the corpus callosum in patients with dup (X)(q28) involving MECP2 is associated with the location of distal breakpoints.

Duplications of Xq28 harboring the methyl CpG binding protein 2 (MECP2) gene explain approximately 1% of X-linked intellectual disability (XLID). The common clinical features observed in patients with dup(X)(q28) are severe ID, infantile hypotonia, mild dysmorphic features and a history of recurrent infections, and MECP2 duplication syndrome is now recognized as a clinical entity. While some patients with this syndrome have other characteristic phenotypes, the reason for the spectrum of phenotypes has not been clarified. Since dup(X)(q28) rearrangements vary in size and location, genes other than MECP2 might affect the phenotype. We used a high-density oligonucleotide array to carry out precise mapping in eight Japanese families in which dup(X)(q28) was detected using an in-house bacterial artificial chromosome-based microarray to screen cohorts of individuals with multiple congenital anomalies and intellectual disability (MCA/ID) or with XLID. We hypothesized that the size, gene content, and location of dup(X)(q28) may contribute to variable expressively observed in MECP2 duplication syndrome. Genotype-phenotype correlation in our cases together with cases reported in the literature suggested that copy-number gains between two low copy repeats (LCRK1 and LCRL1) are associated with the incidence of hypoplasia of the corpus callosum. Further studies are necessary to understand the mechanism of this association.

Clinical application of array-based comparative genomic hybridization by two-stage screening for 536 patients with mental retardation and multiple congenital anomalies.

Recent advances in the analysis of patients with congenital abnormalities using array-based comparative genome hybridization (aCGH) have uncovered two types of genomic copy-number variants (CNVs); pathogenic CNVs (pCNVs) relevant to congenital disorders and benign CNVs observed also in healthy populations, complicating the screening of disease-associated alterations by aCGH. To apply the aCGH technique to the diagnosis as well as investigation of multiple congenital anomalies and mental retardation (MCA/MR), we constructed a consortium with 23 medical institutes and hospitals in Japan, and recruited 536 patients with clinically uncharacterized MCA/MR, whose karyotypes were normal according to conventional cytogenetics, for two-stage screening using two types of bacterial artificial chromosome-based microarray. The first screening using a targeted array detected pCNV in 54 of 536 cases (10.1%), whereas the second screening of the 349 cases negative in the first screening using a genome-wide high-density array at intervals of approximately 0.7 Mb detected pCNVs in 48 cases (13.8%), including pCNVs relevant to recently established microdeletion or microduplication syndromes, CNVs containing pathogenic genes and recurrent CNVs containing the same region among different patients. The results show the efficient application of aCGH in the clinical setting.

RNase κ promotes robust piRNA production by generating 2',3'-cyclic phosphate-containing precursors.

In animal germlines, PIWI proteins and the associated PIWI-interacting RNAs (piRNAs) protect genome integrity by silencing transposons. Here we report the extensive sequence and quantitative correlations between 2',3'-cyclic phosphate-containing RNAs (cP-RNAs), identified using cP-RNA-seq, and piRNAs in the Bombyx germ cell line and mouse testes. The cP-RNAs containing 5'-phosphate (P-cP-RNAs) identified by P-cP-RNA-seq harbor highly consistent 5'-end positions as the piRNAs and are loaded onto PIWI protein, suggesting their direct utilization as piRNA precursors. We identified Bombyx RNase Kappa (BmRNase κ) as a mitochondria-associated endoribonuclease which produces cP-RNAs during piRNA biogenesis. BmRNase κ-depletion elevated transposon levels and disrupted a piRNA-mediated sex determination in Bombyx embryos, indicating the crucial roles of BmRNase κ in piRNA biogenesis and embryonic development. Our results reveal a BmRNase κ-engaged piRNA biogenesis pathway, in which the generation of cP-RNAs promotes robust piRNA production.

Copy-number variations on the X chromosome in Japanese patients with mental retardation detected by array-based comparative genomic hybridization analysis.

X-linked mental retardation (XLMR) is a common, clinically complex and genetically heterogeneous disease arising from many mutations along the X chromosome. Although research during the past decade has identified >90 XLMR genes, many more remain uncharacterized. In this study, copy-number variations (CNVs) were screened in individuals with MR from 144 families by array-based comparative genomic hybridization (aCGH) using a bacterial artificial chromosome-based X-tiling array. Candidate pathogenic CNVs (pCNVs) were detected in 10 families (6.9%). Five of the families had pCNVs involving known XLMR genes, duplication of Xq28 containing MECP2 in three families, duplication of Xp11.22-p11.23 containing FTSJ1 and PQBP1 in one family, and deletion of Xp11.22 bearing SHROOM4 in one family. New candidate pCNVs were detected in five families as follows: identical complex pCNVs involved in dup(X)(p22.2) and dup(X)(p21.3) containing part of REPS2, NHS and IL1RAPL1 in two unrelated families, duplication of Xp22.2 including part of FRMPD4, duplication of Xq21.1 including HDX and deletion of Xq24 noncoding region in one family, respectively. Both parents and only mother samples were available in six and three families, respectively, and pCNVs were inherited from each of their mothers in those families other than a family of the proband with deletion of SHROOM4. This study should help to identify the novel XLMR genes and mechanisms leading to MR and reveal the clinical conditions and genomic background of XLMR.

Novel deletion at Xq24 including the UBE2A gene in a patient with X-linked mental retardation.

By using an in-house bacterial artificial chromosome-based X-tilling array, we detected a 0.4 Mb novel deletion at Xq24 that included UBE2A in a 4-year-old and 10-month-old boy with mental retardation and various other characteristics inherited from his mother; for example, marked developmental delay, synophrys, ocular hypertelorism, esotropia, low nasal bridge, marked generalized hirsutism and seizure. Although additional nine transcripts around UBE2A were also defective, a phenotypic similarity with a recently reported X-linked familial case involving a novel X-linked mental retardation syndrome and a nonsense mutation of UBE2A indicates a functional defect of UBE2A to be responsible for most of the abnormalities in these cases. Because some characteristics, such as congenital heart disease and proximal placement of the thumb, were not described in the family reported previously, suggesting genes other than UBE2A within the deleted region to be responsible for those abnormalities.

Dumbbell-PCR for Discriminative Quantification of a Small RNA Variant.

Cellular RNAs are often expressed as multiple isoforms of complex heterogeneity in both length and terminal sequences. IsomiRs, the isoforms of microRNAs, are such an example. Distinct quantification of each RNA variant is necessary to unravel the biogenesis mechanism and biological significance of heterogenetic RNA expression. Here we describe Dumbbell-PCR (Db-PCR), a TaqMan RT-PCR-based method that distinctively quantifies a specific small RNA variant with single-nucleotide resolution at terminal sequences. Db-PCR enables the quantitative analysis of RNA terminal heterogeneity without performing Next-Generation Sequencing.

Increasing cell density globally enhances the biogenesis of Piwi-interacting RNAs in Bombyx mori germ cells.

Piwi proteins and their bound Piwi-interacting RNAs (piRNAs) are predominantly expressed in the germline and play crucial roles in germline development by silencing transposons and other targets. Bombyx mori BmN4 cells are culturable germ cells that equip the piRNA pathway. Because of the scarcity of piRNA-expressing culturable cells, BmN4 cells are being utilized for the analyses of piRNA biogenesis. We here report that the piRNA biogenesis in BmN4 cells is regulated by cell density. As cell density increased, the abundance of Piwi proteins and piRNA biogenesis factors was commonly upregulated, resulting in an increased number of perinuclear nuage-like granules where Piwi proteins localize. Along with these phenomena, the abundance of mature piRNAs also globally increased, whereas levels of long piRNA precursor and transposons decreased, suggesting that increasing cell density promotes piRNA biogenesis pathway and that the resultant accumulation of mature piRNAs is functionally significant for transposon silencing. Our study reveals a previously uncharacterized link between cell density and piRNA biogenesis, designates cell density as a critical variable in piRNA studies using BmN4 cell system, and suggests the alteration of cell density as a useful tool to monitor piRNA biogenesis and function.

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.