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.

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.