Transcriptomic responses of Antarctic clam Laternula elliptica to nanoparticles, at single and combined exposures reveal ecologically relevant biomarkers

In recent years micro- and nanoplastics and metal-oxide nanomaterials have been found in several environmental compartments. The Antarctic soft clam Laternula elliptica is an endemic Antarctic species having a wide distribution in the Southern Ocean. Being a filter-feeder, it could act as suitable bioindicator of pollution from nanoparticles also considering its sensitivity to various sources of stress. The present study aims to assess the impact of polystyrene nanoparticles (PS-NP) and the nanometal titanium-dioxide (n-TiO2) on genome-wide transcript expression of L. elliptica either alone and in combination and at two toxicological relevant concentrations (5 and 50 µg/L) during 96 h exposure. Transcript-target qRT-PCR was performed with the aim to identify suitable biomarkers of exposure and effects. As expected, at the highest concentration tested, the clustering was clearer between control and exposed clams. A total of 221 genes resulted differentially expressed in exposed clams and control ones, and 21 of them had functional annotation such as ribosomal proteins, antioxidant, ion transport (osmoregulation), acid-base balance, immunity, lipid metabolism, cell adhesion, cytoskeleton, apoptosis, chromatin condensation and cell signaling. At functional level, relevant transcripts were shared among some treatments and could be considered as general stress due to nanoparticle exposure. After applying transcript-target approach duplicating the number of clam samples, four ecologically relevant transcripts were revealed as biomarkers for PS-NP, n-TiO2 and their combination at 50 µg/L, that could be used for monitoring clams' health status in different Antarctic localities.

Non-concordant epigenetic and transcriptional responses to acute thermal stress in western mosquitofish (Gambusia affinis)

Climate change is intensifying the frequency and severity of extreme temperatures. Understanding the molecular mechanisms underlying the ability to cope with acute thermal stress is key for predicting species' responses to extreme temperature events. While many studies have focused on the individual roles of gene expression, post-transcriptional processes and epigenetic modifications in response to acute thermal stress, the relative contribution of these molecular mechanisms remains unclear. The wide range of thermal limits of western mosquitofish (Gambusia affinis) provides an opportunity to explore this interplay. Here, we quantified changes in gene expression, alternative splicing, DNA methylation and microRNA (miRNA) expression in muscle tissue dissected from mosquitofish immediately after reaching high (CTmax) or low thermal limit (CTmin). Although the numbers of genes showing expression and splicing changes in response to acute temperature stress were small, we found a possibly larger and non-redundant role of splicing compared to gene expression, with more genes being differentially spliced (DSGs) than differentially expressed (DEGs), and little overlap between DSGs and DEGs. We also identified a small proportion of CpGs showing significant methylation change (i.e. differentially methylated cytosines, DMCs) in fish at thermal limits; however, there was no overlap between DEGs and genes annotated with DMCs in both CTmax and CTmin experiments. The weak interplay between epigenetic modifications and gene expression was further supported by our discoveries of no differentially expressed miRNAs. These findings provide novel insights into the relative role of different molecular mechanisms underlying immediate responses to extreme temperatures and demonstrate non-concordant responses of epigenetic and transcriptional mechanisms to acute temperature stress.

Genomics for monitoring and understanding species responses to global climate change

All life forms across the globe are experiencing drastic changes in environmental conditions as a result of global climate change. These environmental changes are happening rapidly, incur substantial socioeconomic costs, pose threats to biodiversity and diminish a species' potential to adapt to future environments. Understanding and monitoring how organisms respond to human-driven climate change is therefore a major priority for the conservation of biodiversity in a rapidly changing environment. Recent developments in genomic, transcriptomic and epigenomic technologies are enabling unprecedented insights into the evolutionary processes and molecular bases of adaptation. This Review summarizes methods that apply and integrate omics tools to experimentally investigate, monitor and predict how species and communities in the wild cope with global climate change, which is by genetically adapting to new environmental conditions, through range shifts or through phenotypic plasticity. We identify advantages and limitations of each method and discuss future research avenues that would improve our understanding of species' evolutionary responses to global climate change, highlighting the need for holistic, multi-omics approaches to ecosystem monitoring during global climate change.

Smaller herring larval size-at-stage in response to environmental changes is associated with ontogenic processes and stress response

Global change puts coastal systems under pressure, affecting the ecology and physiology of marine organisms. In particular, fish larvae are sensitive to environmental conditions, and their fitness is an important determinant of fish stock recruitment and fluctuations. To assess the combined effects of warming, acidification and change in food quality, herring larvae were reared in a control scenario (11°C*pH 8.0) and a scenario predicted for 2100 (14°C*pH 7.6) crossed with two feeding treatments (enriched in phosphorus and docosahexaenoic acid or not). The experiment lasted from hatching to the beginning of the post-flexion stage (i.e. all fins present) corresponding to 47 days post-hatch (dph) at 14°C and 60 dph at 11°C. Length and stage development were monitored throughout the experiment and the expression of genes involved in growth, metabolic pathways and stress responses were analysed for stage 3 larvae (flexion of the notochord). Although the growth rate was unaffected by acidification and temperature changes, the development was accelerated in the 2100 scenario, where larvae reached the last developmental stage at a smaller size (-8%). We observed no mortality related to treatments and no effect of food quality on the development of herring larvae. However, gene expression analyses revealed that heat shock transcripts expression was higher in the warmer and more acidic treatment. Our findings suggest that the predicted warming and acidification environment are stressful for herring larvae, inducing a decrease in size-at-stage at a precise period of ontogeny. This could either negatively affect survival and recruitment via the extension of the predation window or positively increase the survival by reducing the larval stage duration.

How genomics can help biodiversity conservation

The availability of public genomic resources can greatly assist biodiversity assessment, conservation, and restoration efforts by providing evidence for scientifically informed management decisions. Here we survey the main approaches and applications in biodiversity and conservation genomics, considering practical factors, such as cost, time, prerequisite skills, and current shortcomings of applications. Most approaches perform best in combination with reference genomes from the target species or closely related species. We review case studies to illustrate how reference genomes can facilitate biodiversity research and conservation across the tree of life. We conclude that the time is ripe to view reference genomes as fundamental resources and to integrate their use as a best practice in conservation genomics.

ccImpute: an accurate and scalable consensus clustering based algorithm to impute dropout events in the single-cell RNA-seq data

Background

In recent years, the introduction of single-cell RNA sequencing (scRNA-seq) has enabled the analysis of a cell’s transcriptome at an unprecedented granularity and processing speed. The experimental outcome of applying this technology is a \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$M \times N$$\end{document}M×N matrix containing aggregated mRNA expression counts of M genes and N cell samples. From this matrix, scientists can study how cell protein synthesis changes in response to various factors, for example, disease versus non-disease states in response to a treatment protocol. This technology’s critical challenge is detecting and accurately recording lowly expressed genes. As a result, low expression levels tend to be missed and recorded as zero - an event known as dropout. This makes the lowly expressed genes indistinguishable from true zero expression and different than the low expression present in cells of the same type. This issue makes any subsequent downstream analysis difficult.

Results

To address this problem, we propose an approach to measure cell similarity using consensus clustering and demonstrate an effective and efficient algorithm that takes advantage of this new similarity measure to impute the most probable dropout events in the scRNA-seq datasets. We demonstrate that our approach exceeds the performance of existing imputation approaches while introducing the least amount of new noise as measured by clustering performance characteristics on datasets with known cell identities.

Conclusions

ccImpute is an effective algorithm to correct for dropout events and thus improve downstream analysis of scRNA-seq data. ccImpute is implemented in R and is available at https://github.com/khazum/ccImpute.

Warming Accelerates the Onset of the Molecular Stress Response and Increases Mortality of Larval Atlantic Cod

Temperature profoundly affects ectotherm physiology. Although differential thermal responses influence fitness, thus driving population dynamics and species distributions, our understanding of the molecular architecture underlying these responses is limited, especially during the critical larval stage. Here, using RNA-sequencing of laboratory-reared Atlantic cod (Gadus morhua) larvae of wild origin, we find changes in gene expression in thousands of transcripts consistent with a severe cellular stress response at both ambient and projected (+2°C and +4°C) temperatures. In addition, specific responses to stress, heat, and hypoxia were commonly identified in gene ontology enrichment analyses and 33 of the 44 genes comprising the minimum stress proteome of all organisms were upregulated. Earlier onset of the stress response was evident at higher temperatures; concomitant increased growth and mortality suggests a reduction in fitness. Temporal differences in gene expression levels do not correspond to differences in growing degree days, suggesting negative physiological consequences of warming beyond accelerated development. Because gene expression is costly, we infer that the upregulation of thousands of transcripts in response to warming in larval cod might act as an energetic drain. We hypothesize that the energetically costly stress response, coupled with increased growth rate at warmer temperatures, leads to faster depletion of energy reserves and increased risk of mortality in larval cod. As sea surface temperatures continue to rise over the next century, reduced fitness of Atlantic cod larvae might lead to population declines in this ecologically and socioeconomically important species. Further, our findings expand our understanding of transcriptomic responses to temperature by ectothermic vertebrate larvae beyond the critical first-feeding stage, a time when organisms begin balancing the energetic demands of growth, foraging, development, and maintenance. Linking the molecular basis of a thermal response to key fitness-related traits is fundamentally important to predicting how global warming will affect ectotherms.

Recognizing salinity threats in the climate crisis

Climate change is causing habitat salinity to transform at unprecedented rates across the globe. While much of the research on climate change has focused on rapid shifts in temperature, far less attention has focused on the effects of changes in environmental salinity. Consequently, predictive studies on the physiological, evolutionary, and migratory responses of organisms and populations to the threats of salinity change are relatively lacking. This omission represents a major oversight, given that salinity is among the most important factors that define biogeographic boundaries in aquatic habitats. In this perspective, we briefly touch on responses of organisms and populations to rapid changes in salinity occurring on contemporary time scales. We then discuss factors that might confer resilience to certain taxa, enabling them to survive rapid salinity shifts. Next, we consider approaches for predicting how geographic distributions will shift in response to salinity change. Finally, we identify additional data that are needed to make better predictions in the future. Future studies on climate change should account for the multiple environmental factors that are rapidly changing, especially habitat salinity.

Genomic reaction norms inform predictions of plastic and adaptive responses to climate change

Genomic reaction norms represent the range of gene expression phenotypes (usually mRNA transcript levels) expressed by a genotype along an environmental gradient. Reaction norms derived from common‐garden experiments are powerful approaches for disentangling plastic and adaptive responses to environmental change in natural populations. By treating gene expression as a phenotype in itself, genomic reaction norms represent invaluable tools for exploring causal mechanisms underlying organismal responses to climate change across multiple levels of biodiversity.

Our goal is to provide the context, framework and motivation for applying genomic reaction norms to study the responses of natural populations to climate change.

Here, we describe the utility of integrating genomics with common‐garden‐gradient experiments under a reaction norm analytical framework to answer fundamental questions about phenotypic plasticity, local adaptation, their interaction (i.e. genetic variation in plasticity) and future adaptive potential.

An experimental and analytical framework for constructing and analysing genomic reaction norms is presented within the context of polygenic climate change responses of structured populations with gene flow. Intended for a broad eco‐evo readership, we first briefly review adaptation with gene flow and the importance of understanding the genomic basis and spatial scale of adaptation for conservation and management of structured populations under anthropogenic change. Then, within a high‐dimensional reaction norm framework, we illustrate how to distinguish plastic, differentially expressed (difference in reaction norm intercepts) and differentially plastic (difference in reaction norm slopes) genes, highlighting the areas of opportunity for applying these concepts.

We conclude by discussing how genomic reaction norms can be incorporated into a holistic framework to understand the eco‐evolutionary dynamics of climate change responses from molecules to ecosystems. We aim to inspire researchers to integrate gene expression measurements into common‐garden experimental designs to investigate the genomics of climate change responses as sequencing costs become increasingly accessible.

Physiological Acclimatization in High-Latitude Zooplankton

How individual organisms adapt to non-optimal conditions through physiological acclimatization is central to predicting the consequences of unusual abiotic and biotic conditions such as those produced by marine heat waves. The Northeast Pacific, including the Gulf of Alaska experienced an extreme warming event (2014-2016, "The Blob") that affected all trophic levels leading to large-scale changes in the community. The marine copepod Neocalanus flemingeri is one key member of the subarctic Pacific pelagic ecosystem. During the spring phytoplankton bloom this copepod builds substantial lipid stores as it prepares for its non-feeding adult phase. A three-year comparison of gene expression profiles of copepods collected in Prince William Sound in the Gulf of Alaska between 2015 and 2017 included two high-temperature years (2015 and 2016) and one year with very low phytoplankton abundances (2016). The largest differences in gene expression were between high and low chlorophyll years, and not between warm and cool years. The observed gene expression patterns were indicative of physiological acclimatization. The predominant signal in 2016 was the down-regulation of genes involved in glycolysis and its incoming pathways, consistent with the modulation of metabolic rates in response to prolonged low food conditions. Despite the down-regulation of genes involved in metabolism, there was no evidence of suppression of protein synthesis based on gene expression or behavioral activity. Genes involved in muscle function were up-regulated, and the copepods were actively swimming and responsive to stimuli at collection. However, genes involved in fatty acid metabolism were down-regulated in 2016, suggesting reduced lipid accumulation.

The role of alternative splicing in adaptation and evolution

Regulation of gene expression plays a central role in adaptive divergence and evolution. Although the role of gene regulation in microevolutionary processes is gaining wide acceptance, most studies have only investigated the evolution of transcript levels, ignoring the potentially significant role of transcript structures. We argue that variation in alternative splicing plays an important and widely unexplored role in adaptation (e.g., by increasing transcriptome and/or proteome diversity, or buffering potentially deleterious genetic variation). New studies increasingly highlight the potential for independent evolution in alternative splicing and transcript level, providing alternative paths for selection to act upon. We propose that alternative splicing and transcript levels can provide contrasting, nonredundant mechanisms of equal importance for adaptive diversification of gene function and regulation.

Mediterranean Aquaculture in a Changing Climate: Temperature Effects on Pathogens and Diseases of Three Farmed Fish Species

Climate change is expected to have a drastic effect on aquaculture worldwide. As we move forward with the agenda to increase and diversify aquaculture production, rising temperatures will have a progressively relevant impact on fish farming, linked to a multitude of issues associated with fish welfare. Temperature affects the physiology of both fish and pathogens, and has the potential to lead to significant increases in disease outbreaks within aquaculture systems, resulting in severe financial impacts. Significant shifts in future temperature regimes are projected for the Mediterranean Sea. We therefore aim to review and discuss the existing knowledge relating to disease outbreaks in the context of climate change in Mediterranean finfish aquaculture. The objective is to describe the effects of temperature on the physiology of both fish and pathogens, and moreover to list and discuss the principal diseases of the three main fish species farmed in the Mediterranean, namely gilthead seabream (Sparus aurata), European seabass (Dicentrarchus labrax), and meagre (Argyrosomus regius). We will attempt to link the pathology of each disease to a specific temperature range, while discussing potential future disease threats associated with the available climate change trends for the Mediterranean Sea.

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.

Environmental RNA: A Revolution in Ecological Resolution?

Current advancements in environmental RNA (eRNA) exploit its relatively fast turnover rate relative to environmental DNA (eDNA) to assess 'metabolically active' or temporally/spatially recent community diversity. However, this focus significantly underutilizes the trove of potential ecological information encrypted in eRNA. Here, we argue for pushing beyond current species-level eDNA detection capabilities by using eRNA to detect any organisms with unique eRNA profiles, potentially including different life-history stages, sexes, or even specific phenotypes within a species. We also discuss the future of eRNA as a means of assessing the physiological status of organisms and the ecological health of populations and communities, reflecting ecosystem-level conditions. We posit that eRNA has the potential to significantly improve the resolution of organism detection, biological monitoring, and biomonitoring applications in ecology.

Comparative transcriptomics of spotted seatrout (Cynoscion nebulosus) populations to cold and heat stress

Resilience to climate change depends on a species' adaptive potential and phenotypic plasticity. The latter can enhance survival of individual organisms during short periods of extreme environmental perturbations, allowing genetic adaptation to take place over generations. Along the U.S. East Coast, estuarine-dependent spotted seatrout (Cynoscion nebulosus) populations span a steep temperature gradient that provides an ideal opportunity to explore the molecular basis of phenotypic plasticity. Genetically distinct spotted seatrout sampled from a northern and a southern population were exposed to acute cold and heat stress (5 biological replicates in each treatment and control group), and their transcriptomic responses were compared using RNA-sequencing (RNA-seq). The southern population showed a larger transcriptomic response to acute cold stress, whereas the northern population showed a larger transcriptomic response to acute heat stress compared with their respective population controls. Shared transcripts showing significant differences in expression levels were predominantly enriched in pathways that included metabolism, transcriptional regulation, and immune response. In response to heat stress, only the northern population significantly upregulated genes in the apoptosis pathway, which could suggest greater vulnerability to future heat waves in this population as compared to the southern population. Genes showing population-specific patterns of expression, including hpt, acot, hspa5, and hsc71, are candidates for future studies aiming to monitor intraspecific differences in temperature stress responses in spotted seatrout. Our findings contribute to the current understanding of phenotypic plasticity and provide a basis for predicting the response of a eurythermal fish species to future extreme temperatures.

Transcriptomic responses to heat stress in gill and liver of endangered Brachymystax lenok tsinlingensis

Global warming significantly affects fish, particularly cold-water fish, because increased temperature adversely impacts their abilities to grow or reproduce, and eventually influences their fitness or even causes death. To survive, fish may alter their distribution or behavior to avoid the stress, and perhaps acclimate or evolve resistance to the elevated temperature. Brachymystax lenok tsinlingensis is an endangered cold-water species in China, and it has been found to alter the altitudinal distribution, decrease swimming efficiency and develop resistance under heat exposure, which badly impact the continuing conservation work. To better protect them, it is essential to understand how they respond to thermal stress behaviorally and physiologically. Therefore, the fish were exposed to 24.5 °C and based on the time taken for them to lose equilibrium, they were separately sampled as sensitive and tolerant groups. Both gill and liver tissues were collected from both groups for transcriptome sequencing. Sequencing results demonstrated that control and tolerant groups were similar in transcriptomic patterns and sensitive groups differentially expressed more genes than tolerant ones, suggesting the gene expression of tolerant groups may return to base levels as exposure time increased. Tissue differences were the major factor affecting gene expression, and they also displayed different physiological responses to heat stress. Consistent with other studies, heat shock response, immune response, metabolic adjustment and ion transport were found to be triggered after exposed to elevated temperature. The findings would contribute to a better understanding of responding mechanisms of fish to thermal stress and provide guidance for future conservation programs.

Footprints of global change in marine life: Inferring past environment based on DNA methylation and gene expression marks

Ocean global warming affects the distribution, life history and physiology of marine life. Extreme events, like marine heatwaves, are increasing in frequency and intensity. During sensitive stages of early fish development, the consequences may be long-lasting and mediated by epigenetic mechanisms. Here, we used European sea bass as a model to study the possible long-lasting effects of a marine heatwave during early development. We measured DNA methylation and gene expression in four tissues (brain, muscle, liver and testis) and detected differentially methylated regions (DMRs). Six genes were differentially expressed and contained DMRs three years after exposure to increased temperature, indicating direct phenotypic consequences and representing persistent changes. Interestingly, nine genes contained DMRs around the same genomic regions across tissues, therefore consisting of common footprints of developmental temperature in environmentally responsive loci. These loci are, to our knowledge, the first metastable epialleles (MEs) described in fish. MEs may serve as biomarkers to infer past life history events linked with persistent consequences. These results highlight the importance of subtle phenotypic changes mediated by epigenetics to extreme weather events during sensitive life stages. Also, to our knowledge, it is the first time the molecular effects of a marine heatwave during the lifetime of individuals are assessed. MEs could be used in surveillance programs aimed at determining the footprints of climate change on marine life. Our study paves the way for the identification of conserved MEs that respond equally to environmental perturbations across species. Conserved MEs would constitute a tool of assessment of global change effects in marine life at a large scale.

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.

Comprehensive phenotyping and multi-omic profiling in the toxicity assessment of nanopolystyrene with different surface properties

There is a growing concern regarding the toxic effects of terrestrial nanoplastic contaminants. However, an all-encompassing phenotyping- and omics-based strategy for the toxicity assessment of nanoplastics with different surface properties on soil living organisms remains to be established. Herein, we devised a comprehensive phenotyping and multi-omic profiling method to examine the molecular disturbance of nanopolystyrene (PS)-exposed Caenorhabditis elegans. The exposure time was 24 h with either 1 μg/mL or 10 μg/mL of PS. We found that PS considerably affected the reproduction and locomotion, as well as increased the oxidative stress of worms regardless of their surface properties. Nevertheless, each type of PS affected the metabolome and lipidome of the nematodes differently. Uncharged PS (PS-N) triggered significant metabolic disturbances, whereas the metabolic influences from PS-NH₂ and PS-COOH were subtle. The dysregulated transcriptome profiles of PS-N were strongly associated with the metabolic pathways. Besides, the altered expression of several genes associated with autophagy and longevity was observed. Collectively, we demonstrated that comprehensive phenotyping and omics-based profiling establish a practical framework that allows us to gain deeper insights into the maladaptive consequences of PS in nematodes. It can be utilized for the evaluation of other environmental contaminants in the terrestrial ecosystem.

Transcriptomic evidences of local thermal adaptation for the native fish Colossoma macropomum (Cuvier, 1818)

Brazil has five climatically distinct regions, with an annual average temperature difference up to 14 ºC between the northern and southern extremes. Environmental variation of this magnitude can lead to new genetic patterns among farmed fish populations. Genetically differentiated populations of tambaqui (Colossoma macropomum Cuvier, 1818), an important freshwater fish for Brazilian continental aquaculture, may be associated with regional adaptation. In this study, we selected tambaquis raised in two thermally distinct regions, belonging to different latitudes, to test this hypothesis. De novo transcriptome analysis was performed to compare the significant differences of genes expressed in the liver of juvenile tambaqui from a northern population (Balbina) and a southeastern population (Brumado). In total, 2,410 genes were differentially expressed (1,196 in Balbina and 1,214 in Brumado). Many of the genes are involved in a multitude of biological functions such as biosynthetic processes, homeostasis, biorhythm, immunity, cell signaling, ribosome biogenesis, modification of proteins, intracellular transport, structure/cytoskeleton, and catalytic activity. Enrichment analysis based on biological networks showed a different protein interaction profile for each population, whose encoding genes may play potential functions in local thermal adaptation of fish to their respective farming environments.

Cytosine methylation patterns suggest a role of methylation in plastic and adaptive responses to temperature in European grayling (Thymallus thymallus) populations

Temperature is a key environmental parameter affecting both the phenotypes and distributions of organisms, particularly ectotherms. Rapid organismal responses to thermal environmental changes have been described for several ectotherms; however, the underlying molecular mechanisms often remain unclear. Here, we studied whole genome cytosine methylation patterns of European grayling (Thymallus thymallus) embryos from five populations with contemporary adaptations of early life history traits at either 'colder' or 'warmer' spawning grounds. We reared fish embryos in a common garden experiment using two temperatures that resembled the 'colder' and 'warmer' conditions of the natal natural environments. Genome-wide methylation patterns were similar in populations originating from colder thermal origin subpopulations, whereas single nucleotide polymorphisms uncovered from the same data identified strong population structure among isolated populations, but limited structure among interconnected populations. This was surprising because the previously studied gene expression response among populations was mostly plastic, and mainly influenced by the developmental temperature. These findings support the hypothesis of the magnified role of epigenetic mechanisms in modulating plasticity. The abundance of consistently changing methylation loci between two warmer-to-colder thermal origin population pairs suggests that local adaptation has shaped the observed methylation patterns. The dynamic nature of the methylomes was further highlighted by genome-wide and site-specific plastic responses. Our findings support both the presence of a plastic response in a subset of CpG loci, and the evolutionary role of methylation divergence between populations adapting to contrasting thermal environments.

How will farmed populations of freshwater fish deal with the extreme climate scenario in 2100? Transcriptional responses of Colossoma macropomum from two Brazilian climate regions

Tambaqui (Colossoma macropomum Cuvier, 1818) is an endemic fish of the Amazon and Orinoco basins, and it is the most economically important native species in Brazil being raised in five climatically distinct regions. In the face of current global warming, environmental variations in farm ponds represent additional challenges that may drive new adaptive regional genetic variations among broodstocks of tambaqui. In an experimental context based on the high-emission scenario of the 5th Intergovernmental Panel on Climate Change (IPCC) report, we used two farmed tambaqui populations to test this hypothesis. RNA-seq transcriptome analysis was performed in the liver of juvenile tambaqui from northern (Balbina Experimental Station, Balbina, AM) and southeastern (Brumado Fish Farming, Mogi Mirim, SP) Brazilian regions kept for 30 days in artificial environmental rooms mimicking the current and extreme climate scenarios. Three Illumina MiSeq runs produced close to 120 million 500 bp paired-end reads; 191,139 contigs were assembled with N50 = 1595. 355 genes were differentially expressed for both populations in response to the extreme scenario. After enrichment analysis, each population presented a core set of genes to cope with climate change. Northern fish induced genes related to the cellular response to stress, activation of MAPK activity, response to unfolded protein, protein metabolism and cellular response to DNA damage stimuli. Genes biologically involved in regulating cell proliferation, protein stabilisation and protein ubiquitination for degradation through the ubiquitin-proteasome system were downregulated. Genes associated with biological processes, including the cellular response to stress, MAPK cascade activation, homeostatic processes and positive regulation of immune responses were upregulated in southeastern fish. The downregulated genes were related to cytoskeleton organisation, energy metabolism, and the regulation of transcription and biological rhythms. Our findings reveal the signatures of promising candidate genes involved in the regional plasticity of each population of tambaqui in dealing with upcoming climate changes.

Transcriptome profiling reveals exposure to predicted end-of-century ocean acidification as a stealth stressor for Atlantic cod larvae

Ocean acidification (OA), a direct consequence of increasing atmospheric CO₂ concentration dissolving in ocean waters, is impacting many fish species. Little is known about the molecular mechanisms underlying the observed physiological impacts in fish. We used RNAseq to characterize the transcriptome of 3 different larval stages of Atlantic cod (Gadus morhua) exposed to simulated OA at levels (1179 µatm CO₂) representing end-of-century predictions compared to controls (503 µatm CO₂), which were shown to induce tissue damage and elevated mortality in G. morhua. Only few genes were differentially expressed in 6 and 13 days-post-hatching (dph) (3 and 16 genes, respectively), during a period when maximal mortality as a response to elevated pCO₂ occurred. At 36 dph, 1413 genes were differentially expressed, most likely caused by developmental asynchrony between the treatment groups, with individuals under OA growing faster. A target gene analysis revealed only few genes of the universal and well-defined cellular stress response to be differentially expressed. We thus suggest that predicted ocean acidification levels constitute a "stealth stress" for early Atlantic cod larvae, with a rapid breakdown of cellular homeostasis leading to organismal death that was missed even with an 8-fold replication implemented in this study.

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.

Genomics of Developmental Plasticity in Animals

Developmental plasticity refers to the property by which the same genotype produces distinct phenotypes depending on the environmental conditions under which development takes place. By allowing organisms to produce phenotypes adjusted to the conditions that adults will experience, developmental plasticity can provide the means to cope with environmental heterogeneity. Developmental plasticity can be adaptive and its evolution can be shaped by natural selection. It has also been suggested that developmental plasticity can facilitate adaptation and promote diversification. Here, we summarize current knowledge on the evolution of plasticity and on the impact of plasticity on adaptive evolution, and we identify recent advances and important open questions about the genomics of developmental plasticity in animals. We give special attention to studies using transcriptomics to identify genes whose expression changes across developmental environments and studies using genetic mapping to identify loci that contribute to variation in plasticity and can fuel its evolution.

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.

Transcriptomic responses to heat stress in rainbow trout Oncorhynchus mykiss head kidney

Rainbow trout (Oncorhynchus mykiss) are widely cultured throughout the word for commercial aquaculture. However, as a cold-water species, rainbow trout are highly susceptible to heat stress, which may cause pathological signs or diseases by alleviating the immune roles and then lead to mass mortality. Understanding the molecular mechanisms that occur in the rainbow trout in response to heat stress will be useful to decrease heat stress-related morbidity and mortality in trout aquaculture. In the present study, we conducted transcriptome analysis of head kidney tissue in rainbow trout under heat-stress (24 °C) and control (18 °C) conditions, to identify heat stress-induced genes and pathways. More than 281 million clean reads were generated from six head kidney libraries. Using an adjusted P-value of P < 0.05 as the threshold, a total of 443 differentially expressed genes (DEGs) were identified, including members of the HSP90, HSP70, HSP60, and HSP40 family and several cofactors or cochaperones. The RNA-seq results were confirmed by RT-qPCR. Gene ontology and Kyoto Encyclopedia of Genes and Genomes enrichment analysis of DEGs were performed. Many genes involved in maintaining homeostasis or adapting to stress and stimuli were highly induced in response to high temperature. The most significantly enriched pathway was "Protein processing in endoplasmic reticulum (ER)", a quality control system that ensures correct protein folding or degradation of misfolded polypeptides by ER-associated degradation. Other signaling pathways involved in regulation of immune system and post-transcriptional regulation of spliceosome were also critical for thermal adaptation. These findings improve our understanding of the molecular mechanisms of heat stress responses and are useful to develop strategies for the improvement of rainbow trout survival rate during summer high-temperature period.

RNA sequencing, de novo assembly and differential analysis of the gill transcriptome of freshwater climbing perch Anabas testudineus after 6 days of seawater exposure

To obtain transcriptomic insights into branchial responses to salinity challenge in Anabas testudineus, this study employed RNA sequencing (RNA-Seq) to analyse the gill transcriptome of A. testudineus exposed to seawater (SW) for 6 days compared with the freshwater (FW) control group. A combined FW and SW gill transcriptome was de novo assembled from 169.9 million 101 bp paired-end reads. In silico validation employing 17 A. testudineus Sanger full-length coding sequences showed that 15/17 of them had greater than 80% of their sequences aligned to the de novo assembled contigs where 5/17 had their full-length (100%) aligned and 9/17 had greater than 90% of their sequences aligned. The combined FW and SW gill transcriptome was mapped to 13,780 unique human identifiers at E-value ≤1.0E-20 while 952 and 886 identifiers were determined as up and down-regulated by 1.5 fold, respectively, in the gills of A. testudineus in SW when compared with FW. These genes were found to be associated with at least 23 biological processes. A larger proportion of genes encoding enzymes and transporters associated with molecular transport, energy production, metabolisms were up-regulated, while a larger proportion of genes encoding transmembrane receptors, G-protein coupled receptors, kinases and transcription regulators associated with cell cycle, growth, development, signalling, morphology and gene expression were relatively lower in the gills of A. testudineus in SW when compared with FW. High correlation (R = 0.99) was observed between RNA-Seq data and real-time quantitative PCR validation for 13 selected genes. The transcriptomic sequence information will facilitate development of molecular resources and tools while the findings will provide insights for future studies into branchial iono-osmoregulation and related cellular processes in A. testudineus.

Liver transcriptome analysis reveals extensive transcriptional plasticity during acclimation to low salinity in Cynoglossus semilaevis

BACKGROUND

Salinity is an important abiotic stress that influences the physiological and metabolic activity, reproduction, growth and development of marine fish. It has been suggested that half-smooth tongue sole (Cynoglossus semilaevis), a euryhaline fish species, uses a large amount of energy to maintain osmotic pressure balance when exposed to fluctuations in salinity. To delineate the molecular response of C. semilaevis to different levels of salinity, we performed RNA-seq analysis of the liver to identify the genes and molecular and biological processes involved in responding to salinity changes.

RESULTS

The present study yielded 330.4 million clean reads, of which 83.9% were successfully mapped to the reference genome of C. semilaevis. One hundred twenty-eight differentially expressed genes (DEGs), including 43 up-regulated genes and 85 down-regulated genes, were identified. These DEGs were highly represented in metabolic pathways, steroid biosynthesis, terpenoid backbone biosynthesis, butanoate metabolism, glycerolipid metabolism and the 2-oxocarboxylic acid metabolism pathway. In addition, genes involved in metabolism, osmoregulation and ion transport, signal transduction, immune response and stress response, and cytoskeleton remodeling were affected during acclimation to low salinity. Genes acat2, fdps, hmgcr, hmgcs1, mvk, pmvk, ebp, lss, dhcr7, and dhcr24 were up-regulated and abat, ddc, acy1 were down-regulated in metabolic pathways. Genes aqp10 and slc6a6 were down-regulated in osmoregulation and ion transport. Genes abat, fdps, hmgcs1, mvk, pmvk and dhcr7 were first reported to be associated with salinity adaptation in teleosts.

CONCLUSIONS

Our results revealed that metabolic pathways, especially lipid metabolism were important for salinity adaptation. The candidate genes identified from this study provide a basis for further studies to investigate the molecular mechanism of salinity adaptation and transcriptional plasticity in marine fish.

Similarities in temperature-dependent gene expression plasticity across timescales in threespine stickleback (Gasterosteus aculeatus)

Phenotypic plasticity occurs at a variety of timescales, but little is known about the degree to which plastic responses at different timescales are associated with similar underlying molecular processes, which is critical for assessing the effects of plasticity on evolutionary trajectories. To address this issue, we identified differential gene expression in response to developmental temperature in the muscle transcriptome of adult threespine stickleback (Gasterosteus aculeatus) exposed to 12, 18 and 24°C until hatch and then held at 18°C for 9 months and compared these results to differential gene expression in response to adult thermal acclimation in stickleback developed at 18°C and then acclimated to 5 and 25°C as adults. Adult thermal acclimation affected the expression of 7,940 and 7,015 genes in response to cold and warm acclimation, respectively, and 4,851 of these genes responded in both treatments. In contrast, the expression of only 33 and 29 genes was affected by cold and warm development, respectively. The majority of the genes affected by developmental temperature were also affected by adult acclimation temperature. Many genes that were differentially expressed as a result of adult acclimation were associated with previously identified temperature-dependent effects on DNA methylation patterns, suggesting a role of epigenetic mechanisms in regulating gene expression plasticity during acclimation. Taken together, these results demonstrate similarities between the persistent effects of developmental plasticity on gene expression and the effects of adult thermal acclimation, emphasizing the potential for mechanistic links between plasticity acting at these different life stages.