Phytochrome B Negatively Affects Cold Tolerance by Regulating OsDREB1 Gene Expression through Phytochrome Interacting Factor-Like Protein OsPIL16 in Rice

Cross talk between light signaling and cold signaling has been elucidated in the model plant Arabidopsis and tomato, but little is known about their relationship in rice. Here, we report that phytochrome B (phyB) mutants exhibit improved cold tolerance compared with wild type (WT) rice (Oryza sativa L. cv. Nipponbare). The phyB mutants had a lower electrolyte leakage index and malondialdehyde concentration than the WT, suggesting that they had greater cell membrane integrity and less lipid peroxidation. Real-time PCR analysis revealed that the expression levels of dehydration-responsive element binding protein 1 (OsDREB1) family genes, which functions in the cold stress response in rice, were increased in the phyB mutant under normal and cold stress conditions. PIFs are central players in phytochrome-mediated light signaling networks. To explore the relationship between rice PIFs and OsDREB1 gene expression, we produced overexpression lines of rice PIF genes. OsDREB1 family genes were up-regulated in OsPIL16-overexpression lines, which had improved cold tolerance relative to the WT. Chromatin immunoprecipitation (ChIP)-qPCR assay revealed that OsPIL16 can bind to the N-box region of OsDREB1B promoter. Expression pattern analyses revealed that OsPIL16 transcripts were induced by cold stress and was significantly higher in the phyB mutant than in the WT. Moreover, yeast two-hybrid assay showed that OsPIL16 can bind to rice PHYB. Based on these results, we propose that phyB deficiency positively regulates OsDREB1 expression through OsPIL16 to enhance cell membrane integrity and to reduce the malondialdehyde concentration, resulting in the improved cold tolerance of the phyB mutants.


INTRODUCTION
Light is a major environmental signal influencing a multitude of steps in plant development such as seed germination, carbon assimilation, stem elongation, leaf morphology, and flowering (Carvalho et al., 2011). Plants use phytochromes to sense red light-depleted (shade) and red light-enriched (full sun) conditions (Rockwell and Lagarias, 2006). Phytochromes are light-absorbing photoreceptors that exist in two fundamental forms: the red (R)-light-absorbing form, designated Pr, and the FR-light-absorbing form, designated Pfr. The inactive Pr form is converted to the active Pfr form by exposure to R-light, and is converted back to the inactive Pr form by exposure to FR-light or through dark reversion. Phytochromes interact with members of the basic bHLH family of PIFs to regulate the expression of a large number of light-responsive genes and thus influence many photomorphogenic events (Franklin and Quail, 2010;Nagatani, 2010;Gu et al., 2011). There are five phytochrome genes (PHYA to PHYE) and eight PIF genes (PIF1, PIF3, PIF4, PIF5, PIF6, PIF7, PIF8, and PIL1) in Arabidopsis thaliana, and three phytochrome genes (PHYA to PHYC) and six PIL genes (PIL11 to PIL16) in Oryza sativa (Nakamura et al., 2007;Leivar and Quail, 2011;Jeong and Choi, 2013;Luo et al., 2014).
In addition to their roles in plant photomorphogenesis, cross-talk between phytochrome-mediated light signals and cold signaling pathways has been identified in the model plant Arabidopsis and tomato (Williams et al., 1972;Kim et al., 2002;Wang et al., 2016). As early as 1972, Williams et al. noticed that phytochrome could mediate short day enhancement of cold acclimation, but the molecular basis for phytochrome regulation of cold tolerance remained unknown until recent years (Williams et al., 1972). CBF/DREB1 genes, which are rapidly and transiently induced by low temperatures, play an important role in cold stress responses (Gilmour et al., 2004;Heidarvand and Amiri, 2010;Mao and Chen, 2012). CBF/DREB1s can up-regulate COR genes that contain a C-repeat/drought-responsive promoter element, and thus cause metabolic changes to enhance cold tolerance (Thomashow, 2001). A low R/FR light signal increases CBF gene expression in Arabidopsis in a circadian clock-dependent manner. The transcript abundance of COR15a in a phytochromedeficient mutant analysis revealed that phyB and phyD repress the CBF regulon in high R/FR light in a non-redundant manner, and the phyD mutant showed enhanced cold tolerance (Franklin and Whitelam, 2007). Kidokoro et al. (2009) found that PIF7 specifically binds to the G-box (CACGTG) of the DREB1C and DREB1B promoters in Arabidopsis. Transactivation analysis using mesophyll protoplasts revealed that PIF7 functions as a transcriptional repressor of DREB1B and DREB1C expression under circadian control, and its activity is regulated by the PIF7interacting factors TOC1 and phyB (Kidokoro et al., 2009). Genetic analysis indicated that the CBF pathway is repressed by phyB, PIF4, and PIF7 under a warm long-day (LD) growing season, and this repression is relieved by short-day conditions. As a result, the freezing tolerance of short-day plants was increased relative to that of LD plants in Arabidopsis (Lee and Thomashow, 2012). Wang et al. (2016) found that FR and R light perceived by phyA and phyB positively and negatively regulated cold tolerance, respectively, in tomato (Solanum lycopersicum). FR light-induced activation of phyA triggers ABA signaling and, subsequently, JA signaling, leading to activation of the CBF pathway and a cold response in tomato plants (Wang et al., 2016).
Cold stress is an important factor limiting rice yield in many areas of high latitude and altitude. There are 10 putative DREB1 homologs (DREB1A to DREB1I) in rice, six of which (OsDREB1A, OsDREB1B, OsDREB1C, OsDREB1E, OsDREB1F, and OsDREB1G) were similarly expressed in response to chilling acclimation and cold stress. Moreover, these genes were coexpressed with genes involved in cold signaling, suggesting that they might be involved in the cold response in rice (Mao and Chen, 2012). Functional analysis revealed that over-expression of OsDREB1A, OsDREB1B, or OsDREB1F in rice and overexpression of OsDREB1D or OsDREB1F in Arabidopsis conferred enhanced cold tolerance in transgenic plants (Dubouzet et al., 2003;Wang et al., 2008;Zhang et al., 2009). PhyB deficiency alleviated chilling-induced photoinhibition in rice, probably through a more stabilized chloroplast structure and higher unsaturated fatty acid content in membrane lipids in the phyB mutant (Yang et al., 2013). Cordeiro et al. (2016) observed that OsPIF14 could bind to the OsDREB1B promoter through two N-boxes (CACG(A/C)G) and that the flanking regions of the hexameric core were essential for protein-DNA interaction and stability. Transactivation assays using Arabidopsis protoplasts and rice protoplasts showed that OsPIF14 downregulates OsDREB1B gene expression (Cordeiro et al., 2016). Furthermore, yeast two-hybrid and Co-immunoprecipitation analyses revealed that OsPIF14 preferentially binds to the active Pfr form of rice phyB (Cordeiro et al., 2016).
In this work, we observed that phyB mutants exhibited enhanced cold tolerance compared with the WT, suggesting that phyB may be involved in the regulation of tolerance to cold stress in rice. To investigate the mechanism by which phyB regulates cold tolerance, we analyzed the ELI and MDA content in WT and phyB mutant plants. Rice plants deficient in phyB exhibited reduced ELI and MDA values, presumably as a result of increased membrane integrity. We also analyzed the expression of DREB1 genes related to cold responses in the WT and phyB mutant. Genes acting upstream of the OsDREB1 gene family were further dissected.

Plant Materials and Stress Treatment
OsPIL16 overexpression, phyB1 mutant, phyB2 mutant and WT (O. sativa L., cv. Nipponbare) rice plants were used in this study. The genetic background of the phyB1 and phyB2 mutants was O. sativa cv. Nipponbare, which was described previously (Takano et al., 2001(Takano et al., , 2005. Seeds were surface sterilized in 75% (v/v) ethanol for 30 s and then in 5% NaClO (v/v) for 20 min. The seeds were then rinsed six times in sterile doubledistilled water. The sterilized seeds were soaked in water at room temperature for 2 days and then germinated for 1 day at 28 • C. Seedlings were grown hydroponically in Yoshida's culture solution as described previously (Zhou et al., 2013). Plants were cultured in a growth chamber at 28 • C/28 • C (day/night) under a 15-h-light/9-h-dark photoperiod. For cold stress, three-leaf stage seedlings were placed in a climatic chamber at 4 • C. All treatments were repeated three times. For R light treatment, seeds from WT and phyB were de-husked, surface-sterilized, and grown in 0.4% agar. After incubating at 4 • C overnight, the seeds were transferred to darkness for 7 days at 28 • C, and then transferred to R light for another 24 h. The above-ground parts were harvested for RNA isolation.

Cell Membrane Permeability
The cell membrane permeability of plantlets was assessed by the ELI of their tissues according to Sutinen et al. (1992) with the following modifications. The second leaves were collected from three individual plants at the three-leaf stage grown under normal conditions or treated at 4 • C for different time periods. After rinsing them three times in deionized water, the samples were cut into 1 cm long sections and placed in deionized water. After the leaf sections were placed in a vacuum for 10 min and kept at room temperature for 20 min, the electrical conductivity value (named C1) was determined. Then, the samples were boiled for 30 min, and the total conductivity value (named C2) was determined after the solution was cooled to room temperature. The ELI was expressed as ELI (%) = (C1/C2) × 100.

MDA Concentration Analysis
The MDA concentration was measured as described previously with the following modifications (Kim and Tai, 2011). The leaves were weighed and homogenized in 5 ml of 10% trichloroacetic acid solution. The homogenate was centrifuged, and 2 ml of the supernatant was added to 2 ml 0.67% thiobarbituric acid. The mixture was incubated in boiling water for 30 min, and the reaction was stopped in an ice bath. The MDA concentration was expressed as mol g −1 fresh mass (FM).

Quantitative PCR
Surface-sterilized seeds of the WT and phyB1 mutant were incubated under controlled photoperiodic conditions (15 h light, 28 • C/9 h dark, 28 • C) until the three-leaf stage. Samples from unstressed plants were collected as controls (0 h). Samples were taken throughout the cold stress treatment at 1, 4, 12, and 24 h after exposing the seedlings to cold stress. The collected samples were immediately flash frozen in liquid nitrogen and stored at −80 • C for further analysis. Total RNA was isolated from the third leaves using the RNAiso reagent (TaKaRa, Dalian, China). First-strand cDNAs were synthesized from total RNA using a PrimeScript TM RT reagent Kit with gDNA Eraser (Perfect Real Time), DNase treatment was applied before cDNAs synthesis to reduce DNA contamination according to the manufacturer's instructions. Quantitative PCR was performed on the Thermal Cycler Dice TM Real Time System (TaKaRa) using SYBR Premix Ex Taq TM (TaKaRa). Each reaction contained 10 µl of SYBR Premix Ex Taq TM (TaKaRa), 2 µl of cDNA sample, and 0.2 µl 20 µM gene-specific primer pairs in a final volume of 20 µl. The PCR thermal cycle used was as follows: denaturation at 95 • C for 30 s and 40 cycles of 95 • C for 5 s and 60 • C for 30 s. As an internal control, the rice elongation factor gene (OsEF1α, AK061464) was used to quantify the relative transcript level of each target gene. Three biological replicates were performed. The specific primers used in this study are listed in Supplementary Table 1.

Construction of the Plant Expression Vector and Rice Transformation
To construct the OsPIL16-overexpression (OsPIL16-OX) vector, the ORF region of OsPIL16 was amplified by PCR from cDNA using the primers OsPIL16-F and OsPIL16-R, and cloned into the p1390-Ubi vector between the maize ubiquitin promoter and the nos terminator (Li et al., 2012). To construct the OsPIL16-HA overexpression vector, the ORF region of OsPIL16 was amplified by PCR from cDNA using the primers OsPIL16-F and OsPIL16-HA-R, and cloned in frame at its 3 end with HA tag in the modified p1390-Ubi vector between the maize ubiquitin promoter and the nos terminator. The plasmid was introduced into Agrobacterium tumefaciens strain EHA105 by electroporation. Rice (O. sativa cv. Nipponbare) was transformed via the agroinfection method as described previously (Hood et al., 1993;Hiei et al., 1994).

DNA Extraction and Southern Blot Analysis
Genomic DNA was isolated and purified from the OsPIL16-OX lines and WT plants at the six-leaf stage following Murray and Thompson (Murray and Thompson, 1980). Approximately 100 µg of DNA was digested with HindIII and subjected to electrophoresis on a 0.8% agarose gel. The DNA was then transferred to a nylon membrane (Hybond-N + ; Amersham, Buckinghamshire, UK). PCR primers for amplification of the HptII gene (the selectable marker gene) are listed in Supplementary Table 1. The amplification conditions were one cycle of 3 min at 94 • C, 35 cycles of 30 s at 94 • C, 1 min at 56 • C and 1 min at 72 • C, and finally one cycle of 5 min at 72 • C. The amplified fragment was purified using a PCR clean-up system. DIG-labeled probe preparation and southern blotting were performed according to the DIG High Prime DNA Labeling and Detection Starter Kit I (Roche, Mannheim, Germany) instructions.

Yeast Two-Hybrid Assay
The full coding sequence of OsPIF16 was amplified by PCR from cDNA using the primers OsPIL16-AD-F and OsPIL16-AD-R, and cloned into the vector pGADT7. The C-terminal nonphotoactive region of PHYB (1891-3513 bp) was amplified by PCR from cDNA using the primers C-phyB-F and C-phyB-R, and cloned into pGBKT7. The bait and prey were cotransformed into yeast strain AH109 according to the Yeast Protocols Handbook. Then, positive yeast colonies were plated on SD -Leu -Trp -His media containing 20 mM 3-AT. pGADT7 was used as a prey for the negative control.

Cis-Acting Element Analysis
The promoter sequences of OsDREB1A, OsDREB1B, OsDREB1C, OsDREB1E, OsDREB1F, and OsDREB1G were downloaded from the Rice Genome Annotation Project 1 and used for cis-acting element analysis.

Western Blot
Seedlings of WT and OsPIL16-HA overexpression lines were ground in liquid nitrogen and homogenized in a denaturing buffer (100 mM NaP 2 PO4, 10 mM Tris-Cl, and 8 M urea) by vigorous vortexing. Cell debris was removed by centrifugation at 14, 000 g for 10 min at 4 • C. For western blot analysis, the supernatants were separated on 7.5% SDS-polyacrylamide gel, and the separated proteins were transferred to Immobilon-P-Transfer Membrane (Millipore, Bedford, USA). For detection of OsPIL16-HA, the membrane was incubated with rabbit polyclonal anti-HA antibody (Abcam, Cambridge, UK) in PBS buffer containing 0.05% Tween-20. Bands were visualized with BCIP/NBT Alkaline Phosphatase Color Development Kit (Beyotime, Nanjing, China), according to the manufacturer's instructions.

ChIP -qPCR Assay
ChIP assay was performed using an EpiQuik Plant Chip Kit (Epigentek, San Jose, California, USA). 14-day-old seedlings were incubated in the dark for 12 h after light irradiation and then crosslinked in 1% formaldehyde by vacuum filtration in the dark. ChIP analysis was performed by using an affinity purified anti-HA polyclonal antibody (ab9110, Abcam, Cambridge, UK) and normal IgG was used as negative control. Primers OsDREB1B-G-F and OsDREB1B-G-R, OsDREB1B-N-F and OsDREB1B-N-R that anneal to the G-box and N-box motifs in the OsDREB1B promoters, respectively, were used for qPCR. qPCR was performed using SYBR Green reagent (Takara, Dalian, Japan). Results were presented as percent of input.
All primers used in this study are listed in Supplementary  Table 1.

phyB Mutants Exhibit Enhanced Cold Tolerance
To test their cold tolerance, phyB mutant and WT rice plants were grown to the three-leaf stage in the same tray under well-watered conditions; the phyB1 and phyB2 mutants grew and developed normally (Figure 1A, Control). After 4 • C treatment for 4 days, both the WT and phyB mutant plants showed significant wilting symptoms (Figure 1A, Cold). The cold-stressed seedlings were then removed to normal growth conditions for 7 days. The phyB1 and phyB2 plants recovered almost completely, whereas the WT plants did not (Figure 1A, Recovered); 95 and 94% of the phyB1 and phyB2 plants, respectively, survived and grew new leaves, while 0% of the WT plants grew new leaves 7 days after being removed to normal conditions ( Figure 1B). These results unambiguously demonstrated that the phyB-deficient mutants had improved cold tolerance. To investigate the mechanism of cold tolerance, we monitored changes in the ELI, an indicator of membrane integrity. We first monitored the changes in the ELI over 24 h of cold stress in the WT and phyB1 mutant. Under normal growth conditions, the ELI values of the WT and phyB1 mutant were 4.5 and 2.8%, respectively ( Figure 1C). After cold stress treatment, ELI values increased in both the phyB1 mutant and the WT (Figure 1C). However, the phyB1 mutant exhibited an obviously lower ELI at 6 and 24 h after cold treatment compared with the WT (Figure 1C). We speculated that phyB deficiency probably increased the ability to maintain membrane integrity under cold stress, which is likely to be one of the factors enabling the phyB mutants to cope with cold stress.
Malondialdehyde is the end product of lipid peroxidation that is generated from the oxidation of polyunsaturated fatty acids by ROS and other free radicals in membranes, and acts as an indicator of lipid peroxidation (Liu et al., 2013). We measured the MDA concentration in the WT and the phyB1 mutant after cold stress treatment. The MDA concentration increased to about twofold in the phyB1 mutant and to about fourfold in the WT after 24 h of cold stress (Figure 1D), consistent with the changes in ELI. The lower MDA concentration in the phyB mutant indicated reduced membrane lipid peroxidation relative to the WT, which probably contributed to the enhanced membrane integrity under cold stress.

phyB Altered OsDREB1 Expression Patterns
Dehydration-responsive element binding protein 1 play an important role in cold responses, as indicated by the observation that overexpression of OsDREB1s can enhance cold tolerance in rice (Mao and Chen, 2012). To explore whether OsDREB1s were involved in the enhanced cold tolerance of the phyB mutant, we analyzed the expression patterns of OsDEREB1s in WT and phyB mutant plants grown under normal conditions and cold stress. Six of the nine OsDREB1 family members, including OsDREB1A, OsDREB1B, OsDREB1C, OsDREB1E, OsDREB1F, and OsDREB1G, were rapidly induced by cold stress in both the WT and phyB mutant, but their expression levels were relatively higher in the phyB mutant compared with the WT under normal conditions and cold stress (Figure 2). Although OsDREB1D was induced by cold treatment in the WT but not the phyB mutant, its expression was about twofold higher in the phyB1 mutant relative to the WT under normal conditions (Figure 2). OsDREB1J and OsDREB1I expression was not obviously induced by cold stress, but their transcription levels were clearly higher in the phyB mutant than in the WT. These results demonstrated that phyB deficiency promoted OsDREB1 family gene expression. Considering that overexpression of OsDREB1A, OsDREB1B, and OsDREB1F enhanced the cold tolerance of rice and that overexpression of OsDREB1F or OsDREB1D in Arabidopsis conferred enhanced cold tolerance (Ito et al., 2006;Wang et al., 2008), we speculated that increased expression of OsDREB1 family genes probably contributed to the enhanced cold tolerance of the phyB mutants.

phyB Altered the Expression Pattern of OsPILs
How does phyB regulate OsDREB1 expression? Because PIFs play a pivotal role in phytochrome-mediated light signaling networks (Castillon et al., 2007), we wondered whether OsPIFs were involved in the regulation of the cold response. Nakamura et al. (2007) identified six PIL homologs in rice, designated OsPIL11 through OsPIL16, by extensively evaluating all of the rice databases (Nakamura et al., 2007). To dissect the role of OsPILs in cold responses, we compared their expression patterns in WT and phyB mutant plants grown under normal conditions and cold stress. OsPIL12 expression levels were not affected by cold stress in the WT or phyB mutant, while OsPIL11, OsPIL13, and OsPIL14 were induced by cold stress in the phyB mutant but not the WT (Figure 3). OsPIL15 was induced by cold stress in both the WT and phyB mutant, whereas cold stress induced OsPIL16 expression only in the WT, not in the phyB mutant (Figure 3). The expression levels of OsPIF15 and OsPIF16 were much higher in the phyB mutant than in the WT under normal conditions and cold stress (Figure 3). In particular, OsPIL16 transcript levels were about 400-fold higher in the phyB mutant than in the WT. To understand how phyB regulates OsPIL16 gene, we compared the transcript levels between WT and phyB mutant seedlings that were grown in the dark, or grown in the dark but irradiated with 24 h of continuous R light before harvest. OsPIL16 transcript level was decreased by R light in WT but not in phyB (Supplementary Figure 1), suggesting that phyB perceive R light to inhibit OsPIL16 gene expression at the transcriptional level. Based on these results, we deduce that the expression of multiple OsPILs was regulated by cold stress, and that phyB was involved in modulating the OsPIL gene expression.

OsPIL16-Overexpression Lines Have Enhanced Cold Tolerance and Increased Expression Levels of OsDREB1s
To investigate the role of OsPILs in regulating cold tolerance in rice, we produced the overexpression lines of OsPIL15 and OsPIL16 because of the significant difference in their transcript levels between the WT and phyB (Figure 3). The OsPIL15-overexpression (OsPIL15-OX) lines were reported in our previous study (Zhou et al., 2014). In this study, we produced transgenic rice lines overexpressing OsPIL16 (Figure 4A). Two  independent homozygous lines of the T 4 progeny (#1 and #13) were used to analyze the role of OsPIL16 in rice based on their high expression levels (Figures 4B,C). Seedlings of the OsPIL16-OX lines and WT at the three-leaf stage were hydroponically grown under cold conditions for 4 days and then allowed to recover for 7 days. The OsPIL16-OX lines exhibited enhanced cold tolerance, as more seedlings survived compared with the WT (Figure 4D). However, OsPIL15-OX lines exhibited the same sensitivity as the WT (data not shown). These results suggested that OsPIL16 was involved in the cold response. In Arabidopsis, PIF4 and PIF7 down regulate the CBF pathway and freezing tolerance under LD conditions (Lee and Real-time PCR was performed using cDNA derived from leaves. Error bars indicate SE (n = 3). (D) Representative photographs of the WT and OsPIL16-OX lines under cold stress. Three-leaf stage seedlings were treated with cold (4 • C) for 4 days and then allowed to recover at 28 • C for an additional 7 days. (E) OsDREB1 expression levels in the OsPIL16-OX lines and WT. Real-time PCR was performed using cDNA derived from the leaves of three-leaf stage seedlings. Error bars indicate SE (n = 3). Thomashow, 2012). To clarify the relationship between the OsPILs and OsDREB1s, we compared the transcript levels of the OsDREB1s between the WT and OsPIL16-OX or OsPIL15-OX lines. As shown in Figure 4E, the expression levels of OsDREB1A, OsDREB1B, OsDREB1C, OsDREB1E, OsDREB1F, and OsDREB1G were increased in the OsPIL16-OX lines compared with the WT. These results suggested that OsPIL16 enhanced cold tolerance by up regulating OsDREB1 expression. However, the transcript levels of the OsDREB1 genes were similar in the WT and OsPIL15-OX lines (Supplementary Figure 2). The expression patterns of the OsDREB1s were consistent with the cold tolerance of the OsPIL16-OX and OsPIL15-OX lines.
A previous study revealed that PIFs can bind to the G-box, PBE (CACATG) or N-box (CACGCG or CACGAG) (Cordeiro et al., 2016;Kim et al., 2016). Promoter analysis of OsDREB1A, OsDREB1B, OsDREB1C, OsDREB1D, OsDREB1E, OsDREB1F, and OsDREB1G revealed that they all contained a G-box, PBE, or N-box (Table 1). This result suggested that OsPIL16 probably upregulates the expression of OsDREB1 genes by directly binding to their promoter regions to enhance cold tolerance. To test this hypothesis, we produced the OsPIL16-HA overexpression (OsPIL16-HA OX) lines and detected OsPIL16-HA protein using western blotting (Supplementary Figure 3). ChIP-qPCR assays using primers that anneal to the G-box and N-box motifs in the OsDREB1B promoters were performed. The results show that N-box regions were enriched in the chromatin fractions by anti-HA antibody compared to that by control IgG (Figure 5).

phyB Interacts with OsPIL16 In vivo
Recent evidence suggests that phytochrome signaling is initiated by the direct interaction of the biologically active forms of phytochromes with PIF members in Arabidopsis (Castillon et al., 2007). To examine the physical interaction between phyB and OsPIL16, we performed a yeast two-hybrid assay, using the coding region of OsPIF16 as prey and the C-terminal nonphotoactive coding region of phyB as bait. Three independent clones were used for each sample. As shown in Figure 6 and Supplementary Figure 4, OsPIL16 interacted with the OsPHYB C-terminal domain. Thus, we speculated that phyB affects cold responses through an OsPIL16-mediated pathway.

DISCUSSION
phyB plays an important role in rice photomorphogenesis (Takano et al., 2005;Jumtee et al., 2009). In this study, we demonstrated that phyB deficiency resulted in improved cold tolerance compared with the WT though the regulation of OsDREB1 gene expression via an OsPIL16-mediated pathway (Figure 7). The N-box, PBE (PIF binding E-box), and G-box are marked by capital letters.
FIGURE 5 | Chromatin immunoprecipitation assays of 14-day-old seedlings grown under long-day conditions and shifted to dark for 12 h. Relative enrichment in the G-box region and N-box region of OsDREB1B promoter was quantified by qPCR.

The Enhanced Cold Tolerance of phyB Mutants Is Mainly Attributable to Increased Membrane Integrity and Reduced MDA Content
Cold stress in rice seedlings can negatively affect growth and development or even cause death . In this study, phyB deficiency enhanced cold tolerance compared with the WT. One of the major forms of cold damage is severe membrane injury. ELI is an indicator of membrane integrity (Heidarvand and Maali-Amiri, 2013). The ELI was clearly decreased in the phyB mutants compared with the WT (Figure 1C), indicating that phyB deficiency increased membrane integrity under cold stress. This result is consistent with our previous report that a phyB mutant had a more stable chloroplast structure under cold stress (Yang et al., 2013). MDA is the degradation product of polyunsaturated lipids, and the increased membrane integrity of the phyB mutant can be partly attributed to reduced lipid peroxidation as shown by the reduced MDA content ( Figure 1D). However, Yang et al. (2013) revealed that phyB mutants had higher unsaturated fatty acid contents in their membrane lipids, which also partly contributes to the increased membrane integrity of the phyB mutant (Yang et al., 2013). Therefore, we speculate that phyB deficiency improves membrane integrity in rice, which is an important factor for enhanced cold tolerance in phyB mutants. Additionally, MDA is a reactive aldehyde that initiates toxic stress in cells and subsequently causes cellular dysfunction and tissue damage . Thus, we speculate that the reduced MDA content under cold stress caused by phyB deficiency is likely one of the factors contributing to cold tolerance in the phyB mutants.
phyB Regulates OsDREB1 Expression through OsPIFs to Affect Rice Cold Tolerance C repeat binding factor/DREB1s play an important role in cold signal transduction (Cook et al., 2004;Zhao et al., 2016).
Overexpression of rice or Arabidopsis DREB1s can enhance cold tolerance (Jagloottosen et al., 1998;Dubouzet et al., 2003;Zhang et al., 2009). Interestingly, we found that the phyB mutant had higher expression levels of OsDREB1s (Figure 2), which is probably related to the enhanced cold tolerance in the phyB mutant. Phytochromes have been reported to regulate DREB1 gene expression in Arabidopsis. One study identified phyB as a positive regulator of cold-inducible DRE-controlled expression (Kim et al., 2002). However, the results of a subsequent study indicated a repressive role for phyB and phyD in the control of the CBF-responsive regulon (Franklin and Whitelam, 2007). Kidokoro et al. (2009) also showed that phyB-mediated light signaling had a repressive effect on DREB1 gene expression in Arabidopsis, consistent with our results in this study (Kidokoro et al., 2009). Xu et al. (2011) found that overexpression of ZmCBF3 in rice reduced the MDA concentration (Xu et al., 2011).
In this context, we speculate that the reduced MDA concentration in the phyB mutant probably contributed to the upregulation of OsDREB1 genes in our study. How does phytochrome-mediated pathway regulate DREB1 family gene expression? PIFs, as central players in phytochromemediated light signaling networks, have been reported to regulate DREB1 gene expression in Arabidopsis. PIF7 can bind to the promoter of DREB1C, and the pif7 null mutant was shown to have altered transcript levels of DREB1B and DREB1C (Kidokoro et al., 2009). PIF4 can downregulate the CBF pathway and freezing tolerance under LD conditions (Lee and Thomashow, 2012). There are six OsPIFs (OsPIL11-OsPIL16) in rice (Nakamura et al., 2007). The functions of OsPIL13 and OsPIL15 in rice growth and development have been reported. Overexpression of OsPIL13 in transgenic rice plants promoted internode elongation (Todaka et al., 2012). OsPIL15 represses FIGURE 6 | Analysis of the interaction between OsPIF16 and OsPHYB using a yeast two-hybrid assay. (A) Protein constructs used in the yeast two-hybrid assay. OsPIF16 was used as a prey in a translational fusion with the GAL4 AD. The C-terminal non-photoactive region of rice PHYB was fused with the GAL4 BD to use as a bait. (B) Analysis of protein-protein interactions in yeast growing in medium with (left) or without (right) histidine. The negative control was the interaction between OsPIF16 and GBD. seedling growth in the dark, and when exposed to R or FR light, OsPIL15-OX lines can relieve growth retardation and promote seedling elongation (Zhou et al., 2014). However, to our knowledge, the roles of OsPIL13 and OsPIL15 in regulating OsDREB1 gene expression have not been reported. In this study, our results showed that OsPIL16 upregulated the expression of multiple OsDREB1 genes (Figure 4E), which is consistent with the enhanced cold tolerance of OsPIL16-OX transgenic lines ( Figure 4D). OsPIL14 was reported to downregulate OsDREB1B gene expression in rice protoplasts by binding to the N-box motif of the OsDREB1B promoter (Cordeiro et al., 2016). How does OsPIL16, as a bHLH transcription factor, regulate the expression of DREB1s? bHLH proteins are able to bind to hexameric E-box (CANNTG) or N-box motifs, depending on their transcriptional activity (Cordeiro et al., 2016). A previous study revealed that PIFs not only bind to the G-box in vitro and in vivo, and also bind weakly in vitro to PBE (Kim et al., 2008(Kim et al., , 2016Zhang et al., 2013;Pfeiffer et al., 2014). We found that there were N-boxes, PBE or G-boxes in the promoter regions of the OsDREB1 genes ( Table 1). ChIP-qPCR assay confirmed that OsPIL16 can bind to the N-box of OsDREB1B promoter (Supplementary Figure 3), thus to regulate the expression of OsDREB1B.
To dissect the relationship between phyB and OsPILs, we compared the transcript levels of OsPILs in the WT and phyB mutant. The transcript levels of both OsPIL15 and OsPIL16 were significantly higher in the phyB mutant than in the WT (Figure 2). However, OsPIL15-OX did not show higher OsDREB1 transcript levels (Supplementary Figure 2), which is probably consistent with the similar cold response of the OsPIL15-OX line to the WT. Unlike OsPIL15-OX, OsPIL16-OX had higher OsDREB1 transcript levels (Figure 4D), which contributed to the enhanced cold tolerance in the OsPIL16-OX lines ( Figure 4C). To explain how phyB affects OsDREB1s family gene expression through OsPIL16, it is necessary to clarify the relationship between phyB and OsPIL16. Therefore, we examined the regulation of phyB on OsPIL16 transcript levels and physical interaction between phyB and OsPIL16 proteins. As shown in the results (Figure 3; Supplementary Figure 1), phyB perceives the R light to inhibit OsPIL16 expression at the transcriptional level. Physical interaction between OsPIL16 and OsPHYB were also confirmed by Y2H (Figure 6; Supplementary Figure 4). There is a well-known model of PIF function in phytochrome signaling pathways in Arabidopsis (Castillon et al., 2007;Leivar and Quail, 2011). Light signals induce photoconversion of phytochromes to the active Pfr forms before nuclear migration. In the nucleus, phytochromes physically interact with PIFs, which results in phosphorylation of PIFs. The phosphorylated forms of PIFs are subsequently degraded by the 26S proteasome. The lightinduced proteolytic removal of PIFs results in relieving the negative regulation of photomorphogenesis. However, Park et al. (2012) reported that phyB inhibits the regulatory activity of PIF1 and PIF3 by two different modes of action: by releasing them from their target promoters and by mediating their degradation. In this study, OsPIL16 induced OsDREB1s gene expression probably through binding to the N-box of their promoters. Thus, OsPIL16 is the positive regulator of OsDREB1 gene expression and rice cold tolerance. However, the physical interaction between phyB and OsPIL16 proteins inhibits the OsDREB1 gene expression by either degradation of OsPIL16 protein or by inhibiting the bind of OsPIL16 to N-box of OsDREB1 gene promoter. This hypothesis can explain the higher transcript levels of OsDREB1s in the phyB mutants than in the WT.
In summary, we determined the function of phyB in the cold tolerance of rice using genetic and physiological approaches. Although our findings suggest that phyB affect rice cold tolerance through OsPIL16 pathway in this study, some other pathways, for example ABA pathway and ICE pathway, are probably also involved in the phyB-regulated cold responses in rice. ABA has been shown to increase levels of CBF1 transcript and protein (Knight et al., 2004). Our previous result revealed that light signals mediated by phytochrome B affect ABA pathway in rice (Gu et al., 2012). To test this possibility, we must produce the null mutants of pil16 and phyBpil16 double mutant in the future experiments. In addition, further work is needed to decipher the mechanisms by which OsPIL16 regulates OsDREB1s.

AUTHOR CONTRIBUTIONS
YH and XX conceived the idea, led the study design, data analysis, and manuscript writing. YL performed southern blot, assisted in quantitative PCR and yeast two-hybrid. CZ and LX performed vector construction and rice transformation. LC and GZ assisted in the measurement of ELI and MDA. LC and YH performed ChIP-qPCR. JZ assisted in rice growing.