Signals of heterogeneous selection at an MHC locus in geographically proximate ecotypes of sockeye salmon

Recent selection created distinctive variability patterns on MHC class II loci in three dolphin species from the Mediterranean Sea

The major histocompatibility complex (MHC) includes highly polymorphic genes involved in antigen presentation, which is crucial for adaptive immune response. They represent fitness related genetic markers particularly informative for populations exposed to environmental challenges. Here we analyse the diversity and evolutionary traits of MHC class II DQA and DQB genes in the dolphins Stenella coeruleoalba and Grampus griseus from the Mediterranean Sea. We found substantial nucleotide and functional diversity, as well as strong evidence of balancing selection indicated by allele and supertype frequencies, Tajima's D statistics and dN/dS tests. The Risso's dolphin, considered the least abundant in the region, showed the effect of divergent allele advantage at the nucleotide and functional-peptide levels. An outstanding polymorphism was found in the striped dolphin, particularly intriguing in the DQA gene where the Ewens-Watterson test detected a selection sweep that occurred in recent history. We hypothesize that morbillivirus, which has recurrently invaded Mediterranean populations over the last decades, exerted the detected selective pressure.

Conserved islands of divergence associated with adaptive variation in sockeye salmon are maintained by multiple mechanisms

Local adaptation is facilitated by loci clustered in relatively few regions of the genome, termed genomic islands of divergence. The mechanisms that create and maintain these islands and how they contribute to adaptive divergence is an active research topic. Here, we use sockeye salmon as a model to investigate both the mechanisms responsible for creating islands of divergence and the patterns of differentiation at these islands. Previous research suggested that multiple islands contributed to adaptive radiation of sockeye salmon. However, the low-density genomic methods used by these studies made it difficult to fully elucidate the mechanisms responsible for islands and connect genotypes to adaptive variation. We used whole genome resequencing to genotype millions of loci to investigate patterns of genetic variation at islands and the mechanisms that potentially created them. We discovered 64 islands, including 16 clustered in four genomic regions shared between two isolated populations. Characterisation of these four regions suggested that three were likely created by structural variation, while one was created by processes not involving structural variation. All four regions were small (< 600 kb), suggesting low recombination regions do not have to span megabases to be important for adaptive divergence. Differentiation at islands was not consistently associated with established population attributes. In sum, the landscape of adaptive divergence and the mechanisms that create it are complex; this complexity likely helps to facilitate fine-scale local adaptation unique to each population.

The sockeye salmon genome, transcriptome, and analyses identifying population defining regions of the genome

Sockeye salmon (Oncorhynchus nerka) is a commercially and culturally important species to the people that live along the northern Pacific Ocean coast. There are two main sockeye salmon ecotypes—the ocean-going (anadromous) ecotype and the fresh-water ecotype known as kokanee. The goal of this study was to better understand the population structure of sockeye salmon and identify possible genomic differences among populations and between the two ecotypes. In pursuit of this goal, we generated the first reference sockeye salmon genome assembly and an RNA-seq transcriptome data set to better annotate features of the assembly. Resequenced whole-genomes of 140 sockeye salmon and kokanee were analyzed to understand population structure and identify genomic differences between ecotypes. Three distinct geographic and genetic groups were identified from analyses of the resequencing data. Nucleotide variants in an immunoglobulin heavy chain variable gene cluster on chromosome 26 were found to differentiate the northwestern group from the southern and upper Columbia River groups. Several candidate genes were found to be associated with the kokanee ecotype. Many of these genes were related to ammonia tolerance or vision. Finally, the sex chromosomes of this species were better characterized, and an alternative sex-determination mechanism was identified in a subset of upper Columbia River kokanee.

Climate vulnerability assessment for Pacific salmon and steelhead in the California Current Large Marine Ecosystem

Major ecological realignments are already occurring in response to climate change. To be successful, conservation strategies now need to account for geographical patterns in traits sensitive to climate change, as well as climate threats to species-level diversity. As part of an effort to provide such information, we conducted a climate vulnerability assessment that included all anadromous Pacific salmon and steelhead (Oncorhynchus spp.) population units listed under the U.S. Endangered Species Act. Using an expert-based scoring system, we ranked 20 attributes for the 28 listed units and 5 additional units. Attributes captured biological sensitivity, or the strength of linkages between each listing unit and the present climate; climate exposure, or the magnitude of projected change in local environmental conditions; and adaptive capacity, or the ability to modify phenotypes to cope with new climatic conditions. Each listing unit was then assigned one of four vulnerability categories. Units ranked most vulnerable overall were Chinook (O. tshawytscha) in the California Central Valley, coho (O. kisutch) in California and southern Oregon, sockeye (O. nerka) in the Snake River Basin, and spring-run Chinook in the interior Columbia and Willamette River Basins. We identified units with similar vulnerability profiles using a hierarchical cluster analysis. Life history characteristics, especially freshwater and estuary residence times, interplayed with gradations in exposure from south to north and from coastal to interior regions to generate landscape-level patterns within each species. Nearly all listing units faced high exposures to projected increases in stream temperature, sea surface temperature, and ocean acidification, but other aspects of exposure peaked in particular regions. Anthropogenic factors, especially migration barriers, habitat degradation, and hatchery influence, have reduced the adaptive capacity of most steelhead and salmon populations. Enhancing adaptive capacity is essential to mitigate for the increasing threat of climate change. Collectively, these results provide a framework to support recovery planning that considers climate impacts on the majority of West Coast anadromous salmonids.

Pathogens in space: Advancing understanding of pathogen dynamics and disease ecology through landscape genetics

Landscape genetics has provided many insights into how heterogeneous landscape features drive processes influencing spatial genetic variation in free-living organisms. This rapidly developing field has focused heavily on vertebrates, and expansion of this scope to the study of infectious diseases holds great potential for landscape geneticists and disease ecologists alike. The potential application of landscape genetics to infectious agents has garnered attention at formative stages in the development of landscape genetics, but systematic examination is lacking. We comprehensively review how landscape genetics is being used to better understand pathogen dynamics. We characterize the field and evaluate the types of questions addressed, approaches used and systems studied. We also review the now established landscape genetic methods and their realized and potential applications to disease ecology. Lastly, we identify emerging frontiers in the landscape genetic study of infectious agents, including recent phylogeographic approaches and frameworks for studying complex multihost and host-vector systems. Our review emphasizes the expanding utility of landscape genetic methods available for elucidating key pathogen dynamics (particularly transmission and spread) and also how landscape genetic studies of pathogens can provide insight into host population dynamics. Through this review, we convey how increasing awareness of the complementarity of landscape genetics and disease ecology among practitioners of each field promises to drive important cross-disciplinary advances.

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.

Associations of stream geomorphic conditions and prevalence of alternative reproductive tactics among sockeye salmon populations

In many species, males may exhibit alternative life histories to circumvent the costs of intrasexual competition and female courtship. While the evolution and underlying genetic and physiological mechanisms behind alternative reproductive tactics are well studied, there has been less consideration of the ecological factors that regulate their prevalence. Here, we examine six decades of age composition records from thirty-six populations of sockeye salmon (Oncorhynchus nerka) to quantify associations between spawning habitat characteristics and the prevalence of precocious sneakers known as 'jacks'. Jack prevalence was independent of neutral genetic structure among stream populations, but varied among habitat types and as a function of continuous geomorphic characteristics. Jacks were more common in streams relative to beaches and rivers, and their prevalence was negatively associated with stream width, depth, elevation, slope and area, but positively related to bank cover. Behavioural observations showed that jacks made greater use of banks, wood and shallows than guard males, indicating that their reproductive success depends on the availability of such refuges. Our results emphasize the role of the physical habitat in shaping reproductive tactic frequencies among populations, likely through local adaptation in response to variable fitness expectations under different geomorphic conditions.

Thermal adaptation and phenotypic plasticity in a warming world: Insights from common garden experiments on Alaskan sockeye salmon

An important unresolved question is how populations of coldwater-dependent fishes will respond to rapidly warming water temperatures. For example, the culturally and economically important group, Pacific salmon (Oncorhynchus spp.), experience site-specific thermal regimes during early development that could be disrupted by warming. To test for thermal local adaptation and heritable phenotypic plasticity in Pacific salmon embryos, we measured the developmental rate, survival, and body size at hatching in two populations of sockeye salmon (Oncorhynchus nerka) that overlap in timing of spawning but incubate in contrasting natural thermal regimes. Using a split half-sibling design, we exposed embryos of 10 families from each of two populations to variable and constant thermal regimes. These represented both experienced temperatures by each population, and predicted temperatures under plausible future conditions based on a warming scenario from the downscaled global climate model (MIROC A1B scenario). We did not find evidence of thermal local adaptation during the embryonic stage for developmental rate or survival. Within treatments, populations hatched within 1 day of each other, on average, and among treatments, did not differ in survival in response to temperature. We did detect plasticity to temperature; embryos developed 2.5 times longer (189 days) in the coolest regime compared to the warmest regime (74 days). We also detected variation in developmental rates among families within and among temperature regimes, indicating heritable plasticity. Families exhibited a strong positive relationship between thermal variability and phenotypic variability in developmental rate but body length and mass at hatching were largely insensitive to temperature. Overall, our results indicated a lack of thermal local adaptation, but a presence of plasticity in populations experiencing contrasting conditions, as well as family-specific heritable plasticity that could facilitate adaptive change.

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.

Genomic islands of divergence linked to ecotypic variation in sockeye salmon

Regions of the genome displaying elevated differentiation (genomic islands of divergence) are thought to play an important role in local adaptation, especially in populations experiencing high gene flow. However, the characteristics of these islands as well as the functional significance of genes located within them remain largely unknown. Here, we used data from thousands of SNPs aligned to a linkage map to investigate genomic islands of divergence in three ecotypes of sockeye salmon (Oncorhynchus nerka) from a single drainage in southwestern Alaska. We found ten islands displaying high differentiation among ecotypes. Conversely, neutral structure observed throughout the rest of the genome was low and not partitioned by ecotype. One island on linkage group So13 was particularly large and contained six SNPs with F_ST  > 0.14 (average F_ST of neutral SNPs = 0.01). Functional annotation revealed that the peak of this island contained a nonsynonymous mutation in a gene involved in growth in other species (TULP4). The islands that we discovered were relatively small (80-402 Kb), loci found in islands did not show reduced levels of diversity, and loci in islands displayed slightly elevated linkage disequilibrium. These attributes suggest that the islands discovered here were likely generated by divergence hitchhiking; however, we cannot rule out the possibility that other mechanisms may have produced them. Our results suggest that islands of divergence serve an important role in local adaptation with gene flow and represent a significant advance towards understanding the genetic basis of ecotypic differentiation.

Major histocompatibility complex diversity is positively associated with stream water temperatures in proximate populations of sockeye salmon

Local adaptation to heterogeneous environments generates population diversity within species, significantly increasing ecosystem stability and flows of ecosystem services. However, few studies have isolated the specific mechanisms that create and maintain this diversity. Here, we examined the relationship between water temperature in streams used for spawning and genetic diversity at a gene involved in immune function [the major histocompatibility complex (MHC)] in 14 populations of sockeye salmon (Oncorhynchus nerka) sampled across the Wood River basin in south-western Alaska. The largest influence on MHC diversity was lake basin, but we also found a significant positive correlation between average water temperature and MHC diversity. This positive relationship between temperature and MHC diversity appears to have been produced by natural selection at very local scales rather than neutral processes, as no correlation was observed between temperature and genetic diversity at 90 neutral markers. Additionally, no significant relationship was observed between temperature variability and MHC diversity. Although lake basin was the largest driver of differences in MHC diversity, our results also demonstrate that fine-scale differences in water temperature may generate variable selection regimes in populations that spawn in habitats separated by as little as 1 km. Additionally, our results indicated that some populations may be reaching a maximum level of MHC diversity. We postulate that salmon from populations in warm streams may delay spawning until late summer to avoid thermal stress as well as the elevated levels of pathogen prevalence and virulence associated with warm temperatures earlier in the summer.

The signature of fine scale local adaptation in Atlantic salmon revealed from common garden experiments in nature

Understanding the extent, scale and genetic basis of local adaptation (LA) is important for conservation and management. Its relevance in salmonids at microgeographic scales, where dispersal (and hence potential gene flow) can be substantial, has however been questioned. Here, we compare the fitness of communally reared offspring of local and foreign Atlantic salmon Salmo salar from adjacent Irish rivers and reciprocal F₁ hybrid crosses between them, in the wild ‘home’ environment of the local population. Experimental groups did not differ in wild smolt output but a catastrophic flood event may have limited our ability to detect freshwater performance differences, which were evident in a previous study. Foreign parr exhibited higher, and hybrids intermediate, emigration rates from the natal stream relative to local parr, consistent with genetically based behavioural differences. Adult return rates were lower for the foreign compared to the local group. Overall lifetime success of foreigners and hybrids relative to locals was estimated at 31% and 40% (mean of both hybrid groups), respectively. The results imply a genetic basis to fitness differences among populations separated by only 50 km, driven largely by variation in smolt to adult return rates. Hence even if supplementary stocking programs obtain broodstock from neighbouring rivers, the risk of extrinsic outbreeding depression may be high.