Evolutionary meta-analysis of solanaceous resistance gene and solanum resistance gene analog sequences and a practical framework for cross-species comparisons

Cross-species comparative genomics approaches have been employed to map and clone many important disease resistance (R) genes from Solanum species-especially wild relatives of potato and tomato. These efforts will increase with the recent release of potato genome sequence and the impending release of tomato genome sequence. Most R genes belong to the prominent nucleotide binding site-leucine rich repeat (NBS-LRR) class and conserved NBS-LRR protein motifs enable survey of the R gene space of a plant genome by generation of resistance gene analogs (RGA), polymerase chain reaction fragments derived from R genes. We generated a collection of 97 RGA from the disease-resistant wild potato S. bulbocastanum, complementing smaller collections from other Solanum species. To further comparative genomics approaches, we combined all known Solanum RGA and cloned solanaceous NBS-LRR gene sequences, nearly 800 sequences in total, into a single meta-analysis. We defined R gene diversity bins that reflect both evolutionary relationships and DNA cross-hybridization results. The resulting framework is amendable and expandable, providing the research community with a common vocabulary for present and future study of R gene lineages. Through a series of sequence and hybridization experiments, we demonstrate that all tested R gene lineages are of ancient origin, are shared between Solanum species, and can be successfully accessed via comparative genomics approaches.     

Co-Variation Among Major Classes of LRR-Encoding Genes in Two Pairs of Plant Species

NBS-LRR (nucleotide-binding site-leucine-rich repeat), LRR-RLK (LRR-receptor-like kinase), and LRR-only are the three major LRR-encoding genes. Owing to the crucial role played by them in plant resistance, development, and growth, extensive studies have been performed on the NBS-LRR and LRR-RLK genes. However, few studies have focused on these genes collectively; they may co-vary as all of them contain LRR motifs. To investigate their common evolutionary patterns, all major classes of LRR-encoding genes were identified in 12 plant species, and particularly compared in two pairs of close relatives, Arabidopsis thaliana-A. lyrata (At-Al) and Zea mays-Sorghum bicolor. Our results showed that these genes co-vary significantly in terms of their numbers between species and that the genes with certain evolutionary parameters are most likely to have similar functions. The development-related genes have clear orthologous relationships between closely related species, as well as lower nucleotide divergence, and Ka/Ks ratio. In contrast, resistance-related genes have exactly opposite characteristics and favor 11-15 LRRs per gene. This association could be very useful in predicting the function of LRR-encoding genes. The presence of co-variation suggests that LRRs, combined with other domains, can work better in some common functions. In order to cooperate efficiently, there should be balanced gene numbers among the different gene classes.     

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.     

Heterogeneous evolutionary rates of Pi2/9 homologs in rice

ABSTRACT: BACKGROUND: The Pi2/9 locus contains multiple nucleotide binding site--leucine-rich repeat (NBS-LRR) genes in the rice genome. Although three functional R-genes have been cloned from this locus, little is known about the origin and evolutionary history of these genes. Herein, an extensive genome-wide survey of Pi2/9 homologs in rice, sorghum, Brachypodium and Arabidopsis, was conducted to explore this theme. RESULTS: In our study, 1, 1, 5 and 156 Pi2/9 homologs were detected in Arabidopsis, Brachypodium, sorghum and rice genomes, respectively. Two distinct evolutionary patterns of Pi2/9 homologs, Type I and Type II, were observed in rice lines. Type I Pi2/9 homologs showed evidence of rapid gene diversification, including substantial copy number variations, obscured orthologous relationships, high levels of nucleotide diversity or/and divergence, frequent sequence exchanges and strong positive selection, whereas Type II Pi2/9 homologs exhibited a fairly slow evolutionary rate. Interestingly, the three cloned R-genes from the Pi2/9 locus all belonged to the Type I genes. CONCLUSIONS: Our data show that the Pi2/9 locus had an ancient origin predating the common ancestor of gramineous species. The existence of two types of Pi2/9 homologs suggest that diversifying evolution should be an important strategy of rice to cope with different types of pathogens. The relationship of cloned Pi2/9 genes and Type I genes also suggests that rapid gene diversification might facilitate rice to adapt quickly to the changing spectrum of the fungal pathogen M. grisea. Based on these criteria, other potential candidate genes that might confer novel resistance specificities to rice blast could be predicted.     

Uncovering the dynamic evolution of nucleotide‐binding site‐leucine‐rich repeat (NBS‐LRR) genes in Brassicaceae

Plant genomes harbor dozens to hundreds of nucleotide-binding site-leucine-rich repeat(NBS-LRR) genes;however,the long-term evolutionary history of these resistance genes has not been fully understood. This study focuses on five Brassicaceae genomes and the Carica papaya genome to explore changes in NBS-LRR genes that have taken place in this Rosid II lineage during the past 72 million years. Various numbers of NBS-LRR genes were identified from Arabidopsis lyrata(198),A. thaliana(165),Brassica rapa(204),Capsella rubella(127),Thellungiella salsuginea(88),and C. papaya(51). In each genome,the identified NBS-LRR genes were found to be unevenly distributed among chromosomes and most of them were clustered together.Phylogenetic analysis revealed that,before and after Brassicaceae speciation events,both toll/interleukin-1receptor-NBS-LRR(TNL) genes and non-toll/interleukin-1receptor-NBS-LRR(n TNL) genes exhibited a pattern of first expansion and then contraction,suggesting that both subclasses of NBS-LRR genes were responding to pathogen pressures synchronically. Further,by examining the gain/loss of TNL and n TNL genes at different evolutionary nodes,this study revealed that both events often occurred more drastically in TNL genes. Finally,the phylogeny of n TNL genes suggested that this NBS-LRR subclass is composed o two separate ancient gene types: RPW8-NBS-LRR and Coiled-coil-NBS-LRR.     

Genome duplication and the origin of angiosperms

Despite intensive research, little is known about the origin of the angiosperms and their rise to ecological dominance during the Early Cretaceous. Based on whole-genome analyses of Arabidopsis thaliana, there is compelling evidence that angiosperms underwent two whole-genome duplication events early during their evolutionary history. Recent studies have shown that these events were crucial for the creation of many important developmental and regulatory genes found in extant angiosperm genomes. Here, we argue that these ancient polyploidy events might have also had an important role in the origin and diversification of the angiosperms.     

Tracing the origin and evolution of plant TIR-encoding genes

Toll-interleukin-1 receptor (TIR)-encoding proteins represent one of the most important families of disease resistance genes in plants. Studies that have explored the functional details of these genes tended to focus on only a few limited groups; the origin and evolutionary history of these genes were therefore unclear. In this study, focusing on the four principal groups of TIR-encoding genes, we conducted an extensive genome-wide survey of 32 fully sequenced plant genomes and Expressed Sequence Tags (ESTs) from the gymnosperm Pinus taeda and explored the origins and evolution of these genes. Through the identification of the TIR-encoding genes, the analysis of chromosome positions, the identification and analysis of conserved motifs, and sequence alignment and phylogenetic reconstruction, our results showed that the genes of the TIR-X family (TXs) had an earlier origin and a wider distribution than the genes from the other three groups. TIR-encoding genes experienced large-scale gene duplications during evolution. A skeleton motif pattern of the TIR domain was present in all spermatophytes, and the genes with this skeleton pattern exhibited a conserved and independent evolutionary history in all spermatophytes, including monocots, that followed their gymnosperm origin. This study used comparative genomics to explore the origin and evolutionary history of the four main groups of TIR-encoding genes. Additionally, we unraveled the mechanism behind the uneven distribution of TIR-encoding genes in dicots and monocots.     

Red and green algal origin of diatom membrane transporters: insights into environmental adaptation and cell evolution

Membrane transporters (MTs) facilitate the movement of molecules between cellular compartments. The evolutionary history of these key components of eukaryote genomes remains unclear. Many photosynthetic microbial eukaryotes (e.g., diatoms, haptophytes, and dinoflagellates) appear to have undergone serial endosymbiosis and thereby recruited foreign genes through endosymbiotic/horizontal gene transfer (E/HGT). Here we used the diatoms Thalassiosira pseudonana and Phaeodactylum tricornutum as models to examine the evolutionary origin of MTs in this important group of marine primary producers. Using phylogenomics, we used 1,014 diatom MTs as query against a broadly sampled protein sequence database that includes novel genome data from the mesophilic red algae Porphyridium cruentum and Calliarthron tuberculosum, and the stramenopile Ectocarpus siliculosus. Our conservative approach resulted in 879 maximum likelihood trees of which 399 genes show a non-lineal history between diatoms and other eukaryotes and prokaryotes (at the bootstrap value ≥70%). Of the eukaryote-derived MTs, 172 (ca. 25% of 697 examined phylogenies) have members of both red/green algae as sister groups, with 103 putatively arising from green algae, 19 from red algae, and 50 have an unresolved affiliation to red and/or green algae. We used topology tests to analyze the most convincing cases of non-lineal gene history in which red and/or green algae were nested within stramenopiles. This analysis showed that ca. 6% of all trees (our most conservative estimate) support an algal origin of MTs in stramenopiles with the majority derived from green algae. Our findings demonstrate the complex evolutionary history of photosynthetic eukaryotes and indicate a reticulate origin of MT genes in diatoms. We postulate that the algal-derived MTs acquired via E/HGT provided diatoms and other related microbial eukaryotes the ability to persist under conditions of fluctuating ocean chemistry, likely contributing to their great success in marine environments.     

The origin of nucleus: rebuild from the prokaryotic ancestors of ribosome export factors

Various hypotheses have been proposed on the evolutionary origin of eukaryotic nucleus. Because one of the major cargoes in the nucleocytoplasmic export in the eukaryotic cell is the ribosome, its stimulating proteins called Ribosome Export Factors (REFs) might have an evolutionary history of inscribing the origin of eukaryotic nucleus. With the aim of understanding the evolutionary origin of the nucleus, here we employed the yeast REFs and searched for their evolutionary origin in more than 500 genomes of archaea and eubacteria by the PSI-BLAST search. Our results showed that the non-membranous REFs (non-mREFs) originated exclusively from eubacterial proteins, whereas the membranous REFs (mREFs) are from both archaeal and eubacterial proteins. Since the non-mREFs just work inside the nucleus while the mREFs shuttle between the nucleus and the cytoplasm, these results suggest that the extant REFs working inside the nucleus have derived exclusively from eubacterial proteins, implying that the nucleus arose in a cell that contained chromosomes possessing a substantial fraction of eubacterial genes, in line with the predictions of several models entailing endosymbiosis at eukaryote origins.     

Genome-wide identification of NBS genes in<Emphasis Type="Italic"> japonica</Emphasis> rice reveals significant expansion of divergent non-TIR NBS-LRR genes

A complete set of candidate disease resistance ( R) genes encoding nucleotide-binding sites (NBSs) was identified in the genome sequence of japonica rice ( Oryza sativa L. var. Nipponbare). These putative R genes were characterized with respect to structural diversity, phylogenetic relationships and chromosomal distribution, and compared with those in Arabidopsis thaliana. We found 535 NBS-coding sequences, including 480 non-TIR (Toll/IL-1 receptor) NBS-LRR (Leucine Rich Repeat) genes. TIR NBS-LRR genes, which are common in A. thaliana, have not been identified in the rice genome. The number of non-TIR NBS-LRR genes in rice is 8.7 times higher than that in A. thaliana, and they account for about 1% of all of predicted ORFs in the rice genome. Some 76% of the NBS genes were located in 44 gene clusters or in 57 tandem arrays, and 16 apparent gene duplications were detected in these regions. Phylogenetic analyses based both NBS and N-terminal regions classified the genes into about 200 groups, but no deep clades were detected, in contrast to the two distinct clusters found in A. thaliana. The structural and genetic diversity that exists among NBS-LRR proteins in rice is remarkable, and suggests that diversifying selection has played an important role in the evolution of R genes in this agronomically important species. (Supplemental material is available online at .)Communicated by R. HagemannThe first three authors contributed equally to this work     

The genomic dynamics and evolutionary mechanism of the Pi2/9 locus in rice

The Pi2/9 locus contains at least four resistance specificities to Magnaporthe grisea and belongs to a gene complex comprised of multiple genes that encode highly homologous nucleotide binding site (NBS) and leucine rich repeat (LRR) proteins. To investigate the genetic events involved in the evolution of the Pi2/9 locus, we analyzed the Pi2/9 locus at the inter- and intralocus levels in five rice cultivars. The NBS-LRR genes in the five cultivars belong to the same phylogenetic clade among rice NBS-LRR genes, and all have a phase-2 intron at the N-terminus. However, the paralogs within each haplotype show a significant sequence divergence and their N-terminal intron and 5' regulatory regions are very different. On the contrary, the orthologs from different haplotypes are highly similar, indicating an obvious orthologous relationship has been maintained during the evolution of the Pi2/9 locus. These results suggest that sequence diversification in the 5' regulatory regions and N-terminal introns of the paralogs may have led to suppression of meiotic recombination between the paralogs within each haplotype, facilitating the maintenance of the orthologous relationship among rice cultivars. Our observations provide valuable insight into the genomic dynamics and evolutionary mechanism of an NBS-LRR resistance-gene complex in rice.     

Mode of amplification and reorganization of resistance genes during recent Arabidopsis thaliana evolution.

The NBS-LRR (nucleotide-binding site plus leucine-rich repeat) genes represent the major class of disease resistance genes in flowering plants and comprise 166 genes in the ecotype Col-0 of Arabidopsis thaliana. NBS-LRR genes are organized in single-gene loci, clusters, and superclusters. Phylogenetic analysis reveals nine monophyletic clades and a few phylogenetic orphans. Most clusters contain only genes from the same phylogenetic lineage, reflecting their origin from the exchange of sequence blocks as a result of intralocus recombination. Multiple duplications increased the number of NBS-LRR genes in the progenitors of Arabidopsis, suggesting that the present complexity in Col-0 may derive from as few as 17 progenitors. The combination of physical and phylogenetic analyses of the NBS-LRR genes makes it possible to detect relatively recent gene rearrangements, which increased the number of NBS-LRR genes by about 50, but which are almost never associated with large segmental duplications. The identification of 10 heterogeneous clusters containing members from different clades demonstrates that sequence sampling between different resistance gene loci and clades has occurred. Such events may have taken place early during flowering plant evolution, but they generated modules that have been duplicated and remobilized also more recently.     

Functional characterization and signal transduction ability of nucleotide-binding site-leucine-rich repeat resistance genes in plants

Pathogen infection in plants is often limited by a multifaceted defense response triggered by resistance genes. The most prevalent class of resistance proteins includes those that contain a nucleotide-binding site-leucine-rich repeat (NBS-LRR) domain. Over the past 15 years, more than 50 novel NBS-LRR class resistance genes have been isolated and characterized; they play a significant role in activating conserved defense-signaling networks. Recent molecular research on NBS-LRR resistance proteins and their signaling networks has the potential to broaden the use of resistance genes for disease control. Various transgenic approaches have been tested to broaden the disease resistance spectrum using NBS-LRR genes. This review highlights the recent progress in understanding the structure, function, signal transduction ability of NBS-LRR resistance genes in different host-pathogen systems and suggests new strategies for engineering pathogen resistance in crop plants.     

RPP13 is a simple locus in Arabidopsis thaliana for alleles that specify downy mildew resistance to different avirulence determinants in Peronospora parasitica

Disease resistance (R) genes are found in plants as either simple (single allelic series) loci, or more frequently as complex loci of tandemly repeated genes. These different loci are likely to be under similar evolutionary forces from pathogens, but the contrast between them suggests important differences in mechanisms associated with DNA structure and recombination that generate and maintain R gene diversity. The RPP13 locus in Arabidopsis represents an important paradigm for studying the evolution of an R gene at a simple locus. The RPP13 allele from the accession Nd-1, designated RPP13-Nd, confers resistance to five different isolates of the biotrophic oomycete, Peronospora parasitica (causal agent of downy mildew), and encodes an NBS-LRR type R protein with a putative amino-terminal leucine zipper. The RPP13-Rld allele, cloned from the accession Rld-2, encodes a different specificity. Comparison of three RPP13 alleles revealed a high rate of amino acid divergence within the LRR domain, less than 80% identity overall, compared to the remainder of the protein (> 95% identity). We also found evidence for positive selection in the LRR domain for amino acid diversification outside the core conserved beta-strand/beta-turn motif, suggesting that more of the LRR structure is available for interaction with target molecules than has previously been reported for other R gene products. Furthermore, an amino acid sequence (LLRVLDL) identical in an LRR among RPP13 alleles is conserved in other LZ NBS-LRR type R proteins, suggesting functional significance.     

Isolation of a family of resistance gene analogue sequences of the nucleotide binding site (NBS) type from Lens species.

Most known plant disease-resistance genes (R genes) include in their encoded products domains such as a nucleotide-binding site (NBS) or leucine-rich repeats (LRRs). Sequences with unknown function, but encoding these conserved domains, have been defined as resistance gene analogues (RGAs). The conserved motifs within plant NBS domains make it possible to use degenerate primers and PCR to isolate RGAs. We used degenerate primers deduced from conserved motifs in the NBS domain of NBS-LRR resistance proteins to amplify genomic sequences from Lens species. Fragments from approximately 500-850 bp were obtained. The nucleotide sequence analysis of these fragments revealed 32 different RGA sequences in Lens species with a high similarity (up to 91%) to RGAs from other plants. The predicted amino acid sequences showed that lentil sequences contain all the conserved motifs (P-loop, kinase-2, kinase-3a, GLPL, and MHD) present in the majority of other known plant NBS-LRR resistance genes. Phylogenetic analyses grouped the Lens NBS sequences with the Toll and interleukin-1 receptor (TIR) subclass of NBS-LRR genes, as well as with RGA sequences isolated from other legume species. Using inverse PCR on one putative RGA of lentil, we were able to amplify the flanking regions of this sequence, which contained features found in R proteins.     

Recently duplicated maize R2R3 Myb genes provide evidence for distinct mechanisms of evolutionary divergence after duplication

R2R3 Myb genes are widely distributed in the higher plants and comprise one of the largest known families of regulatory proteins. Here, we provide an evolutionary framework that helps explain the origin of the plant-specific R2R3 Myb genes from widely distributed R1R2R3 Myb genes, through a series of well-established steps. To understand the routes of sequence divergence that followed Myb gene duplication, we supplemented the information available on recently duplicated maize (Zea mays) R2R3 Myb genes (C1/Pl1 and P1/P2) by cloning and characterizing ZmMyb-IF35 and ZmMyb-IF25. These two genes correspond to the recently expanded P-to-A group of maize R2R3 Myb genes. Although the origins of C1/Pl1 and ZmMyb-IF35/ZmMyb-IF25 are associated with the segmental allotetraploid origin of the maize genome, other gene duplication events also shaped the P-to-A clade. Our analyses indicate that some recently duplicated Myb gene pairs display substantial differences in the numbers of synonymous substitutions that have accumulated in the conserved MYB domain and the divergent C-terminal regions. Thus, differences in the accumulation of substitutions during evolution can explain in part the rapid divergence of C-terminal regions for these proteins in some cases. Contrary to previous studies, we show that the divergent C termini of these R2R3 MYB proteins are subject to purifying selection. Our results provide an in-depth analysis of the sequence divergence for some recently duplicated R2R3 Myb genes, yielding important information on general patterns of evolution for this large family of plant regulatory genes.     

Species-specific duplications driving the recent expansion of NBS-LRR genes in five Rosaceae species

Background

Disease resistance (R) genes from different Rosaceae species have been identified by map-based cloning for resistance breeding. However, there are few reports describing the pattern of R-gene evolution in Rosaceae species because several Rosaceae genome sequences have only recently become available.

Results

Since most disease resistance genes encode NBS-LRR proteins, we performed a systematic genome-wide survey of NBS-LRR genes between five Rosaceae species, namely Fragaria vesca (strawberry), Malus × domestica (apple), Pyrus bretschneideri (pear), Prunus persica (peach) and Prunus mume (mei) which contained 144, 748, 469, 354 and 352 NBS-LRR genes, respectively. A high proportion of multi-genes and similar Ks peaks (Ks = 0.1- 0.2) of gene families in the four woody genomes were detected. A total of 385 species-specific duplicate clades were observed in the phylogenetic tree constructed using all 2067 NBS-LRR genes. High percentages of NBS-LRR genes derived from species-specific duplication were found among the five genomes (61.81% in strawberry, 66.04% in apple, 48.61% in pear, 37.01% in peach and 40.05% in mei). Furthermore, the Ks and Ka/Ks values of TIR-NBS-LRR genes (TNLs) were significantly greater than those of non-TIR-NBS-LRR genes (non-TNLs), and most of the NBS-LRRs had Ka/Ks ratios less than 1, suggesting that they were evolving under a subfunctionalization model driven by purifying selection.

Conclusions

Our results indicate that recent duplications played an important role in the evolution of NBS-LRR genes in the four woody perennial Rosaceae species. Based on the phylogenetic tree produced, it could be inferred that species-specific duplication has mainly contributed to the expansion of NBS-LRR genes in the five Rosaceae species. In addition, the Ks and Ka/Ks ratios suggest that the rapidly evolved TNLs have different evolutionary patterns to adapt to different pathogens compared with non-TNL resistant genes.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1291-0) contains supplementary material, which is available to authorized users.     

Frequent,independent transfers of a catabolic gene from bacteria to contrasted filamentous eukaryotes

Even genetically distant prokaryotes can exchange genes between them, and these horizontal gene transfer events play a central role in adaptation and evolution. While this was long thought to be restricted to prokaryotes, certain eukaryotes have acquired genes of bacterial origin. However, gene acquisitions in eukaryotes are thought to be much less important in magnitude than in prokaryotes. Here, we describe the complex evolutionary history of a bacterial catabolic gene that has been transferred repeatedly from different bacterial phyla to stramenopiles and fungi. Indeed, phylogenomic analysis pointed to multiple acquisitions of the gene in these filamentous eukaryotes—as many as 15 different events for 65 microeukaryotes. Furthermore, once transferred, this gene acquired introns and was found expressed in mRNA databases for most recipients. Our results show that effective inter-domain transfers and subsequent adaptation of a prokaryotic gene in eukaryotic cells can happen at an unprecedented magnitude.     

Phylogenetic and evolutionary analysis of NBS-encoding genes in Rutaceae fruit crops

The nucleotide-binding site leucine-rich repeat (NBS-LRR) genes are the largest class of disease resistance genes in plants. However, our understanding of the evolution of NBS-LRR genes in Rutaceae fruit crops is rather limited. We report an evolutionary study of 103 NBS-encoding genes isolated from Poncirus trifoliata (trifoliate orange), Citrus reticulata (tangerine) and their F₁ progeny. In all, 58 of the sequences contained a continuous open reading frame. Phylogenetic analysis classified the 58 NBS genes into nine clades, eight of which were genus specific. This was taken to imply that most of the ancestors of these NBS genes evolved after the genus split. The motif pattern of the 58 NBS-encoding genes was consistent with their phylogenetic profile. An extended phylogenetic analysis, incorporating citrus NBS genes from the public database, classified 95 citrus NBS genes into six clades, half of which were genus specific. RFLP analysis showed that citrus NBS-encoding genes have been evolving rapidly, and that they are unstable when passed through an intergeneric cross. Of 32 NBS-encoding genes tracked by gene-specific PCR, 24 showed segregation distortion among a set of 94 F₁ individuals. This study provides new insight into the evolution of Rutaceae NBS genes and their behaviour following an intergeneric cross.