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

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.

Exposure to salinity induces oxidative damage and changes in the expression of genes related to appetite regulation in Nile tilapia (Oreochromis niloticus)

Variations in water salinity and other extrinsic factors have been shown to induce changes in feeding rhythms and growth in fish. However, it is unknown whether appetite-related hormones mediate these changes in Nile tilapia (Oreochromis niloticus), an important species for aquaculture in several countries. This study aimed to evaluate the expression of genes responsible for appetite regulation and genes related to metabolic and physiological changes in tilapia exposed to different salinities. Moreover, the study proposed to sequence and to characterize the cart, cck, and pyy genes, and to quantify their expression in the brain and intestine of the fish by quantitative polymerase chain reaction (qPCR). The animals were exposed to three salinities: 0, 6, and 12 parts per thousand (ppt) of salt for 21 days. Furthermore, lipid peroxidation, reactive oxygen species, DNA damage, and membrane fluidity in blood cells were quantified by flow cytometry. The results indicated an increased expression of cart, pyy, and cck and a decreased expression of npy in the brain, and the same with cck and npy in the intestine of fish treated with 12 ppt. This modulation and other adaptive responses may have contributed to the decrease in weight gain, specific growth rate, and final weight. In addition, we showed oxidative damage in blood cells resulting from increasing salinity. These results provide essential data on O. niloticus when exposed to high salinities that have never been described before and generate knowledge necessary for developing biotechnologies that may help improve the production of economically important farmed fish.

Regulating Strategies of Transcription and Alternative Splicing for Cold Tolerance Harpadon nehereus Fish

In recent years, Harpadon nehereus gradually become a dominant species with great potential for exploitation in the East China Sea, and it is worth investigating whether H. nehereus would tolerate cold stress to continue to expand into the colder northern waters. The molecular regulation level is favorable evidence to explore the cold tolerance of H. nehereus, a total of 6,650, 1,936, and 2,772 differentially expressed genes (DEGs) in transcription regulation, and 4,409, 1,250, and 2,303 differential alternative splicing genes (DASGs) in alternative splicing regulation were identified in H. nehereus at 13, 15, and 17°C, respectively, importantly, 47 genes were identified as the key candidate genes for cold tolerance in H. nehereus. In transcription regulation, up-regulated DEGs were enriched in metabolic process terms and ribosome, spliceosome pathway, etc., while down-regulated DEGs were enriched in signal transduction terms, focal adhesion, proteoglycans in cancer pathway, etc., at 13, 15, and 17°C, respectively. In alternative splicing regulation, spliceosome, mRNA surveillance pathway, etc., were significantly enriched in DASGs. In a word, H. nehereus adapts to cold environments mainly through transcription and translation, transmembrane transport, protein modification, etc., while cold stress may also induce some diseases in H. nehereus.

Different transcriptomic architecture of the gill epithelia in Nile and Mozambique tilapia after salinity challenge

Tilapiine fishes of the genus Oreochromis vary in their euryhaline capabilities, therefore inhabiting aquatic environments of different salinities across the African continent. We analyzed the differential gene expression in the gills before and after 6 weeks salinity challenge between the highly tolerant Mozambique tilapia (Oreochromis mossambicus) and the less tolerant Nile tilapia (O. niloticus). The pathways triggered by salinity in both tilapia species reveal immune and cell stress responses as well as turnover of ionocytes. Nevertheless, the actual differential expressed genes vary between these two species, pointing at differential transcriptomic architecture, which likely contribute to the species osmoregulation capabilities in elevated salinities.

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.

Chemical niches and ionoregulatory traits: applying ionoregulatory physiology to the conservation management of freshwater fishes

Alterations in water chemistry can challenge resident fish species. More specifically, chemical changes that disrupt ion balance will negatively affect fish health and impact physiological and ecological performance. However, our understanding of which species and populations are at risk from ionoregulatory disturbances in response to changing freshwater environments is currently unclear. Therefore, we propose a novel framework for incorporating ionoregulatory physiology into conservation management of inland fishes. This framework introduces the concepts of fundamental chemical niche, which is the tolerable range of chemical conditions for a given species based on laboratory experiments, and realized chemical niche, which is the range of chemical conditions in which a species resides based on distribution surveys. By comparing these two niches, populations that may be at risk from ionoregulatory disturbances and thus require additional conservation considerations can be identified. We highlight the potential for commonly measured ionoregulatory traits to predict fundamental and realized chemical niches but caution that some traits may not serve as accurate predictors despite being important for understanding ionoregulatory mechanisms. As a sample application of our framework, the minimum pH distribution (realized niche) and survival limit pH (fundamental niche) of several North American fishes were determined by systematic review and were compared. We demonstrate that ionoregulatory capacity is significantly correlated with a realized niche for many species, highlighting the influence of ionoregulatory physiology on fish distribution patterns along chemical gradients. Our aim is that this framework will stimulate further research in this field and result in a broader integration of physiological data into conservation management decisions for inland waters.

Increased polyamine levels and maintenance of γ-aminobutyric acid (Gaba) homeostasis in the gills is indicative of osmotic plasticity in killifish

The Fundulus genus of killifish includes species that inhabit marshes along the U.S. Atlantic coast and the Gulf of Mexico, but differ in their ability to adjust rapidly to fluctuations in salinity. Previous work suggests that euryhaline killifish stimulate polyamine biosynthesis and accumulate putrescine in the gills during acute hypoosmotic challenge. Despite evidence that polyamines have an osmoregulatory role in euryhaline killifish species, their function in marine species is unknown. Furthermore, the consequences of hypoosmotic-induced changes in polyamine synthesis on downstream pathways, such as ƴ-aminobutyric acid (Gaba) production, have yet to be explored. Here, we examined the effects of acute hypoosmotic exposure on polyamine, glutamate, and Gaba levels in the gills of a marine (F. majalis) and two euryhaline killifish species (F. heteroclitus and F. grandis). Fish acclimated to 32 ppt or 12 ppt water were transferred to fresh water, and concentrations of glutamate (Glu), Gaba, and the polyamines putrescine (Put), spermidine (Spd), and spermine (Spm) were measured in the gills using high-performance liquid chromatography. F. heteroclitus and F. grandis exhibited an increase in gill Put concentration, but showed no change in Glu or Gaba levels following freshwater transfer. F. heteroclitus also accumulated Spd in the gills, whereas F. grandis showed transient increases in Spd and Spm levels. In contrast, gill Put, Spm, Glu, and Gaba levels decreased in F. majalis following freshwater transfer. Together, these findings suggest that increasing polyamine levels and maintaining Glu and Gaba levels in the gills may enable euryhaline teleosts to acclimate to shifts in environmental salinity.

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.

Reductionist approaches to the study of ionoregulation in fishes

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

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.

Interspecific variation in hypoxia tolerance and hypoxia acclimation responses in killifish from the family Fundulidae

Hypoxia is a pervasive stressor in aquatic environments, and both phenotypic plasticity and evolutionary adaptation could shape the ability to cope with hypoxia. We investigated evolved variation in hypoxia tolerance and the hypoxia acclimation response across fundulid killifishes that naturally experience different patterns of hypoxia exposure. We compared resting O2 consumption rate (ṀO2 ), and various indices of hypoxia tolerance [critical O2 tension (Pcrit), regulation index (RI), O2 tension (PO2 ) at loss of equilibrium (PLOE) and time to LOE (tLOE) at 0.6 kPa O2] in Fundulus confluentus, Fundulus diaphanus, Fundulus heteroclitus, Fundulus rathbuni, Lucania goodei and Lucania parva We examined the effects of chronic (28 days) exposure to constant hypoxia (2 kPa) or nocturnal intermittent hypoxia (12 h normoxia:12 h hypoxia) in a subset of species. Some species exhibited a two-breakpoint model in ṀO2  caused by early, modest declines in ṀO2  in moderate hypoxia. We found that hypoxia tolerance varied appreciably across species: F. confluentus was the most tolerant (lowest PLOE and Pcrit, longest tLOE), whereas F. rathbuni and F. diaphanus were the least tolerant. However, there was not a consistent pattern of interspecific variation for different indices of hypoxia tolerance, with or without taking phylogenetic relatedness into account, probably because these different indices are underlain by partially distinct mechanisms. Hypoxia acclimation generally improved hypoxia tolerance, but the magnitude of plasticity and responsiveness to different hypoxia patterns varied interspecifically. Our results therefore suggest that hypoxia tolerance is a complex trait that is best appreciated by considering multiple indices of tolerance.

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

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

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.

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.

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.

Gene expression differs in susceptible and resistant amphibians exposed to Batrachochytrium dendrobatidis

Chytridiomycosis, the disease caused by the fungal pathogen Batrachochytrium dendrobatidis (Bd), has devastated global amphibian biodiversity. Nevertheless, some hosts avoid disease after Bd exposure even as others experience near-complete extirpation. It remains unclear whether the amphibian adaptive immune system plays a role in Bd defence. Here, we describe gene expression in two host species-one susceptible to chytridiomycosis and one resistant-following exposure to two Bd isolates that differ in virulence. Susceptible wood frogs (Rana sylvatica) had high infection loads and mortality when exposed to the more virulent Bd isolate but lower infection loads and no fatal disease when exposed to the less virulent isolate. Resistant American bullfrogs (R. catesbeiana) had high survival across treatments and rapidly cleared Bd infection or avoided infection entirely. We found widespread upregulation of adaptive immune genes and downregulation of important metabolic and cellular maintenance components in wood frogs after Bd exposure, whereas American bullfrogs showed little gene expression change and no evidence of an adaptive immune response. Wood frog responses suggest that adaptive immune defences may be ineffective against virulent Bd isolates that can cause rapid physiological dysfunction. By contrast, American bullfrogs exhibited robust resistance to Bd that is likely attributable, at least in part, to their continued upkeep of metabolic and skin integrity pathways as well as greater antimicrobial peptide expression compared to wood frogs, regardless of exposure. Greater understanding of these defences will ultimately help conservationists manage chytridiomycosis.

Parallel habitat acclimatization is realized by the expression of different genes in two closely related salamander species (genus Salamandra)

The utilization of similar habitats by different species provides an ideal opportunity to identify genes underlying adaptation and acclimatization. Here, we analysed the gene expression of two closely related salamander species: Salamandra salamandra in Central Europe and Salamandra infraimmaculata in the Near East. These species inhabit similar habitat types: 'temporary ponds' and 'permanent streams' during larval development. We developed two species-specific gene expression microarrays, each targeting over 12 000 transcripts, including an overlapping subset of 8331 orthologues. Gene expression was examined for systematic differences between temporary ponds and permanent streams in larvae from both salamander species to establish gene sets and functions associated with these two habitat types. Only 20 orthologues were associated with a habitat in both species, but these orthologues did not show parallel expression patterns across species more than expected by chance. Functional annotation of a set of 106 genes with the highest effect size for a habitat suggested four putative gene function categories associated with a habitat in both species: cell proliferation, neural development, oxygen responses and muscle capacity. Among these high effect size genes was a single orthologue (14-3-3 protein zeta/YWHAZ) that was downregulated in temporary ponds in both species. The emergence of four gene function categories combined with a lack of parallel expression of orthologues (except 14-3-3 protein zeta) suggests that parallel habitat adaptation or acclimatization by larvae from S. salamandra and S. infraimmaculata to temporary ponds and permanent streams is mainly realized by different genes with a converging functionality.

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.

Kidney morphology and candidate gene expression shows plasticity in sticklebacks adapted to divergent osmotic environments

Novel physiological challenges in different environments can promote the evolution of divergent phenotypes, either through plastic or genetic changes. Environmental salinity serves as a key barrier to the distribution of nearly all aquatic organisms, and species diversification is likely to be enabled by adaptation to alternative osmotic environments. The threespine stickleback (Gasterosteus aculeatus) is a euryhaline species with populations found both in marine and freshwater environments. It has evolved both highly plastic and locally adapted phenotypes due to salinity-derived selection, but the physiological and genetic basis of adaptation to salinity is not fully understood. We integrated comparative cellular morphology of the kidney, a key organ for osmoregulation, and candidate gene expression to explore the underpinnings of evolved variation in osmotic plasticity within two populations of sticklebacks from distinct salinity zones in the Baltic Sea: the high salinity Kattegat, representative of the ancestral marine habitat, and the low salinity Bay of Bothnia. A common-garden experiment revealed that kidney morphology in the ancestral high salinity population had a highly plastic response to salinity conditions, whereas this plastic response was reduced in the low salinity population. Candidate gene expression in kidney tissue revealed a similar pattern of population-specific differences, with a higher degree of plasticity in the native high salinity population. Together these results suggest that renal cellular morphology has become canalized to low salinity, and that these structural differences may have functional implications for osmoregulation.

Transcriptomic differences between euryhaline and stenohaline malaria vector sibling species in response to salinity stress

Evolution of osmoregulatory systems is a key factor in the transition of species between fresh- and saltwater habitats. Anopheles coluzzii and Anopheles merus are stenohaline and euryhaline malaria vector mosquitoes belonging to a larger group of sibling species, the Anopheles gambiae complex, which radiated in Africa within the last 2 million years. Comparative ecological genomics of these vector species can provide insight into the mechanisms that permitted the rapid radiation of this species complex into habitats of contrasting salinity. Here, we use RNA-Seq to investigate gene expression differences between An. coluzzii and An. merus after briefly exposing both young and old larval instars of each species to either saltwater (SW) or freshwater (FW). Our study aims to identify candidate genes and pathways responsible for the greater SW tolerance of An. merus. Our results are congruent with the ability of gene induction to mediate salinity tolerance, with both species showing increasing amounts of differential gene expression between SW and FW as salt concentrations increase. Besides ion transporters such as AgAE2 that may serve as effectors for osmoregulation, we also find mitogen-activated protein kinases that may serve in a phosphorylation signalling pathway responding to salinity, and report potential cross-talk between the mosquito immune response and osmoregulation. This study provides a key step towards applying the growing molecular knowledge of these malaria vectors to improve understanding of their ecological tolerances and habitat occupancy.

Considerations for the use of transcriptomics in identifying the ‘genes that matter’ for environmental adaptation

Transcriptomics has emerged as a powerful approach for exploring physiological responses to the environment. However, like any other experimental approach, transcriptomics has its limitations. Transcriptomics has been criticized as an inappropriate method to identify genes with large impacts on adaptive responses to the environment because: (1) genes with large impacts on fitness are rare; (2) a large change in gene expression does not necessarily equate to a large effect on fitness; and (3) protein activity is most relevant to fitness, and mRNA abundance is an unreliable indicator of protein activity. In this review, these criticisms are re-evaluated in the context of recent systems-level experiments that provide new insight into the relationship between gene expression and fitness during environmental stress. In general, these criticisms remain valid today, and indicate that exclusively using transcriptomics to screen for genes that underlie environmental adaptation will overlook constitutively expressed regulatory genes that play major roles in setting tolerance limits. Standard practices in transcriptomic data analysis pipelines may also be limiting insight by prioritizing highly differentially expressed and conserved genes over those genes that undergo moderate fold-changes and cannot be annotated. While these data certainly do not undermine the continued and widespread use of transcriptomics within environmental physiology, they do highlight the types of research questions for which transcriptomics is best suited and the need for more gene functional analyses. Such information is pertinent at a time when transcriptomics has become increasingly tractable and many researchers may be contemplating integrating transcriptomics into their research programs.

Development and application of the adverse outcome pathway framework for understanding and predicting chronic toxicity: I. Challenges and research needs in ecotoxicology

To elucidate the effects of chemicals on populations of different species in the environment, efficient testing and modeling approaches are needed that consider multiple stressors and allow reliable extrapolation of responses across species. An adverse outcome pathway (AOP) is a concept that provides a framework for organizing knowledge about the progression of toxicity events across scales of biological organization that lead to adverse outcomes relevant for risk assessment. In this paper, we focus on exploring how the AOP concept can be used to guide research aimed at improving both our understanding of chronic toxicity, including delayed toxicity as well as epigenetic and transgenerational effects of chemicals, and our ability to predict adverse outcomes. A better understanding of the influence of subtle toxicity on individual and population fitness would support a broader integration of sublethal endpoints into risk assessment frameworks. Detailed mechanistic knowledge would facilitate the development of alternative testing methods as well as help prioritize higher tier toxicity testing. We argue that targeted development of AOPs supports both of these aspects by promoting the elucidation of molecular mechanisms and their contribution to relevant toxicity outcomes across biological scales. We further discuss information requirements and challenges in application of AOPs for chemical- and site-specific risk assessment and for extrapolation across species. We provide recommendations for potential extension of the AOP framework to incorporate information on exposure, toxicokinetics and situation-specific ecological contexts, and discuss common interfaces that can be employed to couple AOPs with computational modeling approaches and with evolutionary life history theory. The extended AOP framework can serve as a venue for integration of knowledge derived from various sources, including empirical data as well as molecular, quantitative and evolutionary-based models describing species responses to toxicants. This will allow a more efficient application of AOP knowledge for quantitative chemical- and site-specific risk assessment as well as for extrapolation across species in the future.

The effects of salinity on swimming performance of two estuarine fishes, Fundulus heteroclitus and Fundulus majalis

Prolonged and high-speed swimming performance measurements were used to explore the swimming abilities of two species of estuarine fishes, the mummichog Fundulus heteroclitus and the striped killifish Fundulus majalis, under different salinities. Critical swimming performance was significantly higher for F. majalis in high salinity than in low salinity, but no difference was observed in brief constant acceleration swimming trials in this species; however, the swimming performance of F. heteroclitus was not significantly affected by salinity changes, indicating that this species is well adapted to regular estuarine salinity oscillations. Fundulus majalis displayed higher swimming speeds than F. heteroclitus in both high and low salinities, and while this cannot be explained by their respective salinity preferences, the specific habitat preferences of F. majalis for sandy subtidal habitats and F. heteroclitus for vegetated marshes could explain the better swimming performance of F. majalis.

Ten years of transcriptomics in wild populations: what have we learned about their ecology and evolution?

Molecular ecology has moved beyond the use of a relatively small number of markers, often noncoding, and it is now possible to use whole-genome measures of gene expression with microarrays and RNAseq (i.e. transcriptomics) to capture molecular response to environmental challenges. While transcriptome studies are shedding light on the mechanistic basis of traits as complex as personality or physiological response to catastrophic events, these approaches are still challenging because of the required technical expertise, difficulties with analysis and cost. Still, we found that in the last 10 years, 575 studies used microarrays or RNAseq in ecology. These studies broadly address three questions that reflect the progression of the field: (i) How much variation in gene expression is there and how is it structured? (ii) How do environmental stimuli affect gene expression? (iii) How does gene expression affect phenotype? We discuss technical aspects of RNAseq and microarray technology, and a framework that leverages the advantages of both. Further, we highlight future directions of research, particularly related to moving beyond correlation and the development of additional annotation resources. Measuring gene expression across an array of taxa in ecological settings promises to enrich our understanding of ecology and genome function.

Reciprocal osmotic challenges reveal mechanisms of divergence in phenotypic plasticity in the killifish Fundulus heteroclitus

The killifish Fundulus heteroclitus is an estuarine species with broad physiological plasticity enabling acclimation to diverse stressors. Previous work suggests freshwater populations expanded their physiology to accommodate low salinity environments, however, it is unknown if this compromises their tolerance to high salinity. We employed a comparative approach to investigate the mechanisms of a derived freshwater phenotype and the fate of an ancestral euryhaline phenotype after invasion of a freshwater environment. We compared physiological and transcriptomic responses to high and low salinity stress in fresh and brackish water populations and found an enhanced plasticity to low salinity in the freshwater population coupled with a reduced ability to acclimate to high salinity. Transcriptomic data identified genes with a conserved common response, a conserved salinity dependent response, and responses associated with population divergence. Conserved common acclimation responses revealed stress responses and alterations in cell-cycle regulation as important mechanisms in the general osmotic response. Salinity-specific responses included the regulation of genes involved in ion transport, intracellular calcium, energetic processes, and cellular remodeling. Genes diverged between populations were primarily those showing salinity-specific expression and included those regulating polyamine homeostasis and cell cycle. Additionally, when populations were matched with their native salinity, expression patterns were consistent with the concept of “transcriptomic resilience,” suggesting local adaptation. These findings provide insight into the fate of a plastic phenotype after a shift in environmental salinity and help to reveal mechanisms allowing for euryhalinity.

Evolution of salinity tolerance from transcriptome to physiological system

The relationship between microevolution and macroevolution is a topic of fundamental importance in evolutionary biology. The increasing accessibility of genomic tools is making the hunt for genes that underlie evolutionary divergence more tractable and, when combined with physiological approaches, provides exceptional power to elucidate the causal mechanisms of the relationship. In this issue of Molecular Ecology, Whitehead et al. (2013) employ this strategy to show that common physiological and genomic mechanisms lead to divergence in salinity tolerance across micro- and macroevolutionary timescales. They compare two killifish species from the genus Fundulus, F. majalis, which inhabits primarily marine and brackish environments and represents the ancestral state of the genus, and F. heteroclitus, which has derived an osmotic niche that expands into freshwater. Corresponding to the differences in osmotic niche, the species differ strikingly in how the structure of the ion-transporting epithelium and the transcriptome of the gills respond to osmotic challenge. These inter-specific differences were similar to but more pronounced than the differences associated with the more subtle intra-specific variation in osmotic niche within each species. It appears that a progression of the same functional adjustments first allowed expansion of the osmotic niche of F. heteroclitus into freshwater and then further expanded the niche of select F. heteroclitus populations towards more dilute freshwater environments. The work of Whitehead et al. therefore emphasizes how the mechanisms of adaptive divergence between populations can be expanded over time to produce the more complex differences that can exist between species.

Proteome variance differences within populations of European whitefish (Coregonus lavaretus) originating from contrasting salinity environments

Variation in gene expression is an important component of the phenotypic differences observed in nature. Gene expression variance across biological groups and environmental conditions has been studied extensively and has revealed specific genes and molecular mechanisms of interest. However, little is known regarding the importance of within-population gene expression variation to environmental adaptation. To address this issue, we quantified the proteomes of individuals of European whitefish (Coregonus lavaretus) from populations that have previously been shown to have adapted during early development to freshwater and brackishwater salinity environments. Using MS-based label-free proteomics, we studied 955 proteins in eight hatch-stage fish embryos from each population that had been reared in either freshwater or brackishwater salinity conditions. By comparing the levels of within-population protein expression variance over individuals and per protein between populations, we found that fish embryos from the population less affected by salinity level had also markedly higher levels of expression variance. Gene Ontologies and molecular pathways associated with osmoregulation showed the most significant difference of within-population proteome variance between populations. Several new candidate genes for salinity adaptation were identified, emphasising the added value of combining assessments of within-population gene expression variation with standard gene expression analysis practices for better understanding the mechanisms of environmental adaptation.

BIOLOGICAL SIGNIFICANCE

We demonstrate the benefits of studying within-population gene expression variance together with more typical methods of gene expression profiling. Proteome variance differences within European whitefish populations originating from different salinity environments allowed us to identify several new candidate genes for salinity adaptation in teleost fish and generate many further hypotheses to be tested. This article is part of a Special Issue entitled: Proteomics of non-model organisms.

Conservation physiology across scales: insights from the marine realm

As the field of conservation physiology develops and becomes increasingly integrated with ecology and conservation science, the fundamental concept of scale is being recognized as important, particularly for ensuring that physiological knowledge is contextualized in a manner most relevant to policy makers, conservation practitioners and stakeholders. Failure to consider the importance of scale in conservation physiology-both the challenges and the opportunities that it creates-will impede the ability of this discipline to generate the scientific understanding needed to contribute to meaningful conservation outcomes. Here, we have focused on five aspects of scale: biological, spatial, temporal, allometric and phylogenetic. We also considered the scale of policy and policy application relevant to those five types of scale as well as the merits of upscaling and downscaling to explore and address conservation problems. Although relevant to all systems (e.g. freshwater, terrestrial) we have used examples from the marine realm, with a particular emphasis on fishes, given the fact that there is existing discourse regarding scale and its relevance for marine conservation and management. Our synthesis revealed that all five aspects of scale are relevant to conservation physiology, with many aspects inherently linked. It is apparent that there are both opportunities and challenges afforded by working across scales but, to understand mechanisms underlying conservation problems, it is essential to consider scale of all sorts and to work across scales to the greatest extent possible. Moreover, given that the scales in biological processes will often not match policy and management scales, conservation physiology needs to show how it is relevant to aspects at different policy/management scales, change the scales at which policy/management intervention is applied or be prepared to be ignored.

Functional and population genomic divergence within and between two species of killifish adapted to different osmotic niches

Adaptation to salinity affects species distributions, promotes speciation, and guides many evolutionary patterns in fishes. To uncover the basis of a complex trait like osmoregulation, genome-level analyses are sensible. We combine population genomic scans with genome expression profiling to discover candidate genes and pathways associated with divergence between osmotic environments. We compared transcriptome sequence divergence between multiple freshwater and saltwater populations of the rainwater killifish, Lucania parva. We also compared sequence divergence between L. parva and its sister species, Lucania goodei, a freshwater specialist. We found highly differentiated single nucleotide polymorphisms (SNPs) between freshwater and saltwater L. parva populations in cell junction and ion transport genes, including V-type H(+) ATPase. Between species, we found divergence in reproduction and osmotic stress genes. Genes that were differentially expressed between species during osmotic acclimation included genes involved in ion transport and cell volume regulation. Gene sets that were divergent in coding sequence and divergent in expression did not overlap, although they did converge in function. Like many studies using genomic scans, our approach may miss some loci that contribute to adaptation but have complicated patterns of allelic variation. Our study suggests that gene expression and coding sequence may evolve independently as populations adapt to a complex physiological challenge.