Virulence and molecular genetic diversity, variation, and evolution of the Puccinia triticina population in Hebei Province of China from 2001 to 2010

Wheat leaf rust, caused by Puccinia triticina, is one of the most important fungal diseases of wheat in China. However, little is known about the dynamic changes of population structure and genetic diversity of P. triticina during a period of time. In this study, 247 isolates of P. triticina collected from Hebei Province from 2001 to 2010 were tested on 36 Thatcher near-isogenic lines for virulence diversity and detected by 21 pairs of Expressed Sequence Tag derived Simple Sequence Repeat (EST-SSR) primers for genetic diversity. A total of 204 isolates were successfully identified as 164 races, and THTT, THST, PHRT, THTS, and PHTT were the most common races in Hebei Province from 2001 to 2010. The cluster analysis based on virulence showed that P. triticina has a rich virulence polymorphism, which had a certain correlation with the years, while the cluster analysis based on EST-SSR showed that the genetic diversity of the P. triticina population was significantly different between years in Hebei Province from 2001 to 2010. In addition, the population structure of P. triticina may have changed greatly in 2007 and 2009, which was significantly different from that of 2001–2006 on either virulence or genetic characteristics. The variation frequency of the population structure had an increasing trend during this period. From 2001 to 2010, there was a certain degree of gene flow among the P. triticina populations. No significant correlation was found between virulence and molecular polymorphism. The genetic differentiation analysis of the 10 tested populations (each year as a population) showed that the coefficient of genetic differentiation (Gst) was 0.27, indicating that there was a certain genetic differentiation among or within populations of P. triticina in Hebei Province. The genetic variation within populations (73.08%) was higher than that among populations (26.92%), which indicated that the genetic variations were mainly found within populations. Our study provides the foundation for a better understanding of the population structure change and genetic diversity of P. triticina over a period in Hebei Province of China.


Introduction
Wheat leaf rust, caused by Puccinia triticina Erikss., is the most common fungal disease of common wheat (Triticum aestivum L.) in major wheat-growing regions worldwide (Bolton et al., 2008). The occurrence and epidemic of leaf rust have seriously affected the yield and quality of wheat, and it is one of the main factors affecting global wheat production (Huerta-Espino et al., 2011). Due to the difference in wheat cultivars and disease period, the yield loss may be 7%-30%, even more than 50% in severe cases (Huerta-Espino et al., 2011;Zhang et al., 2022). China is the world's largest wheatproducing country, with more than 23 million hectares of wheatplanting area and more than 138 million tons yield in 2022 (USDA-FAS). A serious incidence of leaf rust with severe epidemics in 1969, 1973, 1975, 1979, 2012, 2013, and 2015 caused significant yield losses in major wheat-producing regions in China (Zhang et al., 2020a;Zhang et al., 2020b). Especially in 2015, wheat leaf rust was widespread in Hebei, Shandong, Henan, Jiangsu, Anhui, Hubei, and other provinces of China Zhang et al., 2020a;Zhang et al., 2020b;Zhang et al., 2022). Hebei Province is one of the main wheat-producing regions in China, and wheat in this region is affected by leaf rust every year. The promotion and application of resistant cultivars are the most economical, effective, and environment-friendly methods to control the disease (Pink, 2002). However, P. triticina has a wide range of temperature tolerance and highly variable virulence. These features easily lead to the loss of the effectiveness of leaf rust resistance genes, which makes leaf rust difficult to control. Therefore, it is of great significance to understand the population virulence and genetic structure of P. triticina for controlling the disease and providing guidance for breeders in wheat resistance breeding programs.
In order to monitor the virulence dynamics and genetic diversity of P. triticina, studies on race surveys and virulence identification of P. triticina have been carried out in different main wheat-producing regions worldwide since the 1920s. The annual survey of race and virulence of P. triticina in the United States began in 1926 (Mains andJackson, 1926;Johnston et al., 1968). Kolmer and Fajolu (2022) conducted a comprehensive and systematic analysis of the race and virulence dynamics of P. triticina in the United States from 2000 to 2020. Similar surveys and studies have been carried out in Canada since 1931 and from 1997 to 2010 (Johnson, 1956;Kolmer, 1991;Kolmer, 1999, Kolmer, 2001McCallum et al., 2016). Moreover, similar studies have also been carried out in France (Goyeau et al., 2006), South America (Germań et al., 2007), Ukraine (Elyasi-Gomari and Panteleev, 2006), India (Bhardwaj et al., 2010;Bhardwaj et al., 2019), and Egypt (El-Orabey et al., 2018). Wang (1947) first reported three Chinese P. triticina races that were identified based on eight international identification wheat cultivars chosen by Johnston and Mains (1932). Then, based on the Thatcher near-isogenic lines, the race and virulence studies of P. triticina since 1986 were carried out in China (Chen et al., 1993;Chen et al., 1998;Hu and Roelfs, 1989;Yang et al., 2004;Liu and Chen, 2012;Kolmer, 2015;Zhang et al., 2020a;Zhang et al., 2020b). Investigations of races and virulence of P. triticina can promote the study of genetic structure and virulence variation of P. triticina populations.
In addition to the methods of the race and virulence identification of P. triticina, many molecular marker technologies have also been applied to the genetic diversity analysis of P. triticina populations in recent years. Park et al. (2000) analyzed the population structure of P. triticina in Western Europe by assessing the variability in pathogenicity and randomly amplified polymorphic DNA (RAPD) among 61 single-uredinial isolates, and little evidence was found for robust distinct clusters among the isolates. In 2001, Kolmer (2001) found a correlation between virule nce phenotype and amplifi ed fragment le ngth polymorphism (AFLP) genotype for 69 isolates of P. triticina from Canada. In 2007, Kolmer and Ordoñez (2007) identified and analyzed the virulence and molecular genotype of isolates of P. triticina collected from common wheat in Central Asia and the Caucasus regions based on near-isogenic lines of Thatcher and SSR molecular markers and found that all populations of P. triticina from Central Asia and the Caucasus regions were significantly differentiated from that of North America. In addition,  reported a high genetic diversity in P. triticina populations collected in Hebei Province of China based on universally primed PCR (UP-PCR) primers. The Expressed Sequence Tag derived Simple Sequence Repeat (EST-SSR) markers of P. triticina have been developed and used in genetic diversity analysis for P. triticina in Canada and China (Wang et al., 2010a;Wang et al., 2010b;Zhao et al., 2014;Zhang et al., 2018). The EST-SSRs were polymorphic and informative in determining intraspecific genetic diversity (Wang et al. 2010b). P. triticina isolates collected from Hebei Province of China in 2008 and 2010 showed high diversity and differences between two populations on EST-SSR (Zhang et al., 2018). However, it only revealed the population differences between 2 years in Hebei Province, which did not fully reveal the P. triticina population dynamic patterns in this region. In this study, we analyzed the population genetic diversity of P. triticina in Hebei Province for 10 years (2001)(2002)(2003)(2004)(2005)(2006)(2007)(2008)(2009)(2010) based on the results of Zhang et al. (2018), which can provide a basis for further exploring the race variation patterns of P. triticina and the molecular mechanism of race virulence evolution.

Leaf rust sample collection and multiplication
Leaf rust samples were collected from 10 regions of Hebei Province from 2001 to 2010 during the wheat filling period by researchers from the College of Plant Protection, Hebei Agricultural University (Figure 1; Supplementary Table S1). All samples were collected from fields under natural infection. Wheat leaves with uredinia of P. triticina sampled from a single plant or cultivar were treated as a single sample. The sample treatment and singleuredinial isolate multiplication were performed as described by Zhang et al. (2020a). Urediniospores collected from each sample were inoculated to the highly susceptible cultivar 'Zhengzhou 5389' to obtain a single uredinium. Urediniospores from the single uredinium were increased by inoculating the new seedlings of 'Zhengzhou 5389' using the same inoculation procedure. Approximately 10-14 days after inoculation, urediniospores were collected by shaking the leaf, put into glass tubes for drying and lyophilization, and stored at -20°C.

DNA extraction and PCR amplification
The genomic DNA of each isolate sample was extracted from the urediniospores according to the modified Cetyltrimethylammonium Bromide (CTAB) method described by Gao (2003). The 21 pairs of EST-SSR primers (Supplementary Table S2) developed by Wang et al. (2010aWang et al. ( , 2010b were used for PCR reaction, and the total reaction system and amplification procedure were performed as described by Zhang et al. (2018). PCR products were separated by 10% (w/v) polyacrylamide gel in 0.5× TBE buffer (45 mM Tris, 45 mM boric acid, and 1 mM EDTA), visualized by silver staining.

Data analysis
To analyze the genetic relationships among selected P. triticina isolates, two binary matrices based on virulence ITs and EST-SSR data were constructed respectively as described by Zhang et al. (2018). Briefly, the virulence binary matrix of 0 and 1 was constructed based on the high (1) or low (0) ITs, and the EST-SSR binary matrix of 0 and 1 was generated based on the presence (1) or absence (0) of PCR amplification bands. The software program NTSYS-pc version 2.10 (Rohlf, 2000) was used to construct the dendrograms, the twodimensional (2D) principal coordinates analysis (PCA), and the correlation between the virulence and EST-SSR. The software POPGENE version 1.32 (Yeh et al., 1999) was used to assess the population diversity of P. triticina. To evaluate population genetic differentiations, the analysis of molecular variance (AMOVA) was performed using the AMOVA model of Arlequin3.ll software (http:// cmpg.unibe.ch/software/arlequin3).

Genetic polymorphism analysis
Two hundred forty-seven isolates of P. triticina were analyzed by 21 pairs of EST-SSR primers, and only 17 pairs of primers can amplify clear and stable bands (Supplementary Figure S2). In addition, three of the 17 pairs of primers, PtSSR0536, PtSSR3145, and PtSSR0801, did not amplify the polymorphism bands in the tested isolates. A total of 54 alleles were amplified from the remaining 14 pairs of primers, of which 42 alleles were  polymorphism alleles, and the percentage of polymorphism was 77.78%. The number of alleles and polymorphic alleles of these 14 pairs of primers was 2-9 and 1-8, respectively (Table 3). The genetic similarity coefficient based on EST-SSR analysis of the 247 isolates was 0.67-0.99 (Supplementary Figure S3), which was slightly higher than the virulence similarity coefficient. All isolates were clustered into four groups (S 1 -S 4 ) at a similarity coefficient of 0.74. Cluster S1 contained 214 isolates (86.64%), including two distinct subclusters (S 1-1 and S 1-2 ) at a similarity coefficient of 0.758. The isolates of subcluster S 1-1 and subsubcluster S 1-2-1 mainly came from 2001 to 2006, those of S 1-2-2 were from 2009, those of S 1-2-3 were from 2007, those of S 1-2-4 and S 1-2-6 were from 2008, and those of S 1-2-5 were from 2010. While the isolates of clusters S 2 -S 4 were from 2007 to 2009. These results indicated that the molecular genetic structure of P. triticina in Hebei Province from 2001 to 2006 was relatively consistent, and the genetic structure of P. triticina after 2007 had significantly changed compared with 2001-2006. Most of the isolates originating from the same or neighboring years were clustered into the same group or several neighboring subclusters with a highlevel similarity, indicating that the genetic structure in these isolates might be related to their occurrence year.    relatively macro level. In addition, the correlation coefficient of virulence and EST-SSR polymorphisms was only 0.006 by using MXCOMP module of NTSYS-pc, which indicated that the correlation between virulence and EST-SSR polymorphisms was very low. In the 2D-PCA plot analysis of EST-SSR, all isolates were mainly distributed in two areas (A and B) (Figure 4). The isolates from 2001 to 2006 were mainly concentrated in area A, which indicated that these isolates have similar genetic loci. While the isolates in 2007-2010 were mainly scattered and distributed in area B, which may suggest that the genetic diversity of P. triticina populations in 2007-2010 was relatively rich. Moreover, there was a partial overlap between area A and area B, which indicated that some isolates in P. triticina populations had a similar genetic structure in some years. The isolates of P. triticina in 2007-2010 were significantly different from those of other years, which may be due to the large changes in the genetic structure of P. triticina populations in these years.

Virulence and EST-SSR polymorphism analysis of predominant races
In order to further determine the relations of P. triticina pathotypes with virulence polymorphism or EST-SSR polymorphism, the cluster analysis was carried out on 49 isolates of five predominant races (THTT, THST, THTS, PHRT, and PHTT). In the polymorphism analysis based on virulence phenotypes, most of the same races tended to group in the same or adjacent clade ( Figure 5). While in the polymorphism analysis based on EST-SSR, the races were clustered in different branches by year (Figure 6), which indicated that there was a certain genetic differentiation among or within the P. triticina population in different years. In addition, the correlation coefficient between pathotypes and virulence polymorphisms was 0.574 (>0.3), and the correlation coefficient between pathotypes and EST-SSR polymorphisms was 0.025 (<0.3), which was analyzed by the MXCOMP module of NTSYS-pc. These results indicated a correlation between pathotypes and virulence polymorphism, but there was no correlation between pathotypes and EST-SSR polymorphism.

Genetic diversity analysis of P. triticina populations
The P. triticina isolates tested from 10 years (2001)(2002)(2003)(2004)(2005)(2006)(2007)(2008)(2009)(2010) in Hebei Province have been designated as 10 corresponding populations by the year. The H of these P. triticina populations ranged from 0.13 to 0.26, and I ranged from 0.20 to 0.38, while the highest and lowest levels of genetic diversity were found in 2009 and 2010, respectively (Table 4) In addition, the analysis results of POPGENE software showed that the coefficient of genetic differentiation (Gst) was 0.27, of which genetic diversity within the population accounts for 73.08% of the total genetic diversity, and genetic diversity among the populations accounts for 26.92% of the total population genetic diversity. There were also certain genetic differentiations among or within P. triticina populations analyzed by the Arlequin3.ll software, and the genetic variation within populations accounted for 72.54% of the total variation, and the genetic variation among the populations accounted for 27.46% (Table 5). So, the variations were mainly found within populations of P. triticina. Moreover, a certain level of  To clarify the genetic relationship between different P. triticina populations from 2001 to 2010, the population genetic distance was estimated using the POPGENE software (Table 6). A dendrogram (UPGMA, unweighted average method) was then constructed by MEGA 6.0 software based on the Nei's genetic distances (Figure 7). The 10 populations of P. triticina could be clustered into three main branches (I-III), 2001-2006, 2007-2008, and 2009-2010

Discussion
The most common method to study population dynamics of P. triticina is to monitor the virulence changes based on the infection types of isolates to a set of Thatcher near-isogenic lines or currently planted cultivars, which can directly reflect the changes of pathotypes and virulence of P. triticina. While molecular marker technology can reveal the genetic structure, genetic differentiation, and genetic diversity of pathogens at the molecular level, which can reflect the difference between individuals of pathogens based on the variation of gene or genotype (Park et al., 2000;Xu et al., 2002;Mebrate et al., 2006;Kolmer and Ordoñez, 2007;Mantovani et al., 2010;Ordoñez et al., 2010). In addition, molecular marker technology is also widely used in the virulence variation, origin evolution, and migration of pathogens (Chen et al., 1999;Xu et al., 2009). In this study, the genetic diversity, evolution, and variation of P. triticina from 2001 to 2010 in Hebei Province were researched by combining virulence identification with EST-SSR analysis.
A relatively rich genetic diversity was found among the P.  (Figure 7). In general, the factors affecting the genetic structure of P. triticina populations are mainly due to the mutation of the pathogen, the alternation of host cultivars, and the climate environment changes. According to our , and so on were introduced into Hebei Province and have a wide planting area for high yield characterization. The alternation of wheat cultivars in this periods may produce affections on P. triticina population structure changes. In addition, the precipitations, temperatures, and altitude can also affect the development of P. triticina populations . According to the Climate Bulletin of Hebei Province, the annual average temperature in Hebei Province was 0.6°C-1.3°C higher than normal in 2007-2010, which were the warmest years since meteorological observations, and the winter temperature was significantly higher, which was the second warmest winter since 1951. Warm winter will be beneficial to the overwintering and survival rate of P. triticina, which may be conducive to the emergence of new races and the development of the population size of P. triticina. Since the period of large changes in the population structure of P. triticina was consistent with the changes in climate conditions, the warm winter may be another reason for the changes in the population structure of P. triticina in 2007-2010. Based on the above analysis, we speculated that the changes in P. triticina population structure during 2001-2010 were caused by the combined effects of host cultivars and climate conditions. Kosman et al. (Kosman et al., 2014;Kosman et al., 2019) and Kolmer et al. (2017) thought that migration was also likely an important factor in the genetic structure change of P. triticina populations. Wright (1990) proposed that when parameter N m (gene flow) >1, there was gene flow between populations, which could prevent genetic differentiations between populations. Gene flow has the role of homogenization, and the stronger the gene flow, the higher the similarity and the smaller the differentiation. On the contrary, the lower the gene flow, the lower the similarity and the larger the differentiation. The gene flow generated by the spread and diffusion of spores in different regions may change the gene frequency of the population. In addition, the interannual gene flow was mainly caused by asexual spores overwintering and oversummering, which was generally relatively stable, but low temperature or drastic environmental changes may also cause variations. The north and northwest of Hebei Province are located in mountainous areas, which have fewer wheat-planting areas, so no samples were collected for identification. The main wheat-growing region of Hebei Province is in Phylogenetic tree of the 10 Puccinia triticina populations based on EST-SSR. The isolates of each year are regarded as a population. the northern part of the North China Plain (only including Hebei Province), which is narrow compared with the southern part of the North China Plain (including Henan Province, Shandong Province, Anhui Province, and North part of Jiangsu Province). So, the plain is very conducive to the migration of P. triticina by the monsoon, and the gene flow caused by the migration of P. triticina among regions may reduce the diversity of P. triticina in the same year. Therefore, the difference of P. triticina population in the same year may come from the difference of wheat cultivars planted in different regions. Wang et al. (Wang et al., 2010a, Wang et al. 2010b) developed 21 pairs of EST-SSR primers, which were effective in analyzing the molecular diversity of P. triticina isolates in Canada. However, similar to the study of Zhao (2013), four pairs of primers in the present study cannot amplify clear and stable bands. In addition, although three pairs of primers can amplify clear and stable bands, there are no polymorphic bands, which cannot distinguish the isolates of P. triticina. So, only 14 pairs of EST-SSR primers in this study were adapted to the genetic diversity analysis of P. triticina in China. Due to the significant difference in virulence of P. triticina between North America and China, we speculated that the genetic difference of P. triticina between the two regions may be the main reason for this phenomenon.
Many studies on the genetic diversity of P. triticina populations showed that the genetic variations were mainly found within populations, which was higher than among populations Ordoñez and Kolmer, 2007;Ordoñez et al., 2010;Dadrezaie et al., 2013;Xu et al., 2013). Moreover, similar results were reported in the Puccinia striiformis f. sp. tritici (Wan et al., 2015). In this study, we found a certain degree of genetic differentiation among populations or within populations. The variation within populations was higher than that among populations of P. triticina in Hebei Province from 2001 to 2010. These results were consistent with the above studies. This indicated that the differentiation within populations in the same year played a more important role than the differentiation among populations, revealing that year separation was also apparently an important limiting factor to gene flow and population structure. Because the overwintering survival rate of P. triticina in Hebei Province is relatively low, the source of the P. triticina inoculum in Hebei Province, particularly north of Shijiazhuang, may have migrated primarily from other major wheat-growing regions that are conducive to P. triticina overwintering, such as Henan, Shandong, and Jiangsu Province, which are located in the middle-south of the North China Plain. The main wheat cultivars in these provinces are different. For example, the wheat cultivars of Henan Province are mainly 'Zhengmai' series and 'Zhoumai' series, the wheat cultivars of Shandong Province are mainly 'Jimai' series, 'Lumai' series, and 'Yanmai' series, and the wheat cultivars of Jiangsu Province are mainly 'Sumai' series and 'Yangmai' series, which also easily lead to the diversity of P. triticina races. Therefore, the genetic variation of P. triticina in Hebei Province from 2001 to 2010 was mainly found within populations, which may be caused by the differences of races from other regions outside Hebei Province mentioned above. In addition, according to the study of Dadrezaie et al. (2013), although the high molecular variability within P. triticina populations may indicate the colonial behavior of P. triticina in different years, the similarity between P. triticina populations also indicates the gene flow. So the population similarity between years may be the result of gene flow.
According to the virulence and genetic polymorphism analysis, we found that the diversity of P. triticina in Hebei Province was rich, and the virulence structure was very complex. Both virulence and genetic polymorphism analyses showed that the clustering of some P. triticina isolates in Hebei Province had a certain correlation with the year, which indicated that the genotypes of the P. triticina isolates in the same year had more similarity than that in different years. Both analyses showed significant changes in the population structure of P. triticina in Hebei Province after 2007 compared with 2001-2006. With the change in climate, the occurrence of wheat leaf rust in China has been increasing yearly, especially after 2008, such as -2009-2013(Zhang et al., 2020aZhang et al., 2020b). The increase in P. triticina population size will increase the possibility of survival of individual isolates with favorable mutations for virulence or at genetic loci (Kolmer, 2015), which will further enrich the population diversity of P. triticina. In addition, we found that there was no significant correlation between these isolates and the regions, which further indicated the migration of P. triticina isolates among these regions.
Moreover, there is a contradiction in the correlation between virulence phenotypes and molecular genotypes of P. triticina in the world. Many studies on genetic diversity have reported the significant correlation between the virulence phenotypes and molecular genotypes of P. triticina in Pakistan, Europe, Central Asia, and the Caucasus regions (Halkett et al., 2005;Kolmer and Ordoñez, 2007;Mantovani et al., 2010;Kolmer et al., 2012;Kolmer et al., 2017). However, Park et al. (2000); Prasad et al. (2017), and Gultyaeva et al. (2018) thought that there was no significant correlation between virulence polymorphism and molecular polymorphism of the P. triticina population diversity in Western Europe, South Asia, or Russia and Kazakhstan. There was also no significant correlation between virulence polymorphism and molecular polymorphism of the P. triticina population in China (Xu et al., 2002;Pu et al., 2004;Xu et al., 2013;Zhao, 2013), which was consistent with the results of this study. In addition, the same conclusion had been reported in the studies of population genetics of P. striiformis f. sp. tritici in the United States and Northwestern China (Chen et al., 2010;Wan et al., 2015). Whereas Kolmer (2015) reported that the SSR genotypes of 100 P. triticina isolates from seven provinces (the samples provided by our lab) in 2009 and 2010 in China had a significant correlation with the virulence phenotypes, which may be related to the sample size or the regional source differences of the collections of P. triticina. The greater regional disparity, the greater possibility of virulence and genetic variation of P. triticina populations.
In this study, the results of polymorphism analysis of the predominant pathotypes showed a correlation between pathotype and virulence polymorphism. However, the different pathotypes may also be grouped into the same cluster. For example, THTT and PHRT or PHTT and PHRT, which differed in virulence on Lr2a and/or Lr17, were clustered into the same cluster. Similar phenomena have been found in other studies. Wang et al. (2010a) found that TDBG and TDBJ, which differed in virulence on Lr14a, were clustered into the same cluster. While some of the same races did not completely cluster into the same group, which indicated that these isolates have genetic differences, although they were of the same pathotype. Zhao and Wang (Zhao and Wang, 1986;Wang, 1988) found the virulence heterogeneity in the same pathotype of P. triticina, and the pathogenicity of different P. triticina isolates of the same pathotype was not completely consistent. The more differential hosts were used, the more pathotypes of P. triticina were identified. In this study, in addition to 16 standardized differentials that were used for race naming, 20 extra differentials were used to test the virulence of P. triticina isolates. Therefore, although the same races had the same infection types to those of 16 standard hosts, there may be differences in the infection types to additional extra differentials. For example, isolate THTT 10-2 was virulent to Lr15 and avirulent to Lr28 and Lr36, but THTT 09-25 was avirulent to Lr15 and virulent to Lr28 and Lr36. THTT 09-2 was avirulent to Lr28 and Lr33+34 and virulent to Lr15 and Lr36, but THTT 08-27 was virulent to Lr28 and Lr33+34 and avirulent to Lr15 and Lr36. Therefore, although these isolates are of the same pathotype according to the 4-letter nomenclature, some of them were different in the virulence to other Lr genes. With the increase of the differential hosts for pathotype identification, the diversity of isolates may be more abundant.

Conclusion
In the present study, we identified the virulence diversity of 247 P. triticina isolates collected in Hebei Province from 2001 to 2010 on 36 Thatcher near-isogenic lines and detected the EST-SSR genetic diversity. The cluster analysis based on virulence phenotype and EST-SSR genotype showed that P. triticina populations in Hebei Province have rich virulence and genetic diversity. Due to the alternation of wheat cultivars and climatic condition changes, the population structure of P. triticina has changed significantly since 2007. Moreover, there was a certain genetic differentiation among or within populations of P. triticina in Hebei Province, and the genetic variation of within populations was the main origin of virulence phenotype or genetic genotype variations in P. triticina populations in Hebei Province. Our study will provide a theoretical basis for a better understanding of the population structure and genetic variations of P. triticina over a period in Hebei Province of China.

Data availability statement
The original contributions presented in the study are included in the article/Supplementary Material. Further inquiries can be directed to the corresponding authors.