Temperature stress differentially modulates transcription in meiotic anthers of heat-tolerant and heat-sensitive tomato plants
© Bita et al; licensee BioMed Central Ltd. 2011
Received: 7 February 2011
Accepted: 31 July 2011
Published: 31 July 2011
Fluctuations in temperature occur naturally during plant growth and reproduction. However, in the hot summers this variation may become stressful and damaging for the molecular mechanisms involved in proper cell growth, impairing thus plant development and particularly fruit-set in many crop plants. Tolerance to such a stress can be achieved by constitutive gene expression or by rapid changes in gene expression, which ultimately leads to protection against thermal damage. We have used cDNA-AFLP and microarray analyses to compare the early response of the tomato meiotic anther transcriptome to moderate heat stress conditions (32°C) in a heat-tolerant and a heat-sensitive tomato genotype. In the light of the expected global temperature increases, elucidating such protective mechanisms and identifying candidate tolerance genes can be used to improve breeding strategies for crop tolerance to heat stress.
The cDNA-AFLP analysis shows that 30 h of moderate heat stress (MHS) alter the expression of approximately 1% of the studied transcript-derived fragments in a heat-sensitive genotype. The major effect is gene down-regulation after the first 2 h of stress. The microarray analysis subsequently applied to elucidate early responses of a heat-tolerant and a heat-sensitive tomato genotype, also shows about 1% of the genes having significant changes in expression after the 2 h of stress. The tolerant genotype not only reacts with moderate transcriptomic changes but also exhibits constitutively higher expression levels of genes involved in protection and thermotolerance.
In contrast to the heat-sensitive genotype, the heat-tolerant genotype exhibits moderate transcriptional changes under moderate heat stress. Moreover, the heat-tolerant genotype also shows a different constitutive gene expression profile compared to the heat-sensitive genotype, indicating genetic differences in adaptation to increased temperatures. In the heat-tolerant genotype, the majority of changes in gene expression is represented by up-regulation, while in the heat-sensitive genotype there is a general trend to down-regulate gene expression upon MHS. The putative functions associated with the genes identified by cDNA-AFLP or microarray indicate the involvement of heat shock, metabolism, antioxidant and development pathways. Based on the observed differences in response to MHS and on literature sources, we identified a number of candidate transcripts involved in heat-tolerance.
Abiotic stresses such as extreme temperatures, drought, flooding or chemical toxicity, pose serious threats to agricultural production. A rapid adaptation or an innate tolerance mechanism can protected the further development of the plant and importantly for yield, secure successful fruit-set. Gamete development in angiosperms takes place within two floral organs; the male stamen and the female pistil [1, 2]. The male gamete development starts with the differentiation of the reproductive tissues of the anther. After meiosis of the pollen mother cell and mitotic divisions, microspore maturation follows, resulting in the mature pollen grain. After initiation, highly specialized anther tissues will acquire non-reproductive (e.g., the tapetum for support, stomium for dehiscence) or reproductive functions (pollen - formation). Both tapetum and microspore development are essential for male fertility as documented by numerous studies on male sterile mutants [3–9]. In Petunia, a cDNA-AFLP study revealed that during meiosis in anthers, under standard conditions, besides the typical meiotic genes, other genes are also important for the development of the pollen grains. Thus, during meiosis, tapetum differentiation genes, serine-proteinases, hormone metabolism genes, genes involved in cell wall biosynthesis and even ribonucleases and polyamine biosynthesis are all modulated in expression . These genes represent members of physiological and metabolic pathways that are naturally differentially expressed during pollen biogenesis under normal conditions.
Tomato (Solanum lycopersicum L.) displays marked responses to heat, similar to other crop species including pepper, potato, melon, cowpea, wheat, common bean, rice and barley [11–13]. Hot summers in many agricultural regions can negatively affect the vegetative and reproductive growth phases of such crops  and can result in up to 70% losses in tomato yield . However, heat stress has numerous specific effects depending on the genotype. Physiological observations both under field and greenhouse conditions show a variable degree of tolerance between different genotypes. A wide range of heat stress phenotypes has been described [16, 17]. For example, of five tomato cultivars grown under moderate heat stress (MHS) conditions (32°C day and 26°C night), only one set fruit . The differences in pollen grain development among the tolerant genotypes are most critical factors to determine fruit set under heat stress. Comparing the effects of heat stress on a heat-tolerant and a heat-sensitive tomato cultivar showed that temperature stress affected mainly the development of pollen grains, where reduced viability was more pronounced in the heat-sensitive cultivar .
Under heat stress, it is the pollen grain development and particularly meiosis that shows the highest susceptibility, followed by germination and pollen tube growth, which in more severe cases is also significantly affected [20–22]. Recently, an analysis of maturing tomato microspores (mitotic anthers) exposed to heat-shock conditions has shown no differential gene expression between heat-tolerant and heat-sensitive genotypes but it seems that the capacity for thermotolerance may be achieved by modulating the expression levels of such 'responsive' genes prior to heat stress exposure .
In this paper we present an analysis of the different changes in gene expression in developing meiotic anthers (the first and most sensitive anther developmental stage) in response to MHS of tolerant and sensitive tomato genotypes. Assuming the protective mechanisms are initiated shortly after the commencement of the stress period, we decided to study gene expression changes in the first hours of MHS. Both cDNA-AFLP and Combimatrix microarray technology were applied to obtain a general overview of molecular mechanisms that participate in the response to MHS of anthers from a heat-tolerant and two heat-sensitive genotypes of tomato. We describe a set of candidate genes and pathways that open up the possibility of understanding and modulating male heat tolerance in tomato and other commercially important Solanaceae crops.
Results and discussion
In order to analyse the dynamics of transcriptional responses to moderate heat stress (MHS) in meiotic anthers of tolerant and sensitive genotypes, the temperature range for the experiment was chosen based on agronomically relevant temperatures shown to have a significant effect on pollen grain viability , rather than using classical heat shock conditions of 42°-45°C . We focused our analysis of gene expression on whole meiotic anther cones isolated from flower buds of plants that were exposed to MHS (32°C/26°C, day/night) for up to 30 h, and compared this to standard temperatures (26°C/18°C day/night). To analyse the heat stress response we used several genotypes characterised by plant breeders as relatively tolerant and sensitive.
Identification of differentially expressed transcripts during MHS by cDNA-AFLP
The most marked transcriptional responses for the heat-sensitive genotype were found at 2 h after initiating MHS: 96 TDFs showed a clear differential expression pattern. Of these, 41 showed a down-regulated pattern, 26 TDFs were transiently induced, 16 were up-regulated and 13 were transiently repressed. No further significant changes in gene expression were found after 6, 16 or 30 h of MHS. As a next step, a selection of TDFs exhibiting differential gene expression were sequenced in order to gain insight into the processes involved in the response(s) to MHS.
Characterisation of differentially expressed TDFs in MM from the cDNA-AFLP analysis
TC number and homology
e-value and putative function
similar to DNA-binding PD1-like protein; Pisum sativum, partial (35%)
6.2 e-36 metabolism
similar to 70 kD heat shock protein; Arabidopsis thaliana partial (25%)
heat shock response
homologue to cytochrome c oxidase subunit 1; Amoebidium
parasiticum partial (3%)
1.1 e-13 oxidation reduction
homologue to calcium/calmodulin-dependent protein kinase CaMK3;
Nicotiana tabacum, partial (47%)
heat shock response
similar to putative peptide transporter; Arabidopsis thaliana, partial (13%)
7.9 e-81 transport
similar to amino acid transporter-like protein; Arabidopsis
thaliana, partial (57%)
2.7 e-32 transport
similar to 4Fe-4S ferredoxin iron-sulfur binding domain protein;
Alkalilimnicola ehrlichel, partial (7%)
9.0 e-0.7 transport
Phosphosulfolactate synthase-related protein; Solanum lycopersicum, complete
heat shock response
similar to OSJNBb0015D13.13 protein; Oryza sativa, partial (87%)
Metabolism, cell wall
homologue to 40S ribosomal protein S8-like protein; Solanum
2.6 e-28 translation
weakly similar to TGF-beta receptor type l/lI extracellular region;
Medicago truncatula, partial (8%)
2.4 e-41 transport
similar to uncharacterized protein; Vitis vinifera partial (44%)
8.8 e-12 unknown
similar to Os02g0159700 protein; Oryza sativa Japonica
Group, partial (44%)
8.6 e-08 oxidation reduction
similar to P70 protein; Nicotiana tabacum
heat shock response
Similar to cytosolic ascorbate peroxidase 1; Solanum
lycopersicum, partial (32%)
1.1 e-19 oxidation reduction
similar to 3-deoxy-D-arabino-heptulosonate 7-phosphate synthase;
Vitis vinifera, partial (66%)
homologue to succinyl-CoA ligase alpha 1 subunit; Solanum
lycopersicum, partial (7%)
similar to DnaJ protein; Arabidopsis thaliana, partial (24%)
heat shock response
weakly similar to RNA polymerase Rpb4; Medicago truncatula, partial (52%)
2.3 e-38 translation
homologue to ribosomal protein S6-like protein; Solanum
3.6 e-41 Translation
similar to elicitor-inducible protein EIG-J7; Capsicum annuum, partial (84%)
6.0e-13 stress response
similar to 4Fe-4S ferredoxin iron-sulfur binding domain protein;
Alkalilimnicola ehrlichei, partial (7%)
2.8 e-34 transport
similar to chromosome chrl scaffold_5 whole genome shotgun
sequence; Wis vinifera, Rep:, partial (13%)
Similar to Zea mays 18S ribosomal RNA gene partial sequence,
3.1 e-80 translation
similar to chromosome chr15 scaffold_40 whole genome shotgun
sequence; Vitis vinifera, partial (10%)
3.1 e-27 unknown
Validation of cDNA-AFLP expression patterns
We conclude that the cDNA-AFLP experiments showed the effect of MHS on the meiotic anther transcriptome to be visible as of 2 h of treatment with different effects in heat sensitive and tolerant genotypes. In response to MHS, a differential modulation of approximately 1% of the anther transcriptome was observed in the MM genotype. In an effort to identify components involved in the heat stress response, TDFs showing differential gene expression during the heat treatment have been sequenced and analysed further. The identified genes are significantly changed in expression and indicate that meiotic anther development is sensitive to heat stress as early as 2 h of MHS.
Microarray analysis complements cDNA-AFLP analysis and indicates differential responses to heat stress in heat-tolerant and heat-sensitive genotypes
In tomato, cDNA-AFLP and microarray technology have been previously used to dissect stress responses and already enabled identification of candidate genes for tolerance to both abiotic and biotic stresses such as high temperature, low temperature, salt stress and pathogen attack [23, 29–31]. We used the Combimatrix Tomato Array 1.0 to compare and analyse the response to heat in a tolerant and in a sensitive genotype. Modulation of gene expression in the two genotypes was analysed using multiclass, paired, and un-paired Significance Analysis of Microarray; SAM .
To evaluate the relationship between genome-wide expression profiles of the two genotypes in relation to the heat treatment, SAM multiclass analysis was performed for three consecutive time points, 0 h-2 h-6 h (Additional file 1, SAM multiclass). Using this statistical approach we were able to identify 56 genes in HS1 and 75 genes in FR whose mean expression was significantly changed across 2 time points or over the entire duration of the heat treatment, at a false discovery rate (FDR) of 5%. Among these genes, 45 are shared between the two genotypes. Transient up-regulation at the 2 h time point in both genotypes was the most common profile, and comparable to that of cDNA-AFLP experiment. Interestingly, up-regulation showed a marked genotypic difference with 12 genes up-regulated in HS1, in comparison to only 5 in FR. Down-regulation of gene expression was the second most frequently observed modulation for FR (16 genes) but not for HS1 (9 genes). These modulation types also show a good correlation with the cDNA-AFLP data, where the sensitive MM genotype is also characterised by a considerable down-regulation of gene expression. These results suggest that the difference in tolerance to MHS between the two genotypes is mainly associated with a moderated transcriptomic response in the heat-tolerant genotype with respect to the number of differentially expressed genes and the level of expression changes.
A further analysis of the heat response per genotype was achieved by looking at differentially expressed genes between three pair wise comparisons with a SAM paired analysis: t0-t2, t2-t6, and t0-t6 (Additional file 1, SAM paired). The statistical nature of the SAM paired analysis, in which samples are assigned to one group and there is also a one-to-one pairing between biological replicates, reveals additional information on subtle differences in the response to heat between the two genotypes, particularly in the sensitive genotype. We detected 74 genes in HS1 and 137 in FR, which exhibited significant changes in expression (FDR = 5%). As in the multiclass analysis, the majority were significantly different for the 2 h time point (t0-t2) with 38 genes significantly changed in HS1 compared to 72 genes in FR. The SAM paired analysis also shows that, in the heat-tolerant HS1 genotype, less transcriptional changes are invested into reacting to the heat stress than in the heat-sensitive genotype FR. Overall, the paired analysis shows that HS1 up-regulates 56 and down-regulates 49 genes, while FR up-regulates 72 genes and down-regulates 109, clearly more when compared to the number of down-regulated genes in the heat-tolerant genotype. The SAM paired analysis allows a more in-depth analysis compared to the SAM multiclass analysis, suggesting not only that the difference in tolerance between the two genotypes is indeed associated with a lower transcriptional response of the heat-tolerant genotype but also with more functional classes and members affected in the sensitive genotype.
To detect if significant differences can also be found between the two genotypes prior to the heat treatment, we used the unpaired SAM analysis (Additional file 1, SAM unpaired). By using this method, where biological replicates are not paired each other between two samples, 24 genes were identified as significantly different between HS1 and FR (FDR = 5%). Of these, 4 are significantly higher in the heat-tolerant genotype, while the other 20 genes are at significantly higher levels in the heat-sensitive genotype. The 4 genes more highly expressed in the tolerant genotype and some of the other 20 more highly expressed in the sensitive genotype are further discussed below with regards to their functional classification. A previous study in salt cress proposed that the salt tolerant plant had a smaller number of salt-regulated genes because of the constitutive high level of expression of stress protection-related genes even under standard growth conditions . Similarly, Frank et al.  showed that constitutive expression of a heat shock transcription factor and several heat shock proteins mark these genes as candidates for taking part in microspore thermotolerance. In our case the genes with constitutively higher expression in the tolerant genotype are the heat shock Hsp82 (TC170030), a gene coding for mitochondrial small heat shock protein (msHsp) LeMtHSP (TC187014), a cathepsin B-like cysteine proteinase (TC171192) and the fructose-1,6-bisphosphate aldolase gene (TC176475). We suggest that in HS1, the differentially expressed genes prior to heat treatment may represent a constitutive tolerance to MHS.
The statistical analysis of gene expression during MHS reveals a differential response to heat in the two genotypes, in terms of number of significantly differentially expressed transcripts, their levels of expression and the functional classes they belong to. Moreover, the analysis also shows that significant differences in gene expression are present prior to the commencement of the MHS treatment.
Unlike the microarray analysis, the cDNA-AFLP identified genes involved in the response to heat, which may be interesting as candidates for MHS tolerance. Furthermore, of the 25 sequenced and identified TDFs, 19 were also present on the microarray and 3 even appeared as significantly changed upon MHS on the microarray (TC170123/TC211882; a hsp70 like fragment, TC169993/TC191544; the hsa32 also involved in heat shock response and TC179740/TC213983; a fragment of unknown function). The alternate TC numbering arose from the updated Gene Index Release version 8 to version 11 by the DFCI database (included here for clarity).
Main components determining gene expression modulation under MHS
Transcript profile clustering under MHS
Functional groups modulated by MHS
Protection and repair transcripts represent the majority of genes affected by MHS
From the SAM paired analysis, the highest induced transcripts during the first 2 h of MHS in the tolerant genotype belong to the protection and repair genes and are similar to: AthHsp22.3 (TC185802, fold induction 1500), a class I heat shock protein (TC184669, fold induction 170), a class III sHsp (TC173974, fold induction 106), the cytosolic ascorbate peroxidase (TC182989, fold induction 36.6) and a class I sHsp (TC178169, fold induction 35.8). The highest induced transcripts in the sensitive genotype during the first two hours of MHS are a class I heat shock protein (TC184669, fold induction 73), a DNA-J (176514, fold induction 55), a cytosolic ascorbate peroxidase (TC182929, fold induction 34) and mitochondrial sHsp (TC187014, fold induction 32). However, the highest induction over the entire duration of MHS was recorded in both genotypes and consists of a similar group in each genotype: class I sHsp 17.6 (TC178169), class II sHsp 17.6 (TC170069), and mitochondrial sHsp (TC187014). Of these, only the class I sHsp 17.6 is more highly induced in the tolerant genotype (fold induction 36 as compared to 26 in FR). A single transcript corresponding to a DNA-J like protein (TC176727) was down-regulated in FR over the 6 h of MHS, but during the last 4 h of MHS, 15 heat shock proteins were down-regulated in expression in FR and 18 in HS1.
An additional set of genes that are differentially expressed exclusively in the heat-sensitive genotype, include heat shock proteins Hsa32 (TC191544), Rof1 (TC187143) and ROC7 (TC175253) genes that are required for protein folding. It is evident that oxido-reduction processes are affected by MHS in the heat-sensitive genotype from the differential expression of a carotenoid cleavage dioxygenase 1 (TC181718), a steroid alpha-reductase, (TC170471), a flavonol synthase (TC172800) and a flavoprotein (TC83140). The stress response of the heat-sensitive genotype can be seen from repression of transcripts similar to pepsin A (TC171322) and subtilisine (TC181342). Unlike the heat-sensitive genotype, the heat-tolerant genotype HS1 shows fewer genes significantly induced in the protection and repair group (Figure 6).
In conclusion, the heat-tolerant genotype HS1 reacts to MHS by enhancing the expression of heat shock proteins, oxido-reductive, transport and stress protective transcripts. In stark contrast, the heat-sensitive FR shows a more extensive, wider transcriptomic response characterised first by a complementary modulation pattern as compared to the heat-tolerant genotype. For example, transport and signalling genes are down-regulated in the heat-sensitive and up-regulated in the heat-tolerant. Secondly, the stress response in the heat-sensitive genotype is characterised by an amplified expression of a more diverse range of functional classes including transcription, photosynthesis, hormone related and protection genes which are increased, while transport and carbohydrate metabolism are reduced in expression. Thirdly, although the two genotypes involve similar functional classes in the response to heat, each genotype modulates "common" and "specific" genes, as described above. The genotype-specific modulated transcripts are of particular interest for the present work as their annotation may provide clues to the putative mechanism of heat tolerance.
The results obtained in the present study are similar to those obtained by Frank et al. , where the heat-tolerant genotype is shown to modulate fewer genes in response to heat, but also similar to results on other abiotic stresses [36–38]. The most marked differences between the two genotypes can be seen in the heat shock, metabolism and oxido-reduction groups: 32% (HS1) respectively 20% (FR) of the modulated genes belong to the heat shock group, 12% (HS1) respectively 20% (FR) belong to metabolic genes group and 8% (HS1) respectively 10% (FR) belong to the oxido-reduction group. It appears that the heat-tolerant genotype modulates mainly genes in the heat shock group while the heat-sensitive genotype modulates genes of the heat shock and metabolic groups in equal fractions; processes which may underpin the better performance of the heat-tolerant genotype under elevated temperatures.
Metabolism and development transcripts are the second most affected by MHS
Interestingly, in comparison to the heat-tolerant genotype, all development-related genes significantly changed in expression by MHS were very highly expressed in the heat-sensitive genotype FR. In the metabolism group it is interesting to note that FR exhibits increased expression levels of carbohydrate metabolism genes as well, for example a transcript coding for beta-amylase (TC178310) is induced 3 fold in FR during the first 2 h of MHS but then it is repressed again. ß-amylase induction and the resultant maltose accumulation may function as a compatible-solute stabilizing factor in the chloroplast stroma in response to acute temperature stress .
Transcription and signalling transcripts are also affected by MHS
Unlike the heat-tolerant genotype HS1, the heat-sensitive FR shows a different response with regard to transcription, translation and signalling. The ethylene-responsive transcriptional co-activator (TC171720) is transiently induced at 2 h by heat only in the sensitive genotype, which may also reflect the sensitivity of FR to MHS. It is known that the A. thaliana ortholog, AtMBF1c, enhances tolerance to heat and osmotic stress when over-expressed in A. thaliana. Unlike the heat-tolerant genotype HS1, FR reacts with a different set of signalling partners and transcriptional regulators, down-regulating after 2 h of MHS a photoperiod responsive protein (TC177921) and an ARF5-like transcript (TC182088). In Arabidopsis hypocotyls, high temperature (29°C) causes an increase in free IAA, and this mediates a marked cell elongation response. Therefore, temperature signals may modify auxin synthesis or distribution in the plant, and this could represent a growth-regulating mechanism . The heat-tolerant genotype reacts to heat by continuously increasing the expression of an Adagio 3-like transcript (TC188199) while another Bel1 transcriptional regulator (TC175335) is decreased at 6 h of MHS.
In addition, the heat-tolerant genotype not only reacted differently from the heat-sensitive during heat stress but also exhibited a constitutive and specific gene expression pattern, characterised by very high expression levels of protection and repair genes. The heat shock Hsp82 (TC170030) and the gene coding for mitochondrial small heat shock protein (msHsp) LeMtHSP (TC187014) are highly expressed in HS1 at time point 0 h prior to the MHS. Hsp82 is an hsp90 chaperone required for cell growth, is highly conserved among eukaryotes and the presence of at least one of the HSP90 gene product family members is essential for viability in yeast, Drosophila, and humans. In yeast, Hsp82 is part of the Hsp90 chaperone, which directly interacts with Hsf1p, the heat-shock transcription factor from yeast . The reduction in tomato mitochondrial small hsp LEMTSHP expression, also known as Hsp23.8  seems to lead directly to susceptibility to heat stress. In addition, a transcript similar to a cathepsin B-like cysteine proteinase (TC171192) is also constitutively higher expressed in the tolerant genotype. Cysteine proteinases have a positive role in plant growth, development, senescence and programmed cell death, but also in storage protein mobilization . All these 3 genes are significantly increased after 2 h of MHS. The heat-tolerant genotype also exhibited a constitutive and specific gene expression pattern, characterised by differential expression levels of a carbohydrate metabolism gene, the fructose-1,6-bisphosphate aldolase (TC176475), which is highly expressed in HS1 at time point 0 h prior to the MHS. The fructose-1,6-bisphosphate aldolase is an enzyme of the glycolytic pathway. Reports from several laboratories have suggested that increased rates of glycolysis play an essential role in the initiation of DNA synthesis and may be involved in maintenance of replication and protein activity at high temperature .
It is interesting to note that in general, the heat-tolerant genotype HS1 appears to have an enhanced innate heat protection system exemplified by the hsp82 and hsp90 levels before and during the treatment. Furthermore, HS1 seems to have a lower rate of metabolism as exemplified by the low number of genes significantly changed in comparison to the sensitive genotype, not only before the treatment but also during the 6 h of MHS.
Unlike HS1, the sensitive genotype FR exhibits an increased number of genes (20) with a constitutively higher expression prior to the MHS, which may also indicate the sensitivity of this genotype to adverse conditions. For example, the highest constitutive expression in the sensitive genotype compared to the tolerant genotype is of genes that belong to the stress response pathway (TC189289; pepsin A, TC181342; subtilisin like protease), metabolism (TC176104; a transferase) and interestingly, in development (TC182608; histone H3, TC171121; male sterility 2 gene) or transcription (TC182088; ARF5). Of the genes constitutively higher expressed in the sensitive genotype, 19 are down-regulated within the first 2 h of MHS and only 1 is not affected at all by the MHS treatment (TC172148; an aquaporin involved in water and neutral solute transport).
The analysis of the stress response reveals that the tolerant and sensitive genotypes react largely by involving similar functional classes: heat shock, oxido-reduction, photosynthesis, stress response, protection, development, transcription and transport. The nature of the genes and intensity of their expression level is very genotype specific as has been described above. However, the group of transcript significantly changed in both genotypes shows probably the general response to heat and comprises in majority genes from the heat shock, metabolism and oxido-reduction functional classes.
The majority of "common" transcripts generally show the same induction pattern in both genotypes, but very often at different expression levels. However, the question whether they are relevant for tolerance is still open. For example, expression of a transcript similar to hsp 17.6 (TC170069), a class I shsp is higher and transiently induced at two hours of MHS in FR. A transcript similar to Hsc70 (TC176566) is also more highly expressed in FR as an result of heat. Other transcripts similar to Class I and Class III heat shock proteins are expressed more highly in the tolerant genotype in the first 2 h of stress. Oxidation and reduction genes are also induced over the first 2 h of heat and they are similarly induced in both genotypes, for example a cytosolic ascorbate peroxidase (TC170369); a central component of the reactive oxygen gene network in Arabidopsis , is more highly induced in the tolerant genotype than in the sensitive genotype. A metallothionein-like protein type 2 B (TC179817) involved in cellular response to stress  is expressed and transiently induced at a higher expression level in the sensitive genotype. TC170658, corresponding to a chalcone synthase 1 B involved in flavonoids biosynthesis, well known for its antioxidant properties , continuously decreases in expression in the two genotypes but the decrease is more pronounced in the sensitive genotype during the first 2 h of MHS. Photosynthesis is also affected, for example RUBISCO activase (TC176226) is equally increased at 2 h of MHS then decreased at the end of the heat treatment in both genotypes, however with a higher level in the sensitive genotype. Interestingly, expression of a transcript similar to the Zeatin O-glucosyltransferase ZOG (TC187640) decreases in the tolerant genotype over the entire duration of MHS but transiently decreasing in the sensitive genotype. The significance of this particular expression pattern may lie in the fact that ZOG regulates responses to water deficits and thus regulates stress protection by increasing active cytokinin levels .
It is clear that numerous biochemical pathways and metabolic routes are affected by MHS and the range of responses is generally broader in the sensitive genotype than in the tolerant. Thus it is not yet clear which genes are most directly involved in tolerance or in sensitivity to heat stress.
Validation of microarray expression data
Our transcript profiling experiments and those from others  demonstrate that many genes, involved in seemingly unrelated processes, are modulated by moderate heat stress. Expression of many of these genes is also modified in sensitive as well as in relative tolerant lines, differing only in amplitude. This suggests that tolerance is based on fine-tuning of quantitative expression of many genes. It is difficult to say whether the differences observed in our experiments are the basis of heat tolerance or solely a consequence of a better performing plant under stress conditions independent of a causative link. Furthermore, it is possible that different mechanisms exist in tomatoes, which have been selected during cultivation in different part of the world. The selection of tomato in cultivation is based on a small set of domestication genes causing little genetic variation in the tomato germplasm . Based on this, an unintended selection during cultivation for epigenetic factors causing differential gene expression could partly explain differences in stress adaptation.
Plant material and heat stress conditions
Anther isolation and RNA extraction
We focused our analysis of gene expression on whole meiotic anther cones isolated from flower buds of plants that were exposed to MHS (32°C/26°C, day/night) for up to 30 h, and compared this to standard temperatures (26°C/18°C day/night). For the cDNA-AFLP and for the microarray experiment and in an attempt to reduce the biological variation, a pool of six meiotic flower buds from 3 plants were harvested (according to ). Pilot experiments indicated the third and fourth flower clusters to give the most consistent phenotypic response to heat stress (data not shown). From these, anther cones ranging in size from 2-4 mm were isolated, immediately frozen in liquid N2 and pooled for RNA extraction. Thus, tissues from pooled samples from 3 biological repeats were homogenized using glass rods and RNA was extracted using a Plant RNA isolation kit (Plant RNeasy Mini Kit, Qiagen Benelux BV, Venlo the Netherlands) according to the manufacturer's recommendations.
Transcript profiling by cDNA-AFLP
For cDNA-AFLP, material was harvested at 0 h (ST) and at 2 h, 6 h, 16 h and 30 h from both ST and MHS treatments. The isolated RNA was subjected to cDNA-AFLP template preparation as described previously . RNA fingerprinting was carried out using 92 primer combinations with 2 selective nucleotides and gave rise to an average of 80 transcript-derived fragments per primer combination (TDFs). Primer sequences were as described in Bachem et al.  (Additional file 2). The majority of bands showed no change in intensity in response to heat. Changes in the intensity of individual bands did not affect others in the same lane, indicating that product accumulation was not affected by the concentration of individual substrates in the reaction. Inconsistent bands were mainly observed in the region of the gel with very small DNA fragments (< 70 bp) and were therefore considered as aspecific amplification products. Bands corresponding to differentially expressed genes were cut out from the gel and the eluted DNA was re-amplified under the same conditions as for the selective amplification. Fragments were subsequently ligated in a T-tailed EcoRV digested phagemid (pBlueScriptII SK(+), Stratagene, La Jolla, CA, USA) and sequenced (CEQ™ DTCS Quick Start Kit and CEQ2000 DNA Analysis System, Beckman Coulter, Fullerton, CA, USA). Fragments that did not have the expected size, based on the height on the acryl amide gel from which they were isolated, were discarded. All cDNA-AFLP expression patterns displayed in the Results section, was confirmed by real-time quantitative PCR and using the microarray data.
Microarray hybridization and data analysis
The 90 K Custom TomatoArray 1.0 chip (Combimatrix microarray platform, http://ddlab.sci.univr.it/FunctionalGenomics/) consists of 20200 unique probes derived from Solanum lycopersicum transcripts and various controls, and was produced by the Plant Functional Genomics Center, University of Verona. The gene specific probes (oligonucleotide of 35-40-mer) randomly distributed in quadruplicate across the array, were designed using the program design OligoArray 2.1 . The sequences represented on the chip correspond to 20115 S. lycopersicum TCs (Gene Index Release 11.0, 21 June, 2006) and to 85 technical controls (negative and positive spiking controls). From the 20115 tomato probes represented on the chip, 17018 probes were found to be expressed above background and consistent among the technical repeats within each chip. The quality of the biological replicates was evaluated by Pearson coefficient, which ranged from 0.87 to 0.98. For the microarray experiment, 3 biological replicas for each sample were used. Material was harvested at 0 h (ST) and at 2 h and 6 h from both genotypes (FR and HS1). Amino allyl-RNA synthesis (aRNA) and labelling with the Alexa647 dye were performed from 1 μg of total RNA, with the "SuperScript TM Indirect RNA Amplification System" (Invitrogen, Carlsbad, CA) according to manufacturer's recommendations. Labelled aRNA was quantified by spectrophotometer and the efficiency of Cy5 dye incorporation (DOL) was calculated. 4 μg of labelled RNA, with DOL value ranging from 2.0 to 3.0, were first fragmented and hybridized to the array as indicated by the manufacturer (http://www.combimatrix.com) Pre-hybridization, hybridization, washings and imaging were performed according to the protocols given by CombiMatrix. The array was scanned with a ScanArray 4000 × L microarray scanner (Perkin Elmer, Waltham, MA USA). Tiff images were exported to the Microarray Imager 5.8 (Combimatrix) for the densitometry analysis of the spots.
The microarray data was normalized by median scaling and the quality of the biological replicates was assessed by means of Pearson's coefficient. Only those genes that were up or down-regulated at least twofold, relative to the average of the three negative controls, for at least one time point, were included in the analysis. The two-fold change in expression as threshold for consideration is a convention employed in various transcript profiling studies  and therefore is used here to assess the findings from the presented work in comparison with those of earlier reports. All microarray expression data are available at GEO under the series entry GSE24805. The resulting set of significantly changed expression values has been taken into further analysis with paired, unpaired and multiclass SAM analysis, to detect the effect of 'treatment' on expression levels. Genes with similar expression patterns were grouped according to a hierarchical clustering algorithm using Euclidean Distance coefficients.
The paired and unpaired SAM statistical analyses were carried out to determine the significant differences in gene expression between two time-points for each genotype and the significant differences in gene expression between individual time-points observed in both genotypes, respectively. These analyses were performed by T-Mev with a FDR = 5%. For each analysis, a new input dataset was created by filtering the data of the two samples that were compared on the basis of the C.V. values.
In order to determine the significant differences in gene expression accumulated over the entire duration of the experiment all data were filtered on the basis of the CV value (CV < 0.5) among spot replicates present on the chip and the restricted dataset obtained, of 17018 genes, was suited to run a multiclass comparison method of Significance Analysis of Microarray (SAM), with false discovery rate (FDR) = 5% .
The main components of the changes in gene expression in response to MHS were determined using the set of differentially expressed genes obtained by SAM multiclass analysis, by performing an O2PLS analysis with SIMCA P+ (Umetrics, Umea, Sweden). OPLS combines the existing theory of partial least squares (PLS) regression  and orthogonal signal correction (OSC) . The unique property of O2PLS is its capacity to identify joint variation between two datasets, while acknowledging also systematic variation that is unique to each dataset. The expression profile clustering for each genotype was carried out with the Hierarchical Clustering (HCL) method with Euclidean Distance using T-Mev software (Version 4.3).
Validation of gene expression by q-PCR
Total DNA-free RNA was isolated from tomato anthers (n > 5) using a RNA isolation kit (Qiagen RNeasy mini kit with the RNase-Free DNase Set, Qiagen). A PCR reaction with RNA- and DNA-specific primers (Additional file 2) based on an intron in the S. lycopersicum actin gene Tom51 (GenBank accession number U60481) was performed to ensure the absence of contaminating genomic DNA. Total DNA-free RNA (1 μg) was used for cDNA synthesis (iScript™ cDNA Synthesis Kit; Bio-Rad Laboratories, Hercules, CA, USA) in a total volume of 25 μl. PCR reactions were carried out in 25 μl containing 0.125 μl of cDNA synthesis reaction mixture, 400 nm of each primer and 12.5 μl of iQ SYBR Green Supermix (Bio-Rad Laboratories). PCRs were performed in a 96-well Bio-Rad iCycler (Bio-Rad Laboratories) using a temperature program starting with 3 min at 95°C followed by 40 cycles consisting of 15 s at 95°C and 45 s at 57°C, and finally the melting temperature of the amplified product was determined to verify the presence of a specific product. Fold changes in expression levels were calculated using the MS-Excel macros (Biorad) where the lowest expression value was set to one.
In addition, a fraction of the PCR mixture was analysed on a 1% agarose/ethidium bromide gel to check the size of the amplified DNA fragment. The primers that were used for the real-time quantitative PCR reactions were designed using a computer program (Beacon Designer 5.01; Premier Biosoft International, Palo Alto, CA, USA) to obtain primers that have close to identical melting temperatures and do not form secondary structures with each other in the given PCR conditions. Primer sequences are listed in Additional file 2. In addition, to enhance primer efficiency, primer-binding sites were chosen such that secondary structures of the template were avoided. All reactions were performed on a pool of 6 anther cones collected from 3 different plants. The pooling has been done in such a way as to minimize the eventual biological variation.
Functional classification criteria
For the microarray analysis, the 20115 tomato probes were designed from sequences deposited in the tomato TGI database . Thus, the annotations of the gene sets on the cDNA-AFLP or on the microarray were taken from the DFCI Tomato Gene Index (LGI), which integrates research data from all international tomato gene research projects. Where biological ontology was not clear, the UniProt Knowledgebase provided further functional information.
All microarray expression data are available at GEO under the series entry GSE24805.
List of abbreviations
cDNA Amplified Fragment Length Polymorphism.
Acknowledgements and Funding
The authors thank Nunhems BV for making germplasm (HS1 and FR) available. The authors would like to mention that the significance analysis of microarray (SAM) has been exclusively performed by the Italian collaborator at the Dipartimento di Biotecnologie, Universitá degli Studi di Verona, Italy.
This work was supported by the EU-EST Marie Curie [ADONIS, 2000327 to C.E.B] and additionally by a grant [personal, to C.E.B] from the Institute of Water and Wetland Research Nijmegen.
- Drews GN, Goldberg RB: Genetic control of flower development. Trends in Genetics. 1989, 5 (8): 256-261.PubMedView ArticleGoogle Scholar
- Goldberg RB: Plants: novel developmental processes. Science. 1988, 240 (4858): 1460-1467. 10.1126/science.3287622.PubMedView ArticleGoogle Scholar
- Aarts MG, Hodge R, Kalantidis K, Florack D, Wilson ZA, Mulligan BJ, Stiekema WJ, Scott R, Pereira A: The Arabidopsis MALE STERILITY 2 protein shares similarity with reductases in elongation/condensation complexes. The Plant Journal. 1997, 12 (3): 615-623. 10.1046/j.1365-313X.1997.d01-8.x.PubMedView ArticleGoogle Scholar
- Ariizumi T, Hatakeyama K, Hinata K, Sato S, Kato T, Tabata S, Toriyama K: A novel male-sterile mutant of Arabidopsis thaliana, faceless pollen-1, produces pollen with a smooth surface and an acetolysis-sensitive exine. Plant Molecular Biology. 2003, 53 (1-2): 107-116.PubMedView ArticleGoogle Scholar
- Chen YC, McCormick S: sidecar pollen, an Arabidopsis thaliana male gametophytic mutant with aberrant cell divisions during pollen development. Development. 1996, 122 (10): 3243-3253.PubMedGoogle Scholar
- Hulskamp M, Parekh NS, Grini P, Schneitz K, Zimmermann I, Lolle SJ, Pruitt RE: The STUD gene is required for male-specific cytokinesis after telophase II of meiosis in Arabidopsis thaliana. Dev Biol. 1997, 187 (1): 114-124. 10.1006/dbio.1997.8554.PubMedView ArticleGoogle Scholar
- Taylor PE, Glover JA, Lavithis M, Craig S, Singh MB, Knox RB, Dennis ES, Chaudhury AM: Genetic control of male fertility in Arabidopsis thaliana: structural analyses of postmeiotic developmental mutants. Planta. 1998, 205 (4): 492-505. 10.1007/s004250050348.PubMedView ArticleGoogle Scholar
- van der Meer IM, Stam ME, van Tunen AJ, Mol JN, Stuitje AR: Antisense inhibition of flavonoid biosynthesis in petunia anthers results in male sterility. The Plant cell. 1992, 4 (3): 253-262.PubMed CentralPubMedView ArticleGoogle Scholar
- Kaul MLH: Male sterility in higher plants. 1988, Berline: Springer-Verlag, 10:View ArticleGoogle Scholar
- Cnudde F, Hedatale V, de Jong H, Pierson ES, Rainey DY, Zabeau M, Weterings K, Gerats T, Peters JL: Changes in gene expression during male meiosis in Petunia hybrida. Chromosome Res. 2006, 14 (8): 919-932. 10.1007/s10577-006-1099-5.PubMedView ArticleGoogle Scholar
- Abiko M, Akibayashi K, Sakata T, Kimura M, Kihara M, Itoh K, Asamizu E, Sato S, Takahashi H, Higashitani A: High-temperature induction of male sterility during barley (Hordeum vulgare L.) anther development is mediated by transcriptional inhibition. Sexual Plant Reproduction. 2005, 18 (2): 91-100. 10.1007/s00497-005-0004-2.View ArticleGoogle Scholar
- Lindquist S, Craig EA: The heat-shock proteins. Annu Rev Genet. 1988, 22: 631-677. 10.1146/annurev.ge.22.120188.003215.PubMedView ArticleGoogle Scholar
- Porch TG, Jahn M: Effects of high-temperature stress on microsporogenesis in heat-sensitive and heat-tolerant genotypes of Phaseolus vulgaris. Plant Cell and Environment. 2001, 24 (7): 723-731. 10.1046/j.1365-3040.2001.00716.x.View ArticleGoogle Scholar
- Pressman E, Peet MM, Pharr DM: The effect of heat stress on tomato pollen characteristics is associated with changes in carbohydrate concentration in the developing anthers. Ann Bot. 2002, 90 (5): 631-636. 10.1093/aob/mcf240.PubMed CentralPubMedView ArticleGoogle Scholar
- Sato S, Peet MM, Gardner RG: Altered flower retention and development patterns in 9 tomato cultivars under elevated temperature. Hort Sci. 2004, 101: 95-101. 10.1016/j.scienta.2003.10.008.View ArticleGoogle Scholar
- Sato S, Peet MM: Effects of moderately elivated temperature stress on the timing of pollen release and its germination in Tomato (Lycopersicon exculentum Mill.). J Hort Sci Biotech. 2005, 80 (1): 23-28.Google Scholar
- Kamel A, Qian Y, Kolbe E, Stafford C: Development and validation of a multiresidue method for the determination of neonicotinoid and macrocyclic lactone pesticide residues in milk, fruits, and vegetables by ultra-performance liquid chromatography/MS/MS. J AOAC Int. 2010, 93 (2): 389-399.PubMedGoogle Scholar
- Sato Y, Murakami T, Funatsuki H, Matsuba S, Saruyama H, Tanida M: Heat shock-mediated APX gene expression and protection against chilling injury in rice seedlings. J Exp Bot. 2000, 52 (354): 145-151.View ArticleGoogle Scholar
- Levy A, Rabinowitch HD, Kedar N: Morphological and physiological characters affecting flower drop and fruit set of tomatoes at high temperatures. Euphytica. 1978, 27 (1): 211-218. 10.1007/BF00039137.View ArticleGoogle Scholar
- Sato Y, Aoki Y, Matsuoka M: A loss-of-function mutation in the rice KNOX type homeobox gene, OSH3. Plant Cell Physiol. 2002, 43 (1): 44-51. 10.1093/pcp/pcf003.PubMedView ArticleGoogle Scholar
- Kim JC, Lee SH, Cheong YH, Yoo CM, Lee SI, Chun HJ, Yun DJ, Hong JC, Lee SY, Lim CO, Cho MJ: A novel cold-inducible zinc finger protein from soybean, SCOF-1, enhances cold tolerance in transgenic plants. Plant J. 2001, 25 (3): 247-259. 10.1046/j.1365-313x.2001.00947.x.PubMedView ArticleGoogle Scholar
- Kakani VG, Reddy KR, Koti S, Wallace TP, Prasad PV, Reddy VR, Zhao D: Differences in in vitro pollen germination and pollen tube growth of cotton cultivars in response to high temperature. Ann Bot. 2005, 96 (1): 59-67. 10.1093/aob/mci149.PubMed CentralPubMedView ArticleGoogle Scholar
- Frank G, Pressman E, Ophir R, Althan L, Shaked R, Freedman M, Shen S, Firon N: Transcriptional profiling of maturing tomato (Solanum lycopersicum L.) microspores reveals the involvement of heat shock proteins, ROS scavengers, hormones, and sugars in the heat stress response. J Exp Bot. 2009, 60 (13): 3891-3908. 10.1093/jxb/erp234.PubMed CentralPubMedView ArticleGoogle Scholar
- Sato S, Peet MM, Thomas JF: Determining critical pre- and post-anthesis periods and physiological processes in Lycopersicon esculentum Mill. exposed to moderately elevated temperatures. J Exp Bot. 2002, 53 (371): 1187-1195. 10.1093/jexbot/53.371.1187.PubMedView ArticleGoogle Scholar
- Bachem CWB, Vriezen WH, Bita CE, Fernandez del Carmen A: AFLP-Based RNA Fingerprinting: Novel Variants and Applications. The Handbook of Plant Funcitonal Genomics. Edited by: Kahl G, Meksem K. 2008, Weinheim: Wiley-Blackwell, 21-37.View ArticleGoogle Scholar
- Solomon S, Plattner GK, Knutti R, Friedlingstein P: Irreversible climate change due to carbon dioxide emissions. Proc Natl Acad Sci USA. 2009, 106 (6): 1704-1709. 10.1073/pnas.0812721106.PubMed CentralPubMedView ArticleGoogle Scholar
- Mittler R, Kim Y, Song L, Coutu J, Coutu A, Ciftci-Yilmaz S, Lee H, Stevenson B, Zhu JK: Gain- and loss-of-function mutations in Zat10 enhance the tolerance of plants to abiotic stress. FEBS Lett. 2006, 580 (28-29): 6537-6542. 10.1016/j.febslet.2006.11.002.PubMed CentralPubMedView ArticleGoogle Scholar
- Sun W, Xu X, Zhu H, Liu A, Liu L, Li J, Hua X: Comparative transcriptomic profiling of a salt-tolerant wild tomato species and a salt-sensitive tomato cultivar. Plant Cell Physiol. 2010, 51 (6): 997-1006. 10.1093/pcp/pcq056.PubMedView ArticleGoogle Scholar
- Gabriels SH, Takken FL, Vossen JH, de Jong CF, Liu Q, Turk SC, Wachowski LK, Peters J, Witsenboer HM, de Wit PJ: CDNA-AFLP combined with functional analysis reveals novel genes involved in the hypersensitive response. Mol Plant Microbe Interact. 2006, 19 (6): 567-576. 10.1094/MPMI-19-0567.PubMedView ArticleGoogle Scholar
- Zhang S, Scheller HV: Photoinhibition of photosystem I at chilling temperature and subsequent recovery in Arabidopsis thaliana. Plant Cell Physiol. 2004, 45 (11): 1595-1602. 10.1093/pcp/pch180.PubMedView ArticleGoogle Scholar
- Zhu JK: Plant salt tolerance. Trends Plant Sci. 2001, 6 (2): 66-71. 10.1016/S1360-1385(00)01838-0.PubMedView ArticleGoogle Scholar
- Tusher VG, Tibshirani R, Chu G: Significance analysis of microarrays applied to the ionizing radiation response. Proc Natl Acad Sci USA. 2001, 98 (9): 5116-5121. 10.1073/pnas.091062498.PubMed CentralPubMedView ArticleGoogle Scholar
- Taji T, Seki M, Satou M, Sakurai T, Kobayashi M, Ishiyama K, Narusaka Y, Narusaka M, Zhu JK, Shinozaki K: Comparative genomics in salt tolerance between Arabidopsis and aRabidopsis-related halophyte salt cress using Arabidopsis microarray. Plant Physiol. 2004, 135 (3): 1697-1709. 10.1104/pp.104.039909.PubMed CentralPubMedView ArticleGoogle Scholar
- Bylesjo M, Eriksson D, Kusano M, Moritz T, Trygg J: Data integration in plant biology: the O2PLS method for combined modeling of transcript and metabolite data. The Plant journal: for cell and molecular biology. 2007, 52 (6): 1181-1191.View ArticleGoogle Scholar
- Qin D, Wu H, Peng H, Yao Y, Ni Z, Li Z, Zhou C, Sun Q: Heat stress-responsive transcriptome analysis in heat susceptible and tolerant wheat (Triticum aestivum L.) by using Wheat Genome Array. BMC Genomics. 2008, 9: 432-10.1186/1471-2164-9-432.PubMed CentralPubMedView ArticleGoogle Scholar
- Page D, Gouble B, Valot B, Bouchet JP, Callot C, Kretzschmar A, Causse M, Renard CM, Faurobert M: Protective proteins are differentially expressed in tomato genotypes differing for their tolerance to low-temperature storage. Planta. 2010, 232 (2): 483-500. 10.1007/s00425-010-1184-z.PubMedView ArticleGoogle Scholar
- Swindell WR, Huebner M, Weber AP: Transcriptional profiling of Arabidopsis heat shock proteins and transcription factors reveals extensive overlap between heat and non-heat stress response pathways. BMC Genomics. 2007, 8: 125-10.1186/1471-2164-8-125.PubMed CentralPubMedView ArticleGoogle Scholar
- Zhou N, Robinson SJ, Huebert T, Bate NJ, Parkin IA: Comparative genome organization reveals a single copy of CBF in the freezing tolerant crucifer Thlaspi arvense. Plant Mol Biol. 2007, 65 (5): 693-705. 10.1007/s11103-007-9235-y.PubMedView ArticleGoogle Scholar
- Kaplan F, Guy CL: beta-Amylase induction and the protective role of maltose during temperature shock. Plant Physiol. 2004, 135 (3): 1674-1684. 10.1104/pp.104.040808.PubMed CentralPubMedView ArticleGoogle Scholar
- Suzuki N, Bajad S, Shuman J, Shulaev V, Mittler R: The transcriptional co-activator MBF1c is a key regulator of thermotolerance in Arabidopsis thaliana. J Biol Chem. 2008, 283 (14): 9269-9275. 10.1074/jbc.M709187200.PubMedView ArticleGoogle Scholar
- Gray WM, Ostin A, Sandberg G, Romano CP, Estelle M: High temperature promotes auxin-mediated hypocotyl elongation in Arabidopsis. Proc Natl Acad Sci USA. 1998, 95 (12): 7197-7202. 10.1073/pnas.95.12.7197.PubMed CentralPubMedView ArticleGoogle Scholar
- Borkovich KA, Farrelly FW, Finkelstein DB, Taulien J, Lindquist S: Hsp82 Is an Essential Protein That Is Required in Higher Concentrations for Growth of Cells at Higher Temperatures. Mol Cell Biol. 1989, 9 (9): 3919-3930.PubMed CentralPubMedView ArticleGoogle Scholar
- Liu JA, Shono M: Characterization of mitochondria-located small heat shock protein from tomato (Lycopersicon esculentum). Plant and Cell Physiology. 1999, 40 (12): 1297-1304.PubMedView ArticleGoogle Scholar
- Grudkowska M, Zagdanska B: Multifunctional role of plant cysteine proteinases. Acta Biochim Pol. 2004, 51 (3): 609-624.PubMedGoogle Scholar
- Janniere L, Canceill D, Suski C, Kanga S, Dalmais B, Lestini R, Monnier AF, Chapuis J, Bolotin A, Marina T, Le Chatelier E, Ehrlich SD: Genetic Evidence for a Link Between Glycolysis and DNA Replication. Plos One. 2007, 2 (5):
- Davletova S, Rizhsky L, Liang H, Shengqiang Z, Oliver DJ, Coutu J, Shulaev V, Schlauch K, Mittler R: Cytosolic ascorbate peroxidase 1 is a central component of the reactive oxygen gene network of Arabidopsis. The Plant cell. 2005, 17 (1): 268-281. 10.1105/tpc.104.026971.PubMed CentralPubMedView ArticleGoogle Scholar
- Yeh SC, Hsieh HM, Huang PC: Transcripts of metallothionein genes in Arabidopsis thaliana. DNA Seq. 1995, 5 (3): 141-144.PubMedGoogle Scholar
- Pietta PG: Flavonoids as antioxidants. J Nat Prod. 2000, 63 (7): 1035-1042. 10.1021/np9904509.PubMedView ArticleGoogle Scholar
- Polanska L, Vicankova A, Novakova M, Malbeck J, Dobrev PI, Brzobohaty B, Vankova R, Machackova I: Altered cytokinin metabolism affects cytokinin, auxin, and abscisic acid contents in leaves and chloroplasts, and chloroplast ultrastructure in transgenic tobacco. Journal of experimental botany. 2007, 58 (3): 637-649.PubMedView ArticleGoogle Scholar
- Paran I, van der Knaap E: Genetic and molecular regulation of fruit and plant domestication traits in tomato and pepper. J Exp Bot. 2007, 58 (14): 3841-3852. 10.1093/jxb/erm257.PubMedView ArticleGoogle Scholar
- Giorno F, Wolters-Arts M, Grillo S, Scharf KD, Vriezen WH, Mariani C: Developmental and heat stress-regulated expression of HsfA2 and small heat shock proteins in tomato anthers. J Exp Bot. 2010, 61 (2): 453-462. 10.1093/jxb/erp316.PubMed CentralPubMedView ArticleGoogle Scholar
- Brukhin V, Hernould M, Gonzalez N, Chevalier C, Mouras A: Flower development schedule in tomato Lycopersicon esculentum cv. sweet cherry. Sexual Plant Reproduction. 2003, 15 (6): 311-320.Google Scholar
- Levy C, Donnadieu M, Halpern GM: Demonstration of specific serum IgG with an "in vitro" radioimmunoassay (RAST IgG) (author's transl). Pathol Biol (Paris). 1977, 25 (7): 505-508.Google Scholar
- Wold S, Sjostrom M, Eriksson L: PLS-regression: a basic tool of chemometrics. Chemometrics and Intelligent Laboratory Systems. 2001, 58 (2): 109-130. 10.1016/S0169-7439(01)00155-1.View ArticleGoogle Scholar
- Wold S, Antti H, Lindgren F, Ohman J: Orthogonal signal correction of near-infrared spectra. Chemometrics and Intelligent Laboratory Systems. 1998, 44 (1-2): 175-185. 10.1016/S0169-7439(98)00109-9.View ArticleGoogle Scholar
- van der Hoeven R, Ronning C, Giovannoni J, Martin G, Tanksley S: Deductions about the number, organization, and evolution of genes in the tomato genome based on analysis of a large expressed sequence tag collection and selective genomic sequencing. The Plant cell. 2002, 14 (7): 1441-1456. 10.1105/tpc.010478.PubMed CentralPubMedView ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.