Genome Wide Association Mapping of Seedling and Adult Plant Resistance to Barley Stripe Rust (Puccinia striiformis f. sp. hordei) in India

Barley stripe rust is caused by Puccinia striiformis f.sp. hordei, (Psh), occurs worldwide, and is a major disease in South Asia. The aim of this work was to identify and estimate effects of loci underlying quantitative resistance to rust at seedling and adult plant stages. HI-AM panel of 261 barley genotypes consisting of released cultivars from North and South America, Europe, Australia, advanced breeding lines, and local landraces from ICARDA barley program were screened at seedling and adult plant stages for resistance to Psh. Seedling resistance was evaluated with the five prevalent Psh races in India. Screening for the adult plant stage resistance was also performed in two different locations by inoculating with a mixture of the five races used for seedling screeing. The panel was genotyped using DaRT-Seq high-throughput genotyping platform. The genome-wide association mapping (GWAM) showed a total of 45 QTL located across the seven barley chromosomes for seedling resistance to the five races and 18 QTL for adult plant stage resistance. Common QTL for different races at seedling stage were found on all chromosomes except on chromosome 1H. Four common QTL associated with seedling and adult plant stage resistance were found on chromosomes 2, 5, and 6H. Moreover, one of the QTL located on the long arm of chromosome 5H showed stable effects across environments for adult plant stage resistance. Several QTL identified in this study were also reported before in bi-parental and association mapping populations studies validating current GWAM. However 15 new QTL were found at adult plant stage on all chromosomes except the 4H, explaining up to 36.79% of the variance. The promising QTL detected at both stages, once validated, can be used for MAS in Psh resistance breeding program globally.

Barley stripe rust is caused by Puccinia striiformis f.sp. hordei, (Psh), occurs worldwide, and is a major disease in South Asia. The aim of this work was to identify and estimate effects of loci underlying quantitative resistance to rust at seedling and adult plant stages. HI-AM panel of 261 barley genotypes consisting of released cultivars from North and South America, Europe, Australia, advanced breeding lines, and local landraces from ICARDA barley program were screened at seedling and adult plant stages for resistance to Psh. Seedling resistance was evaluated with the five prevalent Psh races in India. Screening for the adult plant stage resistance was also performed in two different locations by inoculating with a mixture of the five races used for seedling screeing. The panel was genotyped using DaRT-Seq high-throughput genotyping platform. The genome-wide association mapping (GWAM) showed a total of 45 QTL located across the seven barley chromosomes for seedling resistance to the five races and 18 QTL for adult plant stage resistance. Common QTL for different races at seedling stage were found on all chromosomes except on chromosome 1H. Four common QTL associated with seedling and adult plant stage resistance were found on chromosomes 2, 5, and 6H. Moreover, one of the QTL located on the long arm of chromosome 5H showed stable effects across environments for adult plant stage resistance. Several QTL identified in this study were also reported before in bi-parental and association mapping populations studies validating current GWAM. However 15 new QTL were found at adult plant stage on all chromosomes except the 4H, explaining up to 36.79% of the variance. The promising QTL detected at both stages, once validated, can be used for MAS in Psh resistance breeding program globally.
Keywords: Puccinia striiformis f. sp. hordei, GWAM, seedling stage, barley (Hordeum vulgare L.), adult plant resistance (APR), QTL INTRODUCTION Barley (Hordeum vulgare L.) is one of the most important cereal crops in the word and especially in the dry areas where often is the only crop that can be grown under extreme drought conditions (Ceccarelli, 1994;Li et al., 2006). On the other hand, barley managed by irrigation and high rainfall is common in South Asia and East Africa as well as other regions where rusts and foliar blights are important production constraints. The stripe rust caused by Puccinia striiformis f. sp. hordei, (Psh) is the major constraint in South Asia, East Africa, and Central and North America affecting both quantity and quality of barley produced (Luthra and Chopra, 1990;Roelfs and Huerta-Espino, 1994). Though, Psh can be effectively managed with fungicides, but the use of resistant varieties is considered the most sustainable option for both environmental and economic reasons. Incorporation of qualitative and quantitative resistance is important to obtain cultivars with durable resistance to stripe rust (St. Clair, 2010). Quantitative rust resistance is mediated by quantitative trait loci (QTL) conferring partial/non race specific resistance or "slow rusting" type of resistance. Many QTL reduce the disease severity by increasing rust incubation or latent periods (Mundt, 2014). Usually slow rusting genes have small or intermediate effects when present alone but a higher degree of durable resistance can be achieved by combining 4 or 5 such genes (Singh et al., 2000;Herrera-Foessel, 2011). Quantitative resistance is usually considered race-non-specific but as reported by Poland et al. (2009), the biological mechanisms underlying it is poorly understood and a wide variation of mechanisms is expected. Rust pathogens have intrinsic characteristic like windaided migration, ability to easily increase population size and to mutate and acquire new virulence to resistance genes. Due to these characteristics quantitative resistance is preferable to qualitative resistance for long lasting cultivation of new resistant cultivars. Furthermore, with more than 70 Psh races identified (Hovmøller, 2002) a better understanding of the genetic control of quantitative resistance is of crucial importance and breeding would be more effective if based on extensive knowledge of the resistance genes/QTL (Gutiérrez et al., 2015). A previuous study performed by Verma et al. (2018), using the same germplasm, identified accessions carrying already know resistance genes and Psh resistant genotypes to all the five races that may possess novel resistance genes. The current study was undertaken to identify QTL effective against Psh individual races at seedling stage and QTL for quantitative resistance in field globally and especially for South Asia where annual recurrence of stripe rust from the Himalayas is a big challenge.

Plant Material and Yellow Rust Races
Two hundred sixty-one spring barley genotypes (172 two-row and 89 six-row types), including released cultivars, advanced breeding lines, and landraces were used in this study. The set is named as HI-AM (High Input Association Mapping) panel as mostly barley genotypes were from ICARDA breeding program targeted toward optimum management conditions. Out of the 261 genotypes, 124 were from ICARDA's barley breeding program (50 two-row and 74 six-row type), 32 from Europe (28 two-row and 4 six-row type), 34 North America (28 two-row and 6 six-row type), 67 from South America (62 two-row and 5 six-row type) and 4 from Australia (only two-row type). The seedling resistance test (SRT) for HI-AM was done against five prevalent Psh races in India, namely, Q (5S0), 24 (0S0-1), 57 (0S0), M (1S0), and G (4S0) under glasshouse at Regional Station, Indian Institute of Wheat and Barley Research (IIWBR), Shimla, India. Among these races, 57, 24, and G were old races while Q and M were identified and characterized more recently (Nayar et al., 1997;Prashar et al., 2007). The barley geneotypes were also screened at their Adult-plant stage (APS) at two locations using mixture of the five races received from Regional Station, IIWBR, Shimla, India.
Screening for Seedling Resistance to Five Races of Psh Seeds (5 seeds per hole per genotype) were sown in aluminum trays filled with a mixture of fine loam soil and farmyard manure. A total of 18 genotypes and a susceptible check "Bilara-2" were planted in each tray. Seedlings were inoculated 1 week when primary leaves were fully expanded. Inoculation was done using a glass atomizer with 100 mg of rust spores of each race suspended in 10 ml light grade mineral oil, Soltrol 170, (Chevron Phillips Chemicals Asia Pvt. Ltd., Singapore). Inoculated seedlings were misted with sterile water and placed for 48 h in dew chambers at 16 ± 2 • C with >90% relative humidity and 12 h day/night cycle. Seedlings were then transferred to glasshouse and incubated at 16 ± 2 • C with >70% relative humidity, illuminated at about 15,000 lux for 12 h. Leaves were treated with fine elemental sulfur to prevent powdery mildew infection, without affecting rust infections. Infection types (ITs), were recorded 16-18 days after inoculation following the modified methods of Nayar et al. (1997) and Stakman et al. (1962). Infection types 0 to 2 were considered resistant and 3 to 33 + as susceptible (Table 1A). Genotypes with resiant recation were tested again against the same races to confirm their reaction. The experiment was repeated once with resistant genotypes only to ascertain the consistency of the ITs.

Screening for Adult-Plant Stage Resistance
All genotypes used in the seedling stage resistance screening were also evalauted for APS under artificial inoculation conditions at RARI, Durgapura (75 • 47 ′ E, 26 • 51 ′ N), Rajasthan, in crop seasons 2012-13 and 2013-14 and at ICAR-IIWBR, Karnal (76 • 98 ′ E, 29 • 69 ′ N),Haryana during 2013-2014, in India. The experiment was laid out in augmented design with yellow rust susceptible check, "Bilara-2, " after every 20 test genotypes in both locations and all seasons. Seeds were sown in 1-m rows with 25 cm row to row distance for each genotype in the first fortnight of November each year at Durgapura and Karnal. The infector rows of Bilara-2 were also sown perpendicular to the test material rows all around the blocks. In order to multiply the initial inoculum six rows of Bilara-2 (about 50 M long) were sown close in the same field about 15 days before the sowing of experimental materials. The five races Q (5S0), 24 (0S0-1), 57 (0S0), M (1S0), and G (4S0) received from Regional Station, IIWBR, Shimla were multiplied on this large plot of Bilara 2. The infector rows were syringe inoculated at seedling stages (Zadoks GS 20) with the mixed inoculum of stripe rust races taken from early sown Bilara 2 plot, followed by repeated sprays of inoculum collected from infector rows as well as from early sown Bilara 2 plot on the test material. The field was given extra irrigations (sprinkler) to maintain an appropriate humid microclimate for better disease devlopment. Stripe rust severity was recorded at the early to late flowering stages (Zadoks GS 60-69) when maximum disease severity reached on the susceptible check rows, sown after every 20 rows of test material.
The modified Cobb's scale (Peterson et al., 1948) was used for classifying geneotypes to different reaction groups (Table 1B). Genotypes were classified as moderately resistant (up to 10 MR); moderately susceptible (5S to 20 MS); susceptible (>20MS-60S) and highly susceptible (>60S) based on stripe rust severity over two seasons and two locations.

Genotyping With DaRT-Seq
Total genomic DNA extraction was performed on a single plant from each genotype from the HI-AM panel using a Biosprint 96 DNA Plant Kit (Qiagen, Hilden, Germany). DNA sample were processed for DaRT-Seq (Diversity Array Technology Pty Ltd, DaRT P/L) by a series of digestion/ligation reactions (Cruz et al., 2013;Ames et al., 2015). The system combines complexity reduction methods with next-generation sequencing platforms, targeting primarily genic regions (Carling et al., 2015). It produces two types of markers, SilicoDaRT markers characterized by presence/absence variation (PAVs) and classical SNPs present in the sequenced fragments (http://www. diversityarrays.com/dartapplication-dartseq-data-types). The PAV/SNP markers were subsequently aligned by using sequence information available at ftp://ftpmips.helmholtz-muenchen. de/plants/barley/public_data/anchoring. The Thresholds for minimum base id of 90% and E-value of 5 −10 were imposed to declare positive matches against the available datasets of the physical map (ftp://ftpmips.helmholtz-muenchen.de/plants/ barley/public_data/). Markers quality control of the initial dataset was conducted by removing heterozygous and monomorphic markers and markers with minor allele frequencies (MAF) < 5% and markers with missing data > 10%.

Population Structure and Linkage Disequilibrium
The genetic structure of the 261 genotypes of the HI-AM panel was investigated using 105 and 101 unlinked markers from the PAV set and from the SNP set, respectively, distributed across the 7 barley chromosomes. Markers subsets used for population structure assessment were obtained by selecting one marker every 10 cM within both PAVs and SNPs sets, in order to avoid to enclose linked markers in the subsets. Population structure was firstly determined using STRUCTURE version 2.3.4 (Pritchard et al., 2000). The admixture model option was run using a burnin length of 10 10 cycles, in order to minimize the effect of starting configuration, and a simulation of 10 6 cycles was applied. Cluster values (k) from 2 to 10 were chosen and 5 independent runs for each k were chosen to obtain consistent results. Additionally, the adegenet package for R statistical software (The R Development core team) was used to confirm the number of sub-populations by the Bayesian Information Criterion (BIC). Finally, on the base of PCA results, genotypes were assigned to subgroups or considered admixed on the base of 80% membership criterion. A principal component analysis was also used to determine population structure and used as covariate in the subsequent GWAM study. Linkage disequilibrium (LD) was estimated with Tassel software V 5.2.32 (Bradbury et al., 2007) using a subset of 1,577 polymorphic markers with known position selected from the original SNP marker-set. Linkage disequilibrium was calculated, using the SNPs marker set, separately for locus pairs within the same chromosomes and between chromosomes. LD was estimated as the squared allele frequency correlations (R 2 ) with only p-values ≤ 0.01 for each pair of loci considered as significant. The nlstools package for R Statistical Software (The R Development core team) was used to estimate the extent of LD by non-linear regression analysis on the basis of intrachromosomial r 2 values (Hill and Weir, 1988;Remington et al., 2001).

Genome Wide Association Mapping
Disease severity scores at seedling and adult plant stages and the genotypic data were used to perform GWAM using Tassel V 5.2.32. GWAM was performed using both General Linear Model (GLM) and Mixed Linear Model (MLM) methods. The general equations for GLM and MLM are the followings: y = Xa + e and y = Xa + Qb + Ku + e, respectively. X denotes at the marker while Q is the Q-matrix obtained by STRUCTURE software and K is the kinship matrix (Q+K). The vector for phenotypes is indicated as y while a is the effect of marker fixed effects and b represent the vector of fixed effects, while u is an unknown vector of random additive genetic effects. Association analysis using the GLM model was performed incorporating as covariate population structure derived from PCA analysis (GLM+PCA model) or the Q-matrix (GLM+Q model) in order to avoid type I errors. The MLM model consider the familiar relatedness (the K model) and takes into consideration both population structure and familiar relatedness (Q+K and PCA+K models). The kinship matrix (K) was estimated using Tassel V 5.2.32 from the whole set of DaRT markers with unique position. For comparison we also conducted GWAM without any correction for population substructure. For all scans threshold of (−log 10 p ≥ 3) was set for identifying significant marker-trait associations. Significant markers mapping within 5 cM of each other were considered as being linked to the same QTL and the marker with the highest p-value was chosen as representing the QTL. We firstly determined the critical p-value for the significance of marker-trait association using the false discovery rate (FDR). Since the FDR was found to be highly stringent and considering the stringency of the model used for accounting for population structure, in which most of the false positives were inherently controlled. Markers were declared significant at the p = 0.0001 [−log(p) = 3] with the selected models according to a liberal approach (Chan et al., 2010).

QTL Alignment and Candidate Genes
To align detected QTL with those previously reported in different barley germplasm, we checked for the presence of common markers and/or positions in the barleymap database (http:// floresta.eead.csic.es/barleymap/; Cantalapiedra et al., 2015), that allows to search the position of barley genetic markers on the Barley Physical Map (IBSC, 2012). The position of the marker representative of the QTL was compared with those of markers at QTL peaks reported in previous studies. Candidate genes search was done using the PGSB database 1 that provides access to the barley gene annotation described by the IBSC (2012). The markers at the QTL peaks, were used to search, in the genomic region encompassing the QTL, for functional domains or genes functionally related with disease resistance mechanisms.

Seedling and Adult Plant Resistance
The statistical analysis of the phenotypic data collected has been already reported by Verma et al. (2018), while the reaction type of barley genotypes at SRT and APS are summarized in Table 1. The results showed that the genotypes of HI-AM panel were resistant to races 24 (54.8%) and M (62 %), followed by race G (47.5%), while for the races 57 (31.4%), and Q (16.5%). Races Q and 57 were the most virulent races where over 84 and 68% of genotypes were susceptible (Table 1A).
In the APS screening, 46.7 and 27.9% genotypes were resistant at Durgapura during 2013 and 2014, respectively and 77% resistant at Karnal during 2014 (Table 1B). Durgapura represents an optimal site for Psh due to climatic conditions (temperature and humidity) that favors Psh development while Karnal is characterized by the occurrence of severe winters that may limit the development and sporulation of the rust for the secondary spread. However, Karnal also represents the location in Haryana, which is also prone to stripe rust losses in the wider region. We observed a variation in the range of field reactions between Dg13 and Dg14, as well as between Durgapura and Karnal reactions. As expected data collected in Durgapura indicated a wide range of reaction types when compared with data collected in Karnal. In both locations the susceptible check showed the maximum level of susceptibility (100S) indicating that the variations observed in reactions of the test genotypes are due to climatic factors and not because of escape or inoculum load.

Marker Statistics, Population Structure, and Linkage Disequilibrium
Markers selected after filtering steps have been subsequently used to estimate population structure, linkage disequilibrium (LD) and to perform GWAM for both SRT and APS. The final sets of markers comprise 13.182 PAV and 6.311 SNPs. Population structure analysis performed with Structure (Pritchard et al., 2000) showed three subpopulations (k = 3) on the base of the k parameter and according to Evanno et al. (2005), the same number of k was also confirmed by the BIC estimation perfomed with R Statistical Software (The R Development core team). The population structure is shown in Figures 1A,B. The first subpopulation (Q1) is mainly composed by ICARDA germplasm (70%). The second group (Q2), located alone in the left side of the PCA chart ( Figure 1A) shows the higher degree of diversity and the highest number of entries from South America (37%), followed by ICARDA germplasm (22%), North America (21%), Europe (13%), and Australia (3%). Q3 was again mainly composed by ICARDA germplasm (70%). The 6-row genotypes are spread across the three subgroups representing the 53% (Q1), 26% (Q2), and 45% (Q3) of the total number of genotypes for ech subpopulation, respectively. There were 153872 (24.4%) inter-chromosomal pairs of loci showing significant LD (P < 0.01), 3847 (2.5%) of which had R 2 > 0.2. Of the intra-chromosomal locus pairs, 34227 (31.5%) had a significant LD of which 3022 (8.89%) had R 2 > 0.2 Intra-chromosomal locus pairs have a higher mean R 2 value (0.10) than inter-chromosomal locus pairs (0.02). The scatter plots of LD (R 2 ) as a function of the inter-marker distance (cM) within the same chromosome for all genotypes indicated a clear LD decay at 4 cM (R 2 = 0.18) with genetic distance (Figure 2).

Genome Wide Association Mapping
Using the PAVs marker set, the best model fitting was the MLM using Q+K model when analyzing SRT data; while PCA+K was the more suitable for APS data analyses. Using the SNP dataset the best fitting model for SRT and APS was the GLM model using PCA for accounting population structure and relatedness. The QQ plots for GWAM are available as Supplementary Materials (Figures S1-S4). The GWAM analyses at SRT identified 45 QTL located across the seven barley chromosomes ( Table 2). The marker R 2 ranged from 4.25% to 6.56% Table 2). The race specific QTL detected for SRT explained together the 41.77 % (Race Q), 50.1% (Race 24), 36.42% (Race 57), 53.0% (Race G), and 49.84% (Race M) of phenotypic variance, respectively. GWAM for APS showed 18 QTL using phenotypic data from two locations and during two seasons ( Table 3). The marker R 2 ranged from 4.54 to 8.11% for APS (Table 3). A QTL located on chromosome 5H was found consistantly stable across seasons and environments. Phenotypic variance explained by QTL detected in case of APS was 15.35% for Dg13, 36.79% for Dg14 and 45.82% for Kr14.    Among the QTL detected, 8 QTL were significant for two races ( Table 4A). The QTL on chromosome 5H (cM 137.08), overlaps with a QTL for resistance at SRT for the race M. Furthermore, other 2 QTL detected for APS are located on chromosome 2H at the same position of QTL detected for resistance at SRT for race G (Table 4B). Similarly another QTL for APS located on chromosiome 6H overlaps with OTL for SRT for race Q.

Candidate Genes for Resistance to Psh
Out of the 45 QTL identified for resistance at SRT, 15 were coincident with prior reports Those QTL were reported from different barley germplasm and different Psh races (Thomas et al., 1995;Toojinda et al., 2000;Vales et al., 2005;Rao et al., 2007;Verhoeven et al., 2011;Gutiérrez et al., 2015;Dracatos et al., 2016;Esvelt Klos et al., 2016;Belcher et al., 2018). For Psh resistance at the APS stage using different germplasm and races, only 3 out of the 18 QTL detected, shown in Table 5, are coincident with previous reports (Vales et al., 2005;Rao et al., 2007;Verhoeven et al., 2011;Gutiérrez et al., 2015;Dracatos et al., 2016;Belcher et al., 2018). Candidate genes (CG) were found for both QTL at the SRT and APS (Table 5). QTL detected at both SRT and APS are located in genomic regions enriched in genes or functional domains that according to their annotations can be considered indicative of common R gene products like: nucleotide binding site (NBS), leucine rich repeat (LLR) and disease resistance protein, genes involved in b-glucan biosynthesis. We identified 27 (CG) for Psh resistance at SRT and 10 for APS. As expected most of the CGs belongs to the NBS-LLR disease resistance protein family, others are classified as generic disease resistance proteins and also as genes involved in β-glucans biosynthesis. We also find NAC protein overlapping with the QTL SRT _RG_4 on chromosome 4H and Lr21 for the QTL SRT _R24_1 located on the chromosome 1H.

DISCUSSION
Barley stripe rust specialization and race structure are poorly defined as compared with wheat rusts and only few studies on the genetic control of the resistance are available in barley. The lack of barley genetic stocks for resistance represent one of the main limiting factors in the identification of genetic determinats of Psh resistance (Dracatos et al., 2016). In order to identify genomic regions controlling resistance to five prevalent races of Psh in India at both SRT and APS stages, the study was taken up on the barley genotypes from much diverse sources. The SRT under artificial inoculation for the five races of the HI-AM panel showed that the most virulent races were Q and 57, respectively. For APS resistance, we observed variation in rust severity in the two seasons in Durgapura (Dg13 vs. Dg14), however, the range of field reaction was higher at Durgapura as compared to Karnal (Kr14). Durgapura is an optimal location for rust screening due to relatively less severe winter favoring faster stripe rust development (Verma et al., 2018), while, Karnal station with severe winters sometimes experience delayed strip rust development with less secondary spread. However, the susceptible check Bilara-2 showed high severy with susceptible reaction at both locations indicating that there is no escape and it was possible to select resistant genotypes at both locations. Beside climate variation, others factors also contribute to the unmanageable variation in infection response in field trials. Inoculum composition, sequential infection, and plant phenology are others main factors that can modify plant response to Psh and others pathogens (Hickey et al., 2011;Gutiérrez et al., 2015). Therefore testing genotypes in different environments is important to confirm the resistance of barley genotypes to Psh. As mentioned before the HI-AM comprises 261 spring genotypes from different breeding programs in different continents and landraces from germplasm collections. Ear type and origin of germplasm are often the main determinants of population subgrouping but in our case seems that both traits do not have any effect on population structure. The absence of subgrouping based on ear type may be due to the fact that the genotypes from breeding programs represent outcome of 2× 6-row hybridization being frequently used for germplasm improvement. Furthermore out of 89 six-row genotypes present in the panel, 73 were bred and selected at ICARDA, where crosses between 2 and 6-row genotypes are routinary. In the case of germplasm origin, the absence of subgrouping is most probably due to the extensive germplasm exchange between ICARDA and breeders in North and South American. The distribution of genotypes across the three different subgroups seems to be due to the breeding history of individuals, and further analysis based on pedigrees might be helpful to understand relationship of genotypes belonging to different groups.
Normally seedling resistance is considered qualitative and based on gene for gene interaction between host and pathogen while quantitative resistance has been defined as a non-race specific resistance expressed only in adult plants (Milus and Line, 1986;Richardson et al., 2006). APR is generally best expressed at adult phase and usually involve additive and/or epistatic effects of multiple genes that confers a durable partial resistance (Hickey et al., 2011). Those QTL may represent an interesting source of quantitative resistance and, if validated they can be introgressesed in breeding materials, through MAS, to combine both qualitative and quantitative resistance. Qualitative resistance mechanisms have been extensively studied in terms of genomic location and specificity (Giese et al., 1993;Thomas et al., 1995;Graner and Tekauz, 1996) while mechanisms underlying quantitative resistance still to be clarified. As reviewed by Richardson et al. (2006) quantitative resistance may be controlled by uncharacterized classes of R genes or by alternative alleles at   Thomas et al., 1995;Qi et al., 1998;Hayes et al., 2000). We detected several race specific QTL at SRT stage that were already reported as QTL for resistance at both SRT and APS in prior reports. The fact that most of them were detected at APS and using different races support the hypothesis that both qualitative and quantitative resistance genes may be located at the same loci. Interestingly the phenotypic variance explained by the QTL detected for the two most virulent races at SRT stage was the lowest (41.77% for race Q and 36.42% for race 57, respectively) when comparated with less virulent races. We detected 18 QTL for APS using a mixture of the five races, out of which only 3 were reported before. APS_Dg14_1 detected at APS (Verhoeven et al., 2011), while APS_Kr14_1 and APS_Kr_14_4 were detected by Dracatos et al. (2016) at SRT. Furthermore, we found that eight of the QTL detected for SRT were significant for two different races, furthermore QTL SRT _R0_2/SRT _RG_2, SRT _R24_6/SRT _RM_3 and SRT _R24_9/SRT _R0_5 shows opposite effects within races. We also found that 13 out of 18 QTL detected at APS were not coincident with those detected at SRT, therefore they can be considered QTL for Adult Plant Resistance (APR). Among the QTL detected the most interesting are located on the long arm of chromosome 5H at cM 137.08 (APS_Dg13_3, APS_Dg14_5, and APS_Kr_7), near the telomeric region, it was significant across environments and represents the best candidate for validation detected at APS stage. Commons QTL for SRT and APS resistance are located on chromosomes 2HS, 2HL, 5H, and 6HL respectively. It is noteworthy to mention that the QTL for APS on chromosome 5H at 137.08 cM is overlapping with QTL SRT_RM_6. The QTL SRT_RM_6 is located 3.39 cM (133.69 cM) from the QTL for APS and is position is coincident with another QTL previously reported for APS by Vales et al. (2005) and by Gutiérrez et al. (2015) located at 129.44 cM in the IBSC 2012 2 barley genetic map.
Again this could be explained by presence of clusters of R genes where qualitative genes may regulate the response in SRT, while quantitative genes control the non-race specific response at APS. In case of the overlapping race specific QTL, the opposite effects detected between different races may be due to different race specific resistance genes at the same loci or to the same resistance QTL/gene that responds in a race specific way. Since, as reported by Park (2008) the use of both SRT and APR is desirable for increasing durability of rust resistance in cereals those QTL also represent a valuable source of resistance to Psh. Furthermore, is noteworthy that the most part of QTL detected at APS stage shows negative effect on rust pathogenesis and if validated, they can be used in MAS in stripe rust resistance breeding in barley, especially in India and South Asia.
For several of QTL detected at both SRT and APS we were able to identify CGs; preference was given to genes predicted to be members of the most common family of R genes and other genes predicted to be relevant to stripe rust pathogenesis. QTL genomic region encompass genes annotated for functions or for domains related to disease resistance thus these genes can be considered putative candidate genes for the corresponding QTL. Those genes are part of most represented resistance gene family that encodes cytoplasmatic proteins with nucleotide-binding sites and several leucine rich repeats (Halterman et al., 2001). Furthermore we also found several candidate genes annotated as part of the in β-glucans biosynthesis pathway. β-glucans are component of cell wall and glucanases have a role in early plant response to fungal pathogens. In fact glucanases are secreted by cell walls that are directed toward degradation of fungal walls (Thomma, 2003;Veronese et al., 2003;Vorwerk et al., 2004). The CG associated with race specific QTL SRT_RG_4 was a protein with a NAC domain; proteins encoded by NAC gene family constitute a large family of specific transcription factors, involved in both abiotic and biotic stress response (Al Abdallat et al., 2014;McGrann et al., 2015). For instance, the wheat transcription factor TaNAC4 expression is induced in leaves by the infection of stripe rust (Xia et al., 2010). Many other genes implicated in plant defense response process may be co-located in the same regions. For instance Esvelt Klos et al. (2016) reported several CGs related with cellular reactive oxygen species that are known to play an important role in plant defense mechanisms. Different hypothesis, based on experimental evidences, have been postulated on the mechanisms underlying the quantitative disease resistance and further studies are required to decipher how this mechanisms confers resistance to disease like Psh.

CONCLUSION
Expanding the catalog of mapped QTL for stripe rust resistance and its validation represent an important step toward the application of MAS for the introgression and pyramiding of resistance genes in new barley cultivars. In this work, novel QTL for Psh resistance at SRT and adult plant stages were identified which could be helpful in dissection the resistance mechanism to this pathogen. New QTL need to be validated for their diversity, effectiveness in different genetic background and with more races of Psh existing in other regions of the world to ensure their use for introgression in barley germplasm or for MAS globally.