Candidate genes associated with bud dormancy release in blackcurrant (Ribes nigrumL.)
© Hedley et al; licensee BioMed Central Ltd. 2010
Received: 12 March 2010
Accepted: 14 September 2010
Published: 14 September 2010
The detrimental effects of mild winter temperatures on the consistency of cropping of blackcurrant (Ribes nigrum L.) in parts of Europe have led to increasing interest in the genetic control of dormancy release in this species. This study examined patterns of gene expression in leaf buds of blackcurrant to identify key differential changes in these profiles around the time of budbreak.
Using leaf bud tissue of blackcurrant, a cDNA library was generated as a source of blackcurrant ESTs for construction of a custom microarray, which was used to identify differential gene expression during dormancy release. Gene activity was lowest in early stages of dormancy, increasing to reach a maximum around the time of budbreak. Genes with significantly changing expression profiles were clustered and evidence is provided for the transient activity of genes previously associated with dormancy processes in other species. Expression profiling identified candidate genes which were mapped onto a blackcurrant genetic linkage map containing budbreak-related QTL. Three genes, which putatively encode calmodulin-binding protein, beta tubulin and acetyl CoA carboxylase respectively, were found to co-localise with budbreak QTL.
This study provides insight into the genetic control of dormancy transition in blackcurrant, identifying key changes in gene expression around budbreak. Genetic mapping of ESTs enabled the identification of genes which co-localise with previously-characterised blackcurrant QTL, and it is concluded that these genes have probable roles in release of dormancy and can therefore provide a basis for the development of genetic markers for future breeding deployment.
Blackcurrant (Ribes nigrum L.) is grown widely throughout temperate areas of Europe and New Zealand, with an annual production of 160,000 tonnes in Europe and 185,000 tonnes globally . The centres of genetic diversity for blackcurrant are situated in northern Scandinavia and Russia, although species are also found in North and South America, Asia and northwest Africa . The fruit has very high nutritional value, in terms of the content of ascorbic acid and other antioxidants , and its primary commercial use is processing for juice. In the UK, varieties are chosen to give a spread of cropping seasons. Timing of leaf budbreak is an important physiological factor in the developmental processes leading to cropping of blackcurrant, as uniformity of development is essential for both successful harvesting and for optimum fruit quality and yield.
The control of budbreak in many woody perennial plants, including berry fruits such as Ribes and Rubus, is dependent on exposure to chilling temperatures during the winter for a sufficient duration to release dormancy, followed by appropriate warmer temperatures in the spring to induce growth . An insufficient amount of chilling during the dormant period can seriously impact key phenological traits, notably time of budbreak, time and duration of flowering and fruit quality at harvest . The chilling requirement, in terms of both duration and range of suitable temperatures, varies between species, between cultivars within species [6, 7], and even between different buds of the same plant [8, 9]. Blackcurrant has a relatively high chilling requirement for dormancy break, from ca. 1,300 h in the case of New Zealand-selected germplasm  to over 2,000 h for some late-flowering types such as 'Ben Lomond' from Scotland and most Nordic cultivars . As warmer winters are likely to become increasingly prevalent , there is a risk of insufficient chilling becoming more widespread in cultivars requiring high levels of winter chilling, so breeding strategies to address this using low-chill parental material are now required .
However, these must be balanced against the risk of spring frost damage at flowering in cultivars that break bud too early after minimal chilling and whilst the incidence of spring frosts is declining across Europe , there can still be occasional damage in some areas. Also, the risk of spring frost damage at flowering may depend on temperature variance as well as on specific mean temperatures, and climate projections are predicting a future rise in such variance .
Dormancy in woody perennials is an adaptive mechanism for the survival of winter conditions  and was classified by Lang  into three major types: (i) ecodormancy, where growth inactivity is caused by unfavourable environmental conditions, with a resumption of growth when conditions improve; (ii) endodormancy, caused by endogenous factors within the dormant buds that cannot be overcome even by favourable environmental conditions; and (iii) paradormancy, or correlative inhibition , where the dormancy status of buds is influenced by physiological factors in other parts of the plant (eg. apical dominance). The effects of winter chilling on dormancy break at bud-burst are in endodormant tissues, and further investigation of winter chilling in blackcurrant is structured around improved phenotyping of diverse germplasm  coupled with genetic analysis of Quantitative Trait Loci (QTL) and genes associated with dormancy-related processes. Putative QTL linked to key developmental traits, including time of budbreak, have recently been identified by Brennan et al. .
Dormancy break and subsequent physiological events are controlled through the coordinated action of large numbers of genes in woody plants  and whilst there have been many studies of the molecular genetics of dormancy and budbreak, many aspects remain unclear. Microarray studies to examine global gene expression during dormancy and budbreak have been used in a range of woody genera, including Vitis, [16–18], Quercus  and Rubus . Most of these studies attempted functional classification of the differentially expressed genes from the transcriptome related to budbreak or, in the case of Mathiason et al. , to the fulfilment of chilling requirement. Some studies, eg. Mazzitelli et al.  and Derory et al. , have identified candidate genes for dormancy-related metabolism, and several of these genes were previously implicated in these processes in other species.
The aim of this study was to examine patterns of gene expression in leaf buds of blackcurrant and to identify the key differential changes in these profiles around the time of budbreak. A custom cDNA microarray was utilised which incorporated blackcurrant Expressed Sequence Tags (ESTs) derived from buds at various stages of dormancy. Expression profiling identified candidate genes associated with budbreak which were mapped onto a blackcurrant genetic linkage map. Initial steps have been taken to develop associated markers that can be deployed in the characterisation of diverse germplasm for inclusion in downstream breeding programmes.
Blackcurrant expressed sequence tags
RNA was isolated from blackcurrant leaf bud tissue sampled at five time-points prior to budbreak and, following pooling, was used to construct a targeted cDNA library to maximise collation of relevant ESTs. Single-pass sequencing of approximately 9,000 clones generated 7,400 high-quality ESTs. Following identification of a low-redundancy set of blackcurrant clones, single ESTs (3,633 in total) representing each contiguous sequence were selected and whole cDNA inserts used to fabricate microarrays as described.
Gene expression profiling during dormancy transition
Leaf buds were sampled from blackcurrant plants throughout the different dormancy phases, from mid-winter following onset of dormancy through to spring when bud-burst was prevalent. RNA was isolated from bud samples from winter 2005 to spring 2006: week 0 (early December), week 2, week 4, week 8, week 12, week 14 and week 16 (end of March). For the microarray experimental design, RNAs from consecutive time-points were hybridised on the same arrays in a two-channel loop-design (design and complete dataset available at ArrayExpress  accession E-TABM-807). Data were extracted from each array and quality filtering was applied to leave 2,320 reliable replicated expressed probes.
In order to identify and characterise variation within the entire reliable dataset Principal Component Analysis (PCA) was used (Figure 1a). Principal component 1 accounted for 66% of the variance, whereas principal component 2 comprised 12% of the variance, based upon the filtered genelist of 2,320 probes. Principal component 1 separates the samples of the earlier weeks (0, 2, 4 and 8) from those of the later weeks (12, 14 and 16). Principal component 2 in particular separates the week 8 data from the other time points. To analyse this further, an overview of gene activity in blackcurrant leaf bud tissue was obtained using volcano plots to identify numbers of differentially expressed genes between consecutive sample time points (Figure 1b). Early winter (week 2 to week 4) shows the lowest levels of differential gene expression, and in the first two weeks of sampling the majority of genes are down-regulated. Differential gene expression reaches a peak between week 8 and week 12, with similar numbers of up-and down-regulated genes detected. Subsequently, the majority of genes from week 12 to week 14 are reduced in expression levels.
Significant differentially expressed genes were identified over the entire period of sampling using ANalysis Of VAriance (ANOVA). False Discovery Rate (FDR) was applied at high stringency multiple testing correction; identifying 1,040 significantly changing probes (Additional File 1), of which 842 unigenes have significant homology to known Arabidopsis gene loci. The high proportion of significantly changing genes reflects our generation and use of a microarray targeted to dormancy processes in Ribes. Unigenes considered to have significantly changing gene expression profiles following ANOVA analysis were subjected to K-means clustering to group similarly regulated genes together (Figure 2). For simplicity, average gene expression profiles are shown for each group to indicate general trends. The four groups obtained (Set 1-Set 4, Figure 2) consist of 260 (227 annotated), 181 (160), 124 (107) and 475 (348) unigenes respectively.
Validation of microarray gene expression profiling
Genetic mapping and association of candidate genes with dormancy traits
Quantitative trait loci interval mapping previously identified a QTL for budbreak in 2002 on linkage group 3 , a region which also contains SERB09L03 (GenBank accession GT022654) and SERB03C05 (GT027101). The recombination frequency between these SNP-markers was zero. The maximum LOD (logarithm (base 10) of odds) score in this region was 3.4. In 2003, the LOD score for budbreak in this region was less than the usual threshold of 3.0, but the QTL effect was the same sign . A more powerful QTL by environment (year) analysis of the combined 2003 and 2004 data  showed that SERB09L03 (GT022654) and SERB03C05 (GT027101) were significantly associated with budbreak (p-value < 0.001) and that there was no significant QTL by year interaction. The mean budbreak for the heterozygous genotype at these SNPs was 4.0 days (se 2.35) earlier than the homozygote in 2002 and 1.6 (se 0.84) earlier in 2003. These two ESTs encode proteins with homology to an acetyl-CoA carboxylase (blastx 2e-47) and a calmodulin-binding protein (blastx 2e-17), respectively. Another EST, SERB02K12 (GT026930) on linkage group 4, was found using the QTL by environment analysis to be associated with budbreak (p-value < 0.001). The mean budbreak for the heterozygous genotype at this SNP was 2.5 days (se 1.38) later than the homozygote in 2002 and 2.1 (se 1.44) later in 2003. This has derived amino acid homology to a beta-tubulin (blastx 1e-94). Specific gene expression profiles, including quantitative RT-PCR validation, is presented in Figure 3 and indicates that the genes encoding acetyl-CoA carboxylase and beta-tubulin have similar patterns, with relatively constant expression to week 8, followed by rapid accumulation of transcripts later in the time course. The gene encoding the calmodulin-binding protein however, initially dips in activity at week 2, recovers and gradually declines through the remainder of the sampling period.
The need for blackcurrant varieties requiring lower levels of chilling during the dormant period is likely to increase in the future , although current unpredictability of weather events in many of the growing regions means that utilising a range of varieties with varying chilling requirements is probably advisable for most growers. For example, varieties that break bud too early may be vulnerable to damage by spring frosts during the flowering period in some northern areas, notably Scandinavia . However, the selection of varieties based on budbreak characteristics currently depends largely on phenotypic assessments in the field using prevailing weather conditions. To move beyond this into a more focused selection procedure requires detailed knowledge about the genetic control of budbreak and subsequent developmental events in blackcurrant, which will enable development of molecular markers directly linked to traits of interest.
General observations of gene activity in blackcurrant bud tissue
This study provides insight into the genetic components associated with dormancy transition in blackcurrant. A targeted blackcurrant cDNA microarray was constructed and utilised to identify sets of genes with significant differential expression profiles during dormancy release and budbreak. Utilising PCA and identifying differentially expressed probes between sampling stages clearly indicated maximum gene activity between week 8 and week 12 (Figure 1), during the month of February. This was also demonstrated in all four of the primary gene expression profiles determined by K-means clustering (Figure 2), with rapid induction or repression after week 8. Up to this point, the cumulative chilling (temperatures below 7°C) was 1,723 h (see Methods), similar levels to which have been demonstrated under controlled environments in Rubus to develop 50-100% budburst in whole canes or isolated nodes respectively . This was slightly lower than that estimated by Fraser  using cuttings, rather than whole plants, which indicated the chilling requirements of commercial UK cultivars of blackcurrant, including 'Ben Gairn' and 'Ben Hope', to be > 2000 h. In Vitis riparia, almost 100% budburst was achieved in canes after 1,500 h of chilling . Accumulated hours of chilling in our study would indicate that our material received sufficient chilling for budbreak between week 8 and week 12 and, combined with the general observations in gene expression profiles, suggest that the blackcurrant samples include bud material in transition from the endodormant to ecodormant phases .
Putative roles of the significantly changing genes during dormancy release from this study remain speculative and care must be taken in interpreting the data, however a number of gene homologues identified have been previously described in several related studies. Some of these common genes of interest are briefly discussed in the context of dormancy transition and budbreak.
Association of dormancy transition with signalling, transcription factors and secondary metabolism
Almost half of the significant differentially expressed blackcurrant genes clustered into expression profile Set 4 (Figure 2), whereby gene activity in general was highest in early winter, subsequently declining throughout the remaining sampling period into spring-time. This set was overrepresented by GO terms associated with transcription factors (TFs; Additional File 1). Of note within this group are genes encoding TFs (see Additional File 4), of which there are several zinc finger family members (B-box, CCCH, C2H2, and C3HC4/RING types) which are known to play important roles in signal transduction pathways of plants . A zinc finger protein from Quercus  has been implicated in budburst, which is related to a TF involved in seed germination from Arabidopsis . Another class of TF which shows gradual down-regulation through the time course in Ribes encodes a Short Vegetative Phase (SVP) type MADS box transcription factor (Additional File 4). Homologues of SVP have been identified in Japanese apricot  and raspberry , where the gene has similar expression patterns to that observed in blackcurrant under controlled conditions leading to dormancy release. SVP has been shown in Arabidopsis to act as a direct repressor of flowering  and appears to have elevated expression under short-day environments . A TF family highly related to SVP, termed Dormancy Associated MADS-box (DAM), has been described in several other perennial species, including poplar , leafy spurge [32, 30] and peach . Expression of two leafy spurge genes encoding DAM proteins reflect the dormancy status within the leaf meristem, with both genes induced by cold temperatures and differentially expressed according to day length. MADS-box TFs are clearly important components in signalling and dormancy transition in several species.
Oxidative stress has been reported to be an important factor in controlling budbreak, and induction of dormancy release by artificial stimuli, such as heat stress and hydrogen cyanamide, leads to the differential expression of several genes associated with this process in grape . From the blackcurrant array data, there is evidence of transient expression of several genes whose products are coupled with this mechanism (Additional File 4), including ascorbate peroxidase (Set 1), calcium-dependent protein kinases (Set 1 & Set 4), glutathione S-transferases (GST; Set 1 & Set 4), catalase (Set 4) and sucrose synthase (Set 4). Similar mechanisms appear to also be involved in dormancy release of raspberry leaf buds, with transient induction of genes encoding ascorbate peroxidase, catalase and GST . Other signalling systems, demonstrated by artificial dormancy control in grape, involve ethylene and abscisic acid (ABA), leading to the hypothesis that there is complex interaction between these mechanisms to control dormancy release . Within the blackcurrant data, there is evidence of ethylene signalling during transition with transient expression of several ethylene response factors (ERFs; Set 4), which have also been observed in leafy spurge . Many of the blackcurrant genes strongly induced later in the time-course (Set 3, Figure 2) are associated with secondary metabolism, specifically those leading to flavonoid biosynthesis. These include genes encoding chalcone isomerase, chalcone synthase, dihydroflavonol reductase and flavonol synthase (Additional File 4), some of which have also been reported to be significantly regulated during dormancy release in raspberry  and leafy spurge [14, 37]. Flavonoids provide protection of tissues against UV damage (reviewed in ) and are known to control auxin transport in plants, subsequently altering growth of tissues such as axillary buds .
Candidate genes associated with dormancy traits
Genetic mapping of significantly expressed ESTs enabled identification of three genes which co-localise with previously characterised blackcurrant budbreak QTL , two of which map with high significance to a single locus on linkage group 3 and another which maps to linkage group 4 (Figure 4). These associations provide independent evidence that these genes may be implicated in dormancy processes, however it should be noted that they may also be simply co-located with the causal genes. Each of these ESTs is discussed on the basis of their homologies to known genes and gene products, and in the context of previous work on dormancy in other species.
Acetyl CoA carboxylase (ACCase, encoded by GT022654) functions to provide malonyl-CoA as a substrate for biosynthesis of fatty acids  and secondary metabolites, such as flavonoids and anthocyanins  which are essential for plant development and protection. It is known that the homomeric form is induced by UV-B radiation for flavonoid production . The gene expression profile derived from this EST (Figure 3) is a member of Set 3 of the clustered groups (Figure 2), which has overrepresentation of terms associated with secondary metabolism and phenylpropanoid biosynthesis. ACCase is associated with rapidly dividing tissues and early stages of cell growth and development and this gene is clearly induced in expression during early spring. In addition, an ACCase has been shown to be induced in floral buds of common bean .
Clone GT027101 encodes a gene product with homology to calmodulin-binding protein. Calmodulin itself binds to calcium, which is common regulator of many different protein targets and a transducer of secondary messenger signals . As previously mentioned, oxidative stress is also known to induce calcium signalling in plant cells . Calmodulin-binding proteins exist in many different forms with a range of associated functions and have been demonstrated to control calcium signalling in woody perennial species such as grape  and blueberry . The expression profile of the gene in blackcurrant (Figure 3) indicates repression from week 4 and is present within cluster Set 4 (Figure 2), which contains other genes with GO terms associated with transcription factors. The polypeptide encoded by this gene also shows similarity to a putative pheromone receptor protein AR781 in Arabidopsis, whose EST profiles show high levels of expression in bud and seed tissue (data from dbEST ; Id 5922760).
A third blackcurrant gene (GT026930), which maps to a budbreak QTL on linkage group 4, encodes a protein with high homology to beta-tubulin, which is an essential component of microtubules within the cytoskeleton. High levels of beta-tubulin expression is a clear indicator of rapidly dividing and expanding cells and the gene expression profile shows rapid induction after week 8 (Figure 2). This profile is a member of Set 1 (Figure 3) which contains many terms representative of increased protein metabolism and biosynthesis. From previous studies, beta-tubulin is consistently induced on dormancy release in buds of a range of woody perennial species, including Malus and Rosa , and is proposed as a genetic marker for dormancy status in such species. It has also been demonstrated that a beta-tubulin gene is induced in expression during bud-burst in poplar .
Extending this study to look for association of these genes with dormancy-related traits in a wide blackcurrant germplasm collection is planned and will further substantiate their roles in physiological processes related to budbreak. However, it is clear that they have excellent potential as genetic markers to assist future breeding strategies. Further work to confirm the genetic map locations for both phenological and fruit quality traits is currently in progress with an extended blackcurrant mapping population and appropriate phenotyping.
This study provides insight into the genetic control of dormancy transition in blackcurrant, and the use of microarrays over a time course leading up to budbreak shows key changes in gene expression at budbreak. Further mapping of ESTs with significant differential changes in expression around budbreak enabled the identification of three genes which co-localise with previously-characterised blackcurrant QTL, and it is concluded that these genes have probable roles in the dormancy and budbreak processes, and can therefore provide a basis for the development of genetic markers for future breeding deployment.
Blackcurrant (Ribes nigrum L. cv. Ben Hope) plants were grown in the field at SCRI, Dundee, Scotland (56° 27'N, 3° 04'W) to four years of age. Leaf buds were sampled at eight time-points: week 0, 7/12/05; week 2, 21/12/05; week 4, 04/01/06; week 8, 01/02/06; week 12, 01/03/06; week 14, 16/03/06; week 16, 29/03/06; week 43, 16/10/06. Air temperature measurements were taken from October and cumulative hours of chilling (below 4°C) determined: week 0, 610 h; week 2, 837 h; week 4, 1,140 h; week 8, 1,723 h; week 12, 2,262 h; week 14, 2,622 h; week 16, 2,847 h. Buds were taken from three biological replicate pools at random positions within the bush, flash frozen in liquid nitrogen and stored at -80°C. Candidate genes were mapped using the reference blackcurrant population designated 9328, described by Brennan et al .
Total RNA was extracted from frozen blackcurrant bud material (100 mg) using the Plant RNeasy Mini Extraction Kit (Qiagen) according to the manufacturer's recommendations (substituting buffer RLC for RLT and including 10% v/v RNA Isolation Aid (Ambion) and 1% v/v β-mercaptoethanol). RNA quality was checked by spectrophotometry and integrity assessed using a Bioanalyzer (Agilent Technologies).
cDNA library construction and EST sequencing
Only very limited numbers of publicly accessible gene sequences are available for blackcurrant (less than 30 sequences present in NCBI , 2009) and closely related woody perennial species. Therefore a cDNA library was generated from blackcurrant leaf bud tissue, which was used as a source of probes for construction of a custom microarray. Utilisation of these expressed gene sequences maximised the chances of identifying putative candidate genes associated with dormancy traits and budbreak in blackcurrant. Total RNA from blackcurrant buds (equal quantities of each) were pooled from five time-points (week 0, week 2, week 4, week 8 and week 12) prior to and potentially including budbreak. Messenger RNA was purified from 100 μg total RNA using Dynabeads mRNA Purification Kit (Invitrogen) as recommended. A directional cDNA library was constructed in pSPORT1 from pooled mRNA using the Superscript Plasmid System (Invitrogen) according to the manufacturer's instructions. Randomly selected clones (~9,000) were picked and stored as glycerol stocks. From these, plasmids were prepared using Multiscreen 96-well Filter Plates (Millipore) as recommended. Sequencing reactions were carried out using M13 reverse primer and BigDye version 3.1 Kit (Applied Biosystems) according to the manufacturer's instructions. Completed sequencing reactions were purified using ethanol precipitation prior to analysis on an ABI PRISM 3730 DNA Sequencer (Applied Biosystems). Sequences were trimmed according to chromatogram quality criteria, using the base-calling program Phred  at a score >20, and vector-derived sequences were removed using Crossmatch software .
In order to estimate the redundancy of the ESTs, related cDNA sequences were assembled as contiguous sequences using the CAP3 program  with default settings of all parameters . ESTs were defined as redundant when they exhibited more than 95% identity over aligned regions and were assembled into a single contiguous sequence.
Similarities to previously identified sequences were obtained by searching public databases using the BLASTn and BLASTx algorithms  for nucleotide and deduced amino acid sequences respectively. Local databases containing non-redundant nucleotide and protein sequences obtained from the National Center for Biotechnology Information  were searched. Matches were considered non-significant (no match) when e-values were greater than 0.01. All EST sequences used to fabricate microarrays have been submitted to dbEST at GenBank .
Blackcurrant cDNA inserts were amplified by PCR using plasmid DNA template and M13 forward and reverse primers that span the multiple cloning site of the vector. Each reaction was performed in 100 μl containing 50 ng of plasmid, M13 primers (0.5 μM each), 0.2 mM dNTPs, 2 mM dNTPs, 2 mM MgCl2 and 0.10 U/μl Taq DNA Polymerase (Promega) in 1 × PCR buffer. PCR conditions were 94°C for 3 min for 1 cycle and then 94°C for 30 s, 54°C for 30 s, 72°C for 2 min for 38 cycles, 72°C for 7 min for 1 cycle. PCR products were purified using the MinElute 96 UF PCR Purification Kit (Qiagen) following the manufacturer's recommendations. Purified PCR products were standardised to 400 μg/μl.
A total of 3,633 cDNA amplicons from blackcurrant library clones were prepared for spotting onto modified glass slides (Nexterion Slide A, Schott Glass). Purified amplified cDNA (7.5 μl at 400 μg/μl) was mixed with 2.5 μl of 99% dimethylsulfoxide (DMSO), and transferred to 384-well print plates (Genetix). Probes were printed in triplicate with a defined random pattern onto each slide using a 24-pin robotic system (Q-Array Mini, Genetix) at constant relative humidity (50%) at 23°C. Arrays were allowed to air-dry for 10 min and the cDNA spots were immobilised following the Nexterion protocol. The quality of the printing process was checked using the SpotCheck Microarray Slide QC kit (Genetix). Details of the blackcurrant dormancy array design can be found at ArrayExpress (accession A-MEXP-1694; )
Target total RNA (2 μg) was labelled using the 3DNA Array 900 Kit (Genisphere) as recommended. RNAs from consecutive time-points were hybridised on the same arrays in a two-colour loop-type design (see ArrayExpress accession E-TABM-807 ) incorporating three biological replicates for each time point (8 time points, 24 samples, 12 arrays total).
Arrays were scanned using an ArrayWoRx Auto scanner (Applied Precision) at appropriate exposure settings for Cy3 (595 nm) and Cy5 (685 nm) at 9.75 μm resolution, generating separate TIFF images. Exposure levels were adjusted to compensate for slight variations in labelling efficiencies. Data were acquired from images using GenePix Pro software (Molecular Devices) and median signal and background intensities were determined for the Cy3 and Cy5 channels for each spot on each microarray. Background-subtracted intensity values were imported into GeneSpring (v.7.3; Agilent Technologies), whereby data from replicate spots within each array were averaged.
Microarray data analysis
Data sets for each array were normalised in GeneSpring using the LOWESS (LOcally WEighted polynomial regreSSion) algorithm to minimize differences in dye incorporation efficiency , prior to re-importing the data as single colour intensity values. Subsequent analysis was carried out using only the first seven sampling points (week 0, week 2, week 4, week 8, week 12, week 14 and week 16) representing winter-spring seasons. Scaling of the data was performed using default GeneSpring settings. Data were filtered by expression level (>100 in 3/21 samples) to remove unreliable, low intensity data points, leaving 2,320 genes for downstream analysis. Data from biological replicates were combined in GeneSpring for statistical analysis: (i) volcano plots (p-value <0.05, fold-change >2 ×) were used to identify differentially expressed genes between sample time points; (ii) 1-way ANOVA was used to estimate differential response between time points (default parameters: parametric test, variances not assumed equal (Welch ANOVA)) and probes were selected (1,040 genes) that were judged to be statistically significant (p-value <0.01) with Benjamini & Hochberg multiple testing correction .
Clustering of gene expression profiles across the developmental time series was performed in GeneSpring using the K-means algorithm. Default parameters (100 iterations, Pearson measure as similarity correlation) were used to generate 4 cluster sets from the selected ANOVA gene-list.
Analysis of Gene Ontology terms
Blackcurrant unigenes used to construct the microarray were annotated with the most significant BLAST homologue to Arabidopsis gene loci. These homologues from selected gene-lists were subsequently analysed for enrichment of Gene Ontology (GO) terms using the Term Enrichment Tool  at AmiGO (version 1.7; ). The whole microarray gene set was used as the defined background and thresholds of p-value 0.05 and a minimum of 1 gene product were used.
Real time RT-PCR validation
Contiguous EST sequences containing probes selected following microarray analysis were used to design Universal Probe Library (UPL, Roche) assays using default parameters . Primer and probe sequences are detailed in Additional File 3 and were designed to the three candidate genes described and also fifteen additional randomly selected unigenes with contrasting expression patterns. Total RNA (5 μg) was reverse transcribed using oligo d(T)18 and You-Prime RT Beads (GE Healthcare) as recommended. cDNA was diluted to 50 μl and 2 μl of each sample used as template for UPL assays with FastStart TaqMan Probe Master ROX (Roche) as recommended, with standard cycling and data acquisition on 7500 Fast Start (Applied Biosystems) equipment. Assay efficiencies were evaluated with standard curves for reference (eIF4A) and each test gene, which determined if ΔΔCt or Pfaffl  calculations were used to estimate relative gene expression profiles.
Candidate gene mapping
Inserts from 112 selected EST clones (Additional File 5) identified by ANOVA filtering were completely sequenced using Sanger technology as above, independently with both M13 forward and reverse priming. Sequences were assembled using Sequencher software (v. 4.9, Genecodes) with default parameters. Oligonucleotide primers (Additional File 5) were designed using the Primer 3 programme . If the lengths of sequenced clones were greater than 1,000 bp, two sets of primers were designed. Genomic DNA was isolated using DNeasy Plant Mini Kit (Qiagen) as recommended, from the seed parent (S36-1-100), pollen parent (S10-2-27/28) and six segregating progeny. PCR was used to amplify gene fragments from each of the lines in a 20 μl reaction volume containing 10-50 ng genomic DNA, 1.0 U of Roche Taq Polymerase, 1 × reaction buffer containing 1.5 mM MgCl2 and supplemented with 0.2 mM dNTPs and 1.0 μM each primer. Thermocycling consisted of 5 min at 94°C; 7 cycles of 30 s at 94°C, 30 s at 65°C, and 30 s at 72°C decreasing to 58°C at 1°C per cycle, followed by 25 cycles of 30 s at 94°C, 30 s at 58°C and 30 s at 72°C, followed by 7 min at 72°C. Single PCR products were gel purified from each reaction using the Qiagen MinElute Gel Extraction Kit according to the manufacturer's specifications and were sequenced using the forward primer used in PCR amplification as above. SNPs were identified by visual inspection in Sequencher v4.6 (Gene Codes Corporation) and Pyrosequencing assays were designed using the Pyrosequencing Assay Design Software v1 (Qiagen).
SNPs which were identified by direct sequencing as being heterozygous in one or both parents were subsequently mapped in the segregating mapping population using the Pyrosequencing SNP detection platform. Forward and reverse primers, and internal sequencing primers within a few bases of the identified SNP, were designed using the Pyrosequencing Primer Design Software. One of the external primers was biotinylated and 25 μl standard PCR reactions were performed as recommended, with thermocycling consisting of 5 min at 94°C, 40 cycles of 30 s at 94°C, 30 s at 58°C and 30 s at 72°C, followed by 10 min at 72°C. Biotinylated products were immobilised onto Sepharose beads and the samples were analysed using a PSQ 96MA system together with SNP reagent kits and SNP software (Qiagen). Two mapping parents and 93 segregating progeny were scored and mapped using the methodology developed by Brennan et al. .
Simple sequence repeat identification and analysis
SSRs were identified using the Sputnik programme and primers were designed using Primer 3 software . Each 10 μl PCR reaction contained 20 ng genomic DNA, 1 μM each forward and reverse primers (one of which was labelled with 6-FAM), 0.2 mM dNTPs, 1.0 U of Roche Taq Polymerase, in 1 × reaction buffer containing 1.5 mM MgCl2. Thermocycling conditions were 5 min at 94°C, 7 cycles of 30 s at 94°C, 30 s at 65°C, and 30 s at 72°C decreasing to 58°C at 1°C per cycle, followed by 25 cycles of 30 s at 94°C, 30 s at 58°C and 30 s at 72°C, followed by 7 min at 72°C. Following PCR, fragments (along with 500 ROX-labelled size standard) were separated on an ABI Prism 3730 DNA Analyzer according to the manufacturer's instructions (Applied Biosystems). Genotypes were identified using GeneMapper v3.7 software (Applied Biosystems).
Genetic mapping and QTL analysis
As discussed in Brennan et al. , the mapping population S36-1-100 × S10-2-27/28 contains both sibs from these parents and selfs from SCRI S36-1-100. The new markers were added to the existing linkage map using JoinMap 3  to combine data from the selfs and sibs, as described . QTL analysis of the budbreak data was carried out with the new map and the sibs from this cross using MapQTL 5 , as described in Brennan et al. . The new markers were also tested using a more powerful mixed model analysis that combines the sib and self phenotypic data in a QTL by environment (year) analysis of the budbreak data from 2002 to 2003, as described .
The authors are very grateful to Jim McNicol (BioSS, Dundee) for providing advice on microarray experimental design and GO analysis, and Mark Young (SCRI) for supplying meteorological data. Thanks also to Dr Louise Shepherd & Dr Mary Woodhead for their constructive comments on the manuscript. This work was funded under the Scottish Government RERAD Programme.
- Jensen S: Judgement on the harvest of the year 2009. Proceedings of the 15th European Blackcurrant Conference, Nyborg, Denmark. 2009, [http://www.gartneriraadgivningen.dk/Nyheder/Solbaer/Presentations/Presentations.html]Google Scholar
- Brennan RM: Currants and gooseberries. Fruit Breeding Volume II:Vine and Small Fruits Crops. Edited by: Janick J, Moore JN. New York, J Wiley & Sons Inc, 1996, 191-295.Google Scholar
- Brennan RM, Graham J: Improving fruit quality in Rubus and Ribesthrough breeding. Functional Plant Science and Biotechnology. 2009,Google Scholar
- Martin GC: Bud dormancy in deciduous fruit trees. Plant physiology: a treatise. Edited by: Steward FC. New York Academic Press Inc, 1991, 183-225.Google Scholar
- Atkinson CJ, Sunley RJ, Jones HG, Brennan R, Darby P: Winter chill in fruit. UK Department of Food, Environment and Rural Affairs Report NoCTC026. 2004Google Scholar
- Fuchigami LH, Wisniewski ME: Quantifying bud dormancy: physiological approaches. HortScience. 1997, 32: 618-623.Google Scholar
- Jones HG, Brennan RM: Potential impacts of climate change on soft fruit production: the example of winter chill in Ribes. Acta Horticulturae. 2009, 838: 28-32.Google Scholar
- Faust M, Liu D, Wang SY, Stutte GW: Involvement of apical dominance in winter dormancy of apple buds. Acta Horticulturae. 1995, 395: 47-56.View ArticleGoogle Scholar
- Rohde A, Howe GT, Olsen JE, Moritz T, Van Montagu M, Juntilla O, Boerjan W: Molecular aspects of bud dormancy in trees. MolecularBiology of Woody Plants. Edited by: Jain S, Minocha SC. Dordrecht, Netherlands, Kluwer Academic, 2000, 1: 89-134.View ArticleGoogle Scholar
- Snelling C, Langford G: The development of low chill blackcurrants from the New Zealand breeding programme. Acta Horticulturae. 2002, 585: 167-169.View ArticleGoogle Scholar
- UK Climate Impact Projections. 2009, [http://www.ukcip.org.uk]
- Rigby JR, Porporato A: Spring frost risk in a changing climate. Geophysical Research Letters. 2008, 35: L12703-10.1029/2008GL033955.Google Scholar
- Lang GA: Dormancy: a new universal terminology. HortScience. 1987, 22: 817-820.Google Scholar
- Horvath DP, Anderson JV, Chao WS, Foley ME: Knowing when to grow: signals regulating bud dormancy. Trends in Plant Science. 2003, 8: 534-540. 10.1016/j.tplants.2003.09.013.PubMedView ArticleGoogle Scholar
- Brennan R, Jorgensen L, Hackett C, Woodhead M, Gordon SL, Russell J: The development of a genetic linkage map of blackcurrant (Ribes nigrum L.) and the identification of regions associated with key fruit quality and agronomic traits. Euphytica. 2008, 161: 19-34. 10.1007/s10681-007-9412-8.View ArticleGoogle Scholar
- Pacey-Miller T, Scott K, Ablett E, Tingey S, Ching A, Henry R: Genes associated with the end of dormancy in grapes. Functional Integrative Genomics. 2003, 3: 144-152. 10.1007/s10142-003-0094-6.PubMedView ArticleGoogle Scholar
- Tattershall EAR, Grimplet J, DeLuc L, Wheatley MD, Vincent D, Osborne C, Ergül A, Lomen E, Blank RR, Schlauch KA, Cushman JC, Cramer GR: Transcript abundance profiles reveal larger and more complex responses of grapevine to chilling compared to osmotic and salinity stress. Functional Integrative Genomics. 2007, 7: 317-333. 10.1007/s10142-007-0051-x.View ArticleGoogle Scholar
- Mathiason K, He D, Grimplet J, Venkateswari J, Galbraith DW, Or E, Fennell A: Transcript profiling in Vitis riparia during chilling requirement fulfilment reveals coordination of gene expression patterns with optimized bud break. Functional Integrative Genomics. 2009, 9: 81-96. 10.1007/s10142-008-0090-y.PubMedView ArticleGoogle Scholar
- Derory J, Léger P, Garcia V, Schaeffer J, Hauser M-T, Salin F, Luschnig C, Plomion C, Glössl J, Kremer A: Transcriptome analysis of bud burst in sessile oak (Quercus petraea). New Phytologist. 2006, 170: 723-728. 10.1111/j.1469-8137.2006.01721.x.PubMedView ArticleGoogle Scholar
- Mazzitelli L, Hancock RD, Haupt S, Walker PG, Pont SDA, McNicol J, Cardle L, Morris J, Viola R, Brennan R, Hedley PE, Taylor MA: Coordinated gene expression during phases of dormancy release in raspberry (Rubus idaeus L.) buds. Journal Experimental Botany. 2007, 58: 1035-1045. 10.1093/jxb/erl266.View ArticleGoogle Scholar
- ArrayExpress. [http://www.ebi.ac.uk/microarray-as/ae/]
- AmiGO. [http://amigo.geneontology.org/cgi-bin/amigo/term_enrichment]
- Hackett CA, Russell JR, Jorgenson L, Gordon S, Brennan RM: Multi-environment QTL mapping in blackcurrant (Ribes nigrumL.) using mixed models. Theoretical and Applied Genetics. 2010,Google Scholar
- Lantin B: The chilling requirements of buds of blackcurrant (Ribes nigrum L.) and of some redcurrants (Ribes spp.). Ann Amél Plantes. 1973, 23: 27-44.Google Scholar
- Fraser GA: An analysis of dormancy and chilling of Ribes nigrum L. PhD thesis. University of Reading, UK: 2005,Google Scholar
- Ciftci-Yilmaz S, Mittler R: The zinc finger network of plants. Cellular and Molecular Life Sciences. 2008, 65: 1150-1160. 10.1007/s00018-007-7473-4.PubMedView ArticleGoogle Scholar
- Gualberti G, Papi M, Bellucci L, Ricci I, Bouchez D, Camilleri C, Costantino P, Vittorioso P: Mutations in the Dof zinc finger genes DAG2 and DAG1 influence with opposite effects the germination of Arabidopsis seeds. Plant Cell. 2002, 14: 1253-63. 10.1105/tpc.010491.PubMedPubMed CentralView ArticleGoogle Scholar
- Yamane H, Kashiwa Y, Ooka T, Tao R, Yonemori K: Suppression subtractive hybridization and differential screening reveals endodormancy-associated expression of an SVP/AGL24-type MADS-box gene in lateral vegetative buds of Japanese apricot. Journal American Society Horticultural Science. 2008, 33: 708-716.Google Scholar
- Yant L, Mathieu J, Schmid M: Just say no: floral repressors help Arabidopsis bide the time. Current Opinions Plant Biology. 2009, 12: 580-586. 10.1016/j.pbi.2009.07.006.View ArticleGoogle Scholar
- Horvath DP, Sung S, Kim D, Chao W, Anderson J: Characterization, expression and function of DORMANCY ASSOCIATED MADS-BOX genes from leafy spurge. Plant Molecular Biology. 2010,Google Scholar
- Ruttink T, Arend M, Morreel K, Storme V, Rombauts S, Fromm J, Bhalerao RP, Boerjan W, Rohdea A: A molecular timetable for apical bud formation and dormancy induction in poplar. Plant Cell. 2007, 19: 2370-2390. 10.1105/tpc.107.052811.PubMedPubMed CentralView ArticleGoogle Scholar
- Horvath DP, Chao WS, Suttle JC, Thimmapuram J, Anderson JV: Transcriptome analysis identifies novel responses and potential regulatory genes involved in seasonal dormancy transitions of leafy spurge (Euphorbia esula L.). BMC Genomics. 2008, 9: 536-554. 10.1186/1471-2164-9-536.PubMedPubMed CentralView ArticleGoogle Scholar
- Jiminez S, Reighard GL, Bielenberg DG: Gene expression of DAM5 and DAM6 is suppressed by chilling temperatures and inversely correlated with bud break rate. Plant Molecular Biology. 2010, 73: 157-67. 10.1007/s11103-010-9608-5.View ArticleGoogle Scholar
- Halaly T, Pang X, Batikoff T, Crane O, Keren A, Venkateswari J, Ogrodovitch A, Sadka A, Lavee S, Or E: Similar mechanisms might be triggered by alternative external stimuli that induce dormancy release in grape buds. Planta. 2008, 228: 79-88. 10.1007/s00425-008-0720-6.PubMedView ArticleGoogle Scholar
- Ophir R, Pang X, Halaly T, Venkateswari J, Lavee S, Galbraith D, Or E: Gene-expression profiling of grape bud response to two alternative dormancy-release stimuli expose possible links between impaired mitochondrial activity, hypoxia, ethylene-ABA interplay and cell enlargement. Plant Molecular Biology. 2009, 71: 403-423. 10.1007/s11103-009-9531-9.PubMedView ArticleGoogle Scholar
- Dogramaci M, Horvath DP, Chao WS, Foley ME, Christoffers MJ, Anderson JV: Low temperatures impact dormancy status, flowering competence, and transcript profiles in crown buds of leafy spurge. Plant Molecular Biology. 2010, 73: 207-226. 10.1007/s11103-010-9621-8.PubMedView ArticleGoogle Scholar
- Horvath DP, Anderson JV, Soto-Suarez M, Chao WS: Transcriptome analysis of leafy spurge (Euphorbia esula) crown buds during shifts in well-defined phases of dormancy. Weed Science. 2006, 54: 821-827. 10.1614/WS-06-06R1.1.View ArticleGoogle Scholar
- Ferrer J-L, Austin MB, Stewart C, Noel JP: Structure and function of enzymes involved in the biosynthesis of phenylpropanoids. Plant Physiology Biochemistry. 2008, 46: 356-370. 10.1016/j.plaphy.2007.12.009.PubMedPubMed CentralView ArticleGoogle Scholar
- Lazar G, Goodman HM: MAX1, a regulator of the flavonoid pathway, controls vegetative axillary bud outgrowth in Arabidopsis. Proceedings National Academy Sciences. 2006, 103: 472-476. 10.1073/pnas.0509463102.View ArticleGoogle Scholar
- Harwood JL: Fatty acid metabolism. Annual Review Plant Physiology Plant Molecular Biology. 1988, 39: 101-138. 10.1146/annurev.pp.39.060188.000533.View ArticleGoogle Scholar
- Sasaki Y, Nagano Y: Plant acetyl-CoA carboxylase: structure, biosynthesis, regulation and gene manipulation for plant breeding. Bioscience Biotechnology Biochemistry. 2004, 68: 1175-1184. 10.1271/bbb.68.1175.View ArticleGoogle Scholar
- Figueroa-Balderas RE, Garcia-Ponce B, Rocha-Sosa M: Hormonal and stress induction of the gene encoding common bean acetyl-coenzyme A carboxylase. Plant Physiology. 2006, 142: 609-619. 10.1104/pp.106.085597.PubMedPubMed CentralView ArticleGoogle Scholar
- Zielinski RE: Calmodulin and calmodulin-binding proteins in plants. Annual Review Plant Physiology. 1998, 49: 697-725. 10.1146/annurev.arplant.49.1.697.View ArticleGoogle Scholar
- Pang X, Halaly T, Crane O, Keilin T, Keren-Keiserman A, Ogrodovitch A, Galbraith D, Or E: Involvement of calcium signalling in dormancy release of grape buds. Journal Experimental Botany. 2007, 58: 3249-3262. 10.1093/jxb/erm172.View ArticleGoogle Scholar
- Dhanaraj AL, Alkharouf NW, Beard HS, Chouikha IB, Matthews BF, Wei H, Arora R, Rowland LJ: Major differences observed in transcript profiles of blueberry during cold acclimation under field and cold room conditions. Planta. 2007, 225: 735-751. 10.1007/s00425-006-0382-1.PubMedView ArticleGoogle Scholar
- National Center for Biotechnology Information Expressed Sequence Tags database. [http://www.ncbi.nlm.nih.gov/dbEST/]
- Bergervoet JHW, Jing HC, van den Hout JWE, de Castro RD, Kunneman BPAM, Bino RJ, Groot SPC: Expression of beta-tubulin during dormancy induction and release in apical and axillary buds of five woody species. Physiologia Plantarum. 1999, 106: 238-245. 10.1034/j.1399-3054.1999.106214.x.View ArticleGoogle Scholar
- Rohde A, Ruttink T, Hostyn V, Sterck L, Van Driessche K, Boerjan W: Gene expression during the induction, maintenance, and release of dormancy in apical buds of poplar. Journal Experimental Botany. 2007, 58: 4047-4060. 10.1093/jxb/erm261.View ArticleGoogle Scholar
- National Center for Biotechnology Information Entrez Nucleotide database. [http://www.ncbi.nlm.nih.gov/nuccore]
- Phrap software. [http://www.phrap.org/]
- Huang X, Madan A: CAP3: A DNA sequence assembly program. Genome Research. 1999, 9: 868-77. 10.1101/gr.9.9.868.PubMedPubMed CentralView ArticleGoogle Scholar
- CAP3 Sequence Assembly Program. [http://seq.cs.iastate.edu/cap3.html]
- Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ: Basic local alignment search tool. Journal Molecular Biology. 1990, 5: 403-10.View ArticleGoogle Scholar
- Yang YH, Dudoit S, Luu P, Lin DM, Peng V, Ngai J, Speed TP: Normalization for cDNA microarray data: a robust composite method addressing single and multiple slide systematic variation. Nucleic Acids Research. 2002, 30: e15-10.1093/nar/30.4.e15.PubMedPubMed CentralView ArticleGoogle Scholar
- Benjamini Y, Hochberg Y: Controlling the false discovery rate: a practical and powerful approach to multiple testing. Journal Royal Statistical Society Series B. 1995, 1: 289-300.Google Scholar
- Boyle EI, Weng S, Gollub J, Jin H, Botstein D, Cherry JM, Sherlock G: GO::TermFinder--open source software for accessing Gene Ontology information and finding significantly enriched Gene Ontology terms associated with a list of genes. Bioinformatics. 2004, 20: 3710-5. 10.1093/bioinformatics/bth456.PubMedPubMed CentralView ArticleGoogle Scholar
- Universal Probe Library Assay Design Center. [https://www.roche-applied-science.com/sis/rtpcr/upl/index.jsp?id=UP030000]
- Tichopad A, Dilger M, Schwarz G, Pfaffl MW: Standardized determination of real-time PCR efficiency from a single reaction set-up. Nucleic Acids Research. 2003, 31: 20-10.1093/nar/gng122. e122View ArticleGoogle Scholar
- Rozen S, Skaletsky HJ: Primer3 on the WWW for general users and for biologist programmers. Bioinformatics Methods and Protocols: Methods in Molecular Biology. Edited by: Krawetz S, Misener S. Humana Press, Totowa, NJ: 2000. 365-386.Google Scholar
- Van Ooijen JW, Voorrips RE: JoinMap® 3.0, Software for the calculation of genetic linkage maps. Plant Research International, Wageningen, Netherlands. 2001Google Scholar
- Van Ooijen JW: MapQTL ® 5, Software for the mapping of quantitative trait loci in experimental populations. Kyazma BV, Wageningen, Netherlands. 2004Google 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.