A Comprehensive Map of Intron Branchpoints and Lariat RNAs in Plants.

Lariats are formed by excised introns, when the 5' splice site joins with the branchpoint (BP) during splicing. Although lariat RNAs are usually degraded by RNA debranching enzyme 1, recent findings in animals detected many lariat RNAs under physiological conditions. By contrast, the features of BPs and to what extent lariat RNAs accumulate naturally are largely unexplored in plants. Here, we analyzed 948 RNA sequencing data sets to document plant BPs and lariat RNAs on a genome-wide scale. In total, we identified 13,872, 5199, 29,582, and 13,478 BPs in Arabidopsis (Arabidopsis thaliana), tomato (Solanum lycopersicum), rice (Oryza sativa), and maize (Zea mays), respectively. Features of plant BPs are highly similar to those in yeast and human, in that BPs are adenine-preferred and flanked by uracil-enriched sequences. Intriguingly, ∼20% of introns harbor multiple BPs, and BP usage is tissue-specific. Furthermore, 10,580 lariat RNAs accumulate in wild-type Arabidopsis plants, and most of these lariat RNAs originate from longer or retroelement-depleted introns. Moreover, the expression of these lariat RNAs is accompanied by the incidence of back-splicing of parent exons. Collectively, our results provide a comprehensive map of intron BPs and lariat RNAs in four plant species and uncover a link between lariat turnover and splicing.

Anaphase-promoting complex/cyclosome regulates RdDM activity by degrading DMS3 in Arabidopsis.

During RNA-directed DNA methylation (RdDM), the DDR complex, composed of DRD1, DMS3, and RDM1, is responsible for recruiting DNA polymerase V (Pol V) to silence transposable elements (TEs) in plants. However, how the DDR complex is regulated remains unexplored. Here, we show that the anaphase-promoting complex/cyclosome (APC/C) regulates the assembly of the DDR complex by targeting DMS3 for degradation. We found that a substantial set of RdDM loci was commonly de-repressed in apc/c and pol v mutants, and that the defects in RdDM activity resulted from up-regulated DMS3 protein levels, which finally caused reduced Pol V recruitment. DMS3 was ubiquitinated by APC/C for degradation in a D box-dependent manner. Competitive binding assays and gel filtration analyses showed that a proper level of DMS3 is critical for the assembly of the DDR complex. Consistent with the importance of the level of DMS3, overaccumulation of DMS3 caused defective RdDM activity, phenocopying the apc/c and dms3 mutants. Moreover, DMS3 is expressed in a cell cycle-dependent manner. Collectively, these findings provide direct evidence as to how the assembly of the DDR complex is regulated and uncover a safeguarding role of APC/C in the regulation of RdDM activity.

Transcriptional dynamics of transposable elements when converting fibroblast cells of Macaca mulatta to neuroepithelial stem cells.

BACKGROUND: Transposable elements (TE) account for more than 50% of human genome. It has been reported that some types of TEs are dynamically regulated in the reprogramming of human cell lines. However, it is largely unknown whether some TEs in Macaca mulatta are also regulated during the reprogramming of cell lines of monkey. RESULTS: Here, we systematically examined the transcriptional activities of TEs during the conversion of Macaca mulatta fibroblast cells to neuroepithelial stem cells (NESCs). Hundreds of TEs were dynamically regulated during the reprogramming of Macaca mulatta fibroblast cells. Furthermore, 48 Long Terminal Repeats (LTRs), as well as some integrase elements, of Macaca endogenous retrovirus 3 (MacERV3) were transiently activated during the early stages of the conversion process, some of which were further confirmed with PCR experiments. These LTRs were potentially bound by critical transcription factors for reprogramming, such as KLF4 and ETV5. CONCLUSION: These results suggest that the transcription of TEs are delicately regulated during the reprogramming of Macaca mulatta fibroblast cells. Although the family of ERVs activated during the reprogramming of fibroblast cells in Macaca mulatta is different from those in the reprogramming of human fibroblast cells, our results suggest that the activation of some ERVs is a conserved mechanism in primates for converting fibroblast cells to stem cells.

Genome-wide identification and comprehensive analysis of microRNAs and phased small interfering RNAs in watermelon.

BACKGROUND: MicroRNAs (miRNAs) are a class of endogenous small non-coding RNAs involved in the post-transcriptional gene regulation and play a critical role in plant growth, development and stress responses. Watermelon (Citrullus lanatus L.) is one of the important agricultural crops worldwide. However, the watermelon miRNAs and phasiRNAs and their functions are not well explored. RESULTS: Here we carried out computational and experimental analysis of miRNAs and phased small interfering RNAs (phasiRNAs) in watermelon by analyzing 14 small RNA profiles from roots, leaves, androecium, petals, and fruits, and one published small RNA profile of mixed tissues. To identify the targets of miRNAs and phasiRNAs, we generated a degradome profile for watermelon leaf which is analyzed using the SeqTar algorithm. We identified 97 conserved pre-miRNAs, of which 58 have not been reported previously and 348 conserved mature miRNAs without precursors. We also found 9 novel pre-miRNAs encoding 18 mature miRNAs. One hundred and one 21 nucleotide (nt) PHAS loci, and two hundred and forty one 24 nt PHAS loci were also identified. We identified 127 conserved targets of the conserved miRNAs and TAS3-derived tasiRNAs by analyzing a degradome profile of watermelon leaf. CONCLUSIONS: The presented results provide a comprehensive view of small regulatory RNAs and their targets in watermelon.

Phased secondary small interfering RNAs in Panaxnotoginseng.

BACKGROUND: Recent results demonstrated that either non-coding or coding genes generate phased secondary small interfering RNAs (phasiRNAs) guided by specific miRNAs. Till now, there is no studies for phasiRNAs in Panax notoginseng (Burk.) F.H. Chen (P. notoginseng), an important traditional Chinese herbal medicinal plant species. METHODS: Here we performed a genome-wide discovery of phasiRNAs and its host PHAS loci in P. notoginseng by analyzing small RNA sequencing profiles. Degradome sequencing profile was used to identify the trigger miRNAs of these phasiRNAs and potential targets of phasiRNAs. We also used RLM 5'-RACE to validate some of the identified phasiRNA targets. RESULTS: After analyzing 24 small RNA sequencing profiles of P. notoginseng, 204 and 90 PHAS loci that encoded 21 and 24 nucleotide (nt) phasiRNAs, respectively, were identified. Furthermore, we found that phasiRNAs produced from some pentatricopeptide repeat-contain (PPR) genes target another layer of PPR genes as validated by both the degradome sequencing profile and RLM 5'-RACE analysis. We also found that miR171 with 21 nt triggers the generations of 21 nt phasiRNAs from its conserved targets. CONCLUSIONS: We validated that some phasiRNAs generated from PPRs and TASL genes are functional by targeting other PPRs in trans. These results provide the first set of PHAS loci and phasiRNAs in P. notoginseng, and enhance our understanding of PHAS in plants.

Accurate detection for a wide range of mutation and editing sites of microRNAs from small RNA high-throughput sequencing profiles.

Various types of mutation and editing (M/E) events in microRNAs (miRNAs) can change the stabilities of pre-miRNAs and/or complementarities between miRNAs and their targets. Small RNA (sRNA) high-throughput sequencing (HTS) profiles can contain many mutated and edited miRNAs. Systematic detection of miRNA mutation and editing sites from the huge volume of sRNA HTS profiles is computationally difficult, as high sensitivity and low false positive rate (FPR) are both required. We propose a novel method (named MiRME) for an accurate and fast detection of miRNA M/E sites using a progressive sequence alignment approach which refines sensitivity and improves FPR step-by-step. From 70 sRNA HTS profiles with over 1.3 billion reads, MiRME has detected thousands of statistically significant M/E sites, including 3'-editing sites, 57 A-to-I editing sites (of which 32 are novel), as well as some putative non-canonical editing sites. We demonstrated that a few non-canonical editing sites were not resulted from mutations in genome by integrating the analysis of genome HTS profiles of two human cell lines, suggesting the existence of new editing types to further diversify the functions of miRNAs. Compared with six existing studies or methods, MiRME has shown much superior performance for the identification and visualization of the M/E sites of miRNAs from the ever-increasing sRNA HTS profiles.

A repertoire of intronic lariat RNAs reveals tissue-specific regulation and target mimicry potential in plants.

Lariat RNA is concomitantly produced by excised intron during RNA splicing, which is usually debranched by DBR1, an RNA debranching enzyme. However, increasing evidence showed that some lariat RNA could escape debranching. Little is known about how and why these lariat RNAs could be retained. By comparing the atlas of lariat RNAs between the non-dividing cell (mature pollen) and three actively dividing tissues (young shoot apex, young seeds, and young roots), we identified hundreds to thousands of lariat RNA naturally retained in each tissue, and the incidence of lariat RNA retention is much less in shoot apex while much more in pollen. Many lariat RNAs derived from the same intron or different lariat RNAs from the same pre-mRNA could be retained in one tissue while degraded in the other tissues. By deciphering lariat RNA sequences, we identified an AG-rich (RAAAAVAAAR) motif and a UC-rich (UCUCUYUCUC) motif for pollen-specific and the other three tissues-retained lariat RNAs, respectively. Reconstitution of the pollen-specific AG-rich motif indeed enhanced lariat RNA retention in plants. Biologically, hundreds of lariat RNAs harbored miRNA binding sites, and dual-luciferase reporter assay showed that these natural lariat RNAs had the potential to protect expression of miRNA target genes. Collectively, our results uncover that selective retention of lariat RNA is an actively regulatory process, and provide new insights into understanding how lariat RNA metabolism may impact miRNA activity.

Identification of Intronic Lariat-Derived Circular RNAs in Arabidopsis by RNA Deep Sequencing.

Lariat RNAs are well-known by-products of pre-mRNA splicing in eukaryotes, which are produced by the excised introns when the 5' splice site (5' ss) joins with the branchpoint (BP) during splicing. In general, most of lariat RNAs are usually linearized by RNA debranching enzyme 1 (DBR1), followed by degradation for intron turnover. However, with the high-throughput RNA sequencing technology and bioinformatics methods, increasing evidences have shown that many lariat RNAs can stably accumulate under physiological conditions in both animals and plants. Here, we describe a large-scale analysis to systematically identify the lariat RNAs (i.e., intronic circular RNAs) in Arabidopsis by utilizing the RNA-sequencing data.

The chromatin remodeling complex imitation of switch controls stamen filament elongation by promoting jasmonic acid biosynthesis in Arabidopsis.

Plant reproduction requires the coordinated development of both male and female reproductive organs. Jasmonic acid (JA) plays an essential role in stamen filament elongation. However, the mechanism by which the JA biosynthesis genes are regulated to promote stamen elongation remains unclear. Here, we show that the chromatin remodeling complex Imitation of Switch (ISWI) promotes stamen filament elongation by regulating JA biosynthesis. We show that AT-Rich Interacting Domain 5 (ARID5) interacts with CHR11, CHR17, and RLT1, several known subunits of ISWI. Mutations in ARID5 and RLTs caused a reduced seed set due to greatly shortened stamen filaments. RNA-seq analyses reveal that the expression of key genes responsible for JA biosynthesis is significantly down-regulated in the arid5 and rlt mutants. Consistently, the JA levels are drastically decreased in both arid5 and rlt mutants. Chromatin immunoprecipitation-quantitative PCR analyses further show that ARID5 is recruited to the chromatin of JA biosynthesis genes. Importantly, exogenous JA treatments can fully rescue the defects of stamen filament elongation in both arid5 and rlt mutants, leading to the partial recovery of fertility. Our results provide a clue how JA biosynthesisis positively regulated by the chromatin remodeling complex ISWI, thereby promoting stamen filament elongation in Arabidopsis.

Identifying microRNAs and Their Editing Sites in Macaca mulatta.

MicroRNAs (miRNAs) are small non-coding RNAs that are critical in post-transcriptional regulation. Macaca mulatta is an important nonhuman primate that is often used in basic and translational researches. However, the annotation of miRNAs in Macaca mulatta is far from complete, and there are no reports of miRNA editing events in Macaca mulatta, although editing may affect the biogenesis or functions of the miRNAs. To improve miRNA annotation and to reveal editing events of miRNAs in Macaca mulatta, we generated 12 small RNA profiles from eight tissues and performed comprehensive analysis of these profiles. We identified 479 conserved pre-miRNAs that have not been reported in Macaca mulatta and 17 species specific miRNAs. Furthermore, we identified 3386 editing sites with significant editing levels from 471 pre-miRNAs after analyzing the 12 self-generated and 58 additional published sRNA-seq profiles from 17 different types of organs or tissues. In addition to 16 conserved A-to-I editing sites, we identified five conserved C-to-U editing sites in miRNAs of Macaca mulatta and Homo sapiens. We also identified 11 SNPs in the miRNAs of Macaca mulatta. The analysis of the potential targets of 69 miRNAs with editing or mutation events in their seed regions suggest that these editing or mutation events severely changed their targets and their potential functions. These results significantly increase our understanding of miRNAs and their mutation/editing events in Macaca mulatta.

A lariat-derived circular RNA is required for plant development in Arabidopsis.

Lariat RNA is produced during pre-mRNA splicing, and it is traditionally thought as by-products, due to the quick turnover by debranching followed by degradation. However, recent findings identified many lariat RNAs accumulate with a circular form in higher eukaryotes. Although the remarkable accumulation, biological consequence of lariat-derived circular RNAs (here we name laciRNAs) remains largely unknown. Here, we report that a specific laciRNA from At5g37720 plays an essential role in plant development by regulating gene expression globally. We focus on 17 laciRNAs with accumulation in wild type plants by circular RNA sequencing in Arabidopsis. To determine biological functions of these laciRNAs, we constructed one pair of transgenic plants for each laciRNA, in which the local gene with or without introns was over-expressed in wild type plants, respectively. By comparing morphological phenotypes and transcriptomic profiles between two classes of transgenic plants, we show that over-expression of the laciRNA derived from the 1st intron of At5g37720 causes pleiotropic phenotypes, including curly and clustered leaf, late flowering, reduced fertility, and accompanied with altered expression of approximately 800 genes. Our results provide another example that a specific plant circular RNA regulates gene expression in a similar manner to that of other non-coding RNAs under physiological conditions.

Small RNA profiles from Panax notoginseng roots differing in sizes reveal correlation between miR156 abundances and root biomass levels.

Plant genomes encode several classes of small regulatory RNAs (sRNAs) that play critical roles in both development and stress responses. Panax notoginseng (Burk.) F.H. Chen (P. notoginseng) is an important traditional Chinese herbal medicinal plant species for its haemostatic effects. Therefore, the root yield of P. notoginseng is a major economically important trait since the roots of P. notoginseng are the parts used to produce medicine. To identify sRNAs that are critical for the root biomass of P. notoginseng, we performed a comprehensive study of miRNA transcriptomes from P. notoginseng roots of different biomasses. We identified 675 conserved miRNAs, of which 180 pre-miRNAs are also identified, and three TAS3 loci in P. notoginseng. By using degradome sequencing, we identified 79 conserved miRNA:target or tasiRNA:target interactions, of which eight were further confirmed with the RLM 5'-RACE experiments. More importantly, our results revealed that a member of miR156 family and one of its SPL target genes have inverse expression levels, which is tightly correlated with greater root biomass contents. These results not only contributes to overall understanding of post-transcriptional gene regulation in roots of P. notoginseng but also could serve as markers for breeding P. notoginseng with greater root yield.

Small RNA profiles in soybean primary root tips under water deficit.

BACKGROUND: Soybean (Glycine max) production is significantly hampered by frequent droughts in many regions of the world including the United States. Identifying microRNA (miRNA)-controlled posttranscriptional gene regulation under drought will enhance our understanding of molecular basis of drought tolerance in this important cash crop. Indeed, miRNA profiles in soybean exposed to drought were studied but not from the primary root tips, which is not only a main zone of water uptake but also critical for water stress sensing and signaling. METHODS: Here we report miRNA profiles specifically from well-watered and water-stressed primary root tips (0 to 8 mm from the root apex) of soybean. Small RNA sequencing confirmed the expression of vastly diverse miRNA (303 individual miRNAs) population, and, importantly several conserved miRNAs were abundantly expressed in primary root tips. RESULTS: Notably, 12 highly conserved miRNA families were differentially regulated in response to water-deficit; six were upregulated while six others were downregulated at least by one fold (log2) change. Differentially regulated soybean miRNAs are targeting genes include auxin response factors, Cu/Zn Superoxide dismutases, laccases and plantacyanin and several others. CONCLUSIONS: These results highlighted the importance of miRNAs in primary root tips both under control and water-deficit conditions; under control conditions, miRNAs could be important for cell division, cell elongation and maintenance of the root apical meristem activity including quiescent centre whereas under water stress differentially regulated miRNAs could decrease auxin signaling and oxidative stress as well as other metabolic processes that save energy and water.