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.

Diurnal oscillations of epigenetic modifications are associated with variation in rhythmic expression of homoeologous genes in Brassica napus.

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.

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.

Transcriptional and translational landscape fine-tune genome annotation and explores translation control in cotton.

INTRODUCTION: The unavailability of intergenic region annotation in whole genome sequencing and pan-genomics hinders efforts to enhance crop improvement. OBJECTIVES: Despite advances in research, the impact of post-transcriptional regulation on fiber development and translatome profiling at different stages of fiber growth in cotton (G. hirsutum) remains unexplored. METHODS: We utilized a combination of reference-guided de novo transcriptome assembly and ribosome profiling techniques to uncover the hidden mechanisms of translational control in eight distinct tissues of upland cotton. RESULTS: Our study identified P-site distribution at three-nucleotide periodicity and dominant ribosome footprint at 27 nucleotides. Specifically, we have detected 1,589 small open reading frames (sORFs), including 1,376 upstream ORFs (uORFs) and 213 downstream ORFs (dORFs), as well as 552 long non-coding RNAs (lncRNAs) with potential coding functions, which fine-tune the annotation of the cotton genome. Further, we have identified novel genes and lncRNAs with strong translation efficiency (TE), while sORFs were found to affect mRNA transcription levels during fiber elongation. The reliability of these findings was confirmed by the high consistency in correlation and synergetic fold change between RNA-sequencing (RNA-seq) and Ribosome-sequencing (Ribo-seq) analyses. Additionally, integrated omics analysis of the normal fiber ZM24 and short fiber pag1 cotton mutant revealed several differentially expressed genes (DEGs), and fiber-specific expressed (high/low) genes associated with sORFs (uORFs and dORFs). These findings were further supported by the overexpression and knockdown of GhKCS6, a gene associated with sORFs in cotton, and demonstrated the potential regulation of the mechanism governing fiber elongation on both the transcriptional and post-transcriptional levels. CONCLUSION: Reference-guided transcriptome assembly and the identification of novel transcripts fine-tune the annotation of the cotton genome and predicted the landscape of fiber development. Our approach provided a high-throughput method, based on multi-omics, for discovering unannotated ORFs, hidden translational control, and complex regulatory mechanisms in crop plants.

Diurnal RNAPII-tethered chromatin interactions are associated with rhythmic gene expression in rice.

BACKGROUND: The daily cycling of plant physiological processes is speculated to arise from the coordinated rhythms of gene expression. However, the dynamics of diurnal 3D genome architecture and their potential functions underlying the rhythmic gene expression remain unclear. RESULTS: Here, we reveal the genome-wide rhythmic occupancy of RNA polymerase II (RNAPII), which precedes mRNA accumulation by approximately 2 h. Rhythmic RNAPII binding dynamically correlates with RNAPII-mediated chromatin architecture remodeling at the genomic level of chromatin interactions, spatial clusters, and chromatin connectivity maps, which are associated with the circadian rhythm of gene expression. Rhythmically expressed genes within the same peak phases of expression are preferentially tethered by RNAPII for coordinated transcription. RNAPII-associated chromatin spatial clusters (CSCs) show high plasticity during the circadian cycle, and rhythmically expressed genes in the morning phase and non-rhythmically expressed genes in the evening phase tend to be enriched in RNAPII-associated CSCs to orchestrate expression. Core circadian clock genes are associated with RNAPII-mediated highly connected chromatin connectivity networks in the morning in contrast to the scattered, sporadic spatial chromatin connectivity in the evening; this indicates that they are transcribed within physical proximity to each other during the AM circadian window and are located in discrete "transcriptional factory" foci in the evening, linking chromatin architecture to coordinated transcription outputs. CONCLUSION: Our findings uncover fundamental diurnal genome folding principles in plants and reveal a distinct higher-order chromosome organization that is crucial for coordinating diurnal dynamics of transcriptional regulation.