Next Article in Journal
Importance of Clinical Isolates in Cryptococcus neoformans Research
Next Article in Special Issue
Special Issue “Genomics of Fungal Plant Pathogens”
Previous Article in Journal
Advances and Challenges in CRISPR/Cas-Based Fungal Genome Engineering for Secondary Metabolite Production: A Review
Previous Article in Special Issue
The Chromosome-Scale Genomes of Exserohilum rostratum and Bipolaris zeicola Pathogenic Fungi Causing Rice Spikelet Rot Disease
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Virulence and Genetic Types of Blumeria graminis f. sp. hordei in Tibet and Surrounding Areas

1
Research Institude of Tibet Plateau Ecology, Plant Sciences College, Tibet Agricultural and Animal Husbandry University, Nyingchi 860000, China
2
Institute of Plant Protection, Sichuan Academy of Agricultural Sciences/MOA Key Laboratory for Integrated Management of Pest on Crops in Southwest China, Chengdu 610000, China
*
Authors to whom correspondence should be addressed.
J. Fungi 2023, 9(3), 363; https://doi.org/10.3390/jof9030363
Submission received: 4 November 2022 / Revised: 10 March 2023 / Accepted: 14 March 2023 / Published: 16 March 2023
(This article belongs to the Special Issue Genomics of Fungal Plant Pathogens)

Abstract

:
Barley (Hordeum vulgare L.) is the most important cereal crop in the Qinghai-Tibet Plateau, and the yield has been seriously threatened by Blumeria graminis f. sp. hordei (Bgh) in recent years. To understand the virulence and genetic traits of different Bgh populations, 229 isolates of Bgh were collected from Tibet, Sichuan, Gansu and Yunnan provinces of China during 2020 and 2021, and their pathogenicity to 21 barley lines of different genotypes was assessed. A total of 132 virulent types were identified. The Bgh isolates from Yunnan showed the highest diversity in terms of virulence complexity (Rci) and genetic diversity (KWm), followed by those from Sichuan, Gansu, and Tibet, in that order. Single nucleotide polymorphism (SNP) in genes coding for alternative oxidase (AOX), protein kinase A (PKA), and protein phosphatase type 2A (PPA) were detected at seven polymorphic sites. Nine haplotypes (H1–H9) with an average haplotype diversity (Hd) and nucleotide diversity π of 0.564 and 0.00034, respectively, were observed. Of these, haplotypes H1 and H4 accounted for 88.8% of the isolates, and H4 was predominant in Tibet. Genetic diversity analysis using the STRUCTURE (K = 2) and AMOVE indicated that the inter-group variation accounted for 54.68%, and inter- and intra-population genotypic heterogeneity accounted for 23.90% and 21.42%, respectively. The results revealed the recent expansion of the Bgh population in Tibet, accompanied by an increase in virulence and a loss of genetic diversity.

1. Introduction

Barley (Hordeum vulgare L.) is the fourth largest cereal crop in the world [1,2]. In China, the area for barley cultivation averaged 8 Mha during 1914 to 1918, dropped to 0.85 Mha in 2004, and further declined to 0.50 Mha in 2007 [3]. In 1997~2000, the barley production was mainly concentrated in the middle and lower reaches of the Yangtze River and North China. Then, the main barley production areas translocated to Southwest and Northwest China, and Jiangsu, Yunnan, Sichuan, Gansu and Inner Mongolia became the top five provinces for barley production in China [4]. Accounting for nearly 40% of the barley cultivation area in China, Yunnan, Sichuan and Tibet in the southwest are neighbored by Qinghai and Gansu in the northwest, geographically segregated from other regions of scaled barley production in China [5]. In the Tibet-Qinghai plateau and the surrounding Sichuan, Gansu, and Yunnan, the production areas covered approximately 0.03 Mha [6,7]. Unlike other types of barley mainly used for fodder and in the brewery industry, the hull-less barley (H. vulgare L. var. nudum), known as Qingke, is more suitable to the harsh environment in the cool, high altitude regions of the Qinghai-Tibet Plateau, and is consequently cultivated as the staple food of people living there [8].
Barley powdery mildew is an airborne fungal disease caused by Blumeria graminis f. sp. hordei (Bgh), which frequently occurs in Europe, North Africa and China [9,10,11]. It damages barley quality and results in a yield loss of generally 20% and more than 30% in outbreaks [12,13,14]. In contrast to many other plant pathogens, Bgh is extraordinarily adaptable, and some commonly recommended strategies involving genetic resistance may be ineffective [15]. As one of the consequences of global climate change and warming, powdery mildew has become more prevalent in the Tibet-Qinghai plateau and surrounding areas, and it urgently needs to be controlled [16,17]. Among the feasible strategies, breeding and deploying resistance varieties is the most efficient, economic and environment-friendly way to control the disease. Over 50 genes, including alleles, have been identified, and are available for resistance breeding [18]. However, the variation in the composition of the pathogen population, especially in its virulence to different genes can frustrate such effort, and the breakdown in resistance caused by such changes has caused the bottleneck of resistance breeding and deployment. To obtain virulence and genetic information of the Bgh populations in Tibet and the surrounding areas, 12 disease nurseries were set up in 2020 and 2021. Together with the observation of the resistance expression of barley varieties and lines bearing different resistance genes (data published elsewhere), 321 Bgh isolates derived from single colonies from disease samples were collected and characterized for their virulence and SNP polymorphism of genes coding for alternative oxidase (AOX), protein kinase A (PKA), and protein phosphatase type 2A (PPA) [19].

2. Materials and Methods

2.1. Lines and Varieties of Barley

Seeds of the susceptible control Hua 30 and 30 near-isogenic lines (NILs) harboring known resistance genes to Bgh (Table A1) were provided by Dr. Lin at the Institute of Plant Protection, China Academy of Agricultural Sciences. The susceptible variety QB01 (Zangqing 23) was a hull-less barley, provided by the Tibet Academy of Agriculture and Animal Husbandry Sciences for the isolation, purification and multiplication of the Bgh isolates.

2.2. Disease Samples and Bgh Isolates

Twelve disease nurseries including five in Tibet, three in Sichuan, two in Qinghai, and one in both Gansu and Yunnan Province (details shown in Figure A1), were set up. Seeds of each line or variety were sown in a 1.5 m long line, without replicates. The susceptible variety Hua 30 was sown to surround the disease nurseries, as the spreading line. Symptoms of powdery mildew developed at nine disease nurseries, at Qamdo (30°56′33″ N, 97°21′42″ E, 3160 m asl), Linzhi (29°43′9″ N, 94°34′46″ E, 2998 m asl), Shannan (29°10′9″ N, 91°45′55″ E, 3595 m asl), Lhasa (29°16′32″ N, 90°23′57″ E, 3662 m asl), and Rikaze (29°10′17″ N, 89°5′15″ E, 3820 m asl) in Tibet (5/5), Hezuo (35°5′53″ N, 102°54′29″ E, 2730 m asl) in Gansu (1/1), Chengdu (30°47′55″ N, 103°54′2″ E, 549 m) and Xichang (27°56′59″ N, 102°10′36″ E, 1527 m) in Sichuan (2/3) and Diqing (27°48′39″ N, 99°39′28″ E, 3280 m asl) in Yunnan (1/1). Diseased leaves were sampled from the spreading lines in the disease nurseries, and a total of 229 single-colony Bgh isolates were obtained and kept on the seedlings of the barley variety QB01 in test tubes of 5 cm in diameter containing sterilized soil, following the methods of Xu et al. [20].

2.3. Spore Preparation

To multiply the Bgh isolates, conidia of each isolates from the above seedlings were spread onto the seedlings of QB01 grown in pots containing 400 mL sterilized soil and then incubated in a growth chamber at 18 ± 2 °C, with continuous illumination. To prevent cross contamination, seedlings were completely encased by a glass cylinder of 10 cm in diameter, with five layers of gauze on top. When white mycelia emerged 3–5 d post inoculation (dpi), leaf segments of 5 cm in length were cut off the inoculated leaves, placed face-up onto 1% agar plates amended with 60 mg/L benzimidazole, and grown at 16 ± 2 °C, with a 16/8 h (light/dark) illumination regime. After 5–7 d of incubation, the conidia were collected onto parchment paper in the laminar flow and placed into 2.0 mL centrifuge tubes. The fresh conidia were used for the virulence assay, or first frozen in liquid nitrogen and then conserved at −80 °C for further DNA extraction [21].

2.4. Virulence Typing

Clumps of five germinated seeds of QB-01 and Hua 30, as well as those of 30 barley NILs, were sown in 7 rows in a 26 cm × 18 cm × 2 cm enamel tray containing sterilized soil without replicates. At the 2-leaf stage, each tray was transferred into steel cylinders of 0.4 m in diameter and 1 m in height, and inoculated with the conidia of one Bgh isolate by gently shaking the sporulating leaf segments. The tray was then encased in a transparent plastic bag supported by an iron framework, and incubated for 7–10 d at 18 ± 2 °C with continued illumination until susceptible QB-01 and Hua 30 showed symptoms. The infection types (ITs) of each barley NIL from different Bgh isolates were scored on a scale of 0 to 4, according to Jensen et al. [22], and resistance/susceptibility responses were obtained, i.e., 0–2 for resistant (R) and 3–4 for susceptible (S). The 21 selected NILs, shown in Figure 1, were clustered into groups using an octal numbering system. These groups were ranked according to the sum of three values of each group, and the resulting number (reverse-octal notation) was used to designate the pathotype of the isolates [23,24].

2.5. DNA Extraction and PCR Amplification

Genomic DNA of Bgh isolates was extracted using a Plant Genomic DNA kit DP305 (Tiangen Biotech, Beijing, China), according to the manufacturer’s instructions. Three primer sets specific for alternative oxidase (AOX; F/R-CGCATAGCCCTTTACTTAG/ATGGATTTGGGTCCTCGTT), protein kinase A (PKA; F/R-ATTTCGGTAGGGTTCATCTGG/TACCGTTCCGTCTCTTCAGG), and protein phosphatase type 2A (PPA; F/R-TAGATGGGTGGATTGAGAAC/ATCGTCAGGATCAGACCATA) were used to study the single nucleotide polymorphism (SNP) [25]. The 25 μL PCR reaction system contained 2 × Taq PCR MasterMix 12.5 μL, template DNA (10 ng/μL) 2.5 μL, with each primer (10 μmol/L) 1.0 μL, and ddH2O 8.0 μL. The amplification program consisted of a pre-denaturation at 95 °C for 5 min, followed by 34 cycles of denaturation at 95 °C for 25 s, annealing at 56 °C for 25 s, extension at 72 °C for 45 s, and a final extension at 72 °C for 7 min. After that, 5 μL of the PCR products were detected using 1% agarose gel electrophoresis. The PCR products with clear and bright bands and of the expected size were sent to Tsingke Biotechnology (Beijing, China) for sequencing.

2.6. Data Processing

The QD-SHAN in NTSYSpc 2.10s package was used for the clustering analysis of the 21 differential lines of pallas, based on the virulence assessment data of 229 Bgh isolates. Principal coordinate analysis (PCoA) was performed using an online tool (https://hiplot.com.cn accessed on 4 November 2022). VAT software was used to identify the Bgh pathotype based on the parameters related to the virulence diversity, including average relative-virulence complexity (Rci), gene diversity (Hs), genotypic diversity (Sh) and genetic diversity (KWm) [26]. DNAMAN was used to assemble the nucleotide sequences of 321 isolates, using previously published sequences as references, and the AOX-PKA-PPA sequences of the isolates was aligned and compared using the MEGA7.0. Using the Euclidean squared distance matrix in the Arlequin 3.5 package, a gene haplotype network was constructed using the PopART software, and the geographic distribution of haplotypes was processed using the ArcGIS10.2. STRUCTURE 2.3.4 under K = 1–9 was used to predict population structure, and the ΔK value was used to determine the optimal number of populations [27]. Finally, a histogram of the population-structure-related data was generated, using the Distruct software. DnaSP 5.0 [28] was used to compute diversity indices including haplotype diversity (Hd) and nucleotide diversity (Pi), to conduct neutrality tests including Tajima’s D as well as Fu and Li’F* tests, and to infer population expansion. AMOVA (Analysis of Molecular Variance) in the Arlequin 3.5 package [29,30] was used to compute the genetic differentiation fixation coefficient FCT/FST of the inter-geographic groups/inter-populations, and to conduct a neutrality test with 1000 random permutations.

3. Results

3.1. Virulence Diversity of Bgh Populations

3.1.1. Virulence Frequencies of Bgh

The virulence frequencies to different NILs as well as Hua 30 of the Bgh isolates from different locations were obtained as seen in Table 1. The virulence frequencies of all of the sampled Bgh population to P01 (Mla1 + MlaAl2), P03 (Mla6 + Mla14), P04A (Mla7 + Mlk), P04B (Mla7 + MlaNo3), P06 (Mla7 + M1Lg2), P07 (Mla9 + Mlk), and P08B (Mla9) was 0. Significant differences in virulence to the other NILs were found among different Bgh populations. Virulent isolates to P08A (Mla9 + Mlk) and P10 (Mla12 + MlaEm2) emerged only at 1 or 2 locations, while virulent isolates to P09 (Mla10 + MlaDu2), P12 (Mla22), P15 (MlRu2), P18 (Mlnn) P20 (Mlat + Mla8), P23 (MlLa) and P30 were present in all of the populations. The virulence frequencies to P02 (Mla3), P05, P13 (Mla23), P16 (Mlk) and mlo5 were less than 20%. The virulence frequency to the other eight resistance genes ranged from 25% to 50%, indicating that the NILs possessing these genes had moderate susceptibility (MS) or moderate resistance (MR). The genes in P09 (Mla10 + MlaDu2), P12 (Mla22), P15 (MlRu2), P18 (mlnn), P30 and pallas were no longer efficient. Based on the resistance frequencies of different NIL lines and the clarity of their resistance genes, 21 NIls were selected as the differential lines for virulence typing (Figure 1).

3.1.2. Virulence Types of Bgh

Based on the virulence spectra of 229 Bgh isolates to 21 barley NILs, a total of 132 pathotypes were designated (Table A2), with an average of 1.7 isolates per pathotype. Of these, 16 pathotypes had frequencies of higher than 2, and most of the isolates had one specific pathotype. Pathotypes 7.4.0.0.0.0.0, 7.3.0.0.0.0.0, 7.2.0.0.0.0.0, 7.1.0.0.0.0.0, 7.0.0.0.0.0.0, 5.4.0.0.0.0.0, 3.5.0.0.0.0.0, and 3.3.0.0.0.0.0 were common in the Bgh isolates from Tibet, and the pathotype 7.3.0.0.0.0.0 had the highest frequency, of 5. Pathotypes 7.3.6.5.0.0.0, 7.0.1.0.0.0.0, 6.4.2.0.1.0.0, and 6.1.2.0.0.0.0 were popular in the Bgh populations of Sichuan. The PCoA of the Bgh (Figure 2) was performed to represent the distance between the population of isolates from the nine disease nurseries, indicating that the Bgh populations of Diqing in Yunnan province, Hezuo in Gansu province and Chengdu and Xichang in Sichuan province had similar virulence spectra, which were different from those of the four Tibetan Bgh populations.

3.1.3. Virulence Diversity of Bgh

The virulence diversity of the Bgh populations from nine disease nurseries in Tibet, Sichuan, Yunnan, and Gansu were evaluated. As shown in Table 2, the Bgh population at Diqing in Yunnan had the highest virulence complexity, with an Rci value of 0.363, whereas that of Rikaze in Tibet was the lowest (0.160). The Bgh populations at Diqing in Yunnan and Xichang in Sichuan had genetic diversity value of KWm at 0.385, while the populations at Shannan in Tibet had the lowest KWm (0.143). The relationship between the populations was analyzed based on the genetic distance Nei’s (N) and mean character difference (MCD) (Table 3). The genetic distance between the populations at Diqing and Shannan was greatest for N (0.192), and the MCD value (0.292) between these two population was also the highest. The N and MCD values between the populations at Lhasa, Shannan, and Qamdo in Tibet were the lowest, and were 0.017 and 0.082, respectively.

3.2. Molecular Genetic Structure of the Bgh Populations

3.2.1. Haplotype Composition of Different Bgh Population

After the sequences of the genes of AOX, PKA and PPA in the genome of 321 isolates were organized through alignment analysis and assembled, the length of the sequence containing polymorphic sites was 2196 bp for all the isolates. Excluding gaps and three deletion sites, a total of seven polymorphic sites, two singleton variable sites (at positions 394 and 1883), and five parsimony informative sites (at positions 238, 738, 1279, 1539, and 1569) were found in the sequences. In total, nine haplotypes (H1–H9) were detected (Table 3). The haplotype network and the geographical distribution of the isolates (Figure 3 and Figure 4) indicated that, among the nine haplotypes, two haplotypes (H1 and H4) were dominant and accounted for 29.6% and 59.2% of the isolates, respectively. H1 existed in populations in Sichuan, Gansu, and Yunnan and the population at Qamdo in Tibet, while H4 haplotype existed in all of the five Tibetan populations and in the population at Chengdu in Sichuan.
The genetic structure of the Bgh populations was analyzed employing STRUCTURE software, and the K values were set from 1 to 9, to predict the optimal population structure. When K = 2, the ΔΚ value was the largest, suggesting that the optimal number of sub-populations was 2 for the grouping (Figure 5). The nine Bgh populations were classified into two sub-populations, of which subpopulation 1 involved all of the Tibetan populations, and subpopulation 2 involved the populations from the surrounding Sichuan, Gansu, and Yunnan. Populations from closer geographic locations had higher genetic-structure similarity.

3.2.2. Genotype Diversity of Bgh Populations

The total haplotype diversity (represented by Hd) of the investigated Bgh populations was 0.564, and the average nucleotide diversity (represented by π) was 0.00034. The Hd value ranged from 0 to 0.6485, and the π value ranged from 0 to 0.00049, illustrating the large variation between the populations. The population at Qamdo in Tibet had the highest Hd (0.64848), followed by the population at Hezuo in Gansu and Xichang in Sichuan, which were 0.34762 and 0.31169, respectively. The populations at Lhasa, Linzhi, and Shannan in Tibet and Diqing in Yunnan showed the lowest haplotype and nucleotide diversity, and only one haplotype was detected in each of these populations There were two haplotypes in the population at Rikaze (Table 4 and Table 5).
The AMOVE analysis showed (Table 6) that inter-group genetic variation accounted for 54.68%, inter-population genetic variation accounted for 23.90% and intra-population genetic variation accounted for 21.42%. The inter-group genetic differentiation coefficient, FCT, was 0.54681, and the average inter-population genetic differentiation coefficient, FST, was 0.69714, indicating that larger genetic differentiation of Bgh populations mainly existed between groups. There was a higher level of genetic differentiation and a lower level of gene exchange between the populations. The neutrality test showed that the Fu’s Fs value was −2.992, and the Tajima’s D value was −0.65769. These were non-significant (p > 0.10), indicating that the Bgh populations followed the neutral model of evolution, and there were many loci with low-frequency alleles. The Fu and Li’ D (−0.89307) and Fu and Li’F (−0.97079) were also non-significant (p > 0.10), indicating that the Bgh populations were affected by selection and recombination.

4. Discussion

The differentiation of Bgh into different populations of distinct virulence and genetic composition has been reported globally [3,11,12,13,14,15]. In China, it has been found that there has been a recent change in the virulence frequency of Bgh to resistance genes and or gene combinations, which shows an overall trend toward higher pathotype diversity and virulence complexity [12,21,31]. The Bgh populations in central Europe also appear to exhibit a similar trend [32,33,34]. The isolates from Yunnan and Zhejiang had similar virulence profiles, but differed from those identified in Tibet [21]. In this research, disease samples were mainly collected from the cultivation areas of hull-less barley, except for Xichang and Chengdu in Sichuan. There were still 11 NILs with a virulence frequency of less than 10.0%. Eight out of the seventeen NILs expressing high resistance in the research of Wang et al. [21] were no longer highly resistant. This may be due to the differences in the sampling sites and possible breakdown of resistance in Tibet. From the result of the virulence assay of the Bgh population and disease incidence at different disease nurseries (unpublished data), resistance genotypes Mla6 + Mla14, Mla7 + Mlk, Mla7 + Mla (No3), Mla7 + M1 (Lg2), Mla9, Mla9 + Mlk and Mla12 + MlaEm2 are highly resistant to Bgh in Tibet and the surrounding areas. Inside Tibet, Mlra, Mla23, Mlk1, Mlg + Ml(CP) and the broad-spectrum resistance gene, mlo5, are also still resistant to powdery mildew. In order to obtain higher and more durable resistance, it would be valuable to pyramid such genes in combination with Mlo, together with the continuous monitoring of the changes in the virulence frequencies of different Bgh populations.
While barley is newly cultivated in Xichang and rare in Chengdu and Xining, the hull-less barley has long been introduced, first into the southern and eastern areas of the Qinghai-Tibet plateau, and has experienced a process of adaptation to the environment of varying altitudes [35]. Although it is in the recent two decades that powdery mildew has been reported to cause severe loss in Tibet [21], the incidence of the disease was documented as early as in 1987 [36]. Powdery mildew disease has not been observed in the disease nurseries set up at Daofu in Sichuan, Haibei and Xining in Qinghai. Severe powdery mildew was observed in the disease nursery at Hezuo in Gansu in 2019, but only a small number of colonies were observed in late September 2020, during 2020–2022. Obviously, the conclusion after the genetic analysis in this research, that the Bgh population had experienced a recent expansion, is congruent with the recent disease epidemic in the Tibet-Qinghai Plateau.
Both geographic isolation and long-distance migration had impacts on the population diversity of Blumeria graminis, and [21] has revealed the difference in Bgh genotypes between Tibet, Yunnnan, Sichuan in the southwest and Zhejiang in the southeast of China. The Bgh genotypes in this research revealed two centers of genetic diversity, i.e., Qamdo and Hezuo. The prevalence of H1 haplotypes at Qamdo and the locations surrounding Tibet suggested the possibility of the pathogen isolates being introduced into Tibet from other areas. When populations inside Tibet were looked at closely, Bgh seemed to have originated in the areas of Qamdo and expanded to the western areas of higher altitudes and lower temperature and humidity at Linzhi, Lhasa, Shannan and Rikaze. The introduction of barley into the Tibet-Qinghai plateau has resulted in the loss of most of the genetic diversity of the crop [35]. It is reasonable to postulate that the loss of resistance genes also happened during such a process due to the lack of disease pressure on the high and arid plateau. The results of virulence typing both in this research and in that of Wang et al. [21] revealed the weaker virulence of the Tibet Bgh population compared with that in the other areas of China. In contrast to the less diversified genetic structure than that of Bgh at Qamdo, the virulence become stronger in the pathogen population at the western sites in Tibet. This could be explained by the more intensive cultivation of barley with newly introduced resistance genes in commercial cultivars in the latter locations.

Author Contributions

Y.W., Y.P. and M.W. conceived the project and drafted the manuscript; Q.Z. and Z.X. collected the isolates and performed the pathogenicity test. All authors have read and agreed to the published version of the manuscript.

Funding

This work was financially supported by the Research on Key Technology of Ecological Disease and Pest Control for Barley (ZX202001ZY0043N), the Special Fund Project of the Central Government for Supporting the Development of Local Colleges and Universities (KY2022ZY-02), and the 2019 Innovation Plan for Postgraduates of Tibet Agricultural and Animal Husbandry University (YJS2019-22).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

All data supporting the findings of this study are available within the paper.

Acknowledgments

We would like to thank Ruiming Lin of the Institute of Plant Sciences, Chinese Academy of Agricultural Sciences for providing barley seeds, Qing Xue of the Department of Plant Pathology, Nanjing Agricultural University for providing technical assistance. We also thank the Agricultural Research Institute, Tibet Academy of Agriculture and Animal Husbandry Sciences, Changdu Academy of Agricultural Science, Shannan Academy of Agricultural Science, Rikaze Academy of Agricultural Science, Institute of Gannan Agricultural Science, Institute of Diqing Agricultural Science, and Qinghai Academy of Agricultural and Forestry Sciences for the kind assistance in the disease nurseries and with isolate sampling.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Figure A1. Geographical locations of the variation of barley powdery mildew in observational field. The base drawing is from the website of the Ministry of Natural Resources (http://www.mnr.gov.cn accessed on 9 November 2022), and the review drawing number is GS (2020) 4621.
Figure A1. Geographical locations of the variation of barley powdery mildew in observational field. The base drawing is from the website of the Ministry of Natural Resources (http://www.mnr.gov.cn accessed on 9 November 2022), and the review drawing number is GS (2020) 4621.
Jof 09 00363 g0a1
Table A1. The 30 barley genotypes with resistance allele(s) to Blumeria. graminis f. sp. hordei.
Table A1. The 30 barley genotypes with resistance allele(s) to Blumeria. graminis f. sp. hordei.
NumberVarietyResistance GenesOrigin
1P01Mla1, Mla(Al2)Denmark
2P02Mla3Denmark
3P03Mla6, Mla14Denmark
4P04AMla7,MlkDenmark
5P04BMla7, Mla(No3)Denmark
6P05/Denmark
7P06Mla7Ml(Lg2)Denmark
8P07Mla9, MlkDenmark
9P08AMla9, MlkDenmark
10P08BMla9Denmark
11P09Mla10, MlaDu2Denmark
12P10Mla12, MlaEm2Denmark
13P11Mla13, MlRu3, MlaRu4Denmark
14P12Mla22Denmark
15P13Mla23Denmark
16P14MlraDenmark
17P15Ml(Ru2)Denmark
18P16MlkDenmark
19P17Mlk1Denmark
20P18MlnnDenmark
21P19Mlp1Denmark
22P20Mlat, Mla8Denmark
23P21Mlg, Ml(CP)Denmark
24P22mlo5Denmark
25P23MlLaDenmark
26P24MlhDenmark
27P29/Denmark
28P30/Denmark
29P31/Denmark
30PallasMla8Sweden
Table A2. Numerical designations of the pathotypes showing virulence to resistance genes present in 21 differential varieties. Those with a frequency of >1 are marked in bold.
Table A2. Numerical designations of the pathotypes showing virulence to resistance genes present in 21 differential varieties. Those with a frequency of >1 are marked in bold.
Pathotypes
1.0.1.0.0.0.03.0.0.0.0.0.03.4.3.2.3.0.05.6.6.1.1.0.06.5.1.5.6.0.07.2.0.0.0.0.0
1.0.1.4.0.0.03.0.0.4.0.0.03.5.0.0.0.0.05.7.0.0.0.0.06.5.4.0.0.0.07.2.1.0.0.0.0
1.2.0.0.0.0.03.0.2.0.0.0.03.5.1.0.1.0.06.0.0.0.0.0.06.6.4.0.0.0.07.2.4.0.0.0.0
1.2.0.0.1.0.03.0.3.0.0.0.03.5.1.2.2.0.06.0.1.3.0.0.06.7.3.3.2.0.07.2.5.5.4.0.0
1.2.2.0.1.0.03.0.3.2.0.0.03.5.1.4.5.0.06.0.2.0.0.0.07.0.0.0.0.0.07.2.6.5.0.0.0
1.3.0.0.0.0.03.0.3.4.1.0.03.5.3.2.2.0.06.0.4.0.0.0.07.0.0.2.0.0.07.3.0.0.0.0.0
1.3.0.0.2.0.03.0.3.5.1.0.04.0.3.4.2.0.06.1.2.0.0.0.07.0.1.0.0.0.07.3.3.0.0.0.0
1.3.1.6.0.0.03.0.3.6.1.0.04.0.4.0.0.0.06.1.5.0.0.0.07.0.1.0.4.0.07.3.4.0.0.0.0
1.4.0.0.0.0.03.0.3.7.1.0.04.4.2.0.0.0.06.1.7.1.4.0.07.0.1.1.4.0.07.3.4.1.0.0.0
1.5.4.0.0.0.03.1.0.0.0.0.04.6.0.0.0.0.06.2.5.0.0.0.07.0.1.3.2.0.07.3.6.5.0.0.0
1.5.4.1.0.0.03.1.1.0.0.0.04.7.0.0.0.0.06.2.5.3.1.0.07.0.2.0.0.0.07.4.0.0.0.0.0
1.7.0.0.0.0.03.1.1.2.3.0.05.0.0.0.5.1.06.2.5.7.0.0.07.0.3.2.0.0.07.4.1.3.0.0.0
2.0.0.0.0.0.03.1.3.2.2.0.05.0.1.0.4.0.06.2.7.3.1.0.07.0.3.4.1.0.07.4.1.3.2.0.0
2.0.4.6.0.0.03.2.3.3.1.0.05.0.2.0.4.0.06.3.2.2.5.0.07.0.4.0.0.0.07.4.2.0.0.0.0
2.0.5.6.0.0.03.2.4.0.0.0.05.1.2.0.0.0.06.3.4.2.2.0.07.1.0.0.0.0.07.4.3.2.0.0.0
2.0.7.0.0.0.03.2.4.0.0.1.05.1.2.0.4.0.06.3.5.7.0.0.07.1.0.0.4.0.07.4.3.6.0.0.0
2.0.7.2.0.0.03.3.0.0.0.0.05.1.2.0.5.1.06.4.2.0.1.0.07.1.1.0.0.0.07.5.0.0.0.0.0
2.1.0.0.1.0.03.3.3.5.1.0.05.1.3.0.0.0.06.4.4.0.0.0.07.1.1.0.4.0.07.5.0.0.1.0.0
2.2.2.0.0.0.03.3.3.7.0.0.05.3.0.0.0.0.06.4.6.0.1.0.07.1.2.0.0.0.07.5.0.0.3.0.0
2.5.0.0.0.0.03.3.3.7.1.0.05.3.0.1.0.0.06.5.0.0.0.0.07.1.2.0.1.0.07.6.0.0.0.0.0
2.6.4.0.0.0.03.4.0.0.0.0.05.4.0.0.0.0.06.5.0.1.0.0.07.1.3.1.1.0.07.7.0.0.0.0.0
2.7.0.0.1.0.03.4.3.0.3.0.05.6.2.1.0.0.06.5.1.4.0.0.07.1.4.6.0.0.07.7.0.0.1.0.0

References

  1. IGC. 2019. Available online: https://www.igc.int/en/default.aspx (accessed on 5 August 2022).
  2. Çelik, O.; Karakaya, A. Genetic diversity of barley foliar fungal pathogens. Agronomy 2021, 11, 434. [Google Scholar] [CrossRef]
  3. Dreiseitl, A.; Yang, J. Powdery mildew resistance in a collection of Chinese barley varieties. Genet. Resour. Crop Evol. 2007, 54, 2. [Google Scholar] [CrossRef]
  4. Jia, X.; Sun, Z.; Li, X. Layout of barley production and its comparative Advantage in China. Agric. Outlook 2017, 13, 40–46. [Google Scholar]
  5. Zeng, Y.; Pu, X.; Zhang, J. Synthetic research and utilization on industrial development of barley in Southwestern China. J. Agric. Sci. Technol. 2013, 15, 48–56. [Google Scholar]
  6. Qiang, X.; Chi, D.; Feng, J. Status of the production and development of highland barley in the Qinghai-Tibet plateau region. Tibet Sci. Technol. 2008, 3, 11–17. [Google Scholar]
  7. Nyima, T.; Yu, D.; Tang, Y.; Bian, B.; Zeng, X. Great improvement of the yield of highland barley to ensure the safety of highland barley. Tibet Sci. Technol. 2013, 2, 5–7. [Google Scholar]
  8. Yuan, H.; Zeng, X.; Yang, Q.; Xu, Q.; Wang, Y.; Jabu, D.; Sang, Z.; Tashi, N. Gene coexpression network analysis combined with metabonomics reveals the resistance responses to powdery mildew in Tibetan hulless barley. Sci. Rep. 2018, 8, 14928. [Google Scholar] [CrossRef] [Green Version]
  9. Mathre, D.E. Compendium of Barley Disease, 2nd ed.; American Phytopathological Society: St. Paul, MN, USA, 1997. [Google Scholar]
  10. Yahyaoui, A.; Reinhold, M.; Scharen, A. Virulence spectrum in populations of the barley powdery mildew pathogen Erysiphe graminis f. sp. hordei in Tunisia and Morocco in 1992. Plant Pathol. 1997, 46, 139–146. [Google Scholar]
  11. Dreiseitl, A. Pathogenic divergence of central European and Australian populations of Blumeria graminis f. sp. hordei. Ann. Appl. Biol. 2014, 165, 364–372. [Google Scholar] [CrossRef]
  12. Zhu, J.; Zhou, Y.; Yang, J. Advances in the study of genetics of Blumeria graminis f. sp. hordei. J. Plant Prot. 2014, 41, 109–117. [Google Scholar]
  13. Chaure, P.; Gurr, S.; Spanu, P. Stable transformation of Erysiphe graminis an obligate biotrophic pathogen of barley. Nat. Biotechnol. 2000, 18, 205–207. [Google Scholar] [CrossRef]
  14. Jørgensen, J.; Jensen, H. Powdery mildew resistance in barley landrace material. I. Screening for resistance. Euphytica 1997, 97, 227–233. [Google Scholar] [CrossRef]
  15. Dreiseitl, A. Specific resistance of barley to powdery mildew, Its use and beyond. A concise critical review. Genes 2020, 11, 971. [Google Scholar] [CrossRef] [PubMed]
  16. Wan, Y.; Li, Y.; Gao, Q.; Qin, X.; Ma, X.; Liu, G. Climate change trend and its impact on yield of highland barley in Tibet, China. J. Agric. Resour. Environ. 2018, 35, 374–380. [Google Scholar]
  17. Liu, G. Effect of climate change on Tibetan agricultural production. Tibet J. Agric. Sci. 2019, 41, 49–54. [Google Scholar]
  18. Reinstädler, A.; Müller, J.; Czembor, J.; Piffanelli, P.; Panstruga, R. Novel induced mlo mutant alleles in combination with site-directed mutagenesis reveal functionally important domains in the heptahelical barley Mlo protein. BMC Plant Biol. 2010, 10, 31. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  19. Parks, R.; Carbone, I.; Murphy, J.; Cowger, C. Population genetic analysis of an eastern U.S. wheat powdery mildew population reveals geographic subdivision and recent common ancestry with U.K. and Israeli populations. Phytopathology 2009, 99, 840–849. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  20. Xu, Z.; Duan, X.; Zhou, Y. Population genetic analysis of Blumeria graminis f. sp. tritici in Qinghai Province, China. J. Integr. Agric. 2014, 13, 1952–1961. [Google Scholar] [CrossRef]
  21. Wang, Y.; Zhang, G.; Mu, W.; Lin, R. Virulence variability and genetic diversity in Blumeria graminis f. sp. hoedei in southeastern and southwestern China. Plant Dis. 2022, 8, 10. [Google Scholar]
  22. Jensen, H.; Christensen, E.; Jørgensen, J. Powdery Mildew Resistance Genes in 127 Northwest European Spring Barley Varieties. Plant Breed. 1992, 108, 210–228. [Google Scholar] [CrossRef]
  23. Gilmour, J. Octal notation for designating physiologic races of plant pathogens. Nature 1973, 242, 620. [Google Scholar] [CrossRef]
  24. Limpert, E.; Müller, K. Designation of pathotypes of plant pathogens. J. Phytopathol. 1994, 140, 346–358. [Google Scholar] [CrossRef]
  25. Cowger, C.; Parks, R.; Kosman, E. Structure and migration in U.S. Blumeria graminis f. sp. tritici populations. Phytopathology 2016, 106, 295–304. [Google Scholar] [PubMed] [Green Version]
  26. Schachtel, G.; 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] [Green Version]
  27. Earl, D.A.; vonHoldt, B.M. Structure Harvester: A website and program for visualizing STRUCTURE output and implementing the Evanno method. Conserv. Genet. Resour. 2012, 4, 359–361. [Google Scholar] [CrossRef]
  28. Librado, P.; Rozas, J. DnaSP v5: A software for comprehensive analysis of DNA polymorphism data. Bioinformatics 2009, 25, 1451–1452. [Google Scholar] [CrossRef] [Green Version]
  29. Excoffier, L.; Lischer, H. Arlequin suite ver 3.5: A new series of programs to perform population genetics analyses under Linux and Windows. Mol. Ecol. Resour. 2010, 10, 564–567. [Google Scholar] [CrossRef]
  30. Dupanloup, I.; Schneider, S.; Excoffier, L. A simulated annealing approach to define the genetic structure of populations. Mol. Ecol. 2002, 11, 2571–2581. [Google Scholar] [CrossRef]
  31. Wang, J.; Dreiseitl, A. Virulence and diversity of Blumeria graminis f. sp. hordei in East China. Eur. J. Plant Pathol. 2007, 117, 357–368. [Google Scholar]
  32. Dreiseitl, A. Great pathotype diversity and reduced virulence complexity in a central European population of Blumeria graminis f. sp. hordei in 2015–2017. Eur. J. Plant Pathol. 2019, 153, 801–811. [Google Scholar] [CrossRef]
  33. Dreiseitl, A. Genotype heterogeneity in accessions of a winter barley core collection assessed on postulated specific powdery mildew resistance genes. Agronomy 2021, 11, 513. [Google Scholar] [CrossRef]
  34. Helen, R.; Dreiseitl, A.; Sadiki, M.; Daniel, J. High diversity, low spatial structure and rapid pathotype evolution in Moroccan populations of Blumeria graminis f. sp. hordei. Eur. J. Plant Pathol. 2013, 136, 323–336. [Google Scholar]
  35. Zeng, X.; Guo, Y.; Xu, Q.; Mascher, M.; Guo, G.; Li, S.; Mao, L.; Liu, Q.; Xia, Z.; Zhou, J.; et al. Origin and evolution of qingke barley in Tibet. Nat. Commun. 2018, 9, 5433. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  36. Zhang, S. Diseases, Pests and Weeds in Tibet Agriculture; The Tibet People’s Publisher: Lhasa, China, 1987; pp. 393–394. [Google Scholar]
Figure 1. Dendrogram based on the cluster analysis of infection types of 21NILs of barley variety pallas from 229 Bgh isolates, using the QD-SAHN program in NTSYSpc 2.10s.
Figure 1. Dendrogram based on the cluster analysis of infection types of 21NILs of barley variety pallas from 229 Bgh isolates, using the QD-SAHN program in NTSYSpc 2.10s.
Jof 09 00363 g001
Figure 2. PCoA of the virulence variability of Bgh isolates from nine disease nurseries in Tibet, Sichuan, Yunnan, and Gansu.
Figure 2. PCoA of the virulence variability of Bgh isolates from nine disease nurseries in Tibet, Sichuan, Yunnan, and Gansu.
Jof 09 00363 g002
Figure 3. PopART Network of nine haplotypes of Blumeria. graminis f. sp. hordei.
Figure 3. PopART Network of nine haplotypes of Blumeria. graminis f. sp. hordei.
Jof 09 00363 g003
Figure 4. Distribution of the test isolates for haplotypes of Blumeria. graminis f. sp. hordei. (Note: The base drawing is from the website of the Ministry of Natural Resources (http://www.mnr.gov.cn accessed on 9 November 2022), the review drawing number is GS (2020) 4621, and the base drawing has not been modified).
Figure 4. Distribution of the test isolates for haplotypes of Blumeria. graminis f. sp. hordei. (Note: The base drawing is from the website of the Ministry of Natural Resources (http://www.mnr.gov.cn accessed on 9 November 2022), the review drawing number is GS (2020) 4621, and the base drawing has not been modified).
Jof 09 00363 g004
Figure 5. Genetic structure map inferred by STRUCTURE analysis consisting of Blumeria. graminis f. sp. hordei from nine populations.
Figure 5. Genetic structure map inferred by STRUCTURE analysis consisting of Blumeria. graminis f. sp. hordei from nine populations.
Jof 09 00363 g005
Table 1. Virulence frequencies of Blumeria. graminis f. sp. hordei from nine populations isolates on 30 barley NILs.
Table 1. Virulence frequencies of Blumeria. graminis f. sp. hordei from nine populations isolates on 30 barley NILs.
VarietyResistance Gene(s)Virulence Frequency (%)
QD aLZSNLSRKZHZCDXCDQTotal
P01Mla1, Mla(Al2)0000000000
P02Mla3035.7100016.6711.5425.0009.88
P03Mla6, Mla140000000000
P04AMla7, Mlk0000000000
P04BMla7, Mla(No3)0000000000
P05/021.4300016.677.6941.6758.3316.20
P06Mla7, Ml(Lg2)0000000000
P07Mla9, Mlk0000000000
P08AMla9, Mlk0004.55000000.50
P08BMla90000000000
P09Mla10, MlaDu260.0010066.6710070.0091.6776.9210010085.03
P10Mla12, MlaEm26.67000000000.74
P11Mla13, Ru3, Ru440.000054.543.70023.0854.0033.3323.18
P12Mla2253.3310083.3363.6477.7875.0065.3833.0058.3367.75
P-13Mla23021.4316.6713.64012.5019.234.1754.1715.76
P-14Mlra021.4300054.1765.3846.0066.6728.18
P15Ml(Ru2)60.0078.5783.3381.8244.4483.3376.9292.0091.6776.90
P16Mlk021.43007.41015.3833.006.259.27
P17Mlk107.14003.7062.50037.5075.0020.65
P18Mlnn60.0064.2971.0072.7348.1583.3384.6250.0041.6763.98
P19Mlp1075.00007.418.757.6971.0091.6729.06
P20Mlat, Mla873.3321.4366.6772.7314.8125.0023.0816.6716.6736.71
P21Mlg, Ml(CP)00018.18054.1723.0833.0058.3320.75
P22mlo50004.5522.2216.67021.0012.508.55
P23MlLa46.6771.4391.6750.0040.7429.1715.3854.0033.3348.04
P24Mlh33.3327.4366.6736.3659.2612.530.77058.3336.07
P29/0000000000
P30/73.3392.8691.6781.8285.1910092.3110010090.80
P31/0000000000
PallasMla8100100100100100100100100100100
a QD, Qamdo; LZ, Linzhi; SN, Shannan; LS, Lasha; RKZ, Rikaze; CD, Chengdu; XC, Xichang; DQ, Diqing.
Table 2. Virulence diversity parameters of nine Bgh populations.
Table 2. Virulence diversity parameters of nine Bgh populations.
QDLZSNLSRKZHZCDXCDQ
Number of isolates302824222724262424
Number of pathotypes131810181516212219
Average virulence complexity per isolate (Ci)3.7335.3935.0004.6823.3706.0424.8085.6257.625
Relative virulence complexity per isolate (Rci)0.1780.2570.2380.2230.1600.2880.2290.2680.363
Genotypic diversity (sh)0.7420.8460.6700.9030.7970.8360.9120.9640.891
Gene diversity (HS)0.1620.1650.1070.1730.1650.2090.2050.2610.256
Genetic diversity (KWm)0.2480.2480.1430.2640.2500.2980.2670.3850.385
Table 3. Nei’s (N) standard genetic distances and mean character difference (MCD) between nine Blumeria graminis f. sp. hordei populations.
Table 3. Nei’s (N) standard genetic distances and mean character difference (MCD) between nine Blumeria graminis f. sp. hordei populations.
PopulationsQDLZSNLSRKZHZCDXCDQ
QD-0.184 20.1030.1840.1110.2370.1440.180.282
LZ0.105 1-0.1470.1820.1570.1510.1620.1440.234
SN0.0440.08-0.0820.1190.2200.1700.2400.292
LS0.0170.1020.041-0.1390.1510.1620.1690.234
RKZ0.0480.0740.0570.062-0.2010.1560.1840.237
HZ0.140.070.1430.1190.114-0.1570.1270.151
CD0.0720.0850.0980.060.0710.083-0.1520.206
XC0.1010.0610.1340.0860.0930.0520.061-0.155
DQ0.1850.1190.1920.1580.1510.070.1120.057-
1 Nei’s (N) values: below diagonal; 2 Mean character difference (MCD) values: above diagonal.
Table 4. Distribution of Blumeria. graminis f. sp. hordei single-nucleotide polymorphism haplotypes.
Table 4. Distribution of Blumeria. graminis f. sp. hordei single-nucleotide polymorphism haplotypes.
HaplotypeAOX aPKAPPAQDLZSNLSRKZHZCDXCDQ
H1A bGAGACC6 20232224
H2 G 21 2
H3 G G 3
H4 A 1543454439 1
H5 G A 1
H6T 3
H7 G 1
H8 A 1
H9 A 4
Total 464345444224242624
a AOX = alternative oxidase, PKA = protein kinase A, and PPA = protein phosphatase type 2A. b Most frequent or consensus nucleotide for each variable position.
Table 5. Genetic diversity of multi-locus sequence haplotypes of Blumeria. graminis f. sp. hordei isolates from nine populations.
Table 5. Genetic diversity of multi-locus sequence haplotypes of Blumeria. graminis f. sp. hordei isolates from nine populations.
PopulationsNumber of IsolatesHaplotypes Diversityjof (Hd)Nucleotide Diversity (π)Average Number of Difference (K)
QD460.64850.000491.00735
LZ430.0000.000000.00000
SN450.0000.000000.00000
LS440.0000.000000.00000
RKZ450.11510.000010.23020
HZ240.34770.000170.37143
CD240.31170.000070.07407
XC260.14530.000140.31169
DQ240.00000.000000.00000
Total3210.57800.000340.74034
Table 6. Analysis of molecular variance (AMOVA) for Blumeria. graminis f. sp. hordei isolates from nine populations.
Table 6. Analysis of molecular variance (AMOVA) for Blumeria. graminis f. sp. hordei isolates from nine populations.
Source of Varianced.fSum of SquaresVariance ComponentsPercentage of Variation (%)Fixation Indices
Among groups145.2680.31240 Va54.68FCT = 0.54681
Among populations835.6260.13655 Vb23.90FST = 0.69714
Within populations31338.3020.12237 Vc21.42
Total322119.1950.57132
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

Wang, Y.; Zhuoma, Q.; Xu, Z.; Peng, Y.; Wang, M. Virulence and Genetic Types of Blumeria graminis f. sp. hordei in Tibet and Surrounding Areas. J. Fungi 2023, 9, 363. https://doi.org/10.3390/jof9030363

AMA Style

Wang Y, Zhuoma Q, Xu Z, Peng Y, Wang M. Virulence and Genetic Types of Blumeria graminis f. sp. hordei in Tibet and Surrounding Areas. Journal of Fungi. 2023; 9(3):363. https://doi.org/10.3390/jof9030363

Chicago/Turabian Style

Wang, Yunjing, Qucuo Zhuoma, Zhi Xu, Yunliang Peng, and Mu Wang. 2023. "Virulence and Genetic Types of Blumeria graminis f. sp. hordei in Tibet and Surrounding Areas" Journal of Fungi 9, no. 3: 363. https://doi.org/10.3390/jof9030363

APA Style

Wang, Y., Zhuoma, Q., Xu, Z., Peng, Y., & Wang, M. (2023). Virulence and Genetic Types of Blumeria graminis f. sp. hordei in Tibet and Surrounding Areas. Journal of Fungi, 9(3), 363. https://doi.org/10.3390/jof9030363

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