Functional genomics of physiological plasticity and local adaptation in killifish

Effect of climate change on hematotoxicity/hepatoxicity oxidative stress, Oncorhynchus mykiss, under controlled conditions

Described as the 'main ecological factor', temperature, strongly affects the physiological stress responses of fish. In order to evaluate the effects of temperature variations on fish culture and food value chain, the present study was designed as a climate change model. Furthermore, the present study provides a theoretical basis for a better understanding of the mechanisms of the environmentally induced changes. In this direction, we examined the blood physiology and oxidative stress responses induced by temperature variation in the rainbow trout, a temperature-sensitive cold-water fish. The obtained results showed that climate changes promoted the inhibited activities' expressions and the development of potential tissue and hematological defense mechanisms against temperature-induced toxic damage. This study showed that climate change could be a subset of the studies on the stress physiology in aquaculture, which can be developed for new experimental designs and research collaborations. Furthermore, it highlights knowledge gaps to guide future research in this emerging field.

Disturbance of oxidation/antioxidant status and histopathological damage in tsinling lenok trout under acute thermal stress

The tsinling lenok trout (Brachymystax Lenok tsinlingensis) is a typical land-locked cold salmon. In this study, through the acute high temperature stress (16, 24, and 26°C), samples were taken at multiple temperature points to analyze the dynamic changes of serum non-specific immune indexes and histopathological changes of tissues of tsinling lenok trout. The histopathological investigation of different studied tissues revealed an increase of histological lesions' frequency and severity with increasing water temperature. The activity of T-SOD, GSH-Px, CAT, ACP, and LZM and MDA concentration are all impacted by high temperature stress. The activities of T-SOD, GSH-Px, and ACP are significantly lower in temperatures higher than 16°C (P<0.05). However, with the increase of water temperature, MDA content increased significantly. The activities of CAT and LZM were the highest at 24°C, which were significantly higher than those at 26°C (P<0.05). The above results indicate that 24°C is a "critical high temperature point" for tsinling lenok trout under heat stress, and this temperature point may be the critical point for tsinling lenok trout to enter "damaged" from adaptive adjustment. Our results can provide a theoretical basis for the development of genetic breeding, improvement, and control measures of heat stress in tsinling lenok trout in the future.

Transcriptomic Response of the Liver Tissue in Trachinotus ovatus to Acute Heat Stress

Trachinotus ovatus is a major economically important cultured marine fish in the South China Sea. However, extreme weather and increased culture density result in uncontrollable problems, such as increases in water temperature and a decline in dissolved oxygen (DO), hindering the high-quality development of aquaculture. In this study, liver transcriptional profiles of T. ovatus were investigated under acute high-temperature stress (31 °C and 34 °C) and normal water temperature (27 °C) using RNA sequencing (RNA-Seq) technology. Differential expression analysis and STEM analysis showed that 1347 differentially expressed genes (DEGs) and four significant profiles (profiles 0, 3, 4, and 7) were screened, respectively. Of these DEGs, some genes involved in heat shock protein (HSPs), hypoxic adaptation, and glycolysis were up-regulated, while some genes involved in the ubiquitin-proteasome system (UPS) and fatty acid metabolism were down-regulated. Our results suggest that protein dynamic balance and function, hypoxia adaptation, and energy metabolism transformation are crucial in response to acute high-temperature stress. Our findings contribute to understanding the molecular response mechanism of T. ovatus under acute heat stress, which may provide some reference for studying the molecular mechanisms of other fish in response to heat stress.

Abiotic factors and heavy metals defining eco-physiological niche in fish

Ecosystem dynamics undergoing alterations in structure and function highlights the need to look into the relations between ecological parameters and organismal fitness and tolerance. Ecophysiological studies are used to understand how organisms adapt to and cope up with environmental stress. Current study uses a process-based approach to model physiochemical parameters regarding seven different fish species. Species respond to climatic variations via acclimation or adaptation through physiological plasticity. Four sites are differentiated into two types based on the water quality parameters and metal contamination. Seven fish species are clustered into two groups, each group depicting separate pattern of response in similar habitat. In this manner, biomarkers from three different physiological axes- stress, reproduction, and neurology were taken to determine the organism's ecological niche. Cortisol, Testosterone, Estradiol, and AChE are the signature molecules estimated for the said physiological axes. The ordination technique, nonmetric multidimensional scaling, has been utilized to visualize the differentiated physiological response to changing environmental conditions. Then, Bayesian Model Averaging (BMA) was used to identify the factors that play a key role in refining the stress physiology and determining the niche. Current study confirms different species belonging to similar habitats respond to various environmental and physiological factors in a different manner as various biomarkers respond in a species-specific pattern that induces the choice of habitat preference controlling its ecophysiological niche. In the present study, it is quite apparent that adaptive mechanism of fish to environmental stress is achieved through modification of physiological mechanisms through a panel of biochemical markers. These markers organize a cascade of physiological event at various levels including reproduction.

Comparative transcriptome reveals the thermal stress response differences between Heilongjiang population and Xinjiang population of Lota lota

Some cold-water fishes are particularly sensitive to the water temperature increasing caused by current global warming. However, the alterations in the physiology and behavior of infraspecific populations living in heterogeneous landscapes in response to water temperature increasing were significantly different. Consequently, understanding the impact of temperature increasing on different populations may be crucial for the conservation of cold-water fishes in the context of global warming. The burbot is the only freshwater specie in Gadiformes. To better understand the differences of different populations of burbot under similar thermal stress, Lota lota was selected as the research objects. Firstly, RNA-seq was applied to identify the transcriptomic responses of Heilongjiang population exposed to three temperature gradients (0 °C, 18 °C and 28 °C). Compared with 0 °C, 4216 and 12,657 genes were significantly differentially expressed at 18 °C and 28 °C, respectively. Meanwhile, 49 genes were significantly differentially expressed in three temperature pairs and these genes were presumed to involve in stress response process, immunologic process, reproductive process, development process, material metabolism process, signal transduction process, spermatogenesis process and cell apoptosis process. The response differences of two L. lota populations to similar thermal stress were compared and the results showed that they have different gene expression responses (the number of differentially expression genes and biological processes). The lower annual temperature of the Heilongjiang River might make it more sensitive to temperature increasing. Based on the comparative transcriptome analyses, 12 orthologous genes were considered as the potential regulators of L. lota preference for cold-water environment and these genes are potentially related to the immunologic process, reproductive process, development process, signal transduction process, and cell apoptosis process. Those results can provide basic information for the rational development of conservation strategies of different L. lota populations under the background of global warming.

Transcriptome analysis reveals fluid shear stress (FSS) and atherosclerosis pathway as a candidate molecular mechanism of short-term low salinity stress tolerance in abalone

Background

Transcriptome sequencing is an effective tool to reveal the essential genes and pathways underlying countless biotic and abiotic stress adaptation mechanisms. Although severely challenged by diverse environmental conditions, the Pacific abalone Haliotis discus hannai remains a high-value aquaculture mollusk and a Chinese predominantly cultured abalone species. Salinity is one of such environmental factors whose fluctuation could significantly affect the abalone’s cellular and molecular immune responses and result in high mortality and reduced growth rate during prolonged exposure. Meanwhile, hybrids have shown superiority in tolerating diverse environmental stresses over their purebred counterparts and have gained admiration in the Chinese abalone aquaculture industry.

The objective of this study was to investigate the molecular and cellular mechanisms of low salinity adaptation in abalone. Therefore, this study used transcriptome analysis of the gill tissues and flow cytometric analysis of hemolymph of H. discus hannai (DD) and interspecific hybrid H. discus hannai ♀ x H. fulgens ♂ (DF) during low salinity exposure. Also, the survival and growth rate of the species under various salinities were assessed.

Results

The transcriptome data revealed that the differentially expressed genes (DEGs) were significantly enriched on the fluid shear stress and atherosclerosis (FSS) pathway. Meanwhile, the expression profiles of some essential genes involved in this pathway suggest that abalone significantly up-regulated calmodulin-4 (CaM-4) and heat-shock protein90 (HSP90), and significantly down-regulated tumor necrosis factor (TNF), bone morphogenetic protein-4 (BMP-4), and nuclear factor kappa B (NF-kB).

Also, the hybrid DF showed significantly higher and sustained expression of CaM and HSP90, significantly higher phagocytosis, significantly lower hemocyte mortality, and significantly higher survival at low salinity, suggesting a more active molecular and hemocyte-mediated immune response and a more efficient capacity to tolerate low salinity than DD.

Conclusions

Our study argues that the abalone CaM gene might be necessary to maintain ion equilibrium while HSP90 can offset the adverse changes caused by low salinity, thereby preventing damage to gill epithelial cells (ECs).

The data reveal a potential molecular mechanism by which abalone responds to low salinity and confirms that hybridization could be a method for breeding more stress-resilient aquatic species.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08611-8.

Nature versus nurture in heat stress induced learning between inbred and outbred populations of Lymnaea stagnalis

Changing environmental conditions often lead to microevolution of traits that are adaptive under the current selection pressure. Currently, one of the major selection pressures is the rise in temperatures globally that has a severe impact on the behavioral ecology of animals. However, the role of thermal stress on neuronal plasticity and memory formation is not well understood. Thermal tolerance and sensitivity to heat stress show variation across populations of the same species experiencing different thermal regimes. We used two populations of the pond snail Lymnaea stagnalis: one lab-bred W-snails and the other wild Delta snails to test heat shock induced learning and memory formation for the Garcia effect learning paradigm. In Garcia effect, a single pairing of a heat stressor (30 °C for 1h) with a novel taste results in a taste-specific negative hedonic shift lasting 24h as long-term memory (LTM) in lab bred W-snails. In this study we used a repeated heat stress procedure to test for increased or decreased sensitivity to the heat before testing for the Garcia effect. We found that lab-bred W-snails show increased sensitivity to heat stress after repeated heat exposure for 7days, leading to enhanced LTM for Garcia effect with only 15min of heat exposure instead of standard 1h. Surprisingly, the freshly collected wild snails do not show Garcia effect. Additionally, F1 generation of wild snails raised and maintained under laboratory conditions still retain their heat stress tolerance similar to their parents and do not show a Garcia effect under standard learning paradigm or even after repeated heat stressor. Thus, we found a differential effect of heat stress on memory formation in wild and lab bred snails. Most interestingly we also show that local environmental (temperature) conditions for one generation is not enough to alter thermal sensitivity in a wild population of L. stagnalis.

Metal Accumulation and Ion Regulation in the Fish Hyphessobrycon luetkenii Living in a Site Chronically Contaminated by Copper: Insights from Translocation Experiments

Fish living in the João Dias creek (southern Brazil) have to deal with trace-metal contamination in the long-term basis, as this aquatic environment has been historically impacted by copper mining activities. In order to survive in this harsh environment, the local biota had to develop adaptations related to pollution tolerance. The aim of this study was to test if biochemical mechanisms related to osmoregulation were among these adaptations, using translocation experiments. Water ionic and trace-metal compositions were measured in a nonmetal impacted site (NMIS) and in a metal impacted site (MIS) of this creek. Also, whole-body metal accumulation, ion concentration and branchial enzyme activity (Na,K-ATPase and carbonic anhydrase) were evaluated in Hyphessobrycon luetkenii. In both NMIS and MIS, fish were collected and immediately stored, kept caged or translocated from sites. The result shows that waterborne Cu was 3.4-fold higher at the MIS. Accordingly, animals that had contact with this site showed elevated whole-body Cu levels. Moreover, both translocated groups showed elevated Na,K-ATPase activity. Additionally, fish translocated from the NMIS to the MIS showed lower carbonic anhydrase activity. These findings indicate that H. luetkenii chronically or acutely exposed to naturally elevated waterborne Cu showed a rapid Cu bioaccumulation but was unable to readily excrete it. Moreover, classic Cu osmoregulatory toxicity related to Na,K-ATPase inhibition was not observed. Conversely, impacts in ammonia excretion related to carbonic anhydrase inhibition may have occurred.

Climate differently influences the genomic patterns of two sympatric marine fish species

Climate influences population genetic variation in marine species. Capturing these impacts remains challenging for marine fishes which disperse over large geographical scales spanning steep environmental gradients. It requires the extensive spatial sampling of individuals or populations, representative of seascape heterogeneity, combined with a set of highly informative molecular markers capable of revealing climatic-associated genetic variations. We explored how space, dispersal and environment shape the genomic patterns of two sympatric fish species in the Mediterranean Sea, which ranks among the oceanic basins most affected by climate change and human pressure. We hypothesized that the population structure and climate-associated genomic signatures of selection would be stronger in the less mobile species, as restricted gene flow tends to facilitate the fixation of locally adapted alleles. To test our hypothesis, we genotyped two species with contrasting dispersal abilities: the white seabream Diplodus sargus and the striped red mullet Mullus surmuletus. We collected 823 individuals and used genotyping by sequencing (GBS) to detect 8,206 single nucleotide polymorphisms (SNPs) for the seabream and 2,794 for the mullet. For each species, we identified highly differentiated genomic regions (i.e. outliers) and disentangled the relative contribution of space, dispersal and environmental variables (climate, marine primary productivity) on the outliers' genetic structure to test the prevalence of gene flow and local adaptation. We observed contrasting patterns of gene flow and adaptive genetic variation between the two species. The seabream showed a distinct Alboran sea population and panmixia across the Mediterranean Sea. The mullet revealed additional differentiation within the Mediterranean Sea that was significantly correlated to summer and winter temperatures, as well as marine primary productivity. Functional annotation of the climate-associated outlier SNPs then identified candidate genes involved in heat tolerance that could be examined to further predict species' responses to climate change. Our results illustrate the key steps of a comparative seascape genomics study aiming to unravel the evolutionary processes at play in marine species, to better anticipate their response to climate change. Defining population adaptation capacities and environmental niches can then serve to incorporate evolutionary processes into species conservation planning.

Reductionist approaches to the study of ionoregulation in fishes

The mechanisms underlying ionoregulation in fishes have been studied for nearly a century, and reductionist methods have been applied at all levels of biological organization in this field of research. The complex nature of ionoregulatory systems in fishes makes them ideally suited to reductionist methods and our collective understanding has been dramatically shaped by their use. This review provides an overview of the broad suite of techniques used to elucidate ionoregulatory mechanisms in fishes, from the whole-animal level down to the gene, discussing some of the advantages and disadvantages of these methods. We provide a roadmap for understanding and appreciating the work that has formed the current models of organismal, endocrine, cellular, molecular, and genetic regulation of ion balance in fishes and highlight the contribution that reductionist techniques have made to some of the fundamental leaps forward in the field throughout its history.

Molecular mechanisms of local adaptation for salt-tolerance in a treefrog

Salinization is a global phenomenon affecting ecosystems and forcing freshwater organisms to deal with increasing levels of ionic stress. However, our understanding of mechanisms that permit salt tolerance in amphibians is limited. This study investigates mechanisms of salt tolerance in locally adapted, coastal populations of a treefrog, Hyla cinerea. Using a common garden experiment, we (i) determine the extent that environment (i.e., embryonic and larval saltwater exposure) or genotype (i.e., coastal vs. inland) affects developmental benchmarks and transcriptome expression, and (ii) identify genes that may underpin differences in saltwater tolerance. Differences in gene expression, survival, and plasma osmolality were most strongly associated with genotype. Population genetic analyses on expressed genes also delineated coastal and inland groups based on genetic similarity. Coastal populations differentially expressed osmoregulatory genes including ion transporters (atp1b1, atp6V1g2, slc26a), cellular adhesion components (cdh26, cldn1, gjb3, ocln), and cytoskeletal components (odc1-a, tgm3). Several of these genes are the same genes expressed by euryhaline fish after exposure to freshwater, which is a novel finding for North American amphibians and suggests that these genes may be associated with local salinity adaptation. Coastal populations also highly expressed glycerol-3-phosphate dehydrogenase 1 (gpd1), which indicates they use glycerol as a compatible osmolyte to reduce water loss - another mechanism of saltwater tolerance previously unknown in frogs. These data signify that Hyla cinerea inhabiting coastal, brackish wetlands have evolved a salt-tolerant ecotype, and highlights novel candidate pathways that can lead to salt tolerance in freshwater organisms facing habitat salinization.

Transcriptome analysis of the Larimichthys polyactis under heat and cold stress

The small yellow croaker (Larimichthys polyactis) is an economically important marine fish that is widely distributed in the East Sea, Yellow Sea and Bohai of China. However, the wild populations of L. polyactis are severely depleted, and there is currently a developing large-scale artificial propagation of this fish for aquaculture. However, the current variety of L. polyactis that is cultivated is not capable to coping with large fluctuations in temperature. Therefore, it is important to understand the molecular mechanisms that are activated in response to temperature stress in the L. polyactis. Here, we conducted transcriptomic analysis of the liver of L. polyactis under heat and cold stress. A total of 270,844,888, 265,727,006 and 259,666,218 clean reads were generated from high temperature group, low temperature group and control group, respectively, and comparing expression of genes in these transcriptomes, 10,878 unigenes that were differential expressed were identified. Seventeen of the differentially expressed unigenes were validated by qRT-PCR. Pathway enrichment analysis identified that the ER pathway, immune signaling pathway and metabolic response pathway were affected by temperature stress. The results of this study provide a comprehensive overview of temperature stress-induced transcriptional patterns in liver tissues of the L. polyactis. In addition, these results can guide future molecular studies of heat and cold stress response in this species for improving the stock used for aquaculture.

Is aquaporin-3 involved in water-permeability changes in the killifish during hypoxia and normoxic recovery, in freshwater or seawater?

Aquaporins are the predominant water-transporting proteins in vertebrates, but only a handful of studies have investigated aquaporin function in fish, particularly in mediating water permeability during salinity challenges. Even less is known about aquaporin function in hypoxia (low oxygen), which can profoundly affect gill function. Fish deprived of oxygen typically enlarge gill surface area and shrink the water-to-blood diffusion distance, to facilitate oxygen uptake into the bloodstream. However, these alterations to gill morphology can result in unfavorable water and ion fluxes. Thus, there exists an osmorespiratory compromise, whereby fish must try to balance high branchial gas exchange with low ion and water permeability. Furthermore, the gills of seawater and freshwater teleosts have substantially different functions with respect to osmotic and ion fluxes; consequently, hypoxia can have very different effects according to the salinity of the environment. The purpose of this study was to determine what role aquaporins play in water permeability in the hypoxia-tolerant euryhaline common killifish (Fundulus heteroclitus), in two important osmoregulatory organs-the gills and intestine. Using immunofluorescence, we localized aquaporin-3 (AQP3) protein to the basolateral and apical membranes of ionocytes and enterocytes, respectively. Although hypoxia increased branchial AQP3 messenger-RNA expression in seawater and freshwater, protein abundance did not correlate. Indeed, hypoxia did not alter AQP3 protein abundance in seawater and reduced it in the cell membranes of freshwater gills. Together, these observations suggest killifish AQP3 contributes to reduced diffusive water flux during hypoxia and normoxic recovery in freshwater and facilitates intestinal permeability in seawater and freshwater.

Resilience to multiple stressors in an aquatic plant and its microbiome

PREMISE

Outcomes of species interactions, especially mutualisms, are notoriously dependent on environmental context, and environments are changing rapidly. Studies have investigated how mutualisms respond to or ameliorate anthropogenic environmental changes, but most have focused on nutrient pollution or climate change and tested stressors one at a time. Relatively little is known about how mutualisms may be altered by or buffer the effects of multiple chemical contaminants, which differ fundamentally from nutrient or climate stressors and are especially widespread in aquatic habitats.

METHODS

We investigated the impacts of two contaminants on interactions between the duckweed Lemna minor and its microbiome. Sodium chloride (salt) and benzotriazole (a corrosion inhibitor) often co-occur in runoff to water bodies where duckweeds reside. We tested three L. minor genotypes with and without the culturable portion of their microbiome across field-realistic gradients of salt (3 levels) and benzotriazole (4 levels) in a fully factorial experiment (24 treatments, tested on each genotype) and measured plant and microbial growth.

RESULTS

Stressors had conditional effects. Salt decreased both plant and microbial growth and decreased plant survival more as benzotriazole concentrations increased. In contrast, benzotriazole did not affect microbial abundance and even benefited plants when salt and microbes were absent, perhaps due to biotransformation into growth-promoting compounds. Microbes did not ameliorate duckweed stressors; microbial inoculation increased plant growth, but not at high salt concentrations.

CONCLUSIONS

Our results suggest that multiple stressors matter when predicting responses of mutualisms to global change and that beneficial microbes may not always buffer hosts against stress.

Histopathology and immunohistochemistry of gills of Van fish (Alburnus tarichi Güldenstädt, 1814) infected with myxosporean parasites

The Van fish are a cyprinid species endemic to Turkey's largest soda lake, Lake Van, and have great economic value because they are a food source. Once a year, the fish take part in reproductive migration to the fresh waters flowing into the lake. The fish migrate from an extreme environment with high salinity (2.2%) and high pH (9.8). These fish are unable to reproduce in this alkaline environment and must migrate to fresh water during their breeding season. The aim of the present study is to report the presence of the myxosporean parasites on the gills and the pathological changes. Changes in gill histopathology, mucocytes, mitochondria-rich cells, expression of Heat Shock Protein 70 (Hsp70), and ATPase (NKA) were observed in the gill tissue. As a result of the histopathological changes in gills, infected fish had abundant plasmodia with different sizes. Plasmodia were found on gill filaments inside white ovoid-shaped structures. It was observed that plasmodia were contained on the primary filament which changed the histological structure of the gill tissue to a large extent. It was determined that the density and size of mucocytes in the infected areas of the gill tissue increased, whereas the number of mitochondria-rich cells decreased. Hsp70, an indicator of stress, was not different between normal and infected fish.

Draft genome assembly of Tenualosa ilisha, Hilsa shad, provides resource for osmoregulation studies

This study provides the first high-quality draft genome assembly (762.5 Mb) of Tenualosa ilisha that is highly contiguous and nearly complete. We observed a total of 2,864 contigs, with 96.4% completeness with N₅₀ of 2.65 Mbp and the largest contig length of 17.4 Mbp, along with a complete mitochondrial genome of 16,745 bases. A total number of 33,042 protein coding genes were predicted, among these, 512 genes were classified under 61 Gene Ontology (GO) terms, associated with various homeostasis processes. Highest number of genes belongs to cellular calcium ion homeostasis, followed by tissue homeostasis. A total of 97 genes were identified, with 16 GO terms related to water homeostasis. Claudins, Aquaporins, Connexins/Gap junctions, Adenylate cyclase, Solute carriers and Voltage gated potassium channel genes were observed to be higher in number in T. ilisha, as compared to that in other teleost species. Seven novel gene variants, in addition to claudin gene (CLDZ), were found in T. ilisha. The present study also identified two putative novel genes, NKAIN3 and L4AM1, for the first time in fish, for which further studies are required for pinpointing their functions in fish. In addition, 1.6 million simple sequence repeats were mined from draft genome assembly. The study provides a valuable genomic resource for the anadromous Hilsa. It will form a basis for future studies, pertaining to its adaptation mechanisms to different salinity levels during migration, which in turn would facilitate in its domestication.

The osmorespiratory compromise in the euryhaline killifish: water regulation during hypoxia

Freshwater- and seawater-acclimated Fundulus heteroclitus were exposed to acute hypoxia (10% air saturation, 3 h), followed by normoxic recovery (3 h). In both salinities, ventilation increased and heart rate fell in the classic manner, while Ṁ _O2  initially declined by ∼50%, with partial restoration by 3 h of hypoxia, and no O₂ debt repayment during recovery. Gill paracellular permeability (measured with [¹⁴C] PEG-4000) was 1.4-fold higher in seawater, and declined by 50% during hypoxia with post-exposure overshoot to 188%. A similar pattern with smaller changes occurred in freshwater. Drinking rate (also measured with [¹⁴C] PEG-4000) was 8-fold higher in seawater fish, but declined by ∼90% during hypoxia in both groups, with post-exposure overshoots to ∼270%. Gill diffusive water flux (measured with ³H₂O) was 1.9-fold higher in freshwater fish, and exhibited a ∼35% decrease during hypoxia, which persisted throughout recovery, but was unchanged during hypoxia in seawater fish. Nevertheless, freshwater killifish gained mass while seawater fish lost mass during hypoxia, and these changes were not corrected during normoxic recovery. We conclude that this hypoxia-tolerant teleost beneficially reduces gill water permeability in a salinity-dependent fashion during hypoxia, despite attempting to simultaneously improve Ṁ _O2 , but nevertheless incurs a net water balance penalty in both freshwater and seawater.

Transcriptome analyses reveal alterations in muscle metabolism, immune responses and reproductive behavior of Japanese mantis shrimp (Oratosquilla oratoria) at different cold temperature

Low temperature reduction is thought to cause widespread effects on the physical and behavioral traits of marine organisms, which include metabolic processes, immune responses, and reproductive behavior. Crustaceans are generally considered sensitive to temperature reduction due to the lack of efficient regulators. To better understand the molecular regulatory mechanisms of crustacean exposure to cold stress, Japanese mantis shrimp (Oratosquilla oratoria) was chosen as a representative crustacean. Transcriptomic responses in O. oratoria from five temperatures (25 °C, 22 °C, 19 °C, 16 °C, and 13 °C) were studied using RNA-seq. A total of 64.91 Gb of clean transcriptomic data were generated in 10 libraries and then spliced into 52,107 unigenes with an average length of 1089 bp and an N50 length of 1872 bp. A total of 14,841 unigenes was annotated in at least one database using Blastx alignment. Compared with the control temperature (25 °C), 7, 21, 58, and 236 unigenes were significantly differentially expressed at 22 °C, 19 °C, 16 °C, and 13 °C, respectively. GO analysis showed that 6, 20, 27, and 35 terms were significantly enriched at 22 °C, 19 °C, 16 °C, and 13 °C, respectively. In addition, 2, 5, 2, and 10 significant pathways were presented at 22 °C, 19 °C, 16 °C, and 13 °C, respectively. Combining NR, GO, and KEGG annotation information, many genes significantly differentially expressed at low temperatures may be associated with metabolic processes, immune response, and reproductive behavior. Additionally, we reconstructed the phylogenetic relationship based on 366 orthologous genes and the predicted differentiation time of O. oratoria and P. vannamei range from 212.82 to 365.30 Mya. Furthermore, 16 orthologous genes were identified as PSGs and 30 orthologous genes were identified as FEGs and these adaptive genes were associated with energy metabolism, stress response and immunity, and multiple cellular processing. These results provide fundamental information about molecular mechanisms regulating cold stress response of O. oratoria.

Comparative transcriptomics of 3 high-altitude passerine birds and their low-altitude relatives

High-altitude environments present strong stresses for living organisms, which have driven striking phenotypic and genetic adaptations. While previous studies have revealed multiple genetic adaptations in high-altitude species, how evolutionary history (i.e., phylogenetic background) contributes to similarity in genetic adaptations to high-altitude environments is largely unknown, in particular in a group of birds. We explored this in 3 high-altitude passerine birds from the Qinghai-Tibet Plateau and their low-altitude relatives in lowland eastern China. We generated transcriptomic data for 5 tissues across these species and compared sequence changes and expression shifts between high- and low-altitude pairs. Sequence comparison revealed that similarity in all 3 high-altitude species was high for genes under positive selection (218 genes) but low in amino acid substitutions (only 4 genes sharing identical amino acid substitutions). Expression profiles for all genes identified a tissue-specific expression pattern (i.e., all species clustered by tissue). By contrast, an altitude-related pattern was observed in genes differentially expressed between all 3 species pairs and genes associated with altitude, suggesting that the high-altitude environment may drive similar expression shifts in the 3 high-altitude species. Gene expression level, gene connectivity, and the interactions of these 2 factors with altitude were correlated with evolutionary rates. Our results provide evidence for how gene sequence changes and expression shifts work in a concerted way in a group of high-altitude birds, leading to similar evolution routes in response to high-altitude environmental stresses.

Evolutionary Toxicogenomics of the Striped Killifish (Fundulus majalis) in the New Bedford Harbor (Massachusetts, USA)

In this paper, we used a Genotyping-by-Sequencing (GBS) approach to find and genotype more than 4000 genome-wide SNPs (Single Nucleotide Polymorphisms) from striped killifish exposed to a variety of polychlorinated biphenyls (PCBs) and other aromatic pollutants in New Bedford Harbor (NBH, Massachusetts, USA). The aims of this study were to identify the genetic consequences of exposure to aquatic pollutants and detect genes that may be under selection. Low genetic diversity (H_E and π) was found in the site exposed to the highest pollution level, but the pattern of genetic diversity did not match the pollution levels. Extensive connectivity was detected among sampling sites, which suggests that balanced gene flow may explain the lack of genetic variation in response to pollution levels. Tests for selection identified 539 candidate outliers, but many of the candidate outliers were not shared among tests. Differences among test results likely reflect different test assumptions and the complex pollutant mixture. Potentially, selectively important loci are associated with 151 SNPs, and enrichment analysis suggests a likely involvement of these genes with pollutants that occur in NBH. This result suggests that selective processes at genes targeted by pollutants may be occurring, even at a small geographical scale, and may allow the local striped killifish to resist the high pollution levels.

Netting the Stress Responses in Fish

In the last decade, the concept of animal stress has been stressed thin to accommodate the effects of short-term changes in cell and tissue physiology, major behavioral syndromes in individuals and ecological disturbances in populations. Seyle's definition of stress as "the nonspecific (common) result of any demand upon the body" now encompasses homeostasis in a broader sense, including all the hierarchical levels in a networked biological system. The heterogeneity of stress responses thus varies within individuals, and stressors become multimodal in terms of typology, source and effects, as well as the responses that each individual elicits to cope with the disturbance. In fish, the time course of changes after stress strongly depends on several factors, including the stressful experiences in early life, the vertical transmission of stressful-prone phenotypes, the degree of individual phenotypic plasticity, the robustness and variety of the epigenetic network related to environmentally induced changes, and the intrinsic behavioral responses (individuality/personality) of each individual. The hierarchical heterogeneity of stress responses demands a code that may decrypt and simplify the analysis of both proximate and evolutionary causes of a particular stress phenotype. We propose an analytical framework, the stressotope, defined as an adaptive scenario dominated by common environmental selective pressures that elicit common multilevel acute stress-induced responses and produce a measurable allostatic load in the organism. The stressotope may constitute a blueprint of embedded interactions between stress-related variations in cell states, molecular mediators and systemic networks, a map of circuits that reflect the inherited and acquired stress responses in an ever-changing, microorganismal-loaded medium. Several features of the proposed model are discussed as a starting point to pin down the maximum common stress responses across immune-neuroendocrine relevant physiological levels and scenarios, including the characterization of behavioral responses, in fish.

Dynamics of Gene Expression Responses for Ion Transport Proteins and Aquaporins in the Gill of a Euryhaline Pupfish during Freshwater and High-Salinity Acclimation

Pupfishes (genus Cyprinodon) evolved some of the broadest salinity tolerances of teleost fishes, with some taxa surviving in conditions from freshwater to nearly 160 ppt. In this study, we examined transcriptional dynamics of ion transporters and aquaporins in the gill of the desert Amargosa pupfish (Cyprinodon nevadensis amargosae) during rapid salinity change. Pupfish acclimated to 7.5 ppt were exposed to freshwater (0.3 ppt), seawater (35 ppt), or hypersaline (55 ppt) conditions over 4 h and sampled at these salinities over 14 d. Plasma osmolality and Cl^- concentration became elevated 8 h after the start of exposure to 35 or 55 ppt but returned to baseline levels after 14 d. Osmolality recovery was paralleled by increased gill Na⁺/K⁺-ATPase activity and higher relative levels of messenger RNAs (mRNAs) encoding cystic fibrosis transmembrane conductance regulator (cftr) and Na⁺/K⁺/2Cl^- cotransporter-1 (nkcc1). Transcripts encoding one Na⁺-HCO₃^- cotransporter-1 isoform (nbce1.1) also increased in the gills at higher salinities, while a second isoform (nbce1.2) increased expression in freshwater. Pupfish in freshwater also had lower osmolality and elevated gill mRNAs for Na⁺/H⁺ exchanger isoform-2a (nhe2a) and V-type H⁺-ATPase within 8 h, followed by increases in Na⁺/H⁺ exchanger-3 (nhe3), carbonic anhydrase 2 (ca2), and aquaporin-3 (aqp3) within 1 d. Gill mRNAs for Na⁺/Cl^- cotransporter-2 (ncc2) also were elevated 14 d after exposure to 0.3 ppt. These results offer insights into how coordinated transcriptional responses for ion transporters in the gill facilitate reestablishment of osmotic homeostasis after changes in environmental salinity and provide evidence that the teleost gill expresses two Na⁺-HCO₃^- cotransporter-1 isoforms with different roles in freshwater and seawater acclimation.

Adaptation to environmental stress at different timescales

Environments are changing rapidly, and to cope with these changes, organisms have to adapt. Adaptation can take many shapes and occur at different speeds, depending on the type of response, the trait, the population, and the environmental conditions. The biodiversity crisis that we are currently facing illustrates that numerous species and populations are not capable of adapting with sufficient speed to ongoing environmental changes. Here, we discuss current knowledge on the ability of animals and plants to adapt to environmental stress on different timescales, mainly focusing on thermal stress and ectotherms. We discuss within-generation responses that can be fast and induced within minutes or hours, evolutionary adaptations that are often slow and take several generations, and mechanisms that lay somewhere in between and that include epigenetic transgenerational effects. To understand and predict the impacts of environmental change and stress on biodiversity, we suggest that future studies should (1) have an increased focus on understanding the type and speed of responses to fast environmental changes; (2) focus on the importance of environmental fluctuations and the predictability of environmental conditions on adaptive capabilities, preferably in field studies encompassing several fitness components; and (3) look at ecosystem responses to environmental stress and their resilience when disturbed.

Metal accumulation and expression of genes encoding for metallothionein and copper transporters in a chronically exposed wild population of the fish Hyphessobrycon luetkenii

In the present study, metal (As, Cd, Cu, Fe, Mn, Pb and Zn) accumulation and expression of genes involved in metal metabolism (metallothioneins, ATP7A and CTR1) were evaluated in gills and muscle of the fish Hyphessobrycon luetkenii living in the João Dias creek, a site historically (~1870-1996) contaminated with a metal mixture associated with copper mining (Minas do Camaquã, southern Brazil). Fish were collected in a metal impacted site of the João Dias creek and kept in a cage at this site (PP fish) or translocated to a non-metal impacted reach of this creek (PC fish). Gill metal concentrations and metallothionein gene expression were lower in PC fish than in PP fish at any experimental time (24, 48 and 72 h). In muscle, no significant changes were observed. These findings indicate that metal accumulated in gills of wild fish chronically exposed to the metal mixture are more easily excreted than those accumulated in the muscle. In this case, expression of gene encoding for metallothionein is shown to play a key role in the regulation of metal accumulation in gills of H. luetkenii living in an area historically contaminated with a metal mixture associated with copper mining.

Arsenic Reduces Gene Expression Response to Changing Salinity in Killifish

Toxicogenomic approaches can detect and classify adverse interactions between environmental toxicants and other environmental stressors but require more complex experimental designs and analytical approaches. Here we use novel toxicogenomic techniques to analyze the effect of arsenic exposure in wild killifish populations acclimating to changing salinity. Fish from three populations were acclimated to full strength seawater and transferred to fresh water for 1 or 24 h. Linear models of gene expression in gill tissue identified 31 genes that responded to osmotic shock at 1 h and 178 genes that responded at 24 h. Arsenic exposure (100 μg/L) diminished the responses (reaction norms) of these genes by 22% at 1 h ( p = 1.0 × 10^-6) and by 10% at 24 h ( p = 3.0 × 10^-10). Arsenic also significantly reduced gene coregulation in gene regulatory networks ( p = 0.002, paired Levene's test), and interactions between arsenic and salinity acclimation were uniformly antagonistic at the biological pathway level ( p < 0.05, binomial test). Arsenic's systematic interference with gene expression reaction norms was validated in a mouse multistressor experiment, demonstrating the ability of these toxicogenomic approaches to identify biologically relevant adverse interactions between environmental toxicants and other environmental stressors.

Expression of genes related to metal metabolism in the freshwater fish Hyphessobrycon luetkenii living in a historically contaminated area associated with copper mining

Copper (Cu) mining in Minas do Camaquã-Brazil, released significant amounts of metals into the João Dias creek, where Hyphessobrycon luetkenii inhabit. Because the involvement of Cu in biological processes its concentration and availability is regulated by molecules as the metal regulatory transcription factor (MTF-1), metallothionein (MT) and transporters (ATP7A and CTR1). These genes were whole sequenced and their expression (GE) evaluated in gills, liver and intestine. Were collected fish in non-contaminated and contaminated (Cu 3.4-fold higher) sites of the creek (CC and PP) and respectively translocated (CP and PC) for 96 h. The GE of the non-translocated groups evidenced that MT, MTF-1 and CTR1 have organ specific differences between both communities. Additionally the translocation allowed to identify organ specific changes associated with the activation/inactivation of protective mechanisms. These findings indicate that MTF-1, MT and CTR-1 GE play an important role in the tolerance of H. luetkenii to Cu contamination.

The utility of transcriptomics in fish conservation

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

Evidence of the Complexity of Gene Expression Analysis in Fish Wild Populations

The present work examines the induction of the band 3 anion transport protein, mitogen-activated protein kinase, and lactate dehydrogenase, respectively related to osmolyte transport, cell volume regulation, and energy production in the gills of two tilapia strains exposed to either freshwater or hypersaline water. Overall, genes showed similar expression patterns between strains. However, a wild population survey across a range of natural habitats and salinities did not reveal the expected patterns. Although significant, the correlations between gene expression and salinity were slightly ambiguous and did not show any link with phenotypic differences in life history traits previously reported between the same populations. The differential expression was also not associated with the population genetic structure inferred from neutral markers. The results suggest that the differential expression observed is not the result of evolutionary forces such as genetic drift or adaptation by natural selection. Instead, it can be speculated that genes responded to various abiotic and biotic stressors, including factors intrinsic to animals. This study provides clear evidence of the complexity of gene expression analysis in wild populations and shows that more attention needs to be paid when selecting candidates as potential biomarkers for monitoring adaptive responses to a specific environmental perturbation.

Transcriptomic differentiation underlying marine-to-freshwater transitions in the South American silversides Odontesthes argentinensis and O. bonariensis (Atheriniformes)

Salinity gradients are critical habitat determinants for freshwater organisms. Silverside fishes in the genus Odontesthes have recently and repeatedly transitioned from marine to freshwater habitats, overcoming a strong ecological barrier. Genomic and transcriptomic changes involved in this kind of transition are only known for a few model species. We present new data and analyses of gene expression and microbiome composition in the gills of two closely related silverside species, marine O. argentinensis and freshwater O. bonariensis and find more than three thousand transcripts differentially expressed, with osmoregulatory/ion transport genes and immune genes showing very different expression patterns across species. Interspecific differences also involve more than one thousand transcripts with nonsynonymous SNPs in the coding sequences, most of which were not differentially expressed. In addition to characterizing gill transcriptomes from wild‐caught marine and freshwater fishes, we test experimentally the response to salinity increases by O. bonariensis collected from freshwater habitats. Patterns of expression in gill transcriptomes of O. bonariensis exposed to high salinity do not resemble O. argentinensis mRNA expression, suggesting lack of plasticity for adaptation to marine conditions in this species. The diversity of functions associated with both the differentially expressed set of transcripts and those with sequence divergence plus marked microbiome differences suggest that multiple abiotic and biotic factors in marine and freshwater habitats are driving transcriptomic differences between these species.

Long-Term Overgrazing-Induced Memory Decreases Photosynthesis of Clonal Offspring in a Perennial Grassland Plant

Previous studies of transgenerational plasticity have demonstrated that long-term overgrazing experienced by Leymus chinensis, an ecologically dominant, rhizomatous grass species in eastern Eurasian temperate grassland, significantly affects its clonal growth in subsequent generations. However, there is a dearth of information on the reasons underlying this overgrazing-induced memory effect in plant morphological plasticity. We characterized the relationship between a dwarf phenotype and photosynthesis function decline of L. chinensis from the perspective of leaf photosynthesis by using both field measurement and rhizome buds culture cultivated in a greenhouse. Leaf photosynthetic functions (net photosynthetic rate, stomatal conductance, intercellular carbon dioxide concentration, and transpiration rate) were significantly decreased in smaller L. chinensis individuals that were induced to have a dwarf phenotype by being heavily grazed in the field. This decreased photosynthetic function was maintained a generation after greenhouse tests in which grazing was excluded. Both the response of L. chinensis morphological traits and photosynthetic functions in greenhouse were deceased relative to those in the field experiment. Further, there were significant decreases in leaf chlorophyll content and Rubisco enzyme activities of leaves between bud-cultured dwarf and non-dwarf L. chinensis in the greenhouse. Moreover, gene expression patterns showed that the bud-cultured dwarf L. chinensis significantly down-regulated (by 1.86- to 5.33-fold) a series of key genes that regulate photosynthetic efficiency, stomata opening, and chloroplast development compared with the non-dwarf L. chinensis. This is among the first studies revealing a linkage between long-term overgrazing affecting the transgenerational morphological plasticity of clonal plants and physiologically adaptive photosynthesis function. Overall, clonal transgenerational effects in L. chinensis phenotypic traits heavily involve photosynthetic plasticity.

Identifying selectively important amino acid positions associated with alternative habitat environments in fish mitochondrial genomes

Fish species inhabitating seawater (SW) or freshwater (FW) habitats have to develop genetic adaptations to alternative environment factors, especially salinity. Functional consequences of the protein variations associated with habitat environments in fish mitochondrial genomes have not yet received much attention. We analyzed 829 complete fish mitochondrial genomes and compared the amino acid differences of 13 mitochondrial protein families between FW and SW fish groups. We identified 47 specificity determining sites (SDS) that associated with FW or SW environments from 12 mitochondrial protein families. Thirty-two (68%) of the SDS sites are hydrophobic, 13 (28%) are neutral, and the remaining sites are acidic or basic. Seven of those SDS from ND1, ND2 and ND5 were scored as probably damaging to the protein structures. Furthermore, phylogenetic tree based Bayes Empirical Bayes analysis also detected 63 positive sites associated with alternative habitat environments across ten mtDNA proteins. These signatures could be important for studying mitochondrial genetic variation relevant to fish physiology and ecology.

The landscape of extreme genomic variation in the highly adaptable Atlantic killifish

Understanding and predicting the fate of populations in changing environments require knowledge about the mechanisms that support phenotypic plasticity and the adaptive value and evolutionary fate of genetic variation within populations. Atlantic killifish (Fundulus heteroclitus) exhibit extensive phenotypic plasticity that supports large population sizes in highly fluctuating estuarine environments. Populations have also evolved diverse local adaptations. To yield insights into the genomic variation that supports their adaptability, we sequenced a reference genome and 48 additional whole genomes from a wild population. Evolution of genes associated with cell cycle regulation and apoptosis is accelerated along the killifish lineage, which is likely tied to adaptations for life in highly variable estuarine environments. Genome-wide standing genetic variation, including nucleotide diversity and copy number variation, is extremely high. The highest diversity genes are those associated with immune function and olfaction, whereas genes under greatest evolutionary constraint are those associated with neurological, developmental, and cytoskeletal functions. Reduced genetic variation is detected for tight junction proteins, which in killifish regulate paracellular permeability that supports their extreme physiological flexibility. Low-diversity genes engage in more regulatory interactions than high-diversity genes, consistent with the influence of pleiotropic constraint on molecular evolution. High genetic variation is crucial for continued persistence of species given the pace of contemporary environmental change. Killifish populations harbor among the highest levels of nucleotide diversity yet reported for a vertebrate species, and thus may serve as a useful model system for studying evolutionary potential in variable and changing environments.

Is Exposure to Macondo Oil Reflected in the Otolith Chemistry of Marsh-Resident Fish?

Genomic and physiological responses in Gulf killifish (Fundulus grandis) in the northern Gulf of Mexico have confirmed oil exposure of resident marsh fish following the Macondo blowout in 2010. Using these same fish, we evaluated otolith microchemistry as a method for assessing oil exposure history. Laser-ablation inductively-coupled-plasma mass spectrometry was used to analyze the chemical composition of sagittal otoliths to assess whether a trace metal signature could be detected in the otoliths of F. grandis collected from a Macondo-oil impacted site in 2010, post-spill relative to pre-spill, as well as versus fish from areas not impacted by the spill. We found no evidence of increased concentrations of two elements associated with oil contamination (nickel and vanadium) in F. grandis otoliths regardless of Macondo oil exposure history. One potential explanation for this is that Macondo oil is relatively depleted of those metals compared to other crude oils globally. During and after the spill, however, elevated levels of barium, lead, and to a lesser degree, copper were detected in killifish otoliths at the oil-impacted collection site in coastal Louisiana. This may reflect oil contact or other environmental perturbations that occurred concomitant with oiling. For example, increases in barium in otoliths from oil-exposed fish followed (temporally) freshwater diversions in Louisiana in 2010. This implicates (but does not conclusively demonstrate) freshwater diversions from the Mississippi River (with previously recorded higher concentrations of lead and copper), designed to halt the ingress of oil, as a mechanism for elevated elemental uptake in otoliths of Louisiana marsh fishes. These results highlight the potentially complex and indirect effects of the Macondo oil spill and human responses to it on Gulf of Mexico ecosystems, and emphasize the need to consider the multiple stressors acting simultaneously on inshore fish communities.

Hormonal regulation of aquaporin 3: opposing actions of prolactin and cortisol in tilapia gill

Aquaporins (Aqps) are expressed within key osmoregulatory tissues where they mediate the movement of water and selected solutes across cell membranes. We leveraged the functional plasticity of Mozambique tilapia (Oreochromis mossambicus) gill epithelium to examine how Aqp3, an aquaglyceroporin, is regulated in response to osmoregulatory demands. Particular attention was paid to the actions of critical osmoregulatory hormones, namely, prolactin (Prl), growth hormone and cortisol. Branchial aqp3 mRNA levels were modulated following changes in environmental salinity, with enhanced aqp3 mRNA expression upon transfer from seawater to freshwater (FW). Accordingly, extensive Aqp3 immunoreactivity was localized to cell membranes of branchial epithelium in FW-acclimated animals. Upon transferring hypophysectomized tilapia to FW, we identified that a pituitary factor(s) is required for Aqp3 expression in FW. Replacement with ovine Prl (oPrl) was sufficient to stimulate Aqp3 expression in hypophysectomized animals held in FW, an effect blocked by coinjection with cortisol. Both oPrl and native tilapia Prls (tPrl177 and tPrl188) stimulated aqp3 in incubated gill filaments in a concentration-related manner. Consistent with in vivo responses, coincubation with cortisol blocked oPrl-stimulated aqp3 expression in vitro Our data indicate that Prl and cortisol act directly upon branchial epithelium to regulate Aqp3 in tilapia. Thus, within the context of the diverse actions of Prl on hydromineral balance in vertebrates, we define a new role for Prl as a regulator of Aqp expression.

Transcriptional response of yellow perch to changes in ambient metal concentrations-A reciprocal field transplantation experiment

Recent local adaptation to pollution has been evidenced in several organisms inhabiting environments heavily contaminated by metals. Nevertheless, the molecular mechanisms underlying adaptation to high metal concentrations are poorly understood, especially in fishes. Yellow perch (Perca flavescens) populations from lakes in the mining area of Rouyn-Noranda (QC, Canada) have been faced with metal contamination for about 90 years. Here, we examine gene transcription patterns of fish reciprocally transplanted between a reference and a metal-contaminated lake and also fish caged in their native lake. After four weeks, 111 genes were differentially transcribed in metal-naïve fish transferred to the metal-contaminated lake, revealing a plastic response to metal exposure. Genes involved in the citric cycle and beta-oxidation pathways were under-transcribed, suggesting a potential strategy to mitigate the effects of metal stress by reducing energy turnover. However, metal-contaminated fish transplanted to the reference lake did not show any transcriptomic response, indicating a reduced plastic response capability to sudden reduction in metal concentrations. Moreover, the transcription of other genes, especially ones involved in energy metabolism, was affected by caging. Overall, our results highlight environmental stress response mechanisms in yellow perch at the transcriptomic level and support a rapid adaptive response to metal exposure through genetic assimilation.

Cytogenetic variation of repetitive DNA elements in Hoplias malabaricus (Characiformes - Erythrinidae) from white, black and clear water rivers of the Amazon basin

Hoplias malabaricus is a common fish species occurring in white, black and clear water rivers of the Amazon basin. Its large distribution across distinct aquatic environments can pose stressful conditions for dispersal and creates possibilities for the emergence of local adaptive profiles. We investigated the chromosomal localization of repetitive DNA markers (constitutive heterochromatin, rDNA and the transposable element REX-3) in populations from the Amazonas river (white water), the Negro river (black water) and the Tapajós river (clear water), in order to address the variation/association of cytogenomic features and environmental conditions. We found a conserved karyotypic macrostructure with a diploid number of 40 chromosomes (20 metacentrics + 20 submetacentrics) in all the samples. Heteromorphism in pair 14 was detected as evidence for the initial differentiation of an XX/XY system. Minor differences detected in the amount of repetitive DNA markers are interpreted as possible signatures of local adaptations to distinct aquatic environments.

Effects of Anthropogenic Pollution on the Oxidative Phosphorylation Pathway of Hepatocytes from Natural Populations of Fundulus heteroclitus

Persistent organic pollutants (POPs), including polycyclic aromatic hydrocarbons (PAHs) and polychlorinated biphenyls (PCBs), potentially target mitochondria and cause toxicity. We compared the effects of POPs on mitochondrial respiration by measuring oxidative phosphorylation (OxPhos) metabolism in hepatocytes isolated from lab-depurated Fundulus heteroclitus from a Superfund site contaminated with PAHs (Elizabeth River VA, USA) relative to OxPhos metabolism in individuals from a relatively clean, reference population (King's Creek VA, USA). In individuals from the polluted Elizabeth River population, OxPhos metabolism displayed lower LEAK and lower activities in complex III, complex IV, and E State, but higher activity in complex I compared to individuals from the reference King's Creek population. To test the supposition that these differences were due to or related to the chronic PAH contamination history of the Elizabeth River population, we compared the OxPhos functions of undosed individuals from the polluted and reference populations to individuals from these populations dosed with a PAH {benzo [α] pyrene (BaP)} or a PCB {PCB126 (3,3',4,4',5-pentachlorobiphenyl)}, respectively. Exposure to PAH or PCB affected OxPhos in the reference King's Creek population but had no detectable effects on the polluted Elizabeth River population. Thus, PAH exposure significantly increased LEAK, and exposure to PCB126 significantly decreased State 3, E state and complex I activity in the reference King's Creek population. These data strongly implicate an evolved tolerance in the Elizabeth River fish where dosed fish are not affected by PAH exposure and undosed fish show decreased LEAK and increased State 3 and E state.

Initial Evidence for Adaptive Selection on the NADH Subunit Two of Freshwater Dolphins by Analyses of Mitochondrial Genomes

A small number of cetaceans have adapted to an entirely freshwater environment, having colonized rivers in Asia and South America from an ancestral origin in the marine environment. This includes the ‘river dolphins’, early divergence from the odontocete lineage, and two species of true dolphins (Family Delphinidae). Successful adaptation to the freshwater environment may have required increased demands in energy involved in processes such as the mitochondrial osmotic balance. For this reason, riverine odontocetes provide a compelling natural experiment in adaptation of mammals from marine to freshwater habitats. Here we present initial evidence of positive selection in the NADH dehydrogenase subunit 2 of riverine odontocetes by analyses of full mitochondrial genomes, using tests of selection and protein structure modeling. The codon model with highest statistical support corresponds to three discrete categories for amino acid sites, those under positive, neutral, and purifying selection. With this model we found positive selection at site 297 of the NADH dehydrogenase subunit 2 (d_N/d_S>1.0,) leading to a substitution of an Ala or Val from the ancestral state of Thr. A phylogenetic reconstruction of 27 cetacean mitogenomes showed that an Ala substitution has evolved at least four times in cetaceans, once or more in the three ‘river dolphins’ (Families Pontoporidae, Lipotidae and Inidae), once in the riverine Sotalia fluviatilis (but not in its marine sister taxa), once in the riverine Orcaella brevirostris from the Mekong River (but not in its marine sister taxa) and once in two other related marine dolphins. We located the position of this amino acid substitution in an alpha-helix channel in the trans-membrane domain in both the E. coli structure and Sotalia fluviatilis model. In E. coli this position is located in a helix implicated in a proton translocation channel of respiratory complex 1 and may have a similar role in the NADH dehydrogenases of cetaceans.

Comparative renal gene expression in response to abrupt hypoosmotic shock in spotted scat (Scatophagus argus)

Scatophagus argus, a euryhaline fish, is notable for its ability to tolerate a wide range of environmental salinities and especially for its tolerance to a rapid, marked reduction in salinity. Therefore, S. argus is a good model for studying the molecular mechanisms mediating abrupt hyperosmoregulation. The serum osmotic pressure decreased steeply within one hour after transferring S. argus from seawater (SW) to freshwater (FW) and remained at new balance throughout the duration of one week. To explain this phenomenon and understand the molecular responses to an abrupt hypoosmotic shock, hypoosmotic stress responsive genes were identified by constructing two suppression subtractive hybridization (SSH) cDNA libraries from the kidneys of S. argus that had been transferred from SW to FW. After trimming and blasting, 52 ESTs were picked out from the subtractive library. Among them, 11 genes were significantly up-regulated (p < 0.05). The kinetics studies of gene expression levels were conducted for 1 week after the transfer using quantitative real-time PCR. A significant variation in the expression of these genes occurred within 12h after the hypoosmotic shock, except for growth hormone (GH) and polyadenylate binding protein 1 (PBP1), which were significantly up-regulated 2 days post-transfer. Our results suggest different functional roles for these genes in response to hypoosmotic stress during the stress response phase (1 hpt-12 hpt) and stable phase (12 hpt-7 dpt). Furthermore, the plasma growth hormone level was detected to be significantly elevated at 1 hpt and 24 hpt following abrupt hypoosmotic shock. Meanwhile, several hematological parameters, hemoglobin (HGB), red blood cell (RBC) and mean cellular hemoglobin concentration (MCHC), were observed to be significantly increased at 12 hpt and 2 dpt compared with that of control group. Our results provide a solid basis from which to conduct future studies on the osmoregulatory mechanisms in the euryhaline fish.

Paracellular pathway remodeling enhances sodium secretion by teleost fish in hypersaline environments

In vertebrate salt-secreting epithelia, Na(+) moves passively down an electrochemical gradient via a paracellular pathway. We assessed how this pathway is modified to allow Na(+) secretion in hypersaline environments. Mummichogs (Fundulus heteroclitus) acclimated to hypersaline [2× seawater (2SW), 64‰] for 30 days developed invasive projections of accessory cells with an increased area of tight junctions, detected by punctate distribution of CFTR (cystic fibrosis transmembrane conductance regulator) immunofluorescence and transmission electron miscroscopy of the opercular epithelia, which form a gill-like tissue rich in ionocytes. Distribution of CFTR was not explained by membrane raft organization, because chlorpromazine (50 μmol l(-1)) and filipin (1.5 μmol l(-1)) did not affect opercular epithelia electrophysiology. Isolated opercular epithelia bathed in SW on the mucosal side had a transepithelial potential (Vt) of +40.1±0.9 mV (N=24), sufficient for passive Na(+) secretion (Nernst equilibrium voltage≡ENa=+24.11 mV). Opercular epithelia from fish acclimated to 2SW and bathed in 2SW had higher Vt of +45.1±1.2 mV (N=24), sufficient for passive Na(+) secretion (ENa=+40.74 mV), but with diminished net driving force. Bumetanide block of Cl(-) secretion reduced Vt by 45% and 29% in SW and 2SW, respectively, a decrease in the driving force for Na(+) extrusion. Estimates of shunt conductance from epithelial conductance (Gt) versus short-circuit current (Isc) plots (extrapolation to zero Isc) suggested a reduction in total epithelial shunt conductance in 2SW-acclimated fish. In contrast, the morphological elaboration of tight junctions, leading to an increase in accessory-cell-ionocyte contact points, suggests an increase in local paracellular conductance, compensating for the diminished net driving force for Na(+) and allowing salt secretion, even in extreme salinities.

Transcriptomic response to heat stress among ecologically divergent populations of redband trout

Background

As ectothermic organisms have evolved to differing aquatic climates, the molecular basis of thermal adaptation is a key area of research. In this study, we tested for differential transcriptional response of ecologically divergent populations of redband trout (Oncorhynchus mykiss gairdneri) that have evolved in desert and montane climates. Each pure strain and their F1 cross were reared in a common garden environment and exposed over four weeks to diel water temperatures that were similar to those experienced in desert climates within the species’ range. Gill tissues were collected from the three strains of fish (desert, montane, F1 crosses) at the peak of heat stress and tested for mRNA expression differences across the transcriptome with RNA-seq.

Results

Strong differences in transcriptomic response to heat stress were observed across strains confirming that fish from desert environments have evolved diverse mechanisms to cope with stressful environments. As expected, a large number of total transcripts (12,814) were differentially expressed in the study (FDR ≤ 0.05) with 2310 transcripts in common for all three strains, but the desert strain had a larger number of unique differentially expressed transcripts (2875) than the montane (1982) or the F1 (2355) strain. Strongly differentiated genes (>4 fold change and FDR ≤ 0.05) were particularly abundant in the desert strain (824 unique contigs) relative to the other two strains (montane = 58; F1 = 192).

Conclusions

This study demonstrated patterns of acclimation (i.e., phenotypic plasticity) within strains and evolutionary adaptation among strains in numerous genes throughout the transcriptome. Key stress response genes such as molecular chaperones (i.e., heat shock proteins) had adaptive patterns of gene expression among strains, but also a much higher number of metabolic and cellular process genes were differentially expressed in the desert strain demonstrating these biological pathways are critical for thermal adaptation to warm aquatic climates. The results of this study further elucidate the molecular basis for thermal adaptation in aquatic ecosystems and extend the potential for identifying genes that may be critical for adaptation to changing climates.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1246-5) contains supplementary material, which is available to authorized users.

Genome-wide single-generation signatures of local selection in the panmictic European eel

Next-generation sequencing and the collection of genome-wide data allow identifying adaptive variation and footprints of directional selection. Using a large SNP data set from 259 RAD-sequenced European eel individuals (glass eels) from eight locations between 34 and 64(o) N, we examined the patterns of genome-wide genetic diversity across locations. We tested for local selection by searching for increased population differentiation using F(ST) -based outlier tests and by testing for significant associations between allele frequencies and environmental variables. The overall low genetic differentiation found (F(ST) = 0.0007) indicates that most of the genome is homogenized by gene flow, providing further evidence for genomic panmixia in the European eel. The lack of genetic substructuring was consistent at both nuclear and mitochondrial SNPs. Using an extensive number of diagnostic SNPs, results showed a low occurrence of hybrids between European and American eel, mainly limited to Iceland (5.9%), although individuals with signatures of introgression several generations back in time were found in mainland Europe. Despite panmixia, a small set of SNPs showed high genetic differentiation consistent with single-generation signatures of spatially varying selection acting on glass eels. After screening 50 354 SNPs, a total of 754 potentially locally selected SNPs were identified. Candidate genes for local selection constituted a wide array of functions, including calcium signalling, neuroactive ligand-receptor interaction and circadian rhythm. Remarkably, one of the candidate genes identified is PERIOD, possibly related to differences in local photoperiod associated with the >30° difference in latitude between locations. Genes under selection were spread across the genome, and there were no large regions of increased differentiation as expected when selection occurs within just a single generation due to panmixia. This supports the conclusion that most of the genome is homogenized by gene flow that removes any effects of diversifying selection from each new generation.

Transcriptomics profiling and steroid production in mummichog (Fundulus heteroclitus) testes after treatment with 5α-dihydrotestosterone

5α-Dihydrotestosterone (DHT) is a potent androgen in mammals with multiple roles; however the physiological actions of DHT in male fishes are not well known. To address this knowledge gap, male mummichog (Fundulus heteroclitus) were continuously exposed to 0, 5, and 50 μg/L DHT for 21 days. Following exposure, testes were separated for histology, ex vivo incubation to measure steroidogenic capacity, and gene expression analyses (real-time PCR and microarray). DHT significantly decreased ex vivo 11-ketotestosterone (11KT) production in males exposed to 50 μg/L DHT but not 5 μg/L DHT, and DHT exposure did not affect ex vivo testosterone production. Histological examination revealed that the amount of interlobular and connective tissue present in the testes was increased in the 50 μg/L DHT treatment. Despite reductions in the production of 11KT, DHT did not affect the expression of targeted genes in the steroidogenic pathway such as steroidogenic acute regulatory protein (star), P450 side chain cleavage (cyp11a1) and 11β-hydroxysteroid dehydrogenase (hsd11b3). Microarray analysis in the testes of individuals from control and 50 μg/L DHT revealed that males exposed to 50 μg/L DHT showed regulated transcriptional sub-networks that were related to immunity, regulation of blood flow, lipids and xenobiotic clearance, suggesting that DHT may be involved in the physiological regulation of these processes in the fish testes. A second objective of this study was to determine the feasibility of measuring mRNA levels in tissues used for ex vivo steroid production by comparing RNA integrity and transcript levels in testes of both immediately flash frozen tissue and incubated tissue. There was no significant difference in RNA quality between the two time points, indicating RNA integrity can remain intact for at least 18 h in ex vivo assays, thereby providing a viable option for researchers assessing multi-level biological reproductive endpoints when limited tissue is available. While the gene expression levels of actb, efla, rps12, rps18, star, and hsd11b3 remained unchanged, esr2a (esrba), esr2b (esrbb) and cyp11a1 were significantly lower in incubated tissue compared to flash frozen tissue. Therefore caution must be used as the steady-state levels of select genes may change over time. This study improves our understanding of DHT action in the teleostean testis and generates new hypotheses regarding cell processes that are regulated by this underexplored and potent androgen.

Common functional targets of adaptive micro- and macro-evolutionary divergence in killifish

Environmental salinity presents a key barrier to dispersal for most aquatic organisms, and adaptation to alternate osmotic environments likely enables species diversification. Little is known of the functional basis for derived tolerance to environmental salinity. We integrate comparative physiology and functional genomics to explore the mechanistic underpinnings of evolved variation in osmotic plasticity within and among two species of killifish; Fundulus majalis harbours the ancestral mainly salt-tolerant phenotype, whereas Fundulus heteroclitus harbours a derived physiology that retains extreme salt tolerance but with expanded osmotic plasticity towards the freshwater end of the osmotic continuum. Common-garden comparative hypo-osmotic challenge experiments show that F. heteroclitus is capable of remodelling gill epithelia more quickly and at more extreme osmotic challenge than F. majalis. We detect an unusual pattern of baseline transcriptome divergence, where neutral evolutionary processes appear to govern expression divergence within species, but patterns of divergence for these genes between species do not follow neutral expectations. During acclimation, genome expression profiling identifies mechanisms of acclimation-associated response that are conserved within the genus including regulation of paracellular permeability. In contrast, several responses vary among species including those putatively associated with cell volume regulation, and these same mechanisms are targets for adaptive physiological divergence along osmotic gradients within F. heteroclitus. As such, the genomic and physiological mechanisms that are associated with adaptive fine-tuning within species also contribute to macro-evolutionary divergence as species diversify across osmotic niches.

Genomics and the Evolution of Phenotypic Traits

Evolutionary genetics has entered an unprecedented era of discovery, catalyzed in large part by the development of technologies that provide information about genome sequence and function. An important benefit is the ability to move beyond a handful of model organisms in lab settings to identify the genetic basis for evolutionarily interesting traits in many organisms in natural settings. Other benefits are the abilities to identify causal mutations and validate their phenotypic consequences more readily and in many more species. Genomic technologies have reinvigorated interest in some of the most fundamental and persistent questions in evolutionary genetics, revealed previously unsuspected evolutionary phenomena, and opened the door to a wide range of new questions.