Identification and Optimization of Novel Inhibitors of the Polyketide Synthase 13 Thioesterase Domain with Antitubercular Activity.

There is an urgent need for new tuberculosis (TB) treatments, with novel modes of action, to reduce the incidence/mortality of TB and to combat resistance to current treatments. Through both chemical and genetic methodologies, polyketide synthase 13 (Pks13) has been validated as essential for mycobacterial survival and as an attractive target for Mycobacterium tuberculosis growth inhibitors. A benzofuran series of inhibitors that targeted the Pks13 thioesterase domain, failed to progress to preclinical development due to concerns over cardiotoxicity. Herein, we report the identification of a novel oxadiazole series of Pks13 inhibitors, derived from a high-throughput screening hit and structure-guided optimization. This new series binds in the Pks13 thioesterase domain, with a distinct binding mode compared to the benzofuran series. Through iterative rounds of design, assisted by structural information, lead compounds were identified with improved antitubercular potencies (MIC < 1 µM) and in vitro ADMET profiles.

Lysyl-tRNA synthetase, a target for urgently needed M. tuberculosis drugs.

Tuberculosis is a major global cause of both mortality and financial burden mainly in low and middle-income countries. Given the significant and ongoing rise of drug-resistant strains of Mycobacterium tuberculosis within the clinical setting, there is an urgent need for the development of new, safe and effective treatments. Here the development of a drug-like series based on a fused dihydropyrrolidino-pyrimidine scaffold is described. The series has been developed against M. tuberculosis lysyl-tRNA synthetase (LysRS) and cellular studies support this mechanism of action. DDD02049209, the lead compound, is efficacious in mouse models of acute and chronic tuberculosis and has suitable physicochemical, pharmacokinetic properties and an in vitro safety profile that supports further development. Importantly, preliminary analysis using clinical resistant strains shows no pre-existing clinical resistance towards this scaffold.

Toolkit of Approaches To Support Target-Focused Drug Discovery for Plasmodium falciparum Lysyl tRNA Synthetase.

There is a pressing need for new medicines to prevent and treat malaria. Most antimalarial drug discovery is reliant upon phenotypic screening. However, with the development of improved target validation strategies, target-focused approaches are now being utilized. Here, we describe the development of a toolkit to support the therapeutic exploitation of a promising target, lysyl tRNA synthetase (PfKRS). The toolkit includes resistant mutants to probe resistance mechanisms and on-target engagement for specific chemotypes; a hybrid KRS protein capable of producing crystals suitable for ligand soaking, thus providing high-resolution structural information to guide compound optimization; chemical probes to facilitate pulldown studies aimed at revealing the full range of specifically interacting proteins and thermal proteome profiling (TPP); as well as streamlined isothermal TPP methods to provide unbiased confirmation of on-target engagement within a biologically relevant milieu. This combination of tools and methodologies acts as a template for the development of future target-enabling packages.

Neutral and adaptive drivers of genomic change in introduced brook trout (Salvelinus fontinalis) populations revealed by pooled sequencing.

Understanding the drivers of successful species invasions is important for conserving native biodiversity and for mitigating the economic impacts of introduced species. However, whole-genome resolution investigations of the underlying contributions of neutral and adaptive genetic variation in successful introductions are rare. Increased propagule pressure should result in greater neutral genetic variation, while environmental differences should elicit selective pressures on introduced populations, leading to adaptive differentiation. We investigated neutral and adaptive variation among nine introduced brook trout (Salvelinus fontinalis) populations using whole-genome pooled sequencing. The populations inhabit isolated alpine lakes in western Canada and descend from a common source, with an average of ~19 (range of 7-41) generations since introduction. We found some evidence of bottlenecks without recovery, no strong evidence of purifying selection, and little support that varying propagule pressure or differences in local environments shaped observed neutral genetic variation differences. Putative adaptive loci analysis revealed nonconvergent patterns of adaptive differentiation among lakes with minimal putatively adaptive loci (0.001%-0.15%) that did not correspond with tested environmental variables. Our results suggest that (i) introduction success is not always strongly influenced by genetic load; (ii) observed differentiation among introduced populations can be idiosyncratic, population-specific, or stochastic; and (iii) conservatively, in some introduced species, colonization barriers may be overcome by support through one aspect of propagule pressure or benign environmental conditions.

Assessing conservation status with extensive but low-resolution data: Application of frequentist and Bayesian models to endangered Athabasca River rainbow trout.

Use of extensive but low-resolution abundance data is common in the assessment of species at-risk status based on quantitative decline criteria under International Union for Conservation of Nature (IUCN) and national endangered species legislation. Such data can be problematic for 3 reasons. First, statistical power to reject the null hypothesis of no change is often low because of small sample size and high sampling uncertainty leading to a high frequency of type II errors. Second, range-wide assessments composed of multiple site-specific observations do not effectively weight site-specific trends into global trends. Third, uncertainty in site-specific temporal trends and relative abundance are not propagated at the appropriate spatial scale. A common result is the propensity to underestimate the magnitude of declines and therefore fail to identify the appropriate at-risk status for a species. We used 3 statistical approaches, from simple to more complex, to estimate temporal decline rates for a designatable unit (DU) of rainbow trout in the Athabasca River watershed in western Canada. This DU is considered a native species for purposes of listing because of its genetic composition characterized as >0.95 indigenous origin in the face of continuing introgressive hybridization with introduced populations in the watershed. Analysis of abundance trends from 57 time series with a fixed-effects model identified 33 sites with negative trends, but only 2 were statistically significant. By contrast, a hierarchical linear mixed model weighted by site-specific abundance provided a DU-wide decline estimate of 16.4% per year and a 3-generation decline of 93.2%. A hierarchical Bayesian mixed model yielded a similar 3-generation decline trend of 91.3% and the posterior distribution showed that the estimate had a >99% probability of exceeding thresholds for an endangered listing. We conclude that the Bayesian approach was the most useful because it provided a probabilistic statement of threshold exceedance in support of an at-risk status recommendation.

El uso de datos extensivos, pero de baja resolución, de la abundancia es una práctica común en la evaluación del estado de riesgo de una especie con base en los criterios cuantitativos de declinación establecidos por la Unión Internacional para la Conservación de la Naturaleza (UICN) y la legislación nacional sobre especies en peligro extinción. Dicha información puede ser problemática por tres razones: primero, el poder estadístico para rechazar la hipótesis nula de ningún cambio es frecuentemente bajo debido a un tamaño pequeño de la muestra y a la elevada incertidumbre del muestreo, lo que resulta en una frecuencia elevada de errores de tipo II; segundo, las evaluaciones de amplia variedad compuestas de varias observaciones específicas de sitio no sopesan efectivamente las tendencias específicas de sitio dentro de las tendencias globales; y tercero, la incertidumbre en las tendencias temporales específicas de sitio y en la abundancia relativa no se propagan a la escala espacial apropiada. Un resultado común del uso de esta información es la propensión a subestimar la magnitud de las declinaciones, y por lo tanto equivocarse en la identificación del estado de riesgo apropiado para la especie. Usamos tres estrategias estadísticas, de simples a más complejas, para estimar las tasas de declinación temporal para una unidad designable (UD) de trucha arcoíris en la cuenca del río Athabasca al oeste de Canadá. Esta UD es considerada una especie nativa por razones de listado debido a su composición genética, caracterizada como >0-95 de origen nativo de frente a la continua hibridación introgresiva con poblaciones introducidas a la cuenca. El análisis de las tendencias de abundancia de 57 series de tiempo con un modelo de efectos fijos identificó 33 sitios con tendencias negativas, pero sólo dos fueron estadísticamente significativas. En contraste, un modelo lineal mixto de jerarquías sopesado por abundancia específica de sitio proporcionó una estimación de declinación en toda la UD de 16.4% año−1 y una declinación a tres generaciones de 93.2%. Un modelo bayesiano de jerarquías produjo una tendencia de declinación a tres generaciones de 91.3% y la distribución posterior mostró que el estimado tuvo una probabilidad >99% de exceder los umbrales para la categorización como especie en peligro. Concluimos que la estrategia bayesiana fue la más útil porque proporcionó una afirmación probabilística de la superación del umbral a favor de una recomendación de categorizar el estado como en riesgo.

Repurposed floxacins targeting RSK4 prevent chemoresistance and metastasis in lung and bladder cancer.

Lung and bladder cancers are mostly incurable because of the early development of drug resistance and metastatic dissemination. Hence, improved therapies that tackle these two processes are urgently needed to improve clinical outcome. We have identified RSK4 as a promoter of drug resistance and metastasis in lung and bladder cancer cells. Silencing this kinase, through either RNA interference or CRISPR, sensitized tumor cells to chemotherapy and hindered metastasis in vitro and in vivo in a tail vein injection model. Drug screening revealed several floxacin antibiotics as potent RSK4 activation inhibitors, and trovafloxacin reproduced all effects of RSK4 silencing in vitro and in/ex vivo using lung cancer xenograft and genetically engineered mouse models and bladder tumor explants. Through x-ray structure determination and Markov transient and Deuterium exchange analyses, we identified the allosteric binding site and revealed how this compound blocks RSK4 kinase activation through binding to an allosteric site and mimicking a kinase autoinhibitory mechanism involving the RSK4's hydrophobic motif. Last, we show that patients undergoing chemotherapy and adhering to prophylactic levofloxacin in the large placebo-controlled randomized phase 3 SIGNIFICANT trial had significantly increased (P = 0.048) long-term overall survival times. Hence, we suggest that RSK4 inhibition may represent an effective therapeutic strategy for treating lung and bladder cancer.

Targeting Mycobacterium tuberculosis CoaBC through Chemical Inhibition of 4'-Phosphopantothenoyl-l-cysteine Synthetase (CoaB) Activity.

Coenzyme A (CoA) is a ubiquitous cofactor present in all living cells and estimated to be required for up to 9% of intracellular enzymatic reactions. Mycobacterium tuberculosis (Mtb) relies on its own ability to biosynthesize CoA to meet the needs of the myriad enzymatic reactions that depend on this cofactor for activity. As such, the pathway to CoA biosynthesis is recognized as a potential source of novel tuberculosis drug targets. In prior work, we genetically validated CoaBC as a bactericidal drug target in Mtb in vitro and in vivo. Here, we describe the identification of compound 1f, a small molecule inhibitor of the 4'-phosphopantothenoyl-l-cysteine synthetase (PPCS; CoaB) domain of the bifunctional Mtb CoaBC, and show that this compound displays on-target activity in Mtb. Compound 1f was found to inhibit CoaBC uncompetitively with respect to 4'-phosphopantothenate, the substrate for the CoaB-catalyzed reaction. Furthermore, metabolomic profiling of wild-type Mtb H37Rv following exposure to compound 1f produced a signature consistent with perturbations in pantothenate and CoA biosynthesis. As the first report of a direct small molecule inhibitor of Mtb CoaBC displaying target-selective whole-cell activity, this study confirms the druggability of CoaBC and chemically validates this target.

Compartmentation of ß2 -adrenoceptor stimulated cAMP responses by phosphodiesterase types 2 and 3 in cardiac ventricular myocytes.

BACKGROUND AND PURPOSE: In cardiac myocytes, cyclic AMP (cAMP) produced by both ß1 - and ß2 -adrenoceptors increases L-type Ca2+ channel activity and myocyte contraction. However, only cAMP produced by ß1 -adrenoceptors enhances myocyte relaxation through phospholamban-dependent regulation of the sarco/endoplasmic reticulum Ca2+ ATPase 2 (SERCA2). Here we have tested the hypothesis that stimulation of ß2 -adrenoceptors produces a cAMP signal that is unable to reach SERCA2 and determine what role, if any, phosphodiesterase (PDE) activity plays in this compartmentation. EXPERIMENTAL APPROACH: The cAMP responses produced by ß1 -and ß2 -adrenoceptor stimulation were studied in adult rat ventricular myocytes using two different fluorescence resonance energy transfer (FRET)-based biosensors, the Epac2-camps, which is expressed uniformly throughout the cytoplasm of the entire cell and the Epac2-αKAP, which is targeted to the SERCA2 signalling complex. KEY RESULTS: Selective activation of ß1 - or ß2 -adrenoceptors produced cAMP responses detected by Epac2-camps. However, only stimulation of ß1 -adrenoceptors produced a cAMP response detected by Epac2-αKAP. Yet, stimulation of ß2 -adrenoceptors was able to produce a cAMP signal detected by Epac2-αKAP in the presence of selective inhibitors of PDE2 or PDE3, but not PDE4. CONCLUSION AND IMPLICATIONS: These results support the conclusion that cAMP produced by ß2 -adrenoceptor stimulation was not able to reach subcellular locations where the SERCA2 pump is located. Furthermore, this compartmentalized response is due at least in part to PDE2 and PDE3 activity. This discovery could lead to novel PDE-based therapeutic treatments aimed at correcting cardiac relaxation defects associated with certain forms of heart failure.

Inhibiting Mycobacterium tuberculosis CoaBC by targeting an allosteric site.

Coenzyme A (CoA) is a fundamental co-factor for all life, involved in numerous metabolic pathways and cellular processes, and its biosynthetic pathway has raised substantial interest as a drug target against multiple pathogens including Mycobacterium tuberculosis. The biosynthesis of CoA is performed in five steps, with the second and third steps being catalysed in the vast majority of prokaryotes, including M. tuberculosis, by a single bifunctional protein, CoaBC. Depletion of CoaBC was found to be bactericidal in M. tuberculosis. Here we report the first structure of a full-length CoaBC, from the model organism Mycobacterium smegmatis, describe how it is organised as a dodecamer and regulated by CoA thioesters. A high-throughput biochemical screen focusing on CoaB identified two inhibitors with different chemical scaffolds. Hit expansion led to the discovery of potent and selective inhibitors of M. tuberculosis CoaB, which we show to bind to a cryptic allosteric site within CoaB.

Discovery of an Allosteric Binding Site in Kinetoplastid Methionyl-tRNA Synthetase.

Methionyl-tRNA synthetase (MetRS) is a chemically validated drug target in kinetoplastid parasites Trypanosoma brucei and Leishmania donovani. To date, all kinetoplastid MetRS inhibitors described bind in a similar way to an expanded methionine pocket and an adjacent, auxiliary pocket. In the current study, we have identified a structurally novel class of inhibitors containing a 4,6-diamino-substituted pyrazolopyrimidine core (the MetRS02 series). Crystallographic studies revealed that MetRS02 compounds bind to an allosteric pocket in L. major MetRS not previously described, and enzymatic studies demonstrated a noncompetitive mode of inhibition. Homology modeling of the Trypanosoma cruzi MetRS enzyme revealed key differences in the allosteric pocket between the T. cruzi and Leishmania enzymes. These provide a likely explanation for the lower MetRS02 potencies that we observed for the T. cruzi enzyme compared to the Leishmania enzyme. The identification of a new series of MetRS inhibitors and the discovery of a new binding site in kinetoplastid MetRS enzymes provide a novel strategy in the search for new therapeutics for kinetoplastid diseases.

Lysyl-tRNA synthetase as a drug target in malaria and cryptosporidiosis.

Malaria and cryptosporidiosis, caused by apicomplexan parasites, remain major drivers of global child mortality. New drugs for the treatment of malaria and cryptosporidiosis, in particular, are of high priority; however, there are few chemically validated targets. The natural product cladosporin is active against blood- and liver-stage Plasmodium falciparum and Cryptosporidium parvum in cell-culture studies. Target deconvolution in P. falciparum has shown that cladosporin inhibits lysyl-tRNA synthetase (PfKRS1). Here, we report the identification of a series of selective inhibitors of apicomplexan KRSs. Following a biochemical screen, a small-molecule hit was identified and then optimized by using a structure-based approach, supported by structures of both PfKRS1 and C. parvum KRS (CpKRS). In vivo proof of concept was established in an SCID mouse model of malaria, after oral administration (ED90 = 1.5 mg/kg, once a day for 4 d). Furthermore, we successfully identified an opportunity for pathogen hopping based on the structural homology between PfKRS1 and CpKRS. This series of compounds inhibit CpKRS and C. parvum and Cryptosporidium hominis in culture, and our lead compound shows oral efficacy in two cryptosporidiosis mouse models. X-ray crystallography and molecular dynamics simulations have provided a model to rationalize the selectivity of our compounds for PfKRS1 and CpKRS vs. (human) HsKRS. Our work validates apicomplexan KRSs as promising targets for the development of drugs for malaria and cryptosporidiosis.

Identification of Leishmania major UDP-Sugar Pyrophosphorylase Inhibitors Using Biosensor-Based Small Molecule Fragment Library Screening.

Leishmaniasis is a neglected disease that is caused by different species of the protozoan parasite Leishmania, and it currently affects 12 million people worldwide. The antileishmanial therapeutic arsenal remains very limited in number and efficacy, and there is no vaccine for this parasitic disease. One pathway that has been genetically validated as an antileishmanial drug target is the biosynthesis of uridine diphosphate-glucose (UDP-Glc), and its direct derivative UDP-galactose (UDP-Gal). De novo biosynthesis of these two nucleotide sugars is controlled by the specific UDP-glucose pyrophosphorylase (UGP). Leishmania parasites additionally express a UDP-sugar pyrophosphorylase (USP) responsible for monosaccharides salvage that is able to generate both UDP-Gal and UDP-Glc. The inactivation of the two parasite pyrophosphorylases UGP and USP, results in parasite death. The present study reports on the identification of structurally diverse scaffolds for the development of USP inhibitors by fragment library screening. Based on this screening, we selected a small set of commercially available compounds, and identified molecules that inhibit both Leishmania major USP and UGP, with a half-maximal inhibitory concentration in the 100 µM range. The inhibitors were predicted to bind at allosteric regulation sites, which were validated by mutagenesis studies. This study sets the stage for the development of potent USP inhibitors.

Life-history variation along environmental and harvest clines of a northern freshwater fish: Plasticity and adaptation.

Plasticity, local adaptation and evolutionary trade-offs drive clinal variation in traits associated with lifetime growth. Disentangling the processes and determinants that cause these traits to vary helps to understand species' responses to changing environments. This is particularly urgent for exploited populations, where size-selective harvest can induce life-history evolution. Lake trout (Salvelinus namaycush) are an exploited fish with a life history adapted to low-productivity freshwaters of northern North America, which makes them highly vulnerable to ecosystem changes and overfishing. We characterized life-history variation across a broad and diverse landscape for this iconic northern freshwater fish and evaluated whether clinal variation was consistent with hypotheses for local adaptation or growth plasticity. We estimated growth-associated traits for 90 populations exposed to a diversity of environments using a Bayesian multivariate hierarchical model. We tested for clinal variation in their somatic growth, size at maturity and reproductive allocation along environmental gradients of lake productivity, climate, prey and exploitation clines under competing hypotheses of plasticity and local adaptation. Clinal life-history variation was consistent with growth plasticity and local adaptations but not harvest-induced evolution. Variation in somatic growth was explained by exploitation, climate and prey fish occurrence. Increased exploitation, from pristine to fully exploited conditions, led to increased somatic growth (from 32 to 45 mm/year) and adult life spans, and reduced age at maturity (from 11 to 8 years). Variation in size at maturity was explained by climate and, less certainly, prey fish occurrence, while reproductive allocation was explained by evolutionary trade-offs with mortality and other traits, but not environment. Lake trout life-history variation within this range was as wide as that observed across dozens of other freshwater species. Lake trout life histories resulted from evolutionary trade-offs, growth plasticity and local adaptations along several environmental clines. Presuming a plastic response, we documented ~1.4-fold growth compensation to exploitation-lower growth compensation than observed in many freshwater fishes. These results suggest that harvested species exposed to spatially structured and diverse environments may have substantial clinal variation on different traits, but due to different processes, and this has implications for their resilience and successful management.

Metabolic Scope as a Proximate Constraint on Individual Behavioral Variation: Effects on Personality, Plasticity, and Predictability.

Behavioral ecologists have hypothesized that among-individual differences in resting metabolic rate (RMR) may predict consistent individual differences in mean values for costly behaviors or for behaviors that affect energy intake rate. This hypothesis has empirical support and presently attracts considerable attention, but, notably, it does not provide predictions for individual differences in (a) behavioral plasticity or (b) unexplained variation (residual variation from mean individual behavior, here termed predictability). We outline how consideration of aerobic maximum metabolic rate (MMR) and particularly aerobic scope (= MMR - RMR) can be used to simultaneously make predictions about mean and among- and within-individual variation in behavior. We predict that while RMR should be proportional to an individual's mean level of sustained behavioral activity (one aspect of its personality), individuals with greater aerobic scope will also have greater scope to express behavioral plasticity and/or greater unpredictability in behavior (=greater residual variation). As a first step toward testing these predictions, we analyze existing activity data from selectively bred lines of mice that differ in both daily activity and aerobic scope. We find that replicate high-scope mice are more active on average and show greater among-individual variation in activity, greater among-individual variation in plasticity, and greater unpredictability. These data provide some tentative first support for our hypothesis, suggesting that further research on this topic would be valuable.

Habitat availability and ontogenetic shifts alter bottlenecks in size-structured fish populations.

For species that utilize different habitats throughout their life cycle, the habitat limitation at a given stage can act as a bottleneck on population abundance, impacting density-dependent processes such as individual growth and survival. We explore the influence of habitat limitation on population dynamics by developing a multi-stage population model based on lake-dwelling rainbow trout (Oncorhynchus mykiss) populations where adults occupy the lake habitat but use tributaries for spawning and juvenile rearing. The model details density-dependent ecological processes and ontogenetic habitat shifts, harvest mortality, and the impact of climate on growth. We ran model simulations using a range of early life stage habitat availabilities and climatic conditions representative of the native range of rainbow trout in Canada and compared the results to empirical data. The results suggest that (1) increases in early life stage habitat leads to increases in population abundance but, due to density-dependent processes, also results in slower growing stunted populations; (2) population bottlenecks can occur at any life stage, even at the adult stage if spawning and rearing habitats are abundant; (3) when the level of competition for early life stages is increased, inter-cohort competition can lead to population cycles. The model's conclusions are further reinforced by empirical data showing a similar trend in the relationship between fish density and maximum size and providing evidence that limited early life stage habitat leads to lower fish densities and larger fish size. We provide a model that links environmental conditions to population dynamics and is useful for fisheries management and habitat conservation decisions.

Interaction of ecological and angler processes: experimental stocking in an open access, spatially structured fishery.

Effective management of socioecological systems requires an understanding of the complex interactions between people and the environment. In recreational fisheries, which are prime examples of socioecological systems, anglers are analogous to mobile predators in natural predator-prey systems, and individual fisheries in lakes across a region are analogous to a spatially structured landscape of prey patches. Hence, effective management of recreational fisheries across large spatial scales requires an understanding of the dynamic interactions among ecological density dependent processes, landscape-level characteristics, and angler behaviors. We focused on the stocked component of the open access rainbow trout (Oncorhynchus mykiss) fishery in British Columbia (BC), and we used an experimental approach wherein we manipulated stocking densities in a subset of 34 lakes in which we monitored angler effort, fish abundance, and fish size for up to seven consecutive years. We used an empirically derived relationship between fish abundance and fish size across rainbow trout populations in BC to provide a measure of catch-based fishing quality that accounts for the size-abundance trade off in this system. We replicated our experimental manipulation in two regions known to have different angler populations and broad-scale access costs. We hypothesized that angler effort would respond to variation in stocking density, resulting in spatial heterogeneity in angler effort but homogeneity in catch-based fishing quality within regions. We found that there is an intermediate stocking density for a given lake or region at which angler effort is maximized (i.e., an optimal stocking density), and that this stocking density depends on latent effort and lake accessibility. Furthermore, we found no clear effect of stocking density on our measure of catch-based fishing quality, suggesting that angler effort homogenizes catch-related attributes leading to an eroded relationship between stocking density and catch-based fishing quality at the timescale of annual surveys. We conclude that declines in fishing quality resulting from understocking (due to declines in catch rate with low fish abundance) and overstocking (due to suppressed growth and limited recruitment at high density) give an optimal stocking rate that depends on accessibility and latent effort.

Supply-demand equilibria and the size-number trade-off in spatially structured recreational fisheries.

Recreational fishing effort varies across complex inland landscapes (e.g., lake-districts) and appears influenced by both angler preferences and qualities of the fishery resource, like fish size and abundance. However, fish size and abundance have an ecological trade-off within a population, thereby structuring equal-quality isopleths expressing this trade-off across the fishing landscape. Since expressed preferences of recreational anglers (i.e., site-selection of high-quality fishing opportunities among many lakes) can be analogous to optimal foraging strategies of natural predators, adopting such concepts can aid in understanding scale-dependence in fish-angler interactions and impacts of fishing across broad landscapes. Here, we assumed a fish supply-angler demand equilibria and adapted a novel bivariate measure of fishing quality based on fish size and catch rates to assess how recreational anglers influence fishing quality among a complex inland landscape. We then applied this metric to evaluate (1) angler preferences for caught and released fish compared to harvested fish, (2) the nonlinear size-numbers trade-off with uncertainty in both traits, and (3) the spatial-scale of the equilibria across 62 lakes and four independent management regions in British Columbia's (BC) rainbow trout Oncorhynchus mykiss fishery. We found anglers had low preference for caught and released fish (~10% of the value compared to harvested fish), which modified anglers' perception of fishing quality. Hence, fishing quality and angler effort was not influenced simply by total fish caught, but largely by harvested fish catch rates. Fishing quality varied from BC's northern regions (larger fish and more abundant) compared to southern regions (smaller fish and less abundant) directly associated with a 2.5 times increase in annual fishing effort in southern regions, suggesting that latent fishing pressure can structure the size-numbers trade-off in rainbow trout populations. The presence of two different equal-quality isopleths suggests at least two effective landscapes support co-occurring ideal free distributions of recreational fishing effort in BC's rainbow fishery. Anglers' expressed preferences among lakes interacted with density dependent growth and survival within lakes to structure a size-numbers trade-off influencing how anglers perceive fishing quality and, ultimately, distribute across complex inland landscapes.

Rainbow trout in seasonal environments: phenotypic trade-offs across a gradient in winter duration.

Survival through periods of resource scarcity depends on the balance between metabolic demands and energy storage. The opposing effects of predation and starvation mortality are predicted to result in trade-offs between traits that optimize fitness during periods of resource plenty (e.g., during the growing season) and those that optimize fitness during periods of resource scarcity (e.g., during the winter). We conducted a common environment experiment with two genetically distinct strains of rainbow trout to investigate trade-offs due to (1) the balance of growth and predation risk related to foraging rate during the growing season and (2) the allocation of energy to body size prior to the winter. Fry (age 0) from both strains were stocked into replicate natural lakes at low and high elevation that differed in winter duration (i.e., ice cover) by 59 days. Overwinter survival was lowest in the high-elevation lakes for both strains. Activity rate and growth rate were highest at high elevation, but growing season survival did not differ between strains or between environments. Hence, we did not observe a trade-off between growth and predation risk related to foraging rate. Growth rate also differed significantly between the strains across both environments, which suggests that growth rate is involved in local adaptation. There was not, however, a difference between strains or between environments in energy storage. Hence, we did not observe a trade-off between growth and storage. Our findings suggest that intrinsic metabolic rate, which affects a trade-off between growth rate and overwinter survival, may influence local adaptation in organisms that experience particularly harsh winter conditions (e.g., extended periods trapped beneath the ice in high-elevation lakes) in some parts of their range.

Energy allocation strategy modifies growth-survival trade-offs in juvenile fish across ecological and environmental gradients.

In young temperate zone fishes, conflicting energy demands lead to variability in growing season and winter survival. Growing season survival is driven by size-dependent predation risk whereas winter survival is constrained by autumn body size, energy storage and winter duration. We developed a model of the seasonality of energetics coupled to empirical measures of resource availability, size-dependent predation and temperature seasonality for rainbow trout (Oncorhynchus mykiss) in two sets of lakes in British Columbia, Canada, representing endpoints of a gradient of temperature, growing season duration and winter duration. This model was used to determine the energy allocation strategy which maximized first-year survival across these gradients. Survival was sensitive to the timing of the switch from somatic to storage strategies in cold, short growing season, low resource environments. A broader range of energy allocation strategies were viable in warmer, longer growing season and higher resource lakes. We used empirical observations of autumn energy storage and our modeled values for size-dependent minimal lipid levels needed to survive winter in each system to estimate winter survival for juvenile rainbow trout. Winter survival estimates were 6% in cold lakes with low resources, 82% in warm, lakes with low resources and 100% in warm lakes with high resources. Fish in warm lakes with ample resources allocated substantially more to storage than the minimum required to survive winter generated from our model, suggesting additional selection pressures for increased storage when there was ample surplus energy. We concluded that growth-survival trade-offs, modified by seasonality of the environment, influenced the growing season energy allocation strategies for young-of-the-year fish, and suggested this may be important for understanding population viability across environmental gradients.