Divergent sensory and immune gene evolution in sea turtles with contrasting demographic and life histories.

Sea turtles represent an ancient lineage of marine vertebrates that evolved from terrestrial ancestors over 100 Mya. The genomic basis of the unique physiological and ecological traits enabling these species to thrive in diverse marine habitats remains largely unknown. Additionally, many populations have drastically declined due to anthropogenic activities over the past two centuries, and their recovery is a high global conservation priority. We generated and analyzed high-quality reference genomes for the leatherback (Dermochelys coriacea) and green (Chelonia mydas) turtles, representing the two extant sea turtle families. These genomes are highly syntenic and homologous, but localized regions of noncollinearity were associated with higher copy numbers of immune, zinc-finger, and olfactory receptor (OR) genes in green turtles, with ORs related to waterborne odorants greatly expanded in green turtles. Our findings suggest that divergent evolution of these key gene families may underlie immunological and sensory adaptations assisting navigation, occupancy of neritic versus pelagic environments, and diet specialization. Reduced collinearity was especially prevalent in microchromosomes, with greater gene content, heterozygosity, and genetic distances between species, supporting their critical role in vertebrate evolutionary adaptation. Finally, diversity and demographic histories starkly contrasted between species, indicating that leatherback turtles have had a low yet stable effective population size, exhibit extremely low diversity compared with other reptiles, and harbor a higher genetic load compared with green turtles, reinforcing concern over their persistence under future climate scenarios. These genomes provide invaluable resources for advancing our understanding of evolution and conservation best practices in an imperiled vertebrate lineage.

Smaller body size under warming is not due to gill-oxygen limitation in a cold-water salmonid.

Declining body size in fishes and other aquatic ectotherms associated with anthropogenic climate warming has significant implications for future fisheries yields, stock assessments and aquatic ecosystem stability. One proposed mechanism seeking to explain such body-size reductions, known as the gill oxygen limitation (GOL) hypothesis, has recently been used to model future impacts of climate warming on fisheries but has not been robustly empirically tested. We used brook trout (Salvelinus fontinalis), a fast-growing, cold-water salmonid species of broad economic, conservation and ecological value, to examine the GOL hypothesis in a long-term experiment quantifying effects of temperature on growth, resting metabolic rate (RMR), maximum metabolic rate (MMR) and gill surface area (GSA). Despite significantly reduced growth and body size at an elevated temperature, allometric slopes of GSA were not significantly different than 1.0 and were above those for RMR and MMR at both temperature treatments (15°C and 20°C), contrary to GOL expectations. We also found that the effect of temperature on RMR was time-dependent, contradicting the prediction that heightened temperatures increase metabolic rates and reinforcing the importance of longer-term exposures (e.g. >6 months) to fully understand the influence of acclimation on temperature-metabolic rate relationships. Our results indicate that although oxygen limitation may be important in some aspects of temperature-body size relationships and constraints on metabolic supply may contribute to reduced growth in some cases, it is unlikely that GOL is a universal mechanism explaining temperature-body size relationships in aquatic ectotherms. We suggest future research focus on alternative mechanisms underlying temperature-body size relationships, and that projections of climate change impacts on fisheries yields using models based on GOL assumptions be interpreted with caution.

Species and population specific gene expression in blood transcriptomes of marine turtles.

BACKGROUND: Transcriptomic data has demonstrated utility to advance the study of physiological diversity and organisms' responses to environmental stressors. However, a lack of genomic resources and challenges associated with collecting high-quality RNA can limit its application for many wild populations. Minimally invasive blood sampling combined with de novo transcriptomic approaches has great potential to alleviate these barriers. Here, we advance these goals for marine turtles by generating high quality de novo blood transcriptome assemblies to characterize functional diversity and compare global transcriptional profiles between tissues, species, and foraging aggregations. RESULTS: We generated high quality blood transcriptome assemblies for hawksbill (Eretmochelys imbricata), loggerhead (Caretta caretta), green (Chelonia mydas), and leatherback (Dermochelys coriacea) turtles. The functional diversity in assembled blood transcriptomes was comparable to those from more traditionally sampled tissues. A total of 31.3% of orthogroups identified were present in all four species, representing a core set of conserved genes expressed in blood and shared across marine turtle species. We observed strong species-specific expression of these genes, as well as distinct transcriptomic profiles between green turtle foraging aggregations that inhabit areas of greater or lesser anthropogenic disturbance. CONCLUSIONS: Obtaining global gene expression data through non-lethal, minimally invasive sampling can greatly expand the applications of RNA-sequencing in protected long-lived species such as marine turtles. The distinct differences in gene expression signatures between species and foraging aggregations provide insight into the functional genomics underlying the diversity in this ancient vertebrate lineage. The transcriptomic resources generated here can be used in further studies examining the evolutionary ecology and anthropogenic impacts on marine turtles.

Environment and phenology shape local adaptation in thermal performance.

Populations within species often exhibit variation in traits that reflect local adaptation and further shape existing adaptive potential for species to respond to climate change. However, our mechanistic understanding of how the environment shapes trait variation remains poor. Here, we used common garden experiments to quantify thermal performance in eight populations of the marine snail Urosalpinx cinerea across thermal gradients on the Atlantic and the Pacific coasts of North America. We then evaluated the relationship between thermal performance and environmental metrics derived from time-series data. Our results reveal a novel pattern of 'mixed' trait performance adaptation, where thermal optima were positively correlated with spawning temperature (cogradient variation), while maximum trait performance was negatively correlated with season length (countergradient variation). This counterintuitive pattern probably arises because of phenological shifts in the spawning season, whereby 'cold' populations delay spawning until later in the year when temperatures are warmer compared to 'warm' populations that spawn earlier in the year when temperatures are cooler. Our results show that variation in thermal performance can be shaped by multiple facets of the environment and are linked to organismal phenology and natural history. Understanding the impacts of climate change on organisms, therefore, requires the knowledge of how climate change will alter different aspects of the thermal environment.

The utility of transcriptomics in fish conservation.

There is growing recognition of the need to understand the mechanisms underlying organismal resilience (i.e. tolerance, acclimatization) to environmental change to support the conservation management of sensitive and economically important species. Here, we discuss how functional genomics can be used in conservation biology to provide a cellular-level understanding of organismal responses to environmental conditions. In particular, the integration of transcriptomics with physiological and ecological research is increasingly playing an important role in identifying functional physiological thresholds predictive of compensatory responses and detrimental outcomes, transforming the way we can study issues in conservation biology. Notably, with technological advances in RNA sequencing, transcriptome-wide approaches can now be applied to species where no prior genomic sequence information is available to develop species-specific tools and investigate sublethal impacts that can contribute to population declines over generations and undermine prospects for long-term conservation success. Here, we examine the use of transcriptomics as a means of determining organismal responses to environmental stressors and use key study examples of conservation concern in fishes to highlight the added value of transcriptome-wide data to the identification of functional response pathways. Finally, we discuss the gaps between the core science and policy frameworks and how thresholds identified through transcriptomic evaluations provide evidence that can be more readily used by resource managers.

Effects of high temperatures on threatened estuarine fishes during periods of extreme drought.

Climate change and associated increases in water temperatures may impact physiological performance in ectotherms and exacerbate endangered species declines. We used an integrative approach to assess the impact of elevated water temperature on two fishes of immediate conservation concern in a large estuary system, the threatened longfin smelt (Spirinchus thaleichthys) and endangered delta smelt (Hypomesus transpacificus). Abundances have reached record lows in California, USA, and these populations are at imminent risk of extirpation. California is currently impacted by a severe drought, resulting in high water temperatures, conditions that will become more common as a result of climate change. We exposed fish to environmentally relevant temperatures (14°C and 20°C) and used RNA sequencing to examine the transcriptome-wide responses to elevated water temperature in both species. Consistent with having a lower temperature tolerance, longfin smelt exhibited a pronounced cellular stress response, with an upregulation of heat shock proteins, after exposure to 20°C that was not observed in delta smelt. We detected an increase in metabolic rate in delta smelt at 20°C and increased expression of genes involved in metabolic processes and protein synthesis, patterns not observed in longfin smelt. Through examination of responses across multiple levels of biological organization, and by linking these responses to habitat distributions in the wild, we demonstrate that longfin smelt may be more susceptible than delta smelt to increases in temperatures, and they have little room to tolerate future warming in California. Understanding the species-specific physiological responses of sensitive species to environmental stressors is crucial for conservation efforts and managing aquatic systems globally.

Linking transcriptional responses to organismal tolerance reveals mechanisms of thermal sensitivity in a mesothermal endangered fish.

Forecasting species' responses to climate change requires understanding the underlying mechanisms governing environmental stress tolerance, including acclimation capacity and acute stress responses. Current knowledge of these physiological processes in aquatic ectotherms is largely drawn from eurythermal or extreme stenothermal species. Yet many species of conservation concern exhibit tolerance windows and acclimation capacities in between these extremes. We linked transcriptome profiles to organismal tolerance in a mesothermal endangered fish, the delta smelt (Hypomesus transpacificus), to quantify the cellular processes, sublethal thresholds and effects of thermal acclimation on acute stress responses. Delta smelt initiated rapid molecular changes in line with expectations of theoretical thermal limitation models, but also exhibited diminished capacity to modify the expression of some genes and cellular mechanisms key to coping with acute thermal stress found in eurytherms. Sublethal critical thresholds occurred 4-6 °C below their upper tolerance limits, and thermal acclimation shifted the onset of acute thermal stress and tolerance as predicted. However, we found evidence that delta smelt's limited thermal plasticity may be partially due to an inability of individuals to effectively make physiological adjustments to truly achieve new homoeostasis under heightened temperatures, resulting in chronic thermal stress. These findings provide insight into the physiological basis of the diverse patterns of thermal tolerances observed in nature. Moreover, understanding how underlying molecular mechanisms shape thermal acclimation capacity, acute stress responses and ultimately differential phenotypes contributes to a predictive framework to deduce species' responses in situ to changes in selective pressures due to climate change.

A Scientist's Guide to Achieving Broader Impacts through K-12 STEM Collaboration.

The National Science Foundation and other funding agencies are increasingly requiring broader impacts in grant applications to encourage US scientists to contribute to science education and society. Concurrently, national science education standards are using more inquiry-based learning (IBL) to increase students' capacity for abstract, conceptual thinking applicable to real-world problems. Scientists are particularly well suited to engage in broader impacts via science inquiry outreach, because scientific research is inherently an inquiry-based process. We provide a practical guide to help scientists overcome obstacles that inhibit their engagement in K-12 IBL outreach and to attain the accrued benefits. Strategies to overcome these challenges include scaling outreach projects to the time available, building collaborations in which scientists' research overlaps with curriculum, employing backward planning to target specific learning objectives, encouraging scientists to share their passion, as well as their expertise with students, and transforming institutional incentives to support scientists engaging in educational outreach.

Accounting for unobserved population dynamics and aging error in close-kin mark-recapture assessments.

Obtaining robust estimates of population abundance is a central challenge hindering the conservation and management of many threatened and exploited species. Close-kin mark-recapture (CKMR) is a genetics-based approach that has strong potential to improve the monitoring of data-limited species by enabling estimates of abundance, survival, and other parameters for populations that are challenging to assess. However, CKMR models have received limited sensitivity testing under realistic population dynamics and sampling scenarios, impeding the application of the method in population monitoring programs and stock assessments. Here, we use individual-based simulation to examine how unmodeled population dynamics and aging uncertainty affect the accuracy and precision of CKMR parameter estimates under different sampling strategies. We then present adapted models that correct the biases that arise from model misspecification. Our results demonstrate that a simple base-case CKMR model produces robust estimates of population abundance with stable populations that breed annually; however, if a population trend or non-annual breeding dynamics are present, or if year-specific estimates of abundance are desired, a more complex CKMR model must be constructed. In addition, we show that CKMR can generate reliable abundance estimates for adults from a variety of sampling strategies, including juvenile-focused sampling where adults are never directly observed (and aging error is minimal). Finally, we apply a CKMR model that has been adapted for population growth and intermittent breeding to two decades of genetic data from juvenile lemon sharks (Negaprion brevirostris) in Bimini, Bahamas, to demonstrate how application of CKMR to samples drawn solely from juveniles can contribute to monitoring efforts for highly mobile populations. Overall, this study expands our understanding of the biological factors and sampling decisions that cause bias in CKMR models, identifies key areas for future inquiry, and provides recommendations that can aid biologists in planning and implementing an effective CKMR study, particularly for long-lived data-limited species.

Sources of bias in applying close-kin mark-recapture to terrestrial game species with different life histories.

Close-kin mark-recapture (CKMR) is a method analogous to traditional mark-recapture but without requiring recapture of individuals. Instead, multilocus genotypes (genetic marks) are used to identify related individuals in one or more sampling occasions, which enables the opportunistic use of samples from harvested wildlife. To apply the method accurately, it is important to build appropriate CKMR models that do not violate assumptions linked to the species' and population's biology and sampling methods. In this study, we evaluated the implications of fitting overly simplistic CKMR models to populations with complex reproductive success dynamics or selective sampling. We used forward-in-time, individual-based simulations to evaluate the accuracy and precision of CKMR abundance and survival estimates in species with different longevities, mating systems, and sampling strategies. Simulated populations approximated a range of life histories among game species of North America with lethal sampling to evaluate the potential of using harvested samples to estimate population size. Our simulations show that CKMR can yield nontrivial biases in both survival and abundance estimates, unless influential life history traits and selective sampling are explicitly accounted for in the modeling framework. The number of kin pairs observed in the sample, in combination with the type of kinship used in the model (parent-offspring pairs and/or half-sibling pairs), can affect the precision and/or accuracy of the estimates. CKMR is a promising method that will likely see an increasing number of applications in the field as costs of genetic analysis continue to decline. Our work highlights the importance of applying population-specific CKMR models that consider relevant demographic parameters, individual covariates, and the protocol through which individuals were sampled.

Benchmarking ultra-high molecular weight DNA preservation methods for long-read and long-range sequencing.

BACKGROUND: Studies in vertebrate genomics require sampling from a broad range of tissue types, taxa, and localities. Recent advancements in long-read and long-range genome sequencing have made it possible to produce high-quality chromosome-level genome assemblies for almost any organism. However, adequate tissue preservation for the requisite ultra-high molecular weight DNA (uHMW DNA) remains a major challenge. Here we present a comparative study of preservation methods for field and laboratory tissue sampling, across vertebrate classes and different tissue types. RESULTS: We find that storage temperature was the strongest predictor of uHMW fragment lengths. While immediate flash-freezing remains the sample preservation gold standard, samples preserved in 95% EtOH or 20-25% DMSO-EDTA showed little fragment length degradation when stored at 4°C for 6 hours. Samples in 95% EtOH or 20-25% DMSO-EDTA kept at 4°C for 1 week after dissection still yielded adequate amounts of uHMW DNA for most applications. Tissue type was a significant predictor of total DNA yield but not fragment length. Preservation solution had a smaller but significant influence on both fragment length and DNA yield. CONCLUSION: We provide sample preservation guidelines that ensure sufficient DNA integrity and amount required for use with long-read and long-range sequencing technologies across vertebrates. Our best practices generated the uHMW DNA needed for the high-quality reference genomes for phase 1 of the Vertebrate Genomes Project, whose ultimate mission is to generate chromosome-level reference genome assemblies of all ∼70,000 extant vertebrate species.

Diminished warming tolerance and plasticity in low-latitude populations of a marine gastropod.

Models of species response to climate change often assume that physiological traits are invariant across populations. Neglecting potential intraspecific variation may overlook the possibility that some populations are more resilient or susceptible than others, creating inaccurate predictions of climate impacts. In addition, phenotypic plasticity can contribute to trait variation and may mediate sensitivity to climate. Quantifying such forms of intraspecific variation can improve our understanding of how climate can affect ecologically important species, such as invasive predators. Here, we quantified thermal performance (tolerance, acclimation capacity, developmental traits) across seven populations of the predatory marine snail (Urosalpinx cinerea) from native Atlantic and non-native Pacific coast populations in the USA. Using common garden experiments, we assessed the effects of source population and developmental acclimation on thermal tolerance and developmental traits of F1 snails. We then estimated climate sensitivity by calculating warming tolerance (thermal tolerance - habitat temperature), using field environmental data. We report that low-latitude populations had greater thermal tolerance than their high latitude counterparts. However, these same low-latitude populations exhibited decreased thermal tolerance when exposed to environmentally realistic higher acclimation temperatures. Low-latitude native populations had the greatest climate sensitivity (habitat temperatures near thermal limits). In contrast, invasive Pacific snails had the lowest climate sensitivity, suggesting that these populations are likely to persist and drive negative impacts on native biodiversity. Developmental rate significantly increased in embryos sourced from populations with greater habitat temperature but had variable effects on clutch size and hatching success. Thus, warming can produce widely divergent responses within the same species, resulting in enhanced impacts in the non-native range and extirpation in the native range. Broadly, our results highlight how intraspecific variation can alter management decisions, as this may clarify whether management efforts should be focused on many or only a few populations.

Plasma proteomics of green turtles (Chelonia mydas) reveals pathway shifts and potential biomarker candidates associated with health and disease.

Evaluating sea turtle health can be challenging due to an incomplete understanding of pathophysiologic responses in these species. Proteome characterization of clinical plasma samples can provide insights into disease progression and prospective biomarker targets. A TMT-10-plex-LC-MS/MS platform was used to characterize the plasma proteome of five, juvenile, green turtles (Chelonia mydas) and compare qualitative and quantitative protein changes during moribund and recovered states. The 10 plasma samples yielded a total of 670 unique proteins. Using ≥1.2-fold change in protein abundance as a benchmark for physiologic upregulation or downregulation, 233 (34.8%) were differentially regulated in at least one turtle between moribund and recovered states. Forty-six proteins (6.9%) were differentially regulated in all five turtles with two proteins (0.3%) demonstrating a statistically significant change. A principle component analysis showed protein abundance loosely clustered between moribund samples or recovered samples and for turtles that presented with trauma (n = 3) or as intestinal floaters (n = 2). Gene Ontology terms demonstrated that moribund samples were represented by a higher number of proteins associated with blood coagulation, adaptive immune responses and acute phase response, while recovered turtle samples included a relatively higher number of proteins associated with metabolic processes and response to nutrients. Abundance levels of 48 proteins (7.2%) in moribund samples significantly correlated with total protein, albumin and/or globulin levels quantified by biochemical analysis. Differentially regulated proteins identified with immunologic and physiologic functions are discussed for their possible role in the green turtle pathophysiologic response and for their potential use as diagnostic biomarkers. These findings enhance our ability to interpret sea turtle health and further progress conservation, research and rehabilitation programs for these ecologically important species.

Testing the Stability of Plasma Protein and Whole Blood RNA in Archived Blood of Loggerhead Sea Turtles, Caretta caretta.

Sample storage conditions can affect accuracy and reproducibility of biological measurements. Storing samples rapidly at the lowest available temperatures is considered ideal but is not always feasible when sampling in remote and logistically challenging field conditions, as is often the case with sea turtles. The objective of this study was to examine the stability of plasma proteins and quality of whole blood RNA from loggerhead sea turtle samples collected as part of an eighteen-year-long curated specimen collection. These biological variables are often used to assess sea turtle health; therefore, it is necessary to maintain the integrity of these components during storage. Protein electrophoresis was conducted on heparinized plasma from individual turtles collected in 2018 (n = 3), 2008 (n = 3), and 2001 (n = 3). Plasma was also pooled from four turtles sampled in 2018 and subjected to various storage temperatures. Whole blood was collected in blood collection tubes containing sodium heparin or PAXgene tubes with an RNA preservative. These were subjected to different storage treatments that can possibly occur during logistically difficult field sampling. Following various treatments, plasma proteins showed minor differences across collection years and no differences among storage treatments were observed, even when exposed to 38°C for three hours. RNA quality was assessed from whole blood using an RNA integrity number (RIN). RINs were poor from sodium heparin tubes that were frozen and from PAXgene tubes after an extended thaw. High-quality RNA was obtained from sodium heparin tubes that were never frozen and from PAXgene tubes with freezing delayed by up to 11 days. Overall, these results indicate that plasma proteins remain stable over time and when exposed to undesirable storage conditions, and RNA degrades rapidly in sea turtle blood after freezing and when not properly preserved. These aspects are important to consider when planning sampling protocols and logistics for optimal long-term sample preservation.

Persistent organic pollutants in green sea turtles (Chelonia mydas) inhabiting two urbanized Southern California habitats.

Within Southern California, east Pacific green sea turtles (Chelonia mydas) forage year-round, taking advantage of diverse food resources, including seagrass, marine algae, and invertebrates. Assessing persistent organic pollutants (POP) in green turtle aggregations in the Seal Beach National Wildlife Refuge (SBNWR, n = 17) and San Diego Bay (SDB, n = 25) can help quantify contamination risks for these populations. Blood plasma was analyzed for polychlorinated biphenyls (PCBs), organochlorinated pesticides (OCPs), and polybrominated diphenyl ethers (PBDEs). PCBs and body size explained much of the separation of turtles by foraging aggregation in a principal component analysis. Turtles from SDB had significantly (p < 0.001) higher total PCBs than SBNWR turtles. Most PCBs detected in turtles were non-dioxin-like PCB congeners (153, 138, 99) that are associated with neurotoxicity. Recaptured turtles' POP levels changed significantly over time indicating significant variation in POP levels through time and space, even among adjacent foraging locations.

A versatile Rapture (RAD-Capture) platform for genotyping marine turtles.

Advances in high-throughput sequencing (HTS) technologies coupled with increased interdisciplinary collaboration are rapidly expanding capacity in the scope and scale of wildlife genetic studies. While existing HTS methods can be directly applied to address some evolutionary and ecological questions, certain research goals necessitate tailoring methods to specific study organisms, such as high-throughput genotyping of the same loci that are comparable over large spatial and temporal scales. These needs are particularly common for studies of highly mobile species of conservation concern like marine turtles, where life history traits, limited financial resources and other constraints require affordable, adaptable methods for HTS genotyping to meet a variety of study goals. Here, we present a versatile marine turtle HTS targeted enrichment platform adapted from the recently developed Rapture (RAD-Capture) method specifically designed to meet these research needs. Our results demonstrate consistent enrichment of targeted regions throughout the genome and discovery of candidate variants in all species examined for use in various conservation genetics applications. Accurate species identification confirmed the ability of our platform to genotype over 1,000 multiplexed samples and identified areas for future methodological improvement such as optimization for low initial concentration samples. Finally, analyses within green turtles supported the ability of this platform to identify informative SNPs for stock structure, population assignment and other applications over a broad geographic range of interest to management. This platform provides an additional tool for marine turtle genetic studies and broadens capacity for future large-scale initiatives such as collaborative global marine turtle genetic databases.

Trace metals in green sea turtles (Chelonia mydas) inhabiting two southern California coastal estuaries.

Foraging aggregations of east Pacific green sea turtles (Chelonia mydas) inhabit the Seal Beach National Wildlife Refuge (SBNWR) and San Diego Bay (SDB), two habitats in southern California, USA, located near urbanized areas. Both juvenile and adult green turtles forage in these areas and exhibit high site fidelity, which potentially exposes green turtles to anthropogenic contaminants. We assessed 21 trace metals (TM) bioaccumulated in green turtle scute and red blood cell (RBC) samples collected from SBNWR (n = 16 turtles) and SDB (n = 20 turtles) using acid digestion and inductively coupled plasma mass spectrometry. Principal component analyses of TM composition indicate that SBNWR and SDB turtles have location-specific contaminant signatures, characterized by differences in cadmium and selenium concentrations: SBNWR turtles had significantly more cadmium and selenium in RBC and more selenium in scute samples, than SDB turtles. Cadmium and selenium concentrations in RBC had a strong positive relationship, regardless of location. SBNWR turtles had higher selenium in RBCs than previously measured in other green turtle populations globally. Due to different retention times in blood vs. scute, these results suggest that SBNWR turtles have high long- and short-term selenium exposure. Turtles from SBNWR and SDB had higher trace metal concentrations than documented in green turtle populations that inhabit non-urbanized areas, supporting the hypothesis that coastal cities can increase trace metal exposure to local green turtles. Our study finds evidence that green turtle TM concentrations can differ between urbanized habitats and that long-term monitoring of these green turtles may be necessary.

Juvenile rockfish show resilience to CO2-acidification and hypoxia across multiple biological scales.

California's coastal ecosystems are forecasted to undergo shifting ocean conditions due to climate change, some of which may negatively impact recreational and commercial fish populations. To understand if fish populations have the capacity to respond to multiple stressors, it is critical to examine interactive effects across multiple biological scales, from cellular metabolism to species interactions. This study examined the effects of CO2-acidification and hypoxia on two naturally co-occurring species, juvenile rockfish (genus Sebastes) and a known predator, cabezon (Scorpaenichthys marmoratus). Fishes were exposed to two PCO2 levels at two dissolved oxygen (DO) levels: ~600 (ambient) and ~1600 (high) µatm PCO2 and 8.0 (normoxic) and 4.5 mg l-1 DO (hypoxic) and assessments of cellular metabolism, prey behavior and predation mortality rates were quantified after 1 and 3 weeks. Physiologically, rockfish showed acute alterations in cellular metabolic enzyme activity after 1 week of acclimation to elevated PCO2 and hypoxia that were not evident in cabezon. Alterations in rockfish energy metabolism were driven by increases in anaerobic LDH activity, and adjustments in enzyme activity ratios of cytochrome c oxidase and citrate synthase and LDH:CS. Correlated changes in rockfish behavior were also apparent after 1 week of acclimation to elevated PCO2 and hypoxia. Exploration behavior increased in rockfish exposed to elevated PCO2 and spatial analysis of activity indicated short-term interference with anti-predator responses. Predation rate after 1 week increased with elevated PCO2; however, no mortality was observed under the multiple-stressor treatment suggesting negative effects on cabezon predators. Most noteworthy, metabolic and behavioral changes were moderately compensated after 3 weeks of acclimation, and predation mortality rates also decreased suggesting that these rockfish may be resilient to changes in environmental stressors predicted by climate models. Linking physiological and behavioral responses to multiple stressors is vital to understand impacts on populations and community dynamics.

Assessments at multiple levels of biological organization allow for an integrative determination of physiological tolerances to turbidity in an endangered fish species.

Turbidity can influence trophic levels by altering species composition and can potentially affect fish feeding strategies and predator-prey interactions. The estuarine turbidity maximum, described as an area of increased suspended particles, phytoplankton and zooplankton, generally represents a zone with higher turbidity and enhanced food sources important for successful feeding and growth in many fish species. The delta smelt (Hypomesus transpacificus) is an endangered, pelagic fish species endemic to the San Francisco Estuary and Sacramento-San Joaquin River Delta, USA, where it is associated with turbid waters. Turbidity is known to play an important role for the completion of the species' life cycle; however, turbidity ranges in the Delta are broad, and specific requirements for this fish species are still unknown. To evaluate turbidity requirements for early life stages, late-larval delta smelt were maintained at environmentally relevant turbidity levels ranging from 5 to 250 nephelometric turbidity units (NTU) for 24 h, after which a combination of physiological endpoints (molecular biomarkers and cortisol), behavioural indices (feeding) and whole-organism measures (survival) were determined. All endpoints delivered consistent results and identified turbidities between 25 and 80 NTU as preferential. Delta smelt survival rates were highest between 12 and 80 NTU and feeding rates were highest between 25 and 80 NTU. Cortisol levels indicated minimal stress between 35 and 80 NTU and were elevated at low turbidities (5, 12 and 25 NTU). Expression of stress-related genes indicated significant responses for gst, hsp70 and glut2 in high turbidities (250 NTU), and principal component analysis on all measured genes revealed a clustering of 25, 35, 50 and 80 NTU separating the medium-turbidity treatments from low- and high-turbidity treatments. Taken together, these data demonstrate that turbidity levels that are either too low or too high affect delta smelt physiological performance, causing significant effects on overall stress, food intake and mortality. They also highlight the need for turbidity to be considered in habitat and water management decisions.