Signatures of balancing selection are maintained at disease resistance loci following mating system evolution and a population bottleneck in the genus Capsella

ABSTRACT: BACKGROUND: Population bottlenecks can lead to a loss of variation at disease resistance loci, which could have important consequences for the ability of populations to adapt to pathogen pressure. Alternatively, current or past balancing selection could maintain high diversity, creating a strong heterogeneity in the retention of polymorphism across the genome of bottlenecked populations. We sequenced part of the LRR region of 9 NBS-LRR disease resistance genes in the outcrossing Capsella grandiflora and the recently derived, bottlenecked selfing species Capsella rubella, and compared levels and patterns of nucleotide diversity and divergence with genome-wide reference loci. RESULTS: In strong contrast with reference loci, average diversity at resistance loci was comparable between C. rubella and C. grandiflora, primarily due to two loci with highly elevated diversity indicative of past or present balancing selection. Average between-species differentiation was also reduced at the set of R-genes compared with reference loci, which is consistent with the maintenance of ancestral polymorphism. CONCLUSIONS: Historical or ongoing balancing selection on plant disease resistance genes is a likely contributor to the retention of ancestral polymorphism in some regions of the bottlenecked Capella rubella genome.     

MHC-mediated local adaptation in reciprocally translocated Chinook salmon

Most Pacific salmonid populations have faced significant population declines over the past 30 years. In order to effectively conserve and manage these populations, knowledge of the evolutionary adaptive state of individuals and the scale of adaptation across populations is needed. The vertebrate major histocompatibility complex (MHC) represents an important adaptation to parasites, and genes encoding for the MHC are widely held to be undergoing balancing selection. However, the generality of balancing selection across populations at MHC loci is not well documented. Using Chinook salmon (Oncorhynchus tshawytscha) from two populations, we follow the survival of full-sib family replicates reared in their natal river and reciprocally transplanted to a foreign river to examine selection and local adaptation at the MHC class I and II loci. In both populations, we found evidence of a survivorship advantage associated with nucleotide diversity at the MHC class I locus. In contrast, we found evidence that MHC class II diversity was disadvantageous in one population. There was no evidence that these effects occurred in translocated families, suggesting some degree of local adaptation at the MHC loci. Thus, our results implicate balancing selection at the MHC class I but potentially differing selection across populations at the class II locus.     

Balancing Selection for Pathogen Resistance Reveals an Intercontinental Signature of Red Queen Coevolution

The link between long-term host–parasite coevolution and genetic diversity is key to understanding genetic epidemiology and the evolution of resistance. The model of Red Queen host–parasite coevolution posits that high genetic diversity is maintained when rare host resistance variants have a selective advantage, which is believed to be the mechanistic basis for the extraordinarily high levels of diversity at disease-related genes such as the major histocompatibility complex in jawed vertebrates and R-genes in plants. The parasites that drive long-term coevolution are, however, often elusive. Here we present evidence for long-term balancing selection at the phenotypic (variation in resistance) and genomic (resistance locus) level in a particular host–parasite system: the planktonic crustacean Daphnia magna and the bacterium Pasteuria ramosa. The host shows widespread polymorphisms for pathogen resistance regardless of geographic distance, even though there is a clear genome-wide pattern of isolation by distance at other sites. In the genomic region of a previously identified resistance supergene, we observed consistent molecular signals of balancing selection, including higher genetic diversity, older coalescence times, and lower differentiation between populations, which set this region apart from the rest of the genome. We propose that specific long-term coevolution by negative-frequency-dependent selection drives this elevated diversity at the host''s resistance loci on an intercontinental scale and provide an example of a direct link between the host’s resistance to a virulent pathogen and the large-scale diversity of its underlying genes.     

Geographical distribution of selected and putatively neutral SNPs in Southeast Asian malaria parasites

Loci targeted by directional selection are expected to show elevated geographical population structure relative to neutral loci, and a flurry of recent papers have used this rationale to search for genome regions involved in adaptation. Studies of functional mutations that are known to be under selection are particularly useful for assessing the utility of this approach. Antimalarial drug treatment regimes vary considerably between countries in Southeast Asia selecting for local adaptation at parasite loci underlying resistance. We compared the population structure revealed by 10 nonsynonymous mutations (nonsynonymous single-nucleotide polymorphisms [nsSNPs]) in four loci that are known to be involved in antimalarial drug resistance, with patterns revealed by 10 synonymous mutations (synonymous single-nucleotide polymorphisms [sSNPs]) in housekeeping genes or genes of unknown function in 755 Plasmodium falciparum infections collected from 13 populations in six Southeast Asian countries. Allele frequencies at known nsSNPs underlying resistance varied markedly between locations (F(ST) = 0.18-0.66), with the highest frequencies on the Thailand-Burma border and the lowest frequencies in neighboring Lao PDR. In contrast, we found weak but significant geographic structure (F(ST) = 0-0.14) for 8 of 10 sSNPs. Importantly, all 10 nsSNPs showed significantly higher F(ST) (P < 8 x 10(-5)) than simulated neutral expectations based on observed F(ST) values in the putatively neutral sSNPs. This result was unaffected by the methods used to estimate allele frequencies or the number of populations used in the simulations. Given that dense single-nucleotide polymorphism (SNP) maps and rapid SNP assay methods are now available for P. falciparum, comparing genetic differentiation across the genome may provide a valuable aid to identifying parasite loci underlying local adaptation to drug treatment regimes or other selective forces. However, the high proportion of polymorphic sites that appear to be under balancing selection (or linked to selected sites) in the P. falciparum genome violates the central assumption that selected sites are rare, which complicates identification of outlier loci, and suggests that caution is needed when using this approach.     

Geographic heterogeneity in natural selection on an MHC locus in sockeye salmon

Balancing selection maintains high levels of polymorphism and heterozygosity in genes of the MHC (major histocompatibility complex) of vertebrate organisms, and promotes long evolutionary persistence of individual alleles and strongly differentiated allelic lineages. In this study, genetic variation at the MHC class II DAB-beta1 locus was examined in 31 populations of sockeye salmon (Oncorhynchus nerka) inhabiting the Fraser River drainage of British Columbia, Canada. Twenty-five percent of variation at the locus was partitioned among sockeye populations, as compared with 5% at neutral genetic markers. Geographic heterogeneity of balancing selection was detected among four regions in the Fraser River drainage and among lake systems within regions. High levels of beta1 allelic diversity and heterozygosity, as well as distributions of alleles and allelic lineages that were more even than expected for a neutral locus, indicated the presence of balancing selection in populations throughout much of the interior Fraser drainage. However, proximate populations in the upper Fraser region, and four of six populations from the lower Fraser drainage, exhibited much lower levels of genetic diversity and had beta1 allele frequency distributions in conformance with those expected for a neutral locus, or a locus under directional selection. Pair-wise FST values for beta1 averaged 0.19 and tended to exceed the corresponding values estimated for neutral loci at all levels of population structure, although they were lower among populations experiencing balancing selection than among other populations. The apparent heterogeneity in selection resulted in strong genetic differentiation between geographically proximate populations with and without detectable levels of balancing selection, in stark contrast to observations at neutral loci. The strong partitioning and complex structure of beta1 diversity within and among sockeye populations on a small geographic scale illustrates the value of incorporating adaptive variation into conservation planning for the species.     

Comparative analysis of population structure across environments and geographical scales at major histocompatibility complex and microsatellite loci in Atlantic salmon (Salmo salar)

Evidence of selection acting on major histocompatibility complex (MHC) genes has been illustrated with the analysis of their nucleotide sequences and allele frequency distribution. Comparing the patterns of population differentiation at neutral markers and MHC genes in the wild may provide further insights about the relative role of selection and neutrality in shaping their diversity. In this study, we combine both methods to assess the role of selection on a MHC gene in Atlantic salmon. We compare variation at a MHC class II B locus and microsatellites among 14 samples from seven different rivers and seven subpopulations within a single river system covering a variety of habitats and different geographical scales. We show that diversifying selection is acting on the sites involved in antigen presentation and that balancing selection maintains a high level of polymorphism within populations. Despite important differences in habitat type, the comparison of the population structure at MHC and microsatellites on large geographical scales reveals a correlation between patterns of differentiation, indicating that drift and migration have been more important than selection in shaping population differentiation at the MHC locus. In contrast, strong discrepancies between patterns of population differentiation at the two types of markers provides support for the role of selection in shaping population structure within rivers. Together, these results confirm that natural selection is influencing MHC gene diversity in wild Atlantic salmon although neutral forces may also be important in their evolution.     

A temporal analysis shows major histocompatibility complex loci in the Scandinavian wolf population are consistent with neutral evolution

The major histocompatibility complex (MHC) has an integral role in the immune system, and hence diversity at its genes may be of particular importance for the health of populations. In large populations, balancing selection maintains diversity in MHC genes, but theoretical expectations indicate that this form of selection is absent or inefficient in small populations. We examine the level of diversity at three MHC class II loci in the wolf population of Scandinavia, a population naturally recolonized with a genetic contribution from as few as three founders, and in four neighbouring wolf populations. In the Scandinavian wolf population, two alleles were found for each locus and the distribution of alleles is compatible with their linkage into two haplotypes. Changes in the level of heterozygosity over time since recolonization demonstrate the effects of the proposed arrival of an immigrant wolf. The maintenance of diversity is shown to be compatible with a neutral, random allocation of alleles, in conjunction with crossing between packs. A total of 15 DRB1, seven DQA and 10 DQB1 alleles are found in four neighbouring wolf populations, with substantial sharing across populations. Even in these larger populations, bottlenecks and fragmentation with consequent genetic drift are likely to have resulted in few indicators for balancing selection and significant differentiation of populations.     

Evolving disease resistance genes

Defenses against most specialized plant pathogens are often initiated by a plant disease resistance gene. Plant genomes encode several classes of genes that can function as resistance genes. Many of the mechanisms that drive the molecular evolution of these genes are now becoming clear. The processes that contribute to the diversity of R genes include tandem and segmental gene duplications, recombination, unequal crossing-over, point mutations, and diversifying selection. Diversity within populations is maintained by balancing selection. Analyses of whole-genome sequences have and will continue to provide new insight into the dynamics of resistance gene evolution.     

Local Adaptation in the Flowering-Time Gene Network of Balsam Poplar, Populus balsamifera L

Identifying the signature and targets of local adaptation is an increasingly important goal in empirical population genetics. Using data from 443 balsam poplar Populus balsamifera trees sampled from 31 populations, we tested for evidence of geographically variable selection shaping diversity at 27 homologues of the Arabidopsis flowering-time network. These genes are implicated in the control of seasonal phenology, an important determinant of fitness. Using 335 candidate and 412 reference single nucleotide polymorphisms (SNPs), we tested for evidence of local adaptation by searching for elevated population differentiation using F(ST)-based outlier analyses implemented in BayeScan or a Hierarchical Model in Arelquin and by testing for significant associations between allele frequency and environmental variables using BAYENV. A total of 46 SNPs from 14 candidate genes had signatures of local adaptation-either significantly greater population differentiation or significant covariance with one or more environmental variable relative to reference SNP distributions. Only 11 SNPs from two genes exhibited both elevated population differentiation and covariance with one or more environmental variables. Several genes including the abscisic acid gene ABI1B and the circadian clock genes ELF3 and GI5 harbored a large number of SNPs with signatures of local adaptation-with SNPs in GI5 strongly covarying with both latitude and precipitation and SNPs in ABI1B strongly covarying with temperature. In contrast to several other systems, we find little evidence that photoreceptors, including phytochromes, play an important role in local adaptation. Our results additionally show that detecting local adaptation is sensitive to the analytical approaches used and that model-based significance thresholds should be viewed with caution.     

Unraveling the Effects of Selection and Demography on Immune Gene Variation in Free-Ranging Plains Zebra (Equus quagga) Populations

Demography, migration and natural selection are predominant processes affecting the distribution of genetic variation among natural populations. Many studies use neutral genetic markers to make inferences about population history. However, the investigation of functional coding loci, which directly reflect fitness, is critical to our understanding of species'' ecology and evolution. Immune genes, such as those of the Major Histocompatibility Complex (MHC), play an important role in pathogen recognition and provide a potent model system for studying selection. We contrasted diversity patterns of neutral data with MHC loci, ELA-DRA and -DQA, in two southern African plains zebra (Equus quagga) populations: Etosha National Park, Namibia, and Kruger National Park, South Africa. Results from neutrality tests, along with observations of elevated diversity and low differentiation across populations, supported previous genus-level evidence for balancing selection at these loci. Despite being low, MHC divergence across populations was significant and may be attributed to drift effects typical of geographically separated populations experiencing little to no gene flow, or alternatively to shifting allele frequency distributions driven by spatially variable and fluctuating pathogen communities. At the DRA, zebra exhibited geographic differentiation concordant with microsatellites and reduced levels of diversity in Etosha due to highly skewed allele frequencies that could not be explained by demography, suggestive of spatially heterogeneous selection and local adaptation. This study highlights the complexity in which selection affects immune gene diversity and warrants the need for further research on the ecological mechanisms shaping patterns of adaptive variation among natural populations.     

Balancing selection is the main force shaping the evolution of innate immunity genes

The evolutionarily recent geographic expansion of humans, and the even more recent development of large, relatively dense human settlements, has exposed our species to new pathogenic environments. Potentially lethal pathogens are likely to have exerted important selective pressures on our genome, so immunity genes can be expected to show molecular signatures of the adaptation of human populations to these recent conditions. While genes related to the acquired immunity system have indeed been reported to show traces of local adaptation, little is known about the response of the innate immunity system. In this study, we analyze the variability patterns in different human populations of fifteen genes related to innate immunity. We have used both single nucleotide polymorphism and sequence data, and through the analysis of interpopulation differentiation, the linkage disequilibrium pattern, and intrapopulation diversity, we have discovered some signatures of positive and especially balancing selection in these genes, thus confirming the importance of the immune system genetic plasticity in the evolutionary adaptive process. Interestingly, the strongest evidence is found in three TLR genes and CD14. These innate immunity genes play a pivotal role, being involved in the primary recognition of pathogens. In general, more evidences of selection appear in the European populations, in some case possibly related to severe population specific pressures. However, we also describe evidence from African populations, which may reflect parallel or long-term selective forces acting in different geographic areas.     

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.     

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

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

Diversity and linkage of genes in the self-incompatibility gene family in Arabidopsis lyrata

We report studies of seven members of the S-domain gene family in Arabidopsis lyrata, a member of the Brassicaceae that has a sporophytic self-incompatibility (SI) system. Orthologs for five loci are identifiable in the self-compatible relative A. thaliana. Like the Brassica stigmatic incompatibility protein locus (SRK), some of these genes have kinase domains. We show that several of these genes are unlinked to the putative A. lyrata SRK, Aly13. These genes have much lower nonsynonymous and synonymous polymorphism than Aly13 in the S-domains within natural populations, and differentiation between populations is higher, consistent with balancing selection at the Aly13 locus. One gene (Aly8) is linked to Aly13 and has high diversity. No departures from neutrality were detected for any of the loci. Comparing different loci within A. lyrata, sites corresponding to hypervariable regions in the Brassica S-loci (SLG and SRK) and in comparable regions of Aly13 have greater replacement site divergence than the rest of the S-domain. This suggests that the high polymorphism in these regions of incompatibility loci is due to balancing selection acting on sites within or near these regions, combined with low selective constraints.     

Genetic variation affecting host-parasite interactions: major-effect quantitative trait loci affect the transmission of sigma virus in Drosophila melanogaster

In natural populations, genetic variation affects resistance to disease. Whether that genetic variation comprises lots of small-effect polymorphisms or a small number of large-effect polymorphisms has implications for adaptation, selection and how genetic variation is maintained in populations. Furthermore, how much genetic variation there is, and the genes that underlie this variation, affects models of co-evolution between parasites and their hosts. We are studying the genetic variation that affects the resistance of Drosophila melanogaster to its natural pathogen — the vertically transmitted sigma virus. We have carried out three separate quantitative trait locus mapping analyses to map gene variants on the second chromosome that cause variation in the rate at which males transmit the infection to their offspring. All three crosses identified a locus in a similar chromosomal location that causes a large drop in the rate at which the virus is transmitted. We also found evidence for an additional smaller-effect quantitative trait locus elsewhere on the chromosome. Our data, together with previous experiments on the sigma virus and parasitoid wasps, indicate that the resistance of D. melanogaster to co-evolved pathogens is controlled by a limited number of major-effect polymorphisms.     

Genetic analysis of differentiation among breeding ponds reveals a candidate gene for local adaptation in Rana arvalis

One of the main questions in evolutionary and conservation biology is how geographical and environmental features of the landscape shape neutral and adaptive genetic variation in natural populations. The identification of genomic polymorphisms that account for adaptive variation can aid in finding candidate loci for local adaptation. Consequently, a comparison of spatial patterns in neutral markers and loci under selection may help disentangle the effects of gene flow, genetic drift and selection at the landscape scale. Many amphibians breed in wetlands, which differ in environmental conditions and in the degree of isolation, enhancing the potential for local adaptation. We used microsatellite markers to measure genetic differentiation among 17 local populations of Rana arvalis breeding in a network of wetlands. We found that locus RC08604 deviated from neutral expectations, suggesting that it is a good candidate for directional selection. We used a genetic network analysis to show that the allele distribution in this locus is correlated with habitat characteristics, whereas this was not the case at neutral markers that displayed a different allele distribution and population network in the study area. The graph approach illustrated the genomic heterogeneity (neutral loci vs. the candidate locus for directional selection) of gene exchange and genetic divergence among populations under directional selection. Limited gene flow between wetlands was only observed at the candidate genomic region under directional selection. RC08604 is partially located inside an up‐regulated thyroid‐hormone receptor (TRβ) gene coordinating the expression of other genes during metamorphosis and appears to be linked with variation in larval life‐history traits found among R. arvalis populations. We suggest that directional selection on genes coding larval life‐history traits is strong enough to maintain the divergence in these genomic regions, reducing the effective recombination of locally adapted alleles but not in other regions of the genome. Integrating this knowledge into conservation plans at the landscape scale will improve the design of management strategies to preserve adaptive genetic diversity in wetland networks.     

Distinguishing between Selective Sweeps from Standing Variation and from a De Novo Mutation

An outstanding question in human genetics has been the degree to which adaptation occurs from standing genetic variation or from de novo mutations. Here, we combine several common statistics used to detect selection in an Approximate Bayesian Computation (ABC) framework, with the goal of discriminating between models of selection and providing estimates of the age of selected alleles and the selection coefficients acting on them. We use simulations to assess the power and accuracy of our method and apply it to seven of the strongest sweeps currently known in humans. We identify two genes, ASPM and PSCA, that are most likely affected by selection on standing variation; and we find three genes, ADH1B, LCT, and EDAR, in which the adaptive alleles seem to have swept from a new mutation. We also confirm evidence of selection for one further gene, TRPV6. In one gene, G6PD, neither neutral models nor models of selective sweeps fit the data, presumably because this locus has been subject to balancing selection.     

Hitch-hiking to a locus under balancing selection: high sequence diversity and low population subdivision at the S-locus genomic region in Arabidopsis halleri.

Hitch-hiking to a site under balancing selection is expected to produce a local increase in nucleotide polymorphism and a decrease in population differentiation compared with the background genomic level, but empirical evidence supporting these predictions is scarce. We surveyed molecular diversity at four genes flanking the region controlling self-incompatibility (the S-locus) in samples from six populations of the herbaceous plant Arabidopsis halleri, and compared their polymorphism with sequences from five control genes unlinked to the S-locus. As a preliminary verification, the S-locus flanking genes were shown to co-segregate with SRK, the gene involved in the self-incompatibility reaction at the pistil level. In agreement with theory, our results demonstrated a significant peak of nucleotide diversity around the S-locus as well as a significant decrease in population genetic structure in the S-locus region compared with both control genes and a set of seven unlinked microsatellite markers. This is consistent with the theoretical expectation that balancing selection is increasing the effective migration rate in subdivided populations. Although only four S-locus flanking genes were investigated, our results suggest that these two signatures of the hitch-hiking effect are localized in a very narrow genomic region.     

Balancing selection at the tomato RCR3 Guardee gene family maintains variation in strength of pathogen defense

Coevolution between hosts and pathogens is thought to occur between interacting molecules of both species. This results in the maintenance of genetic diversity at pathogen antigens (or so-called effectors) and host resistance genes such as the major histocompatibility complex (MHC) in mammals or resistance (R) genes in plants. In plant-pathogen interactions, the current paradigm posits that a specific defense response is activated upon recognition of pathogen effectors via interaction with their corresponding R proteins. According to the "Guard-Hypothesis," R proteins (the "guards") can sense modification of target molecules in the host (the "guardees") by pathogen effectors and subsequently trigger the defense response. Multiple studies have reported high genetic diversity at R genes maintained by balancing selection. In contrast, little is known about the evolutionary mechanisms shaping the guardee, which may be subject to contrasting evolutionary forces. Here we show that the evolution of the guardee RCR3 is characterized by gene duplication, frequent gene conversion, and balancing selection in the wild tomato species Solanum peruvianum. Investigating the functional characteristics of 54 natural variants through in vitro and in planta assays, we detected differences in recognition of the pathogen effector through interaction with the guardee, as well as substantial variation in the strength of the defense response. This variation is maintained by balancing selection at each copy of the RCR3 gene. Our analyses pinpoint three amino acid polymorphisms with key functional consequences for the coevolution between the guardee (RCR3) and its guard (Cf-2). We conclude that, in addition to coevolution at the "guardee-effector" interface for pathogen recognition, natural selection acts on the "guard-guardee" interface. Guardee evolution may be governed by a counterbalance between improved activation in the presence and prevention of auto-immune responses in the absence of the corresponding pathogen.