The effect of RNA polymerase V on 24-nt siRNA accumulation depends on DNA methylation contexts and histone modifications in rice

Upstream regulator of genomic imprinting in rice endosperm is a small RNA-associated chromatin remodeler

Genomic imprinting is observed in endosperm, a placenta-like seed tissue, where transposable elements (TEs) and repeat-derived small RNAs (sRNAs) mediate epigenetic changes in plants. In imprinting, uniparental gene expression arises due to parent-specific epigenetic marks on one allele but not on the other. The importance of sRNAs and their regulation in endosperm development or in imprinting is poorly understood in crops. Here we show that a previously uncharacterized CLASSY (CLSY)-family chromatin remodeler named OsCLSY3 is essential for rice endosperm development and imprinting, acting as an upstream player in the sRNA pathway. Comparative transcriptome and genetic analysis indicated its endosperm-preferred expression and its likely paternal imprinted nature. These important features are modulated by RNA-directed DNA methylation (RdDM) of tandemly arranged TEs in its promoter. Upon perturbation of OsCLSY3 in transgenic lines, we observe defects in endosperm development and a loss of around 70% of all sRNAs. Interestingly, well-conserved endosperm-specific sRNAs (siren) that are vital for reproductive fitness in angiosperms are also dependent on OsCLSY3. We observed that many imprinted genes and seed development-associated genes are under the control of OsCLSY3. These results support an essential role of OsCLSY3 in rice endosperm development and imprinting, and propose similar regulatory strategies involving CLSY3 homologs among other cereals.

Candidate Genes and Favorable Haplotypes Associated with Iron Toxicity Tolerance in Rice

Iron (Fe) toxicity is a major issue adversely affecting rice production worldwide. Unfortunately, the physiological and genetic mechanisms underlying Fe toxicity tolerance in rice remain relatively unknown. In this study, we conducted a genome-wide association study using a diverse panel consisting of 551 rice accessions to identify genetic mechanisms and candidate genes associated with Fe toxicity tolerance. Of the 29 quantitative trait loci (QTL) for Fe toxicity tolerance detected on chromosomes 1, 2, 5, and 12, five (qSH_Fe5, qSFW_Fe2.3, qRRL5.1, qRSFW1.1, and qRSFW12) were selected to identify candidate genes according to haplotype and bioinformatics analyses. The following five genes were revealed as promising candidates: LOC_Os05g40160, LOC_Os05g40180, LOC_Os12g36890, LOC_Os12g36900, and LOC_Os12g36940. The physiological characteristics of rice accessions with contrasting Fe toxicity tolerance reflected the importance of reactive oxygen species-scavenging antioxidant enzymes and Fe homeostasis for mitigating the negative effects of Fe toxicity on rice. Our findings have clarified the genetic and physiological mechanisms underlying Fe toxicity tolerance in rice. Furthermore, we identified valuable genetic resources for future functional analyses and the development of Fe toxicity-tolerant rice varieties via marker-assisted selection.

Mating system is associated with seed phenotypes upon loss of RNA-directed DNA methylation in Brassicaceae

In plants, de novo DNA methylation is guided by 24-nt short interfering (si)RNAs in a process called RNA-directed DNA methylation (RdDM). Primarily targeted at transposons, RdDM causes transcriptional silencing and can indirectly influence expression of neighboring genes. During reproduction, a small number of siRNA loci are dramatically upregulated in the maternally derived seed coat, suggesting that RdDM might have a special function during reproduction. However, the developmental consequence of RdDM has been difficult to dissect because disruption of RdDM does not result in overt phenotypes in Arabidopsis (Arabidopsis thaliana), where the pathway has been most thoroughly studied. In contrast, Brassica rapa mutants lacking RdDM have a severe seed production defect, which is determined by the maternal sporophytic genotype. To explore the factors that underlie the different phenotypes of these species, we produced RdDM mutations in 3 additional members of the Brassicaceae family: Camelina sativa, Capsella rubella, and Capsella grandiflora. Among these 3 species, only mutations in the obligate outcrosser, C. grandiflora, displayed a seed production defect similar to Brassica rapa mutants, suggesting that mating system is a key determinant for reproductive phenotypes in RdDM mutants.

Rice requires a chromatin remodeler for Polymerase IV-small interfering RNA production and genomic immunity

Transgenes are often spontaneously silenced, which hinders the application of genetic modifications to crop breeding. While gene silencing has been extensively studied in Arabidopsis (Arabidopsis thaliana), the molecular mechanism of transgene silencing remains elusive in crop plants. We used rice (Oryza sativa) plants silenced for a 35S::OsGA2ox1 (Gibberellin 2-oxidase 1) transgene to isolate five elements mountain (fem) mutants showing restoration of transgene expression. In this study, we isolated multiple fem2 mutants defective in a homolog of Required to Maintain Repression 1 (RMR1) of maize (Zea mays) and CLASSY (CLSY) of Arabidopsis. In addition to failing to maintain transgene silencing, as occurs in fem3, in which mutation occurs in NUCLEAR RNA POLYMERASE E1 (OsNRPE1), the fem2 mutant failed to establish transgene silencing of 35S::OsGA2ox1. Mutation in FEM2 eliminated all RNA POLYMERASE IV (Pol-IV)-FEM1/OsRDR2 (RNA-DEPENDENT RNA POLYMERASE 2)-dependent small interfering RNAs (siRNAs), reduced DNA methylation on genome-wide scale in rice seedlings, caused pleiotropic developmental defects, and increased disease resistance. Simultaneous mutation in 2 FEM2 homologous genes, FEM2-Like 1 (FEL1) and FEL2, however, did not affect DNA methylation and rice development and disease resistance. The predominant expression of FEM2 over FEL1 and FEL2 in various tissues was likely caused by epigenetic states. Overexpression of FEL1 but not FEL2 partially rescued hypomethylation of fem2, indicating that FEL1 maintains the cryptic function. In summary, FEM2 is essential for establishing and maintaining gene silencing; moreover, FEM2 is solely required for Pol IV-FEM1 siRNA biosynthesis and de novo DNA methylation.

Barley AGO4 proteins show overlapping functionality with distinct small RNA-binding properties in heterologous complementation

KEY MESSAGE

Barley AGO4 proteins complement expressional changes of epigenetically regulated genes in Arabidopsis ago4-3 mutant and show a distinct affinity for the 5' terminal nucleotide of small RNAs, demonstrating functional conservation and divergence. The function of Argonaute 4 (AGO4) in Arabidopsis thaliana has been extensively characterized; however, its role in monocots, which have large genomes abundantly supplemented with transposable elements (TEs), remains elusive. The study of barley AGO4 proteins can provide insights into the conserved aspects of RNA-directed DNA methylation (RdDM) and could also have further applications in the field of epigenetics or crop improvement. Bioinformatic analysis of RNA sequencing data identified two active AGO4 genes in barley, HvAGO4a and HvAGO4b. These genes function similar to AtAGO4 in an Arabidopsis heterologous complementation system, primarily binding to 24-nucleotide long small RNAs (sRNAs) and triggering methylation at specific target loci. Like AtAGO4, HvAGO4B exhibits a preference for binding sRNAs with 5' adenine residue, while also accepting 5' guanine, uracil, and cytosine residues. In contrast, HvAGO4A selectively binds only sRNAs with a 5' adenine residue. The diverse binding capacity of barley AGO4 proteins is reflected in TE-derived sRNAs and in their varying abundance. Both barley AGO4 proteins effectively restore the levels of extrachromosomal DNA and transcript abundancy of the heat-activated ONSEN retrotransposon to those observed in wild-type Arabidopsis plants. Our study provides insight into the distinct binding specificities and involvement in TE regulation of barley AGO4 proteins in Arabidopsis by heterologous complementation.

DNA-dependent RNA polymerases in plants

DNA-dependent RNA polymerases (Pols) transfer the genetic information stored in genomic DNA to RNA in all organisms. In eukaryotes, the typical products of nuclear Pol I, Pol II, and Pol III are ribosomal RNAs, mRNAs, and transfer RNAs, respectively. Intriguingly, plants possess two additional Pols, Pol IV and Pol V, which produce small RNAs and long noncoding RNAs, respectively, mainly for silencing transposable elements. The five plant Pols share some subunits, but their distinct functions stem from unique subunits that interact with specific regulatory factors in their transcription cycles. Here, we summarize recent advances in our understanding of plant nucleus-localized Pols, including their evolution, function, structures, and transcription cycles.

Small RNA-mediated DNA methylation during plant reproduction

Reproductive tissues are a rich source of small RNAs, including several classes of short interfering (si)RNAs that are restricted to this stage of development. In addition to RNA polymerase IV-dependent 24-nt siRNAs that trigger canonical RNA-directed DNA methylation, abundant reproductive-specific siRNAs are produced from companion cells adjacent to the developing germ line or zygote and may move intercellularly before inducing methylation. In some cases, these siRNAs are produced via non-canonical biosynthesis mechanisms or from sequences with little similarity to transposons. While the precise role of these siRNAs and the methylation they trigger is unclear, they have been implicated in specifying a single megaspore mother cell, silencing transposons in the male germ line, mediating parental dosage conflict to ensure proper endosperm development, hypermethylation of mature embryos, and trans-chromosomal methylation in hybrids. In this review, we summarize the current knowledge of reproductive siRNAs, including their biosynthesis, transport, and function.

Plant polymerase IV sensitizes chromatin through histone modifications to preclude spread of silencing into protein-coding domains

Across eukaryotes, gene regulation is manifested via chromatin states roughly distinguished as heterochromatin and euchromatin. The establishment, maintenance, and modulation of the chromatin states is mediated using several factors including chromatin modifiers. However, factors that avoid the intrusion of silencing signals into protein-coding genes are poorly understood. Here we show that a plant specific paralog of RNA polymerase (Pol) II, named Pol IV, is involved in avoidance of facultative heterochromatic marks in protein-coding genes, in addition to its well-established functions in silencing repeats and transposons. In its absence, H3K27 trimethylation (me3) mark intruded the protein-coding genes, more profoundly in genes embedded with repeats. In a subset of genes, spurious transcriptional activity resulted in small(s) RNA production, leading to post-transcriptional gene silencing. We show that such effects are significantly pronounced in rice, a plant with a larger genome with distributed heterochromatin compared with Arabidopsis Our results indicate the division of labor among plant-specific polymerases, not just in establishing effective silencing via sRNAs and DNA methylation but also in influencing chromatin boundaries.

An OsPRMT5-OsAGO2/miR1875-OsHXK1 module regulates rice immunity to blast disease

Rice ARGONAUTE2 (OsAGO2) is a core component of the rice RNA-induced silencing complex (RISC), which is repressed by Magnaporthe oryzae (M. oryzae) infection. Whether and how OsAGO2-mediated gene silencing plays a role in rice blast resistance and which sRNAs participate in this process are unknown. Our results indicate that OsAGO2 is a key immune player that manipulates rice defense responses against blast disease. OsAGO2 associates with the 24-nt miR1875 and binds to the promoter region of HEXOKINASE1 (OsHXK1), which causes DNA methylation and leads to gene silencing. Our multiple genetic evidence showed that, without M. oryzae infection, OsAGO2/miR1875 RISC promoted OsHXK1 promoter DNA methylation and OsHXK1 silencing; after M. oryzae infection, the reduced OsAGO2/miR1875 led to a relatively activated OsHXK1 expression. OsHXK1 acts as a positive regulator of blast disease resistance that OsHXK1-OE rice exhibited enhanced resistance, whereas Cas9-Oshxk1 rice showed reduced resistance against M. oryzae infection. OsHXK1 may function through its sugar sensor activity as glucose induced defense-related gene expression and reactive oxygen species (ROS) accumulation in Nipponbare and OsHXK1-OE but not in Cas9-Oshxk1 rice. OsAGO2 itself is delicately regulated by OsPRMT5, which senses M. oryzae infection and attenuates OsAGO2-mediated gene silencing through OsAGO2 arginine methylation. Our study reveals an OsPRMT5-OsAGO2/miR1875-OsHXK1 regulatory module that fine tunes the rice defense response to blast disease.

Factor of DNA methylation 1 affects woodland strawberry plant stature and organ size via DNA methylation

RNA-directed DNA methylation (RdDM) is an epigenetic process that directs silencing to specific genomic regions and loci. The biological functions of RdDM are not well studied in horticultural plants. Here, we isolated the ethyl methane-sulfonate-induced mutant reduced organ size (ros) producing small leaves, flowers, and fruits in woodland strawberry (Fragaria vesca) due to reduced cell numbers compared with that in the wild-type (WT). The candidate mutation causes a premature stop codon in FvH4_6g28780, which shares high similarity to Arabidopsis (Arabidopsis thaliana) Factor of DNA Methylation1 (FDM1) encoding an RdDM pathway component and was named FveFDM1. Consistently, the fvefdm1CR mutants generated by CRISPR/Cas9 also produced smaller organs. Overexpressing FveFDM1 in an Arabidopsis fdm1-1 fdm2-1 double mutant restored DNA methylation at the RdDM target loci. FveFDM1 acts in a protein complex with its homolog Involved in De Novo 2 (FveIDN2). Furthermore, whole-genome bisulfite sequencing revealed that DNA methylation, especially in the CHH context, was remarkably reduced throughout the genome in fvefdm1. Common and specific differentially expressed genes were identified in different tissues of fvefdm1 compared to in WT tissues. DNA methylation and expression levels of several gibberellic acid (GA) biosynthesis and cell cycle genes were validated. Moreover, the contents of GA and auxin were substantially reduced in the young leaves of fvefdm1 compared to in the WT. However, exogenous application of GA and auxin could not recover the organ size of fvefdm1. In addition, expression levels of FveFDM1, FveIDN2, Nuclear RNA Polymerase D1 (FveNRPD1), Domains Rearranged Methylase 2 (FveDRM2), and cell cycle genes were greatly induced by GA treatment. Overall, our work demonstrated the critical roles of FveFDM1 in plant growth and development via RdDM-mediated DNA methylation in horticultural crops.

Expansion and contraction of small RNA and methylation machinery throughout plant evolution

The revolution in sequencing has created a wealth of plant genomes that can be mined to understand the evolution of biological complexity. Complexity is often driven by gene duplication, which allows paralogs to specialize in an activity of the ancestral gene or acquire novel functions. Angiosperms encode a variety of gene silencing pathways that share related machinery for small RNA biosynthesis and function. Recent phylogenetic analysis of these gene families plots the expansion, specialization, and occasional contraction of this core machinery. This analysis reveals the ancient origin of RNA-directed DNA Methylation in early land plants, or possibly their algal ancestors, as well as ongoing duplications that evolve novel small RNA pathways.

A null allele of the pol IV second subunit impacts stature and reproductive development in Oryza sativa

All eukaryotes possess three DNA-dependent RNA polymerases, Pols I-III, while land plants possess two additional polymerases, Pol IV and Pol V. Derived through duplication of Pol II subunits, Pol IV produces 24-nt short interfering RNAs that interact with Pol V transcripts to target de novo DNA methylation and silence transcription of transposons. Members of the grass family encode additional duplicated subunits of Pol IV and V, raising questions regarding the function of each paralog. In this study, we identify a null allele of the putative Pol IV second subunit, NRPD2, and demonstrate that NRPD2 is the sole subunit functioning with NRPD1 in small RNA production and CHH methylation in leaves. Homozygous nrpd2 mutants have neither gametophytic defects nor embryo lethality, although adult plants are dwarf and sterile.

The Rise and Fall of Billionaire siRNAs during Reproductive Development in Rice

RNA polymerase IV-dependent siRNAs, usually 24 nt in length, function in the RNA-directed DNA methylation that is responsible for de novo methylation in plants. We analyzed 24 nt siRNAs in inflorescences and found that among the 20,200 24 nt siRNA clusters, the top 0.81% highly expressed clusters accounted for more than 68% of the 24 nt siRNA reads in inflorescences. We named the highly expressed siRNAs as billionaire siRNAs (bill-siRNAs) and the less-expressed siRNAs as pauper siRNAs (pau-siRNAs). The bill-siRNAs in inflorescences are mainly derived from the ovary. Female gametes produced more bill-siRNAs than male gametes. In embryos and seedlings developed from fertilized egg cells, the bill-siRNAs from gametes disappeared. The endosperm, which develops from the fertilized central cell, also contained no bill-siRNAs from gametes but did contain newly and highly expressed siRNAs produced in different regions. In contrast, bill-siRNAs from the ovaries were maintained in the seed coat. The biosynthesis of bill-siRNAs in various tissues and cells is dependent on OsRDR2 (RNA-dependent RNA polymerase 2) and Pol IV (DNA-dependent RNA polymerase IV). Similar to the pau-siRNAs, the first base of bill-siRNAs is enriched at adenine, and bill-siRNAs can direct DNA methylation in various tissues.

A rice single cell transcriptomic atlas defines the developmental trajectories of rice floret and inflorescence meristems

Rice inflorescence development determines yield and relies on the activity of axillary meristems (AMs); however, high-resolution analysis of its early development is lacking. Here, we have used high-throughput single-cell RNA sequencing to profile 37 571 rice inflorescence cells and constructed a genome-scale gene expression resource covering the inflorescence-to-floret transition during early reproductive development. The differentiation trajectories of florets and AMs were reconstructed, and discrete cell types and groups of regulators in the highly heterogeneous young inflorescence were identified and then validated by in situ hybridization and with fluorescent marker lines. Our data demonstrate that a WOX transcription factor, DWARF TILLER1, regulates flower meristem activity, and provide evidence for the role of auxin in rice inflorescence branching by exploring the expression and biological role of the auxin importer OsAUX1. Our comprehensive transcriptomic atlas of early rice inflorescence development, supported by genetic evidence, provides single-cell-level insights into AM differentiation and floret development.

Reinforcement of CHH methylation through RNA-directed DNA methylation ensures sexual reproduction in rice

DNA methylation is an important epigenetic mark that regulates the expression of genes and transposons. RNA-directed DNA methylation (RdDM) is the main molecular pathway responsible for de novo DNA methylation in plants. Although the mechanism of RdDM has been well studied in Arabidopsis (Arabidopsis thaliana), most mutations in RdDM genes cause no remarkable developmental defects in Arabidopsis. Here, we isolated and cloned Five Elements Mountain 1 (FEM1), which encodes RNA-dependent RNA polymerase 2 (OsRDR2) in rice (Oryza sativa). Mutation in OsRDR2 abolished the accumulation of 24-nt small interfering RNAs, and consequently substantially decreased genome-wide CHH (H = A, C, or T) methylation. Moreover, male and female reproductive development was disturbed, which led to sterility in osrdr2 mutants. We discovered that OsRDR2-dependent DNA methylation may regulate the expression of multiple key genes involved in stamen development, meiosis, and pollen viability. In wild-type (WT) plants but not in osrdr2 mutants, genome-wide CHH methylation levels were greater in panicles, stamens, and pistils than in seedlings. The global increase of CHH methylation in reproductive organs of the WT was mainly explained by the enhancement of RdDM activity, which includes OsRDR2 activity. Our results, which revealed a global increase in CHH methylation through enhancement of RdDM activity in reproductive organs, suggest a crucial role for OsRDR2 in the sexual reproduction of rice.