Considerations in the analysis of plant chromatin accessibility data

Evolutionary Dynamics of Chromatin Structure and Duplicate Gene Expression in Diploid and Allopolyploid Cotton

Polyploidy is a prominent mechanism of plant speciation and adaptation, yet the mechanistic understandings of duplicated gene regulation remain elusive. Chromatin structure dynamics are suggested to govern gene regulatory control. Here, we characterized genome-wide nucleosome organization and chromatin accessibility in allotetraploid cotton, Gossypium hirsutum (AADD, 2n = 4X = 52), relative to its two diploid parents (AA or DD genome) and their synthetic diploid hybrid (AD), using DNS-seq. The larger A-genome exhibited wider average nucleosome spacing in diploids, and this intergenomic difference diminished in the allopolyploid but not hybrid. Allopolyploidization also exhibited increased accessibility at promoters genome-wide and synchronized cis-regulatory motifs between subgenomes. A prominent cis-acting control was inferred for chromatin dynamics and demonstrated by transposable element removal from promoters. Linking accessibility to gene expression patterns, we found distinct regulatory effects for hybridization and later allopolyploid stages, including nuanced establishment of homoeolog expression bias and expression level dominance. Histone gene expression and nucleosome organization are coordinated through chromatin accessibility. Our study demonstrates the capability to track high-resolution chromatin structure dynamics and reveals their role in the evolution of cis-regulatory landscapes and duplicate gene expression in polyploids, illuminating regulatory ties to subgenomic asymmetry and dominance.

A review of artificial intelligence-assisted omics techniques in plant defense: current trends and future directions

Plants intricately deploy defense systems to counter diverse biotic and abiotic stresses. Omics technologies, spanning genomics, transcriptomics, proteomics, and metabolomics, have revolutionized the exploration of plant defense mechanisms, unraveling molecular intricacies in response to various stressors. However, the complexity and scale of omics data necessitate sophisticated analytical tools for meaningful insights. This review delves into the application of artificial intelligence algorithms, particularly machine learning and deep learning, as promising approaches for deciphering complex omics data in plant defense research. The overview encompasses key omics techniques and addresses the challenges and limitations inherent in current AI-assisted omics approaches. Moreover, it contemplates potential future directions in this dynamic field. In summary, AI-assisted omics techniques present a robust toolkit, enabling a profound understanding of the molecular foundations of plant defense and paving the way for more effective crop protection strategies amidst climate change and emerging diseases.

Selective penetration of fullerenol through pea seed coats mitigates osmosis-repressed germination via chromatin remodeling and transcriptional reprograming

BACKGROUND

The alteration of chromatin accessibility plays an important role in plant responses to abiotic stress. Carbon-based nanomaterials (CBNMs) have attracted increasing interest in agriculture due to their potential impact on crop productivity, showcasing effects on plant biological processes at transcriptional levels; however, their impact on chromatin accessibility remains unknown.

RESULTS

This study found that fullerenol can penetrate the seed coat of pea to mitigate the reduction of seed germination caused by osmotic stress. RNA sequencing (RNA-seq) revealed that the application of fullerenol caused the high expression of genes related to oxidoreduction to return to a normal level. Assay for transposase accessible chromatin sequencing (ATAC-seq) confirmed that fullerenol application reduced the overall levels of chromatin accessibility of numerous genes, including those related to environmental signaling, transcriptional regulation, and metabolism.

CONCLUSION

This study suggests that fullerenol alleviates osmotic stress on various fronts, encompassing antioxidant, transcriptional, and epigenetic levels. This advances knowledge of the working mechanism of this nanomaterial within plant cells. © 2024 Society of Chemical Industry.

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.

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.

Optimization of ATAC-seq in wheat seedling roots using INTACT-isolated nuclei

BACKGROUND

The genetic information contained in the genome of an organism is organized in genes and regulatory elements that control gene expression. The genomes of multiple plants species have already been sequenced and the gene repertory have been annotated, however, cis-regulatory elements remain less characterized, limiting our understanding of genome functionality. These elements act as open platforms for recruiting both positive- and negative-acting transcription factors, and as such, chromatin accessibility is an important signature for their identification.

RESULTS

In this work we developed a transgenic INTACT [isolation of nuclei tagged in specific cell types] system in tetraploid wheat for nuclei purifications. Then, we combined the INTACT system together with the assay for transposase-accessible chromatin with sequencing [ATAC-seq] to identify open chromatin regions in wheat root tip samples. Our ATAC-seq results showed a large enrichment of open chromatin regions in intergenic and promoter regions, which is expected for regulatory elements and that is similar to ATAC-seq results obtained in other plant species. In addition, root ATAC-seq peaks showed a significant overlap with a previously published ATAC-seq data from wheat leaf protoplast, indicating a high reproducibility between the two experiments and a large overlap between open chromatin regions in root and leaf tissues. Importantly, we observed overlap between ATAC-seq peaks and cis-regulatory elements that have been functionally validated in wheat, and a good correlation between normalized accessibility and gene expression levels.

CONCLUSIONS

We have developed and validated an INTACT system in tetraploid wheat that allows rapid and high-quality nuclei purification from root tips. Those nuclei were successfully used to performed ATAC-seq experiments that revealed open chromatin regions in the wheat genome that will be useful to identify cis-regulatory elements. The INTACT system presented here will facilitate the development of ATAC-seq datasets in other tissues, growth stages, and under different growing conditions to generate a more complete landscape of the accessible DNA regions in the wheat genome.

Machine learning bridges omics sciences and plant breeding

Some of the biological knowledge obtained from fundamental research will be implemented in applied plant breeding. To bridge basic research and breeding practice, machine learning (ML) holds great promise to translate biological knowledge and omics data into precision-designed plant breeding. Here, we review ML for multi-omics analysis in plants, including data dimensionality reduction, inference of gene-regulation networks, and gene discovery and prioritization. These applications will facilitate understanding trait regulation mechanisms and identifying target genes potentially applicable to knowledge-driven molecular design breeding. We also highlight applications of deep learning in plant phenomics and ML in genomic selection-assisted breeding, such as various ML algorithms that model the correlations among genotypes (genes), phenotypes (traits), and environments, to ultimately achieve data-driven genomic design breeding.

Chromatin dynamics associated with seed desiccation tolerance/sensitivity at early germination in Medicago truncatula

Desiccation tolerance (DT) has contributed greatly to the adaptation of land plants to severe water-deficient conditions. DT is mostly observed in reproductive parts in flowering plants such as seeds. The seed DT is lost at early post germination stage but is temporally re-inducible in 1 mm radicles during the so-called DT window following a PEG treatment before being permanently silenced in 5 mm radicles of germinating seeds. The molecular mechanisms that activate/reactivate/silence DT in developing and germinating seeds have not yet been elucidated. Here, we analyzed chromatin dynamics related to re-inducibility of DT before and after the DT window at early germination in Medicago truncatula radicles to determine if DT-associated genes were transcriptionally regulated at the chromatin levels. Comparative transcriptome analysis of these radicles identified 948 genes as DT re-induction-related genes, positively correlated with DT re-induction. ATAC-Seq analyses revealed that the chromatin state of genomic regions containing these genes was clearly modulated by PEG treatment and affected by growth stages with opened chromatin in 1 mm radicles with PEG (R1P); intermediate openness in 1 mm radicles without PEG (R1); and condensed chromatin in 5 mm radicles without PEG (R5). In contrast, we also showed that the 103 genes negatively correlated with the re-induction of DT did not show any transcriptional regulation at the chromatin level. Additionally, ChIP-Seq analyses for repressive marks H2AK119ub and H3K27me3 detected a prominent signal of H3K27me3 on the DT re-induction-related gene sequences at R5 but not in R1 and R1P. Moreover, no clear H2AK119ub marks was observed on the DT re-induction-related gene sequences at both developmental radicle stages, suggesting that silencing of DT process after germination will be mainly due to H3K27me3 marks by the action of the PRC2 complex, without involvement of PRC1 complex. The dynamic of chromatin changes associated with H3K27me3 were also confirmed on seed-specific genes encoding potential DT-related proteins such as LEAs, oleosins and transcriptional factors. However, several transcriptional factors did not show a clear link between their decrease of chromatin openness and H3K27me3 levels, suggesting that their accessibility may also be regulated by additional factors, such as other histone modifications. Finally, in order to make these comprehensive genome-wide analyses of transcript and chromatin dynamics useful to the scientific community working on early germination and DT, we generated a dedicated genome browser containing all these data and publicly available at https://iris.angers.inrae.fr/mtseedepiatlas/jbrowse/?data=Mtruncatula.

Genome-wide chromatin accessibility analysis unveils open chromatin convergent evolution during polyploidization in cotton

Allopolyploidization, resulting in divergent genomes in the same cell, is believed to trigger a “genome shock”, leading to broad genetic and epigenetic changes. However, little is understood about chromatin and gene-expression dynamics as underlying driving forces during allopolyploidization. Here, we examined the genome-wide DNase I-hypersensitive site (DHS) and its variations in domesticated allotetraploid cotton ( Gossypium hirsutum and Gossypium barbadense , AADD) and its extant AA ( Gossypium arboreum ) and DD ( Gossypium raimondii ) progenitors. We observed distinct DHS distributions between G. arboreum and G. raimondii . In contrast, the DHSs of the two subgenomes of G. hirsutum and G. barbadense showed a convergent distribution. This convergent distribution of DHS was also present in the wild allotetraploids Gossypium darwinii and G. hirsutum var. yucatanense , but absent from a resynthesized hybrid of G. arboreum and G. raimondii , suggesting that it may be a common feature in polyploids, and not a consequence of domestication after polyploidization. We revealed that putative cis -regulatory elements (CREs) derived from polyploidization-related DHSs were dominated by several families, including Dof, ERF48, and BPC1. Strikingly, 56.6% of polyploidization-related DHSs were derived from transposable elements (TEs). Moreover, we observed positive correlations between DHS accessibility and the histone marks H3K4me3, H3K27me3, H3K36me3, H3K27ac, and H3K9ac, indicating that coordinated interplay among histone modifications, TEs, and CREs drives the DHS landscape dynamics under polyploidization. Collectively, these findings advance our understanding of the regulatory architecture in plants and underscore the complexity of regulome evolution during polyploidization.

Quality control and evaluation of plant epigenomics data

Epigenomics is the study of molecular signatures associated with discrete regions within genomes, many of which are important for a wide range of nuclear processes. The ability to profile the epigenomic landscape associated with genes, repetitive regions, transposons, transcription, differential expression, cis-regulatory elements, and 3D chromatin interactions has vastly improved our understanding of plant genomes. However, many epigenomic and single-cell genomic assays are challenging to perform in plants, leading to a wide range of data quality issues; thus, the data require rigorous evaluation prior to downstream analyses and interpretation. In this commentary, we provide considerations for the evaluation of plant epigenomics and single-cell genomics data quality with the aim of improving the quality and utility of studies using those data across diverse plant species.

Genome accessibility dynamics in response to phosphate limitation is controlled by the PHR1 family of transcription factors in Arabidopsis

As phosphorus is one of the most limiting nutrients in many natural and agricultural ecosystems, plants have evolved strategies that cope with its scarcity. Genetic approaches have facilitated the identification of several molecular elements that regulate the phosphate (Pi) starvation response (PSR) of plants, including the master regulator of the transcriptional response to phosphate starvation PHOSPHATE STARVATION RESPONSE1 (PHR1). However, the chromatin modifications underlying the plant transcriptional response to phosphate scarcity remain largely unknown. Here, we present a detailed analysis of changes in chromatin accessibility during phosphate starvation in Arabidopsis thaliana root cells. Root cells undergo a genome-wide remodeling of chromatin accessibility in response to Pi starvation that is often associated with changes in the transcription of neighboring genes. Analysis of chromatin accessibility in the phr1 phl2 double mutant revealed that the transcription factors PHR1 and PHL2 play a key role in remodeling chromatin accessibility in response to Pi limitation. We also discovered that PHR1 and PHL2 play an important role in determining chromatin accessibility and the associated transcription of many genes under optimal Pi conditions, including genes involved in the PSR. We propose that a set of transcription factors directly activated by PHR1 in Pi-starved root cells trigger a second wave of epigenetic changes required for the transcriptional activation of the complete set of low-Pi-responsive genes.

The native cistrome and sequence motif families of the maize ear

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

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

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.

Network-based approaches for understanding gene regulation and function in plants

Expression reprogramming directed by transcription factors is a primary gene regulation underlying most aspects of the biology of any organism. Our views of how gene regulation is coordinated are dramatically changing thanks to the advent and constant improvement of high-throughput profiling and transcriptional network inference methods: from activities of individual genes to functional interactions across genes. These technical and analytical advances can reveal the topology of transcriptional networks in which hundreds of genes are hierarchically regulated by multiple transcription factors at systems level. Here we review the state of the art of experimental and computational methods used in plant biology research to obtain large-scale datasets and model transcriptional networks. Examples of direct use of these network models and perspectives on their limitations and future directions are also discussed.