Clinal variation in MHC diversity with temperature: evidence for the role of host-pathogen interaction on local adaptation in Atlantic salmon

The interaction of climate, plant, and soil factors drives putative soil fungal pathogen diversity and community structure in dry grasslands

Soil pathogens play important roles in shaping soil microbial diversity and controlling ecosystem functions. Though climate and local environmental factors and their influences on fungal pathogen communities have been examined separately, few studies explore the relative contributions of these factors. This is particularly crucial in eco-fragile regions, which are more sensitive to environmental changes. Herein we investigated the diversity and community structure of putative soil fungal pathogens in cold and dry grasslands on the Tibetan Plateau, using high-throughput sequencing. The results showed that steppe soils had the highest diversity of all pathogens and plant pathogens; contrastingly, meadow soils had the highest animal pathogen diversity. Structural equation modelling revealed that climate, plant, and soil had similar levels of influence on putative soil fungal pathogen diversity, with total effects ranging from 52% to 59% (all p < 0.001), with precipitation exhibiting a stronger direct effect than plant and soil factors. Putative soil fungal pathogen community structure gradually changed with desert, steppe, and meadow, and was primarily controlled by the interactions of climate, plant, and soil factors rather than by distinct factors individually. This finding contrasts with most studies of soil bacterial and fungal community structure, which generally report dominant roles of individual environmental factors.

Multifaceted framework for defining conservation units: An example from Atlantic salmon (Salmo salar) in Canada

Conservation units represent important components of intraspecific diversity that can aid in prioritizing and protecting at-risk populations, while also safeguarding unique diversity that can contribute to species resilience. In Canada, identification and assessments of conservation units is done by the Committee on the Status of Endangered Wildlife in Canada (COSEWIC). COSEWIC can recognize conservation units below the species level (termed "designatable units"; DUs) if the unit has attributes that make it both discrete and evolutionarily significant. There are various ways in which a DU can meet criteria of discreteness and significance, and increasing access to "big data" is providing unprecedented information that can directly inform both criteria. Specifically, the incorporation of genomic data for an increasing number of non-model species is informing more COSEWIC assessments; thus, a repeatable, robust framework is needed for integrating these data into DU characterization. Here, we develop a framework that uses a multifaceted, weight of evidence approach to incorporate multiple data types, including genetic and genomic data, to inform COSEWIC DUs. We apply this framework to delineate DUs of Atlantic salmon (Salmo salar, L.), an economically, culturally, and ecologically significant species, that is also characterized by complex hierarchical population structure. Specifically, we focus on an in-depth example of how our approach was applied to a previously data limited region of northern Canada that was defined by a single large DU. Application of our framework with newly available genetic and genomic data led to subdividing this DU into three new DUs. Although our approach was developed to meet criteria of COSEWIC, it is widely applicable given similarities in the definitions of a conservation unit.

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.

Tenets in Microbial Endocrinology: A New Vista in Teleost Reproduction

Climate vulnerability and induced changes in physico-chemical properties of aquatic environment can bring impairment in metabolism, physiology and reproduction in teleost. Variation in environmental stimuli mainly acts on reproduction by interfering with steroidogenesis, gametogenesis and embryogenesis. The control on reproductive function in captivity is essential for the sustainability of aquaculture production. There are more than 3,000 teleost species across the globe having commercial importance; however, adequate quality and quantity of seed production have been the biggest bottleneck. Probiotics are widely used in aquaculture as a growth promoter, stress tolerance, pathogen inhibition, nutrient digestibility and metabolism, reproductive performance and gamete quality. As the gut microbiota exerts various effects on the intestinal milieu which influences distant organs and pathways, therefore it is considered to be a full-fledged endocrine organ. Researches on Gut-Brain-Gonad axis (GBG axis) and its importance on physiology and reproduction have already been highlighted for higher mammals; however, the study on fish physiology and reproduction is limited. While looking into the paucity of information, we have attempted to review the present status of microbiome and its interaction between the brain and gut. This review will address a process of the microbiome physiological mechanism involved in fish reproduction. The gut microbiota influences the BPG axis through a wide variety of compounds, including neuropeptides, neurotransmitter homologs and transmitters. Currently, research is being conducted to determine the precise process by which gut microbial composition influences brain function in fish. The gut-brain bidirectional interaction can influence brain biochemistry such as GABA, serotonin and tryptophan metabolites which play significant roles in CNS regulation. This review summarizes the fact, how microbes from gut, skin and other parts of the body influence fish reproduction through the Gut-Brain-Gonad axis.

Selection and demography drive range-wide patterns of MHC-DRB variation in mule deer

BACKGROUND

Standing genetic variation is important especially in immune response-related genes because of threats to wild populations like the emergence of novel pathogens. Genetic variation at the major histocompatibility complex (MHC), which is crucial in activating the adaptive immune response, is influenced by both natural selection and historical population demography, and their relative roles can be difficult to disentangle. To provide insight into the influences of natural selection and demography on MHC evolution in large populations, we analyzed geographic patterns of variation at the MHC class II DRB exon 2 locus in mule deer (Odocoileus hemionus) using sequence data collected across their entire broad range.

RESULTS

We identified 31 new MHC-DRB alleles which were phylogenetically similar to other cervid MHC alleles, and one allele that was shared with white-tailed deer (Odocoileus virginianus). We found evidence for selection on the MHC including high dN/dS ratios, positive neutrality tests, deviations from Hardy-Weinberg Equilibrium (HWE) and a stronger pattern of isolation-by-distance (IBD) than expected under neutrality. Historical demography also shaped variation at the MHC, as indicated by similar spatial patterns of variation between MHC and microsatellite loci and a lack of association between genetic variation at either locus type and environmental variables.

CONCLUSIONS

Our results show that both natural selection and historical demography are important drivers in the evolution of the MHC in mule deer and work together to shape functional variation and the evolution of the adaptive immune response in large, well-connected populations.

Microevolution of Pieris butterfly genes involved in host-plant adaptation along a host-plant community cline

Herbivorous insects have evolved counteradaptations to overcome the chemical defenses of their host plants. Several of these counteradaptations have been elucidated at the molecular level, in particular for insects specialized on cruciferous host plants. While the importance of these counteradaptations for host plant colonization is well established, little is known about their microevolutionary dynamics in the field. In particular, it is not known whether and how host plant diversity shapes diversity in insect counteradaptations. In this study, we examine patterns of host plant use and insect counteradaptation in three Pieris butterfly species across Japan. The larvae of these butterflies express nitrile-specifier protein (NSP) and its paralog major allergen (MA) in their gut to overcome the highly diversified glucosinolate-myrosinase defense system of their cruciferous host plants. Pieris napi and Pieris melete colonize wild Brassicaceae whereas Pieris rapae typically uses cultivated Brassica as a host, regardless of the local composition of wild crucifers. As expected, NSP and MA diversity was independent of the local composition of wild Brassicaceae in P. rapae. In contrast, NSP diversity correlated with local host plant diversity in both species that preferred wild Brassicaceae. P. melete and P. napi both revealed two distinct major NSP alleles, which shaped diversity among local populations, albeit with different evolutionary trajectories. In comparison, MA showed no indication for local adaptation. Altogether, MA appeared to be evolutionary more conserved than NSP, suggesting that both genes play different roles in diverting host plant chemical defense.

MHC class I-α population differentiation in a commercial fish, the European sea bass (Dicentrarchus labrax)

Identifying population structuring in highly fecund marine species with high dispersal rates is challenging, but critical for conservation and stock delimitation for fisheries management. European sea bass (Dicentrarchus labrax) is a commercial species of fisheries and aquaculture relevance whose stocks are declining in the North Atlantic, despite management measures to protect them and identifying their fine population structure is needed for managing their exploitation. As for other marine fishes, neutral genetic markers indicate that eastern Atlantic sea bass form a panmictic population and is currently managed as arbitrarily divided stocks. The genes of the major histocompatibility complex (MHC) are key components of the adaptive immune system and ideal candidates to assess fine structuring arising from local selective pressures. We used Illumina sequencing to characterise allelic composition and signatures of selection at the MHC class I-α region of six D. labrax populations across the Atlantic range. We found high allelic diversity driven by positive selection, corresponding to moderate supertype diversity, with 131 alleles clustering into four to eight supertypes, depending on the Bayesian information criterion threshold applied, and a mean number of 13 alleles per individual. Alleles could not be assigned to particular loci, but private alleles allowed us to detect regional genetic structuring not found previously using neutral markers. Our results suggest that MHC markers can be used to detect cryptic population structuring in marine species where neutral markers fail to identify differentiation. This is particularly critical for fisheries management, and of potential use for selective breeding or identifying escapes from sea farms.

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.

Adaptation and cryptic pseudogenization in penguin Toll-like Receptors

Penguins (Sphenisciformes) are an iconic order of flightless, diving seabirds distributed across a large latitudinal range in the Southern Hemisphere. The extensive area over which penguins are endemic is likely to have fostered variation in pathogen pressure, which in turn will have imposed differential selective pressures on the penguin immune system. At the front line of pathogen detection and response, the Toll-like receptors (TLRs) provide insight into host evolution in the face of microbial challenge. TLRs respond to conserved pathogen-associated molecular patterns and are frequently found to be under positive selection, despite retaining specificity for defined agonist classes. We undertook a comparative immunogenetics analysis of TLRs for all penguin species and found evidence of adaptive evolution that was largely restricted to the cell surface-expressed TLRs, with evidence of positive selection at, or near, key agonist-binding sites in TLR1B, TLR4, and TLR5. Intriguingly, TLR15, which is activated by fungal products, appeared to have been pseudogenized multiple times in the Eudyptes spp., but a full-length form was present as a rare haplotype at the population level. However, in vitro analysis revealed that even the full-length form of Eudyptes TLR15 was nonfunctional, indicating an ancestral cryptic pseudogenization prior to its eventual disruption multiple times in the Eudyptes lineage. This unusual pseudogenization event could provide an insight into immune adaptation to fungal pathogens such as Aspergillus, which is responsible for significant mortality in wild and captive bird populations.

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.

Genetic structure of immunologically associated candidate genes suggests arctic rabies variants exert differential selection in arctic fox populations

Patterns of local adaptation can emerge in response to the selective pressures diseases exert on host populations as reflected in increased frequencies of respective, advantageous genotypes. Elucidating patterns of local adaptation enhance our understanding of mechanisms of disease spread and the capacity for species to adapt in context of rapidly changing environments such as the Arctic. Arctic rabies is a lethal disease that largely persists in northern climates and overlaps with the distribution of its natural host, arctic fox. Arctic fox populations display little neutral genetic structure across their North American range, whereas phylogenetically unique arctic rabies variants are restricted in their geographic distributions. It remains unknown if arctic rabies variants impose differential selection upon host populations, nor what role different rabies variants play in the maintenance and spread of this disease. Using a targeted, genotyping-by-sequencing assay, we assessed correlations of arctic fox immunogenetic variation with arctic rabies variants to gain further insight into the epidemiology of this disease. Corroborating past research, we found no neutral genetic structure between sampled regions, but did find moderate immunogenetic structuring between foxes predominated by different arctic rabies variants. F_ST outliers associated with host immunogenetic structure included SNPs within interleukin and Toll-like receptor coding regions (IL12B, IL5, TLR3 and NFKB1); genes known to mediate host responses to rabies. While these data do not necessarily reflect causation, nor a direct link to arctic rabies, the contrasting genetic structure of immunologically associated candidate genes with neutral loci is suggestive of differential selection and patterns of local adaptation in this system. These data are somewhat unexpected given the long-lived nature and dispersal capacities of arctic fox; traits expected to undermine local adaptation. Overall, these data contribute to our understanding of the co-evolutionary relationships between arctic rabies and their primary host and provide data relevant to the management of this disease.

Urbanization processes drive divergence at the major histocompatibility complex in a common waterbird

Urban sprawl is one of the most common landscape alterations occurring worldwide, and there is a growing list of species that are recognised to have adapted to urban life. To be successful, processes of urban colonization by wildlife require a broad spectrum of phenotypic (e.g., behavioural or physiological) adjustments, but evidence for genetic adaptations is much scarcer. One hypothesis proposes that different pathogen-driven selective pressures between urban and non-urban landscapes leads to adaptations in host immune genes. Here, we examined urbanization-related differentiation at the key pathogen-recognition genes of vertebrate adaptive immunity-the major histocompatibility complex (MHC)-in a common waterbird, the Eurasian coot (Fulica atra). Samples were collected from an old urban population (established before the 1950s), a new urban population (established in the 2000s), and two rural populations from central Poland. We found strong significant divergence (as measured with Jost's D) at the MHC class II between the old urban population and the remaining (new urban and rural) populations. Also, there was a moderate, but significant divergence at the MHC between the new urban population and two rural populations, while no divergence was found between the two rural populations. The total number of MHC alleles and the number of private (population-specific) MHC alleles was lower in old urban populations, as compared to the rural ones. These patterns of differentiation at the MHC were not consistent with patterns found for neutral genetic markers (microsatellites), which showed few differences between the populations. Our results indicate that MHC allele composition depended on the level of anthropogenic disturbance and the time which passed since urban colonization, possibly due to the processes of genotype sorting and local adaptation. As such, our study contributes to the understanding of genetic mechanisms associated with urbanization processes in wildlife.

Artificial Rearing of Atlantic Salmon Juveniles for Supportive Breeding Programs Induces Long-Term Effects on Gut Microbiota after Stocking

In supportive breeding programs for wild salmon populations, stocked parr experience higher mortality rates than wild ones. Among other aspects of phenotype, the gut microbiota of artificially raised parr differs from that of wild parr before stocking. Early steps of microbiota ontogeny are tightly dependent upon environmental conditions, both of which exert long-term effects on host physiology. Therefore, our objective was to assess to what extent the resilience capacity of the microbiota of stocked salmon may prevent taxonomic convergence with that of their wild congeners after two months in the same natural environment. Using the 16S SSU rRNA marker gene, we tested the general hypothesis that environmental conditions during the very first steps of microbiota ontogeny imprint a permanent effect on later stages of microbiota recruitment. Our results first showed that gut microbiota composition of stocked and wild parr from the same genetic population, and sharing the same environment, was dependent on the early rearing environment. In contrast, skin microbiota in stocked individuals converged to that of wild individuals. Taxonomic composition and co-occurrence network analyses suggest an impairment of wild bacteria recruitment and a higher instability for the gut microbiota of stocked parr. This study is the first to demonstrate the long-term effect of early microbiota ontogeny in artificial rearing for natural population conservation programs, raising the need to implement microbial ecology.

Clinal and seasonal changes are correlated in Drosophila melanogaster natural populations

Spatial and seasonal variations in the environment are ubiquitous. Environmental heterogeneity can affect natural populations and lead to covariation between environment and allele frequencies. Drosophila melanogaster is known to harbor polymorphisms that change both with latitude and seasons. Identifying the role of selection in driving these changes is not trivial, because nonadaptive processes can cause similar patterns. Given the environment changes in similar ways across seasons and along the latitudinal gradient, one promising approach may be to look for parallelism between clinal and seasonal changes. Here, we test whether there is a genome-wide correlation between clinal and seasonal changes, and whether the pattern is consistent with selection. Allele frequency estimates were obtained from pooled samples from seven different locations along the east coast of the United States, and across seasons within Pennsylvania. We show that there is a genome-wide correlation between clinal and seasonal variations, which cannot be explained by linked selection alone. This pattern is stronger in genomic regions with higher functional content, consistent with natural selection. We derive a way to biologically interpret these correlations and show that around 3.7% of the common, autosomal variants could be under parallel seasonal and spatial selection. Our results highlight the contribution of natural selection in driving fluctuations in allele frequencies in natural fly populations and point to a shared genomic basis to climate adaptation that happens over space and time in D. melanogaster.

Summer Is Coming! Tackling Ocean Warming in Atlantic Salmon Cage Farming

Atlantic salmon (Salmo salar) cage farming has traditionally been located at higher latitudes where cold seawater temperatures favor this practice. However, these regions can be impacted by ocean warming and heat waves that push seawater temperature beyond the thermo-tolerance limits of this species. As more mass mortality events are reported every year due to abnormal sea temperatures, the Atlantic salmon cage aquaculture industry acknowledges the need to adapt to a changing ocean. This paper reviews adult Atlantic salmon thermal tolerance limits, as well as the deleterious eco-physiological consequences of heat stress, with emphasis on how it negatively affects sea cage aquaculture production cycles. Biotechnological solutions targeting the phenotypic plasticity of Atlantic salmon and its genetic diversity, particularly that of its southernmost populations at the limit of its natural zoogeographic distribution, are discussed. Some of these solutions include selective breeding programs, which may play a key role in this quest for a more thermo-tolerant strain of Atlantic salmon that may help the cage aquaculture industry to adapt to climate uncertainties more rapidly, without compromising profitability. Omics technologies and precision breeding, along with cryopreservation breakthroughs, are also part of the available toolbox that includes other solutions that can allow cage farmers to continue to produce Atlantic salmon in the warmer waters of the oceans of tomorrow.

Positive selection on the MHC class II DLA-DQA1 gene in golden jackals (Canis aureus) from their recent expansion range in Europe and its effect on their body mass index

Background

In Europe, golden jackals (Canis aureus) have been expanding their range out of the southern and southeastern Balkans towards central Europe continually since the 1960s. Here, we investigated the level of functional diversity at the MHC class II DLA-DQA1 exon 2 in golden jackal populations from Bulgaria, Serbia, and Hungary. Specifically, we tested for positive selection on and geographic variation at that locus due to adaptation to supposedly regionally varying pathogenic landscapes. To test for potential fitness effects of different protein variants on individual body condition, we used linear modeling of individual body mass indexes (bmi) and accounted for possible age, sex, geographical, and climatic effects. The latter approach was performed, however, only on Serbian individuals with appropriate data.

Results

Only three different DLA-DQA1 alleles were detected, all coding for different amino-acid sequences. The neutrality tests revealed no significant but positive values; there was no signal of spatial structuring and no deviation from the Hardy–Weinberg equilibrium across the studied range of expansion. However, we found a signal of trans-species polymorphism and significant test results for positive selection on three codons. Our information-theory based linear modeling results indicated an effect of ambient temperature on the occurrence of individual DLA-DQA1 genotypes in individuals from across the studied expansion range, independent from geographical position. Our linear modeling results of individual bmi values indicated that yearlings homozygous for DLA-DQA1*03001 reached values typical for adults contrary to yearlings carrying other genotypes (protein combinations). This suggested better growth rates and thus a possible fitness advantage of yearlings homozygous for DLA-DQA1*03001.

Conclusions

Our results indicate a demographic (stochastic) signal of reduced DLA-DQA1 exon 2 variation, in line with the documented historical demographic bottleneck. At the same time, however, allelic variation was also affected by positive selection and adaptation to varying ambient temperature, supposedly reflecting geographic variation in the pathogenic landscape. Moreover, an allele effect on body mass index values of yearlings suggested differential fitness associated with growth rates. Overall, a combination of a stochastic effect and positive selection has shaped and is still shaping the variation at the studied MHC locus.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12862-021-01856-z.

MHC Reflects Fine-Scale Habitat Structure in White-Tailed Eagles, Haliaeetus albicilla

The major histocompatibility complex (MHC) genes code for key immune receptors responsible for recognition of intra- and extracellular pathogens (MHC class I and class II, respectively). It was hypothesized that MHC polymorphism can be maintained via fluctuating selection resulting from between-habitat variation in pathogen regimes. We examined associations between MHC class I and class II genes and habitat structure in an apex avian predator, the white-tailed eagle, Haliaeetus albicilla. We genotyped MHC class I and class II genes in ca. 150 white-tailed eagle chicks from nearly 100 nesting territories distributed across 3 distinct populations in Poland. Habitat structure was quantified at the level of foraging territories and directly at the nest sites. We found strong support for associations of habitat traits with diversity and allelic composition at the MHC class II. Forest area within territory and forest productivity were identified as the major habitat predictors of MHC class II polymorphism, whereas other habitat traits (distance to nearest open water, grassland, and water area within territory or understory presence) showed fewer associations with class II alleles. In contrast, there was little support for associations between MHC class I genes and habitat structure. All significant associations were apparent at the within-population level rather than between populations. Our results suggest that extracellular (rather than intracellular) pathogens may exert much stronger selective pressure on the white-tailed eagle. Associations of habitat structure with MHC class II may reflect fluctuating (balancing) selection, which maintains MHC diversity within populations.

The role of a mechanistic host in maintaining arctic rabies variant distributions: Assessment of functional genetic diversity in Alaskan red fox (Vulpes vulpes)

Populations are exposed to different types and strains of pathogens across heterogeneous landscapes, where local interactions between host and pathogen may present reciprocal selective forces leading to correlated patterns of spatial genetic structure. Understanding these coevolutionary patterns provides insight into mechanisms of disease spread and maintenance. Arctic rabies (AR) is a lethal disease with viral variants that occupy distinct geographic distributions across North America and Europe. Red fox (Vulpes vulpes) are a highly susceptible AR host, whose range overlaps both geographically distinct AR strains and regions where AR is absent. It is unclear if genetic structure exists among red fox populations relative to the presence/absence of AR or the spatial distribution of AR variants. Acquiring these data may enhance our understanding of the role of red fox in AR maintenance/spread and inform disease control strategies. Using a genotyping-by-sequencing assay targeting 116 genomic regions of immunogenetic relevance, we screened for sequence variation among red fox populations from Alaska and an outgroup from Ontario, including areas with different AR variants, and regions where the disease was absent. Presumed neutral SNP data from the assay found negligible levels of neutral genetic structure among Alaskan populations. The immunogenetically-associated data identified 30 outlier SNPs supporting weak to moderate genetic structure between regions with and without AR in Alaska. The outliers included SNPs with the potential to cause missense mutations within several toll-like receptor genes that have been associated with AR outcome. In contrast, there was a lack of genetic structure between regions with different AR variants. Combined, we interpret these data to suggest red fox populations respond differently to the presence of AR, but not AR variants. This research increases our understanding of AR dynamics in the Arctic, where host/disease patterns are undergoing flux in a rapidly changing Arctic landscape, including the continued northward expansion of red fox into regions previously predominated by the arctic fox (Vulpes lagopus).

Genetic drift and bottleneck do not influence diversity in Toll-like receptor genes at a small spatial scale in a Himalayan passerine

Genetic diversity is important for long-term viability of a population. Low genetic diversity reduces persistence and survival of populations and increases susceptibility to diseases. Comparisons of the neutral markers with functional loci such as immune genes [Toll-like receptors; TLR] can provide useful insights into evolutionary potential of a species and how the diversity of pathogens and selection pressures on their hosts are directly linked to their environment. In this study, we compare genetic diversity in neutral (eleven microsatellite loci) and adaptive (seven TLR loci) loci to determine genetic variation in a nonmigratory western Himalayan passerine, the black-throated tit (Aegithalos concinnus), distributed across an elevation gradient with varying degree of pathogen-mediated selection pressure. We further compare the diversity in TLR loci with a high-elevation sister species, the white-throated tit (Aegithalos niveogularis). Our results indicate a lack of population genetic structure in the black-throated tit and signatures of a past bottleneck. In contrast, we found high diversity in TLR loci and locus-specific (TLR7) signatures of pathogen-mediated selection, which was comparable to diversity in the white-throated tit. Levels of diversity at TLR5 locus corresponded very closely with neutral microsatellite variation. We found evidence of positive selection in TLR1LA, TLR5, and TLR7 loci highlighting the importance in pathogen recognition. Our finding demonstrates that reduction in neutral variation does not necessarily lead to reduction in functional genetic diversity and probably helps in revival of population in a widespread species.

A continuously changing selective context on microbial communities associated with fish, from egg to fork

Fast increase of fish aquaculture production to meet consumer demands is accompanied by important ecological concerns such as disease outbreaks. Meanwhile, food waste is an important concern with fish products since they are highly perishable. Recent aquaculture and fish product microbiology, and more recently, microbiota research, paved the way to a highly integrated approach to understand complex relationships between host fish, product and their associated microbial communities at health/disease and preservation/spoilage frontiers. Microbial manipulation strategies are increasingly validated as promising tools either to replace or to complement traditional veterinary and preservation methods. In this review, we consider evolutionary forces driving fish microbiota assembly, in particular the changes in the selective context along the production chain. We summarize the current knowledge concerning factors governing assembly and dynamics of fish hosts and food microbial communities. Then, we discuss the current microbial community manipulation strategies from an evolutionary standpoint to provide a perspective on the potential for risks, conflict and opportunities. Finally, we conclude that to harness evolutionary forces in the development of sustainable microbiota manipulation applications in the fish industry, an integrated knowledge of the controlling abiotic and especially biotic factors is required.

Evidence of Pathogen-Induced Immunogenetic Selection across the Large Geographic Range of a Wild Seabird

Over evolutionary time, pathogen challenge shapes the immune phenotype of the host to better respond to an incipient threat. The extent and direction of this selection pressure depend on the local pathogen composition, which is in turn determined by biotic and abiotic features of the environment. However, little is known about adaptation to local pathogen threats in wild animals. The Gentoo penguin (Pygoscelis papua) is a species complex that lends itself to the study of immune adaptation because of its circumpolar distribution over a large latitudinal range, with little or no admixture between different clades. In this study, we examine the diversity in a key family of innate immune genes-the Toll-like receptors (TLRs)-across the range of the Gentoo penguin. The three TLRs that we investigated present varying levels of diversity, with TLR4 and TLR5 greatly exceeding the diversity of TLR7. We present evidence of positive selection in TLR4 and TLR5, which points to pathogen-driven adaptation to the local pathogen milieu. Finally, we demonstrate that two positively selected cosegregating sites in TLR5 are sufficient to alter the responsiveness of the receptor to its bacterial ligand, flagellin. Taken together, these results suggest that Gentoo penguins have experienced distinct pathogen-driven selection pressures in different environments, which may be important given the role of the Gentoo penguin as a sentinel species in some of the world's most rapidly changing environments.

Northern areas as refugia for temperate species under current climate warming: Atlantic salmon (Salmo salar L.) as a model in northern Europe

In this work, patterns of geographical genetic diversity in Atlantic salmon Salmo salar were studied across the whole Atlantic Arc; whether these patterns (and thus genetic population structure) were affected by water temperatures was also evaluated. Salmo salar populations were characterized using microsatellite loci and then analysed with reference to ocean surface temperature data from across the region. Analysis showed the presence of a latitudinal cline of genetic variability (higher in northern areas) and water temperatures (sea surface temperatures) determining genetic population structure (the latter in combination with genetic drift in southern populations). Under the current global change scenario, northern areas of Europe would constitute refugia for diversity in the future. This is effectively the inverse of what appears to have happened in glacial refugia during the last glacial maximum. From this perspective, the still abundant and large northern populations S. salar should be considered as precious as the small almost relict southern ones and given appropriate protection. Careful management of the species, coordinated across countries and latitudes, is needed in order to avoid its extinction in Europe.

Population genetic analysis reveals a geographically limited transition zone between two genetically distinct Atlantic salmon lineages in Norway

Atlantic salmon is characterized by a high degree of population genetic structure throughout its native range. However, while populations inhabiting rivers in Norway and Russia make up a significant proportion of salmon in the Atlantic, thus far, genetic studies in this region have only encompassed low to modest numbers of populations. Here, we provide the first "in-depth" investigation of population genetic structuring in the species in this region. Analysis of 18 microsatellites on >9,000 fish from 115 rivers revealed highly significant population genetic structure throughout, following a hierarchical pattern. The highest and clearest level of division separated populations north and south of the Lofoten region in northern Norway. In this region, only a few populations displayed intermediate genetic profiles, strongly indicating a geographically limited transition zone. This was further supported by a dedicated cline analysis. Population genetic structure was also characterized by a pattern of isolation by distance. A decline in overall genetic diversity was observed from the south to the north, and two of the microsatellites showed a clear decrease in number of alleles across the observed transition zone. Together, these analyses support results from previous studies, that salmon in Norway originate from two main genetic lineages, one from the Barents-White Sea refugium that recolonized northern Norwegian and adjacent Russian rivers, and one from the eastern Atlantic that recolonized the rest of Norway. Furthermore, our results indicate that local conditions in the limited geographic transition zone between the two observed lineages, characterized by open coastline with no obvious barriers to gene flow, are strong enough to maintain the genetic differentiation between them.

Divergent transcriptomic signatures in response to salinity exposure in two populations of an estuarine fish

In estuary and coastal systems, human demand for freshwater, climate change-driven precipitation variability, and extreme weather impact salinity levels, reducing connectivity between mesohaline coastal fish populations and potentially contributing to genomic divergence. We examined gill transcriptome responses to salinity in wild-caught juveniles from two populations of Sacramento splittail (Pogonichthys macrolepidotus), a species of conservation concern that is endemic to the San Francisco Estuary, USA, and the lower reaches of its tributaries. Recent extreme droughts have led to salinities above the tolerance limits for this species, creating a migration barrier between these populations, which potentially contributed to population divergence. We identified transcripts involved in a conserved response to salinity; however, the more salinity-tolerant San Pablo population had greater transcriptome plasticity (3.6-fold more transcripts responded than the Central Valley population) and a response consistent with gill remodeling after 168 hr of exposure to elevated salinity. The reorganization of the gill in response to changing osmotic gradients is a process critical for acclimation and would facilitate enhanced salinity tolerance. We detected an upregulation of receptors that control the Wnt (wingless-type) cell signaling pathway that may be required for an adaptive response to increases in salinity, patterns not observed in the relatively salinity-sensitive Central Valley population. We detected 62 single nucleotide polymorphisms (SNPs) in coding regions of 26 transcripts that differed between the populations. Eight transcripts that contained SNPs were associated with immune responses, highlighting the importance of diversity in immune gene sequences as a defining characteristic of genomic divergence between these populations. Our data demonstrate that these populations have divergent transcriptomic responses to salinity, which is consistent with observed physiological differences in salinity tolerance.

Molecular and Structural Characterization of MHC Class II β Genes Reveals High Diversity in the Cold-Adapted Icefish Chionodraco hamatus

This study reports the presence of two distinct MHC class II β genes in the Antarctic icefish Chionodraco hamatus, belonging to the classical (ChhaDAB) and nonclassical (ChhaDBB) evolutionary lineages. By the application of targeted sequencing approach, a remarkable molecular diversity in the exon 2 sequence of the highly expressed gene ChhaDAB has been observed, resulting in an estimate of 92 different variants translated in 87 different peptides from 54 analysed icefish individuals. A highly conservative estimate, based on a 95% sequence identity threshold clustering, translate this variability in 41 different peptide clusters belonging to four different clades and showing the signature of different kinds of selection. In stark contrast, the poorly expressed ChhaDBB gene displayed a very low level of molecular diversity within exon 2, in agreement with expectations for a nonclassical MHC class II β gene.

Natural Selection and Origin of a Melanistic Allele in North American Gray Wolves

Pigmentation is often used to understand how natural selection affects genetic variation in wild populations since it can have a simple genetic basis, and can affect a variety of fitness-related traits (e.g., camouflage, thermoregulation, and sexual display). In gray wolves, the K locus, a β-defensin gene, causes black coat color via a dominantly inherited KB allele. The allele is derived from dog-wolf hybridization and is at high frequency in North American wolf populations. We designed a DNA capture array to probe the geographic origin, age, and number of introgression events of the KB allele in a panel of 331 wolves and 20 dogs. We found low diversity in KB, but not ancestral ky, wolf haplotypes consistent with a selective sweep of the black haplotype across North America. Further, North American wolf KB haplotypes are monophyletic, suggesting that a single adaptive introgression from dogs to wolves most likely occurred in the Northwest Territories or Yukon. We use a new analytical approach to date the origin of the KB allele in Yukon wolves to between 1,598 and 7,248 years ago, suggesting that introgression with early Native American dogs was the source. Using population genetic simulations, we show that the K locus is undergoing natural selection in four wolf populations. We find evidence for balancing selection, specifically in Yellowstone wolves, which could be a result of selection for enhanced immunity in response to distemper. With these data, we demonstrate how the spread of an adaptive variant may have occurred across a species' geographic range.

Rapid seasonal evolution in innate immunity of wild Drosophila melanogaster

Understanding the rate of evolutionary change and the genetic architecture that facilitates rapid adaptation is a current challenge in evolutionary biology. Comparative studies show that genes with immune function are among the most rapidly evolving genes across a range of taxa. Here, we use immune defence in natural populations of Drosophila melanogaster to understand the rate of evolution in natural populations and the genetics underlying rapid change. We probed the immune system using the natural pathogens Enterococcus faecalis and Providencia rettgeri to measure post-infection survival and bacterial load of wild D. melanogaster populations collected across seasonal time along a latitudinal transect along eastern North America (Massachusetts, Pennsylvania and Virginia). There are pronounced and repeatable changes in the immune response over the approximately 10 generations between spring and autumn collections, with a significant but less distinct difference observed among geographical locations. Genes with known immune function are not enriched among alleles that cycle with seasonal time, but the immune function of a subset of seasonally cycling alleles in immune genes was tested using reconstructed outbred populations. We find that flies containing seasonal alleles in Thioester-containing protein 3 (Tep3) have different functional responses to infection and that epistatic interactions among seasonal Tep3 and Drosomycin-like 6 (Dro6) alleles underlie the immune phenotypes observed in natural populations. This rapid, cyclic response to seasonal environmental pressure broadens our understanding of the complex ecological and genetic interactions determining the evolution of immune defence in natural populations.

Chromosomal fusion and life history-associated genomic variation contribute to within-river local adaptation of Atlantic salmon

Chromosomal inversions have been implicated in facilitating adaptation in the face of high levels of gene flow, but whether chromosomal fusions also have similar potential remains poorly understood. Atlantic salmon are usually characterized by population structure at multiple spatial scales; however, this is not the case for tributaries of the Miramichi River in North America. To resolve genetic relationships between populations in this system and the potential for known chromosomal fusions to contribute to adaptation, we genotyped 728 juvenile salmon using a 50 K SNP array. Consistent with previous work, we report extremely weak overall population structuring (Global F_ST  = 0.004) and failed to support hierarchical structure between the river's two main branches. We provide the first genomic characterization of a previously described polymorphic fusion between chromosomes 8 and 29. Fusion genomic characteristics included high LD, reduced heterozygosity in the fused homokaryotes, and strong divergence between the fused and the unfused rearrangement. Population structure based on fusion karyotype was five times stronger than neutral variation (F_ST  = 0.019), and the frequency of the fusion was associated with summer precipitation supporting a hypothesis that this rearrangement may contribute local adaptation despite weak neutral differentiation. Additionally, both outlier variation among populations and a polygenic framework for characterizing adaptive variation in relation to climate identified a 250-Kb region of chromosome 9, including the gene six6 that has previously been linked to age-at-maturity and run-timing for this species. Overall, our results indicate that adaptive processes, independent of major river branching, are more important than neutral processes for structuring these populations.

Structural and compositional mismatch between captive and wild Atlantic salmon (Salmo salar) parrs' gut microbiota highlights the relevance of integrating molecular ecology for management and conservation methods

Stocking methods are used in the Province of Quebec to restore Salmo salar populations. However, Atlantic salmon stocked juveniles show higher mortality rates than wild ones when introduced into nature. Hatchery environment, which greatly differs from the natural environment, is identified as the main driver of the phenotypic mismatch between captive and wild parrs. The latter is also suspected to impact the gut microbiota composition, which can be associated with essential metabolic functions for their host. We hypothesized that hatchery-raised parrs potentially recruit gut microbial communities that are different from those recruited in the wild. This study evaluated the impacts of artificial rearing on gut microbiota composition in 0+ parrs meant for stocking in two distinct Canadian rivers: Rimouski and Malbaie (Quebec, Canada). Striking differences between hatchery and wild-born parrs' gut microbiota suggest that microbiota could be another factor that could impact their survival in the targeted river, because the microbiome is narrowly related to host physiology. For instance, major commensals belonging to Enterobacteriaceae and Clostridiacea from wild parrs' gut microbiota were substituted in captive parrs by lactic acid bacteria from the Lactobacillaceae family. Overall, captive parrs host a generalist bacterial community whereas wild parrs' microbiota is much more specialized. This is the very first study demonstrating extensive impact of captive rearing on intestinal microbiota composition in Atlantic salmon intended for wild population stocking. Our results strongly suggest the need to implement microbial ecology concepts into conservation management of endangered salmon stocks supplemented with hatchery-reared parrs.

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.

Genomewide evidence of environmentally mediated secondary contact of European green crab (Carcinus maenas) lineages in eastern North America

Genetic-environment associations are increasingly revealed through population genomic data and can occur through a number of processes, including secondary contact, divergent natural selection, or isolation by distance. Here, we investigate the influence of the environment, including seasonal temperature and salinity, on the population structure of the invasive European green crab (Carcinus maenas) in eastern North America. Green crab populations in eastern North America are associated with two independent invasions, previously shown to consist of distinct northern and southern ecotypes, with a contact zone in southern Nova Scotia, Canada. Using a RAD-seq panel of 9,137 genomewide SNPs, we detected 41 SNPs (0.49%) whose allele frequencies were highly correlated with environmental data. A principal components analysis of 25 environmental variables differentiated populations into northern, southern, and admixed sites in concordance with the observed genomic spatial structure. Furthermore, a spatial principal components analysis conducted on genomic and geographic data revealed a high degree of global structure (p < .0001) partitioning a northern and southern ecotype. Redundancy and partial redundancy analyses revealed that among the environmental variables tested, winter sea surface temperature had the strongest association with spatial structuring, suggesting that it is an important factor defining range and expansion limits of each ecotype. Understanding environmental thresholds associated with intraspecific diversity will facilitate the ability to manage current and predict future distributions of this aquatic invasive species.

SNP hot-spots in the clam parasite QPX

Background

Quahog Parasite Unknown (QPX) is an opportunistic protistan pathogen of the clam Mercenaria mercenaria. Infections with QPX have caused significant economic losses in the Northeastern United States. Previous research demonstrated a geographic gradient for disease prevalence and intensity, but little information is available on the genetic diversity of the parasite throughout its distribution range. Also, QPX virulence factors are not well understood. This study addresses the occurrence of QPX genetic variants with a particular focus on functions involved in virulence and adaptation to environmental conditions.

Results

Analyses were performed using transcriptome-wide single-nucleotide polymorphism (SNP) of four QPX isolates cultured from infected clams collected from disparate locations along the Northeastern United States. For contig assembly and mapping, two different genome builds and four transcriptomes of the parasite were examined. Genomic variants appeared at a differential rate relative to sequenced transcripts at 20.18 and 22.55% occurrence under 1000 base pairs upstream and downstream protein domains respectively and at 57.26% rate in protein domain coding sequences. QPX strains shared 30.50% of the mutations and exhibited a preferential nucleotide substitution towards thymine. Sequence identity suggested relatedness between different QPX strains, with the parasite being possibly introduced to Virginia from the Massachusetts region during clam trading, while QPX could have been naturally present in New York. Diversity in virulence, temperature, and salinity domains suggested a common variability between strains, but with a preferential higher variation in local adaptation genes. This could explain differences in disease prevalence noted in different regions. Overall, the results supported views that this opportunistic parasite might be able to adapt to varying environmental conditions.

Conclusion

Relatedness and mutations between the four QPX strains suggested that variability in environmental-related functions favors parasite survival, potentially promoting resilience against stressful conditions. These findings are in agreement with the widespread presence of QPX in the environment. Although QPX levels are enzootic in most areas, an increase in disease outbreaks were often associated with seasonal changes in environmental conditions. A selection mediated by the parasitic life of QPX remains possible, but the effect of the environment on the biology of the parasite appears more obvious.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4866-8) contains supplementary material, which is available to authorized users.

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.

Pathogen richness and abundance predict patterns of adaptive major histocompatibility complex variation in insular amphibians

The identification of the factors responsible for genetic variation and differentiation at adaptive loci can provide important insights into the evolutionary process and is crucial for the effective management of threatened species. We studied the impact of environmental viral richness and abundance on functional diversity and differentiation of the MHC class Ia locus in populations of the black-spotted pond frog (Pelophylax nigromaculatus), an IUCN-listed species, on 24 land-bridge islands of the Zhoushan Archipelago and three nearby mainland sites. We found a high proportion of private MHC alleles in mainland and insular populations, corresponding to 32 distinct functional supertypes, and strong positive selection on MHC antigen-binding sites in all populations. Viral pathogen diversity and abundance were reduced at island sites relative to the mainland, and islands housed distinctive viral communities. Standardized MHC diversity at island sites exceeded that found at neutral microsatellites, and the representation of key functional supertypes was positively correlated with the abundance of specific viruses in the environment (Frog virus 3 and Ambystoma tigrinum virus). These results indicate that pathogen-driven diversifying selection can play an important role in maintaining functionally important MHC variation following island isolation, highlighting the importance of considering functionally important genetic variation and host-pathogen associations in conservation planning and management.

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.

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.

Testing mate choice and overdominance at MH in natural families of Atlantic salmon Salmo salar

This study aimed to test mate choice and selection during early life stages on major histocompatibility (MH) genotype in natural families of Atlantic salmon Salmo salar spawners and juveniles, using nine microsatellites to reconstruct families, one microsatellite linked to an MH class I gene and one minisatellite linked to an MH class II gene. MH-based mate choice was only detected for the class I locus on the first year, with lower expected heterozygosity in the offspring of actually mated pairs than predicted under random mating. The genotype frequencies of MH-linked loci observed in the juveniles were compared with frequencies expected from Mendelian inheritance of parental alleles to detect selection during early life stages. No selection was detected on the locus linked to class I gene. For the locus linked to class II gene, observed heterozygosity was higher than expected in the first year and lower in the second year, suggesting overdominance and underdominance, respectively. Within family, juveniles' body size was linked to heterozygosity at the same locus, with longer heterozygotes in the first year and longer homozygotes in the second year. Selection therefore seems to differ from one locus to the other and from year to year.

The evolution of the major histocompatibility complex in upstream versus downstream river populations of the longnose dace

Populations in upstream versus downstream river locations can be exposed to vastly different environmental and ecological conditions and can thus harbor different genetic resources due to selection and neutral processes. An interesting question is how upstream–downstream directionality in rivers affects the evolution of immune response genes. We used next‐generation amplicon sequencing to identify eight alleles of the major histocompatibility complex (MHC) class II β exon 2 in the cyprinid longnose dace (Rhinichthys cataractae) from three rivers in Alberta, upstream and downstream of municipal and agricultural areas along contaminant gradients. We used these data to test for directional and balancing selection on the MHC. We also genotyped microsatellite loci to examine neutral population processes in this system. We found evidence for balancing selection on the MHC in the form of increased nonsynonymous variation relative to neutral expectations, and selection occurred at more amino acid residues upstream than downstream in two rivers. We found this pattern despite no population structure or isolation by distance, based on microsatellite data, at these sites. Overall, our results suggest that MHC evolution is driven by upstream–downstream directionality in fish inhabiting this system.

Outlier analyses to test for local adaptation to breeding grounds in a migratory arctic seabird

Investigating the extent (or the existence) of local adaptation is crucial to understanding how populations adapt. When experiments or fitness measurements are difficult or impossible to perform in natural populations, genomic techniques allow us to investigate local adaptation through the comparison of allele frequencies and outlier loci along environmental clines. The thick‐billed murre (Uria lomvia) is a highly philopatric colonial arctic seabird that occupies a significant environmental gradient, shows marked phenotypic differences among colonies, and has large effective population sizes. To test whether thick‐billed murres from five colonies along the eastern Canadian Arctic coast show genomic signatures of local adaptation to their breeding grounds, we analyzed geographic variation in genome‐wide markers mapped to a newly assembled thick‐billed murre reference genome. We used outlier analyses to detect loci putatively under selection, and clustering analyses to investigate patterns of differentiation based on 2220 genomewide single nucleotide polymorphisms (SNPs) and 137 outlier SNPs. We found no evidence of population structure among colonies using all loci but found population structure based on outliers only, where birds from the two northernmost colonies (Minarets and Prince Leopold) grouped with birds from the southernmost colony (Gannet), and birds from Coats and Akpatok were distinct from all other colonies. Although results from our analyses did not support local adaptation along the latitudinal cline of breeding colonies, outlier loci grouped birds from different colonies according to their non‐breeding distributions, suggesting that outliers may be informative about adaptation and/or demographic connectivity associated with their migration patterns or nonbreeding grounds.

Variation in Major Histocompatibility Complex diversity in invasive cane toad populations

Context The cane toad (Rhinella marina), a native species of central and southern America, was introduced to Australia in 1935 as a biocontrol agent after a complex history of prior introductions. The population rapidly expanded and has since spread through much of the Australian landmass, with severe impacts on the endemic wildlife, primarily via toxicity to predators. The invasion process has taken its toll on the cane toad, with changes in the immunological capacity across the Australian invasive population. Aims To investigate the immunogenetic underpinnings of these changes, we studied the diversity of the Major Histocompatiblity Complex (MHC) genes in introduced cane toad populations. Methods We studied the diversity of two MHC genes (the classical class I UA locus and a class II DAB locus) and compared these with neutral microsatellite markers in toads from the Australian site of introduction and the Australian invasion front. We also included toads from Hawai’i, the original source of the Australian toads, to infer founder effect. Key results Diversity across all markers was low across Australian and Hawai’ian samples, consistent with a reduction in genetic diversity through multiple founder effects during the course of the successive translocations. In Australia, allelic diversity at the microsatellite markers and the UA locus was reduced at the invasion front, whereas all three alleles at the DAB locus were maintained in the invasion-front toads. Conclusions Loss of allelic diversity observed at the microsatellite markers and the UA locus could be the result of drift and bottlenecking along the invasion process, however, the persistence of DAB diversity warrants further investigation to disentangle the evolutionary forces influencing this locus. Implications Through the use of different molecular markers, we provide a preliminary description of the adaptive genetic processes occurring in this invasive population. The extremely limited MHC diversity may represent low immunogenetic competence across the Australian population, which could be exploited for invasive species management.

Coloniality and migration are related to selection on MHC genes in birds

The major histocompatibility complex (MHC) plays a key role in pathogen recognition as a part of the vertebrate adaptive immune system. The great diversity of MHC genes in natural populations is maintained by different forms of balancing selection and its strength should correlate with the diversity of pathogens to which a population is exposed and the rate of exposure. Despite this prediction, little is known about how life-history characteristics affect selection at the MHC. Here, we examined whether the strength of balancing selection on MHC class II genes in birds (as measured with nonsynonymous nucleotide substitutions, dN) was related to their social or migratory behavior, two life-history characteristics correlated with pathogen exposure. Our comparative analysis indicated that the rate of nonsynonymous substitutions was higher in colonial and migratory species than solitary and resident species, suggesting that the strength of balancing selection increases with coloniality and migratory status. These patterns could be attributed to: (1) elevated transmission rates of pathogens in species that breed in dense aggregations, or (2) exposure to a more diverse fauna of pathogens and parasites in migratory species. Our study suggests that differences in social structure and basic ecological traits influence MHC diversity in natural vertebrate populations.