Skip to main content
  • Research article
  • Open access
  • Published:

DNA polymorphisms and haplotype patterns of transcription factors involved in barley endosperm development are associated with key agronomic traits

Abstract

Background

Association mapping is receiving considerable attention in plant genetics for its potential to fine map quantitative trait loci (QTL), validate candidate genes, and identify alleles of interest. In the present study association mapping in barley (Hordeum vulgare L.) is investigated by associating DNA polymorphisms with variation in grain quality traits, plant height, and flowering time to gain further understanding of gene functions involved in the control of these traits. We focused on the four loci BLZ1, BLZ2, BPBF and HvGAMYB that play a role in the regulation of B-hordein expression, the major fraction of the barley storage protein. The association was tested in a collection of 224 spring barley accessions using a two-stage mixed model approach.

Results

Within the sequenced fragments of four candidate genes we observed different levels of nucleotide diversity. The effect of selection on the candidate genes was tested by Tajima's D which revealed significant values for BLZ1, BLZ2, and BPBF in the subset of two-rowed barleys. Pair-wise LD estimates between the detected SNPs within each candidate gene revealed different intra-genic linkage patterns. On the basis of a more extensive examination of genomic regions surrounding the four candidate genes we found a sharp decrease of LD (r2<0.2 within 1 cM) in all but one flanking regions.

Significant marker-trait associations between SNP sites within BLZ1 and flowering time, BPBF and crude protein content and BPBF and starch content were detected. Most haplotypes occurred at frequencies <0.05 and therefore were rejected from the association analysis. Based on haplotype information, BPBF was associated to crude protein content and starch content, BLZ2 showed association to thousand-grain weight and BLZ1 was found to be associated with flowering time and plant height.

Conclusions

Differences in nucleotide diversity and LD pattern within the candidate genes BLZ1, BLZ2, BPBF, and HvGAMYB reflect the impact of selection on the nucleotide sequence of the four candidate loci.

Despite significant associations, the analysed candidate genes only explained a minor part of the total genetic variation although they are known to be important factors influencing the expression of seed quality traits. Therefore, we assume that grain quality as well as plant height and flowering time are influenced by many factors each contributing a small part to the expression of the phenotype. A genome-wide association analysis could provide a more comprehensive picture of loci involved in the regulation of grain quality, thousand grain weight and the other agronomic traits that were analyzed in this study. However, despite available high-throughput genotyping arrays the marker density along the barely genome is still insufficient to cover all associations in a whole genome scan. Therefore, the candidate gene-based approach will further play an important role in barley association studies.

Background

Association mapping is receiving considerable attention in plant genetics for its potential to fine map quantitative trait loci (QTL), validate candidate genes, and identify alleles of interest. Association mapping has several advantages over linkage mapping: First, a potentially larger number of alleles per locus can be surveyed simultaneously [1]. Second, results refer to a more representative genetic background. Third, the resolution of association mapping is increased because all recombination events accumulated in the population history are taken into consideration [2]. There are two ways to identify DNA-markers for QTL via association mapping: whole genome association mapping and re-sequencing of candidate genes. In whole genome association mapping populations are genotyped with a genome-wide set of closely linked and evenly distributed markers. This essentially requires a large number of markers and is therefore expensive and statistically complex [3]. The number of markers to be employed depends on the genome size and the extent of LD along the chromosomes. In a candidate gene-based approach, genotyping is targeted to functional and positional candidate genes for the trait under consideration [4]. This approach is assisted by (i) plant genomics resources such as expressed sequence tag (EST) databases, (ii) available knowledge on gene function in model organisms, and (iii) referenced information on physiology, biochemistry, and molecular genetics available for the trait of interest. In the present study we applied a candidate gene-based approach to find marker-trait associations for agronomic important traits in a spring barley collection.

The improvements of grain yield and quality, either for food or for feed, are paramount targets in any barley breeding program. It is known that transcription factors play an important role in controlling expression during seed development. Genetic differences in the synthesis of storage proteins can already be observed at the transcriptional level [5–7]. In barley, B-hordein represents the largest fraction of the storage protein. Functional analysis of the promoters of genes specifically expressed in the cereal endosperm, such as those encoding B-hordein (e.g. Hor2), has demonstrated the existence of cis-acting motifs capable of interacting with nuclear proteins that are putatively responsible for their tissue specificity and temporal regulation [8–10]. The endosperm box is a conserved cis-acting element, which contains two distinct protein binding sites: the prolamin-box (PB) and the GCN4-like motif (GLM). Four transcription factors (TFs) are the gibberellin-regulated Myb factor (GAMYB), the barley leucine zippers 1 and 2 (BLZ1, BLZ2), and the barley prolamin box binding factor (BPBF) that were shown to be involved in the transcription of B-hordeins encoded by the Hor2 locus.

BLZ1 mRNA is detected during early endosperm development. The single copy gene is a transcriptional activator that interacts with endosperm-specific gene promoters (Figure 1). Vicente-Carbajosa et al. [11] demonstrated the involvement of BLZ1 in the regulation of hordein gene expression through binding to the GLM. BLZ1 protein functions as a transcriptional activator and is able to form either homodimers or heterodimers with BLZ2 [12]. The BLZ2 mRNA expression is restricted to the endosperm and its protein specifically binds to the GLM [12]. As indicated by its designation, the BPBF has been shown to activate hordein genes through binding to the PB [13, 14]. Transient expression experiments in developing barley endosperms demonstrate that BPBF trans-activates transcription from the PB element of a native Hor2 promoter [14]. Positive regulatory interaction was observed between BPBF and HvGAMYB in the control of endosperm gene expression during seed development [13]. In developing seeds abundant expression of the transcription factor HvGAMYB is induced by gibberellic acid. Its mRNA can be detected in the starchy endosperm and other grain tissues [13]. The protein trans-actives transcription from the native Hor2 promoter through binding to a third motif (5'-AACA/TA-3') that is present in endosperm-specific genes. Thus, HvGAMYB represents a key regulator of genes specifically expressed in the endosperm during seed development [13]. In addition to seed tissue, HvGAMYB also plays a role in other aspects of plant growth and development [15] and BLZ1 expression was also detected in leaves and roots [11].

Figure 1
figure 1

Interplay between the candidate genes and the promoter region of a target gene (e.g. Hor2 ). Transcription start is displayed as ATG. Arrows to both sides show known interactions between the four transcription factors BLZ1, BLZ2, BPBF and HvGAMYB. Grey boxes indicate cis-regulatory motifs named as mentioned in the black boxes.

A phenotypically well characterized spring barley collection was recently established by Haseneyer et al. [16] as resource for this association study. Information about morphological properties of the accessions is available and population structure was determined with 45 EST-derived SSR markers. In the current paper we report on the analysis of nucleotide diversity parameters for the above mentioned candidate genes. Results are presented on the association between sequence polymorphisms within BLZ1 (chromosome 5 H), BLZ2 (chromosome 1 H), BPBF (chromosome 5 H), and HvGAMYB (chromosome 3 H) and the phenotypic variation of the five agronomic traits thousand-grain weight, starch content, protein content, plant height and flowering time.

Methods

Plant material and phenotypic analyses

The above mentioned collection of spring barleys selected from the Barley Core Collection (BCC) and the Federal ex situ Genebank (HOR) was used in this study (Additional file 1). The germplasm set consists of 128 two-rowed and 96 six-rowed accessions originating from Europe (N = 109), East Asia (N = 40), America (N = 30), and West Asia and North Africa (N = 45). Eighteen accessions were classified as "breeding/research material", 55 accessions as landraces/traditional cultivars while the remaining accessions represent advanced breeding lines and cultivars. Accessions were phenotypically evaluated at Stuttgart-Hohenheim (South Germany), Irlbach (South Germany) and Bergen-Wohlde (North Germany) in 2004 and 2005. Each trial was arranged in microplots in a 25 × 15 lattice design with three replicates. Thousand-grain weight (TGW), flowering time (FT), and plant height (PH) were recorded. Grain quality (crude protein content (CPC) and starch content (STR)) was assessed by near infrared reflectance spectroscopy (NIRS, for further details see [16]).

Population structure

All 224 accessions were genotyped with 45 simple sequence repeat (SSR) markers that are evenly distributed across the barley genome [17]. A population structure with K = 2 subgroups was inferred from the SSR data by using the STRUCTURE 2.0 software package [18, 19]. The individual steps of analysis were described in detail by Haseneyer et al. [16].

Genotyping and genetic mapping

Eight seeds from each accession were grown in the greenhouse and leaves from 2-week-old seedlings were harvested and bulked for genomic DNA extraction using the method described in Stein et al. [20]. PCR-primers were designed using the software Primer3 [21]. Primer sequences and the fragment-specific PCR profile conditions are given in additional file 2. PCR for single nucleotide polymorphism (SNP) analysis by DNA sequencing was performed as described in full detail by [22]. In preparation for DNA sequencing, we purified the PCR amplicons in 384-well plates and adjusted to similar molarity. 10 ng PCR product was used as template for cycle sequencing. DNA sequences were determined using ABI BigDye Terminator 3.1 chemistry and 96-capillary sequencer systems (ABI 3730 × l). Forward and reverse PCR primers were used as sequencing primers (Additional file 2). DNA sequence ladders were processed for quality scoring using a software package based on the poly-phred system [23]. We applied the program Sequencher™ Version 4.5 (Gene Codes Cooperation) for sequence alignment and editing. All positions given in the text correspond to the positions in the haplotype sequence alignments related to the start codon (Additional file 3).

BLZ1, BLZ2 and BPBF were genetically mapped in the Oregon Wolfe Barley (OWB) mapping population developed by Costa et al. [24]. Positions were determined on an updated OWB map [25]. Therefore, we designed cleaved amplified polymorphic sequence (CAPS) markers that require the use of the restriction enzymes Nci I (BLZ1, SNP 1733), Ssp I (BLZ2, SNP 2161), and Sty I (BPBF, SNP -210). HvGAMYB was mapped earlier by Haseneyer et al. [26].

Diversity and association analysis

The candidate genes' DNA fragments were sequenced for each accession of the collection. DnaSP Version 4.10 [27] was applied for the statistical sequence analysis. This software does not take into account the alignment gaps that may lead to underestimated diversity values. To avoid potential bias, insertion-deletion events (indels) were treated as single sites. Nucleotide diversity estimated as Pi (Ï€) [28], haplotype diversity (Hd), and Tajima's D [29] were computed. Diversity values of gene fragments showing no sequence overlap were calculated fragment-wise and then the arithmetic average was computed.

LD between pairs of polymorphic sites (minor allele frequency, MAF ≥ 0.05) was estimated by TASSEL software, version 1.9.3 [30]. LD is expressed by r2 [31] and the statistical significance (P-value) of the observed LD is estimated by Monte-Carlo approximation of Fisher's exact test [32], with 1,000 permutations. In order to estimate the local decay of LD, additional markers flanking the candidate genes at increasing distances were investigated in the entire collection. The expected value of r2 is E(r2) = 1/(1+C), where C = 4 Nc, N is the effective population size, and c is the recombination fraction between sites [33]. This model was employed in nonlinear regression of r2 on c, treating N as a parameter to be estimated, using PROC NLIN of the SAS System for Windows (Version 9.1.3.)

Combined analyses of phenotypic and genotypic data were performed using Version 9.1.3 of the SAS System for Windows. We followed a two-stage mixed model approach [34, 35] where in the first stage adjusted entry means and weights were computed for each trial, which were then subjected to a mixed model analysis combined over trials in the second stage. Our analysis is based on the assumption that genotypes are a random sample from the world collection of barley genotypes. In order to compute adjusted means for single trials, however, we formally took genotypes as fixed in the first stage, fitting a linear model with fixed effects for genotypes and replicate and random effects for block and error. Thus, adjusted means were unbiased estimates of the genotypes' performances in the different environments, which allowed formulating a mixed model for adjusted means in the second stage. Note that taking genotypes random, and hence computing best linear unbiased predictors (BLUPs) of genotype performances, in the first stage would have caused biases that would have been difficult to account for in stage two [35]. In the second stage, the following model terms were fitted: overall mean (fixed), trial main effects (fixed), genotype main effect (random), genotype-by-trial interaction (random). In addition, spike morphology and geographic origin were modelled by fixed effects for 'row number', which had two levels, and 'origin', which had four levels. Population structure was modelled by fixed-effects regression on a Q matrix of membership probabilities of N genotypes in each of K subgroups. The Q matrix was computed using the Bayesian approach of Pritchard et al. [19]. Associations of haplotypes and SNP markers were tested by adding a haplotype or SNP marker covariate to the fixed part of the model. Tests of fixed effects were based on variance estimates using the restricted maximum likelihood (REML) method and denominator degrees of freedom approximated by the method of Kenward and Roger [36]. The genetic variance explained by a fixed effect was computed by the relative reduction in genetic variance when the fixed term was added. Weights to model the error variance of adjusted means in stage two were computed based on the diagonal elements of the inverse of the asymptotic variance-covariance matrix of adjusted means [35]. All variance components were estimated by the REML method. Adjusted means were compared by Wald t-tests [37]. As the haplotype means were not variance balanced, we used the method of Piepho [38] to generate a letter display showing the significance of comparisons. Type I error rate was controlled by the Bonferroni-Holm procedure [39].

Results

Sequence diversity and haplotype analysis

The polymorphism density ranged from 1 polymorphism/31 bp (BLZ2), 1 polymorphism/42 bp (BPBF), 1 polymorphism/55 bp (BLZ1) to 1 polymorphism/74 bp (HvGAMYB). Nucleotide diversities (π) were determined for BLZ1 (1,113 bp), BLZ2 (2,232 bp), BPBF (1,119 bp) and HvGAMYB (3,337 bp) for the entire germplasm set and the geographical and morphological subsets individually (Table 1). Diversity estimates for the entire collection ranged from π = 2.4 × 10-3 (HvGAMYB) to π = 8.1 × 10-3 (BPBF). In most cases individual subgroups showed a similar range of nucleotide diversities for all candidate genes. An exception was only noted for BLZ2 where the two-rowed subset displayed a highly reduced π-value, whereas a high diversity was observed for the East Asian accessions.

Table 1 Estimates of nucleotide and haplotype diversity for the candidate genes BLZ1, BLZ2, BPBF and HvGAMYB

Haplotype analysis indicated a similar diversity at most gene loci and for all subpopulations, although the number of haplotypes per locus ranged from 8 (BLZ1) to 21 (BPBF). The haplotype diversity at BLZ2, BPBF and HvGAMYB was mainly caused by the six-rowed accessions that were particularly frequent in the American, East Asian and West Asian and North African subsets. The two-rowed subset, that primarily included European genotypes, revealed the lowest estimates for all loci considered, especially for the BLZ2 gene.

Linkage disequilibrium

The pairwise LD values revealed different patterns for the genes studied (Figure 2). BLZ1 and HvGAMYB showed strong LD (r2>0.8, P < 0.0001) only between a few polymorphic sites. At the BLZ1 locus two blocks of polymorphism (positions 1740 to 1890 and 2520 to 2774) displayed significant LD estimates higher than r2 = 0.5 (P < 0.0001). BLZ2 and BPBF showed significant LD across the entire sequence. Even beyond the gap of 482 bp between the two sequenced fragments of BPBF (positions -368 to 62 and 579 to 1129) LD persisted at a high level (r2>0.4, P < 0.0001). The sites 2316 and 2361 at the BLZ2 gene and 870 at the BPBF locus segregated separately from the remaining polymorphic sites.

Figure 2
figure 2

Linkage disequilibrium between the polymorphic sites (MAF>0.05) within the candidate loci BLZ1, BLZ2, BPBF, HvGAMYB. Asterisk indicates transcription start, dashed lines indicate regions that were sequenced and "i" and "e" column indicates polymorphisms in introns and exons, respectively. MAF = minor allele frequency.

The results of the extended LD study of markers flanking the four candidate genes showed that LD remained significant at distances up to 19 cM. However, individual r2 values sharply decreased to r2<0.1 within 1 cM in the surrounding regions of all four candidate genes (Figure 3). Only in the proximal region of BLZ2 sustained levels of LD were observed up to 10 cM (Additional file 4).

Figure 3
figure 3

LD decay plot in the surrounding regions of the four candidate genes as a function of genetic distance (in cM). Dots indicate pairwise comparisons between SNP alleles with minor allele frequency larger 0.05. The curve shows nonlinear regression of r2 on genetic distance.

The impact of selection on the four candidate genes was tested by calculating Tajima's D. Significant deviations from the mutation-drift equilibrium were observed for BLZ1 and BPBF for the entire collection (Table 2). Within the two-rowed subset BLZ1, BLZ2, and BPBF were significant, while in the six-rowed subset only BPBF revealed a significant Tajima D-value. No significant values were observed for HvGAMYB.

Table 2 Tajima's D for the candidate genes BLZ1, BLZ2, BPBF, and HvGAMYB

Marker-trait association

For all association analyses the model including population structure (two subgroups referred to as 'K2'), 'row number' and 'origin' was applied. Several SNPs within the candidate gene BLZ1 being in high LD with one another were significantly associated with flowering time (Table 3, Additional file 5). They explained between 6.5 to 7.5% of the genetic variation and phenotypic means of the respective SNP alleles revealed a significant two-day difference in flowering time.

Table 3 Percentage explained variance (%Var), phenotypic means of SNP alleles, and significant (P < 0.05) marker-trait associations

Within the BPBF gene polymorphisms in the exonic and 5'- and 3'-untranslated regions were associated with crude protein content. Phenotypic means of the SNP alleles revealed a significant difference of 5.9% in crude protein content. One SNP (pos. 579) explained 12.4% of the genetic variation. Estimates for pairwise LD were significant for these sites with r2>0.5 (P < 0.0001). A portion (5'-untranslated region, and SNP at pos. 62) of these sites showed association to starch content revealing a significant difference between phenotypic means of the contrasting alleles.

Haplotype-trait association

In accordance with the handling of SNP data, haplotype-trait associations were restricted to haplotype classes which were more frequent than 0.05. Applying this frequency threshold, three haplotype classes were detected for BLZ1, BLZ2, and BPBF and four haplotype classes for HvGAMYB (Additional file 5) that were entered in the association analysis. BLZ1 was significantly (P-value < 0.05) associated with flowering time and plant height, and explained 3.3% and 3.1% of the genetic variation, respectively (Table 4). A weak association of BLZ1 with crude protein content was observed explaining 2.7% of the genetic variation. BLZ2 haplotypes were associated with thousand-grain weight and explained 4.0% of the genetic variance (Table 4). Haplotypes of the candidate gene BPBF were significantly associated to crude protein content and starch content and explained 8.2% and 6.0% of the genetic variation, respectively.

Table 4 Haplotype-trait associations (P = 0.05) and percentage explained genetic variance (%Var) of the candidate genes' haplotypes

Discussion

In this study, a worldwide collection of spring barley accessions was used to perform marker-trait association analyses based on a set of four candidate genes for grain quality.

Different patterns of sequence diversity, haplotype diversity and LD were observed for the candidate genes BLZ1, BLZ2, BPBF, and HvGAMYB. A similar variability of LD patterns was found for different members of the CBF (C-repeat binding factor) transcription factor family [40]. In the present study LD within genes was weak for BLZ1 and HvGAMYB but strong for the other two genes. The high number of sequence polymorphisms detected at the BLZ2 locus is in accordance with observations on the homologous gene Opaque 2 in maize [41]. Compared to the remaining members of the bzip class of regulatory genes, BLZ2 and its homologues seem to be characterized by exceptionally high levels of polymorphism. The high SNP frequency in BLZ2 is not reflected in a high nucleotide or haplotype diversity since diversity in this gene is caused by only few frequent and many rare SNPs. HvGAMYB showed the lowest and BPBF the highest values of nucleotide diversity, whereas the opposite was found for the haplotype diversity. This pattern is due to the high pairwise LD at the BPBF locus resulting in few frequent and many rare haplotypes. The low level of LD, which was observed at the HvGAMYB locus, might be due to a low selection pressure on this gene during its domestication and breeding history [26].

Malting barley is characterized by a low protein and high starch content [42]. In this regard, two-rowed barley is preferred by European brewing industry due to the favourable protein to starch relation. A strong selection for these two negatively correlated traits might have had a bearing on nucleotide diversity in the underlying candidate genes. This is apparent for the BLZ2 locus where the reduced diversity in the European subset corresponds with a high proportion of two-rowed genotypes in this geographic subset. The observed reduction in sequence variation might be a consequence of purifying selection [43]. The negative Tajima D value might indicate such kind of selection for BLZ2 in the two-rowed subset caused by the elimination of deleterious alleles and leaving only one major haplotype which is common to 95 of the 108 two-rowed accessions.

It is well known that selection in autogamous organisms leads to an increase in LD [44]. In this context, selection may affect the regulatory regions of genes, or target regulatory loci rather than the protein-coding region of genes [45]. In Zea mays L. the ear underwent dramatic morphological alteration upon domestication and has been a continuing target of selection for grain yield [46]. Therefore, Hufford et al. [46] hypothesize that genes targeted by selection are more likely to be expressed in tissues that experienced high levels of morphological divergence during crop improvement. One such tissue in barley is the endosperm since its characteristics are the determinants of malting quality [47]. Since expression of BLZ2 and BPBF is restricted to the endosperm [12, 14] the selection and corresponding enrichment of only a few favourable alleles at these loci entails an increase in LD. Determining the nucleotide diversity of these two genes in wild barley would allow verification of this hypothesis.

The tentative appraisal about the impact of selection on the four candidate genes was investigated by calculating Tajima's D. A significant deviation from the mutation-drift-equilibrium, especially in the two-rowed subgroup, was observed for the three candidate genes that were found to be associated to the target traits. In Europe, two-rowed barley is the main target for the improvement of seed quality parameters. This is in accordance with the significant Tajima D values obtained for the three loci in this subgroup indicating footprints of selection on BLZ1, BLZ2 and BPBF. However, selection might act in different ways: In case of BLZ2 selection resulted in the accumulation of a large number of low frequency SNP alleles as 61% of the recorded SNPs have a MAF < 5%. In conjunction with the extended LD across this gene, this results in the presence of only one major haplotype for this gene which is present in 54% of the accessions. Within the subset of two-rowed barleys, this haplotype is even more dominant showing a frequency of 88% (see previous pragraph). In case of BLZ1 and BPBF, 11% and 23% of the SNPs show a MAF < 5%. Hence, selection was effective in the elimination of rare SNP alleles and the accumulation of moderate frequent SNP alleles was promoted. The indication that these two genes are targeted by balancing selection is supported by significant Tajima D values.

The detected marker-trait associations, even for polymorphisms explaining only a minor portion of the trait variation, are attributed to the high statistical power achieved by (i) extensive and precise phenotyping of the target traits as reflected by high heritability estimates [16], (ii) considering the population structure of the collection and (iii) the high phenotypic variability of the worldwide collection and the large nucleotide diversity within the selected candidate genes. However, the power to detect an association also depends on the number of accessions in the individual haplotype classes on which the analysis is based. In the analysed collection the high degree of diversity resulted in prevalence of rare haplotypes that occurred in less than 5% of accessions and thus were excluded from the analysis to avoid spurious associations. Interestingly, most of the phenotypic differences were found between those rare haplotype classes. Hence, a considerably larger collection size or the selective enrichment of haplotype classes would be needed to warrant a proper sample size for rare haplotypes as well.

The observed haplotype associations of BLZ1 with flowering time and plant height corroborate the hypothesis of Vicente-Carbajosa et al. [11] that this gene is involved in developmental processes and photoperiodic response. Pleiotropic effects of a single gene as observed for BLZ1 lead to overlapping QTL position estimates for different traits providing a basis for enhancing the effectiveness of marker-assisted selection [48]. Thus, candidate gene-based association studies for two or more traits might substantially contribute to cultivar improvement. However, in the present study, we could not identify an advantageous haplotype or SNP sites in the investigated candidate genes comparable to the ones found in the sh4-d gene in rice, the Q-gene in wheat and the ppd-H1 gene in barley [49–51]. As the present candidate genes were described as trans-active regulators for hordein encoding genes [11, 12, 14, 52], we hypothesize that they influence both grain protein composition and protein content and thus are of importance not only for malting [53] but also for nutritional quality [54].

Both, marker-trait and haplotype-trait associations yielded comparable results. In both approaches significant associations of BLZ1 with flowering time and BPBF with crude protein and starch content were found. Using haplotypes instead of SNP alleles revealed a higher number of associations. This shows the higher sensitivity and statistical power of haplotype-trait associations [55, 56] as here accessions are divided in several classes whereas in marker-trait association only two classes, representing the two SNP alleles, are considered. The portion of explained genetic variance by SNP sites was in reasonable agreement with the explained genetic variance by haplotypes. As would be expected for a quantitative trait, only a small part of the entire genetic variation could be explained by the variation occurring at the candidate loci. It follows that the remaining variation is due to additional loci that also influence the expression of crude protein content, starch content, thousand-grain weight, plant height, and flowering time.

With the increasing availability of high-throughput genotyping platforms for barley (DArT array [57], oligonucleotide pool assay [58]), estimation of genome-wide LD decay and whole genome association studies become a feasible alternative to the analysis of candidate genes. LD studies based on such genotyping data that were retrieved for a collection of genotypes resulted in a decay of intrachromosomal LD below r2<0.2 within 2.6 cM [59], r2<0.15 within 3.2 cM [60] and r2<0.5 within 3.9 cM [58], respectively. Complementary to the decrease in genetic diversity, LD has been shown to increase from wild barley via landraces to modern cultivars [58, 61]. Notwithstanding this observation, LD within cultivated barley is also population dependent so that comparison of genome-wide LD between collections composed of accessions with different origins is difficult. In our world-wide collection the extent of genome-wide LD decreases more rapidly than in geographically restricted collections of domesticated barley germplasm [58–60]. The chromosomal regions surrounding the four candidate genes display a rapid LD decay. However, genome wide DNA fingerprinting of the present population would significantly increase the knowledge about LD structure in the present collection and facilitate comparisons to other mapping panels regarding local LD patterns and trait associations.

Conclusions

Nucleotide diversity and LD patterns of BLZ1, BLZ2, BPBF, and HvGAMYB revealed differences between the candidate genes and between geographical and morphological subsets of the collection. This reflects the impact of selection on the nucleotide sequence of these four candidate loci.

According to literature, the four candidate genes represent transcriptional key regulators in barley. However, only three of the four selected candidate genes could be confirmed by haplotype-trait association studies. We conclude that there is still an incomplete knowledge about the expression and interaction of genes controlling the quantitative traits crude protein content, starch content, thousand-grain weight, plant height, and flowering time in barley. Additionally, both haplotypes and SNPs only explained a part of the genetic variation. Therefore, and in accordance with their quantitative inheritance, we assume that the investigated seed traits, plant height, and flowering time are influenced by many additional hitherto unknown factors each contributing a small part to the expression of the phenotype.

Although genome-wide association mapping could provide a more comprehensive picture of loci involved in the regulation of crude protein content, starch content, thousand-grain weight, flowering time, and plant height there is a risk of overlooking an association in genome-wide association studies. As has been demonstrated in the present study, a gene may contain SNPs that are associated and others that are not associated with the trait under consideration. If only one or two SNPs per locus (e.g. EST) would be interrogated as is presently the case with many SNP marker arrays used for whole genome scans, it is possible that the "right" SNP was not included in the array. On the other hand, a candidate gene-based approach might suffer from the limited knowledge about candidates for a given trait and hence only a part of the genetic variation for this trait is captured. Further verification of the observed associations is difficult owing to the quantitative nature of the target trait. Moreover, LD decay and hence genetic resolution of the present population is still insufficient to preclude that the observed association is not due to the presence of a physically linked gene being in LD with the candidate gene. Notwithstanding this fact, future candidate gene-based approaches will greatly benefit from the continuous accumulation of knowledge on gene function and regulation. Because of this and due to the still insufficient marker coverage of the barley genome, the candidate gene-based association mapping will continue to play an important role in barley.

References

  1. Flint-Garcia SA, Thuillet AC, Yu J, Pressoir G, Romero SM, Mitchell SE, Doebley J, Kresovich S, Goodman MM, Buckler ES: Maize association population: a high-resolution platform for quantitative trait locus dissection. Plant J. 2005, 44: 1054-1064. 10.1111/j.1365-313X.2005.02591.x.

    Article  PubMed  CAS  Google Scholar 

  2. Buckler ES, Thornsberry JM: Plant molecular diversity and applications to genomics. Curr Opin Plant Biol. 2002, 5: 107-111. 10.1016/S1369-5266(02)00238-8.

    Article  PubMed  CAS  Google Scholar 

  3. Hirschhorn JN, Daly MJ: Genome-wide association studies for common diseases and complex traits. Nat Rev Genet. 2005, 6: 95-108. 10.1038/nrg1521.

    Article  PubMed  CAS  Google Scholar 

  4. Pflieger S, Lefebvre V, Causse M: The candidate gene approach in plant genetics: a review. Mol Breed. 2001, 7: 275-291. 10.1023/A:1011605013259.

    Article  CAS  Google Scholar 

  5. Bartels D, Thompson RD: Synthesis of messenger-RNAs coding for abundant endosperm proteins during wheat-grain development. Plant Sci. 1986, 46: 117-125. 10.1016/0168-9452(86)90118-4.

    Article  CAS  Google Scholar 

  6. Soave C, Salamini F: Organization and regulation of zein genes in maize endosperm. Philos T Roy Soc B. 1984, 304: 341-347. 10.1098/rstb.1984.0029.

    Article  CAS  Google Scholar 

  7. Sorensen MB, Cameronmills V, Brandt A: Transcriptional and post-transcriptional regulation of gene-expression in developing barley endosperm. Mol Gen Genet. 1989, 217: 195-201. 10.1007/BF02464881.

    Article  Google Scholar 

  8. Hammond-Kosack MC, Holdsworth MJ, Bevan MW: In vivo footprinting of a low molecular weight glutenin gene (LMWG-1 D1) in wheat endosperm. Embo J. 1993, 12: 545-554.

    PubMed  CAS  PubMed Central  Google Scholar 

  9. Vicente-Carbajosa J, Beritashvili DR, Kraev AS, Skryabin KG: Conserved structure and organization of B hordein genes in the Hor 2 locus of barley. Plant Mol Biol. 1992, 18: 453-458. 10.1007/BF00040661.

    Article  PubMed  CAS  Google Scholar 

  10. Zhao Y, Leisy DJ, Okita TW: Tissue-specific expression and temporal regulation of the rice glutelin Gt3 gene are conferred by at least two spatially separated cis-regulatory elements. Plant Mol Biol. 1994, 25: 429-436. 10.1007/BF00043871.

    Article  PubMed  CAS  Google Scholar 

  11. Vicente-Carbajosa J, Onate L, Lara P, Diaz I, Carbonero P: Barley BLZ1: a bZIP transcriptional activator that interacts with endosperm-specific gene promoters. Plant J. 1998, 13: 629-640. 10.1111/j.1365-313X.1998.00068.x.

    Article  PubMed  CAS  Google Scholar 

  12. Onate L, Vicente-Carbajosa J, Lara P, Diaz I, Carbonero P: Barley BLZ2, a seed-specific bZIP protein that interacts with BLZ1 in vivo and activates transcription from the GCN4-like motif of B-hordein promoters in barley endosperm. J Biol Chem. 1999, 274: 9175-9182. 10.1074/jbc.274.14.9175.

    Article  PubMed  CAS  Google Scholar 

  13. Diaz I, Vicente-Carbajosa J, Abraham Z, Martinez M, Isabel-La Moneda I, Carbonero P: The GAMYB protein from barley interacts with the DOF transcription factor BPBF and activates endosperm-specific genes during seed development. Plant J. 2002, 29: 453-464. 10.1046/j.0960-7412.2001.01230.x.

    Article  PubMed  CAS  Google Scholar 

  14. Mena M, Vicente-Carbajosa J, Schmidt RJ, Carbonero P: An endosperm-specific DOF protein from barley, highly conserved in wheat, binds to and activates transcription from the prolamin-box of a native B-hordein promoter in barley endosperm. Plant J. 1998, 16: 53-62. 10.1046/j.1365-313x.1998.00275.x.

    Article  PubMed  CAS  Google Scholar 

  15. Woodger FJ, Gubler F, Pogson BJ, Jacobsen JV: A Mak-like kinase is a repressor of GAMYB in barley aleurone. Plant J. 2003, 33: 707-717. 10.1046/j.1365-313X.2003.01663.x.

    Article  PubMed  CAS  Google Scholar 

  16. Haseneyer G, Stracke S, Paul C, Einfeldt C, Broda A, Piepho H-P, Graner A, Geiger HH: Population structure and phenotypic variation of a spring barley world collection set up for association studies. Plant Breeding. 2009

    Google Scholar 

  17. Thiel T, Michalek W, Varshney RK, Graner A: Exploiting EST databases for the development and characterization of gene-derived SSR-markers in barley (Hordeum vulgare L.). Theor Applied Genet. 2003, 106: 411-422.

    CAS  Google Scholar 

  18. Falush D, Stephens M, Pritchard JK: Inference of population structure using multilocus genotype data: linked loci and correlated allele frequencies. Genetics. 2003, 164: 1567-1587.

    PubMed  CAS  PubMed Central  Google Scholar 

  19. Pritchard JK, Stephens M, Donnelly P: Inference of population structure using multilocus genotype data. Genetics. 2000, 155: 945-959.

    PubMed  CAS  PubMed Central  Google Scholar 

  20. Stein N, Herren G, Keller B: A new DNA extraction method for high-throughput marker analysis in a large-genome species such as Triticum aestivum. Plant Breeding. 2001, 120: 354-356. 10.1046/j.1439-0523.2001.00615.x.

    Article  CAS  Google Scholar 

  21. Rozen S, Skaletsky HJ: Primer3 on the WWW for general users and for biologist programmers. Bioinformatics Methods and Protocols: Methods in Molecular Biology. Edited by: Krawetz S, Misener S. 2002, Humana Press Totowa, NJ, 365-386.

    Google Scholar 

  22. Sauer S, Reinhardt R, Lehrach H, Gut IG: Single-nucleotide polymorphisms: analysis by mass spectrometry. Nat Protoc. 2006, 1: 1761-1771. 10.1038/nprot.2006.257.

    Article  PubMed  CAS  Google Scholar 

  23. Nickerson DA, Tobe VO, Taylor SL: PolyPhred: automating the detection and genotyping of single nucleotide substitutions using fluorescence-based resequencing. Nucleic Acids Res. 1997, 25: 2745-2751. 10.1093/nar/25.14.2745.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  24. Costa JM, Corey A, Hayes PM, Jobet C, Kleinhofs A, Kopisch-Obusch A, Kramer SF, Kudrna D, Li M, Riera-Lizarazu O, Sato K, Szücs P, Toojinda T, Vales MI, Wolfe RI: Molecular mapping of the Oregon Wolfe Barleys: a phenotypically polymorphic doubled-haploid population. Theor Appl Genet. 2001, 103: 415-424. 10.1007/s001220100622.

    Article  CAS  Google Scholar 

  25. Stein N, Prasad M, Scholz U, Thiel T, Zhang H, Wolf M, Kota R, Varshney RK, Perovic D, Grosse I, Graner A: A 1,000-loci transcript map of the barley genome: new anchoring points for integrative grass genomics. Theor Appl Genet. 2007, 114: 823-839. 10.1007/s00122-006-0480-2.

    Article  PubMed  CAS  Google Scholar 

  26. Haseneyer G, Ravel C, Dardevet M, Balfourier F, Sourdille P, Charmet G, Brunel D, Sauer S, Geiger HH, Graner A, Stracke S: High level of conservation between genes coding for the GAMYB transcription factor in barley (Hordeum vulgare L.) and bread wheat (Triticum aestivum L.) collections. Theor Appl Genet. 2008, 117: 321-331. 10.1007/s00122-008-0777-4.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  27. Rozas J, Sanchez-DelBarrio JC, Messeguer X, Rozas R: DnaSP, DNA polymorphism analyses by the coalescent and other methods. Bioinformatics. 2003, 19: 2496-2497. 10.1093/bioinformatics/btg359.

    Article  PubMed  CAS  Google Scholar 

  28. Nei M: Molecular evolutionary genetics. 1987, Columbia University Press

    Google Scholar 

  29. Tajima F: Statistical method for testing the neutral mutation hypothesis by DNA polymorphism. Genetics. 1989, 123: 585-595.

    PubMed  CAS  PubMed Central  Google Scholar 

  30. MaizeGenetics.net - Website of the Edward Buckler lab. [http://www.maizegenetics.net/]

  31. Hill WG, Robertson A: Linkage disequilibrium in finite populations. Theor Appl Genet. 1968, 38: 226-231. 10.1007/BF01245622.

    Article  PubMed  CAS  Google Scholar 

  32. Weir BS: Genetic data analysis II. 1996, Massachusetts, Sinauer

    Google Scholar 

  33. Sved JA: Linkage disequilibrium and homozygosity of chromosome segments in finite populations. Theor Popul Biol. 1971, 2: 125-141. 10.1016/0040-5809(71)90011-6.

    Article  PubMed  CAS  Google Scholar 

  34. Piepho HP: Stability analysis using the SAS system. Agron J. 1999, 91: 154-160. 10.2134/agronj1999.00021962009100010024x.

    Article  Google Scholar 

  35. Smith AB, Cullis BR, Gilmour AR: Analysing variety by environment data using multiplicative mixed models and adjustments for spatial field trend. Biometrics. 2001, 57: 1138-1147. 10.1111/j.0006-341X.2001.01138.x.

    Article  PubMed  CAS  Google Scholar 

  36. Kenward MG, Roger JH: Small sample inference for fixed effects from restricted maximum likelihood. Biometrics. 1997, 53: 983-997. 10.2307/2533558.

    Article  PubMed  CAS  Google Scholar 

  37. Wald A: Tests of statistical hypotheses concerning several parameters when the number of observations is large. T Am Math Soc. 1943, 54: 426-482. 10.2307/1990256.

    Article  Google Scholar 

  38. Piepho HP: algorithm for a letter-based representation of all-pairwise comparisons. J Comput Graph Stat. 2004, 13: 456-466. 10.1198/1061860043515.

    Article  Google Scholar 

  39. Holm S: A simple sequentially rejective Bonferroni test procedure. Scand J Stat. 1979, 6: 65-70.

    Google Scholar 

  40. 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. Theor Appl Genet. 2009, 119: 1335-48. 10.1007/s00122-009-1138-7.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  41. Henry AM, Manicacci D, Falque M, Damerval C: Molecular evolution of the Opaque-2 gene in Zea mays L. J Mol Evol. 2005, 61: 551-558. 10.1007/s00239-005-0003-9.

    Article  PubMed  CAS  Google Scholar 

  42. Eurequa Poland - Principles of barley breeding. [http://www.eurequa.pl/eng/I.2.htm]

  43. Charlesworth B, Morgan M, Charlesworth D: The effect of deleterious mutations on neutral molecular variation. Genetics. 1993, 134: 1289-1303.

    PubMed  CAS  PubMed Central  Google Scholar 

  44. Gupta PK, Rustgi S, Kulwal PL: Linkage disequilibrium and association studies in higher plants: present status and future prospects. Plant Mol Biol. 2005, 57: 461-485. 10.1007/s11103-005-0257-z.

    Article  PubMed  CAS  Google Scholar 

  45. Zeder MA, Emshwiller E, Smith BD, Bradley DG: Documenting domestication: the intersection of genetics and archaeology. Trends Genet. 2006, 22: 139-155. 10.1016/j.tig.2006.01.007.

    Article  PubMed  CAS  Google Scholar 

  46. Hufford KM, Canaran P, Ware DH, McMullen MD, Gaut BS: Patterns of selection and tissue-specific expression among maize domestication and crop improvement loci. Plant Physiol. 2007, 144: 1642-1653. 10.1104/pp.107.098988.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  47. Brennan CS, Harris N, Smith D, Shewry PR: Structural differences in the mature endosperms of good and poor malting barley cultivars. J Cereal Sci. 1996, 24: 171-177. 10.1006/jcrs.1996.0050.

    Article  Google Scholar 

  48. Han F, Romagosa I, Ullrich SE, Jones BL, Hayes PM, Wesenberg DM: Molecular marker-assisted selection for malting quality traits in barley. Mol Breed. 1997, 3: 427-437. 10.1023/A:1009608312385.

    Article  CAS  Google Scholar 

  49. Li C, Zhou A, Sang T: Rice domestication by reducing shattering. Science. 2006, 311: 1936-1939. 10.1126/science.1123604.

    Article  PubMed  CAS  Google Scholar 

  50. Simons KJ, Fellers JP, Trick HN, Zhang Z, Tai YS, Gill BS, Faris JD: Molecular characterization of the major wheat domestication gene Q. Genetics. 2006, 172: 547-555. 10.1534/genetics.105.044727.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  51. Turner A, Beales J, Faure S, Dunford RP, Laurie DA: The pseudo-response regulator Ppd-H1 provides adaptation to photoperiod in barley. Science. 2005, 310: 1031-1034. 10.1126/science.1117619.

    Article  PubMed  CAS  Google Scholar 

  52. Gubler F, Watts RJ, Kalla R, Jacobsen JV: GAMyb: A transcription factor mediating gibberellin-regulated gene expression in aleurone cells of barley. Plant Physiol. 1997, 114: 1493-1493.

    Google Scholar 

  53. Howard KA, Gayler KR, Eagles HA, Halloran GM: The relationship between D hordein and malting quality in barley. J Cereal Sci. 1996, 24: 47-53. 10.1006/jcrs.1996.0036.

    Article  CAS  Google Scholar 

  54. Shewry PR, Tatham AS, Halford NG, Barker JH, Hannappel U, Gallois P, Thomas M, Kreis M: Opportunities for manipulating the seed protein composition of wheat and barley in order to improve quality. Transgenic Res. 1994, 3: 3-12. 10.1007/BF01976021.

    Article  PubMed  CAS  Google Scholar 

  55. Cardon LR, Abecasis GR: Using haplotype blocks to map human complex trait loci. Trends Genet. 2003, 19: 135-140. 10.1016/S0168-9525(03)00022-2.

    Article  PubMed  CAS  Google Scholar 

  56. Fallin D, Cohen A, Essioux L, Chumakov I, Blumenfeld M, Cohen D, Schork NJ: Genetic analysis of case/control data using estimated haplotype frequencies: application to APOE locus variation and Alzheimer's disease. Genome Res. 2001, 11: 143-151. 10.1101/gr.148401.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  57. Wenzl P, Carling J, Kudrna D, Jaccoud D, Huttner E, Kleinhofs A, Kilian A: Diversity arrays technology (DArT) for whole-genome profiling of barley. Proc Natl Acad Sci USA. 2004, 101: 9915-9920. 10.1073/pnas.0401076101.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  58. Rostoks N, Ramsay L, MacKenzie K, Cardle L, Bhat PR, Roose ML, Svensson JT, Stein N, Varshney RK, Marshall DF, Graner A, Close TJ, Waugh R: Recent history of artificial outcrossing facilitates whole-genome association mapping in elite inbred crop varieties. Proc Natl Acad Sci USA. 2006, 103: 18656-18661. 10.1073/pnas.0606133103.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  59. Zhang LY, Marchand S, Tinker NA, Belzile F: Population structure and linkage disequilibrium in barley assessed by DArT markers. Theor Appl Genet. 2009, 119: 43-52. 10.1007/s00122-009-1015-4.

    Article  PubMed  CAS  Google Scholar 

  60. Comadran J, Thomas WTB, van Eeuwijk FA, Ceccarelli S, Grando S, Stanca AM, Pecchioni N, Akar T, Al-Yassin A, Benbelkacem A, Ouabbou H, Bort J, Romagosa I, Hackett CA, Russell JR: Patterns of genetic diversity and linkage disequilibrium in a highly structured Hordeum vulgare association-mapping population for the Mediterranean basin. Theor Appl Genet. 2009, 119: 175-187. 10.1007/s00122-009-1027-0.

    Article  PubMed  CAS  Google Scholar 

  61. Morrell PL, Toleno DM, Lundy KE, Clegg MT: Low levels of linkage disequilibrium in wild barley (Hordeum vulgare ssp. spontaneum) despite high rates of self-fertilization. Proc Natl Acad Sci USA. 2005, 102: 2442-2447. 10.1073/pnas.0409804102.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

Download references

Acknowledgements

The authors thank Dr. Christian Paul and Merle Alex (Julius-Kuehn-Institute Braunschweig-Voelkenrode, Germany) for laboratory and advisory assistance in the determination of crude protein and starch content by NIRS. We acknowledge Raj Kishore Pasam and Dr. Benjamin Kilian (IPK Gatersleben, Germany) for his help and support in the extended LD study. This work was financed by the German Ministry of Education and Research (BMBF, PTJ-BIO/0313098).

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Andreas Graner.

Additional information

Authors' contributions

GH carried out the molecular genetic studies, the sequence alignment and analyses, the statistical association analyses, and drafted the manuscript. SSt participated in the design and coordination of the study. HPP developed the concept for the statistical analysis. SSa carried out sequencing of the candidate genes. HHG and AG participated in the design and coordination of the study, interpretation of the data and the development of the manuscript. All authors read and approved the final manuscript.

Electronic supplementary material

12870_2009_499_MOESM1_ESM.XLS

Additional file 1: Accessions under study. Information about origin, row number, biological status and haplotypes observed for the candidate genes BLZ1, BLZ2, BPBF and HvGAMYB are given. (XLS 59 KB)

12870_2009_499_MOESM2_ESM.PDF

Additional file 2: Primer sequences for PCR and sequencing of the candidate genes, PCR conditions and fragment range. 1: numbers indicate positions in the nucleotide sequence alignment of the candidate genes' haplotypes given in additional file 3. (PDF 12 KB)

12870_2009_499_MOESM3_ESM.PDF

Additional file 3: Nucleotide sequence alignments of the candidate gene fragments. Description: Haplotype sequences of BLZ1 (reference = [GenBank:X80068.1]), BLZ2 (reference = [GenBank:Y10834.1]), BPBF (reference = [GenBank:AJ000991.1]) and HvGAMYB (reference = [GenBank:AY008692.1]). The alignment position is relative to the ATG and gaps are counted. Abbreviations: hpt = haplotype, cds = coding sequence, gene = genomic sequence (if available). (PDF 175 KB)

12870_2009_499_MOESM4_ESM.PDF

Additional file 4: Linkage disequilibrium in the surrounding region of the candidate genes BLZ1 (A), BLZ2 (B), BPBF (C), and HvGAMYB (D). The position 0.0 cM refers to the candidate gene. The symbols × and ◆ indicate significant (P = 0.05) and non significant pairwise comparisons, respectively. (PDF 48 KB)

12870_2009_499_MOESM5_ESM.XLS

Additional file 5: Haplotype sequence and marker-trait associations detected in the four candidate genes BLZ1 , BLZ2 , BPBF , HvGAMYB. Significant associations (P = 0.05) are indicated by 'x'. The traits crude protein content (CPC), starch content (STR), thousand-grain weight (TGW), plant height (PH), and flowering time (FT) were considered. Haplotype frequencies (in %) and minor allele frequencies (MAF, in %) are given. (XLS 97 KB)

Authors’ original submitted files for images

Rights and permissions

Open Access This article is published under license to BioMed Central Ltd. This is an Open Access article is distributed under the terms of the Creative Commons Attribution License ( https://creativecommons.org/licenses/by/2.0 ), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and permissions

About this article

Cite this article

Haseneyer, G., Stracke, S., Piepho, HP. et al. DNA polymorphisms and haplotype patterns of transcription factors involved in barley endosperm development are associated with key agronomic traits. BMC Plant Biol 10, 5 (2010). https://doi.org/10.1186/1471-2229-10-5

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1471-2229-10-5

Keywords