Physiological analysis and transcriptome sequencing reveal the effects of combined cold and drought on tomato leaf
BMC Plant Biology volume 19, Article number: 377 (2019)
Co-occurrence of cold and drought stress can alter the response of plants at morphological, physiological and molecular levels, which finally affect crop production, more than individual stress. Understanding the responses of crop to combined stress is necessary to improve tolerance and maintain crop production especially in the field where combined stress frequently occurs. We aimed to clarify the underlying leaf physiological and molecular mechanisms of tomato by imposing combining cold and drought on one popular tomato cultivar ‘Jinlingmeiyu’ as an example.
The physiological and genetic responses were identified in tomatoes after 42 h exposure to control, cold, drought and combined treatments. As compared with control, water loss rate at the three stresses including cold, drought and combined stress significantly decreased until 40 min after taking samples from the plants. The content of H2O2, zeatin riboside (ZR) and melatonin in all stress treatments were significantly higher than the control. Drought stress alone and combined stress induced the accumulation of abscisic acid (ABA) and auxin (IAA) as compared with control. The individual cold and combined stress significantly decreased the maximum quantum efficiency of PSII (Fv/Fm), quantum yield of PSII (Fq′/Fm′) and electron transport rate (ETR). In total, 7141, 1850 and 7841 genes were involved in the stress response to cold, drought and their combination. Functional analysis of the stress-inducible genes provided more insights concerning the complex regulatory mechanisms that were involved in combined stress. The expression level of 12 genes were validated by quantitative real-time PCR (qRT-PCR).
We found that the expression of stress-specific genes changed with physiological variation, indicating the close crosstalk between physiological and genetic response especially under combined stress. This study provides new knowledge on the complex regulatory mechanism genes in tomato (‘Jinlingmeiyu’) leaf to abiotic stresses.
One sentence summary
Combination of cold and drought is a new state stress with cold as a dominant factor in tomato.
Plants suffer low temperature and water shortage in cold- and drought-prone regions in the world [1,2,3]. Cold stress has been reported to negatively affect plant growth and development in various crops, such as rice , wheat  and tomato [6, 7]. Likewise, drought stress is a natural challenge for crops including tomato with an adverse effect on production [1, 8, 9]. Due to the rapid and dynamic global environmental changes, plants grown under the field conditions more and more frequently face a combination of various abiotic stresses , such as cold and drought. Studies involving cold and drought stress normally address the stress individually, while only few studies combine cold and drought [10, 11].
Tomato (Solanum lycopersicum L.) is an economically important crop worldwide, which is sensitive to a series of abiotic stresses, especially extreme temperature and drought. Cold and drought have been reported to decrease the tomato production, the damage of which was aggravated as a consequence of global climate change [12, 13]. Simultaneous cold and drought negatively affected crop growth and hampered productivity [11, 14], leading to economy loss. Tomato seedlings could easily suffer low temperature and water deficit in the beginning of spring and the late of winter in many regions of China. For instance, the promotion of cultivar ‘Jinlingmeiyu’ in Jiangsu Province were restricted due to the cold and drought stress. Therefore, it is crucial to determine how tomatoes such as ‘Jinlingmeiyu’ respond to combined cold and drought stresses.
Relative water content (RWC) is a quantitative indicator of plant water status . Cold stress (3 °C for 16 h) and drought stress (without irrigation for 5 days) induced a significant drop in the RWC of tomato [6, 15], which was correlated with the water loss. Higher leaf RWC of tomato under cold stress was promoted by exogenous H2O2 pretreatment . Tomato plants can maintain acquired tolerance though increasing NADPH oxidase-mediated H2O2 to temperature stress such as heat stress . Drought stress resulted in an increase in reactive oxygen species (ROS) triggering a response by signal transduction pathways with H2O2 as secondary messenger, which was closely related to ABA (abscisic acid) . Phytohormones such as ABA, IAA (auxin), GA3 (gibberellin), ZR (zeatin riboside) and melatonin played a dynamic and important role in increasing the tolerance of plants to abiotic stress . For instance, both cold and drought stress induced ABA accumulation in plants [17, 18]. Previous studies clearly demonstrated the significance of IAA in cold−/drought-mediated plant growth, regulated by polar deployment and trafficking of auxin carriers [19, 20]. The reduction of GA biosynthesis in plant under abiotic stress contributed to growth restraint mediated by DELLA proteins . On the contrary, enhanced cold tolerance together with high ZR content was observed in transgenic tall fescue containing the Agrobacterium tumefaciens ipt gene . Melatonin is naturally an antioxidant in plants that can efficiently remove ROS . In addition, due to the high sensitivity of the photosynthesis to abiotic stress, lower net photosynthesis rate and decreased Fv/Fm (maximum quantum efficiency of PSII) were observed in tomato under cold (3 ± 2 °C for 16 h)  and drought without irrigation for 5 days . Generally, the plants under cold and drought stress share a number of common responses in the aspect of phytohormones being related to ROS, antioxidant response and signaling transduction. . However, changes in phytohormones content of tomato under combined cold and drought as compared with individual stress remained unclear and how the changes affect tomato photosynthesis and growth need to clarify.
As a consequence of the increased efficiency of RNA sequencing (RNA-seq), transcriptome analysis has identified many genes involved in cold or drought response in plants, such as rice , maize  and tomato [26, 27]. Profound changes in the gene expression in plants were induced by cold  and drought . For example, genes involved in signal transduction (e.g. mitogen-activated protein kinase or MAPK) and osmotic stress signal perception (e.g. MYB29) were both responsive to cold and drought conditions . Furthermore, WRKY transcription factors has been reported to play important roles in response to abiotic stresses including cold and drought [29, 30]. Some of the genes were induced only by cold, while some genes are drought-inducible only . A key characteristic of both cold and drought stress is the accumulation of ABA that can promote a series of plant response such as ABA-dependent gene expression [18, 31].
Crop physiology and genomics have provided useful knowledge to improve plant tolerance (Tuberosa & Salvi, 2006). To the best of our knowledge, the physiological and molecular response of tomatoes under combined cold and drought stress are scarce. This study focuses on the responsive mechanisms of tomato under cold, drought and their combination in relation to leaf water status, photosystem II (PSII) activity, antioxidants, phytohormones and key genes by transcriptome analyses. We aim to clarify the effect of cold, drought and their combination ranging from leaf physiological and molecular levels and expound the relationship between individual and combined stress response with tomato cultivar ‘Jinlingmeiyu’ as an example. Our hypotheses are that 1) combined cold and drought stress create a different stress response compared to individual stress with cold as a dominant stress factor; 2) combined cold and drought stress negatively affected tomato leaf physiological response regulated by genes involved in signal transduction and photosynthetic electron transport. This study provides overview responses of tomato ‘Jinlingmeiyu’ to combined stresses by integrating plant physiological and transcriptional approaches.
Tomato exposed to the three stresses for 42 h showed slower growth than the control (Fig. 1a). The leaf RWC after drought stress was significantly lower than in the control, while that after cold stress was significantly higher than drought and combined stress (Fig. 1b). The water loss rate after all stress exposures was significantly lower than control from 10 min to 40 min after sampling, but this trend disappeared in the late stage (from 120 min to 240 min) (Fig. 1c). In addition, all stress treatments significantly increased the H2O2 content in the leaves with the highest H2O2 content induced by cold and combined stress (Fig. 2).
The phytohormones were regulated very differently by the various stress treatments. The content of ABA and IAA were only upregulated when the tomato plants were exposed to drought alone or combined with cold (Fig. 3a, b), and for ABA more after the combined stress than after single drought stress. The content of GA3 was only upregulated when exposed to cold alone and combination stress (Fig. 3c), while ZR and melatonin were upregulated by all stress treatments but more if drought was involved alone or in combination with cold (Fig. 3d, e).
The Fv/Fm and Fq′/Fm′ (quantum yield of PSII) after cold and combined stress were significantly lower than control and drought stress (Fig. 4a, b). The qL (fraction of open PSII centers) showed the opposite pattern from NPQ (non-photochemical quenching) with cold increasing and drought decreasing qL but here the variation prevented any establishment of significant differences (Fig. 4c). After 42 h stress, drought increased and cold decreased NPQ making the cold and drought treatments significantly different with the control and combined stress in between (Fig. 4d). The electron transport rate (ETR) in the cold and the combined stress treatments decreased significantly compared with control at PPFD (photosynthetic photon flux density) from 0 to 822 μmol m− 2 s− 1 (Fig. 4e). By comparison, in the drought treatment ETR was significantly lower than control only at the highest PPFD (Fig. 4e).
The four cDNA libraries of tomato treated by control, cold, drought and their combination generated 25,082,239, 29,815,250, 28,754,112 and 27,195,304 clean reads with 43.3, 43.0, 42.8 and 43.3% GC percentage, respectively. The respective Q30 percentage (sequencing error rates < 0.1%) was 92.9, 93.0, 92.5 and 91.9% in tomato in response to control, cold, drought and their combination, indicating that the sequencing data was qualified for further analysis. In total, there were 50,162,479, 59,630,499, 57,508,225 and 54,390,609 total reads identified in tomato in response to the four treatments (Additional file 3: Table S1). De novo assembly predicted 755 new genes that have not appeared in tomato genome database annotated by COG (Clusters of Orthologous Groups, http://www.ncbi.nlm.nih.gov/COG/), GO (Gene Ontology Consortium, http://www.geneontology.org/), KEGG (Kyoto Encyclopedia of Genes and Genomes, http://www.genome.jp/kegg/), KOG (Eukaryotic Orthologous Groups, http://www.ncbi.nlm.nih.gov/KOG/), Pfam (Homologous Protein Family, http://pfam.xfam.org/), Swissprot (A manually annotated, non-redundant protein sequence database, http://www.uniprot.org/), eggNOG (evolutionary genealogy of genes: Non-supervised Orthologous Groups, http://eggnog.embl.de/) and nr (non-redundant protein sequence database, ftp://ftp.ncbi.nih.gov/blast/db/) (Additional file 1: Figure S1). The eggNOG function classification of consensus sequence showed that the function of the majority of the genes remained unknown and many genes played roles in the section of ‘L’ (replication, recombination and repair), ‘Q’ (secondary metabolites biosynthesis, transport and catabolism) and ‘T’ (signal transduction mechanisms) (Additional file 1: Figure S1).
RNA-seq provided an overview of genes differentially expressed in tomato during the different stresses (Fig. 5). Compared with control, 7141, 1850 and 7841 genes were significantly differently expressed in tomato after cold, drought and combined stress, respectively (Fig. 5a). Combined stress induced significantly different expression levels of 44 and 7886 genes in comparison with individual cold and drought, respectively (Fig. 5a). Among the significantly DEGs (differentially expressed genes), 1165 genes existed only in control vs combined stress as compared with control vs cold and control vs drought (Fig. 5b). Overall, the expression pattern of the genes in tomato during cold were similar to combined stress (Fig. 5c).
The genes with significantly different levels in tomato after combined stress had a role in carbon metabolism, biosynthesis of amino acids, plant-pathogen interaction and so on (Additional file 2: Figure S2). When compared to cold stress, genes expressed in response to combined stress played roles in phytohormone signal transduction, plant-pathogen interaction etc. (Fig. 6a). The genes responsive to combined stress differed from that of drought with a role in ribosomes, carbon metabolism, biosynthesis of amino acids, plant-pathogen interaction, glycolysis/gluconeogenesis etc. (Fig. 6b).
The expression levels of 12 genes with significantly changed expression levels based on RNA-seq were detected by qRT-PCR (quantitative real-time PCR validation). The drought stress significantly increased the expression level of Solyc12g099390.1 in comparison with control (Fig. 7a). The expression level of Solyc07g044840.2 in tomato after drought and combined stress was significantly higher than control (Fig. 7b). The cold, drought and combined stress significantly decreased the expression level of Solyc02g077990.2, Solyc05g052600.2, Solyc09g011810.2, Solyc12g009600.1, Solyc06g051400.2 and Solyc01g103100.2, but significantly increased the expression level of Solyc09g009020.2 (Fig. 7c-i). Compared with the control, the expression level of Solyc01g028810.2, Solyc11g069380.1 and Solanum_lycopersicum_newGene_4540 (Solyc_newGene_4540) in tomato dropped significantly after cold and combined stress, while the expression level of Solyc11g069380.1 and Solyc_newGene_4540 in tomato increased significantly after drought stress (Fig. 7j-l). The primary functions of the 12 verified genes were signal transduction, carbohydrate transport and metabolism, translation as well as posttranslational modification responding to drought and cold (Fig. 7). The details of the 12 genes were shown Additional file 3: Table S2.
Plants respond to cold and drought stress through a series of biochemical, physiological and genetic processes [11, 32]. Molecular and genomic studies have shown that the genes with a wide variety of functions are induced by cold and drought stresses . Understanding the mechanisms of physiology and genes involved in stress signaling networks and stress tolerance is the key for crop improvement and production [3, 11].
The physiological responses were closely related to key genes in tomato during stress
The common phenotypic symptoms of low temperature and water shortage in plants are leaf wilting, stunted growth and even death during the exposure [3, 17]. In accordance with our previous studies [15, 33], drought stress decreased the RWC and water loss rate in tomato (Fig. 1). Plants exposed to cold stress often show water-stress symptoms due to cold-induced inhibition of water uptake and loss [6, 17]. Thereby, leaf wilting and curling in tomato in response to three stresses were correlated with low RWC and decreased water loss rate, especially in the early stage after sampling.
In addition to changes in water status, increased H2O2 content was observed in tomato responding to three stresses (Fig. 2), indicating that low temperature and water shortage conditions caused ROS generation. With the chosen stress levels, the chilling temperature created more H2O2 than the drought. We found that the expression level of genes being involved in the ROS response were significantly changed by the three stresses (Additional file 3: Table S3). For example, as compared with control, the expression of Solyc02g086050.2 and Solyc12g070270.1 significantly decreased by 1.76 and 1.48, while that of Solyc02g078360.2 significantly increased by 1.29 in tomato leaf under cold stress. Similarly, the three genes involved in ROS responding to cold stress showed the same trend in expression levels to combined stress. Besides, the expression level of Solyc10g044550.1 was down-regulated (− 1.61) in tomato leaf under drought condition in comparison with control. The H2O2 may be part of the signal to turn on the stress-responsive genes in tomato under cold and/or drought stress as indicated by Cruz de Carvalho (2008)  and Yadav (2011) . The H2O2 can act as an oxidative burst signal, activating the stress defense system in plants . In this study, 829, 175 and 907 genes with roles in signal transduction mechanisms were induced by cold, drought and combined stress, respectively (Additional file 3: Table S4). Numerous transcription-factor genes were stress-inducible and different transcriptional regulatory mechanisms were induced by cold and/or drought stress in signal transduction pathways [35, 36]. Accordingly, in this study, cold, drought and combined stress induced transcription factors such as members of the basic-domain leucine zipper (bZIP) family, the MYB family and the WRKY family. These transcription factors can regulate many stress-inducible genes and establish gene networks in response to various stress conditions in tomato.
Moreover, phytohormones are crucial components in various signal pathways , and for instance, ABA can stimulate the synthesis of H2O2 in guard cell and in turn H2O2 mediated ABA-dependent stomatal closure . A significant crosstalk among the drought and ABA responses have been shown since more than half of the drought-inducible genes are also induced by ABA treatments . Taking account into the increased ABA content in tomato after drought and combined stress (Fig. 3), it is suggested that there was a positive relationship between the water deficiency and ABA enhancement. A similar response of ABA and IAA (growth promoting hormone) was observed in tomato in all stresses, which was consistent with the result in cold-shock treated wheat . Reduction of GA content and signaling adversely affected plant growth under abiotic stress such as cold . However, it is opposite with our finding that cold induced the GA3 accumulation regardless of dourght stress, which could be due to the effect of other hormonal responses. ZR contents in ipt-transgenic lines with enhanced cold tolerance were 65–400% higher than non-transformed control , suggesting that the studied tomato actively responded to the stress conditions by increasing ZR content. Previous study found that melatonin was involved in the regulation of ROS and ABA in Malus species when exposed to drought . Increased melatonin content in tomato during drought stress help to alleviate the ROS damage due to comparatively lower H2O2 content as compared with cold and combined stress. Therefore, it is suggested that inducing a signal transduction such as the production of H2O2 and phytohormones is a key mechanism for plants to respond to abiotic stress. ABA synthesis is one of the quickest responses of plants under abiotic stress that trigger the expression of ABA-inducible gene . We found that 120, 43 and 125 genes playing roles in the phytohormone signal transduction were significantly regulated by cold, drought and combined stress, respectively (Additional file 3: Table S5). This can be partly explained by the significant changes in phytohormone contents in tomato under the stress conditions (Fig. 3).
The tomato exhibited significant variation in PSII activity especially after exposure to cold and combined stress. According to Baker and Rosenqvist (2004) , as expected, drought alone did not affect neither Fv/Fm nor the fluorescence quenching parameters measured at the growth PPFD (Fig. 4). Any stress that involved 10 °C temperature, however, significantly decreased Fv/Fm and the operating efficiency, Fq′/Fm′, at the growth PPFD. Neither the oxidation state of PSII, qL, and the light controlled heat dissipation in the antenna, NPQ, were significantly affected at the growth PPFD (Fig. 4). However, in light of the increasing difference between the treatments with increasing PPFD seen in the rapid light response curve of ETR, it implied that the tendencies seen at intermediate PPFD become significantly different when approaching light saturation. Meanwhile, we found that 36, four and 38 genes with significantly different expression levels in tomato under cold, drought and combined stress, respectively, were involved in the photosynthetic electron transport (Additional file 3: Table S6). Such physiological and genetic responses in tomato under cold or combined stress in turn lead to ROS generation and oxidative damage if protective mechanisms could not dissipate the excessive energy, which partially explained why tomato under cold and combined stress showed higher H2O2 content than drought.
Combined cold and drought was a new state of stress condition with cold as a dominant factor
Deng et al. (2012) found that the physiological changes induced by drought are quite different from that induced by cold . This was confirmed in our study by the physiological parameters such as RWC, water loss rate, H2O2 content, phytohormone content and chlorophyll fluorescence. Sales et al. (2013) found that drought stress alone (11 days without irrigation) severely affected ETR in sugarcane . We found that ETR was more affected by cold and combined stress than drought, since an increase in photorespiration might have protected ETR from damage caused by drought. More importantly, in most cases here, the physiological responses were similar between cold and combined stress, suggesting the dominant role of low temperature when cold was combined with drought. The irrigation was stopped three days before temperature stress since water deficit is a gradual process that took time to happen. The possible reason for the dominant role of cold was that the response of tomato to drought was slow due to gradual water deficit, while the response of tomato to low temperature was very fast.
Only 10% of the drought-inducible genes were also induced by cold , being consistent with our results that a few genes were shared in tomato by individual cold and drought condition. Apart from the physiological response, 4111 and 3457 genes were up- and down- regulated by drought as compared with cold (Fig. 5), indicating different molecular strategies in their reaction to the two stresses. Cold and drought induced common stress-inducible genes, while one of the stresses specifically induced some genes . By contrast, a few genes with different expression levels between cold and combined stress were also found (Fig. 5). For example, the genes involved in ROS response showed the similar pattern in expression level when the plants were exposed to cold and combined stress, but exhibited a different pattern under drought and combined stress (Additional file 3: Table S3). As compared with individual cold and drought, there were three and 852 specific genes involved in signal transduction mechanisms responding to combined stress, respectively (Additional file 3: Table S4). This explained why a few genes played roles in tomato under combined stress in comparison with cold (Fig. 6). In our previous study, drought stress played a dominant role when combined with heat stress . Obviously, it is not the same story when drought was combined with cold as indicated by the complex physiological and molecular networks responding to combined cold and drought. Combined cold and drought stress induced dehydration and declined tomato photosynthetic electron transport followed by ROS metabolism disorder and increased hormone level, during which the activation of stress-responsive genes and transcription factors happened (Fig. 8). However, the photosynthetic electron transport was less affected after drought than cold and combined stress due to an increasing photorespiration, which was accompanied with lower H2O2 and GA3 content (Fig. 8). The fact that the overlapping genes in response to individual and combined stress in tomato as shown in Fig. 8 highlights their importance in enhancing the tolerance of plants to abiotic stress. By comparison, specific genes responding to individual and combined stress, respectively, indicated that combined stress is a new state of condition with the activation of specifically responsive genes (Fig. 8).
In general, combining cold and drought stress was a new state of stress condition rather than a simple addition of individual stresses on tomato ‘Jinlingmeiyu’ seedlings. The crosstalk between leaf physiology and gene expressions involved in ROS, phytohormone, signal transduction, photosynthetic electron transport and carbohydrate transport in tomato at individual and combined stress induced the disorder in carbohydrate metabolism, which contributed to delayed growth and development in plant under abiotic stresses (Fig. 8). This study provides us insights into the leaf physiological and molecular mechanisms of tolerance and adaption to cold, drought and their combination in plants with tomato ‘Jinlingmeiyu’ as an example.
Plant material and growth condition
Single seeds of the popular commercial tomato cultivar ‘Jinlingmeiyu’ bred by Jiangsu Province Academy of Agricultural Sciences were sown in plug trays. The seedlings were grown at 26/18 °C (14 h/10 h, day/night) and 50–60% relative humidity in climate chambers (RDN-560E-4, Dongnan Instrument Co, Ltd., Ningbo, China). The light intensity was 300 μmol m− 2 s− 1 PPFD provided by white LED light (Dongnan Instrument Co, Ltd., Ningbo, China) at the level of the plants. The plants were irrigated by water every five days after sowing and irrigated by nutrient solution every three days 15 days after sowing. The nutrient solution was made up based on the Japanese Garden test formula [Ca(NO3)2•4H2O, 945 mg L− 1; KNO3, 809 mg L− 1; NH4H2PO4, 153 mg L− 1; MgSO4•7H2O, 493 mg L− 1; FeSO4•7H2O, 13.9 mg L− 1; Na2-EDTA, 18.6 mg L− 1; H3BO3, 2.86 mg L− 1; MnSO4•4H2O, 2.13 mg L− 1; ZnSO4•7H2O, 0.22 mg L− 1; CuSO4•5H2O, 0.08 mg L− 1; (NH4)6Mo7O24•4H2O, 0.02 mg L− 1]. The 21-days-old uniform plants were transferred to plastic pots (11 cm diameter, 9 cm height) with one plant per pot. Both the plug trays and plastic pots were filled with sphagnum substrate (Pindstrup, Denmark). After three days, we stopped irrigating half of the seedlings. The 27-day-old irrigated plants were treated at continuous 26 °C and 10 °C for 42 h at 300 μmol m− 2 s− 1 as control and cold stress treated plants, respectively. The 27-day-old non-irrigated plants were treated at 26 °C and 10 °C for 42 h at 300 μmol m− 2 s− 1 as drought and combined stress treated plants, respectively.
Leaf relative water content and water loss rate
The last fully expanded leaf from the top was cut from the plants after 42 h of the treatments with four replications. The fresh weight (FW) of the leaf was measured immediately. The leaf was immersed in dd-H2O and incubated at room temperature for four hours. Turgid weight (TW) of the leaf was measured after taking the leaf out form dd-H2O and blotting water from the leaf surface. Dry weight (DW) of the leaf was measured by drying the leaf at 80 °C for 24 h. The relative water content (RWC) of the leaf in % was calculated as RWC = [(FW - DW)/(TW - DW)] * 100. The FW of the leaf immediately after cutting (0 min), 10 min, 20 min, 30 min, 40 min, 60 min, 90 min, 2 h, 3 h and 4 h after cutting were measured at room temperature. Water loss rate of the leaf x min after cutting in % = [(FW0 min - FWx min)/FW0 min] * 100.
Leaf H2O2 content
The last fully expanded leaf from three plants per treatment were taken for the measurements as three replications. The absorbance was detected at 560 nm using a UV spectrophotometer (TU-1810, Beijing purkinje general instrument, Beijing, China) using ferrous oxidation xylenol orange assay. In detail, 2 g leaf sample was grounded at 4 °C in 2 mL pre-cool acetone and centrifuged at 10,000 r min− 1 at 4 °C for 10 min. Then, 1 mL supernatant was well mixed with 2.25 mL carbon tetrachloride and 0.75 mL trichloromethane. Afterwards, 5 mL ddH2O was added, well mixed and centrifuged at 5000 r min− 1 for 1 min. After that, 1 mL supernatant was well mixed with 2 mL mixture of reagent A and B (1:10, volume ratio). The reagent A contained 3.3 mmol L− 1 ferrous sulfate, 3.3 mmol L− 1 ammonia sulfate and 412.5 mmol L− 1 sulfuric acid. The reagent B contained 165 μmol L− 1 xylenol orange and 165 mmol L− 1 sorbitol.
Leaf phytohormone content
The last fully expanded leaf from the top were collected with three replications and immediately frozen in liquid nitrogen for the measurements of ABA, IAA, GA3, ZR and melatonin content using ELISA technique. As described by He (1993)  and Yang et al. (2001) , 2 mL pre-cold 80% (v/v) methanol with 1 mmol L− 1 BHT (butylated hydroxytoluene) were added. The samples were homogenized and centrifuged at 10000 r min− 1 at 4 °C. Afterwards, the supernatant were passed through a C-18 solid phase extraction column (Waters, Milford, MA) and dried in N2. The samples were dissolved in 0.01 mol L− 1 PBS (pH 7.4) to determine the levels of phytohormone content. By adding known amounts of standard hormone, the percentage recovery of the ABA, IAA, GA3 and ZR were calculated. The samples for melatonin content was homogenized and determined according to the plant MT ELISA kit (Lanpai Bio, Shanghai, China). In detail, the standard and samples were added in the wells, after which 100 μL HRP horseradish peroxidase were added. The samples were covered by adhesive strip and incubated for 1 h at 37 °C. The wells were aspirated and washed, which was repeated five times. Subsequently 50 μL substrate A and 50 μL substrate B were added to the wells, which were incubated for 15 min at 37 °C in darkness. Afterwards, 50 μL stop buffer was added to the wells. The optical density (O.D.) at 450 nm was recorded using a microtiter plate reader within 15 min. The melatonin content in the samples was calculated by comparing the O.D. of the samples from a standard curve.
Leaf chlorophyll fluorescence
The last fully expanded leaf from the top of four plants per treatment were used for the measurements after 42 h treatment. The seedlings were put in darkness for 20 min before quenching analysis. The Fv/Fm, Fq′/Fm′, qL and NPQ under a PPFD of 285 μmol m− 2 s− 1 with internal light source were measured using MINI-PAM-II (Walz, Eiffeltrich, Germany) and calculated according to the summary by Murchie and Lawson (2013) . Rapid light-response curves of ETR were obtained by exposing the leaf sample for 30 s to each of 10 increasing actinic light levels from 0 to 1496 μmol m− 2 s− 1 PPFD (0, 24, 45, 65, 90, 126, 190, 285, 420, 630 and 822 μmol m− 2 s− 1) using MINI-PAM-II. The main vein of the leaf was avoided during the measurements.
Library construction and RNA-seq of leaf samples
The first fully expanded leaf from the top of the plants was collected with three replicates after 42 h per treatment. The 12 samples were immediately frozen in liquid nitrogen and stored at − 80 °C before sequencing. Total RNA was extracted using the Total Plant RNA Extraction Kit (Karroten, Beijing, China). The RNA concentration and integrity of the RNA was detected using NanoDrop 2000 (Thermo Scientific, DE, USA) and RNA Nano 6000 of Agilent Bioanalyzer 2100 system (Agilent Technologies, CA, USA) to ensure the qualified sample for RNA-seq. The libraries were constructed using NEBNext UltraTM RNA Library Prep Kit for Illumina (NEB, USA). In detail, the mRNA was captured poly-T oligo-attached magnetic beads (NEB, USA). The mRNAs were randomly fragmented using NEBNext First Strand Synthesis Reaction Buffer (NEB, USA). The first strand of cDNA was synthesized using random hexamers and the second strand of cDNA was synthesized by adding buffer, dNTPs, RNase H and DNA polymerase I with the mRNA as a template. The cDNAs were purified using AMPure XP system (Beckman Coulter, Beverly, USA), after which the distal ends were repaired and the tails as well as adapters were linked. To select the fragment size, the AMPure XP beads were used and cDNA library was obtained by PCR enrichment.
The quality of the cDNA library were determined using Agilent Bioanalyzer 2100 (Agilent Technologies, Palo Alto, CA) and qPCR (BioMarker Technologies Co. Ltd., Beijing) to make sure the cDNA library was qualified for sequencing. Then, the qualified cDNA library was sequenced using the HiSeq X-ten (BioMarker Technologies Co. Ltd., Beijing).
Identification of differentially expressed genes (DEGs) and functional annotation in leaf of tomato ‘Jinlingmeiyu’
RNA-Seq data for samples of plants treated with the four treatments were obtained from three biological replicates. Empty reads, adapter sequences and low-quality sequences were removed from raw reads to obtain clean reads. The transcript abundances of genes in the samples were calculated by fragments per kilobase of exon per million fragments mapped (FPKM). The DEGs in tomatoes in response to the four treatments were obtained using bioinformatics methods. The criteria for DEGs selection was fold change ≥2 and false discovery rate (FDR) < 0.05. The function of the DEGs were conducted based on GO and KEGG database.
Total RNAs were extracted using Trizol reagent (Invitrogen, CA, USA). Using Prime Script RT reagent Kit (TaKaRa, Dalian, China), the qualified total RNAs were reverse transcribed to cDNA at 42 °C for 60 min, 70 °C for 15 min and kept on ice for 3 min in an Eppendorf Mastercycler Gradient (Mastercycler®ep realplex, Hamburg, Germany). Using SYBR Premix Ex TaqTM (Takara, Dalian, China), the cDNA was amplified at 95 °C for 2 min, 40 cycles of 95 °C for 15 s, 60 °C for 15 s and 72 °C for 20 s in an Eppendorf real-time PCR machine (Mastercycler®ep realplex, Hamburg, Germany). The reactions were repeated three times as three technical repetitions with three biological repetitions. The expression levels were calculated by the 2−ΔΔCt method with EF1α as reference gene. The primers for qRT-PCR were designed using primer 5.0 (Primer-E Ltd., Plymouth, UK) shown in Additional file 3: Table S7.
Data analysis and data access
Analysis of variance (ANOVA) between the physiological parameters of plants at the control, cold, drought and combined stress were performed using SPSS 16.0 (SPSS Inc. Chicago, IL, USA). The qRT-PCR data of the 12 genes in tomatoes at the four treatments were analyzed by ANOVA. The KEGG pathway enrichment analysis of the DEGs was performed using hypergeometric distribution (BioMarker Technologies Co. Ltd., Beijing). The significant differences were considered at P < 0.05. The sequencing data has been submitted to NCBI under the accession number of SRP156535 at https://www.ncbi.nlm.nih.gov/sra/SRP156535.
Availability of data and materials
The sequencing data has been submitted to NCBI under the accession number of SRP156535 at https://www.ncbi.nlm.nih.gov/sra/SRP156535.
The seeds, plant materials and datasets during or analyzed during the current study available from the corresponding author on reasonable request.
differentially expressed genes
electron transport rate
false discovery rate
fragments per kilobase of exon per million fragments mapped
- Fq ′/Fm ′ :
quantum yield of photosystem II
- Fv/Fm :
maximum potential quantum efficiency of photosystem II
- GA3 :
gene ontology consortium
kyoto encyclopedia of genes and genomes
photosynthetic photon flux density
- qL :
fraction of open PSII centers
quantitative real-time PCR
relative water content
Cattivelli L, Rizza F, Badeck FW, Mazzucotelli E, Mastrangelo AM, Francia E, Marè C, Tondelli A, Stanca AM. Drought tolerance improvement in crop plants: an integrated view from breeding to genomics. Field Crop Res. 2008;105:1–14.
Inman-Bamber NG, Smith DM. Water relations in sugarcane and response to water deficits. Field Crop Res. 2005;92:185–202.
Yadav SK. Cold stress tolerance mechanisms in plants. A review. Agron Sustain Dev. 2011;30:605–20.
Hussain S, Khan F, Cao W, Wu L, Geng M. Seed priming alters the production and detoxification of reactive oxygen intermediates in rice seedlings grown under sub-optimal temperature and nutrient supply. Front Plant Sci. 2016;7:439.
Li X, Cai J, Liu F, Zhou Q, Dai T, Cao W, Jiang D. Wheat plants exposed to winter warming are more susceptible to low temperature stress in the spring. Plant Growth Regul. 2015;77:11–9.
İşeri ÖD, Körpe DA, Sahin FI, Haberal M. Hydrogen peroxide pretreatment of roots enhanced oxidative stress response of tomato under cold stress. Acta Physiol Plant. 2013;35:1905–13.
Starck Z, Niemyska B, Bogdan J, Tawalbeh RNA. Response of tomato plants to chilling stress in association with nutrient or phosphorus starvation. Plant Soil. 2000;226:99–106.
Parry M, Flexas J, Medrano H. Prospects for crop production under drought: research priorities and future directions. Ann Appl Biol. 2005;147:211–26.
Petrozza A, Santaniello A, Summerer S, Di Tommaso G, Di Tommaso D, Paparelli E, Piaggesi A, Perata P, Cellini F. Physiological responses to Megafol® treatments in tomato plants under drought stress: a phenomic and molecular approach. Sci Hortic. 2014;174:185–92.
Mittler R. Abiotic stress, the field environment and stress combination. Trends Plant Sci. 2006;11:15–9.
Hussain HA, Hussain S, Khaliq A, Ashraf U, Anjum SA, Men S, Wang L. Chilling and drought stresses in crop plants: implications, cross talk, and potential management opportunities. Front Plant Sci. 2018;9:393.
Kalloo G. Breeding for environmental stress resistance in tomato. In: Kalloo G, editor. Genetic improvement of tomato. Berlin: Springer; 1991. p. 153–65.
Park EJ, Jeknić Z, Sakamoto A, DeNoma J, Yuwansiri R, Murata N, Chen THH. Genetic engineering of glycinebetaine synthesis in tomato protects seeds, plants, and flowers from chilling damage. Plant J. 2004;40:474–87.
Tommasini L, Svensson JT, Rodriguez EM, Wahid A, Malatrasi M, Kato K, Wanamaker S, Resnik J, Close TJ. Dehydrin gene expression provides an indicator of low temperature and drought stress: transcriptome-based analysis of barley (Hordeum vulgare L.). Functional and Integrative Genomics. 2008;8:387–405.
Zhou R, Yu X, Ottosen CO, Rosenqvist E, Zhao L, Wang Y, Yu W, Zhao T, Wu Z. Drought stress had a predominant effect over heat stress on three tomato cultivars subjected to combined stress. BMC Plant Biol. 2017;17:24.
Sun M, Jiang F, Cen B, Wen J, Zhou Y, Wu ZJPC. Environment: respiratory burst oxidase homologue-dependent H2O2 and chloroplast H2O2 are essential for the maintenance of acquired thermotolerance during recovery after acclimation. Plant Cell and Environment. 2018;41:2372–89.
Cruz de Carvalho MH. Drought stress and reactive oxygen species: production, scavenging and signaling. Plant Signal and Behavior. 2008;3:156–65.
Liu Q, Kasuga M, Sakuma Y, Abe H, Miura S, Yamaguchi-shinozaki K, Shinozaki K. Two transcription factors, DREB1 and DREB2, with an EREBP/AP2 DNA binding domain separate two cellular signal transduction pathways in drought- and low-temperature-responsive gene expression, respectively, in Arabidopsis. Plant Cell. 1998;10:1391–406.
Rahman A. Auxin: a regulator of cold stress response. Physiol Plant. 2013;147:28–35.
Shi H, Chen L, Ye T, Liu X, Ding K, Chan Z. Modulation of auxin content in Arabidopsis confers improved drought stress resistance. Plant Physiol Biochem. 2014;82:209–17.
Colebrook EH, Thomas SG, Phillips AL, Hedden P. The role of gibberellin signalling in plant responses to abiotic stress. J Exp Biol. 2014;217:67–75.
Hu Y, Jia W, Wang J, Zhang Y, Yang L, Lin Z. Transgenic tall fescue containing the Agrobacterium tumefaciens ipt gene shows enhanced cold tolerance. Plant Cell Rep. 2005;23:705–9.
Tan D-X, Hardeland R, Manchester LC, Korkmaz A, Ma S, Rosales-Corral S, Reiter RJ. Functional roles of melatonin in plants, and perspectives in nutritional and agricultural science. J Exp Bot. 2012;63:577–97.
Shen C, Li D, He R, Fang Z, Xia Y, Gao J, Shen H, Cao M. Comparative transcriptome analysis of RNA-seq data for cold-tolerant and cold-sensitive rice genotypes under cold stress. J Plant Biol. 2014;57:337–48.
Song K, Kim HC, Shin S, Kim K-H, Moon JC, Kim JY, Lee BM. Transcriptome analysis of flowering time genes under drought stress in maize leaves. Front Plant Sci. 2017;8:267.
Chen H, Chen X, Chen D, Li J, Zhang Y, Wang A. A comparison of the low temperature transcriptomes of two tomato genotypes that differ in freezing tolerance: Solanum lycopersicum and Solanum habrochaites. BMC Plant Biol. 2015;15:132.
Zhu M, Meng X, Cai J, Li G, Dong T, Li Z. Basic leucine zipper transcription factor SlbZIP1 mediates salt and drought stress tolerance in tomato. BMC Plant Biol. 2018;18:83.
Kaur G, Asthir B. Molecular responses to drought stress in plants. Biol Plant. 2017;61:201–9.
Chen H, Lai Z, Shi J, Xiao Y, Chen Z, Xu X. Roles of Arabidopsis WRKY18, WRKY40 and WRKY60 transcription factors in plant responses to abscisic acid and abiotic stress. BMC Plant Biol. 2010;10:1–15.
Chen L, Song Y, Li S, Zhang L, Zou C, Yu D. The role of WRKY transcription factors in plant abiotic stresses. BBA - Gene Regulatory Mechanisms. 1819;2012:120–8.
Zhu JK. Salt and drought stress signal transduction in plants. Annu Rev Plant Biol. 2002;53:247–73.
Shinozaki K, Yamaguchi-shinozaki K, Seki M. Regulatory network of gene expression in the drought and cold stress responses. Curr Opin Plant Biol. 2003;6:410–7.
Zhou R, Kong L, Wu Z, Rosenqvist E, Wang Y, Zhao L, Zhao T, Ottosen CO. Physiological response of tomatoes at drought, heat and their combination followed by recovery. Physiol Plant. 2019;165:144–54.
Mittler R, Vanderauwera S, Suzuki N, Miller G, Tognetti VB, Vandepoele K, Gollery M, Shulaev V, Van Breusegem F. ROS signaling: the new wave? Trends Plant Sci. 2011;16:300–9.
Kreps JA, Wu Y, Chang H-S, Zhu T, Wang X, Harper J. Transcriptome changes for Arabidopsis in response to salt, osmotic, and cold stress. Plant Physiol. 2002;130:2129–41.
Seki M, Narusaka M, Ishida J, Nanjo T, Fujita M, Oono Y, Kamiya A, Nakajima M, Enju A, Sakurai K. Monitoring the expression profiles of 7000 Arabidopsis genes under drought, cold, and high-salinity stresses using a full-length cDNA microarray. Plant J. 2002;31:279–92.
Han Y, Li A, Li F, Zhao M, Wang W. Characterization of a wheat (Triticum aestivum L.) expansin gene, TaEXPB23, involved in the abiotic stress response and phytohormone regulation. Plant Physiol Biochem. 2012;54:49–58.
Li C, Tan DX, Liang D, Chang C, Jia D, Ma F. Melatonin mediates the regulation of ABA metabolism, free-radical scavenging, and stomatal behaviour in two Malus species under drought stress. J Exp Bot. 2015;66:669–80.
Kalapos B, Dobrev P, Nagy T, Vítámvás P, Györgyey J, Kocsy G, Marincs F, Galiba G. Transcript and hormone analyses reveal the involvement of ABA-signalling, hormone crosstalk and genotype-specific biological processes in cold-shock response in wheat. Plant Sci. 2016;253:86–97.
Yamaguchi-Shinozaki K, Shinozaki K. Transcriptional regulatory networks in cellular responses and toletance to dehydration and cold stresses. Annu Rev Plant Biol. 2006;57:781–803.
Baker NR, Rosenqvist E. Applications of chlorophyll fluorescence can improve crop production strategies: an examination of future possibilities. J Exp Bot. 2004;55:1607–21.
Deng B, Du W, Liu C, Sun W, Tian S, Dong H. Antioxidant response to drought, cold and nutrient stress in two ploidy levels of tobacco plants: low resource requirement confers polytolerance in polyploids? Plant Growth Regul. 2012;66:37–47.
Sales CRG, Ribeiro RV, Silveira JAG, Machado EC, Martins MO, Lagôa A.M. Superoxide dismutase and ascorbate peroxidase improve the recovery of photosynthesis in sugarcane plants subjected to water deficit and low substrate temperature. Plant Physiol Biochem 2013;73:326.
Rabbani MA, Maruyama K, Abe H, Khan MA, Katsura K, Ito Y, Yoshiwara K, Seki M, Shinozaki K, Yamaguchi-Shinozaki K. Monitoring expression profiles of rice genes under cold, drought, and high-salinity stresses and abscisic acid application using cDNA microarray and RNA gel-blot analyses. Plant Physiol. 2003;133:1755–67.
He Z. A laboratory guide to chemical control technology on field crop. In: He Z, editor. Beijing: Beijing Agricultural University Press; 1993. p. 60–8.
Yang YM, Xu CN, Wang BM, Jia JZ. Effects of plant growth regulators on secondary wall thickening of cotton fibres. Plant Growth Regul. 2001;35:233–7.
Murchie EH, Lawson T. Chlorophyll fluorescence analysis: a guide to good practice and understanding some new applications. J Exp Bot. 2013;64:3983–98.
The authors thank the assistance from Dr. Yinlei Wang, Ms. Liping Zhao and Mr. Duanyue Huang during the experiment.
The funders had no role in the experiment design, data analysis, decision to publish or preparation of the manuscript. The authors acknowledge the funding from Natural Science Foundation of China (31601745), Natural Science Foundation of Youth in Jiangsu Province (BK20160579) and Aarhus University Research Foundation (AUFF grant, 30379).
Ethics approval and consent to participate
The seeds and plant materials of tomato cultivar ‘Jinlingmeiyu’ used in this study were donated by Jiangsu Province Academy of Agricultural Sciences, Nanjing, China. The experimental research on plants comply with institutional, national, or international guidelines.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Figure S1. (A) Number of new genes predicted by de novo assembly, (B) distribution of Nr homologous species and (C) eggNOG function classification of consensus sequence. The Nr indicated non-redundant protein sequence database, the website of which was ftp://ftp.ncbi.nih.gov/blast/db/. The eggNOG indicated evolutionary genealogy of genes: Non-supervised Orthologous Groups, the website of which was http://eggnog.embl.de/. (JPG 158 kb)
Figure S2. The KEGG (Kyoto Encyclopedia of Genes and Genomes, http://www.genome.jp/kegg) pathway enrichment analysis of the DEGs (differentially expressed genes) between (A) control vs cold, (B) control vs drought, (C) control vs combined stress and (D) cold vs drought. (JPG 230 kb)
Table S1. Reads overview of the RNA-seq data in tomato at control, cold, drought and their combination. Table S2. The detailed information of the 12 differentially expressed genes (DEGs) for qRT-PCR validation. Table S3. Differentially expressed genes (DEGs) involved in reactive oxygen species (ROS) in tomato at control, cold, drought and their combination. Table S4. Differentially expressed genes (DEGs) involved in signal transduction in tomato at control, cold, drought and their combination. Table S5. Differentially expressed genes (DEGs) involved in phytohormone signal transduction in tomato at control, cold, drought and their combination. Table S6. Differentially expressed genes (DEGs) involved in photosynthetic electron transport in tomato at control, cold, drought and their combination. Table S7. Primers for qRT-PCR validation of differentially expressed genes (DEGs). (XLSX 714 kb)
About this article
Cite this article
Zhou, R., Yu, X., Zhao, T. et al. Physiological analysis and transcriptome sequencing reveal the effects of combined cold and drought on tomato leaf. BMC Plant Biol 19, 377 (2019). https://doi.org/10.1186/s12870-019-1982-9