Transposable element-associated microRNA hairpins produce 21-nt sRNAs integrated into typical microRNA pathways in rice.

microRNAs (miRNAs) are a class of small RNAs (sRNAs) of ~21 nucleotides (nt) in length processed from foldback hairpins by dicer-like1 (DCL1) or DCL4. They regulate the expression of target mRNAs by base pairing through RNA-induced silencing complex (RISC). In the RISC, Argonaute1 (AGO1) is the key protein that cleaves miRNA targets at position ten of a miRNA:target duplex. The authenticity of many annotated rice miRNA hairpins is under debate because of their homology to repeat sequences. Some of them, like miR1884b, have been removed from the current release of miRBase based on incomplete information. In this study, we investigated the association of transposable element (TE)-derived miRNAs with typical miRNA pathways (DCL1/4- and AGO1-dependent) using publicly available deep sequencing datasets. Seven miRNA hairpins with 13 unique sRNAs were specifically enriched in AGO1 immunoprecipitation samples and relatively reduced in DCL1/4 knockdown genotypes. Interestingly, these species are ~21-nt long, instead of 24-nt as annotated in miRBase and the literature. Their expression profiles meet current criteria for functional annotation of miRNAs. In addition, diagnostic cleavage tags were found in degradome datasets for predicted target mRNAs. Most of these miRNA hairpins share significant homology with miniature inverted-repeat transposable elements, one type of abundant DNA transposons in rice. Finally, the root-specific production of a 24-nt miRNA-like sRNA was confirmed by RNA blot for a novel EST that maps to the 3'-UTR of a candidate pseudogene showing extensive sequence homology to miR1884b hairpin. Our data are consistent with the hypothesis that TEs can serve as a driving force for the evolution of some MIRNAs, where co-opting of DICER-LIKE1/4 processing and integration into AGO1 could exapt transcribed TE-associated hairpins into typical miRNA pathways.

An autoregulatory feedback loop involving PAP1 and TAS4 in response to sugars in Arabidopsis.

miR828 in Arabidopsis triggers the cleavage of Trans-Acting SiRNA Gene 4 (TAS4) transcripts and production of small interfering RNAs (ta-siRNAs). One siRNA, TAS4-siRNA81(-), targets a set of MYB transcription factors including PAP1, PAP2, and MYB113 which regulate the anthocyanin biosynthesis pathway. Interestingly, miR828 also targets MYB113, suggesting a close relationship between these MYBs, miR828, and TAS4, but their evolutionary origins are unknown. We found that PAP1, PAP2, and TAS4 expression is induced specifically by exogenous treatment with sucrose and glucose in seedlings. The induction is attenuated in abscisic acid (ABA) pathway mutants, especially in abi3-1 and abi5-1 for PAP1 or PAP2, while no such effect is observed for TAS4. PAP1 is under regulation by TAS4, demonstrated by the accumulation of PAP1 transcripts and anthocyanin in ta-siRNA biogenesis pathway mutants. TAS4-siR81(-) expression is induced by physiological concentrations of Suc and Glc and in pap1-D, an activation-tagged line, indicating a feedback regulatory loop exists between PAP1 and TAS4. Bioinformatic analysis revealed MIR828 homologues in dicots and gymnosperms, but only in one basal monocot, whereas TAS4 is only found in dicots. Consistent with this observation, PAP1, PAP2, and MYB113 dicot paralogs show peptide and nucleotide footprints for the TAS4-siR81(-) binding site, providing evidence for purifying selection in contrast to monocots. Extended sequence similarities between MIR828, MYBs, and TAS4 support an inverted duplication model for the evolution of MIR828 from an ancestral gymnosperm MYB gene and subsequent formation of TAS4 by duplication of the miR828* arm. We obtained evidence by modified 5'-RACE for a MYB mRNA cleavage product guided by miR828 in Pinus resinosa. Taken together, our results suggest that regulation of anthocyanin biosynthesis by TAS4 and miR828 in higher plants is evolutionarily significant and consistent with the evolution of TAS4 since the dicot-monocot divergence.

Evidence for antisense transcription associated with microRNA target mRNAs in Arabidopsis.

Antisense transcription is a pervasive phenomenon, but its source and functional significance is largely unknown. We took an expression-based approach to explore microRNA (miRNA)-related antisense transcription by computational analyses of published whole-genome tiling microarray transcriptome and deep sequencing small RNA (smRNA) data. Statistical support for greater abundance of antisense transcription signatures and smRNAs was observed for miRNA targets than for paralogous genes with no miRNA cleavage site. Antisense smRNAs were also found associated with MIRNA genes. This suggests that miRNA-associated "transitivity" (production of small interfering RNAs through antisense transcription) is more common than previously reported. High-resolution (3 nt) custom tiling microarray transcriptome analysis was performed with probes 400 bp 5' upstream and 3' downstream of the miRNA cleavage sites (direction relative to the mRNA) for 22 select miRNA target genes. We hybridized RNAs labeled from the smRNA pathway mutants, including hen1-1, dcl1-7, hyl1-2, rdr6-15, and sgs3-14. Results showed that antisense transcripts associated with miRNA targets were mainly elevated in hen1-1 and sgs3-14 to a lesser extent, and somewhat reduced in dcl11-7, hyl11-2, or rdr6-15 mutants. This was corroborated by semi-quantitative reverse transcription PCR; however, a direct correlation of antisense transcript abundance in MIR164 gene knockouts was not observed. Our overall analysis reveals a more widespread role for miRNA-associated transitivity with implications for functions of antisense transcription in gene regulation. HEN1 and SGS3 may be links for miRNA target entry into different RNA processing pathways.

RNA structure probing reveals the structural basis of Dicer binding and cleavage.

It is known that an RNA's structure determines its biological function, yet current RNA structure probing methods only capture partial structure information. The ability to measure intact (i.e., full length) RNA structures will facilitate investigations of the functions and regulation mechanisms of small RNAs and identify short fragments of functional sites. Here, we present icSHAPE-MaP, an approach combining in vivo selective 2'-hydroxyl acylation and mutational profiling to probe intact RNA structures. We further showcase the RNA structural landscape of substrates bound by human Dicer based on the combination of RNA immunoprecipitation pull-down and icSHAPE-MaP small RNA structural profiling. We discover distinct structural categories of Dicer substrates in correlation to both their binding affinity and cleavage efficiency. And by tertiary structural modeling constrained by icSHAPE-MaP RNA structural data, we find the spatial distance measuring as an influential parameter for Dicer cleavage-site selection.

Using a barcoded AAV capsid library to select for clinically relevant gene therapy vectors.

While gene transfer using recombinant adeno-associated viral (rAAV) vectors has shown success in some clinical trials, there remain many tissues that are not well transduced. Because of the recent success in reprogramming islet-derived cells into functional ß cells in animal models, we constructed 2 highly complex barcoded replication competent capsid shuffled libraries and selected for high-transducing variants on primary human islets. We describe the generation of a chimeric AAV capsid (AAV-KP1) that facilitates transduction of primary human islet cells and human embryonic stem cell-derived ß cells with up to 10-fold higher efficiency compared with previously studied best-in-class AAV vectors. Remarkably, this chimeric capsid also enabled transduction of both mouse and human hepatocytes at very high levels in a humanized chimeric mouse model, thus providing a versatile vector that has the potential to be used in both preclinical testing and human clinical trials for liver-based diseases and diabetes.

Identification of 20 microRNAs from Oryza sativa.

MicroRNAs (miRNAs) are present in both plant and animal kingdoms and represents a growing family of non-coding RNAs. These tiny RNAs act as small guides and direct negative regulations usually in the process of development through sequence complementarity to target mRNAs. Although a large number of miRNAs have been identified from various animals, so far plant miRNA studies have focused mainly on Arabidopsis. Here we describe the identification of 20 miRNAs from a rice cDNA library. All the miRNAs were presumably processed from precursors with stem-loop structures and were positively detected in rice cells from at least one tissue, some of which showed tissue-specific expression. Twenty-three unique rice genes were identified to be feasible targets for seven rice miRNAs, including four members of Scarecrow-like transcription factor, the targets of miR-39 that had been characterized in Arabidopsis. Lacking long complementarity, the regulatory targets of 13 miRNAs remain to be further investigated. A possible mechanism of translational repressor for plant miRNAs that lack perfect complementarity to target mRNAs is discussed.

Identification of a novel box C/D snoRNA from mouse nucleolar cDNA library.

By construction and screen of mouse nucleolar cDNA library, a novel mammalian small nucleolar RNAs (snoRNA) was identified. The novel snoRNA, 70 nt in length, displays structural features typical of C/D box snoRNA family. The snoRNA possesses an 11-nt-long rRNA antisense element and is predicted to guide the 2'-O-methylation of mouse 28S rRNA at G4043, a site unknown so far to be modified in vertebrates. The comparison of functional element of snoRNA guides among eukaryotes reveals that the novel snoRNA is a mammalian counterpart of yeast snR38 despite highly divergent sequence between them. Mouse and human snR38 and other cognates in distant vertebrates were positively detected with slight length variability. As expected, the rRNA ribose-methylation site predicted by mouse snR38 was precisely mapped by specific-primer extension assay. Furthermore, our analyses show that mouse and human snR38 gene have multiple variants and are nested in the introns of different host genes with unknown function. Thus, snR38 is a phylogenetically conserved methylation guide but exhibits different genomic organization in eukaryotes.

Regulation of microRNA-mediated gene silencing by microRNA precursors.

Processing of microRNAs (miRNAs) from their precursors to their biologically active mature forms is regulated during development and cancer. We show that mouse pri- or pre-miR-151 can bind to and compete with mature miR-151-5p and miR-151-3p for binding sites contained within the complementary regions of the E2f6 mRNA 3' untranslated region (UTR). E2f6 mRNA levels were directly regulated by pri- or pre-miR-151. Conversely, miR-151-mediated repression of ARHGDIA mRNA was dependent on the level of mature miR-151 because only the mature miRNA binds the 3' UTR. Thus, processing of miR-151 can have different effects on separate mRNA targets within a cell. A bioinformatics pipeline revealed additional candidate regions where precursor miRNAs can compete with their mature miRNA counterparts. We validated this experimentally for miR-124 and the SNAI2 3' UTR. Hence, miRNA precursors can serve as post-transcriptional regulators of miRNA activity and are not mere biogenesis intermediates.

Analysis of antisense expression by whole genome tiling microarrays and siRNAs suggests mis-annotation of Arabidopsis orphan protein-coding genes.

BACKGROUND: MicroRNAs (miRNAs) and trans-acting small-interfering RNAs (tasi-RNAs) are small (20-22 nt long) RNAs (smRNAs) generated from hairpin secondary structures or antisense transcripts, respectively, that regulate gene expression by Watson-Crick pairing to a target mRNA and altering expression by mechanisms related to RNA interference. The high sequence homology of plant miRNAs to their targets has been the mainstay of miRNA prediction algorithms, which are limited in their predictive power for other kingdoms because miRNA complementarity is less conserved yet transitive processes (production of antisense smRNAs) are active in eukaryotes. We hypothesize that antisense transcription and associated smRNAs are biomarkers which can be computationally modeled for gene discovery. PRINCIPAL FINDINGS: We explored rice (Oryza sativa) sense and antisense gene expression in publicly available whole genome tiling array transcriptome data and sequenced smRNA libraries (as well as C. elegans) and found evidence of transitivity of MIRNA genes similar to that found in Arabidopsis. Statistical analysis of antisense transcript abundances, presence of antisense ESTs, and association with smRNAs suggests several hundred Arabidopsis 'orphan' hypothetical genes are non-coding RNAs. Consistent with this hypothesis, we found novel Arabidopsis homologues of some MIRNA genes on the antisense strand of previously annotated protein-coding genes. A Support Vector Machine (SVM) was applied using thermodynamic energy of binding plus novel expression features of sense/antisense transcription topology and siRNA abundances to build a prediction model of miRNA targets. The SVM when trained on targets could predict the "ancient" (deeply conserved) class of validated Arabidopsis MIRNA genes with an accuracy of 84%, and 76% for "new" rapidly-evolving MIRNA genes. CONCLUSIONS: Antisense and smRNA expression features and computational methods may identify novel MIRNA genes and other non-coding RNAs in plants and potentially other kingdoms, which can provide insight into antisense transcription, miRNA evolution, and post-transcriptional gene regulation.