Origin and maintenance of a broad-spectrum disease resistance locus in Arabidopsis

The broad-spectrum mildew resistance genes RPW8.1 and RPW8.2 define a unique type of plant disease resistance (R) gene, and so far homologous sequences have been found in Arabidopsis thaliana only, which suggests a recent origin. In addition to RPW8.1 and RPW8.2, the RPW8 locus contains three homologs of RPW8, HR1, HR2, and HR3, which do not contribute to powdery mildew resistance. To investigate whether RPW8 has originated recently, and if so the processes involved, we have isolated and analyzed the syntenic RPW8 loci from Arabidopsis lyrata, and from Brassica rapa and B. oleracea. The A. lyrata locus contains four genes orthologous to HR1, HR2, HR3, and RPW8.2, respectively. Two syntenic loci have been characterized in Brassica; one locus contains three genes and is present in both B. oleracea and B. rapa, and the other locus contains a single gene and is detected in B. rapa only. The Brassica homologs have highest similarity to HR3. Sequence analyses suggested that the RPW8 gene family in Brassicaceae originated from an HR3-like ancestor gene through a series of duplications and that RPW8.1 and RPW8.2 evolved from functional diversification through positive selection several MYA. Examination of the sequence polymorphism of 32 A. thaliana accessions at the RPW8 locus and their disease reaction phenotypes revealed that the polymorphic RPW8 locus defines a major source of resistance to powdery mildew diseases. A possible evolutionary mechanism by which functional polymorphism at the AtRPW8 locus has been maintained in contemporary populations of A. thaliana is discussed.     

Intraspecific genetic variations, fitness cost and benefit of RPW8, a disease resistance locus in Arabidopsis thaliana

The RPW8 locus of Arabidopsis thaliana confers broad-spectrum resistance to powdery mildew pathogens. In many A. thaliana accessions, this locus contains two homologous genes, RPW8.1 and RPW8.2. In some susceptible accessions, however, these two genes are replaced by HR4, a homolog of RPW8.1. Here, we show that RPW8.2 from A. lyrata conferred powdery mildew resistance in A. thaliana, suggesting that RPW8.2 might have gained the resistance function before the speciation of A. thaliana and A. lyrata. To investigate how RPW8 has been maintained in A. thaliana, we examined the nucleotide sequence polymorphisms in RPW8 from 51 A. thaliana accessions, related disease reaction phenotypes to the evolutionary history of RPW8.1 and RPW8.2, and identified mutations that confer phenotypic variations. The average nucleotide diversities were high at RPW8.1 and RPW8.2, showing no sign of selective sweep. Moreover, we found that expression of RPW8 incurs fitness benefits and costs on A. thaliana in the presence and absence of the pathogens, respectively. Our results suggest that polymorphisms at the RPW8 locus in A. thaliana may have been maintained by complex selective forces, including those from the fitness benefits and costs both associated with RPW8.     

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.     

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.     

Linkage disequilibrium between incompatibility locus region genes in the plant Arabidopsis lyrata

We have studied diversity in Arabidopsis lyrata of sequences orthologous to the ARK3 gene of A. thaliana. Our main goal was to test for recombination in the S-locus region. In A. thaliana, the single-copy ARK3 gene is closely linked to the non-functional copies of the self-incompatibility loci, and the ortholog in A. lyrata (a self-incompatible species) is in the homologous genome region and is known as Aly8. It is thus of interest to test whether Aly8 sequence diversity is elevated due to close linkage to the highly polymorphic incompatibility locus, as is theoretically predicted. However, Aly8 is not a single-copy gene, and the presence of paralogs could also lead to the appearance of elevated diversity. We established a typing approach based on different lengths of Aly8 PCR products and show that most A. lyrata haplotypes have a single copy, but some have two gene copies, both closely linked to the incompatibility locus, one being a pseudogene. We determined the phase of multiple haplotypes in families of plants from Icelandic and other populations. Different Aly8 sequence types are associated with different SRK alleles, while haplotypes with the same SRK sequences tend to have the same Aly8 sequence. There is evidence of some exchange of sequences between different Aly8 sequences, making it difficult to determine which ones are allelic or to estimate the diversity. However, the homogeneity of the Aly8 sequences of each S-haplotype suggests that recombination between the loci has been very infrequent over the evolutionary history of these populations. Overall, the results suggest that recombination rarely occurs in the interval between the S-loci and Aly8 and that linkage to the S-loci can probably account for the observed high Aly8 diversity.     

A cluster of disease resistance genes in Arabidopsis is coordinately regulated by transcriptional activation and RNA silencing

The RPP5 (for recognition of Peronospora parasitica 5) locus in the Arabidopsis thaliana Columbia strain contains a cluster of paralogous disease Resistance (R) genes that play important roles in innate immunity. Among the R genes in this locus, RPP4 confers resistance to two races of the fungal pathogen Hyaloperonospora parasitica, while activation of SNC1 (for suppressor of npr1-1, constitutive 1) results in the resistance to another race of H. parasitica and to pathovars of the bacterial pathogen Pseudomonas syringae through the accumulation of salicylic acid (SA). Here, we demonstrate that other Columbia RPP5 locus R genes can be induced by transgenic overexpression of SNC1, which itself is regulated by a positive amplification loop involving SA accumulation. We also show that small RNA species that can target RPP5 locus R genes are produced in wild-type plants and that these R genes can be cosuppressed in transgenic plants overexpressing SNC1. Steady state expression levels of SNC1 increase in some mutants (dcl4-4, ago1-36, and upf1-5) defective in RNA silencing as well as in transgenic plants expressing the P1/Helper Component-Protease viral suppressor of RNA silencing. However, steady state levels of small RNA species do not change in mutants that upregulate SNC1. These data indicate many Columbia RPP5 locus R genes can be coordinately regulated both positively and negatively and suggest that the RPP5 locus is poised to respond to pathogens that disturb RNA silencing.     

Genome evolution among cruciferous plants: a lecture from the comparison of the genetic maps of three diploid species--Capsella rubella, Arabidopsis lyrata subsp. petraea, and A. thaliana

Comparative mapping in cruciferous plants is ongoing, and recently two additional genetic maps of diploid Capsella and Arabidopsis lyrata subsp. petraea have been presented. We compared both maps with each other using the sequence map and genomic data resources from Arabidopsis thaliana as a reference. The ancestors of the species pair Capsella-Arabidopsis diverged from one another approximately 10-14 million years ago (mya), whereas Arabidopsis thaliana and Arabidopsis lyrata have been separated since roughly 5-6 mya. Our analysis indicated that among diploid Capsella and Arabidopsis lyrata all eight genetic linkage groups are totally colinear to each other, with only two inversions significantly differentiating these two species.By minimizing the number of chromosomal rearrangements during genome evolution, we presented a model of chromosome evolution involving all three species. From this scenario, it is obvious that Arabidopsis thaliana underwent a dramatic genome reconstruction, with a base chromosome number reduction from five to eight and with approximately 1.3 chromosomal rearrangements per million years. In contrast, the terminal lineage leading to Capsella has only undergone less than 0.09 rearrangements per million years. This is the same rate as calculated for Arabidopsis lyrata since its separation from the Capsella lineage 10-14 mya. These results are in strong contrast to all overestimated rates calculated from comparisons of the systems Arabidopsis thaliana and Brassica, and our data demonstrate the problematic nature of both model systems.     

Maintenance of genetic variation in plants and pathogens involves complex networks of gene-for-gene interactions

The RPP13 [recognition of Hyaloperonospora arabidopsidis (previously known as Peronospora parasitica )] resistance ( R ) gene in Arabidopsis thaliana exhibits the highest reported level of sequence diversity among known R genes. Consistent with a co-evolutionary model, the matching effector protein ATR13 ( A. thaliana -recognized) from H. arabidopsidis reveals extreme levels of allelic diversity. We isolated 23 new RPP13 sequences from a UK metapopulation, giving a total of 47 when combined with previous studies. We used these in functional studies of the A. thaliana accessions for their resistance response to 16 isolates of H. arabidopsidis . We characterized the molecular basis of recognition by the expression of the corresponding ATR13 genes from these 16 isolates in these host accessions. This allowed the determination of which alleles of RPP13 were responsible for pathogen recognition and whether recognition was dependent on the RPP13/ATR13 combination. Linking our functional studies with phylogenetic analysis, we determined that: (i) the recognition of ATR13 is mediated by alleles in just a single RPP13 clade; (ii) RPP13 alleles in other clades have evolved the ability to detect other pathogen ATR protein(s); and (iii) at least one gene, unlinked to RPP13 in A. thaliana, detects a different subgroup of ATR13 alleles.     

Contrasting patterns of nucleotide polymorphism at the alcohol dehydrogenase locus in the outcrossing Arabidopsis lyrata and the selfing Arabidopsis thaliana

Nucleotide variation at the alcohol dehydrogenase locus (Adh) was studied in the outcrossing Arabidopsis lyrata, a close relative of the selfing Arabidopsis thaliana. Overall, estimated nucleotide diversity in the North American ssp. lyrata and two European ssp. petraea populations was 0.0038, lower than the corresponding specieswide estimate for A. thaliana at the same set of nucleotide sites. The distribution of segregating sites across the gene differed between the two species. Estimated sequence diversity within an A. lyrata population with a large sample size (0.0023) was much higher than has previously been observed for A. thaliana. This North American population has an excess of sites at intermediate frequencies compared with neutral expectation (Tajima's D = 2.3, P < 0.005), suggestive of linked balancing selection or a recent population bottleneck. In contrast, an excess of rare polymorphisms has been found in A. thaliana. Polymorphism within A. lyrata and divergence from A. thaliana appear to be correlated across the Adh gene sequence. The geographic distribution of polymorphism was quite different from that of A. thaliana, for which earlier studies of several genes found low within-population nucleotide site polymorphism and no overall continental differentiation of variation despite large differences in site frequencies between local populations. Differences between the outcrossing A. lyrata and the selfing A. thaliana reflect the impact of differences in mating system and the influence of bottlenecks in A. thaliana during rapid colonization on DNA sequence polymorphism. The influence of additional variability-reducing mechanisms, such as background selection or hitchhiking, may not be discernible.     

Members of the Arabidopsis HRT/RPP8 family of resistance genes confer resistance to both viral and oomycete pathogens

Turnip crinkle virus (TCV) inoculation onto TCV-resistant Arabidopsis leads to a hypersensitive response (HR) controlled by the dominant gene HRT. HRT is a member of the class of resistance (R) genes that contain a leucine zipper, a nucleotide binding site, and leucine-rich repeats. The chromosomal position of HRT and its homology to resistance gene RPP8 and two RPP8 homologs indicate that unequal crossing over and gene conversion may have contributed to HRT evolution. RPP8 confers resistance to an oomycete pathogen, Peronospora parasitica. Despite very strong similarities within the HRT/RPP8 family, HRT and RPP8 are specific for the respective pathogens they detect. Hence, the HRT/RPP8 family provides molecular evidence that sequence changes between closely related members of multigene families can generate novel specificities for radically different pathogens. Transgenic plants expressing HRT developed an HR but generally remained susceptible to TCV because of a second gene, RRT, that regulates resistance to TCV. However, several transgenic plants that overexpressed HRT produced micro-HRs or no HR when inoculated with TCV and were resistant to infection. Expression of the TCV coat protein gene in seedlings containing HRT resulted in massive necrosis and death, indicating that the avirulence factor detected by the HRT-encoded protein is the TCV coat protein.     

Further insights into the phylogeny of Arabidopsis (Brassicaceae) from nuclear Atmyb2 flanking sequence

Arabidopsis thaliana is the preeminent plant model organism. However, significant advances in evolution and ecology are being made by expanding the scope of research beyond this single species into the broader genus Arabidopsis. Surprisingly, few studies have rigorously investigated phylogenetic relationships between the nine Arabidopsis species, and this study evaluates both these and hypotheses related to two instances of intra-generic hybridization. DNA sequences from the 5' flanking region of the nuclear Atmyb2 gene from 12 of the 14 Arabidopsis taxa were used to reconstruct the generic phylogeny. The strict consensus tree was highly concordant with previous studies, identifying lineages corresponding to widespread species but exhibiting a large basal polytomy. Our data indicates that the paternal parent of the allopolyploid A. suecica is A. neglecta rather than A. arenosa s.l., although the need for a detailed phylogeographical study of these three species is noted. Finally, our data provided additional phylogenetic evidence of hybridization between Arabidopsis lyrata s.l. and A. halleri s.l. Taken together, the well-defined lineages within the genus and the potential for hybridization between them highlight Arabidopsis as a promising group for comparative and experimental studies of hybridization.     

White rust (Albugo candida) resistance loci on three Arabidopsis chromosomes are closely linked to downy mildew (Peronospora parasitica) resistance loci

Two accessions of Arabidopsis thaliana (Ksk-1 and Ksk-2) were used to identify and map three loci ( RAC1 , RAC2 and RAC3 ) of genes that confer R esistance to A lbugo c andida (white rust). The phenotypes associated with these genes were classified as either FN (necrotic f lecks on upper surface of cotyledons and n o blisters) for RAC2 and RAC3, or FYN ( f lecks surrounded by y ellowing and n o blisters) for RAC1 . Both phenotypes exhibited rapid death of host cells penetrated by the parasite (hypersensitive response), with callose deposition commonly encasing the haustorium. F₆ recombinant inbred lines were produced specifically for the purpose of mapping each RAC locus relative to molecular markers. Dominant resistance at the locus RAC1 in Ksk-1 was previously mapped to chromosome 1 between RFLP markers m253 and m254, and co-segregating with a downy mildew resistance specificity RPP9 in the accession Wei-0. We report here a fine-scale map interval and co-segregating markers for this locus, which in turn enabled mapping of a previously unnoticed source of resistance in Ksk-1 designated RAC3 that exhibits an FN phenotype hyperstatic to the FYN phenotype of RAC1. RAC3 is closely linked to the RPP8/HRT on chromosome 5, a locus which contains specificities for resistance to downy mildew and turnip crinkle virus. Recombinant inbreds also enabled mapping of recessive resistance at RAC2 in Ksk-2 to the bottom arm of chromosome 3, in the 6 cM interval between two downy mildew resistance loci ( RPP1 and RPP13 ) .     

Self-incompatibility in the genus Arabidopsis: characterization of the S locus in the outcrossing A. lyrata and its autogamous relative A. thaliana

As a starting point for a phylogenetic study of self-incompatibility (SI) in crucifers and to elucidate the genetic basis of transitions between outcrossing and self-fertilizing mating systems in this family, we investigated the SI system of Arabidopsis lyrata. A. lyrata is an outcrossing close relative of the self-fertile A. thaliana and is thought to have diverged from A. thaliana approximately 5 million years ago and from Brassica spp 15 to 20 million years ago. Analysis of two S (sterility) locus haplotypes demonstrates that the A. lyrata S locus contains tightly linked orthologs of the S locus receptor kinase (SRK) gene and the S locus cysteine-rich protein (SCR) gene, which are the determinants of SI specificity in stigma and pollen, respectively, but lacks an S locus glycoprotein gene. As described previously in Brassica, the S haplotypes of A. lyrata differ by the rearranged order of their genes and by their variable physical sizes. Comparative mapping of the A. lyrata and Brassica S loci indicates that the S locus of crucifers is a dynamic locus that has undergone several duplication events since the Arabidopsis--Brassica split and was translocated as a unit between two distant chromosomal locations during diversification of the two taxa. Furthermore, comparative analysis of the S locus region of A. lyrata and its homeolog in self-fertile A. thaliana identified orthologs of the SRK and SCR genes and demonstrated that self-compatibility in this species is associated with inactivation of SI specificity genes.     

Comparative gene mapping in Arabidopsis lyrata chromosomes 1 and 2 and the corresponding A. thaliana chromosome 1: recombination rates, rearrangements and centromere location

To add detail to the genetic map of Arabidopsis lyrata, and compare it with that of A. thaliana, we have developed many additional markers in the A. lyrata linkage groups, LG1 and LG2, corresponding to A. thaliana chromosome 1. We used a newly developed method for marker development for single nucleotide polymorphisms present in gene sequences, plus length differences, to map genes in an A. lyrata family, including variants in several genes close to the A. thaliana centromere 1, providing the first data on the location of an A. lyrata centromere; we discuss the implications for the evolution of chromosome 1 of A. thaliana. With our larger marker density, large rearrangements between the two Arabidopsis species are excluded, except for a large inversion on LG2. This was previously known in Capsella; its presence in A. lyrata suggests that, like most other rearrangements, it probably arose in the A. thaliana lineage. Knowing that marker orders are similar, we can now compare homologous, non-rearranged map distances to test the prediction of more frequent crossing-over in the more inbreeding species. Our results support the previous conclusion of similar distances in the two species for A. lyrata LG1 markers. For LG2 markers, the distances were consistently, but non-significantly, larger in A. lyrata. Given the two species' large DNA content difference, the similarity of map lengths, particularly for LG1, suggests that crossing-over is more frequent across comparable physical distances in the inbreeder, A. thaliana, as predicted.     

The transition to self-compatibility in Arabidopsis thaliana and evolution within S-haplotypes over 10 Myr

A recent investigation found evidence that the transition of Arabidopsis thaliana from ancestral self-incompatibility (SI) to full self-compatibility occurred very recently and suggested that this occurred through a selective fixation of a nonfunctional allele (PsiSCR1) at the SCR gene, which determines pollen specificity in the incompatibility response. The main evidence is the lack of polymorphism at the SCR locus in A. thaliana. However, the nearby SRK gene, which determines stigma specificity in self-incompatible Brassicaceae species, has extremely high sequence diversity, with 3 very divergent SRK haplotypes, 2 of them present in multiple strains. Such high diversity is extremely unusual in this species, and it suggests the possibility that multiple, different SRK haplotypes may have been preserved from A. thaliana's self-incompatible ancestor. To study the evolution of S-haplotypes in the A. thaliana lineage, we searched the 2 most closely related Arabidopsis species Arabidopsis lyrata and Arabidopsis halleri, in which most populations have retained SI, and found SRK sequences corresponding to all 3 A. thaliana haplogroup sequences. Our molecular evolutionary analyses of these 3 S-haplotypes provide an independent estimate of the timing of the breakdown of SI and again exclude an ancient transition to selfing in A. thaliana. Comparing sequences of each of the 3 haplogroups between species, we find that 2 of the 3 SRK sequences (haplogroups A and B) are similar throughout their length, suggesting that little or no recombination with other SRK alleles has occurred since these species diverged. The diversity difference between the SCR and SRK loci in A. thaliana, however, suggests crossing-over, either within SRK or between the SCR and SRK loci. If the loss of SI involved fixation of the PsiSCR1 sequence, the exchange must have occurred during its fixation. Divergence between the species is much lower at the S-locus, compared with reference loci, and we discuss two contributory possibilities. Introgression may have occurred between A. lyrata and A. halleri and between their ancestral lineage and A. thaliana, at least for some period after their split. In addition, the coalescence times of sequences of individual S-haplogroups are expected to be less than those of alleles at non-S-loci.     

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.