High levels of nucleotide diversity and fast decline of linkage disequilibrium in rye (Secale cerealeL.) genes involved in frost response
© Li et al; licensee BioMed Central Ltd. 2011
Received: 12 August 2010
Accepted: 10 January 2011
Published: 10 January 2011
Rye (Secale cereale L.) is the most frost tolerant cereal species. As an outcrossing species, rye exhibits high levels of intraspecific diversity, which makes it well-suited for allele mining in genes involved in the frost responsive network. For investigating genetic diversity and the extent of linkage disequilibrium (LD) we analyzed eleven candidate genes and 37 microsatellite markers in 201 lines from five Eastern and Middle European rye populations.
A total of 147 single nucleotide polymorphisms (SNPs) and nine insertion-deletion polymorphisms were found within 7,639 bp of DNA sequence from eleven candidate genes, resulting in an average SNP frequency of 1 SNP/52 bp. Nucleotide and haplotype diversity of candidate genes were high with average values π = 5.6 × 10-3 and Hd = 0.59, respectively. According to an analysis of molecular variance (AMOVA), most of the genetic variation was found between individuals within populations. Haplotype frequencies varied markedly between the candidate genes. ScCbf14, ScVrn1, and ScDhn1 were dominated by a single haplotype, while the other 8 genes (ScCbf2, ScCbf6, ScCbf9b, ScCbf11, ScCbf12, ScCbf15, ScIce2, and ScDhn3) had a more balanced haplotype frequency distribution. Intra-genic LD decayed rapidly, within approximately 520 bp on average. Genome-wide LD based on microsatellites was low.
The Middle European population did not differ substantially from the four Eastern European populations in terms of haplotype frequencies or in the level of nucleotide diversity. The low LD in rye compared to self-pollinating species promises a high resolution in genome-wide association mapping. SNPs discovered in the promoters or coding regions, which attribute to non-synonymous substitutions, are suitable candidates for association mapping.
Rye (Secale cereale L.) is a cross-pollinated cereal with a diploid genome. It is grown on approximately 6 million hectares in Europe for bread-making, animal feed, forage feeding, and vodka production (FAO, 2010). As the most frost tolerant small grain cereal  it is well-suited for investigations of frost tolerance. Findings in rye are of interest for less frost tolerant cereals such as wheat and barley.
Cold and frost stress, namely chilling injury at temperatures lower than 10°C and freezing injury at temperatures lower than 0°C, adversely affect plant growth and productivity via cellular damage, dehydration and metabolic reaction slow-down. A major focus of this study was to investigate candidate genes with a putative role in frost tolerance. Frost tolerance has a polygenic inheritance. Many genes involved in the cold/frost responsive network have been identified in Arabidopsis via quantitative trait loci (QTL) mapping, microarray analysis and transgenic expression [2, 3]. These genes are mainly involved in stress signalling, transcriptional regulation, and direct response to cold/frost, including cellular membrane stabilization. The gene I nducer of C bf E xpression 2 (Ice2) is a basic helix-loop-helix transcription factor that binds to promoters of the C -repeat B inding F actor (Cbf) gene family and activates their transcription under frost stress in hexaploid wheat . Over-expression of Arabidopsis Ice2  results in increased tolerance to deep freezing stress at a temperature of -20C° after cold acclimation. The Cbf gene family belongs to the family of APETALA2 transcription factors. In barley, diploid and hexaploid wheat several cereal Cbf homologs have been cloned and mapped to the Fr2 locus on homoeologous group 5, which coincides with a major QTL for frost tolerance [6–8]. Using wheat-rye addition lines, Campoli et al.  assigned twelve members of the Cbf gene family to the long arm of chromosome 5R in rye. Several studies in Arabidopsis provide evidence that allelic variation in the Cbf gene family forms the molecular basis for the freezing tolerance QTL [10, 11]. Cbf transcription factors activate Co ld R esponsive (COR) genes through binding to cis-elements in the promoters of COR genes under cold stress in Arabidopsis . More than 70 proteins encoded by COR genes are involved in direct response to cold/frost. Dehydrins, also known as Late Embryogenesis Abundant II (LEA II), are among the proteins that protect other proteins and membranes from cellular damage caused by dehydration . In barley, 13 dehydrin genes (Dhn 1-13) have been identified . Transcripts of Dhn1, Dhn2, Dhn3, Dhn4, Dhn7, and Dhn9 were detected in plants subjected to cold acclimation at 4°C followed by mild frost at -2°C or -4°C . Dhn1 and Dhn3 were mapped in barley to chromosome 5H near a QTL for winter hardiness and on chromosome 6H, respectively . Recent studies showed that cold/frost regulation and vernalization are interconnected [16, 17]. Winter cereals require long exposure to cold in winter, the so-called vernalization, to accelerate flowering in the next spring. This process prevents the early transition of winter cereals into the less cold-tolerant reproductive phase. Vrn1 has been mapped to the second locus conferring frost tolerance, Fr1, on the long arm of homoeologous group 5 near the Fr2 locus . Transcript levels of all cold-induced Cbf genes at the frost tolerance locus Fr-H2 in barley are significantly higher in lines harbouring the vrn1 winter allele than in lines harbouring the Vrn1 spring allele . It remains unknown how the Cbf family members interact with Vrn1 under frost stress.
To unveil genetic diversity among candidate genes involved in the frost response network in rye, one Middle European and four Eastern European populations were studied. Cultivated rye shows a wide range of diversity, reflecting adaptation to various environments and selection pressures . Middle European populations are well-adapted to the more moderate Middle European climate which is in the transition zone between temperate and continental climate, whereas Eastern European populations show good adaptation to a continental climate with severe winters. Thus, differences between Middle and Eastern European populations in allele number and/or frequencies of frost-related candidate genes are expected. Several studies have investigated genome-wide genetic diversity in rye based on molecular markers, including isoenzymes  and simple sequence repeats (SSRs) . None, however, have investigated locus-specific genetic diversity at the gene level.
Linkage disequilibrium (LD), the non-random combination of alleles at different loci, determines the marker density required for marker-based studies, such as association mapping or genomic selection . Studies on the extent of LD in various crops, such as Triticum durum , Zea mays [25, 26], and Sorghum bicolor , indicate large variation in the extent of LD. The effect of germplasm on LD is clearly observed in barley, where LD decays within 0.4 kb in wild material and extends up to 212 kb in elite lines . LD decay can also vary considerably from locus to locus due to different recombination rates and selection pressures at different regions of the genome. In addition, higher levels of LD are observed in self-pollinating species compared to outcrossing species, indicating that mating systems play a role . Since rye is an outcrossing species, a low level of LD with a rapid decay is expected. To the best of our knowledge there is no prior study on the pattern of LD within and between rye genes.
The objectives of this study were to investigate nucleotide and haplotype diversity, the extent and pattern of LD, and population differences among eleven candidate genes (ScCbf2, ScCbf6, ScCbf9b, ScCbf11, ScCbf12, ScCbf14, ScCbf15, ScVrn1, ScIce2, ScDhn1, and ScDhn3) involved in the frost tolerance network in five winter rye populations from Belarus, Germany and Poland.
Plant material and DNA extraction
Plant material was derived from five open-pollinated winter rye breeding populations, four from Eastern Europe, PR 2733 (Belarus), EKOAGRO (Poland), SMH2502 (Poland), ROM103 (Poland), and one from Middle Europe, Petkus (Germany). For convenience, they will be referred to as PR, EKO, SMH, ROM, and Petkus, respectively. The Petkus population has undergone several cycles of recurrent selection, while the breeding history of the four Eastern European populations is unknown. Since rye is an outcrossing species, it is highly heterozygous, which leads to difficulties in determining haplotype phase. To address this problem, gamete capture was performed. Between 15 and 68 heterozygous plants from each of the five populations were crossed with the self-fertile inbred line Lo152 resulting in 201 heterozygous S0 plants, each with one gamete known. The plants were grown in a growth chamber and DNA was extracted from leaves according to Rogowsky et al. .
Candidate gene selection and primer design
Eleven candidate genes, ScCbf2, ScCbf6, ScCbf9b, ScCbf11, ScCbf12, ScCbf14, ScCbf15, ScVrn1, ScIce2, ScDhn1, and ScDhn3, were selected based on their association with frost tolerance in closely related species. Individual Cbf genes were selected based on an expression study in rye  and linkage mapping in barley and diploid wheat [6, 8], Vrn1 based on linkage mapping and a real-time PCR expression study in wheat [18, 31], Ice2 based on an expression study in wheat , and Dhn1 and Dhn3 based on an expression study in barley . We followed the Cbf nomenclature proposed by Skinner et al. , whereby names with the same number followed by different letters describe highly identical but distinct genes, for example, the highly identical Cbf9a and Cbf9b genes first identified by Jaglo et al. . Primers for all genes were designed using Primer-BLAST from the NCBI database (http://www.ncbi.nlm.nih.gov/tools/primer-blast/) based on sequences available in GenBank; information can be found in Additional file 1. Due to limited information on rye DNA sequences in GenBank, primers for ScVrn1, ScIce2, ScDhn1 and ScDhn3 were designed based on homologous genes in H. vulgare, T. aestivum and T. monococcum. Despite lack of homology in non-coding regions, putative functional regions of the candidate genes could be amplified. A 250 bp fragment of the promoter and first exon of ScVrn1 was amplified since there is evidence that this region is one of the determinants of winter/spring growth habit in barley and wheat [34, 35].
Amplification of candidate genes and DNA sequencing
Fourteen fragments of eleven candidate genes were amplified by PCR in 10 μl reaction volumes containing 10 ng DNA, 150 nM of each primer, 1x Taq DNA polymerase reaction buffer, 1.5 or 2.0 mM MgCl2, 0.2 mM of each dNTP, and 0.5 U Taq DNA polymerase. After an initial denaturation at 96°C for 10 min, 35 cycles were conducted at 96°C for 1 min, primer-specific annealing temperatures at 52-66°C for 1 min, 72°C for 1 min, and a final extension step at 72°C for 15 min. Details on candidate gene amplification were described in Additional file 1. The PCR products were purified in 96-well MultiScreen PCR plates (Millipore Corporation, Billerica, MA, USA) and directly sequenced through the QIAGEN sequencing service (QIAGEN, Hilden, Germany). Amplicons of each S0 plant were sequenced with both forward and reverse PCR primers. Sequence data were assembled into contigs and SNPs were detected using the software Variant Reporter™ V1.0 (Applied Biosystems, Foster City, CA, USA). The DNA sequence of Lo152, a homozygous inbred line, was used as the reference sequence, and alleles of this common parent were subtracted from all sequences to determine the haplotype phase. Heterozygous insertion and deletion events were detected manually by checking sequences from both strands. The web-based program Indelligent v1.2 (http://ctap.inhs.uiuc.edu/dmitriev) was used to resolve heterozygous insertion-deletion events (Indels). In case of large Indels, for example, 200 bp in ScCbf2, which Indelligent could not resolve, amplicons from the respective lines were sub-cloned using the TOPO TA Cloning Kit (Invitrogen, Carlsbad, CA, USA). At least five clones were sequenced to resolve heterozygous Indels. Sequences of the Lo152 reference alleles from the eleven candidate genes were submitted to GenBank under accession numbers HQ730763-HQ730773.
The actual numbers of successful PCR amplification of the 201 lines differed from gene to gene ranging from 128 lines (64%) in ScCbf11 to 198 (98%) in ScVrn1. Missing amplification products in individual lines were most likely the result of SNPs/Indels in the primer binding sites. However, absence of some Cbf genes in particular lines, as has recently been reported in barley and wheat [36, 37] cannot be excluded as an alternative explanation.
Sequence polymorphisms were deduced from sequence comparisons in gene-wise sequence alignments. For convenience, polymorphic sites along the sequence were numbered starting with "SNP1". Lo152 alleles were excluded from all analyses. Haplotypes and haplotype frequencies were determined within each candidate gene using DnaSP v5.10  and Arlequin v3.1 , respectively.
Nucleotide diversity (π) was calculated as the average number of nucleotide differences per site between two sequences for both, the complete sequences and restricted to exons, and haplotype diversity (Hd) as the probability that two randomly chosen haplotypes from a given population were different . Analyses of nucleotide and haplotype diversity were performed separately for each population as well as for all populations grouped together using the software DnaSP v5.10. DnaSP v5.10 does not take into account alignment gaps that may lead to underestimated diversity values. Hence, to avoid potential bias, Indels were treated as single polymorphic sites. Average nucleotide diversity (π) over all genes was calculated using concatenated sequences in software TASSEL v2.1 (http://www.maizegenetics.net/).
To test for selection Tajima's D was calculated as the difference between the mean pairwise nucleotide differences (π) and the number of segregating sites (S) relative to their standard error using the software DnaSP v5.10. The statistical significance of Tajima's D was obtained assuming that D follows the beta distribution . The rate ratio of non-synonymous to synonymous substitutions (d N /d S ) was calculated according to the method introduced by Yang and Nielsen  implemented in the program YN00 of software package PAML v4.4c . Significant departure from the standard neutral model, i.e. d N /d S = 1, was assessed by the likelihood ratio test implemented in the CODEML program of PAML v4.4c.
SSR genotyping and genetic diversity analyses
where N was the number of observed genotypes (201 or less), N i the number of genotypes with type i alleles among the N genotypes, n the number of genotypes in each population, and S was the total number of alleles among the N genotypes. To visualize the degree of variation within and between populations, principal co-ordinate analysis (PCoA) was performed using NTSYSpc v2.2 (Applied Biostatistics Inc., Setauket, NY, USA) based on DICE similarity coefficients for SSRs and haplotypes of candidate genes . Analysis of molecular variance (AMOVA)  was performed based on SSRs using Arlequin v3.1  with 15,000 permutations of the data to estimate statistical significance at P < 0.001 for each variance component in Fisher's exact test. The Lo152 alleles were excluded from all analyses.
where Γ = 4Nc and n is the number of lines compared. The LD decay curve was estimated using a non-linear least-squares estimate of Γ fit by the nls function in the R software package, http://www.r-project.org, separately for each population and for all populations pooled together. The approach of Breseghello and Sorrells  was used to determine threshold values of r2 that indicated significant LD. Briefly, r 2 values were estimated from 37 unlinked SSR markers and square root transformed so that they would be better approximated by a Normal distribution. The 95th percentile from the empirical distribution of all pairwise r (n = 666) derived from the 37 unlinked SSR markers was selected as the threshold value, with the rationale that any values above the threshold could in high likelihood be attributable to genetic linkage. Threshold values were calculated separately for each population and for all populations pooled together. The extent of LD was estimated as the point where the LD decay curve passed below the threshold.
DNA sequence polymorphisms
Summary information of candidate gene (CG) sequences: Analyzed fragment length, gene coverage, number of lines, number of SNPs, rate ratio of non-synonymous to synonymous substitutions (d N /d S ), number of Indels and haplotypes, haplotype (Hd) and nucleotide diversity (π), Tajima's D, and linkage disequilibrium (LD)
Fragment length (bp)
Gene coverage a
No. of lines b
No. of SNPs c(non-synonymous)
d N /d S
No. of Indels
No. of haplotypes
π± SD × 10-3(only exon)
Intra-genic LD (r 2 )
0.67 ± 0.02
1.5 ± 0.1 (1.4 ± 0.1)
0.44 ± 0.04
3.6 ± 0.3
0.98 ± 0.03
7.1 ± 0.3 (11.5 ± 0.2)
0.65 ± 0.02
14.5 ± 0. 9
0.89 ± 0.02
8.8 ± 1.0 (7.7 ± 0.1)
0.17 ± 0.04
1.5 ± 0.3
0.68 ± 0.04
3.0 ± 0.2
0.33 ± 0.05
2.7 ± 0.5 (4.4 ± 0.1)
0.73 ± 0.03
8.1 ± 0.6 (8.9 ± 0.1)
0.80 ± 0.02
11.2 ± 0.6 (0)
0.11 ± 0.03
0.4 ± 0.1 (0)
Locus-wise and genome-wide genetic diversity
Nucleotide diversity (π) ranged from 0.4 × 10-3 in ScVrn1 to 14.5 × 10-3 in ScCbf11, and when restricted to exons, from 0 in ScIce2 and ScVrn1 to 14.5 × 10-3 in ScCbf11 (Table 1). The biggest difference between analyses of π for the whole gene compared to restriction to exons occurred in ScIce2 where π decreased from 11.2 to 0 due to absence of SNPs in the exon. Haplotype diversity (Hd) ranged from 0.11 in ScVrn1 to 0.98 in ScCbf9b. A significant positive Tajima's D value was observed over all populations for ScCbf15 and ScIce2, whereas a significant negative value was observed in ScDhn1. Rate ratios of non-synonymous to synonymous substitutions (d N /d S ) were < 1 for ScCbf2, ScCbf6, ScCbf9b, ScCbf11, ScCbf12, ScCbf14, ScDhn1, and ScDhn3. ScCbf15 was the only gene with a d N /d S ratio > 1. d N /d S was significant for ScCbf9b, ScCbf12, ScCbf14, ScCbf15, ScDhn1, and ScDhn3. Due to lack of polymorphisms in their coding sequences d N /d S was not calculated for ScIce2 and ScVrn1.
Genetic diversity within populations based on 37 SSR markers
No. of lines
No. of private alleles a(%)
Average no. of alleles (range)
PIC b± SD
Allelic richness c
0.50 ± 0.02
0.49 ± 0.03
0.46 ± 0.03
0.51 ± 0.01
0.37 ± 0.02
Genetic variation within and between populations
High level of nucleotide and haplotype diversity in rye
We investigated the genetic diversity of five winter rye populations from Middle and Eastern Europe. SNP frequency and nucleotide diversity are affected by several factors, including selection, mutation, mating system, effective population size, and demography . SNP frequency observed in the 5 rye populations under study was on average 1 SNP every 52 bp and the average nucleotide diversity (π) ranged from 0.4 × 10-3 to 14.5 × 10-3 with an average value of π = 5.6 × 10-3. These values are as high as those reported in maize landraces, where one study reported a rate of one SNP per 62 bp, a range of π from 0.1 × 10-3 to 13.3 × 10-3 and an average value of π equal to 4.0 × 10-3 . Some studies have suggested that comparisons among different species should be restricted to homologous genes . Nucleotide diversities of three Cbf homologs (AtCbf1, AtCbf2 and AtCbf3) in 34 Arabidopsis ecotypes ranged from π = 2.6 × 10-3 to 6.9 × 10-3 , a smaller range compared to this study (π = 1.5 × 10-3 to 14.5 × 10-3), which is likely due to the different mating system. In addition, the Cbf gene family in rye encompasses more members than in Arabidopsis, which could result in less selection pressure on individual genes with complementary function in the frost tolerance network and consequently in a higher nucleotide diversity. The buffering effect induced by a large number of duplicated genes leads to a higher variation in individual duplicated genes, a phenomenon also observed in polyploid plants . It is worth re-iterating that inference concerning the nucleotide diversity of ScVrn1 was restrained since only a partial fragment of the gene, 30% of the coding region, could be amplified due to limited available rye sequences for primer design. Observed haplotype diversities of HvCbf9b in Hordeum spontaneum, old cultivars and modern cultivars of H. vulgare were 0.48, 0.18, and 0.06, respectively, which is much lower than that of ScCbf9b in this study (0.98 ± 0.03) .
A reduced genetic diversity was observed in five of the eleven genes. One possible explanation is that directional selection on the loci responsible for fitness related traits such as frost tolerance might reduce diversity within locally adapted populations due to an increase in the frequency of alleles contributing to adaptation . ScCbf15 and ScIce2 showed significant positive values of Tajima's D (2.14 and 2.34, respectively; P < 0.05) over all populations, indicating balancing selection, whereby genotypes carrying alleles with intermediate frequency are favored. Positive Tajima's D values can also be observed if a population was formed from a recent admixture of two different populations, which cannot be excluded in this study. Dhn1 showed a significant negative value of Tajima's D (P < 0.05), indicating purifying selection, whereby an excess of polymorphisms with low frequencies was observed. However, population growth can also result in significant negative values of Tajima's D. Interestingly, Dhn1 in Scots pine has also been described as subject to positive selection , implying that Dhn1 is possibly a target of selection in different species. ScCbf9b, ScCbf12, ScCbf14, ScDhn1, and ScDhn3 had a d N /d S ratio significantly smaller than 1 (P < 0.01 or P < 0.001), whereas ScCbf15 had a d N /d S ratio significantly greater than 1 (P < 0.001). These findings can be interpreted as indication for purifying and positive selection, respectively . However, it was pointed out that inferring selection pressure based on the d N /d S ratio is difficult from within-species data where segregating polymorphisms rather than fixed substitutions are observed [53, 58, 59].
Based on haplotype frequencies in the eleven candidate genes, the single Middle European population did not distinctly differ from the Eastern European ones (Figure 1). Since we have no information on the temporal breeding history of the Eastern European populations, it is beyond the scope of this study to make inferences on the selection pressure due to contrasting winter temperatures in these Middle and Eastern European populations. One possible explanation for a lack of differentiation might be seed exchanges between them. However, little is known about these processes, since pedigrees of the four Eastern European populations were not accessible.
Genetic variation within and among populations
Assessment of genetic diversity based on genome-wide SSRs and locus-specific candidate genes are complementary investigations, the former providing a global view of the rye genome and the latter restricted to genes involved in the frost tolerance network. Genome-wide assessment of diversity using SSR markers revealed a higher genetic diversity for the Eastern European populations PR, EKO, SMH, and ROM compared to the Middle European Petkus population. One reason for this finding might be a bottleneck effect due to a higher selection pressure in the Petkus population, whereby it could be assumed that many "unfavourable" minor alleles were eliminated to pave the way for plants with desirable traits. The Petkus population, one of the major heterotic groups in rye, has systematically been improved by more than 5 cycles of full sib recurrent selection and a reduction in allele diversity of SSR markers due to hitchhiking with linked loci which were targets of selection is probable. The reduction of genetic diversity due to human-induced selection has been well documented in barley and maize [36, 60, 61]. By contrast, the Eastern European populations experienced a lower selection pressure by mass or half sib selection in the breeding programs, where introgression of foreign material was common in order to keep genetic variability on a high level. Interestingly, no reduction of genetic diversity was observed in the Petkus population based on candidate genes. One possible explanation is that at the time where selection took place, winters in Germany, the provenance of the Petkus population, were harsh enough to form a similar selection pressure on the Petkus population compared to Eastern European populations under Eastern European winters. It must be stated however, that Petkus is the only representative for the Middle European rye populations in our study and thus our conclusions on population differences must be limited to the Petkus population. Another reason could be that frost tolerance is a complex quantitative trait involving large gene networks comprising individual genes contributing only small effects, thereby making it difficult to detect selection signatures, such as reduction of genetic diversity in candidate genes. PCoA based on both candidate genes and SSRs showed high genetic variation between individuals within populations and limited clustering of lines from the same population, findings in accordance with previously reported investigations of 26 rye populations based on isoenzyme markers  and 12 rye populations based on RFLP markers . Similar results have also been reported in other outcrossing species, including white clover  and perennial ryegrass , probably a consequence of the obligate cross-pollinated reproductive behaviour of outcrossing species. On the contrary, investigations in the self-pollinated species rice have revealed larger variation between populations .
Rapid decay of linkage disequilibrium in rye
The extent of LD in rye across all eleven candidate genes and over all populations was approximately 520 bp using r 2 = 0.16 as a critical threshold estimated from a separate analysis of 37 unlinked SSR markers. This rapid decay of LD could be expected, because compared to self-pollinated species, cross-pollinated rye has a higher effective recombination rate, which leads to a rapid decay of LD . LD decays rapidly in other cross-pollinated species, including douglas fir, maize and ryegrass [25, 53, 65]. However, in self-pollinated species LD can extend up to 10-30 kb as in Arabidopsis [66, 67] and 212 kb in cultivated barley . Pairwise LD measured by r2 based on SSRs was very low (mean r 2 = 0.01), which was expected since the 37 SSRs have an average marker interval of 21 cM according to the integrated consensus map of Gustafson et al. .
LD results from the interplay of many factors. Selection, which causes locus-specific bottlenecks, is one of the factors that increases LD between selected alleles at linked loci. Homologs of ScCbf (except ScCbf11 in this case) were closely linked and located in the Fr-H2/Fr-A m 2 frost locus spanning approximately 0.8 cM in the genetic maps of barley and diploid wheat on homoeologous group 5 [6, 8, 17]. The order of Cbf genes in the genetic map is consistent in both species . The Cbf gene family is a large regulatory gene family with more than 20 members in barley, diploid and hexaploid wheat [6, 8, 69], sharing a high sequence similarity and induced under frost stress. It has been suggested that the members of the Cbf gene family have slightly different functions in the frost responsive network [8, 9]. In this study, a large variation of mean r 2 in seven Cbf genes (0.13 to 0.92) was observed, indicating that the family has probably undergone diverse selection history. LD can be increased by selection, for instance, by selective sweeps in which the alleles at flanking loci of a locus under selection are rapidly swept to high frequency or fixation . Arabidopsis' AtCbf2 was implicated as subject to selection, resulting in functional divergence from AtCbf1 and AtCbf3 after Cbf gene duplication . In this study an observed strong LD block and low nucleotide diversity in ScCbf14 indicated a selective sweep. Among Cbf family members, TaCbf14 has been mapped to the highest peak of the frost tolerance QTL in hexaploid wheat . Two of the SNPs in HvCbf14 were statistically associated with frost tolerance in a European germplasm collection of spring and winter barley . It remains to be demonstrated that the LD block of ScCbf14 found in this study has an influence on frost tolerance in rye.
Genetic diversity is vital to crop improvement. This study of eleven candidate genes with a putative role in frost response and 37 genome-wide SSRs demonstrated high genetic diversity among five winter rye populations from Middle and Eastern Europe. Most of the diversity was observed within populations. The Middle European Petkus population differed neither in terms of haplotype frequencies nor in nucleotide diversities in eleven candidate genes from the four Eastern European populations. LD within candidate genes decayed rapidly, falling below r 2 = 0.16 within approximately 520 bp. In contrast to selfing species, such as Arabidopsis or barley, low LD in rye promises a higher resolution in genome-wide association mapping. A challenge, however, is that many more markers are required for covering the whole genome. Given the huge genome size of rye, (~8,100 Mb) and until high-density genotyping arrays for rye become available, candidate gene based association mapping remains the most appropriate strategy for gene identification. The SNPs discovered in the promoter or coding regions of the genes investigated in this study, which cause non-synonymous substitutions, are suitable candidates for association mapping and will be studied in more detail with respect to their role in the expression of frost tolerance in rye.
We would like to thank Susanne Schrack and Tobias Dreser for technical assistance and Carmen Berlanas for sequencing ScCbf6 in her master thesis. We acknowledge Andreas Böck and Valentin Wimmer for their help in determining the LD extent. We thank the two anonymous reviewers for their constructive comments. The first author gratefully acknowledges the support of the Graduate School at the Technische Universität München, München, Germany. The project GABI RYE-FROST is funded by the German Federal Ministry of Education and Research (Grant numbers 0315062A and 0315062B).
- Fowler DB, Limin AE: Exploitable genetic variability for cold tolerance in commercially grown cereals. Can J Plant Sci. 1987, 67 (1): 278-278.Google Scholar
- Thomashow MF: Plant cold acclimation: Freezing tolerance genes and regulatory mechanisms. Annu Rev Plant Physiol Plant Mol Biol. 1999, 50: 571-599. 10.1146/annurev.arplant.50.1.571.PubMedView ArticleGoogle Scholar
- Yamaguchi-Shinozaki K, Shinozaki K: Transcriptional regulatory networks in cellular responses and tolerance to dehydration and cold stresses. Annu Rev Plant Biol. 2006, 57: 781-803. 10.1146/annurev.arplant.57.032905.105444.PubMedView ArticleGoogle Scholar
- Badawi M, Reddy YV, Agharbaoui Z, Tominaga Y, Danyluk J, Sarhan F, Houde M: Structure and functional analysis of wheat ICE (inducer of CBF expression) genes. Plant Cell Physiol. 2008, 49 (8): 1237-1249. 10.1093/pcp/pcn100.PubMedView ArticleGoogle Scholar
- Fursova OV, Pogorelko GV, Tarasov VA: Identification of ICE2, a gene involved in cold acclimation which determines freezing tolerance in Arabidopsis thaliana. Gene. 2009, 429 (1-2): 98-103.PubMedView ArticleGoogle Scholar
- Francia E, Barabaschi D, Tondelli A, Laido G, Rizza F, Stanca AM, Busconi M, Fogher C, Stockinger EJ, Pecchioni N: Fine mapping of a HvCBF gene cluster at the frost resistance locus Fr-H2 in barley. Theor Appl Genet. 2007, 115 (8): 1083-1091. 10.1007/s00122-007-0634-x.PubMedView ArticleGoogle Scholar
- Baga M, Chodaparambil SV, Limin AE, Pecar M, Fowler DB, Chibbar RN: Identification of quantitative trait loci and associated candidate genes for low-temperature tolerance in cold-hardy winter wheat. Funct Integr Genomics. 2007, 7 (1): 53-68. 10.1007/s10142-006-0030-7.PubMedView ArticleGoogle Scholar
- Knox AK, Li CX, Vagujfalvi A, Galilba G, Stockinger EJ, Dubcovsky J: Identification of candidate CBF genes for the frost tolerance locus Fr-A(m)2 in Triticum monococcum. Plant Mol Biol. 2008, 67 (3): 257-270. 10.1007/s11103-008-9316-6.PubMedView ArticleGoogle Scholar
- Campoli C, Matus-Cadiz MA, Pozniak CJ, Cattivelli L, Fowler DB: Comparative expression of Cbf genes in the Triticeae under different acclimation induction temperatures. Mol Genet Genomics. 2009, 282 (2): 141-152. 10.1007/s00438-009-0451-9.PubMedPubMed CentralView ArticleGoogle Scholar
- Alonso-Blanco C, Gomez-Mena C, Llorente F, Koornneef M, Salinas J, Martinez-Zapater JM: Genetic and molecular analyses of natural variation indicate CBF2 as a candidate gene for underlying a freezing tolerance quantitative trait locus in Arabidopsis. Plant Physiol. 2005, 139 (3): 1304-1312. 10.1104/pp.105.068510.PubMedPubMed CentralView ArticleGoogle Scholar
- McKhann HI, Gery C, Berard A, Leveque S, Zuther E, Hincha DK, De Mita S, Brunel D, Teoule E: Natural variation in CBF gene sequence, gene expression and freezing tolerance in the Versailles core collection of Arabidopsis thaliana. BMC Plant Biol. 2008, 8: 105-10.1186/1471-2229-8-105.PubMedPubMed CentralView ArticleGoogle Scholar
- Chinnusamy V, Zhu J, Zhu JK: Cold stress regulation of gene expression in plants. Trends Plant Sci. 2007, 12 (10): 444-451. 10.1016/j.tplants.2007.07.002.PubMedView ArticleGoogle Scholar
- Kosova K, Vitamvas P, Prasil IT: The role of dehydrins in plant response to cold. Biol Plant. 2007, 51 (4): 601-617. 10.1007/s10535-007-0133-6.View ArticleGoogle Scholar
- Choi DW, Zhu B, Close TJ: The barley (Hordeum vulgare L.) dehydrin multigene family: sequences, allele types, chromosome assignments, and expression characteristics of 11 Dhn genes of cv Dicktoo. Theor Appl Genet. 1999, 98 (8): 1234-1247. 10.1007/s001220051189.View ArticleGoogle Scholar
- Zhu B, Choi DW, Fenton R, Close TJ: Expression of the barley dehydrin multigene family and the development of freezing tolerance. Mol Gen Genet. 2000, 264 (1-2): 145-153. 10.1007/s004380000299.PubMedView ArticleGoogle Scholar
- Kosova K, Prasil IT, Vitamvas P: The relationship between vernalization-and photoperiodically-regulated genes and the development of frost tolerance in wheat and barley. Biol Plant. 2008, 52 (4): 601-615. 10.1007/s10535-008-0120-6.View ArticleGoogle Scholar
- Galiba G, Vagujfalvi A, Li CX, Soltesz A, Dubcovsky J: Regulatory genes involved in the determination of frost tolerance in temperate cereals. Plant Sci. 2009, 176 (1): 12-19. 10.1016/j.plantsci.2008.09.016.View ArticleGoogle Scholar
- Galiba G, Quarrie SA, Sutka J, Morgounov A, Snape JW: RFLP mapping of the vernalization (Vrn1) and frost resistance (Fr1) genes on chromosome 5A of wheat. Theor Appl Genet. 1995, 90 (7-8): 1174-1179. 10.1007/BF00222940.PubMedView ArticleGoogle Scholar
- Stockinger EJ, Skinner JS, Gardner KG, Francia E, Pecchioni N: Expression levels of barley Cbf genes at the Frost resistance-H2 locus are dependent upon alleles at Fr-H1 and Fr-H2. Plant J. 2007, 51 (2): 308-321. 10.1111/j.1365-313X.2007.0141.x.PubMedView ArticleGoogle Scholar
- Persson K, von Bothmer R, Gullord M, Gunnarsson E: Phenotypic variation and relationships in landraces and improved varieties of rye (Secale cereale L.) from northern Europe. Genet Resour Crop Evol. 2006, 53 (4): 857-866. 10.1007/s10722-004-6694-8.View ArticleGoogle Scholar
- Persson K, Von Bothmer R: Genetic diversity amongst landraces of rye (Secale cereale L.) from northern Europe. Hereditas. 2002, 136 (1): 29-38. 10.1034/j.1601-5223.2002.1360105.x.PubMedView ArticleGoogle Scholar
- Shang HY, Wei YM, Wang XR, Zheng YL: Genetic diversity and phylogenetic relationships in the rye genus Secale L. (rye) based on Secale cereale microsatellite markers. Genet Mol Biol. 2006, 29 (4): 685-691. 10.1590/S1415-47572006000400018.View ArticleGoogle Scholar
- Flint-Garcia SA, Thornsberry JM, Buckler ES: Structure of linkage disequilibrium in plants. Annu Rev Plant Biol. 2003, 54: 357-374. 10.1146/annurev.arplant.54.031902.134907.PubMedView ArticleGoogle Scholar
- Maccaferri M, Sanguineti MC, Noli E, Tuberosa R: Population structure and long-range linkage disequilibrium in a durum wheat elite collection. Mol Breed. 2005, 15 (3): 271-289. 10.1007/s11032-004-7012-z.View ArticleGoogle Scholar
- Tenaillon MI, Sawkins MC, Long AD, Gaut RL, Doebley JF, Gaut BS: Patterns of DNA sequence polymorphism along chromosome 1 of maize (Zea mays ssp mays L.). Proc Natl Acad Sci USA. 2001, 98 (16): 9161-9166. 10.1073/pnas.151244298.PubMedPubMed CentralView ArticleGoogle Scholar
- Ching A, Caldwell KS, Jung M, Dolan M, Smith OS, Tingey S, Morgante M, Rafalski AJ: SNP frequency, haplotype structure and linkage disequilibrium in elite maize inbred lines. BMC Genet. 2002, 3: 19-10.1186/1471-2156-3-19.PubMedPubMed CentralView ArticleGoogle Scholar
- Hamblin MT, Fernandez MGS, Casa AM, Mitchell SE, Paterson AH, Kresovich S: Equilibrium processes cannot explain high levels of short- and medium-range linkage disequilibrium in the domesticated grass Sorghum bicolor. Genetics. 2005, 171 (3): 1247-1256. 10.1534/genetics.105.041566.PubMedPubMed CentralView ArticleGoogle Scholar
- Caldwell KS, Russell J, Langridge P, Powell W: Extreme population-dependent linkage disequilibrium detected in an inbreeding plant species, Hordeum vulgare. Genetics. 2006, 172 (1): 557-567. 10.1534/genetics.104.038489.PubMedPubMed CentralView ArticleGoogle Scholar
- Rogowsky PM, Guidet FLY, Langridge P, Shepherd KW, Koebner RMD: Isolation and characterisation of wheat-rye recombinants involving chromosome arm 1DS of wheat. Theor Appl Genet. 1991, 82 (5): 537-544. 10.1007/BF00226788.PubMedView ArticleGoogle Scholar
- Campoli C, Matus-Cadiz MA, Pozniak CJ, Cattivelli L, Fowler DB: Comparative expression of Cbf genes in the Triticeae under different acclimation induction temperatures. Molecular Genetics and Genomics. 2009, 282 (2): 141-152. 10.1007/s00438-009-0451-9.PubMedPubMed CentralView ArticleGoogle Scholar
- Sutton F, Chen DG, Ge XJ, Kenefick D: Cbf genes of the Fr-A2 allele are differentially regulated between long-term cold acclimated crown tissue of freeze-resistant and - susceptible, winter wheat mutant lines. BMC Plant Biol. 2009, 9: 34-10.1186/1471-2229-9-34.PubMedPubMed CentralView ArticleGoogle Scholar
- Skinner JS, von Zitzewitz J, Szucs P, Marquez-Cedillo L, Filichkin T, Amundsen K, Stockinger EJ, Thomashow MF, Chen THH, Hayes PM: Structural, functional, and phylogenetic characterization of a large CBF gene family in barley. Plant Mol Biol. 2005, 59 (4): 533-551. 10.1007/s11103-005-2498-2.PubMedView ArticleGoogle Scholar
- Jaglo KR, Kleff S, Amundsen KL, Zhang X, Haake V, Zhang JZ, Deits T, Thomashow MF: Components of the Arabidopsis C-repeat/dehydration-responsive element binding factor cold-response pathway are conserved in Brassica napus and other plant species. Plant Physiol. 2001, 127 (3): 910-917. 10.1104/pp.010548.PubMedPubMed CentralView ArticleGoogle Scholar
- Beales J, Laurie DA, Devos KM: Allelic variation at the linked AP1 and PhyC loci in hexaploid wheat is associated but not perfectly correlated with vernalization response. Theor Appl Genet. 2005, 110 (6): 1099-1107. 10.1007/s00122-005-1938-3.PubMedView ArticleGoogle Scholar
- Yan L, Loukoianov A, Tranquilli G, Helguera M, Fahima T, Dubcovsky J: Positional cloning of the wheat vernalization gene VRN1. Proc Natl Acad Sci USA. 2003, 100 (10): 6263-6268. 10.1073/pnas.0937399100.PubMedPubMed CentralView ArticleGoogle Scholar
- Fricano A, Rizza F, Faccioli P, Pagani D, Pavan P, Stella A, Rossini L, Piffanelli P, Cattivelli L: Genetic variants of HvCbf14 are statistically associated with frost tolerance in a European germplasm collection of Hordeum vulgare. Theoretical and Applied Genetics. 2009, 119 (7): 1335-1348. 10.1007/s00122-009-1138-7.PubMedPubMed CentralView ArticleGoogle Scholar
- Knox AK, Dhillon T, Cheng HM, Tondelli A, Pecchioni N, Stockinger EJ: CBF gene copy number variation at Frost Resistance-2 is associated with levels of freezing tolerance in temperate-climate cereals. Theor Appl Genet. 2010, 121 (1): 21-35. 10.1007/s00122-010-1288-7.PubMedView ArticleGoogle Scholar
- Yang ZH: PAML 4: Phylogenetic analysis by maximum likelihood. Mol Biol Evol. 2007, 24 (8): 1586-1591. 10.1093/molbev/msm088.PubMedView ArticleGoogle Scholar
- Khlestkina EK, Ma HMT, Pestsova EG, Roder MS, Malyshev SV, Korzun V, Borner A: Mapping of 99 new microsatellite-derived loci in rye (Secale cereale L.) including 39 expressed sequence tags. Theor Appl Genet. 2004, 109 (4): 725-732. 10.1007/s00122-004-1659-z.PubMedView ArticleGoogle Scholar
- Hackauf B, Wehling P: Identification of microsatellite polymorphisms in an expressed portion of the rye genome. Plant Breed. 2002, 121 (1): 17-25. 10.1046/j.1439-0523.2002.00649.x.View ArticleGoogle Scholar
- Liu KJ, Muse SV: PowerMarker: an integrated analysis environment for genetic marker analysis. Bioinformatics. 2005, 21 (9): 2128-2129. 10.1093/bioinformatics/bti282.PubMedView ArticleGoogle Scholar
- Petit RJ, El Mousadik A, Pons O: Identifying populations for conservation on the basis of genetic markers. Conserv Biol. 1998, 12 (4): 844-855. 10.1046/j.1523-1739.1998.96489.x.View ArticleGoogle Scholar
- Goudet J: FSTAT (Version 1.2): A computer program to calculate F-statistics. J Hered. 1995, 86 (6): 485-486.Google Scholar
- Dice LR: Measure of the amount of ecologic association between species. Ecology. 1945, 26 (3): 297-302. 10.2307/1932409.View ArticleGoogle Scholar
- Excoffier L, Smouse PE, Quattro JM: Analysis of molecular variance inferred from metric distances among DNA haplotypes: application to human mitochondrial DNA restriction data. Genetics. 1992, 131 (2): 479-491.PubMedPubMed CentralGoogle Scholar
- Excoffier L, Laval G, Schneider S: Arlequin (version 3.0): An integrated software package for population genetics data analysis. Evolutionary Bioinformatics. 2005, 1: 47-50.Google Scholar
- Hill WG, Robertson A: Linkage disequilibrium in finite populations. Theor Appl Genet. 1968, 38: 226-231. 10.1007/BF01245622.PubMedView ArticleGoogle Scholar
- Weir B: Genetic data analysis II. Sunderland, MA: Sinauer Associates; 1996.Google Scholar
- Hill WG, Weir BS: Variances and covariance of squared linkage disequilibrium in finite populations. Theor Popul Biol. 1988, 33 (1): 54-78. 10.1016/0040-5809(88)90004-4.PubMedView ArticleGoogle Scholar
- Breseghello F, Sorrells ME: Association mapping of kernel size and milling quality in wheat (Triticum aestivum L.) cultivars. Genetics. 2006, 172 (2): 1165-1177. 10.1534/genetics.105.044586.PubMedPubMed CentralView ArticleGoogle Scholar
- Nei M: Molecular Evolutionary Genetics. NY: Columbia University Press;1987.Google Scholar
- Yamasaki M, Schroeder S, Sanchez-Villeda H, Gaut B, McMullen MD: Empirical analysis of selection screens for domestication and improvement loci in maize by extended DNA sequencing. Plant Genome. 2008, 1 (1): 33-43. 10.3835/plantgenome2007.09.0530.View ArticleGoogle Scholar
- Krutovsky KV, Neale DB: Nucleotide diversity and linkage disequilibrium in cold-hardiness- and wood quality-related candidate genes in Douglas fir. Genetics. 2005, 171 (4): 2029-2041. 10.1534/genetics.105.044420.PubMedPubMed CentralView ArticleGoogle Scholar
- Lin YH, Hwang SY, Hsu PY, Chiang YC, Huang CL, Wang CN, Lin TP: Molecular population genetics and gene expression analysis of duplicated CBF genes of Arabidopsis thaliana. BMC Plant Biol. 2008, 8: 16-10.1186/1471-2229-8-16.View ArticleGoogle Scholar
- Lawton-Rauh A: Evolutionary dynamics of duplicated genes in plants. Mol Phylogenet Evol. 2003, 29 (3): 396-409. 10.1016/j.ympev.2003.07.004.PubMedView ArticleGoogle Scholar
- Biswas S, Akey JM: Genomic insights into positive selection. Trends Genet. 2006, 22 (8): 437-446. 10.1016/j.tig.2006.06.005.PubMedView ArticleGoogle Scholar
- Wachowiak W, Balk PA, Savolainen O: Search for nucleotide diversity patterns of local adaptation in dehydrins and other cold-related candidate genes in Scots pine (Pinus sylvestris L.). Tree Genet Genomes. 2009, 5 (1): 117-132. 10.1007/s11295-008-0188-3.View ArticleGoogle Scholar
- Yang Z: Adaptive molecular evolution. Handbook of Statistical Genetics volume 1. Edited by: Balding DJ, Bishop M, Cannings C. West Sussex: John Wiley & Sons; 2003: 229-250.Google Scholar
- Kryazhimskiy S, Plotkin JB: The Population Genetics of dN/dS. Plos Genet. 2008, 4 (12): e1000304-10.1371/journal.pgen.1000304.PubMedPubMed CentralView ArticleGoogle Scholar
- Kilian B, Ozkan H, Kohl J, von Haeseler A, Barale F, Deusch O, Brandolini A, Yucel C, Martin W, Salamini F: Haplotype structure at seven barley genes: relevance to gene pool bottlenecks, phylogeny of ear type and site of barley domestication. Mol Genet Genomics. 2006, 276 (3): 230-241. 10.1007/s00438-006-0136-6.PubMedView ArticleGoogle Scholar
- Tian F, Stevens NM, Buckler ES: Tracking footprints of maize domestication and evidence for a massive selective sweep on chromosome 10. Proc Natl Acad Sci USA. 2009, 106: 9979-9986. 10.1073/pnas.0901122106.PubMedPubMed CentralView ArticleGoogle Scholar
- George J, Dobrowolski MP, de Jong EV, Cogan NOI, Smith KF, Forster JW: Assessment of genetic diversity in cultivars of white clover (Trifolium repens L.) detected by SSR polymorphisms. Genome. 2006, 49 (8): 919-930. 10.1139/G06-079.PubMedView ArticleGoogle Scholar
- Bolaric S, Barth S, Melchinger AE, Posselt UK: Genetic diversity in European perennial ryegrass cultivars investigated with RAPD markers. Plant Breed. 2005, 124 (2): 161-166. 10.1111/j.1439-0523.2004.01032.x.View ArticleGoogle Scholar
- Zhang DL, Zhang HL, Wang MX, Sun JL, Qi YW, Wang FM, Wei XH, Han LZ, Wang XK, Li ZC: Genetic structure and differentiation of Oryza sativa L. in China revealed by microsatellites. Theor Appl Genet. 2009, 119 (6): 1105-1117. 10.1007/s00122-009-1112-4.PubMedView ArticleGoogle Scholar
- Xing Y, Frei U, Schejbel B, Asp T, Lubberstedt T: Nucleotide diversity and linkage disequilibrium in 11 expressed resistance candidate genes in Lolium perenne. BMC Plant Biol. 2007, 7: 43-10.1186/1471-2229-7-43.PubMedPubMed CentralView ArticleGoogle Scholar
- Kim S: Recombination and linkage disequilibrium in Arabidopsis thaliana. Nat Genet. 2007, 39: 1151-1155. 10.1038/ng2115.PubMedView ArticleGoogle Scholar
- Nordborg M, Hu TT, Ishino Y, Jhaveri J, Toomajian C, Zheng HG, Bakker E, Calabrese P, Gladstone J, Goyal R, Jakobsson M, Kim S, Morozov Y, Padhukasahasram B, Plagnol V, Rosenberg NA, Shah C, Wall JD, Wang J, Zhao KY, Kalbfleisch T, Schulz V, Kreitman M, Bergelson J: The pattern of polymorphism in Arabidopsis thaliana. PLoS Biol. 2005, 3 (7): 1289-1299. 10.1371/journal.pbio.0030196.View ArticleGoogle Scholar
- Gustafson JP, Ma XF, Korzun V, Snape JW: A consensus map of rye integrating mapping data from five mapping populations. Theor Appl Genet. 2009, 118 (4): 793-800. 10.1007/s00122-008-0939-4.PubMedView ArticleGoogle Scholar
- Badawi M, Danyluk J, Boucho B, Houde M, Sarhan F: The CBF gene family in hexaploid wheat and its relationship to the phylogenetic complexity of cereal CBFs. Mol Genet Genomics. 2007, 277 (5): 533-554. 10.1007/s00438-006-0206-9.PubMedPubMed CentralView ArticleGoogle Scholar
- Kraakman ATW, Niks RE, Van den Berg P, Stam P, Van Eeuwijk FA: Linkage disequilibrium mapping of yield and yield stability in modern spring barley cultivars. Genetics. 2004, 168 (1): 435-446. 10.1534/genetics.104.026831.PubMedPubMed CentralView ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.