Understanding the early cold response mechanism in IR64 indica rice variety through comparative transcriptome analysis
BMC Genomics volume 21, Article number: 425 (2020)
Cellular reprogramming in response to environmental stress involves alteration of gene expression, changes in the protein and metabolite profile for ensuring better stress management in plants. Similar to other plant species originating in tropical and sub-tropical areas, indica rice is highly sensitive to low temperature that adversely affects its growth and grain productivity. Substantial work has been done to understand cold induced changes in gene expression in rice plants. However, adequate information is not available for early gene expression, especially in indica variety. Therefore, a transcriptome profile was generated for cold shock treated seedlings of IR64 variety to identify early responsive genes.
The functional annotation of early DEGs shows enrichment of genes involved in altered membrane rigidity and electrolytic leakage, the onset of calcium signaling, ROS generation and activation of stress responsive transcription factors in IR64. Gene regulatory network suggests that cold shock induced Ca2+ signaling activates DREB/CBF pathway and other groups of transcription factors such as MYB, NAC and ZFP; for activating various cold-responsive genes. The analysis also indicates that cold induced signaling proteins like RLKs, RLCKs, CDPKs and MAPKK and ROS signaling proteins. Further, several late-embryogenesis-abundant (LEA), dehydrins and low temperature-induced-genes were upregulated under early cold shock condition, indicating the onset of water-deficit conditions. Expression profiling in different high yielding cultivars shows high expression of cold-responsive genes in Heera and CB1 indica varieties. These varieties show low levels of cold induced ROS production, electrolytic leakage and high germination rate post-cold stress, compared to IR36 and IR64. Collectively, these results suggest that these varieties may have improved adaptability to cold stress.
The results of this study provide insights about early responsive events in Oryza sativa l.ssp. indica cv IR64 in response to cold stress. Our data shows the onset of cold response is associated with upregulation of stress responsive TFs, hydrophilic proteins and signaling molecules, whereas, the genes coding for cellular biosynthetic enzymes, cell cycle control and growth-related TFs are downregulated. This study reports that the generation of ROS is integral to the early response to trigger the ROS mediated signaling events during later stages.
Oryza sativa L. ssp. Indica being a tropical crop is highly sensitive to low-temperature stress leading to impaired growth and massive losses in grain productivity. Reports suggest that rice plants are more susceptible to cold stress during the seedling, tillering, panicle development and flowering stages. Rice grown below the ambient temperature leads to a lower rate of germination, retarded seedling emergence, delayed vegetative growth, reduced rates of photosynthesis. Continued exposure to cold stress causes tissue necrosis and ultimately, cellular death [1,2,3,4]. Thus, the agronomic productivity of rice plants is heavily affected under low-temperature conditions, especially in the elevated regions where cold mountain water is used for irrigation.
After perceiving the cold stress, cells undergo increase in their membrane rigidity due to a reduction in the plasma membrane fluidity . This rise in the membrane rigidity causes increased electrolytic leakage from the cell, which acts as the primary signal for triggering cold response via activation of the cold-responsive gene expression [6,7,8,9]. Studies have shown that increase in membrane rigidity activates early cytoplasmic signals, such as triggering the MAPK signaling cascade, and the influx of cytoplasmic Ca2+ via mechano-sensitive Ca2+ channel or ligand-gated Ca2+ channel . The increase of cytoplasmic Ca2+ activates a myriad of downstream signaling pathways, mainly via calcium decoders such as calcium-dependent protein kinases (CDPKs), to further activate the transcription of cold-responsive transcription factor belonging to C-repeat binding factor (CBF)/ Dehydration responsive element binding (DREB) family . The DREB transcription factor binds to DRE site in the promoter region, thereby activating the expression of many cold-responsive (COR) genes (such as LT1, KIN, RAB, ERD genes). These DREB regulon genes play an important role in stabilizing membrane structure, activating ROS scavengers, and promoting the production of osmoprotectants to protect both the membrane and organelle damage during cold stress [12,13,14,15,16]. Microarray analysis has identified other transcription factors such as HSF1C, ZAT12, ZF, ZAT10 and SZF2 that are co-expressed with CBF and can positively regulate COR gene expression to impart cold stress tolerance . Other transcription factors have been reported in Arabidopsis such as Eskimo1 and HOS9, which also participate in freezing tolerance and together constitute the CBF independent regulation [18, 19].
High-throughput RNA sequencing data have contributed significantly to understanding the molecular mechanism of cold response in rice. Owing to the diverse growing conditions and availability of various rice cultivable varieties worldwide; it is integral to continue with the high-throughput study of different varieties. Such studies provide a better understanding of the complexity of cold signaling that greatly enhances growth and grain productivity. Moreover, the response to cold stress in rice varies with tissue type, as well as varying developmental stages. Further, more pronounced effects are observed during the seedling stage and flower development in rice plants, when subjected to cold stress conditions. Among the two major subspecies of rice, the japonica varieties, usually grown at higher altitudes, are more tolerant to cold stress compared to indica varieties which are typically grown in tropical regions [20,21,22]. Previous studies have shown that cold-responsive genes can be clustered under two major groups, i.e., regulatory protein-coding genes that perceive the signal and functional protein-coding genes which initiate the abiotic stress response [23, 24]. The genes coding for regulatory proteins includes signaling molecules, such as kinases, phosphatases, calcium-binding messenger molecules, transcription factors, micro RNAs, and Two-Component Systems [25,26,27]. The second cluster of functional protein-coding genes comprise antioxidants, players of ROS removal, compatible solutes, and other hydrophilic molecules crucial for maintaining the osmotic balance [27,28,29,30]. However, the expression of this group of genes majorly depends on the stress exposure time and tissue type.
Previous work has identified several cold-tolerant wild rice varieties, such as Dongxiang common wild rice, Chaling wild rice, and Guangxi wild rice [31,32,33]. These varieties have been reported to withstand temperatures as low as − 9 °C to − 12 °C. Shen et al., 2014, suggested that Dongxiang common wild rice is an ideal germplasm source for the generation of cold-resistance breeding [34, 35]. In the modern era, various approaches like hybridizations among wild rice varieties and novel genetic manipulation methods have been used to improve the crop’s growth and productivity under cold conditions. Cold tolerant introgressed lines, generated from indica variety, with japonica rice have shown increased tolerance to cold stress compared to cold-sensitive indica varieties . Although strong selection pressure may be attributed to the evolution of such cold-tolerant genotypes, the molecular basis of tolerance in such intrinsically tolerant rice varieties is still a less explored field. In the present scenario, the knowledge of cold-responsive genes, and the cold-tolerant QTLs triggered at various stages, during cold stress can be exploited for breeding cold-tolerant rice varieties. Indica rice varieties with better adaptability to low-temperature conditions need to be identified. Furthermore, improving high yielding tropical indica varieties, such that they are better suited for lower temperatures, can provide a solution to the loss of yield in the paddy fields of high-altitude terrain. This study was aimed at identifying genes that undergo differential expression in response to early cold stress in the IR64 rice variety and further extended to different indica cultivars to understand the early signaling events associated with the cold stress response.
Cold shock induces differential gene expression in IR64 indica rice
This work aimed at studying the early cold stress response in rice (Oryza sativa L. ssp. indica) involved the treatment of 14 days old rice seedlings to cold shock (at 4 °C) for 2 h. cDNA libraries generated from IR64 seedlings grown under the control (28 °C ± 1 °C) and 2 h cold shock (4 °C) conditions are denoted as CT replicates and CS replicates respectively. A total of 58.53 million (CT1), 33.61 million (CT2), 60.85 million (CS1), 39.82 million reads (CS2) were generated for the IR64 rice cultivars. Statistics of cleaned reads (approx. 63 nucleotides) were assessed with FastQC, which revealed that all reads were of fairly good quality and without adapters. The read mapping was carried out using HiSat2, resulting in 95.5% reads mapping to reference Japonica genome Oryza sativa japonica (Os-Nipponbare-Reference-IRGSP-1.0 (IRGSP-1.0) (Table 1). Normalized expression profiling was done on the aligned reads resulting in the identification of 32,161 transcripts expressed in at least one of the four samples profiled. As shown in Fig. 1a, 24,988 transcripts were found in both the control data sets, whereas 24,836 transcripts were common for both cold shock replicates. When compared between control and cold data sets, 72% (23232) of the transcripts were expressed in both, indicating a basal level of expression. Further, 539 transcripts were found exclusively expressed in control, whereas 931 transcripts were expressed only under cold shock condition.
We used Deseq2 package for differential expression analysis of the genes. A filter with a p-value cut-off of < 0.05 and log2fold change ≥1.5 and ≤ − 1.5 was set as the criteria to identify the differentially expressed genes (DEGs). These DEGs were visualized using volcano plot (Fig. 1b) to understand the distribution of up and downregulated genes. For this analysis, FPKM of > = 0.1 for a transcript was considered as expressed supported by a median read count of at least five reads per transcript covering 100% of the sequence. Among these DEGs, 380 genes were upregulated, whereas, 136 genes were downregulated in cold shock (CS) Vs control condition (CT) seedlings. (Additional file 1) Analysis of unsupervised hierarchical clustering of differentially expressed transcripts shows distinct gene expression patterns of up and down-regulation levels during cold shock treated (CS) (Fig. 1c). The differentially expressed genes could be categorized into five different clusters, based on their expression patterns. The largest group, Cluster V comprises of the differentially expressed genes, upregulated under cold shock, whereas, Clusters II and IV constitute the downregulated genes under cold shock. The clustering of all these replicates exhibited high sample reproducibility.
Functional annotation of cold shock-induced genes indicate a significant increase in cold-responsive TF and ROS activity
To understand the biological function of the cold induced differentially expressed genes (DEGs), GO enrichment was performed using an FDR adjusted p-value of ≤0.05 as the cut-off. The Blast2GO analysis for 516 DEGs featured 234 GO term annotation in biological process, 273 in molecular function (MF), and 262 for cellular component. Comparative analysis of the upregulated and downregulated GO terms indicates cell division (GO:0051301), proliferation (GO:0008283), developmental processes (GO:0032502) and growth (GO:0040007) were specific for downregulated genes. GO terms such as transport (GO:0006810), and homeostatic process (GO:0042592) were majorly associated with the upregulated genes (Fig. 2a).
GO analysis identified oxidation-reduction process, processes related to water stress and lignin metabolism were significantly enriched during cold shock treatment, in addition to generic terms such as cellular biosynthetic processes and transcription regulation (Fig. 2b). Under stress response, significant enrichment for GO-terms such as response to alcohol (GO:0097305), response to temperature stimulus (GO:0009266), response to abscisic acid (GO:0009737), response to lipid (GO:0033993), response to acid chemical (GO:0001101), response to osmotic stress (GO:0006970), and cold acclimation (GO:0009631) was observed (Fig. 2c). For the oxidation-reduction process, response to Oxygen-containing compound (GO: 1901700) and lignin catabolic process (GO:0046274) were significantly enriched under cold shock. GO-molecular function (MF) terms comparison indicates that metal ion binding (GO:0005488), oxidoreductase activity (GO:0016491) and transcription regulator activity (GO:0140100) were enriched for upregulated genes (Fig. 2c). GO enrichment analysis was also performed for the downregulated genes, but no significantly enriched terms were detected.
The biological pathway associated (KEGG pathway) with cold shock response was analyzed using BLASTKOALA (24.9% of input sequences). The analysis revealed that most genes were assigned to metabolism (40) of carbohydrate, amino acid, lipid and secondary metabolites; environmental information processing (15), genetic information processing (5) like transcription, translation and protein processes (Additional file 6C). Further, KEGG-BRITE reconstruction revealed that compared to control, a higher number of genes were assigned to ko01000 Enzymes (42), ko02000 Transporters (8), ko01003 Glycosyltransferases (6), ko03000 Transcription factors (4) and, ko04147 exosomes (4) in the cold shock treated sample (Additional file 2). The Interpro domain search data indicated that DNA binding domain and cytochromes were most abundant in the upregulated genes (Additional file 6B).
Gene regulatory network induced during early cold stress
Analysis of the DEGs showed the presence of a milieu of stress responsive genes upregulated, which included heat shock protein genes (Os02g0758000, Os03g0266900, Os06g0253100) Terpene synthases (OsTPS1, OsTPS31), Dehydrins, LEA and RAB group of proteins coding genes (OsDHN1, OsLEA28, OsRAB16, OsLEA14), pathogen-related proteins and chitinase and glucanases. Among signaling molecules, calcium-calmodulin molecules and receptor kinases and protein phosphatases, along with several redox homeostasis proteins were induced under cold shock conditions. Besides hydrophilic proteins and signaling protein, our data set indicate the presence of 38 upregulated TFs (10% of total upregulated DEGs) and 9 downregulated TFs in cold shock transcriptome (Fig. 3a and b; Additional file 3). Gene network of these upregulated DEGs shows three major clusters that are highly interconnected. Cluster I comprise Zinc finger and NAC transcription factors and signaling proteins such as calmodulin and kinases. Cluster II represents DREB/AP2 and MYB transcription factors as major nodes. Cluster III contain proteins that mostly belong to osmoprotectants activated in response to dehydration stress (Fig. 3c).
Cell wall modification and ROS generation are crucial to stress perception during early cold shock in IR64
Differentially expressed gene set unique to this study (Additional file 4), has a significant number of genes responsible for cell wall modification and ROS generation (Fig. 3a and b). The genes coding for redox molecules comprised majorly of lignin catabolic laccase genes (OsLAC10: Os02g0749700, OsLAC17: Os10g0346300, OsLAC23: Os11g0641500, and OsLAC29: Os12g0258700), the germin-like oxalate oxidases, and other ROS generating enzymes. These findings suggest that the generation of ROS occurs during early cold shock and is essential for activating the redox signaling at the later stages of the stress response. Other stress responsive genes, unique to 2 h cold shock include the terpene biosynthesis genes (OsTPS1, OsTPS10), salt stress responsive lectin proteins (Os01g0348800, Os01g0348900) and receptor-like kinases (Os11g0672200, Os04g0540900). Other genes that were induced within 2 h include cell wall degrading enzymes, such as, chitinases (Os05g0399300, Os11g0701200, Os11g0702100), cellulases (Os01g0946600, Os01g0946700), pectin methylesterase (Os04g0458900), and a group of membrane transporter genes.
Cold shock-induced TFs upregulated under cold shock constitute major gene regulatory networks
Sequence analysis suggests that around 5–7% of coding sequences in plant genomes constitute transcription factors [37, 38]. In plants, the role of AP2/EREBP, bZIP, NAC/NAM (ATAF and CUC), MYC/MYB, and WRKY transcription factor families has been elucidated in the abiotic stress response that regulates stress-responsive gene expression via ABA-dependent or independent pathways [23, 39, 40]. The upregulated transcription factors in this study belong to various families; the majority being the AP2/ERF, DREB and MYB group of TFs (Fig. 4a). Downregulated genes included nine transcription factors needed for growth and development of the plant, which consisted primarily of bHLH TFs, and growth-related TFs such as OsGIF3 (Os03g0733600), OsGRAS1(Os01g0646300) and OsGRF7 (Os12g0484900) (Additional file 3). A gene regulatory network analysis using STRING shows that these upregulated transcription factors constitute a major network consisting of 15 nodes and a second network with 3 nodes. Further, the search suggested MYB2, MYB4, DREB1B, ZFP37, DREB1E and, DREB1G were highly connected and formed the central cluster (Cluster I, Fig. 4b). The cluster II consists of NAC39, which connects other TF like DERF5 (ERF103), NAC077 and ERF71. The heat shock transcription factors, HSF21, HSF11, are connected among themselves and also connected to central cluster via DREB1E, DREB1G and MYB2. HOX transcription factors constitute a secondary network which may contribute to cold stress response (Fig. 4b).
Differential expression of DREB1 regulon genes is integral to early cold stress response
The RNA-seq results from cold shock and control condition were validated using quantitative real-time PCR (qRT-PCR) for 25 genes. This included 15 genes from upregulated DEGs and 10 genes from downregulated DEGs of the cold shock transcriptome. Figure 5a shows the fold change obtained from three biological replicates for these 25 genes and their corresponding log2fold change obtained from RNA-seq data. As seen in the figure, 80% of validated genes exhibited differential expression with significant p-value (p < 0.05) that matches with the genes’ expression profiles from RNA-seq data.
In this study, we observed that the DREB1 group of genes viz., OsDREB1A, OsDREB1B, OsDREB1C, OsDREB1G, OsDREB1E, OsDREB1H were significantly upregulated under cold shock conditions. This indicates that the transcription factors which recognize the DRE motif are triggered early during cold shock, to induce downstream regulators for mounting the entire cold stress response in rice plants. Previous studies have indicated that DREB/CBF dependent regulation is considered as the major pathway in cold acclimation and is highly conserved in various plant species . Among the 516 DEGs reported in this study, 27.7% of genes were identified to have at least one DRE-binding motif at the upstream region (Additional file 5). The upregulated DREB1 regulon genes identified in this study may be grouped into four major categories: the hydrophilic proteins, LEA, DHN1 (COR410) proteins, stress associated proteins (SAP17), osmotin precursors, photoperiod sensitive and transporter proteins; signal transducing molecules, such as phosphatases, membrane kinases, Ca2+-CaM proteins, and various cold-responsive ubiquitin ligases; multiple classes of enzymes including catabolic chitinases, exocyst complex proteins and oxidoreductases, laccases; and the final group of other zinc fingers, HOX, and AP2/ERF group of transcription factors. Twenty genes from the 107 upregulated DREB regulon genes were validated, among which expression of the genes coding for AP2/ERF transcription factor OsERF102 (Os09g0457900), Calcium-binding protein OsCCD1 (Os06g0683400), and the exocyst subunit EXO70 family protein, OsEXO70FX14 (Os01g0905300) exhibited high levels of significant upregulation (Fig. 5b).
Expression of cold-responsive genes in different unexplored indica cultivars
We extended our study to ten different indica rice cultivars that were never characterized for cold stress response. These varieties are mostly high yielding varieties (HYVs), including both hybrids and field-selected varieties, cultivated all over India (Table 2). The expression profiles of 18 genes (14 upregulated and 4 downregulated, as validated in IR64) from the above-validated list were tested for this study. Our data shows upregulation of transcript levels of OsDREB1b, ONAC039, OsCML31, ERD15, OsLEA14, OsCIPK7 and the DREB regulon genes Quinone Oxidoreductase, OsFbox208, Transcription factor IIS, Pyridoxal phosphate dependent transferase, OsERF102, OsCCD1, Pectin-glucuronyltransferase, and OsEXO70FX14 genes post cold shock treatment in the majority of the varieties (Fig. 6a). In the case of the downregulated genes, OsGRF7, OsPRP1, OsbHLH155 showed consistent down-regulation in all varieties in response to cold shock treatment (Fig. 6b). Our data indicated that expression of majority of the cold-responsive genes was significantly upregulated in the CB1, and Heera cultivars under control conditions, compared to the other varieties. Interestingly, the DREB regulon genes OsEXO70FX14 (Os01g0905300), Transcription Factor IIS (Os06g0693800), Pyridoxal-phosphate dependent transferase (Os04g0614500), OsFbox208 (Os04g0414500) which exhibited highly significant upregulation under cold shock conditions in most varieties, was highly upregulated under control condition in CB1 rice variety (Fig. 6c). This higher level of expression of cold-responsive genes under control condition in CB1 and Heera varieties prompted us to examine the physiological response of the indica rice varieties used in this study.
The germination experiments performed in this study revealed that the indica rice variety CB1 remains nearly unaffected by cold stress (72 h) prior to germination, whereas, the variety IR36 shows the highest sensitivity to the same. As illustrated by Fig. 6d, the decrease in percent germination of CB1 with increasing cold stress is not significant, compared to the other varieties. The rice variety Heera also exhibits a slight reduction in the percent germination, which indicates that this variety may also have better resilience to the low-temperature stress. IR64 and IR36 varieties show a significant drop in their percent germination, which suggests their high level of sensitivity to cold stress conditions. All the other indica varieties used in this study exhibit low to moderate sensitivity to cold stress as shown by their percent germination (Fig. 6d).
The electrolyte leakage of the indica rice varieties IR64, IR36, Hamsahamas, CB1, and Heera under control and cold stress conditions (4 °C treatment for 2 h, 4 h, 24 h, 48 h, and 72 h) was studied by determining the relative electrical conductivity (REC%) of leaf tissue samples of 14 days old seedlings. As shown in Fig. 6e, REC% shows a significant increase under cold stress in IR64, IR36, and Hamsahamas varieties. In contrast, a delayed increase of REC% was observed for CB1, and Heera varieties (48 h and 72 h cold stress treatments). There was no significant change in REC% in CB1 and Heera during initial cold stress time points.
The ROS generation of the five indica varieties, under control and cold stress conditions (4 °C treatment for 2 h, 4 h, 8 h, and 24 h) was quantitated using H2DCFDA. Our data shows that the ROS generation is induced significantly under cold shock conditions in IR64, IR36 and Hamsahamas varieties, with higher levels of increase in IR64 and IR36 varieties. As exhibited in Fig. 6f, the ROS levels in these three varieties show a significant increase with increasing stress duration, up to 8 h. However, the ROS generation decreases after 24 h of cold treatment in IR64, IR36 and Hamsahamas, suggesting the induction of the ROS scavengers during later stages of cold stress. The tolerant varieties CB1, and Heera show no significant changes in the ROS levels during cold stress.
Studies have implicated cold stress as a major threat for indica rice plants in upland areas. The low temperature of the water used for irrigation in such area results in a lower rate of germination, delayed seedling emergence, increased electrolytic leakage, changes in chlorophyll fluorescence, retarded growth, lower spikelet fertility, reduced tillering and high death rate [1,2,3,4]. Therefore, despite many high yielding varieties (HYVs) being developed, the challenge remains to grow such HYVs under low-temperature conditions.
This study was designed to identify the early changes in the expression associated with cold stress response in IR64 indica rice variety. Functional annotation of upregulated DEGs generated from this study exhibited a high level of significant enrichment for GO-terms, such as transcription factor activity, response to temperature stimulus, response to osmotic stress, and response to redox changes. Analysis of genes unique to this study is centralized around stress perception and early signal relays, such as altered membrane rigidity and electrolytic leakage, the onset of calcium signaling, ROS generation and activation of stress responsive transcription factors. Interestingly, when compared with previous reports of cold stress responsive genes in indica varieties, IR64 early cold response is associated with a higher number of upregulated genes unlike late chilling stress response in 93–11 variety (4 °C for 72 h), (2298 upregulated and 3992 downregulated DEGs, ).
The role of transcription factors in stress response has been extensively studied in rice plants. Reports suggest that manipulating the expression levels of various stress induced TFs alters the expression of downstream target genes that are involved in enhancing stress tolerance mechanism in plants. In rice, while the role of WRKY family of transcription factors has been more prevalent in plant-pathogen response; AP2/EREBP and NAC/NAMTFs have gained more considerable significance in abiotic stress response studies [43,44,45]. Our data elucidates the upregulation of 38 transcription factors (10% of total upregulated DEGs) in the cold transcriptome, where AP2/EREBP family represent the dominant group. The role of DREB1 genes, a subtype of AP2/EREBP TF family, in cold stress response via the ABA-independent pathway has been well established in rice plants [4, 46,47,48]. Studies have reported that OsDREB1A and OsDREB1B gene expression was induced within 40 min after cold exposure . Interestingly, our results show similar upregulation of the DREB1 genes: OsDREB1B, OsDREB1A, OsDREB1C, OsDREB1E, OsDREB1G, and OsDREB1H genes in indica rice IR64 within 2 h of cold shock. A large number of other AP2/ERF transcription factors, such as OsERF141, OsDERF5, EREBP139, OsERF102, and OsDERF8 were also upregulated; further indicating the importance of the AP2/EREBP TFs in early cold stress response . The R2R3-type Myb TF in rice, OsMYB2 has been reported to be a master regulator for abiotic stress response conferring salt, cold and dehydration stress tolerance when overexpressed in rice . Overexpression of OsMYB4 has imparted significant tolerance to chilling and freezing stress in transgenic Arabidopsis [51,52,53]. Furthermore, Dai et al. , suggested that the R1R2R3 Myb factor, OsMYB3R-2 was upregulated under cold stress (72 h) and activated DREB/CBF pathway to increase the tolerance of the plant against freezing, drought and salinity stress . Our study also indicates upregulation of OsMYB2 and OsMYB4 genes within 2 h of cold shock, thus implying their involvement in early signaling events. ZFP182 (Os03g0820300) gene, coding for a TFIIIA-type zinc finger protein type transcription factor known to be involved in multiple abiotic stress tolerance mechanisms in rice [56, 57] was significantly upregulated under cold 2 h shock in this study (Additional file 3).
Sensing of the low temperature occurs via changes in the membrane fluidity or by sensor proteins such as RLKs, phospholipases, and Calcium channels. This perception leads to cytoskeletal reorganization due to Ca2+ influx, triggering several signaling events, thereby mounting the entire cellular response during stress [4, 58, 59]. The signaling proteins upregulated in this study involves components in calcium signaling such as Calmodulin/CBL, Ca2+ transporter and Ca2+ decoder & various other kinases, such as RLKs, RLCK, MAPKKK and Phospholipase A2 (Additional file 3). The role of these proteins is to transduce the cold signal in the cytosol and activate stress responsive transcription factors and various hydrophilic polypeptides such as LEA and Dehydrin, to stabilize against cold-stress injuries in plant cells. Increased electrolytic leakage in IR64 rice variety was evident from a significant increase in the relative electrical conductivity (REC%) for plants subjected to cold 2 h stress treatment at 4 °C (Fig. 6e).
The Late Embryogenesis Abundant (LEA) group of proteins, as the name suggests, are synthesized at the later stage of embryogenesis, prior to seed desiccation . Reports in barley, indicate that the expression of LEA genes is induced under higher ABA concentration and water-deficit conditions, asserting their role as a dehydration responsive gene . Several LEA, dehydrins and low temperature-induced  genes were upregulated in our data set (Additional file 3), suggesting the onset of dehydration response during initial cold stress conditions. Among this, OsLEA14 or wsi18 which was previously reported to have a role in chilling stress [63, 64] was also found to be present as upregulated LEA under 2 h cold shock condition.
Low-temperature stress triggers the generation of ROS like singlet oxygen species and H2O2 . This is metabolized by ROS scavenger like cytochrome P450 to prevent cellular damages, ultimately leading to redox homeostasis. Studies suggest that production of ROS varies during abiotic stress response between sensitive and tolerant varieties of rice [66,67,68]. Zhang et al. have shown that rice can better adapt to the chilling stress condition when ROS-mediated signaling genes were upregulated. Our result shows upregulation of Os71Z6, CYP78A9, OsKO4, OsABAox2, CYP701A8, and OsNR1 genes of the cytochrome family in response to cold shock. Reports suggest that the germin-like-oxalate oxidases that generate H2O2 in the apoplast  may have a role in plant defense response. Upregulation of these oxalate oxidases during cold stress response in our study agrees with the finding that ROS generation is induced as a primary abiotic stress response . Further, the significant increase in ROS production was observed under cold shock condition in IR64 variety, as observed by increased H2DCFDA fluorescence for cold 2 h treated plants as opposed to control plants (Fig. 6f). Similar to the observations from Zhang et al. study, our data shows pathway enrichment of Plant-Pathogen interaction and Diterpenoid Biosynthesis (adjusted p-values < 0.014 and 0.011 respectively) when IR64 plants were subjected to cold shock conditions (Additional file 6A). Studies in Arabidopsis report that Aquaporins or TIP1 and TIP2 are required for the transport of H2O2 species from the chloroplast and peroxisomes to the cytoplasm for regulating ROS signaling . Interestingly, our results show upregulation of both ROS transporter Aquaporin (OsTIP1) and heat shock transcription factors (HSFs) as molecular sensors of ROS, in early cold stress response to regulate the oxidation stress responsive genes  (Additional file 5). Various other cold stress response genes established by previous studies, such as OsTPP2 , OsTPS1 , OsPHS1 , were significantly upregulated during cold 2 h shock in this study.
The Dehydration Responsive Element (DRE) cis-acting elements are involved in both osmotic and cold stress induced gene expression. While osmotic stress response involves changes in DREB2 regulon gene expression, DREB1 genes are majorly responsible for mounting the cold stress response in plants . Further, several studies have reported that overexpression of the DREB1A, DREB1B, and DREB1C genes lead to certain biochemical alterations that are associated with the phenomenon of cold acclimation, and freezing stress tolerance in Arabidopsis [77,78,79,80,81]. Reports suggest that on perceiving the cold stress signal, the DREB1 cassette genes are promptly and transiently expressed, which then activates a milieu of downstream stress responsive genes in both dicotyledons and monocotyledons . In this study, significant upregulation of the DREB1 gene cassette prompted us to focus on the DREB regulon genes. Genes with the A/GCCGAC (DRE core motif) site in their 1 kb upstream sequences were screened from among the differentially expressed gene set. Upregulated DREB regulons (107 genes), obtained from this analysis include abiotic stress responsive TFs, calcium-binding proteins, RLCKs, redox- signaling molecules, and other hydrophilic water-deficit responsive proteins. The down-regulation of certain DREB regulons (36 genes) during cold shock suggests that they may be associated with normal growth and development of the plant (Additional file 5).
Cold stress, like any other abiotic stress, has an adverse effect on normal growth and metabolism in rice plants. The validated genes from the downregulated differentially expressed list are mostly related to cell cycle, protein kinases, and growth-promoting transcription factors. UDP-glucosyltransferase N-acetyltransferase (Os01g0686300) and Protein kinase (Os07g0145400) are genes encoding proteins involved in metabolic processes, such as transferring of hexosyl groups, and phosphate, respectively. These genes were significantly downregulated under cold shock, indicating hindrance to metabolism in rice plants. N-acetyltransferase ESCO1 (Os04g0498900) which codes for an acetyltransferase having a role in meiotic chromosome segregation and sister chromatid cohesion was found to be downregulated. This suggests that under cold shock inhibition of chromosomal segregation occurs that may lead to inhibition of the cell cycle. Cell cycle genes such as OsSDS (Os03g0225200) that codes for Cyclin A/B/C/D domain-containing protein and endosulphine family protein Lg106 (Os01g0249300) involved in the initiation of G0 program were both downregulated; further emphasizing the negative effect of cold stress on the cell cycle. Transcription factors significantly downregulated under cold shock included Basic helix-loop-helix dimerization region bHLH domain-containing protein, OsbHLH155 (Os06g0724800) and Transcription factor-GRAS domain-containing OsSLRL1(Os01g0646300), responsible for gibberellin (GA) signaling in plant growth and development. Transcription activator, OsGRF7 (Os12g0484900), growth-regulating factor 7, responsible for growth and development, was also downregulated under cold shock condition.
Germination experiments were employed to identify indica rice varieties with contrasting physiological response to cold stress. CB1 and Heera showed better germination profile under cold treatment, which was highly contrasted to the sensitive varieties, IR36 and IR64 exhibiting reduced rates of germination under cold stress condition (Fig. 6d and g). These varieties, along with one with semi-cold sensitive variety Hamsahamas were then subjected to relative electrical conductivity and ROS generation studies. IR64 and IR36 rice varieties exhibited a significant increase in REC% and ROS production under cold shock conditions. This indicated that electrolytic leakage and ROS generation are triggered early during cold stress response in the sensitive rice varieties, as opposed to the tolerant lines, CB1 and Heera. The low level of REC% and ROS production, together with the expression profile of the cold-responsive genes in CB1 and Heera rice varieties further validate their high level of tolerance to the low-temperature condition (Fig. 6e and f). Furthermore, transcription of some of the cold stress responsive genes and cold induced DREB regulons were found to be higher in CB1 and Heera varieties under control conditions suggesting that these varieties may be primed to mount the stress response during cold shock. However, all the DEGs that were profiled in other indica varieties show upregulation under cold shock condition (Fig. 6a to c).
Studying the cold induced change in transcription may be one approach to decipher the complexity of cold signaling in plants. However, deciphering the relationship between gene expression and dynamicity of the epigenome, changes in metabolites, alteration in protein modifications together, can lead to a better understanding of signaling network associated with cold stress response especially in cereal plants like rice.
This study was aimed at identifying some of the key responsive genes activated early under cold stress (2 h) in IR64 indica rice variety. The extensively studied cold-responsive DREB subfamily members were induced as an early event, along with other EREBP, MYB, NAC, HSF transcription factors families. A large number of Ca2+ binding proteins and kinases indicate that calcium signalling as an essential cellular mediator of the cold signal perception and response. Triggering of the ROS generation is also evident as numerous ROS producing oxidases were upregulated in the cold shock transcriptome. The growth promoter factors and transcription factors were downregulated as an initial response to the cold condition. Taken together, this study indicates that the Ca2+ and ROS mediated pathways are early cold shock-induced events which prime the cells for the later response. Physiological and expression studies with different indica rice cultivars suggest that CB1 and Heera indica rice varieties are better suited for low-temperature conditions, as opposed to the highly sensitive IR36 and IR64 varieties.
Plant growth conditions
Oryza sativa L. ssp. indica rice genotypes used in this study include IR64, CB1, Heera, Hamsahamas, Ratna, Rasi, IR36, GB1, GB3, IET4786 (Miniket), and Rajendra Bhagawati (RB) (Table 2). The seeds were surface sterilized with 0.1% (w/v) HgCl2 for 15 min, washed several times with sterile water, following which they were germinated over water-soaked sterile gauge placed in trays at 28 °C ± 1 °C in the dark for 3 days. The germinated seedlings were transferred to fresh water-soaked sterile gauge in trays, in the presence of 0.25X Murashige and Skoog complete media at 28 °C ± 1 °C in 16 h light and 8 h dark photoperiodic cycle with 50% relative humidity and 700 lmol photons m− 2 s− 1 in a plant growth chamber. For cold shock treatment, the 14-days-old seedlings (~ 100 seedlings for each experimental set) were transferred to 4 °C, whereas the control plants (~ 100 seedlings) were maintained at 28 °C ± 1 °C (marked as CS and CT respectively). For cold shock, seedling samples were collected from plants incubated at 4 °C for 2 h.
RNA was isolated from ~ 200 mg of leaf tissue of rice seedlings (n ~ 6 seedlings) for qRT-PCR from each treatment (control, CT and cold shock, CS) of IR64 cultivar (considered as a single replicate). Total RNA from each replicate was extracted using RNASure® Mini Kit (Nucleopore-Genetix), according to manufacturer’s protocol. RNA samples were treated with DNaseI to remove DNA contamination. For RNA-seq, RNA was pooled from ~ 500 mg seedlings (n ~ 15seedlings) and considered as one biological replicate. Two such replicates for each control and cold shock (CS) samples were used for sequencing. The RNA concentration was determined using the Qubit Fluorometer. An aliquot of the samples was run on an (Agilent) RNA Bioanalyzer chip to check for integrity.
cDNA library preparation and Illumina sequencing
The RNA quality check, quantification, cDNA library preparation and sequencing were done at Genotypic Technology, India using Illumina NextSeq500 platform generating 400 million paired-end reads with an average size of 75 bp.
Raw sequence processing and differential gene expression
Raw reads were assessed using FastQC . The adapters attached with the raw reads were subsequently removed using BBDuK  tool and quality was confirmed using FastQC. Processed reads were corrected using Rcorrector . HiSat2 index file was built using japonica reference genome and gtf files (gene transfer format) containing the feature list from RapDB (https://rapdb.dna.affrc.go.jp/download/archive/irgsp1/IRGSP-1.0_genome.fasta.gz) and (https://rapdb.dna.affrc.go.jp/download/archive/irgsp1/IRGSP-1.0_representative_2019-03-22.tar.gz). Better and complete annotation of the Japonica rice was the rationale behind using it as the reference genome . HiSat2  was used to align corrected reads with the reference genome. The alignment output file was in SAM format and was converted into BAM format using SAMtools . Qualimap  analysis was done to assess the quality of alignment on sorted BAM files. FeatureCounts  was used to convert sorted BAM file to count reads. Later, DeSeq2  that takes count files as input was used for calculating differentially expressed genes. In order to filter out differentially expressed genes, threshold fold-change was set at 1.5 with a p-value cut off of ≤0.05. Subsequently, unsupervised hierarchical clustering and analysis were done on analyzed output.
The BLAST2GO [92, 93] program was used against NR database for the GO annotation analysis of DEGs (Differentially Expressed Genes) in terms of Biological Process, Molecular Function and Cellular Components ontologies. WEGO tool was used to visualize the GO annotation, where both the upregulated and downregulated GO data sets were uploaded to obtain a comparative GO analysis. KO (KEGG Orthology) analysis was performed by submitting the upregulated gene list in BLASTKOALA . PANTHER  classification system was used to determine the protein family classification. Significantly enriched GO and KEGG Pathways of the differentially expressed genes was carried out using ShinyGO v0.60  online tool and Cluster profiler package .
Further, biologically relevant and statistically significant enriched genes were identified, and regulatory network modelling analysis was carried out using in-house scripts, and the output was visualized using Cytoscape V2.8.3 . The p-value calculated by the hypergeometric test and was corrected by FDR. The FDR value of ≤0.05 was used as the threshold to identify the significant functional categories and metabolic pathways. Upregulated genes including transcription factors, kinases, calcium signaling components, redox components and other stress responsive genes were filtered out from the upregulated gene list and submitted to STRING  database for network analysis.
DRE motif searching in upstream of genes
To search for the DRE (Dehydration Responsive Element) motifs, DNA sequences were fetched within 1 KB upstream (from the transcription start site) of the coding genes using getfasta program (bedtools) . The extracted region sequences of the differentially regulated genes were and examined for the DRE core-motif (A/GCCGAC) . The Find Individual Motif Occurrences (FIMO) tool of the MEME package was used for motif search .
For qRT-PCR, RNA was isolated from different samples using TRIzol reagent (Invitrogen) as described in manufacturers’ protocol. cDNA samples were generated using 5 μg of total RNA from three biological replicates. The cDNA thus generated was used for subsequent validation experiments. Forty-five genes (35 upregulated and 10 downregulated) were selected from the differentially regulated gene list, based on the potential role in cold stress response, presence of DRE sites and literature study. The primers used for this study are listed in Additional file 7. OsActin1(Os03g0718100) and OsUbq5 (Os01g0328400) genes were selected as endogenous control genes. All reactions were performed in three independent biological replicates, and the expression levels for each sample were calculated using the ΔCt method. Two-tailed paired t-test was performed to determine the level of significance. In case of the expression studies in different indica varieties (Fig. 6a to c), the 2-ΔCt values were used to generate heatmaps for each gene, under control and 2 h cold shock conditions.
Percent germination under cold stress condition
Surface sterilized seeds of the nine indica rice varieties (IR64, IR36, Hamsahamas, Rasi, Ratna, Heera, CB1, GB1, and GB3) were incubated at 4 °C for varying time points: 24 h, 48 h, and 72 h under moist conditions. Following each cold stress time points, the seeds were brought to 28 °C ± 1 °C and incubated for 3 days to mimic the control germination condition. Finally, ten seeds of each variety were plated in separate sterile Petri plates on adequately moist 90 mm filter paper disc. The germination rate, root length, and shoot length of each seed was recorded for the next 5 days. Three biological replicates were generated for each variety, and two-tailed paired t-test was performed to determine the level of significance. The percent germination of the varieties was calculated as under:
Determining the electrolyte leakage (using REC assay) and ROS production in indica varieties
The relative electrical conductivity (REC) of five indica rice varieties (IR36, IR64, Hamsahamas, CB1, and Heera) was measured under cold stress condition, to determine the variation in the physiological response of these varieties. For this, 14 days old seedlings of each variety were subjected to 4 °C treatment for 2 h, 4 h, 6 h, 24 h, 48 h, and 72 h respectively. Seedlings maintained at 28 °C ± 1 °C were used as the control set. Two hundred milligrams of leaf tissue from each set (after stress treatment, and control) was harvested and immersed in tubes containing 20 mL of distilled water. The tubes were then incubated at room temperature (25 °C) with constant shaking at 120 rpm  for 2 h. After incubation, the initial electrical conductivity of the solution (EC1) was measured, which represents the ion leakage from the leaf tissue samples. Following measurement of EC1, the solution was heated to a temperature of 100 °C for 30mins, cooled to room temperature, and the electrical conductivity of the solution (EC2) was measured at room temperature. The electrical conductivity of distilled water, ECw1 and ECw2 were measured for normalization purpose . The relative electrolyte content (REC) was calculated as under:
The ROS production of IR36, IR64, Hamsahamas, CB1, and Heera indica rice varieties was determined under control and cold stress treatment (at 4 °C for 2 h, 4 h, 8 h, and 24 h) conditions. Ground tissue powder from 100 mg leaf tissue samples of 14 days old seedlings of each variety was suspended in 10 mM Tris-Cl (pH = 7.2). After removing the cellular debris, the plant extract (was diluted?) was subjected to ROS measurement using H2DCFDA (2′7’- Dichlorofluorescein diacetate: 100 mM solution in DMSO). The protein concentration of each sample was determined using the Bradford reagent. BSA standard curve was plotted to determine the protein concentration. The ROS generation of the samples was expressed as fluorescence units/mg of protein .
For both the physiological experiments, three biological replicates were generated for each variety, and two-tailed paired t-test was performed to determine the level of significance.
Availability of data and materials
The datasets generated and/or analysed during this current study are available in the NCBI Sequence Read Archive repository (https://www.ncbi.nlm.nih.gov/bioproject/?term=PRJNA506503) under the following accession numbers: SRX5055383 (IR64_cold 2 h_replicate_1), SRX5055384 (IR64_cold 2 h_replicate_2), SRX5055385 (IR64_control_replicate_1), and SRX5055386 (IR64_control_replicate_2).
Cold shock treated seedlings
Differentially expressed genes
Dehydration responsive element
Find Individual Motif Occurrences
High yielding varieties
Reactive oxygen species
Quantitative real-time PCR
Kaneda C. Response of indica japonica rice hybrids to low temperature. SABRAO J. 1974;6:17–32.
Mackill DJ, Lei X. Genetic variation for traits related to temperate adaptation of rice cultivars. Crop Sci. 1997;37(4):1340–6.
Nakagahra M, Okuno K, Vaughan D. Rice genetic resources: history, conservation, investigative characterization and use in Japan. Plantmolecular biology; 1997;35(1-2):69–77.
Zhang Q, Chen Q, Wang S, Hong Y, Wang Z. Rice and cold stress: methods for its evaluation and summary of cold tolerance-related quantitative trait loci. Rice. 2014;7(1):24.
Los DA, Murata N. Membrane fluidity and its roles in the perception of environmental signals. Biochim Biophys Acta. 2004;1666(1–2):142–57.
Huang J, Sun S-J, Xu D-Q, Yang X, Bao Y-M, Wang Z-F, Tang H-J, Zhang H. Increased tolerance of rice to cold, drought and oxidative stresses mediated by the overexpression of a gene that encodes the zinc finger protein ZFP245. Biochem Biophys Res Commun. 2009;389(3):556–61.
Yun K-Y, Park MR, Mohanty B, Herath V, Xu F, Mauleon R, Wijaya E, Bajic VB, Bruskiewich R, de los Reyes BG. Transcriptional regulatory network triggered by oxidative signals configures the early response mechanisms of japonica rice to chilling stress. BMC Plant Biol. 2010;10(1):16.
Tian Y, Zhang H, Pan X, Chen X, Zhang Z, Lu X, Huang R. Overexpression of ethylene response factor TERF2 confers cold tolerance in rice seedlings. Transgenic Res. 2011;20(4):857–66.
Zhang Q, Jiang N, Wang G-L, Hong Y, Wang Z. Advances in understanding cold sensing and the cold-responsive network in rice. Adv Crop Sci Tech. 2013;1(1):104.
Chinnusamy V, Zhu J, Zhu JK. Gene regulation during cold acclimation in plants. Physiol Plant. 2006;126(1):52–61.
Reddy AS, Ali GS, Celesnik H, Day IS. Coping with stresses: roles of calcium-and calcium/calmodulin-regulated gene expression. Plant Cell. 2011;23(6):2010–32.
Liu J, Shi Y, Yang S. Insights into the regulation of C-repeat binding factors in plant cold signaling. J Integr Plant Biol. 2018;60(9):780–95.
Cook D, Fowler S, Fiehn O, Thomashow MF. A prominent role for the CBF cold response pathway in configuring the low-temperature metabolome of Arabidopsis. Proc Natl Acad Sci. 2004;101(42):15243–8.
Guy CL. Cold acclimation and freezing stress tolerance: role of protein metabolism. Annu Rev Plant Biol. 1990;41(1):187–223.
Hannah MA, Heyer AG, Hincha DK. A global survey of gene regulation during cold acclimation in Arabidopsis thaliana. PLoS Genet. 2005;1(2):e26.
Maruyama K, Takeda M, Kidokoro S, Yamada K, Sakuma Y, Urano K, Fujita M, Yoshiwara K, Matsukura S, Morishita Y. Metabolic pathways involved in cold acclimation identified by integrated analysis of metabolites and transcripts regulated by DREB1A and DREB2A. Plant Physiol. 2009;150(4):1972–80.
Park S, Lee CM, Doherty CJ, Gilmour SJ, Kim Y, Thomashow MF. Regulation of the Arabidopsis CBF regulon by a complex low-temperature regulatory network. Plant J. 2015;82(2):193–207.
Kim SH, Kim HS, Bahk S, An J, Yoo Y, Kim J-Y, Chung WS. Phosphorylation of the transcriptional repressor MYB15 by mitogen-activated protein kinase 6 is required for freezing tolerance in Arabidopsis. Nucleic Acids Res. 2017;45(11):6613–27.
Zhu J, Shi H, Lee B-H, Damsz B, Cheng S, Stirm V, Zhu J-K, Hasegawa PM, Bressan RA. An Arabidopsis homeodomain transcription factor gene, HOS9, mediates cold tolerance through a CBF-independent pathway. Proc Natl Acad Sci. 2004;101(26):9873–8.
Lv Y, Guo Z, Li X, Ye H, Li X, Xiong L. New insights into the genetic basis of natural chilling and cold shock tolerance in rice by genome-wide association analysis. Plant Cell Environ. 2016;39(3):556–70.
Cheng C, Yun K-Y, Ressom HW, Mohanty B, Bajic VB, Jia Y, Yun SJ, de los Reyes BG. An early response regulatory cluster induced by low temperature and hydrogen peroxide in seedlings of chilling-tolerant japonica rice. BMC Genomics. 2007;8(1):175.
Zhang, Q., Chen, Q., Wang, S., Hong, Y., & Wang, Z. Rice and cold stress: methods for its evaluation and summary of cold tolerance-relatedquantitative trait loci. New York: Rice; 2014;7(1):24. https://doi.org/10.1186/s12284-014-0024-3.
Lan Thi Hoang X, Du Nhi NH, Binh Anh Thu N, Phuong Thao N, Phan Tran L-S. Transcription factors and their roles in signal transduction in plants under abiotic stresses. Curr Genomics. 2017;18(6):483–97.
Agarwal PK, Agarwal P, Reddy M, Sopory SK. Role of DREB transcription factors in abiotic and biotic stress tolerance in plants. Plant Cell Rep. 2006;25(12):1263–74.
Iuchi S, Kobayashi M, Taji T, Naramoto M, Seki M, Kato T, Tabata S, Kakubari Y, Yamaguchi-Shinozaki K, Shinozaki K. Regulation of drought tolerance by gene manipulation of 9-cis-epoxycarotenoid dioxygenase, a key enzyme in abscisic acid biosynthesis in Arabidopsis. Plant J. 2001;27(4):325–33.
Kulik A, Wawer I, Krzywińska E, Bucholc M, Dobrowolska G. SnRK2 protein kinases—key regulators of plant response to abiotic stresses. Omics. 2011;15(12):859–72.
Jeong D-H, Green PJ. The role of rice microRNAs in abiotic stress responses. J Plant Biol. 2013;56(4):187–97.
Yamada M, Morishita H, Urano K, Shiozaki N, Yamaguchi-Shinozaki K, Shinozaki K, Yoshiba Y. Effects of free proline accumulation in petunias under drought stress. J Exp Bot. 2005;56(417):1975–81.
Hussain Wani S, Brajendra Singh N, Haribhushan A, Iqbal Mir J. Compatible solute engineering in plants for abiotic stress tolerance-role of glycine betaine. Curr Genomics. 2013;14(3):157–65.
Slama I, Abdelly C, Bouchereau A, Flowers T, Savoure A. Diversity, distribution and roles of osmoprotective compounds accumulated in halophytes under abiotic stress. Ann Bot. 2015;115(3):433–47.
Guo H, Chen C, Zhang X, Qin X, Feng R. Study on cold tolerance identification and genetic homozygosity of wild rice in Guangxi. Southwest China J Agric Sci. 2017;30(6):1245–50.
Mao D, Yu L, Chen D, Li L, Zhu Y, Xiao Y, Zhang D, Chen C. Multiple cold resistance loci confer the high cold tolerance adaptation of Dongxiang wild rice (Oryza rufipogon) to its high-latitude habitat. Theor Appl Genet. 2015;128(7):1359–71.
Liu J, Cheng J, Gong J, Kang M, Hong Y. Cloning of the DREB-like transcription factor from Chaling wild rice and the construction of plant expression vector. J Hunan Agric Univ. 2010;36(6):630–3.
Dazhou C, Pingan Z, Yeqing X. Identification of QTLs for cold tolerance at seedling stage in Dongxiang wild rice (Oryza rufipogon Griff.) by SSR markers. China Jiangxi Nongxue Daoxue Xuebao . 2004;24(6):753–56.
Li F, Guo S, Zhao Y, Chen D, Chong K, Xu Y. Overexpression of a homopeptide repeat-containing bHLH protein gene (OrbHLH001) from Dongxiang wild rice confers freezing and salt tolerance in transgenic Arabidopsis. Plant Cell Rep. 2010;29(9):977–86.
Zhang F, Huang L, Wang W, Zhao X, Zhu L, Fu B, Li Z. Genome-wide gene expression profiling of introgressed indica rice alleles associated with seedling cold tolerance improvement in a japonica rice background. BMC Genomics. 2012;13(1):461.
Udvardi MK, Kakar K, Wandrey M, Montanari O, Murray J, Andriankaja A, Zhang J-Y, Benedito V, Hofer JM, Chueng F. Legume transcription factors: global regulators of plant development and response to the environment. Plant Physiol. 2007;144(2):538–49.
Riechmann JL, Heard J, Martin G, Reuber L, Jiang C-Z, Keddie J, Adam L, Pineda O, Ratcliffe O, Samaha R. Arabidopsis transcription factors: genome-wide comparative analysis among eukaryotes. Science. 2000;290(5499):2105–10.
Lata C, Yadav A, Prasad M. Role of plant transcription factors in abiotic stress tolerance. Abiotic Stress Response in Plants, Physiological, Biochemical and Genetic Perspectives, Arun Shanker and B. Venkateswarlu, INTECH Open Access Publishers; 2011. https://doi.org/10.5772/23172.
Dubouzet JG, Sakuma Y, Ito Y, Kasuga M, Dubouzet EG, Miura S, Seki M, Shinozaki K, Yamaguchi-Shinozaki K. OsDREB genes in rice, Oryza sativa L., encode transcription activators that function in drought-, high-salt-and cold-responsive gene expression. Plant J. 2003;33(4):751–63.
Lata C, Prasad M. Role of DREBs in regulation of abiotic stress responses in plants. J Exp Bot. 2011;62(14):4731–48.
Shen C, Li D, He R, Fang Z, Xia Y, Gao J, Shen H, Cao M. Comparative transcriptome analysis of RNA-seq data for cold-tolerant and cold-sensitive rice genotypes under cold stress. J Plant Biol. 2014;57(6):337–48.
Khong GN, Richaud F, Coudert Y, Pati PK, Santi C, Périn C, Breitler J-C, Meynard D, Vinh DN, Guiderdoni E. Modulating rice stress tolerance by transcription factors. Biotechnol Genet Eng Rev. 2008;25(1):381–404.
Pandey SP, Somssich IE. The role of WRKY transcription factors in plant immunity. Plant Physiol. 2009;150(4):1648–55.
Viana VE, Busanello C, da Maia LC, Pegoraro C, de Oliveira AC. Activation of rice WRKY transcription factors: an army of stress fighting soldiers? Curr Opin Plant Biol. 2018;45:268–75.
Liu Q, Zhao N, Yamaguch-Shinozaki K, Shinozaki K. Regulatory role of DREB transcription factors in plant drought, salt and cold tolerance. Chin Sci Bull. 2000;45(11):970–5.
Lee S-C, Huh K-W, An K, An G, Kim S-R. Ectopic expression of a cold-inducible transcription factor, CBF1/DREB1b, in transgenic rice (Oryza sativa L.). Mol Cells (Springer Science & Business Media BV). 2004;18(1):107–14.
Roy D, Paul A, Roy A, Ghosh R, Ganguly P, Chaudhuri S. Differential acetylation of histone H3 at the regulatory region of OsDREB1b promoter facilitates chromatin remodelling and transcription activation during cold stress. PLoS One. 2014;9(6):e100343.
Mizoi J, Shinozaki K, Yamaguchi-Shinozaki K. AP2/ERF family transcription factors in plant abiotic stress responses. Biochim Biophys Acta. 2012;1819(2):86–96.
Yang A, Dai X, Zhang W-H. A R2R3-type MYB gene, OsMYB2, is involved in salt, cold, and dehydration tolerance in rice. J Exp Bot. 2012;63(7):2541–56.
Vannini C, Locatelli F, Bracale M, Magnani E, Marsoni M, Osnato M, Mattana M, Baldoni E, Coraggio I. Overexpression of the rice Osmyb4 gene increases chilling and freezing tolerance of Arabidopsis thaliana plants. Plant J. 2004;37(1):115–27.
PARK MR, YUN KY, Mohanty B, Herath V, Xu F, Wijaya E, Bajic VB, YUN SJ, De Los Reyes BG. Supra-optimal expression of the cold-regulated OsMyb4 transcription factor in transgenic rice changes the complexity of transcriptional network with major effects on stress tolerance and panicle development. Plant Cell Environ. 2010;33(12):2209–30.
Soltész A, Vágújfalvi A, Rizza F, Kerepesi I, Galiba G, Cattivelli L, Coraggio I, Crosatti C. The rice Osmyb4 gene enhances tolerance to frost and improves germination under unfavourable conditions in transgenic barley plants. J Appl Genet. 2012;53(2):133–43.
Dai X, Xu Y, Ma Q, Xu W, Wang T, Xue Y, Chong K. Overexpression of an R1R2R3 MYB gene, OsMYB3R-2, increases tolerance to freezing, drought, and salt stress in transgenic Arabidopsis. Plant Physiol. 2007;143(4):1739–51.
Ma Q, Dai X, Xu Y, Guo J, Liu Y, Chen N, Xiao J, Zhang D, Xu Z, Zhang X. Enhanced tolerance to chilling stress in OsMYB3R-2 transgenic rice is mediated by alteration in cell cycle and ectopic expression of stress genes. Plant Physiol. 2009;150(1):244–56.
Huang J, Yang X, Wang M-M, Tang H-J, Ding L-Y, Shen Y, Zhang H-S. A novel rice C2H2-type zinc finger protein lacking DLN-box/EAR-motif plays a role in salt tolerance. Biochim Biophys Acta. 2007;1769(4):220–7.
Huang J, Sun S, Xu D, Lan H, Sun H, Wang Z, Bao Y, Wang J, Tang H, Zhang H. A TFIIIA-type zinc finger protein confers multiple abiotic stress tolerances in transgenic rice (Oryza sativa L.). Plant Mol Biol. 2012;80(3):337–50.
Solanke AU, Sharma AK. Signal transduction during cold stress in plants. Physiol Mol Biol Plants. 2008;14(1–2):69–79.
Guo X, Liu D, Chong K. Cold signaling in plants: insights into mechanisms and regulation. J Integr Plant Biol. 2018;60(9):745–56.
Ingram J, Bartels D. The molecular basis of dehydration tolerance in plants. Annu Rev Plant Biol. 1996;47(1):377–403.
Hong B, Uknes SJ, Ho T-hD. Cloning and characterization of a cDNA encoding a mRNA rapidly-induced by ABA in barley aleurone layers. Plant Mol Biol. 1988;11(4):495–506.
Vergnolle C, Vaultier M-N, Taconnat L, Renou J-P, Kader J-C, Zachowski A, Ruelland E. The cold-induced early activation of phospholipase C and D pathways determines the response of two distinct clusters of genes in Arabidopsis cell suspensions. Plant Physiol. 2005;139(3):1217–33.
Takahashi R, Joshee N, Kitagawa Y. Induction of chilling resistance by water stress, and cDNA sequence analysis and expression of water stress-regulated genes in rice. Plant Mol Biol. 1994;26(1):339–52.
Joshee N, Kisaka H, Kitagawa Y. Isolation and characterization of a water stress-specific genomic gene, pwsi 18, from rice. Plant Cell Physiol. 1998;39(1):64–72.
Habibi G. Hydrogen peroxide (H2O2) generation, scavenging and signaling in plants. In: Oxidative damage to plants. AcademicPress; 2014. p. 557–84.
Kaur N, Dhawan M, Sharma I, Pati PK. Interdependency of reactive oxygen species generating and scavenging system in salt sensitive and salt tolerant cultivars of rice. BMC Plant Biol. 2016;16(1):131.
Zhang J, Luo W, Zhao Y, Xu Y, Song S, Chong K. Comparative metabolomic analysis reveals a reactive oxygen species-dominated dynamic model underlying chilling environment adaptation and tolerance in rice. New Phytol. 2016;211(4):1295–310.
Ijaz B, Formentin E, Ronci B, Locato V, Barizza E, Hyder MZ, Schiavo FL, Yasmin T. Salt tolerance in indica rice cell cultures depends on a fine tuning of ROS signalling and homeostasis. PLoS One. 2019;14(4):e0213986.
Mittler R. Oxidative stress, antioxidants and stress tolerance. Trends Plant Sci. 2002;7(9):405–10.
Choudhury FK, Rivero RM, Blumwald E, Mittler R. Reactive oxygen species, abiotic stress and stress combination. Plant J. 2017;90(5):856–67.
Bienert GP, Møller AL, Kristiansen KA, Schulz A, Møller IM, Schjoerring JK, Jahn TP. Specific aquaporins facilitate the diffusion of hydrogen peroxide across membranes. J Biol Chem. 2007;282(2):1183–92.
Miller G, Mittler R. Could heat shock transcription factors function as hydrogen peroxide sensors in plants? Ann Bot. 2006;98(2):279–88.
Song S-Y, Chen Y, Chen J, Dai X-Y, Zhang W-H. Physiological mechanisms underlying OsNAC5-dependent tolerance of rice plants to abiotic stress. Planta. 2011;234(2):331–45.
Li H-W, Zang B-S, Deng X-W, Wang X-P. Overexpression of the trehalose-6-phosphate synthase gene OsTPS1 enhances abiotic stress tolerance in rice. Planta. 2011;234(5):1007–18.
Du H, Wu N, Chang Y, Li X, Xiao J, Xiong L. Carotenoid deficiency impairs ABA and IAA biosynthesis and differentially affects drought and cold tolerance in rice. Plant Mol Biol. 2013;83(4–5):475–88.
Yamaguchi-Shinozaki K, Shinozaki K. DREB regulons in abiotic-stress-responsive gene expression in plants. In Molecular breeding of forage and turf. New York: Springer; 2009. p. 15–28.
Gilmour SJ, Fowler SG, Thomashow MF. Arabidopsis transcriptional activators CBF1, CBF2, and CBF3 have matching functional activities. Plant Mol Biol. 2004;54(5):767–81.
Gilmour SJ, Sebolt AM, Salazar MP, Everard JD, Thomashow MF. Overexpression of the Arabidopsis CBF3transcriptional activator mimics multiple biochemical changes associated with cold acclimation. Plant Physiol. 2000;124(4):1854–65.
Jaglo-Ottosen KR, Gilmour SJ, Zarka DG, Schabenberger O, Thomashow MF. Arabidopsis CBF1 overexpression induces COR genes and enhances freezing tolerance. Science. 1998;280(5360):104–6.
Kasuga M, Liu Q, Miura S, Yamaguchi-Shinozaki K, Shinozaki K. Improving plant drought, salt, and freezing tolerance by gene transfer of a single stress-inducible transcription factor. Nat Biotechnol. 1999;17(3):287–91.
Liu Q, Kasuga M, Sakuma Y, Abe H, Miura S, Yamaguchi-Shinozaki K, Shinozaki K. Two transcription factors, DREB1 and DREB2, with an EREBP/AP2 DNA binding domain separate two cellular signal transduction pathways in drought-and low-temperature-responsive gene expression, respectively, in Arabidopsis. Plant Cell. 1998;10(8):1391–406.
Nakashima K, Yamaguchi-Shinozaki K. Regulons involved in osmotic stress-responsive and cold stress-responsive gene expression in plants. Physiol Plant. 2006;126(1):62–71.
Andrews S. FastQC: a quality control tool for high throughput sequence data. In: Babraham Bioinformatics. Cambridge: Babraham Institute; 2010.
Bushnell B. BBMap: a fast, accurate, splice-aware aligner. Berkeley: Lawrence Berkeley National Lab.(LBNL); 2014.
Song L, Florea L. Rcorrector: efficient and accurate error correction for Illumina RNA-seq reads. GigaScience. 2015;4(1):48.
Sinha S, Sevanthi V, Chaudhary S, Tyagi P, Venkadesan S, Rani M, Mandal P. Transcriptome analysis of two rice varieties contrasting for nitrogen use efficiency under chronic N starvation reveals differences in chloroplast and starch metabolism-related genes. Genes. 2018;9(4):206.
Kim D, Langmead B, Salzberg SL. HISAT: a fast spliced aligner with low memory requirements. Nat Methods. 2015;12(4):357.
Li H, Handsaker B, Wysoker A, Fennell T, Ruan J, Homer N, Marth G, Abecasis G, Durbin R. The sequence alignment/map format and SAMtools. Bioinformatics. 2009;25(16):2078–9.
Okonechnikov K, Conesa A, García-Alcalde F. Qualimap 2: advanced multi-sample quality control for high-throughput sequencing data. Bioinformatics. 2015;32(2):292–4.
Liao Y, Smyth GK, Shi W. The R package Rsubread is easier, faster, cheaper and better for alignment and quantification of RNA sequencing reads. Nucleic Acids Res. 2019;47(8):e47.
Love MI, Huber W, Anders S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 2014;15(12):550.
Conesa A, Götz S. Blast2GO: a comprehensive suite for functional analysis in plant genomics. Int J Plant Genomics. 2008;2008:619832.
Conesa A, Götz S, García-Gómez JM, Terol J, Talón M, Robles M. Blast2GO: a universal tool for annotation, visualization and analysis in functional genomics research. Bioinformatics. 2005;21(18):3674–6.
Kanehisa M, Sato Y, Morishima K. BlastKOALA and GhostKOALA: KEGG tools for functional characterization of genome and metagenome sequences. J Mol Biol. 2016;428(4):726–31.
Mi H, Huang X, Muruganujan A, Tang H, Mills C, Kang D, Thomas PD. PANTHER version 11: expanded annotation data from Gene Ontology and Reactome pathways, and data analysis tool enhancements. Nucleic Acids Res. 2016;45(D1):D183–9.
Ge S, Jung D. ShinyGO: a graphical enrichment tool for animals and plants. bioRxiv. 2018:315150.
Yu G, Wang L-G, Han Y, He Q-Y. clusterProfiler: an R package for comparing biological themes among gene clusters. Omics. 2012;16(5):284–7.
Smoot ME, Ono K, Ruscheinski J, Wang P-L, Ideker T. Cytoscape 2.8: new features for data integration and network visualization. Bioinformatics. 2010;27(3):431–2.
Szklarczyk D, Gable AL, Lyon D, Junge A, Wyder S, Huerta-Cepas J, Simonovic M, Doncheva NT, Morris JH, Bork P. STRING v11: protein–protein association networks with increased coverage, supporting functional discovery in genome-wide experimental datasets. Nucleic Acids Res. 2018;47(D1):D607–13.
Quinlan AR, Hall IM. BEDTools: a flexible suite of utilities for comparing genomic features. Bioinformatics. 2010;26(6):841–2.
Yamaguchi-Shinozaki K, Shinozaki K. A novel cis-acting element in an Arabidopsis gene is involved in responsiveness to drought, low-temperature, or high-salt stress. Plant Cell. 1994;6(2):251–64.
Bailey TL, Boden M, Buske FA, Frith M, Grant CE, Clementi L, Ren J, Li WW, Noble WS. MEME SUITE: tools for motif discovery and searching. Nucleic Acids Res. 2009;37(suppl_2):W202–8.
Guan S, Xu Q, Ma D, Zhang W, Xu Z, Zhao M, Guo Z. Transcriptomics profiling in response to cold stress in cultivated rice and weedy rice. Gene. 2019;685:96–105.
Jambunathan N. Determination and detection of reactive oxygen species (ROS), lipid peroxidation, and electrolyte leakage in plants. In Plant stress tolerance. Humana Press; 2010. p. 291–97.
The authors extend gratitude to Dr. Somnath Bhattacharya (Professor, Genetics and Plant Breeding, BCKV, Mohanpur, West Bengal, India) and Dr. Sukanta Dasgupta (Deputy Director of Agriculture (WBP), Government of West Bengal, India) for kindly supplying the indica varieties, used in this study. Dr. Tripathy and Abhishek Das would like to acknowledge CSIR-IICB for infrastructure for data analysis.
This work was supported by SERB, Department of Science and Technology, Government of India, (SERB/SB/SO/PS/80/2013). Pratiti Dasgupta sincerely acknowledges CSIR; Government of India, for her fellowship [09/015(0490)/2016-EMR-I]. The authors sincerely acknowledge Bose Institute for institutional support. Abhishek Das sincerely acknowledges ICMR; Government of India, for his fellowship. Contents of this manuscript are solely the responsibility of the authors and do not necessarily represent the official views of the funding agency. The funders have no role in designing or data collection and analysis, decision to publish, or preparation of the manuscript.
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Table containing detailed analysis of cold shock (2 h at 4 °C) significantly upregulated and downregulated genes, compared to control conditions, in IR64 seedlings.
Sheet 1 contains a table for GO enrichment, genes for each functional category for upregulated DEGs. Sheet 2 contains a table showing KEGG Brite analysis for cold induced DEGs.
Gene details for differentially regulated Transcription factors, components of Redox pathway, and other cold-responsive genes.
Differentially expressed genes unique to IR64 cold shock 2 h transcriptome data.
Differentially regulated gene list with DRE-core motif in 1 kb upstream of the transcription start site.
Functional annotation for differentially regulated genes. (A) KEGG pathway enrichment data for upregulated DEGs (B) Top ten abundant domains present in the upregulated genes, obtained using Blast2GO (C) the pathway reconstruction result of upregulated differentially expressed genes, based on KEGG, generated using BLASTKOALA.
Table containing sequences of primers used for qRT-PCR assays.)
About this article
Cite this article
Dasgupta, P., Das, A., Datta, S. et al. Understanding the early cold response mechanism in IR64 indica rice variety through comparative transcriptome analysis. BMC Genomics 21, 425 (2020). https://doi.org/10.1186/s12864-020-06841-2