State-dependent life history models in a changing (and regulated) environment: steelhead in the California Central Valley

Age at release affects developmental physiology and sex-specific phenotypic diversity of hatchery steelhead trout (Oncorhynchus mykiss)

Most steelhead trout hatcheries increase growth rate during rearing to produce and release yearling smolts for harvest augmentation, but natural steelhead exhibit variable age of smoltification, so this common rearing practice may not be ideal for programs focused on recovering imperiled wild stocks; therefore, it is important to investigate and compare alternative hatchery rearing methods that promote life history diversity. Over six consecutive years, the Winthrop National Fish Hatchery on the Methow River, WA reared and released paired groups of age-1 (S1) and age-2 (S2) steelhead smolts. To understand how the two rearing methods affected developmental ontogeny and life-history, fish were sampled prior to hatchery release for factors associated with smoltification (size, gill Na+/K+ ATPase activity, and a qualitative smolt phenotype) and sexual maturation (sex, pituitary and testis mRNA transcripts, gonadosomatic index, and plasma 11-ketotestosterone). Our objectives were to quantify levels of smoltification and male maturation during hatchery rearing, combine metrics to estimate residualism (failure to migrate upon release), and compare the treatments by sex. Overall, S2 rearing produced 7.8% more smolts and 44-fold (4.4 vs. 0.1%) more precociously mature males than S1 rearing. Conversely, S1 rearing produced 31.6% more residuals than S2 rearing. While the proportion of total male residuals was comparable between treatments, the S1 treatment produced approximately five-fold more female residuals (20.6 vs. 4.2%). Because residuals contribute minimally to adult returns and the number of returning adult females is critical to the success of salmonid supplementation efforts, developing rearing techniques that maximize migration in females is a management priority. Physiological assessments are useful for characterizing and quantifying the effects and risks of different hatchery rearing regimes on steelhead life-history, in addition to providing sex-specific guidance to inform and optimize conservation management goals in supplementation programs.

Genetic parentage reveals the (un)natural history of Central Valley hatchery steelhead

Populations composed of individuals descended from multiple distinct genetic lineages often feature significant differences in phenotypic frequencies. We considered hatchery production of steelhead, the migratory anadromous form of the salmonid species Oncorhynchus mykiss, and investigated how differences among genetic lineages and environmental variation impacted life history traits. We genotyped 23,670 steelhead returning to the four California Central Valley hatcheries over 9 years from 2011 to 2019, confidently assigning parentage to 13,576 individuals to determine age and date of spawning and rates of iteroparity and repeat spawning within each year. We found steelhead from different genetic lineages showed significant differences in adult life history traits despite inhabiting similar environments. Differences between coastal and Central Valley steelhead lineages contributed to significant differences in age at return, timing of spawning, and rates of iteroparity among programs. In addition, adaptive genomic variation associated with life history development in this species varied among hatchery programs and was associated with the age of steelhead spawners only in the coastal lineage population. Environmental variation likely contributed to variations in phenotypic patterns observed over time, as our study period spanned both a marine heatwave and a serious drought in California. Our results highlight evidence of a strong genetic component underlying known phenotypic differences in life history traits between two steelhead lineages.

Environmentally triggered shifts in steelhead migration behavior and consequences for survival in the mid-Columbia River

The majority of Columbia River summer-run steelhead encounter high river temperatures (near or > 20°C) during their spawning migration. While some steelhead pass through the mid-Columbia River in a matter of days, others use tributary habitats as temperature refuges for periods that can last months. Using PIT tag detection data from adult return years 2004-2016, we fit 3-component mixture models to differentiate between "fast", "slow", and "overwintering" migration behaviors in five aggregated population groups. Fast fish migrated straight through the reach on average in ~7-9 days while slow fish delayed their migration for weeks to months, and overwintering fish generally took ~150-250 days. We then fit covariate models to examine what factors contributed to the probability of migration delay during summer months (slow or overwintering behaviors), and to explore how migration delay related to mortality. Finally, to account for the impact of extended residence times in the reach for fish that delayed, we compared patterns in estimated average daily rates of mortality between migration behaviors and across population groups. Results suggest that migration delay was primarily triggered by high river temperatures but temperature thresholds for delay were lowest just before the seasonal peak in river temperatures. While all populations groups demonstrated these general patterns, we documented substantial variability in temperature thresholds and length of average delays across population groups. Although migration delay was related to higher reach mortality, it was also related to lower average daily mortality rates due to the proportional increase in reach passage duration being larger than the associated increase in mortality. Lower daily mortality rates suggest that migration delay could help mitigate the impacts of harsh migration conditions, presumably through the use of thermal refuges, despite prolonged exposure to local fisheries. Future studies tracking individual populations from their migration through reproduction could help illuminate the full extent of the tradeoffs between different migration behaviors.

Eco-evolutionary dynamics may show an irreversible regime shift, illustrated by salmonids facing climate change

The enhanced or reduced growth of juvenile masu salmon (Oncorhynchus masou masou) may result from climate changes to their environment and thus impact on the eco-evolutionary dynamics of their life-history choices. Male juveniles with status, i.e., if their body size is larger than a threshold, stay in the stream and become resident males reproducing for multiple years, while those with smaller status, i.e., their body size is below the threshold, migrate to the ocean and return to the stream one year later to reproduce only once. Since juvenile growth is suppressed by the density of resident males, the fraction of resident males may stay in equilibrium or fluctuate wildly over a 2-year period. When the threshold value evolves, the convergence stable strategy may generate either an equilibrium or large fluctuations of male residents. If environmental changes occur faster than the rate of evolutionary adaptation, the eco-evolutionary dynamics exhibit a qualitative shift in the population dynamics. We also investigated the relative assessment models, in which individual life-history choices are made based on the individual’s relative status within the juvenile population. The eco-evolutionary dynamics are very different from the absolute assessment model, demonstrating the importance of understanding the mechanisms of life history choices when predicting the impacts of climate change.

Biological and statistical interpretation of size-at-age, mixed-effects models of growth

The differences in life-history traits and processes between organisms living in the same or different populations contribute to their ecological and evolutionary dynamics. We developed mixed-effect model formulations of the popular size-at-age von Bertalanffy and Gompertz growth functions to estimate individual and group variation in body growth, using as a model system four freshwater fish populations, where tagged individuals were sampled for more than 10 years. We used the software Template Model Builder to estimate the parameters of the mixed-effect growth models. Tests on data that were not used to estimate model parameters showed good predictions of individual growth trajectories using the mixed-effects models and starting from one single observation of body size early in life; the best models had R 2 > 0.80 over more than 500 predictions. Estimates of asymptotic size from the Gompertz and von Bertalanffy models were not significantly correlated, but their predictions of size-at-age of individuals were strongly correlated (r > 0.99), which suggests that choosing between the best models of the two growth functions would have negligible effects on the predictions of size-at-age of individuals. Model results pointed to size ranks that are largely maintained throughout the lifetime of individuals in all populations.

Anadromy, potamodromy and residency in brown trout Salmo trutta: the role of genes and the environment

Brown trout Salmo trutta is endemic to Europe, western Asia and north-western Africa; it is a prominent member of freshwater and coastal marine fish faunas. The species shows two resident (river-resident, lake-resident) and three main facultative migratory life histories (downstream-upstream within a river system, fluvial-adfluvial potamodromous; to and from a lake, lacustrine-adfluvial (inlet) or allacustrine (outlet) potamodromous; to and from the sea, anadromous). River-residency v. migration is a balance between enhanced feeding and thus growth advantages of migration to a particular habitat v. the costs of potentially greater mortality and energy expenditure. Fluvial-adfluvial migration usually has less feeding improvement, but less mortality risk, than lacustrine-adfluvial or allacustrine and anadromous, but the latter vary among catchments as to which is favoured. Indirect evidence suggests that around 50% of the variability in S. trutta migration v. residency, among individuals within a population, is due to genetic variance. This dichotomous decision can best be explained by the threshold-trait model of quantitative genetics. Thus, an individual's physiological condition (e.g., energy status) as regulated by environmental factors, genes and non-genetic parental effects, acts as the cue. The magnitude of this cue relative to a genetically predetermined individual threshold, governs whether it will migrate or sexually mature as a river-resident. This decision threshold occurs early in life and, if the choice is to migrate, a second threshold probably follows determining the age and timing of migration. Migration destination (mainstem river, lake, or sea) also appears to be genetically programmed. Decisions to migrate and ultimate destination result in a number of subsequent consequential changes such as parr-smolt transformation, sexual maturity and return migration. Strong associations with one or a few genes have been found for most aspects of the migratory syndrome and indirect evidence supports genetic involvement in all parts. Thus, migratory and resident life histories potentially evolve as a result of natural and anthropogenic environmental changes, which alter relative survival and reproduction. Knowledge of genetic determinants of the various components of migration in S. trutta lags substantially behind that of Oncorhynchus mykiss and other salmonines. Identification of genetic markers linked to migration components and especially to the migration-residency decision, is a prerequisite for facilitating detailed empirical studies. In order to predict effectively, through modelling, the effects of environmental changes, quantification of the relative fitness of different migratory traits and of their heritabilities, across a range of environmental conditions, is also urgently required in the face of the increasing pace of such changes.

Size-conditional smolting and the response of Carmel River steelhead to two decades of conservation efforts

Threshold effects are common in ecosystems and can generate counterintuitive outcomes in management interventions. A threshold effect proposed for steelhead trout (Oncorhynchus mykiss) is size-conditional smolting and marine survival. Steelhead are anadromous, maturing in the ocean but migrating to freshwater to spawn, where their offspring reside for one or more years before smolting—physiologically transforming to a saltwater form—and migrating to the ocean. In conditional smolting, juveniles transform only if growth exceeds a threshold body size prior to migration season, and subsequent marine survival correlates with size at ocean entry. Conditional smolting suggests that efforts to improve freshwater survival of juveniles may reduce smolt success if they increase competition and reduce growth. Using model-selection techniques, we asked if this effect explained declining numbers of adult Carmel River steelhead. This threatened population has been the focus of two decades of habitat restoration, as well as active translocation and captive-rearing of juveniles stranded in seasonally dewatered channels. In the top-ranked model selected by information-theoretic criteria, adult decline was linked to reduced juvenile growth rates in the lower river, consistent with the conditional smolting hypothesis. According to model inference, since 2005 most returning adult steelhead were captively-reared. However, a lower-ranked model without conditional smolting also had modest support, and suggested a negative effect of captive rearing. Translocations of juvenile fish to perennial reaches may have reduced the steelhead run slightly by raising competition, but this effect is confounded in the data with effects of river flow on growth. Efforts to recover Carmel River steelhead will probably be more successful if they focus on conditions promoting rapid growth in the river. Our analysis clearly favored a role for size-conditional smolting and marine survival in the decline of the population, but did not definitively rule out alternative explanations.

Evolutionary restoration potential evaluated through the use of a trait-linked genetic marker

Human‐driven evolution can impact the ecological role and conservation value of impacted populations. Most evolutionary restoration approaches focus on manipulating gene flow, but an alternative approach is to manipulate the selection regime to restore historical or desired trait values. Here we examined the potential utility of this approach to restore anadromous migratory behavior in coastal California steelhead trout (Oncorhynchus mykiss) populations. We evaluated the effects of natural and anthropogenic environmental variables on the observed frequency of alleles at a genomic marker tightly associated with migratory behavior across 39 steelhead populations from across California, USA. We then modeled the potential for evolutionary restoration at sites that have been impacted by anthropogenic barriers. We found that complete barriers such as dams are associated with major reductions in the frequency of anadromy‐associated alleles. The removal of dams is therefore expected to restore anadromy significantly. Interestingly, accumulations of large numbers of partial barriers (passable under at least some flow conditions) were also associated with significant reductions in migratory allele frequencies. Restoration involving the removal of partial barriers could be evaluated alongside dam removal and fishway construction as a cost‐effective tool to restore anadromous fish migrations. Results encourage broader consideration of in situ evolution during the development of habitat restoration projects.

The puzzle of partial migration: Adaptive dynamics and evolutionary game theory perspectives

We consider the phenomenon of partial migration which is exhibited by populations in which some individuals migrate between habitats during their lifetime, but others do not. First, using an adaptive dynamics approach, we show that partial migration can be explained on the basis of negative density dependence in the per capita fertilities alone, provided that this density dependence is attenuated for increasing abundances of the subtypes that make up the population. We present an exact formula for the optimal proportion of migrants which is expressed in terms of the vital rates of migrant and non-migrant subtypes only. We show that this allocation strategy is both an evolutionary stable strategy (ESS) as well as a convergence stable strategy (CSS). To establish the former, we generalize the classical notion of an ESS because it is based on invasion exponents obtained from linearization arguments, which fail to capture the stabilizing effects of the nonlinear density dependence. These results clarify precisely when the notion of a "weak ESS", as proposed in Lundberg (2013) for a related model, is a genuine ESS. Secondly, we use an evolutionary game theory approach, and confirm, once again, that partial migration can be attributed to negative density dependence alone. In this context, the result holds even when density dependence is not attenuated. In this case, the optimal allocation strategy towards migrants is the same as the ESS stemming from the analysis based on the adaptive dynamics. The key feature of the population models considered here is that they are monotone dynamical systems, which enables a rather comprehensive mathematical analysis.

Shifting Thresholds: Rapid Evolution of Migratory Life Histories in Steelhead/Rainbow Trout, Oncorhynchus mykiss

Expression of phenotypic plasticity depends on reaction norms adapted to historic selective regimes; anthropogenic changes in these selection regimes necessitate contemporary evolution or declines in productivity and possibly extinction. Adaptation of conditional strategies following a change in the selection regime requires evolution of either the environmentally influenced cue (e.g., size-at-age) or the state (e.g., size threshold) at which an individual switches between alternative tactics. Using a population of steelhead (Oncorhynchus mykiss) introduced above a barrier waterfall in 1910, we evaluate how the conditional strategy to migrate evolves in response to selection against migration. We created 9 families and 917 offspring from 14 parents collected from the above- and below-barrier populations. After 1 year of common garden-rearing above-barrier offspring were 11% smaller and 32% lighter than below-barrier offspring. Using a novel analytical approach, we estimate that the mean size at which above-barrier fish switch between the resident and migrant tactic is 43% larger than below-barrier fish. As a result, above-barrier fish were 26% less likely to express the migratory tactic. Our results demonstrate how rapid and opposing changes in size-at-age and threshold size contribute to the contemporary evolution of a conditional strategy and indicate that migratory barriers may elicit rapid evolution toward the resident life history on timescales relevant for conservation and management of conditionally migratory species.

Effects of pyrethroid insecticides in urban runoff on Chinook salmon, steelhead trout, and their invertebrate prey

Pyrethroid insecticides can affect salmonids either indirectly through toxicity to their prey or directly by toxicity to the fish themselves. In support of a study on pyrethroid impacts to Chinook salmon and steelhead trout in the American River (Sacramento, California, USA), 96-h median effective concentration (EC50) and median lethal concentration (LC50) values for the pyrethroid bifenthrin were determined for taxa not traditionally used for toxicity testing but of interest as salmonid prey, including a chironomid, caddisflies, mayflies, and stoneflies. A laboratory was constructed on the banks of the American River to expose macroinvertebrates, Chinook salmon, and steelhead trout to flow-through river water containing urban runoff during storm events. Bifenthrin from urban runoff was found in river water following 5 rain events, reaching 14.6 ng/L. Mortality to the exposed salmonids was not observed, and sublethal effects were not seen in vitellogenin or sex steroid levels. Indirect effects via toxicity to salmonid prey are possible. Mortality to Hyalella azteca, a potential prey, was observed in every event tested, and peak bifenthrin concentrations were comparable to the 96-h EC50 of the caddisfly, Hydropsyche sp., the most important prey species on a biomass basis for American River Chinook salmon. The other invertebrates tested had EC50s exceeding bifenthrin concentrations seen in the American River, though could potentially be at risk at concentrations previously reported in smaller urban tributaries. Environ Toxicol Chem 2015;34:649-657. © 2014 SETAC.

Climate warming causes life-history evolution in a model for Atlantic cod (Gadus morhua)

Climate change influences the marine environment, with ocean warming being the foremost driving factor governing changes in the physiology and ecology of fish. At the individual level, increasing temperature influences bioenergetics and numerous physiological and life-history processes, which have consequences for the population level and beyond. We provide a state-dependent energy allocation model that predicts temperature-induced adaptations for life histories and behaviour for the North-East Arctic stock (NEA) of Atlantic cod (Gadus morhua) in response to climate warming. The key constraint is temperature-dependent respiratory physiology, and the model includes a number of trade-offs that reflect key physiological and ecological processes. Dynamic programming is used to find an evolutionarily optimal strategy of foraging and energy allocation that maximizes expected lifetime reproductive output given constraints from physiology and ecology. The optimal strategy is then simulated in a population, where survival, foraging behaviour, growth, maturation and reproduction emerge. Using current forcing, the model reproduces patterns of growth, size-at-age, maturation, gonad production and natural mortality for NEA cod. The predicted climate responses are positive for this stock; under a 2°C warming, the model predicted increased growth rates and a larger asymptotic size. Maturation age was unaffected, but gonad weight was predicted to more than double. Predictions for a wider range of temperatures, from 2 to 7°C, show that temperature responses were gradual; fish were predicted to grow faster and increase reproductive investment at higher temperatures. An emergent pattern of higher risk acceptance and increased foraging behaviour was also predicted. Our results provide important insight into the effects of climate warming on NEA cod by revealing the underlying mechanisms and drivers of change. We show how temperature-induced adaptations of behaviour and several life-history traits are not only mediated by physiology but also by trade-offs with survival, which has consequences for conservation physiology.

Stochastic dynamic programming illuminates the link between environment, physiology, and evolution

I describe how stochastic dynamic programming (SDP), a method for stochastic optimization that evolved from the work of Hamilton and Jacobi on variational problems, allows us to connect the physiological state of organisms, the environment in which they live, and how evolution by natural selection acts on trade-offs that all organisms face. I first derive the two canonical equations of SDP. These are valuable because although they apply to no system in particular, they share commonalities with many systems (as do frictionless springs). After that, I show how we used SDP in insect behavioral ecology. I describe the puzzles that needed to be solved, the SDP equations we used to solve the puzzles, and the experiments that we used to test the predictions of the models. I then briefly describe two other applications of SDP in biology: first, understanding the developmental pathways followed by steelhead trout in California and second skipped spawning by Norwegian cod. In both cases, modeling and empirical work were closely connected. I close with lessons learned and advice for the young mathematical biologists.

Using scale characteristics and water temperature to reconstruct growth rates of juvenile steelhead Oncorhynchus mykiss

Juvenile steelhead Oncorhynchus mykiss from a northern California Central Valley population were reared in a controlled laboratory experiment. Significantly different rates of growth were observed among fish reared under two ration treatments and three temperature treatments (8, 14 and 20°C). Wider circulus spacing and faster deposition was associated with faster growth. For the same growth rate, however, circulus spacing was two-fold wider and deposited 36% less frequently in the cold compared to the hot temperature treatment. In a multiple linear regression, median circulus spacing and water temperature accounted for 68% of the variation in observed O. mykiss growth. These results corroborate previous research on scale characteristics and growth, while providing novel evidence that highlights the importance of water temperature in these relationships. Thus, this study establishes the utility of using scale analysis as a relatively non-invasive method for inferring growth in salmonids.

Separating intrinsic and environmental contributions to growth and their population consequences

Among-individual heterogeneity in growth is a commonly observed phenomenon that has clear consequences for population and community dynamics yet has proved difficult to quantify in practice. In particular, observed among-individual variation in growth can be difficult to link to any given mechanism. Here, we develop a Bayesian state-space framework for modeling growth that bridges the complexity of bioenergetic models and the statistical simplicity of phenomenological growth models. The model allows for intrinsic individual variation in traits, a shared environment, process stochasticity, and measurement error. We apply the model to two populations of steelhead trout (Oncorhynchus mykiss) grown under common but temporally varying food conditions. Models allowing for individual variation match available data better than models that assume a single shared trait for all individuals. Estimated individual variation translated into a roughly twofold range in realized growth rates within populations. Comparisons between populations showed strong differences in trait means, trait variability, and responses to a shared environment. Together, individual- and population-level variation have substantial implications for variation in size and growth rates among and within populations. State-dependent life-history models predict that this variation can lead to differences in individual life-history expression, lifetime reproductive output, and population life-history diversity.

Contrasting effects of climate change in continental vs. oceanic environments on population persistence and microevolution of Atlantic salmon

Facing climate change (CC), species are prone to multiple modifications in their environment that can lead to extinction, migration or adaptation. Identifying the role and interplay of different potential stressors becomes a key question. Anadromous fishes will be exposed to both river and oceanic habitat changes. For Atlantic salmon, the river water temperature, river flow and oceanic growth conditions appear as three main stressing factors. They could act on population dynamics or as selective forces on life-history pathways. Using an individual-based demo-genetic model, we assessed the effects of these factors (1) to compare risks of extinction resulting from CC in river and ocean, and (2) to assess CC effects on life-history pathways including the evolution of underlying genetic control of phenotypic plasticity. We focused on Atlantic salmon populations from Southern Europe for a time horizon of three decades. We showed that CC in river alone should not lead to extinction of Southern European salmon populations. In contrast, the reduced oceanic growth appeared as a significant threat for population persistence. An increase in river flow amplitude increased the risk of local extinction in synergy with the oceanic effects, but river temperature rise reduced this risk. In terms of life-history modifications, the reduced oceanic growth increased the age of return of individuals through plastic and genetic responses. The river temperature rise increased the proportion of sexually mature parr, but the genetic evolution of the maturation threshold lowered the maturation rate of male parr. This was identified as a case of environmentally driven plastic response that masked an underlying evolutionary response of plasticity going in the opposite direction. We concluded that to counteract oceanic effects, river flow management represented the sole potential force to reduce the extinction probability of Atlantic salmon populations in Southern Europe, although this might not impede changes in migration life history.

Assessing adaptive phenotypic plasticity by means of conditional strategies from empirical data: the latent environmental threshold model

Conditional strategies are the most common form of discrete phenotypic plasticity. In a conditional strategy, the phenotype expressed by an organism is determined by the difference between an environmental cue and a threshold, both of which may vary among individuals. The environmental threshold model (ETM) has been proposed as a mean to understand the evolution of conditional strategies, but has been surprisingly seldom applied to empirical studies. A hindrance for the application of the ETM is that often, the proximate cue triggering the phenotypic expression and the individual threshold are not measurable, and can only be assessed using a related observable cue. We describe a new statistical model that can be applied in this common situation. The Latent ETM (LETM) allows for a measurement error in the phenotypic expression of the individual environmental cue and a purely genetically determined threshold. We show that coupling our model with quantitative genetic methods allows an evolutionary approach including an estimation of the heritability of conditional strategies. We evaluate the performance of the LETM with a simulation study and illustrate its utility by applying it to empirical data on the size-dependent smolting process for stream-dwelling Atlantic salmon juveniles.