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

Two antagonistic gene regulatory networks drive Arabidopsis root hair growth at low temperature linked to a low-nutrient environment

Root hair (RH) cells can elongate to several hundred times their initial size, and are an ideal model system for investigating cell size control. Their development is influenced by both endogenous and external signals, which are combined to form an integrative response. Surprisingly, a low-temperature condition of 10°C causes increased RH growth in Arabidopsis and in several monocots, even when the development of the rest of the plant is halted. Previously, we demonstrated a strong correlation between RH growth response and a significant decrease in nutrient availability in the growth medium under low-temperature conditions. However, the molecular basis responsible for receiving and transmitting signals related to the availability of nutrients in the soil, and their relation to plant development, remain largely unknown. We have discovered two antagonic gene regulatory networks (GRNs) controlling RH early transcriptome responses to low temperature. One GNR enhances RH growth and it is commanded by the transcription factors (TFs) ROOT HAIR DEFECTIVE 6 (RHD6), HAIR DEFECTIVE 6-LIKE 2 and 4 (RSL2-RSL4) and a member of the homeodomain leucine zipper (HD-Zip I) group I 16 (AtHB16). On the other hand, a second GRN was identified as a negative regulator of RH growth at low temperature and it is composed by the trihelix TF GT2-LIKE1 (GTL1) and the associated DF1, a previously unidentified MYB-like TF (AT2G01060) and several members of HD-Zip I group (AtHB3, AtHB13, AtHB20, AtHB23). Functional analysis of both GRNs highlights a complex regulation of RH growth response to low temperature, and more importantly, these discoveries enhance our comprehension of how plants synchronize RH growth in response to variations in temperature at the cellular level.

Comprehensive transcriptional regulatory networks in potato through chromatin accessibility and transcriptome under drought and salt stresses

Drought and high salt stress have a great negative impact on potato growth and development. However, the molecular mechanisms by which different tissues and organs of potato plants respond to drought and high salt stress at different stress times lack definition. In this study, we mapped the whole genome of THSs in potato in response to different stresses using RNA-seq and ATAC-seq technologies and constructed the unique and shared transcriptional regulatory networks (TRNs) under different stresses, stress time points, and tissues in potato. The results showed opposite trends for changes in chromatin accessibility and expression of genes under drought and salt stresses. Forty-eight hours and root-specific TRNs were more complex than those of 3 h and leaf, and there were genes with inconsistent drought- and salt-stress-regulated expression only in root-shared TRNs, demonstrating the variability of potato's response to stresses under different tissues and treatment times. These results provide a basis for elucidating the transcriptional mechanisms underlying the specific response of potato to drought or salt stress and the common response to salt and drought stress.

The role of DNA content in shaping chromatin architecture and gene expression

Whole-genome duplication is an evolutionary force that drives speciation in all living kingdoms and is notably prevalent in plants. The evolutionary history of plants involved at least two genomic duplications that significantly expanded the plant morphology and physiology spectrum. Many important crops are polyploids, showing valuable features relative to morphological and stress response traits. After genome duplication, diploidization processes facilitate genomic adjustments to restore disomic inheritance. However, little is known about the chromatin changes triggered by nuclear DNA content alterations. Here, we report that synthetically induced genome duplication leads to chromatinization and significant changes in gene expression, resulting in a transcriptional landscape resembling a natural tetraploid. Interestingly, synthetic diploidization elicits only minor alterations in transcriptional activity and chromatin accessibility compared to the more pronounced effects of tetraploidization. We identified epigenetic factors, including specific histone variants, that showed increased expression following genome duplication and decreased expression after genome reduction. These changes may play a key role in the epigenetic mechanisms underlying the phenotypic complexity after tetraploidization in plants. Our findings shed light on the mechanisms that modulate chromatin accessibility remodeling and gene transcription regulation underlying plant genome adaptation in response to changes in genome size.

A dominant role of transcriptional regulation during the evolution of C4 photosynthesis in Flaveria species

C4 photosynthesis exemplifies convergent evolution of complex traits. Herein, we construct chromosome-scale genome assemblies and perform multi-omics analysis for five Flaveria species, which represent evolutionary stages from C3 to C4 photosynthesis. Chromosome-scale genome sequence analyses reveal a gradual increase in genome size during the evolution of C4 photosynthesis attributed to the expansion of transposable elements. Systematic annotation of genes encoding C4 enzymes and transporters identify additional copies of three C4 enzyme genes through retrotranspositions in C4 species. C4 genes exhibit elevated mRNA and protein abundances, reduced protein-to-RNA ratios, and comparable translation efficiencies in C4 species, highlighting a critical role of transcriptional regulation in C4 evolution. Furthermore, we observe an increased abundance of ethylene response factor (ERF) transcription factors and cognate cis-regulatory elements associated with C4 genes regulation. Altogether, our study provides valuable genomic resources for the Flaveria genus and sheds lights on evolutionary and regulatory mechanisms underlying C4 photosynthesis.

A mobile miR160-triggered transcriptional axis controls root stem cell niche maintenance and regeneration in Arabidopsis

In multicellular organisms, communication between cells is vital for their fate determination. In plants, the quiescent center (QC) signals to adjacent stem cells to maintain them undifferentiated. However, how surrounding stem cells instruct the QC remains poorly understood. Here, we show that in the Arabidopsis root, microRNA160 (miR160) moves from stele stem cells (SSCs) to the QC, where it degrades the mRNAs of two auxin response factors, ARF10 and ARF17. This degradation relieves BRAVO from direct transcriptional repression, maintaining QC quiescence. We further identify that blocking miR160 movement due to DNA damage-induced SSC death and restricted symplastic transport reduces BRAVO and WOX5 expression, leading to QC division to replenish damaged stem cells during root regeneration. Together, our results demonstrate that a transcriptional axis initiated by mobile miR160 regulates the QC and stem cell behavior, advancing our understanding of the communication between stem cells and their surrounding cellular environment.

Interspecific hybridization in Brassica species leads to changes in agronomic traits through the regulation of gene expression by chromatin accessibility and DNA methylation

Interspecific hybridization is a common method in plant breeding to combine traits from different species, resulting in allopolyploidization and significant genetic and epigenetic changes. However, our understanding of genome-wide chromatin and gene expression dynamics during allopolyploidization remains limited. This study generated two Brassica allotriploid hybrids via interspecific hybridization. We observed that accessible chromatin regions (ACRs) and DNA methylation interact to regulates gene expression after interspecific hybridization, ultimately influencing the agronomic traits of the hybrids. In total, 234,649 ACRs were identified in the parental lines and hybrids; the hybridization process induces changes in the distribution and abundance of their accessible chromatin regions, particularly in gene regions and their proximity. Genes associated with proximal ACRs were more highly expressed than those associated with distal and genic ACRs. More than half of novel ACRs drove transgressive gene expression in the hybrids, and the transgressive upregulated genes showed significant enrichment in metal ion binding, especially magnesium ion, calcium ion, and potassium ion binding. We also identified Bna.bZIP11 in the single-parent activation ACR, which binds to BnaA06.UF3GT to promote anthocyanin accumulation in F1 hybrids. DNA methylation plays a role in repressing gene expression, and unmethylated ACRs are more transcriptionally active. Additionally, the A-subgenome ACRs were associated with genome dosage rather than DNA methylation. The interplay among DNA methylation, transposable elements, and sRNA contributes to the dynamic landscape of ACRs during interspecific hybridization, resulting in distinct gene expression patterns on the genome.

Enhanced pollen tube performance at high temperature contributes to thermotolerant fruit and seed production in tomato

Rising temperature extremes during critical reproductive periods threaten the yield of major grain and fruit crops. Flowering plant reproduction depends on the ability of pollen grains to generate a pollen tube, which elongates through the pistil to deliver sperm cells to female gametes for double fertilization. We used tomato as a model fruit crop to determine how high temperature affects the pollen tube growth phase, taking advantage of cultivars noted for fruit production in exceptionally hot growing seasons. We found that exposure to high temperature solely during the pollen tube growth phase limits fruit biomass and seed set more significantly in thermosensitive cultivars than in thermotolerant cultivars. Importantly, we found that pollen tubes from the thermotolerant Tamaulipas cultivar have enhanced growth in vivo and in vitro under high temperature. Analysis of the pollen tube transcriptome's response to high temperature allowed us to define two response modes (enhanced induction of stress responses and higher basal levels of growth pathways repressed by heat stress) associated with reproductive thermotolerance. Importantly, we define key components of the pollen tube stress response, identifying enhanced reactive oxygen species (ROS) homeostasis and pollen tube callose synthesis and deposition as important components of reproductive thermotolerance in Tamaulipas. Our work identifies the pollen tube growth phase as a viable target to enhance reproductive thermotolerance and delineates key pathways that are altered in crop varieties capable of fruiting under high-temperature conditions.

Heat stress causes chromatin accessibility and related gene expression changes in crown tissues of barley (Hordeum vulgare)

Plant responses to stress caused by high temperatures involve changes occurring at the molecular, metabolic, and physiological levels. Understanding the mechanisms by which plants recognize signals to activate this response is a prerequisite for identifying key genes and signaling pathways and for obtaining heat-tolerant plants. We demonstrated the first implementation of an assay for transposase-accessible chromatin to identify open chromatin regions (OCRs) in crown tissues of barley using three genotypes carrying different allelic forms of the sdw1 gene encoding gibberellin 20-oxidase subjected to elevated temperatures. In parallel, we performed gene expression analysis, which allowed us to relate changes in chromatin state to changes in transcriptional activity. The obtained data revealed that the hypersensitive chromatin regions within the genes were more repeatable than those outside the gene intervals. We observed that prolonged exposure to high temperatures increased chromatin accessibility. Genes with OCRs in their regulatory regions were involved in stress signaling and tolerance, including calcium-dependent protein kinase, mitogen-activated protein kinase (MAPK3), receptor-like cytoplasmic kinase (RLK), TIFY domain-containing transcriptional regulator, bZIP transcription factor, and regulatory protein NPR1. The effect of genotype on gene expression was not as pronounced as that of temperature. By combining results from the differential analysis of chromatin accessibility and expression profiles, we identified genes with high temperature-induced changes in chromatin accessibility associated with expression alterations. Importantly, our data revealed a relationship between the loss of chromatin accessibility in response to heat and the downregulation of genes related to gibberellin signaling.

Evolution of plant cell-type-specific cis-regulatory elements

Cis-regulatory elements (CREs) are critical in regulating gene expression, and yet understanding of CRE evolution remains challenging. Here, we constructed a comprehensive single-cell atlas of chromatin accessibility in Oryza sativa, integrating data from 103,911 nuclei representing 126 discrete cell states across nine distinct organs. We used comparative genomics to compare cell-type resolved chromatin accessibility between O. sativa and 57,552 nuclei from four additional grass species (Zea mays, Sorghum bicolor, Panicum miliaceum, and Urochloa fusca). Accessible chromatin regions (ACRs) had different levels of conservation depending on the degree of cell-type specificity. We found a complex relationship between ACRs with conserved noncoding sequences, cell-type specificity, conservation, and tissue-specific switching. Additionally, we found that epidermal ACRs were less conserved compared to other cell types, potentially indicating that more rapid regulatory evolution has occurred in the L1-derived epidermal layer of these species. Finally, we identified and characterized a conserved subset of ACRs that overlapped the repressive histone modification H3K27me3, implicating them as potentially silencer-like CREs maintained by evolution. Collectively, this comparative genomics approach highlights the dynamics of plant cell-type-specific CRE evolution.

Chromatin Accessibility and Gene Expression Vary Between a New and Evolved Autopolyploid of Arabidopsis arenosa

Polyploids arise from whole-genome duplication (WGD) events, which have played important roles in genome evolution across eukaryotes. WGD can increase genome complexity, yield phenotypic novelty, and influence adaptation. Neo-polyploids have been reported to often show seemingly stochastic epigenetic and transcriptional changes, but this leaves open the question whether these changes persist in evolved polyploids. A powerful approach to address this is to compare diploids, neo-polyploids, and evolved polyploids of the same species. Arabidopsis arenosa is a species that allows us to do this-natural diploid and autotetraploid populations exist, while neo-tetraploids can be artificially generated. Here, we use ATAC-seq to assay local chromatin accessibility, and RNA-seq to study gene expression on matched leaf and petal samples from diploid, neo-tetraploid and evolved tetraploid A. arenosa. We found over 8,000 differentially accessible chromatin regions across all samples. These are largely tissue specific and show distinct trends across cytotypes, with roughly 70% arising upon WGD. Interestingly, only a small proportion is associated with expression changes in nearby genes. However, accessibility variation across cytotypes associates strongly with the number of nearby transposable elements. Relatively few genes were differentially expressed upon genome duplication, and ∼60% of these reverted to near-diploid levels in the evolved tetraploid, suggesting that most initial perturbations do not last. Our results provide new insights into how epigenomic and transcriptional mechanisms jointly respond to genome duplication and subsequent evolution of autopolyploids, and importantly, show that one cannot be directly predicted from the other.

Investigating the cis-regulatory basis of C3 and C4 photosynthesis in grasses at single-cell resolution

While considerable knowledge exists about the enzymes pivotal for C4 photosynthesis, much less is known about the cis-regulation important for specifying their expression in distinct cell types. Here, we use single-cell-indexed ATAC-seq to identify cell-type-specific accessible chromatin regions (ACRs) associated with C4 enzymes for five different grass species. This study spans four C4 species, covering three distinct photosynthetic subtypes: Zea mays and Sorghum bicolor (NADP-dependent malic enzyme), Panicum miliaceum (NAD-dependent malic enzyme), Urochloa fusca (phosphoenolpyruvate carboxykinase), along with the C3 outgroup Oryza sativa. We studied the cis-regulatory landscape of enzymes essential across all C4 species and those unique to C4 subtypes, measuring cell-type-specific biases for C4 enzymes using chromatin accessibility data. Integrating these data with phylogenetics revealed diverse co-option of gene family members between species, showcasing the various paths of C4 evolution. Besides promoter proximal ACRs, we found that, on average, C4 genes have two to three distal cell-type-specific ACRs, highlighting the complexity and divergent nature of C4 evolution. Examining the evolutionary history of these cell-type-specific ACRs revealed a spectrum of conserved and novel ACRs, even among closely related species, indicating ongoing evolution of cis-regulation at these C4 loci. This study illuminates the dynamic and complex nature of cis-regulatory elements evolution in C4 photosynthesis, particularly highlighting the intricate cis-regulatory evolution of key loci. Our findings offer a valuable resource for future investigations, potentially aiding in the optimization of C3 crop performance under changing climatic conditions.

Light-responsive transcription factors VvHYH and VvGATA24 mediate wax terpenoid biosynthesis in Vitis vinifera

The cuticular wax that covers the surfaces of plants is the first barrier against environmental stresses and increasingly accumulates with light exposure. However, the molecular basis of light-responsive wax biosynthesis remains elusive. In grape (Vitis vinifera), light exposure resulted in higher wax terpenoid content and lower decay and abscission rates than controls kept in darkness. Assay for transposase-accessible chromatin with high-throughput sequencing and RNA-seq data were integrated to draw the chromatin accessibility and cis-elements regulatory map to identify the potential action sites. Terpenoid synthase 12 (VvTPS12) and 3-hydroxy-3-methylglutaryl-CoA reductase 2 (VvHMGR2) were identified as grape wax biosynthesis targets, while VvHYH and VvGATA24 were identified as terpenoid biosynthesis activators, as more abundant wax crystals and higher wax terpenoid content were observed in transiently overexpressed grape berries and Nicotiana benthamiana leaves. The interaction between VvHYH and the open chromatin of VvTPS12 was confirmed qualitatively using a dual luciferase assay and quantitatively using surface plasma resonance, with an equilibrium dissociation constant of 2.81 nm identified via the latter approach. Molecular docking simulation implied the structural nature of this interaction, indicating that 24 amino acid residues of VvHYH, including Arg106A, could bind to the VvTPS12 G-box cis-element. VvGATA24 directly bound to the open chromatin of VvHMGR2, with an equilibrium dissociation constant of 8.59 nm. Twelve amino acid residues of VvGATA24, including Pro218B, interacted with the VvHMGR2 GATA-box cis-element. Our work characterizes the mechanism underlying light-mediated wax terpenoid biosynthesis and provides gene targets for future molecular breeding.

Genome-wide epigenetic dynamics of tea leaves under mechanical wounding stress during oolong tea postharvest processing

Understanding the epigenetic responses to mechanical wounding stress during the postharvest processing of oolong tea provides insight into the reprogramming of the tea genome and its impact on tea quality. Here, we characterized the 5mC DNA methylation and chromatin accessibility landscapes of tea leaves subjected to mechanical wounding stress during the postharvest processing of oolong tea. Analysis of the differentially methylated regions and preferentially accessible promoters revealed many overrepresented TF-binding motifs, highlighting sets of TFs that are likely important for the quality of oolong tea. Within these sets, we constructed a chromatin accessibility-mediated gene regulatory network specific to mechanical wounding stress. In combination with the results of the TF-centred yeast one-hybrid assay, we identified potential binding sites of CsMYC2 and constructed a gene regulatory network centred on CsMYC2, clarifying the potential regulatory role of CsMYC2 in the postharvest processing of oolong tea. Interestingly, highly accessible chromatin and hypomethylated cytosine were found to coexist in the promoter region of the indole biosynthesis gene (tryptophan synthase β-subunit, CsTSB) under wounding stress, which indicates that these two important epigenetic regulatory mechanisms are jointly involved in regulating the synthesis of indole during the postharvest processing of oolong tea. These findings improve our understanding of the epigenetic regulatory mechanisms involved in quality formation during the postharvest processing of oolong tea.

Single-Cell Transcriptomics Applied in Plants

Single-cell RNA sequencing (scRNA-seq) is a high-tech method for characterizing the expression patterns of heterogeneous cells in the same tissue and has changed our evaluation of biological systems by increasing the number of individual cells analyzed. However, the full potential of scRNA-seq, particularly in plant science, has not yet been elucidated. To explore the utilization of scRNA-seq technology in plants, we firstly conducted a comprehensive review of significant scRNA-seq findings in the past few years. Secondly, we introduced the research and applications of scRNA-seq technology to plant tissues in recent years, primarily focusing on model plants, crops, and wood. We then offered five databases that could facilitate the identification of distinct expression marker genes for various cell types. Finally, we analyzed the potential problems, challenges, and directions for applying scRNA-seq in plants, with the aim of providing a theoretical foundation for the better use of this technique in future plant research.

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.

Functional Characterization of Accessible Chromatin in Common Wheat

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

Enhanced pollen tube performance at high temperature contributes to thermotolerant fruit production in tomato

Rising temperature extremes during critical reproductive periods threaten the yield of major grain and fruit crops. Flowering plant reproduction depends on development of sufficient numbers of pollen grains and on their ability to generate a cellular extension, the pollen tube, which elongates through the pistil to deliver sperm cells to female gametes for double fertilization. These critical phases of the life cycle are sensitive to temperature and limit productivity under high temperature (HT). Previous studies have investigated the effects of HT on pollen development, but little is known about how HT applied during the pollen tube growth phase affects fertility. Here, we used tomato as a model fruit crop to determine how HT affects the pollen tube growth phase, taking advantage of cultivars noted for fruit production in exceptionally hot growing seasons. We found that exposure to HT solely during the pollen tube growth phase limits fruit biomass and seed set more significantly in thermosensitive cultivars than in thermotolerant cultivars. Importantly, we found that pollen tubes from the thermotolerant Tamaulipas cultivar have enhanced growth in vivo and in vitro under HT. Analysis of the pollen tube transcriptome's response to HT allowed us to develop hypotheses for the molecular basis of cellular thermotolerance in the pollen tube and we define two response modes (enhanced induction of stress responses, and higher basal levels of growth pathways repressed by heat stress) associated with reproductive thermotolerance. Importantly, we define key components of the pollen tube stress response identifying enhanced ROS homeostasis and pollen tube callose synthesis and deposition as important components of reproductive thermotolerance in Tamaulipas. Our work identifies the pollen tube growth phase as a viable target to enhance reproductive thermotolerance and delineates key pathways that are altered in crop varieties capable of fruiting under HT conditions.

A foundational large language model for edible plant genomes

Significant progress has been made in the field of plant genomics, as demonstrated by the increased use of high-throughput methodologies that enable the characterization of multiple genome-wide molecular phenotypes. These findings have provided valuable insights into plant traits and their underlying genetic mechanisms, particularly in model plant species. Nonetheless, effectively leveraging them to make accurate predictions represents a critical step in crop genomic improvement. We present AgroNT, a foundational large language model trained on genomes from 48 plant species with a predominant focus on crop species. We show that AgroNT can obtain state-of-the-art predictions for regulatory annotations, promoter/terminator strength, tissue-specific gene expression, and prioritize functional variants. We conduct a large-scale in silico saturation mutagenesis analysis on cassava to evaluate the regulatory impact of over 10 million mutations and provide their predicted effects as a resource for variant characterization. Finally, we propose the use of the diverse datasets compiled here as the Plants Genomic Benchmark (PGB), providing a comprehensive benchmark for deep learning-based methods in plant genomic research. The pre-trained AgroNT model is publicly available on HuggingFace at https://huggingface.co/InstaDeepAI/agro-nucleotide-transformer-1b  for future research purposes.

Evolutionary insights into the organization of chromatin structure and landscape of transcriptional regulation in plants

Development of complex traits necessitates the functioning and coordination of intricate regulatory networks involving multiple genes. Understanding 3D chromatin structure can facilitate insight into the regulation of gene expression by regulatory elements. This potential, of visualizing the role of chromatin organization in the evolution and function of regulatory elements, remains largely unexplored. Here, we describe new perspectives that arise from the dual considerations of sequence variation of regulatory elements and chromatin structure, with a special focus on whole-genome doubling or polyploidy. We underscore the significance of hierarchical chromatin organization in gene regulation during evolution. In addition, we describe strategies for exploring chromatin organization in future investigations of regulatory evolution in plants, enabling insights into the evolutionary influence of regulatory elements on gene expression and, hence, phenotypes.

Identification and functional characterization of conserved cis-regulatory elements responsible for early fruit development in cucurbit crops

A number of cis-regulatory elements (CREs) conserved during evolution have been found to be responsible for phenotypic novelty and variation. Cucurbit crops such as cucumber (Cucumis sativus), watermelon (Citrullus lanatus), melon (Cucumis melo), and squash (Cucurbita maxima) develop fruits from an inferior ovary and share some similar biological processes during fruit development. Whether conserved regulatory sequences play critical roles in fruit development of cucurbit crops remains to be explored. In six well-studied cucurbit species, we identified 392,438 conserved noncoding sequences (CNSs), including 82,756 that are specific to cucurbits, by comparative genomics. Genome-wide profiling of accessible chromatin regions (ACRs) and gene expression patterns mapped 20,865 to 43,204 ACRs and their potential target genes for two fruit tissues at two key developmental stages in six cucurbits. Integrated analysis of CNSs and ACRs revealed 4,431 syntenic orthologous CNSs, including 1,687 cucurbit-specific CNSs that overlap with ACRs that are present in all six cucurbit crops and that may regulate the expression of 757 adjacent orthologous genes. CRISPR mutations targeting two CNSs present in the 1,687 cucurbit-specific sequences resulted in substantially altered fruit shape and gene expression patterns of adjacent NAC1 (NAM, ATAF1/2, and CUC2) and EXT-like (EXTENSIN-like) genes, validating the regulatory roles of these CNSs in fruit development. These results not only provide a number of target CREs for cucurbit crop improvement, but also provide insight into the roles of CREs in plant biology and during evolution.

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.

Predmoter-cross-species prediction of plant promoter and enhancer regions

Motivation

Identifying cis-regulatory elements (CREs) is crucial for analyzing gene regulatory networks. Next generation sequencing methods were developed to identify CREs but represent a considerable expenditure for targeted analysis of few genomic loci. Thus, predicting the outputs of these methods would significantly cut costs and time investment.

Results

We present Predmoter, a deep neural network that predicts base-wise Assay for Transposase Accessible Chromatin using sequencing (ATAC-seq) and histone Chromatin immunoprecipitation DNA-sequencing (ChIP-seq) read coverage for plant genomes. Predmoter uses only the DNA sequence as input. We trained our final model on 21 species for 13 of which ATAC-seq data and for 17 of which ChIP-seq data was publicly available. We evaluated our models on Arabidopsis thaliana and Oryza sativa. Our best models showed accurate predictions in peak position and pattern for ATAC- and histone ChIP-seq. Annotating putatively accessible chromatin regions provides valuable input for the identification of CREs. In conjunction with other in silico data, this can significantly reduce the search space for experimentally verifiable DNA-protein interaction pairs.

Availability and implementation

The source code for Predmoter is available at: https://github.com/weberlab-hhu/Predmoter. Predmoter takes a fasta file as input and outputs h5, and optionally bigWig and bedGraph files.

Investigating the cis-Regulatory Basis of C3 and C4 Photosynthesis in Grasses at Single-Cell Resolution

While considerable knowledge exists about the enzymes pivotal for C4 photosynthesis, much less is known about the cis-regulation important for specifying their expression in distinct cell types. Here, we use single-cell-indexed ATAC-seq to identify cell-type-specific accessible chromatin regions (ACRs) associated with C4 enzymes for five different grass species. This study spans four C4 species, covering three distinct photosynthetic subtypes: Zea mays and Sorghum bicolor (NADP-ME), Panicum miliaceum (NAD-ME), Urochloa fusca (PEPCK), along with the C3 outgroup Oryza sativa. We studied the cis-regulatory landscape of enzymes essential across all C4 species and those unique to C4 subtypes, measuring cell-type-specific biases for C4 enzymes using chromatin accessibility data. Integrating these data with phylogenetics revealed diverse co-option of gene family members between species, showcasing the various paths of C4 evolution. Besides promoter proximal ACRs, we found that, on average, C4 genes have two to three distal cell-type-specific ACRs, highlighting the complexity and divergent nature of C4 evolution. Examining the evolutionary history of these cell-type-specific ACRs revealed a spectrum of conserved and novel ACRs, even among closely related species, indicating ongoing evolution of cis-regulation at these C4 loci. This study illuminates the dynamic and complex nature of CRE evolution in C4 photosynthesis, particularly highlighting the intricate cis-regulatory evolution of key loci. Our findings offer a valuable resource for future investigations, potentially aiding in the optimization of C3 crop performance under changing climatic conditions.

Deep learning the cis-regulatory code for gene expression in selected model plants

Elucidating the relationship between non-coding regulatory element sequences and gene expression is crucial for understanding gene regulation and genetic variation. We explored this link with the training of interpretable deep learning models predicting gene expression profiles from gene flanking regions of the plant species Arabidopsis thaliana, Solanum lycopersicum, Sorghum bicolor, and Zea mays. With over 80% accuracy, our models enabled predictive feature selection, highlighting e.g. the significant role of UTR regions in determining gene expression levels. The models demonstrated remarkable cross-species performance, effectively identifying both conserved and species-specific regulatory sequence features and their predictive power for gene expression. We illustrated the application of our approach by revealing causal links between genetic variation and gene expression changes across fourteen tomato genomes. Lastly, our models efficiently predicted genotype-specific expression of key functional gene groups, exemplified by underscoring known phenotypic and metabolic differences between Solanum lycopersicum and its wild, drought-resistant relative, Solanum pennellii.

Protocol to isolate nuclei from Chlamydomonas reinhardtii for ATAC sequencing

The isolation of sufficient amounts of intact nuclei is essential to obtain high-resolution maps of chromatin accessibility via assay for transposase-accessible chromatin using sequencing (ATAC-seq). Here, we present a protocol for tag-free isolation of nuclei from both cell walled and cell wall-deficient strains of the green model alga Chlamydomonas reinhardtii at a suitable quality for ATAC-seq. We describe steps for nuclei isolation, quantification, and downstream ATAC-seq. This protocol is optimized to shorten the time of isolation and quantification of nuclei.

Differences in transcription initiation directionality underlie distinctions between plants and animals in chromatin modification patterns at genes and cis-regulatory elements

Transcriptional initiation is among the first regulated steps controlling eukaryotic gene expression. High-throughput profiling of fungal and animal genomes has revealed that RNA Polymerase II often initiates transcription in both directions at the promoter transcription start site, but generally only elongates productively into the gene body. Additionally, Pol II can initiate transcription in both directions at cis-regulatory elements such as enhancers. These bidirectional RNA Polymerase II initiation events can be observed directly with methods that capture nascent transcripts, and they are also revealed indirectly by the presence of transcription-associated histone modifications on both sides of the transcription start site or cis-regulatory elements. Previous studies have shown that nascent RNAs and transcription-associated histone modifications in the model plant Arabidopsis thaliana accumulate mainly in the gene body, suggesting that transcription does not initiate widely in the upstream direction from genes in this plant. We compared transcription-associated histone modifications and nascent transcripts at both transcription start sites and cis-regulatory elements in A. thaliana, Drosophila melanogaster, and Homo sapiens. Our results provide evidence for mostly unidirectional RNA Polymerase II initiation at both promoters and gene-proximal cis-regulatory elements of A. thaliana, whereas bidirectional transcription initiation is observed widely at promoters in both D. melanogaster and H. sapiens, as well as cis-regulatory elements in Drosophila. Furthermore, the distribution of transcription-associated histone modifications around transcription start sites in the Oryza sativa (rice) and Glycine max (soybean) genomes suggests that unidirectional transcription initiation is the norm in these genomes as well. These results suggest that there are fundamental differences in transcriptional initiation directionality between flowering plant and metazoan genomes, which are manifested as distinct patterns of chromatin modifications around RNA polymerase initiation sites.

PlantFUNCO: Integrative Functional Genomics Database Reveals Clues into Duplicates Divergence Evolution

Evolutionary epigenomics and, more generally, evolutionary functional genomics, are emerging fields that study how non-DNA-encoded alterations in gene expression regulation are an important form of plasticity and adaptation. Previous evidence analyzing plants' comparative functional genomics has mostly focused on comparing same assay-matched experiments, missing the power of heterogeneous datasets for conservation inference. To fill this gap, we developed PlantFUN(ctional)CO(nservation) database, which is constituted by several tools and two main resources: interspecies chromatin states and functional genomics conservation scores, presented and analyzed in this work for three well-established plant models (Arabidopsis thaliana, Oryza sativa, and Zea mays). Overall, PlantFUNCO elucidated evolutionary information in terms of cross-species functional agreement. Therefore, providing a new complementary comparative-genomics source for assessing evolutionary studies. To illustrate the potential applications of this database, we replicated two previously published models predicting genetic redundancy in A. thaliana and found that chromatin states are a determinant of paralogs degree of functional divergence. These predictions were validated based on the phenotypes of mitochondrial alternative oxidase knockout mutants under two different stressors. Taking all the above into account, PlantFUNCO aim to leverage data diversity and extrapolate molecular mechanisms findings from different model organisms to determine the extent of functional conservation, thus, deepening our understanding of how plants epigenome and functional noncoding genome have evolved. PlantFUNCO is available at https://rocesv.github.io/PlantFUNCO.

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.

BCL7A and BCL7B potentiate SWI/SNF-complex-mediated chromatin accessibility to regulate gene expression and vegetative phase transition in plants

Switch defective/sucrose non-fermentable (SWI/SNF) chromatin remodeling complexes are multi-subunit machineries that establish and maintain chromatin accessibility and gene expression by regulating chromatin structure. However, how the remodeling activities of SWI/SNF complexes are regulated in eukaryotes remains elusive. B-cell lymphoma/leukemia protein 7 A/B/C (BCL7A/B/C) have been reported as subunits of SWI/SNF complexes for decades in animals and recently in plants; however, the role of BCL7 subunits in SWI/SNF function remains undefined. Here, we identify a unique role for plant BCL7A and BCL7B homologous subunits in potentiating the genome-wide chromatin remodeling activities of SWI/SNF complexes in plants. BCL7A/B require the catalytic ATPase BRAHMA (BRM) to assemble with the signature subunits of the BRM-Associated SWI/SNF complexes (BAS) and for genomic binding at a subset of target genes. Loss of BCL7A and BCL7B diminishes BAS-mediated genome-wide chromatin accessibility without changing the stability and genomic targeting of the BAS complex, highlighting the specialized role of BCL7A/B in regulating remodeling activity. We further show that BCL7A/B fine-tune the remodeling activity of BAS complexes to generate accessible chromatin at the juvenility resetting region (JRR) of the microRNAs MIR156A/C for plant juvenile identity maintenance. In summary, our work uncovers the function of previously elusive SWI/SNF subunits in multicellular eukaryotes and provides insights into the mechanisms whereby plants memorize the juvenile identity through SWI/SNF-mediated control of chromatin accessibility.

Application of Single-Cell Assay for Transposase-Accessible Chromatin with High Throughput Sequencing in Plant Science: Advances, Technical Challenges, and Prospects

The Single-cell Assay for Transposase-Accessible Chromatin with high throughput sequencing (scATAC-seq) has gained increasing popularity in recent years, allowing for chromatin accessibility to be deciphered and gene regulatory networks (GRNs) to be inferred at single-cell resolution. This cutting-edge technology now enables the genome-wide profiling of chromatin accessibility at the cellular level and the capturing of cell-type-specific cis-regulatory elements (CREs) that are masked by cellular heterogeneity in bulk assays. Additionally, it can also facilitate the identification of rare and new cell types based on differences in chromatin accessibility and the charting of cellular developmental trajectories within lineage-related cell clusters. Due to technical challenges and limitations, the data generated from scATAC-seq exhibit unique features, often characterized by high sparsity and noise, even within the same cell type. To address these challenges, various bioinformatic tools have been developed. Furthermore, the application of scATAC-seq in plant science is still in its infancy, with most research focusing on root tissues and model plant species. In this review, we provide an overview of recent progress in scATAC-seq and its application across various fields. We first conduct scATAC-seq in plant science. Next, we highlight the current challenges of scATAC-seq in plant science and major strategies for cell type annotation. Finally, we outline several future directions to exploit scATAC-seq technologies to address critical challenges in plant science, ranging from plant ENCODE(The Encyclopedia of DNA Elements) project construction to GRN inference, to deepen our understanding of the roles of CREs in plant biology.

Genome-wide DNase I-hypersensitive site assay reveals distinct genomic distributions and functional features of open chromatin in autopolyploid sugarcane

The characterization of cis-regulatory DNA elements (CREs) is essential for deciphering the regulation of gene expression in eukaryotes. Although there have been endeavors to identify CREs in plants, the properties of CREs in polyploid genomes are still largely unknown. Here, we conducted the genome-wide identification of DNase I-hypersensitive sites (DHSs) in leaf and stem tissues of the auto-octoploid species Saccharum officinarum. We revealed that DHSs showed highly similar distributions in the genomes of these two S. officinarum tissues. Notably, we observed that approximately 74% of DHSs were located in distal intergenic regions, suggesting considerable differences in the abundance of distal CREs between S. officinarum and other plants. Leaf- and stem-dependent transcriptional regulatory networks were also developed by mining the binding motifs of transcription factors (TFs) from tissue-specific DHSs. Four TEOSINTE BRANCHED 1, CYCLOIDEA, and PCF1 (TCP) TFs (TCP2, TCP4, TCP7, and TCP14) and two ethylene-responsive factors (ERFs) (ERF109 and ERF03) showed strong causal connections with short binding distances from each other, pointing to their possible roles in the regulatory networks of leaf and stem development. Through functional validation in transiently transgenic protoplasts, we isolate a set of tissue-specific promoters. Overall, the DHS maps presented here offer a global view of the potential transcriptional regulatory elements in polyploid sugarcane and can be expected to serve as a valuable resource for both transcriptional network elucidation and genome editing in sugarcane breeding.

Transcription factor NAC1 activates expression of peptidase-encoding AtCEPs in roots to limit root hair growth

Plant genomes encode a unique group of papain-type Cysteine EndoPeptidases (CysEPs) containing a KDEL endoplasmic reticulum (ER) retention signal (KDEL-CysEPs or CEPs). CEPs process the cell-wall scaffolding EXTENSIN (EXT) proteins that regulate de novo cell-wall formation and cell expansion. Since CEPs cleave EXTs and EXT-related proteins, acting as cell-wall-weakening agents, they may play a role in cell elongation. The Arabidopsis (Arabidopsis thaliana) genome encodes 3 CEPs (AtCPE1-AtCEP3). Here, we report that the genes encoding these 3 Arabidopsis CEPs are highly expressed in root-hair (RH) cell files. Single mutants have no evident abnormal RH phenotype, but atcep1-3 atcep3-2 and atcep1-3 atcep2-2 double mutants have longer RHs than wild-type (Wt) plants, suggesting that expression of AtCEPs in root trichoblasts restrains polar elongation of the RH. We provide evidence that the transcription factor NAC1 (petunia NAM and Arabidopsis ATAF1, ATAF2, and CUC2) activates AtCEPs expression in roots to limit RH growth. Chromatin immunoprecipitation indicates that NAC1 binds to the promoter of AtCEP1, AtCEP2, and, to a lower extent, AtCEP3 and may directly regulate their expression. Inducible NAC1 overexpression increases AtCEP1 and AtCEP2 transcript levels in roots and leads to reduced RH growth while the loss of function nac1-2 mutation reduces AtCEP1-AtCEP3 gene expression and enhances RH growth. Likewise, expression of a dominant chimeric NAC1-SRDX repressor construct leads to increased RH length. Finally, we show that RH cell walls in the atcep1-3 atcep3-2 double mutant have reduced levels of EXT deposition, suggesting that the defects in RH elongation are linked to alterations in EXT processing and accumulation. Our results support the involvement of AtCEPs in controlling RH polar growth through EXT processing and insolubilization at the cell wall.

Distinct guard cell-specific remodeling of chromatin accessibility during abscisic acid- and CO2-dependent stomatal regulation

In plants, epidermal guard cells integrate and respond to numerous environmental signals to control stomatal pore apertures, thereby regulating gas exchange. Chromatin structure controls transcription factor (TF) access to the genome, but whether large-scale chromatin remodeling occurs in guard cells during stomatal movements, and in response to the hormone abscisic acid (ABA) in general, remains unknown. Here, we isolate guard cell nuclei from Arabidopsis thaliana plants to examine whether the physiological signals, ABA and CO2 (carbon dioxide), regulate guard cell chromatin during stomatal movements. Our cell type-specific analyses uncover patterns of chromatin accessibility specific to guard cells and define cis-regulatory sequences supporting guard cell-specific gene expression. We find that ABA triggers extensive and dynamic chromatin remodeling in guard cells, roots, and mesophyll cells with clear patterns of cell type specificity. DNA motif analyses uncover binding sites for distinct TFs enriched in ABA-induced and ABA-repressed chromatin. We identify the Abscisic Acid Response Element (ABRE) Binding Factor (ABF) bZIP-type TFs that are required for ABA-triggered chromatin opening in guard cells and roots and implicate the inhibition of a clade of bHLH-type TFs in controlling ABA-repressed chromatin. Moreover, we demonstrate that ABA and CO2 induce distinct programs of chromatin remodeling, whereby elevated atmospheric CO2 had only minimal impact on chromatin dynamics. We provide insight into the control of guard cell chromatin dynamics and propose that ABA-induced chromatin remodeling primes the genome for abiotic stress resistance.

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.

Dynamics of chromatin accessibility and genome wide control of desiccation tolerance in the resurrection plant Haberlea rhodopensis

BACKGROUND

Drought is one of the main consequences of global climate change and this problem is expected to intensify in the future. Resurrection plants evolved the ability to withstand the negative impact of long periods of almost complete desiccation and to recover at rewatering. In this respect, many physiological, transcriptomic, proteomic and genomic investigations have been performed in recent years, however, few epigenetic control studies have been performed on these valuable desiccation-tolerant plants so far.

RESULTS

In the present study, for the first time for resurrection plants we provide evidences about the differential chromatin accessibility of Haberlea rhodopensis during desiccation stress by ATAC-seq (Assay for Transposase Accessible Chromatin with high-throughput sequencing). Based on gene similarity between species, we used the available genome of the closely related resurrection plant Dorcoceras hygrometricum to identify approximately nine hundred transposase hypersensitive sites (THSs) in H. rhodopensis. The majority of them corresponds to proximal and distal regulatory elements of different genes involved in photosynthesis, carbon metabolism, synthesis of secondary metabolites, cell signalling and transcriptional regulation, cell growth, cell wall, stomata conditioning, chaperons, oxidative stress, autophagy and others. Various types of binding motifs recognized by several families of transcription factors have been enriched from the THSs found in different stages of drought. Further, we used the previously published RNA-seq data from H. rhodopensis to evaluate the expression of transcription factors putatively interacting with the enriched motifs, and the potential correlation between the identified THS and the expression of their corresponding genes.

CONCLUSIONS

These results provide a blueprint for investigating the epigenetic regulation of desiccation tolerance in resurrection plant H. rhodopensis and comparative genomics between resurrection and non-resurrection species with available genome information.

Changes in epigenetic features in legumes under abiotic stresses

Legume crops are rich in nutritional value for human and livestock consumption. With global climate change, developing stress-resilient crops is crucial for ensuring global food security. Because of their nitrogen-fixing ability, legumes are also important for sustainable agriculture. Various abiotic stresses, such as salt, drought, and elevated temperatures, are known to adversely affect legume production. The responses of plants to abiotic stresses involve complicated cellular processes including stress hormone signaling, metabolic adjustments, and transcriptional regulations. Epigenetic mechanisms play a key role in regulating gene expressions at both transcriptional and posttranscriptional levels. Increasing evidence suggests the importance of epigenetic regulations of abiotic stress responses in legumes, and recent investigations have extended the scope to the epigenomic level using next-generation sequencing technologies. In this review, the current knowledge on the involvement of epigenetic features, including DNA methylation, histone modification, and noncoding RNAs, in abiotic stress responses in legumes is summarized and discussed. Since most of the available information focuses on a single aspect of these epigenetic features, integrative analyses involving omics data in multiple layers are needed for a better understanding of the dynamic chromatin statuses and their roles in transcriptional regulation. The inheritability of epigenetic modifications should also be assessed in future studies for their applications in improving stress tolerance in legumes through the stable epigenetic optimization of gene expressions.

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.

Differences in transcription initiation directionality underlie distinctions between plants and animals in chromatin modification patterns at genes and cis-regulatory elements

Transcriptional initiation is among the first regulated steps controlling eukaryotic gene expression. High-throughput profiling of fungal and animal genomes has revealed that RNA Polymerase II (Pol II) often initiates transcription in both directions at the promoter transcription start site (TSS), but generally only elongates productively into the gene body. Additionally, Pol II can initiate transcription in both directions at cis-regulatory elements (CREs) such as enhancers. These bidirectional Pol II initiation events can be observed directly with methods that capture nascent transcripts, and they are also revealed indirectly by the presence of transcription-associated histone modifications on both sides of the TSS or CRE. Previous studies have shown that nascent RNAs and transcription-associated histone modifications in the model plant Arabidopsis thaliana accumulate mainly in the gene body, suggesting that transcription does not initiate widely in the upstream direction from genes in this plant. We compared transcription-associated histone modifications and nascent transcripts at both TSSs and CREs in Arabidopsis thaliana, Drosophila melanogaster, and Homo sapiens. Our results provide evidence for mostly unidirectional Pol II initiation at both promoters and gene-proximal CREs of Arabidopsis thaliana, whereas bidirectional transcription initiation is observed widely at promoters in both Drosophila melanogaster and Homo sapiens, as well as CREs in Drosophila. Furthermore, the distribution of transcription-associated histone modifications around TSSs in the Oryza sativa (rice) and Glycine max (soybean) genomes suggests that unidirectional transcription initiation is the norm in these genomes as well. These results suggest that there are fundamental differences in transcriptional initiation directionality between flowering plant and metazoan genomes, which are manifested as distinct patterns of chromatin modifications around RNA polymerase initiation sites.

The eINTACT method for studying nuclear changes in host plant cells targeted by bacterial effectors in native infection contexts

Type-III effector proteins are major virulence determinants that most gram-negative bacteria inject into host cells to manipulate cellular processes for infection. Because effector-targeted cells are embedded and underrepresented in infected plant tissues, it is technically challenging to isolate them for focused studies of effector-induced cellular changes. This protocol describes a novel technique, effector-inducible isolation of nuclei tagged in specific cell types (eINTACT), for isolating biotin-labeled nuclei from Arabidopsis plant cells that have received Xanthomonas bacterial effectors by using streptavidin-coated magnetic beads. This protocol is an extension of the existing Nature Protocols Protocol of the INTACT method for the affinity-based purification of nuclei of specific cell types in the context of developmental biology. In a phytopathology scenario, our protocol addresses how to obtain eINTACT transgenic lines and compatible bacterial mutants, verify the eINTACT system and purify nuclei of bacterial effector-recipient cells from infected tissues. Differential analyses of purified nuclei from plants infected by bacteria expressing the effector of interest and those from plants infected by effector-deletion bacterial mutants will reveal the effector-dependent nuclear changes in targeted host cells. Provided that the eINTACT system is available, the infection experiment takes 5 d, and the procedures, from collecting bacteria-infected leaves to obtaining nuclei of effector-targeted cells, can be completed in 4 h. eINTACT is a unique method for isolating high-quality nuclei from bacterial effector-targeted host cells in native infection contexts. This method is adaptable to study the functions of type-III effectors from numerous gram-negative bacteria in host plants that are amenable to transformation.

Distinct guard cell specific remodeling of chromatin accessibility during abscisic acid and CO2 dependent stomatal regulation

In plants, epidermal guard cells integrate and respond to numerous environmental signals to control stomatal pore apertures thereby regulating gas exchange. Chromatin structure controls transcription factor access to the genome, but whether large-scale chromatin remodeling occurs in guard cells during stomatal movements, and in response to the hormone abscisic acid (ABA) in general, remain unknown. Here we isolate guard cell nuclei from Arabidopsis thaliana plants to examine whether the physiological signals, ABA and CO2, regulate guard cell chromatin during stomatal movements. Our cell type specific analyses uncover patterns of chromatin accessibility specific to guard cells and define novel cis-regulatory sequences supporting guard cell specific gene expression. We find that ABA triggers extensive and dynamic chromatin remodeling in guard cells, roots, and mesophyll cells with clear patterns of cell-type specificity. DNA motif analyses uncover binding sites for distinct transcription factors enriched in ABA-induced and ABA-repressed chromatin. We identify the ABF/AREB bZIP-type transcription factors that are required for ABA-triggered chromatin opening in guard cells and implicate the inhibition of a set of bHLH-type transcription factors in controlling ABA-repressed chromatin. Moreover, we demonstrate that ABA and CO2 induce distinct programs of chromatin remodeling. We provide insight into the control of guard cell chromatin dynamics and propose that ABA-induced chromatin remodeling primes the genome for abiotic stress resistance.

Engineering custom morpho- and chemotypes of Populus for sustainable production of biofuels, bioproducts, and biomaterials

Humans have been modifying plant traits for thousands of years, first through selection (i.e., domestication) then modern breeding, and in the last 30 years, through biotechnology. These modifications have resulted in increased yield, more efficient agronomic practices, and enhanced quality traits. Precision knowledge of gene regulation and function through high-resolution single-cell omics technologies, coupled with the ability to engineer plant genomes at the DNA sequence, chromatin accessibility, and gene expression levels, can enable engineering of complex and complementary traits at the biosystem level. Populus spp., the primary genetic model system for woody perennials, are among the fastest growing trees in temperate zones and are important for both carbon sequestration and global carbon cycling. Ample genomic and transcriptomic resources for poplar are available including emerging single-cell omics datasets. To expand use of poplar outside of valorization of woody biomass, chassis with novel morphotypes in which stem branching and tree height are modified can be fabricated thereby leading to trees with altered leaf to wood ratios. These morphotypes can then be engineered into customized chemotypes that produce high value biofuels, bioproducts, and biomaterials not only in specific organs but also in a cell-type-specific manner. For example, the recent discovery of triterpene production in poplar leaf trichomes can be exploited using cell-type specific regulatory sequences to synthesize high value terpenes such as the jet fuel precursor bisabolene specifically in the trichomes. By spatially and temporally controlling expression, not only can pools of abundant precursors be exploited but engineered molecules can be sequestered in discrete cell structures in the leaf. The structural diversity of the hemicellulose xylan is a barrier to fully utilizing lignocellulose in biomaterial production and by leveraging cell-type-specific omics data, cell wall composition can be modified in a tailored and targeted specific manner to generate poplar wood with novel chemical features that are amenable for processing or advanced manufacturing. Precision engineering poplar as a multi-purpose sustainable feedstock highlights how genome engineering can be used to re-imagine a crop species.

Transcription factors KANADI 1, MYB DOMAIN PROTEIN 44, and PHYTOCHROME INTERACTING FACTOR 4 regulate long intergenic noncoding RNAs expressed in Arabidopsis roots

Thousands of long intergenic noncoding RNAs (lincRNAs) have been identified in plant genomes. While some lincRNAs have been characterized as important regulators in different biological processes, little is known about the transcriptional regulation for most plant lincRNAs. Through the integration of 8 annotation resources, we defined 6,599 high-confidence lincRNA loci in Arabidopsis (Arabidopsis thaliana). For lincRNAs belonging to different evolutionary age categories, we identified major differences in sequence and chromatin features, as well as in the level of conservation and purifying selection acting during evolution. Spatiotemporal gene expression profiles combined with transcription factor (TF) chromatin immunoprecipitation (ChIP) data were used to construct a TF-lincRNA regulatory network containing 2,659 lincRNAs and 15,686 interactions. We found that properties characterizing lincRNA expression, conservation, and regulation differ between plants and animals. Experimental validation confirmed the role of 3 TFs, KANADI 1, MYB DOMAIN PROTEIN 44, and PHYTOCHROME INTERACTING FACTOR 4, as key regulators controlling root-specific lincRNA expression, demonstrating the predictive power of our network. Furthermore, we identified 58 lincRNAs, regulated by these TFs, showing strong root cell type-specific expression or chromatin accessibility, which are linked with genome-wide association studies genetic associations related to root system development and growth. The multilevel genome-wide characterization covering chromatin state information, promoter conservation, and chromatin immunoprecipitation-based TF binding, for all detectable lincRNAs across 769 expression samples, permits rapidly defining the biological context and relevance of Arabidopsis lincRNAs through regulatory networks.

Integrating ATAC-seq and RNA-seq to identify differentially expressed genes with chromatin-accessible changes during photosynthetic establishment in Populus leaves

Leaves are the primary photosynthetic organs, providing essential substances for tree growth. It is important to obtain an anatomical understanding and regulatory network analysis of leaf development. Here, we studied leaf development in Populus Nanlin895 along a development gradient from the newly emerged leaf from the shoot apex to the sixth leaf (L1 to L6) using anatomical observations and RNA-seq analysis. It indicated that mesophyll cells possess obvious vascular, palisade, and spongy tissue with distinct intercellular spaces after L3. Additionally, vacuoles fuse while epidermal cells expand to form pavement cells. RNA-seq analysis indicated that genes highly expressed in L1 and L2 were related to cell division and differentiation, while those highly expressed in L3 were enriched in photosynthesis. Therefore, we selected L1 and L3 to integrate ATAC-seq and RNA-seq and identified 735 differentially expressed genes (DEGs) with changes in chromatin accessibility regions within their promoters, of which 87 were transcription factors (TFs), such as ABI3VP1, AP-EREBP, MYB, NAC, and GRF. Motif enrichment analysis revealed potential regulatory functions for the DEGs through upstream TFs including TCP, bZIP, HD-ZIP, Dof, BBR-BPC, and MYB. Overall, our research provides a potential molecular foundation for regulatory network exploration in leaf development during photosynthesis establishment.

Spatially resolved gene regulatory networks in Asian rice (Oryza sativa cv. Nipponbare) leaves

Transcriptome profiles in plants are heterogenous at every level of morphological organization. Even within organs, cells of the same type can have different patterns of gene expression depending on where they are positioned within tissues. This heterogeneity is associated with non-uniform distribution of biological processes within organs. The regulatory mechanisms that establish and sustain the spatial heterogeneity are unknown. Here, we identify regulatory modules that support functional specialization of different parts of Oryza sativa cv. Nipponbare leaves by leveraging transcriptome data, transcription factor binding motifs and global gene regulatory network prediction algorithms. We generated a global gene regulatory network in which we identified six regulatory modules that were active in different parts of the leaf. The regulatory modules were enriched for genes involved in spatially relevant biological processes, such as cell wall deposition, environmental sensing and photosynthesis. Strikingly, more than 86.9% of genes in the network were regulated by members of only five transcription factor families. We also generated targeted regulatory networks for the large MYB and bZIP/bHLH families to identify interactions that were masked in the global prediction. This analysis will provide a baseline for future single cell and array-based spatial transcriptome studies and for studying responses to environmental stress and demonstrates the extent to which seven coarse spatial transcriptome analysis can provide insight into the regulatory mechanisms supporting functional specialization within leaves.

Exploring plant cis-regulatory elements at single-cell resolution: overcoming biological and computational challenges to advance plant research

Cis-regulatory elements (CREs) are important sequences for gene expression and for plant biological processes such as development, evolution, domestication, and stress response. However, studying CREs in plant genomes has been challenging. The totipotent nature of plant cells, coupled with the inability to maintain plant cell types in culture and the inherent technical challenges posed by the cell wall has limited our understanding of how plant cell types acquire and maintain their identities and respond to the environment via CRE usage. Advances in single-cell epigenomics have revolutionized the field of identifying cell-type-specific CREs. These new technologies have the potential to significantly advance our understanding of plant CRE biology, and shed light on how the regulatory genome gives rise to diverse plant phenomena. However, there are significant biological and computational challenges associated with analyzing single-cell epigenomic datasets. In this review, we discuss the historical and foundational underpinnings of plant single-cell research, challenges, and common pitfalls in the analysis of plant single-cell epigenomic data, and highlight biological challenges unique to plants. Additionally, we discuss how the application of single-cell epigenomic data in various contexts stands to transform our understanding of the importance of CREs in plant genomes.

Low temperature-induced regulatory network rewiring via WRKY regulators during banana peel browning

Banana (Musa spp.) fruits, as typical tropical fruits, are cold sensitive, and lower temperatures can disrupt cellular compartmentalization and lead to severe browning. How tropical fruits respond to low temperature compared to the cold response mechanisms of model plants remains unknown. Here, we systematically characterized the changes in chromatin accessibility, histone modifications, distal cis-regulatory elements, transcription factor binding, and gene expression levels in banana peels in response to low temperature. Dynamic patterns of cold-induced transcripts were generally accompanied by concordant chromatin accessibility and histone modification changes. These upregulated genes were enriched for WRKY binding sites in their promoters and/or active enhancers. Compared to banana peel at room temperature, large amounts of banana WRKYs were specifically induced by cold and mediated enhancer-promoter interactions regulating critical browning pathways, including phospholipid degradation, oxidation, and cold tolerance. This hypothesis was supported by DNA affinity purification sequencing, luciferase reporter assays, and transient expression assay. Together, our findings highlight widespread transcriptional reprogramming via WRKYs during banana peel browning at low temperature and provide an extensive resource for studying gene regulation in tropical plants in response to cold stress, as well as potential targets for improving cold tolerance and shelf life of tropical fruits.

Gene expression and expression quantitative trait loci analyses uncover natural variations underlying the improvement of important agronomic traits during modern maize breeding

Maize (Zea mays L.) is a major staple crop worldwide, and during modern maize breeding, cultivars with increased tolerance to high-density planting and higher yield per plant have contributed significantly to the increased yield per unit land area. Systematically identifying key agronomic traits and their associated genomic changes during modern maize breeding remains a significant challenge because of the complexity of genetic regulation and the interactions of the various agronomic traits, with most of them being controlled by numerous small-effect quantitative trait loci (QTLs). Here, we performed phenotypic and gene expression analyses for a set of 137 elite inbred lines of maize from different breeding eras in China. We found four yield-related traits are significantly improved during modern maize breeding. Through gene-clustering analyses, we identified four groups of expressed genes with distinct trends of expression pattern change across the historical breeding eras. In combination with weighted gene co-expression network analysis, we identified several candidate genes regulating various plant architecture- and yield-related agronomic traits, such as ZmARF16, ZmARF34, ZmTCP40, ZmPIN7, ZmPYL10, ZmJMJ10, ZmARF1, ZmSWEET15b, ZmGLN6 and Zm00001d019150. Further, by combining expression quantitative trait loci (eQTLs) analyses, correlation coefficient analyses and population genetics, we identified a set of candidate genes that might have been under selection and contributed to the genetic improvement of various agronomic traits during modern maize breeding, including a number of known key regulators of plant architecture, flowering time and yield-related traits, such as ZmPIF3.3, ZAG1, ZFL2 and ZmBES1. Lastly, we validated the functional variations in GL15, ZmPHYB2 and ZmPYL10 that influence kernel row number, flowering time, plant height and ear height, respectively. Our results demonstrates the effectiveness of our combined approaches for uncovering key candidate regulatory genes and functional variation underlying the improvement of important agronomic traits during modern maize breeding, and provide a valuable genetic resource for the molecular breeding of maize cultivars with tolerance for high-density planting.

Genome-wide identification of accessible chromatin regions by ATAC-seq upon induction of the transcription factor bZIP11 in Arabidopsis

Basic leucine zipper 11 (bZIP11) is a transcription factor that is activated under low energy conditions in plants and plays a crucial role in enabling plants to adapt to starvation situations. Although previous results indicate that bZIP11 regulates chromatin accessibility based on evidence obtained from single genomic loci, to what extent this transcription factor regulates the chromatin landscape at the whole genome level remains unknown. Here we addressed this by performing an ATAC-seq (Assay for Transposase-Accessible Chromatin with high-throughput sequencing) on Arabidopsis thaliana (Arabidopsis) leaf protoplasts to obtain a profile of chromatin patterning in response upon bZIP11 induction. We identified, on average, 10,000 differentially accessible regions upon bZIP11 induction, corresponding to over 8,420 different genes out of the 25,000 genes present in the Arabidopsis genome. Our study provides a resource for understanding how bZIP11 regulates the genome at the chromatin level and provides an example of the impact of a single transcription factor on a whole plant genome.

Plant Promoters and Terminators for High-Precision Bioengineering

High-precision bioengineering and synthetic biology require fine-tuning gene expression at both transcriptional and posttranscriptional levels. Gene transcription is tightly regulated by promoters and terminators. Promoters determine the timing, tissues and cells, and levels of the expression of genes. Terminators mediate transcription termination of genes and affect mRNA levels posttranscriptionally, e.g., the 3'-end processing, stability, translation efficiency, and nuclear to cytoplasmic export of mRNAs. The promoter and terminator combination affects gene expression. In the present article, we review the function and features of plant core promoters, proximal and distal promoters, and terminators, and their effects on and benchmarking strategies for regulating gene expression.

Genome-wide mapping of DNase I hypersensitive sites in pineapple leaves

Pineapple [Ananas comosus (L.) Merr.] is the most economically important crop possessing crassulacean acid metabolism (CAM) photosynthesis which has a higher water use efficiency by control of nocturnal opening and diurnal closure of stomata. To provide novel insights into the diel regulatory landscape in pineapple leaves, we performed genome-wide mapping of DNase I hypersensitive sites (DHSs) in pineapple leaves at day (2a.m.) and night (10a.m.) using a simplified DNase-seq method. As a result, totally 33340 and 28753 DHSs were found in green-tip tissue, and 29597 and 40068 were identified in white-base tissue at 2a.m. and 10a.m., respectively. We observed that majority of the pineapple genes occupied less than two DHSs with length shorter than 1 kb, and the promotor DHSs showed a proximal trend to the transcription start site (>77% promotor DHSs within 1 kb). In addition, more intergenic DHSs were identified around transcription factors or transcription co-regulators (TFs/TCs) than other functional genes, indicating complex regulatory contexts around TFs/TCs. Through combined analysis of tissue preferential DHSs and genes, we respectively found 839 and 888 coordinately changed genes in green-tip at 2a.m. and 10a.m. (AcG2 and AcG10). Furthermore, AcG2-specific, AcG10-specific and common accessible DHSs were dissected from the total photosynthetic preferential DHSs, and the regulatory networks indicated dynamic regulations with multiple cis-regulatory elements occurred to genes preferentially expressed in photosynthetic tissues. Interestingly, binding motifs of several cycling TFs were identified in the DHSs of key CAM genes, revealing a circadian regulation to CAM coordinately diurnal expression. Our results provide a chromatin regulatory landscape in pineapple leaves during the day and night. This will provide important information to assist with deciphering the circadian regulation of CAM photosynthesis.