Abstract
Cold acclimation (CA) and deacclimation (DA), which are often accompanied by changes in freezing tolerance (FT), carbohydrates and hormones, are crucial for winter survival, especially under global warming. Plants with weak CA and premature DA caused by warm winters and/or unseasonal warm spells can be easily injured by adverse reactions to cold. Thus, understanding the molecular mechanisms of FT is imperative. In this study, we used high-throughput RNA-seq to profile the CA and DA of leaves of overwintering Rhododendron "Miyo-no-Sakae" over time; these leaves do not undergo dormancy but do undergo photoprotection during CA, and they do not grow during DA. Using Mfuzz and weighted gene coexpression network analysis, we identified specific transcriptional characteristics in each phase of CA and DA and proposed networks involving coexpressed genes and physiological traits. In particular, we discovered that the circadian rhythm is critical for obtaining the strongest FT, and high expression of circadian rhythm-related genes might be linked to sugar accumulation during winter. Furthermore, evergreen leaves exhibited robust photoprotection during winter, as revealed by high values of nonphotochemical quenching, high expression of transcripts annotated as "early light-induced proteins", loss of granum stacks and destacking of thylakoids, all of which were alleviated during DA. The strong requirement of photoprotection could be the reason for decreased abscisic acid (ABA) and jasmonic acid (JA) contents during CA, and decreases in ABA and JA contents may contribute to decreases in lignin content. Our data suggest that the molecular mechanisms of FT in overwintering leaves are unique, which may be due to the high requirements for photoprotection during winter.