Mining the Vavilov wheat diversity panel for new sources of adult plant resistance to stripe rust

doi:10.21203/rs.3.rs-1001609/v1

Download PDF

Research Article

Mining the Vavilov wheat diversity panel for new sources of adult plant resistance to stripe rust

https://doi.org/10.21203/rs.3.rs-1001609/v1

This work is licensed under a CC BY 4.0 License

Journal Publication

published 03 Feb, 2022

Read the published version in Theoretical and Applied Genetics →

You are reading this latest preprint version

Wheat stripe rust (YR) caused by Puccinia striiformis f. sp. tritici (Pst), poses a significant threat to global food security. Resistance genes commonly found in many wheat varieties have been rendered ineffective due to the rapid evolution of the pathogen. To identify novel sources of adult plant resistance (APR), 292 accessions from the N.I. Vavilov Institute of Plant Genetic Resources, Saint Petersburg, Russia, were screened for known APR genes (i.e. Yr18, Yr29, Yr46, Yr33, Yr39 and Yr59) using linked polymerase chain reaction (PCR) molecular markers. Accessions were evaluated against Pst (pathotype 134 E16 A+ Yr17+Yr27) at seedling and adult plant stages across multiple years (2014, 2015 and 2016) in Australia. Phenotypic analyses identified 132 lines that potentially carry novel sources of APR to YR. Genome-wide association studies (GWAS) identified 68 significant marker-trait associations (P<0.001) for YR resistance, representing 47 independent quantitative trait loci (QTL) regions. Fourteen genomic regions overlapped with previously reported Yr genes, including Yr29, Yr56, Yr5, Yr43, Yr57, Yr30, Yr46, Yr47, Yr35, Yr36, Yrxy1, Yr59, Yr52 and YrYL. In total, seven QTL (positioned on chromosomes 1D, 2A, 3A, 3D, 5D, 7B and 7D) did not collocate with previously reported genes or QTLs, indicating the presence of promising novel resistance factors. Overall, the Vavilov diversity panel provides a rich source of new alleles which could be used to broaden the genetic bases of YR resistance in modern wheat varieties.

Agronomy

Molecular Genetics

Triticum aestivum

Puccinia striiformis

disease resistance

genome-wide association study

landraces

Multi-year evaluation of the Vavilov wheat diversity panel identified new sources of adult plant resistance to stripe rust. Genome-wide association studies revealed the key genomic regions influencing resistance, including seven novel loci.

Wheat is one of the most widely cultivated cereal crops worldwide, with an annual production of approximately 761 million tons (Crop Prospects and Food situation, 2020). Global wheat yields are threatened by climate change (Asseng et al., 2015, Barlow et al., 2015) and rapidly evolving pathogens, including rust diseases (Chaves et al., 2013). Among the rusts, stripe rust (YR) caused by Puccinia striiformis f.sp. tritici (Pst) is an economically important disease that has caused several major epidemics worldwide, resulting in significant production losses (Sanin and Nazarova 2010, Hovmoller et al., 2011, Ellis et al., 2014b, Xia et al., 2016a, Rahmatov, 2016, Ali et al., 2017). Historically, YR was prevalent in cooler climates, however, the majority of wheat growing areas in the world have now become prone to YR (Ali et al., 2014b, Hubbard et al., 2015). Currently, 88% of global wheat production is under threat to YR, which accounts for annual losses of more than one billion US dollars (Beddow et al., 2015a). Deployment of YR resistant varieties is the preferred method of rust disease management because it is cost-effective and reduces the reliance on fungicides (Chen, 2005a).

Genetic resistance to YR is broadly categorized into two major classes: all stage resistance (ASR) (e.g. R-genes) and adult plant resistance (APR-genes). ASR is often underpinned by a single gene with a large effect that provides effective resistance at all stages of plant growth. The R-gene interacts with the pathogen in a gene-for-gene relationship (Flor, 1971), and therefore is commonly referred to as “race-specific” resistance (Ellis et al., 2014a). When deployed in a variety grown on a large scale, strong selection pressure is exerted on the pathogen population to select for mutations that overcome the resistance mechanism. Therefore, this type of resistance is often rendered ineffective within just 3-5 years. In contrast, APR is typically controlled by multiple genes, each with minor or partial effect and are usually best expressed at adult growth stages. Most of the well-characterised APR genes are race nonspecific and usually confer a “slow rusting” phenotype, which is known to be more durable (Lagudah, 2011a, Ellis et al., 2014a, Mundt, 2014). Importantly, APR genes can contribute to high levels of resistance through additive or epistatic effects (Sorensen et al., 2014). Some of the genes, such as pleiotropic APRs and high temperature adult plant (HTAP) resistance genes, are highly valuable in breeding programs. For instance, APR genes Yr18, Yr29, Yr30 and Yr46 confer pleiotropic resistance to YR, leaf rust, stem rust and powdery mildew of wheat (Lan et al., 2015), while Yr18, Yr29, Yr36, Yr39 and Yr52 exhibit HTAP resistance (Chen, 2013b).

To date, 83 Yr resistance genes (Yr1–83) have been officially designated, along with 47 genes that have been temporarily named (McIntosh et al., 2019). According to the available information, most are classified as R-genes, whereas only 18 are classified as APR genes (Wu et al., 2016). Notably, only three APR genes Yr18/Lr34 (Krattinger et al., 2009), Yr36 (Fu et al., 2009), and Yr46 (Moore et al., 2015a) have been cloned to date (Liu et al., 2015). The majority of the R-genes, which have been deployed in various wheat varieties, are no longer effective due to the emergence of virulent pathotypes (Hovmoller et al., 2011). For instance, a large number of predominant races evolved from the year 2000 onwards and displayed added virulence to numerous resistance genes such as Yr2, Yr6, Yr7, Yr8, Yr9, Yr10, Yr17, Yr27, Yr43 and Yr44 (Wan and Chen, 2014). Virulence was also reported for some of the APR genes (Sorensen et al., 2014). Therefore, an additional level of durable genetic resistance could be achieved by pyramiding both seedling and APR genes in future varieties (Mundt, 2014). Hence, the discovery of new sources of genetic resistance is a priority for wheat research and the successful integration of new technologies in crop improvement programs is important to achieve long-term rust control.

Plant genetic resources that are stored in gene banks worldwide are a valuable source of genetic diversity for biotic and abiotic stresses (Rao, 2004b). Among these, wheat landraces and wild relatives are valuable sources of novel alleles for YR resistance (Sthapit et al., 2014, Manickavelu et al., 2016, Sthapit Kandel et al., 2017, Pasam et al., 2017). In previous studies, three important Yr genes (Yr47, Yr51 and Yr57) have been successfully characterized from wheat landraces in the Watkins collection (Bansal et al., 2011, Randhawa et al., 2015). Another historically important germplasm collection is preserved at the N.I. Vavilov Institute of Plant Genetic Resources (VIR) in Saint Petersburg, Russia, which holds ~38,430 wheat accessions. The collection comprises 76% bread wheat, 16% durum wheat and 7.9% wild and primitive wheats from diverse geographical origins (i.e., Africa, east and west Asia, USA, Canada, Central and South America, Europe) and of diverse biological status, including wild forms, local cultivars, breeding lines, mutants and artificial allopolyploids (Mitrofanova, 2012). While the VIR collection is yet to be explored for YR resistance, previous studies have highlighted the genetic variation for a number of biotic stresses (Mitrofanova, 2012, Sadovaya et al., 2015, Riaz et al., 2016a, Riaz et al., 2018).

To unravel the genetic architecture of rust resistance, bi-parental linkage mapping studies are traditionally performed (Yang et al,2017). However, high costs associated with population development, poor mapping resolution due to low recombination and the constraint of low allelic diversity, are some of the limitations associated with the linkage mapping approach (Flint-Garcia, 2013). Alternatively, a GWAS approach can be applied to a collection of accessions or a natural population. It offers broader allele coverage and higher mapping resolution due to historical recombination events among the panel of lines (Brachi et al., 2011a). This helps to localise the association signals to smaller regions within the chromosome and supports more efficient identification of candidate gene(s). However, applying GWAS to germplasm collections can be challenging because of population structure, which can result in spurious correlations between markers and traits (Gupta, Kulwal et al. 2019, Yang et al,2017). Another drawback is the low detection power of rare alleles with larger effects or multi-allelic variants with minor effects. Thus, an ideal GWAS analysis requires a large population size, high marker density and a mixed linear model to detect true genotype-phenotype associations (Bulli et al., 2016a). GWAS was initially successfully implemented for rust resistance in hexaploid wheat by Crossa et al. (2007b). Since then, it has been widely used in several studies to detect genomic regions associated with YR resistance in spring, winter, synthetic wheat germplasm, and landraces (Zegeye et al., 2014, Jighly et al., 2015, Maccaferri et al., 2015b, Naruoka et al., 2015,Bulli et al., 2016a, Godoy et al., 2017, Liu et al., 2017).

In this study, we evaluated 292 hexaploid bread wheat accessions from the VIR against Pst to identify novel sources of YR resistance. We applied GWAS to identify key genomic regions that could support the development of new cultivars incorporating durable resistance to YR.

Plant material

This study examined the Vavilov wheat diversity panel for resistance to YR. The composition of the panel was previously described by Riaz et al. (2017). Briefly, it includes 292 bread wheat accessions, including 136 landraces, 32 cultivars, 10 breeding lines and 115 lines with unknown cultivation classification, which were collected over a 70-year period (Online Resource 1). A total of 206 lines have origin information, whereas the origin is unknown for the remaining 89 lines. Notably, the diversity panel is both morphologically and genetically more diverse compared to Australian and CIMMYT elite wheat materials (Riaz et al., 2017).

Pathogen

A highly virulent and widespread Pst pathotype in Australia, 134 E16 A+ Yr17+Yr27, was used for screening the diversity panel. The pathotype 134 E16 A+ was first detected in Western Australia in 2002 where it initially displayed virulence for Yr6, Yr7, Yr8, Yr9, and YrA (Wellings et al., 2003) and gradually acquired virulence for Yr10 and Yr17 (Wellings, 2007). Virulence for Yr27 was first reported in 2008 and combined virulence for Yr17 and Yr27 in pathotype E16A+Yr17+Yr27+ was detected in 2010 (Randhawa et al., 2015).

PCR marker screening for known APR genes

A subset of 283 wheat accessions were screened for known APR genes. The accessions were screened by a PCR reaction, using markers corresponding to Yr18, Yr29, Yr46, Yr33, Yr39 and Yr59. Marker screening for Yr18, Yr29 and Yr46 (Riaz et al., 2016b) was performed with gene-specific cleaved amplified polymorphic sequence (CAPS) markers cssfr5 (Lagudah et al., 2009), csLV46 (unpublished) and a gene-specific single-nucleotide polymorphic marker SNP1-TM4 (Moore et al., 2015b), respectively. The genes Yr33 and Yr59 were evaluated, each with two closely linked markers to detect their presence or absence. SSR markers gwm111 and gwm437 were used to identify Yr33 (Zahravi et al., 2003), whereas the SSR marker barc32 and Resistance Gene Analog Polymorphism (RGAP) marker wmc557 were used to identify Yr59 (Zhou et al., 2014). The Yr39 gene was evaluated using the RGAP marker wgp45 (Lin and Chen, 2007). PCR was carried out in a 20 µl reaction containing upto 100 ng of genomic DNA, 1X GoTaq Flexi buffer, 1.5 mM MgCl₂, 200 µM dNTP, 200 nM of both forward and reverse primers and 1U of Taq polymerase (M829B, Promega, USA). A touch down PCR cycling condition was used as follows: denaturation at 94°C for 30 seconds; annealing at 62°C for 30 seconds, decreasing by 0.5°C per cycle; extension at 72°C for 1 minute followed by repeating these steps for 10 cycles; after enrichment, the program continued for 29 cycles as follows: 94°C for 30 seconds, 58°C for 30 seconds and 72°C for 40 seconds. Amplified products were evaluated using agarose gel electrophoresis.

Phenotyping for seedling resistance

All accessions were evaluated for Pst resistance at the seedling stage in a glasshouse (GH) at the University of Queensland, St Lucia, Queensland, Australia. The experiment was conducted using a complete randomised design with three replicates. Three to four seeds per accession were sown at four different positions around the circumference of 140 mm ANOVApots® filled with potting media. Seedlings were raised under glasshouse conditions maintained at day/night temperature of 22/17°C using a natural 12 h diurnal photoperiod. Eleven-day-old seedlings (i.e., two-leaf stage) were inoculated with Pst pathotype 134E16 A+ Yr17+Yr27 according to the air-brush protocol reported by Hickey et al. (2012). Seedlings were evaluated after 14 days post-inoculation based on 0–4 Stakman et al. (1962) scale which contains both numbers (e.g., 0, 1…4) and symbols (e.g., +). The numbers represent the disease score: highly resistant (HR, 0–1), moderately resistant (MR, 2), moderately susceptible (MS, 3), and susceptible (S, 4).

Phenotyping for adult plant resistance

A total of 292 accessions of the Vavilov wheat diversity panel were evaluated across three years (2014, 2015 and 2016) at The University of Queensland Gatton Research Station, Queensland, Australia. Wheat lines were sown in the field as un-replicated hill plots. The hill plots were 0.5 m apart along the direction of sowing and in two rows positioned between two spreader rows of the very susceptible genotype, Morocco (Figure 1a). To establish YR in the nursery, wheat seedlings that were raised in the glasshouse, infected with Pst, and were transplanted among the spreader rows approximately one month after sowing the nursery. The reaction to YR was evaluated using a modified Cobb’s scale (Peterson et al., 1948) and grouped on a 0‒9 scale in 2014, a 1‒9 scale (Bariana et al., 2007) was used in 2015 and 2016, the lower score values on both scales indicated increased resistance (Figure 1b). Disease assessments were performed at different time-points from heading (Zadoks 50) to grain filling stage (Zadoks 85). Specifically, data was collected at three time-points in 2014 and two time-points in 2015 and 2016. The seven phenotypic datasets were referred to as: Field_2014_1, Field_2014_2, Field_2014_3, Field_2015_1, Field_2015_2, Field_2016_1and Field_2016_2.

Genotyping

DNA of each accession was extracted according to the recommended DArT protocol (www.diversityarrays.com) and genotyped using the genotyping-by-sequencing (GBS) Diversity Array Technology (DArTseq) platform (Li et al., 2015). A total of 56,306 raw silicoDArT markers were returned, of which 34,311 were polymorphic in the panel. Only markers positioned on the current DArTseq consensus map (i.e., 14,228 polymorphic markers) were selected and filtered based on 10% and 15% threshold for missing data for markers and accessions, respectively. Markers with a minor allele frequency (MAF) <0.05 were also excluded. Additionally, three markers specific for known APR genes (Yr18, Yr29 and Yr46) were included. In total, 13,934 high quality markers and 292 accessions were retained and used for GWAS. These markers were ordered based on the genetic map positions in a high-resolution DArTseq consensus map (version 4.0), developed by Diversity Array Technology Inc., Canberra, ACT, Australia.

Data analysis

Prior to analysis, seedling resistance scores (Stakman 0-4) were converted to the 0-9 scale where, 0 = immune and 9 = very susceptible, using a conversion table (Ziems et al., 2014). The IT were converted as follows: 0;, ;n, ;, 1−, 1, 1+, 2−, 2, 2+, 2++, 3−, 3, 3+, 3++ and 4 were coded as 0, 0.5, 1, 2.5, 3, 3.5, 4, 5, 6, 6.5, 7, 8, 8.5 and 9 respectively. The accessions were grouped based on the phenotypic data (IT – Infection Type) of the seedlings; highly resistant (IT = 0‒3), intermediate (IT = 4‒6) and highly susceptible (IT = 7‒9). Accessions representing landraces, cultivars, breeding lines and the unclassified group displaying susceptible reactions at the seedling stage (IT = 7‒9) were selected (132 lines) for further analysis to identify lines carrying potentially novel APR genes.

Data analyses were performed using R (R Core Team, 2019). To investigate groups of accessions within the Vavilov diversity panel, hierarchical cluster analysis was performed using passport data (including biological status and origin) and the YR disease response data. Squared Euclidean distances were used to calculate the dissimilarity matrix and the tree was created according to Ward agglomerative clustering criterion. To compare the results across the seedling and adult plant experiments, correlation analyses were performed using mean disease severity (IT) scores. Principal component analysis (PCA) was performed using all datasets and results were displayed as a biplot to investigate the correlation of disease responses across experiments.

Genome-wide association studies

The genetic diversity and population structure of the Vavilov wheat diversity panel was previously described by Riaz et al. (2016). GWAS was performed using a compressed mixed linear model (MLM) (Yu et al., 2006) implemented in the R package known as Genome Association and Prediction Integrated Tool (GAPIT) (Lipka et al., 2012). To minimise spurious associations due to population structure and relatedness, the Q matrix was used as a fixed effect and the K matrix was used to fit a random genetic effect (Gupta et al., 2014). To decrease the Type II error rate, a relaxed significance threshold for marker-trait associations (MTA) of –log₁₀(P) > 3.0 was used, as previously employed in other GWAS studies examining the Vavilov panel (Riaz et al. 2018, Dinglasan et al. 2019). Markers detected in different assays, but in the same chromosomal position within a 1cM window were considered as the same QTL. To confirm the independence of closely positioned QTL, the local LD (linkage disequilibrium) value for associated markers was calculated in R using the ‘genetics’ package (Warnes et al., 2012). For each significantly associated marker, the resistance allele was assigned based on the direction of the allele effect on the resistance score value (Riaz et al., 2018). For example, if the “1” marker allele had a negative effect (e.g., -1.2), it was considered the resistance allele and the “0” marker allele was considered the susceptibility allele.

Comparison of QTL with previously reported Yr genes (QTL) and identification of putative candidate genes

To compare the position of the 47 QTL detected in this study with previously reported Yr genes and QTL, we projected the genetic positions on to the integrated consensus map developed by Maccaferri et al. (2015b) using the MapChart software version 3.2 (Voorrips, 2002). The integrated consensus map includes markers from the 9K (Cavanagh et al., 2013) and 90K (Wang et al., 2014b) consensus maps, the tetraploid consensus map (Maccaferri et al., 2015a), the Synthetic x Opata DH GBS map (Saintenac et al., 2013a), the Diversity Array Technology (DArT) integrated map (http://www.diversityarrays.com/search/node/wheat%20DArT%20map), SSR consensus map (Somers et al., 2004) and the Synthetic x Opata ITMI BARC SSR map (Song et al., 2005). In addition, recently reported genes were also projected onto the integrated map. All QTL regions were compared with the genes/QTLs reported by Chen (2017). QTL were considered novel when they were positioned more than 10 cM away from previously reported genes or QTL. The flanking sequences were BLAST searched against the Ensemble Plant genome-centric portal to find the physical position based on the wheat genome assembly IWGSC Ref-Sequence v1.1 (IWGSC 2018). For five DArTSeq markers (1255550, 981525, 1236960, 3026338,1159261), the 10+ genome assembly browsers at GrainGenes (http://www.graingenes.org) were used to retrieve the physical position. For each QTL, gene models and functional annotation was retrieved from the wheat genome assembly IWGSC Ref-Sequence v1.1(IWGSC 2018) using Ensemble Plant genome-centric portal.

Assigning and stacking of alleles associated with novel APR QTL

Seven newly discovered APR QTL were further analysed to determine their effect on disease response. To remove the masking effect of major resistance genes, only accessions displaying susceptibility at the seedling stage (i.e., IT scores >7) were selected for these analyses. The difference between the mean disease responses for groups of accessions carrying different combinations of resistance alleles was tested for significance based on Tukey’s test for multiple comparisons with a family-wise error rate of 5%. To investigate the trend between disease response and the accumulation of resistance alleles for the seven APR QTL, a relative disease index (YRi) was calculated across field trials according to:

$$YRi=\sum _{k}^{n}\frac{individual Dis. Score \left[k\right]}{mean Dis.Score \left[k\right]}$$

where YRi is the disease index of a line, in experiment k in relation to the population mean over n field experiments. A disease index below 0 reflects a high level of resistance to YR across the trials, whereas values above 0 imply susceptibility.

Analysis of stripe rust response at seedling and adult plant stages

In total, 292 accessions were tested for their IT at the seedling stage and the results revealed that 39% of the seedlings showed highly resistant reactions (IT = 0–3), 16% moderately resistant reactions (IT = 4–6) and 45% were susceptible (IT = 7–9). Most accessions displaying susceptibility were of unknown origin (n = 50; Online Resource 1), followed by accessions from Russia (n = 27), India (n = 17) and Pakistan (n = 11). In total, 132 lines displayed susceptible reactions to YR at the seedling stage and therefore these accessions were considered potential candidates for carrying uncharacterised APR genes.

Under field conditions, the accessions displayed a wide range of responses across the three years (Figure1a,b; 2). In 2014 and 2015, more than 80% of the tested lines (n = 284) displayed resistant to moderately resistant reactions, whereas in 2016, only 58% of the tested lines (n = 247) showed resistant reactions of which 43% displayed intermediate resistant phenotypes (Figure 2).

Weak positive correlations were observed between seedling and most field-based assessments; for example Field_2014_3 (r = 0.33, P < 0.001) and Field_2015_2 (r = 0.30, P < 0.001). The seedling disease response showed a moderate correlation with Field_2016_2 (r = 0.53, P < 0.001). Conversely, stronger correlations were observed between the field assessments; for example Field_2015_1 and Field_2016_2 (r = 0.81, P < 0.001). In Figure 3, the accessions displayed in black are seedling-susceptible and lack known genes, thus represent accessions that could carry uncharacterised APR.

Cluster analysis using phenotype and passport data

Ward’s hierarchical agglomerative cluster analysis was performed for 238 wheat accessions using trait data, including biological status (Unclassified, Classified – Breeding line, Cultivar, Landrace), origin of collection and YR disease response in the glasshouse and field experiments. The analysis revealed three major clusters divided into eight subgroups (Figure 4 and Online Resource 1). Cluster I comprised accessions belonging only to the unclassified group (i.e., unknown geographical information and unknown biological status). All the accessions in this cluster were grouped under the APR and moderate resistance (MR) category (Subgroups 2,3 and 5). The two other major clusters (Cluster II and III) comprised accessions belonging to a combination of classified and unclassified groups. Cluster II comprised accessions carrying APR (Subgroups 6 and 8) and accessions that were generally susceptible (Subgroup 4). Accessions in Sub-group 6 were more resistant and stable across the years compared to accessions belonging to Sub-group 8. Cluster III comprised accessions carrying mainly APR (Subgroup 7) and ASR (Subgroup 1).

Identification of novel sources of APR

Based on the results from YR assessments and PCR marker screening, 48 lines were deemed to carry potentially novel sources of APR to YR (Online resource 1). This included lines classified as unknown origin (n = 18), landraces (n = 14), cultivars (n = 10) and breeding lines (n = 6). The landraces were mainly from India (n = 3), Spain (n = 2), Russia (n = 2), and Kazakhstan (n = 2), although landraces from Pakistan, Turkey, Azerbaijan, Tajikistan and Armenia were also included in the set. Most cultivars (n = 8) and breeding lines (n = 5) that carried potentially novel sources of APR were from Russia. The 10 most promising lines were selected for crossing to elite cultivars. These showed a high degree of susceptibility at the seedling stage and intermediate levels of resistance at the adult plant stage (Figure 5). The majority of these lines were from the unclassified group (n = 6), followed by a landrace from Pakistan (WLA–043), a breeding line from Russia (WLA–302) and two cultivars from Chile (WLA–300) and Russia (WLA–087).

Marker‒trait associations

A total of 68 significant marker-trait associations were detected (-log₁₀ (p-value) > 3) which represented 47 different QTL regions (Table 1, Online Resource 2). Fifty-seven markers were associated with resistance at the adult plant stage, and many (n = 26) were detected only under the environmental conditions experienced in 2014. Nine markers were associated with resistance at the seedling stage while only one marker was detected at both seedling and adult plant stage. Based on the chromosomal position and LD between adjacent markers, the QTLs represented 41 APR loci and 6 seedling resistance loci. Comparison of these QTL regions with the previously reported Yr genes and QTL identified seven genomic regions as potentially novel QTL associated with YR resistance at the adult plant stage (Online Resource 2, 3). Notably, 14 QTL overlapped with the positions of previously reported Yr genes: Yr29, Yr5, Yr43, Yr57, Yr30, Yr51, Yr48, Yr47, Yr35, Yr36, Yrxy1, Yr59, Yr52 and YrYL on chromosome 1B, 2B, 3B, 4A, 5B, 6B, 7A, 7B and 7D, suggesting they could be linked to, or represent, those known resistance factors (Online resource 4).

Table 1

Summary of the stripe rust resistance QTL identified at both the seedling and adult plant stage in the Vavilov wheat diversity panel.
QTL name	Peak Marker ^a	Chr.	Position (cM)^b	Position (Mb)^c	Phenotypic dataset(s)^d	-log₁₀(P)	R. allele ^e	Effect on trait	Gene/QTL	References
qNV.Yr-1B-1	1219437	1B	241.91	660.6	Field_2_2016	3.14	1	-0.2	Previously reported	Maccaferri et al., 2015
qNV.Yr-1B-2	1119286	1B	269.26	677.4	Field_2_2016	3.34	0	-0.2	Yr29	Lan et al., 2014,
										Maccaferri et al., 2015
qNV.Yr-1D-1	1862252	1D	29.72	21	Field_3_2014	3.17	1	-0.8	New
qNV.Yr-1D-2	1229892	1D	47.7	35.7	Field_2_2016	3.23	1	-0.22	Previously reported	Maccaferri et al., 2015
qNV.Yr-2A.1	2295954	2A	9.93	12.2	Field_1_2015	3.2	1	0.05	Previously reported	Maccaferri et al., 2015
qNV.Yr-2A.2	1255550	2A	24.56	33.6	Seedling	3.5	1	0.05	Yr56	Bansal et al., 2014,
										Basnet et al., 2014
qNV.Yr-2A.3	1093102	2A	46.78	49.5	Field_2_2015	3.78	1	-0.2	Previously reported	Naruoka et al., 2015
qNV.Yr-2A.4	1122504	2A	60.52	87	Field_1_2015	3.59	1	0.08	Previously reported	Bulli et al., 2016
qNV.Yr-2A.5	1862406	2A	68.57	256.2	Field_2_2015	3.45	1	0.11	New
qNV.Yr-2A.6	1183380	2A	124.25	761.1	Field_2_2014	3.6	0	0.23	Previously reported	Dedryver et al., 2009
					Field_3_2014	3.76		0.79
qNV.Yr-2B.1	2289349	2B	78.31	583.7	Field_2_2016	3.45	1	0.29	Yr5	Godoy et al., 2017
qNV.Yr-2B.2	1235357		88.54	763	Field _2_2016	3.19	0	0.37	Yr43	Maccaferri et al., 2015,
	1233396		88.92	763.4	Seedling	3.43	0	-0.51		Xu et al., 2013
	3064531		88.98	763.5	Seedling	4.11	0	0.07
	2322232		89.52	765.6	Seedling	3.06	1	0.001
	981525		89.52	764.6	Seedling	3.23	1	-0.07
	1046333		89.52	759.7	Seedling	3.21	1	-0.06
	2322248		89.88	766.2	Field_1_2015	3.45	1	0.1
qNV.Yr-2B.3	1236960	2B	107.02	783.9	Field_2_2015	4.33	0	-0.33	Previously reported	Crossa et al., 2007,
	1117771			784.4	Field _2_2015	4.41	0	-0.4		Buerstamayr et al., 2014
					Field_2_2014	3.44		1.31
	1337290			790.6	Field_2_2014	4.11	1	-0.57
	1296890			790.6	Field_2_2015	3.01	0	-0.12
	1297955			790.9	Field_2_2015	3.63	0	-0.6
	4663965			794.9	Field_2_2015	4.15	0	-1.08
	1298829			800.1	Field_2_2015	3.24	0	-0.41
qNV.Yr-2D.1	2246033	2D	70.85	207.1	Field_2_2014	3.35	1	0.97	Previously reported	Mallard et al., 2005,
	2338436			337.4	Field_2_2015	3.08	1	0.15		Ren et al., 2012
qNV.Yr-2D.2	2243695	2D	136.45	619.6	Seedling	3.04	0	0.45	Previously reported	Suenaga et al., 2003,
										Melicher et al., 2008
qNV.Yr-3A.1	1123104	3A	44.75	170.1	Seedling	3.22	1	0.33	Previously reported	Zegeye et al., 2014
qNV.Yr-3A.2	2253468	3A	49.23	574.3	Field_1_2015	3.41	0	0.06	Previously reported	Pasam et al., 2017
qNV.Yr-3A.3	2295584	3A	53.26	600.3	Field_3_2014	3.09	1	1.5	Previously reported	Hou et al., 2015
qNV.Yr-3A.4	1116501	3A	109.54	100.4	Field_2_2015	3.23	0	-0.11	Previously reported	Rosewarne et al., 2013
qNV.Yr-3A.5	1237203	3A	142.13	738.9	Field_1_2014	4.38	0	0.6	New
qNV.Yr-3B	2307351	3B	4.97	5.3	Field_1_2016	4.28	0	-0.23	Yr57	Lan et al., 2014,
					Field_2_2016	3.03		-0.13	Yr30	Basnet et al., 2013
	1239212		5.25	5.3	Field_1_2016	3.83	0	-0.02
					Field_2_2016	3.24		-0.006
					Seedling	3.27		-0.07
qNV.Yr-3D.1	992483	3D	59.49	374	Field_1_2014	3.88	1	-0.41	New
qNV.Yr-3D.2	1228106	3D	123.84	582.5	Field_2_2014	3.45	0	3.73	Previously reported	Basnet et al., 2013
qNV.Yr-4A.1	1120960	4A	96.08	668.4	Field_2_2014	3.76	1	1.02	Previously reported	Bulli et al., 2016,
										Zegeye et al,2014
qNV.Yr-4A.2	1131408	4A	117.26	252.7	Field_2-2016	3.79	0	0.42	Previously reported	Vazques et al., 2012,
										Chen et al., 2012
qNV.Yr-4D.1	3021792	4D	43.44	158.4	Field_2_2016	3.97	1	-0.38	Yr46	Herrera-Foessel et al., 2011
					Field_2_2014	3.74		0.3
					Field_1_2016	3.4		-0.07
qNV.Yr-4D.2	1117919	4D	54.74	439.3	Field_1_2014	3.08	1	1.19	Previously reported	Ren et al., 2012
qNV.Yr-5A	1127003	5A	58.39	524.8	Seedling	4.8	0	-0.03	Previously reported	Manickavelu et al., 2016,
										Lan et al., 2014
qNV.Yr-5B.1	1088238	5B	12.41	15.3	Field_2_2015	3.26	0	NA	Yr47	Bansal et al., 2013
qNV.Yr-5B.2	1125706	5B	22.69	42.3	Field_2_2016	3.07	1	-0.01	Previously reported	Hao et al., 2011
qNV.Yr-5B.3	1104944	5B	27.83	54.2	Field_2_2015	3.07	0	0.04	Previously reported	Yang et al., 2013,
	1128792		29.32	53	Field_2_2015	3.03	0	0.01		Lu et al., 2014
qNV.Yr-5B.4	1107669	5B	46.54	590	Field_2_2014	3.15	1	0.44	Previously reported	Lu et al., 2009,
					Field_3_2014	3.41		1.64		Vazques et al., 2015
qNV.Yr-5D	991465	5D	58.61	424	Field_1_2014	3.62	0	0.14	New
qNV.Yr-6A	3022417	6A	100.76	616.3	Field_2_2016	3.61	0	-0.01	Previously reported	Vazquez et al., 2012,
										Bulli et al., 2016
qNV.Yr-6B.1	1302727	6B	2.41	9.8	Field_2_2014	3.64	0	-1.78	Yr35	Prins et al., 2011,
	1097180		2.41	0.3	Field_1_2015	3.19	0	-0.2		Bansal et al., 2013
	1121581		4.36	18	Field_1_2014	3.89	1	1.06
qNV.Yr-6B.2	1229248	6B	10.89	20.7	Field_3_2014	3.56	0	8.01	Previously reported	Prins et al., 2011, Bansal et al., 2013
qNV.Yr-6B.3	2290856	6B	23.19	117.1	Field_1_2015	3.63	1	0.12	Yr36	Pasam et al., 2017,
										Santra et al., 2008
qNV.Yr-6B.4	1058394	6B	35.34	576.5	Field_2_2016	3.55	0	0.09	Previously reported	Klarqyist et al.,.2016
qNV.Yr-6B.5	2293763	6B	79.73	710.6	Field_3_2014	3.07	1	-3.44	Previously reported	Basnet et al., 2013,
	1109194		82.95	715.8	Field_1_2015	3.19	0	-0.37		Rosewarne et al., 2013
					Field_2_2016	3.28		-0.31
qNV.Yr-7A.1	1244188	7A	96.07	648.3	Field_2_2015	3.34	0	0.15	Yrxy1	Zhou et al., 2011,
										Singh et al., 2014,
										Zegeye et al., 2014
qNV.Yr-7A.2	1090077	7A	149.56	712.3	Field_2_2015	3.85	1	-0.04	Previously reported	Rosewarne et al., 2013
										Vazques et al., 2012
qNV.Yr-7B.1	1126816	7B	20.97	170.5	Field_2_2014	3.26	1	2.31	New
qNV.Yr-7B.2	1059624	7B	38.66	138.7	Field_2_2014	3.23	1	0.12	Previously reported	Suenaga et al., 2003
	1259637		38.65	129.5	Field_1_2014	3.17		-0.01
qNV.Yr-7B.3	1298605	7B	58.63	607.6	Field_3_2014	3.21	0	-0.11	Previously reported	Pasam et al., 2017
	3026338		59.74	614.1	Field_3_2014	3.33	0	-2.98
	1159261		60.75	614.1	Field_1_2014	3.03	0	0.01
					Field_3_2014	3.24		0.6
qNV.Yr-7B.4	1066421	7B	118.47	718.1	Field_1_2014	3.36	1	0.14	Yr52	Ren et al., 2012,
									Yr59	Zhou et al., 2014
qNV.Yr-7D.1	1105401	7D	56.58	65.9	Field_2_2015	3.3	1	-0.27	YrYL	Wu et al., 2016
qNV.Yr-7D.2	2303363	7D	72.12	103.6	Field_1_2016	3.03	1	-0.04	New
^a Most strongly associated marker for the QTL
^b Chromosome position (cM) of the QTL on the DArTseq wheat consensus map (v 4.0)
^c Chromosome position (Mb) of the QTL on the physical map (Ref-Sequence v1.1 (IWGSC 2018), 10+ genome assembly browsers - GrainGenes)
^d The dataset for which the QTL was detected
^e Allele associated with resistance

Six QTL regions defined by 10 markers (2A, 2B, 2D, 3A, 3B and 5A) were associated with seedling resistance. Of these QTL, only four were detected at the seedling stage, while the remaining two QTL (qNV.Yr-2B.2 and qNV.Yr-3B) were detected at both seedling and adult plant growth stages. The QTL qNV.Yr-2B.2 included seven markers within the 88.5-89.9 cM region, which were associated with resistance across test environments (i.e., seedling, field trials in 2015 and 2016). It also co-located with the seedling resistance gene Yr43. The remaining three QTL were detected only in one adult plant environment (besides the seedling stage). The QTL with the strongest association was qNV.Yr-3B (-log₁₀ P ≥ 4) with larger marker effects on disease response (-0.23). This QTL co-located with partial resistance (and pleiotropic) gene Yr30 and the seedling resistance gene Yr57, which provides resistance against the Australian pathotype 134E16 A+Yr17+Yr27. However, due to the narrow genomic interval between these genes on chromosome 3B, we were unable to determine whether Yr30, Yr57 or both were present in this QTL.

A total of 41 QTL (57 markers) were considered “APR QTL” and were significantly associated with IT in one or multiple environments (Table 1). A total of 26, 23 and 16 markers were detected in 2014, 2015 and 2016, respectively, and clustered into different QTL regions. Comparison of the detected loci with the previously reported Yr genes and QTL (using the integrated map from Maccaferri et al., 2015b), revealed that 41 genomic regions likely represented previously reported genes and/or QTL (Table 1, Online Resource 2, 3 and 4). Interestingly, genomic regions harbouring the HTAP resistance genes Yr36, Yr59 and Yr52 were detected in the 2014 and 2015 field environments, whereas APR genes Yr29 and Yr30 were detected in 2016. Notably, high average maximum temperatures (28.5°C ‒ 31.4°C) were experienced during the YR assessment periods in 2014 and 2015 compared to 2016.

Seven QTL (chromosomes 1D, 2A, 3A, 3D, 5D, 7B and 7D) were detected at the adult plant stage and appeared to be novel (Table 1). The QTL qNV-Yr-7D.2 was detected under environmental conditions in 2016, whereas all other QTL in the D genome were detected in 2014. The QTL qNV.Yr-2A.5 and qNV.Yr-5D were significant in multiple field environments, while qNV.Yr-1D-1, qNV.Yr-3A.5, qNV.Yr-3D.1, qNV.Yr-7B.1, and QYr.uq-7D.2 were detected only in a single field environment.

Identification of candidate genes and their functional annotation

Candidate genes for the seven novel QTL had annotated protein functions that were related to disease resistance signaling pathways (Online Resource 5). For instance, qNV.Yr-1D-1 was positioned adjacent to candidate gene TraesCS1D01G042500, an Ankyrin repeat, which interact with nucleotide-binding, leucine-rich repeat (NLR) proteins in mediating plant immunity response (Wang et al., 2020). The qNV.Yr-2A5 was collocated with TraesCS2A02G229000, which encodes an RNA recognition motif that plays a key role in regulating defence-related signaling pathway (Woloshen et al., 2011). TraesCS3A02G524800 and TraesCS3D02G269300 were identified as gene models for qNV.Yr-3A.5 and qNV.Yr-3D.1, which encode Zinc finger protein and Spear family protein, respectively, and are known to paly roles in plant immunity and stress tolerance (Cui et al., 2018, Yang et al., 2018). The qNV.Yr-7B.1 colocated with TraesCS7B02G137200 which encodes a multi antimicrobial extrusion protein involved in regulating disease tolerance (Julião et al., 2020). Finally, TraesCS5D02G334700 and TraesCS7D02G155300 were considered candidate genes for qNV.Yr-5D and QYr.uq-7D.2, respectively, which are known to encode enzyme in the Cytochrome P450 superfamily which play a role in hormone signaling under stress conditions (Pandian et al., 2020).

Allele stacking effects and geographical distribution of the novel APR QTL

The seven novel APR QTL were investigated for their potential additive effect to reduce disease severity. A total of 132 accessions that were seedling susceptible were selected for the analysis. To further verify allelic effects in the subset, the mean resistance level of lines carrying alleles associated with resistance were compared with those carrying alleles for susceptibility at each locus. Significant allelic effects were verified through t-tests for qNV.Yr-2A.5, qNV.Yr-1D-1, qNV.Yr-5D and qNV.Yr-7B.1, whereas significant effects were not detected for qNV.Yr-3A.5, qNV.Yr-3D.1 and qNV.Yr-7B.1(Figure 6). To investigate the accumulation effect of the seven novel APR loci for YR resistance, lines were grouped based on the number of QTL present in each line.

A clear decreasing trend in the field relative disease index was observed as the number of resistance alleles increased from one to three (Figure 7a). None of the Vavilov lines combined more than three resistance alleles. Notably, most of the lines carrying three novel resistance loci (n = 23) showed a disease index lower than zero. The novel alleles occurred singly or in combination in accessions from India, Kazakhstan, Ukraine and Russia (Figure 7b). Further, countries from central Asia to west Asia such as Uzbekistan, Tajikistan, Armenia and Georgia also carried combinations of 2–3 novel alleles.

Vavilov wheat accessions offer diverse sources of stripe rust resistance

In this study, the Vavilov wheat accessions displayed a wide range of YR responses in both seedling and adult plant stages across multiple environments. At the seedling stage, about 55% of the accessions displayed high levels of resistance (39%, IT = 0-3) to intermediate levels of resistance (16%, IT = 4-6). Based on field assessments in 2014, 83% of the lines also showed resistance at the adult plant stage, indicating the presence of strong seedling resistance or combination of effective seedling and APR to YR. Alternatively, these resistant accessions could carry multiple APR genes, as the presence of just 2 or 3 minor genes can confer “near-immune” levels of resistance (Bariana et al., 2010, Singh et al., 2014b). Some genes are known to have a synergistic effect on disease resistance, such as Lr34/Yr18 (German and Kolmer, 1992), which directly contributes resistance to leaf rust, YR and powdery mildew (Rinaldo et al., 2017), while also contributing resistance to stem rust due to additive and epistatic interactions (Hiebert et al., 2016). These types of genes are more valuable when breeding for rust resistance because they are less likely to be selected and deployed alone, and some genes can enhance resistance levels to more than one economically important disease.

A total of 45% of the accessions (n = 132) displayed susceptibility at the seedling stage, but resistance at adult growth stages. Considering the diversity of the germplasm, it suggests a high frequency of multiple APR genes. Despite using the same pathotype each year, the correlation between disease response varied across the three years of field evaluations. This highlights the role of genotype × environment interactions influencing the expression of resistance. The low temperature and high rainfall in 2016 compared to 2014 and 2015 likely favoured disease development and could be the reason for increased susceptibility in 2016. Although, this does not rule out other factors related to seasonal variation, such as differences in physiological development across seasons (Wang et al., 2014) and abiotic stress (Guo et al., 2008). Such effects are well-described for the Lr34/Yr18 gene, where studies have found that gene expression is highly influenced by the genetic background, environmental factors and level of infection (Risk et al., 2012).

Screening for known APR genes (Yr18, Yr29, Yr46, Yr33, Yr39, Yr59) using linked molecular markers revealed some interesting relationships when accession passport information was explored. For instance, the presence of one or more gene combinations of Yr59 and Yr46 in almost all accessions from India and Pakistan, suggests that these APR genes are common in germplasm originating from these regions. A total of 65% of the 88 lines that displayed strong resistant responses in both seedling and adult plant stages were classified as landraces. Previous studies have also found that wheat landraces offer a rich source of disease resistance loci (Tala et al., 2011, Pasam et al. 2017).

By combining the field screening data with results from marker screening for known APR genes, we identified accessions that carried potentially novel resistance alleles. Interestingly, a number of lines originating from European and Asian countries such as Russia, Kazakhstan, Armenia, Ukraine, Tajikistan, Spain, Uzbekistan, Turkey, Mongolia and China showed an absence of known APR genes and clustered in the APR category, which suggests they likely carry novel resistance alleles. The top 10 lines deemed to carry potentially novel sources of APR were from Russia, Pakistan, Chile and the group of unknown origin. The Russian accessions WLA–302 (1978) and WLA–087 (1990) displayed intermediate resistance across the three years of evaluation. The landrace WLA–043 from Pakistan collected in 1936 was found to be a good source of APR to YR across the years. Pakistan is a diversity hotspot for wheat due to its wide agroecological zones, where many landraces that are genetically diverse with rare alleles have been found, particularly in the Himalayan region (Hirano et al., 2008). Evaluation of 5,700 accessions from USDA‒ARS National Small Grain Collection (NSGS) showed that Chile is another important geographic center of rust resistance with the highest number of accessions (28%) for stem rust resistance determined at the adult plant stage (Bonman et al., 2007). In this study, the cultivar WLA–300 collected from Chile in 1963 was identified as one of the best accessions for APR to YR. The top 10 accessions identified in this study should be proiritised for bi-parental mapping studies to characterise the underlying resistance loci.

GWAS reveals 47 genomic regions associated with resistance

This study identified 68 DArTseq markers, representing 47 QTL regions, that were significantly (P<0.001) associated with seedling and adult plant resistance to YR. Notably, a larger number of marker-trait associations (n = 26) were detected under field conditions in 2014, which represented 21 QTL regions. In contrast, the markers detected in 2016 (n = 15) represented 14 different genomic regions and 13 of them were unique to the field environment in 2016. Therefore, the field conditions experienced in 2016 resulted in the detection of many loci that were either not expressed or weakly expressed in the 2014 and 2015 field environments because GWAS was unable to detect them. This highlights the importance of testing material for rust resistance across multiple environments.

Among the 6 QTL detected at the seedling stage, a major haplotype block comprising seven markers was identified on chromosome 2B (qNV.Yr-2B.2). These markers were significantly associated with resistance at seedling, adult plant or both stages, suggesting the possible presence of more than one locus within this region. Previously, four different genes (Yr5, Yr44, Yr53 and Yr43) have all been mapped in this region within 35.1cM (Xu et al., 2013). Based on the positions of key markers, qNV.Yr-2B.2 is most likely Yr43. It should be noted that virulence for Yr43 and Yr44 occurs in many Pst races that emerged from 2000 onwards in the USA (Chen 2005, Wan and Chen 2014), while Yr5 and Yr53 have remained effective against all races. Therefore, it is important to study this QTL in detail to understand the specific gene/genes which confer resistance to the Australian pathotype 134E16 A+Yr17+Yr27. Such information would be very useful for breeding programs.

Three well-characterised APR multi-resistance loci were detected in this study, including Lr46/Yr29/Sr58/Pm39/Ltn2 on chromosome 1BS, Sr2/Yr30/Lr27/Pm on 3BL and Lr67/Yr46/Pm46/Sr55/Ltn3 on 4DS. These are the most commonly introduced and selected genes in global wheat breeding programs. Interestingly, effects associated with these genes were detected in the 2016 field environment that experienced relatively lower temperatures and higher rainfall, in comparison to the 2014 and 2015 seasons. Apart from these genes, only a few QTLs were detected in 2016 compared to the other years. The QTL (QYr.ucw-1B and QYr.ucw-1D) reported by Maccaferri et al. (2015a) and QTL (QYrst.orr-4A and QYrst.orr-6A) reported by Dolores Vazquez et al., (2012) were also detected in the present study in 2016. The resistance conferred by these QTL may be further studied in detail by analysing two seedling-susceptible landraces from Armenia (WLA–249) and Russia (WLA–315). Further, the HTAP resistance genes such as Yr36 (Uay et al, 2005), Yr52 (Ren etal,2012) and Yr59 (Zhou et al., 2014) coincided with QTL detected in 2014 and 2015 field environments which were much warmer than 2016, which further suggests these regions could be underpinned by HTAP resistance loci.

Other important QTL detected in this study were qNV.Yr-3B and qNV.Yr-7B.4 (Table 1). The first QTL, qNV.Yr-3B was detected at the seedling stage; however, this QTL was also detected at the adult stage under field conditions in 2016. This QTL co-located with seedling resistance gene Yr57 and adult plant pleiotropic gene Yr30. Yr57 is a broadly effective against Australian pre-and post-2000 Pst pathotypes 104 E137A+, 108E141A+, 110 E143A+, 134 E16A+, 134 E16A+Yr17+, 134 E16A+ Yr17+Yr27+ and 150 E16A+ (Randawa,2015), whereas Yr30 is closely linked to Sr2 and confers partial resistance to YR. These two genes are closely linked (< 5 cM apart), which could be exploited for gene pyramiding (Randawa, 2015). The second QTL, qNV.Yr-7B.4 was represented by DArT marker 1066421 in the field in 2014. Notably, this region on chromosome 7B harbours two HTAP genes, Yr52 (Ren et al, 2012) and Yr59 (Zhou et al, 2014), and two temporally designated genes YrC591 (Li et al, 2008) and YrZH84 (Li et al, 2006). Yr52 confers HTAP resistance, which is highly consistent and partial in nature, making it an attractive target for breeding programs targeting durable rust resistance (Ren et al, 2012). Yr59 is another important HTAP resistance gene, which is linked to YrZH84 (Zhou et al., 2014) and is effective against all known North American Pst races, thus it is considered highly valuable for resistance breeding.

Novel QTL for adult plant resistance

Seven new QTL for YR resistance were identified on chromosomes 1D, 2A, 3A, 3D, 5D, 7B and 7D. These genomic regions were associated with different levels of resistance across the three years of field assessments. qNV.Yr-2A.5 (chromosome 2A) and qNV.Yr-5D (chromosome 5D) were associated with resistance in 2014 and 2015, while other QTL were detected only in a single environment. Therefore, these QTL may be environment-specific or convey smaller effects. Notably, the majority of novel QTL detected in this study were located in the D-genome. Aegilops tauschii, the D-genome progenitor of hexaploid wheat, is considered a valuable source of YR resistance (Liu et al, 2013). However, until now, only a few resistance loci have been identified in the D-genome. For instance, Yr28 (Singh et al, 2000), slow rusting APR gene Yr46 (Herrera foessel, 2011) and YrAs2388R (Zhang et al, 2019) have been mapped on chromosome 4D. Therefore, the novel QTL detected in the current study provides valuable sources for APR. Accession numbers WLA-106 (Ukrain), WLA-124 (Russia) for qNV.Yr-1D-1, WLA-43 (Pakistan) for qNV.Yr-3D.1 and qNV.Yr-5D, and WLA-114 (Azerbaijan) for QYr.uq-7D.2 represent valuable donor lines that can be used for pre-breeding or future research seeking to fine-map or clone the underlying gene(s).

Analysis of the annotated genes with the novel QTL regions revealed key candidates with roles related to plant defence (e.g. ankyrin repeat, RNA recognition motif, Zinc finger protein, Spear family protein,multi antimicrobial extrusion protein and Cytochrome P450 superfamily). While these candidate genes occur in the reference genome, their presence should be confirmed in the resistant accessions prior to transgenic efforts to validate gene function.

Novel APR QTL confer additive effects to reduce disease

We demonstrated that pyramiding of novel alleles associated with APR can significantly reduce YR disease. As a group, APR lines that carry a combination of 2‒3 resistance alleles show stronger levels of resistance compared to lines carrying just one resistance allele. However, allele effects are variable depending on the QTL, the QTL combination, genetic background and environmental conditions. For example, WLA-043 from Pakistan carries three QTL (3A, 3D and 5D) and showed stable resistance across years in comparison with WLA-279 from Kazakhstan which also carries three QTL (1D, 2A and 3A). It has been well described that some resistance genes (such as Yr18, Yr30 and Yr39) can increase the effectiveness of other R genes or weak APR genes (Chen et al, 2013, Ellis et al, 2014). While Yr18 and Yr49 confer additive effects for YR resistance, increased stem rust severity was detected in Chinese Spring due to the inactivation of Lr34 gene in the presence of Sr2 (Ellis et al, 2014). Therefore, it is important to understand gene interactions to incorporate effective disease resistance into wheat varieties.

This is the first study to document YR resistance in diverse wheat accessions sourced from the VIR. GWAS identified 47 genomic regions including 68 markers that were significantly associated with resistance. Over 90% of the QTL detected in this study aligned with previously reported genes and QTL, which highlights the power of association mapping and accuracy of resistance loci that have been identified in the present study. The molecular markers linked to QTL identified in this study will be valuable for further validation of resistance loci and potentially future marker-assisted selection breeding approaches. Further, insight into the distribution of resistance alleles and their frequencies in relation to geographical origin could aid selection of specific sources of resistance in breeding programs. The significant effect of pyramiding novel APR alleles identified in this study underscores the potential for breeding wheat varieties with improved YR resistance.

Author’s contribution:DJ, AR and LTH conceived and designed the study; DJ, AR, NA and WLN collected the data; DJ wrote the manuscript; DJ, EGD, KVF, SA, LZ, OA, EA, SKP and IG contributed to data analyses and interpretation of results, GP guided the design and establishment of rust screening nurseries in Queensland, Australia; all authors read and revised the manuscript.

Acknowledgements

DJ was supported by an Australian Postgraduate Award. LTH was supported by a Discovery Early Career Researcher Award (DECRA) Fellowship from the Australian Research Council (ARC) and an Early Career Researcher grant from UQ. We give thanks to the past and present members of the N.I. Vavilov Institute of Plant Genetic Resources in St Petersburg, Russia, for preserving and sharing the historic wheat germplasm. We are grateful to the technical staff at UQ Central Glasshouse Facilities at St. Lucia and the UQ Gatton Research Farm for assistance with establishing and maintaining experiments.

Ethical standard

This article does not contain any studies with human participants or animals performed by any of the authors.

Conflict of interest

The authors declare that they have no conflict of interest.

Data availability

The datasets generated during and/or analysed during the current study (i.e. phenotypes, genotypes, and mapping data) are available from the corresponding author on reasonable request.

Ali S, Gladieux P, Rahman H, Saqib MS, Fiaz M, Ahmad H, Leconte M, Gautier A, Justesen AF, Hovmoller MS, Enjalbert J, De Vallavieille-Pope C (2014) Inferring the contribution of sexual reproduction, migration and off-season survival to the temporal maintenance of microbial populations: a case study on the wheat fungal pathogen Puccinia striiformis f. sp. tritici. Mol Ecol 23:603–617. doi.org/10.1111/mec.12629
Ali S, Rodriguez-Algaba J, Thach T, Sorensen CK, HansenJG, Lassen P, NazariK, HodsonDP JustesenAF, Hovmoller MS (2017) Yellow rust epidemics worldwide were caused by pathogen races from divergent genetic lineages. Front Plant Sci. doi.org/10.3389/fpls.2017.01057
IWGSC (2018) IWGSC whole-genome assembly principal investigators; Pozniak CJ et al. Shifting the limits in wheat research and breeding using a fully annotated reference genome. Science, 361 DOI:10.1126/science.aar7191
Asseng S, Ewert F, Martre P, Rotter R, Lobel D, Cammarano D, Kimball B, Ottman MJ, Wall G, White J et al (2015) Rising temperatures reduce global wheat production. Nat Clim Chang 5:143–147. doi:10.1038/nclimate2470
Bansal UK, Arief VN, Delacy IH, Bariana HS (2013) Exploring wheat landraces for rust resistance using a single marker scan. Euphytica 194:219–233. doi:10.1007/s10681-013-0940-0
BansalUK, Forrest KL, Hayden MJ, Miah H, Singh D, Bariana HS (2011) Characterisation of a new stripe rust resistance gene Yr47 and its genetic association with the leaf rust resistance gene Lr52. Theor Appl Genet 122:1461–4666. doi:10.1007/s00122-011-1545-4
Bariana HS, Bansal UK, Schmidt A, Lehmensiek A, Kaur J, Miah H, Howes N, Mcintyre CL (2010) Molecular mapping of adult plant stripe rust resistance in wheat and identification of pyramided QTL genotypes. Euphytica 176:251–260. doi.org/10.1007/s10681-010-0240-x
Bariana HS, Miah H, Brown GN, Willey N, Lehmensiek A (2007) Molecular mapping of durable rust resistance in wheat and its implication in breeding. Dordrecht. Springer Netherlands, 723–728
Barlow KM, Christy BP, O'Leary GJ, Riffkin PA, Nuttall JG (2015) Simulating the impact of extreme heat and frost events on wheat crop production: A review. Field Crops Res 171:109–119. doi.org/10.1016/j.fcr.2014.11.010
Basnet BR, Ibrahim AM, Chen X, Singh RP, Mason ER, Bowden RL, Liu S, Hays DB, Devkota RN, Subramanian NK (2014b) Molecular mapping of stripe rust resistance in hard red winter wheat TAM 111 adapted to the US High Plains. Crop Sci 54:1361–1373. doi.org/10.2135/cropsci2013.09.0625
Basnet B, Singh R, Ibrahim A, Herrera-Foessel S, Huerta-Espino J, Lan C, Rudd J (2014a) Characterization of Yr54 and other genes associated with adult plant resistance to yellow rust and leaf rust in common wheat Quaiu 3. Mol breed 33:385–399. doi:10.1007/s11032-013-9957-2
Beddow JM, Pardey PG, Chai Y, Hurley TM, Kriticos DJ, Braun HJ, Park RF, Cuddy WS, Yonow T (2015) Research investment implications of shifts in the global geography of wheat stripe rust. Nat Plants 1:15132. doi:10.1038/nplants.2015.132
Bonman JM, Bockelman HE, Jin Y, Hijmans RJ, Gironella AIN (2007) Geographic distribution of stem rust resistance in wheat landraces. Crop Sci 47:1955–1963. doi.org/10.2135/cropsci2007.01.0028
Brachi B, Morris GP, Borevitz JO (2011) Genome-wide association studies in plants: the missing heritability is in the field. Genome Bio 12:232. doi.org/10.1186/gb-2011-12-10-232
Bulli P, Zhang J, Chao S, Chen X, Pumphrey M (2016a) Genetic architecture of resistance to stripe rust in a global winter wheat germplasm collection. G3 9:6(8): 2237–2253. doi: 10.1534/g3.116.028407
Chaves MS, Martinelli JA, Wesp-Guterres C, Graichen FAS, Brammer SP, Scagliusi SM, Da Silva PR, Wietholter P, Torres GAM, Lau EY (2013) The importance for food security of maintaining rust resistance in wheat. Food secur 5:157–176. doi 10.1007/s12571-013-0248-x
Chen X (2005) Epidemiology and control of stripe rust [Puccinia striiformis f. sp. tritici] on wheat. Can J Plant Patho 27:314–337. doi.org/10.1080/07060660509507230
Chen X (2013) Review article: High-temperature adult-plant resistance, key for sustainable control of stripe rust. Am J Plant Sci 4:605–627. doi:10.4236/ajps.2013.43080
Chen X (2017) Stripe Rust Epidemiology. In: Chen, Kan Z (eds) Stripe Rust. Springer Netherlands, Dordrecht
Crossa J, Burgueno J, Dreisigacker S, Vargas M, Herrera-Foessel SA, Lillemo M, Singh RP, Trethowan R, Warburton M, Franco J (2007) Association analysis of historical bread wheat germplasm using additive genetic covariance of relatives and population structure. Genetics 177:1889–1913. doi.org/10.1534/genetics.107.078659
Cui B, Pan Q, Clarke D et al (2018) S-nitrosylation of the zinc finger protein SRG1 regulates plant immunity. Nat Commun 9:4226. doi.org/10.1038/s41467-018-06578-3
Dinglasan EG, Singh D, Shankar M, Afanasenko O, Platz G, Godwin ID, Voss-Fels KP, Hickey LT (2019) Discovering new alleles for yellow spot resistance in the Vavilov wheat collection. Theor Appl Genet 132:149–162. doi.org/10.1007/s00122-018-3204-5
Dedryver F, Paillard S, Mallard S, Robert O, Trottet M, Negre S, Verplancke G, Jahier J (2009) Characterization of genetic components involved in durable resistance to stripe rust in the bread wheat 'Renan'. Phytopathology 99:968–973. doi:10.1094/PHYTO-99-8-0968
Dolores Vazquez M, James Peterson C, Riera-Lizarazu O, Chen X, Heesacker A, Ammar K, Crossa J, Mundt C (2012) Genetic analysis of adult plant, quantitative resistance to stripe rust in wheat cultivar 'Stephens' in multi-environment trials. Theor Appl Genet 124:1–11. doi:10.1007/s00122-011-1681-x
Ellis JG, Lagudah ES, Spielmeyer W, Dodds PN (2014) The past, present and future of breeding rust resistant wheat. Front Plant Sci 5:641. doi.org/10.3389/fpls.2014.00641
Crop Prospects and Food Situation (2020) http://www.fao.org/3/ca9803en/CA9803EN.pdf. doi.org/10.4060/ca9803en
Flint-Garcia SA (2013) Genetics and consequences of crop domestication. J Agric Food Chem 61:8267–8276. doi:10.1021/jf305511d
Flor HH (1971) Current status of the gene-for-gene concept. Ann Rev of Phytopathol 9:275–296. doi.org/10.1146/annurev.py.09.090171.001423
Fu D, Uauy C, Distelfeld A, Blechl A, Epstein L, Chen X, Sela H, Fahima T, Dubcovsky J (2009) A kinase-START gene confers temperature-dependent resistance to wheat stripe rust. Science 323:1357–1360. doi:10.1126/science.1166289
German SE, Kolmer JA (1992) Effect of gene Lr34 in the enhancement of resistance to leaf rust of wheat. Theor Appl Genet 84:97–105. doi:10.1007/BF00223987
Godoy JG, Rynearson S, Chen X, Pumphrey M (2017) Genome-wide association mapping of loci for resistance to stripe rust in North American elite spring wheat germplasm. Phytopathology 108:234–245. doi.org/10.1094/PHYTO-06-17-0195-R
Guo Q, Zhang ZJ, Xu YB, Li GH, Feng J, Zhou Y (2008) Quantitative trait loci for high-temperature adult-plant and slow-rusting resistance to Puccinia striiformis f. sp. tritici in wheat cultivars. Phytopathology 98:803–809. doi.org/10.1094/PHYTO-98-7-0803
Gupta PK, Kulwal PL, Jaiswal V (2014) Association mapping in crop plants: opportunities and challenges. Adv Genet 85:109–147. doi:10.1016/B978-0-12-800271-1.00002-0
Gupta PK, Kulwal PL, Jaiswal V (2019) Association mapping in plants in the post-GWAS genomics era. Adv Genet 104:75–154. doi:10.1016/bs.adgen.2018.12.001
Hao Y, Chen Z, Wang Y, Bland D, Buck J, Brown-Guedira G, Johnson J (2011) Characterization of a major QTL for adult plant resistance to stripe rust in US soft red winter wheat. Theor Appl Genet 123:1401–1411. doi:10.1007/s00122-011-1675-8
Herrera-Foessel SA, Lagudah ES, Huerta-Espino J, Hayden MJ, Bariana HS, Singh D, Singh RP (2011) New slow-rusting leaf rust and stripe rust resistance genes Lr67 and Yr46 in wheat are pleiotropic or closely linked. Theor Appl Genet 122:239–249. doi:10.1007/s00122-010-1439-x
Hickey LT, Wilkinson PM, Knight CR, Godwin ID, Kravchuk OY, Aitken EAB, Bansal UK, Bariana HS, Delacy IH, Dieters MJ (2012) Rapid phenotyping for adult-plant resistance to stripe rust in wheat. Plant Breed 131:54–61. doi.org/10.1111/j.1439-0523.2011.01925.x
Hiebert CW, Kolmer JA, Mccartney CA, Briggs J, Fetch T, Bariana H, Choulet F, Rouse MN, Spielmeyer W (2016) Major gene for field stem rust resistance co-locates with resistance gene Sr12 in 'Thatcher' wheat PLOS ONE, 11. doi.org/10.1371/journal.pone.0157029
Hirano R, Kikuchi A, Kawase M, Watanabe KN (2008) Evaluation of genetic diversity of bread wheat landraces from Pakistan by AFLP and implications for a future collection strategy. Genet Resour Crop Evol 55:1007–1015. doi:10.1007/s10722-008-9308-z
Hou L, Chen X, Wang M, See DR, Chao S, Bulli P, Jing J (2015) Mapping a large number of qtl for durable resistance to stripe rust in winter wheat druchamp using SSR and SNP markers. PLOS ONE, 10, doi.org/10.1371/journal.pone.0126794
Hovmoller MS, Sorensen CK, Walter S, Justesen AF (2011) Diversity of Puccinia striiformis on cereals and grasses. Annu Rev Phytopathol 49:197–217. doi:10.1146/annurev-phyto-072910-095230
Hubbard A, Lewis CM, Yoshida K, Ramirez-Gonzalez RH, De Vallavieille-Pope C, Thomas J, Kamoun S, Bayles R, Uauy C, Saunders DG (2015) Field pathogenomics reveals the emergence of a diverse wheat yellow rust population. Genome Biol 16:23. doi.org/10.1186/s13059-015-0590-8
Jighly A, Oyiga BC, Makdis F, Nazari K, Youssef O, Tadesse W, Abdalla O, Ogbonnaya FC (2015) Genome-wide DArT and SNP scan for QTL associated with resistance to stripe rust (Puccinia striiformis f. sp. tritici) in elite ICARDA wheat (Triticum aestivum L.) germplasm. Theor Appl Genet 128:1277–1295. doi:10.1007/s00122-015-2504-2
Julião MHM, Silva SR, Ferro JA, Varani AM (2020) A genomic and transcriptomic overview of MATE, ABC, and MFS transporters in Citrus sinensis interaction with Xanthomonas citri subsp. citri. Plants 9:794. doi:10.3390/plants9060794
Klarquist FE, Chen MX, Carter HA (2016) Novel QTL for stripe rust resistance on chromosomes 4A and 6B in soft white winter wheat cultivars. Agronomy 6:4. doi.org/10.3390/agronomy6010004
Krattinger SG, Lagudah ES, Spielmeyer W, Singh RP, Huerta-Espino J, Mcfadden H, BossolinI E, Selter LL, Keller B (2009) A putative ABC transporter confers durable resistance to multiple fungal pathogens in wheat. Science 323:1360–13633. doi:10.1126/science.1166453
Lagudah ES (2011) Molecular genetics of race non-specific rust resistance in wheat. Euphytica 179:81-91.doi.org/10.1007/s10681-010-0336-3
Lagudah ES, Krattinger SG, Herrera-Foessel S, Singh RP, Huerta-Espino J, Spielmeyer W, Brown-Guedira G, Selter LL, Keller B (2009) Gene-specific markers for the wheat gene Lr34/Yr18/Pm38 which confers resistance to multiple fungal pathogens. Theor Appl Genet 119:889–898. doi.org/10.1007/s00122-009-1097-z
Lan C, Rosewarne GM, Singh RP, Herrera-Foessel SA, Huerta-Espino J, Basnet BR, Zhang Y, Yang E (2014) QTL characterization of resistance to leaf rust and stripe rust in the spring wheat line Francolin#1. Mol Breed 34:789–803. doi:10.1007/s11032-014-0075-6
Lan C, Zhang Y, Herrera-Foessel SA, Basnet BR, Huerta-Espino J, Lagudah ES, Singh RP (2015) Identification and characterization of pleiotropic and co-located resistance loci to leaf rust and stripe rust in bread wheat cultivar Sujata. Theor Appl Genet 128:549–561. doi:10.1007/s00122-015-2454-8
Li Y, Niu YC, Chen XM (2008) Mapping a stripe rust resistance gene YrC591 in wheat variety C591 with SSR and AFLP markers. Theor Applied Genet 118:339–346. doi:10.1007/s00122-008-0903-3
Li ZF, Zheng TC, He ZH, Li GQ, Xu SC, Li XP, Yang GY, Singh RP, Xia XC (2006) Molecular tagging of stripe rust resistance gene YrZH84 in Chinese wheat line Zhou 8425B. Theor Appl Genet 112:1098–1103. doi:10.1007/s00122-006-0211-8
Li H, Vikram P, Singh RP, Kilian A, Carling J, Song J, Burgueno-Ferreira JA, Bhavani S, Huerta-Espino J, Payne T, Sehgal D, Wenzl P, Singh S (2015) A high density GBS map of bread wheat and its application for dissecting complex disease resistance traits. BMC genom 16:216. doi.org/10.1186/s12864-015-1424-5
Lin F, Chen XM (2007) Genetics and molecular mapping of genes for race-specific all-stage resistance and non-race-specific high-temperature adult-plant resistance to stripe rust in spring wheat cultivar Alpowa. Theor Appl Genet 114:1277–1287. doi:10.1007/s00122-007-0518-0
Lipka AE, Tian F, Wang Q, Peiffer J, Li M, Bradbury PJ, Gore MA, Buckler ES, Zhang Z (2012) GAPIT: genome association and prediction integrated tool. Bioinform 28:2397–2399. doi:10.1093/bioinformatics/bts444
Liu M, Zhang C, Yuan C, Zhang L, Huang L, Wu J, Wang J, Zheng Y, Zhang H, Liu D, Fu D (2013) Stripe rust resistance in Aegilops tauschii germplasm. Crop Sci 53:2014. doi.org/10.2135/cropsci2013.01.0008
Liu J, He Z, Wu L, Bai B, Wen W, Xie C, Xia X (2015) Genome-wide linkage mapping of qtl for adult-plant resistance to stripe rust in a Chinese wheat population Linmai 2 x Zhong 892. PLoS One, 10, doi.org/10.1371/journal.pone.0145462
Liu W, Maccaferri M, Bulli P, Rynearson S, Tuberosa R, Chen X, Pumphrey M (2017) Genome-wide association mapping for seedling and field resistance to Puccinia striiformis f. sp. tritici in elite durum wheat. Theor Appl Genet 130:649–667. doi:10.1007/s00122-016-2841-9
Lu Y, Lan C, Liang S, Zhou X, Liu D, Zhou G, Lu Q, Jing J, Wang M, Xia X, He Z (2009) QTL mapping for adult-plant resistance to stripe rust in Italian common wheat cultivars Libellula and Strampelli. Theor Appl Genet 119::1349–1359. doi:10.1007/s00122-009-1139-6
Lu Y, Wang M, Chen X, See D, Chao S, Jing J (2014) Mapping of Yr62 and a small-effect QTL for high-temperature adult-plant resistance to stripe rust in spring wheat PI 192252. Theor Appl Genet 127:1449–1459. doi.org/10.1007/s00122-014-2312-0
Maccaferri M, Ricci A, Salvi S, Milner SG, Noli E, Martelli PL, Casadio R, Akhunov E, Scalabrin S, Vendramin V, Ammar K, Blanco A, Desiderio F, Distelfeld A, Dubcovsky J, Fahima T, Faris J, Korol A, Massi A, Mastrangelo AM, Morgante M, Pozniak C, N'Diaye A, Xu S, Tuberosa R (2015a) A high-density, SNP-based consensus map of tetraploid wheat as a bridge to integrate durum and bread wheat genomics and breeding. Plant Biotechnol J 13:648–663. doi:10.1111/pbi.12288
Maccaferri M, Zhang J, Bulli P, Abate Z, Chao S, Cantu D, Bossolini E, Chen X, Pumphrey M, Dubcovsky J (2015b) A genome-wide association study of resistance to stripe rust (Puccinia striiformis f. sp. tritici) in a worldwide collection of hexaploid spring wheat (Triticum aestivum L.). G3 (Bethesda), 5:449-65.doi: 10.1534/g3.114.014563
McIntosh RA, Dubcovsky J, Roger WJ, Xia XC, Raupp WJ (2019) Catalogue of Gene Symbols for Wheat. http://www.shigen.nig.ac.jp/wheat/komugi/top/top.jsp
Mallard S, Gaudet D, Aldeia A, Abelard C, Besnard AL, Sourdille P, Dedryver F (2005) Genetic analysis of durable resistance to yellow rust in bread wheat. Theor Appl Genet 110:1401–1409. doi:10.1007/s00122-005-1954-3
Manickavelu A, Joukhadar R, Jighly A, Lan C, Huerta-Espino J, Stanikzai AS, Kilian A, Singh RP, Ban T (2016) Genome wide association mapping of stripe rust resistance in Afghan wheat landraces. Plant Sci 252:222–229. doi:10.1016/j.plantsci.2016.07.018
Melichar JP, Berry S, Newell C, Maccormack R, Boyd LA (2008) QTL identification and microphenotype characterisation of the developmentally regulated yellow rust resistance in the UK wheat cultivar Guardian. Theor Appl Genet 117:391–399. doi:10.1007/s00122-008-0783-6
Mitrofanova O (2012) Wheat genetic resources in Russia: Current status and prebreeding studies. Russ J Genet Appl Res 2:277–285. doi:10.1134/S2079059712040077
Moore JW, Herrera-Foessel S, Lan C, Schnippenkoetter W, Ayliffe M, Huerta-Espino J, Lillemo M, Viccars L, Milne R, Periyannan S, Kong X, Spielmeyer W, Talbot M, Bariana H, Patrick JW, Dodds P, Singh R, Lagudah E (2015) A recently evolved hexose transporter variant confers resistance to multiple pathogens in wheat. Nat Genet 47:1494–1498. doi:10.1038/ng.3439
Mundt CC (2014) Durable resistance: a key to sustainable management of pathogens and pests. Infect Genet Evol 27:446–455. doi:10.1016/j.meegid.2014.01.011
Naruoka Y, Garland-Campbell KA, Carter AH (2015) Genome-wide association mapping for stripe rust (Puccinia striiformis F. sp. tritici) in US Pacific Northwest winter wheat (Triticum aestivum L.). Theor Appl Genet 128:1083–1101. doi:10.1007/s00122-015-2492-2
Pandian BA, Sathishraj R, Djanaguiraman M, Prasad PVV, Jugulam M (2020) Role of cytochrome P450 enzymes in plant stress response. Antioxidants 9:454. doi:10.3390/antiox9050454
Pasam RK, Bansal U, Daetwyler HD, Forrest KL, Wong D, Petkowski J, Willey N, Randhawa M, Chhetri M, Miah H, Tibbits J, Bariana H, Hayden MJ (2017) Detection and validation of genomic regions associated with resistance to rust diseases in a worldwide hexaploid wheat landrace collection using BayesR and mixed linear model approaches. Theor Appl Genet 130:777–793. https://doi.org/10.1007/s00122-016-2851-7
Peterson RF, Campbell AB, Hannah AE (1948) A diagrammatic scale for estimating rust intensity on leaves and stems of cereals. Can J of Res 26:496–500. .doi.org/10.1139/cjr48c-033
Prins R, Pretorius Z, Bender C, Martin A (2011) QTL mapping of stripe, leaf and stem rust resistance genes in a Kariega × Avocet S doubled haploid wheat population. Mol Breed 27::259–270. doi:10.1007/s11032-010-9428-y
R Core Team (2019) R: A language and environment for statistical computing. R Foundation for Statistical Computing. Vienna, Austria. http://www.R–project.org/
Rahmatov M (2016) Genetic characterization of novel resistance to stem rust and stripe rust in wheat-alien introgression lines. Dissertation, Acta Universitatis Agriculturae Sueciae
Randhawa MS, Bariana HS, Mago R, Bansal UK (2015) Mapping of a new stripe rust resistance locus Yr57 on chromosome 3BS of wheat. Mol Breed 35. doi:10.1007/s11032-015-0270-0
Rao NK (2004) Plant genetic resources: advancing conservation and use through biotechnology. Afr J biotechnol 3:136–145. doi:10.5897/AJB2004.000-2025
Ren Y, Li Z, He Z, Wu L, Bai B, Lan C, Wang C, Zhou G, Zhu H, Xia X (2012) QTL mapping of adult-plant resistances to stripe rust and leaf rust in Chinese wheat cultivar Bainong 64. Theor Appl Genet 125::1253–1262. doi:10.1007/s00122-012-1910-y
Ren RS, Wang MN, Chen XM, Zhang ZJ (2012) Characterization and molecular mapping of Yr52 for high temperature adult-plant resistance to stripe rust in spring wheat germplasm PI 183527. Theor Appl Genet 125:847–857. doi:10.1007/s00122-012-1877-8
Riaz A, Athiyannan N, Periyannan SK, Afanasenko O, Mitrofanova OP, Platz GJ, Aitken EAB, Snowdon RJ, Lagudah ES, Hickey LT, Voss-Fels KP (2018) Unlocking new alleles for leaf rust resistance in the Vavilov wheat collection. Theor Appl Genet 131:127–144. doi:10.1007/s00122-017-2990-5
Riaz A, Athiyannan N, Periyannan S, Afanasenko O, Mitrofanova O, Aitken EAB, Lagudah ET, Hickey L (2016) Mining Vavilov's treasure chest of wheat diversity for adult plant resistance to Puccinia triticina. Plant Dis 101:317–323. doi:10.1094/PDIS-05-16-0614-RE
Riaz A, Hathorn A, Dinglasan E, Ziems L, Richard C, Singh D, Mitrofanova O, Afanasenko O, Aitken E, Godwin I, Hickey L (2016) Into the vault of the Vavilov wheats: old diversity for new alleles. Genet Resour Crop Evol 64:531–544. doi.org/10.1007/s10722-016-0380-5
Rinaldo A, Gilbert B, Boni R, Krattinger SG, Singh D, Park RF, Lagudah E, Ayliffe M (2017) The Lr34 adult plant rust resistance gene provides seedling resistance in durum wheat without senescence. Plant Biotechnol J 15:894–905. doi:10.1111/pbi.12684
Risk JM, Selter LL, Krattinger SG, Viccars LA, Richardson TM, Buesing G, Herren G, Lagudah ES, Keller B (2012) Functional variability of the Lr34 durable resistance gene in transgenic wheat. Plant biotechnol J 10:477–487. doi.org/10.1111/j.1467-7652.2012.00683.x
Rosewarne GM, Herrera-Foessel SA, Singh RP, Huerta-Espino J, Lan CX, He ZH (2013) Quantitative trait loci of stripe rust resistance in wheat. Theor Appl Genet 126:2427–2449. doi.org/10.1007/s00122-013-2159-9
Sadovaya AS, Gultyaeva EI, Mitrofanova OP, Shaidayuk EL, Hakimova AG, Zuev EV (2015) Leaf rust resistance in common wheat varieties and lines from the collection of the Vavilov Plant Industry Institute carrying alien genetic material. Russ J Genet Appl Res 5:233–241. doi:10.1134/S2079059715030144
Saintenac C, Jiang D, Wang S, Akhunov E (2013) Sequence-based mapping of the polyploid wheat genome. G3 3:1105–1114. doi: 10.1534/g3.113.005819
Sanin SS, Nazarova LN (2010) Phytosanitary situation on wheat crops in the Russian Federation (1991-2008). Plant Protection Quarantine 2:69–80 (in Russ.)
Santra DK, Chen XM, Santra M, Campbell KG, Kidwell KK (2008) Identification and mapping QTL for high-temperature adult-plant resistance to stripe rust in winter wheat (Triticum aestivum L.) cultivar 'Stephens'. Theor Appl Genet 117:793–802. doi:10.1007/s00122-008-0820-5
Singh A, Knox RE, Depauw RM, Singh AK, Cuthbert RD, Campbell HL, Shorter S, Bhavani S (2014) Stripe rust and leaf rust resistance QTL mapping, epistatic interactions, and co-localization with stem rust resistance loci in spring wheat evaluated over three continents. Theor Appl Genet 127:2465–2477. doi:10.1007/s00122-014-2390-z
Singh R, Nelson J. Sorrells M (2000) Mapping Yr28 and other genes for resistance to stripe rust in wheat. Crop Sci 40:1148–1155. .doi.org/10.2135/cropsci2000.4041148x
Singh R, Sukhwinder S, Bhavani S, Basenet BR (2014) Progress towards genetics and breeding for minor genes based resistance to Ug99 and other rust in CIMMYT high yeilding. J Integ Agric. doi.org/10.1016/S2095-3119(13)60649-8
Somers DJ, Isaac P, Edwards K (2004) A high-density microsatellite consensus map for bread wheat (Triticum aestivum L.). Theo Appl Genet 109:1105–1114. doi.org/10.1007/s00122-004-1740-7
Song QJ, Shi JR, Singh S, Fickus EW, Costa JM, Lewis J, Gill BS, Ward R, Cregan PB (2005) Development and mapping of microsatellite (SSR) markers in wheat. Theor Appl Genet 110:550–560. doi:10.1007/s00122-004-1871-x
Sorensen CK, Hovmoller MS, Leconte M, Dedryver F, DeVallavieille-Pope C (2014) New races of Puccinia striiformis found in europe reveal race specificity of long-term effective adult plant resistance in wheat. Phytopathology 104:1042–1051. doi.org/10.1094/PHYTO-12-13-0337-R
Stakman EC, Stewart DM, Leogering WQ (1962) Identification of physiological races of Puccinia graminis VAR tritici. United States Department of Agriculture, ARS E617, p53
Sthapit Kandel J, Krishnan V, Jiwan D, Chen X, Skinner DZ, See DR (2017) Mapping genes for resistance to stripe rust in spring wheat landrace PI 480035. PLoS One 12. doi.org/10.1371/journal.pone.0177898
Sthapit J, Newcomb M, Bonman JM, Chen X, See DR (2014) Genetic diversity for stripe rust resistance in wheat landraces and identification of accessions with resistance to stem rust and stripe rust. Crop Sci 54:2131. doi:10.2135/cropsci2013.07.0438
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. doi:10.1094/PHYTO.2003.93.7.881
Talas F, Longin F, Miedaner T (2011) Sources of resistance to Fusarium head blight within Syrian durum wheat landraces. Plant Breed 130:398–400. doi.org/10.1111/j.1439-0523.2011.01867.x
Uauy C, Brevis JC, Chen XM, Khan I, Jackson L, Chicaiza O, Distenfeld A, Fahima T, Dubcovsky J (2005) High-temperature adult-plant stripe rust resistance gene Yr36 from Triticum turgidum ssp. dicoccoides is closely linked to the grain protein content locus Gpc-B1. Theor Appl Genet 112:97–105. doi: 10.1007/s00122-005-0109-x
Vazquez MD, Zemetra R, Peterson CJ, Chen XM, Heesacker A, Mundt CC (2015) Multi-location wheat stripe rust QTL analysis: genetic background and epistatic interactions. Theor Appl Genet 128::1307–1318. doi:10.1007/s00122-015-2507-z
Voorrips RE (2002) )MapChart: Software for the graphical presentation of linkage maps and QTLs. J Hered 93:77–78. doi.org/10.1093/jhered/93.1.77
Wan A, Chen X (2014) Virulence characterization of Puccinia striiformis f. sp. tritici using a new set of Yr single-gene line differentials in the United States in 2010. Plant Dis 98:534–1542. doi:10.1094/PDIS-01-14-0071-RE
Wang S, Wong D, Forrest K, Allen A, Chao S, Huang BE, Maccaferri M, Salvi S, Milner SG, Cattivelli L (2014a) Characterization of polyploid wheat genomic diversity using a high_density 90 000 single nucleotide polymorphism array. Plant biotechnol J 12:787–796. doi:10.1111/pbi.12183
Wang H, Zou S, Li Y et al (2020) An ankyrin-repeat and WRKY-domain-containing immune receptor confers stripe rust resistance in wheat. Nat Commun 11:1353. https://doi.org/10.1038/s41467-020-15139-6
Warnes G, Gorjanc G, LeischF, Man M (2012) Population Genetics.R Package Version 1.3.8.1. Available from: http://CRAN.R-project.org/package=genetics
Wellings CR (2007) Puccinia striiformis in Australia: a review of the incursion, evolution, and adaptation of stripe rust in the period 1979-2006. Aust J Agric Res 58:567–575. doi.org/10.1071/AR07130
Wellings CR, Wright DG, Keiper F, Loughman R (2003) First detection of wheat stripe rust in Western Australia: evidence for a foreign incursion. Australas Plant Pathol 32:321–322. doi.org/10.1071/AP03023
Woloshen V, Huang S, Li X (2011) RNA-binding proteins in plant immunity. J Pathog doiorg/. 10.4061/2011/278697
Wu XL, Wang JW, Cheng YK, Ye XL, Li W, Pu ZE, Jiang QT, WeI YM, Deng M, Zheng YL, Chen GY (2016) Inheritance and molecular mapping of an all-stage stripe rust resistance gene derived from the Chinese common wheat landrace "Yilongtuomai". J Hered 107:463–470. doi:10.1093/jhered/esw032
Xia C, Wan A, Wang M, Jiwan DA, See DR, Chen X (2016) Secreted protein gene derived-single nucleotide polymorphisms (SP-SNPs) reveal population diversity and differentiation of Puccinia striiformis f. sp. tritici in the United States. Fungal Biol 120:729–744. doi:10.1016/j.funbio.2016.02.007
Xu LS, Wang MN, Cheng P, Kang Z, Hulbert S, Chen X (2013) Molecular mapping of Yr53, a new gene for stripe rust resistance in durum wheat accession PI 480148 and its transfer to common wheat. Theor Appl Genet 126:523–533. doi:10.1007/s00122-012-1998-0
Yang EN, Rosewarne GM, Herrera-Foessel SA, Huerta-Espino J, Tang ZX, Sun CF, Ren ZL, Singh RP (2013a) QTL analysis of the spring wheat "Chapio" identifies stable stripe rust resistance despite inter-continental genotype x environment interactions. Theor Appl Genet 126:1721–1732. doi:10.1007/s00122-013-2087-8
Yang J, Liu Y, Yan H, Tian T, You Q, Zang L, Xu W, Su Z (2018) Plant EAR: Functional analysis platform for Plant EAR motif-containing proteins. Front Genet 9:590. doi=10.3389/fgene.2018.00590
Yang X, Pengcheng L, Zefeng Y, Chenwu X (2017) Genetic mapping of quantitative trait loci in crops
Crop J 5:175–184. doi.org/10.1016/j.cj.2016.06.003
Yu J, Pressoir G, Briggs WH, Vroh Bi I, Yamasaki M, Doebley JF, McmulleN MD, Gaut BS, Nielsen DM, Holland JB, Kresovich S, Buckler ES (2006) A unified mixed-model method for association mapping that accounts for multiple levels of relatedness. Nat Genet 38:203–208. doi:10.1038/ng1702
Zahravi M, Balakrishna PV, Banks PM, Bariana HS, Ghannadha MR, Pogna NE, Romano M, Pogna EA, Galterio G, Shariflou MR (2003) Bulk segregant analysis of stripe rust resistance in wheat (Triticum aestivum) using microsatellite markers. Instituto Sperimentale per la Cerealicoltura, Rome, pp 861–863
Zegeye H, Rasheed A, Makdis F, Badebo A, Ogbonnaya FC (2014) Genome-wide association mapping for seedling and adult plant resistance to stripe rust in synthetic hexaploid wheat. PLoS One 9:e105593. .doi.org/10.1371/journal.pone.0105593
Zhang C, Huang L, Zhang H et al(2019)An ancestral NB-LRR with duplicated 3′UTRs confers stripe rust resistance in wheat and barley. Nat Commun 10, 4023. doi.org/10.1038/s41467-019-11872-9
Zhou XL, Wang WL, Wang LL, Hou DY, Jing JX, Wang Y, Xu ZQ, Yao Q, Yin JL, Ma DF (2011) Genetics and molecular mapping of genes for high-temperature resistance to stripe rust in wheat cultivar Xiaoyan 54. Theor Appl Genet 123:431–438. doi:10.1007/s00122-011-1595-7
Zhou XL, Wang MN, Chen XM, Lu Y, Kang ZS, Jing JX (2014) Identification of Yr59 conferring high-temperature adult-plant resistance to stripe rust in wheat germplasm PI 178759. Theor Appl Genet 127:935–945. doi:10.1007/s00122-014-2269-z
Ziems LA, Hickey LT, Hunt CH, Mace ES, Platz GJ, Franckowiak JD, Jordan DR (2014) Association mapping of resistance to Puccinia hordei in Australian barley breeding germplasm. Theor Appl Genet 127:1199–1212. doi:10.1007/s00122-014-2291-1

Download PDF

Journal Publication

published 03 Feb, 2022

Read the published version in Theoretical and Applied Genetics →

Editorial decision: Minor revisions
24 Nov, 2021
Reviews received at journal
21 Oct, 2021
Reviewers invited by journal
21 Oct, 2021
Editor assigned by journal
21 Oct, 2021
First submitted to journal
20 Oct, 2021

You are reading this latest preprint version

Mining the Vavilov wheat diversity panel for new sources of adult plant resistance to stripe rust

Status:

Journal Publication

Version 1

Abstract

Figures

Key message

Introduction

Materials And Methods

Plant material

Pathogen

PCR marker screening for known APR genes

Phenotyping for seedling resistance

Phenotyping for adult plant resistance

Genotyping

Data analysis

Genome-wide association studies

Comparison of QTL with previously reported Yr genes (QTL) and identification of putative candidate genes

Assigning and stacking of alleles associated with novel APR QTL

Results

Analysis of stripe rust response at seedling and adult plant stages

Cluster analysis using phenotype and passport data

Identification of novel sources of APR

Marker‒trait associations

Identification of candidate genes and their functional annotation

Allele stacking effects and geographical distribution of the novel APR QTL

Discussion

Vavilov wheat accessions offer diverse sources of stripe rust resistance

GWAS reveals 47 genomic regions associated with resistance

Novel QTL for adult plant resistance

Novel APR QTL confer additive effects to reduce disease

Conclusion

Declarations

References

Supplementary Files

Status:

Journal Publication

Version 1