Next Article in Journal
Complete Chloroplast Genome of Alternanthera sessilis and Comparative Analysis with Its Congeneric Invasive Weed Alternanthera philoxeroides
Next Article in Special Issue
A Genome-Wide Association Study Approach to Identify Novel Major-Effect Quantitative Trait Loci for End-Use Quality Traits in Soft Red Winter Wheat
Previous Article in Journal
Downregulation of Histone H3.3 Induces p53-Dependent Cellular Senescence in Human Diploid Fibroblasts
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Genes
Volume 15
Issue 5
10.3390/genes15050542
genes-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Comparative Analysis of Virulence and Molecular Diversity of Puccinia striiformis f. sp. tritici Isolates Collected in 2016 and 2023 in the Western Region of China

by
Tesfay Gebrekirstos Gebremariam
1,2,
Fengtao Wang
3,
Ruiming Lin
3,* and
Hongjie Li
1,4,*
1
The National Engineering Laboratory of Crop Molecular Breeding, Institute of Crop Sciences, Chinese Academy of Agricultural Sciences, Beijing 100081, China
2
Tigray Agricultural Research Institute, Mekelle P.O. Box 492, Ethiopia
3
State Key Laboratory for Biology of Diseases and Insect Pests, Institute of Plant Protection, Chinese Academy of Agricultural Sciences, Beijing 100193, China
4
Institute of Biotechnology, Xianghu Laboratory, Hangzhou 311231, China
*
Authors to whom correspondence should be addressed.
Genes 2024, 15(5), 542; https://doi.org/10.3390/genes15050542
Submission received: 8 March 2024 / Revised: 15 April 2024 / Accepted: 20 April 2024 / Published: 25 April 2024
(This article belongs to the Special Issue Quality Gene Mining and Breeding of Wheat)
Browse Figures
Versions Notes

Abstract

:
Puccinia striiformis f. sp. tritici (Pst) is adept at overcoming resistance in wheat cultivars, through variations in virulence in the western provinces of China. To apply disease management strategies, it is essential to understand the temporal and spatial dynamics of Pst populations. This study aimed to evaluate the virulence and molecular diversity of 84 old Pst isolates, in comparison to 59 newer ones. By using 19 Chinese wheat differentials, we identified 98 pathotypes, showing virulence complexity ranging from 0 to 16. Associations between 23 Yr gene pairs showed linkage disequilibrium and have the potential for gene pyramiding. The new Pst isolates had a higher number of polymorphic alleles (1.97), while the older isolates had a slightly higher number of effective alleles, Shannon’s information, and diversity. The Gansu Pst population had the highest diversity (uh = 0.35), while the Guizhou population was the least diverse. Analysis of molecular variance revealed that 94% of the observed variation occurred within Pst populations across the four provinces, while 6% was attributed to differences among populations. Overall, Pst populations displayed a higher pathotypic diversity of H > 2.5 and a genotypic diversity of 96%. This underscores the need to develop gene-pyramided cultivars to enhance the durability of resistance.

1. Introduction

Resistance to rust diseases in bread wheat cultivars (Triticum aestivum L., 2n = 6× = 42, AABBDD) is a common challenge due to the constant evolution of virulence in rust pathogens. Pst causes stripe rust (also known as yellow rust) on more than 20 grass and crop species including wheat, barley (Hordeum vulgare L.), triticale (×Triticosecale Wittmack), and rye (Secale cereale L.) [1]. Even though the extent of yield loss caused by stripe rust depends on factors such as pathogen aggressiveness, cultivar resistance level, and environmental variables, it has been reported that the global annual yield loss is estimated to be 1%, equivalent to USD 1 billion [2,3]. Globally, stripe rust is highly prevalent, with an occurrence rate of 88%. Moreover, the majority of wheat cultivars exhibit compatibility with this disease due to the presence of susceptible genes [3,4]. The prevalence of the pathogen is closely associated with the availability of suitable host plants and favorable weather conditions [5].
Wheat stripe rust can be managed through the utilization of cultivars carrying resistance genes, the application of fungicides, and the implementation of improved agronomic practices. Among these strategies, the deployment of resistant Yr genes stands out as the most effective approach due to its environmental friendliness, cost-effectiveness, and ease of implementation [6]. At present, over 80 Yr genes have been identified from various domesticated and wild species of wheat, comprising 54 genes conferring all-stage resistance (ASR) and 24 genes providing adult-plant resistance (APR) or high-temperature adult-plant (HTAP) resistance [7].
The diverse genetic makeup of Pst populations, resulting from mechanisms such as sexual recombination, mutation, and somatic hybridization, is responsible for virulence variation within the pathogen of stripe rust [8,9,10]. The wind-mediated dispersal of Pst urediniospores across wheat-growing environments also contributes to the genetic diversity of the pathogen [11]. Over time, these diversifications facilitate the emergence or re-emergence of highly virulent races of Pst, which are capable of overcoming genetic resistance in wheat cultivars. Subsequently, this results in recurrent disease epidemics [3,9,12]. For instance, the emergence of new races within the Pst population has led to disease outbreak, due to the breakdown of resistance for genes Yr7, Yr8, Yr9, Yr10, Yr17, and Yr27 in various countries [12,13,14]. High levels of virulence have been reported in numerous characterized Yr genes, with the exceptions of Yr5 and Yr15 [7,15]. The genetic variabilities and virulence spectra of Pst populations have been extensively studied across countries, including Australia [11], South Africa [16], Canada [14], Russia [1], the USA [15,17,18], China [19], and Turkey [20].
From 1950 to 2000, China witnessed seven significant epidemics of stripe rust [21]. Since then, numerous research studies have been undertaken to monitor the virulence dynamics and genetic diversity of Pst populations. The analysis of virulence variation in Pst populations utilized a set of 19 differentials of wheat cultivars established in 2010, along with Yr single-gene differentials [7]. Chen et al. [22] identified 41 races of Pst from a total of 4714 isolates and they also observed changes in the virulence of CYR33 (Chinese yellow rust) on YrSu. Additionally, the race composition of a Pst population across various provinces of the country was determined [23,24]. In the hotspot provinces of China, more than 750 Pst races/pathotypes have been identified. These include various predominant and virulent races, namely CYR1, CYR8, CYR10, CYR13, CYR16, CYR17, CYR18, CYR19, CYR23, CYR25, CYR28, CYR29, CYR30, CYR33, and the L13 and SY11 pathotype groups [7]. The implementation of molecular markers for genotyping Pst populations has been widely employed in China. Using SSR markers, Chen et al. [24] identified shared multi-locus genotypes with the Pst populations in the provinces of Yunnan and Guizhou, which demonstrated a significant genetic resemblance between these geographical areas. Hu et al. [25] conducted a study to assess the genetic differentiation among 961 Pst isolates, obtained from various wheat-growing regions in the provinces of Gansu, Shaanxi, Sichuan, and Tibet. Their research findings explicitly indicated the existence of gene flow between Pst populations. Similarly, Yan et al. [26] confirmed genetic differentiation within individuals of isolates, while a small proportion of genetic variation was observed between the Pst populations of Yunnan and the Yangtze river. Furthermore, the simple nucleotide polymorphism (SNP) and amplified fragment length polymorphism (AFLP) techniques are commonly utilized in China for genotyping Pst populations [27,28]. In China, the provinces of Gansu, Sichuan, Guizhou, and Yunnan have been identified as epidemiological zones, exhibiting frequent variations in the virulence of stripe rust [21].
Understanding the temporal and spatial dynamics of Pst populations is essential for the implementation of timely and effective disease management strategies in these high-risk regions. In this study, we characterized and compared the virulence patterns of two groups of Pst isolates collected during the epidemic seasons of 2016 and 2023. The specific objectives of this study were to (1) analyze the pathotypes of Pst collections based on the Chinese differentials, (2) determine the frequency of virulence and effective genes using Yr single-gene differentials, and (3) assess the genetic diversity among Pst populations using molecular markers.

2. Materials and Methods

2.1. Collection of Pst Isolates

Wheat leaves exhibiting symptoms of wheat stripe rust were systematically collected from experimental plots, nurseries, and commercial fields located across the provinces of Yunnan, Guizhou, Sichuan, and Gansu. To obtain an adequate quantity of urediniospores and a representative sample of isolates, a total of five infected leaves exhibiting Pst uredinia from a single cultivar or field were collected and considered as a single sample. Two rounds of field surveys were conducted during the cropping seasons of 2023 and 2016. For the comparative analysis of virulence and genotypic variation, isolates collected in the 2016 cropping season were designated as old isolates of Yunnan (YNOI), Guizhou (GZOI), Sichuan (SCOI), and Gansu (GSOI), while those collected in the 2023 cropping season were designated new isolates of Yunnan (YNNI), Guizhou (GZNI), and Sichuan (SCNI). Leaf samples were dried with moisture absorbent for 24 h and stored at 4 °C until further use.

2.2. Isolate Identification and Multiplication

To obtain fresh urediniospores, 3–5 cm segments of infected leaves from each sample were placed on wet absorbent paper in Petri dishes and were incubated in a dark room at 10 °C for 24 h [29,30]. Fresh urediniospores produced by a single uredinium were inoculated onto 10-day-old susceptible seedlings of Mingxian 169 (MX169), using a sterilized needle for isolate multiplication. The MX169 seedlings were treated with a 0.33% maleic hydrazide solution to inhibit further growth [31]. Urediniospores were harvested 14–18 days after inoculation and were dried in a desiccator at 5 °C for 24–48 h. Urediniospores produced from a single uredinium on each leaf were considered as a single isolate. In total, 143 isolates of Pst were recovered from 160 samples. Pst isolates were temporarily stored at 4 °C, or at −40 °C, for months. Long-term stored isolates were heat-shocked at 42 °C for 5 min to reactivate and were then stored at 4 °C for 12 h.

2.3. Inoculation of Chinese Differentials and Yr Single Gene Lines

Nineteen Chinese differentials and 31 single Yr gene differentials were planted in pots (10 × 10 × 7 cm3) and boxes (10 × 25 × 50 cm3), respectively (Supplementary Table S1). MX169 was included as the susceptible check. Growth conditions and inoculation procedures were previously described [29]. Seedlings were grown in a growth chamber, programmed at 15 ± 1 °C/12 ± 1 °C with a photoperiod of 16 h light/8 h dark, maintaining a relative humidity of 70%. Seedlings at the two-leaf stage were sprayed with 15–20 mg of urediniospores suspended in 1 mL of Novec 7100 engineered fluid (3M, Maplewood, MN, USA). Inoculated seedlings were then placed in a dark dew chamber at 10 °C for 24 h, before being transferred to growth chambers automated with a diurnal temperature cycle ranging from 4 °C at 2:00 AM to 20 °C at 2:00 PM, 8 h dark and 16 h light from 10 PM to 6 AM and from 6 AM to 10 PM, respectively. Stringent measures were taken to prevent cross-contamination during the inoculation, incubation, and handling of isolates.

2.4. Characterization of Virulence

The infection types (ITs) produced by Pst isolates on seedlings of Chinese differentials and Yr single-gene lines were assessed after complete sporulation of MX169, using a previously described 0–4 scale [22], where 0 = absence of visible symptoms; 0; = presence of necrotic and/or chlorotic flecks, without sporulation; 1 = presence of necrotic and/or chlorotic blotches or stripes, with trace sporulation; 2 = presence of necrotic and/or chlorotic blotches or stripes, with moderate sporulation; 3 = presence of chlorotic blotches or stripes, with moderate to abundant sporulation; and 4 = absence of necrosis or chlorosis, with abundant sporulation. ITs 0–2 were considered as pathogen avirulence, while ITs 3–4 were deemed as pathogen virulence. Isolates lacking clear virulence patterns on the testing materials were subjected to re-inoculation to verify their virulence patterns.

2.5. Pathotype Designation and Effectiveness of Yr Genes

Pst pathotypes were identified and characterized based on the virulence formulae of the 19 Chinese differentials, using the methodology outlined by Draz [32]. For instance, pathotype V2, 4, 8, 14 exhibited virulence on the corresponding differentials of Fulhard (unknown), Mentana (unknown), Funo (YrA+), and Suwon 11 (YrSu). Pst isolates displaying similar virulence patterns on the differentials were grouped together as a single pathotype, while those displaying unique virulence patterns were classified as distinct pathotypes. The efficacy levels of the Yr genes were classified as high, moderate, or low, based on their respective virulence frequencies of 0–20%, 20–50%, and 50–100%, as previously described [15].

2.6. Association of Yr Genes

Pairwise virulence/avirulence associations of Yr loci were empirically assessed according to Parks et al. [33]. Four pairwise groups of alleles at two loci were identified, as follows: isolates virulent to both Yr genes (VV); isolates avirulent to both Yr genes (AA); isolates virulent to one Yr gene, but avirulent to the other (VA); and isolates avirulent to one Yr gene, but virulent to the other (AV). Positive associations were determined if the AA or VV groups represented more than 60% of the tested isolates, while negative associations were determined if AV or VA accounted for more than 60% of the isolates. If the combined frequency of AA and VV fell within the range of 40–60%, the association was deemed unclear.

2.7. Molecular Marker Analysis

DNA was extracted from 20 mg of dried urediniospores from each isolate, using the cetyl trimethyl ammonium bromide method [34]. Prior to PCR amplification, the quality and quantity of the extracted DNA were determined using an ND-1000 spectrophotometer (Thermo Scientific, Waltham, MA, USA), followed by dilution with ddH2O to a concentration of 50 ng/μL. The genotyping of DNA samples was conducted using eleven previously designed SSR markers including RJ3 and RJ22 [35]; RJ3N and RJ13N [36]; CPS08, CPS09, CPS13, and CPS27 [37]; and PstP002, PstP003, and PstP006 [38]. These markers were chosen for their high informativeness. Primers of these markers were synthesized by Sangon Biotech Inc. (Shanghai, China, http://www.sangon.com/, accessed on 10 February 2024).
A reaction mixture consisted of 25 μL, comprising 2 μL of DNA (50 ng/μL), 2.5 μL of 10× reaction buffer (Mg2+-free), 2 μL of Mg2+ (25 mM), 2 μL of dNTPs (2.5 mM), 1 μL of each primer (10 mM), 0.2 μL of Taq DNA polymerase (5 U/μL), and 14.3 μL of ddH2O. The amplification of DNA was carried out in a thermocycler (Bio-Rad, Los Angelas, CA, USA) under the following conditions: initial denaturation at 94 °C for 5 min, 30 cycles of denaturation at 94 °C for 30 s, annealing at 50–58 °C for 30 s, extension at 72 °C for 1 min, and a final extension at 72 °C for 10 min. The amplified products were resolved using capillary electrophoresis technology [39].

2.8. Data Analysis

The diversity of pathotypes in the Pst population was assessed using the Shannon–Wiener diversity index (H) [40]. According to this index, low, moderate, and high diversities were indicated with values of H < 1.5, 1.5 < H < 2.5, and H > 2.5, respectively.
The data scored using a 0–4 scale were transformed into a coding system, where a value of 1 represents ratings of 3–4, indicating virulence, while a value of 0 represents ratings of 0–2, indicating avirulence. The virulence frequency (VF) indicates the percentage of virulent isolates of Pst in a population to an individual Yr gene line. The virulence complexity (VC) of each isolate indicates the isolate’s potential to overcome host Yr genes. VF and VC were analyzed using a Jaccard dissimilarity coefficient in the virulence analysis tool (VAT) [41]. The statistical significance of virulence frequencies among subpopulations of Pst was assessed using the t-test (https://www.scribbr.com/statistics/t-test/, accessed on 2 March 2024). The unique virulence patterns (UVPs), average virulence complexity (AVC), relative virulence complexity (RVC), and virulence uniformity (VU) of the Pst isolates were also computed using VAT [41].
A dendrogram showing the virulence patterns of both old and new Pst isolates was generated using the unweighted pair group method with arithmetic mean (UPGMA), in the SAHN program of NTSYSpc 2.2 [42]. Genetic diversity parameters, including number of alleles (Na), number of effective alleles (Ne), Shannon index (I), unbiased diversity (uh), percentage of polymorphic loci (PPL), and analysis of molecular variance (AMOVA) with 999 permutations, were analyzed using GenAlEx 6.5 [43]. The Nei’s genetic distance and principal coordinate analysis (PCoA) plots were generated using GenAlEx 6.5 [43].

3. Results

3.1. Identification of Pathotypes and Diversities in the New and Old Isolates of Pst

The analysis of virulence patterns among 143 Pst isolates on the 19 Chinese wheat differentials for stripe rust revealed 98 distinct pathotypes (Supplementary Table S2). Nineteen Pst pathotypes were detected more than once, with virulence frequencies ranging from 1.4% to 5.0% (Table 1). The most frequent pathotype was V2, 4, 8, 14, with a pathotype frequency of 5.0% (7 out of 143) in the Guizhou, Sichuan, and Yunnan provinces. This pathotype exhibited virulence on the differential lines Fulhard, Mentana, Funo, and Suwon 11, and was avirulent on all other differentials including Trigo Eureka, Lutescens 128, Virgilio, Abbondanza, Early Premium, Danish 1, Jubilejina 2, Fengchan 3, Lovrin 13, Kangyin 655, Zhong 4, Lovrin 10, Hybrid 46, T. spelta var album, and Guinong 22. The second most frequent pathotype, V2, 4, 8, 12, 14, 16, had a frequency of 4.2% in the Sichuan and Gansu Pst populations. Pathotypes V2, 8, 14, and V0 were each detected five times (3.5%). The former was found in both new and old isolate groups of Yunnan, Guizhou, and Sichuan, while the latter was detected in the old isolate groups of the Yunnan, Guizhou, and Gansu provinces. Three pathotypes (V2, 14; V2, 8, 14, 19; and V2, 8) each had frequencies of 2.8% across the old isolates of the provinces. Pathotypes V2, 7, 8, 14, 19; V2, 8, 19; and V2 were detected, each at 1.4%, across old isolates of the provinces. Pathotypes V4, 14, 16 and V2, 4, 14, 16 had virulence frequencies of 2.8% and 2.1% in the new isolate groups of the Guizhou and Sichuan provinces, respectively. Similarly, V2, 3, 4, 6, 7, 8, 10, 11, 12, 14, 16, 19; V2 10; V2, 4, 8, 11, 14, 16; and V2, 8, 14, 16 were detected two times (1.4%) each across new isolates of the provinces. Six pathotypes (V2, 4, 8, 14; V2, 8, 14; V2, 4, 8, 10, 14, 19; V3, 7, 10, 19; V2, 4, 8, 12, 14, 16; and V2, 4, 7, 8, 14) with frequencies ranging from 1.4% to 5.0% were detected in both isolate groups. The remaining 79 pathotypes were each detected only once (Supplementary Table S2).
Pathotype diversity was comparable in Yunnan and Gansu (each with H = 3.6), as well as in Sichuan (H = 3.5). However, pathotype diversity was lower in the Guizhou Pst population, with an H value of 2.8 (Table 1). These data provide strong evidence for the higher diversity and uniform distribution of pathotypes across the Pst populations.

3.2. Virulence Complexities and Frequencies on the Chinese Differentials

The virulence complexities (VCs) of old pathotypes (derived from old isolates) ranged from 0 to 16, with a mean of 5.2 (Figure 1). Pathotype V0 displayed avirulence to all differential sets. Conversely, pathotype V1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 16, 17 was virulent to all differentials, except Zhong 4, T. spelta var album, and Guinong 22. In new pathotypes (derived from new isolates), the VCs ranged from 1 to 12, with a mean of 7. Pathotype V12 was solely virulent to Lovrin 13, while pathotype V2, 3, 4, 5, 7, 8, 9, 10, 11, 12, 14, 16 showed virulence to Fulhard, Lutescens 128, Mentana, Virgilio, Early Premium, Funo, Danish 1, Jubilejina 2, Fengchan 3, Lovrin 13, Suwon 11, and Lovrin 10. In both pathotype groups, narrow virulence (e.g., 0–1) and wide virulence (e.g., 11–16) were limited. Instead, virulence between 2 and 10 differentials tended to be frequently observed among the Pst pathotypes (Supplementary Table S2 and Figure 1).
Virulence frequencies in the new isolates were higher, ranging from 56 to 100% for Fulhard, 50 to 80% for Mentana, 56 to 83% for Funo, and 82 to 90% for Suwon11, compared to 78–86%, 40–54%, 60–80%, and 54–71% in the old isolates, respectively (Figure 2). Virulence to Lutescens 128 and Lovrin 10 was as high as 70% and 50% in the new isolates, respectively, while it was as low as 39% and 43% in the old isolates. Virulence was moderate, though higher in the new isolates for Abbondanza (9–28%), Early Premium (33–46%), Jubilejina 2 (17–39%), Lovrin 13 (28–44%), and Hybrid 46 (6–20%), compared to 9–18%, 18–32%, 20–32%, 9–28%, and 4–10% in the old isolates, respectively. Only Guinong 22 exhibited a higher virulence (50%) to the old isolates than to the new ones (33%). Virulence to Trigo Eureka, Virgilio, Danish 1, Fengchan 3, and Kangyin 655 was less than 20% across isolate groups. Neither isolate group exhibited virulence to Zhong 4 and T. spelta var album.

3.3. Virulence Frequencies and Complexities on Yr Single-Gene Differentials

The frequencies of virulence to 31 Yr single-gene differentials ranged from 0 to 93% (Table 2). All Pst isolates were avirulent to Yr5 and Yr15. Low frequencies of virulence were observed for Yr24 (10–20%) and Yr26 (0–20%); moderate frequencies to Yr76 (18–50%), Yr50 (0–42%), Yr41 (0–35%), Yr32 (0–47%), Yr10 (0–41%), Yr17 (13–50%), YrSp (0–40%), and YrTr1 (0–32%); and high frequencies to Yr45 (63–92%), YrSu (58–82%), Yr43 (57–75%), Yr25 (50–90%), YrJu4 (55–88%), and YrRes (58–93%) across the provinces. Some Yr genes showed province-specific effectiveness. Frequency of virulence to Yr21 exceeded 55% in the Yunnan, Guizhou, and Gansu Pst populations, but remained below 37% in Sichuan. Virulence to Yr1 was less than 30% in Guizhou, but ranged from 30% to 55% in the other provinces. Virulence to Yr40 was less than 40% in the Guizhou and Gansu provinces, whereas it reached as high as 72% in Sichuan and 68% in Yunnan. Virulence to Yr64 was 80% in Guizhou, but below 50% in the other provinces.
Yr genes, such as Yr50, Yr41, Yr10, Yr17, YrSp, YrTr1, Yr26, and Yr29, exhibited significantly higher frequencies of virulence to the new isolates than the old ones (p < 0.05). Virulence to Yr50 ranged between 10% and 42% in the new isolates, while an immune response was observed in the old isolates. Virulence to Yr41, Yr10, Yr17, YrSp, YrTr1, Yr26, and Yr29 was between 0 and 21%, 0 and 11%, 13 and 25%, 0 and 25%, 0 and 8%, 0 and 11%, and 30 and 55% in the old isolates. However, it increased in the new isolates, with the virulence to the corresponding genes of 30–35%, 10–41%, 26–50%, 21–40%, 20–32%, 17–20%, and 55–79%, respectively. Virulence to Yr76, Yr40, YrA, Yr1, Yr6, Yr3c, Yr7, Yr9, Yr32, Yr64, and YrRes was also substantially increased in the new isolates, despite being statistically non-significant. In contrast, Yr45, Yr25, Yr8, Yr44, and YrKy2 exhibited higher frequencies of virulence in the old isolates of Pst.
Each Pst population displayed 100% UVP (Table 3). New isolates from the Guizhou and Sichuan provinces exhibited the highest AVC, rendering each 14.4 out of 31 Yr genes ineffective. The AVC was comparable, with values of 13.6 and 13.7 observed in old isolates from Gansu and new isolates from Yunnan, respectively. AVC values of 11.5 and 12.4 were noted in the old isolates from the Guizhou and Yunnan provinces, respectively. The AVC decreased to 10.6 in the old isolates from Sichuan. A similar trend between AVC and RVC was noted among the Pst populations. The highest VU (0.36) was computed from Sichuan’s old isolates. Guizhou’s new isolate and Gansu’s old isolates exhibited the lowest VU values of 0.23 and 0.24, respectively.

3.4. Associations of Alleles at the Avirulence Loci

The pairwise association between virulence and avirulence of Pst isolates at Yr loci are presented in Table 4. Twenty-three Yr gene pairs showed a positive association. However, the association between Yr50 and YrJu4, Yr50 and YrKy2, and Yr41 and YrKy2 remained unclear, as the proportion of AA plus VV fell within the range of 40–60%, despite statistical significance.

3.5. Cluster Analysis of Virulence

Cluster analysis was performed based on the virulence data of the Pst populations using two sets of differentials (Figure 3). At a similarity coefficient of 0.14, the Pst populations were grouped into five clusters, based on the 19 sets of Chinese differentials (Figure 3A). New isolates from the Guizhou and Yunnan provinces, as well as old isolates from the Guizhou and Sichuan provinces, were clustered separately. Old isolates from Yunnan and Gansu, as well as new isolates from Sichuan were clustered independently. Old isolates from Gansu and new isolates from Sichuan are the most distantly clustered isolates. Additionally, based on the virulence frequencies of 31 single Yr genes, the Pst populations were further grouped into six clusters (Figure 3B). With the exception of the old isolates of the Guizhou and Gansu provinces, the remaining populations were independently clustered.

3.6. Molecular Diversity

The diversity parameters of the Pst populations, as determined using SSR markers, are presented in Table 5. The number of alleles (Na) was higher in the new isolates (Na = 1.97) compared to the old isolates (Na = 1.94). On the other hand, the number of effective alleles (Ne = 1.58), Shannon’s information index (I = 0.51), and unbiased diversity (uh = 0.34) were higher in the old isolates compared to the new ones. A provincial comparison revealed that the highest Na value (2.00) was observed in the Sichuan and Yunnan populations. The maximum Ne (1.56), I (0.51), and uh (0.35) values were observed in the Gansu population, followed by the Yunnan population with Ne, I, and uh values of 1.56, 0.51, and 0.34, respectively. In the Sichuan population, these values were 1.54, 0.49, and 0.33, respectively. The Guizhou population exhibited the lowest values of genetic parameters Na (1.81), Ne (1.54), I (0.46), uh (0.32), and percentage of polymorphic loci (86%). The average polymorphism value among the Pst subpopulations was 96%.
The AMOVA was used to partition the molecular variance among four provinces of the Pst populations (Yunnan, Guizhou, Sichuan, and Gansu), as well as between the old and new isolate groups. The results showed that 94% of the genetic variance existed within populations, while 6% was among the Pst populations in the four provinces. Within the isolate groups, we observed 89% of genetic variance, whereas 11% of the genetic variance was attributed to the differences between the old and new isolates of Pst (p = 0.001) (Table 6).

3.7. Genetic Distance

The genetic distance among the Pst populations in the provinces of Yunnan (YN), Guizhou (GZ), Sichuan (SC), and Gansu (GS) are presented in Table 7 and Figure 4. The Nei’s distance ranged from 0.03 to 0.10. There was a decrease in genetic distance between the populations of GS-SC, GS-YN, and SC-YN, implying that the SC-YN populations were genetically closer to each other than the other populations. The genetic distance was comparable between population pairs SC-GZ, YN-GZ, and GS-GZ, each with a value of 0.06. However, the populations of GS-SC were found to be the most genetically distant from the others. The results are consistent with the cluster analysis of virulence (Figure 3).
A principal coordinate analysis (PCoA) plot was generated based on the Nei’s distances computed between every pair of SSR genotypes. The majority of populations did not show distinct separation, except for the Guizhou population, which was predominantly clustered within the third and the fourth quadrants (Figure 4). PC1 (axis1) and PC2 (axis 2) accounted for 41% and 31% of the total genetic variance, respectively.

4. Discussion

The frequent outbreaks of wheat stripe rust in China led to the constant replacements of cultivars [21,44]. The northwestern (including Gansu) and the southwestern regions (including Yunnan, Guizhou, and Sichuan) of China are among the largest over-summering and over-wintering regions of Pst [45]. Monitoring the phenotypic and genetic diversity of Pst in these regions is essential to understand the pathogen’s population structure and to identify effective Yr genes for breeding programs. The virulence profiles of the Pst isolates on the Chinese differentials revealed a wide distribution of pathotypes across provinces. Approximately 19% of pathotypes (44% of isolates) were detected more than once within or among provinces. This phenomenon is likely attributed to the migration of spores between provinces. Pst pathotypes, including V2, 8, 14; V3, 7, 10, 19; V2, 4, 8, 14; and V2, 4, 8, 10, 14, 19, were consistently found in both old and new isolates. The presence of ancestral isolates in the current populations of Pst indicates their capacity to adapt to diverse environmental factors in the past eight years. In contrast, pathotypes such as V0; V2, 8; V2, 8, 19; and V2, 8, 14, 19 were prevalent in old isolates, but were not detected in new ones, indicating a potential genetic drift in these specific pathotypes. This result is consistent with the observation of the previous report [21].
It was observed that the newly identified isolates exhibited a higher frequency of virulence to the Chinese and Yr single-gene differentials compared to the old isolates, thereby suggesting a shift in the virulence. Two possible explanations can be derived. First, mutations, as a mechanism of virulence variation, give rise to new isolates through evolution [19]. In most cases, mutations increase the magnitude of virulence [46]. Second, the newly emerged isolates of Pst were found to be more aggressive than the old isolates [17]. However, the high virulence uniformity observed in the old isolates from Sichuan province should not be underestimated. This uniformity in virulence could increase the risk of stripe rust outbreaks, since the virulence profiles on Yr genes remain consistent. It was observed that the Pst populations collected from different provinces during the 2016 or 2023 epidemic seasons demonstrate a higher virulence similarity within their respective populations. Conversely, a lower degree of virulence similarity was observed between the Pst populations collected during the two epidemic seasons. This suggests that the changes introduced between 2016 and 2023 have a greater impact on Pst virulence dynamics compared to the environmental variations among the provinces at a given time. Previous studies have revealed that the populations of Pst from Yunnan and Guizhou are closely related, while those from Sichuan and Gansu exhibit distinct clustering patterns [47], which is consistent with our findings. Similarly, in Europe, Pst isolates collected before 2011 from the United Kingdom and France showed minimal genetic diversity. Nevertheless, these isolates were genetically distinct from those collected in 2013 [48].
With the exception of trace sporulation (IT 1) on a 0–4 scale [22], incompatibility was observed between the Pst isolates and Zhong 4 (unknown gene), Yr5, and Yr15. This suggests their effectiveness against Pst isolates. The effectiveness of Zhong 4, Yr5, and Yr15 has been extensively documented [15,29,44,49]. However, instances of virulence to Yr15 in Afghanistan [50] and Yr5 in India, Australia, Tajikistan [11], Turkey [51], and Syria [52] have been reported. In China, virulence on Yr5 was detected in the Shaanxi and Qinghai provinces [53]. Due to their close geographical proximity, the incursion of Yr5 virulent isolates through wind from Shaanxi and Qinghai to Sichuan, Guizhou, and Gansu provinces is likely to occur. Thus, continuously monitoring the Pst population structures is a crucial task. Virulence to Yr76, Yr50, Yr32, Yr24, Yr10, Yr17, YrSp, YrTr1, and Yr26 was reached at a moderate level in the studied provinces. The moderate virulence signals a potential risk of complete resistance loss under high selection pressure. Other Yr genes displayed a province-specific effectiveness. For instance, Yr6 and Yr21 exhibited a moderate level of effectiveness against the Pst population in the Sichuan province. However, their resistance was overcome by the isolates found in the remaining provinces. Yr40 was effective against the Guizhou isolates, but ineffective to the Sichuan and Yunnan isolates. This location-specific efficacy suggests the limited applicability of resistance genes across provinces. High selection pressure on Yr45, YrSu, Yr43, YrA, Yr25, Yr7, Yr9, Yr29, YrJu4, and YrRes was observed across the Pst isolate groups. Among them, Yr7, Yr9, Yr29, YrSu, and YrA have been deployed widely in Chinese wheat cultivars [54]. However, the continuous cultivation of these susceptible cultivars may lead to an accumulation of inoculum, thereby increasing disease pressure over time. Therefore, it is essential to replace these cultivars at the regional level with alternative cultivars that possess Yr5 and Yr15. Compared to the 2016 survey conducted by Chen et al. [55], dynamics in virulence were observed among Yr genes. Consistency with our results was noted for virulence to Yr5, Yr7, Yr43, and Yr44. While increases in virulence were observed for Yr8, Yr17, Yr24, Yr32, and YrTr1, Yr6, Yr9, YrSp, and Yr76 exhibited a clear decline in virulence. Likewise, our results showed a relatively lower frequency of virulence compared to previous studies in the Yunnan and Guizhou provinces [24]. This dynamic in frequency of virulence could be attributed to factors such as mutation, selection, sexual recombination, and migration [56].
Associations among alleles at pairs of avirulence loci enhance the durability of resistance against stripe rust. Several Yr gene pairs exhibited positive association, with a VV proportion of less than 10%. For instance, there was a significant association between avirulence to Yr24 and Yr26 in the Pst isolates. These Yr genes have been frequently used as sources of resistance genes in wheat breeding programs in China [57]. However, the emergence of virulent races to Yr24 and Yr26 has been reported in the Sichuan, Gansu, Shaanxi, Ningxia, and Qinghai provinces, posing a threat to wheat production [55,58]. This underscores an urgent need for pyramiding Yr24 and Yr26, as the co-selection likelihood of these Yr genes by the pathogen is low. Similarly, 22 Yr gene pairs were found to be positively associated with a lower proportion of virulence compared to their individual virulence frequency. This demonstrates that pyramiding these genes could serve as a viable strategy to extend the duration of resistant cultivars in production.
Both the phenotypic and genotypic results revealed the presence of genetic diversity in the Pst population. The comparison among provinces showed that the Pst population of Gansu exhibit a higher genetic diversity. This can be attributed to favorable environmental conditions that support over-summering and over-wintering, along with the presence of alternative hosts like Berberis spp., which encourage sexual recombination within the population [59]. Previous studies have confirmed the existence of genotypic diversity within the Pst population in Gansu province [27,60,61]. Moreover, the Yunnan Pst population exhibited a high level of pathotypic and genotypic diversity. The cultivation of various wheat genotypes throughout the year in Yunnan likely increases the genetic diversity of Pst [61]. The cultivation of various wheat cultivars frequently leads to the deployment of diverse Yr genes, thereby augmenting the overall diversity of the pathogen population. A diverse pathogen population increases the likelihood of the emergence of more virulent pathotypes through mechanisms such as mutation, recombination, and host selection [8,9,10]. Consequently, this can lead to the occurrence of novel disease epidemics. To mitigate this risk, regions with high pathogen diversity should implement a broader range of breeding programs, aiming to develop resistant cultivars that can effectively combat the disease. On the contrary, the Pst isolates from Guizhou province exhibited the lowest diversity, presumably due to unfavorable weather conditions, narrow host range, geographic barriers, and a small sample size.
The results of AMOVA revealed that the majority of genetic differentiation occurred within populations, despite limited genetic variation observed among populations. This is consistent with a previous study [25]. Several factors influence the level of genetic differentiation in Pst populations. Gene flow increases genetic variation within a population, while decreasing genetic differentiation between populations [25]. The lack of distinct clustering among Pst populations in the PCoA signifies the likelihood of gene flow among provinces. It is well-known that mutation and somatic recombination can augment genetic variation within populations [9]. Similarly, studies have confirmed that sexual recombination significantly increases the degree of genetic variability within Pst populations [8,25]. However, natural selection and genetic drift contribute to genetic differentiation between populations [19]. This suggests that the collective effect of these factors could determine the genetic variation among Pst populations.
Despite limited genetic variation observed among provinces, the Pst populations of Gansu and Sichuan were found to be genetically distant. The urediniospores of Pst exhibit a migratory behavior towards the mountainous regions of Gansu during the summer season, to undergo over-summering in favorable conditions. Similarly, during the winter season, these spores move to the lower-elevation areas of Sichuan for over-wintering, taking advantage of the cold conditions prevalent in that region [62]. Gradually, the Pst populations could potentially limit or even stop such migratory behavior, undergoing location-specific evolutionary changes to cope with temperature fluctuations. Consequently, this could result in a divergence in their genetic makeup. Recent evidence has confirmed distinct clustering patterns among Pst isolates from Sichuan and Gansu [47].
In conclusion, it can be deduced that the newly evolved isolates of Pst demonstrate a higher virulence complexity, leading to the ineffectiveness of 14.2 out of 31 Yr genes, compared to older isolates, which rendered only 12 Yr genes ineffective. All Pst isolates were avirulent to lines with Yr5, Yr15, and Zhong 4. These resistant Yr genes can be recommended for stripe rust breeding programs in China. The Pst population exhibited higher pathotypic and genotypic diversities across the studied provinces of China. These findings underscore the necessity to develop multiple resistant cultivars or employing gene-pyramiding strategies to mitigate the risks associated with the diversity of the pathogen.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/genes15050542/s1, Table S1: List of two sets of wheat differentials used to differentiate Pst pathotypes in China; Table S2: Distribution of Pst pathotypes in western provinces of China.

Author Contributions

Conceptualization, T.G.G. and R.L.; methodology, T.G.G. and F.W.; formal analysis, T.G.G.; data curation, T.G.G.; writing manuscript, T.G.G.; review and editing, R.L. and H.L.; visualization, T.G.G.; supervision, H.L.; funding acquisition, H.L. and R.L. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the Zhejiang Provincial Leading Talents Program (2024SSYS0099), the National Key Research and Development Program of China (2021YFD1401000), and the Agricultural Science and Technology Innovation Program of the Chinese Academy of Agricultural Sciences (CAAS-ASTIP-ZDRW202002).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data utilized and/or analyzed in the present study can be obtained from the corresponding author upon reasonable request.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Gultyaeva, E.; Shaydayuk, E.; Kosman, E. Virulence diversity of Puccinia striiformis f. sp. tritici in common wheat in Russian regions in 2019–2021. Agriculture 2022, 12, 1957. [Google Scholar] [CrossRef]
  2. Savary, S.; Willocquet, L.; Pethybridge, S.J.; Esker, P.; McRoberts, N.; Nelson, A. The global burden of pathogens and pests on major food crops. Nat. Ecol. Evol. 2019, 3, 430–443. [Google Scholar] [CrossRef] [PubMed]
  3. Chen, X.M. Pathogens which threaten food security: Puccinia striiformis, the wheat stripe rust pathogen. Food Secur. 2020, 12, 239–251. [Google Scholar] [CrossRef]
  4. Beddow, J.M.; Pardey, P.G.; Chai, Y.; Hurley, T.M.; Kriticos, D.J.; Braun, H.J.; Park, R.F.; Cuddy, W.S.; Yonow, T. Research investment implications of shifts in the global geography of wheat stripe rust. Nat. Plants 2015, 1, 15132. [Google Scholar] [CrossRef] [PubMed]
  5. Hovmøller, M.S.; Sørensen, C.K.; Walter, S.; Justesen, A.F. Diversity of Puccinia striiformis on cereals and grasses. Annu. Rev. Phytopathol. 2011, 49, 197–217. [Google Scholar] [CrossRef] [PubMed]
  6. Chen, X.M. Review article: High-temperature adult-plant resistance, key for sustainable control of stripe rust. Am. J. Plant Sci. 2013, 4, 608–627. [Google Scholar] [CrossRef]
  7. Wan, A.M.; Wang, X.J.; Kang, Z.S.; Chen, X.M. Variability of the stripe rust pathogen. In Stripe Rust; Chen, X.M., Kang, Z.S., Eds.; Springer: Dordrecht, The Netherlands, 2017; pp. 35–155. [Google Scholar] [CrossRef]
  8. Mboup, M.; Leconte, M.; Gautier, A.; Wan, A.M.; Chen, W.; de Vallavieille-Pope, C.; Enjalbert, J. Evidence of genetic recombination in wheat yellow rust populations of a Chinese oversummering area. Fungal Genet. Biol. 2009, 46, 299–307. [Google Scholar] [CrossRef] [PubMed]
  9. Park, R.F.; Wellings, C.R. Somatic hybridization in the uredinales. Annu. Rev. Phytopathol. 2012, 50, 219–239. [Google Scholar] [CrossRef] [PubMed]
  10. Schwessinger, B. Fundamental wheat stripe rust research in the 21st century. New Phytol. 2017, 213, 1625–1631. [Google Scholar] [CrossRef] [PubMed]
  11. Wellings, C.R.; Mclntosh, R.A. Puccinia striiformis f. sp. tritici in Australasia: Pathogenic changes during the first 10 years. Plant Pathol. 1990, 39, 316–325. [Google Scholar] [CrossRef]
  12. Bayles, R.A.; Flath, K.; Hovmoller, M.S.; de Vallavieille-Pope, C. Breakdown of the resistance to yellow rust of wheat in northern Europe. Agronomie 2000, 20, 805–811. [Google Scholar] [CrossRef]
  13. Chen, X.M.; Moore, M.; Milus, E.A.; Long, D.L.; Line, R.F.; Marshall, D.; Jackson, L. Wheat stripe rust epidemics and races of Puccinia striiformis f. sp. tritici in the United States in 2000. Plant Dis. 2002, 86, 39–46. [Google Scholar] [CrossRef] [PubMed]
  14. Ghanbarnia, K.; Gourlie, R.; Amundsen, E.; Aboukhaddour, R. The changing virulence of stripe rust in Canada from 1984 to 2017. Phytopathology 2021, 111, 1840–1850. [Google Scholar] [CrossRef] [PubMed]
  15. Wan, A.M.; Chen, X.M. 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. 2014, 98, 1534–1542. [Google Scholar] [CrossRef] [PubMed]
  16. Boshoff, W.H.P.; Pretorius, Z.A.; van Niekerk, B.D. Establishment, distribution, and pathogenicity of Puccinia striiformis f. sp. tritici in South Africa. Plant Dis. 2002, 86, 485–492. [Google Scholar] [CrossRef] [PubMed]
  17. Milus, E.A.; Seyran, E.; McNew, R. Aggressiveness of Puccinia striiformis f. sp. tritici isolates in the south-central United States. Plant Dis. 2006, 90, 847–852. [Google Scholar] [CrossRef]
  18. Wan, A.M.; Chen, X.M.; Yuen, J. Races of Puccinia striiformis f. sp. tritici in the United States in 2011 and 2012 and comparison with races in 2010. Plant Dis. 2016, 100, 966–975. [Google Scholar] [CrossRef] [PubMed]
  19. Wang, C.H.; Li, Y.X.; Wang, B.T.; Hu, X.P. Genetic analysis reveals relationships among populations of Puccinia striiformis f. sp. tritici from the Longnan, Longdong, and central Shaanxi regions of China. Phytopathology 2022, 112, 278–289. [Google Scholar] [CrossRef] [PubMed]
  20. Cat, A.; Tekin, M.; Akan, K.; Akar, T.; Catal, M. Races of Puccinia striiformis f. sp. tritici identified from the coastal areas of Turkey. Can. J. Plant Pathol. 2021, 43, 323–332. [Google Scholar] [CrossRef]
  21. Wan, A.M.; Chen, X.M.; He, Z.H. Wheat stripe rust in China. Aust. J. Agric. Res. 2007, 58, 605–619. [Google Scholar] [CrossRef]
  22. Chen, W.Q.; Wu, L.R.; Liu, T.G.; Xu, S.C.; Jin, S.L.; Peng, Y.L.; Wang, B.T. Race dynamics, diversity, and virulence evolution in Puccinia striiformis f. sp. tritici, the causal agent of wheat stripe rust in china from 2003 to 2007. Plant Dis. 2009, 93, 1093–1101. [Google Scholar] [CrossRef] [PubMed]
  23. Chen, W.Q.; Wellings, C.; Chen, X.M.; Kang, Z.S.; Liu, T.G. Wheat stripe (yellow) rust caused by Puccinia striiformis f. sp. tritici. Mol. Plant Pathol. 2014, 15, 433–446. [Google Scholar] [CrossRef] [PubMed]
  24. Chen, W.; Zhang, Z.D.; Tian, Y.; Kang, Z.S.; Zhao, J. Race composition and genetic diversity of a Puccinia striiformis f. sp. tritici population from Yunnan and Guizhou epidemiological regions in China in 2018. J. Plant Pathol. 2023, 42, 253–267. [Google Scholar] [CrossRef]
  25. Hu, X.P.; Ma, L.J.; Liu, T.G.; Wang, C.H.; Peng, Y.L.; Pu, Q.; Xu, X.M. Population genetic analysis of Puccinia striiformis f. sp. tritici suggests two distinct populations in Tibet and the other regions of China. Plant Dis. 2017, 101, 288–296. [Google Scholar] [CrossRef] [PubMed]
  26. Yan, D.X.; Huang, L.; Xing, Y.; Li, H.F.; Xia, C.J.; Wang, A.; Gao, L.; Liu, T.G.; Chen, W.Q. Analysis of genetic diversity and population structure of Puccinia striiformis f. sp. tritici infers inoculum relationships from Yunnan to the middle and lower reaches of the Yangtze river. Plant Pathol. 2023, 72, 1517–1527. [Google Scholar] [CrossRef]
  27. Liang, J.M.; Wan, Q.; Luo, Y.; Ma, Z.H. Population genetic structures of Puccinia striiformis in Ningxia and Gansu provinces of China. Plant Dis. 2013, 97, 501–509. [Google Scholar] [CrossRef] [PubMed]
  28. Alam, M.A.; Li, H.X.; Hossain, A.; Li, M.J. Genetic diversity of wheat stripe rust fungus Puccinia striiformis f. sp. tritici in Yunnan, China. Plants 2021, 10, 1735. [Google Scholar] [CrossRef] [PubMed]
  29. Chen, X.M.; Line, R.F. Ientification of stripe rust resistance genes in wheat genotypes used to differentiate North American races of Puccinia striiformis. Phytopathology 1992, 82, 1428–1434. [Google Scholar] [CrossRef]
  30. Chen, X.M.; Line, R.F. Inheritance of stripe rust resistance in wheat cultivars used to differentiate races of Puccinia striiformis in North America. Phytopathology 1992, 82, 633–637. [Google Scholar] [CrossRef]
  31. Hebblethwaite, P.D.; Burbidge, A. The effect of maleic hydrazide and chlorocholine chloride on the growth, seed yield components and seed yield of 23 ryegrass. J. Agric. Sci. 1976, 86, 343–353. [Google Scholar] [CrossRef]
  32. Draz, I.S. Pathotypic and molecular evolution of contemporary population of Puccinia striiformis f. sp. tritici in Egypt during 2016–2018. J. Phytopathol. 2019, 167, 26–34. [Google Scholar] [CrossRef]
  33. Parks, R.; Carbone, I.; Murphy, J.P.; Marshall, D.; Cowger, C. Virulence structure of the eastern U.S. wheat powdery mildew population. Plant Dis. 2008, 92, 1074–1082. [Google Scholar] [CrossRef] [PubMed]
  34. Chen, X.M.; Ronald, F.L.; Leung, H. Relationship between virulence and DNA Polymorphism in Puccinia striiformis. Phytopathology 1993, 83, 1489–1497. [Google Scholar] [CrossRef]
  35. Enjalbert, J.; Duan, X.; Giraud, T.; Vautrin, D.; de Vallavieille-Pope, C.; Solignac, M. Isolation of twelve microsatellite loci, using an enrichment protocol, in the phytopathogenic fungus Puccinia striiformis f. sp. tritici. Mol. Ecol. Notes 2002, 2, 563–565. [Google Scholar] [CrossRef]
  36. Bahri, B.; Leconte, M.; de Vallavieille-Pope, C.; Enjalbert, J. Isolation of ten microsatellite loci in an EST library of the phytopathogenic fungus Puccinia striiformis f. sp. tritici. Conserv. Genet. 2009, 10, 1425–1428. [Google Scholar] [CrossRef]
  37. Chen, C.Q.; Zheng, W.M.; Buchenauer, H.; Huang, L.L.; Lu, N.H.; Kang, Z.S. Isolation of microsatellite loci from expressed sequence tag library of Puccinia striiformis f. sp. tritici. Mol. Ecol. Resour. 2009, 9, 236–238. [Google Scholar] [CrossRef] [PubMed]
  38. Cheng, P.; Chen, X.M.; Xu, L.S.; See, D.R. Development and characterization of expressed sequence tag-derived microsatellite markers for the wheat stripe rust fungus Puccinia striiformis f. sp. tritici. Mol. Ecol. Resour. 2012, 12, 779–781. [Google Scholar] [CrossRef] [PubMed]
  39. Horká, M.; Horký, J.; Matoušková, H.; Šlais, K. Separation of plant pathogens from different hosts and tissues by capillary electromigration techniques. Anal. Chem. 2007, 79, 9539–9546. [Google Scholar] [CrossRef] [PubMed]
  40. Shannon, C.E. A mathematical theory of communication. Bell Syst. Tech. J. 1948, 27, 623–656. [Google Scholar] [CrossRef]
  41. Schachtel, G.A.; Dinoor, A.; Herrmann, A.; Kosman, E. Comprehensive evaluation of virulence and resistance data: A new analysis tool. Plant Dis. 2012, 96, 1060–1063. [Google Scholar] [CrossRef] [PubMed]
  42. Kamvar, Z.N.; Brooks, J.C.; Grünwald, N.J. Novel R tools for analysis of genome-wide population genetic data with emphasis on clonality. Front. Genet. 2015, 6, 208. [Google Scholar] [CrossRef] [PubMed]
  43. Peakall, R.; Smouse, P.E. GenALEx 6.5: Genetic analysis in excel. Population genetic software for teaching and research an update. Bioinformatics 2012, 28, 2537–2539. [Google Scholar] [CrossRef] [PubMed]
  44. Li, C.; Awais, M.; Yang, H.; Shen, Y.Y.; Li, G.K.; Guo, H.F.; Ma, J.B. Races CYR34 and Suwon11-1 of Puccinia striiformis f. sp. tritici played an important role in causing the stripe rust epidemic in winter wheat in Yili, Xinjiang, China. J. Fungi 2023, 9, 436. [Google Scholar]
  45. Zeng, S.M.; Luo, Y. Long-distance spread and interregional epidemics of wheat stripe rust in China. Plant Dis. 2006, 90, 980–988. [Google Scholar] [CrossRef]
  46. Hovmøller, M.S.; Justesen, A.F. Rates of evolution of avirulence phenotypes and DNA markers in a northwest European population of Puccinia striiformis f. sp. tritici. Mol. Ecol. 2007, 16, 4637–4647. [Google Scholar] [CrossRef] [PubMed]
  47. Zhou, A.H.; Wang, J.; Chen, X.M.; Xia, M.H.; Feng, Y.X.; Ji, F.; Huang, L.L. Virulence characterization of Puccinia striiformis f. sp. tritici in China using the Chinese and Yr single-gene differentials. Plant Dis. 2021, 108, 671–683. [Google Scholar]
  48. Hubbard, A.; Lewis, C.M.; Yoshida, K.; Ramirez-Gonzalez, R.H.; de Vallavieille-Pope, C.; Thomas, J.; Kamoun, S.; Bayles, R.; Uauy, C.; Saunders, D.G. Field pathogenomics reveals the emergence of adiverse wheat yellow rust population. Genome Biol. 2015, 16, 23. [Google Scholar] [CrossRef] [PubMed]
  49. Zhan, G.M.; Wang, F.P.; Wan, C.P.; Han, Q.M.; Huang, L.L.; Kang, Z.S.; Chen, X.M. Virulence and molecular diversity of the Puccinia striiformis f. sp. tritici population in Xinjiang in relation to other regions of western China. Plant Dis. 2016, 100, 99–107. [Google Scholar] [CrossRef]
  50. van Silfhout, C.H.; Adriana, G.; Gerechter-amitai, Z.K.; Frida, K. Identification and characterization of resistance to yellow rust and powdery mildew in wild emmer wheat and their transfer to bread wheat. Neth. J. Plant Pathol. 1989, 95, 73–78. [Google Scholar]
  51. Tekin, M.; Cat, A.; Akan, K.; Catal, M.; Akar, T. A new virulent race of wheat stripe rust pathogen (Puccinia striiformis f. sp. tritici) on the resistance gene Yr5 in Turkey. Plant Dis. 2021, 105, 3292. [Google Scholar] [CrossRef]
  52. Kharouf, S.H.; Hamzeh, S.H.; Azmeh, M.F. Races identification of wheat rusts in Syria during the 2019 growing season. Arab J. Plant Prot. 2021, 39, 1–13. [Google Scholar] [CrossRef]
  53. Zhang, G.S.; Zhao, Y.Y.; Kang, Z.S.; Zhao, J. First report of a Puccinia striiformis f. sp. tritici race virulent to wheat stripe rust resistance gene Yr5 in China. Plant Dis. 2019, 104, 284. [Google Scholar] [CrossRef]
  54. Zheng, S.G.; Li, Y.F.; Lu, L.; Liu, Z.; Zhang, C.H.; Ao, D.H.; Li, L.R.; Zhang, C.Y.; Liu, R.; Luo, C.P.; et al. Evaluating the contribution of Yr genes to stripe rust resistance breeding through marker-assisted detection in wheat. Euphytica 2017, 213, 50. [Google Scholar] [CrossRef]
  55. Chen, X.M.; Wang, M.N.; Wan, A.M.; Bai, Q.; Li, M.J.; López, P.F.; Maccaferr, M.; Mastrangelo, A.M.; Barnes, W.C.; Cruz, D.F.C.; et al. Virulence characterization of Puccinia striiformis f. sp. tritici collections from six countries in 2013 to 2020. Can. J. Plant Pathol. 2021, 43, 308–322. [Google Scholar] [CrossRef]
  56. Burdon, J.J.; Silk, J. Sources and patterns of diversity in plant-pathogenic fungi. Phytopathology 1997, 87, 664–669. [Google Scholar] [CrossRef] [PubMed]
  57. Li, G.Q.; Li, Z.F.; Yang, W.Y.; Zhang, Y.; He, Z.H.; Xu, S.C.; Singh, R.P.; Qu, Y.Y.; Xia, X.C. Molecular mapping of stripe rust resistance gene YrCH42 in Chinese wheat cultivar Chuanmai 42 and its allelism with Yr24 and Yr26. Theor. Appl. Genet. 2006, 112, 1434–1440. [Google Scholar] [CrossRef]
  58. Han, D.J.; Wang, Q.L.; Chen, X.M. Emerging Yr26 virulent races of Puccinia striiformis f. sp. tritici are threatening wheat production in the Sichuan Basin, China. Plant Dis. 2015, 99, 754–760. [Google Scholar] [CrossRef]
  59. Wang, C.C.; Zhang, R.; Chu, B.Y.; Wu, B.M.; Ma, Z.H. Population genetic structure of Puccinia striiformis f. sp. tritici at the junction of Gansu, Sichuan and Shaanxi provinces in China. Phytopathol. Res. 2019, 1, 25. [Google Scholar] [CrossRef]
  60. Liu, X.F.; Huang, C.; Sun, Z.Y.; Liang, J.M.; Luo, Y.; Ma, Z.H. Analysis of population structure of Puccinia striiformis in Yunnan province of China by using AFLP. Eur. J. Plant Pathol. 2011, 129, 43–55. [Google Scholar] [CrossRef]
  61. Lu, N.H.; Wang, J.F.; Chen, X.M.; Zhan, G.M.; Chen, C.Q.; Huang, L.L.; Kang, Z.S. Spatial genetic diversity and interregional spread of Puccinia striiformis f. sp. tritici in northwest China. Eur. J. Plant Pathol. 2011, 131, 685–693. [Google Scholar] [CrossRef]
  62. Liang, J.M.; Liu, X.F.; Tsui, C.K.M.; Ma, Z.H.; Luo, Y. Genetic structure and asymmetric migration of wheat stripe rust pathogen in western epidemic areas of China. Phytopathology 2021, 111, 1252–1260. [Google Scholar] [CrossRef] [PubMed]
Genes 15 00542 g001
Figure 1. Distributions of Pst pathotypes according to their virulence complexity.
Figure 1. Distributions of Pst pathotypes according to their virulence complexity.
Genes 15 00542 g001
Genes 15 00542 g002
Figure 2. Frequencies of virulence in both new and old Pst isolates on 19 Chinese differentials across provinces.
Figure 2. Frequencies of virulence in both new and old Pst isolates on 19 Chinese differentials across provinces.
Genes 15 00542 g002
Genes 15 00542 g003
Figure 3. Dendrograms showing virulence similarity among Pst isolates defined by provinces and isolate types, using the 19 Chinese differentials (A) and the 31 single Yr gene lines (B).
Figure 3. Dendrograms showing virulence similarity among Pst isolates defined by provinces and isolate types, using the 19 Chinese differentials (A) and the 31 single Yr gene lines (B).
Genes 15 00542 g003
Genes 15 00542 g004
Figure 4. Principal coordinates analysis (PCoA) showing four Pst populations collected from the Yunnan, Guizhou, Sichuan, and Gansu provinces.
Figure 4. Principal coordinates analysis (PCoA) showing four Pst populations collected from the Yunnan, Guizhou, Sichuan, and Gansu provinces.
Genes 15 00542 g004
Table 1. Distribution and diversity of Pst pathotypes across the four provinces in China.
Table 1. Distribution and diversity of Pst pathotypes across the four provinces in China.
No.Pathotype a YunnanGuizhou Sichuan Gansu No. Isolates/Mean Frequency (%)
1V2, 4, 8, 14(2)1 (2) b1 (1) 75.0
2V2, 4, 8, 12, 14, 16 5(1)64.2
3V2, 8, 141 (1)11 (1) 53.5
5V0(1)(2) (2)53.5
4V4, 14, 16 13 42.8
6V2, 14(2)(1)(1) 42.8
7V2, 8, 14, 19(1)(1)(2) 42.8
8V2, 4, 8, 10, 14, 192(1)(1) 42.8
9V2, 8 (2)(2)42.8
10V2, 4, 14, 16 3 32.1
11V2, 3, 4, 6, 7, 8, 10, 11, 12, 14, 16, 191 1 21.4
12V2, 101 1 21.4
13V3, 7, 10, 19 1(1) 21.4
14V2, 4, 8, 11, 14, 1611 21.4
15V2, 8, 14, 162 21.4
16V2, 8, 19 (2) 21.4
17V2, 7, 8, 14, 19(1) (1)21.4
18V2, 4, 7, 8, 14(1)1 21.4
19V2 (1) (1)21.4
H c3.62.83.53.63.4
a Pathotypes detected more than once are listed with their corresponding virulence (V1–V19) to the 19 Chinese differentials: V0 = avirulent to all differentials; V1 = Trigo Eureka (Yr6); V2 = Fulhard (unknown); V3 = Lutescens 128 (unknown); V4 = Mentana (unknown); V5 = Virgilio (YrVir1 and YrVir2); V6 = Abbondanza (unknown); V7 = Early Premium (unknown); V8 = Funo (YrA,+); V9 = Danish 1 (Yr3); V10 = Jubilejina 2 (YrJu1, YrJu2, YrJu3, and YrJu4); V11 = Fengchan 3 (Yr1); V12 = Lovrin 13 (Yr9,+); V13 = Kangyin 655 (Yr1, YrKy1, and YrKy2); V14 = Suwon 11 (YrSu); V15 = Zhong 4 (unknown); V16 = Lovrin 10 (Yr9); V17 = Hybrid 46 (Yr3b and Yr4b); V18 = T. spelta var album (Yr5); and V19 = Guinong 22 (Yr26). b Old and new isolates are shown inside and outside parentheses, respectively. c Shannon–Wiener pathotype diversity index (H) for the 143 Pst isolates.
Table 2. Virulence frequencies of 143 Pst isolates on 31 Yr genes of wheat lines at the seedling stage.
Table 2. Virulence frequencies of 143 Pst isolates on 31 Yr genes of wheat lines at the seedling stage.
Yr Genes New IsolatesOld Isolatesp-Value
SCNI (29) YNNI (19)GZNI (11)MeanSCOI (22)YNOI (24)GZOI (10)GSOI (28)Mean
Yr2124.157.964.048.736.455.260.057.151.90.400
Yr7631.047.450.0 42.818.220.844.422.226.40.059
Yr5037.942.110.030.00.000.000.000.000.000.008
Yr4072.468.430.056.936.450.040.032.239.70.107
Yr4134.531.630.032.013.64.200.0021.49.800.006
Yr4569.063.291.074.481.891.790.082.186.40.090
YrSu72.457.982.070.863.658.370.075.967.00.314
Yr4362.157.973.064.372.775.070.057.168.70.251
YrA58.652.646.052.431.854.250.060.749.20.351
Yr2558.650.482.063.763.683.390.075.078.00.114
Yr155.247.430.044.231.830.820.046.431.80.265
Yr50.000.000.000.00.000.000.000.000.00
Yr645.268.446.056.518.254.260.053.646.50.230
Yr3c44.857.960.054.231.841.760.046.445.00.150
Yr782.855.691.076.545.558.370.071.461.30.121
Yr827.647.430.035.059.162.570.064.364.00.002
Yr958.652.660.057.127.362.530.057.144.20.146
Yr2420.020.010.016.713.620.020.020.020.00.195
Yr3213.847.420.027.113.64.200.0025.010.70.096
Yr6448.331.682.054.050.041.755.036.045.70.637
Yr1041.410.540.030.64.504.200.0011.04.900.017
Yr150.000.000.000.000.000.000.000.000.00
Yr1744.826.350.042.013.612.520.025.017.80.036
Yr4458.642.150.050.268.287.550.064.367.50.072
YrSp31.021.140.030.713.616.70.0025.013.80.039
YrTr131.031.620.027.50.008.300.000.002.100.001
Yr2617.216.820.018.04.508.300.0011.16.000.005
Yr2979.355.073.069.155.050.030.046.445.40.022
YrJu483.368.470.073.986.487.555.085.777.40.387
YrKy230.052.660.047.554.554.250.064.355.80.185
YrRes93.378.970.080.768.258.370.082.169.70.115
Table 3. Unique virulence patterns (UVPs), average virulence complexity (AVC), relative virulence complexity (RVC), and virulence uniformity (VU) among populations of Pst in the western provinces of China.
Table 3. Unique virulence patterns (UVPs), average virulence complexity (AVC), relative virulence complexity (RVC), and virulence uniformity (VU) among populations of Pst in the western provinces of China.
Population Total Isolates UVP (%)AVCRVC VU
SCNI 2910014.4 0.460.31
SCOI2210010.60.340.36
YNNI1910013.70.440.30
YNOI2410012.40.400.26
GZNI1110014.40.460.23
GZOI1010011.50.370.27
GSOI2810013.60.440.24
Table 4. Pairwise associations between virulence and avirulence in Pst isolates, at respective Yr loci.
Table 4. Pairwise associations between virulence and avirulence in Pst isolates, at respective Yr loci.
Yr Gene PairPathogen Alleles aTotal VV b ProportionAssociation Type c p-Value d
AAAVVAVV
Yr76 41881011341430.24+<0.001
Yr76 329363591430.06+0.016
Yr76 1092734101430.07+0.012
Yr76 Sp861430131430.09+0.038
Yr76 Tr191834101430.07+0.026
Yr50 4111917341430.03+0.009
Yr50 51302741430.03+<0.001
Yr50 2411026341430.03+0.035
Yr50 1012511341430.03+0.002
Yr50 Sp11422341430.03+0.021
Yr50 Tr112511341430.03+0.002
Yr50 2612115431430.02+0.043
Yr50 Ju428060551430.38+/−<0.001
Yr50 Ky26967071430.05+/−0.014
Yr41 1011481471430.05+0.002
Yr41 1798247141430.10+<0.001
Yr41 Ky264585161430.11+/−0.018
Yr24 32101273121430.08+<0.001
Yr24 1010582371430.05+0.017
Yr24 26105820101430.07+<0.001
Yr32 2611711871430.05+<0.001
Yr17 Sp921325131430.09+0.006
Yr17 2698727111430.08+0.001
YrSp Tr110981971430.05+0.007
YrSp 26107101881430.06+0.005
YrTr1 2611711871430.05+<0.001
a Number of virulent isolates (V) or avirulent (A), at respective Yr genes. b Proportion of isolates virulent to both Yr genes. c +, a positive association with AA or VV; +/−, an unclear association. d Pairs significantly different at p < 0.05 based on Fisher’s exact test are listed.
Table 5. Population genetic parameters among Pst populations.
Table 5. Population genetic parameters among Pst populations.
Pst Population NNaNeIuhPPL (%)
Old population841.941.580.510.3497
New population591.971.530.480.3297
Guizhou population 211.811.540.460.3286
Sichuan population 512.001.540.490.33100
Yunnan population 432.001.560.510.34100
Gansu population 281.941.560.510.3597
Mean 47.671.941.550.490.3396
N, number of samples; Na, number of polymorphic alleles; Ne, number of effective alleles; I, Shannon information index; uh, unbiased diversity; PPL, percentage of polymorphic loci.
Table 6. Analysis of molecular variance (AMOVA) among provinces (Sichuan, Yunnan, Gansu, and Guizhou) and isolate types (old and new isolates) of Pst populations.
Table 6. Analysis of molecular variance (AMOVA) among provinces (Sichuan, Yunnan, Gansu, and Guizhou) and isolate types (old and new isolates) of Pst populations.
Source of VariationDegree of FreedomSum of SquaresComponent of VariationsPercentage of
Variation (%)		p-Value
Among provinces356.60.460.001
Within provinces139839.66.0940.001
Total142896.36.4100
Between isolate groups 157.10.7110.001
Within isolate groups141848.06.0890.001
Total142905.16.8100
Table 7. Pairwise matrix showing the Nei’s genetic distance among provinces of the Pst populations.
Table 7. Pairwise matrix showing the Nei’s genetic distance among provinces of the Pst populations.
YNGZ SCGS
YN
GZ0.06
SC0.030.06
GS 0.040.060.10
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Gebremariam, T.G.; Wang, F.; Lin, R.; Li, H. Comparative Analysis of Virulence and Molecular Diversity of Puccinia striiformis f. sp. tritici Isolates Collected in 2016 and 2023 in the Western Region of China. Genes 2024, 15, 542. https://doi.org/10.3390/genes15050542

AMA Style

Gebremariam TG, Wang F, Lin R, Li H. Comparative Analysis of Virulence and Molecular Diversity of Puccinia striiformis f. sp. tritici Isolates Collected in 2016 and 2023 in the Western Region of China. Genes. 2024; 15(5):542. https://doi.org/10.3390/genes15050542

Chicago/Turabian Style

Gebremariam, Tesfay Gebrekirstos, Fengtao Wang, Ruiming Lin, and Hongjie Li. 2024. "Comparative Analysis of Virulence and Molecular Diversity of Puccinia striiformis f. sp. tritici Isolates Collected in 2016 and 2023 in the Western Region of China" Genes 15, no. 5: 542. https://doi.org/10.3390/genes15050542

APA Style

Gebremariam, T. G., Wang, F., Lin, R., & Li, H. (2024). Comparative Analysis of Virulence and Molecular Diversity of Puccinia striiformis f. sp. tritici Isolates Collected in 2016 and 2023 in the Western Region of China. Genes, 15(5), 542. https://doi.org/10.3390/genes15050542

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop