Comparative analysis of NBS domain sequences of NBS-LRR disease resistance genes from sunflower, lettuce, and chicory

Plant resistance to many types of pathogens and pests can be achieved by the presence of disease resistance (R) genes. The nucleotide binding site-leucine rich repeat (NBS-LRR) class of R-genes is the most commonly isolated class of R-genes and makes up a super-family, which is often arranged in the genome as large multi-gene clusters. The NBS domain of these genes can be targeted by polymerase chain reaction (PCR) amplification using degenerate primers. Previous studies have used PCR derived NBS sequences to investigate both ancient R-gene evolution and recent evolution within specific plant families. However, comparative studies with the Asteraceae family have largely been ignored. In this study, we address recent evolution of NBS sequences within the Asteraceae and extend the comparison to the Arabidopsis thaliana genome. Using multiple sets of primers, NBS fragments were amplified from genomic DNA of three species from the family Asteraceae: Helianthus annuus (sunflower), Lactuca sativa (lettuce), and Cichorium intybus (chicory). Analysis suggests that Asteraceae species share distinct families of R-genes, composed of genes related to both coiled-coil (CC) and toll-interleukin-receptor homology (TIR) domain containing NBS-LRR R-genes. Between the most closely related species, (lettuce and chicory) a striking similarity of CC subfamily composition was identified, while sunflower showed less similarity in structure. These sequences were also compared to the A. thaliana genome. Asteraceae NBS gene subfamilies appear to be distinct from Arabidopsis gene clades. These data suggest that NBS families in the Asteraceae family are ancient, but also that gene duplication and gene loss events occur and change the composition of these gene subfamilies over time.     

The genetic architecture of resistance

Plant resistance genes (R genes), especially the nucleotide binding site leucine-rich repeat (NBS-LRR) family of sequences, have been extensively studied in terms of structural organization, sequence evolution and genome distribution. These studies indicate that NBS-LRR sequences can be split into two related groups that have distinct amino-acid motif organizations, evolutionary histories and signal transduction pathways. One NBS-LRR group, characterized by the presence of a Toll/interleukin receptor domain at the amino-terminal end, seems to be absent from the Poaceae. Phylogenetic analysis suggests that a small number of NBS-LRR sequences existed among ancient Angiosperms and that these ancestral sequences diversified after the separation into distinct taxonomic families. There are probably hundreds, perhaps thousands, of NBS-LRR sequences and other types of R gene-like sequences within a typical plant genome. These sequences frequently reside in 'mega-clusters' consisting of smaller clusters with several members each, all localized within a few million base pairs of one another. The organization of R-gene clusters highlights a tension between diversifying and conservative selection that may be relevant to gene families that are unrelated to disease resistance.     

Putting knowledge of plant disease resistance genes to work.

Plant disease resistance genes trigger defence mechanisms upon recognition of pathogen compatibility factors, which are encoded by avirulence genes. Isolation of the barley powdery mildew resistance gene Mla opens the door to understanding the extensive allelic diversity of this locus. Completion of the Arabidopsis genome sequence enables the analysis of the complete set of R-gene homologues in a flowering plant. A new R gene, RPW8, conferring resistance in Arabidopsis to powdery mildew, reveals a new class of protein associated with pathogen recognition. New prospects for using R-gene polymorphism in agriculture are becoming apparent.     

Characterization and comparative sequence analysis of the DNA mismatch repair MSH2 and MSH7 genes from tomato

DNA mismatch repair proteins play an essential role in maintaining genomic integrity during replication and genetic recombination. We successfully isolated a full length MSH2 and partial MSH7 cDNAs from tomato, based on sequence similarity between MutS and plant MSH homologues. Semi-quantitative RT-PCR reveals higher levels of mRNA expression of both genes in young leaves and floral buds. Genetic mapping placed MSH2 and MSH7 on chromosomes 6 and 7, respectively, and indicates that these genes exist as single copies in the tomato genome. Analysis of protein sequences and phylogeny of the plant MSH gene family show that these proteins are evolutionarily conserved, and follow the classical model of asymmetric protein evolution. Genetic manipulation of the expression of these MSH genes in tomato will provide a potentially useful tool for modifying genetic recombination and hybrid fertility between wide crosses.     

Structure, function and evolution of plant disease resistance genes

Gene-for-gene plant disease resistance involves two basic processes: perception of pathogen attack, followed by responses to limit disease. Perception involves receptors with high degrees of specificity for pathogen strains, which are encoded by disease resistance genes. Large repertoires of distantly related resistance (R) genes with diverse recognitional specificities are found within a single plant species. The generation of R-gene polymorphism involves gene duplication, followed by DNA-sequence divergence by point mutation, and by deletion and duplication of intragenic DNA repeats encoding blocks of leucine-rich elements. Recombination between related genes reassorts this variation to further diversify gene sequences. Pathogen pressure selects functional resistance specificities and results in the maintenance of R-gene diversity. Recent genome-sequence data reveal that the NBS-LRR (i.e. nucleotide-binding site-leucine-rich repeat) class of R genes represents as much as 1% of the Arabidopsis genome. Experimental data have shown that the LRR has a role in determination of specificity. Mutation experiments, in which R-gene signaling has been dissociated from specificity in constitutive signal mutants, have provided the potential for non-specific resistance to be expressed from pathogen-infection-induced promoters in transgenic plants.     

The R1 resistance gene cluster contains three groups of independently evolving, type I R1 homologues and shows substantial structural variation among haplotypes of Solanum demissum

Cultivated and wild potatoes contain a major disease-resistance cluster on the short arm of chromosome V, including the R1 resistance (R) gene against potato late blight. To explore the functional and evolutionary significance of clustering in the generation of novel disease-resistance genes, we constructed three approximately 1 Mb physical maps in the R1 gene region, one for each of the three genomes (haplotypes) of allohexaploid Solanum demissum, the wild potato progenitor of the R1 locus. Totals of 691, 919 and 559 kb were sequenced for each haplotype, and three distinct resistance-gene families were identified, one homologous to the potato R1 gene and two others homologous to either the Prf or the Bs4 R-gene of tomato. The regions with R1 homologues are highly divergent among the three haplotypes, in contrast to the conserved flanking non-resistance gene regions. The R1 locus shows dramatic variation in overall length and R1 homologue number among the three haplotypes. Sequence comparisons of the R1 homologues show that they form three distinct clades in a distance tree. Frequent sequence exchanges were detected among R1 homologues within each clade, but not among those in different clades. These frequent sequence exchanges homogenized the intron sequences of homologues within each clade, but did not homogenize the coding sequences. Our results suggest that the R1 homologues represent three independent groups of fast-evolving type I resistance genes, characterized by chimeric structures resulting from frequent sequence exchanges among group members. Such genes were first identified among clustered RGC2 genes in lettuce, where they were distinguished from slow-evolving type II R-genes. Our findings at the R1 locus in S. demissum may indicate that a common or similar mechanism underlies the previously reported differentiation of type I and type II R-genes and the differentiation of type I R-genes into distinct groups, identified here.     

Structure and evolution of plant disease resistance genes

This article reviews recent advances that shed light on plant disease resistance genes, beginning with a brief overview of their structure, followed by their genomic organization and evolution. Plant disease resistance genes have been exhaustively investigated in terms of their structural organization, sequence evolution and genome distribution. There are probably hundreds of NBS-LRR sequences and other types of R-gene-like sequences within a typical plant genome. Recent studies revealed positive selection and selective maintenance of variation in plant resistance and defence-related genes. Plant resistance genes are highly polymorphic and have diverse recognition specificities. R-genes occur as members of clustered gene families that have evolved through duplication and diversification. These genes appear to evolve more rapidly than other regions of the genome, and domains such as the leucine-rich repeat, are subject to adaptive selection     

Characterization and linkage mapping of R-gene analogous DNA sequences in pea (Pisum sativum L.)

Pea (Pisum sativum L.) sequences that are analogous to the conserved nucleotide binding site (NBS) domain found in a number of plant disease resistance genes (R-genes) were cloned. Using redundant oligonucleotide primers and the polymerase chain reaction (PCR), we amplified nine pea sequences and characterised their sequences. The pea R-gene analog (RGA)- deduced amino acid sequences demonstrated significant sequence similarity with known R-gene sequences lodged in public databases. The genomic locations of eight of the pea RGAs were determined by linkage mapping. The eight RGAs identified ten loci that mapped to six linkage groups. In addition, the genomic organization of the RGAs was inferred. Both single-copy and multicopy sequence families were present among the RGAs, and the multicopy families occurred most often as tightly linked clusters of related sequences. Intraspecific copy number variability was observed in three of the RGA sequence families, suggesting that these sequence families are evolving rapidly. The genomic locations of the pea RGAs were compared with the locations of known pea R-genes and sym genes involved in the pea-rhizobia symbiosis. Two pea RGAs mapped in the genomic region containing a pea R-gene, Fw, and four pea RGAs mapped in regions of the genome containing sym genes. Received: 4 August 1999 / Accepted: 11 November 1999     

Origin, diversity and evolution of NBS-type disease-resistance gene homologues in coffee trees (Coffea L.)

The majority of plant disease-resistance genes (R-genes) isolated so far encode a predicted nucleotide-binding site (NBS) domain. NBS domains related to R-genes show a highly conserved backbone of amino acid motifs, which makes it possible to isolate resistance gene analogues (RGAs) by PCR with degenerate primers. Multiple combinations of primers with low degeneracy, designed from two conserved motifs in the NBS regions of R-genes of various plants, were used on genomic DNA from coffee trees, an important perennial tropical crop. Nine distinct classes of RGAs of the NBS-like type, representing a highly diverse sample, were isolated from Coffea arabica and C. canephora species. The analysis of one coffee RGA family suggested point mutations as the primary source of diversity. With one exception, coffee RGA families appeared to be closely related in sequence to at least one cloned R-gene. In addition, deduced amino acid sequences of coffee RGAs were identified that showed strong sequence similarity to almost all known non-TIR (Toll/Interleukin 1 Receptor)-type R-genes. The high degree of similarity between particular coffee RGAs and R-genes isolated from other angiosperm species, such as Arabidopsis, tomato and rice, indicates an ancestral relationship and the existence of common ancestors. The data obtained from coffee species suggests that the evolution of NBS-encoding sequences involves the gradual accumulation of mutations and slow rates of divergence within distinct R-gene families, rather than being a rapid process. Functional inferences drawn from the suggested pattern of evolution of NBS-type R-genes is also discussed.     

Genome-wide analysis and identification of stress-responsive genes of the CCCH zinc finger family in Solanum lycopersicum

Zinc finger genes comprise a large and diverse gene family. Based on their individual finger structures and spacing, zinc finger proteins are further divided into different families according to their specific molecular functions. Genes in the CCCH family encode zinc finger proteins containing a motif with three cysteines and one histidine. They play important roles in plant growth and development, and in response to biotic and abiotic stresses. However, the limited analysis of the genome sequence has meant that there is no detailed information concerning the CCCH zinc finger family in tomato (Solanum lycopersicum). Here, we identified 80 CCCH zinc finger protein genes in the tomato genome. A complete overview of this gene family in tomato was presented, including the chromosome locations, gene duplications, phylogeny, gene structures and protein motifs. Promoter sequences and expression profiles of putative stress-responsive members were also investigated. These results revealed that, with the exception of four genes, the 80 CCCH genes are distributed over all 12 chromosomes with different densities, and include six segmental duplication events. The CCCH family in tomato could be divided into 12 groups based on their different CCCH motifs and into eight subfamilies by phylogenetic analysis. Analysis showed that almost all CCCH genes contain putative stress-responsive cis-elements in their promoter regions. Nine CCCH genes chosen for further quantitative real-time PCR analysis showed differential expression patterns in three representative tomato tissues. In addition, their expression levels indicated that these genes are mostly involved in the response to mannitol, heat, salicylic acid, ethylene or methyl jasmonate treatments. To the best of our knowledge, this is the first report of a genome-wide analysis of the tomato CCCH zinc finger family. Our data provided valuable information on tomato CCCH proteins and form a foundation for future studies of these proteins, especially for those members that may play important roles in stress responses.     

番茄PPR基因家族的鉴定与生物信息学分析

PPR(Pentatricopeptide repeats)基因家族在植物中广泛存在, 其在植物生长发育过程中至关重要。文章采用生物信息学方法, 利用Pfam已鉴定的PPR保守结构域序列检索番茄(Solanum lycopersicum L.)基因组计划注释的蛋白序列, 最终确定了番茄中可能存在的471个PPR编码基因; 根据拟南芥(Arabidopsis thaliana L.)中鉴定的各个结构域的特点对其进行了蛋白结构分析、分类和保守序列分析, 并对番茄PPR基因家族进行了系统进化树构建、染色体定位、亚细胞定位预测、表达和GO分析等。结果表明：番茄PPR基因家族分为P和PLS两个亚家族, 各占序列数目的一半, PLS亚家族又分为PLS、E、E+和DYW四类, 且在进化树中形成不同的分支; 各个结构域在植物中非常保守; PPR基因家族分布在番茄12条染色体上, 且多数无内含子结构; 大部分PPR蛋白具有线粒体或叶绿体定位序列, GO分析表明PPR蛋白参与RNA相关的生物学过程     

Genome-level evolution of resistance genes in Arabidopsis thaliana

Pathogen resistance genes represent some of the most abundant and diverse gene families found within plant genomes. However, evolutionary mechanisms generating resistance gene diversity at the genome level are not well understood. We used the complete Arabidopsis thaliana genome sequence to show that most duplication of individual NBS-LRR sequences occurs at close physical proximity to the parent sequence and generates clusters of closely related NBS-LRR sequences. Deploying the statistical strength of phylogeographic approaches and using chromosomal location as a proxy for spatial location, we show that apparent duplication of NBS-LRR genes to ectopic chromosomal locations is largely the consequence of segmental chromosome duplication and rearrangement, rather than the independent duplication of individual sequences. Although accounting for a smaller fraction of NBS-LRR gene duplications, segmental chromosome duplication and rearrangement events have a large impact on the evolution of this multigene family. Intergenic exchange is dramatically lower between NBS-LRR sequences located in different chromosome regions as compared to exchange between sequences within the same chromosome region. Consequently, once translocated to new chromosome locations, NBS-LRR gene copies have a greater likelihood of escaping intergenic exchange and adopting new functions than do gene copies located within the same chromosomal region. We propose an evolutionary model that relates processes of genome evolution to mechanisms of evolution for the large, diverse, NBS-LRR gene family.     

Arabidopsis thaliana: a source of candidate disease-resistance genes for Brassica napus.

Common structural and amino acid motifs among cloned plant disease-resistance genes (R genes), have made it possible to identify putative disease-resistance sequences based on DNA sequence identity. Mapping of such R-gene homologues will identify candidate disease-resistance loci to expedite map-based cloning strategies in complex crop genomes. Arabidopsis thaliana expressed sequence tags (ESTs) with homology to cloned plant R genes (R-ESTs), were mapped in both A. thaliana and Brassica napus to identify candidate R-gene loci and investigate intergenomic collinearity. Brassica R-gene homologous sequences were also mapped in B. napus. In total, 103 R-EST loci and 36 Brassica R-gene homologous loci were positioned on the N-fo-61-9 B. napus genetic map, and 48 R-EST loci positioned on the Columbia x Landsberg A. thaliana map. The mapped loci identified collinear regions between Arabidopsis and Brassica which had been observed in previous comparative mapping studies; the detection of syntenic genomic regions indicated that there was no apparent rapid divergence of the identified genomic regions housing the R-EST loci.     

Characterisation and genetic mapping of resistance and defence gene analogs in cocoa (Theobroma cacao L.)

Disease resistance and defence gene analog (RGA/DGA) sequences were isolated in cocoa using a PCR approach with degenerate primers designed from conserved domains of plant resistance and defence genes: the NBS (nucleotide binding site) motif present in a number of resistance genes such as the tobacco N, sub-domains of plant serine/threonine kinases such as the Pto tomato gene, and conserved domains of two defence gene families: pathogenesis-related proteins (PR) of classes 2 and 5. Nucleotide identity between thirty six sequences isolated from cocoa and known resistance or defence genes varied from 58 to 80%. Amino acid sequences translated from corresponding coding sequences produced sequences without stop codons, except for one NBS –like sequence. Most of the RGAs could be mapped on the cocoa genome and three clusters of genes could be observed : NBS-like sequences clustered in two regions located on chromosomes 7 and 10, Pto-like sequences mapped in five genome regions of which one, located on chromosome 4, corresponded to a cluster of five different sequences. PR2-like sequences mapped in two regions located on chromosome 5 and 9 respectively. An enrichment of the genetic map with microsatellite markers allowed us to identify several co-localisations of RGAs, DGAs and QTL for resistance to Phytophthora detected in several progenies, particularly on chromosome 4 where a cluster of Pto-like sequences and 4 QTL for resistance to Phytophthora were observed. Many other serious diseases affect cocoa and the candidate genes, isolated in this study, could be of broader interest in cocoa disease management.     

Analysis of sequence, map position, and gene expression reveals conserved essential genes for iron uptake in Arabidopsis and tomato

Arabidopsis (Arabidopsis thaliana) and tomato (Lycopersicon esculentum) show similar physiological responses to iron deficiency, suggesting that homologous genes are involved. Essential gene functions are generally considered to be carried out by orthologs that have remained conserved in sequence and map position in evolutionarily related species. This assumption has not yet been proven for plant genomes that underwent large genome rearrangements. We addressed this question in an attempt to deduce functional gene pairs for iron reduction, iron transport, and iron regulation between Arabidopsis and tomato. Iron uptake processes are essential for plant growth. We investigated iron uptake gene pairs from tomato and Arabidopsis, namely sequence, conserved gene content of the regions containing iron uptake homologs based on conserved orthologous set marker analysis, gene expression patterns, and, in two cases, genetic data. Compared to tomato, the Arabidopsis genome revealed more and larger gene families coding for the iron uptake functions. The number of possible homologous pairs was reduced if functional expression data were taken into account in addition to sequence and map position. We predict novel homologous as well as partially redundant functions of ferric reductase-like and iron-regulated transporter-like genes in Arabidopsis and tomato. Arabidopsis nicotianamine synthase genes encode a partially redundant family. In this study, Arabidopsis gene redundancy generally reflected the presumed genome duplication structure. In some cases, statistical analysis of conserved gene regions between tomato and Arabidopsis suggested a common evolutionary origin. Although involvement of conserved genes in iron uptake was found, these essential genes seem to be of paralogous rather than orthologous origin in tomato and Arabidopsis.     

The root-knot nematode resistance gene (Mi) in tomato: construction of a molecular linkage map and identification of dominant cDNA markers in resistant genotypes

A dominant allele at the Mi locus on chromosome 6 of tomato (Lycopersicon esculentum Mill) confers resistance to three species of root-knot nematodes (Meloidogyne). The resistance, which is associated with a localized necrotic response, was originally introduced into tomato from the wild species Lycopersicon peruvianum. As a step towards the molecular cloning of Mi, we have identified closely linked DNA markers from both cDNA and genomic DNA libraries as restriction fragment length polymorphisms (RFLPs). DNA from tomato populations segregating for nematode resistance was analyzed to generate a high-resolution genetic map of this region. Additional information on gene order was obtained by comparing the size of the introgressed L. peruvianum chromosomal segment within a collection of nematode-resistant tomato lines. Among the four cDNA markers that are tightly linked to Mi, three are dominant, i.e. L. peruvianum-specific. One cDNA marker corresponds to a gene family comprising 20-30 members, one of which is diagnostic for all nematode-resistant genotypes tested. The presence of non-homologous sequences around the Mi gene may contribute to the suppression of recombination in this region of the genome in crosses heterozygous for Mi. The potential of 'walking' from closely linked markers to Mi is discussed.     

High-Throughput Genomics Enhances Tomato Breeding Efficiency

Tomato (Solanum lycopersicum) is considered a model plant species for a group of economically important crops, such as potato, pepper, eggplant, since it exhibits a reduced genomic size (950 Mb), a short generation time, and routine transformation technologies. Moreover, it shares with the other Solanaceous plants the same haploid chromosome number and a high level of conserved genomic organization. Finally, many genomic and genetic resources are actually available for tomato, and the sequencing of its genome is in progress. These features make tomato an ideal species for theoretical studies and practical applications in the genomics field. The present review describes how structural genomics assist the selection of new varieties resistant to pathogens that cause damage to this crop. Many molecular markers highly linked to resistance genes and cloned resistance genes are available and could be used for a high-throughput screening of multiresistant varieties. Moreover, a new genomics-assisted breeding approach for improving fruit quality is presented and discussed. It relies on the identification of genetic mechanisms controlling the trait of interest through functional genomics tools. Following this approach, polymorphisms in major gene sequences responsible for variability in the expression of the trait under study are then exploited for tracking simultaneously favourable allele combinations in breeding programs using high-throughput genomic technologies. This aims at pyramiding in the genetic background of commercial cultivars alleles that increase their performances. In conclusion, tomato breeding strategies supported by advanced technologies are expected to target increased productivity and lower costs of improved genotypes even for complex traits.Key Words: Solanum lycopersicum, genetic and genomic resources, molecular markers, microarray, resistance to pathogens, fruit quality.     

The R1 gene for potato resistance to late blight (Phytophthora infestans) belongs to the leucine zipper/NBS/LRR class of plant resistance genes

Late blight caused by the oomycete Phytophthora infestans is the most destructive disease in potato cultivation worldwide. New, more virulent P. infestans strains have evolved which overcome the genetic resistance that has been introgressed by conventional breeding from wild potato species into commercial varieties. R genes (for single-gene resistance) and genes for quantitative resistance to late blight are present in the germplasm of wild and cultivated potato. The molecular basis of single-gene and quantitative resistance to late blight is unknown. We have cloned R1, the first gene for resistance to late blight, by combining positional cloning with a candidate gene approach. The R1 gene is member of a gene family. It encodes a protein of 1293 amino acids with a molecular mass of 149.4 kDa. The R1 gene belongs to the class of plant genes for pathogen resistance that have a leucine zipper motif, a putative nucleotide binding domain and a leucine-rich repeat domain. The most closely related plant resistance gene (36% identity) is the Prf gene for resistance to Pseudomonas syringae of tomato. R1 is located within a hot spot for pathogen resistance on potato chromosome V. In comparison to the susceptibility allele, the resistance allele at the R1 locus represents a large insertion of a functional R gene.