Transcriptional responses of winter barley to cold indicate nucleosome remodelling as a specific feature of crown tissues

Integrated Transcriptomics and Metabolomics Analysis of the Fructan Metabolism Response to Low-Temperature Stress in Garlic

As the main reserve carbohydrate in garlic, fructan contributes to garlic's yield and quality formation. Numerous studies have shown that plant fructan metabolism induces a stress response to adverse environments. However, the transcriptional regulation mechanism of garlic fructan in low-temperature environments is still unknown. In this study, the fructan metabolism of garlic seedlings under low-temperature stress was revealed by transcriptome and metabolome approaches. With the extension of stress time, the number of differentially expressed genes and metabolites increased. Using weighted gene co-expression network analysis (WGCNA), three key enzyme genes related to fructan metabolism were screened (a total of 12 transcripts): sucrose: sucrose 1-fructosyltransferase (1-SST) gene; fructan: fructan 6G fructosyltransferase (6G-FFT) gene; and fructan 1-exohydrolase (1-FEH) gene. Finally, two hub genes were obtained, namely Cluster-4573.161559 (6G-FFT) and Cluster-4573.153574 (1-FEH). The correlation network and metabolic heat map analysis between fructan genes and carbohydrate metabolites indicate that the expression of key enzyme genes in fructan metabolism plays a positive promoting role in the fructan response to low temperatures in garlic. The number of genes associated with the key enzyme of fructan metabolism in trehalose 6-phosphate was the highest, and the accumulation of trehalose 6-phosphate content may mainly depend on the key enzyme genes of fructan metabolism rather than the enzyme genes in its own synthesis pathway. This study not only obtained the key genes of fructan metabolism in garlic seedlings responding to low temperatures but also preliminarily analyzed its regulatory mechanism, providing an important theoretical basis for further elucidating the cold resistance mechanism of garlic fructan metabolism.

Field plus lab experiments help identify freezing tolerance and associated genes in subtropical evergreen broadleaf trees: A case study of Camellia oleifera

The molecular mechanisms of freezing tolerance are unresolved in the perennial trees that can survive under much lower freezing temperatures than annual herbs. Since natural conditions involve many factors and temperature usually cannot be controlled, field experiments alone cannot directly identify the effects of freezing stress. Lab experiments are insufficient for trees to complete cold acclimation and cannot reflect natural freezing-stress responses. In this study, a new method was proposed using field plus lab experiments to identify freezing tolerance and associated genes in subtropical evergreen broadleaf trees using Camellia oleifera as a case. Cultivated C. oleifera is the dominant woody oil crop in China. Wild C. oleifera at the high-elevation site in Lu Mountain could survive below -30°C, providing a valuable genetic resource for the breeding of freezing tolerance. In the field experiment, air temperature was monitored from autumn to winter on wild C. oleifera at the high-elevation site in Lu Mountain. Leave samples were taken from wild C. oleifera before cold acclimation, during cold acclimation and under freezing temperature. Leaf transcriptome analyses indicated that the gene functions and expression patterns were very different during cold acclimation and under freezing temperature. In the lab experiments, leaves samples from wild C. oleifera after cold acclimation were placed under -10°C in climate chambers. A cultivated C. oleifera variety "Ganwu 1" was used as a control. According to relative conductivity changes of leaves, wild C. oleifera showed more freezing-tolerant than cultivated C. oleifera. Leaf transcriptome analyses showed that the gene expression patterns were very different between wild and cultivated C. oleifera in the lab experiment. Combing transcriptome results in both of the field and lab experiments, the common genes associated with freezing-stress responses were identified. Key genes of the flg22, Ca²⁺ and gibberellin signal transduction pathways and the lignin biosynthesis pathway may be involved in the freezing-stress responses. Most of the genes had the highest expression levels under freezing temperature in the field experiment and showed higher expression in wild C. oleifera with stronger freezing tolerance in the lab experiment. Our study may help identify freezing tolerance and underlying molecular mechanisms in trees.

Molecular breeding of barley for quality traits and resilience to climate change

Barley grains are a rich source of compounds, such as resistant starch, beta-glucans and anthocyanins, that can be explored in order to develop various products to support human health, while lignocellulose in straw can be optimised for feed in husbandry, bioconversion into bioethanol or as a starting material for new compounds. Existing natural variations of these compounds can be used to breed improved cultivars or integrated with a large number of mutant lines. The technical demands can be in opposition depending on barley's end use as feed or food or as a source of biofuel. For example beta-glucans are beneficial in human diets but can lead to issues in brewing and poultry feed. Barley breeders have taken action to integrate new technologies, such as induced mutations, transgenics, marker-assisted selection, genomic selection, site-directed mutagenesis and lastly machine learning, in order to improve quality traits. Although only a limited number of cultivars with new quality traits have so far reached the market, research has provided valuable knowledge and inspiration for future design and a combination of methodologies to achieve the desired traits. The changes in climate is expected to affect the quality of the harvested grain and it is already a challenge to mitigate the unpredictable seasonal and annual variations in temperature and precipitation under elevated [CO₂] by breeding. This paper presents the mutants and encoded proteins, with a particular focus on anthocyanins and lignocellulose, that have been identified and characterised in detail and can provide inspiration for continued breeding to achieve desired grain and straw qualities.

The maize low-lignin brown midrib3 mutant shows pleiotropic effects on photosynthetic and cell wall metabolisms in response to chilling

Maize (Zea mays L.) is one of the major cereal crops in the world and is highly sensitive to low temperature. Here, changes in photosynthetic and cell wall metabolisms were investigated during a long chilling exposure in inbred line F2 and a low-lignin near-isogenic brown midrib3 mutant (F2bm3), which has a mutation in the caffeic acid O-methyltransferase (COMT) gene. Results revealed that the plant biomass was reduced, and this was more pronounced in F2bm3. Photosynthesis was altered in both lines with distinct changes in photosynthetic pigment content between F2bm3 and F2, indicating an alternative photoprotection mechanism between lines under chilling. Starch remobilization was observed in F2bm3 while concentrations of sucrose, fructose and starch increased in F2, suggesting a reduced sugar partitioning in F2. The cell wall was altered upon chilling, resulting in changes in the composition of glucuronorabinoxylan and a reduced cellulose level in F2. Chilling shifted lignin subunit composition in F2bm3 mutant to a higher proportion of p-hydroxyphenyl (H) units, whereas it resulted in lignin with a higher proportion of syringyl (S) residues in F2. On average, the total cell wall ferulic acid (FA) content increased in both genotypes, with an increase in ether-linked FA in F2bm3, suggesting a greater degree of cross-linking to lignin. The reinforcement of the cell wall with lignin enriched in H-units and a higher concentration in cell-wall-bound FA observed in F2bm3 as a response to chilling, could be a strategy to protect the photosystems.

Metabolomics Reveals 5-Aminolevulinic Acid Improved the Ability of Tea Leaves (Camellia sinensis L.) against Cold Stress

Tea is an important woody crop whose cultivation is severely limited by cold stress. Although 5-aminolevulinic acid (ALA) is known to be effective in alleviating abiotic stresses in plants, knowledge of the detailed metabolic response of tea plants to exogenous ALA-induced cold resistance is still limited—a lack which restricts our ability to protect tea plants from cold stress. In the present study, we performed an in-depth metabolomics analysis to elucidate the metabolic responses of tea plants to cold stress and explore the role of ALA in improving tea plants’ cold-resistance capability. Metabolic profiles showed that cold stress altered various metabolisms in tea plants, especially galactose composition and flavonoid contents. Furthermore, exogenous ALA application altered a series of metabolisms associated with cold stress. Importantly, increases in metabolites, including catechin, 3,4-dihydroxyphenylacetic acid and procyanidin B2, involved in the mechanisms of ALA improved tea plants’ cold resistance. Overall, our study deciphered detailed metabolic responses of tea plants to cold stress and elucidated the mechanisms of ALA in enhancing cold resistance through rebuilding compositions of soluble carbohydrates and flavonoids. Therefore, we have provided a basis for exogenous usage of ALA to protect tea plants from cold stress.

High-throughput sequencing data revealed genotype-specific changes evoked by heat stress in crown tissue of barley sdw1 near-isogenic lines

BACKGROUND

High temperature shock is becoming increasingly common in our climate, affecting plant growth and productivity. The ability of a plant to survive stress is a complex phenomenon. One of the essential tissues for plant performance under various environmental stimuli is the crown. However, the molecular characterization of this region remains poorly investigated. Gibberellins play a fundamental role in whole-plant stature formation. This study identified plant stature modifications and crown-specific transcriptome re-modeling in gibberellin-deficient barley sdw1.a (BW827) and sdw1.d (BW828) mutants exposed to increased temperature.

RESULTS

The deletion around the sdw1 gene in BW827 was found to encompass at least 13 genes with primarily regulatory functions. A bigger genetic polymorphism of BW828 than of BW827 in relation to wild type was revealed. Transcriptome-wide sequencing (RNA-seq) revealed several differentially expressed genes involved in gibberellin metabolism and heat response located outside of introgression regions. It was found that HvGA20ox4, a paralogue of the HvGA20ox2 gene, was upregulated in BW828 relative to other genotypes, which manifested as basal internode elongation. The transcriptome response to elevated temperature differed in the crown of sdw1.a and sdw1.d mutants; it was most contrasting for HvHsf genes upregulated under elevated temperature in BW828, whereas those specific to BW827 were downregulated. In-depth examination of sdw1 mutants revealed also some differences in their phenotypes and physiology.

CONCLUSIONS

We concluded that despite the studied sdw1 mutants being genetically related, their heat response seemed to be genotype-specific and observed differences resulted from genetic background diversity rather than single gene mutation, multiple gene deletion, or allele-specific expression of the HvGA20ox2 gene. Differences in the expressional reaction of genes to heat in different sdw1 mutants, found to be independent of the polymorphism, could be further explained by in-depth studies of the regulatory factors acting in the studied system. Our findings are particularly important in genetic research area since molecular response of crown tissue has been marginally investigated, and can be useful for wide genetic research of crops since barley has become a model plant for them.

Plant Proteoforms Under Environmental Stress: Functional Proteins Arising From a Single Gene

Proteins are directly involved in plant phenotypic response to ever changing environmental conditions. The ability to produce multiple mature functional proteins, i.e., proteoforms, from a single gene sequence represents an efficient tool ensuring the diversification of protein biological functions underlying the diversity of plant phenotypic responses to environmental stresses. Basically, two major kinds of proteoforms can be distinguished: protein isoforms, i.e., alterations at protein sequence level arising from posttranscriptional modifications of a single pre-mRNA by alternative splicing or editing, and protein posttranslational modifications (PTMs), i.e., enzymatically catalyzed or spontaneous modifications of certain amino acid residues resulting in altered biological functions (or loss of biological functions, such as in non-functional proteins that raised as a product of spontaneous protein modification by reactive molecular species, RMS). Modulation of protein final sequences resulting in different protein isoforms as well as modulation of chemical properties of key amino acid residues by different PTMs (such as phosphorylation, N- and O-glycosylation, methylation, acylation, S-glutathionylation, ubiquitinylation, sumoylation, and modifications by RMS), thus, represents an efficient means to ensure the flexible modulation of protein biological functions in response to ever changing environmental conditions. The aim of this review is to provide a basic overview of the structural and functional diversity of proteoforms derived from a single gene in the context of plant evolutional adaptations underlying plant responses to the variability of environmental stresses, i.e., adverse cues mobilizing plant adaptive mechanisms to diminish their harmful effects.

Membrane-Enriched Proteomics Link Ribosome Accumulation and Proteome Reprogramming With Cold Acclimation in Barley Root Meristems

Due to their sessile nature, plants rely on root systems to mediate many biotic and abiotic cues. To overcome these challenges, the root proteome is shaped to specific responses. Proteome-wide reprogramming events are magnified in meristems due to their active protein production. Using meristems as a test system, here, we study the major rewiring that plants undergo during cold acclimation. We performed tandem mass tag-based bottom-up quantitative proteomics of two consecutive segments of barley seminal root apexes subjected to suboptimal temperatures. After comparing changes in total and ribosomal protein (RP) fraction-enriched contents with shifts in individual protein abundances, we report ribosome accumulation accompanied by an intricate translational reprogramming in the distal apex zone. Reprogramming ranges from increases in ribosome biogenesis to protein folding factors and suggests roles for cold-specific RP paralogs. Ribosome biogenesis is the largest cellular investment; thus, the vast accumulation of ribosomes and specific translation-related proteins during cold acclimation could imply a divergent ribosomal population that would lead to a proteome shift across the root. Consequently, beyond the translational reprogramming, we report a proteome rewiring. First, triggered protein accumulation includes spliceosome activity in the root tip and a ubiquitous upregulation of glutathione production and S-glutathionylation (S-GSH) assemblage machineries in both root zones. Second, triggered protein depletion includes intrinsically enriched proteins in the tip-adjacent zone, which comprise the plant immune system. In summary, ribosome and translation-related protein accumulation happens concomitantly to a proteome reprogramming in barley root meristems during cold acclimation. The cold-accumulated proteome is functionally implicated in feedbacking transcript to protein translation at both ends and could guide cold acclimation.

COR/LEA Proteins as Indicators of Frost Tolerance in Triticeae: A Comparison of Controlled versus Field Conditions

Low temperatures in the autumn induce enhanced expression/relative accumulation of several cold-inducible transcripts/proteins with protective functions from Late-embryogenesis-abundant (LEA) superfamily including dehydrins. Several studies dealing with plants grown under controlled conditions revealed a correlation (significant quantitative relationship) between dehydrin transcript/protein relative accumulation and plant frost tolerance. However, to apply these results in breeding, field experiments are necessary. The aim of the review is to provide a summary of the studies dealing with the relationships between plant acquired frost tolerance and COR/LEA transcripts/proteins relative accumulation in cereals grown in controlled and field conditions. The impacts of cold acclimation and vernalisation processes on the ability of winter-type Triticeae to accumulate COR/LEA proteins are discussed. The factors determining dehydrin relative accumulation under controlled cold acclimation treatments versus field trials during winter seasons are discussed. In conclusion, it can be stated that dehydrins could be used as suitable indicators of winter survival in field-grown winter cereals but only in plant prior to the fulfilment of vernalisation requirement.

ETHYLENE RESPONSE FACTOR39-MYB8 complex regulates low-temperature-induced lignification of loquat fruit

Flesh lignification is a specific chilling response that causes deterioration in the quality of stored red-fleshed loquat fruit (Eribotrya japonica) and is one aspect of wider chilling injury. APETALA2/ETHLENE RESPONSIVE FACTOR (AP2/ERF) transcription factors are important regulators of plant low-temperature responses and lignin biosynthesis. In this study, the expression and action of 27 AP2/ERF genes from the red-fleshed loquat cultivar 'Luoyangqing' were investigated in order to identify transcription factors regulating low-temperature-induced lignification. EjERF27, EjERF30, EjERF36, and EjERF39 were significantly induced by storage at 0 °C but inhibited by a low-temperature conditioning treatment (pre-storage at 5 °C for 6 days before storage at 0 °C, which reduces low-temperature-induced lignification), and their transcript levels positively correlated with flesh lignification. A dual-luciferase assay indicated that EjERF39 could transactivate the promoter of the lignin biosynthetic gene Ej4CL1, and an electrophoretic mobility shift assay confirmed that EjERF39 recognizes the DRE element in the promoter region of Ej4CL1. Furthermore, the combination of EjERF39 and the previously characterized EjMYB8 synergistically transactivated the Ej4CL1 promoter, and both transcription factors showed expression patterns correlated with lignification in postharvest treatments and red-fleshed 'Luoyangqing' and white-fleshed 'Ninghaibai' cultivars with different lignification responses. Bimolecular fluorescence complementation and luciferase complementation imaging assays confirmed direct protein-protein interaction between EjERF39 and EjMYB8. These data indicate that EjERF39 is a novel cold-responsive transcriptional activator of Ej4CL1 that forms a synergistic activator complex with EjMYB8 and contributes to loquat fruit lignification at low temperatures.

Physiological and transcriptomic analyses reveal the mechanisms underlying the salt tolerance of Zoysia japonica Steud

BACKGROUND

Areas with saline soils are sparsely populated and have fragile ecosystems, which severely restricts the sustainable development of local economies. Zoysia grasses are recognized as excellent warm-season turfgrasses worldwide, with high salt tolerance and superior growth in saline-alkali soils. However, the mechanism underlying the salt tolerance of Zoysia species remains unknown.

RESULTS

The phenotypic and physiological responses of two contrasting materials, Zoysia japonica Steud. Z004 (salt sensitive) and Z011 (salt tolerant) in response to salt stress were studied. The results show that Z011 was more salt tolerant than was Z004, with the former presenting greater K+/Na+ ratios in both its leaves and roots. To study the molecular mechanisms underlying salt tolerance further, we compared the transcriptomes of the two materials at different time points (0 h, 1 h, 24 h, and 72 h) and from different tissues (leaves and roots) under salt treatment. The 24-h time point and the roots might make significant contributions to the salt tolerance. Moreover, GO and KEGG analyses of different comparisons revealed that the key DEGs participating in the salt-stress response belonged to the hormone pathway, various TF families and the DUF family.

CONCLUSIONS

Zoysia salt treatment transcriptome shows the 24-h and roots may make significant contributions to the salt tolerance. The auxin signal transduction family, ABA signal transduction family, WRKY TF family and bHLH TF family may be the most important families in Zoysia salt-stress regulation.

Overexpression of the wheat trehalose 6-phosphate synthase 11 gene enhances cold tolerance in Arabidopsis thaliana

Low temperature is a key stress factor for the growth and development of wheat (Triticum aestivum L.), and glycometabolism plays an important role in plant cold tolerance. Our previous study identified trehalose 6-phosphate synthase 11 gene (TaTPS11), which had a significantly different expression pattern between a high freezing-tolerant wheat cultivar and a low freezing-tolerant wheat cultivar. In this study, TaTPS11 was isolated from a winter-hardy wheat cultivar (D1) and overexpressed in Arabidopsis thaliana to study its effect on cold tolerance in plants. Transgenic plants expressing TaTPS11 had lower sucrose content, higher starch content, and higher activity of key enzyme (sucrose phosphate synthase, sucrose synthase, and invertase) involved in sucrose metabolism. In addition, the expression level of sucrose non-fermenting 1-related kinase 1 (SnRK1), which catalyzes the sucrose in plants, increased in the TaTPS11-overexpressed plants. These results indicated that heterologous expression of TaTPS11 influenced carbohydrate metabolism in Arabidopsis plants. The resultant plants had a significantly higher survival rate after -5 °C treatment for 2 h and exhibited enhanced cold tolerance without unfavorable phenotypes compared to wild-type. Our findings indicated that manipulation of TaTPS11 improved cold tolerance in plants and TaTPS11 had potential values in wheat cold-tolerance breeding.

On a Cold Night: Transcriptomics of Grapevine Flower Unveils Signal Transduction and Impacted Metabolism

Low temperature is a critical environmental factor limiting plant productivity, especially in northern vineyards. To clarify the impact of this stress on grapevine flower, we used the Vitis array based on Roche-NimbleGen technology to investigate the gene expression of flowers submitted to a cold night. Our objectives were to identify modifications in the transcript levels after stress and during recovery. Consequently, our results confirmed some mechanisms known in grapes or other plants in response to cold stress, notably, (1) the pivotal role of calcium/calmodulin-mediated signaling; (2) the over-expression of sugar transporters and some genes involved in plant defense (especially in carbon metabolism), and (3) the down-regulation of genes encoding galactinol synthase (GOLS), pectate lyases, or polygalacturonases. We also identified some mechanisms not yet known to be involved in the response to cold stress, i.e., (1) the up-regulation of genes encoding G-type lectin S-receptor-like serine threonine-protein kinase, pathogen recognition receptor (PRR5), or heat-shock factors among others; (2) the down-regulation of Myeloblastosis (MYB)-related transcription factors and the Constans-like zinc finger family; and (3) the down-regulation of some genes encoding Pathogen-Related (PR)-proteins. Taken together, our results revealed interesting features and potentially valuable traits associated with stress responses in the grapevine flower. From a long-term perspective, our study provides useful starting points for future investigation.

Cytokinin at the Crossroads of Abiotic Stress Signalling Pathways

Cytokinin is a multifaceted plant hormone that plays major roles not only in diverse plant growth and development processes, but also stress responses. We summarize knowledge of the roles of its metabolism, transport, and signalling in responses to changes in levels of both macronutrients (nitrogen, phosphorus, potassium, sulphur) and micronutrients (boron, iron, silicon, selenium). We comment on cytokinin's effects on plants' xenobiotic resistance, and its interactions with light, temperature, drought, and salinity signals. Further, we have compiled a list of abiotic stress-related genes and demonstrate that their expression patterns overlap with those of cytokinin metabolism and signalling genes.

Heterologous Expression of a Novel Zoysia japonica C2H2 Zinc Finger Gene, ZjZFN1, Improved Salt Tolerance in Arabidopsis

Growing evidence indicates that some grass species are more tolerant to various abiotic and biotic stresses than many crops. Zinc finger proteins play important roles in plant abiotic and biotic stresses. Although genes coding for these proteins have been cloned and identified in various plants, their function and underlying transcriptional mechanisms in the halophyte Zoysia japonica are barely known. In the present study, ZjZFN1 was isolated from Z. japonica using RACE method. Quantitative real time PCR results revealed that the expression of ZjZFN1 was much higher in leaf than in root and stem tissues, and induced by salt, cold or ABA treatment. The subcellular localization assay demonstrated that ZjZFN1 was localized to the nucleus. Expression of the ZjZFN1 in Arabidopsis thaliana improved seed germination and enhanced plant adaption to salinity stress with improved percentage of green cotyledons and growth status under salinity stress. Physiological and transcriptional analyses suggested that ZjZFN1 might, at least in part, influence reactive oxygen species accumulation and regulate the transcription of salinity responsive genes. Furthermore, RNA-sequencing analysis of ZjZFN1-overexpressing plants revealed that ZjZFN1 may serve as a transcriptional activator in the regulation of stress responsive pathways, including phenylalanine metabolism, α-linolenic acid metabolism and phenylpropanoid biosynthesis pathways. Taken together, these results provide evidence that ZjZFN1 is a potential key player in plants' tolerance to salt stress, and it could be a valuable gene in Z. japonica breeding projects.

Transcriptome analysis of Brassica juncea var. tumida Tsen responses to Plasmodiophora brassicae primed by the biocontrol strain Zhihengliuella aestuarii

Mustard clubroot, caused by Plasmodiophora brassicae, is a serious disease that affects Brassica juncea var. tumida Tsen, a mustard plant that is the raw material for a traditional fermented food manufactured in Chongqing, China. In our laboratory, we screened the antagonistic bacteria Zhihengliuella aestuarii against P. brassicae. To better understand the biocontrol mechanism, three transcriptome analyses of B. juncea var. tumida Tsen were conducted using Illumina HiSeq 4000, one from B. juncea only inoculated with P. brassicae (P), one inoculated with P. brassica and the biocontrol agent Z. aestuarii at the same time (P + B), and the other was the control (H), in which P. brassicae was replaced by sterile water. A total of 19.94 Gb was generated by Illumina HiSeq sequencing. The sequence data were de novo assembled, and 107,617 unigenes were obtained. In total, 5629 differentially expressed genes between biocontrol-treated (P + B) and infected (P) samples were assigned to 126 KEGG pathways. Using multiple testing corrections, 20 pathways were significantly enriched with Qvalue ≤ 0.05. The resistance-related genes, involved in the production of pathogenesis-related proteins, pathogen-associated molecular pattern-triggered immunity, and effector-triggered immunity signaling pathways, calcium influx, salicylic acid pathway, reactive oxygen intermediates, and mitogen-activated protein kinase cascades, and cell wall modification, were obtained. The various defense responses induced by the biocontrol strain combatted the P. brassicae infection. The genes and pathways involved in plant resistance were induced by a biocontrol strain. The transcriptome data explained the molecular mechanism of the potential biocontrol strain against P. brassicae. The data will also serve as an important public information platform to study B. juncea var. tumida Tsen and will be useful for breeding mustard plants resistant to P. brassicae.

Plant Abiotic Stress Proteomics: The Major Factors Determining Alterations in Cellular Proteome

HIGHLIGHTS: Major environmental and genetic factors determining stress-related protein abundance are discussed.Major aspects of protein biological function including protein isoforms and PTMs, cellular localization and protein interactions are discussed.Functional diversity of protein isoforms and PTMs is discussed. Abiotic stresses reveal profound impacts on plant proteomes including alterations in protein relative abundance, cellular localization, post-transcriptional and post-translational modifications (PTMs), protein interactions with other protein partners, and, finally, protein biological functions. The main aim of the present review is to discuss the major factors determining stress-related protein accumulation and their final biological functions. A dynamics of stress response including stress acclimation to altered ambient conditions and recovery after the stress treatment is discussed. The results of proteomic studies aimed at a comparison of stress response in plant genotypes differing in stress adaptability reveal constitutively enhanced levels of several stress-related proteins (protective proteins, chaperones, ROS scavenging- and detoxification-related enzymes) in the tolerant genotypes with respect to the susceptible ones. Tolerant genotypes can efficiently adjust energy metabolism to enhanced needs during stress acclimation. Stress tolerance vs. stress susceptibility are relative terms which can reflect different stress-coping strategies depending on the given stress treatment. The role of differential protein isoforms and PTMs with respect to their biological functions in different physiological constraints (cellular compartments and interacting partners) is discussed. The importance of protein functional studies following high-throughput proteome analyses is presented in a broader context of plant biology. In summary, the manuscript tries to provide an overview of the major factors which have to be considered when interpreting data from proteomic studies on stress-treated plants.

Integrated analysis of transcriptomic and metabolomic data reveals critical metabolic pathways involved in polyphenol biosynthesis in Nicotiana tabacum under chilling stress

Chilling stress increases the amount of polyphenols, especially lignin, which protects tobacco (Nicotiana tabacum L. cv. k326) from chilling stress. To clarify the molecular biosynthesis mechanism of the key representative compounds, specifically lignin, RNA sequencing and ultra-high pressure liquid chromatography coupled to quadrupole-time of flight mass spectrometry technologies were used to construct transcriptomic and metabolomic libraries from the leaves of tobacco plants subjected to normal (25°C) and chilling (4°C) temperature treatments. Transcriptomic libraries from the different samples were sequenced, generating more than 40million raw reads. Among nine samples, metabolomic analysis identified a total of 97 encoding enzymes that function in the key steps of pathways related to polyphenol biosynthesis, where 42 metabolites were also located. An integrated analysis of metabolic and transcriptomic data revealed that most of the intermediate metabolites and enzymes related to lignin biosynthesis were synthesised in the leaves under chilling stress, which suggests that the biosynthesis of lignin plays an important role in the response of tobacco leaves to cold temperatures. In addition, the cold insensitivity of chalcone synthase genes might be considered to be an important rate-limiting factor in the process of precursor substance flow to flavonoid biosynthesis under chilling stress. Furthermore, the upregulated expression of phenylalanine ammonia lyase (PAL), hydroxycinnamoyl transferase (HCT) and cinnamyl-alcohol dehydrogenase (CAD) under chilling stress is the key to an increase in lignin synthesis. This study provides a hypothetical basis for the screening of new active metabolites and the metabolic engineering of polyphenols in tobacco.

Differential Proteomic Analysis Reveals the Effect of Calcium on Malus baccata Borkh. Leaves under Temperature Stress

In the cool apple-producing areas of northern China, air temperature during early spring changes in a rapid and dramatic manner, which affects the growth and development of apple trees at the early stage of the growing season. Previous studies have shown that the treatment of calcium can increase the cold tolerance of Malus baccata Borkh., a widely-used rootstock apple tree in northern China. To better understand the physiological function of calcium in the response of M. baccata to temperature stress, we analyzed the effect of calcium treatment (2% CaCl₂) on M. baccata leaves under temperature stress. Physiological analysis showed that temperature stress aggravated membrane lipid peroxidation, reduced chlorophyll content and induced photo-inhibition in leaves, whereas these indicators of stress injuries were alleviated by the application of calcium. An isobaric tags for relative and absolute quantitation (iTRAQ)-based proteomics approach was used in this study. Among the 2114 proteins that were detected in M. baccata leaves, 41, 25, and 34 proteins were differentially regulated by the increasing, decreasing, and changing temperature treatments, respectively. Calcium treatment induced 9 and 15 proteins after increasing and decreasing temperature, respectively, in comparison with non-treated plants. These calcium-responsive proteins were mainly related to catalytic activity, binding, and structural molecule activity. Hierarchical cluster analysis indicated that the changes in abundance of the proteins under increasing temperature and changing temperature treatments were similar, and the changes in protein abundance under decreasing temperature and increasing temperature with calcium treatment were similar. The findings of this study will allow a better understanding of the mechanisms underlying the role of calcium in M. baccata leaves under temperature stress.

Discovery of novel cold-induced CISP genes encoding small RNA-binding proteins related to cold adaptation in barley

To adapt to cold conditions, barley plants rely on specific mechanisms, which have not been fully understood. In this study, we characterized a novel barley cold-induced gene identified using a PCR-based high coverage gene expression profiling method. The identified gene encodes a small protein that we named CISP1 (Cold-induced Small Protein 1). Homology searches of sequence databases revealed that CISP1 homologs (CISP2 and CISP3) exist in barley genome. Further database analyses showed that the CISP1 homologs were widely distributed in cold-tolerant plants such as wheat and rye. Quantitative reverse transcription PCR analyses indicated that the expression of barley CISP genes was markedly increased in roots exposed to cold conditions. In situ hybridization analyses showed that the CISP1 transcripts were localized in the root tip and lateral root primordium. We also demonstrated that the CISP1 protein bound to RNA. Taken together, these findings indicate that CISP1 and its homologs encoding small RNA-binding proteins may serve as RNA chaperones playing a vital role in the cold adaptation of barley root. This is the first report describing the likely close relationship between root-specific genes and the cold adaptation process, as well as the potential function of the identified genes.

Global Scale Transcriptional Profiling of Two Contrasting Barley Genotypes Exposed to Moderate Drought Conditions: Contribution of Leaves and Crowns to Water Shortage Coping Strategies

Drought is a serious threat for sustainable agriculture. Barley represents a species well adapted to environmental stresses including drought. To elucidate the adaptive mechanism of barley on transcriptional level we evaluated transcriptomic changes of two contrasting barley cultivars upon drought using the microarray technique on the level of leaves and crowns. Using bioinformatic tools, differentially expressed genes in treated vs. non-treated plants were identified. Both genotypes revealed tissue dehydration under drought conditions as shown at water saturation deficit and osmotic potential data; however, dehydration was more severe in Amulet than in drought-resistant Tadmor under the same ambient conditions. Performed analysis showed that Amulet enhanced expression of genes related to active plant growth and development, while Tadmor regarding the stimulated genes revealed conservative, water saving strategy. Common reactions of both genotypes and tissues included an induction of genes encoding several stress-responsive signaling proteins, transcription factors as well as effector genes encoding proteins directly involved in stress acclimation. In leaf, tolerant cultivar effectively stimulated mainly the expression of genes encoding proteins and enzymes involved in protein folding, sulfur metabolism, ROS detoxification or lipid biosynthesis and transport. The crown specific reaction of tolerant cultivar was an enhanced expression of genes encoding proteins and enzymes involved in cell wall lignification, ABRE-dependent abscisic acid (ABA) signaling, nucleosome remodeling, along with genes for numerous jasmonate induced proteins.

Proteomic Response of Hordeum vulgare cv. Tadmor and Hordeum marinum to Salinity Stress: Similarities and Differences between a Glycophyte and a Halophyte

Response to a high salinity treatment of 300 mM NaCl was studied in a cultivated barley Hordeum vulgare Syrian cultivar Tadmor and in a halophytic wild barley H. marinum. Differential salinity tolerance of H. marinum and H. vulgare is underlied by qualitative and quantitative differences in proteins involved in a variety of biological processes. The major aim was to identify proteins underlying differential salinity tolerance between the two barley species. Analyses of plant water content, osmotic potential and accumulation of proline and dehydrin proteins under high salinity revealed a relatively higher water saturation deficit in H. marinum than in H. vulgare while H. vulgare had lower osmotic potential corresponding with high levels of proline and dehydrins. Analysis of proteins soluble upon boiling isolated from control and salt-treated crown tissues revealed similarities as well as differences between H. marinum and H. vulgare. The similar salinity responses of both barley species lie in enhanced levels of stress-protective proteins such as defense-related proteins from late-embryogenesis abundant family, several chaperones from heat shock protein family, and others such as GrpE. However, there have also been found significant differences between H. marinum and H. vulgare salinity response indicating an active stress acclimation in H. marinum while stress damage in H. vulgare. An active acclimation to high salinity in H. marinum is underlined by enhanced levels of several stress-responsive transcription factors from basic leucine zipper and nascent polypeptide-associated complex families. In salt-treated H. marinum, enhanced levels of proteins involved in energy metabolism such as glycolysis, ATP metabolism, and photosynthesis-related proteins indicate an active acclimation to enhanced energy requirements during an establishment of novel plant homeostasis. In contrast, changes at proteome level in salt-treated H. vulgare indicate plant tissue damage as revealed by enhanced levels of proteins involved in proteasome-dependent protein degradation and proteins related to apoptosis. The results of proteomic analysis clearly indicate differential responses to high salinity and provide more profound insight into biological mechanisms underlying salinity response between two barley species with contrasting salinity tolerance.

24-Epibrassinoslide enhances plant tolerance to stress from low temperatures and poor light intensities in tomato (Lycopersicon esculentum Mill.)

Brassinosteroids (Brs) are a newly recognized group of active steroidal hormones that occur at low concentrations in all plant parts and one of the active and stable forms is 24-epibrassinolide (EBR). We investigated the effect of EBR on tomato (Lycopersicon esculentum Mill.) and its mechanism when seedlings were exposed to low temperature and poor light stress conditions. Leaves of stress-tolerant 'Zhongza9' and stress-sensitive 'Zhongshu4' cultivars were pre-treated with spray solutions containing either 0.1 μM EBR or no EBR (control). The plants were then transferred to chambers where they were exposed to low temperatures of 12 °C/6 °C (day/night) under a low light (LL) level of 80 μmol · m(-2) · s(-1). Exogenous application of EBR significantly increased the antioxidant activity of superoxide dismutase, catalase and peroxidase, and decreased the rate of O2 · (-) formation and H2O2 and malondialdehyde contents. Additionally, the ATP synthase β subunit content was increased by exogenous hormone application. Based on these results, we conclude that exogenous EBR can elicit synergism between the antioxidant enzyme systems and the ATP synthase β subunit so that scavenging of reactive oxygen species becomes more efficient. These activities enable plants to cope better under combined low temperature and poor light stresses.

Low temperature tolerance in plants: Changes at the protein level

Low temperature (LT) is one of several important environmental stresses influencing plant performance and distribution. Adaptation to LT is a highly dynamic stress-response phenomenon and involves complex cross-talk between different regulatory levels. Although plants differ in their sensitivity to LT, in temperate species low nonfreezing temperatures cause noticeable alterations in various biochemical and physiological processes that can potentially improve freezing tolerance. This adaptation is associated with changes in the expression pattern of genes and their protein products. Proteins are the major players in most cellular events and are directly involved in plant LT responses, thereby proteome analysis could help uncover additional novel proteins associated with LT tolerance. Proteomics is recommended as an appropriate strategy for complementing transcriptome level changes and characterizing translational and post-translational regulations. In this review, we considered alterations in the expression and accumulation of proteins in response to LT stress in the three major cereal crops produced worldwide (wheat, barley, and rice). LT stress down-regulates many photosynthesis-related proteins. On the contrary, pathways/protein sets that are up-regulated by LT include carbohydrate metabolism (ATP formation), ROS scavenging, redox adjustment, cell wall remodelling, cytoskeletal rearrangements, cryoprotection, defence/detoxification. These modifications are common adaptation reactions also observed in the plant model Arabidopsis, thus representing key potential biomarkers and critical intervention points for improving LT tolerance of crop plants in cold regions with short summers. We believe that an assessment of the proteome within a broad time frame and during the different phenological stages may disclose the molecular mechanisms related to the developmental regulation of LT tolerance and facilitate the progress of genetically engineered stress-resistant plant varieties.

Comparative Transcriptomics of Sijung and Jumli Marshi Rice during Early Chilling Stress Imply Multiple Protective Mechanisms

INTRODUCTION

Low temperature is one of the major environmental factors that adversely affect plant growth and yield. Many cereal crops from tropical regions, such as rice, are chilling sensitive and, therefore, are affected already at <10 °C. Interestingly, it has been demonstrated that chilling susceptibility varies greatly among rice varieties, which indicates differences in the underlying molecular responses. Understanding these differences is vital for continued development of rational breeding and transgenic strategies for more tolerant varieties. Thus, in this study, we conducted a comparative global gene expression profiling analysis of the chilling tolerant varieties Sijung and Jumli Marshi (spp. Japonica) during early chilling stress (<24 h, 10 °C).

METHODS AND RESULTS

Global gene expression experiments were conducted with Agilent Rice Gene Expression Microarray 4 x 44 K. The analysed results showed that there was a relatively low (percentage or number) overlap in differentially expressed genes in the two varieties and that substantially more genes were up-regulated in Jumli Marshi than in Sijung but the number of down-regulated genes were higher in Sijung. In broad GO annotation terms, the activated response pathways in Sijung and Jumli Marshi were coherent, as a majority of the genes belonged to the catalytic, transcription regulator or transporter activity categories. However, a more detailed analysis revealed essential differences. For example, in Sijung, activation of calcium and phosphorylation signaling pathways, as well as of lipid transporters and exocytosis-related proteins take place very early in the stress response. Such responses can be coupled to processes aimed at strengthening the cell wall and plasma membrane against disruption. On the contrary, in Jumli Marshi, sugar production, detoxification, ROS scavenging, protection of chloroplast translation, and plausibly the activation of the jasmonic acid pathway were the very first response activities. These can instead be coupled to detoxification processes.

CONCLUSIONS

Based on the results inferred from this study, we conclude that different, but overlapping, strategies are undertaken by the two varieties to cope with the chilling stress; in Sijung the initial molecular responses seem to be mainly targeted at strengthening the cell wall and plasma membrane, whereas in Jumli Marshi the protection of chloroplast translation and detoxification is prioritized.

Transcriptome profiling of Vitis amurensis, an extremely cold-tolerant Chinese wild Vitis species, reveals candidate genes and events that potentially connected to cold stress

Vitis amurensis Rupr. is an exceptional wild-growing Vitis (grape) species that can safely survive a wide range of cold conditions, but the underlying cold-adaptive mechanism associated with gene regulation is poorly investigated. We have analyzed the physiochemical and transcriptomic changes caused by cold stress in a cold-tolerant accession, 'Heilongjiang seedling', of Chinese wild V. amurensis. We statistically determined that a total of 6,850 cold-regulated transcripts were involved in cold regulation, including 3,676 up-regulated and 3,174 down-regulated transcripts. A global survey of messenger RNA revealed that skipped exon is the most prevalent form of alternative spicing event. Importantly, we found that the total splicing events increased with the prolonged cold stress. We also identified thirty-eight major TF families that were involved in cold regulation, some of which were previously unknown. Moreover, a large number of candidate pathways for the metabolism or biosynthesis of secondary metabolites were found to be regulated by cold, which is of potential importance in coordinating cold tolerance with growth and development. Several heat shock proteins and heat shock factors were also detected to be intensively cold-regulated. Furthermore, we validated the expression profiles of 16 candidates using qRT-PCR to further confirm the accuracy of the RNA-seq data. Our results provide a genome-wide view of the dynamic changes in the transcriptome of V. amurensis, in which it is evident that various structural and regulatory genes are crucial for cold tolerance/adaptation. Moreover, our robust dataset advances our knowledge of the genes involved in the complex regulatory networks of cold stress and leads to a better understanding of cold tolerance mechanisms in this extremely cold-tolerant Vitis species.

The up-regulation of elongation factors in the barley leaf and the down-regulation of nucleosome assembly genes in the crown are both associated with the expression of frost tolerance

We report a series of microarray-based leaf and crown transcriptome comparisons involving three barley cultivars (cvs. Luxor, Igri and Atlas 68) which express differing degrees of frost tolerance. The transcripts were obtained following the exposure of seedlings to low (above and below zero) temperatures, aiming to identify those genes and signalling/metabolic pathways which are associated with frost tolerance. Both the leaves and the crowns responded to low temperature by the up-regulation of a suite of abscisic acid (ABA)-responsive genes, most of which have already been recognized as components of the plant low temperature response. The inter-cultivar comparison indicated that genes involved in maintaining the leaf's capacity to synthesize protein and to retain chloroplast activity were important for the expression of frost tolerance. In the crown, the repression of genes associated with nucleosome assembly and transposon regulation were the most relevant transcriptional changes associated with frost tolerance, highlighting the role of gene repression in the cold acclimation response.

Analysis of global gene expression in Brachypodium distachyon reveals extensive network plasticity in response to abiotic stress

Brachypodium distachyon is a close relative of many important cereal crops. Abiotic stress tolerance has a significant impact on productivity of agriculturally important food and feedstock crops. Analysis of the transcriptome of Brachypodium after chilling, high-salinity, drought, and heat stresses revealed diverse differential expression of many transcripts. Weighted Gene Co-Expression Network Analysis revealed 22 distinct gene modules with specific profiles of expression under each stress. Promoter analysis implicated short DNA sequences directly upstream of module members in the regulation of 21 of 22 modules. Functional analysis of module members revealed enrichment in functional terms for 10 of 22 network modules. Analysis of condition-specific correlations between differentially expressed gene pairs revealed extensive plasticity in the expression relationships of gene pairs. Photosynthesis, cell cycle, and cell wall expression modules were down-regulated by all abiotic stresses. Modules which were up-regulated by each abiotic stress fell into diverse and unique gene ontology GO categories. This study provides genomics resources and improves our understanding of abiotic stress responses of Brachypodium.

Global expression profiling of low temperature induced genes in the chilling tolerant japonica rice Jumli Marshi

Low temperature is a key factor that limits growth and productivity of many important agronomical crops worldwide. Rice (Oryza sativa L.) is negatively affected already at temperatures below +10°C and is therefore denoted as chilling sensitive. However, chilling tolerant rice cultivars exist and can be commercially cultivated at altitudes up to 3,050 meters with temperatures reaching as low as +4°C. In this work, the global transcriptional response to cold stress (+4°C) was studied in the Nepalese highland variety Jumli Marshi (spp. japonica) and 4,636 genes were identified as significantly differentially expressed within 24 hours of cold stress. Comparison with previously published microarray data from one chilling tolerant and two sensitive rice cultivars identified 182 genes differentially expressed (DE) upon cold stress in all four rice cultivars and 511 genes DE only in the chilling tolerant rice. Promoter analysis of the 182 genes suggests a complex cross-talk between ABRE and CBF regulons. Promoter analysis of the 511 genes identified over-represented ABRE motifs but not DRE motifs, suggesting a role for ABA signaling in cold tolerance. Moreover, 2,101 genes were DE in Jumli Marshi alone. By chromosomal localization analysis, 473 of these cold responsive genes were located within 13 different QTLs previously identified as cold associated.

Adaptation to seasonality and the winter freeze

Flowering plants initially diversified during the Mesozoic era at least 140 million years ago in regions of the world where temperate seasonal environments were not encountered. Since then several cooling events resulted in the contraction of warm and wet environments and the establishment of novel temperate zones in both hemispheres. In response, less than half of modern angiosperm families have members that evolved specific adaptations to cold seasonal climates, including cold acclimation, freezing tolerance, endodormancy, and vernalization responsiveness. Despite compelling evidence for multiple independent origins, the level of genetic constraint on the evolution of adaptations to seasonal cold is not well understood. However, the recent increase in molecular genetic studies examining the response of model and crop species to seasonal cold offers new insight into the evolutionary lability of these traits. This insight has major implications for our understanding of complex trait evolution, and the potential role of local adaptation in response to past and future climate change. In this review, we discuss the biochemical, morphological, and developmental basis of adaptations to seasonal cold, and synthesize recent literature on the genetic basis of these traits in a phylogenomic context. We find evidence for multiple genetic links between distinct physiological responses to cold, possibly reinforcing the coordinated expression of these traits. Furthermore, repeated recruitment of the same or similar ancestral pathways suggests that land plants might be somewhat pre-adapted to dealing with temperature stress, perhaps making inducible cold traits relatively easy to evolve.

Proteins involved in distinct phases of cold hardening process in frost resistant winter barley (Hordeum vulgare L.) cv Luxor

Winter barley is an economically important cereal crop grown in higher latitudes and altitudes where low temperatures represent an important environmental constraint limiting crop productivity. In this study changes in proteome of leaves and crowns in a frost tolerant winter barley cv. Luxor in relation to short and long term periods of cold followed by a brief frost treatment were studied in order to disclose proteins responsible for the cold hardening process in distinct plant tissues. The mentioned changes have been monitored using two dimensional difference gel electrophoresis (2D-DIGE) with subsequent peptide-mapping protein identification. Regarding approximately 600–700 distinct protein spots detected on 2D gels, there has been found at least a two-fold change after exposure to low temperatures in about 10% of proteins in leaves and 13% of proteins in crowns. Protein and nitrogen metabolic processes have been influenced by low temperature to a similar extent in both tissues while catabolism, carbohydrate metabolism and proteins involved in stress response have been more affected in crowns than in leaves. The range of changes in protein abundance was generally higher in leaves and chloroplast proteins were frequently affected which suggests a priority to protect photosynthetic apparatus. Overall, our data proved existence of slightly different response strategies to low temperature stress in crowns and leaves, i.e., tissues with different biological role. Moreover, there have been found several proteins with large increase in accumulation, e.g., 33 kDa oxygen evolving protein of photosystem II in leaves and “enhanced disease susceptibility 1” in crowns; these proteins might have potential to indicate an enhanced level of frost tolerance in barley.

Cell wall compositional modifications of Miscanthus ecotypes in response to cold acclimation

Miscanthus, a potential energy crop grass, can be damaged by late frost when shoots emerge too early in the spring and during the first winter after planting. The effects of cold acclimation on cell wall composition were investigated in a frost-sensitive clone of Miscanthus x giganteus compared to frost-tolerant clone, Miscanthus sinensis August Feder, and an intermediate frost-tolerant clone, M. sinensis Goliath. Cellulose and lignin contents were higher in M. x giganteus than in the M. sinensis genotypes. In ambient temperature controls, each clone displayed different glucuronoarabinoxylan (GAX) contents and degree of arabinose substitution on the xylan backbone. During cold acclimation, an increase in (1→3),(1→4)-β-D-glucan content was observed in all genotypes. Uronic acid level increased in the frost sensitive genotype but decreased in the frost tolerant genotypes in response to cold. In all clones, major changes in cell wall composition were observed with modifications in phenylalanine ammonia-lyase (PAL) and cinnamyl alcohol dehydrogenase (CAD) activities in both non- and cold-acclimated experiments. A large increase in CAD activity under cold stress was displayed in each clone, but it was largest in the frost-tolerant clone, M. sinensis August Feder. The marked increase in PAL activity observed in the frost-tolerant clones under cold acclimation, suggests a reorientation of the products towards the phenylpropanoid pathway or aromatic synthesis. How changes in cell wall physical properties can impact frost tolerance is discussed.

Induction of DREB2A pathway with repression of E2F, jasmonic acid biosynthetic and photosynthesis pathways in cold acclimation-specific freeze-resistant wheat crown

Winter wheat lines can achieve cold acclimation (development of tolerance to freezing temperatures) and vernalization (delay in transition from vegetative to reproductive phase) in response to low non-freezing temperatures. To describe cold-acclimation-specific processes and pathways, we utilized cold acclimation transcriptomic data from two lines varying in freeze survival but not vernalization. These lines, designated freeze-resistant (FR) and freeze-susceptible (FS), were the source of crown tissue RNA. Well-annotated differentially expressed genes (p ≤ 0.005 and fold change ≥ 2 in response to 4 weeks cold acclimation) were used for gene ontology and pathway analysis. "Abiotic stimuli" was identified as the most enriched and unique for FR. Unique to FS was "cytoplasmic components." Pathway analysis revealed the "triacylglycerol degradation" pathway as significantly downregulated and common to both FR and FS. The most enriched of FR pathways was "neighbors of DREB2A," with the highest positive median fold change. The "13-LOX and 13-HPL" and the "E2F" pathways were enriched in FR only with a negative median fold change. The "jasmonic acid biosynthesis" pathway and four "photosynthetic-associated" pathways were enriched in both FR and FS but with a more negative median fold change in FR than in FS. A pathway unique to FS was "binding partners of LHCA1," which was enriched only in FS with a significant negative median fold change. We propose that the DREB2A, E2F, jasmonic acid biosynthesis, and photosynthetic pathways are critical for discrimination between cold-acclimated lines varying in freeze survival.

Differential modulation of photosynthesis, signaling, and transcriptional regulation between tolerant and sensitive tomato genotypes under cold stress

The wild species Solanum habrochaites is more cold tolerant than the cultivated tomato (S. lycopersicum). To explore the mechanisms underlying cold tolerance of S. habrochaites, seedlings of S. habrochaites LA1777 introgression lines (ILs), as well as the two parents, were evaluated under low temperature (4°C). The IL LA3969 and its donor parent LA1777 were found to be more cold tolerant than the recurrent parent S. lycopersicum LA4024. The differences in physiology and global gene expression between cold-tolerant (LA1777 and LA3969) and -sensitive (LA4024) genotypes under cold stress were further investigated. Comparative transcriptome analysis identified 1613, 1456, and 1523 cold-responsive genes in LA1777, LA3969, and LA4024, respectively. Gene ontology (GO) term enrichment analysis revealed that more GO biological process terms were significantly enriched among the up-regulated genes in the two tolerant genotypes, whereas more biological processes were significantly repressed by cold stress in the sensitive one. A total of 92 genes with significant differential expression between tolerant and sensitive genotypes under cold stress were identified. Among these, many stress-related GO terms were significantly enriched, such as 'response to stimulus' and 'response to stress'. Moreover, GO terms 'response to hormone stimulus', 'response to reactive oxygen species (ROS)', and 'calcium-mediated signaling' were also overrepresented. Several transcripts involved in hormone or ROS homeostasis were also differentially expressed. ROS, hormones, and calcium as signaling molecules may play important roles in regulating gene expression in response to cold stress. Moreover, the expression of various transcription factors, post-translational proteins, metabolic enzymes, and photosynthesis-related genes was also specifically modulated. These specific modifications may play pivotal roles in conferring cold tolerance in tomato. These results not only provide new insights into the molecular mechanisms of cold tolerance in tomato, but also provide potential candidate genes for genetic improvement.