Founding events influence genetic population structure of sockeye salmon (Oncorhynchus nerka) in Lake Clark, Alaska

Bottlenecks can have lasting effects on genetic population structure that obscure patterns of contemporary gene flow and drift. Sockeye salmon are vulnerable to bottleneck effects because they are a highly structured species with excellent colonizing abilities and often occupy geologically young habitats. We describe genetic divergence among and genetic variation within spawning populations of sockeye salmon throughout the Lake Clark area of Alaska. Fin tissue was collected from sockeye salmon representing 15 spawning populations of Lake Clark, Six-mile Lake, and Lake Iliamna. Allele frequencies differed significantly at 11 microsatellite loci in 96 of 105 pairwise population comparisons. Pairwise estimates of FST ranged from zero to 0.089. Six-mile Lake and Lake Clark populations have historically been grouped together for management purposes and are geographically proximate. However, Six-mile Lake populations are genetically similar to Lake Iliamna populations and are divergent from Lake Clark populations. The reduced allelic diversity and strong divergence of Lake Clark populations relative to Six-mile Lake and Lake Iliamna populations suggest a bottleneck associated with the colonization of Lake Clark by sockeye salmon. Geographic distance and spawning habitat differences apparently do not contribute to isolation and divergence among populations. However, temporal isolation based on spawning time and founder effects associated with ongoing glacial retreat and colonization of new spawning habitats contribute to the genetic population structure of Lake Clark sockeye salmon. Nonequilibrium conditions and the strong influence of genetic drift caution against using estimates of divergence to estimate gene flow among populations of Lake Clark sockeye salmon.     

Lack of parallel genetic patterns underlying the repeated ecological divergence of beach and stream‐spawning kokanee salmon

Recent progress in methods for detecting adaptive population divergence in situ shows promise for elucidating the conditions under which selection acts to generate intraspecific diversity. Rapid ecological diversification is common in fishes; however, the role of phenotypic plasticity and adaptation to local environments is poorly understood. It is now possible to investigate genetic patterns to make inferences regarding phenotypic traits under selection and possible mechanisms underlying ecotype divergence, particularly where similar novel phenotypes have arisen in multiple independent populations. Here, we employed a bottom‐up approach to test for signatures of directional selection associated with divergence of beach‐ and stream‐spawning kokanee, the obligate freshwater form of sockeye salmon (Oncorhynchus nerka). Beach‐ and stream‐spawners co‐exist in many post‐glacial lakes and exhibit distinct reproductive behaviours, life‐history traits and spawning habitat preferences. Replicate ecotype pairs across five lakes in British Columbia, Canada were genotyped at 57 expressed sequence tag‐linked and anonymous microsatellite loci identified in a previous genome scan. Fifteen loci exhibited signatures of directional selection (high F_ST outliers), four of which were identified in multiple lakes. However, the lack of parallel genetic patterns across all lakes may be a result of: 1) an inability to detect loci truly under selection; 2) alternative genetic pathways underlying ecotype divergence in this system; and/or 3) phenotypic plasticity playing a formative role in driving kokanee spawning habitat differences. Gene annotations for detected outliers suggest pathogen resistance and energy metabolism as potential mechanisms contributing to the divergence of beach‐ and stream‐spawning kokanee, but further study is required.     

Contrasting patterns of morphological and neutral genetic divergence among geographically proximate populations of sockeye salmon Oncorhynchus nerka in Lake Aleknagik,Alaska

Levels of differentiation in morphological traits (age at maturity, body length at age, egg mass and body depth) and spawning time were examined in sockeye salmon Oncorhynchus nerka from three geographically proximate but physically distinct creeks in Lake Aleknagik, Alaska. Happy Creek fish had significantly greater values for most measured morphological traits, and Eagle Creek fish spawned significantly later than fish in the other creeks. Phenotypic differentiation between creeks, measured as P_ST, was then compared with microsatellite marker differentiation between creeks, measured as F_ST. No correlations were apparent between P_ST and F_ST values, and P_ST values were generally significantly larger than zero (P_ST= 0·0018–0·31) whereas F_ST values were not (F_ST=?0·0004 to 0·0016). The insignificant pair‐wise F_ST values between creek samples indicated that gene flow occurs between creeks, assuming the creek populations have reached migration–drift equilibrium. However, the strong homing behaviour of sockeye salmon precludes a scenario in which fish from the three creeks constitute a single population that segregates by body size. Rather, significant phenotypic differentiation suggests that strong divergent selection occurs on the phenotypic traits despite the homogenizing effects of gene flow.     

MOLECULAR GENETIC EVIDENCE FOR PARALLEL LIFE-HISTORY EVOLUTION WITHIN A PACIFIC SALMON (SOCKEYE SALMON AND KOKANEE,ONCORHYNCHUS NERKA)

The Pacific salmon Oncorhynchus nerka typically occurs as a sea-run form (sockeye salmon) or may reside permanently in lakes (kokanee) thoughout its native North Pacific. We tested whether such geographically extensive ecotypic variation resulted from parallel evolutionary divergence thoughout the North Pacific or whether the two forms are monophyletic groups by examining allelic variation between sockeye salmon and kokanee at two minisatellite DNA repeat loci and in mitochondrial DNA (mtDNA) Bgl II restriction sites. Our examination of over 750 fish from 24 populations, ranging from Kamchatka to the Columbia River, identified two major genetic groups of North Pacific O. nerka: a “northwestern” group consisting of fish from Kamchatka, western Alaska, and northwestern British Columbia, and a “southern” group consisting of sockeye salmon and kokanee populations from the Fraser and Columbia River systems. Maximum-likelihood analysis accompanied by bootstrapping provided strong support for these two genetic groups of O. nerka; the populations did not cluster by migratory form, but genetic affinities were organized more strongly by geographic proximity. The two major genetic groups resolved in our study probably stem from historical isolation and dispersal of O. nerka from two major Wisconsinan glacial refugia in the North Pacific. There were significant minisatellite DNA allele frequency differences between sockeye salmon and kokanee populations from different parts of the same watershed, between populations spawning in different tributaries of the same lake, and also between sympatric populations spawning in the same stream at the same time. MtDNA Bgl II restriction site variation was significant between sockeye salmon and kokanee spawning in different parts of the same major watershed but not between forms spawning in closer degrees of reproductive sympatry. Patterns of genetic affinity and allele sharing suggested that kokanee have arisen from sea-run sockeye salmon several times independently in the North Pacific. We conclude that sockeye salmon and kokanee are para- and polyphyletic, respectively, and that the present geographic distribution of the ecotypes results from parallel evolutionary origins of kokanee from sockeye (divergences between them) thoughout the North Pacific.     

Potential adaptive divergence between subspecies and populations of snapdragon plants inferred from QST–FST comparisons

Phenotypic divergence among natural populations can be explained by natural selection or by neutral processes such as drift. Many examples in the literature compare putatively neutral (F_ST) and quantitative genetic (Q_ST) differentiation in multiple populations to assess their evolutionary signature and identify candidate traits involved with local adaptation. Investigating these signatures in closely related or recently diversified species has the potential to shed light on the divergence processes acting at the interspecific level. Here, we conducted this comparison in two subspecies of snapdragon plants (eight populations of Antirrhinum majus pseudomajus and five populations of A. m. striatum) in a common garden experiment. We also tested whether altitude was involved with population phenotypic divergence. Our results identified candidate phenological and morphological traits involved with local adaptation. Most of these traits were identified in one subspecies but not the other. Phenotypic divergence increased with altitude for a few biomass‐related traits, but only in A. m. striatum. These traits therefore potentially reflect A. m. striatum adaptation to altitude. Our findings imply that adaptive processes potentially differ at the scale of A. majus subspecies.     

Resurrecting an extinct salmon evolutionarily significant unit: archived scales, historical DNA and implications for restoration

Archival scales from 603 sockeye salmon (Oncorhynchus nerka), sampled from May to July 1924 in the lower Columbia River, were analysed for genetic variability at 12 microsatellite loci and compared to 17 present‐day O. nerka populations—exhibiting either anadromous (sockeye salmon) or nonanadromous (kokanee) life histories—from throughout the Columbia River Basin, including areas upstream of impassable dams built subsequent to 1924. Statistical analyses identified four major genetic assemblages of sockeye salmon in the 1924 samples. Two of these putative historical groupings were found to be genetically similar to extant evolutionarily significant units (ESUs) in the Okanogan and Wenatchee Rivers (pairwise F_ST = 0.004 and 0.002, respectively), and assignment tests were able to allocate 77% of the fish in these two historical groupings to the contemporary Okanogan River and Lake Wenatchee ESUs. A third historical genetic grouping was most closely aligned with contemporary sockeye salmon in Redfish Lake, Idaho, although the association was less robust (pairwise F_ST = 0.060). However, a fourth genetic grouping did not appear to be related to any contemporary sockeye salmon or kokanee population, assigned poorly to the O. nerka baseline, and had distinctive early return migration timing, suggesting that this group represents a historical ESU originating in headwater lakes in British Columbia that was probably extirpated sometime after 1924. The lack of a contemporary O. nerka population possessing the genetic legacy of this extinct ESU indicates that efforts to reestablish early‐migrating sockeye salmon to the headwater lakes region of the Columbia River will be difficult.     

Gene flow does not prevent personality and morphological differentiation between two blue tit populations

Understanding the causes and consequences of population phenotypic divergence is a central goal in ecology and evolution. Phenotypic divergence among populations can result from genetic divergence, phenotypic plasticity or a combination of the two. However, few studies have deciphered these mechanisms for populations geographically close and connected by gene flow, especially in the case of personality traits. In this study, we used a common garden experiment to explore the genetic basis of the phenotypic divergence observed between two blue tit (Cyanistes caeruleus) populations inhabiting contrasting habitats separated by 25 km, for two personality traits (exploration speed and handling aggression), one physiological trait (heart rate during restraint) and two morphological traits (tarsus length and body mass). Blue tit nestlings were removed from their population and raised in a common garden for up to 5 years. We then compared adult phenotypes between the two populations, as well as trait‐specific Q_st and F_st. Our results revealed differences between populations similar to those found in the wild, suggesting a genetic divergence for all traits. Q_st–F_st comparisons revealed that the trait divergences likely result from dissimilar selection patterns rather than from genetic drift. Our study is one of the first to report a Q_st–F_st comparison for personality traits and adds to the growing body of evidence that population genetic divergence is possible at a small scale for a variety of traits including behavioural traits.     

Size differentiation in Finnish house sparrows follows Bergmann's rule with evidence of local adaptation

Bergmann's rule predicts that individuals are larger in more poleward populations and that this size gradient has an adaptive basis. Hence, phenotypic divergence in size traits between populations (P_ST) is expected to exceed the level of divergence by drift alone (F_ST). We measured 16 skeletal traits, body mass and wing length in 409 male and 296 female house sparrows Passer domesticus sampled in 12 populations throughout Finland, where the species has its northernmost European distributional margin. Morphometric differentiation across populations (P_ST) was compared with differentiation in 13 microsatellites (F_ST). We find that twelve traits phenotypically diverged more than F_ST in both sexes, and an additional two traits diverged in males. The phenotypic divergence exceeded F_ST in several traits to such a degree that findings were robust also to strong between‐population environmental effects. Divergence was particularly strong in dimensions of the bill, making it a strong candidate for the study of adaptive molecular genetic divergence. Divergent traits increased in size in more northern populations. We conclude that house sparrows show evidence of an adaptive latitudinal size gradient consistent with Bergmann's rule on the modest spatial scale of ca. 600 km.     

Morphometric differentiation across House Sparrow Passer domesticus populations in Finland in comparison with the neutral expectation for divergence

To understand the biology of organisms it is important to take into account the evolutionary forces that have acted on their constituent populations. Neutral genetic variation is often assumed to reflect variation in quantitative traits under selection, though with even low neutral divergence there can be substantial differentiation in quantitative genetic variation associated with locally adapted phenotypes. To study the relative roles of natural selection and genetic drift in shaping phenotypic variation, the levels of quantitative divergence based on phenotypes (P_ST) and neutral genetic divergence (F_ST) can be compared. Such a comparison was made between 10 populations of Finnish House Sparrows (n = 238 individuals) collected in 2009 across the whole country. Phenotypic variation in tarsus‐length, wing‐length, bill‐depth, bill‐length and body mass were considered and 13 polymorphic microsatellite loci were analysed to quantify neutral genetic variation. Calculations of P_ST were based on Markov‐Chain Monte Carlo Bayesian estimates of phenotypic variances across and within populations. The robustness of the conclusions of the P_ST–F_ST comparison was evaluated by varying the proportion of variation due to additive genetic effects within and across populations. Our results suggest that body mass is under directional selection, whereas the divergence in other traits does not differ from neutral expectations. These findings suggest candidate traits for considering gene‐based studies of local adaptation. The recognition of locally adapted populations may be of value in the conservation of this declining species.     

Low levels of quantitative and molecular genetic differentiation among natural populations of Medicago ciliaris Kroch. (Fabaceae) of different Tunisian eco-geographical origin

In this paper, we analyze the genetic variability in four Tunisian natural populations of Medicago ciliaris using 19 quantitative traits and six polymorphic microsatellite loci. We investigated the amplification transferability of 30 microsatellites developed in the model legume M. truncatula to M. ciliaris. Results revealed that about 56.66% of analyzed markers are valuable genetic markers for M. ciliaris. The most genetic diversity at quantitative traits and microsatellite loci was found to occur within populations (>80%). Low differentiations among populations at quantitative traits Q _ST  = 0.146 and molecular markers F _ST  = 0.18 were found. The majority of measured traits exhibited no significant difference in the level of Q _ST and F _ST . Furthermore, significant correlations established between these traits and eco-geographical factors suggested that natural selection should be invoked to explain the level of phenotypic divergence among populations rather than drift. There was no significant correlation between population differentiation at quantitative traits and molecular markers. Significant spatial genetic structure consistent with models of isolation by distance was detected within all studied populations. The site-of-origin environmental factors explain about 9.07% of total phenotypic genetic variation among populations. The eco-geographical factors that influence more the variation of measured traits among populations are the soil texture and altitude. Nevertheless, there were no consistent pattern of associations between gene diversity (He) and environmental factors.     

Post‐glacial colonization routes coincide with a life‐history breakpoint along a latitudinal gradient

Although adaptive divergence along environmental gradients has repeatedly been demonstrated, the role of post‐glacial colonization routes in determining phenotypic variation along gradients has received little attention. Here, we used a hierarchical Q_ST–F_ST approach to separate the roles of adaptive and neutral processes in shaping phenotypic variation in moor frog (Rana arvalis) larval life histories along a 1,700 km latitudinal gradient across northern Europe. This species has colonized Scandinavia via two routes with a contact zone in northern Sweden. By using neutral SNP and common garden phenotypic data from 13 populations at two temperatures, we showed that most of the variation along the gradient occurred between the two colonizing lineages. We found little phenotypic divergence within the lineages; however, all phenotypic traits were strongly diverged between the southern and northern colonization routes, with higher growth and development rates and larger body size in the north. The Q_ST estimates between the colonization routes were four times higher than F_ST, indicating a prominent role for natural selection. Q_ST within the colonization routes did not generally differ from F_ST, but we found temperature‐dependent adaptive divergence close to the contact zone. These results indicate that lineage‐specific variation can account for much of the adaptive divergence along a latitudinal gradient.     

GENETIC AND PHENOTYPIC DIVERGENCE BETWEEN LOW‐ AND HIGH‐ALTITUDE POPULATIONS OF TWO RECENTLY DIVERGED CINNAMON TEAL SUBSPECIES

Spatial variation in the environment can lead to divergent selection between populations occupying different parts of a species’ range, and ultimately lead to population divergence. The colonization of new areas can thus facilitate divergence in beneficial traits, yet with little differentiation at neutral genetic markers. We investigated genetic and phenotypic patterns of divergence between low‐ and high‐altitude populations of cinnamon teal inhabiting normoxic and hypoxic regions in the Andes and adjacent lowlands of South America. Cinnamon teal showed strong divergence in body size (PC1; P_ST= 0.56) and exhibited significant frequency differences in a single nonsynonymous α‐hemoglobin amino acid polymorphism (Asn/Ser‐α9; F_ST= 0.60) between environmental extremes, despite considerable admixture of mtDNA and intron loci (F_ST= 0.004–0.168). Inferences of strong population segregation were further supported by the observation of few mismatched individuals in either environmental extreme. Coalescent analyses indicated that the highlands were most likely colonized from lowland regions but following divergence, gene flow has been asymmetric from the highlands into the lowlands. Multiple selection pressures associated with high‐altitude habitats, including cold and hypoxia, have likely shaped morphological and genetic divergence within South American cinnamon teal populations.     

Divergence at neutral and non-neutral loci in Drosophila buzzatii populations and their hybrids

The impact of intraspecific hybridisation on fitness and morphological traits depends on the history of natural selection and genetic drift, which may have led to differently coadapted gene-complexes in the parental populations. The divergence at neutral and non-neutral loci between populations can be evaluated by estimating F_ST and Q_ST respectively, and hence give an estimate of drift and selection in the populations. Here we investigate (1) whether divergence between populations in quantitative traits (wing size and shape) can be attributed to selection or drift alone, (2) The impact of intraspecific hybridisation on estimators for divergence at neutral (F_ST) and non-neutral loci (Q_ST) in hybrids, (3) If measurement of shape is more informative than size in order to detect divergence in quantitative traits between populations. The aims were addressed by performing two hybridisations between three populations of Drosophila buzzatii, one between populations from Argentina and the Canary Islands (separated for 200 years), and the other between populations from Argentina and Australia (separated for 80 years). We observed the highest divergence at neutral loci between the Argentinean and Canary Island populations, but highest morphological divergence between the Argentinean and Australian populations, indicating that natural selection is acting on the wings. Divergence based on Q_ST measures in the hybrids was sensitive towards increased phenotypic variance (σ²p) within groups and should be used with care when σ²p of populations differ. Our results indicate that measures of shape give a better estimate of divergence at the underlying quantitative traits loci than measures of size.     

Genetic and Morphological Divergence in Three Strains of Brook Trout Salvelinus fontinalis Commonly Stocked in Lake Superior

Fitness related traits often show spatial variation across populations of widely distributed species. Comparisons of genetic variation among populations in putatively neutral DNA markers and in phenotypic traits susceptible to selection (Q_ST F_ST analysis) can be used to determine to what degree differentiation among populations can be attributed to selection or genetic drift. Traditionally, Q_ST F_ST analyses require a large number of populations to achieve sufficient statistical power; however, new methods have been developed that allow Q_ST F_ST comparisons to be conducted on as few as two populations if their pedigrees are informative. This study compared genetic and morphological divergence in three strains of brook trout Salvelinus fontinalis that were historically or currently used for stocking in the Lake Superior Basin. Herein we examined if morphological divergence among populations showed temporal variation, and if divergence could be attributed to selection or was indistinguishable from genetic drift. Multivariate Q_ST F_ST analysis showed evidence for divergent selection between populations. Univariate analyses suggests that the pattern observed in the multivariate analyses was largely driven by divergent selection for length and weight, and moreover by divergence between the Assinica strain and each of the Iron River and Siskiwit strains rather than divergent selection between each population pair. While it could not be determined if divergence was due to natural selection or inadvertent artificial selection in hatcheries, selected differences were consistent with patterns of domestication commonly found in salmonids.     

Population genetic structure and a possible role for selection in driving phenotypic divergence in a rainbowfish (Melanotaeniidae)

Quantitative genetic divergence may be driven by drift or selection. The rainbowfish Melanotaenia australis exhibits phenotypic divergence among populations in Western Australia, although the mechanisms driving this divergence are unknown. We used microsatellites to assess neutral genetic divergence (F_ST), and found a hierarchical pattern of subdivision consistent with low divergence between upstream and downstream populations (within drainages), moderate divergence between drainages (within regions), and high divergence between regions. Using a common‐garden approach, we measured quantitative genetic divergence in phenotypic traits (Q_ST). By comparing this to expectations from neutral processes (F_ST), we concluded that the effect of selection varies depending on the spatial scale considered. Within drainages, selection may be causing divergence between upstream and downstream phenotypes but, between regions, selection appears to homogenize phenotypes. This highlights the importance of spatial scale in studies of this kind, and suggests that, because variance in selection pressures can drive speciation, polymorphism in M. australis may represent speciation in action. © 2010 The Linnean Society of London, Biological Journal of the Linnean Society, 2011, 102 , 144–160.     

Phenotypic divergence of the common toad (Bufo bufo) along an altitudinal gradient: evidence for local adaptation

Variation in the environment can induce different patterns of genetic and phenotypic differentiation among populations. Both neutral processes and selection can influence phenotypic differentiation. Altitudinal phenotypic variation is of particular interest in disentangling the interplay between neutral processes and selection in the dynamics of local adaptation processes but remains little explored. We conducted a common garden experiment to study the phenotypic divergence in larval life-history traits among nine populations of the common toad (Bufo bufo) along an altitudinal gradient in France. We further used correlation among population pairwise estimates of quantitative trait (Q_ST) and neutral genetic divergence (F_ST from neutral microsatellite markers), as well as altitudinal difference, to estimate the relative role of divergent selection and neutral genetic processes in phenotypic divergence. We provided evidence for a neutral genetic differentiation resulting from both isolation by distance and difference in altitude. We found evidence for phenotypic divergence along the altitudinal gradient (faster development, lower growth rate and smaller metamorphic size). The correlation between pairwise Q_STs–F_STs and altitude differences suggested that this phenotypic differentiation was most likely driven by altitude-mediated selection rather than by neutral genetic processes. Moreover, we found different divergence patterns for larval traits, suggesting that different selective agents may act on these traits and/or selection on one trait may constrain the evolution on another through genetic correlation. Our study highlighted the need to design more integrative studies on the common toad to unravel the underlying processes of phenotypic divergence and its selective agents in the context of environmental clines.     

Divergence of mitochondrial DNA in a population of sockeye salmon (Oncorhynchus nerka Walbaum) from Azabache'e lake (Kamchatka)

Variability of three PCR-amplified mtDNA regions was examined in five populations of sockeye salmon from Azabach'e Lake. Eighteen haplotypes were detected in 144 fish. Significant differences were found between seasonal races of sockeye salmon spawning in the lake. The short time of independent divergence between the seasonal races indicates that these races formed independently in each spawning region. No difference in mtDNA between lake samples of early sockeye salmon (subisolates) was revealed, which confirms the existence of gene flow between them. A high level of differences between the sockeye salmon spawning in the lake and spawning in the tributaries of the lake, the Bushuev and Lotnaya rivers, suggests that there were no migration between them during many generations and that the nature of spawning grounds (lake or river) is essential for within species differentiation in this species.     

The role of selection and historical factors in driving population differentiation along an elevational gradient in an island bird

Adaptation to local environmental conditions and the range dynamics of populations can influence evolutionary divergence along environmental gradients. Thus, it is important to investigate patterns of both phenotypic and genetic variations among populations to reveal the respective roles of these two types of factors in driving population differentiation. Here, we test for evidence of phenotypic and genetic structure across populations of a passerine bird (Zosterops borbonicus) distributed along a steep elevational gradient on the island of Réunion. Using 11 microsatellite loci screened in 401 individuals from 18 localities distributed along the gradient, we found that genetic differentiation occurred at two spatial levels: (i) between two main population groups corresponding to highland and lowland areas, respectively, and (ii) within each of these two groups. In contrast, several morphological traits varied gradually along the gradient. Comparison of neutral genetic differentiation (F_ST) and phenotypic differentiation (P_ST) showed that P_ST largely exceeds F_ST at several morphological traits, which is consistent with a role for local adaptation in driving morphological divergence along the gradient. Overall, our results revealed an area of secondary contact midway up the gradient between two major, cryptic, population groups likely diverged in allopatry. Remarkably, local adaptation has shaped phenotypic differentiation irrespective of population history, resulting in different patterns of variation along the elevational gradient. Our findings underscore the importance of understanding both historical and selective factors when trying to explain variation along environmental gradients.     

The effect of the infectious hematopoietic necrosis virus on the population of sockeye salmon <Emphasis Type="Italic">Oncorhynchus nerka</Emphasis> (Salmoniformes,Salmonidae) from Lake Nachikinskoe

The results of virological examination of adult sockeye salmon Oncorhynchus nerka in Lake Nachikinskoe from 2003 to 2008 demonstrated high values of prevalence of the infectious hematopoietic necrosis virus (IHNV) in the fish population (up to 100%). In the studied water reservoir there are optimal conditions for the development and revealing of the epizootic of infectious hematopoietic necrosis (IHN). In Lake Nachikinskoe, two natural epizootics of IHN in juvenile sockeye salmon were recorded in 2003 and 2006. The tendencies for change of prevalence of IHNV by years were analyzed, and their relation to the numbers of adult sockeye salmon at the spawning grounds of Lake Nachikinskoe was determined. The correlation of IHNV prevalence with the numbers of adult sockeye salmon in the studied water reservoir is direct but insignificant (r _s = 0.695, p < 0.05). It was shown that the epizootic of 2003 could affect the fact that fish return in 2007 and 2008 was minimum for the considered time interval while the numbers of their parents were, on the other hand, very high.