MHC standing genetic variation and pathogen resistance in wild Atlantic salmon

Mélanie Dionne, Kristina M. Miller, Julian J. Dodson, Louis Bernatchez


Pathogens are increasingly emerging in human-altered environments as a serious threat to biodiversity. In this context of rapid environmental changes, improving our knowledge on the interaction between ecology and evolution is critical. The objective of this study was to evaluate the influence of an immunocompetence gene, the major histocompatibility complex (MHC) class IIβ, on the pathogen infection levels in wild Atlantic salmon populations, Salmo salar, and identify selective agents involved in contemporary coevolution. MHC variability and bacterial infection rate were determined throughout the summer in juvenile salmon from six rivers belonging to different genetic and ecological regions in Québec, Canada. A total of 13 different pathogens were identified in kidney by DNA sequence analysis, including a predominant myxozoa, most probably recently introduced in North America. Infection rates were the highest in southern rivers at the beginning of the summer (average 47.6±6.3% infected fish). One MHC allele conferred a 2.9 times greater chance of being resistant to myxozoa, while another allele increased susceptibility by 3.4 times. The decrease in frequency of the susceptibility allele but not other MHC or microsatellite alleles during summer was suggestive of a mortality event from myxozoa infection. These results supported the hypothesis of pathogen-driven selection in the wild by means of frequency-dependent selection or change in selection through time and space rather than heterozygous advantage, and underline the importance of MHC standing genetic variation for facing pathogens in a changing environment.

1. Introduction

Recent evidence suggests that parasites and infectious diseases are increasingly emerging in wild populations (Harvell et al. 1999; Daszak et al. 2000; Dobson & Foufopoulos 2001). Emerging diseases are due, in part, to increased human activity and most particularly to pollution, environmental degradation, climate change, habitat fragmentation, animal transfer and hatcheries (Harvell et al. 1999; Daszak et al. 2000; Dobson & Foufopoulos 2001). In the aquatic environment and in fish populations, in particular, several diseases such as whirling disease, furunculosis and infectious salmon anaemia have been associated with aquaculture-related activities (Dobson & Foufopoulos 2001). For instance, louse infestations associated with salmon farms are blamed for depressed wild pink salmon populations, potentially causing over 80 per cent mortality in some cases (Krkošek et al. 2007). These emerging pathogens have the potential to significantly reduce the abundance of threatened species (Dobson & Foufopoulos 2001; Lafferty et al. 2004) and as such, represent a serious threat to which wild populations have to adapt rapidly. In such a context of rapid environmental changes, linking ecological and evolutionary processes then becomes critical to evaluate the influence of standing genetic variation on contemporary adaptation (see also Porlier et al. 2009; Zheng et al. 2009).

Adaptation can occur over a very short time scale (Thompson 1998), suggesting that wild populations could potentially adapt to the increasing number of novel pathogens emerging in the environment. In vertebrates, genes of the major histocompatibility complex (MHC) are involved in pathogen resistance by encoding cell-surface proteins that bind peptide fragments derived from pathogens and present them to T cells that activate a specific immune response (Potts & Wakeland 1990). These genes are among the most polymorphic in vertebrates and most of the variability is concentrated at the peptide-binding region (PBR) involved in pathogen binding (Hughes & Yeager 1998). Theory predicts that, at the population level, high MHC diversity should confer resistance to a diverse array of pathogens (reviewed in Bernatchez & Landry 2003; Sommer 2005; Piertney & Oliver 2006). Recent studies conducted at a large spatial scale supported a positive association between MHC diversity and pathogen diversity in the wild, independent of historical colonization processes, in humans (Prugnolle et al. 2005) and Atlantic salmon (Dionne et al. 2007). Yet, the exact underlying mechanism involved in pathogen resistance at the individual level and how this translates into the maintenance of MHC diversity at the population level is not well understood (Sommer 2005; Piertney & Oliver 2006).

At the population level, MHC diversity is thought to be maintained by balancing selection through (i) heterozygote advantage, (ii) frequency-dependent selection, or (iii) variable selection in time and space (Nei & Hughes 1991; Hedrick 2002). The heterozygote advantage (or overdominance) hypothesis suggests that the fitness of heterozygote individuals is greater than that of both homozygotes because of their ability to present a broader spectrum of antigens to T cells (Hughes & Nei 1988). MHC heterozygote advantage was observed in humans for resistance to hepatitis B and HIV infections (Thursz et al. 1997; Carrington et al. 1999) and in mice for resistance to Salmonella infections (Penn et al. 2002; McClelland et al. 2003). On the other hand, frequency-dependent selection occurs when the advantage of possessing an MHC allele depends on its frequency in the population and is usually referred to as negative frequency-dependent selection or rare allele advantage (Takahata & Nei 1990). The variable selection in time and space hypothesis reasons that MHC polymorphism is maintained by shifting pathogen compositions that occur through time and space (Hedrick 2002). Support for the latter two hypotheses is provided by the association of certain MHC alleles with resistance to specific diseases, such as those found for malaria and HIV in humans (Hill et al. 1991; Carrington et al. 1999), Marek's disease in chicken (Briles et al. 1977) and furunculosis, infectious salmon anaemia and infectious haemopoietic necrosis virus in salmonids (Langefors et al. 2001; Grimholt et al. 2003; Miller et al. 2004). Further support for the latter hypothesis was also found in a temporal study conducted on the great reed warbler, Acrocephalus arundinaceus, where change in frequency of two MHC class I alleles contrasted with the stochastic variability of microsatellite markers (Westerdahl et al. 2004). When contrasting the heterozygote advantage and the frequency-dependent selection hypotheses in pathogen resistance, the latter has received more support in the literature (Sommer 2005). However, the exact pathogen resistance mechanism occurring at the individual level in the wild is yet to be clarified. Moreover, although the importance of pre-existing genetic variation, referred to as ‘standing genetic variation’, has been emphasized for adaptation and rapid evolution (Barrett & Schulter 2008), its importance for pathogen resistance in the wild is yet to be evaluated. Finally, although associations between pathogenic infection and population abundance have been observed in the wild (Krkošek et al. 2007; Miller & Vincent 2008), investigating demographic consequences of immune maladaptation is still in its infancy.

Anadromous populations of Atlantic salmon, Salmo salar, reproduce and spend the first years of their life in rivers along the North American and European Atlantic coasts, where a major decline in abundance has been observed during the last decades (Caron et al. 2005). After a feeding period of 1–3 years at sea, the majority of Atlantic salmon returns to spawn in natal rivers (Stabell 1984) and evidence suggests that they adapt to local conditions (reviewed in Taylor 1991; García de Leániz et al. 2007). Recently, adaptation to local temperature regime and bacterial community prevailing in the natal river has been supported through large-scale variations in MHC class IIβ diversity specifically at the PBR, suggesting the influence of pathogen-driven balancing selection on MHC diversity in Atlantic salmon (Dionne et al. 2007). Also, indirect evidence was found for lower fitness of immigrants coming from genetically and thermally distinct regions, suggesting a potential influence of local adaptation in limiting gene flow (Dionne et al. 2008). Evidence of adaptation through disassortative mating at the PBR in MHC class II was also found, possibly providing better pathogen resistance for the progeny (Landry et al. 2001). Evidence of selection at MHC class II was also found at the juvenile stage. Indeed, observed genotype frequencies of hatchery juveniles introduced in a river were significantly different after six months than the expected frequencies from parental crosses (de Eyto et al. 2007). However, uncertainty still persists whether early juvenile stages are able to mount specific immune responses (Bakke & Harris 1998). More importantly, information about selective agents, namely pathogens infecting Atlantic salmon in the wild along with the associated infection rates, is lacking. This lack of knowledge is even more pronounced for juvenile stages, as most pathogen resistance studies have been conducted on adults in laboratory conditions (Langefors et al. 2001; Grimholt et al. 2003; Miller et al. 2004).

The objectives of this study were, first, to document infection rates and identify pathogens infecting juvenile Atlantic salmon inhabiting rivers from distinct genetic groups and ecological regions; second to evaluate the influence of MHC class IIβ alleles in conferring resistance or susceptibility to pathogen infection in the wild; and third, to test which of the heterozygote advantage, frequency-dependent selection or variable selection in time and space hypotheses are prevailing at the individual level. This was achieved by testing whether (i) heterozygotes at MHC class II are less infected than homozygotes and whether (ii) there is an association between specific MHC alleles and the resistance or susceptibility of juveniles to the most prevalent infections. A final objective consisted of an indirect test for pathogen-induced mortality by following MHC allele frequency changes during the summer for alleles found to be associated with the most prevalent infections in the wild. In order to dissociate the influence of selection from other evolutionary forces, MHC variability was compared with that observed at putative neutral microsatellite loci.

2. Material and methods

(a) Sampling

Sampling was conducted in six rivers from Québec, Canada: the Trinité (TR), Laval (LA), and Sainte-Marguerite (SM) rivers on the north shore, and the Causapscal (CA), Dartmouth (DA), and Saint-Jean (SJ) rivers on the south shore of the St Lawrence Estuary (figure 1). These rivers belong to three different Atlantic salmon genetic regions according to a Bayesian clustering analysis based on neutral microsatellite markers: TR and LA in region 4; SM in region 5; and all southern rivers in region 6 (Dionne et al. 2008). These genetic regions are also ecologically differentiated, mainly in terms of their historical summer temperature index (average temperature regime for multiple rivers in regions 4, 5 and 6, respectively: 1066±132; 1542±160; and 1343±77 degree days; Dionne et al. 2008). Rivers draining the north shore of the St Lawrence Estuary are characterized by lower bacterial diversity in the water than those draining the south shore (2.5±2.4 and 38.9±12.8 ×103 intergrated optical density (IOD) mm−2, respectively, see details in Dionne et al. 2007). Each river was sampled on three occasions (June, July and August) during summer 2005. On each occasion, we sampled on average 50 parr per river (age 1+ and 2+ in an approximate ratio of 1 : 1 per sample), for a total of 898 juvenile salmon collected. Fish were handled using gloves to avoid external contamination. Each fish was individually inserted in a sterile sampling bag (Whirl-Pak, Nasco, Fort Atkinson, WI, USA), measured for fork length and quick-frozen on dry ice in order to preserve tissues for DNA and pathogen analyses. In the laboratory and under sterile conditions, the adipose fin was clipped and stored in 95 per cent ethanol for subsequent DNA analysis while the whole fish was kept at −80°C.

Figure 1

Location of the six Atlantic salmon rivers where sampling took place in Québec, Canada. The Sainte-Marguerite, Laval and Trinité rivers are located on the north shore while the Causapscal, Saint-Jean and Dartmouth rivers are located on the south shore of the St Lawrence Estuary.

(b) Microsatellite and MHC class IIβ genotyping

DNA was extracted from fin clips using the Qiagen DNeasy Tissue Kit (Mississauga, Ontario, Canada) following the guidelines of the manufacturer. Microsatellite polymorphism was quantified at 13 loci as detailed in the appendix of Dionne et al. (2007): Ssa85; Ssa202; Ssa197 (O'Reilly et al. 1996); Ssosl417 (Slettan et al. 1995); SsaD85 (T. King, unpublished data); SsaD71; SsaD144 (King et al. 2005); MST-3 (Presa & Guyomard 1996); Sssp1605; Sssp2210; Sssp2215; Sssp2216; and SsspG7 (Paterson et al. 2004). MHC class IIβ exon 2 was amplified using fluorescently labelled primers (Dionne et al. 2007) and alleles were separated via denaturing gradient gel electrophoresis (DGGE), as detailed in the appendix of Dionne et al. (2007).

(c) Bacterial infection analyses

Kidney tissue was selected to evaluate bacterial pathogenic infection in juvenile salmon as this organ is affected by multiple well-known fish diseases such as furunculosis, bacterial kidney disease, proliferative kidney disease and vibriosis (Bakke & Harris 1998). Also, a healthy kidney is normally free of bacteria, so the presence of bacteria implies pathogenic infection (Uhland et al. 2000). A total of 421 age 2+ parr were selected for infection analyses. The kidney of each fish was carefully dissected under a fume hood and under sterile conditions, with approximately 25 mg of tissue used for DNA extraction. Total DNA was extracted from the kidney using the QIAamp DNA Mini Kit following the guidelines of the manufacturer for the bacteria DNA protocol (Qiagen). Bacterial DNA was then specifically amplified using 16S rDNA universal primers (Crump et al. 2003) in a nested PCR protocol, as detailed in the analysis of bacterial diversity in water samples (Appendix 2 in Dionne et al. 2007). The 16S rDNA is a conserved region that is widely used to amplify and identify all bacteria present in a given environment (Woo et al. 2008). Distinct amplified bacterial sequences were separated using DGGE (BioRad, Hercules, CA, USA) under the conditions outlined in Appendix 2 of Dionne et al. (2007). Eleven different 16S rDNA fragments, amplified from 11 different bacteria species/strains, were included as standards on each DGGE gel to allow band comparison based on the positions of the standard sequences. Negative and positive controls (infected salmon with a known bacterium) were included in each PCR and were run on each gel to detect and avoid potential contamination and ensure proper amplification, which proved to be adequate. Each bacterial band was cut from the gel and sequenced using the 16S-357F and 16S-519R primers and the Big Dye Terminator cycle sequencing kit (Applied Biosystems, Foster City, CA, USA). Sequencing reactions were purified using the DyeEx 2.0 Spin Kit (Qiagen) and electrophoresis was conducted on a 3730 Capillary DNA Sequencer (Applied Biosystems). Nucleotide sequences were edited in BioEdit (Hall 1999) and identified using the Basic Local Alignment Search Tool in the National Center for Biotechnology Information.

(d) Statistical analyses

(i) Microsatellite basic analyses

The potential occurrence of null alleles and scoring errors due to stuttering or large allele dropout in the dataset was assessed using the software Micro-Checker (van Oosterhout et al. 2004). Deviation from Hardy–Weinberg expectations was tested for each locus and river sampled, and linkage disequilibrium was tested for each pair of loci using Genepop 3.4 (Raymond & Rousset 1995). The LnRH neutrality test was performed on each microsatellite locus to verify the neutrality of each marker (Schlötterer 2002; see Dionne et al. (2007) for calculation details). Temporal stability of microsatellite allele frequency distributions per river across the three sampling periods was assessed using the Fisher's exact probability test in Genepop 3.4 (Raymond & Rousset 1995).

(ii) Infection analyses

In order to evaluate the overall pattern of variation in infection rate of wild juvenile salmon between regions and sampling periods, a two-factor logistic regression was conducted on infection status of the fish (infected or uninfected) using SAS 9.1 (SAS Institute). A second general analysis, a two-factor Poisson regression analysis, was conducted using SAS to evaluate the regional and temporal variations in pathogen diversity (number of different pathogen strains found in each fish). It was not possible to include individual rivers as a fixed factor in these analyses because of a lack of power due to the limited number of infected fish in certain rivers.

(iii) Infection and MHC polymorphism

Levels of heterozygosity were calculated for infected and uninfected fish using Genepop 3.4 (Raymond & Rousset 1995) and compared between MHC and microsatellite markers. Under the heterozygote advantage hypothesis, heterozygosity is expected to be lower for infected compared with uninfected fish at MHC but not at microsatellite markers. To test this hypothesis further, MHC p-distance for each fish was calculated using Mega 3.1 (Kumar et al. 2000) for both the overall polymorphic sites and for the PBR only. This p-distance represents a genetic measure of allelic divergence within individuals, which incorporates the proportion of different amino acids between two MHC alleles of any given fish (Nei & Kumar 2000). The PBR was identified based on codons involved in pathogen binding in humans (Brown et al. 1988, 1993) and corresponded to those presented in Landry & Bernatchez (2001) and Dionne et al. (2007). For comparison, a comparable individual distance measure for microsatellites, representing the difference in mutational steps, was estimated by calculating mean difference in the number of repeats between two microsatellite alleles per locus for each individual (hereafter referred to as ‘repeat distance’). A two-factor multivariate analysis of variance (MANOVA) was then conducted using SAS to evaluate the potential implication of the MHC p-distance (overall and PBR), the microsatellite repeat distance and the fork length in differentiating between infected and uninfected salmon among the different sampling periods. Fork length was used here as a phenotypic index representing general body condition, which might influence indirectly the probability for a fish to be infected in the wild.

To evaluate the potential association between specific MHC alleles and either resistance or susceptibility to the most prevalent infection, a Fisher's exact probability test was first conducted using SAS and odds ratios were calculated for each allele. A high odds ratio would indicate a high probability for a fish to be infected when having a given allele. Secondly, a multivariate logistic regression was conducted on fish infection status (infected or uninfected with the most prevalent infection), with region, sampling period and all MHC alleles as factors integrated in a stepwise process using SAS. The Fisher's exact probability test and the multivariate logistic regression were also conducted at the amino acid level to identify specific codon positions potentially associated with resistance or susceptibility to the most prevalent infection. Finally, frequency of the resistance and susceptibility alleles identified in the previous analyses was compared among sampling periods, along with the other MHC and microsatellite alleles, in order to test for differential mortality caused by the predominant infection during summer. Mantel–Haenszel Χ2-tests were conducted on all juveniles sampled to evaluate trends in allele frequency changes using SAS and sequential Bonferroni correction was applied to evaluate significance (Rice 1989). All analyses in SAS were conducted in SAS software 9.1 (SAS Institute).

3. Results

(a) Microsatellite polymorphism

No evidence for null alleles or scoring errors due to stuttering or large allele dropout was found in the overall dataset. Deviations from Hardy–Weinberg equilibrium were identified in 13 out of 234 comparisons (all heterozygote deficits), a number close to the 12 significant tests expected by chance at α=0.05. There was no evidence of significant deviations associated with either particular loci or populations, which translated into small global FIS values ranging from −0.010 to 0.018 depending on the population. No significant linkage disequilibrium was detected within populations or between locus pairs across populations. The LnRH neutrality test identified microsatellite Sssp1605 as a locus potentially under selection, which was thus removed from further analyses (Schlötterer 2002, see Dionne et al. (2007) for calculation details).

(b) Infection types and prevalence

Overall, 12.1 per cent of the juvenile salmon analysed were infected with one or more pathogens (51/421 salmon). The observed infection rates in rivers of both north and south shores were highest in June and declined rapidly thereafter according to the logistic regression and the post hoc comparisons (p<0.0001 between June and the two other summer periods; table 1). Indeed, a proportion of 20.6±6.9 and 47.6±6.3 per cent of juvenile salmon were found to be infected at the beginning of the summer in north and south shore rivers, respectively, while infection rates dropped, respectively, to 3.2±1.9 and 2.3±1.6 per cent by the end of the summer (table 2). There was also a nearly significant trend for infection rates to be higher in St Lawrence south shore than in north shore rivers (p=0.058). Pathogen diversity, in terms of number of different pathogen species per individual, also followed the same pattern with higher pathogen diversity in June (p<0.0001) and a tendency for higher pathogen diversity in south shore when compared with the north shore rivers, although the latter was not significant (p=0.23). Note that values of pathogen diversity are close to zero due to the high number of uninfected fish. Sex ratio was slightly male biased in infected salmon (male/female=1.27) compared with uninfected individuals (male/female=1.09), although both were not departing significantly from equal sex ratios (Χ2-tests: p=0.40 and 0.39 for infected and uninfected groups).

View this table:
Table 1

Logistic regression analysis on infection status (infected or uninfected) of juvenile Atlantic salmon and Poisson regression analysis on pathogen diversity (number of different pathogens) infecting each individual according to the region (north or south shore) and the sampling period during summer (June, July or August).

View this table:
Table 2

Mean infection rate and pathogen diversity infecting juvenile Atlantic salmon from north and south shore rivers of the St Lawrence Estuary. (s.d., standard deviation.)

A total of 13 different pathogens were identified based on their DNA sequences of 151 bp on average, as indicated by the best match with known sequences from GenBank (table 3). The nearest taxonomic neighbour represented the genus (or order for myxozoa) corresponding to all best matches in GenBank with the unknown sequence (homology above 80%) and the proportion of sequence similarity represented highest homology with species in that genus group. A total of 92.2 per cent of the infected juvenile salmon had a single pathogen in the kidney (47/51 juvenile salmon). Bacterial kidney infections were present in 29.4 per cent of the infected fish (15/51 juvenile salmon) and 3.6 per cent of the salmon sampled (15/421 juvenile salmon). Each identified bacteria had low prevalence in the populations studied, infecting on an average 0.36±0.81 per cent juvenile salmon per river over all time periods (range: 0.0–4.4%, table 3). However, one predominant pathogen was present in 76.5 per cent of the infected fish (39/51 juvenile salmon) and was found to infect 9.3 per cent of the juvenile salmon sampled (39/421 juvenile salmon). This represented an average prevalence of 8.9±6.4 per cent infected juvenile salmon per river during summer (range: 0.0–16.4%), with higher infection rates in the St Lawrence south shore than in north shore rivers (south shore prevalence=14.2±2.3%, north shore prevalence=3.7±3.9%; t-test, p=0.011). This predominant pathogen was identified as a myxozoa, class Myxosporea, order Bivalvulidae based on GenBank known sequences (85% similarity with available species in the Bivalvulidae order). This pathogen was also compared with unpublished myxozoa sequences available from C. Whipps (2008; SUNY College of Environmental Science & Forestry, Syracuse, NY) and fell within the Myxobolus clade as defined by Fiala (2006), yielding matches (78% similarity) to Myxobolus species infecting salmonids, i.e. M. arcticus and M. kisutchi (C. Whipps 2008, personal communication). Histological observations of infected kidneys yielded no myxospore to allow further identification (C. Whipps 2008, personal communication).

View this table:
Table 3

Pathogen identity found in the kidney of juvenile Atlantic salmon from six rivers in Québec: Sainte-Marguerite (SM), Laval (LA), Trinité (TR), Causapscal (CA), Dartmouth (DA) and Saint-Jean (SJ). (GenBank sequence accession number for the myxozoa and each bacterium is indicated. Nearest taxonomic neighbour (genus or order for myxozoa) and proportion of sequence similarity with GenBank available species for each associated genus is indicated. Prevalence indicates the proportion of salmon infected per river over all time periods during summer.)

(c) Infection and MHC polymorphism

A total of 13 different MHC class IIβ alleles were found among all sampled juvenile salmon (figure 2). The number of different amino acids between MHC alleles averaged 9.3±2.7 (range: 4–17), while the proportion of different amino acids averaged 0.12±0.04 (range: 0.03–0.22). The observed heterozygosity at microsatellite loci was 0.83 and 0.86 for infected and uninfected fish, respectively, and both groups exhibited a slight yet significant deficit in heterozygotes compared with the expected values (HE: 0.86 and 0.88, respectively; p<0.0001 for both groups). A deficit in heterozygotes was expected as each group represented a pool of mixed individuals from different populations (Wahlund effect; all rivers pooled because of low number of infected fish). A deficit in heterozygotes was also observed at MHC class IIβ for uninfected fish (HO: 0.84, HE: 0.90, p<0.0001). However, in infected fish, an excess in heterozygotes was observed at MHC class IIβ (HO: 0.96, HE: 0.90, p=0.048), contrasting with neutral expectations. Overall, the observed microsatellite heterozygosity levels were similar for infected and uninfected juvenile salmon (0.83±0.08 and 0.86±0.08, respectively; t-test, p=0.45). However, MHC heterozygosity was higher for infected than for uninfected salmon (0.96 and 0.84, respectively). The MANOVA revealed an overall significant difference in MHC p-distance, microsatellite repeat distance or fork length between infected and uninfected juvenile salmon (p=0.027) and across the different time periods (p=0.029), but there was no significant interaction between infection status and time period (p=0.20). In the univariate tests, MHC p-distance, calculated both using all polymorphic sites or PBR only, was higher for infected than for uninfected juvenile salmon, although this trend was only significant within the PBR (p=0.068 and 0.019, respectively; table 4). On the other hand, no significant difference in microsatellite repeat distance was observed between infected and uninfected juvenile salmon (p=0.829). Univariate tests also revealed that fork length was the factor varying across time, with juvenile salmon becoming longer with time during summer as expected following summer growth (p=0.008). There was also a non-significant tendency for infected juvenile salmon to be smaller than the uninfected salmon (p=0.07).

Figure 2

Allele frequency of each MHC class IIβ allele found in infected (black bars) and uninfected (grey bars) juvenile Atlantic salmon sampled in six rivers in the province of Québec, Canada. (Asterisks indicate alleles significantly associated with infection probability according to the Fisher's exact probability test.)

View this table:
Table 4

Two-factor MANOVA comparing individual genetic distances (MHC p-distance and microsatellite repeat distance) and fork lengths between the different infection status of juvenile Atlantic salmon (infected or uninfected) and among the different sampling periods during summer (June, July and August). (MHC p-distance was calculated either on all polymorphic sites or on PBR only. The p-values indicate significance in the univariate tests for each dependent variable compared between infected and uninfected fish (infection status) and among sampling periods (period). s.d., standard deviation. Values in bold represent significant tests.)

An association was uncovered between the occurrence of specific MHC alleles and the probability of juvenile salmon to be infected by myxozoa, the most prevalent pathogen, according to the Fisher's exact probability test (table 5). Allele 7 had a low odds ratio and was therefore identified as a potential resistance allele decreasing the risk of myxozoa infection (p=0.02, odds ratio=0.28). Conversely, allele 8 had a high odds ratio and represented a potential susceptibility allele increasing the chance of myxozoa infection (p=0.03, odds ratio=2.57). In order to take into account multiple testing in this analysis, 200 odds ratios per allele were calculated based on 200 simulated datasets generated by randomly permuting MHC individual genotypes with infection status (infected or uninfected by myxozoa) using Poptool 3.0.3 (Hood 2008). Simulations indicated that the observed odds ratios for alleles 7 and 8 were outside the theoretical expectation (p=0.020 and 0.025, respectively), with four and five simulated odds ratios out of 200 respectively lower and higher than the observed values for alleles 7 and 8. The Fisher's exact probability test at the amino acid level identified a valine in position 65 and a leucine in position 68 as amino acid positions associated with resistance of allele 7 and susceptibility of allele 8 to myxozoa infection, respectively (p<0.05). Phenylalanine or tyrosine and phenylalanine or histidine were otherwise found at positions 65 and 68, respectively. The multivariate logistic regression gave concordant results and revealed that the infection status of salmon depended on the region of origin, the sampling period and the presence of allele 8 (p=0.018, p<0.001 and p=0.019, respectively). Indeed, in concordance with the previous global logistic regression on infection rate, myxozoa infection was predominant in St Lawrence south shore relative to north shore rivers and was highest in June and decreased thereafter. Also, salmon with allele 8 had 3.4 times more chance to be infected with myxozoa than those not having this allele, when controlling for spatial and temporal variations in infection rate. Allele 7 was only retained by the model when the factor region was excluded from the analysis (p=0.031), indicating that the frequency of allele 7 was correlated with the regional variation, being slightly more frequent and less variable among north shore than south shore rivers (0.12±0.05 and 0.11±0.11, respectively). Still, according to the logistic regression and when controlling for spatial and temporal variations in infection rate, salmon carrying allele 7 had 2.9 times greater chance of being uninfected by myxozoa than salmon not carrying this allele. A comparison of allele frequencies for each microsatellite locus across time indicated no significant temporal change during summer (p>0.05). The frequency of MHC allele 7 increased slightly through time but this trend was not significant (p=0.30), while the frequency of MHC allele 8 decreased significantly through summer (p=0.003; figure 3). All other MHC alleles did not show any trend in frequency change during the summer (p>0.05).

View this table:
Table 5

Fisher's exact probability test on myxozoa infection status of juvenile Atlantic salmon as a function of MHC class IIβ alleles. (Values in bold represent alleles significantly associated with the probability of being infected by myxozoa in the wild according to the Fisher's exact probability test and the permutation test of significance.)

Figure 3

Average allele frequency of the identified MHC resistance allele 7 (filled circles) and MHC susceptibility allele 8 (open circles) to myxozoa across the summer in juvenile Atlantic salmon from six rivers in Québec. (Error bars represent standard deviation associated with the different populations. Asterisk indicates significance while n.s. indicates a non-significant change in allele frequency according to the Mantel–Haenszel Χ2-tests.)

4. Discussion

The goal of this study was to clarify the mechanism of pathogen resistance at the individual level and evaluate the spatio-temporal intensity of immune challenges in the wild by screening bacterial infections in juvenile Atlantic salmon and relating these to MHC and microsatellite variability. An overall kidney infection rate of 12.1 per cent was found in wild juvenile Atlantic salmon within the six rivers studied. Infection rate and pathogen diversity were greater at the beginning of the summer when an average proportion as high as one in every two juveniles was found to be infected in St Lawrence south shore rivers and one in every five in north shore rivers. The sex ratio of the infected individuals was slightly male biased, albeit not significantly, which suggested that both sexes had similar infection rates in the wild. A myxozoa parasite of the genus Myxobolus, most probably recently introduced in North America, was discovered and showed an average prevalence of 14.2 per cent in south shore rivers and 3.7 per cent in north shore rivers during summer. One MHC class IIβ allele was associated with resistance and another with susceptibility to the myxozoa, suggesting a role of MHC standing genetic variation in influencing infection levels towards an emerging parasite. As no evidence for heterozygote advantage was found, these results better support the frequency-dependent selection or the variable selection in time and space hypotheses. In addition, a significant decrease in frequency of the susceptibility allele but not other MHC or microsatellite alleles during summer suggested the occurrence of a mortality event due to myxozoa infection. Overall, these findings suggest that the myxozoa parasite can exert a selective pressure on its host, but that standing genetic variation at MHC in host populations could contribute to promoting adaptation over time.

Spatial variation in infectivity suggested a trend for higher infection rates and pathogen diversity in juvenile salmon from south shore compared with north shore rivers in the St Lawrence system. This result was concordant with the evidence of a higher pathogen selection pressure in southern warm than northern cold habitats found in a large-scale study conducted on wild Atlantic salmon populations (Dionne et al. 2007). Indeed, the bacterial community in each of the six rivers was analysed during the same year and bacterial diversity was higher in south shore than in north shore rivers (38.9±12.8 and 2.5±2.4×103 IOD mm−2, respectively; see Dionne et al. (2007) for bacterial community analyses). Temporal variation in infectivity also suggested higher infection rates in June when water temperature rise rapidly from approximately 6°C to over 18°C (Ministère des Ressources Naturelles et de la Faune du Québec). Pathogenic bacteria and other parasites may increase in abundance and virulence once a temperature threshold is reached, as is the case for myxozoa infections that mostly occur at more than 15°C (Uhland et al. 2000), and as previously observed in the laboratory for other fish pathogens (e.g. Larsen et al. 2004). The observed decrease in infection rate thereafter suggested that either juveniles lost their infection with time, following the pathogen natural life cycle or through immune resistance, or mortality occurred following infection. Although these hypotheses are not mutually exclusive, our data are more consistent with the latter hypothesis, as a significant decrease in frequency of a MHC susceptibility allele towards myxozoa infection, as well as a trend towards the increase in a MHC resistance allele, were observed during the summer while no significant change was detected for other MHC and microsatellite alleles. Overall, we found evidence for spatio-temporal variation in infection rate concordant with variation in pathogen selection pressure in the environment, and also with the immune competence associated with specific MHC alleles (see below).

Twelve species of potentially pathogenic bacteria were found in kidneys of wild juvenile salmon, all of which had a low prevalence in the populations. Some of them were known pathogens, namely Mycoplasma infecting farmed cod and other vertebrates (Messick et al. 2002; Nylund et al. 2006), the pathogenic bacteria Rickettsiella mainly infecting arthropods (Cordeaux et al. 2007), the opportunistic fish pathogen Pseudomonas causing septicaemia in stressful conditions (Thune et al. 1993), Shewanella responsible for ulcer disease in some molluscs and fish (Cai et al. 2006) and Piscirickettsia causing epizootic diseases in salmonids (Fryer & Hedrick 2003). Other bacteria, namely Sphingomonas, Methylobacterium, Bosea, Lactobacillus and Burkholderia, have been previously found in soil and freshwater environments and some such as Burkholderia can become opportunistic pathogens (Berriatua et al. 2001). Although prevalence was low for these bacteria, some of them such as Piscirickettsia are qualified as epizootic, affecting a small number of individuals but, under favourable environmental conditions, possessing the capacity to spread rapidly in the population (Fryer & Hedrick 2003). As such, these could constitute latent pathogens representing a potential health risk for salmon populations in the future. The most prevalent pathogen infecting juvenile salmon was a myxozoa of the Myxobolus clade, a group of fish parasites first discovered in farmed salmonids in Europe (Hoffman 1990). Myxobolus infections, as well as other documented myxozoa infections, have been detected in the USA since the 1950s and epidemiological and molecular evidence suggests that these parasites were initially introduced through fish transported from Europe (Hedrick et al. 1993; Andree et al. 1999; Kent et al. 2001; Miller & Vincent 2008). Myxobolus infections have been recently reported in eastern Canada in some fish species such as the log-perch, Percina caprodes, and the banded killifish, Fundulus diaphanus (Cone & Marcogliese 2004; Cone et al. 2006). However, to our knowledge, this is the first published report of a myxozoa infection in wild Atlantic salmon in eastern Canada. As such, this most probably represents an emerging salmon parasite in this system. Myxozoa parasites are responsible for multiple diseases in salmonids, including the well-documented whirling disease caused by Myxobolus cerebralis, which can severely compromise survival in some species (Kent et al. 2001). Even though MHC allele 7 covaried with regional variability, MHC alleles 7 and 8 were associated with resistance and susceptibility to myxozoa infection, respectively, after controlling for spatio-temporal variability in infection rate. Previous studies conducted in the laboratory also found an association between an MHC allele and incidence of diseases, namely furunculosis, infectious salmon anaemia and infectious haemopoietic necrosis virus in Atlantic salmon (Langefors et al. 2001; Grimholt et al. 2003; Miller et al. 2004). However, very few studies have attempted to test for an association between MHC alleles and parasite resistance in wild populations (but see Harf & Sommer (2005) and references therein on small mammals and Madsen & Ujvari (2006) on pythons). In addition to identifying an association between MHC alleles and susceptibility, we have also attempted to disentangle the influence of spatio-temporal variability in pathogen selection pressure from immune resistance processes in determining infection rates in the wild, which has always been a challenge in the past (Wegner et al. 2003). The significant shift in the frequency of the susceptibility allele 8 over the summer supports the hypothesis of ‘real-time’ myxozoa selection pressure on juvenile Atlantic salmon causing mortality over the summer. An increasing number of studies have demonstrated the effect of natural selection on MHC genes over an evolutionary time scale in the wild (Miller et al. 2001; Prugnolle et al. 2005; van Oosterhout et al. 2006; Dionne et al. 2007), but few have identified differential mortality associated with MHC allelic composition during a single generation (Piertney & Oliver 2006; but see de Eyto et al. 2007). Assuming that this susceptibility allele frequency change is strictly associated with mortality, we can speculate that the deleterious consequences of maladaptation in these populations resulted roughly in a 5 per cent reduction in population abundance from June to August (allele frequency changes from 0.088 to 0.040, figure 3). In another system, the introduction of M. cerebralis was associated with a drastic decline in rainbow trout (Oncorhynchus mykiss) recruitment over 18 years (Miller & Vincent 2008). Altogether, these results then suggest that Myxobolus parasites could have demographic consequences on infected populations, and underline the importance of immune adaptation for population persistence. Admittedly, more in-depth studies will be necessary to rigorously assess the deleterious consequences of immune maladaptation in salmon populations. Such studies will become increasingly important in the face of climate change, which may cause rapid variations in pathogen communities and selection pressures.

The heterozygote advantage hypothesis was not supported by our data, which might not be surprising in the context of a predominant pathogen, as suggested by Wegner et al. (2004). Admittedly, not all possible pathogens were identified in this study, which might have influenced the relationship observed between heterozygosity and infection rate. However, what was interesting was the detection of the reverse pattern, with significantly more MHC heterozygotes and higher MHC amino acid genetic distance in infected than in uninfected fish. A similar trend was also observed in Atlantic salmon with furunculosis in the laboratory (Langefors et al. 2001). One possible explanation could be the low frequency of the MHC susceptibility allele 8 (overall frequency less than 0.05), which would result in this allele being mostly restricted to heterozygous individuals, with the result that heterozygotes appear more infected than expected. In such a case, MHC amino acid genetic distance should follow this pattern only if the susceptibility allele is more differentiated in terms of amino acids compared with other alleles, which was indeed found at PBR but not at overall polymorphic sites (t-test comparing p-distances associated with allele 8 with p-distances associated with all other allele comparisons: p=0.029 and 0.39, respectively), in concordance with the results of the MANOVA. Additionally, higher MHC heterozygosity in infected fish could also be partially explained by spatial variability, as a tendency for higher MHC heterozygosity was observed in south shore rivers where infections were more frequent (HO south shore, 0.90±0.05; north shore, 0.81±0.05; t-test: p=0.085). This suggested that the heterozygote advantage hypothesis should be tested with care in the wild, as frequency and identity of MHC alleles conferring susceptibility as well as spatial variability in MHC heterozygosity could influence the relationship between MHC heterozygosity (or other MHC individual genetic distance) and infection rate.

Our results, showing an association between MHC alleles and infection level, are more compatible with the negative frequency-dependent selection and the variable selection in time and space hypotheses (Takahata & Nei 1990; Hedrick 2002), but possibly better support the latter. Although rare alleles were not found to be more advantageous than the others, this study suggested that at the individual level, when facing an emerging pathogen, possessing the right alleles to initiate an appropriate immune response may be important for contemporary adaptation. As the pathogen community is changing, different alleles could be selected for through time, maintaining MHC polymorphism at the population level over an evolutionary time frame, as suggested by the variable selection in time and space hypothesis (Hedrick 2002) and as evidenced in the temporal study of Westerdahl et al. (2004). Indeed, at the population level, balancing selection maintaining MHC diversity could represent the combined influence of multiple directional selection pressures associated with each pathogen over time. However, allele frequency and infection rates would need to be followed across generations to find a stronger support for these hypotheses and evaluate whether a resistance process will develop over an evolutionary time frame in these myxozoa-exposed populations.

Overall, this study represents one of the very few attempts to test for an association between MHC alleles and infection level in wild populations by differentiating the influence of spatio-temporal pathogen selection pressure from individual immune competence. In doing so, we documented the prevalence and type of pathogen representing a contemporary or a potential future selection pressure for Atlantic salmon in the wild. Moreover, this study identified a potential MHC-selective mortality event, and an emerging parasite was identified as the most probable selective agent, possibly representing the first steps of a coevolution process between a host and its pathogen in the wild. Finally, these results suggest the implication of MHC standing genetic variation in facing pathogen challenges in the wild and underline the importance of maintaining diversity at MHC genes in natural populations to increase chances of host contemporary adaptation in a changing environment.


Special thanks to F. Caron, D. Fournier, M. Dorais, J.-P. Lebel, R. Isabel, M. Talbot and M. Valentine from ‘Ministère des Ressources Naturelles et de la Faune’ (MRNF) and A. Boivin and J.-F. Bourque from the ‘Centre Inter-universitaire de Recherche sur le Saumon Atlantique’ (CIRSA) for their assistance in the field. R. Firth, J. Roy, R. Lamy and C. Bernier are also gratefully acknowledged for their help during sampling. We also thank K. Giguère, L. Papillon, V. Albert, M. Gauthier, L. Landry, R. Martel and T. Ming for their laboratory assistance. A special thanks to B. Wilson for his great help with the pathogen detection protocol and to C. Whipps for his generous help in identifying the myxozoa parasite. Funding for this project was provided by Natural Sciences and Engineering Research Council of Canada (NSERC, collaborative project) to L.B., J.J.D. and K.M.M. and from MRNF. M.D. was financially supported by an NSERC postgraduate scholarship. Diversitas is also acknowledged for its financial support to M.D. This is a contribution to the research programmes of Québec-Océan and CIRSA.


  • One contribution of 14 to a Theme Issue ‘Eco-evolutionary dynamics’.


View Abstract