Spatially and temporally fluctuating selection at non-MHC immune genes: evidence from TAP polymorphism in populations of brown trout (Salmo trutta, L.)

ddRAD-seq reveals the genetic structure and detects signals of selection in Italian brown trout

Background

Brown trout is one of the most widespread fresh-water fish species in Europe. The evolutionary history of and phylogenetic relationships between brown trout populations are complex, and this is especially true for Italian populations, which are heavily influenced in different ways by stocking practices. The characterization of the genetic structure of Italian brown trout populations may give information on the risk of losing endemic Italian populations due to lack of genetic diversity or to admixture with stocking populations. The identification of signatures of selection, and the information deriving from dense genotyping data will help genotype-informed breeding programs. We used a ddRAD-seq approach to obtain more than 100,000 single nucleotide polymorphisms (SNPs), and to characterize the population structure and signatures of selection in 90 brown trout samples.

Results

Italian brown trout populations are genetically differentiated, although the stocking practices have introduced strong admixture in endemic Italian trout, especially with the Atlantic lineage. Most of the analysed populations showed high levels of kinship and inbreeding. We detected putative signatures of selection using different approaches, and investigated if the regions were enriched for functional categories. Several regions putatively under selection and characterized by a reduction in heterozygosity across all the studied populations are enriched for genes involved in the response to viral infections.

Conclusions

Our results, which show evidence of admixture with the Atlantic lineage (commonly used for stocking), confirm the need for controlling stocking practices, in order to avoid the erosion of the endemic gene pool; given the apparently high levels of kinship and inbreeding in local populations, our results also show the need to take action for increasing gene diversity. In addition, we used the genetically-distinct lineages to detect signatures of selection and we identified putative signatures of selection in several regions associated with resistance to infectious diseases. These constitute candidate regions for the study of resistance to infections in wild and farmed trout.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12711-022-00698-7.

Footprints of local adaptation span hundreds of linked genes in the Atlantic silverside genome

The study of local adaptation in the presence of ongoing gene flow is the study of natural selection in action, revealing the functional genetic diversity most relevant to contemporary pressures. In addition to individual genes, genome‐wide architecture can itself evolve to enable adaptation. Distributed across a steep thermal gradient along the east coast of North America, Atlantic silversides (Menidia menidia) exhibit an extraordinary degree of local adaptation in a suite of traits, and the capacity for rapid adaptation from standing genetic variation, but we know little about the patterns of genomic variation across the species range that enable this remarkable adaptability. Here, we use low‐coverage, whole‐transcriptome sequencing of Atlantic silversides sampled along an environmental cline to show marked signatures of divergent selection across a gradient of neutral differentiation. Atlantic silversides sampled across 1371 km of the southern section of its distribution have very low genome‐wide differentiation (median F_ST = 0.006 across 1.9 million variants), consistent with historical connectivity and observations of recent migrants. Yet almost 14,000 single nucleotide polymorphisms (SNPs) are nearly fixed (F_ST > 0.95) for alternate alleles. Highly differentiated SNPs cluster into four tight linkage disequilibrium (LD) blocks that span hundreds of genes and several megabases. Variants in these LD blocks are disproportionately nonsynonymous and concentrated in genes enriched for multiple functions related to known adaptations in silversides, including variation in lipid storage, metabolic rate, and spawning behavior. Elevated levels of absolute divergence and demographic modeling suggest selection maintaining divergence across these blocks under gene flow. These findings represent an extreme case of heterogeneity in levels of differentiation across the genome, and highlight how gene flow shapes genomic architecture in continuous populations. Locally adapted alleles may be common features of populations distributed along environmental gradients, and will likely be key to conserving variation to enable future responses to environmental change.

Genetic variability and ontogeny predict microbiome structure in a disease-challenged montane amphibian

Amphibian populations worldwide are at risk of extinction from infectious diseases, including chytridiomycosis caused by the fungal pathogen Batrachochytrium dendrobatidis (Bd). Amphibian cutaneous microbiomes interact with Bd and can confer protective benefits to the host. The composition of the microbiome itself is influenced by many environment- and host-related factors. However, little is known about the interacting effects of host population structure, genetic variation and developmental stage on microbiome composition and Bd prevalence across multiple sites. Here we explore these questions in Amietia hymenopus, a disease-affected frog in southern Africa. We use microsatellite genotyping and 16S amplicon sequencing to show that the microbiome associated with tadpole mouthparts is structured spatially, and is influenced by host genotype and developmental stage. We observed strong genetic structure in host populations based on rivers and geographic distances, but this did not correspond to spatial patterns in microbiome composition. These results indicate that demographic and host genetic factors affect microbiome composition within sites, but different factors are responsible for host population structure and microbiome structure at the between-site level. Our results help to elucidate complex within- and among- population drivers of microbiome structure in amphibian populations. That there is a genetic basis to microbiome composition in amphibians could help to inform amphibian conservation efforts against infectious diseases.

Drift, not selection, shapes toll-like receptor variation among oceanic island populations

Understanding the relative role of different evolutionary forces in shaping the level and distribution of functional genetic diversity among natural populations is a key issue in evolutionary and conservation biology. To do so accurately genetic data must be analysed in conjunction with an unambiguous understanding of the historical processes that have acted upon the populations. Here, we focused on diversity at toll‐like receptor (TLR) loci, which play a key role in the vertebrate innate immune system and, therefore, are expected to be under pathogen‐mediated selection. We assessed TLR variation within and among 13 island populations (grouped into three archipelagos) of Berthelot's pipit, Anthus berthelotii, for which detailed population history has previously been ascertained. We also compared the variation observed with that found in its widespread sister species, the tawny pipit, Anthus campestris. We found strong evidence for positive selection at specific codons in TLR1LA, TLR3 and TLR4. Despite this, we found that at the allele frequency level, demographic history has played the major role in shaping patterns of TLR variation in Berthelot's pipit. Levels of diversity and differentiation within and across archipelagos at all TLR loci corresponded very closely with neutral microsatellite variation and with the severity of the bottlenecks that occurred during colonization. Our study shows that despite the importance of TLRs in combating pathogens, demography can be the main driver of immune gene variation within and across populations, resulting in patterns of functional variation that can persist over evolutionary timescales.

Isolation over 35 years in a heated biotest basin causes selection on MHC class IIß genes in the European perch (Perca fluviatilis L.)

Genes that play key roles in host immunity such as the major histocompatibility complex (MHC) in vertebrates are expected to be major targets of selection. It is well known that environmental conditions can have an effect on host–parasite interactions and may thus influence the selection on MHC. We analyzed MHC class IIß variability over 35 years in a population of perch (Perca fluviatilis) from the Baltic Sea that was split into two populations separated from each other. One population was subjected to heating from cooling water of a nuclear power plant and was isolated from the surrounding environment in an artificial lake, while the other population was not subjected to any change in water temperature (control). The isolated population experienced a change of the allelic composition and a decrease in allelic richness of MHC genes compared to the control population. The two most common MHC alleles showed cyclic patterns indicating ongoing parasite–host coevolution in both populations, but the alleles that showed a cyclic behavior differed between the two populations. No such patterns were observed at alleles from nine microsatellite loci, and no genetic differentiation was found between populations. We found no indications for a genetic bottleneck in the isolated population during the 35 years. Additionally, differences in parasitism of the current perch populations suggest that a change of the parasite communities has occurred over the isolation period, although the evidence in form of in-depth knowledge of the change of the parasite community over time is lacking. Our results are consistent with the hypothesis of a selective sweep imposed by a change in the parasite community.

Applications in the search for genomic selection signatures in fish

Selection signatures are genomic regions harboring DNA sequences functionally involved in the genetic variation of traits subject to selection. Selection signatures have been intensively studied in recent years because of their relevance to evolutionary biology and their potential association with genes that control phenotypes of interest in wild and domestic populations. Selection signature research in fish has been confined to a smaller scale, due in part to the relatively recent domestication of fish species and limited genomic resources such as molecular markers, genetic mapping, DNA sequences, and reference genomes. However, recent genomic technology advances are paving the way for more studies that may contribute to the knowledge of genomic regions underlying phenotypes of biological and productive interest in fish.

Immunogenetic factors affecting susceptibility of humans and rodents to hantaviruses and the clinical course of hantaviral disease in humans

We reviewed the associations of immunity-related genes with susceptibility of humans and rodents to hantaviruses, and with severity of hantaviral diseases in humans. Several class I and class II HLA haplotypes were linked with severe or benign hantavirus infections, and these haplotypes varied among localities and hantaviruses. The polymorphism of other immunity-related genes including the C4A gene and a high-producing genotype of TNF gene associated with severe PUUV infection. Additional genes that may contribute to disease or to PUUV infection severity include non-carriage of the interleukin-1 receptor antagonist (IL-1RA) allele 2 and IL-1β (-511) allele 2, polymorphisms of plasminogen activator inhibitor (PAI-1) and platelet GP1a. In addition, immunogenetic studies have been conducted to identify mechanisms that could be linked with the persistence/clearance of hantaviruses in reservoirs. Persistence was associated during experimental infections with an upregulation of anti-inflammatory responses. Using natural rodent population samples, polymorphisms and/or expression levels of several genes have been analyzed. These genes were selected based on the literature of rodent or human/hantavirus interactions (some Mhc class II genes, Tnf promoter, and genes encoding the proteins TLR4, TLR7, Mx2 and β3 integrin). The comparison of genetic differentiation estimated between bank vole populations sampled over Europe, at neutral and candidate genes, has allowed to evidence signatures of selection for Tnf, Mx2 and the Drb Mhc class II genes. Altogether, these results corroborated the hypothesis of an evolution of tolerance strategies in rodents. We finally discuss the importance of these results from the medical and epidemiological perspectives.

Genetic variation at the MHC in a population of introduced wild turkeys

Genetic variation in the major histocompatibility complex (MHC) is known to affect disease resistance in many species. Investigations of MHC diversity in populations of wild species have focused on the antigen presenting class IIβ molecules due to the known polymorphic nature of these genes and the role these molecules play in pathogen recognition. Studies of MHC haplotype variation in the turkey ( Meleagris gallopavo ) are limited. This study was designed to examine MHC diversity in a group of Eastern wild turkeys ( Meleagris gallopavo silvestris ) collected during population expansion following reintroduction of the species in southern Wisconsin, USA. Southern blotting with BG and class IIβ probes and single nucleotide polymorphism (SNP) genotyping was used to measure MHC variation. SNP analysis focused on single copy MHC genes flanking the highly polymorphic class IIβ genes. Southern blotting identified 27 class IIβ phenotypes, whereas SNP analysis identified 13 SNP haplotypes occurring in 28 combined genotypes. Results show that genetic diversity estimates based on RFLP (Southern blot) analysis underestimate the level of variation detected by SNP analysis. Sequence analysis of the mitochondrial D-loop identified 7 mitochondrial haplotypes (mitotypes) in the sampled birds. Results show that wild turkeys located in southern Wisconsin have a genetically diverse MHC and originate from several maternal lineages.

Parasite-mediated selection of major histocompatibility complex variability in wild brandt's voles (Lasiopodomys brandtii) from Inner Mongolia, China

Background

Genes of the major histocompatibility complex (MHC) exhibit high levels of variability, which is believed to have arisen through pathogen-mediated selection. We investigated the relationship between parasite load and genetic diversity at selectively neutral, non-coding markers (microsatellites) and adaptive genetic variation at a functionally important part of the MHC in six independent natural populations of Brandt’s voles (Lasiopodomys brandtii) from two regions of the Xilingol Grassland area of Inner Mongolia.

Results

Two-hundred and fifty-two voles were screened for gastrointestinal parasites, and were assessed for genetic variation. Parasite screening was done through non-invasive fecal egg counts, while allelic diversity was determined via single-stranded conformation polymorphism and DNA sequencing. We detected eight distinct helminth egg morphotypes. A total of 10 microsatellite loci were genotyped and 19 unique MHC class II B alleles were isolated. The rate of nonsynonymous substitutions (d_N) exceeded the rate of synonymous substitutions (d_S) at putative antigen binding sites of DRB. Neutral and adaptive genetic diversity differed between the six vole populations. To test the main pathogen-driven selection hypotheses for the maintenance of host MHC diversity and parasite species-specific co-evolutionary effects, multivariate approaches (generalized linear mixed models) were used to test for associations between the MHC class II DRB genotype and infections with nematodes. We found no evidence for heterozygote advantage, and overall heterozygosity was lower than expected in the MHC alleles. We identified an association between the parasite load and specific MHC alleles in the voles, and this pattern varied between geographic regions.

Conclusions

The results suggest that MHC variability in Brandt’s voles is maintained by rare allele advantage and fluctuating selection, but the data failed to show any heterozygote advantage effect. Our results add to a growing body of evidence showing that the mode and relative strength of pathogen-driven selection acting on MHC diversity varies within specific wild populations. In addition, our study contributes to the understanding of what maintains MHC diversity, of host-pathogen coevolution and of how genetic diversity is maintained in voles.

Inference of selection based on temporal genetic differentiation in the study of highly polymorphic multigene families

The co-evolutionary arms race between host immune genes and parasite virulence genes is known as Red Queen dynamics. Temporal fluctuations in allele frequencies, or the ‘turnover’ of alleles at immune genes, are concordant with predictions of the Red Queen hypothesis. Such observations are often taken as evidence of host-parasite co-evolution. Here, we use computer simulations of the Major Histocompatibility Complex (MHC) of guppies (Poecilia reticulata) to study the turnover rate of alleles (temporal genetic differentiation, G'_ST). Temporal fluctuations in MHC allele frequencies can be ≥≤order of magnitude larger than changes observed at neutral loci. Although such large fluctuations in the MHC are consistent with Red Queen dynamics, simulations show that other demographic and population genetic processes can account for this observation, these include: (1) overdominant selection, (2) fluctuating population size within a metapopulation, and (3) the number of novel MHC alleles introduced by immigrants when there are multiple duplicated genes. Synergy between these forces combined with migration rate and the effective population size can drive the rapid turnover in MHC alleles. We posit that rapid allelic turnover is an inherent property of highly polymorphic multigene families and that it cannot be taken as evidence of Red Queen dynamics. Furthermore, combining temporal samples in spatial F_ST outlier analysis may obscure the signal of selection.

Exploring neutral and adaptive processes in expanding populations of gilthead sea bream, Sparus aurata L., in the North-East Atlantic

Recent studies in empirical population genetics have highlighted the importance of taking into account both neutral and adaptive genetic variation in characterizing microevolutionary dynamics. Here, we explore the genetic population structure and the footprints of selection in four populations of the warm-temperate coastal fish, the gilthead sea bream (Sparus aurata), whose recent northward expansion has been linked to climate change. Samples were collected at four Atlantic locations, including Spain, Portugal, France and the South of Ireland, and genetically assayed using a suite of species-specific markers, including 15 putatively neutral microsatellites and 23 expressed sequence tag-linked markers, as well as a portion of the mitochondrial DNA (mtDNA) control region. Two of the putatively neutral markers, Bld-10 and Ad-10, bore signatures of strong directional selection, particularly in the newly established Irish population, although the potential 'surfing effect' of rare alleles at the edge of the expansion front was also considered. Analyses after the removal of these loci suggest low but significant population structure likely affected by some degree of gene flow counteracting random genetic drift. No signal of historic divergence was detected at mtDNA. BLAST searches conducted with all 38 markers used failed to identify specific genomic regions associated to adaptive functions. However, the availability of genomic resources for this commercially valuable species is rapidly increasing, bringing us closer to the understanding of the interplay between selective and neutral evolutionary forces, shaping population divergence of an expanding species in a heterogeneous milieu.

Contrasted patterns of selection on MHC-linked microsatellites in natural populations of the Malagasy plague reservoir

Plague (Yersinia pestis infection) is a highly virulent rodent disease that persists in many natural ecosystems. The black rat (Rattus rattus) is the main host involved in the plague focus of the central highlands of Madagascar. Black rat populations from this area are highly resistant to plague, whereas those from areas in which the disease is absent (low altitude zones of Madagascar) are susceptible. Various lines of evidence suggest a role for the Major Histocompatibility Complex (MHC) in plague resistance. We therefore used the MHC region as a candidate for detecting signatures of plague-mediated selection in Malagasy black rats, by comparing population genetic structures for five MHC-linked microsatellites and neutral markers in two sampling designs. We first compared four pairs of populations, each pair including one population from the plague focus and one from the disease-free zone. Plague-mediated selection was expected to result in greater genetic differentiation between the two zones than expected under neutrality and this was observed for one MHC-class I-linked locus (D20Img2). For this marker as well as for four other MHC-linked loci, a geographic pattern of genetic structure was found at local scale within the plague focus. This pattern would be expected if plague selection pressures were spatially variable. Finally, another MHC-class I-linked locus (D20Rat21) showed evidences of balancing selection, but it seems more likely that this selection would be related to unknown pathogens more widely distributed in Madagascar than plague.

Evidence for selection at cytokine loci in a natural population of field voles (Microtus agrestis)

Individuals in natural populations are frequently exposed to a wide range of pathogens. Given the diverse profile of gene products involved in responses to different types of pathogen, this potentially results in complex pathogen-specific selection pressures acting on a broad spectrum of immune system genes in wild animals. Thus far, studies into the evolution of immune genes in natural populations have focused almost exclusively on the Major Histocompatibility Complex (MHC). However, the MHC represents only a fraction of the immune system and there is a need to broaden research in wild species to include other immune genes. Here, we examine the evidence for natural selection in a range of non-MHC genes in a natural population of field voles (Microtus agrestis). We concentrate primarily on genes encoding cytokines, signalling molecules critical in eliciting and mediating immune responses and identify signatures of natural selection acting on several of these genes. In particular, genetic diversity within Interleukin 1 beta and Interleukin 2 appears to have been maintained through balancing selection. Taken together with previous findings that polymorphism within these genes is associated with variation in resistance to multiple pathogens, this suggests that pathogen-mediated selection may be an important force driving genetic diversity at cytokine loci in voles and other natural populations. These results also suggest that, along with the MHC, preservation of genetic variation within cytokine genes should be a priority for the conservation genetics of threatened wildlife populations.

Drift and selection influence geographic variation at immune loci of prairie-chickens

Previous studies of immunity in wild populations have focused primarily on genes of the major histocompatibility complex (MHC); however, studies of model species have identified additional immune-related genes that also affect fitness. In this study, we sequenced five non-MHC immune genes in six greater prairie-chicken (Tympanuchus cupido) populations that have experienced varying degrees of genetic drift as a consequence of population bottlenecks and fragmentation. We compared patterns of geographic variation at the immune genes with six neutral microsatellite markers to investigate the relative effects of selection and genetic drift. Global F(ST) outlier tests identified positive selection on just one of five immune genes (IAP-1) in one population. In contrast, at other immune genes, standardized G'(ST) values were lower than those at microsatellites for a majority of pairwise population comparisons, consistent with balancing selection or with species-wide positive or purifying selection resulting in similar haplotype frequencies across populations. The effects of genetic drift were also evident as summary statistics (e.g., Tajima's D) did not differ from neutrality for the majority of cases, and immune gene diversity (number of haplotypes per gene) was correlated positively with population size. In summary, we found that both genetic drift and selection shaped variation at the five immune genes, and the strength and type of selection varied among genes. Our results caution that neutral forces, such as drift, can make it difficult to detect current selection on genes.

Evidence of neutral and adaptive genetic divergence between European trout populations sampled along altitudinal gradients

Species with a wide geographical distribution are often composed of distinct subgroups which may be adapted to their local environment. European trout (Salmo trutta species complex) provide an example of such a complex consisting of several genetically and ecologically distinct forms. However, trout populations are strongly influenced by human activities, and it is unclear to what extent neutral and adaptive genetic differences have persisted. We sampled 30 Swiss trout populations from heterogeneous environments along replicated altitudinal gradients in three major European drainages. More than 850 individuals were genotyped at 18 microsatellite loci which included loci diagnostic for evolutionary lineages and candidate markers associated with temperature tolerance, reproductive timing and immune defence. We find that the phylogeographic structure of Swiss trout populations has not been completely erased by stocking. Distinct genetic clusters corresponding to the different drainages could be identified, although nonindigenous alleles were clearly present, especially in the two Mediterranean drainages. We also still detected neutral genetic differentiation within rivers which was often associated with the geographical distance between populations. Five loci showed evidence of divergent selection between populations with several drainage-specific patterns. Lineage-diagnostic markers, a marker linked to a quantitative trait locus for upper temperature tolerance in other salmonids and a marker linked to the major histocompatibility class I gene were implicated in local adaptation and some patterns were associated with altitude. In contrast, tentative evidence suggests a signal of balancing selection at a second immune relevant gene (TAP2). Our results confirm the persistence of both neutral and potentially adaptive genetic differences between trout populations in the face of massive human-mediated dispersal.

Extent and scale of local adaptation in salmonid fishes: review and meta-analysis

What is the extent and scale of local adaptation (LA)? How quickly does LA arise? And what is its underlying molecular basis? Our review and meta-analysis on salmonid fishes estimates the frequency of LA to be ∼55-70%, with local populations having a 1.2 times average fitness advantage relative to foreign populations or to their performance in new environments. Salmonid LA is evident at a variety of spatial scales (for example, few km to>1000 km) and can manifest itself quickly (6-30 generations). As the geographic scale between populations increases, LA is generally more frequent and stronger. Yet the extent of LA in salmonids does not appear to differ from that in other assessed taxa. Moreover, the frequency with which foreign salmonid populations outperform local populations (∼23-35%) suggests that drift, gene flow and plasticity often limit or mediate LA. The relatively few studies based on candidate gene and genomewide analyses have identified footprints of selection at both small and large geographical scales, likely reflecting the specific functional properties of loci and the associated selection regimes (for example, local niche partitioning, pathogens, parasites, photoperiodicity and seasonal timing). The molecular basis of LA in salmonids is still largely unknown, but differential expression at the same few genes is implicated in the convergent evolution of certain phenotypes. Collectively, future research will benefit from an integration of classical and molecular approaches to understand: (i) species differences and how they originate, (ii) variation in adaptation across scales, life stages, population sizes and environmental gradients, and (iii) evolutionary responses to human activities.

Identifying footprints of selection in stocked brown trout populations: a spatio-temporal approach

Studies of interactions between farmed and wild salmonid fishes have suggested reduced fitness of farmed strains in the wild, but evidence for selection at the genic level is lacking. We studied three brown trout populations in Denmark which have been significantly admixed with stocked hatchery trout (19-64%), along with two hatchery strains used for stocking. The wild populations were represented by contemporary samples (2000-2006) and two of them by historical samples (1943-1956). We analysed 61 microsatellite loci, nine of which showed putative functional relationships [expressed sequence tag (EST)-linked or quantitative trait loci]. F(ST)-based outlier tests provided support for diversifying selection at chromosome regions marked by three loci, two anonymous and one EST-linked. Patterns of differentiation suggested that the loci were candidates for being under diversifying hitch-hiking selection in hatchery vs. wild environments. Analysis of hatchery strain admixture proportions showed that in one wild population, two of the loci showed significantly lower admixture proportions than the putatively neutral loci, implying contemporary selection against alleles introduced by hatchery strain trout. In the most strongly admixed population, however, there was no evidence for selection, possibly because of immigration by stocked trout overcoming selection against hatchery-derived alleles or supportive breeding practices allowing hatchery strain trout to escape natural selection. To our knowledge, this is the first study demonstrating footprints of selection in wild salmonid populations subject to spawning intrusion by farmed fish.

Beyond MHC: signals of elevated selection pressure on Atlantic salmon (Salmo salar) immune-relevant loci

Using Atlantic salmon (Salmo salar) as a model system, we investigated whether 18 microsatellites tightly linked to immune-relevant genes have experienced different selection pressures than 76 loci with no obvious association with immune function. Immune-relevant loci were identified as outliers by two outlier tests significantly more often than nonimmune linked loci (22% vs. 1.6%). In addition, the allele frequencies of immune relevant markers were more often correlated with latitude and temperature. Combined, these results support the hypothesis that immune-relevant loci more frequently exhibit footprints of selection than other loci. They also indicate that the correlation between immune-relevant loci and latitude may be due to temperature-induced differences in pathogen-driven selection or some other environmental factor correlated with latitude.

Temporal variation at the MHC class IIb in wild populations of the guppy (Poecilia reticulata)

Understanding genetic diversity in natural populations is a fundamental objective of evolutionary biology. The immune genes of the major histocompatibility complex (MHC) are excellent candidates to study such diversity because they are highly polymorphic in populations. Although balancing selection may be responsible for maintaining diversity at these functionally important loci, temporal variation in selection pressure has rarely been examined. We examine temporal variation in MHC class IIB diversity in nine guppy (Poecilia reticulata) populations over two years. We found that five of the populations changed significantly more at the MHC than at neutral (microsatellite) loci as measured by F(ST), which suggests that the change at the MHC was due to selection and not neutral processes. Additionally, pairwise population differentiation measures at the MHC were higher in 2007 than in 2006, with the signature of selection changing from homogenizing to diversifying selection or neutral evolution. Interestingly, within the populations the magnitude of the change at the MHC between years was related to the change in the proportion of individuals infected by a common parasite, indicating a link between genetic structure and the parasite. Our data thereby implicate temporal variation in selective pressure as an important mechanism maintaining diversity at the MHC in wild populations.

Evidence for directional selection at a novel major histocompatibility class I marker in wild common frogs (Rana temporaria) exposed to a viral pathogen (Ranavirus)

Whilst the Major Histocompatibility Complex (MHC) is well characterized in the anuran Xenopus, this region has not previously been studied in another popular model species, the common frog (Rana temporaria). Nor, to date, have there been any studies of MHC in wild amphibian host-pathogen systems. We characterise an MHC class I locus in the common frog, and present primers to amplify both the whole region, and specifically the antigen binding region. As no more than two expressed haplotypes were found in over 400 clones from 66 individuals, it is likely that there is a single class I locus in this species. This finding is consistent with the single class I locus in Xenopus, but contrasts with the multiple loci identified in axolotls, providing evidence that the diversification of MHC class I into multiple loci likely occurred after the Caudata/Anura divergence (approximately 350 million years ago) but before the Ranidae/Pipidae divergence (approximately 230 mya). We use this locus to compare wild populations of common frogs that have been infected with a viral pathogen (Ranavirus) with those that have no history of infection. We demonstrate that certain MHC supertypes are associated with infection status (even after accounting for shared ancestry), and that the diseased populations have more similar supertype frequencies (lower F_ST) than the uninfected. These patterns were not seen in a suite of putatively neutral microsatellite loci. We interpret this pattern at the MHC locus to indicate that the disease has imposed selection for particular haplotypes, and hence that common frogs may be adapting to the presence of Ranavirus, which currently kills tens of thousands of amphibians in the UK each year.

Balancing and directional selection at exon-2 of the MHC DQB1 locus among populations of odontocete cetaceans

The diversity of exon-2 (peptide-binding region) of the DQB1 locus (Class II, major histocompatibility complex, MHC) was investigated on an extended sample of populations of three focal cetacean species (two sibling delphinid species and another in the same family). We tested the hypothesis that dolphin populations with a worldwide distribution across different habitats and geographic regions will be under differential selective pressure by comparing DQB1 variation with variation at neutral markers and by investigating putative functional residues within the exon-2 sequence at the population level. Variation at the DQB1 locus was not correlated to neutral differentiation (assessed by comparison with microsatellite DNA markers), and overall F(ST) values were significantly lower for the MHC locus, consistent with expectations for balancing selection. Measures of heterozygosity and d(n)/d(s) ratios were also consistent with balancing selection. However, outliers in the F(ST) comparisons and the analysis of putative functional residues suggested incidences of directional selection in local populations.

Balancing selection at a frog antimicrobial peptide locus: fluctuating immune effector alleles?

Balancing selection is common on many defense genes, but it has rarely been reported for immune effector proteins such as antimicrobial peptides (AMPs). We describe genetic diversity at a brevinin-1 AMP locus in three species of leopard frogs (Rana pipiens, Rana blairi, and Rana palustris). Several highly divergent allelic lineages are segregating at this locus. That this unusual pattern results from balancing selection is demonstrated by multiple lines of evidence, including a ratio of nonsynonymous/synonymous polymorphism significantly higher than 1, the ZnS test, incongruence between the number of segregating sites and haplotype diversity, and significant Tajima's D values. Our data are more consistent with a model of fluctuating selection in which alleles change frequencies over time than with a model of stable balancing selection such as overdominance. Evidence for fluctuating selection includes skewed allele frequencies, low levels of synonymous variation, nonneutral values of Tajima's D within allelic lineages, an inverse relationship between the frequency of an allelic lineage and its degree of polymorphism, and divergent allele frequencies among populations. AMP loci could be important sites of adaptive genetic diversity, with consequences for host–pathogen coevolution and the ability of species to resist disease epidemics.