Mapping quantitative trait Loci for resistance to rice blast

Refinement of rice blast disease resistance QTLs and gene networks through meta-QTL analysis

Rice blast disease is the most devastating disease constraining crop productivity. Vertical resistance to blast disease is widely studied despite its instability. Clusters of genes or QTLs conferring blast resistance that offer durable horizontal resistance are important in resistance breeding. In this study, we aimed to refine the reported QTLs and identify stable meta-QTLs (MQTLs) associated with rice blast resistance. A total of 435 QTLs were used to project 71 MQTLs across all the rice chromosomes. As many as 199 putative rice blast resistance genes were identified within 53 MQTL regions. The genes included 48 characterized resistance gene analogs and related proteins, such as NBS-LRR type, LRR receptor-like kinase, NB-ARC domain, pathogenesis-related TF/ERF domain, elicitor-induced defense and proteins involved in defense signaling. MQTL regions with clusters of RGA were also identified. Fifteen highly significant MQTLs included 29 candidate genes and genes characterized for blast resistance, such as Piz, Nbs-Pi9, pi55-1, pi55-2, Pi3/Pi5-1, Pi3/Pi5-2, Pikh, Pi54, Pik/Pikm/Pikp, Pb1 and Pb2. Furthermore, the candidate genes (42) were associated with differential expression (in silico) in compatible and incompatible reactions upon disease infection. Moreover, nearly half of the genes within the MQTL regions were orthologous to those in O. sativa indica, Z. mays and A. thaliana, which confirmed their significance. The peak markers within three significant MQTLs differentiated blast-resistant and susceptible lines and serve as potential surrogates for the selection of blast-resistant lines. These MQTLs are potential candidates for durable and broad-spectrum rice blast resistance and could be utilized in blast resistance breeding.

Advancement in the Breeding, Biotechnological and Genomic Tools towards Development of Durable Genetic Resistance against the Rice Blast Disease

Rice production needs to be sustained in the coming decades, as the changeable climatic conditions are becoming more conducive to disease outbreaks. The majority of rice diseases cause enormous economic damage and yield instability. Among them, rice blast caused by Magnaportheoryzae is a serious fungal disease and is considered one of the major threats to world rice production. This pathogen can infect the above-ground tissues of rice plants at any growth stage and causes complete crop failure under favorable conditions. Therefore, management of blast disease is essentially required to sustain global food production. When looking at the drawback of chemical management strategy, the development of durable, resistant varieties is one of the most sustainable, economic, and environment-friendly approaches to counter the outbreaks of rice blasts. Interestingly, several blast-resistant rice cultivars have been developed with the help of breeding and biotechnological methods. In addition, 146 R genes have been identified, and 37 among them have been molecularly characterized to date. Further, more than 500 loci have been identified for blast resistance which enhances the resources for developing blast resistance through marker-assisted selection (MAS), marker-assisted backcross breeding (MABB), and genome editing tools. Apart from these, a better understanding of rice blast pathogens, the infection process of the pathogen, and the genetics of the immune response of the host plant are very important for the effective management of the blast disease. Further, high throughput phenotyping and disease screening protocols have played significant roles in easy comprehension of the mechanism of disease spread. The present review critically emphasizes the pathogenesis, pathogenomics, screening techniques, traditional and molecular breeding approaches, and transgenic and genome editing tools to develop a broad spectrum and durable resistance against blast disease in rice. The updated and comprehensive information presented in this review would be definitely helpful for the researchers, breeders, and students in the planning and execution of a resistance breeding program in rice against this pathogen.

Mapping Blast Resistance Genes in Rice Varieties 'Minghui 63' and 'M-202'

Rice blast caused by the fungus Magnaporthe oryzae (syn. Magnaporthe grisea) is one of the most lethal diseases for sustainable rice production worldwide. Blast resistance mediated by major resistance genes is often broken down after a short period of deployment, while minor blast resistance genes, each providing a small effect on disease reactions, are more durable. In the present study, we first evaluated disease reactions of two rice breeding parents 'Minghui 63' and 'M-202' with 11 blast races, IA45, IB1, IB45, IB49, IB54, IC1, IC17, ID1, IE1, IG1, and IH1, commonly present in the United States, under greenhouse conditions using a category disease rating resembling infection types under field conditions. 'Minghui 63' exhibited differential resistance responses in comparison with those of 'M-202' to the tested blast races. A recombinant inbred line (RIL) population of 275 lines from a cross between 'Minghui 63' and 'M-202' was also evaluated with the above-mentioned blast races. The population was genotyped with 156 simple sequence repeat (SSR) and insertion and deletion (Indel) markers. A linkage map with a genetic distance of 1,022.84 cM was constructed using inclusive composite interval mapping (ICIM) software. A total of 10 resistance QTLs, eight from 'Minghui 63' and two from 'M-202', were identified. One major QTL, qBLAST2 on chromosome 2, was identified by seven races/isolates. The remaining nine minor resistance QTLs were mapped on chromosomes 1, 3, 6, 9, 10, 11, and 12. These findings provide useful genetic markers and resources to tag minor blast resistance genes for marker-assisted selection in rice breeding program and for further studies of underlying genes.

A Meta-Analysis of Quantitative Trait Loci Associated with Multiple Disease Resistance in Rice (Oryza sativa L.)

Rice blast, sheath blight and bacterial leaf blight are major rice diseases found worldwide. The development of resistant cultivars is generally perceived as the most effective way to combat these diseases. Plant disease resistance is a polygenic trait where a combinatorial effect of major and minor genes affects this trait. To locate the source of this trait, various quantitative trait loci (QTL) mapping studies have been performed in the past two decades. However, investigating the congruency between the reported QTL is a daunting task due to the heterogeneity amongst the QTLs studied. Hence, the aim of our study is to integrate the reported QTLs for resistance against rice blast, sheath blight and bacterial leaf blight and objectively analyze and consolidate the location of QTL clusters in the chromosomes, reducing the QTL intervals and thus identifying candidate genes within the selected meta-QTL. A total of twenty-seven studies for resistance QTLs to rice blast (8), sheath blight (15) and bacterial leaf blight (4) was compiled for QTL projection and analyses. Cumulatively, 333 QTLs associated with rice blast (114), sheath blight (151) and bacterial leaf blight (68) resistance were compiled, where 303 QTLs could be projected onto a consensus map saturated with 7633 loci. Meta-QTL analysis on 294 QTLs yielded 48 meta-QTLs, where QTLs with membership probability lower than 60% were excluded, reducing the number of QTLs within the meta-QTL to 274. Further, three meta-QTL regions (MQTL2.5, MQTL8.1 and MQTL9.1) were selected for functional analysis on the basis that MQTL2.5 harbors the highest number of QTLs; meanwhile, MQTL8.1 and MQTL9.1 have QTLs associated with all three diseases mentioned above. The functional analysis allows for determination of enriched gene ontology and resistance gene analogs (RGAs) and other defense-related genes. To summarize, MQTL2.5, MQTL8.1 and MQTL9.1 have a considerable number of R-genes that account for 10.21%, 4.08% and 6.42% of the total genes found in these meta-QTLs, respectively. Defense genes constitute around 3.70%, 8.16% and 6.42% of the total number of genes in MQTL2.5, MQTL8.1 and MQTL9.1, respectively. This frequency is higher than the total frequency of defense genes in the rice genome, which is 0.0096% (167 defense genes/17,272 total genes). The integration of the QTLs facilitates the identification of QTL hotspots for rice blast, sheath blight and bacterial blight resistance with reduced intervals, which helps to reduce linkage drag in breeding. The candidate genes within the promising regions could be utilized for improvement through genetical engineering.

Identification of Novel QTL Conferring Sheath Blight Resistance in Two Weedy Rice Mapping Populations

BACKGROUND

Rice sheath blight (ShB) disease, caused by the pathogenic fungus Rhizoctonia solani, causes significant yield losses globally. US weedy rice populations, which are de-domesticated forms of indica and aus cultivated rice, appear to be more resistant to ShB than local japonica cultivated rice. We mapped quantitative trait loci (QTL) associated with ShB resistance using two F₈ recombinant inbred line populations generated from crosses of an indica crop variety, Dee-Geo-Woo-Gen (DGWG), with individuals representing the two major US weed biotypes, straw hull (SH) and black hull awned (BHA).

RESULTS

We identified nine ShB resistance QTL across both mapping populations. Five were attributable to alleles that affect plant height (PH) and heading date (HD), two growth traits that are known to be highly correlated with ShB resistance. By utilizing an approach that treated growth traits as covariates in the mapping model, we were able to infer that the remaining four QTL are involved in ShB resistance. Two of these, qShB1-2 and qShB4, are different from previously identified ShB QTL and represent new candidates for further study.

CONCLUSION

Our findings suggest that ShB resistance can be improved through favorable plant growth traits and the combined effects of small to moderate-effect resistance QTL. Additionally, we show that including PH and HD as covariates in QTL mapping models is a powerful way to identify new ShB resistance QTL.

Pi5 and Pii Paired NLRs Are Functionally Exchangeable and Confer Similar Disease Resistance Specificity

Effector-triggered immunity (ETI) is an effective layer of plant defense initiated upon recognition of avirulence (Avr) effectors from pathogens by cognate plant disease resistance (R) proteins. In rice, a large number of R genes have been characterized from various cultivars and have greatly contributed to breeding programs to improve resistance against the rice blast pathogen Magnaporthe oryzae. The extreme diversity of R gene repertoires is thought to be a result of co-evolutionary history between rice and its pathogens including M. oryzae. Here we show that Pii is an allele of Pi5 by DNA sequence characterization and complementation analysis. Pii-1 and Pii-2 cDNAs were cloned by reverse transcription polymerase chain reaction from the Pii -carrying cultivar Fujisaka5 . The complementation test in susceptible rice cultivar Dongjin demonstrated that the rice blast resistance mediated by Pii , similar to Pi5 , requires the presence of two nucleotide-binding leucine-rich repeat genes, Pii-1 and Pii-2 . Consistent with our hypothesis that Pi5 and Pii are functionally indistinguishable, the replacement of Pii-1 by Pi5-1 and Pii-2 by Pi5-2 , respectively, does not change the level of disease resistance to M. oryzae carrying AVR-Pii. Surprisingly, Exo70F3, required for Pii-mediated resistance, is dispensable for Pi5-mediated resistance. Based on our results, despite similarities observed between Pi5 and Pii, we hypothesize that Pi5 and Pii pairs require partially distinct mechanisms to function.

Functional Interactions Between Major Rice Blast Resistance Genes, Pi-ta and Pi-b, and Minor Blast Resistance Quantitative Trait Loci

Major blast resistance (R) genes confer resistance in a gene-for-gene manner. However, little information is available on interactions between R genes. In this study, interactions between two rice blast R genes, Pi-ta and Pi-b, and other minor blast resistance quantitative trait loci (QTLs) were investigated in a recombinant inbred line (RIL) population comprising 243 RILs from a Cybonnet (CYBT) × Saber (SB) cross. CYBT has the R gene Pi-ta and SB has Pi-b. Ten differential isolates of four Magnaporthe oryzae races (IB-1, IB-17, IB-49, and IE-1K) were used to evaluate disease reactions of the 243 RILs under greenhouse conditions. Five resistance QTLs were mapped on chromosomes 2, 3, 8, 9, and 12 with a linkage map of 179 single nucleotide polymorphism markers. Among them, qBR12 (Q1), was mapped at the Pi-ta locus and accounted for 45.41% of phenotypic variation while qBR2 (Q2) was located at the Pi-b locus and accounted for 24.81% of disease reactions. The additive-by-additive epistatic interaction between Q1 (Pi-ta) and Q2 (Pi-b) was detected; they can enhance the disease resistance by an additive 0.93 using the 0 to 9 standard phenotyping method. These results suggest that Pi-ta interacts synergistically with Pi-b.

Joint parameter estimation in the QTL mapping of ordinal traits

With the rapid development of statistical genetics, the deep researches of ordinal traits have been gradually emphasized. The data of these traits bear relatively less information than those of continuous phenotypes, therefore it is more complex to map the quantitative trait loci (QTL) of ordinal traits. In this paper, the multiple-interval mapping method is considered in the genetic mapping of ordinal traits. By combining threshold model and statistical model, we build a cumulative logistic regression model to express the relationship between the ordinal data and the QTL genotypes. In order to make the interval mapping more straightforward, we treat the recombination rates as unknown parameters, and then simultaneously obtain the estimates of QTL positions, QTL effects and threshold parameters via the EM algorithm. We perform simulation experiments to investigate and compare the proposed method. We also present a real example to test the reasonableness of the considered model and estimate both model parameters and QTL parameters. Both results of simulations and example show that the method we proposed is reasonable and effective.

Dissection of the genetic architecture of rice resistance to the blast fungus Magnaporthe oryzae

Resistance in rice cultivars to the rice blast fungus Magnaporthe oryzae is complex and is controlled by both major genes and quantitative trait loci (QTLs). We undertook a genome-wide association study (GWAS) using the rice diversity panel 1 (RDP1) that was genotyped using a high-density (700 000 single nucleotide polymorphisms) array and inoculated with five diverse M. oryzae isolates. We identified 97 loci associated with blast resistance (LABRs). Among them, 82 were new regions and 15 co-localized with known blast resistance loci. The top 72 LABRs explained up to 98% of the phenotypic variation. The candidate genes in the LABRs encode nucleotide-binding site leucine-rich repeat (NBS-LRR) resistance proteins, receptor-like protein kinases, transcription factors and defence-related proteins. Among them, LABR_64 was strongly associated with resistance to all five isolates. We analysed the function of candidate genes underlying LABR_64 using RNA interference (RNAi) technology and identified two new resistance alleles at the Pi5 locus. We demonstrate an efficient strategy for rapid allele discovery using the power of GWAS, coupled with RNAi technology, for the dissection of complex blast resistance in rice.

Coevolutionary Dynamics of Rice Blast Resistance Gene Pi-ta and Magnaporthe oryzae Avirulence Gene AVR-Pita 1

The Pi-ta gene in rice is effective in preventing infections by Magnaporthe oryzae strains that contain the corresponding avirulence gene, AVR-Pita1. Diverse haplotypes of AVR-Pita1 have been identified from isolates of M. oryzae from rice production areas in the United States and worldwide. DNA sequencing and mapping studies have revealed that AVR-Pita1 is highly unstable, while expression analysis and quantitative resistance loci mapping of the Pi-ta locus revealed complex evolutionary mechanisms of Pi-ta-mediated resistance. Among these studies, several Pi-ta transcripts were identified, most of which are probably derived from alternative splicing and exon skipping, which could produce functional resistance proteins that support a new concept of coevolution of Pi-ta and AVR-Pita1. User-friendly DNA markers for Pi-ta have been developed to support marker-assisted selection, and development of new rice varieties with the Pi-ta markers. Genome-wide association studies revealed a link between Pi-ta-mediated resistance and yield components suggesting that rice has evolved a complicated defense mechanism against the blast fungus. In this review, we detail the current understanding of Pi-ta allelic variation, its linkage with rice productivity, AVR-Pita allelic variation, and the coevolution of Pi-ta and AVR-Pita in Oryza species and M. oryzae populations, respectively. We also review the genetic and molecular basis of Pi-ta and AVR-Pita interaction, and its value in marker-assisted selection and engineering resistance.

QTL Analysis for Resistance to Blast Disease in U.S. Weedy Rice

Understanding the genetic architecture of adaptation is of great importance in evolutionary biology. U.S. weedy rice is well adapted to the local conditions in U.S. rice fields. Rice blast disease is one of the most destructive diseases of cultivated rice worldwide. However, information about resistance to blast in weedy rice is limited. Here, we evaluated the disease reactions of 60 U.S. weedy rice accessions with 14 blast races, and investigated the quantitative trait loci (QTL) associated with blast resistance in two major ecotypes of U.S. weedy rice. Our results revealed that U.S. weedy rice exhibited a broad resistance spectrum. Using genotyping by sequencing, we identified 28 resistance QTL in two U.S. weedy rice ecotypes. The resistance QTL with relatively large and small effects suggest that U.S. weedy rice groups have adapted to blast disease using two methods, both major resistance (R) genes and QTL. Three genomic loci shared by some of the resistance QTL indicated that these loci may contribute to no-race-specific resistance in weedy rice. Comparing with known blast disease R genes, we found that the R genes at these resistance QTL are novel, suggesting that U.S. weedy rice is a potential source of novel blast R genes for resistant breeding.

Combination Patterns of Major R Genes Determine the Level of Resistance to the M. oryzae in Rice (Oryza sativa L.)

Rice blast caused by Magnaporthe oryzae is the most devastating disease of rice and poses a serious threat to world food security. In this study, the distribution and effectiveness of 18 R genes in 277 accessions were investigated based on pathogenicity assays and molecular markers. The results showed that most of the accessions exhibited some degree of resistance (resistance frequency, RF >50%). Accordingly, most of the accessions were observed to harbor two or more R genes, and the number of R genes harbored in accessions was significantly positively correlated with RF. Some R genes were demonstrated to be specifically distributed in the genomes of rice sub-species, such as Pigm, Pi9, Pi5 and Pi1, which were only detected in indica-type accessions, and Pik and Piz, which were just harbored in japonica-type accessions. By analyzing the relationship between R genes and RF using a multiple stepwise regression model, the R genes Pid3, Pi5, Pi9, Pi54, Pigm and Pit were found to show the main effects against M. oryzae in indica-type accessions, while Pita, Pb1, Pik, Pizt and Pia were indicated to exhibit the main effects against M. oryzae in japonica-type accessions. Principal component analysis (PCA) and cluster analysis revealed that combination patterns of major R genes were the main factors determining the resistance of rice varieties to M. oryzae, such as 'Pi9+Pi54', 'Pid3+Pigm', 'Pi5+Pid3+Pigm', 'Pi5+Pi54+Pid3+Pigm', 'Pi5+Pid3' and 'Pi5+Pit+Pid3' in indica-type accessions and 'Pik+Pib', 'Pik+Pita', 'Pik+Pb1', 'Pizt+Pia' and 'Pizt+Pita' in japonica-type accessions, which were able to confer effective resistance against M. oryzae. The above results provide good theoretical support for the rational utilization of combinations of major R genes in developing rice cultivars with broad-spectrum resistance.

A multiple gene complex on rice chromosome 4 is involved in durable resistance to rice blast

Quantitative trait loci (QTLs) for resistance to rice blast offer a potential source of durable disease resistance in rice. However, few QTLs have been validated in progeny testing, on account of their small phenotypic effects. To understand the genetic basis for QTL-mediated resistance to blast, we dissected a resistance QTL, qBR4-2, using advanced backcross progeny derived from a chromosome segment substitution line in which a 30- to 34-Mb region of chromosome 4 from the resistant cultivar Owarihatamochi was substituted into the genetic background of the highly susceptible Aichiasahi. The analysis resolved qBR4-2 into three loci, designated qBR4-2a, qBR4-2b, and qBR4-2c. The sequences of qBR4-2a and qBR4-2b, which lie 181 kb apart from each other and measure, 113 and 32 kb, respectively, appear to encode proteins with a putative nucleotide-binding site (NBS) and leucine-rich repeats (LRRs). Sequence analysis of the donor allele of qBR4-2a, the region with the largest effect among the three, revealed sequence variations in the NBS-LRR region. The effect of qBR4-2c was smallest among the three, but its combination with the donor alleles of qBR4-2a and qBR4-2b significantly enhanced blast resistance. qBR4-2 comprises three tightly linked QTLs that control blast resistance in a complex manner, and thus gene pyramiding or haplotype selection is the recommended strategy for improving QTL-mediated resistance to blast disease through the use of this chromosomal region.

Molecular mapping of the blast resistance genes Pi2-1 and Pi51(t) in the durably resistant rice 'Tianjingyeshengdao'

Tianjingyeshengdao' (TY) is a rice cultivar with durable resistance to populations of Magnaporthe oryzae (the causal agent of blast) in China. To understand the genetic basis of its resistance to blast, we developed a population of recombinant inbred lines from a cross between TY and the highly susceptible 'CO39' for gene mapping analysis. In total, 22 quantitative trait loci (QTLs) controlling rice blast resistance were identified on chromosomes 1, 3, 4, 5, 6, 9, 11, and 12 from the evaluation of four disease parameters in both greenhouse and blast nursery conditions. Among these QTLs, 19 were contributed by TY and three by CO39. Two QTL clusters on chromosome 6 and 12 were named Pi2-1 and Pi51(t), respectively. Pi2-1 was detected under both growth chamber and natural blast nursery conditions, and explained 31.24 to 59.73% of the phenotypic variation. Pi51(t) was only detected in the natural blast nursery and explained 3.67 to 10.37% of the phenotypic variation. Our results demonstrate that the durable resistance in TY is controlled by two major and seven minor genes. Identification of the markers linked to both Pi2-1 and Pi51(t) in this study should be useful for marker-aided selection in rice breeding programs as well as for molecular cloning of the identified resistance genes.