Single Nucleotide Polymorphisms in HSP17.8 and Their Association with Agronomic Traits in Barley

Abstract

Small heat shock protein 17.8 (HSP17.8) is produced abundantly in plant cells under heat and other stress conditions and may play an important role in plant tolerance to stress environments. However, HSP17.8 may be differentially expressed in different accessions of a crop species exposed to identical stress conditions. The ability of different genotypes to adapt to various stress conditions resides in their genetic diversity. Allelic variations are the most common forms of genetic variation in natural populations. In this study, single nucleotide polymorphisms (SNPs) of the HSP17.8 gene were investigated across 210 barley accessions collected from 30 countries using EcoTILLING technology. Eleven SNPs including 10 from the coding region of HSP17.8 were detected, which form nine distinguishable haplotypes in the barley collection. Among the 10 SNPs in the coding region, six are missense mutations and four are synonymous nucleotide changes. Five of the six missense changes are predicted to be deleterious to HSP17.8 function. The accessions from Middle East Asia showed the higher nucleotide diversity of HSP17.8 than those from other regions and wild barley (H. spontaneum) accessions exhibited greater diversity than the cultivated barley (H. vulgare) accessions. Four SNPs in HSP17.8 were found associated with at least one of the agronomic traits evaluated except for spike length, namely number of grains per spike, thousand kernel weight, plant height, flag leaf area and leaf color. The association between SNP and these agronomic traits may provide new insight for study of the gene's potential contribution to drought tolerance of barley.

About the Authors:

Yanshi Xia

Contributed equally to this work with: Yanshi Xia, Ronghua Li

Affiliations International Crop Research Center for Stress Resistance, College of Life Sciences, Guangzhou University, Guangzhou, China, College of Light Industry and Food Science, South China University of Technology, Guangzhou, China

Ronghua Li

Contributed equally to this work with: Yanshi Xia, Ronghua Li

Affiliation: International Crop Research Center for Stress Resistance, College of Life Sciences, Guangzhou University, Guangzhou, China

Zhengxiang Ning

Affiliation: College of Light Industry and Food Science, South China University of Technology, Guangzhou, China

Guihua Bai

Affiliation: Hard Winter Wheat Genetics Research Unit, United States Department of Agriculture - Agricultural Research Service, Manhattan, Kansas, United States of America

Kadambot H. M. Siddique

Affiliation: The Institute of Agriculture, The University of Western Australia, Crawley, Perth, Australia

Guijun Yan

Affiliation: The Institute of Agriculture, The University of Western Australia, Crawley, Perth, Australia

Michael Baum

Affiliation: International Center for Agricultural Research in the Dry Areas, Aleppo, Syria

Rajeev K. Varshney

Affiliation: International Crops Research Institute for the Semi-Arid Tropics, Patancheru, Greater Hyderabad, India

Peiguo Guo

* E-mail: [email protected]

Affiliation: International Crop Research Center for Stress Resistance, College of Life Sciences, Guangzhou University, Guangzhou, China

Introduction

Small heat shock proteins (sHSPs) with a molecular weight of 15 to 42 kDa are produced ubiquitously in prokaryotic and eukaryotic cells under heat and other stress conditions at various growth stages [1]–[5]. Most sHSPs have strong cytoprotective effects by keeping functional conformations of proteins or refolding denatured proteins, relieving protein from aggregation and removing harmful polypeptides under stress conditions [6]. These sHSPs are classified into several subgroups based on DNA sequence similarity, immunological cross-reactivity and intracellular localization [7].

Heat shock protein 17.8 (HSP17.8), a member of class I cytosolic sHSPs, presents as a dimer under normal physiological conditions and is converted to high oligomeric complexes, ranging from 240 kDa to >480 kDa, after heat shock [3]. Several studies have postulated that HSP17.8 plays a role in tolerance to heat [8], freezing [9], [10] and drought [11], [12]. HSP17.8 showed chaperone-like activity to protect citrate synthase from thermal aggregation at 43°C in the cyanobacterium Anabaena sp. PCC 7120 [8], and played a role in membrane protein targeting to the chloroplast outer membrane in Arabidopsis [3]. At the freezing temperature, HSP17.8 showed at least 1.5 times more expression in freezing-tolerant transgenic maize than in non-transgenic maize [10]. Similarly, HSP17.8 transcripts were present at a higher level in drought-tolerant control plants of Katya cultivar than in non-tolerant Sadovo cultivar of wheat [11]. In the case of barley, HSP17.8 showed expression exclusively in drought-tolerant barley genotypes e.g. Martin and Hordeum spontaneum 41-1, and not in drought-sensitive genotype e.g. Moroc9-75, under drought stress [12].

The ability of different accessions of a crop species to adapt to stress is a result of sequence variation in genes in the accessions of the given crop species. Single nucleotide polymorphisms (SNPs) and small insertions and deletions (indels) are the most common forms of genetic variation in natural populations, and is a reflection of evolution and adaptation that play a prominent role in the heritability of phenotypes [13]. EcoTILLING is the application of TILLING (Targeting Induced Local Lesions IN Genome) for discovery of SNPs and small indels in natural populations. Initially EcoTILLING was used to detect nucleotide variation in five genes in a natural population of Arabidopsis [14]. In the case of poplar (Populus trichocarpa), EcoTILLING identified 63 SNPs in nine different genes in 41 different populations [15]. In brief, several reports have indicated that EcoTILLING is an accurate, high-throughput and low-cost technique for the discovery and evaluation of nucleotide variations in candidate genes associated with target traits [16]–[18]. Recently, association analysis emerged as a powerful approach to identify the role of genetic polymorphism in phenotypic variations in response to environmental stresses [19], [20]. For instance, by using association analysis, one SNP of the HvCBF4 gene was significantly (P<0.001) associated with salt tolerance in 188 Tibetan barley accessions [21], and five SNPs in Lhcb1 gene were significantly (P<0.01) associated with several agronomic traits in 292 barley accessions [22]. However, no effort has been made so far to assay allelic variation in HSP17.8 and its association with agronomic traits in barley.

In this study, EcoTILLING approach was used to detect genetic variations of HSP17.8 in a natural barley population comprising of 210 accessions collected from diverse geographical origins. Distribution of SNPs across different geographic origins (Africa, Middle East Asia, North East Asia, Arabian Peninsula and Europe), row types (two-row or six-row) and other possible categories (such as wild versus landrace) were investigated. Population parameters were estimated using SNPs found in different barley populations. An attempt was also made to assess potential effect of SNPs on protein function and their association with six agronomic traits in barley.

Materials and Methods

Plant Material

Seeds of 210 barley (Hordeum vulgare L.) accessions containing 171 H. vulgare landraces and 39 wild relatives of Hordeum vulgare ssp. spontaneum (hereafter named H. spontaneum) were obtained from the International Centre for Agricultural Research in the Dry Areas (Table S1). Of the 210 accessions, 164 originated from 19 countries in Asia, 40 accessions from six African countries, and 6 accessions from five European countries. For the purpose of comparison, accessions were divided into five adjacent geographic regions as described in Varshney et al. [23]. The geographical distribution of investigated accessions is given in Table 1.

[Figure omitted. See PDF.]

Table 1. Summary of the geographic origins of barley accessions used for allele mining of HSP17.8.

https://doi.org/10.1371/journal.pone.0056816.t001

Phenotypic data

The germplasm collection (210 accessions) was sown with three replications during two growing seasons (2009/2010 and 2010/2011) at the Experimental Station of Guangzhou University, Guangzhou, Guangdong Province, China (23° 16' N; 113°23'E, elevation 16 m asl). Eleven seeds of each accession were planted 30 cm apart in 1.5 m long, single-row plots. All accessions were evaluated for six agronomic traits—flag leaf area (FLA in cm2), spike length (SL in cm), number of grains per spike (NGS), leaf color (LC in SPAD), plant height (PH in cm) and thousand kernel weight (TKW in g)—in replicated field experiments. Three randomly selected plants of each accession from each plot were characterized for the six traits using the following methods:

1. Flag leaf area (FLA; cm2) of the uppermost, fully expanded leaf of the main tiller at flowering. FLA = leaf length×leaf width×0.75 [24].

2. Spike length (SL; cm) was measured from the base of each main spike to the top of the spike, excluding awns, at maturity.

3. Number of grains per spike (NGS) was counted at maturity and the average number of seeds from three spikes used for analysis.

4. Leaf color (LC; SPAD) was determined at heading, before any symptoms of senescence were visible. The middle section of a randomly-selected flag leaf was evaluated between 09:00 and 12:00 h using a chlorophyll meter (SPAD-502, Minolta, Japan) and averaged over two measurements for analysis.

5. Plant height (PH; cm) was measured from ground level to the base of spike at maturity.

6. Thousand kernel weight (TKW; g) was calculated based on a sample of 250 seeds per plot.

Designing of primers

For amplification of a target fragment of HSP17.8 as described in Wienholds et al. [25], gene-specific primers were designed according to the published mRNA sequence of HSP17.8 from GeneBank (accession no. AK368988.1) with melting temperatures around 60°C by using Primer 5.0 software (Premier Biosoft International, Palo Alto, CA, USA) (Table S2). The region covered by the primer pair is 600 base pairs, which includes the whole open reading frame (483 bp) of HSP17.8. Sequence of universal primers M13 was added to 5' ends of HSPforward and HSPreverse as adaptors (Table S2). M13 forward primers labeled with IRDye800 at 5′-end and M13 reverse primers labeled with IRDye700 at 5′-end were synthesized by LI-COR Inc.

DNA extraction and PCR

Genomic DNA of barley accessions was extracted from 200 mg young leaf tissue following Guo et al. [26]. DNA from all samples was quantified using a spectrophotometer and normalized to a concentration of 20 ng/ µl. DNA from each accession was mixed in a 1∶1 ratio with reference DNA (ICARDA IG: 26727) to generate heteroduplexes for point mutation detection.

For nucleotide polymorphism screening with EcoTILLING, the target region of HSP17.8 was amplified by a nested PCR described by Wang et al. [27]. The first PCR was performed with a gene-specific primer pair. After the first PCR reaction, samples were diluted in 90 µl water; 1 µl diluted PCR was used as a template for the second nested PCR reaction. This reaction contained a mixture of four primers: 0.08 µM IRD800-labeled M13 forward primers, 0.02 µM M13HSPforward primers, 0.06 µM IRD700-labeled M13 reverse primers and 0.04 µM M13HSPreverse primers. After the nested PCR, heteroduplexes formation was performed by incubating the reaction mix at 99°C for 10 min, and slow renaturation with 70 cycles of 20 sec at 70°C with a decrement of 0.3°C per cycle.

The PCR products were digested by CEL I as described by Raghavan et al. [28]. Celery juice extract (CEL I) was produced following Guo and Li [29]. For digestion, 10 µl of solution containing 10 mM HEPES (pH 7.5), 10 mM MgSO4, 0.002% (w/v) Triton X-100, 0.2 µg/ml of bovine serum albumin, and 0.4 µl CEL I enzyme solution, was added to 10 µl of heteroduplex DNAs, and incubated at 45°C for 15 min. Digestion was stopped by adding 5 µl 0.25 M EDTA (pH 8.0), mixing thoroughly and then incubating on ice. Digested products were purified using isopropanol as described by Raghavan et al. [28]. After samples purification, 5 µl formamide loading dye and 10 µl ddH2O were added and heated at 85°C for concentrating sample to about 3 µl. Samples were loaded (0.8 µl/well) on denaturing 6.5% polyacrylamide gels on LI-COR 4300 DNA Analyzer. Two electronic image files were produced per gel run, one in the IRD700 channel and the other in IRD800 channel. Tiff images were manually scored using the GelBuddy program [30]. The appearance of cleavage products in both channels at reciprocal size that add up to the full length PCR product, were considered a polymorphic site (Fig. 1). Data summary reports generated by GelBuddy were imported to Microsoft Excel for further analysis. Samples were grouped into putative haplotype categories based on the cleaved banding pattern in evaluated gel-frames.

[Figure omitted. See PDF.]

Figure 1. Detection of polymorphisms for a targeted region of the HSP17.8 gene by EcoTILLING.

Sampled images of the IRD 800 and IRD 700 channels are shown at left and right, respectively. The specific cleavage products appear as intense dark bands between 127 to 512 bp with molecular weights listed to the left in each channel image by arrows. Complementary fragments in corresponding lanes between the IRD 700 and IRD 800 channel images labeled with the same box pattern (including rectangle, oval and diamond). The sizes of complementary fragments in the IRD 700 labeled and the IRD 800 labeled add up to the size of PCR fragment (639 bp). Several intense dark bands near the bottom of the gel in both channels result from random mispriming. Molecular weights are provided by the GelBuddy program. The size of the DNA ladder is listed to the right of the IRD700 image.

https://doi.org/10.1371/journal.pone.0056816.g001

DNA sequencing and statistical analysis

One representative genotype for each unique haplotype was reamplified by gene-specific primers using 40 ng of genomic DNA. The resulted PCR fragment was directly sequenced from both directions using an ABI 3730xl DNA Analyzer by a commercial company (Sangon Biotech Co., Ltd., China) to confirm the polymorphisms. Each polymorphic site was sequenced in more than one accession to confirm only two alleles segregating at any specific site. Sequences were analyzed using ClustalW software (http://www.ebi.ac.uk/tools). The SIFT (Sorting Intolerant from Tolerant) and PARSESNP (Position-Specific Scoring Matrix) programs were used to predict the impact of missense mutations on protein function [31], [32]. Population genetics parameters, including nucleotide diversity (π), haplotype diversity (HD) and Tajima's D [33], were analyzed using DnaSP v5.0 [34].

Association analysis between SNPs and agronomic traits

Association between markers and traits was evaluated using a General Linear Model (GLM_Q) in the TASSEL v3.0 software (http://www.maizegenetics.net/tassel), where the SNP being tested was considered as a fixed effect, and the factor and matrix of subpopulation membership (Q matrix) were used as cofactors to account for population structure. Possible population structure was depicted using genotypic data of the 210 barley accessions and 21 genome-wide SSR markers (3 SSRs for each chromosome) (Table S3) by Structure software version 2 [35]. Three groups were identified (unpublished). Permutations of 1,000 runs were performed to calculate the significant p value for F-test. The association between a marker and a trait is represented by its R2 value, an estimate of the percentage of variance explained by the marker.

Results

Variability of phenotypic traits

Three developmental traits—flag leaf area (FLA in cm2), leaf color (LC in SPAD) and plant height (PH in cm) were measured for all 210 accessions (Table 2). Three yield-related traits—spike length (SL in cm), number of grains per spike (NGS) and 1000-kernel weight (TKW in g) were measured for only 192 and 195 barley accessions in the 2009/2010 and 2010/2011 growing seasons, because 18 and 15 barley accessions did not head, respectively. Large phenotypic variation was observed for all traits, and significant correlations between the various phenotypic traits were found among the barley accessions (Table 2). Two-row and 6-row barley differed significantly in all phenotypic measurements except for TKW. In addition, H. spontaneum accessions had significantly lower scores in FLA, NGS, LC and TKW than H. vulgare landraces.

[Figure omitted. See PDF.]

Table 2. Phenotypic scores from field trials.

https://doi.org/10.1371/journal.pone.0056816.t002

Nucleotide polymorphisms

By using EcoTILLING, 13 polymorphic sites were detected in the targeted region of HSP17.8 across the 210 accessions (Table 3, Fig. 1). The frequency of polymorphic sites in 210 accessions ranged from 0.005 to 0.143 with an average 0.030. To determine the precise position and nature of these polymorphic sites, several samples containing each of these sites were randomly sequenced. Accessions with the same polymorphic sites in EcoTILLING exhibited the same nucleotide changes in sequence with only one exception that one sample showing two polymorphic sites in EcoTILLING did not show nucleotide variation in sequencing. Thus, sequencing confirmed 11 SNPs in the targeted region of HSP17.8 at a frequency of 1 SNP per 54.5 bp. Among the 11 SNPs, 10 were from coding regions and one from non-coding region, and seven were transitions (C-T and A-G) and four were transversions (A-C, A-T, C-G and G-T). Among the ten coding SNPs, six were nonsynonymous mutations s and four were silent synonymous mutations. Five of the missense mutations were predicted to have a severe effect on the HSP17.8 protein function.

[Figure omitted. See PDF.]

Table 3. List of nucleotide polymorphisms in HSP17.8 and their effects on codon frequencies.

https://doi.org/10.1371/journal.pone.0056816.t003

The nucleotide diversity in the sequenced region of HSP17.8 as measured by π (pairwise nucleotide diversity) was 0.00188 among the 210 accessions. Across different geographic regions, the range of π values spanned from 0.00049 for African accessions (40 accessions) to 0.00212 for Middle East Asian accessions (44 accessions). Similarly, π for H. spontaneum was higher than for H. vulgare landrace accessions. In addition, 2-row and 6-row barley had similar nucleotide diversity (Table 4). To test whether the SNPs in the targeted region of HSP17.8 were neutral mutations, Tajima's D statistics (Tajima 1989) were estimated. Tajima's D values were negative for all sub-populations except for Africa population. However none of the values was statistically significant (P<0.05) (Table 4), thus, nucleotide variations in HSP17.8 gene could be the results from the standard neutral selection.

[Figure omitted. See PDF.]

Table 4. Frequency of HSP17.8 haplotypes, nucleotide diversity (π), haplotype diversity (HD) and Tajima's D test for different barley populations.

https://doi.org/10.1371/journal.pone.0056816.t004

Haplotype diversity analysis

For the 11 SNPs confirmed by sequencing, nine distinguishable haplotypes were detected across 210 accessions (Table 5). The level of haplotype diversity was 0.526, with the frequency of each haplotype shown in Table 4. Among all the haplotypes, three major haplotypes were detected in 210 accessions, with the cumulative frequency of the first three haplotypes being 0.9. Haplotype H1 was found in more than two-thirds of the accessions screened. Haplotype H2 was observed in 30 accessions (14.3%), and haplotype H3 observed in 21 accessions (10.0%). The frequency of the other haplotypes (H4 to H9) was low, between 1.0% and 2.4% (Table 4).

[Figure omitted. See PDF.]

Table 5. Distribution of polymorphic SNPs across nine HSP17.8 haplotypes.

https://doi.org/10.1371/journal.pone.0056816.t005

The frequencies of HSP17.8 haplotypes differed markedly across geographical regions (Table 4 and Fig. 2). Haplotype diversity ranged from 0 in Europe (6 accessions) to 0.722 in the Middle East Asia (44 accessions). This was particularly evident for the haplotype H2 which was absent in Africa, Middle East Asia and Europe, but most frequent in North East Asia (0.264). The less frequently identified haplotypes were confined to specific geographic regions. Of the six haplotypes present in <10% of accessions sampled, five were unique to one region (Asia); three of which were exclusive to accessions from Middle East Asia; and two from North East Asia. In addition, H. spontaneum and H. vulgare landraces were widely separated into six haplotypes with only three common between the two groups (Table 4). When comparing 2-row and 6-row barley, more haplotypes and higher diversity were observed in 2-row barley (Table 4).

[Figure omitted. See PDF.]

Figure 2. Composition of HSP17.8 haplotype in accessions of different geographic region.

AFR: Africa, APS: Arabian Peninsula, EUR: Europe, MEA: Middle East Asia, NEA: North East Asia. H (H1–H9) represents haplotype as described in Table 5.

https://doi.org/10.1371/journal.pone.0056816.g002

Association between SNPs and phenotypic traits

Because two SNPs (positions 196 bp and 483 bp in HSP17.8) were rare alleles (frequency <1%), and three SNPs (positions 104 bp, 267 bp and 300 bp) or two SNPs (positions 525 bp and 565 bp) were in complete linkage disequilibrium, five SNPs (positions 196 bp, 483 bp, 104 bp, 267 bp and 565 bp) were excluded from further analysis. Although two SNPs (position 469 and 582 bp) were synonymous variations, they may be in complete linkage disequilibrium with other non-synonymous mutations, which could result in biological changes in the organism. Therefore, these two SNPs were further analyzed. Thus only six SNPs (positions 300 bp, 428 bp, 469 bp, 525 bp, 582 bp and 599 bp) were used in association analysis. Significant association (P<0.05) was observed for four distinct SNPs with at least one of the evaluated traits except for spike length, and four of the associations were highly significant (P<0.01) (Table 6). These associated SNPs explained 2.0% to 4.6% of the variation for individual traits. FLA was significantly associated (P<0.01) with one SNP (position 469 bp in HSP17.8) in 2009/2010, which explained 3.9% of the phenotypic variation. NGS was significantly associated (P<0.01) with one SNP at position 599 bp, which explained 3.3% and 4.6% phenotypic variation in 2009/2010 and 2010/2011, respectively. One SNP (position 300 in HSP17.8) significantly associated (P<0.05) with PH was observed, accounting for approximately 2% phenotypic variation in 2009/2010 and 2010/2011. LC showed significant (P<0.05) association with three SNPs (position 469 bp, 525 bp and 599 bp in HSP17.8), which explained 2.8%, 3.9% and 2.4% of the phenotypic variation, respectively. TGW was significantly associated (P<0.05) with one SNP at position 300 bp in 2009/2010, which explained 2.6% phenotypic variation.

[Figure omitted. See PDF.]

Table 6. SNPs of HSP17.8 associated with agronomic traits of barley using a significance level corresponding to α = 0.05.

https://doi.org/10.1371/journal.pone.0056816.t006

Discussion

Natural variation in the barley population

To characterize genetic variation in HSP17.8, 171 H. vulgare landraces and 39 H. spontaneum accessions from 30 countries in three continents were characterized for allele diversity. EcoTILLING revealed 11 SNPs and nine haplotypes after analysis of 126,000 bp sequences in HSP17.8. Of the 11 unique SNPs, ten were in the coding region, which included four silent synonymous mutations and six missense mutations. Of the six missense mutations, four SNPs were predicted by the SIFT program to severely affect protein function, and one mutation was predicted by PSSM program to severely damage the function of the predicted protein. Because there is a 20% false-positive error in SIFT [31], some mutations predicted to be deleterious may be functionally neutral. However, these scores may be useful in prioritizing mutations for further study and analysis of possible contributions and roles of HSP17.8 in stress tolerance.

Tajima's D is used to measure deviation from neutral evolution by comparing diversity estimates based on nucleotide diversity (θ) and average pairwise nucleotide diversity (π). Parameters θ and π are affected differently by natural selection. In this study, Tajima's D neutrality test revealed no evidence of natural selection for HSP17.8 (Table 4) but under some kind of purifying selection as revealed by a high negative value of Tajima's D. This insignificant result may be attributed to the low number of SNPs observed, which weakens the neutrality test. This result agrees with previous reports on other functional genes in barley [36] and CPsHSP-2 in Machilus kusano [37]. However, the Tajima's D value in Africa was positive, which may result from balancing selection or bottleneck effect.

Compared to previous reports on barley [36], [38], nucleotide diversity (π = 0.00118) of HSP17.8 was lower, and the level of haplotype diversity (0.526) was higher in this study. The average frequency of SNPs was 1 per 54.5 bp, which was similar to that found in Chen et al. (1 SNP/53.8 bp) [39], lower than Zeng et al.(1 SNP/9.8 bp) [40], and higher than Rostoks et al. (1 SNP/200 bp) [41]. The discrepancy in SNP frequency among studies may be due to differences in genomic regions assayed, and number, content and geographic origins of germplasm used [38], [42], [43]. As predicted, H. spontaneum had a higher nucleotide diversity (π) and haplotype diversity than H. vulgare landraces in this study, which agrees with previous reports [38], [44], [45]. It is likely that H. vulgare landraces would have gone through a population bottleneck during domestication, which resulted in a reduction in genetic diversity[38], [45]. Comparing genetic diversity across geographical regions, accessions from Middle East Asia exhibited the highest nucleotide diversity (π), which agrees with studies of Malysheva-Otto et al. [46] and Varshney et al. [23]. One plausible reason was that the Middle East's ‘fertile crescent’ region was the center of origin of barley and the main distribution region of wild barley [47], [48]. This important result may indicate that variations in HSP17.8 in wild barleys and Middle East landraces might offer elite alleles for the improvement of stress tolerance in barley. To date, a number of wild barley and Middle East types have been used in barley breeding programs [49].

Different ecotypes of barley germplasm show their specific traits in terms of ecology, morphology and physiology [50]. Different ecotype populations from several different geographic regions may include abundant genetic variation associated with phenotypic traits. As extremely low nucleotide diversity could be observed in one ecotype, only a small number of ecotype populations may be sufficient to capture already adequate genetic variability [51]. Since wild and cultivated barley are cross-compatible, it is possible to increase the genetic diversity of barley using wild ancestors as a parent in crosses with cultivated barley [52]. There is evidence that the introgression of chromosome segments from wild progenitors improves agronomic performance of an elite cultivar [52]–[54].

Association analysis between SNPs and phenotypic traits

In our previous studies, HSP17.8 in five barley accessions, namely Martin, HS4-1, Moroc9-75, Tadmor and WI2291, was up-regulated under drought stress with the expression fold changes of 15.67, 19.43, 23.92, 2.48 and 1.15, respectively [12], [55]. Simple linear regression analysis between expression fold changes and six agronomic traits was performed in five barley accessions, and significant correlations were observed for PH, NGS, TKW, FLA and LC with the R2 values of 0.138, 0.163, 0.169, 0.134 and 0.218, respectively, suggesting the high expression level of HSP17.8 is likely associated with PH, NGS, TKW, FLA and LC. In addition, BLAST search in GenBank (http://www.ncbi.nlm.nih.gov/genbank/) showed that HSP17.8 exhibited nucleotide sequence homology to a barley EST (GenBank accession AL510041) that was located in 63.4 cM on chromosome 4H [56]. Several previous studies also reported chromosome 4H had QTLs for PH [57]–[59], TGW [57]–[59], NGS [58]–[60], LC [61] and FLA [62]. Therefore, HSP17.8 may have cis affects for the QTLs of PH, TGW, NGS, LC, and FLA located on chromosome 4H.

Association mapping has recently emerged as an alternative approach to mapping QTL and genes associated with quantitative traits using a diverse collection of germplasm lines or breeding materials [63]. Compared to traditional QTL mapping, association mapping is faster and provides greater capacity and power for QTL/gene detection [19], [20]. In addition, candidate-gene association analysis is more precise than genome-wide association analysis. Because association analysis links specific nucleotide polymorphisms to trait variations, one may be able to associate SNPs with specific biological effects [64]. Many agronomically important traits are controlled by QTLs [65]. By association analysis of 816 genome-wide markers, Varshney et al. [66] identified one to eight significant QTLs for nine agronomic traits in 223 barley accessions. With 204 polymorphisms in 24 transcription families, Yu et al. [67] found three genes associated with the drought tolerance index and five genes with the drought tolerance level in 95 diverse rice landraces. In the present study, four SNPs in the HSP17.8 gene were significantly associated with at least one agronomic trait investigated. Three SNPs at positions 469 bp, 525 bp and 599 bp of HSP17.8 were significantly associated with LC, one SNP at position 469 bp with FLA, one SNP at position 599 bp with NGS, and one SNP at position 300 bp with PH and TGW. Of the four SNPs associated with phenotypic traits, three at positions 300 bp, 525 bp and 599 bp in HSP17.8 were missense mutations which, according to SIFT, were predicted to severely affect protein function. Association analysis showed that the three missense mutations were deleterious mutations. Due to low minor allele frequency for two SNPs (positions 300 and 525), their association results should be interpreted with caution. However, one SNP at position 469 associated with phenotypic traits was a synonymous mutation that may not change the structure of the gene product. A plausible reason for this result is the hitchhiking effort of a locus in positive selection [68] or a false positive association. Thus, these SNPs identified by association analysis need to be validated for individual cultivars involved in crosses before they can be applied to marker-assisted selection in the progeny [69], [70]. The findings from this study indicate that further research of these newly detected SNPs in HSP17.8 is necessary to evaluate their possible influence on agronomic traits and usefulness as future selection markers in modern barley breeding program.

Conclusions

In this study, 11 SNPs were detected and nine unique haplotypes were identified in HSP17.8 among 210 barley accessions collected from 30 countries by EcoTILLING technology, and 5 SNPs with missense changes are predicted to be deleterious to protein function. Four SNPs significantly associated with agronomic traits identified in this study can be used as DNA markers in marker-assisted selection to improve these agronomic traits after further validation of their functions in individual cultivars.

Supporting Information

[Figure omitted. See PDF.]

Table S1.

Details on barley accessions used in this study.

https://doi.org/10.1371/journal.pone.0056816.s001

(DOC)

Table S2.

Primer sequences used for PCR amplification of HSP17.8 .

https://doi.org/10.1371/journal.pone.0056816.s002

(DOC)

Table S3.

SSR marker information used in evaluation of population structure.

https://doi.org/10.1371/journal.pone.0056816.s003

(DOC)

Author Contributions

Conceived and designed the experiments: PG RL YX. Performed the experiments: YX RL ZN PG. Analyzed the data: YX RL ZN GB. Contributed reagents/materials/analysis tools: PG RL ZN GB RKV MB. Wrote the paper: YX GB KHMS PG RKV GY MB.

Citation: Xia Y, Li R, Ning Z, Bai G, Siddique KHM, Yan G, et al. (2013) Single Nucleotide Polymorphisms in HSP17.8 and Their Association with Agronomic Traits in Barley. PLoS ONE 8(2): e56816. https://doi.org/10.1371/journal.pone.0056816

References

1. Sun W, Montagu MV, Verbruggen N (2002) Small heat shock proteins and stress tolerance in plants. Biochim Biophys Acta 1577: 1–9.

2. Hu W, Hu G, Han B (2009) Genome-wide survey and expression profiling of heat shock proteins and heat shock factors revealed overlapped and stress specific response under abiotic stresses in rice. Plant Sci 176: 583–590.

3. Kim DH, Xu ZY, Na YJ, Yoo YJ, Lee J, et al. (2011) Small heat shock protein HSP17.8 functions as an AKR2A cofactor in the targeting of chloroplast outer membrane proteins in Arabidopsis. Plant Physiol 157: 132–146.

4. Sato Y, Yokoya S (2008) Enhanced tolerance to drought stress in transgenic rice plants overexpressing a small heat-shock protein, sHSP17.7. Plant Cell Rep 27: 329–334.

5. Song NH, Ahn YJ (2011) DcHSP17.7, a small heat shock protein in carrot, is tissue-specifically expressed under salt stress and confers tolerance to salinity. New Biotechnol 28: 6.

6. Timperio AM, Egidi MG, Zolla L (2008) Proteomics applied on plant abiotic stresses: role of heat shock proteins (HSP). J Proteom 71: 391–411.

7. Waters ER, Lee GJ, Vierling E (1996) Evolution, structure and function of the small heat shock proteins in plants. J Exp Bot 47: 325–338.

8. Liu X, Huang W, Li M, Wu Q (2005) Purification and characterization of two small heat shock proteins from Anabaena sp. PCC 7120. IUBMB Life 57: 449–454.

9. Kant P, Gordon M, Kant S, Zolla G, Davydov O, et al. (2008) Functional-genomics-based identification of genes that regulate Arabidopsis responses to multiple abiotic stresses. Plant Cell Environ 31: 697–714.

10. Shou H, Bordallo P, Fan JB, Yeakley JM, Bibikova M, et al. (2004) Expression of an active tobacco mitogen-activated protein kinase kinase kinase enhances freezing tolerance in transgenic maize. PNAS 101: 3298–3303.

11. Grigorova B, Vaseva II, Demirevska K, Feller U (2011) Expression of selected heat shock proteins after individually applied and combined drought and heat stress. Acta Physiol Plant 33: 2041–2049.

12. Guo P, Baum M, Grando S, Ceccarelli S, Bai G, et al. (2009) Differentially expressed genes between drought-tolerant and drought-sensitive barley genotypes in response to drought stress during the reproductive stage. J Exp Bot 12: 3531–3544.

13. Witkin SS, Linhares IM (2010) Heat shock proteins, genital tract infections and reproductive outcome. Prokaryotic and Eukaryotic Heat Shock Proteins in Infectious Disease, Heat Shock Proteins 4: :241–256.

14. Comai L, Young K, Till BJ, Reynolds SH, Greene EA, et al. (2004) Efficient discovery of DNA polymorphisms in natural populations by EcoTILLING. Plant J 37: 778–786.

15. Gilchrist EJ, Haughn GW, Ying CC, Otto SP, Zhuang J, et al. (2006) Use of EcoTILLING as an efficient SNP discovery tool to survey genetic variation in wild populations of Populus trichocarpa. Mol Ecol 15: 1367–1378.

16. Nieto C, Piron F, Dalmais M, Marco CF, Moriones E, et al. (2007) EcoTILLING for the identification of allelic variants of melon eIF4E, a factor that controls virus susceptibility. BMC Plant Biology 7: 34.

17. Till BJ, Zerr T, Comai L, Henikoff S (2006) A protocol for TILLING and EcoTILLING in plants and animals. Nat Protoc 5: 2465–2477.

18. Wang N, Shi L, Tian F, Ning H, Wu X, et al. (2010) Assessment of FAE1 polymorphisms in three Brassica species using EcoTILLING and their association with differences in seed erucic acid contents. BMC Plant Biology 10: 137.

19. Yu J, Buckler ES (2006) Genetic association mapping and genome organization of maize. Curr Opin Biotech 17: 155–160.

20. Zhu C, Gore M, Buckler ES, Yu J (2008) Status and prospects of association mapping in plants. The Plant Genome 1: 5–20.

21. Wu D, Qiu L, Xu L, Ye L, Chen M, et al. (2011) Genetic variation of HvCBF genes and their association with salinity tolerance in Tibetan annual wild barley. PLoS one 6(7): e22938 Doi: 10.1371/journal.Pone.0022938.

22. Xia Y, Ning Z, Bai G, Li R, Yan G, et al. (2012) Allelic Variations of a Light Harvesting Chlorophyll A/B-Binding Protein Gene (Lhcb1) Associated with Agronomic Traits in Barley. PLoS ONE 7(5): e37573 .

23. Varshney RK, Baum M, Guo P, Grando S, Ceccarelli S, et al. (2010) Features of SNP and SSR diversity in a set of ICARDA barley germplasm collection. Mol Breeding 26: 229–242.

24. Gupta R, Guhey A, Jadhav A, Ahad I (2011) Physiological boost to improve the yield of rice germplasm under different water regimes. Research Journal of Agricultural Science 2: 87–90.

25. Wienholds E, van Eeden F, Kosters M, Mudde J, Plasterk RH, et al. (2003) Efficient target-selected mutagenesis in zebrafish. Genome Res 12: 2700–2707.

26. Guo P, Bai G, Li R, Carver B, Baum M (2007) Molecular characterization of atlas 66-derived wheat near-isogenic lines contrasting in aluminum (Al) tolerance. Agr Sci China 5: 522–528.

27. Wang GX, Tan MK, Rakshit SR, Saitoh H, Terauchi R, et al. (2007) Discovery of single-nucleotide mutations in acetolactate synthase genes by Ecotilling. Pestic Biochem Phys 88: 143–148.

28. Raghavan C, Naredo MEB, Wang HH, Atienza G, Liu B, et al. (2007) Rapid method for detecting SNPs on agarose gels and its application in candidate gene mapping. Mol Breeding 19: 87–101.

29. Guo PG, Li RH (2010) A method for isolation of CEL I enzyme from celery. Chinese Patent, ZL200910038725.9, 2010-9-29.

30. Zerr T, Henikoff S (2005) Automated band mapping in electrophoretic gel image using background information. Nucleic Acids Res 33: 2806–2812.

31. Ng PC, Henikoff S (2003) SIFT: predicting amino acid changes that affect protein function. Nucleic Acids Res 31: 3812–3814.

32. Taylor NE, Greene EA (2003) PARSESNP: a tool for the analysis of nucleotide polymorphisms. Nucleic Acids Res 31: 3808–3811.

33. Tajima F (1989) Statistical method for testing the neutral mutation hypothesis by DNA polymorphism. Genetics 123: 585–595.

34. Librado P, Rozas J (2009) DnaSP v5: a software for comprehensive analysis of DNA polymorphism data. Bioinformatics 25: 1451–1452.

35. Pritchard JK, Stephens M, Donnelly P (2000) Inference of population structure using multilocus genotype data. Genetics 155: 945–959.

36. Cseri A, Cserháti M, von Korff M, Nagy B, Horváth GV, et al. (2011) Allele mining and haplotype discovery in barley candidate genes for drought tolerance. Euphytica 181: 341–356.

37. Chung JD, Lin TP, Cheng YP, Hwang SY (2010) Geographical variation and differential selection modes of paralogs of chloroplast small heat shock protein genes in Machilus kusanoi (Lauraceae). Bot Stud 51: 27–36.

38. Russell J, Booth A, Fuller J, Harrower B, Hedley P, et al. (2004) A comparison of sequence-based polymorphism and haplotype content in transcribed and anonymous regions of the barley genome. Genome 47: 389–398.

39. Chen G, Yang S, Li W, Li W, Wei Y, et al. (2009) Variations of single nucleotide polymorphism (SNP) of Hordoindoline a gene (hina) and its relationships with kernel hardness in 92 North American barley varieties (Hordeum vulgare L.). Agr Sci China 12: 1440–1447.

40. Zeng X, Wei Y, Jiang Q, Qi P, Zheng Y (2009) SNP analysis and haplotype identification in chymotrypsin inhibitor-2 (CI-2) gene of barley. Agr Sci China 8: 8–14.

41. Rostoks N, Borevitz JO, Hedley PE, Russell J, Mudies S, et al. (2005) Single-feature polymorphism discovery in the barley transcriptome. Genome Biol 6: R54.

42. Bundock PC, Henry RJ (2004) Single nucleotide polymorphism, haplotype diversity and recombination in the Isa gene of barley. Theor Appl Genet 109: 543–551.

43. Duran C, Appleby N, Vardy M, Imelfort M, Edwards D, et al. (2009) Single nucleotide polymorphism discover in barley using autoSNPdb. Plant Biotechnol J 7: 326–333.

44. Kilian B, Ozkan H, Kohl J, von Haeseler A, Barale F, et al. (2006) Haplotype structure at seven barley genes: relevance to gene pool bottlenecks, phylogeny of ear type and site of barley domestication. Mol Genet Genomics 276: 230–241.

45. Tanksley SD, McCouch SR (1997) Seed banks and molecular maps: unlocking genetic potential from the wild. Science 277: 1063–1066.

46. Malysheva-Otto LV, Ganal MW, Röder MS (2006) Analysis of molecular diversity, population structure and linkage disequilibrium in a worldwide survey of cultivated barley germplasm (Hordeum vulgare L.). BMC Genetics 7: 6.

47. Nevo E (1992) Origin, evolution, population genetics and resources for breeding of wild barley, Hordeum spontaneum in the Fertile Crescent. In: Shewry PR(ed) Barley: genetics, biochemistry, molecular biology and biotechnology. C.A.B. International, Wallingford, Oxford, pp19–43.

48. Ceccarelli S, Grando S, van Leur JAG (1995) Barley landraces of the fertile crescent offer new breeding options for stress environments. Diversity 11: 112–113.

49. Ellis RP, Forster BP, Robinson D, Handley LL, Gordon DC, et al. (2000) Wild barley: a source of genes for crop improvement in the 21st century? J Exp Bot 51: 9–17.

50. Abdellaoui R, M'Hamed HC, Naceur MB, Bettaïeb-Kaab L, Hanida JB (2007) Morpho-physiological and molecular characterization of some Tunisian barley ecotypes. Asian Journal of Plant Sciences 6: 261–268.

51. Peter-Schmid MKI, Boller B, Kölliker R (2008) Habitat and management affect genetic structure of Festuca pratensis but not Lolium multiflorum ecotype populations. Plant Breeding 127: 510–517.

52. Inostroza L, del Pozo A, Matus I, Castillo D, Hayes P, et al. (2009) Association mapping of plant height, yield, and yield stability in recombinant chromosome substitution lines (RCSLs) using Hordeum vulgare subsp. spontaneum as a source of donor alleles in a Hordeum vulgare subsp. vulgare background. Mol Breeding 23: 365–376.

53. Lakew B, Eglinton J, Henry RJ, Baum M, Grando S, et al. (2011) The potential contribution of wild barley (Hordeum vulgare ssp. spontaneum) germplasm to drought tolerance of cultivated barley (H. vulgare ssp. vulgare). Field Crops Research 120: 161–168.

54. Schmalenbach I, Léon J, Pillen K (2008) Identification and verification of QTLs for agronomic traits using wild barley introgression lines. Theor Appl Genet 18: 483–497.

55. Guo PG, Baum M, Li RH, Grando S, Varshney RK, et al. (2007) Transcriptional analysis of barley genes in response to drought stress at the reproductive growth stage using affymetrix Barley 1 genechip. J Guangzhou Univ (Nat Sci Ed) 5: 32–36.

56. Stein N, Prasad M, Scholz U, Thiel T, Zhang H, et al. (2007) A 1000 loci transcript map of the barley genome: new anchoring points for integrative grass genomics. Theor Appl Genet 114: 823–839.

57. Li JZ, Huang XQ, Heinrichs F, Ganal MW, Röder MS (2006) Analysis of QTLs for yield components, agronomic traits, and disease resistance in an advanced backcross population of spring barley. Genome 49: 454–466.

58. Teulat B, Merah O, Souyris I, This D (2001) QTLs for agronomic traits from a Mediterranean barley progeny grown in several environments. Theor Appl Genet 103: 774–787.

59. Schmalenbach I, Léon J, Pillen K (2009) Identification and verification of QTLs for agronomic traits using wild barley introgression lines. Theor Appl Genet 118: 483–497.

60. Sun D, Ren W, Sun G, Peng J (2011) Molecular diversity and association mapping of quantitative traits in Tibetan wild and worldwide originated barley (Hordeum vulgare L.) germplasm. Euphytica 178: 31–43.

61. Guo P, Baum M, Varshney RK, Graner A, Grando S, et al. (2008) QTLs for chlorophyll and chlorophyll fluorescence parameters in barley under post-flowering drought. Euphytica 163: 203–214.

62. Mohammed KAH (2004) Improving crop varieties of spring barley for drought and heat tolerance with AB-QTL analysis. PhD Thesis, Bonn University, Germany.

63. Myles S, Peiffer J, Brown PJ, Ersoz ES, Zhang Z, et al. (2009) Association mapping: critical considerations shift from genotyping to experimental design. Plant Cell 21: 2194–2202.

64. Szalma SJ, Buckler ES, Snook ME, McMullen MD (2005) Association analysis of candidate genes for maysin and chlorogenic acid accumulation in maize silks. Theor Appl Genet 110: 1324–1333.

65. Paterson AH, Lander ES, Hewitt JD, Peterson S, Lincoln SE, et al. (1988) Resolution of quantitative traits into Mendelian factors by using a complete linkage map of restriction fragment length polymorphisms. Nature 335: 721–726.

66. Varshney RK, Paulo MJ, Grando S, van Eeuwijk FA, Keizer LCP, et al. (2012) Genome wide association analyses for drought tolerance related traits in barley (Hordeum vulgare L.). Field Crops Research 126: 171–180.

67. Yu S, Liao F, Wang F, Wen W, Li J, et al. (2012) Identification of rice transcription factors associated with drought tolerance using the ecotilling method. PLoS ONE 7(2): e30765.

68. Thornton KR, Jensen JD, Becquet C, Andolfatto P (2007) Progress and prospects in mapping recent selection in the genome. Heredity 98: 340–348.

69. Breseghello F, Sorrells ME (2006) Association mapping of kernel size and milling quality in wheat (Triticum aestivum L.) cultivars. Genetics 172: 1165–1177.

70. Peng JH, Bai Y, Haley SD, Lapitan NLV (2009) Microsatellite-based molecular diversity of bread wheat germplasm and association mapping of wheat resistance to the Russian wheat aphid. Genetica 135: 95–122.