The utility of transcriptomics in fish conservation

Integrated analysis of intestinal microbiota and transcriptome reveals that a coordinated interaction of the endocrine, immune system and gut microbiota response to heat stress in Litopenaeus vannamei

Due to the ongoing global warming, the risk of heatwaves in the oceans is continuously increasing while our understanding of the physiological response of Litopenaeus vannamei under extreme temperature conditions remains limited. Therefore, this study aimed to evaluate the physiological responses of L. vannamei under heat stress. Our results indicated that as temperature rose, the structure of intestinal and hepatopancreatic tissues was damaged sequentially. Activity of immune-related enzymes (acid phosphatase/alkaline phosphatase) initially increased before decreased, while antioxidant enzymes (superoxide dismutase and glutathione-S transferase) activity and malondialdehyde content increased with rising temperature. In addition, the total antioxidant capacity decreased with rising temperature. With the rising temperature, there was a significant increase in the expression of caspase-3, heat shock protein 70, lipopolysaccharide-induced tumor necrosis factor-α, transcriptional enhanced associate domain and yorkie in intestinal and hepatopancreatic tissues. Following heat stress, the number of potentially beneficial bacteria (Rhodobacteraceae and Gemmonbacter) increased which maintain balance and promote vitamin synthesis. Intestinal transcriptome analysis revealed 852 differentially expressed genes in the heat stress group compared with the control group. KEGG functional annotation results showed that the endocrine system was the most abundant in Organismal systems followed by the immune system. These results indicated that heat stress leads to tissue damage in shrimp, however the shrimp may respond to stress through a coordinated interaction strategy of the endocrine system, immune system and gut microbiota. This study revealed the response mechanism of L. vannamei to acute heat stress and potentially provided a theoretical foundation for future research on shrimp environmental adaptations.

Unraveling the molecular mechanisms of fish physiological response to freshwater salinization: A comparative multi-tissue transcriptomic study in a river polluted by potash mining

Freshwater salinization is an escalating global environmental issue that threatens freshwater biodiversity, including fish populations. This study aims to uncover the molecular basis of salinity physiological responses in a non-native minnow species (Phoxinus septimaniae x P. dragarum) exposed to saline effluents from potash mines in the Llobregat River, Barcelona, Spain. Employing high-throughput mRNA sequencing and differential gene expression analyses, brain, gills, and liver tissues collected from fish at two stations (upstream and downstream of saline effluent discharge) were examined. Salinization markedly influenced global gene expression profiles, with the brain exhibiting the most differentially expressed genes, emphasizing its unique sensitivity to salinity fluctuations. Pathway analyses revealed the expected enrichment of ion transport and osmoregulation pathways across all tissues. Furthermore, tissue-specific pathways associated with stress, reproduction, growth, immune responses, methylation, and neurological development were identified in the context of salinization. Rigorous validation of RNA-seq data through quantitative PCR (qPCR) underscored the robustness and consistency of our findings across platforms. This investigation unveils intricate molecular mechanisms steering salinity physiological response in non-native minnows confronting diverse environmental stressors. This comprehensive analysis sheds light on the underlying genetic and physiological mechanisms governing fish physiological response in salinity-stressed environments, offering essential knowledge for the conservation and management of freshwater ecosystems facing salinization.

A multi-tissue de novo transcriptome assembly and relative gene expression of the vulnerable freshwater salmonid Thymallus ligericus

Freshwater ecosystems are among the most endangered ecosystems worldwide. While numerous taxa are on the verge of extinction as a result of global changes and direct or indirect anthropogenic activity, genomic and transcriptomic resources represent a key tool for comprehending species' adaptability and serve as the foundation for conservation initiatives. The Loire grayling, Thymallus ligericus, is a freshwater European salmonid endemic to the upper Loire River basin. The species is comprised of fragmented populations that are dispersed over a small area and it has been identified as a vulnerable species. Here, we provide a multi-tissue de novo transcriptome assembly of T. ligericus. The completeness and integrity of the transcriptome were assessed before and after redundancy removal with lineage-specific libraries from Eukaryota, Metazoa, Vertebrata, and Actinopterygii. Relative gene expression was assessed for each of the analyzed tissues, using the de novo assembled transcriptome and a genome-based analysis using the available T. thymallus genome as a reference. The final assembly, with a contig N50 of 1221 and Benchmarking Universal Single-Copy Orthologs (BUSCO) scores above 94%, is made accessible along with structural and functional annotations and relative gene expression of the five tissues (NCBI SRA and FigShare databases). This is the first transcriptomic resource for this species, which provides a foundation for future research on this and other salmonid species that are increasingly exposed to environmental stressors.

Environmental RNA can distinguish life stages in amphibian populations

Applications of environmental DNA (eDNA) analysis methods for biomonitoring have grown exponentially over the last decade and provide a wealth of new information on the distribution of species. However, eDNA methods have limited application for estimating population-level metrics. Environmental RNA (eRNA) has the potential to address ecological questions by gathering population demographic information from environmental media but may be challenging to detect and analyze. We developed gene-specific eRNA assays targeting keratin-associated genes in two focal species, American bullfrogs (Lithobates catesbeianus) and tiger salamanders (Ambystoma mavortium) to answer an important question in amphibian management: whether species detections represent breeding populations versus transitory adults. We performed an extensive laboratory validation with amphibians housed across development stages, where we collected 95 and 127 environmental samples for bullfrogs and salamanders, respectively. Both assays were highly specific to the larval stage and amplified with high sensitivity (90% in bullfrog and 88.4% in tiger salamander samples). We then applied our validated assays to multiple natural systems. When larvae were present, we found 74.1% overall detection in bullfrog field samples and 70.8% and 48.5% overall detection in field samples from ponds with A. macrodactylum and A. californiense larvae, correlating with eDNA detection rates. When only adults were present, we did not detect larvae-specific eRNA in A. macrodactylum ponds, despite high eDNA detection rates. Although much work is ahead for optimizing assay design, sampling and filtering methods, we demonstrate that eRNA can successfully be used to discern life stages with direct application for ecology and conservation management.

A whole-body transcriptome assembly of the annelid worm Hediste diversicolor

The Annelida phylum is composed of a myriad of species exhibiting key phenotypic adaptations. They occupy key ecological niches in a variety of marine, freshwater and terrestrial ecosystems. Importantly, the increment of omic resources is rapidly modifying the taxonomic landscape and knowledge of species belonging to this phylum. Here, we comprehensively characterised and annotated a transcriptome of the common ragworm, Hediste diversicolor (OF Müller). This species belongs to the family Nereididae and inhabits estuarine and lagoon areas on the Atlantic coasts of Europe and North America. Ecologically, H. diversicolor plays an important role in benthic food webs. Given its commercial value, H. diversicolor is a promising candidate for aquaculture development and production in farming facilities, under a circular economy framework. We used Illumina next-generation sequencing technology, to produce a total of 105 million (M) paired-end (PE) raw reads and generate the first whole-body transcriptome assembly of H. diversicolor species. This high-quality transcriptome contains 69,335 transcripts with an N50 transcript length of 2313 bp and achieved a BUSCO gene completeness of 97.7% and 96% in Eukaryota and Metazoa lineage-specific profile libraries. Our findings offer a valuable resource for multiple biological applications using this species.

Navigating noisy waters: A review of field studies examining anthropogenic noise effects on wild fisha)

Anthropogenic noise is globally increasing in aquatic ecosystems, and there is concern that it may have adverse consequences in many fish species, yet the effects of noise in field settings are not well understood. Concern over the applicability of laboratory-conducted bioacoustic experiments has led to a call for, and a recent increase in, field-based studies, but the results have been mixed, perhaps due to the wide variety of techniques used and species studied. Previous reviews have explored the behavioral, physiological, and/or anatomical costs of fish exposed to anthropogenic noise, but few, if any, have focused on the field techniques and sound sources themselves. This review, therefore, aims to summarize, quantify, and interpret field-based literature, highlight novel approaches, and provide recommendations for future research into the effects of noise on fish.

Anemonefishes: A model system for evolutionary genomics

Anemonefishes are an iconic group of coral reef fish particularly known for their mutualistic relationship with sea anemones. This mutualism is especially intriguing as it likely prompted the rapid diversification of anemonefish. Understanding the genomic architecture underlying this process has indeed become one of the holy grails of evolutionary research in these fishes. Recently, anemonefishes have also been used as a model system to study the molecular basis of highly complex traits such as color patterning, social sex change, larval dispersal and life span. Extensive genomic resources including several high-quality reference genomes, a linkage map, and various genetic tools have indeed enabled the identification of genomic features controlling some of these fascinating attributes, but also provided insights into the molecular mechanisms underlying adaptive responses to changing environments. Here, we review the latest findings and new avenues of research that have led to this group of fish being regarded as a model for evolutionary genomics.

Environmental transcriptomics under heat stress: Can environmental RNA reveal changes in gene expression of aquatic organisms?

To safeguard biodiversity in a changing climate, taxonomic information about species turnover and insights into the health of organisms are required. Environmental DNA approaches are increasingly used for species identification, but cannot provide functional insights. Transcriptomic methods reveal the physiological states of macroorganisms, but are currently species-specific and require tissue sampling or animal sacrifice, making community-wide assessments challenging. Here, we test whether broad functional information (expression level of the transcribed genes) can be harnessed from environmental RNA (eRNA), which includes extra-organismal RNA from macroorganisms along with whole microorganisms. We exposed Daphnia pulex as well as phytoplankton prey and microorganism colonizers to control (20°C) and heat stress (28°C) conditions for 7 days. We sequenced eRNA from tank water (after complete removal of Daphnia) as well as RNA from Daphnia tissue, enabling comparisons of extra-organismal and organismal RNA-based gene expression profiles. Both RNA types detected similar heat stress responses of Daphnia. Using eRNA, we identified 32 Daphnia genes to be differentially expressed following heat stress. Of these, 17 were also differentially expressed and exhibited similar levels of relative expression in organismal RNA. In addition to the extra-organismal Daphnia response, eRNA detected community-wide heat stress responses consisting of distinct functional profiles and 121 differentially expressed genes across eight taxa. Our study demonstrates that environmental transcriptomics based on extra-organismal eRNA can noninvasively reveal gene expression responses of macroorganisms following environmental changes, with broad potential implications for the biomonitoring of health across the trophic chain.

Tissue-specific transcriptomes reveal potential mechanisms of microbiome heterogeneity in an ancient fish

The lake sturgeon (Acipenser fulvescens) is an ancient, octoploid fish faced with conservation challenges across its range in North America, but a lack of genomic resources has hindered molecular research in the species. To support such research, we created a transcriptomic database from 13 tissues: brain, esophagus, gill, head kidney, heart, white muscle, liver, glandular stomach, muscular stomach, anterior intestine, pyloric cecum, spiral valve and rectum. The transcriptomes for each tissue were sequenced and assembled individually from a mean of 98.3 million (±38.9 million SD) reads each. In addition, an overall transcriptome was assembled and annotated with all data used for each tissue-specific transcriptome. All assembled transcriptomes and their annotations were made publicly available as a scientific resource. The non-gut transcriptomes provide important resources for many research avenues. However, we focused our analysis on messenger ribonucleic acid (mRNA) observations in the gut because the gut represents a compartmentalized organ system with compartmentalized functions, and seven of the sequenced tissues were from each of these portions. These gut-specific analyses were used to probe evidence of microbiome regulation by studying heterogeneity in microbial genes and genera identified from mRNA annotations. Gene set enrichment analyses were used to reveal the presence of photoperiod and circadian-related transcripts in the pyloric cecum, which may support periodicity in lake sturgeon digestion. Similar analyses were used to identify different types of innate immune regulation across the gut, while analyses of unique transcripts annotated to microbes revealed heterogeneous genera and genes among different gut tissues. The present results provide a scientific resource and information about the mechanisms of compartmentalized function across gut tissues in a phylogenetically ancient vertebrate. Database URL: https://figshare.com/projects/Lake_Sturgeon_Transcriptomes/133143.

Harnessing the omics revolution to address the global biodiversity crisis

Human disturbances are altering global biodiversity in unprecedented ways. We identify three fundamental challenges underpinning our understanding of global biodiversity (namely discovery, loss, and preservation), and discuss how the omics revolution (e.g. genomics, transcriptomics, proteomics, metabolomics, and meta-omics) can help address these challenges. We also discuss how omics tools can illuminate the major drivers of biodiversity loss, including invasive species, pollution, urbanization, overexploitation, and climate change, with a special focus on highly diverse tropical environments. Although omics tools are transforming the traditional toolkit of biodiversity research, their application to addressing the current biodiversity crisis remains limited and may not suffice to offset current rates of biodiversity loss. Despite technical and logistical challenges, omics tools need to be fully integrated into global biodiversity research, and better strategies are needed to improve their translation into biodiversity policy and practice. It is also important to recognize that although the omics revolution can be considered the biologist's dream, socioeconomic disparity limits their application in biodiversity research.

Transcriptomic impacts and potential routes of detoxification in a lampricide-tolerant teleost exposed to TFM and niclosamide

Sea lamprey (Petromyzon marinus) control in the Laurentian Great Lakes of North America often relies on the application of 3-trifluoromethyl-4-nitrophenol (TFM) and niclosamide mixtures to kill larval sea lamprey. Selectivity of TFM against lampreys appears to be due to differential detoxification ability in these jawless fishes compared to bony fishes, particularly teleosts. However, the proximate mechanisms of tolerance to the TFM and niclosamide mixture and the mechanisms of niclosamide toxicity on its own are poorly understood, especially among non-target fishes. Here, we used RNA sequencing to identify specific mRNA transcripts and functional processes that responded to niclosamide or a TFM:niclosamide mixture in bluegill (Lepomis macrochirus). Bluegill were exposed to niclosamide or TFM:niclosamide mixture, along with a time-matched control group, and gill and liver tissues were sampled at 6, 12, and 24 h. We summarized the whole-transcriptome patterns through gene ontology (GO) term enrichment and through differential expression of detoxification genes. The niclosamide treatment resulted in an upregulation of several transcripts associated with detoxification (cyp, ugt, sult, gst), which may help explain the relatively high detoxification capacity in bluegill. Conversely, the TFM:niclosamide mixture resulted in an enrichment of processes related to arrested cell cycle and growth, and cell death alongside a diverse detoxification gene response. Detoxification of both lampricides likely involves the use of phase I and II biotransformation genes. Our findings strongly suggest that the unusually high tolerance of bluegill to lampricides is due to these animals having an inherently high capacity and flexible detoxification response to such compounds.

Influence of RNA-Seq library construction, sampling methods, and tissue harvesting time on gene expression estimation

RNA sequencing (RNA-Seq) is popular for measuring gene expression in non-model organisms, including wild populations. While RNA-Seq can detect gene expression variation among wild-caught individuals and yield important insights into biological function, sampling methods can also affect gene expression estimates. We examined the influence of multiple technical variables on estimated gene expression in a non-model fish, the westslope cutthroat trout (Oncorhynchus clarkii lewisi), using two RNA-Seq library types: 3' RNA-Seq (QuantSeq) and whole mRNA-Seq (NEB). We evaluated effects of dip netting versus electrofishing, and of harvesting tissue immediately versus 5 min after euthanasia on estimated gene expression in blood, gill, and muscle. We found no significant differences in gene expression between sampling methods or tissue collection times with either library type. When library types were compared using the same blood samples, 58% of genes detected by both NEB and QuantSeq showed significantly different expression between library types, and NEB detected 31% more genes than QuantSeq. Although the two library types recovered different numbers of genes and expression levels, results with NEB and QuantSeq were consistent in that neither library type showed differences in gene expression between sampling methods and tissue harvesting times. Our study suggests that researchers can safely rely on different fish sampling strategies in the field. In addition, while QuantSeq is more cost effective, NEB detects more expressed genes. Therefore, when it is crucial to detect as many genes as possible (especially low expressed genes), when alternative splicing is of interest, or when working with an organism lacking good genomic resources, whole mRNA-Seq is more powerful.

Blood-based gene expression as non-lethal tool for inferring salinity-habitat history of European eel (Anguilla anguilla)

The European eel is a facultative catadromous species, meaning that it can skip the freshwater phase or move between marine and freshwater habitats during its continental life stage. Otolith microchemistry, used to determine the habitat use of eel or its salinity history, requires the sacrifice of animals. In this context, blood-based gene expression may represent a non-lethal alternative. In this work, we tested the ability of blood transcriptional profiling to identify the different salinity-habitat histories of European eel. Eels collected from different locations in Norway were classified through otolith microchemistry as freshwater residents (FWR), seawater residents (SWR) or inter-habitat shifters (IHS). We detected 3451 differentially expressed genes from blood by comparing FWR and SWR groups, and then used that subset of genes in a machine learning approach (i.e., random forest) to the extended FWR, SWR, and IHS group. Random forest correctly classified 100% of FWR and SWR and 83% of the IHS using a minimum of 30 genes. The implementation of this non-lethal approach may replace otolith-based microchemistry analysis for the general assessment of life-history tactics in European eels. Overall, this approach is promising for the replacement or reduction of other lethal analyses in determining certain fish traits.

aTAP: automated transcriptome analysis platform for processing RNA-seq data by de novo assembly

RNA-seq is a sequencing technique that uses next-generation sequencing (NGS) to explore and study the entire transcriptome of a biological sample. NGS-based analyses are mostly performed via command-line interfaces, which is an obstacle for molecular biologists and researchers. Therefore, the higher throughputs from NGS can only be accessed with the help of bioinformatics and computer science expertise. As the cost of sequencing is continuously falling, the use of RNA-seq seems certain to increase. To minimize the problems encountered by biologists and researchers in RNA-seq data analysis, we propose an automated platform with a web application that integrates various bioinformatics pipelines. The platform is intended to enable academic users to more easily analyze transcriptome datasets. Our automated Transcriptome Analysis Platform (aTAP) offers comprehensive bioinformatics workflows, including quality control of raw reads, trimming of low-quality reads, de novo transcriptome assembly, transcript expression quantification, differential expression analysis, and transcript annotation. aTAP has a user-friendly graphical interface, allowing researchers to interact with and visualize results in the web browser. This project offers an alternative way to analyze transcriptome data, by integrating efficient and well-known tools, that is simpler and more accessible to research communities. aTAP is freely available to academic users at https://atap.psu.ac.th/.

Acute thermal stress elicits interactions between gene expression and alternative splicing in a fish of conservation concern

Transcriptomic research provides a mechanistic understanding of an organism's response to environmental challenges such as increasing temperatures, which can provide key insights into the threats posed by thermal challenges associated with urbanization and climate change. Differential gene expression and alternative splicing are two elements of the transcriptomic stress response that may work in tandem, but relatively few studies have investigated these interactions in fishes of conservation concern. We studied the imperilled redside dace (Clinostomus elongatus) as thermal stress is hypothesised to be an important cause of population declines. We tested the hypothesis that gene expression-splicing interactions contribute to the thermal stress response. Wild fish exposed to acute thermal stress were compared with both handling controls and fish sampled directly from a river. Liver tissue was sampled to study the transcriptomic stress response. With a gene set enrichment analysis, we found that thermally stressed fish showed a transcriptional response related to transcription regulation and responses to unfolded proteins, and alternatively spliced genes related to gene expression regulation and metabolism. One splicing factor, prpf38b, was upregulated in the thermally stressed group compared to the other treatments. This splicing factor may have a role in the Jun/AP-1 cellular stress response, a pathway with wide-ranging and context-dependent effects. Given large gene interaction networks and the context-dependent nature of transcriptional responses, our results highlight the importance of understanding interactions between gene expression and splicing for understanding transcriptomic responses to thermal stress. Our results also reveal transcriptional pathways that can inform conservation breeding, translocation, and reintroduction programs for redside dace and other imperilled species by identifying appropriate source populations.

Transcriptome profiles revealed high- and low-salinity water altered gill homeostasis in half-smooth tongue sole (Cynoglossus semilaevis)

Salinity is an important environmental factor that affects fish growth, development, and reproduction. As euryhaline fish, half-smooth tongue sole (Cynoglossus semilaevis) are a suitable species for deciphering the salinity adaptation mechanism of fish; however, the molecular mechanisms underlying low- and high-salinity responses remain unclear. In this study, RNA-seq was applied to characterize the genes and regulatory pathways involved in C. semilaevis gill responses to high- (32 ppt), low- (8 ppt), and control-salinity (24 ppt) water. Gills were rich in mitochondria-rich cells (MRCs) in high salinity. Compared with control, 2137 and 218 differentially expressed genes (DEGs) were identified in low and high salinity, respectively. The enriched functions of most DEGs were metabolism, ion transport, regulation of cell cycle, and immune response. The DEGs involved in oxidative phosphorylation, citrate cycle, and fatty acid metabolism were down-regulated in low salinity. For ion transport, high and low salinity significantly altered the expressions of prlr, ca12, and cftr. In cell cycle arrest and cellular repair, gadd45b, igfbp5, and igfbp2 were significantly upregulated in high and low salinity. For immune response, il10, il34, il12b, and crp increased in high and low salinity. Our findings suggested that alterations in material and energy metabolism, ions transport, cell cycle arrest, cellular repair, and immune response, are required to maintain C. semilaevis gill homeostasis under high and low salinity. This study provides insight into the divergence of C. semilaevis osmoregulation mechanisms acclimating to high and low salinity, which will serve as reference for the healthy culture of C. semilaevis.

Non-Lethal Sampling Supports Integrative Movement Research in Freshwater Fish

Freshwater ecosystems and fishes are enormous resources for human uses and biodiversity worldwide. However, anthropogenic climate change and factors such as dams and environmental contaminants threaten these freshwater systems. One way that researchers can address conservation issues in freshwater fishes is via integrative non-lethal movement research. We review different methods for studying movement, such as with acoustic telemetry. Methods for connecting movement and physiology are then reviewed, by using non-lethal tissue biopsies to assay environmental contaminants, isotope composition, protein metabolism, and gene expression. Methods for connecting movement and genetics are reviewed as well, such as by using population genetics or quantitative genetics and genome-wide association studies. We present further considerations for collecting molecular data, the ethical foundations of non-lethal sampling, integrative approaches to research, and management decisions. Ultimately, we argue that non-lethal sampling is effective for conducting integrative, movement-oriented research in freshwater fishes. This research has the potential for addressing critical issues in freshwater systems in the future.

Understanding risks and consequences of pathogen infections on the physiological performance of outmigrating Chinook salmon

The greatest concentration of at-risk anadromous salmonids is found in California (USA)-the populations that have been negatively impacted by the degradation of freshwater ecosystems. While climate-driven environmental changes threaten salmonids directly, they also change the life cycle dynamics and geographic distribution of pathogens, their resulting host-pathogen interactions and potential for disease progression. Recent studies have established the correlation between pathogen detection and salmonid smolt mortality during their migration to the ocean. The objective of the present study was to screen for up to 47 pathogens in juvenile Chinook salmon (Oncorhynchus tshawytscha) that were held in cages at two key sites of the Sacramento River (CA, USA) and measure potential consequences on fish health. To do so, we used a combination of transcriptomic analysis, enzymatic assays for energy metabolism and hypoxia and thermal tolerance measures. Results revealed that fish were infected by two myxozoan parasites: Ceratonova shasta and Parvicapsula minibicornis within a 2-week deployment. Compared to the control fish maintained in our rearing facility, infected fish displayed reduced body mass, depleted hepatic glycogen stores and differential regulation of genes involved in the immune and general stress responses. This suggests that infected fish would have lower chances of migration success. In contrast, hypoxia and upper thermal tolerances were not affected by infection, suggesting that infection did not impair their capacity to cope with acute abiotic stressors tested in this study. An evaluation of long-term consequences of the observed reduced body mass and hepatic glycogen depletion is needed to establish a causal relationship between salmon parasitic infection and their migration success. This study highlights that to assess the potential sublethal effects of a stressor, or to determine a suitable management action for fish, studies need to consider a combination of endpoints from the molecular to the organismal level.

Host gene expression in wildlife disease: making sense of species-level responses

Emerging infectious diseases are significant threats to wildlife conservation, yet the impacts of pathogen exposure and infection can vary widely among host species. As such, conservation biologists and disease ecologists have increasingly aimed to understand species-specific host susceptibility using molecular methods. In particular, comparative gene expression assays have been used to contrast the transcriptomic responses of disease-resistant and disease-susceptible hosts to pathogen exposure. This work usually assumes that the gene expression responses of disease-resistant species will reveal the activation of molecular pathways contributing to host defence. However, results often show that disease-resistant hosts undergo little gene expression change following pathogen challenge. Here, we discuss the mechanistic implications of these "null" findings and offer methodological suggestions for future molecular studies of wildlife disease. First, we highlight that muted transcriptomic responses with minimal immune system recruitment may indeed be protective for nonsusceptible hosts if they limit immunopathology and promote pathogen tolerance in systems where susceptible hosts suffer from genetic dysregulation. Second, we argue that overly narrow investigation of responses to pathogen exposure may overlook important, constitutively active molecular pathways that underlie species-specific defences. Finally, we outline alternative study designs and approaches that complement interspecific transcriptomic comparisons, including intraspecific gene expression studies and genomic methods to detect signatures of selection. Collectively, these insights will help ecologists extract maximal information from conservation-relevant transcriptomic data sets, leading to a deeper understanding of host defences and, ultimately, the implementation of successful conservation interventions.

Molecular and physiological responses predict acclimation limits in juvenile brook trout (Salvelinus fontinalis)

Understanding the resilience of ectotherms to high temperatures is essential because of the influence of climate change on aquatic ecosystems. The ability of species to acclimate to high temperatures may determine whether populations can persist in their native ranges. We examined physiological and molecular responses of juvenile brook trout (Salvelinus fontinalis) to six acclimation temperatures (5, 10, 15, 20, 23 and 25°C) that span the thermal distribution of the species to predict acclimation limits. Brook trout exhibited an upregulation of stress-related mRNA transcripts (heat shock protein 90-beta, heat shock cognate 71 kDa protein, glutathione peroxidase 1) and downregulation of transcription factors and osmoregulation-related transcripts (nuclear protein 1, Na+/K+/2Cl- co-transporter-1-a) at temperatures ≥20°C. We then examined the effects of acclimation temperature on metabolic rate (MR) and physiological parameters in fish exposed to an acute exhaustive exercise and air exposure stress. Fish acclimated to temperatures ≥20°C exhibited elevated plasma cortisol and glucose, and muscle lactate after exposure to the acute stress. Fish exhibited longer MR recovery times at 15 and 20°C compared with the 5 and 10°C groups; however, cortisol levels remained elevated at temperatures ≥20°C after 24 h. Oxygen consumption in fish acclimated to 23°C recovered quickest after exposure to acute stress. Standard MR was highest and factorial aerobic scope was lowest for fish held at temperatures ≥20°C. Our findings demonstrate how molecular and physiological responses predict acclimation limits in a freshwater fish as the brook trout in the present study had a limited ability to acclimate to temperatures beyond 20°C.

Nonlinear effects of environmental salinity on the gill transcriptome versus proteome of Oreochromis niloticus modulate epithelial cell turnover

A data-independent acquisition (DIA) assay library for targeted quantitation of thousands of Oreochromis niloticus gill proteins using a label- and gel-free workflow was generated and used to compare protein and mRNA abundances. This approach generated complimentary rather than redundant data for 1899 unique genes in gills of tilapia acclimated to freshwater and brackish water. Functional enrichment analyses identified mitochondrial energy metabolism, serine protease and immunity-related functions, and cytoskeleton/ extracellular matrix organization as major processes controlled by salinity in O. niloticus gills. Non-linearity in salinity-dependent transcriptome versus proteome regulation was revealed for specific functional groups of genes. The relationship was more linear for other molecular functions/ cellular processes, suggesting that the salinity-dependent regulation of O. niloticus gill function relies on post-transcriptional mechanisms for some functions/ processes more than others. This integrative systems biology approach can be adopted for other tissues and organisms to study cellular dynamics for many changing ecological contexts.

Conservation insights from wild bee genetic studies: Geographic differences, susceptibility to inbreeding, and signs of local adaptation

Conserving bees are critical both ecologically and economically. Genetic tools are valuable for monitoring these vital pollinators since tracking these small, fast-flying insects by traditional means is difficult. By surveying the current state of the literature, this review discusses how recent advances in landscape genetic and genomic research are elucidating how wild bees respond to anthropogenic threats. Current literature suggests that there may be geographic differences in the vulnerability of bee species to landscape changes. Populations of temperate bee species are becoming more isolated and more genetically depauperate as their landscape becomes more fragmented, but tropical bee species appear unaffected. These differences may be an artifact of historical differences in land-use, or it suggests that different management plans are needed for temperate and tropical bee species. Encouragingly, genetic studies on invasive bee species indicate that low levels of genetic diversity may not lead to rapid extinction in bees as once predicted. Additionally, next-generation sequencing has given researchers the power to identify potential genes under selection, which are likely critical to species' survival in their rapidly changing environment. While genetic studies provide insights into wild bee biology, more studies focusing on a greater phylogenetic and life-history breadth of species are needed. Therefore, caution should be taken when making broad conservation decisions based on the currently few species examined.

The use of non-lethal sampling for transcriptomics to assess the physiological status of wild fishes

Fishes respond to different abiotic and biotic stressors through changes in gene expression as a part of an integrated physiological response. Transcriptomics approaches have been used to quantify gene expression patterns as a reductionist approach to understand responses to environmental stressors in animal physiology and have become more commonly used to study wild fishes. We argue that non-lethal sampling for transcriptomics should become the norm for assessing the physiological status of wild fishes, especially when there are conservation implications. Processes at the level of the transcriptome provide a "snapshot" of the cellular conditions at a given time; however, by using a non-lethal sampling protocol, researchers can connect the transcriptome profile with fitness-relevant ecological endpoints such as reproduction, movement patterns and survival. Furthermore, telemetry is a widely used approach in fisheries to understand movement patterns in the wild, and when combined with transcriptional profiling, provides arguably the most powerful use of non-lethal sampling for transcriptomics in wild fishes. In this review, we discuss the different tissues that can be successfully incorporated into non-lethal sampling strategies, which is particularly useful in the context of the emerging field of conservation transcriptomics. We briefly describe different methods for transcriptional profiling in fishes from high-throughput qPCR to whole transcriptome approaches. Further, we discuss strategies and the limitations of using transcriptomics for non-lethally studying fishes. Lastly, as 'omics' technology continues to advance, transcriptomics paired with different omics approaches to study wild fishes will provide insight into the factors that regulate phenotypic variation and the physiological responses to changing environmental conditions in the future.

RNA sequencing describes both population structure and plasticity-selection dynamics in a non-model fish

Background

Messenger RNA sequencing is becoming more common in studies of non-model species and is most often used for gene expression-based investigations. However, the method holds potential for numerous other applications as well—including analyses of alternative splicing, population structure, and signatures of selection. To maximize the utility of mRNA data sets, distinct analyses may be combined such as by exploring dynamics between gene expression with signatures of selection in the context of population structure. Here, we compare two published data sets describing two populations of a minnow species endemic to the San Francisco Estuary (Sacramento splittail, Pogonichthys macrolepidotus): a microsatellite data set showing population structure, and an mRNA whole transcriptome data set obtained after the two populations were exposed to a salinity challenge. We compared measures of population structure and genetic variation using single nucleotide polymorphisms (SNPs) called from mRNA from the whole transcriptome sequencing study with those patterns determined from microsatellites. For investigating plasticity and evolution, intra- and inter-population transcriptome plasticity was investigated with differential gene expression, differential exon usage, and gene expression variation. Outlier SNP analysis was also performed on the mRNA data set and signatures of selection and phenotypic plasticity were investigated on an individual-gene basis.

Results

We found that mRNA sequencing revealed patterns of population structure consistent with those found with microsatellites, but with lower magnitudes of genetic variation and population differentiation consistent with widespread purifying selection expected when using mRNA. In addition, within individual genes, phenotypic plasticity or signatures of selection were found in almost mutual exclusion (except heatr6, nfu1, slc22a6, sya, and mmp13).

Conclusions

These results show that an mRNA sequencing data set may have multiple uses, including describing population structure and for investigating the mechanistic interplay of evolution and plasticity in adaptation. MRNA sequencing thus complements traditional sequencing methods used for population genetics, in addition to its utility for describing phenotypic plasticity.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07592-4.

Infected juvenile salmon can experience increased predation during freshwater migration

Predation risk for animal migrants can be impacted by physical condition. Although size- or condition-based selection is often observed, observing infection-based predation is rare due to the difficulties in assessing infectious agents in predated samples. We examined predation of outmigrating sockeye salmon (Oncorhynchus nerka) smolts by bull trout (Salvelinus confluentus) in south-central British Columbia, Canada. We used a high-throughput quantitative polymerase chain reaction (qPCR) platform to screen for the presence of 17 infectious agents found in salmon and assess 14 host genes associated with viral responses. In one (2014) of the two years assessed (2014 and 2015), the presence of infectious haematopoietic necrosis virus (IHNv) resulted in 15-26 times greater chance of predation; in 2015 IHNv was absent among all samples, predated or not. Thus, we provide further evidence that infection can impact predation risk in migrants. Some smolts with high IHNv loads also exhibited gene expression profiles consistent with a virus-induced disease state. Nine other infectious agents were observed between the two years, none of which were associated with increased selection by bull trout. In 2014, richness of infectious agents was also associated with greater predation risk. This is a rare demonstration of predator consumption resulting in selection for prey that carry infectious agents. The mechanism by which this selection occurs is not yet determined. By culling infectious agents from migrant populations, fish predators could provide an ecological benefit to prey.

Applying a gene-suite approach to examine the physiological status of wild-caught walleye (Sander vitreus)

Molecular techniques have been increasingly used in a conservation physiology framework to provide valuable information regarding the mechanisms underlying responses of wild organisms to environmental and anthropogenic stressors. In the present study, we developed a reference gill transcriptome for walleye (Sander vitreus), allowing us to pair a gene-suite approach (i.e. multiple genes across multiple cellular processes) with multivariate statistics to examine the physiological status of wild-caught walleye. For molecular analyses of wild fish, the gill is a useful target for conservation studies, not only because of its importance as an indicator of the physiological status of fish but also because it can be biopsied non-lethally. Walleye were non-lethally sampled following short- (~1.5 months) and long-term (~3.5 months) confinement in the Delta Marsh, which is located south of Lake Manitoba in Manitoba, Canada. Large-bodied walleye are confined in the Delta Marsh from late April to early August by exclusion screens used to protect the marsh from invasive common carp (Cyprinus carpio), exposing fish to potentially stressful water quality conditions. Principal components analysis revealed patterns of transcript abundance consistent with exposure of fish to increasingly high temperature and low oxygen conditions with longer holding in the marsh. For example, longer-term confinement in the marsh was associated with increases in the mRNA levels of heat shock proteins and a shift in the mRNA abundance of aerobic to anaerobic metabolic genes. Overall, the results of the present study suggest that walleye confined in the Delta Marsh may be exhibiting sub-lethal responses to high temperature and low oxygen conditions. These results provide valuable information for managers invested in mediating impacts to a local species of conservation concern. More broadly, we highlight the usefulness of pairing transcriptomic techniques with multivariate statistics to address potential confounding factors that can affect measured physiological responses of wild-caught fish.

Omics approaches for conservation biology research on the bivalve Chamelea gallina

The striped venus (Chamelea gallina) is an important economic resource in the Mediterranean Basin; this species has exhibited a strong quantitative decline in the Adriatic Sea. The aim of this work was to provide a comprehensive view of the biological status of C. gallina to elucidate the bioecological characteristics and genetic diversity of wild populations. To the best of our knowledge, this investigation is the first to perform a multidisciplinary study on C. gallina based on two omics approaches integrated with histological, ecotoxicological, and chemical analyses and with the assessment of environmental parameters. The results obtained through RNA sequencing indicated that the striped venus has a notable ability to adapt to different environmental conditions. Moreover, the stock reduction exhibited by this species in the last 2 decades seems not to have negatively affected its genetic diversity. Indeed, the high level of genetic diversity that emerged from our ddRAD dataset analyses is ascribable to the high larval dispersal rate, which might have played a “compensatory role” on local fluctuations, conferring to this species a good adaptive potential to face the environmental perturbations. These findings may facilitate the efforts of conservation biologists to adopt ad hoc management plans for this fishery resource.

Transcriptomic response to elevated water temperatures in adult migrating Yukon River Chinook salmon (Oncorhynchus tshawytscha)

Chinook salmon (Oncorhynchus tshawytscha) declines are widespread and may be attributed, at least in part, to warming river temperatures. Water temperatures in the Yukon River and tributaries often exceed 18°C, a threshold commonly associated with heat stress and elevated mortality in Pacific salmon. Untangling the complex web of direct and indirect physiological effects of heat stress on salmon is difficult in a natural setting with innumerable system challenges but is necessary to increase our understanding of both lethal and sublethal impacts of heat stress on populations. The goal of this study was to characterize the cellular stress response in multiple Chinook salmon tissues after acute elevated temperature challenges. We conducted a controlled 4-hour temperature exposure (control, 18°C and 21°C) experiment on the bank of the Yukon River followed by gene expression (GE) profiling using a 3'-Tag-RNA-Seq protocol. The full transcriptome was analysed for 22 Chinook salmon in muscle, gill and liver tissue. Both the 21°C and 18°C treatments induced greater activity in genes associated with protein folding (e.g. HSP70, HSP90 mRNA) processes in all tissues. Global GE patterns indicate that transcriptomic responses to heat stress were highly tissue-specific, underscoring the importance of analyzing multiple tissues for determination of physiological effect. Primary superclusters (i.e. groupings of loosely related terms) of altered biological processes were identified in each tissue type, including regulation of DNA damage response (gill), regulation by host of viral transcription (liver) and regulation of the force of heart contraction (muscle) in the 21°C treatment. This study provides insight into mechanisms potentially affecting adult Chinook salmon as they encounter warm water during their spawning migration in the Yukon River and suggests that both basic and more specialized cellular functions may be disrupted.

Genomic signals found using RNA sequencing show signatures of selection and subtle population differentiation in walleye (Sander vitreus) in a large freshwater ecosystem

RNA sequencing is an effective approach for studying aquatic species yielding both physiological and genomic data. However, its population genetic applications are not well-characterized. We investigate this possible role for RNA sequencing for population genomics in Lake Winnipeg, Manitoba, Canada, walleye (Sander vitreus). Lake Winnipeg walleye represent the largest component of the second-largest freshwater fishery in Canada. In the present study, large female walleye were sampled via nonlethal gill biopsy over two years at three spawning sites representing a latitudinal gradient in the lake. Genetic variation from sequenced mRNA was analyzed for neutral and adaptive markers to investigate population structure and possible adaptive variation. We find low population divergence (F ST = 0.0095), possible northward gene flow, and outlier loci that vary latitudinally in transcripts associated with cell membrane proteins and cytoskeletal function. These results indicate that Lake Winnipeg walleye may be effectively managed as a single demographically connected metapopulation with contributing subpopulations and suggest genomic differences possibly underlying observed phenotypic differences. Despite its high cost relative to other genotyping methods, RNA sequencing data can yield physiological in addition to genetic information discussed here. We therefore argue that it is useful for addressing diverse molecular questions in the conservation of freshwater species.

Causes and consequences of intraspecific variation in animal responses to anthropogenic noise

Anthropogenic noise is a recognized global pollutant, affecting a wide range of nonhuman animals. However, most research considers only whether noise pollution has an impact, ignoring that individuals within a species or population exhibit substantial variation in responses to stress. Here, we first outline how intrinsic characteristics (e.g., body size, condition, sex, and personality) and extrinsic factors (e.g., environmental context, repeated exposure, prior experience, and multiple stressors) can affect responses to environmental stressors. We then present the results of a systematic search of the anthropogenic-noise literature, identifying articles that investigated intraspecific variation in the responses of nonhuman animals to noise. This reveals that fewer than 10% of articles (51 of 589) examining impacts of noise test experimentally for intraspecific variation in responses; of those that do, more than 75% report significant effects. We assess these existing studies to determine the current scope of research and findings to-date, and to provide suggestions for good practice in the design, implementation, and reporting of robust experiments in this field. We close by explaining how understanding intraspecific variation in responses to anthropogenic noise is crucial for improving how we manage captive animals, monitor wild populations, model species responses, and mitigate effects of noise pollution on wildlife. Our aim is to stimulate greater knowledge and more effective management of the harmful consequences of this global pollutant.

Differential regulation of select osmoregulatory genes and Na+/K+-ATPase paralogs may contribute to population differences in salinity tolerance in a semi-anadromous fish

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.

Transcriptomics of Tasmanian Devil (Sarcophilus Harrisii) Ear Tissue Reveals Homogeneous Gene Expression Patterns across a Heterogeneous Landscape

In an era of unprecedented global change, exploring patterns of gene expression among wild populations across their geographic range is crucial for characterizing adaptive potential. RNA-sequencing studies have successfully characterized gene expression differences among populations experiencing divergent environmental conditions in a wide variety of taxa. However, few of these studies have identified transcriptomic signatures to multivariate, environmental stimuli among populations in their natural environments. Herein, we aim to identify environmental and sex-driven patterns of gene expression in the Tasmanian devil (Sarcophilus harrisii), a critically endangered species that occupies a heterogeneous environment. We performed RNA-sequencing on ear tissue biopsies from adult male and female devils from three populations at the extremes of their geographic range. There were no transcriptome-wide patterns of differential gene expression that would be suggestive of significant, environmentally-driven transcriptomic responses. The general lack of transcriptome-wide variation in gene expression levels across the devil's geographic range is consistent with previous studies that documented low levels of genetic variation in the species. However, genes previously implicated in local adaptation to abiotic environment in devils were enriched for differentially expressed genes. Additionally, three modules of co-expressed genes were significantly associated with either population of origin or sex.

Physiological status of silver carp (Hypophthalmichthys molitrix) in the Illinois River: An assessment of fish at the leading edge of the invasion front

Silver carp (Hypophthalmichthys molitrix) are invasive to North America, and their range has expanded within the Mississippi River Basin, seemingly unchecked, since their introduction in the late 1970s, with the exception of the upper reaches of the Illinois River. With the imminent threat of their movement into the Great Lakes, the goal of the present study was to assess whether differences in the physiological status between silver carp at the leading edge of their invasion front and core population sites could explain their lack of expansion upstream toward Lake Michigan over the past decade. A transcriptomic approach using RNA sequencing and analysis of plasma variables were used to quantify differences among fish at the leading edge and two downstream core population sites. Leading-edge fish exhibited upregulation of genes associated with xenobiotic defense (e.g., ATP-binding cassette C1 [abcc1], abcc2, abcc6), decreased cell integrity (i.e., macroautophagy and apoptosis; autophagy-related protein 9A [atg9a], caspase 3b [casp3b]), and cholesterol metabolism (e.g., abca1, apolipoprotein A1 [apoa1], sterol O-acyltransferase [soat1]) and downregulation of genes associated with DNA repair (e.g., tumor suppressor p53-binding protein 1 [tp53bp1]) compared to core population sites. Transcriptomic profiles of leading-edge fish were consistent with fish inhabiting a polluted environment and suggest that poorer water quality conditions upstream of the leading edge may represent a non-permanent barrier to silver carp range expansion. The present study provides potential molecular targets for monitoring the physiological status of silver carp over time and in response to future improvements in water quality upstream of their leading edge.

Salmonid gene expression biomarkers indicative of physiological responses to changes in salinity and temperature, but not dissolved oxygen

An organism's ability to respond effectively to environmental change is critical to its survival. Yet, life stage and overall condition can dictate tolerance thresholds to heightened environmental stressors, such that stress may not be equally felt across individuals and at all times. Also, the transcriptional responses induced by environmental changes can reflect both generalized responses as well as others that are highly specific to the type of change being experienced. Thus, if transcriptional biomarkers specific to a stressor, even under multi-stressor conditions, can be identified, the biomarkers could then be applied in natural environments to determine when and where an individual experiences such a stressor. Here, we experimentally challenged juvenile Chinook salmon (Oncorhynchus tshawytscha) to validate candidate gill gene expression biomarkers. A sophisticated experimental design manipulated salinity (freshwater, brackish water and seawater), temperature (10, 14 and 18°C) and dissolved oxygen (normoxia and hypoxia) in all 18 possible combinations for 6 days using separate trials for three smolt statuses (pre-smolt, smolt and de-smolt). In addition, changes in juvenile behaviour, plasma variables, gill Na+/K+-ATPase activity, body size, body morphology and skin pigmentation supplemented the gene expression responses. We identified biomarkers specific to salinity and temperature that transcended the multiple stressors, smolt status and mortality (live, dead and moribund). Similar biomarkers for dissolved oxygen were not identified. This work demonstrates the unique power of gene expression biomarkers to identify a specific stressor even under multi-stressor conditions, and we discuss our next steps for hypoxia biomarkers using an RNA-seq study.

Development of SNVs and indels markers mined out of the first multi-organ transcriptomes from Hypancistrus zebra (Loricariidae), an endangered Amazonian catfish

Hypancistrus zebra is a catfish, endemic from the Xingu River, threatened with extinction due to the impacts of Belo Monte dam, of its illegal capture, of gold mining activities and of climate change. Currently, there are three nucleotide sequences from this species in GenBank, what impedes the development of genetic markers to assist on its conservation. A total of 217 million RNA-Seq reads from seven organs were sequenced and used to assemble 566,607 transcripts, including 98% of BUSCO vertebrates orthologs, 11,321 transcripts with SNVs and 1,724 transcripts with indels. Three transcripts with SNVs and five transcripts with indels were validated as the best candidate markers to conservation practices. This work illustrates the use of transcriptomics in conservation, by the development of a bigger toolbox for an endangered fish, and shall further contribute to studies on this and others related species reproduction, physiology, and adaptability to environmental changes.

Divergent transcriptomic signatures in response to salinity exposure in two populations of an estuarine fish

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.

Developing specific molecular biomarkers for thermal stress in salmonids

Background

Pacific salmon (Oncorhynchus spp.) serve as good biological indicators of the breadth of climate warming effects on fish because their anadromous life cycle exposes them to environmental challenges in both marine and freshwater environments. Our study sought to mine the extensive functional genomic studies in fishes to identify robust thermally-responsive biomarkers that could monitor molecular physiological signatures of chronic thermal stress in fish using non-lethal sampling of gill tissue.

Results

Candidate thermal stress biomarkers for gill tissue were identified using comparisons among microarray datasets produced in the Molecular Genetics Laboratory, Pacific Biological Station, Nanaimo, BC, six external, published microarray studies on chronic and acute temperature stress in salmon, and a comparison of significant genes across published studies in multiple fishes using deep literature mining. Eighty-two microarray features related to 39 unique gene IDs were selected as candidate chronic thermal stress biomarkers. Most of these genes were identified both in the meta-analysis of salmon microarray data and in the literature mining for thermal stress markers in salmonids and other fishes. Quantitative reverse transcription PCR (qRT-PCR) assays for 32 unique genes with good efficiencies across salmon species were developed, and their activity in response to thermally challenged sockeye salmon (O. nerka) and Chinook salmon (O. tshawytscha) (cool, 13–14 °C and warm temperatures 18–19 °C) over 5–7 days was assessed. Eight genes, including two transcripts of each SERPINH1 and HSP90AA1, FKBP10, MAP3K14, SFRS2, and EEF2 showed strong and robust chronic temperature stress response consistently in the discovery analysis and both sockeye and Chinook salmon validation studies.

Conclusions

The results of both discovery analysis and gene expression showed that a panel of genes involved in chaperoning and protein rescue, oxidative stress, and protein biosynthesis were differentially activated in gill tissue of Pacific salmon in response to elevated temperatures. While individually, some of these biomarkers may also respond to other stressors or biological processes, when expressed in concert, we argue that a biomarker panel comprised of some or all of these genes could provide a reliable means to specifically detect thermal stress in field-caught salmon.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-5108-9) contains supplementary material, which is available to authorized users.

Transcriptomic, lipid, and histological profiles suggest changes in health in fish from a pesticide hot spot

Barramundi (Lates calcarifer) were collected at the beginning (1st sampling) and end (2nd sampling) of the wet season from Sandy Creek, an agriculturally impacted catchment in the Mackay Whitsundays region of the Great Barrier Reef catchment area, and from Repulse Creek, located approximately 100 km north in Conway National Park, to assess the impacts of pesticide exposure. Gill and liver histology, lipid class composition in muscle, and the hepatic transcriptome were examined. The first sample of Repulse Creek fish showed little tissue damage and low transcript levels of xenobiotic metabolism enzymes. Sandy Creek fish showed altered transcriptomic patterns, including those that regulate lipid metabolism, xenobiotic metabolism, and immune response; gross histological alterations including lipidosis; and differences in some lipid classes. The second sampling of Repulse Creek fish showed similar alterations in hepatic transcriptome and tissue structure as fish from Sandy Creek. These changes may indicate a decrease in health of pesticide exposed fish.

Contaminant exposure effects in a changing climate: how multiple stressors can multiply exposure effects in the amphipod Hyalella azteca

Global climate change (GCC) is likely to intensify the synergistic effects between altered physicochemical parameters [of changing habitats] and other anthropogenic threats, such as water pollution, posing increased risks to aquatic biodiversity. As such, it is critical to understand how organisms will respond to changes in water temperature and salinity in the presence of contaminants. We exposed the epibenthic amphipod Hyalella azteca to a 3 × 3 factorial treatment design of three temperatures and three salinities ranging from 12 to 18 °C and 0 to 8 parts per thousand (ppt), respectively, in combination with a low-level environmentally relevant concentration of the pyrethroid insecticide bifenthrin (1 ng/L). Effects on survival and swimming behavior were evaluated after 96 h exposure. Transcription of a select suite of genes was monitored at 24, 48, and 96 h using quantitative polymerase chain reaction (qPCR). Our results not only demonstrate that the changes in salinity and temperature result in negative effects to invertebrate survival, behavior, and gene response, but that the effects were significantly more pronounced in the presence of bifenthrin. This is particularly important since greater thermal fluctuations, changes in timing and extent of glacial melt, and changes in precipitation, could result in H. azteca experiencing lower temperatures at times that coincide with increased spraying of pyrethroids. These environmentally relevant exposures using the standard test species H. azteca provide essential information for understanding effects caused by GCC in conjunction with increasing pesticide use, further highlighting the need to incorporate GCC impacts into risk assessments of contaminants of concern.