Next Article in Journal
Relationship between Key Environmental Factors and the Architecture of Fruit Shape and Size in Near-Isogenic Lines of Cucumber (Cucumis sativus L.)
Next Article in Special Issue
Identification of Elite R-Gene Combinations against Blast Disease in Geng Rice Varieties
Previous Article in Journal
Characterization of the Secreted Acid Phosphatase SapS Reveals a Novel Virulence Factor of Staphylococcus aureus That Contributes to Survival and Virulence in Mice
Previous Article in Special Issue
Nanosheet-Facilitated Spray Delivery of dsRNAs Represents a Potential Tool to Control Rhizoctonia solani Infection
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
IJMS
Volume 23
Issue 22
10.3390/ijms232214032
ijms-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

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

Genome-Wide Association Study Identifies a Rice Panicle Blast Resistance Gene Pb3 Encoding NLR Protein

by
Lu Ma
1,†,
Yao Yu
1,†,
Changqing Li
1,
Panting Wang
2,
Kunquan Liu
1,
Wenjing Ma
3,
Wei Wang
1,
Yunxin Fan
1,
Ziwei Xiong
1,
Tingting Jiang
1,
Jingran Zhang
1,
Zhixue Wang
1,
Jianfei Wang
1,
Hongsheng Zhang
1 and
Yongmei Bao
1,*
1
State Key Laboratory of Crop Genetics and Germplasm Enhancement, College of Agriculture, Jiangsu Collaborative Innovation Center for Modern Crop Production, Cyrus Tang Innovation Center for Crop Seed Industry, Jiangsu Province Engineering Research Center of Seed Industry Science and Technology, Nanjing Agricultural University, Nanjing 210095, China
2
School of Life Sciences and Biotechnology, Shanghai Jiao Tong University, Shanghai 200240, China
3
School of Agriculture and Biology, Shanghai Jiao Tong University, Shanghai 200240, China
*
Author to whom correspondence should be addressed.
†
These authors contributed equally to this work.
Int. J. Mol. Sci. 2022, 23(22), 14032; https://doi.org/10.3390/ijms232214032
Submission received: 17 October 2022 / Revised: 9 November 2022 / Accepted: 11 November 2022 / Published: 14 November 2022
(This article belongs to the Special Issue Molecular Genetics of Rice Disease Resistance)
Browse Figures
Review Reports Versions Notes

Abstract

:
Rice blast is a worldwide fungal disease that seriously affects the yield and quality of rice. Identification of resistance genes against rice blast disease is one of the effective ways to control this disease. However, panicle blast resistance genes, which are useful in the fields, have rarely been studied due to the difficulty in phenotypic identification and the environmental influences. Here, panicle blast resistance-3 (Pb3) was identified by a genome-wide association study (GWAS) based on the panicle blast resistance phenotypes of 230 Rice Diversity Panel I (RDP-I) accessions with 700,000 single-nucleotide polymorphism (SNP) markers. A total of 16 panicle blast resistance loci (PBRLs) within three years including one repeated locus PBRL3 located in chromosome 11 were identified. In addition, 7 genes in PBRL3 were identified as candidate genes by haplotype analysis, which showed significant differences between resistant and susceptible varieties. Among them, one nucleotide-binding domain and Leucine-rich Repeat (NLR) gene Pb3 was highly conserved in multiple resistant rice cultivars, and its expression was significantly induced after rice blast inoculation. Evolutionary analysis showed that Pb3 was a typical disease resistance gene containing coiled-coil, NB-ARC, and LRR domains. T-DNA insertion mutants and CRISPR lines of Pb3 showed significantly reduced panicle blast resistance. These results indicate that Pb3 is a panicle blast resistance gene and GWAS is a rapid method for identifying panicle blast resistance in rice.

1. Introduction

Rice blast, caused by the fungi Magnaporthe oryzae, is one of the most serious and recurrent diseases in rice-growing regions worldwide [1]. It is estimated that every year, the disease-destroyed rice could feed 60 million people [2,3]. Breeding new resistant varieties with resistance (R) genes are the most effective way to solve the problem [4]. To date, more than 100 rice blast resistance (R) genes have been identified and 35 of them have been cloned [5]. However, panicle blast resistance genes have rarely been studied due to the difficulty in phenotypic identification and the environmental influences [6]. Panicle blast resistance-1 (Pb1) was the first cloned panicle blast resistance gene encoding a coiled-coil-nucleotide binding site-leucine rich repeat (CC-NBS-LRR) protein, which was isolated from the indica cultivar Modan and conferred durable and broad-spectrum resistance to rice blast [7]. Pi25 was derived from the indica rice variety Gumei 2 and provided both panicle blast and leaf blast resistance [8]. The Pigm contains two NBS-LRR proteins, PigmR and PigmS, which endow rice with blast resistance and yield balance [9]. Pb2 is an NLR gene located at the end of chromosome 11 identified by GWAS and complementation transgenic plants phenotype analysis confirmed that it has resistance to panicle blast and leaf blast resistance [10].
Among the 35 cloned resistance genes, 31 genes encode the nucleotide binding-leucine rich repeats (NBS–LRR) protein [11]. The NBS domain contains three highly conserved kinase structures that bind to phosphoric acid to obtain energy for defense against pathogens, and the C-terminal LRR participates in protein–protein interactions of R and Avr-genes [11]. Pib is the first NBS-LRR-type rice blast resistance gene that was cloned in 1999 [12]. The NLR gene is often repeated in tandem to regulate rice blast resistance [13]. For example, Pi2, Pi9, Pi50, Pigm, and Piz-t are all located in the pericentric region of chromosome 6, belonging to the same genetic locus [5]. In addition, Pi-d2, pi21, Ptr, and bsr-k1 encode receptor-like kinase, proline-rich protein, atypical four Armadillo repeats protein, and tetratricopeptide repeats protein, respectively [14,15,16,17]. In the field cultivation, Blast resistance mediated by single or a few major blast resistance genes (R genes) is easily overcome by M. oryzae. However, partial resistance maintains continuous field resistance, which is controlled by minor polygenes or the quantitative trait locus (QTL) [18]. Therefore, identifying new resistant QTLs from rice germplasms is an effective way to control rice blast [19].
Genome-wide association study (GWAS) is a genome-wide genetic marker detection technology, which has been developed to be an important method to mine the inheritance of complex diseases [20]. Compared with traditional bi-parental linkage analysis, GWAS has a high efficiency to identify QTLs in natural varieties [21,22]. In 2005, GWAS was first applied to human disease age-related macular degeneration [23]. With the advent of next-generation sequencing technology, a large amount of genotype data has been obtained in plants, including single-nucleotide polymorphisms (SNPs), insertions and deletions (INDEL), and structure variation (SV), and these played an important role in the development of the GWAS method. In the past 10 years, GWAS has grown into a popular and powerful tool to mine genes/QTL for complex agricultural traits [24]. In wheat, 160 resistance loci associated with stem rust resistance including 21 newly discovered loci were identified by GWAS with 283 durum wheat varieties and 26,439 SNP markers [25]. In maize, 39 SNP loci associated with fumonisin resistance were identified by GWAS by utilizing GBS (Genotyping-By-Sequencing) technology to genotype 270 inbred lines [26]. In soybean, 30 loci associated with flowering and maturation time, and plant height and yield were detected by GWAS based on the phenotype data of 113 germplasms and 4442 high-quality SNP markers [27].
GWAS has also been widely used to study rice complex traits such as yield and disease resistance. With 28,445 SNP markers by resequencing 115 rice japonica cultivars, 16 loci associated with amylose content and resistant starch were identified by GWAS [28]. GWAS analysis of 1,345,417 SNPs on the whole genome of 950 widely collected rice varieties detected 10 QTLs associated with yield traits [29]. With 2,888,332 high-confidence SNPs of 259 rice varieties, 63 BLB (bacterial leaf blight) resistance loci containing 954 genes were identified by GWAS [30]. Ten genomic regions (GRs) with 147 single-nucleotide polymorphisms associated with Rice black-streaked dwarf virus (RBSDV) resistance were identified by GWAS with 1956 accessions from the 3000 Rice Genomes Project (3K RGP) [31]. Rice diversity panel 1 (RDP-1), consisting of 413 O. sativa cultivars collected from 82 countries, is publicly available and widely applied to identify resistance-related loci by GWAS [32,33]. By using 162 rice cultivars of RDP-1, 31 QTLs associated with blast resistance were identified by inoculating eight rice blast isolates from four African counties [34]. Totals of 373, 356, and 336 rice cultivars in RDP-1 were evaluated in the field in three major rice production areas of China, and 74 candidate genes in the region of 16 LAFBRs were identified by GWAS [33]. Combining 362 rice cultivars of RDP-1 and 700K high-density genetic markers, 97 loci were identified associated with blast resistance (LABRs) against five M. oryzae isolates [35]. The Rice Diversity Panel II (RDP-II), which contains 584 rice accessions and is genotyped with 700K SNPs, was inoculated with three M. oryzae isolates. Then, 27 loci associated with rice blast resistance (LABRs) were detected by GWAS, and a partial resistance gene PiPR1 was identified [36]. Pb2 encodes an NBS-LRR protein, identified in the GWAS of 230 cultivars of RDP-1, which enhanced the resistance of panicle blast and leaf blast without affecting agronomic traits and can be used for rice blast resistance breeding [10]. GWAS analysis was carried out by using broad-spectrum resistant rice variety Digu and 66 no broad-spectrum resistant accessions, and found that the natural variation of bsr-d1 promoter altered resistance to rice blast [37]. bsr-d1 is a C2H2-type transcription factor that is regulated by the MYB transcription factor and enhances blast resistance by inhibiting H2O2 degradation, and has no obvious effect on rice yield and quality. Therefore, bsr-d1 has a broad application prospect in rice breeding for blast resistance.
In this study, we evaluated the panicle blast resistance of 230 cultivars derived from RDP-1 by field injection inoculation with the M. oryzae strain 2014-290-B27. By GWAS analysis of population phenotypes and 700K SNPs, 16 panicle blast resistance loci (PBRL) within three years including one repeated locus PBRL3 located in chromosome 11 were identified. In addition, 7 genes in PBRL3 were identified as candidate genes by haplotype analysis, which showed a significant difference between resistant and susceptible varieties. Among them, one Nucleotide-binding domain and Leucine-rich Repeat (NLR) gene Pb3 was highly conserved in multiple resistant rice cultivars, and its expression was significantly induced after rice blast inoculation. Compared with the wild-type, the T-DNA insertion mutant lines and CRISPR knockout transgenic lines showed significantly reduced resistance to leaf blast and panicle blast. In summary, this study identified a Nucleotide-binding domain and Leucine-rich Repeat (NLR) panicle blast resistance gene Pb3 by GWAS analysis, which can significantly improve rice panicle blast and leaf blast resistance.

2. Results

2.1. Phenotype Characteristics of Panicle Blast Resistance in the Natural Population

The natural population used in this study consisted of 230 rice accessions of RDP-1, including 43 indica (IND), 51 tropical japonica (TRJ), 51 temperate japonica (TEJ), 7 aromatic (ARO), 45 aus (AUS), and 34 area mixture (ADM) (Figure 1A). A neighbor-joining analysis was constructed to understand the genetic diversity of 230 accessions used in this study and classified the population into six groups (Figure 1B). To evaluate the panicle blast resistance of the 230 accessions, we inoculated the population with blast strain 2014-290-B27 and evaluated the percentage of diseased grains for each accession in the field in 2017, 2018, and 2019 (Figure 2A–C). Based on the phenotype in three experiments, we analyzed the differences in panicle blast resistance among sub-populations of RDP-1. In 2017, the percentages of diseased grains were 35.57%, 17.10%, 51.08%, 26.00%, 56.15%, and 26.43% in IND, AUS, ADM, ARO, TEJ, and TRJ, respectively (Figure 2D). In 2018, the percentages of diseased grains of IND, AUS, ADM, ARO, TEJ, and TRJ were 11.67%, 6.68%, 35.12%, 4.14%, and 40.87%, respectively (Figure 2E). In 2019, the percentages of diseased grains of IND, AUS, ADM, ARO, TEJ, and TRJ were 28.44%, 16.64%, 44.71%, 14.14%, and 55.59%, respectively (Figure 2F).

2.2. Identification of PBRLs by Using GWAS

We performed GWAS for the three experiments with the Efficient Mixed-Model Association eXpedited (EMMAX), which accounts for both population structure and kinship [38]. The quantile–quantile plots showed that the model of the three environments could be used to identify association signals (Supplementary Figure S1). A total of 16 panicle blast resistance loci (PBRLs) mainly distributed on chromosomes 6 and 10 were identified by GWAS with high-quality genetic markers filtered from 700K SNPs (Figure 3A–C). Among the 16 PBRLs, 9 loci were co-localized with previously mapped or cloned R genes or loci (Table 1). PBRL1 was located at the terminal of chromosome 2 and co-localized with the previously identified Pig(t) [39]. PBRL6 was located in the 20 Mb region of chromosome 4 and contained two reported resistance genes Pi-21 or PBRL-14 locus [10,15]. Among the 16 PBRLs, PBRL3 was the only repeated locus detected both in 2017 and 2019, and it was located in the Pik locus at the end of chromosome 11 [40].

2.3. Analysis of Candidate Genes for Panicle Blast Resistance

To identify candidate resistance genes underlying the stable PBRL3 locus, we investigated the haplotype block structure based on locus within 500 kb intervals; however, it was difficult to find the LD block in the interval (Figure 4A–C). According to genome-wide LD, decay rates of indica and japonica were estimated at ∼200 kb [41], and a 200 kb interval (Chr11:27.1 Mb–27.3 Mb) was selected to predict candidate genes. With reference to the Nipponbare genome sequence and gene annotation, 15 candidate genes with SNPs on the promoter or exons region were obtained. Haplotype analysis results showed that 7 candidate genes with different haplotypes had significant phenotypic differences in the natural population and were selected as potential candidate genes for further analysis (Figure 5A). Among them, LOC_Os11g44990 and LOC_Os11g45090 encode CC-NBS-LRR proteins, LOC_Os11g44910 encodes a DEAD-box RNA helicase, LOC_Os11g44930 encodes pentatricopeptide repeat domain protein, LOC_Os11g44950 encodes a glycosyl hydrolase, and LOC_Os11g44890 and LOC_Os11g45030 encode expressed proteins (Table 2). To investigate expression patterns of candidate genes after M.oryzae infection, we performed quantitative real-time PCR (qRT-PCR) analysis of resistant japonica landrace Bodao after seedling blast inoculation. The results showed that LOC_Os11g44910, LOC_Os11g449300, LOC_Os11g44950, LOC_Os11g44990, and LOC_Os11g45090 were significantly induced after inoculation. Among them, LOC_Os11g45090, which was most significantly induced by rice blast fungi and encodes NBS-LRR domain protein, was reported to be involved in rice blast resistance [42]. We concluded that LOC_Os11g45090 was the candidate gene in PBRL3 and then was named Pb3 (Figure 5B).

2.4. Pb3 Is an R Gene That Regulates Panicle Blast Resistance

Pb3 encodes a protein with a CC-NBS-LRR domain. In the high-quality 700K SNPs dataset, six SNPs causing amino acid changes were found in exons. In addition, SNP11-27282578 was located on the coiled-coil domain, SNP11-27282896 was located on the predicted NB-ARC domain, and SNP11-27284006 and SNP11-27285004 were located on the LRR domains. Moreover, the variation in the SNP11-27285004 in natural populations (CAG to TAG) leads to premature termination of Pb3 protein translation, suggesting that the terminal LRR domain may be involved in resistance to rice blast. Haplotype analysis showed that there were 5 haplotypes of Pb3 in the natural population (Figure 6A,B). The Hap.A haplotype was significantly more resistant to panicle blast than Hap.C and Hap.D (Figure 6C,D). The resistance haplotype A mainly contains TEJ cultivars, while the susceptible haplotypes C and D include 20% indica and 60% japonica (Figure 6E). To study the genetic diversity of Pb3 in natural populations, we calculated the nucleotide diversity (π value) of the chromosome 11 region of 230 varieties by VCF tools, and the results showed that the π value of the Pb3 locus was 0.0002, which was significantly lower than the average nucleotide diversity of the chromosome 11 (π value = 0.0003), illustrating the sequence conservation of Pb3 in natural populations (Figure 6F). We carried out evolutionary analysis on the cloned rice blast R gene, and the results showed that Pb3 was a typical NLR gene, which was highly homologous to Pit and was predicted to play similar functions (Figure 7A,B).
To investigate whether Pb3 was involved in rice blast resistance, we obtained T-DNA insertion lines (pb3) from Salk Institute Genomic Analysis Laboratory (Figure 8A). PCR analysis confirmed T-DNA insertions in the exon region (Supplementary Figure S2). Then, we inoculated mutant and wild-type Dongjin at the booting and seedling stages, respectively. Blast vaccination results showed that the average lesion length and number of Pb3 were significantly higher than those of the wild type (Figure 8B,C,E). Moreover, 5 mutants were more susceptible to the M. oryzae strain than WT by blast inoculation at the booting stage (Figure 8D,F). To further verify the rice blast resistance of Pb3, we knocked out this gene in the Bodao background with the CRISPR/Cas9 gene editing system. A total of 3 homozygous CRISPR lines (type1, type2, and type3) were obtained (Figure 9A,B). The mutation target was located on the LRR domain at the end of the gene. Inoculation of the booting stage results showed that 5 mutants of 2 lines were more susceptible than WT to the M. oryzae strain (Figure 9C,D). It showed that the LRR domain confers resistance to rice blast in Pb3, which was consistent with the results of haplotype analysis. These results further validated that Pb3 was a new R gene in the Pik locus that regulates both leaf blast and panicle blast resistance in rice.

3. Discussion

Rice blast is one of the critical rice diseases, and leaf blast and panicle blast are the most common types of this disease. Compared with leaf blast, panicle blast is considered to be more destructive, which can cause direct yield losses up to 70% or even 100% in fields by affecting grain sterility, rotting the branch and neck, and even losing the entire panicle [43,44]. Isolating the panicle blast resistance gene and introducing it into the susceptible elite cultivars is the most economical and effective way to control rice blast disease. Blast resistance is generally classified into complete and partial resistance [45]. Complete resistance to blast is qualitative and race-specific and is usually controlled by a major gene [4]. However, this type of resistance is easily overcome by the variety of races in the field. Partial resistance controlled by minor polygenes or the quantitative trait locus (QTL) confers broad-spectrum and durable resistance to M. oryzae.
The Rice Diversity Panel (RDP-1) contains substantial genetic and phenotypic diversity and has been used for GWAS analysis of various traits in rice [46]. We selected 230 cultivars in RDP-1 for panicle blast resistance identification in the field. Three experimental field panicle blast inoculations were performed on these cultivars. The average diseased grains percentages were 23.33%, 12.53%, 50.37%, and 38.70% in IND, AUS, TEJ, and TRJ, respectively. These results are consistent with previous studies where the TRJ and IND subgroups of rice were more resistant than the TEJ subgroups [33]. We chose the EMMAX model for association analysis and identified a repetitively mapped locus PBRL3 within the Pik locus at the end of chromosome 11. Haplotype analysis and expression pattern analysis identified 5 candidate genes that were significantly induced by M. oryzae. Among them, LOC_Os11g44910 encodes a DEAD-box RNA helicase, which contains a conserved amino acid sequence Asp-Glu-Ala-Asp (DEAD). DEAD box RNA helicase family proteins can hydrolyze ATP to obtain double-stranded energy unwinding RNA, and are reported to be involved in plant leaf growth and the cold stress response [47,48]. LOC_Os11g44930 encodes a pentatricopeptide-repeat (PPR)-domain-containing protein. Plant PPR proteins are usually located in mitochondria or chloroplasts, which edit and modify the RNA transcripts to regulate plant growth and development [49,50]. LOC_Os11g44950 is a member of the glycosyl hydrolase family 3. Glycosyl hydrolase can effectively hydrolyze sugar-based substances and are involved in cell wall remodeling [51]. LOC_Os11g44990 and LOC_Os11g45090 are NBS-LRR-type genes that are widely reported to be involved in rice blast resistance [52]. Furthermore, after inoculation with M. oryzae, the expression of LOC_Os11g45090 significantly increased at 24 h by 4.07-fold. These results indicate that LOC_Os11g45090 (Pb3) is the most likely candidate gene in PBRL3.
In the 230 natural populations, five SNPs causing amino acid changes were detected in the coding region of Pb3. Among them, SNP11-27282578 fell on the CC domain, SNP11-27282896 was located in the NB-ARC domain, and SNP27284006 and SNP11-27285004 fell on the LRR domain. Combining the haplotype analysis results, we found that the difference between resistant haplotype Hap. A and susceptible haplotypes Hap. C and Hap. D contained these four SNPs that were located in the functional domain of Pb3. In NBS-LRR proteins, the CC domain mainly mediates the interaction with transcription factors and structural proteins [53]. NBS has the function of ATP or GTP binding and hydrolysis, which can resist the invasion of pathogens by gaining energy [11]. The easily altered structure of LRR may be related to its ability to specifically recognize pathogen-effector proteins [54]. To verify the function of these domains, we obtained the transgenic plants of Pb3. We first obtained a T-DNA insertion mutant of Pb3, and the insertion position was located in the double coiled-coil region. The susceptible phenotype of the mutant line verified that Pb3 was indeed involved in rice blast resistance. The genotype data revealed an SNP11-27285004 on the LRR domain in the end of Pb3 in natural populations, which can prematurely terminate the translation of the protein, resulting in the deletion of the LRR domain. Therefore, we specifically constructed a CRISPR-Cas9 mutation target located on the terminal LRR domain, and CRISPR showed that it was more susceptible than wild-type thin rice, indicating the importance of the LRR domain for Pb3 to play the role of panicle blast resistance.
Most of the cloned rice blast resistance genes encode NBS-LRR domain proteins. Therefore, we performed phylogenetic tree analysis on Pb3 and cloned R genes and found that Pit has the highest homology with PB3. Pit interacts with OsRac1 through its NB-ARC domain and activates OsRac1, which plays an important role in Pit-mediated immune responses [55]. Therefore, Pb3 may also mediate resistance by binding to a key immune regulatory factor. In the following studies, we can verify the interaction between Pb3 and OsRac1, and can also search for other interacting proteins in combination with yeast two-hybrid assays to research the specific resistance mechanism of Pb3. On the other hand, we verified the function of the LRR domain at the end of Pb3 through the CRISPR plants. The deletion of the LRR domain may affect the recognition of effector proteins of rice blast fungi, further inhibit the activation of the rice defense response, and reduce rice resistance to rice blast. In summary, we identified a repeat site PBRL3 on chromosome 11 by GWAS. Through functional verification of candidate genes, we found that Pb3, which encodes NBS-LRR protein, is a causal gene in PBRL3, which can regulate the resistance of rice to leaf blast and panicle blast. In this work, we carried out GWAS analysis to quickly and accurately obtain the rice blast resistance R gene, which provides a theoretical basis for mining a large number of panicle blast resistance genes. With the continuous development of sequencing technology, it will be more accurate to locate panicle blast resistance genes through GWAS in the future. However, due to the influence of environmental factors, identification of resistance to panicle blast in the field is still a difficult problem, which is crucial to the results of GWAS analysis.

4. Materials and Methods

4.1. Plant Materials and Fungal Growth

RDP-1, which contains 413 O. sativa accessions from 82 countries, was provided by Cornell University [46]. According to the growth period of these accessions, we selected 230 accessions from the panel to ensure the growth period consistency of rice. The whole panel was divided into six subpopulations, including Indica, Aus, Tropical japonica, Temperate japonica, Aromatic, and Admixtures [46]. These germplasms were planted in an experimental field in Nanjing for panicle blast inoculation. The seedlings of these germplasms were cultivated in 96-hole plugs in the greenhouse for seedling leaf blast inoculation. The rice blast isolate 2014-290B27 was provided by the Jiangsu Academy of Agricultural Sciences. The blast fungi were grown on corn rice straw agar plates at 28 °C for seven days and then transferred in black light (20 W) in an incubator at 28 °C for sporulation and culture for 7–10 days to promote spore production. They were then washed with sterilized distilled water to a spore suspension with a concentration of 1 × 105 /mL [16,56].

4.2. Inoculation and Phenotypic Survey

In 2017 and 2018, natural population materials were planted in Nanjing, and inoculation was carried out at the booting stage. Each rice accession was inoculated with 9 rice panicles. The ear was sampled 15 days after inoculation, and the resistance survey was measured by calculating the percentage of the number of diseased grains in the whole panicle of rice (0–100%) [57]. For the leaf blast resistance identification of T-DNA insertion mutants and CRISPR plants of Pb3, the seedlings were planted in the greenhouse of Nanjing Agricultural University to the two-leaf and one-heart stage for spray inoculation [58]. The length and number of lesions were recorded after six days. Six plants for each line and three replicates were measured, and the average was calculated for the resistance phenotype. Suyunuo, a highly susceptible japonica rice variety, was planted around the tray as a control.

4.3. Population Structure Analysis and Genome-Wide Association Study

The phylogenetic tree of the 230 accessions was constructed using the MEGA7 program [59]. Multiple sequence alignment of 230 accessions was performed with the ClustalW program with the standard setting. The neighbor-joining (NJ) with maximum composite likelihood method was used to construct phylogenetic trees with a bootstrap value of 1000 replicates in MEGA 7. The beautification of the phylogenetic tree was conducted with iTOL online tools [60]. Based on the mixed linear model, GWAS was used to detect resistance-related loci [61]. The data of 400K SNPs (MAF ≥ 5%, missing rate < 25%) corresponding to the study accessions were filtered from the 700K high-density genetic markers of RDP-1 through the TASSEL 5.0 software. The MLM model was combined with the EMMAX online script (http://genetics.cs.ucla.edu/emmax/install.html, accessed on 20 December 2019) and R packages for analysis. Manhattan plots were drawn by using the R package qqman. p-value < 10−4 was used as the significance threshold to obtain the loci associated with leaf blast.

4.4. Candidate Genes and Haplotype Analysis

The LD heatmap of the associated loci in the GWAS was constructed using the Haploview4.2 [62]. Haplotype blocks were constructed using the confidence interval (r2 > 0.6) [21]. Candidate genes were predicted according to the MSU Rice Genome Annotation Project (http://rice.plantbiology.msu.edu, accessed on 10 June 2020). We analyzed the polymorphisms in the candidate region that were located on the promoter and caused amino acid changes. Haplotype analysis of these candidate genes was performed to screen out significant difference genes by using the Haploview4.2.

4.5. Genomic Nucleotide Diversity and Phylogenetic Analysis

To clarify the Pb3 genetic differentiation of the natural population, the PLINK online script was used to extract SNPs on chromosome 11 in natural populations [63]. Then, nucleotide diversity (π) was calculated using VCFtools software [64]. Sliding 100 kb windows were used during the calculation with a 10 kb sliding step. The multiple sequences of R genes were aligned with MEGA7, then the phylogenetic trees were calculated and constructed using the neighbor-joining method within 1000 replicates. The conserved domain of the R gene was predicted by the Pfam online website (http://pfam.xfam.org/, accessed on 25 November 2021), and the gene structure was drawn by the IBS1.0.3 software [65].

4.6. RNA Extraction and qRT-PCR Analysis

Variety leaves infected with M. oryzae strains were collected at different time points after inoculation. The samples were immediately frozen with liquid nitrogen and stored in the refrigerator at −80 °C. The total RNA was extracted from leaves using TRIpure Reagent (Wuxi, China, Bio Teke Corporation, http://www.bioteke.com, accessed on 5 June 2021). First-strand cDNA synthesis was performed using HiScript II Q-RT Super Mix (Nanjing, China, Vazyme, http://www.vazyme.com, accessed on 6 June 2021). qRT-PCR was performed using the AceQ® qPCR SYBR Green Master Mix. The primers in Supplemental Table S1 were amplified for the genes’ expression analysis. The rice actin gene (LOC_Os03g50885) was used as the internal control. Each experiment was performed with three biological samples and each sample was assayed with three technical replications.

4.7. Generation of Transgenic Plants

The CRISPR-Cas9 gene mutation system was obtained from the Chen Qijun Laboratory, College of Biology, China Agricultural University, Beijing China [66]. The target primers were designed using the online website (http://skl.scau.edu.cn/primerdesign/, accessed on 12 September 2021). The CRISPR vectors pBUE411-2gR containing the Pb3-specific silencing target were transformed into japonica rice Bodao by an agrobacterium-mediated transformation method [67]. We obtained six CRISPR mutants of Pb3 containing three editing types: type1:KO-1/ KO-2; type2: KO-3/KO-4/KO-5; type3: KO-6. The T-DNA insertion mutant of the Dongjin background of Pb3 was purchased from Salk Institute Genomic Analysis Laboratory (http://signal.salk.edu/index.html, accessed on 11 July 2021). Sequencing results showed that an 8275 bp T-DNA was inserted into the exon of Pb3 in the mutant. We obtained 5 homozygous mutant plants (pb3-1 to pb3-5) from 10 T0 generation seeds of the mutant. The primers in Supplemental Table S1 were amplified for positive verification.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/ijms232214032/s1.

Author Contributions

Experimental design, L.M., Y.Y., C.L. and Y.B.; field phenotyping, L.M., Y.Y., C.L., P.W., K.L. and Y.B.; data analysis, L.M. and Y.Y.; gene cloning and experimental validation, L.M., Y.Y., C.L., P.W., K.L., W.M., W.W. and Y.F.; percentage diseased grains investigation, L.M., Y.Y., C.L., Z.X., T.J. and J.Z.; writing—original draft preparation, L.M. and Y.B.; writing—review and editing, L.M., Z.W., J.W., H.Z. and Y.B. All authors have read and agreed to the published version of the manuscript.

Funding

The project was supported by the Guidance Foundation, the Sanya Institute of Nanjing Agricultural University (NAUSY-MS06), the National Key Project for Jiangsu Agriculture Science and Technology Innovation Fund (CX(2019)1008), and grants from the National Key Research and Development Program of China (2020YFE0202900), the Natural Science Foundation of China (32072033, 31871602, 31171516), the Natural Science Foundation of Jiangsu Province (BK20211212), and the Fundamental Research Funds for the Central Universities (KYZZ2022003).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

We thank Yongfeng Liu for providing the rice blast isolates.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Dean, R.; Van Kan, J.A.; Pretorius, Z.A.; Hammond-Kosack, K.E.; Di Pietro, A.; Spanu, P.D.; Rudd, J.J.; Dickman, M.; Kahmann, R.; Ellis, J.; et al. The Top 10 fungal pathogens in molecular plant pathology. Mol. Plant Pathol. 2012, 13, 414–430. [Google Scholar] [CrossRef] [Green Version]
  2. Skamnioti, P.; Gurr, S.J. Against the grain: Safeguarding rice from rice blast disease. Trends Biotechnol. 2009, 27, 141–150. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Pennisi, E. Armed and dangerous. Science 2010, 327, 804–805. [Google Scholar] [PubMed] [Green Version]
  4. Wang, C.; Yang, Y.; Yuan, X.; Xu, Q.; Feng, Y.; Yu, H.; Wang, Y.; Wei, X. Genome-wide association study of blast resistance in indica rice. BMC Plant Biol. 2014, 14, 311. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Wang, B.H.; Ebbole, D.J.; Wang, Z.H. The arms race between Magnaporthe oryzae and rice: Diversity and interaction of Avr and R genes. J. Integr. Agr. 2017, 16, 2746–2760. [Google Scholar] [CrossRef]
  6. Liu, Q.; Yang, J.; Zhang, S.; Zhao, J.; Feng, A.; Yang, T.; Wang, X.; Mao, X.; Dong, J.; Zhu, X.; et al. OsGF14b Positively Regulates Panicle Blast Resistance but Negatively Regulates Leaf Blast Resistance in Rice. Mol. Plant Microbe Interact 2016, 29, 46–56. [Google Scholar] [CrossRef] [Green Version]
  7. Hayashi, N.; Inoue, H.; Kato, T.; Funao, T.; Shirota, M.; Shimizu, T.; Kanamori, H.; Yamane, H.; Hayano-Saito, Y.; Matsumoto, T.; et al. Durable panicle blast-resistance gene Pb1 encodes an atypical CC-NBS-LRR protein and was generated by acquiring a promoter through local genome duplication. Plant J. 2010, 64, 498–510. [Google Scholar] [CrossRef]
  8. Chen, J.; Shi, Y.; Liu, W.; Chai, R.; Fu, Y.; Zhuang, J.; Wu, J. A Pid3 allele from rice cultivar Gumei2 confers resistance to Magnaporthe oryzae. J. Genet. Genom. 2011, 38, 209–216. [Google Scholar] [CrossRef]
  9. Deng, Y.; Zhai, K.; Xie, Z.; Yang, D.; Zhu, X.; Liu, J.; Wang, X.; Qin, P.; Yang, Y.; Zhang, G.; et al. Epigenetic regulation of antagonistic receptors confers rice blast resistance with yield balance. Science 2017, 355, 962–965. [Google Scholar] [CrossRef]
  10. Yu, Y.; Ma, L.; Wang, X.; Zhao, Z.; Wang, W.; Fan, Y.; Liu, K.; Jiang, T.; Xiong, Z.; Song, Q.; et al. Genome-Wide Association Study Identifies a Rice Panicle Blast Resistance Gene, Pb2, Encoding NLR Protein. Int. J. Mol. Sci. 2022, 23, 5668. [Google Scholar] [CrossRef]
  11. McHale, L.; Tan, X.; Koehl, P.; Michelmore, R.W. Plant NBS-LRR proteins: Adaptable guards. Genome Biol. 2006, 7, 212. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  12. Wang, Z.X.; Yano, M.; Yamanouchi, U.; Iwamoto, M.; Monna, L.; Hayasaka, H.; Katayose, Y.; Sasaki, T. The Pib gene for rice blast resistance belongs to the nucleotide binding and leucine-rich repeat class of plant disease resistance genes. Plant J. 1999, 19, 55–64. [Google Scholar] [CrossRef] [PubMed]
  13. Tanweer, F.A.; Rafii, M.Y.; Sijam, K.; Rahim, H.A.; Ahmed, F.; Latif, M.A. Current advance methods for the identification of blast resistance genes in rice. Comptes Rendus Biol. 2015, 338, 321–334. [Google Scholar] [CrossRef]
  14. Chen, X.; Shang, J.; Chen, D.; Lei, C.; Zou, Y.; Zhai, W.; Liu, G.; Xu, J.; Ling, Z.; Cao, G.; et al. A B-lectin receptor kinase gene conferring rice blast resistance. Plant J. 2006, 46, 794–804. [Google Scholar] [CrossRef] [PubMed]
  15. Fukuoka, S.; Saka, N.; Koga, H.; Ono, K.; Shimizu, T.; Ebana, K.; Hayashi, N.; Takahashi, A.; Hirochika, H.; Okuno, K.; et al. Loss of Function of a Proline-Containing Protein Confers Durable Disease Resistance in Rice. Science 2009, 325, 998–1001. [Google Scholar] [CrossRef]
  16. Zhao, H.; Wang, X.; Jia, Y.; Minkenberg, B.; Wheatley, M.; Fan, J.; Jia, M.H.; Famoso, A.; Edwards, J.D.; Wamishe, Y.; et al. The rice blast resistance gene Ptr encodes an atypical protein required for broad-spectrum disease resistance. Nat. Commun. 2018, 9, 2039. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  17. Zhou, X.; Liao, H.; Chern, M.; Yin, J.; Chen, Y.; Wang, J.; Zhu, X.; Chen, Z.; Yuan, C.; Zhao, W.; et al. Loss of function of a rice TPR-domain RNA-binding protein confers broad-spectrum disease resistance. Proc. Natl. Acad. Sci. USA 2018, 115, 3174–3179. [Google Scholar] [CrossRef] [Green Version]
  18. Nelson, R.; Wiesner-Hanks, T.; Wisser, R.; Balint-Kurti, P. Navigating complexity to breed disease-resistant crops. Nat. Rev. Genet. 2018, 19, 21–33. [Google Scholar] [CrossRef]
  19. Xiao, N.; Wu, Y.; Pan, C.; Yu, L.; Chen, Y.; Liu, G.; Li, Y.; Zhang, X.; Wang, Z.; Dai, Z.; et al. Improving of Rice Blast Resistances in Japonica by Pyramiding Major R Genes. Front. Plant Sci. 2016, 7, 1918. [Google Scholar] [CrossRef] [Green Version]
  20. Burghardt, L.T.; Young, N.D.; Tiffin, P. A Guide to Genome-Wide Association Mapping in Plants. Curr. Protoc. Plant Biol. 2017, 2, 22–38. [Google Scholar] [CrossRef]
  21. Yano, K.; Yamamoto, E.; Aya, K.; Takeuchi, H.; Lo, P.C.; Hu, L.; Yamasaki, M.; Yoshida, S.; Kitano, H.; Hirano, K.; et al. Genome-wide association study using whole-genome sequencing rapidly identifies new genes influencing agronomic traits in rice. Nat. Genet. 2016, 48, 927–934. [Google Scholar] [CrossRef] [PubMed]
  22. Tibbs Cortes, L.; Zhang, Z.; Yu, J. Status and prospects of genome-wide association studies in plants. Plant Genome 2021, 14, e20077. [Google Scholar] [CrossRef] [PubMed]
  23. Klein, R.J.; Zeiss, C.; Chew, E.Y.; Tsai, J.Y.; Sackler, R.S.; Haynes, C.; Henning, A.K.; SanGiovanni, J.P.; Mane, S.M.; Mayne, S.T.; et al. Complement factor H polymorphism in age-related macular degeneration. Science 2005, 308, 385–389. [Google Scholar] [CrossRef] [PubMed]
  24. Lu, J.; Wang, C.; Zeng, D.; Li, J.; Shi, X.; Shi, Y.; Zhou, Y. Genome-Wide Association Study Dissects Resistance Loci against Bacterial Blight in a Diverse Rice Panel from the 3000 Rice Genomes Project. Rice 2021, 14, 22. [Google Scholar] [CrossRef]
  25. Megerssa, S.H.; Ammar, K.; Acevedo, M.; Brown-Guedira, G.; Ward, B.; Degete, A.G.; Randhawa, M.S.; Sorrells, M.E. Multiple-Race Stem Rust Resistance Loci Identified in Durum Wheat Using Genome-Wide Association Mapping. Front. Plant Sci. 2020, 11, 598509. [Google Scholar] [CrossRef]
  26. Samayoa, L.F.; Cao, A.; Santiago, R.; Malvar, R.A.; Butron, A. Genome-wide association analysis for fumonisin content in maize kernels. BMC Plant Biol. 2019, 19, 166. [Google Scholar] [CrossRef] [Green Version]
  27. Zatybekov, A.; Abugalieva, S.; Didorenko, S.; Gerasimova, Y.; Sidorik, I.; Anuarbek, S.; Turuspekov, Y. GWAS of agronomic traits in soybean collection included in breeding pool in Kazakhstan. BMC Plant Biol. 2017, 17 (Suppl. 1), 179. [Google Scholar] [CrossRef] [Green Version]
  28. Biselli, C.; Volante, A.; Desiderio, F.; Tondelli, A.; Gianinetti, A.; Finocchiaro, F.; Taddei, F.; Gazza, L.; Sgrulletta, D.; Cattivelli, L.; et al. GWAS for Starch-Related Parameters in Japonica Rice (Oryza sativa L.). Plants 2019, 8, 292. [Google Scholar] [CrossRef] [Green Version]
  29. Huang, X.; Zhao, Y.; Wei, X.; Li, C.; Wang, A.; Zhao, Q.; Li, W.; Guo, Y.; Deng, L.; Zhu, C.; et al. Genome-wide association study of flowering time and grain yield traits in a worldwide collection of rice germplasm. Nat. Genet. 2011, 44, 32–39. [Google Scholar] [CrossRef]
  30. Shu, X.; Wang, A.; Jiang, B.; Jiang, Y.; Xiang, X.; Yi, X.; Li, S.; Deng, Q.; Wang, S.; Zhu, J.; et al. Genome-wide association study and transcriptome analysis discover new genes for bacterial leaf blight resistance in rice (Oryza sativa L.). BMC Plant Biol. 2021, 21, 255. [Google Scholar] [CrossRef]
  31. Liu, Q.; Lan, G.; Zhu, Y.; Chen, K.; Shen, C.; Zhao, X.; Zhang, F.; Xu, J.; Li, Z. Genome-Wide Association Study on Resistance to Rice Black-Streaked Dwarf Disease Caused by Rice black-streaked dwarf virus. Plant Dis. 2021, 105, 607–615. [Google Scholar] [CrossRef] [PubMed]
  32. McCouch, S.R.; Wright, M.H.; Tung, C.W.; Maron, L.G.; McNally, K.L.; Fitzgerald, M.; Singh, N.; De Clerck, G.; Agosto-Perez, F.; Korniliev, P.; et al. Open access resources for genome-wide association mapping in rice. Nat. Commun. 2016, 7, 10532. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. Zhu, D.; Kang, H.; Li, Z.; Liu, M.; Zhu, X.; Wang, Y.; Wang, D.; Wang, Z.; Liu, W.; Wang, G.L. A Genome-Wide Association Study of Field Resistance to Magnaporthe Oryzae in Rice. Rice 2016, 9, 44. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. Mgonja, E.M.; Balimponya, E.G.; Kang, H.; Bellizzi, M.; Park, C.H.; Li, Y.; Mabagala, R.; Sneller, C.; Correll, J.; Opiyo, S.; et al. Genome-Wide Association Mapping of Rice Resistance Genes Against Magnaporthe oryzae Isolates from Four African Countries. Phytopathology 2016, 106, 1359–1365. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  35. Kang, H.; Wang, Y.; Peng, S.; Zhang, Y.; Xiao, Y.; Wang, D.; Qu, S.; Li, Z.; Yan, S.; Wang, Z.; et al. Dissection of the genetic architecture of rice resistance to the blast fungus Magnaporthe oryzae. Mol. Plant Pathol. 2016, 17, 959–972. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  36. Liu, M.H.; Kang, H.; Xu, Y.; Peng, Y.; Wang, D.; Gao, L.; Wang, X.; Ning, Y.; Wu, J.; Liu, W.; et al. Genome-wide association study identifies an NLR gene that confers partial resistance to Magnaporthe oryzae in rice. Plant Biotechnol. J. 2019, 18, 1376–1383. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  37. Li, W.; Zhu, Z.; Chern, M.; Yin, J.; Yang, C.; Ran, L.; Cheng, M.; He, M.; Wang, K.; Wang, J.; et al. A Natural Allele of a Transcription Factor in Rice Confers Broad-Spectrum Blast Resistance. Cell 2017, 170, 114–126.e15. [Google Scholar] [CrossRef] [Green Version]
  38. Kang, H.M.; Sul, J.H.; Service, S.K.; Zaitlen, N.A.; Kong, S.Y.; Freimer, N.B.; Sabatti, C.; Eskin, E. Variance component model to account for sample structure in genome-wide association studies. Nat. Genet. 2010, 42, 348–354. [Google Scholar] [CrossRef] [Green Version]
  39. Zhou, J.H.; Wang, J.L.; Xu, J.C.; Lei, C.L.; Ling, Z.Z. Identification and mapping of a rice blast resistance gene Pi-g(t) in the cultivar Guangchangzhan. Plant Pathol. 2004, 53, 191–196. [Google Scholar] [CrossRef]
  40. Ashikawa, I.; Hayashi, N.; Abe, F.; Wu, J.Z.; Matsumoto, T. Characterization of the rice blast resistance gene Pik cloned from Kanto51. Mol. Breed. 2012, 30, 485–494. [Google Scholar] [CrossRef]
  41. Huang, X.; Wei, X.; Sang, T.; Zhao, Q.; Feng, Q.; Zhao, Y.; Li, C.; Zhu, C.; Lu, T.; Zhang, Z.; et al. Genome-wide association studies of 14 agronomic traits in rice landraces. Nat. Genet. 2010, 42, 961–967. [Google Scholar] [CrossRef] [PubMed]
  42. Tabien, R.E.; Li, Z.; Paterson, A.H.; Marchetti, M.A.; Stansel, J.W.; Pinson, S.R.M. Mapping of four major rice blast resistance genes from ‘Lemont’ and ‘Teqing’ and evaluation of their combinatorial effect for field resistance. Theor. Appl. Genet. 2000, 101, 1215–1225. [Google Scholar] [CrossRef]
  43. Khan, M.A.; Bhuiyan, M.R.; Hossain, M.S.; Sen, P.P.; Ara, A.; Siddique, M.A.; Ali, M.A. Neck blast disease influences grain yield and quality traits of aromatic rice. Comptes Rendus Biol. 2014, 337, 635–641. [Google Scholar] [CrossRef]
  44. Fang, N.; Wei, X.; Shen, L.; Yu, Y.; Li, M.; Yin, C.; He, W.; Guan, C.; Chen, H.; Zhang, H.; et al. Fine mapping of a panicle blast resistance gene Pb-bd1 in Japonica landrace Bodao and its application in rice breeding. Rice 2019, 12, 18. [Google Scholar] [CrossRef] [PubMed]
  45. Wang, G.L.; Mackill, D.J.; Bonman, J.M.; McCouch, S.R.; Champoux, M.C.; Nelson, R.J. RFLP mapping of genes conferring complete and partial resistance to blast in a durably resistant rice cultivar. Genetics 1994, 136, 1421–1434. [Google Scholar] [CrossRef]
  46. Zhao, K.; Tung, C.W.; Eizenga, G.C.; Wright, M.H.; Ali, M.L.; Price, A.H.; Norton, G.J.; Islam, M.R.; Reynolds, A.; Mezey, J.; et al. Genome-wide association mapping reveals a rich genetic architecture of complex traits in Oryza sativa. Nat. Commun. 2011, 2, 467. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  47. Nishimura, K.; Ashida, H.; Ogawa, T.; Yokota, A. A DEAD box protein is required for formation of a hidden break in Arabidopsis chloroplast 23S rRNA. Plant J. 2010, 63, 766–777. [Google Scholar] [CrossRef] [PubMed]
  48. Braud, C.; Zheng, W.; Xiao, W. LONO1 encoding a nucleoporin is required for embryogenesis and seed viability in Arabidopsis. Plant Physiol. 2012, 160, 823–836. [Google Scholar] [CrossRef] [Green Version]
  49. Cao, Z.L.; Yu, Q.B.; Sun, Y.; Lu, Y.; Cui, Y.L.; Yang, Z.N. A point mutation in the pentatricopeptide repeat motif of the AtECB2 protein causes delayed chloroplast development. J. Integr. Plant Biol. 2011, 53, 258–269. [Google Scholar] [CrossRef]
  50. Haili, N.; Planchard, N.; Arnal, N.; Quadrado, M.; Vrielynck, N.; Dahan, J.; des Francs-Small, C.C.; Mireau, H. The MTL1 Pentatricopeptide Repeat Protein Is Required for Both Translation and Splicing of the Mitochondrial NADH DEHYDROGENASE SUBUNIT7 mRNA in Arabidopsis. Plant Physiol. 2016, 170, 354–366. [Google Scholar] [CrossRef]
  51. Urbanowicz, B.R.; Bennett, A.B.; Del Campillo, E.; Catala, C.; Hayashi, T.; Henrissat, B.; Hofte, H.; McQueen-Mason, S.J.; Patterson, S.E.; Shoseyov, O.; et al. Structural organization and a standardized nomenclature for plant endo-1,4-beta-glucanases (cellulases) of glycosyl hydrolase family 9. Plant Physiol. 2007, 144, 1693–1696. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  52. Sharma, T.R.; Rai, A.K.; Gupta, S.K.; Vijayan, J.; Devanna, B.N.; Ray, S. Rice Blast Management Through Host-Plant Resistance: Retrospect and Prospects. Agric. Res. 2012, 1, 37–52. [Google Scholar] [CrossRef] [Green Version]
  53. Stirnweis, D.; Milani, S.D.; Jordan, T.; Keller, B.; Brunner, S. Substitutions of two amino acids in the nucleotide-binding site domain of a resistance protein enhance the hypersensitive response and enlarge the PM3F resistance spectrum in wheat. Mol. Plant Microbe Interact 2014, 27, 265–276. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  54. 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, 4004–4014. [Google Scholar] [CrossRef]
  55. Kawano, Y.; Akamatsu, A.; Hayashi, K.; Housen, Y.; Okuda, J.; Yao, A.; Nakashima, A.; Takahashi, H.; Yoshida, H.; Wong, H.L.; et al. Activation of a Rac GTPase by the NLR family disease resistance protein Pit plays a critical role in rice innate immunity. Cell Host Microbe 2010, 7, 362–375. [Google Scholar] [CrossRef] [Green Version]
  56. Guo, M.; Chen, Y.; Du, Y.; Dong, Y.; Guo, W.; Zhai, S.; Zhang, H.; Dong, S.; Zhang, Z.; Wang, Y.; et al. The bZIP transcription factor MoAP1 mediates the oxidative stress response and is critical for pathogenicity of the rice blast fungus Magnaporthe oryzae. PLoS Pathog. 2011, 7, e1001302. [Google Scholar] [CrossRef] [Green Version]
  57. Wang, R.; Fang, N.; Guan, C.; He, W.; Bao, Y.; Zhang, H. Characterization and Fine Mapping of a Blast Resistant Gene Pi-jnw1 from the japonica Rice Landrace Jiangnanwan. PLoS ONE 2016, 11, e0169417. [Google Scholar] [CrossRef] [Green Version]
  58. Park, C.H.; Chen, S.; Shirsekar, G.; Zhou, B.; Khang, C.H.; Songkumarn, P.; Afzal, A.J.; Ning, Y.; Wang, R.; Bellizzi, M.; et al. The Magnaporthe oryzae effector AvrPiz-t targets the RING E3 ubiquitin ligase APIP6 to suppress pathogen-associated molecular pattern-triggered immunity in rice. Plant Cell 2012, 24, 4748–4762. [Google Scholar] [CrossRef] [Green Version]
  59. Tamura, K.; Peterson, D.; Peterson, N.; Stecher, G.; Nei, M.; Kumar, S. MEGA5: Molecular evolutionary genetics analysis using maximum likelihood, evolutionary distance, and maximum parsimony methods. Mol. Biol. Evol. 2011, 28, 2731–2739. [Google Scholar] [CrossRef] [Green Version]
  60. Letunic, I.; Bork, P. Interactive Tree Of Life (iTOL) v5: An online tool for phylogenetic tree display and annotation. Nucleic Acids Res. 2021, 49, W293–W296. [Google Scholar] [CrossRef]
  61. Zhang, Z.; Ersoz, E.; Lai, C.Q.; Todhunter, R.J.; Tiwari, H.K.; Gore, M.A.; Bradbury, P.J.; Yu, J.; Arnett, D.K.; Ordovas, J.M.; et al. Mixed linear model approach adapted for genome-wide association studies. Nat. Genet. 2010, 42, 355–360. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  62. Barrett, J.C.; Fry, B.; Maller, J.; Daly, M.J. Haploview: Analysis and visualization of LD and haplotype maps. Bioinformatics 2005, 21, 263–265. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  63. Purcell, S.; Neale, B.; Todd-Brown, K.; Thomas, L.; Ferreira, M.A.; Bender, D.; Maller, J.; Sklar, P.; de Bakker, P.I.; Daly, M.J.; et al. PLINK: A tool set for whole-genome association and population-based linkage analyses. Am. J. Hum. Genet. 2007, 81, 559–575. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  64. Danecek, P.; Auton, A.; Abecasis, G.; Albers, C.A.; Banks, E.; De Pristo, M.A.; Handsaker, R.E.; Lunter, G.; Marth, G.T.; Sherry, S.T. The variant call format and VCFtools. Bioinformatics 2011, 27, 2156–2158. [Google Scholar] [CrossRef] [PubMed]
  65. Liu, W.; Xie, Y.; Ma, J.; Luo, X.; Nie, P.; Zuo, Z.; Lahrmann, U.; Zhao, Q.; Zheng, Y.; Zhao, Y.; et al. IBS: An illustrator for the presentation and visualization of biological sequences. Bioinformatics 2015, 31, 3359–3361. [Google Scholar] [CrossRef] [Green Version]
  66. Zhou, H.; Liu, B.; Weeks, D.P.; Spalding, M.H.; Yang, B. Large chromosomal deletions and heritable small genetic changes induced by CRISPR/Cas9 in rice. Nucleic Acids Res. 2014, 42, 10903–10914. [Google Scholar] [CrossRef] [PubMed]
  67. Wang, M.; Chen, C.; Xu, Y.Y.; Jiang, R.X.; Han, Y.; Xu, Z.H.; Chong, K. A practical vector for efficient knockdown of gene expression in rice (Oryza sativa L.). Plant Mol. Biol. Rep. 2004, 22, 409–417. [Google Scholar] [CrossRef]
Ijms 23 14032 g001 550
Figure 1. Population structure in O. sativa. (A) The large pie chart summarizes the distribution of subpopulations in the 230 O. sativa samples in our diversity panel, and the smaller pie charts on the world map correspond to the country-specific distribution of subpopulations sampled. (B) Phylogenetic tree of 230 accessions from the RDP1 in this study. ADM = admixture; ARO = aromatic; AUS = aus; IND = indica; TEJ = temperate japonica; TRJ = tropic japonica.
Figure 1. Population structure in O. sativa. (A) The large pie chart summarizes the distribution of subpopulations in the 230 O. sativa samples in our diversity panel, and the smaller pie charts on the world map correspond to the country-specific distribution of subpopulations sampled. (B) Phylogenetic tree of 230 accessions from the RDP1 in this study. ADM = admixture; ARO = aromatic; AUS = aus; IND = indica; TEJ = temperate japonica; TRJ = tropic japonica.
Ijms 23 14032 g001
Ijms 23 14032 g002 550
Figure 2. Phenotypic distribution of the 230 cultivars’ panicle blast resistance: (A–C) Diseased grains rate of the 230 cultivars by inoculating isolate 2014-290B27 in 2017, 2018, and 2019 in Nanjing, respectively. (D–F) Distribution of average diseased grains of the population in the six sub-populations.
Figure 2. Phenotypic distribution of the 230 cultivars’ panicle blast resistance: (A–C) Diseased grains rate of the 230 cultivars by inoculating isolate 2014-290B27 in 2017, 2018, and 2019 in Nanjing, respectively. (D–F) Distribution of average diseased grains of the population in the six sub-populations.
Ijms 23 14032 g002
Ijms 23 14032 g003 550
Figure 3. GWAS of the 230 cultivars’ inoculation and the distribution of rice-blast-associated loci. (A–C) Manhattan plots for panicle blast resistance in 2017, 2018, and 2019, respectively. The x-axis shows the single-nucleotide polymorphisms (SNPs) along each chromosome; the y-axis is the −log 10 P for the association. The red arrows indicate the significant SNPs above the threshold, and the black arrows indicate the panicle blast resistance loci. (D) Summary of the identified loci of panicle blast. Each bar represents an associated locus. Blue and red lines indicate the associated loci and repeated loci, respectively.
Figure 3. GWAS of the 230 cultivars’ inoculation and the distribution of rice-blast-associated loci. (A–C) Manhattan plots for panicle blast resistance in 2017, 2018, and 2019, respectively. The x-axis shows the single-nucleotide polymorphisms (SNPs) along each chromosome; the y-axis is the −log 10 P for the association. The red arrows indicate the significant SNPs above the threshold, and the black arrows indicate the panicle blast resistance loci. (D) Summary of the identified loci of panicle blast. Each bar represents an associated locus. Blue and red lines indicate the associated loci and repeated loci, respectively.
Ijms 23 14032 g003
Ijms 23 14032 g004 550
Figure 4. GWAS for leaf blast by inoculating the isolate 2014-290-B27 and LD block analysis for the peak on chromosome 11. (A,B) Manhattan plots from GWAS for rice blast resistance in June 2017 and 2019. The genomic coordinates are displayed along the x-axis and the logarithm of the odds (LOD) score for SNP is displayed on the y-axis. Different colored dots represent significant SNP obtained in different years. The red dashed box indicates the repeat site PBRL3. (C) Local Manhattan plot (top) and LD heatmap (bottom) surrounding the peak on chromosome 11. The red dot indicates the lead SNP at PBRL3.
Figure 4. GWAS for leaf blast by inoculating the isolate 2014-290-B27 and LD block analysis for the peak on chromosome 11. (A,B) Manhattan plots from GWAS for rice blast resistance in June 2017 and 2019. The genomic coordinates are displayed along the x-axis and the logarithm of the odds (LOD) score for SNP is displayed on the y-axis. Different colored dots represent significant SNP obtained in different years. The red dashed box indicates the repeat site PBRL3. (C) Local Manhattan plot (top) and LD heatmap (bottom) surrounding the peak on chromosome 11. The red dot indicates the lead SNP at PBRL3.
Ijms 23 14032 g004
Ijms 23 14032 g005 550
Figure 5. Candidate gene analysis of PBRL3 locus. (A) Haplotype analysis for candidate genes of PBRL3. (B) q-RT PCR analysis of nine potential candidate genes after inoculation. Data represent means ± s.d. (n = 3).
Figure 5. Candidate gene analysis of PBRL3 locus. (A) Haplotype analysis for candidate genes of PBRL3. (B) q-RT PCR analysis of nine potential candidate genes after inoculation. Data represent means ± s.d. (n = 3).
Ijms 23 14032 g005
Ijms 23 14032 g006 550
Figure 6. The haplotype of Pb3 in the natural 230 population of RDP-1. (A) Gene structure of the candidate gene Pb3 (bottom) and conserved domain in Pb3 (top). The blue, red, and green boxes represent different domains of PB3, respectively. The grey boxes represent the coding sequence of Pb3. (B) Haplotype analysis of the Pb3 gene region based on the polymorphic sites shown in (A). (C,D) The distribution of phenotypes corresponding to different haplotypes in 2017 and 2019. (E) Frequency distribution of the five subpopulations in the five haplotypes of Pb3. (F) Nucleotide diversity of chromosome 11. The Pb3 locus is indicated by a green arrow.
Figure 6. The haplotype of Pb3 in the natural 230 population of RDP-1. (A) Gene structure of the candidate gene Pb3 (bottom) and conserved domain in Pb3 (top). The blue, red, and green boxes represent different domains of PB3, respectively. The grey boxes represent the coding sequence of Pb3. (B) Haplotype analysis of the Pb3 gene region based on the polymorphic sites shown in (A). (C,D) The distribution of phenotypes corresponding to different haplotypes in 2017 and 2019. (E) Frequency distribution of the five subpopulations in the five haplotypes of Pb3. (F) Nucleotide diversity of chromosome 11. The Pb3 locus is indicated by a green arrow.
Ijms 23 14032 g006
Ijms 23 14032 g007 550
Figure 7. Evolutionary analysis of Pb3: (A) Phylogenetic tree was constructed based on the protein sequence of R genes by MAGE7.0. Pb3 is shaded in red. (B) The gene structure of R genes. Conserved domains were predicted via the Pfam online site. Amino acid sequences are represented by black lines and domains are represented by colorful boxes.
Figure 7. Evolutionary analysis of Pb3: (A) Phylogenetic tree was constructed based on the protein sequence of R genes by MAGE7.0. Pb3 is shaded in red. (B) The gene structure of R genes. Conserved domains were predicted via the Pfam online site. Amino acid sequences are represented by black lines and domains are represented by colorful boxes.
Ijms 23 14032 g007
Ijms 23 14032 g008 550
Figure 8. The phenotype of T-DNA insertion mutation lines of Pb3. (A) Schematic presentation of the gene structure of Pb3 and T-DNA insertion sites. (B,C) The average lesion number and length in of wild-type Dongjin and five mutants inoculated with RO1-1 at the seedling stage. (D) Percentage of diseased grains of wild-type Dongjin and four mutant lines. Each material is inoculated with 9 panicles. (E) Spray inoculation of the four mutants and wild-type plants using strain RO1-1. (F) Injection inoculation at the booting stage of the four mutants and wild-type plants using strain RO1-1. * p < 0.05, ** p < 0.01.
Figure 8. The phenotype of T-DNA insertion mutation lines of Pb3. (A) Schematic presentation of the gene structure of Pb3 and T-DNA insertion sites. (B,C) The average lesion number and length in of wild-type Dongjin and five mutants inoculated with RO1-1 at the seedling stage. (D) Percentage of diseased grains of wild-type Dongjin and four mutant lines. Each material is inoculated with 9 panicles. (E) Spray inoculation of the four mutants and wild-type plants using strain RO1-1. (F) Injection inoculation at the booting stage of the four mutants and wild-type plants using strain RO1-1. * p < 0.05, ** p < 0.01.
Ijms 23 14032 g008
Ijms 23 14032 g009 550
Figure 9. The Phenotype of CRISPR lines of Pb3. (A) Schematic presentation of the gene structure of Pb3 and CRISPR/CAS9 mutation sites. (B) Mutations in the three knockout lines of Pb3 in the Bodao background. The red arrow indicates the position of the CDs of Pb3. The bases marked red represent the inserted sequence in the genome. The red line indicates the missing sequence on the genome. (C) Percentage of diseased grains of the wild-type Bodao and three mutant lines. Different lines are distinguished by different colors. Each material is inoculated with 9 panicles. (D) Injection inoculation at the booting stage of the six mutants and wild-type plants using strain 2015-248A3. BO: Bodao. * p < 0.05.
Figure 9. The Phenotype of CRISPR lines of Pb3. (A) Schematic presentation of the gene structure of Pb3 and CRISPR/CAS9 mutation sites. (B) Mutations in the three knockout lines of Pb3 in the Bodao background. The red arrow indicates the position of the CDs of Pb3. The bases marked red represent the inserted sequence in the genome. The red line indicates the missing sequence on the genome. (C) Percentage of diseased grains of the wild-type Bodao and three mutant lines. Different lines are distinguished by different colors. Each material is inoculated with 9 panicles. (D) Injection inoculation at the booting stage of the six mutants and wild-type plants using strain 2015-248A3. BO: Bodao. * p < 0.05.
Ijms 23 14032 g009
Table
Table 1. The PBRLs that were associated with panicle blast in three experiments by inoculating 2014-290-B27.
Table 1. The PBRLs that were associated with panicle blast in three experiments by inoculating 2014-290-B27.
LocusChrPositionTop SNPp-ValueYearsLocus Reference
PBRL1233611914–33651599336119144.95 × 10−52017Pig(t)
PBRL2104834530–484426548345307.73 × 10−52017
PBRL31127114833–27391298271183703.58 × 10−62017/2019Pik
PBRL41226529664–26995595269542477.56 × 10−52017
PBRL5140778316–41195405411954056.32 × 10−62018LABR_19
PBRL6420342333–20624109203423338.79 × 10−72018PBRL-14 Pi-21
PBRL7613298661–13477930134052161.88 × 10−62018
PBRL8620225421–20855898203671248.57 × 10−62018
PBRL9621936925–22049270219996471.59 × 10−52018LABR_49
PBRL10103525543–363129535829611.45 × 10−72018
PBRL11125009847–25590336254720016.05 × 10−72019 LABR_11
PBRL12321526136–22360496215484376.22 × 10−72019LABR_33
PBRL13325404248–25800205254042489.70 × 10−72019LABR_34
PBRL14613982991–14520587142040453.33 × 10−6 2019
PBRL15910788559–10895765108098433.13 × 10−6 2019LABR_67
PBRL161015103798–15256763152515372.26 × 10−62019
Table
Table 2. The candidate genes of PBRL3 that were significantly different by haplotype analysis.
Table 2. The candidate genes of PBRL3 that were significantly different by haplotype analysis.
LocusChrGene IDSNP LocationGene AnnotationHaplotype Analysis
PBRL311LOC_Os11g44890CDsexpressed proteinp < 0.01
PBRL311LOC_Os11g44910Promoter CDsDEAD-box ATP-dependent RNA helicase, putative, expressedp < 0.05
PBRL311LOC_Os11g44930CDspentatricopeptide repeat domain containing protein, putative, expressedp < 0.05
PBRL311LOC_Os11g44950Promoter CDsglycosyl hydrolase family 3 protein, putative, expressedp < 0.05
PBRL311LOC_Os11g44990CDsNB-ARC domain containing protein, expressedp < 0.05
PBRL311LOC_Os11g45030CDsexpressed proteinp < 0.05
PBRL311LOC_Os11g45090Promoter CDsNB-ARC domain containing protein, expressedp < 0.05
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Ma, L.; Yu, Y.; Li, C.; Wang, P.; Liu, K.; Ma, W.; Wang, W.; Fan, Y.; Xiong, Z.; Jiang, T.; et al. Genome-Wide Association Study Identifies a Rice Panicle Blast Resistance Gene Pb3 Encoding NLR Protein. Int. J. Mol. Sci. 2022, 23, 14032. https://doi.org/10.3390/ijms232214032

AMA Style

Ma L, Yu Y, Li C, Wang P, Liu K, Ma W, Wang W, Fan Y, Xiong Z, Jiang T, et al. Genome-Wide Association Study Identifies a Rice Panicle Blast Resistance Gene Pb3 Encoding NLR Protein. International Journal of Molecular Sciences. 2022; 23(22):14032. https://doi.org/10.3390/ijms232214032

Chicago/Turabian Style

Ma, Lu, Yao Yu, Changqing Li, Panting Wang, Kunquan Liu, Wenjing Ma, Wei Wang, Yunxin Fan, Ziwei Xiong, Tingting Jiang, and et al. 2022. "Genome-Wide Association Study Identifies a Rice Panicle Blast Resistance Gene Pb3 Encoding NLR Protein" International Journal of Molecular Sciences 23, no. 22: 14032. https://doi.org/10.3390/ijms232214032

APA Style

Ma, L., Yu, Y., Li, C., Wang, P., Liu, K., Ma, W., Wang, W., Fan, Y., Xiong, Z., Jiang, T., Zhang, J., Wang, Z., Wang, J., Zhang, H., & Bao, Y. (2022). Genome-Wide Association Study Identifies a Rice Panicle Blast Resistance Gene Pb3 Encoding NLR Protein. International Journal of Molecular Sciences, 23(22), 14032. https://doi.org/10.3390/ijms232214032

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