Temperature influences the distribution, range, and phenology of plants. The key transcriptional activators of heat shock response in eukaryotes, the heat shock factors (HSFs), have undergone large-scale gene amplification in plants. While HSFs are central in heat stress responses, their role in the response to ambient temperature changes is less well understood. We show here that the warm ambient temperature transcriptome is dependent upon the HSFA1 clade of Arabidopsis HSFs, which cause a rapid and dynamic eviction of H2A.Z nucleosomes at target genes. A transcriptional cascade results in the activation of multiple downstream stress-responsive transcription factors, triggering large-scale changes to the transcriptome in response to elevated temperature. H2A.Z nucleosomes are enriched at temperature-responsive genes at non-inducible temperature, and thus likely confer inducibility of gene expression and higher responsive dynamics. We propose that the antagonistic effects of H2A.Z and HSF1 provide a mechanism to activate gene expression rapidly and precisely in response to temperature, while preventing leaky transcription in the absence of an activation signal.

Plants have significantly more transcription factor (TF) families than animals and fungi, and plant TF families tend to contain more genes; these expansions are linked to adaptation to environmental stressors. Many TF family members bind to similar or identical sequence motifs, such as G-boxes (CACGTG), so it is difficult to predict regulatory relationships. We determined that the flanking sequences near G-boxes help determine in vitro specificity but that this is insufficient to predict the transcription pattern of genes near G-boxes. Therefore, we constructed a gene regulatory network that identifies the set of bZIPs and bHLHs that are most predictive of the expression of genes downstream of perfect G-boxes. This network accurately predicts transcriptional patterns and reconstructs known regulatory subnetworks. Finally, we present Ara-BOX-cis (araboxcis.org), a Web site that provides interactive visualizations of the G-box regulatory network, a useful resource for generating predictions for gene regulatory relations.

Retrotransposons have played an important role in the evolution of host genomes 1 , 2 . Their impact is mainly deduced from the composition of DNA sequences that have been fixed over evolutionary time 2 . Such studies provide important “snapshots” reflecting the historical activities of transposons but do not predict current transposition potential. We previously reported Sequence-Independent Retrotransposon Trapping (SIRT) as a method that, by identification of extrachromosomal linear DNA (eclDNA), revealed the presence of active LTR retrotransposons in Arabidopsis 3 . However, SIRT cannot be applied to large and transposon-rich genomes, as found in crop plants. We have developed an alternative approach named ALE-seq ( a mplification of L TR of e clDNAs followed by seq uencing) for such situations. ALE-seq reveals sequences of 5’ LTRs of eclDNAs after two-step amplification: in vitro transcription and subsequent reverse transcription. Using ALE-seq in rice, we detected eclDNAs for a novel Copia family LTR retrotransposon, Go-on , which is activated by heat stress. Sequencing of rice accessions revealed that Go-on has preferentially accumulated in indica rice grown at higher temperatures. Furthermore, ALE-seq applied to tomato fruits identified a developmentally regulated Gypsy family of retrotransposons. A bioinformatic pipeline adapted for ALE-seq data analyses is used for the direct and reference-free annotation of new, active retroelements. This pipeline allows assessment of LTR retrotransposon activities in organisms for which genomic sequences and/or reference genomes are either unavailable or of low quality.