Mapping a Large Number of QTL for Durable Resistance to Stripe Rust in Winter Wheat Druchamp Using SSR and SNP Markers

Winter wheat Druchamp has both high-temperature adult-plant (HTAP) resistance and all-stage resistance to stripe rust caused by Puccinia striiformis f. sp. tritici (Pst). The HTAP resistance in Druchamp is durable as the variety has been resistant in adult-plant stage since it was introduced from France to the United States in late 1940s. To map the quantitative trait loci (QTL) for stripe rust resistance, an F8 recombinant inbred line (RIL) population from cross Druchamp × Michigan Amber was phenotyped for stripe rust response in multiple years in fields under natural infection and with selected Pst races under controlled greenhouse conditions, and genotyped with simple sequence repeat (SSR) and single nucleotide polymorphism (SNP) markers. Composite interval mapping (CIM) identified eight HTAP resistance QTL and three all-stage resistance QTL. Among the eight HTAP resistance QTL, QYrdr.wgp-1BL.2 (explaining 2.36-31.04% variation), QYrdr.wgp-2BL (2.81–15.65%), QYrdr.wgp-5AL (2.27–17.22%) and QYrdr.wgp-5BL.2 (2.42–15.13%) were significant in all tests; and QYrdr.wgp-1BL.1 (1.94–10.19%), QYrdr.wgp-1DS (2.04–27.24%), QYrdr.wgp-3AL (1.78–13.85%) and QYrdr.wgp-6BL.2 (1.69–33.71%) were significant in some of the tests. The three all-stage resistance QTL, QYrdr.wgp-5BL.1 (5.47–36.04%), QYrdr.wgp-5DL (9.27–11.94%) and QYrdr.wgp-6BL.1 (13.07-20.36%), were detected based on reactions in the seedlings tested with certain Pst races. Among the eleven QTL detected in Druchamp, at least three (QYrdr.wgp-5DL for race-specific all-stage resistance and QYrdr.wgp-3AL and QYrdr.wgp-6BL.2 for race non-specific HTAP resistance) are new. All these QTL, especially those for durable HTAP resistance, and their closely linked molecular markers could be useful for developing wheat cultivars with durable resistance to stripe rust.


Introduction
Wheat (Triticum aestivum L.) is one of the most important crops worldwide, and one of its most damaging diseases is stripe rust, caused by the fungus Puccinia striiformis f. sp. tritici (Pst) [1][2][3]. Infection by Pst on wheat plants can occur anytime from the one-leaf stage to prior to plant maturity [1]. Fungicide application can control stripe rust. However the use of fungicides adds extra cost for wheat production; can have adverse effects on the environment; and may not be fully effective if not applied on time. Growing cultivars with genetic resistance is the most effective, economical, easy to use and environment-friendly way to control the disease [4], [5].
Stripe rust resistance can be classified as race specific resistance and race non-specific resistance based on specificity, and also can be classified as all-stage resistance (also known as seedling resistance) and adult plant resistance (APR) based on the plant growth stage [1], [4], [6]. All-stage resistance, usually race specific and complete, can be detected at seedling stage and protects plants throughout the entire growth cycle when effective against races. However, new virulent races may overcome race specific resistance. In contrast, APR, usually non-race specific, expresses at adult-plant stage. High-temperature adult-plant (HTAP) resistance is further characterized by effectiveness at late stages of plant growth when the weather becomes warm. It is important to know if the APR in a particular cultivar is sensitive to temperature or not, as the information is useful for determining where and under which weather conditions the APR can be effective or not effective, and also for understanding if it is race specific or race nonspecific. The level of HTAP resistance is often incomplete and is affected by the plant growth stage, temperature, humidity and inoculum load [1], [4], [7], [8]. Cultivars with HTAP resistance are susceptible at seedling stage if they do not have effective all-stage resistance, but express resistance at adult-plant stage usually after jointing stage. Many wheat cultivars have all-stage and/or HTAP resistance. Identification and mapping of genes conferring HTAP and all-stage resistance can allow breeders to combine different types of resistance into single genetic backgrounds to achieve complete and long-lasting protection of cultivars from stripe rust [1], [4], [5].
Druchamp (PI 174622) is a soft white winter wheat cultivar developed in Ville-de-Paris, France in 1940 and introduced to the US Pacific Northwest for production in 1949 (http:// www.ars-grin.gov/cgi-bin/npgs/acc/display.pl?1151600). Although Druchamp is no longer grown for commercial production because other cultivars have higher yields, it has been used to monitor virulence in the Pst population. Because of its race specific all-stage resistance, Duchamp has been used in seedling stage to differentiate Pst races since 1969 [9]. Later, Druchamp was also found to have high level HTAP resistance [10][11][12][13]. A previous monosomic study located three genes for race-specific all-stage resistance to chromosomes 1B, 5B and 6A [14], but the chromosomal locations of these genes have not been determined. Using a biometric approach, three genes were estimated to confer HTAP resistance in Druchamp and significant additive and dominant components and complex gene interactions were observed among resistance gene loci [12], [13]. However, none of the HTAP resistance genes in Druchamp had been mapped to chromosomes prior to the present study.
The objectives of the present study were to map genes or quantitative trait loci (QTL) for either all-stage or HTAP resistance using molecular markers and to assess the specificity or stability of the effects of identified QTL for different types of resistance across multiple environments or tested with different races. This study is expected to generate information vital for understanding durability of HTAP resistance and identify markers that can be used to incorporate HTAP resistance genes from Druchamp into new wheat cultivars.

Ethics statement
No permits were necessary to conduct reported field experiments, because Pst is a naturally occurring plant pathogen in the reported environments and no exotic cultures of the pathogen were used. Research was conducted on land owned by Washington State University. No protected species were sampled. No animal subjects were used in described research. All experiments reported in this manuscript comply with all federal, state and university rules and regulations.

Plant and pathogen materials
Druchamp (resistant parent) used as the female parent was crossed with susceptible winter wheat Michigan Amber, and 94 F 8 recombinant inbred lines (RILs) were obtained from 94 F 2 plants of a single F 1 plant of the cross through single-seed descent [12], [13]. The RILs were phenotyped in fields under natural Pst infection and in greenhouses with selected races, and also genotyped with DNA markers. Seven Pst races, PST-25, PST-29, PST-35, PST-45, PST-100, PST-114 and PST-127, were chosen based on the reactions of Druchamp to these races and their virulence formulae and predominance in different periods over the past 40 years [1], [9], [15], [16].

Greenhouse tests for race specific all-stage resistance
Seedling tests were conducted under controlled conditions in a greenhouse as previously described [10], [11]. About 10 seeds of each line were planted in a 7×7×7 cm pot filled with soil mixture and grown in a rust-free greenhouse. Two-leaf stage seedlings were inoculated with urediniospores of a selected race. Five (PST-29, PST-35, PST-45, PST-100 and PST-114) of the seven selected races were used in the seedling tests as they are avirulent to Druchamp and virulent to Michigan Amber [1], [9], [16]. The inoculated plants were kept in a dew chamber for 24 h at 10°C without light, and then grown in a growth chamber using a low diurnal temperature cycle gradually changing between 4°C at 2:00 am and 20°C at 2:00 pm with 16 h light/8 h dark [10], [11]. A set of wheat varieties that were used to differentiate Pst races was included in each race test to confirm the identity of the race [15], [16]. Infection type (IT) based on the 0-9 scale [9] was scored for each line 18 to 21 days after inoculation when stripe rust was fully developed on Michigan Amber.

Greenhouse tests for HTAP resistance
Adult-plants of the F 8 RILs and parents were evaluated for studying HTAP resistance at high temperatures in the greenhouse with two (PST-25 and PST-127) of the seven selected races, as the two races are virulent on seedlings of Druchamp and both seedling and adult plants of Michigan Amber, but not virulent to adult plants of Druchamp. In addition, PST-25 and PST-127 represent predominant races in 1980s and recent years, respectively [9], [16]. One-leaf stage seedlings were vernalized in a growth chamber at 2-5°C for 40 days. After vernalization, twelve seedlings for each line and each race test were transplanted in three pots of 15-cm in diameter filled with soil mixture. The pots were arranged using a completely randomized block design and grown in a rust-free greenhouse. Plants were inoculated at the heading stage, kept in a dew chamber for 24 h at 10°C without light, and then grown in a growth chamber under conditions similar to the seedling tests, except at a higher diurnal temperature cycle gradually changing between 10°C at 2:00 am and 30°C at 2:00 pm [4], [12]. Data of IT and disease severity (DS, percentage of diseased foliage) were scored as average for each plant 20-22 days after inoculation.

Field tests
The RIL population and parents were evaluated for stripe rust response in the field nurseries at The two locations are about 500 km apart and have different Pst race compositions and climatic conditions. In each field experiment, the F 8 RIL population and parents were planted in a randomized complete block design with three replications. About 30 seeds for each RIL or parent were planted in a 60-cm row with 20 cm between rows. Susceptible variety PS 279 were planted every 20 rows throughout the field to increase the speed and uniformity of stripe rust development. The cultural practices commonly used in wheat production of these regions were performed for fertilization and weed control. All field experiments were conducted under natural infection of Pst to allow the evaluation of responses to various races since the natural occurrence of the disease is adequate [1], [17]. Data of IT and DS were visually recorded for each row three times at the heading, flowering and milk stages at Pullman and at the jointing, heading and milk stages at the Mt. Vernon location when Michigan Amber had approximately 50, 80 and >95% DS, respectively. The three-time DS values were used to calculate the value of relative area under the disease progress curve (rAUDPC) for each RIL as previously described [12], [18]. Both the rAUDPC and IT data were used in QTL mapping.

DNA extraction
Genomic DNA was extracted from the leaf samples using the CTAB method [19]. DNA concentrations were determined using agarose gel electrophoresis [20] and spectrophotometry (NanoDrop ND-1000, Thermo Scientific, Wilmington, DE, USA).

SSR marker analysis
More than 700 SSR markers were screened to identify those polymorphic between the parents. The polymorphic markers were used to genotype the RILs. The distributions of the markers on the 21 wheat chromosomes were determined based on Somers et al. [21] and the GrainGenes database (http://www.wheat.pw.usda.gov). The M13 tail (5'-CACGACGTTGTAAAACGAC) was added to the 5' end of each forward primer to detect polymerase chain reaction (PCR) products through direct labeling [22]. The M13 universal primers were labeled with one of the fluorescent dyes FAM (blue), VIC (green), NED (yellow) and PET (red) (Applied Biosystems, Foster City, CA, USA) for detecting different PCR products. Each 12-μl PCR reaction mix con- PCR products with four different fluorescent dyes were pooled together including 3, 3, 4 and 6 μl of FAM, VIC, NET and PET, respectively and added ddH 2 O to 25 μl, from which 3 μl was transferred into a new tube. A volume of 9 μl Hi-Di Formamide and 1 μl 445-bp Cassul DNA ladder (Applied Biosystems) were added to the tube, giving a total of 13 μl. After denaturing at 95°C for 5 min, the 13 μl mixture was subjected to capillary electrophoresis using an ABI3730 DNA Analyzer (Applied Biosystems, Foster City, CA). Marker alleles were scored using software GeneMapper v1.5 (Softgenetics, State College, PA, USA).

SNP marker analysis
The RIL population and parents were genotyped with single nucleotide polymorphism (SNP) markers using the Illumina Infinium assay and the Wheat SNP 9K iSelect BeadChips developed by the International Wheat SNP Consortium [23]. SNP genotyping was performed on BeadStation and iScan instruments at the USDA-ARS Biosciences Research Laboratory, Fargo, ND, USA. The raw SNP data were processed with the Illumina GenomeStudio v2011.1 software (Ilumina Inc, San Diego, CA).

Map construction and QTL analysis
The SSR and SNP data were used to construct the linkage groups using JoinMap version 4.0 [24]. Genetic distances were calculated using the Kosambi mapping function [25]. For each linkage group, the SSR marker order and the assignment to chromosomes were based on the wheat maps published by Somers et al. [21], Sourdille et al. [26] and the GrainGenes database (http://wheat.pw.usda.gov). The chromosomal positions of SNP markers were determined using the genetic maps developed by Cavanagh et al. [23]. The MapChart computer program [27] was used to draw the linkage maps.
For the QTL analysis, different locations and different years of the same location were considered as different environments as the weather conditions, the time and speed of stripe rust development and the race compositions could be different from location to location and from year to years. QTL analysis was performed using mean rAUDPC and IT values of each environment, three-year (2006, 2010 and 2011) means of Pullman, three-year (2005, 2010 and 2011) means of Mt. Vernon experiments and also means of all six environments. QTL mapping was conducted using the composite the interval mapping (CIM) program [28], [29] in the WinQTL Cartographer v2.5 software [30]. The likelihood odds (LOD) thresholds for determining statistically significant QTL were calculated by 1,000 permutations [31]. Based on the permutation tests, LOD 3.0 was set as the threshold to determine HTAP resistance QTL and 5.0 to determine all-stage resistance QTL. A walk speed of 0.5 cM was used for all QTL detections. LOD, additive effects (a) and phenotypic coefficients of determination (R 2 ) for individual QTL were calculated using CIM [32].

Statistical analysis of phenotypic data
Analysis of variance (ANOVA) was performed using the rAUDPC and IT data to determine the effects of genetic and environmental factors and their interaction using the SAS statistics package (SAS Institute, Inc., Cary, NC, USA). The PROC GLM procedure was used to test lines as a fixed effect, and environments including combination of locations and years and replicates as random effects. The variance components were determined based on ANOVA for a random model generated from PROC GLM. The broad-sense heritability (H 2 ) was estimated based on the formula H 2 = σ 2 g / σ 2 p , where σ 2 g is the genetic variance and σ 2 p represents the phenotypic variance. The genetic variance (σ 2 g ) was calculated from (σ 2 L -σ 2 E )/r, where σ 2 L is the mean variance of the RILs, σ 2 E is the error variance and r equals the number of replications [33]. The correlation coefficients were calculated for pairwise comparison of the population responses in the six environments.

Race specific all-stage resistance
When tested with the five races avirulent on Duchamp at the seedling stage and low temperature cycle, Druchamp was resistant (IT 2) and Michigan Amber was susceptible (IT 8) ( Table 1). The F 8 RIL population showed a continuous segregation with various IT. Arbitrary classification of the lines into resistant (IT 2-5) and susceptible (IT 6-8) groups suggested the presence of one gene for resistance in Druchamp to races PST-29 and PST-114, and two genes for resistance to PST-35, PST-45 and PST-100 (Table 1). In each test, not all of the RILs were highly resistant (IT 2) or highly susceptible (IT 8), indicating that the race-specific all-stage resistance was mainly quantitative. Thus, QTL analysis was suitable for mapping the genes conferring the race-specific all-stage resistance in Druchamp. The ANOVA results showed significant (P < 0.0001) genetic variations among RILs for both rAUDPC and IT in the field and greenhouse experiments ( Table 2). No significant variation was detected among the replications within each experiment (P = 0.16-0.79). The estimated broad-sense heritability values based on all data sets for rAUDPC and IT ranged from 0.66 to 0.94 (Table 2). Correlation coefficients, ranging from 0.43 to 0.95, for either rAUDPC or IT Table 1. Seedling infection types (ITs) of the parents and number of F 8 recombinant inbred lines (RILs) from the Druchamp × Michigan Amber-derived recombinant inbred lines (RILs) tested in the greenhouse with races of Puccinia striiformis f. sp. tritici, theoretical segregation ratios of resistant and susceptible lines, and χ 2 and P values of chi-squared tests for goodness of fit of the observed numbers and expected ratios.

IT b
No. among the six field environments and two race tests in the greenhouse were all significant (P < 0.001) ( Table 3), suggesting that the expression of HTAP resistance was consistent across the different environments and against different races.

Construction of linkage maps
Of the 768 SSR markers screened, 240 (31.3%) were found to be polymorphic between the two parents and were used to test the RIL population. Based on the chi-squared tests, 156 were considered reliable for mapping as they fit the expected 1:1 ratio and thus used in constructing linkage groups. Additionally, a total of 9,000 SNP markers were evaluated on the population, of which 2,535 were found to be polymorphic among the RILs and suitable for linkage construction. Using the 156 SSR and 2,535 SNP markers, 32 linkage groups consisting of 132 SSR and 2,300 SNP markers were constructed, and the remaining 24 SSR markers and 235 SNPs were unlinked based on the LOD 3.0 set by the permutation test. The 32 linkage groups were assigned to 18 wheat chromosomes. Chromosomes 1D, 2B, 3A, 4B, 5A, 5B and 7B each    (Table 4). All three QTL were contributed by Druchamp. The numbers of the QTL was consistent with those determined through classic genetic analysis by arbitrarily classifying the RILs based on their IT for all five race tests ( Table 1). The various infection types observed among the RILs indicated effects of mostly additive and some epistatic interactions. The QTL analysis detected more significant effects of additive (P < 0.001) then epistatic interactions (P = 0.05) when two QTL were detected.   QYrdr.wgp-5BL.1, most closely linked to SNP marker IWA8581, explained 36.04, 30.57 and 5.47% of the phenotypic variations in the tests with PST-29, PST-45 and PST-100, respectively (Fig 3; Table 4). QYrdr.wgp-5DL, close to IWA8331, explained 11.94% in the test with PST-35
Among the eight QTL detected for HTAP resistance, QYrdr.wgp-1BL.2 was the most consistent and provided the highest level of resistance. It was detected in all six field and two greenhouse experiments, and also detected with Pullman and Mt. Vernon mean data and the overall mean data of all six field experiments (Tables 5 and 6). This QTL explained 2.36-31.04% of the phenotypic variance, depending upon the experiment. The closest marker to the QTL was SNP  marker IWA8581, and SSR markers Xwmc830 and Xgwm259 flanked the locus by approximately 4.6 and 2.0 cM, respectively (Fig 4).
QYrdr.wgp-6BL.2, most closely associated with SNP marker IWA6420 (Fig 4) and explained 1.69-33.71% of the phenotypic variation, was detected in all experiments, but only the LOD values from the experiments of Pullman 2010, Pullman 2011, Mt. Vernon 2011 and race PST-127 in the greenhouse and the Pullman three-year mean data were greater than the threshold value of 3.0, indicating that relatively low stability compared to the other QTL.
QYrdr.wgp-1BL.1 was detected in four of the six field experiments, the PST-25 test in the greenhouse and with the Mt. Vernon three-year mean data as well as the mean data of all six- field experiments (Tables 5 and 6). The QTL explained 1.94-10.19% of the observed phenotypic variation. QYrdr.wgp-1BL.1 was significantly associated with SSR marker Xgwm131 and flanked by Xwmc611 and Xwmc694 at approximate distances of 3.9 and 2.1 cM, respectively (Fig 4). The nearest SNP marker to QYrdr.wgp-1DS was IWA2268 (Fig 4). This QTL was detected in five of the six field experiments and the two tests with races PST-25 and PST-127 in the greenhouse, and with the means from the Mt. Vernon three-year data and the six field-experiment data. It explained 2.04-27.24% of the phenotypic variance (Tables 5 and 6). QYrdr. wgp-3AL was detected in four of the six field experiments and in the tests with races PST-25 and PST-127 in the greenhouse, and with means of the data from Pullman, Mt. Vernon and both Pullman and Mt. Vernon. This QTL, explaining 1.78-13.85% of the observed phenotypic variation, was closely associated with IWA6834 (Fig 4; Tables 5 and 6).

Effect of the number of QTL in combination
To determine the effects of QTL in various combinations on HTAP resistance, the 94 RILs were classified into genotypic groups based on the presence of markers closely associated with the eight QTL (S1 Table). These genotypes were further grouped into eight groups based on the number of potential QTL for HTAP resistance. Fig 5 shows the differences in the mean IT and rAUDPC values of the eight groups. In general, RILs with more resistance QTL had lower IT and rAUDPC values. This observation is an indication of additive effects of the QTL, with resistance increasing as the number of QTL increases. This is supported by the significant additive effects (P < 0.01) obtained by the QTL analysis.

Interactions between the QTL and environments
The interactions to different environments were determined for each of the HTAP QTL. Significant interactions were detected for QYrdr.wgp-1BL.2 (P = 0.01 for IT and 0.06 for rAUDPC), QYrdr.wgp-1DS (P = 0.0001 for IT and 0.05 for rAUDPC) and QYrdr.wgp-6BL.1 (P = 0.06 for IT and 0.02 for rAUDPC). For the other five QTL, the environmental interactions were not significant for either the IT or rAUDPC data (P = 0.16 to 0.99). The lack of significant interaction with different environments for the majority of the QTL may explain the durability and consistency of the HTAP resistance against the diverse races that are virulent to the seedlings of Druchamp over the years and in the different locations.

Discussion
Druchamp has been resistant to stripe rust since its introduction from France to the US more than 60 years ago. Even after the first Druchamp seedling-virulent race, PST-6, was detected in 1974 and subsequently many virulent races were predominant [9], [16], the cultivar has never become susceptible in the field nurseries. Druchamp was determined to have both race-specific all-stage and race non-specific HTAP resistance [10][11][12][13][14]. In the present study, we mapped three QTL for the all-stage resistance and eight QTL for the HTAP resistance in Druchamp. Compared to the previous QTL mapping studies [4], [34], this study reports the highest number of QTL for resistance to stripe rust in a single wheat cultivar.
QYrdr.wgp-1BL.  [41]. QYr.jirc-1B, a minor QTL derived from 'Fukuho-komugi', with the closed marker Xwmc320, was located in the region near the centromere [35]. QYr.cimmyt-1BL was mapped between SSR markers Xgwm140 and Xgwm259 [36]. Based on the closely linked markers and chromosomal positions, the two QTL for HTAP resistance identified in the present study are likely different from these genes. QPst.jic-1B was considered the same as Yr29, as it mapped to the distal end of 1BL between markers Xgwm818 and Xgwm259, the same position as Yr29 on 1BL [38]. Yr29 is flanked by the SSR loci Xwmc44 and Xgwm259 [37]. These two markers are also flanking QYrdr.wgp-1BL.2 (Fig 4), suggesting that QYrdr.wgp-1BL.2 is likely the same gene as Yr29 or a different allele at that locus. The chromosomal region containing QYrdr.wgp-1BL.1 and QYrdr.wgp-1BL.2 on 1BL appears to be a hot spot for resistance genes against stripe rust and other diseases.
Druchamp was reported to have all-stage resistance gene Yr3a, which was reported on chromosome 1B using monosomic analysis [14]. That gene provides resistance to Pst races CDL-1, CDL-21 and CDL-29, which were later renamed as PST-1, PST-21 and PST-29, respectively [1]. In the present study, the gene for resistance to PST-29 was mapped to 5BL (QYrdr.wgp-5BL.1). The two QTL on 1BL, QYrdr.wgp-1BL.1 and QYrdr.wgp-1BL.2, were only detected in adult-plant stage, and no all-stage resistance QTL was detected on 1B with the five tested races including PST-29. The present study corrected the chromosomal location of the gene in Druchamp for resistance to PST-29 to be in 5BL, instead of 1B.

QYrdr.wgp-1DS
Yr25 was reported on chromosome 1D (Calonnec and Johnson 1998), but its specific chromosomal location is unknown. Yr25 is a major gene for resistance to some races in other countries [44], but largely ineffective to the Pst population in the US [45]. Ren et al. [46] mapped adult plant resistance QTL QYr.caas-1DS in wheat cultivar Naxos on chromosome 1DS flanked by SSR markers Xgwm353 and Xgdm33b, and Xwmc432 is close to the QTL region. In the present study, QYrdr.wgp-1DS.1 is about 5 cM away from Xwmc432 in (Fig 4). Thus, these two QTL could be the same or different alleles in the same chromosomal region.

QYrdr.wgp-3AL
On chromosome 3A, only three stripe rust resistance genes have been reported. YrTr2 is a gene for race-specific all-stage resistance, and its chromosomal arm is unknown [53]. Lillemo et al. [33] identified a QTL (QRYr3A.1) on chromosome 3AS linked to SSR marker Xbarc310 in wheat cultivar Saar. YrQ2 conferring slow-rusting in wheat line Xichang 76-9 was mapped on chromosome 3AS and linked to SSR markers Xwmc11 and Xbarc57 within 12.1 cM [54]. Because QYrdr.wgp-3AL confers HTAP resistance and is on chromosome 3AL, it should be a new gene for stripe rust resistance.

QYrdr.wgp-5DL
QYrdr.wgp-5DL confers an all-stage resistance in Druchamp to old races like PST-35 and PST-45. Chen et al. [53] located YrDa2 to chromosome 5D in US cultivar Daws using monosomic analysis, but its specific chromosomal location is unknown. Because Daws is susceptible to PST-35 and PST-45, while QYrdr.wgp-5DL is effective against the two races, the two genes should be different, and QYrdr.wgp-5DL should be a new gene, since no any other stripe rust resistance QTL have been reported on 5DL [4], [34].

Non-race specificity of HTAP resistance in Druchamp
Druchamp and Michigan Amber were included in a stripe rust monitoring nursery that has been evaluated every year at Pullman and Mt. Vernon in the State of Washington and other locations across the US for more than 40 years, and Druchamp has always shown a high level of resistance. In the present study, the F 8 RIL population was evaluated at the two locations for a period spanning six years from 2005 to 2011, during which 6-21 races were detected each year in these two locations with the presence of races virulent on the seedlings of Druchamp [16], [45], [69]. During the field experiments at the two locations in 2010 and 2011, stripe rust samples were collected from some of the RILs, and all samples were identified as races virulent on the seedlings of Druchamp (data not shown). The virulence data of the Pst populations indicated that the resistant phenotypes observed in the fields were controlled by QTL for the HTAP resistance, but not for the all-stage resistance in Druchamp. This is supported by the findings that none of the three race-specific all-stage resistance QTL was detected in the field experiments and that none of the eight QTL detected in the field experiments contributed to any of the observed seedling resistance in the greenhouse experiments. Most of the QTL for HTAP resistance identified in this study were detected in all field experiments and five of the eight HTAP QTL did not show significant interactions with the environments, indicating the race non-specificity nature of the resistance. The observed variations in individual gene effects in the experiments were mostly due to environmental effects as HTAP resistance is affected by temperature, moisture and disease pressure [1], [4].

Conclusions
This study demonstrates the durability of HTAP resistance, and that the high level of HTAP resistance in Druchamp is controlled by a large number of QTL with various degrees of effectiveness. Among the total of 11 QTL detected in Druchamp, at least three (QYrdr.wgp-5DL for race-specific all-stage resistance and QYrdr.wgp-3AL and QYrdr.wgp-6BL.2 for race non-specific HTAP resistance) are new. The study establishes the presence of Yr29 in Druchamp. The QTL for either all-stage resistance or HTAP resistance are mostly recessive as the distribution of stripe rust responses of the F 8 RILs were skewed more toward susceptibility. Because the HTAP QTL mostly had additive effects, RILs with more QTL displayed higher levels of resistance (Fig 5). The results of the present study and previous studies indicate that high level of HTAP resistance can be controlled by one, few or a large number of QTL [4]. Although the durability of stripe rust resistance in some other wheat cultivars is more dependent on the type of resistance and less dependent on the number of genes [4], [70], the durability of the resistance in Druchamp may be due to both the HTAP type and the large number of QTL. More studies are needed to elucidate the molecular mechanisms of durable type resistance. Practically, QTL for resistance that has been proven to be durable like the HTAP resistance in Druchamp should be valuable for developing wheat cultivars with durable resistance to stripe rust. Molecular markers linked to the QTL can be useful for marker-assisted selection. Multiple markers should be used as most of the closely linked markers for the QTL identified in this study are SNPs.
Supporting Information S1 Dataset. Infection type data of seedling tests in the greenhouse. Infection types produced by Pst races on F 8 lines derived from the Druchamp × Michigan Amber cross tested in seedling stage under low temperature (4-20°C) in the greenhouse. (XLSX) S2 Dataset. Infection type and severity data of adult-plants tested in the greenhouse. Infection type and severity (%) produced by Pst races PST-25 (