Dynamics of disease resistance polymorphism at the Rpm1 locus of Arabidopsis

Ancient diversity of splicing motifs and protein surfaces in the wild emmer wheat (Triticum dicoccoides) LR10 coiled coil (CC) and leucine-rich repeat (LRR) domains

In this study, we explore the diversity and its distribution along the wheat leaf rust resistance protein LR10 three-dimensional structure. Lr10 is a leaf rust resistance gene encoding a coiled coil-nucleotide-binding site-leucine-rich repeat (CC-NBS-LRR) class of protein. Lr10 was cloned and sequenced from 58 accessions representing diverse habitats of wild emmer wheat in Israel. Nucleotide diversity was very high relative to other wild emmer wheat genes (π= 0.029). The CC domain was found to be the most diverse domain and subject to positive selection. Superimposition of the diversity on the CC three-dimensional structure showed that some of the variable and positively selected residues were solvent exposed and may interact with other proteins. The LRR domain was relatively conserved, but showed a hotspot of amino acid variation between two haplotypes in the ninth repeat. This repeat was longer than the other LRRs, and three-dimensional modelling suggested that an extensive α helix structure was formed in this region. The two haplotypes also differed in splicing regulation motifs. In genotypes with one haplotype, an intron was alternatively spliced in this region, whereas, in genotypes with the other haplotype, this intron did not splice at all. The two haplotypes are proposed to be ancient and maintained by balancing selection.

Genetic structure, linkage disequilibrium and signature of selection in Sorghum: lessons from physically anchored DArT markers

Population structure, extent of linkage disequilibrium (LD) as well as signatures of selection were investigated in sorghum using a core sample representative of worldwide diversity. A total of 177 accessions were genotyped with 1122 informative physically anchored DArT markers. The properties of DArTs to describe sorghum genetic structure were compared to those of SSRs and of previously published RFLP markers. Model-based (STRUCTURE software) and Neighbor-Joining diversity analyses led to the identification of 6 groups and confirmed previous evolutionary hypotheses. Results were globally consistent between the different marker systems. However, DArTs appeared more robust in terms of data resolution and bayesian group assignment. Whole genome linkage disequilibrium as measured by mean r(2) decreased from 0.18 (between 0 to 10 kb) to 0.03 (between 100 kb to 1 Mb), stabilizing at 0.03 after 1 Mb. Effects on LD estimations of sample size and genetic structure were tested using i. random sampling, ii. the Maximum Length SubTree algorithm (MLST), and iii. structure groups. Optimizing population composition by the MLST reduced the biases in small samples and seemed to be an efficient way of selecting samples to make the best use of LD as a genome mapping approach in structured populations. These results also suggested that more than 100,000 markers may be required to perform genome-wide association studies in collections covering worldwide sorghum diversity. Analysis of DArT markers differentiation between the identified genetic groups pointed out outlier loci potentially linked to genes controlling traits of interest, including disease resistance genes for which evidence of selection had already been reported. In addition, evidence of selection near a homologous locus of FAR1 concurred with sorghum phenotypic diversity for sensitivity to photoperiod.

Flagellin perception varies quantitatively in Arabidopsis thaliana and its relatives

Much is known about the evolution of plant immunity components directed against specific pathogen strains: They show pervasive functional variation and have the potential to coevolve with pathogen populations. However, plants are effectively protected against most microbes by generalist immunity components that detect conserved pathogen-associated molecular patterns (PAMPs) and control the onset of PAMP-triggered immunity. In Arabidopsis thaliana, the receptor kinase flagellin sensing 2 (FLS2) confers recognition of bacterial flagellin (flg22) and activates a manifold defense response. To decipher the evolution of this system, we performed functional assays across a large set of A. thaliana genotypes and Brassicaceae relatives. We reveal extensive variation in flg22 perception, most of which results from changes in protein abundance. The observed variation correlates with both the severity of elicited defense responses and bacterial proliferation. We analyzed nucleotide variation segregating at FLS2 in A. thaliana and detected a pattern of variation suggestive of the rapid fixation of a novel adaptive allele. However, our study also shows that evolution at the receptor locus alone does not explain the evolution of flagellin perception; instead, components common to pathways downstream of PAMP perception likely contribute to the observed quantitative variation. Within and among close relatives, PAMP perception evolves quantitatively, which contrasts with the changes in recognition typically associated with the evolution of R genes.

Parallel Loss-of-Function at the RPM1 Bacterial Resistance Locus in Arabidopsis thaliana

Dimorphism at the Resistance to Pseudomonas syringae pv. maculicola 1 (RPM1) locus is well documented in natural populations of Arabidopsis thaliana and has been portrayed as a long-term balanced polymorphism. The haplotype from resistant plants contains the RPM1 gene, which enables these plants to recognize at least two structurally unrelated bacterial effector proteins (AvrB and AvrRpm1) from bacterial crop pathogens. A complete deletion of the RPM1 coding sequence has been interpreted as a single event resulting in susceptibility in these individuals. Consequently, the ability to revert to resistance or for alternative R-gene specificities to evolve at this locus has also been lost in these individuals. Our survey of variation at the RPM1 locus in a large species-wide sample of A. thaliana has revealed four new loss-of-function alleles that contain most of the intervening sequence of the RPM1 open reading frame. Multiple loss-of-function alleles may have originated due to the reported intrinsic cost to plants expressing the RPM1 protein. The frequency and geographic distribution of rpm1 alleles observed in our survey indicate the parallel origin and maintenance of these loss-of-function mutations and reveal a more complex history of natural selection at this locus than previously thought.

Molecular evolution of a family of resistance gene analogs of nucleotide-binding site sequences in Solanum lycopersicum

Nucleotide-binding site-leucine-rich repeats (NBS-LRR) gene families are one of the major plant resistance genes. Genomic NBS evolution was studied in many plant species for diverse arrays of NBS gene families. In this study, we focused on one family of NBS sequences in an attempt to understand how closely related NBS sequences evolved in the light of selection in domesticated plant species. A phylogenetic analysis revealed five major clades (A-E) and five subclades (A1-A5) within clade A of cloned NBS sequences. Positive selection was only detected in newly evolved NBS lineages in subclades of clade A. Positively selected codon sites were found among NBS sequences of clade A. A sliding-window analysis revealed that regions with Ka/Ks ratios of >1 were in the inter-motifs when paired clades were compared, but regions with Ka/Ks ratios of >1 were found across NBS sequences when subclades of clade A were compared. Our results based on a family of closely related NBS sequences showed that positive selection was first exerted on specific lineages across all NBS sequences after selective constraints. Subsequently, sequences with mutations in commonly conserved motifs were scrutinized by purifying selection. In the long term, conserved high frequency alleles in commonly conserved motifs and changes in inter-motifs were maintained in the investigated family of NBS sequences. Moreover, codons identified to be under positive selection in the inter-motifs were mainly located in regions involved in functions of ATP binding or hydrolysis.

Population genetics of duplicated alternatively spliced exons of the Dscam gene in Daphnia and Drosophila

In insects and crustaceans, the Down syndrome cell adhesion molecule (Dscam) occurs in many different isoforms. These are produced by mutually exclusive alternative splicing of dozens of tandem duplicated exons coding for parts or whole immunoglobulin (Ig) domains of the Dscam protein. This diversity plays a role in the development of the nervous system and also in the immune system. Structural analysis of the protein suggested candidate epitopes where binding to pathogens could occur. These epitopes are coded by regions of the duplicated exons and are therefore diverse within individuals. Here we apply molecular population genetics and molecular evolution analyses using Daphnia magna and several Drosophila species to investigate the potential role of natural selection in the divergence between orthologs of these duplicated exons among species, as well as between paralogous exons within species. We found no evidence for a role of positive selection in the divergence of these paralogous exons. However, the power of this test was low, and the fact that no signs of gene conversion between paralogous exons were found suggests that paralog diversity may nonetheless be maintained by selection. The analysis of orthologous exons in Drosophila and in Daphnia revealed an excess of non-synonymous polymorphisms in the epitopes putatively involved in pathogen binding. This may be a sign of balancing selection. Indeed, in Dr. melanogaster the same derived non-synonymous alleles segregate in several populations around the world. Yet other hallmarks of balancing selection were not found. Hence, we cannot rule out that the excess of non-synonymous polymorphisms is caused by segregating slightly deleterious alleles, thus potentially indicating reduced selective constraints in the putative pathogen binding epitopes of Dscam.

Genome-wide patterns of genetic variation in sweet and grain sorghum (Sorghum bicolor)

Background

Sorghum (Sorghum bicolor) is globally produced as a source of food, feed, fiber and fuel. Grain and sweet sorghums differ in a number of important traits, including stem sugar and juice accumulation, plant height as well as grain and biomass production. The first whole genome sequence of a grain sorghum is available, but additional genome sequences are required to study genome-wide and intraspecific variation for dissecting the genetic basis of these important traits and for tailor-designed breeding of this important C₄ crop.

Results

We resequenced two sweet and one grain sorghum inbred lines, and identified a set of nearly 1,500 genes differentiating sweet and grain sorghum. These genes fall into ten major metabolic pathways involved in sugar and starch metabolisms, lignin and coumarin biosynthesis, nucleic acid metabolism, stress responses and DNA damage repair. In addition, we uncovered 1,057,018 SNPs, 99,948 indels of 1 to 10 bp in length and 16,487 presence/absence variations as well as 17,111 copy number variations. The majority of the large-effect SNPs, indels and presence/absence variations resided in the genes containing leucine rich repeats, PPR repeats and disease resistance R genes possessing diverse biological functions or under diversifying selection, but were absent in genes that are essential for life.

Conclusions

This is a first report of the identification of genome-wide patterns of genetic variation in sorghum. High-density SNP and indel markers reported here will be a valuable resource for future gene-phenotype studies and the molecular breeding of this important crop and related species.

Spatial heterogeneity, frequency-dependent selection and polymorphism in host-parasite interactions

BACKGROUND

Genomic and pathology analysis has revealed enormous diversity in genes involved in disease, including those encoding host resistance and parasite effectors (also known in plant pathology as avirulence genes). It has been proposed that such variation may persist when an organism exists in a spatially structured metapopulation, following the geographic mosaic of coevolution. Here, we study gene-for-gene relationships governing the outcome of plant-parasite interactions in a spatially structured system and, in particular, investigate the population genetic processes which maintain balanced polymorphism in both species.

RESULTS

Following previous theory on the effect of heterogeneous environments on maintenance of polymorphism, we analysed a model with two demes in which the demes have different environments and are coupled by gene flow. Environmental variation is manifested by different coefficients of natural selection, the costs to the host of resistance and to the parasite of virulence, the cost to the host of being diseased and the cost to an avirulent parasite of unsuccessfully attacking a resistant host. We show that migration generates negative direct frequency-dependent selection, a condition for maintenance of stable polymorphism in each deme. Balanced polymorphism occurs preferentially if there is heterogeneity for costs of resistance and virulence alleles among populations and to a lesser extent if there is variation in the cost to the host of being diseased. We show that the four fitness costs control the natural frequency of oscillation of host resistance and parasite avirulence alleles. If demes have different costs, their frequencies of oscillation differ and when coupled by gene flow, there is amplitude death of the oscillations in each deme. Numerical simulations show that for a multiple deme island model, costs of resistance and virulence need not to be present in each deme for stable polymorphism to occur.

CONCLUSIONS

Our theoretical results confirm the importance of empirical studies for measuring the environmental heterogeneity for genetic costs of resistance and virulence alleles. We suggest that such studies should be developed to investigate the generality of this mechanism for the long-term maintenance of genetic diversity at host and parasite genes.

Successive increases in the resistance of Drosophila to viral infection through a transposon insertion followed by a Duplication

To understand the molecular basis of how hosts evolve resistance to their parasites, we have investigated the genes that cause variation in the susceptibility of Drosophila melanogaster to viral infection. Using a host-specific pathogen of D. melanogaster called the sigma virus (Rhabdoviridae), we mapped a major-effect polymorphism to a region containing two paralogous genes called CHKov1 and CHKov2. In a panel of inbred fly lines, we found that a transposable element insertion in the protein coding sequence of CHKov1 is associated with increased resistance to infection. Previous research has shown that this insertion results in a truncated messenger RNA that encodes a far shorter protein than the susceptible allele. This resistant allele has rapidly increased in frequency under directional selection and is now the commonest form of the gene in natural populations. Using genetic mapping and site-specific recombination, we identified a third genotype with considerably greater resistance that is currently rare in the wild. In these flies there have been two duplications, resulting in three copies of both the truncated allele of CHKov1 and CHKov2 (one of which is also truncated). Remarkably, the truncated allele of CHKov1 has previously been found to confer resistance to organophosphate insecticides. As estimates of the age of this allele predate the use of insecticides, it is likely that this allele initially functioned as a defence against viruses and fortuitously "pre-adapted" flies to insecticides. These results demonstrate that strong selection by parasites for increased host resistance can result in major genetic changes and rapid shifts in allele frequencies; and, contrary to the prevailing view that resistance to pathogens can be a costly trait to evolve, the pleiotropic effects of these changes can have unexpected benefits.

Molecular evolution under increasing transposable element burden in Drosophila: a speed limit on the evolutionary arms race

BACKGROUND

Genome architecture is profoundly influenced by transposable elements (TEs), and natural selection against their harmful effects is a critical factor limiting their spread. Genome defense by the piRNA silencing pathway also plays a crucial role in limiting TE proliferation. How these two forces jointly determine TE abundance is not well understood. To shed light on the nature of factors that predict TE success, we test three distinct hypotheses in the Drosophila genus. First, we determine whether TE abundance and relaxed genome-wide purifying selection on protein sequences are positively correlated. This serves to test the hypothesis that variation in TE abundance in the Drosophila genus can be explained by the strength of natural selection, relative to drift, acting in parallel against mildly deleterious non-synonymous mutations. Second, we test whether increasing TE abundance is correlated with an increased rate of amino-acid evolution in genes encoding the piRNA machinery, as might be predicted by an evolutionary arms race model. Third, we test whether increasing TE abundance is correlated with greater codon bias in genes of the piRNA machinery. This is predicted if increasing TE abundance selects for increased efficiency in the machinery of genome defense.

RESULTS

Surprisingly, we find neither of the first two hypotheses to be true. Specifically, we found that genome-wide levels of purifying selection, measured by the ratio of non-synonymous to synonymous substitution rates (ω), were greater in species with greater TE abundance. In addition, species with greater TE abundance have greater levels of purifying selection in the piRNA machinery. In contrast, it appears that increasing TE abundance has primarily driven adaptation in the piRNA machinery by increasing codon bias.

CONCLUSIONS

These results indicate that within the Drosophila genus, a historically reduced strength of selection relative to drift is unlikely to explain patterns of increased TE success across species. Other factors, such as ecological exposure, are likely to contribute to variation in TE abundances within species. Furthermore, constraints on the piRNA machinery may temper the evolutionary arms race that would drive increasing rates of evolution at the amino acid level. In the face of these constraints, selection may act primarily by improving the translational efficiency of the machinery of genome defense through efficient codon usage.

Molecular adaptation in flowering and symbiotic recognition pathways: insights from patterns of polymorphism in the legume Medicago truncatula

Background

We studied patterns of molecular adaptation in the wild Mediterranean legume Medicago truncatula. We focused on two phenotypic traits that are not functionally linked: flowering time and perception of symbiotic microbes. Phenology is an important fitness component, especially for annual plants, and many instances of molecular adaptation have been reported for genes involved in flowering pathways. While perception of symbiotic microbes is also integral to adaptation in many plant species, very few reports of molecular adaptation exist for symbiotic genes. Here we used data from 57 individuals and 53 gene fragments to quantify the overall strength of both positive and purifying selection in M. truncatula and asked if footprints of positive selection can be detected at key genes of rhizobia recognition pathways.

Results

We examined nucleotide variation among 57 accessions from natural populations in 53 gene fragments: 5 genes involved in nitrogen-fixing bacteria recognition, 11 genes involved in flowering, and 37 genes used as control loci. We detected 1757 polymorphic sites yielding an average nucleotide diversity (pi) of 0.003 per site. Non-synonymous variation is under sizable purifying selection with 90% of amino-acid changing mutations being strongly selected against. Accessions were structured in two groups consistent with geographical origins. Each of these two groups harboured an excess of rare alleles, relative to expectations of a constant-sized population, suggesting recent population expansion. Using coalescent simulations and an approximate Bayesian computation framework we detected several instances of genes departing from selective neutrality within each group and showed that the polymorphism of two nodulation and four flowering genes has probably been shaped by recent positive selection.

Conclusion

We quantify the intensity of purifying selection in the M. truncatula genome and show that putative footprints of natural selection can be detected at different time scales in both flowering and symbiotic pathways.

Source verification of mis-identified Arabidopsis thaliana accessions

A major strength of Arabidopsis thaliana as a model lies in the availability of a large number of naturally occurring inbred lines. Recent studies of A. thaliana population structure, using thousands of accessions from stock center and natural collections, have revealed a robust pattern of isolation by distance at several spatial scales, such that genetically identical individuals are generally found close to each other. However, some individual accessions deviate from this pattern. While some of these may be the products of rare long-distance dispersal events, many deviations may be the result of mis-identification, in the sense that the data regarding location of origin data are incorrect. Here, we aim to identify such discrepancies. Of the 5965 accessions examined, we conclude that 286 deserve special attention as being potentially mis-identified. We describe these suspicious accessions and their possible origins, and advise caution with regard to their use in experiments in which accurate information on geographic origin is important. Finally, we discuss possibilities for maintaining the integrity of stock lines.

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.

Unusual signatures of highly adaptable R-loci in closely-related Arabidopsis species

Plant resistance genes (R-genes) evolve rapidly in response to changing environments. What are the most remarkable signatures of fast adaptive genes, besides the commonly revealed rapid divergence and high non-synonymous substitution rate? Here we investigated these changes in five R-loci following recent differentiation between Arabidopsis thaliana and Arabidopsis lyrata. Extreme differences in evolutionary rates were observed: e.g., an overall 5.46-9.83-fold different nucleotide diversity at two R-loci between species, ten-fold higher non-synonymous substitution rates within one species versus the other, significantly different Ka/Ks ratios between species for the same R-gene, and high interspecific divergence at one R-locus. Particularly, we observed an elevated level of trans-specific polymorphism at one R-locus and a differentially maintained presence/absence polymorphism at another. The high frequency of ancestral polymorphisms amongst R-genes suggests that the persistence of some functional variation is an important evolutionary mechanism shaping genetic variation in R-genes, while the variation of presence/absence polymorphisms provides a potential mechanism for malleable activation of adaptive resistance response pathways. The distinct patterns among R-genes suggest that the same R-gene ortholog can be quickly shaped by different evolutionary processes, e.g., purifying selection in one species but positive selection in a closely-related species.

A multifaceted genomics approach allows the isolation of the rice Pia-blast resistance gene consisting of two adjacent NBS-LRR protein genes

The Oryza sativa (rice) resistance gene Pia confers resistance to the blast fungus Magnaporthe oryzae carrying the AVR-Pia avirulence gene. To clone Pia, we employed a multifaceted genomics approach. First, we selected 12 R-gene analog (RGA) genes encoding nucleotide binding site-leucine rich repeats (NBS-LRRs) proteins from a region on chromosome 11 that shows linkage to Pia. By using seven rice accessions, we examined the association between Pia phenotypes and DNA polymorphisms in the 10 genes, which revealed three genes (Os11gRGA3-Os11gRGA5) exhibiting a perfect association with the Pia phenotypes. We also screened ethyl methane sulfonate (EMS)-treated mutant lines of the rice cultivar 'Sasanishiki' harboring Pia, and isolated two mutants that lost the Pia phenotype. DNA sequencing of Os11gRGA3-Os11gRGA5 from the two mutant lines identified independent mutations of major effects in Os11gRGA4. The wild-type 'Sasanishiki' allele of Os11gRGA4 (SasRGA4) complemented Pia function in both mutants, suggesting that SasRGA4 is necessary for Pia function. However, when the rice cultivar 'Himenomochi' lacking Pia was transfected with SasRGA4, the Pia phenotype was not recovered. An additional complementation study revealed that the two NBS-LRR-type R genes, SasRGA4 and SasRGA5, that are located next to each other and oriented in the opposite direction are necessary for Pia function. A population genetics analysis of SasRGA4 and SasRGA5 suggests that the two genes are under long-term balancing selection.

Genomic and resistance gene homolog diversity of the dominant tallgrass prairie species across the U.S. Great Plains precipitation gradient

BACKGROUND

Environmental variables such as moisture availability are often important in determining species prevalence and intraspecific diversity. The population genetic structure of dominant plant species in response to a cline of these variables has rarely been addressed. We evaluated the spatial genetic structure and diversity of Andropogon gerardii populations across the U.S. Great Plains precipitation gradient, ranging from approximately 48 cm/year to 105 cm/year.

METHODOLOGY/PRINCIPAL FINDINGS

Genomic diversity was evaluated with AFLP markers and diversity of a disease resistance gene homolog was evaluated by PCR-amplification and digestion with restriction enzymes. We determined the degree of spatial genetic structure using Mantel tests. Genomic and resistance gene homolog diversity were evaluated across prairies using Shannon's index and by averaging haplotype dissimilarity. Trends in diversity across prairies were determined using linear regression of diversity on average precipitation for each prairie. We identified significant spatial genetic structure, with genomic similarity decreasing as a function of distance between samples. However, our data indicated that genome-wide diversity did not vary consistently across the precipitation gradient. In contrast, we found that disease resistance gene homolog diversity was positively correlated with precipitation.

SIGNIFICANCE

Prairie remnants differ in the genetic resources they maintain. Selection and evolution in this disease resistance homolog is environmentally dependent. Overall, we found that, though this environmental gradient may not predict genomic diversity, individual traits such as disease resistance genes may vary significantly.

Faster evolution of Z-linked duplicate genes in chicken

It has been shown that duplicate genes on the X chromosome evolve much faster than duplicate genes on autosomes in Drosophila melanogaster. However, whether this phenomenon is general and can be applied to other species is not known. Here we examined this issue in chicken that have heterogametic females (females have ZW sex chromosome). We compared sequence divergence of duplicate genes on the Z chromosome with those on autosomes. We found that duplications on the Z chromosome indeed evolved faster than those on autosomes and show distinct patterns of molecular evolution from autosomal duplications. Examination of the expression of duplicate genes revealed an enrichment of duplications on the Z chromosome having male-biased expression and an enrichment of duplications on the autosomes having female-biased expression. These results suggest an evolutionary trend of the recruitment of duplicate genes towards reproduction-specific function. The faster evolution of duplications on Z than on the autosomes is most likely contributed by the selective forces driving the fixation of adaptive mutations on Z. Therefore, the common phenomena observed in both flies and chicken suggest that duplicate genes on sex chromosomes have distinct dynamics and are more influenced by natural selection than autosomal duplications, regardless of the kind of sex determination systems.

In vitro tests of natural allelic variation of innate immune genes (avian β-defensins) reveal functional differences in microbial inhibition

Allelic variation in immune genes might result from, and contribute to, host-pathogen evolution. Functional allelic variation in the innate immune system has received little attention. Here, we investigate whether naturally occurring allelic variation within the avian innate immune system (β-defensins) is associated with variation in antimicrobial activity. We tested differences in in vitro antimicrobial properties of the synthesized products of two alleles of avian β-defensin 7, both of which occur at high frequency in natural populations of the great tit (Parus major). Only one allele strongly inhibited the growth of the gram-positive bacterium Staphylococcus aureus, but both alleles strongly inhibited growth of the gram-negative bacterium Escherechia coli. Our data demonstrate functional allelic variation in natural defensin genes, and we discuss how differences in efficacy against microbial species might contribute to maintaining this variation.

Plant-parasite coevolution: bridging the gap between genetics and ecology

We review current ideas about coevolution of plants and parasites, particularly processes that generate genetic diversity. Frequencies of host resistance and parasite virulence alleles that interact in gene-for-gene (GFG) relationships coevolve in the familiar boom-and-bust cycle, in which resistance is selected when virulence is rare, and virulence is selected when resistance is common. The cycle can result in stable polymorphism when diverse ecological and epidemiological factors cause negative direct frequency-dependent selection (ndFDS) on host resistance, parasite virulence, or both, such that the benefit of a trait to fitness declines as its frequency increases. Polymorphism can also be stabilized by overdominance, when heterozygous hosts have greater resistance than homozygotes to diverse pathogens. Genetic diversity can also persist in the form of statistical polymorphism, sustained by random processes acting on gene frequencies and population size. Stable polymorphism allows alleles to be long-lived and genetic variation to be detectable in natural populations. In agriculture, many of the factors promoting stability in host-parasite interactions have been lost, leading to arms races of host defenses and parasite effectors.

The isolation and characterization of Pik, a rice blast resistance gene which emerged after rice domestication

• The rice-rice blast pathosystem is of great interest, not only because of the damaging potential of rice blast to the rice crop, but also because both the pathogen and its host are experimentally amenable. The rice blast resistance gene Pik, which is one of the five classical alleles located at the Pik locus on the long arm of chromosome 11, confers high and stable resistance to many Chinese rice blast isolates. • The isolation and functional characterization of Pik were performed in the present study through genetic and genomic approaches. • A combination of Pik-1 and Pik-2 is required for the expression of Pik resistance. Both Pik-1 and Pik-2 encode coiled-coil nucleotide binding site leucine-rich repeat (NBS-LRR) proteins, and each shares a very high level of protein identity with corresponding proteins encoded by the Pik-m and Pik-p alleles. Pik could be distinguished from other Pik alleles, including Pik-m and Pik-p, by the allele-specific, single-nucleotide polymorphism T1-2944G. • The coupled genes probably did not evolve as a result of a duplication event, and are far from any NBS-LRR R gene characterized. Pik is a younger allele at the locus that probably emerged after rice domestication.

Rapid linkage disequilibrium decay in the Lr10 gene in wild emmer wheat (Triticum dicoccoides) populations

INTRODUCTION

Recombination is a key evolutionary factor enhancing diversity. However, the effect of recombination on diversity in inbreeding species is expected to be low. To estimate this effect, recombination and diversity patterns of Lr10 gene were studied in natural populations of the inbreeder species, wild emmer wheat (Triticum dicoccoides). Wild emmer wheat is the progenitor of most cultivated wheats and it harbors rich genetic resources for disease resistance. Lr10 is a leaf rust resistance gene encoding three domains: a coiled-coil, nucleotide-binding site, and leucine-rich repeat (CC-NBS-LRR).

RESULTS

Lr10 was sequenced from 58 accessions representing 12 diverse habitats in Israel. Diversity analysis revealed a high rate of synonymous and non-synonymous substitutions (d (S) = 0.029, d (N) = 0.018, respectively) in the NBS-LRR domains. Moreover, in contrast to other resistance genes, in Lr10 the CC domain was more diverse than the NBS-LRR domains (d (S) = 0.069 vs. 0.029, d (N) = 0.094 vs. 0.018) and was subjected to positive selection in some of the populations. Seventeen recombination events were detected between haplotypes, especially in the CC domain. Linkage disequilibrium (LD) analysis has shown a rapid decay from r (2) = 0.5 to r (2) = 0.1 within a 2-kb span.

CONCLUSION

These results suggest that recombination is a diversifying force for the R-gene, Lr10, in the selfing species T. dicoccoides. This is the first report of a short-range LD decay in wild emmer wheat.

Towards identifying genes underlying ecologically relevant traits in Arabidopsis thaliana

A major challenge in evolutionary biology and plant breeding is to identify the genetic basis of complex quantitative traits, including those that contribute to adaptive variation. Here we review the development of new methods and resources to fine-map intraspecific genetic variation that underlies natural phenotypic variation in plants. In particular, the analysis of 107 quantitative traits reported in the first genome-wide association mapping study in Arabidopsis thaliana sets the stage for an exciting time in our understanding of plant adaptation. We also argue for the need to place phenotype-genotype association studies in an ecological context if one is to predict the evolutionary trajectories of plant species.

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.

The CC-NB-LRR-type Rdg2a resistance gene confers immunity to the seed-borne barley leaf stripe pathogen in the absence of hypersensitive cell death

Background

Leaf stripe disease on barley (Hordeum vulgare) is caused by the seed-transmitted hemi-biotrophic fungus Pyrenophora graminea. Race-specific resistance to leaf stripe is controlled by two known Rdg (Resistance to Drechslera graminea) genes: the H. spontaneum-derived Rdg1a and Rdg2a, identified in H. vulgare. The aim of the present work was to isolate the Rdg2a leaf stripe resistance gene, to characterize the Rdg2a locus organization and evolution and to elucidate the histological bases of Rdg2a-based leaf stripe resistance.

Principal Findings

We describe here the positional cloning and functional characterization of the leaf stripe resistance gene Rdg2a. At the Rdg2a locus, three sequence-related coiled-coil, nucleotide-binding site, and leucine-rich repeat (CC-NB-LRR) encoding genes were identified. Sequence comparisons suggested that paralogs of this resistance locus evolved through recent gene duplication, and were subjected to frequent sequence exchange. Transformation of the leaf stripe susceptible cv. Golden Promise with two Rdg2a-candidates under the control of their native 5′ regulatory sequences identified a member of the CC-NB-LRR gene family that conferred resistance against the Dg2 leaf stripe isolate, against which the Rdg2a-gene is effective. Histological analysis demonstrated that Rdg2a-mediated leaf stripe resistance involves autofluorescing cells and prevents pathogen colonization in the embryos without any detectable hypersensitive cell death response, supporting a cell wall reinforcement-based resistance mechanism.

Conclusions

This work reports about the cloning of a resistance gene effective against a seed borne disease. We observed that Rdg2a was subjected to diversifying selection which is consistent with a model in which the R gene co-evolves with a pathogen effector(s) gene. We propose that inducible responses giving rise to physical and chemical barriers to infection in the cell walls and intercellular spaces of the barley embryo tissues represent mechanisms by which the CC-NB-LRR-encoding Rdg2a gene mediates resistance to leaf stripe in the absence of hypersensitive cell death.

Diversity, distribution, and evolution of Solanum bulbocastanum late blight resistance genes

Knowledge on the evolution and distribution of late blight resistance genes is important for a better understanding of the dynamics of these genes in nature. We analyzed the presence and allelic diversity of the late blight resistance genes Rpi-blb1, Rpi-blb2, and Rpi-blb3, originating from Solanum bulbocastanum, in a set of tuber-bearing Solanum species comprising 196 different taxa. The three genes were only present in some Mexican diploid as well as polyploid species closely related to S. bulbocastanum. Sequence analysis of the fragments obtained from the Rpi-blb1 and Rpi-blb3 genes suggests an evolution through recombinations and point mutations. For Rpi-blb2, only sequences identical to the cloned gene were found in S. bulbocastanum accessions, suggesting that it has emerged recently. The three resistance genes occurred in different combinations and frequencies in S. bulbocastanum accessions and their spread is confined to Central America. A selected set of genotypes was tested for their response to the avirulence effectors IPIO-2, Avr-blb2, and Pi-Avr2, which interact with Rpi-blb1, Rpi-blb2, and Rpi-blb3, respectively, as well as by disease assays with a diverse set of isolates. Using this approach, some accessions could be identified that contain novel, as yet unknown, late blight resistance factors in addition to the Rpi-blb1, Rpi-blb2, and Rpi-blb3 genes.

Towards population genomics of effector-effector target interactions

Pathogen-plant host coevolutionary interactions exert strong natural selection on both organisms, specifically on the genes coding for effectors (pathogens), as well as on those coding for effector targets and R proteins (plant hosts). Natural selection leaves behind DNA sequence signatures on such genes and on linked genomic regions. These signatures can readily be detected by studying the patterns of intraspecies polymorphisms and interspecies divergence of the DNA sequences. Recent developments in DNA sequencing technology have made whole-genome studies on patterns of DNA polymorphisms : divergence possible. This type of analysis, called 'population genomics', appears to be powerful enough to identify novel effector-effector target genes. Here, we provide an overview of the statistical tools used for population genomics and their applications. This is followed by a brief review of evolutionary studies on plant genes involved in host-pathogen interactions. Finally we provide an example from our study on Magnaporthe oryzae.

In search of decoy/guardee to R genes: deciphering the role of sugars in defense against Fusarium wilt in chickpea

Plant responses are coordinately controlled by both external and internal signals. Apt perception of pathogen attack and its appropriate conversion to internal signals ultimately determine the outcome of innate immunity. The present review predicts the involvement of unconventional 'Guard/Decoy Model' in chickpea-Fusarium encounter. Rapid alkalinization factor is predicted to act as initial 'Gatekeeper decoy' counteracting fungal entry. Phospholipases and cystatins probably function as 'Guardees' being shielded by R gene(s). Serine Threonine Kinases decodes external pathogenic signals to in planta defense alarms. 14.3.3 provides clues to the wilt mechanism. The versatile sugars serve as signal generators and transmitters maintaining intra and inter cellular connectivity during stress.

Human IRGM gene "to be or not to be"

The immunity-related GTPases (IRG proteins) are one of the strongest early resistance systems against intracellular pathogens. The IRG gene family contains 21 copies arranged as tandem gene clusters on two chromosomes in the C57BL/6 mouse genome but has been reduced to only two copies in humans: IRGC and IRGM. IRGC is not involved in immunity, but the human IRGM gene plays a role in autophagy-targeted destruction of Mycobacterium tuberculosis (BCG) and Salmonella typhimurium. Variant IRGM haplotypes have been associated with increased risk for Crohn's disease and correlated with differential expression of IRGM transcripts. This article reviews in detail the studies performed on human samples, in vitro, and in sequence analyses that provide evidence for the unusual evolutionary history of the IRGM locus and the important role of the IRGM gene in autophagy and Crohn's disease in response to pathogenesis.

Structure, allelic diversity and selection of Asr genes, candidate for drought tolerance, in Oryza sativa L. and wild relatives

Asr (ABA, stress, ripening) genes represent a small gene family potentially involved in drought tolerance in several plant species. To analyze their interest for rice breeding for water-limited environments, this gene family was characterized further. Genomic organization of the gene family reveals six members located on four different chromosomes and with the same exon-intron structure. The maintenance of six members of the Asr gene family, which are the result of combination between tandem duplication and whole genome duplication, and their differential regulation under water stress, involves probably some sub-functionalization. The polymorphism of four members was studied in a worldwide collection of 204 accessions of Oryza sativa L. and 14 accessions of wild relatives (O. rufipogon and O. nivara). The nucleotide diversity of the Asr genes was globally low, but contrasted for the different genes, leading to different shapes of haplotype networks. Statistical tests for neutrality were used and compared to their distribution in a set of 111 reference genes spread across the genome, derived from another published study. Asr3 diversity exhibited a pattern concordant with a balancing selection at the species level and with a directional selection in the tropical japonica sub-group. This study provides a thorough description of the organization of the Asr family, and the nucleotide and haplotype diversity of four Asr in Oryza sativa species. Asr3 stood out as the best potential candidate. The polymorphism detected here represents a first step towards an association study between genetic polymorphisms of this gene family and variation in drought tolerance traits.

Natural diversity in flowering responses of Arabidopsis thaliana caused by variation in a tandem gene array

Tandemly arrayed genes that belong to gene families characterize genomes of many organisms. Gene duplication and subsequent relaxation of selection can lead to the establishment of paralogous cluster members that may evolve along different trajectories. Here, we report on the structural variation in MADS AFFECTING FLOWERING 2 (MAF2) gene, one member of the tandemly duplicated cluster of MADS-box-containing transcription factors in Arabidopsis thaliana. The altered gene structure at the MAF2 locus is present as a moderate-frequency polymorphism in Arabidopsis and leads to the extensive diversity in transcript patterns due to alternative splicing. Rearrangements at the MAF2 locus are associated with an early flowering phenotype in BC(5) lines. The lack of suppression of flowering time in a MAF2-insertion line expressing the MAF2-specific artificial miRNA suggests that these MAF2 variants are behaving as loss-of-function alleles. The variation in gene architecture is also associated with segregation distortion, which may have facilitated the spread and the establishment of the corresponding alleles throughout the Eurasian range of the A. thaliana population.

Single feature polymorphism (SFP)-based selective sweep identification and association mapping of growth-related metabolic traits in Arabidopsis thaliana

Background

Natural accessions of Arabidopsis thaliana are characterized by a high level of phenotypic variation that can be used to investigate the extent and mode of selection on the primary metabolic traits. A collection of 54 A. thaliana natural accession-derived lines were subjected to deep genotyping through Single Feature Polymorphism (SFP) detection via genomic DNA hybridization to Arabidopsis Tiling 1.0 Arrays for the detection of selective sweeps, and identification of associations between sweep regions and growth-related metabolic traits.

Results

A total of 1,072,557 high-quality SFPs were detected and indications for 3,943 deletions and 1,007 duplications were obtained. A significantly lower than expected SFP frequency was observed in protein-, rRNA-, and tRNA-coding regions and in non-repetitive intergenic regions, while pseudogenes, transposons, and non-coding RNA genes are enriched with SFPs. Gene families involved in plant defence or in signalling were identified as highly polymorphic, while several other families including transcription factors are depleted of SFPs. 198 significant associations between metabolic genes and 9 metabolic and growth-related phenotypic traits were detected with annotation hinting at the nature of the relationship. Five significant selective sweep regions were also detected of which one associated significantly with a metabolic trait.

Conclusions

We generated a high density polymorphism map for 54 A. thaliana accessions that highlights the variability of resistance genes across geographic ranges and used it to identify selective sweeps and associations between metabolic genes and metabolic phenotypes. Several associations show a clear biological relationship, while many remain requiring further investigation.

Arabidopsis and the plant immune system

Understanding the fundamental mechanisms of plant disease resistance is of central importance to sustainable agriculture and human health. Use of the model plant Arabidopsis thaliana has resulted in an explosion of information regarding both disease resistance and susceptibility to pathogens. The last 20 years of research have demonstrated the commonalities between Arabidopsis and crop species. In this review, commemorating the 10th anniversary of the sequencing of the Arabidopsis genome, we will address some of the insights derived from the use of Arabidopsis as a model plant pathology system.

Remarkably ancient balanced polymorphisms in a multi-locus gene network

Local adaptations within species are often governed by several interacting genes scattered throughout the genome. Single-locus models of selection cannot explain the maintenance of such complex variation because recombination separates co-adapted alleles. Here we report a novel type of intraspecific multi-locus genetic variation that has been maintained over a vast period of time. The galactose (GAL) utilization gene network of the brewer’s yeast relative Saccharomyces kudriavzevii exists in two distinct states: a functional gene network in Portuguese strains and, in Japanese strains, a non-functional gene network of allelic pseudogenes. Genome sequencing of all available S. kudriavzevii strains revealed that none of the functional GAL genes were acquired from other species. Rather, these polymorphisms have been maintained for nearly the entire history of the species, despite more recent gene flow genome-wide. Experimental evidence suggests that inactivation of the GAL3 and GAL80 regulatory genes facilitated the origin and long-term maintenance of the two gene network states. This striking example of a balanced unlinked gene network polymorphism introduces a remarkable type of intraspecific variation that may be widespread.

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.

Clustering, haplotype diversity and locations of MIC-3: a unique root-specific defense-related gene family in Upland cotton (Gossypium hirsutum L.)

MIC-3 is a recently identified gene family shown to exhibit increased root-specific expression following nematode infection of cotton plants that are resistant to root-knot nematode. Here, we cloned and sequenced MIC-3 genes from selected diploid and tetraploid cotton species to reveal sequence differences at the molecular level and identify chromosomal locations of MIC-3 genes in Gossypium species. Detailed sequence analysis and phylogenetic clustering of MIC-3 genes indicated the presence of multiple MIC-3 gene members in Gossypium species. Haplotypes of a MIC-3 gene family member were discovered by comparative analysis among consensus sequences across genotypes within an individual clade in the phylogram to overcome the problem of duplicated loci in the tetraploid cotton. Deficiency tests of the SNPs delimited six A(t)-genome members of the MIC-3 family clustered to chromosome arm 4sh, and one D(t)-genome member to chromosome 19. Clustering was confirmed by long-PCR amplification of the intergenic regions using A(t)-genome-specific MIC-3 primer pairs. The clustered distribution may have been favored by selection for responsiveness to evolving disease and/or pest pressures, because large variants of the MIC-3 gene family may have been recovered from small physical areas by recombination. This could give a buffer against selection pressure from a broad range of pest and pathogens in the future. To our knowledge, these are the first results on the evolution of clustering and genome-specific haplotype members of a unique cotton gene family associated with resistant response against a major pathogen.

Targets of balancing selection in the human genome

Balancing selection is potentially an important biological force for maintaining advantageous genetic diversity in populations, including variation that is responsible for long-term adaptation to the environment. By serving as a means to maintain genetic variation, it may be particularly relevant to maintaining phenotypic variation in natural populations. Nevertheless, its prevalence and specific targets in the human genome remain largely unknown. We have analyzed the patterns of diversity and divergence of 13,400 genes in two human populations using an unbiased single-nucleotide polymorphism data set, a genome-wide approach, and a method that incorporates demography in neutrality tests. We identified an unbiased catalog of genes with signatures of long-term balancing selection, which includes immunity genes as well as genes encoding keratins and membrane channels; the catalog also shows enrichment in functional categories involved in cellular structure. Patterns are mostly concordant in the two populations, with a small fraction of genes showing population-specific signatures of selection. Power considerations indicate that our findings represent a subset of all targets in the genome, suggesting that although balancing selection may not have an obvious impact on a large proportion of human genes, it is a key force affecting the evolution of a number of genes in humans.

Quantitative fitness effects of infection in a gene-for-gene system

* It is often assumed that pathogen infection decreases plant fitness, thereby driving the evolution of plant resistance (R) genes. However, the impact of bacterial pathogens on fitness has been shown to be relatively subtle, ranging from positive to negative. * In this study, we focus on the Rps5-mediated resistance in Arabidopsis thaliana and examine the fitness effects of resistance by experimentally infecting resistant (R) and susceptible (S) plants with a natural avirulent Pseudomonas syringae strain at each of three initial infection dosage levels. Our methodology ensured control of the plant genetic backgrounds; within each of two natural accessions we created isolines varying in the presence or absence of Rps5. * In terms of lifetime fitness, R plants outperformed their S controls by 9.6-32% when infected by a pathogen carrying an associated Avr gene, depending on the initial dosage levels and genetic backgrounds. * We also found that the naturally R line, Col-0, is more tolerant than the naturally S accession, Ga-0. The negative impact of infection on fitness was 20% less in Col-0 than Ga-0. We found no effect of Rps5 itself on the tolerance of either accession. We therefore failed to find evidence for a trade-off between tolerance and resistance.

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.

The diversity of C-type lectins in the genome of a basal metazoan, Nematostella vectensis

C-type lectins (CTLs) are involved in cell-cell adhesion, recognition, and innate immunity in higher vertebrates, but little is known about CTLs in basal metazoans. The recent sequencing of the cnidarian Nematostella vectensis genome allowed us to explore the CTL-like gene family at the base of metazoan evolution. Sixty-seven predicted CTLs, with a total of 92 putative C-type lectin domains (CTLDs), were classified according to number of CTLDs present and their association with other protein domains in the CTL. Conserved residues in the glycan-binding pocket suggest that approximately half of the CTLDs retain glycan-binding function. Phylogenetic analysis of N. vectensis CTLDs with respect to other model invertebrates and humans indicates N. vectensis CTLD sequences more closely resemble vertebrate CTLDs. This study provides a N. vectensis CTL database that can be used for further research on the evolution of cnidarian CTLs and the role of CTLs in cnidarian innate immunity.

Continua of specificity and virulence in plant host-pathogen interactions: causes and consequences

Ecological, evolutionary and molecular models of interactions between plant hosts and microbial pathogens are largely based around a concept of tightly coupled interactions between species pairs. However, highly pathogenic and obligate associations between host and pathogen species represent only a fraction of the diversity encountered in natural and managed systems. Instead, many pathogens can infect a wide range of hosts, and most hosts are exposed to more than one pathogen species, often simultaneously. Furthermore, outcomes of pathogen infection vary widely because host plants vary in resistance and tolerance to infection, while pathogens are also variable in their ability to grow on or within hosts. Environmental heterogeneity further increases the potential for variation in plant host-pathogen interactions by influencing the degree and fitness consequences of infection. Here, we describe these continua of specificity and virulence inherent within plant host-pathogen interactions. Using this framework, we describe and contrast the genetic and environmental mechanisms that underlie this variation, outline consequences for epidemiology and community structure, explore likely ecological and evolutionary drivers, and highlight several key areas for future research.

The zig-zag-zig in oomycete-plant interactions

In addition to a range of preformed barriers, plants defend themselves against microbial invasion by detecting conserved, secreted molecules, called pathogen-associated molecular patterns (PAMPs). PAMP-triggered immunity (PTI) is the first inducible layer of plant defence that microbial pathogens must navigate by the delivery of effector proteins that act to suppress or otherwise manipulate key components of resistance. Effectors may themselves be targeted by a further layer of defence, effector-triggered immunity (ETI), as their presence inside or outside host cells may be detected by resistance proteins. This 'zig-zag-zig' of tightly co-evolving molecular interactions determines the outcome of attempted infection. In this article, we consider the complex molecular interplay between plants and plant pathogenic oomycetes, drawing on recent literature to illustrate what is known about oomycete PAMPs and elicitors of defence responses, the effectors they utilize to suppress PTI, and the phenomenal molecular 'battle' between effector and resistance (R) genes that dictates the establishment or evasion of ETI.

In the trenches of plant pathogen recognition: Role of NB-LRR proteins

As in nearly every discipline of plant biology, new insights are constantly changing our understanding of plant immunity. It is now clear that plant immunity is controlled by two layers of inducible responses: basal responses triggered by conserved microbial features and specific responses triggered by gene-for-gene recognition of pathogen effector proteins by host resistance (R) proteins. The nucleotide-binding domain leucine-rich repeat (NB-LRR) class of R proteins plays a major role in the combat against a wide range of plant pathogens. The variation that has been generated and is maintained within these conserved proteins has diversified their specificity, subcellular localisations, activation and recognition mechanisms, allowing them to specifically adapt to different plant-pathogen interaction systems. This review addresses recent advances in the molecular role of NB-LRR proteins in pathogen recognition and activation of plant defence responses.

Identification and analyses of candidate genes for rpp4-mediated resistance to Asian soybean rust in soybean

Asian soybean rust is a formidable threat to soybean (Glycine max) production in many areas of the world, including the United States. Only five sources of resistance have been identified (Resistance to Phakopsora pachyrhizi1 [Rpp1], Rpp2, Rpp3, Rpp4, and Rpp5). Rpp4 was previously identified in the resistant genotype PI459025B and mapped within 2 centimorgans of Satt288 on soybean chromosome 18 (linkage group G). Using simple sequence repeat markers, we developed a bacterial artificial chromosome contig for the Rpp4 locus in the susceptible cv Williams82 (Wm82). Sequencing within this region identified three Rpp4 candidate disease resistance genes (Rpp4C1-Rpp4C3 [Wm82]) with greatest similarity to the lettuce (Lactuca sativa) RGC2 family of coiled coil-nucleotide binding site-leucine rich repeat disease resistance genes. Constructs containing regions of the Wm82 Rpp4 candidate genes were used for virus-induced gene silencing experiments to silence resistance in PI459025B, confirming that orthologous genes confer resistance. Using primers developed from conserved sequences in the Wm82 Rpp4 candidate genes, we identified five Rpp4 candidate genes (Rpp4C1-Rpp4C5 [PI459025B]) from the resistant genotype. Additional markers developed from the Wm82 Rpp4 bacterial artificial chromosome contig further defined the region containing Rpp4 and eliminated Rpp4C1 (PI459025B) and Rpp4C3 (PI459025B) as candidate genes. Sequencing of reverse transcription-polymerase chain reaction products revealed that Rpp4C4 (PI459025B) was highly expressed in the resistant genotype, while expression of the other candidate genes was nearly undetectable. These data support Rpp4C4 (PI459025B) as the single candidate gene for Rpp4-mediated resistance to Asian soybean rust.

Functional and geographical differentiation of candidate balanced polymorphisms in Arabidopsis thaliana

Molecular population genetic analysis of three chromosomal regions in Arabidopsis thaliana suggested that balancing selection might operate to maintain variation at three novel candidate adaptive trait genes, including SOLUBLE STARCH SYNTHASE I (SSI), PLASTID TRANSCRIPTIONALLY ACTIVE 7(PTAC7), and BELL-LIKE HOMEODOMAIN 10 (BLH10). If balanced polymorphisms are indeed maintained at these loci, then we would expect to observe functional variation underlying the previously detected signatures of selection. We observe multiple replacement polymorphisms within and in the 32 amino acids just upstream of the protein-protein interacting BELL domain at the BLH10 locus. While no clear protein sequence differences are found between allele types in SSI and PTAC7, these two genes show evidence for allele-specific variation in expression levels. Geographical patterns of allelic differentiation seem consistent with population stratification in this species and a significant longitudinal cline was observed at all three candidate loci. These data support a hypothesis of balancing selection at all three candidate loci and provide a basis for more detailed functional work by identifying possible functional differences that might be selectively maintained.

Evolution of new disease specificity at a simple resistance locus in a crop-weed complex: reconstitution of the Lr21 gene in wheat

The wheat leaf-rust resistance gene Lr21 was first identified in an Iranian accession of goatgrass, Aegilops tauschii Coss., the D-genome donor of hexaploid bread wheat, and was introgressed into modern wheat cultivars by breeding. To elucidate the origin of the gene, we analyzed sequences of Lr21 and lr21 alleles from 24 wheat cultivars and 25 accessions of Ae. tauschii collected along the Caspian Sea in Iran and Azerbaijan. Three basic nonfunctional lr21 haplotypes, H1, H2, and H3, were identified. Lr21 was found to be a chimera of H1 and H2, which were found only in wheat. We attempted to reconstitute a functional Lr21 allele by crossing the cultivars Fielder (H1) and Wichita (H2). Rust inoculation of 5876 F(2) progeny revealed a single resistant plant that proved to carry the H1H2 haplotype, a result attributed to intragenic recombination. These findings reflect how plants balance the penalty and the necessity of a resistance gene and suggest that plants can reuse "dead" alleles to generate new disease-resistance specificity, leading to a "death-recycle" model of plant-resistance gene evolution at simple loci. We suggest that selection pressure in crop-weed complexes contributes to this process.

Possible diversifying selection in the imprinted gene, MEDEA, in Arabidopsis

Coevolutionary conflict among imprinted genes that influence traits such as offspring growth may arise when maternal and paternal genomes have different evolutionary optima. This conflict is expected in outcrossing taxa with multiple paternity, but not self-fertilizing taxa. MEDEA (MEA) is an imprinted plant gene that influences seed growth. Disagreement exists regarding the type of selection acting on this gene. We present new data and analyses of sequence diversity of MEA in self-fertilizing and outcrossing Arabidopsis and its relatives, to help clarify the form of selection acting on this gene. Codon-based branch analysis among taxa (PAML) suggests that selection on the coding region is changing over time, and nonsynonymous substitution is elevated in at least one outcrossing branch. Codon-based analysis of diversity within outcrossing Arabidopsis lyrata ssp. petraea (OmegaMap) suggests that diversifying selection is acting on a portion of the gene, to cause elevated nonsynonymous polymorphism. Providing further support for balancing selection in A. lyrata, Hudson, Kreitman and Aguadé analysis indicates that diversity/divergence at silent sites in the MEA promoter and genic region is elevated relative to reference genes, and there are deviations from the neutral frequency spectrum. This combination of positive selection as well as balancing and diversifying selection in outcrossing lineages is consistent with other genes influence by evolutionary conflict, such as disease resistance genes. Consistent with predictions that conflict would be eliminated in self-fertilizing taxa, we found no evidence of positive, balancing, or diversifying selection in A. thaliana promoter or genic region.