Stable unmethylated DNA demarcates expressed genes and their cis-regulatory space in plant genomes

Potent pollen gene regulation by DNA glycosylases in maize

Although DNA methylation primarily represses TEs, it also represses select genes that are methylated in plant body tissues but demethylated by DNA glycosylases (DNGs) in endosperm or pollen. Activity of either one of two DNGs, MDR1 or DNG102, is essential for pollen viability in maize. Using single-pollen mRNA sequencing on pollen segregating mutations in both genes, we identified 58 candidate DNG target genes that account for 11.1% of the wild-type transcriptome but are silent or barely detectable in the plant body (sporophyte). They are unusual in their tendency to lack introns but even more so in their having TE-like methylation in their CDS. The majority have predicted functions in cell wall modification, and they likely support the rapid tip growth characteristic of pollen tubes. These results suggest a critical role for DNA methylation and demethylation in regulating maize genes with potential for extremely high expression in pollen but constitutive silencing elsewhere.

Computational Reconstruction of the Transcription Factor Regulatory Network Induced by Auxin in Arabidopsis thaliana L

In plant hormone signaling, transcription factor regulatory networks (TFRNs), which link the master transcription factors to the biological processes under their control, remain insufficiently characterized despite their crucial function. Here, we identify a TFRN involved in the response to the key plant hormone auxin and define its impact on auxin-driven biological processes. To reconstruct the TFRN, we developed a three-step procedure, which is based on the integrated analysis of differentially expressed gene lists and a representative collection of transcription factor binding profiles. Its implementation is available as a part of the CisCross web server. With the new method, we distinguished two transcription factor subnetworks. The first operates before auxin treatment and is switched off upon hormone application, the second is switched on by the hormone. Moreover, we characterized the functioning of the auxin-regulated TFRN in control of chlorophyll and lignin biosynthesis, abscisic acid signaling, and ribosome biogenesis.

Genome-wide chromatin accessibility reveals transcriptional regulation of heterosis in inter-subspecific hybrid rice

The utilization of rice heterosis is essential for ensuring global food security; however, its molecular mechanism remains unclear. In this study, comprehensive analyses of accessible chromatin regions (ACRs), DNA methylation, and gene expression in inter-subspecific hybrid and its parents were performed to determine the potential role of chromatin accessibility in rice heterosis. The hybrid exhibited abundant ACRs, in which the gene ACRs and proximal ACRs were directly related to transcriptional activation rather than the distal ACRs. Regarding the dynamic accessibility contribution of the parents, paternal ZHF1015 transmitted a greater number of ACRs to the hybrid. Accessible genotype-specific target genes were enriched with overrepresented transcription factors, indicating a unique regulatory network of genes in the hybrid. Compared with its parents, the differentially accessible chromatin regions with upregulated chromatin accessibility were much greater than those with downregulated chromatin accessibility, reflecting a stronger regulation in the hybrid. Furthermore, DNA methylation levels were negatively correlated with ACR intensity, and genes were strongly affected by CHH methylation in the hybrid. Chromatin accessibility positively regulated the overall expression level of each genotype. ACR-related genes with maternal Z04A-bias allele-specific expression tended to be enriched during carotenoid biosynthesis, whereas paternal ZHF1015-bias genes were more active in carbohydrate metabolism. Our findings provide a new perspective on the mechanism of heterosis based on chromatin accessibility in inter-subspecific hybrid rice.

DNA methylation variations underlie lettuce domestication and divergence

BACKGROUND

Lettuce (Lactuca sativa L.) is an economically important vegetable crop worldwide. Lettuce is believed to be domesticated from a single wild ancestor Lactuca serriola and subsequently diverged into two major morphologically distinct vegetable types: leafy lettuce and stem lettuce. However, the role of epigenetic variation in lettuce domestication and divergence remains largely unknown.

RESULTS

To understand the genetic and epigenetic basis underlying lettuce domestication and divergence, we generate single-base resolution DNA methylomes from 52 Lactuca accessions, including major lettuce cultivars and wild relatives. We find a significant increase of DNA methylation during lettuce domestication and uncover abundant epigenetic variations associated with lettuce domestication and divergence. Interestingly, DNA methylation variations specifically associated with leafy and stem lettuce are related to regulation and metabolic processes, respectively, while those associated with both types are enriched in stress responses. Moreover, we reveal that domestication-induced DNA methylation changes could influence expression levels of nearby and distal genes possibly through affecting chromatin accessibility and chromatin loop.

CONCLUSION

Our study provides population epigenomic insights into crop domestication and divergence and valuable resources for further domestication for diversity and epigenetic breeding to boost crop improvement.

Transcription factor binding site divergence across maize inbred lines drives transcriptional and phenotypic variation

Regulatory elements are important constituents of plant genomes that have shaped ancient and modern crops. Their identification, function, and diversity in crop genomes however are poorly characterized, thus limiting our ability to harness their power for further agricultural advances using induced or natural variation. Here, we use DNA affinity purification-sequencing (DAP-seq) to map transcription factor (TF) binding events for 200 maize TFs belonging to 30 distinct families and heterodimer pairs in two distinct inbred lines historically used for maize hybrid plant production, providing empirical binding site annotation for 5.3% of the maize genome. TF binding site comparison in B73 and Mo17 inbreds reveals widespread differences, driven largely by structural variation, that correlate with gene expression changes. TF binding site presence-absence variation helps clarify complex QTL such as vgt1, an important determinant of maize flowering time, and DICE, a distal enhancer involved in herbivore resistance. Modification of TF binding regions via CRISPR-Cas9 mediated editing alters target gene expression and phenotype. Our functional catalog of maize TF binding events enables collective and comparative TF binding analysis, and highlights its value for agricultural improvement.

Establishing an optimized ATAC-seq protocol for the maize

The advent of next-generation sequencing in crop improvement offers unprecedented insights into the chromatin landscape closely linked to gene activity governing key traits in plant development and adaptation. Particularly in maize, its dynamic chromatin structure is found to collaborate with massive transcriptional variations across tissues and developmental stages, implying intricate regulatory mechanisms, which highlights the importance of integrating chromatin information into breeding strategies for precise gene controls. The depiction of maize chromatin architecture using Assay for Transposase Accessible Chromatin with high-throughput sequencing (ATAC-seq) provides great opportunities to investigate cis-regulatory elements, which is crucial for crop improvement. In this context, we developed an easy-to-implement ATAC-seq protocol for maize with fewer nuclei and simple equipment. We demonstrate a streamlined ATAC-seq protocol with four key steps for maize in which nuclei purification can be achieved without cell sorting and using only a standard bench-top centrifuge. Our protocol, coupled with the bioinformatic analysis, including validation by read length periodicity, key metrics, and correlation with transcript abundance, provides a precise and efficient assessment of the maize chromatin landscape. Beyond its application to maize, our testing design holds the potential to be applied to other crops or other tissues, especially for those with limited size and amount, establishing a robust foundation for chromatin structure studies in diverse crop species.

Three-year study of DNA cytosine methylation dynamics in transplanted Malbec grapevines

DNA cytosine methylation, an epigenetic mechanism involved in gene regulation and genome stability, remains poorly understood in terms of its role under changing environmental conditions. Previous research using methylation-sensitive amplified polymorphism (MSAP) markers in a Vitis vinifera L. cv. Malbec clone showed vineyard-specific DNA methylation polymorphism, but no change in overall methylation levels. To complement these findings, the present study investigates the intra-seasonal epigenetic dynamics between genetically identical plants grown in different vineyards through a transplanting experiment. Cuttings of the same clone, showing differential methylation patterns imposed by the vineyard of origin (Agrelo and Gualtallary), were cultivated in a common vineyard (Lunlunta). Using high-performance liquid chromatography-ultraviolet detection, the quantification of global DNA 5-methylcytosine (5-mC) levels revealed relatively low overall 5-mC percentages in grapevines, with higher levels in Agrelo (5.8%) compared to Gualtallary plants (3.7%). The transplanted plants maintained the 5-mC levels differences between vineyards (9.8% vs 6.2%), which equalized in subsequent seasons (7.5% vs 7%). Additionally, the study examined 5-mC polymorphism using MSAP markers in Lunlunta transplanted plants over three seasons. The observed differences between vineyards in MSAP patterns during the initial growing season gradually diminished, suggesting a reprogramming of the hemimethylated pattern following implantation in the common vineyard. In contrast, the non-methylated pattern exhibited greater stability, indicating a potential memory effect. Overall, this study provides valuable insights into the dynamic nature of DNA methylation in grapevines under changing environmental conditions, with potential implications for crop management and breeding strategies.

Extreme restructuring of cis-regulatory regions controlling a deeply conserved plant stem cell regulator

A striking paradox is that genes with conserved protein sequence, function and expression pattern over deep time often exhibit extremely divergent cis-regulatory sequences. It remains unclear how such drastic cis-regulatory evolution across species allows preservation of gene function, and to what extent these differences influence how cis-regulatory variation arising within species impacts phenotypic change. Here, we investigated these questions using a plant stem cell regulator conserved in expression pattern and function over ~125 million years. Using in-vivo genome editing in two distantly related models, Arabidopsis thaliana (Arabidopsis) and Solanum lycopersicum (tomato), we generated over 70 deletion alleles in the upstream and downstream regions of the stem cell repressor gene CLAVATA3 (CLV3) and compared their individual and combined effects on a shared phenotype, the number of carpels that make fruits. We found that sequences upstream of tomato CLV3 are highly sensitive to even small perturbations compared to its downstream region. In contrast, Arabidopsis CLV3 function is tolerant to severe disruptions both upstream and downstream of the coding sequence. Combining upstream and downstream deletions also revealed a different regulatory outcome. Whereas phenotypic enhancement from adding downstream mutations was predominantly weak and additive in tomato, mutating both regions of Arabidopsis CLV3 caused substantial and synergistic effects, demonstrating distinct distribution and redundancy of functional cis-regulatory sequences. Our results demonstrate remarkable malleability in cis-regulatory structural organization of a deeply conserved plant stem cell regulator and suggest that major reconfiguration of cis-regulatory sequence space is a common yet cryptic evolutionary force altering genotype-to-phenotype relationships from regulatory variation in conserved genes. Finally, our findings underscore the need for lineage-specific dissection of the spatial architecture of cis-regulation to effectively engineer trait variation from conserved productivity genes in crops.

Gapless genome assembly and epigenetic profiles reveal gene regulation of whole-genome triplication in lettuce

BACKGROUND

Lettuce, an important member of the Asteraceae family, is a globally cultivated cash vegetable crop. With a highly complex genome (∼2.5 Gb; 2n = 18) rich in repeat sequences, current lettuce reference genomes exhibit thousands of gaps, impeding a comprehensive understanding of the lettuce genome.

FINDINGS

Here, we present a near-complete gapless reference genome for cutting lettuce with high transformability, using long-read PacBio HiFi and Nanopore sequencing data. In comparison to stem lettuce genome, we identify 127,681 structural variations (SVs, present in 0.41 Gb of sequence), reflecting the divergence of leafy and stem lettuce. Interestingly, these SVs are related to transposons and DNA methylation states. Furthermore, we identify 4,612 whole-genome triplication genes exhibiting high expression levels associated with low DNA methylation levels and high N6-methyladenosine RNA modifications. DNA methylation changes are also associated with activation of genes involved in callus formation.

CONCLUSIONS

Our gapless lettuce genome assembly, an unprecedented achievement in the Asteraceae family, establishes a solid foundation for functional genomics, epigenomics, and crop breeding and sheds new light on understanding the complexity of gene regulation associated with the dynamics of DNA and RNA epigenetics in genome evolution.

Extreme restructuring of cis -regulatory regions controlling a deeply conserved plant stem cell regulator

A striking paradox is that genes with conserved protein sequence, function and expression pattern over deep time often exhibit extremely divergent cis -regulatory sequences. It remains unclear how such drastic cis -regulatory evolution across species allows preservation of gene function, and to what extent these differences influence how cis- regulatory variation arising within species impacts phenotypic change. Here, we investigated these questions using a plant stem cell regulator conserved in expression pattern and function over ∼125 million years. Using in-vivo genome editing in two distantly related models, Arabidopsis thaliana (Arabidopsis) and Solanum lycopersicum (tomato), we generated over 70 deletion alleles in the upstream and downstream regions of the stem cell repressor gene CLAVATA3 ( CLV3 ) and compared their individual and combined effects on a shared phenotype, the number of carpels that make fruits. We found that sequences upstream of tomato CLV3 are highly sensitive to even small perturbations compared to its downstream region. In contrast, Arabidopsis CLV3 function is tolerant to severe disruptions both upstream and downstream of the coding sequence. Combining upstream and downstream deletions also revealed a different regulatory outcome. Whereas phenotypic enhancement from adding downstream mutations was predominantly weak and additive in tomato, mutating both regions of Arabidopsis CLV3 caused substantial and synergistic effects, demonstrating distinct distribution and redundancy of functional cis -regulatory sequences. Our results demonstrate remarkable malleability in cis -regulatory structural organization of a deeply conserved plant stem cell regulator and suggest that major reconfiguration of cis -regulatory sequence space is a common yet cryptic evolutionary force altering genotype-to-phenotype relationships from regulatory variation in conserved genes. Finally, our findings underscore the need for lineage-specific dissection of the spatial architecture of cis -regulation to effectively engineer trait variation from conserved productivity genes in crops.

Author summary

We investigated the evolution of cis -regulatory elements (CREs) and their interactions in the regulation of a plant stem cell regulator gene, CLAVATA3 (CLV3) , in Arabidopsis and tomato. Despite diverging ∼125 million years ago, the function and expression of CLV3 is conserved in these species; however, cis -regulatory sequences upstream and downstream have drastically diverged, preventing identification of conserved non-coding sequences between them. We used CRISPR-Cas9 to engineer dozens of mutations within the cis -regulatory regions of Arabidopsis and tomato CLV3. In tomato, our results show that tomato CLV3 function primarily relies on interactions among CREs in the 5' non-coding region, unlike Arabidopsis CLV3 , which depends on a more balanced distribution of functional CREs between the 5' and 3' regions. Therefore, despite a high degree of functional conservation, our study demonstrates divergent regulatory strategies between two distantly related CLV3 orthologs, with substantial alterations in regulatory sequences, their spatial arrangement, and their relative effects on CLV3 regulation. These results suggest that regulatory regions are not only extremely robust to mutagenesis, but also that the sequences underlying this robustness can be lineage-specific for conserved genes, due to the complex and often redundant interactions among CREs that ensure proper gene function amidst large-scale sequence turnover.

Mining the Utricularia gibba genome for insulator-like elements for genetic engineering

Introduction

Gene expression is often controlled via cis-regulatory elements (CREs) that modulate the production of transcripts. For multi-gene genetic engineering and synthetic biology, precise control of transcription is crucial, both to insulate the transgenes from unwanted native regulation and to prevent readthrough or cross-regulation of transgenes within a multi-gene cassette. To prevent this activity, insulator-like elements, more properly referred to as transcriptional blockers, could be inserted to separate the transgenes so that they are independently regulated. However, only a few validated insulator-like elements are available for plants, and they tend to be larger than ideal.

Methods

To identify additional potential insulator-like sequences, we conducted a genome-wide analysis of Utricularia gibba (humped bladderwort), one of the smallest known plant genomes, with genes that are naturally close together. The 10 best insulator-like candidates were evaluated in vivo for insulator-like activity.

Results

We identified a total of 4,656 intergenic regions with expression profiles suggesting insulator-like activity. Comparisons of these regions across 45 other plant species (representing Monocots, Asterids, and Rosids) show low levels of syntenic conservation of these regions. Genome-wide analysis of unmethylated regions (UMRs) indicates ~87% of the targeted regions are unmethylated; however, interpretation of this is complicated because U. gibba has remarkably low levels of methylation across the genome, so that large UMRs frequently extend over multiple genes and intergenic spaces. We also could not identify any conserved motifs among our selected intergenic regions or shared with existing insulator-like elements for plants. Despite this lack of conservation, however, testing of 10 selected intergenic regions for insulator-like activity found two elements on par with a previously published element (EXOB) while being significantly smaller.

Discussion

Given the small number of insulator-like elements currently available for plants, our results make a significant addition to available tools. The high hit rate (2 out of 10) also implies that more useful sequences are likely present in our selected intergenic regions; additional validation work will be required to identify which will be most useful for plant genetic engineering.

Domains Rearranged Methylase 2 maintains DNA methylation at large DNA hypomethylated shores and long-range chromatin interactions in rice

DNA methylation plays an important role in gene regulation and genomic stability. However, large DNA hypomethylated regions known as DNA methylation valleys (DMVs) or canyons have also been suggested to serve unique regulatory functions, largely unknown in rice (Oryza sativa). Here, we describe the DMVs in rice seedlings, which were highly enriched with developmental and transcription regulatory genes. Further detailed analysis indicated that grand DMVs (gDMVs) might be derived from nuclear integrants of organelle DNA (NORGs). Furthermore, Domains Rearranged Methylase 2 (OsDRM2) maintained DNA methylation at short DMV (sDMV) shores. Epigenetic maps indicated that sDMVs were marked with H3K4me3 and/or H3K27me3, although the loss of DNA methylation had a negligible effect on histone modification within these regions. In addition, we constructed H3K27me3-associated interaction maps for homozygous T-DNA insertion mutant of the gene (osdrm2) and wild type (WT). From a global perspective, most (90%) compartments were stable between osdrm2 and WT plants. At a high resolution, we observed a dramatic loss of long-range chromatin loops in osdrm2, which suffered an extensive loss of non-CG (CHG and CHH, H = A, T, or C) methylation. From another viewpoint, the loss of non-CG methylation at sDMV shores in osdrm2 could disrupt H3K27me3-mediated chromatin interaction networks. Overall, our results demonstrated that DMVs are a key genomic feature in rice and are precisely regulated by epigenetic modifications, including DNA methylation and histone modifications. OsDRM2 maintained DNA methylation at sDMV shores, while OsDRM2 deficiency strongly affected three-dimensional (3D) genome architectures.

Natural methylation epialleles correlate with gene expression in maize

DNA methylation in plants is depleted from cis-regulatory elements in and near genes but is present in some gene bodies, including exons. Methylation in exons solely in the CG context is called gene body methylation (gbM). Methylation in exons in both CG and non-CG contexts is called TE-like methylation (teM). Assigning functions to both forms of methylation in genes has proven to be challenging. Toward that end, we utilized recent genome assemblies, gene annotations, transcription data, and methylome data to quantify common patterns of gene methylation and their relations to gene expression in maize. We found that gbM genes exist in a continuum of CG methylation levels without a clear demarcation between unmethylated genes and gbM genes. Analysis of expression levels across diverse maize stocks and tissues revealed a weak but highly significant positive correlation between gbM and gene expression except in endosperm. gbM epialleles were associated with an approximately 3% increase in steady-state expression level relative to unmethylated epialleles. In contrast to gbM genes, which were conserved and were broadly expressed across tissues, we found that teM genes, which make up about 12% of genes, are mainly silent, are poorly conserved, and exhibit evidence of annotation errors. We used these data to flag teM genes in the 26 NAM founder genome assemblies. While some teM genes are likely functional, these data suggest that the majority are not, and their inclusion can confound the interpretation of whole-genome studies.

Asymmetric variation in DNA methylation during domestication and de-domestication of rice

Hundreds of plant species have been domesticated to feed human civilization, while some crops have undergone de-domestication into agricultural weeds, threatening global food security. To understand the genetic and epigenetic basis of crop domestication and de-domestication, we generated DNA methylomes from 95 accessions of wild rice (Oryza rufipogon L.), cultivated rice (Oryza sativa L.) and weedy rice (O. sativa f. spontanea). We detected a significant decrease in DNA methylation over the course of rice domestication but observed an unexpected increase in DNA methylation through de-domestication. Notably, DNA methylation changes occurred in distinct genomic regions for these 2 opposite stages. Variation in DNA methylation altered the expression of nearby and distal genes through affecting chromatin accessibility, histone modifications, transcription factor binding, and the formation of chromatin loops, which may contribute to morphological changes during domestication and de-domestication of rice. These insights into population epigenomics underlying rice domestication and de-domestication provide resources and tools for epigenetic breeding and sustainable agriculture.

Promoter editing for the genetic improvement of crops

Gene expression plays a fundamental role in the regulation of agronomically important traits in crop plants. The genetic manipulation of plant promoters through genome editing has emerged as an effective strategy to create favorable traits in crops by altering the expression pattern of the pertinent genes. Promoter editing can be applied in a directed manner, where nucleotide sequences associated with favorable traits are precisely generated. Alternatively, promoter editing can also be exploited as a random mutagenic approach to generate novel genetic variations within a designated promoter, from which elite alleles are selected based on their phenotypic effects. Pioneering studies have demonstrated the potential of promoter editing in engineering agronomically important traits as well as in mining novel promoter alleles valuable for plant breeding. In this review, we provide an update on the application of promoter editing in crops for increased yield, enhanced tolerance to biotic and abiotic stresses, and improved quality. We also discuss several remaining technical bottlenecks and how this strategy may be better employed for the genetic improvement of crops in the future.

Exploring the crop epigenome: a comparison of DNA methylation profiling techniques

Epigenetic modifications play a vital role in the preservation of genome integrity and in the regulation of gene expression. DNA methylation, one of the key mechanisms of epigenetic control, impacts growth, development, stress response and adaptability of all organisms, including plants. The detection of DNA methylation marks is crucial for understanding the mechanisms underlying these processes and for developing strategies to improve productivity and stress resistance of crop plants. There are different methods for detecting plant DNA methylation, such as bisulfite sequencing, methylation-sensitive amplified polymorphism, genome-wide DNA methylation analysis, methylated DNA immunoprecipitation sequencing, reduced representation bisulfite sequencing, MS and immuno-based techniques. These profiling approaches vary in many aspects, including DNA input, resolution, genomic region coverage, and bioinformatics analysis. Selecting an appropriate methylation screening approach requires an understanding of all these techniques. This review provides an overview of DNA methylation profiling methods in crop plants, along with comparisons of the efficacy of these techniques between model and crop plants. The strengths and limitations of each methodological approach are outlined, and the importance of considering both technical and biological factors are highlighted. Additionally, methods for modulating DNA methylation in model and crop species are presented. Overall, this review will assist scientists in making informed decisions when selecting an appropriate DNA methylation profiling method.

Transgressive and parental dominant gene expression and cytosine methylation during seed development in Brassica napus hybrids

KEY MESSAGE

Transcriptomic and epigenomic profiling of gene expression and small RNAs during seed and seedling development reveals expression and methylation dominance levels with implications on early stage heterosis in oilseed rape. The enhanced performance of hybrids through heterosis remains a key aspect in plant breeding; however, the underlying mechanisms are still not fully elucidated. To investigate the potential role of transcriptomic and epigenomic patterns in early expression of hybrid vigor, we investigated gene expression, small RNA abundance and genome-wide methylation in hybrids from two distant Brassica napus ecotypes during seed and seedling developmental stages using next-generation sequencing. A total of 31117, 344, 36229 and 7399 differentially expressed genes, microRNAs, small interfering RNAs and differentially methylated regions were identified, respectively. Approximately 70% of the differentially expressed or methylated features displayed parental dominance levels where the hybrid followed the same patterns as the parents. Via gene ontology enrichment and microRNA-target association analyses during seed development, we found copies of reproductive, developmental and meiotic genes with transgressive and paternal dominance patterns. Interestingly, maternal dominance was more prominent in hypermethylated and downregulated features during seed formation, contrasting to the general maternal gamete demethylation reported during gametogenesis in angiosperms. Associations between methylation and gene expression allowed identification of putative epialleles with diverse pivotal biological functions during seed formation. Furthermore, most differentially methylated regions, differentially expressed siRNAs and transposable elements were in regions that flanked genes without differential expression. This suggests that differential expression and methylation of epigenomic features may help maintain expression of pivotal genes in a hybrid context. Differential expression and methylation patterns during seed formation in an F₁ hybrid provide novel insights into genes and mechanisms with potential roles in early heterosis.

Charting plant gene functions in the multi-omics and single-cell era

Despite the increased access to high-quality plant genome sequences, the set of genes with a known function remains far from complete. With the advent of novel bulk and single-cell omics profiling methods, we are entering a new era where advanced and highly integrative functional annotation strategies are being developed to elucidate the functions of all plant genes. Here, we review different multi-omics approaches to improve functional and regulatory gene characterization and highlight the power of machine learning and network biology to fully exploit the complementary information embedded in different omics layers. Finally, we discuss the potential of emerging single-cell methods and algorithms to further increase the resolution, allowing generation of functional insights about plant biology.

Characterization of Gene Regulatory Networks in Plants Using New Methods and Data Types

A major question in plant biology is to understand how plant growth, development, and environmental responses are controlled and coordinated by the activities of regulatory factors. Gene regulatory network (GRN) analyses require integrated approaches that combine experimental approaches with computational analyses. A wide range of experimental approaches and tools are now available, such as targeted perturbation of gene activities, quantitative and cell-type specific measurements of dynamic gene activities, and systematic analysis of the molecular 'hard-wiring' of the systems. At the computational level, different tools and databases are available to study regulatory sequences, including intuitive visualizations to explore data-driven gene regulatory networks in different plant species. Furthermore, advanced data integration approaches have recently been developed to efficiently leverage complementary regulatory data types and learn context-specific networks.

Canalization of genome-wide transcriptional activity in Arabidopsis thaliana accessions by MET1-dependent CG methylation

BACKGROUND

Despite its conserved role on gene expression and transposable element (TE) silencing, genome-wide CG methylation differs substantially between wild Arabidopsis thaliana accessions.

RESULTS

To test our hypothesis that global reduction of CG methylation would reduce epigenomic, transcriptomic, and phenotypic diversity in A. thaliana accessions, we knock out MET1, which is required for CG methylation, in 18 early-flowering accessions. Homozygous met1 mutants in all accessions suffer from common developmental defects such as dwarfism and delayed flowering, in addition to accession-specific abnormalities in rosette leaf architecture, silique morphology, and fertility. Integrated analysis of genome-wide methylation, chromatin accessibility, and transcriptomes confirms that MET1 inactivation greatly reduces CG methylation and alters chromatin accessibility at thousands of loci. While the effects on TE activation are similarly drastic in all accessions, the quantitative effects on non-TE genes vary greatly. The global expression profiles of accessions become considerably more divergent from each other after genome-wide removal of CG methylation, although a few genes with diverse expression profiles across wild-type accessions tend to become more similar in mutants. Most differentially expressed genes do not exhibit altered chromatin accessibility or CG methylation in cis, suggesting that absence of MET1 can have profound indirect effects on gene expression and that these effects vary substantially between accessions.

CONCLUSIONS

Systematic analysis of MET1 requirement in different A. thaliana accessions reveals a dual role for CG methylation: for many genes, CG methylation appears to canalize expression levels, with methylation masking regulatory divergence. However, for a smaller subset of genes, CG methylation increases expression diversity beyond genetically encoded differences.

Dynamic DNA methylation changes reveal tissue-specific gene expression in sugarcane

DNA methylation is an important mechanism for the dynamic regulation of gene expression and silencing of transposons during plant developmental processes. Here, we analyzed genome-wide methylation patterns in sugarcane (Saccharum officinarum) leaves, roots, rinds, and piths at single-base resolution. DNA methylation patterns were similar among the different sugarcane tissues, whereas DNA methylation levels differed. We also found that DNA methylation in different genic regions or sequence contexts plays different roles in gene expression. Differences in methylation among tissues resulted in many differentially methylated regions (DMRs) between tissues, particularly CHH DMRs. Genes overlapping with DMRs tended to be differentially expressed (DEGs) between tissues, and these DMR-associated DEGs were enriched in biological pathways related to tissue function, such as photosynthesis, sucrose synthesis, stress response, transport, and metabolism. Moreover, we observed many DNA methylation valleys (DMVs), which always overlapped with transcription factors (TFs) and sucrose-related genes, such as WRKY, bZIP, WOX, SPS, and FBPase. Collectively, these findings provide significant insights into the complicated interplay between DNA methylation and gene expression and shed light on the epigenetic regulation of sucrose-related genes in sugarcane.

The evolution and function of transposons in epigenetic regulation in response to the environment

Transposable elements (TEs) make up a major proportion of plant genomes. Despite their prevalence genome-wide, TEs are often tossed aside as "junk DNA" since they rarely cause phenotypes, and epigenetic mechanisms silence TEs to prevent them from causing deleterious mutations through movement. While this bleak picture of TEs in genomes is true on average, a growing number of examples across many plant species point to TEs as drivers of phenotypic diversity and novel stress responses. Examples of TE-influenced phenotypes illustrate the many ways that novel transposition events can alter local gene expression and how this relates to potential variation in plant responses to environmental stress. Since TE families and insertions at the locus level lack evolutionary conservation, advancements in the field will require TE experts across diverse species to identify and utilize TE variation in their own systems as a means of crop improvement.

Epigenetic features drastically impact CRISPR-Cas9 efficacy in plants

CRISPR-Cas9-mediated genome editing has been widely adopted for basic and applied biological research in eukaryotic systems. While many studies consider DNA sequences of CRISPR target sites as the primary determinant for CRISPR mutagenesis efficiency and mutation profiles, increasing evidence reveals the substantial role of chromatin context. Nonetheless, most prior studies are limited by the lack of sufficient epigenetic resources and/or by only transiently expressing CRISPR-Cas9 in a short time window. In this study, we leveraged the wealth of high-resolution epigenomic resources in Arabidopsis (Arabidopsis thaliana) to address the impact of chromatin features on CRISPR-Cas9 mutagenesis using stable transgenic plants. Our results indicated that DNA methylation and chromatin features could lead to substantial variations in mutagenesis efficiency by up to 250-fold. Low mutagenesis efficiencies were mostly associated with repressive heterochromatic features. This repressive effect appeared to persist through cell divisions but could be alleviated through substantial reduction of DNA methylation at CRISPR target sites. Moreover, specific chromatin features, such as H3K4me1, H3.3, and H3.1, appear to be associated with significant variation in CRISPR-Cas9 mutation profiles mediated by the non-homologous end joining repair pathway. Our findings provide strong evidence that specific chromatin features could have substantial and lasting impacts on both CRISPR-Cas9 mutagenesis efficiency and DNA double-strand break repair outcomes.

The maize gene maternal derepression of r1 encodes a DNA glycosylase that demethylates DNA and reduces siRNA expression in the endosperm

Demethylation of transposons can activate the expression of nearby genes and cause imprinted gene expression in the endosperm; this demethylation is hypothesized to lead to expression of transposon small interfering RNAs (siRNAs) that reinforce silencing in the next generation through transfer either into egg or embryo. Here we describe maize (Zea mays) maternal derepression of r1 (mdr1), which encodes a DNA glycosylase with homology to Arabidopsis thaliana DEMETER and which is partially responsible for demethylation of thousands of regions in endosperm. Instead of promoting siRNA expression in endosperm, MDR1 activity inhibits it. Methylation of most repetitive DNA elements in endosperm is not significantly affected by MDR1, with an exception of Helitrons. While maternally-expressed imprinted genes preferentially overlap with MDR1 demethylated regions, the majority of genes that overlap demethylated regions are not imprinted. Double mutant megagametophytes lacking both MDR1 and its close homolog DNG102 result in early seed failure, and double mutant microgametophytes fail pre-fertilization. These data establish DNA demethylation by glycosylases as essential in maize endosperm and pollen and suggest that neither transposon repression nor genomic imprinting is its main function in endosperm.

Modeling chromatin state from sequence across angiosperms using recurrent convolutional neural networks

Accessible chromatin regions are critical components of gene regulation but modeling them directly from sequence remains challenging, especially within plants, whose mechanisms of chromatin remodeling are less understood than in animals. We trained an existing deep-learning architecture, DanQ, on data from 12 angiosperm species to predict the chromatin accessibility in leaf of sequence windows within and across species. We also trained DanQ on DNA methylation data from 10 angiosperms because unmethylated regions have been shown to overlap significantly with ACRs in some plants. The across-species models have comparable or even superior performance to a model trained within species, suggesting strong conservation of chromatin mechanisms across angiosperms. Testing a maize (Zea mays L.) held-out model on a multi-tissue chromatin accessibility panel revealed our models are best at predicting constitutively accessible chromatin regions, with diminishing performance as cell-type specificity increases. Using a combination of interpretation methods, we ranked JASPAR motifs by their importance to each model and saw that the TCP and AP2/ERF transcription factor (TF) families consistently ranked highly. We embedded the top three JASPAR motifs for each model at all possible positions on both strands in our sequence window and observed position- and strand-specific patterns in their importance to the model. With our publicly available across-species 'a2z' model it is now feasible to predict the chromatin accessibility and methylation landscape of any angiosperm genome.

Characterization of Transposon-Derived Accessible Chromatin Regions in Rice (Oryza Sativa)

Growing evidence indicates that transposons or transposable elements (TEs)-derived accessible chromatin regions (ACRs) play essential roles in multiple biological processes by interacting with trans-acting factors. However, the function of TE-derived ACRs in the regulation of gene expression in the rice genome has not been well characterized. In this study, we examined the chromatin dynamics in six types of rice tissues and found that ~8% of ACRs were derived from TEs and exhibited distinct levels of accessibility and conservation as compared to those without TEs. TEs exhibited a TE subtype-dependent impact on ACR formation, which can be mediated by changes in the underlying DNA methylation levels. Moreover, we found that tissue-specific TE-derived ACRs might function in the tissue development through the modulation of nearby gene expression. Interestingly, many genes in domestication sweeps were found to overlap with TE-derived ACRs, suggesting their potential functions in the rice domestication. In addition, we found that the expression divergence of 1070 duplicate gene pairs were associated with TE-derived ACRs and had distinct distributions of TEs and ACRs around the transcription start sites (TSSs), which may experience different selection pressures. Thus, our study provides some insights into the biological implications of TE-derived ACRs in the rice genome. Our results imply that these ACRs are likely involved in the regulation of tissue development, rice domestication and functional divergence of duplicated genes.

Elucidating the biology of transcription factor-DNA interaction for accurate identification of cis-regulatory elements

Transcription factors (TFs) play a critical role in determining cell fate decisions by integrating developmental and environmental signals through binding to specific cis-regulatory modules and regulating spatio-temporal specificity of gene expression patterns. Precise identification of functional TF binding sites in time and space not only will revolutionize our understanding of regulatory networks governing cell fate decisions but is also instrumental to uncover how genetic variations cause morphological diversity or disease. In this review, we discuss recent advances in mapping TF binding sites and characterizing the various parameters underlying the complexity of binding site recognition by TFs.

Genome-wide loss of CHH methylation with limited transcriptome changes in Setaria viridis DOMAINS REARRANGED METHYLTRANSFERASE (DRM) mutants

The DOMAINS REARRANGED METHYLTRANSFERASEs (DRMs) are crucial for RNA-directed DNA methylation (RdDM) in plant species. Setaria viridis is a model monocot species with a relatively compact genome that has limited transposable element (TE) content. CRISPR-based genome editing approaches were used to create loss-of-function alleles for the two putative functional DRM genes in S. viridis to probe the role of RdDM. Double mutant (drm1ab) plants exhibit some morphological abnormalities but are fully viable. Whole-genome methylation profiling provided evidence for the widespread loss of methylation in CHH sequence contexts, particularly in regions with high CHH methylation in wild-type plants. Evidence was also found for the locus-specific loss of CG and CHG methylation, even in some regions that lack CHH methylation. Transcriptome profiling identified genes with altered expression in the drm1ab mutants. However, the majority of genes with high levels of CHH methylation directly surrounding the transcription start site or in nearby promoter regions in wild-type plants do not have altered expression in the drm1ab mutant, even when this methylation is lost, suggesting limited regulation of gene expression by RdDM. Detailed analysis of the expression of TEs identified several transposons that are transcriptionally activated in drm1ab mutants. These transposons are likely to require active RdDM for the maintenance of transcriptional repression.

Genome Assembly and Analysis of the Flavonoid and Phenylpropanoid Biosynthetic Pathways in Fingerroot Ginger (Boesenbergia rotunda)

Boesenbergia rotunda (Zingiberaceae), is a high-value culinary and ethno-medicinal plant of Southeast Asia. The rhizomes of this herb have a high flavanone and chalcone content. Here we report the genome analysis of B. rotunda together with a complete genome sequence as a hybrid assembly. B. rotunda has an estimated genome size of 2.4 Gb which is assembled as 27,491 contigs with an N50 size of 12.386 Mb. The highly heterozygous genome encodes 71,072 protein-coding genes and has a 72% repeat content, with class I TEs occupying ~67% of the assembled genome. Fluorescence in situ hybridization of the 18 chromosome pairs at the metaphase showed six sites of 45S rDNA and two sites of 5S rDNA. An SSR analysis identified 238,441 gSSRs and 4604 EST-SSRs with 49 SSR markers common among related species. Genome-wide methylation percentages ranged from 73% CpG, 36% CHG and 34% CHH in the leaf to 53% CpG, 18% CHG and 25% CHH in the embryogenic callus. Panduratin A biosynthetic unigenes were most highly expressed in the watery callus. B rotunda has a relatively large genome with a high heterozygosity and TE content. This assembly and data (PRJNA71294) comprise a source for further research on the functional genomics of B. rotunda, the evolution of the ginger plant family and the potential genetic selection or improvement of gingers.

High-quality genome and methylomes illustrate features underlying evolutionary success of oaks

The genus Quercus, which emerged ∼55 million years ago during globally warm temperatures, diversified into ∼450 extant species. We present a high-quality de novo genome assembly of a California endemic oak, Quercus lobata, revealing features consistent with oak evolutionary success. Effective population size remained large throughout history despite declining since early Miocene. Analysis of 39,373 mapped protein-coding genes outlined copious duplications consistent with genetic and phenotypic diversity, both by retention of genes created during the ancient γ whole genome hexaploid duplication event and by tandem duplication within families, including numerous resistance genes and a very large block of duplicated DUF247 genes, which have been found to be associated with self-incompatibility in grasses. An additional surprising finding is that subcontext-specific patterns of DNA methylation associated with transposable elements reveal broadly-distributed heterochromatin in intergenic regions, similar to grasses. Collectively, these features promote genetic and phenotypic variation that would facilitate adaptability to changing environments.

Epigenetic regulation of seed-specific gene expression by DNA methylation valleys in castor bean

BACKGROUND

Understanding the processes governing angiosperm seed growth and development is essential both for fundamental plant biology and for agronomic purposes. Master regulators of angiosperm seed development are expressed in a seed-specific manner. However, it is unclear how this seed specificity of transcription is established. In some vertebrates, DNA methylation valleys (DMVs) are highly conserved and strongly associated with key developmental genes, but comparable studies in plants are limited to Arabidopsis and soybean. Castor bean (Ricinus communis) is a valuable model system for the study of seed biology in dicots and source of economically important castor oil. Unlike other dicots such as Arabidopsis and soybean, castor bean seeds have a relatively large and persistent endosperm throughout seed development, representing substantial structural differences in mature seeds. Here, we performed an integrated analysis of RNA-seq, whole-genome bisulfite sequencing, and ChIP-seq for various histone marks in the castor bean.

RESULTS

We present a gene expression atlas covering 16 representative tissues and identified 1162 seed-specific genes in castor bean (Ricinus communis), a valuable model for the study of seed biology in dicots. Upon whole-genome DNA methylation analyses, we detected 32,567 DMVs across five tissues, covering ~33% of the castor bean genome. These DMVs are highly hypomethylated during development and conserved across plant species. We found that DMVs have the potential to activate transcription, especially that of tissue-specific genes. Focusing on seed development, we found that many key developmental regulators of seed/endosperm development, including AGL61, AGL62, LEC1, LEC2, ABI3, and WRI1, were located within DMVs. ChIP-seq for five histone modifications in leaves and seeds clearly showed that the vast majority of histone modification peaks were enriched within DMVs, and their remodeling within DMVs has a critical role in the regulation of seed-specific gene expression. Importantly, further experiment analysis revealed that distal DMVs may act as cis-regulatory elements, like enhancers, to activate downstream gene expression.

CONCLUSIONS

Our results point to the importance of DMVs and special distal DMVs behaving like enhancers, in the regulation of seed-specific genes, via the reprogramming of histone modifications within DMVs. Furthermore, these results provide a comprehensive understanding of the epigenetic regulator roles in seed development in castor bean and other important crops.

Decoding the sorghum methylome: understanding epigenetic contributions to agronomic traits

DNA methylation is a chromatin modification that plays an essential role in regulating gene expression and genome stability and it is typically associated with gene silencing and heterochromatin. Owing to its heritability, alterations in the patterns of DNA methylation have the potential to provide for epigenetic inheritance of traits. Contemporary epigenomic technologies provide information beyond sequence variation and could supply alternative sources of trait variation for improvement in crops such as sorghum. Yet, compared with other species such as maize and rice, the sorghum DNA methylome is far less well understood. The distribution of CG, CHG, and CHH methylation in the genome is different compared with other species. CG and CHG methylation levels peak around centromeric segments in the sorghum genome and are far more depleted in the gene dense chromosome arms. The genes regulating DNA methylation in sorghum are also yet to be functionally characterised; better understanding of their identity and functional analysis of DNA methylation machinery mutants in diverse genotypes will be important to better characterise the sorghum methylome. Here, we catalogue homologous genes encoding methylation regulatory enzymes in sorghum based on genes in Arabidopsis, maize, and rice. Discovering variation in the methylome may uncover epialleles that provide extra information to explain trait variation and has the potential to be applied in epigenome-wide association studies or genomic prediction. DNA methylation can also improve genome annotations and discover regulatory elements underlying traits. Thus, improving our knowledge of the sorghum methylome can enhance our understanding of the molecular basis of traits and may be useful to improve sorghum performance.

Cis-regulatory sequences in plants: Their importance, discovery, and future challenges

The identification and characterization of cis-regulatory DNA sequences and how they function to coordinate responses to developmental and environmental cues is of paramount importance to plant biology. Key to these regulatory processes are cis-regulatory modules (CRMs), which include enhancers and silencers. Despite the extraordinary advances in high-quality sequence assemblies and genome annotations, the identification and understanding of CRMs, and how they regulate gene expression, lag significantly behind. This is especially true for their distinguishing characteristics and activity states. Here, we review the current knowledge on CRMs and breakthrough technologies enabling identification, characterization, and validation of CRMs; we compare the genomic distributions of CRMs with respect to their target genes between different plant species, and discuss the role of transposable elements harboring CRMs in the evolution of gene expression. This is an exciting time to study cis-regulomes in plants; however, significant existing challenges need to be overcome to fully understand and appreciate the role of CRMs in plant biology and in crop improvement.

Chromatin and regulatory differentiation between bundle sheath and mesophyll cells in maize

C₄ plants partition photosynthesis enzymes between the bundle sheath (BS) and the mesophyll (M) cells for the better delivery of CO₂ to RuBisCO and to reduce photorespiration. To better understand how C₄ photosynthesis is regulated at the transcriptional level, we performed RNA-seq, ATAC-seq, ChIP-seq and Bisulfite-seq (BS-seq) on BS and M cells isolated from maize leaves. By integrating differentially expressed genes with chromatin features, we found that chromatin accessibility coordinates with epigenetic features, especially H3K27me3 modification and CHH methylation, to regulate cell type-preferentially enriched gene expression. Not only the chromatin-accessible regions (ACRs) proximal to the genes (pACRs) but also the distal ACRs (dACRs) are determinants of cell type-preferentially enriched expression. We further identified cell type-preferentially enriched motifs, e.g. AAAG for BS cells and TGACC/T for M cells, and determined their corresponding transcription factors: DOFs and WRKYs. The complex interaction between cis and trans factors in the preferential expression of C₄ genes was also observed. Interestingly, cell type-preferentially enriched gene expression can be fine-tuned by the coordination of multiple chromatin features. Such coordination may be critical in ensuring the cell type-specific function of key C₄ genes. Based on the observed cell type-preferentially enriched expression pattern and coordinated chromatin features, we predicted a set of functionally unknown genes, e.g. Zm00001d042050 and Zm00001d040659, to be potential key C₄ genes. Our findings provide deep insight into the architectures associated with C₄ gene expression and could serve as a valuable resource to further identify the regulatory mechanisms present in C₄ species.

Epigenome guided crop improvement: current progress and future opportunities

Epigenomics encompasses a broad field of study, including the investigation of chromatin states, chromatin modifications and their impact on gene regulation; as well as the phenomena of epigenetic inheritance. The epigenome is a multi-modal layer of information superimposed on DNA sequences, instructing their usage in gene expression. As such, it is an emerging focus of efforts to improve crop performance. Broadly, this might be divided into avenues that leverage chromatin information to better annotate and decode plant genomes, and into complementary strategies that aim to identify and select for heritable epialleles that control crop traits independent of underlying genotype. In this review, we focus on the first approach, which we term ‘epigenome guided’ improvement. This encompasses the use of chromatin profiles to enhance our understanding of the composition and structure of complex crop genomes. We discuss the current progress and future prospects towards integrating this epigenomic information into crop improvement strategies; in particular for CRISPR/Cas9 gene editing and precision genome engineering. We also highlight some specific opportunities and challenges for grain and horticultural crops.

Prediction of conserved and variable heat and cold stress response in maize using cis-regulatory information

Changes in gene expression are important for responses to abiotic stress. Transcriptome profiling of heat- or cold-stressed maize genotypes identifies many changes in transcript abundance. We used comparisons of expression responses in multiple genotypes to identify alleles with variable responses to heat or cold stress and to distinguish examples of cis- or trans-regulatory variation for stress-responsive expression changes. We used motifs enriched near the transcription start sites (TSSs) for thermal stress-responsive genes to develop predictive models of gene expression responses. Prediction accuracies can be improved by focusing only on motifs within unmethylated regions near the TSS and vary for genes with different dynamic responses to stress. Models trained on expression responses in a single genotype and promoter sequences provided lower performance when applied to other genotypes but this could be improved by using models trained on data from all three genotypes tested. The analysis of genes with cis-regulatory variation provides evidence for structural variants that result in presence/absence of transcription factor binding sites in creating variable responses. This study provides insights into cis-regulatory motifs for heat- and cold-responsive gene expression and defines a framework for developing models to predict expression responses across multiple genotypes.

Arabidopsis bioinformatics: tools and strategies

The sequencing of the Arabidopsis thaliana genome 21 years ago ushered in the genomics era for plant research. Since then, an incredible variety of bioinformatic tools permit easy access to large repositories of genomic, transcriptomic, proteomic, epigenomic and other '-omic' data. In this review, we cover some more recent tools (and highlight the 'classics') for exploring such data in order to help formulate quality, testable hypotheses, often without having to generate new experimental data. We cover tools for examining gene expression and co-expression patterns, undertaking promoter analyses and gene set enrichment analyses, and exploring protein-protein and protein-DNA interactions. We will touch on tools that integrate different data sets at the end of the article.

The arches and spandrels of maize domestication, adaptation, and improvement

People living in the Balsas River basin in southwest México domesticated maize from the bushy grass teosinte. Nine thousand years later, in 2021, Ms. Deb Haaland - a member of the Pueblo of Laguna tribe of New Mexico - wore a dress adorned with a cornstalk when she was sworn in as the Secretary of Interior of the United States of America. This choice of garment highlights the importance of the coevolution of maize and the farmers who, through careful selection over thousands of years, domesticated maize and adapted the physiology and shoot architecture of maize to fit local environments and growth habits. Some traits such as tillering were directly selected on (arches), and others such as tassel size are the by-products (spandrels) of maize evolution. Here, we review current knowledge of the underlying cellular, developmental, physiological, and metabolic processes that were selected by farmers and breeders, which have positioned maize as a top global staple crop.

Intergenic spaces: a new frontier to improving plant health

To more sustainably mitigate the impact of crop diseases on plant health and productivity, there is a need for broader spectrum, long-lasting resistance traits. Defense response (DR) genes, located throughout the genome, participate in cellular and system-wide defense mechanisms to stave off infection by diverse pathogens. This multigenic resistance avoids rapid evolution of a pathogen to overcome host resistance. DR genes reside within resistance-associated quantitative trait loci (QTL), and alleles of DR genes in resistant varieties are more active during pathogen attack relative to susceptible haplotypes. Differential expression of DR genes results from polymorphisms in their regulatory regions, that includes cis-regulatory elements such as transcription factor binding sites as well as features that influence epigenetic structural changes to modulate chromatin accessibility during infection. Many of these elements are found in clusters, known as cis-regulatory modules (CRMs), which are distributed throughout the host genome. Regulatory regions involved in plant-pathogen interactions may also contain pathogen effector binding elements that regulate DR gene expression, and that, when mutated, result in a change in the plants' response. We posit that CRMs and the multiple regulatory elements that comprise them are potential targets for marker-assisted breeding for broad-spectrum, durable disease resistance.

Transcriptional regulation in plants: Using omics data to crack the cis-regulatory code

Innovative omics technologies, advanced bioinformatics, and machine learning methods are rapidly becoming integral tools for plant functional genomics, with tremendous recent advances made in this field. In transcriptional regulation, an initial lag in the accumulation of plant omics data relative to that of animals stimulated the development of computational methods capable of extracting maximum information from the available data sets. Recent comprehensive studies of transcription factor-binding profiles in Arabidopsis and maize and the accumulation of uniformly processed omics data in public databases have brought plant biologists into the big leagues, with many cutting-edge methods available. Here, we summarize the state-of-the-art bioinformatics approaches used to predict or infer the cis-regulatory code behind transcriptional gene regulation, focusing on their plant research applications.

The native cistrome and sequence motif families of the maize ear

Elucidating the transcriptional regulatory networks that underlie growth and development requires robust ways to define the complete set of transcription factor (TF) binding sites. Although TF-binding sites are known to be generally located within accessible chromatin regions (ACRs), pinpointing these DNA regulatory elements globally remains challenging. Current approaches primarily identify binding sites for a single TF (e.g. ChIP-seq), or globally detect ACRs but lack the resolution to consistently define TF-binding sites (e.g. DNAse-seq, ATAC-seq). To address this challenge, we developed MNase-defined cistrome-Occupancy Analysis (MOA-seq), a high-resolution (< 30 bp), high-throughput, and genome-wide strategy to globally identify putative TF-binding sites within ACRs. We used MOA-seq on developing maize ears as a proof of concept, able to define a cistrome of 145,000 MOA footprints (MFs). While a substantial majority (76%) of the known ATAC-seq ACRs intersected with the MFs, only a minority of MFs overlapped with the ATAC peaks, indicating that the majority of MFs were novel and not detected by ATAC-seq. MFs were associated with promoters and significantly enriched for TF-binding and long-range chromatin interaction sites, including for the well-characterized FASCIATED EAR4, KNOTTED1, and TEOSINTE BRANCHED1. Importantly, the MOA-seq strategy improved the spatial resolution of TF-binding prediction and allowed us to identify 215 motif families collectively distributed over more than 100,000 non-overlapping, putatively-occupied binding sites across the genome. Our study presents a simple, efficient, and high-resolution approach to identify putative TF footprints and binding motifs genome-wide, to ultimately define a native cistrome atlas.

Understanding gene regulation remains a central goal of modern biology. Delineating the full set of regulatory DNA elements that orchestrate this regulation requires information at two scales; the broad landscape of accessible chromatin, and the site-specific binding of transcription factors (TFs) at discrete cis-regulatory DNA elements. Here we describe a single assay that uses micrococcal nuclease (MNase) as a structural probe to simultaneously reveal regions of accessible chromatin in addition to high-resolution footprints with signatures of TF-occupied cis-elements. We have used maize developing ear tissue as proof of concept, showing the method detects known TF-binding sites. This genome-wide assay not only defines chromatin landscapes, but crucially enables global discovery and mapping of sequence motifs underlying small footprints of ~30 bp to produce an atlas of candidate TF occupancy.

Stability of DNA methylation and chromatin accessibility in structurally diverse maize genomes

Accessible chromatin and unmethylated DNA are associated with many genes and cis-regulatory elements. Attempts to understand natural variation for accessible chromatin regions (ACRs) and unmethylated regions (UMRs) often rely upon alignments to a single reference genome. This limits the ability to assess regions that are absent in the reference genome assembly and monitor how nearby structural variants influence variation in chromatin state. In this study, de novo genome assemblies for four maize inbreds (B73, Mo17, Oh43, and W22) are utilized to assess chromatin accessibility and DNA methylation patterns in a pan-genome context. A more complete set of UMRs and ACRs can be identified when chromatin data are aligned to the matched genome rather than a single reference genome. While there are UMRs and ACRs present within genomic regions that are not shared between genotypes, these features are 6- to 12-fold enriched within regions between genomes. Characterization of UMRs present within shared genomic regions reveals that most UMRs maintain the unmethylated state in other genotypes with only ∼5% being polymorphic between genotypes. However, the majority (71%) of UMRs that are shared between genotypes only exhibit partial overlaps suggesting that the boundaries between methylated and unmethylated DNA are dynamic. This instability is not solely due to sequence variation as these partially overlapping UMRs are frequently found within genomic regions that lack sequence variation. The ability to compare chromatin properties among individuals with structural variation enables pan-epigenome analyses to study the sources of variation for accessible chromatin and unmethylated DNA.

De novo assembly, annotation, and comparative analysis of 26 diverse maize genomes

We report de novo genome assemblies, transcriptomes, annotations, and methylomes for the 26 inbreds that serve as the founders for the maize nested association mapping population. The number of pan-genes in these diverse genomes exceeds 103,000, with approximately a third found across all genotypes. The results demonstrate that the ancient tetraploid character of maize continues to degrade by fractionation to the present day. Excellent contiguity over repeat arrays and complete annotation of centromeres revealed additional variation in major cytological landmarks. We show that combining structural variation with single-nucleotide polymorphisms can improve the power of quantitative mapping studies. We also document variation at the level of DNA methylation and demonstrate that unmethylated regions are enriched for cis-regulatory elements that contribute to phenotypic variation.

Beyond the gene: epigenetic and cis-regulatory targets offer new breeding potential for the future

For millennia, natural and artificial selection has combined favourable alleles for desirable traits in crop species. While modern plant breeding has achieved steady increases in crop yields over the last century, on the current trajectory we will simply not meet demand by 2045. Novel breeding strategies and sources of genetic variation will be required to sustainably fill predicted yield gaps and meet new consumer preferences. Here, we highlight that stepping up to meet this grand challenge will increasingly require thinking 'beyond the gene'. Significant progress has been made in understanding the contributions of both epigenetic variation and cis-regulatory variation to plant traits. This non-genic variation has great potential in future breeding, synthetic biology and biotechnology applications.

Anno genominis XX: 20 years of Arabidopsis genomics

Twenty years ago, the Arabidopsis thaliana genome sequence was published. This was an important moment as it was the first sequenced plant genome and explicitly brought plant science into the genomics era. At the time, this was not only an outstanding technological achievement, but it was characterized by a superb global collaboration. The Arabidopsis genome was the seed for plant genomic research. Here, we review the development of numerous resources based on the genome that have enabled discoveries across plant species, which has enhanced our understanding of how plants function and interact with their environments.

Genetic and epigenetic variation in transposable element expression responses to abiotic stress in maize

Transposable elements (TEs) pervade most eukaryotic genomes. The repetitive nature of TEs complicates the analysis of their expression. Evaluation of the expression of both TE families (using unique and multi-mapping reads) and specific elements (using uniquely mapping reads) in leaf tissue of three maize (Zea mays) inbred lines subjected to heat or cold stress reveals no evidence for genome-wide activation of TEs; however, some specific TE families generate transcripts only in stress conditions. There is substantial variation for which TE families exhibit stress-responsive expression in the different genotypes. In order to understand the factors that drive expression of TEs, we focused on a subset of families in which we could monitor expression of individual elements. The stress-responsive activation of a TE family can often be attributed to a small number of elements in the family that contains regions lacking DNA methylation. Comparisons of the expression of TEs in different genotypes revealed both genetic and epigenetic variation. Many of the specific TEs that are activated in stress in one inbred are not present in the other inbred, explaining the lack of activation. Among the elements that are shared in both genomes but only expressed in one genotype, we found that many exhibit differences in DNA methylation such that the genotype without expression is fully methylated. This study provides insights into the regulation of expression of TEs in normal and stress conditions and highlights the role of chromatin variation between elements in a family or between genotypes for contributing to expression variation. The highly repetitive nature of many TEs complicates the analysis of their expression. Although most TEs are not expressed, some exhibits expression in certain tissues or conditions. We monitored the expression of both TE families (using unique and multi-mapping reads) and specific elements (using uniquely mapping reads) in leaf tissue of three maize (Zea mays) inbred lines subjected to heat or cold stress. While genome-wide activation of TEs did not occur, some TE families generated transcripts only in stress conditions with variation by genotype. To better understand the factors that drive expression of TEs, we focused on a subset of families in which we could monitor expression of individual elements. In most cases, stress-responsive activation of a TE family was attributed to a small number of elements in the family. The elements that contained small regions lacking DNA methylation regions showed enriched expression while fully methylated elements were rarely expressed in control or stress conditions. The cause of varied expression in the different genotypes was due to both genetic and epigenetic variation. Many specific TEs activated by stress in one inbred were not present in the other inbred. Among the elements shared in both genomes, full methylation inhibited expression in one of the genotypes. This study provides insights into the regulation of TE expression in normal and stress conditions and highlights the role of chromatin variation between elements in a family or between genotypes for contributing to expression.

A cis-regulatory atlas in maize at single-cell resolution

cis-regulatory elements (CREs) encode the genomic blueprints of spatiotemporal gene expression programs enabling highly specialized cell functions. Using single-cell genomics in six maize organs, we determined the cis- and trans-regulatory factors defining diverse cell identities and coordinating chromatin organization by profiling transcription factor (TF) combinatorics, identifying TFs with non-cell-autonomous activity, and uncovering TFs underlying higher-order chromatin interactions. Cell-type-specific CREs were enriched for enhancer activity and within unmethylated long terminal repeat retrotransposons. Moreover, we found cell-type-specific CREs are hotspots for phenotype-associated genetic variants and were targeted by selection during modern maize breeding, highlighting the biological implications of this CRE atlas. Through comparison of maize and Arabidopsis thaliana developmental trajectories, we identified TFs and CREs with conserved and divergent chromatin dynamics, showcasing extensive evolution of gene regulatory networks. In addition to this rich dataset, we developed single-cell analysis software, Socrates, which can be used to understand cis-regulatory variation in any species.

A cis-regulatory atlas in maize at single-cell resolution

cis-regulatory elements (CREs) encode the genomic blueprints of spatiotemporal gene expression programs enabling highly specialized cell functions. Using single-cell genomics in six maize organs, we determined the cis- and trans-regulatory factors defining diverse cell identities and coordinating chromatin organization by profiling transcription factor (TF) combinatorics, identifying TFs with non-cell-autonomous activity, and uncovering TFs underlying higher-order chromatin interactions. Cell-type-specific CREs were enriched for enhancer activity and within unmethylated long terminal repeat retrotransposons. Moreover, we found cell-type-specific CREs are hotspots for phenotype-associated genetic variants and were targeted by selection during modern maize breeding, highlighting the biological implications of this CRE atlas. Through comparison of maize and Arabidopsis thaliana developmental trajectories, we identified TFs and CREs with conserved and divergent chromatin dynamics, showcasing extensive evolution of gene regulatory networks. In addition to this rich dataset, we developed single-cell analysis software, Socrates, which can be used to understand cis-regulatory variation in any species.

Dissecting cis-regulatory control of quantitative trait variation in a plant stem cell circuit

Cis-regulatory mutations underlie important crop domestication and improvement traits^1,2. However, limited allelic diversity has hindered functional dissection of the large number of cis-regulatory elements and their potential interactions, thereby precluding a deeper understanding of how cis-regulatory variation impacts traits quantitatively. Here, we engineered over 60 promoter alleles in two tomato fruit size genes^3,4 to characterize cis-regulatory sequences and study their functional relationships. We found that targeted mutations in conserved promoter sequences of SlCLV3, a repressor of stem cell proliferation^5,6, have a weak impact on fruit locule number. Pairwise combinations of these mutations mildly enhance this phenotype, revealing additive and synergistic relationships between conserved regions and further suggesting even higher-order cis-regulatory interactions within the SlCLV3 promoter. In contrast, SlWUS, a positive regulator of stem cell proliferation repressed by SlCLV3 (refs. ^5,6), is more tolerant to promoter perturbations. Our results show that complex interplay among cis-regulatory variants can shape quantitative variation, and suggest that empirical dissections of this hidden complexity can guide promoter engineering to predictably modify crop traits.