Transcriptomic imprints of adaptation to fresh water: parallel evolution of osmoregulatory gene expression in the Alewife

A key evolutionary step determining osmoregulatory ability for freshwater colonisation in early life stages of fish

Colonisation of freshwater habitats by marine animals is a remarkable evolutionary event that has enriched biodiversity in freshwater ecosystems. The acquisition of tolerance to hypotonic stress during early life stages is presumed to be essential for their successful freshwater colonisation, but very little empirical evidence has been obtained to support this idea. This study aimed to comprehend the evolutionary changes in osmoregulatory mechanisms that enhance larval freshwater tolerance in amphidromous fishes, which typically spend their larval period in marine (ancestral) habitats and the rest of their life history stages in freshwater (derived) habitats. We compared the life history patterns and changes in larval survivorship and gene expression depending on salinity among three congeneric marine-originated amphidromous goby species (Gymnogobius), which had been suggested to differ in their larval dependence on freshwater habitats. An otolith microchemical analysis and laboratory-rearing experiment confirmed the presence of freshwater residents only in G. urotaenia and higher larval survivorship of this species in the freshwater condition than in the obligate amphidromous G. petschiliensis and G. opperiens. Larval whole-body transcriptome analysis revealed that G. urotaenia from both amphidromous and freshwater-resident populations exhibited the greatest differences in expression levels of several osmoregulatory genes, including aqp3, which is critical for water discharge from their body during early fish development. The present results consistently support the importance of enhanced freshwater tolerance and osmoregulatory plasticity in larval fish to establish freshwater forms, and further identified key candidate genes for larval freshwater adaptation and colonisation in the goby group.

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.

Salinity-induced transcriptome profiles in marine and freshwater threespine stickleback after an abrupt 6-hour exposure

Saltwater and freshwater environments have opposing physiological challenges, yet, there are fish species that are able to enter both habitats during short time spans, and as individuals they must therefore adjust quickly to osmoregulatory contrasts. In this study, we conducted an experiment to test for plastic responses to abrupt salinity changes in two populations of threespine stickleback, Gasterosteus aculeatus, representing two ecotypes (freshwater and ancestral saltwater). We exposed both ecotypes to abrupt native (control treatment) and non‐native salinities (0‰ and 30‰) and sampled gill tissue for transcriptomic analyses after 6 h of exposure. To investigate genomic responses to salinity, we analyzed four different comparisons; one for each ecotype (in their control and exposure salinity; (1) and (2), one between ecotypes in their control salinity (3), and the fourth comparison included all transcripts identified in (3) that did not show any expressional changes within ecotype in either the control or the exposed salinity (4)). Abrupt salinity transfer affected the expression of 10 and 1530 transcripts for the saltwater and freshwater ecotype, respectively, and 1314 were differentially expressed between the controls, including 502 that were not affected by salinity within ecotype (fixed expression). In total, these results indicate that factors other than genomic expressional plasticity are important for osmoregulation in stickleback, due to the need for opposite physiological pathways to survive the abrupt change in salinity.

Repeated Genetic Targets of Natural Selection Underlying Adaptation of Fishes to Changing Salinity

Ecological transitions across salinity boundaries have led to some of the most important diversification events in the animal kingdom, especially among fishes. Adaptations accompanying such transitions include changes in morphology, diet, whole-organism performance, and osmoregulatory function, which may be particularly prominent since divergent salinity regimes make opposing demands on systems that maintain ion and water balance. Research in the last decade has focused on the genetic targets underlying such adaptations, most notably by comparing populations of species that are distributed across salinity boundaries. Here, we synthesize research on the targets of natural selection using whole-genome approaches, with a particular emphasis on the osmoregulatory system. Given the complex, integrated and polygenic nature of this system, we expected that signatures of natural selection would span numerous genes across functional levels of osmoregulation, especially salinity sensing, hormonal control, and cellular ion exchange mechanisms. We find support for this prediction: genes coding for V-type, Ca2+, and Na+/K+-ATPases, which are key cellular ion exchange enzymes, are especially common targets of selection in species from six orders of fishes. This indicates that while polygenic selection contributes to adaptation across salinity boundaries, changes in ATPase enzymes may be of particular importance in supporting such transitions.

Physiological, morphological and transcriptomic responses of Tibetan naked carps (Gymnocypris przewalskii) to salinity variations

Gymnocypris przewalskii is a native cyprinid fish that dwells in the Lake Qinghai with salinity of 12-13‰. It migrates annually to the freshwater rivers for spawning, experiencing the significant changes in salinity. In the present study, we performed the physiological, morphological and transcriptomic analyses to understand the osmoregulation in G. przewalskii. The physiological assay showed that the osmotic pressure of G. przewalskii was almost isosmotic to the brackish lake water. The low salinity reduced its ionic concentrations and osmotic pressure. The plasticity of gill microstructure was linked to the salinity variations, including the presence of mucus and intact tight junctions in brackish water and the development of the mitochondria-rich cells and the loosened tight junctions in freshwater. RNA-seq analysis identified 1926 differentially expressed genes, including 710 and 1216 down- and up-regulated genes in freshwater, which were enriched in ion transport, cell-cell adhesion, and mucus secretion. Genes in ion uptake were activated in low salinity, and mucus pathways and tight junction showed the higher transcription in brackish water. The isosmoticity between the body fluid and the environment suggested G. przewalskii was in the metabolic-saving condition in the brackish water. The decreased salinity disrupted this balance, which activated the ion uptake in freshwater to maintain osmotic homeostasis. The gill remodeling was involved in this process through the development of the mitochondria-rich cells to enhance ion uptake. The current finding provided insights into the potential mechanisms of G. przewalskii to cope with salinity alteration.

Recognizing salinity threats in the climate crisis

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

The genetics of adaptation in freshwater Eurasian shad (Alosa)

Studying the genetics of phenotypic convergence can yield important insights into adaptive evolution. Here, we conducted a comparative genomic study of four lineages (species and subspecies) of anadromous shad (Alosa) that have independently evolved life cycles entirely completed in freshwater. Three naturally diverged (A. fallax lacustris, A. f. killarnensis, and A. macedonica), and the fourth (A. alosa) was artificially landlocked during the last century. To conduct this analysis, we assembled and annotated a draft of the A. alosa genome and generated whole‐genome sequencing for 16 anadromous and freshwater populations of shad. Widespread evidence for parallel genetic changes in freshwater populations within lineages was found. In freshwater A. alosa, which have only been diverging for tens of generations, this shows that parallel adaptive evolution can rapidly occur. However, parallel genetic changes across lineages were comparatively rare. The degree of genetic parallelism was not strongly related to the number of shared polymorphisms between lineages, thus suggesting that other factors such as divergence among ancestral populations or environmental variation may influence genetic parallelism across these lineages. These overall patterns were exemplified by genetic differentiation involving a paralog of ATPase‐α1 that appears to be under selection in just two of the more distantly related lineages studied, A. f. lacustris and A. alosa. Our findings provide insights into the genetic architecture of adaptation and parallel evolution along a continuum of population divergence.

Population genomics provides insights into the population structure and temperature-driven adaptation of Collichthys lucidus

BACKGROUND

Understanding the genetic structure and local adaptive evolutionary mechanisms of marine organisms is crucial for the management of biological resources. As the ecologically and commercially important small-sized shallow-sea fish, Collichthys lucidus plays a vital role in the structure and functioning of marine ecosystem processes. C. lucidus has been shown to have an obvious population structure. Therefore, it is an ideal candidate for investigating population differentiation and local adaptation under heterogeneous environmental pressure.

RESULTS

A total of 184,708 high-quality single nucleotide polymorphisms (SNPs) were identified and applied to elucidate the fine-scale genetic structure and local thermal adaptation of 8 C. lucidus populations. Population structure analysis based on all SNPs indicated that the northern group and southern group of C. lucidus have a strong differentiation. Moreover, 314 SNPs were found to be significantly associated with temperature variation, and annotations of genes containing temperature-related SNPs suggested that genes were involved in material (protein, lipid, and carbohydrate) metabolism and immune responses.

CONCLUSION

The high genetic differentiation of 8 C. lucidus populations may have been caused by long-term geographic isolation during the glacial period. Moreover, we suspected that variation in these genes associated with material (protein, lipid, and carbohydrate) metabolism and immune responses was critical for adaptation to spatially heterogeneous temperatures in natural C. lucidus populations. In conclusion, this study could help us determine how C. lucidus populations will respond to future ocean temperature rising.

Effects of Bicarbonate Stress on Serum Ions and Gill Transporters in Alkali and Freshwater Forms of Amur Ide (Leuciscus waleckii)

The Amur ide (Leuciscus waleckii) is a fish in the Cyprinidae family. Compared with other Amur ide living in freshwater ecosystems, the Amur ide population in Lake Dali Nor of China is famous for its high tolerance to the alkaline conditions of 54 mM (pH 9.6). Yet, surprisingly, the ionoregulatory mechanism responsible for this remarkable alkaline adaptation remains unclear. Therefore, this study sought to investigate how bicarbonate affects the acid-base balancing and ionoregulatory responses of this animal. Here, using a comparative approach, the alkali form of Amur ide and its ancestral freshwater form living in other freshwater basins were each exposed to 50 mM (pH 9.59 ± 0.09), a level close to the alkalinity of Lake Dali Nor, and their physiological (AE1) adjustment of ions and acid-base regulation were investigated. This study highlighted differences in blood pH and serum ions (e.g., Na⁺, K⁺, Cl⁻, and Ca²⁺), Na⁺/K⁺ ATPase (NKA) activity and its mRNA level, and mRNA expression of gill transporters (Na⁺/H⁺ exchanger member 2 and/or 3, Na⁺/HCO3- cotransporter (NBC1), Cl⁻/HCO3- exchanger, Na⁺/Cl⁻ cotransporter (NCC), Na⁺/K⁺/2Cl⁻ (NKCC1), SLC26A5, and SLC26A6) for alkalinity adaptation between the two forms of Amur ide differing in alkalinity tolerance. Specifically, close relationships among the serum Na⁺ and mRNA levels of NCC, NKCC1, and NHE, and also NKA and NBC1, in addition to serum Cl⁻ and bicarbonate transporters (e.g., SLC26A5 and SLC26A6), characterized the alkali form of Amur ide. We propose that this ecotype can ensure its transepithelial Cl⁻ and Na⁺ uptake/base secretions are highly functional, by its basolateral NKA with NBC1 and apical ionic transporters, and especially NCC incorporated with other transporters (e.g., SLC26). This suggests an evolved strong ability to maintain an ion osmotic and acid-base balance for more effectively facilitating its adaptability to the high alkaline environment. This study provides new insights into the physiological responses of the alkaline form of the Amur ide fish for adapting to extreme alkaline conditions. This information could be used as a reference to cultivating alkaline-tolerant fish species in abandoned alkaline waters.

Physiological plasticity and acclimatory responses to salinity stress are ion-specific in the mayfly, Neocloeon triangulifer

Freshwater salinization is a rapidly emerging ecological issue and is correlated with significant declines in aquatic biodiversity. It remains unclear how changing salinity regimes affect the physiology of sensitive aquatic insects. We used the parthenogenetic mayfly, Neocloeon triangulifer, to ask how ionic exposure history alters physiological processes and responses to subsequent major ion exposures. Using radiotracers (²²Na, ³⁵SO₄, and ⁴⁵Ca), we observed that mayflies chronically reared in elevated sodium or sulfate (157 mg L^-1 Na or 667 mg L^-1 SO₄) had 2-fold (p < 0.0001) and 8-fold (p < 0.0001) lower ion uptake rates than mayflies reared in dilute control water (16 mg L^-1 Na and 23 mg L^-1 SO₄) and subsequently transferred to elevated salinities, respectively. These acclimatory ion transport changes provided protection in 96-h toxicity bioassays for sodium, but not sulfate. Interestingly, calcium uptake was uniformly much lower and minimally influenced by exposure history, but was poorly tolerated in the toxicity bioassays. With qRT-PCR, we observed that the expression of many ion transporter genes in mayflies was influenced by elevated salinity in an ion-specific manner (general upregulation in response to sulfate, downregulation in response to calcium). Elevated sodium exposure had minimal influence on the same genes. Finally, we provide novel light microscopic evidence of histomorphological changes within the epithelium of the Malpighian tubules (insect primary excretory system) that undergoes cellular degeneration and necrosis secondary to calcium toxicity. We conclude that physiological plasticity to salinity stress is ion-specific and provide evidence for ion-specific toxicity mechanisms in N. triangulifer.

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.

Comprehensive phylogeny of Konosirus punctatus (Clupeiformes: Clupeidae) based on transcriptomic data

Konosirus punctatus is an important species for the structure of marine ecosystems. Meanwhile, it is a native species in the northwest Pacific Ocean and supports important commercial fishery. In the present study, we generated the whole transcriptome of K. punctatus from combined tissues (muscle, liver, gill, heart, kidney, swim bladder and sexual gland) using Illumina RNA-seq technology and a total of 46087110 clean reads were obtained, corresponding to 6531521430 nucleotides. Meanwhile, 10000 clean reads were randomly selected and compared with NT database to examine the possible contamination. Results showed that 6754 clean reads were distributed among some species closely related with K. punctatus, indicating no-pollution. De novo assembly was performed and all clean reads were assembled to produce 71610 longest unigenes with an N50 of 906 bp. Among all the unigenes, 43974 unigenes were annotated in at least one database and 3172 unigenes were annotated in all databases. All unigenes were further analyzed to predict the gene structure and we have obtained a total of 54864 coding sequences (CDSs) and 17326 simple sequence repeats (SSRs). Saturation analyses were applied to evaluate the accuracy of gene expression and we hypothesized that the detection of gene expression might be effective. Finally, single-copy orthologous genes were applied to construct the phylogenetic relationship of K. punctatus. Results showed that K. punctatus diverged from the common ancestor with Alosa alosa, Alosa pseudoharengus and Sardina pilchardus at approx. 61.16–92.52 MYA. The present study will provide a foundational molecular information for the biological research of K. punctatus.

Using seasonal genomic changes to understand historical adaptation to new environments: Parallel selection on stickleback in highly-variable estuaries

Parallel evolution is considered strong evidence for natural selection. However, few studies have investigated the process of parallel selection as it plays out in real time. The common approach is to study historical signatures of selection in populations already well adapted to different environments. Here, to document selection under natural conditions, we study six populations of threespine stickleback (Gasterosteus aculeatus) inhabiting bar-built estuaries that undergo seasonal cycles of environmental changes. Estuaries are periodically isolated from the ocean due to sandbar formation during dry summer months, with concurrent environmental shifts that resemble the long-term changes associated with postglacial colonization of freshwater habitats by marine populations. We used pooled whole-genome sequencing to track seasonal allele frequency changes in six of these populations and search for signatures of natural selection. We found consistent changes in allele frequency across estuaries, suggesting a potential role for parallel selection. Functional enrichment among candidate genes included transmembrane ion transport and calcium binding, which are important for osmoregulation and ion balance. The genomic changes that occur in threespine stickleback from bar-built estuaries could provide a glimpse into the early stages of adaptation that have occurred in many historical marine to freshwater transitions.

Parallelism in gene expression between foothill and alpine ecotypes in Arabidopsis arenosa

Parallel adaptation results from the independent evolution of similar traits between closely related lineages and allows us to test to what extent evolution is repeatable. Similar gene expression changes are often detected but the identity of genes shaped by parallel selection and the causes of expression parallelism remain largely unknown. By comparing genomes and transcriptomes of four distinct foothill-alpine population pairs across four treatments, we addressed the genetic underpinnings, plasticity and functional consequences of gene expression parallelism in alpine adaptation. Seeds of eight populations of Arabidopsis arenosa were raised under four treatments that differed in temperature and irradiance, factors varying strongly with elevation. Parallelism in differential gene expression between the foothill and alpine ecotypes was quantified by RNA-seq in leaves of young plants. By manipulating temperature and irradiance, we also tested for parallelism in plasticity (i.e., gene-environment interaction, GEI). In spite of global non-parallel patterns transcriptome wide, we found significant parallelism in gene expression at the level of individual loci with an over-representation of genes involved in biotic stress response. In addition, we demonstrated significant parallelism in GEI, indicating a shared differential response of the originally foothill versus alpine populations to environmental variation across mountain regions. A fraction of genes showing expression parallelism also encompassed parallel outliers for genomic differentiation, with greater enrichment of such variants in cis-regulatory elements in some mountain regions. In summary, our results suggest frequent evolutionary repeatability in gene expression changes associated with the colonization of a challenging environment that combines constitutive expression differences and plastic interaction with the surrounding environment.

Investigating Diadromy in Fishes and Its Loss in an -Omics Era

Diadromy, the predictable movements of individuals between marine and freshwater environments, is biogeographically and phylogenetically widespread across fishes. Thus, despite the high energetic and potential fitness costs involved in moving between distinct environments, diadromy appears to be an effective life history strategy. Yet, the origin and molecular mechanisms that underpin this migratory behavior are not fully understood. In this review, we aim first to summarize what is known about diadromy in fishes; this includes the phylogenetic relationship among diadromous species, a description of the main hypotheses regarding its origin, and a discussion of the presence of non-migratory populations within diadromous species. Second, we discuss how recent research based on -omics approaches (chiefly genomics, transcriptomics, and epigenomics) is beginning to provide answers to questions on the genetic bases and origin(s) of diadromy. Finally, we suggest future directions for -omics research that can help tackle questions on the evolution of diadromy.

Capacity for freshwater acclimation and differences in the transcription of ion transporter genes underlying different migratory life histories of Takifugu fish

The genus Takifugu is a group of approximately 20 species of puffer fishes living in a wide range of salinity environments around East Asian countries. This group presents a broad spectrum of evolutionary stages adapted to anadromy as a result of speciation that occurred a short time (2-5 million years) ago on an evolutionary timescale. This group thus can be considered as a model for studying the evolutionary mechanisms of anadromy. We firstly conducted a transfer experiment from seawater to low-salinity waters on five Takifugu species: two anadromous species T. obscurus and T. ocellatus, two euryhaline wanderer marine species T. rubripes and T. niphobles, and a strictly marine species T. snyderi, and confirmed that the capacity for acclimation to hypotonic environments was associated with their life history strategies. Next, transcriptomes of the gill and intestine of these species in hypotonic condition were compared to those under hypertonic condition for each species using RNA-Sequencing so as to determine possible candidate transporters playing an important role on freshwater adaptation. As this analysis suggested that cftr, encoding an important ion transporter for seawater acclimation in the gill, and ncc, encoding a transporter that is suggested to play important osmoregulatory roles in the intestine, are important candidates, their expression was validated by quantitative real-time PCR analysis. Expression of cftr was downregulated in the gills of the four euryhaline species under the hypotonic condition, but no change was detected in the gill of stenohaline T. snyderi, which may be one reason for the poor hypotonic acclimation capacity of T. snyderi. Expression of ncc was clearly upregulated in the intestines of the two anadromous species under the hypotonic condition, but not in other three species. Different ion transporter expression patterns between the five species indicate that the transcriptional regulation of cftr in the gill and ncc in the intestine may be important for the improvement of hypotonic acclimation capacity and evolution of anadromy in the Takifugu species.

Echolocation call frequency variation in horseshoe bats: molecular basis revealed by comparative transcriptomics

Recently diverged taxa with contrasting phenotypes offer opportunities for unravelling the genetic basis of phenotypic variation in nature. Horseshoe bats are a speciose group that exhibit a derived form of high-duty cycle echolocation in which the inner ear is finely tuned to echoes of the narrowband call frequency. Here, by focusing on three recently diverged subspecies of the intermediate horseshoe bat (Rhinolophus affinis) that display divergent echolocation call frequencies, we aim to identify candidate loci putatively involved in hearing frequency variation. We used de novo transcriptome sequencing of two mainland taxa (himalayanus and macrurus) and one island taxon (hainanus) to compare expression profiles of thousands of genes. By comparing taxa with divergent call frequencies (around 15 kHz difference), we identified 252 differentially expressed genes, of which six have been shown to be involved in hearing or deafness in human/mouse. To obtain further validation of these results, we applied quantitative reverse transcription-PCR to the candidate gene FBXL15 and found a broad association between the level of expression and call frequency across taxa. The genes identified here represent strong candidate loci associated with hearing frequency variation in bats.

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.

Error, noise and bias in de novo transcriptome assemblies

De novo transcriptome assembly is a powerful tool, and has been widely used over the last decade for making evolutionary inferences. However, it relies on two implicit assumptions: that the assembled transcriptome is an unbiased representation of the underlying expressed transcriptome, and that expression estimates from the assembly are good, if noisy approximations of the relative abundance of expressed transcripts. Using publicly available data for model organisms, we demonstrate that, across assembly algorithms and data sets, these assumptions are consistently violated. Bias exists at the nucleotide level, with genotyping error rates ranging from 30% to 83%. As a result, diversity is underestimated in transcriptome assemblies, with consistent underestimation of heterozygosity in all but the most inbred samples. Even at the gene level, expression estimates show wide deviations from map-to-reference estimates, and positive bias at lower expression levels. Standard filtering of transcriptome assemblies improves the robustness of gene expression estimates but leads to the loss of a meaningful number of protein-coding genes, including many that are highly expressed. We demonstrate a computational method, length-rescaled CPM, to partly alleviate noise and bias in expression estimates. Researchers should consider ways to minimize the impact of bias in transcriptome assemblies.

Restoration-mediated secondary contact leads to introgression of alewife ecotypes separated by a colonial-era dam

Secondary contact may have important implications for ecological and evolutionary processes; however, few studies have tracked the outcomes of secondary contact from its onset in natural ecosystems. We evaluated an anadromous alewife (Alosa pseudoharengus ) reintroduction project in Rogers Lake (Connecticut, USA), which contains a landlocked alewife population that was isolated as a result of colonial-era damming. After access to the ocean was restored, adult anadromous alewife were stocked into the lake. We assessed anadromous juvenile production, the magnitude and direction of introgression, and the potential for competition between ecotypes. We obtained fin clips from all adult alewife stocked into the lake during the restoration and a sample of juveniles produced in the lake two years after the stocking began. We assessed the ancestry of juveniles using categorical assignment and pedigree reconstruction with newly developed microhaplotype genetic markers. Anadromous alewives successfully spawned in the lake and hybridized with the landlocked population. Parentage assignments revealed that male and female anadromous fish contributed equally to juvenile F1 hybrids. The presence of landlocked backcrosses shows that some hybrids were produced within the first two years of secondary contact, matured in the lake, and reproduced. Therefore, introgression appears directional, from anadromous into landlocked, in the lake environment. Differences in estimated abundance of juveniles of different ecotypes in different habitats were also detected, which may reduce competition between ecotypes as the restoration continues. Our results illustrate the utility of restoration projects to study the outcomes of secondary contact in real ecosystems.

Parallel selection on ecologically relevant gene functions in the transcriptomes of highly diversifying salmonids

Background

Salmonid fishes are characterised by a very high level of variation in trophic, ecological, physiological, and life history adaptations. Some salmonid taxa show exceptional potential for fast, within-lake diversification into morphologically and ecologically distinct variants, often in parallel; these are the lake-resident charr and whitefish (several species in the genera Salvelinus and Coregonus). To identify selection on genes and gene categories associated with such predictable diversifications, we analysed 2702 orthogroups (4.82 Mbp total; average 4.77 genes/orthogroup; average 1783 bp/orthogroup). We did so in two charr and two whitefish species and compared to five other salmonid lineages, which do not evolve in such ecologically predictable ways, and one non-salmonid outgroup.

Results

All selection analyses are based on Coregonus and Salvelinus compared to non-diversifying taxa. We found more orthogroups were affected by relaxed selection than intensified selection. Of those, 122 were under significant relaxed selection, with trends of an overrepresentation of serine family amino acid metabolism and transcriptional regulation, and significant enrichment of behaviour-associated gene functions. Seventy-eight orthogroups were under significant intensified selection and were enriched for signalling process and transcriptional regulation gene ontology terms and actin filament and lipid metabolism gene sets. Ninety-two orthogroups were under diversifying/positive selection. These were enriched for signal transduction, transmembrane transport, and pyruvate metabolism gene ontology terms and often contained genes involved in transcriptional regulation and development. Several orthogroups showed signs of multiple types of selection. For example, orthogroups under relaxed and diversifying selection contained genes such as ap1m2, involved in immunity and development, and slc6a8, playing an important role in muscle and brain creatine uptake. Orthogroups under intensified and diversifying selection were also found, such as genes syn3, with a role in neural processes, and ctsk, involved in bone remodelling.

Conclusions

Our approach pinpointed relevant genomic targets by distinguishing among different kinds of selection. We found that relaxed, intensified, and diversifying selection affect orthogroups and gene functions of ecological relevance in salmonids. Because they were found consistently and robustly across charr and whitefish and not other salmonid lineages, we propose these genes have a potential role in the replicated ecological diversifications.

Effect of salinity fluctuation on the transcriptome of the Japanese mantis shrimp Oratosquilla oratoria

Salinity fluctuation may detrimentally affect the composition and biological processes of crustaceans. As a euryhaline crustacean, Oratosquilla oratoria can survive at salinities ranging from 20 psu to 40 psu. Therefore, we designed five salinity gradients (20, 25, 30, 35, and 40 psu) and 66.39 Gb clean transcriptome data were obtained after O. oratorias were exposed to each gradient for 24 h. All clean data were spliced into 50,482 unigenes and 17,035 unigenes were annotated in at least one database. Compared with 30 psu, 1010, 851, 1733 and 2188 differentially expressed genes were obtained at 20, 25, 35 and 40 psu, respectively. Results also showed that the osmoregulation of O. oratoria is primarily regulated by lipid and amino acid metabolism, amongst others. No significant up-regulated pathways were enriched at 25 psu and 35 psu, although more significant down-regulated pathways were obtained at 35 psu. Therefore, we assumed that the optimum survival salinity of O. oratoria may range from 25 psu to 35 psu. However, 35 psu may be more suitable for O. oratoria. In addition, 55 unigenes that encode putative inorganic ion exchanges were identified. This study aims to provide fundamental information for understanding the osmoregulation mechanisms of crustaceans.

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).

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.

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.

Predominance of cis-regulatory changes in parallel expression divergence of sticklebacks

Regulation of gene expression is thought to play a major role in adaptation, but the relative importance of cis- and trans- regulatory mechanisms in the early stages of adaptive divergence is unclear. Using RNAseq of threespine stickleback fish gill tissue from four independent marine-freshwater ecotype pairs and their F1 hybrids, we show that cis-acting (allele-specific) regulation consistently predominates gene expression divergence. Genes showing parallel marine-freshwater expression divergence are found near to adaptive genomic regions, show signatures of natural selection around their transcription start sites and are enriched for cis-regulatory control. For genes with parallel increased expression among freshwater fish, the quantitative degree of cis- and trans-regulation is also highly correlated across populations, suggesting a shared genetic basis. Compared to other forms of regulation, cis-regulation tends to show greater additivity and stability across different genetic and environmental contexts, making it a fertile substrate for the early stages of adaptive evolution.

The evolutionary consequences for seawater performance and its hormonal control when anadromous Atlantic salmon become landlocked

Populations of anadromous fish have become landlocked in relatively recent geological history (<14,000 years), but the evolutionary impacts of this altered lifecycle on traits underlying seawater performance have not been established. In order to examine the effects of relaxed selection on seawater traits, anadromous and landlocked Atlantic salmon were reared under identical conditions and examined for differences in seawater performance and its underlying physiological and endocrine control during the time of spring downstream migration. Salinity tolerance, survival and initial growth in seawater were greater in anadromous than in landlocked salmon. Abundance of the seawater isoform of gill Na⁺/K⁺-ATPase increased in spring in both strains but was greater in anadromous salmon. Hormones associated with seawater acclimation (adrenocorticotropic hormone, cortisol and growth hormone) increased in spring in both strains but were higher in anadromous salmon, whereas plasma thyroid hormones did not differ. Hypothalamic urotensin I mRNA levels also increased in spring and were higher in the anadromous strain. The results provide evidence that salinity tolerance and associated physiological traits are regulated by seasonal stimulation of the hypothalamic-pituitary-interrenal axis, and that relaxed selection on seawater entry traits has decreased this stimulation in landlocked salmon.

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.

(Non)Parallel Evolution

Parallel evolution across replicate populations has provided evolutionary biologists with iconic examples of adaptation. When multiple populations colonize seemingly similar habitats, they may evolve similar genes, traits, or functions. Yet, replicated evolution in nature or in the laboratory often yields inconsistent outcomes: Some replicate populations evolve along highly similar trajectories, whereas other replicate populations evolve to different extents or in distinct directions. To understand these heterogeneous outcomes, biologists are increasingly treating parallel evolution not as a binary phenomenon but rather as a quantitative continuum ranging from parallel to nonparallel. By measuring replicate populations’ positions along this (non)parallel continuum, we can test hypotheses about evolutionary and ecological factors that influence the extent of repeatable evolution. We review evidence regarding the manifestation of (non)parallel evolution in the laboratory, in natural populations, and in applied contexts such as cancer. We enumerate the many genetic, ecological, and evolutionary processes that contribute to variation in the extent of parallel evolution.

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.

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.

Transcriptomic responses to environmental change in fishes: Insights from RNA sequencing

The need to better understand how plasticity and evolution affect organismal responses to environmental variability is paramount in the face of global climate change. The potential for using RNA sequencing (RNA-seq) to study complex responses by non-model organisms to the environment is evident in a rapidly growing body of literature. This is particularly true of fishes for which research has been motivated by their ecological importance, socioeconomic value, and increased use as model species for medical and genetic research. Here, we review studies that have used RNA-seq to study transcriptomic responses to continuous abiotic variables to which fishes have likely evolved a response and that are predicted to be affected by climate change (e.g., salinity, temperature, dissolved oxygen concentration, and pH). Field and laboratory experiments demonstrate the potential for individuals to respond plastically to short- and long-term environmental stress and reveal molecular mechanisms underlying developmental and transgenerational plasticity, as well as adaptation to different environmental regimes. We discuss experimental, analytical, and conceptual issues that have arisen from this work and suggest avenues for future study.