Direct hybridization facilitates the simultaneous identification and introgression of QTL for adult plant resistance to the Ug99 stem rust lineage from tetraploid Khorasan wheat to common wheat

Direct hybridization facilitates the simultaneous identification and introgression of QTL for adult plant resistance to the Ug99 stem rust lineage from tetraploid Khorasan wheat to common wheat | Research Square window.SnipcartSettings = { analytics: { enabled: false } }; (function() { var accessVector = localStorage.getItem('access_vector') || ''; window.dataLayer = window.dataLayer || []; if (accessVector) { window.dataLayer.push({ user: { profile: { profileInfo: { snid: accessVector } } } }); } })(); (function(w,d,s,l,i){w[l]=w[l]||[];w[l].push({'gtm.start':new Date().getTime(),event:'gtm.js'});var f=d.getElementsByTagName(s)[0],j=d.createElement(s),dl=l!='dataLayer'?'&l='+l:'';j.async=true;j.src='https://www.googletagmanager.com/gtm.js?id='+i+dl;f.parentNode.insertBefore(j,f);})(window,document,'script','dataLayer','GTM-K279D39R'); Browse Preprints In Review Journals COVID-19 Preprints AJE Video Bytes Research Tools Research Promotion AJE Professional Editing AJE Rubriq About Preprint Platform In Review Editorial Policies Our Team Advisory Board Help Center Sign In Submit a Preprint Cite Share Download PDF Research Article Direct hybridization facilitates the simultaneous identification and introgression of QTL for adult plant resistance to the Ug99 stem rust lineage from tetraploid Khorasan wheat to common wheat Max Fraser, Emily Conley, Zennah Kosgey, Ashenafi Gemechu Degete, and 2 more This is a preprint; it has not been peer reviewed by a journal. https://doi.org/ 10.21203/rs.3.rs-2958205/v1 This work is licensed under a CC BY 4.0 License Status: Under Revision Version 1 posted 4 You are reading this latest preprint version Abstract The Puccinia graminis f. sp. tritici ( Pgt ) Ug99 race group presents a major challenge to global wheat production. Satisfying current and future demands hinges on the discovery of new sources of resistance. It is imperative that the durability and diversity of Ug99 resistance is improved by identifying and deploying novel resistance genes. Progenitor species and wild relatives of common wheat ( Triticum aestivum ) have proven to be rich sources of genetic diversity. The Khorasan wheat ( Triticum turgidum ssp. turanicum ) accession CItr 11390 displays adult plant resistance (APR) to Ug99 races. 121 BC 1 F 5 -derived recombinant inbred lines were developed from a cross between CItr 11390 and MN07098-6 to map and introgress resistance loci from CItr 11390. The population was evaluated in Kenya and Ethiopia in 2019 and 2020. Two APR QTL from CItr 11390 were detected in multiple environments. QSr.umn-2A is believed to be the APR gene Sr63 on chromosome 2AL. QSr.umn-6BL was identified on 6BL upstream from Sr11 . The distance from Sr11 and lack of APR QTL reported on 6BL suggest QSr.umn-6BL is a novel locus. Additional QTL were mapped to chromosomes 1AS, 3AL, 3BL, 5AL, and 6BS in single environments. The population segregates for TKTTF seedling resistance conferred by Sr7a and a novel locus, QSr.umn-5A.1 . The population consists of the first hexaploid wheat lines to pyramid Sr7a , Sr57 / Lr34 / Yr18 , Sr63 , and QSr.umn-6BL . This study is the first report of Pgt resistance QTL from Khorasan wheat, and it demonstrates the feasibility of simultaneously identifying and transferring resistance QTL from tetraploid to hexaploid wheat. Khorasan wheat Triticum turgidum ssp. turanicum Stem rust Ug99 Introgression Figures Figure 1 Figure 2 Figure 3 Figure 4 Figure 5 Key message Multiple quantitative trait loci for adult plant resistance to Pgt races endemic to Kenya and Ethiopia were simultaneously identified and introgressed from tetraploid wheat accession CItr 11390 into common wheat. Introduction Worldwide, common wheat ( Triticum aestivum L.) is the most widely grown crop with approximately 220 million hectares harvested annually. Global production of ~ 780 million tons of wheat provides 20% of the world’s calories and protein (FAO 2022 ; Braun et al. 2010 ). In order to sustain this demand, wheat yields must increase 60% by 2050 (Godfray et al. 2010 ). Unfortunately, trends in wheat yields are inconsistent across major production regions and increasingly threatened by numerous biotic and abiotic factors (Grassini et al. 2013 ). Among these biotic factors is wheat stem rust caused by the fungus Puccinia graminis f. sp. tritici Eriks. and E. Henn. ( Pgt ). Capable of destroying entire fields and reducing regional grain yields by up to 70% during epidemics, stem rust has been among the most damaging and feared diseases of wheat throughout history. Historical and archeological records of Pgt infection and stem rust epidemics extend across all major growing regions and date as far back as 3,300 BC (Kislev 1982 ; Roelfs 1985 ). However, outbreaks have greatly diminished since the early 20th century. Widespread deployment of genetic resistance to Pgt , disease control programs such as the Barberry Eradication Program in the United States, and the development of chemical fungicides effectively contained stem rust outbreaks to manageable levels throughout much of the 20th century (Singh et al. 2011 ). During this period, stem rust disease resistance was largely qualitative in nature. Race-specific, all-stage resistance (ASR), also referred to as seedling resistance, was the cornerstone of stem rust resistance breeding efforts. The release of the hard red spring wheat variety ‘Thatcher’ by the University of Minnesota in 1934 was in direct response to a series of epidemics that afflicted the growing region in the preceding decades (Hayes et al. 1936 ). The rise of ‘Thatcher’ as the predominant variety grown throughout the Northern Great Plains of the United States and Canada into the 1960’s was driven by a stem rust resistance package built around the ASR gene Sr12 (Knott 2000 ). Beyond North America, stem rust was effectively controlled through the introgression and deployment of another ASR gene, Sr31 , from the rye ( Secale cereale L.) variety ‘Petkus’ beginning in the mid-1960s (McIntosh et al. 1995 ). Genetic control of stem rust using ASR genes such as Sr12 and Sr31 remained effective throughout much of the mid to late 20th century. Subsequently, breeding efforts and funding for fundamental and applied stem rust research were shifted towards the control of other pathogens (Singh et al. 2011 ). Stem rust reemerged as a central topic of wheat research and breeding following the discovery of Pgt isolate Ug99 in Uganda in 1998 (Pretorius et al. 2000 ) and subsequent wheat stem rust epidemics in Kenya in 2004 and 2005. Later given the race designation TTKSK (Roelfs and Martens 1988 ), it was virulent to Sr31 and several other resistance genes from common wheat and related species. Since the initial discovery of race TTKSK, the pathogen has undergone rapid genetic mutation and moved as far as South Africa, Yemen, and Iran (Terefe et al. 2018; Pretorius et al. 2012 ). Further migration of the Ug99 lineage of Pgt poses a major threat to global wheat production, as a study found that 85–90% of wheat varieties were susceptible to Ug99 (Singh et al. 2015 ). The Ug99 lineage of stem rust is not the only threat to wheat production within Sub-Saharan Africa and beyond. Recent epidemics in Kenya and Ethiopia are attributed to races TKTTF, TTRTF, and JRCQC (Letta et al. 2013 ; Olivera et al. 2015 ; Patpour et al. 2020 ). The movement of these races outside the regions of their initial detection and the outbreaks they have caused highlight the need for constant pathogen monitoring, discovery of novel wheat stem rust resistance genes, and an emphasis on breeding for adult plant resistance. Adult plant resistance (APR) is characterized by a reduction in both urediniospore size and production, and it is generally race nonspecific unlike like ASR (Stuthman et al. 2007 ). APR is believed to be a more durable form of resistance, as it is less prone to the boom and bust cycle associated with the successive deployment of individual ASR genes (Singh et al. 2011 ). Unfortunately, only six of the 63 named stem rust resistance genes confer APR: Sr2 (McFadden 1930 ), Sr55 / Lr67 / Yr46 (Herrera-Foessel et al. 2014 ), Sr56 (Bansal et al. 2014 ), Sr57 / Lr34 / Yr18 (Singh and Huerta-Espino 2003 ), Sr58 / Lr46 / Yr29 (Martínez et al. 2001 ), and Sr63 (Mago et al. 2022 ). Singh et al. ( 2011 ) hypothesized that pyramiding four to five APR genes could provide near immunity to Pgt . However, the expression and effectiveness of resistance genes is dependent on genetic background and environmental conditions (Kolmer 1996 ). Therefore, discovery and characterization of novel genes are essential to optimize gene pyramids and overcome the ability of Pgt to defeat deployed resistance genes. Mapping and dissection of stem rust resistance genes has not been limited to hexaploid common wheat (AABBDD, 2n = 6x = 42). The progenitor species and wild relatives of common wheat have proven to be rich sources of diversity for numerous traits including stem rust resistance (Rowland and Kerber 1974 ). In total, stem rust resistance genes have been identified in 14 relatives of hexaploid wheat. The resistance genes identified in these relatives account for a majority of the genes that remain effective against the races of the Ug99 lineage (Singh et al. 2015 ). However, most of the introgressions are not used by breeding programs (Singh et al. 2015 ; Yu et al. 2014 ). Introgressions from species in the secondary and tertiary genes pools including T. timopheevii Zhuk. (AAGG, 2n = 4x = 28), Secale cereale L., and members of Thinopyrum and Aegilops are carried on Robertsonian translocations, which often harbor unfavorable alleles at genes influencing agronomic traits (Singh et al. 2011 ). As opposed to members of the secondary and tertiary gene pools, the diploid and tetraploid progenitors of common wheat carry homologous chromosomes to those of common wheat capable of undergoing pairing and recombination, albeit with varying difficulty (Gill and Raupp 1987 ). In the case of T. monococcum (AA, 2n = 2x = 14), gene transfer is most often conducted through a bridge cross with durum wheat ( T. turgidum , AABB, 2n = 4x = 28) (Sharma and Gill 1983 ). Direct hybridization with T. monococcum is complicated by the high rates of male and female sterility and lethality in F 1 hybrids (The 1973 ). However, direct introgression from T. monococcum is possible through extensive backcrossing as demonstrated by the transfer of Sr21 and Sr35 into hexaploid germplasm (The 1973 ; Zhang 2010). Introgressions from Aegilops tauschii (DD, 2n = 2x = 14) have historically been facilitated by first producing synthetic allohexaploid wheat through chromosome doubling of durum x Ae. tauschii F 1 hybrids. Direct crossing between Ae. tauschii and common wheat is a faster approach than developing synthetic allohexaploid lines despite the necessity of performing embryo rescue, and it also avoids the disruption of adaptive allele combinations that have evolved in the A and B genomes (Gill and Raupp 1987 ). Olson et al. ( 2013a , b ) demonstrated the feasibility of direct crossing with Ae. tauschii by simultaneously mapping and introgressing stem rust resistance loci SrTA1662 , Sr10171 , and Sr10187 in populations developed after two rounds of backcrossing to T. aestivum recurrent parents. Direct introgression of stem rust resistance genes from emmer and durum wheat (AABB, 2n = 4x = 28) is well documented and does not require specialized techniques (Lanning et al 2008 ; Laugerotte et al. 2022 ). Examples include Sr2 introgressed from T. turgidum ssp. dicoccum variety ‘Yaroslav’ (McFadden 1930 ), Sr11 from T. turgidum ssp. durum variety ‘Gaza’ (Green et al. 1960 ), Sr12 from T. turgidum ssp. durum variety ‘Iumillo’ (Knott 2000 ), and Sr13 and Sr14 from T. turgidum ssp. dicoccum variety ‘Khapli’ (Knott 1962 ). Similar to the approach described by Olson et al., ( 2013a , b ), QTL influencing yield component traits and Yr53 , a gene for stripe rust (caused by P. striiformis f. sp. tritici ) resistance, were simultaneously mapped and transferred into common wheat from tetraploid wheat accessions (Kalous et al. 2015 ; Xu et al. 2013 ). The primary objective of this study was to identify the underlying loci conferring APR to Ug99 lineage races of Pgt in the Khorasan wheat accession ( Triticum turgidum ssp. turanicum ) CItr 11390 (‘Sun Ray’). In an effort to simultaneously transfer these loci to hexaploid wheat lines while reducing the effects of poor agronomic alleles from CItr 11390 and recovering euploid progeny, a BC 1 F 5 -derived recombinant inbred line (RIL) population was developed. Materials and methods Plant materials A mapping population of 121 BC 1 F 5 -derived RILs was developed from an interspecific cross between tetraploid accession CItr 11390 and the hexaploid recurrent parent MN07098-6. CItr 11390 (‘Sun Ray’) is a Khorasan wheat ( Triticum turgidum ssp. turanicum ) accession of unknown improvement status within the USDA National Small Grains Collection. MN07098-6 (SD3696/MN97803-3BS) is a hard red spring wheat breeding line susceptible to the Ug99 lineage from the University of Minnesota. The population was developed from a single BC 1 F 1 plant. The BC 1 F 2 progeny were advanced in the greenhouse using the single seed descent method, and a seed increase was conducted at the BC 1 F 5 generation to obtain enough seed for field experiments. Seedling stage stem rust phenotyping The RILs and parents were evaluated for seedling infection types to Pgt races TTKSK (isolate 04KEN156/04; Ug99), TTKTT (isolate 14KEN58-1; Ug99 with Sr24 and SrTmp virulence), TKTTF (isolate 13ETH18-1), and TTRTF (isolate 14GEO189-1) at the United States Department of Agriculture – Agricultural Research Service (USDA-ARS) Cereal Disease Lab in Saint Paul, Minnesota. Inoculation and phenotyping procedures were conducted as described by Rouse et al. ( 2011 ). All lines were evaluated in four replications. Infection types (ITs) were recorded on a 0–4 scale described by Stakman et al. ( 1962 ). IT scores were converted to a linearized scale for QTL mapping according to Gao et al. ( 2016 ). Field-based stem rust evaluation The BC 1 F 5:6 and BC 1 F 5:7 populations, as well as the parents, were evaluated for field response to African Pgt races during the off-seasons (January – May) of 2019 and 2020 in Njoro, Kenya (hereafter referred to as Ken19 and Ken20, respectively) and Debre Zeit, Ethiopia (hereafter referred to as Eth19 and Eth20, respectively). In all environments, the lines were grown in a randomized complete block design with two replications. In Kenya, each line was grown in two 0.7 m rows spaced 0.3 m apart. The susceptible check CItr 6255 (‘Red Bobs’) was planted after every fifty entries. Twin spreader rows consisting of susceptible varieties ‘Cacucke’, ‘Eagle 10’, ‘Robin’, and six CIMMYT experimental lines were grown alongside and between each plot. Additionally, a border of spreader rows encompassed the nursery. Starting at the booting growth stage through heading (Zadok’s growth stage 37–60), spreader rows were inoculated with a bulk inoculum of locally collected urediniospores. Pgt races TTKSK, TTKST, TTKTK, and TTKTT were among the races present in the Kenya nurseries (Mago et al. 2022 ). In Ethiopia, lines were grown in 1.0 m double rows flanked by spreader rows consisting of susceptible varieties ‘Arendeto’, ‘Digalu’, ‘Local Red’, ‘Morocco’, and ‘PBW343’ (Kosgey et al. 2021 ). Spreader rows were inoculated at Zadok’s growth stages 37–60 using bulk of locally collected urediniospores. Stem rust severity was assessed once spreaders and the check line CItr 6255 reached 50% severity in both Kenya and Ethiopia. Stem rust severity was scored on a whole plot basis using the modified Cobb scale of 0 to 100%, where 0 is equal to complete immunity or absence of infection symptoms and 100% is equal to complete susceptibility (Peterson et al., 1948 ). Infection response (IR) was recorded based on pustule size and prevalence of chlorosis and necrosis. Four classifications were assigned for IR: S, susceptible; MS, moderately susceptible; MR, moderately resistant; and R, resistant (Roelfs et al. 1992 ). Lines with heterogeneous infection responses were denoted using a hyphen (e.g. MS-MR) with the most prevalent IR category listed first. Linearized infection response (LIR) and coefficient of infection (CI) were calculated from the field observations using a custom R script (Gao et al. 2016 ). LIR is infection response scaled from 0 to 1.0. CI was then calculated as the product of LIR and severity (Stubbs et al. 1986 ). Two observations were recorded in both years in Kenya, while three observations were made in Ethiopia. However, only the CI value from the final observation in each environment was considered during data analysis and QTL mapping. Genotyping Single nucleotide polymorphism markers (SNPs) were called using genotyping-by-sequencing (GBS; Elshire et al. 2011 ). DNA was extracted from juvenile leaves collected from the parents and BC 1 F 5:6 RILs using BioSprint 96 DNA Plant Kits (QIAGEN, Hilden, Germany). Extracted DNA was quantified using Quant-iT PicoGreen dsDNA assay kits (Invitrogen, Waltham, Massachusetts, USA) and normalized to 20ng/µL. GBS libraries were constructed using a two-enzyme approach as described by Poland et al. ( 2012 ) with two modifications: 1) two unique barcodes were ligated to each sample in an effort to reduce sequencing bias of certain barcodes and 2) barcode and common adapter concentrations were increased to 0.1 µM and 50 µM, respectively. Two 64-plex libraries were generated including four biological replicates of parent MN07098-6 and three biological replicates of parent CItr 11390. Each library was sequenced on one Illumina HiSeq 2000 lane to produce single-end 100 bp reads. Marker calling, filtering, imputation, and distortion Raw sequencing reads were demultiplexed using Sabre ( https://github.com/najoshi/sabre ). Adapters and low quality bases (Q < 30) were trimmed using Cutadapt (Martin 2011 ). Trimmed reads were aligned to the T. aestivum cv. Chinese Spring reference assembly, RefSeq v2.1 (Zhu et al. 2021 ) using the BWA-MEM algorithm with default parameters (Li 2013). Aligned reads were filtered (MAPQ < 40) and sorted with SAMtools, and SNP calling was done with BCFtools (Li 2011 ; Li et al. 2009 ). Genotypes within a site were set to missing data based on read depth (DP < 3). Sites with ≥ 20% missing data and QUAL < 20 were then discarded as well as monomorphic and multiallelic sites. Sites were further filtered by removing those with a minor allele frequency (MAF) 10%. Missing marker data were imputed in TASSEL 5 using LD-kNNi, an imputation method based on the k -nearest neighbor algorithm (Bradbury et al. 2007 ; Money et al. 2015 ). Genotypes were converted from a nucleotide-based format to a parent-based format, and consensus genotypes were called for both parents from the replicated samples using the GenosToABHPlugin in TASSEL 5 (Bradbury et al. 2007 ). Further marker analysis and filtering was performed in R using the ‘qtl’ and ‘ASMap’ packages (Broman et al. 2003 ; R Core Team 2022; Taylor and Butler 2017 ). Redundant markers with identical genotype patterns across all individuals were identified, and all but one representative marker were removed from the data set. Using the filtered dataset, all pairwise marker similarities were calculated among the individuals. Individuals with identical genotypes at ≥ 95% markers were identified, and the line with a greater amount of missing data among the pair was removed from analysis. Lastly, markers were tested for segregation distortion. Molecular marker screening The parents and population were screened with a series of markers previously described as linked to resistance loci: Sr7a , Sr11 , Sr12 , Sr13 , Sr21 , Sr38/Lr37/Yr17 , Sr63 , Sr8155B1 , and SrTm4 (Table 1 ). Additionally, shortly after population development concluded, MN07098-6 was found to be heterogeneous for Sr57/Lr34/Yr18 on chromosome 7D. The plant used for the backcross was not the selfed progeny of the original MN07098-6 crossing parent. Therefore, the lines used for the initial and backcross were potentially near-isogenic for Sr57/Lr34/Yr18 . The population was screened for Sr57/Lr34/Yr18 to determine if the population was segregating at this locus using a perfect marker (Table 1 ). Aneuploidy and chromosomal deletion detection In silico detection of aneuploidy and chromosomal deletions among the progeny lines was conducted using the GBS data and a method adapted from Singh et al. ( 2020 ). MN07098-6 and CItr 11390 were used as euploid allohexaploid and allotetraploid controls, respectively. Each chromosome was divided in 100 Mb bins with a sliding window step size of 50 Mb. Bins < 50 Mb in length at the end of chromosomes were discarded as they nested completely within another bin. BEDtools multicov (Quinlan and Hall 2010 ) was used to count the number of reads uniquely aligned within each bin. Counts were normalized across the genome for all bins and individuals to remove bias due to differential sequencing. Karyotypes for each sample were plotted using base R functions and the ‘chromPlot’ package (Oróstica and Verdugo 2016 ; R Core Team 2022). Bins were assigned a copy number based on normalized values. Due to artifacts in the assembly leading to a small number of reads aligning to D genome chromosomes for CItr 11390, the mean normalized value across all bins deviated slightly from one. Therefore, the upper bin value limit for each copy number was determined using the following equation: $${bin}_{max }=\frac{copy number* \mu }{2}+2\sigma$$ where µ and σ are the mean normalized values and standard deviation across all bins, respectively. Karyotype plots were visually assessed to detect chromosomal deletions and aneuploidy. Individuals with signatures for any form of aneuploidy involving chromosomes of either the A or B genome were not included in genetic and phenotypic analysis. Additionally, RILs with all forms of aneuploidy or deletions other than monosomy in the D genome were excluding from further analysis. RILs with monosomy affecting a single D genome chromosome were included in initial data analysis as there is representation of all 21 chromosomes. Linkage map construction Linkage map construction was performed in R using the MSTmap algorithm implemented in the package ‘ASMap’ (R Core Team 2022; Taylor and Butler 2017 ; Wu et al. 2008 ). Due to the population being developed from a single BC 1 F 1 plant, it was treated as a self-pollinated RIL population despite the generation of backcrossing. Linkage groups were independently constructed for each chromosome based on the physical positions assigned to the GBS markers. Map construction and marker ordering was conducted using the Kosambi mapping function and a significance threshold for marker clustering of p = 10 − 6 (Kosambi 1943 ). Final marker order was set by jittering the positions of markers mapped to the same position. Data analysis Seedling-stage stem rust resistance Mean values calculated from the converted infection type scores across the four replications for Pgt races observed to have variable virulence to the parents were used for analysis and QTL mapping. Lines with a linearized score of < 6 were considered resistant whereas those with scores ≥ 6 were considered susceptible. A chi-square (χ 2 ) test was conducted to determine the genetic architecture of resistance to the screened Pgt races within the population. Field-based stem rust The population was determined to be segregating for Sr57/Lr34/Yr18 . To determine the significance of the effects of Sr57/Lr34/Yr18 and ploidy on stem rust CI, mixed linear models were fit in R using the ‘ lme4 ’ package (Bates et al. 2015 ). Sr57/Lr34/Yr18 , ploidy state, and environments treated as fixed effects. Genotypes, genotype-environment interactions, and blocks within environments were treated as random effects. Backwards elimination of all effects in the across-environment model was performed using the ‘ lmerTest ’ package in R (Kuznetsova et al. 2017). Significance values for the fixed and random effects were calculated based on Sattethwaite’s approximation and likelihood ratio tests, respectively. The interaction between genotype and environment was determined to be significant, which warranted the separate analysis of each environment. CI values within environments were then adjusted for the effect of Sr57/Lr34/Yr18 by fitting mixed linear models for each of the four environments separately. Sr57 / Lr34 / Yr18 was significant in the across-environment model, while ploidy was determined to be insignificant. Therefore, Sr57 / Lr34 / Yr18 was the only fixed effect retained in the individual environment models. Genotype and block were included as random effects. In addition to visual inspection of diagnostic plots for each single environment model, residuals were tested for normality and homogeneity with Sharpiro-Wilk and Levene’s tests, respectively. If necessary, data transformation was performed if assumptions of normality and residual homogeneity were violated. Log and square root transformations were tested for each environment. Additionally, an arcsine transformation of the Eth20 data was tested after dividing the data by 100. The model selection approach above was repeated to assess the effectiveness of each transformation. Coefficients from the selected models were extracted to correct phenotypic values for the effect of Sr57/Lr34/Yr18 . Pearson correlation coefficients between environments were calculated in R using the adjusted datasets. QTL Mapping QTL detection within individual environments was conducted using multiple QTL mapping (MQM) in R using the ‘ qtl ’ package (Broman et al. 2003 ). Multiple QTL mapping was conducted using an automated forward/backward analysis, and the final models were selected based on penalized LOD score (pLOD; Manichaikul et al. 2009 ). First, a single locus model was fitted for each environment and race with segregating resistance using multiple imputation and the scanone function. A 95% significance threshold for each single locus scan was determined by 1000 permutations across the entire genome. Significant markers (α = 0.05) from the single locus models were used as the initial set of cofactors for MQM analysis. Next, a two-dimensional genome scan was performed with 1000 permutations to calculate the threshold for significant LOD scores (α = 0.05) described by Broman et al. ( 2003 ) using the scantwo function. These thresholds were then used to calculate the pLOD scores to be used for multiple QTL mapping. Using the loci detected by single-QTL scans and the calculated pLOD thresholds, Multiple QTL models with maximal pLOD scores were selected and refined using the stepwiseqtl command and multiple imputation. Within the stepwiseqtl function, forward selection was allowed to proceed with a maximum of 10 QTL. The reported LOD support intervals, reported as 95% Bayes credible intervals, were calculated using the command bayesint. Assessment of QTL interactions and combinations All pairwise QTL x QTL interactions were simultaneously tested during multiple QTL model selection using stepwiseqtl, and models were refined in the event of a significant interaction between QTL. Significance and percent of phenotypic variance explained of the final models including interactions was determined in R using the fitql command from the ‘qtl’ package. Mean CI values of RILs carrying all combinations of QTL detected in multiple environments, as well as Sr57 / Lr34 / Yr18 , were compared within each environments. Mean CI values of each class of RILs were calculated based on the allele at the peak marker of the QTL using the unadjusted and untransformed data from each environment. Tukey’s HSD tests were performed to determine whether the means of each class of RIL were significantly different within each environment. Molecular marker development PCR Allele Competitive Extension (PACE™) SNP genotyping assays (3CR Bioscience, Harlow, UK) were developed for novel QTL detected in multiple environments. Primers were designed at GBS SNPs throughout the LOD intervals. Primers targeting resistance alleles were labelled with the HEX fluorophore, and those for susceptible alleles were labelled with the FAM fluorophore. Each PACE™ assay mix contained 60 µL HEX-labeled primer (100 µM), 60 µL FAM-labeled primer (100 µM), 120 µL common primer (100 µM), and 230 µL molecular-grade water (Sigma Cat. No. 4502). For each assay, genotyping was conducted using standard conditions for PACE. 8 µL reactions consisted of 2 µL template DNA (20 ng/µL), 0.11 µL assay mix, 4 µL 2x PACE™ master mix, and 1.89 µL molecular grade water. Genotyping was performed on a Roche Lightcycler® 480 (Roche Diagnostics, Indianapolis, IN, USA) using the following reaction conditions: hot-start activation at 95°C for 15 minutes, 10 touchdown cycles at 94°C for 20 seconds, touchdown at 61°C for 1 minute with a 0.6°C reduction per cycle until a secondary target of 55°C, and 48 cycles at 94°C for 20 seconds followed by 55°C for 1 minute/cycle. An endpoint read was taken following a 2-minute cooldown cycle at 37°C. Alleles were called using the Roche Lightcycler® 480 Endpoint Genotyping Software Module. Results Genotyping and Aneuploid Detection Genotyping-by-sequencing (GBS) resulted in a total of 269,278,521 uniquely mapped reads across all samples. The distribution of read counts per sample was relatively normally distributed, and the median number of reads per sample was ~ 2.2 million (Fig. S1 A). Total read counts were first normalized across all samples. Counts in each bin were then normalized by the median count for each bin across all samples. As expected, the distribution of normalized read counts per bin had a normal distribution (Fig. S1 B). Small shoulder peaks were observed in the distribution centered around 0, 0.5, and 2.0, which represent deleted, monosomic, and duplicated bins, respectively. Over 98% of bins had a normalized read count around 1.0, representative of the disomic state. The analysis accurately assigned a copy number of 0 to all D genome chromosomes of the tetraploid parent, CItr 11390, while all 21 chromosomes in MN07098-6 were determined to be disomic and fully intact as expected (Fig. 1 ). Of the 121 RILs, 30 were determined to be aneuploids (Fig. S2 ). Of the 30 aneuploid RILs, 24 had an abnormality with chromosomes in homoeologous groups 4 or 6, with 18 possessing chromosome 6D abnormalities (Table S1 ). The most common form of aneuploidy was monosomy (15 out of 30 lines). Notes made in observation nurseries in 2019 in Saint Paul, Minnesota and Crookston, Minnesota corroborate the in silico analysis. Several of the RILs identified as aneuploids were noted in both locations to have poor vigor, short stature, and/or late maturity. The 30 RILs identified as aneuploids, as well as one RIL of a pair with > 95% identical SNPs, were initially retained during variant calling. However, the 16 RILs with forms of aneuploidy or deletions involving the A or B genome, as well as 1 euploid RIL of the pair found to have > 95% identical SNPs, were removed during SNP filtering. The remaining 14 aneuploid RILs (11 M6D and 3 M4D) were included in phenotypic data analysis to determine whether monosomy involving either chromosome had a significant effect on stem rust infection. Linkage map construction After filtering, 2,287 SNPs were obtained from GBS across the A and B genomes. No SNPs in the D genome remained after filtering as expected. In addition to the GBS SNPs, markers linked to previously mapped stem rust resistance genes were used for linkage mapping. Of the molecular markers assayed, those linked to Sr7a , Sr11 , Sr63 , and Sr8155B1 were polymorphic between the parents and added to the GBS marker dataset (Table 1 ). Additionally, the population was segregating for Sr57/Lr34/Yr18 . However, as no GBS SNPs were retained on chromosome 7D after filtering, Sr57/Lr34/Yr18 was not included in linkage mapping. Table 1 Molecular marker assays for known stem rust resistance genes Parent allele Gene Chromosome Marker MN07098-6 CItr 11390 Reference Sr7a 4AL KASP_IWB12146 KASP_IWB47019 KASP_IWB68386 G b A a Bajgain et al. 2015a ; Olivera et al. 2022 Sr11 6BL KASP_IWB72471 T a G b Nirmala et al. 2016 Sr12 3BL NB-LRR3 G b G b Hiebert et al. 2016 Sr13 6AL Sr13R/S_SNP A a Null maswheat.ucdavis.edu/protocols/Sr13 Sr21 2AL Sr21TRYF5R5 Null Null Chen et al. 2018 Sr38/Lr37/Yr17 2NS:2AS Lr37-Yr17-Sr38_GBG-KASP A a A a Helguera et al. 2003 ; Mustahsan 2022 Sr57/Lr34/Yr18 7DS wMAS000003 T/A c A a maswheat.ucdavis.edu/protocols/Lr34 Sr63 2AL KASP_IWB32429 T a G b Mago et al. 2022 Sr8155B1 6AS KASP_IWB10558 A a G b Nirmala et al. 2017 SrTm4 2AL BQ421276 DR732348 Null Null Briggs et al. 2015 a Reportedly predictive of susceptibility b Reportedly predictive of resistance c MN07098-6 is heterogeneous for Sr57 / Lr34 / Yr18 The 2,291 GBS and molecular markers used for linkage mapping were assigned to 21 linkage groups representing the 14 chromosomes of the A and B genomes. Chromosomes 1B, 2A, 2B, 3B, 5A, 5B, and 6A were represented by a single linkage group, while the other seven chromosomes were split into two linkage groups each. Following marker ordering, markers with identical genotypes were identified and all but one representative marker for each group of markers were removed from the map. The final map consisted of 1,212 markers with a total length of 1,230.7 cM and an average spacing of ~ 1.0 cM between markers. Seedling-stage disease evaluation Parental infection types to the four Pgt races screened are shown in Table 2 . MN07098-6 and CItr 11390 were both susceptible at the seedling stage (IT = 33 + c to 3+) to the two Ug99 lineage races screened (TTKSK and TTKTT). The results of the TTKTT screen were unexpected, as CItr 11390 was found to carry the allele at KASP_IWB10558 reported to be predictive of resistance provided by Sr8155B1 (Nirmala et al. 2017 ). Note that Nirmala et al. ( 2017 ) reported the same isolate of TTKTT used in this study, 14KEN58-1, to be avirulent to Sr8155B1 with infection type ‘0;’. However, the same study reported the 04KEN156/04 isolate to be virulent to Sr8155B1 with infection type ‘3+’. Initial inspection of the data indicated CItr 11390 was moderately resistant to TTRTF (IT = 2). However, analysis of the RILs found only four resistant lines (mean linearized IT < 6) to 100 susceptible lines (mean linearized IT ≥ 6). Additionally, the upper bound of the 95% confidence interval for all four RILs with a mean linearized IT < 6, was greater than the threshold for resistance (Fig. S4). Therefore, TTRTF was not included in QTL mapping. Table 2 Isolate designation, origin, parental infection type, and virulence phenotype of Puccinia graminis f. sp. tritici races used to evaluate seedling resistance in the population. Race Isolate Origin MN07098-6 CItr 11390 TTKSK 04KEN156/04 Kenya 3+ 3+ TTKTT 14KEN58-1 Kenya 3+ 33 + c TKTTF 13ETH18-1 Ethiopia 12- 3+ TTRTF 14GEO189-1 Georgia 3- 2 Variation between the parents for seedling resistance to TKTTF was observed. MN07098-6 had a mean linearized IT of 2.7. This finding was as expected given MN07098-6 carries Sr7a , a gene known to confer resistance to TKTTF (Edae et al. 2018 ; Turner et al. 2016 ). CItr 11390 was susceptible to TKTTF with a mean linearized IT of 8.8 (Fig. S4). The susceptibility of CItr 11390 was notable as it was found to carry the allele reported to be predictive of resistance at KASP_IWB72471, a marker linked to Sr11 (Nirmala et al. 2016 ). Analysis of the RILs found an approximately 1:1 segregation ratio of resistance with 57 resistant RILs and 47 susceptible RILs supporting the hypothesis that resistance to TKTTF in the population is provided by a single resistance gene (Table 3 ). A chi-square test resulted in the acceptance of the null hypothesis that resistance to TKTTF in the population is monogenic (χ 2 = 1.0, p > 0.32; Table 3 ). Table 3 Chi-square test for segregation of stem rust reaction to TKTTF in the 15xR012 population Race No. of resistant RILs (IT ≤ 6) No. of susceptible RILs (IT > 6) Observed segregation χ 2 Number of expected genes TKTTF 57 47 1:1 (1, N = 104) = 1.0, p > 0.33 1 gene Linkage mapping for seedling resistance to TKTTF MQM analysis was performed to validate the effectiveness of Sr7a against TKTTF. A large-effect QTL was detected on 4AL (LOD = 29.2; R 2 = 69.8%; Table 4 ; Fig. S5) that colocated with the Sr7a -linked KASP markers screened on the population. Unexpectedly, a second resistance loci was retained in the final model on the long arm of chromosome 5A (LOD = 29.2; R 2 = 4.6%; Table 4 ; Fig. S5). Like Sr7a , the resistance allele at this locus is donated by the hexaploid parent, MN07098-6. Table 4 Quantitative trait loci (QTL) for seedling resistance to TKTTF in the 15xR012 population detected by multiple QTL mapping (MQM). All QTL are significant at α < 0.05 Race QTL a Chr Flanking markers b Pos (cM) LOD c R 2 d Add e Left Right TKTTF QSr.umn-4A 4AL S4A_720268179 S4A_728943337 7.0 29.2 69.8 2.7 QSr.umn-5A.1 5AL S5A_553400727 S5A_605408042 109.2 3.6 4.6 0.6 a Named according to McIntosh et al. ( 2020 ) b 95% Bayes credible interval c Peak logarithm of odds d Phenotypic variance explained e Estimated additive effect of QTL; negative value indicates that the allele donor is CItr 11390 Field-based disease evaluation CItr 11390 displayed high levels of resistance in Kenya in both 2019 and 2020, but the accession was far less resistant in both Ethiopian nurseries (Table 5 ). MN07098-6 was susceptible in all four environments with mean severity ranging from 60 to 75% within individual environments (Table S2 ). Disease pressure was especially high in Ethiopia in 2019, and scoring was conducted later than the optimal window. Both parents and the RIL population displayed completely susceptible infection responses and almost no variation in severity was observed (Tables S2 and S3). Raw stem rust CI distributions were continuous and positively-skewed in Kenya in 2019 and 2020 (Fig. 2 ). Visual assessment of a Q-Q plot indicated the Eth20 data had a slightly negative skew. Residuals for Ken19 and Ken20 were in violation of the assumption of normality based on Shapiro-Wilk tests ( p < 0.05). Only the residuals for Ken20 were heteroscedastic based on Levene’s test ( p < 0.05). Table 5 Mean stem rust coefficient of infection values, standard deviations, and ranges in the 15xR012 mapping population in four field environments. Environments a Parent mean coefficient of infection Population CItr 11390 MN07098-6 Mean ± SD Range Eth19 60.0 60.0 58.3 ± 3.1 45.0–60.0 Eth20 40.0 65.0 59.2 ± 12.3 24.0–80.0 Ken19 2.5 60.0 35.5 ± 22.4 2.0–90.0 Ken20 8.3 75.0 41.0 ± 21.8 8.0-100 a Eth19 = Debre Zeit, Ethiopia 2019; Eth20 = Debre Zeit, Ethiopia 2020; Ken19 = Njoro, Kenya 2019; Ken20 = Njoro, Kenya 2020 Square root transformation of the Ken19 values resulted in the restoration of normality. Square root transformation of the Ken20 values also resulted in the restoration of normality and homoscedasticity, as well as marked improvement in the Q-Q plot. Therefore, the transformed data were used for Ken19 and Ken20.The Eth20 data were divided by 100 followed by an arcsine transformation. While a Shapiro-Wilk test indicated that the transformed data were normally distributed, little to no differences were viewed between the Q-Q plots of the untransformed and transformed data fitted with the same model. Given these observation, as well as the homogeneity of the untransformed residuals and the increased difficulty of interpreting QTL mapping results from transformed datasets, the untransformed data from Eth20 were used for further analyses. The 15xR012 population was found to segregate for Sr57/Lr34/Yr18 due to near-isogenic MN07098-6 lines being used for population development. Sr57 / Lr34 / Yr18 was found to have a significant effect on CI in a mixed model using the data from all environments. Due to a complete absence of additional markers linked to Sr57/Lr34/Yr18 on chromosome 7D, the effect of the gene was adjusted using mixed models fitted within each environment. The estimated marginal means of CI for the resistance allele at Sr57 / Lr34 / Yr18 was 30.0%, 46.5%, and 11.1% lower than the susceptible allele in Ken19, Ken20, and Eth20, respectively. After adjustment, mean CI increased slightly in Ken19, Ken20, and Eth20 (Fig. 3 ). Adjusted CI values of the RILs between all environments were significantly correlated (Table 6 ). Table 6 Pearson correlation coefficients ( r ) of adjusted stem rust coefficient of infection values observed in the 15xR012 population in four field environments in Ethiopia and Kenya in 2019 and 2020 Eth19 Eth20 Ken19 Ken20 Eth20 0.37 *** Ken19 0.23 * 0.33 *** Ken20 0.20 * 0.39 *** 0.66 *** Eth19 = Debre Zeit, Ethiopia 2019; Eth20 = Debre Zeit, Ethiopia 2020; Ken19 = Njoro, Kenya 2019; Ken20 = Njoro, Kenya 2020 * Correlation is significant at α < 0.05 ** Correlation is significant at α < 0.01 *** Correlation is significant at α < 0.001 Linkage mapping of stem rust resistance Multiple QTL mapping (MQM) detected significant QTL associated with APR in three of the four African nurseries (Table 7 ; Fig. S6). No QTL were detected in Eth19. In total, seven QTL were detected in at least one environment. The LOD scores of the detected QTL ranged from 3.3 to 20.9. Of the seven QTL, two were detected in more than one environment. QSr.umn-6BL was detected in Eth20 and Ken20. It explained 10.2% of the observed phenotypic variation in Eth20 and 10.5% in Ken20. While QSr.umn-6BL was not retained in the multiple QTL model, it was detected by single-QTL and two dimensional genome scans (Table S4). The second multi-environment QTL, QSr.umn-2A , was detected in Kenya in 2019 and 2020. QSr.umn-2A explained 42.8% of the phenotypic variation in 2019 and 44.6% of variation in 2020. The peaks for QSr.umn-2A were sharply centered around KASP marker Sr63_IWB32429 for Ken19 and Ken20, indicating that QSr.umn-2A could be Sr63 . QSr.umn-2A and QSr.umn-6BL were each contributed by the tetraploid parent, CItr 11390. The five other QTL were only detected in a single environment and explained 5.0-21.5% of the observed phenotypic variation within their respective environments. Notably, the interval of QSr.umn-5A.2 co-located with the TKTTF ASR-associated QTL, QSr.umn-5A.1 . The resistance allele was contributed by MN07098-6 at both QSr.umn-5A.1 and QSr.umn-5A.2 . Of the remaining single-environment QTL, three were contributed by MN07098-6 and one was contributed by CItr 11390 (Table 7 ). Table 7 Quantitative trait loci (QTL) for field-based stem rust resistance in the 15xR012 population detected by multiple QTL mapping (MQM). All QTL are significant at α < 0.05 Environment a QTL b Chr Flanking markers c Pos (cM) LOD d R 2 e Add f Left Right Eth20 QSr.umn-1A 1AS S1A_5914856 S1A_15479164 10.0 6.4 21.5 1.6 QSr.umn-6BL 6BL S6B_467576061 S6B_651836748 13.0 3.3 10.2 -1.1 Ken19 QSr.umn-2A 2AL S2A_681051011 S2A_702504398 139.0 14.1 42.8 -0.8 g QSr.umn-5A.2 5AL S5A_488255473 S5A_659345920 77.0 4.2 10.1 0.4 g Ken20 QSr.umn-2A 2AL S2A_681051011 Sr63_IWB32429 136.0 20.9 44.6 -0.67 g QSr.umn-3A 3AL S3A_603255953 S3A_654993190 105.0 3.6 5.0 0.22 g QSr.umn-3B 3BL S3B_429246336 S3B_490054756 33.1 6.6 10.0 0.32 g QSr.umn-6BS 6BS S6B_17383541 S6B_36319780 16.7 10.0 16.2 -0.22 g QSr.umn-6BL 6BL S6B_572334604 S6B_602404199 31.7 6.9 10.5 -0.31 g a Eth20 = Debre Zeit, Ethiopia 2020; Ken19 = Njoro, Kenya 2019; Ken20 = Njoro, Kenya 2020 b Named according to McIntosh et al. ( 2020 ) c 95% Bayes credible interval d Peak logarithm of odds e Phenotypic variance explained f Estimated additive effect of QTL; negative value indicates that the allele donor is CItr 11390 g Value is reported as the square root of coefficient of infection Assessment of QTL interactions and combinations Tests for epistasis while fitting the multiple QTL models detected one significant interaction within Ken20. The interaction between QSr.umn-2A and QSr.umn-6BS explained 10.9% of the observed phenotypic variation in Ken20 ( p < 0.001). The resistance allele was donated by CItr 11390 at both QTL. The mean CI of lines carrying only the resistance allele at QSr.umn-2A was significantly lower than those carrying only the resistance allele at QSr.umn-6BS . However, the difference in mean CI between lines only carrying the resistance allele at QSr.umn-2A and those carrying resistance alleles at both loci was not significant (Fig. 4 a, 4 b, Table S5). QSr.umn-6BL was detected in two environments by MQM while QSr.umn-2A was detected in Kenya in both years. In addition to the detected QTL, the effect of Sr57/Lr34 / Yr18 was significant on CI in all three environments. Multiple RILs in each of the eight possible combinations of these three loci were found among the population: No QTL, Sr57 , QSr.umn-2A (hereafter referred to as Sr63 ), QSr.umn-6BL , Sr57 + Sr63 , Sr57 + QSr.umn-6BL , Sr63 + QSr.umn-6BL , and Sr57 + Sr63 + QSr.umn-6B . Significant differences were observed among the allelic combinations in all environments. In Eth20, the mean CI of RILs carrying a single gene was lower than lines lacking all three genes. However, these differences were not significant. All combinations of two or all three resistance alleles resulted in a significant reduction in mean CI (Fig. 5 , Table S6). The mean unadjusted CI of lines carrying all three genes was 52.0 compared to 69.3 for lines lacking all three genes. Within the Kenyan environments, many of the pairwise comparisons were not statistically significant. However, this is possibly due to the low sample size of RILs constituting each of the eight groups of gene combinations (Fig. 5 , Table S6). Despite this finding, additive relationships were observed among the three genes. The mean unadjusted CI values of RILs carrying no resistance alleles at the three loci (Ken19: 61.6; Ken20: 66.9) were over three times greater than those carrying the resistance allele each of the three genes in both environments (Ken19: 18.2; Ken20: 18.5). Molecular marker development PACE ™ SNP genotyping assays were developed for QSr.umn-6BL at loci approximately flanking the QTL peaks detected by single-QTL genome scans in Eth20, Ken19, and Ken20. The markers with the clearest results were determined to be at loci correspond to positions 655,795,233 and 675,524,532 on chromosome 6B of RefSeq v2.1 assembly (Table S7). Genotyping calls for both assays were found to be concordant with the GBS markers (Table S8). Discussion The discovery and deployment of new stem rust resistance genes from diverse sources is vital to the sustained production of wheat on both regional and global scales. Since the initial discovery of Pgt race TTKSK (Ug99) in Uganda in 1998 (Pretorius et al. 2000 ), an additional fourteen variants have been identified with unique virulence profiles. The Ug99 race group has moved well beyond Uganda, with variants detected in fourteen countries spanning a geographic range from Iraq to South Africa (Terefe et al. 2018; Pretorius et al. 2012 ). The Ug99 race group is not the only threat posed by Pgt . The need for continued monitoring and resistance breeding for other Pgt races is evidenced by the detection and geographic expansion of races TKTTF, TTRTF, JRCQC, and TTTTF, some of which are virulent to widely deployed resistance genes including SrTmp and Sr13b (Letta et al. 2013 ; Olivera et al. 2015 ; Patpour et al. 2020 ). The challenge of discovering novel sources of stem rust resistance for wheat is aided by the rich gene pools available to wheat breeders and pathologists (Rowland and Kerber 1974 ). However, as demonstrated by this study, the integration of ancestral wheat species and other members of the Aegilops tribe into downstream breeding is not a straightforward process (Leigh et al. 2022 ). Of the 121 BC 1 F5-derived RILs developed, 30 were aneuploids and removed from further phenotypic analysis and QTL mapping (Table S1 ). Stem rust infection was much greater in both Ethiopian environments than either Kenyan environment. It is possible this difference was driven by the warmer environment in Ethiopia, as temperature has been found to significantly impact the effectiveness of ASR and APR genes (Chen et al. 2018 ; Gao et al. 2019 ; McIntosh et al. 1995 ; Zhang et al. 2017 ). Alternatively, the differences could be driven by the presence of additional races in Ethiopia. Races TRTTF and JRCQC have previously been identified from isolates collected in the Debre Zeit nursery (Olivera et al. 2012 ). While the 15xR012 population was not screened for resistance to JRCQC, all RILs were susceptible to TRTTF at the seedling stage (Fig. S3 ). In addition to TRTTF, the population was susceptible to Ug99 races TTKSK and TTKTT. Of the Pgt isolates screened, only seedling resistance to TKTTF was segregating within the population. The results from our seedling stage screenings contradict previous reports of a pair of molecular markers, KASP_IWB72471 and KASP_IWB10558, being predictive of ASR genes Sr11 and Sr8155B1 , respectively (Nirmala et al. 2016 ; Nirmala et al. 2017 ). Sr11 , residing on chromosome 6BL, confers resistance to Pgt race TKTTF (Green et al. 1960 ; Nirmala et al. 2016 ). However, the CItr 11390 was susceptible to TKTTF despite carrying the allele at KASP_IWB72471 reported to be predictive of Sr11 resistance. Additionally, KASP_IWB72471 was found to have no effect on seedling resistance to TKTTF in single marker analysis, and Sr11 was not detected by MQM for TKTTF or APR resistance despite substantial marker coverage in the region. Similarly, CItr 11390 was found to carry the allele reported to be predictive of Sr8155B1 resistance at KASP_IWB10558 on chromosome 6AS. However, CItr 11390, as well as MN07098-6 and the entire RIL population, were susceptible to race TTKTT, which is avirulent to Sr8155B1. These results indicate that KASP_IWB72471 and KASP_IWB10558 are not predictive of resistance at their respective loci in this population. Initial analysis of the phenotypic data from seedling testing to race TKTTF indicated resistance was conferred by a single gene (Table 3 ). This was expected given the population segregates for Sr7a . However, MQM detected a second locus stretching approximately 50Mb on the long arm of chromosome 5A. Like Sr7a , the resistance allele at QSr.umn-5A.1 was donated by the common wheat parent, MN07098-6. No known ASR genes are located on chromosome 5AL. Association mapping of North American wheat breeding germplasm by Bajgain et al. ( 2015a ) detected three significant MTAs associated with seedling resistance to TKTTF on chromosome 5AL. Each of the three SNPs fall within the 95% Bayes credible interval of QSr.umn-5A.1 , and one is located less than 3Mb from the peak of QSr.umn-5A.1 . MN07098-6 was included among the association mapping panel used by Bajgain et al. ( 2015a ), and the line carries the resistance allele at each of the three SNPs. The detection of co-locating QTL and MTAs sourced from common germplasm suggests the interval on chromosome 5AL may contain a novel source of resistance to Pgt race TKTTF. Multiple QTL mapping detected a total of seven QTL for APR to Pgt in three of the four test environments. Among these QTL, two were detected in multiple environments that were both donated by CItr 11390. The largest effect QTL of the three, QSr.umn-2A , was detected in Kenya in 2019 and 2020. It is likely QSr.umn-2A is Sr63 , an APR locus first reported by Mago et al. ( 2022 ) from the durum wheat ( Triticum turgidum ssp. durum ) variety ‘Glossy Huguenot’. A KASP marker linked to Sr63 , KASP_IWB32429, was located at the peak of QSr.umn-2A in both environments. Additionally, the physical positions of the QSr.umn-2A interval (681–702 Mb) were nearly identical to that of Sr63 reported by Mago el al. (2022) in the RefSeq v2.1 assembly (683–696 Mb). The detected interval is distinct from the location of Sr21 (713 Mb), and molecular marker screening confirmed the lack of the resistance allele at Sr21 in CItr 11390. The detection of QSr.umn-2A and Sr63 in a common environment (Kenya 2019) further supports the validation of Sr63 as a unique APR locus on chromosome 2AL. Given the population’s apparent seedling susceptibility to Ug99 lineage races (Table 2 ), Sr63 on chromosome 2AL could be an important target for resistance breeding to the Ug99 race group given its relatively large effect on disease reduction at the adult growth stages. The lack of detection of Sr63 in Ethiopia indicates the gene may not be effective against other races present in Ethiopia, such as TRTTF and JRCQC. Screening the population in single-race nurseries could provide greater detail on efficacy of the gene to specific races. Additionally, APR genes have been found to provide pleiotropic resistance to the other cereal rusts, as well as other plant pathogens (Herrera-Foessel et al. 2014 ; Hiebert et al. 2010 ; Krattinger et al. 2009 ; Martinez et al. 2004; Rinaldo et al. 2017 ). Further disease resistance characterization of the 15xR012 RIL population to other pathogens, namely stripe rust (caused by P. striiformis f. sp. tritici ) and leaf rust (caused by P. triticina f. sp. tritici ) may be worth pursuing. QSr.umn-6BL on the long arm of chromosome 6B was also detected in more than one environment. While the QTL was not retained in the multiple QTL model for Ken19, it was initially detected by both single-QTL and two-dimensional genome scans for Ken19, Ken20, and Eth20. The ASR gene Sr11 resides on 6BL (Green et al. 1960 ), but molecular marker screening suggests that neither MN07098-6 nor CItr 11390 carry the resistance allele at Sr11 . Additionally, the nearest flanking marker of QSr.umn-6BL was more than 40 Mb from Sr11 . Bajgain et al. ( 2015a ) previously reported a marker trait association (MTA) on chromosome 6BL at IWB35697 for APR in both Kenya and Ethiopia. MTAs in close proximity to IWB35697 were also detected in winter and durum wheat diversity panels in multiple years in Kenya and Ethiopia (Megeressa et al. 2020; Yu et al. 2012 ). However, local alignment of the significant markers from each of the aforementioned studies, including IWB35697, to RefSeq v2.1 revealed they were outside of the interval of QSr.umn-6BL and close to Sr11 . The distance from Sr11 , susceptibility of CItr 11390 to TKTTF, and the lack of APR QTL reported in this region suggest QSr.umn-6BL is a novel APR locus (Yu et al. 2014 ). Given the potential novelty of QSr.umn-6BL , refinement of the interval and the development of near isogenic lines (NILs) should be pursued to determine whether the QTL is a valuable breeding target for stem rust adult plant resistance. Efforts to do so can be initiated using the genotyping assays developed in this study (Table S7). Additional assays that more closely approximate the intervals detected during MQM are needed. In addition to the two QTL identified in multiple environments, five QTL were detected in single environments. Of these, four were contributed by the hexaploid parent MN07098-6: QSr.umn-1A (Eth20), QSr.umn-3A (Ken20), QSr.umn-3B (Ken20), and QSr.umn-5A.2 (Ken19). QSr.umn-1A explained a large portion of the phenotypic variation observed in Ethiopia in 2020 ( R 2 = 21.5). However, this is likely inflated by the small population size, relatively little phenotypic variation in the Eth20 environment, the lack of effectiveness of other loci such as Sr63 , and the absence of ASR loci in the population besides Sr7a and QSr.umn-5A.1 . While no known APR gene has been reported on 1AS, other QTL within the first ~ 16 Mb of chromosome 1AS have been detected in multiple studies (Bajgain et al. 2015b ; Bansal et al. 2008 ; Bhavani et al. 2011 ; Yu et al. 2012 ). Notably, Bajgain et al. ( 2015b ) detected a QTL in another University of Minnesota hard red spring wheat line for field resistance in 2013 in Kenya that co-locates with QSr.umn-1A . QSr.umn-3A ( R 2 = 5.0) was detected in Ken20 on the long arm of chromosome 3A with resistance conferred by the allele from MN07098-6. Sr35 is the most well-known stem rust resistance loci on chromosome 3AL (Zhang et al. 2010 ). However, it is unlikely that QSr.umn-3A is Sr35 as MN07098-6 is susceptible to TTKSK, the nearest flanking marker of QSr.umn-3A is ~ 11Mb from Sr35 , and the pedigree of MN07098-6 lacks T. monococcum introgressions. QSr.umn-3A is also unlikely to be Sr27 , which is carried on a 3AL.3RS translocation derived from the rye ( Secale cereale L.) variety ‘Imperial’ (Marais 2001 ). Sr27 provides strong seedling resistance to TTKSK that likely would have been detected during screening. Additionally, the deployment of the 3AL.3RS translocation has been almost exclusive to triticale (Singh et al. 2011 ). An MTA within 5 Mb of QSr.umn-3A was reported by Letta et al. ( 2013 ). However, this association was detected in a durum wheat panel grown in an Ethiopian nursery. In common wheat, the only APR QTL reported to partially overlap with QSr.umn-3A was detected in a CIMMYT biparental population grown in Njoro, Kenya during the 2011 main season (Singh et al. 2013 ). Further testing is needed to verify QSr.umn-3A given the limited reports of stem rust APR QTL within or near the interval. QSr.umn-3B ( R 2 = 6.6) was detected in Ken20 in the centromeric region of the long arm of chromosome 3B. The interval is in close proximity to a KASP marker linked to an NB-LRR motif within the Sr12 locus that co-segregates with seedling stage resistance (Hiebert et al. 2016 ). The population is fixed for the allele linked to resistance at this Sr12 (Table 1 ). However, adult plant resistance to Ug99 lineage and North American Pgt races also co-segregates with Sr12 (Hiebert et al. 2016 ; Rouse et al. 2014 ). The co-localization of QSr.umn-3B with the aforementioned loci suggests it may be the same source of APR. Hiebert et al. ( 2016 ) also determined that the APR locus within Sr12 interval interacted with Sr57/Lr34/Yr18 . Due to the inability to include Sr57 / Lr34 / Yr18 in QTL mapping, such interactions could not be evaluated while fitting multiple QTL models. Development and testing of additional genetic resource, such as families of NILs carrying various allelic combinations at the two loci, is necessary to further evaluate the effect of QSr.umn-3B and determine whether the locus interacts with Sr57/Lr34/Yr18 . QSr.umn-5A.2 ( R 2 = 8.6) was detected in Ken19 on the long arm of chromosome 5A. During the 2013 off-season in Ethiopia, Bajgain et al. ( 2015a ) detected a single MTA associated with APR that is ~ 13Mb from the peak marker of QSr.umn-5A lines. QTL on the chromosome 5AL associated with APR in African stem rust nurseries have also been reported in CIMMYT germplasm (Bhavani et al. 2011 ). Additionally, Edae et al. ( 2018 ) detected MTAs near the MTA reported by Bajgain et al. ( 2015a ) for both infection response and disease severity to race QTHJC across three growing seasons in Minnesota using the same panel of spring wheat lines. As previously discussed, a QTL providing ASR to TKTTF was also detected in this study that co-localized with QSr.umn-5A.1 . In the case of both QTL, the resistance allele was provided by MN07098-6. Given QSr.umn.5A.1 was only detected in a single environment, further testing is needed to validate the QTL and also determine the mechanism for co-localization of ASR and APR much like the Sr12 locus. The last QTL detected in this study, QSr.umn-6BS ( R 2 = 16.2), was detected only in Ken20 and contributed by CItr 11390. The effect of QSr.umn-6BS was relative small compared to the other QTL detected in Ken20. Additionally, an interaction between this QTL and Sr63 was detected during MQM. The mean CI value for lines carrying the resistance allele at both QSr.umn-6BS and Sr63 was not significantly different from those only carrying the resistance allele at Sr63 . This suggests that the presence of the resistance allele at Sr63 masks the effect of QSr.umn-6BS and is sufficient to reduce disease symptoms. Given the observed masking effect, as well as the substantially larger effect and detection in multiple environments of Sr63 , there is little reason to prioritize further investigations and selection of QSr.umn-6BS . Stacking multiple APR genes is an effective approach to developing high levels of resistance to stem rust (Pretorius et al. 2017 ). The mean CI of RILs carrying three resistance alleles at Sr57/Lr34/Yr18, Sr63 , and QSr.umn-6BL was markedly lower across environments compared to those carrying resistance alleles at zero, one, or two of the loci (Fig. 5 , Table S6). However, developing pyramids of favorable stem rust resistance alleles from alien sources that are accessible for commercial breeding is challenging due to issues associated with interspecific hybridization such as linkage drag and chromosome pairing. Sources for introgression are limited by the frequency of the resistance allele within a donor species. For example, of the 155 CIMMYT durum lines Mago et al. ( 2022 ) screened for Sr63 , two were found to carry the resistance allele, making CItr 11390 only the third wheat accession known to carry the resistance allele at this locus. The RILs of the 15xR012 population are therefore the first hexaploid wheat lines known to carry Sr63 . This is also likely the case for QSr.umn-6BL , as this is the first report of QTL for stem rust APR within this region. Multiple euploid RILs carry the resistance allele at Sr7a and Sr57/Lr34/Yr18 in addition to Sr63 and QSr.umn-6BL (Table S9). The development and discovery of these lines demonstrates the feasibility of simultaneously mapping and transferring QTL from Khorasan wheat into hexaploid breeding materials. The RILs and SNP genotyping assays developed in this study provide value not only to targeted breeding efforts in areas where Ug99 stem rust races are present, but also facilitate the rapid deployment of newly validated and discovered loci such as Sr63 and QSr.umn-6BL to the greater wheat breeding community. Declarations Funding This project was supported by the USDA National Institute of Food and Agriculture competitive grant no. 2022-68013-36439, USDA-ARS appropriated project “Surveillance, Pathogen Biology, and Host Resistance of Cereal Rusts”, and the USDA-ARS National Plant Disease Recovery System. Author Contributions MF coordinated sequencing, constructed the genetic map, performed data analysis, and wrote the final manuscript. EC conducted molecular marker genotyping, aided in sequencing library preparation, and contributed to data analysis. ZK and AG coordinated field trials and oversaw data collection. MR conceived the study, developed the 15xR012 mapping population, facilitated seedling disease screenings, and aided in trial design. JA contributed to data analysis and provided suggestions to experimental design. All authors reviewed the manuscript. Acknowledgments This project was supported by the USDA National Institute of Food and Agriculture competitive grant no. 2022-68013-36439, USDA-ARS appropriated project “Surveillance, Pathogen Biology, and Host Resistance of Cereal Rusts”, and the USDA-ARS National Plant Disease Recovery System. Data Availability Phenotypic and genotypic data used for analysis are included in this article and supplementary files. References Bajgain P, Rouse M, Bulli P, et al (2015a) Association mapping of North American spring wheat breeding germplasm reveals loci conferring resistance to Ug99 and other African stem rust races. BMC Plant Biology 15 Bajgain P, Rouse MN, Bhavani S, Anderson JA (2015b) QTL mapping of adult plant resistance to Ug99 stem rust in the spring wheat population RB07/MN06113-8. Molecular Breeding 35 Bansal U, Bariana H, Wong D, et al (2014) Molecular mapping of an adult plant stem rust resistance gene Sr56 in winter wheat cultivar Arina. Theoretical and Applied Genetics 127:1441–1448 Bansal UK, Bossolini E, Miah H, et al (2008) Genetic mapping of seedling and adult plant stem rust resistance in two European winter wheat cultivars. Euphytica 164:821–828 Bates D, Mächler M, Bolker B, Walker S (2015) Fitting Linear Mixed-Effects Models Using lme4. Journal of Statistical Software 67:1–48 Bhavani S, Singh R, Argillier O, et al (2011) Mapping durable adult plant stem rust resistance to the race Ug99 group in six CIMMYT wheats. 2011 Borlaug Global Rust Initiative Technical Workshop 43–53 Bradbury PJ, Zhang Z, Kroon DE, et al (2007) TASSEL: software for association mapping of complex traits in diverse samples. Bioinformatics 23:2633–2635 Braun HJ, Atlin G, Payne T (2010) Multi-location testing as a tool to identify plant response to global climate change Briggs J, Chen S, Zhang W, et al (2015) Mapping of SrTm4 , a Recessive Stem Rust Resistance Gene from Diploid Wheat Effective to Ug99. Phytopathology® 105:1347–1354 Broman KW, Wu H, Sen Ś, Churchill GA (2003) R/qtl: QTL mapping in experimental crosses. Bioinformatics 19:889–890 Chen S, Zhang W, Bolus S, et al (2018) Identification and characterization of wheat stem rust resistance gene Sr21 effective against the Ug99 race group at high temperature. PLOS Genetics 14:e1007287 Edae EA, Pumphrey MO, Rouse MN (2018) A Genome-Wide Association Study of Field and Seedling Response to Individual Stem Rust Pathogen Races Reveals Combinations of Race-Specific Genes in North American Spring Wheat. Frontiers in Plant Science 9 Elshire RJ, Glaubitz JC, Sun Q, et al (2011) A Robust, Simple Genotyping-by-Sequencing (GBS) Approach for High Diversity Species. PLoS ONE 6:e19379 FAO (2022) Crop Prospects and Food Situation. Food and Agriculture Organization of the United Nations Gao L, Babiker EM, Nava IC, et al (2019) Temperature‐sensitive wheat stem rust resistance gene Sr15 is effective against Puccinia graminis f. sp. tritici race TTKSK. Plant Pathology 68:143–151 Gao L, Turner MK, Chao S, et al (2016) Genome Wide Association Study of Seedling and Adult Plant Leaf Rust Resistance in Elite Spring Wheat Breeding Lines. PLOS ONE 11:e0148671 Gill BS, Friebe B, Endo TR (1991) Standard karyotype and nomenclature system for description of chromosome bands and structural aberrations in wheat (Triticum aestivum). Genome 34:830–839 Gill BS, Raupp WJ (1987) Direct Genetic Transfers from Aegilops squarrosa L. to Hexaploid Wheat. Crop Science 27:445–450 Godfray HCJ, Beddington JR, Crute IR, et al (2010) Food Security: The Challenge of Feeding 9 Billion People. Science 327:812–818 Grassini P, Eskridge KM, Cassman KG (2013) Distinguishing between yield advances and yield plateaus in historical crop production trends. Nature Communications 4 Green GJ, Knott DR, Watson IA, Pugsley AT (1960) Seedling reactions to stem rust of lines of Marquis wheat with substituted genes for rust resistance. Canadian Journal of Plant Science 40:524–538 Hayes HK, Ausemus ER, Stakman EC, et al (1936) Thatcher Wheat Helguera M, Khan IA, Kolmer J, et al (2003) PCR Assays for the Lr37‐Yr17‐Sr38 Cluster of Rust Resistance Genes and Their Use to Develop Isogenic Hard Red Spring Wheat Lines. Crop Science 43:1839–1847 Herrera-Foessel SA, Singh RP, Lillemo M, et al (2014) Lr67/Yr46 confers adult plant resistance to stem rust and powdery mildew in wheat. Theoretical and Applied Genetics 127:781–789 Hiebert CW, Kolmer JA, Mccartney CA, et al (2016) Major Gene for Field Stem Rust Resistance Co-Locates with Resistance Gene Sr12 in ‘Thatcher’ Wheat. PLOS ONE 11:e0157029 Hiebert CW, Thomas JB, Mccallum BD, et al (2010) An introgression on wheat chromosome 4DL in RL6077 (Thatcher*6/PI 250413) confers adult plant resistance to stripe rust and leaf rust (Lr67). Theoretical and Applied Genetics 121:1083–1091 Kalous JR, Martin JM, Sherman JD, et al (2015) Impact of the D genome and quantitative trait loci on quantitative traits in a spring durum by spring bread wheat cross. Theoretical and Applied Genetics 128:1799–1811 Kimber G, Sears ER (1968) Nomenclature for the description of aneuploids in the Triticinae. 3rd International Wheat Genetics Symposium 468–473 Kislev ME (1982) Stem Rust of Wheat 3300 Years Old Found in Israel. Science 216:993–994 Knott DR (1962) The inheritance of rust resistance: IX. The inheritance of resistance to races 15B and 56 of stem rust in the wheat variety Khapstein. Canadian Journal of Plant Science 42:415–419 Knott DR (2000) The inheritance of stem rust resistance in Thatcher wheat. Canadian Journal of Plant Science 80:53–63 Kolmer JA (1996) Genetics of resistance to wheat leaf rust. Annual Review of Phytopathology 34:435–455 Kosambi DD (1943) The estimation of map distances from recombination values. Annals of Eugenics 12:172–175 Kosgey ZC, Edae EA, Dill-Macky R, et al (2021) Mapping and Validation of Stem Rust Resistance Loci in Spring Wheat Line CI 14275. Frontiers in Plant Science 11 Krattinger SG, Lagudah ES, Spielmeyer W, et al (2009) A Putative ABC Transporter Confers Durable Resistance to Multiple Fungal Pathogens in Wheat. Science 323:1360–1363 Kumar S, Fetch TG, Knox RE, et al (2021) Mapping of Ug99 stem rust resistance in Canadian durum wheat. Canadian Journal of Plant Pathology 43:599–611 Lanning SP, Blake NK, Sherman JD, Talbert LE (2008) Variable Production of Tetraploid and Hexaploid Progeny Lines from Spring Wheat by Durum Wheat Crosses. Crop Science 48:199–202 Laugerotte J, Baumann U, Sourdille P (2022) Genetic control of compatibility in crosses between wheat and its wild or cultivated relatives. Plant Biotechnology Journal 20:812–832 Leigh FJ, Wright TIC, Horsnell RA, et al (2022) Progenitor species hold untapped diversity for potential climate-responsive traits for use in wheat breeding and crop improvement. Heredity 128:291–303 Letta T, Maccaferri M, Badebo A, et al (2013) Searching for novel sources of field resistance to Ug99 and Ethiopian stem rust races in durum wheat via association mapping. Theoretical and Applied Genetics 126:1237–1256 Li H (2011) A statistical framework for SNP calling, mutation discovery, association mapping and population genetical parameter estimation from sequencing data. Bioinformatics 27:2987–2993 Li H, Handsaker B, Wysoker A, et al (2009) The Sequence Alignment/Map format and SAMtools. Bioinformatics 25:2078–2079 Mago R, Chen C, Xia X, et al (2022) Adult plant stem rust resistance in durum wheat Glossy Huguenot: mapping, marker development and validation. Theoretical and Applied Genetics 135:1541–1550 Manichaikul A, Moon JY, Sen S, et al (2009) A Model Selection Approach for the Identification of Quantitative Trait Loci in Experimental Crosses, Allowing Epistasi. Genetics 181:1077–1086 Marais GF (2001) An evaluation of three Sr27-carrying wheat × rye translocations. South African Journal of Plant and Soil 18:135–136 Martin M (2011) Cutadapt removes adapter sequences from high-throughput sequencing reads. EMBnet.journal 17:10 Martínez F, Niks RE, Singh RP, Rubiales D (2001) Characterization of Lr46, a Gene Conferring Partial Resistance to Wheat Leaf Rust. Hereditas 135:111–114 McFadden ES (1930) A successful transfer of emmer characters to vulgare wheat. Journal of the American Society of Agronomy 22:1020–1034 McIntosh RA, Dubcovsky J, Rogers WJ, et al (2020) Catalogue of Gene Symbols for Wheat: 2020 Supplement. Annual Wheat Newsletter 66 McIntosh RA, Wellings CR, Park RF (1995) Wheat Rusts: An Atlas of Resistance Genes Megerssa SH, Ammar K, Acevedo M, et al (2020) Multiple-Race Stem Rust Resistance Loci Identified in Durum Wheat Using Genome-Wide Association Mapping. Frontiers in Plant Science 11 Money D, Gardner K, Migicovsky Z, et al (2015) LinkImpute: Fast and Accurate Genotype Imputation for Nonmodel Organisms. G3 Genes|Genomes|Genetics 5:2383–2390 Mustahsan W (2022) Characterization of stripe rust resistance and yield-related traits in an ‘Overley’ x ‘Overland’ population. Thesis Nirmala J, Chao S, Olivera P, et al (2016) Markers Linked to Wheat Stem Rust Resistance Gene Sr11 Effective to Puccinia graminis f. sp. tritici Race TKTTF. Phytopathology 106:1352–1358 Nirmala J, Saini J, Newcomb M, et al (2017) Discovery of a Novel Stem Rust Resistance Allele in Durum Wheat that Exhibits Differential Reactions to Ug99 Isolates. G3 Genes|Genomes|Genetics 7:3481–3490 Olivera P, Newcomb M, Szabo LJ, et al (2015) Phenotypic and Genotypic Characterization of Race TKTTF of Puccinia graminis f. sp. tritici that Caused a Wheat Stem Rust Epidemic in Southern Ethiopia in 2013–14. Phytopathology® 105:917–928 Olivera PD, Jin Y, Rouse M, et al (2012) Races of Puccinia graminis f. sp. tritici with Combined Virulence to Sr13 and Sr9e in a Field Stem Rust Screening Nursery in Ethiopia. Plant Disease 96:623–628 Olivera PD, Villegas D, Cantero‐Martínez C, et al (2022) A unique race of the wheat stem rust pathogen with virulence on Sr31 identified in Spain and reaction of wheat and durum cultivars to this race. Plant Pathology 71:873–889 Olson EL, Rouse MN, Pumphrey MO, et al (2013a) Introgression of stem rust resistance genes SrTA10187 and SrTA10171 from Aegilops tauschii to wheat. Theoretical and Applied Genetics 126:2477–2484 Olson EL, Rouse MN, Pumphrey MO, et al (2013b) Simultaneous transfer, introgression, and genomic localization of genes for resistance to stem rust race TTKSK (Ug99) from Aegilops tauschii to wheat. Theoretical and Applied Genetics 126:1179–1188 Oróstica KY, Verdugo RA (2016) chromPlot: visualization of genomic data in chromosomal context. Bioinformatics 32:2366–2368 Patpour M, Justesen AF, Tecle AW, et al (2020) First Report of Race TTRTF of Wheat Stem Rust (Puccinia graminis f. sp. tritici) in Eritrea. Plant Disease 104:973 Peterson RF, Campbell AB, Hannah AE (1948) A diagrammatic scale for estimating rust intensity on leaves and stems of cereals. Canadian Journal of Research 26c:496–500 Poland JA, Brown PJ, Sorrells ME, Jannink J-L (2012) Development of High-Density Genetic Maps for Barley and Wheat Using a Novel Two-Enzyme Genotyping-by-Sequencing Approach. PLoS ONE 7:e32253 Pretorius ZA, Ayliffe M, Bowden RL, et al (2017) Advances in control of wheat rusts. Achieving sustainable cultivation of wheat 1:295–343 Pretorius ZA, Jin Y, Bender CM, et al (2012) Seedling resistance to stem rust race Ug99 and marker analysis for Sr2, Sr24 and Sr31 in South African wheat cultivars and lines. Euphytica 186:15–23 Pretorius ZA, Singh RP, Wagoire WW, Payne TS (2000) Detection of Virulence to Wheat Stem Rust Resistance Gene Sr31 in Puccinia graminis f. sp. tritici in Uganda. Plant Disease 84:203–203 Quinlan AR, Hall IM (2010) BEDTools: a flexible suite of utilities for comparing genomic features. Bioinformatics 26:841–842 Rinaldo A, Gilbert B, Boni R, et al (2017) The Lr34 adult plant rust resistance gene provides seedling resistance in durum wheat without senescence. Plant Biotechnology Journal 15:894–905 Roelfs AP (1985) Wheat and rye stem rust. The Cereal Rusts, Vol. 2 2:3–37 Roelfs AP, Martens JWM (1988) An international system of nomenclature for Puccinia graminis f.sp. tritici. Phytopathology 78:526–533 Roelfs AP, Singh RP, Saari EE (1992) Rust Diseases of Wheat: Concepts and Methods of Disease Management. Rouse MN, Talbert LE, Singh D, Sherman JD (2014) Complementary epistasis involving Sr12 explains adult plant resistance to stem rust in Thatcher wheat (Triticum aestivum L.). Theoretical and Applied Genetics 127:1549–1559 Rouse MN, Wanyera R, Njau P, Jin Y (2011) Sources of Resistance to Stem Rust Race Ug99 in Spring Wheat Germplasm. Plant Disease 95:762–766 Rowland GG, Kerber ER (1974) Telocentric mapping in hexaploid wheat of genes for leaf rust resistance and other characters derived from Aegilops squarrosa. Canadian Journal of Genetics and Cytology 16:137–144 Sharma HC, Gill BS (1983) Current status of wide hybridization in wheat. Euphytica 32:17–31 Singh N, Raupp J, Koo D-H, et al (2020) In-silico detection of aneuploidy and chromosomal deletions in wheat using genotyping-by-sequencing. Plant Methods 16 Singh RP, Hodson DP, Huerta-Espino J, et al (2011) The Emergence of Ug99 Races of the Stem Rust Fungus is a Threat to World Wheat Production. Annual Review of Phytopathology 49:465–481 Singh RP, Hodson DP, Jin Y, et al (2015) Emergence and Spread of New Races of Wheat Stem Rust Fungus: Continued Threat to Food Security and Prospects of Genetic Control. Phytopathology® 105:872–884 Singh RP, Huerta-Espino J (2003) Effect of leaf rust resistance gene Lr34 on components of slow rusting at seven growth stages in wheat. Euphytica 129:371–376 Singh S, Singh RP, Bhavani S, et al (2013) QTL mapping of slow-rusting, adult plant resistance to race Ug99 of stem rust fungus in PBW343/Muu RIL population. Theoretical and Applied Genetics 126:1367–1375 Stakman EC, Stewart DM, Loegering WQ (1962) Identification of physiologic races of Puccinia graminis var. tritici. E617 Stubbs RW, Prescott JM, Saari EE, Dubin HJ (1986) Cereal disease methodology manual Stuthman DD, Leonard KJ, Miller‐Garvin J (2007) Breeding Crops for Durable Resistance to Disease. 95:319–367 Suenaga K, Singh RP, Huerta-Espino J, William HM (2003) Microsatellite Markers for Genes Lr34/Yr18 and Other Quantitative Trait Loci for Leaf Rust and Stripe Rust Resistance in Bread Wheat. Phytopathology 93:881–890 Taylor J, Butler D (2017) R Package ASMap: Efficient Genetic Linkage Map Construction and Diagnosis. Journal of Statistical Software 79:1–29 Terefe TG, Visser B, and Pretorius ZA (2016) Variation in Puccinia graminis f. sp. tritici detected on wheat and triticale in South Africa from 2009 to 2013. Crop Protection 86:9–16 The TT (1973) Chromosome Location of Genes Conditioning Stem Rust Resistance Transferred from Diploid to Hexaploid Wheat. Nature New Biology 241:256–256 Turner MK, Jin Y, Rouse MN, Anderson JA (2016) Stem Rust Resistance in ‘Jagger’ Winter Wheat. Crop Science 56:1719–1725 Wu Y, Bhat PR, Close TJ, Lonardi S (2008) Efficient and Accurate Construction of Genetic Linkage Maps from the Minimum Spanning Tree of a Graph. PLoS Genetics 4:e1000212 Xu LS, Wang MN, Cheng P, et al (2013) Molecular mapping of Yr53, a new gene for stripe rust resistance in durum wheat accession PI 480148 and its transfer to common wheat. Theoretical and Applied Genetics 126:523–533 Yu L-X, Barbier H, Rouse MN, et al (2014) A consensus map for Ug99 stem rust resistance loci in wheat. Theoretical and Applied Genetics 127:1561–1581 Yu L-X, Morgounov A, Wanyera R, et al (2012) Identification of Ug99 stem rust resistance loci in winter wheat germplasm using genome-wide association analysis. Theoretical and Applied Genetics 125:749–758 Zhang W, Chen S, Abate Z, et al (2017) Identification and characterization of Sr13, a tetraploid wheat gene that confers resistance to the Ug99 stem rust race group. Proceedings of the National Academy of Sciences 114:E9483–E9492 Zhang W, Olson E, Saintenac C, et al (2010) Genetic Maps of Stem Rust Resistance Gene Sr35 in Diploid and Hexaploid Wheat. Crop Science 50:2464–2474 Zhu T, Wang L, Rimbert H, et al (2021) Optical maps refine the bread wheat Triticum aestivum cv. Chinese Spring genome assembly. The Plant Journal 107:303–314 Supplementary Files FigureS2.pdf UMN15xR012AfricanRustGeneticMap.xlsx UMN15xR012AfricanRustPhenotypes.xlsx UMN15xR012AfricanRustSupplementalsFINAL.docx Cite Share Download PDF Status: Under Revision Version 1 posted Reviewers agreed at journal 23 May, 2023 Reviewers invited by journal 23 May, 2023 Editor assigned by journal 23 May, 2023 First submitted to journal 22 May, 2023 You are reading this latest preprint version Research Square lets you share your work early, gain feedback from the community, and start making changes to your manuscript prior to peer review in a journal. As a division of Research Square Company, we’re committed to making research communication faster, fairer, and more useful. We do this by developing innovative software and high quality services for the global research community. Our growing team is made up of researchers and industry professionals working together to solve the most critical problems facing scientific publishing. Also discoverable on Platform About Our Team In Review Editorial Policies Advisory Board Help Center Resources Author Services Accessibility API Access RSS feed Manage Cookie Preferences © Research Square 2026 | ISSN 2693-5015 (online) Privacy Policy Terms of Service Do Not Sell My Personal Information {"props":{"pageProps":{"initialData":{"identity":"rs-2958205","acceptedTermsAndConditions":true,"allowDirectSubmit":false,"archivedVersions":[],"articleType":"Research Article","associatedPublications":[],"authors":[{"id":203244885,"identity":"1ec295fc-2c2a-4156-a7b9-ad06aa8b42e7","order_by":0,"name":"Max Fraser","email":"data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAAAZAAAAAyAQMAAABI0h/eAAAABlBMVEX///8AAABVwtN+AAAACXBIWXMAAA7EAAAOxAGVKw4bAAAAuElEQVRIiWNgGAWjYBAC9gYQaWBjxw/hMxPWwnMARFakJUs2kKblzGHGDQeI1iJ9+PGHj22HmY1vZKc9YKiwTmwgqIUvzUxyZls6n9mN3O0GDGfSCWux52EwY+Zts2YGatkmwdh2mAhbeNg/f/7bxsy4eQZIyz+itPAYSDOccWbcIAHS0kCcljLJHmAgS5x5u00i4Vi6MTEO2/zhBygq24G2fKixliWoBRUkkKZ8FIyCUTAKRgEuAAAMbzs84EiEAwAAAABJRU5ErkJggg==","orcid":"https://orcid.org/0009-0001-8090-336X","institution":"University of Minnesota Twin Cities College of Food Agricultural and Natural Resource Sciences","correspondingAuthor":true,"prefix":"","firstName":"Max","middleName":"","lastName":"Fraser","suffix":""},{"id":203244886,"identity":"98e68c3d-cd7e-4c9a-9888-597ff5893de1","order_by":1,"name":"Emily Conley","email":"","orcid":"","institution":"University of Minnesota Twin Cities College of Food Agricultural and Natural Resource Sciences","correspondingAuthor":false,"prefix":"","firstName":"Emily","middleName":"","lastName":"Conley","suffix":""},{"id":203244887,"identity":"de33f6f6-a034-4c2d-aaf8-2aae47c82c94","order_by":2,"name":"Zennah Kosgey","email":"","orcid":"","institution":"Kenya Agricultural \u0026 Livestock Research Organization","correspondingAuthor":false,"prefix":"","firstName":"Zennah","middleName":"","lastName":"Kosgey","suffix":""},{"id":203244888,"identity":"76e455fb-3151-47d1-89ce-f03fdaf689cc","order_by":3,"name":"Ashenafi Gemechu Degete","email":"","orcid":"","institution":"Ethiopian Institute of Agricultural Research","correspondingAuthor":false,"prefix":"","firstName":"Ashenafi","middleName":"Gemechu","lastName":"Degete","suffix":""},{"id":203244889,"identity":"9c55a812-76ce-4850-a87e-101d84b9899f","order_by":4,"name":"Matthew Rouse","email":"","orcid":"","institution":"USDA-ARS CDL: USDA-ARS Cereal Disease Laboratory","correspondingAuthor":false,"prefix":"","firstName":"Matthew","middleName":"","lastName":"Rouse","suffix":""},{"id":203244890,"identity":"b22caed3-8c77-4e5b-a147-4478d606b8cb","order_by":5,"name":"James Anderson","email":"","orcid":"https://orcid.org/0000-0003-4655-6517","institution":"University of Minnesota Twin Cities College of Food Agricultural and Natural Resource Sciences","correspondingAuthor":false,"prefix":"","firstName":"James","middleName":"","lastName":"Anderson","suffix":""}],"badges":[],"createdAt":"2023-05-19 21:49:09","currentVersionCode":1,"declarations":"","doi":"10.21203/rs.3.rs-2958205/v1","doiUrl":"https://doi.org/10.21203/rs.3.rs-2958205/v1","draftVersion":[],"editorialEvents":[],"editorialNote":"","failedWorkflow":false,"files":[{"id":37494299,"identity":"76c6b882-e392-4944-ba8b-b9339426f683","added_by":"auto","created_at":"2023-05-25 14:47:29","extension":"jpg","order_by":1,"title":"Figure 1","display":"","copyAsset":false,"role":"figure","size":104167,"visible":true,"origin":"","legend":"\u003cp\u003eKaryotypes of MN07098-6, CItr 11390, and four representative RILs from the 15xR012 population with various ploidy states: \u003cstrong\u003eA)\u003c/strong\u003ehexaploid parent MN07098-6, \u003cstrong\u003eB)\u003c/strong\u003etetraploid parent CItr 11390, \u003cstrong\u003eC)\u003c/strong\u003editelosomic 6AS 15xR012-1, \u003cstrong\u003eD) \u003c/strong\u003ehexaploid 15xR012-2, \u003cstrong\u003eE) \u003c/strong\u003enullisomic-tetrasomic 6D/6A 15xR012-13, and \u003cstrong\u003eF) \u003c/strong\u003enullisomic 6D 15xR012-84. Centromeres are represented by the black circles. Chromosome copy number is indicated by color.\u003c/p\u003e","description":"","filename":"1.jpg","url":"https://assets-eu.researchsquare.com/files/rs-2958205/v1/0f99508aa31c0862bae6fe7c.jpg"},{"id":37494300,"identity":"c8321241-0528-48b3-80bb-d9020a8bcfac","added_by":"auto","created_at":"2023-05-25 14:47:29","extension":"jpg","order_by":2,"title":"Figure 2","display":"","copyAsset":false,"role":"figure","size":80320,"visible":true,"origin":"","legend":"\u003cp\u003eDistribution of stem rust coefficient of infection values (CI) for the 15xR012 population within \u003cstrong\u003ea) \u003c/strong\u003eEthiopia 2019 (Eth19), \u003cstrong\u003eb) \u003c/strong\u003eEthiopia 2020 (Eth20) ,\u003cstrong\u003ec) \u003c/strong\u003eKenya 2019 (Ken19), and \u003cstrong\u003ed) \u003c/strong\u003eKenya 2020 (Ken20). Data presented are mean CI values for each RIL before and after adjustment for the effect of \u003cem\u003eSr57/Lr34/Yr18\u003c/em\u003e.\u003c/p\u003e","description":"","filename":"2.jpg","url":"https://assets-eu.researchsquare.com/files/rs-2958205/v1/ab5e9970838bb716d1f7cedf.jpg"},{"id":37494795,"identity":"a6eea825-194f-46a3-8030-ebb6c6feefdc","added_by":"auto","created_at":"2023-05-25 14:55:29","extension":"jpg","order_by":3,"title":"Figure 3","display":"","copyAsset":false,"role":"figure","size":96388,"visible":true,"origin":"","legend":"\u003cp\u003eDistribution of stem rust coefficient of infection values (CI) for the 15xR012 population within \u003cstrong\u003ea) \u003c/strong\u003eKenya 2019 (Ken19) and \u003cstrong\u003eb) \u003c/strong\u003eKenya 2020 (Ken20). Data presented are mean CI values for each RIL after square root transformation and adjustment for the effect of \u003cem\u003eSr57/Lr34/Yr18\u003c/em\u003e\u003c/p\u003e","description":"","filename":"3.jpg","url":"https://assets-eu.researchsquare.com/files/rs-2958205/v1/6c539a90d894c35bb0fc5f44.jpg"},{"id":37494302,"identity":"e574bea3-7011-4a27-aa2d-d49d1e08b9a8","added_by":"auto","created_at":"2023-05-25 14:47:29","extension":"jpg","order_by":4,"title":"Figure 4","display":"","copyAsset":false,"role":"figure","size":74851,"visible":true,"origin":"","legend":"\u003cp\u003eEpistatic interaction of QTL for stem rust coefficient of infection (CI) in the 15xR012 population. \u003cstrong\u003ea) \u003c/strong\u003eDistribution of the square root transformed CI values in Ken20 for each allelic combination of \u003cem\u003eQSr.umn-2A \u003c/em\u003eand \u003cem\u003eQSr.umn-6BS \u003c/em\u003eamong the RIL population. \u003cstrong\u003eb) \u003c/strong\u003eEffect plot between \u003cem\u003eQSr.umn-2A \u003c/em\u003eand \u003cem\u003eQSr.umn-6BS\u003c/em\u003e in Ken20. The A allele is contributed by hexaploid parent MN07098-6, and the B allele is contributed by tetraploid parent CItr 11390. Black diamonds in panels A and C indicate mean values. Groups that do not share letters above the plot are significantly different (\u003cem\u003ep \u003c/em\u003e\u0026lt; 0.05; Tukey’s HSD). Percentages after the QTL names and the interaction term are the percent of variation explained. Error bars indicate the 95% confidence interval\u003c/p\u003e","description":"","filename":"4.jpg","url":"https://assets-eu.researchsquare.com/files/rs-2958205/v1/9a6ba7260ff323cf4919dcbd.jpg"},{"id":37494301,"identity":"d1925cef-6fa1-418c-af0d-5f7e871290a8","added_by":"auto","created_at":"2023-05-25 14:47:29","extension":"jpg","order_by":5,"title":"Figure 5","display":"","copyAsset":false,"role":"figure","size":126900,"visible":true,"origin":"","legend":"\u003cp\u003eUnadjusted stem rust coefficient of infection among the 15xR012 RILs carrying each combination of \u003cem\u003eSr57/Lr34/Yr18\u003c/em\u003e,\u003cem\u003e Sr63\u003c/em\u003e, and \u003cem\u003eQSr.umn-6BL\u003c/em\u003e alleles within the three environments \u003cem\u003eQSr.umn-6BL\u003c/em\u003e was detected. White diamonds indicate mean CI of each QTL combination within an environment. Groups within an environment that do not share letters above the plot are significantly different (\u003cem\u003ep \u003c/em\u003e\u0026lt; 0.05; Tukey’s HSD).\u003c/p\u003e","description":"","filename":"5.jpg","url":"https://assets-eu.researchsquare.com/files/rs-2958205/v1/692e84f3e82c4cb61a3b1ed4.jpg"},{"id":37496168,"identity":"98e49361-b0df-48e3-aa87-2419b9086172","added_by":"auto","created_at":"2023-05-25 15:03:36","extension":"pdf","order_by":0,"title":"","display":"","copyAsset":false,"role":"manuscript-pdf","size":921155,"visible":true,"origin":"","legend":"","description":"","filename":"manuscript.pdf","url":"https://assets-eu.researchsquare.com/files/rs-2958205/v1/5183e01f-dfd5-4c54-ab50-f38f9351c6a3.pdf"},{"id":37494797,"identity":"1899e894-95f2-47c0-9f4a-532a627eb4c7","added_by":"auto","created_at":"2023-05-25 14:55:29","extension":"pdf","order_by":1,"title":"","display":"","copyAsset":false,"role":"supplement","size":140562,"visible":true,"origin":"","legend":"","description":"","filename":"FigureS2.pdf","url":"https://assets-eu.researchsquare.com/files/rs-2958205/v1/1763ff0116510727d0200cbb.pdf"},{"id":37496162,"identity":"bed2a67a-3a18-481d-92a2-f2868b354c13","added_by":"auto","created_at":"2023-05-25 15:03:29","extension":"xlsx","order_by":2,"title":"","display":"","copyAsset":false,"role":"supplement","size":430675,"visible":true,"origin":"","legend":"","description":"","filename":"UMN15xR012AfricanRustGeneticMap.xlsx","url":"https://assets-eu.researchsquare.com/files/rs-2958205/v1/93b4a14d05bc6bb079d42975.xlsx"},{"id":37496163,"identity":"478c0988-fe71-4d18-924e-23524940f136","added_by":"auto","created_at":"2023-05-25 15:03:29","extension":"xlsx","order_by":3,"title":"","display":"","copyAsset":false,"role":"supplement","size":21542,"visible":true,"origin":"","legend":"","description":"","filename":"UMN15xR012AfricanRustPhenotypes.xlsx","url":"https://assets-eu.researchsquare.com/files/rs-2958205/v1/1113a11c4a6257f7bf651ae6.xlsx"},{"id":37494305,"identity":"2f235ed4-10f9-47d3-81ac-d88e112828c1","added_by":"auto","created_at":"2023-05-25 14:47:29","extension":"docx","order_by":4,"title":"","display":"","copyAsset":false,"role":"supplement","size":159066,"visible":true,"origin":"","legend":"","description":"","filename":"UMN15xR012AfricanRustSupplementalsFINAL.docx","url":"https://assets-eu.researchsquare.com/files/rs-2958205/v1/9036dbd5fdc5aa169a7084d5.docx"}],"financialInterests":"","formattedTitle":"Direct hybridization facilitates the simultaneous identification and introgression of QTL for adult plant resistance to the Ug99 stem rust lineage from tetraploid Khorasan wheat to common wheat","fulltext":[{"header":"Key message","content":"\u003cp\u003eMultiple quantitative trait loci for adult plant resistance to \u003cem\u003ePgt\u0026nbsp;\u003c/em\u003eraces endemic to Kenya and Ethiopia were simultaneously identified and introgressed from tetraploid wheat accession CItr 11390 into common wheat.\u003c/p\u003e"},{"header":"Introduction","content":"\u003cp\u003eWorldwide, common wheat (\u003cem\u003eTriticum aestivum\u003c/em\u003e L.) is the most widely grown crop with approximately 220\u0026nbsp;million hectares harvested annually. Global production of ~\u0026thinsp;780\u0026nbsp;million tons of wheat provides 20% of the world\u0026rsquo;s calories and protein (FAO \u003cspan citationid=\"CR14\" class=\"CitationRef\"\u003e2022\u003c/span\u003e; Braun et al. \u003cspan citationid=\"CR8\" class=\"CitationRef\"\u003e2010\u003c/span\u003e). In order to sustain this demand, wheat yields must increase 60% by 2050 (Godfray et al. \u003cspan citationid=\"CR19\" class=\"CitationRef\"\u003e2010\u003c/span\u003e). Unfortunately, trends in wheat yields are inconsistent across major production regions and increasingly threatened by numerous biotic and abiotic factors (Grassini et al. \u003cspan citationid=\"CR20\" class=\"CitationRef\"\u003e2013\u003c/span\u003e). Among these biotic factors is wheat stem rust caused by the fungus \u003cem\u003ePuccinia graminis\u003c/em\u003e f. sp. \u003cem\u003etritici\u003c/em\u003e Eriks. and E. Henn. (\u003cem\u003ePgt\u003c/em\u003e).\u003c/p\u003e \u003cp\u003eCapable of destroying entire fields and reducing regional grain yields by up to 70% during epidemics, stem rust has been among the most damaging and feared diseases of wheat throughout history. Historical and archeological records of \u003cem\u003ePgt\u003c/em\u003e infection and stem rust epidemics extend across all major growing regions and date as far back as 3,300 BC (Kislev \u003cspan citationid=\"CR29\" class=\"CitationRef\"\u003e1982\u003c/span\u003e; Roelfs \u003cspan citationid=\"CR70\" class=\"CitationRef\"\u003e1985\u003c/span\u003e). However, outbreaks have greatly diminished since the early 20th century. Widespread deployment of genetic resistance to \u003cem\u003ePgt\u003c/em\u003e, disease control programs such as the Barberry Eradication Program in the United States, and the development of chemical fungicides effectively contained stem rust outbreaks to manageable levels throughout much of the 20th century (Singh et al. \u003cspan citationid=\"CR78\" class=\"CitationRef\"\u003e2011\u003c/span\u003e).\u003c/p\u003e \u003cp\u003eDuring this period, stem rust disease resistance was largely qualitative in nature. Race-specific, all-stage resistance (ASR), also referred to as seedling resistance, was the cornerstone of stem rust resistance breeding efforts. The release of the hard red spring wheat variety \u0026lsquo;Thatcher\u0026rsquo; by the University of Minnesota in 1934 was in direct response to a series of epidemics that afflicted the growing region in the preceding decades (Hayes et al. \u003cspan citationid=\"CR22\" class=\"CitationRef\"\u003e1936\u003c/span\u003e). The rise of \u0026lsquo;Thatcher\u0026rsquo; as the predominant variety grown throughout the Northern Great Plains of the United States and Canada into the 1960\u0026rsquo;s was driven by a stem rust resistance package built around the ASR gene \u003cem\u003eSr12\u003c/em\u003e (Knott \u003cspan citationid=\"CR31\" class=\"CitationRef\"\u003e2000\u003c/span\u003e). Beyond North America, stem rust was effectively controlled through the introgression and deployment of another ASR gene, \u003cem\u003eSr31\u003c/em\u003e, from the rye (\u003cem\u003eSecale cereale\u003c/em\u003e L.) variety \u0026lsquo;Petkus\u0026rsquo; beginning in the mid-1960s (McIntosh et al. \u003cspan citationid=\"CR50\" class=\"CitationRef\"\u003e1995\u003c/span\u003e). Genetic control of stem rust using ASR genes such as \u003cem\u003eSr12\u003c/em\u003e and \u003cem\u003eSr31\u003c/em\u003e remained effective throughout much of the mid to late 20th century. Subsequently, breeding efforts and funding for fundamental and applied stem rust research were shifted towards the control of other pathogens (Singh et al. \u003cspan citationid=\"CR78\" class=\"CitationRef\"\u003e2011\u003c/span\u003e).\u003c/p\u003e \u003cp\u003eStem rust reemerged as a central topic of wheat research and breeding following the discovery of \u003cem\u003ePgt\u003c/em\u003e isolate Ug99 in Uganda in 1998 (Pretorius et al. \u003cspan citationid=\"CR67\" class=\"CitationRef\"\u003e2000\u003c/span\u003e) and subsequent wheat stem rust epidemics in Kenya in 2004 and 2005. Later given the race designation TTKSK (Roelfs and Martens \u003cspan citationid=\"CR71\" class=\"CitationRef\"\u003e1988\u003c/span\u003e), it was virulent to \u003cem\u003eSr31\u003c/em\u003e and several other resistance genes from common wheat and related species. Since the initial discovery of race TTKSK, the pathogen has undergone rapid genetic mutation and moved as far as South Africa, Yemen, and Iran (Terefe et al. 2018; Pretorius et al. \u003cspan citationid=\"CR66\" class=\"CitationRef\"\u003e2012\u003c/span\u003e). Further migration of the Ug99 lineage of \u003cem\u003ePgt\u003c/em\u003e poses a major threat to global wheat production, as a study found that 85\u0026ndash;90% of wheat varieties were susceptible to Ug99 (Singh et al. \u003cspan citationid=\"CR79\" class=\"CitationRef\"\u003e2015\u003c/span\u003e).\u003c/p\u003e \u003cp\u003eThe Ug99 lineage of stem rust is not the only threat to wheat production within Sub-Saharan Africa and beyond. Recent epidemics in Kenya and Ethiopia are attributed to races TKTTF, TTRTF, and JRCQC (Letta et al. \u003cspan citationid=\"CR40\" class=\"CitationRef\"\u003e2013\u003c/span\u003e; Olivera et al. \u003cspan citationid=\"CR56\" class=\"CitationRef\"\u003e2015\u003c/span\u003e; Patpour et al. \u003cspan citationid=\"CR62\" class=\"CitationRef\"\u003e2020\u003c/span\u003e). The movement of these races outside the regions of their initial detection and the outbreaks they have caused highlight the need for constant pathogen monitoring, discovery of novel wheat stem rust resistance genes, and an emphasis on breeding for adult plant resistance.\u003c/p\u003e \u003cp\u003eAdult plant resistance (APR) is characterized by a reduction in both urediniospore size and production, and it is generally race nonspecific unlike like ASR (Stuthman et al. \u003cspan citationid=\"CR84\" class=\"CitationRef\"\u003e2007\u003c/span\u003e). APR is believed to be a more durable form of resistance, as it is less prone to the boom and bust cycle associated with the successive deployment of individual ASR genes (Singh et al. \u003cspan citationid=\"CR78\" class=\"CitationRef\"\u003e2011\u003c/span\u003e). Unfortunately, only six of the 63 named stem rust resistance genes confer APR: \u003cem\u003eSr2\u003c/em\u003e (McFadden \u003cspan citationid=\"CR48\" class=\"CitationRef\"\u003e1930\u003c/span\u003e), \u003cem\u003eSr55\u003c/em\u003e/\u003cem\u003eLr67\u003c/em\u003e/\u003cem\u003eYr46\u003c/em\u003e (Herrera-Foessel et al. \u003cspan citationid=\"CR24\" class=\"CitationRef\"\u003e2014\u003c/span\u003e), \u003cem\u003eSr56\u003c/em\u003e (Bansal et al. \u003cspan citationid=\"CR3\" class=\"CitationRef\"\u003e2014\u003c/span\u003e), \u003cem\u003eSr57\u003c/em\u003e/\u003cem\u003eLr34\u003c/em\u003e/\u003cem\u003eYr18\u003c/em\u003e (Singh and Huerta-Espino \u003cspan citationid=\"CR80\" class=\"CitationRef\"\u003e2003\u003c/span\u003e), \u003cem\u003eSr58\u003c/em\u003e/\u003cem\u003eLr46\u003c/em\u003e/\u003cem\u003eYr29\u003c/em\u003e (Mart\u0026iacute;nez et al. \u003cspan citationid=\"CR47\" class=\"CitationRef\"\u003e2001\u003c/span\u003e), and \u003cem\u003eSr63\u003c/em\u003e (Mago et al. \u003cspan citationid=\"CR43\" class=\"CitationRef\"\u003e2022\u003c/span\u003e). Singh et al. (\u003cspan citationid=\"CR78\" class=\"CitationRef\"\u003e2011\u003c/span\u003e) hypothesized that pyramiding four to five APR genes could provide near immunity to \u003cem\u003ePgt\u003c/em\u003e. However, the expression and effectiveness of resistance genes is dependent on genetic background and environmental conditions (Kolmer \u003cspan citationid=\"CR32\" class=\"CitationRef\"\u003e1996\u003c/span\u003e). Therefore, discovery and characterization of novel genes are essential to optimize gene pyramids and overcome the ability of \u003cem\u003ePgt\u003c/em\u003e to defeat deployed resistance genes.\u003c/p\u003e \u003cp\u003eMapping and dissection of stem rust resistance genes has not been limited to hexaploid common wheat (AABBDD, 2n\u0026thinsp;=\u0026thinsp;6x\u0026thinsp;=\u0026thinsp;42). The progenitor species and wild relatives of common wheat have proven to be rich sources of diversity for numerous traits including stem rust resistance (Rowland and Kerber \u003cspan citationid=\"CR75\" class=\"CitationRef\"\u003e1974\u003c/span\u003e). In total, stem rust resistance genes have been identified in 14 relatives of hexaploid wheat. The resistance genes identified in these relatives account for a majority of the genes that remain effective against the races of the Ug99 lineage (Singh et al. \u003cspan citationid=\"CR79\" class=\"CitationRef\"\u003e2015\u003c/span\u003e). However, most of the introgressions are not used by breeding programs (Singh et al. \u003cspan citationid=\"CR79\" class=\"CitationRef\"\u003e2015\u003c/span\u003e; Yu et al. \u003cspan citationid=\"CR92\" class=\"CitationRef\"\u003e2014\u003c/span\u003e). Introgressions from species in the secondary and tertiary genes pools including \u003cem\u003eT. timopheevii\u003c/em\u003e Zhuk. (AAGG, 2n\u0026thinsp;=\u0026thinsp;4x\u0026thinsp;=\u0026thinsp;28), \u003cem\u003eSecale cereale\u003c/em\u003e L., and members of \u003cem\u003eThinopyrum\u003c/em\u003e and \u003cem\u003eAegilops\u003c/em\u003e are carried on Robertsonian translocations, which often harbor unfavorable alleles at genes influencing agronomic traits (Singh et al. \u003cspan citationid=\"CR78\" class=\"CitationRef\"\u003e2011\u003c/span\u003e).\u003c/p\u003e \u003cp\u003eAs opposed to members of the secondary and tertiary gene pools, the diploid and tetraploid progenitors of common wheat carry homologous chromosomes to those of common wheat capable of undergoing pairing and recombination, albeit with varying difficulty (Gill and Raupp \u003cspan citationid=\"CR18\" class=\"CitationRef\"\u003e1987\u003c/span\u003e). In the case of \u003cem\u003eT. monococcum\u003c/em\u003e (AA, 2n\u0026thinsp;=\u0026thinsp;2x\u0026thinsp;=\u0026thinsp;14), gene transfer is most often conducted through a bridge cross with durum wheat (\u003cem\u003eT. turgidum\u003c/em\u003e, AABB, 2n\u0026thinsp;=\u0026thinsp;4x\u0026thinsp;=\u0026thinsp;28) (Sharma and Gill \u003cspan citationid=\"CR76\" class=\"CitationRef\"\u003e1983\u003c/span\u003e). Direct hybridization with \u003cem\u003eT. monococcum\u003c/em\u003e is complicated by the high rates of male and female sterility and lethality in F\u003csub\u003e1\u003c/sub\u003e hybrids (The \u003cspan citationid=\"CR88\" class=\"CitationRef\"\u003e1973\u003c/span\u003e). However, direct introgression from \u003cem\u003eT. monococcum\u003c/em\u003e is possible through extensive backcrossing as demonstrated by the transfer of \u003cem\u003eSr21\u003c/em\u003e and \u003cem\u003eSr35\u003c/em\u003e into hexaploid germplasm (The \u003cspan citationid=\"CR88\" class=\"CitationRef\"\u003e1973\u003c/span\u003e; Zhang 2010). Introgressions from \u003cem\u003eAegilops tauschii\u003c/em\u003e (DD, 2n\u0026thinsp;=\u0026thinsp;2x\u0026thinsp;=\u0026thinsp;14) have historically been facilitated by first producing synthetic allohexaploid wheat through chromosome doubling of durum x \u003cem\u003eAe. tauschii\u003c/em\u003e F\u003csub\u003e1\u003c/sub\u003e hybrids. Direct crossing between \u003cem\u003eAe. tauschii\u003c/em\u003e and common wheat is a faster approach than developing synthetic allohexaploid lines despite the necessity of performing embryo rescue, and it also avoids the disruption of adaptive allele combinations that have evolved in the A and B genomes (Gill and Raupp \u003cspan citationid=\"CR18\" class=\"CitationRef\"\u003e1987\u003c/span\u003e). Olson et al. (\u003cspan citationid=\"CR59\" class=\"CitationRef\"\u003e2013a\u003c/span\u003e, \u003cspan citationid=\"CR60\" class=\"CitationRef\"\u003eb\u003c/span\u003e) demonstrated the feasibility of direct crossing with \u003cem\u003eAe. tauschii\u003c/em\u003e by simultaneously mapping and introgressing stem rust resistance loci \u003cem\u003eSrTA1662\u003c/em\u003e, \u003cem\u003eSr10171\u003c/em\u003e, and \u003cem\u003eSr10187\u003c/em\u003e in populations developed after two rounds of backcrossing to \u003cem\u003eT. aestivum\u003c/em\u003e recurrent parents.\u003c/p\u003e \u003cp\u003eDirect introgression of stem rust resistance genes from emmer and durum wheat (AABB, 2n\u0026thinsp;=\u0026thinsp;4x\u0026thinsp;=\u0026thinsp;28) is well documented and does not require specialized techniques (Lanning et al \u003cspan citationid=\"CR37\" class=\"CitationRef\"\u003e2008\u003c/span\u003e; Laugerotte et al. \u003cspan citationid=\"CR38\" class=\"CitationRef\"\u003e2022\u003c/span\u003e). Examples include \u003cem\u003eSr2\u003c/em\u003e introgressed from \u003cem\u003eT. turgidum\u003c/em\u003e ssp. \u003cem\u003edicoccum\u003c/em\u003e variety \u0026lsquo;Yaroslav\u0026rsquo; (McFadden \u003cspan citationid=\"CR48\" class=\"CitationRef\"\u003e1930\u003c/span\u003e), \u003cem\u003eSr11\u003c/em\u003e from \u003cem\u003eT. turgidum\u003c/em\u003e ssp. \u003cem\u003edurum\u003c/em\u003e variety \u0026lsquo;Gaza\u0026rsquo; (Green et al. \u003cspan citationid=\"CR21\" class=\"CitationRef\"\u003e1960\u003c/span\u003e), \u003cem\u003eSr12\u003c/em\u003e from \u003cem\u003eT. turgidum\u003c/em\u003e ssp. \u003cem\u003edurum\u003c/em\u003e variety \u0026lsquo;Iumillo\u0026rsquo; (Knott \u003cspan citationid=\"CR31\" class=\"CitationRef\"\u003e2000\u003c/span\u003e), and \u003cem\u003eSr13\u003c/em\u003e and \u003cem\u003eSr14\u003c/em\u003e from \u003cem\u003eT. turgidum\u003c/em\u003e ssp. \u003cem\u003edicoccum\u003c/em\u003e variety \u0026lsquo;Khapli\u0026rsquo; (Knott \u003cspan citationid=\"CR30\" class=\"CitationRef\"\u003e1962\u003c/span\u003e). Similar to the approach described by Olson et al., (\u003cspan citationid=\"CR59\" class=\"CitationRef\"\u003e2013a\u003c/span\u003e, \u003cspan citationid=\"CR60\" class=\"CitationRef\"\u003eb\u003c/span\u003e), QTL influencing yield component traits and \u003cem\u003eYr53\u003c/em\u003e, a gene for stripe rust (caused by \u003cem\u003eP. striiformis\u003c/em\u003e f. sp. \u003cem\u003etritici\u003c/em\u003e) resistance, were simultaneously mapped and transferred into common wheat from tetraploid wheat accessions (Kalous et al. \u003cspan citationid=\"CR27\" class=\"CitationRef\"\u003e2015\u003c/span\u003e; Xu et al. \u003cspan citationid=\"CR91\" class=\"CitationRef\"\u003e2013\u003c/span\u003e).\u003c/p\u003e \u003cp\u003eThe primary objective of this study was to identify the underlying loci conferring APR to Ug99 lineage races of \u003cem\u003ePgt\u003c/em\u003e in the Khorasan wheat accession (\u003cem\u003eTriticum turgidum\u003c/em\u003e ssp. \u003cem\u003eturanicum\u003c/em\u003e) CItr 11390 (\u0026lsquo;Sun Ray\u0026rsquo;). In an effort to simultaneously transfer these loci to hexaploid wheat lines while reducing the effects of poor agronomic alleles from CItr 11390 and recovering euploid progeny, a BC\u003csub\u003e1\u003c/sub\u003eF\u003csub\u003e5\u003c/sub\u003e-derived recombinant inbred line (RIL) population was developed.\u003c/p\u003e"},{"header":"Materials and methods","content":"\u003cdiv id=\"Sec3\" class=\"Section2\"\u003e \u003ch2\u003ePlant materials\u003c/h2\u003e \u003cp\u003eA mapping population of 121 BC\u003csub\u003e1\u003c/sub\u003eF\u003csub\u003e5\u003c/sub\u003e-derived RILs was developed from an interspecific cross between tetraploid accession CItr 11390 and the hexaploid recurrent parent MN07098-6. CItr 11390 (\u0026lsquo;Sun Ray\u0026rsquo;) is a Khorasan wheat (\u003cem\u003eTriticum turgidum\u003c/em\u003e ssp. \u003cem\u003eturanicum\u003c/em\u003e) accession of unknown improvement status within the USDA National Small Grains Collection. MN07098-6 (SD3696/MN97803-3BS) is a hard red spring wheat breeding line susceptible to the Ug99 lineage from the University of Minnesota. The population was developed from a single BC\u003csub\u003e1\u003c/sub\u003eF\u003csub\u003e1\u003c/sub\u003e plant. The BC\u003csub\u003e1\u003c/sub\u003eF\u003csub\u003e2\u003c/sub\u003e progeny were advanced in the greenhouse using the single seed descent method, and a seed increase was conducted at the BC\u003csub\u003e1\u003c/sub\u003eF\u003csub\u003e5\u003c/sub\u003e generation to obtain enough seed for field experiments.\u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec4\" class=\"Section2\"\u003e \u003ch2\u003eSeedling stage stem rust phenotyping\u003c/h2\u003e \u003cp\u003eThe RILs and parents were evaluated for seedling infection types to \u003cem\u003ePgt\u003c/em\u003e races TTKSK (isolate 04KEN156/04; Ug99), TTKTT (isolate 14KEN58-1; Ug99 with \u003cem\u003eSr24\u003c/em\u003e and \u003cem\u003eSrTmp\u003c/em\u003e virulence), TKTTF (isolate 13ETH18-1), and TTRTF (isolate 14GEO189-1) at the United States Department of Agriculture \u0026ndash; Agricultural Research Service (USDA-ARS) Cereal Disease Lab in Saint Paul, Minnesota. Inoculation and phenotyping procedures were conducted as described by Rouse et al. (\u003cspan citationid=\"CR74\" class=\"CitationRef\"\u003e2011\u003c/span\u003e). All lines were evaluated in four replications. Infection types (ITs) were recorded on a 0\u0026ndash;4 scale described by Stakman et al. (\u003cspan citationid=\"CR82\" class=\"CitationRef\"\u003e1962\u003c/span\u003e). IT scores were converted to a linearized scale for QTL mapping according to Gao et al. (\u003cspan citationid=\"CR16\" class=\"CitationRef\"\u003e2016\u003c/span\u003e).\u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec5\" class=\"Section2\"\u003e \u003ch2\u003eField-based stem rust evaluation\u003c/h2\u003e \u003cp\u003eThe BC\u003csub\u003e1\u003c/sub\u003eF\u003csub\u003e5:6\u003c/sub\u003e and BC\u003csub\u003e1\u003c/sub\u003eF\u003csub\u003e5:7\u003c/sub\u003e populations, as well as the parents, were evaluated for field response to African \u003cem\u003ePgt\u003c/em\u003e races during the off-seasons (January \u0026ndash; May) of 2019 and 2020 in Njoro, Kenya (hereafter referred to as Ken19 and Ken20, respectively) and Debre Zeit, Ethiopia (hereafter referred to as Eth19 and Eth20, respectively).\u003c/p\u003e \u003cp\u003eIn all environments, the lines were grown in a randomized complete block design with two replications. In Kenya, each line was grown in two 0.7 m rows spaced 0.3 m apart. The susceptible check CItr 6255 (\u0026lsquo;Red Bobs\u0026rsquo;) was planted after every fifty entries. Twin spreader rows consisting of susceptible varieties \u0026lsquo;Cacucke\u0026rsquo;, \u0026lsquo;Eagle 10\u0026rsquo;, \u0026lsquo;Robin\u0026rsquo;, and six CIMMYT experimental lines were grown alongside and between each plot. Additionally, a border of spreader rows encompassed the nursery. Starting at the booting growth stage through heading (Zadok\u0026rsquo;s growth stage 37\u0026ndash;60), spreader rows were inoculated with a bulk inoculum of locally collected urediniospores. \u003cem\u003ePgt\u003c/em\u003e races TTKSK, TTKST, TTKTK, and TTKTT were among the races present in the Kenya nurseries (Mago et al. \u003cspan citationid=\"CR43\" class=\"CitationRef\"\u003e2022\u003c/span\u003e). In Ethiopia, lines were grown in 1.0 m double rows flanked by spreader rows consisting of susceptible varieties \u0026lsquo;Arendeto\u0026rsquo;, \u0026lsquo;Digalu\u0026rsquo;, \u0026lsquo;Local Red\u0026rsquo;, \u0026lsquo;Morocco\u0026rsquo;, and \u0026lsquo;PBW343\u0026rsquo; (Kosgey et al. \u003cspan citationid=\"CR34\" class=\"CitationRef\"\u003e2021\u003c/span\u003e). Spreader rows were inoculated at Zadok\u0026rsquo;s growth stages 37\u0026ndash;60 using bulk of locally collected urediniospores.\u003c/p\u003e \u003cp\u003eStem rust severity was assessed once spreaders and the check line CItr 6255 reached 50% severity in both Kenya and Ethiopia. Stem rust severity was scored on a whole plot basis using the modified Cobb scale of 0 to 100%, where 0 is equal to complete immunity or absence of infection symptoms and 100% is equal to complete susceptibility (Peterson et al., \u003cspan citationid=\"CR63\" class=\"CitationRef\"\u003e1948\u003c/span\u003e). Infection response (IR) was recorded based on pustule size and prevalence of chlorosis and necrosis. Four classifications were assigned for IR: S, susceptible; MS, moderately susceptible; MR, moderately resistant; and R, resistant (Roelfs et al. \u003cspan citationid=\"CR72\" class=\"CitationRef\"\u003e1992\u003c/span\u003e). Lines with heterogeneous infection responses were denoted using a hyphen (e.g. MS-MR) with the most prevalent IR category listed first. Linearized infection response (LIR) and coefficient of infection (CI) were calculated from the field observations using a custom R script (Gao et al. \u003cspan citationid=\"CR16\" class=\"CitationRef\"\u003e2016\u003c/span\u003e). LIR is infection response scaled from 0 to 1.0. CI was then calculated as the product of LIR and severity (Stubbs et al. \u003cspan citationid=\"CR83\" class=\"CitationRef\"\u003e1986\u003c/span\u003e). Two observations were recorded in both years in Kenya, while three observations were made in Ethiopia. However, only the CI value from the final observation in each environment was considered during data analysis and QTL mapping.\u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec6\" class=\"Section2\"\u003e \u003ch2\u003eGenotyping\u003c/h2\u003e \u003cp\u003eSingle nucleotide polymorphism markers (SNPs) were called using genotyping-by-sequencing (GBS; Elshire et al. \u003cspan citationid=\"CR13\" class=\"CitationRef\"\u003e2011\u003c/span\u003e). DNA was extracted from juvenile leaves collected from the parents and BC\u003csub\u003e1\u003c/sub\u003eF\u003csub\u003e5:6\u003c/sub\u003e RILs using BioSprint 96 DNA Plant Kits (QIAGEN, Hilden, Germany). Extracted DNA was quantified using Quant-iT PicoGreen dsDNA assay kits (Invitrogen, Waltham, Massachusetts, USA) and normalized to 20ng/\u0026micro;L. GBS libraries were constructed using a two-enzyme approach as described by Poland et al. (\u003cspan citationid=\"CR64\" class=\"CitationRef\"\u003e2012\u003c/span\u003e) with two modifications: 1) two unique barcodes were ligated to each sample in an effort to reduce sequencing bias of certain barcodes and 2) barcode and common adapter concentrations were increased to 0.1 \u0026micro;M and 50 \u0026micro;M, respectively. Two 64-plex libraries were generated including four biological replicates of parent MN07098-6 and three biological replicates of parent CItr 11390. Each library was sequenced on one Illumina HiSeq 2000 lane to produce single-end 100 bp reads.\u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec7\" class=\"Section2\"\u003e \u003ch2\u003eMarker calling, filtering, imputation, and distortion\u003c/h2\u003e \u003cp\u003eRaw sequencing reads were demultiplexed using Sabre (\u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://github.com/najoshi/sabre\u003c/span\u003e\u003cspan address=\"https://github.com/najoshi/sabre\" targettype=\"URL\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e). Adapters and low quality bases (Q\u0026thinsp;\u0026lt;\u0026thinsp;30) were trimmed using Cutadapt (Martin \u003cspan citationid=\"CR46\" class=\"CitationRef\"\u003e2011\u003c/span\u003e). Trimmed reads were aligned to the \u003cem\u003eT. aestivum\u003c/em\u003e cv. Chinese Spring reference assembly, RefSeq v2.1 (Zhu et al. \u003cspan citationid=\"CR96\" class=\"CitationRef\"\u003e2021\u003c/span\u003e) using the BWA-MEM algorithm with default parameters (Li 2013). Aligned reads were filtered (MAPQ\u0026thinsp;\u0026lt;\u0026thinsp;40) and sorted with SAMtools, and SNP calling was done with BCFtools (Li \u003cspan citationid=\"CR41\" class=\"CitationRef\"\u003e2011\u003c/span\u003e; Li et al. \u003cspan citationid=\"CR42\" class=\"CitationRef\"\u003e2009\u003c/span\u003e). Genotypes within a site were set to missing data based on read depth (DP\u0026thinsp;\u0026lt;\u0026thinsp;3). Sites with \u0026ge;\u0026thinsp;20% missing data and QUAL\u0026thinsp;\u0026lt;\u0026thinsp;20 were then discarded as well as monomorphic and multiallelic sites. Sites were further filtered by removing those with a minor allele frequency (MAF)\u0026thinsp;\u0026lt;\u0026thinsp;40% or heterozygosity\u0026thinsp;\u0026gt;\u0026thinsp;10%. Missing marker data were imputed in TASSEL 5 using LD-kNNi, an imputation method based on the \u003cem\u003ek\u003c/em\u003e-nearest neighbor algorithm (Bradbury et al. \u003cspan citationid=\"CR7\" class=\"CitationRef\"\u003e2007\u003c/span\u003e; Money et al. \u003cspan citationid=\"CR52\" class=\"CitationRef\"\u003e2015\u003c/span\u003e). Genotypes were converted from a nucleotide-based format to a parent-based format, and consensus genotypes were called for both parents from the replicated samples using the GenosToABHPlugin in TASSEL 5 (Bradbury et al. \u003cspan citationid=\"CR7\" class=\"CitationRef\"\u003e2007\u003c/span\u003e).\u003c/p\u003e \u003cp\u003eFurther marker analysis and filtering was performed in R using the \u0026lsquo;qtl\u0026rsquo; and \u0026lsquo;ASMap\u0026rsquo; packages (Broman et al. \u003cspan citationid=\"CR10\" class=\"CitationRef\"\u003e2003\u003c/span\u003e; R Core Team 2022; Taylor and Butler \u003cspan citationid=\"CR86\" class=\"CitationRef\"\u003e2017\u003c/span\u003e). Redundant markers with identical genotype patterns across all individuals were identified, and all but one representative marker were removed from the data set. Using the filtered dataset, all pairwise marker similarities were calculated among the individuals. Individuals with identical genotypes at \u0026ge;\u0026thinsp;95% markers were identified, and the line with a greater amount of missing data among the pair was removed from analysis. Lastly, markers were tested for segregation distortion.\u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec8\" class=\"Section2\"\u003e \u003ch2\u003eMolecular marker screening\u003c/h2\u003e \u003cp\u003eThe parents and population were screened with a series of markers previously described as linked to resistance loci: \u003cem\u003eSr7a\u003c/em\u003e, \u003cem\u003eSr11\u003c/em\u003e, \u003cem\u003eSr12\u003c/em\u003e, \u003cem\u003eSr13\u003c/em\u003e, \u003cem\u003eSr21\u003c/em\u003e, \u003cem\u003eSr38/Lr37/Yr17\u003c/em\u003e, \u003cem\u003eSr63\u003c/em\u003e, \u003cem\u003eSr8155B1\u003c/em\u003e, and \u003cem\u003eSrTm4\u003c/em\u003e (Table\u0026nbsp;\u003cspan refid=\"Tab1\" class=\"InternalRef\"\u003e1\u003c/span\u003e). Additionally, shortly after population development concluded, MN07098-6 was found to be heterogeneous for \u003cem\u003eSr57/Lr34/Yr18\u003c/em\u003e on chromosome 7D. The plant used for the backcross was not the selfed progeny of the original MN07098-6 crossing parent. Therefore, the lines used for the initial and backcross were potentially near-isogenic for \u003cem\u003eSr57/Lr34/Yr18\u003c/em\u003e. The population was screened for \u003cem\u003eSr57/Lr34/Yr18\u003c/em\u003e to determine if the population was segregating at this locus using a perfect marker (Table\u0026nbsp;\u003cspan refid=\"Tab1\" class=\"InternalRef\"\u003e1\u003c/span\u003e).\u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec9\" class=\"Section2\"\u003e \u003ch2\u003eAneuploidy and chromosomal deletion detection\u003c/h2\u003e \u003cp\u003eIn silico detection of aneuploidy and chromosomal deletions among the progeny lines was conducted using the GBS data and a method adapted from Singh et al. (\u003cspan citationid=\"CR77\" class=\"CitationRef\"\u003e2020\u003c/span\u003e). MN07098-6 and CItr 11390 were used as euploid allohexaploid and allotetraploid controls, respectively. Each chromosome was divided in 100 Mb bins with a sliding window step size of 50 Mb. Bins\u0026thinsp;\u0026lt;\u0026thinsp;50 Mb in length at the end of chromosomes were discarded as they nested completely within another bin. BEDtools multicov (Quinlan and Hall \u003cspan citationid=\"CR68\" class=\"CitationRef\"\u003e2010\u003c/span\u003e) was used to count the number of reads uniquely aligned within each bin. Counts were normalized across the genome for all bins and individuals to remove bias due to differential sequencing. Karyotypes for each sample were plotted using base R functions and the \u0026lsquo;chromPlot\u0026rsquo; package (Or\u0026oacute;stica and Verdugo \u003cspan citationid=\"CR61\" class=\"CitationRef\"\u003e2016\u003c/span\u003e; R Core Team 2022). Bins were assigned a copy number based on normalized values. Due to artifacts in the assembly leading to a small number of reads aligning to D genome chromosomes for CItr 11390, the mean normalized value across all bins deviated slightly from one. Therefore, the upper bin value limit for each copy number was determined using the following equation:\u003cdiv id=\"Equa\" class=\"Equation\"\u003e\u003cdiv format=\"TEX\" class=\"mathdisplay\" id=\"FileID_Equa\" name=\"EquationSource\"\u003e\n$${bin}_{max }=\\frac{copy number* \\mu }{2}+2\\sigma$$\u003c/div\u003e\u003c/div\u003e\u003c/p\u003e \u003cp\u003ewhere \u0026micro; and σ are the mean normalized values and standard deviation across all bins, respectively. Karyotype plots were visually assessed to detect chromosomal deletions and aneuploidy. Individuals with signatures for any form of aneuploidy involving chromosomes of either the A or B genome were not included in genetic and phenotypic analysis. Additionally, RILs with all forms of aneuploidy or deletions other than monosomy in the D genome were excluding from further analysis. RILs with monosomy affecting a single D genome chromosome were included in initial data analysis as there is representation of all 21 chromosomes.\u003c/p\u003e \u003cdiv id=\"Sec10\" class=\"Section3\"\u003e \u003ch2\u003eLinkage map construction\u003c/h2\u003e \u003cp\u003eLinkage map construction was performed in R using the MSTmap algorithm implemented in the package \u0026lsquo;ASMap\u0026rsquo; (R Core Team 2022; Taylor and Butler \u003cspan citationid=\"CR86\" class=\"CitationRef\"\u003e2017\u003c/span\u003e; Wu et al. \u003cspan citationid=\"CR90\" class=\"CitationRef\"\u003e2008\u003c/span\u003e). Due to the population being developed from a single BC\u003csub\u003e1\u003c/sub\u003eF\u003csub\u003e1\u003c/sub\u003e plant, it was treated as a self-pollinated RIL population despite the generation of backcrossing. Linkage groups were independently constructed for each chromosome based on the physical positions assigned to the GBS markers. Map construction and marker ordering was conducted using the Kosambi mapping function and a significance threshold for marker clustering of \u003cem\u003ep\u003c/em\u003e\u0026thinsp;=\u0026thinsp;10\u003csup\u003e\u0026minus;\u0026thinsp;6\u003c/sup\u003e (Kosambi \u003cspan citationid=\"CR33\" class=\"CitationRef\"\u003e1943\u003c/span\u003e). Final marker order was set by jittering the positions of markers mapped to the same position.\u003c/p\u003e \u003c/div\u003e \u003c/div\u003e \u003cdiv id=\"Sec11\" class=\"Section2\"\u003e \u003ch2\u003eData analysis\u003c/h2\u003e \u003cdiv id=\"Sec12\" class=\"Section3\"\u003e \u003ch2\u003eSeedling-stage stem rust resistance\u003c/h2\u003e \u003cp\u003eMean values calculated from the converted infection type scores across the four replications for \u003cem\u003ePgt\u003c/em\u003e races observed to have variable virulence to the parents were used for analysis and QTL mapping. Lines with a linearized score of \u0026lt;\u0026thinsp;6 were considered resistant whereas those with scores\u0026thinsp;\u0026ge;\u0026thinsp;6 were considered susceptible. A chi-square (χ\u003csup\u003e2\u003c/sup\u003e) test was conducted to determine the genetic architecture of resistance to the screened \u003cem\u003ePgt\u003c/em\u003e races within the population.\u003c/p\u003e \u003c/div\u003e \u003c/div\u003e \u003cdiv id=\"Sec13\" class=\"Section2\"\u003e \u003ch2\u003eField-based stem rust\u003c/h2\u003e \u003cp\u003eThe population was determined to be segregating for \u003cem\u003eSr57/Lr34/Yr18\u003c/em\u003e. To determine the significance of the effects of \u003cem\u003eSr57/Lr34/Yr18\u003c/em\u003e and ploidy on stem rust CI, mixed linear models were fit in R using the \u0026lsquo;\u003cem\u003elme4\u003c/em\u003e\u0026rsquo; package (Bates et al. \u003cspan citationid=\"CR5\" class=\"CitationRef\"\u003e2015\u003c/span\u003e). \u003cem\u003eSr57/Lr34/Yr18\u003c/em\u003e, ploidy state, and environments treated as fixed effects. Genotypes, genotype-environment interactions, and blocks within environments were treated as random effects. Backwards elimination of all effects in the across-environment model was performed using the \u0026lsquo;\u003cem\u003elmerTest\u003c/em\u003e\u0026rsquo; package in R (Kuznetsova et al. 2017). Significance values for the fixed and random effects were calculated based on Sattethwaite\u0026rsquo;s approximation and likelihood ratio tests, respectively.\u003c/p\u003e \u003cp\u003eThe interaction between genotype and environment was determined to be significant, which warranted the separate analysis of each environment. CI values within environments were then adjusted for the effect of \u003cem\u003eSr57/Lr34/Yr18\u003c/em\u003e by fitting mixed linear models for each of the four environments separately. \u003cem\u003eSr57\u003c/em\u003e/\u003cem\u003eLr34\u003c/em\u003e/\u003cem\u003eYr18\u003c/em\u003e was significant in the across-environment model, while ploidy was determined to be insignificant. Therefore, \u003cem\u003eSr57\u003c/em\u003e/\u003cem\u003eLr34\u003c/em\u003e/\u003cem\u003eYr18\u003c/em\u003e was the only fixed effect retained in the individual environment models. Genotype and block were included as random effects.\u003c/p\u003e \u003cp\u003eIn addition to visual inspection of diagnostic plots for each single environment model, residuals were tested for normality and homogeneity with Sharpiro-Wilk and Levene\u0026rsquo;s tests, respectively. If necessary, data transformation was performed if assumptions of normality and residual homogeneity were violated. Log and square root transformations were tested for each environment. Additionally, an arcsine transformation of the Eth20 data was tested after dividing the data by 100. The model selection approach above was repeated to assess the effectiveness of each transformation. Coefficients from the selected models were extracted to correct phenotypic values for the effect of \u003cem\u003eSr57/Lr34/Yr18\u003c/em\u003e. Pearson correlation coefficients between environments were calculated in R using the adjusted datasets.\u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec14\" class=\"Section2\"\u003e \u003ch2\u003eQTL Mapping\u003c/h2\u003e \u003cp\u003eQTL detection within individual environments was conducted using multiple QTL mapping (MQM) in R using the \u0026lsquo;\u003cem\u003eqtl\u003c/em\u003e\u0026rsquo; package (Broman et al. \u003cspan citationid=\"CR10\" class=\"CitationRef\"\u003e2003\u003c/span\u003e). Multiple QTL mapping was conducted using an automated forward/backward analysis, and the final models were selected based on penalized LOD score (pLOD; Manichaikul et al. \u003cspan citationid=\"CR44\" class=\"CitationRef\"\u003e2009\u003c/span\u003e). First, a single locus model was fitted for each environment and race with segregating resistance using multiple imputation and the scanone function. A 95% significance threshold for each single locus scan was determined by 1000 permutations across the entire genome. Significant markers (α\u0026thinsp;=\u0026thinsp;0.05) from the single locus models were used as the initial set of cofactors for MQM analysis. Next, a two-dimensional genome scan was performed with 1000 permutations to calculate the threshold for significant LOD scores (α\u0026thinsp;=\u0026thinsp;0.05) described by Broman et al. (\u003cspan citationid=\"CR10\" class=\"CitationRef\"\u003e2003\u003c/span\u003e) using the scantwo function. These thresholds were then used to calculate the pLOD scores to be used for multiple QTL mapping. Using the loci detected by single-QTL scans and the calculated pLOD thresholds, Multiple QTL models with maximal pLOD scores were selected and refined using the stepwiseqtl command and multiple imputation. Within the stepwiseqtl function, forward selection was allowed to proceed with a maximum of 10 QTL. The reported LOD support intervals, reported as 95% Bayes credible intervals, were calculated using the command bayesint.\u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec15\" class=\"Section2\"\u003e \u003ch2\u003eAssessment of QTL interactions and combinations\u003c/h2\u003e \u003cp\u003eAll pairwise QTL x QTL interactions were simultaneously tested during multiple QTL model selection using stepwiseqtl, and models were refined in the event of a significant interaction between QTL. Significance and percent of phenotypic variance explained of the final models including interactions was determined in R using the fitql command from the \u0026lsquo;qtl\u0026rsquo; package.\u003c/p\u003e \u003cp\u003eMean CI values of RILs carrying all combinations of QTL detected in multiple environments, as well as \u003cem\u003eSr57\u003c/em\u003e/\u003cem\u003eLr34\u003c/em\u003e/\u003cem\u003eYr18\u003c/em\u003e, were compared within each environments. Mean CI values of each class of RILs were calculated based on the allele at the peak marker of the QTL using the unadjusted and untransformed data from each environment. Tukey\u0026rsquo;s HSD tests were performed to determine whether the means of each class of RIL were significantly different within each environment.\u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec16\" class=\"Section2\"\u003e \u003ch2\u003eMolecular marker development\u003c/h2\u003e \u003cp\u003ePCR Allele Competitive Extension (PACE\u0026trade;) SNP genotyping assays (3CR Bioscience, Harlow, UK) were developed for novel QTL detected in multiple environments. Primers were designed at GBS SNPs throughout the LOD intervals. Primers targeting resistance alleles were labelled with the HEX fluorophore, and those for susceptible alleles were labelled with the FAM fluorophore. Each PACE\u0026trade; assay mix contained 60 \u0026micro;L HEX-labeled primer (100 \u0026micro;M), 60 \u0026micro;L FAM-labeled primer (100 \u0026micro;M), 120 \u0026micro;L common primer (100 \u0026micro;M), and 230 \u0026micro;L molecular-grade water (Sigma Cat. No. 4502). For each assay, genotyping was conducted using standard conditions for PACE. 8 \u0026micro;L reactions consisted of 2 \u0026micro;L template DNA (20 ng/\u0026micro;L), 0.11 \u0026micro;L assay mix, 4 \u0026micro;L 2x PACE\u0026trade; master mix, and 1.89 \u0026micro;L molecular grade water. Genotyping was performed on a Roche Lightcycler\u0026reg; 480 (Roche Diagnostics, Indianapolis, IN, USA) using the following reaction conditions: hot-start activation at 95\u0026deg;C for 15 minutes, 10 touchdown cycles at 94\u0026deg;C for 20 seconds, touchdown at 61\u0026deg;C for 1 minute with a 0.6\u0026deg;C reduction per cycle until a secondary target of 55\u0026deg;C, and 48 cycles at 94\u0026deg;C for 20 seconds followed by 55\u0026deg;C for 1 minute/cycle. An endpoint read was taken following a 2-minute cooldown cycle at 37\u0026deg;C. Alleles were called using the Roche Lightcycler\u0026reg; 480 Endpoint Genotyping Software Module.\u003c/p\u003e \u003c/div\u003e"},{"header":"Results","content":"\u003cdiv id=\"Sec18\" class=\"Section2\"\u003e \u003ch2\u003eGenotyping and Aneuploid Detection\u003c/h2\u003e \u003cp\u003eGenotyping-by-sequencing (GBS) resulted in a total of 269,278,521 uniquely mapped reads across all samples. The distribution of read counts per sample was relatively normally distributed, and the median number of reads per sample was ~\u0026thinsp;2.2\u0026nbsp;million (Fig. \u003cspan refid=\"MOESM1\" class=\"InternalRef\"\u003eS1\u003c/span\u003eA). Total read counts were first normalized across all samples. Counts in each bin were then normalized by the median count for each bin across all samples. As expected, the distribution of normalized read counts per bin had a normal distribution (Fig. \u003cspan refid=\"MOESM1\" class=\"InternalRef\"\u003eS1\u003c/span\u003eB). Small shoulder peaks were observed in the distribution centered around 0, 0.5, and 2.0, which represent deleted, monosomic, and duplicated bins, respectively. Over 98% of bins had a normalized read count around 1.0, representative of the disomic state.\u003c/p\u003e \u003cp\u003eThe analysis accurately assigned a copy number of 0 to all D genome chromosomes of the tetraploid parent, CItr 11390, while all 21 chromosomes in MN07098-6 were determined to be disomic and fully intact as expected (Fig.\u0026nbsp;\u003cspan refid=\"Fig1\" class=\"InternalRef\"\u003e1\u003c/span\u003e). Of the 121 RILs, 30 were determined to be aneuploids (Fig. \u003cspan refid=\"MOESM2\" class=\"InternalRef\"\u003eS2\u003c/span\u003e). Of the 30 aneuploid RILs, 24 had an abnormality with chromosomes in homoeologous groups 4 or 6, with 18 possessing chromosome 6D abnormalities (Table \u003cspan refid=\"MOESM1\" class=\"InternalRef\"\u003eS1\u003c/span\u003e). The most common form of aneuploidy was monosomy (15 out of 30 lines). Notes made in observation nurseries in 2019 in Saint Paul, Minnesota and Crookston, Minnesota corroborate the \u003cem\u003ein silico\u003c/em\u003e analysis. Several of the RILs identified as aneuploids were noted in both locations to have poor vigor, short stature, and/or late maturity. The 30 RILs identified as aneuploids, as well as one RIL of a pair with \u0026gt;\u0026thinsp;95% identical SNPs, were initially retained during variant calling. However, the 16 RILs with forms of aneuploidy or deletions involving the A or B genome, as well as 1 euploid RIL of the pair found to have \u0026gt;\u0026thinsp;95% identical SNPs, were removed during SNP filtering. The remaining 14 aneuploid RILs (11 M6D and 3 M4D) were included in phenotypic data analysis to determine whether monosomy involving either chromosome had a significant effect on stem rust infection.\u003c/p\u003e \u003cp\u003e \u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec19\" class=\"Section2\"\u003e \u003ch2\u003eLinkage map construction\u003c/h2\u003e \u003cp\u003eAfter filtering, 2,287 SNPs were obtained from GBS across the A and B genomes. No SNPs in the D genome remained after filtering as expected. In addition to the GBS SNPs, markers linked to previously mapped stem rust resistance genes were used for linkage mapping. Of the molecular markers assayed, those linked to \u003cem\u003eSr7a\u003c/em\u003e, \u003cem\u003eSr11\u003c/em\u003e, \u003cem\u003eSr63\u003c/em\u003e, and \u003cem\u003eSr8155B1\u003c/em\u003e were polymorphic between the parents and added to the GBS marker dataset (Table\u0026nbsp;\u003cspan refid=\"Tab1\" class=\"InternalRef\"\u003e1\u003c/span\u003e). Additionally, the population was segregating for \u003cem\u003eSr57/Lr34/Yr18\u003c/em\u003e. However, as no GBS SNPs were retained on chromosome 7D after filtering, \u003cem\u003eSr57/Lr34/Yr18\u003c/em\u003e was not included in linkage mapping.\u003c/p\u003e \u003cp\u003e \u003cdiv class=\"gridtable\"\u003e\u003ctable float=\"Yes\" id=\"Tab1\" border=\"1\"\u003e \u003ccaption language=\"En\"\u003e \u003cdiv class=\"CaptionNumber\"\u003eTable 1\u003c/div\u003e \u003cdiv class=\"CaptionContent\"\u003e \u003cp\u003eMolecular marker assays for known stem rust resistance genes\u003c/p\u003e \u003c/div\u003e \u003c/caption\u003e \u003ccolgroup cols=\"6\"\u003e \u003cdiv align=\"left\" class=\"colspec\" colname=\"c1\" colnum=\"1\"\u003e\u003c/div\u003e \u003cdiv align=\"left\" class=\"colspec\" colname=\"c2\" colnum=\"2\"\u003e\u003c/div\u003e \u003cdiv align=\"left\" class=\"colspec\" colname=\"c3\" colnum=\"3\"\u003e\u003c/div\u003e \u003cdiv align=\"left\" class=\"colspec\" colname=\"c4\" colnum=\"4\"\u003e\u003c/div\u003e \u003cdiv align=\"left\" class=\"colspec\" colname=\"c5\" colnum=\"5\"\u003e\u003c/div\u003e \u003cdiv align=\"left\" class=\"colspec\" colname=\"c6\" colnum=\"6\"\u003e\u003c/div\u003e \u003cthead\u003e \u003ctr\u003e \u003cth align=\"left\" colname=\"c1\"\u003e\u0026nbsp;\u003c/th\u003e \u003cth align=\"left\" colname=\"c2\"\u003e\u0026nbsp;\u003c/th\u003e \u003cth align=\"left\" colname=\"c3\"\u003e\u0026nbsp;\u003c/th\u003e \u003cth align=\"left\" colspan=\"2\" nameend=\"c5\" namest=\"c4\"\u003e \u003cp\u003eParent allele\u003c/p\u003e \u003c/th\u003e \u003cth align=\"left\" colname=\"c6\"\u003e\u0026nbsp;\u003c/th\u003e \u003c/tr\u003e \u003c/thead\u003e \u003ctbody\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003eGene\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003eChromosome\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c3\"\u003e \u003cp\u003eMarker\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c4\"\u003e \u003cp\u003eMN07098-6\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c5\"\u003e \u003cp\u003eCItr 11390\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c6\"\u003e \u003cp\u003eReference\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003e\u003cem\u003eSr7a\u003c/em\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003e4AL\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c3\"\u003e \u003cp\u003eKASP_IWB12146\u003c/p\u003e \u003cp\u003eKASP_IWB47019\u003c/p\u003e \u003cp\u003eKASP_IWB68386\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c4\"\u003e \u003cp\u003eG\u003csup\u003eb\u003c/sup\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c5\"\u003e \u003cp\u003eA\u003csup\u003ea\u003c/sup\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c6\"\u003e \u003cp\u003eBajgain et al. \u003cspan citationid=\"CR1\" class=\"CitationRef\"\u003e2015a\u003c/span\u003e;\u003c/p\u003e \u003cp\u003eOlivera et al. \u003cspan citationid=\"CR58\" class=\"CitationRef\"\u003e2022\u003c/span\u003e\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003e\u003cem\u003eSr11\u003c/em\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003e6BL\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c3\"\u003e \u003cp\u003eKASP_IWB72471\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c4\"\u003e \u003cp\u003eT\u003csup\u003ea\u003c/sup\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c5\"\u003e \u003cp\u003eG\u003csup\u003eb\u003c/sup\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c6\"\u003e \u003cp\u003eNirmala et al. \u003cspan citationid=\"CR54\" class=\"CitationRef\"\u003e2016\u003c/span\u003e\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003e\u003cem\u003eSr12\u003c/em\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003e3BL\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c3\"\u003e \u003cp\u003eNB-LRR3\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c4\"\u003e \u003cp\u003eG\u003csup\u003eb\u003c/sup\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c5\"\u003e \u003cp\u003eG\u003csup\u003eb\u003c/sup\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c6\"\u003e \u003cp\u003eHiebert et al. \u003cspan citationid=\"CR25\" class=\"CitationRef\"\u003e2016\u003c/span\u003e\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003e\u003cem\u003eSr13\u003c/em\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003e6AL\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c3\"\u003e \u003cp\u003eSr13R/S_SNP\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c4\"\u003e \u003cp\u003eA\u003csup\u003ea\u003c/sup\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c5\"\u003e \u003cp\u003eNull\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c6\"\u003e \u003cp\u003emaswheat.ucdavis.edu/protocols/Sr13\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003e\u003cem\u003eSr21\u003c/em\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003e2AL\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c3\"\u003e \u003cp\u003eSr21TRYF5R5\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c4\"\u003e \u003cp\u003eNull\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c5\"\u003e \u003cp\u003eNull\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c6\"\u003e \u003cp\u003eChen et al. \u003cspan citationid=\"CR11\" class=\"CitationRef\"\u003e2018\u003c/span\u003e\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003e\u003cem\u003eSr38/Lr37/Yr17\u003c/em\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003e2NS:2AS\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c3\"\u003e \u003cp\u003eLr37-Yr17-Sr38_GBG-KASP\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c4\"\u003e \u003cp\u003eA\u003csup\u003ea\u003c/sup\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c5\"\u003e \u003cp\u003eA\u003csup\u003ea\u003c/sup\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c6\"\u003e \u003cp\u003eHelguera et al. \u003cspan citationid=\"CR23\" class=\"CitationRef\"\u003e2003\u003c/span\u003e;\u003c/p\u003e \u003cp\u003eMustahsan \u003cspan citationid=\"CR53\" class=\"CitationRef\"\u003e2022\u003c/span\u003e\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003e\u003cem\u003eSr57/Lr34/Yr18\u003c/em\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003e7DS\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c3\"\u003e \u003cp\u003ewMAS000003\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c4\"\u003e \u003cp\u003eT/A\u003csup\u003ec\u003c/sup\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c5\"\u003e \u003cp\u003eA\u003csup\u003ea\u003c/sup\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c6\"\u003e \u003cp\u003emaswheat.ucdavis.edu/protocols/Lr34\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003e\u003cem\u003eSr63\u003c/em\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003e2AL\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c3\"\u003e \u003cp\u003eKASP_IWB32429\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c4\"\u003e \u003cp\u003eT\u003csup\u003ea\u003c/sup\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c5\"\u003e \u003cp\u003eG\u003csup\u003eb\u003c/sup\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c6\"\u003e \u003cp\u003eMago et al. \u003cspan citationid=\"CR43\" class=\"CitationRef\"\u003e2022\u003c/span\u003e\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003e\u003cem\u003eSr8155B1\u003c/em\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003e6AS\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c3\"\u003e \u003cp\u003eKASP_IWB10558\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c4\"\u003e \u003cp\u003eA\u003csup\u003ea\u003c/sup\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c5\"\u003e \u003cp\u003eG\u003csup\u003eb\u003c/sup\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c6\"\u003e \u003cp\u003eNirmala et al. \u003cspan citationid=\"CR55\" class=\"CitationRef\"\u003e2017\u003c/span\u003e\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003e\u003cem\u003eSrTm4\u003c/em\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003e2AL\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c3\"\u003e \u003cp\u003eBQ421276\u003c/p\u003e \u003cp\u003eDR732348\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c4\"\u003e \u003cp\u003eNull\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c5\"\u003e \u003cp\u003eNull\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c6\"\u003e \u003cp\u003eBriggs et al. \u003cspan citationid=\"CR9\" class=\"CitationRef\"\u003e2015\u003c/span\u003e\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003c/tbody\u003e \u003c/colgroup\u003e \u003ctfoot\u003e \u003ctr\u003e\u003ctd colspan=\"6\"\u003e\u003csup\u003ea\u003c/sup\u003e Reportedly predictive of susceptibility\u003c/td\u003e\u003c/tr\u003e \u003ctr\u003e\u003ctd colspan=\"6\"\u003e\u003csup\u003eb\u003c/sup\u003e Reportedly predictive of resistance\u003c/td\u003e\u003c/tr\u003e \u003ctr\u003e\u003ctd colspan=\"6\"\u003e\u003csup\u003ec\u003c/sup\u003e MN07098-6 is heterogeneous for \u003cem\u003eSr57\u003c/em\u003e/\u003cem\u003eLr34\u003c/em\u003e/\u003cem\u003eYr18\u003c/em\u003e\u003c/td\u003e\u003c/tr\u003e \u003c/tfoot\u003e \u003c/table\u003e\u003c/div\u003e \u003c/p\u003e \u003cp\u003eThe 2,291 GBS and molecular markers used for linkage mapping were assigned to 21 linkage groups representing the 14 chromosomes of the A and B genomes. Chromosomes 1B, 2A, 2B, 3B, 5A, 5B, and 6A were represented by a single linkage group, while the other seven chromosomes were split into two linkage groups each. Following marker ordering, markers with identical genotypes were identified and all but one representative marker for each group of markers were removed from the map. The final map consisted of 1,212 markers with a total length of 1,230.7 cM and an average spacing of ~\u0026thinsp;1.0 cM between markers.\u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec20\" class=\"Section2\"\u003e \u003ch2\u003eSeedling-stage disease evaluation\u003c/h2\u003e \u003cp\u003eParental infection types to the four \u003cem\u003ePgt\u003c/em\u003e races screened are shown in Table\u0026nbsp;\u003cspan refid=\"Tab2\" class=\"InternalRef\"\u003e2\u003c/span\u003e. MN07098-6 and CItr 11390 were both susceptible at the seedling stage (IT\u0026thinsp;=\u0026thinsp;33\u0026thinsp;+\u0026thinsp;c to 3+) to the two Ug99 lineage races screened (TTKSK and TTKTT). The results of the TTKTT screen were unexpected, as CItr 11390 was found to carry the allele at KASP_IWB10558 reported to be predictive of resistance provided by \u003cem\u003eSr8155B1\u003c/em\u003e (Nirmala et al. \u003cspan citationid=\"CR55\" class=\"CitationRef\"\u003e2017\u003c/span\u003e). Note that Nirmala et al. (\u003cspan citationid=\"CR55\" class=\"CitationRef\"\u003e2017\u003c/span\u003e) reported the same isolate of TTKTT used in this study, 14KEN58-1, to be avirulent to \u003cem\u003eSr8155B1\u003c/em\u003e with infection type \u0026lsquo;0;\u0026rsquo;. However, the same study reported the 04KEN156/04 isolate to be virulent to \u003cem\u003eSr8155B1\u003c/em\u003e with infection type \u0026lsquo;3+\u0026rsquo;.\u003c/p\u003e \u003cp\u003eInitial inspection of the data indicated CItr 11390 was moderately resistant to TTRTF (IT\u0026thinsp;=\u0026thinsp;2). However, analysis of the RILs found only four resistant lines (mean linearized IT\u0026thinsp;\u0026lt;\u0026thinsp;6) to 100 susceptible lines (mean linearized IT\u0026thinsp;\u0026ge;\u0026thinsp;6). Additionally, the upper bound of the 95% confidence interval for all four RILs with a mean linearized IT\u0026thinsp;\u0026lt;\u0026thinsp;6, was greater than the threshold for resistance (Fig. S4). Therefore, TTRTF was not included in QTL mapping.\u003c/p\u003e \u003cp\u003e \u003cdiv class=\"gridtable\"\u003e\u003ctable float=\"Yes\" id=\"Tab2\" border=\"1\"\u003e \u003ccaption language=\"En\"\u003e \u003cdiv class=\"CaptionNumber\"\u003eTable 2\u003c/div\u003e \u003cdiv class=\"CaptionContent\"\u003e \u003cp\u003eIsolate designation, origin, parental infection type, and virulence phenotype of \u003cem\u003ePuccinia graminis\u003c/em\u003e f. sp. \u003cem\u003etritici\u003c/em\u003e races used to evaluate seedling resistance in the population.\u003c/p\u003e \u003c/div\u003e \u003c/caption\u003e \u003ccolgroup cols=\"5\"\u003e \u003cdiv align=\"left\" class=\"colspec\" colname=\"c1\" colnum=\"1\"\u003e\u003c/div\u003e \u003cdiv align=\"left\" class=\"colspec\" colname=\"c2\" colnum=\"2\"\u003e\u003c/div\u003e \u003cdiv align=\"left\" class=\"colspec\" colname=\"c3\" colnum=\"3\"\u003e\u003c/div\u003e \u003cdiv align=\"left\" class=\"colspec\" colname=\"c4\" colnum=\"4\"\u003e\u003c/div\u003e \u003cdiv align=\"left\" class=\"colspec\" colname=\"c5\" colnum=\"5\"\u003e\u003c/div\u003e \u003cthead\u003e \u003ctr\u003e \u003cth align=\"left\" colname=\"c1\"\u003e \u003cp\u003eRace\u003c/p\u003e \u003c/th\u003e \u003cth align=\"left\" colname=\"c2\"\u003e \u003cp\u003eIsolate\u003c/p\u003e \u003c/th\u003e \u003cth align=\"left\" colname=\"c3\"\u003e \u003cp\u003eOrigin\u003c/p\u003e \u003c/th\u003e \u003cth align=\"left\" colname=\"c4\"\u003e \u003cp\u003eMN07098-6\u003c/p\u003e \u003c/th\u003e \u003cth align=\"left\" colname=\"c5\"\u003e \u003cp\u003eCItr 11390\u003c/p\u003e \u003c/th\u003e \u003c/tr\u003e \u003c/thead\u003e \u003ctbody\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003eTTKSK\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003e04KEN156/04\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c3\"\u003e \u003cp\u003eKenya\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c4\"\u003e \u003cp\u003e3+\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c5\"\u003e \u003cp\u003e3+\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003eTTKTT\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003e14KEN58-1\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c3\"\u003e \u003cp\u003eKenya\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c4\"\u003e \u003cp\u003e3+\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c5\"\u003e \u003cp\u003e33\u0026thinsp;+\u0026thinsp;c\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003eTKTTF\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003e13ETH18-1\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c3\"\u003e \u003cp\u003eEthiopia\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c4\"\u003e \u003cp\u003e12-\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c5\"\u003e \u003cp\u003e3+\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003eTTRTF\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003e14GEO189-1\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c3\"\u003e \u003cp\u003eGeorgia\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c4\"\u003e \u003cp\u003e3-\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c5\"\u003e \u003cp\u003e2\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003c/tbody\u003e \u003c/colgroup\u003e \u003c/table\u003e\u003c/div\u003e \u003c/p\u003e \u003cp\u003eVariation between the parents for seedling resistance to TKTTF was observed. MN07098-6 had a mean linearized IT of 2.7. This finding was as expected given MN07098-6 carries \u003cem\u003eSr7a\u003c/em\u003e, a gene known to confer resistance to TKTTF (Edae et al. \u003cspan citationid=\"CR12\" class=\"CitationRef\"\u003e2018\u003c/span\u003e; Turner et al. \u003cspan citationid=\"CR89\" class=\"CitationRef\"\u003e2016\u003c/span\u003e). CItr 11390 was susceptible to TKTTF with a mean linearized IT of 8.8 (Fig. S4). The susceptibility of CItr 11390 was notable as it was found to carry the allele reported to be predictive of resistance at KASP_IWB72471, a marker linked to \u003cem\u003eSr11\u003c/em\u003e (Nirmala et al. \u003cspan citationid=\"CR54\" class=\"CitationRef\"\u003e2016\u003c/span\u003e). Analysis of the RILs found an approximately 1:1 segregation ratio of resistance with 57 resistant RILs and 47 susceptible RILs supporting the hypothesis that resistance to TKTTF in the population is provided by a single resistance gene (Table\u0026nbsp;\u003cspan refid=\"Tab3\" class=\"InternalRef\"\u003e3\u003c/span\u003e). A chi-square test resulted in the acceptance of the null hypothesis that resistance to TKTTF in the population is monogenic (χ\u003csup\u003e2\u003c/sup\u003e\u0026thinsp;=\u0026thinsp;1.0, \u003cem\u003ep\u003c/em\u003e\u0026thinsp;\u0026gt;\u0026thinsp;0.32; Table\u0026nbsp;\u003cspan refid=\"Tab3\" class=\"InternalRef\"\u003e3\u003c/span\u003e).\u003c/p\u003e \u003cp\u003e \u003cdiv class=\"gridtable\"\u003e\u003ctable float=\"Yes\" id=\"Tab3\" border=\"1\"\u003e \u003ccaption language=\"En\"\u003e \u003cdiv class=\"CaptionNumber\"\u003eTable 3\u003c/div\u003e \u003cdiv class=\"CaptionContent\"\u003e \u003cp\u003eChi-square test for segregation of stem rust reaction to TKTTF in the 15xR012 population\u003c/p\u003e \u003c/div\u003e \u003c/caption\u003e \u003ccolgroup cols=\"6\"\u003e \u003cdiv align=\"left\" class=\"colspec\" colname=\"c1\" colnum=\"1\"\u003e\u003c/div\u003e \u003cdiv align=\"left\" class=\"colspec\" colname=\"c2\" colnum=\"2\"\u003e\u003c/div\u003e \u003cdiv align=\"left\" class=\"colspec\" colname=\"c3\" colnum=\"3\"\u003e\u003c/div\u003e \u003cdiv align=\"left\" class=\"colspec\" colname=\"c4\" colnum=\"4\"\u003e\u003c/div\u003e \u003cdiv align=\"left\" class=\"colspec\" colname=\"c5\" colnum=\"5\"\u003e\u003c/div\u003e \u003cdiv align=\"left\" class=\"colspec\" colname=\"c6\" colnum=\"6\"\u003e\u003c/div\u003e \u003cthead\u003e \u003ctr\u003e \u003cth align=\"left\" colname=\"c1\"\u003e \u003cp\u003eRace\u003c/p\u003e \u003c/th\u003e \u003cth align=\"left\" colname=\"c2\"\u003e \u003cp\u003eNo. of resistant RILs (IT\u0026thinsp;\u0026le;\u0026thinsp;6)\u003c/p\u003e \u003c/th\u003e \u003cth align=\"left\" colname=\"c3\"\u003e \u003cp\u003eNo. of susceptible RILs (IT\u0026thinsp;\u0026gt;\u0026thinsp;6)\u003c/p\u003e \u003c/th\u003e \u003cth align=\"left\" colname=\"c4\"\u003e \u003cp\u003eObserved segregation\u003c/p\u003e \u003c/th\u003e \u003cth align=\"left\" colname=\"c5\"\u003e \u003cp\u003eχ\u003csup\u003e2\u003c/sup\u003e\u003c/p\u003e \u003c/th\u003e \u003cth align=\"left\" colname=\"c6\"\u003e \u003cp\u003eNumber of expected genes\u003c/p\u003e \u003c/th\u003e \u003c/tr\u003e \u003c/thead\u003e \u003ctbody\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003eTKTTF\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003e57\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c3\"\u003e \u003cp\u003e47\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c4\"\u003e \u003cp\u003e1:1\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c5\"\u003e \u003cp\u003e(1, N\u0026thinsp;=\u0026thinsp;104)\u0026thinsp;=\u0026thinsp;1.0, \u003cem\u003ep\u003c/em\u003e\u0026thinsp;\u0026gt;\u0026thinsp;0.33\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c6\"\u003e \u003cp\u003e1 gene\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003c/tbody\u003e \u003c/colgroup\u003e \u003c/table\u003e\u003c/div\u003e \u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec21\" class=\"Section2\"\u003e \u003ch2\u003eLinkage mapping for seedling resistance to TKTTF\u003c/h2\u003e \u003cp\u003eMQM analysis was performed to validate the effectiveness of \u003cem\u003eSr7a\u003c/em\u003e against TKTTF. A large-effect QTL was detected on 4AL (LOD\u0026thinsp;=\u0026thinsp;29.2; \u003cem\u003eR\u003c/em\u003e\u003csup\u003e\u003cem\u003e2\u003c/em\u003e\u003c/sup\u003e\u0026thinsp;=\u0026thinsp;69.8%; Table\u0026nbsp;\u003cspan refid=\"Tab4\" class=\"InternalRef\"\u003e4\u003c/span\u003e; Fig. S5) that colocated with the \u003cem\u003eSr7a\u003c/em\u003e-linked KASP markers screened on the population. Unexpectedly, a second resistance loci was retained in the final model on the long arm of chromosome 5A (LOD\u0026thinsp;=\u0026thinsp;29.2; \u003cem\u003eR\u003c/em\u003e\u003csup\u003e\u003cem\u003e2\u003c/em\u003e\u003c/sup\u003e\u0026thinsp;=\u0026thinsp;4.6%; Table\u0026nbsp;\u003cspan refid=\"Tab4\" class=\"InternalRef\"\u003e4\u003c/span\u003e; Fig. S5). Like \u003cem\u003eSr7a\u003c/em\u003e, the resistance allele at this locus is donated by the hexaploid parent, MN07098-6.\u003c/p\u003e \u003cp\u003e \u003cdiv class=\"gridtable\"\u003e\u003ctable float=\"Yes\" id=\"Tab4\" border=\"1\"\u003e \u003ccaption language=\"En\"\u003e \u003cdiv class=\"CaptionNumber\"\u003eTable 4\u003c/div\u003e \u003cdiv class=\"CaptionContent\"\u003e \u003cp\u003eQuantitative trait loci (QTL) for seedling resistance to TKTTF in the 15xR012 population detected by multiple QTL mapping (MQM). All QTL are significant at α\u0026thinsp;\u0026lt;\u0026thinsp;0.05\u003c/p\u003e \u003c/div\u003e \u003c/caption\u003e \u003ccolgroup cols=\"9\"\u003e \u003cdiv align=\"left\" class=\"colspec\" colname=\"c1\" colnum=\"1\"\u003e\u003c/div\u003e \u003cdiv align=\"left\" class=\"colspec\" colname=\"c2\" colnum=\"2\"\u003e\u003c/div\u003e \u003cdiv align=\"left\" class=\"colspec\" colname=\"c3\" colnum=\"3\"\u003e\u003c/div\u003e \u003cdiv align=\"left\" class=\"colspec\" colname=\"c4\" colnum=\"4\"\u003e\u003c/div\u003e \u003cdiv align=\"left\" class=\"colspec\" colname=\"c5\" colnum=\"5\"\u003e\u003c/div\u003e \u003cdiv align=\"char\" char=\".\" class=\"colspec\" colname=\"c6\" colnum=\"6\"\u003e\u003c/div\u003e \u003cdiv align=\"char\" char=\".\" class=\"colspec\" colname=\"c7\" colnum=\"7\"\u003e\u003c/div\u003e \u003cdiv align=\"char\" char=\".\" class=\"colspec\" colname=\"c8\" colnum=\"8\"\u003e\u003c/div\u003e \u003cdiv align=\"char\" char=\".\" class=\"colspec\" colname=\"c9\" colnum=\"9\"\u003e\u003c/div\u003e \u003cthead\u003e \u003ctr\u003e \u003cth align=\"left\" colname=\"c1\"\u003e \u003cp\u003eRace\u003c/p\u003e \u003c/th\u003e \u003cth align=\"left\" colname=\"c2\"\u003e \u003cp\u003eQTL\u003csup\u003ea\u003c/sup\u003e\u003c/p\u003e \u003c/th\u003e \u003cth align=\"left\" colname=\"c3\"\u003e \u003cp\u003eChr\u003c/p\u003e \u003c/th\u003e \u003cth align=\"left\" colspan=\"2\" nameend=\"c5\" namest=\"c4\"\u003e \u003cp\u003eFlanking markers\u003csup\u003eb\u003c/sup\u003e\u003c/p\u003e \u003c/th\u003e \u003cth align=\"left\" colname=\"c6\"\u003e \u003cp\u003ePos (cM)\u003c/p\u003e \u003c/th\u003e \u003cth align=\"left\" colname=\"c7\"\u003e \u003cp\u003eLOD\u003csup\u003ec\u003c/sup\u003e\u003c/p\u003e \u003c/th\u003e \u003cth align=\"left\" colname=\"c8\"\u003e \u003cp\u003e\u003cem\u003eR\u003c/em\u003e\u003csup\u003e\u003cem\u003e2 d\u003c/em\u003e\u003c/sup\u003e\u003c/p\u003e \u003c/th\u003e \u003cth align=\"left\" colname=\"c9\"\u003e \u003cp\u003eAdd\u003csup\u003ee\u003c/sup\u003e\u003c/p\u003e \u003c/th\u003e \u003c/tr\u003e \u003c/thead\u003e \u003ctbody\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e\u0026nbsp;\u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e\u0026nbsp;\u003c/td\u003e \u003ctd align=\"left\" colname=\"c3\"\u003e\u0026nbsp;\u003c/td\u003e \u003ctd align=\"left\" colname=\"c4\"\u003e \u003cp\u003eLeft\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c5\"\u003e \u003cp\u003eRight\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c6\"\u003e\u0026nbsp;\u003c/td\u003e \u003ctd align=\"left\" colname=\"c7\"\u003e\u0026nbsp;\u003c/td\u003e \u003ctd align=\"left\" colname=\"c8\"\u003e\u0026nbsp;\u003c/td\u003e \u003ctd align=\"left\" colname=\"c9\"\u003e\u0026nbsp;\u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003eTKTTF\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003e\u003cem\u003eQSr.umn-4A\u003c/em\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c3\"\u003e \u003cp\u003e4AL\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c4\"\u003e \u003cp\u003eS4A_720268179\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c5\"\u003e \u003cp\u003eS4A_728943337\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c6\"\u003e \u003cp\u003e7.0\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c7\"\u003e \u003cp\u003e29.2\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c8\"\u003e \u003cp\u003e69.8\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c9\"\u003e \u003cp\u003e2.7\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e\u0026nbsp;\u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003e\u003cem\u003eQSr.umn-5A.1\u003c/em\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c3\"\u003e \u003cp\u003e5AL\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c4\"\u003e \u003cp\u003eS5A_553400727\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c5\"\u003e \u003cp\u003eS5A_605408042\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c6\"\u003e \u003cp\u003e109.2\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c7\"\u003e \u003cp\u003e3.6\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c8\"\u003e \u003cp\u003e4.6\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c9\"\u003e \u003cp\u003e0.6\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003c/tbody\u003e \u003c/colgroup\u003e \u003ctfoot\u003e \u003ctr\u003e\u003ctd colspan=\"9\"\u003e\u003csup\u003ea\u003c/sup\u003e Named according to McIntosh et al. (\u003cspan citationid=\"CR49\" class=\"CitationRef\"\u003e2020\u003c/span\u003e)\u003c/td\u003e\u003c/tr\u003e \u003ctr\u003e\u003ctd colspan=\"9\"\u003e\u003csup\u003eb\u003c/sup\u003e 95% Bayes credible interval\u003c/td\u003e\u003c/tr\u003e \u003ctr\u003e\u003ctd colspan=\"9\"\u003e\u003csup\u003ec\u003c/sup\u003e Peak logarithm of odds\u003c/td\u003e\u003c/tr\u003e \u003ctr\u003e\u003ctd colspan=\"9\"\u003e\u003csup\u003ed\u003c/sup\u003e Phenotypic variance explained\u003c/td\u003e\u003c/tr\u003e \u003ctr\u003e\u003ctd colspan=\"9\"\u003e\u003csup\u003ee\u003c/sup\u003e Estimated additive effect of QTL; negative value indicates that the allele donor is CItr 11390\u003c/td\u003e\u003c/tr\u003e \u003c/tfoot\u003e \u003c/table\u003e\u003c/div\u003e \u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec22\" class=\"Section2\"\u003e \u003ch2\u003eField-based disease evaluation\u003c/h2\u003e \u003cp\u003eCItr 11390 displayed high levels of resistance in Kenya in both 2019 and 2020, but the accession was far less resistant in both Ethiopian nurseries (Table\u0026nbsp;\u003cspan refid=\"Tab5\" class=\"InternalRef\"\u003e5\u003c/span\u003e). MN07098-6 was susceptible in all four environments with mean severity ranging from 60 to 75% within individual environments (Table \u003cspan refid=\"MOESM2\" class=\"InternalRef\"\u003eS2\u003c/span\u003e). Disease pressure was especially high in Ethiopia in 2019, and scoring was conducted later than the optimal window. Both parents and the RIL population displayed completely susceptible infection responses and almost no variation in severity was observed (Tables S2 and S3). Raw stem rust CI distributions were continuous and positively-skewed in Kenya in 2019 and 2020 (Fig.\u0026nbsp;\u003cspan refid=\"Fig2\" class=\"InternalRef\"\u003e2\u003c/span\u003e). Visual assessment of a Q-Q plot indicated the Eth20 data had a slightly negative skew. Residuals for Ken19 and Ken20 were in violation of the assumption of normality based on Shapiro-Wilk tests (\u003cem\u003ep\u003c/em\u003e\u0026thinsp;\u0026lt;\u0026thinsp;0.05). Only the residuals for Ken20 were heteroscedastic based on Levene\u0026rsquo;s test (\u003cem\u003ep\u003c/em\u003e\u0026thinsp;\u0026lt;\u0026thinsp;0.05).\u003c/p\u003e \u003cp\u003e \u003c/p\u003e \u003cp\u003e \u003cdiv class=\"gridtable\"\u003e\u003ctable float=\"Yes\" id=\"Tab5\" border=\"1\"\u003e \u003ccaption language=\"En\"\u003e \u003cdiv class=\"CaptionNumber\"\u003eTable 5\u003c/div\u003e \u003cdiv class=\"CaptionContent\"\u003e \u003cp\u003eMean stem rust coefficient of infection values, standard deviations, and ranges in the 15xR012 mapping population in four field environments.\u003c/p\u003e \u003c/div\u003e \u003c/caption\u003e \u003ccolgroup cols=\"5\"\u003e \u003cdiv align=\"left\" class=\"colspec\" colname=\"c1\" colnum=\"1\"\u003e\u003c/div\u003e \u003cdiv align=\"left\" class=\"colspec\" colname=\"c2\" colnum=\"2\"\u003e\u003c/div\u003e \u003cdiv align=\"left\" class=\"colspec\" colname=\"c3\" colnum=\"3\"\u003e\u003c/div\u003e \u003cdiv align=\"left\" class=\"colspec\" colname=\"c4\" colnum=\"4\"\u003e\u003c/div\u003e \u003cdiv align=\"left\" class=\"colspec\" colname=\"c5\" colnum=\"5\"\u003e\u003c/div\u003e \u003cthead\u003e \u003ctr\u003e \u003cth align=\"left\" colname=\"c1\"\u003e \u003cp\u003eEnvironments\u003csup\u003ea\u003c/sup\u003e\u003c/p\u003e \u003c/th\u003e \u003cth align=\"left\" colspan=\"2\" nameend=\"c3\" namest=\"c2\"\u003e \u003cp\u003eParent mean coefficient of infection\u003c/p\u003e \u003c/th\u003e \u003cth align=\"left\" colspan=\"2\" nameend=\"c5\" namest=\"c4\"\u003e \u003cp\u003ePopulation\u003c/p\u003e \u003c/th\u003e \u003c/tr\u003e \u003c/thead\u003e \u003ctbody\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e\u0026nbsp;\u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003eCItr 11390\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c3\"\u003e \u003cp\u003eMN07098-6\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c4\"\u003e \u003cp\u003eMean\u0026thinsp;\u0026plusmn;\u0026thinsp;SD\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c5\"\u003e \u003cp\u003eRange\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003eEth19\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003e60.0\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c3\"\u003e \u003cp\u003e60.0\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c4\"\u003e \u003cp\u003e58.3\u0026thinsp;\u0026plusmn;\u0026thinsp;3.1\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c5\"\u003e \u003cp\u003e45.0\u0026ndash;60.0\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003eEth20\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003e40.0\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c3\"\u003e \u003cp\u003e65.0\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c4\"\u003e \u003cp\u003e59.2\u0026thinsp;\u0026plusmn;\u0026thinsp;12.3\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c5\"\u003e \u003cp\u003e24.0\u0026ndash;80.0\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003eKen19\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003e2.5\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c3\"\u003e \u003cp\u003e60.0\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c4\"\u003e \u003cp\u003e35.5\u0026thinsp;\u0026plusmn;\u0026thinsp;22.4\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c5\"\u003e \u003cp\u003e2.0\u0026ndash;90.0\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003eKen20\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003e8.3\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c3\"\u003e \u003cp\u003e75.0\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c4\"\u003e \u003cp\u003e41.0\u0026thinsp;\u0026plusmn;\u0026thinsp;21.8\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c5\"\u003e \u003cp\u003e8.0-100\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003c/tbody\u003e \u003c/colgroup\u003e \u003ctfoot\u003e \u003ctr\u003e\u003ctd colspan=\"5\"\u003e\u003csup\u003ea\u003c/sup\u003e Eth19\u0026thinsp;=\u0026thinsp;Debre Zeit, Ethiopia 2019; Eth20\u0026thinsp;=\u0026thinsp;Debre Zeit, Ethiopia 2020; Ken19\u0026thinsp;=\u0026thinsp;Njoro, Kenya 2019; Ken20\u0026thinsp;=\u0026thinsp;Njoro, Kenya 2020\u003c/td\u003e\u003c/tr\u003e \u003c/tfoot\u003e \u003c/table\u003e\u003c/div\u003e \u003c/p\u003e \u003cp\u003eSquare root transformation of the Ken19 values resulted in the restoration of normality. Square root transformation of the Ken20 values also resulted in the restoration of normality and homoscedasticity, as well as marked improvement in the Q-Q plot. Therefore, the transformed data were used for Ken19 and Ken20.The Eth20 data were divided by 100 followed by an arcsine transformation. While a Shapiro-Wilk test indicated that the transformed data were normally distributed, little to no differences were viewed between the Q-Q plots of the untransformed and transformed data fitted with the same model. Given these observation, as well as the homogeneity of the untransformed residuals and the increased difficulty of interpreting QTL mapping results from transformed datasets, the untransformed data from Eth20 were used for further analyses.\u003c/p\u003e \u003cp\u003eThe 15xR012 population was found to segregate for \u003cem\u003eSr57/Lr34/Yr18\u003c/em\u003e due to near-isogenic MN07098-6 lines being used for population development. \u003cem\u003eSr57\u003c/em\u003e/\u003cem\u003eLr34\u003c/em\u003e/\u003cem\u003eYr18\u003c/em\u003e was found to have a significant effect on CI in a mixed model using the data from all environments. Due to a complete absence of additional markers linked to \u003cem\u003eSr57/Lr34/Yr18\u003c/em\u003e on chromosome 7D, the effect of the gene was adjusted using mixed models fitted within each environment. The estimated marginal means of CI for the resistance allele at \u003cem\u003eSr57\u003c/em\u003e/\u003cem\u003eLr34\u003c/em\u003e/\u003cem\u003eYr18\u003c/em\u003e was 30.0%, 46.5%, and 11.1% lower than the susceptible allele in Ken19, Ken20, and Eth20, respectively. After adjustment, mean CI increased slightly in Ken19, Ken20, and Eth20 (Fig.\u0026nbsp;\u003cspan refid=\"Fig3\" class=\"InternalRef\"\u003e3\u003c/span\u003e). Adjusted CI values of the RILs between all environments were significantly correlated (Table\u0026nbsp;\u003cspan refid=\"Tab6\" class=\"InternalRef\"\u003e6\u003c/span\u003e).\u003c/p\u003e \u003cp\u003e \u003c/p\u003e \u003cp\u003e \u003cdiv class=\"gridtable\"\u003e\u003ctable float=\"Yes\" id=\"Tab6\" border=\"1\"\u003e \u003ccaption language=\"En\"\u003e \u003cdiv class=\"CaptionNumber\"\u003eTable 6\u003c/div\u003e \u003cdiv class=\"CaptionContent\"\u003e \u003cp\u003ePearson correlation coefficients (\u003cem\u003er\u003c/em\u003e) of adjusted stem rust coefficient of infection values observed in the 15xR012 population in four field environments in Ethiopia and Kenya in 2019 and 2020\u003c/p\u003e \u003c/div\u003e \u003c/caption\u003e \u003ccolgroup cols=\"5\"\u003e \u003cdiv align=\"left\" class=\"colspec\" colname=\"c1\" colnum=\"1\"\u003e\u003c/div\u003e \u003cdiv align=\"char\" char=\".\" class=\"colspec\" colname=\"c2\" colnum=\"2\"\u003e\u003c/div\u003e \u003cdiv align=\"char\" char=\".\" class=\"colspec\" colname=\"c3\" colnum=\"3\"\u003e\u003c/div\u003e \u003cdiv align=\"char\" char=\".\" class=\"colspec\" colname=\"c4\" colnum=\"4\"\u003e\u003c/div\u003e \u003cdiv align=\"left\" class=\"colspec\" colname=\"c5\" colnum=\"5\"\u003e\u003c/div\u003e \u003cthead\u003e \u003ctr\u003e \u003cth align=\"left\" colname=\"c1\"\u003e\u0026nbsp;\u003c/th\u003e \u003cth align=\"left\" colname=\"c2\"\u003e \u003cp\u003eEth19\u003c/p\u003e \u003c/th\u003e \u003cth align=\"left\" colname=\"c3\"\u003e \u003cp\u003eEth20\u003c/p\u003e \u003c/th\u003e \u003cth align=\"left\" colname=\"c4\"\u003e \u003cp\u003eKen19\u003c/p\u003e \u003c/th\u003e \u003cth align=\"left\" colname=\"c5\"\u003e \u003cp\u003eKen20\u003c/p\u003e \u003c/th\u003e \u003c/tr\u003e \u003c/thead\u003e \u003ctbody\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003eEth20\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c2\"\u003e \u003cp\u003e0.37\u003csup\u003e***\u003c/sup\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c3\"\u003e\u0026nbsp;\u003c/td\u003e \u003ctd align=\"left\" colname=\"c4\"\u003e\u0026nbsp;\u003c/td\u003e \u003ctd align=\"left\" colname=\"c5\"\u003e\u0026nbsp;\u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003eKen19\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c2\"\u003e \u003cp\u003e0.23\u003csup\u003e*\u003c/sup\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c3\"\u003e \u003cp\u003e0.33\u003csup\u003e***\u003c/sup\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c4\"\u003e\u0026nbsp;\u003c/td\u003e \u003ctd align=\"left\" colname=\"c5\"\u003e\u0026nbsp;\u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003eKen20\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c2\"\u003e \u003cp\u003e0.20\u003csup\u003e*\u003c/sup\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c3\"\u003e \u003cp\u003e0.39\u003csup\u003e***\u003c/sup\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c4\"\u003e \u003cp\u003e0.66\u003csup\u003e***\u003c/sup\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c5\"\u003e\u0026nbsp;\u003c/td\u003e \u003c/tr\u003e \u003c/tbody\u003e \u003c/colgroup\u003e \u003ctfoot\u003e \u003ctr\u003e\u003ctd colspan=\"5\"\u003eEth19\u0026thinsp;=\u0026thinsp;Debre Zeit, Ethiopia 2019; Eth20\u0026thinsp;=\u0026thinsp;Debre Zeit, Ethiopia 2020; Ken19\u0026thinsp;=\u0026thinsp;Njoro, Kenya 2019; Ken20\u0026thinsp;=\u0026thinsp;Njoro, Kenya 2020\u003c/td\u003e\u003c/tr\u003e \u003ctr\u003e\u003ctd colspan=\"5\"\u003e\u003csup\u003e*\u003c/sup\u003e Correlation is significant at α\u0026thinsp;\u0026lt;\u0026thinsp;0.05\u003c/td\u003e\u003c/tr\u003e \u003ctr\u003e\u003ctd colspan=\"5\"\u003e\u003csup\u003e**\u003c/sup\u003e Correlation is significant at α\u0026thinsp;\u0026lt;\u0026thinsp;0.01\u003c/td\u003e\u003c/tr\u003e \u003ctr\u003e\u003ctd colspan=\"5\"\u003e\u003csup\u003e***\u003c/sup\u003e Correlation is significant at α\u0026thinsp;\u0026lt;\u0026thinsp;0.001\u003c/td\u003e\u003c/tr\u003e \u003c/tfoot\u003e \u003c/table\u003e\u003c/div\u003e \u003c/p\u003e \u003cdiv id=\"Sec23\" class=\"Section3\"\u003e \u003ch2\u003eLinkage mapping of stem rust resistance\u003c/h2\u003e \u003cp\u003eMultiple QTL mapping (MQM) detected significant QTL associated with APR in three of the four African nurseries (Table\u0026nbsp;\u003cspan refid=\"Tab7\" class=\"InternalRef\"\u003e7\u003c/span\u003e; Fig. S6). No QTL were detected in Eth19. In total, seven QTL were detected in at least one environment. The LOD scores of the detected QTL ranged from 3.3 to 20.9. Of the seven QTL, two were detected in more than one environment. \u003cem\u003eQSr.umn-6BL\u003c/em\u003e was detected in Eth20 and Ken20. It explained 10.2% of the observed phenotypic variation in Eth20 and 10.5% in Ken20. While \u003cem\u003eQSr.umn-6BL\u003c/em\u003e was not retained in the multiple QTL model, it was detected by single-QTL and two dimensional genome scans (Table S4).\u003c/p\u003e \u003cp\u003eThe second multi-environment QTL, \u003cem\u003eQSr.umn-2A\u003c/em\u003e, was detected in Kenya in 2019 and 2020. \u003cem\u003eQSr.umn-2A\u003c/em\u003e explained 42.8% of the phenotypic variation in 2019 and 44.6% of variation in 2020. The peaks for \u003cem\u003eQSr.umn-2A\u003c/em\u003e were sharply centered around KASP marker Sr63_IWB32429 for Ken19 and Ken20, indicating that \u003cem\u003eQSr.umn-2A\u003c/em\u003e could be \u003cem\u003eSr63\u003c/em\u003e.\u003c/p\u003e \u003cp\u003e \u003cem\u003eQSr.umn-2A\u003c/em\u003e and \u003cem\u003eQSr.umn-6BL\u003c/em\u003e were each contributed by the tetraploid parent, CItr 11390. The five other QTL were only detected in a single environment and explained 5.0-21.5% of the observed phenotypic variation within their respective environments. Notably, the interval of \u003cem\u003eQSr.umn-5A.2\u003c/em\u003e co-located with the TKTTF ASR-associated QTL, \u003cem\u003eQSr.umn-5A.1\u003c/em\u003e. The resistance allele was contributed by MN07098-6 at both \u003cem\u003eQSr.umn-5A.1\u003c/em\u003e and \u003cem\u003eQSr.umn-5A.2\u003c/em\u003e. Of the remaining single-environment QTL, three were contributed by MN07098-6 and one was contributed by CItr 11390 (Table\u0026nbsp;\u003cspan refid=\"Tab7\" class=\"InternalRef\"\u003e7\u003c/span\u003e).\u003c/p\u003e \u003cp\u003e \u003cdiv class=\"gridtable\"\u003e\u003ctable float=\"Yes\" id=\"Tab7\" border=\"1\"\u003e \u003ccaption language=\"En\"\u003e \u003cdiv class=\"CaptionNumber\"\u003eTable 7\u003c/div\u003e \u003cdiv class=\"CaptionContent\"\u003e \u003cp\u003eQuantitative trait loci (QTL) for field-based stem rust resistance in the 15xR012 population detected by multiple QTL mapping (MQM). All QTL are significant at α\u0026thinsp;\u0026lt;\u0026thinsp;0.05\u003c/p\u003e \u003c/div\u003e \u003c/caption\u003e \u003ccolgroup cols=\"9\"\u003e \u003cdiv align=\"left\" class=\"colspec\" colname=\"c1\" colnum=\"1\"\u003e\u003c/div\u003e \u003cdiv align=\"left\" class=\"colspec\" colname=\"c2\" colnum=\"2\"\u003e\u003c/div\u003e \u003cdiv align=\"left\" class=\"colspec\" colname=\"c3\" colnum=\"3\"\u003e\u003c/div\u003e \u003cdiv align=\"left\" class=\"colspec\" colname=\"c4\" colnum=\"4\"\u003e\u003c/div\u003e \u003cdiv align=\"left\" class=\"colspec\" colname=\"c5\" colnum=\"5\"\u003e\u003c/div\u003e \u003cdiv align=\"char\" char=\".\" class=\"colspec\" colname=\"c6\" colnum=\"6\"\u003e\u003c/div\u003e \u003cdiv align=\"char\" char=\".\" class=\"colspec\" colname=\"c7\" colnum=\"7\"\u003e\u003c/div\u003e \u003cdiv align=\"char\" char=\".\" class=\"colspec\" colname=\"c8\" colnum=\"8\"\u003e\u003c/div\u003e \u003cdiv align=\"left\" class=\"colspec\" colname=\"c9\" colnum=\"9\"\u003e\u003c/div\u003e \u003cthead\u003e \u003ctr\u003e \u003cth align=\"left\" colname=\"c1\"\u003e \u003cp\u003eEnvironment\u003csup\u003ea\u003c/sup\u003e\u003c/p\u003e \u003c/th\u003e \u003cth align=\"left\" colname=\"c2\"\u003e \u003cp\u003eQTL\u003csup\u003eb\u003c/sup\u003e\u003c/p\u003e \u003c/th\u003e \u003cth align=\"left\" colname=\"c3\"\u003e \u003cp\u003eChr\u003c/p\u003e \u003c/th\u003e \u003cth align=\"left\" colspan=\"2\" nameend=\"c5\" namest=\"c4\"\u003e \u003cp\u003eFlanking markers\u003csup\u003ec\u003c/sup\u003e\u003c/p\u003e \u003c/th\u003e \u003cth align=\"left\" colname=\"c6\"\u003e \u003cp\u003ePos (cM)\u003c/p\u003e \u003c/th\u003e \u003cth align=\"left\" colname=\"c7\"\u003e \u003cp\u003eLOD\u003csup\u003ed\u003c/sup\u003e\u003c/p\u003e \u003c/th\u003e \u003cth align=\"left\" colname=\"c8\"\u003e \u003cp\u003e\u003cem\u003eR\u003c/em\u003e\u003csup\u003e\u003cem\u003e2 e\u003c/em\u003e\u003c/sup\u003e\u003c/p\u003e \u003c/th\u003e \u003cth align=\"left\" colname=\"c9\"\u003e \u003cp\u003eAdd\u003csup\u003ef\u003c/sup\u003e\u003c/p\u003e \u003c/th\u003e \u003c/tr\u003e \u003c/thead\u003e \u003ctbody\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e\u0026nbsp;\u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e\u0026nbsp;\u003c/td\u003e \u003ctd align=\"left\" colname=\"c3\"\u003e\u0026nbsp;\u003c/td\u003e \u003ctd align=\"left\" colname=\"c4\"\u003e \u003cp\u003eLeft\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c5\"\u003e \u003cp\u003eRight\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c6\"\u003e\u0026nbsp;\u003c/td\u003e \u003ctd align=\"left\" colname=\"c7\"\u003e\u0026nbsp;\u003c/td\u003e \u003ctd align=\"left\" colname=\"c8\"\u003e\u0026nbsp;\u003c/td\u003e \u003ctd align=\"left\" colname=\"c9\"\u003e\u0026nbsp;\u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003eEth20\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003e\u003cem\u003eQSr.umn-1A\u003c/em\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c3\"\u003e \u003cp\u003e1AS\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c4\"\u003e \u003cp\u003eS1A_5914856\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c5\"\u003e \u003cp\u003eS1A_15479164\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c6\"\u003e \u003cp\u003e10.0\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c7\"\u003e \u003cp\u003e6.4\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c8\"\u003e \u003cp\u003e21.5\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c9\"\u003e \u003cp\u003e1.6\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e\u0026nbsp;\u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003e\u003cem\u003eQSr.umn-6BL\u003c/em\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c3\"\u003e \u003cp\u003e6BL\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c4\"\u003e \u003cp\u003eS6B_467576061\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c5\"\u003e \u003cp\u003eS6B_651836748\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c6\"\u003e \u003cp\u003e13.0\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c7\"\u003e \u003cp\u003e3.3\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c8\"\u003e \u003cp\u003e10.2\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c9\"\u003e \u003cp\u003e-1.1\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003eKen19\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003e\u003cem\u003eQSr.umn-2A\u003c/em\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c3\"\u003e \u003cp\u003e2AL\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c4\"\u003e \u003cp\u003eS2A_681051011\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c5\"\u003e \u003cp\u003eS2A_702504398\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c6\"\u003e \u003cp\u003e139.0\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c7\"\u003e \u003cp\u003e14.1\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c8\"\u003e \u003cp\u003e42.8\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c9\"\u003e \u003cp\u003e-0.8\u003csup\u003eg\u003c/sup\u003e\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e\u0026nbsp;\u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003e\u003cem\u003eQSr.umn-5A.2\u003c/em\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c3\"\u003e \u003cp\u003e5AL\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c4\"\u003e \u003cp\u003eS5A_488255473\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c5\"\u003e \u003cp\u003eS5A_659345920\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c6\"\u003e \u003cp\u003e77.0\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c7\"\u003e \u003cp\u003e4.2\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c8\"\u003e \u003cp\u003e10.1\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c9\"\u003e \u003cp\u003e0.4\u003csup\u003eg\u003c/sup\u003e\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003eKen20\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003e\u003cem\u003eQSr.umn-2A\u003c/em\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c3\"\u003e \u003cp\u003e2AL\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c4\"\u003e \u003cp\u003eS2A_681051011\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c5\"\u003e \u003cp\u003eSr63_IWB32429\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c6\"\u003e \u003cp\u003e136.0\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c7\"\u003e \u003cp\u003e20.9\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c8\"\u003e \u003cp\u003e44.6\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c9\"\u003e \u003cp\u003e-0.67\u003csup\u003eg\u003c/sup\u003e\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e\u0026nbsp;\u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003e\u003cem\u003eQSr.umn-3A\u003c/em\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c3\"\u003e \u003cp\u003e3AL\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c4\"\u003e \u003cp\u003eS3A_603255953\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c5\"\u003e \u003cp\u003eS3A_654993190\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c6\"\u003e \u003cp\u003e105.0\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c7\"\u003e \u003cp\u003e3.6\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c8\"\u003e \u003cp\u003e5.0\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c9\"\u003e \u003cp\u003e0.22\u003csup\u003eg\u003c/sup\u003e\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e\u0026nbsp;\u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003e\u003cem\u003eQSr.umn-3B\u003c/em\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c3\"\u003e \u003cp\u003e3BL\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c4\"\u003e \u003cp\u003eS3B_429246336\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c5\"\u003e \u003cp\u003eS3B_490054756\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c6\"\u003e \u003cp\u003e33.1\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c7\"\u003e \u003cp\u003e6.6\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c8\"\u003e \u003cp\u003e10.0\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c9\"\u003e \u003cp\u003e0.32\u003csup\u003eg\u003c/sup\u003e\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e\u0026nbsp;\u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003e\u003cem\u003eQSr.umn-6BS\u003c/em\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c3\"\u003e \u003cp\u003e6BS\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c4\"\u003e \u003cp\u003eS6B_17383541\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c5\"\u003e \u003cp\u003eS6B_36319780\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c6\"\u003e \u003cp\u003e16.7\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c7\"\u003e \u003cp\u003e10.0\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c8\"\u003e \u003cp\u003e16.2\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c9\"\u003e \u003cp\u003e-0.22\u003csup\u003eg\u003c/sup\u003e\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e\u0026nbsp;\u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003e\u003cem\u003eQSr.umn-6BL\u003c/em\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c3\"\u003e \u003cp\u003e6BL\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c4\"\u003e \u003cp\u003eS6B_572334604\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c5\"\u003e \u003cp\u003eS6B_602404199\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c6\"\u003e \u003cp\u003e31.7\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c7\"\u003e \u003cp\u003e6.9\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c8\"\u003e \u003cp\u003e10.5\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c9\"\u003e \u003cp\u003e-0.31\u003csup\u003eg\u003c/sup\u003e\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003c/tbody\u003e \u003c/colgroup\u003e \u003ctfoot\u003e \u003ctr\u003e\u003ctd colspan=\"9\"\u003e\u003csup\u003ea\u003c/sup\u003e Eth20\u0026thinsp;=\u0026thinsp;Debre Zeit, Ethiopia 2020; Ken19\u0026thinsp;=\u0026thinsp;Njoro, Kenya 2019; Ken20\u0026thinsp;=\u0026thinsp;Njoro, Kenya 2020\u003c/td\u003e\u003c/tr\u003e \u003ctr\u003e\u003ctd colspan=\"9\"\u003e\u003csup\u003eb\u003c/sup\u003e Named according to McIntosh et al. (\u003cspan citationid=\"CR49\" class=\"CitationRef\"\u003e2020\u003c/span\u003e)\u003c/td\u003e\u003c/tr\u003e \u003ctr\u003e\u003ctd colspan=\"9\"\u003e\u003csup\u003ec\u003c/sup\u003e 95% Bayes credible interval\u003c/td\u003e\u003c/tr\u003e \u003ctr\u003e\u003ctd colspan=\"9\"\u003e\u003csup\u003ed\u003c/sup\u003e Peak logarithm of odds\u003c/td\u003e\u003c/tr\u003e \u003ctr\u003e\u003ctd colspan=\"9\"\u003e\u003csup\u003ee\u003c/sup\u003e Phenotypic variance explained\u003c/td\u003e\u003c/tr\u003e \u003ctr\u003e\u003ctd colspan=\"9\"\u003e\u003csup\u003ef\u003c/sup\u003e Estimated additive effect of QTL; negative value indicates that the allele donor is CItr 11390\u003c/td\u003e\u003c/tr\u003e \u003ctr\u003e\u003ctd colspan=\"9\"\u003e\u003csup\u003eg\u003c/sup\u003e Value is reported as the square root of coefficient of infection\u003c/td\u003e\u003c/tr\u003e \u003c/tfoot\u003e \u003c/table\u003e\u003c/div\u003e \u003c/p\u003e \u003c/div\u003e \u003c/div\u003e \u003cdiv id=\"Sec24\" class=\"Section2\"\u003e \u003ch2\u003eAssessment of QTL interactions and combinations\u003c/h2\u003e \u003cp\u003eTests for epistasis while fitting the multiple QTL models detected one significant interaction within Ken20. The interaction between \u003cem\u003eQSr.umn-2A\u003c/em\u003e and \u003cem\u003eQSr.umn-6BS\u003c/em\u003e explained 10.9% of the observed phenotypic variation in Ken20 (\u003cem\u003ep\u003c/em\u003e\u0026thinsp;\u0026lt;\u0026thinsp;0.001). The resistance allele was donated by CItr 11390 at both QTL. The mean CI of lines carrying only the resistance allele at \u003cem\u003eQSr.umn-2A\u003c/em\u003e was significantly lower than those carrying only the resistance allele at \u003cem\u003eQSr.umn-6BS\u003c/em\u003e. However, the difference in mean CI between lines only carrying the resistance allele at \u003cem\u003eQSr.umn-2A\u003c/em\u003e and those carrying resistance alleles at both loci was not significant (Fig.\u0026nbsp;\u003cspan refid=\"Fig4\" class=\"InternalRef\"\u003e4\u003c/span\u003ea, \u003cspan refid=\"Fig4\" class=\"InternalRef\"\u003e4\u003c/span\u003eb, Table S5).\u003c/p\u003e \u003cp\u003e \u003c/p\u003e \u003cp\u003e \u003cem\u003eQSr.umn-6BL\u003c/em\u003e was detected in two environments by MQM while \u003cem\u003eQSr.umn-2A\u003c/em\u003e was detected in Kenya in both years. In addition to the detected QTL, the effect of \u003cem\u003eSr57/Lr34\u003c/em\u003e/\u003cem\u003eYr18\u003c/em\u003e was significant on CI in all three environments. Multiple RILs in each of the eight possible combinations of these three loci were found among the population: No QTL, \u003cem\u003eSr57\u003c/em\u003e, \u003cem\u003eQSr.umn-2A\u003c/em\u003e (hereafter referred to as \u003cem\u003eSr63\u003c/em\u003e), \u003cem\u003eQSr.umn-6BL\u003c/em\u003e, \u003cem\u003eSr57\u003c/em\u003e\u0026thinsp;+\u0026thinsp;\u003cem\u003eSr63\u003c/em\u003e, \u003cem\u003eSr57\u003c/em\u003e\u0026thinsp;+\u0026thinsp;\u003cem\u003eQSr.umn-6BL\u003c/em\u003e, \u003cem\u003eSr63\u003c/em\u003e\u0026thinsp;+\u0026thinsp;\u003cem\u003eQSr.umn-6BL\u003c/em\u003e, and \u003cem\u003eSr57\u003c/em\u003e\u0026thinsp;+\u0026thinsp;\u003cem\u003eSr63\u003c/em\u003e\u0026thinsp;+\u0026thinsp;\u003cem\u003eQSr.umn-6B\u003c/em\u003e.\u003c/p\u003e \u003cp\u003eSignificant differences were observed among the allelic combinations in all environments. In Eth20, the mean CI of RILs carrying a single gene was lower than lines lacking all three genes. However, these differences were not significant. All combinations of two or all three resistance alleles resulted in a significant reduction in mean CI (Fig.\u0026nbsp;\u003cspan refid=\"Fig5\" class=\"InternalRef\"\u003e5\u003c/span\u003e, Table S6). The mean unadjusted CI of lines carrying all three genes was 52.0 compared to 69.3 for lines lacking all three genes.\u003c/p\u003e \u003cp\u003eWithin the Kenyan environments, many of the pairwise comparisons were not statistically significant. However, this is possibly due to the low sample size of RILs constituting each of the eight groups of gene combinations (Fig.\u0026nbsp;\u003cspan refid=\"Fig5\" class=\"InternalRef\"\u003e5\u003c/span\u003e, Table S6). Despite this finding, additive relationships were observed among the three genes. The mean unadjusted CI values of RILs carrying no resistance alleles at the three loci (Ken19: 61.6; Ken20: 66.9) were over three times greater than those carrying the resistance allele each of the three genes in both environments (Ken19: 18.2; Ken20: 18.5).\u003c/p\u003e \u003cp\u003e \u003c/p\u003e \u003cdiv id=\"Sec25\" class=\"Section3\"\u003e \u003ch2\u003eMolecular marker development\u003c/h2\u003e \u003cp\u003ePACE \u0026trade; SNP genotyping assays were developed for \u003cem\u003eQSr.umn-6BL\u003c/em\u003e at loci approximately flanking the QTL peaks detected by single-QTL genome scans in Eth20, Ken19, and Ken20. The markers with the clearest results were determined to be at loci correspond to positions 655,795,233 and 675,524,532 on chromosome 6B of RefSeq v2.1 assembly (Table S7). Genotyping calls for both assays were found to be concordant with the GBS markers (Table S8).\u003c/p\u003e \u003c/div\u003e \u003c/div\u003e"},{"header":"Discussion","content":"\u003cp\u003eThe discovery and deployment of new stem rust resistance genes from diverse sources is vital to the sustained production of wheat on both regional and global scales. Since the initial discovery of \u003cem\u003ePgt\u003c/em\u003e race TTKSK (Ug99) in Uganda in 1998 (Pretorius et al. \u003cspan citationid=\"CR67\" class=\"CitationRef\"\u003e2000\u003c/span\u003e), an additional fourteen variants have been identified with unique virulence profiles. The Ug99 race group has moved well beyond Uganda, with variants detected in fourteen countries spanning a geographic range from Iraq to South Africa (Terefe et al. 2018; Pretorius et al. \u003cspan citationid=\"CR66\" class=\"CitationRef\"\u003e2012\u003c/span\u003e). The Ug99 race group is not the only threat posed by \u003cem\u003ePgt\u003c/em\u003e. The need for continued monitoring and resistance breeding for other \u003cem\u003ePgt\u003c/em\u003e races is evidenced by the detection and geographic expansion of races TKTTF, TTRTF, JRCQC, and TTTTF, some of which are virulent to widely deployed resistance genes including \u003cem\u003eSrTmp\u003c/em\u003e and \u003cem\u003eSr13b\u003c/em\u003e (Letta et al. \u003cspan citationid=\"CR40\" class=\"CitationRef\"\u003e2013\u003c/span\u003e; Olivera et al. \u003cspan citationid=\"CR56\" class=\"CitationRef\"\u003e2015\u003c/span\u003e; Patpour et al. \u003cspan citationid=\"CR62\" class=\"CitationRef\"\u003e2020\u003c/span\u003e).\u003c/p\u003e \u003cp\u003eThe challenge of discovering novel sources of stem rust resistance for wheat is aided by the rich gene pools available to wheat breeders and pathologists (Rowland and Kerber \u003cspan citationid=\"CR75\" class=\"CitationRef\"\u003e1974\u003c/span\u003e). However, as demonstrated by this study, the integration of ancestral wheat species and other members of the \u003cem\u003eAegilops\u003c/em\u003e tribe into downstream breeding is not a straightforward process (Leigh et al. \u003cspan citationid=\"CR39\" class=\"CitationRef\"\u003e2022\u003c/span\u003e). Of the 121 BC\u003csub\u003e1\u003c/sub\u003eF5-derived RILs developed, 30 were aneuploids and removed from further phenotypic analysis and QTL mapping (Table \u003cspan refid=\"MOESM1\" class=\"InternalRef\"\u003eS1\u003c/span\u003e).\u003c/p\u003e \u003cp\u003eStem rust infection was much greater in both Ethiopian environments than either Kenyan environment. It is possible this difference was driven by the warmer environment in Ethiopia, as temperature has been found to significantly impact the effectiveness of ASR and APR genes (Chen et al. \u003cspan citationid=\"CR11\" class=\"CitationRef\"\u003e2018\u003c/span\u003e; Gao et al. \u003cspan citationid=\"CR15\" class=\"CitationRef\"\u003e2019\u003c/span\u003e; McIntosh et al. \u003cspan citationid=\"CR50\" class=\"CitationRef\"\u003e1995\u003c/span\u003e; Zhang et al. \u003cspan citationid=\"CR94\" class=\"CitationRef\"\u003e2017\u003c/span\u003e). Alternatively, the differences could be driven by the presence of additional races in Ethiopia. Races TRTTF and JRCQC have previously been identified from isolates collected in the Debre Zeit nursery (Olivera et al. \u003cspan citationid=\"CR57\" class=\"CitationRef\"\u003e2012\u003c/span\u003e). While the 15xR012 population was not screened for resistance to JRCQC, all RILs were susceptible to TRTTF at the seedling stage (Fig. \u003cspan refid=\"MOESM3\" class=\"InternalRef\"\u003eS3\u003c/span\u003e).\u003c/p\u003e \u003cp\u003eIn addition to TRTTF, the population was susceptible to Ug99 races TTKSK and TTKTT. Of the \u003cem\u003ePgt\u003c/em\u003e isolates screened, only seedling resistance to TKTTF was segregating within the population. The results from our seedling stage screenings contradict previous reports of a pair of molecular markers, KASP_IWB72471 and KASP_IWB10558, being predictive of ASR genes \u003cem\u003eSr11\u003c/em\u003e and \u003cem\u003eSr8155B1\u003c/em\u003e, respectively (Nirmala et al. \u003cspan citationid=\"CR54\" class=\"CitationRef\"\u003e2016\u003c/span\u003e; Nirmala et al. \u003cspan citationid=\"CR55\" class=\"CitationRef\"\u003e2017\u003c/span\u003e). \u003cem\u003eSr11\u003c/em\u003e, residing on chromosome 6BL, confers resistance to \u003cem\u003ePgt\u003c/em\u003e race TKTTF (Green et al. \u003cspan citationid=\"CR21\" class=\"CitationRef\"\u003e1960\u003c/span\u003e; Nirmala et al. \u003cspan citationid=\"CR54\" class=\"CitationRef\"\u003e2016\u003c/span\u003e). However, the CItr 11390 was susceptible to TKTTF despite carrying the allele at KASP_IWB72471 reported to be predictive of \u003cem\u003eSr11\u003c/em\u003e resistance. Additionally, KASP_IWB72471 was found to have no effect on seedling resistance to TKTTF in single marker analysis, and \u003cem\u003eSr11\u003c/em\u003e was not detected by MQM for TKTTF or APR resistance despite substantial marker coverage in the region. Similarly, CItr 11390 was found to carry the allele reported to be predictive of \u003cem\u003eSr8155B1\u003c/em\u003e resistance at KASP_IWB10558 on chromosome 6AS. However, CItr 11390, as well as MN07098-6 and the entire RIL population, were susceptible to race TTKTT, which is avirulent to Sr8155B1. These results indicate that KASP_IWB72471 and KASP_IWB10558 are not predictive of resistance at their respective loci in this population.\u003c/p\u003e \u003cp\u003eInitial analysis of the phenotypic data from seedling testing to race TKTTF indicated resistance was conferred by a single gene (Table\u0026nbsp;\u003cspan refid=\"Tab3\" class=\"InternalRef\"\u003e3\u003c/span\u003e). This was expected given the population segregates for \u003cem\u003eSr7a\u003c/em\u003e. However, MQM detected a second locus stretching approximately 50Mb on the long arm of chromosome 5A. Like \u003cem\u003eSr7a\u003c/em\u003e, the resistance allele at \u003cem\u003eQSr.umn-5A.1\u003c/em\u003e was donated by the common wheat parent, MN07098-6. No known ASR genes are located on chromosome 5AL. Association mapping of North American wheat breeding germplasm by Bajgain et al. (\u003cspan citationid=\"CR1\" class=\"CitationRef\"\u003e2015a\u003c/span\u003e) detected three significant MTAs associated with seedling resistance to TKTTF on chromosome 5AL. Each of the three SNPs fall within the 95% Bayes credible interval of \u003cem\u003eQSr.umn-5A.1\u003c/em\u003e, and one is located less than 3Mb from the peak of \u003cem\u003eQSr.umn-5A.1\u003c/em\u003e. MN07098-6 was included among the association mapping panel used by Bajgain et al. (\u003cspan citationid=\"CR1\" class=\"CitationRef\"\u003e2015a\u003c/span\u003e), and the line carries the resistance allele at each of the three SNPs. The detection of co-locating QTL and MTAs sourced from common germplasm suggests the interval on chromosome 5AL may contain a novel source of resistance to \u003cem\u003ePgt\u003c/em\u003e race TKTTF.\u003c/p\u003e \u003cp\u003eMultiple QTL mapping detected a total of seven QTL for APR to \u003cem\u003ePgt\u003c/em\u003e in three of the four test environments. Among these QTL, two were detected in multiple environments that were both donated by CItr 11390. The largest effect QTL of the three, \u003cem\u003eQSr.umn-2A\u003c/em\u003e, was detected in Kenya in 2019 and 2020. It is likely \u003cem\u003eQSr.umn-2A\u003c/em\u003e is \u003cem\u003eSr63\u003c/em\u003e, an APR locus first reported by Mago et al. (\u003cspan citationid=\"CR43\" class=\"CitationRef\"\u003e2022\u003c/span\u003e) from the durum wheat (\u003cem\u003eTriticum turgidum\u003c/em\u003e ssp. \u003cem\u003edurum\u003c/em\u003e) variety \u0026lsquo;Glossy Huguenot\u0026rsquo;. A KASP marker linked to \u003cem\u003eSr63\u003c/em\u003e, KASP_IWB32429, was located at the peak of \u003cem\u003eQSr.umn-2A\u003c/em\u003e in both environments. Additionally, the physical positions of the \u003cem\u003eQSr.umn-2A\u003c/em\u003e interval (681\u0026ndash;702 Mb) were nearly identical to that of \u003cem\u003eSr63\u003c/em\u003e reported by Mago el al. (2022) in the RefSeq v2.1 assembly (683\u0026ndash;696 Mb). The detected interval is distinct from the location of \u003cem\u003eSr21\u003c/em\u003e (713 Mb), and molecular marker screening confirmed the lack of the resistance allele at \u003cem\u003eSr21\u003c/em\u003e in CItr 11390. The detection of \u003cem\u003eQSr.umn-2A\u003c/em\u003e and \u003cem\u003eSr63\u003c/em\u003e in a common environment (Kenya 2019) further supports the validation of \u003cem\u003eSr63\u003c/em\u003e as a unique APR locus on chromosome 2AL.\u003c/p\u003e \u003cp\u003eGiven the population\u0026rsquo;s apparent seedling susceptibility to Ug99 lineage races (Table\u0026nbsp;\u003cspan refid=\"Tab2\" class=\"InternalRef\"\u003e2\u003c/span\u003e), \u003cem\u003eSr63\u003c/em\u003e on chromosome 2AL could be an important target for resistance breeding to the Ug99 race group given its relatively large effect on disease reduction at the adult growth stages. The lack of detection of \u003cem\u003eSr63\u003c/em\u003e in Ethiopia indicates the gene may not be effective against other races present in Ethiopia, such as TRTTF and JRCQC. Screening the population in single-race nurseries could provide greater detail on efficacy of the gene to specific races. Additionally, APR genes have been found to provide pleiotropic resistance to the other cereal rusts, as well as other plant pathogens (Herrera-Foessel et al. \u003cspan citationid=\"CR24\" class=\"CitationRef\"\u003e2014\u003c/span\u003e; Hiebert et al. \u003cspan citationid=\"CR26\" class=\"CitationRef\"\u003e2010\u003c/span\u003e; Krattinger et al. \u003cspan citationid=\"CR35\" class=\"CitationRef\"\u003e2009\u003c/span\u003e; Martinez et al. 2004; Rinaldo et al. \u003cspan citationid=\"CR69\" class=\"CitationRef\"\u003e2017\u003c/span\u003e). Further disease resistance characterization of the 15xR012 RIL population to other pathogens, namely stripe rust (caused by \u003cem\u003eP. striiformis\u003c/em\u003e f. sp. \u003cem\u003etritici\u003c/em\u003e) and leaf rust (caused by \u003cem\u003eP. triticina\u003c/em\u003e f. sp. \u003cem\u003etritici\u003c/em\u003e) may be worth pursuing.\u003c/p\u003e \u003cp\u003e \u003cem\u003eQSr.umn-6BL\u003c/em\u003e on the long arm of chromosome 6B was also detected in more than one environment. While the QTL was not retained in the multiple QTL model for Ken19, it was initially detected by both single-QTL and two-dimensional genome scans for Ken19, Ken20, and Eth20. The ASR gene \u003cem\u003eSr11\u003c/em\u003e resides on 6BL (Green et al. \u003cspan citationid=\"CR21\" class=\"CitationRef\"\u003e1960\u003c/span\u003e), but molecular marker screening suggests that neither MN07098-6 nor CItr 11390 carry the resistance allele at \u003cem\u003eSr11\u003c/em\u003e. Additionally, the nearest flanking marker of \u003cem\u003eQSr.umn-6BL\u003c/em\u003e was more than 40 Mb from \u003cem\u003eSr11\u003c/em\u003e. Bajgain et al. (\u003cspan citationid=\"CR1\" class=\"CitationRef\"\u003e2015a\u003c/span\u003e) previously reported a marker trait association (MTA) on chromosome 6BL at IWB35697 for APR in both Kenya and Ethiopia. MTAs in close proximity to IWB35697 were also detected in winter and durum wheat diversity panels in multiple years in Kenya and Ethiopia (Megeressa et al. 2020; Yu et al. \u003cspan citationid=\"CR93\" class=\"CitationRef\"\u003e2012\u003c/span\u003e). However, local alignment of the significant markers from each of the aforementioned studies, including IWB35697, to RefSeq v2.1 revealed they were outside of the interval of \u003cem\u003eQSr.umn-6BL\u003c/em\u003e and close to \u003cem\u003eSr11\u003c/em\u003e. The distance from \u003cem\u003eSr11\u003c/em\u003e, susceptibility of CItr 11390 to TKTTF, and the lack of APR QTL reported in this region suggest \u003cem\u003eQSr.umn-6BL\u003c/em\u003e is a novel APR locus (Yu et al. \u003cspan citationid=\"CR92\" class=\"CitationRef\"\u003e2014\u003c/span\u003e). Given the potential novelty of \u003cem\u003eQSr.umn-6BL\u003c/em\u003e, refinement of the interval and the development of near isogenic lines (NILs) should be pursued to determine whether the QTL is a valuable breeding target for stem rust adult plant resistance. Efforts to do so can be initiated using the genotyping assays developed in this study (Table S7). Additional assays that more closely approximate the intervals detected during MQM are needed.\u003c/p\u003e \u003cp\u003eIn addition to the two QTL identified in multiple environments, five QTL were detected in single environments. Of these, four were contributed by the hexaploid parent MN07098-6: \u003cem\u003eQSr.umn-1A\u003c/em\u003e (Eth20), \u003cem\u003eQSr.umn-3A\u003c/em\u003e (Ken20), \u003cem\u003eQSr.umn-3B\u003c/em\u003e (Ken20), and \u003cem\u003eQSr.umn-5A.2\u003c/em\u003e (Ken19). \u003cem\u003eQSr.umn-1A\u003c/em\u003e explained a large portion of the phenotypic variation observed in Ethiopia in 2020 (\u003cem\u003eR\u003c/em\u003e\u003csup\u003e\u003cem\u003e2\u003c/em\u003e\u003c/sup\u003e\u0026thinsp;=\u0026thinsp;21.5). However, this is likely inflated by the small population size, relatively little phenotypic variation in the Eth20 environment, the lack of effectiveness of other loci such as \u003cem\u003eSr63\u003c/em\u003e, and the absence of ASR loci in the population besides \u003cem\u003eSr7a\u003c/em\u003e and \u003cem\u003eQSr.umn-5A.1\u003c/em\u003e. While no known APR gene has been reported on 1AS, other QTL within the first\u0026thinsp;~\u0026thinsp;16 Mb of chromosome 1AS have been detected in multiple studies (Bajgain et al. \u003cspan citationid=\"CR2\" class=\"CitationRef\"\u003e2015b\u003c/span\u003e; Bansal et al. \u003cspan citationid=\"CR4\" class=\"CitationRef\"\u003e2008\u003c/span\u003e; Bhavani et al. \u003cspan citationid=\"CR6\" class=\"CitationRef\"\u003e2011\u003c/span\u003e; Yu et al. \u003cspan citationid=\"CR93\" class=\"CitationRef\"\u003e2012\u003c/span\u003e). Notably, Bajgain et al. (\u003cspan citationid=\"CR2\" class=\"CitationRef\"\u003e2015b\u003c/span\u003e) detected a QTL in another University of Minnesota hard red spring wheat line for field resistance in 2013 in Kenya that co-locates with \u003cem\u003eQSr.umn-1A\u003c/em\u003e.\u003c/p\u003e \u003cp\u003e \u003cem\u003eQSr.umn-3A\u003c/em\u003e (\u003cem\u003eR\u003c/em\u003e\u003csup\u003e\u003cem\u003e2\u003c/em\u003e\u003c/sup\u003e\u0026thinsp;=\u0026thinsp;5.0) was detected in Ken20 on the long arm of chromosome 3A with resistance conferred by the allele from MN07098-6. \u003cem\u003eSr35\u003c/em\u003e is the most well-known stem rust resistance loci on chromosome 3AL (Zhang et al. \u003cspan citationid=\"CR95\" class=\"CitationRef\"\u003e2010\u003c/span\u003e). However, it is unlikely that \u003cem\u003eQSr.umn-3A\u003c/em\u003e is \u003cem\u003eSr35\u003c/em\u003e as MN07098-6 is susceptible to TTKSK, the nearest flanking marker of \u003cem\u003eQSr.umn-3A\u003c/em\u003e is ~\u0026thinsp;11Mb from \u003cem\u003eSr35\u003c/em\u003e, and the pedigree of MN07098-6 lacks \u003cem\u003eT. monococcum\u003c/em\u003e introgressions. \u003cem\u003eQSr.umn-3A\u003c/em\u003e is also unlikely to be \u003cem\u003eSr27\u003c/em\u003e, which is carried on a 3AL.3RS translocation derived from the rye (\u003cem\u003eSecale cereale\u003c/em\u003e L.) variety \u0026lsquo;Imperial\u0026rsquo; (Marais \u003cspan citationid=\"CR45\" class=\"CitationRef\"\u003e2001\u003c/span\u003e). \u003cem\u003eSr27\u003c/em\u003e provides strong seedling resistance to TTKSK that likely would have been detected during screening. Additionally, the deployment of the 3AL.3RS translocation has been almost exclusive to triticale (Singh et al. \u003cspan citationid=\"CR78\" class=\"CitationRef\"\u003e2011\u003c/span\u003e). An MTA within 5 Mb of \u003cem\u003eQSr.umn-3A\u003c/em\u003e was reported by Letta et al. (\u003cspan citationid=\"CR40\" class=\"CitationRef\"\u003e2013\u003c/span\u003e). However, this association was detected in a durum wheat panel grown in an Ethiopian nursery. In common wheat, the only APR QTL reported to partially overlap with \u003cem\u003eQSr.umn-3A\u003c/em\u003e was detected in a CIMMYT biparental population grown in Njoro, Kenya during the 2011 main season (Singh et al. \u003cspan citationid=\"CR81\" class=\"CitationRef\"\u003e2013\u003c/span\u003e). Further testing is needed to verify \u003cem\u003eQSr.umn-3A\u003c/em\u003e given the limited reports of stem rust APR QTL within or near the interval.\u003c/p\u003e \u003cp\u003e \u003cem\u003eQSr.umn-3B\u003c/em\u003e (\u003cem\u003eR\u003c/em\u003e\u003csup\u003e\u003cem\u003e2\u003c/em\u003e\u003c/sup\u003e\u0026thinsp;=\u0026thinsp;6.6) was detected in Ken20 in the centromeric region of the long arm of chromosome 3B. The interval is in close proximity to a KASP marker linked to an NB-LRR motif within the \u003cem\u003eSr12\u003c/em\u003e locus that co-segregates with seedling stage resistance (Hiebert et al. \u003cspan citationid=\"CR25\" class=\"CitationRef\"\u003e2016\u003c/span\u003e). The population is fixed for the allele linked to resistance at this \u003cem\u003eSr12\u003c/em\u003e (Table\u0026nbsp;\u003cspan refid=\"Tab1\" class=\"InternalRef\"\u003e1\u003c/span\u003e). However, adult plant resistance to Ug99 lineage and North American \u003cem\u003ePgt\u003c/em\u003e races also co-segregates with \u003cem\u003eSr12\u003c/em\u003e (Hiebert et al. \u003cspan citationid=\"CR25\" class=\"CitationRef\"\u003e2016\u003c/span\u003e; Rouse et al. \u003cspan citationid=\"CR73\" class=\"CitationRef\"\u003e2014\u003c/span\u003e). The co-localization of \u003cem\u003eQSr.umn-3B\u003c/em\u003e with the aforementioned loci suggests it may be the same source of APR. Hiebert et al. (\u003cspan citationid=\"CR25\" class=\"CitationRef\"\u003e2016\u003c/span\u003e) also determined that the APR locus within \u003cem\u003eSr12\u003c/em\u003e interval interacted with \u003cem\u003eSr57/Lr34/Yr18\u003c/em\u003e. Due to the inability to include \u003cem\u003eSr57\u003c/em\u003e/\u003cem\u003eLr34\u003c/em\u003e/\u003cem\u003eYr18\u003c/em\u003e in QTL mapping, such interactions could not be evaluated while fitting multiple QTL models. Development and testing of additional genetic resource, such as families of NILs carrying various allelic combinations at the two loci, is necessary to further evaluate the effect of \u003cem\u003eQSr.umn-3B\u003c/em\u003e and determine whether the locus interacts with \u003cem\u003eSr57/Lr34/Yr18\u003c/em\u003e.\u003c/p\u003e \u003cp\u003e \u003cem\u003eQSr.umn-5A.2\u003c/em\u003e (\u003cem\u003eR\u003c/em\u003e\u003csup\u003e\u003cem\u003e2\u003c/em\u003e\u003c/sup\u003e\u0026thinsp;=\u0026thinsp;8.6) was detected in Ken19 on the long arm of chromosome 5A. During the 2013 off-season in Ethiopia, Bajgain et al. (\u003cspan citationid=\"CR1\" class=\"CitationRef\"\u003e2015a\u003c/span\u003e) detected a single MTA associated with APR that is ~\u0026thinsp;13Mb from the peak marker of \u003cem\u003eQSr.umn-5A\u003c/em\u003e lines. QTL on the chromosome 5AL associated with APR in African stem rust nurseries have also been reported in CIMMYT germplasm (Bhavani et al. \u003cspan citationid=\"CR6\" class=\"CitationRef\"\u003e2011\u003c/span\u003e). Additionally, Edae et al. (\u003cspan citationid=\"CR12\" class=\"CitationRef\"\u003e2018\u003c/span\u003e) detected MTAs near the MTA reported by Bajgain et al. (\u003cspan citationid=\"CR1\" class=\"CitationRef\"\u003e2015a\u003c/span\u003e) for both infection response and disease severity to race QTHJC across three growing seasons in Minnesota using the same panel of spring wheat lines. As previously discussed, a QTL providing ASR to TKTTF was also detected in this study that co-localized with \u003cem\u003eQSr.umn-5A.1\u003c/em\u003e. In the case of both QTL, the resistance allele was provided by MN07098-6. Given \u003cem\u003eQSr.umn.5A.1\u003c/em\u003e was only detected in a single environment, further testing is needed to validate the QTL and also determine the mechanism for co-localization of ASR and APR much like the \u003cem\u003eSr12\u003c/em\u003e locus.\u003c/p\u003e \u003cp\u003eThe last QTL detected in this study, \u003cem\u003eQSr.umn-6BS\u003c/em\u003e (\u003cem\u003eR\u003c/em\u003e\u003csup\u003e\u003cem\u003e2\u003c/em\u003e\u003c/sup\u003e\u0026thinsp;=\u0026thinsp;16.2), was detected only in Ken20 and contributed by CItr 11390. The effect of \u003cem\u003eQSr.umn-6BS\u003c/em\u003e was relative small compared to the other QTL detected in Ken20. Additionally, an interaction between this QTL and \u003cem\u003eSr63\u003c/em\u003e was detected during MQM. The mean CI value for lines carrying the resistance allele at both \u003cem\u003eQSr.umn-6BS\u003c/em\u003e and \u003cem\u003eSr63\u003c/em\u003e was not significantly different from those only carrying the resistance allele at \u003cem\u003eSr63\u003c/em\u003e. This suggests that the presence of the resistance allele at \u003cem\u003eSr63\u003c/em\u003e masks the effect of \u003cem\u003eQSr.umn-6BS\u003c/em\u003e and is sufficient to reduce disease symptoms. Given the observed masking effect, as well as the substantially larger effect and detection in multiple environments of \u003cem\u003eSr63\u003c/em\u003e, there is little reason to prioritize further investigations and selection of \u003cem\u003eQSr.umn-6BS\u003c/em\u003e.\u003c/p\u003e \u003cp\u003eStacking multiple APR genes is an effective approach to developing high levels of resistance to stem rust (Pretorius et al. \u003cspan citationid=\"CR65\" class=\"CitationRef\"\u003e2017\u003c/span\u003e). The mean CI of RILs carrying three resistance alleles at \u003cem\u003eSr57/Lr34/Yr18, Sr63\u003c/em\u003e, and \u003cem\u003eQSr.umn-6BL\u003c/em\u003e was markedly lower across environments compared to those carrying resistance alleles at zero, one, or two of the loci (Fig.\u0026nbsp;\u003cspan refid=\"Fig5\" class=\"InternalRef\"\u003e5\u003c/span\u003e, Table S6). However, developing pyramids of favorable stem rust resistance alleles from alien sources that are accessible for commercial breeding is challenging due to issues associated with interspecific hybridization such as linkage drag and chromosome pairing.\u003c/p\u003e \u003cp\u003eSources for introgression are limited by the frequency of the resistance allele within a donor species. For example, of the 155 CIMMYT durum lines Mago et al. (\u003cspan citationid=\"CR43\" class=\"CitationRef\"\u003e2022\u003c/span\u003e) screened for \u003cem\u003eSr63\u003c/em\u003e, two were found to carry the resistance allele, making CItr 11390 only the third wheat accession known to carry the resistance allele at this locus. The RILs of the 15xR012 population are therefore the first hexaploid wheat lines known to carry \u003cem\u003eSr63\u003c/em\u003e. This is also likely the case for \u003cem\u003eQSr.umn-6BL\u003c/em\u003e, as this is the first report of QTL for stem rust APR within this region. Multiple euploid RILs carry the resistance allele at \u003cem\u003eSr7a\u003c/em\u003e and \u003cem\u003eSr57/Lr34/Yr18\u003c/em\u003e in addition to \u003cem\u003eSr63\u003c/em\u003e and \u003cem\u003eQSr.umn-6BL\u003c/em\u003e (Table S9).\u003c/p\u003e \u003cp\u003eThe development and discovery of these lines demonstrates the feasibility of simultaneously mapping and transferring QTL from Khorasan wheat into hexaploid breeding materials. The RILs and SNP genotyping assays developed in this study provide value not only to targeted breeding efforts in areas where Ug99 stem rust races are present, but also facilitate the rapid deployment of newly validated and discovered loci such as \u003cem\u003eSr63\u003c/em\u003e and \u003cem\u003eQSr.umn-6BL\u003c/em\u003e to the greater wheat breeding community.\u003c/p\u003e"},{"header":"Declarations","content":"\u003ch2\u003eFunding\u003c/h2\u003e \u003cp\u003eThis project was supported by the USDA National Institute of Food and Agriculture competitive grant no. 2022-68013-36439, USDA-ARS appropriated project \u0026ldquo;Surveillance, Pathogen Biology, and Host Resistance of Cereal Rusts\u0026rdquo;, and the USDA-ARS National Plant Disease Recovery System.\u003c/p\u003e\u003ch2\u003eAuthor Contributions\u003c/h2\u003e \u003cp\u003eMF coordinated sequencing, constructed the genetic map, performed data analysis, and wrote the final manuscript. EC conducted molecular marker genotyping, aided in sequencing library preparation, and contributed to data analysis. ZK and AG coordinated field trials and oversaw data collection. MR conceived the study, developed the 15xR012 mapping population, facilitated seedling disease screenings, and aided in trial design. JA contributed to data analysis and provided suggestions to experimental design. All authors reviewed the manuscript.\u003c/p\u003e\u003ch2\u003eAcknowledgments\u003c/h2\u003e \u003cp\u003eThis project was supported by the USDA National Institute of Food and Agriculture competitive grant no. 2022-68013-36439, USDA-ARS appropriated project \u0026ldquo;Surveillance, Pathogen Biology, and Host Resistance of Cereal Rusts\u0026rdquo;, and the USDA-ARS National Plant Disease Recovery System.\u003c/p\u003e\u003ch2\u003eData Availability\u003c/h2\u003e \u003cp\u003ePhenotypic and genotypic data used for analysis are included in this article and supplementary files.\u003c/p\u003e"},{"header":"References","content":"\u003col\u003e\n\u003cli\u003eBajgain P, Rouse M, Bulli P, et al (2015a) Association mapping of North American spring wheat breeding germplasm reveals loci conferring resistance to Ug99 and other African stem rust races. BMC Plant Biology 15\u003c/li\u003e\n\u003cli\u003eBajgain P, Rouse MN, Bhavani S, Anderson JA (2015b) QTL mapping of adult plant resistance to Ug99 stem rust in the spring wheat population RB07/MN06113-8. Molecular Breeding 35\u003c/li\u003e\n\u003cli\u003eBansal U, Bariana H, Wong D, et al (2014) Molecular mapping of an adult plant stem rust resistance gene Sr56 in winter wheat cultivar Arina. Theoretical and Applied Genetics 127:1441\u0026ndash;1448\u003c/li\u003e\n\u003cli\u003eBansal UK, Bossolini E, Miah H, et al (2008) Genetic mapping of seedling and adult plant stem rust resistance in two European winter wheat cultivars. Euphytica 164:821\u0026ndash;828\u003c/li\u003e\n\u003cli\u003eBates D, M\u0026auml;chler M, Bolker B, Walker S (2015) Fitting Linear Mixed-Effects Models Using lme4. Journal of Statistical Software 67:1\u0026ndash;48\u003c/li\u003e\n\u003cli\u003eBhavani S, Singh R, Argillier O, et al (2011) Mapping durable adult plant stem rust resistance to the race Ug99 group in six CIMMYT wheats. 2011 Borlaug Global Rust Initiative Technical Workshop 43\u0026ndash;53\u003c/li\u003e\n\u003cli\u003eBradbury PJ, Zhang Z, Kroon DE, et al (2007) TASSEL: software for association mapping of complex traits in diverse samples. Bioinformatics 23:2633\u0026ndash;2635\u003c/li\u003e\n\u003cli\u003eBraun HJ, Atlin G, Payne T (2010) Multi-location testing as a tool to identify plant response to global climate change\u003c/li\u003e\n\u003cli\u003eBriggs J, Chen S, Zhang W, et al (2015) Mapping of \u003cem\u003eSrTm4\u003c/em\u003e, a Recessive Stem Rust Resistance Gene from Diploid Wheat Effective to Ug99. Phytopathology\u0026reg; 105:1347\u0026ndash;1354\u003c/li\u003e\n\u003cli\u003eBroman KW, Wu H, Sen Ś, Churchill GA (2003) R/qtl: QTL mapping in experimental crosses. Bioinformatics 19:889\u0026ndash;890\u003c/li\u003e\n\u003cli\u003eChen S, Zhang W, Bolus S, et al (2018) Identification and characterization of wheat stem rust resistance gene Sr21 effective against the Ug99 race group at high temperature. PLOS Genetics 14:e1007287\u003c/li\u003e\n\u003cli\u003eEdae EA, Pumphrey MO, Rouse MN (2018) A Genome-Wide Association Study of Field and Seedling Response to Individual Stem Rust Pathogen Races Reveals Combinations of Race-Specific Genes in North American Spring Wheat. Frontiers in Plant Science 9\u003c/li\u003e\n\u003cli\u003eElshire RJ, Glaubitz JC, Sun Q, et al (2011) A Robust, Simple Genotyping-by-Sequencing (GBS) Approach for High Diversity Species. PLoS ONE 6:e19379\u003c/li\u003e\n\u003cli\u003eFAO (2022) Crop Prospects and Food Situation. Food and Agriculture Organization of the United Nations\u003c/li\u003e\n\u003cli\u003eGao L, Babiker EM, Nava IC, et al (2019) Temperature‐sensitive wheat stem rust resistance gene Sr15 is effective against Puccinia graminis f. sp. tritici race TTKSK. Plant Pathology 68:143\u0026ndash;151\u003c/li\u003e\n\u003cli\u003eGao L, Turner MK, Chao S, et al (2016) Genome Wide Association Study of Seedling and Adult Plant Leaf Rust Resistance in Elite Spring Wheat Breeding Lines. PLOS ONE 11:e0148671\u003c/li\u003e\n\u003cli\u003eGill BS, Friebe B, Endo TR (1991) Standard karyotype and nomenclature system for description of chromosome bands and structural aberrations in wheat (Triticum aestivum). Genome 34:830\u0026ndash;839\u003c/li\u003e\n\u003cli\u003eGill BS, Raupp WJ (1987) Direct Genetic Transfers from Aegilops squarrosa L. to Hexaploid Wheat. Crop Science 27:445\u0026ndash;450\u003c/li\u003e\n\u003cli\u003eGodfray HCJ, Beddington JR, Crute IR, et al (2010) Food Security: The Challenge of Feeding 9 Billion People. Science 327:812\u0026ndash;818\u003c/li\u003e\n\u003cli\u003eGrassini P, Eskridge KM, Cassman KG (2013) Distinguishing between yield advances and yield plateaus in historical crop production trends. Nature Communications 4\u003c/li\u003e\n\u003cli\u003eGreen GJ, Knott DR, Watson IA, Pugsley AT (1960) Seedling reactions to stem rust of lines of Marquis wheat with substituted genes for rust resistance. Canadian Journal of Plant Science 40:524\u0026ndash;538\u003c/li\u003e\n\u003cli\u003eHayes HK, Ausemus ER, Stakman EC, et al (1936) Thatcher Wheat\u003c/li\u003e\n\u003cli\u003eHelguera M, Khan IA, Kolmer J, et al (2003) PCR Assays for the \u003cem\u003eLr37‐Yr17‐Sr38\u003c/em\u003e Cluster of Rust Resistance Genes and Their Use to Develop Isogenic Hard Red Spring Wheat Lines. Crop Science 43:1839\u0026ndash;1847\u003c/li\u003e\n\u003cli\u003eHerrera-Foessel SA, Singh RP, Lillemo M, et al (2014) Lr67/Yr46 confers adult plant resistance to stem rust and powdery mildew in wheat. Theoretical and Applied Genetics 127:781\u0026ndash;789\u003c/li\u003e\n\u003cli\u003eHiebert CW, Kolmer JA, Mccartney CA, et al (2016) Major Gene for Field Stem Rust Resistance Co-Locates with Resistance Gene Sr12 in \u0026lsquo;Thatcher\u0026rsquo; Wheat. PLOS ONE 11:e0157029\u003c/li\u003e\n\u003cli\u003eHiebert CW, Thomas JB, Mccallum BD, et al (2010) An introgression on wheat chromosome 4DL in RL6077 (Thatcher*6/PI 250413) confers adult plant resistance to stripe rust and leaf rust (Lr67). Theoretical and Applied Genetics 121:1083\u0026ndash;1091\u003c/li\u003e\n\u003cli\u003eKalous JR, Martin JM, Sherman JD, et al (2015) Impact of the D genome and quantitative trait loci on quantitative traits in a spring durum by spring bread wheat cross. Theoretical and Applied Genetics 128:1799\u0026ndash;1811\u003c/li\u003e\n\u003cli\u003eKimber G, Sears ER (1968) Nomenclature for the description of aneuploids in the Triticinae. 3rd International Wheat Genetics Symposium 468\u0026ndash;473\u003c/li\u003e\n\u003cli\u003eKislev ME (1982) Stem Rust of Wheat 3300 Years Old Found in Israel. Science 216:993\u0026ndash;994\u003c/li\u003e\n\u003cli\u003eKnott DR (1962) The inheritance of rust resistance: IX. The inheritance of resistance to races 15B and 56 of stem rust in the wheat variety Khapstein. Canadian Journal of Plant Science 42:415\u0026ndash;419\u003c/li\u003e\n\u003cli\u003eKnott DR (2000) The inheritance of stem rust resistance in Thatcher wheat. Canadian Journal of Plant Science 80:53\u0026ndash;63\u003c/li\u003e\n\u003cli\u003eKolmer JA (1996) Genetics of resistance to wheat leaf rust. Annual Review of Phytopathology 34:435\u0026ndash;455\u003c/li\u003e\n\u003cli\u003eKosambi DD (1943) The estimation of map distances from recombination values. Annals of Eugenics 12:172\u0026ndash;175\u003c/li\u003e\n\u003cli\u003eKosgey ZC, Edae EA, Dill-Macky R, et al (2021) Mapping and Validation of Stem Rust Resistance Loci in Spring Wheat Line CI 14275. Frontiers in Plant Science 11\u003c/li\u003e\n\u003cli\u003eKrattinger SG, Lagudah ES, Spielmeyer W, et al (2009) A Putative ABC Transporter Confers Durable Resistance to Multiple Fungal Pathogens in Wheat. Science 323:1360\u0026ndash;1363\u003c/li\u003e\n\u003cli\u003eKumar S, Fetch TG, Knox RE, et al (2021) Mapping of Ug99 stem rust resistance in Canadian durum wheat. Canadian Journal of Plant Pathology 43:599\u0026ndash;611\u003c/li\u003e\n\u003cli\u003eLanning SP, Blake NK, Sherman JD, Talbert LE (2008) Variable Production of Tetraploid and Hexaploid Progeny Lines from Spring Wheat by Durum Wheat Crosses. Crop Science 48:199\u0026ndash;202\u003c/li\u003e\n\u003cli\u003eLaugerotte J, Baumann U, Sourdille P (2022) Genetic control of compatibility in crosses between wheat and its wild or cultivated relatives. Plant Biotechnology Journal 20:812\u0026ndash;832\u003c/li\u003e\n\u003cli\u003eLeigh FJ, Wright TIC, Horsnell RA, et al (2022) Progenitor species hold untapped diversity for potential climate-responsive traits for use in wheat breeding and crop improvement. Heredity 128:291\u0026ndash;303\u003c/li\u003e\n\u003cli\u003eLetta T, Maccaferri M, Badebo A, et al (2013) Searching for novel sources of field resistance to Ug99 and Ethiopian stem rust races in durum wheat via association mapping. Theoretical and Applied Genetics 126:1237\u0026ndash;1256\u003c/li\u003e\n\u003cli\u003eLi H (2011) A statistical framework for SNP calling, mutation discovery, association mapping and population genetical parameter estimation from sequencing data. Bioinformatics 27:2987\u0026ndash;2993\u003c/li\u003e\n\u003cli\u003eLi H, Handsaker B, Wysoker A, et al (2009) The Sequence Alignment/Map format and SAMtools. Bioinformatics 25:2078\u0026ndash;2079\u003c/li\u003e\n\u003cli\u003eMago R, Chen C, Xia X, et al (2022) Adult plant stem rust resistance in durum wheat Glossy Huguenot: mapping, marker development and validation. Theoretical and Applied Genetics 135:1541\u0026ndash;1550\u003c/li\u003e\n\u003cli\u003eManichaikul A, Moon JY, Sen S, et al (2009) A Model Selection Approach for the Identification of Quantitative Trait Loci in Experimental Crosses, Allowing Epistasi. Genetics 181:1077\u0026ndash;1086\u003c/li\u003e\n\u003cli\u003eMarais GF (2001) An evaluation of three Sr27-carrying wheat \u0026times; rye translocations. South African Journal of Plant and Soil 18:135\u0026ndash;136\u003c/li\u003e\n\u003cli\u003eMartin M (2011) Cutadapt removes adapter sequences from high-throughput sequencing reads. EMBnet.journal 17:10\u003c/li\u003e\n\u003cli\u003eMart\u0026iacute;nez F, Niks RE, Singh RP, Rubiales D (2001) Characterization of Lr46, a Gene Conferring Partial Resistance to Wheat Leaf Rust. Hereditas 135:111\u0026ndash;114\u003c/li\u003e\n\u003cli\u003eMcFadden ES (1930) A successful transfer of emmer characters to vulgare wheat. Journal of the American Society of Agronomy 22:1020\u0026ndash;1034\u003c/li\u003e\n\u003cli\u003eMcIntosh RA, Dubcovsky J, Rogers WJ, et al (2020) Catalogue of Gene Symbols for Wheat: 2020 Supplement. Annual Wheat Newsletter 66\u003c/li\u003e\n\u003cli\u003eMcIntosh RA, Wellings CR, Park RF (1995) Wheat Rusts: An Atlas of Resistance Genes\u003c/li\u003e\n\u003cli\u003eMegerssa SH, Ammar K, Acevedo M, et al (2020) Multiple-Race Stem Rust Resistance Loci Identified in Durum Wheat Using Genome-Wide Association Mapping. Frontiers in Plant Science 11\u003c/li\u003e\n\u003cli\u003eMoney D, Gardner K, Migicovsky Z, et al (2015) LinkImpute: Fast and Accurate Genotype Imputation for Nonmodel Organisms. G3 Genes|Genomes|Genetics 5:2383\u0026ndash;2390\u003c/li\u003e\n\u003cli\u003eMustahsan W (2022) Characterization of stripe rust resistance and yield-related traits in an \u0026lsquo;Overley\u0026rsquo; x \u0026lsquo;Overland\u0026rsquo; population. Thesis\u003c/li\u003e\n\u003cli\u003eNirmala J, Chao S, Olivera P, et al (2016) Markers Linked to Wheat Stem Rust Resistance Gene Sr11 Effective to Puccinia graminis f. sp. tritici Race TKTTF. Phytopathology 106:1352\u0026ndash;1358\u003c/li\u003e\n\u003cli\u003eNirmala J, Saini J, Newcomb M, et al (2017) Discovery of a Novel Stem Rust Resistance Allele in Durum Wheat that Exhibits Differential Reactions to Ug99 Isolates. G3 Genes|Genomes|Genetics 7:3481\u0026ndash;3490\u003c/li\u003e\n\u003cli\u003eOlivera P, Newcomb M, Szabo LJ, et al (2015) Phenotypic and Genotypic Characterization of Race TKTTF of Puccinia graminis f. sp. tritici that Caused a Wheat Stem Rust Epidemic in Southern Ethiopia in 2013\u0026ndash;14. Phytopathology\u0026reg; 105:917\u0026ndash;928\u003c/li\u003e\n\u003cli\u003eOlivera PD, Jin Y, Rouse M, et al (2012) Races of Puccinia graminis f. sp. tritici with Combined Virulence to Sr13 and Sr9e in a Field Stem Rust Screening Nursery in Ethiopia. Plant Disease 96:623\u0026ndash;628\u003c/li\u003e\n\u003cli\u003eOlivera PD, Villegas D, Cantero‐Mart\u0026iacute;nez C, et al (2022) A unique race of the wheat stem rust pathogen with virulence on \u003cem\u003eSr31\u003c/em\u003e identified in Spain and reaction of wheat and durum cultivars to this race. Plant Pathology 71:873\u0026ndash;889\u003c/li\u003e\n\u003cli\u003eOlson EL, Rouse MN, Pumphrey MO, et al (2013a) Introgression of stem rust resistance genes SrTA10187 and SrTA10171 from Aegilops tauschii to wheat. Theoretical and Applied Genetics 126:2477\u0026ndash;2484\u003c/li\u003e\n\u003cli\u003eOlson EL, Rouse MN, Pumphrey MO, et al (2013b) Simultaneous transfer, introgression, and genomic localization of genes for resistance to stem rust race TTKSK (Ug99) from Aegilops tauschii to wheat. Theoretical and Applied Genetics 126:1179\u0026ndash;1188\u003c/li\u003e\n\u003cli\u003eOr\u0026oacute;stica KY, Verdugo RA (2016) chromPlot: visualization of genomic data in chromosomal context. Bioinformatics 32:2366\u0026ndash;2368\u003c/li\u003e\n\u003cli\u003ePatpour M, Justesen AF, Tecle AW, et al (2020) First Report of Race TTRTF of Wheat Stem Rust (Puccinia graminis f. sp. tritici) in Eritrea. Plant Disease 104:973\u003c/li\u003e\n\u003cli\u003ePeterson RF, Campbell AB, Hannah AE (1948) A diagrammatic scale for estimating rust intensity on leaves and stems of cereals. Canadian Journal of Research 26c:496\u0026ndash;500\u003c/li\u003e\n\u003cli\u003ePoland JA, Brown PJ, Sorrells ME, Jannink J-L (2012) Development of High-Density Genetic Maps for Barley and Wheat Using a Novel Two-Enzyme Genotyping-by-Sequencing Approach. PLoS ONE 7:e32253\u003c/li\u003e\n\u003cli\u003ePretorius ZA, Ayliffe M, Bowden RL, et al (2017) Advances in control of wheat rusts. Achieving sustainable cultivation of wheat 1:295\u0026ndash;343\u003c/li\u003e\n\u003cli\u003ePretorius ZA, Jin Y, Bender CM, et al (2012) Seedling resistance to stem rust race Ug99 and marker analysis for Sr2, Sr24 and Sr31 in South African wheat cultivars and lines. Euphytica 186:15\u0026ndash;23\u003c/li\u003e\n\u003cli\u003ePretorius ZA, Singh RP, Wagoire WW, Payne TS (2000) Detection of Virulence to Wheat Stem Rust Resistance Gene Sr31 in Puccinia graminis f. sp. tritici in Uganda. Plant Disease 84:203\u0026ndash;203\u003c/li\u003e\n\u003cli\u003eQuinlan AR, Hall IM (2010) BEDTools: a flexible suite of utilities for comparing genomic features. Bioinformatics 26:841\u0026ndash;842\u003c/li\u003e\n\u003cli\u003eRinaldo A, Gilbert B, Boni R, et al (2017) The Lr34 adult plant rust resistance gene provides seedling resistance in durum wheat without senescence. Plant Biotechnology Journal 15:894\u0026ndash;905\u003c/li\u003e\n\u003cli\u003eRoelfs AP (1985) Wheat and rye stem rust. The Cereal Rusts, Vol. 2 2:3\u0026ndash;37\u003c/li\u003e\n\u003cli\u003eRoelfs AP, Martens JWM (1988) An international system of nomenclature for Puccinia graminis f.sp. tritici. Phytopathology 78:526\u0026ndash;533\u003c/li\u003e\n\u003cli\u003eRoelfs AP, Singh RP, Saari EE (1992) Rust Diseases of Wheat: Concepts and Methods of Disease Management.\u003c/li\u003e\n\u003cli\u003eRouse MN, Talbert LE, Singh D, Sherman JD (2014) Complementary epistasis involving Sr12 explains adult plant resistance to stem rust in Thatcher wheat (Triticum aestivum L.). Theoretical and Applied Genetics 127:1549\u0026ndash;1559\u003c/li\u003e\n\u003cli\u003eRouse MN, Wanyera R, Njau P, Jin Y (2011) Sources of Resistance to Stem Rust Race Ug99 in Spring Wheat Germplasm. Plant Disease 95:762\u0026ndash;766\u003c/li\u003e\n\u003cli\u003eRowland GG, Kerber ER (1974) Telocentric mapping in hexaploid wheat of genes for leaf rust resistance and other characters derived from Aegilops squarrosa. Canadian Journal of Genetics and Cytology 16:137\u0026ndash;144\u003c/li\u003e\n\u003cli\u003eSharma HC, Gill BS (1983) Current status of wide hybridization in wheat. Euphytica 32:17\u0026ndash;31\u003c/li\u003e\n\u003cli\u003eSingh N, Raupp J, Koo D-H, et al (2020) In-silico detection of aneuploidy and chromosomal deletions in wheat using genotyping-by-sequencing. Plant Methods 16\u003c/li\u003e\n\u003cli\u003eSingh RP, Hodson DP, Huerta-Espino J, et al (2011) The Emergence of Ug99 Races of the Stem Rust Fungus is a Threat to World Wheat Production. Annual Review of Phytopathology 49:465\u0026ndash;481\u003c/li\u003e\n\u003cli\u003eSingh RP, Hodson DP, Jin Y, et al (2015) Emergence and Spread of New Races of Wheat Stem Rust Fungus: Continued Threat to Food Security and Prospects of Genetic Control. Phytopathology\u0026reg; 105:872\u0026ndash;884\u003c/li\u003e\n\u003cli\u003eSingh RP, Huerta-Espino J (2003) Effect of leaf rust resistance gene Lr34 on components of slow rusting at seven growth stages in wheat. Euphytica 129:371\u0026ndash;376\u003c/li\u003e\n\u003cli\u003eSingh S, Singh RP, Bhavani S, et al (2013) QTL mapping of slow-rusting, adult plant resistance to race Ug99 of stem rust fungus in PBW343/Muu RIL population. Theoretical and Applied Genetics 126:1367\u0026ndash;1375\u003c/li\u003e\n\u003cli\u003eStakman EC, Stewart DM, Loegering WQ (1962) Identification of physiologic races of Puccinia graminis var. tritici. E617\u003c/li\u003e\n\u003cli\u003eStubbs RW, Prescott JM, Saari EE, Dubin HJ (1986) Cereal disease methodology manual\u003c/li\u003e\n\u003cli\u003eStuthman DD, Leonard KJ, Miller‐Garvin J (2007) Breeding Crops for Durable Resistance to Disease. 95:319\u0026ndash;367\u003c/li\u003e\n\u003cli\u003eSuenaga K, Singh RP, Huerta-Espino J, William HM (2003) Microsatellite Markers for Genes Lr34/Yr18 and Other Quantitative Trait Loci for Leaf Rust and Stripe Rust Resistance in Bread Wheat. Phytopathology 93:881\u0026ndash;890\u003c/li\u003e\n\u003cli\u003eTaylor J, Butler D (2017) R Package ASMap: Efficient Genetic Linkage Map Construction and Diagnosis. Journal of Statistical Software 79:1\u0026ndash;29\u003c/li\u003e\n\u003cli\u003eTerefe TG, Visser B, and Pretorius ZA (2016) Variation in Puccinia graminis f. sp. tritici detected on wheat and triticale in South Africa from 2009 to 2013. Crop Protection 86:9\u0026ndash;16\u003c/li\u003e\n\u003cli\u003eThe TT (1973) Chromosome Location of Genes Conditioning Stem Rust Resistance Transferred from Diploid to Hexaploid Wheat. Nature New Biology 241:256\u0026ndash;256\u003c/li\u003e\n\u003cli\u003eTurner MK, Jin Y, Rouse MN, Anderson JA (2016) Stem Rust Resistance in \u0026lsquo;Jagger\u0026rsquo; Winter Wheat. Crop Science 56:1719\u0026ndash;1725\u003c/li\u003e\n\u003cli\u003eWu Y, Bhat PR, Close TJ, Lonardi S (2008) Efficient and Accurate Construction of Genetic Linkage Maps from the Minimum Spanning Tree of a Graph. PLoS Genetics 4:e1000212\u003c/li\u003e\n\u003cli\u003eXu LS, Wang MN, Cheng P, et al (2013) Molecular mapping of Yr53, a new gene for stripe rust resistance in durum wheat accession PI 480148 and its transfer to common wheat. Theoretical and Applied Genetics 126:523\u0026ndash;533\u003c/li\u003e\n\u003cli\u003eYu L-X, Barbier H, Rouse MN, et al (2014) A consensus map for Ug99 stem rust resistance loci in wheat. Theoretical and Applied Genetics 127:1561\u0026ndash;1581\u003c/li\u003e\n\u003cli\u003eYu L-X, Morgounov A, Wanyera R, et al (2012) Identification of Ug99 stem rust resistance loci in winter wheat germplasm using genome-wide association analysis. Theoretical and Applied Genetics 125:749\u0026ndash;758\u003c/li\u003e\n\u003cli\u003eZhang W, Chen S, Abate Z, et al (2017) Identification and characterization of Sr13, a tetraploid wheat gene that confers resistance to the Ug99 stem rust race group. Proceedings of the National Academy of Sciences 114:E9483\u0026ndash;E9492\u003c/li\u003e\n\u003cli\u003eZhang W, Olson E, Saintenac C, et al (2010) Genetic Maps of Stem Rust Resistance Gene Sr35 in Diploid and Hexaploid Wheat. Crop Science 50:2464\u0026ndash;2474\u003c/li\u003e\n\u003cli\u003eZhu T, Wang L, Rimbert H, et al (2021) Optical maps refine the bread wheat Triticum aestivum cv. Chinese Spring genome assembly. The Plant Journal 107:303\u0026ndash;314\u003c/li\u003e\n\u003c/ol\u003e"}],"fulltextSource":"","fullText":"","funders":[],"hasAdminPriorityOnWorkflow":false,"hasManuscriptDocX":true,"hasOptedInToPreprint":true,"hasPassedJournalQc":"","hasAnyPriority":false,"hideJournal":false,"highlight":"","institution":"","isAcceptedByJournal":false,"isAuthorSuppliedPdf":false,"isDeskRejected":"","isHiddenFromSearch":false,"isInQc":false,"isInWorkflow":true,"isPdf":false,"isPdfUpToDate":true,"isWithdrawnOrRetracted":false,"journal":{"display":true,"email":"[email protected]","identity":"theoretical-and-applied-genetics","isNatureJournal":false,"hasQc":true,"allowDirectSubmit":false,"externalIdentity":"taag","sideBox":"Learn more about [Theoretical and Applied Genetics](https://www.springer.com/journal/122)","snPcode":"122","submissionUrl":"https://submission.nature.com/new-submission/122/3","title":"Theoretical and Applied Genetics","twitterHandle":"","acdcEnabled":true,"dfaEnabled":true,"editorialSystem":"em","reportingPortfolio":"Springer Hybrid","inReviewEnabled":true,"inReviewRevisionsEnabled":false},"keywords":"Khorasan wheat, Triticum turgidum ssp. turanicum, Stem rust, Ug99, Introgression ","lastPublishedDoi":"10.21203/rs.3.rs-2958205/v1","lastPublishedDoiUrl":"https://doi.org/10.21203/rs.3.rs-2958205/v1","license":{"name":"CC BY 4.0","url":"https://creativecommons.org/licenses/by/4.0/"},"manuscriptAbstract":"\u003cp\u003eThe \u003cem\u003ePuccinia graminis\u003c/em\u003e\u0026nbsp;f. sp.\u0026nbsp;\u003cem\u003etritici\u003c/em\u003e\u0026nbsp;(\u003cem\u003ePgt\u003c/em\u003e) Ug99 race group presents a major challenge to global wheat production. Satisfying current and future demands hinges on the discovery of new sources of resistance. It is imperative that the durability and diversity of Ug99 resistance is improved by identifying and deploying novel resistance genes. Progenitor species and wild relatives of common wheat (\u003cem\u003eTriticum aestivum\u003c/em\u003e) have proven to be rich sources of genetic diversity. The Khorasan wheat (\u003cem\u003eTriticum turgidum \u003c/em\u003essp.\u003cem\u003e turanicum\u003c/em\u003e) accession CItr 11390 displays adult plant resistance (APR) to Ug99 races. 121 BC\u003csub\u003e1\u003c/sub\u003eF\u003csub\u003e5\u003c/sub\u003e-derived recombinant inbred lines were developed from a cross between CItr 11390 and MN07098-6 to map and introgress resistance loci from CItr 11390. The population was evaluated in Kenya and Ethiopia in 2019 and 2020. Two APR QTL from CItr 11390 were detected in multiple environments. \u003cem\u003eQSr.umn-2A\u003c/em\u003e is believed to be the APR gene \u003cem\u003eSr63 \u003c/em\u003eon chromosome 2AL. \u003cem\u003eQSr.umn-6BL \u003c/em\u003ewas identified on 6BL upstream from \u003cem\u003eSr11\u003c/em\u003e. The distance from \u003cem\u003eSr11\u003c/em\u003e and lack of APR QTL reported on 6BL\u003cem\u003e \u003c/em\u003esuggest \u003cem\u003eQSr.umn-6BL\u003c/em\u003e is a novel locus. Additional QTL were mapped to chromosomes 1AS, 3AL, 3BL, 5AL, and 6BS in single environments. The population segregates for TKTTF seedling resistance conferred by \u003cem\u003eSr7a \u003c/em\u003eand a novel locus, \u003cem\u003eQSr.umn-5A.1\u003c/em\u003e. The population consists of the first hexaploid wheat lines to pyramid \u003cem\u003eSr7a\u003c/em\u003e, \u003cem\u003eSr57\u003c/em\u003e/\u003cem\u003eLr34\u003c/em\u003e/\u003cem\u003eYr18\u003c/em\u003e, \u003cem\u003eSr63\u003c/em\u003e, and \u003cem\u003eQSr.umn-6BL\u003c/em\u003e. This study is the first report of \u003cem\u003ePgt \u003c/em\u003eresistance QTL from Khorasan wheat, and it demonstrates the feasibility of simultaneously identifying and transferring resistance QTL from tetraploid to hexaploid wheat.\u003c/p\u003e","manuscriptTitle":"Direct hybridization facilitates the simultaneous identification and introgression of QTL for adult plant resistance to the Ug99 stem rust lineage from tetraploid Khorasan wheat to common wheat","msid":"","msnumber":"","nonDraftVersions":[{"code":1,"date":"2023-05-25 14:47:24","doi":"10.21203/rs.3.rs-2958205/v1","editorialEvents":[{"type":"communityComments","content":0},{"type":"reviewerAgreed","content":"","date":"2023-05-23T23:25:45+00:00","index":0,"fulltext":""},{"type":"reviewersInvited","content":"","date":"2023-05-23T23:15:34+00:00","index":"","fulltext":""},{"type":"editorAssigned","content":"","date":"2023-05-23T07:55:08+00:00","index":"","fulltext":""},{"type":"submitted","content":"Theoretical and Applied Genetics","date":"2023-05-22T14:44:49+00:00","index":"","fulltext":""}],"status":"published","journal":{"display":true,"email":"[email protected]","identity":"theoretical-and-applied-genetics","isNatureJournal":false,"hasQc":true,"allowDirectSubmit":false,"externalIdentity":"taag","sideBox":"Learn more about [Theoretical and Applied Genetics](https://www.springer.com/journal/122)","snPcode":"122","submissionUrl":"https://submission.nature.com/new-submission/122/3","title":"Theoretical and Applied Genetics","twitterHandle":"","acdcEnabled":true,"dfaEnabled":true,"editorialSystem":"em","reportingPortfolio":"Springer Hybrid","inReviewEnabled":true,"inReviewRevisionsEnabled":false}}],"origin":"","ownerIdentity":"a77b55ce-6543-4e64-a27a-b961470e50d4","owner":[],"postedDate":"May 25th, 2023","published":true,"recentEditorialEvents":[],"rejectedJournal":[],"revision":"","amendment":"","status":"in-revision","subjectAreas":[],"tags":[],"updatedAt":"2023-07-24T00:43:16+00:00","versionOfRecord":[],"versionCreatedAt":"2023-05-25 14:47:24","video":"","vorDoi":"","vorDoiUrl":"","workflowStages":[]},"version":"v1","identity":"rs-2958205","journalConfig":"researchsquare"},"__N_SSP":true},"page":"/article/[identity]/[[...version]]","query":{"redirect":"/article/rs-2958205","identity":"rs-2958205","version":["v1"]},"buildId":"J0_U0BvcaRcwD8yVFaRlm","isFallback":false,"isExperimentalCompile":false,"dynamicIds":[84888],"gssp":true,"scriptLoader":[]}