Whole-Genome Sequencing Reveals That Regulatory and Low Pleiotropy Variants Underlie Local Adaptation to Environmental Variability in Purple Sea Urchins

AbstractOrganisms experience environments that vary across both space and time. Such environmental heterogeneity shapes standing genetic variation and may influence species' capacity to adapt to rapid environmental change. However, we know little about the kind of genetic variation that is involved in local adaptation to environmental variability. To address this gap, we sequenced the whole genomes of 140 purple sea urchins (Strongylocentrotus purpuratus) from seven populations that vary in their degree of pH variability. Despite no evidence of global population structure, we found a suite of single-nucleotide polymorphisms (SNPs) tightly correlated with local pH variability (outlier SNPs), which were overrepresented in regions putatively involved in gene regulation (long noncoding RNA and enhancers), supporting the idea that variation in regulatory regions is important for local adaptation to variability. In addition, outliers in genes were found to be (i) enriched for biomineralization and ion homeostasis functions related to low pH response, (ii) less central to the protein-protein interaction network, and (iii) underrepresented among genes highly expressed during early development. Taken together, these results suggest that loci that underlie local adaptation to pH variability in purple sea urchins fall in regions with potentially low pleiotropic effects (based on analyses involving regulatory regions, network centrality, and expression time) involved in low pH response (based on functional enrichment).

Transcriptome Analysis Revealed Osmoregulation Related Regulatory Networks and Hub Genes in the Gills of Hilsa shad, Tenualosa ilisha, during the Migratory Osmotic Stress

Tenualosa ilisha (Hilsa shad), an anadromous fish, usually inhabits coastal and estuarine waters, and migrates to freshwater for spawning. In this study, large-scale gill transcriptome analyses from three salinity regions, i.e., fresh, brackish and marine water, revealed 3277 differentially expressed genes (DEGs), out of which 232 were found to be common between marine vs freshwater and brackish vs freshwater. These genes were mapped into 54 KEGG Pathways, and the most significant of these were focal adhesion, adherens junction, tight junction, and PI3K-Akt signaling pathways. A total of 24 osmoregulatory genes were found to be differentially expressed in different habitats. The gene members of slc16 and slc2 families showed a dissimilar pattern of expressions, while two claudin genes (cldn11 & cldn10), transmembrane tm56b, and voltage-gated potassium channel gene kcna10 were downregulated in freshwater samples, as compared to that of brackish and marine environment. Protein-protein interaction (PPI) network analysis of 232 DEGs showed 101 genes to be involved in PPI, while fn1 gene was found to be interacting with the highest number of genes (36). Twenty-five hub genes belonged to 12 functional groups, with muscle structure development with seven genes, forming the major group. These results provided valuable information about the genes, potentially involved in the molecular mechanisms regulating water homeostasis in gills, during migration for spawning and low-salinity adaptation in Hilsa shad. These genes may form the basis for the bio-marker development for adaptation to the stress levied by major environmental changes, due to hatchery/culture conditions.

Transcriptomic Analysis in Marine Medaka Gill Reveals That the Hypo-Osmotic Stress Could Alter the Immune Response via the IL17 Signaling Pathway

Fish gills are the major osmoregulatory tissue that contact the external water environment and have developed an effective osmoregulatory mechanism to maintain cellular function. Marine medaka (Oryzias melastigma) has the ability to live in both seawater and fresh water environments. The present study performed a seawater (SW) to 50% seawater (SFW) transfer, and the gill samples were used for comparative transcriptomic analysis to study the alteration of hypo-osmotic stress on immune responsive genes in this model organism. The result identified 518 differentiated expressed genes (DEGs) after the SW to SFW transfer. Various pathways such as p53 signaling, forkhead box O signaling, and the cell cycle were enriched. Moreover, the immune system was highlighted as one of the top altered biological processes in the enrichment analysis. Various cytokines, chemokines, and inflammatory genes that participate in the IL-17 signaling pathway were suppressed after the SW to SFW transfer. On the other hand, some immunoglobulin-related genes were up-regulated. The results were further validated by real-time qPCR. Taken together, our study provides additional gill transcriptome information in marine medaka; it also supports the notion that osmotic stress could influence the immune responses in fish gills.

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.

Loss of transcriptional plasticity but sustained adaptive capacity after adaptation to global change conditions in a marine copepod

Adaptive evolution and phenotypic plasticity will fuel resilience in the geologically unprecedented warming and acidification of the earth’s oceans, however, we have much to learn about the interactions and costs of these mechanisms of resilience. Here, using 20 generations of experimental evolution followed by three generations of reciprocal transplants, we investigated the relationship between adaptation and plasticity in the marine copepod, Acartia tonsa, in future global change conditions (high temperature and high CO₂). We found parallel adaptation to global change conditions in genes related to stress response, gene expression regulation, actin regulation, developmental processes, and energy production. However, reciprocal transplantation showed that adaptation resulted in a loss of transcriptional plasticity, reduced fecundity, and reduced population growth when global change-adapted animals were returned to ambient conditions or reared in low food conditions. However, after three successive transplant generations, global change-adapted animals were able to match the ambient-adaptive transcriptional profile. Concurrent changes in allele frequencies and erosion of nucleotide diversity suggest that this recovery occurred via adaptation back to ancestral conditions. These results demonstrate that while plasticity facilitated initial survival in global change conditions, it eroded after 20 generations as populations adapted, limiting resilience to new stressors and previously benign environments.

Rapid adaptation will facilitate species resilience under global climate change, but its effects on plasticity are less commonly investigated. This study shows that 20 generations of experimental adaptation in a marine copepod drives a rapid loss of plasticity that carries costs and might have impacts on future resilience to environmental change.

The time course of molecular acclimation to seawater in a euryhaline fish

The Arabian pupfish, Aphanius dispar, is a euryhaline fish inhabiting both inland nearly-freshwater desert ponds and highly saline Red Sea coastal lagoons of the Arabian Peninsula. Desert ponds and coastal lagoons, located respectively upstream and at the mouths of dry riverbeds (“wadies”), have been found to potentially become connected during periods of intense rainfall, which could allow the fish to migrate between these different habitats. Flash floods would therefore flush Arabian pupfish out to sea, requiring a rapid acclimation to a greater than 40 ppt change in salinity. To investigate the molecular pathways of salinity acclimation during such events, a Red Sea coastal lagoon and a desert pond population were sampled, with the latter exposed to a rapid increase in water salinity. Changes in branchial gene expression were investigated via genome-wide transcriptome measurements over time from 6 h to 21 days. The two natural populations displayed basal differences in genes related to ion transport, osmoregulation and immune system functions. These mechanisms were also differentially regulated in seawater transferred fish, revealing their crucial role in long-term adaptation. Other processes were only transiently activated shortly after the salinity exposure, including cellular stress response mechanisms, such as molecular chaperone synthesis and apoptosis. Tissue remodelling processes were also identified as transient, but took place later in the timeline, suggesting their importance to long-term acclimation as they likely equip the fish with lasting adaptations to their new environment. The alterations in branchial functional pathways displayed by Arabian pupfish in response to salinity increases are diverse. These reveal a large toolkit of molecular processes important for adaptation to hyperosmolarity that allow for successful colonization to a wide variety of different habitats.

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.

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.

Evolutionary responses of marine organisms to urbanized seascapes

Many of the world's major cities are located in coastal zones, resulting in urban and industrial impacts on adjacent marine ecosystems. These pressures, which include pollutants, sewage, runoff and debris, temperature increases, hardened shorelines/structures, and light and acoustic pollution, have resulted in new evolutionary landscapes for coastal marine organisms. Marine environmental changes influenced by urbanization may create new selective regimes or may influence neutral evolution via impacts on gene flow or partitioning of genetic diversity across seascapes. While some urban selective pressures, such as hardened surfaces, are similar to those experienced by terrestrial species, others, such as oxidative stress, are specific to aquatic environments. Moreover, spatial and temporal scales of evolutionary responses may differ in the ocean due to the spatial extent of selective pressures and greater capacity for dispersal/gene flow. Here, we present a conceptual framework and synthesis of current research on evolutionary responses of marine organisms to urban pressures. We review urban impacts on genetic diversity and gene flow and examine evidence that marine species are adapting, or are predicted to adapt, to urbanization over rapid evolutionary time frames. Our findings indicate that in the majority of studies, urban stressors are correlated with reduced genetic diversity. Genetic structure is often increased in urbanized settings, but artificial structures can also act as stepping stones for some hard-surface specialists, promoting range expansion. Most evidence for rapid adaptation to urban stressors comes from studies of heritable tolerance to pollutants in a relatively small number of species; however, the majority of marine ecotoxicology studies do not test directly for heritability. Finally, we highlight current gaps in our understanding of evolutionary processes in marine urban environments and present a framework for future research to address these gaps.

Evolution of Plasticity in Response to Ethanol between Sister Species with Different Ecological Histories (Drosophila melanogaster and D. simulans)

AbstractWhen populations evolve adaptive reaction norms in response to novel environments, it can occur through a process termed genetic accommodation. Under this model, the initial response to the environment is widely variable between genotypes as a result of cryptic genetic variation, which is then refined by selection to a single adaptive response. Here, I empirically test these predictions from genetic accommodation by measuring reaction norms in individual genotypes and across several time points. I compare two species of Drosophila that differ in their adaptation to ethanol (D. melanogaster and D. simulans). Both species are human commensals with a recent cosmopolitan expansion, but only D. melanogaster is adapted to ethanol exposure. Using gene expression as a phenotype and an approach that combines information about expression and alternative splicing, I find that D. simulans exhibits cryptic genetic variation in the response to ethanol, while D. melanogaster has almost no genotype-specific variation in reaction norm. This is evidence for adaptation to ethanol through genetic accommodation, suggesting that the evolution of phenotypic plasticity could be an important contributor to the ability to exploit novel resources.

Comparative transcriptomic analysis reveals mechanisms of divergence in osmotic regulation of the turbot Scophthalmus maximus

The turbot Scophthalmus maximus has evolved extensive physiological ability to adapt to multiple environmental salinities. The morphological changes of the kidney indicated the adaptability difference and similarity of turbot to salinity stress. Identify transcriptome-wide differences between low-salinity seawater (LSW, salinity 5)- and high-salinity seawater (HSW, salinity 50)-acclimated kidneys of turbot to decipher the osmotic regulation mechanism. We identified 688 differentially expressed genes (DEGs) in the LSW-acclimated kidneys and 2441 DEGs in the HSW-acclimated kidneys of turbot compared with seawater-acclimated kidneys, respectively. We investigated three patterns of gene regulation to salinity stress that involved in ion channels and transporters, functions of calcium regulation, organic osmolytes, energy demand, cell cycle regulation, and cell protection. Additionally, protein-protein interaction (PPI) analysis of DEGs suggested the presence of a frequent functional interaction pattern and that crucial genes in the PPI network are involved in hyper-osmotic regulation. Based on the analysis of comparative transcriptome data and related literature reports, we conclude that the mechanisms responsible for osmotic regulation and its divergence in turbot are related to various genes that are involved in canonical physiological functions. These findings provide insight into the divergence in osmoregulation of turbot and valuable information about osmoregulation mechanisms that will benefit other studies in this field.

Adaptation of plasticity to projected maximum temperatures and across climatically defined bioregions

Resilience to environmental stressors due to climate warming is influenced by local adaptations, including plastic responses. The recent literature has focused on genomic signatures of climatic adaptation, but little is known about how plastic capacity may be influenced by biogeographic and evolutionary processes. We investigate phenotypic plasticity as a target of climatic selection, hypothesizing that lineages that evolved in warmer climates will exhibit greater plastic adaptive resilience to upper thermal stress. This was experimentally tested by comparing transcriptomic responses within and among temperate, subtropical, and desert ecotypes of Australian rainbowfish subjected to contemporary and projected summer temperatures. Critical thermal maxima were estimated, and ecological niches delineated using bioclimatic modeling. A comparative phylogenetic expression variance and evolution model was used to assess plastic and evolved changes in gene expression. Although 82% of all expressed genes were found in the three ecotypes, they shared expression patterns in only 5 out of 236 genes that responded to the climate change experiment. A total of 532 genes showed signals of adaptive (i.e., genetic-based) plasticity due to ecotype-specific directional selection, and 23 of those responded to projected summer temperatures. Network analyses demonstrated centrality of these genes in thermal response pathways. The greatest adaptive resilience to upper thermal stress was shown by the subtropical ecotype, followed by the desert and temperate ecotypes. Our findings indicate that vulnerability to climate change will be highly influenced by biogeographic factors, emphasizing the value of integrative assessments of climatic adaptive traits for accurate estimation of population and ecosystem responses.

Using asexual vertebrates to study genome evolution and animal physiology: Banded (Fundulus diaphanus) x Common Killifish (F. heteroclitus) hybrid lineages as a model system

Wild, asexual, vertebrate hybrids have many characteristics that make them good model systems for studying how genomes evolve and epigenetic modifications influence animal physiology. In particular, the formation of asexual hybrid lineages is a form of reproductive incompatibility, but we know little about the genetic and genomic mechanisms by which this mode of reproductive isolation proceeds in animals. Asexual lineages also provide researchers with the ability to produce genetically identical individuals, enabling the study of autonomous epigenetic modifications without the confounds of genetic variation. Here, we briefly review the cellular and molecular mechanisms leading to asexual reproduction in vertebrates and the known genetic and epigenetic consequences of the loss of sex. We then specifically discuss what is known about asexual lineages of Fundulus diaphanus x F. heteroclitus to highlight gaps in our knowledge of the biology of these clones. Our preliminary studies of F. diaphanus and F. heteroclitus karyotypes from Porter's Lake (Nova Scotia, Canada) agree with data from other populations, suggesting a conserved interspecific chromosomal arrangement. In addition, genetic analyses suggest that: (a) the same major clonal lineage (Clone A) of F. diaphanus x F. heteroclitus has remained dominant over the past decade, (b) some minor clones have also persisted, (c) new clones may have recently formed, and iv) wild clones still mainly descend from F. diaphanus ♀ x F. heteroclitus ♂ crosses (96% in 2017-2018). These data suggest that clone formation may be a relatively rare, but continuous process, and there are persistent environmental or genetic factors causing a bias in cross direction. We end by describing our current research on the genomic causes and consequences of a transition to asexuality and the potential physiological consequences of epigenetic variation.

Transcriptomic Analysis of Gill and Kidney from Asian Seabass (Lates calcarifer) Acclimated to Different Salinities Reveals Pathways Involved with Euryhalinity

Asian seabass (or commonly known as barramundi), Lates calcarifer, is a bony euryhaline teleost from the Family Latidae, inhabiting nearshore, estuarine, and marine connected freshwaters throughout the tropical Indo-West Pacific region. The species is catadromous, whereby adults spawn in salinities between 28 and 34 ppt at the mouth of estuaries, with resultant juveniles usually moving into brackish and freshwater systems to mature, before returning to the sea to spawn again as adults. The species lives in both marine and freshwater habitats and can move quickly between the two; thus, the species' ability to tolerate changes in salinity makes it a good candidate for studying the salinity acclimation response in teleosts. In this study, the transcriptome of two major osmoregulatory organs (gills and kidneys) of young juvenile Asian seabass reared in freshwater and seawater were compared. The euryhaline nature of Asian seabass was found to be highly pliable and the moldability of the trait was further confirmed by histological analyses of gills and kidneys. Differences in major expression pathways were observed, with differentially expressed genes including those related to osmoregulation, tissue/organ morphogenesis, and cell volume regulation as central to the osmo-adaptive response. Additionally, genes coding for mucins were upregulated specifically under saline conditions, whereas several genes important for growth and development, as well as circadian entrainment were specifically enriched in fish reared in freshwater. Routing of the circadian rhythm mediated by salinity changes could be the initial step in salinity acclimation and possibly migration in euryhaline fish species such as the Asian seabass.

Choice consequences: salinity preferences and hatchling survival in the mangrove rivulus (Kryptolebias marmoratus)

In heterogeneous environments, mobile species should occupy habitats in which their fitness is maximized. Mangrove rivulus fish inhabit mangrove ecosystems where salinities range from 0 to 65 ppt, but are most often collected from areas with salinities of ∼25 ppt. We examined the salinity preference of mangrove rivulus in a lateral salinity gradient, in the absence of predators and competitors. Fish could swim freely for 8 h throughout the gradient with chambers containing salinities ranging from 5 to 45 ppt (or 25 ppt throughout in the control). We defined preference as the salinity in which the fish spent most of their time, and also measured preference strength, latency to begin exploring the arena, and number of transitions between chambers. To determine whether these traits were repeatable, each fish experienced three trials. Mangrove rivulus spent a greater proportion of time in salinities lower (5-15 ppt) than they occupy in the wild. Significant among-individual variation in the (multivariate) behavioral phenotype emerged when animals experienced the gradient, indicating strong potential for selection to drive behavioral evolution in areas with diverse salinity microhabitats. We also showed that mangrove rivulus had a significantly greater probability of laying eggs in low salinities compared with control or high salinities. Eggs laid in lower salinities also had higher hatching success compared with those laid in higher salinities. Thus, although mangrove rivulus can tolerate a wide range of salinities, they prefer low salinities. These results raise questions about factors that prevent mangrove rivulus from occupying lower salinities in the wild, whether higher salinities impose energetic costs, and whether fitness changes as a function of salinity.

Comparative transcriptome profiling of selected osmotic regulatory proteins in the gill during seawater acclimation of chum salmon (Oncorhynchus keta) fry

Salmonid fishes, chum salmon (Oncorhynchus keta) have the developed adaptive strategy to withstand wide salinity changes from the early life stage. This study investigated gene expression patterns of cell membrane proteins in the gill of chum salmon fry on the transcriptome level by tracking the salinity acclimation of the fish in changing environments ranging from freshwater (0 ppt) to brackish water (17.5 ppt) to seawater (35 ppt). Using GO analysis of DEGs, the known osmoregulatory genes and their functional groups such as ion transport, transmembrane transporter activity and metal ion binding were identified. The expression patterns of membrane protein genes, including pump-mediated protein (NKA, CFTR), carrier-mediated protein (NKCC, NHE3) and channel-mediated protein (AQP) were similar to those of other salmonid fishes in the smolt or adult stages. Based on the protein-protein interaction analysis between transmembrane proteins and other related genes, we identified osmotic-related genes expressed with salinity changes and analyzed their expression patterns. The findings of this study may facilitate the disentangling of the genetic basis of chum salmon and better able an understanding of the osmophysiology of the species.

Transcriptomic response to three osmotic stresses in gills of hybrid tilapia (Oreochromis mossambicus female × O. urolepis hornorum male)

Background

Osmotic stress is a widespread phenomenon in aquatic animal. The ability to cope with salinity stress and alkaline stress is quite important for the survival of aquatic species under natural conditions. Tilapia is an important commercial euryhaline fish species. What’s more tilapia is a good experimental material for osmotic stress regulation research, but the molecular regulation mechanism underlying different osmotic pressure of tilapia is still unexplored.

Results

To elucidate the osmoregulation strategy behind its hyper salinity, alkalinity and salinity-alkalinity stress of tilapia, the transcriptomes of gills in hybrid tilapia (Oreochromis mossambicus ♀ × O. urolepis hornorum ♂) under salinity stress (S: 25‰), alkalinity stress(A: 4‰) and salinity-alkalinity stress (SA: S: 15‰, A: 4‰) were sequenced using deep-sequencing platform Illumina/HiSeq-2000 and differential expression genes (DEGs) were identified. A total of 1958, 1472 and 1315 upregulated and 1824, 1940 and 1735 downregulated genes (P-value < 0.05) were identified in the salt stress, alkali stress and saline-alkali stress groups, respectively, compared with those in the control group. Furthermore, Kyoto Encyclopedia of Genes and Genomes pathway analyses were conducted in the significant different expression genes. In all significant DEGs, some of the typical genes involved in osmoregulation, including carbonic anhydrase (CA), calcium/calmodulin-dependent protein kinase (CaM kinase) II (CAMK2), aquaporin-1(AQP1), sodium bicarbonate cotransporter (SLC4A4/NBC1), chloride channel 2(CLCN2), sodium/potassium/chloride transporter (SLC12A2 / NKCC1) and other osmoregulation genes were also identified. RNA-seq results were validated with quantitative real-time PCR (qPCR), the 17 random selected genes showed a consistent direction in both RNA-Seq and qPCR analysis, demonstrated that the results of RNA-seq were reliable.

Conclusions

The present results would be helpful to elucidate the osmoregulation mechanism of aquatic animals adapting to saline-alkali challenge. This study provides a global overview of gene expression patterns and pathways that related to osmoregulation in hybrid tilapia, and could contribute to a better understanding of the molecular regulation mechanism in different osmotic stresses.

Salinity-dependent expression of ncc2 in opercular epithelium and gill of mummichog (Fundulus heteroclitus)

Mummichogs (Fundulus heteroclitus) can tolerate abrupt changes in environmental salinity because of their ability to rapidly adjust the activities of ionocytes in branchial and opercular epithelia. In turn, the concerted expression of sub-cellular effectors of ion transport underlies adaptive responses to fluctuating salinities. Exposure to seawater (SW) stimulates the expression of Na⁺/K⁺/2Cl^- cotransporter 1 (nkcc1) and cystic fibrosis transmembrane regulator (cftr) mRNAs in support of ion extrusion by SW-type ionocytes. Given the incomplete understanding of how freshwater (FW)-type ionocytes actually operate in mummichogs, the transcriptional responses essential for ion absorption in FW environments remain unresolved. In a subset of species, a 'fish-specific' Na⁺/Cl^- cotransporter denoted Ncc2 (Slc12a10) is responsible for the uptake of Na⁺ and Cl^- across the apical surface of FW-type ionocytes. In the current study, we identified an ncc2 transcript that is highly expressed in gill filaments and opercular epithelium of FW-acclimated mummichogs. Within 1 day of transfer from SW to FW, ncc2 levels in both tissues increased in parallel with reductions in nkcc1 and cftr. Conversely, mummichogs transferred from FW to SW exhibited marked reductions in ncc2 concurrent with increases in nkcc1 and cftr. Immunohistochemical analyses employing a homologous antibody revealed apical Ncc2-immunoreactivity in Na⁺/K⁺-ATPase-immunoreactive ionocytes of FW-acclimated animals. Our combined observations suggest that Ncc2/ncc2-expressing ionocytes support the capacity of mummichogs to inhabit FW environments.

Osmoregulation by the myo-inositol biosynthesis pathway in turbot Scophthalmus maximus and its regulation by anabolite and c-Myc

The induction of the myo-inositol biosynthesis (MIB) pathway in euryhaline fishes is an important component of the cellular response to osmotic challenge. The MIPS and IMPA1 genes were sequenced in turbot and found to be highly conserved in phylogenetic evolution, especially within the fish species tested. Under salinity stress in turbot, both MIPS and IMPA1 showed adaptive expression, a turning point in the level of expression occurred at 12 h in all tissues tested. We performed an RNAi assay mediated by long fragment dsRNA prepared by transcription in vitro. The findings demonstrated that knockdown of the MIB pathway weakened the function of gill osmotic regulation, and may induce a genetic compensation response in the kidney and gill to maintain physiological function. Even though the gill and kidney conducted stress reactions or compensatory responses to salinity stress, this inadequately addressed the consequences of MIB knockdown. Therefore, the survival time of turbot under salinity stress after knockdown was obviously less than that under seawater, especially under low salt stress. Pearson's correlation analysis between gene expression and dietary myo-inositol concentration indicated that the MIB pathway had a remarkable negative feedback control, and the dynamic equilibrium mediated by negative feedback on the MIB pathway played a crucial role in osmoregulation in turbot. An RNAi assay with c-Myc in vivo and the use of a c-Myc inhibitor (10058-F4) in vitro demonstrated that c-Myc was likely to positively regulate the MIB pathway in turbot.

Developing a High-Quality Linkage Map for the Atlantic Killifish Fundulus heteroclitus

Killifish (Fundulus heteroclitus) are widely distributed among different aquatic environments where they demonstrate an impressive range of highly-plastic and locally adaptive phenotypes. High-throughput sequencing has begun to unravel the mechanisms and evolutionary history of these interesting features by establishing relationships in the genotype-phenotype map. However, some genotype-phenotype analyses require a higher order of contiguity than what initial scaffolded (fragmented genome assembly where contigs have been assemble into scaffolds) genome assemblies can provide. Here, we used 5,685 high-quality RAD-Seq markers from a single mapping family to order 84% of the scaffolded genome assembly to 24 chromosomes. This serves to: 1) expand the killifish genomic toolkit, 2) estimate genome-wide recombination rates, and 3) compare genome synteny to humans and other fishes. After initially building our map, we found that the selection of thresholds for sequence data filtration highly impacted scaffold placement in the map. We outline each step of the approach that dramatically improved our map to help guide others toward more effective linkage mapping for genome assembly. Our final map supports strong conservation of genomic synteny among closely related fish species and reveals previously described chromosomal rearrangements between more distantly related clades. However, we also commonly found minor scaffold misorientations in F. heteroclitus and in other assemblies, suggesting that further mapping (such as optical mapping) is necessary for finer scale resolution of genome structure. Lastly, we discuss the problems that would be expected from misoriented/unplaced scaffolds and stress the importance of a quality mapped genome as a key feature for further investigating population and comparative genomic questions with F. heteroclitus and other taxa.

Effects of oxidative stress on sex-specific gene expression in the copepod Tigriopus californicus revealed by single individual RNA-seq

Oxidative stress reflects the imbalance of pro-oxidants and antioxidants. Prolonged oxidative stress can induce cellular damage, diseases and aging, and the effects may be sex-specific. Tigriopus californicus has recently been proposed as an alternative model system for sex-specific studies due to the absence of sex chromosomes. In this study, we used comparative transcriptomic analyses to assess sex-specific transcriptional responses to oxidative stress. Male and female individuals were maintained separately in one of three treatments: 1) control conditions with an algae diet, 2) pro-oxidant (H₂O₂) conditions with an algae diet or 3) decreased antioxidant conditions (reduced carotenoids due to a yeast diet). Single individual RNA-seq was then conducted for twenty-four libraries using Ligation Mediated RNA sequencing (LM-Seq). Variance in gene expression was partitioned into 62.3% between sexes, 26.85% among individuals and 10.85% among treatments. Within each of the three treatments, expression was biased toward females. However, compared to the control treatment, males in both pro-oxidant and decreased antioxidant treatments differentially expressed more genes while females differentially expressed fewer genes but with a greater magnitude of fold change. As the first study of copepods to apply single individual RNA-seq, the findings will contribute to a better understanding of transcriptomic variation among individuals as well as sex-specific response mechanisms to oxidative stress in the absence of sex chromosomes.

The Molecular Basis of Freshwater Adaptation in Prawns: Insights from Comparative Transcriptomics of Three Macrobrachium Species

Elucidating the molecular basis of adaptation to different environmental conditions is important because adaptive ability of a species can shape its distribution, influence speciation, and also drive a variety of evolutionary processes. For crustaceans, colonization of freshwater habitats has significantly impacted diversity, but the molecular basis of this process is poorly understood. In the current study, we examined three prawn species from the genus Macrobrachium (M. australiense, M. tolmerum, and M. novaehollandiae) to better understand the molecular basis of freshwater adaptation using a comparative transcriptomics approach. Each of these species naturally inhabit environments with different salinity levels; here, we exposed them to the same experimental salinity conditions (0‰ and 15‰), to compare expression patterns of candidate genes that previously have been shown to influence phenotypic traits associated with freshwater adaptation (e.g., genes associated with osmoregulation). Differential gene expression analysis revealed 876, 861, and 925 differentially expressed transcripts under the two salinities for M. australiense, M. tolmerum, and M. novaehollandiae, respectively. Of these, 16 were found to be unannotated novel transcripts and may be taxonomically restricted or orphan genes. Functional enrichment and molecular pathway mapping revealed 13 functionally enriched categories and 11 enriched molecular pathways that were common to the three Macrobrachium species. Pattern of selection analysis revealed 26 genes with signatures of positive selection among pairwise species comparisons. Overall, our results indicate that the same key genes and similar molecular pathways are likely to be involved with freshwater adaptation widely across this decapod group; with nonoverlapping sets of genes showing differential expression (mainly osmoregulatory genes) and signatures of positive selection (genes involved with different life history traits).

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

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

Profiling microRNA expression in Atlantic killifish (Fundulus heteroclitus) gill and responses to arsenic and hyperosmotic stress

The Atlantic killifish (Fundulus heteroclitus), native to estuarine areas of the Atlantic coast of the United States, has become a valuable ecotoxicological model as a result of its ability to acclimate to rapid environmental changes and adapt to polluted habitats. MicroRNAs (miRNAs) are highly conserved small RNAs that regulate gene expression and play critical roles in stress responses in a variety of organisms. Global miRNA expression in killifish and the potential roles miRNA have in environmental acclimation have yet to be characterized. Accordingly, we profiled miRNA expression in killifish gill for the first time and identified a small group of highly expressed, well-conserved miRNAs as well as 16 novel miRNAs not yet identified in other organisms. Killifish respond to large fluctuations in salinity with rapid changes in gene expression and protein trafficking to maintain osmotic balance, followed by a secondary phase of gene and protein expression changes that enable remodeling of the gills. Arsenic, a major environmental toxicant, was previously shown to inhibit gene expression responses in killifish gill, as well the ability of killifish to acclimate to a rapid increase in salinity. Thus, we examined the individual and combined effects of salinity and arsenic on miRNA expression in killifish gill. Using small RNA sequencing, we identified 270 miRNAs expressed in killifish, and found that miR-135b was differentially expressed in response to arsenic and at 24 h following transfer to salt water. Predicted targets of miR-135b are involved in ion transport, cell motility and migration, GTPase mediated signal transduction and organelle assembly. Consistent with previous studies of these two environmental stressors, we found a significant interaction (i.e., arsenic dependent salinity effect), whereby killifish exposed to arsenic exhibited an opposite response in miR-135b expression at 24 h post hyperosmotic challenge compared to controls. By examining mRNA expression of predicted miRNA targets during salinity acclimation and arsenic exposure, we found that miR-135b targets were significantly more likely to decrease during salinity acclimation than non-targets. Our identification of a significant interaction effect of arsenic and salinity on miR-135b expression supports the hypothesis that arsenic alters upstream regulators of stress response networks, which may adversely affect the killifish response to osmotic stress. The characterization of miRNAs in this ecotoxicological model will be a valuable resource for future studies investigating the role of miRNAs in response to environmental stress.

Insect communities in saline waters consist of realized but not fundamental niche specialists

Considering how organisms adapt to stress is essential if we are to anticipate biological responses to global change in ecosystems. Communities in stressful environments can potentially be assembled by specialists (i.e. species that only occur in a limited range of environmental conditions) and/or generalist species with wider environmental tolerances. We review the existing literature on the salinity tolerance of aquatic insects previously identified as saline specialists because they were exclusively found in saline habitats, and explore if these saline realized niche specialists are also specialists in their fundamental niches or on the contrary are fundamental niche generalist species confined to the highest salinities they can tolerate. The results suggest that species inhabiting saline waters are generalists in their fundamental niches, with a predominant pattern of high survival in freshwater-low salinity conditions, where their fitness tends to be similar or even higher than in saline waters. Additionally, their performance in freshwater tends to be similar to related strictly freshwater species, so no apparent trade-off of generalization is shown. These results are discussed in the framework of the ecological and evolutionary processes driving community assembly across the osmotic stress gradient, and their potential implications for predicting impacts from saline dilution and freshwater salinization.This article is part of the theme issue 'Salt in freshwaters: causes, ecological consequences and future prospects'.

Integrative Population and Physiological Genomics Reveals Mechanisms of Adaptation in Killifish

Adaptive divergence between marine and freshwater (FW) environments is important in generating phyletic diversity within fishes, but the genetic basis of this process remains poorly understood. Genome selection scans can identify adaptive loci, but incomplete knowledge of genotype-phenotype connections makes interpreting their significance difficult. In contrast, association mapping (genome-wide association mapping [GWAS], random forest [RF] analyses) links genotype to phenotype, but offer limited insight into the evolutionary forces shaping variation. Here, we combined GWAS, RF, and selection scans to identify loci important in adaptation to FW environments. We utilized FW-native and brackish water (BW)-native populations of Atlantic killifish (Fundulus heteroclitus) as well as a naturally admixed population between the two. We measured morphology and multiple physiological traits that differ between populations and may contribute to osmotic adaptation (salinity tolerance, hypoxia tolerance, metabolic rate, body shape) and used a reduced representation approach for genome-wide genotyping. Our results show patterns of population divergence in physiological capabilities that are consistent with local adaptation. Population genomic scans between BW-native and FW-native populations identified genomic regions evolving by natural selection, whereas association mapping revealed loci that contribute to variation for each trait. There was substantial overlap in the genomic regions putatively under selection and loci associated with phenotypic traits, particularly for salinity tolerance, suggesting that these regions and genes are important for adaptive divergence between BW and FW environments. Together, these data provide insight into the mechanisms that enable diversification of fishes across osmotic boundaries.

Reduced Swimming Performance Repeatedly Evolves on Loss of Migration in Landlocked Populations of Alewife

Whole-organism performance tasks are accomplished by the integration of morphological traits and physiological functions. Understanding how evolutionary change in morphology and physiology influences whole-organism performance will yield insight into the factors that shape its own evolution. We demonstrate that nonmigratory populations of alewife (Alosa pseudoharengus) have evolved reduced swimming performance in parallel, compared with their migratory ancestor. In contrast to theoretically and empirically based predictions, poor swimming among nonmigratory populations is unrelated to the evolution of osmoregulation and occurs despite the fact that nonmigratory alewives have a more fusiform (torpedo-like) body shape than their ancestor. Our results suggest that elimination of long-distance migration from the life cycle has shaped performance more than changes in body shape and physiological regulatory capacity.

Transcriptomic responses to extreme low salinity among locally adapted populations of Olympia oyster (Ostrea lurida)

The Olympia oyster (Ostrea lurida) is a foundation species inhabiting estuaries along the North American west coast. In California estuaries, O. lurida is adapted to local salinity regimes and populations differ in low salinity tolerance. In this study, oysters from three California populations were reared for two generations in a laboratory common garden and subsequently exposed to low salinity seawater. Comparative transcriptomics was then used to understand species-level responses to hyposmotic stress and population-level mechanisms underlying divergent salinity tolerances. Gene expression patterns indicate Olympia oysters are sensitive to hyposmotic stress: All populations respond to low salinity by up-regulating transcripts indicative of protein unfolding, DNA damage and cell cycle arrest after sub-lethal exposure. Among O. lurida populations, transcriptomic profiles differed constitutively and in response to low salinity. Despite two generations in common-garden conditions, transcripts encoding apoptosis modulators were constitutively expressed at significantly different levels in the most tolerant population. Expression of cell death regulators may facilitate cell fate decisions when salinity declines. Following low salinity exposure, oysters from the more tolerant population expressed a small number of mRNAs at significantly higher levels than less tolerant populations. Proteins encoded by these transcripts regulate ciliary activity within the mantle cavity and may function to prolong valve closure and reduce mortality in low salinity seawater. Collectively, gene expression patterns suggest sub-lethal impacts of hyposmotic stress in Olympia oysters are considerable and that even oysters with greater low salinity tolerance may be vulnerable to future freshwater flooding events.

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.

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.

Phenotypic and transcriptomic responses to salinity stress across genetically and geographically divergent Tigriopus californicus populations

Species inhabiting the North American west coast intertidal must tolerate an extremely variable environment, with large fluctuations in both temperature and salinity. Uncovering the mechanisms for this tolerance is key to understanding species' persistence. We tested for differences in salinity tolerance between populations of Tigriopus californicus copepods from locations in northern (Bodega Reserve) and southern (San Diego) California known to differ in temperature, precipitation and humidity. We also tested for differences between populations in their transcriptomic responses to salinity. Although these two populations have ~20% mtDNA sequence divergence and differ strongly in other phenotypic traits, we observed similarities in their phenotypic and transcriptomic responses to low and high salinity stress. Salinity significantly affected respiration rate (increased under low salinity and reduced under high salinity), but we found no significant effect of population on respiration or a population by salinity interaction. Under high salinity, there was no population difference in knock-down response, but northern copepods had a smaller knock-down under low salinity stress, corroborating previous results for T. californicus. Northern and southern populations had a similar transcriptomic response to salinity based on a principle components analysis, although differential gene expression under high salinity stress was three times lower in the northern population compared to the southern population. Transcripts differentially regulated under salinity stress were enriched for "amino acid transport" and "ion transport" annotation categories, supporting previous work demonstrating that the accumulation of free amino acids is important for osmotic regulation in T. californicus.

Transcriptome response of Atlantic salmon (Salmo salar) to competition with ecologically similar non-native species

Non-native species may be introduced either intentionally or unintentionally, and their impact can range from benign to highly disruptive. Non-native salmonids were introduced into Lake Ontario, Canada, to provide recreational fishing opportunities; however, the establishment of those species has been proposed as a significant barrier to the reintroduction of native Atlantic salmon (Salmo salar) due to intense interspecific competition. In this study, we compared population differences of Atlantic salmon in transcriptome response to interspecific competition. We reared Atlantic salmon from two populations (LaHave River and Sebago Lake) with fish of each of three non-native salmonids (Chinook salmon Oncorhynchus tshawytscha, rainbow trout O. mykiss, and brown trout S. trutta) in artificial streams. We used RNA-seq to assess transcriptome differences between the Atlantic salmon populations and the responses of these populations to the interspecific competition treatments after 10 months of competition in the stream tanks. We found that population differences in gene expression were generally greater than the effects of interspecific competition. Interestingly, we found that the two Atlantic salmon populations exhibited similar responses to interspecific competition based on functional gene ontologies, but the specific genes within those ontologies were different. Our transcriptome analyses suggest that the most stressful competitor (as measured by the highest number of differentially expressed genes) differs between the two study populations. Our transcriptome characterization highlights the importance of source population selection for conservation applications, as organisms with different evolutionary histories can possess different transcriptional responses to the same biotic stressors. The results also indicate that generalized predictions of the response of native species to interactions with introduced species may not be appropriate without incorporating potential population-specific response to introduced species.

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.

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.

NFAT5 genes are part of the osmotic regulatory system in Atlantic salmon (Salmo salar)

The anadromous Atlantic salmon utilizes both fresh and salt water (FW and SW) habitats during its life cycle. The parr-smolt transformation (PST) is an important developmental transition from a FW adapted juvenile parr to a SW adapted smolt. Physiological changes in osmoregulatory tissues, particularly the gill, are key in maintaining effective ion regulation during PST. Changes are initiated prior to SW exposure (preparative phase), and are completed when smolts enter the sea (activational phase) where osmotic stress may directly stimulate changes in gene expression. In this paper we identify 4 nuclear factor of activated T cells (NFAT5, an osmotic stress transcription factor) paralogues in Atlantic salmon, which showed strong homology in characterized functional domains with those identified in other vertebrates. Two of the identified paralogues (NFAT5b1 and NFAT5b2) showed increased expression following transfer from FW to SW. This effect was largest in parr that were maintained under short day photoperiod, and showed the highest increases in chloride ion levels in response to SW exposure. The results of this study suggest that NFAT5 is involved in the osmotic stress response of Atlantic salmon.

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

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