Search

Abstract Background In response to drought stress (DS), plants undergo complex processes that entail significant transcriptome reprogramming. However, the intricate relationship between the dynamic alterations in the three-dimensional (3D) genome and the modulation of gene co-expression in drought responses remains a relatively unexplored area. Results In this study, we reconstruct high-resolution 3D genome maps based on genomic regions marked by H3K9ac, an active histone modification that dynamically responds to soil water variations in rice. We discover a genome-wide disconnection of 3D genome contact upon DS with over 10,000 chromatin loops lost, which are partially recovered in the subsequent re-watering. Loops integrating promoter–promoter interactions (PPI) contribute to gene expression in addition to basal H3K9ac modifications. Moreover, H3K9ac-marked promoter regions with high affinities in mediating PPIs, termed as super-promoter regions (SPRs), integrate spatially clustered PPIs in a super-enhancer-like manner. Interestingly, the knockout mutation of OsbZIP23, a well-defined DS-responsive transcription factor, leads to the disassociation of over 80% DS-specific PPIs and decreased expression of the corresponding genes under DS. As a case study, we show how OsbZIP23 integrates the PPI cluster formation and the co-expression of four dehydrin genes, RAB16A–D, through targeting the RAB16C SPR in a stress signaling-dependent manner. Conclusions Our high-resolution 3D genome maps unveil the principles and details of dynamic genome folding in response to water supply variations and illustrate OsbZIP23 as an indispensable integrator of the yet unique 3D genome organization that is essential for gene co-expression under DS in rice.

Abstract Background Epigenetic modifications that exhibit circadian oscillations also promote circadian oscillations of gene expression. Brassica napus is a heterozygous polyploid species that has undergone distant hybridization and genome doubling events and has a young and distinct species origin. Studies incorporating circadian rhythm analysis of epigenetic modifications can offer new insights into differences in diurnal oscillation behavior among subgenomes and the regulation of diverse expressions of homologous gene rhythms in biological clocks. Results In this study, we created a high-resolution and multioscillatory gene expression dataset, active histone modification (H3K4me3, H3K9ac), and RNAPII recruitment in Brassica napus. We also conducted the pioneering characterization of the diurnal rhythm of transcription and epigenetic modifications in an allopolyploid species. We compared the evolution of diurnal rhythms between subgenomes and observed that the Cn subgenome had higher diurnal oscillation activity in both transcription and active histone modifications than the An subgenome. Compared to the A subgenome in Brassica rapa, the An subgenome of Brassica napus displayed significant changes in diurnal oscillation characteristics of transcription. Homologous gene pairs exhibited a higher proportion of diurnal oscillation in transcription than subgenome-specific genes, attributed to higher chromatin accessibility and abundance of active epigenetic modification types. We found that the diurnal expression of homologous genes displayed diversity, and the redundancy of the circadian system resulted in extensive changes in the diurnal rhythm characteristics of clock genes after distant hybridization and genome duplication events. Epigenetic modifications influenced the differences in the diurnal rhythm of homologous gene expression, and the diurnal oscillation of homologous gene expression was affected by the combination of multiple histone modifications. Conclusions Herein, we presented, for the first time, a characterization of the diurnal rhythm characteristics of gene expression and its epigenetic modifications in an allopolyploid species. Our discoveries shed light on the epigenetic factors responsible for the diurnal oscillation activity imbalance between subgenomes and homologous genes’ rhythmic expression differences. The comprehensive time-series dataset we generated for gene expression and epigenetic modifications provides a valuable resource for future investigations into the regulatory mechanisms of protein-coding genes in Brassica napus.

Summary: 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.

Comprehensive epigenomic maps of various rice varieties are still unavailable. Here, the authors report the development of eChIP as a fast and low-input upgrade of regular plant ChIP-seq protocol for epigenome analysis of 20 rice varieties and annotate over 80% of the genome with different epigenome properties for transcriptional regulation.

Three-dimensional genome organization and its effect on transcription remain elusive in rice. Here, the authors map promoter–promoter interactions and heterochromatin interactions using ChIA-PET and reveal spatial correlation between the genetic regulation of eQTLs and e-traits.

Spatial organization of regulatory elements and its impact on gene expression in plants remain unclear. Here, the authors construct maize chromatin interaction maps using chromatin interaction analysis by paired-end tag sequencing (ChIA-PET) and show their associations with gene expression and agronomic traits.

Characterizing genome-wide histone posttranscriptional modifications and transcriptional factor occupancy is crucial for deciphering their biological functions. Chromatin immunoprecipitation followed by sequencing (ChIP-seq) is a powerful method for genome-wide profiling of histone modifications and transcriptional factor-binding sites. However, the current ChIP-seq experimental procedure in plants requires significant material and several days for completion. CUT&Tag is an alternative method of ChIP-seq for low-sample and single-cell epigenomic profiling using protein A-Tn5 transposase fusion proteins (PAT). In this study, we developed a nucleus CUT&Tag (nCUT&Tag) protocol based on the live-cell CUT&Tag technology. Our results indicate that nCUT&Tag could be used for histone modifications profiling in both monocot rice and dicot rapeseed using crosslinked or fresh tissues. In addition, both active and repressive histone marks such as H3K4me3 and H3K9me2 can be identified using our nCUT&Tag. More importantly, all the steps in nCUT&Tag can be finished in only 1 day, and the assay can be performed with as little as 0.01 g of plant tissue as starting materials. Therefore, our results demonstrate that nCUT&Tag is an efficient alternative strategy for plant epigenomic studies.