The impact of fishing-induced mortality on the evolution of alternative life-history tactics in brook charr.

Although contemporary trends indicative of evolutionary change have been detected in the life-history traits of exploited populations, it is not known to what extent fishing influences the evolution of alternative life-history tactics in migratory species such as salmonids. Here, we build a model to predict the evolution of anadromy and residency in an exploited population of brook charr, Salvelinus fontinalis. Our model allows for both phenotypic plasticity and genetic change in the age and size at migration by including migration reaction norms. Using this model, we predict that fishing of anadromous individuals over the course of 100 years causes evolution in the migration reaction norm, resulting in a decrease in average probabilities of migration with increasing harvest rate. Moreover, we show that differences in natural mortalities in freshwater greatly influence the magnitude and rate of evolutionary change. The fishing-induced changes in migration predicted by our model alter population abundances and reproductive output and should be accounted for in the sustainable management of salmonids.

Introduction

Fishing is now acknowledged as a potential evolutionary force, described as a ‘massive uncontrolled experiment in evolutionary selection’ (Stokes and Law 2000). Whenever individuals with certain characteristics are more likely to survive harvest or to produce more viable offspring than others, fishing can induce evolutionary changes in life-history traits (Law and Grey 1989; Jørgensen 1990; Sutherland 1990; Law 1991, 2000; Smith 1994; Haugen and Vøllestad 2001; Conover and Munch 2002; Barot et al. 2004; Olsen et al. 2004; Reznick and Ghalambor 2005). That fishing can generate substantial selection differentials on phenotypic traits that are influenced by additive genetic variation is beyond doubt (Heino and Godø 2002; Swain et al. 2007). Yet, the rate of these changes and their consequences for stock viability, stability, yield, and recovery are less clear (Law 2000; Hutchings and Fraser 2008). Survivors of the harvesting process are likely to be genotypes with traits that confer relatively high fitness under fishing selection, but may be less than optimal with respect to natural selection (Conover 2000; Carlson et al. 2007). This may lead to slow recovery when fishing mortality is relaxed. Moreover, because cessation of fishing does not automatically produce equal selection pressures in the opposite direction, paying off this ‘Darwinian debt’ (Cookson 2004) may take a long time (Conover 2000; Law 2000).

Salmonids are well known for their diversity of life-history forms, with alternative mating tactics such as early maturing jacks in coho salmon, Oncorhynchus kisutch (Gross 1985), precocious parr in Atlantic Salmon, Salmo salar (Hutchings and Myers 1988), or various benthic and pelagic morphs in Artic charr, Salvelinus alpinus, and brook charr, Salvelinus fontinalis (Skúlason et al. 1996; Proulx and Magnan 2004). We follow Gross (1996) and Gross and Repka (1997, 1998) in using the term life-history tactics to refer to outcomes of life-history strategies, or decision rules, that determine how somatic and reproductive effort is allocated among alternative phenotypes. A common feature of many salmonid systems is the presence in sympatry of both anadromous (sea-run) and resident males and females, with resident fish completing their entire life cycle without migrating to sea (Jonsson and Jonsson 1993). Accumulating evidence suggests that these two forms may occur as alternative tactics within a single breeding population (Nordeng 1983; Morita et al. 2000; Olsson and Greenberg 2004; Thériault et al. 2007a). Individuals are understood to adopt a particular migration tactic by following a conditional life-history strategy involving energy thresholds, and various components of the energetic state of individuals (growth, lipid deposition, and metabolic rate) have been implicated in this process (Thorpe 1986; Bohlin et al. 1990; Hutchings and Myers 1994; Thorpe et al. 1998; Forseth et al. 1999; Morinville and Rasmussen 2003). Although influenced by environmental conditions (e.g. Olsson et al. 2006), the adoption of alternative life-history tactics in salmonids involves significant additive genetic variation, which has been demonstrated both in the laboratory (Silverstein and Hershberger 1992; Heath et al. 1994; Wild et al. 1994) and in the field (Garant et al. 2003; Thériault et al. 2007b). Moreover, whether an individual migrates or not will have critical consequences for its growth, survival, maturation, and reproduction. Survival is elevated in freshwater, but growth rates are reduced and resident individuals attain a smaller size at maturation (Gross 1987). As reproductive success is linked to body size in females (Fleming 1996; Morita and Takashima 1998; Thériault et al. 2007a), resident females experience decreased reproductive success relative to the bigger anadromous females. Reproductive success of males seems to be less affected by smaller size, as resident males employ alternative reproductive tactics, such as sneaking, to get access to mating opportunities (Hutchings and Myers 1988; Fleming 1996).

Owing to its size selectivity and temporally variable nature (Ricker 1995; Quinn et al. 2007), commercial fishing in salmonids has been shown to impact several life-history traits such as growth, age and size at maturation, and run timing. However, despite wide commercial and recreational interests in salmonids, evidence of evolutionary change caused by salmonid fisheries is still mostly circumstantial (Myers et al. 1986; Fukuwaka and Morita 2008, Hard et al. 2008). A fishery that targets only the migrant part of a population will inevitably be selective with respect to life-history tactics such as anadromy and residency. However, the consequences of such differential fishery-induced selection on the evolution of alternative life-history tactics have, to our knowledge, never been rigorously investigated.

Here, we used a recently developed modeling approach to predict the consequences of fishery-induced mortality on the evolution of anadromy and residency. The modeling approach (hereafter termed ‘eco-genetic’) incorporates both ecological and quantitative genetic processes, providing a mechanistically rich framework in which to predict the rate of evolutionary change on ecological timescales (Dunlop et al. 2007). In particular, our modeling approach enables distinguishing between plastic and evolutionary responses to fishing.

In the wild, salmonids show phenotypic plasticity in the age and size at migration. To account for such plasticity in the process of migration, we adopted a reaction norm approach. Reaction norms in the narrow sense describe how a single genotype is translated into different phenotypes depending on environmental conditions (Stearns 1992), while estimations of reaction norms in field studies must typically rely on the broader notion of population-level reaction norms (Sarkar and Fuller 2003; Hutchings 2004). Alternative tactics in salmonids have previously been described by reaction norms, based on the idea that the adoption of a particular tactic is governed by thresholds in growth rate (Myers and Hutchings 1986; Thorpe 1986; Bohlin et al. 1990; Hazel et al. 1990). Here, we extend this approach and consider the probability for the adoption of a particular migration tactic (anadromy or residency) as a function of size-at-age, where size-at-age has the helpful feature of integrating all environmental factors affecting growth. Such an approach has previously been used to model the evolution of maturation reaction norms (Ernande et al. 2004; Dunlop et al. 2007) and to tease apart phenotypically plastic responses from possible genetic changes in the age and size at maturation (Heino et al. 2002a; Grift et al. 2003; Barot et al. 2004, 2005; Olsen et al. 2004, 2005; Dunlop et al. 2005; Dieckmann and Heino 2007). Our study represents an extension of the maturation reaction norm approach so as to account for phenotypic plasticity in another fundamental life-history transition in the study of exploited populations.

We used data from a well-studied brook charr population in Québec, Canada, to parameterize our model. Recreational fishermen in the region are increasingly exploiting the sea-run components of this species as a result of the decline in Atlantic salmon stocks. Yet, anadromous populations of brook charr are not rigorously managed in many systems. Here, we examine the impact of various exploitation rates on the evolution of migration reaction norms, as well as on ecological and demographic characteristics of the population. We chose to model dynamics over a 100-year time horizon (approximately 30 generations) as this timeframe is commonly viewed as a manageable window from a conservation standpoint (Frankham et al. 2002).

The main purpose of our study was to address the following two questions: (i) Is fishing expected to induce evolutionary changes in the conditional migration strategy of salmonids? (ii) In a population with an evolving migration strategy, what are the effects of fishing on fecundity, abundance, and fishery yield? We also explored whether different freshwater mortality rates counteract or exacerbate the impact of fishing in saltwater on the evolution of anadromy and residency.

Methods

We constructed an individual-based eco-genetic model similar to that developed by Dunlop et al. (2007) to evaluate the effects of selective fishing mortality on the evolution of anadromy and residency within a conditional strategy framework. The model was built to reflect the life history of a sympatric population of brook charr in which anadromous and resident migration tactics coexist (Fig. 1) inhabiting a small tributary of the Ste-Marguerite River in Québec, Canada, named Morin Creek. The behavior and life history of brook charr in this system are well studied (Morinville and Rasmussen 2003, 2006, 2007; Thériault and Dodson 2003; Lenormand et al. 2004; Thériault et al. 2007a,b) and ample data from the years 1998–2004 were available to parameterize the model (Fig. 2, Table 1). The model follows evolution of the migration reaction norm, a quantitative trait that is passed on at the individual-level from parents to offspring. We assumed a closed population such that no new genetic variance was introduced by immigration. Model simulations were run for a total duration of 100 years in discrete, 1-year time steps and each simulation was repeated 30 times. As we consider the rate and amount of evolutionary change over the course of 100 years, our predictions represent transients and not evolutionary endpoints. Each year, individuals had the chance to experience the processes of migration to and from saltwater, growth, maturation, reproduction, and mortality (Fig. 1).

Figure 1

Schematic overview of the life cycle of brook charr, showing the sequence of events in the eco-genetic model.

Figure 2

Empirically derived functions used in the model. (A) Probabilistic migration reaction norms estimated for Morin creek (data from Thériault and Dodson 2003). (B) Probabilistic maturation reaction norms estimated for anadromous individuals (thick lines) and resident individuals (thin lines), in each case showing the 50%, 1%, and 99% probability percentile curves. (C) Relationship between fecundity and body length for anadromous individuals (open circles) and resident individuals (filled circles) (data from Lenormand 2003). (D) Stock-recruitment relationship (based on Elliott 1993). (E) Selectivity curves for different maximal harvest probabilities increasing from 5% to 100% in increments of 10%.

Table 1

Model parameters and their values.

Symbol	Description	Equation	Source	Value
–	Initial mean body length (mm)	–	1	82.98
–	Initial standard deviation of body length (mm)	–	1	13.68
–	Mean emergence size (mm)	–	1	31.71
–	Standard deviation of emergence size (mm)	–	1	5.52
h2	Initial heritability of evolving PMigRN traits	–	–	0.5
CV	Coefficient of variation for all evolving traits	–	–	0.08
c0	Evolving PMigRN trait	1	1	−12.36
c1	Evolving PMigRN trait (mm−1)	1	1	0.11
c2	Evolving PMigRN trait (year−1)	1	1	9.69
c3	Evolving PMigRN trait (mm−1 year−1)	1	1	−0.08
gs	Mean annual growth in saltwater (mm)	2	2	95.53
gf	Mean annual growth in freshwater (mm)	2	1	35.35
GSI	Mean gonado-somatic index	2	2	0.147
Sf	PMRN slope of resident morph (mm year−1)	3	1	−32.41
if	PMRN intercept of resident morph (mm)	3	1	259.72
wf	PMRN width of resident morph (mm)	3	1	114.53
Ss	PMRN slope of anadromous morph (mm year−1)	3	2	−177.58
is	PMRN intercept of anadromous morph (mm)	3	2	843.38
ws	PMRN width of anadromous morph (mm)	3	2	532.9
H1	Constant in fecundity function (mm−1)	4	2	0.04
H2	Constant in fecundity function	4	2	2.86
r	Constant in stock-recruitment function	5	3	25.0
b	Constant in stock-recruitment function	5	3	0.0027
mi,r	Immature natural mortality probability for resident morph under default conditions	–	2	0.60
mm,r	Mature natural mortality probability for resident morph under default conditions	–	2	0.88
mi,a1	Immature natural mortality probability in 1st year for anadromous morph under default conditions	–	2	0.80
mi,a2	Immature natural mortality probability in 2nd year for anadromous morph under default conditions	–	2	0.60
mm,a	Mature natural mortality probability for anadromous morph under default conditions	–	2	0.55
mi,r,poor	Immature natural mortality probability for resident morph under poor conditions	–	2	0.80
mm,r,poor	Mature natural mortality probability for resident morph under poor conditions	–	2	0.95
mi,r,good	Immature natural mortality probability for resident morph under good conditions	–	2	0.20
mm,r,good	Mature natural mortality probability for resident morph under good conditions	–	2	0.70

PMRN, probabilistic maturation reaction norm; PMigRN, probabilistic migration reaction norm.

Data sources: (1) Morin Creek data from Thériault (2001) and Thériault and Dodson (2003), (2) Ste-Marguerite River data for anadromous fish and Morin Creek data for resident fish from Lenormand (2003), (3) based on Elliott (1993).

Migration

The migration reaction norm was represented by a logistic function, describing the probability P of migrating as a function of age a and body length l,

where . This form was chosen because it allows the probability of migrating as a function of size to change its slope with age (Fig. 2A). Each individual was thus characterized by the four evolving parameters c0, c1, c2, and c3 describing its probabilistic migration reaction norm (PMigRN), which, together with its age and length, in turn determined its probability of migrating in a given year. An individual could migrate at either age 1 or 2 years only: if a fish did not migrate by age 2, we assumed that it would be a freshwater resident for all its life. This understanding is corroborated by field observations on this system (Thériault and Dodson 2003).

Somatic growth

Individuals grew according to the growth model introduced by Lester et al. (2004). Newborns in the model were given a random size at emergence, in accordance with the empirical mean and standard deviation estimated from back-calculations of the 1998-, 1999-, and 2000-year classes of fish captured in Morin creek (Thériault 2001; Thériault and Dodson 2003). Prior to maturation, individuals grew with an annual phenotypic growth increment g determined by the environment in which they resided during that year (freshwater, e = f, or saltwater, e = s). The maximal growth increment is expressed in saltwater, whereas individuals living in freshwater grow slower due to the poorer growing environment they experience. The environment-specific growth rates gf and gs were empirically derived from immature individuals of the Ste-Marguerite River system and Morin Creek (Lenormand 2003; Thériault and Dodson 2003; Table 1).

Immediately following maturation, individuals devoted a proportion of energy to reproductive tissues, so that the body length l at age a + 1 was given by

where GSI was the gonado-somatic index (gonadic mass divided by somatic mass) estimated for anadromous females. The gonado-somatic index was assumed to be similar and constant for all mature individuals in the population for simplicity. Growth rates were assumed to be density-independent, both in saltwater and in freshwater, to keep predictions simple. This simplifying assumption was further motivated by the following two reasons. First, in view of the high productivity of marine habitats, the small population sizes modeled here, and the importance of density-dependent predation mortality in saltwater, the density dependence of growth rates at sea must be expected to be weak. Second, results gathered from a creek adjacent and similar to Morin Creek failed to detect any density-dependence in freshwater growth for brook charr of age 0 and older (Centre Interuniversitaire de Recherche sur le Saumon Atlantique, CIRSA, unpublished data). When we tested the effects of relaxing a simplifying model assumption by adding density-dependent freshwater growth, there was little impact on the probability of migrating at age 1, but the probability of migrating at age 2 did evolve to be higher when the strength of density dependence was increased (Supplementary material, Section 3). We also evaluated the sensitivity of our results with regard to relaxing the simplifying assumption that growth rates were not evolving (Supplementary material, Section 2). We could thus confirm the robustness of our results under the incorporation of growth evolution.

Maturation

In any given year, an immature individual had a probability to mature during the upcoming year that was based on its environment (freshwater, e = f; or saltwater, e = s) and on its age a and body length l. These probabilities Pm were given by probabilistic maturation reaction norms (PMRNs, Fig. 2B, Heino et al. 2002a; Dieckmann and Heino 2007) with logistic length dependence, linear age dependence, and constant width,

where i, s, and w were the environment-specific PMRN intercepts, slopes, and widths, respectively, and k = logit(99%) − logit(1%). To keep the model simple, and to focus on the evolution of migration, maturation tendency was not considered as an evolving trait in our model.

A logistic regression was applied to age-specific length distributions of immature and mature fish to provide an approximate estimation of the population’s PMRNs (Heino et al. 2002a; Dieckmann and Heino 2007). For anadromous fish, we used data gathered from the whole Ste-Marguerite River system (pooled from 1998 to 2001), whereas for resident fish, we used data gathered from Morin creek (pooled from 1998 to 2002). A linear regression of lengths at ages 2 and 3 (the two age classes for which sufficient data were available) at which the probability to mature in the next year was 50% was used to estimate the slopes and intercepts of linear PMRNs (Fig. 2B). We then used the 1% and 99% maturation probability percentiles of 3-year-old individuals to determine the PMRN widths. Our estimated PMRNs are only an approximation because, due to sampling constraints, we included fish that could have matured in a previous year in our mature length distribution. However, maturation occurs over such a narrow range in this population (the majority of individuals mature between ages 2 and 3), that the effect of including previously matured fish should be minimal relative to a population in which the range of maturation ages is larger. Owing to the uncertainty in our approximated PMRN, we performed sensitivity analyses and found that varying the PMRN slope and intercept had little impact on our model results (Supplementary material, Section 4).

Reproduction

There was no sex-structure in our model and reproduction occurred annually in freshwater between random pairs of mature individuals. The largest individual in the reproductive pair was chosen to be the mother, so as to account for frequently observed mating between anadromous females (bigger) and resident males (smaller) and the apparent absence of the reverse (big anadromous males are not expected to mate with small resident females, Thériault et al. 2007a). The number of eggs produced by a reproductive pair was estimated from the body length l of the mother according to an empirically derived relationship between fecundity and body length (Fig. 2C):

with allometric constants H1 and H2. The number of new individuals recruiting to the population at age 1 was determined from a Ricker stock-recruitment function (Fig. 2D):

where S is the number of adults, and r and b are constants. As the necessary data to derive such a stock-recruitment function specifically for our system were not available, we used constants estimated for brown trout (Elliott 1993), a species with life-history characteristics very similar to brook charr, and scaled them so as to yield realistic estimates of recruitment and spawner abundance in our system (Table 1).

Inheritance and expression

Genotype determination

Inheritance of the PMigRN was described by the infinitesimal model of quantitative genetics (Cavalli-Sforza and Feldman 1976). We assumed that phenotypic plasticity for migration was heritable by modeling genetically based reaction norms that were passed from parents to offspring (e.g. Brommer et al. 2005; Nussey et al. 2005; Dunlop et al. 2007). The four parameters c0, c1, c2, and c3 describing the PMigRN were thus considered as evolving traits. The genetic trait values of an offspring were drawn at random from normal distributions with mean values given by the mid-parental genetic trait values and variances that equaled half the corresponding genetic variances in the initial population. Modeling offspring variance in this way assumes equal variances of maternal and paternal traits and that the segregation and recombination of genes during reproduction introduce a constant amount of variation into the population (Roughgarden 1979).

Phenotype determination

The phenotypically expressed values of an individual’s four PMigRN traits were drawn randomly in each year from normal distributions with mean values given by the individual’s genetic trait values and variances that equaled the assumed environmental variances. The latter were calculated based on an assumed initial heritability h (see section on Initial population structure below) for the PMigRN traits and on an assumed initial genetic coefficient of variation. Based on the definition of heritability, h2 = VA/VP with VP = VA + VE (VA is the additive genetic variance; VP, total phenotypic variance; and VE, environmental variance), it is possible to calculate VE for each trait from the initial values of h2 and VA. While the environmental variances for each trait were kept constant in the model, the corresponding values of VA, and thus of h2, were free to evolve after the initial year (Supplementary material, Section 1).

Natural mortality

Default natural mortalities

Age-specific annual mortality probabilities were estimated from Morin Creek data for resident individuals, and from a larger mark-recapture experiment in the whole Ste-Marguerite River system for anadromous individuals (Lenormand 2003). Immature (mi) and mature (mm) mortality probabilities were applied annually to resident and anadromous individuals (Table 1). For anadromous fish, the mortality probability of an immature individual varied depending on whether it was the first or second year the individual spent in saltwater (mi,a1 and mi,a2, respectively, Table 1).

Alternative natural mortalities

In addition to the default values representing ‘normal’ freshwater mortality probabilities, we also simulated ‘poor’ freshwater survival conditions and ‘good’ freshwater survival conditions (Table 1), while keeping natural saltwater mortalities unchanged.

Fishing mortality

We applied fishing to anadromous individuals only, as they are the only targets of the recreational fishery. Length-dependent annual harvest probabilities for these fish were derived based on the observed sizes of fish caught and on data quantifying overall annual exploitation rates (Lenormand 2003; CIRSA, unpublished data, Fig. 2E). Medium-sized anadromous fish (with body lengths between 200 and 350 mm) are most likely to be caught, because they are abundant and, during the upstream migration of immature anadromous brook charr in early fall (Lenormand et al. 2004), concentrated in the river’s estuary, where their exploitation is little regulated. Smaller brook charr (with lengths between 110 and 200 mm) are not attractive to fishermen, whereas the bigger, mostly mature charr (with lengths larger than 350 mm) are under spatial and temporal regulations that prevent high fishing pressures on these larger fish. Size-selective fishing mortality was applied to individuals regardless of their maturation status. We varied the maximal harvest probabilities in the selectivity curves of anadromous fish between 0% and 100% in increments of 10% (Fig. 2E).

Initial population structure

Initial length distribution

The initial population in the model consisted of 5000 age-1 individuals with initial body lengths following a normal distribution with mean and standard deviation estimated for the 1998- to 2000-year classes of fish captured in Morin creek (Thériault and Dodson 2003; Table 1).

Initial migration reaction norms

The initial population-level PMigRN was estimated using data on size and age at migration from Morin Creek. Data on fish of ages 1 and 2 were analyzed for the years 1998–2000, as migration occurs almost exclusively at these two ages (Thériault and Dodson 2003). All individuals from the initial population were assigned genetic values for the four evolving traits c0, c1, c2, and c3 following normal distributions with mean values given by the trait values implied by the initial population-level PMigRN (Fig. 2A) and standard deviations given by the assumed initial genetic coefficient of variation.

Initial heritabilities

The initial heritability of each trait describing the PMigRN was assumed as 0.5. We do not know the actual value of heritability of plasticity for anadromy and residency in this system, but genetic variation and heritability have been demonstrated for plasticity in general (Scheiner 1993; Nussey et al. 2005) and have been assumed for migration in salmonids in particular (Hazel et al. 1990; Hutchings and Myers 1994; see also the review by Hutchings 2004). We chose the value 0.5 because the heritability of threshold traits that influence the adoption of alternative migration tactics varies between 0.52 and 0.56 for brook charr in this system (Thériault et al. 2007b), and between 0.12 and 0.98, with a mean of 0.53 for various other threshold traits (heritability of binomial threshold traits are reported on a ‘liability’ scale; Roff 1996). After initialization, heritabilities, genetic variances, and genetic covariances were free to evolve, and can thus be regarded as emerging properties of the model. Even though heritabilities directly scale the speed of evolution, so that we must expect slower or faster changes in reaction norms if we assume lower or higher heritabilities, the nature of predicted evolutionary changes remains unchanged as heritabilities are jointly increased or decreased (see, e.g., Dunlop et al. 2007).

Results

Impact of different fishing mortalities

Increasing the probability of harvest causes an evolutionary shift in the migration reaction norms for both age-1 and age-2 individuals: for an individual of the same size, the probability of migration is lowered as the maximal harvest probability increases (Fig. 3A,B). This translates into an overall probability of migration that is decreasing with increasing maximal harvest probability (Fig. 3C). The absolute number of fish that migrate decreases as the maximal harvest probability increases, and this trend is more pronounced for age-1 individuals than for age-2 individuals (Fig. 3D). Mean age at migration increases with maximal harvest probability (Fig. 3E), primarily reflecting the fact that the proportion of fish migrating at age 2 increases with harvest probability. Mean age at maturation did not change for residents, but decreased for anadromous individuals (Fig. 3F). Mean individual fecundity, highly dependent on size, decreased with the maximal harvest probability for anadromous fish, but showed no change for resident fish (Fig. 3G). Overall abundance of the population shows little change with increasing maximal harvest probability, because the number of fish in freshwater rose while the number in saltwater decreased to almost zero (Fig. 3H). The heritability of the migration reaction norm traits varied though time, but did not show a significant increase or decrease, either at low or at high maximal harvest probability (Fig. S1).

Figure 3

Model results after 100 years of fishing with different maximal harvest probabilities. Panels (A) and (B) show the resultant age-specific migration reaction norms (line thickness increases with increasing maximal harvest probability between 0% and 100% in increments of 10%). Panels (C–H) show how maximal harvest probabilities affect age-specific migration probabilities, age-specific numbers of migrating fish, ages at migration and maturation, individual fecundity, and population abundances. Results are averaged over 30 independent model runs.

Impact of different natural mortalities

Survival conditions in freshwater influenced the evolution of the migration reaction norm. After 100 years of fishing, low survival in freshwater associated with poor conditions leads to the evolution of a migration reaction norm that implies a higher probability of migrating for a given size than the reaction norm evolved under normal freshwater survival (Fig. 4A,B,C). In contrast, good survival conditions in freshwater leads to a lower probability of migrating for a given size than under normal conditions (Fig. 4A,B,C). Poor survival in freshwater thus offsets the evolutionary effect of fishing by increasing the probability of migrating (Fig. 4C) and the number of migrants (Fig. 4D), whereas good survival in freshwater had the opposite effect. The population’s abundance as a function of the maximal harvest probability is higher for poor than for normal freshwater survival conditions, whereas for good conditions, a nonmonotonic dependence on maximal harvest probability is found (Fig. 4E). The cumulative catch shows similar dome-shaped relationships for the three survival conditions in freshwater, but peaks at higher catches and higher maximal harvest probabilities as the survival conditions in freshwater are worsened (Fig. 4F).

Figure 4

Model results after 100 years of fishing and different survival conditions in freshwater (poor, normal, or good). The maximum harvest probability was 0.5. Panels (A) and (B) show the resultant age-specific migration reaction norms. Panels (C–F) show how maximal harvest probabilities affect migration probabilities (for ages 1 and 2 combined), numbers of migrating fish (for ages 1 and 2 combined), total population abundances, and cumulative catches. Results are averaged over 30 independent model runs.

Discussion

Using an eco-genetic modeling approach, we explored the impact of recreational fishing on the evolution of anadromy and residency of a small population of brook charr. During a 100 years of fishing of anadromous individuals, we predicted evolution in this population’s migration reaction norm, with the average probability of migration decreasing with increasing harvest rate. These changes were accompanied by increases in the proportion of fish migrating at age 2, resulting in higher mean ages at migration. Our findings suggest that selective harvesting of anadromous fish results in a higher tendency for residency, through an increased fitness advantage of staying longer in freshwater and delaying migration. Shifts in the maturation reaction norms of several commercially important marine species, toward younger ages and smaller sizes at maturation, have been reported in the wake of heavy fishing pressures that selected against genotypes predisposing fish to mature later and larger (Atlantic cod, Gadus morhua, Heino et al. 2002b; Barot et al. 2004; Olsen et al. 2004, 2005; Baulier et al. 2006; plaice Pleuronectes platessa, Grift et al. 2003, 2007; American plaice, Hippoglossoides platessoides, Barot et al. 2005; sole, Solea solea, Mollet et al. 2006). Here, we have shown that evolutionary shifts in the reaction norm of another fundamental ontogenetic process – i.e., migration, an important life-history characteristic in salmonids – are also expected to result from elevated fishing mortality.

By changing the distribution of heritable traits, harvesting by humans can unintentionally select against the most desirable phenotypes (i.e. bigger individuals), which increase harvestable biomass (Law and Grey 1989; Coltman 2008; Hutchings and Fraser 2008). In a terrestrial context, trophy hunting of bighorn sheep over a period of 30 years has generated an undesired evolutionary response in weight and horn length, shown by the decrease in the mean breeding values of both of these traits (Coltman et al. 2003). Here, we have shown that the intensive harvesting of anadromous fish evolutionarily reduces the probability of migration. This ultimately leads to a reduced number of fish in saltwater, and thus to a situation in which less fish are available to the recreational fishery. We also found changes in the age at maturation for the anadromous part of the population. These changes, however, must not be interpreted as genetic changes, as the maturation tendency was not allowed to evolve in the model. Instead, these changes reflect the fact that at high harvest rates only the smallest and youngest anadromous fish are escaping harvest (as a result of the size-dependent harvest probabilities shown in Fig. 2E). The removal of large fish by the recreational fishery also causes a decline in the mean individual fecundity of anadromous fish at high harvest probabilities.

Genetic changes induced by fisheries, as any other selective force, are potentially reversible, given that sufficient heritable genetic variation remains and that adequate selection differentials in the opposite direction are generated once fishing is relaxed or stopped (Law 2000; Reznick and Ghalambor 2005). According to the body of work on commercially exploited marine species, however, it appears that such reversal is a difficult and slow process (Barot et al. 2004; Olsen et al. 2004; Swain et al. 2007; but see Fukuwaka and Morita 2008, who show an increase in the maturation threshold following cessation of fishing). Although we have not explored the extent of trait reversal following the cessation of fishing, our results suggest that heritable additive genetic variance in the migration reaction norm is preserved even at high harvest rates. The maintenance of genetic variation even under strong directional selection has been demonstrated for various traits (Houle 1992), including threshold traits (Roff 1994). In the latter case, variation remains ‘hidden’ by virtue of the threshold nature of the trait. This preservation of genetic variation may imply that evolutionary changes in threshold traits could be reversed after the return of favorable conditions. Lost alternative life-history tactics could be restored provided that the mechanisms for their phenotypic expression have not degenerated during a period of disuse (West-Eberhard 2003). Indeed, cases have been documented of nonanadromous salmonid fish stocks that maintained their capacity to migrate or that again gave rise to anadromous morphs after transplantation (Staurnes et al. 1992; Thrower et al. 2004).

Our study also aimed at assessing how variations in natural survival in freshwater could influence the evolution of migration. We found that shifts in migration reaction norms caused by selective harvesting were either impeded or exacerbated by low or high survival in freshwater, respectively. Fluctuating selection on body size has also been hypothesized as a factor favoring threshold variation, and thus ultimately maintaining alternative male life cycles, in Atlantic salmon (Aubin-Horth et al. 2005). Fluctuations in the direction of evolution have been demonstrated on the time scale of decades in Darwin finches, where selection on body size and beak shape has changed direction in time (Grant and Grant 2002). In the face of high temporal variability in ecological conditions, typical of northern temperate rivers inhabited by salmonids, predicting the rate and magnitude of fishing-induced life-history evolution thus becomes difficult over the long term.

This study deliberately focused on fishing-induced changes in the migration reaction norm of salmonids, as this phenomenon has received little, if any, attention in the published literature. To keep predictions simple, we did not allow for concomitant evolution of other life-history traits, such as growth capacity, the maturation reaction norm, or reproductive investment. While these other traits are also likely to evolve, most models published to date have focused, like our study here, on fishing-induced evolution in only one such character, such as maturation age or size, or the maturation reaction norm (Heino 1998; Ernande et al. 2004; Gårdmark and Dieckmann 2006; Dunlop et al. 2007). With this study being the first to address the fishing-induced evolution of a PMigRN, the inclusion of additional evolving traits would have unduly complicated the model, at least at this stage of investigation. We also confirmed that including growth capacity as an evolving trait did not significantly alter our predictions about migration evolution (Fig. S2). Future extensions of our model could explore the evolution of additional traits.

The effects of correlated responses to selection would also merit further investigation. Fishing-induced selection on one trait could generate responses in other genetically correlated traits: this could either diminish or amplify rates of evolutionary change, depending on the sign and magnitude of the correlation (Lynch 1999; Walsh et al. 2006; Dieckmann and Heino 2007; Hutchings and Fraser 2008). For example, smolting, maturation, and growth have all been shown to be genetically correlated, to varying degrees, in rainbow trout (Oncorhynchus mykiss,Thrower et al. 2004). Furthermore, the dynamic interactions of these traits with season-specific growth rates have been hypothesized to be a key factor maintaining genetic variation in smolting, despite complete selection against this phenotypic expression of migration (Thrower et al. 2004). Genes associated with smolting are thus conserved in the population, through selection for, or against, other genetically correlated traits. In brook charr, the probability of migrating at age 1 was found to be positively genetically correlated with body size at age 1 (Thériault et al. 2007b). One might expect that such a correlation would exacerbate the effect of fishing: selecting against anadromous individuals migrating at age 1 might indirectly select for smaller size at age 1, which in turn is associated with a higher probability of staying resident. However, a smaller size at age 1 could also translate into a smaller size at age 2, implying a higher probability of migrating. The effects of genetic correlations among traits related to the life-history tactics adopted in the face of fishing-induced selection remain to be explored.

Despite their important commercial and recreational implications, the effects of fishing on salmonids have rarely been demonstrated beyond the immediate consequences for abundance (see review by Hard et al. 2008). Genetic responses to commercial fisheries have been proposed for explaining a decrease in size of Pacific salmon (Ricker 1995), as well as changes in age and size at maturation in European grayling, Thymallus thymallus (Haugen and Vøllestad 2001), lake whitefish, Coregonus clupeaformis (Handford et al. 1977), and Atlantic salmon, Salmo salar (Bielak and Power 1986). However, these studies relied on indirect evidence as genetic changes were not demonstrated directly. By illustrating the impact of a sport fishery on the evolution of life-history tactics in salmonids, along with the associated ecological and demographic consequences, our present study contributes to the emerging view that evolution can be rapid enough to be an integral part of ecological interactions (Hendry and Kinnison 1999; Reznick and Ghalambor 2005; Jørgensen et al. 2007). The high harvest probabilities that caused the most significant demographic and evolutionary changes in this study are more common in commercial salmonid fisheries (such as in Atlantic salmon; Dempson et al. 2001) than in recreational fisheries (such as in the brook charr system studied here, where the mean harvest probability equaled 0.3; CIRSA, unpublished data). Recreational practices, however, still have the potential to contribute to the decline of fisheries (Cooke and Cowx 2006). We therefore recommend that managers and policy-makers increasingly complement traditional methods of fisheries management based on population dynamics alone with modeling approaches of the kind presented here.