Open Access

EcoTILLING in Beta vulgaris reveals polymorphisms in the FLC-like gene BvFL1that are associated with annuality and winter hardiness

  • Sebastian LM Frerichmann1,
  • Martin Kirchhoff1, 2,
  • Andreas E Müller1, 3,
  • Axel J Scheidig4,
  • Christian Jung1 and
  • Friedrich J Kopisch-Obuch1Email author
BMC Plant Biology201313:52

DOI: 10.1186/1471-2229-13-52

Received: 12 September 2012

Accepted: 21 February 2013

Published: 25 March 2013

Abstract

Background

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.

Results

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.

Conclusions

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.

Keywords

Beta vulgaris EcoTILLING Diversity Flowering Winter hardiness BvFL1

Background

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 [1]. 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) [2], Musa spp. (various banana species) [3], Populus trichocarpa (black cottonwood) [4], Phaseolus vulgaris (common bean) [5], and Vigna radiata (mung bean) [6]. Furthermore, it has been used for candidate gene-based detection of new alleles conferring resistance to biotic and abiotic stress in Hordeum vulgare (barley) [7], Oryza sativa (rice) [8, 9], Solanum tuberosum (potato) [10]. Cucumis spp. (including cucumber) [11] and Solanum lycopersicum (tomato) [12]. EcoTILLING has not been reported for sugar beet (Beta vulgaris ssp. vulgaris L.) which contributes to 22% of the world production of white sugar [13].

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) [14]. 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 [15]. 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 [16]. 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 [16].

In order to avoid an untimely transition to the extremely cold-sensitive generative phase [19] 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 [17]. Similar to FLC, BvFT1 is gradually down-regulated during the prolonged cold of winter [17]. 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[18]. 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 [18]. 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 [18].

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 [24]. 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 [2527]. 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 [28] and elevated levels of FLC expression [29]. Deng et al. [30] 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[31].

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 [32]. 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.

Results

Phenotyping and model-based analyses of population structure

In fall 2009, 41 out of the 268 accessions sown in the field in July or August at four different locations, had at least one plant which had started bolting before the first frost. These accessions included two sugar beets (2.2% of accessions tested), four garden beets (6.9%), four fodder beets (10.0%), 20 leaf beets (42.3%) and eleven wild sea beets (B. vulgaris ssp. maritima) (31.4%). Across all four environments, bolting rates for the 41 accessions ranged from 0.05 to 0.75 (Figure 1). Variation for survival rate survival rate after winter across eight environments in 2008/09 and 2009/10 was described in detail by Kirchhoff et al. [32], and ranged from 0.07 to 0.66. On average, sugar beet accessions performed best (0.39) while fodder beet and garden beet performed worst (0.24 and 0.19, respectively). The largest variation for survival rate was found in B. vulgaris ssp. maritima followed by leaf beets, whereas sugar beets showed the smallest variation. Population structure was analysed by a model-based method using the genotypic data for 40 polymorphic AFLPs detected among the 268 accessions. An independent calculation of k was repeated six times for each value of k from k = 1 to k = 8. The log probability L(K) increased sharply from k = 1 to k = 3, but only slowly after k = 3 (Additional file 1). When k is approaching a true value, L(K) plateaus or continues to increase slightly [33]. Therefore the structure analysis suggested the presence of three subgroups (k = 3), where most of the sugar beets (0.81) fell in the first group, most of the fodder beets (0.61) and garden beets (0.86) in the second group, and most of the leaf beets and B. vulgaris ssp. maritima in the third group (Figure 2). To further increase confidence in the k value estimate, we calculated ΔK and obtained the highest ΔK value (60.04) for k = 3 (Additional file 2).
https://static-content.springer.com/image/art%3A10.1186%2F1471-2229-13-52/MediaObjects/12870_2012_Article_1270_Fig1_HTML.jpg
Figure 1

Distribution of bolting rate before winter. Distribution of bolting rate before winter among a subset of 41 B. vulgaris accessions tested in four different environments. BVM = B. vulgaris ssp. maritima.

https://static-content.springer.com/image/art%3A10.1186%2F1471-2229-13-52/MediaObjects/12870_2012_Article_1270_Fig2_HTML.jpg
Figure 2

Overview of population structure analysis. Population structure of 268 B. vulgaris accessions based on 40 AFLP markers under the assumption of k = 3 subpopulations. Each B. vulgaris accession is represented by one bar that is divided in up to k segments, each proportional to the inferred subpopulation. Accessions are grouped by their respective B. vulgaris group. BVM = B. vulgaris ssp. maritima.

Amplification of candidate genes

For the three candidate genes BvFL1, BvFT1 and BTC1 we designed 24 primer pairs, of which five primer pairs passed the primer “crash-test” and therefore were suitable for EcoTILLING. The genomic sequence of BvFL1 [GenBank: DQ189214.1 and DQ189215.1] was taken from Reeves et al. [24] while the genomic sequence of BvFT1 [GenBank: HM448909.1] and BTC1 [GenBank: HQ709091.1] were taken from Pin et al. [17, 18]. We adopted a primer pre-screen as described by Weil and Monde [34] to test for the occurrence of unwanted amplification from single primers prior to the costly synthesis of labelled primers. Amplification from the five successfully tested primer pairs resulted in a total amplicon length of 4,234 bp (Table 1). The remaining primer combinations were not suitable for EcoTILLING due to miss-priming and single primer amplification revealed in the primer crash-test (data not shown).
Table 1

Overview of candidate genes investigated with EcoTILLING

Gene

Genomic size

Protein domain(s)

Amplicon size(s)

Primer name and sequence (5’-3’)

Length and portion of genomic sequence covered by EcoTILLING

Portion of reading frame covered by EcoTILLING

BvFL1

8.62 kb

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%)

62%

BvFT1

7.53 kb

PEBP

FT1a: 916 bp FT1b: 713 bp

FT1a-fw tggtacgtgtatgaaacagaagctg FT1a-rv catcaactccatatttggggtg FT1b-fw acccatctatacttgtcgatgacc FT1b-rv caatggggaagtggttcacact

1,629 bp (22%)

92%

BTC1

11.27 kb

REC, CCT

BTC1: 996 bp

BTC1-fw cagctgtaggatgttatcgtgctgag BTC1-rv agtaggtgataaggacaagacattgc

992 bp (9%)

15%

The gene names, genomic sizes in kilo base pairs (kb), protein domains, sizes of amplicons in base pairs (bp) per gene, primer names and sequences, and the sizes of the amplicons in base pairs (bp) are given. Genomic size was defined here as the size of the genomic sequence of a gene from the start to the stop codon plus 1 kb each upstream and downstream. The genomic size given for BvFL1 refers to the known portion of the sequence. A large part of intron 1 was not sequenced (Reeves et al. [24]) and is not considered here. Furthermore, the length of the investigated genomic sequence and the percentage of the open reading frame which was investigated via EcoTILLING are given.

For BvFL1 we could amplify two regions (Figure 3). The first amplicon (‘FL1a’) covered 977 bp of the promoter, exon 1 (including the 5'-UTR and the 5' region of the coding sequence), and part of intron 1. The second amplicon (‘FL1b’) spanned a region of 632 bp in exons 3 and 4 and the intervening intron. These regions were chosen because they contain the promoter region, the region encoding the MADS box domain as well as a TGTGAT sequence motif (K box) which is associated with transcription factor binding activity. Thus, 62% of the BvFL1 ORF was covered. Furthermore, we targeted the first intron, which is known to include a number of regulatory regions in Arabidopsis[35]. For BvFT1 we amplified a 916 bp region (‘FT1a’) that extends from the 5'-UTR to the 5' region of the coding sequence in exon 1. A second amplicon (‘FT1b’) spanned 713 bp and was located in exon 4 and the 3'-UTR (Figure 3). Both regions together cover 92% of the ORF and were chosen because they contain parts of the promoter and the PEB domain (Table 1). In BTC1 we amplified a 992 bp fragment (‘BTC1’) of the conserved response regulator receiver (REC) domain region (Figure 3), which covered 15% of the ORF (Table 1).
https://static-content.springer.com/image/art%3A10.1186%2F1471-2229-13-52/MediaObjects/12870_2012_Article_1270_Fig3_HTML.jpg
Figure 3

Flowering time genes BvFL1, BvFT1 and BTC1 and distribution of polymorphisms. Shown are graphical outputs from PARSESNP. a), c), e) Genomic sequences with intron/exon structure. b), d), f) Coding sequences. Red dotted arrowed lines indicate the amplicons investigated by EcoTILLING. Blue striped boxes indicate exons. White arrowheads that point downwards indicate changes in non-coding regions. Black arrowheads that point upwards indicate changes that induce missense mutations in the protein products. White arrowheads that point upwards indicate silent changes. Green boxes mark important protein domains.

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.

The mean non-reference nucleotide frequencies (NNFs; s. Materials and Methods) for BvFL1 (0.18) and BvFT1 (0.17) are similar across all accessions tested, but varied between individual B. vulgaris forms, ranging from 0.12 in garden beets to 0.23 in B. vulgaris ssp. maritima for BvFT1, and from 0.04 in garden beets to 0.55 in B. vulgaris ssp. maritima for BvFL1 (Additional file 4). The mean NNF for BTC1 over all B. vulgaris forms was 0.07, and ranged from 0.03 in fodder beets to 0.30 in B. vulgaris ssp. maritima. We identified 55 haplotypes across all five amplified candidate regions. The numbers of haplotypes ranged from four in amplicon BTC1 to 18 in amplicon FL1b. Sixty per cent of the detected haplotypes were rare, occurring at frequencies below 0.05. The reference haplotype (H0) was the most common in each amplicon, with frequencies ranging from 0.49 for FT1b_H0 to 0.87 for BTC1_H0 (Figure 4). The non-reference haplotype frequencies (NHF) ranged from 0.01 to 0.38 across the amplicons within the distinct B. vulgaris forms (Figure 4). The highest NHF was observed in B. vulgaris ssp. maritima, while garden and sugar beets had the lowest NHF. Gene diversity (Ht) for each amplicon and within each cultivar group is displayed in Additional file 5. Ht range was lowest in the amplicon FT1b (0.16 to 0.30) and highest in the amplicon FL1a (0.02 to 0.51). The highest and lowest diversity for the B. vulgaris forms was observed in amplicon FL1a for B. vulgaris ssp. maritima (0.51) and garden beets (0.02), respectively. Subdividing the 88 investigated sugar beet accessions into 49 accessions of elite breeding material (SBEBM) provided by Strube GmbH & Co. KG (Söllingen, Germany) and the remaining 39 accessions of sugar beet germplasm (SBGP, mostly composed of various gene bank accessions) revealed a trend towards lower diversity in the amplicons FL1a, FL1b and BTC1, while the diversity increased in FT1a and FT1b (Additional file 6).
https://static-content.springer.com/image/art%3A10.1186%2F1471-2229-13-52/MediaObjects/12870_2012_Article_1270_Fig4_HTML.jpg
Figure 4

Haplotype diversity diagram. Shown is the haplotype distribution for the amplicons FL1a, FL1b, FT1a, FT1b and BTC1. The percentage of haplotypes for each amplicon is given for each B. vulgaris form. Haplotypes are colour coded. The reference haplotype H0 is always shown in light blue, while non-reference haplotypes are shown in various colours. The B. vulgaris form is given on the x-axis. SBEBM = Sugar beet elite breeding material, SBGP = Sugar beet germplasm, BVM = B. vulgaris ssp. maritima.

BvFL1sequence variations are associated with bolting and survival rate

Based on the Q matrix for k = 3, associations with bolting rate, survival rate, and survival rate with bolting rate as cofactor were each significant (P ≤ 0.05) for the amplicon FL1a and for the amplicon FL1b (Table 2). Dunnet comparisons revealed that B. vulgaris ssp. maritima accessions with haplotype FL1a_H6, FL1b_H5 and FL1b_H10 had a significantly higher bolting rate of 55% (P < 0.0001), 75% (P < 0.0001) and 11% (P < 0.05), respectively, compared to 1% for FL1a_H0 and 2% for FL1b_H0 (Additional file 6). Furthermore, garden beet accessions with haplotype FL1b_H6 bolted before winter with a bolting rate of 6% (P < 0.0001) compared to 1% for accessions with the reference haplotype FL1b_H0 (Figure 5). Dunnet comparisons for survival rate revealed that B. vulgaris ssp. maritima accessions with haplotype FL1a_H6 had a significantly lower survival rate of 13% (P = 0.015) compared to 39% observed for accessions with the reference haplotype FL1a_H0 (Additional file 7). By contrast, leaf beet accessions with the haplotype FL1b_H3 had a significantly higher survival rate of 37% (P = 0.012) compared to 19% of accessions with the reference haplotype FL1b_H0 (Figure 6). DNA sequences of significant haplotypes of BvFL1 are shown in Additional file 8.
https://static-content.springer.com/image/art%3A10.1186%2F1471-2229-13-52/MediaObjects/12870_2012_Article_1270_Fig5_HTML.jpg
Figure 5

BvFL1 haplotype effect on bolting rate before winter. Graphic representation of the effect of BvFL1 haplotypes of the amplicons FL1a and FL1b on bolting rate before winter in five B. vulgaris forms. * or *** mark non-reference haplotypes with a significantly different bolting rate from reference haplotypes FL1a_H0 or FL1b_H0 (P < 0.05 or P < 0.0001, respectively); BVM = B. vulgaris ssp. maritima.

https://static-content.springer.com/image/art%3A10.1186%2F1471-2229-13-52/MediaObjects/12870_2012_Article_1270_Fig6_HTML.jpg
Figure 6

BvFL1 haplotype effect on survival rate after winter. Graphic representation of the effect of BvFL1 haplotypes of the amplicons FL1a and FL1b on survival rate after winter in five B. vulgaris forms. * marks non-reference haplotypes with a significantly different survival rate from the reference haplotypes FL1a_H0 or FL1b_H0 (alpha = 0.05); BVM = B. vulgaris ssp. maritima.

Table 2

Significant associations of bolting and survival rate with amplicon

Trait

Amplicon

Fd)

Pe)

R2 f)

BRa)

FL1_a

9.04**

6.47E-11

0.22

BR

FL1_b

11.28**

2.65E-23

0.44

SRb)

FL1_a

7.17**

1.44E-08

0.13

SR

FL1_b

2.91**

7.72E-05

0.13

SR/BRc)

FL1_a

6.54**

9.52E-08

0.11

SR/BR

FL1_b

2.90**

1.14E-04

0.12

a) BR = bolting rate, b) SR = survival rate c) SR with BR as cofactor, d) F value from the F test on marker, e) Bonferroni corrected p value, f) R2 is the fraction of the total variation explained by the amplicon.

**(P < 0.001) indicates the amplicon is highly significantly associated with trait.

Discussion

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. [36] and Galeano et al. [5] 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 [39] and banana [3] 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. [40] and by McGrath et al. [41] 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. [40] and Fénart et al. [42] 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 [4345]. 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. [46] and McGrath et al. [41]. 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 [17]. 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[47], the role of FLC-like genes outside the Brassicaceae is not well understood [16], 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 [18], 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 [17], 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 [50]. 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 [19]. 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. [29] 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 [31] 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[50] 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.

Conclusions

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.

Methods

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. [32]. 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 [53]. 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 [34]. 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. [38]. 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. [38] 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. [38]. 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 [54]. 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 [38]. 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 [55].

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 [56] and SIFT [57] 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. [58]. 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)]) [60] by the web based interface of STRUCTURE HARVESTER [61].

Statistical analysis

Association mapping was conducted using the general linear model (GLM) in TASSEL v. 3.0 [62]. 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 [63] separately for each B. vulgaris form.

Declarations

Acknowledgements

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.

Authors’ Affiliations

(1)
Plant Breeding Institute, Christian-Albrechts-University of Kiel
(2)
Nordsaat Saatzucht GmbH, Böhnshauser Straße
(3)
Strube Research GmbH & Co. KG
(4)
Zoological Institute, Department of Structural Biology, Christian-Albrechts-University of Kiel

References

  1. 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
  2. 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
  3. 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
  4. 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
  5. 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
  6. 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
  7. 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
  8. 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
  9. 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
  10. 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
  11. 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
  12. 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
  13. FAOSTAT - food and agricultural commodities production world wide. 2010, http://faostat.fao.org/site/339/default.aspx.
  14. 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
  15. 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
  16. 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
  17. 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
  18. 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
  19. 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
  20. Van Dijk H, Boudry P: Genetic variation for life histories in Beta maritima. Int Board Plant Genetic Res. 1991, 7: 44-55.Google Scholar
  21. 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
  22. 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
  23. 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
  24. 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
  25. 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
  26. 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
  27. 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
  28. 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
  29. 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
  30. 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
  31. 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
  32. 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
  33. Pritchard JK, Stephens M, Donnelly P: Inference of population structure using multilocus genotype data. Genetics. 2000, 155: 945-959.PubMedPubMed CentralGoogle Scholar
  34. Weil CF, Monde RA: Getting the point - mutations in maize. Crop Sci. 2007, 47: S60-S67.View ArticleGoogle Scholar
  35. 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
  36. 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
  37. 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
  38. 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
  39. 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
  40. 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
  41. 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
  42. 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
  43. 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
  44. 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
  45. 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
  46. 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
  47. 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
  48. 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
  49. 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
  50. 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
  51. 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
  52. 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
  53. 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
  54. 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
  55. 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
  56. 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
  57. 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
  58. 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
  59. 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
  60. 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
  61. 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
  62. 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
  63. 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

Copyright

© Frerichmann et al.; licensee BioMed Central Ltd. 2013

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.