Background Bread wheat is an allohexaploid species with a 16-Gb genome that has large intergenic regions, which presents a big challenge for pinpointing regulatory elements and further revealing the transcriptional regulatory mechanisms. Chromatin profiling to characterize the combinatorial patterns of chromatin signatures is a powerful means to detect functional elements and clarify regulatory activities in human studies. Results In the present study, through comprehensive analyses of the open chromatin, DNA methylome, seven major chromatin marks, and transcriptomic data generated for seedlings of allohexaploid wheat, we detected distinct chromatin architectural features surrounding various functional elements, including genes, promoters, enhancer-like elements, and transposons. Thousands of new genic regions and cis-regulatory elements are identified based on the combinatorial pattern of chromatin features. Roughly 1.5% of the genome encodes a subset of active regulatory elements, including promoters and enhancer-like elements, which are characterized by a high degree of chromatin openness and histone acetylation, an abundance of CpG islands, and low DNA methylation levels. A comparison across sub-genomes reveals that evolutionary selection on gene regulation is targeted at the sequence and chromatin feature levels. The divergent enrichment of cis-elements between enhancer-like sequences and promoters implies these functional elements are targeted by different transcription factors. Conclusions We herein present a systematic epigenomic map for the annotation of cis-regulatory elements in the bread wheat genome, which provides new insights into the connections between chromatin modifications and cis-regulatory activities in allohexaploid wheat. Electronic supplementary material The online version of this article (10.1186/s13059-019-1746-8) contains supplementary material, which is available to authorized users.

Wheat (Triticum aestivum) has a large allohexaploid genome. Subgenome-divergent regulation contributed to genome plasticity and the domestication of polyploid wheat. However, the specificity encoded in the wheat genome determining subgenome-divergent spatio-temporal regulation has been largely unexplored. The considerable size and complexity of the genome are major obstacles to dissecting the regulatory specificity. Here, we compared the epigenomes and transcriptomes from a large set of samples under diverse developmental and environmental conditions. Thousands of distal epigenetic regulatory elements (distal-epiREs) were specifically linked to their target promoters with coordinated epigenomic changes. We revealed that subgenome-divergent activity of homologous regulatory elements are affected by specific epigenetic signatures. Subgenome-divergent epiRE regulation of tissue specificity is associated with dynamic modulation of H3K27me3 mediated by Polycomb complex and demethylases. Furthermore, quantitative epigenomic approaches detected key stress responsive cis- and trans-acting factors validated by DNA Affinity Purification and sequencing (DAP-seq), and demonstrated the coordinated interplay between epiRE sequence contexts, epigenetic factors, and transcription factors in regulating subgenome divergent transcriptional responses to external changes. Thus, this study provides a wealth of resources for elucidating the epiRE regulomics and subgenome-divergent regulation in hexaploid wheat, and gives new clues for interpreting genetic and epigenetic interplay in regulating the benefits of polyploid wheat.

The elucidation of epigenetic responses of salt-responsive genes facilitates understanding of the underlying mechanisms that confer salt tolerance in rice. However, it is still largely unknown how epigenetic mechanisms are associated with the expression of salt-responsive genes in rice and other crops. In this study, we reported tissue-specific gene expression and tissue-specific changes in chromatin modifications or signatures between seedlings and roots in response to salt treatment. Our study indicated that among six of individual mark examined (H3K4me3, H3K27me3, H4K12ac, H3K9ac, H3K27ac and H3K36me3), a positive association between salt-related changes in histone marks and the expression of differentially expressed genes (DEGs) was observed only for H3K9ac and H4K12ac in seedlings and H3K36me3 in roots. In contrast, chromatin states (CSs) with combinations of six histone modification marks played crucial roles in the differential expression of salt-responsive genes between seedlings and roots. Most importantly, CS7 containing the bivalent marks H3K4me3 and H3K27me3, with a mutual exclusion of functions with each other, displayed distinct functions in the expression of DEGs in both tissues. Specifically, H3K27me3 in CS7 mainly suppressed the expression of DEGs in roots, while H3K4me3 affected the expression of down- and up-regulated genes, possibly by antagonizing the repressive role of H3K27me3 in seedlings. Our findings indicate distinct impacts of the CSs on the differential expression of salt-responsive genes between seedlings and roots in rice, which provides an important background for understanding chromatin-based epigenetic mechanisms that might confer salt tolerance in plants.

We thank the Bioinformatics Center, Nanjing Agricultural University for providing computing facilities for data processing and analyses. This research was partially supported by the China Scholarship Council (grant no. 201706850071), which awarded a fellowship to Y.F. Author contributions: W.Z. conceived and designed the study; W.Z. developed the DRIP-seq methodology and generated DRIPseq and NBS-seq data sets. Y.F. analyzed the data; P.Z. analyzed the 6mA data. L.C. prepared all transgenic rice lines and conducted G4-ChIP-qPCR and histone-related western blotting experiments. K.L. grew the material and performed the 6mA and 5mC dot blotting assays. Y.L.F. prepared the genomic DNA and conducted the G4-ChIP experiments. J.S. proofread the manuscript. X.P. performed the qPCR and dot blotting assays. C.C. assisted in the NBS-seq library sequencing and provided 6mA-related cas9 rice transgenic lines. X.W. helped with DNA fragmentation. Z.S. helped with rice callus preparation. Y.W. and H.W. supervised the data analysis. W.Z., Y.F., and H.W. interpreted the results. W.Z. and H.W. wrote the paper. All authors read and approved the final version of the manuscript.

More than 80% of the wheat genome consists of transposable elements (TEs), which act as one major driver of wheat genome evolution. However, their contributions to the regulatory evolution of wheat adaptations remain largely unclear. Here, we created genome-binding maps for 53 transcription factors (TFs) underlying environmental responses by leveraging DAP-seq in Triticum urartu, together with epigenomic profiles. Most TF-binding sites (TFBS) located distally from genes are embedded in TEs, whose functional relevance is supported by purifying selection and active epigenomic features. About 24% of the non-TE TFBS share significantly high sequence similarity with TE-embedded TFBS. These non-TE TFBS have almost no homologous sequences in non-Triticeae species and are potentially derived from Triticeae-specific TEs. The expansion of TE-derived TFBS linked to wheat-specific gene responses, suggesting TEs are an important driving force for regulatory innovations. Altogether, TEs have been significantly and continuously shaping regulatory networks related to wheat genome evolution and adaptation.