Introduction {#Sec1}
============

Wheat (*Triticum aestivum* L.) is one of the staple food crops and serves nearly 35% of the world populations^[@CR1]^. Powdery mildew (PM), caused by *Blumeria graminis* f. sp. *tritici*, is one of the most destructive wheat diseases worldwide^[@CR2],[@CR3]^. PM occurs in wheat fields where high moisture is available^[@CR4],[@CR5]^ and its epidemics can cause a significant grain yield reduction up to 34%^[@CR6]--[@CR8]^. As high demand for wheat yield increase to meet world-growing population, PM is becoming more severe because denser plant canopy of modern cultivars under high nitrogen fertilizer supply favors PM development. Use of resistant cultivars is the most economical environmental friendly approach to reduce the yield loss due to the disease^[@CR9]^. However, rapid development of new races of the pathogen can quickly overcome the host resistance and result in PM outbursts^[@CR10]^.

Wheat PM resistance is conditioned by both race specific and non-race specific resistance genes. To date, 55 resistant genes have been formally designated (*Pm1*-*Pm55*) and mapped to various chromosomes^[@CR11],[@CR12]^. *Pm1a* was the first PM resistance gene identified on the long arm of chromosome 7A from a Canadian wheat cultivar Axminster^[@CR13]--[@CR17]^. Since then, other genes have been identified and tightly linked DNA markers have been reported. Some PM resistance genes, such as *Pm2* on 5DS^[@CR13],[@CR14]^, *Pm3* on 1AS^[@CR14],[@CR18],[@CR19]^, *Pm18* on 7A^[@CR13]^, were identified from bread wheat, whereas others were introgressed into wheat from wheat close relatives, including *Pm4a* on 2AL from *T. monococcum, Pm4b* on 2AL from *T. carthlicum* ^[@CR14]^, *Pm12* on 6BS from *Aegilops speltoides* ^[@CR20]^, *Pm13* on 3B and 3D from *Ae. longissima* ^[@CR21]^, and *Pm25* on 1A from *T. monococcum* ^[@CR22]^. Some *Pm* genes were identified from non-progenitor species of wheat, including *Pm7*, *Pm8*, *Pm17*, *Pm20* and *PmCn17* from *Secale cereal*, *Pm21* from *Haynaldia villosa*, and *Pm40* and *Pm43* from *Thinopyrum intermedium* ^[@CR23]--[@CR30]^. Most of these genes have been used in breeding for PM resistance in wheat. However, new PM pathogen races may easily overcome the resistance conferred by these major genes^[@CR31]^. Therefore, it is very important to continuously explore new resistance genes to diversify resistance sources against rapidly evolving pathogen races.

Molecular markers have been successfully used to tag genes or quantitative trait loci (QTLs) and to estimate their effects^[@CR32]--[@CR35]^. Traditionally, bi-parental mapping populations are used to determine the locations of resistance genes or QTLs in one or two cultivars^[@CR36]--[@CR40]^. However, markers linked to a QTL identified from a specific mapping population may not be useful for marker-assisted breeding in other breeding populations^[@CR41]^. More recently, association mapping (AM) has been used for identification and dissection of disease resistance genes or QTLs in many plant species^[@CR42]^. Unlike bi-parental QTL mapping, AM does not require development of populations, and can quickly assemble a population by collecting a set of diversity germplasm^[@CR43]^. Tightly linked markers to QTLs identified from AM study have high potential to be used for selection of the identified genes or QTLs in breeding.

Population structure displayed a systematic difference in allele frequencies between subpopulations^[@CR44]^. The unequal distribution of alleles among the subpopulations may result in an increase in false association^[@CR42],[@CR45]^. Population structure (Q) and genetic relatedness (Kingship) can be integrated together into a statistical model for association analysis to reduce such false association^[@CR46],[@CR47]^.

In this study, we used the AM approach to study PM resistance in a set of elite breeding lines from U.S. winter wheat breeding programs. The objectives of this study were to 1) determine the genes or QTLs for PM resistance in a panel of US winter wheat germplasm, and 2) identify SNP and SSR markers associated with the QTLs or genes for marker-assisted selection.

Results {#Sec2}
=======

Powdery mildew resistance in the association mapping population {#Sec3}
---------------------------------------------------------------

The PM scores of the 185 wheat accessions in the association mapping population ranged from 0 to 90%. In seedling stage, about 18% accessions were resistant when \<20% severity was classified resistant; whereas in adult plant stage, 43% and 66% of accessions, respectively, were resistant in the two greenhouse experiments and 71% of the tested accessions showed resistance in the field experiment (Supplemental Table [S1](#MOESM1){ref-type="media"}, Fig. [1](#Fig1){ref-type="fig"}). Thirty-three lines including 20 soft wheat lines and 13 hard wheat lines showed immune response in adult plant stage for all three experiments and they are good sources of resistance to PM. About 50 lines showed trace PM at adult stage in one of the three experiments. The results suggest that more tested accessions, especially soft winter wheat, had adult plant resistance than seedling resistance. For adult plant resistance, more resistant accessions were observed in the 2013 field and the 2011 greenhouse experiments than the 2010 greenhouse experiments. The correlation coefficients of PM scores were all significant among the four experiments and was the highest between 2010 and 2011 greenhouse experiments (r = 0.762) and the lowest between 2013 field and 2013 greenhouse (r = 0.466) experiments, suggesting that some adult plant resistance genes expressed in the field experiment might be different from these expressed in seedling stage in the greenhouse experiments.Figure 1Percentage of powdery mildew resistant and susceptible accessions in the association mapping population evaluated at adult plant stage in fall 2010 (2010FA) and spring 2011 (2011SA) greenhouse experiments and field 2013 experiment (2013_Field), and at seedling stage in spring 2013 greenhouse experiment (2013SS).

Population structure {#Sec4}
--------------------

Structure analysis indicated that the population could be divided into three groups (Fig. [2](#Fig2){ref-type="fig"}). Group I (102 accessions) and Group II (23 accessions) are mainly hard winter wheat, whereas Group III (60 accessions) is mainly soft winter wheat. One quarter of accessions in the Group I are hard white winter wheat (HWWW) and others are hard red winter wheat (HRWW). Group II consists of all HRWW with most of accessions having Jagger (13 accessions) in their pedigrees. Principal coordinate analysis (PCoA) and a similarity matrix heat-map (Supplementary Fig. [S1](#MOESM1){ref-type="media"}) confirmed the three groups derived from structure analysis (Fig. [3](#Fig3){ref-type="fig"}).Figure 2Structure analysis devided the population of 185 U.S. wheat accessions into three groups (GI, GII and GIII). Figure 3Principal coordinate analysis (PCoA) separated the population of 185 U.S. wheat accessions into three groups that correspond to the three groups derived from structure analysis.

AMOVA showed that individuals within groups accounted for 89% of the genetic variation, whereas only 11% was explained by the variation among the groups (Table [1](#Tab1){ref-type="table"}).Table 1Analysis of molecular variance on the association mapping population of 185 winter wheat accessions using SNP data.SourcesdfSSMSEst. var.%*P* valueAmong pops21939.733969.86615.82611%0.001Within pops18222579.348124.062124.06289%Total18424519.081139.889100%

Markers significantly associate with PM resistance {#Sec5}
--------------------------------------------------

The optimal model was selected based on the observed *P*-values and the expected *P*-values for each trait (Q--Q plots). The MLM model including population structure (Q) and kingship (K) showed a better fit than GLM, and thus was selected for further association analysis (Supplementary Fig. [S2](#MOESM1){ref-type="media"}). A total of 37 SNPs showed significant associations with PM resistance (Table [2](#Tab2){ref-type="table"}, Fig. [4](#Fig4){ref-type="fig"}) and were mapped to nine regions on eight chromosomes (Fig. [5](#Fig5){ref-type="fig"}). SNPs on three chromosomes (6A, 1B and 1D) were significantly associated with PM resistance in at least two experiments. Other five chromosome regions showed a significant association with PM resistance in only one experiment. On chromosome 6A, three SNPs (wsnp_Ex_c11621_18716254, RFL_Contig5693_646, RAC875_c16962_288) were significant in the 2010 greenhouse experiment and two SNPs (Tdurum_contig9612_80, Excalibur_c23748_1050) were significant in the 2011 greenhouse experiment. The five SNPs were located within 1 cM (Table [2](#Tab2){ref-type="table"}) and mapped in the same LD group (Fig. [6](#Fig6){ref-type="fig"}), and thus they are tightly linked markers and most likely associated with the same gene that showed a significant effect on PM resistance in the two experiments. Four SNPs on chromosome 1B were significant in the 2010 and 2011 greenhouse experiments. One of them was mapped in the chromosome position 43.66 cM and other three were mapped in about 19 cM away on the reference map (Table [2](#Tab2){ref-type="table"}, Fig. [5](#Fig5){ref-type="fig"})^[@CR48]^. LD analysis indicated they were in the same LD (Fig. [6](#Fig6){ref-type="fig"}) and likely the same QTL for PM resistance. Three SNPs that were mapped within a 5 cM region on the chromosome 1D were also significantly associated with PM resistance in the 2010 and 2011 greenhouse experiments (Table [2](#Tab2){ref-type="table"}), and they are more likely linked to the same QTL on the chromosome 1D.Table 2Significant SNP markers associated with wheat powdery mildew resistance evaluated in the greenhouse experiments of fall 2010 (2010FA) and spring 2011 (2011SA) and the field experiment of 2013 (2013_Field) for adult plant resistance, spring 2013 greenhouse for seedling resistance (2013SS).Marker nameChromosomePosition^\*^ (cM)2010FA2011SA2013_Field2013SSResistance alleleSensitive allele*P* value*R* Squared*P* value*R* Squared*P* value*R* Squared*P* value*R* SquaredKukri_c46010_1461A56.394.31E-050.09------------------CTwsnp_Ex_c31983_407096071A56.394.31E-050.09------------------CTRa_c5696_15561A70.12.85E-040.07------------------TCwsnp_Ex_c12101_193602131A70.13.71E-040.06------------------TCBobWhite_c34661_2081A70.14.81E-040.06------------------GTwsnp_Ex_c31525_403027471A70.15.90E-040.06------------------TCKukri_c4099_3211A70.15.90E-040.06------------------CTJagger_c2047_3622A96.63------------1.76E-040.07------GAKu_c57425_4132A96.63------------1.76E-040.07------GACAP8_c4813_2522A97.51------------7.72E-050.07------ACKukri_rep_c73477_8882A97.51-----------1.81E-040.06------TCwsnp_Ex_c2887_53304262A97.51------------1.81E-040.06------TCwsnp_Ex_c2887_53307872A99.28------------2.29E-040.06------TCBS00068178_515A70.3------------------1.91E-040.08TCKukri_c14889_1165A70.3------------------1.91E-040.08TCGENE-3189_3775A70.3------------------2.73E-040.07CTBS00098207_515A70.3------------------2.73E-040.07TCwsnp_Ex_c11621_187162546A42.957.32E-040.06------------------TCTdurum_contig9612_806A43.1------5.12E-040.05------------TCRFL_Contig5693_6466A43.17.32E-040.06------------------AGRAC875_c16962_2886A43.17.32E-040.06------------------CTExcalibur_c23748_10506A43.1------7.92E-040.05------------AGBobWhite_c14258_4341B43.665.67E-040.069.12E-040.05------------ACBS00093945_511B62.585.39E-040.06------------------AGJD_c12243_3601B64.14.47E-040.068.45E-040.05------------CTBobWhite_c26130_4751B64.15.67E-040.069.12E-040.05------------CTBobWhite_c19725_13293B70.09------------------8.87E-050.09GAwsnp_JD_c19725_177325263B70.09------------------9.96E-050.09AGExcalibur_c5298_1713B70.09------------------1.14E-040.08CTExcalibur_c5416_84622-------------------6.18E-040.07ATKu_c71049_118022-------------------6.13E-040.07TCwsnp_Ex_c13496_212431675B40.55------------7.03E-040.05------AGKukri_c18702_1325B41.35------------5.89E-040.05------AGwsnp_Ex_c20988_301076095B41.35------------8.10E-040.05------GABS00077837_511D3.55.67E-040.069.12E-040.05------------AGIACX65861D3.55.67E-040.069.12E-040.05------------TCKukri_c3181_23401D8.715.67E-040.069.12E-040.05------------TGNote: "---" denotes not significant. ^\*^The position of markers on wheat consensus SNP map. Figure 4Manhattan plots for SNPs associated with powdery mildew resistance in wheat. Figure 5Linkage maps to show chromosomal locations of the significant QTLs for powdery mildew resistance from association mapping using 185 winter wheat accessions. Figure 6Linkage Disequilibrium (LD) of some significant SNP markers.

In five other chromosome regions that were significantly associated with PM resistance in a single experiment (Table [2](#Tab2){ref-type="table"}), two SNPs, Kukri_c46010_146 and wsnp_Ex_c31983_40709607, which were mapped at 56.39 cM on chromosome 1A (Table [2](#Tab2){ref-type="table"}), were significant in the 2010 greenhouse experiment. Five other significant SNPs (Ra_c5696_1556, wsnp_Ex_c12101_19360213, BobWhite_c34661_208, wsnp_Ex_c31525_40302747, and Kukri_c4099_321) were mapped at 14 cM away on the same chromosome (Table [2](#Tab2){ref-type="table"}). The two sets of markers showed no LD and more likely associated with two different QTLs for PM resistance (Fig. [6](#Fig6){ref-type="fig"}). Six SNPs were mapped within a 3 cM region on the chromosome 2A and all significant in 2013 field experiment (*R* ^2^ = 0.073, *P* \< 0.000077), suggesting one PM resistance QTL is more likely in the region. Three SNPs mapped within a 1 cM region on 5B were also significant in 2013 field experiment (*R* ^2^ = 0.078, *P* \< 0.0003). Four SNPs on chromosome 5A (*P* \< 0.0003) and five SNPs on chromosome 3B (*P* \< 0.0007) were significant in the 2013 greenhouse experiment. The four SNPs on 5A were mapped together, suggesting that one QTL for seedling resistance may link to these SNPs (Table [2](#Tab2){ref-type="table"}). SNPs on 3B either were mapped together or shared the same LD, suggesting they more likely link to a single gene for seedling resistance on the chromosome 3B (Table [2](#Tab2){ref-type="table"}, Fig. [6](#Fig6){ref-type="fig"}).

Among 457 SSR markers from all 21 chromosomes screened, three were significant for PM resistance. Xscm0009 on chromosome 1A was significantly associated with PM resistance in 2010 greenhouse experiments (*R* ^2^ = 0.08, *P* \< 0.0007, Table [3](#Tab3){ref-type="table"}). Two other markers, Xcfd9-2 on 3D and Xcfd95 on 6D, were significant in the 2013 field experiment.Table 3Significant SSR markers associated with wheat powdery mildew resistance evaluated in the greenhouse experiments of fall 2010 (2010FA) and spring 2011 (2011SA) and the field experiment of 2013 (2013_Field) for adult plant resistance, spring 2013 greenhouse for seedling resistance (2013SS).Test MethodsMarker nameChromosome2010FA2011SA2013_Field2013SS*P* value*R* Squared*P* value*R* Squared*P* value*R* Squared*P* value*R* SquaredMLMXscm00091 A7.09E-040.08------------------Xcfd9-23D------------5.92E-040.09------Xcfd956D-----------3.78E-040.07------Note: "--" denotes not significant.

The adult plant PM resistance level of each line is highly correlated with number of QTLs in the lines (r = 0.6083). These QTLs showed obvious additive effects. When a line carried four or more QTLs, it usually showed immune response to PM infection.

Discussion {#Sec6}
==========

Population structure may result in false associations between traits and markers if it is not properly treated during analysis^[@CR49]^. In this study, both structure analysis and PCoA stratified the AM panel of 185 accessions into three groups and AMOVA also indicated significant population differentiation (*P* \< 0.001), demonstrating the presence of obvious structure in the population. Wheat cultivars from NRPN were mostly clustered to Group I and wheat cultivars from RGON and SRPN were mostly clustered to Group I and Group II, whereas accessions in USSRWWN and UESRWWN showed higher PM resistance and clustered into Group III. The result indicated that PM resistance was partly influenced by the geographic distribution.

A significantly higher level of variation was detected within a population (89%) than among populations (11%) by AMOVA, indicating that a high level of genetic diversity, but not genetic divergence, was observed in US winter wheat, therefore, breeding selection has played a significant role in maintaining genetic diversity in the breeding populations.

Using the AM panel in this study, we identified SNP markers linked to nine QTLs for PM resistance and SSR markers linked to three QTLs for PM resistance on the chromosomes 1A, 2A, 5A, 6A, 1B, 3B, 5B, 1D, 3D and 6D. Among them, SNPs closely linked to the three genes on the chromosomes 6A, 1B, 1D were significant in adult plant stage in two experiments, suggesting these genes may have stable adult plant resistance to multiple races presented in both environments.

When significant QTLs identified in the current study were compared with the QTLs reported previously, we found that several QTLs have been mapped to the similar positions where PM resistance genes were reported previously. The seven SNPs that were significantly linked to PM resistance QTL in the 2010 field experiment were mapped to two locations of the chromosome 1A at 15 cM apart, and LD analysis suggested that they were linked to two different genes. *Pm3* has been previously reported on the chromosome 1A^[@CR50]--[@CR54]^; *Pm 17* on chromosome 1RS.1AL was mapped about 1.5 cM away from the IAG95-CA (RFLP)/CT355 (AFLP) markers^[@CR55]^; and *Pm25* on chromosome 1A was mapped at about 12.8 cM from a RAPD marker OPA04950^[@CR22]^. Five SNP-associated QTLs on the chromosome 1A are more than 40 cM away from *Pm3*, therefore neither of them are *Pm3*. However, this QTL could be the same or a closely linked gene to *QPm.caas-1A* between SSR marker Xbarc148 and Xwmc550^[@CR56]^ according to the marker sequence position in the W7984 reference map (Supplemental Table [S2](#MOESM1){ref-type="media"}). The SSR marker, Xscm0009 was associated with PM resistance in the greenhouse experiments of 2010. The banding pattern of Xscm0009 showed that the marker is in 1A/1 R translocation, which suggests the linked gene is most likely *Pm17*. For other two SNPs on chromosome 1A, known linked genes can not be found based on available information, and they may link to a novel gene. However further research is needed to determine its identity.

*Pm4a*, *Pm4b*, *Pm4c* (or *Pm23*), *Pm4d*, *PmDR147*, *PmPS5A* and *PmLK906* were identified on the chromosome 2A in previous reports^[@CR14],[@CR57]--[@CR60],[@CR61],[@CR62]^. Among them, *Pm4a* was about 1.5 cM away from the marker Xbcd292^[@CR14]^ and *QPm.inra-2A* is linked to the SSR marker *Xgwm275* that is also in the vicinity of *Pm4a* ^[@CR63]^. In this study, the QTL on 2 A was mapped to a chromosome region near *QPm.inra-2A* according to the marker sequence position in the W7984 reference map (Supplemental Table [S2](#MOESM1){ref-type="media"}) therefore it may be the same QTL as *QPm.inra-2A*.

*Pm2026* was mapped to the distal portion of chromosome 5 AL and flanked by SSR markers Xcfd39 and Xgwm126^[@CR64]^. The QTLs *QPm.sfr-5A.1*, *QPm.ttu-5A*, *QPm.sfr-5A.2*, *QPm.sfr-5A.3* and *QPm.nuls-5A* were also reported on the chromosome 5A^[@CR65]--[@CR67]^. The QTLs *QPm.sfr-5A.1*, *QPm.sfr-5A.2*, *QPm.sfr-5A.3* were linked to RFLP markers^[@CR66]^. *QPm.ttu-5A* and *QPm.nuls-5A* were mapped in the marker intervals Xgwm186-Xgwm415 and Xgwm617b-Xwmc327, respectively^[@CR65],[@CR67]^. In the current study, the 5A QTL is likely different from the QTLs *QPm.ttu-5A*, *QPm.nuls-5A* and *Pm2026* based on its marker position in the W7984 reference (Supplemental Table [S2](#MOESM1){ref-type="media"}). Whether the 5A QTL is the one of the previously reported genes remains to be determined due to lack of common markers among the QTLs reported in different studies.

One QTL was identified on chromosome 3B in this study. This QTL is likely *PmHNK*, a gene was previously mapped at 3.8 cM away from the SSR marker Xwmc291^[@CR68]^ according to the W7984 reference. The map locations of SNPs, Excalibur_c5416_846 and Ku_c71049_1180, are not available, but a strong LD between the two SNPs and the markers mapped on the chromosome 3B located the two SNPs on the same chromosome (Fig. [6](#Fig6){ref-type="fig"}).

*Pm30, Pm36, PmAS846* and *MIVE29* were all mapped on chromosome 5B in previous studies. Liu *et al*.^[@CR69]^ mapped *Pm30* on the chromosome 5B about 6 cM from the SSR marker Xgwm159/460. Three significant SNPs (wsnp_Ex_c13496_21243167, Kukri_c18702_132, wsnp_Ex_c20988_30107609) mapped on the chromosome 5B in the current study are close to *Pm30* based on the W7984 reference map (Supplemental Table [S2](#MOESM1){ref-type="media"}), therefore, they are likely the same gene.

*Pm10*, *Pm22*, *Pm24a* and *Pm24b* were previously reported on chromosomes 1D^[@CR37],[@CR70]--[@CR72]^. Tightly linked markers were identified for *Pm24a* ^[@CR37]^ and *Pm24b* ^[@CR72]^. Other QTLs have also been reported in this chromosome including *QPm.inra-1D.1* and *QPm.sfr-1D* ^[@CR63],[@CR66]^. The markers flanking *Pm24a* were Xgwm789/Xgwm603 and Xbarc229 at 2.4 and 3.6 cM to the gene, respectively. *Pm24b* was flanked by Xgwm337 and Xbarc229 at the genetic distances of 3.7 and 1.0 cM, respectively. Here, we identified three significant SNPs linked to a gene on chromosome 1D (Table [2](#Tab2){ref-type="table"}, Fig. [5](#Fig5){ref-type="fig"}), which is 2.2 cM away from SSR marker Xbarc229 based on W7984 reference map, and thus the QTL identified in this study is more likely *Pm24*.

Among these nine PM resistance QTLs/genes identified by SNP markers in our study, some are located in the chromosome locations where PM resistance genes have not been reported before, and thus they are likely novel genes for PM resistance (Table [2](#Tab2){ref-type="table"}). Resistance genes *Pm21, Pm31, PmY39-2* and *MIRE* were mapped on chromosome 6A^[@CR73]--[@CR76]^. However, the QTL on 6A identified in the current study does not link to any of the genes and is more likely a new QTL for PM resistance.

*Pm8*, *Pm32* and *Pm39* were previously reported on chromosomes 1B^[@CR28],[@CR67],[@CR77]^. However, *Pm8* and *Pm32* have not been mapped to date. *Pm39* ^[@CR67]^ was closely linked to SSR marker Xwmc719, but it is about 56.9 cM away from the 1B QTL identified in this study (Table [2](#Tab2){ref-type="table"}, Fig. [5](#Fig5){ref-type="fig"}) according to the W7984 reference map. Several other QTLs have been reported in 1B including *QPm.sfr-1B*, *QPm.ttu-1B*, *QPm.vt-1B*, *QPm.vt-1BL*, *QPm.vt-1B*, *QPm.osu-1B* ^[@CR65],[@CR66],[@CR78]--[@CR80]^, but we can not determine the relationship between the newly identified QTL from this study and previous reported ones due to lack of common markers between these studies.

The SSR markers, Xcfd9-2 on chromosome 3D and Xcfd95 on chromosome 6D, also showed significant associations with PM resistance. *Pm45* was mapped on chromosome 6D flanked by Xcfd80, Xmag6139 and Xmag6140^[@CR81]^. The QTL detected in this study can be mapped near *Pm45* based on the W7984 reference sequence. The QTL detected on 3D is close to *QPm.inra-3D* that was flanked by Xcfd152 and Xgwm707^[@CR63]^, hence, they may be the same QTL.

In this study, soft wheat carries at least two resistance QTLs that were identified by SNP markers, and showed a higher level of resistance than hard wheat (Supplemental Table [1](#Tab1){ref-type="table"}). For adult plant, a high level correlation between number of QTLs and PM resistance was observed (r = 0.6083). The PM resistance QTLs showed obvious additive effects. The result indicated that four or more QTLs together produce immune response in adult plant. Therefore, pyramiding four or more of the QTLs using marker-assisted selection can obtain PM immune wheat cultivars.

Our study demonstrates that most of US winter wheat germplasm, especially soft winter wheat, have high levels of PM resistance and carry several QTLs for adult plant resistance to PM. Association mapping effectively identified these genomic regions associated with PM resistance and associated markers linked to these QTLs. They showed obvious additive effect on adult plant PM resistance. The accessions carrying multiple resistance QTLs could be excellent sources of PM resistance for US winter wheat breeding because they are either elite breeding lines or locally adapted cultivars. Future work will be to develop bi-parental populations of the accessions to validate the resistance loci identified in this study and develop user-friendly markers that can be used to accelerate the incorporation of these resistance QTLs into new cultivars.

Materials and Methods {#Sec7}
=====================

Plant materials {#Sec8}
---------------

A set of 185 elite breeding lines and cultivars were selected from the 2008 U.S Southern (SRPN) and Northern (NRPN) Hard Winter Wheat Regional Performance Nurseries, Regional Germplasm Observation Nurseries (RGON), U.S. Uniform Eastern Soft Red Winter Wheat Nursery (UESRWWN), Uniform Southern Soft Red Winter Wheat Nursery (USSRWWN), and elite breeding lines from Oklahoma State University by removal of sibling lines (Supplemental Table [S1](#MOESM1){ref-type="media"}). Among these accessions, 130 are hard winter wheat and 55 are soft winter wheat.

DNA extraction and marker analysis {#Sec9}
----------------------------------

Leaf tissue at the two-leaf stage was collected into 1.1-ml deep-well plates and dried for 2 d in a freeze-dryer (Thermo Fisher Scientific Inc., Waltham, MA, USA) for DNA isolation. The plates containing dried tissue and a 3.2-mm stainless steel bead in each well were shaken at 25 times per sec for 5 min in a Mixer Mill (Retsch GmbH, Germany). Genomic DNA was extracted using the cetyl trimethyl ammonium bromide method and SSR markers were analyzed using an M13-tailed primer as described by Li *et al*.^[@CR82]^. All PCR products were separated on an ABI PRISM 3730 DNA Analyzer (Applied Biosystems, Foster City, CA, USA), and marker data were scored using GeneMarker version 1.6 (Soft Genetics LLC, State College, PA, USA) and manually checked twice for accuracy.

Wheat accessions were genotyped using the wheat 90 K SNP assay developed by Illumina Inc. (San Diego, CA, USA). The assay was designed under the protocols of the International Wheat SNP Consortium^[@CR83]^ and conducted at the USDA Small Grains Genotyping Laboratory in Fargo, ND. SNP genotypes were called using GenomeStudio v2011.1 software (Illumina Inc.). SNPs with minor allele frequency less than 5% or with missing data more than 15% were removed. A total of 21,600 SNPs were scored and used for association analysis.

Disease evaluation {#Sec10}
------------------

PM was evaluated in both greenhouse and field experiments at Kansas State University, Manhattan, KS from 2010 to 2013. In the greenhouse experiments, five plants were transplanted into 12.5 × 12.5 cm Tara pots (Hummert International, Topeka, KS) after seven weeks of vernalization at 6 °C. The greenhouse temperatures were set at 22 ± 5 °C during a day with supplemental light of 12 h, and at 17 ± 2 °C during a night. The cultivar, Wesley, was used as the susceptible control in all trials. Field experiments were carried out in Rocky Ford Wheat Disease Nursery in Manhattan, KS. Thirty seeds per accession were sown in a 1.3-m row using a randomized complete block design with two replicates. Naturally occurred *B. graminis* f. sp. *tritici* was used as inoculum for both field and greenhouse experiments.

PM severity was visually estimated as overall percentage of infected leaf area two weeks after anthesis when disease reached maximum levels following Chen and Xu^[@CR10]^. In each experiment, two replications were evaluated for each accession and a mean value of two replications was used for association analysis. Resistance classification followed Liu *et al*. ^[@CR31]^, with some modification. In brief, plants without PM symptom (0) or low PM coverage (≤20%) were rated as resistant, whereas plants with \>20% PM coverage were rated as susceptible.

Association analysis {#Sec11}
--------------------

The population structure (Q matrix) was analyzed using the program STRUCTURE version 2.2^[@CR84]^, with a burn-in length of 10,000 and a total of 10,000 Markov chain Monte Carlo iterations for each *k*. Ten independent runs were carried out for each *k* value. The maximum likelihood of each *k* value, the variance among 10 runs, and the pedigree information of each line were weighted to determine the optimal number of groups. The relative kinship (K) matrix was calculated using SPAGeDi^[@CR85],[@CR86]^. Components of genetic variances among and within groups were estimated by analysis of molecular variance (AMOVA) with 1000 permutations using GenAlEx 6.501^[@CR87]^.

Association analysis was carried out using both generalized linear model (GLM) and mixed linear model (MLM). GLM includes Q matrix for fixed effects, whereas MLM includes both a Q matrix for fixed effects and a kinship matrix for random effects. The observed *P* values and expected *P* values for each trait (Q-Q plot) were used for model comparison to select the best model^[@CR88]--[@CR90]^. Association analysis between SSR and PM severity was conducted using TASSEL 2.1^[@CR47],[@CR91]--[@CR93]^. Associations between SNP markers and PM severity were determined using the Genome Association and Prediction Integrated Tool (GAPIT)^[@CR94]^, an R package for genome wide association study (GWAS) and genome prediction (<http://www.r-project.org>). A threshold of *p* \< 0.001 was set up to claim significant associations between SSR or SNP markers and the traits. The genetic positions (cM) of SNP markers on chromosomes were determined based on the 2015 wheat consensus map^[@CR48]^. The marker-trait associations were cross-referenced against all reported QTLs in the literature and the GrainGenes database (<https://wheat.pw.usda.gov/GG3/>)^[@CR95]^. Sequences that harbored significant SNPs were further blasted against the W7984 reference sequences to estimate their putative chromosome positions.

To estimate linkage between significant SNPs, Haploview 4.2 (<http://www.broad.mit.edu/mpg/haploview/>) was used to calculate the linkage disequilibrium (LD) among all significant markers^[@CR96]^. Markers in close vicinity with strong LDs were considered to represent the same gene.

Electronic supplementary material
=================================

 {#Sec12}

Supplementary info Supplementary table S1

**Electronic supplementary material**

**Supplementary information** accompanies this paper at doi:10.1038/s41598-017-11230-z

**Publisher\'s note:** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

This is contribution number 17--J from the Kansas Agricultural Experiment Station. This study was funded by the National Research Initiative Competitive Grants 2011-68002-30029 and 2017-67007-25939, 2017-67007-25929 from the USDA National Institute of Food and Agriculture. Mention of trade names or commercial products in this publication is solely for the purpose of providing specific information and does not imply recommendation or endorsement by the US Department of Agriculture. USDA is an equal opportunity provider and employer.

G.B. conceived and designed this study, N.L. and M.L. performed the experiments, N.L., G.B., M.L., X.X. and W.Z. analyzed the data, N.L. and G.B. wrote the paper. All authors revised the manuscript and approved the final version of the manuscript.

Competing Interests {#FPar1}
===================

The authors declare that they have no competing interests.
