利用二代测序技术高通量分析水稻品种间抗病性差异的分子基础

稻瘟病和白叶枯病是由稻温病菌(Magnaporthe oryzae)和白叶枯病菌(Xanthomonas oryzae pv.oryzae)引起的两种主要水稻病害,也是制约中国水稻生产的主要病害。为了从DNA水平探索造成水稻感病品种‘丽江新团黑谷’(LTH)和高抗品种‘特特普’(Tetep,TTP)间抗病性差异的分子基础,该研究对其已知的3个抗稻瘟病基因和3个抗白叶枯基因所在DNA区段分别进行PCR扩增,将等量混合的PCR产物再与基因组重测序样品按Ct值差值(ΔCt)~10的比例混合,采用二代测序技术进行一次性测序和比较分析,并对有差异的基因区域进行常规传统测序验证,以确定这2个品种中抗性基因(R基因)的数目和结构与品种抗病或感病表型的关联性。实验结果表明,二代测序能够快速并准确地寻找到2个不同水稻品种中多个特定基因的序列差异,且差异位点与常规测序结果相符。从LTH和TTP这2种抗性不同水稻品种在多个抗性基因的DNA水平差异来看,有差异的抗性基因位点在高抗品种TTP中大都与原始抗性基因序列相同,而对应的普感品种LTH的抗性基因往往多表现为氨基酸突变,这些序列差异很可能就是导致TTP与LTH抗性差异的分子基础。     

Molecular Mapping of Rice Blast Resistance Gene Pi-k h in the Rice Variety Tetep

We have used rice line Tetep as a resistant donor with the aim of mapping a durable blast resistance gene Pi-k ^h using RAPD and AFLP techniques in conjunction with bulk segregant analysis. An F₂ mapping population consisting of 205 plants was generated by crossing Tetep with HP2216, a highly susceptible cultivar. Inoculation with specific isolate (PLP-1) of Magnaporthe grisea at seeding stage showed that the Pi-k ^h gene inherited as a single dominant gene in F2 population. RAPD analysis was performed with 240 primers to detect polymorphism between resistant and susceptible parents. Of these, 48 primers produced polymorphic banding pattern between resistant and susceptible parents. Bulk segregant analysis was performed with 48 primers of which 5 showed polymorphism between resistant and susceptible bulks. A 700 bp DNA band was obtained in resistant F2 plants with primer 5-129 indicating its linkage to the resistance gene. Out of 64 AFLP primer combinations used for polymorphism survey between HP 2216 and Tetep, 11 AFLP primer combinations were able to distinguish the resistant and susceptible bulks. An AFLP band of 75 bp obtained with primer combination, E-TAlM-CTC co-segregated with the resistance gene. The RAPD marker 5-129₇₀₀ and AFLP₇₅ were placed on the linkage map at a distance of 2.1 eM and 15.1 eM flanking to Pi-k ^hgene, respectively. The RAPD band closely linked to Pi-k ^h gene was sequenced and used for the development of CAPs markers which also co-segregated with resistant phenotype in the mapping population. On sequence analysis and homology search of RAPD fragment with whole rice genome sequence database and the information available on physical, genetic and sequence maps of rice, the co-segregating CAPs marker was placed at long arm of rice chromosome 11. CAPs marker developed in this study showed polymorphism in different rice cultivars grown in North-Western Himalayan region and is being used for the pyramiding of Pi-k ^h gene along with other blast resistance genes using marker-assisted selection.     

Allele-mining of rice blast resistance genes at AC134922 locus

The AC134922 locus is one of the most rapidly evolving nucleotide binding site-leucine-rich repeat (NBS-LRR) gene family in rice genome. Six rice blast resistance (R) genes have been cloned from this locus and other two resistance candidate genes, Pi34 and Pi47, are also mapped to this complex locus. Therefore, it seems that more functional R genes could be identified from this locus. In this study, we cloned 22 genes from 12 cultivars based on allele-mining strategy at this locus and identified 6 rice blast R genes with 4 of them recognizing more than one isolates. Our result suggests that gene stacking might be the evolutionary strategy for complex gene locus to interact with rapidly evolving pathogens, which might provide a potential way for the cloning of durable resistance genes. Moreover, the mosaic structure and ambiguous ortholog/paralog relationships of these homologous genes, caused by frequent recombination and gene conversion, indicate that multiple alleles of this complex locus may serve as a reservoir for the evolutionary novelty of these R genes.     

Functional complementation of rice blast resistance gene Pi-k h (Pi54) conferring resistance to diverse strains of Magnaporthe oryzae

Blast disease of rice, caused by Magnaporthe oryzae is an explosive disease that can spread rapidly in conducive conditions. R-gene mediated resistance offers an environmentally sustainable solution for management of this important disease of rice. We have earlier identified a unique R-gene of rice, on chromosome 11 of Oryza sativa ssp. indica cultivar Tetep. In this study we report functional validation of the Pi-k ^h (Pi54) gene using complementation assay. The blast resistance candidate gene Pi-k ^h (Pi54) was cloned into a plant transformation vector and the construct was used to transform a japonica cultivar of rice Taipei 309, which is susceptible to M. oryzae. Transgenic lines containing Pi-k ^h (Pi54) gene were found to confer high degree of resistance to diverse isolates of M. oryzae. The callose deposition was analyzed and compared between the transgenic and non-transgenic rice plants and widespread deposition was observed at the infection sites in plants showing incompatible interaction. Successful complementation of Pi-k ^h (Pi54) gene confirmed that the gene is responsible for resistance to M. oryzae in transgenic lines developed during this study. Expression analysis of the gene in resistant plants revealed that the gene is pathogen inducible in nature and is not expressed constitutively. Detection of callose deposition in resistant plants containing Pi-k ^h (Pi54) gene implicates its involvement in the initiation of defense response cascade.     

Transposon-based high sequence diversity in Avr-Pita alleles increases the potential for pathogenicity of Magnaporthe oryzae populations

Magnaporthe oryzae causes rice blast that is one of the most devastating diseases of rice worldwide. Highly variable nature of this fungus has evolved itself against major resistance genes in newly released rice varieties. Understanding the population structure of this fungus is essential for proper utilization of the rice blast resistance genes in rice crop plants. In the present study, we analyzed 133 isolates of M. oryzae from ten countries to find the allelic variation of Avr-Pita gene that is triggering Pita-mediated resistance in rice plant. The diversity analysis of these alleles showed higher level of nucleotide variation in the coding regions than the noncoding regions. Evolutionary analysis of these alleles indicates that Avr-Pita gene is under purifying selection to favor its major alleles in 133 isolates analyzed in this study. We hypothesize that the selection of favorable Avr-Pita allele in these isolates may occur through a genetic mechanism known as recurrent selective sweeps. A total of 22 functional Avr-Pita protein variants were identified in this study. Insertion of Pot3 transposable element into the promoter of Avr-Pita gene was identified in virulent isolates and was suggested that mobility of repeat elements in avirulence genes of M. oryzae seems to help in emergence of new virulent types of the pathogen. Allele-specific markers developed in this study will be helpful to identify a particular type of Avr-Pita allele from M. oryzae population which can form the basis for the deployment of Pita gene in different epidemiological regions.     

Molecular diversity in rice blast resistance gene Pi-ta makes it highly effective against dynamic population of Magnaporthe oryzae

Rice blast is one of the important diseases of rice which can be effectively managed by the deployment of resistance genes. Pi-ta is one of the major blast resistant genes effective against pathogen populations in different parts of India. We analysed allelic variants of Pi-ta from 48 rice lines selected after phenotyping of 529 rice landraces across three eco-geographical blast hot spot regions. Besides, Pi-ta orthologue sequences of 220 rice accessions belonging to wild and cultivated species of rice were also included in the study for a better evo–devo perspective of the diversity present in the gene and the selection pressures acting on this locus. We obtained high nucleotide variations (SNPs and insertion–deletions) in the intronic region. We also identified 64 haplotypes based on nucleotide polymorphism in these alleles. Pi-ta orthologues of Indian landraces were scattered in eight major haplotypes indicating its heterogenous nature. We identified a total of 47 different Pi-ta protein variants on the basis of deduced amino acid residues amongst the orthologues. Five unique and novel Pi-ta variants were identified for the first time in rice landraces exhibiting different reaction types against the Magnaporthe oryzae population. A high value of Pi_non/syn was observed only in the leucine-rich domain of the alleles cloned from Indian landraces, indicating strong selective forces acting on this region. The detailed molecular analysis of the Pi-ta orthologues provides insights to a high degree of inter- and intraspecific relationships amongst the Oryza species. We identified rice landraces possessing the effective alleles of this resistance gene which can be used in future blast resistance breeding programmes.     

High-resolution mapping, cloning and molecular characterization of the Pi-k h gene of rice, which confers resistance to Magnaporthe grisea

In order to understand the molecular mechanisms involved in the gene-for-gene type of pathogen resistance, high-resolution genetic and physical mapping of resistance loci is required to facilitate map-based cloning of resistance genes. Here, we report the molecular mapping and cloning of a dominant gene (Pi-k ^h ) present in the rice line Tetep, which is associated with resistance to rice blast disease caused by Magnaporthe grisea. This gene is effective against M. grisea populations prevalent in the Northwestern Himalayan region of India. Using 178 sequence tagged microsatellite, sequence-tagged site, expressed sequence tag and simple sequence repeat (SSR) markers to genotype a population of 208 F₂ individuals, we mapped the Pi-k ^h gene between two SSR markers (TRS26 and TRS33) which are 0.7 and 0.5 cM away, respectively, and can be used in marker-assisted-selection for blast-resistant rice cultivars. We used the markers to identify the homologous region in the genomic sequence of Oryza sativa cv. Nipponbare, and a physical map consisting of two overlapping bacterial artificial chromosome and P1 artificial chromosome clones was assembled, spanning a region of 143,537 bp on the long arm of chromosome 11. Using bioinformatic analyses, we then identified a candidate blast-resistance gene in the region, and cloned the homologous sequence from Tetep. The putative Pi-k ^h gene cloned from Tetep is 1.5 kbp long with a single ORF, and belongs to the nucleotide binding site-leucine rich repeat class of disease resistance genes. Structural and expression analysis of the Pi-k ^h gene revealed that its expression is pathogen inducible.     

Sequence variation of avirulence gene AVR-Pita1 in rice blast fungus, Magnaporthe oryzae

The interaction between rice, Oryza sativa, and rice blast fungus, Magnaporthe oryzae, is triggered by an interaction between the protein products of the host resistant gene, and the pathogen avirulence gene. This interaction follows the ‘gene-for-gene' concept. The resistant gene has effectively protected rice plants from rice blast infection. However, the resistant genes usually break down several years after the release of the resistant rice varieties because the fungus has evolved to new races. The objective of this study is to investigate the nucleotide sequence variation of the AVR-Pita1 gene that influences the adaption of rice blast fungus to overcome the resistant gene, Pi-ta. Thirty rice blast fungus isolates were collected in 2005 and 2010 from infected rice plants in northern and northeastern Thailand. The nucleotide sequences of AVR-Pita1 were amplified and analyzed. Phylogenetic analysis was conducted using the MEGA 5.0 program. The results showed a high level of nucleotide sequence polymorphisms and the positive genetic selection pressure in Thai rice blast isolates. The details of sequence variation analysis were described in this article. The information from this study can be used for rice blast resistant breeding program in the future.     

Genome‐wide association study identifies an NLR gene that confers partial resistance to Magnaporthe oryzae in rice

Because of the frequent breakdown of major resistance (R) genes, identification of new partial R genes against rice blast disease is an important goal of rice breeding. In this study, we used a core collection of the Rice Diversity Panel II (C‐RDP‐II), which contains 584 rice accessions and are genotyped with 700 000 single‐nucleotide polymorphism (SNP) markers. The C‐RDP‐II accessions were inoculated with three blast strains collected from different rice‐growing regions in China. Genome‐wide association study identified 27 loci associated with rice blast resistance (LABRs). Among them, 22 LABRs were not associated with any known blast R genes or QTLs. Interestingly, a nucleotide‐binding site leucine‐rich repeat (NLR) gene cluster exists in the LABR12 region on chromosome 4. One of the NLR genes is highly conserved in multiple partially resistant rice cultivars, and its expression is significantly up‐regulated at the early stages of rice blast infection. Knockout of this gene via CRISPR‐Cas9 in transgenic plants partially reduced blast resistance to four blast strains. The identification of this new non‐strain specific partial R gene, tentatively named rice blast Partial Resistance gene 1 (PiPR1), provides genetic material that will be useful for understanding the partial resistance mechanism and for breeding durably resistant cultivars against blast disease of rice.     

SSRs for Marker-Assisted Selection for Blast Resistance in Rice (Oryza sativa L.)

Rice blast caused by the fungus Magnaporthe oryzae is one of the most devastating diseases of rice in nearly all rice growing areas of the world including Malaysia. To develop cultivars with resistance against different races of M. oryzae, availability of molecular markers along with marker-assisted selection strategies are essential. In this study, 11 polymorphic simple sequence repeat (SSR) markers with good fit of 1:2:1 ratio for single gene model in F₂ population derived from the cross of Pongsu seribu 2 (Resistant) and Mahsuri (Susceptible) rice cultivars were analysed in 296 F₃ families derived from individual F₂ plants to investigate association with Pi gene conferring resistance to M. oryzae pathotype. Parents and progeny were grouped into two phenotypic classes based on their blast reactions. Chi-square test for the segregation of resistance and susceptibility in F₃ generation fitted a ratio of approximately 3:1. Association of SSR markers with phenotypic trait in F₃ families was identified by statistical analysis. Four SSR markers (RM413, RM5961, RM1233 and RM8225) were significantly associated with blast resistance to pathotype 7.2 of M. oryzae in rice (p ≤ 0.01). These four markers accounted for about 20% of total phenotypic variation. So, these markers were confirmed as suitable markers for use in marker-assisted selection and confirmation of blast resistance genes to develop rice cultivars with durable blast resistance in Malaysian rice breeding programmes.     

Regions outside the leucine-rich repeat domain determine the distinct resistance specificities of the rice blast resistance genes Pik and Pik-m

Two alleles of the rice blast resistance (R) Pik locus, Pik-m and Pik, are each composed of a pair of nucleotide-binding site–leucine-rich repeat (NBS–LRR) genes, referred to as the first gene and the second gene. Pik-m and Pik are unique in that many of the amino acid substitutions between them are distributed in or near the N-terminal coiled-coil (CC) domain of the first gene, suggesting that the CC domain of the first gene plays an important role in determinating their R specificity. To examine this hypothesis, I investigated resistance phenotypes of transgenic plants carrying each of two kinds of domain-swapped Pik-m-based recombinant first genes. Replacement of the LRR domain of Pik-m with the equivalent region of Pik did not change the Pik-m-type specificity, indicating that regions outside the LRR domain are responsible for differentiating the R specificity of Pik-m from Pik. In contrast, replacement of both the NBS and LRR domains of Pik-m with the corresponding region of Pik resulted in loss of blast resistance, suggesting that co-adaptation of polymorphisms in the CC and NBS domains is necessary to maintain resistance.     

Gene interactions and genetics of blast resistance and yield attributes in rice (Oryza sativa L.)

Blast disease caused by the pathogen Pyricularia oryzae is a serious threat to rice production. Six generations viz., P₁, P₂, F₁, F₂, B₁ and B₂ of a cross between blast susceptible high-yielding rice cultivar ADT 43 and resistant near isogenic line (NIL) CT13432-3R, carrying four blast resistance genes Pi1, Pi2, Pi33 and Pi54 in combination were used to study the nature and magnitude of gene action for disease resistance and yield attributes. The epistatic interaction model was found adequate to explain the gene action in most of the traits. The interaction was complementary for number of productive tillers, economic yield, lesion number, infected leaf area and potential disease incidence but duplicate epistasis was observed for the remaining traits. Among the genotypes tested under epiphytotic conditions, gene pyramided lines were highly resistant to blast compared to individuals with single genes indicating that the nonallelic genes have a complementary effect when present together. The information on genetics of various contributing traits of resistance will further aid plant breeders in choosing appropriate breeding strategy for blast resistance and yield enhancement in rice.     

Identification and fine mapping of a major R gene to Magnaporthe oryzae in a broad-spectrum resistant germplasm in rice

Identification of R genes and development of associated molecular markers will facilitate their application in the development of crop cultivars resistant to disease. We evaluated the resistance of a resistant germplasm ??D69??, 10 monogenic lines, and model cultivar ??Nipponbare?? to 56 M. oryzae isolates of blast disease in rice. The results demonstrated that only D69 exhibited full-spectrum resistance among the 12 investigated materials. Resistance inheritance in D69 was analyzed using a stable isolate GD08T13 with strong pathogenicity, collected from diseased panicles. A single dominant R gene was revealed and designated as Pi51(t). Through linkage analysis and the development of new markers, Pi51(t) was subsequently delimited to an interval of ~100.8?kb flanked by markers Ind306 and RM19818, where Pi2, Pi9, Piz, Piz-t, Pigm(t), and Pi40(t) reside. Different genotypes identified by linked markers pB8, Pi9-2, zt56591, and T845, and different pathotypes to the same set of isolates, distinguished Pi51(t) from Pi2, Pi9, Piz, and Piz-t. The origin of Pi40(t) in wild rice suggests that Pi51(t) and Pi40(t) are different. Comparison of resistance spectra suggests multiple R genes in D69, making its resistance durable and valuable in breeding programs. The results of this work will facilitate future studies on cloning and functional analysis of blast resistance genes for rice improvement.     

Genomic structure and evolution of the Pi2/9 locus in wild rice species

Rice blast, caused by the fungal pathogen Magnaporthe oryzae, is a devastating disease of rice worldwide. Among the 85 mapped resistance (R) genes against blast, 13 have been cloned and characterized. However, how these genes originated and how they evolved in the Oryza genus remains unclear. We previously cloned the rice blast R-genes Pi2, Pi9, and Piz-t, and analyzed their genomic structure and evolution in cultivated rice. In this study, we determined the genomic sequences of the Pi2/9 locus in four wild Oryza species representing three genomes (AA, BB and CC). The number of Pi2/9 family members in the four wild species ranges from two copies to 12 copies. Although these genes are conserved in structure and categorized into the same subfamily, sequence duplications and subsequent inversions or uneven crossing overs were observed, suggesting that the locus in different wild species has undergone dynamic changes. Positive selection was found in the leucine-rich repeat region of most members, especially in the largest clade where Pi9 is included. We also provide evidence that the Pi9 gene is more related to its homologues in the recurrent line and other rice cultivars than to those in its alleged donor species O. minuta, indicating a possible origin of the Pi9 gene from O. sativa. Comparative sequence analysis between the four wild Oryza species and the previously established reference sequences in cultivated rice species at the Pi2/9 locus has provided extensive and unique information on the genomic structure and evolution of a complex R-gene cluster in the Oryza genus.     

Gene diagnosis and targeted breeding for blast-resistant Kongyu 131 without changing regional adaptability

The fungus Magnaporthe oryzae threatens the rice production of Kongyu 131 (KY131), a leading japonica variety in Northeast China. In this study, two rice lines, KP1 and KP2-Hd1, were obtained by introgressing the blast resistance genes Pi1 and Pi2 into KY131, respectively. However, both lines headed later than KY131. RICE60K SNP array analysis showed that Hd1 closely linked to Pi2 was introgressed into KP2-Hd1, and the linkage drag of Hd1 was broken by recombination. On the other hand, no known flowering genes were introgressed into KP1. Gene diagnosis by resequencing six flowering genes showed that KP1 carried functional Hd16 and Ghd8 alleles. Due to its suppression role in heading under long-day conditions, Ghd8 was chosen as the target for gene editing to disrupt its function. Four sgRNAs targeting different sites within Ghd8 were utilized to induce large-deletion mutations, which were easy to detect via agarose gel electrophoresis. All the ghd8-mutated KP1 lines were resistant to rice blast disease and headed earlier than the control KP1, even than KY131, under natural long-day conditions, which ensures its growth in Northeast China. This study confirmed that a combination of gene diagnosis and targeted gene editing is a highly efficient way to quickly eliminate undesired traits in a breeding line.     

RFLP- and physical mapping of resistance gene homologues in rice (O. sativa) and Barley (H. vulgare)

 The deduced peptide sequences of 25 gene fragments of NBS-LRR resistance (R) gene homologues from rice and barley and of characterized R genes were compared, revealing a string of six conserved motifs. Mapping of the R-gene candidates in rice showed linkage to genes conferring race-specific resistance to rice blast (Pi-k, Pi-f and Pi-1) and bacterial blight disease (Xa-1, Xa-3 and Xa-4), in barley to powdery mildew (Mla) and the rust fungus (Rpg1). In rice four mixed clusters were detected, each harboring at least two highly dissimilar NBS-LRR genes. A YAC-contig was established for one of these mixed clusters. YAC fragmentation experiments revealed the presence of at least five NBS-LRR genes within 200 kb in head-to-tail orientation. Received: 24 July 1998 / Accepted: 14 August 1998