Skip to main content
  • Research article
  • Open access
  • Published:

Transcriptome analysis of postharvest blueberries (Vaccinium corymbosum ‘Duke’) in response to cold stress



Blueberry (Vaccinium spp.) is a small berry with high economic value. Although cold storage can extend the storage time of blueberry to more than 60 days, it leads to chilling injury (CI) displaying as pedicle pits; and the samples of 0 °C-30 days was the critical point of CI. However, little is known about the mechanism and the molecular basis response to cold stress in blueberry have not been explained definitely. To comprehensively reveal the CI mechanisms in response to cold stress, we performed high-throughput RNA Seq analysis to investigate the gene regulation network in 0d (control) and 30d chilled blueberry. At the same time, the pitting and decay rate, electrolyte leakage (EL), malondialdehyde (MDA) proline content and GSH content were measured.


Two cDNA libraries from 0d (control) and 30d chilled samples were constructed and sequenced, generating a total of 35,060 unigenes with an N50 length of 1348 bp. Of these, 1852 were differentially expressed, with 1167 upregulated and 685 downregulated. Forty-five cold-induced transcription factor (TF) families containing 1023 TFs were identified. The DEGs indicated biological processes such as stress responses; cell wall metabolism; abscisic acid, gibberellin, membrane lipid, energy metabolism, cellular components, and molecular functions were significantly responsed to cold storage. The transcriptional level of 40 DEGs were verified by qRT-PCR.


The postharvest cold storage leads serious CI in blueberry, which substantially decreases the quality, storability and consumer acceptance. The MDA content, proline content, EL increased and the GSH content decreased in this chilled process. The biological processes such as stress responses, hormone metabolic processes were significantly affected by CI. Overall, the results obtained here are valuable for preventing CI under cold storage and could help to perfect the lack of the genetic information of non-model plant species.


Both production and consumption of blueberry (Vaccinium spp.) have increased sharply worldwide in recent years at least partly due to their known nutritional, economic value and health benefits. However, it can only be stored for 5–10 days at room temperature (RT) due to the rapid softening. Since how to prevent the decay and prolonging the shelf-life became important topics in postharvest research. Storage at low temperature (LT) is advantageous to preventing softening and prolonging the postharvest life; however, cold storage triggers the pedicle pitting, pericarp and pulp adhesion in postharvest blueberry; which also accompanied with abnormal changes of stress-related enzymes, destruction of outer cell wall structure, and the decrease of fiber filaments [1, 2]. The development of this CI reduces consumer acceptance of the fruit, thus limiting its storage life. And it has been confirmed that cold stress accelerated the energy consumption, membrane lipid peroxidation and affected normal active oxygen metabolism in blueberry.

Chilling injury is one of the main problems affecting the market value of horticultural fruits and vegetables; which is also one of the principal abiotic stress. Cold stress at low temperatures above zero (0–10 °C) and freezing stress at subzero temperatures lead to decreased cell membrane fluidity, water potential, osmotic stress causing irreversible damage. In previous research, cold stress (4 °C) has been found to damage the cellular structure of peppers (Capsicum annuum L.) [3,4,5]; and the most typical CI symptom in peppers is water immersion sunken spots [6, 7]; in banana, the peel would pit and discolor at low temperatures [8]; in peach, CI symptoms include lack of juiciness [9], poor flavor, and the brown-red discoloration of flesh [10]; in mango, CI symptoms include scalding, softening, internal browning and electrolyte leakage [11]; in litchi, fruit quality, membrane permeability, enzyme activities and energy changed during cold storage [12]; in pear, the aroma was less under or after long term storage at 0 °C, seriously affecting its quality [13,14,15]; in harvested cucumber (Cucumis sativus L.) fruit, disorder characterized by surface pitting and dark watery patches was confirmed when when held at 7 °C or below [16]. In conclusion, it has been identified when the fruits are subjected to cold stress, the cell structure changed, the membranes disrupted, the malondialdehyde (MDA) content increased, and the reactive oxygen (ROS) accumulated [17, 18], finally led to the decline of fruit quality. Although the CI symptoms in different fruits and vegetables are distinct, the nutritional and economic values are adversely affected. And the mechanisms involved in the CI of blueberry are still not clear. Therefore there also remains a need to study the chilling mechanism and develop more effective techniques for blueberry fruit storage and transport at LT. Moreover, with the rapid development of high-throughput sequencing, great progress has been made towards the understanding the cold-response mechanism underlying CI.

Among these, RNA-sequencing (RNA-Seq) has become a powerful technology for characterizing molecular regulators in numerous postharvest fruits, such as peppers, tomatoes, and ‘Nanguo’ pear [14, 19]. Additionally, (Zhao et al. 2019) [20] identified four PbBAM genes associated with stress in pear using RNA-Seq; (Lou et al. 2018) [21] combined the isobaric tags for relative and absolute quantification (iTRAQ) and RNA-Seq illustrated the mechanism of cold tolerance in the loquat. Previous studies have also identified and confirmed the expression pattern of some key genes and transcription factor (TFs) are in response to cold stress. Among these, bZIP, MYB, CBF, and AP2/ERF are the most well-known TFs and are regarded as important regulators of cold-responsive genes in plants [22,23,24,25,26,27,28]. However, CI studies on postharvest blueberry are scarce and the genetic characterization of the underlying regulatory mechanisms of pitting in blueberry has rarely been investigated.

In the present study, in order to characterize the CI mechanism at the molecular level, enrich available transcriptome data, identify genes, important pathways and regulation network potentially involved in fruit response to cold storage, a comprehensive transcriptome profiling analysis and physiological experiments were performed on 0d (control) and 30d chilled blueberry (pitted). The results obtained here provide a theoretical and a practical basis for preventing blueberry pitting; it is also the basis for studying the functions of genes and transcription factors affecting the blueberry fruit pitting and provide a reference transcriptome for future studies involving blueberry fruit.


Changes in phenotype and physiology at low temperature

Changes in blueberry fruit phenotype at low temperature

The physiological performance of control blueberry and 30d chilled blueberry was characterized. The blueberry stored at 20 °C did not show CI symptoms (Fig. 1a), whereas the fruits stored at 0 °C for 30 days showed typical CI symptoms as pedicle pitting (Fig. 1b); these symptoms were maintained during the shelf-life.

Fig. 1
figure 1

The appearance of blueberry (a) without pitting, (b) with pitting, 30d chilled berries

Changes in pitting rate and decaying rate at low temperature

After the measurements of the pitting rate in blueberry at shelf-life after different storage times at 0 °C; we found no signs of pitting were observed in the blueberry stored directly at 20 °C or for those stored at 0 °C for 15 days beforehand (Fig. 2). However, signs of pitting were observed in blueberry stored at 0 °C for 30 days beforehand. The pitting rate was 5.2% after 30 days 0 °C storage, after then rapidly increased to 21.9% at 2d shelf-life, and then 28.6% at 4d, 35.8% at 6d and 45.2% at 8d. For blueberry stored at 0 °C for 45 days and 60 days, the pitting rates increased to 65.3 and 81% when the fruits were removed from 0 °C and kept at 8d shelf-life. The decay rate had a similar trend; both increased gradually, but the decay rate was slightly higher in long-term low temperature stored blueberry than that in the short-term low temperature stored (Fig. 2).

Fig. 2
figure 2

Changes of pitting and decay rates at shelf-life (20 °C) after different storage time at 0 °C. The lines indicate the pitting rate and the column diagram indicate the decay rate. Each value is the mean of three replicates of 150 fruits (P < 0.05). The error bars show the SD of the means (n = 3; three biological replicates, each with three technical replicates)

Changes in physiology at low temperature

Membrane plays a key role in fruit. Membrane lipid peroxidation often occurs during CI and MDA is one of the products of membrane lipid peroxidation. As shown in (Fig. 3a), there was an upward trend in the MDA content, and the MDA content of blueberry fruit stored at 0 °C for 30 days (9.2 mmol kg− 1) was higher than that of fruit stored at 20 °C (6.5 mmol kg− 1). Meanwhile, the permeability of the cell membrane can reflect the stability and the injury degree of the cell membrane; which can reflect by the cell membrane relative conductivity. It can be seen from (Fig. 3b) that with the prolongation of storage, the permeability of cell membrane of blueberry fruit increased continuously. After 30 days storage at 0 °C, membrane structure, phase change and permeability changed; the relative electrolytic leakage of the blueberry stored at 0 °C for 30 days (42%) was higher than that of blueberry stored at 20 °C (36%), then reached to 81%. Meanwhile, with the occurrence of pedicle pitting, the relative conductivity of fruit cell membrane increased sharply. The results showed that the GSH content in blueberry fruit stored at 0 °C for 30 days (89.2%) was higher than that of fruit stored at 20 °C (75%); the GSH content accumulated with the prolongation of time, especially in the fruits stored at 20 °C-6 days (92.4%) (Fig. 3c). After 30 days’ storage at 0 °C, the GSH content of blueberry fruit decreased rapidly from 89.2%, and then maintained a low stable state (64.6%). The proline content of blueberry fruit stored at 0 °C for 30 days (2.25 g kg− 1) was higher than that of fruit stored at 20 °C (0.6 g kg− 1) (Fig. 3d), and the proline content increased significantly after the fruits were removed from cold storage and reached 3.95 g kg− 1 at 8d shelf-life. The results indicated that the cell membrane structure and membrane function were seriously damaged during cold storage.

Fig. 3
figure 3

Changes in (a) MDA content, (b) relative electrolyte leakage (EL), (c) GSH content, (d) proline content at shelf-life and shelf-life after 30 days’ 0 °C storage. The error bars show the SD of the means (n = 3; three biological replicates, each with three technical replicates)

Summary of the sequencing data

To study the effect of cold storage on postharvest blueberries, total RNAs from six samples (control and chilled) were used for deep sequencing. To evaluate the quality of the sequencing data, randomization test of mRNA fragmentation saturation, length test of insert fragment and saturation test were conducted. These tests determined the number of obtained reads was enough to cover most of the expressed genes. The total 50.74Gb of clean data were generated (Table 1). Average clean data were above 6.25Gb and Q30 values were obtained for more than 93.93% of the data; 76.38 to 79.53% of the clean reads were successfully mapped to the unigene database. Meanwhile the transcription group data detection has a high sensitivity, with the protein encoding gene expression level FPKM (fragments per kilobase million) having values across 10− 2 to 10− 4 six orders of magnitude. The dispersion degree of sample gene expression level distribution was average, and the overall gene expression abundance in different samples was good (Fig. 4a). According to the PCA analysis, the three biological replicates of control and chilled groups were clustered together respectively (Fig. 4c), showing a Pearson correlation above 0.94 (Fig. 4b).

Table 1 Summary of RNA-Seq data and sequence assembly
Fig. 4
figure 4

(a) Box plot representation of expression range from six libraries; (b) Pearson correlation between samples and (c) PCA analysis among T01-T06 samples; Different colored squares represent the degree of correlation between two samples; T01-T03 represented 0d control samples and T04-T06 represented 30d-chilled samples

Comprehensive profiling of transcript expression analysis

Number of DEGs

The analysis of DEGs revealed 561 up-regulated and 158 down-regulated genes for FDR ≤ 0.01, and 1167 up-regulated and 685 down-regulated genes for FDR < 0.05. To assess the diversity of DEGs, the Venn map was constructed (Fig. 5).

Fig. 5
figure 5

Venn diagram of differentially expressed genes between each two samples

Heatmap of all DEGs and the selection of significantly changed DEGs

The identified DEGs were analyzed by hierarchical clustering to organize genes with the same or similar expression behaviors to show the different expression patterns of gene sets under different experimental conditions (Fig. 6). In the blueberry kept at 20 °C, most genes were significantly down-regulated; in blueberry kept at 0 °C for 30 days, most of these genes were up-regulated. We clustered these genes into 15 categories according to their functional pathways. These included GABA receptor activity, extracellular-glutamate-gated ion channel activity, and the G-protein coupled receptor signaling pathway. There were 30, 27, 37, 92, 14, 351, 217, 193, 315, 1, 161, 409, 2, 1, 2 genes in each cluster via heatmap analysis, and there were more up-regulated genes in the G4, G5, G6, G7, G8, and G12 clusters, while the down-regulated genes mostly belonged to the G1, G2, G3, G9, G10, G11, G13, G14 and G15 clusters (Fig. 6).

Fig. 6
figure 6

Thermal maps of DEGs and 15 clusters produced for thermographic analysis. Different columns represent different samples and different rows represent different genes. Each small square represents a gene and the color indicates the level of expression; red represents up-regulation and green represents down-regulation

Meantime, the 10 DEGs most significantly changed among treatment groups were selected to assess the important genes and pathways involved in CI (Table 2). Among these DEGs, c127608 was up-regulated and involved in six pathways and c120356 was down-regulated and involved in three pathways of the blueberry kept at 0 °C for 30 days.

Table 2 Genes significantly up-regulated and down-regulated in 30d-chilled blueberry

Functional annotations and classifications

In our research, the N50 was 1348 bp, and there were 21,649 unigenes longer than 1Kb. The assembly integrity was high among the 84,260 unigenes. To predict and analyze the functions of the assembled unigenes, the NR, Swiss-Prot, GO, COG, KOG, Pfam, KEGG, KOG, and eggNOG databases were used as the basic local alignment search tool: 1142 (18.4%) unigenes were annotated in NR; 834 (13.4%) in Swiss-Prot; 685 (11.2%) in GO; 394 (6.4%) unigenes in COG; 922 (14.8%) in Pfam; 413 (6.7%) in KEGG; 638 (10.4%) in KOG; and 1099 (17.7%) in eggNOG.

Then the GO enrichment analysis was performed to obtain functional information for the DEGs mentioned above. Significantly enriched GO terms, identified based on corrected P-value< 0.05, were found for DEGs involved in biological processes, cellular components, and molecular functions (Fig. 7a). The DEGs found in chilled blueberry were mainly enriched in biological processes, metabolic processes, cellular processes, single organism processes, responses to stimulus, biological regulation, localization, cellular components, developmental processes, and 12 other different biological processes. The DEGs included in the metabolic process, cellular process, and single-organism process categories were not expressed in blueberry that did not pit.

Fig. 7
figure 7

Functional annotations and classifications of unigenes in 30d-chilled blueberry according to (a) GO, (b) KEGG and (c) COG databases

To further identify the functions of DEGs in response to cold stress, the KEGG database was used to classify and characterize the DEGs into corresponding pathways. By comparing the samples, we found that 413 DEGs were enriched into 50 pathways; three of these 50 pathways were related to cellular processes, two pathways belonged to environmental information processing, 12 pathways contained 98 DEGs involved in genetic information processing, 32 pathways included 174 DEGs involved in metabolism, and the remaining pathway was related to organism systems (Fig. 7b). The results of the KEGG analysis showed that, in blueberry under cold stress, 26 DEGs were involved in plant hormone signal transduction pathways, including auxin (Aux), cytokinin (CK), gibberellin (GA), abscisic acid (ABA), ethylene (ET), brassinosteroid (BR), jasmonic acid (JA), and salicylic acid (SA) signaling pathways, each gathering different numbers of DEGs. Moreover, carotenoid biosynthesis, which is involved in the biosynthesis of ABA, comprised five DEGs. These results indicated that not only fatty acid metabolism, but also multiple plant hormone signal transduction pathways were activated under cold stress leading to blueberry fruit pitting. Based on enrichment results, the six most significantly changed pathways in response to cold stress were plant hormone transduction, carotenoid biosynthesis, GSH metabolism, starch and sucrose metabolism, protein processing in endoplasmic reticulum, and chlorophyll metabolism (Table 3).

Table 3 The top 6 enriched pathways of DEGs in 30d-chilled blueberry

In addition, 394 DEGs (25classes) were annotated in the COG database (Fig. 7c). The top three COG terms were: general function prediction only (95DEGs); replication, recombination and repair (50 DEGs). Energy production and conversion, amino acid transport and metabolism of carbohydrate transport, inorganic ion transport, and metabolic biosynthesis transport and catabolism were also important.

Identification and analysis of the DEGs under cold stress

Cold-response DEGs involved in membrane lipid and energy metabolism

In this study, six highly differentially expressed pathways of the 129 pathways including membrane lipid metabolism, including glycerolipid metabolism (Ko00561), glycerophospholipid metabolism (Ko00564), ether lipid metabolism (Ko00565), α-linolenic acid metabolism (Ko00592), sphingolipid metabolism (Ko00600), and phosphatidylinositol signaling system (Ko04070) were significantly changed in chilled blueberry. There was one down-regulated DEG in Ko00561; and the Ko00564; Ko00592, Ko00565, Ko00564, Ko00600, and Ko04070 comprised one, one, four, two, and one highly differentially up-regulated expressed structural genes, respectively. Among these six pathways, Ko00564 and Ko00592 were strongly induced in plants under cold stress; c107765.graphc0(Phosphocholine CytidylyltransferaseA) was down-regulated; the c100748.graphc0 (glycerol-3-phosphate dehydrogenase1), c131228.graphc1 (phospholipaseD1/2), c125240.graphc0 (N-myristoyltransferase), c126471.graphc0 (Triacylglycerol lipase4), and c130733.graphc0 (phosphatidylserine synthase2) were up-regulated in Ko00564, while c130733.graphc0 was up-regulated in Ko00592.

The changes in cell membrane lipid compositions are closely related to fatty acid metabolism under cold stress. Six of the 129 pathways significantly changed in chilled blueberry were related to fatty acid metabolism, which included fatty acid biosynthesis (Ko00061), fatty acid elongation (Ko00062), fatty acid degradation (Ko00071), linoleic acid metabolism (Ko00591), biosynthesis of unsaturated fatty acids (Ko01040), and fatty acid metabolism (Ko01212). There were three, one, and two DEGs down-regulated in Ko00061, Ko00071, and Ko01040, respectively; there were one, one, five, one, and two DEGs up-regulated in Ko00062, Ko00071, Ko00591, Ko01040, and Ko01212, respectively. Genes c126515.graphc0 (Fas-associated protein with death domain), c109637.graphc0 (fatty acyl-ACP thioesteraseB), c126471.graphc0 (doxycycline) were down-regulated in Ko00061 (Fig. 8a).

Fig. 8
figure 8

Transcripts involved in (a) liquid related pathways, (b) proline, glutathione, and flavonoid metabolism, (c) brassinosteroid biosynthesis, carotenoid biosynthesis, and zeatin biosynthesis in chilled blueberry

In the Ko00230 purine metabolism pathway, the DEGs related to energy metabolism affected by cold stress. Gene c121781.graphc0, which was up-regulated, participates in the pentose phosphate pathway. This is involved in several biosynthesis pathways of microbial metabolism in diverse environments and in carbon metabolism biosynthesis; c123543.graphc0 was expressed in chilled blueberry and its expression was lower in this group than control fruits. In the Ko00280 valine, leucine, and isoleucine metabolism pathway, c122233graphc0 was up-regulated; 3-Hydroxyisobutyryl-CoA Hydrolase was co-regulated by two genes, c120356.graph and c122037.graphc0, which were significantly down-regulated (fold-change = − 5.71) and are involved in the beta-alanine metabolic and carbon metabolism pathways.

Cold-response genes involved in proline, glutathione, and flavonoid metabolism

Additionally, dehydration stress is often accompanied by osmotic adjustment. Plants need to accumulate or reduce organic or inorganic substances such as proline, glutathione, and flavonoid to maintain cellular moisture holding capacity. According to pathway enrichment analysis, arginine and proline metabolism is activated by cold-induced dehydration. In our study, c126860.graphc0 (NADP-specific glutamate dehydrogenase) also participates in alanine, aspartate, glutamate metabolism, D-glutamate metabolism, and carbon metabolism, was up-regulated; c123304.graphc0 (polyamine oxidase2/3/4) was up-regulated during spermidine synthesis and it is also involved in beta-alanine metabolism. In our validation experiments, we found that the relative expression of c126860.graphc0 in chilled blueberry was significantly up-regulated. The expression of c123304.graphc0 was also significantly up-regulated in chilled blueberry.

In the present study, c123222.graphc0 [Isocitrate dehydrogenase1/2 (IDH1/2)] plays a role in intermediary metabolism and energy production, might tightly associate or interact with the pyruvate dehydrogenase complex in the citric acid cycle (TCA cycle). Some genes related to glutathione were also identified. GPX and glutathione transferase (GST) are also important scavengers of ROS that participate in many kinds of environmental stress responses to adverse conditions. Glutathione peroxidase participates in arachidonic acid metabolism and, together with c112252.graphc0 [glutathione s-transferase (GST)], c119253.graphc0 (GST), c102135.graphc0 (GST), c105671.graphc0 (GST), c105671.graphc1 (GST), c122443.graphc0 (GST), c110117.graphc0, and c120425.graphc0 coordinates the GST metabolism of xenobiotics by cytochrome P450. In our validation, GST and IDH1 were down-regulated in chilled blueberry compared to the control blueberry. The expression of GST1 was significantly higher in chilled blueberry; the expressions of GST2 and GST3 were significantly lower in chilled blueberry (Fig. 8b).

The gene c116187.graphc0 [Hydroxycinnamoyl-CoA shikimate (HCT)] encodes shikimate O-hydroxycinn amoyltransferase (EC:, which participates in metabolic pathways and secondary metabolite biosynthesis; the leucoanthocyanidin reductases (LAR) encoded by c119207.graphc1 and c110095.graphc0 participate in secondary metabolite biosynthesis. Our qRT-PCR results showed that the expressions of HCT and LAR (c110095.graphc0) were down-regulated and that of LAR (c119207.graphc1) was up-regulated in chilled blueberry. (Fig. 8b).

Cold-response genes involved in hormone biosynthesis and signal transduction

Brassinosteroid biosynthesis (Ko00905), carotenoid biosynthesis (Ko00906), and zeatin biosynthesis (Ko00908).

In the zeatin biosynthesis pathway, c111548.graphc0 (Cytokinin oxidase/dehydrogenase) was significantly down-regulated to − 4.5 in chilled blueberry; c126830.graphc0 (tRNA isopentenyltransferase1) was up-regulated during the production of uridine phosphorylase and cis-zeatin glucoside. Genes c112010.graphc0 [Cytochrome-P450-90A1-(CYP90A1)] and c115136.graphc0 [Cytochrome-P450-85A2-(CYP85A2)] participate in secondary metabolite biosynthesis and metabolic pathways, including BR biosynthesis. These two genes had higher expression in chilled blueberry (Fig. 8c).

In the carotenoid biosynthesis pathway, c99806.graphc0 (beta-carotene 3H-hydroxylase), the precursor of ABA synthesis, was significantly down-regulated, while the synthesis and metabolism of ABA, also included in this pathway, and related genes were significantly up-regulated in chilled blueberry. This finding was in agreement with the fact that stress can increase ABA biosynthesis and accumulation as part of the plant defense response. 9-Cis-epoxycarotenoid dioxygenase [NCED; c126206.graphc0) is the most critical rate-limiting enzyme in ABA synthesis and it was significantly up-regulated after storage at 0 °C. Xanthoxin dehydrogenase/ABA2 (c125901.graphc0) is another rate-limiting enzyme in the ABA synthesis pathway, and genes c120652.graphc1 and c122315.graphc1 co-regulate ABA 8′-hydroxylase. The results of our qRT-PCR showed that the expression of beta-carotene 3-hydroxylase was up-regulated in chilled blueberry. The ABA2 genes were significantly up-regulated in chilled blueberry. Expression of Cytochrome P450 707A (CYP707A), which encodes a cytochrome P450 that participates in plant metabolism of terpenoids and polyketides in the ABA metabolic pathway, was significantly up-regulated in chilled blueberry and was 25 times that measured in control blueberry. Thus, cold storage stimulated the expression of genes that encode ABA synthase-related enzymes and are involved in ABA synthesis, but also promoted its ABA metabolism to dihydroxy phaseic acid. Pitting in blueberry after a 30 days’ storage at 0 °C might be affected by these genes; however, the specific mechanism needs further study (Fig. 8c).

Plant hormone signal transduction analysis

The results of transcriptome analysis indicated a large number of DEGs were enriched in the plant hormone signal transduction pathways, especially in IAA signal transduction, ABA signal transduction, SA signal transduction and GA3 signal transduction (Fig. 9).

Fig. 9
figure 9

Transcripts involved in hormone related pathways in 30d-chilled blueberry

Our results indicated that stress-response pathways interact with auxin regulation through the transcription of the Aux/IAA, IAA5, IAA6, and IAA19 are required for stress tolerance; one gene was down-regulated and two genes were up-regulated in the auxin signaling pathway in chilled blueberries. Gene c113524.graphc0 encodes transport inhibitor response1 was down-regulated in chilled blueberry. Genes c128076.graphc0 (IAA), c122133.graphc1 (IAA), and c97770graphc0 coordinate the regulation of IAA, an auxin-responsive protein that participates in plant hormone signal transduction; c111137.graphc1, a member of the GH3 gene family, showed higher expression in chilled blueberry. Genes c102439.graphc0 and c113051.graphc0 belong to the small auxin-upregulated RNA gene family, and are regulated by auxin and environmental factors; these genes were similarly expressed with IAA, which was significantly up-regulated in chilled blueberry. Pyrabactin resistance 1-like (PYL) is an ABA receptor, PYL1, PYL4, PYL5, PYL6, PYL7, PYL10, PYL11, PYL13, and PYL15 have been isolated and identified from Arabidopsis thaliana. In the present study, the expression of c121430.graphc0, an ABA receptor of the PYR/PYL family involved in mitogen-activated protein kinase signaling, was significantly higher in chilled blueberry fruits. The upregulation of both genes indicated that ABA biosynthesis and catabolism were activated by low temperature. The ABA response element binding factor gene (c112990.graphc0) was down-regulated in our transcriptome results, while the brassinosteroid- signaling kinase (c122496.graphc0) was up-regulated in chilled blueberry fruits which was eight times higher than that in control blueberry.

In the SA pathway, gene expression of the plant-pathogen interaction genes c129050.graphc0 (Pathogenesis related protein1) and c126145.graphc0 (TGA) were down-regulated in chilled blueberries. The expression pattern of both genes was similar. Moreover, the treatment of horticultural crops and fruits with GA3 is known to alleviate CI symptoms after treatments at 0 °C [29]. In the GA3 pathway, the expression of DELLA, a negative regulator of GA3 signaling, was up-regulated after 30 days’ storage at 0 °C. Three genes were up-regulated in the CK pathway, a CK receptor and the 2-component response regulators ARR-B or ARR-A. The expression levels of c107755.graphc0 (ARR-B) and c107605.graphc0 (ARR-A) in chilled blueberries were significantly up-regulated, and 18 and 59 times that in control blueberries, respectively. Therefore, genes in hormone signal transduction pathways were significantly affected by cold storage at 0 °C, especially those involved in CK and Aux regulation and metabolism.

TFs in response to cold stress

The different gene expression patterns across the 30d-chilled blueberry indicated that multiple structural genes have contributed to blueberry fruit pitting. In the present study, we screened our assembled transcripts and predicted a total of 1023 TFs from 45 families and identified 738 protein kinase, and 327 transcriptional regulators (TRs); the expressions of most of them in chilled blueberry fruits were changed. The 1023 TFs comprised 42 categories of TFs including 92 C2H2, 87 MYB 68, 74 Ap2/erf-erf, 56 bHLH, 53 C2C2, 51 bZIP, 51 C3H, 45 FAR1, 43WRKY, 39 NAC (Fig. 10).

Fig. 10
figure 10

The numbers and classification of differentially expressed TFs in chilled blueberries. Categories of transcription factors less than 1% of the total are not marked on the pie chart

Validation of the RNA-Seq results by qRT-PCR

To ensure the reliability of the RNA-Seq data, the expression patterns of 40 random DEGs were evaluated by qRT-PCR. (Figs. 11 and 12). The genes represented various functional categories or pathways, including liquid related, defense systems, flavonoid metabolism, brassinosteroid biosynthesis, carotenoid biosynthesis, zeatin biosynthesis and hormone signal transduction pathways. The linear regression showed that the results from RNA-Seq and qRT-PCR were highly relevant (Pearson’s r = 0.8624), despite the difference in the absolute FC between the two methods. The results showed the expression patterns derived by both methods were consistent, confirming the reliability of the transcriptome results. These consistent expression patterns further confirmed the transcriptome data reported in this study constitutes a valuable supplement to the available blueberry genomic and transcriptome information.

Fig. 11
figure 11

The correlation analysis between qRT-PCR of RNA-Seq results of 40 randomly genes

Fig. 12
figure 12

qRT-PCR verification of expression pattern obtained by RNA sequencing. The relative transcription level was calculated according to the 2−ΔΔCT method with actin reference genes as control. Each value is the mean ± SD; P < 0.05 (n = 3, three biological replicates, each with three technical replicates)


Although LT storage is the major method to prevent most fruits and vegetables from decaying, it may cause CI symptoms, such as the surface pitting and dark watery patches in harvested cucumber [16], peel pitting in banana [8], and less aroma in pear [13]. In our present study, the blueberry were stored at low temperature (0 ± 0.5 °C) for 15, 30, 45, and 60 days. This method can preserve the fruit quality and extend their postharvest life. However, pedicle pitting developed due to LT storage and the cold stress; the longer the blueberries are stored, the more serious the pitting is. And the 0 °C-30 days storage was the critical point of pitting. With the occurance of pitting, the phenotype of fruit changed accompanied with the wrinkled cell membrane and the destroyed structure [30]. Meantime, the physiology indices of EL, MDA, proline content, and GSH content, which were considered to reflect physiological state of plant exposure to cold stress [31], were also determined in our study. And the significant increase in EL, content of MDA and proline, and the decrease in GSH content in pitted blueberry indicated that membrane lipid metabolism plays an important role in cold stress response. However, this is only the changes in physiological.

Although cold stress responses and physiological changes have been extensively studied [6, 7, 9, 10], the its molecular and genetic mechanisms remain poorly characterized in blueberry. In order to explore the direct cause of fruit pitting and the gene regulation net under cold stress at the molecular level, we performed deep transcriptome sequencing of the chilled blueberry and control blueberry samples. According to the RNA-Seq data, the three independent biological replicates were highly consistent (0.9462 < R2 < 0.9998), which indicated that the RNA-Seq data was reliable. The transcriptome sequences were assembled into 35,060 transcripts, of which 73.7% were annotated. The annotations provide a valuable resource for investigating specific processes, functions, and pathways during blueberry research. For further verification, we used qRT-PCR, and the expression patterns of 40 significantly changed DEGs corresponded (Pearson’s r = 0.8624) in both methods. Under cold stress, numerous DEGs participating in diverse cell, molecular, and biological pathways were identified as putatively involved in pitting. We mainly focused on the genes involved in lipid metabolism, energy metabolism, antioxidant defense system, hormone signaling networks and TFs based on enrichment results (Tables. 456 and Figs. 89 and S1). These results demonstrate that our transcriptome database is a rich resource for mining cold stress-responsive genes. These transcriptome analysis results provide key points to explain the mechanism underlying the metabolic changes in pitted blueberry in response to cold stress at the molecular level.

Table 4 Genes involved in lipid related pathways in in response to low temperature
Table 5 Genes involved in proline, glutathione, and flavonoid in response to low temperature
Table 6 Genes involved in plant hormone signal transduction in response to low temperature

Lipid metabolism and energy metabolism under cold stress

When the environmental temperature decreases, physical phase and cell arrangement change first; then membrane permeability, lipids’ composition and content begin to degrade. Electrolyte leakage is an important indicator of cell membrane permeability and is the most frequently used method to evaluate plant tissue injury under severe stress conditions. The substantial increase in EL is commonly caused by changes in membrane structure and lipid composition and considered as an index reflecting cell membrane function and integrity [32, 33]. In our present study, the severity of pitting in blueberries was related to relatively higher EL. In addition to electrical leakage, the lipid peroxidation are also evaluated in studies of fruit mechanisms under cold stress [34]. Since MDA content is a direct measure of lipid peroxidation, it can be evaluated to assess the degree of cell damage [35]. In the present study, we observed a significant increase in MDA content during cold storage, especially during the shelf life after 0 °C storage. These clearly indicated that blueberry membrane lipid peroxidation occurs during LT storage and it is intensified during shelf life after 0 °C storage.

Meantime, the accumulation of proline is also related to the cold tolerance in most plant [36]. And the proline has been confirmed to accumulate when plants are exposed to salt or low temperature stress [37]. In our present study, the proline content of blueberries after 30 days’ cold storage was significantly higher than that of other blueberries during the 8-days shelf life; additionally, the proline content of blueberry fruits after 30 days’ cold storage was more than 2 times higher than the blueberries in the other group. Obviously, the increase of proline provided strong evidence that LT storage caused cold stress and blueberries responded to cold stress. In the meantime, we also observed that the GSH in the cold storage group had a downward trend while that in the other group had an upward trend. In conclusion, membrane lipid metabolism plays an important role in blueberry response to cold stress and may be crucially related to the tolerance of fruit to cold stress. This was consistent with some studies in other fruits. Li et al. (2012) [38] confirmed relative EL had a remarkable increase in the mango fruit stored at 5 °C. Kong et al. (2017) [4] also indicated that LT (4 °C) caused serious membrane damage in peppers and the MGDG, PC, PE and PA changed in response to cold. Wang et al. (2019) [30] also identified the blueberries had a higher level of DGDG after LT storage. Additionally, changes of fatty acids were also involved in membrane lipid metabolism. And the change of the lipid not only occurs to the fruits and vegetables in the post-harvest low-temperature storage process, but also occur to the fruits and vegetables or the plants during the cold acclimation [39].

Additionally, genetic and molecular evidence shows that pitting is a complex phenomenon involving the alteration of metabolism with synthesis of specific metabolites, lipids, energy and other pathways. And the changes in gene expression underlie some of the biochemical and physiological changes that occur during cold storage. In the present study, the genes regulating membrane lipid components and fatty acids also significantly changed under cold stress. Six highly differentially expressed pathways of membrane lipid metabolism, including glycerolipid metabolism (Ko00561), glycerophospholipid metabolism (Ko00564), ether lipid metabolism (Ko00565), α-linolenic acid metabolism (Ko00592), sphingolipid metabolism (Ko00600), and phosphatidylinositol signaling system (Ko04070) were significantly changed in chilled blueberries. Additionally, six highly differentially expressed pathways related to fatty acid metabolism including fatty acid biosynthesis (Ko00061), fatty acid elongation (Ko00062), fatty acid degradation (Ko00071), linoleic acid metabolism (Ko00591), biosynthesis of unsaturated fatty acids (Ko01040), and fatty acid metabolism (Ko01212) were also significantly changed in chilled blueberries. Among these genes, the fadD, c109637.graph_c0 encoding palmitoyl-acyl carrier protein was the most obvious up-regulated gene, its expression was 2.2 times higher than that in the control group; the FATB, c126515.graph_c0 encoding acyl-CoA synthetase 1 was the most obvious down-regulated gene, its expression was 2.5 times lower than that the control group. These results suggested the pathways related to membrane lipid had a strong response to cold stress, which was consistent withe the results in loquat [21] and these two genes may be the most ideal genes to further study the regulation of membrane lipid metabolism response to fruit pitting. The study of Die and Rowland. 2014 [40] also confirmed the pathways associated with lipid metabolism, especially regarding enzymes involved in fatty acid elongation, glycerophospholipid metabolism and -linolenic acid metabolism changed significantly in response to cold stress. Moreover, Die and Rowland. 2014 [40] identified the genes encoding phospholipaseC, histidine kinase, 5 lipid-transfer protein (LTP) genes and 2 fatty acid desaturases with potential roles in altering the composition of proteins and lipids.

In fact, the occurrence of CI is not just the results in detrimental changes in membrane structure but also related to the changes of energy status [2, 41]. Published data on other fruits also suggest a close association between membrane integrity regulation is cellular energy status. It has been suggested that the damage of cell membranes was associated with lack of energy status, and ATP played important roles in synthesis of fatty acid and repair of membranes [42]. Recently, many evidences have shown that development of CI in postharvest fruit is partly attributed to limited availability of energy or low energy production, whereas acquaintance of higher levels of ATP and energy charge alleviate CI [16, 43]. Moreover, energy metabolism is important to sustain plant life and plant resistance to environmental stresses [14, 44, 45]. In the present study, there were numerous down-regulated DEGs involved in energy metabolism, which may suggest that the available energy status and the stable enzymatic system in blueberry collectively contribute to improve chilling tolerance by alleviating pitting and maintaining the quality of blueberry fruits during storage at 0 °C [1]. Energy metabolism is linked to adenosine triphosphate (ATP) production and it is the major determinant of cell function and viability, as indicated by plants with higher energy and better maintenance of biological membrane structure being less susceptible to CI [46, 47]. Conversely, long-term storage at 0 °C could result in the decrease of ATP levels in postharvest fruit, causing a disruption in their energy metabolism [38]. In particular, it has been shown that ATP and energy content in papaya decreased after fruit storage at 0 °C, indicating that CI has a close relationship with energy metabolism [45]. Consistent with these findings, our results suggest that energy metabolism might be a main reason underlying the pitting symptoms and therefore exogenous ATP treatment might reduce pitting of blueberry fruits. And it has been proved that the high ATP level can be maintained by increasing the D-pyrroline-5-carboxylatesynthetase (P5CS) activity and then reduce the CI [44]. And suppressing reactive oxygen species [16] modulating proline accumulation [48] may can also enhance the chilling tolerance.

Hormone signals under cold stress

Previous studies have demonstrated that several plant hormones are involved in modulating response and adaptation to a changing environment in plants [49]. As secondary signals, they can initiate a series of signal events (cascade reactions) that ultimately induce stress responsive genes. They may highly correlated with stress responses under LT after harvest and the changes in plant hormone contents can also affect fruit quality. BRs can confer resistance of plants to various abiotic and biotic stresses [50]. Mango fruits treated with Br had higher cold stress tolerance through the regulation of lipids [44]. And exogenous BL can markedly decreased CI incidence of mango fruit. In addition, other plant hormones, including ABA, GA3, JA and ethylene [51, 52], have been also implicated to modulate abiotic and biotic stresses. MeJA treatment can effectively inhibit the CI in postharvest loquat fruit [53, 54] and exhibit higher tolerance to cold stress in postharvest peach fruit by inducing enzyme activities related to energy metabolism and maintaining high levels of ATP and energy charge [55]. Recent studies revealed that exogenous application of SA was able to significantly remit chilling symptoms in harvested loquat fruit [56] and mango fruit. And treatment with ETH or BR, at an effective concentration, could also alleviate CI of tomato fruit [57, 58]. And the role of ABA in the stress-resistant process of the plant is also essential. The content of ABA increases during plant defense responses and this phytohormone plays multiple roles in plant stress responses to drought and cold [59]. ABA is also involved in regulation of cold acclimation. The ABA-insensitive 5 (ABI5) protein was one of the 13 functionally and structurally-related bZIP proteins found to play crucial roles in mediating changes in stress responses [60]. Thus, plant hormones have an important role in low temperature tolerance in postharvest fruit.

In the present study, we found numerous plant-hormone-related genes were responsive to CI. The numbers and expression patterns of these genes are shown in (Fig. 9); in particular, genes involved in BR, CK, Aux, and SA pathways were up-regulated while those involved in ET and methyl jasmonates pathways were not affected by CI. But it may be not consistent with Ding et al., 2015 [51], they revealed that, during cold storage, endogenous levels of GA3 and the expression level of the key GA metabolic genes were lower in chilled fruit compared to that of fruit stored at room temperature. Perhaps because of the different species, the trend of hormone content change under low temperature stress is also different. Additionally, genes involved in Aux and SA regulation were similarly affected CI, probably because the expression of receptor genes decreased, but the expression of PR1 and GH3 increased. Overall, ABA, GA3, BRs, and CK played important roles in pitting caused by cold storage by regulating the expression of numerous downstream genes in blueberries. These results indicate that hormone signaling pathways may play unique functions under cold stress conditions, and may be forming a complex antioxidant defense system in blueberry. However, whether and how ABA, GA3, BRs, or CK mediated signaling participates in cold stress responses in blueberry need to be further determined.

TFs involved in blueberry pitting responses to cold stress

In many biological processes, defense responses also require the regulation of specific TFs, which comprise one of the complex regulatory networks in plants [14]. TFs may play essential roles in stress responses by regulating their target genes through specific binding to cis-acting elements in their promoters. Although the TFs families play diverse roles in plant developmental processes and environmental responses, most of them have been reported to be linked to cold stress resistance in plants. For example, six BrbZIP genes have been identified as putative key factors in cold stress response [23]. NAC genes are also reported to be induced by cold in Arabidopsis [61], chrysanthemum [62], rice [63] and sugarcane [64]. Additionally, the NAC TFs can also interact with CBF1 to regulate cold tolerance in apples [65] and bananas [66]. The bHLH family, the second largest family [67], can also increase cold tolerance in plants by regulating the expression of ROS clearance-related and stress-responsive genes. Meanwhile, the first WRKY transcription factor was found in sweet potato under cold stress [68]; in cucumber, CsWRKY46 confers cold tolerance and positively regulates ABA-dependent cold signaling pathways [69].

Some less common TFs like C2H2-zinc finger proteins are also essential in plant stress responses, although their transcriptional regulatory mechanisms remain largely unclear [70]. Overall, our results revealed that the 1023 TFs comprised 42 categories of TFs including 92 C2H2, 87 MYB 68, 74 Ap2/erf-erf, 56 bHLH, 53 C2C2, 51 bZIP, 51 C3H, 45 FAR1, 43WRKY, 39 NAC were changed in the pitting process; this was a little similarity with Die and Rowland (2014) [40]; they found the seven most highly represented TF families in response to cold stress were WRKY, ARF, C3H, AP2/ERF, bHLH, C2H2, and NAC. Based on our current results, the bHLH family was significantly up-regulated in chilled blueberries, but Ap2, ZIP, and WRKY families were significantly down-regulated under this condition. In addition, Plant hormones and transcription factors may also interact during cold stress. For instance, NAC TFs have been demonstrated to regulate the expression of ABA-related genes during abiotic stress responses [71]. And ethylene was shown to negatively regulate cold signaling at least partially through the direct transcriptional control of CBFs [52]. Since our current research focus on further confirming whether the C2H2, MYB, and NAC TFs play an important role in the response to LT, and how they regulate the expression of downstream genes in pitted blueberry.


Overall, the harvest longterm cold storage leads serious CI in blueberries, which substantially decreases the quality, storability and consumer acceptance. Our results identified the occurrence of pitting, and proved 30 days’ cold storage is the critical point of pitting. We found higher MDA content, proline content, EL and lower GSH content in 30d-chilled berries. Additionally, a comprehensive transcriptome profile of blueberry under cold stress were explored by RNA-seq this time in order to characterize the CI mechanism at the molecular level. We identified and summarized the genes in response to cold stress as follows: (1) transcription factors; (2) membrane lipid and energy metabolism; (3) defensive; and (4) hormone biosynthesis and signal transduction, including the genes involved in membrane lipids, proline, glutathione, flavonoids, brassinosteroid, carotenoid, and zeatin biosynthesis pathway. These DEGs and the biological processes related to the stress responses, lipid and hormone metabolic processes almost all accumulated and up-regulated in the CI process. And may play a crucial role in the response to cold stress in blueberry. These results provides novel insights into a series of molecular mechanisms underlying physiological metabolism and defense. Further, the enormous amount of blueberry transcriptome data generated here will serve as the foundation in finding available effective way in alleviating CI and also help to perfect the lack of the genetic information of non-model plant species. In future work, we aim to focus on the function verification of special structural genes fadD (c109637.graph_c0) related to membrane lipid and fatty acid metabolism and 8′-hydroxylas (c122315.graph_c0) related to plant hormone and explore how they work in chilling process.


Plant materials and temperature treatments

‘Duke’ blueberry (Vaccinium corymbosum L.) fruits were harvested at commercial maturity from an blueberry base in Shenyang (N41°39′55.82″, E123°05′12.46″), Liaoning Province, China, from April to July in 2016 to 2018. The blueberries were manually picked to prevent any mechanical damage. After 2 h pre-cooling at 20 °C, blueberry fruits without physical injury, disease, or rot were enclosed in plastic boxes; then stored at different storage temperatures (20 °C and 0 °C). Each box contained 125 g blueberry fruits and the relative humidity was kept at 85%. The fleshy tissues were collected, quickly frozen in liquid nitrogen, and stored at − 80 °C [72]. The blueberry samples non-pitted (T01, T02, and T03) were defined as control blueberry and the pitted blueberry stored at 0 °C for 30 days (T04, T05, and T06) were defined as 30d-chilled blueberry. They were collected for RNA-Seq and the RNA-Seq was did by Biomarker Technologies (Beijing, China). During this period, we ground 20 different blueberry fruits of each samples (T01–06), and in each samples (T01–06), we also perform three technical replicates.

Evaluation of pitting rate, decay rate, electrolyte leakage, MDA and proline

Pitting and decay rate were calculated for fruits at 0, 2, 4, 6, and 8 days shelf-life, with or without previous storage at 0 °C, using the following equations:

$$ \mathrm{Pitting}\ \mathrm{rate}\ \left(\%\right)=\mathrm{An}/\mathrm{Am}\times 100 $$
$$ \mathrm{Decay}\ \mathrm{rate}\ \left(\%\right)=\mathrm{Ai}/\mathrm{Am}\times 100 $$

In these equations, An is the number of blueberry fruits with pitting, Ai is the number of blueberry fruits cankered, juicy, or mildew, and Am is the total number of blueberry fruits. Three independent measurements were performed at each time point for 150 fruits per replicate sample.

The cell membrane permeability determined by electrical conductivity meter according to the method described in Mao et al. (2007) [73]. MDA content was determined according to the method described in Zhao et al. (2005) [74]. The proline (Pro) content and reduced glutathione (GSH) content were determined according to the method described in Zhao et al. (2009) [75].

Total RNA isolation, cDNA library construction, and sequencing

Total RNAs from all samples were isolated using the OminiPlant RNA kit (CWBio, Beijing, China). The integrity was monitored on 1% agarose gels and the purity was checked using the NanoDrop spectrophotometer (Thermo Fisher Scientific, Inc.) according to Niu et al. (2017) [76]. Sequencing libraries were generated using NEBNext®Ultra™ RNA Library Prep Kit (Illumina®NEB, USA) and index codes were added to attribute sequences to each sample. The library quality was assessed on the Agilent Bioanalyzer 2100 system. Finally, the cDNA libraries were sequenced on the Illumina® HiSeq 2000 platform (Illumina Inc.) and paired-end reads were generated at Beijing Biomarker Technology Company (Beijing, China) [77]. All the downstream analyses were based on clean data with high quality. And all the sequencing data generated in this study have been deposited in the NCBI with the link of, under the accession number SAMN13049031, the Temporary Submission ID SUB643610 and the BioProject ID PRJNA578101.

Functional annotation and differential expression analysis of genes

Because of the absence of a accurate reference-grade genome for blueberry, a transcriptome was de novo assembled from this RNAseq data and other publicly available transcriptome data from blueberry downloaded from the National Center for Biotechnology Information (NCBI). And the non-redundant unigenes were annotated based on NCBI non-redundant (NR); Protein family (Pfam); Swiss-Prot, Clusters of Orthologous Groups (KOG/COG/eggNOG); Kyoto Encyclopedia of Genes and Genomes (KEGG) and Gene Ontology (GO) databases according to Liu et al. (2015) [77].

DESeq R package (1.10.1) was used to identify the DEGs between control and chilled groups with the standard P-value< 0.05 and |log2 (fold change)| ≥ 1. The GO enrichment and KEGG analysis of the DEGs was performed to identify significantly enriched biological pathways of DEGs, with an FDR adjusted P-value< 0.05.

Validation of the RNA-Seq data and expression analysis by qRT-PCR

Quantitative real-time PCR was performed to validate the expression patterns obtained from the RNA-Seq data. Gene specific primers were designed in Primer 5.0 and synthesized by Genewiz Biotechnology Synthesis Lab (Jiangsu, China) (Table S1 Primers used in qRT-PCR). The program and reaction conditions for the 2 × ULtraSYBR Mixture (Low ROX) Kit (ComWin Biotech, Beijing, China) using QuantStudio 6 Flex (Life Technologies) were as described in the manufacturer protocols. Template cDNA was diluted five-fold, and 2 μl of the diluted cDNA was then added to each 20-μl PCR mixture. The qRT-PCR was performed 2-step method. Each sample was carried out with three independent biological replicates and three technical replicates. Relative expression levels were calculated by the comparative 2−ΔΔCT method using the expression of actins as the internal control.

Statistical analyses

One-way analysis of variance (ANOVA) with the least significant difference (LSD) was performed for all data using SPSS 20.0 software (BM Corp, Armonk, NY, USA). Each experiment was repeated three times independently and the represent significance was considered to P < 0.05. Origin 8.1 ( and R3.5.1 ( were used to display the data obtained from the experiments.

Availability of data and materials

The data generated and analysed in this study are available from the corresponding author on reasonable request. And all the sequencing data generated in this study have been deposited in the NCBI with the link of, under the accession number SAMN13049031, the Temporary Submission ID SUB643610 and the BioProject ID PRJNA578101.



Abscisic acid


Adenosine triphosphate






Basic leucine-zipper


C-Repeat binding factor


Chilling injury




Differentially expressed gene


Electrolyte leakage




False discovery rate


Fragments per kilobase million




Gene ontology




Glutathione s-transferase


Hydroxycinnamoyl-CoA shikimate


Indole acetic acid


Jasmonic acid


Leuco anthocyanidin reductase




Low temperature




Non-redundant protein sequences


Phenylalanine ammonia lyase


Protein family


Pyrabactin resistance 1-like


Pyrabactin resistance


Reactive oxygen species


Room temperature


Salicylic acid


Transcription factors


Transcriptional regulator


  1. Zhou Q, Zhang CL, Cheng SC, Wei BD, Liu XY, Ji SJ. Changes in energy metabolism accompanying pitting in blueberries stored at low temperature. Food Chem. 2014;164:493–501.

    Article  CAS  PubMed  Google Scholar 

  2. Zhou Q, Ma C, Cheng SH, Wei BD, Liu XY, Ji SJ. Changes in antioxidative metabolism accompanying pitting development in stored blueberry fruit. Postharvest Biol Technol. 2014:88–95.

    Article  CAS  Google Scholar 

  3. Ge WY, Kong XM, Zhao YB, Wei BD, Zhou Q, Ji SJ. Insights into the metabolism of membrane lipid fatty acids associated with chilling injury in post-harvest bell peppers. Food Chem. 2019;295:26–35.

    Article  CAS  PubMed  Google Scholar 

  4. Kong XM, Wei BD, Gao Z, Zhou Y, Shi F, Zhou X, Ji SJ. Changes in membrane lipid composition and function accompanying chilling injury in bell peppers. Plant Cell Physiol. 2017:167–78.

    Article  CAS  Google Scholar 

  5. Purvis AC. Diphenylamine reduces chilling injury of green bell pepper fruit. Postharvest Biol Technol. 2002;25:41–8.

    Article  CAS  Google Scholar 

  6. Lim CS, Kang SM, Cho JL. Bell pepper (Capsicum annuum L) fruits are susceptible to chilling injury at the breaker stage of ripeness. HortScience. 2007;42:1659–64.

    Article  CAS  Google Scholar 

  7. Wang Q, Ding T, Zuo JH, Gao LP, Fan LL. Amelioration of postharvest chilling injury in sweet pepper by glycine betaine. Postharvest Biol Technol. 2016;112:114–20.

    Article  CAS  Google Scholar 

  8. Chen JY, He LH, Jiang YM, Wang Y, Joyce DC, Lu WJ. Role of phenylalanine ammonia-lyase in heat pretreatment-induced chilling tolerance in banana fruit. Physiol Plant. 2008;132:318–28.

    Article  CAS  PubMed  Google Scholar 

  9. Lurie S, Crisosto CH. Chilling injury in peach and nectarine. Postharvest Biol Technol. 2005;37(3):195–208.

    Article  Google Scholar 

  10. Crisosto CH, Mitchell FG, Ju Z. Susceptibility to chilling injury of peach nectarine and plum cultivars grown in California. HortScience. 1999;34:1116–8.

    Article  Google Scholar 

  11. Rémy R, Mathieu L, Claudie DM. Antioxidant and enzymatic responses to oxidative stress induced by cold temperature storage and ripening in mango (Mangifera indica L , cv , ‘Cogshall’) in relation to carotenoid content. J Plant Physiol. 2018;224–225:75–85.

    Google Scholar 

  12. Liu H, Song L, You Y, Li YB, Duan XW, Jiang YM, Lu WJ. Cold storage duration affects litchi fruit quality membrane permeability enzyme activities and energy charge during shelf time at ambient temperature. Postharvest Biol Technol. 2011;60:24–30.

    Article  CAS  Google Scholar 

  13. Shi F, Zhou X, Zhou Q, Tan Z, Yao MM, Ji SJ. Membrane lipid metabolism changes and aroma ester loss in low temperature stored ‘Nanguo’ pears. Food Chem. 2018;245:446–53.

    Article  CAS  PubMed  Google Scholar 

  14. Shi F, Zhou X, Zhou Q, Tan Z, Yao MM, Ji SJ. Transcriptome analyses provide new possible mechanisms of aroma ester weakening of ‘Nanguo’ pear after cold storage. Sci Hortic. 2018;237:247–56.

    Article  CAS  Google Scholar 

  15. Zhou X, Dong L, Li R, Zhou Q, Wang JW, Ji SJ. Low temperature conditioning prevents loss of aroma-related esters from ‘Nanguo’ pears during ripening at room temperature. Postharvest Biol Technol. 2014;100:23–32.

    Article  CAS  Google Scholar 

  16. Chen B, Yang H. 6-Benzylaminopurine alleviates chilling injury of postharvest cucumber fruit through modulating antioxidant system and energy status. J Sci Food Agric. 2013;93:1915–21.

    Article  CAS  PubMed  Google Scholar 

  17. Kratsch HA, Wise RR. The ultrastructure of chilling stress. Plant Cell Environ. 2000;23:337–50.

    Article  CAS  Google Scholar 

  18. Singh Z, Singh RK, Sane VA, Nath P. Mango postharvest biology and biotechnology. Crit Rev Plant Sci. 2013;32:217–36.

    Article  CAS  Google Scholar 

  19. Kong XM, Zhou Q, Luo F, Wei BD, Wang YJ, Zhou X, Ji SJ. Transcriptome analysis of harvested bell peppers (Capsicum annuum L,) in response to cold stress. Plant Physiol Biochem. 2019;139:314–24.

    Article  CAS  PubMed  Google Scholar 

  20. Zhao LY, Gong X, Gao JZ, Dong HZ, Zhang SL, Huang XS. Transcriptomic and evolutionary analyses of white pear (Pyrus bretschneideri) β-amylase genes reveals their importance for cold and drought stress responses. Gene. 2019;689:102–13.

    Article  CAS  PubMed  Google Scholar 

  21. Lou XM, Wang KH, Ni XP, Gao ZH, Iqbal S. Integrating proteomic and transcriptomic analyses of loquat (Eriobotrya japonica Lindl ,) in response to cold stress. Gene. 2018; 677: 57–65.

    Article  CAS  PubMed  Google Scholar 

  22. Baloglu MC, Eldem V, Hajyzadeh M, Unver T. Genome-wide analysis of the bZIP transcription factors in cucumber. PLoS One. 2014;9(4):e96014.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  23. Dubos C, Stracke R, Grotewold E, Weisshaar B, Martin C, Lepiniec L. MYB transcription factors in Arabidopsis. Trends Plant Sci. 2010;15:573–81.

    Article  CAS  PubMed  Google Scholar 

  24. Fowler S, Thomashow M. Arabidopsis transcriptome profiling indicates that multiple regulatory pathways are activated during cold acclimation in addition to the CBF cold response pathway. Plant Cell. 2005;14:1675–90.

    Article  CAS  Google Scholar 

  25. Jolma A, Nitta KR, Dave K, Popov A, Taipale M, Taipale J. DNA-dependent formation of transcription factor pairs alters their binding specificity. Nature. 2015;527:384–388.

    Article  CAS  PubMed  Google Scholar 

  26. Kauzo S, Kauzo YS, Seki M. Regulatory network of gene expression in the drought and cold stress responses. Curr Opin Plant Biol. 2003;6:410–7.

    Article  CAS  Google Scholar 

  27. Liu L, Cao XL, Bai R, Yao N, Li LB, He CF. Isolation and characterization of the cold-induced Phyllostachys edulis AP2/ERF family transcription factor peDREB1. Plant Mol Biol. 2012;30:679–89.

    Article  CAS  Google Scholar 

  28. Thomashow MF, Gilmour SJ, Stockinger EJ, Zarka DG. Role of the Arabidopsis CBF transcriptional activators in cold acclimation. Physiol Plant. 2001;112:171–5.

    Article  CAS  Google Scholar 

  29. Tian WN, Lv YC, Cao JK, Jiang WB. Retention of iceberg lettuce quality by low temperature storage and postharvest application of 1-methylcyclopropene or gibberellic acid. J Food Sci Technol. 2014;51(5):943–9.

    Article  CAS  PubMed  Google Scholar 

  30. Wang YJ, Ji SJ, Dai HY, Kong XM, Hao J, Wang SY, Zhou Q. Changes in membrane lipid metabolism accompany pitting in blueberry during refrigeration and subsequent storage at room temperature. Front Plant Sci. 2019;10:1–12.

    Article  Google Scholar 

  31. Zhao DY, Shen L, Fan B, Liu KL, Yu MM, Zheng Y, Sheng JP. Physiological and genetic properties of tomato fruits from 2 cultivars differing in chilling tolerance at cold storage. J Food Sci. 2009;74:348–52.

    Article  CAS  Google Scholar 

  32. Kwon SY, Jeong YJ, Lee HS, Kim JS, Cho KY, Allen RD, Kwak SS. Enhanced tolerance of transgenic tobacco plants expressing both superoxide dismutase and ascorbate peroxidase in chloroplasts against methyl viologen mediated oxidative stress. Plant Cell Environ. 2002;25:873–82.

    Article  Google Scholar 

  33. Yabuta Y, Motoki T, Yoshimura K, Takeda T, Ishikawa T, Shigeoka S. Thylakoid membrane-bound ascorbate peroxidase is a limiting factor of antioxidative systems under photo-oxidative stress. Plant J. 2002;32:915–25.

    Article  CAS  PubMed  Google Scholar 

  34. Zhang XH, Shen L, Li FJ, Zhang YX, Meng DM, Sheng JP. Up-regulating arginase contributes to amelioration of chilling stress and the antioxidant system in cherry tomato fruits. J Sci Food Agric. 2010;90:2195–202.

    Article  CAS  PubMed  Google Scholar 

  35. Xu XZ, Guo RR, Cheng CX, Zhang HJ, Zhang YC, Wang XP. Overexpression of ALDH2B8 an aldehyde dehydrogenase gene from grapevine sustains Arabidopsis growth upon salt stress and protects plants against oxidative stress. Plant Cell Tissue Organ Cult. 2013;114:187–96.

    Article  CAS  Google Scholar 

  36. Chen WP, Li PH. Membrane stabilization by abscisic acid under cold aids proline in alleviating chilling injury in maize (Zea mays L,) cultured cells. Plant Cell Environ. 2002;25:955–62.

    Article  CAS  Google Scholar 

  37. Iba K. Acclimative response to temperature stress in higher plants: approaches of gene engineering for temperature tolerance. Annu Rev Plant Biol. 2002;53:225–45.

    Article  CAS  PubMed  Google Scholar 

  38. Li B, Zhang C, Cao B, Qin GZ, Wang WH, Tian SP. Brassinolide enhances cold stress tolerance of fruit by regulating plasma membrane proteins and lipids. Amino Acids. 2012;43:2469–80.

    Article  CAS  PubMed  Google Scholar 

  39. Moellering ER, Muthan B, Benning C. Freezng tolerance in plants requires lipid remodeling at the outer chloroplast membrane. Science. 2010;330:227–30.

    Article  CAS  Google Scholar 

  40. Die JV, Rowland LJ. Elucidating cold acclimation pathway in blueberry by transcriptome profiling. Environ Exp Bot. 2014;06:87–98.

    Article  CAS  Google Scholar 

  41. Lee DH, Lee CB. Chilling stress-induced changes of antioxidant enzymes in the leaves of cucumber: in gel enzyme activity assays. Plant Sci. 2002;159:75–85.

    Article  Google Scholar 

  42. Rawyler A, Pavelic D, Gianinazzi C, Oberson J, Braendle R. Membrane lipid integrity relies on a threshold of ATP production rate in potato cell cultures submitted to anoxia. Plant Physiol. 1999;120:293–300.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  43. Jiang YM, Jiang YL, Qu HX, Duan XW, Luo YB, Jiang WB. Energy aspects in ripening and senescence of harvested horticultural crops. Stewart Postharvest Rev. 2007;4:1–5.

    Google Scholar 

  44. Li P, Zheng XL, Liu Y, Zhu YY. Storage application of oxalic acid alleviates chilling injury in mango fruit by modulating proline metabolism and energy status under chilling stress. Food Chem. 2014;142:72–8.

    Article  CAS  PubMed  Google Scholar 

  45. Jin P, Zhu H, Wang L, Shan T, Zheng YH. Oxalic acid alleviates chilling injury in peach fruit by regulating energy metabolism and fatty acid contents. Food Chem. 2014;161:87–93.

    Article  CAS  PubMed  Google Scholar 

  46. Dong CJ, Li L, Shang QM, Liu XY, Zhang ZG. Endogenous salicylic acid accumulation is required for chilling tolerance in cucumber (Cucumissativus L,) seedlings. Planta. 2014;240:687–700.

    Article  CAS  PubMed  Google Scholar 

  47. Wang Z, Ma L, Zhang XF, Xu LM, Cao JK, Jiang WB. The effect of exogenous salicylic acid on antioxidant activity bioactive compounds and antioxidant system in apricot fruit. Sci Hortic. 2015;181:113–20.

    Article  CAS  Google Scholar 

  48. Cao SF, Cai YT, Yang ZF, Zheng YH. MeJA induces chilling tolerance in loquat fruit by regulating proline and c-aminobutyric acid contents. Food Chem. 2012;133:1466–70.

    Article  CAS  Google Scholar 

  49. Belkhadir Y, Chory J. Brassinosteroid signaling: a paradigm for steroid hormone signaling from the cell surface. Science. 2006;314:1410–1.

    Article  CAS  PubMed  Google Scholar 

  50. Bajguz A, Hayat S. Effects of brassinosteroids on the plant responses to environmental stresses. Plant Physiol Biochem. 2009;47:1–8.

    Article  CAS  PubMed  Google Scholar 

  51. Ding Y, Sheng JP, Li SY, Nie Y, Zhao JH, Tang XM. The role of gibberellins in the mitigation of chilling injury in cherry tomato (Solanum lycopersicum L) fruit. Postharvest Biol Technol. 2015;101:88–95.

    Article  CAS  Google Scholar 

  52. Shi Y, Tian S, Hou L, Huang X, Zhang X, Guo H, Yang S. Ethylene signaling negatively regulates freezing tolerance by repressing expression of CBF and type-a ARR genes in Arabidopsis. Plant Cell. 2012;24:2578–95.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  53. Cao SF, Zheng YH, Yang ZF, Tang SS, Jin P, Wang KT, Wang XM. Effect of methyl jasmonate on the inhibition of Colletotrichum acutatum infection in loquat fruit and the possible mechanisms. Postharvest Biol Technol. 2008;49:301–7.

    Article  CAS  Google Scholar 

  54. Cao SF, Zheng YH, Wang KT, Rui HJ. Effect of methyl jasmonate on cell wall modification of loquat fruit in relation to chilling injury after harvest. Food Chem. 2010;118:641–7.

    Article  CAS  Google Scholar 

  55. Meng XH, Han J, Wang Q, Tian SP. Changes in physiology and quality of peach fruit treated by methyl jasmonate under low temperature stress. Food Chem. 2009;106:501–8.

    Article  CAS  Google Scholar 

  56. Cai C, Li X, Chen KS. Acetylsalicylic acid alleviates chilling injury of postharvest loquat (Eriobotrya japonica Lindl,) fruit. Eur Food Res Technol. 2006;223:533–9.

    Article  CAS  Google Scholar 

  57. Zhao DY, Shen L, Fan B, Yu MM, Zheng Y, Lv SN, Sheng JP. Ethylene and cold participate in the regulation of LeCBF1 gene expression in postharvest tomato fruits. FEBS Lett. 2009;583:3329–34.

    Article  CAS  PubMed  Google Scholar 

  58. Aghdam MS, Mohammadkhani N. Enhancement of chilling stress tolerance of tomato fruit by postharvest brassinolide treatment. Food Bioprocess Technol. 2014;7:909–14.

    Article  CAS  Google Scholar 

  59. Qi XN, Y Y, Fan ZQ, Chen JY, Lu WJ, Kuang JF. A banana fruit transcriptional repressor MaERF10 interacts with MaJAZ3 to strengthen the repression of JA biosynthetic genes involved in MeJA-mediated cold tolerance. Postharvest Biol Technol. 2016; 120: 222–231.

    Article  CAS  Google Scholar 

  60. Fujita Y, Fujita M, Shinozaki K, Kazuko YS. ABA-mediated transcriptional regulation in response to osmotic stress in plants. J Plant Res. 2010;124:509–25.

    Article  CAS  Google Scholar 

  61. Yang X, Wang X, Ji L, Yi ZL, Fu CX, Hu RB, Zhou GH. Over expression of a Miscanthus lutarioriparius NAC gene MlNAC5 confers enhanced drought and cold tolerance in Arabidopsis. Plant Cell Rep. 2015;34:943.

    Article  CAS  PubMed  Google Scholar 

  62. Huang H, Wang Y, Wang S, Wu X, Yang K, Niu Y, Dai S. Transcriptome-wide survey and expression analysis of stress responsive NAC genes in Chrysanthemum lavandulifolium. Plant Sci. 2012;193–194:18–27.

    Article  PubMed  CAS  Google Scholar 

  63. Hu HH, You J, Fang YJ, Zhu XY, Qi ZY, Xiong LZ. Characterization of transcription factor gene SNAC2 conferring cold and salt tolerance in rice. Plant Mol Biol. 2008;67:169–81.

    Article  CAS  PubMed  Google Scholar 

  64. Nogueira FTS, Schlögl PS, Camargo SR, Fernandez JH, RVE D Jr, Pompermayer P, Arruda P. SsNAC23 a member of the NAC domain protein famiy is associated with cold herbivory and water stress in sugarcane. Plant Sci. 2005;169:93–106.

    Article  CAS  Google Scholar 

  65. An JP, Li R, Qu FJ, You CX, Wang FX, Hao YJ. An apple NAC transcription factor negatively regulates cold tolerance via CBF-dependent pathway. J Plant Physiol. 2018;221:74–80.

    Article  CAS  PubMed  Google Scholar 

  66. Shan W, Kuang JF, Lu WJ, Chen JY. Banana fruit NAC transcription factor MaNAC1 is a direct target of MaICE1 and involved in cold stress through interacting with MaCBF1. Plant Cell Environ. 2014;37:2116–27.

    Article  CAS  PubMed  Google Scholar 

  67. Feller A, Machemer K, Braun EL, Grotewold E. Evolutionary and comparative analysis of MYB and bHLH plant transcription factors. Plant J. 2011;66(1):94–116.

    Article  CAS  PubMed  Google Scholar 

  68. Ishiguro S, Nakamura K. Characterization of a cDNA encoding a novel DNA binding protein SPF1 that recognizes SP8 sequences in the 50 upstream regions of genes coding for sporamin and b-amylase from sweet potato. Mol Gen Genet. 1994;244:563–71.

    Article  CAS  PubMed  Google Scholar 

  69. Zhang Y, Yu HJ, Yang XY, Li Q, Ling J, Wang H, Jiang WJ. CsWRKY46 a WRKY transcription factor from cucumber confers cold resistance in transgenic-plant by regulating a set of cold-stress responsive genes in an ABA-dependent manner. Plant Physiol Biochem. 2016;108:478–87.

    Article  CAS  PubMed  Google Scholar 

  70. Han YC, Fu CC, Kuang JF, Chen JY, Lu WJ. Two banana fruit ripening-related C2H2 zinc finger proteins are transcriptional repressors of ethylene biosynthetic genes. Postharvest Biol Technol. 2016;116:8–15.

    Article  CAS  Google Scholar 

  71. Puranik S, Sahu PP, Srivastava PS. Prasad. NAC proteins: regulation and role in stress tolerance. Trends Plant Sci. 2012;17:369–81.

    Article  CAS  PubMed  Google Scholar 

  72. Sheng L, Zhou X, Liu ZY, Wang JW, Zhou Q, Zhang Q, Ji SJ. Changed activities of enzymes crucial to membrane lipid metabolism accompany pericarp browning in ‘Nanguo’ pears during cold storage and subsequent shelf-life at room temperature. Postharvest Biol Technol. 2016;117:1–8.

    Article  CAS  Google Scholar 

  73. Mao LC, Wang GZ, Zhu CG, Pang HQ. Involvement of phospholipase D and lipoxygenase in response to chilling stress in postharvest cucumber fruits. Plant Sci. 2007;172:400–5.

    Article  CAS  Google Scholar 

  74. Zhao YF, Lin HT, Lin JF, Chen SJ, Xi YF. Changes of respiration intensity membrane permeability and quality of ‘lonyan’ fruit after harvest. J Fujian Agric Forestry Univ. 2005;34(2):263–7.

    Google Scholar 

  75. Zhao MG, Chen L, Zhang LL, Zhang WH. Nitric reductase dependent nitric oxide production is involved in cold acclimation and freezing tolerance in Arabidopsis. Plant Physiol. 2009;151:755–67.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  76. Niu XS, Wang YY, Li M, Zhang XR, Wu YT. Transcriptome analysis of avian reovirus-mediated changes in gene expression of normal chicken fibroblast DF-1 cells. BMC Genomics. 2017:18–911.

  77. Liu MM, Xing YM, Zhang DW, Guo SX. Transcriptome analysis of genes involved in defence response in Polyporus umbellatus with Armillaria mellea infection. Sci Rep. 2015;5:16075.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

Download references


We thank Professor Shujuan Ji, Department of Food, Shenyang Agricultural University for her guidance.


Material collection, transcriptional sequencing and data analyses in this study were financially supported by the “National Natural Science Foundation of China (31501534); the gene verification and physiological experiments were financially supported by the 64th China Postdoctoral Science Foundation (2018 M640260).

Author information

Authors and Affiliations



ZQ and JSJ designed this experiment; WBD, CSC, ZF, WYJ, HJ, WSY researched literature; ZF did experiments (RNA isolation, qRT-PCR), data analysis and manuscript editing; ZQ and JSJ revised this manuscript. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Qian Zhou.

Ethics declarations

Ethics approval and consent to participate

Not Applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

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

Supplementary information

Additional file 1: Table S1.

Primers used in qRT-PCR. Figure S1. Graphical Abstract.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (, which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Zhang, F., Ji, S., Wei, B. et al. Transcriptome analysis of postharvest blueberries (Vaccinium corymbosum ‘Duke’) in response to cold stress. BMC Plant Biol 20, 80 (2020).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: