Transcriptomic identification of candidate genes involved in sunflower responses to chilling and salt stresses based on cDNA microarray analysis
© Fernandez et al; licensee BioMed Central Ltd. 2008
Received: 30 May 2007
Accepted: 26 January 2008
Published: 26 January 2008
Considering that sunflower production is expanding to arid regions, tolerance to abiotic stresses as drought, low temperatures and salinity arises as one of the main constrains nowadays. Differential organ-specific sunflower ESTs (expressed sequence tags) were previously generated by a subtractive hybridization method that included a considerable number of putative abiotic stress associated sequences. The objective of this work is to analyze concerted gene expression profiles of organ-specific ESTs by fluorescence microarray assay, in response to high sodium chloride concentration and chilling treatments with the aim to identify and follow up candidate genes for early responses to abiotic stress in sunflower.
Abiotic-related expressed genes were the target of this characterization through a gene expression analysis using an organ-specific cDNA fluorescence microarray approach in response to high salinity and low temperatures. The experiment included three independent replicates from leaf samples. We analyzed 317 unigenes previously isolated from differential organ-specific cDNA libraries from leaf, stem and flower at R1 and R4 developmental stage. A statistical analysis based on mean comparison by ANOVA and ordination by Principal Component Analysis allowed the detection of 80 candidate genes for either salinity and/or chilling stresses. Out of them, 50 genes were up or down regulated under both stresses, supporting common regulatory mechanisms and general responses to chilling and salinity. Interestingly 15 and 12 sequences were up regulated or down regulated specifically in one stress but not in the other, respectively. These genes are potentially involved in different regulatory mechanisms including transcription/translation/protein degradation/protein folding/ROS production or ROS-scavenging. Differential gene expression patterns were confirmed by qRT-PCR for 12.5% of the microarray candidate sequences.
Eighty genes isolated from organ-specific cDNA libraries were identified as candidate genes for sunflower early response to low temperatures and salinity. Microarray profiling of chilling and NaCl-treated sunflower leaves revealed dynamic changes in transcript abundance, including transcription factors, defense/stress related proteins, and effectors of homeostasis, all of which highlight the complexity of both stress responses. This study not only allowed the identification of common transcriptional changes to both stress conditions but also lead to the detection of stress-specific genes not previously reported in sunflower. This is the first organ-specific cDNA fluorescence microarray study addressing a simultaneous evaluation of concerted transcriptional changes in response to chilling and salinity stress in cultivated sunflower.
Sunflower (Helianthus annuus L.) is the third most important source of edible vegetable oil worldwide which is also thought to become an efficient source of biodiesel (Sunflower Statistics NSA 2007, USA) . Considering that sunflower production is expanding to arid regions in the Mediterranean area, North America, India and Argentina, tolerance to drought and salinity arises as important issues for breeding programs [2–4]. In addition, requirements of early sow to maximize the growing season and to escape to drought stress have increased the need of better chilling tolerance, particularly at early stages of development. Molecular mechanisms involved in response to these stresses have been extensively studied in model species like Arabidopsis thaliana [5–7] and in important crop species like rice . The expression of a number of plant genes is regulated by abiotic environmental stresses including drought, high salinity and cold [9–12]. Transcriptome analysis using microarrays have proven to be a powerful tool for discovery of many stressed-induced genes involved in stress response and tolerance. Macro and microarray studies of abiotic stress responses in Arabidopsis and Oryza sativa allowed the identification of genes involving both functional and regulatory proteins [6, 8, 13–23]. The first group comprises membrane transporters and water channel proteins, key enzymes for osmolite biosynthesis; detoxification enzymes and macromolecules protection proteins. The second group comprises transcription factors (TFs) (i.e. bZIP, MYC, MYB, CREB/CBF, HD-ZIP), protein kinases and proteinases involved in the regulation of signal transduction and gene expression. These regulatory systems have been reported either as dependent or independent on abscisic acid (ABA) which indicate the existence of complex regulatory mechanisms between perception of abiotic stress signals and gene expression [20, 21].
Cross talk signaling cascades among drought, cold and salinity has been reported for A. thaliana and large number of stress-inducible genes were isolated and characterized including osmotic response genes as rd22BP1, AtMYBB2, DREB1A and DREB2A, signaling molecules that activates effectors as SOS3 (Ca2+ binding protein), SOS2 (Ca2+ dependent kinase), SOS1 (Na+/H+ membrane antiporter) [8, 22]. Regarding responses to low temperatures, cold-induced genes were reported in many species as lucerne , Arabidopsis [24–26], barley  and wheat [28, 29]. Many of these genes encode for proteins of unknown function, being some of them described as LEA proteins (Late Embryogenesis Abundant) . Other cold-resistance genes (COR) as LT1, KIN, RD and ERD have been isolated mainly from Arabidopsis. These genes present the CRT/DRE (C-repeat/dehydratation-responsive elements) in the promoter region that bind CBF and DREB (C-repeat binding factors/dehydration-responsive elements binding proteins) TFs. While cis elements in cold response genes bind DREB1/CBF TFs, regulatory regions of drought response genes bind TFs belonging to DREB 2 type protein [5, 31–33].
As mentioned, tolerance to a combination of different abiotic stresses is a well-known breeding target for sunflower as well as for other crops. Studies of simultaneous stress exposure were documented in various plant systems [8, 13, 34–41]. Nevertheless, little is known about the comparative molecular mechanisms underlying the acclimation responses of plants to a combination of different stresses .
Gene expression databases are increasing exponentially and the resulting information is stored and classified. While powerful software algorithms allow structural sequence similarity comparisons between species, difficulties arise to predict molecular function based on comparisons with homologous genes identified as determinant for a specific trait in different species. Identification of true orthologous among species is a powerful tool for candidate gene detection, particularly when comparing species having their full genome sequenced and those based on EST sequencing projects. In the case of Asteraceae species (including sunflower), only small syntenic fragments with Arabidopsis could be identified and their evolution involving major chromosomal rearrangements makes orthologous gene pairs difficult to identify . Even for other plant taxa comparative functional transcriptomic studies among crop plant genomes is relatively scarce . A microarray analysis in A. thaliana to identify simultaneously conserved and differentially expressed genes in oat, poplar and Euphorbia esula L. ("leafy spurge") was recently reported .
Sunflower was described as normally susceptible to low temperatures and salinity [45, 46], however, available information on gene expression in response to abiotic stresses are still limited to few studies [35, 45–49]. Recently the detection of a large number of down-regulated genes in plants exposed to extensive periods of low temperatures was reported , indicating that acclimation to chilling temperatures does not occur in sunflower. Meanwhile, transcriptional profiles in drought-tolerant and non-tolerant sunflower genotypes in response to water-stress allowed the identification of differential gene expression related to amino acid and carbohydrate metabolisms and signal-transduction processes . In this work we report for the first time a concerted study on gene expression in early responses to chilling and salinity using a fluorescence microarray assay based on organ-specific unigenes in sunflower. The aim of this work was to detect candidate genes associated to regulatory and stress-response pathways common to both stresses and at the same time identify those genes exclusively expressed in response to each kind of stress conditions.
Results and Discussion
Data analysis and accuracy of biological replicates
Differential expression of organ-specific sunflower sequences previously obtained by suppression subtractive hybridisation (SSH)  was evaluated with regard to the response to chilling and salt stresses using cDNA fluorescence microarray hybridization followed by Northern blot and qRT-PCR validation. Thus, three biological replicates were evaluated for chilling and salinity stresses as well as for control plants. The stress treatments were designed considering that sunflower has been described as normally susceptible to low temperatures and salinity [45, 46]. Regarding chilling tolerance, there is only one report studying sunflower response to low temperatures . However, in that study only long-term acclimatizing was evaluated, thus meaning that detected changes in gene expression reflect mainly plant metabolism adaptation to grow under suboptimal conditions more than short term responses to chilling. No previous studies on concerted gene expression of cultivated sunflower to salinity were reported before. In the present work, an early response to chilling and salinity is evaluated in order to detect early transcriptional changes in genes induced at the onset of the tolerance process.
Candidate gene selection
Data normalization - Normalization within microarrays
Principal axis 1, which retained 75.6% of the total variability, sorts genes according their fold change in over (on the right) and under (on the left) expression conditions independently of the stress condition applied. Principal axis 2 retaining 17.2% of total variability emphasizes the differences between cold and salinity treatments.
In order to experimentally validate differentially expressed genes derived from microarray analysis, Northern blot analysis and qRT-PCR were performed for candidate gene T411 [GenBank: BU671801], not only to validate the differential expression of this gene but also to set the cut-off significance criteria in the Bi-plot analysis within those genes retained by the p value criteria. This gene is up-regulated under cold and salinity stresses (see Additional files 2 and 3) and its position in the distance-to-the-origin distribution corresponds to the 70th percentile, as previously described.
Oligonucleotides used for qRT-PCR validations
GenBank (dbEST) Accession Number
Forward primer 5'-3'
Reverse primer 3'-5'
Comparison of gene expression levels obtained by cDNA microarray and qRT-PCR analysis for 10 differentially expressed genes
GenBank (dbEST) Accession Number
Fold change qRT-PCR
Fold change microarray
Differences in transcriptional changes detected in organ-specific cDNA libraries
Differences between cold and salinity treatments
Specific gene expression patterns in plants exposed to different treatments are evidenced when profile graphs by gene are analyzed (see Additional file 2). Drought, salinity and cold stress reduce water availability decreasing cell water potential. In order to avoid dehydration, plant responses include solute accumulation, cell wall components modification, and synthesis of protective proteins, avoiding or repairing cellular damage . The activation of these responses requires a complex signaling network being many of them shared by various abiotic stresses as those involving the DREB/CBF pathway [9, 54] and other ones typically from a determined abiotic stress [7, 55].
Although heatmaps are largely used in microarray analysis as a tool to visualize and detect differentially expressed genes  (see Additional file 3), an alternative and more reliable tool is the evaluation of individual gene transcription profile (see Additional file 2). Indeed, several reports on the usefulness of both methodologies to analyze microarray results have been recently published .
The availability of microarray technology allowing the comparison of large numbers of genes that are regulated under certain condition represents a powerful tool for many researchers. However, caution should be exercised when interpreting the outcome results. For example, salt stress treatment results in a rapid decline in photosynthesis within minutes . Therefore, genes that are regulated by photosynthetic activities maybe affected but they are not regulated by salt stress per se. In another study the expression of 150 genes in response to wounding and insect feeding was carried out , finding that some of the responding genes, for which the inducing stimulus was unclear in these processes, were also characterized as induced by drought . Thus, several factors could contribute to a complex pattern of transcript levels in which the interplay of stimuli that control gene expression can override each another. The stress intensity and the time course of gene induction are also important factors to be considered . The primary and the secondary stresses may induce genes in different time frames.
In Arabidopsis studies on concerted genes expression revealed a large number of cold- acclimation and freezing tolerance genes [6, 60–62], but there are only few studies on acclimation to chilling temperature . The ability of sunflower plants to gain a frost tolerance after exposure to a period of low temperature is still poorly understood. Although, recent results suggest that sunflower plants are non-acclimating plants under two conditions of low-term long-temperature at 15°C and 7°C . Thus, in the present study the response of sunflower to low temperatures has been focused at primary responses of young plants to a 24 hs treatment at 10°C with the aim to detect regulatory mechanisms induced at this early stage. Tolerance to low temperatures arise as an important trait considering that sunflower productive area is expanding to marginal regions with suboptimal growing conditions and the increasing requirement of early sow to maximize the growing season in many countries. Regarding tolerance to high NaCl contents, there is not much information available for sunflower but expansion of crops to more arid regions is associated to an increasing problem of soil salinity.
Low-temperature exposure in combination with high irradiance causes rapid inhibition of photosynthesis in a broad range of plants including tomato, cucumber and maize. Several elements contributing to this inhibition have been identified . Damage to the reducing side of photosystem II is well documented [67, 68] and, for moderately sensitive species such as maize, it may be the major cause of impaired whole plant photosynthesis following chilling. However, in the most chill-sensitive species, such as tomato, impaired reductive activation of the stromal biphosphatases appears to be the dominating factor limiting carbon assimilation following chilling in the light . Low temperature at night can also cause severe reductions in CO2 fixation on the day after chilling. In this work, genes that encode products with predicted functions related to energy metabolism were down-regulated under both stresses in sunflower. Among them, many genes potentially encoding components involved in photosynthesis, such as photosystem proteins, chlorophyll-binding proteins, Rubisco and light harvesting proteins showed differential expression patterns in this study. Considering sunflower as a medium-tolerant plant to chilling sensitivity , it has been suggested that decreasing energy metabolism is one of the cellular processes associated with the sunflower response to low temperatures . Here, we show that this process is not only down-regulated under chilling stress but also under salinity stress in H. annuus. Down-regulation of fructose-1,6 biphosfatase under drought stress in sunflower was reported as mainly associated to the stomata closure produced during water deficit . Yet, it was observed that NaCl reduces photosynthetic activity in Phaseolus vulgaris independent of stomata closure and by reducing the RuBP pool size through an effect on the RuBP regeneration potential . Thus, decline in photosynthesis in response to salinity has been attributed to the salinity effect on both stomatal and non-stomatal controls . However, the same enzyme from the halophytic wild rice, Portesia coarctata, decreases its catalytic activity by salt and may have intrinsic structural properties to withstand such decline .
Calvin cycle is mainly affected under dehydration and salinity stress in stress-sensitive plants as sunflower. However, the effect of a cold stress on this pathway is still unknown. It has been previously reported that during cold acclimation there is an increase in the availability of Pi and phosphorylated intermediates in both the pathway for Suc synthesis and the Calvin cycle, and increased capacities of enzymes in both pathways . However, one consequence of these long-term changes in cytosolic Pi availability and the capacity for Suc synthesis could be to pull too much carbon out of the chloroplast via the triose-phosphate transporter. This would, in turn, reduce the capacity of the Calvin cycle to regenerate RuBP and inhibit photosynthesis. However, cold-tolerant species such as A. thaliana and winter cereals are able to recover their photosynthetic capacity and resume growth after several days to weeks at low temperature through cold acclimation process [72–74], being actually little evidence available about sunflower response during the days after cold stress.
There are many reports on defense signals induced by different stresses including cold, salinity and drought [8, 35, 75–77]. In this category we report ten genes with differential expression patterns (see Additional file 1). One of them is up-regulated under both stresses (EST T411, similar to a plastidic aldolase) [GenBank: BU671801], a second one is down-regulated under chilling stress (EST T340, similar to a chloroplastic glutamine synthetase) [GenBank: BU671843] whereas a third one is specifically down-regulated under salt stress conditions. Plastidic aldolase genes characterized in Nicotiana plants can be grouped in two sub-families: AldP1 and AldP2. It was first reported that AldP2 was up-regulated by salt stress, whereas AldP1 was suppressed by salt stress . Thus, EST T411 identified in this work as up-regulated in both stresses would be hypothetically similar to an AldP2 type due to the up-regulation observed under salinity stress. On the other hand, down regulation of the transcriptional profile of EST H136 (similar to a chloroplast drought-induced stress protein) [GenBank: O04002] under chilling and salinity was observed. Typically, pathways leading to CO2 fixation and light harvesting are suppressed under abiotic stresses; although there is evidence of an over expression of glutamine synthase to enhance salinity tolerance in plants .
Signaling and transcription machinery
Regulatory proteins as TFs (bZIP, MYC, MYB and DREB) as well as protein kinases and proteases are involved in transcriptional changes under abiotic stresses . The activation of the transcriptional machinery in regulation of salt-dependent gene expression requires the induction of specific TFs as well as RNA polymerases . Many regulatory proteins, mainly identified in A. thaliana , showed changes in TFs under environmental stresses. In this analysis, and in agreement with a previous report , an EST encoding a protein with similarity to a zinc finger family protein was identified as up-regulated under low temperature, although this transcription factor (TF) does not show significant similarity to the one previously reported. A zinc finger protein associated to saline stress in Arabidopsis was recently reported  which is different from the one previously identified as DREB1A .
Many candidate genes had been identified as TFs or as sensing receptors for calcium signaling by in silico analysis, showing the relative abundance of transcriptional machinery related genes in the organ-specific cDNA libraries developed by our group . In the present work, two DNA binding proteins isolated by stress organ-specific cDNA libraries were detected as up-regulated specifically under salt (T187, T454) [GenBank: BU671817] [GenBank: BU671860] and a zinc finger protein specifically induced under chilling treatment. The large amount of TFs identified in Arabidopsis, indicating the complexity in the secondary metabolism of the plants , could explain the dramatic implication of those proteins in the abiotic stress responses beyond the critical interaction between plants and the environment and the level of duplications found in the Arabidopsis genome . A. thaliana TFs involved in stress response are traditional classified in ABA-dependent and ABA-independent regulatory pathways. According to microarray analysis in this species there are several independently responses to abiotic stress, one of them involving the DREB/CBF regulon . While DREB1 genes are specifically induced by cold, DREB2 genes are induced by dehydration and high salt but not by cold [7, 32, 84, 85]. This response was also reported in rice .
Recent studies have reported the importance of HD-ZIP TFs in response to drought in an independent pathway respect to the DREB transcriptional cascade . In Arabidopsis, ESTs libraries analysed by digital northern represented a 13% of signaling associated genes . The up-regulation of an ADP-ribosylation factor described in this work (EF127) [GenBank: BU671885] would explain the regulation of the intracellular traffic through vesicles . In addition, another gene under-regulated in both stresses (T234) [GenBank: BU671830] highly similar to an extracellular Ca2+ sensor was detected. That receptor was recently identified in A. thaliana . The authors demonstrated that Ca2+ extracellular level regulates Ca2+-dependent intracellular signaling through specific sensors. In this way those receptors would modulate calcium-dependant kinases previously described as enzymes highly expressed under chilling acclimation mainly. This knowledge was reported for Medicago sativa plants evaluated after 10 minutes exposition at low temperatures [90, 7]. The role of calcium-dependent protein-kinases under different environmental stresses was also reported in A. thaliana  and rice .
In general, genes involved in this cellular process were up-regulated under abiotic stress as a protective mechanism against key enzymes activity (see Additional file 1). Regulation of the translational machinery is considered an integral component in the cellular stress response [77, 93]. It has been indicated that ribosomal proteins are not only central to translational efficiency but have extra-ribosomal functions . In the present assay, this is confirmed by up-regulation of ribosomal proteins under both stresses. Besides, cDNAs encoding elongation factors were detected as salt induced, as previously reported for stress-associated genes in several systems [8, 95].
Protein turnover/folding/protein interactions
Protein degradation during stress is a highly conserved and regulated phenomenon in all the organisms reported so far . In this analysis, EST F443 [GenBank: BU671999] similar to a copper chaperone from tomato, was down-regulated under both stresses, as previously reported by a MALDI-TOF analysis of cold stress induction in rice seedlings . This tomato's chaperone seems to play a role in copper mobilization from decaying organs towards reproductive structures, contributing to growth in other parts of the plant . In addition, EST F231 [GenBank: BU671928] similar to a cyclophilin, was also down-regulated under chilling and salinity. Three peptidylprolyl isomerases (PPI) were detected in Arabidopsis plants treated with NaCl being one of them similar to a cyclophilin down-regulated under salinity stress .
By contrast, EST T124 [GenBank: BU671806] similar to a heat shock protein was up-regulated in both stresses. Heat shock proteins (HSPs), often called the stress proteins, are now recognized as important to a range of physiological and cellular functions under both normal growth conditions and in response to stresses other than heat shock . Starting in the mid-1980s, the concept of molecular chaperones evolved from the work of biochemists and cell biologists, and several HSPs were soon recognized as having such chaperoning functions. Proteomic analysis allowed the identification of several HSP's up-regulated in poplar under chilling stress, being one of them strongly similar to one of H. annuus detected in this work .
Cold stress, salinity and drought, combined with high light conditions, result in enhanced production of ROS by the photosynthetic apparatus because these conditions limit the availability of CO2 for the dark reaction, leaving oxygen as one of the main reductive products of photosynthesis . Drought, salt, and cold stress all induce the accumulation of ROS such as super oxide, hydrogen peroxide, and hydroxyl radicals . H2O2 is generated in peroxisomes by the enzymatic activity of glycolate oxidase . In this study an EST with homology to a glycolate oxidase (T120) [GenBank: BU671805] was up-regulated under both stresses, probably involved in a general generation of ROS under different abiotic stresses. These ROS may be signals inducing ROS scavengers and other protective mechanisms, as well as damaging agents contributing to stress injury in plants . Many ESTs from leaf and stem cDNA libraries encoding peroxidases, thioredoxins, catalases and oxygen-evolving enhancer proteins showed transcriptional changes in response to the studied stresses. Most of those proteins were up-regulated in both stresses due to the accumulation of these products along the oxidative stress. However, a NADH-plastoquinone reductase and a catalytic hydrolase were down-regulated. Genes encoding proteins associated with cellular homeostasis (respiration, cellular biogenesis and DNA repair) showed a distinct decline under abiotic stresses .
Ion homeostasis during salt stress is affected by sodium fluxes, transport and compartmentalization . Abundant transport-related genes have been described in a differential gene expression study that involve hybrid sunflower species, as preferentially expressed in Helianthus deserticola, a xerophytic species restricted to desert habitat . These genes seem to be important in the extreme environment of desert soil, functioning as both osmotic sensors and ionic regulators to prevent desiccation [103, 104].
In this work, EST H360 [GenBank: BU672086] with similarity to an ATP synthase was up-regulated in both stresses, as happens with EST F557 [GenBank: BU672042] similar to a putative carrier protein. These transcriptional changes could take place as a result of the large activity of ion transporters during salt tolerance and potassium nutrition . Potassium transporters may function in the transport of K+, which is an essential cofactor for many enzymes ; or control K+ uptake and regulate Na2+ uptake, which can be an important determinant of salinity tolerance [12, 106]. Moreover, carrier proteins, water-channel proteins and sugar transporters are thought to function through plasma membranes and tonoplast to adjust the osmotic pressure under stress conditions . On the other hand, lipid transfer proteins (LTPs) are another group of transport-related proteins associated to fatty-acid metabolism which may have a function in repairing stress-induced damage in membranes or changes in the lipid composition of the membranes, perhaps to regulate permeability to toxic ions and the fluidity of the membrane [107, 108]. Many LTPs have been shown to affect cell wall extensibility or to be secreted in response to NaCl stress . Nearly half of the detectable LTP transcripts in Arabidopsis root microarray were down-regulated by NaCl treatment . Moreover, LTPs with similarity to Arabidopsis' LTPs were also detected as down-regulated in sunflower under chilling stress . In the present work we identified an EST (EF502) highly similar to an LTP protein [GenBank: BU671910] being up-regulated under chilling stress and down-regulated under saline environment.
EST H123 [GenBank: BU672069], which shows a high identity to a myo-inositol phosphate synthase (MIPS protein, isomerase involved in inositol metabolism) [IUBMB enzyme nomenclature: EC 220.127.116.11.] was down-regulated in chilling and salinity. Inositol is a natural cell wall osmoprotector subcellular synthesized into phosphatidyl-inositol as part of a complex process and then recycled into phosphatidyl-inositol cycle as a complex signaling mechanism under abiotic stress conditions. Transgenic tobacco tolerant to salinity mediated by a MIPS gene product has been also described in P. coarctata . By contrast, in sesame, down-regulation by salt stress in seeds during germination was reported for the SeMIPS1 gene . In addition, transcription of the MIPS gene was found to be affected by salinity during biosynthesis of myo-inositol and its derivatives [112, 113], whereas evidence was reported that expression of the MIPS gene is up-regulated during salinity stress in salt-tolerant plants, while its transcriptional levels are reduced in salt-sensitive A. thaliana . The down-regulation of EST H123 [GenBank: BU672069] in sunflower (salt-sensitive crop) reported here is strongly similar to a MIPS gene product in sesame [GO Term GO: 0004512] .
This work presents the first cDNA sunflower fluorescence microarray analysed by a combined statistical method, studying transcriptional changes in early responses to chilling and salinity. The statistical approach to select candidate genes combining a classical hypothesis test of equal mean across treatment, i.e ANOVA, and an ordination technique based on principal component analysis appears as an efficient methodology to identify differentially expressed genes revealing a total of 80 candidate genes either under chilling and/or salt stress. Even when this represents a high percentage (28%) of differentially expressed genes from the initial number of organ-specific sequences evaluated, this is a lower proportion than that identified by the ANOVA which was about 50%. The reduction on the number of proposed differentially expressed genes by the combined selection criteria mentioned above is useful to prevent against false discoveries. Ten candidate genes from 12 selected ESTs representing different expression patterns were successfully validated by Real-Time quantitative PCR. Out of the 80 candidate genes, 50 genes were found up or down-regulated under abiotic stresses, supporting common regulatory mechanisms and general responses to low temperature and salinity. In addition, 15 and 12 genes were up or down-regulated respectively under one specific stress whereas the three remaining genes showed a contrasting transcriptional profile, being induced under one stress and suppressed under the other. Interestingly, almost half of the differentially expressed genes (39) detected in this study correspond to genes with unknown predicted function. This result indicates that even though ESTs database for Compositae plants comprises a large number of sequences (509,904), many of them do not have an assigned putative function. The difficulty in finding orthologous pairs in sequence comparisons with a fully sequenced genome species as Arabidopsis which is highly divergent from the Compositae species, could explain the large number of unknown or unclassified ESTs in expression studies involving partially sequenced genomes as sunflower.
There are many efforts from different research groups reporting the use of ESTs in microarrays analysis to contribute to the identification of those candidate genes whose expression level changed in presence of abiotic stresses [21, 35, 37, 61, 62, 115]. However, many EST collections are not complete and are derived from cDNA libraries of plants grown under non-stressed conditions, being defense/stress ESTs under represented. In the case of sunflower, functional genomic studies targeting different forms of water-deficit stress [49, 89] have been conducted using low scale thematic microarrays, though there are an important number of ESTs available for this species. In this sense, we consider that this work makes an essential contribution to the knowledge of an oil crop plant genome by its transcriptome characterization under two abiotic stresses: chilling and salinity. It represents the first work studying concerted gene expression of sunflower in response to salinity, allowing the identification of genes involved in common regulatory mechanisms to both stresses from those specifically associated to either chilling or salinity. Further studies exploring profile expression of these candidate genes under different low temperatures and salinity treatments combining different stress intensities and different stress extension periods will help to understand their role at different points of the complex regulatory mechanisms associated to stress response. Selected candidate genes could be used ultimate to test their molecular functions using expression studies for over expression or suppression of target genes in transgenic plants. These candidate sequences constitute at the same time a valuable tool to develop functional molecular markers based on SNPS/IndDels to assist selection in breeding programs.
Sunflower (H. annuus L.) plants belonging to inbred line HA89 provided by sunflower breeding program from EEA INTA Balcarce, Argentina were grown under standard conditions in greenhouse (16 h photoperiod and 20–24°C temperature) in pots of 1 liter (volume) containing composite soil. Three pots, representing three biological samples, each of them containing 4 seedlings, were grown for each treatment including control plants, watered daily with tap water and fertilized weekly with Hakaphos (COMPO®) 18-18-18 (NPK) at a final concentration of 100 ppm (0.55 W/V) during 2 weeks.
Chilling and salinity treatments
For high salinity treatment, 2-week old seedlings (2-full expanded leaves) were watered with 150 mM salt solution for three days (adapted from Liu and Baird) . Control seedlings were watered with tap water and grown under the same conditions. For chilling stress, sunflower seedlings at the two-leaf stage were subjected to 10°C (adapted from Huang et al.) . All plants were grown in growth chamber (Conviron®) during 24 hour with a 16 photoperiod provided by daylight fluorescent tubes (Philips, Argentina). Leaves from control and stressed seedlings were collected separately per biological replicate, immediately frozen in liquid nitrogen and stored at -80°C until processed for RNA extraction. Three biological samples represented by three pots, each one composed by four seedlings, were processed and evaluated independently for each stress treatment: chilling, salinity and control.
Amplification and preparation of cDNAs for microarrays construction
Sunflower EST clones derived from different organ-specific cDNA libraries were grouped in contigs using BioPipeline . Running BLASTN and BLASTX routines , 319 sunflower sequences with significant similarities (E value < 1.0 E-10 or BLAST score > 80) to already known or predicted genes involved in main stress/defense responses, primary metabolic pathways, gene expression or signal transduction were selected for microarray construction. BLAST-based GO term prediction application were used to assign functional categories to the defined unigene collection considered all three GO categories at the most specific term as described in a previous work . GO annotation updates were done through AmiGO .
EST libraries amplification
Recombinant plasmid preparations from different organ-specific cDNA libraries  that were previously used for 3' sequencing were used again for PCR amplification. LacZ1 and 2 forward and reverse primers were used to amplify the cloned inserts in 96 well plates using an Eppendorf thermocycler.
(LacZ1) 5'-3' sequence: GCT TCC GGC TCG TAT GTT GTG TG 5'-3'
(LacZ2) 5'-3'sequence: AAA GGG GGA TGT GCT GCA AGG CG.
Four PCR reaction plates were prepared, one from each cryogenic glycerol tube containing the cloned EST. Long Expand Template PCR System kit (Roche Diagnostics, Inc.) was used in a 50 μl final volume master mix consisting of final concentrations of 1× PCR buffer 3, dNTPs 0.2 mM, 0.25 mM Lac Z primers, Taq Pol Mix 0.75 u. Forty eight μl from this volume were aliquoted into each well of a 96 well PCR reaction plate (MJ Research). A 2 μl aliquot of an undiluted plasmid template DNA was aliquoted into the 48 μl of master mix. The plates were briefly centrifuged for 30 sec at 1500 rpm and placed into a Eppendorf thermocycler at denaturing conditions (94°C) for 1 min, followed by 32 cycles of 94°C for 30 sec; annealing temperatures were programmed at a descendant ramp from 60°C to 55°C and 72°C for 30 sec and a final extension of 72°C for 5 minutes. A typical yield from the PCR was about 50–150 ng/μl of amplified DNA.
The PCR products were loaded into a 96-well plate (Millipore #MANU3050) and vacuum filtered at 15 psi for about 10 min until the wells were completely empty. Purified products (QIAquick 96 PCR Purification Kit, Qiagen, Germany) were eluted in sterile water. One aliquot from each sample was evaluated by electrophoresis on 1% agarose electrophoresis gels to confirm amplification quality and quantity. Low DNA Mass Ladder (Invitrogen 10068013, Invitrogen, Argentina) was simultaneously run in the same electrophoresis gels to get a reliable PCR product quantification by fluorescence comparison using a Typhoon (GE Healthcare Life Sciences, Argentina) digitalization machine and software.
Spotting, microarray construction and post-print processing
Five μl form each of the purified PCR cDNA products (317 in total) were transferred into a 384 well plate containing 5 μl of DMSO 100% and spotted onto coated glass-slides (Ultragap II, Corning Systems, USA) by Gentron Genomic Services (Gentron, Buenos Aires, Argentina), using VersArray Chip Writer (BioRad, USA). The array design consisted in 4 supergrids, each containing 6 subgrids of 64 spots each (8 × 8), being each cDNA printed in quadruplicate. Three clones corresponding to house mouse (Mus musculus) were used as negative controls [GenBank: NM009060, NM008690, NM019476] while actin and rRNA sequences from sunflower (H. annuus) [GenBank: AAF82805] were used as a positive control for expression analysis in all microarray slides. The printed arrays were cross-linked to the slide by UV irradiation at 250 mJ using UV Stratalinker 2400 (Stratagene, USA). The slides were stored in a dessicator chamber until use. All the slides were hybridized with a pooled control RNA used as reference hybridization.
The microarray derived data platform was entered in The Gene Expression Omnibus database [117, 118] from which a platform accession number was assigned [GEO: GPL 4366]. Thus, complete tables of sequence identifiers and organ-specific unigenes accession numbers printed on arrays are available .
RNA, extraction, purification, amplification and labeling
Total RNA was extracted from approximately 2 g of leaf tissue using TRIzol® reagent following manufacturer recommendations (Invitrogen, Argentina). RNA integrity was analyzed by checking its electrophoretic mobility on 1.5% agarose gels in ME buffer (400 mM MOPS, 100 mM Na acetate, 10 mM EDTA pH 8.0, in diethyl-pyrocarbonate treated water). RNA was further purified by use of RNeasy Mini columns (Qiagen, Germany) according to manufacturer's instructions. To control biological variation between individuals, three biological samples from the same tissue were pooled on one sample prior to probe preparation. The reference (control) sample consisted of pooled RNA extracted from sunflower seedlings growing under unaltered environmental greenhouse conditions, whereas chilling and salinity samples were RNA extracted from sunflower seedlings growing in greenhouse under those stressed conditions.
The RNA (800 ng) samples were labeled by using SuperScript Indirect RNA Amplification System Kit (Invitrogen, Argentina) based on the method previously reported . Following RNA amplification (with the incorporation of UTP aminoallyl), labeled product was achieved by incubating with Cy3 or Cy5 esters in alkaline media.
Microarray hybridization reactions
The microarray slides were used in order to quantify the relative expression of ESTs in control and treated leaves by Cy3 and Cy5 hybridization technique. Dye-swaps were used to correct for differences in incorporation and fluorescent properties of both dyes, generating a number of 9 slides per experiment (three slides for control and three slides for each treatment) with a total number of 18 slides considering dye-swaps hybridizations. The microarray slides were prehybridized by incubation in 5× SSC, 0.1% SDS, 1% BSA at 42°C. In the next day, the cover slip was removed and the slide was washed once in 1× SSC, 0.2% SDS (prewarmed to 42°C); once in 0.2× SSC, 0.2% SDS at room temperature; and once in 0.1× SSC at room temperature. Washes were conduced with gentle shaking at 100 rpm for 5 minutes. Slides were subjected to low speed centrifugation for 2 min at 500 rpm to dry them.
Slide scanning and signal quantitation
The hybridized slides were scanned using a VersArray Chip Reader (BioRad, USA) scanner (two different channels for the two different dyes were used) at three different detector sensitivities. Image analysis and signal quantification were performed using free open source software Spotfinder , quantifying signal intensity for each spot. Then, data integration from multiple scanning processes was achieved.
Data normalization – Normalization within microarrays
Background subtraction was performed before calculating ratios. The elements with either printing or hybridization artifacts were flagged and discarded before analysis. Only spots with an intensity of at least 1.5 times above the local background in both channels were used for subsequent analysis. The outcoming data from each slide were then log transformed (using log base 2) and normalized using 3-D normalization (depending on spot intensity and it's location in the array) (Alvarez et al., unpublished data) using "The R statistical language" . Potential artifacts and false positives were eliminated and only those clones that exhibited similar expression patterns between the original hybridization and the corresponding dye swaps were selected for further analysis . A gene expression matrix was generated and its analysis was focused on differentially expressed genes.
Normalization between microarrays
Methodology used among biological replicates hypothesize that most of genes do not change their expression level among treatments. In this context, quantiles equalization or other normalization tools will not substantially modify the change to detect patterns of different expression levels. However, according to exploratory data analysis it was determined that an important fraction of the organ-specific ESTs did show some expression level differences among treatments. So, no additional normalization was applied under the risk of increasing the false positive or negative identification discovery rates. In order to have gene profiles ranged according to their experimental error, the gene expression matrix was scaled in a gene by gene basis, dividing by the common within-treatments standard deviation, thus generating the Gene Expression Matrix (as it is referred in this work).
Gene expression matrix analysis
The whole analysis related to gene expression matrix was performed using software Infostat 2006® . A two step procedure was used for candidate gene identification. First an analysis of variance for every gene was performed and only those genes having p-values lower that 5% were retained for complementary analysis. The analysis of variance was run for a fixed effect model under a complete random design. In a second step, the location of genes in the space spanned by the two first principal components of the gene expression matrix was used to filter genes and reduce the rate of positive false discoveries. The rationality behind this procedure is that the farther the gene is to the origin of the space, the larger the fold change expression level. Under this assumption, those genes located in the periphery of the resulting bi-plot should be the most dramatically involved in the responses to stress. The cut-off criterion was set as the distance-to-the-origin of the EST T411 (contig of the EST T111, AN: BU671801). The main basis of using this EST as cutting point was that it had been already experimentally validated as differentially expressed by Northern-blot and qRT-PCR. The position of EST T411 in the distance-to-the-origin distribution corresponds to the 70th percentile.
A graphical representation of the gene expressions among treatment conditions is presented as a heatmap plot (see Additional file 2). The average linkage method with Euclidean distance was used to generate the clustering relationships using the heatmap function in "The R statistical language" [, #241].
For northern blot analyses total RNA (20 ug) from leaves was fractionated on 1.5% agarosa-MOPS 1× gel and blotted onto nylon membranes (Hybond-N+, GE Healthcare Life Sciences, Amersham, Argentina). In all cases, membranes were cross-linked by UV illumination. Probes used for northern hybridization were prepared by random priming of the purified PCR products corresponding to EST candidate gene clone using [32P]-dCTP (NEN Perkin Elmer, USA). Hybridizations were performed at 42°C in the presence of 50% formamide (Ambion ULTRAhyb® Ultrasensitive Hybridization Buffer, Ambion, USA) and washes were also done at 42°C, according to provider's instructions. Exposed to BIOMAX MR Kodak X-ray films (KODAK, SIGMA Argentina) were scanned and analyzed with the TN-image program in a Typhoon device to calculate the relative signal intensities standardized with respect to rRNA and actin sequence from sunflower [GenBank: AAF82805] signals.
To confirm the results obtained from microarrays experiments, the transcript abundances of 10 differentially expressed ESTs were tested. Gene-specific primers were designed using Primer 3 . Oligonucleotide primer sequences are shown in Table 1. First-strand cDNA was reverse transcribed from 500 ng of DNase treated RNA according to manufacturer instructions (Invitrogen, Argentina). The reaction was performed in a 30-ul volume containing 15 ul QuantiTect™ SYBR® Green PCR (Qiagen, Germany), 300 nm of each primer and 1 μl of cDNA derived from RT product. The PCR reactions were run in an ABI PRISM 7000 HT Sequence Detection System (Applied Biosystems, USA) using the following program: 50°C for 2 min, 95°C for 10 min and 40 cycles of 95°C for 15 sec and 60°C for 1 min. Following PCR amplification, the reactions were subjected to a temperature ramp to create the dissociation curve, measured as changes in fluorescence readings as a function of temperature, allowing the detection of non-specific products. The dissociation program was 95°C for 15 sec, 60°C for 15 sec, followed by 20 min of slow ramp from 60°C to 95°C. Three replicates of each reaction were performed and actin sequence from sunflower [GenBank: AAF82805] was used as an internal control to normalize gene expression level. Negative control reactions using untranscribed RNA were also run to confirm absence of genomic DNA. Quantifying the relative changes in gene expression was performed using the 2-ΔΔCT method . Comparative results between qRT-PCR and microarray fold changes are presented in Table 2.
This research was supported by the ANPCyT/FONCYT; BID 1728 AC/AR PID 267 and PAV 137 and INTA-PE 243.540. Dr. P. Fernandez holds a post-doctoral fellowship from INTA, Luis Fernandez holds a technical position at INTA, MSc. Julio Di Rienzo is Professor of Statistics and Biometry at National University of Cordoba, Dr. R. Heinz and Dr. N. Paniego are career members of the Consejo Nacional de Investigaciones Científicas y Técnicas (CONICET, Argentina) and Dr. H.E. Hopp is a career member of the Comisión de Investigaciones Científicas de la Provincia de Buenos Aires (CIC) and Professor at the Facultad de Ciencias Exactas y Naturales, University of Buenos Aires (UBA).
- National Sunflower Association. [http://www.sunflowernsa.com/stats/table.asp?contentID=109&htmlID=74&submit170=View&submit.x=57&submit.y=12)]
- Connor D, Hall A: Sunflower physiology. Monograph No. 35. Madison, WI: ASA, CSSA, SSSA. 1997Google Scholar
- Miller J, Gulya T: Registration of four maintainer (HA 382 to HA 385) and four restorer (RHA 386 to RHA 389) sunflower germplasm lines. Crop Science. 1995, 35: 286-Google Scholar
- Paniego N, Heinz R, Fernandez P, Talia P, Nishinakamasu V, Hopp H: Sunflower. Genome Mapping and Molecular Breeding in Plants 2nd edition. Edited by: Kole C. Berlin Heidelberg: Springer-Verlag; 2007,: 153-177.Google Scholar
- Pradeep KA, Parinita A, Reddy M, Sopory Sudhir K: Role of DREB transcription factors in abiotic and biotic stress tolerance in plants. Plant Cell Reports. 2006, 25 (12): 1263-10.1007/s00299-006-0204-8.Google Scholar
- Seki M, Narusaka M, Abe H, Kasuga M, Yamaguchi-Shinozaki K, Carninci P, Hayashizaki Y, Shinozaki K: Monitoring the expression pattern of 1300 Arabidopsis genes under drought and cold stresses by using a full-length cDNA microarray. Plant Cell. 2001, 13 (1): 61-72. 10.1105/tpc.13.1.61.PubMedPubMed CentralGoogle Scholar
- Zhu JK: Plant salt tolerance. Trends Plant Sci. 2001, 6 (2): 66-71. 10.1016/S1360-1385(00)01838-0.PubMedGoogle Scholar
- Sahi C, Singh A, Kumar K, Blumwald E, Grover A: Salt stress response in rice: genetics, molecular biology, and comparative genomics. Funct Integr Genomics. 2006, 6 (4): 263-284. 10.1007/s10142-006-0032-5.PubMedGoogle Scholar
- Bray E: Plant responses to water deficit. Trends in Plant Science. 1997, 2: 48-54. 10.1016/S1360-1385(97)82562-9.Google Scholar
- Shinozaki K, Yamaguchi-Shinozaki K, Seki M: Regulatory network of gene expression in the drought and cold stress responses. Curr Opin Plant Biol. 2003, 6 (5): 410-417. 10.1016/S1369-5266(03)00092-X.PubMedGoogle Scholar
- Yamaguchi-Shinozaki K, Shinozaki K: Organization of cis- acting regulatory elements in osmotic- and cold stress-responsive promoters. Trends in Plant Science. 2005, 10: 1360-1385. 10.1016/j.tplants.2005.10.002.Google Scholar
- 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.PubMedGoogle Scholar
- Kanesaki Y, Suzuki I, Allakhverdiev SI, Mikami K, Murata N: Salt stress and hyperosmotic stress regulate the expression of different sets of genes in Synechocystis sp. PCC 6803. Biochem Biophys Res Commun. 2002, 290 (1): 339-348. 10.1006/bbrc.2001.6201.PubMedGoogle Scholar
- Sahi C, Agarwal M, Reddy MK, Sopory SK, Grover A: Isolation and expression analysis of salt stress-associated ESTs from contrasting rice cultivars using a PCR-based subtraction method. Theor Appl Genet. 2003, 106 (4): 620-628.PubMedGoogle Scholar
- Kore-eda S, Cushman MA, Akselrod I, Bufford D, Fredrickson M, Clark E, Cushman JC: Transcript profiling of salinity stress responses by large-scale expressed sequence tag analysis in Mesembryanthemum crystallinum. Gene. 2004, 341: 83-92. 10.1016/j.gene.2004.06.037.PubMedGoogle Scholar
- Marin K, Kanesaki Y, Los DA, Murata N, Suzuki I, Hagemann M: Gene expression profiling reflects physiological processes in salt acclimation of Synechocystis sp. strain PCC 6803. Plant Physiol. 2004, 136 (2): 3290-3300. 10.1104/pp.104.045047.PubMedPubMed CentralGoogle Scholar
- Sottosanto JB, Gelli A, Blumwald E: DNA array analysis of Arabidopsis thaliana lacking a vacuolar Na+/H+ antiporter: impact of AtNHX1 on gene expression. Plant J. 2004, 40 (5): 752-771. 10.1111/j.1365-313X.2004.02253.x.PubMedGoogle 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.PubMedPubMed CentralGoogle Scholar
- Shiozaki N, Yamada M, Yoshiba Y: Analysis of salt-stress-inducible ESTs isolated by PCR-subtraction in salt-tolerant rice. Theor Appl Genet. 2005, 110 (7): 1177-1186. 10.1007/s00122-005-1931-x.PubMedGoogle Scholar
- 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.PubMedGoogle Scholar
- Seki M, Narusaka M, Ishida J, Nanjo T, Fujita M, Oono Y, Kamiya A, Nakajima M, Enju A, Sakurai T, Satou M, Akiyama K, Taji T, Yamaguchi-Shinozaki K, Carninci P, Kawai J, Hayashizaki Y, Shinozaki K: Monitoring the expression profiles of 7000 Arabidopsis genes under drought, cold and high-salinity stresses using a full-length cDNA microarray. Plant J. 2002, 31 (3): 279-292. 10.1046/j.1365-313X.2002.01359.x.PubMedGoogle Scholar
- Chinnusamy V, Jagendorf A, Zhu J: Understanding and improving salt tolerance in plants. Crop Science. 2005, 45: 437-448.Google Scholar
- Wolfraim LA, Langis R, Tyson H, Dhindsa RS: cDNA sequence, expression, and transcript stability of a cold acclimation-specific gene, cas18, of alfalfa (Medicago falcata) cells. Plant Physiol. 1993, 101 (4): 1275-1282. 10.1104/pp.101.4.1275.PubMedPubMed CentralGoogle Scholar
- Gilmour SJ, Artus NN, Thomashow MF: cDNA sequence analysis and expression of two cold-regulated genes of Arabidopsis thaliana. Plant Mol Biol. 1992, 18 (1): 13-21. 10.1007/BF00018452.PubMedGoogle Scholar
- Nordin K, Heino P, Palva ET: Separate signal pathways regulate the expression of a low-temperature-induced gene in Arabidopsis thaliana (L.) Heynh. Plant Mol Biol. 1991, 16 (6): 1061-1071. 10.1007/BF00016077.PubMedGoogle Scholar
- Yamaguchi-Shinozaki K, Shinozaki K: Characterization of the expression of a desiccation-responsive rd29 gene of Arabidopsis thaliana and analysis of its promoter in transgenic plants. Mol Gen Genet. 1993, 236 (2–3): 331-340. 10.1007/BF00277130.PubMedGoogle Scholar
- Dunn MA, Hughes MA, Zhang L, Pearce RS, Quigley AS, Jack PL: Nucleotide sequence and molecular analysis of the low temperature induced cereal gene, BLT4. Mol Gen Genet. 1991, 229 (3): 389-394. 10.1007/BF00267460.PubMedGoogle Scholar
- Houde M, Dhindsa RS, Sarhan F: A molecular marker to select for freezing tolerance in Gramineae. Mol Gen Genet. 1992, 234 (1): 43-48.PubMedGoogle Scholar
- Oullet F, Vazquez-Tello A, Sarhan F: The wheat wcs120 promoter is cold-inducible in both monocotyledonous and dicotyledonous species. FEBS Letters. 1998, 423: 324-328. 10.1016/S0014-5793(98)00116-1.Google Scholar
- Ingram J, Bartels D: The Molecular Basis of Dehydration Tolerance in Plants. Annu Rev Plant Physiol Plant Mol Biol. 1996, 47: 377-403. 10.1146/annurev.arplant.47.1.377.PubMedGoogle Scholar
- Viswanathan C, Zhu JK: Molecular genetic analysis of cold-regulated gene transcription. Philos Trans R Soc Lond B Biol Sci. 2002, 357 (1423): 877-886. 10.1098/rstb.2002.1076.PubMedPubMed CentralGoogle Scholar
- Maruyama K, Sakuma Y, Kasuga M, Ito Y, Seki M, Goda H, Shimada Y, Yoshida S, Shinozaki K, Yamaguchi-Shinozaki K: Identification of cold-inducible downstream genes of the Arabidopsis DREB1A/CBF3 transcriptional factor using two microarray systems. Plant J. 2004, 38 (6): 982-993. 10.1111/j.1365-313X.2004.02100.x.PubMedGoogle Scholar
- Dubouzet JG, Sakuma Y, Ito Y, Kasuga M, Dubouzet EG, Miura S, Seki M, Shinozaki K, Yamaguchi-Shinozaki K: OsDREB genes in rice, Oryza sativa L., encode transcription activators that function in drought, high-salt and cold-responsive gene expression. Plant J. 2003, 33 (4): 751-763. 10.1046/j.1365-313X.2003.01661.x.PubMedGoogle Scholar
- 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.PubMedPubMed CentralGoogle Scholar
- Hewezi T, Leger M, El Kayal W, Gentzbittel L: Transcriptional profiling of sunflower plants growing under low temperatures reveals an extensive down-regulation of gene expression associated with chilling sensitivity. J Exp Bot. 2006, 57 (12): 3109-3122. 10.1093/jxb/erl080.PubMedGoogle Scholar
- Karrenberg S, Edelist C, Lexer C, Rieseberg L: Response to salinity in the homoploid hybrid species Helianthus paradoxus and its progenitors H. annuus and H. petiolaris. New Phytol. 2006, 170 (3): 615-629. 10.1111/j.1469-8137.2006.01687.x.PubMedPubMed CentralGoogle Scholar
- Kawasaki S, Borchert C, Deyholos M, Wang H, Brazille S, Kawai K, Galbraith D, Bohnert HJ: Gene expression profiles during the initial phase of salt stress in rice. Plant Cell. 2001, 13 (4): 889-905. 10.1105/tpc.13.4.889.PubMedPubMed CentralGoogle Scholar
- Renaut J, Lutts S, Hoffmann L, Hausman J: Responses of poplar to chilling temperatures: proteomic and physiological aspects. BMC Plant Biology. 2004, 6: 81-90.Google Scholar
- Cramer G, Ergül A, Grimplet J, Tillett R, Tattersall E, Bohlman M, Vincent D, Sonderegger J, Evans J, Osborne C, Quilici D, Schlauch K, Schooley D, Cushman J: Water and salinity stress in grapevines: early and late changes in transcript and metabolite profiles. Funct Integ Genomics. 2007, 7: 111-134. 10.1007/s10142-006-0039-y.Google Scholar
- Ozturk Z, Talamé V, Deyholos M, Michalowski C, Galbraith D, Gozukirmizi N, Tuberosa R, Bohnert H: Monitoring large-scale changes in transcript abundance in drought- and salt stressed barley. Plant Molecular Biology. 2002, 48: 551-573. 10.1023/A:1014875215580.Google Scholar
- Ueda A, Kathiresan A, Inada M, Narita Y, Nakamura T, Shi W, Takabe T, Bennett J: Osmotic stress in barley regulates expression of a different set of genes than salt stress does. Journal of Experimental Botany. 2005, 55: 2213-2218. 10.1093/jxb/erh242.Google Scholar
- Mittler R: Abiotic stress, the field environment and stress combination. Trends in Plant Science. 2006, 11 (1): 15-19. 10.1016/j.tplants.2005.11.002.PubMedGoogle Scholar
- Timms L, Jimenez R, Chase M, Lavelle D, McHale L, Kozik A, Lai Z, Heesacker A, Knapp S, Rieseberg L, Michelmore R, Kesseli R: Analyses of synteny between Arabidopsis thaliana and species in the Asteraceae reveal a complex network of small syntenic segments and major chromosomal rearrangements. Genetics. 2006, 173 (4): 2227-2235. 10.1534/genetics.105.049205.PubMedPubMed CentralGoogle Scholar
- Horvath D, Anderson J, Soto Suarez M, Chao W: Transcriptome analysis of "leafy spurge" (Euphorbia esula) crown buds during shifts in well-defined phases of dormancy. Weed Science. 2006, 54 (5): 821-827. 10.1614/WS-06-06R1.1.Google Scholar
- Kratsch H, Wise R: The ultrastructure of chilling stress. Plant Cell and Environment. 2000, 23: 337-350. 10.1046/j.1365-3040.2000.00560.x.Google Scholar
- Huang L, Ye Z, Bell RW, Dell B: Boron nutrition and chilling tolerance of warm climate crop species. Ann Bot (Lond). 2005, 96 (5): 755-767. 10.1093/aob/mci228.Google Scholar
- Lai Z, Livingstone K, Zou Y, Church SA, Knapp SJ, Andrews J, Rieseberg LH: Identification and mapping of SNPs from ESTs in sunflower. Theor Appl Genet. 2005, 111 (8): 1532-1544. 10.1007/s00122-005-0082-4.PubMedPubMed CentralGoogle Scholar
- Liu X, Baird W: Differential expression of genes regulated in response to drought or salinity stress in sunflower. Crop Science. 2003, 43: 678-687.Google Scholar
- Roche J, Hewezi T, Bouniols A, Gentzbittel L: Transcriptional profiles of primary metabolism and signal transduction-related genes in response to water stress in field-grown sunflower genotypes using a thematic cDNA microarray. Planta. 2007, 10.1007/s00425-007-0508-0.Google Scholar
- Fernandez P, Paniego N, Lew S, Hopp HE, Heinz RA: Differential representation of sunflower ESTs in enriched organ-specific cDNA libraries in a small scale sequencing project. BMC Genomics. 2003, 4 (1): 40-10.1186/1471-2164-4-40.PubMedPubMed CentralGoogle Scholar
- Livak K, Schmittgen T: Analysis of relative gene expression data using real-time quantitative PCR and the 2-ΔΔCt method. Methods. 2001, 25: 402-408. 10.1006/meth.2001.1262.PubMedGoogle Scholar
- Forment J, Gadea J, Huerta L, Abizanda J, Agusti S, Alamar E, Alos F, Andres R, Arribas J, Beltran A, Berbel M, Blazquez J, Brumos L, Canas M, Cercos J, Colmenero-Flores A, Conesa B, Estables M, Gandia J, Garcia-Martinez J, Gimeno A, Gisbert G, Gomez L, Gonzalez-Candelas A, Granell J, Guerri M, F L, Madueno J, Marcos M, Marques F, Martinez M, Martinez-Godoy S, Miralles P, Moreno L, Navarro V, Pallas M, Perez-Amador J, Perez-Valle C, Pons I, Rodrigo P, Rodriguez C, Royo R, Serrano G, Soler F, Tadeo M, Talon J, Terol M, Trenor L, Vaello O, Vicente C, Vidal L, Zacarias L, Conejero V: Development of a citrus genome-wide EST collection and cDNA microarray as resources for genomic studies. Plant Molecular Biology. 2005, 57: 375-391. 10.1007/s11103-004-7926-1.PubMedGoogle Scholar
- Verslues PE, Bray EA: Role of abscisic acid (ABA) and Arabidopsis thaliana ABA-insensitive loci in low water potential-induced ABA and proline accumulation. J Exp Bot. 2006, 57 (1): 201-212. 10.1093/jxb/erj026.PubMedGoogle Scholar
- Chinnusamy V, Schumaker K, Zhu JK: Molecular genetic perspectives on cross-talk and specificity in abiotic stress signaling in plants. J Exp Bot. 2004, 55 (395): 225-236. 10.1093/jxb/erh005.PubMedGoogle Scholar
- 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.PubMedPubMed CentralGoogle Scholar
- Allison D, Xiangquin C, Grier P, Sabripour M: Microarray data analysis: from disarray to consolidation and consensus. Nature Reviews Genetics. 2006, 7: 55-65. 10.1038/nrg1749.PubMedGoogle Scholar
- Lemaire-Chamley M, Petit J, Garcia V, Just D, Baldet P, Germain V, Fagard M, Mouassite M, Cheniclet C, Rothan C: Changes in transcriptional profiles are associated with early fruit tissue specialization in tomato. Plant Physiol. 2005, 139 (2): 750-769. 10.1104/pp.105.063719.PubMedPubMed CentralGoogle Scholar
- Reymond P, Weber H, Damond M, Farmer EE: Differential gene expression in response to mechanical wounding and insect feeding in Arabidopsis. Plant Cell. 2000, 12 (5): 707-720. 10.1105/tpc.12.5.707.PubMedPubMed CentralGoogle Scholar
- Xiong L, Schumaker KS, Zhu JK: Cell signaling during cold, drought, and salt stress. Plant Cell. 2002, S165-183. 14 SupplPubMedPubMed CentralGoogle Scholar
- Ma S, Gong Q, Bohnert J: Dissecting salt stress pathways. Journal of Experimental Botany. 2005, 57 (5): 1097-1107. 10.1093/jxb/erj098.Google Scholar
- Jiang Y, Deyholos MK: Comprehensive transcriptional profiling of NaCl-stressed Arabidopsis roots reveals novel classes of responsive genes. BMC Plant Biol. 2006, 6: 25-10.1186/1471-2229-6-25.PubMedPubMed CentralGoogle Scholar
- Kreps JA, Wu Y, Chang HS, Zhu T, Wang X, Harper JF: Transcriptome changes for Arabidopsis in response to salt, osmotic, and cold stress. Plant Physiol. 2002, 130 (4): 2129-2141. 10.1104/pp.008532.PubMedPubMed CentralGoogle Scholar
- Lee J, Prasad V, Yang P, Wu J, David Ho T, Charng Y, Chan M: Expression of Arabidopsis CBF1 regulated by an ABA/stress inducible promoter in transgenic tomato confers stress tolerance without affecting yield. Plant Cell and Environment. 2003, 26 (7): 1181-10.1046/j.1365-3040.2003.01048.x.Google Scholar
- Camon E, Magrane M, Barrell D, Binns D, Fleischmann W, Kersey P, Mulder N, Oinn T, Maslen J, Cox A, Apweiler R: The Gene Ontology Annotation (GOA) project: implementation of GO in SWISS-PROT, TrEMBL, and InterPro. Genome Res. 2003, 13 (4): 662-672. 10.1101/gr.461403.PubMedPubMed CentralGoogle Scholar
- Ashburner M, Ball CA, Blake JA, Botstein D, Buttler H, Cherry JM, Davis AP, Dolinsky K, Dwight SS, Eppig JT: Gene Ontology tool for the unification of biology. Nature Genetics. 2000, 25: 25-29. 10.1038/75556.PubMedPubMed CentralGoogle Scholar
- Martino-Catt S, Ort DR: Low temperature interrupts circadian regulation of transcriptional activity in chilling-sensitive plants. Proc Natl Acad Sci U S A. 1992, 89 (9): 3731-3735. 10.1073/pnas.89.9.3731.PubMedPubMed CentralGoogle Scholar
- Kee S, Martin B, Ort R: The effects of chilling in the dark and in the light on photosynthesis of tomato: electron transfer reactions. Photosynthesis research. 1986, 8: 41-51. 10.1007/BF00028475.PubMedGoogle Scholar
- Percival M, Bradbury M, Hayden D, Baker N: Modification of the photochemical apparatus in maize by photoinhibitory stress at low temperature4th edition. Dordrecht, The Netherlands: Martinus Nijhoff; 1987,:Google Scholar
- Sassenrath G, Ort D, Portis A: Impaired reductive activation of stromal biphosphatases in tomato leaves following low-temperature exposure at high light. Arch Biochem Biophys. 1991, 202: 302-308.Google Scholar
- Walker R, Torofalvy E, Subton W: Photosynthetic responses of Citrus varieties Rangpur lime and Etrog Citron to salt treatment. Australian Journal of Plant Pathology. 1982, 9: 783-790.Google Scholar
- Ghosh S, Bagchi S, Majumder A: Chloroplast fructose-1,6-biphosphatase from Oryza differs in salt tolerance property from the Porteresia enzyme and is protected by osmolytes. Plant Science. 2001, 160: 1171-1181. 10.1016/S0168-9452(01)00361-2.PubMedGoogle Scholar
- Hurry V, Hunner N: Low growth temperature effects a differential inhibition of photosynthesis in spring and winter-wheat. Plant Physiology. 1991, 93: 491-497.Google Scholar
- Levitt J: Chilling, freezing, and high temperature stresses. New York. 1980Google Scholar
- Strand Å, Hurry V, Gustafsson P, Gardeström P: Development of Arabidopsis thaliana leaves at low temperatures releases the suppression of photosynthesis and photosynthetic gene expression despite the accumulation of soluble carbohydrates. Plant Journal. 1997, 12: 605-614. 10.1046/j.1365-313X.1997.00605.x.PubMedGoogle Scholar
- Hasegawa P, Bressan R, Zhu J, Bohnert H: Plant cellular and molecular responses to high salinity. Annual Review of Plant Physiology and Plant Molecular Biology. 2000, 51: 463-499. 10.1146/annurev.arplant.51.1.463.PubMedGoogle Scholar
- Chen W, Provart N, Glazebrook J, Katagiri F, Chang H, Eulgem T, Mauch F, Luan S, Zou G, Whitham S, Budworth P, Tao Y, Xie Z, Chen X, Lam S, Kreps J, Harper J, Si-Ammour A, Mauch-Mani B, Heinlein M, Kobayash iK, Hohn T, Dangl J, Wang X, Zhu T: Expression profile matrix of Arabidopsis transcription factor genes suggests their putative functions in response to environmental stresses. The Plant Cell. 2002, 14: 559-574. 10.1105/tpc.010410.PubMedPubMed CentralGoogle Scholar
- Dhanaraj A, Alkharouf N, Beard H, Chouikha I, Matthews B, Wei H, Arora R, Rowland L: Major differences observed in transcript profiles of blueberry during cold acclimation under field and cold room conditions. Planta. 2006, 10.1007/s00425-006-0382-1.Google Scholar
- Yamada S, Toshiyuki K, Akiko H, Kuwata S, Hidemasa I, Tomoaki K: Differential expression of plastidic aldolase genes in Nicotiana plants under salt stress. Plant Science. 2000, 154: 61-69. 10.1016/S0168-9452(00)00188-6.PubMedGoogle Scholar
- Hoshida H, Tanaka Y, Hibino T, Hayashi Y, Tanaka A, Takabe T: Enhanced tolerance to salt stress in transgenic rice that overexpresses chloroplast glutamine synthetase. Plant Molecular Biology. 2000, 43 (103-111):PubMedGoogle Scholar
- Estruch F: Stress-controlled transcription factors, stress-induced genes and stress tolerance in budding yeast. FEMS Microbiol Rev. 2000, 25: 469-486. 10.1111/j.1574-6976.2000.tb00551.x.Google Scholar
- Rodriguez Milla M, Townsend J, Feng Chang I, Cushman J: The Arabidopsis AtDi19 gene family encodes a novel type of Cys2/His2 zinc-finger protein implicated in ABA-independent dehydration, high-salinity stress and light signaling pathways. Plant Molecular Biology. 2006, 61: 13-30. 10.1007/s11103-005-5798-7.Google Scholar
- Szathmáry E, Jordán F, Pál C: Can genes explain biological complexity?. Science. 2001, 292: 1315-1316. 10.1126/science.1060852.PubMedGoogle Scholar
- Arabidopsis GI: Analysis of the genome sequence of the flowering plant Arabidopsis thaliana. Nature. 2000, 408: 796-815. 10.1038/35048692.Google Scholar
- Liu J, Zhu JK: A calcium sensor homolog required for plant salt tolerance. Science. 1998, 280 (5371): 1943-1945. 10.1126/science.280.5371.1943.PubMedGoogle Scholar
- Nagashima A, Hanaoka M, Motohashi R, Seki M, Shinozaki K, Kanamaru K, Takahashi H, Tanaka K: DNA microarray analysis of plastid gene expression in an Arabidopsis mutant deficient in a plastid transcription factor sigma, SIG2. Biosci Biotechnol Biochem. 2004, 68 (3): 694-704. 10.1271/bbb.68.694.PubMedGoogle Scholar
- Manavella PA, Arce AL, Dezar CA, Bitton F, Renou JP, Crespi M, Chan RL: Cross-talk between ethylene and drought signaling pathways is mediated by the sunflower Hahb-4 transcription factor. Plant J. 2006, 48 (1): 125-137. 10.1111/j.1365-313X.2006.02865.x.PubMedGoogle Scholar
- Mahalingam R, Gomez-Buitrago A, Eckardt N, Shah N, Guevara-Garcia A, Day P, Raina R, Fedoroff Nina V: Characterizing the stress/defense transcriptome of Arabidopsis. Genome Biology. 2003, 4: R20-10.1186/gb-2003-4-3-r20.PubMedPubMed CentralGoogle Scholar
- Song X, Yang C, Liu J, Yang W: RPA, a class II ARFGAP protein, activates ARF1 and U5 and plays a role in root hair development in Arabidopsis. Plant Physiology. 2006, 141 (3):
- Han S, Tang R, Anderson L, Woerner T, Pei Z: A cell surface receptor mediates extracellular Ca2+ sensing in guard cells. Nature. 2003, 425: 196-200. 10.1038/nature01932.PubMedGoogle Scholar
- 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.PubMedPubMed CentralGoogle Scholar
- Choi H, Park HJ, Park JH, Kim S, Im MY, Seo HH, Kim Y-W, Hwang I, Kim SY: Arabidopsis Calcium-Dependent Protein Kinase AtCPK32 interacts with ABF4, a transcriptional regulator of abscisic acid-responsive gene expression, and modulates its activity. Plant Physiology. 2005, 139 (4):
- Wan B, Lin Y, Mou T: Expression of rice Ca(2+)-dependent protein kinases (CDPKs) genes under different environmental stresses. FEBS Letters. 2007, 581 (6): 1179-1189. 10.1016/j.febslet.2007.02.030.PubMedGoogle Scholar
- Wood A, Oliver M: Translational control in plant stress: formation of messenger ribonucleoprotein complexes (mRNPs) in Tortula ruralis in response to desiccation. Plant Journal. 1999, 18: 359-370. 10.1046/j.1365-313X.1999.00458.x.Google Scholar
- Wool I: Extraribosomal functions of ribosomal proteins. Trends in Biochemical Sciences. 1996, 21: 164-165. 10.1016/0968-0004(96)20011-8.PubMedGoogle Scholar
- Harding H, Calfon M, Urano F, Novoa I, Ron D: Transcriptional and translational control in the mammalian unfolded protein response. Annu Rev Cell Dev Biol. 2002, 18: 575-599. 10.1146/annurev.cellbio.18.011402.160624.PubMedGoogle Scholar
- Vierstra R, Callis J: Polypeptide tags, ubiquitous modifiers for plant protein regulation. Plant Molecular Biology. 1999, 41: 435-442. 10.1023/A:1006323317890.PubMedGoogle Scholar
- Cui S, Huang F, Wang J, Ma X, Cheng Y, Liu J: A proteomic analysis of cold stress response in rice seedlings. Proteomics. 2005, 5: 3162-3172. 10.1002/pmic.200401148.PubMedGoogle Scholar
- Mira H, Martínez-García F, Peñarrubia L: Evidence for the plant-specific intercellular transport of the Arabidopsis copper chaperone CCH. Plant Journal. 2001, 25: 521-528. 10.1046/j.1365-313x.2001.00985.x.PubMedGoogle Scholar
- Grover A: Molecular biology of stress responses. Cell Stress and Chaperones. 2002, 1: 1-5. 10.1379/1466-1268(2002)007<0001:MBOSR>2.0.CO;2.Google Scholar
- Mittler R: Oxidative stress, antioxidants, and stress tolerance. Plant Science. 2002, 7: 405-410.Google Scholar
- Prasad T, Anderson M, Martin B, Stewart C: Evidence for chilling-induced oxidative stress in maize seedlings and a regulatory role for hydrogen peroxide. The Plant Cell. 1994, 6 (1): 65-74. 10.2307/3869675.PubMedPubMed CentralGoogle Scholar
- Lai Z, Gross BL, Zou Y, Andrews J, Rieseberg LH: Microarray analysis reveals differential gene expression in hybrid sunflower species. Mol Ecol. 2006, 15 (5): 1213-1227. 10.1111/j.1365-294X.2006.02775.x.PubMedPubMed CentralGoogle Scholar
- Jang J, Kim D, Kim Y, Kim J, Kang H: An expression analysis of a gene family encoding plasma membrane aquaporins in response to abiotic stresses in Arabidopsis thaliana. Plant Molecular Biology. 2004, 54: 713-725. 10.1023/B:PLAN.0000040900.61345.a6.PubMedGoogle Scholar
- Zhu J: Salt and drought stress signal transduction in plants. Annu Rev Plant Biol. 2002, 53: 247-273. 10.1146/annurev.arplant.53.091401.143329.PubMedPubMed CentralGoogle Scholar
- Blumwald E, Grover A: Salt tolerance. UK: John Wiley and Sons Ltd; 2006.Google Scholar
- Walia E, Wilson C, Zeng L, Ismail A, Condamine P, Close T: Genome-wide transcriptional analysis of salinity stressed japonica and indica rice genotypes during panicle initiation stage. Plant Molecular Biology. 2007, 63: 609-623. 10.1007/s11103-006-9112-0.PubMedPubMed CentralGoogle Scholar
- Torres-Schumann S, Godoy J, Pintor-Toro J: A probable lipid transfer protein gene is induced by NaCl stems of tomato plants. Plant Molecular Biology. 1992, 18: 749-757. 10.1007/BF00020016.PubMedGoogle Scholar
- Holmberg N, Bülow L: Improving stress tolerance in plants by gene transfer. Trends in Plant Science. 1998, 3: 61-66. 10.1016/S1360-1385(97)01163-1.Google Scholar
- Nieuwland J, Feron R, Huisman B, Fasolino A, Hilbers C, Derksen J, Mariani C: Lipid transfer proteins enhance cell wall extension in tobacco. The Plant Cell. 2005, 7: 2007-2019.Google Scholar
- Das-Chatterjee A, Goswami L, Maitra S, Ghosh Dastidar K, Ray S, Majumder A: Introgression of a novel salt-tolerant L-myo-inositol 1-phosphate synthase from Porteresia coarctata (Roxb.) Tateoka (PcINO1) confers salt tolerance to evolutionary diverse organisms. FEBS Letters. 2006, 580: 3980-3988. 10.1016/j.febslet.2006.06.033.PubMedGoogle Scholar
- Chun J-A, Jin U-H, Lee J-W, Yi Y-B, Hyung N-I, Kang M-H, Pyee J-H, Suh M, Kang C-W, Seo H-Y, Lee S-W, Chung C-H: Isolation and characterization of a myo-inositol 1-phosphate synthase cDNA from developing sesame (Sesamum indicum L.) seeds: functional and differential expression, and salt-induced transcription during germination. Planta. 2003, 216: 874-880.PubMedGoogle Scholar
- Nelson D, Rammesmayer G, Bohnert H: Regulation of cell-specific inositol metabolism and transport in plant salinity tolerance. The Plant Cell. 1998, 10: 753-764. 10.2307/3870662.PubMedPubMed CentralGoogle Scholar
- Nelson D, Koukoumanos M, Bohnert H: Myo-inositol-dependent sodium uptake in ice plant. Plant Physiology. 1999, 119: 165-172. 10.1104/pp.119.1.165.PubMedPubMed CentralGoogle Scholar
- Ishitani M, Majumder A, Bornhouser A, Michalowski C, Jensen R, Bohnert H: Coordinate transcriptional induction of myo-inositol metabolism during environmental stress. Plant Journal. 1996, 9: 537-548. 10.1046/j.1365-313X.1996.09040537.x.PubMedGoogle Scholar
- Terol J, Conesa A, Colmenero J, Cercos M, Tadeo F, Agustí J, Alós E, Andres F, Soler G, Brumos J, Iglesias D, Stefan G, Legaz F, Argout X, Courtois B, Ollitrault P, Dossat C, Wincker P, Morillon R, Talon M: Analysis of 13,000 unique Citrus clusters associated with fruit quality, production and salinity tolerance. BMC Genomics. 2007, 8 (31):
- Altschul S, Gish W, Miller W, Myers E, Lipman D: Basic local alignment search tool. Journal of Molecular Biology. 1990, 215: 403-410.PubMedGoogle Scholar
- Edgar R, Domrachev M, Lash A: Gene Expression Omnibus: NCBI gene expression and hybridization array data repository. Nucleic Acid Research. 2002, 30 (1): 207-210. 10.1093/nar/30.1.207.Google Scholar
- Barrett T, Troup D, Wilhite S, Ledoux P, Rudnev D, Evangelista C, Kim I, Soboleva A, Tomashevsky M, Edgar R: NCBI GEO: mining tens of millions of expression profiles-database and tools update. Nucleic Acid Research. 2006, 760-765. 35 DatabaseGoogle Scholar
- Eberwine J, Spencer C, Miyashiro K, Mackler S, Finnell R: Complementary DNA synthesis in situ: methods and applications. 1992Google Scholar
- Saeed A, Sharov V, White J, Li J, Liang W, Bhagabati N, Braisted J, Klapa M, Currier T, Thiagarajan M, Sturn A, Snuffin M, Rezantsev A, Popov D, Ryltsov A, Kostukovich E, Borisovsky I, Liu Z, Vinsavich A, Trush V, Quackenbush J: TM4: a free, open-source system for microarray data management and analysis. Biotechniques. 2003, 34 (2): 374-378.PubMedGoogle Scholar
- "The R statistical language". [http://cran.r-project.org/faqs.html]
- Yang Y, Dudoit S, Luu P, Lin D, Peng V, Ngai J, Speed T: Normalization for cDNA microarray data: a robust composite method addressing single and multiple slide systematic variation. Nucleic Acid Research. 2002, 30 (4): e15-10.1093/nar/30.4.e15.Google Scholar
- Infostat 2006®. [http://www.infostat.com.ar]
- Rozen S, Skaletsky H: Primer3 on the WWW for general users and for biologist programmers. Totowa, NJ: Humana Press; 2000.Google Scholar