Skip to main content

Advertisement

Isolation and functional characterization of cold-regulated promoters, by digitally identifying peach fruit cold-induced genes from a large EST dataset

Abstract

Background

Cold acclimation is the process by which plants adapt to the low, non freezing temperatures that naturally occur during late autumn or early winter. This process enables the plants to resist the freezing temperatures of winter. Temperatures similar to those associated with cold acclimation are also used by the fruit industry to delay fruit ripening in peaches. However, peaches that are subjected to long periods of cold storage may develop chilling injury symptoms (woolliness and internal breakdown). In order to better understand the relationship between cold acclimation and chilling injury in peaches, we isolated and functionally characterized cold-regulated promoters from cold-inducible genes identified by digitally analyzing a large EST dataset.

Results

Digital expression analyses of EST datasets, revealed 164 cold-induced peach genes, several of which show similarities to genes associated with cold acclimation and cold stress responses. The promoters of three of these cold-inducible genes (Ppbec1, Ppxero2 and Pptha1) were fused to the GUS reporter gene and characterized for cold-inducibility using both transient transformation assays in peach fruits (in fruta) and stable transformation in Arabidopsis thaliana. These assays demonstrate that the promoter Pptha1 is not cold-inducible, whereas the Ppbec1 and Ppxero2 promoter constructs are cold-inducible.

Conclusion

This work demonstrates that during cold storage, peach fruits differentially express genes that are associated with cold acclimation. Functional characterization of these promoters in transient transformation assays in fruta as well as stable transformation in Arabidopsis, demonstrate that the isolated Ppbec1 and Ppxero2 promoters are cold-inducible promoters, whereas the isolated Pptha1 promoter is not cold-inducible. Additionally, the cold-inducible activity of the Ppbec1 and Ppxero2 promoters suggest that there is a conserved heterologous cold-inducible regulation of these promoters in peach and Arabidopsis. These results reveal that digital expression analyses may be used in non-model species to identify candidate genes whose promoters are differentially expressed in response to exogenous stimuli.

Background

Cold temperature is an environmental factor that plays an important role in plant growth and development. Temperate plants have developed mechanisms to adapt to periods of low non-freezing temperatures, enabling these plants to survive subsequent freezing temperatures. This process is called cold acclimation [1]. Cold acclimation is a complex process that involves physiological, biochemical and molecular modifications [24]. Hundreds of genes have been shown to have altered expression levels during cold acclimation [5]. These alterations enable the plant to withstand freezing by creating a chronic response that protects the integrity of the cellular membranes, enhances anti-oxidative mechanisms and accumulates molecular cryoprotectants [6].

Under normal conditions, cold acclimation is initiated by the cold temperatures of late fall and early winter, when fruit trees lack fruits. Similar cold temperatures have been used in the fruit industry to store fruits for prolonged periods of time. These temperatures inhibit fruit ripening, thereby extending fruit postharvest life. Despite the benefits, peaches that are subjected to long periods of cold storage can develop chilling injury symptoms (i.e. woolliness and internal breakdown) which reduce the postharvest quality of these fruits and results in significant economical losses [79].

Most of the efforts directed towards understanding the molecular basis of cold acclimation have been performed in the model plant A. thaliana [14]. Little is known about what occurs under low, non-freezing temperatures in fruits or fruit trees. Since chilling injury occurs in fruits that have undergone long-term cold storage, perhaps cold acclimation processes are associated with this injury. A better understanding of cold acclimation and cold-responsive genes in peach trees and fruits may provide clues about the association of cold acclimation and chilling injury.

Several transcription factors associated with cold acclimation have been shown to regulate the expression of cold-inducible genes containing conserved ABRE (abscisic acid response elements) and/or DRE (dehydration-responsive) elements in their promoters [1013]. The regulation of cold-inducible promoters in peaches may be mediated by the interaction between promoters containing these types of cis-elements and orthologous transcription factors. However, the identification and functional characterization of these types of promoters in fruit trees is lacking.

We have demonstrated previously that there is a conserved heterologous regulation of the wheat putative high-affinity Pi transporter, TaPT2 in both monocots (wheat) and dicots (Arabidopsis) [14]. These findings demonstrate that Arabidopsis may be used as a heterologous system to test the functionality of promoters. However, this type of heterologous regulation may not exist for all promoters and may not be conserved among all plant species. An alternative to functional analyses in heterologous systems is transient transformation of fruits using agro-infiltration. Agro-infiltration of fruits have been performed to test the activity of the 35S CaMV promoter fused to reporter genes such as GUS or luciferase in tomatoes, apples, pears, peaches, strawberries and oranges [15, 16]. However, to our knowledge, it has not been used to determine the activity of cold-inducible promoters within the fruit (in fruta).

To identify cold-responsive genes expressed in peach fruits, digital expression analyses of ESTs from fruits exposed to four different postharvest conditions were analyzed [17]. Isolation of the promoter regions of three genes highly expressed in fruits that have undergone long-term cold storage, allowed us to identify common regulatory elements present in these promoters. Functional characterization of these promoters (stably in A. thaliana and transiently in peach fruits) demonstrates that these are peach cold-inducible promoters and that there is a conserved heterologous regulation of these promoters in peach and Arabidopsis.

Methods

Digital expression analyses

We have previously described the contigs used in this work [17]. The ESTs that make up these contigs represent transcripts from peach fruit mesocarp at four different postharvest conditions. The post-harvest conditions include: fruits processed in a packing plant (E1: non-ripe; no long term cold storage); packing followed by a shelf-life at 20°C for 2-6 days (E2: Ripe; no long term cold storage; juicy fruits); packing followed by cold storage at 4°C for 21 days (E3: non-ripe; long term cold storage) and packing followed by cold storage at 4°C for 21 days and shelf-life at 20°C for 2-6 days (E4: Ripe; long term cold storage; woolly fruits).

As we described in Vizoso et al [17], the contigs that represent differentially expressed genes were identified using the Winflat program that submits the sequence data to a rigorous statistical analysis described by Audic and Claverie [18] http://igs-server.cnrs-mrs.fr. This analysis calculates the probability that a gene is equally expressed in two different conditions by observing the distribution of tag counts (number of ESTs). Therefore, small probability values (p-values) are associated with non-symmetrical distributions, characteristic of differentially expressed genes [18, 19].

To analyze the co-expression of differentially expressed genes, contigs were clustered using the Pearson linear correlation coefficient [19, 20]. Briefly, contigs with at least five ESTs were selected to make the expression profile matrix, which consisted of 1,402 rows (the contigs) and 4 columns (four cDNA libraries). The similarity between clusters and libraries was estimated using an un-centered Pearson's correlation coefficient in the Cluster 3.0 program [20] http://rana.lbl.gov/EisenSoftware.htm. Pearson correlation coefficients > 0.85 (zero values indicate no association and a coefficient equal to 1 indicate a fully correlated pattern) are indicated by an asterisk in Additional File 1. Dendrograms were constructed from the pair wise distances using the UPGMA algorithm. The results were visualized and analyzed using the Java TreeView program http://jtreeview.sourceforge.net.

Gene Ontology molecular function and biological process annotations of the contigs are described in Vizoso et al [17]. Each annotation and contig assembly was manually corrected, when necessary.

mRNA isolation and reverse transcriptase (RT)-PCR

The kit Oligotex™ mRNA Spin-Column (Qiagen, New York, USA) was used to purify mRNA. The mRNA was purified from pools of total RNA obtained from peach fruit mesocarp (O'Henry var.) representing the stages E1, E2, E3 and E4 as described previously [17, 21]. The mRNA was quantified using the Poly (A) mRNA Detection System™ (Promega, Madison, USA). First strand cDNA was synthesized from 5 ng of the mRNA in a 20 μl final volume. The reaction mix was prepared using the ImProm-II™ reverse Transcription System (Promega, Madison, USA) and anchored oligo (dT) of 18-mers, according to the manufacturer's instructions. As an internal control for normalization, heterologous mRNA (1.2 kb mRNA coding for Kanamycin) was added to each mRNA sample. To control for genomic DNA contamination, PCR amplification was performed on template RNA that was not reverse transcribed. To confirm that the amplified fragments correspond to the cDNAs of interest, these fragments were cloned in pBluescript and sequenced (Macrogen, Korea). The primer sequences used to amplify the internal regions of the basic endochitinase Ppbec1 (BEC226F and BEC576R), dehydrin Ppxero2 (DX-82F and DX176R), thaumatin Pptha1 (THA30F and THA382R), lipoxygenase Pplox1 (LOX982F and LOX1267R) and the actin Ppact7 (ACT-F and ACT-R) genes are shown in Table 1. Primers used to amplify a 323 bp fragment of the cDNA from the Kanamycin mRNA control are: "Upstream Control Primer" (5'-gCCATTCTCACCggATTCAgTCgTC-3') and "Downstream Control Primer" (5'-AgCCgCCgTCCCgTCAAgTCAg-3'). PCR reactions were performed by diluting the cDNAs a 100 fold and using 1 μl of each dilution as a template in a final reaction volume of 20 μl, containing 0.5 μM primers; 0.2 mM dNTPs; 1.5 mM MgCl2; 5U Taq polymerase and 1× buffer. The PCR conditions were: 93°C for 5 min and then a variable number of cycles (26 to 34) at 93°C for 30 sec, 1 min at 55°C, and 1 min at 72°C. The PCR reaction was with a final step at 72°C for 10 min.

Table 1 Primers used in this study

Cloning of the promoters

Genomic DNA was isolated from peach leaves (Prunus persica var. persica (L.) Batch cv. O'Henry) as described in Manubens et al [22]. The Universal Genome Walker™ Kit (Clontech Laboratories, Inc., Palo Alto, CA, USA) was used to isolate the promoters regions of Ppbec1, Ppxero2, Pptha1 and Pplox1. The isolated genomic DNA was digested with four restriction enzymes (EcoRV, PvuII, SspI, and MlsI). DNA fragments containing adaptors at both ends were used as a template for amplifying the promoter regions. GSP1 and GSP2 gene specific primers were designed to isolate the promoters (Table 1). For the first group of PCR reactions, a specific adaptor primer (AP1, 5'-ggATCCTAATACgACTCACTATAgggC-3') and the GSP1 primers specific for each gene were used. The final primer concentration in the PCR reaction was 0.2 μM in a final volume of 50 μL. Manual Hot Start was performed using 5 U of the Synergy DNA polymerase (Genecraft, Münster, Germany). The conditions for this first round of amplifications was: 1 cycle at 93°C for 10 min, 7 cycles of 93°C for 30 sec, 72°C for 15 min, followed by 37 cycles of 93°C for 30 sec, 67°C for 15 min. For the nested PCR, the specific adaptor primer 2 (AP2, 5'-ACTATAgggCACgCgTggT-3') and the gene specific GSP2 primers were used. As a DNA template in these reactions, 1 μL of a 50 fold dilution of end-product of the first round of amplifications was used. The conditions for the second round of amplification were: 1 cycle at 93°C for 10 min, 5 cycles (7 cycles in the case of Ppxero2) of 93°C for 30 sec, 72°C for 15 min, followed by 20 cycles (30 cycles in the case of Ppxero2) of 93°C for 30 sec, 67°C for 15 min. The amplified products were cloned in pGEM-T vector and sequenced (Macrogen, Korea). The Ppbec1 and Ppxero2 promoters were subsequently amplified from the pGEM-T clones using the AP2 and BEC-32BamHI or DX24BamHI primers, respectively (Table 1). The products of this amplification were also cloned in the pGEM-T vector and re-sequenced (Macrogen, Korea). The promoter fragments were extracted from the pGEM-T vector (including the Pptha1 promoter), with a BamHI-SalI sequential digestion, and transcriptionally fused to the uidA reporter gene in the promoterless binary vector pBI101.1 [23]. The binary vector was introduced into A. tumefaciens (GV3101) for subsequent Arabidopsis and peach fruit transformations.

Promoter sequences analysis

Analysis of putative transcription factor binding sites was carried out using the database PLACE http://www.dna.affrc.go.jp/htdocs/PLACE/ [24] coupled with visual analyses. To identify predicted conserved motifs, the promoter sequences were analyzed using the YMF 3.0 program [25] http://wingless.cs.washington.edu/YMF/YMFWeb/YMFInput.pl. Only the statistically significant motifs (Z score value > 6.5) were selected [26].

Growth, transformation and cold treatments of A. thaliana

Wild-type and transgenic A. thaliana (ecotype Columbia) were grown in a mixture of soil-vermiculite (3:1) in a growth chamber with a 16-h light cycle (140 μmol m-2 s-1) at 22°C. Alternatively, seeds were surface sterilized as described in Gonzalez et al [27], plated on Murashige-Skoog (1 × MS) media containing 0.8% agar, 0.1% sucrose and 50 mg/l Kanamycin for transgenic lines and grown under the same conditions as the soil-grown plants.

Transgenic Arabidopsis was obtained by using the GV3101 A. tumefaciens-mediated floral dip method [28]. A. tumefaciens previously transformed with the binary vector pBI101.3 harboring the promoter::uidA fusions: Ppbec1::uidA (PBIPpbec1); Ppxero2::uidA (pBIPpxero2); Pptha1::uidA (pBIPptha1), or the control vectors pBI121 (containing the 35S CaMV promoter) and pBI101.3 (promoterless), were used. In cold treatments, T3 homozygous transgenic Arabidopsis seedlings were grown on plates containing 1× MS media, 0.8% agar, and 0.1% sucrose in a growth chamber with a 16-h light cycle (140 μmol m-2 s-1) at 24°C for two weeks, and then transferred to 4°C for 7 days. A minimum of three independent transgenic lines were used for each construct.

Peach fruit transient transformation and cold treatments

A. tumefaciens transformed with the vectors pBIPpbec1, pBIPpxero2, pBIPptha1, pBI121 or pBI101.3 were grown in LB medium supplemented with Kanamycin (100 μg/ml), Rifampicin (10 μg/ml) and Gentamycin (100 μg/ml). The cultures were grown for two days at 28°C until they reached an OD600 between 0.6 and 0.8. The culture was then centrifuged and the pellet re-suspended in MMA medium (1× MS, MES 10 mM (pH 5.6), 20 g/l sucrose, and 200 μM acetosyringone) to reach an OD600 of 2.4. Approximately 0.7 mL of this bacterial suspension was used to infiltrate mature fruits from O'Henry, Elegant Lady and Florida King varieties of peach as described by Spolaore et al [15].

To analyze the promoter activity at 20°C, the fruits infiltrated with the different constructs, were stored in a dark growth chamber for five days. To analyze the cold-responsive promoter activity, the infiltrated fruits were stored 2 days post-infiltration (dpi) in a dark growth chamber at 4°C for 10 days. After the growth chamber incubation time, the infiltrated region of the fruit was extracted with a cork bore and stained for GUS activity as described by Tittarelli et al [14].

GUS activity measurement

Histochemical staining of Arabidopsis seedlings for β-glucuronidase (GUS) activity was performed as described by Jefferson et al [23], with the following modifications: transgenic Arabidopsis seedlings used in the cold-treatments described earlier were vacuum infiltrated in 50 mM NaH2PO4, pH 7.0; 0.1 mM X-Gluc; 10 mM EDTA and 0.1% Triton X-100. These samples were incubated in the dark at 37°C for 24-72 h. Samples that did not develop color after 72 h were considered negative for GUS activity. Plant material was subsequently fixed in 0.04% formaldehyde, 0.04% acetic acid and 0.285% ethanol for 30 min, followed by an ethanol dilution series to remove chlorophyll from the plant tissue (70% ethanol for 1 h, 100% ethanol for 1 h, 70% ethanol for 1 h and distilled water).

Slices (2 mm) of transiently transformed peaches were imbibed in the GUS staining solution (0.72 M K2HPO4; 0.17 M KH2PO4; 0.5 mM K3Fe(CN)6; 0.5 mM K4Fe(CN)6; 1× Triton X-100; 12.7 mM EDTA; 20% (v/v) methanol and 0.5 mM X-Gluc) [15]. Samples were vacuum-infiltrated for 30 min at room-temperature and then incubated overnight at 37°C. Fluorometric GUS assays were performed as described by Jefferson et al [23]. The Arabidopsis seedlings were ground in a mortar using liquid nitrogen, and the tissue powder was transferred to a microtube. One ml of the extraction buffer (50 mM NaH2PO4, pH 7.0; 1 mM EDTA; 0.1% Triton X-100; 0.1% (w/v) sodium laurylsarcosine and 5 mM dithiothreitol) was added. Samples were centrifuged for 10 min at 12,000 g at 4°C and the supernatant was transferred to a new microtube. The fluorogenic reaction was carried out in 2 ml volume containing 1 mM 4-methyl umbelliferyl glucuronide (MUG) in an extraction buffer supplemented with a 50 μL aliquot of the protein extract supernatants. The protein quantity of the sample extracts was determined as described previously [29], using bovine serum albumin (BSA) as a standard.

Results

Identification of peach cold-regulated genes by digital expression analyses of EST datasets

Coordinated gene expression analyses of peach fruit ESTs datasets revealed 10 major hierarchical clusters (Additional File 1), containing unique contigs. We identified 164 contigs with preferential expression in fruits stored at 4°C (E3: non-ripe; long term cold storage). Table 2 contains a complete list of these contigs together with their annotations, GO biological process annotations and the origin of the ESTs in each contig. Contigs with statistically differential expression, in E3 compared to the other stages are also indicated.

Table 2 Putative function of 164 genes preferentially expressed in cold stored peach fruits.

Approximately 95% of the 164 cold-induced peach genes share significant identity with sequences in Arabidopsis, suggesting that these may be putative orthologs. The putative Arabidopsis orthologs that are induced or repressed by cold, based on ColdArrayDB analyses http://cold.stanford.edu/cold/cgi-bin/data.cgi are shown in Table 2. Only 29 contigs (18% of the 164 cold-induced genes) share significant sequence identity with genes of unknown function. Approximately 38% of these contigs (11 contigs) share significant sequence identity with plant gene sequences annotated as expressed proteins. Six of the contigs with unknown function do not share sequence identity with any sequences in the public databases, suggesting that these are novel genes.

Annotation frequency comparative analyses of cold-induced (164 contigs), cold-repressed (138 contigs) or contigs unrelated to cold (1,238 contigs), revealed an overrepresentation of stress response genes and an underrepresentation of genes related to energy metabolism in fruits that were stored in the cold (Figure 1). Among the genes related to stress response we identified four contigs that are similar to thaumatin-like proteins: C1708, C2177, C2317 and C2147 (98%, 99%, 98% and 93% amino acid identity with P. persica thaumatin-like protein 1 precursor, respectively, GenBank accession number: P83332). Three of the stress response genes are similar to chitinases: C910 (76% amino acid identity with Malus domestica class III acidic endochitinase, GenBank accession number: ABC47924); C2131 (74% amino acid identity with Galega orientalis class Ib basic endochitinase, GenBank accession number: AAP03087) and C2441 (72% amino acid identity with A. thaliana class IV chitinase, GenBank accession number: NP_191010). Two of the stress response genes are similar to dehydrins: C254 (97% amino acid identity with P. persica Ppdhn1, GenBank accession number: AAC49658) and C304, 100% amino acid identity with P. persica type II SK2 dehydrin Ppdhn3 (Genbank accession number: AAZ83586).

Figure 1
figure1

Annotation frequency comparison of cold-induced, cold-repressed or unrelated to cold-induction contigs. The frequency of contigs that are associated with a specific Gene Ontology are expressed as the percentage of the total annotations for each analyzed group (164 for the cold-induced, 138 for the cold-repressed and 1,238 for unrelated to cold-induction). The numbers of contigs in each group, belonging to each biological process classification, are show at the top of each bar. The category "others process" are: cell adhesion (GO: 0007155, 1 contig); cell communication (GO: 0007154, 1 contig); cell cycle (GO: 0007049, 5 contigs); cell death (GO: 0008219, 1 contig); cell homeostasis (GO: 0019725, 4 contigs); organism physiological process (GO: 0050874; 1 contig); regulation of GTPase activity (GO: 0043087; 1 contig); response to stimulus (GO: 0050896; 10 contigs) and viral life cycle (GO: 0016032; 1 contig).

Cold-induced expression of Ppbec1, Ppxero2 and Pptha1

We evaluated the expression levels of three cold-induced candidate genes by RT-PCR: a basic endochitinase (C2131, Ppbec1), a dehydrin (C254, Ppxero2) and a thaumatin-like protein (C2317, Pptha1). These genes were chosen due to the high number of ESTs in cold-stored fruits (E3), as revealed by the digital expression analyses (Figure 2). The expression level of a contig similar to lipoxygenase (C3336, Pplox1) that does not express preferentially in cold stored fruits (E3) as well as the expression level of a contig (C407, Ppact7) that does not significantly change expression under the different post-harvest conditions, were analyzed (Figure 2). Interestingly, all five genes analyzed showed an expression pattern significantly similar to the ones predicted by the digital expression analyses (Figure 2). The genes Ppbec1, Ppxero2 and Pptha1 have an increased expression in cold-stored fruits, whereas the Pplox1 gene increased expression in woolly fruits rather than cold-stored fruits.

Figure 2
figure2

Evaluation of the accuracy of the predicted expression patterns of selected genes by RT-PCR. (A) RT-PCR analysis of RNA expression of three cold-induced genes: Ppbec1, Ppxero2, and Pptha1 under different post-harvest conditions. These post-harvest conditions include: fruits processed in a packing plant (E1: non-ripe; no long term cold storage); packing followed by a shelf-life at 20°C for 2-6 days (E2: Ripe; no long term cold storage; juicy fruits); packing followed by cold storage at 4°C for 21 days (E3: non-ripe; long term cold storage) and packing followed by cold storage at 4°C for 21 days and shelf-life at 20°C for 2-6 days (E4: Ripe; long term cold storage; woolly fruits). The expression level of Pplox1 was analyzed as a control for genes that do not express preferentially in cold stored fruits (E3). Ppact7 was analyzed as a control for genes that do not significantly change expression levels between the four post-harvest conditions analyzed. The two arrows associated with each gel represent 500 bp (upper) and 300 bp (lower). The number of ESTs associated with each contig and library source is indicated. (B) Densitometry quantification of the expression level obtained by RT-PCR, the figure shows the bands intensities for each gene relative to Ppact7 intensity.

Identification of conserved motifs in the promoters of cold-inducible genes Ppbec1, Ppxero2 and Pptha1

We cloned 826 bp, 1,348 bp and 1,559 bp fragments corresponding to the regions upstream of the translation start codons of Ppbec1, Ppxero2 and Pptha1, respectively. The sequences of these promoter regions as well as the cDNA of their corresponding genes are shown in the Additional Files 2, 3 and 4.

The high sequence identity between the Ppxero2 contig with the coding region of Ppdhn1[30] was also observed within the promoter sequences of these two genes. Only one nucleotide difference at position -469 was found, suggesting that Ppxero2 and Ppdhn1 may be the same gene (Additional File 3). However, the promoter isolated in this work is about 230 bp longer (at the 5' end) than the previously published promoter [30].

Cis-element regulatory motifs related to cold gene expression regulation such as ABRE [13], MYCR [31, 32], MYBR [31, 33] and DRE/CRT [34] were identified in all three promoters of these cold-inducible genes (Figure 3). In addition, three statistically significant predicted motifs were present in the promoters of these cold-inducible genes (TACGTSGS, TGTGTGYS and CTAGAASY (Figure 3). These motifs were not found in the Pplox1 promoter identified in this work (Additional File 5).

Figure 3
figure3

Putative cis -regulatory elements identified in Ppbec1 , Ppxero2 and Pptha1 promoter sequences. Topologies of the Ppbec1 (A), Ppxero2 (B) and Pptha1 (C) promoters are shown. The promoters are draw proportionally (the bar correspond to 100 bp). Boxed regions: predicted 5' UTR region. Black arrow shows the position of different cis-regulatory elements related to low temperature responses: ABRE, DRE/CRT, MYBR and MYCR. The putative cis-regulatory elements identified by the motif prediction program YMF3.0 are shown as grey triangle, black circle and asterisk. The sequences, the symbol and the significance score (Zscore) of the motifs, are shown in the upper left corner. The degenerate bases allowed in the motifs are S (C or G) and Y (C or T). Note: in order to ensure at the legibility of the figure, not all cis-elements are marked in (B) and (C). However, the complete sequences of these promoters are available in Additional Files 3 and 4.

Cold-induced Ppbec1 and Ppxero2 promoters in transiently transformed peach fruits and stably transformed Arabidopsis

Transient transformation assays of peach fruits revealed that all three cloned promoters (pBIPpbec1, pBIPxero2 and pBIPptha1) were able to activate GUS (uidA) expression (Figure 4). However, only the pBIPpbec1 and pBIPxero2 promoter constructs showed cold-inducible increases in GUS activity (Figure 4). The pBIPtha1 construct was expressed at both 20°C and 4°C. Comparable results were seen in fruits from three different peach varieties (data not shown).

Figure 4
figure4

Cold-inducible peach Ppbec1 and Ppxero2 promoters in transiently transformed peach fruits. (A) Structure of the binary vector constructs used for functional analysis of the Ppbec1, Ppxero2 and Pptha1 promoter-uidA fusions. LB and RB: left and right T-DNA border. (B) Histochemical GUS staining of fruit slices from agro-infiltrated peaches stored at 20°C for 5 days post-inoculation or 4°C for 10 days. These images correspond to the transient transformation of O'Henry variety fruits. However, similar results were seen in all varieties assayed (data not shown).

Similar results were seen when these promoter-GUS constructs were analyzed in stably transformed Arabidopsis. All three constructs were able to activate GUS expression, but only the Ppbec1 and Ppxero2 promoters (pBIPpbec1 and pBIPxero2, respectively) induced expression in response to cold (Figure 5). As observed with the fruit transient transformation assays, the Pptha1 promoter (pBIPtha1) expressed GUS under all conditions analyzed.

Figure 5
figure5

Conserved heterologous regulation of the cold-inducible peach Ppbec1 and Ppxero2 promoters in transgenic Arabidopsis plants. The upper panel shows histochemical GUS staining of representative transgenic Arabidopsis lines carrying the Ppbec1 promoter-uidA fusion, Ppxero2 promoter-uidA fusion and Pptha1 promoter-uidA fusion. The lower panel shows the results of fluorometric GUS-assays of three independent Arabidopsis transgenic lines (L1, L2 and L3) containing the Ppbec1 promoter-uidA fusion, Ppxero2 promoter-uidA fusion or Pptha1 promoter-uidA fusion. Homozygous T3 plants were grown for 14 days in MS plates with 0.8% agar at 24°C (white bars) and then transfer to 4°C for 7 days (blacks bars). The asterisk above each bar represents those samples that have a statistically significant increase in GUS activity in the cold treated plants when compared to the untreated plants. Bars represent the mean ± standard deviation, n = 5. t-student * p < 0.01.

Discussion and Conclusion

Digital expression analyses of EST datasets have permitted us to identify a large diversity of cold-inducible genes in peach fruits, three of which were chosen for further analyses (Ppbec1, Ppxero2 y Pptha1). Both digital expression analyses and RT-PCR suggest that the Ppbec1, Ppxero2 and Pptha1 are cold-inducible genes. The promoters of these cold-inducible genes were isolated and characterized using both transient transformation assays in peach fruits and stable transformation in Arabidopsis. These analyses have revealed that the isolated Ppbec1 and Ppxero2 promoters are cold-inducible promoters, whereas the isolated Pptha1 promoter was not cold-inducible. These results, therefore, demonstrate that the isolated Ppbec1 and Ppxero2 promoters are sufficient for cold-induced gene expression. Furthermore, these results suggest that there is a conserved heterologous cold-inducible regulation of these promoters in peach and Arabidopsis.

Plants respond to cold temperatures by modifying the transcription and translation levels of hundreds of genes [35, 36]. These acute molecular changes are related to plant cell physiological and biochemical modifications (cold acclimation) that lead to stress tolerance and cold adaptation (a chronic response). In peach fruits, cold temperatures induce chilling injury, possibly due to global transcriptome changes [37]. With the exception of studies in the model organism A. thaliana [4] and work published recently [17, 38], little is known about the peach global transcriptional response to cold. Using the Pearson correlation coefficient, we analyze the coordinated gene expression of 1,402 contigs. This analysis revealed 164 genes preferentially expressed in peach fruits, of which digital expression analyses [18] revealed 45 of these genes (27%) with statistically significant cold-induction. A large proportion of the contigs preferentially expressed at 4°C (around 74% of the total) do not exhibited significant sequence homology (e-value < e-10) with the rest of the analyzed contigs (Table 2). This result could suggest that these contigs represent genes with non-redundant functions that will have a special importance during the exposure of the fruits to low temperatures.

Among the highly expressed genes in cold stored fruits, we found genes related to stress response in plants, including three dehydrins (C30, C254 and C304), three chitinases (C910, C2131 and C2441), four thaumatin-like proteins (C1708, C2177, C2317 and C2147), and polygalacturonase inhibiting protein (C2988), similar to what was reported by Ogundiwin et al [38]. Dehydrins are hydrophilic proteins that belong to the subgroup D-11 of the LEA ("late-embryogenesis-abundant") proteins [39]. There is some evidence that suggests that dehydrins protect macromolecules such as membranes and proteins against the damages associated with water deficiency [4042]. In peach, these genes are induced during cold acclimation and in cold-stored fruits [30, 38]. It has been observed that pathogenesis-related (PR) proteins such as chitinases and thaumatins are accumulated in the apoplastic space in winter rye during cold acclimation. These proteins also may have antifreeze properties that will protect the integrity of the plant cell avoiding the formation of ice [43, 44]. It has also been observed that these types of proteins retain their enzymatic activity under low temperatures, and may form part of a general response mechanism associated with unfavorable conditions, by providing protection from opportunist pathogen attack whilst the plant is in a weakened state [4547]. A similar role is shared by polygalacturonase inhibiting proteins in different plants models [48, 49].

We also found some genes related to protein folding and degradation, such as heat shock proteins, BiP-1 and DJ-1 family proteins (Table 2). These processes are very active when plants face low temperatures, chemical and oxidative stress. These proteins participate in the prevention and repair of damage produced by cold, through the stabilization of protein structure and the degradation of proteins that are not folded correctly [50, 51].

In this work we were interested in isolating and functionally characterizing promoters of cold-inducible peach genes. To date, only a few inducible promoters have been identified in crop plants. The Pptha1, Ppbec1 and Ppxero2 genes were chosen for promoter cloning and characterization based on the up-regulation that these genes showed in the in silico analysis and RT-PCR. The promoter sequences of these genes contain several cis-regulatory elements such as DRE/CRT, ABRE, MYCR (MYC recognition site) and MYBR (MYB recognition site) [13, 3134] that are related to stress response, specifically to cold/dehydration. These cis-regulatory elements are conserved in several plant species [52]. The presence of these conserved motifs suggests that these promoters may respond to the cold. Using transient transformation in peach fruit we confirmed that the promoters isolated from Ppbec1 and Ppxero2 are induced during low temperature storage, but not at room temperature. On the other hand, the Pptha1 promoter is active under all the temperatures analyzed. This could indicate that the Pptha1 promoter sequence might not contain all the elements needed to regulate expression in a cold-inducible manner. Alternatively, the agro-infiltration technique may induce stress signals that will activate this promoter. However, this last possibility is not likely because the activation of the Pptha1 promoter at all analyzed temperatures is also seen in the stably transformed transgenic Arabidopsis plants. The promoters Ppbec1 and Ppxero2, however, are cold-induced both in Arabidopsis transgenic plants as well as transient expressing fruits, suggesting that the Ppbec1 and Ppxero2 promoters are cold-inducible peach promoters. The cold-inducibility of these promoters in A. thaliana also suggests that this model plant may be used to functionally analyze peach cold-induced genes as well as their corresponding cis-elements and trans-acting factors.

The identification of these fruit tree cold-inducible promoters as well as the conserved heterologous regulation of these promoters in peach and Arabidopsis, demonstrates that these two transformation assays may be used to molecularly define the cis-elements and trans-acting regulatory factors that are associated with cold-responsive genes. By better understanding the regulatory mechanisms associated with cold-responsive genes, we may better understand the molecular differences and similarities between cold acclimation and chilling injury as well as the role these processes play in fruit tree growth and fruit quality.

References

  1. 1.

    Thomashow MF: So what's new in the field of plant cold acclimation? Lots!. Plant Physiol. 2001, 125 (1): 89-93. 10.1104/pp.125.1.89.

  2. 2.

    Sharma P, Sharma N, Deswal R: The molecular biology of the low-temperature response in plants. Bioessays. 2005, 27 (10): 1048-1059. 10.1002/bies.20307.

  3. 3.

    Thomashow MF: PLANT COLD ACCLIMATION: Freezing Tolerance Genes and Regulatory Mechanisms. Annu Rev Plant Physiol Plant Mol Biol. 1999, 50: 571-599. 10.1146/annurev.arplant.50.1.571.

  4. 4.

    Zhu J, Dong CH, Zhu JK: Interplay between cold-responsive gene regulation, metabolism and RNA processing during plant cold acclimation. Curr Opin Plant Biol. 2007, 10 (3): 290-295. 10.1016/j.pbi.2007.04.010.

  5. 5.

    Fowler S, Thomashow MF: Arabidopsis transcriptome profiling indicates that multiple regulatory pathways are activated during cold acclimation in addition to the CBF cold response pathway. Plant Cell. 2002, 14 (8): 1675-1690. 10.1105/tpc.003483.

  6. 6.

    Mahajan S, Tuteja N: Cold, salinity and drought stresses: an overview. Arch Biochem Biophys. 2005, 444 (2): 139-158. 10.1016/j.abb.2005.10.018.

  7. 7.

    Crisosto C, Mitchell F, Ju Z: Susceptibility to chilling injury of peach, nectarine, and plum cultivars grown in California. HortScience. 1999, 34: 1116-1118.

  8. 8.

    Crisosto CH, Mitchell FG, Johnson S: Factors in fresh market stone fruit quality. Postharvest News and Information. 1995, 6 (2): 17-21.

  9. 9.

    Retamales J, Cooper T, Streif J, Kania J: Preventing cold storage disorders in nectarines. J Hort Sci. 1992, 67: 619-626.

  10. 10.

    Bonetta D, McCourt P: Genetic analysis of ABA signal transduction pathways. Trends in Plant Science. 1998, 3 (6): 231-235. 10.1016/S1360-1385(98)01241-2.

  11. 11.

    Kasuga M, Liu Q, Miura S, Yamaguchi-Shinozaki K, Shinozaki K: Improving plant drought, salt, and freezing tolerance by gene transfer of a single stress-inducible transcription factor. Nat Biotechnol. 1999, 17 (3): 287-291. 10.1038/7036.

  12. 12.

    Kyonoshin Maruyama YS, Mie Kasuga, Yusuke Ito, Motoaki Seki, Hideki Goda, Yukihisa Shimada, Shigeo Yoshida, Kazuo Shinozaki, Kazuko Yamaguchi-Shinozaki: Identification of cold-inducible downstream genes of the Arabidopsis DREB1A/CBF3 transcriptional factor using two microarray systems. The Plant Journal. 2004, 38 (6): 982-993. 10.1111/j.1365-313X.2004.02100.x.

  13. 13.

    Yamaguchi-Shinozaki K, Shinozaki K: A novel cis-acting element in an Arabidopsis gene is involved in responsiveness to drought, low-temperature, or high-salt stress. Plant Cell. 1994, 6 (2): 251-264. 10.1105/tpc.6.2.251.

  14. 14.

    Tittarelli A, Milla L, Vargas F, Morales A, Neupert C, Meisel LA, Salvo GH, Penaloza E, Munoz G, Corcuera LJ, et al: Isolation and comparative analysis of the wheat TaPT2 promoter: identification in silico of new putative regulatory motifs conserved between monocots and dicots. J Exp Bot. 2007, 58 (10): 2573-2582. 10.1093/jxb/erm123.

  15. 15.

    Spolaore S, Trainotti L, Casadoro G: A simple protocol for transient gene expression in ripe fleshy fruit mediated by Agrobacterium. J Exp Bot. 2001, 52 (357): 845-850.

  16. 16.

    Ahmad M, Mirza B: An efficient protocol for transient transformation of intact fruit and transgene expression in Citrus. Plant Molecular Biology Reporter. 2005, 23 (4): 419-420. 10.1007/BF02788891.

  17. 17.

    Vizoso P, Meisel L, Tittarelli A, Latorre M, Saba J, Caroca R, Maldonado J, Cambiazo V, Campos-Vargas R, Gonzalez M, et al: Comparative EST transcript profiling of peach fruits under different postharvest conditions reveals candidate genes associated with peach fruit quality. BMC Genomics. 2009, 10: 421-10.1186/1471-2164-10-423.

  18. 18.

    Audic S, Claverie JM: The significance of digital gene expression profiles. Genome Res. 1997, 7 (10): 986-995.

  19. 19.

    Ewing RM, Ben Kahla A, Poirot O, Lopez F, Audic S, Claverie JM: Large-scale statistical analyses of rice ESTs reveal correlated patterns of gene expression. Genome Res. 1999, 9 (10): 950-959. 10.1101/gr.9.10.950.

  20. 20.

    Eisen MB, Spellman PT, Brown PO, Botstein D: Cluster analysis and display of genome-wide expression patterns. Proc Natl Acad Sci USA. 1998, 95 (25): 14863-14868. 10.1073/pnas.95.25.14863.

  21. 21.

    Meisel L, Fonseca B, Gonzalez S, Baeza-Yates R, Cambiazo V, Campos R, Gonzalez M, Orellana A, Retamales J, Silva H: A rapid and efficient method for purifying high quality total RNA from peaches (Prunus persica) for functional genomics analyses. Biol Res. 2005, 38 (1): 83-88. 10.4067/S0716-97602005000100010.

  22. 22.

    Manubens A, Lobos S, Jadue Y, Toro M, Messina R, Lladser M, Seelenfreund D: DNA Isolation and AFLP Fingerprinting of Nectarine and Peach Varieties (Prunus persica). Plant Molecular Biology Reporter. 1999, 17 (3): 255-267. 10.1023/A:1007656110444.

  23. 23.

    Jefferson RA, Kavanagh TA, Bevan MW: GUS fusions: beta-glucuronidase as a sensitive and versatile gene fusion marker in higher plants. EMBO J. 1987, 6 (13): 3901-3907.

  24. 24.

    Higo K, Ugawa Y, Iwamoto M, Korenaga T: Plant cis-acting regulatory DNA elements (PLACE) database: 1999. Nucleic Acids Res. 1999, 27 (1): 297-300. 10.1093/nar/27.1.297.

  25. 25.

    Blanchette M, Sinha S: Separating real motifs from their artifacts. Bioinformatics. 2001, 17 (Suppl 1): S30-38.

  26. 26.

    Sinha S, Tompa M: Discovery of novel transcription factor binding sites by statistical overrepresentation. Nucleic Acids Res. 2002, 30 (24): 5549-5560. 10.1093/nar/gkf669.

  27. 27.

    Gonzalez J, Reyes F, Salas C, Santiag M, Codriansky Y, Coliheuque N, Silva H: Arabidopsis thaliana: a model host plant to study plant-pathogen interaction using Chilean field isolates of Botrytis cinerea. Biol Res. 2006, 39 (2): 221-228. 10.4067/S0716-97602006000200004.

  28. 28.

    Clough SJ, Bent AF: Floral dip: a simplified method for Agrobacterium-mediated transformation of Arabidopsis thaliana. Plant J. 1998, 16 (6): 735-743. 10.1046/j.1365-313x.1998.00343.x.

  29. 29.

    Bradford MM: A rapid and sensitive method for the quantitation of microgram quantities of protein utilizing the principle of protein-dye binding. Anal Biochem. 1976, 72: 248-254. 10.1016/0003-2697(76)90527-3.

  30. 30.

    Wisniewski ME, Bassett CL, Renaut J, Farrell R jr, Tworkoski T, Artlip TS: Differential regulation of two dehydrin genes from peach (Prunus persica) by photoperiod, low temperature and water deficit. Tree Physiol. 2006, 26 (5): 575-584.

  31. 31.

    Abe H, Urao T, Ito T, Seki M, Shinozaki K, Yamaguchi-Shinozaki K: Arabidopsis AtMYC2 (bHLH) and AtMYB2 (MYB) function as transcriptional activators in abscisic acid signaling. Plant Cell. 2003, 15 (1): 63-78. 10.1105/tpc.006130.

  32. 32.

    Chinnusamy V, Ohta M, Kanrar S, Lee BH, Hong X, Agarwal M, Zhu JK: ICE1: a regulator of cold-induced transcriptome and freezing tolerance in Arabidopsis. Genes Dev. 2003, 17 (8): 1043-1054. 10.1101/gad.1077503.

  33. 33.

    Solano R, Nieto C, Avila J, Canas L, Diaz I, Paz-Ares J: Dual DNA binding specificity of a petal epidermis-specific MYB transcription factor (MYB.Ph3) from Petunia hybrida. EMBO J. 1995, 14 (8): 1773-1784.

  34. 34.

    Shinozaki K, Yamaguchi-Shinozaki K: Molecular responses to dehydration and low temperature: differences and cross-talk between two stress signaling pathways. Curr Opin Plant Biol. 2000, 3 (3): 217-223.

  35. 35.

    Umezawa T, Fujita M, Fujita Y, Yamaguchi-Shinozaki K, Shinozaki K: Engineering drought tolerance in plants: discovering and tailoring genes to unlock the future. Curr Opin Biotechnol. 2006, 17 (2): 113-122.

  36. 36.

    Yamaguchi-Shinozaki K, Shinozaki K: Organization of cis-acting regulatory elements in osmotic- and cold-stress-responsive promoters. Trends Plant Sci. 2005, 10 (2): 88-94. 10.1016/j.tplants.2004.12.012.

  37. 37.

    Gonzalez-Aguero M, Pavez L, Ibanez F, Pacheco I, Campos-Vargas R, Meisel LA, Orellana A, Retamales J, Silva H, Gonzalez M, et al: Identification of woolliness response genes in peach fruit after post-harvest treatments. J Exp Bot. 2008, 59 (8): 1973-1986. 10.1093/jxb/ern069.

  38. 38.

    Ogundiwin EA, Marti C, Forment J, Pons C, Granell A, Gradziel TM, Peace CP, Crisosto CH: Development of ChillPeach genomic tools and identification of cold-responsive genes in peach fruit. Plant Mol Biol. 2008, 68 (4-5): 379-397. 10.1007/s11103-008-9378-5.

  39. 39.

    Close TJ: Dehydrins: A commonalty in the response of plants to dehydration and low temperature. Physiologia Plantarum. 1997, 100: 291-296. 10.1111/j.1399-3054.1997.tb04785.x.

  40. 40.

    Koag MC, Fenton RD, Wilkens S, Close TJ: The binding of maize DHN1 to lipid vesicles. Gain of structure and lipid specificity. Plant Physiol. 2003, 131 (1): 309-316. 10.1104/pp.011171.

  41. 41.

    Rorat T: Plant dehydrins--tissue location, structure and function. Cell Mol Biol Lett. 2006, 11 (4): 536-556. 10.2478/s11658-006-0044-0.

  42. 42.

    Sanchez-Ballesta MT, Rodrigo MJ, Lafuente MT, Granell A, Zacarias L: Dehydrin from citrus, which confers in vitro dehydration and freezing protection activity, is constitutive and highly expressed in the flavedo of fruit but responsive to cold and water stress in leaves. J Agric Food Chem. 2004, 52 (7): 1950-1957. 10.1021/jf035216+.

  43. 43.

    Hon WC, Griffith M, Mlynarz A, Kwok YC, Yang DS: Antifreeze proteins in winter rye are similar to pathogenesis-related proteins. Plant Physiol. 1995, 109 (3): 879-889. 10.1104/pp.109.3.879.

  44. 44.

    Yu XM, Griffith M, Wiseman SB: Ethylene induces antifreeze activity in winter rye leaves. Plant Physiol. 2001, 126 (3): 1232-1240. 10.1104/pp.126.3.1232.

  45. 45.

    Dani V, Simon WJ, Duranti M, Croy RR: Changes in the tobacco leaf apoplast proteome in response to salt stress. Proteomics. 2005, 5 (3): 737-745. 10.1002/pmic.200401119.

  46. 46.

    Diaz-Vivancos P, Rubio M, Mesonero V, Periago PM, Barcelo AR, Martinez-Gomez P, Hernandez JA: The apoplastic antioxidant system in Prunus: response to long-term plum pox virus infection. J Exp Bot. 2006, 57 (14): 3813-3824. 10.1093/jxb/erl138.

  47. 47.

    Hiilovaara-Teijo M, Hannukkala A, Griffith M, Yu XM, Pihakaski-Maunsbach K: Snow-mold-induced apoplastic proteins in winter rye leaves lack antifreeze activity. Plant Physiol. 1999, 121 (2): 665-674. 10.1104/pp.121.2.665.

  48. 48.

    Meyer K, Keil M, Naldrett MJ: A leucine-rich repeat protein of carrot that exhibits antifreeze activity. FEBS Lett. 1999, 447 (2-3): 171-178. 10.1016/S0014-5793(99)00280-X.

  49. 49.

    Yao C, Conway WS, Ren R, Smith D, Ross GS, Sams CE: Gene encoding polygalacturonase inhibitor in apple fruit is developmentally regulated and activated by wounding and fungal infection. Plant Mol Biol. 1999, 39 (6): 1231-1241. 10.1023/A:1006155723059.

  50. 50.

    Iba K: Acclimative response to temperature stress in higher plants: approaches of gene engineering for temperature tolerance. Annu Rev Plant Biol. 2002, 53: 225-245. 10.1146/annurev.arplant.53.100201.160729.

  51. 51.

    Sabehat A, Lurie S, Weiss D: Expression of small heat-shock proteins at low temperatures. A possible role in protecting against chilling injuries. Plant Physiol. 1998, 117 (2): 651-658. 10.1104/pp.117.2.651.

  52. 52.

    Chen WJ, Zhu T: Networks of transcription factors with roles in environmental stress response. Trends Plant Sci. 2004, 9 (12): 591-596. 10.1016/j.tplants.2004.10.007.

Download references

Acknowledgements

This work was supported by ICM P06-065-F; FDI G02P1001 (Chilean Genome Initiative) with funding from the Chilean government as well as ASOEX (Asociación de Exportadores de Chile A.G.), FDF (Fundación para el Desarrollo Frutícola) and Fundación Chile; Proyecto Consorcio BIOFRUTALES S.A.; PBCT R11 and CONICYT Fellowship D-21080654 to AM.

Author information

Correspondence to Herman Silva.

Additional information

Authors' contributions

AT: identified and cloned the promoters. AT, MS, LM and HS drafted the manuscript. AT and MS: performed the digital expression analysis. AM and AT: performed the construction of Arabidopsis transgenic plants as well as the transient assay. HS: conceived, supervised and participated in all the analysis. All authors read and approved the manuscript.

Electronic supplementary material

Additional file 1: Identification of fruit cold-induced contigs using correlated expression analysis of peach ESTs. The data provided represents the co-expression analysis of differentially expressed genes. The contigs were clustered using the Pearson linear correlation coefficient. (DOC 140 KB)

Additional file 2: Sequence of the Ppbec1 promoter and open reading frame. The data provided represents the sequences of the Ppbec1 promoter and open reading frame. (DOC 784 KB)

Additional file 3: Sequence of the Ppxero2 promoter and open reading frame. The data provided represents the sequences of the Ppxero2 promoter and open reading frame. (DOC 1 MB)

Additional file 4: Sequence of the Pptha1 promoter and open reading frame. The data provided represents the sequences of the Pptha1 promoter and open reading frame. (DOC 864 KB)

Additional file 5: Sequence of the Pplox1 promoter and open reading frame. The data provided represents the sequences of the Pplox1 promoter and open reading frame. (DOC 1 MB)

Authors’ original submitted files for images

Below are the links to the authors’ original submitted files for images.

Authors’ original file for figure 1

Authors’ original file for figure 2

Authors’ original file for figure 3

Authors’ original file for figure 4

Authors’ original file for figure 5

Rights and permissions

Reprints and Permissions

About this article

Keywords

  • Cold Storage
  • Cold Acclimation
  • Peach Fruit
  • Transgenic Arabidopsis Seedling
  • Chilling Injury Symptom