Parallel signatures of selection in temporally isolated lineages of pink salmon

Studying the effect of similar environments on diverse genetic backgrounds has long been a goal of evolutionary biologists with studies typically relying on experimental approaches. Pink salmon, a highly abundant and widely ranging salmonid, provide a naturally occurring opportunity to study the effects of similar environments on divergent genetic backgrounds due to a strict two‐year semelparous life history. The species is composed of two reproductively isolated lineages with overlapping ranges that share the same spawning and rearing environments in alternate years. We used restriction‐site‐associated DNA (RAD) sequencing to discover and genotype approximately 8000 SNP loci in three population pairs of even‐ and odd‐year pink salmon along a latitudinal gradient in North America. We found greater differentiation within the odd‐year than within the even‐year lineage and greater differentiation in the southern pair from Puget Sound than in the northern Alaskan population pairs. We identified 15 SNPs reflecting signatures of parallel selection using both a differentiation‐based method (BAYESCAN) and an environmental correlation method (BAYENV). These SNPs represent genomic regions that may be particularly informative in understanding adaptive evolution in pink salmon and exploring how differing genetic backgrounds within a species respond to selection from the same natural environment.     

Climate warming causes phenological shift in Pink Salmon, Oncorhynchus gorbuscha, behavior at Auke Creek, Alaska

Thirty‐four years (1972–2005) of water temperature data and extensive biological observations at Auke Creek, Alaska indicate a general warming trend that affected the native pink salmon (Oncorhynchus gorbuscha) population. Serial environmental records at nearby Auke Bay, Alaska over 46 years show trends of increasing air and sea surface temperatures. Trends of increased total precipitation and earlier date of ice out on nearby Auke Lake also occurred, but not at significant rates. Average water temperatures during the incubation of pink salmon in Auke Creek increased at a rate of 0.03 °C yr^?1 over the 34‐year period. For the 1972–2005 broods, midpoints of fry migrations from Auke Creek ranged between April 2 and May 7, and there was a trend of earlier migration of pink salmon fry at a rate of ? 0.5 days yr^?1. The migration timing of adult salmon into Auke Creek also showed a trend toward earlier timing. The earlier adult migration combined with warmer incubation temperatures are related to earlier migration of pink salmon fry. If the observed warming trend continues, Auke Creek may become unsuitable habitat for pink salmon. Given the trend for salmon fry to migrate earlier, a larger portion of the population may become mismatched with optimum environmental conditions during their early marine life history. If salmon adults continue to migrate into the creek earlier when water temperatures are commonly high, it will result in increased prespawning mortality.     

Anadromous red-spotted masu salmon (Oncorhynchus masou ishikawae), a southernmost sea-migration form of salmonid,displays low variation in both age at seaward migration and sea age

We examined variations in age at seaward migration and sea age for the anadromous form of red-spotted masu salmon (Oncorhynchus masou ishikawae) in two Japanese rivers. The anadromous form of red-spotted masu salmon expressed only two sea migration patterns in the two rivers: (a) the majority of the salmon (95%, n = 81) were of age-0, and age-1 migrants were rare (n = 4); and (b) all the salmon examined (n = 22) made a return migration within a year, with 23% of the salmon exhibiting potamodromy in the river. Owing to low variation in their sea migratory patterns, the anadromous form of red-spotted masu salmon is likely vulnerable to environmental fluctuations.     

Genetic monitoring of northern sea of Okhotsk populations of pink salmon (<Emphasis Type="Italic">Oncorhynchus gorbuscha</Emphasis>)

Results from the 1993–2004 genetic monitoring of pink salmon populations reproducing in the rivers of Tauy Bay on the Sea of Okhotsk are analyzed. A statistically significant heterogeneity of samples as determined by gene frequencies is found only in the pink salmon generations of even years. The genetic differentiation of samples from even years (G_ST = 1.39 ± 0.41) is higher than that of odd years (G_ST = 0.740.09). The pattern for the indicator of genetic variability (heterozygosity) is exactly the opposite (0.076 ± 0.02564 vs. 0.8760 ± 0.01950). Consequently, the lower-heterozygosity samples of lines from even years are on average more genetically distinct than the analogous indicator for odd years. In addition, the interpopulation ratio in the general value of genetic diversity is almost always smaller than both the intraannual and interannual ratios, leading to a low level of interpopulation genetic differences. Cluster analysis reveals that most 2001–2004 samples are grouped separately from samples collected prior to 2000. In our opinion, the reason for this could be the turnover of a numerically dominant generation of northern Sea of Okhtosk pink salmon and the change in gene frequencies accompanying it.     

Population Dynamics of the Pink Salmon Oncorhynchus gorbuschain the Basin of the Sea of Okhotsk in the 1990s

In the 1990s, an extensive body of data was gathered on the size of the Oncorhynchus gorbuschapink salmon populations of the Sea of Okhotsk at all the main developmental stages. A significant increase in numbers was found for juvenile pink salmon migrating into the offshore regions of the Sea of Okhotsk and the Pacific waters around the Kurils: from 250–450 million in 1990–1991 to 807–1570 million fish in 1993–1999. The overall number of migrating pink salmon in even years sharply increased in 1994 up to 215 million fish. After 1994, this estimate exceeded the number of migrating pink salmon in preceding odd years. Ocean survival of juvenile pink salmon gradually declined. This may be related to changes in the North Pacific pelagic ecosystems.     

Range-wide migration corridors and non-breeding areas of a northward expanding Afro-Palaearctic migrant,the European Bee-eater Merops apiaster

Across their ranges, different populations of migratory species often use separate routes to migrate between breeding and non-breeding grounds. Recent changes in climate and land-use have led to breeding range expansions in many species but it is unclear whether these populations also establish new migratory routes, non-breeding sites and migration phenology. Thus, we compared the migration patterns of European Bee-eaters Merops apiaster from two established western (n = 5) and eastern (n = 6) breeding populations in Europe, with those from a newly founded northern population (n = 19). We aimed to relate the breeding populations to the two known non-breeding clusters in Africa, and to test for similarities of migration routes and timing between the old and new populations. Western Bee-eaters used the western flyway to destinations in West Africa; the eastern birds uniformly headed south to southern African non-breeding sites, confirming a complete separation in time and space between these long-established populations. The recently founded northern population, however, also used a western corridor, but crossed the Mediterranean further east than the western population and overwintered mainly in a new non-breeding area in southern Congo/northern Angola. The migration routes and the new non-breeding range overlapped only slightly with the western, but not with the eastern, population. In contrast, migration phenology appeared to differ between the western and both the northern and the eastern populations, with tracked birds from the western population migrating 2–4 weeks earlier. The northern population thus shares some spatial traits with western Bee-eaters, but similar phenology only with eastern population. This divergence highlights the adjustments in the timing of migration to local environmental conditions in newly founded populations, and a parallel establishment of new breeding and non-breeding sites.     

Life-history divergence in Chinook salmon: historic contingency and parallel evolution

By jointly considering patterns of genetic and life-history diversity in over 100 populations of Chinook salmon from California to British Columbia, we demonstrate the importance of two different mechanisms for life-history evolution. Mapping adult run timing (the life-history trait most commonly used to characterize salmon populations) onto a tree based on the genetic data shows that the same run-time phenotypes exist in many different genetic lineages. In a hierarchical gene diversity analysis, differences among major geographic and ecological provinces explained the majority (62%) of the overall G(ST), whereas run-time differences explained only 10%. Collectively, these results indicate that run-timing diversity has developed independently by a process of parallel evolution in many different coastal areas. However, genetic differences between coastal populations with different run timing from the same basin are very modest (G(ST) < 0.02), indicating that evolutionary divergence of this trait linked to reproductive isolation has not led to parallel speciation, probably because of ongoing gene flow. A strikingly different pattern is seen in the interior Columbia River Basin, where run timing and other correlated life-history traits map cleanly onto two divergent genetic lineages (G(ST) approximately 0.15), indicating that some patterns of life-history diversity have a much older origin. Indeed, genetic data indicate that in the interior Columbia Basin, the two divergent lineages behave essentially as separate biological species, showing little evidence of genetic contact in spite of the fact that they comigrate through large areas of the river and ocean and in some locations spawn in nearly adjacent areas.     

Geographical and temporal mitochondrial DNA variability in populations of pink salmon

Pink salmon Oncorhynchus gorbuscha from odd and even year generations in rivers of Sakhalin Island, Kuril Island, Kamchatka Peninsula, and Alaska were investigated with five informative restriction endonucleases for mtDNA variation. The odd and even generations from the same rivers of South Sakhalin differed greatly. The time of divergence between the two broodlines was estimated at 0.9-1.1 Myr. The variability of mtDNA in odd year generations was higher than in even year generations and may have been due to' founder' and/or' bottleneck' effects. The differences among river populations within the Sakhalin region in 1991-1993 were not significant and this confirms the highly migratory nature of pink compared with other Pacific salmon. The mtDNA samples revealed statistically significant differences between regions. The northern populations (Kamchatka, Alaska) were less diverse in number and frequency of haplotypes than the southern populations (Sakhalin). This suggests that pink salmon originated in the Sakhalin-Kuril region and that a founder effect during the spread of this species may have restricted the mtDNA variability in other regions.     

Run timing of pelagic fishes in Gulf of St Lawrence: area and species effects

Run timings of four ubiquitous pelagic fishes in the Gulf of St Lawrence were compared for area and year effects. Run timing of alewife, Alosa pseudoharengus (Wilson), smelt, Osmerus mordax (Mitchill), and Atlantic salmon, Salmo salar L., was estimated at index rivers, and time of arrival of herring, Clupea harengus L., on spawning grounds was estimated from gillnet fisheries located around the Gulf, 1978–87. Time of pelagic fish migrations had strong species and area effects, but weak year effects. The species had their own unique migration times and durations, and could be divided into two groups: herring and smelt, which arrived in all areas of the Gulf at the same time; and alewife and salmon, which had significant area effects. Environmental data were not correlated with run timing. There was little annual variation in run timing among early-run species; only alewife had significant year effects. There was annual variation in the timing of late-run species, but it was not correlated with proximate environmental factors. Because run timing is an easily defined stock characteristic, observations over large areas and many years might be useful for understanding the impact of large-scale environmental changes on natural populations.     

Basin‐scale phenology and effects of climate variability on global timing of initial seaward migration of Atlantic salmon (Salmo salar)

Migrations between different habitats are key events in the lives of many organisms. Such movements involve annually recurring travel over long distances usually triggered by seasonal changes in the environment. Often, the migration is associated with travel to or from reproduction areas to regions of growth. Young anadromous Atlantic salmon (Salmo salar) emigrate from freshwater nursery areas during spring and early summer to feed and grow in the North Atlantic Ocean. The transition from the freshwater (‘parr’) stage to the migratory stage where they descend streams and enter salt water (‘smolt’) is characterized by morphological, physiological and behavioural changes where the timing of this parr‐smolt transition is cued by photoperiod and water temperature. Environmental conditions in the freshwater habitat control the downstream migration and contribute to within‐ and among‐river variation in migratory timing. Moreover, the timing of the freshwater emigration has likely evolved to meet environmental conditions in the ocean as these affect growth and survival of the post‐smolts. Using generalized additive mixed‐effects modelling, we analysed spatio‐temporal variations in the dates of downstream smolt migration in 67 rivers throughout the North Atlantic during the last five decades and found that migrations were earlier in populations in the east than the west. After accounting for this spatial effect, the initiation of the downstream migration among rivers was positively associated with freshwater temperatures, up to about 10 °C and levelling off at higher values, and with sea‐surface temperatures. Earlier migration occurred when river discharge levels were low but increasing. On average, the initiation of the smolt seaward migration has occurred 2.5 days earlier per decade throughout the basin of the North Atlantic. This shift in phenology matches changes in air, river, and ocean temperatures, suggesting that Atlantic salmon emigration is responding to the current global climate changes.     

Experimental microevolution: transplantation of pink salmon into the European North

Human-mediated translocations of species beyond their native ranges can enhance evolutionary processes in populations introduced to novel environments. We studied such processes in several generations of pink salmon Oncorhynchus gorbuscha introduced to the European North of Russia using a set of morphological and life-history traits as well as molecular genetic markers with different selective values: protein-coding loci, mtDNA, microsatellites, and MHC. The introduction of reproductively isolated pink salmon broodlines of odd and even years yielded different results. The odd-year broodline established self-reproducing local populations in many rivers of new range, but sustainable changes in external morphology, reproduction, and life-history, as well as the impoverishment of the gene pool occurred. Their successful colonisation of the new range resulted in specialisation manifested in the rapid directional shifts in some highly heritable phenotypic traits accompanied by increased homozygosity at molecular markers as a consequence of genetic drift and selective processes. The returns of transplanted pink salmon of even-year broodline decreased sharply already in the second generation, but there was no marked reduction of genetic diversity. Our data, as well as the analysis of the history of all pink salmon transplantations beyond the species range, demonstrate comparatively greater success of introduced odd-year broodline and permit to assume different adaptive plasticity of the even- and odd-year broodlines in pink salmon, what is most likely determined by differences in their evolutionary histories. Population genetic data suggest that the even-year broodline probably diverged from the odd-year broodline relatively recently and, due to the founder effect, may have lost a part of its genetic variation with which adaptive plasticity potential is associated.     

Comparative migration strategies of wild and captive‐bred Asian Houbara Chlamydotis macqueenii

For migratory species, the success of population reintroduction or reinforcement through captive‐bred released individuals depends on survivors undertaking appropriate migrations. We assess whether captive‐bred Asian Houbara Chlamydotis macqueenii from a breeding programme established with locally sourced individuals and released into suitable habitat during spring or summer undertake similar migrations to those of wild birds. Using satellite telemetry, we compare the migrations of 29 captive‐bred juveniles, 10 wild juveniles and 39 wild adults (including three birds first tracked as juveniles), examining migratory propensity (proportion migrating), timing, direction, stopover duration and frequency, efficiency (route deviation), and wintering and breeding season locations. Captive‐bred birds initiated autumn migration an average of 20.6 (±4.6 se) days later and wintered 470.8 km (±76.4) closer to the breeding grounds, mainly in Turkmenistan, northern Iran and Afghanistan, than wild birds, which migrated 1217.8 km (±76.4), predominantly wintering in southern Iran and Pakistan (juveniles and adults were similar). Wintering locations of four surviving captive‐bred birds were similar in subsequent years (median distance to first wintering site = 70.8 km, range 6.56–221.6 km), suggesting that individual captive‐bred birds (but not necessarily their progeny) remain faithful to their first wintering latitude. The migratory performance of captive‐bred birds was otherwise similar to that of wild juveniles. Although the long‐term fitness consequences for captive‐bred birds establishing wintering sites at the northern edge of those occupied by wild birds remain to be quantified, it is clear that the pattern of wild migrations established by long‐term selection is not replicated. If the shorter migration distance of young captive‐bred birds has a physiological rather than a genetic basis, then their progeny may still exhibit wild‐type migration. However, as there is a considerable genetic component to migration, captive breeding management must respect migratory population structure as well as natal and release‐site fidelity.     

Gene expression of neuronal soluble N-ethylmaleimide-sensitive factor attachment protein receptor complex components in the olfactory organ and brain during seaward and homeward migration by pink salmon (Oncorhynchus gorbuscha)

The expression of synaptic vesicle exocytosis-regulator SNARE complex component genes (snap25, stx1 and vamp2) was examined in the olfactory nervous system during seaward and homeward migration by pink salmon (Oncorhynchus gorbuscha). The expression levels of snares in the olfactory organ were higher in seaward fry than in feeding and homeward adults, reflecting the development of the olfactory nervous system. The expression of snap25a, b and stx1a was upregulated or stable in the adult olfactory bulb and telencephalon. This upregulated expression suggested alterations in olfactory neuronal plasticity that may be related to the discrimination of natal rivers. The expression of stx1b was downregulated in the adult olfactory bulb, but remained stable in the adult telencephalon. The expression of vamp2 was initially strong in seaward fry, but was downregulated in adults in both the olfactory bulb and telencephalon. Pink salmon has the lowest diversity of maturation age, the largest population, and the most evolutional position in Pacific salmon (genus Oncorhynchus). The expression of snares in the olfactory center of pink salmon reflected the timing of sexual maturation and homeward migration. The present results and our previous studies indicate that snares show distinct expression patterns between two salmon species that depend on physiological and ecological features of migration.     

Parallel adaptation to climate above the 35th parallel

Independent or parallel evolution of similar traits is key to understanding the genetics and limitations of adaptation. Adaptation from the same genetic changes in different populations defines parallel evolution. Such genetic changes can derive from standing ancestral variation or de novo mutations and excludes instances of adaptive introgression. In this issue of Molecular Ecology, Walden et al.(2020) investigate the scale of parallel climate adaptation from standing genetic variation between two North American Arabidopsis lyrata lineages, each formed by a distinct evolutionary history during the last glacial cycle. By identifying adaptive variants correlated with three ecologically significant climatic gradients, they show that instead of the same genetic variants or even genes, parallel evolution is only observed at the level of biological processes. The evolution of independent adaptive variants to climate in two genetically close lineages is explained by their different post‐glacial demographic histories. Separate glacial refugia and strong population bottlenecks were probably sufficient to change the landscape of shared allele frequencies, hindering the possibility of parallel evolution.     

Fine‐scale temporal adaptation within a salmonid population: mechanism and consequences

We demonstrate a clear example of local adaptation of seasonal timing of spawning and embryo development. The consequence is a population of pink salmon that is segmented into spawning groups that use the same limited habitat. We synthesize published observations with results of new analyses to demonstrate that genetic variation of these traits results in survival differentials related to that variation, and that density‐dependent embryo mortality and seasonally variable juvenile mortality are a mechanism of selection. Most examples of local adaptation in natural systems depend on observed correlations between environments and fitness traits, but do not fully demonstrate local adaptation: that the trait is genetically determined, exhibits different fitness in common environments or across different environments, and its variation is mechanistically connected to fitness differences. The geographic or temporal scales of local adaptation often remain obscure. Here, we show that heritable, fine‐scale differences of timing of reproductive migration in a pink salmon (Oncorhynchus gorbuscha) resulted in temporal structure that persisted several generations; the differences enable a density‐dependent population to pack more spawners into limited spawning habitat, that is, enhance its fitness. A balanced trade‐off of survivals results because embryos from early‐migrating fish have a lower freshwater survival (harsh early physical conditions and disturbance by late spawners), but emigrant fry from late‐migrating fish have lower marine survivals (timing of their vernal emergence into the estuarine environment). Such fine‐scale local adaptations increase the genetic portfolio of the populations and may provide a buffer against the impacts of climate change.     

Integration of Random Forest with population‐based outlier analyses provides insight on the genomic basis and evolution of run timing in Chinook salmon (Oncorhynchus tshawytscha)

Anadromous Chinook salmon populations vary in the period of river entry at the initiation of adult freshwater migration, facilitating optimal arrival at natal spawning. Run timing is a polygenic trait that shows evidence of rapid parallel evolution in some lineages, signifying a key role for this phenotype in the ecological divergence between populations. Studying the genetic basis of local adaptation in quantitative traits is often impractical in wild populations. Therefore, we used a novel approach, Random Forest, to detect markers linked to run timing across 14 populations from contrasting environments in the Columbia River and Puget Sound, USA. The approach permits detection of loci of small effect on the phenotype. Divergence between populations at these loci was then examined using both principle component analysis and F_ST outlier analyses, to determine whether shared genetic changes resulted in similar phenotypes across different lineages. Sequencing of 9107 RAD markers in 414 individuals identified 33 predictor loci explaining 79.2% of trait variance. Discriminant analysis of principal components of the predictors revealed both shared and unique evolutionary pathways in the trait across different lineages, characterized by minor allele frequency changes. However, genome mapping of predictor loci also identified positional overlap with two genomic outlier regions, consistent with selection on loci of large effect. Therefore, the results suggest selective sweeps on few loci and minor changes in loci that were detected by this study. Use of a polygenic framework has provided initial insight into how divergence in a trait has occurred in the wild.     

On possible re-distribution of pink salmon <Emphasis Type="Italic">Onchorhynchus gorbuscha</Emphasis> between the reproduction areas of different stocks in the Sakhalin-Kuril region

The scale structure of pink salmon Oncorhynchus gorbuscha inhabiting southern Sakhalin and Iturup Island were studied in 2014. These locations are characterized by the unusual dynamics of catches and biological parameters of the fish that come to spawn in different areas in the Sakhalin-Kuril region. The pink salmon that has originated from Iturup Island appeared en masse in the waters of southern Sakhalin, according to the analysis of the scleritogram fragments that reflect the fish growth during the first months of life. This is the first confirmation of the hypothesis of the fluctuating stocks of pink salmon approved by the ichthyological methods. It is suggested that wide-scale straying of pink salmon takes place in the years of dominance shift between the generations of even and odd spawning years.     

Distribution of genetic variation underlying adult migration timing in steelhead of the Columbia River basin

Fish migrations are energetically costly, especially when moving between freshwater and saltwater, but are a viable strategy for Pacific salmon and trout (Oncorhynchus spp.) due to the advantageous resources available at various life stages. Anadromous steelhead (O. mykiss) migrate vast distances and exhibit variation for adult migration phenotypes that have a genetic basis at candidate genes known as greb1L and rock1. We examined the distribution of genetic variation at 13 candidate markers spanning greb1L, intergenic, and rock1 regions versus 226 neutral markers for 113 populations (n = 9,471) of steelhead from inland and coastal lineages in the Columbia River. Patterns of population structure with neutral markers reflected genetic similarity by geographic region as demonstrated in previous studies, but candidate markers clustered populations by genetic variation associated with adult migration timing. Mature alleles for late migration had the highest frequency overall in steelhead populations throughout the Columbia River, with only 9 of 113 populations that had a higher frequency of premature alleles for early migration. While a single haplotype block was evident for the coastal lineage, we identified multiple haplotype blocks for the inland lineage. The inland lineage had one haplotype block that corresponded to candidate markers within the greb1L gene and immediately upstream in the intergenic region, and the second block only contained candidate markers from the intergenic region. Haplotype frequencies had similar patterns of geographic distribution as single markers, but there were distinct differences in frequency between the two haplotype blocks for the inland lineage. This may represent multiple recombination events that differed between lineages where phenotypic differences exist between freshwater entry versus arrival timing as indicated by Micheletti et al. (2018a). Redundancy analyses were used to model environmental effects on allelic frequencies of candidate markers, and significant variables were migration distance, temperature, isothermality, and annual precipitation. This study improves our understanding of the spatial distribution of genetic variation underlying adult migration timing in steelhead as well as associated environmental factors and has direct conservation and management implications.     

Environmental adaptation in Chinook salmon (Oncorhynchus tshawytscha) throughout their North American range

Landscape genomics is a rapidly growing field with recent advances in both genotyping efficiency and statistical analyses that provide insight towards local adaptation of populations under varying environmental and selective pressure. Chinook salmon (Oncorhynchus tshawytscha) are a broadly distributed Pacific salmon species, occupying a diversity of habitats throughout the northeastern Pacific with pronounced variation in environmental and climate features but little is understood regarding local adaptation in this species. We used a multivariate method, redundancy analysis (RDA), to identify polygenic correlations between 19 703 SNP loci and a suite of environmental variables in 46 collections of Chinook salmon (1956 total individuals) distributed throughout much of its North American range. Models in RDA were conducted on both rangewide and regional scales by hierarchical partitioning of the populations into three distinct genetic lineages. Our results indicate that between 5.8 and 21.8% of genomic variation can be accounted for by environmental features, and 566 putatively adaptive loci were identified as targets of environmental adaptation. The most influential drivers of adaptive divergence included precipitation in the driest quarter of the year (Rangewide and North Coastal Lineage, anova P = 0.002 and 0.01, respectively), precipitation in the wettest quarter of the year (Interior Columbia River Stream‐Type Lineage, anova P = 0.03), variation in mean diurnal range in temperature (South Coastal Lineage, anova P = 0.005), and migration distance (Rangewide, anova P = 0.001). Our results indicate that environmental features are strong drivers of adaptive genomic divergence in this species, and provide a foundation to investigate how Chinook salmon might respond to global environmental change.