EcoTILLING in Beta vulgaris reveals polymorphisms in the FLC-like gene BvFL1that are associated with annuality and winter hardiness
© Frerichmann et al.; licensee BioMed Central Ltd. 2013
Received: 12 September 2012
Accepted: 21 February 2013
Published: 25 March 2013
Sugar beet (Beta vulgaris ssp. vulgaris L.) is an important crop for sugar and biomass production in temperate climate regions. Currently sugar beets are sown in spring and harvested in autumn. Autumn-sown sugar beets that are grown for a full year have been regarded as a cropping system to increase the productivity of sugar beet cultivation. However, for the development of these “winter beets” sufficient winter hardiness and a system for bolting control is needed. Both require a thorough understanding of the underlying genetics and its natural variation.
We screened a diversity panel of 268 B. vulgaris accessions for three flowering time genes via EcoTILLING. This panel had been tested in the field for bolting behaviour and winter hardiness. EcoTILLING identified 20 silent SNPs and one non-synonymous SNP within the genes BTC1, BvFL1 and BvFT1, resulting in 55 haplotypes. Further, we detected associations of nucleotide polymorphisms in BvFL1 with bolting before winter as well as winter hardiness.
These data provide the first genetic indication for the function of the FLC homolog BvFL1 in beet. Further, it demonstrates for the first time that EcoTILLING is a powerful method for exploring genetic diversity and allele mining in B. vulgaris.
KeywordsBeta vulgaris EcoTILLING Diversity Flowering Winter hardiness BvFL1
EcoTILLING is a fast and easy method to detect rare SNPs or small indels in target genes in natural populations. It is an adaptation of the TILLING (Targeting Induced Local Lesion In Genomes) technique that is used to detect point mutations in mutant populations . In EcoTILLING, endonucleases such as CEL I are used to cut mismatched sites in the heteroduplex DNA formed by hybridization of different genotypes in a test panel. It is a cost effective technology as sequencing is limited to individual genotypes each representing a different haplotype. EcoTILLING has been used for the characterization of the genetic variability in Arabidopsis thaliana (thale cress) , Musa spp. (various banana species) , Populus trichocarpa (black cottonwood) , Phaseolus vulgaris (common bean) , and Vigna radiata (mung bean) . Furthermore, it has been used for candidate gene-based detection of new alleles conferring resistance to biotic and abiotic stress in Hordeum vulgare (barley) , Oryza sativa (rice) [8, 9], Solanum tuberosum (potato) . Cucumis spp. (including cucumber)  and Solanum lycopersicum (tomato) . EcoTILLING has not been reported for sugar beet (Beta vulgaris ssp. vulgaris L.) which contributes to 22% of the world production of white sugar .
Sugar beets are herbaceous, dicotyledonous plants that belong to the Amaranthaceae family (formerly Chenopodiaceae). The genus Beta is divided into the two sections Corollinae and Beta, the latter of which is further divided into cultivated beets (B. vulgaris ssp. vulgaris), wild sea beets (B. vulgaris ssp. maritima L.) and wild beets (B. vulgaris ssp. adanensis) . Within B. vulgaris ssp. vulgaris, four cultivated groups can be distinguished: fodder beet, leaf beet, garden beet and sugar beet. While leaf beets and garden beets show an annual or biennial life cycle, sugar beets and also fodder beets are biennial plants that stay in the vegetative phase in their first year, forming a storage root with a high sucrose concentration of up to 20%. Both vernalization and long days are required for stem elongation (bolting) and flowering to occur in the second year of growth. Vernalization in sugar beet is achieved by exposure to cold temperatures for ten to 14 weeks.
Currently, sugar beets are cultivated as a spring sown crop in cool temperate climate regions. Seeds are sown in April and the roots are harvested starting in September. The late formation of a closed leaf canopy in spring is regarded as the main factor limiting beet yield . One strategy to overcome this is the production of autumn sown winter beets which develop a closed canopy earlier in spring. However, breeding of autumn sown winter beets requires sufficient winter hardiness to survive the winter and a system for bolting control which allows bolting for seed production but represses bolting after winter during crop production . With key regulators of flowering and bolting in B. vulgaris recently having been identified [17, 18], bolting control may be achieved by genetic modification which on the one hand allows suppression of bolting after winter for cultivation of beets, but on the other hand enables bolting for seed production .
In order to avoid an untimely transition to the extremely cold-sensitive generative phase  before or during winter, and to facilitate the accumulation of sufficient resources for reproduction, winter-annual and biennial plants growing in temperate zones require vernalization for induction of flowering. Cultivated beets are biennials, whereas annual beets without a requirement for vernalization are frequently observed in wild beet populations [20, 21]. The vernalization response in biennial beets is mediated by the FLOWERING LOCUS T (FT) homolog BvFT1, which in contrast to the promotive action of FT in Arabidopsis functions as a repressor of flowering . Similar to FLC, BvFT1 is gradually down-regulated during the prolonged cold of winter . In annual beets, BvFT1 is not expressed even in the absence of vernalization and was shown to be negatively regulated by the pseudo-response regulator gene BOLTING TIME CONTROL 1 (BTC1), formerly referred to as BvBTC1. This gene is located at the bolting locus B and is a major determinant of the annual growth habit in beet. The dominant BTC1 allele promotes bolting in annuals in response to long days, whereas biennials carry a partial-loss-of-function allele which is not able to mediate the promotive effect of long days without prior vernalization . All cultivated (biennial) beet accessions tested were found to carry the same haplotype whereas the vast majority of wild sea beets harbour haplotypes which resemble the BTC1 allele found in annual reference accessions .
Several other genes in beet have been identified on the basis of homology to floral transition genes in Arabidopsis, including the central regulator of vernalization requirement and response in this species, FLOWERING LOCUS C (FLC) [22, 23]. The FLC-LIKE 1 gene BvFL1 is gradually down-regulated during a prolonged exposure to cold under continuous light . Constitutive expression of BvFL1 in an FLC null mutant of Arabidopsis significantly delayed flowering, suggesting at least partial evolutionary conservation of function between FLC homologs in Arabidopsis and beet.
Interestingly, flowering time control genes also seem to affect frost tolerance, which is the most important factor contributing to winter hardiness [25–27]. Plants can further increase their frost tolerance by a gradual adaptation of the metabolism during a hardening process that occurs at non-freezing temperatures below 10°C. In Arabidopsis, frost tolerance is regulated by the C-REPEAT BINDING FACTOR (CBF) transcription factor family, with plants constitutively overexpressing CBF genes showing an increase in frost tolerance  and elevated levels of FLC expression . Deng et al.  reported that FLC plays a dual role in flowering time control and cold stress response. Interestingly, a recent study suggested that the recruitment of a repressive chromatin complex at the FLC locus involves the cold-induced expression of a long non-coding RNA, termed COLDAIR, from intron 1 of FLC.
In the present study, we established EcoTILLING in B. vulgaris to survey a large panel of cultivated and wild beets for allelic variants of candidate genes for regulators of vernalization requirement and/or winter hardiness. This panel had been phenotyped before for variation in the occurrence of bolting before winter (i.e. in the absence of vernalization) and survival rates after winter . As candidate genes we chose (i) BTC1 and (ii) BvFT1, because of their known functions in the regulation of vernalization requirement and response in beet, and (iii) BvFL1, because of the regulatory role of its homolog FLC in both vernalization and cold stress response in Arabidopsis. We found that haplotype variation at the BvFL1 locus was associated with variation in bolting rate before winter and survival rate after winter. These data provide the first genetic indication for the function of the FLC homolog BvFL1 in beet, and are relevant for sugar beet breeding and our understanding of the bolting time control network in Beta.
Phenotyping and model-based analyses of population structure
Amplification of candidate genes
Overview of candidate genes investigated with EcoTILLING
Primer name and sequence (5’-3’)
Length and portion of genomic sequence covered by EcoTILLING
Portion of reading frame covered by EcoTILLING
MADS box/K box
FL1a: 977 bp FL1b: 632 bp
FL1a-fw tcggactttccctataagct FL1a-rv cacgtgaatcgttacagaca FL1b-fw gctgatagtctgtcccttttgtc FL1b-rv tgactccaacaccacgatgca
1,609 bp (19%)
FT1a: 916 bp FT1b: 713 bp
FT1a-fw tggtacgtgtatgaaacagaagctg FT1a-rv catcaactccatatttggggtg FT1b-fw acccatctatacttgtcgatgacc FT1b-rv caatggggaagtggttcacact
1,629 bp (22%)
BTC1: 996 bp
BTC1-fw cagctgtaggatgttatcgtgctgag BTC1-rv agtaggtgataaggacaagacattgc
992 bp (9%)
Genetic diversity, SNP densities and haplotype frequencies
Across the five amplicons (FL1a, FL1b, FT1a, FT1b and BTC1) a total of 21 SNPs were identified among the 268 accessions tested by LI-COR analyses. Eighteen SNPs were located in introns, while the remaining three SNPs were located in exon 4 of BvFT1 and in exon 3 and exon 4 of BvFL1, the latter SNP of which was non-synonymous. The number of polymorphisms varied from gene to gene and had an overall density of 5.3 SNP/kb. The lowest SNP density was found in BTC1 (2.01 SNP/kb), whereas the highest SNP density occurred in FL1b (9.82 SNP/kb) (Additional file 3). The SNP allocations and gene structures are shown in Figure 3. To evaluate the efficiency of EcoTILLING in B. vulgaris we estimated the rate of false negatives by sequencing all amplicons in four selected accessions, which resulted in a false negative rate of 5%. Sequencing of FT1a, FT1b, FL1a and BTC1 did not reveal additional SNPs which had not been already identified by LI-COR analyses. For the amplicon FL1b one additional SNP was identified after sequencing in a single accession. This SNP was not detected for any of the 268 accession by EcoTILLING on the LI-COR.
BvFL1sequence variations are associated with bolting and survival rate
Significant associations of bolting and survival rate with amplicon
This is the first report of EcoTILLING applied to B. vulgaris. We established EcoTILLING based on a panel of 268 accessions representing the wild and domesticated gene pool of B. vulgaris. In this panel we successfully screened the allelic variation in three genes that are candidates for regulators of vernalization requirement and/or winter hardiness. As a result we were able to provide a snapshot of the species-wide diversity within these genes. Further, we identified haplotypes that are associated with bolting rate before winter and with survival rate, which in turn might be useful for improvement of winter hardiness in sugar beets. Our results show that EcoTILLING is a suitable and cost effective method for allele mining in B. vulgaris.
In most EcoTILLING protocols heteroduplexed DNA is digested by purified CEL I endonuclease. Instead of the purified enzyme, Till et al.  and Galeano et al.  used celery juice obtained from salted out or dialyzed crude celery extract for EcoTILLING screens of Arabidopsis thaliana and Phaseolus vulgaris, respectively. We went a step further and used the crude celery extract (CCE) without further processing, and observed the same activity as compared to the commercial CEL I enzyme Surveyor® (data not shown). Also, CCE was very stable and kept its activity for weeks even when stored at 4°C. As using CCE eliminates the need for special enzyme purification steps like chromatography and specialized laboratory equipment, this increases the cost efficiency of EcoTILLING. We were further able to show that once suitable primers are designed, EcoTILLING provides a high throughput method for the analysis of natural nucleotide diversity in B. vulgaris. Also, EcoTILLING is a rather cost effective method. When evaluating LI-COR gels, signals can be grouped according to size and pattern, and only a limited number of samples per group need to be sequenced to break down the detected variation to the nucleotide level. This drastically reduces the sequencing costs, in our case by 1/3. If only SNPs/haplotypes with effect on the phenotype of interest are sequenced, costs can be further reduced. However, it has to be considered that EcoTILLING is prone to false negative detection because some fragment sizes are masked by background “noise”, due to miss-priming, or because of weaker fluorescence toward the top of each lane and increasing fluorescence “noise” toward the bottom [3, 37, 38]. The false negative rate in our case was 5% which is similar to rates reported in human  and banana  EcoTILLING.
The population structure analysis by AFLP markers indicates that the B. vulgaris accessions can be grouped into three groups (k = 3). Under consideration of the phenotypic classification of the panel (gene bank information as well as our field observations on plant habitus), the three groups can be referred to as a sugar beet group, a fodder beet and garden beet group, and a group comprising leaf beets and wild sea beets. This reflects the evolutionary history of Beta and the selection intensity during the past 200 years of beet breeding. A similar structure has also been described by Jung et al.  and by McGrath et al.  after genotyping with completely different marker systems. Both groups report that sugar beets can be clearly distinguished from B. vulgaris ssp. maritima. In our study, a few accessions were classified differently by genotype (according to AFLP analysis) than by phenotype (see Figure 2). Regarding B.vulgaris ssp. maritima this could hint at gene flow from cultivated forms into wild material, either in their natural habitat or during propagation by gene banks. At the same time, classification by phenotype was sometimes ambiguous. For instance ‘Patak’ accessions from India (PI 116809 and PI 121838), although cultivated, showed a plant habitus whose classification into B.vulgaris ssp. maritima seems more reasonable than classification into any of the cultivated forms. Interestingly, both approaches to account for population structure resulted in significant associations of the same amplicons, hinting at the robustness of the results. Nevertheless, to account for genotyping errors as a source for putative misclassification by STRUCTURE, genotypic outliers were removed from the dataset for a further analysis by TASSEL. These outliers were sugar beet, fodder beet and table beet accessions with an estimated portion of non-cultivated beet genome > 50% and B.vulgaris ssp. maritima accessions with an estimated portion of cultivated beet genome > 50% (see Figure 2). In this analysis previous associations were still highly significant (data not shown).
Comparing all five B. vulgaris forms we observed the highest genetic diversity for the investigated genes in B. vulgaris ssp. maritima. This is indicated by a mean NNF of 0.36 compared to 0.19 in leaf beets followed by fodder beets (0.12), sugar beets (0.08) and garden beets (0.07). The same trend is observed when looking at the average NHF and Ht. Our findings are in accordance with Jung et al.  and Fénart et al.  who reported a higher genetic diversity in B. vulgaris ssp. maritima compared to sugar beets. As selection results in a loss of genetic diversity, it is not surprising that the genetic diversity in B. vulgaris ssp. maritima appears to be higher not only compared to sugar beet but also in comparison to all four cultivar groups taken together. Crop evolution is best understood for sugar beet which has been affected by founder effects as it was derived from a single fodder beet population and also by genetic bottlenecks through introgression of a series of traits from a limited number of genetic resources [43–45]. This explains why sugar beet together with red table beet showed the lowest diversity.
In our study the genetic diversity in sugar beet based on Ht ranged from 0.03 to 0.28 for the single amplicons with an average of 0.17. These estimates are likely to be upward biased as we could not distinguish between the occurrence of non-reference nucleotides in the heterozygous or homozygous state. Nevertheless, gene diversity in our study is lower compared to Li et al.  and McGrath et al. . However, in contrast to Li et al. and McGrath et al. we estimated the genetic diversity for nucleotide polymorphisms in three genes that may have been under selective pressure. This is especially the case for BvFL1 where we estimated Ht values of 0.03 and 0.12 for FL1a and FL1b, respectively, and also for BTC1 (Ht = 0.06). This could be the effect of selection for bolting resistance to prevent bolting caused by late frosts after sowing in spring. Comparing sugar beet elite breeding material with sugar beet germplasm, the genetic diversity turned out to have been further decreased by selection for BTC1 and BvFL1 (see also Additional file 5). At the same time BvFT1 showed even more diversity in SBEBM indicating that this gene is obviously not under selective pressure. This is somehow surprising, as BvFT1 was shown to play a key role in bolting suppression under non-inductive conditions . Still, these data have to be interpreted with care, as sample sizes are moderate and the SBEBM material represents only one breeding company.
For BvFL1, we were able to detect an association with bolting. Four haplotypes of this gene (FL1a_H6, FL1b_H5, FL1b_H6 and FL1b_H10) had a significant effect on bolting rate before winter in B. vulgaris ssp. maritima and/or garden beets. Although variation in FLC is known to affect flowering time in A. thaliana, the role of FLC-like genes outside the Brassicaceae is not well understood , and a functional analysis of BvFL1 in B. vulgaris, e.g. through mutational or transgenic approaches, is still lacking. Effects on bolting rate were not observed for all B. vulgaris forms, which in part may be due to the absence of the divergent haplotypes that affect bolting rate in B. vulgaris ssp. maritima or garden beet. The complete absence of these haplotypes in sugar beet may reflect the breeding history of sugar beet, during which breeders strongly selected against bolting before vernalization [48, 49].
As BTC1 is known to be a major factor controlling bolting without prior vernalization in beets , we expected an effect of BTC1 sequence variations on bolting before winter. However, this was not observed here. This may be due in part to an underrepresentation of annual BTC1 alleles in our panel or the fact that the current analysis was limited to a relatively small portion of the coding sequence of BTC1 (15%; Table 1), and did not include the promoter. Although BvFT1 is known to respond to vernalization and is down-regulated by cold temperatures, which in turn enables induction of flowering , we also did not detect an effect of haplotype variation in this gene on bolting rate. While other reasons for this cannot be excluded, as discussed for BvFL1 and BTC1, it is also conceivable that possible phenotypic effects of haplotype variation at BvFT1 are difficult to detect under the environmental conditions present in the current study.
Besides significant effects on bolting rate, plants with haplotypes FL1a_H6 and FL1b_H3 showed a significant impact on survival rate after winter in B. vulgaris ssp. maritima and leaf beets, respectively. A similar effect has been shown before for A. thaliana ecotypes, where a SNP in intron 1 of FLC led to a 1.6-fold increase in winter survival rates in genotypes carrying a functional FRI allele . The authors suggested that survival after winter is associated with time to bolting. Similarly, we found that the survival rate of truly biennial (vernalization requiring) leaf beet accessions with the FL1b_H3 haplotype was higher (by 20%) when compared to the reference haplotype FL1_H0 (p = 0.006, data not shown). Hence, winter hardiness in sugar beet might be improved by introgressing FL1b_H3 from leaf beet. Interestingly, B. vulgaris ssp. maritima accessions with FL1a_H6, which had a lower survival rate than accessions with the reference haplotype, also showed an increased bolting rate. Furthermore, after removal of B. vulgaris ssp. maritima accessions which bolted before winter, the accessions with haplotype FL1a_H6 did not have a significant effect anymore. Therefore, the lower survival rate observed for this haplotype might be a direct physiological effect of bolting before winter as plants in the generative phase are more frost sensitive . However, by using bolting rate as a cofactor in a further TASSEL analysis, we can exclude that increased frost sensitivity in the generative phase is the mere cause for association of BvFL1 with survival rate since this association stayed significant. Interestingly, Seo et al.  reported that transient cold temperatures and overexpression of CBFs lead to elevated FLC expression and delayed flowering, suggesting a possible role of FLC in cold stress response in A. thaliana and, by analogy, a possible explanation for the detected effect of BvFL1 haplotypes on survival rate in B. vulgaris. Effects on survival rate could not be observed consistently for both haplotypes throughout all B. vulgaris forms. This may in part be due to the absence of the two haplotypes in some of the other B. vulgaris forms (see Additional file 7). Similar to bolting rate, the absence of haplotype effects on survival rate in some B. vulgaris forms might also be due to the polygenic inheritance of survival rate.
Among the SNPs underlying haplotypes FL1a_H6, FL1b_H3, FL1b_H5, FL1b_H6, and FL1b_H10, two are located in an exon. The SNP in exon 3 is synonymous, whereas the SNP in exon 4 is non-synonymous, leading to an amino acid substitution from valine to isoleucine. The other SNPs identified in BvFL1 are silent as they are located in introns, including intron 1. These SNPs might influence gene function by affecting the transcriptional regulation of BvFL1, as was reported for intronic polymorphisms in FLC in Arabidopsis [51, 52]. Also, Heo and Sung  reported that the regulatory non-coding RNA COLDAIR is expressed from intron 1 of FLC. Finally, the increase in winter survival rates observed for an allelic variant of FLC was also associated with polymorphisms in intron 1. As with association studies in general, it cannot be excluded that the functional polymorphisms for the traits investigated are located outside the amplified gene regions and that the SNPs detected here are merely linked to these polymorphisms.
In conclusion, our study demonstrates that EcoTILLING can be successfully employed in B. vulgaris to survey a large panel of plant accessions for allelic variants in different candidate genes. Our data also provide the first genetic indication that an FLC homolog indeed may also affect flowering time (and winter survival) in a species which is only distantly related to A. thaliana. The above described panel of diverse B. vulgaris forms is an excellent resource to identify allelic variation in additional flowering time control genes such as BvFT2 or candidate genes for agronomic traits such as stress response and plant architecture. Allelic variants identified by EcoTILLING can be used to introduce new genetic variation into elite beet breeding material.
Plant material and phenotypic data
Phenotypic data for bolting before winter and winter hardiness were taken from a recent study described in detail by Kirchhoff et al. . In short, a panel of 396 B. vulgaris accessions covering a wide range of genetic diversity was tested for winter hardiness in a replicated overwintering field experiment in eight environments at five different locations in Germany and Belarus in the winters of 2008/09 and 2009/10. Survival rates were determined as the fraction of surviving individuals among all plants of a given accession ranging from 0 to 1, where 0 means no plants survived and 1 means all plants from one accessions survived. The mean survival rates were estimated as best linear unbiased predictors (BLUPs) for each accession across all environments. Accordingly, bolting rates before winter were determined in the 2009/10 environments as the fraction of bolting individuals among all plants of a given accession ranging from 0 to 1, where 0 means none of the plants bolted and 1 means all plants of a given accession bolted. Recording time was before the first frost (2 December 2009). To avoid unbalanced data, we reduced the data set to a subpanel of 268 accessions that were tested in all environments. These comprise the four cultivar groups fodder beet (40), leaf beet (47), garden beet (58) and sugar beet (88), as well as 35 B. vulgaris ssp. maritima accessions. The 88 sugar beets can be further subdivided into 49 elite accessions (sugar beet elite breeding material, SBEBM) provided by Strube GmbH & Co. KG (Söllingen, Germany) and 39 mostly gene bank accessions of various origins (sugar beet germplasm, SBGP).
DNA isolation and screening for polymorphisms
DNA was isolated from freeze dried leaf samples taken from up to eight plants per accession. This was done with a NucleoSpin® 96 Plant II Kit (MACHEREY-NAGEL GmbH & Co. KG, Düren, Germany) as recommended by the manufacturer on a Tecan “Freedom Evo” Robot. DNA concentration was measured via the Tecan Robot using a photometer and SYBR® Green (Invitrogen GmbH, Darmstadt, Germany) and normalized with DNase free water to a final concentration of 10 ng/μl in a total volume of 160 μl. The 268 DNA samples representing the 268 B.vulgaris accessions of the test panel were each pooled 1:1 with DNA of the biennial sugar beet 93161P as reference type and stored in 96 well plates. 93161P is an inbred line homozygous for the investigated candidate genes and was provided by Saatzucht Dieckmann.
Oligonucleotide primers amplifying conserved domains of the genes BvFT1, BvFL1 and BTC1 were designed from genomic sequences with FastPCR . Regions were chosen after analyses of genomic DNA sequence with CODDLE (Codons Optimized to Discover Deleterious Lesions; http://www.proweb.org/coddle/). The primers were pre-screened before labelling in a so called “crash-test” adapted from Weil and Monde . Forward and reverse primers were end dye labelled with Dyomics fluorescent tags DY-681 (700 nm absorption) and DY-781 (800 nm absorption), respectively. PCR amplification was done in a 20 μl volume containing 1 ng pooled DNA, 1 × Taq buffer, 1.5 mM MgCl2, 2 mM dNTPs (Invitrogen, Darmstadt, Germany), 0.2 units recombinant TAQ DNA polymerase (Invitrogen, Darmstadt, Germany) and 0.8 pmol primer (biomers.net, Ulm, Germany). The primers were used in a labelled versus non labelled ratio of 3:2 for DY-681 and 4:1 for DY-781 according to Till et al. . PCR was performed on a DNA Engine DYAD thermal cycler (MJ Research Inc., Waltham, MA, USA). PCR steps for the amplification were as follows: an initial denaturation step at 95°C for 5 min followed by 30 cycles of a 30 sec denaturation step at 95°C, 30 sec annealing at 60°C and 60 sec elongation at 72°C. The PCR was concluded with 5 min elongation at 72°C. The crude celery extract (CCE) was extracted as described by Till et al.  with a slight modification. We did not dialyze and re-buffer our celery juice as recommended. Instead we used the crude extract for enzymatic mismatch cleavage and tested it against commercial products. The results (data not shown) were identical with those obtained by using Surveyor® endonuclease. For SNP evaluation we only used the 700 nm channel of the LI-COR 4300, because the 800 nm channel did not provide additional information. For heteroduplex formation the PCR product was denatured at 95°C for 10 min and slowly re-annealed by cooling down to 85°C by 2° per sec and further cooling down to 25°C by 0.5°C per sec. The re-annealed PCR product was digested at 42°C for 15 min with crude celery extract (CCE) containing 0.6 μl CCE and 5.4 μl CCE buffer for each 20 μl reaction. The CCE buffer was prepared according to Till et al. . The reaction was stopped with 4 μl 200 mM EDTA. PCR products were cleaned up after endonuclease digestion by Sephadex purification.
Fragment analysis was performed on a LI-COR 4300 DNA analyser using a 6.5% KB Plus gel matrix (LI-COR®, Bad Homburg, Germany). The gel run was performed at 1,500 V, 40 mA and 40 W for 2 hours and 30 minutes. Acquired data were analysed visually using the software Gelbuddy . For each gel run, an analysis window smaller than the target amplicon size was manually chosen based on image quality and the absence of PCR mispriming artefacts that can occur near the primer binding region . For considering gel bands as digestion fragments, bands in the 700 nm channel were scored and a binary matrix was generated reflecting the presence (1) or absence (0) of all different fragment sizes for each sample.
Analysis of polymorphisms and haplotypes
For simplification, each unique fragment visible after acrylamide gel electrophoresis was considered as a SNP despite the fact that fragments could also be caused by small indels. SNP densities were calculated as the number of polymorphic SNP loci divided by the total length of screened sequence in kb. Non-reference nucleotide frequencies (NNFs) were calculated for each SNP locus as the number of accessions with the SNP allele deviating from the reference allele of 93161P divided by the number of screened accessions. Average heterozygosity Ht (i.e. gene diversity) for each SNP was calculated with the genetic distance and phylogenetic analysis package DISPAN .
Accessions with identical SNP pattern were assigned to the same haplotype. Accessions with no restriction bands on a LI-COR gel were assigned to the reference haplotype H0 (93161P). A haplotype with a frequency of less than 5% was declared rare. Non-reference haplotype frequencies (NHF) were calculated for each haplotype as the number of accessions with a haplotype deviating from the respective reference haplotype (FT1a_H0, FT1b_H0, FL1a_H0, FL1b_H0 or BTC1_H0) divided by the total number of accessions screened per B. vulgaris form.
Selected accessions with significant haplotypes associated with bolting rate and survival rate were sequenced via Sanger sequencing with the respective primer combination. To predict the functional impact of the SNPs characteristic of these haplotypes, the web based tools PARSESNP  and SIFT  were used.
AFLP analysis and population structure analysis
The population structure of the 268 B. vulgaris accessions was analysed with the AFLP (amplified fragment length polymorphism) technique essentially as described by Vos et al. . The following modification was applied: restriction of DNA was carried out with PstI instead of EcoRI.
Following AFLP marker analysis the population structure was calculated with the software package STRUCTURE version 2.3.4 [33, 59]. The optimum number of populations (k) was selected after six independent runs with a burn-in of 50,000 iterations followed by 100,000 iterations for each value of k (testing from k = 1 to k = 8). As program parameters for the investigation of the whole panel, the no-admixture model with the correlated allele frequency model was chosen. The most likely value for k was determined on the basis of the following criteria: (1) comparison of values for L(K) of each k; (2) stability of grouping patterns across five runs, and (3) value of ΔK calculated based on the second order rate of change of the likelihood (ΔK = m(|L’(K)|)/s[L(K)])  by the web based interface of STRUCTURE HARVESTER .
Association mapping was conducted using the general linear model (GLM) in TASSEL v. 3.0 . An association of a given amplicon with bolting rate, survival rate, or survival rate with bolting rate as cofactor was claimed at an experiment wise alpha level of 0.05 (Bonferroni correction). In case of a significant association of a given amplicon with bolting rate or survival rate, a Dunnett’s post hoc test for multiple comparisons of all haplotypes against the reference haplotype H0 was performed. This was performed with the statistical software R  separately for each B. vulgaris form.
The authors thank Michaela Jahn for her excellent technical assistance, and Monika Bruisch, and Erwin Danklefsen for substantial help and support in field and greenhouse work. Financial support of the project was provided by the Federal Ministry of Education and Research (BMBF) within the GABI program (FKZ: 0315052C and 0315058A). We thank Dr. Axel Schechert from Strube Research GmbH & Co.KG for providing seed material, and Gretel Schulze-Buxloh and Sebastian Vogt for providing the full-length genomic sequence of BvFT1. We are indebted to Sonja Hollmer from the Zoological Institute, Department of Structural Biology for preparing the crude celery extract, and Dr. Mario Hassler for support in statistical analysis. We also acknowledge the helpful comments of four unknown reviewers.
- McCallum CM, Comai L, Greene EA, Henikoff S: Targeting induced local lesions in genomes (TILLING) for plant functional genomics. Plant Physiol. 2000, 123: 439-442. 10.1104/pp.123.2.439.PubMedPubMed CentralView ArticleGoogle Scholar
- Comai L, Young K, Till BJ, Reynolds SH, Greene EA, Codomo CA, Enns LC, Johnson JE, Burtner C, Odden AR, Henikoff S: Efficient discovery of DNA polymorphisms in natural populations by EcoTILLING. Plant J. 2004, 37: 778-786. 10.1111/j.0960-7412.2003.01999.x.PubMedView ArticleGoogle Scholar
- Till BJ, Jankowicz-Cieslak J, Sígi L, Huynh O, Utsushi H, Swennen R, Terauchi R, Mba C: Discovery of nucleotide polymorphisms in the Musa gene pool by EcoTILLING. Theor Appl Genet. 2010, 121: 1381-1389. 10.1007/s00122-010-1395-5.PubMedPubMed CentralView ArticleGoogle Scholar
- Gilchrist EJ, Haughn GW, Ying CC, Otto SP, Zhuang J, Cheung D, Hamberger B, Aboutorabi F, Kalynyak T, Johnson L, Bohlmann J, Ellis B, Douglas CJ, Bronk QCB: Use of Ecotilling as an efficient SNP discovery tool to survey genetic variation in wild populations of Populus trichocarpa. Mol Ecol. 2006, 15: 1367-1378. 10.1111/j.1365-294X.2006.02885.x.PubMedView ArticleGoogle Scholar
- Galeano CH, Gomez M, Rodriguez LM, Blair MW: CEL I nuclease digestion for SNP discovery and marker development in common bean (Phaseolus vulgaris L). Crop Sci. 2009, 49: 381-394. 10.2135/cropsci2008.07.0413.View ArticleGoogle Scholar
- Barkley NA, Wang M, Gillaspie A, Dean R, Pederson G, Jenkins T: Discovering and verifying DNA polymorphisms in a mung bean [V. radiata (L.) R. Wilczek] collection by EcoTILLING and sequencing. BMC Res Notes. 2008, 1: 28. 10.1186/1756-0500-1-28.PubMedPubMed CentralView ArticleGoogle Scholar
- Mejlhede N, Kyjovska Z, Backes G, Burhenne K, Rasmussen SK, Jahoor A: EcoTILLING for the identification of allelic variation in the powdery mildew resistance genes mlo and Mla of barley. Plant Breeding. 2006, 125: 461-467. 10.1111/j.1439-0523.2006.01226.x.View ArticleGoogle Scholar
- Kadaru S, Yadav A, Fjellstrom R, Oard J: Alternative ecotilling protocol for rapid, cost-effective single-nucleotide polymorphism discovery and genotyping in rice (Oryza sativa L). Plant Mol Biol Rep. 2006, 24: 3-22. 10.1007/BF02914042.View ArticleGoogle Scholar
- Negrao S, Almadanim C, Pires I, McNally K, Oliveira M: Use of EcoTILLING to identify natural allelic variants of rice candidate genes involved in salinity tolerance. Plant Genetic Res. 2011, 9: 300-304. 10.1017/S1479262111000566.View ArticleGoogle Scholar
- Elias R, Till BJ, Mba C, Al-Safadi B: Optimizing TILLING and EcoTILLING techniques for potato (Solanum tuberosum L). BMC Res Notes. 2009, 2: 1-5. 10.1186/1756-0500-2-1.View ArticleGoogle Scholar
- Nieto C, Piron F, Dalmais M, Marco CF, Moriones E, Gomez-Guillamon ML, Truniger V, Gomez P, Garcia-Mas J, Aranda MA, Bendahmane A: EcoTILLING for the identification of allelic variants of melon eIF4E, a factor that controls virus susceptibility. BMC Plant Biol. 2007, 7: 1-9. 10.1186/1471-2229-7-1.View ArticleGoogle Scholar
- Rigola D, van Oeveren J, Janssen A, Bonne A, Schneiders H, van der Poel HJ, van Orsouw NJ, Hogers RC, de Both MT, van Eijk MJ: High-throughput detection of induced mutations and natural variation using keypoint technology. PLoS One. 2009, 4: e4761. 10.1371/journal.pone.0004761.PubMedPubMed CentralView ArticleGoogle Scholar
- FAOSTAT - food and agricultural commodities production world wide. 2010, http://faostat.fao.org/site/339/default.aspx.
- Letschert J, Lange W, Frese L, Van Den Berg R: Taxonomy of Beta section Beta.Journal of sugar beet research. 1994, 31 (1-2): 69-85.View ArticleGoogle Scholar
- Hoffmann CM, Kluge-Severin S: Light absorption and radiation use efficiency of autumn and spring sown sugar beets. Field Crops Res. 2010, 119: 238-244. 10.1016/j.fcr.2010.07.014.View ArticleGoogle Scholar
- Jung C, Müller AE: Flowering time control and applications in plant breeding. Trends Plant Sci. 2009, 14: 563-573. 10.1016/j.tplants.2009.07.005.PubMedView ArticleGoogle Scholar
- Pin PA, Benlloch R, Bonnet D, Wremerth-Weich E, Kraft T, Gielen JJ, Nilsson O: An antagonistic pair of FT homologs mediates the control of flowering time in sugar beet. Science. 2010, 330: 1397-1400. 10.1126/science.1197004.PubMedView ArticleGoogle Scholar
- Pin PA, Zhang W, Vogt SH, Dally N, Büttner B, Schulze-Buxloh G, Jelly NS, Chia TYP, Mutasa-Göttgens ES, Dohm JC, Himmelbauer H, Weisshaar B, Kraus J, Gielen JJL, Lommel M, Weyens G, Wahl B, Schechert A, Nilsson O, Jung C, Kraft T, Müller AE: The role of a pseudo-response regulator gene in life cycle adaptation and domestication of beet. Curr Biol. 2012, 22: 1095-1101. 10.1016/j.cub.2012.04.007.PubMedView ArticleGoogle Scholar
- Fowler DB, Limin AE, Robertson AJ, Gusta LV: International crop science I. Breeding for low-temperature tolerance in field crops. Ames, Iowa: Crops Science Society of America; 1993: 357-362.Google Scholar
- Van Dijk H, Boudry P: Genetic variation for life histories in Beta maritima. Int Board Plant Genetic Res. 1991, 7: 44-55.Google Scholar
- Van Dijk H, Boudry P, McCombre H, Vernet P: Flowering time in wild beet (Beta vulgaris ssp. maritima) along a latitudinal cline. Acta Oecologica. 1997, 18: 47-60. 10.1016/S1146-609X(97)80080-X.View ArticleGoogle Scholar
- Michaels SD, Amasino RM: FLOWERING LOCUS C encodes a novel MADS domain protein that acts as a repressor of flowering. Plant Cell. 1999, 11: 949-956.PubMedPubMed CentralView ArticleGoogle Scholar
- Sheldon CC, Burn JE, Perez PP, Metzger J, Edwards JA, Peacock WJ, Dennis ES: The FLF MADS box bene: a repressor of flowering in arabidopsis regulated by vernalization and methylation. Plant Cell. 1999, 11: 445-458.PubMedPubMed CentralView ArticleGoogle Scholar
- Reeves PA, He YH, Schmitz RJ, Amasino RM, Panella LW, Richards CM: Evolutionary conservation of the flowering locus C-mediated vernalization response: evidence from the sugar beet (Beta vulgaris). Genetics. 2007, 176: 295-307. 10.1534/genetics.106.069336.PubMedPubMed CentralView ArticleGoogle Scholar
- Bieniawska Z, Espinoza C, Schlereth A, Sulpice R, Hincha DK, Hannah MA: Disruption of the Arabidopsis circadian clock is responsible for extensive variation in the cold-responsive transcriptome. Plant Physiol. 2008, 147: 263-279. 10.1104/pp.108.118059.PubMedPubMed CentralView ArticleGoogle Scholar
- Eriksson ME, Webb AAR: Plant cell responses to cold are all about timing. Current Opinion Plant Biol. 2011, 14: 731-737. 10.1016/j.pbi.2011.08.005.View ArticleGoogle Scholar
- Franklin KA, Whitelam GC: Light-quality regulation of freezing tolerance in Arabidopsis thaliana. Nat Genet. 2007, 39: 1410-1413. 10.1038/ng.2007.3.PubMedView ArticleGoogle Scholar
- Gilmour SJ, Fowler SG, Thomashow MF: Arabidopsis transcriptional activators CBF1, CBF2, and CBF3 have matching functional activities. Plant Mol Biol. 2004, 54: 767-781.PubMedView ArticleGoogle Scholar
- Seo E, Lee H, Jeon J, Park H, Kim J, Noh YS, Lee I: Crosstalk between cold response and flowering in arabidopsis is mediated through the flowering-time gene SOC1 and its upstream negative regulator FLC. Plant Cell. 2009, 21: 3185-3197. 10.1105/tpc.108.063883.PubMedPubMed CentralView ArticleGoogle Scholar
- Deng W, Ying H, Helliwell CA, Taylor JM, Peacock WJ, Dennis ES: Flowering locus C (FLC) regulates development pathways throughout the life cycle of Arabidopsis. Proc Natl Acad Sci. 2011, 108: 6680-6685. 10.1073/pnas.1103175108.PubMedPubMed CentralView ArticleGoogle Scholar
- Heo JB, Sung S: Vernalization-mediated epigenetic silencing by a long intronic noncoding RNA. Science. 2011, 331: 76-79. 10.1126/science.1197349.PubMedView ArticleGoogle Scholar
- Kirchhoff M, Svirshchevskaya A, Hoffmann C, Schechert A, Jung C, Kopisch-Obuch F: High degree of genetic variation of winter hardiness in a panel of Beta vulgaris L. Crop Sci. 2012, 52: 179-188. 10.2135/cropsci2011.04.0185.View ArticleGoogle Scholar
- Pritchard JK, Stephens M, Donnelly P: Inference of population structure using multilocus genotype data. Genetics. 2000, 155: 945-959.PubMedPubMed CentralGoogle Scholar
- Weil CF, Monde RA: Getting the point - mutations in maize. Crop Sci. 2007, 47: S60-S67.View ArticleGoogle Scholar
- Schmitz RJ, Amasino RM: Vernalization: a model for investigating epigenetics and eukaryotic gene regulation in plants. Biochimica et Biophysica Acta (BBA) Gene Struct Expression. 2007, 1769: 269-275. 10.1016/j.bbaexp.2007.02.003.View ArticleGoogle Scholar
- Till BJ, Colbert T, Codomo C, Enns L, Johnson J, Reynolds SH, Henikoff JG, Greene EA, Steine MN, Comai L, Henikoff S: High-throughput tilling for Arabidopsis. Methods Mol Biol. 2006, 323: 127-135.PubMedGoogle Scholar
- Greene EA, Codomo CA, Taylor NE, Henikoff JG, Till BJ, Reynolds SH, Enns LC, Burtner C, Johnson JE, Odden AR, Comai L, Henikoff S: Spectrum of chemically induced mutations from a large-scale reverse-genetic screen in arabidopsis. Genetics. 2003, 164: 731-740.PubMedPubMed CentralGoogle Scholar
- Till BJ, Zerr T, Comai L, Henikoff S: A protocol for tilling and Ecotilling in plants and animals. Nat Protoc. 2006, 1: 2465-2477. 10.1038/nprot.2006.329.PubMedView ArticleGoogle Scholar
- Till BJ, Zerr T, Bowers E, Greene EA, Comai L, Henikoff S: High-throughput discovery of rare human nucleotide polymorphisms by EcoTILLING. Nucleic Acids Res. 2006, 34: e99. 10.1093/nar/gkl479. 01–12PubMedPubMed CentralView ArticleGoogle Scholar
- Jung C, Pillen K, Frese L, Melchinger A: Phylogenetic relationships between cultivated and wild species of the genus Beta revealed by DNA “fingerprinting”. Theor Appl Genet. 1993, 86: 449-457. 10.1007/BF00838560.PubMedView ArticleGoogle Scholar
- McGrath JM, Derrico CA, Yu Y: Genetic diversity in selected, historical US sugarbeet germplasm and Beta vulgaris ssp. maritima. Theor Appl Genet. 1999, 98: 968-976. 10.1007/s001220051157.View ArticleGoogle Scholar
- Fenart S, Arnaud JF, De Cauwer I, Cuguen J: Nuclear and cytoplasmic genetic diversity in weed beet and sugar beet accessions compared to wild relatives: new insights into the genetic relationships within the Beta vulgaris complex species. Theor Appl Genet. 2008, 116: 1063-1077. 10.1007/s00122-008-0735-1.PubMedView ArticleGoogle Scholar
- Lewellen RT: Use of plant introductions to improve populations and hybrids of sugar beet. Use of plant introductions in cultivar development. Edited by: Shands HL. 1992, Madison, WI: Crop Sci. Soc. Am, 117-135. Volume 2Google Scholar
- Savitsky VF: Proceedings american society of sugar beet technologists. A genetic study of monogerm and multigerm characters in beets. Salt Lake City, Utah, USA: The American Society of Sugar Beet Technologists; 1952: 824.Google Scholar
- Biancardi E, Lewellen RT, De Biaggi M, Erichsen A, Stevanato P: The origin of rhizomania resistance in sugar beet. Euphytica. 2002, 127: 383-397. 10.1023/A:1020310718166.View ArticleGoogle Scholar
- Li J, Luhmann A-K, Weißleder K, Stich B: Genome-wide distribution of genetic diversity and linkage disequilibrium in elite sugar beet germplasm. BMC Genomics. 2011, 12: 1-10. 10.1186/1471-2164-12-1.PubMedPubMed CentralView ArticleGoogle Scholar
- Michaels SD, He YH, Scortecci KC, Amasino RM: Attenuation of flowering locus C activity as a mechanism for the evolution of summer-annual flowering behavior in Arabidopsis. Proc Natl Acad Sci. 2003, 100: 10102-10107. 10.1073/pnas.1531467100.PubMedPubMed CentralView ArticleGoogle Scholar
- Biancardi E, McGrath JM, Panella L, Lewellen RT, Stevanato P: Sugar beet. Handbook of plant breeding: root and tuber crops. Edited by: Bradshaw JE. New York: Springer; 2010: 173-219. Prohens J, Nuez F, Carena MJ (Series Editor): HANDBOOK OF PLANT BREEDINGView ArticleGoogle Scholar
- McGrath JM, Saccomani M, Stevanato P, Biancardi E: Beet. Genome mapping and molecular breeding in plants. Edited by: Kole C. Berlin Heidelberg: Springer; 2007: 191-202.Google Scholar
- Korves TM, Schmid KJ, Caicedo AL, Mays C, Stinchcombe JR, Purugganan MD, Schmitt J: Fitness effects associated with the major flowering time gene FRIGIDA in Arabidopsis thaliana in the field. Am Nat. 2007, 169: E141-E157. 10.1086/513111.PubMedView ArticleGoogle Scholar
- Sheldon CC, Conn AB, Dennis ES, Peacock WJ: Different regulatory regions are required for the vernalization-induced repression of flowering locus C and for the epigenetic maintenance of repression. Plant Cell. 2002, 14: 2527-2537. 10.1105/tpc.004564.PubMedPubMed CentralView ArticleGoogle Scholar
- Choi K, Kim J, Hwang H-J, Kim SY, Park C, Kim SY, Lee I: The FRIGIDA complex activates transcription of FLC, a strong flowering repressor in arabidopsis, by recruiting chromatin modification factors. Plant Cell. 2011, 23: 289-303. 10.1105/tpc.110.075911.PubMedPubMed CentralView ArticleGoogle Scholar
- Kalendar R, Lee D, Schulman AH: FastPCR software for PCR primer and probe design and repeat search. Genes, Genomes Genomics. 2009, 3: 1-14.Google Scholar
- Zerr T, Henikoff S: Automated band mapping in electrophoretic gel images using background information. Nucleic Acids Res. 2005, 33: 2806-2812. 10.1093/nar/gki580.PubMedPubMed CentralView ArticleGoogle Scholar
- Ota T: DISPAN: genetic distance and phylogenetic analysis. Book DISPAN: genetic distance and phylogenetic analysis. Edited by: City: Institute of Molecular Evolutionary Genetics. University Park, PA, USA: Pennsylvania State University; 1993: 10Google Scholar
- Taylor NE, Greene EA: PARSESNP: a tool for the analysis of nucleotide polymorphisms. Nucleic Acids Res. 2003, 31: 3808-3811. 10.1093/nar/gkg574.PubMedPubMed CentralView ArticleGoogle Scholar
- Ng PC, Henikoff S: SIFT: predicting amino acid changes that affect protein function. Nucleic Acids Res. 2003, 31: 3812-3814. 10.1093/nar/gkg509.PubMedPubMed CentralView ArticleGoogle Scholar
- Vos P, Hogers R, Bleeker M, Reijans M, Lee T, Hornes M, Friters A, Pot J, Paleman J, Kuiper M, Zabeau M: AFLP: a new technique for DNA fingerprinting. Nucleic Acids Res. 1995, 23: 4407-4414. 10.1093/nar/23.21.4407.PubMedPubMed CentralView ArticleGoogle Scholar
- Falush D, Stephens M, Pritchard JK: Inference of population structure using multilocus genotype data: dominant markers and null alleles. Mol Ecol Notes. 2007, 7: 574-578. 10.1111/j.1471-8286.2007.01758.x.PubMedPubMed CentralView ArticleGoogle Scholar
- Evanno G, Regnaut S, Goudet J: Detecting the number of clusters of individuals using the software structure: a simulation study. Mol Ecol. 2005, 14: 2611-2620. 10.1111/j.1365-294X.2005.02553.x.PubMedView ArticleGoogle Scholar
- Earl D, von Holdt B: Structure harvester: a website and program for visualizing structure output and implementing the evanno method. Conserv Genet Resour. 2012, 4: 359-361. 10.1007/s12686-011-9548-7.View ArticleGoogle Scholar
- Bradbury PJ, Zhang Z, Kroon DE, Casstevens TM, Ramdoss Y, Buckler ES: TASSEL: software for association mapping of complex traits in diverse samples. Bioinformatics. 2007, 23: 2633-2635. 10.1093/bioinformatics/btm308.PubMedView ArticleGoogle Scholar
- R Development Core Team: Book R: a language and environment for statistical computing. R: a language and environment for statistical computing. City: R Foundation for Statistical Computing; 2011: 12Google Scholar