Integrative analysis of reference epigenomes in 20 rice varieties

Genome-wide mapping of DNase I hypersensitive sites revealed differential chromatin accessibility and regulatory DNA elements under drought stress in rice cultivars

Drought stress (DS) is one of the major constraints limiting yield in crop plants including rice. Gene regulation under DS is largely governed by accessibility of the transcription factors (TFs) to their cognate cis-regulatory elements (CREs). In this study, we used DNase I hypersensitive assays followed by sequencing to identify the accessible chromatin regions under DS in a drought-sensitive (IR64) and a drought-tolerant (N22) rice cultivar. Our results indicated that DNase I hypersensitive sites (DHSs) were highly enriched at transcription start sites (TSSs) and numerous DHSs were detected in the promoter regions. DHSs were concurrent with epigenetic marks and the genes harboring DHSs in their TSS and promoter regions were highly expressed. In addition, DS induced changes in DHSs (∆DHSs) in TSS and promoter regions were positively correlated with upregulation of several genes involved in drought/abiotic stress response, those encoding TFs and located within drought-associated quantitative trait loci, much preferentially in the drought-tolerant cultivar. The CREs representing the binding sites of TFs involved in DS response were detected within the ∆DHSs, suggesting differential accessibility of TFs to their cognate sites under DS in different rice cultivars, which may be further deployed for enhancing drought tolerance in rice.

Multiplexed CRISPR-Cas9 mutagenesis of rice PSBS1 noncoding sequences for transgene-free overexpression

Understanding CRISPR-Cas9's capacity to produce native overexpression (OX) alleles would accelerate agronomic gains achievable by gene editing. To generate OX alleles with increased RNA and protein abundance, we leveraged multiplexed CRISPR-Cas9 mutagenesis of noncoding sequences upstream of the rice PSBS1 gene. We isolated 120 gene-edited alleles with varying non-photochemical quenching (NPQ) capacity in vivo-from knockout to overexpression-using a high-throughput screening pipeline. Overexpression increased OsPsbS1 protein abundance two- to threefold, matching fold changes obtained by transgenesis. Increased PsbS protein abundance enhanced NPQ capacity and water-use efficiency. Across our resolved genetic variation, we identify the role of 5'UTR indels and inversions in driving knockout/knockdown and overexpression phenotypes, respectively. Complex structural variants, such as the 252-kb duplication/inversion generated here, evidence the potential of CRISPR-Cas9 to facilitate significant genomic changes with negligible off-target transcriptomic perturbations. Our results may inform future gene-editing strategies for hypermorphic alleles and have advanced the pursuit of gene-edited, non-transgenic rice plants with accelerated relaxation of photoprotection.

Conserved functions of chromatin regulators in basal Archaeplastida

Chromatin is a dynamic network that regulates genome organization and gene expression. Different types of chromatin regulators are highly conserved among Archaeplastida, including unicellular algae, while some chromatin genes are only present in land plant genomes. Here, we review recent advances in understanding the function of conserved chromatin factors in basal land plants and algae. We focus on the role of Polycomb-group genes which mediate H3K27me3-based silencing and play a role in balancing gene dosage and regulating haploid-to-diploid transitions by tissue-specific repression of the transcription factors KNOX and BELL in many representatives of the green lineage. Moreover, H3K27me3 predominantly occupies repetitive elements which can lead to their silencing in a unicellular alga and basal land plants, while it covers mostly protein-coding genes in higher land plants. In addition, we discuss the role of nuclear matrix constituent proteins as putative functional lamin analogs that are highly conserved among land plants and might have an ancestral function in stress response regulation. In summary, our review highlights the importance of studying chromatin regulation in a wide range of organisms in the Archaeplastida.

Transgenerational epigenetic inheritance during plant evolution and breeding

Plants can program and reprogram their genomes to create genetic variation and epigenetic modifications, leading to phenotypic plasticity. Although consequences of genetic changes are comprehensible, the basis for transgenerational inheritance of epigenetic variation is elusive. This review addresses contributions of external (environmental) and internal (genomic) factors to the establishment and maintenance of epigenetic memory during plant evolution, crop domestication, and modern breeding. Dynamic and pervasive changes in DNA methylation and chromatin modifications provide a diverse repertoire of epigenetic variation potentially for transgenerational inheritance. Elucidating and harnessing epigenetic inheritance will help us develop innovative breeding strategies and biotechnological tools to improve crop yield and resilience in the face of environmental challenges. Beyond plants, epigenetic principles are shared across sexually reproducing organisms including humans with relevance to medicine and public health.

Are cereal grasses a single genetic system?

In 1993, a passionate and provocative call to arms urged cereal researchers to consider the taxon they study as a single genetic system and collaborate with each other. Since then, that group of scientists has seen their discipline blossom. In an attempt to understand what unity of genetic systems means and how the notion was borne out by later research, we survey the progress and prospects of cereal genomics: sequence assemblies, population-scale sequencing, resistance gene cloning and domestication genetics. Gene order may not be as extraordinarily well conserved in the grasses as once thought. Still, several recurring themes have emerged. The same ancestral molecular pathways defining plant architecture have been co-opted in the evolution of different cereal crops. Such genetic convergence as much as cross-fertilization of ideas between cereal geneticists has led to a rich harvest of genes that, it is hoped, will lead to improved varieties.

Epigenetic control on transcription of vernalization genes and whole-genome gene expression profile induced by vernalization in common wheat

Vernalization is necessary for winter wheat to flower. However, it is unclear whether vernalization is also required for spring wheat, which is frequently sown in fall, and what molecular mechanisms underlie the vernalization response in wheat varieties. In this study, we examined the molecular mechanisms that regulate vernalization response in winter and spring wheat varieties. For this purpose, we determined how major vernalization genes (VRN1, VRN2, and VRN3) respond to vernalization in these varieties and whether modifications to histones play a role in changes in gene expression. We also identified genes that are differentially regulated in response to vernalization in winter and spring wheat varieties. We found that in winter wheat, but not in spring wheat, VRN1 expression decreases when returned to warm temperature following vernalization. This finding may be associated with differences between spring and winter wheat in the levels of tri-methylation of lysine 27 on histone H3 (H3K27me3) and tri-methylation of lysine 4 on histone H3 (H3K4me3) at the VRN1 gene. Analysis of winter wheat transcriptomes before and after vernalization revealed that vernalization influences the expression of several genes, including those involved in leucine catabolism, cysteine biosynthesis, and flavonoid biosynthesis. These findings provide new candidates for further study on the mechanism of vernalization regulation in wheat.

PlantFUNCO: Integrative Functional Genomics Database Reveals Clues into Duplicates Divergence Evolution

Evolutionary epigenomics and, more generally, evolutionary functional genomics, are emerging fields that study how non-DNA-encoded alterations in gene expression regulation are an important form of plasticity and adaptation. Previous evidence analyzing plants' comparative functional genomics has mostly focused on comparing same assay-matched experiments, missing the power of heterogeneous datasets for conservation inference. To fill this gap, we developed PlantFUN(ctional)CO(nservation) database, which is constituted by several tools and two main resources: interspecies chromatin states and functional genomics conservation scores, presented and analyzed in this work for three well-established plant models (Arabidopsis thaliana, Oryza sativa, and Zea mays). Overall, PlantFUNCO elucidated evolutionary information in terms of cross-species functional agreement. Therefore, providing a new complementary comparative-genomics source for assessing evolutionary studies. To illustrate the potential applications of this database, we replicated two previously published models predicting genetic redundancy in A. thaliana and found that chromatin states are a determinant of paralogs degree of functional divergence. These predictions were validated based on the phenotypes of mitochondrial alternative oxidase knockout mutants under two different stressors. Taking all the above into account, PlantFUNCO aim to leverage data diversity and extrapolate molecular mechanisms findings from different model organisms to determine the extent of functional conservation, thus, deepening our understanding of how plants epigenome and functional noncoding genome have evolved. PlantFUNCO is available at https://rocesv.github.io/PlantFUNCO.

Genome-wide profiling of histone (H3) lysine 4 (K4) tri-methylation (me3) under drought, heat, and combined stresses in switchgrass

BACKGROUND

Switchgrass (Panicum virgatum L.) is a warm-season perennial (C4) grass identified as an important biofuel crop in the United States. It is well adapted to the marginal environment where heat and moisture stresses predominantly affect crop growth. However, the underlying molecular mechanisms associated with heat and drought stress tolerance still need to be fully understood in switchgrass. The methylation of H3K4 is often associated with transcriptional activation of genes, including stress-responsive. Therefore, this study aimed to analyze genome-wide histone H3K4-tri-methylation in switchgrass under heat, drought, and combined stress.

RESULTS

In total, ~ 1.3 million H3K4me3 peaks were identified in this study using SICER. Among them, 7,342; 6,510; and 8,536 peaks responded under drought (DT), drought and heat (DTHT), and heat (HT) stresses, respectively. Most DT and DTHT peaks spanned 0 to + 2000 bases from the transcription start site [TSS]. By comparing differentially marked peaks with RNA-Seq data, we identified peaks associated with genes: 155 DT-responsive peaks with 118 DT-responsive genes, 121 DTHT-responsive peaks with 110 DTHT-responsive genes, and 175 HT-responsive peaks with 136 HT-responsive genes. We have identified various transcription factors involved in DT, DTHT, and HT stresses. Gene Ontology analysis using the AgriGO revealed that most genes belonged to biological processes. Most annotated peaks belonged to metabolite interconversion, RNA metabolism, transporter, protein modifying, defense/immunity, membrane traffic protein, transmembrane signal receptor, and transcriptional regulator protein families. Further, we identified significant peaks associated with TFs, hormones, signaling, fatty acid and carbohydrate metabolism, and secondary metabolites. qRT-PCR analysis revealed the relative expressions of six abiotic stress-responsive genes (transketolase, chromatin remodeling factor-CDH3, fatty-acid desaturase A, transmembrane protein 14C, beta-amylase 1, and integrase-type DNA binding protein genes) that were significantly (P < 0.05) marked during drought, heat, and combined stresses by comparing stress-induced against un-stressed and input controls.

CONCLUSION

Our study provides a comprehensive and reproducible epigenomic analysis of drought, heat, and combined stress responses in switchgrass. Significant enrichment of H3K4me3 peaks downstream of the TSS of protein-coding genes was observed. In addition, the cost-effective experimental design, modified ChIP-Seq approach, and analyses presented here can serve as a prototype for other non-model plant species for conducting stress studies.

Plant pangenomes for crop improvement, biodiversity and evolution

Plant genome sequences catalogue genes and the genetic elements that regulate their expression. Such inventories further research aims as diverse as mapping the molecular basis of trait diversity in domesticated plants or inquiries into the origin of evolutionary innovations in flowering plants millions of years ago. The transformative technological progress of DNA sequencing in the past two decades has enabled researchers to sequence ever more genomes with greater ease. Pangenomes - complete sequences of multiple individuals of a species or higher taxonomic unit - have now entered the geneticists' toolkit. The genomes of crop plants and their wild relatives are being studied with translational applications in breeding in mind. But pangenomes are applicable also in ecological and evolutionary studies, as they help classify and monitor biodiversity across the tree of life, deepen our understanding of how plant species diverged and show how plants adapt to changing environments or new selection pressures exerted by human beings.

Basic helix-loop-helix transcription factor OsbHLH110 positively regulates abscisic acid biosynthesis and salinity tolerance in rice

Salinity is a significant abiotic stress factor affecting plant growth, consequently reducing crop yield. Abscisic acid (ABA), a well-known phytohormone, is crucial in conferring resistance to abiotic stress, thus, understanding the mechanisms underlying ABA biosynthesis is crucial. In rice (Oryza sativa L.), OsABA2, a short-chain dehydrogenase protein, has a pivotal role in modulating ABA biosynthesis and salt tolerance by undergoing phosphorylation at Ser197 through mitogen-activated protein kinase OsMPK1. However, the interaction between OsABA2 and other proteins in regulating ABA biosynthesis remains unclear. We employed OsABA2 as a bait in yeast two-hybrid screening: a basic helix-loop-helix transcription factor interacting with OsABA2, named OsbHLH110, was identified. Our results showed that firefly luciferase complementary imaging, pull-down, and co-immunoprecipitation assays validated the interaction between OsbHLH110 and OsABA2, affirming their interaction in vivo and in vitro. Moreover, the expression of OsbHLH110 significantly increases in response to salt and ABA treatments. Additionally, OsbHLH110 can directly bind to the G-box element in the OsABA2 promoter. This binding enhances luciferase activity controlled by the OsABA2 promoter, thereby increasing the expression of the OsABA2 gene and content of the OsABA2 protein, resulting in an increase in ABA content. OsABA2 enhanced the interaction between OsbHLH110 and OsABA2 promoter. This collaborative effect enhanced the regulation of ABA biosynthesis. Subsequent genetic analysis demonstrated that OsbHLH110 improved the tolerance of rice to salt stress.

Application of Single-Cell Assay for Transposase-Accessible Chromatin with High Throughput Sequencing in Plant Science: Advances, Technical Challenges, and Prospects

The Single-cell Assay for Transposase-Accessible Chromatin with high throughput sequencing (scATAC-seq) has gained increasing popularity in recent years, allowing for chromatin accessibility to be deciphered and gene regulatory networks (GRNs) to be inferred at single-cell resolution. This cutting-edge technology now enables the genome-wide profiling of chromatin accessibility at the cellular level and the capturing of cell-type-specific cis-regulatory elements (CREs) that are masked by cellular heterogeneity in bulk assays. Additionally, it can also facilitate the identification of rare and new cell types based on differences in chromatin accessibility and the charting of cellular developmental trajectories within lineage-related cell clusters. Due to technical challenges and limitations, the data generated from scATAC-seq exhibit unique features, often characterized by high sparsity and noise, even within the same cell type. To address these challenges, various bioinformatic tools have been developed. Furthermore, the application of scATAC-seq in plant science is still in its infancy, with most research focusing on root tissues and model plant species. In this review, we provide an overview of recent progress in scATAC-seq and its application across various fields. We first conduct scATAC-seq in plant science. Next, we highlight the current challenges of scATAC-seq in plant science and major strategies for cell type annotation. Finally, we outline several future directions to exploit scATAC-seq technologies to address critical challenges in plant science, ranging from plant ENCODE(The Encyclopedia of DNA Elements) project construction to GRN inference, to deepen our understanding of the roles of CREs in plant biology.

Dissecting postharvest chilling injuries in pome and stone fruit through integrated omics

Lowering the storage temperature is an effective method to extend the postharvest and shelf life of fruits. Nevertheless, this technique often leads to physiological disorders, commonly known as chilling injuries. Apples and pears are susceptible to chilling injuries, among which superficial scald is the most economically relevant. Superficial scald is due to necrotic lesions of the first layers of hypodermis manifested through skin browning. In peaches and nectarines, chilling injuries are characterized by internal symptoms, such as mealiness. Fruits with these aesthetic or compositional/structural defects are not suitable for fresh consumption. Genetic variation is a key factor in determining fruit susceptibility to chilling injuries; however, physiological, or technical aspects such as harvest maturity and storage conditions also play a role. Multi-omics approaches have been used to provide an integrated explanation of chilling injury development. Metabolomics in pome fruits specifically targets the identification of ethylene, phenols, lipids, and oxidation products. Genomics and transcriptomics have revealed interesting connections with metabolomic datasets, pinpointing specific genes linked to cold stress, wax synthesis, farnesene metabolism, and the metabolic pathways of ascorbate and glutathione. When applied to Prunus species, these cutting-edge approaches have uncovered that the development of mealiness symptoms is linked to ethylene signaling, cell wall synthesis, lipid metabolism, cold stress genes, and increased DNA methylation levels. Emphasizing the findings from multi-omics studies, this review reports how the integration of omics datasets can provide new insights into understanding of chilling injury development. This new information is essential for successfully creating more resilient fruit varieties and developing novel postharvest strategies.

Epigenetic modifications regulate cultivar-specific root development and metabolic adaptation to nitrogen availability in wheat

The breeding of crops with improved nitrogen use efficiency (NUE) is crucial for sustainable agriculture, but the involvement of epigenetic modifications remains unexplored. Here, we analyze the chromatin landscapes of two wheat cultivars (KN9204 and J411) that differ in NUE under varied nitrogen conditions. The expression of nitrogen metabolism genes is closely linked to variation in histone modification instead of differences in DNA sequence. Epigenetic modifications exhibit clear cultivar-specificity, which likely contributes to distinct agronomic traits. Additionally, low nitrogen (LN) induces H3K27ac and H3K27me3 to significantly enhance root growth in KN9204, while remarkably inducing NRT2 in J411. Evidence from histone deacetylase inhibitor treatment and transgenic plants with loss function of H3K27me3 methyltransferase shows that changes in epigenetic modifications could alter the strategy preference for root development or nitrogen uptake in response to LN. Here, we show the importance of epigenetic regulation in mediating cultivar-specific adaptation to LN in wheat.

Genome resources for the elite bread wheat cultivar Aikang 58 and mining of elite homeologous haplotypes for accelerating wheat improvement

Despite recent progress in crop genomics studies, the genomic changes brought about by modern breeding selection are still poorly understood, thus hampering genomics-assisted breeding, especially in polyploid crops with compound genomes such as common wheat (Triticum aestivum). In this work, we constructed genome resources for the modern elite common wheat variety Aikang 58 (AK58). Comparative genomics between AK58 and the landrace cultivar Chinese Spring (CS) shed light on genomic changes that occurred through recent varietal improvement. We also explored subgenome diploidization and divergence in common wheat and developed a homoeologous locus-based genome-wide association study (HGWAS) approach, which was more effective than single homoeolog-based GWAS in unraveling agronomic trait-associated loci. A total of 123 major HGWAS loci were detected using a genetic population derived from AK58 and CS. Elite homoeologous haplotypes (HHs), formed by combinations of subgenomic homoeologs of the associated loci, were found in both parents and progeny, and many could substantially improve wheat yield and related traits. We built a website where users can download genome assembly sequence and annotation data for AK58, perform blast analysis, and run JBrowse. Our work enriches genome resources for wheat, provides new insights into genomic changes during modern wheat improvement, and suggests that efficient mining of elite HHs can make a substantial contribution to genomics-assisted breeding in common wheat and other polyploid crops.

Cotranscriptional demethylation induces global loss of H3K4me2 from active genes in Arabidopsis

Based on studies of animals and yeasts, methylation of histone H3 lysine 4 (H3K4me1/2/3, for mono-, di-, and tri-methylation, respectively) is regarded as the key epigenetic modification of transcriptionally active genes. In plants, however, H3K4me2 correlates negatively with transcription, and the regulatory mechanisms of this counterintuitive H3K4me2 distribution in plants remain largely unexplored. A previous genetic screen for factors regulating plant regeneration identified Arabidopsis LYSINE-SPECIFIC DEMETHYLASE 1-LIKE 3 (LDL3), which is a major H3K4me2 demethylase. Here, we show that LDL3-mediated H3K4me2 demethylation depends on the transcription elongation factor Paf1C and phosphorylation of the C-terminal domain (CTD) of RNA polymerase II (RNAPII). In addition, LDL3 binds to phosphorylated RNAPII. These results suggest that LDL3 is recruited to transcribed genes by binding to elongating RNAPII and demethylates H3K4me2 cotranscriptionally. Importantly, the negative correlation between H3K4me2 and transcription is significantly attenuated in the ldl3 mutant, demonstrating the genome-wide impacts of the transcription-driven LDL3 pathway to control H3K4me2 in plants. Our findings implicate H3K4me2 demethylation in plants as chromatin records of transcriptional activity, which ensures robust gene control.

Sea-ATI unravels novel vocabularies of plant active cistrome

The cistrome consists of all cis-acting regulatory elements recognized by transcription factors (TFs). However, only a portion of the cistrome is active for TF binding in a specific tissue. Resolving the active cistrome in plants remains challenging. In this study, we report the assay sequential extraction assisted-active TF identification (sea-ATI), a low-input method that profiles the DNA sequences recognized by TFs in a target tissue. We applied sea-ATI to seven plant tissues to survey their active cistrome and generated 41 motif models, including 15 new models that represent previously unidentified cis-regulatory vocabularies. ATAC-seq and RNA-seq analyses confirmed the functionality of the cis-elements from the new models, in that they are actively bound in vivo, located near the transcription start site, and influence chromatin accessibility and transcription. Furthermore, comparing dimeric WRKY CREs between sea-ATI and DAP-seq libraries revealed that thermodynamics and genetic drifts cooperatively shaped their evolution. Notably, sea-ATI can identify not only positive but also negative regulatory cis-elements, thereby providing unique insights into the functional non-coding genome of plants.

LHP1-mediated epigenetic buffering of subgenome diversity and defense responses confers genome plasticity and adaptability in allopolyploid wheat

Polyploidization is a major driver of genome diversification and environmental adaptation. However, the merger of different genomes may result in genomic conflicts, raising a major question regarding how genetic diversity is interpreted and regulated to enable environmental plasticity. By analyzing the genome-wide binding of 191 trans-factors in allopolyploid wheat, we identified like heterochromatin protein 1 (LHP1) as a master regulator of subgenome-diversified genes. Transcriptomic and epigenomic analyses of LHP1 mutants reveal its role in buffering the expression of subgenome-diversified defense genes by controlling H3K27me3 homeostasis. Stripe rust infection releases latent subgenomic variations by eliminating H3K27me3-related repression. The simultaneous inactivation of LHP1 homoeologs by CRISPR-Cas9 confers robust stripe rust resistance in wheat seedlings. The conditional repression of subgenome-diversified defenses ensures developmental plasticity to external changes, while also promoting neutral-to-non-neutral selection transitions and adaptive evolution. These findings establish an LHP1-mediated buffering system at the intersection of genotypes, environments, and phenotypes in polyploid wheat. Manipulating the epigenetic buffering capacity offers a tool to harness cryptic subgenomic variations for crop improvement.

Replacement of Arabidopsis H2A.Z with human H2A.Z orthologs reveals extensive functional conservation and limited importance of the N-terminal tail sequence for Arabidopsis development

The incorporation of histone variants, distinct paralogs of core histones, into chromatin affects all DNA-templated processes in the cell, including the regulation of transcription. In recent years, much research has been focused on H2A.Z, an evolutionarily conserved H2A variant found in all eukaryotes. In order to investigate the functional conservation of H2A.Z histones during eukaryotic evolution we transformed h2a.z deficient plants with three human H2A.Z proteins to assess their ability to rescue the mutant defects. We discovered that human H2A.Z.1 and H2A.Z.2.1 fully complement the phenotypic abnormalities of h2a.z plants despite the fact that Arabidopsis and human H2A.Z N-terminal tail sequences are quite divergent. In contrast, the brain-specific splice variant H2A.Z.2.2 has a dominant-negative effect in wild-type plants. Furthermore, H2A.Z.1 almost completely re-establishes normal H2A.Z chromatin occupancy in h2a.z plants and restores the transcript levels of more than 84 % of misexpressed genes. Finally, our hypothesis that the N-terminal tail of Arabidopsis H2A.Z is not crucial for its developmental functions was supported by the ability of N-terminal end truncations of Arabidopsis HTA11 to largely rescue the defects of h2a.z mutants.

Epigenomic divergence correlates with sequence polymorphism in Arabidopsis paralogs

Processes affecting rates of sequence polymorphism are fundamental to the evolution of gene duplicates. The relationship between gene activity and sequence polymorphism can influence the likelihood that functionally redundant gene copies are co-maintained in stable evolutionary equilibria vs other outcomes such as neofunctionalization. Here, we investigate genic variation in epigenome-associated polymorphism rates in Arabidopsis thaliana and consider whether these affect the evolution of gene duplicates. We compared the frequency of sequence polymorphism and patterns of genetic differentiation between genes classified by exon methylation patterns: unmethylated (unM), gene-body methylated (gbM), and transposon-like methylated (teM) states, which reflect divergence in gene expression. We found that the frequency of polymorphism was higher in teM (transcriptionally repressed, tissue-specific) genes and lower in gbM (active, constitutively expressed) genes. Comparisons of gene duplicates were largely consistent with genome-wide patterns - gene copies that exhibit teM accumulate more variation, evolve faster, and are in chromatin states associated with reduced DNA repair. This relationship between expression, the epigenome, and polymorphism may lead to the breakdown of equilibrium states that would otherwise maintain genetic redundancies. Epigenome-mediated polymorphism rate variation may facilitate the evolution of novel gene functions in duplicate paralogs maintained over evolutionary time.

Domains Rearranged Methylase 2 maintains DNA methylation at large DNA hypomethylated shores and long-range chromatin interactions in rice

DNA methylation plays an important role in gene regulation and genomic stability. However, large DNA hypomethylated regions known as DNA methylation valleys (DMVs) or canyons have also been suggested to serve unique regulatory functions, largely unknown in rice (Oryza sativa). Here, we describe the DMVs in rice seedlings, which were highly enriched with developmental and transcription regulatory genes. Further detailed analysis indicated that grand DMVs (gDMVs) might be derived from nuclear integrants of organelle DNA (NORGs). Furthermore, Domains Rearranged Methylase 2 (OsDRM2) maintained DNA methylation at short DMV (sDMV) shores. Epigenetic maps indicated that sDMVs were marked with H3K4me3 and/or H3K27me3, although the loss of DNA methylation had a negligible effect on histone modification within these regions. In addition, we constructed H3K27me3-associated interaction maps for homozygous T-DNA insertion mutant of the gene (osdrm2) and wild type (WT). From a global perspective, most (90%) compartments were stable between osdrm2 and WT plants. At a high resolution, we observed a dramatic loss of long-range chromatin loops in osdrm2, which suffered an extensive loss of non-CG (CHG and CHH, H = A, T, or C) methylation. From another viewpoint, the loss of non-CG methylation at sDMV shores in osdrm2 could disrupt H3K27me3-mediated chromatin interaction networks. Overall, our results demonstrated that DMVs are a key genomic feature in rice and are precisely regulated by epigenetic modifications, including DNA methylation and histone modifications. OsDRM2 maintained DNA methylation at sDMV shores, while OsDRM2 deficiency strongly affected three-dimensional (3D) genome architectures.

Paternal DNA methylation is remodeled to maternal levels in rice zygote

Epigenetic reprogramming occurs during reproduction to reset the genome for early development. In flowering plants, mechanistic details of parental methylation remodeling in zygote remain elusive. Here we analyze allele-specific DNA methylation in rice hybrid zygotes and during early embryo development and show that paternal DNA methylation is predominantly remodeled to match maternal allelic levels upon fertilization, which persists after the first zygotic division. The DNA methylation remodeling pattern supports the predominantly maternal-biased gene expression during zygotic genome activation (ZGA) in rice. However, parental allelic-specific methylations are reestablished at the globular embryo stage and associate with allelic-specific histone modification patterns in hybrids. These results reveal that paternal DNA methylation is remodeled to match the maternal pattern during zygotic genome reprogramming and suggest existence of a chromatin memory allowing parental allelic-specific methylation to be maintained in the hybrid.

Integrated 3D genome, epigenome and transcriptome analyses reveal transcriptional coordination of circadian rhythm in rice

Photoperiods integrate with the circadian clock to coordinate gene expression rhythms and thus ensure plant fitness to the environment. Genome-wide characterization and comparison of rhythmic genes under different light conditions revealed delayed phase under constant darkness (DD) and reduced amplitude under constant light (LL) in rice. Interestingly, ChIP-seq and RNA-seq profiling of rhythmic genes exhibit synchronous circadian oscillation in H3K9ac modifications at their loci and long non-coding RNAs (lncRNAs) expression at proximal loci. To investigate how gene expression rhythm is regulated in rice, we profiled the open chromatin regions and transcription factor (TF) footprints by time-series ATAC-seq. Although open chromatin regions did not show circadian change, a significant number of TFs were identified to rhythmically associate with chromatin and drive gene expression in a time-dependent manner. Further transcriptional regulatory networks mapping uncovered significant correlation between core clock genes and transcription factors involved in light/temperature signaling. In situ Hi-C of ZT8-specific expressed genes displayed highly connected chromatin association at the same time, whereas this ZT8 chromatin connection network dissociates at ZT20, suggesting the circadian control of gene expression by dynamic spatial chromatin conformation. These findings together implicate the existence of a synchronization mechanism between circadian H3K9ac modifications, chromatin association of TF and gene expression, and provides insights into circadian dynamics of spatial chromatin conformation that associate with gene expression rhythms.

Asymmetric variation in DNA methylation during domestication and de-domestication of rice

Hundreds of plant species have been domesticated to feed human civilization, while some crops have undergone de-domestication into agricultural weeds, threatening global food security. To understand the genetic and epigenetic basis of crop domestication and de-domestication, we generated DNA methylomes from 95 accessions of wild rice (Oryza rufipogon L.), cultivated rice (Oryza sativa L.) and weedy rice (O. sativa f. spontanea). We detected a significant decrease in DNA methylation over the course of rice domestication but observed an unexpected increase in DNA methylation through de-domestication. Notably, DNA methylation changes occurred in distinct genomic regions for these 2 opposite stages. Variation in DNA methylation altered the expression of nearby and distal genes through affecting chromatin accessibility, histone modifications, transcription factor binding, and the formation of chromatin loops, which may contribute to morphological changes during domestication and de-domestication of rice. These insights into population epigenomics underlying rice domestication and de-domestication provide resources and tools for epigenetic breeding and sustainable agriculture.

A DNA adenine demethylase impairs PRC2-mediated repression of genes marked by a specific chromatin signature

BACKGROUND

The Fe (II)- and α-ketoglutarate-dependent AlkB family dioxygenases are implicated in nucleotide demethylation. AlkB homolog1 (ALKBH1) is shown to demethylate DNA adenine methylation (6mA) preferentially from single-stranded or unpaired DNA, while its demethylase activity and function in the chromatin context are unclear.

RESULTS

Here, we find that loss-of-function of the rice ALKBH1 gene leads to increased 6mA in the R-loop regions of the genome but has a limited effect on the overall 6mA level. However, in the context of mixed tissues, rather than on individual loci, the ALKBH1 mutation or overexpression mainly affects the expression of genes with a specific combination of chromatin modifications in the body region marked with H3K4me3 and H3K27me3 but depleted of DNA CG methylation. In the similar context of mixed tissues, further analysis reveals that the ALKBH1 protein preferentially binds to genes marked by the chromatin signature and has a function to maintain a high H3K4me3/H3K27me3 ratio by impairing the binding of Polycomb repressive complex 2 (PRC2) to the targets, which is required for both the basal and stress-induced expression of the genes.

CONCLUSION

Our findings unravel a function of ALKBH1 to control the balance between the antagonistic histone methylations for gene activity and provide insight into the regulatory mechanism of PRC2-mediated H3K27me3 deposition within the gene body region.

QTG-Miner aids rapid dissection of the genetic base of tassel branch number in maize

Genetic dissection of agronomic traits is important for crop improvement and global food security. Phenotypic variation of tassel branch number (TBN), a major breeding target, is controlled by many quantitative trait loci (QTLs). The lack of large-scale QTL cloning methodology constrains the systematic dissection of TBN, which hinders modern maize breeding. Here, we devise QTG-Miner, a multi-omics data-based technique for large-scale and rapid cloning of quantitative trait genes (QTGs) in maize. Using QTG-Miner, we clone and verify seven genes underlying seven TBN QTLs. Compared to conventional methods, QTG-Miner performs well for both major- and minor-effect TBN QTLs. Selection analysis indicates that a substantial number of genes and network modules have been subjected to selection during maize improvement. Selection signatures are significantly enriched in multiple biological pathways between female heterotic groups and male heterotic groups. In summary, QTG-Miner provides a large-scale approach for rapid cloning of QTGs in crops and dissects the genetic base of TBN for further maize breeding.

3D organization of regulatory elements for transcriptional regulation in Arabidopsis

BACKGROUND

Although spatial organization of compartments and topologically associating domains at large scale is relatively well studied, the spatial organization of regulatory elements at fine scale is poorly understood in plants.

RESULTS

Here we perform high-resolution chromatin interaction analysis using paired-end tag sequencing approach. We map chromatin interactions tethered with RNA polymerase II and associated with heterochromatic, transcriptionally active, and Polycomb-repressive histone modifications in Arabidopsis. Analysis of the regulatory repertoire shows that distal active cis-regulatory elements are linked to their target genes through long-range chromatin interactions with increased expression of the target genes, while poised cis-regulatory elements are linked to their target genes through long-range chromatin interactions with depressed expression of the target genes. Furthermore, we demonstrate that transcription factor MYC2 is critical for chromatin spatial organization, and propose that MYC2 occupancy and MYC2-mediated chromatin interactions coordinately facilitate transcription within the framework of 3D chromatin architecture. Analysis of functionally related gene-defined chromatin connectivity networks reveals that genes implicated in flowering-time control are functionally compartmentalized into separate subdomains via their spatial activity in the leaf or shoot apical meristem, linking active mark- or Polycomb-repressive mark-associated chromatin conformation to coordinated gene expression.

CONCLUSION

The results reveal that the regulation of gene transcription in Arabidopsis is not only by linear juxtaposition, but also by long-range chromatin interactions. Our study uncovers the fine scale genome organization of Arabidopsis and the potential roles of such organization in orchestrating transcription and development.

Suppressing a phosphohydrolase of cytokinin nucleotide enhances grain yield in rice

One-step and two-step pathways are proposed to synthesize cytokinin in plants. The one-step pathway is mediated by LONELY GUY (LOG) proteins. However, the enzyme for the two-step pathway remains to be identified. Here, we show that quantitative trait locus GY3 may boost grain yield by more than 20% through manipulating a two-step pathway. Locus GY3 encodes a LOG protein that acts as a 5'-ribonucleotide phosphohydrolase by excessively consuming the cytokinin precursors, which contrasts with the activity of canonical LOG members as phosphoribohydrolases in a one-step pathway. The residue S41 of GY3 is crucial for the dephosphorylation of iPRMP to produce iPR. A solo-LTR insertion within the promoter of GY3 suppressed its expression and resulted in a higher content of active cytokinins in young panicles. Introgression of GY302428 increased grain yield per plot by 7.4% to 16.3% in all investigated indica backgrounds, which demonstrates the great value of GY302428 in indica rice production.

Histone Maps in Gossypium darwinii Reveal Epigenetic Regulation Drives Subgenome Divergence and Cotton Domestication

The functional annotation of genomes, including chromatin modifications, is essential to understand the intricate architecture of chromatin and the consequential gene regulation. However, such an annotation remains limited for cotton genomes. Here, we conducted chromatin profiling in a wild allotetraploid cotton Gossypium darwinii (AD genome) by integrating the data of histone modification, transcriptome, and chromatin accessibility. We revealed that the A subgenome showed a higher level of active histone marks and lower level of repressive histone marks than the D subgenome, which was consistent with the expression bias between the two subgenomes. We show that the bias in transcription and histone modification between the A and D subgenomes may be caused by genes unique to the subgenome but not by homoeologous genes. Moreover, we integrate histone marks and open chromatin to define six chromatin states (S1-S6) across the cotton genome, which index different genomic elements including genes, promoters, and transposons, implying distinct biological functions. In comparison to the domesticated cotton species, we observed that 23.2% of genes in the genome exhibit a transition from one chromatin state to another at their promoter. Strikingly, the S2 (devoid of epigenetic marks) to S3 (enriched for the mark of open chromatin) was the largest transition group. These transitions occurred simultaneously with changes in gene expression, which were significantly associated with several domesticated traits in cotton. Collectively, our study provides a useful epigenetic resource for research on allopolyploid plants. The domestication-induced chromatin dynamics and associated genes identified here will aid epigenetic engineering, improving polyploid crops.

The multi-BRCT domain protein DDRM2 promotes the recruitment of RAD51 to DNA damage sites to facilitate homologous recombination

DNA double-strand breaks (DSBs) are the most toxic form of DNA damage in cells. Homologous recombination (HR) is an error-free repair mechanism for DSBs as well as a basis for gene targeting using genome-editing techniques. Despite the importance of HR, the HR mechanism in plants is poorly understood. Through genetic screens for DNA damage response mutants (DDRMs), we find that the Arabidopsis ddrm2 mutant is hypersensitive to DSB-inducing reagents. DDRM2 encodes a protein with four BRCA1 C-terminal (BRCT) domains and is highly conserved in plants including bryophytes, the earliest land plant lineage. The plant-specific transcription factor SOG1 binds to the promoter of DDRM2 and activates its expression. In consistence, the expression of DDRM2 is induced by DSBs in a SOG1-dependent manner. In support, genetic analysis suggests that DDRM2 functions downstream of SOG1. Similar to the sog1 mutant, the ddrm2 mutant shows dramatically reduced HR efficiency. Mechanistically, DDRM2 interacts with the core HR protein RAD51 and is required for the recruitment of RAD51 to DSB sites. Our study reveals that SOG1-DDRM2-RAD51 is a novel module for HR, providing a potential target for improving the efficiency of gene targeting.

Haplotype mapping of H3K27me3-associated chromatin interactions defines topological regulation of gene silencing in rice

Histone modification H3K27me3 is an important chromatin mark that plays vital roles in repressing expression of developmental genes. Here, we construct high-resolution 3D genome maps using long-read chromatin interaction analysis by paired-end tag sequencing (ChIA-PET) and characterize H3K27me3-associated chromatin interactions in an elite rice hybrid, Shanyou 63. We find that many H3K27me3-marked regions may function as silencer-like regulatory elements. The silencer-like elements can come into proximity with distal target genes via forming chromatin loops in 3D space of the nuclei, regulating gene silencing and plant traits. Natural and induced deletion of silencers upregulate expression of distal connected genes. Furthermore, we identify extensive allele-specific chromatin loops. We find that genetic variations alter allelic chromatin topology, thus modulating allelic gene imprinting in rice hybrids. In conclusion, the characterization of silencer-like regulatory elements and haplotype-resolved chromatin interaction maps provide insights into the understanding of molecular mechanisms underlying allelic gene silencing and plant trait controlling.

A VEL3 histone deacetylase complex establishes a maternal epigenetic state controlling progeny seed dormancy

Mother plants play an important role in the control of dormancy and dispersal characters of their progeny. In Arabidopsis seed dormancy is imposed by the embryo-surrounding tissues of the endosperm and seed coat. Here we show that VERNALIZATION5/VIN3-LIKE 3 (VEL3) maintains maternal control over progeny seed dormancy by establishing an epigenetic state in the central cell that primes the depth of primary seed dormancy later established during seed maturation. VEL3 colocalises with MSI1 in the nucleolus and associates with a histone deacetylase complex. Furthermore, VEL3 preferentially associates with pericentromeric chromatin and is required for deacetylation and H3K27me3 deposition established in the central cell. The epigenetic state established by maternal VEL3 is retained in mature seeds, and controls seed dormancy in part through repression of programmed cell death-associated gene ORE1. Our data demonstrates a mechanism by which maternal control of progeny seed physiology persists post-shedding, maintaining parental control of seed behaviour.

The SMC5/6 complex recruits the PAF1 complex to facilitate DNA double-strand break repair in Arabidopsis

DNA double-strand breaks (DSBs) are one of the most toxic forms of DNA damage, which threatens genome stability. Homologous recombination is an error-free DSB repair pathway, in which the evolutionarily conserved SMC5/6 complex (SMC5/6) plays essential roles. The PAF1 complex (PAF1C) is well known to regulate transcription. Here we show that SMC5/6 recruits PAF1C to facilitate DSB repair in plants. In a genetic screen for DNA damage response mutants (DDRMs), we found that the Arabidopsis ddrm4 mutant is hypersensitive to DSB-inducing agents and is defective in homologous recombination. DDRM4 encodes PAF1, a core subunit of PAF1C. Further biochemical and genetic studies reveal that SMC5/6 recruits PAF1C to DSB sites, where PAF1C further recruits the E2 ubiquitin-conjugating enzymes UBC1/2, which interact with the E3 ubiquitin ligases HUB1/2 to mediate the monoubiquitination of histone H2B at DSBs. These results implicate SMC5/6-PAF1C-UBC1/2-HUB1/2 as a new axis for DSB repair through homologous recombination, revealing a new mechanism of SMC5/6 and uncovering a novel function of PAF1C.

Pan-genome inversion index reveals evolutionary insights into the subpopulation structure of Asian rice

Understanding and exploiting genetic diversity is a key factor for the productive and stable production of rice. Here, we utilize 73 high-quality genomes that encompass the subpopulation structure of Asian rice (Oryza sativa), plus the genomes of two wild relatives (O. rufipogon and O. punctata), to build a pan-genome inversion index of 1769 non-redundant inversions that span an average of ~29% of the O. sativa cv. Nipponbare reference genome sequence. Using this index, we estimate an inversion rate of ~700 inversions per million years in Asian rice, which is 16 to 50 times higher than previously estimated for plants. Detailed analyses of these inversions show evidence of their effects on gene expression, recombination rate, and linkage disequilibrium. Our study uncovers the prevalence and scale of large inversions (≥100 bp) across the pan-genome of Asian rice and hints at their largely unexplored role in functional biology and crop performance.

Approaches for Modes of Action Study of Long Non-Coding RNAs: From Single Verification to Genome-Wide Determination

Long non-coding RNAs (lncRNAs) are transcripts longer than 200 nucleotides (nt) that are not translated into known functional proteins. This broad definition covers a large collection of transcripts with diverse genomic origins, biogenesis, and modes of action. Thus, it is very important to choose appropriate research methodologies when investigating lncRNAs with biological significance. Multiple reviews to date have summarized the mechanisms of lncRNA biogenesis, their localization, their functions in gene regulation at multiple levels, and also their potential applications. However, little has been reviewed on the leading strategies for lncRNA research. Here, we generalize a basic and systemic mind map for lncRNA research and discuss the mechanisms and the application scenarios of ‘up-to-date’ techniques as applied to molecular function studies of lncRNAs. Taking advantage of documented lncRNA research paradigms as examples, we aim to provide an overview of the developing techniques for elucidating lncRNA interactions with genomic DNA, proteins, and other RNAs. In the end, we propose the future direction and potential technological challenges of lncRNA studies, focusing on techniques and applications.

Genome-wide analysis of transcriptome and histone modifications in Brassica napus hybrid

Although utilization of heterosis has largely improved the yield of many crops worldwide, the underlying molecular mechanism of heterosis, particularly for allopolyploids, remains unclear. Here, we compared epigenome and transcriptome data of an elite hybrid and its parental lines in three assessed tissues (seedling, flower bud, and silique) to explore their contribution to heterosis in allopolyploid B. napus. Transcriptome analysis illustrated that a small proportion of non-additive genes in the hybrid compared with its parents, as well as parental expression level dominance, might have a significant effect on heterosis. We identified histone modification (H3K4me3 and H3K27me3) variation between the parents and hybrid, most of which resulted from the differences between parents. H3K4me3 variations were positively correlated with gene expression differences among the hybrid and its parents. Furthermore, H3K4me3 and H3K27me3 were rather stable in hybridization and were mainly inherited additively in the B. napus hybrid. Together, our data revealed that transcriptome reprogramming and histone modification remodeling in the hybrid could serve as valuable resources for better understanding heterosis in allopolyploid crops.

Natural variation and artificial selection at the BnaC2.MYB28 locus modulate Brassica napus seed glucosinolate

The degradation products of glucosinolates (GSLs) greatly lower the nutritional value of rapeseed (Brassica napus) meal; thus, reduction of seed GSL content (SGC) has become an important objective of rapeseed breeding. In our previous study, we finely mapped a major QTL (qGSL-C2) for SGC to a 49-kb collinear region on B. rapa chromosome A2. Here, we experimentally validated that BnaC2.MYB28, encoding an R2R3-MYB transcription factor, is the causal gene of qGSL-C2. BnaC2.MYB28 is a nucleus-localized protein mainly expressed in vegetative tissues. Knockout of BnaC2.MYB28 in the high-SGC parent G120 reduced SGC to a value lower than that in the low-SGC parent ZY50, while overexpression of BnaC2.MYB28 in both parental lines (G120 and ZY50) led to extremely high SGC, indicating that BnaC2.MYB28 acts as a positive regulator of SGC in both parents. Molecular characterization revealed that BnaC2.MYB28 forms a homodimer and specifically interacts with BnaMYC3. Moreover, BnaC2.MYB28 can directly activate the expression of GSL biosynthesis genes. Differential expression abundance resulting from the polymorphic promoter sequences, in combination with the different capability in activating downstream genes involved in aliphatic GSL biosynthesis, caused the functional divergence of BnaC2.MYB28 in SGC regulation between the parents. Natural variation of BnaC2.MYB28 was highly associated with SGC in natural germplasm and has undergone artificial selection in modern low-GSL breeding. This study provides important insights into the core function of BnaC2.MYB28 in regulating SGC and a promising strategy for manipulating SGC in rapeseed.

Plant pan-genomics and its applications

Plant genomes are so highly diverse that a substantial proportion of genomic sequences are not shared among individuals. The variable DNA sequences, along with the conserved core sequences, compose the more sophisticated pan-genome that represents the collection of all non-redundant DNA in a species. With rapid progress in genome sequencing technologies, pan-genome research in plants is now accelerating. Here we review recent advances in plant pan-genomics, including major driving forces of structural variations that constitute the variable sequences, methodological innovations for representing the pan-genome, and major successes in constructing plant pan-genomes. We also summarize recent efforts toward decoding the remaining dark matter in telomere-to-telomere or gapless plant genomes. These new genome resources, which have remarkable advantages over numerous previously assembled less-than-perfect genomes, are expected to become new references for genetic studies and plant breeding.

CUT&Tag for Mapping In Vivo Protein-DNA Interactions in Plants

Histone post-translational modifications and transcription factors (TFs) play vital roles in regulating gene expression. A comprehensive understanding of transcriptional regulation requires genome-wide mapping of chromatin features such as histone modifications and TF binding sites. Here, we describe a detailed nucleus CUT&Tag (Cleavage Under Targets and Tagmentation) protocol, which is an antibody-guided in situ protein-DNA interaction mapping method using protein A/G fused Tn5 transposase. Compared with regular ChIP-seq in plants, nucleus CUT&Tag (nCUT&Tag) omits many steps such as sonication and immunoprecipitation, thus saving much time and making it possible to efficiently profile chromatin features from low-input and even single cells with higher signal-to-noise ratio.

Chromatin accessibility illuminates single-cell regulatory dynamics of rice root tips

BACKGROUND

Root development and function have central roles in plant adaptation to the environment. The modification of root traits has additionally been a major driver of crop performance since the green revolution; however, the molecular underpinnings and the regulatory programmes defining root development and response to environmental stress remain largely unknown. Single-cell reconstruction of gene regulatory programmes provides an important tool to understand the cellular phenotypic variation in complex tissues and their response to endogenous and environmental stimuli. While single-cell transcriptomes of several plant organs have been elucidated, the underlying chromatin landscapes associated with cell type-specific gene expression remain largely unexplored.

RESULTS

To comprehensively delineate chromatin accessibility during root development of an important crop, we applied single-cell ATAC-seq (scATAC-seq) to 46,758 cells from rice root tips under normal and heat stress conditions. Our data revealed cell type-specific accessibility variance across most of the major cell types and allowed us to identify sets of transcription factors which associate with accessible chromatin regions (ACRs). Using root hair differentiation as a model, we demonstrate that chromatin and gene expression dynamics during cell type differentiation correlate in pseudotime analyses. In addition to developmental trajectories, we describe chromatin responses to heat and identify cell type-specific accessibility changes to this key environmental stimulus.

CONCLUSIONS

We report chromatin landscapes during rice root development at single-cell resolution. Our work provides a framework for the integrative analysis of regulatory dynamics in this important crop organ at single-cell resolution.

Disease resistance conferred by components of essential chrysanthemum oil and the epigenetic regulation of OsTPS1

The sesquiterpene alpha-bisabolol is the predominant active ingredient in essential oils that are highly valued in the cosmetics industry due to its wound healing, anti-inflammatory, and skin-soothing properties. Alpha-bisabolol was thought to be restricted to Compositae plants. Here we reveal that alpha-bisabolol is also synthesized in rice, a non-Compositae plant, where it acts as a novel sesquiterpene phytoalexin. Overexpressing the gene responsible for the biosynthesis of alpha-bisabolol, OsTPS1, conferred bacterial blight resistance in rice. Phylogenomic analyses revealed that alpha-bisabolol-synthesizing enzymes in rice and Compositae evolved independently. Further experiments demonstrated that the natural variation in the disease resistance level was associated with differential transcription of OsTPS1 due to polymorphisms in its promoter. We demonstrated that OsTPS1 was regulated at the epigenetic level by JMJ705 through the methyl jasmonate pathway. These data reveal the cross-family accumulation and regulatory mechanisms of alpha-bisabolol production.

Transcription elongator SPT6L regulates the occupancies of the SWI2/SNF2 chromatin remodelers SYD/BRM and nucleosomes at transcription start sites in Arabidopsis

Chromatin remodelers have been thought to be crucial in creating an accessible chromatin environment before transcription activation. However, it is still unclear how chromatin remodelers recognize and bind to the active regions. In this study, we found that chromatin remodelers SPLAYED (SYD) and BRAHMA (BRM) interact and co-occupy with Suppressor of Ty6-like (SPT6L), a core subunit of the transcription machinery, at thousands of the transcription start sites (TSS). The association of SYD and BRM to chromatin is dramatically reduced in spt6l and can be restored mainly by SPT6LΔtSH2, which binds to TSS in a RNA polymerase II (Pol II)-independent manner. Furthermore, SPT6L and SYD/BRM are involved in regulating the nucleosome and Pol II occupancy around TSS. The presence of SPT6L is sufficient to restore the association of the chromatin remodeler SYD to chromatin and maintain normal nucleosome occupancy. Our findings suggest that the two chromatin remodelers can form protein complexes with the core subunit of the transcription machinery and regulate nucleosome occupancy in the early transcription stage.

Systematic annotation of conservation states provides insights into regulatory regions in rice

Plant genomes contain a large fraction of noncoding sequences. The discovery and annotation of conserved noncoding sequences (CNSs) in plants is an ongoing challenge. Here we report the application of comparative genomics to systematically identify CNSs in 50 well-annotated Gramineae genomes using rice (Oryza sativa) as the reference. We conduct multiple-way whole-genome alignments to the rice genome. The rice genome is annotated as 20 conservation states (CSs) at single-nucleotide resolution using a multivariate hidden Markov model (ConsHMM) based on the multiple-genome alignments. Different states show distinct enrichments for various genomic features, and the conservation scores of CSs are highly correlated with the level of associated chromatin accessibility. We find that at least 33.5% of the rice genome is highly under selection, with more than 70% of the sequence lying outside of coding regions. A catalog of 855,366 regulatory CNSs is generated, and they significantly overlapped with putative active regulatory elements such as promoters, enhancers, and transcription factor binding sites. Collectively, our study provides a resource for elucidating functional noncoding regions of the rice genome and an evolutionary aspect of regulatory sequences in higher plants.

Genome-wide prediction of activating regulatory elements in rice by combining STARR-seq with FACS

Self-transcribing active regulatory region sequencing (STARR-seq) is widely used to identify enhancers at the whole-genome level. However, whether STARR-seq works as efficiently in plants as in animal systems remains unclear. Here, we determined that the traditional STARR-seq method can be directly applied to rice (Oryza sativa) protoplasts to identify enhancers, though with limited efficiency. Intriguingly, we identified not only enhancers but also constitutive promoters with this technique. To increase the performance of STARR-seq in plants, we optimized two procedures. We coupled fluorescence activating cell sorting (FACS) with STARR-seq to alleviate the effect of background noise, and we minimized PCR cycles and retained duplicates during prediction, which significantly increased the positive rate for activating regulatory elements (AREs). Using this method, we determined that AREs are associated with AT-rich regions and are enriched for a motif that the AP2/ERF family can recognize. Based on GC content preferences, AREs are clustered into two groups corresponding to promoters and enhancers. Either AT- or GC-rich regions within AREs could boost transcription. Additionally, disruption of AREs resulted in abnormal expression of both proximal and distal genes, which suggests that STARR-seq-revealed elements function as enhancers in vivo. In summary, our work provides a promising method to identify AREs in plants.

Plant synthetic epigenomic engineering for crop improvement

Efforts have been directed to redesign crops with increased yield, stress adaptability, and nutritional value through synthetic biology-the application of engineering principles to biology. A recent expansion in our understanding of how epigenetic mechanisms regulate plant development and stress responses has unveiled a new set of resources that can be harnessed to develop improved crops, thus heralding the promise of "synthetic epigenetics." In this review, we summarize the latest advances in epigenetic regulation and highlight how innovative sequencing techniques, epigenetic editing, and deep learning-driven predictive tools can rapidly extend these insights. We also proposed the future directions of synthetic epigenetics for the development of engineered smart crops that can actively monitor and respond to internal and external cues throughout their life cycles.

Integration of multi-omics technologies for crop improvement: Status and prospects

With the rapid development of next-generation sequencing (NGS), multi-omics techniques have been emerging as effective approaches for crop improvement. Here, we focus mainly on addressing the current status and future perspectives toward omics-related technologies and bioinformatic resources with potential applications in crop breeding. Using a large amount of omics-level data from the functional genome, transcriptome, proteome, epigenome, metabolome, and microbiome, clarifying the interaction between gene and phenotype formation will become possible. The integration of multi-omics datasets with pan-omics platforms and systems biology could predict the complex traits of crops and elucidate the regulatory networks for genetic improvement. Different scales of trait predictions and decision-making models will facilitate crop breeding more intelligent. Potential challenges that integrate the multi-omics data with studies of gene function and their network to efficiently select desirable agronomic traits are discussed by proposing some cutting-edge breeding strategies for crop improvement. Multi-omics-integrated approaches together with other artificial intelligence techniques will contribute to broadening and deepening our knowledge of crop precision breeding, resulting in speeding up the breeding process.

DNA methylation underpins the epigenomic landscape regulating genome transcription in Arabidopsis

BACKGROUND

It is challenging to determine the effect of DNA methylation on the epigenetic landscape and the function in higher organisms due to the lack of DNA methylation-free mutants.

RESULTS

Here, the analysis of a recently generated Arabidopsis mutant completely devoid of DNA methylation reveals that DNA methylation underpins the genome-wide landscape of histone modifications. Complete loss of DNA methylation causes an upheaval of the histone modification landscape, including complete loss of H3K9me2 and widespread redistribution of active and H3K27me3 histone marks, mostly owing to the role of DNA methylation in initiating H3K9me2 deposition and excluding active marks and repressive mark H3K27me3; CG and non-CG methylation can act independently at some genomic regions while they act cooperatively at many other regions. The transcriptional reprogramming upon loss of all DNA methylation correlates with the extensive redistribution or switches of the examined histone modifications. Histone modifications retained or gained in the DNA methylation-free mutant serve as DNA methylation-independent transcriptional regulatory signals: active marks promote genome transcription, whereas the repressive mark H3K27me3 compensates for the lack of DNA hypermethylation/H3K9me2 at multiple transposon families.

CONCLUSIONS

Our results show that an intact DNA methylome constitutes the scaffolding of the epigenomic landscape in Arabidopsis and is critical for controlled genome transcription and ultimately for proper growth and development.

Genome-wide binding analysis of transcription factor Rice Indeterminate 1 reveals a complex network controlling rice floral transition

RICE INDETERMINATE 1 (RID1) plays a critical role in controlling floral transition in rice (Oryza sativa). However, the molecular basis for this effect, particularly the target genes and regulatory specificity, remains largely unclear. Here, we performed chromatin immunoprecipitation followed by sequencing (ChIP-seq) in young leaves at the pre-floral-transition stage to identify the target genes of RID1, identifying 2,680 genes associated with RID1 binding sites genome-wide. RID1 binding peaks were highly enriched for TTTGTC, the direct binding motif of the INDETERMINATE DOMAIN protein family that includes RID1. Interestingly, CACGTG and GTGGGCCC, two previously uncharacterized indirect binding motifs, were enriched through the interactions of RID1 with the novel flowering-promoting proteins OsPIL12 and OsTCP11, respectively. Moreover, the ChIP-seq data demonstrated that RID1 bound to numerous rice heading-date genes, such as HEADING DATE 1 (HD1) and FLAVIN-BINDING, KELCH REPEAT, F-BOX 1 (OsFKF1). Notably, transcriptome sequencing (RNA-seq) analysis revealed roles of RID1 in diverse developmental pathways. Genetic analysis combined with genome-wide ChIP-seq and RNA-seq results showed that RID1 directly binds to the promoter of OsERF#136 (a repressor of rice flowering) and negatively regulates its expression. Overall, our findings provide new insights into the molecular and genetic mechanisms underlying rice floral transition and characterize OsERF#136 as a previously unrecognized direct target of RID1.

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.

Characterization of Transposon-Derived Accessible Chromatin Regions in Rice (Oryza Sativa)

Growing evidence indicates that transposons or transposable elements (TEs)-derived accessible chromatin regions (ACRs) play essential roles in multiple biological processes by interacting with trans-acting factors. However, the function of TE-derived ACRs in the regulation of gene expression in the rice genome has not been well characterized. In this study, we examined the chromatin dynamics in six types of rice tissues and found that ~8% of ACRs were derived from TEs and exhibited distinct levels of accessibility and conservation as compared to those without TEs. TEs exhibited a TE subtype-dependent impact on ACR formation, which can be mediated by changes in the underlying DNA methylation levels. Moreover, we found that tissue-specific TE-derived ACRs might function in the tissue development through the modulation of nearby gene expression. Interestingly, many genes in domestication sweeps were found to overlap with TE-derived ACRs, suggesting their potential functions in the rice domestication. In addition, we found that the expression divergence of 1070 duplicate gene pairs were associated with TE-derived ACRs and had distinct distributions of TEs and ACRs around the transcription start sites (TSSs), which may experience different selection pressures. Thus, our study provides some insights into the biological implications of TE-derived ACRs in the rice genome. Our results imply that these ACRs are likely involved in the regulation of tissue development, rice domestication and functional divergence of duplicated genes.