Characterization of the rice blast resistance gene Pik cloned from Kanto51

To study similar, but distinct, plant disease resistance (R) specificities exhibited by allelic genes at the rice blast resistance locus Pik/Pikm, we cloned the Pik gene from rice cultivar Kanto51 and compared its molecular features with those of Pikm and of another Pik gene cloned from cv. Kusabue. Like Pikm, Pik is composed of two adjacent NBS-LRR (nucleotide-binding site, leucine-rich repeat) genes: the first gene, Pik1-KA, and the second gene, Pik2-KA. Pik from Kanto51 and Pik from Kusabue were not identical; although the predicted protein sequences of the second genes were identical, the sequences differed by three amino acids within the NBS domain of the first genes. The Pik proteins from Kanto51 and Kusabue differed from Pikm in eight and seven amino acids, respectively. Most of these substituted amino acids were within the coiled-coil (CC) and NBS domains encoded by the first gene. Of these substitutions, all within the CC domain were conserved between the two Pik proteins, whereas all within the NBS domain differed between them. Comparison of the two Pik proteins and Pikm suggests the importance of the CC domain in determining the resistance specificities of Pik and Pikm. This feature contrasts with that of most allelic or homologous NBS-LRR genes characterized to date, in which the major specificity determinant is believed to lie in the highly diverged LRR domain. In addition, our study revealed high evolutionary flexibility in the genome at the Pik locus, which may be relevant to the generation of new R specificities at this locus.     

Genetic diversity of NBS–LRR class disease-resistance gene analogs in cultivated and wild eggplants

The genes encoding the nucleotide-binding site (NBS) and leucine-rich repeat (LRR) motifs constitute a large gene family in plants and have attracted much interest, because most of the plant disease-resistance genes that have been cloned are from this gene family. In this study, degenerate oligonucleotide primers, designed on the basis of conserved regions of the NBS domains from known plant resistance genes, were used to isolate resistance gene analogs (RGAs) from cultivated and wild eggplants, i.e., S. melongena, S. aethiopicum gr. Gilo, S. linnaeanum, S. integrifolium, S. sisymbriifolium, and S. khasianum. Sequence analysis indicated that the cloned eggplant RGAs belong to the non-TIR–NBS–LRR type, which are very similar to the R genes or the RGAs identified in other plant species, especially Solanaceae plants, suggesting the existence of common ancestors. Wide genetic diversity of eggplant RGAs was observed both in interspecific and intraspecific sequences, and eight distinct families of eggplant RGAs were identified. Further studies revealed a high average ratio of synonymous to non-synonymous substitution and a low level of recombination. These results suggest that NBS-encoding sequences of RGAs in cultivated and wild eggplants are subject to gradual accumulation of mutations leading to purifying selection. This is the first report of NBS–LRR class RGAs in eggplants.     

Two different CC‐NBS‐LRR genes are required for Lr10‐mediated leaf rust resistance in tetraploid and hexaploid wheat

Comparative study of disease resistance genes in crop plants and their relatives provides insight on resistance gene function, evolution and diversity. Here, we studied the allelic diversity of the Lr10 leaf rust resistance gene, a CC‐NBS‐LRR coding gene originally isolated from hexaploid wheat, in 20 diploid and tetraploid wheat lines. Besides a gene in the tetraploid wheat variety ‘Altar’ that is identical to the hexaploid wheat Lr10, two additional, functional resistance alleles showing sequence diversity were identified by virus‐induced gene silencing in tetraploid wheat lines. In contrast to most described NBS‐LRR proteins, the N‐terminal CC domain of LR10 was found to be under strong diversifying selection. A second NBS‐LRR gene at the Lr10 locus, RGA2, was shown through silencing to be essential for Lr10 function. Interestingly, RGA2 showed much less sequence diversity than Lr10. These data demonstrate allelic diversity of functional genes at the Lr10 locus in tetraploid wheat, and these new genes can now be analyzed for agronomic relevance. Lr10‐based resistance is highly unusual both in its dependence on two, only distantly, related CC‐NBS‐LRR proteins, as well as in the pattern of diversifying selection in the N‐terminal domain. This indicates a new and complex molecular mechanism of pathogen detection and signal transduction.     

Isolation of a full-length CC–NBS–LRR resistance gene analog candidate from sugar pine showing low nucleotide diversity

The nucleotide-binding-site and leucine-rich-repeat (NBS–LRR) class of R proteins is abundant and widely distributed in plants. By using degenerate primers designed on the NBS domain in lettuce, we amplified sequences in sugar pine that shared sequence identity with many of the NBS–LRR class resistance genes catalogued in GenBank. The polymerase chain reaction products were used to probe a cDNA library constructed from needle tissue of sugar pine seedlings. A full-length cDNA was obtained that demonstrated high predicted amino acid sequence similarity to the coiled coil (CC)–NBS–LRR subclass of NBS–LRR resistance proteins in GenBank. Sequence analyses of this gene in megagametophytes from two sugar pine trees segregating for the hypersensitive response to white pine blister rust revealed zero nucleotide variation. Moreover, there was no variation found in 24 unrelated sugar pine trees except for three single-nucleotide polymorphisms located in the 3′ untranslated region. Compared to other genes sequenced in Pinaceae, such a low level of sequence variation in unrelated individuals is unusual. Although, numerous studies have reported that plant R genes are under diversifying selection for specificity to evolving pathogens, the resistance gene analog discussed here appears to be under intense purifying selection.An erratum to this article can be found at     

A set of near-isogenic lines of Indica-type rice variety CO 39 as differential varieties for blast resistance

Twenty-seven near-isogenic lines (NILs) with the genetic background of a blast-susceptible variety, CO 39, were developed by repeated backcrossing as a first set of a large number of differential varieties (DVs) with Indica-type genetic background. The NILs included 14 resistance genes—Pish, Pib, Piz-5, Piz-t, Pi5(t), Pik-s, Pik, Pik-h, Pik-m, Pik-p, Pi1, Pi7(t), Pita, and Pita-2—derived from 26 donor varieties. The reaction patterns of NILs against 20 standard isolates from the Philippines were similar to those of blast monogenic lines with the same resistance gene, except for those against two isolates that are avirulent to Pia in the genetic background of CO 39. A genome-wide DNA marker survey revealed that chromosome segments were introgressed in the regions where each resistance gene was previously mapped and most of the other chromosome regions in each NIL were CO 39 type. Segregation analysis of resistance and co-segregation analysis between resistance and DNA markers using F3 populations derived from the crosses between each NIL and the recurrent parent, CO 39, revealed a single-gene control of resistance and association between resistance and target introgressed segments. The morphological characters of each NIL were almost the same as those of the recurrent parent except for some lines, suggesting that these NILs can be used even under tropical conditions where Japonica-type DVs are not suitable for cropping. Thus, these NILs are useful not only as genetic tools for blast resistance study but also as sources of genes for breeding of Indica-type rice varieties.     

Evolution of the number of LRRs in plant disease resistance genes

The largest group of plant resistance (R) genes contain the regions that encode the nucleotide-binding site (NBS) and leucine-rich repeat (LRR) domains (NBS-LRR genes). To gain new resistance, amino acid substitutions and changes in number of the LRRs that recognize the presence of pathogens are considered important. In this study, we focus on the evolution of the number of LRRs and analyze the genome data of five plant species, Arabidopsis thaliana, Oryza sativa, Medicago truncatula, Lotus japonicus and Populus trichocarpa. We first categorized the NBS-LRR genes in each species into groups and subgroups based on the phylogenetic relationships of their NBS domain sequences. Then we estimated the evolutionary rate of the number of LRRs relative to the synonymous divergence in the NBS domain sequences by a maximum likelihood method assuming the single stepwise mutation model. The estimates ranged from 4.5 to 600 and differed between groups in the same species or between species. This indicated different roles played by different groups of the NBS-LRR genes within a species or the effects of various life history characteristics, such as generation time, of the species. We also tested the fit of the model to the data using the variance of number of LRRs in each subgroup. In some subgroups in some plants (16 out of 174 subgroups), the results of simulation using the estimated rates significantly deviated from the observed data. Those subgroups may have undergone different modes of selection from the other subgroups.     

Structure-function analysis of the coiled-coil and leucine-rich repeat domains of the RPS5 disease resistance protein

The Arabidopsis (Arabidopsis thaliana) RESISTANCE TO PSEUDOMONAS SYRINGAE5 (RPS5) disease resistance protein mediates recognition of the Pseudomonas syringae effector protein AvrPphB. RPS5 belongs to the coiled-coil-nucleotide-binding site-leucine-rich repeat (CC-NBS-LRR) family and is activated by AvrPphB-mediated cleavage of the protein kinase PBS1. Here, we present a structure-function analysis of the CC and LRR domains of RPS5 using transient expression assays in Nicotiana benthamiana. We found that substituting the CC domain of RPS2 for the RPS5 CC domain did not alter RPS5 specificity and only moderately reduced its ability to activate programmed cell death, suggesting that the CC domain does not play a direct role in the recognition of PBS1 cleavage. Analysis of an RPS5-super Yellow Fluorescent Protein fusion revealed that RPS5 localizes to the plasma membrane (PM). Alanine substitutions of predicted myristoylation (glycine-2) and palmitoylation (cysteine-4) residues affected RPS5 PM localization, protein stability, and function in an additive manner, indicating that PM localization is essential to RPS5 function. The first 20 amino acids of RPS5 were sufficient for directing super Yellow Fluorescent Protein to the PM. C-terminal truncations of RPS5 revealed that the first four LRR repeats are sufficient for inhibiting RPS5 autoactivation; however, the complete LRR domain was required for the recognition of PBS1 cleavage. Substitution of the RPS2 LRR domain resulted in the autoactivation of RPS5, indicating that the LRR domain must coevolve with the NBS domain. We conclude that the RPS5 LRR domain functions to suppress RPS5 activation in the absence of PBS1 cleavage and promotes RPS5 activation in its presence.     

Six amino acid changes confined to the leucine-rich repeat beta-strand/beta-turn motif determine the difference between the P and P2 rust resistance specificities in flax

At least six rust resistance specificities (P and P1 to P5) map to the complex P locus in flax. The P2 resistance gene was identified by transposon tagging and transgenic expression. P2 is a member of a small multigene family and encodes a protein with nucleotide binding site (NBS) and leucine-rich repeat (LRR) domains and an N-terminal Toll/interleukin-1 receptor (TIR) homology domain, as well as a C-terminal non-LRR (CNL) domain of approximately 150 amino acids. A related CNL domain was detected in almost half of the predicted Arabidopsis TIR-NBS-LRR sequences, including the RPS4 and RPP1 resistance proteins, and in the tobacco N protein, but not in the flax L and M proteins. Presence or absence of this domain defines two subclasses of TIR-NBS-LRR resistance genes. Truncations of the P2 CNL domain cause loss of function, and evidence for diversifying selection was detected in this domain, suggesting a possible role in specificity determination. A spontaneous rust-susceptible mutant of P2 contained a G-->E amino acid substitution in the GLPL motif, which is conserved in the NBS domains of plant resistance proteins and the animal cell death control proteins APAF-1 and CED4, providing direct evidence for the importance of this motif in resistance gene function. A P2 homologous gene isolated from a flax line expressing the P resistance specificity encodes a protein with only 10 amino acid differences from the P2 protein. Chimeric gene constructs indicate that just six of these amino acid changes, all located within the predicted beta-strand/beta-turn motif of four LRR units, are sufficient to alter P2 to the P specificity.     

Expression of two soybean resistance gene candidates shows divergence of paralogous single-copy genes

The cloning of several plant genes directly involved in triggering a disease resistance response has shown that numerous resistance genes in the nucleotide binding site (NBS)/leucine-rich repeat (LRR) class have similar conserved amino acid sequences. In this study, we used a short soybean DNA sequence, previously cloned based on its conserved NBS, as a probe to identify full-length resistance gene candidates. Two homologous, but genetically independent genes were identified. One gene maps to the soybean molecular linkage group (MLG) F and a second is coded on MLG E. The first gene contains a 3,279 nucleotide open reading frame (ORF) sequence and possesses all the functional motifs characteristic of previously cloned NBS/LRR resistance genes. The N-terminal sequence of the deduced gene product is highly characteristic of other resistance genes in the subgroup of NBS/LRR genes which show homology to the Toll/Interleukin-1 receptor genes. The C-terminal region is somewhat more divergent as seen in other cloned disease resistance genes. This region of the F-linked gene contains an LRR region that is characterized by two alternatively spliced products which produce gene products with either a four-repeat or a ten-repeat LRR. The second cloned gene that maps to soybean MLG E contains 1,565 nucleotides of ORF in the N-terminal domain. Despite strong homology, however, the 3′ region of this gene contains several in-frame stop codons and apparent frame shifts compared to the F-linked gene, suggesting that its functionality as a disease resistance gene is questionable. These two disease resistance gene candidates are shown to be closely related to one another and to the members of the NBS/LRR class of disease resistance genes. Received: 29 November 1999 / Accepted: 22 December 1999     

Genome-wide identification and mapping of NBS-encoding resistance genes in Solanum tuberosum group phureja

The majority of disease resistance (R) genes identified to date in plants encode a nucleotide-binding site (NBS) and leucine-rich repeat (LRR) domain containing protein. Additional domains such as coiled-coil (CC) and TOLL/interleukin-1 receptor (TIR) domains can also be present. In the recently sequenced Solanum tuberosum group phureja genome we used HMM models and manual curation to annotate 435 NBS-encoding R gene homologs and 142 NBS-derived genes that lack the NBS domain. Highly similar homologs for most previously documented Solanaceae R genes were identified. A surprising ～41% (179) of the 435 NBS-encoding genes are pseudogenes primarily caused by premature stop codons or frameshift mutations. Alignment of 81.80% of the 577 homologs to S. tuberosum group phureja pseudomolecules revealed non-random distribution of the R-genes; 362 of 470 genes were found in high density clusters on 11 chromosomes.     

Characterization of the Fusarium wilt resistance Fom-2 gene in melon

The melon gene Fom-2, which confers resistance to Fusarium oxysporum f.sp. melonis (Fom) races 0 and 1, has been previously characterized by map-based cloning, and it encodes a protein with a nucleotide binding site (NBS) and leucine-rich repeats (LRRs). Here, we used the primer Fom2-LRR₁₆₃₉ to clone and sequence a partial LRR region of the Fom-2 gene in 11 melon accessions resistant to Fusarium wilt from various geographic regions. Our work revealed that the structure of the partial LRR domain is highly conserved between eight of these resistant accessions and is similar to the resistant allele in the previously characterized PI-161375 line. Conversely, PI-124111 is a unique line that presents the same resistant allele that was previously described in the MR-1 line. The accession Cum-355 presents a protein that differs from that encoded by both the resistant lines PI-161375 and MR-1. This result suggests that Cum-355 has a new resistant allele of Fom-2 that determines the same specificity. Importantly, based on the sequence of the Fom-2 LRR domain, two sequence characterized amplified region (SCAR) markers, Fom2-R₄₀₈ and Fom2-S₃₄₂, were developed for Fom-2 resistant and susceptible alleles, respectively. These allele-specific PCR markers could be used as co-dominant markers when their primer pairs were combined in a multiplex PCR reaction. The specificity of these functional markers (FM) was validated on a set of 27 genotypes representing several melon types. These FM markers are expected to enhance the reliability and cost-effectiveness of marker-assisted selection for the Fom-2 gene in melon.     

Sequence variation at the rice blast resistance gene Pi-km locus: Implications for the development of allele specific markers

The recently cloned blast resistance (R) gene Pi-km protects rice crops against specific races of the fungal pathogen Magnaporthe oryzae in a gene-for-gene manner. The use of blast R genes remains the most cost-effective method for an integrated disease management strategy. To facilitate rice breeding we developed a Pi-km specific DNA marker. For this purpose, we initially explored the existing sequence diversity for alleles of the two genes responsible for the Pi-km specificity. The analysis of 15 rice cultivars revealed that the majority of nucleotide polymorphisms were associated with the Pi-km1 gene. Interestingly, the correspondent amino acid variation was localized within the predicted coiled-coil domain of the putative Pi-km1 protein. In contrast, the sequence of Pi-km2 alleles was highly conserved even within distantly related cultivars. Furthermore, disease reactions of the selected cultivars to five M. oryzae isolates, as well as their determined Pi-km1 allele, showed a good correlation with the known Pi-k genes (-k/-kh/-km/-ks/-kp) historically reported for these cultivars. Based on these findings, specific primer sets have been designed to discriminate among the various Pi-km alleles. The new markers should simplify the introgression of the valuable blast resistance associated with the complex Pi-k locus into rice cultivars.     

Functional markers based molecular characterization and cloning of resistance gene analogs encoding NBS-LRR disease resistance proteins in finger millet (<Emphasis Type="Italic">Eleusine coracana)</Emphasis>

Magnaporthe grisea, the blast fungus is one of the main pathological threats to finger millet crop worldwide. A systematic search for the blast resistance gene analogs was carried out, using functional molecular markers. Three-fourths of the recognition-dependent disease resistance genes (R-genes) identified in plants encodes nucleotide binding site (NBS) leucine-rich repeat (LRR) proteins. NBS-LRR homologs have only been isolated on a limited scale from Eleusine coracana. Genomic DNA sequences sharing homology with NBS region of resistance gene analogs were isolated and characterized from resistant genotypes of finger millet using PCR based approach with primers designed from conserved regions of NBS domain. Attempts were made to identify molecular markers linked to the resistance gene and to differentiate the resistant bulk from the susceptible bulk. A total of 9 NBS-LRR and 11 EST-SSR markers generated 75.6 and 73.5% polymorphism respectively amongst 73 finger millet genotypes. NBS-5, NBS-9, NBS-3 and EST-SSR-04 markers showed a clear polymorphism which differentiated resistant genotypes from susceptible genotypes. By comparing the banding pattern of different resistant and susceptible genotypes, five DNA amplifications of NBS and EST-SSR primers (NBS-05_504, NBS-09₇₁₁, NBS-07₆₈₈, NBS-03₅₀₉ and EST-SSR-04₂₄₁) were identified as markers for the blast resistance in resistant genotypes. Principal coordinate plot and UPGMA analysis formed similar groups of the genotypes and placed most of the resistant genotypes together showing a high level of genetic relatedness and the susceptible genotypes were placed in different groups on the basis of differential disease score. Our results provided a clue for the cloning of finger millet blast resistance gene analogs which not only facilitate the process of plant breeding but also molecular characterization of blast resistance gene analogs from Eleusine coracana.     

Genome-Wide Analysis of Nucleotide-Binding Site (NBS) Disease Resistance (R) Genes in Sacred Lotus (Nelumbo nucifera Gaertn.) Reveals Their Transition Role During Early Evolution of Land Plants

Nucleotide-binding site (NBS) containing genes comprise the largest class in identified plant resistance genes. A total of 137 NBS class resistance genes were identified from the newly sequenced sacred lotus genome (Nelumbo nucifera Gaertn.) through a reiterative computational sequence analysis. Three distinct groups of NBS-encoding genes were classified: 5 with Toll/interleukin-1 receptor homology (TIR) domain at N-terminal (TIR-NBS [-LRR (leucine-rich repeat)]), 37 with CC (coiled coil) domain (CC-NBS [-LRR]), and 95 with neither TIR nor CC at N-terminal (NBS [-LRR]). Sequence analysis revealed high divergence of NBS-LRR genes in sacred lotus. The result of cluster and syntenic analysis of NBS genes suggested a duplication and recombination event, which is consistent with the correspondent result of whole genome analysis. In addition, we also identified 52 NBS genes which have a putative NACHT domain embedded in the NBS domains. This characteristic has only been reported in animals, fungi and bacteria, suggesting that NACHT and NBS domains shared a similar ancient origin; and sacred lotus NBS (NACHT) genes may represent a transition role during the early evolution of disease resistance in land plants.     

Isolation of TIR and nonTIR NBS–LRR resistance gene analogues and identification of molecular markers linked to a powdery mildew resistance locus in chestnut rose (Rosa roxburghii Tratt)

Toll and interleukin-1 receptor (TIR) and nonTIR nucleotide binding site–leucine rich repeat (NBS–LRR) resistance gene analogues (RGAs) were obtained from chestnut rose (Rosa roxburghii Tratt) by two PCR-based amplification strategies (direct amplification and overlap extension amplification) with degenerate primers designed to the conserved P-loop, kinase-2, and Gly-Leu-Pro-Leu (GLPL) motifs within the NBS domain of plant resistance gene (R gene) products. Thirty-four of 65 cloned PCR fragments contained a continuous open reading frame (ORF) and their predicted protein products showed homology to the NBS–LRR class R proteins in the GenBank database. These 34 predicted protein sequences exhibited a wide range (19.5–99.4%) of sequence identity among them and were classified into two distinct groups by phylogenetic analysis. The first group consisted of 23 sequences and seemed to belong to the nonTIR NBS–LRR RGAs, since they contained group specific motifs (RNBS-A-nonTIR motif) that are often present in the coiled-coil domain of the nonTIR NBS–LRR class R genes. The second group comprised 11 sequences that contained motifs found in the TIR domain of TIR NBS–LRR class R genes. Restriction fragment length polymorphic (RFLP) markers were developed from some of the RGAs and used for mapping powdery mildew resistance genes in chestnut rose. Three markers, RGA22C, RGA4A, and RGA7B, were identified to be linked to a resistance gene locus, designated CRPM1 for chestnut rose powdery mildew resistance 1, which accounted for 72% of the variation in powdery mildew resistance phenotype in an F1 segregating population. To our knowledge, this is the first report on isolation, phylogenetic analysis and potential utilization as genetic markers of RGAs in chestnut rose.     

Arms race co‐evolution of Magnaporthe oryzae AVR‐Pik and rice Pik genes driven by their physical interactions

Attack and counter‐attack impose strong reciprocal selection on pathogens and hosts, leading to development of arms race evolutionary dynamics. Here we show that Magnaporthe oryzae avirulence gene AVR‐Pik and the cognate rice resistance (R) gene Pik are highly variable, with multiple alleles in which DNA replacements cause amino acid changes. There is tight recognition specificity of the AVR‐Pik alleles by the various Pik alleles. We found that AVR‐Pik physically binds the N‐terminal coiled‐coil domain of Pik in a yeast two‐hybrid assay as well as in an in planta co‐immunoprecipitation assay. This binding specificity correlates with the recognition specificity between AVR and R genes. We propose that AVR‐Pik and Pik are locked into arms race co‐evolution driven by their direct physical interactions.     

The <Emphasis Type="Italic">Pik</Emphasis>-<Emphasis Type="Italic">p</Emphasis> resistance to <Emphasis Type="Italic">Magnaporthe oryzae</Emphasis> in rice is mediated by a pair of closely linked CC-NBS-LRR genes

The blast resistance gene Pik-p, mapping to the Pik locus on the long arm of rice chromosome 11, was isolated by map-based in silico cloning. Four NBS-LRR genes are present in the target region of cv. Nipponbare, and a presence/absence analysis in the Pik-p carrier cv. K60 excluded two of these as candidates for Pik-p. The other two candidates (KP3 and KP4) were expressed in cv. K60. A loss-of-function experiment by RNAi showed that both KP3 and KP4 are required for Pik-p function, while a gain-of-function experiment by complementation test revealed that neither KP3 nor KP4 on their own can impart resistance, but that resistance was expressed when both were introduced simultaneously. Both Pikp-1 (KP3) and Pikp-2 (KP4) encode coiled-coil NBS-LRR proteins and share, respectively, 95 and 99% peptide identity with the two alleles, Pikm1-TS and Pikm2-TS. The Pikp-1 and Pikp-2 sequences share only limited homology. Their sequence allowed Pik-p to be distinguished from Pik, Pik-s, Pik-m and Pik-h. Both Pikp-1 and Pikp-2 were constitutively expressed in cv. K60 and only marginally induced by blast infection.