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.

The Sacramento splittail (Pogonichthys macrolepidotus) is a species of special concern that is native to the San Francisco Estuary, USA. Two genetically distinct populations exist and differ in maximal salinity tolerances. We examined the expression of 12 genes representative of osmoregulatory functions in the gill over a 14  day time course at two different salinities [11 or 14 PSU (Practical Salinity Units)] and revealed that each population showed distinct patterns of gene expression consistent with population differences in response to osmotic regimes. The relatively more salinity-tolerant San Pablo population significantly upregulated nine out of the 12 transcripts investigated on day 1 of 11 PSU salinity exposure in comparison to the day zero freshwater control. Three transcripts (nka1a, nka1b, and mmp13) were differentially expressed between the populations at 7 and 14 days of salinity exposure, suggesting a reduced ability of the relatively salinity-intolerant Central Valley population to recover. Additionally, a phylogenetic analysis of several Sacramento splittail Na⁺/K⁺-ATPase α1 sequences resulted in grouping by proposed paralog rather than species, suggesting that different paralogs of this gene may exist. These findings, together with prior research conducted on the Sacramento splittail, suggest that the San Pablo population may be able to preferentially regulate select osmoregulatory genes, including different Na⁺/K⁺-ATPase α1 paralogs, to better cope with salinity challenges.

As global change alters multiple environmental conditions, predicting species' responses can be challenging without understanding how each environmental factor influences organismal performance. Approaches quantifying mechanistic relationships can greatly complement correlative field data, strengthening our abilities to forecast global change impacts. Substantial salinity increases are projected in the San Francisco Estuary, California, due to anthropogenic water diversion and climatic changes, where the critically endangered delta smelt (Hypomesus transpacificus) largely occurs in a low-salinity zone (LSZ), despite their ability to tolerate a much broader salinity range. In this study, we combined molecular and organismal measures to quantify the physiological mechanisms and sublethal responses involved in coping with salinity changes. Delta smelt utilize a suite of conserved molecular mechanisms to rapidly adjust their osmoregulatory physiology in response to salinity changes in estuarine environments. However, these responses can be energetically expensive, and delta smelt body condition was reduced at high salinities. Thus, acclimating to salinities outside the LSZ could impose energetic costs that constrain delta smelt's ability to exploit these habitats. By integrating data across biological levels, we provide key insight into the mechanistic relationships contributing to phenotypic plasticity and distribution limitations and advance the understanding of the molecular osmoregulatory responses in nonmodel estuarine fishes.

In estuary and coastal systems, human demand for freshwater, climate change-driven precipitation variability, and extreme weather impact salinity levels, reducing connectivity between mesohaline coastal fish populations and potentially contributing to genomic divergence. We examined gill transcriptome responses to salinity in wild-caught juveniles from two populations of Sacramento splittail (Pogonichthys macrolepidotus), a species of conservation concern that is endemic to the San Francisco Estuary, USA, and the lower reaches of its tributaries. Recent extreme droughts have led to salinities above the tolerance limits for this species, creating a migration barrier between these populations, which potentially contributed to population divergence. We identified transcripts involved in a conserved response to salinity; however, the more salinity-tolerant San Pablo population had greater transcriptome plasticity (3.6-fold more transcripts responded than the Central Valley population) and a response consistent with gill remodeling after 168 hr of exposure to elevated salinity. The reorganization of the gill in response to changing osmotic gradients is a process critical for acclimation and would facilitate enhanced salinity tolerance. We detected an upregulation of receptors that control the Wnt (wingless-type) cell signaling pathway that may be required for an adaptive response to increases in salinity, patterns not observed in the relatively salinity-sensitive Central Valley population. We detected 62 single nucleotide polymorphisms (SNPs) in coding regions of 26 transcripts that differed between the populations. Eight transcripts that contained SNPs were associated with immune responses, highlighting the importance of diversity in immune gene sequences as a defining characteristic of genomic divergence between these populations. Our data demonstrate that these populations have divergent transcriptomic responses to salinity, which is consistent with observed physiological differences in salinity tolerance.

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.

Fishes in estuarine waters are frequently exposed to treated wastewater effluent, among numerous other sources of contaminants, yet the impacts of these anthropogenic chemicals are not well understood in these dynamic and important waterways. Inland silversides (Menidia beryllina) at an early stage of development [12 days posthatch (dph)] were exposed to waters from two estuarine wastewater-treatment outfall locations in a tidal estuary, the Sacramento/San Joaquin Delta (California, USA) that had varied hydrology and input volumes. The genomic response caused by endocrine-disrupting compounds (EDCs) in these waters was determined using quantitative polymerase chain reaction on a suite of hormonally regulated genes. Relative androgenic and estrogenic activities of the waters were measured using CALUX reporter bioassays. The presence of bifenthrin, a pyrethroid pesticide and known EDC, as well as caffeine and the anti-inflammatory pharmaceutical ibuprofen, which were used as markers of wastewater effluent input, were determined using instrumental analysis. Detectable levels of bifenthrin (2.89 ng L(-1)) were found on one of the sampling dates, and caffeine was found on all sampling dates, in water from the Boynton Slough. Neither compound was detected at the Carquinez Strait site, which has a much smaller effluent discharge input volume relative to the receiving water body size compared with Boynton Slough. Water samples from both sites incubated in the CALUX cell line induced estrogenic and androgenic activity in almost all instances, though the estrogenicity was relatively higher than the androgenicity. Changes in the abundance of mRNA transcripts of endocrine-responsive genes and indicators of general chemical stress were observed after a 96-h exposure to waters from both locations. The relative levels of endocrine response, changes in gene transcript abundance, and contaminant concentrations were greater in water from the Boynton Slough site despite those effluents undergoing a more advanced treatment process. The availability of a widely geographically distributed estuarine model species (M. beryllina) now allows for improved assessment of treated effluent impacts across brackish, estuarine, and marine environments.