Genome-Wide Prediction and Validation of Intergenic Enhancers in Arabidopsis Using Open Chromatin Signatures

Plant 3D genomics: the exploration and application of chromatin organization

Eukaryotic genomes are highly folded for packing into higher-order chromatin structures in the nucleus. With the emergence of state-of-the-art chromosome conformation capture methods and microscopic imaging techniques, the spatial organization of chromatin and its functional implications have been interrogated. Our knowledge of 3D chromatin organization in plants has improved dramatically in the past few years, building on the early advances in animal systems. Here, we review recent advances in 3D genome mapping approaches, our understanding of the sophisticated organization of spatial structures, and the application of 3D genomic principles in plants. We also discuss directions for future developments in 3D genomics in plants.

Enhancer-mediated reporter gene expression in Arabidopsis thaliana: a forward genetic screen

Gene expression is controlled and regulated by interactions between cis-regulatory DNA elements (CREs) and regulatory proteins. Enhancers are one of the most important classes of CREs in eukaryotes. Eukaryotic genes, especially those related to development or responses to environmental cues, are often regulated by multiple enhancers in different tissues and/or at different developmental stages. Remarkably, little is known about the molecular mechanisms by which enhancers regulate gene expression in plants. We identified a distal enhancer, CREβ, which regulates the expression of AtDGK7, which encodes a diacylglycerol kinase in Arabidopsis. We developed a transgenic line containing the luciferase reporter gene (LUC) driven by CREβ fused with a minimal cauliflower mosaic virus (CaMV) 35S promoter. The CREβ enhancer was shown to play a role in the response to osmotic pressure of the LUC reporter gene. A forward genetic screen pipeline based on the transgenic line was established to generate mutations associated with altered expression of the LUC reporter gene. We identified a suite of mutants with variable LUC expression levels as well as different segregation patterns of the mutations in populations. We demonstrate that this pipeline will allow us to identify trans-regulatory factors associated with CREβ function as well as those acting in the regulation of the endogenous AtDGK7 gene.

Chromatin accessibility and translational landscapes of tea plants under chilling stress

Plants have evolved regulatory mechanisms at multiple levels to regulate gene expression in order to improve their cold adaptability. However, limited information is available regarding the stress response at the chromatin and translational levels. Here, we characterize the chromatin accessibility, transcriptional, and translational landscapes of tea plants in vivo under chilling stress for the first time. Chilling stress significantly affected both the transcription and translation levels as well as the translation efficiency of tea plants. A total of 3010 genes that underwent rapid and independent translation under chilling stress were observed, and they were significantly enriched in the photosynthesis-antenna protein and phenylpropanoid biosynthesis pathways. A set of genes that were significantly responsive to cold at the transcription and translation levels, including four (+)-neomenthol dehydrogenases (MNDs) and two (E)-nerolidol synthases (NESs) arranged in tandem on the chromosomes, were also found. We detected potential upstream open reading frames (uORFs) on 3082 genes and found that tea plants may inhibit the overall expression of genes by enhancing the translation of uORFs under chilling stress. In addition, we identified distal transposase hypersensitive sites (THSs) and proximal THSs and constructed a transcriptional regulatory network for tea plants under chilling stress. We also identified 13 high-confidence transcription factors (TFs) that may play a crucial role in cold regulation. These results provide valuable information regarding the potential transcriptional regulatory network in plants and help to clarify how plants exhibit flexible responses to chilling stress.

Genomic Features of Open Chromatin Regions (OCRs) in Wild Soybean and Their Effects on Gene Expressions

Transcription activation is tightly associated with the openness of chromatin, which allows direct contact between transcriptional regulators, such as transcription factors, and their targeted DNA for downstream gene activation. However, the annotation of open chromatin regions (OCRs) in the wild soybean (Glycine soja) genome is limited. We performed assay for transposase-accessible chromatin using sequencing (ATAC-seq) and successfully identified 22,333 OCRs in the leaf of W05 (a wild soybean accession). These OCRs were enriched in gene transcription start sites (TSS) and were positively correlated with downstream gene expression. Several known transcription factor (TF)-binding motifs were also enriched at the OCRs. A potential regulatory network was constructed using these transcription factors and the OCR-marked genes. Furthermore, by overlapping the OCR distribution with those of histone modifications from chromatin immunoprecipitation followed by sequencing (ChIP-seq), we found that the distribution of the activation histone mark, H3K4me3, but not that of the repressive H3K27me3 mark, was closely associated with OCRs for gene activation. Several putative enhancer-like distal OCRs were also found to overlap with LincRNA-encoding loci. Moreover, our data suggest that homologous OCRs could potentially influence homologous gene expression. Hence, the duplication of OCRs might be essential for plant genome architecture as well as for regulating gene expression.

Identification of Key Tissue-Specific, Biological Processes by Integrating Enhancer Information in Maize Gene Regulatory Networks

Enhancers are key players in the spatio-temporal coordination of gene expression during numerous crucial processes, including tissue differentiation across development. Characterizing the transcription factors (TFs) and genes they connect, and the molecular functions underpinned is important to better characterize developmental processes. In plants, the recent molecular characterization of enhancers revealed their capacity to activate the expression of several target genes. Nevertheless, identifying these target genes at a genome-wide level is challenging, particularly for large-genome species, where enhancers and target genes can be hundreds of kilobases away. Therefore, the contribution of enhancers to plant regulatory networks remains poorly understood. Here, we investigate the enhancer-driven regulatory network of two maize tissues at different stages: leaves at seedling stage (V2-IST) and husks (bracts) at flowering. Using systems biology, we integrate genomic, epigenomic, and transcriptomic data to model the regulatory relationships between TFs and their potential target genes, and identify regulatory modules specific to husk and V2-IST. We show that leaves at the V2-IST stage are characterized by the response to hormones and macromolecules biogenesis and assembly, which are regulated by the BBR/BPC and AP2/ERF TF families, respectively. In contrast, husks are characterized by cell wall modification and response to abiotic stresses, which are, respectively, orchestrated by the C2C2/DOF and AP2/EREB families. Analysis of the corresponding enhancer sequences reveals that two different transposable element families (TIR transposon Mutator and MITE Pif/Harbinger) have shaped part of the regulatory network in each tissue, and that MITEs have provided potential new TF binding sites involved in husk tissue-specificity.

Rapid Validation of Transcriptional Enhancers Using a Transient Reporter Assay

Enhancers are one of the main classes of cis-regulatory elements (CREs) in the regulation of plant gene expression. Plant enhancers can be predicted based on genomic signatures associated with open chromatin. However, predicted enhancers need to be validated experimentally. We developed an experimental system for rapid enhancer validation. Predicted enhancer candidates are cloned into a vector containing a minimal 35S promoter and a luciferase reporter gene. The construct is then agroinfiltrated into Nicotiana benthamiana leaves followed by bioluminescence signal detection and analysis. Positive bioluminescence signals indicate the enhancer function of each candidate, and the relative signal strength from different enhancers can be quantitatively measured and compared. In summary, we have developed an efficient and rapid plant enhancer validation assay based on a bioluminescent luciferase reporter and agroinfiltration-based N. benthamiana leaf transient expression. This assay can be used for the initial screening of candidate enhancers that are active in leaf tissue. The system can potentially be used to examine the activity of candidate enhancers under different environmental conditions.

An expedient survey and characterization of the soybean JAGGED 1 (GmJAG1) transcription factor binding preference in the soybean genome by modified ChIPmentation on soybean protoplasts

ChIP-seq is widely used for mapping the transcription factor (TF) binding sites throughout the genome in vivo. In this study, we adopted and modified ChIPmentation, a fast, robust, low-input requirement ChIP-seq method, to a transient expression system using soybean protoplasts to expedite the exploration of TF binding sites. To test this new protocol, we expressed a tagged version of a C2H2-type zinc finger TF, JAGGED1 (GmJAG1), in soybean protoplasts and successfully identified its binding sites in the soybean genome. Furthermore, valuable genomic features such as a novel GmJAG1-binding motif, and the epigenetic characteristics as well as an enhancer-like function of GmJBSs were also found via coupling ATAC-seq and H3K27me3 ChIP-seq data. The application of the modified ChIPmentation protocol in this study using soybean protoplasts provided a new approach for rapid elucidation of how a TF binds to the various target genes in the soybean genome, as illustrated here using GmJAG1.

Determinants of correlated expression of transcription factors and their target genes

While transcription factors (TFs) are known to regulate the expression of their target genes (TGs), only a weak correlation of expression between TFs and their TGs has generally been observed. As lack of correlation could be caused by additional layers of regulation, the overall correlation distribution may hide the presence of a subset of regulatory TF-TG pairs with tight expression coupling. Using reported regulatory pairs in the plant Arabidopsis thaliana along with comprehensive gene expression information and testing a wide array of molecular features, we aimed to discern the molecular determinants of high expression correlation of TFs and their TGs. TF-family assignment, stress-response process involvement, short genomic distances of the TF-binding sites to the transcription start site of their TGs, few required protein-protein-interaction connections to establish physical interactions between the TF and polymerase-II, unambiguous TF-binding motifs, increased numbers of miRNA target-sites in TF-mRNAs, and a young evolutionary age of TGs were found particularly indicative of high TF-TG correlation. The modulating roles of post-transcriptional, post-translational processes, and epigenetic factors have been characterized as well. Our study reveals that regulatory pairs with high expression coupling are associated with specific molecular determinants.

Decoding the plant genome: From epigenome to 3D organization

The linear genome of eukaryotes is partitioned into diverse chromatin states and packaged into a three-dimensional (3D) structure, which has functional implications in DNA replication, DNA repair, and transcriptional regulation. Over the past decades, research on plant functional genomics and epigenomics has made great progress, with thousands of genes cloned and molecular mechanisms of diverse biological processes elucidated. Recently, 3D genome research has gradually attracted great attention of many plant researchers. Herein, we briefly review the progress in genomic and epigenomic research in plants, with a focus on Arabidopsis and rice, and summarize the currently used technologies and advances in plant 3D genome organization studies. We also discuss the relationships between one-dimensional linear genome sequences, epigenomic states, and the 3D chromatin architecture. This review provides basis for future research on plant 3D genomics.

Abiotic stress-mediated modulation of the chromatin landscape in Arabidopsis thaliana

Limited information is available on abiotic stress-mediated alterations of chromatin conformation influencing gene expression in plants. In order to characterize the effect of abiotic stresses on changes in chromatin conformation, we employed FAIRE-seq (formaldehyde-assisted isolation of regulatory element sequencing) and DNase-seq to isolate accessible regions of chromatin from Arabidopsis thaliana seedlings exposed to either heat, cold, salt, or drought stress. Approximately 25% of regions in the Arabidopsis genome were captured as open chromatin, the majority of which included promoters and exons. A large proportion of chromatin regions apparently did not change their conformation in response to any of the four stresses. Digital footprints present within these regions had differential enrichment of motifs for binding of 43 different transcription factors. Further, in contrast to drought and salt stress, both high and low temperature treatments resulted in increased accessibility of the chromatin. Also, pseudogenes attained increased chromatin accessibility in response to cold and drought stresses. The highly accessible and inaccessible chromatin regions of seedlings exposed to drought stress correlated with the Ser/Thr protein kinases (MLK1 and MLK2)-mediated reduction and increase in H3 phosphorylation (H3T3Ph), respectively. The presented results provide a deeper understanding of abiotic stress-mediated chromatin modulation in plants.

Genome-Wide Characterization of DNase I-Hypersensitive Sites and Cold Response Regulatory Landscapes in Grasses

Deep sequencing of DNase-I treated chromatin (DNase-seq) can be used to identify DNase I-hypersensitive sites (DHSs) and facilitates genome-scale mining of de novo cis-regulatory DNA elements. Here, we adapted DNase-seq to generate genome-wide maps of DHSs using control and cold-treated leaf, stem, and root tissues of three widely studied grass species: Brachypodium distachyon, foxtail millet (Setaria italica), and sorghum (Sorghum bicolor). Functional validation demonstrated that 12 of 15 DHSs drove reporter gene expression in transiently transgenic B. distachyon protoplasts. DHSs under both normal and cold treatment substantially differed among tissues and species. Intriguingly, the putative DHS-derived transcription factors (TFs) are largely colocated among tissues and species and include 17 ubiquitous motifs covering all grass taxa and all tissues examined in this study. This feature allowed us to reconstruct a regulatory network that responds to cold stress. Ethylene-responsive TFs SHINE3, ERF2, and ERF9 occurred frequently in cold feedback loops in the tissues examined, pointing to their possible roles in the regulatory network. Overall, we provide experimental annotation of 322,713 DHSs and 93 derived cold-response TF binding motifs in multiple grasses, which could serve as a valuable resource for elucidating the transcriptional networks that function in the cold-stress response and other physiological processes.

Global Role of Crop Genomics in the Face of Climate Change

The development of climate change resilient crops is necessary if we are to meet the challenge of feeding the growing world's population. We must be able to increase food production despite the projected decrease in arable land and unpredictable environmental conditions. This review summarizes the technological and conceptual advances that have the potential to transform plant breeding, help overcome the challenges of climate change, and initiate the next plant breeding revolution. Recent developments in genomics in combination with high-throughput and precision phenotyping facilitate the identification of genes controlling critical agronomic traits. The discovery of these genes can now be paired with genome editing techniques to rapidly develop climate change resilient crops, including plants with better biotic and abiotic stress tolerance and enhanced nutritional value. Utilizing the genetic potential of crop wild relatives (CWRs) enables the domestication of new species and the generation of synthetic polyploids. The high-quality crop plant genome assemblies and annotations provide new, exciting research targets, including long non-coding RNAs (lncRNAs) and cis-regulatory regions. Metagenomic studies give insights into plant-microbiome interactions and guide selection of optimal soils for plant cultivation. Together, all these advances will allow breeders to produce improved, resilient crops in relatively short timeframes meeting the demands of the growing population and changing climate.

3D genome architecture coordinates trans and cis regulation of differentially expressed ear and tassel genes in maize

BACKGROUND

Maize ears and tassels are two separate types of inflorescence which are initiated by similar developmental processes but gradually develop distinct architectures. However, coordinated trans and cis regulation of differentially expressed genes determining ear and tassel architecture within the 3D genome context is largely unknown.

RESULTS

We identify 56,055 and 52,633 open chromatin regions (OCRs) in developing maize ear and tassel primordia using ATAC-seq and characterize combinatorial epigenome features around these OCRs using ChIP-seq, Bisulfite-seq, and RNA-seq datasets. Our integrative analysis of coordinated epigenetic modification and transcription factor binding to OCRs highlights the cis and trans regulation of differentially expressed genes in ear and tassel controlling inflorescence architecture. We further systematically map chromatin interactions at high-resolution in corresponding tissues using in situ digestion-ligation-only Hi-C (DLO Hi-C). The extensive chromatin loops connecting OCRs and genes provide a 3D view on cis- and trans-regulatory modules responsible for ear- and tassel-specific gene expression. We find that intergenic SNPs tend to locate in distal OCRs, and our chromatin interaction maps provide a potential mechanism for trait-associated intergenic SNPs that may contribute to phenotypic variation by influencing target gene expression through chromatin loops.

CONCLUSIONS

Our comprehensive epigenome annotations and 3D genome maps serve as valuable resource and provide a deep understanding of the complex regulatory mechanisms of genes underlying developmental and morphological diversities between maize ear and tassel.

Functional Divergence of the Arabidopsis Florigen-Interacting bZIP Transcription Factors FD and FDP

Flowering of many plant species depends on interactions between basic leucine zipper (bZIP) transcription factors and systemically transported florigen proteins. Members of the genus Arabidopsis contain two of these bZIPs, FD and FDP, which we show have largely complementary expression patterns in shoot apices before and during flowering. CRISPR-Cas9-induced null mutants for FDP flower slightly earlier than wild-type, whereas fd mutants are late flowering. Identical G-box sequences are enriched at FD and FDP binding sites, but only FD binds to genes involved in flowering and only fd alters their transcription. However, both proteins bind to genes involved in responses to the phytohormone abscisic acid (ABA), which controls developmental and stress responses. Many of these genes are differentially expressed in both fd and fdp mutant seedlings, which also show reduced ABA sensitivity. Thus, florigen-interacting bZIPs have distinct functions in flowering dependent on their expression patterns and, at earlier stages in development, play common roles in phytohormone signaling.

Integrative analysis of reference epigenomes in 20 rice varieties

Epigenomic modifications are instrumental for transcriptional regulation, but comprehensive reference epigenomes remain unexplored in rice. Here, we develop an enhanced chromatin immunoprecipitation (eChIP) approach for plants, and generate genome-wide profiling of five histone modifications and RNA polymerase II occupancy with it. By integrating chromatin accessibility, DNA methylation, and transcriptome datasets, we construct comprehensive epigenome landscapes across various tissues in 20 representative rice varieties. Approximately 81.8% of rice genomes are annotated with different epigenomic properties. Refinement of promoter regions using open chromatin and H3K4me3-marked regions provides insight into transcriptional regulation. We identify extensive enhancer-like promoters with potential enhancer function on transcriptional regulation through chromatin interactions. Active and repressive histone modifications and the predicted enhancers vary largely across tissues, whereas inactive chromatin states are relatively stable. Together, these datasets constitute a valuable resource for functional element annotation in rice and indicate the central role of epigenomic information in understanding transcriptional regulation.

Understanding epigenomics based on the rice model

The purpose of this paper provides a comprehensive overview of the recent researches on rice epigenomics, including DNA methylation, histone modifications, noncoding RNAs, and three-dimensional genomics. The challenges and perspectives for future research in rice are discussed. Rice as a model plant for epigenomic studies has much progressed current understanding of epigenetics in plants. Recent results on rice epigenome profiling and three-dimensional chromatin structure studies reveal specific features and implication in gene regulation during rice plant development and adaptation to environmental changes. Results on rice chromatin regulator functions shed light on mechanisms of establishment, recognition, and resetting of epigenomic information in plants. Cloning of several rice epialleles associated with important agronomic traits highlights importance of epigenomic variation in rice plant growth, fitness, and yield. In this review, we summarize and analyze recent advances in rice epigenomics and discuss challenges and directions for future research in the field.

UNBRANCHED3 Expression and Inflorescence Development is Mediated by UNBRANCHED2 and the Distal Enhancer, KRN4, in Maize

Enhancers are cis-acting DNA segments with the ability to increase target gene expression. They show high sensitivity to DNase and contain specific DNA elements in an open chromatin state that allows the binding of transcription factors (TFs). While numerous enhancers are annotated in the maize genome, few have been characterized genetically. KERNEL ROW NUMBER4 (KRN4), an intergenic quantitative trait locus for kernel row number, is assumed to be a cis-regulatory element of UNBRANCHED3 (UB3), a key inflorescence gene. However, the mechanism by which KRN4 controls UB3 expression remains unclear. Here, we found that KRN4 exhibits an open chromatin state, harboring sequences that showed high enhancer activity toward the 35S and UB3 promoters. KRN4 is bound by UB2-centered transcription complexes and interacts with the UB3 promoter by three duplex interactions to affect UB3 expression. Sequence variation at KRN4 enhances ub2 and ub3 mutant ear fasciation. Therefore, we suggest that KRN4 functions as a distal enhancer of the UB3 promoter via chromatin interactions and recruitment of UB2-centered transcription complexes for the fine-tuning of UB3 expression in meristems of ear inflorescences. These results provide evidence that an intergenic region helps to finely tune gene expression, providing a new perspective on the genetic control of quantitative traits.

Interacting Genomic Landscapes of REC8-Cohesin, Chromatin, and Meiotic Recombination in Arabidopsis

Meiosis recombines genetic variation and influences eukaryote genome evolution. During meiosis, DNA double-strand breaks (DSBs) enter interhomolog repair to yield crossovers and noncrossovers. DSB repair occurs as replicated sister chromatids are connected to a polymerized axis. Cohesin rings containing the REC8 kleisin subunit bind sister chromatids and anchor chromosomes to the axis. Here, we report the genomic landscape of REC8 using chromatin immunoprecipitation sequencing (ChIP-seq) in Arabidopsis (Arabidopsis thaliana). REC8 associates with regions of high nucleosome occupancy in multiple chromatin states, including histone methylation at H3K4 (expressed genes), H3K27 (silent genes), and H3K9 (silent transposons). REC8 enrichment is associated with suppression of meiotic DSBs and crossovers at the chromosome and fine scales. As REC8 enrichment is greatest in transposon-dense heterochromatin, we repeated ChIP-seq in kyp suvh5 suvh6 H3K9me2 mutants. Surprisingly, REC8 enrichment is maintained in kyp suvh5 suvh6 heterochromatin and no defects in centromeric cohesion were observed. REC8 occupancy within genes anti-correlates with transcription and is reduced in COPIA transposons that reactivate expression in kyp suvh5 suvh6 Abnormal axis structures form in rec8 that recruit DSB-associated protein foci and undergo synapsis, which is followed by chromosome fragmentation. Therefore, REC8 occupancy correlates with multiple chromatin states and is required to organize meiotic chromosome architecture and interhomolog recombination.

Chromatin domains in space and their functional implications

Genome organization displays functional compartmentalization. Many factors, including epigenetic modifications, transcription factors, chromatin remodelers, and RNAs, shape chromatin domains and the three-dimensional genome organization. Various types of chromatin domains with distinct epigenetic and spatial features exhibit different transcriptional activities. As part of the efforts to better understand plant functional genomics, over the past a few years, spatial distribution patterns of plant chromatin domains have been brought to light. In this review, we discuss chromatin domains associated with the nuclear periphery and the nucleolus, as well as chromatin domains staying in proximity and showing physical interactions. The functional implication of these domains is discussed, with a particular focus on the transcriptional regulation and replication timing. Finally, from a biophysical point of view, we discuss potential roles of liquid-liquid phase separation in plant nuclei in the genesis and maintenance of spatial chromatin domains.

Identification and evolution of gene regulatory networks: insights from comparative studies in plants

The availability of genome sequences, genome-wide assays of transcription factor binding, and accessible chromatin maps have unveiled gene regulatory landscapes in plants. This understanding has ushered in comparative gene regulatory network studies that assess network rewiring between species, across time, and between biological tissues. Comparisons of cis-regulatory elements across the plant kingdom have uncovered examples of conserved sequences, but also of divergence, indicating that selective pressures can vary in different plant families. Transcription factor duplication, followed by spatiotemporal expression divergence of the duplicates, also appears to be a key mechanism of network evolution. Here, we review recent literature describing the regulation of gene expression in plants, and how comparative studies provide insights into how these regulatory interactions change and lead to gene regulatory network rewiring.

Anthocyanin Fruit encodes an R2R3-MYB transcription factor, SlAN2-like, activating the transcription of SlMYBATV to fine-tune anthocyanin content in tomato fruit

Anthocyanin fruit (Aft) and atroviolacea (atv) were characterized in wild tomato and can enhance anthocyanin content in tomato fruit. However, the gene underlying the Aft locus and the mechanism by which Aft and atv act remain largely unknown. In this study, the Aft locus was fine-mapped to an approximately 145-kb interval on chromosome 10, excluding SlAN2 (Solyc10g086250), SlANT1 (Solyc10g086260) and SlANT1-like (Solyc10g086270), which have previously been suggested as candidates. Thus, the R2R3-MYB transcription factor SlAN2-like (Solyc10g086290) was considered the best candidate gene for Aft. The CRISPR/Cas9-mediated SlAN2-like mutants show a much lower accumulation of anthocyanins associated with the downregulation of multiple anthocyanin-related genes compared to the wild-type tomato, indicating that SlAN2-like is responsible for the Aft phenotype. The repressive function of SlMYBATV also was confirmed through the CRISPR/Cas9 approach. A yeast-two-hybrid assay revealed that SlMYBATV interacts with the bHLH protein SlJAF13. Furthermore, yeast-one-hybrid and dual-luciferase transient expression assays showed that Aft directly binds to the SlMYBATV promoter and activates its expression. The results herein provide candidate genes to enhance anthocyanin content in tomato fruit. This research also provides insight into a mechanism involving the Aft-SlMYBATV pathway that fine-tunes anthocyanin accumulation in tomato fruit.

Genome-wide MNase hypersensitivity assay unveils distinct classes of open chromatin associated with H3K27me3 and DNA methylation in Arabidopsis thaliana

BACKGROUND

Regulation of transcription depends on interactions between cis-regulatory elements (CREs) and regulatory proteins. Active CREs are imbedded in open chromatin that are accessible to nucleases. Several techniques, including DNase-seq, which is based on nuclease DNase I, and ATAC-seq, which is based on transposase Tn5, have been widely used to identify genomic regions associated with open chromatin. These techniques have played a key role in dissecting the regulatory networks in gene expression in both animal and plant species.

RESULTS

We develop a technique, named MNase hypersensitivity sequencing (MH-seq), to identify genomic regions associated with open chromatin in Arabidopsis thaliana. Genomic regions enriched with MH-seq reads are referred as MNase hypersensitive sites (MHSs). MHSs overlap with the majority (~ 90%) of the open chromatin identified previously by DNase-seq and ATAC-seq. Surprisingly, 22% MHSs are not covered by DNase-seq or ATAC-seq reads, which are referred to "specific MHSs" (sMHSs). sMHSs tend to be located away from promoters, and a substantial portion of sMHSs are derived from transposable elements. Most interestingly, genomic regions containing sMHSs are enriched with epigenetic marks, including H3K27me3 and DNA methylation. In addition, sMHSs show a number of distinct characteristics including association with transcriptional repressors. Thus, sMHSs span distinct classes of open chromatin that may not be accessible to DNase I or Tn5. We hypothesize that the small size of the MNase enzyme relative to DNase I or Tn5 allows its access to relatively more condensed chromatin domains.

CONCLUSION

MNase can be used to identify open chromatin regions that are not accessible to DNase I or Tn5. Thus, MH-seq provides an important tool to identify and catalog all classes of open chromatin in plants.

Evaluation of genetic diversity, agronomic traits, and anthracnose resistance in the NPGS Sudan Sorghum Core collection

Background

The United States Department of Agriculture (USDA) National Plant Germplasm System (NPGS) sorghum core collection contains 3011 accessions randomly selected from 77 countries. Genomic and phenotypic characterization of this core collection is necessary to encourage and facilitate its utilization in breeding programs and to improve conservation efforts. In this study, we examined the genome sequences of 318 accessions belonging to the NPGS Sudan sorghum core set, and characterized their agronomic traits and anthracnose resistance response.

Results

We identified 183,144 single nucleotide polymorphisms (SNPs) located within or in proximity of 25,124 annotated genes using the genotyping-by-sequencing (GBS) approach. The core collection was genetically highly diverse, with an average pairwise genetic distance of 0.76 among accessions. Population structure and cluster analysis revealed five ancestral populations within the Sudan core set, with moderate to high level of genetic differentiation. In total, 171 accessions (54%) were assigned to one of these populations, which covered 96% of the total genomic variation. Genome scan based on Tajima’s D values revealed two populations under balancing selection. Phenotypic analysis showed differences in agronomic traits among the populations, suggesting that these populations belong to different ecogeographical regions. A total of 55 accessions were resistant to anthracnose; these accessions could represent multiple resistance sources. Genome-wide association study based on fixed and random model Circulating Probability (farmCPU) identified genomic regions associated with plant height, flowering time, panicle length and diameter, and anthracnose resistance response. Integrated analysis of the Sudan core set and sorghum association panel indicated that a large portion of the genetic variation in the Sudan core set might be present in breeding programs but remains unexploited within some clusters of accessions.

Conclusions

The NPGS Sudan core collection comprises genetically and phenotypically diverse germplasm with multiple anthracnose resistance sources. Population genomic analysis could be used to improve screening efforts and identify the most valuable germplasm for breeding programs. The new GBS data set generated in this study represents a novel genomic resource for plant breeders interested in mining the genetic diversity of the NPGS sorghum collection.

The prevalence, evolution and chromatin signatures of plant regulatory elements

Chromatin accessibility and modification is a hallmark of regulatory DNA, the study of which led to the discovery of cis-regulatory elements (CREs). Here, we characterize chromatin accessibility, histone modifications and sequence conservation in 13 plant species. We identified thousands of putative CREs and revealed that distal CREs are prevalent in plants, especially in species with large and complex genomes. The majority of distal CREs have been moved away from their target genes by transposable-element (TE) proliferation, but a substantial number of distal CREs also seem to be created by TEs. Finally, plant distal CREs are associated with three major types of chromatin signatures that are distinct from metazoans. Taken together, these results suggest that CREs are prevalent in plants, highly dynamic during evolution and function through distinct chromatin pathways to regulate gene expression.

Novel mRNAs 3' end-associated cis-regulatory elements with epigenomic signatures of mammalian enhancers in the Arabidopsis genome

The precise spatial and temporal control of gene expression requires the coordinated action of genomic cis-regulatory elements (CREs), including transcriptional enhancers. However, our knowledge of enhancers in plants remains rudimentary and only a few plant enhancers have been experimentally defined. Here, we screened the Arabidopsis thaliana genome and identified >1900 unique candidate CREs that carry the genomic signatures of mammalian enhancers. These were termed putative enhancer-like elements (PEs). Nearly all PEs are intragenic and, unexpectedly, most associate with the 3' ends of protein-coding genes. PEs are hotspots for transcription factor binding and harbor motifs resembling cleavage/polyadenylation signals, potentially coupling 3' end processing to the transcriptional regulation of other genes. Hi-C data showed that 24% of PEs are located at regions that can interact intrachromosomally with other protein-coding genes and, surprisingly, many of these target genes interact with PEs through their 3' UTRs. Examination of the genomes of 1135 sequenced Arabidopsis accessions showed that PEs are conserved. Our findings suggest that the identified PEs may serve as transcriptional enhancers and sites for mRNA 3' end processing, and constitute a novel group of CREs in Arabidopsis.

BORDER proteins protect expression of neighboring genes by promoting 3' Pol II pausing in plants

Ensuring that one gene's transcription does not inappropriately affect the expression of its neighbors is a fundamental challenge to gene regulation in a genomic context. In plants, which lack homologs of animal insulator proteins, the mechanisms that prevent transcriptional interference are not well understood. Here we show that BORDER proteins are enriched in intergenic regions and prevent interference between closely spaced genes on the same strand by promoting the 3' pausing of RNA polymerase II at the upstream gene. In the absence of BORDER proteins, 3' pausing associated with the upstream gene is reduced and shifts into the promoter region of the downstream gene. This is consistent with a model in which BORDER proteins inhibit transcriptional interference by preventing RNA polymerase from intruding into the promoters of downstream genes.

The effects of Arabidopsis genome duplication on the chromatin organization and transcriptional regulation

Autopolyploidy is widespread in higher plants and important for agricultural yield and quality. However, the effects of genome duplication on the chromatin organization and transcriptional regulation are largely unknown in plants. Using High-throughput Chromosome Conformation Capture (Hi-C), we showed that autotetraploid Arabidopsis presented more inter-chromosomal interactions and fewer short-range chromatin interactions compared with its diploid progenitor. In addition, genome duplication contributed to the switching of some loose and compact structure domains with altered H3K4me3 and H3K27me3 histone modification status. 539 genes were identified with altered transcriptions and chromatin interactions in autotetraploid Arabidopsis. Especially, we found that genome duplication changed chromatin looping and H3K27me3 histone modification in Flowering Locus C. We propose that genome doubling modulates the transcription genome-wide by changed chromatin interactions and at the specific locus by altered chromatin loops and histone modifications.

Identification and functional evaluation of accessible chromatin associated with wood formation in Eucalyptus grandis

Accessible chromatin changes dynamically during development and harbours functional regulatory regions which are poorly understood in the context of wood development. We explored the importance of accessible chromatin in Eucalyptus grandis in immature xylem generally, and MYB transcription factor-mediated transcriptional programmes specifically. We identified biologically reproducible DNase I Hypersensitive Sites (DHSs) and assessed their functional significance in immature xylem through their associations with gene expression, epigenomic data and DNA sequence conservation. We identified in vitro DNA binding sites for six secondary cell wall-associated Eucalyptus MYB (EgrMYB) transcription factors using DAP-seq, reconstructed protein-DNA networks of predicted targets based on binding sites within or outside DHSs and assessed biological enrichment of these networks with published datasets. 25 319 identified immature xylem DHSs were associated with increased transcription and significantly enriched for various epigenetic signatures (H3K4me3, H3K27me3, RNA pol II), conserved noncoding sequences and depleted single nucleotide variants. Predicted networks built from EgrMYB binding sites located in accessible chromatin were significantly enriched for systems biology datasets relevant to wood formation, whereas those occurring in inaccessible chromatin were not. Our study demonstrates that DHSs in E. grandis immature xylem, most of which are intergenic, are of functional significance to gene regulation in this tissue.

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.

Chromatin Signature and Transcription Factor Binding Provide a Predictive Basis for Understanding Plant Gene Expression

Chromatin accessibility and post-transcriptional histone modifications play important roles in gene expression regulation. However, little is known about the joint effect of multiple chromatin modifications on the gene expression level in plants, despite that the regulatory roles of individual histone marks such as H3K4me3 in gene expression have been well-documented. By using machine-learning methods, we systematically performed gene expression level prediction based on multiple chromatin modifications data in Arabidopsis and rice. We found that as few as four histone modifications were sufficient to yield good prediction performance, and H3K4me3 and H3K36me3 being the top two predictors with known functions related to transcriptional initiation and elongation, respectively. We demonstrated that the predictive powers differed between protein-coding and non-coding genes as well as between CpG-enriched and CpG-depleted genes. We also showed that the predictive model trained in one tissue or species could be applied to another tissue or species, suggesting shared underlying mechanisms. More interestingly, the gene expression levels of conserved orthologs are easier to predict than the species-specific genes. In addition, chromatin state of distal enhancers was moderately correlated to gene expression but was dispensable if given the chromatin features of the proximal regions of genes. We further extended the analysis to transcription factor (TF) binding data. Strikingly, the combinatorial effects of only a few TFs were roughly fit to gene expression levels in Arabidopsis. Overall, by using quantitative modeling, we provide a comprehensive and unbiased perspective on the epigenetic and TF-mediated regulation of gene expression in plants.

Cold stress induces enhanced chromatin accessibility and bivalent histone modifications H3K4me3 and H3K27me3 of active genes in potato

BACKGROUND

Cold stress can greatly affect plant growth and development. Plants have developed special systems to respond to and tolerate cold stress. While plant scientists have discovered numerous genes involved in responses to cold stress, few studies have been dedicated to investigation of genome-wide chromatin dynamics induced by cold or other abiotic stresses.

RESULTS

Genomic regions containing active cis-regulatory DNA elements can be identified as DNase I hypersensitive sites (DHSs). We develop high-resolution DHS maps in potato (Solanum tuberosum) using chromatin isolated from tubers stored under room (22 °C) and cold (4 °C) conditions. We find that cold stress induces a large number of DHSs enriched in genic regions which are frequently associated with differential gene expression in response to temperature variation. Surprisingly, active genes show enhanced chromatin accessibility upon cold stress. A large number of active genes in cold-stored tubers are associated with the bivalent H3K4me3-H3K27me3 mark in gene body regions. Interestingly, upregulated genes associated with the bivalent mark are involved in stress response, whereas downregulated genes with the bivalent mark are involved in developmental processes. In addition, we observe that the bivalent mark-associated genes are more accessible than others upon cold stress.

CONCLUSIONS

Collectively, our results suggest that cold stress induces enhanced chromatin accessibility and bivalent histone modifications of active genes. We hypothesize that in cold-stored tubers, the bivalent H3K4me3-H3K27me3 mark represents a distinct chromatin environment with greater accessibility, which may facilitate the access of regulatory proteins required for gene upregulation or downregulation in response to cold stress.

Global Quantitative Mapping of Enhancers in Rice by STARR-seq

Enhancers activate transcription in a distance-, orientation-, and position-independent manner, which makes them difficult to be identified. Self-transcribing active regulatory region sequencing (STARR-seq) measures the enhancer activity of millions of DNA fragments in parallel. Here we used STARR-seq to generate a quantitative global map of rice enhancers. Most enhancers were mapped within genes, especially at the 5′ untranslated regions (5′UTR) and in coding sequences. Enhancers were also frequently mapped proximal to silent and lowly-expressed genes in transposable element (TE)-rich regions. Analysis of the epigenetic features of enhancers at their endogenous loci revealed that most enhancers do not co-localize with DNase I hypersensitive sites (DHSs) and lack the enhancer mark of histone modification H3K4me1. Clustering analysis of enhancers according to their epigenetic marks revealed that about 40% of identified enhancers carried one or more epigenetic marks. Repressive H3K27me3 was frequently enriched with positive marks, H3K4me3 and/or H3K27ac, which together label enhancers. Intergenic enhancers were also predicted based on the location of DHS regions relative to genes, which overlap poorly with STARR-seq enhancers. In summary, we quantitatively identified enhancers by functional analysis in the genome of rice, an important model plant. This work provides a valuable resource for further mechanistic studies in different biological contexts.

MED25 connects enhancer-promoter looping and MYC2-dependent activation of jasmonate signalling

The lipid-derived hormone jasmonate (JA) regulates plant immunity and adaptive growth by triggering a genome-wide transcriptional programme. In Arabidopsis thaliana, JA-triggered transcriptional programming is largely orchestrated by the master transcription factor MYC2. The function of MYC2 is dependent on its physical interaction with the MED25 subunit of the Mediator transcriptional co-activator complex. Here we report the identification of JA enhancers (JAEs) through profiling the occupancy pattern of MYC2 and MED25. JA regulates the dynamic chromatin looping between JAEs and their promoters in a MED25-dependent manner, while MYC2 auto-regulates itself through JAEs. Interestingly, the JAE of the MYC2 locus (named ME2) positively regulates MYC2 expression during short-term JA responses but negatively regulates it during constant JA responses. We demonstrate that new gene editing tools open up new avenues to elucidate the in vivo function of enhancers. Our work provides a paradigm for functional study of plant enhancers in the regulation of specific physiological processes.

MED25 connects enhancer-promoter looping and MYC2-dependent activation of jasmonate signalling

The lipid-derived hormone jasmonate (JA) regulates plant immunity and adaptive growth by triggering a genome-wide transcriptional programme. In Arabidopsis thaliana, JA-triggered transcriptional programming is largely orchestrated by the master transcription factor MYC2. The function of MYC2 is dependent on its physical interaction with the MED25 subunit of the Mediator transcriptional co-activator complex. Here we report the identification of JA enhancers (JAEs) through profiling the occupancy pattern of MYC2 and MED25. JA regulates the dynamic chromatin looping between JAEs and their promoters in a MED25-dependent manner, while MYC2 auto-regulates itself through JAEs. Interestingly, the JAE of the MYC2 locus (named ME2) positively regulates MYC2 expression during short-term JA responses but negatively regulates it during constant JA responses. We demonstrate that new gene editing tools open up new avenues to elucidate the in vivo function of enhancers. Our work provides a paradigm for functional study of plant enhancers in the regulation of specific physiological processes.

CRISPR/Cas Genome Editing and Precision Plant Breeding in Agriculture

Enhanced agricultural production through innovative breeding technology is urgently needed to increase access to nutritious foods worldwide. Recent advances in CRISPR/Cas genome editing enable efficient targeted modification in most crops, thus promising to accelerate crop improvement. Here, we review advances in CRISPR/Cas9 and its variants and examine their applications in plant genome editing and related manipulations. We highlight base-editing tools that enable targeted nucleotide substitutions and describe the various delivery systems, particularly DNA-free methods, that have linked genome editing with crop breeding. We summarize the applications of genome editing for trait improvement, development of techniques for fine-tuning gene regulation, strategies for breeding virus resistance, and the use of high-throughput mutant libraries. We outline future perspectives for genome editing in plant synthetic biology and domestication, advances in delivery systems, editing specificity, homology-directed repair, and gene drives. Finally, we discuss the challenges and opportunities for precision plant breeding and its bright future in agriculture.

Dynamic control of enhancer activity drives stage-specific gene expression during flower morphogenesis

Enhancers are critical for developmental stage-specific gene expression, but their dynamic regulation in plants remains poorly understood. Here we compare genome-wide localization of H3K27ac, chromatin accessibility and transcriptomic changes during flower development in Arabidopsis. H3K27ac prevalently marks promoter-proximal regions, suggesting that H3K27ac is not a hallmark for enhancers in Arabidopsis. We provide computational and experimental evidence to confirm that distal DNase І hypersensitive sites are predictive of enhancers. The predicted enhancers are highly stage-specific across flower development, significantly associated with SNPs for flowering-related phenotypes, and conserved across crucifer species. Through the integration of genome-wide transcription factor (TF) binding datasets, we find that floral master regulators and stage-specific TFs are largely enriched at developmentally dynamic enhancers. Finally, we show that enhancer clusters and intronic enhancers significantly associate with stage-specific gene regulation by floral master TFs. Our study provides insights into the functional flexibility of enhancers during plant development, as well as hints to annotate plant enhancers.

Enhancer elements can control spatial and temporal patterns of gene expression. Here the authors profile chromatin accessibility, histone modifications and gene expression during Arabidopsis flower development providing evidence for sets of distal enhancers acting in a highly stage-specific manner.

Landscape of gene transposition-duplication within the Brassicaceae family

We developed the CLfinder-OrthNet pipeline that detects co-linearity among multiple closely related genomes, finds orthologous gene groups, and encodes the evolutionary history of each orthologue group into a representative network (OrthNet). Using a search based on network topology, we identified 1,394 OrthNets that included gene transposition-duplication (tr-d) events, out of 17,432 identified in six Brassicaceae genomes. Occurrences of tr-d shared by subsets of Brassicaceae genomes mirrored the divergence times between the genomes and their repeat contents. The majority of tr-d events resulted in truncated open reading frames (ORFs) in the duplicated loci. However, the duplicates with complete ORFs were significantly more frequent than expected from random events. These were derived from older tr-d events and had a higher chance of being expressed. We also found an enrichment of tr-d events with complete loss of intergenic sequence conservation between the original and duplicated loci. Finally, we identified tr-d events uniquely found in two extremophytes among the six Brassicaceae genomes, including tr-d of SALT TOLERANCE 32 and ZINC TRANSPORTER 3 that relate to their adaptive evolution. CLfinder-OrthNet provides a flexible toolkit to compare gene order, visualize evolutionary paths among orthologues as networks, and identify gene loci that share an evolutionary history.

Rapid validation of transcriptional enhancers using agrobacterium-mediated transient assay

BACKGROUND

Enhancers are one of the most important classes of cis-regulatory elements (CREs) and play key roles in regulation of transcription in higher eukaryotes. Enhancers are difficult to identify because they lack positional constraints relative to their cognate genes. Excitingly, several recent studies showed that plant enhancers can be predicted based on their distinct features associated with open chromatin. However, experimental validation is necessary to confirm the predicted enhancer function.

RESULTS

We developed a rapid enhancer validation system based on Nicotiana benthamiana. A set of 12 intergenic and intronic enhancers, identified in Arabidopsis thaliana, were cloned into a vector containing a minimal 35S promoter and a luciferase reporter gene, and were then infiltrated into N. benthamiana leaves mediated by agrobacterium. The enhancer activity of each candidate was quantitatively assayed based on bioluminescence measurement. The data from this luciferase-based validation was correlated with previous data derived from transgenic assays in A. thaliana. In addition, the relative strength of different enhancers for driving the reporter gene can be quantitatively compared. We demonstrate that this system can also be used to map the functional activity of a candidate enhancer under different environmental conditions.

CONCLUSIONS

In summary, we developed a rapid and efficient plant enhancer validation system based on a luciferase reporter and N. benthamiana-based leaf agroinfiltration. This system can be used for initial screening of leaf-specific enhancers and for validating candidate leaf enhancers from dicot species. It can potentially be used to examine the activity of candidate enhancers under different environmental conditions.

Bidirectional Promoter-Based CRISPR-Cas9 Systems for Plant Genome Editing

CRISPR-Cas systems can be expressed in multiple ways, with different capabilities regarding tissue-specific expression, efficiency, and expression levels. Thus far, three expression strategies have been demonstrated in plants: mixed dual promoter systems, dual Pol II promoter systems, and single transcript unit (STU) systems. We explored a fourth strategy to express CRISPR-Cas9 in the model and crop plant, rice, where a bidirectional promoter (BiP) is used to express Cas9 and single guide RNA (sgRNA) in opposite directions. We first tested an engineered BiP system based on double-mini 35S promoter and an Arabidopsis enhancer, which resulted in 20.7% and 52.9% genome editing efficiencies at two target sites in T0 stable transgenic rice plants. We further improved the BiP system drastically by using a rice endogenous BiP, OsBiP1. The endogenous BiP expression system had higher expression strength and led to 75.9-93.3% genome editing efficiencies in rice T0 generation, when the sgRNAs were processed by either tRNA or Csy4. We provided a proof-of-concept study of applying BiP systems for expressing two-component CRISPR-Cas9 genome editing reagents in rice. Our work could promote future research and adoption of BiP systems for CRISPR-Cas-based genome engineering in plants.

Alternative Splicing and Protein Diversity: Plants Versus Animals

Plants, unlike animals, exhibit a very high degree of plasticity in their growth and development and employ diverse strategies to cope with the variations during diurnal cycles and stressful conditions. Plants and animals, despite their remarkable morphological and physiological differences, share many basic cellular processes and regulatory mechanisms. Alternative splicing (AS) is one such gene regulatory mechanism that modulates gene expression in multiple ways. It is now well established that AS is prevalent in all multicellular eukaryotes including plants and humans. Emerging evidence indicates that in plants, as in animals, transcription and splicing are coupled. Here, we reviewed recent evidence in support of co-transcriptional splicing in plants and highlighted similarities and differences between plants and humans. An unsettled question in the field of AS is the extent to which splice isoforms contribute to protein diversity. To take a critical look at this question, we presented a comprehensive summary of the current status of research in this area in both plants and humans, discussed limitations with the currently used approaches and suggested improvements to current methods and alternative approaches. We end with a discussion on the potential role of epigenetic modifications and chromatin state in splicing memory in plants primed with stresses.

Enhancer-Promoter Interaction of SELF PRUNING 5G Shapes Photoperiod Adaptation

Tomato (Solanum lycopersicum) is a major vegetable fruit grown and consumed worldwide. Modern cultivated tomatoes are derived from their wild relative, Solanum pimpinellifolium, a short-day plant that originated from the Andean region of South America. The molecular underpinnings of the regional adaptation and expansion of domesticated tomato remain largely unclear. In this study, we examined flowering time in wild and cultivated tomatoes under both long-day and short-day conditions. Using quantitative trait locus mapping in a recombinant inbred line population, we identified SELF PRUNING 5G (SP5G) as a major locus influencing daylength adaptation in tomato. Genetic diversity analysis revealed that the genomic region harboring SP5G shows signatures of a domestication sweep. We found that a 52-bp sequence within the 3' untranslated region of SP5G is essential for the enhanced expression of this gene, leading to delayed flowering time in tomatoes through a promoter-enhancer interaction that occurs only under long-day conditions. We further demonstrate that the absence of the 52-bp sequence attenuates the promoter-enhancer interaction and reduces SP5G expression in cultivated tomatoes, making their flowering time insensitive to daylength. Our findings demonstrate that cis-regulatory variation at the enhancer region of the SP5G 3' untranslated region confers reduced photoperiodic response in cultivated tomatoes, uncovering a regulatory mechanism that could potentially be used to manipulate flowering time in tomato through novel biotechnological approaches.

A Trihelix Family Transcription Factor Is Associated with Key Genes in Mixed-Linkage Glucan Accumulation

Mixed-linkage glucan (MLG) is a polysaccharide that is highly abundant in grass endosperm cell walls and present at lower amounts in other tissues. Cellulose synthase-like F (CSLF) and cellulose synthase-like H genes synthesize MLG, but it is unknown if other genes participate in the production and restructuring of MLG. Using Brachypodium distachyon transcriptional profiling data, we identified a B distachyon trihelix family transcription factor (BdTHX1) that is highly coexpressed with the B distachyon CSLF6 gene (BdCSLF6), which suggests that BdTHX1 is involved in the regulation of MLG biosynthesis. To determine the genes regulated by this transcription factor, we conducted chromatin immunoprecipitation sequencing (ChIP-seq) experiments using immature B distachyon seeds and an anti-BdTHX1 polyclonal antibody. The ChIP-seq experiment identified the second intron of BdCSLF6 as one of the most enriched sequences. The binding of BdTHX1 to the BdCSLF6 intron sequence was confirmed using electrophoretic mobility shift assays (EMSA). ChIP-seq also showed that a gene encoding a grass-specific glycoside hydrolase family 16 endotransglucosylase/hydrolase (BdXTH8) is bound by BdTHX1, and the binding was confirmed by EMSA. Radiochemical transglucanase assays showed that BdXTH8 exhibits predominantly MLG:xyloglucan endotransglucosylase activity, a hetero-transglycosylation reaction, and can thus produce MLG-xyloglucan covalent bonds; it also has a lower xyloglucan:xyloglucan endotransglucosylase activity. B distachyon shoots regenerated from transformed calli overexpressing BdTHX1 showed an abnormal arrangement of vascular tissue and seedling-lethal phenotypes. These results indicate that the transcription factor BdTHX1 likely plays an important role in MLG biosynthesis and restructuring by regulating the expression of BdCSLF6 and BdXTH8.

Principles and characteristics of the Arabidopsis WRKY regulatory network during early MAMP-triggered immunity

During microbe-associated molecular pattern-triggered immunity more than 5000 Arabidopsis genes are significantly altered in their expression, and the question arises, how such an enormous reprogramming of the transcriptome can be regulated in a safe and robust manner? For the WRKY transcription factors (TFs), which are important regulators of numerous defense responses, it appears that they act in a complex regulatory sub-network rather than in a linear fashion, which would be much more vulnerable to gene function loss either by pathogen-derived effectors or by mutations. In this study we employed RNA-seq, mass spectrometry and chromatin immunoprecipitation-seq to find evidence for and uncover principles and characteristics of this network. Upon flg22-treatment, one can distinguish between two sets of WRKY genes: constitutively expressed and induced WRKY genes. Prior to elicitation the induced WRKY genes appear to be maintained in a repressed state mainly by the constitutively expressed WRKY factors, which themselves appear to be regulated by non-WRKY TFs. Upon elicitation, induced WRKYs rapidly bind to induced WRKY gene promoters and by auto- and cross-regulation build up the regulatory network. Maintenance of this flg22-induced network appears highly robust as removal of three key WRKY factors can be physically and functionally compensated for by other WRKY family members.

Proliferation of Regulatory DNA Elements Derived from Transposable Elements in the Maize Genome

Genomic regions free of nucleosomes, which are hypersensitive to DNase I digestion, are known as DNase I hypersensitive sites (DHSs) and frequently contain cis-regulatory DNA elements. To investigate their prevalence and characteristics in maize (Zea mays), we developed high-resolution genome-wide DHS maps using a modified DNase-seq technique. Maize DHSs exhibit depletion of nucleosomes and low levels of DNA methylation and are enriched with conserved noncoding sequences (CNSs). We developed a protoplast-based transient transformation assay to assess the potential gene expression enhancer and/or promoter functions associated with DHSs, which showed that more than 80% of DHSs overlapping with CNSs showed an enhancer function. Strikingly, nearly 25% of maize DHSs were derived from transposable elements (TEs), including both class I and class II transposons. Interestingly, TE-derived DHSs (teDHSs) homologous to retrotransposons were enriched with sequences related to the intrinsic cis-regulatory elements within the long terminal repeats of retrotransposons. We demonstrate that more than 80% of teDHSs can drive transcription of a reporter gene in protoplast assays. These results reveal the widespread occurrence of TE-derived cis-regulatory sequences and suggest that teDHSs play a major role in transcriptional regulation in maize.

Identification of cis-regulatory elements by chromatin structure

The systematic identification of cis-regulatory elements (CREs) in plant genomes is critically important in understanding transcriptional regulation during development and in response to environmental cues. Several genome-wide structure-based methods have been successfully applied to plant genomes in the past few years. Here, we review recent results on the identification and characterization of CREs in multiple plant species and in different biological processes and discuss future applications of chromatin accessibility data to understand the mechanism, function and evolution of transcriptional regulation networks.

ChIP-ping the branches of the tree: functional genomics and the evolution of eukaryotic gene regulation

Advances in the methods for detecting protein-DNA interactions have played a key role in determining the directions of research into the mechanisms of transcriptional regulation. The most recent major technological transformation happened a decade ago, with the move from using tiling arrays [chromatin immunoprecipitation (ChIP)-on-Chip] to high-throughput sequencing (ChIP-seq) as a readout for ChIP assays. In addition to the numerous other ways in which it is superior to arrays, by eliminating the need to design and manufacture them, sequencing also opened the door to carrying out comparative analyses of genome-wide transcription factor occupancy across species and studying chromatin biology in previously less accessible model and nonmodel organisms, thus allowing us to understand the evolution and diversity of regulatory mechanisms in unprecedented detail. Here, we review the biological insights obtained from such studies in recent years and discuss anticipated future developments in the field.

Profiling of Accessible Chromatin Regions across Multiple Plant Species and Cell Types Reveals Common Gene Regulatory Principles and New Control Modules

The transcriptional regulatory structure of plant genomes remains poorly defined relative to animals. It is unclear how many cis-regulatory elements exist, where these elements lie relative to promoters, and how these features are conserved across plant species. We employed the assay for transposase-accessible chromatin (ATAC-seq) in four plant species (Arabidopsis thaliana, Medicago truncatula, Solanum lycopersicum, and Oryza sativa) to delineate open chromatin regions and transcription factor (TF) binding sites across each genome. Despite 10-fold variation in intergenic space among species, the majority of open chromatin regions lie within 3 kb upstream of a transcription start site in all species. We find a common set of four TFs that appear to regulate conserved gene sets in the root tips of all four species, suggesting that TF-gene networks are generally conserved. Comparative ATAC-seq profiling of Arabidopsis root hair and non-hair cell types revealed extensive similarity as well as many cell-type-specific differences. Analyzing TF binding sites in differentially accessible regions identified a MYB-driven regulatory module unique to the hair cell, which appears to control both cell fate regulators and abiotic stress responses. Our analyses revealed common regulatory principles among species and shed light on the mechanisms producing cell-type-specific transcriptomes during development.

Profiling of Accessible Chromatin Regions across Multiple Plant Species and Cell Types Reveals Common Gene Regulatory Principles and New Control Modules

The transcriptional regulatory structure of plant genomes remains poorly defined relative to animals. It is unclear how many cis-regulatory elements exist, where these elements lie relative to promoters, and how these features are conserved across plant species. We employed the assay for transposase-accessible chromatin (ATAC-seq) in four plant species (Arabidopsis thaliana, Medicago truncatula, Solanum lycopersicum, and Oryza sativa) to delineate open chromatin regions and transcription factor (TF) binding sites across each genome. Despite 10-fold variation in intergenic space among species, the majority of open chromatin regions lie within 3 kb upstream of a transcription start site in all species. We find a common set of four TFs that appear to regulate conserved gene sets in the root tips of all four species, suggesting that TF-gene networks are generally conserved. Comparative ATAC-seq profiling of Arabidopsis root hair and non-hair cell types revealed extensive similarity as well as many cell-type-specific differences. Analyzing TF binding sites in differentially accessible regions identified a MYB-driven regulatory module unique to the hair cell, which appears to control both cell fate regulators and abiotic stress responses. Our analyses revealed common regulatory principles among species and shed light on the mechanisms producing cell-type-specific transcriptomes during development.

Chromatin accessibility dynamics reveal novel functional enhancers in C. elegans

Chromatin accessibility, a crucial component of genome regulation, has primarily been studied in homogeneous and simple systems, such as isolated cell populations or early-development models. Whether chromatin accessibility can be assessed in complex, dynamic systems in vivo with high sensitivity remains largely unexplored. In this study, we use ATAC-seq to identify chromatin accessibility changes in a whole animal, the model organism Caenorhabditis elegans, from embryogenesis to adulthood. Chromatin accessibility changes between developmental stages are highly reproducible, recapitulate histone modification changes, and reveal key regulatory aspects of the epigenomic landscape throughout organismal development. We find that over 5000 distal noncoding regions exhibit dynamic changes in chromatin accessibility between developmental stages and could thereby represent putative enhancers. When tested in vivo, several of these putative enhancers indeed drive novel cell-type- and temporal-specific patterns of expression. Finally, by integrating transcription factor binding motifs in a machine learning framework, we identify EOR-1 as a unique transcription factor that may regulate chromatin dynamics during development. Our study provides a unique resource for C. elegans, a system in which the prevalence and importance of enhancers remains poorly characterized, and demonstrates the power of using whole organism chromatin accessibility to identify novel regulatory regions in complex systems.