Jixian Zhai

Massive Analysis of Rice Small RNAs: Mechanistic Implications of Regulated MicroRNAs and Variants for Differential Target RNA Cleavage

This massive analysis of rice small RNAs evaluated annotated microRNAs (miRNAs) and identified new miRNAs and miRNAs regulated by environmental stresses. Of particular interest are miRNA family members with distinct sequences and organ-preferential expression; some of these guide differential target cleavage and provide new insight about how an agriculturally significant phenotype may be regulated. Small RNAs have a variety of important roles in plant development, stress responses, and other processes. They exert their influence by guiding mRNA cleavage, translational repression, and chromatin modification. To identify previously unknown rice (Oryza sativa) microRNAs (miRNAs) and those regulated by environmental stress, 62 small RNA libraries were constructed from rice plants and used for deep sequencing with Illumina technology. The libraries represent several tissues from control plants and plants subjected to different environmental stress treatments. More than 94 million genome-matched reads were obtained, resulting in more than 16 million distinct small RNA sequences. This allowed an evaluation of ~400 annotated miRNAs with current criteria and the finding that among these, ~150 had small interfering RNA–like characteristics. Seventy-six new miRNAs were found, and miRNAs regulated in response to water stress, nutrient stress, or temperature stress were identified. Among the new examples of miRNA regulation were members of the same miRNA family that were differentially regulated in different organs and had distinct sequences Some of these distinct family members result in differential target cleavage and provide new insight about how an agriculturally important rice phenotype could be regulated in the panicle. This high-resolution analysis of rice miRNAs should be relevant to plant miRNAs in general, particularly in the Poaceae.

miRNAs trigger widespread epigenetically activated siRNAs from transposons in Arabidopsis

In plants, post-transcriptional gene silencing (PTGS) is mediated by DICER-LIKE 1 (DCL1)-dependent microRNAs (miRNAs), which also trigger 21-nucleotide secondary short interfering RNAs (siRNAs) via RNA-DEPENDENT RNA POLYMERASE 6 (RDR6), DCL4 and ARGONAUTE 1 (AGO1), whereas transcriptional gene silencing (TGS) of transposons is mediated by 24-nucleotide heterochromatic (het)siRNAs, RDR2, DCL3 and AGO4 (ref. 4). Transposons can also give rise to abundant 21-nucleotide ‘epigenetically activated’ small interfering RNAs (easiRNAs) in DECREASED DNA METHYLATION 1 (ddm1) and DNA METHYLTRANSFERASE 1 (met1) mutants, as well as in the vegetative nucleus of pollen grains and in dedifferentiated plant cell cultures. Here we show that easiRNAs in Arabidopsis thaliana resemble secondary siRNAs, in that thousands of transposon transcripts are specifically targeted by more than 50 miRNAs for cleavage and processing by RDR6. Loss of RDR6, DCL4 or DCL1 in a ddm1 background results in loss of 21-nucleotide easiRNAs and severe infertility, but 24-nucleotide hetsiRNAs are partially restored, supporting an antagonistic relationship between PTGS and TGS. Thus miRNA-directed easiRNA biogenesis is a latent mechanism that specifically targets transposon transcripts, but only when they are epigenetically reactivated during reprogramming of the germ line. This ancient recognition mechanism may have been retained both by transposons to evade long-term heterochromatic silencing and by their hosts for genome defence.

SDG714, a histone H3K9 methyltransferase, is involved in Tos17 DNA methylation and transposition in rice.

Although the role of H3K9 methylation in rice (Oryza sativa) is unclear, in Arabidopsis thaliana the loss of histone H3K9 methylation by mutation of Kryptonite [also known as SU(VAR)3-9 homolog] reduces genome-wide DNA methylation and increases the transcription of transposable elements. Here, we report that rice SDG714 (for SET Domain Group Protein714) encodes a histone H3K9-specific methyltransferase. The C terminus of SDG714 confers enzymatic activity and substrate specificity, whereas the N terminus localizes it in the nucleus. Loss-of-function mutants of SDG714 (SDG714IR transformants) generated by RNA interference display a mostly glabrous phenotype as a result of the lack of macro trichomes in glumes, leaves, and culms compared with control plants. These mutants also show decreased levels of CpG and CNG cytosine methylation as well as H3K9 methylation at the Tos17 locus, a copia-like retrotransposon widely used for the generation of rice mutants. Most interestingly, loss of function of SDG714 can enhance transcription and cause the transposition of Tos17. Together, these results suggest that histone H3K9 methylation mediated by SDG714 is involved in DNA methylation, the transposition of transposable elements, and genome stability in rice.

Genome-wide Hi-C Analyses in Wild-Type and Mutants Reveal High-Resolution Chromatin Interactions in Arabidopsis

Chromosomes form 3D structures that are critical to the regulation of cellular and genetic processes. Here, we present a study of global chromatin interaction patterns in Arabidopsis thaliana. Our genome-wide approach confirmed interactions that were previously observed by other methods as well as uncovered long-range interactions such as those among small heterochromatic regions embedded in euchromatic arms. We also found that interactions are correlated with various epigenetic marks that are localized in active or silenced chromatin. Arabidopsis chromosomes do not contain large local interactive domains that resemble the topological domains described in animals but, instead, contain relatively small interactive regions scattered around the genome that contain H3K27me3 or H3K9me2. We generated interaction maps in mutants that are defective in specific epigenetic pathways and found altered interaction patterns that correlate with changes in the epigenome. These analyses provide further insights into molecular mechanisms of epigenetic regulation of the genome.

Short-Read Sequencing Technologies for Transcriptional Analyses

The technological advances in DNA sequencing over the past five years have changed our approaches to gene expression analysis, fundamentally altering the basic methods used and in most cases driving a shift from hybridization-based approaches to sequencing-based approaches. Quantitative, tag-based studies of gene expression were one of the earliest applications of these next-generation technologies, but the tremendous depth of sequencing facilitates de novo transcript discovery, which replaces traditional expressed sequence tag (EST) sequencing. In addition, these technologies have created new opportunities for understanding the generation, stability, and decay of RNA and the impacts of chromatin differences on gene expression. As we review the impact of these methods on plant biology, we also mention published studies from animal systems when the methods are broadly applicable. We can anticipate that the published work over the past few years is a harbinger of much broader studies that are yet to be published and are sure to further advance our understanding of plant genomes in a field changing at a dizzying pace.

Small RNA-directed epigenetic natural variation in Arabidopsis thaliana.

Progress in epigenetics has revealed mechanisms that can heritably regulate gene function independent of genetic alterations. Nevertheless, little is known about the role of epigenetics in evolution. This is due in part to scant data on epigenetic variation among natural populations. In plants, small interfering RNA (siRNA) is involved in both the initiation and maintenance of gene silencing by directing DNA methylation and/or histone methylation. Here, we report that, in the model plant Arabidopsis thaliana, a cluster of ∼24 nt siRNAs found at high levels in the ecotype Landsberg erecta (Ler) could direct DNA methylation and heterochromatinization at a hAT element adjacent to the promoter of FLOWERING LOCUS C (FLC), a major repressor of flowering, whereas the same hAT element in ecotype Columbia (Col) with almost identical DNA sequence, generates a set of low abundance siRNAs that do not direct these activities. We have called this hAT element MPF for Methylated region near Promoter of FLC, although de novo methylation triggered by an inverted repeat transgene at this region in Col does not alter its FLC expression. DNA methylation of the Ler allele MPF is dependent on genes in known silencing pathways, and such methylation is transmissible to Col by genetic crosses, although with varying degrees of penetrance. A genome-wide comparison of Ler and Col small RNAs identified at least 68 loci matched by a significant level of ∼24 nt siRNAs present specifically in Ler but not Col, where nearly half of the loci are related to repeat or TE sequences. Methylation analysis revealed that 88% of the examined loci (37 out of 42) were specifically methylated in Ler but not Col, suggesting that small RNA can direct epigenetic differences between two closely related Arabidopsis ecotypes.

Spatiotemporally dynamic, cell-type–dependent premeiotic and meiotic phasiRNAs in maize anthers

Significance By RNA profiling of 10 stages of maize anthers plus mature pollen, we found two distinct classes of phased small-interfering RNAs (phasiRNAs): 21-nt premeiotic phasiRNAs, after germinal and somatic cell specification, and 24-nt meiotic phasiRNAs coordinately accumulated during meiosis and persist into pollen. Sequencing of RNA from five male-sterile, anther developmental mutants—ocl4, mac1, ms23, msca1, and ameiotic1—demonstrated the involvement of specific somatic layers. Premeiotic phasiRNAs require a functional epidermis, whereas meiotic phasiRNAs require a normal tapetum. Mammalian germ cells express “prepachytene” or “pachytene” PIWI-interacting RNAs (piRNAs). Whereas differences in biogenesis indicate independent origins, grass phasiRNAs and mammalian piRNAs share developmental timing, a lack of obvious targets, and an impact on male fertility, suggesting a possible evolutionary convergence. Maize anthers, the male reproductive floral organs, express two classes of phased small-interfering RNAs (phasiRNAs). PhasiRNA precursors are transcribed by RNA polymerase II and map to low-copy, intergenic regions similar to PIWI-interacting RNAs (piRNAs) in mammalian testis. From 10 sequential cohorts of staged maize anthers plus mature pollen we find that 21-nt phased siRNAs from 463 loci appear abruptly after germinal and initial somatic cell fate specification and then diminish, whereas 24-nt phasiRNAs from 176 loci coordinately accumulate during meiosis and persist as anther somatic cells mature and haploid gametophytes differentiate into pollen. Male-sterile ocl4 anthers defective in epidermal signaling lack 21-nt phasiRNAs. Male-sterile mutants with subepidermal defects—mac1 (excess meiocytes), ms23 (defective pretapetal cells), and msca1 (no normal soma or meiocytes)—lack 24-nt phasiRNAs. ameiotic1 mutants (normal soma, no meiosis) accumulate both 21-nt and 24-nt phasiRNAs, ruling out meiotic cells as a source or regulator of phasiRNA biogenesis. By in situ hybridization, miR2118 triggers of 21-nt phasiRNA biogenesis localize to epidermis; however, 21-PHAS precursors and 21-nt phasiRNAs are abundant subepidermally. The miR2275 trigger, 24-PHAS precursors, and 24-nt phasiRNAs all accumulate preferentially in tapetum and meiocytes. Therefore, each phasiRNA type exhibits independent spatiotemporal regulation with 21-nt premeiotic phasiRNAs dependent on epidermal and 24-nt meiotic phasiRNAs dependent on tapetal cell differentiation. Maize phasiRNAs and mammalian piRNAs illustrate putative convergent evolution of small RNAs in male reproduction.

Dicer-like 3 produces transposable element-associated 24-nt siRNAs that control agricultural traits in rice

Significance The functional relationship of transposons and small RNAs remains an important question in the study of gene expression and its effect on agronomic traits. Here, we use deep sequencing of small RNAs to provide the first evidence that the rice Dicer-like 3 homolog OsDCL3a produces 24-nt small interfering RNAs (siRNAs) predominantly associated with miniature inverted repeat transposable elements (MITEs). These 24-nt siRNAs target genes adjacent to MITEs and act as broadly functioning regulators of gene expression. In particular, OsDCL3a directly targets genes involved in homeostasis of the plant hormones gibberellin and brassinosteroid, thus controlling important agricultural traits. This mechanism of fine-tuning gene expression mediated by MITEs may be conserved in organisms with genomes rich in dispersed repeats or transposable elements. Transposable elements (TEs) and repetitive sequences make up over 35% of the rice (Oryza sativa) genome. The host regulates the activity of different TEs by different epigenetic mechanisms, including DNA methylation, histone H3K9 methylation, and histone H3K4 demethylation. TEs can also affect the expression of host genes. For example, miniature inverted repeat TEs (MITEs), dispersed high copy-number DNA TEs, can influence the expression of nearby genes. In plants, 24-nt small interfering RNAs (siRNAs) are mainly derived from repeats and TEs. However, the extent to which TEs, particularly MITEs associated with 24-nt siRNAs, affect gene expression remains elusive. Here, we show that the rice Dicer-like 3 homolog OsDCL3a is primarily responsible for 24-nt siRNA processing. Impairing OsDCL3a expression by RNA interference caused phenotypes affecting important agricultural traits; these phenotypes include dwarfism, larger flag leaf angle, and fewer secondary branches. We used small RNA deep sequencing to identify 535,054 24-nt siRNA clusters. Of these clusters, ∼82% were OsDCL3a-dependent and showed significant enrichment of MITEs. Reduction of OsDCL3a function reduced the 24-nt siRNAs predominantly from MITEs and elevated expression of nearby genes. OsDCL3a directly targets genes involved in gibberellin and brassinosteroid homeostasis; OsDCL3a deficiency may affect these genes, thus causing the phenotypes of dwarfism and enlarged flag leaf angle. Our work identifies OsDCL3a-dependent 24-nt siRNAs derived from MITEs as broadly functioning regulators for fine-tuning gene expression, which may reflect a conserved epigenetic mechanism in higher plants with genomes rich in dispersed repeats or TEs.

Mutations in the Type II Protein Arginine Methyltransferase AtPRMT5 Result in Pleiotropic Developmental Defects in Arabidopsis

Human PROTEIN ARGININE METHYLTRANSFERASE5 (PRMT5) encodes a type II protein arginine (Arg) methyltransferase and its homologs in animals and yeast (Saccharomyces cerevisiae and Schizosaccharomyces pombe) are known to regulate RNA processing, signal transduction, and gene expression. However, PRMT5 homologs in higher plants have not yet been reported and the biological roles of these proteins in plant development remain elusive. Here, using conventional biochemical approaches, we purified a plant histone Arg methyltransferase from cauliflower (Brassica oleracea) that was nearly identical to AtPRMT5, an Arabidopsis (Arabidopsis thaliana) homolog of human PRMT5. AtPRMT5 methylated histone H4, H2A, and myelin basic protein in vitro. Western blot using symmetric dimethyl histone H4 Arg 3-specific antibody and thin-layer chromatography analysis demonstrated that AtPRMT5 is a type II enzyme. Mutations in AtPRMT5 caused pleiotropic developmental defects, including growth retardation, dark green and curled leaves, and FlOWERING LOCUS C (FLC)-dependent delayed flowering. Therefore, the type II protein Arg methyltransferase AtPRMT5 is involved in promotion of vegetative growth and FLC-dependent flowering time regulation in Arabidopsis.

An Atlas of Soybean Small RNAs Identifies Phased siRNAs from Hundreds of Coding Genes

An extensive analysis of small RNAs in soybean identified many miRNAs and phased, secondary siRNA (phasiRNA) loci; some of these miRNAs were the triggers of the phasiRNA loci. Small RNAs are ubiquitous, versatile repressors and include (1) microRNAs (miRNAs), processed from mRNA forming stem-loops; and (2) small interfering RNAs (siRNAs), the latter derived in plants by a process typically requiring an RNA-dependent RNA polymerase. We constructed and analyzed an expression atlas of soybean (Glycine max) small RNAs, identifying over 500 loci generating 21-nucleotide phased siRNAs (phasiRNAs; from PHAS loci), of which 483 overlapped annotated protein-coding genes. Via the integration of miRNAs with parallel analysis of RNA end (PARE) data, 20 miRNA triggers of 127 PHAS loci were detected. The primary class of PHAS loci (208 or 41% of the total) corresponded to NB-LRR genes; some of these small RNAs preferentially accumulate in nodules. Among the PHAS loci, novel representatives of TAS3 and noncanonical phasing patterns were also observed. A noncoding PHAS locus, triggered by miR4392, accumulated preferentially in anthers; the phasiRNAs are predicted to target transposable elements, with their peak abundance during soybean reproductive development. Thus, phasiRNAs show tremendous diversity in dicots. We identified novel miRNAs and assessed the veracity of soybean miRNAs registered in miRBase, substantially improving the soybean miRNA annotation, facilitating an improvement of miRBase annotations and identifying at high stringency novel miRNAs and their targets.