ABSTRACT
Female-biased mortality has been repeatedly reported in Pacific salmon during their upriver migration in both field studies and laboratory holding experiments, especially in the presence of multiple environmental stressors, including thermal stress. Here, we used coho salmon (Oncorhynchus kisutch) to test whether females exposed to elevated water temperatures (18°C) (i) suppress circulating sex hormones (testosterone, 11-ketotestosterone and estradiol), owing to elevated cortisol levels, (ii) have higher activities of enzymes supporting anaerobic metabolism (e.g. lactate dehydrogenase, LDH), (iii) have lower activities of enzymes driving oxidative metabolism (e.g. citrate synthase, CS) in skeletal and cardiac muscle, and (iv) have more oxidative stress damage and reduced capacity for antioxidant defense [lower catalase (CAT) activity]. We found no evidence that a higher susceptibility to oxidative stress contributes to female-biased mortality at warm temperatures. We did, however, find that females had significantly lower cardiac LDH and that 18°C significantly reduced plasma levels of testosterone and estradiol, especially in females. We also found that relative gonad size was significantly lower in the 18°C treatment regardless of sex, whereas relative liver size was significantly lower in females held at 18°C. Further, relative spleen size was significantly elevated in the 18°C treatments across both sexes, with larger warm-induced increases in females. Our results suggest that males may better tolerate bouts of cardiac hypoxia at high temperature, and that thermal stress may also disrupt testosterone- and estradiol-mediated protein catabolism, and the immune response (larger spleens), in migratory female salmon.
INTRODUCTION
Female-biased mortality represents an especially important conservation concern because female numbers typically determine the reproductive potential of a population, and subsequent male biases in operational sex ratios can be disproportionately costly for females (Hardy, 2002; Rankin and Kokko, 2006; Quinn, 1999; Garner et al., 2009). Over the past decade, female-biased mortality has been increasingly reported during upriver spawning migrations of Pacific salmon (e.g. Crossin et al., 2008; Jeffries et al., 2012; Martins et al., 2012; Keefer et al., 2008; Minke-Martin et al., 2018). This trend is especially apparent in the presence of multiple environmental stressors, including temperature (e.g. Jeffries et al., 2014). Although female-biased mortality has been demonstrated in both field (Burnett et al., 2014; Keefer et al., 2008; Martins et al., 2012; Minke-Martin et al., 2018; Roscoe et al., 2011) and laboratory holding studies (Jeffries et al., 2012; Raby et al., 2016), the underlying mechanisms that enhance female vulnerability during this once-in-a-lifetime migration are unknown. Understanding the physiological basis for female-biased mortality represents a crucial goal to help manage these ecologically and economically important species in a changing climate.
Thermal stress arguably represents the most pervasive factor influencing how populations will respond to the effects of climate change (e.g. Pörtner and Farrell, 2008). The current trend has seen a 1–2°C increase in peak summer water temperatures in the Fraser River basin (British Columbia, Canada), with future predictions estimating a 3-fold increase in the number of days above 19°C (a critical salmon threshold) over the next 50–100 years (Ferrari et al., 2007; Hague et al., 2011). Adult Pacific salmon are especially vulnerable to rising summer temperatures because this is when many populations migrate upstream, already encountering river temperatures near their upper thermal limits (Eliason and Farrell, 2016). Increasing water temperatures are known to compromise cardiorespiratory performance (Eliason et al., 2011; Eliason et al., 2013), upregulate the hypothalamic–pituitary–interrenal (HPI) axis (Strange et al., 1977; Raby et al., 2016) and promote the production of reactive oxygen species (ROS; Nakano et al., 2014). The present study aimed to determine the physiological differences that make female Pacific salmon more vulnerable to thermal stress than males.
There is some evidence that sex-specific differences in aerobic capacity may underlie female-biased mortality with rising temperatures. In pink salmon (Oncorhynchus gorbuscha), for instance, aerobic scope (the absolute difference between standard metabolic rate and maximum metabolic rate; AS) is lower in females than in males (Clark et al., 2011). In sockeye salmon (O. nerka), females rely more heavily on anaerobic metabolism during dam passage than males (Burnett et al., 2014), and have higher levels of cardiac lactate following handling stress (Eliason et al., 2020). In another salmonid, the rainbow trout (O. mykiss), females have significantly lower scope for coronary blood flow at high temperatures (Ekström et al., 2017). These combined data suggest that females may be more vulnerable to oxygen supply limitations during thermal stress, especially in the heart. Although upregulation of anaerobic metabolism may help offset oxygen limitations, this compensatory response is less efficient in terms of energy conversion and may carry other costs. Higher rates of glycolytically fuelled burst swimming, for instance, have been associated with upriver migration mortality (Burnett et al., 2014). We therefore hypothesized that skeletal and cardiac muscle of female Pacific salmon are more primed for anaerobic metabolism than their male counterparts.
A second hypothesis is that differences in the neuroendocrine stress response underlie female-biased mortality at high temperatures. Plasma cortisol levels typically peak during upriver migration (Hruska et al., 2007; Carruth et al., 2002), and female Pacific salmon are known to have far higher levels of circulating cortisol than their male counterparts (e.g. Schmidt and Idler, 1962; Hinch et al., 2006). Females are also known to experience elevated and more prolonged cortisol spikes in response to stress (Hruska et al., 2010; Cook et al., 2011; Donaldson et al., 2014; Raby et al., 2016). Increases in cortisol energetically prime fish to respond to acute stressors by increasing blood glucose, liver glycogen and branchial blood flow (Barton and Iwama, 1991; Mommsen et al., 1999). However, chronic elevation of plasma cortisol above baseline levels as a result of stress can disrupt intermediary metabolism (see Van Der Boon et al., 1991) and reproductive endocrine signalling (Carragher et al., 1989; Pankhurst and Van Der Kraak, 2000; McConnachie et al., 2012). This means that a cortisol-driven decrease in testosterone and estradiol can then potentially delay sexual maturation rates (Truscott et al., 1986; Jeffries et al., 2012), and has also been linked to river mortality in Pacific salmon (Cooke et al., 2006). We were therefore interested in how high temperature stress differentially affects glucose metabolism and sex hormone signalling in male and female Pacific salmon.
A third hypothesis is that females may be more susceptible to oxidative damage as temperatures rise. At a molecular level, ROS damage DNA, lipids and proteins (see Beckman and Ames, 1998). Importantly, it has been hypothesized that ROS production indirectly increases costs of reproduction and exercise (Beckman and Ames, 1998; Costantini, 2008). Cessation of feeding during the spawning migration means Pacific salmon only have finite antioxidant stores to protect against ROS (Wilson et al., 2014), and increased metabolic rates as a function of higher river temperatures are also likely to augment ROS production (Bagnyukova et al., 2007). There is some evidence, however, that females may be more vulnerable to ROS than males. High levels of testosterone have been linked to increased production of ROS (see Metcalfe and Alonso-Alvarez, 2010), and high levels of estradiol can compromise antioxidant activity (Thilagam et al., 2010). Plasma testosterone levels are relatively low and stable throughout the upriver migration in male Pacific salmon, where estradiol is virtually absent (Sower and Schreck, 1982; Truscott et al., 1986; Hruska et al., 2007). In females, however, both testosterone and estradiol naturally increase throughout the upriver migration, reaching levels between 2- and 3-fold higher than during their initial entry into the river (Sower and Schreck, 1982; Truscott et al., 1986; Hruska et al., 2007). We might thereby expect increased ROS damage in females, as climbing levels of testosterone and estradiol enhance ROS production and reduce rates of ROS clearance, respectively.
The aim of our study was to explore these three hypotheses that could potentially explain the elevated female mortality phenomenon in Pacific salmon. Specifically, we sought to determine whether differences in anaerobic versus oxidative metabolism, sex hormone signalling and oxidative stress might contribute to female-biased mortality in Pacific salmon exposed to thermal stress. We acclimated mature migratory coho salmon (O. kisutch) to two separate thermal regimes (9°C and 18°C). We then measured differences in metabolic enzyme activities (lactate dehydrogenase, LDH; citrate synthase, CS) in skeletal and cardiac muscles, blood chemistry (lactate, glucose, testosterone, estradiol and electrolytes), ROS damage [thiobarbituric acid reactive substances (TBARS)] and antioxidant enzyme activity (catalase, CAT) between male and females from each temperature treatment. We predicted that relative to males and females at 9°C and males at 18°C, females exposed to 18°C would have: (i) higher LDH but lower CS activities (i.e. females would be more primed for anaerobic exercise), especially in the heart, (ii) higher plasma glucose and lower levels of circulating sex hormones (testosterone and estradiol) and (iii) more ROS damage (higher TBARS concentrations) and reduced capacity for antioxidant defense (lower CAT activity).
MATERIALS AND METHODS
Animal collection and holding
Adult coho salmon [Oncorhynchus kisutch (Walbaum 1792)] were collected from the Chilliwack River Hatchery (British Columbia, Canada; hatchery water temperature 9°C during collection) and transported 23.7 km in a 1250 l tank (8.2–10.4°C; >90% air saturation) to the Fisheries and Oceans Canada Cultus Lake Research Laboratory (British Columbia, Canada). Fish were held in flow-through UV-sterilized and sand-filtered freshwater under natural photoperiod in 8000 l outdoor holding tanks (9°C; >90% air saturation; n≤27 fish per tank, mixed sex) for at least 36 h before experimentation to ensure recovery from handling stress. Fish were then moved to either 9°C or 18°C treatment tanks (1.8 m diameter; 3–6 fish per tank), where water either remained at 9°C, or was increased at a rate of 2–4°C h−1 until 18°C was reached. This represents an ecologically relevant rate of warming in Pacific salmon, which perform rapid vertical migrations both in the ocean and in lakes, during which they can encounter a 5–10°C temperature swing in less than an hour (Eliason and Farrell, 2016). Fish were held at 18°C for 4.77±0.34 days (mean±s.d.; range: 4.51–5.75 days) before they were sampled compared with their 9°C counterparts. We chose 9°C to approximate current river temperatures, whereas 18°C was used to represent a future climate change scenario for this population. It should be noted, however, that some populations across the southern range of the species already experience temperatures as high as 18°C (Welsh et al., 2001; Hayes et al., 2011).
Following these acclimation periods, fish were removed from their tanks and rapidly euthanized by cerebral concussion. A blood sample was immediately collected from the caudal vein and stored in an ice slurry before being centrifuged at 7000 g for 5 min. Plasma was stored at −80°C until analysis. A skeletal muscle sample (<0.5 cm thick), comprising both red and white muscle fibres, was removed from halfway between the dorsal and anal fin, and the ventricle was removed and bisected. The muscle and ventricle samples were blotted dry of blood, rapidly freeze-clamped in liquid nitrogen and stored at −80°C until analysis. Morphometric data (e.g. body mass, length and organ masses) were collected from the carcasses. In general, we acclimated N=12 males and N=12 females to 9°C, and N=12 males and N=12 females to 18°C. However, actual sample sizes vary between the different assays and are further noted in Tables 1 and 2. These sample sizes are known to provide sufficient statistical power to detect effect sizes expected from previous work in Pacific salmon or closely related species. All experimental protocols were approved by the Animal Care Committee at the University of British Columbia in accordance with the Canadian Council on Animal Care (protocol no. A17-0160).
Results and statistical (ANOVA) parameters for blood chemistry, sex hormones, relative organ sizes and thiobarbituric acid reactive substances (TBARS) of wild adult coho salmon acclimated to two test temperatures

Blood chemistry
Plasma metabolites were analyzed according to established methods (Farrell et al., 2001; Richards and Hultin, 2002). Briefly, sodium and potassium were determined using an XP Five-channel Flame Photometer (BWB Technologies, UK), chloride was measured with a Chlorochek Digital Chloridometer (EliTech Group, France), and osmolality was determined using a 3320 Freezing Point Osmometer (Advanced Instruments, USA). Lactate and glucose were measured with a 2300 stat plus glucose and l-lactate analyzer (YSI, USA). All samples were run in duplicate.
Sex hormones
Testosterone, 11-ketotestosterone and 17B-estradiol were extracted using ethyl ether, and quantified using commercially available ELISA kits (Testosterone and 17B-estradiol, Neogen, USA; 11-ketotestosterone, Cayman Chemical, USA) according to manufacturer instructions. All plates were read on a FLUOstar Omega multi-mode microplate reader (BMG Labtech, USA). Measurements for testosterone and 17B-estradiol were made in duplicate, whereas those for 11-ketotestosterone were made in duplicate across at least two dilutions to test for interference. Assays for 11-ketotestosterone were performed approximately 2 years after tissue collection. Some samples therefore were rejected owing to high intra- and inter-assay variation in concentrations.
Enzyme assays
Enzyme activities for LDH (anaerobic), CS (oxidative) and CAT (oxidative repair) were measured in homogenates of ventricle, red muscle and white muscle tissues to determine thermal performance across 9, 14, 18, 22 and 25°C test temperatures according to previously described methods (LDH, Martínez et al., 2006; CS, Moyes et al., 1997; CAT, Ghanizadeh Kazerouni et al., 2016). In short, tissues were homogenized in homogenization buffer (LDH and CS: 0.1% Triton, 50 mmol l−1 Hepes, 1 mmol l−1 EDTA, pH 7.4; CAT: 100 mmol l−1 K phosphate buffer, 100 mmol l−1 KCl, 1 mmol l−1 EDTA, pH 7.4) on ice. All assays were run in triplicate on a SpectraMax iD3 Multi-Mode Microplate Reader (Molecular Devices, USA) using a wavelength of 340 nm to measure the disappearance of NADH for LDH activity, a wavelength of 412 nm to measure the production of 2-nitro-5-thiobenzoic acid as a proxy of CS activity, or a wavelength of 240 nm to measure the disappearance of H2O2 for CAT activity. Absorbance readings were automatically normalized to 1 cm via the integrated PathCheck sensor technology (Molecular Devices, USA). Extinction coefficients were 6.22, 14.1 and 40 l mmol−1 cm−1 for LDH, CS and CAT, respectively.
Oxidative damage
TBARS were quantified separately for tissue from the ventricle, red skeletal muscle and white skeletal muscle to identify levels of lipid peroxidation using a commercially available colorimetric assay kit (Cayman Chemical). Briefly, tissues were homogenized in RIPA buffer and treated with thiobarbituric acid under high temperature (100°C) and acidity. TBARS were then quantified by comparing optical densities at 535 nm to simultaneously determined standard curves for malondialdehyde. All samples were run in duplicate.
Statistical analysis
All statistical analyses were done using the free and open software JAMOVI (Version 0.9; GAMLj module; https://www.jamovi.org). Categorical variables [blood chemistry (lactate, glucose, Cl−, osmolality, Na+, K+), sex hormones (17B-estradiol, testosterone, 11-ketotestosterone), organ masses (gonads, liver, ventricle and spleen)] and TBARS were analyzed with general linear models and ANOVA. Acclimation temperature, sex and their interaction term were used as predictor variables in all models. In models examining organ masses, we included body mass as a covariate. Kinetic enzyme assays were analyzed with mixed models with effect terms: sex, acclimation temperature, test temperature and their interaction terms (sex×acclimation temperature, sex×test temperature, acclimation temperature×test temperature, and a three-way interaction sex×acclimation temperature×test temperature). Fish ID was used as a cluster variable to account for repeated measures across the different test temperatures. Metrics were investigated for normality using quantile-quantile plots. Non-normal data were transformed (log10, square root, Box–Cox) prior to statistical analysis. Heteroskedasticity was subsequently assessed via Levene's test. All data are displayed with untransformed values. Significance levels were set at 0.05, and data are presented as means±s.e.m. in figures and as means±s.d. in tables. Thermal performance curves for enzyme activity are presented with either quadratic or linear curves, depending on best fit (Akaike's information criterion).
RESULTS
Relative organ size
Females had a significantly larger GSI than males, and both sexes in the 18°C treatment had significantly smaller gonads than at 9°C (9°C: males 5.4±1.3%, females 19.8±2.5%; 18°C: males 4.9±1.1%, females: 17.5±0.8%; Table 1). There was no significant effect of acclimation temperature on mean ventricle size relative to body mass (relative ventricular mass; RVM). However, there was a significant effect of sex, where RVM was approximately 12% higher in males (Table 1). There was a significant interaction between acclimation temperature and sex on relative spleen size (splenosomatic index, SSI; Table 1). Specifically, mean female SSI was approximately half that of males at 9°C but increased more drastically with temperature so that it was nearly 20% larger at 18°C. Specifically, SSI was elevated in the 18°C treatment for both sexes (Table 1). Also, a significant interaction existed between sex and acclimation temperature on relative liver size (hepatosomatic index, HSI), with mean HSI being lower in females held at 18°C but no different in males (Table 1).
Blood chemistry
Plasma lactate levels were significantly higher in the 18°C treatments. This effect was independent of sex, as plasma lactate did not significantly differ between males and females at either temperature (Table 1). There were no significant effects of sex or acclimation temperature on plasma osmolality or the plasma concentrations of glucose, Na+, K+ or Cl− (Table 1).
Enzyme assays
In the ventricle, there was a main effect of sex and test temperature on LDH activity (Fig. 1A, Table 2). Males had higher ventricular rates of LDH activity. LDH activity increased with test temperature in both sexes, with a marginally significant interaction between sex and test temperature, where enzyme activity appears to increase more quickly with test temperature in males than in females (Table 2). In contrast, neither sex nor acclimation temperature had significant effects on maximal LDH activity in white or red skeletal muscle, although LDH activity similarly increased with test temperature in both tissues (Fig. 1B,C, Table 2). There were no significant effects of sex or acclimation temperature on CS activity in ventricular, white muscle or red muscle tissues (Fig. 2, Table 2). For the antioxidant enzyme activities, acclimation temperature and test temperature had significant effects on CAT activity in red muscle (Fig. 3C, Table 2). CAT activity was significantly higher in fish acclimated to 9°C, regardless of sex, and increased with test temperature. A significant interaction existed between acclimation temperature and test temperature in CAT activity in white muscle, where CAT activity was significantly higher in 9°C-acclimated fish at the higher test temperatures (Fig. 3B, Table 2). Also, test temperature significantly affected CAT activity in the ventricle, and there was a trend (P=0.097) for ventricular CAT activity to be higher in males (Fig. 3A, Table 2).
Lactate dehydrogenase (LDH) activity in skeletal and cardiac muscle in coho salmon. (A) Cardiac muscle, (B) white muscle and (C) red muscle. A statistically significant (P<0.05) main effect of sex is noted in A. Test temperature represented a significant main effect in all panels. Each sample represents a biological replicate and all samples were run in triplicate. Data are means±s.e.m.
Lactate dehydrogenase (LDH) activity in skeletal and cardiac muscle in coho salmon. (A) Cardiac muscle, (B) white muscle and (C) red muscle. A statistically significant (P<0.05) main effect of sex is noted in A. Test temperature represented a significant main effect in all panels. Each sample represents a biological replicate and all samples were run in triplicate. Data are means±s.e.m.
Citrate synthase (CS) activity in skeletal and cardiac muscle in coho salmon. (A) Cardiac muscle, (B) white muscle and (C) red muscle. A statistically significant (P<0.05) main effect of acclimation temperature (AT) is noted in B. Test temperature represented a significant main effect in all panels. Each sample represents a biological replicate and all samples were run in triplicate. Data are means±s.e.m.
Citrate synthase (CS) activity in skeletal and cardiac muscle in coho salmon. (A) Cardiac muscle, (B) white muscle and (C) red muscle. A statistically significant (P<0.05) main effect of acclimation temperature (AT) is noted in B. Test temperature represented a significant main effect in all panels. Each sample represents a biological replicate and all samples were run in triplicate. Data are means±s.e.m.
Catalase (CAT) activity in skeletal and cardiac muscle in coho salmon. (A) Cardiac muscle, (B) white muscle and (C) red muscle. Statistically significant (P<0.05) effects are represented by AT (main effect of acclimation temperature) and AT×TT (interaction between acclimation temperature and test temperature). Test temperature represented a significant main effect in all cardiac (A) and red muscle (C). Each sample represents a biological replicate and all samples were run in triplicate. Data are means±s.e.m.
Catalase (CAT) activity in skeletal and cardiac muscle in coho salmon. (A) Cardiac muscle, (B) white muscle and (C) red muscle. Statistically significant (P<0.05) effects are represented by AT (main effect of acclimation temperature) and AT×TT (interaction between acclimation temperature and test temperature). Test temperature represented a significant main effect in all cardiac (A) and red muscle (C). Each sample represents a biological replicate and all samples were run in triplicate. Data are means±s.e.m.
Sex hormones
A significant interaction existed between sex and acclimation temperature on plasma testosterone levels (Table 1). Although mean testosterone levels were significantly lower in the 18°C treatment in both sexes, the difference was much more pronounced in females (Table 1). Males had significantly higher levels of 11-ketotestosterone than females, and warmer temperatures significantly decreased 11-ketotestosterone concentrations in both sexes (Table 1). A significant interaction existed between sex and acclimation temperature on plasma estradiol levels (Table 1), with estradiol levels decreasing by nearly 75% in females in the 18°C treatment, but not in males, where the estradiol level was negligible (Table 1).
Oxidative damage
There were no significant effects of sex or acclimation temperature on the concentration of TBARS in red muscle, white muscle or in the ventricle (Table 1). However, the effect of sex on TBARS concentrations approached statistical significance in red muscle (P=0.087), where males tended to have higher levels in red muscle than females (Table 1).
DISCUSSION
Here, we considered potential physiological mechanisms that could contribute to female-biased mortality in response to thermal stress during the spawning migration in Pacific salmon. Although our results render it difficult to pinpoint a specific mechanism, we revealed that (i) sex-specific differences in anaerobic metabolism in the heart, and (ii) changes in sex hormone signalling might play roles in the elevated female mortality rates observed at high temperatures.
Anaerobic metabolism
We found no evidence that females rely more on anaerobic metabolism at rest during high temperature exposure than males. At a tissue level, we predicted that females would be more vulnerable to cardiac oxygen limitations at high temperatures and would therefore upregulate anaerobic pathways. Contrary to this hypothesis, however, ventricular LDH activity was significantly higher in males than in females. Because LDH is needed for a tissue to cope with oxygen limitations, the higher LDH activities in male hearts suggests a greater capacity to cope with hypoxia. This could be due to genetic differences that result in higher constitutive LDH in males or a plastic response to prior hypoxia stress. Increased LDH activities may also indicate that males may have a greater capacity to metabolize lactate (Chippari-Gomes et al., 2005; Crans et al., 2015; Borowiec et al., 2016), which appears to be a preferred extracellular substrate to glucose in salmonids (Lanctin et al., 1980; Farrell et al., 1988). The nearly significant interaction between sex and test temperatures (Table 2) indicates that LDH activities in male ventricles increase more quickly as a function of ambient temperature than in females. In fact, ventricular LDH activities in males were more than 10% higher than for females at the 18°C test temperatures. This represents a departure from previous work in sexually immature hatchery-reared rainbow trout, where female ventricles had higher LDH activities than males (Battiprolu et al., 2007). Moreover, male Pacific salmon are known to have a larger RVM than females (Clark et al., 2008; Sandblom et al., 2009). In the present study, mean ventricular mass in males was approximately 16% bigger than in females (male RVM, 0.199; female RVM, 0.171; Table 1). Earlier work suggested that a cardiac oxygen limitation may contribute to higher mortality in female Pacific salmon during stress (Eliason et al., 2020). Ventricles of female rainbow trout are also more hypoxia tolerant (Battiprolu et al., 2007; Battiprolu and Rodnick, 2014), suggesting that they may experience bouts of hypoxia more routinely – though this idea has not been tested. Our work indicates that a reduced reliance on anaerobic metabolism combined with smaller ventricles in female hearts may undermine cardiac performance as temperatures rise. Collectively, these studies suggest that a cardiac limitation might contribute to heightened female mortality.
Stress response
Migratory female salmon are known to have elevated plasma cortisol (e.g. Schmidt and Idler, 1962; Kubokawa et al., 2001) with higher and more prolonged cortisol spikes in response to stress (e.g. Raby et al., 2016; Hruska et al., 2010; Cook et al., 2011; Donaldson et al., 2014; Eliason et al., 2020). We therefore predicted that downstream pathways (i.e. sex hormone signalling and carbohydrate metabolism) would be disrupted in females exposed to thermal stress (i.e. the 18°C treatment). We did find evidence that female-biased mortality in Pacific salmon may stem from the disruption of sex hormone signalling. As predicted, testosterone levels were markedly higher in females, and the effect of temperature was much more pronounced. Although testosterone drives gonadal development and secondary sex characteristics in male Pacific salmon, it also plays crucial roles in migratory females, including the regulation of oocyte development and final maturation (Truscott et al., 1986). Reductions in testosterone levels may thereby delay rates of maturation in female salmon at high temperatures, which potentially explains the reduced rates of final maturation observed in female pink and sockeye salmon held at 19°C (Jeffries et al., 2012). Supporting this suggestion, we found reduced GSIs in both sexes of the 18°C treatment relative to their 9°C counterparts. We also found reduced mean liver size in females, but relatively no effect in males. Because estradiol signals liver synthesis of vitellogenin and eggshell proteins (Korsgaard et al., 1986; Berg et al., 2004; Larsen et al., 1992), thereby increasing liver size, the relatively low HSI of females at 18°C may also represent a delay in maturation, or an overall reduction in egg mass, as a consequence of reduced estradiol levels. Alternatively, a lower HSI is also consistent with reduced energy reserves (i.e. glycogen; Arndt et al., 1996), potentially indicating substrate limitations in female salmon exposed to high temperatures.
Although the mechanisms remain unknown, lower levels of circulating testosterone have also been linked to upper river mortality in sockeye salmon (Cooke et al., 2006). However, the large decrease in 11-ketotestosterone suggests that males may face their own sex-specific challenges with warming temperatures. 11-Ketotestosterone represents a primary androgen in male salmon (Idler et al., 1961), where it promotes sexual maturation and upstream migratory behaviours (Munakata et al., 2001). This means that maturation may also be delayed in males exposed to high migratory temperatures. Although circulating levels of estradiol were virtually non-detectable in males at either temperature, they plummeted by approximately 75% in females held in the 18°C treatment when compared with the 9°C treatment. Again, this suggests that low estradiol may contribute to delays in female maturation as a function of temperature. These warm temperature effects were seen for other sex hormones, with females appearing to suffer more drastic reductions in sex hormones than do males at warm temperatures. For example, plasma testosterone concentrations, which were more than 3-fold higher in females at 9°C, fell by roughly 80% in females exposed to 18°C and by 70% in males exposed to the same temperature (Table 1). Also, the 75% decrease in estradiol concentration in females during warming is comparable to the 50% decrease 11-ketotestosterone in males during warming. Although it is difficult to reconcile how concentration differences in these three hormones impact each sex comparatively, it is clear that females suffer greater absolute and factorial reductions in sex hormones at high temperatures. Both testosterone and estradiol levels have also been correlated with faster and more successful migrations (Sato et al., 1997; Crossin et al., 2008; Cooke et al., 2008), although mechanisms underlying these effects are not known. However, both testosterone and estradiol increase rates of protein utilization, which likely underlie the repartitioning of white muscle proteins to the gonads (Cleveland and Weber, 2011), and facilitate white muscle as an energy source once lipid reserves become depleted (Ando, 1986). A caveat may be that female Pacific salmon suffer substrate limitations during migration, as high temperatures reduce plasma testosterone and estradiol, potentially compromising these important signalling pathways for protein catabolism.
Contrary to our expectations, there were no significant differences in plasma lactate, glucose or electrolytes (i.e. Na+, K+ or Cl−) between males and females at either acclimation temperature. This reflects a divergence from previous work in pink salmon (O. gorbuscha) and sockeye salmon, where plasma osmolality and Cl− were elevated following exposure to 19°C, and both plasma glucose and lactate were elevated in females (Jeffries et al., 2012). Although a substantial body of work in rainbow trout suggests that cortisol mediates glucose and lactate metabolism (e.g. Eros and Milligan, 1996; Milligan, 1996, 1997; Milligan et al., 2000), the link between circulating cortisol and plasma lactate varies greatly across teleost species (see Mommsen et al., 1999), and can even differ among salmonids (Barton, 2000). Previous work in female pink salmon supplemented with exogenous cortisol during their upriver migration also found no evidence elevated cortisol influences blood glucose levels (McConnachie et al., 2012).
Cortisol is also known to disrupt the immune response (Maule et al., 1987; 1989; Maule and Schreck, 1990), and is positively associated with disease-induced mortality (Pickering and Pottinger, 1989). This means that females may be at greater risk of infection than males. Although we did not set out to test this hypothesis, the higher relative spleen size in males relative to females may represent enhanced immunocompetence. For example, SSI is known to positively correlate with resistance to infections in rainbow trout (Hadidi et al., 2008). Further, the increase in SSI in both sexes in the warm treatment, where aquatic infections likely intensify (Snieszko, 1974; McCullough, 1999; Materna, 2001; Crossin et al., 2008), may signify a mounting immune response. Together, this may mean that the higher cortisol levels and relatively smaller spleens of female Pacific salmon reflect a heightened susceptibility to infection at high temperature. In fact, recent work shows that higher intensities of infection at elevated temperatures at least partly account for higher mortality rates in female coho salmon following handling stress in warm water (Teffer et al., 2019). Differences in spleen size could also reflect differences in the ability to mobilize sequestered red blood cells (Gallaugher et al., 1992), where males have the greater ability.
Oxidative stress
At odds with our hypothesis, we found little evidence that differences in oxidative damage explain female-biased mortality in migratory Pacific salmon despite the fact that oxidative stress markers are regarded as sensitive measures (Birnie-Gauvin et al., 2017). Females did not have elevated evidence of ROS damage (i.e. TBARS), nor did they appear to be investing in antioxidant defences (i.e. CAT activity) more than males in any of the tissues we examined. Red muscle had the highest level of TBARS, presumably owing to its high aerobic capacity (Wilhelm Filho, 2007; Mortelette et al., 2010). Contrary to our prediction, oxidative damage was not more pronounced with acclimation temperature in any of the tissues studied here. However, increased rates of CAT activity in red muscle of fish acclimated to 9°C may indicate a compensatory response to increased ROS challenge at cold temperatures. And it is also possible that reductions in plasma testosterone and estradiol at 18°C (see above) contribute to decreased oxidative defences (Thilagam et al., 2010). Previous work in pink salmon highlights the importance of tissue type and stage of migration when looking for effects of oxidative stress, where oxidative stress in the heart was found to increase over the course of migration, while red and white muscle had no significant changes (Wilson et al., 2014). It should be noted that fish used for our study were collected at their spawning ground (i.e. Chilliwack River hatchery), as opposed to along their migratory route. This means that our sampling may have been biased for the most resilient females by not targeting females at a potentially more vulnerable state earlier in migration.
Conclusions
Our results have helped exclude one hypothesis previously used to explain female-biased mortality at high temperatures. We found no evidence that a higher susceptibility to ROS damage contributes to female-biased mortality at warm temperatures. Our results, however, do bolster other hypotheses. Potential avenues for future work include testing whether high temperatures (i) disproportionately undermine the metabolic capacity of the heart, (ii) disrupt testosterone- and estradiol-mediated protein catabolism and (iii) affect immune responses (larger increases in spleen size) in migratory female salmon. Investigation into these processes, especially earlier along the migratory route, will undoubtedly help to identify whether they contribute to female-biased mortality at high temperatures.
Acknowledgements
We thank Andrew Lotto and John Little for field support, Bryan Smith and the rest of the staff at the Cultus Lake Research Lab, and Jeremy Mothus and the rest of the staff at Chilliwack Hatchery.
Footnotes
Funding
This work was supported by a University of California Santa Barbara Faculty Research Award to E.J.E. A.P.F. and S.J.C. are funded by the Natural Sciences and Engineering Research Council of Canada (NSERC) and the Canada Research Chair program. Funding was also provided by NSERC Discovery and Strategic Grants to S.G.H., S.J.C. and A.P.F. Funded in part by DFO's Aquatic Climate Change Adaptation Services Program to D.A.P.
Author contributions
Conceptualization: A.P.F., S.J.C., D.A.P., S.G.H., E.J.E.; Methodology: A.G.L., E.H., K.K., T.D., T.S.P., B.H., E.J.E.; Validation: E.J.E.; Formal analysis: A.G.L., K.K., E.J.E.; Investigation: A.G.L., E.H., K.K., T.D., T.S.P., B.H., E.J.E.; Resources: J.N.P., A.P.F., S.J.C., D.A.P., S.G.H., E.J.E.; Data curation: E.H.; Writing - original draft: A.G.L., E.J.E.; Writing - review & editing: A.G.L., E.H., K.K., T.D., T.S.P., B.H., J.N.P., A.P.F., S.J.C., D.A.P., S.G.H.; Supervision: A.G.L., J.N.P., A.P.F., S.J.C., D.A.P., S.G.H., E.J.E.; Project administration: E.J.E.; Funding acquisition: J.N.P., A.P.F., S.J.C., D.A.P., S.G.H., E.J.E.
References
Competing interests
The authors declare no competing or financial interests.