Enhancer transcription detected in the nascent transcriptomic landscape of bread wheat

Enhancers in Plant Development, Adaptation and Evolution

Understanding plant responses to developmental and environmental cues is crucial for studying morphological divergence and local adaptation. Gene expression changes, governed by cis-regulatory modules (CRMs) including enhancers, are a major source of plant phenotypic variation. However, while genome-wide approaches have revealed thousands of putative enhancers in mammals, far fewer have been identified and functionally characterized in plants. This review provides an overview of how enhancers function to control gene regulation, methods to predict DNA sequences that may have enhancer activity, methods utilized to functionally validate enhancers and the current knowledge of enhancers in plants, including how they impact plant development, response to environment and evolutionary adaptation.

The Enhancer-Promoter-Mediated Wnt8a Transcription During Neurite Regrowth of Injured Cortical Neurons

Brain injuries can result from accidents, warfare, sports injuries, or brain diseases. Identifying regeneration-associated genes (RAGs) during epigenome remodeling upon brain injury could have a significant impact on reducing neuronal death and subsequent neurodegeneration for patients with brain injury. We previously identified several WNT genes as RAGs involved in the neurite regrowth of injured cortical neurons. Among them, the expression of the Wnt8a gene increased most significantly during neurite regrowth, indicating its potential to promote neuronal regeneration. In this study, we investigated the regulatory mechanism of Wnt8a transcription. An algorithm was developed to predict the novel enhancer regions of candidate genes. By combining active enhancer marks, histone H3 lysine 27 acetylation (H3K27ac), and histone H3 lysine 4 mono-methylation (H3K4me1), we identified a candidate enhancer region for Wnt8a located 1.7 Mb upstream and 0.1 Mb downstream of the Wnt8a gene. This region was organized into enhancers (Ens) 1-15. Enhancer RNA expression from the predicted En1-15 regions, DNA topological dynamics, and the activity of predicted enhancers were analyzed to validate the candidate active enhancers. Our findings showed that the En8, 9, 10, 14, and 15 regions expressed higher eRNAs during neurite regrowth. Notably, the En8-2 and En14-2 subregions showed significantly up-regulated H3K4me1 modification during neurite regrowth. Using chromatin conformation capture assays and enhancer-reporter assays, we delineated that the molecular regulation of Wnt8a transcription during neurite regrowth occurs through looped En8-promoter interplay.

Integrative Multi-Omics Approaches for Identifying and Characterizing Biological Elements in Crop Traits: Current Progress and Future Prospects

The advancement of multi-omics tools has revolutionized the study of complex biological systems, providing comprehensive insights into the molecular mechanisms underlying critical traits across various organisms. By integrating data from genomics, transcriptomics, metabolomics, and other omics platforms, researchers can systematically identify and characterize biological elements that contribute to phenotypic traits. This review delves into recent progress in applying multi-omics approaches to elucidate the genetic, epigenetic, and metabolic networks associated with key traits in plants. We emphasize the potential of these integrative strategies to enhance crop improvement, optimize agricultural practices, and promote sustainable environmental management. Furthermore, we explore future prospects in the field, underscoring the importance of cutting-edge technological advancements and the need for interdisciplinary collaboration to address ongoing challenges. By bridging various omics platforms, this review aims to provide a holistic framework for advancing research in plant biology and agriculture.

Wheat2035: Integrating pan-omics and advanced biotechnology for future wheat design

Wheat (Triticum aestivum) production is vital for global food security, providing energy and protein to millions of people worldwide. Recent advancements in wheat research have led to significant increases in production, fueled by technological and scientific innovation. Here, we summarize the major advancements in wheat research, particularly the integration of biotechnologies and a deeper understanding of wheat biology. The shift from multi-omics to pan-omics approaches in wheat research has greatly enhanced our understanding of the complex genome, genomic variations, and regulatory networks to decode complex traits. We also outline key scientific questions, potential research directions, and technological strategies for improving wheat over the next decade. Since global wheat production is expected to increase by 60% in 2050, continued innovation and collaboration are crucial. Integrating biotechnologies and a deeper understanding of wheat biology will be essential for addressing future challenges in wheat production, ensuring sustainable practices and improved productivity.

Transcriptional regulation mechanism of PARP1 and its application in disease treatment

Poly (ADP-ribose) polymerase 1 (PARP1) is a multifunctional nuclear enzyme that catalyzes poly-ADP ribosylation in eukaryotic cells. In addition to maintaining genomic integrity, this nuclear enzyme is also involved in transcriptional regulation. PARP1 can trigger and maintain changes in the chromatin structure and directly recruit transcription factors. PARP1 also prevents DNA methylation. However, most previous reviews on PARP1 have focused on its involvement in maintaining genome integrity, with less focus on its transcriptional regulatory function. This article comprehensively reviews the transcriptional regulatory function of PARP1 and its application in disease treatment, providing new ideas for targeting PARP1 for the treatment of diseases other than cancer.

Enhancers associated with unstable RNAs are rare in plants

Unstable transcripts have emerged as markers of active enhancers in vertebrates and shown to be involved in many cellular processes and medical disorders. However, their prevalence and role in plants is largely unexplored. Here, we comprehensively captured all actively initiating (nascent) transcripts across diverse crops and other plants using capped small (cs)RNA sequencing. We discovered that unstable transcripts are rare in plants, unlike in vertebrates, and when present, often originate from promoters. In addition, many 'distal' elements in plants initiate tissue-specific stable transcripts and are likely bona fide promoters of as-yet-unannotated genes or non-coding RNAs, cautioning against using reference genome annotations to infer putative enhancer sites. To investigate enhancer function, we integrated data from self-transcribing active regulatory region (STARR) sequencing. We found that annotated promoters and other regions that initiate stable transcripts, but not those marked by unstable or bidirectional unstable transcripts, showed stronger enhancer activity in this assay. Our findings underscore the blurred line between promoters and enhancers and suggest that cis-regulatory elements can encompass diverse structures and mechanisms in eukaryotes, including humans.

Epigenomic identification of vernalization cis-regulatory elements in winter wheat

BACKGROUND

Winter wheat undergoes vernalization, a process activated by prolonged exposure to low temperatures. During this phase, flowering signals are generated and transported to the apical meristems, stimulating the transition to the inflorescence meristem while inhibiting tiller bud elongation. Although some vernalization genes have been identified, the key cis-regulatory elements and precise mechanisms governing this process in wheat remain largely unknown.

RESULTS

In this study, we construct extensive epigenomic and transcriptomic profiling across multiple tissues-leaf, axillary bud, and shoot apex-during the vernalization of winter wheat. Epigenetic modifications play a crucial role in eliciting tissue-specific responses and sub-genome-divergent expressions during vernalization. Notably, we observe that H3K27me3 primarily regulates vernalization-induced genes and has limited influence on vernalization-repressed genes. The integration of these datasets enables the identification of 10,600 putative vernalization-related regulatory elements including distal accessible chromatin regions (ACRs) situated 30Kb upstream of VRN3, contributing to the construction of a comprehensive regulatory network. Furthermore, we discover that TaSPL7/15, integral components of the aging-related flowering pathway, interact with the VRN1 promoter and VRN3 distal regulatory elements. These interactions finely regulate their expressions, consequently impacting the vernalization process and flowering.

CONCLUSIONS

Our study offers critical insights into wheat vernalization's epigenomic dynamics and identifies the putative regulatory elements crucial for developing wheat germplasm with varied vernalization characteristics. It also establishes a vernalization-related transcriptional network, and uncovers that TaSPL7/15 from the aging pathway participates in vernalization by directly binding to the VRN1 promoter and VRN3 distal regulatory elements.

A fine-scale Arabidopsis chromatin landscape reveals chromatin conformation-associated transcriptional dynamics

Plants, as sessile organisms, deploy transcriptional dynamics for adapting to extreme growth conditions such as cold stress. Emerging evidence suggests that chromatin architecture contributes to transcriptional regulation. However, the relationship between chromatin architectural dynamics and transcriptional reprogramming in response to cold stress remains unclear. Here, we apply a chemical-crosslinking assisted proximity capture (CAP-C) method to elucidate the fine-scale chromatin landscape, revealing chromatin interactions within gene bodies closely associated with RNA polymerase II (Pol II) densities across initiation, pausing, and termination sites. We observe dynamic changes in chromatin interactions alongside Pol II activity alterations during cold stress, suggesting local chromatin dynamics may regulate Pol II activity. Notably, cold stress does not affect large-scale chromatin conformations. We further identify a comprehensive promoter-promoter interaction (PPI) network across the genome, potentially facilitating co-regulation of gene expression in response to cold stress. Our study deepens the understanding of chromatin conformation-associated gene regulation in plant response to cold.

The role of histone acetylation in transcriptional regulation and seed development

Histone acetylation is highly conserved across eukaryotes and has been linked to gene activation since its discovery nearly 60 years ago. Over the past decades, histone acetylation has been evidenced to play crucial roles in plant development and response to various environmental cues. Emerging data indicate that histone acetylation is one of the defining features of "open chromatin," while the role of histone acetylation in transcription remains controversial. In this review, we briefly describe the discovery of histone acetylation, the mechanism of histone acetylation regulating transcription in yeast and mammals, and summarize the research progress of plant histone acetylation. Furthermore, we also emphasize the effect of histone acetylation on seed development and its potential use in plant breeding. A comprehensive knowledge of histone acetylation might provide new and more flexible research perspectives to enhance crop yield and stress resistance.

Interplay between coding and non-coding regulation drives the Arabidopsis seed-to-seedling transition

Translation of seed stored mRNAs is essential to trigger germination. However, when RNAPII re-engages RNA synthesis during the seed-to-seedling transition has remained in question. Combining csRNA-seq, ATAC-seq and smFISH in Arabidopsis thaliana we demonstrate that active transcription initiation is detectable during the entire germination process. Features of non-coding regulation such as dynamic changes in chromatin accessible regions, antisense transcription, as well as bidirectional non-coding promoters are widespread throughout the Arabidopsis genome. We show that sensitivity to exogenous ABSCISIC ACID (ABA) during germination depends on proximal promoter accessibility at ABA-responsive genes. Moreover, we provide genetic validation of the existence of divergent transcription in plants. Our results reveal that active enhancer elements are transcribed producing non-coding enhancer RNAs (eRNAs) as widely documented in metazoans. In sum, this study defining the extent and role of coding and non-coding transcription during key stages of germination expands our understanding of transcriptional mechanisms underlying plant developmental transitions.

Mapping nucleosome-resolution chromatin organization and enhancer-promoter loops in plants using Micro-C-XL

Although chromatin organizations in plants have been dissected at the scales of compartments and topologically associating domain (TAD)-like domains, there remains a gap in resolving fine-scale structures. Here, we use Micro-C-XL, a high-throughput chromosome conformation capture (Hi-C)-based technology that involves micrococcal nuclease (instead of restriction enzymes) and long cross-linkers, to dissect single nucleosome-resolution chromatin organization in Arabidopsis. Insulation analysis reveals more than 14,000 boundaries, which mostly include chromatin accessibility, epigenetic modifications, and transcription factors. Micro-C-XL reveals associations between RNA Pols and local chromatin organizations, suggesting that gene transcription substantially contributes to the establishment of local chromatin domains. By perturbing Pol II both genetically and chemically at the gene level, we confirm its function in regulating chromatin organization. Visible loops and stripes are assigned to super-enhancers and their targeted genes, thus providing direct insights for the identification and mechanistic analysis of distal CREs and their working modes in plants. We further investigate possible factors regulating these chromatin loops. Subsequently, we expand Micro-C-XL to soybean and rice. In summary, we use Micro-C-XL for analyses of plants, which reveal fine-scale chromatin organization and enhancer-promoter loops and provide insights regarding three-dimensional genomes in plants.

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.

Comparative analysis of nascent RNA sequencing methods and their applications in studies of cotranscriptional splicing dynamics

High-throughput detection of nascent RNA is critical for studies of transcription and much more challenging than that of mRNA. Recently, several massively parallel nascent RNA sequencing methods were established in eukaryotic cells. Here, we systematically compared 3 classes of methods on the same pure or crude nuclei preparations: GRO-seq for sequence nuclear run-on RNAs, pNET-seq for sequence RNA polymerase II-associated RNAs, and CB RNA-seq for sequence chromatin-bound (CB) RNAs in Arabidopsis (Arabidopsis thaliana). To improve the resolution of CB RNAs, 3'CB RNA-seq was established to sequence the 3' ends of CB RNAs. In addition, we modified pNET-seq to establish the Chromatin Native Elongation Transcript sequencing (ChrNET) method using chromatin as the starting material for RNA immunoprecipitation. Reproducibility, sensitivity and accuracy in detecting nascent transcripts, experimental procedures, and costs were analyzed, which revealed the strengths and weaknesses of each method. We found that pNET and GRO methods best detected active RNA polymerase II. CB RNA-seq is a simple and cost-effective alternative for nascent RNA studies, due to its high correlation with pNET-seq and GRO-seq. Compared with pNET, ChrNET has higher specificity for nascent RNA capture and lower sequencing cost. 3'CB is sensitive to transcription-coupled splicing. Using these methods, we identified 1,404 unknown transcripts, 4,482 unannotated splicing events, and 60 potential recursive splicing events. This comprehensive comparison of different nascent/chromatin RNA sequencing methods highlights the strengths of each method and serves as a guide for researchers aiming to select a method that best meets their study goals.

Transposable element-initiated enhancer-like elements generate the subgenome-biased spike specificity of polyploid wheat

Transposable elements (TEs) comprise ~85% of the common wheat genome, which are highly diverse among subgenomes, possibly contribute to polyploid plasticity, but the causality is only assumed. Here, by integrating data from gene expression cap analysis and epigenome profiling via hidden Markov model in common wheat, we detect a large proportion of enhancer-like elements (ELEs) derived from TEs producing nascent noncoding transcripts, namely ELE-RNAs, which are well indicative of the regulatory activity of ELEs. Quantifying ELE-RNA transcriptome across typical developmental stages reveals that TE-initiated ELE-RNAs are mainly from RLG_famc7.3 specifically expanded in subgenome A. Acquisition of spike-specific transcription factor binding likely confers spike-specific expression of RLG_famc7.3-initiated ELE-RNAs. Knockdown of RLG_famc7.3-initiated ELE-RNAs resulted in global downregulation of spike-specific genes and abnormal spike development. These findings link TE expansion to regulatory specificity and polyploid developmental plasticity, highlighting the functional impact of TE-driven regulatory innovation on polyploid evolution.

Genome-wide enhancer identification by massively parallel reporter assay in Arabidopsis

Enhancers are critical cis-regulatory elements controlling gene expression during cell development and differentiation. However, genome-wide enhancer characterization has been challenging due to the lack of a well-defined relationship between enhancers and genes. Function-based methods are the gold standard for determining the biological function of cis-regulatory elements; however, these methods have not been widely applied to plants. Here, we applied a massively parallel reporter assay on Arabidopsis to measure enhancer activities across the genome. We identified 4327 enhancers with various combinations of epigenetic modifications distinctively different from animal enhancers. Furthermore, we showed that enhancers differ from promoters in their preference for transcription factors. Although some enhancers are not conserved and overlap with transposable elements forming clusters, enhancers are generally conserved across thousand Arabidopsis accessions, suggesting they are selected under evolution pressure and could play critical roles in the regulation of important genes. Moreover, comparison analysis reveals that enhancers identified by different strategies do not overlap, suggesting these methods are complementary in nature. In sum, we systematically investigated the features of enhancers identified by functional assay in A. thaliana, which lays the foundation for further investigation into enhancers' functional mechanisms in plants.

nASAP: A Nascent RNA Profiling Data Analysis Platform

Although nascent RNA profiling data are widely used in transcriptional regulation studies, the development and standardization of data processing pipeline lags far behind RNA-seq. We are filling this gap by establishing the nASAP web server (https://grobase.top/nasap/) to provide practical quality evaluation and comprehensive analysis of nascent RNA datasets. In nASAP, four customized analysis modules are provided, including i) quality assessment, which summarizes the sequencing statistics, mapping ratio, and evaluates RNA integrity and mRNA contamination; ii) quantification analysis for mRNAs, lncRNAs and eRNAs; iii) pausing analysis across the whole genome based on sequencing reads distribution; and iv) network analysis to better understand the gene regulatory mechanism by obtaining annotated enhancer-promoter interactomes. The nASAP is user-friendly and outperforms the existing pipeline for quality control of nascent RNA profiling data. We anticipate that nASAP, which eases both basic and advanced analysis of nascent RNA data, will be extremely useful in various fields.

HSFA1a modulates plant heat stress responses and alters the 3D chromatin organization of enhancer-promoter interactions

The complex and dynamic three-dimensional organization of chromatin within the nucleus makes understanding the control of gene expression challenging, but also opens up possible ways to epigenetically modulate gene expression. Because plants are sessile, they evolved sophisticated ways to rapidly modulate gene expression in response to environmental stress, that are thought to be coordinated by changes in chromatin conformation to mediate specific cellular and physiological responses. However, to what extent and how stress induces dynamic changes in chromatin reorganization remains poorly understood. Here, we comprehensively investigated genome-wide chromatin changes associated with transcriptional reprogramming response to heat stress in tomato. Our data show that heat stress induces rapid changes in chromatin architecture, leading to the transient formation of promoter-enhancer contacts, likely driving the expression of heat-stress responsive genes. Furthermore, we demonstrate that chromatin spatial reorganization requires HSFA1a, a transcription factor (TF) essential for heat stress tolerance in tomato. In light of our findings, we propose that TFs play a key role in controlling dynamic transcriptional responses through 3D reconfiguration of promoter-enhancer contacts.

Identification and functional validation of super-enhancers in Arabidopsis thaliana

Super-enhancers (SEs) are exceptionally large enhancers and are recognized to play prominent roles in cell identity in mammalian species. We surveyed the genomic regions containing large clusters of accessible chromatin regions (ACRs) marked by deoxyribonuclease (DNase) I hypersensitivity in Arabidopsis thaliana. We identified a set of 749 putative SEs, which have a minimum length of 1.5 kilobases and represent the top 2.5% of the largest ACR clusters. We demonstrate that the genomic regions associating with these SEs were more sensitive to DNase I than other nonpromoter ACRs. The SEs were preferentially associated with topologically associating domains. Furthermore, the SEs and their predicted cognate genes were frequently associated with organ development and tissue identity in A. thaliana. Therefore, the A. thaliana SEs and their cognate genes mirror the functional characteristics of those reported in mammalian species. We developed CRISPR/Cas-mediated deletion lines of a 3,578-bp SE associated with the thalianol biosynthetic gene cluster (BGC). Small deletions (131-157 bp) within the SE resulted in distinct phenotypic changes and transcriptional repression of all five thalianol genes. In addition, T-DNA insertions in the SE region resulted in transcriptional alteration of all five thalianol genes. Thus, this SE appears to play a central role in coordinating the operon-like expression pattern of the thalianol BGC.

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.

Protocol for affordable and efficient profiling of nascent RNAs in bread wheat using GRO-seq

Exorbitant sequencing cost is one of the main obstacles limiting the widespread application of Global Run-On sequencing (GRO-seq) to detect transcriptional activity. Here, we describe a more efficient and affordable protocol for GRO-seq that incorporates an rRNA removal step after nuclear RNA isolation and before nascent RNA immunoprecipitation. We have successfully applied this protocol to profile enhancer transcription in allohexaploid bread wheat and increased the proportion of valid data by 20 times.

For complete details on the use and execution of this protocol, please refer to Xie et al. (2022).

•

The detailed GRO-seq protocol integrates rRNA depletion

•

Cost-efficient GRO-seq protocol for detection of eRNA in bread wheat

•

Applicable for any large complex plant or animal genomes

Publisher’s note: Undertaking any experimental protocol requires adherence to local institutional guidelines for laboratory safety and ethics.

Characterization of Expression and Epigenetic Features of Core Genes in Common Wheat

The availability of multiple wheat genome sequences enables us to identify core genes and characterize their genetic and epigenetic features, thereby advancing our understanding of their biological implications within individual plant species. It is, however, largely understudied in wheat. To this end, we reanalyzed genome sequences from 16 different wheat varieties and identified 62,299 core genes. We found that core and non-core genes have different roles in subgenome differentiation. Meanwhile, according to their expression profiles, these core genes can be classified into genes related to tissue development and stress responses, including 3376 genes highly expressed in both spikelets and at high temperatures. After associating with six histone marks and open chromatin, we found that these core genes can be divided into eight sub-clusters with distinct epigenomic features. Furthermore, we found that ca. 51% of the expressed transcription factors (TFs) were marked with both H3K27me3 and H3K4me3, indicative of the bivalency feature, which can be involved in tissue development through the TF-centered regulatory network. Thus, our study provides a valuable resource for the functional characterization of core genes in stress responses and tissue development in wheat.