1. Introduction
About 50% of the world’s population uses rice as a staple food. Increasing rice production is the ultimate goal for meeting the demands of the ever-increasing world population. Although several rice varieties have been introduced over the years, high-yielding cultivars are, ironically, continuously threatened by biotic and abiotic stresses. The fungal pathogen that causes the blast disease (Magnaporthe oryzae) is the most serious worldwide source of disease for rice. It infects several cereal crops such as rice, wheat and rye grass [1,2]. Its infection can cause up to a 30% annual loss in terms of rice grain yield [3]. Studies have shown that the initial infection of the fungus is biotrophic; its hyphae invade and penetrate the plant cells [4]. In later stages, the fungus becomes necrotic and causes the death of the host cells [5]. Rice and its blast disease fungus represent the best pathosystem through which to understand the dynamics of the host–pathogen interaction [6]. Resistance to M. oryzae is controlled by resistance genes (R genes) in rice [7]. It is well known that the avirulence (AVR) genes in M. oryzae correspond to the R genes in rice [8]. Historically, the continuous selection pressure on hosts and their fungal pathogens has led to co-evolution. The selection of R genes and the mutation of pathogen-encoded AVR genes are considered the main driving forces in the arms race hypothesis. There is a continual consensus that the interaction between R and AVR genes lead to their diversification [9]. Previously, a Japonica rice variety, Lijiangxintuanheigu (LTH), was introduced in Yunnan Province [10,11]. LTH is a well-known susceptible variety to thousands of isolates of M. oryzae. Several R genes such as pi60, pi61, pi-g, pi36 and pi19 were cloned based on the genome sequence information obtained from LTH cultivars [12]. In addition, near-isogenic lines that contained single Pi genes were developed to assess the genetic structure and to monitor the new pathotypes of M. oryzae.
The study of R-AVR gene interactions suggests three different modes of host–pathogen interactions [13]. In the first mode, a single R gene directly interacts with its corresponding AVR gene [14]. The second type is characterized by the interaction of one R gene with two unrelated AVR genes [15], while the third type of interaction is represented by the fact that different R genes may interact with different AVR genes. The blast resistance developed by conventional breeding is overcome due to the high genetic diversity of AVR genes [16].
The Pi9 locus in indica rice shows a broad spectrum resistance to a wide range of M. oryzae strains [17]. It is predicted that Pi9 is expressed at high levels in the rice plants and contains a cluster of NBS-LRR [18]. The Pi9-R-gene-based resistance strategy is effective against M. oryzae due to the presence of corresponding AVR-Pi9 genes, which can trigger the plant defense mechanism. AVR-Pi9 is a small gene that codes for only 99 amino acids and is located on chromosome 7 of the M. oryzae genome [19]. The corresponding Pi9 is thought to have been introduced from Oryza minuta [19]. The Pi9 monogenic lines, developed at the International Rice Research Institute (IRRI), show resistance to most rice blast strains, except IE1k [16].
The higher genetic variability in M. oryzae avirulent genes is considered to be an important factor in defeating the R genes encoded by plants. Recently, a diversity analysis of AVR genes from 60 rice blast strains was conducted in Thailand. Only a single amino acid change was detected in the AVR-Pi9 gene. The data suggested that AVR-Pi9 has relatively less genetic diversity compared to other AVR genes [20].
The present study was aimed to demonstrate the genetic diversity and distribution of the AVR-Pi9 gene in M.-oryzae-infected rice and non-rice plants in Yunnan province, China. We collected 328 samples of the blast fungus, and the AVR-Pi9 gene was amplified using specific primers. A sequence comparison indicated the presence of seven different mutants. In this study, the pathogenicity and dynamics of rice blast fungus are discussed.
2. Materials and Methods
Source, isolation and culture of fungi: About 326 isolates of M. oryzae and 2 non-rice fungi were collected from Yunnan Province, China. The samples were collected based on the manifested disease phenotypes in rice plants [21,22,23]. The in planta fungal isolates were collected from the infected rice plants and observed under a light microscope. The isolates were transferred to filter paper at freezing temperatures until they were grown on oatmeal-supplemented agar. All the spores were grown at room temperature under controlled lighting [24]. After confirmation of M. oryzae, the fungal spores were isolated and cultured in the dark at 25 °C for 7 days, as was previously described in [25,26].
DNA isolation, PCR amplification and sequencing: To produce mycelium, the fungal isolates were grown in the dark at 25 °C for 6–8 days. The total genomic DNA was isolated via grown mycelium and plant leaves while using the cetyl trimethylammonium bromide (CTAB) method, as was described previously in [27,28]. The isolated DNA was run on 1% agarose gel and quantified using a spectrophotometer (Nanodrop 2000). Approximately, 50 ng/µL of DNA was used in the PCR amplification reaction. Each PCR was performed in triplicates. Briefly, full-length AVR-Pi9 specific primers were used in a standard 25 µL PCR reaction. AVR-Pi9F atgcagttctctcagatcctcacc and AVR-Pi9R ctaccagtgcgtcttttcgacttg [19,29] were synthesized by Tiangen Biotech Co., Ltd., Beijing, China. All the PCR products were resolved on 1% agarose gel, and the amplification of ~1500 nucleotides (nts) was confirmed (DL2000 DNA Ladder; Tiangen Biotech Co., Ltd., Beijing, China). All the PCR products were sequenced by the specific primers (Shanghai Life Technologies Biotechnology Co., Ltd., Beijing, China).
Sequence assembly and nucleotide diversity analysis: The AVR-Pi9 sequence data were annotated in the SeqMan module of Lasergene-7 (DNA-STAR, Madison, WI, USA). The accurate small reads of ~800 nts were assembled into a single consensus file. To validate the sequencing, the sequences obtained from different isolates were compared, via BLASTn, with the reference sequences already available in the NCBI GenBank. The AVR sequence files of 328 (326 rice and 2 non-rice) isolates were aligned using Clustal-X [30]. For the phylogenetic analysis, the aligned sequence files were exported to MEGA 11.0 software [31]. The phylogenetic tree was constructed using the maximum likelihood method of tree construction with 1000 bootstrap replications. The nucleotide polymorphism, haplotype diversity, (π and Hd, respectively), and DNA haplotype polymorphic sites were calculated using the DnaSP-6 software package [32,33,34]. Population neutrality and selection pressure tests were calculated using Tajima’s D test [35,36]. The DNA sequence-based TCS software was used (
Disease pathogenicity assay: The pathogenicity test was performed on the resistant genotype (IRBL9-W) that harbors the Pi9 gene in the background of LTH. The LTH lines were used as a susceptible control [37]. Briefly, the conidia were randomly collected and maintained on a culture medium, as was described previously in [29,38]. The rice seedlings were grown for 21 days, and the spore suspension that contained conidia was sprayed on the leaves. The spore suspensions of these isolates (105 spores/mL) were sprayed on 3–4 leaf stage seedlings. After inoculation, the plants were transferred to a dark room with more than 90% humidity at 25 °C for 24 h. After 7 days, the plants were photographed and disease symptoms were scored according to a previously described scale [24].
3. Results
3.1. Distribution of M. oryzae Isolates in Yunnan Province, China
A total of 328 samples were collected from the Yunnan province of China during the growing seasons of 2010 to 2020 (Table 1). The main sampling area included the central, northern, southern and eastern areas of the province. In the PCR reaction with specific primers of AVR-Pi9 genes, we obtained a positive amplification in 328 samples (Table 1). To exclude possible contamination, each PCR was repeated three times. The sequencing results also confirmed the presence of M. oryzae in all the samples collected.
The pathogenicity assay indicated that all the isolates were avirulent to the monogenic IRBL9-W line. The frequency of the virulent strains was 100% in the field, with the exception of the western region—where it was observed as 96%. Overall, the frequency of the AVR-Pi9 gene was recorded as 98% (Table 1). The nucleotide alignment of the AVR-Pi9 gene revealed the presence of seven different haplotypes (H1–H7). The nucleotide variation analysis indicated that minor mutations were scattered in the promoter, coding and intronic regions. These findings suggested that although Pi9 is a well-known source of resistance to M. oryzae, virulent strains are also circulating in the region, which may lead to resistance breakdown on a large scale (Table 2).
3.2. Identification of AVR-Pi9 Haplotypes
The AVR-Pi9 gene consisting of 1436 nucleotides was amplified from the 328 isolates by using gene specific primers. Amplified PCR products from 328 isolates (including 326 isolates collected from rice and 2 isolates collected from non-rice) were sequenced and assembled. The gene sequence was compared with the KM004023.1 reference sequence (Table 3). After the sequence annotation and assembly analysis, eight different haplotypes (H1–H8) were identified (Table 3). H1 was identified as the major haplotype with the highest percentage (66.5%) among the other seven haplotypes (Table 3, Figure S1). Sequence analysis revealed the highest sequence identity of the CDS region in the reference sequence. The H2 haplotype was characterized by a deletion of 15 nucleotides at the variant locus. Interestingly, H2, H3 and H4 haplotypes were characterized as deletion mutants in the promotor region of the variant locus. In addition, the H3 haplotype has a substitution (T-nucleotide) when compared to the reference strain. The H5, H6 and H7 haplotypes were found at very low frequencies. However, these haplotypes showed nucleotide substitutions in the CDS region (Table 3). The sequence variation of these seven haplotypes was also evident from the phylogenetic tree (Figure 1). The maximum likelihood phylogenetic tree showed three major groups. One group contained the reference sequence, H1, H2 and H3 haplotypes. Haplotype H4 diverged due to sequence variation in the promoter region (Table 3). The haplotypes H5, H6 and H7 were clustered together due to high-nucleotide identity. However, there were differences in the positions of the nucleotide substitution sites among these three isolates. H8 haplotype contained two isolates collected from Eleusine coracana and Eleusine indica, respectively. The diversity analysis of the haplotypes suggests that AVR-Pi9 may be under either a strong population expansion or a purifying selection due to the introduction of resistant rice germplasms. The results of the polymorphic analysis show that AVR-Pi9 may be under a strong positive selection. The sliding window for the distribution of Ka/Ks values suggests the possibility of potential purifying positive selections (Figure 2). Indeed, the sites of positive selection were clearly present in the exon region of the AVR-Pi9 gene. The haplotype diversity and polymorphic analysis (π) also confirmed the distribution of variation in the CDS region (Figure 2).
3.3. Genetic Variation of the AVR-Pi9 Gene in Non-Rice Hosts
In this study, the AVR-Pi9 gene from three non-rice hosts, including perennial rye grass and wheat (Table 4), is investigated. The AVR-Pi9 sequence from Setaria viridis was also included in this study and was obtained from GenBank. Notably, the non-rice isolates had high genetic similarity to the previously known sources of the AVR-Pi9 gene. Sequence comparison showed that the isolates from ryegrass and Setaria were 98% identical to the reference sequence (AVR-Pi9-KM004023), while the two isolates from wheat, namely PY5033 and PY6045, showed a 97% sequence identity with the reference strain. Close examination revealed the deletions in the nucleotide sequence of the non-rice hosts.
3.4. Selection Pressure of AVR-Pi9 in M. oryzae
To determine the natural selection pressure of the AVR-Pi9 in the M. oryzae from Yunnan Province, the Tajima’s neutrality of the AVR-Pi9 in M. oryzae was tested based on 328 AVR-Pi9 DNA sequences, and the Tajima’s D was −1.85719 (Table 5). The lower value of Tajima’s D indicates the high number of alleles in the rice blast population in Yunnan Province, China. These results also suggest that AVR-Pi9 may be under either strong population expansion or purifying selection due to the introduction of resistant rice germplasms.
3.5. The Phylogenetic Relationship of AVR-Pi9 Haplotypes
To understand the relationship of isolates from the Yunnan province, the phylogenetic tree was constructed from the amplified fragment of AVR-Pi9. The phylogenetic tree that was based on the AVR-Pi9 sequences of 328-blast isolates showed that the majority of them were closely clustered, except H4, H6 and H7 (Figure 3). The close clustering of isolates confirmed their low level of genetic diversity. However, smaller but significant changes in the H4, H6 and H7 isolates resulted in their unique position in the phylogenetic tree. Among the two clusters (A and B), Cluster A included only the isolates collected from rice, while Cluster B contained the isolates from the Eleusine coracana and Eleusine indica of the H8 haplotype.
3.6. Pathogenicity Assay
The previous data showed that the rice varieties harboring the Pi9 gene showed broad-spectrum resistance to M. oryzae. We tested the disease response of the existing rice blast fungus on the landraces collected from Yunnan Province, China. The landraces that were tested against YN700 (H1 haplotype, identical to haplotype KM004023.1) were resistant and showed avirulence. This test confirmed the existing hypothesis that monogenic lines carrying the Pi9 gene are resistant to blast fungus. To confirm the resistance and susceptibility, we inoculated the landraces with a rice blast fungus containing the AVR-Pi9 allele. In the experiment where the haplotype KM004023.1 was used, none of the landraces developed a positive disease response (Figure 4). All the previously tested landraces consistently showed resistance to the existing haplotypes. Nonetheless, the inoculation of CYN201 (H8 isolated from rye grass) on the same ten landraces resulted in a sever disease reaction. All landraces showed sever blast symptoms and eventually died. These data led to the conclusion that the newly identified haplotypes may be lethal to rice production in Yunnan Province, China.
4. Discussion
The introgression of R genes into rice germplasm was an important effort to conduct in order to regain the broad-spectrum resistance of Magnaporthe oryzae [16]. Several monogenic and digenic lines have been developed worldwide to counter the damage caused by rice blast fungus [39]. Despite the large scale high-throughput screening and successful cloning of AVR genes, a complete understanding of them remains a challenge. Although comparative genomics data of field isolates are available, the exact diversity of the AVR locus is still unknown. Usually, after the outbreak of the disease epidemic, the resistant plants become susceptible and, in the following years, a new combination of natural variants in the blast fungus is observed. In Chinese monogenic lines, the resistance to rice blast fungus was introduced in 2000. In fact, five resistant genes were introduced, via the conventional breeding techniques, into the blast susceptible Japonica rice variety named Lijiangxin-tuan-heigu (LTH) [40,41]. At the protein level, Pi9 is a member of the NBS-LRR genes that show broad-spectrum resistance to different rice blast strains.
In this study, we analyzed 326 M. oryzae isolates from the rice plants grown in Yunnan Province, China. We identified small variations in the AVR gene and its regulatory region. According to the gene-for-gene interaction model, the genetic variation in AVR may lead to a disruption of its corresponding R gene. The fungus-encoded AVR genes are targets of the plant-encoded R genes. It is clear that rice varieties harboring the Pi9 gene show immunity to different rice blast strains. It is also well established that there is a strong correlation between Pi9 and AvrPi9. Mechanistically, Pi9 can recognize AvrPi9, which results in a broad-spectrum resistance to diverse rice blast isolates. The AVR genes evolve rapidly, resulting in the breakdown of their corresponding resistance genes. In a recent study on the natural variation of AVR genes in the Sichuan Basin of China, there was a loss of natural resistance against M. oryzae [42]. Therefore, it is believed that closed geographical monitoring of rice lines is important to control the spread of a virulent blast strain. In this study, we analyzed the spread of disease in the Yunnan province of China, where a monogenic resistant line containing the Pi9 (IRBL9-W) R gene was cultivated on a large scale. We also included three non-rice hosts (wheat, ryegrass and Setaria viridis) to determine the natural variation of the AVR-Pi9 genes. The complete sequencing of 328 isolates of M. oryzae was obtained. Interestingly, through sequence analysis, we identified the insertions, deletions and substitution mutations in the different isolates. These data led to the identification of seven unique haplotypes that were not previously known. The H1 haplotype was mostly identical to the existing blast strain sequences. The higher frequency of the H1 haplotype was due to its genetic identity with the existing strain in the coding region of the AVR-Pi9 gene.
The pathogenicity assay suggests that the resistance phenomenon is invalidated by the most prevalent H1 haplotype of AVR-Pi9. The mutations in the promoter region of the AVR-Pi9 gene are intriguing. In the previous studies on AVR genes, it was reported that transposable element insertion sites are located in the promoter region [43,44]. The evolution of mutant magnaporthe races is due to the loss of avirulence genes, which have important interactions with the resistance genes. Even a small point mutation can cause the loss of avirulence. Previous study showed that the deletion of avirulence genes is a common phenomenon in diverse natural populations [45].
It is well established that the nucleotide insertions in AVR genes result in the loss of resistance [46]. Previous studies conducted on a large collection of different cultivars in China suggest that the Pi9 locus in the rice genome confers broad spectrum resistance due to the presence of multiple novel alleles [47]. The allelic variation of the Pi9 gene in the rice genome may be a strong reason for the presence of AVR-Pi9 haplotypes [13]. It is possible that for each allele of Pi9, there is an evolution of AVR-Pi9. Previous studies suggest that there is a direct interaction from the AVR-Pi9 gene with the LRR domain of its cognate Pi9 gene [14]. Currently, more than one hundred R genes have been characterized, and several of them have been incorporated into rice varieties worldwide [13]. Most of these R genes contain nucleotide-binding site leucine-rich repeats (NBS-LRRs), which can potentially interact with pathogen AVR proteins [48]. However, the R genes Pi-d2 and Pi21 contain B lectin receptor kinase and proline-rich metal binding proteins, respectively, instead of NBS-LRRs [45]. In the future, it will be interesting to investigate whether the different haplotypes of AVR-Pi9 are capable of interacting directly with their cognate Pi9 gene.
5. Conclusions
We identified seven new haplotypes of the major AVR-Pi9 gene from the field-grown rice population. Our findings suggest that there may be more variations in the AVR genes than expected. The diversity in AVR sequences can be alarming as the new haplotypes of the fungus can escape the resistance provided by the host genes. Due to the continuous arms race, the evolution of AVR genes may outpace the evolutionary rate of host resistance mechanisms. Therefore, in addition to genetic resistance, the management of fungal pathogens should also be practiced for sustainable crop production.
Conceptualization, J.L. and C.L.; methodology, Q.W.; software, J.L. validation, L.L., Q.W. and Z.S.; formal analysis, Q.W.; investigation, Z.G.; resources, Z.S.; data curation, J.L.; writing—original draft preparation, L.L.; writing—review and editing, J.L.; visualization, Q.W.; supervision, Z.G.; project administration, J.L.; funding acquisition, J.L. All authors have read and agreed to the published version of the manuscript.
Not applicable.
Not applicable.
Not applicable.
We warmly acknowledge Lihua Wang from the Biotechnology and Germplasm Resources institute, Yunnan Academy of Agricultural Sciences of China, as well as Jianbin Zhang from the Chinese Academy of Tropical Agricultural Sciences for providing the plant materials used for the experiments.
The authors declare no conflict of interest.
Footnotes
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.
Enlarge this image.Figure 1. Maximum likelihood phylogenetic tree of the seven haplotypes identified from rice in the Yunnan province of China. The seven haplotypes can be divided into three different clusters as per their sequence variation. The H1 haplotype showed high genetic identity with the reference strain (KM004023.1).
Enlarge this image.Figure 2. The diversification of AVR-Pi9 in avirulent isolates. The distribution of variation in the AVR-Pi9 alleles was analyzed via the sliding window method. The X-axis shows the distribution of variation within the full region, including the promoter region and exons of AVR-Pi9. The lower pane indicates the corresponding schematic presentation of the promoter region in white color, and two exons in red color are for AVR-Pi9. Window length: 10. Step size: 2. The π value corresponds with the level of variation at each site because it is the sum of the pair-wise differences divided by the number of pairs within the population.
Enlarge this image.Figure 3. The phylogenetic tree constructed with the nucleotide sequences of the AVR-Pi9 genes from 328 Magnaporthe oryzae isolates. The maximum likelihood methods of MEGA V5.10. The ID number of KM004023 (GenBank Accession No. KM004023) for AVR-Pi9 was obtained from GenBank.
Enlarge this image.Figure 4. Disease reaction of the identification isolates of YN700 and YNC201 on the landrace accessions. YN700 (which belongs to the H1 haplotype) and YNC201 are avirulent and virulent to IRBL9-W (holding Pi9), respectively. Note: the numbers 14, 22, 62, 68, 69, 132, 143, 148, 241, 249 on top of the photo show the landrace varieties of Qi-He-Da-Hong-Gu, Lao-Leng-Gu, Da-Huang-Pi-Nuo, Bei-Zi-Nuo, Zao-Hong-Gu, Xiao-Hua-Gu, Da-Bo-Gu, Hao-Bu-Ka, Qie-Xie-Gu and San-Bai-Bang, respectively.
Distribution of the AVR-Pi9 genes and the avirulence isolates of the M. oryzae collected from Yunnan, China.
| Host | Locations/Years | No. of Isolates | PCR Detection | Pathogenicity Assay a | ||
|---|---|---|---|---|---|---|
| No. of Isolates with AVR-Pi9 | Frequency (%) | No. of Avirulence Isolates | Frequency (%) | |||
| Rice | Northeastern/2012 | 44 | 44 | 100.0 | 44 | 100.0 |
| Central/2013 | 88 | 88 | 100.0 | 88 | 100.0 | |
| Southeastern/2013 | 21 | 21 | 100.0 | 21 | 100.0 | |
| Western/2013 | 127 | 127 | 100.0 | 122 | 96.1 | |
| Northwestern/2014 | 21 | 21 | 100.0 | 21 | 100.0 | |
| Southwestern/2014 | 25 | 25 | 100.0 | 25 | 100.0 | |
| Total | 326 | 326 | 100.0 | 321 | 98.5 | |
| Non-rice | Western/2014 | 1 | 1 | 100.0 | 1 | 100.0 |
| Central/2014 | 1 | 1 | 100.0 | 1 | 100.0 | |
| Total | 2 | 2 | 100.0 | 2 | 100.0 | |
a Pathogenicity test was performed on the resistant genotype (IRBL9-W) as described previously (35–38).
Distribution of the AVR-Pi9 haplotype in the different rice regions of Yunnan Province.
| Haplotype | No. of Isolates | Frequency (%) | No. of Isolates and Frequency (%) in Each Region | Production a | ||||||
|---|---|---|---|---|---|---|---|---|---|---|
| Northeastern | Central | Southeastern | Western | Northwestern | Southwestern | XI | GJ | |||
| H01 | 218 | 66.9 | 42 (95.5) | 48 (54.5) | 12 (57.1) | 95 (74.8) | 8 (38.1) | 13 (52.0) | 120 (69.4) | 98 (64.1) |
| H02 | 87 | 26.7 | 1 (2.3) | 30 (34.1) | 9 (42.9) | 25 (19.7) | 13 (61.9) | 9 (36.0) | 43 (24.9) | 44 (28.8) |
| H03 | 12 | 3.7 | 0 | 7 (8.0) | 0 | 5 (3.9) | 0 | 0 | 5 (2.9) | 7 (4.6) |
| H04 | 2 | 0.6 | 0 | 0 | 0 | 0 | 0 | 2 (8.0) | 2 (1.2) | 0 |
| H05 | 3 | 0.9 | 0 | 0 | 0 | 2 (1.6) | 0 | 1 (4.0) | 3 (1.7) | 0 |
| H06 | 1 | 0.3 | 1 (2.3) | 0 | 0 | 0 | 0 | 0 | 0 | 1 (0.7) |
| H07 | 3 | 0.9 | 0 | 3 (3.4) | 0 | 0 | 0 | 0 | 0 | 3 (2.0) |
| Total | 326 | 100 | 44 | 88 | 21 | 127 | 21 | 25 | 173 | 153 |
| No. of haplotypes | 3 | 4 | 2 | 4 | 2 | 4 | 5 | 5 | ||
| Index of diversity b | 0.088 | 0.579 | 0.490 | 0.400 | 0.472 | 0.592 | 0.456 | 0.505 | ||
a XI and GJ indicates Xian/Indica and Geng/Japonic, respectively. b Diversity index was calculated as the frequency of the haplotype types in the M. oryzae population following Fontaine’s method. Diversity index = (1 − ∑ni=1pi2) (where pi is the frequency of the haplotype i in a population).
Haplotypes of the AVR-Pi9 loci in the rice blast fungus of Yunnan, China.
| Haplotype | No. of Isolates | Frequency (%) | Variant Locus a | |||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|---|
| Promoter Region | CDS | Intron | ||||||||||
| Between 168–169 | Between 812–827 | 912 | 1139 | 1169 | 1310 | Between 1342–1343 | 1345 | 1374 | 1376 | |||
| KM004023.1 | - | CTCCTACACTGGGGCT | T | C | C | C | - | G | G | C | ||
| H1 | 218 | 66.5 | . | . | . | . | . | . | - | . | . | . |
| H2 | 87 | 26.5 | . | - | . | . | . | . | - | . | . | . |
| H3 | 12 | 3.7 | T | - | . | . | . | . | - | . | . | . |
| H4 | 1 | 0.3 | . | - | C | . | . | . | - | . | . | . |
| H5 | 4 | 1.2 | T | . | . | . | . | . | - | . | . | . |
| H6 | 1 | 0.3 | . | . | . | . | T | . | - | . | . | . |
| H7 | 3 | 0.9 | . | . | . | . | . | . | - | . | . | T |
| H8 * | 2 | 0.6 | . | . | . | T | . | T | GCCCTGTACAATGCTTTTT | T | T | . |
| 328 | 100 | |||||||||||
a “-” indicates absent or deletion. “.” indicates nucleotide identity with KM004023.1 (accession ID) from GenBank. * denotes the haplotypes isolated from non-rice hosts.
The nucleotide identity of the AVR-Pi9 isolated from the three non-rice hosts.
| Host and Isolate | ||||||||
|---|---|---|---|---|---|---|---|---|
| Rice | Rice | Perennial Ryegrass | Perennial Ryegrass | Wheat | Wheat | Setaria viridis (L.) Beauv | ||
| Query | Size (bp) | P131 | PY34 | PgKY | PGPA | PY5033 | PY6045 | US71 |
| AVR-Pi9_KM004023.1 | 3496 letters | 2301/2302 (99%) | 3490/3502 (99%), Gaps = 9/3502 (0%) | 1242/1258 (98%), Gaps = 10/1258 (0%) | 1242/1258 (98%), Gaps = 10/1258 (0%) | 1237/1274 (97%), Gaps = 27/1274 (2%) | 1237/1274 (97%), Gaps = 27/1274 (2%) | 1372/1393 (98%), Gaps = 19/1393 (1%) |
Results from the Tajima’s neutrality test conducted for the AVR-Pi9 in M. oryzae a.
| m | S | π | D |
|---|---|---|---|
| 328 | 7 | 0.00006 | −1.85719 * |
a The analysis involved 328 nucleotide sequences of AVR-Pi9. S indicates the number of segregating sites, π indicates nucleotide diversity, and D is the Tajima test statistic. * statistical significance and p < 0.05.
Supplementary Materials
The following supporting information can be downloaded at:
References
1. Coelho, M.A.D.O.; Torres, G.A.M.; Cecon, P.R.; Santana, F.M. Sowing Date Reduces the Incidence of Wheat Blast Disease. Pesqui. Agropecu. Bras.; 2016; 51, pp. 631-637. [DOI: https://dx.doi.org/10.1590/S0100-204X2016000500025]
2. Farman, M.L. Pyricularia Grisea Isolates Causing Gray Leaf Spot on Perennial Ryegrass (Lolium Perenne) in the United States: Relationship to P. Grisea Isolates from Other Host Plants. Phytopathology; 2002; 92, pp. 245-254. [DOI: https://dx.doi.org/10.1094/PHYTO.2002.92.3.245] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/18943995]
3. Skamnioti, P.; Gurr, S.J. Against the Grain: Safeguarding Rice from Rice Blast Disease. Trends Biotechnol.; 2009; 27, pp. 141-150. [DOI: https://dx.doi.org/10.1016/j.tibtech.2008.12.002] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/19187990]
4. Heath, M.C.; Valent, B.; Howard, R.J.; Chumley, F.G. Interactions of Two Strains of Magnaporthe Grisea with Rice, Goosegrass, and Weeping Lovegrass. Can. J. Bot.; 1990; 68, pp. 1627-1637. [DOI: https://dx.doi.org/10.1139/b90-209]
5. Singh, J.; Gupta, S.K.; Devanna, B.N.; Singh, S.; Upadhyay, A.; Sharma, T.R. Blast Resistance Gene Pi54 Over-Expressed in Rice to Understand Its Cellular and Sub-Cellular Localization and Response to Different Pathogens. Sci. Rep.; 2020; 10, 5243. [DOI: https://dx.doi.org/10.1038/s41598-020-59027-x]
6. Coll, N.S.; Epple, P.; Dangl, J.L. Programmed Cell Death in the Plant Immune System. Cell Death Differ.; 2011; 18, pp. 1247-1256. [DOI: https://dx.doi.org/10.1038/cdd.2011.37]
7. Xing, J.; Jia, Y.; Correll, J.C.; Lee, F.N.; Cartwright, R.; Cao, M.; Yuan, L. Analysis of Genetic and Molecular Identity among Field Isolates of the Rice Blast Fungus with an International Differential System, Rep-PCR, and DNA Sequencing. Plant Dis.; 2013; 97, pp. 491-495. [DOI: https://dx.doi.org/10.1094/PDIS-04-12-0344-RE]
8. Flor, H.H. Current Status of the Gene-for-Gene Concept. Annu. Rev. Phytopathol.; 1971; 9, 275. [DOI: https://dx.doi.org/10.1146/annurev.py.09.090171.001423]
9. Dodds, P.N.; Lawrence, G.J.; Catanzariti, A.-M.; Teh, T.; Wang, C.-I.; Ayliffe, M.A.; Kobe, B.; Ellis, J.G. Direct Protein Interaction Underlies Gene-for-Gene Specificity and Coevolution of the Flax Resistance Genes and Flax Rust Avirulence Genes. Proc. Natl. Acad. Sci. USA; 2006; 103, pp. 8888-8893. [DOI: https://dx.doi.org/10.1073/pnas.0602577103]
10. Telebanco-Yanoria, M.J.; Koide, Y.; Fukuta, Y.; Imbe, T.; Tsunematsu, H.; Kato, H.; Ebron, L.A.; Nguyen, T.M.N.; Kobayashi, N. A Set of Near-Isogenic Lines of Indica-Type Rice Variety CO 39 as Differential Varieties for Blast Resistance. Mol. Breed.; 2011; 27, pp. 357-373. [DOI: https://dx.doi.org/10.1007/s11032-010-9437-x]
11. Fang, W.W.; Liu, C.C.; Zhang, H.W.; Xu, H.; Zhou, S.; Fang, K.X.; Peng, Y.L.; Zhao, W.S. Selection of Differential Isolates of Magnaporthe Oryzae for Postulation of Blast Resistance Genes. Phytopathology; 2018; 108, pp. 878-884. [DOI: https://dx.doi.org/10.1094/PHYTO-09-17-0333-R]
12. Yang, L.; Zhao, M.; Sha, G.; Sun, Q.; Gong, Q.; Yang, Q.; Xie, K.; Yuan, M.; Mortimer, J.C.; Xie, W. The Genome of the Rice Variety LTH Provides Insight into Its Universal Susceptibility Mechanism to Worldwide Rice Blast Fungal Strains. Comput. Struct. Biotechnol. J.; 2022; 20, pp. 1012-1026. [DOI: https://dx.doi.org/10.1016/j.csbj.2022.01.030]
13. Wang, B.; Ebbole, D.J.; Wang, Z. The Arms Race between Magnaporthe Oryzae and Rice: Diversity and Interaction of Avr and R Genes. J. Integr. Agric.; 2017; 16, pp. 2746-2760. [DOI: https://dx.doi.org/10.1016/S2095-3119(17)61746-5]
14. Jia, Y.; McAdams, S.A.; Bryan, G.T.; Hershey, H.P.; Valent, B. Direct Interaction of Resistance Gene and Avirulence Gene Products Confers Rice Blast Resistance. EMBO J.; 2000; 19, pp. 4004-4014. [DOI: https://dx.doi.org/10.1093/emboj/19.15.4004]
15. Cesari, S.; Thilliez, G.; Ribot, C.; Chalvon, V.; Michel, C.; Jauneau, A.; Rivas, S.; Alaux, L.; Kanzaki, H.; Okuyama, Y. The Rice Resistance Protein Pair RGA4/RGA5 Recognizes the Magnaporthe Oryzae Effectors AVR-Pia and AVR1-CO39 by Direct Binding. Plant Cell; 2013; 25, pp. 1463-1481. [DOI: https://dx.doi.org/10.1105/tpc.112.107201] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/23548743]
16. Wang, J.C.; Correll, J.C.; Jia, Y. Characterization of Rice Blast Resistance Genes in Rice Germplasm with Monogenic Lines and Pathogenicity Assays. Crop Prot.; 2015; 72, pp. 132-138. [DOI: https://dx.doi.org/10.1016/j.cropro.2015.03.014]
17. Liu, G.; Lu, G.; Zeng, L.; Wang, G.-L. Two Broad-Spectrum Blast Resistance Genes, Pi9 (t) and Pi2 (t), Are Physically Linked on Rice Chromosome 6. Mol. Genet. Genom.; 2002; 267, pp. 472-480. [DOI: https://dx.doi.org/10.1007/s00438-002-0677-2] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/12111554]
18. Ou, S.; Liu, G.; Zhou, B. The Broad-Spectrum Blast Resistance Gene Pi9 Encodes a Nucleotide-Binding Site Leucine-Rich Repeat Protein and Isa Member of a Multigene Family in Rice. Genetics; 2006; 172, pp. 1901-1914.
19. Wu, J.; Kou, Y.; Bao, J.; Li, Y.; Tang, M.; Zhu, X.; Ponaya, A.; Xiao, G.; Li, J.; Li, C. Comparative Genomics Identifies the Magnaporthe Oryzae Avirulence Effector AvrPi9 That Triggers Pi9-Mediated Blast Resistance in Rice. New Phytol.; 2015; 206, pp. 1463-1475. [DOI: https://dx.doi.org/10.1111/nph.13310]
20. Sutthiphai, T.; Damchuay, K.; Neupane, R.C.; Longya, A.; Sriwongchai, T.; Songkumarn, P.; Parinthawong, N.; Darwell, K.; Jantasuriyarat, C. Genetic Variation of Avirulence Genes (AVR-Pi9, AVR-Pik, AVR-Pita1) and Genetic Diversity of Rice Blast Fungus, Pyricularia Oryzae, in Thailand. Plant Pathol.; 2022; 71, pp. 322-333. [DOI: https://dx.doi.org/10.1111/ppa.13460]
21. Perfect, S.E.; Green, J.R. Infection Structures of Biotrophic and Hemibiotrophic Fungal Plant Pathogens. Mol. Plant Pathol.; 2001; 2, pp. 101-108. [DOI: https://dx.doi.org/10.1046/j.1364-3703.2001.00055.x]
22. Zhang, Q. Strategies for Developing Green Super Rice. Proc. Natl. Acad. Sci. USA; 2007; 104, pp. 16402-16409. [DOI: https://dx.doi.org/10.1073/pnas.0708013104]
23. Scheuermann, K.K.; Jia, Y. Identification of a Pi9-Containing Rice Germplasm with a Newly Developed Robust Marker. Phytopathology; 2016; 106, pp. 871-876. [DOI: https://dx.doi.org/10.1094/PHYTO-02-16-0091-R] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/27050577]
24. Berruyer, R.; Adreit, H.; Milazzo, J.; Gaillard, S.; Berger, A.; Dioh, W.; Lebrun, M.-H.; Tharreau, D. Identification and Fine Mapping of Pi33, the Rice Resistance Gene Corresponding to the Magnaporthe Grisea Avirulence Gene ACE1. Theor. Appl. Genet.; 2003; 107, pp. 1139-1147. [DOI: https://dx.doi.org/10.1007/s00122-003-1349-2] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/12838393]
25. Sirisathaworn, T.; Srirat, T.; Longya, A.; Jantasuriyarat, C. Evaluation of Mating Type Distribution and Genetic Diversity of Three Magnaporthe Oryzae Avirulence Genes, PWL-2, AVR-Pii and Avr-Piz-t, in Thailand Rice Blast Isolates. Agric. Nat. Resour.; 2017; 51, pp. 7-14. [DOI: https://dx.doi.org/10.1016/j.anres.2016.08.005]
26. Belyavskaya, L.G.; Belyavsky, Y.V. Interaction of Modern Soybean Varieties with Biological Preparations of Complex Action and Their Impact on the Yield. Mikrobiolohichnyi Zhurnal; 2016; 78, pp. 61-68. [DOI: https://dx.doi.org/10.15407/microbiolj78.03.061]
27. Tai, T.H.; Tanksley, S.D. A Rapid and Inexpensive Method for Isolation of Total DNA from Dehydrated Plant Tissue. Plant Mol. Biol. Rep.; 1990; 8, pp. 297-303. [DOI: https://dx.doi.org/10.1007/BF02668766]
28. Williams, S.T.; Goodfellow, M.; Alderson, G.; Wellington, E.M.H.; Sneath, P.H.A.; Sackin, M.J. Numerical Classification of Streptomyces and Related Genera. Microbiology; 1983; 129, pp. 1743-1813. [DOI: https://dx.doi.org/10.1099/00221287-129-6-1743]
29. Kumar, S.; Stecher, G.; Tamura, K. MEGA7: Molecular Evolutionary Genetics Analysis Version 7.0 for Bigger Datasets. Mol. Biol. Evol.; 2016; 33, pp. 1870-1874. [DOI: https://dx.doi.org/10.1093/molbev/msw054]
30. Roff, D.A.; Bentzen, P. The Statistical Analysis of Mitochondrial DNA Polymorphisms: Chi 2 and the Problem of Small Samples. Mol. Biol. Evol.; 1989; 6, pp. 539-545.
31. Rozas, J.; Ferrer-Matta, A.; Ramos-Onsins, S.E.; Sanchez-Gracia, A. DnaSP 6: DNA sequence polymorphism analysis of large data sets. Mol. Biol. Evol.; 2017; 32, pp. 3299-3302. [DOI: https://dx.doi.org/10.1093/molbev/msx248]
32. Nei, M.; Miller, J.C. A Simple Method for Estimating Average Number of Nucleotide Substitutions within and between Populations from Restriction Data. Genetics; 1990; 125, pp. 873-879. [DOI: https://dx.doi.org/10.1093/genetics/125.4.873]
33. Korneliussen, T.S.; Albrechtsen, A.; Nielsen, R. ANGSD: Analysis of next Generation Sequencing Data. BMC Bioinform.; 2014; 15, 356. [DOI: https://dx.doi.org/10.1186/s12859-014-0356-4] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/25420514]
34. Tajima, F. The Effect of Change in Population Size on DNA Polymorphism. Genetics; 1989; 123, pp. 597-601. [DOI: https://dx.doi.org/10.1093/genetics/123.3.597]
35. Wang, J.C.; Jia, Y.; Wen, J.W.; Liu, W.P.; Liu, X.M.; Li, L.; Jiang, Z.Y.; Zhang, J.H.; Guo, X.L.; Ren, J.P. Identification of Rice Blast Resistance Genes Using International Monogenic Differentials. Crop Prot.; 2013; 45, pp. 109-116. [DOI: https://dx.doi.org/10.1016/j.cropro.2012.11.020]
36. Li, W.; Wang, B.; Wu, J.; Lu, G.; Hu, Y.; Zhang, X.; Zhang, Z.; Zhao, Q.; Feng, Q.; Zhang, H. The Magnaporthe Oryzae Avirulence Gene AvrPiz-t Encodes a Predicted Secreted Protein That Triggers the Immunity in Rice Mediated by the Blast Resistance Gene Piz-T. Mol. Plant-Microbe Interact.; 2009; 22, pp. 411-420. [DOI: https://dx.doi.org/10.1094/MPMI-22-4-0411] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/19271956]
37. Zhu, M.; Wang, L.; Pan, Q. Identification and Characterization of a New Blast Resistance Gene Located on Rice Chromosome 1 through Linkage and Differential Analyses. Phytopathology; 2004; 94, pp. 515-519. [DOI: https://dx.doi.org/10.1094/PHYTO.2004.94.5.515]
38. Jia, Y.; Valent, B.; Lee, F.N. Determination of Host Responses to Magnaporthe Grisea on Detached Rice Leaves Using a Spot Inoculation Method. Plant Dis.; 2003; 87, pp. 129-133. [DOI: https://dx.doi.org/10.1094/PDIS.2003.87.2.129] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/30812916]
39. Telebanco-Yanoria, M.J.; Koide, Y.; Fukuta, Y.; Imbe, T.; Kato, H.; Tsunematsu, H.; Kobayashi, N. Development of Near-Isogenic Lines of Japonica-Type Rice Variety Lijiangxintuanheigu as Differentials for Blast Resistance. Breed. Sci.; 2010; 60, pp. 629-638. [DOI: https://dx.doi.org/10.1270/jsbbs.60.629]
40. Tian, D.; Guo, X.; Zhang, Z.; Wang, M.; Wang, F. Improving Blast Resistance of the Rice Restorer Line, Hui 316, by Introducing Pi9 or Pi2 with Marker-Assisted Selection. Biotechnol. Biotechnol. Equip.; 2019; 33, pp. 1195-1203. [DOI: https://dx.doi.org/10.1080/13102818.2019.1649095]
41. Ling, Z.; Mew, T.; Wang, J.; Lei, C.; Huang, N. Development of Chinese Near-Isogenic Lines of Rice and Their Differentiating Ability to Pathogenic Races of Pyricularia Grisea. Sci. Agric. Sin.; 2000; 33, pp. 1-8.
42. Hu, Z.-J.; Huang, Y.-Y.; Lin, X.-Y.; Feng, H.; Zhou, S.-X.; Xie, Y.; Liu, X.-X.; Liu, C.; Zhao, R.-M.; Zhao, W.-S. Loss and Natural Variations of Blast Fungal Avirulence Genes Breakdown Rice Resistance Genes in the Sichuan Basin of China. Front. Plant Sci.; 2022; 13, 788876. [DOI: https://dx.doi.org/10.3389/fpls.2022.788876]
43. Chuma, I.; Isobe, C.; Hotta, Y.; Ibaragi, K.; Futamata, N.; Kusaba, M.; Yoshida, K.; Terauchi, R.; Fujita, Y.; Nakayashiki, H. Multiple Translocation of the AVR-Pita Effector Gene among Chromosomes of the Rice Blast Fungus Magnaporthe Oryzae and Related Species. PLoS Pathog.; 2011; 7, e1002147. [DOI: https://dx.doi.org/10.1371/journal.ppat.1002147]
44. Yoshida, K.; Saunders, D.G.O.; Mitsuoka, C.; Natsume, S.; Kosugi, S.; Saitoh, H.; Inoue, Y.; Chuma, I.; Tosa, Y.; Cano, L.M. Host Specialization of the Blast Fungus Magnaporthe Oryzae Is Associated with Dynamic Gain and Loss of Genes Linked to Transposable Elements. BMC Genom.; 2016; 17, 370. [DOI: https://dx.doi.org/10.1186/s12864-016-2690-6] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/27194050]
45. Wang, Q.; Li, J.; Lu, L.; He, C.; Li, C. Novel Variation and Evolution of AvrPiz-t of Magnaporthe Oryzae in Field Isolates. Front. Genet.; 2020; 11, 746. [DOI: https://dx.doi.org/10.3389/fgene.2020.00746]
46. Jeuken, M.J.W.; Zhang, N.W.; McHale, L.K.; Pelgrom, K.; Den Boer, E.; Lindhout, P.; Michelmore, R.W.; Visser, R.G.F.; Niks, R.E. Rin4 Causes Hybrid Necrosis and Race-Specific Resistance in an Interspecific Lettuce Hybrid. Plant Cell; 2009; 21, pp. 3368-3378. [DOI: https://dx.doi.org/10.1105/tpc.109.070334]
47. Chen, X.; Shang, J.; Chen, D.; Lei, C.; Zou, Y.; Zhai, W.; Liu, G.; Xu, J.; Ling, Z.; Cao, G. AB-lectin Receptor Kinase Gene Conferring Rice Blast Resistance. Plant J.; 2006; 46, pp. 794-804. [DOI: https://dx.doi.org/10.1111/j.1365-313X.2006.02739.x] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/16709195]
48. Fukuoka, S.; Saka, N.; Koga, H.; Ono, K.; Shimizu, T.; Ebana, K.; Hayashi, N.; Takahashi, A.; Hirochika, H.; Okuno, K. Loss of Function of a Proline-Containing Protein Confers Durable Disease Resistance in Rice. Science; 2009; 325, pp. 998-1001. [DOI: https://dx.doi.org/10.1126/science.1175550] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/19696351]
© 2023 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/). Notwithstanding the ProQuest Terms and Conditions, you may use this content in accordance with the terms of the License.
Details
; Shi, Zhufeng 2 ; Li, Chengyun 3 ; Guo, Zhixiang 2 ; Li, Jinbin 2
1 Flower Research Institute, Yunnan Academy of Agricultural Sciences, Kunming 650205, China
2 Yunnan Key Laboratory of Green Prevention and Control of Agricultural Transboundary Pests, Agricultural Environment and Resource Research Institute, Yunnan Academy of Agricultural Sciences, Kunming 650205, China
3 The Ministry of Education Key Laboratory for Agricultural Biodiversity and Pest Management, Yunnan Agricultural University, Kunming 650200, China