A Quantitative Computational Framework for Allopolyploid Single-Cell Data Integration and Core Gene Ranking in Development.

Polyploidization drives regulatory and phenotypic innovation. How the merger of different genomes contributes to polyploid development is a fundamental issue in evolutionary developmental biology and breeding research. Clarifying this issue is challenging because of genome complexity and the difficulty in tracking stochastic subgenome divergence during development. Recent single-cell sequencing techniques enabled probing subgenome-divergent regulation in the context of cellular differentiation. However, analyzing single-cell data suffers from high error rates due to high dimensionality, noise, and sparsity, and the errors stack up in polyploid analysis due to the increased dimensionality of comparisons between subgenomes of each cell, hindering deeper mechanistic understandings. In this study, we develop a quantitative computational framework, called "pseudo-genome divergence quantification" (pgDQ), for quantifying and tracking subgenome divergence directly at the cellular level. Further comparing with cellular differentiation trajectories derived from single-cell RNA sequencing data allows for an examination of the relationship between subgenome divergence and the progression of development. pgDQ produces robust results and is insensitive to data dropout and noise, avoiding high error rates due to multiple comparisons of genes, cells, and subgenomes. A statistical diagnostic approach is proposed to identify genes that are central to subgenome divergence during development, which facilitates the integration of different data modalities, enabling the identification of factors and pathways that mediate subgenome-divergent activity during development. Case studies have demonstrated that applying pgDQ to single-cell and bulk tissue transcriptomic data promotes a systematic and deeper understanding of how dynamic subgenome divergence contributes to developmental trajectories in polyploid evolution.

Evolutionary rewiring of the wheat transcriptional regulatory network by lineage-specific transposable elements.

More than 80% of the wheat genome consists of transposable elements (TEs), which act as major drivers 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 (TFBSs) 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 TFBSs share significantly high sequence similarity with TE-embedded TFBSs. These non-TE TFBSs 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.

An atlas of wheat epigenetic regulatory elements reveals subgenome divergence in the regulation of development and stress responses.

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 is 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, and demonstrated the coordinated interplay between epiRE sequence contexts, epigenetic factors, and transcription factors in regulating subgenome divergent transcriptional responses to external changes. Together, 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.

Functional Characterization of Accessible Chromatin in Common Wheat.

Eukaryotic gene transcription is fine-tuned by precise spatiotemporal interactions between cis-regulatory elements (CREs) and trans-acting factors. However, how CREs individually or coordinated with epigenetic marks function in regulating homoeolog bias expression is still largely unknown in wheat. In this study, through comprehensively characterizing open chromatin coupled with DNA methylation in the seedling and spikelet of common wheat, we observed that differential chromatin openness occurred between the seedling and spikelet, which plays important roles in tissue development through regulating the expression of related genes or through the transcription factor (TF)-centered regulatory network. Moreover, we found that CHH methylation may act as a key determinant affecting the differential binding of TFs, thereby resulting in differential expression of target genes. In addition, we found that sequence variations in MNase hypersensitive sites (MHSs) result in the differential expression of key genes responsible for important agronomic traits. Thus, our study provides new insights into the roles of CREs in regulating tissue or homoeolog bias expression, and controlling important agronomic traits in common wheat. It also provides potential CREs for genetic and epigenetic manipulation toward improving desirable traits for wheat molecule breeding.

Characterization of regulatory modules controlling leaf angle in maize.

Leaf angle is an important agronomic trait determining maize (Zea mays) planting density and light penetration into the canopy and contributes to the yield gain in modern maize hybrids. However, little is known about the molecular mechanisms underlying leaf angle beyond the ZmLG1 (liguleless1) and ZmLG2 (Liguleless2) genes. In this study, we found that the transcription factor (TF) ZmBEH1 (BZR1/BES1 homolog gene 1) is targeted by ZmLG2 and regulates leaf angle formation by influencing sclerenchyma cell layers on the adaxial side. ZmBEH1 interacted with the TF ZmBZR1 (Brassinazole Resistant 1), whose gene expression was also directly activated by ZmLG2. Both ZmBEH1 and ZmBZR1 are bound to the promoter of ZmSCL28 (SCARECROW-LIKE 28), a third TF that influences leaf angle. Our study demonstrates regulatory modules controlling leaf angle and provides gene editing targets for creating optimal maize architecture suitable for dense planting.

Characterization of functional relationships of R-loops with gene transcription and epigenetic modifications in rice.

We conducted genome-wide identification of R-loops followed by integrative analyses of R-loops with relation to gene expression and epigenetic signatures in the rice genome. We found that the correlation between gene expression levels and profiled R-loop peak levels was dependent on the positions of R-loops within gene structures (hereafter named "genic position"). Both antisense only (ASO)-R-loops and sense/antisense (S/AS)-R-loops sharply peaked around transcription start sites (TSSs), and these peak levels corresponded positively with transcript levels of overlapping genes. In contrast, sense only (SO)-R-loops were generally spread over the coding regions, and their peak levels corresponded inversely to transcript levels of overlapping genes. In addition, integrative analyses of R-loop data with existing RNA-seq, chromatin immunoprecipitation sequencing (ChIP-seq), DNase I hypersensitive sites sequencing (DNase-seq), and whole-genome bisulfite sequencing (WGBS or BS-seq) data revealed interrelationships and intricate connections among R-loops, gene expression, and epigenetic signatures. Experimental validation provided evidence that the demethylation of both DNA and histone marks can influence R-loop peak levels on a genome-wide scale. This is the first study in plants that reveals novel functional aspects of R-loops, their interrelations with epigenetic methylation, and roles in transcriptional regulation.

Salt-Responsive Genes are Differentially Regulated at the Chromatin Levels Between Seedlings and Roots in Rice.

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.

The bread wheat epigenomic map reveals distinct chromatin architectural and evolutionary features of functional genetic elements.

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.