Figures (4)  Tables (1)
    • Figure 1. 

      Cold transcriptional responses of the three rice varieties. (a) Cold treatment phenotypes of the three rice varieties; the white scale bar represents 2 cm. (b) Principal component analysis (PCA) of the RNA-seq data. (c) Statistics of differentially expressed genes (DEGs). (d) Overlap analysis of DEGs at the 24 h time point. (e) Gene ontology (GO) enrichment analysis of genes with higher expression in KH and LJ compared to YD at the 24 h time point. (f) Selected cold-induced genes with higher expression in KH and LJ compared to YD.

    • Figure 2. 

      Chromatin accessibility and transcription factor binding motif enrichment analysis. (a) Distribution of DNase-seq data relative to gene structures. (b) Distribution of chromatin accessibility signals relative to gene structures. (c) TOBIAS TF footprinting analysis of ACRs in KH at the 24 h time point; the differential binding score in TOBIAS is a normalized statistic summarizing the overall footprint score change of a TF's motifs between two conditions relative to background. (d) TF binding motif enrichment analysis in the promoter ACRs of DEGs: the enrichment ratio is defined as the proportion of transcription factor footprints present in promoter-associated ACRs of a specific gene set divided by that in all promoter-associated ACRs.

    • Figure 3. 

      DAP-seq data analysis of 23 cold-induced transcription factors. (a) Expression profiles of the 23 cold-induced transcription factors. (b) Clustering of DAP-seq sample signals based on the overlap (Szymkiewicz–Simpson) coefficient. (c) Number of peaks in DAP-seq samples and their distribution relative to gene structures. (d) TF motif enrichment analysis based on DAP-seq signals; the most significant motif identified through de novo prediction is used as the representative motif for the TF.

    • Figure 4. 

      DAP-seq defines potential TF-target gene regulatory relationships. (a) Number of potential target genes for each TF defined by DAP-seq and DNase-seq. The top panel shows the number of target genes defined by promoters containing DAP-seq signals, while the bottom panel shows the number of target genes defined by promoters containing overlapping signals from both DAP-seq and DNase-seq. (b) In vitro binding of TFs to the promoters of stress-related genes; the different colored tracks represent distinct transcription factor families.

    • Locus IDTF familySymbolDAP-seq data source
      LOC_Os05g03860bZIPbZIP38This study
      LOC_Os05g34050bZIPbZIP39[27]
      LOC_Os03g55540C2H2ADR3This study
      LOC_Os07g48630EILEIL2This study
      LOC_Os02g45450ERFDREB1GThis study
      LOC_Os06g03670ERFDREB1C[28]
      LOC_Os07g02800G2-likeNIGT2[26]
      LOC_Os01g60270HomeoboxWOX13This study
      LOC_Os05g45410HSFSPL7This study
      LOC_Os11g10130MYBMYB1This study
      LOC_Os06g43090MYBMYB102This study
      LOC_Os10g33810MYBMYB110This study
      LOC_Os11g45740MYBMYB12This study
      LOC_Os03g20090MYBMYB2This study
      LOC_Os02g41510MYBMYB30This study
      LOC_Os04g43680MYBMYB4This study
      LOC_Os04g49450MYB-relatedMYB511This study
      LOC_Os03g04070NACNAC022This study
      LOC_Os03g60080NACNAC19This study
      LOC_Os01g60020NACNAC4This study
      LOC_Os11g08210NACNAC5This study
      LOC_Os05g27730WRKYWRKY53This study
      LOC_Os05g25770WRKYWRKY45This study

      Table 1. 

      Transcription factor genes used for DAP-seq in this study.