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

Background: Pre-exposing plants to abiotic stresses can induce stress memory, which is crucial for adapting to subsequent stress exposure. Although numerous genes involved in salt stress response have been identified, the understanding of memory responses to salt stress remains limited. Results: In this study, we conducted physiological and transcriptional assays on maize plants subjected to recurrent salt stress to characterize salt stress memory. During the second exposure to salt stress, the plants exhibited enhanced salt resistance, as evidenced by increased proline content and higher POD and SOD activity, along with decreased MDA content, indicative of physiological memory behavior. Transcriptional analysis revealed fewer differentially expressed genes and variations in response processes during the second exposure compared to the first, indicative of transcriptional memory behavior. A total of 2,213 salt stress memory genes (SMGs) were identified and categorized into four response patterns. The most prominent group of SMGs consisted of genes with elevated expression during the first exposure to salt stress but reduced expression after recurrent exposure to salt stress, or vice versa ([+ / −] or [− / +]), indicating that a revised response is a crucial process in plant stress memory. Furthermore, nine transcription factors (TFs) (WRKY40, WRKY46, WRKY53, WRKY18, WRKY33, WRKY70, MYB15, KNAT7, and WRKY54) were identified as crucial factors related to salt stress memory. These TFs regulate over 53% of SMGs, underscoring their potential significance in salt stress memory. Conclusions: Our study demonstrates that maize can develop salt stress memory, and the genes identified here will aid in the genetic improvement of maize and other crops.

Cold stress significantly challenges cotton growth and productivity, yet the genetic and molecular mechanisms underlying cold tolerance remain poorly understood. We employed RNA-seq and iterative weighted gene co-expression network analysis (WGCNA) to investigate gene and transposable element (TE) expression changes at six cold stress time points (0 h, 2 h, 4 h, 6 h, 12 h, 24 h). Thousands of differentially expressed genes (DEGs) were identified, exhibiting time-specific patterns that highlight a phase-dependent transcriptional response. While the A and D subgenomes contributed comparably to DEG numbers, numerous homeologous gene pairs showed differential expression, indicating regulatory divergence. Iterative WGCNA uncovered 125 gene co-expression modules, with some enriched in specific chromosomes or chromosomal regions, suggesting localized regulatory hotspots for cold stress response. Notably, transcription factors, including MYB73, ERF017, MYB30, and OBP1, emerged as central regulators within these modules. Analysis of 11 plant hormone-related genes revealed dynamic expression, with ethylene (ETH) and cytokinins (CK) playing significant roles in stress-responsive pathways. Furthermore, we documented over 15,000 expressed TEs, with differentially expressed TEs forming five distinct clusters. TE families, such as LTR/Copia, demonstrated significant enrichment in these expression clusters, suggesting their potential role as modulators of gene expression under cold stress. These findings provide valuable insights into the complex regulatory networks underlying cold stress response in cotton, highlighting key molecular components involved in cold stress regulation. This study provides potential genetic targets for breeding strategies aimed at enhancing cold tolerance in cotton.

Background: The adaptation of plants to cold stress involves changes in gene expression profiles that are associated with epigenetic regulation. Although the three-dimensional (3D) genome architecture is considered an important epigenetic regulator, the role of 3D genome organization in the cold stress response remains unclear. Results: In this study, we developed high-resolution 3D genomic maps using control and cold-treated leaf tissue of the model plant Brachypodium distachyon using Hi-C to determine how cold stress affects the 3D genome architecture. We generated ~ 1.5 kb resolution chromatin interaction maps and showed that cold stress disrupts different levels of chromosome organization, including A/B compartment transition, a reduction in chromatin compartmentalization and the size of topologically associating domains (TADs), and loss of long-range chromatin loops. Integrating RNA-seq information, we identified cold-response genes and revealed that transcription was largely unaffected by the A/B compartment transition. The cold-response genes were predominantly localized in compartment A. In contrast, transcriptional changes are required for TAD reorganization. We demonstrated that dynamic TAD events were associated with H3K27me3 and H3K27ac state alterations. Moreover, a loss of chromatin looping, rather than a gain of looping, coincides with alterations in gene expression, indicating that chromatin loop disruption may play a more important role than loop formation in the cold-stress response. Conclusions: Our study highlights the multiscale 3D genome reprogramming that occurs during cold stress and expands our knowledge of the mechanisms underlying transcriptional regulation in response to cold stress in plants.

Cassava, a crucial tropical crop, faces challenges from cold stress, necessitating an exploration of its molecular response. Here, we investigated the role of DNA methylation in moderating the response to moderate cold stress (10 °C) in cassava. Using whole-genome bisulfite sequencing, we examined DNA methylation patterns in leaf blades and petioles under control conditions, 5 h, and 48 h of cold stress. Tissue-specific responses were observed, with leaf blades exhibiting subtle changes, while petioles displayed a pronounced decrease in methylation levels under cold stress. We identified cold stress-induced differentially methylated regions (DMRs) that demonstrated both tissue and treatment specificity. Importantly, these DMRs were enriched in genes with altered expression, implying functional relevance. The cold-response transcription factor ERF105 associated with DMRs emerged as a significant and conserved regulator across tissues and treatments. Furthermore, we investigated DNA methylation dynamics in transposable elements, emphasizing the sensitivity of MITEs with bHLH binding motifs to cold stress. These findings provide insights into the epigenetic regulation of response to cold stress in cassava, contributing to an understanding of the molecular mechanisms underlying stress adaptation in this tropical plant.