Genome-Wide Mapping of Quantitative Trait Loci Conferring All-Stage and High-Temperature Adult-Plant Resistance to Stripe Rust in Spring Wheat Landrace PI 181410

Stripe rust, caused by Puccinia striiformis f. sp. tritici (Pst), is one of the most destructive diseases of wheat in the world. Genetic resistance is the best strategy for control of the disease. Spring wheat landrace PI 181410 has shown high level resistance to stripe rust. The present study characterized the landrace to have both race-specific all-stage resistance and nonrace-specific high-temperature adult-plant (HTAP) resistance. To map quantitative trait loci (QTL) for the resistance in PI 181410, it was crossed with Avocet S (AvS), from which a recombinant inbred line population was developed. The F5–F8 populations were consecutively phenotyped for stripe rust response in multiple field environments under natural Pst infection, and the F7 population was phenotyped in seedlings at low temperature and in adult-plant stage with selected Pst races in the greenhouse. The F7 population was genotyped using the 90K wheat SNP chip. Three QTL, QYrPI181410.wgp-4AS, QYrPI181410.wgp-4BL, and QYrPI181410.wgp-5BL.1, from PI 181410 for all-stage resistance, were mapped on chromosome arms 4AS, 4BL, and 5BL, respectively. Four QTL, QYrPI181410.wgp-1BL, QYrPI181410.wgp-4BL, QYrPI181410.wgp-5AS, and QYrPI181410.wgp-5BL.2, were identified from PI 181410 for HTAP resistance and mapped to 1BL, 4BL, 5AS, and 5BL, respectively. Two QTL with minor effects on stripe rust response were identified from AvS and mapped to 2BS and 2BL. Four of the QTL from PI 181410 and one from AvS were potentially new. As the 4BL QTL was most effective and likely a new gene for stripe rust resistance, three kompetitive allele specific PCR (KASP) markers were developed for incorporating this gene into new wheat cultivars.


Introduction
Wheat (Triticum aestivum L.) is one of the most important food crops worldwide, as it provides humans most of their daily calories and proteins. With the growing of the human population, wheat demand is expected to increase by 66% in 2040 [1]. However, wheat production has been continually threatened by abiotic and biotic stresses. Stripe rust (also known as yellow rust), caused by fungal pathogen Puccinia striiformis Westend. f. sp. tritici Erikss. (Pst), is especially a devastating disease

Characterization of Stripe Rust Resistance in PI 181410
In the field, PI 181410 has always shown resistance with infection type (IT) 1-2 and severity (SE) 5%-20% to stripe rust at adult-plant stage in both Mount Vernon and Pullman, Washington, USA for 13 years since 2006 ( Figure S1). In contrast, the susceptible check, AvS, has been always highly susceptible (IT 8, SE 90%-100%) in the same nurseries. In the greenhouse tests, the seedlings of PI 181410 were resistant (IT 3) to Pst race PSTv-4, moderately resistant (IT 5) to PSTv-14, and susceptible (IT 8) to races PSTv-37, PSTv-40, and PSTv-51 in the low-temperature (LT) (4-20 • C) tests (Table 1). In the high-temperature (HT) (10-30 • C) tests, PI 181410 seedlings were highly resistant (IT 2) to PSTv-4, PSTv-14, PSTv-40, and PSTv-51 but susceptible (IT 7) to PSTv-37. The adult-plants of PI 181410 showed high resistance (IT 2) to all the five tested races PSTv-4, PSTv-14, PSTv-37, PSTv-40, and PSTv-51 at both LT and HT. For comparison, the seedlings and adult plants of AvS were always highly susceptible Table 1. Stripe rust infection types of PI 181410 to five predominant or most virulent races of Puccinia striiformis f. sp. tritici in the United States tested in seedling and adult-plant stages at low (4-20 • C) and high (10-30 • C) temperature profiles.

Stripe Rust Phenotypes of the Recombinant Inbred Line (RIL) Population
The distributions of seedling IT among the 170 F 7 RILs tested with both PSTv-4 and PSTv-14 were continuous, and the patterns were similar ( Figure 1, Table S1). The F 5 , F 6 , and F 7 RILs were tested under natural infections of Pst in the fields at Mount Vernon and Pullman in 2015, 2016, and 2017, respectively. The IT distributions of adult plants among the 170 RILs in field environments were also continuous but the patterns are different ( Figure S2). Similarly, the relative area under the disease progress curve (rAUDPC) data also had continuous distributions with slightly different patterns ( Figure 2). These results indicated that the stripe rust resistance in PI 181410 was quantitatively inherited. However, the distribution patterns were similar between the races tested in the greenhouse and were different among the year-location environments.  The analysis of variance (ANOVA) showed that the genetic and environmental effects and their interaction were all significant (p < 0.001) ( Table 2). Broad-sense heritability (H 2 ) was estimated as 0.68 for the rAUDPC data and 0.78 for IT data, indicating that the resistance in PI 181410 was highly heritable. The correlation coefficients of either IT or rAUDPC for all possible pairwise environment comparisons were highly significant (p < 0.001), but the values ranged from 0.29 to 0.79 (Table 3). Correlation coefficients were much higher between any two greenhouse tests or field tests than those between a field test and a greenhouse test. The phenotypical analyses indicated that different QTL controlled the resistances observed in different environments and suggested that QTL analysis should be used to map the genes controlling the resistance in the RIL population.

QTL for ASR to Stripe Rust
Using the two sets of IT data of F 7 RILs at the seedling stage tested with races PSTv-4 and PSTv-14 in the greenhouse, we detected three QTL, QYrPI181410.wgp-4AS, QYrPI181410.wgp-4BL, and QYrPI181410.wgp-5BL.1, conferring resistance to both races (Table 4, Figure 3).  QYrPI181410.wgp-4AS was detected by SNP IWB9590 with the logarithm of odds (LOD) of 3.3 and 2.7 for PSTv-4 and PSTv-14, respectively (Table 4). It peaked at the genetic location of 43.0 cM and spanned a relatively small interval of 4.0 cM on the short arm of chromosome 4A. This QTL explained 11.2% and 5.4% of the phenotypic variation to PSTv-4 and PSTv-14, respectively.
QYrPI181410.wgp-4BL was detected by IWB6961 on 4BL using two sets of IT data with both PSTv-4 and PSTv-14 with the LOD ranging from 2.5 to 3.6 ( Table 4). It peaked at the location of 71.0 cM and spanned an interval of 20.0 cM region. This QTL explained 8.6% and 11.0% of the phenotypic   QYrPI181410.wgp-4AS was detected by SNP IWB9590 with the logarithm of odds (LOD) of 3.3 and 2.7 for PSTv-4 and PSTv-14, respectively (Table 4). It peaked at the genetic location of 43.0 cM and spanned a relatively small interval of 4.0 cM on the short arm of chromosome 4A. This QTL explained 11.2% and 5.4% of the phenotypic variation to PSTv-4 and PSTv-14, respectively.
QYrPI181410.wgp-4BL was detected by IWB6961 on 4BL using two sets of IT data with both PSTv-4 and PSTv-14 with the LOD ranging from 2.5 to 3.6 ( Table 4). It peaked at the location of 71.0 cM and spanned an interval of 20.0 cM region. This QTL explained 8.6% and 11.0% of the phenotypic variation to PSTv-4 and PSTv-14, respectively. QYrPI181410.wgp-5BL.1 was detected by IWB865 and peaked at the location of 211.0 cM with an interval of 10.0 cM on 5BL. It explained 15.1% and 14.0% of the phenotypic variation to PSTv-4 and PSTv-14, respectively.
A major QTL, QYrPI181410.wgp-4BL, was consistently detected using all 12 sets of phenotypic data from field tests and two sets of data from greenhouse tests. The LOD value of QYrPI181410.wgp-4BL ranged from 8.5 to 17.4 (Table 5, Figure 4). It peaked around the genetic position of 70.1 cM by SNP markers IWB70599 and IWB6961 within an interval of 3.0 cM (68.0 cM to 71.0 cM) on 4BL and explained 29.2% to 50.0% of the phenotypic variation. However, the phenotypic effect of QYrPI181410.wgp-4BL was lower (29.2 to 32.5) for infection type than for severity in Pullman, WA.
A major QTL, QYrPI181410.wgp-4BL, was consistently detected using all 12 sets of phenotypic data from field tests and two sets of data from greenhouse tests. The LOD value of QYrPI181410.wgp-4BL ranged from 8.5 to 17.4 (Table 5, Figure 4). It peaked around the genetic position of 70.1 cM by SNP markers IWB70599 and IWB6961 within an interval of 3.0 cM (68.0 cM to 71.0 cM) on 4BL and explained 29.2% to 50.0% of the phenotypic variation. However, the phenotypic effect of QYrPI181410.wgp-4BL was lower (29.2 to 32.5) for infection type than for severity in Pullman, WA.  Two minor QTL, QYrAvS.wgp-2BS and QYrAvS.wgp-2BL, were detected on chromosome 2B and the effects were provided by the susceptible parent AvS (Table 5, Figure 3). QYrAvS.wgp-2BS was detected with two sets of IT data from Mount Vernon and explained 5.0%~5.6% of IT variation. QYrAvS.wgp-2BL was only detected with one set of IT data from the Pullman location and explained 7.2% of IT variation. Both QTL did not have any significant effect on rAUDPC.
QYrPI181410.wgp-5AS was identified with two sets of SE data as a minor QTL conferring resistance on reducing severity of stripe rust in PI 181410 (Table 5). This QTL was detected by IWB47184 and located on the telomeric end of 5AS. It explained 3.9% to 4.2% of SE variation.
QYrPI181410.wgp-5BL.2 was detected with five sets of field IT data from Mount Vernon in 2015 and 2016 and Pullman in 2016, 2017, and 2018 (Table 5, Figure 3). It was located within a 7.0 cM genomic interval and peaked around 310.0 cM close to the end of 5BL. This QTL explained 4.8% and 7.9% of the IT variation of data from Mount Vernon. However, it explained 14.6% to 17.3% of IT variation of data from Pullman, WA. This result suggested that the effect of QYrPI181410.wgp-5BL.2 varies among the environments.

Effects of QTL Combinations
To estimate the effects of QTL combinations, the RILs were classified into different groups containing different QTL combinations based on the genotypes spanning the QTL intervals. According to the ASR QTL identified in the greenhouse seedling tests, RILs were classified into four groups ( Figure 5). Although these QTL contributed minor effects to reduce the stripe rust to specific races including PSTv-4 and PSTv-14 in the greenhouse conditions, the RILs containing more QTL showed higher resistance than those containing few QTL.
To estimate the effects of QTL combinations, the RILs were classified into different groups containing different QTL combinations based on the genotypes spanning the QTL intervals. According to the ASR QTL identified in the greenhouse seedling tests, RILs were classified into four groups ( Figure 5). Although these QTL contributed minor effects to reduce the stripe rust to specific races including PSTv-4 and PSTv-14 in the greenhouse conditions, the RILs containing more QTL showed higher resistance than those containing few QTL. The similar trends were also found for HTAP resistance. Four QTL for HTAP resistance were identified in PI 181410. Based on the markers for these QTL, RILs were classified into five groups ( Figure 6). QYrPI181410.wgp-4BL had a higher effect than QYrPI181410.wgp-1BL, QYrPI181410.wgp-5AS, and QYrPI181410.wgp-5BL.2. The resistance level was greatly increased when the major QTL QYrPI181410.wgp-4BL was in combination with minor QTL QYrPI181410.wgp-1BL and/or QYrPI181410.wgp-5BL.2. The RILs containing one of the QTL showed a lower IT and rAUDPC than those without any QTL. The RILs containing more resistance QTL had a higher level of resistance (lower IT and/or rAUDPC values) than those containing a smaller number of QTL ( Figure 6A,B). These results indicated that the combinations of different QTL increased the level of resistance through additive interactions. The similar trends were also found for HTAP resistance. Four QTL for HTAP resistance were identified in PI 181410. Based on the markers for these QTL, RILs were classified into five groups ( Figure 6). QYrPI181410.wgp-4BL had a higher effect than QYrPI181410.wgp-1BL, QYrPI181410.wgp-5AS, and QYrPI181410.wgp-5BL.2. The resistance level was greatly increased when the major QTL QYrPI181410.wgp-4BL was in combination with minor QTL QYrPI181410.wgp-1BL and/or QYrPI181410.wgp-5BL.2. The RILs containing one of the QTL showed a lower IT and rAUDPC than those without any QTL. The RILs containing more resistance QTL had a higher level of resistance (lower IT and/or rAUDPC values) than those containing a smaller number of QTL ( Figure 6A,B). These results indicated that the combinations of different QTL increased the level of resistance through additive interactions.

Kompetitive Allele Specific PCR (KASP) Markers for the Major-Effect QTL
As QYrPI181410.wgp-4BL was the major-effect QTL for HTAP resistance identified in PI 181410 and detected in all greenhouse and field tests, KASP markers were developed for this QTL. Five KASP markers were successfully converted from their original SNP markers ( Table 6). Three of five KASP markers, KIWB81016, KIWB6961, and KIWB80, which were developed from SNP markers IWB81016, IWB6961, and IWB80, respectively, were confirmed to be polymorphic between PI 181410 and AvS. When tested on the RIL population, the KASP markers produced the same genotypes as their original SNP markers among the individual RILs.
The three polymorphic KASP markers, were used to determine their polymorphisms in 155 wheat cultivars and breeding lines from the US (Table S3). These markers had polymorphism rates 45.8%, 69.0%, and 100%, respectively. As expected, none of the 155 wheat accessions had the same marker haplotype of PI 181410 for the three KASP markers. This result indicated that QYrPI181410.wgp-4BL was absent in the collection of US wheat cultivars and breeding lines, and the KASP markers in combination, or individual marker KIWB80, should be useful in MAS for incorporating this QTL into new wheat cultivars.

Kompetitive Allele Specific PCR (KASP) Markers for the Major-Effect QTL
As QYrPI181410.wgp-4BL was the major-effect QTL for HTAP resistance identified in PI 181410 and detected in all greenhouse and field tests, KASP markers were developed for this QTL. Five KASP markers were successfully converted from their original SNP markers ( Table 6). Three of five KASP markers, KIWB81016, KIWB6961, and KIWB80, which were developed from SNP markers IWB81016, IWB6961, and IWB80, respectively, were confirmed to be polymorphic between PI 181410 and AvS. When tested on the RIL population, the KASP markers produced the same genotypes as their original SNP markers among the individual RILs.
The three polymorphic KASP markers, were used to determine their polymorphisms in 155 wheat cultivars and breeding lines from the US (Table S3). These markers had polymorphism rates 45.8%, 69.0%, and 100%, respectively. As expected, none of the 155 wheat accessions had the same marker haplotype of PI 181410 for the three KASP markers. This result indicated that QYrPI181410.wgp-4BL was absent in the collection of US wheat cultivars and breeding lines, and the KASP markers in combination, or individual marker KIWB80, should be useful in MAS for incorporating this QTL into new wheat cultivars.

Discussion
Wheat landrace PI 181410 has been consistently showing high-level resistance to stripe rust in the US Pacific Northwest region since 2006. The four-way tests conducted with predominant and most virulent races of Pst in both seedling and adult-plant stages at both LT and HT profiles showed that PI 181410 has both race-specific ASR and high-level HTAP resistance. Using the high-density genetic maps and multiple environmental tests for QTL mapping, we identified three QTL in PI 181410 for ASR, QYrPI181410.wgp-4AS, QYrPI181410.wgp-4BL, and QYrPI181410.wgp-5BL.1. In addition to the QTL on 4AS, we identified three more QTL, designated as QYrPI181410.wgp-1BL, QYrPI181410.wgp-5AS, and QYrPI181410.wgp-5BL.2 in PI 181410 conferring resistance observed at the adult-plant stage in the fields. Interestingly, we also detected two QTL, designated as QYrAvS.wgp-2BS and QYrAvS.wgp-2BL, from AvS with minor effects only on IT at Mount Vernon and Pullman, respectively. Each of the QTL were located within a relatively small genomic region.
Two minor-effect QTL, QYrAvS.wgp-2BS and QYrAvS.wgp-2BL, were identified from susceptible parent AvS. QYrPI181410.wgp-2BS was mapped close to the distal end of the short arm. According to previous studies, two QTL including QYr.inra-2BS [22] and QYr.wgp-2B.1 [23] were mapped at the same region as QYrAvS.wgp-2BS. QYr.inra-2BS contributed moderate resistance to stripe rust at both seedling and adult-plant stages [22]. Identified in several US Pacific Northwest wheat cultivars by genome-wide association mapping, QYr.wgp-2B.1 contributed adult-plant resistance to stripe rust and was significantly associated with both IT and SE in several environments including Pullman and Mount Vernon, Washington [23]. However, QYrPI181410.wgp-2BS was only associated with IT in Mount Vernon. Thus, QYrAvS.wgp-2BS is likely a locus different from QYr.inra-2BS and QYr.wgp-2B.1. QYrAvS.wgp-2BL was only detected using one set of IT data in a single field condition in Pullman. It was located at the distal end of 2BL. Rosewarne et al. [17] reported a minor QTL affecting adult-plant resistance in AvS flanked by simple sequence repeat (SSR) markers Xgwm1027 and Xgwm619 at the same location on chromosome 2BL as QYrAvS.wgp-2BL. Since their QTL and QYrAvS.wgp-2BL were detected in the same wheat genotype in the same chromosomal location, they should be the same gene.
QYrPI181410.wgp-5AS had a minor effect on reducing disease severity in PI 181410 and was located at the telomere of chromosome 5AS. As no stripe rust resistance QTL have been reported at this location, the 5AS QTL is a new QTL for resistance to stripe rust.
Among the QTL identified in PI 181410, QYrPI181410.wgp-4BL was the most effective. It was consistently detected using all 16 sets of phenotypic data collected from all field tests and both seedling and adult-plant tests in the greenhouse. Its effect was much stronger in the adult-plant stage in the fields and at HT in the greenhouse test than the seedling stage at LT, and therefore, it should be considered as a HTAP resistance gene. QYrPI181410.wgp-4BL was mapped to an interval of 9.0 cM (67.0~76.0 cM) close to the centromeric region of 4BL according to the consensus map [29,30]. The overlap of the genomic regions detected by all 16 data sets narrowed the QTL to 4.0 cM between 68.0 and 72.0 cM positions (Table 5, Figure 4). In previous studies, Yr50 and Yr62 were also mapped to the same location on 4BL as QYrPI181410.wgp-4BL. The dominant resistance gene, Yr50, derived from Thinopyrum intermedium, was mapped on 4BL and flanked by SSR markers Xbarc1096 and Xwmc47, and Yr50 exhibited ASR to all tested Chinese Pst races and eight US predominant Pst races including PSTv-40 [31,32]. In contrast, QYrPI181410.wgp-4BL is from common wheat landrace PI 181410 and confers strong HTAP resistance with only a minor effect at seedling stage. The differences clearly indicate that QYrPI181410.wgp-4BL is different from Yr50. Yr62, from Portuguese spring wheat variety PI 192252, was flanked by SSR markers Xgwm192 and Xgwm251 on 4BL [32]. This gene provides a higher level HTAP resistance (40%-60% of phenotypic variation) than QYrPI181410.wgp-4BL (29.2%-46.9%). Yr62 is not effective at the seedling stage. Thus, QYrPI181410.wgp-4BL is likely different from Yr62. Further studies are needed to determine their relationship using the allelism test approach.
Two major QTL, QPst.jic-4B and QYrva.vt-4BL, were also mapped around the centromeric region on 4BL [21,23]. QPst.jic-4B, identified from a German winter wheat cultivar Alcedo, was flanked by amplified fragment length polymorphism (AFLP) marker S24M37_130 and SSR marker Xcfd39 and peaked at Xwmc692 on 4BL [21]. QYrva.vt-4BL, identified from a soft red winter wheat VA00W-38, was flanked by SSR markers Xwmc692 and Xgwm149 on 4BL [33]. Both QPst.jic-4B and QYrva.vt-4BL explained less effect (17.9%-28.9% and 19.3%, respectively) of the phenotypic variation. Thus, QYrPI181410.wgp-4BL is unlikely the same as any of these two QTL in winter wheat cultivars. As QYrPI181410.wgp-4BL is the most useful among the QTL detected in the present study, we developed three KASP markers for this QTL. The KASP markers performed in the same way as the original SNP markers between the parents and among the RILs, and they were highly polymorphic when tested with 155 wheat cultivars and breeding lines from the US Pacific Northwest. None of the wheat materials have the KIWB80 allele or the PI 181410 haplotype of the three markers. Therefore, the KASP marker KIWB80 or the combination of the three markers can be used in MAS for incorporating the effective QTL into various wheat backgrounds for developing new cultivars with durable, high-level resistance to stripe rust.

Plant Materials
PI 181410, an Afghanistan spring wheat landrace with resistance to stripe rust, was used as the male parent in crossing with Avocet S (AvS), a stripe rust susceptible Australian spring wheat genotype, as a female parent. A recombinant inbred line (RIL) population of the cross was developed using the single-seed decent method. A total of 170 RILs of the F 5 , F 6 , F 7 , and F 8 generations were used for phenotyping the response to stripe rust in different years and locations, and the first 114 F 7 RILs of the F 7 generation were used for SNP genotyping. A collection of 155 US Pacific Northwest spring and winter wheat cultivars and breeding lines was used to determine polymorphisms of the KASP markers developed in this study.  [9,10], following the procedures as previously described [7,34]. Briefly, about 10-15 seeds of each line were planted in one pot in a plastic tray (53.43 × 27.5 × 6.10 cm, containing 72 pots) filled with soil mixture. Seedlings at two-leaf stage [growth stage (GS) 10] [35] were inoculated with urediniospores of a selected race mixed with talc at the ratio of 1:20. The inoculated plants were incubated in a dew chamber at 10 • C in darkness for 24 h for urediniospores to germinate and infect the plants, and then transferred to a growth chamber at a diurnal cycle of low-temperature (LT) gradually changing from 4 • C at 2:00 a.m. to 20 • C at 2:00 p.m. with a 16 h light/8 h dark photoperiod. For high-temperature (HT) tests, after incubating in the dew chamber for 24 h at 10 • C in darkness, the inoculated plants were placed in a growth chamber at a diurnal temperature cycle gradually changing from 10 • C at 2:00 a.m. to 30 • C at 2:00 p.m. with the same photoperiod cycle. For adult-plant tests at both LT and HT profiles, 3 plants grown in a pot (15 × 15 × 18 cm) to the boot stage (GS 45) were inoculated, incubated, and grown under the same conditions as described above. For each race-temperature test, a total of 9 plants in 3 pots were evaluated. Stripe rust infection type (IT) of each plant was scored at 18 to 21 days after inoculation using the 0 to 9 scales [36]. For the seedling tests, ITs of the second leaves were recorded while those of the flag leaves were recorded for the adult-plant tests. Seedlings or adult plants of AvS were included as the susceptible check in each stage-temperature test.

Evaluation of PI 181410 and RILs for Stripe Rust Reaction in the Greenhouse
Races PSTv-4 and PSTv-14 of Pst were used for testing the 170 F 7 RILs at the seedling stage at the LT profile (4-20 • C) with three biological replicates based on their relative avirulence on the seedlings of PI 181410. In these tests, PI 181410 and AvS were included as the resistant and susceptible checks, respectively, and the 18 Yr single-gene differentials [9] were included to ensure the correct race. PSTv-37 and PSTv-40 were used to test the adult plants of 170 RILs at the HT profile (10-30 • C) based on their virulence on seedlings at LT and avirulence on adult plants of PI 181410 at HT. PI 181410 and AvS were also included as the resistant and susceptible checks, respectively. The inoculation, incubation, growth of plants, and IT data record were the same as described above for the seedling LT and adult plant HT tests for the parental lines.

Evaluation for Stripe Rust Responses of PI 181410 and RILs in the Fields
The stripe rust responses of PI 181410, together with other wheat genotypes, were tested under natural Pst infections in field nurseries at Mount Vernon in northwestern Washington and Pullman in southeastern Washington since 2006. The two locations are about 500 km apart with different climate conditions and different Pst races [9,10]. The 170 F 5 , F 6 , F 7 , or F 8 RILs were consecutively phenotyped in both locations in 2015, 2016, 2017, and 2018, respectively. The parental lines and RILs were planted in a randomized block design with three replications. About 20-30 seeds of each line were planted in a 50 cm row with 20 cm between rows. As a susceptible check and spreader, AvS was planted every 20 rows and around the field. Stripe rust IT and severity (SE) were recorded three times at the tillering stage (GS 30-31), booting stage (GS 45), and flowering to grain-filling stages (GS 60-70). Severity was recorded as the percentage of leaf area infected. The SE data were used to calculate the area under the disease progress curve (AUDPC) [37], which was then converted to relative AUDPC (rAUDPC) as percentage of the mean AUDPAC (as 100%) of susceptible parent AvS.
The analysis of variance (ANOVA) was performed for the IT data, SE, and rAUDPC using the aov function of the R program (version 3.5.2, https://www.r-project.org) to determine the effects of genetic and environmental factors and their interactions. Correlation coefficients were calculated using the correlation function in Excel for rAUDPC and IT data set to compare population responses across the four field environments and two greenhouse tests.

Genotyping the Parental Lines and RILs
The parental lines and F 7 RILs were planted in the greenhouse for DNA extraction. Genomic DNA was extracted from the leaves harvested from a single plant of each F 7 RIL using the CTAB (cetyl trimethylammonium bromide) method [38]. DNA quality was determined using 0.8% agarose gel electrophoresis, and DNA concentration was measured using a NanoDrop ND-1000 (Thermo Scientific, Wilmington, DE, USA). Each DNA sample was diluted to 50 ng/µL for use in genotyping. Out of the 170 F 7 RILs, the first 114 RILs and parental lines were genotyped using the 90K Illumina ® iSelect wheat SNP assay [39] at the USDA-ARS Biosciences Research Laboratory, Fargo, ND. The reduced number of RILs genotyped was due to a budget limitation. SNP calling and clustering were determined using the functions of software GenomeStudio v2011.1 (Illumina, San Diego, CA, USA) based on the described genotypic data. Monomorphic SNPs, SNPs with more than 10% missing data, and RILs with 30% missing data were removed before further analyses.

Genetic Map Construction and QTL Mapping
SNP markers polymorphic between the parents were tested for the goodness of the fit to the 1:1 segregation ratio among the RILs using the chi-squared test in JoinMap 4.1 (https://www.kyazma.nl/ index.php/JoinMap/Manual/). SNP markers that did not fit the 1:1 ratio (p < 0.01) were excluded for further analyses.
A genetic map was constructed using software JoinMap 4.1 [40] with the Kosambi mapping function [41]. The threshold value of logarithm of odds (LOD) score of 3.0 was applied to select the genetic linkage groups. The corresponding chromosome of each linkage group was determined based on the chromosome locations of SNPs from the 90K array [39]. The graphic genetic maps were generated using software MapChart v2.3 [42].
The IT data of both seedling stage and adult-plant stage in each year-location environment were used directly for QTL analysis. rAUDPC data were used for QTL mapping if the test had three sets of SE data. If only one set of SE data in the adult-plant stage was collected, as in the case of 2018 Pullman location, the SE data were also directly used for QTL mapping. QTL analysis was performed using both single-marker analysis (SMA) and composite interval mapping (CIM) with the Haley-Knott regression (HK) method in the R/qtl package of RStudio version 1.0.143 [43]. For QTL mapping, the probabilities of association between the genotypes and phenotypes were calculated using the function calc.genoprob with a step size of 1 cM distance and a genotype error of 0.01. SMA was performed using the scanone function and followed by CIM using the cim function. For each set of phenotypic data, a permutation test was performed using the n.perm function (n = 1000) to obtain a genome-wide LOD significance threshold. The resistance QTL interval with 95% confidence was obtained using the bayesint function. SNP markers closest to the QTL peak were identified using find.marker. The effect of each QTL was estimated using the plotPXG function. For the phenotypes with more than one QTL identified, the effect of each QTL was estimated using the multiple QTL mapping method in the R/qtl package [44]. The additive effect for each QTL was estimated using the get.ests = TRUE function. The interaction between two QTL was calculated using the function addint in the multiple QTL mapping method.

Developing KASP Markers
SNP markers flanking or at the peaks of the identified resistance QTL were selected for developing KASP markers using the information provided by PolyMarker (http://polymarker.tgac.ac.uk). Primers of KASP markers were synthesized by Sigma-Aldrich Inc. (Saint Louis, MO, USA). The KASP markers were validated using the parental lines PI 181410 and AvS and the RILs. The PCR reaction, amplification, and detection for KASP markers were conducted following the previously described method [45]. The PCR assay was performed in a LightCycler ® 480 System (Roche Applied Science, Indianapolis, IN, USA) using the endpoint genotyping analysis. To determine the polymorphisms of the KASP markers in various wheat backgrounds, 155 spring and winter wheat cultivars and breeding lines from the US Pacific Northwest were tested with the markers for the major effective QTL.

Conclusions
This study determined that the spring wheat landrace PI 181410, originally from Afghanistan, possesses both race-specific ASR and nonrace-specific HTAP resistance to stripe rust. Three QTL conferring ASR in PI 181410 were mapped to chromosomes 4BS, 4BL, and 5BL; and four QTL (including the 4BL QTL detected in the seedling stage) conferring HTAP resistance were mapped on chromosomes 1BL, 4BL, 5AS, and 5BL. The 1BL QTL was identified as Yr29, and the four other QTL in PI 181410 are potentially new. Two QTL were mapped on 2BS and 2BL from the susceptible parent AvS with minor effects on resistance in the field. The 2BL QTL is the same as a previously mapped QTL in the same variety, while the 2BS QTL is a potentially new locus. As the 4BL QTL was the most effective among the identified QTL and potentially new, three KASP markers (KIWB81016, KIWB6961, and KIWB80) were developed from its associated SNPs. These markers, especially KIWB80, were highly polymorphic among the US Pacific Northwest wheat cultivars and breeding lines. The resistant germplasm, identified QTL, and KASP markers for the most effective QTL are useful for developing wheat cultivars with durable, high-level resistance to stripe rust.