Evolution of plant resistance at the molecular level: ecological context of species interactions

Cucumber mosaic virus infection as a potential selective pressure on Arabidopsis thaliana populations

It has been proposed that in wild ecosystems viruses are often plant mutualists, whereas agroecosystems favour pathogenicity. We seek evidence for virus pathogenicity in wild ecosystems through the analysis of plant-virus coevolution, which requires a negative effect of infection on the host fitness. We focus on the interaction between Arabidopsis thaliana and Cucumber mosaic virus (CMV), which is significant in nature. We studied the genetic diversity of A. thaliana for two defence traits, resistance and tolerance, to CMV. A set of 185 individuals collected in 76 A. thaliana Iberian wild populations were inoculated with different CMV strains. Resistance was estimated from the level of virus multiplication in infected plants, and tolerance from the effect of infection on host progeny production. Resistance and tolerance to CMV showed substantial genetic variation within and between host populations, and depended on the virus x host genotype interaction, two conditions for coevolution. Resistance and tolerance were co-occurring independent traits that have evolved independently from related life-history traits involved in adaptation to climate. The comparison of the genetic structure for resistance and tolerance with that for neutral traits (QST/FST analyses) indicated that both defence traits are likely under uniform selection. These results strongly suggest that CMV infection selects for defence on A. thaliana populations, and support plant-virus coevolution. Thus, we propose that CMV infection reduces host fitness under the field conditions of the wild A. thaliana populations studied.

Assortment of Flowering Time and Immunity Alleles in Natural Arabidopsis thaliana Populations Suggests Immunity and Vegetative Lifespan Strategies Coevolve

The selective impact of pathogen epidemics on host defenses can be strong but remains transient. By contrast, life-history shifts can durably and continuously modify the balance between costs and benefits of immunity, which arbitrates the evolution of host defenses. Their impact on the evolutionary dynamics of host immunity, however, has seldom been documented. Optimal investment into immunity is expected to decrease with shortening lifespan, because a shorter life decreases the probability to encounter pathogens or enemies. Here, we document that in natural populations of Arabidopsis thaliana, the expression levels of immunity genes correlate positively with flowering time, which in annual species is a proxy for lifespan. Using a novel genetic strategy based on bulk-segregants, we partitioned flowering time-dependent from -independent immunity genes and could demonstrate that this positive covariation can be genetically separated. It is therefore not explained by the pleiotropic action of some major regulatory genes controlling both immunity and lifespan. Moreover, we find that immunity genes containing variants reported to impact fitness in natural field conditions are among the genes whose expression covaries most strongly with flowering time. Taken together, these analyses reveal that natural selection has likely assorted alleles promoting lower expression of immunity genes with alleles that decrease the duration of vegetative lifespan in A. thaliana and vice versa. This is the first study documenting a pattern of variation consistent with the impact that selection on flowering time is predicted to have on diversity in host immunity.

Quantitative proteomics identifies 38 proteins that are differentially expressed in cucumber in response to cucumber green mottle mosaic virus infection

BACKGROUND

Since it was first reported in 1935, Cucumber green mottle mosaic virus (CGMMV) has become a serious pathogen in a range of cucurbit crops. The virus is generally transmitted by propagation materials, and to date no effective chemical or cultural methods of control have been developed to combat its spread. The current study presents a preliminary analysis of the pathogenic mechanisms from the perspective of protein expression levels in an infected cucumber host, with the objective of elucidating the infection process and potential strategies to reduce both the economic and yield losses associated with CGMMV.

METHODS

Isobaric tags for relative and absolute quantitation (iTRAQ) technology coupled with liquid chromatography-tandem mass spectrometric (LC-MS/MS) were used to identify the differentially expressed proteins in cucumber plants infected with CGMMV compared with mock-inoculated plants. The functions of the proteins were deduced by functional annotation and their involvement in metabolic processes explored by KEGG pathway analysis to identify their interactions during CGMMV infection, while their in vivo expression was further verified by qPCR.

RESULTS

Infection by CGMMV altered both the expression level and absolute quantity of 38 proteins (fold change >0.6) in cucumber hosts. Of these, 23 were found to be up-regulated, while 15 were down-regulated. Gene ontology (GO) analysis revealed that 22 of the proteins had a combined function and were associated with molecular function (MF), biological process (BP) and cellular component (CC). Several other proteins had a dual function with 1, 7, and 2 proteins being associated with BP/CC, BP/MF, CC/MF, respectively. The remaining 3 proteins were only involved in MF. In addition, Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis identified 18 proteins that were involved in 13 separate metabolic pathways. These pathways were subsequently merged to generate three network diagrams illustrating the interactions between the different pathways, while qPCR was used to track the changes in expression levels of the proteins identified at 3 time points during CGMMV infection. Taken together these results greatly expand our understanding of the relationships between CGMMV and cucumber hosts.

CONCLUSIONS

The results of the study indicate that CGMMV infection significantly changes the physiology of cucumbers, affecting the expression levels of individual proteins as well as entire metabolic pathways. The bioinformatic analysis also identified several pathogenesis-related (PR) proteins that could be useful in the development of disease-resistant plants.

Selection, genome-wide fitness effects and evolutionary rates in the model legume Medicago truncatula

Sequence data for >20 000 annotated genes from 56 accessions of Medicago truncatula were used to identify potential targets of positive selection, the determinants of evolutionary rate variation and the relative importance of positive and purifying selection in shaping nucleotide diversity. Based upon patterns of intraspecific diversity and interspecific divergence, c. 50-75% of nonsynonymous polymorphisms are subject to strong purifying selection and 1% of the sampled genes harbour a signature of positive selection. Combining polymorphism with expression data, we estimated the distribution of fitness effects and found that the proportion of deleterious mutations is significantly greater for expressed genes than for genes with undetected transcripts (nonexpressed) in a previous RNA-seq experiment and greater for broadly expressed genes than those expressed in only a single tissue. Expression level is the strongest correlate of evolutionary rates at nonsynonymous sites, and despite multiple genomic features being significantly correlated with evolutionary rates, they explain less than 20% of the variation in nonsynonymous rates (dN) and <15% of the variation in either synonymous rates (dS) or dN:dS. Among putative targets of selection were genes involved in defence against pathogens and herbivores, genes with roles in mediating the relationship with rhizobial symbionts and one-third of annotated histone-lysine methyltransferases. Adaptive evolution of the methyltransferases suggests that positive selection in gene expression may have occurred through evolution of enzymes involved in epigenetic modification.

Adaptive molecular evolution of a defence gene in sexual but not functionally asexual evening primroses

Theory predicts that sexual reproduction provides evolutionary advantages over asexual reproduction by reducing mutational load and increasing adaptive potential. Here, we test the latter prediction in the context of plant defences against pathogens because pathogens frequently reduce plant fitness and drive the evolution of plant defences. Specifically, we ask whether sexual evening primrose plant lineages (Onagraceae) have faster rates of adaptive molecular evolution and altered gene expression of a class I chitinase, a gene implicated in defence against pathogens, than functionally asexual evening primrose lineages. We found that the ratio of amino acid to silent substitutions (K(a) /K(s) = 0.19 vs. 0.11 for sexual and asexual lineages, respectively), the number of sites identified to be under positive selection (four vs. zero for sexual and asexual lineages, respectively) and the expression of chitinase were all higher in sexual than in asexual lineages. Our results are congruent with the conclusion that a loss of sexual recombination and segregation in the Onagraceae negatively affects adaptive structural and potentially regulatory evolution of a plant defence protein.

Molecular population genetics of elicitor-induced resistance genes in European aspen (Populus tremula L., Salicaceae)

Owing to their long life span and ecological dominance in many communities, forest trees are subject to attack from a diverse array of herbivores throughout their range, and have therefore developed a large number of both constitutive and inducible defenses. We used molecular population genetics methods to examine the evolution of eight genes in European aspen, Populus tremula, that are all associated with defensive responses against pests and/or pathogens, and have earlier been shown to become strongly up-regulated in poplars as a response to wounding and insect herbivory. Our results show that the majority of these defense genes show patterns of intraspecific polymorphism and site-frequency spectra that are consistent with a neutral model of evolution. However, two of the genes, both belonging to a small gene family of polyphenol oxidases, show multiple deviations from the neutral model. The gene PPO1 has a 600 bp region with a highly elevated K(A)/K(S) ratio and reduced synonymous diversity. PPO1 also shows a skew toward intermediate frequency variants in the SFS, and a pronounced fixation of non-synonymous mutations, all pointing to the fact that PPO1 has been subjected to recurrent selective sweeps. The gene PPO2 shows a marked excess of high frequency, derived variants and shows many of the same trends as PPO1 does, even though the pattern is less pronounced, suggesting that PPO2 might have been the target of a recent selective sweep. Our results supports data from both Populus and other species which have found that the the majority of defense-associated genes show few signs of selection but that a number of genes involved in mediating defense against herbivores show signs of adaptive evolution.

Coexistence of trichome variation in a natural plant population: a combined study using ecological and candidate gene approaches

The coexistence of distinct phenotypes within populations has long been investigated in evolutionary ecology. Recent studies have identified the genetic basis of distinct phenotypes, but it is poorly understood how the variation in candidate loci is maintained in natural environments. In this study, we examined fitness consequences and genetic basis of variation in trichome production in a natural population of Arabidopsis halleri subsp. gemmifera. Half of the individuals in the study population produced trichomes while the other half were glabrous, and the leaf beetle Phaedon brassicae imposed intensive damage to both phenotypes. The fitness of hairy and glabrous plants showed no significant differences in the field during two years. A similar result was obtained when sibling hairy and glabrous plants were transplanted at the same field site, whereas a fitness cost of trichome production was detected under a weak herbivory condition. Thus, equivalent fitness of hairy and glabrous plants under natural herbivory allows their coexistence in the contemporary population. The pattern of polymorphism of the candidate trichome gene GLABROUS1 (GL1) showed no evidence of long-term maintenance of trichome variation within the population. Although balancing selection under fluctuating biotic environments is often proposed to explain the maintenance of defense variation, the lack of clear evidence of balancing selection in the study population suggests that other factors such as gene flow and neutral process may have played relatively large roles in shaping trichome variation at least for the single population level.

Effect of hosts on competition among clones and evidence of differential selection between pathogenic and saprophytic phases in experimental populations of the wheat pathogen Phaeosphaeria nodorum

BACKGROUND

Monoculture, multi-cropping and wider use of highly resistant cultivars have been proposed as mechanisms to explain the elevated rate of evolution of plant pathogens in agricultural ecosystems. We used a mark-release-recapture experiment with the wheat pathogen Phaeosphaeria nodorum to evaluate the impact of two of these mechanisms on the evolution of a pathogen population. Nine P. nodorum isolates marked with ten microsatellite markers and one minisatellite were released onto five replicated host populations to initiate epidemics of Stagonospora nodorum leaf blotch. The experiment was carried out over two consecutive host growing seasons and two pathogen collections were made during each season.

RESULTS

A total of 637 pathogen isolates matching the marked inoculants were recovered from inoculated plots over two years. Genetic diversity in the host populations affected the evolution of the corresponding P. nodorum populations. In the cultivar mixture the relative frequencies of inoculants did not change over the course of the experiment and the pathogen exhibited a low variation in selection coefficients.

CONCLUSIONS

Our results support the hypothesis that increasing genetic heterogeneity in host populations may retard the rate of evolution in associated pathogen populations. Our experiment also provides indirect evidence of fitness costs associated with host specialization in P. nodorum as indicated by differential selection during the pathogenic and saprophytic phases.

Spatial variation in disease resistance: from molecules to metapopulations

Variation in disease resistance is a widespread phenomenon in wild plant-pathogen associations. Here, we review current literature on natural plant-pathogen associations to determine how diversity in disease resistance is distributed at different hierarchical levels - within host individuals, within host populations, among host populations at the metapopulation scale and at larger regional scales.We find diversity in resistance across all spatial scales examined. Furthermore, variability seems to be the best counter-defence of plants against their rapidly evolving pathogens. We find that higher diversity of resistance phenotypes also results in higher levels of resistance at the population level.Overall, we find that wild plant populations are more likely to be susceptible than resistant to their pathogens. However, the degree of resistance differs strikingly depending on the origin of the pathogen strains used in experimental inoculation studies. Plant populations are on average 16% more resistant to allopatric pathogen strains than they are to strains that occur within the same population (48 % vs. 32 % respectively).Pathogen dispersal mode affects levels of resistance in natural plant populations with lowest levels detected for hosts of airborne pathogens and highest for waterborne pathogens.Detailed analysis of two model systems, Linum marginale infected by Melampsora lini, and Plantago lanceolata infected by Podosphaera plantaginis, show that the amount of variation in disease resistance declines towards higher spatial scales as we move from individual hosts to metapopulations, but evaluation of multiple spatial scales is needed to fully capture the structure of disease resistance.Synthesis: Variation in disease resistance is ubiquitous in wild plant-pathogen associations. While the debate over whether the resistance structure of plant populations is determined by pathogen-imposed selection versus non-adaptive processes remains unresolved, we do report examples of pathogen-imposed selection on host resistance. Here we highlight the importance of measuring resistance across multiple spatial scales, and of using sympatric strains when looking for signs of coevolution in wild plant-pathogen interactions.

Evolution of disease response genes in loblolly pine: insights from candidate genes

Background

Host-pathogen interactions that may lead to a competitive co-evolution of virulence and resistance mechanisms present an attractive system to study molecular evolution because strong, recent (or even current) selective pressure is expected at many genomic loci. However, it is unclear whether these selective forces would act to preserve existing diversity, promote novel diversity, or reduce linked neutral diversity during rapid fixation of advantageous alleles. In plants, the lack of adaptive immunity places a larger burden on genetic diversity to ensure survival of plant populations. This burden is even greater if the generation time of the plant is much longer than the generation time of the pathogen.

Methodology/Principal Findings

Here, we present nucleotide polymorphism and substitution data for 41 candidate genes from the long-lived forest tree loblolly pine, selected primarily for their prospective influences on host-pathogen interactions. This dataset is analyzed together with 15 drought-tolerance and 13 wood-quality genes from previous studies. A wide range of neutrality tests were performed and tested against expectations from realistic demographic models.

Conclusions/Significance

Collectively, our analyses found that axr (auxin response factor), caf1 (chromatin assembly factor) and gatabp1 (gata binding protein 1) candidate genes carry patterns consistent with directional selection and erd3 (early response to drought 3) displays patterns suggestive of a selective sweep, both of which are consistent with the arm-race model of disease response evolution. Furthermore, we have identified patterns consistent with diversifying selection at erf1-like (ethylene responsive factor 1), ccoaoemt (caffeoyl-CoA-O-methyltransferase), cyp450-like (cytochrome p450-like) and pr4.3 (pathogen response 4.3), expected under the trench-warfare evolution model. Finally, a drought-tolerance candidate related to the plant cell wall, lp5, displayed patterns consistent with balancing selection. In conclusion, both arms-race and trench-warfare models seem compatible with patterns of polymorphism found in different disease-response candidate genes, indicating a mixed strategy of disease tolerance evolution for loblolly pine, a major tree crop in southeastern United States.

Genomic tools development for Aquilegia: construction of a BAC-based physical map

BACKGROUND

The genus Aquilegia, consisting of approximately 70 taxa, is a member of the basal eudicot lineage, Ranuculales, which is evolutionarily intermediate between monocots and core eudicots, and represents a relatively unstudied clade in the angiosperm phylogenetic tree that bridges the gap between these two major plant groups. Aquilegia species are closely related and their distribution covers highly diverse habitats. These provide rich resources to better understand the genetic basis of adaptation to different pollinators and habitats that in turn leads to rapid speciation. To gain insights into the genome structure and facilitate gene identification, comparative genomics and whole-genome shotgun sequencing assembly, BAC-based genomics resources are of crucial importance.

RESULTS

BAC-based genomic resources, including two BAC libraries, a physical map with anchored markers and BAC end sequences, were established from A. formosa. The physical map was composed of a total of 50,155 BAC clones in 832 contigs and 3939 singletons, covering 21X genome equivalents. These contigs spanned a physical length of 689.8 Mb (~2.3X of the genome) suggesting the complex heterozygosity of the genome. A set of 197 markers was developed from ESTs induced by drought-stress, or involved in anthocyanin biosynthesis or floral development, and was integrated into the physical map. Among these were 87 genetically mapped markers that anchored 54 contigs, spanning 76.4 Mb (25.5%) across the genome. Analysis of a selection of 12,086 BAC end sequences (BESs) from the minimal tiling path (MTP) allowed a preview of the Aquilegia genome organization, including identification of transposable elements, simple sequence repeats and gene content. Common repetitive elements previously reported in both monocots and core eudicots were identified in Aquilegia suggesting the value of this genome in connecting the two major plant clades. Comparison with sequenced plant genomes indicated a higher similarity to grapevine (Vitis vinifera) than to rice and Arabidopsis in the transcriptomes.

CONCLUSIONS

The A. formosa BAC-based genomic resources provide valuable tools to study Aquilegia genome. Further integration of other existing genomics resources, such as ESTs, into the physical map should enable better understanding of the molecular mechanisms underlying adaptive radiation and elaboration of floral morphology.

Natural variation in priming of basal resistance: from evolutionary origin to agricultural exploitation

Biotic stress has a major impact on the process of natural selection in plants. As plants have evolved under variable environmental conditions, they have acquired a diverse spectrum of defensive strategies against pathogens and herbivores. Genetic variation in the expression of plant defence offers valuable insights into the evolution of these strategies. The 'zigzag' model, which describes an ongoing arms race between inducible plant defences and their suppression by pathogens, is now a commonly accepted model of plant defence evolution. This review explores additional strategies by which plants have evolved to cope with biotic stress under different selective circumstances. Apart from interactions with plant-beneficial micro-organisms that can antagonize pathogens directly, plants have the ability to prime their immune system in response to selected environmental signals. This defence priming offers disease protection that is effective against a broad spectrum of virulent pathogens, as long as the augmented defence reaction is expressed before the invading pathogen has the opportunity to suppress host defences. Furthermore, priming has been shown to be a cost-efficient defence strategy under relatively hostile environmental conditions. Accordingly, it is possible that selected plant varieties have evolved a constitutively primed immune system to adapt to levels of disease pressure. Here, we examine this hypothesis further by evaluating the evidence for natural variation in the responsiveness of basal defence mechanisms, and discuss how this genetic variation can be exploited in breeding programmes to provide sustainable crop protection against pests and diseases.

Ecological genetics and genomics of plant defenses: Evidence and approaches

Herbivores exert significant selection on plants, and plants have evolved a variety of constitutive and inducible defenses to resist and tolerate herbivory. Assessing the genetic mechanisms that influence defenses against herbivores will deepen our understanding of the evolution of essential phenotypic traits.Ecogenomics is a powerful interdisciplinary approach that can address fundamental questions about the ecology and evolutionary biology of species, such as: which evolutionary forces maintain variation within a population? and What is the genetic architecture of adaptation? This field seeks to identify gene regions that influence ecologically-important traits, assess the fitness consequences under natural conditions of alleles at key quantitative trait loci (QTLs), and test how the abiotic and biotic environment affects gene expression.Here, we review ecogenomics techniques and emphasize how this framework can address long-standing and emerging questions relating to anti-herbivore defenses in plants. For example, ecogenomics tools can be used to investigate: inducible vs. constitutive defenses; tradeoffs between resistance and tolerance; adaptation to the local herbivore community; selection on alleles that confer resistance and tolerance in natural populations; and whether different genes are activated in response to specialist vs. generalist herbivores and to different types of damage.Ecogenomic studies can be conducted with model species, such as Arabidopsis, or their relatives, in which case myriad molecular tools are already available. Burgeoning sequence data will also facilitate ecogenomic studies of non-model species. Throughout this paper, we highlight approaches that are particularly suitable for ecological studies of non-model organisms, discuss the benefits and disadvantages of specific techniques, and review bioinformatic tools for analyzing data.We focus on established and promising techniques, such as QTL mapping with pedigreed populations, genome wide association studies, transcription profiling strategies, population genomics, and transgenic methodologies. Many of these techniques are complementary and can be used jointly to investigate the genetic architecture of defense traits and selection on alleles in nature.

Population-level variation of the preproricin gene contradicts expectation of neutral equilibrium for generalist plant defense toxins

The preproricin gene encodes ricin, the highly toxic, type II ribosome-inactivating protein of castor bean (Ricinus communis L.). As a generalist plant defense gene, preproricin is expected to exhibit population-level variation consistent with the neutral equilibrium model and to comprise few functionally different alleles. We first test the hypothesis that the preproricin gene family should comprise six to eight members by searching the publicly available draft genome sequence of R. communis and analyzing its ricin-like loci. We then test the neutral equilibrium expectation for the preproricin gene by characterizing its allelic variation among 25 geographically diverse castor bean plants. We confirm the presence of six ricin-like loci that share with the preproricin gene 62.9-96.3% nucleotide identity and intact A-chains. DNA sequence variation among the preproricin haplotypes significantly rejects tests of the neutral equilibrium model. Replacement mutations preserve the 12 amino acids known to affect catalytic and electrostatic interactions of the native protein toxin, which suggests functional divergence among alleles has been minimal. Nucleotide polymorphism is maintained by purifying selection (omega < 0.3) yet includes an excess of rare silent mutations greater than predicted by the neutral equilibrium model. Development of robust detection methods for ricin contamination must account for the presence of these other ricin-like molecules and should leverage the specificity provided by the numerous single nucleotide polymorphisms in the preproricin gene.

Impact of initial pathogen density on resistance and tolerance in a polymorphic disease resistance gene system in Arabidopsis thaliana

The evolution of natural enemy defense shapes evolutionary trajectories of natural populations. Although the intensity of selection imposed by enemies clearly varies among natural populations, little is known about the reaction norm of genotypes under a gradient of selective pressure. In this study, we measure the quantitative responses of disease symptoms and plant fitness to a gradient of infection, focusing on the gene-for-gene interaction between the Rpm1 resistance gene in Arabidopsis thaliana and the AvrRpm1 avirulence gene in the bacterial pathogen Pseudomonas syringae. Two complementary sets of plant material were used: resistant (R) and susceptible (S) isogenic lines and a set of six natural accessions, three of which are Rpm1 resistant (R) and three of which are rpm1 susceptible (S). Nine initial pathogen densities were applied to each plant line. Using isogenic lines allows any differences between R and S lines to be attributed directly to the Rpm1 gene, whereas using natural accessions allows the natural variation of resistance and tolerance over a gradient of infection dosages within R and S accessions to be described. For both sets of plant material, increased infection dosage results in more extensive disease symptoms, with a subsequent decrease in seed production. The severity of disease symptoms was reduced in R relative to S subgroups, and the presence of the Rpm1 allele led to an increase in plant fitness. Tolerance, defined as the ability to sustain infection without a reduction in fitness, was directly affected by Rpm1, providing a novel demonstration of an R gene affecting tolerance. Genetic variation for tolerance was also found within the S and R natural accessions, suggesting the potential for selection to act upon this important trait.

DNA polymorphism in the blast disease resistance gene Pita of the wild rice Oryza rufipogon and its related species

Intra-and interspecific DNA variations in the blast resistance gene Pita in wild rice (Oryza rufipogon), cultivated rice (O. sativa), and two other related wild rice species (O. meridionalis and O. officinalis) were analyzed to elucidate the nucleotide polymorphism maintenance mechanisms and evolution of Pita in these species. Nucleotide diversity at silent sites of O. rufipogon Pita was 0.0101, an intermediate value relative to other O. rufipogon nuclear genes. A dimorphic pattern of nucleotide polymorphism was detected in the O. rufipogon Pita region. Inoculation of the blast fungus Magnaporthe oryzae verified that the O. rufipogon Pita gene resides in a dimorphic sequence type. The resistance Pita allele had lower levels of variation than the susceptibility pita allele. A hypothesis of evolutionary relationships indicated that the amino acid mutation in the O. rufipogon Pita protein responsible for the difference between resistance and susceptibility occurred relatively recently. These results suggested that the resistance Pita originated from the susceptibility pita. Nucleotide diversity at replacement sites of the leucine-rich domain (LRD) of both the resistance and susceptibility O. rufipogon pita was low. In tests of neutrality, significantly negative values were detected for the LRD of O. rufipogon susceptibility pita. The low nucleotide diversity at replacement sites of the LRD of the susceptibility pita could be explained by purifying selection. Comparison of Pita between O. rufipogon and O. officinalis revealed an excess of nonconservative amino acid substitutions in the LRD, which could be related to the host-pathogen interaction.

Geographic variation in adaptation at the molecular level: a case study of plant immunity genes

Natural selection imposed by interacting species frequently varies among geographic locations and can lead to local adaptation, where alternative phenotypes are found in different populations. Little is known, however, about whether geographically variable selection acting on traits that mediate species interactions is consistent or strong enough to influence patterns of nucleotide variation at individual loci. To investigate this question, we examined patterns of nucleotide diversity and population structure at 16 plant innate immunity genes, with putative functions in defending plants against pathogens or herbivores, from six populations of teosinte (Zea mays ssp. parviglumis). Specifically, we tested whether patterns of population structure and within-population diversity at immunity genes differed from patterns found at nonimmunity (reference) loci and from neutral expectations derived from coalescent simulations of structured populations. For the majority of genes, we detected no strong evidence of geographically variable selection. However, in the wound-induced serine protease inhibitor (wip1), which inhibits the hydrolysis of dietary proteins in insect herbivores, one population showed unusually high levels of genetic differentiation, very low levels of nucleotide polymorphism, and was fixed for a novel replacement substitution in the active site of the protein. Taken together, these data suggest that wip1 experienced a recent selective sweep in one geographic region; this pattern may reflect local adaptation or an ongoing species-wide sweep. Overall, our results indicate that a signature of local adaptation at the molecular level may be uncommon-particularly for traits that are under complex genetic control.

Molecular evolution of the Pi-ta gene resistant to rice blast in wild rice (Oryza rufipogon)

Rice blast disease resistance to the fungal pathogen Magnaporthe grisea is triggered by a physical interaction between the protein products of the host R (resistance) gene, Pi-ta, and the pathogen Avr (avirulence) gene, AVR-pita. The genotype variation and resistant/susceptible phenotype at the Pi-ta locus of wild rice (Oryza rufipogon), the ancestor of cultivated rice (O. sativa), was surveyed in 36 locations worldwide to study the molecular evolution and functional adaptation of the Pi-ta gene. The low nucleotide polymorphism of the Pi-ta gene of O. rufipogon was similar to that of O. sativa, but greatly differed from what has been reported for other O. rufipogon genes. The haplotypes can be subdivided into two divergent haplogroups named H1 and H2. H1 is derived from H2, with nearly no variation and at a low frequency. H2 is common and is the ancestral form. The leucine-rich repeat (LRR) domain has a high pi(non)/pi(syn) ratio, and the low polymorphism of the Pi-ta gene might have primarily been caused by recurrent selective sweep and constraint by other putative physiological functions. Meanwhile, we provide data to show that the amino acid Ala-918 of H1 in the LRR domain has a close relationship with the resistant phenotype. H1 might have recently arisen during rice domestication and may be associated with the scenario of a blast pathogen-host shift from Italian millet to rice.

The evolution of virulence and pathogenicity in plant pathogen populations

The term virulence has a conflicting history among plant pathologists. Here we define virulence as the degree of damage caused to a host by parasite infection, assumed to be negatively correlated with host fitness, and pathogenicity the qualitative capacity of a parasite to infect and cause disease on a host. Selection may act on both virulence and pathogenicity, and their change in parasite populations can drive parasite evolution and host-parasite co-evolution. Extensive theoretical analyses of the factors that shape the evolution of pathogenicity and virulence have been reported in last three decades. Experimental work has not followed the path of theoretical analyses. Plant pathologists have shown greater interest in pathogenicity than in virulence, and our understanding of the molecular basis of pathogenicity has increased enormously. However, little is known regarding the molecular basis of virulence. It has been proposed that the mechanisms of recognition of parasites by hosts will have consequences for the evolution of pathogenicity, but much experimental work is still needed to test these hypotheses. Much theoretical work has been based on evidence from cellular plant pathogens. We review here the current experimental and observational evidence on which to test theoretical hypotheses or conjectures. We compare evidence from viruses and cellular pathogens, mostly fungi and oomycetes, which differ widely in genomic complexity and in parasitism. Data on the evolution of pathogenicity and virulence from viruses and fungi show important differences, and their comparison is necessary to establish the generality of hypotheses on pathogenicity and virulence evolution.

Genetic diversity in natural populations: a fundamental component of plant-microbe interactions

Genetic diversity for plant defense against microbial pathogens has been studied either by analyzing sequences of defense genes or by testing phenotypic responses to pathogens under experimental conditions. These two approaches give different but complementary information but, till date, only rare attempts at their integration have been made. Here we discuss the advances made, because of the two approaches, in understanding plant-pathogen coevolution and propose ways of integrating the two.

Investigation of the demographic and selective forces shaping the nucleotide diversity of genes involved in nod factor signaling in Medicago truncatula

Symbiotic nitrogen-fixing rhizobia are able to trigger root deformation in their Fabaceae host plants, allowing their intracellular accommodation. They do so by delivering molecules called Nod factors. We analyzed the patterns of nucleotide polymorphism of five genes controlling early Nod factor perception and signaling in the Fabaceae Medicago truncatula to understand the selective forces shaping the evolution of these genes. We used 30 M. truncatula genotypes sampled in a genetically homogeneous region of the species distribution range. We first sequenced 24 independent loci and detected a genomewide departure from the hypothesis of neutrality and demographic equilibrium that suggests a population expansion. These data were used to estimate parameters of a simple demographic model incorporating population expansion. The selective neutrality of genes controlling Nod factor perception was then examined using a combination of two complementary neutrality tests, Tajima's D and Fay and Wu's standardized H. The joint distribution of D and H expected under neutrality was obtained under the fitted population expansion model. Only the gene DMI1, which is expected to regulate the downstream signal, shows a pattern consistent with a putative selective event. In contrast, the receptor-encoding genes NFP and NORK show no significant signatures of selection. Among the genes that we analyzed, only DMI1 should be viewed as a candidate for adaptation in the recent history of M. truncatula.

The Evolution of Resistance and Tolerance to Herbivores

Tolerance and resistance are two different plant defense strategies against herbivores. Empirical evidence in natural populations reveals that individual plants allocate resources simultaneously to both strategies, thus plants exhibit a mixed pattern of defense. In this review we examine the conditions that promote the evolutionary stability of mixed defense strategies in the light of available empirical and theoretical evidence. Given that plant tolerance and resistance are heritable and subject to environmentally dependent selection and genetic constraints, the joint evolution of tolerance and resistance is analyzed, with consideration of multiple species interactions and the plant mating system. The existence of mixed defense strategies in plants makes it necessary to re-explore the coevolutionary process between plants and herbivores, which centered historically on resistance as the only defensive mechanism. In addition, we recognize briefly the potential use of plant tolerance for pest management. Finally, we highlight unresolved issues for future development in this field of evolutionary ecology.

Stability of genetic polymorphism in host-parasite interactions

Allelic diversity is common at host loci involved in parasite recognition, such as the major histocompatibility complex in vertebrates or gene-for-gene relationships in plants, and in corresponding loci encoding antigenic molecules in parasites. Diverse factors have been proposed in models to account for genetic polymorphism in host-parasite recognition. Here, a simple but general theory of host-parasite coevolution is developed. Coevolution implies the existence of indirect frequency-dependent selection (FDS), because natural selection on the host depends on the frequency of a parasite gene, and vice versa. It is shown that polymorphism can be maintained in both organisms only if there is negative, direct FDS, such that the strength of natural selection for the host resistance allele, the parasite virulence allele or both declines with increasing frequency of that allele itself. This condition may be fulfilled if the parasite has more than one generation in the same host individual, a feature which is common to most diseases. It is argued that the general theory encompasses almost all factors previously proposed to account for polymorphism at corresponding host and parasite loci, including those controlling gene-for-gene interactions.

Natural variation in the Pto disease resistance gene within species of wild tomato (Lycopersicon). II. Population genetics of Pto

Disease resistance to the bacterial pathogen Pseudomonas syringae pv. tomato (Pst) in the host species Lycopersicon esculentum, the cultivated tomato, and the closely related L. pimpinellifolium is triggered by the physical interaction between the protein products of the host resistance (R) gene Pto and the pathogen avirulence genes AvrPto and AvrPtoB. Sequence variation at the Pto locus was surveyed in natural populations of seven species of Lycopersicon to test hypotheses of host-parasite coevolution and functional adaptation of the Pto gene. Pto shows significantly higher nonsynonymous polymorphism than 14 other non-R-gene loci in the same samples of Lycopersicon species, while showing no difference in synonymous polymorphism, suggesting that the maintenance of amino acid polymorphism at this locus is mediated by pathogen selection. Also, a larger proportion of ancestral variation is maintained at Pto as compared to these non-R-gene loci. The frequency spectrum of amino acid polymorphisms known to negatively affect Pto function is skewed toward low frequency compared to amino acid polymorphisms that do not affect function or silent polymorphisms. Therefore, the evolution of Pto appears to be influenced by a mixture of both purifying and balancing selection.

Licensed to kill: the lifestyle of a necrotrophic plant pathogen

Necrotrophic plant pathogens have received an increasing amount of attention over the past decade. Initially considered to invade their hosts in a rather unsophisticated manner, necrotrophs are now known to use subtle mechanisms to subdue host plants. The gray mould pathogen Botrytis cinerea is one of the most comprehensively studied necrotrophic fungal plant pathogens. The genome sequences of two strains have been determined. Targeted mutagenesis studies are unraveling the roles played in the infection process by a variety of B. cinerea genes that are required for penetration, host cell killing, plant tissue decomposition or signaling. Our increasing understanding of the tools used by a necrotrophic fungal pathogen to invade plants will be instrumental to designing rational strategies for disease control.

A catalogue of the effector secretome of plant pathogenic oomycetes

The oomycetes form a phylogenetically distinct group of eukaryotic microorganisms that includes some of the most notorious pathogens of plants. Oomycetes accomplish parasitic colonization of plants by modulating host cell defenses through an array of disease effector proteins. The biology of effectors is poorly understood but tremendous progress has been made in recent years. This review classifies and catalogues the effector secretome of oomycetes. Two classes of effectors target distinct sites in the host plant: Apoplastic effectors are secreted into the plant extracellular space, and cytoplasmic effectors are translocated inside the plant cell, where they target different subcellular compartments. Considering that five species are undergoing genome sequencing and annotation, we are rapidly moving toward genome-wide catalogues of oomycete effectors. Already, it is evident that the effector secretome of pathogenic oomycetes is more complex than expected, with perhaps several hundred proteins dedicated to manipulating host cell structure and function.

Genetic diversity and the evolutionary history of plant immunity genes in two species of Zea

Plant pathogenesis-related genes (PR genes) code for enzymes, enzyme inhibitors, and other peptides that confer resistance to pathogens and herbivores. Although several PR genes have been the subject of molecular population genetic analyses, a general understanding of their long-term evolutionary dynamics remains incomplete. Here we analyze sequence data from 17 PR genes from two closely related teosinte species of central Mexico. In addition to testing whether patterns of diversity at individual loci depart from expectations under a neutral model, we compared patterns of diversity at defense genes, as a class, to nondefense genes. In Zea diploperennis, the majority of defense genes have patterns of diversity consistent with neutral expectations while at least two genes showed evidence of recent positive selection consistent with arms-race models of antagonistic coevolution. In Zea mays ssp. parviglumis, by contrast, analyses of both defense and nondefense genes revealed strong and consistent departures from the neutral model, suggestive of nonequilibrium population dynamics or population structure. Nevertheless, we found a significant excess of replacement polymorphism in defense genes compared to nondefense genes. Although we cannot exclude relaxed selective constraint as an explanation, our results are consistent with temporally variable (transient or episodic) selection or geographically variable selection acting on parviglumis defense genes. The different patterns of diversity found in the two Zea species may be explained by parviglumis' greater distribution and population structure together with geographic variation in selection.

Rearrangements in the Cf-9 disease resistance gene cluster of wild tomato have resulted in three genes that mediate Avr9 responsiveness

Cf resistance genes in tomato confer resistance to the fungal leaf pathogen Cladosporium fulvum. Both the well-characterized resistance gene Cf-9 and the related 9DC gene confer resistance to strains of C. fulvum that secrete the Avr9 protein and originate from the wild tomato species Lycopersicon pimpinellifolium. We show that 9DC and Cf-9 are allelic, and we have isolated and sequenced the complete 9DC cluster of L. pimpinellifolium LA1301. This 9DC cluster harbors five full-length Cf homologs, including orthologs of the most distal homologs of the Cf-9 cluster and three central 9DC genes. Two 9DC genes (9DC1 and 9DC2) have an identical coding sequence, whereas 9DC3 differs at its 3' terminus. From a detailed comparison of the 9DC and Cf-9 clusters, we conclude that the Cf-9 and Hcr9-9D genes from the Cf-9 cluster are ancestral to the first 9DC gene and that the three 9DC genes were generated by subsequent intra- and intergenic unequal recombination events. Thus, the 9DC cluster has undergone substantial rearrangements in the central region, but not at the ends. Using transient transformation assays, we show that all three 9DC genes confer Avr9 responsiveness, but that 9DC2 is likely the main determinant of Avr9 recognition in LA1301.

Population genetic evidence for rapid changes in intraspecific diversity and allelic cycling of a specialist defense gene in Zea

Two patterns of plant defense gene evolution are emerging from molecular population genetic surveys. One is that specialist defenses experience stronger selection than generalist defenses. The second is that specialist defenses are more likely to be subject to balancing selection, i.e., evolve in a manner consistent with balanced-polymorphism or trench-warfare models of host-parasite coevolution. Because most of the data of specialist defenses come from Arabidopsis thaliana, we examined the genetic diversity and evolutionary history of three defense genes in two outcrossing species, the autotetraploid Zea perennis and its most closely related extant relative the diploid Z. diploperennis. Intraspecific diversity at two generalist defenses, the protease inhibitors wip1 and mpi, were consistent with a neutral model. Like previously studied genes in these taxa, wip1 and mpi harbored similar levels of diversity in Z. diploperennis and Z. perennis. In contrast, the specialist defense hm2 showed strong although distinctly different departures from a neutral model in the two species. Z. diploperennis appears to have experienced a strong and recent selective sweep. Using a rejection-sampling coalescent method, we estimate the strength of selection on Z. diploperennis hm2 to be approximately 3.0%, which is approximately equal to the strength of selection on tb1 during maize domestication. Z. perennis hm2 harbors three highly diverged alleles, two of which are found at high frequency. The distinctly different patterns of diversity may be due to differences in the phase of host-parasite coevolutionary cycles, although higher hm2 diversity in Z. perennis may also reflect reduced efficacy of selection in the autotetraploid relative to its diploid relative.

Comparative evolutionary histories of chitinase genes in the Genus zea and Family poaceae

Patterns of DNA sequence diversity vary widely among genes encoding proteins that protect plants against pathogens and herbivores. Comparative studies may help determine whether these differences are due to the strength of selection acting on different types of defense, in different evolutionary lineages, or both. I analyzed sequence diversity at three chitinases, a well-studied component of defense, in two species of Zea and several Poaceae taxa. Although the Zea species are closely related and these genes code for proteins with similar biochemical function, patterns of diversity vary widely within and among species. Intraspecific diversity at chiB, chiI, and Z. mays ssp. parviglumis chiA are consistent with a neutral-equilibrium model whereas chiA had no segregating sites within Z. diploperennis--consistent with a recent and strong selective sweep. Codons identified as having diverged among Poaceae taxa in response to positive selection were significantly overrepresented among targets of selection in Arabis, suggesting common responses to selection in distantly related plant taxa. Divergence of the recent duplicates chiA and chiB is consistent with positive selection but relaxed constraint cannot be rejected. Weak evidence for adaptive divergence of these duplicated downstream components of defense contrasts with strong evidence for adaptive divergence of genes involved in pathogen recognition.

Selection on Glycine beta-1,3-endoglucanase genes differentially inhibited by a Phytophthora glucanase inhibitor protein

Plant endo-beta-1,3-glucanases (EGases) degrade the cell wall polysaccharides of attacking pathogens and release elicitors of additional plant defenses. Isozymes EGaseA and EGaseB of soybean differ in susceptibility to a glucanase inhibitor protein (GIP1) produced by Phytophthora sojae, a major soybean pathogen. EGaseA, the major elicitor-releasing isozyme, is a high-affinity ligand for GIP1, which completely inhibits it, whereas EGaseB is unaffected by GIP1. We tested for departures from neutral evolution on the basis of partial sequences of EGaseA and EGaseB from 20 widespread accessions of Glycine soja (the wild progenitor of soybean), from 4 other Glycine species, and across dicotyledonous plants. G. soja exhibited little intraspecific variation at either locus. Phylogeny-based codon evolution models detected strong evidence of positive selection on Glycine EGaseA and weaker evidence for selection on dicot EGases and Glycine EGaseB. Positively selected peptide sites were identified and located on a structural model of EGase bound to GIP1. Positively selected sites and highly variable sites were found disproportionately within 4.5 angstroms of bound GIP1. Low variation within G. soja EGases, coupled with positive selection in both Glycine and dicot lineages and the proximity of rapidly evolving sites to GIP1, suggests an arms race involving repeated adaptation to pathogen attack and inhibition.

Evolution of the innate immune system: the worm perspective

Simple model organisms that are amenable to comprehensive experimental analysis can be used to elucidate the molecular genetic architecture of complex traits. They can thereby enhance our understanding of these traits in other organisms, including humans. Here, we describe the use of the nematode Caenorhabditis elegans as a tractable model system to study innate immunity. We detail our current understanding of the worm's immune system, which seems to be characterized by four main signaling cascades: a p38 mitogen-activated protein kinase, a transforming growth factor-beta-like, a programed cell death, and an insulin-like receptor pathway. Many details, especially regarding pathogen recognition and immune effectors, are only poorly characterized and clearly warrant further investigation. We additionally speculate on the evolution of the C. elegans immune system, taking into special consideration the relationship between immunity, stress responses and digestion, the diversification of the different parts of the immune system in response to multiple and/or coevolving pathogens, and the trade-off between immunity and host life history traits. Using C. elegans to address these different facets of host-pathogen interactions provides a fresh perspective on our understanding of the structure and complexity of innate immune systems in animals and plants.