Next Article in Journal
Aquatic Plants in Ponds at the Brdo Estate (Slovenia) Show Changes in 20 Years
Next Article in Special Issue
YOLOv5s-Based Image Identification of Stripe Rust and Leaf Rust on Wheat at Different Growth Stages
Previous Article in Journal
Landscape, Human Disturbance, and Climate Factors Drive the Species Richness of Alien Invasive Plants on Subtropical Islands
Previous Article in Special Issue
Stem Rust Resistance and Resistance-Associated Genes in 64 Wheat Cultivars from Southern Huanghuai, China
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

The Migration, Diversity, and Evolution of Puccinia triticina in China

1
Biological Control Center of Plant Diseases and Plant Pests of Hebei Province, College of Plant Protection, Hebei Agricultural University, Baoding 071000, China
2
School of Landscape and Ecological Engineering, Hebei Engineering University, Handan 056038, China
*
Authors to whom correspondence should be addressed.
The authors contributed equally to this work.
Plants 2024, 13(17), 2438; https://doi.org/10.3390/plants13172438
Submission received: 17 June 2024 / Revised: 21 August 2024 / Accepted: 21 August 2024 / Published: 31 August 2024
(This article belongs to the Special Issue Plant Pathology and Epidemiology for Grain, Pulses, and Cereal Crops)

Abstract

:
Wheat leaf rust, caused by Puccinia triticina, is one of the most common fungal diseases of wheat in China and occurs widely in various wheat-growing regions. To clarify the epidemic, spread rules, and population structure of P. triticina among different regions, 217 isolates of P. triticina collected from Hebei, Shandong, Sichuan, and Xinjiang in China were tested by 34 Thatcher near-isogenic lines and 21 pairs of EST-SSR primers. A total of 83 races were identified, and THTT, PHTT, THTS, and PHJT were the most predominant races in the four provinces in 2009. We found enriched virulence and genetic diversity in the four P. triticina populations and a significant correlation between genetic polymorphism and geographic regions. However, no significant correlation was found between virulence phenotypes and molecular genotypes. Moreover, a notable high level of gene flow (Nm = 2.82 > 1) among four P. triticina populations was detected. The genetic relationship among Hebei, Shandong, and Sichuan populations was close, possibly due to the spread of P. triticina from Sichuan to Shandong and then to Hebei. In contrast, the Xinjiang population was relatively independent. Genetic differentiation analysis showed some level of differentiation among or within populations of P. triticina in the four provinces, and the genetic variation within populations (74.97%) was higher than across populations (25.03%). Our study provides a basis for a better understanding of the regional migration, epidemic, and population structure of P. triticina in China.

1. Introduction

Puccinia triticina Eriks., the wheat leaf rust fungus of common wheat (Triticum aestivum L.), is prevalent in wheat-growing regions worldwide [1]. Wheat leaf rust is well adapted to the climate of wheat growth, so it is more widely distributed than stem rust or stripe rust [2,3,4,5]. The occurrence and prevalence of leaf rust under favorable climatic conditions can cause significant yield losses in susceptible wheat cultivars [4]. Recent studies predicted that the annual average global food loss attributable to wheat leaf rust will be 8.6 million metric tons from 2000 to 2050 based on a conservative scenario and 18.3 million metric tons based on a high-loss scenario [6]. China is the largest wheat producer in the world, with over 23.63 million hectares of wheat-cultivation area and 136.59 million tons of yield in 2023 (National Bureau of Statistics of China). However, China is also one of the regions most severely impacted by wheat leaf rust, with notable epidemics occurring in 1969, 1973, 1975, 1979, 2012, 2013, 2015, and 2020, leading to serious yield loss in the major wheat-growing regions [7,8]. According to statistics, wheat leaf rust affects approximately 15 million hectares annually in China, resulting in a yield loss of approximately 3 million tons. Resistant cultivars, fungicides, or an improvement of disease management practices are the most commonly used methods to control the disease [6], among which the promotion and application of resistant cultivars is the most economical, effective, and environmentally friendly [9]. However, due to the high variation rate of P. triticina, new races constantly evolve and eventually overcome the resistance of wheat cultivars, threatening wheat production [6]. Therefore, studying the virulence dynamics and population genetic structure of P. triticina will significantly contribute to the control of the occurrence and development of the disease and also guide wheat disease resistance breeding programs.
Since the 1920s, race surveys and virulence identification of P. triticina have been carried out to monitor the virulence dynamics of P. triticina populations. The annual race and virulence survey of P. triticina in the United States began in 1926 and in Canada since 1931 [10,11,12,13]. The comprehensive and systematic analyses of the race and virulence dynamics of P. triticina were conducted in the United States from 2000 to 2020 [14] and in Canada from 1997 to 2010 [15,16,17]. Moreover, similar studies have been reported in France [18], South America [19], Ukraine [20], India [21,22], Egypt [23], and Iran [24]. In China, studies on the virulence dynamics of P. triticina races began in the 1940s [25]. Wheat-growing regions in China are mainly divided into North China, Huang-Huai-Hai, the middle and lower reaches of the Yangtze River, and southwest and northwest regions. In the past, wheat leaf rust mainly occurred in the Southwest China wheat-growing region. However, climate warming has triggered the spread of leaf rust to the north of China. In the 1990s, leaf rust was more severe in all the wheat-growing regions across the country, especially in the Huang-Huai-Hai wheat-growing region. Since 1986, regional race and virulence studies of P. triticina in China have been carried out based on the Thatcher near-isogenic lines [7,8,26,27,28,29,30,31]. These studies on population dynamics of P. triticina based on virulence are of great significance in studying the virulence variation and diversity of P. triticina. Unlike P. striiformis f. sp. tritici, the regional epidemics and spread patterns of P. triticina in China are currently unclear. P. striiformis mainly spread from southern Gansu and northwest Sichuan to other wheat-growing regions, and P. striiformis in the northwest of China, especially in Xinjiang, is different from other regions in terms of population structure and spread [32,33]. However, studies on the regional epidemics and spread rules of P. triticina are still rare and unsystematic, so it is necessary to strengthen the study of Chinese P. triticina in migration, spread, and population structure.
Molecular marker technologies have been widely applied to evaluate the genetic diversity of P. triticina populations in the past 20 years [34]. The development of these technologies enables us to understand the migration, evolution, diversity, and population structure of P. triticina at the molecular level. Random amplified polymorphic DNA (RAPD) [35,36,37], amplified fragment-length polymorphism (AFLP) [38], simple sequence repeat (SSR) [39], and expressed sequence tag-simple sequence repeat (EST-SSR) [40,41,42,43] have been successively developed and applied to study the population structure and genetic diversity of P. triticina. A previous study found that the P. triticina population of China was strongly differentiated from other regions of the world, which may indicate that the Chinese P. triticina population may be endemic in its own region [44]. However, there are few studies on the genetic diversity, migration, spread, and epidemic routes of P. triticina populations among different regions of China. Therefore, in order to clarify the relations of the population of different regions and the rules of the epidemic in different regions, this study analyzed the genetic and virulence diversity of P. triticina populations in Heibei Province, Shandong Province, Sichuan Province, and the Xinjiang Uygur Autonomous Region of China, which can provide further evidence for the regional spread and migration of P. triticina in China.

2. Results

2.1. Dynamics of Races

In total, 217 single-uredinial isolates were selected from Heibei (48 isolates, No. 1–48), Shandong (48 isolates, No. 49–96), Sichuan (44 isolates, No. 97–140), and Xinjiang (77 isolates, No. 141–217) in 2009 (Table 1). These isolates were successfully identified in 83 races (Table 1 and Table 2). THTT, PHTT, THTS, and PHJT were the predominant races with 23.50%, 9.22%, 8.29%, and 5.53% frequencies, respectively (Table 3). However, there were still slight differences in the types of predominant races among different regions (Table 2 and Table 3).

2.2. Virulence Polymorphism Analysis

The similarity coefficient of virulence between 217 isolates was 0.65–1.00 by constructing dendrograms of virulence based on the UPGMA method (Figure 1). All P. triticina isolates were clustered into nine groups (V1–V9) when the similarity coefficient was 0.73. A total of 187 isolates (86.18%) were grouped in cluster V1, and 30 isolates were spread across the other groups (V2–V9). The results showed that these isolates were not clustered into different groups according to their source regions, which indicated that the virulence polymorphism of P. triticinia populations was not closely correlated with geographical origin, although a few isolates with the same geographical origin were clustered in the same group.

2.3. Genetic Polymorphism Analysis

Two hundred and seventeen isolates of P. triticina in 2009 were analyzed by 21 pairs of EST-SSR primers. Only 18 pairs of primers could amplify clear and stable polymorphic bands, and a total of 48 alleles were amplified, of which 30 alleles were polymorphism alleles, and the percentage of polymorphism was 62.50%. The number of alleles and polymorphic loci of these 18 pairs of primers was 2–4 and 1–3, respectively (Table 4).
The genetic similarity coefficient based on an EST-SSR analysis of 217 isolates was 0.62–1.00 (Figure 2). Most isolates were partitioned into two main clusters (S1 and S2) at a similarity coefficient of 0.68. Cluster S1 contained most of the isolates from Hebei, Shandong, and Sichuan. These isolates were clustered into different subclusters according to the source region, such as subclusters S1-1, S1-3, and S1-5 which mainly contained the isolates of Shandong, subcluster S1-2 which mainly contained the isolates of Hebei, and subcluster S1-4 which mainly contained the isolates of Sichuan. Cluster S2 contained three subclusters, S2-1, S2-2, and S2-3, and most of the isolates of these subclusters were mainly from Xinjiang. These results showed that most of the isolates originating from the same region were clustered in the same group or several neighboring subclusters with a high level of similarity, indicating that the genetic structure in these isolates may be related to their source region. Moreover, the clustering results also showed that the isolates from Hebei, Shandong, and Sichuan were closely related, while the isolates from Xinjiang were far from the other three regions.

2.4. Correlation Analysis between Virulence Polymorphism and EST-SSR Polymorphism

By comparing the virulence polymorphism and EST-SSR polymorphism cluster analysis dendrogram of P. triticina, it was found that there was no obvious correlation between the genetic cluster based on EST-SSR and virulence characteristics of P. triticina isolates. The isolates with different phenotypes may have the same or similar DNA fingerprints and have close relationships with each other. Conversely, isolates with the same phenotype may have different DNA fingerprints and distant genetic relationships. In addition, the correlation coefficient of virulence and EST-SSR polymorphisms was 0.08 (<0.3) by using the MXCOMP module of NTSYS-pc, which also indicated that there was no correlation between virulence and EST-SSR polymorphism.

2.5. Two-Dimensional Principal Coordinate Analysis of EST-SSR and Virulence

In order to further determine the correlation between virulence polymorphism or EST-SSR polymorphism and geographical region, all isolates were analyzed by a 2D-PCA plot analysis. In the 2D-PCA plot analysis of virulence, all isolates from four provinces were staggered without obvious regional distribution (Figure 3). In the 2D-PCA plot analysis of EST-SSR, all isolates were also mainly distributed in three areas (A–C) (Figure 4). The isolates from Heibei and Shandong were mainly concentrated in area A, which indicated that these isolates have similar genetic loci. In addition, a few isolates of Hebei were also distributed in area C. The isolates from Sichuan were mainly concentrated in area B, and isolates from Xinjiang were mainly concentrated in area C. These results showed that most of the isolates were obviously clustered according to the geographical region. The population in Xinjiang was relatively independent, while those in the other three regions were relatively close. This further confirmed that the molecular polymorphism based on the EST-SSR analysis was more closely related to the geographical region than the virulence phenotype polymorphism.

2.6. Virulence and EST-SSR Polymorphism Analysis of Predominant Races

In order to further determine whether the pathotypes were related to virulence polymorphism or EST-SSR polymorphism of P. triticina, a cluster analysis of 101 isolates of four predominant races (THTT, PHTT, THTS, and PHJT) was carried out based on virulence or EST-SSR. In the polymorphism analysis based on virulence, most of the same races tended to cluster in the same or adjacent cluster or subcluster (Figure 5), indicating a certain correlation between the pathotypes of the isolates and their virulence polymorphisms. In contrast, the clustering results of most races based on the EST-SSR polymorphism have no obvious rules (Figure 6). The isolates of different races may have the same or similar DNA fingerprints and be grouped into a cluster, while those of the same race with varying fingerprints of DNA can also be grouped into different clusters. The results showed no significant correlation between the phenotype and EST-SSR polymorphism.

2.7. Genetic Diversity Analysis of Puccinia triticina Population

All P. triticina isolates have been classified into four distinct populations based on their regions. Genetic diversity analysis among different populations was further conducted using the POPGEN version 1.32. The percentages of polymorphic loci were 83.67 in Hebei, 67.35 in Shandong, 71.43 in Sichuan, and 81.63 in Xinjiang (Table 5). Nei’s gene diversity index (H) and Shannon’s information index (I) both reflect the degree of variation among different populations with similar trends and the higher the value of H and I, the richer the corresponding gene variations [45]. Our results showed that the Nei’s gene diversity index (H) of each population within the P. triticina populations ranged from 0.28 to 0.36, and Shannon’s information index (I) ranged from 0.40 to 0.52, while the highest and lowest levels of genetic diversity were found in Heibei and Shandong, respectively (Table 5).
In addition, the analysis results of POPGEN software showed that the total genetic diversity (Ht) of all isolates, the genetic diversity within the population (Hs), and the genetic diversity among populations (Dst) were 0.35, 0.30, and 0.05, respectively. Moreover, the coefficient of genetic differentiation (Gst) was 0.15, of which genetic diversity within populations accounted for 85.71% of the total genetic diversity, and genetic diversity among populations accounted for 14.29%. There was also genetic differentiation among or within P. triticina populations per the AMOVA method of Arlequin ver. 3.11. Genetic variation within populations accounted for 74.97% of the total variation, while the genetic variation among the population accounted for 25.03% (Table 6), indicating that variations are mainly found within populations. Moreover, a relatively high level of gene flow (Nm = 2.82 > 1) was detected among the four populations. These results indicated a certain degree of exchange of inoculum sources across populations of different regions, which may play a critical role in keeping the similarity of populations’ genetic structure.
To further clarify the genetic relationship between different P. triticina populations from four regions, the populations’ genetic distance and genetic identities were estimated using the POPGENE software (Table 7), and the results showed that the population in Hebei had the closest genetic distance and the highest genetic identity to that in Shandong.
A dendrogram (UPGMA, unweighted average method) was then constructed by MEGA 7.0 software based on Nei’s genetic distances (Figure 7). The four populations could be clustered into two main branches (I–II) (Figure 7). The genetic relationship among the three populations from Hebei, Shandong, and Sichuan was relatively close, especially Hebei and Shandong, while the populations from Xinjiang had the farthest relationship with other populations. This indicated that the genetic relationship among populations of P. triticina had a certain correlation with geographical distribution.

3. Discussion

Wheat leaf rust is one of the most widespread fungal diseases in China and occurs in all wheat-growing regions across the country. Due to the lack of studies on the spread and migration of P. triticina in China, the regional epidemic of P. triticina is relatively unclear, especially compared with that of P. striiformis f. sp. tritici. So, to clarify the migration route, population diversity, and genetic variation of P. triticina in China, the method of virulence identification combined with EST-SSR molecular marker technology in this study was used for 217 P. triticina isolates from Hebei, Shandong, Sichuan, and Xinjiang in China. The results showed that the analysis methods can reveal the polymorphisms and differences in the genetic structure of P. triticina populations to a certain extent.
In this study, a relatively rich genetic diversity was found among the P. triticina isolates from Hebei, Shandong, Sichuan, and Xinjiang in China, as well as a difference in the genetic structure of P. triticina populations in different regions. In general, the main factors affecting the genetic structure of P. triticina populations are the mutation of the pathogen, the selection effects of host cultivars, and the climate environment. Mutation is likely a recurrent event in P. triticina populations because new virulence phenotypes are often detected shortly after the release of wheat cultivars with race-specific Lr genes [5,46]. So, some level of genetic diversity was maintained mainly by mutation in the highly clonal populations of P. triticina [47]. In the field populations of P. triticina, ultraviolet radiation accounts for most mutations [48]. According to the altitude, latitude, and climate of different regions in China, among the four provinces of Hebei, Shandong, Sichuan, and Xinjiang, the ultraviolet intensity of Sichuan is the weakest, while Xinjiang is the highest. Therefore, the mutation probability may differ in these regions, and the P. triticina isolates are regional specificity in Xinjiang. In addition, the genetic diversity of P. triticina can be affected by wheat cultivars and regional environments [49]. The four provinces of Hebei, Shandong, Sichuan, and Xinjiang belong to three different wheat-growing regions. Hebei and Shandong belong to the North China Plain wheat-growing region, Sichuan belongs to the southwest wheat-growing region, and Xinjiang belongs to the northwest wheat-growing region, and the wheat cultivars planted in these regions are quite different. For example, the main wheat cultivars in Sichuan are ‘Chuanmai’ series cultivars, ‘Xindong’ and ‘Xinchun’ series cultivars are the main wheat cultivars in Xinjiang, the main wheat cultivars in Shandong are ‘Jimai’, ‘Lumai’ and ‘Yanmai’ series cultivars, and ‘Jimai’, ‘Shimai’, ‘Hanmai’, and ‘Liangxing 99’ lines are the main wheat cultivars in Hebei. Therefore, the difference in wheat cultivars in these regions may result in different selection pressures on P. triticina populations, which may cause the virulence or dominant races of P. triticina to evolve in different directions and may also easily lead to the diversity of P. triticina races, resulting in the changes in population structure. Moreover, precipitations, temperatures, and altitude can affect the development of P. triticina populations [50]. The three wheat-growing regions (including the four provinces) have different geographical and climatic conditions. Shandong and Hebei belong to the temperate monsoon climate, Sichuan belongs to the subtropical monsoon climate and plateau mountain climate, and Xinjiang belongs to the temperate continental climate. This significant climatic difference in these three wheat-growing regions may affect the fitness of P. triticina, the size of P. triticina populations, and the composition of dominant races, and may then affect the P. triticina population structure.
Migration was also an important factor in the genetic structure of P. triticina populations [51,52,53]. The spores of P. triticina can spread long distances with airflow, which results in gene flow among populations in different regions, and gene flow could change the gene frequency of the receiving population. We found a relatively high level of gene flow (Nm = 2.82 > 1) among the four populations by using the AMOVA method of Arlequin3.ll software. Nm (gene flow) > 1 means gene flow between populations, which could prevent genetic differentiation between populations [54]. In this study, we found a significant correlation between genetic polymorphism and the geographical origin of four populations, and the populations of Hebei, Shandong, and Sichuan had a closer genetic relationship, while the Xinjiang population was a relatively independent population. This may be related to the geographical environment of the four provinces. Hebei and Shandong are mainly affected by the southeast monsoon from the Pacific Ocean, and the two provinces are adjacent to each other, so P. triticina isolates can easily spread to Hebei from Shandong. The southwest monsoon in the southwest region of China, emanating from the Indian Ocean, affects Yunnan, Chongqing, Sichuan, and even Southeast China. So, P. triticina from India may migrate to China along with the monsoon, which may be the reason that there are similar predominant races in China and India; THTT and PHTT are the predominant races in both regions [22]. Moreover, according to the occurrence and development time of leaf rust in the four provinces, the peak occurrence time of leaf rust in Sichuan is from March to the end of April, in Shandong from mid May to early June, and in Hebei, it is a few days later than that in Shandong. Therefore, the P. triticina isolates in the southwest wheat-growing region are likely to spread towards southeast China with the airflow of southwest wind in May, then spread to Shandong and then Hebei by the southeast wind, resulting in gene flow among the three regions. Compared with the other three provinces, the Xinjiang region is located in the northwest of China and has a relatively high altitude. This region is affected by the west monsoon from the Atlantic Ocean. Hence, P. triticina may spread from west to east along the Tianshan Mountains in Xinjiang. Moreover, leaf rust in the western region of Xinjiang occurs earlier and is more severe than that in the eastern region. However, the wheat-growing regions in Xinjiang are geographically isolated, separated at a long distance from other wheat-growing regions of China by deserts and mountains [33], so the migration of P. triticina isolates is more difficult and most isolates may only spread within Xinjiang like P. striiformis [55], although there was some gene flow between Xinjiang and the other three provinces. Therefore, the above factors led to higher genetic similarity among the P. triticina populations of Hebei, Shandong, and Sichuan, while the P. triticina populations of Xinjiang were relatively independent.
This study found no significant correlation between virulence phenotypes and genotypes of the four P. triticina populations. This was consistent with the results of other studies on the population diversity of P. triticina in China [56,57]. For example, our previous study found no significant correlation between virulence and molecular polymorphism of P. triticina populations in Hebei Province of China from 2001 to 2010 [58]. Similar results were reported in Western Europe, South Asia, Russia, and Kazakhstan [35,59,60]. In addition, the same conclusion had been reported in studies on the population genetics of P. striiformis f. sp. tritici in the United States and Northwestern China [33,61]. However, there were also some different views on the correlation between virulence phenotypes and molecular genotypes of P. triticina. Some studies on genetic diversity have reported a significant correlation between virulence phenotypes and molecular genotypes of P. triticina in Pakistan, Europe, Central Asia, the Caucasus regions, North America, South America, and other worldwide regions [31,39,44,52,62,63,64,65,66,67,68,69], which was controversial with the views of the above studies.
The polymorphism analysis of 101 P. triticina isolates of four predominant pathotypes showed a significant correlation between pathotype and virulence polymorphism, not between phenotype and EST-SSR polymorphism. The majority of the same pathotypes were grouped into one cluster or the same subcluster based on virulence, although the clustering of individual pathotypes had a certain correlation with the regional origin. In the EST-SSR polymorphism analysis, different pathotypes were clustered into one cluster or subcluster, but these pathotypes had a certain correlation with their regional origin. These results indicated that the same pathotypes did not completely cluster into one group. In contrast, the different pathotypes may be clustered, revealing that these isolates were still different or identical although they were of the same or different pathotypes, which also showed that the virulence phenotype was not correlated with the molecular genotype. Similar results have also been reported by Wang et al. [42], who noted that TDBG and TDBJ, which differed in virulence on Lr14a, were clustered in the same cluster. The similarity of virulence might be due to the selection effect of wheat cultivars with different leaf rust resistance genes. Zhao [70] and Wang [71] thought that there was virulence heterogeneity in the same pathotype of P. triticina, and the pathogenicity of different isolates of the same pathotype of P. triticina was not completely consistent.

4. Materials and Methods

4.1. Leaf Rust Sample Collection and Multiplication

Leaf rust samples were collected from Heibei Province, Shandong Province, Sichuan Province, and the Xinjiang Uygur Autonomous Region in 2009 in China by researchers from the College of Plant Protection, Hebei Agricultural University (Table 1). All samples were from fields under natural infection. Wheat leaves with uredinia of P. triticina collected from a single plant or cultivar were treated as a single sample. The sample treatment and single-uredinial isolates multiplication were performed as described by Zhang et al. [7]. Urediniospores collected from each sample were inoculated to the susceptible cultivar ‘Zhengzhou 5389’ to obtain the single uredinium. Urediniospores from the single uredinium were increased by inoculating the new seedlings of ‘Zhengzhou 5389’ using the same inoculation procedure. Twelve days after inoculation, urediniospores were collected by shaking the leaf and putting the spores into glass tubes to dry, after which they were lyophilized and stored at −20 °C.

4.2. Race and Virulence Phenotype Data

Race identification and virulence determination of isolates was performed as described by Zhang et al. [7]. Urediniospores collected from each sample were inoculated on two sets of Thatcher near-isogenic lines, which including 16 standardized differentials (with single resistance genes Lr1, Lr2a, Lr2c, Lr3; Lr9, Lr16, Lr24, Lr26; Lr3Ka, Lr11, Lr17, Lr30; LrB, Lr10, Lr14a, and Lr18) and 18 extra differentials (with resistance genes Lr2b, Lr3bg, Lr14b, Lr15, Lr19, Lr21, Lr23, Lr28, Lr33, Lr33 + 34, Lr36, Lr37, Lr38, Lr39, Lr42, Lr44, Lr47, and Lr50). The susceptible Thatcher was used as a control. When uredinia were fully developed on the leaves of Thatcher about 12 d post-inoculation (dpi), infection types (IT) on the near-isogenic lines were recorded based on a 0 to 4 scale [72]. Races were determined according to Kolmer [38].

4.3. DNA Extraction and PCR Amplification

The genomic DNA of each isolate sample was extracted from the urediniospores according to the modified CTAB method described by Gao [73]. The 21 pairs of EST-SSR primers developed by Wang were used for PCR reaction [43], and the total reaction system and amplification procedure were performed as described by Zhang et al. [40]. PCR products were separated by 10% (w/v) polyacrylamide gel in a 0.5 × TBE buffer (45 mM Tris, 45 mM Boric acid, and 1 mM EDTA) and were visualized by silver staining.

4.4. Data Analysis

To analyze the genetic relationships among the selected P. triticina isolates, two binary matrices based on virulence infection types (IT) and EST-SSR data were constructed, respectively, as described by Zhang et al. [40]. The software program NTSYS-pc version 2.10 [74] was used to construct the dendrograms, the 2D principal coordinate analysis, and the correlation analysis between the virulence and EST-SSR. The software POPGENE version 1.32 [75] was used to assess the population diversity of P. triticina. To evaluate population genetic differentiations, an analysis of molecular variance (AMOVA) was performed using the AMOVA model of the Arlequin ver. 3.11 software.

Author Contributions

H.Y. and Q.M. conceived and designed the research. L.Z. and P.Z. conducted the experiments and analyzed the data. Q.M. participated and assisted in the virulence identification of isolates. L.Z. and H.Y. wrote the manuscript. H.Y. and D.L. reviewed and edited the manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

This research was supported by the Hebei Province Dry Alkali Wheat Industry Technology System (grant number HBCT2024030206); the Natural Science Foundation of Hebei Province (grant number C2024204089, C2024402036); and the Key Research and Development Project of Hebei Province (grant number 21326508D).

Data Availability Statement

The data presented in this study are available upon request from the corresponding author.

Acknowledgments

We thank Vincent Ninkuu from Henan University, Kaifeng, Henan, China, for his assistance in revising this article. We also thank the anonymous reviewers for their constructive comments on this article.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Bolton, M.D.; Kolmer, J.A.; Garvin, D.F. Wheat leaf rust caused by Puccinia triticina. Mol. Plant Pathol. 2008, 9, 563–575. [Google Scholar] [CrossRef] [PubMed]
  2. Roelfs, A.P.; Singh, R.P.; Saari, E.E. Rust Diseases of Wheat: Concepts and Methods of Disease Management; CIMMYT: México, Mexico, 1992. [Google Scholar]
  3. Chester, K.S. The nature and prevention of cereal rust as exemplified in the leaf rust of wheat. Nature 1948, 161, 220–221. [Google Scholar]
  4. Li, H.; Hua, L.; Zhao, S.; Hao, M.; Song, R.; Pang, S.; Liu, Y.; Chen, H.; Zhang, W.; Shen, T.; et al. Cloning of the wheat leaf rust resistance gene Lr47 introgressed from Aegilops speltoides. Nat. Commun. 2023, 14, 6072. [Google Scholar] [CrossRef]
  5. Kolmer, J.; Chen, X.; Jin, Y. Diseases which challenge global wheat production-the wheat rusts. In Wheat Science and Trade; Carver, B.F., Ed.; Wiley-Blackwell: Oxford, UK, 2009; pp. 89–124. [Google Scholar]
  6. Chai, Y.; Pardey, P.G.; Hurley, T.M.; Senay, S.D.; Beddow, J.M. A probabilistic bio-economic assessment of the global consequences of wheat leaf rust. Phytopathology 2020, 110, 1886–1896. [Google Scholar] [CrossRef]
  7. Zhang, L.; Shi, C.; Li, L.; Li, M.; Meng, Q.; Yan, H.; Liu, D. Race and virulence analysis of Puccinia triticina in China in 2014 and 2015. Plant Dis. 2020, 104, 455–464. [Google Scholar] [CrossRef]
  8. Zhang, L.; Xiao, Y.; Gao, Y.; Zhao, N.; An, Y.; Yang, W.; Meng, Q.; Yan, H.; Liu, D. Race and virulence analysis of Puccinia triticina in China during 2011 to 2013. Plant Dis. 2020, 104, 2095–2101. [Google Scholar] [CrossRef]
  9. Pink, D.A.C. Strategies using genes for non-durable disease resistance. Euphytica 2002, 124, 227–236. [Google Scholar] [CrossRef]
  10. Johnson, T. Physiologic races of leaf rust of wheat in Canada 1931 to 1955. Can. J. Agric. Sci. 1956, 36, 371–379. [Google Scholar]
  11. Johnston, C.O.; Caldwell, R.M.; Compton, L.E.; Browder, L.E. Physiologic races of Puccinia recondita f. sp. tritici in the United States from 1926 through 1960; USDA Technical Bulletin; USDA: Washington, DC, USA, 1968; Volume 1393, pp. 1–18. [Google Scholar]
  12. Kolmer, J.A. Phenotypic diversity in two populations of Puccinia recondita f.sp. tritici in Canada during 1931–1987. Phytopathology 1991, 81, 311–315. [Google Scholar] [CrossRef]
  13. Mains, E.B.; Jackson, H.S. Physiologic specialisation in the leaf rust of wheat, Puccinia triticina Erikas. Phytopathology 1926, 16, 89–120. [Google Scholar]
  14. Kolmer, J.A.; Fajolu, O. Virulence phenotypes of the wheat leaf rust pathogen, Puccinia triticina, in the United States from 2018 to 2020. Plant Dis. 2022, 106, 1723–1729. [Google Scholar] [CrossRef]
  15. Kolmer, J.A. Physiologic Specialization of Puccinia triticina in Canada in 1997. Plant Dis. 1999, 83, 194–197. [Google Scholar] [CrossRef]
  16. Kolmer, J.A. Physiologic specialization of Puccinia triticina in Canada in 1998. Plant Dis. 2001, 85, 155–158. [Google Scholar] [CrossRef] [PubMed]
  17. MCcallum, B.; Seto-Goh, P.; Xue, A. Physiologic specialization of Puccinia triticina, the causal agent of wheat leaf rust, in Canada in 2010. Can. J. Plant Pathol. 2016, 38, 440–447. [Google Scholar] [CrossRef]
  18. Goyeau, H.; Park, R.; Schaeffer, B.; Lannou, C. Distribution of pathotypes with regard to host cultivars in French wheat leaf rust populations. Phytopathology 2006, 96, 264–273. [Google Scholar] [CrossRef]
  19. Germán, S.; Labes Barcellos, A.; Chaves, M.; Kohli, M.; Campos, P.; Viedma, L. The situation of common wheat rusts in the Southern Cone of America and perspectives for control. Aust. J. Agric. Res. 2007, 58, 620. [Google Scholar] [CrossRef]
  20. Elyasi-Gomari, S.; Panteleev, V.K. Virulence polymorphism of Puccinia recondita f. sp. tritici and effectiveness of Lr genes for leaf rust resistance of wheat in Ukraine. Plant Dis. 2006, 90, 853–857. [Google Scholar] [CrossRef]
  21. Bhardwaj, S.C.; Gangwar, O.P.; Prasad, P.; Kumar, S.; Khan, H.; Gupta, N. Physiologic specialization and shift in Puccinia triticina pathotypes on wheat in Indian subcontinent during 2013–2016. Ind. J. Phytopathol. 2019, 72, 23–34. [Google Scholar] [CrossRef]
  22. Bhardwaj, S.C.; Prashar, M.; Jain, S.K.; Subodh, K.; Sharma, Y.P. Physiologic specialization of Puccinia triticina on wheat (Triticum species) in India. Indian J. Agric. Sci. 2010, 80, 805–811. [Google Scholar]
  23. El-Orabey, W.M.; Omara, R.I.; Abou-Zeid, M.A. Diversity and virulence dynamics within Puccinia triticina populations in Egypt. J. Plant Protect. Pathol. 2018, 9, 735–745. [Google Scholar] [CrossRef]
  24. Nemati, Z.; Mostowfizadeh-Ghalamfarsa, R.; Dadkhodaie, A.; Mehrabi, R.; Steffenson, B.J. Virulence of leaf rust physiological races in Iran from 2010 to 2017. Plant Dis. 2020, 104, 363–372. [Google Scholar] [CrossRef] [PubMed]
  25. Wang, H.R. Notes on physiologic specialization in leaf rust of wheat in China. Phytopathology 1947, 37, 680–681. [Google Scholar]
  26. Chen, W.Q.; Hu, C.C.; Zhang, S.X. Analysis of virulence genes of Puccinia recondita f. sp. tritici population in China. Sci. Agric. Sin. 1993, 26, 17–23, (In Chinese with English abstract). [Google Scholar]
  27. Chen, W.Q.; Qin, Q.M.; Chen, Y.L.; Yan, S.B. Virulence dynamics of Puccinia recondita f. sp. tritici in China during 1992-1996. Acta Phytopathol. Sin. 1998, 28, 101–106, (In Chinese with English abstract). [Google Scholar]
  28. Hu, C.C.; Roelfs, A.P. Races and virulence of Puccinia recondita f. sp. tritici in China in 1986. Plant Dis. 1989, 73, 499–501. [Google Scholar] [CrossRef]
  29. Liu, T.G.; Chen, W.Q. Race and virulence dynamics of Puccinia triticina in China during 2000–2006. Plant Dis. 2012, 96, 1601–1607. [Google Scholar] [CrossRef]
  30. Li, H.; Huang, L.; Zhang, H.; Liu, B.; Gao, L.; Chen, W.; Liu, T. Race and virulence dynamics of Puccinia triticina in China during 2007 to 2021. Plant Dis. 2024, 108, 256–263. [Google Scholar] [CrossRef]
  31. Kolmer, J.A. Collections of Puccinia triticina in different provinces of China are highly related for virulence and molecular genotype. Phytopathology 2015, 105, 700–706. [Google Scholar] [CrossRef] [PubMed]
  32. Chen, W.Q.; Kang, Z.S.; Ma, Z.H.; Xu, S.C.; Jin, S.L.; Jiang, Y.Y. Integrated management of wheat stripe rust caused by Puccinia striiformis f. sp. tritici in China. Sci. Agric. Sin. 2013, 46, 4254–4262, (In Chinese with English abstract). [Google Scholar]
  33. Wan, Q.; Liang, J.; Luo, Y.; Ma, Z. Population genetic structure of Puccinia striiformis in Northwestern China. Plant Dis. 2015, 99, 1764–1774. [Google Scholar] [CrossRef]
  34. Aoun, M.; Kolmer, J.A.; Breiland, M.; Richards, J.; Brueggeman, R.S.; Szabo, L.J.; Acevedo, M. Genotyping-by-sequencing for the study of genetic diversity in Puccinia triticina. Plant Dis. 2020, 104, 752–760. [Google Scholar] [CrossRef] [PubMed]
  35. Park, R.; Jahoor, A.; Felsenstein, F. Population structure of Puccinia recondita in Western Europe during 1995, as assessed by variability in pathogenicity and molecular markers. J. Phytopathol. 2000, 148, 169–179. [Google Scholar] [CrossRef]
  36. Kolmer, J.A.; Liu, J.Q.; Sies, M. Virulene and molecular polymorphism in Puccinia recondita f. sp. tritici in Canada. Phytopathology 1995, 85, 276–285. [Google Scholar] [CrossRef]
  37. Kolmer, J.A.; Liu, J.Q. Virulence and molecular polymorphism in international collections of the wheat leaf rust fungus Puccinia triticina. Phytopathology 2000, 90, 427–436. [Google Scholar] [CrossRef] [PubMed]
  38. Kolmer, J.A. Molecular polymorphism and virulence phenotypes of the wheat leaf rust fungus Puccinia triticina in Canada. Can. J. Bot. 2001, 79, 917–926. [Google Scholar]
  39. Kolmer, J.A.; Ordoñez, M.E. Genetic differentiation of Puccinia triticina populations in Central Asia and the Caucasus. Phytopathology 2007, 97, 1141–1149. [Google Scholar] [CrossRef] [PubMed]
  40. Zhang, L.Y.; Meng, Q.F.; Yan, H.F.; Liu, D.Q. Virulence and molecular genetic diversity of the Puccinia triticina population in Hebei Province of China in 2008 and 2010. Eur. J. Plant Pathol. 2019, 153, 153–166. [Google Scholar] [CrossRef]
  41. Zhao, P.; Zhang, L.; Jia, S.; Yan, H.; Liu, D. Population genetic diversity of Puccinia triticina based on EST-SSR in Hebei in 2009. Acta Agric. Boreali-Occident. Sin. 2014, 23, 61–67, (In Chinese with English abstract). [Google Scholar]
  42. Wang, X.; Bakkeren, G.; McCallum, B. Virulence and molecular polymorphisms of the wheat leaf rust fungus Puccinia triticina in Canada from 1997 to 2007. Botany 2010, 88, 575–589. [Google Scholar] [CrossRef]
  43. Wang, X.; Mulock, B.; Guus, B.; McCallum, B. Development of EST-derived simple sequence repeat markers for wheat leaf rust fungus, Puccinia triticina Eriks. Can. J. Plant Pathol. 2010, 32, 98–107. [Google Scholar] [CrossRef]
  44. Kolmer, J.A.; Herman, A.; Ordoñez, M.E.; German, S.; Morgounov, A.; Pretorius, Z.; Visser, B.; Anikster, Y.; Acevedo, M. Endemic and panglobal genetic groups, and divergence of host-associated forms in worldwide collections of the wheat leaf rust fungus Puccinia triticina as determined by genotyping by sequencing. Heredity 2020, 124, 397–409. [Google Scholar] [CrossRef] [PubMed]
  45. Zhou, H.; He, X.; Wang, S.; Ma, Q.; Sun, B.; Ding, S.; Chen, L.; Zhang, M.; Li, H. Diversity of the Fusarium pathogens associated with crown rot in the Huanghuai wheat-growing region of China. Environ. Microbiol. 2019, 21, 2740–2754. [Google Scholar] [CrossRef] [PubMed]
  46. Ordoñez, M.E.; Kolmer, J.A. Virulence phenotypes of a worldwide collection of Puccinia triticina from durum wheat. Phytopathology 2007, 97, 344–351. [Google Scholar] [CrossRef]
  47. Goyeau, H.; Halkett, F.; Zapater, M.F.; Carlier, J.; Lannou, C. Clonality and host selection in the wheat pathogenic fungus Puccinia triticina. Fungal Genet. Biol. 2007, 44, 474–483. [Google Scholar] [CrossRef] [PubMed]
  48. Ordoñez, M.E.; Kolmer, J.A. Differentiation of molecular genotypes and virulence phenotypes of Puccinia triticina from common wheat in North America. Phytopathology 2009, 99, 750–758. [Google Scholar] [CrossRef]
  49. Kolmer, J.A.; Long, D.L.; Hughes, M.E. Physiologic specialization of Puccinia triticina on wheat in the United States in 2008. Plant Dis. 2010, 94, 775–780. [Google Scholar] [CrossRef]
  50. Morgounov, A.; Tufan, H.A.; Sharma, R.; Akin, B.; Bagci, A.; Braun, H.; Kaya, Y.; Keser, M.; Payne, T.S.; Sonder, K.; et al. Global incidence of wheat rusts and powdery mildew during 1969-2010 and durability of resistance of winter wheat variety Bezostaya 1. Eur. J. Plant Pathol. 2012, 132, 323–340. [Google Scholar] [CrossRef]
  51. Kolmer, J.A.; Mirza, J.I.; Imtiaz, M.; Shah, S. Genetic differentiation of the wheat leaf rust fungus Puccinia triticina in Pakistan and genetic relationship to other worldwide populations. Phytopathology 2017, 107, 786–790. [Google Scholar] [CrossRef]
  52. Kosman, E.; Ben Yehuda, P.; Manisterski, J. Diversity of virulence phenotypes among annual populations of wheat leaf rust in Israel from 1993 to 2008. Plant Pathol. 2014, 63, 563–571. [Google Scholar] [CrossRef]
  53. Kosman, E.; Ben-Yehuda, P.; Manisterski, J.; Sela, H. Diversity of virulence phenotypes among annual populations of Puccinia triticina originating from common wheat in Israel during the period 2000–2015. Plant Pathol. 2019, 68, 1741–1748. [Google Scholar] [CrossRef]
  54. Wright, S. Evolution in Mendelian populations. 1931. Bull. Math. Biol. 1990, 52, 241–295. [Google Scholar] [CrossRef] [PubMed]
  55. Li, Z.Q.; Zeng, S.M. Wheat Rusts in China; China Agricultural Press: Beijing, China, 2002. [Google Scholar]
  56. Pu, Z.; Liu, T.; Zhang, M.; Chen, W. Development of a molecular detection assay for the physiological race MFR of Puccinia triticina. Acta Phytopathol. Sin. 2004, 34, 449–454, (In Chinese with English abstract). [Google Scholar]
  57. Xu, D.; Huang, G.; Yang, W.; Liu, D. Analysis of the virulence and molecular polymorphism of Puccinia recondita f. sp tritici. J. Agric. Biotechnol. 2002, 10, 41–45, (In Chinese with English abstract). [Google Scholar]
  58. Zhang, L.; Zhang, L.; Meng, Q.; Yan, H.; Liu, D. Virulence and molecular genetic diversity, variation, and evolution of the Puccinia triticina population in Hebei Province of China from 2001 to 2010. Front Plant Sci. 2023, 14, 1095677. [Google Scholar] [CrossRef] [PubMed]
  59. Prasad, P.; Bhardwaj, S.C.; Gangwar, O.P.; Kumar, S.; Khan, H.; Kumar, S.; Rawal, H.; Sharma, T. Population differentiation of wheat leaf rust fungus Puccinia triticina in South Asia. Curr. Sci. 2017, 112, 2073–2083. [Google Scholar] [CrossRef]
  60. Gultyaeva, E.I.; Shaydayuk, E.L.; Kazartsev, I.A.; Akhmetova, A.; Kosman, E. Microsatellite analysis of Puccinia triticina from Triticum and Aegilops hosts. Australas Plant Path. 2018, 47, 163–170. [Google Scholar] [CrossRef]
  61. Chen, X.M.; Penman, L.; Wan, A.M.; Cheng, P. Virulence races of Puccinia striiformis f. sp. tritici in 2006 and 2007 and development of wheat stripe rust and distributions, dynamics, and evolutionary relationships of races from 2000 to 2007 in the United States. Can. J. Plant Pathol. 2010, 32, 315–333. [Google Scholar] [CrossRef]
  62. Kolmer, J.A.; Ordoñez, M.E.; Manisterski, J.; Anikster, Y. Genetic differentiation of Puccinia triticina populations in the Middle East and genetic similarity with populations in Central Asia. Phytopathology 2011, 101, 870–877. [Google Scholar] [CrossRef]
  63. Ordoñez, M.E.; Germán, S.E.; Kolmer, J.A. Genetic differentiation within the Puccinia triticina population in South America and comparison with the North American population suggests common ancestry and intercontinental migration. Phytopathology 2010, 100, 376–383. [Google Scholar] [CrossRef]
  64. Kolmer, J.A.; Hanzalova, A.; Goyeau, H.; Bayles, R.; Morgounov, A. Genetic differentiation of the wheat leaf rust fungus Puccinia triticina in Europe. Plant Pathol. 2012, 62, 21–31. [Google Scholar] [CrossRef]
  65. Halkett, F.; Simon, J.C.; Balloux, F. Tackling the population genetics of clonal and partially clonal organisms. Trends Ecol. Evol. 2005, 20, 194–201. [Google Scholar] [CrossRef] [PubMed]
  66. Kolmer, J.A.; Acevedo, M.A. Genetically divergent types of the wheat leaf fungus Puccinia triticina in Ethiopia, a center of tetraploid wheat diversity. Phytopathology 2016, 106, 380–385. [Google Scholar] [CrossRef]
  67. Kolmer, J.A.; Ordoñez, M.E.; German, S.; Morgounov, A.; Pretorius, Z.; Visser, B.; Goyeau, H.; Anikster, Y.; Acevedo, M. Multilocus genotypes of the wheat leaf rust fungus Puccinia triticina in worldwide regions indicate past and current long-distance migration. Phytopathology 2019, 109, 1453–1463. [Google Scholar] [CrossRef] [PubMed]
  68. Mantovani, P.; Maccaferri, M.; Tuberosa, R.; Kolmer, J. Virulence phenotypes and molecular genotypes in collections of Puccinia triticina from Italy. Plant Dis. 2010, 94, 420–424. [Google Scholar] [CrossRef] [PubMed]
  69. Ordoñez, M.E.; Kolmer, J.A. Simple sequence repeat diversity of a worldwide collection of Puccinia triticina from durum wheat. Phytopathology 2007, 97, 574–583. [Google Scholar] [CrossRef]
  70. Zhao, L.; Wang, H. Heterogeneity of virulence of the so called same race of Puccinia recondita f. sp. tritici. J. Agric. Univ. Hebei. 1986, 9, 57–65, (In Chinese with English abstract). [Google Scholar]
  71. Wang, H.R. Dialectical analysis of host and pathogen interaction in the study of leaf rust of wheat. Sci. Agric. Sin. 1988, 21, 6–11, (In Chinese with English abstract). [Google Scholar]
  72. Roelfs, A.P. Race specificity and methods of study. The Cereal Rusts: Origins, Specificity, Structure, and Physiology; Roelfs, A.P., Bushnell, W.R., Eds.; Academic Press: Orlando, FL, USA, 1984; Volume I, pp. 131–164. [Google Scholar]
  73. Gao, J. Genetic Polymorphismof Wheat Leaf Rust (Puccinia recondita f. sp. tritici) Based on Virulence and AFLP. Ph.D. Dissertation, Hebei Agricultural University, Baoding, China, 2003. [Google Scholar]
  74. Rohlf, F. NTSYS-pc, Numerical Taxonomy and Multivariate Analysis System, version 2.1; Exeter Publications: New York, NY, USA, 2000. [Google Scholar]
  75. Yeh, F.; Yang, R.; Boyle, T.; Ye, Z.H.; Mao, J.X. POPGENE, the User-Friendly Shareware for Population Genetic Analysis; Molecular Biology and Biotechnology Centre, University of Alberta: Edmonton, AB, Canada, 1999. [Google Scholar]
Figure 1. UPGMA dendrogram of Puccinia triticina based on virulence polymorphism.
Figure 1. UPGMA dendrogram of Puccinia triticina based on virulence polymorphism.
Plants 13 02438 g001
Figure 2. UPGMA dendrogram of Puccinia triticina based on EST-SSR polymorphism.
Figure 2. UPGMA dendrogram of Puccinia triticina based on EST-SSR polymorphism.
Plants 13 02438 g002
Figure 3. Two−dimensional plot of principal coordinates analysis of Puccinia triticina based on virulence polymorphism. Note: ○, Hebei; ●, Shandong; ◇, Sichuan; and ◆, Xinjiang.
Figure 3. Two−dimensional plot of principal coordinates analysis of Puccinia triticina based on virulence polymorphism. Note: ○, Hebei; ●, Shandong; ◇, Sichuan; and ◆, Xinjiang.
Plants 13 02438 g003
Figure 4. Two−dimensional plot of principal coordinates analysis of Puccinia triticina based on EST-SSR polymorphism. Note: ○, Hebei; ●, Shandong; ◇, Sichuan; and ◆, Xinjiang. Area A contains most isolates from Heibei and Shandong, area B contains most isolates from Sichuan, and area C contains most isolates from Xinjiang.
Figure 4. Two−dimensional plot of principal coordinates analysis of Puccinia triticina based on EST-SSR polymorphism. Note: ○, Hebei; ●, Shandong; ◇, Sichuan; and ◆, Xinjiang. Area A contains most isolates from Heibei and Shandong, area B contains most isolates from Sichuan, and area C contains most isolates from Xinjiang.
Plants 13 02438 g004
Figure 5. UPGMA dendrogram of P. triticina predominant pathotypes based on virulence.
Figure 5. UPGMA dendrogram of P. triticina predominant pathotypes based on virulence.
Plants 13 02438 g005
Figure 6. UPGMA dendrogram of Puccinia triticina predominant pathotypes based on EST-SSR.
Figure 6. UPGMA dendrogram of Puccinia triticina predominant pathotypes based on EST-SSR.
Plants 13 02438 g006
Figure 7. UPGMA dendrogram of Puccinia triticina populations based on EST-SSR data.
Figure 7. UPGMA dendrogram of Puccinia triticina populations based on EST-SSR data.
Plants 13 02438 g007
Table 1. Isolate codes, location, and race information of Puccinia triticina samples collected from Hebei, Shandong, Sichuan, and Xinjiang.
Table 1. Isolate codes, location, and race information of Puccinia triticina samples collected from Hebei, Shandong, Sichuan, and Xinjiang.
Isolate CodesRegionRacesIsolate CodesRegionRacesIsolate CodesRegionRaces
1Tangshan, HebeiTHTS74Yantai, ShandongPHKT147Shihezi, XinjiangTHKT
2Tangshan, HebeiTHTT75Yantai, ShandongTHTT148Urumqi, XinjiangTHPT
3Tangshan, HebeiPHHP76Jinan, ShandongPHNS149Urumqi, XinjiangTCPT
4Tangshan, HebeiTHTS77Jinan, ShandongPHKT150Urumqi, XinjiangPCTS
5Tangshan, HebeiPHST78Jinan, ShandongPHKT151Urumqi, XinjiangMRTT
6Tangshan, HebeiTHTS79Jinan, ShandongTHKT152Urumqi, XinjiangTHTT
7Tangshan, HebeiPHJT80Jinan, ShandongTHKS153Urumqi, XinjiangTHTT
8Tangshan, HebeiPHTT81Jinan, ShandongPHKT154Kuitun, XinjiangTCRF
9Qinhuangdao, HebeiFHTT82Jinan, ShandongPHKT155Kuitun, XinjiangPRTT
10Qinhuangdao, HebeiRHTT83Jinan, ShandongPHKH156Kuitun, XinjiangPGPT
11Qinhuangdao, HebeiPHKS84Jinan, ShandongTHTT157Kuitun, XinjiangTHTS
12Qinhuangdao, HebeiPGTQ85Jinan, ShandongPHTT158Urumqi, XinjiangMCPS
13Langfang, HebeiTHTS86Yantai, ShandongTHTT159Urumqi, XinjiangMHQS
14Langfang, HebeiPHTT87Yantai, ShandongTHTT160Shihezi, XinjiangPHTT
15Cangzhou, HebeiTHPT88Laiyang, ShandongTHKS161Yining, XinjiangTRRT
16Cangzhou, HebeiTHKT89Tai’an, ShandongTHTT162Yining, XinjiangTHTS
17Cangzhou, HebeiTHTS90Laiyang, ShandongPHKR163Yining, XinjiangTCTT
18Baoding, HebeiPHBT91Laiyang, ShandongPHTS164Shihezi, XinjiangTHTT
19Baoding, HebeiTGKT92Laiyang, ShandongTHTS165Shihezi, XinjiangTGTS
20Baoding, HebeiPHKT93Laiyang, ShandongTHTT166Shihezi, XinjiangPHNS
21Baoding, HebeiPHST94Laiyang, ShandongPHKT167Urumqi, XinjiangTHTT
22Baoding, HebeiTHTT95Laiyang, ShandongPHSS168Shihezi, XinjiangTKFS
23Baoding, HebeiTHTS96Yantai, ShandongTHTT169Urumqi, XinjiangPHTS
24Shijiazhuang, HebeiTKTR97Jiange, SichuanTHRT170Yining, XinjiangTHTS
25Handan, HebeiTHKT98Jiange, SichuanPHFS171Shihezi, XinjiangPHTT
26Handan, HebeiTHTS99Jiange, SichuanTCTT172Urumqi, XinjiangTHTT
27Handan, HebeiTHTT100Jiange, SichuanPHTT173Shihezi, XinjiangPHRT
28Handan, HebeiTHGT101Jiange, SichuanTHTT174Shihezi, XinjiangCHKT
29Handan, HebeiTHHS102Jiange, SichuanPHTT175Yining, XinjiangTHPP
30Xingtai, HebeiTHTT103Jiange, SichuanTHST176Shihezi, XinjiangFHTT
31Tangshan, HebeiPHTT104Jiange, SichuanTHTT177Shihezi, XinjiangTHTS
32Tangshan, HebeiPHKS105Jiange, SichuanPHTS178Shihezi, XinjiangTCTM
33Tangshan, HebeiPHTT106Jiange, SichuanTHTT179Shihezi, XinjiangTHTT
34Tangshan, HebeiTHKT107Jiange, SichuanTKFS180Shihezi, XinjiangPHRT
35Langfang, HebeiPHTS108Zitong, SichuanTCTT181Urumqi, XinjiangPHTT
36Shijiazhuang, HebeiPHKQ109Zitong, SichuanTHTT182Yining, XinjiangTHTP
37Shijiazhuang, HebeiTHTT110Panzhihua, SichuanTHKT183Yining, XinjiangTHPR
38Shijiazhuang, HebeiPHNQ111Panzhihua, SichuanTHTT184Shihezi, XinjiangTCTT
39Shijiazhuang, HebeiTHTT112Panzhihua, SichuanFHTT185Urumqi, XinjiangTHSR
40Handan, HebeiTHTT113Panzhihua, SichuanPHGT186Yining, XinjiangTHTT
41Handan, HebeiPHKS114Ya’an, SichuanTGJP187Shihezi, XinjiangTCKT
42Xingtai, HebeiTHTT115Yibin, SichuanTKKD188Shihezi, XinjiangTHTT
43Xingtai, HebeiTHSS116Zitong, SichuanPHTT189Shihezi, XinjiangTHTS
44Tangshan, HebeiTHTT117Panzhihua, SichuanTKTT190Kuitun, XinjiangPHTS
45Tangshan, HebeiPHTQ118Jiange, SichuanPHFS191Yining, XinjiangTHSR
46Handan, HebeiTHSS119Jiange, SichuanTHTT192Shihezi, XinjiangTHTT
47Langfang, HebeiPHGS120Jiange, SichuanPKPQ193Urumqi, XinjiangTHRR
48Langfang, HebeiTHTT121Jiange, SichuanTHTS194Shihezi, XinjiangTHTT
49Tai’an, ShandongPCJT122Jiange, SichuanPHKN195Kuitun, XinjiangTHKQ
50Tai’an, ShandongPHKT123Zitong, SichuanFHTR196Yining, XinjiangTHTS
51Tai’an, ShandongTHTS124Zitong, SichuanTHTT197Urumqi, XinjiangTHSR
52Tai’an, ShandongPHTT125Zitong, SichuanTHTS198Shihezi, XinjiangPHTT
53Laiyang, ShandongPCHT126Panzhihua, SichuanTKTQ199Yining, XinjiangTHPN
54Laiyang, ShandongTHTT127Panzhihua, SichuanPHTP200Kuitun, XinjiangPGST
55Laiyang, ShandongPHKS128Yibin, SichuanTHTT201Kuitun, XinjiangTHTT
56Laiyang, ShandongTHTT129Yibin, SichuanTHJT202Shihezi, XinjiangTGTQ
57Yantai, ShandongPHKT130Yibin, SichuanTHTT203Urumqi, XinjiangTHTT
58Yantai, ShandongTHRT131Yibin, SichuanTHTJ204Yining, XinjiangKHPQ
59Yantai, ShandongPHRT132Jiange, SichuanFCKT205Kuitun, XinjiangTBTT
60Qingdao, ShandongPHTT133Jiange, SichuanFHRT206Yining, XinjiangPHKQ
61Laiyang, ShandongPHKT134Jiange, SichuanTHKP207Shihezi, XinjiangTHTT
62Laiyang, ShandongTHTT135Panzhihua, SichuanPHST208Urumqi, XinjiangTHTT
63Laiyang, ShandongTHTT136Jiange, SichuanTHTT209Urumqi, XinjiangTGTP
64Laiyang, ShandongPHRS137Jiange, SichuanTHTT210Shihezi, XinjiangTHTT
65Laiyang, ShandongTHRS138Zitong, SichuanPHTT211Yining, XinjiangPGDL
66Laiyang, ShandongPHKT139Panzhihua, SichuanPHST212Shihezi, XinjiangPHTS
67Yantai, ShandongPHKT140Panzhihua, SichuanPHTT213Shihezi, XinjiangPHTT
68Yantai, ShandongTHKT141Shihezi, XinjiangTHTP214Shihezi, XinjiangTHTS
69Yantai, ShandongPHTT142Shihezi, XinjiangMGRR215Yining, XinjiangPHTT
70Yantai, ShandongTHTT143Yining, XinjiangTHTT216Urumqi, XinjiangPHRM
71Yantai, ShandongTHTT144Yining, XinjiangTHTP217Yining, XinjiangTHTP
72Yantai, ShandongTHTT145Yining, XinjiangPHTT
73Yantai, ShandongTHHT146Yining, XinjiangTCCT
Table 2. Races and number of isolates determined from Hebei, Shandong, Sichuan, and Xinjiang in 2009.
Table 2. Races and number of isolates determined from Hebei, Shandong, Sichuan, and Xinjiang in 2009.
RegionsNumber of
Isolates
Number of RacesRaces
Hebei4824THTT(10), THTS(7), PHTT(4), PHJS(3), THJT(3), PHST(2), THJS(2), PHGP, FHTT, RHTT, PGTQ, THPT, PHBT, TGJT, PHKT, PHJT, TKTR, THGT, THGS, PHTS, PHJQ, PHPQ, PHTQ, PHGS
Shandong4819THTT(14), PHJT(11), PHTT(4), THJS(2), THJT(2), THTS(2), PCJT, PCGT, PHJS, THRT, PHRT, PHRS, THRS, THGT, PHPS, PHJH, PHJR, PHTS, PHSS
Sichuan4426THTT(11), PHTT(5), THJT(3), TCTT(2), THTS(2), PHDS(2), PHST(2), THRT, PHTS, TKDS, FHTT, PHGT, TGJP, TKJD, TKTT, PKPQ, PHJN, FHTR, TKTQ, PHTP, THTJ, FCJT, FHRT, THJP
Xinjiang7741THTT(16), THTS(7), PHTT(7), THTP(4), THTR(3), PHTS(3), TCTT(2), PHRT(2), MGRR, TCBT, THJT, THPT, TCPT, PCTS, MRTT, TCRF, PRTT, PGPT, MCPS, MHRS, TRRT, TGTS, PHPS, TKDS, CHJT, THPP, FHTT, TCTM, THPR, TCJT, THRR, THJQ, THPN, PGTT, TGTQ, KHPQ, TBTT, PHJQ, TGTP, PGDL, PHRM
Note: The number in parentheses is for the number of races, while races without parentheses only have one isolate; bold is for the predominant races in each region.
Table 3. Distribution and frequencies of predominant races of Puccinia triticina among single uredinal isolates taken from leaf rust samples collected from wheat fields from Hebei, Shandong, Sichuan, and Xinjiang.
Table 3. Distribution and frequencies of predominant races of Puccinia triticina among single uredinal isolates taken from leaf rust samples collected from wheat fields from Hebei, Shandong, Sichuan, and Xinjiang.
RacesHebeiShandongSichuanXinjiangTotal
N%N%N%N%N%
THTT1020.831429.171125.001620.785123.50
PHTT48.3348.33511.3679.09209.22
THTS714.5824.1724.5579.09188.29
PHKT12.081122.9200.0000.00125.53
THKT36.2524.1712.2711.3073.23
PHTS12.0812.0812.2733.9062.76
PHKS36.2512.0800.0000.0041.84
PHST24.1700.0024.5500.0041.84
TCTT00.0000.0024.5522.6041.84
THTP00.0000.0000.0045.1941.84
FHTT12.0800.0012.2711.3031.38
PHRT00.0012.0800.0022.6031.38
THSR00.0000.0000.0033.9031.38
others1633.331225.001943.183140.267835.94
Total48...48...44...77...217...
Table 4. Numbers of alleles, polymorphic alleles, and percent of polymorphism of 18 EST-SSR primers.
Table 4. Numbers of alleles, polymorphic alleles, and percent of polymorphism of 18 EST-SSR primers.
EST-SSR LocusNumber
of Alleles a
Number of
Polymorphic Alleles b
Percent of
Polymorphic (%)
EST-SSR LocusNumber of Alleles aNumber of
Polymorphic Alleles b
Percent of
Polymorphic (%)
Ptssr00832150.00Ptssr01894250.00
Ptssr00192150.00Ptssr08012150.00
Ptssr56493266.67Ptssr04813266.67
Ptssr00853266.67Ptssr06393266.67
Ptssr29483266.67Ptssr31452150.00
Ptssr05362150.00Ptssr654233100.00
Ptssr68632150.00Ptssr01822150.00
Ptssr02432150.00Ptssr63864375.00
Ptssr01253266.67Average2.671.6762.50
Ptssr55943266.67Total483062.50
a The maximum number of alleles each pair of primers can amplify. It should be noted that not every isolate can amplify all alleles. b The number of polymorphic alleles contained in the above-mentioned alleles (alleles that cannot be amplified in certain isolates).
Table 5. Genetic diversity parameters of Puccinia triticina based on EST-SSR.
Table 5. Genetic diversity parameters of Puccinia triticina based on EST-SSR.
Genetic Polymorphic ParametersHebeiShandongSichuanXinjiangSpecies LevelPopulations Level
Observed number of alleles (Na)1.841.671.711.821.951.76
Effective number of alleles (Ne)1.671.501.491.491.621.54
Nei’s gene diversity index (H)0.360.280.280.280.350.30
Shannon’s information index (I)0.520.400.410.420.520.40
Number of polymorphic loci (NP)41.0033.0035.0040.0047.0037.00
Percent of polymorphic loci (P%)83.6767.3571.4381.6395.9276.02
Table 6. Analysis of molecular variance (AMOVA) across and within populations.
Table 6. Analysis of molecular variance (AMOVA) across and within populations.
Source of VariationDegree of
Freedom
Sum of SquaresPercentage of
Variation
p Value
Among populations3306.1425.03<0.0001
Within populations2131159.0674.97<0.0001
Total2161465.20100.00
Table 7. Nei’s unbiased measures of genetic distances and genetic identities between populations.
Table 7. Nei’s unbiased measures of genetic distances and genetic identities between populations.
PopulationsHebeiShandongSichuanXinjiang
Hebei****0.93370.93180.9229
Shandong0.0686****0.90890.8587
Sichuan0.07070.0955****0.8682
Xinjiang0.08030.15240.1413****
Note: Nei’s genetic identities (above diagonal) and genetic distances (below diagonal).
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

Zhang, L.; Zhao, P.; Meng, Q.; Yan, H.; Liu, D. The Migration, Diversity, and Evolution of Puccinia triticina in China. Plants 2024, 13, 2438. https://doi.org/10.3390/plants13172438

AMA Style

Zhang L, Zhao P, Meng Q, Yan H, Liu D. The Migration, Diversity, and Evolution of Puccinia triticina in China. Plants. 2024; 13(17):2438. https://doi.org/10.3390/plants13172438

Chicago/Turabian Style

Zhang, Lin, Panpan Zhao, Qingfang Meng, Hongfei Yan, and Daqun Liu. 2024. "The Migration, Diversity, and Evolution of Puccinia triticina in China" Plants 13, no. 17: 2438. https://doi.org/10.3390/plants13172438

APA Style

Zhang, L., Zhao, P., Meng, Q., Yan, H., & Liu, D. (2024). The Migration, Diversity, and Evolution of Puccinia triticina in China. Plants, 13(17), 2438. https://doi.org/10.3390/plants13172438

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