Widespread long-range cis-regulatory elements in the maize genome

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.

PlantDeepSEA, a deep learning-based web service to predict the regulatory effects of genomic variants in plants

Characterizing regulatory effects of genomic variants in plants remains a challenge. Although several tools based on deep-learning models and large-scale chromatin-profiling data have been available to predict regulatory elements and variant effects, no dedicated tools or web services have been reported in plants. Here, we present PlantDeepSEA as a deep learning-based web service to predict regulatory effects of genomic variants in multiple tissues of six plant species (including four crops). PlantDeepSEA provides two main functions. One is called Variant Effector, which aims to predict the effects of sequence variants on chromatin accessibility. Another is Sequence Profiler, a utility that performs 'in silico saturated mutagenesis' analysis to discover high-impact sites (e.g., cis-regulatory elements) within a sequence. When validated on independent test sets, the area under receiver operating characteristic curve of deep learning models in PlantDeepSEA ranges from 0.93 to 0.99. We demonstrate the usability of the web service with two examples. PlantDeepSEA could help to prioritize regulatory causal variants and might improve our understanding of their mechanisms of action in different tissues in plants. PlantDeepSEA is available at http://plantdeepsea.ncpgr.cn/.

Genomic editing of intronic enhancers unveils their role in fine-tuning tissue-specific gene expression in Arabidopsis thaliana

Enhancers located in introns are abundant and play a major role in the regulation of gene expression in mammalian species. By contrast, the functions of intronic enhancers in plants have largely been unexplored and only a handful of plant intronic enhancers have been reported. We performed a genome-wide prediction of intronic enhancers in Arabidopsis thaliana using open chromatin signatures based on DNase I sequencing. We identified 941 candidate intronic enhancers associated with 806 genes in seedling tissue and 1,271 intronic enhancers associated with 1,069 genes in floral tissue. We validated the function of 15 of 21 (71%) of the predicted intronic enhancers in transgenic assays using a reporter gene. We also created deletion lines of three intronic enhancers associated with two different genes using CRISPR/Cas. Deletion of these enhancers, which span key transcription factor binding sites, did not abolish gene expression but caused varying levels of transcriptional repression of their cognate genes. Remarkably, the transcriptional repression of the deletion lines occurred at specific developmental stages and resulted in distinct phenotypic effects on plant morphology and development. Clearly, these three intronic enhancers are important in fine-tuning tissue- and development-specific expression of their cognate genes.

High-resolution genome-wide association study pinpoints metal transporter and chelator genes involved in the genetic control of element levels in maize grain

Despite its importance to plant function and human health, the genetics underpinning element levels in maize grain remain largely unknown. Through a genome-wide association study in the maize Ames panel of nearly 2,000 inbred lines that was imputed with ∼7.7 million SNP markers, we investigated the genetic basis of natural variation for the concentration of 11 elements in grain. Novel associations were detected for the metal transporter genes rte2 (rotten ear2) and irt1 (iron-regulated transporter1) with boron and nickel, respectively. We also further resolved loci that were previously found to be associated with one or more of five elements (copper, iron, manganese, molybdenum, and/or zinc), with two metal chelator and five metal transporter candidate causal genes identified. The nas5 (nicotianamine synthase5) gene involved in the synthesis of nicotianamine, a metal chelator, was found associated with both zinc and iron and suggests a common genetic basis controlling the accumulation of these two metals in the grain. Furthermore, moderate predictive abilities were obtained for the 11 elemental grain phenotypes with two whole-genome prediction models: Bayesian Ridge Regression (0.33–0.51) and BayesB (0.33–0.53). Of the two models, BayesB, with its greater emphasis on large-effect loci, showed ∼4–10% higher predictive abilities for nickel, molybdenum, and copper. Altogether, our findings contribute to an improved genotype-phenotype map for grain element accumulation in maize.

Polycomb-dependent differential chromatin compartmentalization determines gene coregulation in Arabidopsis

In animals, distant H3K27me3-marked Polycomb targets can establish physical interactions forming repressive chromatin hubs. In plants, growing evidence suggests that H3K27me3 acts directly or indirectly to regulate chromatin interactions, although how this histone modification modulates 3D chromatin architecture remains elusive. To decipher the impact of the dynamic deposition of H3K27me3 on the Arabidopsis thaliana nuclear interactome, we combined genetics, transcriptomics, and several 3D epigenomic approaches. By analyzing mutants defective for histone H3K27 methylation or demethylation, we uncovered the crucial role of this chromatin mark in short- and previously unnoticed long-range chromatin loop formation. We found that a reduction in H3K27me3 levels led to a decrease in the interactions within Polycomb-associated repressive domains. Regions with lower H3K27me3 levels in the H3K27 methyltransferase clf mutant established new interactions with regions marked with H3K9ac, a histone modification associated with active transcription, indicating that a reduction in H3K27me3 levels induces a global reconfiguration of chromatin architecture. Altogether, our results reveal that the 3D genome organization is tightly linked to reversible histone modifications that govern chromatin interactions. Consequently, nuclear organization dynamics shapes the transcriptional reprogramming during plant development and places H3K27me3 as a key feature in the coregulation of distant genes.

Conserved noncoding sequences provide insights into regulatory sequence and loss of gene expression in maize

Thousands of species will be sequenced in the next few years; however, understanding how their genomes work, without an unlimited budget, requires both molecular and novel evolutionary approaches. We developed a sensitive sequence alignment pipeline to identify conserved noncoding sequences (CNSs) in the Andropogoneae tribe (multiple crop species descended from a common ancestor ∼18 million years ago). The Andropogoneae share similar physiology while being tremendously genomically diverse, harboring a broad range of ploidy levels, structural variation, and transposons. These contribute to the potential of Andropogoneae as a powerful system for studying CNSs and are factors we leverage to understand the function of maize CNSs. We found that 86% of CNSs were comprised of annotated features, including introns, UTRs, putative cis-regulatory elements, chromatin loop anchors, noncoding RNA (ncRNA) genes, and several transposable element superfamilies. CNSs were enriched in active regions of DNA replication in the early S phase of the mitotic cell cycle and showed different DNA methylation ratios compared to the genome-wide background. More than half of putative cis-regulatory sequences (identified via other methods) overlapped with CNSs detected in this study. Variants in CNSs were associated with gene expression levels, and CNS absence contributed to loss of gene expression. Furthermore, the evolution of CNSs was associated with the functional diversification of duplicated genes in the context of maize subgenomes. Our results provide a quantitative understanding of the molecular processes governing the evolution of CNSs in maize.

Single cell gene regulatory networks in plants: Opportunities for enhancing climate change stress resilience

Global warming poses major challenges for plant survival and agricultural productivity. Thus, efforts to enhance stress resilience in plants are key strategies for protecting food security. Gene regulatory networks (GRNs) are a critical mechanism conferring stress resilience. Until recently, predicting GRNs of the individual cells that make up plants and other multicellular organisms was impeded by aggregate population scale measurements of transcriptome and other genome-scale features. With the advancement of high-throughput single cell RNA-seq and other single cell assays, learning GRNs for individual cells is now possible, in principle. In this article, we report on recent advances in experimental and analytical methodologies for single cell sequencing assays especially as they have been applied to the study of plants. We highlight recent advances and ongoing challenges for scGRN prediction, and finally, we highlight the opportunity to use scGRN discovery for studying and ultimately enhancing abiotic stress resilience in plants.

DNA methylation-linked chromatin accessibility affects genomic architecture in Arabidopsis

Plant DNA methylation, which occurs in three sequence contexts (CG, CHG, and CHH, where H refers to A, T, or C), is established and maintained by different mechanisms. In this study, we present genome-wide chromatin accessibility profiles of Arabidopsis mutants that are deficient in CG, CHG, and/or CHH methylation. Through a combination of DNA methylation, chromatin accessibility, and higher-order chromosome conformation profiling of these mutants, we uncover links between DNA methylation, chromatin accessibility, and 3D genome architecture. These results reveal the interplay between CG and non-CG methylation in heterochromatin maintenance and suggest that DNA methylation can directly impact chromatin structure.

DNA methylation is a major epigenetic modification found across species and has a profound impact on many biological processes. However, its influence on chromatin accessibility and higher-order genome organization remains unclear, particularly in plants. Here, we present genome-wide chromatin accessibility profiles of 18 Arabidopsis mutants that are deficient in CG, CHG, or CHH DNA methylation. We find that DNA methylation in all three sequence contexts impacts chromatin accessibility in heterochromatin. Many chromatin regions maintain inaccessibility when DNA methylation is lost in only one or two sequence contexts, and signatures of accessibility are particularly affected when DNA methylation is reduced in all contexts, suggesting an interplay between different types of DNA methylation. In addition, we found that increased chromatin accessibility was not always accompanied by increased transcription, suggesting that DNA methylation can directly impact chromatin structure by other mechanisms. We also observed that an increase in chromatin accessibility was accompanied by enhanced long-range chromatin interactions. Together, these results provide a valuable resource for chromatin architecture and DNA methylation analyses and uncover a pivotal role for methylation in the maintenance of heterochromatin inaccessibility.

Synthetic promoter designs enabled by a comprehensive analysis of plant core promoters

Targeted engineering of plant gene expression holds great promise for ensuring food security and for producing biopharmaceuticals in plants. However, this engineering requires thorough knowledge of cis-regulatory elements to precisely control either endogenous or introduced genes. To generate this knowledge, we used a massively parallel reporter assay to measure the activity of nearly complete sets of promoters from Arabidopsis, maize and sorghum. We demonstrate that core promoter elements-notably the TATA box-as well as promoter GC content and promoter-proximal transcription factor binding sites influence promoter strength. By performing the experiments in two assay systems, leaves of the dicot tobacco and protoplasts of the monocot maize, we detect species-specific differences in the contributions of GC content and transcription factors to promoter strength. Using these observations, we built computational models to predict promoter strength in both assay systems, allowing us to design highly active promoters comparable in activity to the viral 35S minimal promoter. Our results establish a promising experimental approach to optimize native promoter elements and generate synthetic ones with desirable features.

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.

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.

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

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

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.

Advancing crop genomics from lab to field

Crop genomics remains a key element in ensuring scientific progress to secure global food security. It has been two decades since the sequence of the first plant genome, that of Arabidopsis thaliana, was released, and soon after that the draft sequencing of the rice genome was completed. Since then, the genomes of more than 100 crops have been sequenced, plant genome research has expanded across multiple fronts and the next few years promise to bring further advances spurred by the advent of new technologies and approaches. We are likely to see continued innovations in crop genome sequencing, genetic mapping and the acquisition of multiple levels of biological data. There will be exciting opportunities to integrate genome-scale information across multiple scales of biological organization, leading to advances in our mechanistic understanding of crop biological processes, which will, in turn, provide greater impetus for translation of laboratory results to the field.

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.

How to build a crop plant: Defining the cis-regulatory landscape of maize

The functional regulatory elements of agronomically important plant genomes have been long sought after. Marand et. al. generate a comprehensive atlas of cis-regulatory elements at single cell resolution in maize, providing a powerful resource for inquiries into the rules of multicellular development and for precision crop engineering.

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.

A capture-based assay for detection and characterization of transposon polymorphisms in maize

Transposons can create allelic diversity that affects gene expression and phenotypic diversity. The detection of transposon polymorphisms at a genome-wide scale across a large population is difficult. Here we developed a targeted sequencing approach to monitor transposon polymorphisms of interest. This approach can interrogate the presence or absence of transposons reliably across various genotypes. Using this approach, we genotyped a set of 965 transposon-related presence/absence polymorphisms in a diverse panel of 16 maize (Zea mays L.) inbred lines that are representative of the major maize breeding groups. About 70% of the selected regions can be effectively assayed in each genotype. The consistency between the capture-based assay and PCR-based assay are 98.6% based on analysis of 24 randomly selected transposon polymorphisms. By integrating the transposon polymorphisms data with gene expression data, ∼18% of the assayed transposon polymorphisms were found to be associated with variable gene expression levels. A detailed analysis of 18 polymorphisms in a larger association panel confirmed the effects of 10 polymorphisms, with one of them having stronger association with expression than nearby SNP markers. The effects of seven polymorphisms were tested using a luciferase-based expression assay, and one was confirmed. Together, this study demonstrates that the targeted sequencing assay is an effective way to explore transposon function in a high-throughput manner.

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.

Single-cell transcriptome atlas and chromatin accessibility landscape reveal differentiation trajectories in the rice root

Root development relies on the establishment of meristematic tissues that give rise to distinct cell types that differentiate across defined temporal and spatial gradients. Dissection of the developmental trajectories and the transcriptional networks that underlie them could aid understanding of the function of the root apical meristem in both dicots and monocots. Here, we present a single-cell RNA (scRNA) sequencing and chromatin accessibility survey of rice radicles. By temporal profiling of individual root tip cells we reconstruct continuous developmental trajectories of epidermal cells and ground tissues, and elucidate regulatory networks underlying cell fate determination in these cell lineages. We further identify characteristic processes, transcriptome profiles, and marker genes for these cell types and reveal conserved and divergent root developmental pathways between dicots and monocots. Finally, we demonstrate the potential of the platform for functional genetic studies by using spatiotemporal modeling to identify a rice root meristematic mutant from a cell-specific gene cohort.

Asymmetric epigenome maps of subgenomes reveal imbalanced transcription and distinct evolutionary trends in Brassica napus

The complexity of the epigenome landscape and transcriptional regulation is significantly increased during plant polyploidization, which drives genome evolution and contributes to the increased adaptability to diverse environments. However, a comprehensive epigenome map of Brassica napus is still unavailable. In this study, we performed integrative analysis of five histone modifications, RNA polymerase II occupancy, DNA methylation, and transcriptomes in two B. napus lines (2063A and B409), and established global maps of regulatory elements, chromatin states, and their dynamics for the whole genome (including the An and Cn subgenomes) in four tissue types (young leaf, flower bud, silique, and root) of these two lines. Approximately 65.8% of the genome was annotated with different epigenomic signals. Compared with the Cn subgenome, the An subgenome possesses a higher level of active epigenetic marks and lower level of repressive epigenetic marks. Genes from subgenome-unique regions contribute to the major differences between the An and Cn subgenomes. Asymmetric histone modifications between homeologous gene pairs reflect their biased expression patterns. We identified a novel bivalent chromatin state (with H3K4me1 and H3K27me3) in B. napus that is associated with tissue-specific gene expression. Furthermore, we observed that different types of duplicated genes have discrepant patterns of histone modification and DNA methylation levels. Collectively, our findings provide a valuable epigenetic resource for allopolyploid plants.

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.

Technologies enabling rapid crop improvements for sustainable agriculture: example pennycress (Thlaspi arvense L.)

Growing concerns over food insecurity and ecosystems health related to population growth and climate change have challenged scientists to develop new crops, employing revolutionary technologies in combination with traditional methods. In this review, we discuss the domestication of the oilseed-producing cover crop pennycress, which along with the development of other new crops and improvements to farming practices can provide sustainable solutions to address malnutrition and environmental impacts of production agriculture. We highlight some of the new technologies such as bioinformatics-enabled next-generation sequencing and CRISPR genome editing in combination with traditional mutation breeding that has accelerated pennycress development as a new crop and a potential model system. Furthermore, we provide a brief overview of the technologies that can be integrated for improving pennycress and other crops and the status of pennycress development using these technologies.

Developmental genetics of maize vegetative shoot architecture

More than 1.1 billion tonnes of maize grain were harvested across 197 million hectares in 2019 (FAOSTAT 2020). The vast global productivity of maize is largely driven by denser planting practices, higher yield potential per area of land, and increased yield potential per plant. Shoot architecture, the three-dimensional structural arrangement of the above-ground plant body, is critical to maize grain yield and biomass. Structure of the shoot is integral to all aspects of modern agronomic practices. Here, the developmental genetics of the maize vegetative shoot is reviewed. Plant architecture is ultimately determined by meristem activity, developmental patterning, and growth. The following topics are discussed: shoot apical meristem, leaf architecture, axillary meristem and shoot branching, and intercalary meristem and stem activity. Where possible, classical and current studies in maize developmental genetics, as well as recent advances leveraged by "-omics" analyses, are highlighted within these sections.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11032-021-01208-1.

Improving architectural traits of maize inflorescences

The domestication and improvement of maize resulted in radical changes in shoot architecture relative to its wild progenitor teosinte. In particular, critical modifications involved a reduction of branching and an increase in inflorescence size to meet the needs for human consumption and modern agricultural practices. Maize is a major contributor to global agricultural production by providing large and inexpensive quantities of food, animal feed, and ethanol. Maize is also a classic system for studying the genetic regulation of inflorescence formation and its enlarged female inflorescences directly influence seed production and yield. Studies on the molecular and genetic networks regulating meristem proliferation and maintenance, including receptor-ligand interactions, transcription factor regulation, and hormonal control, provide important insights into maize inflorescence development and reveal potential avenues for the targeted modification of specific architectural traits. In this review, we summarize recent findings on the molecular mechanisms controlling inflorescence formation and discuss how this knowledge can be applied to improve maize productivity in the face of present and future environmental challenges.

ConnecTF: A platform to integrate transcription factor-gene interactions and validate regulatory networks

Deciphering gene regulatory networks (GRNs) is both a promise and challenge of systems biology. The promise lies in identifying key transcription factors (TFs) that enable an organism to react to changes in its environment. The challenge lies in validating GRNs that involve hundreds of TFs with hundreds of thousands of interactions with their genome-wide targets experimentally determined by high-throughput sequencing. To address this challenge, we developed ConnecTF, a species-independent, web-based platform that integrates genome-wide studies of TF-target binding, TF-target regulation, and other TF-centric omic datasets and uses these to build and refine validated or inferred GRNs. We demonstrate the functionality of ConnecTF by showing how integration within and across TF-target datasets uncovers biological insights. Case study 1 uses integration of TF-target gene regulation and binding datasets to uncover TF mode-of-action and identify potential TF partners for 14 TFs in abscisic acid signaling. Case study 2 demonstrates how genome-wide TF-target data and automated functions in ConnecTF are used in precision/recall analysis and pruning of an inferred GRN for nitrogen signaling. Case study 3 uses ConnecTF to chart a network path from NLP7, a master TF in nitrogen signaling, to direct secondary TF2s and to its indirect targets in a Network Walking approach. The public version of ConnecTF (https://ConnecTF.org) contains 3,738,278 TF-target interactions for 423 TFs in Arabidopsis, 839,210 TF-target interactions for 139 TFs in maize (Zea mays), and 293,094 TF-target interactions for 26 TFs in rice (Oryza sativa). The database and tools in ConnecTF will advance the exploration of GRNs in plant systems biology applications for model and crop species.

Chromatin Accessibility Is Associated with Artemisinin Biosynthesis Regulation in Artemisia annua

Glandular trichome (GT) is the dominant site for artemisinin production in Artemisia annua. Several critical genes involved in artemisinin biosynthesis are specifically expressed in GT. However, the molecular mechanism of differential gene expression between GT and other tissue types remains elusive. Chromatin accessibility, defined as the degree to which nuclear molecules are able to interact with chromatin DNA, reflects gene expression capacity to a certain extent. Here, we investigated and compared the landscape of chromatin accessibility in Artemisia annua leaf and GT using the Assay for Transposase-Accessible Chromatin using sequencing (ATAC-seq) technique. We identified 5413 GT high accessible and 4045 GT low accessible regions, and these GT high accessible regions may contribute to GT-specific biological functions. Several GT-specific artemisinin biosynthetic genes, such as DBR2 and CYP71AV1, showed higher accessible regions in GT compared to that in leaf, implying that they might be regulated by chromatin accessibility. In addition, transcription factor binding motifs for MYB, bZIP, C2H2, and AP2 were overrepresented in the highly accessible chromatin regions associated with artemisinin biosynthetic genes in glandular trichomes. Finally, we proposed a working model illustrating the chromatin accessibility dynamics in regulating artemisinin biosynthetic gene expression. This work provided new insights into epigenetic regulation of gene expression in GT.

Evolution of Conserved Noncoding Sequences in Arabidopsis thaliana

Recent pangenome studies have revealed a large fraction of the gene content within a species exhibits presence-absence variation (PAV). However, coding regions alone provide an incomplete assessment of functional genomic sequence variation at the species level. Little to no attention has been paid to noncoding regulatory regions in pangenome studies, though these sequences directly modulate gene expression and phenotype. To uncover regulatory genetic variation, we generated chromosome-scale genome assemblies for thirty Arabidopsis thaliana accessions from multiple distinct habitats and characterized species level variation in Conserved Noncoding Sequences (CNS). Our analyses uncovered not only PAV and positional variation (PosV) but that diversity in CNS is nonrandom, with variants shared across different accessions. Using evolutionary analyses and chromatin accessibility data, we provide further evidence supporting roles for conserved and variable CNS in gene regulation. Additionally, our data suggests that transposable elements contribute to CNS variation. Characterizing species-level diversity in all functional genomic sequences may later uncover previously unknown mechanistic links between genotype and phenotype.

Functional classification of plant long noncoding RNAs: a transcript is known by the company it keeps

The extraordinary maturation in high-throughput sequencing technologies has revealed the existence of a complex network of transcripts in eukaryotic organisms, including thousands of long noncoding (lnc) RNAs with little or no protein-coding capacity. Subsequent discoveries have shown that lncRNAs participate in a wide range of molecular processes, controlling gene expression and protein activity though direct interactions with proteins, DNA or other RNA molecules. Although significant advances have been achieved in the understanding of lncRNA biology in the animal kingdom, the functional characterization of plant lncRNAs is still in its infancy and remains a major challenge. In this review, we report emerging functional and mechanistic paradigms of plant lncRNAs and partner molecules, and discuss how cutting-edge technologies may help to identify and classify yet uncharacterized transcripts into functional groups.

Chromatin accessibility profiling methods

Chromatin accessibility, or the physical access to chromatinized DNA, is a widely studied characteristic of the eukaryotic genome. As active regulatory DNA elements are generally 'accessible', the genome-wide profiling of chromatin accessibility can be used to identify candidate regulatory genomic regions in a tissue or cell type. Multiple biochemical methods have been developed to profile chromatin accessibility, both in bulk and at the single-cell level. Depending on the method, enzymatic cleavage, transposition or DNA methyltransferases are used, followed by high-throughput sequencing, providing a view of genome-wide chromatin accessibility. In this Primer, we discuss these biochemical methods, as well as bioinformatics tools for analysing and interpreting the generated data, and insights into the key regulators underlying developmental, evolutionary and disease processes. We outline standards for data quality, reproducibility and deposition used by the genomics community. Although chromatin accessibility profiling is invaluable to study gene regulation, alone it provides only a partial view of this complex process. Orthogonal assays facilitate the interpretation of accessible regions with respect to enhancer-promoter proximity, functional transcription factor binding and regulatory function. We envision that technological improvements including single-molecule, multi-omics and spatial methods will bring further insight into the secrets of genome regulation.

How Stress Facilitates Phenotypic Innovation Through Epigenetic Diversity

Climate adaptation through phenotypic innovation will become the main challenge for plants during global warming. Plants exhibit a plethora of mechanisms to achieve environmental and developmental plasticity by inducing dynamic alterations of gene regulation and by maximizing natural variation through large population sizes. While successful over long evolutionary time scales, most of these mechanisms lack the short-term adaptive responsiveness that global warming will require. Here, we review our current understanding of the epigenetic regulation of plant genomes, with a focus on stress-response mechanisms and transgenerational inheritance. Field and laboratory-scale experiments on plants exposed to stress have revealed a multitude of temporally controlled, mechanistic strategies integrating both genetic and epigenetic changes on the genome level. We analyze inter- and intra-species population diversity to discuss how methylome differences and transposon activation can be harnessed for short-term adaptive efforts to shape co-evolving traits in response to qualitatively new climate conditions and environmental stress.

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.

Accessible chromatin regions and their functional interrelations with gene transcription and epigenetic modifications in sorghum genome

Accessible chromatin regions (ACRs) provide physical scaffolds to recruit transcriptional co-regulators and displace their nearby nucleosomes in multiple plant species. Characterization of ACRs and investigation of their biological effects in Sorghum bicolor has lagged behind. Regulation of gene expression relies on the transcriptional co-regulators that are recruited to ACRs to affect epigenomic modifications of surrounding nucleosomes. In this study, we employed transposase-accessible chromatin sequencing to identify ACRs and decipher how the presence of ACRs affects gene expression and epigenetic signatures in the Sorghum genome. As a result, 21 077 ACRs, which are mapped to 22.9% of genes and 2.7% of repeats, were identified. The profiling of ACRs on gene structures reveals a narrow and sharp peak around the transcription start site, with relatively weak and broad signals covering the entire gene body and an explicit but wide peak from the transcription termination site to its downstream regions. We discovered that the correlations between gene expression levels and profiled ACR densities are dependent on the positions of ACRs. The occurrence of genic ACRs cumulatively enhances the transcriptional activity of intergenic ACR-associated genes. In addition, an intricate crosstalk among ACRs, gene expression, and epigenetic marks has been unveiled by integrating multiple-omics analyses of whole-genome bisulfite sequencing, 6mA immunoprecipitation followed by sequencing, RNA sequencing, chromatin immunoprecipitation sequencing, and DNase I hypersensitive sites sequencing datasets. Our study provides a genome-wide landscape of ACRs in sorghum, decrypts their interrelations with various epigenetic marks, and sheds new light on their roles in transcriptional regulation.

Chromatin loop anchors contain core structural components of the gene expression machinery in maize

BACKGROUND

Three-dimensional chromatin loop structures connect regulatory elements to their target genes in regions known as anchors. In complex plant genomes, such as maize, it has been proposed that loops span heterochromatic regions marked by higher repeat content, but little is known on their spatial organization and genome-wide occurrence in relation to transcriptional activity.

RESULTS

Here, ultra-deep Hi-C sequencing of maize B73 leaf tissue was combined with gene expression and open chromatin sequencing for chromatin loop discovery and correlation with hierarchical topologically-associating domains (TADs) and transcriptional activity. A majority of all anchors are shared between multiple loops from previous public maize high-resolution interactome datasets, suggesting a highly dynamic environment, with a conserved set of anchors involved in multiple interaction networks. Chromatin loop interiors are marked by higher repeat contents than the anchors flanking them. A small fraction of high-resolution interaction anchors, fully embedded in larger chromatin loops, co-locate with active genes and putative protein-binding sites. Combinatorial analyses indicate that all anchors studied here co-locate with at least 81.5% of expressed genes and 74% of open chromatin regions. Approximately 38% of all Hi-C chromatin loops are fully embedded within hierarchical TAD-like domains, while the remaining ones share anchors with domain boundaries or with distinct domains. Those various loop types exhibit specific patterns of overlap for open chromatin regions and expressed genes, but no apparent pattern of gene expression. In addition, up to 63% of all unique variants derived from a prior public maize eQTL dataset overlap with Hi-C loop anchors. Anchor annotation suggests that < 7% of all loops detected here are potentially devoid of any genes or regulatory elements. The overall organization of chromatin loop anchors in the maize genome suggest a loop modeling system hypothesized to resemble phase separation of repeat-rich regions.

CONCLUSIONS

Sets of conserved chromatin loop anchors mapping to hierarchical domains contains core structural components of the gene expression machinery in maize. The data presented here will be a useful reference to further investigate their function in regard to the formation of transcriptional complexes and the regulation of transcriptional activity in the maize genome.

Single-cell RNA sequencing of developing maize ears facilitates functional analysis and trait candidate gene discovery

Crop productivity depends on activity of meristems that produce optimized plant architectures, including that of the maize ear. A comprehensive understanding of development requires insight into the full diversity of cell types and developmental domains and the gene networks required to specify them. Until now, these were identified primarily by morphology and insights from classical genetics, which are limited by genetic redundancy and pleiotropy. Here, we investigated the transcriptional profiles of 12,525 single cells from developing maize ears. The resulting developmental atlas provides a single-cell RNA sequencing (scRNA-seq) map of an inflorescence. We validated our results by mRNA in situ hybridization and by fluorescence-activated cell sorting (FACS) RNA-seq, and we show how these data may facilitate genetic studies by predicting genetic redundancy, integrating transcriptional networks, and identifying candidate genes associated with crop yield traits.

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.

DNA methylation mutants in Physcomitrella patens elucidate individual roles of CG and non-CG methylation in genome regulation

Cytosine (DNA) methylation in plants regulates the expression of genes and transposons. While methylation in plant genomes occurs at CG, CHG, and CHH sequence contexts, the comparative roles of the individual methylation contexts remain elusive. Here, we present Physcomitrella patens as the second plant system, besides Arabidopsis thaliana, with viable mutants with an essentially complete loss of methylation in the CG and non-CG contexts. In contrast to A. thaliana, P. patens has more robust CHH methylation, similar CG and CHG methylation levels, and minimal cross-talk between CG and non-CG methylation, making it possible to study context-specific effects independently. Our data found CHH methylation to act in redundancy with symmetric methylation in silencing transposons and to regulate the expression of CG/CHG-depleted transposons. Specific elimination of CG methylation did not dysregulate transposons or genes. In contrast, exclusive removal of non-CG methylation massively up-regulated transposons and genes. In addition, comparing two exclusively but equally CG- or CHG-methylated genomes, we show that CHG methylation acts as a greater transcriptional regulator than CG methylation. These results disentangle the transcriptional roles of CG and non-CG, as well as symmetric and asymmetric methylation in a plant genome, and point to the crucial role of non-CG methylation in genome regulation.

Unraveling the 3D Genome Architecture in Plants: Present and Future

The eukaryotic genome has a hierarchical three-dimensional (3D) organization with functional implications for DNA replication, DNA repair, and transcriptional regulation. Over the past decade, scientists have endeavored to elucidate the spatial characteristics and functions of plant genome architecture using high-throughput chromatin conformation capturing technologies such as Hi-C, ChIA-PET, and HiChIP. Here, we systematically review current understanding of chromatin organization in plants at multiple scales. We also discuss the emerging opinions and concepts in 3D genome research, focusing on state-of-the-art 3D genome techniques, RNA-chromatin interactions, liquid-liquid phase separation, and dynamic chromatin alterations. We propose the application of single-cell/single-molecule multi-omics, multiway (DNA-DNA, DNA-RNA, and RNA-RNA interactions) chromatin conformation capturing methods, and proximity ligation-independent 3D genome-mapping technologies to explore chromatin organization structure and function in plants. Such methods could reveal the spatial interactions between trait-related SNPs and their target genes at various spatiotemporal resolutions, and elucidate the molecular mechanisms of the interactions among DNA elements, RNA molecules, and protein factors during the formation of key traits in plants.

Rational design of minimal synthetic promoters for plants

Promoters serve a critical role in establishing baseline transcriptional capacity through the recruitment of proteins, including transcription factors. Previously, a paucity of data for cis-regulatory elements in plants meant that it was challenging to determine which sequence elements in plant promoter sequences contributed to transcriptional function. In this study, we have identified functional elements in the promoters of plant genes and plant pathogens that utilize plant transcriptional machinery for gene expression. We have established a quantitative experimental system to investigate transcriptional function, investigating how identity, density and position contribute to regulatory function. We then identified permissive architectures for minimal synthetic plant promoters enabling the computational design of a suite of synthetic promoters of different strengths. These have been used to regulate the relative expression of output genes in simple genetic devices.

Rational design of minimal synthetic promoters for plants

Promoters serve a critical role in establishing baseline transcriptional capacity through the recruitment of proteins, including transcription factors. Previously, a paucity of data for cis-regulatory elements in plants meant that it was challenging to determine which sequence elements in plant promoter sequences contributed to transcriptional function. In this study, we have identified functional elements in the promoters of plant genes and plant pathogens that utilize plant transcriptional machinery for gene expression. We have established a quantitative experimental system to investigate transcriptional function, investigating how identity, density and position contribute to regulatory function. We then identified permissive architectures for minimal synthetic plant promoters enabling the computational design of a suite of synthetic promoters of different strengths. These have been used to regulate the relative expression of output genes in simple genetic devices.

Interspecific introgression and natural selection in the evolution of Japanese apricot (Prunus mume)

Domestication and population differentiation in crops involve considerable phenotypic changes. The logs of these evolutionary paths, including natural/artificial selection, can be found in the genomes of the current populations. However, these profiles have been little studied in tree crops, which have specific characters, such as long generation time and clonal propagation, maintaining high levels of heterozygosity. We conducted exon-targeted resequencing of 129 genomes in the genus Prunus, mainly Japanese apricot (Prunus mume), and apricot (Prunus armeniaca), plum (Prunus salicina), and peach (Prunus persica). Based on their genome-wide single-nucleotide polymorphisms merged with published resequencing data of 79 Chinese P. mume cultivars, we inferred complete and ongoing population differentiation in P. mume. Sliding window characterization of the indexes for genetic differentiation identified interspecific fragment introgressions between P. mume and related species (plum and apricot). These regions often exhibited strong selective sweeps formed in the paths of establishment or formation of substructures of P. mume, suggesting that P. mume has frequently imported advantageous genes from other species in the subgenus Prunus as adaptive evolution. These findings shed light on the complicated nature of adaptive evolution in a tree crop that has undergone interspecific exchange of genome fragments with natural/artificial selections.

Evolutionary and functional genomics of DNA methylation in maize domestication and improvement

DNA methylation is a ubiquitous chromatin feature, present in 25% of cytosines in the maize genome, but variation and evolution of the methylation landscape during maize domestication remain largely unknown. Here, we leverage whole-genome sequencing (WGS) and whole-genome bisulfite sequencing (WGBS) data on populations of modern maize, landrace, and teosinte (Zea mays ssp. parviglumis) to estimate epimutation rates and selection coefficients. We find weak evidence for direct selection on DNA methylation in any context, but thousands of differentially methylated regions (DMRs) are identified population-wide that are correlated with recent selection. For two trait-associated DMRs, vgt1-DMR and tb1-DMR, HiChIP data indicate that the interactive loops between DMRs and respective downstream genes are present in B73, a modern maize line, but absent in teosinte. Our results enable a better understanding of the evolutionary forces acting on patterns of DNA methylation and suggest a role of methylation variation in adaptive evolution.

Variation and evolution of DNA methylation during maize domestication remain largely unknown. Here, the authors generate genome and methylome sequencing data as well as HiChIP-based interactome data to investigate the adaptive and phenotypic consequences of methylation variations in maize.

Mapping Regulatory Determinants in Plants

The domestication and improvement of many plant species have frequently involved modulation of transcriptional outputs and continue to offer much promise for targeted trait engineering. The cis-regulatory elements (CREs) controlling these trait-associated transcriptional variants however reside within non-coding regions that are currently poorly annotated in most plant species. This is particularly true in large crop genomes where regulatory regions constitute only a small fraction of the total genomic space. Furthermore, relatively little is known about how CREs function to modulate transcription in plants. Therefore understanding where regulatory regions are located within a genome, what genes they control, and how they are structured are important factors that could be used to guide both traditional and synthetic plant breeding efforts. Here, we describe classic examples of regulatory instances as well as recent advances in plant regulatory genomics. We highlight valuable molecular tools that are enabling large-scale identification of CREs and offering unprecedented insight into how genes are regulated in diverse plant species. We focus on chromatin environment, transcription factor (TF) binding, the role of transposable elements, and the association between regulatory regions and target genes.

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.

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

The genomic sequences of crops continue to be produced at a frenetic pace. It remains challenging to develop complete annotations of functional genes and regulatory elements in these genomes. Chromatin accessibility assays enable discovery of functional elements; however, to uncover the full portfolio of cis-elements would require profiling of many combinations of cell types, tissues, developmental stages, and environments. Here, we explore the potential to use DNA methylation profiles to develop more complete annotations. Using leaf tissue in maize, we define ∼100,000 unmethylated regions (UMRs) that account for 5.8% of the genome; 33,375 UMRs are found greater than 2 kb from genes. UMRs are highly stable in multiple vegetative tissues, and they capture the vast majority of accessible chromatin regions from leaf tissue. However, many UMRs are not accessible in leaf, and these represent regions with potential to become accessible in specific cell types or developmental stages. These UMRs often occur near genes that are expressed in other tissues and are enriched for binding sites of transcription factors. The leaf-inaccessible UMRs exhibit unique chromatin modification patterns and are enriched for chromatin interactions with nearby genes. The total UMR space in four additional monocots ranges from 80 to 120 megabases, which is remarkably similar considering the range in genome size of 271 megabases to 4.8 gigabases. In summary, based on the profile from a single tissue, DNA methylation signatures provide powerful filters to distill large genomes down to the small fraction of putative functional genes and regulatory elements.

The matrix revolutions: towards the decoding of the plant chromatin three-dimensional reality

In recent years, we have witnessed a significant increase in studies addressing the three-dimensional (3D) chromatin organization of the plant nucleus. Important advances in chromatin conformation capture (3C)-derived and related techniques have allowed the exploration of the nuclear topology of plants with large and complex genomes, including various crops. In addition, the increase in their resolution has permitted the depiction of chromatin compartmentalization and interactions at the gene scale. These studies have revealed the highly complex mechanisms governing plant nuclear architecture and the remarkable knowledge gaps in this field. Here we discuss the state-of-the-art in plant chromosome architecture, including our knowledge of the hierarchical organization of the genome in 3D space and regarding other nuclear components. Furthermore, we highlight the existence in plants of topologically associated domain (TAD)-like structures that display striking differences from their mammalian counterparts, proposing the concept of ICONS-intergenic condensed spacers. Similarly, we explore recent advances in the study of chromatin loops and R-loops, and their implication in the regulation of gene activity. Finally, we address the impact that polyploidization has had on the chromatin topology of modern crops, and how this is related to phenomena such as subgenome dominance and biased gene retention in these organisms.

Plant and animal chromatin three-dimensional organization: similar structures but different functions

Chromatin is the main carrier of genetic information and is non-randomly distributed within the nucleus. Next-generation sequence-based chromatin conformation capture technologies have enabled us to directly examine its three-dimensional organization at an unprecedented scale and resolution. In the best-studied mammalian models, chromatin folding can be broken down into three hierarchical levels, compartment, domains, and loops, which play important roles in transcriptional regulation. Although similar structures have now been identified in plants, they might not possess exactly the same functions as the mammalian ones. Here, we review recent Hi-C studies in plants, compare plant chromatin structures with their mammalian counterparts, and discuss the differences between plants with different genome sizes.