Transcriptome Analysis of the Responses of Rice Leaves to Chilling and Subsequent Recovery

Abstract

Improving chilling tolerance at the seedling stage in rice is essential for agricultural research. We combined a physiological analysis with transcriptomics in a variety Dular subjected to chilling followed by recovery at normal temperature to better understand the chilling tolerance mechanisms of rice. Chilling inhibited the synthesis of chlorophyll and non-structural carbohydrate (NSC) and disrupted the ion balance of the plant, resulting in the impaired function of rice leaves. The recovery treatment can effectively reverse the chilling-related injury. Transcriptome results displayed that 21,970 genes were identified at three different temperatures, and 11,732 genes were differentially expressed. According to KEGG analysis, functional categories for differentially expressed genes (DEGs) mainly included ribosome (8.72%), photosynthesis–antenna proteins (7.38%), phenylpropanoid biosynthesis (11.41%), and linoleic acid metabolism (10.07%). The subcellular localization demonstrated that most proteins were located in the chloroplasts (29.30%), cytosol (10.19%), and nucleus (10.19%). We proposed that some genes involved in photosynthesis, ribosome, phenylpropanoid biosynthesis, and linoleic acid metabolism may play key roles in enhancing rice adaptation to chilling stress and their recovery capacity. These findings provide a foundation for future research into rice chilling tolerance mechanisms.

1. Introduction

Temperature is one of the leading environmental factors affecting plant distribution in terrestrial ecosystems [1]. Therefore, crop losses caused by chilling are an important worldwide problem. The plant damage caused by chilling has become more serious with the continuous deterioration of the global ecological environment. Chilling weakens the cell membrane, resulting in reactive oxygen species (ROS) and other harmful substances, and impacts plant vegetative growth and yield formation [2,3,4]. Several studies have focused on the molecular response of rice to chilling and its physiological and ecological characteristics [5,6,7]. Studies have demonstrated that the growth index and the physiological enzyme activity of rice are significantly inhibited under chilling, although this effect becomes smaller after the chilling pretreatment of rice [8]. This finding indicated that the chilling tolerance of rice could be induced. Because of this, studies have found that chilling environments can induce the expression of multiple genes [9,10,11,12]. The proteins encoded by these genes are involved in the sensing and signal transduction processes of plant tolerance to chilling to promote the synthesis of plant osmotic regulators [13]. They also enhance the activity of antioxidant enzymes [14] and improve cell membrane fluidity [15], thus enhancing the plant’s ability to withstand chilling [16,17]. Therefore, revealing the process and mechanism of sensing, transmitting chilling signals, and stimulating chilling tolerance of plants is of great theoretical and practical significance for further understanding the chilling tolerance mechanism of plants and accelerating the molecular breeding process in crops with chilling tolerance.

Existing studies have revealed that the transcriptional regulation of chilling response genes is an important process in plant responses to chilling injury [18]. Among these, the signal dependent on CBF/DREB, is an essential signaling pathway for Arabidopsis thaliana and rice to produce chilling tolerance ability [7]. The calcium signaling pathway is activated after the plant accumulates chilling tolerance characteristics, and the genes involved are evolutionarily conserved [16]. Through in-depth research, Ma et al. found that COLD1, a key QTL gene in rice that senses chilling, significantly enhanced the chilling tolerance of japonica rice with the overexpression of this gene, whereas functional deletion mutant COLD-1 and antisense gene lines were susceptible to chilling damage [6]. Zhang et al. found that three genes (LOC_Os01g55510, LOC_Os01g55350, and LOC_Os01g55560) in chromosome 1 of rice were closely associated with chilling tolerance [19], indicating the complexity of plant tolerance to chilling. In addition, several transcription factors have been identified, including MYBS3, OsWRKY71, and OsWRKY76, which play important roles in the tolerance of rice to chilling [20,21,22]. Recently, the role of calcium-dependent protein kinases in the chilling tolerance of rice has been widely reported [21,23,24]. These results provided important theoretical references and a scientific basis to further study and understand the molecular mechanisms and regulatory network of rice tolerance to chilling.

Improving the cold tolerance of rice seedlings under low-temperature conditions is very important in promoting rice growth. The use of biostimulants to counteract the effect of abiotic stress is well-recognized [25,26]. However, biostimulants are plant extracts and their bioactive compounds and tolerance mechanisms are unknown [27,28]. The most common components of the biostimulants are mineral elements, humic substances (HSs), vitamins, amino acids, chitin, chitosan, and poly- and oligosaccharides [29]. In addition, biological stimulants positively affect enzyme activity and gene expression in plants’ primary and secondary metabolism [29]. Therefore, biostimulins can improve plant tolerance to extreme temperature, which is an important means and hope for alleviating cold damage in sustainable agriculture [30,31].

Dular was a landrace aus variety that originated from India [32,33]. Dular exhibited significant chilling injury at temperatures below 15 °C [34,35]. The damage to plants under chilling is closely related to the intensity and time of stress. The tolerance of plants to external stress has a limit. The damage caused by short-term stress is reversible, and plants can adapt and return to normal growth. If the stress time is too long or the level of stress is too high, irreversible damage occurs, damaging or killing the plant. The cold spell in late spring is a common but unavoidable climate phenomenon in China. The sudden drop in temperature in a short period causes serious harm to crop growth. In real field production, such temperature changes are often experienced (normal temperature–chilling–restore). In this study, we analyzed rice plants’ physiological response and transcriptional regulation during three continuous and short periods of normal temperature–chilling–recovery to normal temperature. This study confirms the mechanism of hypothermic sensitivity reversal in Dular seedlings and provides knowledge on chilling tolerance mechanisms of rice.

2. Results

2.1. Changes in Chlorophyll and Proline and Non-Structure Carbohydrate (NSC) Contents of Rice Leaves

The chlorophyll and proline contents of the Dular were significantly different between normal-temperature treatment (CT), chilling treatment (TR), and recovery to normal-temperature treatment (RC) (Figure 1A). The chlorophyll a+b and carotenoid content of the Dular were significantly reduced under TR (0.28 mg·g⁻¹ Fresh Weight (FW) and 0.01 mg·g⁻¹ FW) as compared with CT (0.83 mg·g⁻¹ FW and 0.07 mg·g⁻¹ FW). On the one hand, the chlorophyll a+b and carotenoid content in RC (0.37 and 0.03 mg·g⁻¹ FW) was significantly higher than that in TR, although it was lower than in CT (Table S1). On the other hand, the proline content in RC (38.08 μg·g⁻¹·FW) was significantly lower than that in TR (55.60 μg·g⁻¹·FW), although higher than CT (29.13 μg·g⁻¹·FW) (Figure 1B). The NSC, including the soluble sugar content (SSC) and starch content (SC) in RC (8.44 mg·g⁻¹ Dry Weight (DW) and 35.69 mg·g⁻¹ DW) were significantly higher than TR (4.65 mg·g⁻¹ DW and 27.94 mg·g⁻¹ DW), although lower than CT (17.86 mg·g⁻¹ DW and 60.40 mg·g⁻¹ DW). These findings indicated that damage occurred in the leaves under TR and recovered quickly under RC, though they differed from CT.

2.2. Changes in K⁺ and Na⁺ Flux Rate of Rice Roots

The root meristem region of rice was measured in this study. There were significant differences in the ion flow rate of rice roots under three different treatments. K⁺ efflux under TR (31.37 pmol·cm⁻²·s⁻¹) was significantly higher than CT (6.11 pmol·cm⁻²·s⁻¹) and RC (6.43 pmol·cm⁻²·s⁻¹) (Figure 2A). There was no significant difference between RC and CT. The flux rate of Na⁺ revealed an opposite trend. Roots absorbed a small amount of Na⁺ under TR (−5.64 pmol·cm⁻²·s⁻¹), whereas CT (24.34 pmol·cm⁻²·s⁻¹) and RC (140.28 pmol·cm⁻²·s⁻¹) excreted a large amount of Na⁺ externally (Figure 2B). The efflux rate of RC under RC was significantly higher than that of CT. These results indicated that the ion balance in Dular was affected by chilling stress, resulting in the efflux of K⁺ and the absorption of Na⁺, thereby disrupting the K⁺/Na⁺ balance. When the temperature was restored, plants could effectively inhibit K⁺ efflux and excrete Na⁺ in large quantities.

2.3. Gene Identification

A total of 21,970 genes were identified in all CT, TR, and RC samples (Figure 3). Among these, 20,096, 20,177, and 19,209 genes were identified in CT, TR, and RC samples, respectively (Table S2). Furthermore, 839, 733 and 444 genes were identified in CT and TR, TR and RC, and CT and RC, respectively. A total of 17,748 genes were synchronously expressed across all three treatments (CT, TR, and RC).

2.4. Correlation between Different Temperature Treatments in Dular Samples

The principal component analysis (PCA) (Figure 4) of the three treatments depicted the high reproducibility in transcriptomic data. The difference between TR and CT samples was significant, indicating that the transcriptomic of Dular significantly changed under chilling, and the plants were subjected to intense stress. However, the differences between RC and CT samples were reduced than those between TR and CT samples, suggesting that the increase in temperature effectively mitigated the damage of chilling on Dular.

2.5. Transcriptomic Response to Chilling and Normal Temperature Recovery in the Leaves of Dular

A total of 11,732 genes were differentially expressed by expression ratio >2.000 or <0.050 during chilling (TR vs. CT (TR/CT)) and normal-temperature recovery (RC vs. TR (RC/TR)). Among those, 9143 genes were differentially expressed in TR/CT, and 6091 genes were differentially expressed in RC/TR. A total of 3502 genes were differentially expressed in both TR/CT and RC/TR (Figure S1, Table S3). Furthermore, 5524 differentially expressed genes (DEGs) were upregulated and 3619 DEGs were downregulated in TR/CT. Moreover, 2789 DEGs were upregulated and 3302 were downregulated in RC/TR (Figure S2).

2.6. KEGG Pathway Enrichment under Chilling and Subsequent Recovery

KEGG pathway enrichment analysis was performed using an R script for genes in gene concentration. The KEGG pathway was significantly enriched when the p value was < 0.05. In addition, KEGG enrichment analysis was carried out for differentially expressed genes from normal temperature to chilling (TR/CT) (Figure 5A) and their subsequent recovery to normal temperature (RC/TR) (Figure 5B). The results demonstrated that ten metabolic pathways were significantly enriched under chilling and subsequent recovery, including alpha-Linolenic acid metabolism, carbon fixation in photosynthetic organisms, diterpenoid biosynthesis, glyoxylate, and dicarboxylate metabolism, linoleic acid metabolism, phenylpropanoid biosynthesis, photosynthesis-antenna proteins, ribosome, ribosome biogenesis in eukaryotes, and starch and sucrose metabolism (Figure 6).

2.7. DEGs Enriched for KEGG under Chilling and Subsequent Recovery

Enrichment analysis demonstrated that DEGs were mainly enriched in response to the ribosome pathway (map03010), photosynthesis-antenna proteins (map00196), phenylpropanoid biosynthesis (map00940) and linoleic acid metabolism (map00591). There were 103 DEGs (58 upregulated and 45 downregulated) and 144 (71 upregulated and 73 downregulated) in TR/CT and RC/TR (Table S4). There were 5 DEGs (3 upregulated and 2 downregulated) and 46 DEGs (36 upregulated and 10 downregulated) responding to only TR/CT or RC/TR, respectively, whereas there were 97 DEGs in both TR/CT and RC/TR. Among these 97 DEGs, 11 DEGs were upregulated in TR/CT and RC/TR, whereas 19 DEGs were downregulated both in TR/CT and RC/TR. The heatmap illustrates the common genes (Figure 7A) and total genes (Figure S3) under three different temperature treatments. In addition, 44 DEGs were upregulated in RT/CT and downregulated in RC/TR, and 23 DEGs were downregulated in TR/CT and upregulated in RC/TR (Figure 7B). Moreover, 148 DEGs demonstrated 156 subcellular localization results, which indicated that some proteins might locate in multiple organelles. There were 46 (29.30%), 16 (10.19%), 16 (10.19%), 5 (3.18%), 11 (7.01%), 4 (2.55%), 16 (10.19%), 8 (5.73%) and 34 (21.66%) located in the chloroplast, cytoplasm, cytosol, endoplasmic reticulum, membrane, mitochondrial, nucleus, unknown and other regions, respectively (Figure 7C). These results indicated that many chloroplast genes are related to TR/CT and RC/TR.

2.8. DEGs Analysis of Common Responsive Genes between TR/CT and RC/TR

A total of 98 DEGs demonstrated a specific response in TR/CT and RC/TR. There were 44 genes upregulated in TR/CT and downregulated in RC/TR. Twelve of these genes were categorized as being related to phenylpropanoid biosynthesis, including seven peroxidase, two beta-glucosidase, cinnamoyl-CoA reductase 1, cytochrome P450, cinnamyl alcohol dehydrogenase 7. Five genes were categorized as being related to ribosomal, including 60S ribosomal protein L23, 60S ribosomal protein L32-1, syntaxin-61, 60S ribosomal protein L31, and 40S ribosomal protein S23. Moreover, genes related to linoleic acid metabolism included one phospholipase A2 homolog 3.

There were 24 genes downregulated in TR/CT, although they were upregulated in RC/TR. Genes related to linoleic acid metabolism included cytochrome P450 71D8, ent-isokaurene C2-hydroxylase, cytochrome P450 71A1, 3 lipoxygenase (2 lipoxygenases 2.3 and 1 lipoxygenase), linoleate 9S-lipoxygenase 6, and triacylglycerol lipase SDP1. Moreover, genes related to phenylpropanoid biosynthesis included one 3-aminomethylindole N-methyltransferase and one 4-coumarate-CoA ligase-like 2. One nucleolar protein 5-2 was associated with the ribosome.

Only 11 genes were upregulated in both TR/CT and RC/TR. Genes related to linoleic acid metabolism included 1 cytochrome P450 71D7. In addition, the transferase family protein was related to phenylpropanoid biosynthesis. In addition, 19 genes were downregulated in both TR/CT and RC/TR. Genes related to ribosomes included one 50S ribosomal protein L3-2 and one 50S ribosomal protein L25.

2.9. Analysis of Genes Only Responsive in TR/CT or RC/TR

Only five genes were differentially expressed in TR/CT; three were upregulated, and two were downregulated. Chlorophyll a–b binding protein CP24 upregulated was associated with photosynthesis—antenna proteins, and one ent-isokaurene C2-hydroxylase was downregulated due to linoleic acid metabolism.

In total, 36 genes which were only upregulated in RC/TR. Genes related to linoleic acid metabolism included 3 lipoxygenases and 1 ent-cassadiene hydroxylase. One beta-glucosidase 15 was regulated due to phenylpropanoid biosynthesis. Genes related to photosynthesis-antenna proteins included 10 chlorophyll a-b binding proteins.

Ten proteins were downregulated only in RC/TR. Genes related to the ribosome included 60S ribosomal protein L37-1, 60S ribosomal protein L9, RP-L32, ubiquitin-60S ribosomal protein L40-2, 40S ribosomal protein S10-2, and 40S ribosomal protein S13-2.

2.10. The Networking Analysis of Genes Responsive to Chilling and Subsequent Recovery

The Dular leaves were enhanced in response to exogenous chilling stimuli, as evidenced by the downregulation of the pathway of linoleic acid metabolism and the upregulation of the pathways of phenylpropanoid biosynthesis and ribosome. When plants recovered from chilling to normal temperature, leaves changed from tolerance to chilling to growth and development, the linoleic acid metabolism and photosynthesis pathway were upregulated, and the pathway of phenylpropanoid biosynthesis and ribosome were downregulated (Figure 8). qRT-PCR results also confirmed that the gene expression in linoleic acid metabolism and ribosome synthesis were consistent with the transcriptome results (Figure 9). The path of leaves’ responses to chilling and subsequent recovery is displayed in Figure 10.

3. Discussion

3.1. Leaves Responses to Temperature Change

Many plants suffer from physiological dysfunction after being subjected to chilling stress, adversely affecting their growth and development [36,37]. Although plants may acquire chilling tolerance (i.e., chilling adaptation) in their initial response to low but nonfreezing temperatures, continuous chilling exposure is considered catastrophic. Chloroplast plays an important role in rice development. Under chilling, rice leaves revealed symptoms of chilling injury and chlorophyll synthesis was blocked, leading to decreased photosynthesis [38,39,40]. Photosynthesis is an essential factor for rice yield and the most obvious physiological processes in rice under chilling [41]. Compared with chilling sensitive rice, chilling tolerant rice has excellent characteristics of tolerance and adaptation to cold damage caused by low temperature. According to our previous study wild-type Dular rice was sensitive to low temperatures when the ambient temperature was below 15 °C, and the leaves of this rice turned white in a short period [34,35]. Hence, Dular was used for all experiments in this study. This study found that the chlorophyll content of Dular leaves decreased significantly under chilling, whereas it recovered quickly during normal temperature, indicating that the Dular plants were vulnerable regarding chilling tolerance and susceptibility to environmental temperature. Plants produce many reactive oxygen species (ROS) under stress [42]. ROS are important signaling molecules in plant growth, development, and response to various stresses [43]. It has been reported that proline protects enzymes from denaturation and helps stabilize the machinery of protein synthesis, regulates the cytosolic acidity, enhances water-binding capacity, and acts as a reservoir of carbon and nitrogen sources [44]. Under stress, proline biosynthesis was activated in plants, and the proline degradation was inhibited, resulting in a large amount of proline accumulation in plants. Meanwhile, these proline could be decomposed and utilized quickly after the stress is relieved [45]. In this study, the proline content in TR was significantly higher than that in CT, which protected the rice system leaves, stabilized the cell membrane’s function, effectively eliminated ROS damage to the plant body, and improved the chilling tolerance of the plants. The proline content in RC was significantly higher than in TR, indicating a reduced chilling influence on the leaves. Large amounts of proline are converted to carbon and nitrogen sources. Meanwhile, the NSC content was reduced in the TR and rapid increase in RC. This result indicated that rice increased its tolerance under chilling stress by consuming a lot of NSC, and when the stress was relieved, rice recovered rapidly and synthesized a large amount of NSC for the body’s nutritional needs.

Maintaining K⁺/Na⁺ balance is vital for plants to adapt to stressful environments [46,47]. Under stress, maintaining the balance of K⁺/Na⁺ in the cytoplasm mainly includes two aspects: avoiding the excessive accumulation of Na⁺ in the cytoplasm by promoting Na⁺ efflux; limiting the loss of K⁺ [48,49]. In this study, K⁺/Na⁺ in rice roots was balanced under CT, while under TR, a large amount of Na⁺ was absorbed internally, and K⁺ was excreted externally, resulting in disruption of intracellular ion balance and impaired cell membrane function. Under RC, plants effectively inhibit K⁺ efflux and excrete Na⁺ in large quantities, indicating that RC effectively maintains the intracellular ion balance and alleviate the damage caused by low temperature stress.

These results indicated that rice sensed and transmitted signals, stimulating its chilling tolerance function through physiological response. However, the chilling sensitivity of rice is a complex biological process, that involves changes at both physiological and transcriptional regulation levels. Gene expression drives a series of cellular activities and regulating of these requires the interaction of transcription factors (proteins) and transcription factor-binding sites (DNA elements). Transcriptional regulation is not a simple independent process, but a highly interactive gene regulatory network composed of hundreds of transcription factors, target sequences, and co-regulators.

3.2. DEGs Involved in the Pathway of Photosynthesis

This study found that the photosynthetic pathway in the Dular was significantly improved during the transition from chilling to normal-temperature recovery and was enriched in all metabolic pathways, indicating that the activation of the photosynthetic pathway played an important role in the removal of chilling sensitivity in Dular. Photosynthesis is the basis of life on earth. The leaves of higher plants capture light energy mainly through chlorophyll a/b-binding protein, and most of the chlorophyll is located in PSI and PSII natural pigment-protein complexes [50]. Light-trapping regulation is necessary for photosynthesis to maximize the use of light energy under low light exposure and avoid damage to photosynthetic apparatus under high light exposure [51,52]. Chlorophyll a/b-binding proteins play an important role in regulating plant light capture. Sun et al. found that chlorophyll a/b-binding proteins in winter wheat leaves were downregulated under chilling [53], thus affecting photosynthesis. Wang et al. found that under chilling, chlorophyll a/b-binding proteins could induce the upregulated expression of photosynthesis-related proteins, enzymes, and development-related proteins to enhance plant tolerance [54]. All chlorophyll a/b-binding proteins in Dular’s leaves were upregulated at normal temperature, and photosynthesis increased to produce more photosynthetic products for nutritional needs. The chlorophyll content also supported this idea. This result indicated that the temperature recovery could enhance photosynthesis and alleviate the effects of chilling, further confirming the transcriptomic results.

3.3. DEGs Involved in the Pathway of Linoleic Acid Metabolism

Linoleic acid acts as a structural component to maintain the transdermal water barrier of the epidermis membrane fluidity. In addition, when released from membrane phospholipids, it can be enzymatically oxidized to a various derivatives involved in cell signaling [55]. Cytochrome P450 is a metabolic enzyme widely existing in aerobic organisms, which catalyzes and regulates the metabolism and transformation of endogenous substances or the activation and degradation of exogenous compounds in life [56,57]. Lipoxygenases catalyze the oxidation of polyunsaturated fatty acids to produce unsaturated fatty acid and hydroperoxides; this process is known as the hydroperoxidation of lipids [58]. Lipids, primarily stored as triacylglycerols (TAGs) in seed, are the most energy-dense and a common form of seed reserves to support seedling establishment [59]. TAG lipases first hydrolyze TAG stored in seeds to release free fatty acids (FAs) that are transported into peroxisomes and oxidized to acetyl-CoA via multi-ple rounds of β-oxidation cycles [60]. In this study, the expression levels of cytochrome (P450 71D8 and P450 71A1), lipoxygenase, linoleate 9S-lipoxygenase 6 and triacylglycerol lipase SDP1 were downregulated in TR/CT and upregulated in RC/TR, indicating that the plant was highly stressed. The material and energy circulation were obstructed under chilling, whereas the stress was relieved when the temperature subsequent reverted to normal. The linoleic acid metabolism upregulated in RC/TR is beneficial for plants to alleviate the damage caused by chilling.

3.4. DEGs Involved in the Pathway of Phenylpropanoid Biosynthesis

Phenylpropanoid biosynthesis plays a vital role in plant stress tolerance [61], including phenylpropanoid, coumarin, lignin, and flavonoids. Peroxidase is a redox enzyme produced by plants, and it can catalyze many reactions and has the dual function of eliminating the toxicity of hydrogen peroxide and phenols, amines, aldehydes, and benzene [62]. Peroxidase is widely found in plants and is involved in hormone metabolism [63], reactive oxygen signaling molecule metabolism [64], and plant tolerance [65]. β-glucosidase is an important hydrolytic enzyme [66], and it is involved in glucose metabolism and plays an important role in maintaining the normal physiological function of an organism. One study also found that β-glucosidase had biological significance in plants’ virulence, defense, and secondary metabolisms [67]. This study found that phenylpropyl biosynthesis was significantly enhanced in Dular under chilling, in which β-glucosidase and peroxidase were in the highest proportion, and the expression was upregulated overall. The results demonstrated that the tolerance of Dular to chilling was enhanced. The phenylpropane biosynthesis pathway in the Dular was downregulated during the transition from chilling to normal-temperature recovery. This was the exact opposite during the transition from normal-temperature recovery to chilling treatment in Dular. The results revealed that the phenylpropyl biosynthesis pathway was upregulated in Dular in the transition from normal temperature to chilling due to the effect of chilling on the stress tolerance of rice plants. However, when plants recovered from chilling to normal temperature, the phenylpropanoid biosynthesis pathway in the organism was downregulated, and the plant activity changed from resisting chilling to growth and development.

3.5. DEGs Involved in the Pathway of Ribosome

Ribosome synthesis is a key to plant growth and environmental adaptation [68]. For example, Moin et al. found that the ribosomal protein L6 enhanced the salt tolerance of rice [69]. Ribosomes are the oldest and most sophisticated organelles, and their structure and composition remain highly conserved from prokaryotic to eukaryotic. Ribosomes are a vital part of cell replication, responsible for protein synthesis, and play an important role in controlling cell growth, division, and development [70]. A large number of studies have demonstrated that in addition to forming ribosomes and participating in protein biosynthesis [71], ribosomal proteins are involved in chloroplast development [72], translation regulation [73], and other processes. This study further analyzed the downregulated ribosome synthesis pathway with significant enrichment. Ribosome biogenesis involves the transcription of ribosomal DNA (rDNA), precursor-rRNA (pre-rRNA) processing, RNA modifications, and assembly factors [74,75]. The ribosome synthesis in Dular plants under chilling was upregulated, enhancing the synthesis of other proteins and reducing rice damage.

4. Materials and Methods

4.1. Plant Growth and Treatment

In our previous study, wild-type Dular rice was sensitive to low temperatures when the ambient temperature was below 15 °C, and the leaves of this rice turned white in a short period of time [34,35]. Hence, Dular was used for all experiments in this study.

The plump rice seeds were selected, disinfected with 25% sodium hypochlorite solution for 30 min, washed with deionized water, and soaked for 24 h. Then, the seeds were put into an incubator for 24h. After germination, the seeds were evenly planted on the trays filled with matric soil (PINDSTRUP, http://www.pindstrup.cn/, accessed on 15 August 2022), and the trays were placed in the incubator (Yiheng, Shanghai, China). The culture conditions were as follows: 28 °C/day for 14 h, the light intensity of 20,000 Lux, dark culture at 22 °C/night for 10 h, and humidity of 85%. When the rice grew to the stage with 1.5 leaves, three continuous and short periods of treatments were carried out:

Normal-temperature treatment (CT): The culture conditions were 28 °C/day for 14 h, light intensity of 20,000 Lux, and 22 °C/night for 10 h. The treatment period was 24 h.

Chilling treatment (TR): Rice plants of CT were transferred to 15 °C/day for 14 h, with a light intensity of 20,000 Lux, and 10 °C/night for 10 h. The treatment period was 24 h.

Return to normal-temperature treatment (RC): Rice plants of TR were transferred to 28 °C/day for 14 h, with a light intensity of 20,000 Lux, and 22 °C/night for 10 h. The treatment period was 24 h.

The treatment sequence is as follows: normal temperature treatment for 24 h (sampling), then transferred to chilling treatment for 24 h (sampling), and finally transferred to normal temperature treatment for 24 h (sampling). The entire treatment took three days.

4.2. Determination of Chlorophyll, Proline and NSC Content

4.2.1. Determination of Chlorophyll Content

The fresh leaves (0.2 g, three biological repetitions) were grounded with 96% ethanol to a final volume of 25 mL, as proposed by Sadia et al. [35] and Wang [76]. The absorbance value of the extract at the wavelengths of 665 nm, 649 nm and 470 nm was determined using an Infinite M200 PRO NanoQuant (Tecan, Männedorf, Switzerland) microplate reader. As a blank control, 96% ethanol was used. Finally, the content of chlorophyll in plant tissues was calculated according to the following formula: Chlorophyll content (mg/g) = [concentration of chlorophyll × volume of extraction liquid × dilution ratio]/dry weight of the sample.

4.2.2. Determination of Proline Content

The content of the proline was measured using the kit according to the product manual (Cominbio, www.cominbio.com, accessed on 15 August 2022). The contents of proline were expressed as μg·g⁻¹·FW. The determination principle was as follows: proline was extracted with sulfosalicylic acid. After heat treatment, proline reacted with acidic ninhydrin solution to produce a red color. After adding toluene, the absorbance was measured at 520 nm and finally converted to proline content.

4.2.3. Determination of NSC Content

The content of NSC was measured by Wang’s method [76] as follows: the leaves were dried, crushed, and subjected to a 100-mesh sieve. Samples were weighted (0.2 g, three biological repetitions) and placed in glass test tubes. Next, 5 mL of deionized water was added to the test tubes and heated in a 100 °C boiling water bath for 30 min. The solution was filtered into a 50 mL volumetric flask, and the volume was fixed. Then, 0.5 mL extract was transferred from a 50 mL volumetric flask to a 20mL glass test tube, then added to 5mL anthrone-sulfuric acid solution (0.2% w/v anthrone dissolved in sulfuric acid) and heated in 100 °C boiling water for 15 min. After cooling, the absorbance at 630 nm was measured and finally converted into soluble sugar content.

After the soluble sugar extraction, the residue was transferred to a new 50 mL glass test tube. Then, 10 mL of deionized water was added and heated in the 100 °C boiling water for 15 min. Then 5 mL 9.2 M perchloric acid solution was added to a boiling water bath and heated for 15 min. Finally, the volume was fixed to 50 mL. Next, 0.5 mL extract was transferred from a 50 mL volumetric flask to a 20 mL glass test tube, then 5 mL anthrone-sulfuric acid solution (0.2% w/v anthrone dissolved in sulfuric acid) was added and heated in a 100 °C boiling water bath for 15 min. After cooling, the absorbance at 630 nm was measured and finally converted into starch content.

4.3. The Measurement of K⁺ and Na⁺ Fluxes

The non-invasive Micro-test Technology (NMT-YG-100, Younger USA LLC, Amherst, MA, USA) and V2.0 Software (Younger USA LLC, Amherst, MA, USA) was used to measure the flux rates of K⁺ and Na⁺ [77]. The technique was based on the nondestructive microelectrode technology, through the automatic control of a computer and precise motion control system, without touching the sample, the three-dimensional, real-time, and dynamic measurement, the sample in and out of various ion flow rates and their motion direction [78]. The roots were immobilized in a culture dish (diameter: 6 cm) and calibrated with a measuring solution for 10 min. The microelectrode (Younger USA LLC, Amherst, MA, USA) was vibrating between the root surface and measuring solution (two positions, i.e., 30 µm) and kept a perpendicular axis to the root (Figure 2A). The commercially available ionophore (K⁺ and Na⁺) (Younger USA LLC, Amherst, MA, USA) was filled in the microelectrode tip. The measuring solution for K⁺ consisted of 0.1mM CaCl₂, 0.1mM KCl and 0.3mM 2-(N-Morpholino) ethanesulfonic acid (MES), pH 6.0. The measuring solutions for Na⁺ consisted of 0.1mM CaCl₂, 0.1mM KCl, 0.1mM NaCl and 0.3mM MES, pH 6.0. The positive value of flow velocity data refers to efflux, while the negative value refers to the influx.

4.4. RNA Extraction

Total RNA of Dular-CT, Dular-TR and Dular-RC was extracted using TRIzol (TransGen Biotech, Beijing, China). The concentration of RNA was determined by Nanodrop 2000 C (Thermo Fisher Scientific, Waltham, MA, USA). RNA was stored in the Ultra-low temperature Freezer DW-HL678 (Meiling, HeFei, China) at −80 °C.

4.5. Transcriptome Sequencing

Eukaryotic mRNA sequencing was based on Illumina Novaseq 6000 sequencing platform (Meiji, http://www.majorbio.com/, accessed on 15 August 2022), and all mRNA of samples were sequenced. An Illumina Truseq RNA Sample Prep Kit (Meiji, http://www.majorbio.com/, accessed on 15 August 2022) was used to construct the gene pool. The operation flow and instrument reagents were as follows: RNA integrity was detected by agarose gel electrophoresis. The gene pool was constructed based on ensuring the RNA quality of the samples (total content of RNA ≥ 1 μg, content ≥ 35 ng/μL, OD260/280 ≥ 1.8, OD260/230 ≥ 1.0). First, the poly(A) tail structure at the 3’ end of a eukaryotic mRNA was exploited, and the poly(A) mRNA with poly(A) tail was isolated by A-T base pairing using magnetic beads with Oligo (dT) (Meiji, http://www.majorbio.com/, accessed on 15 August 2022). The Illumina Novaseq 6000 platform was designed to sequence short sequences. The enriched mRNA is a complete RNA sequence with an average length of several kb, which requires random interruption. Then mRNA was fragmented by adding fragmentation buffer, and small fragments of 300 bp were screened through magnetic beads to reverse cDNA synthesis. Under reverse transcriptase, random hexamers were added to synthesize one-strand cDNA using mRNA as a template, and then two-strand synthesis was carried out to form a stable double-strand structure. The cDNA structure of the double-strand was the vicious end. In addition, End Repair Mix was added to make the flat end, and then an “A” base was added to the 3’ end to connect the Y-shaped connector. The target strip was then enriched and recovered by the library. After the TBS380 quantitative analysis, the corresponding data were mixed and sequenced on the machine.

4.6. Bioinformatics Analysis of Transcriptomic Data from Oryza sativa L.

Quality control was carried out on the original sequencing data using fastp (http://github.com/OpenGene/fastp, accessed on 15 August 2022) before analysis to obtain high-quality clean data to ensure the accuracy of subsequent analysis results. The genome of Oryza sativa L. ssp. Indica was set as a reference. The mapped data for subsequent transcript assembly and expression volume calculations were obtained using TopHat2 software (http://ccb.jhu.edu/software/tophat/index.shtml, accessed on 15 August 2022). Using StringTie software (http://ccb.jhu.edu/software/stringtie/, accessed on 15 August 2022) to concatenate mapped reads. NCBI_NR database (http://www.ncbi.nlm.nih.gov, accessed on 15 August 2022) and Swiss-PROt database (http://web.expasy.org/docs/swiss-prot_guideline.html, accessed on 15 August 2022) were used for functional annotation and statistics of transcripts. The unigene expression was calculated and normalized to Fragments Per Kilobase of transcript per Million mapped reads (FPKM) by using RSEM software (http://deweylab.github.io/RSEM/, accessed on 15 August 2022). Then Venn and Principal Component Analysis (PCA) were used for sample relationship analysis.

The edgeR software (http://bioconductor.org/packages/stats/bioc/edgeR/, accessed on 15 August 2022) was used for differential gene identification, with the default threshold of p-adjusted < 0.05 and |log2FC| ≥ 1. KEGG pathway (https://www.kegg.jp/, accessed on 15 August 2022) analysis of differential genes was performed using R software (The University of Auckland, Oakland, New Zealand). A p value of 0.05 was used as the threshold for the KEGG pathway enrichment analysis, and the metabolic pathways of the differentially expressed genes were determined.

MajorBio online platform (https://cloud.majorbio.com/page/tools/, accessed on 15 August 2022) was used to analyze the relationship between the different temperature treatments for Dular gene samples.

4.7. qRT-PCR Validation

RNA was reversely transcribed into cDNA using the cDNA Synthesis SuperMix (TransGen Biotech, Beijing, China). The CDS sequences of the target genes were retrieved from the Rice Genome Annotation Project Database (http://rice.uga.edu/index.shtml, accessed on 15 August 2022), and the upstream and downstream primers were designed based on the qRT-PCR primer database (https://biodb.swu.edu.cn/qprimerdb/blast, accessed on 15 August 2022). The qRT-PCR was carried out using the TransScript Green qRT-PCR Supermix kit (TransGen Biotech, Beijing, China) and the Eppendorf Realplex4 (Eppendorf, HAM, DE) instrument was used. The qRT-PCR conditions were as follows: predenaturation at 94 °C for 30 s, denaturation at 94 °C for 5 s, and extension at 60 °C for 30 s (45 cycles). Finally, the relative expression of each candidate mRNA in different samples was calculated by the 2^−ΔΔCt method. The qRT-PCR primers are listed in Table S5.

4.8. Statistical Analysis

All of the data in the paper are the result of three biological replicates. The differences among the treatments were calculated and statistically analyzed using a one-way analysis of variance (ANOVA) by the least-significant-difference multiple-range test (LSD, p < 0.05). The statistical package for OriginPro 8.0 (OriginLab, Northampton, Ma, USA), GraphPad Prism version 5.0 (GraphPad, San Diego, CA, USA), and the Data Processing System (DPS) version 7.05 (Zhejiang University, Hangzhou, China) was used for the statistical analysis.

5. Conclusions

Low temperature severely affects the growth and development of rice. The recovery of temperature after a short period of low temperature ensured the synthesis of chlorophyll and NSC in rice and maintained the ion balance of the plant. The upregulation of photosynthesis, linoleic acid metabolism, and the downregulation of ribosome and phenylpropanoid biosynthesis indicate that plants changed from tolerance to stress to growth and development after temperature recovery. These results suggest that restoring normal temperature treatment can effectively reverse injury caused by chilling. This study revealed the process and mechanism of rice sensing, transmitting chilling signals, and stimulating the function of chilling tolerance, providing a theoretical foundation for further understanding the mechanism of chilling tolerance in plants and accelerating the progress of molecular breeding of chilling tolerant crops.

Author Contributions

Z.L., Y.J., Z.Z. and W.L.: conceptualization; Z.L., Y.J., Z.Z. and W.L.: investigation; Z.L., Y.X., W.Z., W.W. and M.U.K.: formal analysis; Z.L., Z.Z., P.L. and M.U.K.: writing. All authors have read and agreed to the published version of the manuscript.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

All data generated during this study are included in this published article and its supplementary information files, and the raw data used or analyzed during the current study are available from the corresponding author on reasonable request.

Conflicts of Interest

The authors declare no conflict of interest.

Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Figures

Figure 1. (A) The content of chlorophyll a + b under three different treatments. (B) The content of proline under three different treatments. (C) The SSC under three different treatments. (D) The SC under three different treatments. Error bars are the standard error (±SE) of three biological replications. * p < 0.05; ** p < 0.01 by Student’s t-test. CT: Normal-temperature treatment. TR: Chilling treatment. RC: Recovery to normal-temperature treatment. FW: Fresh weight. DW: Dry weight. SSC: Soluble sugar content. SC: Starch content.

Figure 2. (A) The K+ flux rate of root under CT, TR and RC. (B) The Na+ flux rate of root under CT, TR and RC. Columns with different letters are significantly different (LSD test, p < 0.01). Error bars are the standard error (±SE) of three biological replications. CT: Normal-temperature treatment. TR: Chilling treatment. RC: Recovery to normal-temperature treatment. Positive flux rate indicates efflux. Negative flux rate indicates influx.

Figure 3. The Venn diagram between three treatments (CT, TR and RC). CT: Normal-temperature treatment. TR: Low-temperature treatment. RC: Recovery to normal-temperature treatment. The numbers in the figure represent the number of genes identified.

Figure 4. The principal component analysis (PCA) analysis between three samples. After dimensionality reduction analysis, the samples have relative coordinate points on the principal components. The distance of each sample point represents the sample distance, and the closer the distance is, the higher the similarity between samples. The horizontal axis represents the contribution degree of principal component 1 (PC1) in the two-dimensional graph to the distinguished samples, and the vertical axis represents the contribution degree of principal component 2 (PC2) in the two-dimensional graph to the distinguished samples. CT: Normal-temperature treatment. TR: Low-temperature treatment. RC: Recovery to normal-temperature treatment.

Figure 5. (A) The pathways were significantly enriched by KEGG analysis from normal temperature to chilling. (B) The pathways were significantly enriched by KEGG analysis from chilling recovery to normal temperature. The triangles, ellipses, pentagons and rectangles demonstrated the top enriched genes for KEGG in the ribosome, phenylpropanoid biosynthesis, photosynthesis-antenna proteins and linoleic acid metabolism, respectively.

Figure 6. Functional characterization of chilling and recovery-responsive genes under chilling and/or subsequent recovery. CT: Normal-temperature treatment. TR: Chilling treatment. RC: Recovery to normal-temperature treatment. TR/CT: From normal temperature to chilling. RC/TR: From chilling recovery to normal temperature. The red columns show the differential expressed genes under RC/TR. The black columns show differential expressed genes under TR/CT.

Figure 7. (A) Heatmap analysis of the common genes under three different temperature treatments. CT: Normal-temperature treatment. TR: Chilling treatment. RC: Recovery to normal-temperature treatment. (B) Venn diagram of differentially expressed genes that were upregulated or downregulated by chilling stress or recovery. TR/CT: From normal temperature to chilling. RC/TR: From chilling recovery to normal temperature. The green circle on the left shows that the genes are differentially expressed only in TR/CT. The intersection in the middle represents the genes are differentially expressed in both TR/CT and RC/TR. The yellow portion on the right shows the genes that are differentially expressed only in RC/TR. The “+” and “−” indicate upregulated and downregulated genes, respectively. The “+−” indicates genes with upregulated expression under TR/CT and RC/TR treatments. The “+−” indicates genes with upregulated expression under TR/CT and downregulated expression under RC/TR. The “−+” indicates genes with downregulated expression under TR/CT and upregulated expression under RC/TR. The “−−” indicates genes with downregulated expression under TR/CT and RC/TR. (C) Subcellular localization of 149 differentially expressed genes under chilling and/or subsequent recovery.

Figure 8. Heatmap analysis of the genes under three different temperature treatments. CT: Normal-temperature treatment. TR: Chilling treatment. RC: Recovery to normal-temperature treatment. CT1, CT2 and CT3 are the three biological replications of CT. TR1, TR2 and TR3 are the three biological replications of TR. RC1, RC2 and RC3 are the three biological replications of RC. The “+” and “−” indicate upregulated and downregulated pathway by chilling stress (TR/CT) or recovery (RC/TR), respectively.

Figure 9. The qRT-PCR analysis shows the genes expression of Lipoxygenase 2.2, Lipoxygenase, 40S ribosomal protein S23, 60S ribosomal protein L23, 60S ribosomal protein L31 and 60S ribosomal protein L32-1 under three different temperature treatments. Columns with different letters are significant difference (LSD test, p < 0.01). Error bars are the standard error (±SE) of three biological replications. CT: Normal-temperature treatment. TR: Chilling treatment. RC: Recovery to normal-temperature treatment.

Figure 10. The pathway of leaves response to environmental conditions. Red arrows represent upregulation; green arrow represent downregulation.

Supplementary Materials

The supporting information can be downloaded at https://www.mdpi.com/article/10.3390/ijms231810739/s1.

References

1. Zhang, Z.; Li, J.; Li, F.; Liu, H.; Yang, W.; Chong, K.; Xu, Y. OsMAPK3 Phosphorylates OsbHLH002/OsICE1 and Inhibits Its Ubiquitination to Activate OsTPP1 and Enhances Rice Chilling Tolerance. Dev. Cell; 2017; 43, pp. 731-743. [DOI: https://dx.doi.org/10.1016/j.devcel.2017.11.016]

2. Tartoura, K.A.H.; Youssef, S.A. Stimulation of ROS-scavenging systems in squash (Cucurbita pepo L.) plants by compost supplementation under normal and low temperature conditions. Sci. Hortic.; 2011; 130, pp. 862-868. [DOI: https://dx.doi.org/10.1016/j.scienta.2011.08.015]

3. Taghavi, T.; Aghaani, M. Effect of chilling duration on runner production and vegetative growth of ‘Pajaro’ strawberry. Acta Hortic.; 2017; 1156, pp. 505-508. [DOI: https://dx.doi.org/10.17660/ActaHortic.2017.1156.75]

4. Zhang, X.D.; Liu, T.; Zhu, S.H.; Wang, D.; Sun, S.; Xin, L. Short-term hypobaric treatment alleviates chilling injury by regulating membrane fatty acids metabolism in peach fruit. J. Food Biochem.; 2022; 46, e14113. [DOI: https://dx.doi.org/10.1111/jfbc.14113]

5. Mao, D.; Xin, Y.; Tan, Y.; Hu, X.; Bai, J.; Liu, Z.Y.; Yu, Y.; Li, L.; Peng, C.; Fan, T. et al. Natural variation in the HAN1 gene confers chilling tolerance in rice and allowed adaptation to a temperate climate. Proc. Natl. Acad. Sci. USA; 2019; 116, pp. 3494-3501. [DOI: https://dx.doi.org/10.1073/pnas.1819769116]

6. Ma, Y.; Dai, X.; Xu, Y.; Luo, W.; Zheng, X.; Zeng, D.; Pan, Y.; Lin, X.; Liu, H.; Zhang, D. et al. COLD1 confers chilling tolerance in rice. Cell; 2015; 160, pp. 1209-1921. [DOI: https://dx.doi.org/10.1016/j.cell.2015.01.046]

7. Ma, Q.; Dai, X.; Xu, Y.; Guo, J.; Liu, Y.; Chen, N.; Xiao, J.; Zhang, D.; Xu, Z.; Zhang, X. et al. 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, pp. 244-256. [DOI: https://dx.doi.org/10.1104/pp.108.133454]

8. Kuk, Y.; Shin, J.; Burgos, N.; Hwang, T.; Han, O.; Cho, B.; Guh, J. Antioxidative enzymes offer protection from chilling damage in rice plants. Crop Sci.; 2003; 43, pp. 2109-2117. [DOI: https://dx.doi.org/10.2135/cropsci2003.2109]

9. Browse, J.; Xin, Z. Temperature sensing and cold acclimation. Curr. Opin. Plant Biol.; 2001; 4, pp. 241-246. [DOI: https://dx.doi.org/10.1016/S1369-5266(00)00167-9]

10. Thomashow, F. Plant Cold Acclimation: Freezing tolerance genes and regulatory mechanisms. Annu. Rev. Plant Biol.; 1999; 50, pp. 571-599. [DOI: https://dx.doi.org/10.1146/annurev.arplant.50.1.571]

11. Hughes, M.A.; Alison, D.M. The molecular biology of plant acclimation to low temperature. J. Exp. Bot.; 1996; 47, pp. 291-305. [DOI: https://dx.doi.org/10.1093/jxb/47.3.291]

12. Bohnert, J. Adaptations to Environmental Stresses. Plant Cell; 1995; 7, pp. 1099-1111. [DOI: https://dx.doi.org/10.2307/3870060]

13. Chen, T.; Murata, N. Enhancement of tolerance of abiotic stress by metabolic engineering of betaines and other compatible solutes. Curr. Opin. Plant Biol.; 2002; 5, pp. 250-257. [DOI: https://dx.doi.org/10.1016/S1369-5266(02)00255-8]

14. Prasad, T.K.; Anderson, M.D.; Martin, B.A.; Stewart, C.R. Evidence for chilling-induced oxidative stress in maize seedlings and a regulatory role for hydrogen peroxide. Plant Cell; 1994; 6, pp. 65-74. [DOI: https://dx.doi.org/10.2307/3869675]

15. Zhou, A.; Liu, E.; Li, H.; Li, Y.; Feng, S.; Gong, S.; Wang, J. PsCor413pm2, a plasma membrane-localized, cold-regulated protein from phlox subulata, confers low temperature tolerance in Arabidopsis. Int. J. Mol. Sci.; 2018; 19, 2579. [DOI: https://dx.doi.org/10.3390/ijms19092579]

16. Chinnusamy, V.; Zhu, J.; Zhu, J.-K. Cold stress regulation of gene expression in plants. Trends Plant Sci.; 2007; 12, pp. 444-451. [DOI: https://dx.doi.org/10.1016/j.tplants.2007.07.002]

17. Yamaguchi-Shinozaki, K.; Shinozaki, K. Transcriptional regulatory networks in cellular responses and tolerance to dehydration and cold stresses. Annu. Rev. Plant Biol.; 2006; 57, pp. 781-803. [DOI: https://dx.doi.org/10.1146/annurev.arplant.57.032905.105444]

18. Singh, K.B.; Foley, R.C.; Oñate-Sánchez, L. Transcription factors in plant defense and stress responses. Curr. Opin. Plant Biol.; 2002; 5, pp. 430-436. [DOI: https://dx.doi.org/10.1016/S1369-5266(02)00289-3]

19. Zhang, M.; Ye, J.; Xu, Q.; Feng, Y.; Yuan, X.; Yu, H.; Wang, Y.; Wei, X.; Yang, Y. Genome-wide association study of cold tolerance of Chinese indica rice varieties at the bud burst stage. Plant Cell Rep.; 2018; 37, pp. 529-539. [DOI: https://dx.doi.org/10.1007/s00299-017-2247-4]

20. Kim, C.; Vo, K.; Nguyen, C.; Jeong, D.; Lee, S.; Kumar, M.; Kim, S.; Park, S.; Kim, J.; Jeon, J. Functional analysis of a cold-responsive rice WRKY gene, OsWRKY71. Plant Biotechnol. Rep.; 2016; 10, pp. 13-23. [DOI: https://dx.doi.org/10.1007/s11816-015-0383-2]

21. Su, C.; Wang, Y.; Hsieh, T.; Lu, C.; Tseng, T.; Yu, S. A novel MYBS3-dependent pathway confers cold tolerance in rice. Plant physiol.; 2010; 153, pp. 145-158. [DOI: https://dx.doi.org/10.1104/pp.110.153015]

22. Yokotani, N.; Sato, Y.; Tanabe, S.; Chujo, T.; Shimizu, T.; Okada, K.; Yamane, H.; Shimono, M.; Sugano, S.; Takatsuji, H. et al. WRKY76 is a rice transcriptional repressor playing opposite roles in blast disease resistance and cold stress tolerance. J. Exp. Bot.; 2013; 64, pp. 5085-5097. [DOI: https://dx.doi.org/10.1093/jxb/ert298]

23. Liu, Y.; Xu, C.; Zhu, Y.; Zhang, L.; Chen, T.; Zhou, F.; Chen, H.; Lin, Y. The calcium-dependent kinase OsCPK24 functions in cold stress responses in rice. J. Integr. Plant Biol.; 2018; 60, pp. 173-188. [DOI: https://dx.doi.org/10.1111/jipb.12614]

24. Saijo, Y.; Kinoshita, N.; Ishiyama, K.; Hata, S.; Izui, K. A Ca²⁺-dependent protein kinase that endows rice plants with cold- and salt-stress tolerancefunctions in Vascular Bundles. Plant Cell Physiol.; 2001; 42, pp. 1228-1233. [DOI: https://dx.doi.org/10.1093/pcp/pce158] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/11726707]

25. Bulgari, R.; Cocetta, G.; Trivellini, A.; Vernieri, P.; Ferrante, A. Biostimulants and crop responses: A review. Biol. Agric. Hortic.; 2015; 31, pp. 1-17. [DOI: https://dx.doi.org/10.1080/01448765.2014.964649]

26. Sharma, H.; Fleming, C.; Selby, C.; Rao, J.; Martin, T. Plant biostimulants: A review on the processing of macroalgae and use of extracts for crop management to reduce abiotic and biotic stresses. J. Appl. Phycol.; 2014; 26, pp. 465-490. [DOI: https://dx.doi.org/10.1007/s10811-013-0101-9]

27. Light, M.; Daws, M.; Van Staden, J. Smoke-derived butenolide: Towards understanding its biological effects. S. Afr. J. Bot.; 2009; 75, pp. 1-7. [DOI: https://dx.doi.org/10.1016/j.sajb.2008.10.004]

28. Kulkarni, M.; Light, M.; Van Staden, J. Plant-derived smoke: Old technology with possibilities for economic applications in agriculture and horticulture. S. Afr. J. Bot.; 2011; 77, pp. 972-979. [DOI: https://dx.doi.org/10.1016/j.sajb.2011.08.006]

29. Nardi, S.; Pizzeghello, D.; Schiavon, M.; Ertani, A. Plant biostimulants: Physiological responses induced by protein hydrolyzed-based products and humic substances in plant metabolism. Sci. Agric.; 2016; 73, pp. 18-23. [DOI: https://dx.doi.org/10.1590/0103-9016-2015-0006]

30. Yakhin, O.; Lubyanov, A.; Yakhin, I.; Brown, P. Biostimulants in plant science: A global perspective. Front. Plant Sci.; 2017; 7, 2049. [DOI: https://dx.doi.org/10.3389/fpls.2016.02049]

31. Kannan, R.R.; Aderogba, M.A.; Ndhlala, A.R.; Stirk, W.A.; Van Staden, J. Acetylcholinesterase inhibitory activity of phlorotannins isolated from the brown alga, Ecklonia maxima (Osbeck) Papenfuss. Food Res. Int.; 2013; 54, pp. 1250-1254. [DOI: https://dx.doi.org/10.1016/j.foodres.2012.11.017]

32. Mi, J.; Li, G.; Huang, J.; Yu, H.; Zhou, F.; Zhang, Q.; Ouyang, Y.; Mou, T. Stacking S5-n and f5-n to overcome sterility in indica–japonica hybrid rice. Theor. Appl. Genet.; 2016; 129, pp. 563-575. [DOI: https://dx.doi.org/10.1007/s00122-015-2648-0] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/26704419]

33. Wang, J.; Liu, K.; Xu, C.; Li, H.; Zhang, Q. The high level of wide-compatibility of variety ‘Dular’ has a complex genetic basis. Theor. Appl. Genet.; 1998; 97, pp. 407-412. [DOI: https://dx.doi.org/10.1007/s001220050910]

34. Azeem, S.; Li, Z.; Zheng, H.; Lin, W.; Arafat, Y.; Zhang, Z.; Lin, X.; Lin, W. Quantitative proteomics study onLsi1in regulation of rice (Oryza sativa L.) cold resistance. Plant Growth Regul.; 2016; 78, pp. 307-323. [DOI: https://dx.doi.org/10.1007/s10725-015-0094-2]

35. Li, Z.; Feng, S.; Zhan, W.; Xu, L.; Fang, C.; Zhang, Z.; Lin, W. Lsi1 plays an active role in enhancing the chilling tolerance of rice roots. Plant Growth Regul.; 2020; 90, pp. 529-543. [DOI: https://dx.doi.org/10.1007/s10725-020-00577-5]

36. Zhang, Z.; Li, J.; Pan, Y.; Li, J.; Zhou, L.; Shi, H.; Zeng, Y.; Guo, H.; Yang, S.; Zheng, W. Natural variation in CTB4a enhances rice adaptation to cold habitats. Nat. Commun.; 2017; 8, 14788. [DOI: https://dx.doi.org/10.1038/ncomms14788]

37. Zhao, J.; Zhang, S.; Yang, T.; Zeng, Z.; Huang, Z.; Liu, Q.; Wang, X.; Leach, J.; Leung, H.; Liu, B. Global transcriptional profiling of a cold-tolerant rice variety under moderate cold stress reveals different cold stress response mechanisms. Physiol. Plant.; 2015; 154, pp. 381-394. [DOI: https://dx.doi.org/10.1111/ppl.12291]

38. Shimono, H.; Hasegawa, T.; Fujimura, S.; Iwama, K. Responses of leaf photosynthesis and plant water status in rice to low water temperature at different growth stages. Field Crop. Res.; 2004; 89, pp. 71-83. [DOI: https://dx.doi.org/10.1016/j.fcr.2004.01.025]

39. Ariizumi, T.; Kishitani, S.; Inatsugi, R.; Nishida, I.; Murata, N.; Toriyama, K. An increase in unsaturation of fatty acids in phosphatidylglycerol from leaves improves the rates of photosynthesis and growth at low temperatures in transgenic rice seedlings. Plant Cell Physiol.; 2002; 43, pp. 751-758. [DOI: https://dx.doi.org/10.1093/pcp/pcf087]

40. Yahmed, J.; Oliveira, T.; Novillo, P.; Quinones, A.; Forner, M.; Salvador, A.; Froelicher, Y.; Mimoun, M.; Talon, M.; Ollitrault, P. A simple, fast and inexpensive method to assess salt stress tolerance of aerial plant part: Investigations in the mandarin group. J. Plant Physiol.; 2016; 190, pp. 36-43. [DOI: https://dx.doi.org/10.1016/j.jplph.2015.10.008]

41. Jeong, S.; Choi, S.; Lee, D.; Ahn, S.; Hur, Y.; Chow, W.; Park, Y. Differential susceptibility of photosynthesis to light-chilling stress in rice (Oryza sativa L.) depends on the capacity for photochemical dissipation of light. Mol. Cells; 2002; 13, pp. 419-428. [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/12132582]

42. Zhang, K.; Tian, G.; Li, X.; Zhang, Z.; Liu, J.; Li, Y.; Xie, J.; Wang, P. ROS produced via BsRBOHD plays an important role in low temperature-induced anthocyanin biosynthesis in begonia semperflorens. Russ. J. Plant Physiol.; 2020; 67, pp. 250-258. [DOI: https://dx.doi.org/10.1134/S1021443720020181]

43. Apel, K.; Hirt, H. Reactive oxygen species: Metabolism, oxidative stress, and signal transduction. Annu. Rev. Plant Biol.; 2004; 55, pp. 373-399. [DOI: https://dx.doi.org/10.1146/annurev.arplant.55.031903.141701] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/15377225]

44. Liu, K.; Wang, L.; Xu, Y.; Chen, N.; Ma, Q.; Li, F.; Chong, K. Overexpression of OsCOIN, a putative cold inducible zinc finger protein, increased tolerance to chilling, salt and drought, and enhanced proline level in rice. Planta; 2007; 226, pp. 1007-1101. [DOI: https://dx.doi.org/10.1007/s00425-007-0548-5] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/17549515]

45. Trotel-Aziz, P.; Niogret, M.F.; Deleu, C.; Bouchereau, A.; Aziz, A.; Larher, F.R. The control of proline consumption by abscisic acid during osmotic stress recovery of canola leaf discs. Physiol. Plant.; 2003; 117, pp. 213-221. [DOI: https://dx.doi.org/10.1034/j.1399-3054.2003.1170208.x]

46. Sun, J.; Li, L.; Liu, M.; Wang, M.; Ding, M.; Deng, S.; Lu, C.; Zhou, X.; Shen, X.; Zheng, X. Hydrogen peroxide and nitric oxide mediate K⁺/Na⁺ homeostasis and antioxidant defense in NaCl-stressed callus cells of two contrasting poplars. Plant Cell Tissue Organ Cult.; 2010; 103, pp. 205-215. [DOI: https://dx.doi.org/10.1007/s11240-010-9768-7]

47. Sun, J.; Wang, M.; Ding, M.; Deng, S.; Liu, M.; Lu, C.; Zhou, X.; Shen, X.; Zheng, X.; Zhang, Z. H₂O₂ and cytosolic Ca²⁺ signals triggered by the PM H⁺-coupled transport system mediate K⁺/Na⁺ homeostasis in NaCl-stressed Populus euphratica cells. Plant Cell Environ.; 2010; 33, pp. 943-958. [DOI: https://dx.doi.org/10.1111/j.1365-3040.2010.02118.x]

48. Tyerman, S.; Skerrett, I. Root ion channels and salinity. Sci. Hortic.; 1998; 78, pp. 175-235. [DOI: https://dx.doi.org/10.1016/S0304-4238(98)00194-0]

49. Cuin, T.; Betts, S.; Chalmandrier, R.; Shabala, S. A root’s ability to retain K⁺ correlates with salt tolerance in wheat. J. Exp. Bot.; 2008; 59, pp. 2697-2706. [DOI: https://dx.doi.org/10.1093/jxb/ern128]

50. Green, B.; Pichersky, E.; Kloppstech, K. Chlorophyll a/b-binding proteins: An extended family. Trends Biochem. Sci.; 1991; 16, pp. 181-186. [DOI: https://dx.doi.org/10.1016/0968-0004(91)90072-4]

51. Ruban, A.V. Nonphotochemical chlorophyll fluorescence quenching: Mechanism and effectiveness in protecting plants from photodamage. Plant Physiol.; 2016; 170, pp. 1903-1916. [DOI: https://dx.doi.org/10.1104/pp.15.01935] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/26864015]

52. Belgio, E.; Kapitonova, E.; Chmeliov, J.; Duffy, C.; Ungerer, P.; Valkunas, L.; Ruban, A. Economic photoprotection in photosystem II that retains a complete light-harvesting system with slow energy traps. Nat. Commun.; 2014; 5, 4433. [DOI: https://dx.doi.org/10.1038/ncomms5433] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/25014663]

53. Sun, X.; Tan, Q.; Nie, Z.; Hu, C.; An, Y. Differential expression of proteins in response to molybdenum deficiency in winter wheatleaves under low-temperature stress. Plant Mol. Biol. Rep.; 2014; 32, pp. 1057-1069. [DOI: https://dx.doi.org/10.1007/s11105-014-0713-5]

54. Wang, J.; Zhang, H.; Huang, R. Expression analysis of low temperature responsive genes in eupatorium adenophorum spreng using cDNA-AFLP. Plant Mol. Biol. Rep.; 2007; 25, pp. 37-44. [DOI: https://dx.doi.org/10.1007/s11105-007-0005-4]

55. Whelan, J.; Fritsche, K. Linoleic acid. Adv. Nutr.; 2013; 4, pp. 311-312. [DOI: https://dx.doi.org/10.3945/an.113.003772] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/23674797]

56. Yang, X.; Wu, D.; Shi, J.; He, Y.; Pinot, F.; Grausem, B.; Yin, C.; Zhu, L.; Chen, M.; Luo, Z. Rice CYP703A3, a cytochrome P450 hydroxylase, is essential for development of anther cuticle and pollen exine. J. Integr. Plant Biol.; 2014; 56, pp. 979-994. [DOI: https://dx.doi.org/10.1111/jipb.12212] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/24798002]

57. Li, H.; Pinot, F.; Sauveplane, V.; Werck-Reichhart, D.; Diehl, P.; Schreiber, L.; Franke, R.; Zhang, P.; Chen, L.; Gao, Y. Cytochrome P450 family member CYP704B2 catalyzes the ω-hydroxylation of fatty acids and is required for anther cutin biosynthesis and pollen exine formation in rice. Plant Cell; 2010; 22, pp. 173-190. [DOI: https://dx.doi.org/10.1105/tpc.109.070326]

58. Wang, Z.; Benning, C. Chloroplast lipid synthesis and lipid trafficking through ER–plastid membrane contact sites. Biochem. Soc. Trans.; 2012; 40, pp. 457-463. [DOI: https://dx.doi.org/10.1042/BST20110752]

59. Yang, Y.; Benning, C. Functions of triacylglycerols during plant development and stress. Curr. Opin. Biotechnol.; 2018; 49, pp. 191-198. [DOI: https://dx.doi.org/10.1016/j.copbio.2017.09.003]

60. Penfield, S.; Pinfield-Wells, H.; Graham, I. Storage reserve mobilisation and seedling establishment in Arabidopsis. Arab. Book/Am. Soc. Plant Biol.; 2006; 4, e0100. [DOI: https://dx.doi.org/10.1199/tab.0100]

61. Liu, H.; Wang, Z.G.; Xu, W.H.; Zeng, J.; Li, L.X.; Li, S.L.; Gao, Z. Bacillus pumilus LZP02 promotes rice root growth by improving carbohydrate metabolism and phenylpropanoid biosynthesis. Mol. Plant-Microbe Interact.; 2020; 33, pp. 1222-1231. [DOI: https://dx.doi.org/10.1094/MPMI-04-20-0106-R] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/32597697]

62. Rattanawong, K.; Koiso, N.; Toda, E.; Kinoshita, A.; Tanaka, M.; Tsuji, H.; Okamoto, T. Regulatory functions of ROS dynamics via glutathione metabolism and glutathione peroxidase activity in developing rice zygote. Plant J.; 2021; 108, pp. 1097-1115. [DOI: https://dx.doi.org/10.1111/tpj.15497] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/34538012]

63. Susumu, H.; Katsutomo, S.; Hiroyuki, I.; Yuko, O.; Hirokazu, M. A large family of class III plant peroxidases. Plant Cell Physiol.; 2001; 42, pp. 462-468.

64. Almagro, L.; Gómez, R.; Belchi-Navarro, S.; Bru, R.; Ros, B.; Pedreño, M. Class III peroxidases in plant defence reactions. J. Exp. Bot.; 2009; 60, pp. 377-390. [DOI: https://dx.doi.org/10.1093/jxb/ern277]

65. Ma, Y.; Ren, X.; Liang, C. Exogenous Ca²⁺ enhances antioxidant defense in rice to simulated acid rain by regulating ascorbate peroxidase and glutathione reductase. Planta; 2021; 254, 41. [DOI: https://dx.doi.org/10.1007/s00425-021-03679-0]

66. Saloheimo, M.; Kuja-Panula, J.; Ylosmaki, E.; Ward, M.; Penttila, M. Enzymatic Properties and Intracellular Localization of the Novel Trichoderma reesei β-Glucosidase BGLII (Cel1A). Appl. Environ. Microbiol.; 2002; 68, pp. 4546-4553. [DOI: https://dx.doi.org/10.1128/AEM.68.9.4546-4553.2002]

67. Lee, L. Enrichment of bioactive isoflavones in soymilk fermented with β-glucosidase-producing lactic acid bacteria. Food Res. Int.; 2005; 38, pp. 551-559.

68. Hang, R.; Wang, Z.; Deng, X.; Liu, C.; Yan, B.; Yang, C.; Song, X.; Mo, B.; Cao, X. Ribosomal RNA biogenesis and its response to chilling stress in Oryza sativa. Plant Physiol.; 2018; 177, pp. 381-397. [DOI: https://dx.doi.org/10.1104/pp.17.01714]

69. Moin, M.; Saha, A.; Bakshi, A.; Madhav, M.S.; Kirti, P.B. Constitutive expression of Ribosomal Protein L6 modulates salt tolerance in rice transgenic plants. Gene; 2021; 789, 145670. [DOI: https://dx.doi.org/10.1016/j.gene.2021.145670]

70. Barakat, A.; Szick-Miranda, K.; Chang, I.-F.; Guyot, R.; Blanc, G.; Cooke, R.; Delseny, M.; Bailey-Serres, J. The organization of cytoplasmic ribosomal protein genes in the Arabidopsis genome. Plant Physiol.; 2001; 127, pp. 398-415. [DOI: https://dx.doi.org/10.1104/pp.010265]

71. Uzair, M.; Long, H.; Zafar, S.A.; Patil, S.B.; Chun, Y.; Li, L.; Fang, J.; Zhao, J.; Peng, L.; Yuan, S. Narrow leaf21, encoding ribosomal protein RPS3A, controls leaf development in rice. Plant Physiol.; 2021; 186, pp. 497-518. [DOI: https://dx.doi.org/10.1093/plphys/kiab075]

72. Zhou, K.; Zhang, C.; Xia, J.; Yun, P.; Wang, Y.; Ma, T.; Li, Z. Albino seedling lethality 4; Chloroplast 30S Ri-bosomal protein S1 is required for chloroplast ribosome biogenesis and early chloroplast development in rice. Rice; 2021; 14, 47. [DOI: https://dx.doi.org/10.1186/s12284-021-00491-y] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/34046768]

73. Mazahar, M.; Anusree, S.; Achala, B.; Divya, D.; Madhav, M.S.; Kirti, P.B. Study on transcriptional responses and identification of Ri-bosomal protein genes for potential resistance against brown planthopper and gall midge pests in rice. Curr. Genom.; 2021; 22, pp. 98-110.

74. Lin, J.; Lai, S.; Jia, R.; Xu, A.; Zhang, L.; Lu, J.; Ye, K. Structural basis for site-specific ribose methylation by box C/D RNA protein complexes. Nature; 2011; 469, pp. 559-563. [DOI: https://dx.doi.org/10.1038/nature09688] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/21270896]

75. Woolford, J.; Baserga, S. Ribosome biogenesis in the yeast Saccharomyces cerevisiae. Genetics; 2013; 195, pp. 643-681. [DOI: https://dx.doi.org/10.1534/genetics.113.153197]

76. Wang, X. Principles and Techniques of Plant Physiological and Biochemical Experiments; Higher Education Press: Beijing, China, 2006.

77. Zhao, R.; Sun, H.; Zhao, N.; Jing, X.; Shen, X.; Chen, S. The Arabidopsis Ca²⁺-dependent protein kinase CPK27 is required for plant response to salt-stress. Gene; 2015; 563, pp. 203-214. [DOI: https://dx.doi.org/10.1016/j.gene.2015.03.024]

78. Li, Q.; Li, B.; Kronzucker, H.; Shi, W. Root growth inhibition by NH₄⁺ in Arabidopsis is mediated by the root tip and is linked to NH₄⁺ efflux and GMPase activity. Plant Cell Environ.; 2010; 33, pp. 1529-1542. [DOI: https://dx.doi.org/10.1111/j.1365-3040.2010.02162.x]