Skip to main content

A microsatellite diversity analysis and the development of core-set germplasm in a large hulless barley (Hordeum vulgare L.) collection

Abstract

Background

Clarifying genetic diversity in a large germplasm resource plays important roles in experimental designs that provides flexible utility in fundamental research and breeding in crops. However, the work is limited due to small collections of barley that are insufficient representatives.

Results

In the present study, we collected 562 hulless barley (Hordeum vulgare L.) accessions with worldwide geographic origins and evaluated their genetic variability and relatedness based on 93 simple sequence repeat (SSR) markers. In an integrated analysis of the population structure, analysis of molecular variance (AMOVA) and pairwise F ST, the 562 barley accessions exhibited a strong stratification that allowed for them to be divided into two major subpopulations (p1 and p2) and an admixture subpopulation, with 93, 408 and 61 accessions, respectively. In a neutral test, considerable proportions of SSR alleles expressed the strong non-neutrality in specific subpopulations (44 and 37), which are probably responsible for population differentiation. To reduce the diversity redundancy in large barley collections, we delicately selected a core set of 200 barley accessions as a tradeoff between diversity and representativeness in an easily handled population. In comparing the 562 barley accessions, the core barley set accounted for 96.2% of allelic diversity and 93% to 95% of phenotypic variability, whereas it exhibited a significant enhancement in minor allelic frequencies, which probably benefit association mapping in the barley core set.

Conclusions

The results provided additional insight into the genetic structure in a large barley germplasm resource, from which an easily manageable barley core set was identified, demonstrating the great potential for discovering key QTLs and ultimately facilitating barley breeding progress.

Background

Hulless barley (Hordeum vulgare L.) with naked caryopsis is mainly distributed in Tibet and nearby areas. This was believed to be one of the important centers of origin of cultivated barley [1]. The hulless barley has become one of the most important staple crops in the Tibet region for the last 3500―4000 years because it is capable of adapting to the extreme climates and high altitudes. To date, hulless barley has predominantly occupied over 70% of the crop lands in Tibet [2], while it is also rarely distributed in Europe, the USA and Australia, as an important resource to feed animals. In recent years, hulless barley has proved its nutritional and economic value due to its high β-glucan content, which inhibited cholesterol synthesis, and the absence of grain husks probably reduced the processing cost of the barley food industry.

In breeding programs, a priori knowledge of genetic diversity and pairwise relatedness provides beneficial clues for efficiently utilizing a large collection of genetic resources. In the Himalayan region, particularly in Tibet, hulless barley showed a considerable contribution to the diversity of the barley germplasm worldwide. To date, there have been many studies reporting the genetic and phenotypic variability harbored in Himalayan barley accessions. In a survey of Nepalese naked barley germplasm, vastly morphologic [3] and disease-resistant [4] variations were found between and within the landrace populations. Recently, Pandey et al. [5] revealed a complex population structure of 107 hulless barley landraces, possibly due to the geographic differentiation among the barley populations. However, the observations in these studies are limited due to the small collections of barley and their insufficient representativeness.

The draft genomic sequence of the Tibetan hulless landrace “Lasa Goumang” offers opportunities to facilitate genetic improvement by the identification of gene or quantitative trait loci (QTL) of important traits [2]. It is therefore urgently important to build relevant populations for genetic mapping, i.e., linkage mapping and association mapping [6]. Comparatively, natural population-based genome-wide association studies (GWAS) are more powerful and cost-efficient than linkage mapping based on bi-parental crossing populations because the former usually accumulated more recombinant events and showed variation in their genomics and phenomics [7]. Importantly, the natural population confers a straightforward link between identification and breeding applications of key genes responsible for agronomic traits [8]. The trade-off between the experimental cost and mapping power is an advantage of the core germplasm set that harbors the majority (e.g., >90%) of the genetic variability of the total resource by a significantly reduced resource size, which would be an optimal choice in crop genetics. A number of GWAS on a manageable number of accessions (300―500) have been reported to dissect the genetic basis and even isolate single genes underlying key traits in rice [9, 10], maize [11], oilseed rape [12] and peanut [13]. In regards to breeding, we aim to identify the core germplasm that is capable of reducing the laboratory burden in encountering genetically redundant materials, which would increase breeding gain with limited cost.

In this study, we assembled a large-scale germplasm collection of 562 hulless barley accessions worldwide, but mainly in the Tibet region, which were genotyped using 93 simple sequence repeats (SSRs) dispersed across the whole genome. The objectives of this research are to (1) evaluate the genetic variability of the large hulless barley germplasm resource, (2) investigate the population structure and differentiation, and (3) identify and evaluate the utility of an optimal core set in hulless barley.

Methods

A collection of worldwide hulless barley accessions

In this study, a set of 562 naked barley accessions were collected worldwide (Fig. 1a) the majority of which were derived from the Himalayan region, typically in China. This included 426 accessions from Tibet, 29 from Qinghai, and 27 from Sichuan. The remaining accessions were diversely dispersed across different continents (Fig. 1a). All barley accessions involved 98 cultivated varieties, 414 landraces and 50 wild resources. The materials are provided by Chinese Academy of Agricultural Science (CAAS) and accessible in the National Infrastructure of Plant Germplasm Resource in China (http://www.cgris.net/pt/). Detailed information about the 562 hulless barley accessions is listed in Additional file 1.

Fig. 1
figure 1

Geographic coordinates, phylogenic and principal component analyses of the 562 hulless barley accessions. a Colored dots refer to the sampling regions of the barley accessions, while the color gradient corresponds to the number of accessions sampled in each region. b All 562 barley accessions were clustered into two apparent clades, colored red and blue, based on the structural analysis. c The first two principal components (PC) explained over 10% of the molecular variance and discriminated the whole population into three groups, which is largely consistent with the subpopulations inferred by structural analysis

SSR genotyping

For each accession, the genomic DNA was extracted from the young leaf tissues derived from a single plant following the modified cetyltrimethylammonium bromide (CTAB) protocol [14]. The quality of the DNA was evaluated on a 1% agarose gel through comparison with the uncut lambda DNA.

A set of 100 SSR primers were used in this study. Of these, 23 primers were retrieved from the public Graingenes 2.0 database (http://wheat.pw.usda.gov/GG2/) and 77 primers were newly developed by our laboratory based on the hulless barley reference genome [2]. First, 10 randomly selected barley accessions were used to test the polymorphism and amplification efficiency in all 100 SSR primers. Subsequently, all 562 hulless barley accessions were genotyped using 93 high-quality SSR primers, including 22 SSR primers from the public wheat database and 71 newly developed SSRs with the prefix “QK”. The information on the 71 new SSR primers is listed in Additional file 2. The polymerase chain reaction (PCR) analysis was performed according to the previous protocol [15], except that the SSR primers were labelled with fluorescent dyes. The SSR polymorphism was assayed by capillary electrophoresis (ABI 3730 Genetic Analyzer Applied Biosystems). The SSR primer alleles were quantified based on electrophoretic data using the GeneMarker V2.1 software. To avoid genotyping artifacts influenced by the SSR paralogs [16], we recoded the SSR alleles in binary data (1 or 0) to indicate the presence or absence of a certain allele in the whole population. To facilitate data analysis, the binary data per allele was treated as a dominant locus. The statistical value of the SSR primers was defined as an averaged value across all the alleles in a specific SSR.

Genetic diversity

For each SSR primer, we calculated the allelic count, gene diversity (or expected heterozygosity) and polymorphism content (PIC) using PowerMarker version 3.51 [17]. The Shannon index per SSR was calculated using PopGene version 1.31 software [18]. Given the population size affect, the allelic count in a certain population and the allelic richness per SSR were calculated by dividing the total allelic count per SSR by the population size. The differences in the gene diversity, PIC, and allelic richness between the inferred subpopulations was evaluated based on the Wilcoxon’s paired test implemented in SAS 8.02 [19].

Population structure and differentiation analyses

The STRUCTURE v2.2 software [20] based on the Bayesian model was used to explore the population structure of 562 barley accessions based on 93 SSRs. We simulated the number of subpopulations (k) from one to ten, with five replications. For each replication, the posterior probability, LnP(D), of the model was estimated based on 10,000 MCMC (Markov Chain Monte Carlo) iterations before using a burn-in length of 10,000 times with the admixture and related frequency model. The optimal k value was determined by the LnP(D) value with an ad hoc statistic delta k based on the rate of change in the LnP(D) between successive k [21]. Each hulless barley accession was assigned to corresponding subpopulations based on its maximum membership probability. Otherwise, the accession was classified into a “mixed” subpopulation if all the membership probabilities were below 0.65 [22]. The inferred subpopulations based on the STRUCTURE analysis were deployed in the following analysis.

To better illustrate the pairwise relationship across 562 hulless barley accessions, we performed principal component analysis (PCA), which was implemented in SAS 8.02 [19], and a neighbor-joining (N-J) tree based on the Nei’s distance using MEGA 4.0 [23]. To evaluate the dependency between the population structure and several clustering factors, such as geographic origin, eco-type, row-type and evolutionary type, the contingency table test based on an R function “chisq.test” was used (P < 0.05). Based on the subpopulations inferred by structural analysis, we carried out analysis of molecular variance (AMOVA) to assess the population differentiation using Arlequin V3.1 software [24], with 1000 permutations and a molecular distance based on the sum of squared size differences. On the other hand, the pairwise differentiation between subpopulations was evaluated by pairwise F ST, a measurement of expected heterozygosity within subpopulations compared to the total population [25].

To detect the link of population structure with allelic selection, we performed the Ewens-Watterson’s neutrality test [26, 27] for each binary locus (i.e., the SSR allele) in each subpopulation using PopGene version 1.31 software [18]. This was used to test whether the observed allele frequency is significantly higher or lower than the expected frequency by simulation under neutrality expectations [26, 27].

The assembly of hulless barley core set

We used a multi-purpose core selection method to assemble the diverse core subset from the large germplasm collections with genome-covered SSR genotypes, implemented in the R package ‘corehunter’ [28]. To perform a GWAS in the core set, the high power of identifying strong signals necessitates sufficient allele repeats that are linked with interested phenotypes in a given population [6]. In this study, we expected a barley core set of 200 accessions, as the equivalent size has shown sufficient power in identifying major QTLs in self-mating species [10]. We deployed two parameters and weighted them equally to evaluate the hulless barley collection iteratively: (1) average entry-to-nearest-entry distance (EN), the maximization of which increased the average distance between each selected barley and the closest other selected item in the core; (2) gene diversity, which typically reflects the frequency to a balanced extent because the allele count in a specific locus was constantly two in this study. The optimal core set would potentially exhibit a high allelic diversity and low redundancy between any individual pairs. To independently evaluate the reliability of the assembled hulless barley core set, we deployed four agriculturally important traits (plant height, ear length, grain length and thousand seed weight) to assess the ability to maintain the phenotypic variability compared to the original large collection. Finally, we selected several informative SSRs to establish the fingerprints of the hulless barley core set with 200 accessions.

Results

Population structure and genetic diversity in the 562 hulless barley accessions

The population structure of a worldwide collection of 562 hulless barley accessions was investigated using 93 SSR markers. The maximal posterior probability (LnP(D)) of the Bayesian model was estimated at the possible subpopulations (k) from 1 to 10. We found that the most apparent change of the LnP(D) appeared when k increased from 1 to 2 (Additional file 3a). In addition, a sharp peak in delta k appeared at k = 2 (Additional file 3b). Accordingly, the 562 barley accessions were classified into two major subpopulations, designated p1 and p2, as well as a mixed subpopulation (Additional file 3c). The p1 subpopulation contained 93 accessions, of which, 85 accessions were from China (78 from Tibet), 2 were from Australia, 4 were from Canada and 2 were from Mexico. In p1, there were 84 barley landraces, 7 barley cultivars and 2 wild barley resources. The p2 subpopulation contained 408 barley accessions, of which, 344 accessions were from China (296 from Tibet), 14 were from Mexico, 12 were from Canada, 10 were from Japan, 7 were from Australia and 1―5 were from other East Asia countries. In p2, there were 289 barley landraces, including 76 cultivars and 43 wild accessions. The mixed subpopulation contained 61 barley accessions, almost covering all the geographic regions sampled in the present study. The information on the geographic origins, types and inferred subpopulations is listed in detail in Additional file 1. We further employed N-J phylogenetic analysis and principal component analysis (PCA) to determine the genetic relationship across the 562 barley accessions. All 562 barley accessions were clustered into two apparent clades corresponding to the subpopulation determined by structural analysis (Fig. 1b). This pattern was additionally verified by the PCA plot, in which the top two principal components accounted for over 10% of the molecular variance (Fig. 1c). The contingency table test revealed that the inferred population structure was significantly correlated to the evolutionary barley types (P = 0.028). Specifically, the p2 subpopulation was enriched in wild barley relative to the p1 subpopulation (Additional file 1). Therefore, the two major structural populations were used in following analyses.

The genetic diversity in the total barley collection and inferred subpopulations were measured by the allele number, allelic richness, gene diversity, polymorphism information content (PIC) and Shannon index. A summary of the diversity analysis per SSR is listed in Additional file 4. In the total collection, the 93 SSR markers detected 516 alleles, with 5.55 alleles per SSR. The gene diversity, PIC and Shannon index were 0.189 (0.009–0.499), 0.157 (0.009–0.374) and 0.305 (0.029–0.692), respectively (Table 1 and Additional file 5). The p1 and p2 subpopulations contained 456 and 467 alleles, with 4.9 and 5.0 alleles per SSR, respectively. The number of alleles in p2 was slightly higher than that in p1, although the sample size of p2 was over fourfold that of p1. On the other hand, p1 had a similar number of alleles per SSR to p2 but significantly more allelic richness than p2 (P < 0.01). The level of gene diversity, PIC and Shannon index consistently were consistently higher in p1 than in p2, but the extent varied slightly (Table 1 and Additional file 5). The mixed subpopulation had almost one-seventh the sample size in p2 but harbored a comparable number of alleles, which probably led to the highest level of allelic richness in the mixed subpopulation (Table 1).

Table 1 Diversity in barley structurally inferred subpopulations

Genetic differentiation and neutrality between inferred subpopulations

We performed AMOVA and pairwise F ST to investigate the population differentiation between the inferred subpopulations. AMOVA results illustrated that 16.07% (P < 0.001) of the total molecular variation in the 562 hulless barley accessions was attributed to genetic differentiation between inferred subpopulations (Table 2). Pairwise F ST between the p1 and p2 subpopulations was 0.223 (P < 0.001), but the values between the mixed and p1 and the mixed and p2 were 0.079 and 0.075 (P < 0.001), respectively (Additional file 6). This suggested that p1 is significantly divergent from p2, but they were relatively genetically similar to the mixed subpopulation. A similar pattern of genetic differentiation among the inferred subpopulations was confirmed using the results of the pairwise Nei’s minimum distance (Additional file 6). The correlation coefficient between the F ST and Nei’s distance was 0.99 (P < 0.01).

Table 2 Analysis of molecular variance (AMOVA) among the inferred barley subpopulations based on STRUCTURE analysis

To further evaluate the influence of population structure, we deployed the Ewens-Watterson’s neutrality test to determine whether there are specific SSR alleles favored in the subpopulations. In summary, we found that there were 173 SSR alleles that showed significantly higher or lower frequencies than the expected value based on the neutral model (P < 0.01) (Fig. 2). The non-neutrality of these SSR alleles may be attributed to directional selections linked to population differentiation. Of all the non-neutral SSR alleles, there were 77 (44.5%), 69 (39.9%) and 86 (49.7%) alleles that were significantly discovered within p1, p2 and mixed subpopulations, respectively. We found that 9 non-neutral alleles were shared between p1 and p2, while 27 and 26 non-neutral alleles were shared between p1 and mixed and between p2 and mixed, respectively. Moreover, a considerable proportion of the non-neutral alleles that were specific to the subpopulations were identified, with nearly 25.4% and 21.4% of the non-neutral alleles being specific to p1 and p2, respectively (Fig. 2).

Fig. 2
figure 2

Genetic differentiation between the subpopulations was revealed by the neutrality analysis. The number labeled within the different sections of the Venn-plot represent the number of microsatellite alleles that showed significantly non-neutral frequency in specific subpopulations or shared between subpopulations

The assembly and characteristics of barley core resource

Considering the need for GWAS in barley, in this study, we assembled a core set of 200 hulless barley accessions (~1/3 of original population size) to sustain sufficient mapping power but largely reduce genotyping and phenotyping cost. The identifiers of the 200 hulless barley core accessions are indicated in Additional file 1. We found that the core accessions were evenly distributed across the raw collection of 562 accessions that were consistently revealed by the NJ phylogenic and PCA analyses (Additional file 7). Specifically, in the hulless barley core set, 34 accessions were selected from the p1 subpopulation, 148 from the p2 subpopulation and 18 from the mixed subpopulation, accounting for 29.5%–36.6% of the accessions per subpopulation.

To evaluate the reliability of the hulless barley core set, we compared the diversity between the core set and random sets assembled by chance with 1000 repeats. We found that the core set significantly captured more alleles (499 vs. 482; P = 7.8E-4; Fig. 3a) with a higher gene diversity (0.215 vs. 0.195; P = 9.3E-18; Fig. 3b), indicating that the barley core set is a reasonable outcome that is significantly better than random events. Interestingly, the assembled core selection had moved the average pairwise distance between the accessions from 0.17 to 0.26 (Fig. 3c), which potentially eliminated duplications and highly similar accessions, resulting in a highly representative assembled barley core set. In summary, the barley core set contained a total of 499 alleles, accounting for 96.7% of the existing alleles in the original 562 accessions. Typically, the allelic count per accession was comparable between the barley core and raw set in either subpopulation or the whole population (P > 0.05; Fig. 3d). For the gene diversity, the overall diversity value across all the SSRs was 0.221 in the core set, which is significantly beyond the value of 0.189 that was observed in the raw set (P = 0.028). The gene diversity was not significantly different between the barley core and raw set in any subpopulation, but p1 had a marginally higher diversity in the barley core set than in the raw set (Fig. 3d). The allelic frequency exhibited a similar trend with gene diversity, probably explaining that the overall diversity of the barley core set increased compared to the raw set. The summary of the diversity per SSR in the hulless barley core set is listed in Additional file 4.

Fig. 3
figure 3

Diversity analysis of the assembled barley core set. The barley core set of 200 accessions was mathematically selected from all 562 accessions using the R package ‘corehunter’. a The selected barley core set captured 499 microsatellite alleles accounting for 96.7% of the alleles harbored in the 562 accessions, which is significantly more than the 200 random accessions. b Gene diversity of the core set averaged across the SSR primers. It is significantly higher than the raw 562 accessions and 200 random accessions, which is probably due to the high allele maintenance and balanced allelic frequency in the core set. c Comparison of the Nei1972 distance for accession pairs between the barley core set and raw set; d The alleles captured by accession per se (left panel), gene diversity per SSR (middle panel) and minor allele frequency (MAF) of the alleles (right panel) for each subpopulation and whole population in the core (200 accessions) and raw (562 accessions) sets

To evaluate the maintenance of the phenotypic variability in the core set, we compared the hulless barley core and raw set in the four agronomic traits. We found that the ear length, grain length and the thousand seed weight showed a slightly higher coefficient of variation (CV) but plant height had a slightly lower CV in this barley core set than the raw set (Table 3). In addition, the Kolmogorov-Smirnov test revealed that all four traits followed a non-significantly different distribution (i.e., statistically approximately identical) between the core and raw sets (Additional file 8). It demonstrated that the barley core selection based on genetic data largely maintained phenotypic variability.

Table 3 Comparison of the coefficient of variation for the phenotypic value between the core barley set and the raw 562 barley accessions in this study

Regarding the aims in practice, we used ten informative SSR primers that had top10 diversities (see Additional file 4 in details) to construct the fingerprint of each accession in the barley core set (Additional file 9). The fingerprint provided an efficient way to discriminate identities of the barley core set in the breeding program and seed market.

Discussion

Genetic diversity and population structure in the large hulless barley resource

Understanding of the genetic diversity and population structure provides an opportunity to evaluate the applicable potential of a new germplasm resource in the breeding aims and fundamental studies. In this study, we employed 93 SSR primers to evaluate the genetic diversity in a large hulless barley collection of 562 accessions. On average, the SSR amplified 5.55 alleles per SSR, with a range from 2 to 16 alleles per SSR, which was comparable to the value (5.54 in average) of the 44 SSRs in the 107 hulless barley landraces collected from the highlands of Nepal [5]. However, the PIC value was comparatively lower in this study, which may be partly explained by three reasons: (1) the difference in the SSR primers used in two studies; (2) the difference in the PIC calculation algorithm, including the decomposition of the multi-allelic SSRs into the bi-allelic loci, which would unavoidably result in a loss of information that may be proportional to the allele count; (3) the difference in the germplasm resource, including more barley landraces and wild type that significantly enriched the gene library, while the exotic germplasm potentially introduced private or rare alleles (more discussed hereafter).

The Bayesian-model-based structure analysis is widely used a method for the inference of hidden population structure in human and plant species [29, 30]. In this study, two major subpopulations were identified in the collection of 562 hulless barley accessions, the structural pattern of which was verified by the NJ phylogenic and PCA analyses (Fig. 1). Interestingly, we found that the two subpopulations were significantly correlated with the evolutionary types of the barley accessions (P = 0.028). The p2 population contained 10.5% of the wild barley accessions, nearly 5 fold more than p1. This difference may be reflected by the fact that the Himalayan vicinity is one of the centers for barley origin in which the majority of accession regions in p2 were collected. However, the two major subpopulations were not significantly correlated to the geographical origins, row-types and eco-types (Additional file 1). Given that the delta k additionally peaked at a k of 5 and 7, respectively (Additional file 3b), we therefore re-clustered the whole population into 5 and 7 subpopulations. It was probably a refinement of two major subpopulations that was identified in k = 2 (Additional file 10), which enabled these subpopulations to be significantly correlated to the geographic origins (P = 9.38E-7 and 1.78E-11) and row-types (P = 0.032 and 0.017) at k = 5 and 7, respectively (Additional file 11). In the following analysis, we focus on the major two subpopulations that explained a large variance in the population differentiation. We found that the alleles per SSR for the total population and two inferred subpopulations were comparable, but p2 had almost one-fifth the allelic richness, while it had over fourfold the sample size compared to p1 (Table 1). This suggested that, in the p2 subpopulation, the large amounts of exotic barley landraces and wild types may prevalently capture private alleles (the allele that existed specifically only in one accession) or rare alleles. It is reliable because all the diversity measurements (gene diversity, PIC and Shannon index) were consistently significantly lower in p2 compared to p1 (Table 1 and Additional file 4) and the median minor allele frequency (MAF) values were significantly lower in p2 than p1 (Fig. 3d).

Perspectives of the assembled hulless barley core set

In plant breeding, the core collection plays important roles in parent selection for cross breeding and seed conservation. The core collection aims to maximize the genetic diversity of a large germplasm resource in a manageable sized subset, the operation of which will largely reduce redundant labors in the limited breeding cycles. In this study, we exhaustively tested all the possibilities and selected the optimal one (core set), which simultaneously reached a maximum allele amount and genetic diversity. The barley core set captured a total of 96% of the SSR alleles that existed in the raw population but significantly elevated the gene diversity level. Additionally, the barley core set almost mimicked the phenotypic distribution of the four agronomic traits and retained the phenotypic variability of the original 562 hulless barley accessions. These findings demonstrate that the hulless barley core set represented sufficient variation in the genetic and phenomic layers.

GWAS has been increasing in popularity for unraveling the genetic basis of simple and complex traits in plants [31, 32]. Presently, the reference genome of the Tibetan hulless landrace “Lasa Goumang” has been released [2]. It opens an era of GWAS-based gene identification and genomic selection (GS) based on genetic improvement in barley breeding, as tremendous molecular markers would be easily retrieved by resequencing and bioinformatic analysis [33, 34]. The establishment of the multi-ends research population is significantly important in the barley research community. The use of a carefully selected core set of diverse accessions is a routine approach to facilitate GWAS [33], which is probably due to the trade-off between the mapping power and labor- and cost-requirement. However, many researchers argued that the naturally-existing population always harbors numerous rare variations, which probably missed the heritability as given by the low frequency variation, which is difficult to detect statistically by GWAS [35, 36]. Interestingly, to this point, the barley core set provided an optimal mapping population for GWAS, because it significantly balanced the allelic frequencies, elevating the median MAF from 0.077 to 0.11 (all) and from 0.045 to 0.081 (p2) (Fig. 3d), which potentially boosted the mapping power imbedded in this barley core set. In summary, the multi-end use of the core set provided a feasible way to bridge traditional and molecular breeding in hulless barley.

Conclusions

In this manuscript, we collected 562 hulless barley accessions with worldwide geographic origins (the largest collection to date) and evaluated their genetic variability and relatedness based on 93 SSR markers. The joint analyses of PCA, AMOVA, and pairwise-FST and the neutrality test consistently verified the reliability of the population structure. Our results provided new insights into the genetic structure in a large hulless barley germplasm resource. Furthermore, we identified an optimal core set of 200 hulless barley accessions with a high maintenance of the genomic and phenomic diversity. The barley core set enabled the reform of the allelic frequency, which potentially boosted power in discovering key QTLs and ultimately facilitating progress in barley breeding.

Abbreviations

AMOVA:

Analysis of molecular variance

CTAB:

Cetyltrimethylammonium bromide

CV:

Coefficient of variation

EN:

Entry-to-nearest

GS:

Genomic selection

GWAS:

Genome-wide association study

MAF:

Minor allele frequency

MCMC:

Markov Chain Monte Carlo

NJ:

Neighbor joining

PCA:

Principal component analysis

PCR:

Polymerase chain reaction

PIC:

Polymorphism information content

QTL:

Quantitative trait locus

SSR:

Simple sequence repeats

References

  1. Takahashi R. The origin and evolution of cultivated barley. Adv Genet. 1955;7:227–66.

    Google Scholar 

  2. Zeng X, Long H, Wang Z, Zhao S, Tang Y, Huang Z, Wang Y, Xu Q, Mao L, Deng G, et al. The draft genome of Tibetan hulless barley reveals adaptive patterns to the high stressful Tibetan Plateau. Proc Natl Acad Sci U S A. 2015;112(4):1095–100.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  3. Sharma KP, Dahal KR, Basta BK. Genetic diversity of Nepalese naked barley and possibility of yield improvement. In: Proceedings of the 2nd national conference on science and technology, 1994, Kathmandu. Kathmandu: Royal Nepal Acad Sci Tech (RONAST); 1994. p. 231–7.

    Google Scholar 

  4. Baniya BK, Dongol DMS, Riley KW. Characterization of Nepalese barley germplasm. Rachis. 1997;16:16–9.

    Google Scholar 

  5. Pandey M, Wagner C, Friedt W, Ordon F. Genetic relatedness and population differentiation of Himalayan hulless barley (Hordeum vulgare L.) landraces inferred with SSRs. Theor Appl Genet. 2006;113(4):715–29.

    Article  CAS  PubMed  Google Scholar 

  6. Xiao Y, Liu H, Wu L, Warburton M, Yan J. Genome-wide association studies in maize: praise and stargaze. Mol Plant. 2017;10:359–74.

    Article  CAS  PubMed  Google Scholar 

  7. Myles S, Peiffer J, Brown PJ, Ersoz ES, Zhang Z, Costich DE, Buckler ES. Association mapping: critical considerations shift from genotyping to experimental design. Plant Cell. 2009;21:2194–202.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  8. Zhu C, Gore M, Buckler ES, Yu J. Status and Prospects of Association Mapping in Plants. Plant Genome. 2008;1(1):5–20.

    Article  CAS  Google Scholar 

  9. Huang X, Zhao Y, Wei X, Li C, Wang A, Zhao Q, Li W, Guo Y, Deng L, Zhu C, et al. Genome-wide association study of flowering time and grain yield traits in a worldwide collection of rice germplasm. Nat Genet. 2012;44:32–9.

    Article  Google Scholar 

  10. Yano K, Yamamoto E, Aya K, Takeuchi H, Lo PC, Hu L, Yamasaki M, Yoshida S, Kitano H, Hirano K, et al. Genome-wide association study using whole-genome sequencing rapidly identifies new genes influencing agronomic traits in rice. Nat Genet. 2016;48(8):927–34.

    Article  CAS  PubMed  Google Scholar 

  11. Zhang X, Warburton ML, Setter T, Liu H, Xue Y, Yang N, Yan J, Xiao Y. Genome-wide association studies of drought-related metabolic changes in maize using an enlarged SNP panel. Theor Appl Genet. 2016;129:1449–63.

    Article  CAS  PubMed  Google Scholar 

  12. Li F, Chen B, Xu K, Wu J, Song W, Bancroft I, Harper AL, Trick M, Liu S, Gao G, et al. Genome-wide association study dissects the genetic architecture of seed weight and seed quality in rapeseed (Brassica napus L.). DNA Res. 2014;

  13. Jiang H, Huang L, Ren X, Chen Y, Zhou X, Xia Y, Huang J, Lei Y, Yan L, Wan L, et al. Diversity characterization and association analysis of agronomic traits in a Chinese peanut (Arachis hypogaea L.) mini-core collection. J Integr Plant Biol. 2014;56(2):159–69.

    Article  PubMed  Google Scholar 

  14. Doyle J, Doyle J. A rapid DNA isolation procedure for small quantities of fresh leaf tissue. Phytochem Bull. 1987;19(1):11–5.

    Google Scholar 

  15. Ramsay L, Macaulay M, Degli Ivanissevich S, MacLean K, Cardle L, Fuller J, Edwards K, Tuvesson S, Morgante M, Massari A. A simple sequence repeat-based linkage map of barley. Genetics. 2000;156(4):1997–2005.

    CAS  PubMed  PubMed Central  Google Scholar 

  16. Huang L, Jiang H, Ren X, Chen Y, Xiao Y, Zhao X, Tang M, Huang J, Upadhyaya HD, Liao B. Abundant Microsatellite Diversity and Oil Content in Wild Arachis Species. PLoS ONE. 2012;7(11):e50002.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  17. Liu K, Muse SV. PowerMarker: an integrated analysis environment for genetic marker analysis. Bioinformatics. 2005;21(9):2128–9.

    Article  CAS  PubMed  Google Scholar 

  18. Yeh F, Boyle T. Population genetic analysis of codominant and dominant markers and quantitative traits. Belg J Bot. 1997;129:157.

    Google Scholar 

  19. The SAS Institute. SAS/STAT User's Guide, Version 8. Cary: SAS Institute; 1999.

    Google Scholar 

  20. Pritchard J, Stephens M, Donnelly P. Inference of population structure using multilocus genotype data. Genetics. 2000;155(2):945–59.

    CAS  PubMed  PubMed Central  Google Scholar 

  21. Evanno G, Regnaut S, Goudet J. Detecting the number of clusters of individuals using the software structure: a simulation study. Mol Ecol. 2005;14(8):2611–20.

    Article  CAS  PubMed  Google Scholar 

  22. Liu K, Goodman M, Muse S, Smith JS, Buckler E, Doebley J. Genetic structure and diversity among maize inbred lines as inferred from DNA microsatellites. Genetics. 2003;165(4):2117–28.

    CAS  PubMed  PubMed Central  Google Scholar 

  23. Tamura K, Dudley J, Nei M, Kumar S. MEGA4: Molecular Evolutionary Genetics Analysis (MEGA) Software Version 4.0. Mol Biol Evol. 2007;24(8):1596–9.

    Article  CAS  PubMed  Google Scholar 

  24. Excoffier L, Laval G, Schneider S. Arlequin version 3.0: an integrated software package for population genetics data analysis. Evol Bioinformatics Online. 2005;1:47–50.

    CAS  Google Scholar 

  25. Weir B, Cockerham C. Estimating F-statistics for the analysis of population structure. Evolution. 1984;38(6):1358–70.

    CAS  PubMed  Google Scholar 

  26. Ewens WJ. The sampling theory of selectively neutral alleles. Theor Popul Biol. 1972;3(1):87–112.

    Article  CAS  PubMed  Google Scholar 

  27. Watterson GA. The homozygosity test of neutrality. Genetics. 1978;88(2):405–17.

    CAS  PubMed  PubMed Central  Google Scholar 

  28. Beukelaer HD, Davenport G: corehunter: Multi-Purpose Core Subset Selection. R package version 3.0.1. 2016.

    Google Scholar 

  29. The HUGO Pan-Asian SNP Consortium. Mapping Human Genetic Diversity in Asia. Science. 2009;326(5959):1541–5.

    Article  Google Scholar 

  30. Xiao Y, Cai D, Yang W, Ye W, Younas M, Wu J, Liu K. Genetic structure and linkage disequilibrium pattern of a rapeseed (Brassica napus L.) association mapping panel revealed by microsatellites. Theor Appl Genet. 2012;125(3):437–47.

    Article  CAS  PubMed  Google Scholar 

  31. Xiao Y, Tong H, Yang X, Xu S, Pan Q, Qiao F, Raihan MS, Luo Y, Liu H, Zhang X, et al. Genome-wide dissection of the maize ear genetic architecture using multiple populations. New Phytol. 2016;210(3):1095–106.

    Article  CAS  PubMed  Google Scholar 

  32. Li H, Peng Z, Yang X, Wang W, Fu J, Wang J, Han Y, Chai Y, Guo T, Yang N, et al. Genome-wide association study dissects the genetic architecture of oil biosynthesis in maize kernels. Nat Genet. 2013;45(1):43–50.

    Article  CAS  PubMed  Google Scholar 

  33. Lipka AE, Kandianis CB, Hudson ME, Yu J, Drnevich J, Bradbury PJ, Gore MA. From association to prediction: statistical methods for the dissection and selection of complex traits in plants. Curr Opin Plant Biol. 2015;24:110–8.

    Article  PubMed  Google Scholar 

  34. Yu X, Li X, Guo T, Zhu C, Wu Y, Mitchell SE, Roozeboom KL, Wang D, Wang ML, Pederson GA, et al. Genomic prediction contributing to a promising global strategy to turbocharge gene banks. Nat Plants. 2016; doi:10.1038/NPLANTS.2016.150.

  35. Visscher PM, Brown MA, McCarthy MI, Yang J. Five years of GWAS discovery. Am J Hum Genet. 2012;90(1):7–24.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  36. Manolio TA, Collins FS, Cox NJ, Goldstein DB, Hindorff LA, Hunter DJ, McCarthy MI, Ramos EM, Cardon LR, Chakravarti A, et al. Finding the missing heritability of complex diseases. Nature. 2009;461(7265):747–53.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

Download references

Acknowledgments

We thank two anonymous reviewers for providing kind suggestions to improve the manuscript.

Funding

This work was supported by the following funding sources: the Financial Special Fund (2015CZZX001), the Tibet Department of Major Projects (Z2016B01N01), and the National Key Laboratory of Special Research (Z2016D03G01).

Availability of data and materials

The data supporting the results of this article are included within the article and its additional files. The materials are accessible by the corresponding author (nima_zhaxi@sina.com) or the National Infrastructure of Plant Germplasm Resource in China (http://www.cgris.net/pt/) upon requested.

Author information

Authors and Affiliations

Authors

Contributions

QX, YW, ZS performed research; ZL, BL analyzed data; HY provided materials; XZ, NT designed the research and XZ wrote the paper. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Nyima Tashi.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Additional files

Additional file 1: Table S1.

Accession information of 562 hulless barley accessions in the present study. (XLSX 57 kb)

Additional file 2: Table S2.

Information of the 73 SSR primers based on the hulless barley reference genome. (XLSX 14 kb)

Additional file 3: Figure S1.

Structural analysis of the 562 barley accessions based on the Bayesian model. a) the LnP(D) plot against the k series; b) the delta k plot against the k series; c) bar-plot of 562 barley accessions representing the membership probability of the inferred subpopulations. At the probability cutoff of 0.65, the total population was divided into two subpopulations, p1 () and p2 (), appending the mixed subpopulation (). (TIFF 1187 kb)

Additional file 4: Figure S2.

Information of the genetic diversity for the SSR markers used in this study. (XLSX 44 kb)

Additional file 5: Table S3.

Boxplots for the comparison of the genetic diversity across the inferred subpopulations. (TIFF 4562 kb)

Additional file 6: Table S4.

Genetic differentiation coefficient (FST) and Nei1972 distance between the pairwise structural inferred subpopulations. (XLSX 8 kb)

Additional file 7: Figure S3.

Genetic distribution of 200 barley core accessions based on the phylogenic and principal component analysis. a) in the phylogenic tree, the clades were colored by structural results, i.e., red and green colors indicate the accessions from the p1 and p2 subpopulations, respectively, while gray to the accessions from mixed population; the clades labeled with black dots refer to the 200 barley core accessions; b) in the PCA plot, the red points refer to the 200 barley core accession. (TIFF 4117 kb)

Additional file 8: Figure S4.

Phenotypic distribution of the four agriculturally important traits in the core and raw barley sets. (TIFF 88 kb)

Additional file 9: Table S5.

Information on the fingerprint identifiers for the 200 barley core accessions. (XLSX 13 kb)

Additional file 10: Figure S5.

Refinement of the population structure. A1-A2, indicates 2 subpopulations at K = 2; B1-B5, indicates 5 subpopulations at K = 5; C1-C7, indicates 7 subpopulations at K = 7. (TIFF 3398 kb)

Additional file 11: Table S6.

Significance of the correlations between the population structure and several clustering factors. The p values were calculated using the contingency table test. (XLSX 8 kb)

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Xu, Q., Zeng, X., Lin, B. et al. A microsatellite diversity analysis and the development of core-set germplasm in a large hulless barley (Hordeum vulgare L.) collection. BMC Genet 18, 102 (2017). https://doi.org/10.1186/s12863-017-0563-x

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12863-017-0563-x

Keywords