Regulation of Rice Tillering by RNA-Directed DNA Methylation at Miniature Inverted-Repeat Transposable Elements

Comparative and functional analysis unveils the contribution of photoperiod to DNA methylation, sRNA accumulation, and gene expression variations in short-day and long-day grasses

Photoperiod employs complicated networks to regulate various developmental processes in plants, including flowering transition. However, the specific mechanisms by which photoperiod affects epigenetic modifications and gene expression variations in plants remain elusive. In this study, we conducted a comprehensive analysis of DNA methylation, small RNA (sRNA) accumulation, and gene expressions under different daylengths in facultative long-day (LD) grass Brachypodium distachyon and short-day (SD) grass rice. Our results showed that while overall DNA methylation levels were minimally affected by different photoperiods, CHH methylation levels were repressed under their favorable light conditions, particularly in rice. We identified numerous differentially methylated regions (DMRs) that were influenced by photoperiod in both plant species. Apart from differential sRNA clusters, we observed alterations in the expression of key components of the RNA-directed DNA methylation pathway, DNA methyltransferases, and demethylases, which may contribute to the identified photoperiod-influenced CHH DMRs. Furthermore, we identified many differentially expressed genes in response to different daylengths, some of which were associated with the DMRs. Notably, we discovered a photoperiod-responsive gene MYB11 in the transcriptome of B. distachyon, and further demonstrated its role as a flowering inhibitor by repressing FT1 transcription. Together, our comparative and functional analysis sheds light on the effects of daylength on DNA methylation, sRNA accumulation, and gene expression variations in LD and SD plants, thereby facilitating better designing breeding programs aimed at developing high-yield crops that can adapt to local growing seasons.

A DNA demethylase reduces seed size by decreasing the DNA methylation of AT-rich transposable elements in soybean

Understanding how to increase soybean yield is crucial for global food security. The genetic and epigenetic factors influencing seed size, a major crop yield determinant, are not fully understood. We explore the role of DNA demethylase GmDMEa in soybean seed size. Our research indicates that GmDMEa negatively correlates with soybean seed size. Using CRISPR-Cas9, we edited GmDMEa in the Dongnong soybean cultivar, known for small seeds. Modified plants had larger seeds and greater yields without altering plant architecture or seed nutrition. GmDMEa preferentially demethylates AT-rich transposable elements, thus activating genes and transcription factors associated with the abscisic acid pathway, which typically decreases seed size. Chromosomal substitution lines confirm that these modifications are inheritable, suggesting a stable epigenetic method to boost seed size in future breeding. Our findings provide insights into epigenetic seed size control and suggest a strategy for improving crop yields through the epigenetic regulation of crucial genes. This work implies that targeted epigenetic modification has practical agricultural applications, potentially enhancing food production without compromising crop quality.

Transposable elements cause the loss of self-incompatibility in citrus

Self-incompatibility (SI) is a widespread prezygotic mechanism for flowering plants to avoid inbreeding depression and promote genetic diversity. Citrus has an S-RNase-based SI system, which was frequently lost during evolution. We previously identified a single nucleotide mutation in Sm-RNase, which is responsible for the loss of SI in mandarin and its hybrids. However, little is known about other mechanisms responsible for conversion of SI to self-compatibility (SC) and we identify a completely different mechanism widely utilized by citrus. Here, we found a 786-bp miniature inverted-repeat transposable element (MITE) insertion in the promoter region of the FhiS2-RNase in Fortunella hindsii Swingle (a model plant for citrus gene function), which does not contain the Sm-RNase allele but are still SC. We demonstrate that this MITE plays a pivotal role in the loss of SI in citrus, providing evidence that this MITE insertion prevents expression of the S-RNase; moreover, transgenic experiments show that deletion of this 786-bp MITE insertion recovers the expression of FhiS2-RNase and restores SI. This study identifies the first evidence for a role for MITEs at the S-locus affecting the SI phenotype. A family-wide survey of the S-locus revealed that MITE insertions occur frequently adjacent to S-RNase alleles in different citrus genera, but only certain MITEs appear to be responsible for the loss of SI. Our study provides evidence that insertion of MITEs into a promoter region can alter a breeding strategy and suggests that this phenomenon may be broadly responsible for SC in species with the S-RNase system.

High-quality genome assembly enables prediction of allele-specific gene expression in hybrid poplar

Poplar (Populus) is a well-established model system for tree genomics and molecular breeding, and hybrid poplar is widely used in forest plantations. However, distinguishing its diploid homologous chromosomes is difficult, complicating advanced functional studies on specific alleles. In this study, we applied a trio-binning design and PacBio high-fidelity long-read sequencing to obtain haplotype-phased telomere-to-telomere genome assemblies for the 2 parents of the well-studied F1 hybrid "84K" (Populus alba × Populus tremula var. glandulosa). Almost all chromosomes, including the telomeres and centromeres, were completely assembled for each haplotype subgenome apart from 2 small gaps on one chromosome. By incorporating information from these haplotype assemblies and extensive RNA-seq data, we analyzed gene expression patterns between the 2 subgenomes and alleles. Transcription bias at the subgenome level was not uncovered, but extensive-expression differences were detected between alleles. We developed machine-learning (ML) models to predict allele-specific expression (ASE) with high accuracy and identified underlying genome features most highly influencing ASE. One of our models with 15 predictor variables achieved 77% accuracy on the training set and 74% accuracy on the testing set. ML models identified gene body CHG methylation, sequence divergence, and transposon occupancy both upstream and downstream of alleles as important factors for ASE. Our haplotype-phased genome assemblies and ML strategy highlight an avenue for functional studies in Populus and provide additional tools for studying ASE and heterosis in hybrids.

Long noncoding RNAs underlie multiple domestication traits and leafhopper resistance in soybean

The origin and functionality of long noncoding RNA (lncRNA) remain poorly understood. Here, we show that multiple quantitative trait loci modulating distinct domestication traits in soybeans are pleiotropic effects of a locus composed of two tandem lncRNA genes. These lncRNA genes, each containing two inverted repeats, originating from coding sequences of the MYB genes, function in wild soybeans by generating clusters of small RNA (sRNA) species that inhibit the expression of their MYB gene relatives through post-transcriptional regulation. By contrast, the expression of lncRNA genes in cultivated soybeans is severely repressed, and, consequently, the corresponding MYB genes are highly expressed, shaping multiple distinct domestication traits as well as leafhopper resistance. The inverted repeats were formed before the divergence of the Glycine genus from the Phaseolus-Vigna lineage and exhibit strong structure-function constraints. This study exemplifies a type of target for selection during plant domestication and identifies mechanisms of lncRNA formation and action.

Discovery and genome-wide characterization of a novel miniature inverted repeat transposable element reveal genome-specific distribution in Glycine

BACKGROUND

Miniature inverted repeat transposable elements (MITEs) are a dynamic component responsible for genome evolution. Tourist MITEs are derived from and mobilized by elements from the harbinger superfamily.

OBJECTIVE

In this study, a novel family of Tourist-like MITE was characterized in wild soybean species Glycine falcata. The new GftoMITE1 was initially discovered as an insertional polymorphism of the centromere-specific histone H3 (CenH3) gene in G. falcata.

METHODS

Using polymerase chain reaction, cloning and sequencing approaches, we showed a high number of copies of the GftoMITE1 family. Extensive bioinformatic analyses revealed the genome-level distribution and locus-specific mapping of GftoMITE1 members in Glycine species.

RESULTS

Our results provide the first extensive characterization of the GftoMITE1 family and contribute to the understanding of the evolution of MITEs in the Glycine genus. Genome-specific GftoMITE1 was prominent in perennial wild soybean species, but not in annual cultivated soybean (Glycine max) or its progenitor (Glycine soja).

CONCLUSIONS

We discuss that the GftoMITE1 family reveals a single rapid amplification in G. falcata and could have potential implications for gene regulation and soybean breeding as an efficient genetic marker for germplasm utilization in the future.

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.

Transposition of HOPPLA in siRNA-deficient plants suggests a limited effect of the environment on retrotransposon mobility in Brachypodium distachyon

Long terminal repeat retrotransposons (LTR-RTs) are powerful mutagens regarded as a major source of genetic novelty and important drivers of evolution. Yet, the uncontrolled and potentially selfish proliferation of LTR-RTs can lead to deleterious mutations and genome instability, with large fitness costs for their host. While population genomics data suggest that an ongoing LTR-RT mobility is common in many species, the understanding of their dual role in evolution is limited. Here, we harness the genetic diversity of 320 sequenced natural accessions of the Mediterranean grass Brachypodium distachyon to characterize how genetic and environmental factors influence plant LTR-RT dynamics in the wild. When combining a coverage-based approach to estimate global LTR-RT copy number variations with mobilome-sequencing of nine accessions exposed to eight different stresses, we find little evidence for a major role of environmental factors in LTR-RT accumulations in B. distachyon natural accessions. Instead, we show that loss of RNA polymerase IV (Pol IV), which mediates RNA-directed DNA methylation in plants, results in high transcriptional and transpositional activities of RLC_BdisC024 (HOPPLA) LTR-RT family elements, and that these effects are not stress-specific. This work supports findings indicating an ongoing mobility in B. distachyon and reveals that host RNA-directed DNA methylation rather than environmental factors controls their mobility in this wild grass model.

Insertion of a miniature inverted-repeat transposable element into the promoter of OsTCP4 results in more tillers and a lower grain size in rice

A class I PCF type protein, TCP4, was identified as a transcription factor associated with both grain size and tillering through a DNA pull-down-MS assay combined with a genome-wide association study. This transcription factor was found to have a significant role in the variations among the 533 rice accessions, dividing them into two main subspecies. A Tourist-like miniature inverted-repeat transposable element (MITE) was discovered in the promoter of TCP4 in japonica/geng accessions (TCP4M+), which was found to suppress the expression of TCP4 at the transcriptional level. The MITE-deleted haplotype (TCP4M-) was mainly found in indica/xian accessions. ChIP-qPCR and EMSA demonstrated the binding of TCP4 to promoters of grain reservoir genes such as SSIIa and Amy3D in vivo and in vitro, respectively. The introduction of the genomic sequence of TCP4M+ into different TCP4M- cultivars was found to affect the expression of TCP4 in the transgenic rice, resulting in decreased expression of its downstream target gene SSIIa, increased tiller number, and decreased seed length. This study revealed that a Tourist-like MITE contributes to subspecies divergence by regulating the expression of TCP4 in response to environmental pressure, thus influencing source-sink balance by regulating starch biosynthesis in rice.

Genomic asymmetric epigenetic modification of transposable elements is involved in gene expression regulation of allopolyploid Brassica napus

Polyploids are common and have a wide geographical distribution and environmental adaptability. Allopolyploidy may lead to the activation of transposable elements (TE). However, the mechanism of epigenetic modification of TEs in the establishment and evolution of allopolyploids remains to be explored. We focused on the TEs of model allopolyploid Brassica napus (An An Cn Cn ), exploring the TE characteristics of the genome, epigenetic modifications of TEs during allopolyploidization, and regulation of gene expression by TE methylation. In B. napus, approximately 50% of the genome was composed of TEs. TEs increased with proximity to genes, especially DNA transposons. TE methylation levels were negatively correlated with gene expression, and changes in TE methylation levels were able to regulate the expression of neighboring genes related to responses to light intensity and stress, which promoted powerful adaptation of allopolyploids to new environments. TEs can be synergistically regulated by RNA-directed DNA methylation pathways and histone modifications. The epigenetic modification levels of TEs tended to be similar to those of the diploid parents during the genome evolution of B. napus. The TEs of the An subgenome were more likely to be modified, and the imbalance in TE number and epigenetic modification level in the An and Cn subgenomes may lead to the establishment of subgenome dominance. Our study analyzed the characteristics of TE location, DNA methylation, siRNA, and histone modification in B. napus and highlighted the importance of TE epigenetic modifications during the allopolyploidy process, providing support for revealing the mechanism of allopolyploid formation and evolution.

CLSY docking to Pol IV requires a conserved domain critical for small RNA biogenesis and transposon silencing

Eukaryotes must balance the need for gene transcription by RNA polymerase II (Pol II) against the danger of mutations caused by transposable element (TE) proliferation. In plants, these gene expression and TE silencing activities are divided between different RNA polymerases. Specifically, RNA polymerase IV (Pol IV), which evolved from Pol II, transcribes TEs to generate small interfering RNAs (siRNAs) that guide DNA methylation and block TE transcription by Pol II. While the Pol IV complex is recruited to TEs via SNF2-like CLASSY (CLSY) proteins, how Pol IV partners with the CLSYs remains unknown. Here we identified a conserved CYC-YPMF motif that is specific to Pol IV and is positioned on the complex exterior. Furthermore, we found that this motif is essential for the co-purification of all four CLSYs with Pol IV, but that only one CLSY is present in any given Pol IV complex. These findings support a "one CLSY per Pol IV" model where the CYC-YPMF motif acts as a CLSY-docking site. Indeed, mutations in and around this motif phenocopy pol iv null mutants. Together, these findings provide structural and functional insights into a critical protein feature that distinguishes Pol IV from other RNA polymerases, allowing it to promote genome stability by targeting TEs for silencing.

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.

Rice transformation treatments leave specific epigenome changes beyond tissue culture

During transgenic plant production, tissue culture often carries epigenetic, and genetic changes that underlie somaclonal variations, leading to unpredictable phenotypes. Additionally, specific treatments for rice (Oryza sativa) transformation processes may individually or jointly contribute to somaclonal variations, but their specific impacts on rice epigenomes toward transcriptional variations remain unknown. Here, the impact of individual transformation treatments on genome-wide DNA methylation and the transcriptome were examined. In addition to activating stress-responsive genes, individual transformation components targeted different gene expression modules that were enriched in specific functional categories. The transformation treatments strongly impacted DNA methylation and expression; 75% were independent of tissue culture. Furthermore, our genome-wide analysis showed that the transformation treatments consistently resulted in global hypo-CHH methylation enriched at promoters highly associated with downregulation, particularly when the promoters were colocalized with miniature inverted-repeat transposable elements. Our results clearly highlight the specificity of impacts triggered by individual transformation treatments during rice transformation with the potential association between DNA methylation and gene expression. These changes in gene expression and DNA methylation resulting from rice transformation treatments explain a significant portion of somaclonal variations, that is, way beyond the tissue culture effect.

Beyond transcription: compelling open questions in plant RNA biology

The study of RNAs has become one of the most influential research fields in contemporary biology and biomedicine. In the last few years, new sequencing technologies have produced an explosion of new and exciting discoveries in the field but have also given rise to many open questions. Defining these questions, together with old, long-standing gaps in our knowledge, is the spirit of this article. The breadth of topics within RNA biology research is vast, and every aspect of the biology of these molecules contains countless exciting open questions. Here, we asked 12 groups to discuss their most compelling question among some plant RNA biology topics. The following vignettes cover RNA alternative splicing; RNA dynamics; RNA translation; RNA structures; R-loops; epitranscriptomics; long non-coding RNAs; small RNA production and their functions in crops; small RNAs during gametogenesis and in cross-kingdom RNA interference; and RNA-directed DNA methylation. In each section, we will present the current state-of-the-art in plant RNA biology research before asking the questions that will surely motivate future discoveries in the field. We hope this article will spark a debate about the future perspective on RNA biology and provoke novel reflections in the reader.

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.

Geographic-genomic and geographic-phenotypic differentiation of the Aquilegia viridiflora complex

How species diverge into different lineages is a central issue in evolutionary biology. Despite the increasing evidence indicating that such divergences do not need geographic isolation, the correlation between lineage divergence and the adaptive ecological divergence of phenotype corresponding to distribution is still unknown. In addition, gene flow has been widely detected during and through such diverging processes. We used one widely distributed Aquilegia viridiflora complex as a model system to examine genomic differentiation and corresponding phenotypic variations along geographic gradients. Our phenotypic analyses of 20 populations from northwest to northeast China identified two phenotypic groups along the geographic cline. All examined traits are distinct from each other, although a few intermediate individuals occur in their contacting regions. We further sequenced the genomes of representative individuals of each population. However, four distinct genetic lineages were detected based on nuclear genomes. In particular, we recovered numerous genetic hybrids in the contact regions of four lineages. Gene flow is widespread and continuous between four lineages but much higher between contacting lineages than geographically isolated lineages. Gene flow and natural selection might result in inconsistency between heredity and phenotype. Moreover, many genes with fast lineage-specific mutations were identified to be involved in local adaptation. Our results suggest that both geographic isolation and local selection exerted by the environment and pollinators may together create geographic distributions of phenotypic variations as well as the underlying genomic divergences in numerous lineages.

Lost genome segments associate with trait diversity during rice domestication

BACKGROUND

DNA mutations of diverse types provide the raw material required for phenotypic variation and evolution. In the case of crop species, previous research aimed to elucidate the changing patterns of repetitive sequences, single-nucleotide polymorphisms (SNPs), and small InDels during domestication to explain morphological evolution and adaptation to different environments. Additionally, structural variations (SVs) encompassing larger stretches of DNA are more likely to alter gene expression levels leading to phenotypic variation affecting plant phenotypes and stress resistance. Previous studies on SVs in rice were hampered by reliance on short-read sequencing limiting the quantity and quality of SV identification, while SV data are currently only available for cultivated rice, with wild rice largely uncharacterized. Here, we generated two genome assemblies for O. rufipogon using long-read sequencing and provide insights on the evolutionary pattern and effect of SVs on morphological traits during rice domestication.

RESULTS

In this study, we identified 318,589 SVs in cultivated and wild rice populations through a comprehensive analysis of 13 high-quality rice genomes and found that wild rice genomes contain 49% of unique SVs and an average of 1.76% of genes were lost during rice domestication. These SVs were further genotyped for 649 rice accessions, their evolutionary pattern during rice domestication and potential association with the diversity of important agronomic traits were examined. Genome-wide association studies between these SVs and nine agronomic traits identified 413 candidate causal variants, which together affect 361 genes. An 824-bp deletion in japonica rice, which encodes a serine carboxypeptidase family protein, is shown to be associated with grain length.

CONCLUSIONS

We provide relatively accurate and complete SV datasets for cultivated and wild rice accessions, especially in TE-rich regions, by comparing long-read sequencing data for 13 representative varieties. The integrated rice SV map and the identified candidate genes and variants represent valuable resources for future genomic research and breeding in rice.

Polymorphic inverted repeats near coding genes impact chromatin topology and phenotypic traits in Arabidopsis thaliana

Transposons are mobile elements that are commonly silenced to protect eukaryotic genome integrity. In plants, transposable element (TE)-derived inverted repeats (IRs) are commonly found near genes, where they affect host gene expression. However, the molecular mechanisms of such regulation are unclear in most cases. Expression of these IRs is associated with production of 24-nt small RNAs, methylation of the IRs, and drastic changes in local 3D chromatin organization. Notably, many of these IRs differ between Arabidopsis thaliana accessions, causing variation in short-range chromatin interactions and gene expression. CRISPR-Cas9-mediated disruption of two IRs leads to a switch in genome topology and gene expression with phenotypic consequences. Our data show that insertion of an IR near a gene provides an anchor point for chromatin interactions that profoundly impact the activity of neighboring loci. This turns IRs into powerful evolutionary agents that can contribute to rapid adaptation.

Genome-wide analysis of epigenetic and transcriptional changes in the pathogenesis of RGSV in rice

Rice grassy stunt virus (RGSV), a typical negative single-stranded RNA virus, invades rice and generates several disease signs, including dwarfing, tillering, and sterility. Previous research has revealed that RGSV-encoded proteins can force the host's ubiquitin-proteasome system to utilize them for viral pathogenesis. However, most of the studies were limited to a single omics level and lacked multidimensional data collection and correlation analysis on the mechanisms of RGSV-rice interactions. Here, we performed a comprehensive association analysis of genome-wide methylation sequencing, transcriptome sequencing, and histone H3K9me3 modification in RGSV-infested as well as non-infested rice leaves, and the levels of all three cytosine contexts (CG, CHG and CHH) were found to be slightly lower in RGSV-infected rice leaves than in normal rice. Large proportions of DMRs were distributed in the promoter and intergenic regions, and most DMRs were enriched in the CHH context, where the number of CHH hypo-DMRs was almost twice as high as that of hyper-DMRs. Among the genes with down-regulated expression and hypermethylation, we analyzed and identified 11 transcripts involved in fertility, plant height and tillering, and among the transcribed up-regulated and hypermethylated genes, we excavated 7 transcripts related to fertility, plant height and tillering. By analyzing the changes of histone H3K9me3 modification before and after virus infestation, we found that the distribution of H3K9me3 modification in the whole rice genome was prevalent, mainly concentrated in the gene promoter and gene body regions, which was distinctly different from the characteristics of animals. Combined with transcriptomic data, H3K9me3 mark was found to favor targeting highly expressed genes. After RGSV infection, H3K9me3 modifications in several regions of CTK and BR hormone signaling-related genes were altered, providing important targets for subsequent studies.

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.

Evolution of lmiRNAs and their targets from MITEs for rice adaptation

Twenty-four nucleotide long microRNAs (lmiRNAs) direct DNA methylation at target genes and regulate their transcription. The evolutionary origin of lmiRNAs and the range of lmiRNA-mediated regulation remain obscure. Here, we reannotated lmiRNAs and their targets in rice by applying stringent criteria. We found that the majority of lmiRNAs are derived from Miniature Inverted-repeat Transposable Elements (MITEs) and most sites targeted by MITE-derived lmiRNAs reside within MITEs, suggesting co-evolution of lmiRNAs and their targets through MITE amplification. lmiRNAs undergo dynamically changes under stress conditions and the genes targeted by lmiRNAs show an enrichment for stress-responsive genes, suggesting that lmiRNAs are widely involved in plant responses to stresses. We constructed the evolutionary histories of lmiRNAs and their targets. Nearly half of lmiRNAs emerged before or when the AA genome was diverged, while the emergence of lmiRNA targets coincided with or followed the emergence of lmiRNAs. Furthermore, we found that the sequences of a lmiRNA target site underwent variations, coincident with the divergence of rice accessions and the distribution of rice accessions in different geographical locations and climatic conditions. Our findings highlight MITEs as an important origin of lmiRNAs and suggest that the evolution of lmiRNA-target regulatory modules may contribute to rice adaptation to environmental changes.

ZMP recruits and excludes Pol IV-mediated DNA methylation in a site-specific manner

In plants, RNA-directed DNA methylation (RdDM) uses small interfering RNAs (siRNAs) to target transposable elements (TEs) but usually avoids genes. RNA polymerase IV (Pol IV) shapes the landscape of DNA methylation through its pivotal role in siRNA biogenesis. However, how Pol IV is recruited to specific loci, particularly how it avoids genes, is poorly understood. Here, we identified a Pol IV-interacting protein, ZMP (zinc finger, mouse double-minute/switching complex B, Plus-3 protein), which exerts a dual role in regulating siRNA biogenesis and DNA methylation at specific genomic regions. ZMP is required for siRNA biogenesis at some pericentromeric regions and prevents Pol IV from targeting a subset of TEs and genes at euchromatic loci. As a chromatin-associated protein, ZMP prefers regions with depleted histone H3 lysine 4 (H3K4) methylation abutted by regions with H3K4 methylation, probably monitoring changes in local H3K4 methylation status to regulate Pol IV's chromatin occupancy. Our findings uncover a mechanism governing the specificity of RdDM.

Will epigenetics be a key player in crop breeding?

If food and feed production are to keep up with world demand in the face of climate change, continued progress in understanding and utilizing both genetic and epigenetic sources of crop variation is necessary. Progress in plant breeding has traditionally been thought to be due to selection for spontaneous DNA sequence mutations that impart desirable phenotypes. These spontaneous mutations can expand phenotypic diversity, from which breeders can select agronomically useful traits. However, it has become clear that phenotypic diversity can be generated even when the genome sequence is unaltered. Epigenetic gene regulation is a mechanism by which genome expression is regulated without altering the DNA sequence. With the development of high throughput DNA sequencers, it has become possible to analyze the epigenetic state of the whole genome, which is termed the epigenome. These techniques enable us to identify spontaneous epigenetic mutations (epimutations) with high throughput and identify the epimutations that lead to increased phenotypic diversity. These epimutations can create new phenotypes and the causative epimutations can be inherited over generations. There is evidence of selected agronomic traits being conditioned by heritable epimutations, and breeders may have historically selected for epiallele-conditioned agronomic traits. These results imply that not only DNA sequence diversity, but the diversity of epigenetic states can contribute to increased phenotypic diversity. However, since the modes of induction and transmission of epialleles and their stability differ from that of genetic alleles, the importance of inheritance as classically defined also differs. For example, there may be a difference between the types of epigenetic inheritance important to crop breeding and crop production. The former may depend more on longer-term inheritance whereas the latter may simply take advantage of shorter-term phenomena. With the advances in our understanding of epigenetics, epigenetics may bring new perspectives for crop improvement, such as the use of epigenetic variation or epigenome editing in breeding. In this review, we will introduce the role of epigenetic variation in plant breeding, largely focusing on DNA methylation, and conclude by asking to what extent new knowledge of epigenetics in crop breeding has led to documented cases of its successful use.

Methylation of a MITE insertion in the MdRFNR1-1 promoter is positively associated with its allelic expression in apple in response to drought stress

Miniature inverted-repeat transposable elements (MITEs) are widely distributed in the plant genome and can be methylated. However, whether DNA methylation of MITEs is associated with induced allelic expression and drought tolerance is unclear. Here, we identified the drought-inducible MdRFNR1 (root-type ferredoxin-NADP+ oxidoreductase) gene in apple (Malus domestica). MdRFNR1 plays a positive role in drought tolerance by regulating the redox system, including increasing NADP+ accumulation and catalase and peroxidase activities and decreasing NADPH levels. Sequence analysis identified a MITE insertion (MITE-MdRF1) in the promoter of MdRFNR1-1 but not the MdRFNR1-2 allele. MdRFNR1-1 but not MdRFNR1-2 expression was significantly induced by drought stress, which was positively associated with the MITE-MdRF1 insertion and its DNA methylation. The methylated MITE-MdRF1 is recognized by the transcriptional anti-silencing factors MdSUVH1 and MdSUVH3, which recruit the DNAJ domain-containing proteins MdDNAJ1, MdDNAJ2, and MdDNAJ5, thereby activating MdRFNR1-1 expression under drought stress. Finally, we showed that MdSUVH1 and MdDNAJ1 are positive regulators of drought tolerance. These findings illustrate the molecular roles of methylated MITE-MdRF1 (which is recognized by the MdSUVH-MdDNAJ complex) in induced MdRFNR1-1 expression as well as the drought response of apple and shed light on the molecular mechanisms of natural variation in perennial trees.

Dwarf and High Tillering1 represses rice tillering through mediating the splicing of D14 pre-mRNA

Strigolactones (SLs) constitute a class of plant hormones that regulate many aspects of plant development, including repressing tillering in rice (Oryza sativa). However, how SL pathways are regulated is still poorly understood. Here, we describe a rice mutant dwarf and high tillering1 (dht1), which exhibits pleiotropic phenotypes (such as dwarfism and increased tiller numbers) similar to those of mutants defective in SL signaling. We show that DHT1 encodes a monocotyledon-specific hnRNP-like protein that acts as a previously unrecognized intron splicing factor for many precursor mRNAs (pre-mRNAs), including for the SL receptor gene D14. We find that the dht1 (DHT1I232F) mutant protein is impaired in its stability and RNA binding activity, causing defective splicing of D14 pre-mRNA and reduced D14 expression, and consequently leading to the SL signaling-defective phenotypes. Overall, our findings deepen our understanding of the functional diversification of hnRNP-like proteins and establish a connection between posttranscriptional splicing and SL signaling in the regulation of plant development.

A Non-Canonical MITE in the WOX11 Promoter Is Associated with Robust Crown Root development in Rice

The formation of tissues and organs in multicellular organisms is tightly controlled by transcriptional programs determined by temporal and spatial patterns of gene expression. As an important regulator of rice crown root development, WOX11 is essential for crown root formation and its transcript level is positively correlated with crown root biomass. However, how WOX11 is regulated during crown root primordium emergence and outgrowth still remains unknown. In this study, variations of the WOX11 genomic sequence were analyzed, and the highest genetic diversity was found within its promoter, which contained a non-canonical miniature inverted-repeat transposable element (ncMITE) sequence. Analysis of the WOX11 promoter-driven reporter gene GUS (β-glucuronidase) transgenic plants pWOX11(ncMITE+):GUS and pWOX11(ncMITE-):GUS uncovered higher GUS expression levels in crown roots of pWOX11(ncMITE+):GUS plants. Furthermore, pWOX11(ncMITE+):WOX11-FLAG in wox11 background could complement the crown root number and length compared to those of the wild type, while pWOX11(ncMITE-):WOX11-FLAG could not. These results suggested that the ncMITE was positively associated with WOX11 transcripts in rice crown roots. In addition, DNA methylation nearby the ncMITE region attenuated the activation effect of the ncMITE on WOX11 expression, which might also be the cause conferred to the root-specific expression of WOX11. This work provides novel insight into WOX11 expression regulation and reveals a promising target for genetic improvement of root architecture in rice.

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.

Small regulatory RNAs in rice epigenetic regulation

Plant small RNAs (sRNAs) are short non-coding RNAs that are implicated in various regulatory processes involving post-transcriptional gene silencing and epigenetic gene regulation. In epigenetic regulation, sRNAs are primarily involved in RNA-directed DNA methylation (RdDM) pathways. sRNAs in the RdDM pathways play a role not only in the suppression of transposable element (TE) activity but also in gene expression regulation. Although the major components of the RdDM pathways have been well studied in Arabidopsis, recent studies have revealed that the RdDM pathways in rice have important biological functions in stress response and developmental processes. In this review, we summarize and discuss recent literature on sRNA-mediated epigenetic regulation in rice. First, we describe the RdDM mechanisms in plants. We then introduce recent discoveries on the biological roles of rice genes involved in the RdDM pathway and TE-derived sRNAs working at specific genomic loci for epigenetic control in rice.

TaMOR is essential for root initiation and improvement of root system architecture in wheat

Optimal root system architecture is beneficial for water-fertilizer use efficiency, stress tolerance and yield improvement of crops. However, because of the complexity of root traits and difficulty in phenotyping deep roots, the study on mechanisms of root development is rarely reported in wheat (Triticum aestivum L.). In this study, we identified that the LBD (LATERAL ORGAN BOUNDARIES DOMAIN) gene TaMOR (MORE ROOT in wheat) determines wheat crown root initiation. The mor mutants exhibited less or even no crown root, dwarfism, less grain number and lodging caused by few roots. The observation of cross sections showed that crown root initiation is inhibited in the mor mutants. Molecular assays revealed that TaMOR interacts with the auxin response factor ARF5 to directly induce the expression of the auxin transporter gene PIN2 (PIN-FORMED 2) in the root base to regulate crown root initiation. In addition, a 159-bp MITE (miniature inverted-repeat transposable element) insertion causing DNA methylation and lower expression of TaMOR-B was identified in TaMOR-B promoter, which is associated with lower root dry weight and shorter plant height. The results bring new light into regulation mechanisms of crown root initiation and offer a new target for the improvement of root system architecture in wheat.

A MITE variation-associated heat-inducible isoform of a heat-shock factor confers heat tolerance through regulation of JASMONATE ZIM-DOMAIN genes in rice

High temperatures cause huge yield losses in rice. Heat-shock factors (Hsfs) are key transcription factors which regulate the expression of heat stress-responsive genes, but natural variation in and functional characterization of Hsfs have seldom been reported. A significant heat response locus was detected via a genome-wide association study (GWAS) using green leaf area as an indicative trait. A miniature inverted-repeat transposable element (MITE) in the promoter of a candidate gene, HTG3 (heat-tolerance gene on chromosome 3), was found to be significantly associated with heat-induced expression of HTG3 and heat tolerance (HT). The MITE-absent variant has been selected in heat-prone rice-growing regions. HTG3a is an alternatively spliced isoform encoding a functional Hsf, and experiments using overexpression and knockout rice lines showed that HTG3a positively regulates HT at both vegetative and reproductive stages. The HTG3-regulated genes were enriched for heat shock proteins and jasmonic acid signaling. Two heat-responsive JASMONATE ZIM-DOMAIN (JAZ) genes were confirmed to be directly upregulated by HTG3a, and one of them, OsJAZ9, positively regulates HT. We conclude that HTG3 plays an important role in HT through the regulation of JAZs and other heat-responsive genes. The MITE-absent allele may be valuable for HT breeding in rice.

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.

The CLASSY family controls tissue-specific DNA methylation patterns in Arabidopsis

DNA methylation shapes the epigenetic landscape of the genome, plays critical roles in regulating gene expression, and ensures transposon silencing. As is evidenced by the numerous defects associated with aberrant DNA methylation landscapes, establishing proper tissue-specific methylation patterns is critical. Yet, how such differences arise remains a largely open question in both plants and animals. Here we demonstrate that CLASSY1-4 (CLSY1-4), four locus-specific regulators of DNA methylation, also control tissue-specific methylation patterns, with the most striking pattern observed in ovules where CLSY3 and CLSY4 control DNA methylation at loci with a highly conserved DNA motif. On a more global scale, we demonstrate that specific clsy mutants are sufficient to shift the epigenetic landscape between tissues. Together, these findings reveal substantial epigenetic diversity between tissues and assign these changes to specific CLSY proteins, elucidating how locus-specific targeting combined with tissue-specific expression enables the CLSYs to generate epigenetic diversity during plant development.

CLASSY (CLSY) proteins regulate DNA methylation at specific loci in the Arabidopsis genome. Here the authors show that the CLSYs also control tissue-specific DNA methylation, including at siren loci in ovules, and that the lack of an individual CLSYs can shift the epigenetic landscape between tissues.

CG and CHG Methylation Contribute to the Transcriptional Control of OsPRR37-Output Genes in Rice

Plant circadian clock coordinates endogenous transcriptional rhythms with diurnal changes of environmental cues. OsPRR37, a negative component in the rice circadian clock, reportedly regulates transcriptome rhythms, and agronomically important traits. However, the underlying regulatory mechanisms of OsPRR37-output genes remain largely unknown. In this study, whole genome bisulfite sequencing and high-throughput RNA sequencing were applied to verify the role of DNA methylation in the transcriptional control of OsPRR37-output genes. We found that the overexpression of OsPRR37 suppressed rice growth and altered cytosine methylations in CG and CHG sequence contexts in but not the CHH context (H represents A, T, or C). In total, 35 overlapping genes were identified, and 25 of them showed negative correlation between the methylation level and gene expression. The promoter of the hexokinase gene OsHXK1 was hypomethylated at both CG and CHG sites, and the expression of OsHXK1 was significantly increased. Meanwhile, the leaf starch content was consistently lower in OsPRR37 overexpression lines than in the recipient parent Guangluai 4. Further analysis with published data of time-course transcriptomes revealed that most overlapping genes showed peak expression phases from dusk to dawn. The genes involved in DNA methylation, methylation maintenance, and DNA demethylation were found to be actively expressed around dusk. A DNA glycosylase, namely ROS1A/DNG702, was probably the upstream candidate that demethylated the promoter of OsHXK1. Taken together, our results revealed that CG and CHG methylation contribute to the transcriptional regulation of OsPRR37-output genes, and hypomethylation of OsHXK1 leads to decreased starch content and reduced plant growth in rice.

OsRLR4 binds to the OsAUX1 promoter to negatively regulate primary root development in rice

Root architecture is one of the most important agronomic traits that determines rice crop yield. The primary root (PR) absorbs mineral nutrients and provides mechanical support; however, the molecular mechanisms of PR elongation remain unclear in rice. Here, the two loss-of-function T-DNA insertion mutants of root length regulator 4 (OsRLR4), osrlr4-1 and osrlr4-2 with longer PR, and three OsRLR4 overexpression lines, OE-OsRLR4-1/-2/-3 with shorter PR compared to the wild type/Hwayoung (WT/HY), were identified. OsRLR4 is one of five members of the PRAF subfamily of the regulator chromosome condensation 1 (RCC1) family. Phylogenetic analysis of OsRLR4 from wild and cultivated rice indicated that it is under selective sweeps, suggesting its potential role in domestication. OsRLR4 controls PR development by regulating auxin accumulation in the PR tip and thus the root apical meristem activity. A series of biochemical and genetic analyses demonstrated that OsRLR4 functions directly upstream of the auxin transporter OsAUX1. Moreover, OsRLR4 interacts with the TRITHORAX-like protein OsTrx1 to promote H3K4me3 deposition at the OsAUX1 promoter, thus altering its transcription level. This work provides insight into the cooperation of auxin and epigenetic modifications in regulating root architecture and provides a genetic resource for plant architecture breeding.

Intragenic heterochromatin-mediated alternative polyadenylation modulates miRNA and pollen development in rice

Despite a much higher proportion of intragenic heterochromatin-containing genes in crop genomes, the importance of intragenic heterochromatin in crop development remains unclear. Intragenic heterochromatin can be recognised by a protein complex, ASI1-AIPP1-EDM2 (AAE) complex, to regulate alternative polyadenylation. Here, we investigated the impact of rice ASI1 on global poly(A) site usage through poly(A) sequencing and ASI1-dependent regulation on rice development. We found that OsASI1 is essential for rice pollen development and flowering. OsASI1 dysfunction has an important impact on global poly(A) site usage, which is closely related to heterochromatin marks. Intriguingly, OsASI1 interacts with the intronic heterochromatin of OsXRNL, a nuclear XRN family exonuclease gene involved in the processing of an miRNA precursor, to promote the processing of full-length OsXRNL and regulate miRNA abundance. We found that OsASI1-mediated regulation of pollen development partially depends on OsXRNL. Finally, we characterised the rice AAE complex and its involvement in alternative polyadenylation and pollen development. Our findings help to elucidate an epigenetic mechanism governing miRNA abundance and rice development, and provide a valuable resource for studying the epigenetic mechanisms of many important processes in crops.

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

RNA-directed DNA methylation (RdDM) functions in de novo methylation in CG, CHG, and CHH contexts. Here, we performed map-based cloning of OsNRPE1, which encodes the largest subunit of RNA polymerase V (Pol V), a key regulator of gene silencing and reproductive development in rice. We found that rice Pol V is required for CHH methylation on RdDM loci by transcribing long noncoding RNAs. Pol V influences the accumulation of 24-nucleotide small interfering RNAs (24-nt siRNAs) in a locus-specific manner. Biosynthesis of 24-nt siRNAs on loci with high CHH methylation levels and low CG and CHG methylation levels tends to depend on Pol V. In contrast, low methylation levels in the CHH context and high methylation levels in CG and CHG contexts predisposes 24-nt siRNA accumulation to be independent of Pol V. H3K9me1 and H3K9me2 tend to be enriched on Pol V-independent 24-nt siRNA loci, whereas various active histone modifications are enriched on Pol V-dependent 24-nt siRNA loci. DNA methylation is required for 24-nt siRNAs biosynthesis on Pol V-dependent loci but not on Pol V-independent loci. Our results reveal the function of rice Pol V for long noncoding RNA production, DNA methylation, 24-nt siRNA accumulation, and reproductive development.

Amplification dynamics of miniature inverted-repeat transposable elements and their impact on rice trait variability

Transposable elements (TEs) are a rich source of genetic variability. Among TEs, miniature inverted-repeat TEs (MITEs) are of particular interest as they are present in high copy numbers in plant genomes and are closely associated with genes. MITEs are deletion derivatives of class II transposons, and can be mobilized by the transposases encoded by the latter through a typical cut-and-paste mechanism. However, MITEs are typically present at much higher copy numbers than class II transposons. We present here an analysis of 103 109 transposon insertion polymorphisms (TIPs) in 738 Oryza sativa genomes representing the main rice population groups. We show that an important fraction of MITE insertions has been fixed in rice concomitantly with its domestication. However, another fraction of MITE insertions is present at low frequencies. We performed MITE TIP-genome-wide association studies (TIP-GWAS) to study the impact of these elements on agronomically important traits and found that these elements uncover more trait associations than single nucleotide polymorphisms (SNPs) on important phenotypes such as grain width. Finally, using SNP-GWAS and TIP-GWAS we provide evidence of the replicative amplification of MITEs.

A Global Landscape of Miniature Inverted-Repeat Transposable Elements in the Carrot Genome

Miniature inverted-repeat transposable elements (MITEs) are the most abundant group of Class II mobile elements in plant genomes. Their presence in genic regions may alter gene structure and expression, providing a new source of functional diversity. Owing to their small size and lack of coding capacity, the identification of MITEs has been demanding. However, the increasing availability of reference genomes and bioinformatic tools provides better means for the genome-wide identification and analysis of MITEs and for the elucidation of their contribution to the evolution of plant genomes. We mined MITEs in the carrot reference genome DH1 using MITE-hunter and developed a curated carrot MITE repository comprising 428 families. Of the 31,025 MITE copies spanning 10.34 Mbp of the carrot genome, 54% were positioned in genic regions. Stowaways and Tourists were frequently present in the vicinity of genes, while Mutator-like MITEs were relatively more enriched in introns. hAT-like MITEs were relatively more frequently associated with transcribed regions, including untranslated regions (UTRs). Some carrot MITE families were shared with other Apiaceae species. We showed that hAT-like MITEs were involved in the formation of new splice variants of insertion-harboring genes. Thus, carrot MITEs contributed to the accretion of new diversity by altering transcripts and possibly affecting the regulation of many genes.

Modification of cereal plant architecture by genome editing to improve yields

We summarize recent genome editing studies that have focused on the examination (or reexamination) of plant architectural phenotypes in cereals and the modification of these traits for crop improvement. Plant architecture is defined as the three-dimensional organization of the entire plant. Shoot architecture refers to the structure and organization of the aboveground components of a plant, reflecting the developmental patterning of stems, branches, leaves and inflorescences/flowers. Root system architecture is essentially determined by four major shape parameters-growth, branching, surface area and angle. Interest in plant architecture has arisen from the profound impact of many architectural traits on agronomic performance, and the genetic and hormonal regulation of these traits which makes them sensitive to both selective breeding and agronomic practices. This is particularly important in staple crops, and a large body of literature has, therefore, accumulated on the control of architectural phenotypes in cereals, particularly rice due to its twin role as one of the world's most important food crops as well as a model organism in plant biology and biotechnology. These studies have revealed many of the molecular mechanisms involved in the regulation of tiller/axillary branching, stem height, leaf and flower development, root architecture and the grain characteristics that ultimately help to determine yield. The advent of genome editing has made it possible, for the first time, to introduce precise mutations into cereal crops to optimize their architecture and close in on the concept of the ideotype. In this review, we consider recent genome editing studies that have focused on the examination (or reexamination) of plant architectural phenotypes in cereals and the modification of these traits for crop improvement.

Genome-wide analysis of microRNA156 and its targets, the genes encoding SQUAMOSA promoter-binding protein-like (SPL) transcription factors, in the grass family Poaceae

MicroRNAs (miRNAs) are endogenous small non-coding RNAs that play an important role in post-transcriptional gene regulation in plants and animals by targeting messenger RNAs (mRNAs) for cleavage or repressing translation of specific mRNAs. The first miRNA identified in plants, miRNA156 (miR156), targets the SQUAMOSA promoter-binding protein-like (SPL) transcription factors, which play critical roles in plant phase transition, flower and plant architecture, and fruit development. We identified multiple copies of MIR156 and SPL in the rice, Brachypodium, sorghum, maize, and foxtail millet genomes. Sequence and chromosomal synteny analysis showed that both MIR156s and SPLs are conserved across species in the grass family. Analysis of expression data of the SPLs in eleven juvenile and adult rice tissues revealed that four non-miR156-targeted genes were highly expressed and three miR156-targeted genes were only slightly expressed in all tissues/developmental stages. The remaining SPLs were highly expressed in the juvenile stage, but their expression was lower in the adult stage. It has been proposed that under strong selective pressure, non-miR156-targeted mRNA may be able to re-structure to form a miRNA-responsive element. In our analysis, some non-miR156-targeted SPLs (SPL5/8/10) had gene structure and gene expression patterns similar to those of miR156-targeted genes, suggesting that they could diversify into miR156-targeted genes. DNA methylation profiles of SPLs and MIR156s in different rice tissues showed diverse methylation patterns, and hypomethylation of non-CG sites was observed in rice endosperm. Our findings suggested that MIR156s and SPLs had different origination and evolutionary mechanisms: the SPLs appear to have resulted from vertical evolution, whereas MIR156s appear to have resulted from strong evolutionary selection on mature sequences.

Embryo CHH hypermethylation is mediated by RdDM and is autonomously directed in Brassica rapa

BACKGROUND

RNA-directed DNA methylation (RdDM) initiates cytosine methylation in all contexts and maintains asymmetric CHH methylation. Mature plant embryos show one of the highest levels of CHH methylation, and it has been suggested that RdDM is responsible for this hypermethylation. Because loss of RdDM in Brassica rapa causes seed abortion, embryo methylation might play a role in seed development. RdDM is required in the maternal sporophyte, suggesting that small RNAs from the maternal sporophyte might translocate to the developing embryo, triggering DNA methylation that prevents seed abortion. This raises the question of whether embryo hypermethylation is autonomously regulated by the embryo itself or influenced by the maternal sporophyte.

RESULTS

Here, we demonstrate that B. rapa embryos are hypermethylated in both euchromatin and heterochromatin and that this process requires RdDM. Contrary to the current models, B. rapa embryo hypermethylation is not correlated with demethylation of the endosperm. We also show that maternal somatic RdDM is not sufficient for global embryo hypermethylation, and we find no compelling evidence for maternal somatic influence over embryo methylation at any locus. Decoupling of maternal and zygotic RdDM leads to successful seed development despite the loss of embryo CHH hypermethylation.

CONCLUSIONS

We conclude that embryo CHH hypermethylation is conserved, autonomously controlled, and not required for embryo development. Furthermore, maternal somatic RdDM, while required for seed development, does not directly influence embryo methylation patterns.

Single-cell profiling lights different cell trajectories in plants

The molecular mechanism of the maintenance and differentiation of plant stem cells is an eternal theme in studies on plant growth and development. Recent advances in single-cell RNA sequencing (scRNA-seq) methods have completely changed the understanding of cell heterogeneity and cell function, allowing research precision to identify the differentiation trajectory of stem cells maintained and differentiated at the cellular level. This review aimed to mainly discuss the novel insights provided by scRNA-seq for the maintenance and initiation of plant stem cells, cell differentiation, cell response to environmental changes, and improvement strategies for scRNA-seq. In addition, it highlighted additional perspectives beyond scRNA-seq, such as spatial transcriptomes, epigenomes, and single-cell multiomics, for a renewed understanding of stem cell maintenance and cell differentiation, thus providing potential targets and theoretical foundations for crop improvement.

Arms race between rice and viruses: a review of viral and host factors

Much is known about the molecular interactions between positive-strand RNA viruses and dicotyledon plants. However, many important viral pathogens of the monocotyledon rice crop contain negative-strand or double-strand RNA genomes. Recent studies have shown that virus-derived small-interfering RNAs (siRNAs), host microRNAs and phytohormones regulate antiviral responses in rice plants and that rice-infecting RNA viruses encode a diverse repertoire of multifunctional proteins with counter-defensive activities. Moreover, the interactions between viral virulence proteins and host susceptibility factors also shape the virus-rice arms race. This review will focus on these recent advances and discuss strategies and challenges in the translation of discoveries made on molecular virus-rice interactions into practical virus control measures.

Non-coding RNA polymerases that silence transposable elements and reprogram gene expression in plants

Multisubunit RNA polymerase (Pol) complexes are the core machinery for gene expression in eukaryotes. The enzymes Pol I, Pol II and Pol III transcribe distinct subsets of nuclear genes. This family of nuclear RNA polymerases expanded in terrestrial plants by the duplication of Pol II subunit genes. Two Pol II-related enzymes, Pol IV and Pol V, are highly specialized in the production of regulatory, non-coding RNAs. Pol IV and Pol V are the central players of RNA-directed DNA methylation (RdDM), an RNA interference pathway that represses transposable elements (TEs) and selected genes. Genetic and biochemical analyses of Pol IV/V subunits are now revealing how these enzymes evolved from ancestral Pol II to sustain non-coding RNA biogenesis in silent chromatin. Intriguingly, Pol IV-RdDM regulates genes that influence flowering time, reproductive development, stress responses and plant–pathogen interactions. Pol IV target genes vary among closely related taxa, indicating that these regulatory circuits are often species-specific. Data from crops like maize, rice, tomato and Brassicarapa suggest that dynamic repositioning of TEs, accompanied by Pol IV targeting to TE-proximal genes, leads to the reprogramming of plant gene expression over short evolutionary timescales.

A DNA Methylation Reader-Chaperone Regulator-Transcription Factor Complex Activates OsHKT1;5 Expression during Salinity Stress

Irrigated lands are increasingly salinized, which adversely affects agricultural productivity. To respond to high sodium (Na+) concentrations, plants harbor multiple Na+ transport systems. Rice (Oryza sativa) HIGH-AFFINITY POTASSIUM (K+) TRANSPORTER1;5 (OsHKT1;5), a Na+-selective transporter, maintains K+/Na+ homeostasis under salt stress. However, the mechanism regulating OsHKT1;5 expression remains unknown. Here, we present evidence that a protein complex consisting of rice BCL-2-ASSOCIATED ATHANOGENE4 (OsBAG4), OsMYB106, and OsSUVH7 regulates OsHKT1;5 expression in response to salt stress. We isolated a salt stress-sensitive mutant, osbag4-1, that showed significantly reduced OsHKT1;5 expression and reduced K+ and elevated Na+ levels in shoots. Using comparative interactomics, we isolated two OsBAG4-interacting proteins, OsMYB106 (a MYB transcription factor) and OsSUVH7 (a DNA methylation reader), that were crucial for OsHKT1;5 expression. OsMYB106 and OsSUVH7 bound to the MYB binding cis-element (MYBE) and the miniature inverted-repeat transposable element (MITE) upstream of the MYBE, respectively, in the OsHKT1;5 promoter. OsBAG4 functioned as a bridge between OsSUVH7 and OsMYB106 to facilitate OsMYB106 binding to the consensus MYBE in the OsHKT1;5 promoter, thereby activating the OsHKT1;5 expression. Elimination of the MITE or knockout of OsMYB106 or OsSUVH7 decreased OsHKT1;5 expression and increased salt sensitivity. Our findings reveal a transcriptional complex, consisting of a DNA methylation reader, a chaperone regulator, and a transcription factor, that collaboratively regulate OsHKT1;5 expression during salinity stress.

Multifaceted roles of RNA polymerase IV in plant growth and development

We discuss the latest findings on RNA polymerase IV (Pol IV) in plant growth and development, providing new insights and expanding on new ideas for further, more in-depth research on Pol IV.