The USDA Agriculture Research Service National Plant Germplasm System (NPGS) preserves the largest sorghum germplasm collection in the world, which includes 7,217 accessions from the center of diversity in Ethiopia. The characterization of this exotic germplasm at a genome-wide scale will improve conservation efforts and its utilization in research and breeding programs. Therefore, we phenotyped a representative core set of 374 Ethiopian accessions at two locations for agronomic traits and characterized the genomes.
Using genotyping-by-sequencing, we identified 148,476 single-nucleotide polymorphism (SNP) markers distributed across the entire genome. Over half of the alleles were rare (frequency < 0.05). The genetic profile of each accession was unique (i.e., no duplicates), and the average genetic distance among accessions was 0.70. Based on population structure and cluster analyses, we separated the collection into 11 populations with pairwise FST values ranging from 0.11 to 0.47. In total, 198 accessions (53%) were assigned to one of these populations with an ancestry membership coefficient of larger than 0.60; these covered 90% of the total genomic variation. We characterized these populations based on agronomic and seed compositional traits. We performed a cluster analysis with the sorghum association panel based on 26,026 SNPs and determined that nine of the Ethiopian populations expanded the genetic diversity in the panel. Genome-wide association analysis demonstrated that these low-coverage data and the observed population structure could be employed for the genomic dissection of important phenotypes in this core set of Ethiopian sorghum germplasm.
The NPGS Ethiopian sorghum germplasm is a genetically and phenotypically diverse collection comprising 11 populations with high levels of admixture. Genetic associations with agronomic traits can be used to improve the screening of exotic germplasm for selection of specific populations. We detected many rare alleles, suggesting that this germplasm contains potentially useful undiscovered alleles, but their discovery and characterization will require extensive effort. The genotypic data available for these accessions provide a valuable resource for sorghum breeders and geneticists to effectively improve crops.
Sorghum [Sorghum bicolor (L.) Moench] is the fifth most important grain crop after maize, wheat, rice, and barley . Its role as source of grain and cellulosic-based ethanol is of growing importance. The region of domestication and diversity is located in Northeast Africa, an area that extends from Ethiopia to Sudan , but the crop is now grown in an array of environments, ranging from tropical to temperate regions. Sorghum is a highly diverse crop that has experienced multiple re-domestication processes, resulting in five major races differentiated by inflorescence type . These races are associated with particular environments , and ongoing evolution in these regions by recombination and selection has led to a highly diverse crop.
Ex situ sorghum germplasm collections (i.e., gene banks) have been established in several countries as a resource for breeding programs. These germplasm collections provide genetic variation for crop improvement programs. Today, the largest worldwide sorghum collection is maintained by the USDA Agriculture Research Service National Plant Germplasm System (NPGS) and includes >41,860 accessions from 114 countries. Likewise, the International Crops Research Institute for the Semi-arid Tropics (ICRISAT), located in India, has a sorghum collection of 37,904 accessions from different countries. The National Genebank of China maintains a sorghum collection of over 16,000 accessions, including 12,000 Chinese sorghum landraces. To use these resources, breeding programs must have an understanding of the variation within the germplasm collection . The establishment of core reference sets (i.e., core collections) that represent the genetic diversity of the broader collection provides a mechanism to classify these resources for the scientific community. In this regard, ICRISAT has developed a core collection of 2,247 accessions based on phenotypic information [6, 7]; this collection was later reduced to a mini-core of 242 accessions based on genetic profiles . A core collection of 2,438 accessions was established for the NPGS sorghum collection using origin site information (i.e., passport information ).
The selection of landraces for inclusion in a core set based on passport information alone may not adequately represent the genetic diversity owing to limited passport information for many accessions and/or genetic redundancy (i.e., duplicates). Therefore, the use of a molecular genetics approach is imperative for the conservation and utilization of the collection. However, molecular genetic approaches have been limited to specific samples of sorghum from the NPGS [5, 10–13], and a large portion of the genetic diversity remains uncharacterized in the germplasm collection. These studies confirm that modern sorghum breeding materials are from a narrow genetic base, despite high genetic diversity in the germplasm collection. Recently, a sorghum association panel (SAP) comprising 149 U.S. breeding lines and their 228 progenitor accessions was assembled with the purpose of genetically dissecting economically important traits . Certainly, this panel encompasses a wide range of sorghum genetic diversity that has been exploited in breeding programs, but it is not representative of the vast genetic diversity present in the NPGS collection. Further genomic characterization of a broader range of sorghum diversity is required to both improve conservation efforts and provide new knowledge to promote its integration into sorghum breeding programs.
Genome sequence information and next-generation sequencing technology have made important contributions to the development of new genotyping platforms. Today, genotyping-by-sequencing (GBS)  is arguably the most widely used method to genetically characterize plant germplasm. The large number of single-nucleotide polymorphisms (SNPs) identified by GBS makes it possible to integrate into a single research project analyses of genetic diversity and genome-wide association studies (GWAS) , and is invaluable for the identification of allelic diversity and rare genetic variation as well as functional analyses. In maize, an analysis of 2,815 inbred accessions from NPGS resulted in the identification of 681,257 SNP markers, half of which were classified as rare variants (frequency < 0.05) . Likewise, in sorghum, an analysis of 971 worldwide sorghum accessions (including the ICRISAT mini-core and SAP) resulted in the identification of ~265,000 SNP markers , and an analysis of 1,943 georeferenced sorghum landraces (844 and 1,099 from NPGS and ICRISAT, respectively) was used to identify 404,627 SNP markers . These sorghum genetic diversity studies have substantially improved our understanding of the genetic basis of the domestication process as well as the adaptation of various sorghum races within particular environments. Moreover, they substantiate the use of genotypic and environmental information to predict the phenotypes of accessions. Hence, the phenotypes of accessions within large germplasm collections could be predicted based on prior genotypic and phenotypic analyses of a core subset of accessions.
The Ethiopian germplasm maintained by NPGS contains 7,217 accessions, representing 17% of the sorghum collection . Unfortunately, the geographic coordinates of the collection sites are unknown for 83% of these accessions, limiting the establishment of a representative subset based on passport information alone. An analysis of 137 randomly selected Ethiopian accessions with 20 simple sequence repeat markers revealed abundant genetic diversity divided into two main subpopulations . Four of the five sorghum races and wild forms are present in Ethiopia , and previous phenotypic and genotypic analyses identified high levels of genetic diversity in this germplasm [21, 22]. In addition, the geographic distribution of morphological characteristics appears to be determined by topography and climatic variation . Recently, some Ethiopian accessions from the NPGS were genetically characterized for use in GWAS (135 accessions ; 31 accessions ), but the limited number of accessions and/or the selection approach (i.e., breeding lines and/or ecogeographical regions) may bias estimates of the actual genetic diversity in the collection because these samples were not established to represent the Ethiopian collection. In the current study, we established a core subset of 374 NPGS Ethiopian accessions based on the actual frequency distribution of the known sorghum races in this collection. We phenotyped these accessions for several important agronomic traits at two locations and used GBS to 1) evaluate their genetic and phenotype diversity, 2) determine population structure, 3) study the genetic relationships among accessions (and identify duplicates), 4) determine the genetic relationship between the Ethiopia NPGS collection and the SAP, and 5) explore the potential of this germplasm to elucidate the genetic architecture of quantitative traits using GWAS.
GBS of NPGS Ethiopian germplasm
The Ethiopian germplasm core set examined in this experiment comprised 352 exotic accessions and 22 converted tropical lines that are photoperiod-insensitive and dwarf and were produced by crossing exotic lines and modern US cultivars. This core set represented ~5% of the Ethiopian collection, including 16 sorghum races, of which Durra, Durra-bicolor, and unclassified accessions represented 72% of the accessions (Table 1). We pooled tissue from three seedlings for DNA extraction and analyzed the data as a single GBS sample. We obtained approximately 5×sequencing coverage per accession, resulting in the identification of 417,852 SNPs with an average of 34% missing data. We retained SNPs with <20% missing data and minor allele frequencies of >0.01 (i.e., SNPs present in at least three accessions) to obtain a final sample of 148,476 SNPs for further analyses. These SNPs were distributed along the ten sorghum chromosomes with an average of 1 SNP per 5,200 kb; the majority of SNPs were located in sub-telomeric gene-rich regions. In fact, these SNPs tagged 17,196 annotated genes, while another 4,016 genes were within 5 kbp upstream or downstream of a SNP.
Integrity of NPGS Ethiopian collection
The preservation and effective use of the NPGS Ethiopian collection require an adequate knowledge of its genetic diversity to avoid genetic drift and phenotyping duplicate accessions. The majority of these accessions were collected more than 70 years ago and represent multiple independent germplasm collections; likewise, the seed increase process for this enormous collection (i.e., curatorial management) may contribute to problems, such as material duplication by labeling and logistic errors. We calculated the identity-by-state (IBS) genetic distance among all pairs of accessions using 27,306 unlinked SNPs (r2 < 0.5) that were evenly distributed across the sorghum genome. Remarkably, of the 69,564 pairwise IBS genetic distances, only 14 exceeded 0.90 (Fig. 1), and these represented five groups of 11 closely related accessions. Indeed, no genetic profile was duplicated and the average IBS among accessions was 0.70. The homogeneity of each accession is essential for phenotyping analyses. In the GBS analysis, we detected high levels of inbreeding within accessions; the average heterozygosity (0.11) was not larger than the observed heterozygosity in previous studies of sorghum landraces . Even though this study was based on 5% of the NPGS Ethiopian collection, the results indicate that its assembly and maintain has been adequate capturing and preserving its genetic diversity.
Understanding the population structure of the NPGS Ethiopian collection is imperative to establish new strategies to preserve and utilize its genetic diversity and is a prerequisite for GWAS. We inferred that there are 11 ancestral populations based on 10-fold cross-validation using ADMIXTURE  (Fig. 2a; referred to as populations 1 to 11). In total, we assigned 198 accessions (53%) to one of these populations with an ancestry membership coefficient of greater than 0.60, and we categorized the remaining 47% as having evidence of admixture (Fig. 2b). Indeed, the pairwise fixation index (FST) among these populations ranged from 0.11 to 0.47, indicating a relatively high level of genetic differentiation (Table 2). We determined the effective population size (Ne) based on the linkage disequilibrium (LD) of 58,634 SNPs with MAF > 0.05. We obtained an Ne of 199 individuals, which was comparable to the total number of accessions with an ancestry membership coefficient of greater than 0.60. In fact, these 198 accessions covered 90% of the observed allelic diversity (133,874 SNPs); accordingly, additional phenotyping could be limited to or initiated with these accessions. The conservation of these particular accessions is imperative to maintain the genetic diversity of this core set of the collection.
We used an unrooted neighbor-joining tree to determine the genetic relationships among these populations. The results of this analysis were consistent with the previously determined population structure (Fig. 3), with 11 populations belonging to four main clades. The largest clade included populations 1, 2, 4, 6, and 10, suggesting high genetic relatedness among these populations. Additionally, populations 5, 11, and 9 formed a group. The majority of the accessions from populations 3 and 7 belonged to one common clade, while population 8 belonged to a group with three accessions from population 7. We observed accessions with evidence of admixture throughout the tree, although we detected one large cluster that included one accession from population 4. The genetic relationships among populations provide insight into the history of the observed population structure. We inferred that four founder populations may have subsequently subdivided into 11 populations. We did not detect concordance between our phylogenetic and population structure results and the 14 sorghum races present in the Ethiopian germplasm collection. Durra and the intermediate Durra-bicolor are the most frequent races, and their domestication is associated with the Ethiopian region . Accessions classified as these two races were dispersed among the 11 populations, indicating that the observed population structure was most likely the result of adaptation to different Ethiopian environments. The selection of accessions for phenotype evaluations based on this population structure analysis ensures the presence of the majority of the genetic diversity present in the NPGS Ethiopian germplasm collection.
Based on the allele frequency distribution for this NPGS Ethiopia collection core, 60% of the SNPs were rare (MAF < 0.05) (Fig. 4). A high frequency of rare alleles was also observed in a GBS analysis of 2,815 maize inbred accessions from the NPGS collection , confirming the high accuracy of rare allele detection using GBS. In fact, we detected all of these rare alleles in at least three accessions (MAF > 0.01). We found that the 198 representative accessions (including 90% of the SNPs) had a similar proportion of rare alleles. In the four largest populations (3, 6, 7, and 9), which included more than 25 accessions, we found that an average of 55% of alleles were rare. For populations 1, 2, and 10, which included 11 to 19 accessions, an average of 29% of alleles were rare. Since populations 4, 5, and 8 had fewer than ten accessions, we did not estimate the proportion of rare alleles. The large number of admixture accessions in the collection and the distribution of rare alleles among populations suggest these rare alleles could represent recent introgression from neighboring countries. To compare the diversity between these 11 populations, we determined the total number of private alleles in each population. We identified 10,284 private alleles that characterized these 11 populations. Populations 1 and 2 had the most private alleles (1,995 and 2,264, respectively) and these two populations (30 accessions in total) had 66% of the total allelic diversity (98,322 SNPs). Populations six and nine had 1,286 and 1,353 private alleles, respectively, and had 68% of the total allelic diversity (101,681 SNPs) with a total of 62 accessions. In the 176 accessions with evidence of admixture, we only detected 1,473 private alleles. The NPGS Ethiopian collection could be arranged into these 12 groups (11 populations and an admixture group) using a high-throughput genotyping platform based on these private alleles. The organization of this germplasm collection can improve conservation efforts and phenotypic analyses. For instance, the selection of particular populations and/or subsets of accessions that capture the majority of the allelic diversity will optimize the use of phenotyping resources. Allelic richness could be monitored as an indicator of the conservation status of its genetic diversity .
How much of the genetic diversity present in NPGS Ethiopian germplasm has been used in US breeding programs? The US SAP encompasses a wide range of sorghum germplasm that has been exploited in breeding programs and was recently genetically characterized . We identified all of the common alleles between SAP and this Ethiopian core set. We detected 60,962 common alleles, of which 74% had a MAF of greater than 0.05. In fact, 20% (18,213 SNPs) of the rare alleles in the Ethiopia core set were present in SAP, while 57% (10,382) had a MAF of greater than 0.05. In total, 59% of the 10,284 private alleles that characterize the 11 Ethiopian populations were present in SAP. However, only the private alleles from populations 1 and 2 were adequately represented in SAP (87 and 80%, respectively). In fact, the neighbor-joining cluster analysis of SAP and the 11 Ethiopian populations based on 26,026 unlinked SNPs showed that populations 1 and 2 are completely integrated in SAP (Fig. 5), suggesting these two populations have been used extensively in breeding programs. The other nine Ethiopian populations expanded the genetic diversity of SAP accessions originally from Ethiopia and India. Hence, the NPGS Ethiopian collection is a germplasm resource that could be used to support and increase the genetic diversity in SAP. The integration of particular Ethiopian accessions from populations 3 to 11 into SAP will be useful for the discovery of new agronomically important alleles and to advance the use of this germplasm collection.
Genetic diversity and signatures of selection
The management of germplasm collections is intended to maintain the maximum possible genetic diversity of the crop. We measured the levels of genetic diversity and divergence among populations in the NPGS Ethiopia collection based on LD (as estimated by r2) and FST across the genome. In the Ethiopia collection, we observed an average r2 of ≤0.2 at distances of greater than 6 kb (Additional file 1: Figure S1). This decay in LD was more rapid than the decay observed in previous analyses of 330 landraces from different African countries and the SAP, in which LD extended from 20 to 50 kb . The rapid decay in LD among Ethiopian accessions could be explained by the early establishment of these populations after the domestication of sorghum. Moreover, many accessions had evidence of admixture, which reduces the average length of LD blocks, reflecting a mixture of ancient and new Ethiopian landraces. Certainly, the high number of recombination events is associated with the phenotype diversity of this germplasm, and can be useful for the genomic dissection of agronomically important traits by GWAS; however, deep sequencing coverage may be required to cover all genomic blocks.
The population structure of a plant is determined by abiotic and biotic factors that influence selection . The high genetic differentiation index (FST) between these 11 Ethiopian populations (Table 2) indicates selection in particular environmental conditions. The identification of traits that define each of these populations is imperative to improve the use of this germplasm; agronomic traits that are at high frequencies within a particular population can be screened. In this regard, the genetics underlying the observed population structure may reveal genomic regions that are involved in the selection of particular traits. To identify genomic regions that may be under selection pressure, we used the FST outlier method implemented in BayeScan v.2.1 . This approach distinguishes between loci that diverged via random drift and those that diverged via selection. We detected 99 SNPs that diverged by selection, of which 72 SNPs were consistent with diversifying selection and 27 SNPs were consistent with balancing selection (Fig. 6). We observed these signatures of selection across the whole genome, but particularly in sub-telomeric gene-rich regions. We detected diversifying selection on Tannin1  (FST = 0.62). Tannins have diverse biological functions, including protection against pathogenic bacteria and fungi and predation by herbivorous animals . We detected additional FST outlier SNPs with a broad range of gene ontology annotations, including metabolic process (glycosyl transferase), signal transduction (auxin response factor), DNA damage repair, and retrotransposon proteins. Natural selection contributes more substantially to the structure of Ethiopian sorghum populations than selection by farmers.
The NPGS Ethiopian germplasm collection is the largest group of accessions in the sorghum collection; it includes 11 populations based on a genetic analysis. The phenotypic characterization of these populations is necessary to deeply understand the population structure and improve the use of the collection. We characterized these populations with respect to agronomic traits that might be associated with ecogeographical variation and the agricultural system of Ethiopia as well as other traits of potential interest to the sorghum breeding community.
The phenotypic diversity of Ethiopian germplasm was associated with the observed population structure. In fact, among the 11 Ethiopian populations, we observed differences in 14 evaluated traits, with an average of six groups per trait (Tables 3 and 4). These results confirmed the observed diversifying selection in the genetic analysis since contrasting phenotypes characterize these populations (i.e., high, low, and intermediate trait values). In fact, each population had a unique phenotypic profile, suggesting that they represent different environmental conditions affecting selection. For instance, in population 5, we observed evidence of selection for high biomass and small seed size, and we observed higher levels of tannins contentin population 9. Certainly, the population structure could reflect agro-environmental regions in Ethiopia and/or agronomically important traits, and could be used to establish better strategies to select accessions for phenotype screening.
To improve our understanding of the selection process underlying the observed population structure, we used a principal component analysis (PCA) to identify the relative importance of these traits in shaping phenotype diversity. The first four principal components accounted for approximately 60.5% of the total phenotype diversity, and seed compositional traits had the largest loading coefficients (Table 5). In fact, phenol, protein, tannin, and 3-DOA contents had highly similar degrees and directions for PC1 loading. The high diversity of these four seed compositional traits may be explained by opposing selection pressures driven by humans and natural conditions. For instance, high levels of tannins and phenols are associated with grain mold resistance, but these two compounds are unacceptable in seeds for human consumption . Agronomical traits, such as plant height (PH), panicle length and diameter (PL and PD), and flowering, showed positive loading values for PC2, while the dhurrin content (associated with drought tolerance) had loading effects in PC3. These results indicate that seed compositional traits define Ethiopian sorghum cultivars, while agronomical traits have a secondary role owing to the undeveloped agricultural system in the country. Hence, the germplasm includes valuable alleles for seed compositional traits for sorghum breeding programs. In fact, the high genetic and phenotypic diversity within the Ethiopian germplasm is the result of a balance between both artificial and natural selection in each population.
The Ethiopian NPGS germplasm collection is a highly genetically and phenotypically diverse public germplasm with important applications for the improvement of sorghum breeding programs. Despite the large size of this collection, the use of molecular markers associated with important agronomic traits could facilitate the identification of valuable accessions. Therefore, we explored the potential use of these low-coverage data and the population structure analysis for GWAS. We analyzed seed tannin content and PH, two extensively studied traits [16, 29]. The SNP with the strongest association with tannin content (S4_61667908; p < 8.2−10; Additional file 1: Figure S2a) was the same as the one found in an analysis using the SAP  and explained 13.7% of the genetic variance. This SNP was not in complete LD with the tan1 null allele , as observed in the SAP . Additionally, 22 converted dwarf tropical lines in the panel were used to detect dwarfing genes on chromosomes 9 (dw1), 7 (dw3), and 6 (dw2) (Additional file 1: Figure S3); these loci were introgressed during the conversion of these lines  and explained 27% of the genetic variance. Based on Q–Q plots, both traits demonstrated associations when controlling for population structure and genetic relatedness using a compressed mixed linear model, decreasing the identification of spurious associations (Additional file 1: Figure S4).
To identify novel loci in exotic germplasm, we analyzed PH and flowering time (FL) without including the 22 converted dwarf tropical lines. We detected a significant association between PH and the distal region of chromosome 6 (S6_57260888; p < 8.3−7; Fig. 7). This loci explained 6.8% of the genetic variance while the genetic relatedness among accessions (i.e., pseudo-heritability) accounted for 34.5% of the genetic variance indicating that a large portion of the variation is still uncharacterized. Although we did not identify a candidate gene in the region or overlap with a previously detected QTL, the region deserves further study. We did not detect a significant association for FL using the whole subset or the subset lacking the converted lines. In an analysis of other agronomical traits (e.g., PL, width and ratio, seed size, and cold tolerance), we did not identify significant associations. Indeed, in GWAS for some of these traits using the SAP, significant associations have not been detected [33, 34]. The low heritability and phenotype variation for some of these traits, limited number of accessions, low-coverage data, and complexity of the traits could affect the detection of significant associations.
Based on a GWAS of seed compositional traits, we identified novel loci for protein and fat content, and confirmed a previous result for total phenol content. We obtained similar results for total phenol and tannin contents; SNP S4_61667908 demonstrated the most significant association for these traits (p < 1.4−10; Additional file 1: Figure S2b). A significant correlation between these two traits was also observed in the SAP, in which total phenol content was associated with Tannin1 . A genomic region in chromosome 2 was associated with fat and protein content (Fig. 8a and b). Even though we detected a significant correlation between the fat and protein contents (r2 = 0.54; p < 0.0001), we identified two loci in the GWAS that were 57 kbp apart (Fig. 8c). In fact, we detected the strongest associations for fat and protein content with the SNPs S2_57605132 (p < 7.6−8) and S2_57662361 (p < 1.2−6), respectively, which explained 11.5 and 9.6% of the genetic variance, respectively. Since the genetic relatedness among accessions only accounted for 23.3% (fat) and 21.1% (protein) of the genetic variance other important loci are still unidentified. These SNPs were located within the genes Sb02g023740.1 (unknown protein) and Sb02g023770.1 (cytochrome P450), which have not been associated with particular phenotypes in sorghum or Arabidopsis. Additional phenotype and genotype analyses are required to verify these associations. This Ethiopian core set and its genetic characterization could provide a basis for the discovery of important loci in this exotic germplasm.
At present, only a minor portion of the genetic diversity in the NPGS Ethiopian collection is used in sorghum breeding programs. The large size of this particular collection and the lack of genetic and phenotypic data have contributed to the limited use of this germplasm. Therefore, the genomic characterization presented herein should encourage and facilitate its use by breeders and geneticists worldwide. The combination of this genomic information with large-scale phenotyping efforts could be used to mine novel alleles for the improvement of agronomic traits. The development of new cultivars with novel alleles from Ethiopian collections is a promising strategy to translate this genetic diversity into agriculture benefits.
Untapping the genetic variation in the NPGS Ethiopian germplasm for sorghum improvement is challenging and requires precise strategies. Historically, the integration of tropical exotic germplasm into sorghum breeding programs has relied on the use of adapted exotic germplasm (i.e., converted germplasm) for temperate grain production regions via the introgression of dwarfing and early maturity genes . Today, more than 1,000 exotic sorghum lines are adapted to temperate regions, of which approximately 20% originated in Ethiopia . Nevertheless, the rapid fixation of temperate-adapted alleles in early backcross generations reduces genetic variation in multiple genomic regions; in fact, the genetic variation in chromosome 6 is minimal within converted sorghum germplasm . Therefore, breeding schemes that increase recombination in early generations could be used to develop temperate-adapted germplasm. In this regard, a multi-parent advanced generation intercross design  or the employment of various cycles of recurrent genotypic selection should be adequate to exploit the genetic diversity of Ethiopian germplasm. For instance, an initial phenotypic evaluation (i.e., for a trait of interest) may be completed to identify superior and genetically diverse Ethiopian accessions in the core set and these accessions may be crossed with a temperate-adapted germplasm. The progeny can be systematically intercrossed and/or selfed in tropical regions during early generations and genotyped to maintain temperate-adapted alleles (e.g., related to dwarfing and maturity). Superior families can be evaluated in temperate regions for traits of interest and to fix the temperate-adapted alleles. Indeed, previous phenotypic analyses of these types of populations in Arabidopsis , rice , and wheat  identified transgressive segregation, which is valuable for breeding.
Today, the SAP is widely used for the genomic dissection of multiple agronomically important traits [16, 33, 34, 41]. This diversity panel, which includes worldwide sorghum accessions, captures a large portion of sorghum genetic diversity. However, the limited number of Ethiopian accessions (31) does not adequately represent the genetic diversity in the NPGS Ethiopia germplasm collection (Fig. 5). Moreover, these accessions have lower genetic diversity than the original exotic germplasm, since they are temperate-adapted lines developed by the introgression of common genomic regions associated with photoperiod sensitivity and dwarfing . The expanded use of the NPGS collection requires a size reduction to expedite the screening process. Our genomic characterization clearly indicated that the NPGS Ethiopian germplasm collection could be reduced by 50% while maintaining more than 90% of the total genetic diversity. Certainly, a continuous reduction in GBS data makes it possible to generate SNP data on the whole collection at lower coverage (1×) to establish a precise core set.
Germplasm conservation is imperative for continual agricultural improvements. Studies of the genetic relatedness among all ex situ accessions in germplasm collections and genetic redundancy should improve conservation efforts and practical applications. GBS is rapidly becoming the standard choice for genetic studies and germplasm characterization. Hence, in the near future, multiple independent GBS analyses of sorghum accessions from different germplasm collections will likely be integrated into a meta-analysis to identify duplicates and closely related accessions among worldwide germplasm collections. In this regard, a multi-institutional effort is necessary to systematically explore the genetic diversity of all ex situ sorghum collections to broaden the genetic base of breeding programs worldwide. Indeed, the preservation of redundant germplasm in ex situ collections is a waste of limited resources  and does not promote the goal of improving the germplasm.
The NPGS Ethiopian sorghum germplasm is a highly genetically and phenotypically diverse germplasm collection comprising 11 populations. The genetic diversity of the germplasm is largely unused in sorghum breeding programs. Our results provide sorghum breeders and geneticists worldwide with knowledge and tools to utilize and conserve this germplasm. In fact, these results will facilitate resource optimization by improving experimental designs based on population structure and genetic diversity. Moreover, this genomic characterization is the first step toward the development of a NPGS sorghum core collection based on genetic profiles.
The 7,237 Ethiopian sorghum accessions in the current NPGS collection were grouped according to the GRIN database into the 16 known sorghum races and one unknown category (Table 1). Sorghum race classification was selected as the classifiable variable since the information is available for 72% of the accessions . Subsequently, a core set of 374 accessions was established using a stratified sampling method based on these 17 categories (i.e., random selection within each group/category). The final core set included 352 exotic germplasm accessions and 22 converted tropical lines that are photoperiod-insensitive produced by crossing exotic lines and modern US cultivars. In addition, the lines RTx2911 and SC719 were included in the phenotypic analysis as references .
The 374 accessions and two references were evaluated at Isabela and Mayaguez, Puerto Rico during the short day-length season for 2 consecutive years (November 2013 to March 2014 and October 2014 to February 2015). In both locations, a completely randomized design with plots of 1.8 m in length with 0.9 m between rows was used. Plants were maintained using standard management practices, and weeds were controlled with mechanical tillage and hand hoeing.
FL, PH, PL and PD, the ratio of PL/PD, and seed size were obtained for each accession/plot at both locations. FL was defined as the number of days until 50% of the plants within a plot reached anthesis. PH referred to the average distance from the base of the main stalk (i.e., soil) to the top of the panicle at maturity for two representative plants per accession/plot. PL and PD referred to the distance from the first rachis to the top of the panicle, and to the widest section of the panicle , while RA was defined as PL/PD. Seed size was measured based on the volume of 100 seeds. Indeed, the volume of 100 seeds was significantly correlated with seed length (r2 = 0.81), width (r2 = 0.92), diameter (r2 = 0.78), and area (length × width × diameter; r2 = 0.95) (unpublished observations). Thus, seed size referred to the average volume of 100 seed from three representative panicles per accession.
A compositional analysis of whole seeds was performed using near-infrared (NIR) spectroscopy  in 2013. Briefly, 25–50 g of fresh seeds from each accession was scanned using a FOSS XDS NIR (FOSS North America, Eden Praire, MN, USA) spectrometer with a wavelength of range of 400–2,500 nm, and collection spectra data analyzed using ISIscan software (Version 3.10.05933). Based on near-infrared calibration curves, the percent protein, fat, fiber, ash, total phenols (mg gallic acid equivalent/g), condensed tannins (mg catechin equivalents/g), and 3-DOA (absorbance/ml/g) were estimated for each accession.
Early-season cold tolerance was determined using the average cold germination (G72) for fresh seeds from 2013 according to the procedures recommended by the Association of Official Seed Analysts . The lines PI610727 and RTx436 were included in the early-season cold tolerance evaluation as resistant and susceptible references, respectively . Briefly, average G72 values were recorded by sowing ten seeds at 12 °C in polystyrene Petri dishes lined with filter paper moistened with distilled water. The number of germinated seeds after 72 h was determined visually based on the presence or absence of an emerged radicle. Each accession was replicated three times using a completely randomized design. Seeds were allowed to germinate and grow at 12 °C under full light conditions for 3 weeks.
Leaf dhurrin content is a quantitative measure of the level of pre- and post-flowering drought tolerance . In 2013, five representative plants from each accession were punched at anthesis midway between the blade tip and base from the newest fully developed leaf, excluding the flag leaf. The dhurrin content in the leaf samples was extracted and quantified by high-performance liquid chromatography according to Burke et al. .
The combined data over all locations were used to standardize the 16 traits to a mean of zero and a standard deviation of one for multivariate analysis. PCA was performed using the standardized data with a covariance matrix in SAS. Likewise, analysis of variance (ANOVA) were performed using the subset of 198 accessions assigned to the 11 subpopulations. The linear effect model for such ANOVA was the following: Y = μ + L + P + L x P + e; where Y is the trait, μ is the common effect, L is the location effect, P is the population effect, L x P is the year × location effects and e is the plot to plot variation within populations. The ANOVA for traits with one location data was the following Y = μ + P + e. The means of significant different traits (p <0.05) among Ethiopian population was compared using the Fisher’s least significance difference test.
Genomic DNA was isolated from seedlings of the 374 Ethiopian accessions based on the methods of Guillemaut, Marechal-Drouard . DNA was re-purified using the ZR96 DNA Clean & Concentrator-5 (Zymo Research, Irvine, CA, USA). The GBS library was prepared and sequenced at the Institute for Genomic Diversity (University of Cornell, NY; http://www.biotech.cornell.edu/brc/genomic-diversity-facility) using the restriction enzyme ApeKI for digestion. The library had 384 unique barcodes including 4 random blanks (i.e., water). The GBS library was sequenced on four lanes (depth = 5.0) of an Illumina HiSeq 2500, and the GBS analysis pipeline (Tassel Version: 3.0.147), an extension of the Java program TASSEL , was used to call SNPs from the sequenced library by aligning tags to the BTx623 sorghum genome (Sorghum bicolor v.1.0; ). To ensure that sequencing error were not used in our analysis, we retained SNPs with <20% missing data and minor allele frequencies of >0.01 (i.e., SNPs present in at least three accessions) to obtain a final sample of 148,476 SNPs. The numbers of SNPs within genes and exons, and in the proximity of genes (5 kb up- and downstream), were determined using BEDTools .
Population structure, LD, and effective population size
Population structure and admixed ancestry were assessed using a model-based clustering method implemented in ADMIXTURE . Twenty-seven thousand three hundred and six unlinked SNPs (r2 < 0.5) distributed across the genome were identified with PLINK  and used in the analysis. To determine the actual number of populations, ADMIXTURE was run with a 10-fold cross-validation procedure for K ranging from 1 to 20, and the K value with the lowest standard error was selected .
LD between SNPs for each chromosome was estimated in TASSEL 5.0  using a sliding window of 25 bp. These estimates were based on SNPs with a minor allele frequency of >0.05 and R2 between SNPs with a p-value of <0.001. Effective population size was determined based on LD values from SNPs with a minor allele frequency of >0.05 as implemented in SNeP .
Fixation index and outlier detection
The pairwise fixation index among populations was estimated based on the method of Weir and Cockerham  using the R package Hierfstat . To detect signatures of selection, FST outliers were detected among the 27,306 unlinked SNPs using BayeScan v.2.1 . To reduce the identification of false positives, a 50,000-iteration burn-in period and thinning interval size of 10 were used, and the prior odds for the neutral model was up to 1,000. The prior odd threshold to identify FST outliers was determined using a false discovery rate of 0.001 as implemented in the “plot_bayescan” function in R.
Pairwise genetic distances among the 374 Ethiopian accessions were calculated based on IBS as implemented in TASSEL 5.0  using 27,306 previously identified unlinked SNPs. The resulting matrix was subjected to a clustering analysis using neighbor-joining and visualized using Interactive Tree of Life .
Genetic relationships among the NPGS Ethiopia collection and SAP were investigated by the identification of common SNPs in the SAP (265,487 SNPs ) and NPGS Ethiopia collection (148,475 SNPs). In total, 60,962 common SNPs were detected, and a subset of 26,026 unlinked (r2 < 0.5) SNPs were identified in PLINK  for further analysis. Pairwise genetic distances among these 730 lines (374 and 356 from NPGS Ethiopia collection and SAP, respectively) were calculated based on IBS as implemented in TASSEL 5.0 , and the resulting matrix was subjected to a clustering analysis using neighbor-joining and visualized in Interactive Tree of Life .
The presence of population structure and familial relationships among accessions must be adequately accounted for in the GWAS model approach to reduce spurious associations . The association analysis was performed with a compressed mixed linear model  implemented in the R package GAPIT  using 58,634 SNPs with a MAF > 0.5. Kinship was calculated as described by VanRaden , and a co-ancestry matrix from ADMIXTURE was included as a covariate in GAPIT to reduce spurious associations. Log Q–Q plots of p-values were examined to determine how well the models accounted for population structure and familial relatedness. The phenotypic data of PD, RA, and tannin content were normalized prior to the association analysis. Significant associations were determined for each trait using a false discovery rate-adjusted p-value of >0.05  as implemented in GAPIT. Manhattan and Q–Q plots were visualized using the R package qqman .
Average cold germination
Genome Association and Prediction Integrated Tool
Genome-wide association studies
Identity by state
International Crops Research Institute for the Semi-arid Tropics
Ali ML, Rajewski JF, Baenziger PS, Gill KS, Eskridge KM, Dweikat I. Assessment of genetic diversity and relationship among a collection of US sweet sorghum germplasm by SSR markers. Mol Breeding. 2008;21(4):497–509. doi:10.1007/s11032-007-9149-z.
Cuevas HE, Prom LK. Assessment of molecular diversity and population structure of the Ethiopian sorghum [Sorghum bicolor (L.) Moench] germplasm collection maintained by the USDA-ARS National Plant Germplasm Systems using SSR markers. Genet Resour Crop Evol. 2013;60:1817–30.
Wang ML, Chengsong Z, Barkley NA, Chen Z, Erpelding JE, Murray SC, et al. Genetic diversity and population structure analysis of accessions in the US historic sweet sorghum collection. Theor Appl Genet. 2009;120:13–23. doi:10.1007/s00122-009-1155-6.
Morris GP, Ramu P, Deshpande SP, Hash CT, Shah T, Upadhyaya HD, et al. Population genomic and genome-wide association studies of agroclimatic traits in sorghum. Proc Natl Acad Sci U S A. 2013;110(2):453–8.
Teshome A, Baum BR, Fahrig L, Torrance JK, Arnason TJ, Lambert JD. Sorghum [Sorghum bicolor (L.) Moench] landrace variation and classification in north Shewa and south Welo, Ethiopia. Euphytica. 1997;97(3):255–63.
Ayana A, Bekele E. Geographical patterns of morphological variation in sorghum (Sorghum bicolor (L.) Moench) germplasm from Ethiopia and Eritrea: qualitative characters. Hereditas. 1998;129(3):195–205.
Zhang D, Kong W, Robertson J, Goff VH, Epps E, Kerr A, et al. Genetic analysis of inflorescence and plant height components in sorghum (Panicoidae) and comparative genetics with rice (Oryzoidae). BMC Plant Biol. 2015;15:107.
Zhang D, Li JP, Compton RO, Robertson J, Goff VH, Epps E, et al. Comparative genetics of seed size traits in divergent cereal lineages represented by Sorghum (Panicoidae) and rice (Oryzoidae). G3-Genes Genom Genet. 2015;5(6):1117–28.
Bandillo N, Raghavan C, Muyco PA, Sevilla MAL, Lobina IT, Dilla-Ermita CJ, et al. Multi-parent advanced generation inter-cross (MAGIC) populations in rice: progress and potential for genetics research and breeding. Rice. 2013;6:11.
Mackay IJ, Bansept-Basler P, Barber T, Bentley AR, Cockram J, Gosman N, et al. An eight-parent multiparent advanced generation inter-cross population for winter-sown wheat: creation, properties, and validation. G3-Genes Genom Genet. 2014;4(9):1603–10.
Dykes L, Hoffmann L, Portillo-Rodriguez O, Rooney WL, Rooney LW. Prediction of total phenols, condensed tannins, and 3-deoxyanthocyanidins in sorghum grain using near-infrared (NIR) spectroscopy. J Cereal Sci. 2014;60(1):138–42.
Burke JJ, Chen JP, Burow G, Mechref Y, Rosenow D, Payton P, et al. Leaf dhurrin content is a quantitative measure of the level of Pre- and postflowering drought tolerance in sorghum. Crop Sci. 2013;53(3):1056–65.
Barbato M, Orozco-terWengel P, Tapio M, Bruford MW. SNeP: a tool to estimate trends in recent effective population size trajectories using genome-wide SNP data. Front Genet. 2015;6:109. doi:10.3389/fgene.2015.00109.
The authors acknowledge the Sorghum and Millet Germplasm Committee for support with the GBS analysis of the Ethiopian sorghum germplasm.
Availability of data and materials
The seeds of accessions used in this paper are free available through USDA-ARS NPGS system (https://npgsweb.ars-grin.gov/gringlobal/search.aspx). Phenotype data supporting this manuscript is available as Additional file 2, and genotype dataset are available by request to corresponding author.
HEC conceived of the study, collected the phenotype data, analyzed the phenotype and molecular information, and drafted the manuscript. GRV collected the phenotype data and performed DNA extraction. WLR and LH carried out the compositional analysis of the seeds and helped draft the manuscript. All authors read and approved the final manuscript.
The authors declare that they have no competing interests.
Authors and Affiliations
USDA-ARS, Tropical Agriculture Research Station, 2200 Pedro Albizu Campos Avenue, Mayaguez, 00680, Puerto Rico
Hugo E. Cuevas & Giseiry Rosa-Valentin
USDA-ARS, Cropping Systems Research Laboratory, 3810 4th Street, Lubbock, TX, 79415, USA
Chad M. Hayes
Departments of Soil & Crop Sciences, Texas A&M University, College Station, TX, 77843-2474, USA
Linkage disequilibrium curves for the 10 sorghum chromosomes among 374 NPGS Ethiopian sorghum accessions using 58,634 SNPs with a MAF > 0.5. Figure S2. GWAS for tannins and phenols content using core set of 374 Ethiopian accessions from NPGS germplasm collection. (a) GWAS for tannin content. (b) GWAS for total phenol content. Figure S3. GWAS for plant height using a subset of 374 Ethiopian accessions from NPGS germplasm collection. Figure S4. QQ-plots from GWAS using a subset of 374 Ethiopian accessions from NPGS germplasm collection. (a) Tannin content, (b) Plant height using 374 accessions, (c) Plant height using 352 accessions, (d) Total phenol content, (e) Flowering using 352 accessions, (f) Fat content, (g) Protein content. (DOCX 638 kb)
GRIN data base information for 374 sorghum accessions from the U.S. National Plant Germplasm System genetic characterized by genotype-by-sequence analysis (GBS). Table S2. Coancestry coefficient matrix of 374 Ethiopian sorghum accessions from NPGS germplasm collection. Table S3. Phenotype analysis of 374 sorghum accessions from the U.S. National Plant Germplasm System evaluated at Isabela and Mayaguez, Puerto Rico in 2013 and 2014, respectively. Table S4. IBD genetic distance matrix for 374 Ethiopian accessions from NPGS using 27,306 SNPs. Table S5. Kinship Matrix for 374 Ethiopian accessions as decribed by VanRaden  using 58,634 SNPs. (XLSX 3336 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.
Cuevas, H.E., Rosa-Valentin, G., Hayes, C.M. et al. Genomic characterization of a core set of the USDA-NPGS Ethiopian sorghum germplasm collection: implications for germplasm conservation, evaluation, and utilization in crop improvement.
BMC Genomics18, 108 (2017). https://doi.org/10.1186/s12864-016-3475-7