- Research article
- Open Access
Assessment of genetic diversity and yield performance in Jordanian barley (Hordeum vulgare L.) landraces grown under Rainfed conditions
BMC Plant Biologyvolume 17, Article number: 191 (2017)
Barley (Hordeum vulgare L.) is a major cereal crop, which is cultivated under variable environmental conditions and abiotic stresses in marginal areas around the globe. In this study, we evaluated 150 Jordanian landraces obtained from ICARDA Gene Bank and four local checks for yield and yield components related-traits in two locations across Jordan for three growing seasons under rainfed conditions. The study aims to identify superior Jordanian barley genotypes under dry conditions, to understand the genotype × environment (G × E) interactions, to analyze stability parameters and to identify markers associated with yield and yield components under rainfed conditions.
The barley accessions exhibited significant variation for all traits studied. Three accessions with high yield, cultivar superiority and stability under specific environments were identified with accession G69 is the highest yielding and superior for Madaba and overall environments and G144 is the highest yielding at Ramtha. Accession G123 was high yielding in all environments and was stable across different environments. At the genetic level, the Jordanian landraces were found to be diverse with a clustering that was based on row-type. The GWAS analysis identified 77 significant markers-traits associations for multiple traits including grain yield (GY) with three significant QTLs located at 1H, 2H and 7H, which seem important for dry environments.
Utilizing Jordanian barley landraces can effectively improve and adapt the current barley cultivars for cultivation under environmental stresses in dry regions. Utilization of markers associated with important agronomical traits and their incorporation in breeding using marker assisted selection can improve barley tolerance to drought stress.
The domestication of barley (Hordeum vulgare L.) took place prior to 10,000 B.C in the Near East region known as the “Fertile Crescent”, which extends from the southern parts of Jordan to southeastern Turkey to western Iran . The first results of the domestication processes are populations known as “landraces”. Farmers had a crucial role in this crop evolution process and in the maintenance of genetic variation of landraces in their farming environment . Landraces usually have limited geographic range, are diverse within particular types, are adapted to local conditions and are able to tolerate drought, diseases and pests . In the case of barley, landraces are common in marginal, low-input and drought-prone environments in the Fertile Crescent and in other types of stressful environments [1, 4]. In addition to their adaption to harsh environmental conditions, they are popular among farmers for their high feed quality of both grain and straw . In Jordan, barley is cultivated in areas that receive an annual rainfall between 150 and 300 mm. These areas face many challenges that are increasing because of climate change, which already resulted in an increased frequency of dry years. Barley landraces domesticated and evolved in the dry areas of the Fertile Crescent are expected to be adapted to growth-limiting factors as this adaptation contributes to their yield stability.
Molecular markers technology offers the possibility to increasing the efficiency and precision of selection in plant breeding program, especially in QTL identification, which determines the location of genes on chromosomes. Therefore, the identification of markers closely linked to important agronomic traits will allow their utilization in marker-assisted selection and thus increase the efficiency of selection in breeding programs [6, 7]. A genome wide association study (GWAS) is a QTL mapping approach that makes use of natural populations, landraces, collection of cultivars released over years and genotypes with little or no pedigree information . It is used to detect marker-trait relationships based on linkage disequilibrium (LD), considering that associated marker can be utilized in future marker assisted selection. Such approach was used for the identification of novel QTLs for selected traits with practical implications in breeding programs [9,10,11]. For instance, two SSR markers Ebmac415 (chromosome 2H) and GMS21 (chromosome 1H) were associated with increased GY by 19% and 13%, respectively, in most favorable environments for H. spontaneum . Mickelson et al.  detected the regions around markers acat466 and acag135 on chromosome 3H, acag273 on chromosome 4H, TB2122 on chromosome 5H, and acgc132 and acgt517 on chromosome 6H that appear to be related with nitrogen metabolism in barley.
The emergence of new high throughput genotyping platforms enabled the implementation of GWAS in barley [14, 15]). Using GWAS, a natural variant of the barley homolog of Antirrhinum CENTRORADIALIS (HvCEN) was identified in a diverse set of barley germplasm and was found to control flowering and environmental adaptation . In another study, GWAS was used to detect QTL for heading date, plant height, thousand grain weight, starch content and crude protein content in a diverse collection of 224 spring barleys of worldwide origin . GWAS in a large panel of 615 barley cultivars using a mixed linear model identified significant associations for sixteen morphologic and nine agronomic traits . GWAS was applied successfully on a large panel consisting of 2417 accessions using the barley iSelect 9 k SNP assay to dissect hull cover, heading date and spike morphology .
The main objectives of this study were to evaluate the performance of a sample of Jordanian barley landraces accessions under rainfed and drought-prone conditions in Jordan and to identify molecular markers associated with their performance under such conditions. In addition, the tested landraces were evaluated for major agronomic traits in order to identify superior and stable lines with improved adaptation to rainfed conditions. Molecular analysis was carried out to study genetic diversity and population structure and to identify genetic markers associated with high grain yield and yield components under dry conditions by using 9 k Illumina iSelect SNP assay and the GWAS approach.
The material used in this study consists of 150 accessions of barley landraces collected in Jordan and held in the Gene Bank of the International Center for Agricultural Research in the Dry Areas (ICARDA) (Additional file 1). These accessions were collected from a wide range of geographic areas across the country. Four local checks were included in the study, namely Rum (a 6-row barley cultivar developed by CIMMYT and released in 1986), Baladi (a commonly cultivated 2-row barley landrace in Jordan), Yarmouk (a 2-row barley cultivar selected from ICARDA material and released in 2004) and Mutah (a 2-row barley cultivar selected from ICARDA material and released in 2004). The seeds of local checks were obtained from the National Center for Agricultural Research and Extension (NCARE) / Jordan.
A total of six field trials were carried out over three growing seasons (2008–2009, 2009–2010 and 2010–2011) in farmers’ fields at Madaba (with around 350 mm long term average rainfall) and Ramtha (with around 217 mm long term average rainfall). The field trial locations are described further in Additional file 2. All field trials were in accordance with the local legislations under the direct supervision of NCARE/Jordan.
The experimental design was a partially replicated (p-rep) row-column design with 20 columns and 10 rows. The 150 accessions were replicated once, while the four checks (Rum, Baladi, Yarmouk and Mutah) were replicated 13, 13, 12 and 12 times, respectively, which resulted in 200 plots for each field trial. Plot size was 9.75 m2 (6 rows × 6.5 m length × 0.25 m distance between rows) planted at a seeding rate of 100 kg ha−1. The PRDiGGer software (http://nswdpibiom.org/austatgen/software/) was used to randomize the entries and a different randomization was used for each of the six year x location combinations. The software uses a user’s specified blocking sequence to distribute replicated treatments in a balanced way across the design to provide an unbiased estimate of the error variance.
The materials were sown on the 1st week of December in each growing season and grown under rainfed conditions. At each location, experimental plots were managed following the standard agricultural practices including weed control (hand weeding and herbicides against broad-leaved weeds) and pesticide use against main pathogens. Data on precipitation during the growing seasons were obtained from meteorological stations nearby the experimental sites (Additional file 2).
Field data collection
The traits recorded on plot basis were: grain yield (GY in g/m2), as average grain weight of two one m2 sample randomly taken from each plot, measured after threshing and cleaning; biomass yield (BY in g/m2), as average straw and grain weight of two “one-m2” sample randomly taken from each plot; straw yield (StY in g/m2) obtained by subtracting GY from BY; harvest index (HI as %) calculated by dividing GY by multiplied by 100; thousand kernel weight (TKW in g), as the weight of 200 kernels taken from the bulk seed of the plot and multiplied by 5; spike weight (SW in g), as the weight of five randomly selected spikes in the plot; kernels per spike (K/S): as the number of kernels/spike counted on five randomly selected spikes in the plot; spike length (SL in cm), as the length measured from the base to the top of the spike excluding the awns. The following traits were recorded as an average of three plants randomly selected within each plot: plant height (PH in cm): measured at maturity from the ground level to the top of the spike excluding awns; awn length (AL in cm), as the length from the top of the spike to the longest awn; peduncle extrusion (PEX in cm), as the length of the section of the last internode from the ligule of the flag leaf to the base of spike. Negative values of PEX indicate that the spike, is at least partially, inside the boot; peduncle length (PL in cm), as the length from last node to the base of spike.
Analysis of data from the un-replicated trials to estimate genotypic performance of the entries was carried out as described in Singh et al. , and adapted for the p-rep design in rows and columns. The spatial analysis procedure screens the best out of the nine spatial patterns described in terms of local fertility trends and spatially correlated plot errors and formed as all possible combinations of three fertility trends: fixed linear trend along rows (i.e. in column number), random cubic spline in column number, and no trend, and three plot-error structures: independent plot-errors, first order auto-regressive errors along rows, and first order auto-regressive error along rows and along columns). The selection of the most fitted model was done using Akaike information criteria and procedure described in Rollins et al. .
The best linear unbiased estimates (BLUE’s) of the genotype means were obtained for each environment and the trait as the spatial pattern are unique for a field in a given year, and used for assessing the genotypes performance specific to locations and stability over years. For each trial, the field heterogeneity was assessed by coefficient of variation as experimental error standard deviation divided by trial mean × 100, and statistical significance of genotypes effects was assessed by a chi-square based test, called Wald-test, when the genotypes effects were assumed fixed. For each individual trail, heritability in broad sense for each field trial was computed as genotypic variance component estimate divided by sum of genotypic variance component estimate and experimental error variance estimate, in which case the genotypes effects were assumed random as is the usual procedure.
BLUEs were used to examine the distribution of genotypes performance at individual trials using boxplots and using means and stability across location and overall the environments. The static stability of a genotype was obtained as coefficient of variation (CV) of its means (BLUEs) across the environments (CV% = standard deviation of the genotype means across environments divided by overall mean of the genotype × 100 . An index, P, to measure superiority of a genotype was introduced by Lin and Binns , which integrates the genetic effect and G × E interaction in comparison with an ideal genotype achieving the maximum yield in each environment. Since P is at variance scale and to make comparable with the CV measure of stability, we have obtained a standardized coefficient of cultivar superiority (SCCS) for a genotype, which was computed as SCCS = square-root of P divided by mean of the genotype × 100. All the above statistical analyses were carried out using GenStat Edition 15 . The codes for various analyses were organized as separate program files.
For environmental grouping based on location and overall, Restricted Maximum Likelihood Method was used to estimate variances of the genotype × site × year model. For this purpose, analysis of data was carried out based on augmented incomplete block design where contrasts of replicated checks and new un-replicated entries were estimated simultaneously assuming random genotypic values. Thereafter, heritability for each environmental grouping was computed as genotypic variance component estimate divided by sum of genotypic variance component estimate, genotype × year interaction, genotype × site interaction, genotype × year × site interaction and experimental error variance estimate, in which case the genotypes effects were assumed random as is the usual procedure.
The GGE biplot software  was used to exhibit specific adaptation and stability of the genotypes. Location-year combinations were considered one distinct environment as suggested by Costa and Bollero  and Okuyama et al. . The GGE-biplot analysis of tested accessions with the distinct environments was carried out using the GGE biplot software (http://www.ggebiplot.com) for GY.
Genotyping and molecular analysis
For genomic DNA (gDNA) extraction, barley plants were grown in a growth chamber under controlled conditions with temperature range between 22 and 24 °C and short-day photoperiod (12 h light). Total gDNA was extracted from leaf tissue collected from three-week-old plant of each accession using the CTAB (cetyltrimethylammonium bromide) method as described in Winnepenninckx et al. .
For genotyping, a set of 7864 independent, high confidence and gene-based SNPs markers, which are incorporated into the Illumina Infinium iSelect 9 k SNP barley array were used to genotype the accessions as described previously . The genotyping assay was performed by the TraitsGenetics GmbH. (Gatersleben, Germany). Markers with allele frequency < 5% or missing data >10% were removed from further analyses.
Cluster analysis, population structure and linkage disequilibrium
For genetic diversity assessment, SNP marker data was used to generate an un-rooted tree by using a neighbor-joining (NJ) algorithm implemented in DARwin software . The genetic structure (Q) of the 154 barley lines was analyzed using 120 SNP markers distributed across the barley genome (a minimum distance of 10 cM between markers on the same chromosome). A clustering method based on a Bayesian model  and implemented in the STRUCTURE version 2.3.3 algorithm (available from http://web.stanford.edu/group/pritchardlab/structure.html) was employed to determine population structure. Both the length of burn-in period and the number of iterations were set at 100,000 iterations with k value in the range of 1–12. To reach the appropriate k value, the estimated delta k value was calculated for each estimated k according to Evanno et al. .
The LD between markers was estimated using a squared allele frequency correlation (R2) as described by Hedrick  and Weir  and only SNP markers with known chromosomal position were used in the estimation of LD. Linkage disequilibrium statistics were calculated per chromosome and across all chromosomes (Inter- and Intra-chromosomal LD), while the decline of LD with genetic map distance was evaluated by plotting R2 against the distance (cM) between markers and fitting a decay curve using the square root transformation of the equation described by Andreescu et al. . The critical value at 95% of R2 values between pairs of unlinked loci was 0.197. The 2nd degree-loess (local regression) smoothing was calculated and plotted using SigmaPlot V.11 software (http://www.sigmaplot.co.uk/products/sigmaplot/produpdates/prod-updates1.php).
Association between SNP markers and yield and yield components traits from selected individual environments, two locations and overall environments with heritability values ≥5 was analyzed using a mixed-model approach (MLM) to control type I errors, accounting for genetic relatedness or kinship (K). The K matrix is a measure of relative kinship and quantifies the probability that two homologous genes are identical by descent  and it was generated within TASSEL 5 (; http://www.maizegenetics.net/tassel). The analysis was performed with and without taking into account the first three principle components as a covariate (Q matrix) [33, 34]) in TASSEL 5 and only markers that confirmed with both MLM and MLM-Q models were considered as true associations. In order to consider the dependency of markers due to LD, we estimated the significant cutoff for P-value using the method described in Li et al. . The method evaluate the effective number of independent tests from large number of dependent markers based on their LD values to determine the accurate significance threshold. However, associations with P < 0.001 between SNP markers and traits were also considered as suggestive QTL. The map positions of the associated markers were based on previously described map reported in Comadran et al. . Physical positions were optioned from “Morex_2016” map available at https://triticeaetoolbox.org/barley.
Performance of the Jordanian barley landraces under rainfed environments
The rainfall data during the three growing seasons were collected from the nearest meteorological stations. The rainfall amounts received in the six field trials were below the long-term average of the respective locations (Additional file 2). The recorded weather data indicates the presence of variations in rainfall patterns across different locations and years indicating the prevalence of drought conditions at tested sites.
Combined ANOVA showed highly significant differences for location × year × genotype interactions with 4.05% of the total sum of squares attributed to genotypic effects (Additional file 3). The best spatial model screened for GY (Additional file 4) shows that out of six trials, linear trends were present in field along the rows of the layout of the trial at Madaba (2010/2011) and Ramtha (2009/2010) while auto-correlation was found in the plot-errors along rows alone at Madaba, 2008/2009 and along both directions in case of Ramtha 2010/2011. No spatial patterns could be gauzed through the nine models used, although the experimental fields showed a high degree of heterogeneity, where the field coefficient of variation varied in the range 18–52%. The Walt test assuming genotype effects fixed showed significant genotypic differences in all the trials (P < 0.05) except for Ramtha 2008/2009 (P = 0.10). However, when genotype effects were assumed random, the broad sense heritability estimates varied from 5 to 45% in four trials while zero in the remaining trials (Additional file 4). In the latter trials the estimates of variance components were zero using restricted maximum likelihood (REML) method. This may be due to lack of desired replications for the test and the check entries. Distributions of BLUEs of genotype means for GY for each trial are summarized as boxplots in Fig. 1. Of all the 6 trials, the highest means and spread of the BLUEs were observed for 2009/2010 at Madaba while the lower mean levels were observed for Madaba in 2008/2009.
On the another hand, environmental grouping based on location and overall was considered and the Restricted Maximum Likelihood Method was used to calculate their heritability values, which was 20.26% for overall and varied from 8.2 to 12.4% in Madaba and Ramtha, respectively (Additional file 3).
The BLUEs were averaged over the years at each locations and also overall the trials and are presented in Table 1 (full list of genotypes and their statistics are provided in Additional file 5) for top fifteen genotypes along with their measures of stability, CV% and genotype superiority (SCCS). Denoting the genotypes/accessions 1–154 by G1 – G154, a spectacular accession was G69, which was the highest yielding overall (259 g/m2) and at Madaba (389 g/m2). Based on the cultivar superiority index (SCCS), G69 is also ranking first in overall and Madaba (Table 1). The highest yielding genotype at Ramtha was G144 with GY mean value of 195 g/m2. Accession G123 was identified within the top 15 genotypes at Madaba, Ramtha and overall and it was stable at Madaba and overall when compared with Ramtha where it was considerably less stable (Table 1). There is a very high, though negative, correlation between mean and the superiority index SCCS (Additional file 5), which makes SCCS redundant for any additional advantage of the genotype beyond mean over the environments. Correlation between mean and CV%, as a stability measure, was found much lower, ranging from −0.35 to 0.07, and significant at Madaba. This provides a scope for identifying genotypes which are high yielding (or superior) as well as stable.
The broad and specific adaptation of accessions to tested environments for GY was studied by using the GGE biplot analysis based on location-year combinations (Fig. 2). Using the means in each location for available years and considering all the genotypes, the first two principal components (PC) accounted for a total of 71.9% variation in both locations. In Madaba, the GGE biplot analysis showed that the growing seasons were divided into two sections where the 2008/2009 and 2010/2011 growing seasons fall in the same section and the 2009/2010 season was in another section and the tested accessions were scattered into nine sections (Fig. 2a). G25 had the highest GY in Madaba in 2008/2009 and in Madaba in 2010/2011, while G69 had the highest GY in Madaba in 2009/2010 (Fig. 1a). In Ramtha, the GGE biplot showed that the growing seasons were divided into three sections and the accessions into nine sections. G34 had the highest GY in Ramtha in 2008/2009, while G8 had the highest GY in Ramtha in 2009/2010 and G123 (Fig. 2b). The biplot analysis results were in general agreement with the BLUEs values in terms of entries performance across environments.
Association between the agronomic traits
Grain yield showed significant (P ≤ 0.01) correlation with BY and StY in both locations (Madaba and Ramtha) and overall environment (Table 2). TKW was significantly correlated at Ramtha but not at Madaba. The other traits that showed significant correlations with GY were PH at both locations and overall; with PEX and AL at Madaba and under overall (Table 2). Correlation heat maps are given for the individual trials for all possible variables where data were recorded (Additional file 6). In general, GY showed a strong positive correlation with BY and StY across different trials.
For the GWAS analysis, filtering the 9 K SNP markers data resulted in 4956 polymorphic SNPs. Chromosome 5H had the largest number of markers (629), while chromosome 4H carried the smallest number with 253 SNPs. Those SNPs spanned a total genetic distance of 989.9 cM with an average of one marker every 0.3 cM (data not shown). The population structure resulted in two distinct sub-populations (Additional file 7). The first sub-population involved 86 genotypes with Fst value equal to 0.67 while the second one has the rest (68 genotypes) including the four checks with Fst value of 0.29. Using DARwin software, an unrooted NJ tree was built describing the relatedness between Jordanian accessions (Fig. 3). The clustering analysis showed that the panel was diverse and separated into several clades where a large one contained the majority of “2- row” accessions (63) and 12 “6-row” accessions that were highly related. The rest of 6-row barley accessions (29 were identified in distinct clades that included few 2-row accessions. Other subgroups consisted exclusively from 2-row accessions with two distinctive subgroups with the first one including the Baladi check and the second including Mutah and Yarmouk cultivars (Fig. 3).
For LD decay analysis, the LD decayed sharply to an R2 value of about 0.2 in less than 2 cM (Additional file 8). The effective number of independent markers following Li et al.  was 431.2 (~14.7%). This adjustment shift the Bonferroni correction from 1.7 × 10−5 to 1.6 × 10−4. The GWAS analysis revealed eleven marker-trait associations (MTA) significantly associated with yield and yield components across the tested environments with heritability values ≥5 (Table 3) and another 77 suggestive associations at P < 0.001 (Additional file 9). QQ plots that compare the expected and the observed P values for all traits are shown in Additional file 10. Those eleven MTAs represents eight different QTL of which the marker SCRI_RS_1971 on chromosome 5H (168.5 cM) showed association with GY in Madaba 2008/2009 (Table 3). Another interesting QTL located on chromosome 7H (1.6 cM) was associated with K/S in Madaba (Table 3) and K/S in overall and with GY Madaba 2008/2009 environments (Additional file 9). Another important QTL on chromosome 1H also showed association with BY and SL in Madaba 2008/2009 and Madaba 2010/2011, respectively (Table 3). Another QTL was also identified at 1H (133 cM) that was associated with many yield components traits in overall environments (Additional file 9). The other QTLs were found on chromosomes 2H, 4H, 5H and 7H for, “BY and HI”, TKWT, SW and “PH, StY and SL”, respectively (Table 3). Of those, QTL on 2H (~106 cM), which was significantly associated with BY and HI in Madaba 2010/2011and with HI in Ramtha 2009/2010 and, was also associated with other yield components traits across different environments but at lower P-value (Additional file 9). The remaining suggestive MTAs were distributed on all chromosomes and confirmed many of the significant QTLs.
The knowledge of genetic diversity available in germplasm collections or breeding material helps the breeders to plan their programs for specific environments using targeted traits and molecular markers. Hence, this study was conducted to evaluate 150 Jordanian barley landraces and 4 improved varieties used as checks, in different dry areas of Jordan. Significant differences between accessions in yield and yield components were identified across different environments. This indicates the presence of variability in the agro-morphological traits within these accessions, which provide ample scope for selecting superior accessions by plant breeders. This is somehow expected, since Jordan is considered a part of the center of origin for wheat and barley [36, 37]. Similar results were observed by Ceccarelli et al. , who revealed significant genetic variation among 10 traits for 70 barley landraces, 60 of which were collected from Syria and 10 from Jordan. Therefore, the landraces analyzed in this study can be considered as a reservoir of genes that plant breeders need in their barley-breeding program.
Accessions that ranked first for a number of traits were significantly superior to the check varieties under rainfed conditions. The range for GY between accessions varied with location and the mean GY was in the range of 130–389 g/m2 at Madaba and 71–195 g/m2 at Ramtha (Additional file 5). Furthermore, some accessions performed well in some environments where they had the highest value for one trait but the lowest value for other traits (data not shown). Similar results were found by Rodriguez et al. , who found that variations was attributed mainly to environment effect (67.1%), genotypes effect (5.3%) and G × E interaction (26.7%). This is expected in such very dry sites, as rainfall (amounts and its distribution) and temperature differs from year to year in the same location. Similar results were obtained by Rizza et al. , Sudaric et al.  and Mohammadi and Amri . Thus, to reduce the magnitude of the interaction it was suggested by Mohammadi et al.  to divide the testing area into sub-regions to avoid inconsistent environmental conditions prevailed across locations.
In this study, accession G69 gave the highest yield in Madaba and overall environments (Table 1). This is somehow expected knowing that G69 was collected from a semi-humid area in northern parts of Jordan where annual rainfall is higher than 350 mm (Additional file 1). Accession G139 was highly stable and performed well in Ramtha location indicating the suitability of this accession to arid regions in north of Jordan. Accession G123 showed good performance across different environments and was highly stable in Madaba and overall indicating its suitability for dry environments (Table 1). All of these good performing accessions were of 2-row barley type, which are known to be adapted to drought conditions. The “Baladi” traditional local barley landrace in Jordan and in neighboring countries with similar drought conditions have a 2-row type, which is preferred by local farmers. The GGE biplots provided further support for the adaptation of the identified superior genotypes to specific locations (Fig 2).
In this study, SNP markers were used to analyze the level of diversity among and within Jordanian barley landraces. A similar approach was used previously to study genetic diversity in barley lines from different origins [44, 10]. Genetic variation among Jordanian accessions identified to be large and several subgroups were observed based on row-type. A subgroup included the majority of accessions with 2-row barley type that show high level of relatedness, while distal subgroups included the 6-row accessions and distinct 2-row accessions like the Baladi subgroup (Fig. 3). On the other hand, the population structure analysis identified two subgroups that where considered highly structured reflecting the impact of material origin and row type. The separation of Jordanian barley landraces according to head type (two and six rows) is in general agreement with the results of Pasam et al.  and Cuesta-Marcos et al. . The molecular analysis did not reflect the origin of the collection sites for Jordanian barley landraces (data not shown). This might be because these landraces are adapted to a wide range of agro-ecological conditions and the probability of seed exchange between farmers. Russell et al.  found that landraces collected from Jordan and southern Syria clustered together indicating high similarity at the genetic level. In the contrary, Varshney et al.  found that the genetic diversity of barley accessions used in their study was not completely related to the geographic distribution.
To further elaborate on the assumptions extracted from the phenotypic and genotypic data of the tested Jordanian accessions, GWAS was carried out in order to identify QTLs associated with yield and yield components under rainfed conditions. Selecting the best P-value cutoff that control for both type I and type II errors in GWAS analyses is always challenging. Using the naïve false discovery rate (FDR) methods can lead to false negative outcomes as those methods were designed for independent tests while SNPs are dependent due to genetic linkage . To avoid FDR stringency, previous researches used suggestive P-value cutoff [48,49,50]; FDR cutoff at higher values up to 0.25 [51, 52]; or considered clusters of markers with R2 values >0.95 as dependent and can be represented by one marker [53, 54]. In this research, we calculated the number of independent test  followed by Bonferroni correction to accurately control for false negative results.
The GWAS analysis detected eleven significant SNP-trait associations for yield and it components (Table 3), while another 77 suggestive associations were reported in Additional file 9. These numbers are considered relatively close to other GWAS barley studies for yield and yield components under stress conditions [55, 56]. Furthermore, several QTLs identified were previously described in other related studies [47, 55, 56]. For instance, the QTL associated with BY, StY and K/S at overall detected on chromosome 1H (position: 133 cM) coincided with a previously described QTL in Mora et al.  and Varshney et al. . In another case, Wehner et al.  previously detected a QTL for BY on chromosome 5H (167.7 cM) in a region close to the GY QTL in Madaba 2008/2009 while they also detected a QTL associated with BY on 5H (position: 139.1 cM) that is considerably near a GY QTL identified in Ramtha 2009/2010 (additional file 9). The QTL located at 7H (position: 1.6 cM) associated with K/S at Madaba was also reported to be associated with grain number in Ingvordsen et al. . Interestingly, QTLs identified on 2H (Position: 106–108 cM), and associated with several traits across different environments was also reported in Ingvordsen et al.  and von Korff et al. . The QTL located at 5H (position 44 cM) associated with GY in Ramtha and other yield component across all environments was reported in Wehner et al. . This region includes several stress-related genes including an Abscisic acid-inducible protein kinase, which is known to mediate drought responses in cereals . This could indicate the existence of conserved drought tolerance mechanisms since early domestication knowing that wild barley from Jordan is believed to be the primary ancestors in modern cultivated barley .
Nevertheless, QTLs identified in this study did not show clear consistency across the tested environments. This is in general agreement with the results of Varshney et al.  and Maccaferri et al.  and, who attributed this inconsistency to the existence of different genetic systems governing drought tolerance under dry field conditions. In addition, studying complex traits such as yield under dry environments using GWAS and highly structured germplasm seems to be less informative in some instances when compared with bi-parental populations .
In conclusion, we observed a high level of variation for all of the agro-morphological traits measured in a collection of Jordanian barley landraces evaluated for three years in drought prone locations. Such variation was further confirmed at the molecular level using SNP markers analysis, which demonstrates the existence of considerable genetic variability among these accessions. Marker trait association identified highly significant and suggestive major QTLs for yield and yield components under dry environments located at 1H (position: 133 cM), 2H (position: 106 cM) and 5H (position: 44.2 cM). Finally, the superior accessions identified in this study and their respective QTLs can be used in barley breeding programs and marker assisted approach to improve productivity and stability under dry field conditions in specific locations in marginal areas of Jordan and across arid region in the world.
Best linear unbiased estimates
Coefficient of variation
- G × E:
Genotype × Environment
Genome wide association study
Kernels per spike
Quantitative triat loci
Squared allele frequency correlation
Standardized coefficient of cultivar superiority
Thousand kernel weight
Grando S, von Bothmer R, Ceccarelli S. Genetic diversity of barley: use of locally adapted germplasm to enhance yield and yield stability of barley in dry areas. In: Cooper HD, Spillane C, Hodgink T, (eds) Broadening the Genetic Base of crop production. CABI, New York, USA. 2001; 351-372.
Berthaud J, Clément J, Emperaire L, Louette D, Pinton F, Sanou J, Second S. The role of local level gene flow in enhancing and maintaining genetic diversity. In: Cooper HD, Spillane C, Hodgkin T, editors. Broadening the Genetic Base of Crop Production. New York: CABI Publishing in association with FAO and IPGRI; 2001. pp. 81–103.
Brush S. In situ conservation of landraces in centers of crop diversity. Crop Sci. 1995;35:346–54.
Ceccarelli S. Utilization of landraces and H. spontaneum in barley breeding for dry areas. Rachis. 1984;3(2):8–1l.
Grando S, Baum M, Ceccarelli S, Goodchild A, Jaby El-Haramein F, Jahoor A, Backes G. QTLs for straw quality characteristics identified in recombinant inbred lines of a Hordeum vulgare x H. spontaneum cross in a Mediterranean environment. Theor Appl Genet. 2005;110:688–95.
Collard BC, Mackill DJ. Marker-assisted selection: an approach for precision plant breeding in the twenty-first century. Philos Trans R Soc Lond Ser B Biol Sci. 2008;363(1491):557–72.
Ceccarelli S. Efficiency of plant breeding. Crop Sci. 2015;55:87–97.
Bergelson J, Roux F. Towards identifying genes underlying ecologically relevant traits in Arabidopsis Thaliana. Nat Rev Genet. 2010;11(12):867–79.
Comadran J, Kilian B, Russell J, Ramsay L, Stein N, Ganal M, Shaw P, Bayer M, Thomas W, Marshall D, Hedley P, Tondelli A, Pecchioni N, Francia E, Korzun V, Walther A, Waugh R. Natural variation in a homolog of Antirrhinum CENTRORADIALIS contributed to spring growth habit and environmental adaptation in cultivated barley. Nat Genet. 2012;44(12):1388–92.
Pasam R, Sharma R, Malosetti M, van Eeuwijk F, Haseneyer G, Kilian B, Graner A. Genome-wide association studies for agronomical traits in a worldwide spring barley collection. BMC Plant Biol. 2012;12(16):1–22.
Wang M, Jiang N, Jia T, Leach L, Cockram J, Waugh R, Ramsay L, Thomas B, Luo Z. Genome-wide association mapping of agronomic and morphologic traits in highly structured populations of barley cultivars. Theor Appl Genet. 2012;124(2):233–46.
Inostroza L, Del Pozo A, Matus I, Castillo D, Hayes P, Machado S, Corey A. Association mapping of plant height, yield, and yield stability in recombinant chromosome substitution lines (RCSLs) using Hordeum vulgare subsp. Spontaneum as a source of donor alleles in a Hordeum vulgare subsp. vulgare back-ground. Mol Breed. 2009;23:365–76.
Mickelson S, See D, Meyer FD, Garner JP, Foster CR, Blake TK, Fischer AM. Mapping of QTL associated with nitrogen storage and remobilization in barley (Hordeum vulgare L.) leaves. J Exp Bot. 2003;54(383):801–12.
Visioni A, Tondelli A, Francia E, Pswarayi A, Malosetti M, Russell J, Thomas W, Waugh R, Pecchioni N, Romagosa I, Comadran J. Genome-wide association mapping of frost tolerance in barley (Hordeum vulgare L.). BMC Genomics. 2013;14:424–37.
Munoz-Amatriain M, Cuesta-Marco A, Endelman JB, Comadran J, Bonman JM, Bockelman HE, Chao S, Russell J, Waugh R, Hayes PM, Muehlbauer G. The USDA barley Core collection: genetic diversity, population structure, and potential for genome-wide association studies. PLoS One. 2014;9(4):e94688.
Singh M, Malhotra R, Ceccarelli S, Sarker S, Grando S, Erskine W. Spatial variability models to improve dryland field trials. Exp Agric. 2003;39:1–10.
Rollins JA, Drosse B, Mulki MA, Grando S,·Baum M, Singh M, Ceccarelli S,· von Korff M. Variation at the vernalisation genes Vrn-H1 and Vrn-H2 determines growth and yield stability in barley (Hordeum vulgare) grown under dryland conditions in Syria. Theor Appl Genet 2013; 26: 2803–2282.
Francis TR, Kannenberg LW. Yield stability studies in short-season maize. 1. A descriptive method for grouping genotypes. Can J Plant Sci. 1978;58:1029–34.
Lin CS, Binns MR. A superiority measure of cultivar performance for cultivar x location data. Can J Plant Sci. 1988;68:193–8.
Payne R, Harding S, Murray D, Soutar D, Baird D, Welham S, Kane A, Gilmour A, Thompson R, Webster R, Tunnicliffe WG. GenStat for windows. 15th ed. Hemel Hempstead: Introduction. VSN International; 2011.
Yan W. GGEbiplot – a windows application for graphical analysis of multi-environment trial data and other types of two way data. Agron J. 2001;93:1111–8.
Costa JM, Bollero GA. Stability analysis of grain yield in barley (Hordeum Vulgare) in the US mid-Atlantic region. Ann Appl Biol. 2001;139:137.
Okuyama LA, Federizzi LC, Barbosa JF. Grain yield stability of wheat genotypes under irrigated and non-irrigated conditions. Braz Arch Biol Technol. 2005;48(5):697.
Winnepenninckx B, Backeljau T, De Wachter R. Extraction of high molecular weight DNA from mollusks. Trends Genet. 1993;9:407.
Perrier X, Jacquemoud-Collet JP: DARwin software. 2006, [http://darwin.cirad.fr/].
Pritchard J, Stephens M, Donnelly P. Inference of population structure using multilocus genotype data. Genet. 2000;155:945–59.
Evanno G, Regnaut S, Goudet J. Detecting the number of clusters of individuals using the software structure: a simulation study. Mol Ecol. 2005;14:2611–20.
Hedrick PW. Gametic disequilibrium measures: proceed with caution. Genet. 1987;117:331–74.
Weir BS. Genetic data analysis II: methods for discrete population genetic data. Sinauer Associates, Inc. Sunderland, Massachusetts. 1996.
Andreescu C, Avendano S, Brown SR, Hassen A, Lamont SJ, Dekkers JC. Linkage disequilibrium in related breeding lines of chickens. Genet. 2007;177:2161–9.
Massman J, Cooper B, Horsley R, Neate S, Macky RD, Chao S, Dong Y, Schwarz P, Muehlbauer GJ, Smith KP. Genome-wide association mapping of Fusarium head blight resistance in contemporary barley breeding germplasm. Mol Breed. 2011;27:439–54.
Bradbury JC, Zhang Z, Kroon DE, Casstevens TM, Ramdoss Y, Buckler ES. TASSEL: software for association mapping of complex traits in diverse samples. Bioinform. 2007;23:2633–5.
Thornsberry JM, Goodman M, Doebley J, Kresovich S, Nielsen D, Edward ES, Buckler IV. Dwarf8 polymorphisms associate with variation in flowering time. Nat Genet. 2001;28:286–9.
Yu J, Buckler ES. Genetic association mapping and genome organization of maize. Curr Opin Biotechnol. 2006;17:155–60.
Li MX, Yeung JM, Cherny SS, Sham PC. Evaluating the effective numbers of independent tests and significant p-value thresholds in commercial genotyping arrays and public imputation reference datasets. Human Genet. 2012;131(5):747–56.
Harlan J, Zohary D. Distribution of wild wheat and barley. Science. 1966;153:1074–80.
Salamini F, Özkan H, Brandolini A, Schäfer-Pregl R, Martin W. Genetics and geography of wild cereal domestication in the near east. Nat Rev Genet. 2002;3:429–41.
Ceccarelli S, Grando S, van Leur J. Genetic diversity in barley landraces from Syria and Jordan. Euphytica. 1987;36:389–405.
Rodriguez M, Rau D, Papa R, Attene G. Genotype by environment interactions in barley (Hordeum vulgare L.): different responses of landraces, recombinant inbred lines and varieties to Mediterranean environment. Euphytica. 2007;163(2):231–47.
Rizza F, Badeck F, Cattivelli L, Lidestri O, Fonzo N, Stanca A. Use of a water stress index to identify barley genotypes adapted to Rainfed and irrigated conditions. Crop Sci. 2004;44:2127–37.
Sudaric A, Simic D, Vrataric M. Characterization of genotype by environment interactions in soybean breeding programmes of southeast Europe. Plant Breed. 2006;125:191–4.
Mohammadi R. Amri .Genotype x environment interaction for durum wheat grain yield and selection for drought tolerance in irrigated and droughted environments in Iran. J Crop Sci Biotechnol. 2011;14(4):265–74.
Mohammadi R, Mohammadi M, Karimizadeh R, Amri A. Analysis of genotype-by-environment interaction for grain yield of Rainfed durum wheat genotypes in warm winter areas of Iran. J Crop Sci Biotechnol. 2010;13(4):267–74.
Close T, Prasanna R, Lonardi L, Wu Y, Rostoks N, Ramsay L, Druka A, Stein N, Svensson J, Wanamaker S, Bozdag S, Roose M, Moscou M, Chao S, Varshney R, Szűcs P, Sato K, Hayes P, Matthews D, Kleinhofs A, Muehlbauer G, DeYoung J, Marshall D, Madishetty K, Fenton R, Condamine P, Graner A. Waugh. Development and implementation of high-throughput SNP genotyping in barley. BMC Genomics. 2009;10:582.
Cuesta-Marcos A, Szucs P, Close TJ, Filichkin T, Muehlbauer GJ, Smith KP, Hayes PM, Genome-wide SNP. Re-sequencing of growth habit and inflorescence genes in barley: implications for association mapping in germplasm arrays varying in size and structure. BMC Genomics. 2010;15:11–707.
Russell J, Dawson I, Flavell A, Steffenson B, Weltzien E, Booth A, Ceccarelli S, Grando S, Waugh R. Analysis of >1000 single nucleotide polymorphisms in geographically matched samples of landrace and wild barley indicates secondary contact and chromosome level differences in diversity around domestication genes. New Phytol. 2011;191:564–78.
Varshney RK, Paulo MJ, Grando S, Van Eeuwijk FA, Keizer LCP, Guo PP, Ceccarelli S, Kilian A, Baum M, Graner A. Genome wide association analyses for drought tolerance related traits in barley (Hordeum vulgare L.). Field Crops Res. 2012;126:171–80.
Jighly A, Alagu M, Makdis F, Singh M, Singh S, Emebiri LC, Ogbonnaya FC. Genomic regions conferring resistance to multiple fungal pathogens in synthetic hexaploid wheat. Mol Breed. 2016;36:127.
Tadesse W, Ogbonnaya FC, Jighly A, Nazari K, Rajaram S, Baum M. Association mapping of resistance to yellow rust in winter wheat cultivars and elite genotypes. Crop Sci. 2014;54(2):607–16.
Tadesse W, Ogbonnaya FC, Jighly A, Sanchez-Garcia M, Sohail Q, Rajaram S, Baum M. Genome-wide association mapping of yield and grain quality traits in winter wheat genotypes. PLoS One. 2015;10(10):e0141339.
Newell MA, Asoro FG, Scott MP, White PJ, Beavis WD, Jannink J. Genome-wide association study for oat (Avena sativa L.) beta-glucan concentration using germplasm of worldwide origin. Theor Appl Genet. 2012;125:1687–96.
Pozniak CJ, Clarke JM, Clarke FR. Potential for detection of marker–trait associations in durum wheat using unbalanced, historical phenotypic datasets. Mol Breed. 2012;30:1537–50.
Manickavelu A, Joukhadar R, Jighly A, Lan C, Huerta-Espino J, Stanikzai AS, Kilian A, Singh RP, Ban T. Genome wide association mapping of stripe rust resistance in afghan wheat landraces. Plant Sci. 2016;252:222–9.
Jighly A, Oyiga BC, Makdis F, Nazari K, Youssef O, Tadesse W, Abdalla O, Ogbonnaya FC. Genome-wide DArT and SNP scan for QTL associated with resistance to stripe rust (Puccinia striiformis f. Sp. tritici) in elite ICARDA wheat (Triticum aestivum L.) germplasm. Theor Appl Genet. 2015;128(7):1277–95.
Ingvordsen CH, Backes G, Lyngkjær MF, Peltonen-Sainio P, Jahoor A, Mikkelsen TN, Jørgensen RB. Genome-wide association study of production and stability traits in barley cultivated under future climate scenarios. Mol Breed. 2015;35(3):1–14.
Wehner GG, Balko CC, Enders MM, Humbeck KK, Ordon FF. Identification of genomic regions involved in tolerance to drought stress and drought stress induced leaf senescence in juvenile barley. BMC Plant Biol. 2015;15:125.
Mora F, Quitral YA, Matus I, Russell J, Waugh R, del Pozo A. SNP-based QTL mapping of 15 complex traits in barley under rain-fed and well-watered conditions by a mixed modeling approach. Front Plant Sci. 2016;7:909. https://doi.org/10.3389/fpls.2016.00909.
Anderberg RJ, Walker-Simmons M. Isolation of a wheat cDNA clone for an abscisic acid-inducible transcript with homology to protein kinases. Proc Natl Acad Sci U S A. 1992;89(21):10183–7.
von Korff M, Grando S, This D, Baum M, Ceccarelli S. Quantitative trait loci (QTL) associated with agronomic performance of barley under drought. Theor Appl Genet. 2008;117:653–69.
Pankin A, von Korff M. Co-evolution of methods and thoughts in cereal domestication studies: a tale of barley (Hordeum Vulgare). Curr Opin Plant Bio. 2017;36:15–21.
Maccaferri M, Sanguineti MC, Demontis A, El-Ahmed A, Garcia del Moral L, Maalouf F, Nachit M, Nserallah N, Ouabbou H, Rhouma S, Royo C, Villegas D, Tuberosa R. Association mapping in durum wheat grown across a broad range of water regimes. J Exp Bot. 2011;62(2):409–38.
We gratefully acknowledge Miss Shireen Qasrawi for her technical assistance. We sincerely acknowledge Jordanian farmers for their permission to use their fields for experimental work. We gratefully acknowledge NCARE for their assistance in field work. We also acknowledge the statistical support received from the Biometrics and Statistics Section of ICARDA and special thanks to Dr. Murari Singh for his great assistance.
This work was supported in part by a grant from International Development Research Center (IDRC)/Canada and in part by a grant from the Deanship of Scientific Research, The University of Jordan.
Availability of data and materials
The data sets supporting the results of this article will be freely available upon request to corresponding author: firstname.lastname@example.org for non-commercial use only. The genotypic data of the current study can be downloaded from Additional file 11.
Ethics approval and consent to participate
Permission from the farmers to carry out field experiments were obtained under the direct supervision from NCARE/Jordan.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
ICARDA Gene bank information of 154 Jordanian barley genotypes used in this study. (XLSX 21 kb)
Characteristics of field locations, actual total rainfall (mm on monthly basis), number of rainy days and long-term rainfall average in the tested sites over three growing seasons. (XLSX 10 kb)
Combined analysis of variance for grain yield of barley landraces accessions and Restricted Maximum Likelihood Method estimated variances for locations and overall environments. (XLSX 14 kb)
Genotypic coefficient of variation, heritability, coefficient of variation and P-value for grain yield of barley accessions evaluated in eight environments. (XLSX 9 kb)
Means, coefficient of variations (CV%) and standardized cultivar superiority over applicable years of the best linear unbiased estimates (BLUEs) of GY (g/m2) of all tested genotypes. (XLSX 32 kb)
Correlation heat maps between selected triats of all tested genotypes for each indvidual field trial. (XLSX 166 kb)
Population structure of 154 Jordanian barley genotypes used in this study. (XLSX 34 kb)
Decline of Genome-wide LD (as measured by R2) against the genetic distance (in cM). (XLSX 109 kb)
Suggestive markers associated with yield and yield components identified in this study. (XLSX 20 kb)
QQ plots for all tested triats in this study. (PDF 360 kb)
Genotypic data of Jordanian barley landraces and four local checks used in this study. (XLSX 2582 kb)