MHC standing genetic variation and pathogen resistance in wild Atlantic salmon

Changes in capture availability due to infection can lead to detectable biases in population-level infectious disease parameters

Correctly identifying the strength of selection that parasites impose on hosts is key to predicting epidemiological and evolutionary outcomes of host-parasite interactions. However, behavioral changes due to infection can alter the capture probability of infected hosts and thereby make selection difficult to estimate by standard sampling techniques. Mark-recapture approaches, which allow researchers to determine if some groups in a population are less likely to be captured than others, can be used to identify infection-driven capture biases. If a metric of interest directly compares infected and uninfected populations, calculated detection probabilities for both groups may be useful in identifying bias. Here, we use an individual-based simulation to test whether changes in capture rate due to infection can alter estimates of three key metrics: 1) reduction in the reproductive success of infected parents relative to uninfected parents, 2) the relative risk of infection for susceptible genotypes compared to resistant genotypes, and 3) changes in allele frequencies between generations. We explore the direction and underlying causes of the biases that emerge from these simulations. Finally, we argue that short series of mark-recapture sampling bouts, potentially implemented in under a week, can yield key data on detection bias due to infection while not adding a significantly higher burden to disease ecology studies.

Adaptive Evolution of Nearctic Deepwater Fish Vision: Implications for Assessing Functional Variation for Conservation

Intraspecific functional variation is critical for adaptation to rapidly changing environments. For visual opsins, functional variation can be characterized in vitro and often reflects a species' ecological niche but is rarely considered in the context of intraspecific variation or the impact of recent environmental changes on species of cultural or commercial significance. Investigation of adaptation in postglacial lakes can provide key insight into how rapid environmental changes impact functional evolution. Here, we report evidence for molecular adaptation in vision in 2 lineages of Nearctic fishes that are deep lake specialists: ciscoes and deepwater sculpin. We found depth-related variation in the dim-light visual pigment rhodopsin that evolved convergently in these 2 lineages. In vitro characterization of spectral sensitivity of the convergent deepwater rhodopsin alleles revealed blue-shifts compared with other more widely distributed alleles. These blue-shifted rhodopsin alleles were only observed in deep clear postglacial lakes with underwater visual environments enriched in blue light. This provides evidence of remarkably rapid and convergent visual adaptation and intraspecific functional variation in rhodopsin. Intraspecific functional variation has important implications for conservation, and these fishes are of conservation concern and great cultural, commercial, and nutritional importance to Indigenous communities. We collaborated with the Saugeen Ojibway Nation to develop and test a metabarcoding approach that we show is efficient and accurate in recovering the ecological distribution of functionally relevant variation in rhodopsin. Our approach bridges experimental analyses of protein function and genetics-based tools used in large-scale surveys to better understand the ecological extent of adaptive functional variation.

Insights from a chum salmon (Oncorhynchus keta) genome assembly regarding whole-genome duplication and nucleotide variation influencing gene function

Chum salmon are ecologically important to Pacific Ocean ecosystems and commercially important to fisheries. To improve the genetic resources available for this species, we sequenced and assembled the genome of a male chum salmon using Oxford Nanopore read technology and the Flye genome assembly software (contig N50: ∼2 Mbp, complete BUSCOs: ∼98.1%). We also resequenced the genomes of 59 chum salmon from hatchery sources to better characterize the genome assembly and the diversity of nucleotide variants impacting phenotype variation. With genomic sequences from a doubled haploid individual, we were able to identify regions of the genome assembly that have been collapsed due to high sequence similarity between homeologous (duplicated) chromosomes. The homeologous chromosomes are relics of an ancient salmonid-specific genome duplication. These regions were enriched with genes whose functions are related to the immune system and responses to toxins. From analyzing nucleotide variant annotations of the resequenced genomes, we were also able to identify genes that have increased levels of variants thought to moderately impact gene function. Genes related to the immune system and the detection of chemical stimuli (olfaction) had increased levels of these variants based on a gene ontology enrichment analysis. The tandem organization of many of the enriched genes raises the question of why they have this organization.

Understanding the evolution of immune genes in jawed vertebrates

Driven by co-evolution with pathogens, host immunity continuously adapts to optimize defence against pathogens within a given environment. Recent advances in genetics, genomics and transcriptomics have enabled a more detailed investigation into how immunogenetic variation shapes the diversity of immune responses seen across domestic and wild animal species. However, a deeper understanding of the diverse molecular mechanisms that shape immunity within and among species is still needed to gain insight into-and generate evolutionary hypotheses on-the ultimate drivers of immunological differences. Here, we discuss current advances in our understanding of molecular evolution underpinning jawed vertebrate immunity. First, we introduce the immunome concept, a framework for characterizing genes involved in immune defence from a comparative perspective, then we outline how immune genes of interest can be identified. Second, we focus on how different selection modes are observed acting across groups of immune genes and propose hypotheses to explain these differences. We then provide an overview of the approaches used so far to study the evolutionary heterogeneity of immune genes on macro and microevolutionary scales. Finally, we discuss some of the current evidence as to how specific pathogens affect the evolution of different groups of immune genes. This review results from the collective discussion on the current key challenges in evolutionary immunology conducted at the ESEB 2021 Online Satellite Symposium: Molecular evolution of the vertebrate immune system, from the lab to natural populations.

Microsatellites as Molecular Markers with Applications in Exploitation and Conservation of Aquatic Animal Populations

A large number of species and taxa has been studied for genetic polymorphism. Microsatellites have been known as hypervariable neutral molecular markers with the highest resolution power in comparison with any other markers. However, the discovery of a new type of molecular marker—single nucleotide polymorphism (SNP) has put the existing applications of microsatellites to the test. To ensure good resolution power in studies of populations and individuals, a number of microsatellite loci from 14 to 20 was often used, which corresponds to about 200 independent alleles. Recently, these numbers have tended to be increased by the application of genomic sequencing of expressed sequence tags (ESTs), and the choice of the most informative loci for genotyping depends on the aims of research. Examples of successful applications of microsatellite molecular markers in aquaculture, fisheries, and conservation genetics in comparison with SNPs have been summarized in this review. Microsatellites can be considered superior markers in such topics as kinship and parentage analysis in cultured and natural populations, the assessment of gynogenesis, androgenesis and ploidization. Microsatellites can be coupled with SNPs for mapping QTL. Microsatellites will continue to be used in research on genetic diversity in cultured stocks, and also in natural populations as an economically advantageous genotyping technique.

Long term patterns of association between MHC and helminth burdens in the bank vole support Red Queen dynamics

Major histocompatibility complex (MHC) genes encode proteins crucial for adaptive immunity of vertebrates. Negative frequency-dependent selection (NFDS), resulting from adaptation of parasites to common MHC types, has been hypothesized to maintain high, functionally relevant polymorphism of MHC, but demonstration of this relationship has remained elusive. In particular, differentiation of NFDS from fluctuating selection, resulting from changes in parasite communities in time and space (FS), has proved difficult in short-term studies. Here, we used temporal data, accumulated through long-term monitoring of helminths infecting bank voles (Myodes glareolus), to test specific predictions of NFDS on MHC class II. Data were collected in three, moderately genetically differentiated subpopulations in Poland, which were characterized by some stable spatiotemporal helminth communities but also events indicating introduction of new species and loss of others. We found a complex association between individual MHC diversity and species richness, where intermediate numbers of DRB supertypes correlated with lowest species richness, but the opposite was true for DQB supertypes - arguing against universal selection for immunogenetic optimality. We also showed that particular MHC supertypes explain a portion of the variance in prevalence and abundance of helminths, but this effect was subpopulation-specific, which is consistent with both NFDS and FS. Finally, in line with NFDS, we found that certain helminths that have recently colonized or spread in a given subpopulation, more frequently or intensely infected voles with MHC supertypes that have been common in the recent past. Overall, our results highlight complex spatial and temporal patterns of MHC-parasite associations, the latter being consistent with Red Queen coevolutionary dynamics.

Identification of infectious agents in early marine Chinook and Coho salmon associated with cohort survival

Recent decades have seen an increased appreciation for the role infectious diseases can play in mass mortality events across a diversity of marine taxa. At the same time many Pacific salmon populations have declined in abundance as a result of reduced marine survival. However, few studies have explicitly considered the potential role pathogens could play in these declines. Using a multi-year dataset spanning 59 pathogen taxa in Chinook and Coho salmon sampled along the British Columbia coast, we carried out an exploratory analysis to quantify evidence for associations between pathogen prevalence and cohort survival and between pathogen load and body condition. While a variety of pathogens had moderate to strong negative correlations with body condition or survival for one host species in one season, we found that Tenacibaculum maritimum and Piscine orthoreovirus had consistently negative associations with body condition in both host species and seasons and were negatively associated with survival for Chinook salmon collected in the fall and winter. Our analyses, which offer the most comprehensive examination of associations between pathogen prevalence and Pacific salmon survival to date, suggest that pathogens in Pacific salmon warrant further attention, especially those whose distribution and abundance may be influenced by anthropogenic stressors.

The pink salmon genome: Uncovering the genomic consequences of a two-year life cycle

Pink salmon (Oncorhynchus gorbuscha) adults are the smallest of the five Pacific salmon native to the western Pacific Ocean. Pink salmon are also the most abundant of these species and account for a large proportion of the commercial value of the salmon fishery worldwide. A two-year life history of pink salmon generates temporally isolated populations that spawn either in even-years or odd-years. To uncover the influence of this genetic isolation, reference genome assemblies were generated for each year-class and whole genome re-sequencing data was collected from salmon of both year-classes. The salmon were sampled from six Canadian rivers and one Japanese river. At multiple centromeres we identified peaks of Fst between year-classes that were millions of base-pairs long. The largest Fst peak was also associated with a million base-pair chromosomal polymorphism found in the odd-year genome near a centromere. These Fst peaks may be the result of a centromere drive or a combination of reduced recombination and genetic drift, and they could influence speciation. Other regions of the genome influenced by odd-year and even-year temporal isolation and tentatively under selection were mostly associated with genes related to immune function, organ development/maintenance, and behaviour.

Environmentally associated chromosomal structural variation influences fine-scale population structure of Atlantic Salmon (Salmo salar)

Chromosomal rearrangements (e.g., inversions, fusions, and translocations) have long been associated with environmental variation in wild populations. New genomic tools provide the opportunity to examine the role of these structural variants in shaping adaptive differences within and among wild populations of non-model organisms. In Atlantic Salmon (Salmo salar), variations in chromosomal rearrangements exist across the species natural range, yet the role and importance of these structural variants in maintaining adaptive differences among wild populations remains poorly understood. We genotyped Atlantic Salmon (n = 1429) from 26 populations within a highly genetically structured region of southern Newfoundland, Canada with a 220K SNP array. Multivariate analysis, across two independent years, consistently identified variation in a structural variant (translocation between chromosomes Ssa01 and Ssa23), previously associated with evidence of trans-Atlantic secondary contact, as the dominant factor influencing population structure in the region. Redundancy analysis suggested that variation in the Ssa01/Ssa23 chromosomal translocation is strongly correlated with temperature. Our analyses suggest environmentally mediated selection acting on standing genetic variation in genomic architecture introduced through secondary contact may underpin fine-scale local adaptation in Placentia Bay, Newfoundland, Canada, a large and deep embayment, highlighting the importance of chromosomal structural variation as a driver of contemporary adaptive divergence.

Does genetic diversity protect host populations from parasites? A meta-analysis across natural and agricultural systems

If parasites transmit more readily between closely related hosts, then parasite burdens should decrease with increased genetic diversity of host populations. This important hypothesis is often accepted at face value—notorious epidemics of crop monocultures testify to the vulnerability of host populations that have been purged of diversity. Yet the relationship between genetic diversity and parasitism likely varies across contexts, differing between crop and noncrop hosts and between experimental and natural host populations. Here, we used a meta‐analytic approach to ask if host diversity confers protection against parasites over the range of contexts in which it has been tested. We synthesized the results of 102 studies, comprising 2004 effect sizes representing a diversity of approaches and host‐parasite systems. Our results validate a protective effect of genetic diversity, while revealing significant variation in its strength across biological and empirical contexts. In experimental host populations, genetic diversity reduces parasitism by ∼20% for noncrop hosts and by ∼50% for crop hosts. In contrast, observational studies of natural host populations show no consistent relationship between genetic diversity and parasitism, with both strong negative and positive correlations reported. This result supports the idea that, if parasites preferentially attack close relatives, the correlation of genetic diversity with parasitism could be positive or negative depending upon the potential for host populations to evolve in response to parasite selection. Taken together, these results reinforce genetic diversity as a priority for both conservation and agriculture and emphasize the challenges inherent to drawing comparisons between controlled experimental populations and dynamic natural populations.

STAT3/SOCS3 axis contributes to the outcome of salmonid whirling disease

There are differences in disease susceptibility to whirling disease (WD) among strains of rainbow trout. The North American strain Trout Lodge (TL) is highly susceptible, whereas the German Hofer (HO) strain is more resistant. The suppressor of cytokine signaling (SOCS) proteins are key in inhibiting cytokine signaling. Their role in modulating the immune response against whirling disease is not completely clear. This study aimed at investigating the transcriptional response of SOCS1 and SOCS3 genes to Myxobolus cerebralis along with that of several upstream regulators and immune response genes. M. cerebralis induced the expression of SOCS1, the IL-6-dependent SOCS3, the anti-inflammatory cytokine IL-10 and the Treg associated transcription factor FOXP3 in TL fish at multiple time points, which likely caused a restricted STAT1 and STAT3 activity affecting the Th17/Treg17 balance. The expression of SOCS1 and the IL-6-dependent SOCS3 was induced constraining the activation of STAT1 and STAT3 in TL fish, thereby causing Th17/Treg17 imbalance and leaving the fish unable to establish a protective immune response against M. cerebralis or control inflammatory reactions increasing susceptibility to WD. Conversely, in HO fish, the expression of SOCS1 and SOCS3 was restrained, whereas the expression of STAT1 and IL-23-mediated STAT3 was induced potentially enabling more controlled immune responses, accelerating parasite clearance and elevating resistance. The induced expression of STAT1 and IL-23-mediated STAT3 likely maintained a successful Th17/Treg17 balance and enabled fish to promote effective immune responses favouring resistance against WD. The results provide insights into the role of SOCS1 and SOCS3 in regulating the activation and magnitude of host immunity in rainbow trout, which may help us understand the mechanisms that underlie the variation in resistance to WD.

Genome-Wide Association Study Confirms Previous Findings of Major Loci Affecting Resistance to Piscine myocarditis virus in Atlantic Salmon (Salmo salar L.)

Cardiomyopathy syndrome is a viral disease of Atlantic salmon, mostly affecting fish during the late stages of production, resulting in significant losses to the industry. It has been shown that resistance to this disease has a strong genetic component, with quantitative trait loci (QTL) on chromosomes 27 (Ssa27) and Ssa12 to explain most of the additive genetic variance. Here, by analysing animals from a different year-class and a different population, we further aimed to confirm and narrow down the locations of these QTL. The data support the existence of the two QTL and suggest that the causative mutation on Ssa27 is most likely within the 10–10.5 Mbp segment of this chromosome. This region contains a cluster of major histocompatibility complex class I (MHC I) genes with the most strongly associated marker mapped to one of these loci. On Ssa12, the data confirmed the previous finding that the location of the causative mutation is within the 61.3 to 61.7 Mbp region. This segment contains several immune-related genes, but of particular interest are genes related to MHC II. Together, these findings highlight the likely key role of MHC genes in Atlantic salmon following infection with Piscine myocarditis virus (PMCV) and their potential impact on influencing the trajectory of this disease.

Specific pathogens and microbial abundance within liver and kidney tissues of wild marine fish from the Eastern Mediterranean Sea

This study is an initial description and discussion of the kidney and liver microbial communities of five common fish species sampled from four sites along the Eastern Mediterranean Sea shoreline. The goals of the present study were to establish a baseline dataset of microbial communities associated with the tissues of wild marine fish, in order to examine species-specific microbial characteristics and to screen for candidate pathogens. This issue is especially relevant due to the development of mariculture farms and the possible transmission of pathogens from wild to farmed fish and vice versa. Although fish were apparently healthy, 16S rRNA NGS screening identified three potential fish bacterial pathogens: Photobacterium damselae, Vibrio harveyi and Streptococcus iniae. Based on the distribution patterns and relative abundance, 16 samples were classified as potential pathogenic bacteria-infected samples (PPBIS). Hence, PPBIS prevalence was significantly higher in kidneys than in liver samples and variation was found between the fish species. Significant differences were observed between fish species, organs and sites, indicating the importance of the environmental conditions on the fish microbiome. We applied a consistent sampling and analytical method for monitoring in long-term surveys which may be incorporated within other marine fish pathogens surveys around the world.

Divergent Allele Advantage Provides a Quantitative Model for Maintaining Alleles with a Wide Range of Intrinsic Merits

The Major Histocompatibility Complex (MHC) is the most genetically diverse region of the genome in most vertebrates. Some form of balancing selection is necessary to account for the extreme diversity, but the precise mechanism of balancing selection is unknown. Due to the way MHC molecules determine immune recognition, overdominance (also referred to as heterozygote advantage) has been suggested as the main driving force behind this unrivalled diversity. However, both theoretical results and simulation models have shown that overdominance in its classical form cannot maintain large numbers of alleles unless all alleles confer unrealistically similar levels of fitness. There is increasing evidence that heterozygotes containing genetically divergent alleles allow for broader antigen presentation to immune cells, providing a selective mechanism for MHC polymorphism. By framing competing models of overdominance within a general framework, we show that a model based on Divergent Allele Advantage (DAA) provides a superior mechanism for maintaining alleles with a wide range of intrinsic merits, as intrinsically less-fit MHC alleles that are more divergent can survive under DAA. Specifically, our results demonstrate that a quantitative mechanism built from the DAA hypothesis is able to maintain polymorphism in the MHC. Applying such a model to both livestock breeding and conservation could provide a better way of identifying superior heterozygotes, and quantifying the advantages of genetic diversity at the MHC.

Environmental extremes drive population structure at the northern range limit of Atlantic salmon in North America

Conservation of exploited species requires an understanding of both genetic diversity and the dominant structuring forces, particularly near range limits, where climatic variation can drive rapid expansions or contractions of geographic range. Here, we examine population structure and landscape associations in Atlantic salmon (Salmo salar) across a heterogeneous landscape near the northern range limit in Labrador, Canada. Analysis of two amplicon-based data sets containing 101 microsatellites and 376 single nucleotide polymorphisms (SNPs) from 35 locations revealed clear differentiation between populations spawning in rivers flowing into a large marine embayment (Lake Melville) compared to coastal populations. The mechanisms influencing the differentiation of embayment populations were investigated using both multivariate and machine-learning landscape genetic approaches. We identified temperature as the strongest correlate with genetic structure, particularly warm temperature extremes and wider annual temperature ranges. The genomic basis of this divergence was further explored using a subset of locations (n = 17) and a 220K SNP array. SNPs associated with spatial structuring and temperature mapped to a diverse set of genes and molecular pathways, including regulation of gene expression, immune response, and cell development and differentiation. The results spanning molecular marker types and both novel and established methods clearly show climate-associated, fine-scale population structure across an environmental gradient in Atlantic salmon near its range limit in North America, highlighting valuable approaches for predicting population responses to climate change and managing species sustainability.

Neutral variation does not predict immunogenetic variation in the European grayling (Thymallus thymallus)-implications for management

Preservation of genetic diversity is critical to successful conservation, and there is increasing demand for the inclusion of ecologically meaningful genetic information in management decisions. Supportive breeding programmes are increasingly implemented to combat declines in many species, yet their effect on adaptive genetic variation is understudied. This is despite the fact that supportive breeding may interfere with natural evolutionary processes. Here, we assessed the performance of neutral and adaptive markers (major histocompatibility complex; MHC) to inform management of European grayling (Thymallus thymallus), which routinely involves supplementation of natural populations with hatchery-reared fish (stocking). This study is the first to characterize MH II DAA and DAB loci in grayling and to investigate immune genetic variation in relation to management practice in this species. High-throughput Illumina sequencing of "introduced," "stocked native" and "non-stocked native" populations revealed significantly higher levels of allelic richness and heterozygosity for MH markers than microsatellites exclusively in non-stocked native populations. Likewise, significantly lower differentiation at the MH II than for microsatellites was apparent when considering non-stocked native populations, but not stocked populations. We developed a simulation model to test the effects of relaxation of selection during the early life stage within captivity. Dependent on the census population size and stocking intensity, there may be long-term effects of stocking on MH II, but not neutral genetic diversity. This is consistent with our empirical results. This study highlights the necessity for considering adaptive genetic variation in conservation decisions and raises concerns about the efficiency of stocking as a management practice.

Range-wide parallel climate-associated genomic clines in Atlantic salmon

Clinal variation across replicated environmental gradients can reveal evidence of local adaptation, providing insight into the demographic and evolutionary processes that shape intraspecific diversity. Using 1773 genome-wide single nucleotide polymorphisms we evaluated latitudinal variation in allele frequency for 134 populations of North American and European Atlantic salmon (Salmo salar). We detected 84 (4.74%) and 195 (11%) loci showing clinal patterns in North America and Europe, respectively, with 12 clinal loci in common between continents. Clinal single nucleotide polymorphisms were evenly distributed across the salmon genome and logistic regression revealed significant associations with latitude and seasonal temperatures, particularly average spring temperature in both continents. Loci displaying parallel clines were associated with several metabolic and immune functions, suggesting a potential basis for climate-associated adaptive differentiation. These climate-based clines collectively suggest evidence of large-scale environmental associated differences on either side of the North Atlantic. Our results support patterns of parallel evolution on both sides of the North Atlantic, with evidence of both similar and divergent underlying genetic architecture. The identification of climate-associated genomic clines illuminates the role of selection and demographic processes on intraspecific diversity in this species and provides a context in which to evaluate the impacts of climate change.

Do genetic drift and accumulation of deleterious mutations preclude adaptation? Empirical investigation using RADseq in a northern lacustrine fish

Understanding genomic signatures of divergent selection underlying long-term adaptation in populations located in heterogeneous environments is a key goal in evolutionary biology. In this study, we investigated neutral, adaptive and deleterious genetic variation using 7,192 SNPs in 31 Lake Trout (Salvelinus namaycush) populations (n = 673) from Québec, Canada. Average genetic diversity was low, weakly shared among lakes, and positively correlated with lake size, indicating a major role for genetic drift subsequent to lake isolation. Putatively deleterious mutations were on average at lower frequencies than the other SNPs, and their abundance relative to the entire polymorphism in each population was positively correlated with inbreeding, suggesting that the effectiveness of purifying selection was negatively correlated with inbreeding, as predicted from theory. Despite evidence for pronounced genetic drift and inbreeding, several outlier loci were associated with temperature and found in or close to genes with biologically relevant functions notably related to heat stress and immune responses. Outcomes of gene-temperature associations were influenced by the inclusion of the most inbred populations, in which allele frequencies deviated the most from model predictions. This result illustrates challenge in identifying gene-environment associations in cases of high genetic drift and restricted gene flow and suggests limited adaptation in populations experiencing higher inbreeding. We discuss the relevance of these findings for the conservation and management, notably regarding stocking and genetic rescue, of Lake Trout populations and other species inhabiting highly fragmented habitats.

Drift, selection, or migration? Processes affecting genetic differentiation and variation along a latitudinal gradient in an amphibian

BACKGROUND

Past events like fluctuations in population size and post-glacial colonization processes may influence the relative importance of genetic drift, migration and selection when determining the present day patterns of genetic variation. We disentangle how drift, selection and migration shape neutral and adaptive genetic variation in 12 moor frog populations along a 1700 km latitudinal gradient. We studied genetic differentiation and variation at a MHC exon II locus and a set of 18 microsatellites.

RESULTS

Using outlier analyses, we identified the MHC II exon 2 (corresponding to the β-2 domain) locus and one microsatellite locus (RCO8640) to be subject to diversifying selection, while five microsatellite loci showed signals of stabilizing selection among populations. STRUCTURE and DAPC analyses on the neutral microsatellites assigned populations to a northern and a southern cluster, reflecting two different post-glacial colonization routes found in previous studies. Genetic variation overall was lower in the northern cluster. The signature of selection on MHC exon II was weaker in the northern cluster, possibly as a consequence of smaller and more fragmented populations.

CONCLUSION

Our results show that historical demographic processes combined with selection and drift have led to a complex pattern of differentiation along the gradient where some loci are more divergent among populations than predicted from drift expectations due to diversifying selection, while other loci are more uniform among populations due to stabilizing selection. Importantly, both overall and MHC genetic variation are lower at northern latitudes. Due to lower evolutionary potential, the low genetic variation in northern populations may increase the risk of extinction when confronted with emerging pathogens and climate change.

Population Differences at MHC Do Not Explain Enhanced Resistance of Song Sparrows to Local Parasites

Infectious disease represents an emerging threat to natural populations, particularly when hosts are more susceptible to novel parasites (allopatric) than to parasites from the local area (sympatric). This pattern could arise through evolutionary processes (host populations become adapted to their local parasites and genetically differentiated from other populations at immune-related loci) and/or through ecological interactions (host individuals develop resistance to local parasites through previous exposure). The relative importance of these candidate mechanisms remains unclear. In jawed vertebrates, genes of the major histocompatibility complex (MHC) play a fundamental role in immunity and are compelling candidates for spatially varying selection. We recently showed that song sparrows (Melospiza melodia) are more susceptible to allopatric than to sympatric strains of malaria (Plasmodium). In the current study, to determine whether population differences at MHC explain this pattern, we characterized the peptide-binding regions of MHC (classes I and II) of birds that did or did not become infected in the previous experiment. We recovered up to 4 alleles per individual at class I, implying at least 2 loci, and up to 26 alleles per individual at class II, implying at least 13 loci. Individuals with more class I alleles were less likely to become infected by Plasmodium, consistent with parasite-mediated balancing selection. However, we found no evidence for population genetic differentiation at either class of MHC, based on 36 individuals sequenced. Resistance to sympatric parasites previously described for this system likely stems from individuals' prior immune experience, not from population differentiation and locally protective alleles at MHC.

Small-scale intraspecific patterns of adaptive immunogenetic polymorphisms and neutral variation in Lake Superior lake trout

Many fishes express high levels of intraspecific variability, often linked to resource partitioning. Several studies show that a species' evolutionary trajectory of adaptive divergence can undergo reversals caused by changes in its environment. Such a reversal in neutral genetic and morphological variation among lake trout Salvelinus namaycush ecomorphs appears to be underway in Lake Superior. However, a water depth gradient in neutral genetic divergence was found to be associated with intraspecific diversity in the lake. To investigate patterns of adaptive immunogenetic variation among lake trout ecomorphs, we used Illumina high-throughput sequencing. The population's genetic structure of the major histocompatibility complex (MHC Class IIβ exon 2) and 18 microsatellite loci were compared to disentangle neutral and selective processes at a small geographic scale. Both MHC and microsatellite variation were partitioned more by water depth stratum than by ecomorph. Several metrics showed strong clustering by water depth in MHC alleles, but not microsatellites. We report a 75% increase in the number of MHC alleles shared between the predominant shallow and deep water ecomorphs since a previous lake trout MHC study at the same locale (c. 1990s data). This result is consistent with the reverse speciation hypothesis, although adaptive MHC polymorphisms persist along an ecological gradient. Finally, results suggested that the lake trout have multiple copies of the MHC II locus consistent with a historic genomic duplication event. Our findings indicated that conservation approaches for this species could focus on managing various ecological habitats by depth, in addition to regulating the fisheries specific to ecomorphs.

Cheetahs have a stronger constitutive innate immunity than leopards

As a textbook case for the importance of genetics in conservation, absence of genetic variability at the major histocompatibility complex (MHC) is thought to endanger species viability, since it is considered crucial for pathogen resistance. An alternative view of the immune system inspired by life history theory posits that a strong response should evolve in other components of the immune system if there is little variation in the MHC. In contrast to the leopard (Panthera pardus), the cheetah (Acinonyx jubatus) has a relatively low genetic variability at the MHC, yet free-ranging cheetahs are healthy. By comparing the functional competence of the humoral immune system of both species in sympatric populations in Namibia, we demonstrate that cheetahs have a higher constitutive innate but lower induced innate and adaptive immunity than leopards. We conclude (1) immunocompetence of cheetahs is higher than previously thought; (2) studying both innate and adaptive components of immune systems will enrich conservation science.

Spatio-temporal variation in parasite communities maintains diversity at the major histocompatibility complex class IIβ in the endangered Rio Grande silvery minnow

Climate change will strongly impact aquatic ecosystems particularly in arid and semi-arid regions. Fish-parasite interactions will also be affected by predicted altered flow and temperature regimes, and other environmental stressors. Hence, identifying environmental and genetic factors associated with maintaining diversity at immune genes is critical for understanding species' adaptive capacity. Here, we combine genetic (MHC class IIβ and microsatellites), parasitological and ecological data to explore the relationship between these factors in the remnant wild Rio Grande silvery minnow (Hybognathus amarus) population, an endangered species found in the southwestern United States. Infections with multiple parasites on the gills were observed and there was spatio-temporal variation in parasite communities and patterns of infection among individuals. Despite its highly endangered status and chronically low genetic effective size, Rio Grande silvery minnow had high allelic diversity at MHC class IIβ with more alleles recognized at the presumptive DAB1 locus compared to the DAB3 locus. We identified significant associations between specific parasites and MHC alleles against a backdrop of generalist parasite prevalence. We also found that individuals with higher individual neutral heterozygosity and higher amino acid divergence between MHC alleles had lower parasite abundance and diversity. Taken together, these results suggest a role for fluctuating selection imposed by spatio-temporal variation in pathogen communities and divergent allele advantage in maintenance of high MHC polymorphism. Understanding the complex interaction of habitat, pathogens and immunity in protected species will require integrated experimental, genetic and field studies.

Climatic and evolutionary drivers of phase shifts in the plague epidemics of colonial India

Immune heterogeneity in wild host populations indicates that disease-mediated selection is common in nature. However, the underlying dynamic feedbacks involving the ecology of disease transmission, evolutionary processes, and their interaction with environmental drivers have proven challenging to characterize. Plague presents an optimal system for interrogating such couplings: Yersinia pestis transmission exerts intense selective pressure driving the local persistence of disease resistance among its wildlife hosts in endemic areas. Investigations undertaken in colonial India after the introduction of plague in 1896 suggest that, only a decade after plague arrived, a heritable, plague-resistant phenotype had become prevalent among commensal rats of cities undergoing severe plague epidemics. To understand the possible evolutionary basis of these observations, we developed a mathematical model coupling environmentally forced plague dynamics with evolutionary selection of rats, capitalizing on extensive archival data from Indian Plague Commission investigations. Incorporating increased plague resistance among rats as a consequence of intense natural selection permits the model to reproduce observed changes in seasonal epidemic patterns in several cities and capture experimentally observed associations between climate and flea population dynamics in India. Our model results substantiate Victorian era claims of host evolution based on experimental observations of plague resistance and reveal the buffering effect of such evolution against environmental drivers of transmission. Our analysis shows that historical datasets can yield powerful insights into the transmission dynamics of reemerging disease agents with which we have limited contemporary experience to guide quantitative modeling and inference.

Spatial patterns of immunogenetic and neutral variation underscore the conservation value of small, isolated American badger populations

Small and isolated populations often exhibit low genetic diversity due to drift and inbreeding, but may simultaneously harbour adaptive variation. We investigate spatial distributions of immunogenetic variation in American badger subspecies (Taxidea taxus), as a proxy for evaluating their evolutionary potential across the northern extent of the species' range. We compared genetic structure of 20 microsatellites and the major histocompatibility complex (MHC DRB exon 2) to evaluate whether small, isolated populations show low adaptive polymorphism relative to large and well-connected populations. Our results suggest that gene flow plays a prominent role in shaping MHC polymorphism across large spatial scales, while the interplay between gene flow and selection was stronger towards the northern peripheries. The similarity of MHC alleles within subspecies relative to their neutral genetic differentiation suggests that adaptive divergence among subspecies can be maintained despite ongoing gene flow along subspecies boundaries. Neutral genetic diversity was low in small relative to large populations, but MHC diversity within individuals was high in small populations. Despite reduced neutral genetic variation, small and isolated populations harbour functional variation that likely contribute to the species evolutionary potential at the northern range. Our findings suggest that conservation approaches should focus on managing adaptive variation across the species range rather than protecting subspecies per se.

Atlantic salmon populations reveal adaptive divergence of immune related genes - a duplicated genome under selection

Background

Populations of Atlantic salmon display highly significant genetic differences with unresolved molecular basis. These differences may result from separate postglacial colonization patterns, diversifying natural selection and adaptation, or a combination. Adaptation could be influenced or even facilitated by the recent whole genome duplication in the salmonid lineage which resulted in a partly tetraploid species with duplicated genes and regions.

Results

In order to elucidate the genes and genomic regions underlying the genetic differences, we conducted a genome wide association study using whole genome resequencing data from eight populations from Northern and Southern Norway. From a total of ~4.5 million sequencing-derived SNPs, more than 10 % showed significant differentiation between populations from these two regions and ten selective sweeps on chromosomes 5, 10, 11, 13–15, 21, 24 and 25 were identified. These comprised 59 genes, of which 15 had one or more differentiated missense mutation. Our analysis showed that most sweeps have paralogous regions in the partially tetraploid genome, each lacking the high number of significant SNPs found in the sweeps. The most significant sweep was found on Chr 25 and carried several missense mutations in the antiviral mx genes, suggesting that these populations have experienced differing viral pressures. Interestingly the second most significant sweep, found on Chr 5, contains two genes involved in the NF-KB pathway (nkap and nkrf), which is also a known pathogen target that controls a large number of processes in animals.

Conclusion

Our results show that natural selection acting on immune related genes has contributed to genetic divergence between salmon populations in Norway. The differences between populations may have been facilitated by the plasticity of the salmon genome. The observed signatures of selection in duplicated genomic regions suggest that the recently duplicated genome has provided raw material for evolutionary adaptation.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-2867-z) contains supplementary material, which is available to authorized users.

Cryptic female choice enhances fertilization success and embryo survival in chinook salmon

In this study, we investigated two potentially important intersexual postcopulatory gametic interactions in a population of chinook salmon (Oncorhynchus tshawytscha): (i) the effect of female ovarian fluid (OF) on the behaviour of spermatozoa during fertilization and (ii) the effects of multilocus heterozygosity (MLH) (as an index of male quality) and female-male genetic relatedness on sperm behaviour and male fertilization success when there is sperm competition in the presence of that OF. To do this, we conducted a series of in vitro competitive fertilization experiments and found that, when ejaculates from two males are competing for access to a single female's unfertilized eggs, fertilization success was significantly biased towards the male whose sperm swam fastest in the female's OF. Embryo survival--a measure of fitness--was also positively correlated with both sperm swimming speed in OF and male MLH, providing novel evidence that cryptic female choice is adaptive for the female, enhancing the early survival of her offspring and potentially influencing her fitness.

Potential drivers of virulence evolution in aquaculture

Infectious diseases are economically detrimental to aquaculture, and with continued expansion and intensification of aquaculture, the importance of managing infectious diseases will likely increase in the future. Here, we use evolution of virulence theory, along with examples, to identify aquaculture practices that might lead to the evolution of increased pathogen virulence. We identify eight practices common in aquaculture that theory predicts may favor evolution toward higher pathogen virulence. Four are related to intensive aquaculture operations, and four others are related specifically to infectious disease control. Our intention is to make aquaculture managers aware of these risks, such that with increased vigilance, they might be able to detect and prevent the emergence and spread of increasingly troublesome pathogen strains in the future.

Does the parasite-mediated selection drive the MHC class IIB diversity in wild populations of European chub (Squalius cephalus)?

The genes of major histocompatibility complex (MHC) provide an excellent opportunity to study host-parasite relationships because they are expected to evolve in response to parasites and variation in parasite communities. In this study, we investigated the potential role of parasite-mediated selection acting on MHC class IIB (DAB) genes in European chub (Squalius cephalus) natural populations. We found significant differences between populations in metazoan parasites, neutral and adaptive genetic diversities. The analyses based on pairwise data revealed that populations with dissimilar MHC allelic profiles were geographically distant populations with significantly different diversity in microsatellites and a dissimilar composition of parasite communities. The results from the generalized estimating equations method (GEE) on the level of individuals revealed that metazoan parasite load in European chub was influenced by the diversity of DAB alleles as well as by the diversity of neutral genetic markers and host traits reflecting condition and immunocompetence. The multivariate co-inertia analysis showed specific associations between DAB alleles and parasite species. DAB1-like alleles were more involved in associations with ectoparasites, while DAB3-like alleles were positively associated with endoparasites which could suggest potential differences between DAB genes caused by different selection pressure. Our study revealed that parasite-mediated selection is not the only variable affecting MHC diversity in European chub; however, we strongly support the role of neutral processes as the main driver of DAB diversity across populations. In addition, our study contributes to the understanding of the evolution of MHC genes in wild living fish.

Influence of immune-relevant genotype on the reproductive success of a salmonid alternative mating strategy

Major histocompatibility complex (MHC) and immune-relevant gene markers were used to evaluate differences in reproductive success (RS) among naturally spawning coho salmon Oncorhynchus kisutch mate pairs involving an alternative male reproductive phenotype, known as jacks. These mate pairs included both hatchery-reared and wild origin fish such that three classes were evaluated in two consecutive years (2005 and 2006) using a previously constructed multigenerational genetic pedigree: wild × wild (W × W), hatchery × hatchery (H × H) and wild × hatchery (W × H). Oncorhynchus kisutch jack mate pairs mated randomly based on immune-relevant genotype in both years; a result consistent with the opportunistic mating strategy of jacks. An association between greater number of alleles shared at three immune-relevant gene markers and increased RS was found for: W × H mate pairs in 2005 (BHMS429), W × H pairs in 2006 (SsalR016TKU) and W × W pairs in 2006 (OMM3085). No correlation between immune gene diversity and RS was found for H × H pairs in either year. The results suggest that the influence of immune-relevant genotype on mating success may be different for jacks when compared with previous studies of large adult male O. kisutch.

Diversity and relative abundance of the bacterial pathogen, Flavobacterium spp., infecting reproductive ecotypes of kokanee salmon

Background

Understanding the distribution and abundance of pathogens can provide insight into the evolution and ecology of their host species. Previous research in kokanee, the freshwater form of sockeye salmon (Oncorhynchus nerka), found evidence that populations spawning in streams may experience a greater pathogen load compared with populations that spawn on beaches. In this study we tested for differences in the abundance and diversity of the gram-negative bacteria, Flavobacterium spp., infecting tissues of kokanee in both of these spawning habitats (streams and beaches). Molecular assays were carried out using primers designed to amplify a ~200 nucleotide region of the gene encoding the ATP synthase alpha subunit (AtpA) within the genus Flavobacterium. Using a combination of DNA sequencing and quantitative PCR (qPCR) we compared the diversity and relative abundance of Flavobacterium AtpA amplicons present in DNA extracted from tissue samples of kokanee collected from each spawning habitat.

Results

We identified 10 Flavobacterium AtpA haplotypes among the tissues of stream-spawning kokanee and seven haplotypes among the tissues of beach-spawning kokanee, with only two haplotypes shared between spawning habitats. Haplotypes occurring in the same clade as F. psychrophilum were the most prevalent (92% of all reads, 60% of all haplotypes), and occurred in kokanee from both spawning habitats (streams and beaches). Subsequent qPCR assays did not find any significant difference in the relative abundance of Flavobacterium AtpA amplicons between samples from the different spawning habitats.

Conclusions

We confirmed the presence of Flavobacterium spp. in both spawning habitats and found weak evidence for increased Flavobacterium diversity in kokanee sampled from stream-spawning sites. However, the quantity of Flavobacterium DNA did not differ between spawning habitats. We recommend further study aimed at quantifying pathogen diversity and abundance in population-level samples of kokanee combined with environmental sampling to better understand the ecology of pathogen infection in this species.

Electronic supplementary material

The online version of this article (doi:10.1186/1756-0500-7-778) contains supplementary material, which is available to authorized users.

Infectious disease, shifting climates, and opportunistic predators: cumulative factors potentially impacting wild salmon declines

Emerging diseases are impacting animals under high-density culture, yet few studies assess their importance to wild populations. Microparasites selected for enhanced virulence in culture settings should be less successful maintaining infectivity in wild populations, as once the host dies, there are limited opportunities to infect new individuals. Instead, moderately virulent microparasites persisting for long periods across multiple environments are of greatest concern. Evolved resistance to endemic microparasites may reduce susceptibilities, but as barriers to microparasite distributions are weakened, and environments become more stressful, unexposed populations may be impacted and pathogenicity enhanced. We provide an overview of the evolutionary and ecological impacts of infectious diseases in wild salmon and suggest ways in which modern technologies can elucidate the microparasites of greatest potential import. We present four case studies that resolve microparasite impacts on adult salmon migration success, impact of river warming on microparasite replication, and infection status on susceptibility to predation. Future health of wild salmon must be considered in a holistic context that includes the cumulative or synergistic impacts of multiple stressors. These approaches will identify populations at greatest risk, critically needed to manage and potentially ameliorate the shifts in current or future trajectories of wild populations.

Inter individual variations of the fish skin microbiota: host genetics basis of mutualism?

The commensal microbiota of fish skin is suspected to provide a protection against opportunist infections. The skin of fish harbors a complex and diverse microbiota that closely interacts with the surrounding water microbial communities. Up to now there is no clear evidence as to whether the host regulates the recruitment of environmental bacteria to build a specific skin microbiota. To address this question, we detected Quantitative Trait Loci (QTL) associated with the abundance of specific skin microbiota bacterial strains in brook charr (Salvelinus fontinalis), combining 16S RNA tagged-amplicon 454 pyrosequencing with genetic linkage analysis. Skin microbiota analysis revealed high inter-individual variation among 86 F2 fish progeny based upon the relative abundance of bacterial operational taxonomic units (OTUs). Out of those OTUs, the pathogenic strain Flavobacterium psychrophilum and the non-pathogenic strain Methylobacterium rhodesianum explained the majority of inter-individual distances. Furthermore, a strong negative correlation was found between Flavobacterium and Methylobacterium, suggesting a mutually competitive relationship. Finally, after considering a total of 266 markers, genetic linkage analysis highlighted three major QTL associated with the abundance of Lysobacter, Rheinheimera and Methylobacterium. All these three genera are known for their beneficial antibacterial activity. Overall, our results provide evidence that host genotype may regulate the abundance of specific genera among their surface microbiota. They also indicate that Lysobacter, Rheinheimera and Methylobacterium are potentially important genera in providing protection against pathogens.

Neutral and selective processes shape MHC gene diversity and expression in stocked brook charr populations (Salvelinus fontinalis)

The capacity of an individual to battle infection is an important fitness determinant in wild vertebrate populations. The major histocompatibility complex (MHC) genes are crucial for a host's adaptive immune system to detect pathogens. However, anthropogenic activities may disrupt natural cycles of co-evolution between hosts and pathogens. In this study, we investigated the dynamic sequence and expression variation of host parasite interactions in brook charr (Salvelinus fontinalis) in a context of past human disturbance via population supplementation from domestic individuals. To do so, we developed a new method to examine selection shaping MHC diversity within and between populations and found a complex interplay between neutral and selective processes that varied between lakes that were investigated. We provided evidence for a lower introgression rate of domestic alleles and found that parasite infection increased with domestic genomic background of individuals. We also documented an association between individual MHC alleles and parasite taxa. Finally, longer cis-regulatory minisatellites were positively correlated with MHC II down-regulation and domestic admixture, suggesting that inadvertent selection during domestication resulted in a lower immune response capacity, through a trade-off between growth and immunity, which explained the negative selection of domestic alleles at least under certain circumstances.

The roles of standing genetic variation and evolutionary history in determining the evolvability of anti-predator strategies

Standing genetic variation and the historical environment in which that variation arises (evolutionary history) are both potentially significant determinants of a population's capacity for evolutionary response to a changing environment. Using the open-ended digital evolution software Avida, we evaluated the relative importance of these two factors in influencing evolutionary trajectories in the face of sudden environmental change. We examined how historical exposure to predation pressures, different levels of genetic variation, and combinations of the two, affected the evolvability of anti-predator strategies and competitive abilities in the presence or absence of threats from new, invasive predator populations. We show that while standing genetic variation plays some role in determining evolutionary responses, evolutionary history has the greater influence on a population's capacity to evolve anti-predator traits, i.e. traits effective against novel predators. This adaptability likely reflects the relative ease of repurposing existing, relevant genes and traits, and the broader potential value of the generation and maintenance of adaptively flexible traits in evolving populations.

Teleost microbiomes: the state of the art in their characterization, manipulation and importance in aquaculture and fisheries

Indigenous microbiota play a critical role in the lives of their vertebrate hosts. In human and mouse models it is increasingly clear that innate and adaptive immunity develop in close concert with the commensal microbiome. Furthermore, several aspects of digestion and nutrient metabolism are governed by intestinal microbiota. Research on teleosts has responded relatively slowly to the introduction of massively parallel sequencing procedures in microbiomics. Nonetheless, progress has been made in biotic and gnotobiotic zebrafish models, defining a core microbiome and describing its role in development. However, microbiome research in other teleost species, especially those important from an aquaculture perspective, has been relatively slow. In this review, we examine progress in teleost microbiome research to date. We discuss teleost microbiomes in health and disease, microbiome ontogeny, prospects for successful microbiome manipulation (especially in an aquaculture setting) and attempt to identify important future research themes. We predict an explosion in research in this sector in line with the increasing global demand for fish protein, and the need to find sustainable approaches to improve aquaculture yield. The reduced cost and increasing ease of next generation sequencing technologies provides the technological backing, and the next 10 years will be an exciting time for teleost microbiome research.

Parasite-mediated selection drives an immunogenetic trade-off in plains zebras (Equus quagga)

Pathogen evasion of the host immune system is a key force driving extreme polymorphism in genes of the major histocompatibility complex (MHC). Although this gene family is well characterized in structure and function, there is still much debate surrounding the mechanisms by which MHC diversity is selectively maintained. Many studies have investigated relationships between MHC variation and specific pathogens, and have found mixed support for and against the hypotheses of heterozygote advantage, frequency-dependent or fluctuating selection. Few, however, have focused on the selective effects of multiple parasite types on host immunogenetic patterns. Here, we examined relationships between variation in the equine MHC gene, ELA-DRA, and both gastrointestinal (GI) and ectoparasitism in plains zebras (Equus quagga). Specific alleles present at opposing population frequencies had antagonistic effects, with rare alleles associated with increased GI parasitism and common alleles with increased tick burdens. These results support a frequency-dependent mechanism, but are also consistent with fluctuating selection. Maladaptive GI parasite 'susceptibility alleles' were reduced in frequency, suggesting that these parasites may play a greater selective role at this locus. Heterozygote advantage, in terms of allele mutational divergence, also predicted decreased GI parasite burden in genotypes with a common allele. We conclude that an immunogenetic trade-off affects resistance/susceptibility to parasites in this system. Because GI and ectoparasites do not directly interact within hosts, our results uniquely show that antagonistic parasite interactions can be indirectly modulated through the host immune system. This study highlights the importance of investigating the role of multiple parasites in shaping patterns of host immunogenetic variation.

Footprints of directional selection in wild Atlantic salmon populations: evidence for parasite-driven evolution?

Mechanisms of host-parasite co-adaptation have long been of interest in evolutionary biology; however, determining the genetic basis of parasite resistance has been challenging. Current advances in genome technologies provide new opportunities for obtaining a genome-scale view of the action of parasite-driven natural selection in wild populations and thus facilitate the search for specific genomic regions underlying inter-population differences in pathogen response. European populations of Atlantic salmon (Salmo salar L.) exhibit natural variance in susceptibility levels to the ectoparasite Gyrodactylus salaris Malmberg 1957, ranging from resistance to extreme susceptibility, and are therefore a good model for studying the evolution of virulence and resistance. However, distinguishing the molecular signatures of genetic drift and environment-associated selection in small populations such as land-locked Atlantic salmon populations presents a challenge, specifically in the search for pathogen-driven selection. We used a novel genome-scan analysis approach that enabled us to i) identify signals of selection in salmon populations affected by varying levels of genetic drift and ii) separate potentially selected loci into the categories of pathogen (G. salaris)-driven selection and selection acting upon other environmental characteristics. A total of 4631 single nucleotide polymorphisms (SNPs) were screened in Atlantic salmon from 12 different northern European populations. We identified three genomic regions potentially affected by parasite-driven selection, as well as three regions presumably affected by salinity-driven directional selection. Functional annotation of candidate SNPs is consistent with the role of the detected genomic regions in immune defence and, implicitly, in osmoregulation. These results provide new insights into the genetic basis of pathogen susceptibility in Atlantic salmon and will enable future searches for the specific genes involved.

Combining molecular evolution and environmental genomics to unravel adaptive processes of MHC class IIB diversity in European minnows (Phoxinus phoxinus)

Host-pathogen interactions are a major evolutionary force promoting local adaptation. Genes of the major histocompatibility complex (MHC) represent unique candidates to investigate evolutionary processes driving local adaptation to parasite communities. The present study aimed at identifying the relative roles of neutral and adaptive processes driving the evolution of MHC class IIB (MHCIIB) genes in natural populations of European minnows (Phoxinus phoxinus). To this end, we isolated and genotyped exon 2 of two MHCIIB gene duplicates (DAB1 and DAB3) and 1'665 amplified fragment length polymorphism (AFLP) markers in nine populations, and characterized local bacterial communities by 16S rDNA barcoding using 454 amplicon sequencing. Both MHCIIB loci exhibited signs of historical balancing selection. Whereas genetic differentiation exceeded that of neutral markers at both loci, the populations' genetic diversities were positively correlated with local pathogen diversities only at DAB3. Overall, our results suggest pathogen-mediated local adaptation in European minnows at both MHCIIB loci. While at DAB1 selection appears to favor different alleles among populations, this is only partially the case in DAB3, which appears to be locally adapted to pathogen communities in terms of genetic diversity. These results provide new insights into the importance of host-pathogen interactions in driving local adaptation in the European minnow, and highlight that the importance of adaptive processes driving MHCIIB gene evolution may differ among duplicates within species, presumably as a consequence of alternative selective regimes or different genomic context. Using next-generation sequencing, the present manuscript identifies the relative roles of neutral and adaptive processes driving the evolution of MHC class IIB (MHCIIB) genes in natural populations of a cyprinid fish: the European minnow (Phoxinus phoxinus). We highlight that the relative importance of neutral versus adaptive processes in shaping immune competence may differ between duplicates as a consequence of alternative selective regimes or different genomic contexts.