Genome-wide analysis of single non-templated nucleotides in plant endogenous siRNAs and miRNAs

Genome-wide identification of nitrate-responsive microRNAs by small RNA sequencing in the rice restorer cultivar Nanhui 511

Rice productivity relies heavily on nitrogen fertilization, and improving nitrogen use efficiency (NUE) is important for hybrid rice breeding. Reducing nitrogen inputs is the key to achieving sustainable rice production and reducing environmental problems. Here, we analyzed the genome-wide transcriptomic changes in microRNAs (miRNAs) in the indica rice restorer cultivar Nanhui 511 (NH511) under high (HN) and low nitrogen (LN) conditions. The results showed that NH511 is sensitive to nitrogen supplies and HN conditions promoted the growth its lateral roots at the seedling stage. Furthermore, we identified 483 known miRNAs and 128 novel miRNAs by small RNA sequencing in response to nitrogen in NH511. We also detected 100 differentially expressed genes (DEGs), including 75 upregulated and 25 downregulated DEGs, under HN conditions. Among these DEGs, 43 miRNAs that exhibited a 2-fold change in their expression were identified in response to HN conditions, including 28 upregulated and 15 downregulated genes. Additionally, some differentially expressed miRNAs were further validated by qPCR analysis, which showed that miR443, miR1861b, and miR166k-3p were upregulated, whereas miR395v and miR444b.1 were downregulated under HN conditions. Moreover, the degradomes of possible target genes for miR166k-3p and miR444b.1 and expression variations were analyzed by qPCR at different time points under HN conditions. Our findings revealed comprehensive expression profiles of miRNAs responsive to HN treatments in an indica rice restorer cultivar, which advances our understanding of the regulation of nitrogen signaling mediated by miRNAs and provides novel data for high-NUE hybrid rice cultivation.

TENT2, TUT4, and TUT7 selectively regulate miRNA sequence and abundance

TENTs generate miRNA isoforms by 3' tailing. However, little is known about how tailing regulates miRNA function. Here, we generate isogenic HEK293T cell lines in which TENT2, TUT4 and TUT7 are knocked out individually or in combination. Together with rescue experiments, we characterize TENT-specific effects by deep sequencing, Northern blot and in vitro assays. We find that 3' tailing is not random but highly specific. In addition to its known adenylation, TENT2 contributes to guanylation and uridylation on mature miRNAs. TUT4 uridylates most miRNAs whereas TUT7 is dispensable. Removing adenylation has a marginal impact on miRNA levels. By contrast, abolishing uridylation leads to dysregulation of a set of miRNAs. Besides let-7, miR-181b and miR-222 are negatively regulated by TUT4/7 via distinct mechanisms while the miR-888 cluster is upregulated specifically by TUT7. Our results uncover the selective actions of TENTs in generating 3' isomiRs and pave the way to investigate their functions.

A DCL3 dicing code within Pol IV-RDR2 transcripts diversifies the siRNA pool guiding RNA-directed DNA methylation

In plants, selfish genetic elements, including retrotransposons and DNA viruses, are transcriptionally silenced by RNA-directed DNA methylation. Guiding the process are short interfering RNAs (siRNAs) cut by DICER-LIKE 3 (DCL3) from double-stranded precursors of ~30 bp that are synthesized by NUCLEAR RNA POLYMERASE IV (Pol IV) and RNA-DEPENDENT RNA POLYMERASE 2 (RDR2). We show that Pol IV's choice of initiating nucleotide, RDR2's initiation 1-2 nt internal to Pol IV transcript ends and RDR2's terminal transferase activity collectively yield a code that influences which precursor end is diced and whether 24 or 23 nt siRNAs are produced. By diversifying the size, sequence, and strand specificity of siRNAs derived from a given precursor, alternative patterns of DCL3 dicing allow for maximal siRNA coverage at methylated target loci.

An atypical class of non-coding small RNAs is produced in rice leaves upon bacterial infection

Non-coding small RNAs (sRNA) act as mediators of gene silencing and regulate plant growth, development and stress responses. Early insights into plant sRNAs established a role in antiviral defense and they are now extensively studied across plant-microbe interactions. Here, sRNA sequencing discovered a class of sRNA in rice (Oryza sativa) specifically associated with foliar diseases caused by Xanthomonas oryzae bacteria. Xanthomonas-induced small RNAs (xisRNAs) loci were distinctively upregulated in response to diverse virulent strains at an early stage of infection producing a single duplex of 20-22 nt sRNAs. xisRNAs production was dependent on the Type III secretion system, a major bacterial virulence factor for host colonization. xisRNA loci overlap with annotated transcripts sequences, with about half of them encoding protein kinase domain proteins. A number of the corresponding rice cis-genes have documented functions in immune signaling and xisRNA loci predominantly coincide with the coding sequence of a conserved kinase motif. xisRNAs exhibit features of small interfering RNAs and their biosynthesis depend on canonical components OsDCL1 and OsHEN1. xisRNA induction possibly mediates post-transcriptional gene silencing but they do not broadly suppress cis-genes expression on the basis of mRNA-seq data. Overall, our results identify a group of unusual sRNAs with a potential role in plant-microbe interactions.

A tale of non-canonical tails: gene regulation by post-transcriptional RNA tailing

RNA tailing, or the addition of non-templated nucleotides to the 3' end of RNA, is the most frequent and conserved type of RNA modification. The addition of tails and their composition reflect RNA maturation stages and have important roles in determining the fate of the modified RNAs. Apart from canonical poly(A) polymerases, which add poly(A) tails to mRNAs in a transcription-coupled manner, a family of terminal nucleotidyltransferases (TENTs), including terminal uridylyltransferases (TUTs), modify RNAs post-transcriptionally to control RNA stability and activity. The human genome encodes 11 different TENTs with distinct substrate specificity, intracellular localization and tissue distribution. In this Review, we discuss recent advances in our understanding of non-canonical RNA tails, with a focus on the functions of human TENTs, which include uridylation, mixed tailing and post-transcriptional polyadenylation of mRNAs, microRNAs and other types of non-coding RNA.

Reconstitution of siRNA Biogenesis In Vitro: Novel Reaction Mechanisms and RNA Channeling in the RNA-Directed DNA Methylation Pathway

Eukaryotes deploy RNA-mediated gene silencing pathways to guard their genomes against selfish genetic elements, such as transposable elements and invading viruses. In plants, RNA-directed DNA methylation (RdDM) is used to silence selfish elements at the level of transcription. This process involves 24-nt short interfering RNAs (siRNAs) and longer noncoding RNAs to which the siRNAs base-pair. Recently, we showed that 24-nt siRNA biogenesis could be recapitulated in the test tube using purified enzymes, yielding biochemical answers to numerous questions left unresolved by prior genetic and genomic studies. Interestingly, each enzyme has activities that program what happens in the next step, thus channeling the RNAs within the RdDM pathway and restricting their diversion into alternative pathways. However, a similar mechanistic understanding is lacking for other important steps of the RdDM pathway. We discuss some of the steps most in need of biochemical investigation and important questions still in need of answers.

Reaction Mechanisms of Pol IV, RDR2, and DCL3 Drive RNA Channeling in the siRNA-Directed DNA Methylation Pathway

In eukaryotes with multiple small RNA pathways, the mechanisms that channel RNAs within specific pathways are unclear. Here, we reveal the reactions that account for channeling in the small interfering RNA (siRNA) biogenesis phase of the Arabidopsis RNA-directed DNA methylation pathway. The process begins with template DNA transcription by NUCLEAR RNA POLYMERASE IV (Pol IV), whose atypical termination mechanism, induced by nontemplate DNA base-pairing, channels transcripts to the associated RNA-dependent RNA polymerase RDR2. RDR2 converts Pol IV transcripts into double-stranded RNAs and then typically adds an extra untemplated 3' terminal nucleotide to the second strands. The dicer endonuclease DCL3 cuts resulting duplexes to generate 24- and 23-nt siRNAs. The 23-nt RNAs bear the untemplated terminal nucleotide of the RDR2 strand and are underrepresented among ARGONAUTE4-associated siRNAs. Collectively, our results provide mechanistic insights into Pol IV termination, Pol IV-RDR2 coupling, and RNA channeling, from template DNA transcription to siRNA strand discrimination.

Bioinformatic Analysis of Small RNA Sequencing Libraries

Bioinformatic analysis of small RNA sequencing libraries consists of transforming a series of small RNA sequencing experiment fastq files into a table containing small RNA sequences and their abundance. This is achieved by cleaning the reads, aligning the cleaned reads to a reference, and parsing the alignment results. In this protocol we present the most common option, and the rationale, for each of these steps.

RNA uridylation and decay in plants

RNA uridylation consists of the untemplated addition of uridines at the 3' extremity of an RNA molecule. RNA uridylation is catalysed by terminal uridylyltransferases (TUTases), which form a subgroup of the terminal nucleotidyltransferase family, to which poly(A) polymerases also belong. The key role of RNA uridylation is to regulate RNA degradation in a variety of eukaryotes, including fission yeast, plants and animals. In plants, RNA uridylation has been mostly studied in two model species, the green algae Chlamydomonas reinhardtii and the flowering plant Arabidopsis thaliana Plant TUTases target a variety of RNA substrates, differing in size and function. These RNA substrates include microRNAs (miRNAs), small interfering silencing RNAs (siRNAs), ribosomal RNAs (rRNAs), messenger RNAs (mRNAs) and mRNA fragments generated during post-transcriptional gene silencing. Viral RNAs can also get uridylated during plant infection. We describe here the evolutionary history of plant TUTases and we summarize the diverse molecular functions of uridylation during RNA degradation processes in plants. We also outline key points of future research.This article is part of the theme issue '5' and 3' modifications controlling RNA degradation'.

Size Distribution of Small Interfering RNAs in Various Organs at Different Developmental Stages is Primarily Determined by the Dicing Activity of Dicer-Like Proteins in Plants

RNA silencing is a fundamental mechanism to maintain plant growth and development, and regulation of the size distribution of small interfering RNAs (siRNAs) is critical in the control of normal gene expression throughout a plant's life cycle. However, the cause of organ- and developmental stage-specific accumulation of siRNAs has never been reported. Whereas 24 nt siRNAs accumulated about 5.3-fold more than 21 nt siRNAs in Arabidopsis rosette leaves, 21 and 24 nt siRNAs accumulated to similar levels in Arabidopsis pollen grains, rice spikelets and maize anthers. We successfully detected two distinct double-stranded RNA (dsRNA)-cleaving activities that produced 21 and 24 nt RNAs in cell-free extracts prepared from various organs at different developmental stages of A. thaliana, Brassica rapa, rice and maize. Although DCL4 transcript was expressed more than DCL3 transcript in most organs, the 21 nt RNA-producing activity of DCL4 or its orthologs was very low and was 5- to 10-fold lower than the 24 nt RNA-producing activity of DCL3 or its orthologs particularly in leaves, indicating that DCL4 activity is negatively regulated translationally or post-translationally in leaves. High dicing activity of DCL3 and DCL4 was detected in immature inflorescences, developing seeds, germinating embryos and callus, all of which contain actively dividing cells. In various organs at different developmental stages, the size distribution of siRNAs was positively correlated with the dicing activity of two Dicers, DCL3 and DCL4, or their orthologs. Taken together, the size distribution of siRNAs in most organs is primarily determined by the dicing activity of DCL3 and DCL4.

Cis-directed cleavage and nonstoichiometric abundances of 21-nucleotide reproductive phased small interfering RNAs in grasses

Post-transcriptional gene silencing in plants results from independent activities of diverse small RNA types. In anthers of grasses, hundreds of loci yield noncoding RNAs that are processed into 21- and 24-nucleotide (nt) phased small interfering RNAs (phasiRNAs); these are triggered by miR2118 and miR2275. We characterized these 'reproductive phasiRNAs' from rice (Oryza sativa) panicles and anthers across seven developmental stages. Our computational analysis identified characteristics of the 21-nt reproductive phasiRNAs that impact their biogenesis, stability, and potential functions. We demonstrate that 21-nt reproductive phasiRNAs can function in cis to target their own precursors. We observed evidence of this cis regulatory activity in both rice and maize (Zea mays). We validated this activity with evidence of cleavage and a resulting shift in the pattern of phasiRNA production. We characterize biases in phasiRNA biogenesis, demonstrating that the Pol II-derived 'top' strand phasiRNAs are consistently higher in abundance than the bottom strand. The first phasiRNA from each precursor overlaps the miR2118 target site, and this impacts phasiRNA accumulation or stability, evident in the weak accumulation of this phasiRNA position. Additional influences on this first phasiRNA duplex include the sequence composition and length, and we show that these factors impact Argonaute loading.

Cis-directed cleavage and nonstoichiometric abundances of 21-nucleotide reproductive phased small interfering RNAs in grasses

Post-transcriptional gene silencing in plants results from independent activities of diverse small RNA types. In anthers of grasses, hundreds of loci yield noncoding RNAs that are processed into 21- and 24-nucleotide (nt) phased small interfering RNAs (phasiRNAs); these are triggered by miR2118 and miR2275. We characterized these 'reproductive phasiRNAs' from rice (Oryza sativa) panicles and anthers across seven developmental stages. Our computational analysis identified characteristics of the 21-nt reproductive phasiRNAs that impact their biogenesis, stability, and potential functions. We demonstrate that 21-nt reproductive phasiRNAs can function in cis to target their own precursors. We observed evidence of this cis regulatory activity in both rice and maize (Zea mays). We validated this activity with evidence of cleavage and a resulting shift in the pattern of phasiRNA production. We characterize biases in phasiRNA biogenesis, demonstrating that the Pol II-derived 'top' strand phasiRNAs are consistently higher in abundance than the bottom strand. The first phasiRNA from each precursor overlaps the miR2118 target site, and this impacts phasiRNA accumulation or stability, evident in the weak accumulation of this phasiRNA position. Additional influences on this first phasiRNA duplex include the sequence composition and length, and we show that these factors impact Argonaute loading.

Evolution of miRNA Tailing by 3' Terminal Uridylyl Transferases in Metazoa

In bilaterian animals the 3′ ends of microRNAs (miRNAs) are frequently modified by tailing and trimming. These modifications affect miRNA-mediated gene regulation by modulating miRNA stability. Here, we analyzed data from three nonbilaterian animals: two cnidarians (Nematostella vectensis and Hydra magnipapillata) and one poriferan (Amphimedon queenslandica). Our analysis revealed that nonbilaterian miRNAs frequently undergo modifications like the bilaterian counterparts: the majority are expressed as different length isoforms and frequent modifications of the 3′ end by mono U or mono A tailing are observed. Moreover, as the factors regulating miRNA modifications are largely uncharacterized in nonbilaterian animal phyla, in present study, we investigated the evolution of 3′ terminal uridylyl transferases (TUTases) that are known to involved in miRNA 3′ nontemplated modifications in Bilateria. Phylogenetic analysis on TUTases showed that TUTase1 and TUTase6 are a result of duplication in bilaterians and that TUTase7 and TUTase4 are the result of a vertebrate-specific duplication. We also find an unexpected number of Drosophila-specific gene duplications and domain losses in most of the investigated gene families. Overall, our findings shed new light on the evolutionary history of TUTases in Metazoa, as they reveal that this core set of enzymes already existed in the last common ancestor of all animals and was probably involved in modifying small RNAs in a similar fashion to its present activity in bilaterians.

Click Chemical Ligation-Initiated On-Bead DNA Polymerization for the Sensitive Flow Cytometric Detection of 3'-Terminal 2'-O-Methylated Plant MicroRNA

A versatile flow cytometric strategy is developed for the sensitive detection of plant microRNA (miRNA) by coupling the target-templated click nucleic acid ligation (CNAL) with on-bead terminal enzymatic DNA polymerization (TEP). Unlike ligase-catalyzed ligation reaction, the plant miRNA-templated enzyme-free CNAL between two single-stranded DNA (ssDNA) probes, respectively modified with Aza-dibenzocyclooctyne (Aza-DBCO) and N₃, can not only simplify the operation, but also achieve a much higher ligation efficiency. More importantly, the undesirable nonspecific ligation between the Aza-DBCO- and N₃-modified ssDNA, can be effectively eliminated by adding Tween-20, which allows the use of cycling CNAL (CCNAL) in a background-free manner. So each plant miRNA can template many rounds of CNAL reaction to produce numerous ligation products, forming efficient signal amplification. The ligated ssDNA can be anchored on the magnetic beads (MBs) with the 3'-OH termini exposed outside. Then terminal deoxynucleotidyl transferase (TdT), a sequence-independent and template-free polymerase, would specifically catalyze the DNA polymerization along these 3'-OH termini on the MBs, forming poly(T) tails up to thousands of nucleotides long. Each poly(T) tail allows specific binding of numerous 6-carboxyfluorescein (FAM)-labeled poly(A)25 oligonucleotides to accumulate a lot of fluorophores on the MBs, leading to the second step of signal amplification. By integrating the advantages of CCNAL-TEP for highly efficient signal amplification and robust MBs signal readout with powerful flow cytometer, high sensitivity is achieved and the detection limit of plant miRNA has been pushed down to a low level of 5 fM with high specificity to well discriminate even single-base difference between miRNA targets.

Characterization of Conserved and Novel microRNAs in Lilium lancifolium Thunb. by High-Throughput Sequencing

MicroRNAs (miRNAs) are among the class of noncoding small RNA molecules and play a crucial role in post-transcriptional regulation in plants. Although Lilium is one of the most popular ornamental flowers worldwide, however, there is no report on miRNAs identification. In the present study, therefore, miRNAs and their targets were identified from flower, leaf, bulblet and bulb of Lilium lancifolium Thunb. by high-throughput sequencing and bioinformatics analysis. In this study, a total of 38 conserved miRNAs belonging to 17 miRNA families and 44 novel miRNAs were identified. In total, 366 target genes for conserved miRNAs and 415 target genes for novel miRNAs were predicted. The majority of the target genes for conserved miRNAs were transcriptional factors and novel miRNAs targeted mainly protein coding genes. A total of 53 cleavage sites belonging to 6 conserved miRNAs families and 14 novel miRNAs were identified using degradome sequencing. Twenty-three miRNAs were randomly selected, then, their credibility was confirmed using northern blot or stem-loop qRT-PCR. The results from qRT-PCR analysis showed the expression pattern of 4 LL-miRNAs was opposite to their targets. Therefore, our finding provides an important basis to understand the biological functions of miRNAs in Lilium.

AGO4 is specifically required for heterochromatic siRNA accumulation at Pol V-dependent loci in Arabidopsis thaliana

In plants, 24 nucleotide long heterochromatic siRNAs (het-siRNAs) transcriptionally regulate gene expression by RNA-directed DNA methylation (RdDM). The biogenesis of most het-siRNAs depends on the plant-specific RNA polymerase IV (Pol IV), and ARGONAUTE4 (AGO4) is a major het-siRNA effector protein. Through genome-wide analysis of sRNA-seq data sets, we found that AGO4 is required for the accumulation of a small subset of het-siRNAs. The accumulation of AGO4-dependent het-siRNAs also requires several factors known to participate in the effector portion of the RdDM pathway, including RNA POLYMERASE V (POL V), DOMAINS REARRANGED METHYLTRANSFERASE 2 (DRM2) and SAWADEE HOMEODOMAIN HOMOLOGUE 1 (SHH1). Like many AGO proteins, AGO4 is an endonuclease that can 'slice' RNAs. We found that a slicing-defective AGO4 was unable to fully recover AGO4-dependent het-siRNA accumulation from ago4 mutant plants. Collectively, our data suggest that AGO4-dependent siRNAs are secondary siRNAs dependent on the prior activity of the RdDM pathway at certain loci.