Skip to main content

Transcriptomic, proteomic, and physiological comparative analyses of flooding mitigation of the damage induced by low-temperature stress in direct seeded early indica rice at the seedling stage



Low temperature (LT) often occurs at the seedling stage in the early rice-growing season, especially for direct seeded early-season indica rice, and using flooding irrigation can mitigate LT damage in rice seedlings. The molecular mechanism by which flooding mitigates the damage induced by LT stress has not been fully elucidated. Thus, LT stress at 8 °C, LT accompanied by flooding (LTF) and CK (control) treatments were established for 3 days to determine the transcriptomic, proteomic and physiological response in direct seeded rice seedlings at the seedling stage.


LT damaged chloroplasts, and thylakoid lamellae, and increased osmiophilic bodies and starch grains compared to CK, but LTF alleviated the damage to chloroplast structure caused by LT. The physiological characteristics of treated plants showed that compared with LT, LTF significantly increased the contents of rubisco, chlorophyll, PEPCK, ATP and GA3 but significantly decreased soluble protein, MDA and ABA contents. 4D-label-free quantitative proteomic profiling showed that photosynthesis-responsive proteins, such as phytochrome, as well as chlorophyll and the tricarboxylic acid cycle were significantly downregulated in LT/CK and LTF/CK comparison groups. However, compared with LT, phytochrome, chlorophyllide oxygenase activity and the glucan branching enzyme in LTF were significantly upregulated in rice leaves. Transcriptomic and proteomic studies identified 72,818 transcripts and 5639 proteins, and 4983 genes that were identified at both the transcriptome and proteome levels. Differentially expressed genes (DEGs) and differentially expressed proteins (DEPs) were significantly enriched in glycine, serine and threonine metabolism, biosynthesis of secondary metabolites, glycolysis/gluconeogenesis and metabolic pathways.


Through transcriptomic, proteomic and physiological analyses, we determined that a variety of metabolic pathway changes were induced by LT and LTF. GO and KEGG enrichment analyses demonstrated that DEGs and DEPs were associated with photosynthesis pathways, antioxidant enzymes and energy metabolism pathway-related proteins. Our study provided new insights for efforts to reduce the damage to direct seeded rice caused by low-temperature stress and provided a breeding target for low temperature flooding-resistant cultivars. Further analysis of translational regulation and metabolites may help to elucidate the molecular mechanisms by which flooding mitigates low-temperature stress in direct seeded early indica rice at the seedling stage.


Rice (Oryza sativa L.) is the staple food for more than half of the world’s population [1,2,3]. Due to the advances made since the Green Revolution of the 1960s, rice yield has increased considerably. As a typical thermophilic crop, the growth and development of rice are susceptible to changes in temperature, especially decreases in temperature [4, 5]. With the rise of global temperature and frequent occurrence of extreme weather, direct seeded rice is strongly affected by weather compared with traditional transplanting, especially rice seedlings emergence [6, 7]. After direct sowing, heavy rainfall and “cold spell in later spring” disaster make it easy for a large area of rotten seeds and rotten seedlings, resulting in uneven seedling emergence and poor population growth, ultimately reducing the production of direct seeded rice [8, 9]. In addition, due to global warming and changes in production habits, sowing dates have become earlier than before, which increases the probability that direct seeded rice will suffer from low temperature to some extent, especially in the seedling stage of direct seeded early rice, resulting in irreversible cold-tolerant growth of seedlings [10, 11]. In 2008 and 2010, due to the low temperature and severe cold stress in China, the emergence rate of direct seeded early rice decreased by 38–55%, and rice yield notably decreased [12]. At the same time, although rice is a water-loving crop, its growth and development are also affected by long-term flooding stress, resulting in anaerobic respiration of roots and leaves, which produces alcohol toxicity and reduces leaf photosynthesis [13]. Therefore, it is urgently important to perform in-depth research on the response mechanism to low-temperature stress of direct seeded early rice and prevention measures to alleviate this stress and improve the efficiency and stable yield of direct seeded rice.

In direct seeded rice production, when low temperatures occur, farmers often reduce the damage of low temperature to seedlings by irrigating a certain amount of shallow water to protect the seedlings [14]. Through this measure, the survival rate of seedling emergence can effectively increase, and the production risk of direct seeded rice can be reduced [15, 16]. To date, many studies have investigated rice cultivation, physiological traits, genetic mechanisms and other aspects of injuries induced by low-temperature stress [17,18,19] and flooding stress [20, 21]. The effects of flooding on the physiological and ecological characteristics of rice under low temperature have been reported in previous studies, but the conclusions reached by these studies were inconsistent [22, 23]. It has been reported that flooding at the seedling stage significantly increases plant height and internodes and accelerates the growth of germ sheaths, which can effectively reduce the dead seedling rate [24, 25]. The effect of different flooding depths on rice varies significantly. Moderate flooding can stimulate changes in physiological characteristics in plants, thereby promoting the growth of plant height and causing rice to exhibit better adaptability [26]. When the seedlings encounter low-temperature stress, moderate flooding could increase temperature for heat preservation effects, alleviate the accumulation of reactive oxygen species, intensify membrane lipid peroxidation under low temperature conditions, and slow the regulation of endogenous hormones in plants [27]. The direct damage to early rice seedlings caused by low temperature may thus be prevented [28]. At present, most studies on this topic have focused on agronomic traits or related physiological characteristics under different flooding layers [29, 30]. However, the molecular mechanism governing the mitigation effect of shallow flooding irrigation on low-temperature stress in direct seeded early indica rice seedlings has rarely been reported.

With the rapid development of biotechnology, an increasing number of studies on rice in response to different stresses have been analysed in depth by transcriptomic technologies [31, 32]. Transcriptome analysis is rapid and comprehensive and has been constructed and annotated to assist in the identification of differentially expressed genes (DEGs) in different plant populations [33]. However, the analysis of gene expression by measuring mRNA is limited, as mRNA is defined as indirect and temporary messages that transmit information. In contrast, proteins play a direct role in biological processes and are the basis of organisms [34]. Protein is the embodiment and executor of plant functions, which not only regulates plant stress tolerance by changing the catalytic activity of enzymes, but also acts as a transcription factor to regulate the expression of other genes [35,36,37]. Through the combination of the transcriptome and proteome, many differentially expressed proteins (DEPs) have been identified, and metabolic pathways have been found [38]. On this basis, many DEGs related to metabolic pathways have also been identified, providing a molecular mechanism for detecting responses to environmental stress.

At present, many studies have elucidated the mechanisms by which low-temperature stress or flooding stress affect rice seedlings from the aspects of proteomics and transcriptomics [39, 40]. However, the changes in transcriptomics and proteomics associated with low-temperature flooding have not been elucidated. The molecular mechanism of the mitigating effect, rather than superposition effect, of the hypoxic treatment caused by flooding under low temperature is a scientific problem that merits further study. This study combined transcriptomics and 4D-label-free quantitative proteomic analysis to explore the molecular mechanism by which flooding mitigates low-temperature stress on direct seeded early indica rice at the seedling stage. In this study, we identified genes and proteins that were obtained from Illumina-Hiseq and 4D-label-free searching for likely protein identification in Gene Ontology (GO), Kyoto Encyclopedia of Genes and Genomes (KEGG), Klustersof eukaryotic Orthologous Groups (KOG), Swissport and UniProt databases, respectively, and focused on the DEGs and DEPs involved in the flooding-mediated mitigation of low-temperature stress. The results of this study may help to guide the breeding and cultivation of low temperature-tolerant crop cultivars. This research also provides evidence for meteorological disaster mitigation and low temperature-induced damage prevention.


Transmission electron microscopic observation of chloroplast structure

In this study, transmission electron microscopy was employed to compare the differences in chloroplasts structural of early indica rice seedlings after 3 days between LT and LTF groups. The results showed that chloroplasts were regular boat-shaped or spindle-shaped and that thylakoid lamellae were clearly arranged close to the inner wall of cells in CK (Fig. 1c). Compared with CK, chloroplasts began to degrade, exhibiting distorted and loosely structured shapes in LTF (Fig. 1b); however, in LT, chloroplasts were severely degraded, thylakoid lamellae were seriously damaged, and osmiophilic bodies and starch grains increased gradually (Fig. 1a). The damage to chloroplast structure in LTF was less than that observed in LT. These results showed that chloroplasts were affected to some extent by low-temperature stress, and flooding could alleviate low-temperature damage to the chloroplast structure.

Fig. 1

Transmission electron microscope analysis of the top leaves at rice seedlings stage after low temperature and low temperature flooding stress. LT: low temperature, LTF: low temperature flooding, CK: control. Thy: thylakoid lamellae, OB: osmiophilic body, CP: chloroplast, CW: cell wall, MT: mitochondrion, SG: starch grain, NC: nucleus

Analysis of photosynthesis activity and endogenous hormone content

This study showed that the rubisco content of LT was significantly decreased by 26.97% (P < 0.05) compared to that of CK, but there was no significant difference between LTF and CK (Fig. 2a). The PEPCK activity, chlorophyll content and ATP content of LT and LTF were significantly decreased (P < 0.05) compared with CK, and those indexes were lower in LT than in LTF (Fig. 2b, c, d). Compared with CK, the GA3 content of pro-growth hormones in LT and LTF decreased significantly (P < 0.05) (Fig. 2e). In contrast, the ABA content of anti-growth hormones in LT and LTF increased significantly by 35.71 and 16.67% (P < 0.05), respectively (Fig. 2f). There were similar trends in GA3 and ABA between LT and LTF, which reached a significant level. These results indicated that flooding could improve photosynthetic activity and endogenous hormone content under low-temperature stress in direct seeded early indica rice at the seedling stage.

Fig. 2

Analysis photosynthate activity and endogenous hormone content. a rubisco content, b chlorophyll content, c PEPCK activity, d ATP content, e GA3 content. f ABA content. Error bars represent standard deviation (n = 3). Data are mean ± SD. The data were detected by Tukey’s honest significant difference (HSD), and different lowercase letters indicated significant differences at P < 0.05. LT: low temperature, LTF: low temperature flooding, CK: control

Analysis of soluble protein content, antioxidase and osmotic regulatory substances

This study showed that the contents of soluble protein and MDA in LT and LTF increased significantly compared with CK, and LT also significantly increased soluble protein and MDA compared to LTF (P < 0.05) (Fig. 3a, b). In addition, LT and LTF significantly increased the activities of SOD and POD (P < 0.05) compared with CK, and there were no significant differences between LT and LTF, although SOD and POD were higher in LT than in LTF (Fig. 3c, d).

Fig. 3

Analysis of soluble protein content, antioxidase and osmotic regulatory substances. a soluble protein content, b MDA content, c SOD activity, d POD activity. Error bars represent standard deviation (n = 3). Data are mean ± SD. The data were detected by Tukey’s honest significant difference (HSD), and different lowercase letters indicated significant differences at P < 0.05. LT: low temperature, LTF: low temperature flooding, CK: control

Identification of DEPs and DEGs

To evaluate the reliability of the data through proteomic analysis, the Pearson correlation coefficient was calculated for each of three samples, which indicated good reproducibility of the three biological replicates in each treatment (Fig. 4a). In addition, a total of 412,489 spectra were detected, 236,880 of which could be matched to peptides in the database, and 28,934 of the peptides were unique. In total, 5639 proteins could be identified, and 4518 proteins were experimentally quantified (Table 1).

Fig. 4

a Pearson correlation coefficient thermograph of protein quantification. A Pearson coefficient closer to − 1 indicates a negative correlation, a coefficient closer to 1 indicates a positive correlation and a coefficient closer to 0 indicates no correlation. b Summary of up-regulated and down-regulated of DEPs between the treatment groups (LT, LTF) and control group (CK). c Venn diagram the DEPs between the treatment groups (LT, LTF) and control group (CK). d Comparison of transcriptome and proteome identification. LT: low temperature, LTF: low temperature flooding, CK: control

Table 1 MS/MS spectrum database search analysis summary

Proteins with fold change (FC) values > 1.5 or (FC) values < 0.67 (P < 0.05) between the treatment (LT, LTF) and control groups (CK) were regarded as DEPs, and DEPs were hence considered low temperature- and low temperature flooding-responsive proteins at the seedling stage. There were 567 DEPs between LT and CK, 239 DEPs between LTF and CK, and 235 DEPs between LTF and LT. The number of upregulated and downregulated DEPs is shown in Fig. 4b, and the three groups had 16 DEPs in common (Fig. 4c).

In this study, 72,818 transcripts and 5639 proteins were identified by quantitative transcriptome and proteome studies. A total of 4983 genes were identified at both the transcriptome and proteome levels (Fig. 4d). The correlation coefficient between transcripts and proteins in the LT and CK treatment groups was 0.19, that in the LTF and CK treatment groups was 0.25 and that in the LT and LTF treatment groups was 0.22. This finding indicates that the correlation degree of samples in each treatment group is low, and these results are largely consistent with the expected results (Fig. 5).

Fig. 5

The transcript and its corresponding protein expression scatter diagram. LT: low temperature, LTF: low temperature flooding, CK: control

Gene functional description and GO analysis

To annotate the function of low-temperature flooding-responsive proteins, the protein IDs were searched in the NCBI database ( and/or the UniProt database ( For the DEPs between LT and CK, 197 upregulated proteins and 369 downregulated proteins exhibited annotated functions, and 1 downregulated protein remained uncharacterized (Additional file 1: Dataset S1). For the DEPs between LTF and CK, both 114 upregulated proteins and 125 downregulated proteins could be annotated with functions (Additional file 2: Dataset S2). For the DEPs between LTF and LT, both 154 upregulated proteins and 81 downregulated proteins showed annotated functions (Additional file 3: Dataset S3).

To determine the cellular component (CC), molecular function (MF) and biological process (BP) categories of GO for the low temperature- and low temperature flooding- responsive proteins, we searched their protein IDs from the GO database. GO analysis showed that the DEPs were involved in 14 subgroups of BP (Fig. 6a), ten subgroups of CC (Fig. 6b), and ten subgroups of MF (Fig. 6c) between LT and CK. The main biological process categories were metabolic process (30%), cellular process (25%), single-organism process (18%), response to stimulus (7%), localization (7%), biological regulation (5%) and cellular component organization or biogenesis (4%). The cellular component categories were cell (34%), organelle (25%), membrane (24%), and macromolecular complex (12%). The molecular function categories were binding (44%), catalytic activity (42%), transporter activity (5%), structural molecule activity (4%), and antioxidant activity (2%) (Additional file 4: Fig. S1).

Fig. 6

GO classification analysis of DEPs between LT, LTF and CK. a cellular component, b molecular function, c biological process. LT: low temperature, LTF: low temperature flooding, CK: control

GO analysis showed that the DEPs were associated with 13 subgroups of BP (Fig. 6a), nine subgroups of CC (Fig. 6b), and ten subgroups of MF (Fig. 6c) between LTF and CK. The biological process categories were metabolic process (30%), cellular process (26%), single-organism process (21%), and response to stimulus (7%), biological regulation (5%), cellular component organization or biogenesis (4%), and localization (4%). The cellular component categories were cell (36%), membrane (28%), organelle (26%), macromolecular complex (4%), and extracellular region (3%). The molecular function categories were catalytic activity (47%), binding (44%), structural molecule activity (2%), and transporter activity (2%) (Additional file 5: Fig. S2).

GO analysis showed that the DEPs were involved in 14 subgroups of BP (Fig. 6a), eight subgroups of CC (Fig. 6b), and nine subgroups of MF (Fig. 6c) between LTF and LT. The biological process categories were metabolic process (27%), cellular process (19%), single-organism process (17%), localization (10%), response to stimulus (7%), biological regulation (5%), developmental process (3%), multicellular organismal process (3%), cellular component organization or biogenesis (3%), and reproduction (3%). The cellular component categories were membrane (33%), cell (29%), organelle (23%), and macromolecular complex (12%). The molecular function categories were binding (45%), catalytic activity (41%), transporter activity (7%), and structural molecule activity (3%) (Additional file 6: Fig. S3).

Protein–protein interaction

The functional DEPs of all annotations were utilized to analyse protein interactions. This analysis demonstrated that most enzymatic proteins and proteins related to biosynthesis of secondary metabolites, monobactam biosynthesis, metabolic pathways, pentose phosphate pathway, fructose and mannose metabolism, glycolysis/gluconeogenesis, glycine, serine and threonine metabolism, arachidonic acid metabolism, biosynthesis of amino acids, phenylalanine, tyrosine and tryptophan biosynthesis and proteasome-related protein interactions were affected by LT and CK (Additional file 7: Fig. S4). Most enzymatic proteins and metabolic pathways, biosynthesis of secondary metabolites, carotenoid biosynthesis, ribosome biogenesis in eukaryotes, glycolysis/gluconeogenesis, glycine, serine and threonine metabolism photosynthesis and thiamine metabolism were observed for the interaction between LTF and CK (Additional file 8: Fig. S5). Most enzymatic proteins and photosynthesis-antenna proteins and photosynthesis-related protein interactions were affected by LTF and LT (Additional file 9: Fig. S6). In the LT/CK and LTF/CK comparison groups, the relevant DEPs in the metabolic pathway included A2X8P7, A2XLW5, A2XYG6 and B8AYU2, which indicated energy metabolism-related proteins. In the glycine, serine and threonine metabolic pathways, the relevant DEPs had A2YMZ1 and A2YCB9, which indicated photosynthesis-related proteins. In this study, the photosynthesis pathway and energy metabolism pathway were highly enriched under low temperature and low temperature flooding. This result showed that low temperature flooding played an important role in regulating the photosynthetic capacity of rice leaves. Consistent with our GO analysis findings, the majority of proteins were determined to be involved in photosynthesis and metabolic processes. We focused on proteins related to photosynthesis and metabolism at the proteomic level.

KEGG pathway analysis

All of the DEGs and DEPs were analysed for the KEGG over-representation of pathways to obtain functional insights into the differences between LT, LTF and CK treatments. The significantly (P < 0.01) enriched KEGG pathways are shown in Table 2. The KEGG pathways (ordered by rank) were monobactam biosynthesis, glycine, serine and threonine metabolism, biosynthesis of secondary metabolites, pentose phosphate pathway, biosynthesis of amino acids, metabolic pathways, arachidonic acid metabolism, glycolysis/gluconeogenesis, proteasome, phenylalanine, tyrosine and tryptophan biosynthesis, and fructose and mannose metabolism between LT and CK. The KEGG pathways (ordered by rank) were thiamine metabolism, ribosome biogenesis in eukaryotes, carotenoid biosynthesis, biosynthesis of secondary metabolites, metabolic pathways, glycine, serine and threonine metabolism, and glycolysis/gluconeogenesis between LTF and CK. The KEGG pathways (ordered by rank) were photosynthesis and photosynthesis-antenna proteins between LTF and LT. In the LT/CK and LTF/CK comparison groups, the co-enriched KEGG pathways were glycine, serine and threonine metabolism, biosynthesis of secondary metabolites, metabolic pathways, and glycolysis/gluconeogenesis.

Table 2 Enriched KEGG Pathways associated with DEGs and DEPs

Analysis of DEGs by qRT-PCR

To verify the transcriptome results, 18 related genes, including five upregulated genes and 13 downregulated genes, were analysed. The mRNA expression of A2XLW5, B8AYU2, A2X822, B8ASV8 and A2XYC2 was downregulated between LT and CK (Fig. 7a). The mRNA expression levels of A2YM28, A2XKN7, A2WRR8 and A2X8P7 were downregulated between LTF and CK, and B8AZB8, B8BJP8, and B8AS16 were significantly upregulated between LTF and CK (Fig. 7b). In LTF and LT, A2YMZ1, A2YCB9, A2YHC5 and A2YLE6 were determined to be downregulated, and B8B7M5 and A2YP23 were observed to be upregulated (Fig. 7c). These results confirmed that the transcriptome results reflected the relative expression of each gene, in which upregulated or downregulated genes in qRT-PCR were entirely consistent with transcriptome trends. Therefore, the transcriptome results were reliable.

Fig. 7

Comparison of the expressions of RNA-Seq and qRT-PCR results. a: LT/CK, b: LTF/CK, c: LTF/LT. The transcript levels of the selected genes were each normalized to that of the 18 s gene. The relative expression of qRT-PCR transcript was calculated according to the standard curve and normalized to 18 s gene. The calculated data (mean ± SD) of three individual (n = 3). LT: low temperature, LTF: low temperature flooding, CK: control


Effect of flooding on photosynthesis under low-temperature stress

The effects of LT stress on plants are multifaceted, such as growth retardation, physiological metabolic disorders, inhibition of photosynthetic efficiency, and cell membrane damage [41]. Photosynthesis is the most sensitive process for plants in response to LT stress, and chlorophyll is the catalyst for photosynthesis in plants [42, 43]. Our study showed that LT significantly decreased the chlorophyll content, and degraded chloroplasts and membrane lipids, resulting in numerous and large osmiophilic granules, since the expression of chlorophyll-binding proteins (A2YCB9, A2XAH0, A2XN43) induced by LT was significantly downregulated compared with CK (Additional file 1: Dataset S1). LT-induced chlorophyll-associated protein degradation may lead to inhibited photosynthesis that is directly linked with photosynthesis metabolism. This phenomenon could conceivably represent one of the main causes of rice susceptibility to LT stress [44], and our study obtained similar findings (Additional file 3: Dataset S3). Additionally, previous studies found that submergence caused O2 and CO2 deficiency by limiting ambient gas exchange, reduced the light available for photosynthesis and disrupted cellular energy production and ionic homeostasis [45]. However, this study indicated that 5–6 cm flooding had a high chlorophyll content under low-temperature stress, which was attributed to the high expression levels of chlorophyll-related proteins (A2XCH9 and B8B1R3) induced by LTF. In field production, the minimum water temperature is often 3 °C–4 °C higher than the minimum air temperature; therefore, moderate flooding could increase water temperature and alleviate the direct damage to early rice seedlings caused by low temperature, thereby reducing losses in rice production [26, 46].

Chlorophyll is also the main pigment of plant photosynthesis, and its biosynthesis is a series of enzymatic reactions [47]. Low-temperature stress might reduce the activities of key enzymes related to chlorophyll synthesis and affect the process of chlorophyll synthesis, while the effect of flooding on chlorophyll synthesis enzymes is influenced by the duration and depth of flooding [48, 49]. Rubisco and ATP synthases are key enzymes that catalyse the photosynthetic carbon cycle and are highly sensitive to the external environment [50, 51]. ATP synthase β subunit content was significantly reduced at 6 °C in tomato [52], rubisco activity is one of the major limiting factors of photosynthesis, and the photosynthetic capacity of leaves can be enhanced by increasing the rubisco content under low-temperature conditions [53, 54]. Our study had similar results on ATP content, and rubisco and PEPCK activities in rice seedlings, due to the lower expression of ATP synthase-associated proteins (A2XF96 and B8AC70), rubisco-associated protein (A2Z7X6) and pyruvate kinase-associated proteins (A2XEP1 and B8AYC1) (Additional file 1: Dataset S1). However, moderate flooding could promote the growth and development of rice seedlings under low-temperature stress and mitigate the damage caused by low temperature [55, 56]. The present results also demonstrated that LTF significantly increased ATP content, rubisco and PEPCK activities that were associated with high ATP-associated proteins (A2XF96 and B8AMY3), rubisco-associated protein (B8BIM0) and PEPCK-associated proteins expression (B8ACJ0 and B8ACJ0) (Additional file 3: Dataset S3). Flooding could maintain the preferable sheath growth point and root system under low temperature conditions, which improved the ambient temperature tolerance to low temperature of rice, and the increase in rubisco and PEPCK activities significantly improved photosynthetic efficiency within a certain range.

Effect of flooding on the antioxidant system under low-temperature stress

When rice was subjected to LT stress, the balance of intracellular oxygen metabolism was disturbed, resulting in the production of reactive oxygen species (ROS). The accumulation of ROS might cause oxidative stress, but rice can regulate it through antioxidant enzymes (such as POD and SOD) to resist the toxic effects of reactive oxygen free radicals under LT stress [57,58,59]. The accumulation of ROS scavenging enzymes has been determined to be one of the most important responses during LT stress [60]. Moderate flooding could further enhance the antioxidant capacity of plants, inhibit membrane peroxidation and reduce the aggravation of leaf peroxidation damage [61]. Our study found that under LT treatment, the key proteins of the antioxidant pathway, including ROS scavenging (A2XX54 and B8AKJ8), oxidoreductase (A2YAP5) and glutathione peroxidase (A2X822 and B8ASV8), were significantly upregulated (Additional file 1: Dataset S1). Meanwhile, superoxide dismutase activity (A2XAA0), oxidoreductase (A2YAP5) and glutathione peroxidase (B8ASV8) were significantly upregulated between LTF and CK (Additional file 2: Dataset S2), which was consistent with the increase in the activities of antioxidant enzymes, such as POD and SOD, in LT and LTF. The results of this study were in keeping with the findings of a previous study [45]. LT stress upregulated ROS scavenging, redox adjustment, cryoprotection and defence detoxification-related proteins. To survive under LT stress, plants induce the antioxidant enzyme system to remove reactive oxygen groups to prevent these groups from destroying the cell and to enhance their stress tolerance [62,63,64]. However, because the temperature was too low and the damage to rice seedlings was serious, excessive ROS attacked rice, resulting in cell structure damage and metabolic disorder, which might be one of the crucial reasons for the damage to chloroplast structure. In this study, the expression of peroxidase protein (B8ARU3) was significantly downregulated in the LTF and LT groups compared to the LTF group. Flooding plays a role in protecting seedlings from the direct stress of LT because shallow flooding could improve the antioxidant enzyme protection system of direct seeded rice seedlings under LT. The cold resistance of rice may be the result of the coordinated action of antioxidant enzymes and various resistance factors [65].

MDA is the final product of membrane lipid peroxidation, and the change in its content is one of the important indicators of the degree of plasma membrane damage [66]. Our study found that LTF significantly reduced MDA content, as the expression of lipid-associated proteins (A2ZFY8 and A2X7M2) was significantly lower than that of LT (Additional file 3: Dataset S3). This finding indicated that LT induced disintegration of leaf cell membranes and electrolyte leakage in the plants, whereas LTF could contribute to mitigating the involvement of lipoxygenases in membrane lipid peroxidation and reducing damage to the cell membrane system.

Effect of flooding on metabolic regulation under low-temperature stress

Abscisic acid, as an endogenous hormone, is widely found in plants and plays an important role in enhancing plant stress tolerance [67]. The results indicated that the ABA content and the expression of ABA-binding protein (B8APF3) were significantly upregulated in the comparison of the LT and CK groups (Additional file 1: Dataset S1). Ma et al. [68] reported that overexpression of basic leucine zipper could promote the synthesis of ABA and improve the resistance of maize to low-temperature stress, and the results of our study were consistent with this finding. When the receptor cells of leaves and roots sense low-temperature stress, they induce gene expression and protein synthesis related to ABA precursor substances in plants, resulting in a large amount of ABA synthesis. ABA induces the expression of cold-resistance genes and promotes the synthesis of cold-resistance-specific proteins, thereby enhancing plant resistance to low-temperature stress. The dynamic balance regulation mechanism of different ABA contents in plants during growth and development can respond to different degrees of stress [69]. In contrast, the ABA content and ABA-associated protein expression (A2ZKG0) in LTF were significantly lower than those in LT (Additional file 3: Dataset S3). It is possible that flooding has a warming effect on the roots and stems of the aboveground parts. Agural et al. [70] reported that plant hormones, such as growth hormone and cytokinin, also affect plant growth and development during low-temperature stress. Our study showed that the GA3 content and GA3 synthesis-related proteins (A2XJL0 and B8BCQ6) in LT were significantly lower than those in CK (Additional file 1: Dataset S1). Low-temperature stress limited the activities of key enzymes in the GA3 signal transduction pathway and GA3 synthesis process, which decreased the content of GA3 and inhibited the normal growth and development of plants [71]. GA3 is a hormone that efficiently promotes cell division and elongation of roots and leaves, and plays an important role in the entire life cycle of plants [72]. In this study, the GA3 content and cell elongation protein expression (A2Z6Y9) were significantly upregulated in the LTF and LT groups (Additional file 3: Dataset S3). This finding may be observed because flooding promoted the transmission of GA3 signaling molecules in the signaling pathways under low-temperature stress and regulated the expression of GA3 biosynthesis-related genes and the synthesis of specific proteins, thereby facilitating the synthesis of a large number of GA3. In addition, soluble protein is an important osmotic regulator substance [73, 74]. In this study, LTF significantly decreased the soluble protein content compared to LT. This finding might be attributable to the decrease in the low-temperature stress effect after flooding and the response of plant osmoregulatory mechanisms to changes in environmental factors; thus, producing and accumulating less soluble proteins could ensure normal plant growth and development.


In this study, transcriptomic, proteomic and physiological analyses were employed to investigate the molecular mechanism by which flooding mitigates low-temperature stress. The damage that low-temperature stress caused to chloroplast structures in rice leaves was alleviated by flooding. LTF significantly increased the contents of rubisco, chlorophyll, PEPCK, ATP and GA3 and enhanced photosynthesis and energy metabolism compared with LT. GO and KEGG enrichment analysis demonstrated that the DEGs and DEPs were significantly associated with the photosynthesis pathway, antioxidant enzymes and metabolic regulation. Our results may help to elucidate the physiological characteristics, and gene and proteins expression changes by which flooding mitigates low-temperature stress.


Plant materials and growth conditions

Based on previous experiments [75], the early indica rice cultivar Zhongjiazao 17 (ZJZ17) was selected for pot experiments, which is mainly generalized for rice production in Jiangxi Province. Zhongjiazao 17 bred by the China National Rice Research Institute (CNRRI) is inbred indica rice. The pot experiment was conducted at the Key Laboratory of Crop Physiology of Jiangxi Agricultural University (JXAU) of Nanchang City, Jiangxi Province, China (longitude: 115° 50′ E; latitude: 28° 46′ N). Rice was planted in plastic pots measuring 15.0 cm in height, 25.0 cm in length and 23.0 cm in width. Soil samples of the plastic pots were derived from the upper soil layer (0–20 cm) of the rice experiment paddy field at the experimental site located at the science and technology park of JXAU. The soil fertility properties were as the following: soil organic matter 30.35 g·kg− 1, total nitrogen (N) 2.4 g·kg− 1, available phosphorus (P) 25.17 mg·kg− 1, available potassium (K) 84.02 mg·kg− 1, and pH 6.1. The soil was naturally blown-dried, and crushed by a soil grinder (FT-1000A, Changzhou WIK Instrument Manufacturing Co. Ltd., China), and then sifted by a 100-mm mesh. Each pot filled approximately six kilograms of dry soil that was soaked in water for 2 weeks before direct seeding. Three grams compound fertilizer (N-P-K = 15-15-15%) was applied as base fertilizer 1 day before direct seeding. Other management measures were in accordance with local recommendations.

Experimental design

Germinated seeds were selected for direct seeding, and all pots were subsequently placed in an artificial climate chamber with a diurnal temperature of 27/23 °C (day/night). After 20 days of direct seeding, one-third pots of rice seedlings with three leaves continued to be kept in a suitable temperature artificial climate chamber 27/23 °C (day/night) as a control (CK), and the two-thirds pots of rice seedlings were equally moved to another artificial climate chamber for low temperature (LT) and low temperature flooding (LTF) treatment. For low temperature (LT), the diurnal temperature was set to 10/6 °C (day/night), the period of low temperature treatment was 3 days, the light intensity was 100 μmol·m− 2 s− 1, the relative humidity (RH) of the artificial climate chamber was 75%, and potting soil remained moist. During the treatment, the positions of each pot were rotated from time to time to avoid the influence of light on seedling growth. For the low-temperature flooding treatment (LTF), the temperature and treatment period were both the same as those in the LT treatment, and the rice plants were maintained in a 5–6 cm flood water layer. For the control (CK), rice seedlings were grown in an artificial climate chamber in which the diurnal temperature was 27/23 °C (day/night), which kept the soil moist. Each treatment consisted of three replicates, and 20 pots constituted each replication. Under LT, LTF, and CK conditions, the relevant agronomic data were shown in Table S1 (Additional file 10).

Transmission electron microscopic observation

After 3 days of low temperature and low temperature flooding, the middle and upper parts of the leaves were cut into 1 × 1-cm pieces with a blade and immediately placed into the 2.5% glutaric dialdehyde fixative, vacuumized until the sample was completely submerged, and fixed overnight at 4 °C; approximately 10 pieces were taken for each sample. The sample- processing methods were undertaken according to the method described by Huang [76].

Physiological and biochemical characteristic

Enzymatic activity and endogenous hormone contents were measured on the third day of low-temperature and low-temperature flooding treatment. In each replicate experiment, three pots of plants were selected for each treatment group, and a 3-g leaf sample was collected and immediately frozen in liquid N2 and stored at − 80 °C until extraction. PEPCK activity and rubisco, chlorophyll, ATP, ABA, and GA3 contents [77,78,79,80,81] were measured using assay kits according to the manufacturer’s instructions (Cominbio, Suzhou, China). Soluble protein content, peroxidase (POD) activity, superoxide dismutase (SOD) activity and malondialdehyde (MDA) content were determined by configuration of solution [82,83,84].

RNA isolation, library preparation, and sequencing

According to the RNA extraction scheme, the total RNA was separated using TRIzol reagent (Invitrogen, Waltham, MA, USA). RNA purity was checked by a Nano Photometer® spectrophotometer (IMPLEN, CA, USA). The integrity of the RNA was assessed using the RNA 6000 Nano Kit of the Bioanalyzer 2100 system (Agilent Technologies, CA, USA). The quality and quantity of the RNA samples were determined by a Nanodrop 2000 spectrophotometer (Thermo, USA).

The total amount of 3 μg RNA per sample was used as the input material for RNA sample preparation. Sequencing libraries were generated using the NEBNext® Ultra™ RNA Library Prep Kit for Illumina® (NEB, USA) following the manufacturer’s recommendations and index codes were added to attribute sequences to each sample. In brief, mRNA was purified from total RNA using poly-Toligo-attached magnetic beads. Fragmentation was carried out using divalent cations under elevated temperature in NEBNext First Strand Synthesis Reaction Buffer (5X). First strand cDNA was synthesized using random hexamer primers and M-MuLV Reverse Transcriptase (RNase H-). Second strand cDNA synthesis was subsequently performed using DNA polymerase I and RNase H. The remaining overhangs were converted into blunt ends via exonuclease/polymerase activities. After adenylation of the 3′ ends of DNA fragments, NEBNext adaptors with hairpin loop structures were ligated to prepare for hybridization. To preferentially select cDNA fragments 150 ~ 200 bp in length, the library fragments were purified with the AMPure XP system (Beckman Coulter, Beverly, USA). Then, 3 μL USER Enzyme (NEB, USA) was used with size- selected, adaptor-ligated cDNA at 37 °C for 15 min followed by 5 min at 95 °C before PCR. Then PCR was performed with Phusion High-Fidelity DNA polymerase, universal PCR primers and Index (X) Primer. Finally, PCR products were purified (AMPure XP system), and library quality was assessed with the Agilent Bioanalyzer 2100 system. After that step, the library preparation was sequenced on the Illumina HiSeq 2500/X platform, and 125/150-bp paired-end reads were generated. RSEM (version 1.2.15) was used to quantify the mRNA abundance.

Differential expression analysis for RNA-seq data

Expected number of fragments per kilobase of transcript per million fragments mapped (FPKM) values were obtained using Cufflink software (version 2.1.1), which were used as values for normalized gene expression. Differential expression analyses of two comparison groups were performed using the DESeq R package (1.10.1). DESeq provides statistical routines for determining differential expression in digital gene expression data using a model based on the beta negative binomial distribution. The resulting P-values were corrected by the Benjamini-Hochberg method. Genes with an adjusted P-value < 0.01 and fold change > 2 found by DESeq were considered to be differentially expressed.

Protein extraction, trypsin digestion and TMT labelling

The sample was first ground in liquid nitrogen, and then the powder was transferred to a 5-mL centrifuge tube. After that, four volumes of lysis buffer (8 M urea, 1% Protease Inhibitor Cocktail) were added to the cell powder, followed by sonication three times on ice using a high-intensity ultrasonic processor (Scientz, China). The samples were centrifuged at 12,000×g at 4 °C for 10 min and the debris was removed. Finally, the protein concentration was determined with a BCA kit (Beyotime, China) according to the manufacturer’s instructions.

For digestion, the protein solution was reduced with 5 mM dithiothreitol for 30 min at 56 °C and alkylated with 11 mM iodoacetamide for 15 min at room temperature in darkness. The protein sample was then diluted by adding 100 mM TEAB to urea concentrations less than 2 M. Finally, trypsin was added at a 1:50 trypsin-to-protein mass ratio for the first digestion overnight and a 1:100 trypsin-to-protein mass ratio for a second 4 h-digestion.

After trypsin digestion, the peptide was desalted by a Strata X C18 SPE column (Phenomenex, USA) and vacuum-dried. The peptide was reconstituted in 0.5 M TEAB and processed according to the manufacturer’s protocol for the TMT kit. The peptide mixtures were then incubated for 2 h at room temperature and pooled, desalted and dried by vacuum centrifugation.

LC-MS/MS analysis and database search

The tryptic peptides were dissolved in 0.1% formic acid (solvent A) and directly loaded onto a homemade reversed-phase analytical column (15 cm length, 75 μm i.d.). The gradient consisted of an increase from 6 to 23% solvent B (0.1% formic acid in 98% acetonitrile) over 26 min, 23–35% in 8 min and climbing to 80% in 3 min then holding at 80% for the last 3 min, all at a constant flow rate of 400 nL/min on an EASY-nLC 1000 UPLC system (Thermo, Waltham, USA).

The peptides were subjected to an NSI source followed by tandem mass spectrometry (MS/MS) in a Q ExactiveTM Plus (Thermo) coupled online to UPLC. The electrospray voltage applied was 2.0 kV. The m/z scan range was 350 to 1.800 for full scan, and intact peptides were detected in the Orbitrap at a resolution of 70,000. Peptides were then selected for MS/MS using the NCE setting of 28, and the fragments were detected in the Orbitrap at a resolution of 17,500. A data-dependent procedure was employed that alternated between one MS scan followed by 20 MS/MS scans with 15.0 s dynamic exclusion. Automatic gain control (AGC) was set at 5E4. Fixed first mass was set as 100 m/z.

The resulting MS/MS data were processed using the MaxQuant search engine (v. Tandem mass spectra were searched against the human UniProt database concatenated with the reverse decoy database. Trypsin/P was specified as a cleavage enzyme allowing up to 4 missing cleavages. The mass tolerance for precursor ions was set as 20 ppm in the first search and 5 ppm in the main search, and the mass tolerance for fragment ions was set as 0.02 Da. Carbamidomethyl on Cys was specified as a fixed modification and acetylation modification and oxidation on Met were specified as variable modifications. The FDR was adjusted to < 1%, and the minimum score for modified peptides was set as > 40.

Confirmation using qRT-PCR

To verify the accuracy of transcriptome and proteome results, 18 DEGs involved in LT, LTF and CK responses were selected for verification using qRT-PCR. The RNA used for qRT-PCR was the same with those used to construct cDNA library. The primers were designed using Primer 5.0 (Additional file 11: Table S2) and synthesized by Xiamen Life-Int Technology Co., Ltd. and 18 s (GenBank Accession NO. AK059783) acted as the reference gene.

Statistical analysis

The data were analysed by analysis of variance (SAS Institute Inc., Cary, NC), and the means of different treatments were examined by Tukey’s honest significant difference (HSD) test to compare the differences at the probability level of 0.05.

Availability of data and materials

The mass spectrometry proteomics data have been deposited to the ProteomeXchange Consortium ( via the PRIDE [85] partner repository with the dataset identifier PXD024034. RNA sequencing data are available upon request from the corresponding author.





Low temperature


Low temperature flooding


Differentially expressed proteins


Differentially expressed genes


Liquid chromatograph- mass spectrometry


Gene ontology


Kyoto encyclopedia of genes and genomes


Klusters of eukaryotic orthologous groups


Fold change


Biological process


Cellular component


Molecular function


Relative humidity


Ribulose bisphosphate carboxylase oxygenase




Superoxide dismutase




Phosphoenolpyruvate carboxykinase

GA3 :

Gibberellin A3


Abscisic acid


Fresh weight


Honst significant difference


  1. 1.

    Hadiarto T, Tran LSP. Progress studies of drought-responsive genes in rice. Plant Cell Rep. 2011;30(3):297–310.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  2. 2.

    Ruizsanchez M, Aroca R, Munoz Y, Polon R, Ruizlozano JM. The arbuscular mycorrhizal symbiosis enhances the photosynthetic efficiency and the antioxidative response of rice plants subjected to drought stress. J Plant Physiol. 2010;167(11):862–9.

    CAS  Article  Google Scholar 

  3. 3.

    Zeng D, Tian Z, Rao Y, Dong G, Yang Y, Huang L, et al. Rational design of high-yield and superior-quality rice. Nat Plants. 2017;3(4):17031.

    PubMed  Article  PubMed Central  Google Scholar 

  4. 4.

    Clerget B, Bueno CS, Quilty JR, Correa TQ, Sandro JM. Modifications in development and growth of a dual-adapted tropical rice variety grown as either a flooded or an aerobic crop. Field Crops Res. 2014;155:134–43.

    Article  Google Scholar 

  5. 5.

    Sun ZY, Wang XF, Yamamoto H, Zhang JQ, Tani H, Zhong GS, et al. Extraction of rice-planting area and identification of chilling damage by remote sensing technology: a case study of the emerging rice production region in high latitude. Paddy Water Environ. 2016;15(1):181–91.

    Article  Google Scholar 

  6. 6.

    Peyman S, Hashem A. Evaluation eighteen rice genotypes in cold tolerant at germination stage. World Appl Sci J. 2010;11:1476–80.

    Google Scholar 

  7. 7.

    Huang M, Ibrahim MD, Xia B, Zou Y. Significance, progress and prospects for research in simplified cultivation technologies for rice in China. J Agric Sci. 2011;149(4):487–96.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  8. 8.

    Huang M, Jiang LG, Zou YB, Zhang WX. On-farm assessment of effect of low temperature at seedling stage on early-season rice quality. Field Crops Res. 2013;141:63–8.

    Article  Google Scholar 

  9. 9.

    Bonnecarrere V, Borsani O, Diaz P, Capdevielle F, Blanco P, Monza J. Response to photoxidative stress induced by cold in japonica rice is genotype dependent. Plant Sci. 2011;180(5):726–32.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  10. 10.

    Caliskan S, Caliskan ME, Arslan M, Arioglu H. Effects of sowing date and growth duration on growth and yield of groundnut in a Mediterranean-type environment in Turkey. Field Crops Res. 2008;105:131–40.

    Article  Google Scholar 

  11. 11.

    Huang M, Zhang WX, Jiang LG, Zou YB. Impact of temperature changes on early-rice productivity in a subtropical environment of China. Field Crops Res. 2013;146:10–5.

    Article  Google Scholar 

  12. 12.

    Wang WQ, Chen Q, Hussain S, Mei JH, Dong HL, Peng SB, et al. Pre-sowing seed treatments in direct-seeded early rice: consequences for emergence, seedling growth and associated metabolic events under chilling stress. Sci Rep. 2016;6(1):19637.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  13. 13.

    Suralta RR, Kano-Nakata M, Niones JM, Inukai Y, Kameoka E, Tran TT. Root plasticity for maintenance of productivity under abiotic stressed soil environments in rice: progress and prospects. Field Crops Res. 2018;220:57–66.

    Article  Google Scholar 

  14. 14.

    Zhou Z, Vriezen W, Caeneghem WV, Montagu MV, Straeten DVD. Rapid induction of a novel ACC synthase gene in deepwater rice seedlings upon complete submergence. Euphytica. 2001;121(2):137–43.

    CAS  Article  Google Scholar 

  15. 15.

    Tong L, Yoshida T, Maeda T, Kimijima H. Effects of temperature, sowing depth and soil hardness on seedling establishment and yield of Cambodian rice direct-seeded in flood paddy fields. Plant Prod Sci. 2007;10(1):129–35.

    Article  Google Scholar 

  16. 16.

    Mohammad A, Sudhishri S, Das TK, Singh M, Bhattacharyya R, Dass A. Water balance in direct-seeded rice under conservation agriculture in north-western indo-gangetic plains of India. Irrig Sci. 2018;36.

  17. 17.

    Al ME, Cantrill LC, Overall RL, Sutton BG. Mechanism of low-temperature-induced pollen failure in rice. Cell Biol Int. 2013;34(5):469–76.

    Google Scholar 

  18. 18.

    Confalonieri R, Mariani L, Bocchi S. Analysis and modelling of water and near water temperatures in flooded rice. Ecol Model. 2005;183(2):269–80.

    Article  Google Scholar 

  19. 19.

    Zia MS, Salim M, Aslam M, Gill MA, Rahmatullah. Effect of low temperature of irrigation water on rice growth and nutrient uptake. J Agron Crop Sci. 2008;173(1):22–31.

    Article  Google Scholar 

  20. 20.

    Mishra A, Salokhe VM. Flooding stress: the effects of planting pattern and water regime on root morphology, physiology and grain yield of rice. J Agron Crop Sci. 2010;196(5):368–78.

    Article  Google Scholar 

  21. 21.

    Lal B, Gautam P, Mohanty SR, Raja R, Tripathi R, Shahid M, et al. Combined application of silica and nitrogen alleviates the damage of flooding stress in rice. Crop Pasture Sci. 2015;66(7):679–88.

    CAS  Article  Google Scholar 

  22. 22.

    Kuroha T, Nagai K, Gamuyao R, Wang DR, Furuta T, Nakamori M, et al. Ethylene-gibberellin signaling underlies adaptation of rice to periodic flooding. Science. 2018;361(6398):181–6.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  23. 23.

    Li HL, Yu Y, Guo JX, Li XZ, Rensing C, Wang G. Dynamics of the rice rhizosphere microbial community under continuous and intermittent flooding treatment. J Environ Manag. 2019;249:109326.

    CAS  Article  Google Scholar 

  24. 24.

    Hattori Y, Nagai K, Furukawa S, Song XJ, Kawano R, Sakakibara H, et al. The ethylene response factors SNORKEL1 and SNORKEL2 allow rice to adapt to deep water. Nature. 2009;460(7258):1026–30.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  25. 25.

    Kende H. Deepwater rice: a model plant to study stem elongation. Plant Physiol. 1998;118(4):1105–10.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  26. 26.

    Farrell TC, Fox KM, Williams RL, Fukai S. Genotypic variation for cold tolerance during reproductive develop ment in rice: screening with cold air and cold water. Field Crops Res. 2006;98(2):178–94.

    Article  Google Scholar 

  27. 27.

    Jackson MB, Ram PC. Physiological and molecular basis of susceptibility and tolerance of rice plants to complete submergence. Ann Bot. 2003;91(2):227–41.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  28. 28.

    Almeida DM, Almadanim MC, Lourenço T, Abreu IA, Saibo NJM, Oliveira MM. Screening for abiotic stress tolerance in rice: salt, cold, and drought. Methods Mol Biol. 1938;2016:155–82.

    Google Scholar 

  29. 29.

    Lesk C, Rowhani P, Ramankutty N. Influence of extreme weather disasters on global crop production. Nature. 2016;529(7584):84.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  30. 30.

    He H, Ma F, Yang R, Chen L, Jia B, Cui J, et al. Rice performance and water use efficiency under plastic mulching with drip irrigation. PloS One. 2013;8(12):e83103.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  31. 31.

    Cohen SP, Liu HX, Argueso CT, Pereira A, Cruz CC, Verdier V, et al. RNA-Seq analysis reveals insight into enhanced rice Xa7-mediated bacterial blight resistance at high temperature. PLoS One. 2017;12(11):e0187625.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  32. 32.

    Jithesh MN, Shukla PS, Kant P, Joshi J, Critchley AT, Prithiviraj B. Physiological and transcriptomics analyses reveal that ascophyllum nodosum extracts induce salinity tolerance in arabidopsis by regulating the expression of stress responsive genes. J Plant Growth Regul. 2019;38:463–78.

    CAS  Article  Google Scholar 

  33. 33.

    Vera JC, Wheat CW, Fescemyer HW, Frilander MJ, Crawford DL, Hanski I, et al. Rapid transcriptome characterization for a nonmodel organism using 454 pyrosequencing. Mol Ecol. 2008;17(7):1636–47.

    CAS  PubMed  Article  Google Scholar 

  34. 34.

    Marino D, Ariz I, Lasa B, Santamaria E, Fernandez IJ, Gonzalez MC, et al. Quantitative proteomics reveals the importance of N source to control glucosinolate metabolism in Arabidopsis thaliana and Brassica oleracea. J Exp Bot. 2016;67(11):3313–23.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  35. 35.

    Ji ZJ, Zeng YX, Liang Y, Qian Q, Yang CD. Proteomic dissection of the rice-fusarium fujikuroi interaction and the correlation between the proteome and transcriptome under disease stress. BMC Genomics. 2019;20(1):91.

    PubMed  PubMed Central  Article  Google Scholar 

  36. 36.

    Tse WK, Sun J, Zhang H, Law AY, Yeung BH, Chow SC. Transcriptomic and iTRAQ proteomic approaches reveal novel short-term hyperosmotic stress responsive proteins in the gill of the Japanese eel. J Proteome. 2013;89:81–94.

    CAS  Article  Google Scholar 

  37. 37.

    Liang C, Tian J, Liao H. Proteomics dissection of plant responses to mineral nutrient deficiency. Proteomics. 2013;13:624–36.

    CAS  PubMed  Article  Google Scholar 

  38. 38.

    Garcia V, Stevens R, Gil L, Gilbert L, Gest N, Petit J, et al. An integrative genomics approach for deciphering the complex interactions between ascorbate metabolism and fruit growth and composition in tomato. C R Biol. 2009;332(11):1007–21.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  39. 39.

    Yang YW, Chen HC, Jen WF, Liu LY, Chang MC. Comparative transcriptome analysis of shoots and roots of tng67 and tcn1 rice seedlings under cold stress and following subsequent recovery: insights into metabolic pathways, phytohormones, and transcription factors. PLoS One. 2015;10:e0131391.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  40. 40.

    Hussain S, Yin HQ, Peng SB, Khan FA, Khan F, Sameeullah M, et al. Comparative transcriptional profiling of primed and non-primed rice seedlings under submergence stress. Front Plant Sci. 2016;7:1125.

    PubMed  PubMed Central  Google Scholar 

  41. 41.

    Yan SP, Zhang QY, Tang ZC, Su WA, Sun WN. Comparative proteomic analysis provides new insights into chilling stress responses in rice. Mol Cell Proteomics. 2006;5(3):484–96.

    CAS  PubMed  Article  Google Scholar 

  42. 42.

    Gandin A, Koteyeva NK, Voznesenskaya EV, Edwards GE, Cousins AB. The acclimation of photosynthesis and respiration to temperature in the C3-C4 intermediate Salsola divaricata: induction of high respiratory CO2 release under low temperature. Plant Cell Environ. 2014;37(11):2601–12.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  43. 43.

    Sharma P, Sharma N, Deswal R. The molecular biology of the low-temperature response in plants. Bioessays. 2005;27(10):1048–59.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  44. 44.

    Janmohammadi M, Zolla L, Rinalducci S. Low temperature tolerance in plants: changes at the protein level. Phytochemistry. 2015;117:76–89.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  45. 45.

    Minami A, Yano K, Gamuyao R, Nagai K, Kuroha T, Ayano M, et al. Time-course transcriptomics analysis reveals key responses of submerged deepwater rice to flooding. Plant Physiol. 2018;176(4):3081–102.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  46. 46.

    Gunawardena TA, Fukai B. The interaction of nitrogen application and temperature during reproductive stage on spikelet sterility infield-grown rice. Aust J Agric Res. 2005;56:625–36.

    Article  Google Scholar 

  47. 47.

    Li Z, Mo W, Jia L, Xu YC, Tang W, Yang W, et al. Rice FLUORESCENT1 is involved in the regulation of chlorophyll. Plant Cell Physiol. 2019;60(10):2307–18.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  48. 48.

    Mishra A, Salokhe VM. Seedling characteristics and the early growth of transplanted rice under different water regimes. Exp Agr. 2008;44(3):365–83.

    Article  Google Scholar 

  49. 49.

    Liu T, Du QJ, Li SZ, Yang JY, Li XJ, Xu JJ, et al. GSTU43 gene involved in ALA-regulated redox homeostasis, to maintain coordinated chlorophyll synthesis of tomato at low temperature. BMC Plant Biol. 2019;19(1):323.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  50. 50.

    Bowes G. Growth at elevated CO2: photosynthetic responses mediated through rubisco. Plant Cell Environ. 2006;14(8):795–806.

    Article  Google Scholar 

  51. 51.

    Yu J, Chen L, Xu M, Huang B. Effects of elevated CO on physiological responses of tall fescue to elevated temperature, drought stress, and the combined stresses. Crop Sci. 2012;52(4):1848–58.

    CAS  Article  Google Scholar 

  52. 52.

    Cui L, Zou Z, Zhang J, Zhao Y, Yan F. 24-Epibrassinoslide enhances plant tolerance to stress from low temperatures and poor light intensities in tomato. Funct Integr Genomics. 2016;16(1):29–35.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  53. 53.

    Makino A, Sage RF. Temperature response of photosynthesis in transgenic rice transformed with ‘sense’ or ‘antisense’ rbcS. Plant Cell Physiol. 2007;48(10):1472–83.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  54. 54.

    Li Y, Gao Y, Xu X, Shen Q, Guo S. Light-saturated photosynthetic rate in high-nitrogen rice (Oryza sativa L.) leaves is related to chloroplastic CO2 concentration. J Exp Bot. 2009;60(8):2351–60.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  55. 55.

    Lin X, Zhou W, Zhu D, Zhang Y. Effect of SWD irrigation on photosynthesis and grain yield of rice. Field Crops Res. 2005;94(1):67–75.

    Article  Google Scholar 

  56. 56.

    Gautam P, Lal B, Raja R, Baig MJ, Haldar D, Rath L, et al. Post-flood nitrogen and basal phosphorus management affects survival, metabolic changes and anti-oxidant enzyme activities of submerged rice. Funct Plant Biol. 2014;41(12):1284–94.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  57. 57.

    Vighi IL, Benitez LC, Amaral MN, Moraes GP, Auler PA, Rodrigues GS, et al. Functional characterization of the antioxidant enzymes in rice plants exposed to salinity stress. Biol Plant. 2017;61(3):540–50.

    CAS  Article  Google Scholar 

  58. 58.

    Fang CX, Zhang PL, Jian X, Chen WS, Lin HM, Li YZ, et al. Overexpression of Lsi1 in cold-sensitive rice mediates transcriptional regulatory networks and enhances resistance to chilling stress. Plant Sci. 2017;262:115–26.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  59. 59.

    Sanghera GS, Wani SH, Hussain W, Singh NB. Engineering cold stress tolerance in crop plants. Curr Genomics. 2011;12(1):30–43.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  60. 60.

    Cui SX, Huang F, Wang J, Ma X, Cheng YS, Liu JY. A proteomic analysis of cold stress responses in rice seedlings. Proteomics. 2005;5:3162–72.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  61. 61.

    Li YS, Ou SL, Yang CY. The seedlings of different japonica rice varieties exhibit differ physiological properties to modulate plant survival rates under submergence stress. Plants. 2020;9(8):982.

    CAS  PubMed Central  Article  Google Scholar 

  62. 62.

    Gill SS, Tuteja N. Reactive oxygen species and antioxidant machinery in abiotic stress tolerance in crop plants. Plant Physiol Biochem. 2010;48(12):909–30.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  63. 63.

    Ali S, Huang Z, Li HX, Bashir MH, Ren SX. Antioxidant enzyme influences germination, stress tolerance, and virulence of Isaria Fumosorosea. J Basic Microbiol. 2013;53(6):489–97.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  64. 64.

    Liu J, Sun X, Xu F, Zhang Y, Zhang Q, Miao R, Zhang J, Liang J, Xu W. Suppression of OsMDHAR4 enhances heat tolerance by mediating H2O2-induced stomatal closure in rice plants. Rice. 2018;11:38.

    PubMed  PubMed Central  Article  Google Scholar 

  65. 65.

    Chapman JM, Muhlemann JK, Gayomba SR, Muday GK. Rboh-dependent ROS synthesis and ROS scavenging by plant specialized metabolites to modulate plant development and stress responses. Chem Res Toxicol. 2019;32(3):370–96.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  66. 66.

    Liu WY, Yu KM, He TF, Li FF, Zhang DX, Liu JX. The low temperature induced physiological responses of Avena nuda L., a cold-tolerant plant species. Sci World J. 2013:1–7.

  67. 67.

    Wang XQ, Kuang TY, He YK. Conservation between higher plants and the moss physcomitrella patens in response to the phytohormone abscisic acid: a proteomics analysis. BMC Plant Biol. 2010;10(1):192.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  68. 68.

    Ma HZ, Liu C, Li ZX, Ran QJ, Xie GN, Wang BM, et al. ZmbZIP4 contributes to stress resistance in maize by regulating ABA synthesis and root development. Plant Physiol. 2018;178(2):753–70.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  69. 69.

    Ma YL, Cao J, He JH, Chen QQ, Li XF, Yang Y. Molecular mechanism for the regulation of ABA homeostasis during plant development and stress responses. Int J Mol Sci. 2018;19(11):3643.

    PubMed Central  Article  CAS  Google Scholar 

  70. 70.

    Agurla S, Gahir S, Munemasa S, Murata Y, Raghavendra AS. Mechanism of stomatal closure in plants exposed to drought and cold stress: adaptation mechanisms and their applications. Surviv Strat Extreme Cold Desiccation. 2018;1081:215–32.

    CAS  Article  Google Scholar 

  71. 71.

    Wang YL, Cui YT, Hu GH, Wang XD, Chen HZ, Shi QH, et al. Reduced bioactive gibberellin content in rice seeds under low temperature leads to decreased sugar consumption and low seed germination rates. Plant Physiol Biochem. 2018;133:1–10.

    PubMed  Article  CAS  PubMed Central  Google Scholar 

  72. 72.

    Lee ZH, Hirakawa T, Yamaguchi N, Ito T. The roles of plant hormones and their interactions with regulatory genes in determining meristem activity. Int J Mol Sci. 2019;20(16):4065.

    CAS  PubMed Central  Article  PubMed  Google Scholar 

  73. 73.

    Theocharis A, Clement C, Barka EA. Physiological and molecular changes in plants grown at low temperatures. Planta. 2012;235(6):1091–105.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  74. 74.

    Zhu LW, Cao DD, Hu QJ, Guan YJ, Hu WM, Nawaz A, Hu J. Physiological changes and sHSPs genes relative transcription in relation to the acquisition of seed germination during maturation of hybrid rice seed. J Sci Food Agric. 2015;96(5):1764–71.

    PubMed  Article  CAS  Google Scholar 

  75. 75.

    Tang SQ. Effects of low temperature and water-logging stress at bus stage on growth characteristics and yield of early indica rice. Nanchang: Jiangxi Agric Univ. 2019.

  76. 76.

    Huang JL, Lun X, Liu CH. Transmission electron microscopic sample preparation method for rice leaf. J Chin Electron Microsc Soc. 2018;37:386–9.

    Google Scholar 

  77. 77.

    Song HS, Seo HM, Jeon JM, Moon YM, Hong JW, Hong YG, et al. Enhanced isobutanol production from acetate by combinatorial overexpression of acetyl-CoA synthetase and anaplerotic enzymes in engineered Escherichia coli. Biotechnol Bioeng. 2018;115(8):1971–8.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  78. 78.

    Qin LY, Xue YX, Fei Y, Zeng LF, Yang SS, Deng XP. Identification, evolution and expression analyses of ribulose-1,5-bisphosphate carboxylase/oxygenase small subunit gene family in wheat. Acta Physiol Plant. 2018;40(5):85.

    Article  CAS  Google Scholar 

  79. 79.

    Feng B, Huang X, Jiang D, Hua L, Zhuo Y, Wu D. Endoplasmic reticulum stress inducer tunicamycin alters hepatic energy homeostasis in mice. Int J Mol Sci. 2017;18(8):1710.

    PubMed Central  Article  CAS  Google Scholar 

  80. 80.

    Teng YB, Cui HQ, Wang MY, Liu XY. Nitrate reductase is regulated by circadian clock-associated1 in arabidopsis thaliana. Plant Soil. 2017;416(1):477–85.

    CAS  Article  Google Scholar 

  81. 81.

    Bai LQ, Liu YM, Mu Y, Anwar A, He CX, Yan Y, et al. Heterotrimeric G-protein γ subunit CsGG 3.2 positively regulates the expression of CBF genes and chilling tolerance in cucumber. Front. Plant Sci. 2018;9:488.

    Google Scholar 

  82. 82.

    Campion EM, Loughran ST, Walls D. Protein quantitation and analysis of purity. Methods Mol Biol. 2011;681:229–58.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  83. 83.

    Garcia TA, Zenteno ST, Peregrino AB, Yepiz PG. Hypoxia, reoxygenation and cytosolic manganese superoxide dismutase (cMnSOD) silencing in Litopenaeus vannamei: effects on cMnSOD transcripts, superoxide dismutase activity and superoxide anion production capacity. Dev Comp Immunol. 2010;34(11):1230–5.

    Article  CAS  Google Scholar 

  84. 84.

    Castrejon SE, Yatsimirsky AK. Cyclodextrin enhanced fluorimetric determination of malonaldehyde by the thiobarbituric acid method. Talanta. 1997;44:951–7.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  85. 85.

    Ma J, Chen T, Wu S, Yang C, Bai M, Shu K, et al. iProX: an integrated proteome resource. Nucleic Acids Res. 2019;47:1211–7.

    Article  Google Scholar 

Download references


We thank the workers who have been at the forefront of the fight against the COVID-19 pandemic, allowing us to focus on writing this paper.


This work has been supported by the National Natural Science Foundation of China (31760366), the National Key Research and Development Program (2016YFD0300501 and 2017YFND0301605), and Jiangxi Natural Science Foundation of China (20161BAB214171). The funding bodies played no role in the design of the study and collection, analysis, and interpretation of data and in writing the manuscript.

Author information




YZ conceived and designed the experiments. WW, JD, LC conducted experiments and analyzed the data. ZW, YZ and XT contributed new reagents or analytical tools. WW wrote the manuscript. QS, XP revised the manuscript. All authors read and approved the manuscript.

Corresponding authors

Correspondence to Ziming Wu or Yanhua Zeng.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1 : Dataset S1.

Summaries and expression patterns of DEPs between LT and CK.

Additional file 2 : Dataset S2.

Summaries and expression patterns of DEPs between LTF and CK.

Additional file 3 : Dataset S3.

Summaries and expression patterns of DEPs between LTF and LT.

Additional file 4 : Figure S1.

Functional categories of the identified differentially expressed proteins (DEPs) between LT and CK.

Additional file 5 : Figure S2.

Functional categories of the identified differentially expressed proteins (DEPs) between LTF and CK.

Additional file 6 : Figure S3.

Functional categories of the identified differentially expressed proteins (DEPs) between LTF and LT.

Additional file 7 : Figure S4.

DEPs protein–protein interaction analysis between LT and CK.

Additional file 8 : Figure S5.

DEPs protein–protein interaction analysis between LTF and CK.

Additional file 9 : Figure S6.

DEPs protein–protein interaction analysis between LTF and LT.

Additional file 10 : Table S1.

Effects of low-temperature and low temperature flooding on agronomic characters.

Additional file 11 : Table S2.

Primers used for qRT-PCR verification of differently expressed genes.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Wang, W., Du, J., Chen, L. et al. Transcriptomic, proteomic, and physiological comparative analyses of flooding mitigation of the damage induced by low-temperature stress in direct seeded early indica rice at the seedling stage. BMC Genomics 22, 176 (2021).

Download citation


  • Rice
  • Low temperature
  • Flooding
  • Proteome
  • Transcriptome
  • Physiological traits