Genetic Diversity and Signatures of Selection in the Roughskin Sculpin (Trachidermus fasciatus) Revealed by Whole Genome Sequencing

The roughskin sculpin (Trachidermus fasciatus) is an endangered fish species in China. In recent years, artificial breeding technology has made significant progress, and the population of roughskin sculpin has recovered in the natural environment through enhancement programs and the release of juveniles. However, the effects of released roughskin sculpin on the genetic structure and diversity of wild populations remain unclear. Studies on genetic diversity analysis based on different types and numbers of molecular markers have yielded inconsistent results. In this study, we obtained 2,610,157 high-quality SNPs and 494,698 InDels through whole-genome resequencing of two farmed populations and one wild population. Both farmed populations showed consistent levels of genomic polymorphism and a slight increase in linkage compared with wild populations. The population structure of the two farmed populations was distinct from that of the wild population, but the degree of genetic differentiation was low (overall average Fst = 0.015). Selective sweep analysis showed that 523,529 genes were selected in the two farmed populations, and KEGG enrichment analysis showed that the selected genes were related to amino acid metabolism, which might be caused by artificial feeding. The findings of this study provide valuable additions to the existing genomic resources to help conserve roughskin sculpin populations.

Morphological and transcriptomic responses/acclimations of erect-type submerged macrophyte Hydrilla verticillata both at low-light exposure and light recovery phases

Light intensity is a determinant for submerged macrophytes. Little is known about their molecular responses to low-light exposure, despite more informative and responsive than morphological traits. For erect-type submerged macrophytes, the stem is more crucial relative to the leaf in acclimation to low-light stress, but receives less attention. We determined morphological and stem transcriptomic responses/acclimations of Hydrilla verticillata to extremely and mildly low light (7.2 and 36 μmol photons m-2 s-1, respectively), that is, EL and ML, with the radiation intensity of 180 μmol photons m-2 s-1 as the control. Low-light exposure continued for 9 days, followed by a 7-day recovery phase (180 μmol photons m-2 s-1). At the exposure phase, the low-light treatments, in particular the EL, decreased the relative growth ratio, but induced greater height and longer stem internode distance and epidermal cell. Such responses/acclimations continued into the recovery phase, despite more or less changes in the magnitude. Transcriptome showed that the photosynthetic system was inhibited at the exposure phase, but the macrophyte adjusted hormone synthesis relating to cell division and elongation. Moreover, the EL activated cell stress responses such as DNA repair. Following light recovery, the macrophyte exhibited a strong-light response, although energy metabolism enhanced. Especially, the EL enriched the pathways relating to anthocyanin synthesis at such phase, indicating an activation of photoprotective mechanism. Our findings suggest that negative influences of low light occur at both low-light exposure and recovery phases, but submerged macrophytes would acclimate to light environments. Transcriptome can show molecular basis of plant responses/acclimations, including but not limited to morphology. This study establishes a bridge connecting morphological and molecular responses/acclimations.

Metabolic rate and osmoregulation in desert spring amphipods

Desert springs contain many endemic taxa and are of conservation concern due to anthropogenic activities that are expected to increase environmental salinity. Understanding the nature of osmoregulation is necessary to predict how non-vagile aquatic organisms will respond to changes. In the Chihuahuan Desert, the amphipod genus Gammarus Fabricius, 1775 is composed of two lineages. These lineages have species that currently inhabit springs ranging from 0.4 to 7.8 parts per thousand (ppt). All Gammarus in this region are of conservation concern because each is endemic to a single spring system. We exposed individuals of Gammarus colei Walters, Cannizzaro, and Berg in Walters, Cannizzaro, Trujillo and Berg, 2020 and Gammarus seideli Cannizzaro, Walters and Berg, 2018 species occupying low-salinity springs, to the range of salinities found in the Chihuahuan Desert. We measured metabolic rates as [Formula: see text] to examine the energetic cost of osmoregulation. We also measured the hemolymph osmolality of G. colei and compared that with the isosmotic line to determine the degree of osmoregulation. Neither species increased its metabolic rate across increasing salinities. However, G. colei showed an increase in hemolymph osmolality. Despite the divergence (∼66 million years ago) between the two lineages, it appears their physiological tolerances have converged, suggesting that both lineages may be able to persist when exposed to moderate salinity changes.

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.

Why sequence all eukaryotes?

Life on Earth has evolved from initial simplicity to the astounding complexity we experience today. Bacteria and archaea have largely excelled in metabolic diversification, but eukaryotes additionally display abundant morphological innovation. How have these innovations come about and what constraints are there on the origins of novelty and the continuing maintenance of biodiversity on Earth? The history of life and the code for the working parts of cells and systems are written in the genome. The Earth BioGenome Project has proposed that the genomes of all extant, named eukaryotes-about 2 million species-should be sequenced to high quality to produce a digital library of life on Earth, beginning with strategic phylogenetic, ecological, and high-impact priorities. Here we discuss why we should sequence all eukaryotic species, not just a representative few scattered across the many branches of the tree of life. We suggest that many questions of evolutionary and ecological significance will only be addressable when whole-genome data representing divergences at all of the branchings in the tree of life or all species in natural ecosystems are available. We envisage that a genomic tree of life will foster understanding of the ongoing processes of speciation, adaptation, and organismal dependencies within entire ecosystems. These explorations will resolve long-standing problems in phylogenetics, evolution, ecology, conservation, agriculture, bioindustry, and medicine.

Next-Generation Sequencing and Its Impacts on Entomological Research in Ecology and Evolution

The advent of NGS-based methods has been profoundly transforming entomological research. Through continual development and improvement of different methods and sequencing platforms, NGS has promoted mass elucidation of partial or whole genetic materials associated with beneficial insects, pests (of agriculture, forestry and animal, and human health), and species of conservation concern, helping to unravel ecological and evolutionary mechanisms and characterizing survival, trophic interactions, and dispersal. It is shifting the scale of biodiversity and environmental analyses from individuals and biodiversity indicator species to the large-scale study of communities and ecosystems using bulk samples of species or a mixed "soup" of environmental DNA. As the NGS-based methods have become more affordable, complexity demystified, and specificity and sensitivity proven, their use in entomological research has spread widely. This article presents several examples on how NGS-based methods have been used in entomology to provide incentives to apply them when appropriate and to open our minds to the expected advances in entomology that are yet to come.

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

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

The comparative genomic landscape of adaptive radiation in crater lake cichlid fishes

Factors ranging from ecological opportunity to genome composition might explain why only some lineages form adaptive radiations. While being rare, particular systems can provide natural experiments within an identical ecological setting where species numbers and phenotypic divergence in two closely related lineages are notably different. We investigated one such natural experiment using two de novo assembled and 40 resequenced genomes and asked why two closely related Neotropical cichlid fish lineages, the Amphilophus citrinellus species complex (Midas cichlids; radiating) and Archocentrus centrarchus (Flyer cichlid; nonradiating), have resulted in such disparate evolutionary outcomes. Although both lineages inhabit many of the same Nicaraguan lakes, whole-genome inferred demography suggests that priority effects are not likely to be the cause of the dissimilarities. Also, genome-wide levels of selection, transposable element dynamics, gene family expansion, major chromosomal rearrangements and the number of genes under positive selection were not markedly different between the two lineages. To more finely investigate particular subsets of the genome that have undergone adaptive divergence in Midas cichlids, we also examined if there was evidence for 'molecular pre-adaptation' in regions identified by QTL mapping of repeatedly diverging adaptive traits. Although most of our analyses failed to pinpoint substantial genomic differences, we did identify functional categories containing many genes under positive selection that provide candidates for future studies on the propensity of Midas cichlids to radiate. Our results point to a disproportionate role of local, rather than genome-wide factors underlying the propensity for these cichlid fishes to adaptively radiate.

Climate explains population divergence in drought-induced plasticity of functional traits and gene expression in a South African Protea

Long-term environmental variation often drives local adaptation and leads to trait differentiation across populations. Additionally, when traits change in an environment-dependent way through phenotypic plasticity, the genetic variation underlying plasticity will also be under selection. These processes could create a landscape of differentiation across populations in traits and their plasticity. Here, we performed a dry-down experiment under controlled conditions to measure responses in seedlings of a shrub species from the Cape Floristic Region, the common sugarbush (Protea repens). We measured morphological and physiological traits, and sequenced whole transcriptomes of leaf tissues from eight populations that represent both the climatic and the geographical distribution of this species. We found that there is substantial variation in how populations respond to drought, but we also observed common patterns such as reduced leaf size and leaf thickness, and up-regulation of stress-related and down-regulation of growth-related gene groups. Both high environmental heterogeneity and milder source site climates were associated with higher plasticity in various traits and co-expression gene networks. Associations between traits, trait plasticity, gene networks and the source site climate suggest that temperature may play a greater role in shaping these patterns when compared to precipitation, in line with recent changes in the region due to climate change. We also found that traits respond to climatic variation in an environment-dependent manner: some associations between traits and climate were apparent only under certain growing conditions. Together, our results uncover common responses of P. repens populations to drought, and climatic drivers of population differentiation in functional traits, gene expression and their plasticity.

Investigating population-scale allelic differential expression in wild populations of Oithona similis (Cyclopoida, Claus, 1866)

Acclimation allowed by variation in gene or allele expression in natural populations is increasingly understood as a decisive mechanism, as much as adaptation, for species evolution. However, for small eukaryotic organisms, as species from zooplankton, classical methods face numerous challenges. Here, we propose the concept of allelic differential expression at the population-scale (psADE) to investigate the variation in allele expression in natural populations. We developed a novel approach to detect psADE based on metagenomic and metatranscriptomic data from environmental samples. This approach was applied on the widespread marine copepod, Oithona similis, by combining samples collected during the Tara Oceans expedition (2009-2013) and de novo transcriptome assemblies. Among a total of 25,768 single nucleotide variants (SNVs) of O. similis, 572 (2.2%) were affected by psADE in at least one population (FDR < 0.05). The distribution of SNVs under psADE in different populations is significantly shaped by population genomic differentiation (Pearson r = 0.87, p = 5.6 × 10-30), supporting a partial genetic control of psADE. Moreover, a significant amount of SNVs (0.6%) were under both selection and psADE (p < .05), supporting the hypothesis that natural selection and psADE tends to impact common loci. Population-scale allelic differential expression offers new insights into the gene regulation control in populations and its link with natural selection.

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

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

Comparative Transcriptomics Reveals Patterns of Adaptive Evolution Associated with Depth and Age Within Marine Rockfishes (Sebastes)

The genetic underpinnings that contribute to ecological adaptation and speciation are not completely understood, especially within marine ecosystems. These evolutionary processes can be elucidated by studying adaptive radiations, because they provide replicates of divergence within a given environment or time-frame. Marine rockfishes (genus Sebastes) are an adaptive radiation and unique model system for studying adaptive evolution in the marine realm. We investigated molecular evolution associated with ecological (depth) and life history (lifespan) divergence in 2 closely related clades of Sebastes. Brain transcriptomes were sequenced via RNA-Seq from 3 species within the subgenus Pteropodus and a pair of related congeners from the subgenus Sebastosomus in order to identify patterns of adaptive evolution. De novo assemblies from these transcriptomes were used to identify 3867 orthologous clusters, and genes subject to positive selection were identified based on all 5 species, depth, and lifespan. Within all our analyses, we identified hemoglobin subunit α to be under strong positive selection and is associated with the depth of occurrence. In our lifespan analysis we identified immune function genes under positive selection in association with maximum lifespan. This study provides insight on the molecular evolution of rockfishes and these candidate genes may provide a better understanding of how these subgenera radiated within the Northeast Pacific.

Genome-Wide Microarray Analysis Suggests Transcriptomic Response May Not Play a Major Role in High- to Low-Altitude Acclimation in Harvest Mouse (Micromys minutus)

Simple Summary

Micromys minutus is a small rodent species that has a wide range of vertical distribution in Taiwan. By comparing the gene expression profile of the skeletal muscle tissues taken from individuals native to the high-altitude environment and those transferred to the low-altitude captive site, the Tnfrsf12a gene was demonstrated to have a differential expression pattern. Although this finding may be correlated with the altitude acclimation, the observation of only one gene transcript with significant alteration leads us to suggest that genetic response may not play a major role in altitude acclimation in M. minutus. Future comparative functional genomics studies involving other organ systems (in addition to skeletal muscles) and alarger sample size are warranted for better insight into the altitude acclimation of this small rodent species.

Abstract

The harvest mouse (Micromys minutus) is a small rodent species with a wide range of vertical distribution in Taiwan, extending from the sea level to 3100 m altitude. This species has recently suffered from habitat loss in high-altitude areas due to orchard cultivation, which may have resulted in mouse migration from high to low altitude. To investigate whether there is any physiological mechanism involved in altitude acclimation, rat cDNA microarray was used to compare transcriptomic patterns of the skeletal muscle tissues taken from individuals native to the high-altitude environment and those transferred to the low-altitude captive site. Of the 23,188 genes being analyzed, 47 (33 up-regulated and 14 down-regulated) were found to have differential expression (fold change > 4 or < −4, ANOVA p < 0.05). However, after multiple testing correction with a false discovery rate (FDR), only the result for Tnfrsf12a was found to be statistically significant (fold change = 13, FDR p < 0.05). The result was confirmed by quantitative polymerase chain reaction (q-PCR). The expression of Tnfrsf12a possibly relates to the skeletal muscle biology and thus can be correlated with altitude acclimation. However, finding only one gene transcript with significant alteration suggests that transcriptomic response may not play a major role in high- to low-altitude acclimation in harvest mouse.

Selection of suitable reference genes for gene expression analysis in gills and liver of fish under field pollution conditions

To understand the role of gene expression in adaptive variation, it is necessary to examine expression variation in an ecological context. Quantitative real-time PCR (qPCR) is considered the most accurate and reliable technique to measure gene expression and to validate the data obtained by RNA-seq; however, accurate normalization is crucial. In Chile, the freshwater silverside fish Basilichthys microlepidotus inhabits both polluted and nonpolluted areas, showing differential gene expression related to pollution. In this study, we infer the stability of six potential reference genes (tubulin alpha, hypoxanthine-guanine phosphoribosyltransferase, glyceraldehyde-3-phosphate dehydrogenase, beta-actin, 60S ribosomal protein L13, and 60S ribosomal protein L8) in the gills and liver of silverside individuals inhabiting polluted and nonpolluted areas. To validate the reference genes selected, the most and least stable reference genes were used to normalize two target transcripts, one for each organ. The RefFinder tool was used to analyze and identify the most stably expressed genes. The 60S ribosomal protein L8 gene was ranked as the most stable gene for both organs. Our results show that reference gene selection influences the detection of differences in the expression levels of target genes in different organs and, also highlighting candidate reference genes that could be used in field studies.

Variable gene transcription underlies phenotypic convergence of hypoxia tolerance in sculpins

Background

The degree by which mechanisms underlying phenotypic convergence are similar among taxa depends on the number of evolutionary paths available for selection to act upon. Likelihood of convergence will be influenced by an interplay of factors such as genetic architecture, phylogenetic history and population demography. To determine if there is convergence or divergence in mechanisms underlying phenotypic similarity, we assessed whether gene transcription patterns differed among species with similar levels of hypoxia tolerance.

Results

Three species of marine fish from the superfamily Cottoidea (smoothhead sculpin [Artedius lateralis], sailfin sculpin [Nautichthys oculofasciatus] and Pacific staghorn sculpin [Leptocottus armatus]), all of which have previously been shown to share the same level of hypoxia tolerance, were exposed to short-(8 h) and longer-term (72 h) hypoxia and mRNA transcripts were assessed using a custom microarray. We examined hypoxia-induced transcription patterns in metabolic and protein production pathways and found that a high proportion of genes associated with these biological processes showed significant differences among the species. Specifically, the data suggest that the smoothhead sculpin, unlike the sailfin sculpin and the Pacific staghorn sculpin, relied on amino acid degradation rather than glycolysis or fatty acid oxidation to generate ATP during hypoxia exposure. There was also variation across the species in the transcription of genes involved in protein production (e.g. mRNA processing and protein translation), such that it increased in the smoothhead sculpin, decreased in the sailfin sculpin and was variable in the Pacific staghorn sculpin.

Conclusions

Changes in metabolic and protein production pathways are part of the key responses of fishes to exposures to environmental hypoxia. Yet, species with similar overall hypoxia tolerance exhibited different transcriptional responses in these pathways, indicating flexibility and complexity of interactions in the evolution of the mechanisms underlying the hypoxia tolerance phenotype. The variation in the hypoxia-induced transcription of genes across species with similar hypoxia tolerance suggests that similar whole-animal phenotypes can emerge from divergent evolutionary paths that may affect metabolically important functions.

Electronic supplementary material

The online version of this article (10.1186/s12862-018-1275-1) contains supplementary material, which is available to authorized users.

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.

Dissecting the Transcriptional Patterns of Social Dominance across Teleosts

In many species, under varying ecological conditions, social interactions among individuals result in the formation of dominance hierarchies. Despite general similarities, there are robust differences among dominance hierarchies across species, populations, environments, life stages, sexes, and individuals. Understanding the proximate mechanisms underlying the variation is an important step toward understanding the evolution of social behavior. However, physiological changes associated with dominance, such as gonadal maturation and somatic growth, often complicate efforts to identify the specific underlying mechanisms. Traditional gene expression analyses are useful for generating candidate gene lists, but are biased by choice of significance cut-offs and difficult to use for between-study comparisons. In contrast, complementary analysis tools allow one to both test a priori hypotheses and generate new hypotheses. Here we employ a meta-analysis of high-throughput expression profiling experiments to investigate the gene expression patterns that underlie mechanisms and evolution of behavioral social phenotypes. Specifically, we use a collection of datasets on social dominance in fish across social contexts, sex, and species. Using experimental manipulation to produce female dominance hierarchies in the cichlid Astatotilapia burtoni, heralded as a genomic model of social dominance, we generate gene lists, and assess molecular gene modules. In the dominant female gene expression profile, we demonstrate a strong pattern of up-regulation of genes previously identified as having male-biased expression and furthermore, compare expression biases between male and female dominance phenotypes. Using a threshold-free approach to identify correlation throughout ranked gene lists, we query previously published datasets associated with maternal behavior, alternative reproductive tactics, cooperative breeding, and sex-role reversal to describe correlations among these various neural gene expression profiles associated with different instances of social dominance. These complementary approaches capitalize on the high-throughput gene expression profiling from similar behavioral phenotypes in order to address the mechanisms associated with social dominance behavioral phenotypes.

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

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

The roles of plasticity and evolutionary change in shaping gene expression variation in natural populations of extremophile fish

The notorious plasticity of gene expression responses and the complexity of environmental gradients complicate the identification of adaptive differences in gene regulation among populations. We combined transcriptome analyses in nature with common-garden and exposure experiments to establish cause-effect relationships between the presence of a physiochemical stressor and expression differences, as well as to test how evolutionary change and plasticity interact to shape gene expression variation in natural systems. We studied two evolutionarily independent population pairs of an extremophile fish (Poecilia mexicana) living in toxic, hydrogen sulphide (H₂ S)-rich springs and adjacent nontoxic habitats and assessed genomewide expression patterns of wild-caught and common-garden-raised individuals exposed to different concentrations of H₂ S. We found that 7.7% of genes that were differentially expressed between sulphidic and nonsulphidic ecotypes remained differentially expressed in the laboratory, indicating that sources of selection other than H₂ S-or plastic responses to other environmental factors-contribute substantially to gene expression patterns observed in the wild. Concordantly differentially expressed genes in the wild and the laboratory were primarily associated with H₂ S detoxification, sulphur processing and metabolic physiology. While shared, ancestral plasticity played a minor role in shaping gene expression variation observed in nature, we documented evidence for evolved population differences in the constitutive expression as well as the H₂ S inducibility of candidate genes. Mechanisms underlying gene expression variation also varied substantially across the two ecotype pairs. These results provide a springboard for studying evolutionary modifications of gene regulatory mechanisms that underlie expression variation in locally adapted populations.

Complexities of gene expression patterns in natural populations of an extremophile fish (Poecilia mexicana, Poeciliidae)

Variation in gene expression can provide insights into organismal responses to environmental stress and physiological mechanisms mediating adaptation to habitats with contrasting environmental conditions. We performed an RNA-sequencing experiment to quantify gene expression patterns in fish adapted to habitats with different combinations of environmental stressors, including the presence of toxic hydrogen sulphide (H2 S) and the absence of light in caves. We specifically asked how gene expression varies among populations living in different habitats, whether population differences were consistent among organs, and whether there is evidence for shared expression responses in populations exposed to the same stressors. We analysed organ-specific transcriptome-wide data from four ecotypes of Poecilia mexicana (nonsulphidic surface, sulphidic surface, nonsulphidic cave and sulphidic cave). The majority of variation in gene expression was correlated with organ type, and the presence of specific environmental stressors elicited unique expression differences among organs. Shared patterns of gene expression between populations exposed to the same environmental stressors increased with levels of organismal organization (from transcript to gene to physiological pathway). In addition, shared patterns of gene expression were more common between populations from sulphidic than populations from cave habitats, potentially indicating that physiochemical stressors with clear biochemical consequences can constrain the diversity of adaptive solutions that mitigate their adverse effects. Overall, our analyses provided insights into transcriptional variation in a unique system, in which adaptation to H2 S and darkness coincide. Functional annotations of differentially expressed genes provide a springboard for investigating physiological mechanisms putatively underlying adaptation to extreme environments.

Genomic Characterization of the Evolutionary Potential of the Sea Urchin Strongylocentrotus droebachiensis Facing Ocean Acidification

Ocean acidification (OA) is increasing due to anthropogenic CO₂ emissions and poses a threat to marine species and communities worldwide. To better project the effects of acidification on organisms’ health and persistence, an understanding is needed of the 1) mechanisms underlying developmental and physiological tolerance and 2) potential populations have for rapid evolutionary adaptation. This is especially challenging in nonmodel species where targeted assays of metabolism and stress physiology may not be available or economical for large-scale assessments of genetic constraints. We used mRNA sequencing and a quantitative genetics breeding design to study mechanisms underlying genetic variability and tolerance to decreased seawater pH (-0.4 pH units) in larvae of the sea urchin Strongylocentrotus droebachiensis. We used a gene ontology-based approach to integrate expression profiles into indirect measures of cellular and biochemical traits underlying variation in larval performance (i.e., growth rates). Molecular responses to OA were complex, involving changes to several functions such as growth rates, cell division, metabolism, and immune activities. Surprisingly, the magnitude of pH effects on molecular traits tended to be small relative to variation attributable to segregating functional genetic variation in this species. We discuss how the application of transcriptomics and quantitative genetics approaches across diverse species can enrich our understanding of the biological impacts of climate change.

Plasticity in gene transcription explains the differential performance of two invasive fish species

Phenotypic plasticity buffers organisms from environmental change and is hypothesized to aid the initial establishment of nonindigenous species in novel environments and postestablishment range expansion. The genetic mechanisms that underpin phenotypically plastic traits are generally poorly characterized; however, there is strong evidence that modulation of gene transcription is an important component of these responses. Here, we use RNA sequencing to examine the transcriptional basis of temperature tolerance for round and tubenose goby, two nonindigenous fish species that differ dramatically in the extent of their Great Lakes invasions despite similar invasion dates. We used generalized linear models of read count data to compare gene transcription responses of organisms exposed to increased and decreased water temperature from those at ambient conditions. We identify greater response in the magnitude of transcriptional changes for the more successful round goby compared with the less successful tubenose goby. Round goby transcriptional responses reflect alteration of biological function consistent with adaptive responses to maintain or regain homeostatic function in other species. In contrast, tubenose goby transcription patterns indicate a response to stressful conditions, but the pattern of change in biological functions does not match those expected for a return to homeostatic status. Transcriptional plasticity plays an important role in the acute thermal tolerance for these species; however, the impaired response to stress we demonstrate in the tubenose goby may contribute to their limited invasion success relative to the round goby. Transcriptional profiling allows the simultaneous assessment of the magnitude of transcriptional response as well as the biological functions involved in the response to environmental stress and is thus a valuable approach for evaluating invasion potential.

Discovery of Nigri/nox and Panto/pox site-specific recombinase systems facilitates advanced genome engineering

Precise genome engineering is instrumental for biomedical research and holds great promise for future therapeutic applications. Site-specific recombinases (SSRs) are valuable tools for genome engineering due to their exceptional ability to mediate precise excision, integration and inversion of genomic DNA in living systems. The ever-increasing complexity of genome manipulations and the desire to understand the DNA-binding specificity of these enzymes are driving efforts to identify novel SSR systems with unique properties. Here, we describe two novel tyrosine site-specific recombination systems designated Nigri/nox and Panto/pox. Nigri originates from Vibrio nigripulchritudo (plasmid VIBNI_pA) and recombines its target site nox with high efficiency and high target-site selectivity, without recombining target sites of the well established SSRs Cre, Dre, Vika and VCre. Panto, derived from Pantoea sp. aB, is less specific and in addition to its native target site, pox also recombines the target site for Dre recombinase, called rox. This relaxed specificity allowed the identification of residues that are involved in target site selectivity, thereby advancing our understanding of how SSRs recognize their respective DNA targets.

Integrated studies of organismal plasticity through physiological and transcriptomic approaches: examples from marine polar regions

Transcriptomic methods are now widely used in functional genomic research. The vast amount of information received from these studies comes along with the challenge of developing a precise picture of the functional consequences and the characteristic regulatory mechanisms. Here we assess recent studies in marine species and their adaptation to polar (and seasonal) cold and explore how they have been able to draw reliable conclusions from transcriptomic patterns on functional consequences in the organisms. Our analysis indicates that the interpretation of transcriptomic data suffers from insufficient understanding of the consequences for whole organism performance and fitness and comes with the risk of supporting only preliminary and superficial statements.We propose that the functional understanding of transcriptomic data may be improved by their tighter integration into overarching physiological concepts that support the more specific interpretation of the 'omics' data and, at the same time, can be developed further through embedding the transcriptomic phenomena observed. Such possibilities have not been fully exploited.In the context of thermal adaptation and limitation, we explore preliminary evidence that the concept of oxygen and capacity limited thermal tolerance (OCLTT) may provide sufficient complexity to guide the integration of such data and the development of associated functional hypotheses. At the same time, we identify a lack of methodological approaches linking genes and function to higher levels of integration, in terms of organism and ecosystem functioning, at temporal and geographical scales, to support more reliable conclusions and be predictive with respect to the effects of global changes.

Transcriptomic plasticity in brown adipose tissue contributes to an enhanced capacity for nonshivering thermogenesis in deer mice

For small mammals living at high altitude, aerobic heat generation (thermogenesis) is essential for survival during prolonged periods of cold, but is severely impaired under conditions of hypobaric hypoxia. Recent studies in deer mice (Peromyscus maniculatus) reveal adaptive enhancement of thermogenesis in high- compared to low-altitude populations under hypoxic cold stress, an enhancement that is attributable to modifications in the aerobic metabolism of muscles used in shivering. However, because small mammals rely heavily on nonshivering mechanisms for cold acclimatization, we tested for evidence of adaptive divergence in nonshivering thermogenesis (NST) under hypoxia. To do so, we measured NST and characterized transcriptional profiles of brown adipose tissue (BAT) in high- and low-altitude deer mice that were (i) wild-caught and acclimatized to their native altitude, and (ii) born and reared under common garden conditions at low elevation. We found that NST performance under hypoxia is enhanced in wild-caught, high-altitude deer mice, a difference that is associated with increased expression of coregulated genes that influence several physiological traits. These traits include vascularization and O2 supply to BAT, brown adipocyte proliferation and the uncoupling of oxidative phosphorylation from ATP synthesis in the generation of heat. Our results suggest that acclimatization to hypoxic cold stress is facilitated by enhancement of nonshivering heat production, which is driven by regulatory plasticity in a suite of genes that influence intersecting physiological pathways.

Gene expression patterns of two dominant tallgrass prairie species differ in response to warming and altered precipitation

To better understand the mechanisms underlying plant species responses to climate change, we compared transcriptional profiles of the co-dominant C4 grasses, Andropogon gerardii Vitman and Sorghastrum nutans (L.) Nash, in response to increased temperatures and more variable precipitation regimes in a long-term field experiment in native tallgrass prairie. We used microarray probing of a closely related model species (Zea mays) to assess correlations in leaf temperature (Tleaf) and leaf water potential (LWP) and abundance changes of ~10,000 transcripts in leaf tissue collected from individuals of both species. A greater number of transcripts were found to significantly change in abundance levels with Tleaf and LWP in S. nutans than in A. gerardii. S. nutans also was more responsive to short-term drought recovery than A. gerardii. Water flow regulating transcripts associated with stress avoidance (e.g., aquaporins), as well as those involved in the prevention and repair of damage (e.g., antioxidant enzymes, HSPs), were uniquely more abundant in response to increasing Tleaf in S. nutans. The differential transcriptomic responses of the co-dominant C4 grasses suggest that these species may cope with and respond to temperature and water stress at the molecular level in distinct ways, with implications for tallgrass prairie ecosystem function.

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

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

New view of population genetics of zooplankton: RAD-seq analysis reveals population structure of the North Atlantic planktonic copepod Centropages typicus

Detection of population genetic structure of zooplankton at medium-to-small spatial scales in the absence of physical barriers has remained challenging and controversial. The large population sizes and high rates of gene flow characteristic of zooplankton have made resolution of geographical differentiation very difficult, especially when using few genetic markers and assuming equilibrium conditions. Next-generation sequencing now allows simultaneous sampling of hundreds to thousands of genetic markers; new analytical approaches allow studies under nonequilibrium conditions and directional migration. Samples of the North Atlantic Ocean planktonic copepod, Centropages typicus, were analysed using restriction site-associated DNA (RAD) sequencing on a PROTON platform. Although prior studies revealed no genetic differentiation of populations across the geographical range of the species, analysis of RAD tags showed significant structure across the North Atlantic Ocean. We also compared the likelihood for models of connectivity among NW Atlantic populations under various directional flow scenarios that replicate oceanographic conditions of the sampled domain. High-density marker sampling with RAD sequencing markedly outperformed other technical and analytical approaches in detection of population genetic structure and characterization of connectivity of this high geneflow zooplankton species.

Galaxy-M: a Galaxy workflow for processing and analyzing direct infusion and liquid chromatography mass spectrometry-based metabolomics data

BACKGROUND

Metabolomics is increasingly recognized as an invaluable tool in the biological, medical and environmental sciences yet lags behind the methodological maturity of other omics fields. To achieve its full potential, including the integration of multiple omics modalities, the accessibility, standardization and reproducibility of computational metabolomics tools must be improved significantly.

RESULTS

Here we present our end-to-end mass spectrometry metabolomics workflow in the widely used platform, Galaxy. Named Galaxy-M, our workflow has been developed for both direct infusion mass spectrometry (DIMS) and liquid chromatography mass spectrometry (LC-MS) metabolomics. The range of tools presented spans from processing of raw data, e.g. peak picking and alignment, through data cleansing, e.g. missing value imputation, to preparation for statistical analysis, e.g. normalization and scaling, and principal components analysis (PCA) with associated statistical evaluation. We demonstrate the ease of using these Galaxy workflows via the analysis of DIMS and LC-MS datasets, and provide PCA scores and associated statistics to help other users to ensure that they can accurately repeat the processing and analysis of these two datasets. Galaxy and data are all provided pre-installed in a virtual machine (VM) that can be downloaded from the GigaDB repository. Additionally, source code, executables and installation instructions are available from GitHub.

CONCLUSIONS

The Galaxy platform has enabled us to produce an easily accessible and reproducible computational metabolomics workflow. More tools could be added by the community to expand its functionality. We recommend that Galaxy-M workflow files are included within the supplementary information of publications, enabling metabolomics studies to achieve greater reproducibility.

Mechanisms Underlying Adaptation to Life in Hydrogen Sulfide-Rich Environments

Hydrogen sulfide (H2S) is a potent toxicant interfering with oxidative phosphorylation in mitochondria and creating extreme environmental conditions in aquatic ecosystems. The mechanistic basis of adaptation to perpetual exposure to H2S remains poorly understood. We investigated evolutionarily independent lineages of livebearing fishes that have colonized and adapted to springs rich in H2S and compared their genome-wide gene expression patterns with closely related lineages from adjacent, nonsulfidic streams. Significant differences in gene expression were uncovered between all sulfidic and nonsulfidic population pairs. Variation in the number of differentially expressed genes among population pairs corresponded to differences in divergence times and rates of gene flow, which is consistent with neutral drift driving a substantial portion of gene expression variation among populations. Accordingly, there was little evidence for convergent evolution shaping large-scale gene expression patterns among independent sulfide spring populations. Nonetheless, we identified a small number of genes that was consistently differentially expressed in the same direction in all sulfidic and nonsulfidic population pairs. Functional annotation of shared differentially expressed genes indicated upregulation of genes associated with enzymatic H2S detoxification and transport of oxidized sulfur species, oxidative phosphorylation, energy metabolism, and pathways involved in responses to oxidative stress. Overall, our results suggest that modification of processes associated with H2S detoxification and toxicity likely complement each other to mediate elevated H2S tolerance in sulfide spring fishes. Our analyses allow for the development of novel hypotheses about biochemical and physiological mechanisms of adaptation to extreme environments.

Transcriptome profiling of immune tissues reveals habitat-specific gene expression between lake and river sticklebacks

The observation of habitat-specific phenotypes suggests the action of natural selection. The three-spined stickleback (Gasterosteus aculeatus) has repeatedly colonized and adapted to diverse freshwater habitats across the northern hemisphere since the last glaciation, while giving rise to recurring phenotypes associated with specific habitats. Parapatric lake and river populations of sticklebacks harbour distinct parasite communities, a factor proposed to contribute to adaptive differentiation between these ecotypes. However, little is known about the transcriptional response to the distinct parasite pressure of those fish in a natural setting. Here, we sampled wild-caught sticklebacks across four geographical locations from lake and river habitats differing in their parasite load. We compared gene expression profiles between lake and river populations using 77 whole-transcriptome libraries from two immune-relevant tissues, the head kidney and the spleen. Differential expression analyses revealed 139 genes with habitat-specific expression patterns across the sampled population pairs. Among the 139 differentially expressed genes, eight are annotated with an immune function and 42 have been identified as differentially expressed in previous experimental studies in which fish have been immune challenged. Together, these findings reinforce the hypothesis that parasites contribute to adaptation of sticklebacks in lake and river habitats.

Functional Genomic Insights into Regulatory Mechanisms of High-Altitude Adaptation

Recent studies of indigenous human populations at high altitude have provided proof-of-principle that genome scans of DNA polymorphism can be used to identify candidate loci for hypoxia adaptation. When integrated with experimental analyses of physiological phenotypes, genome-wide surveys of DNA polymorphism and tissue-specific transcriptional profiles can provide insights into actual mechanisms of adaptation. It has been suggested that adaptive phenotypic evolution is largely mediated by cis-regulatory changes in genes that are located at integrative control points in regulatory networks. This hypothesis can be tested by conducting transcriptomic analyses of hypoxic signaling pathways in conjunction with experimental measures of vascular oxygen supply and metabolic pathway flux. Such studies may reveal whether the architecture of gene regulatory networks can be used to predict which loci (and which types of loci) are likely to be "hot spots" for adaptive physiological evolution. Functional genomic studies of deer mice (Peromyscus maniculatus) demonstrate how the integrated analysis of variation in tissue-specific transcriptomes, whole-animal physiological performance, and various subordinate traits can yield insights into the mechanistic underpinnings of high-altitude adaptation.

Developmental Transcriptome for a Facultatively Eusocial Bee, Megalopta genalis

Transcriptomes provide excellent foundational resources for mechanistic and evolutionary analyses of complex traits. We present a developmental transcriptome for the facultatively eusocial bee Megalopta genalis, which represents a potential transition point in the evolution of eusociality. A de novo transcriptome assembly of Megalopta genalis was generated using paired-end Illumina sequencing and the Trinity assembler. Males and females of all life stages were aligned to this transcriptome for analysis of gene expression profiles throughout development. Gene Ontology analysis indicates that stage-specific genes are involved in ion transport, cell-cell signaling, and metabolism. A number of distinct biological processes are upregulated in each life stage, and transitions between life stages involve shifts in dominant functional processes, including shifts from transcriptional regulation in embryos to metabolism in larvae, and increased lipid metabolism in adults. We expect that this transcriptome will provide a useful resource for future analyses to better understand the molecular basis of the evolution of eusociality and, more generally, phenotypic plasticity.

Adaptation and acclimation of aerobic exercise physiology in Lake Whitefish ecotypes (Coregonus clupeaformis)

The physiological mechanisms underlying local adaptation in natural populations of animals, and whether the same mechanisms contribute to adaptation and acclimation, are largely unknown. Therefore, we tested for evolutionary divergence in aerobic exercise physiology in laboratory bred, size-matched crosses of ancestral, benthic, normal Lake Whitefish (Coregonus clupeaformis) and derived, limnetic, more actively swimming "dwarf" ecotypes. We acclimated fish to constant swimming (emulating limnetic foraging) and control conditions (emulating normal activity levels) to simultaneously study phenotypic plasticity. We found extensive divergence between ecotypes: dwarf fish generally had constitutively higher values of traits related to oxygen transport (ventricle size) and use by skeletal muscle (percent oxidative muscle, mitochondrial content), and also evolved differential plasticity of mitochondrial function (Complex I activity and flux through Complexes I-IV and IV). The effects of swim training were less pronounced than differences among ecotypes and the traits which had a significant training effect (ventricle protein content, ventricle malate dehydrogenase activity, and muscle Complex V activity) did not differ among ecotypes. Only one trait, ventricle mass, varied in a similar manner with acclimation and adaptation and followed a pattern consistent with genetic accommodation. Overall, the physiological and biochemical mechanisms underlying acclimation and adaptation to swimming activity in Lake Whitefish differ.

Genetic approaches in comparative and evolutionary physiology

Whole animal physiological performance is highly polygenic and highly plastic, and the same is generally true for the many subordinate traits that underlie performance capacities. Quantitative genetics, therefore, provides an appropriate framework for the analysis of physiological phenotypes and can be used to infer the microevolutionary processes that have shaped patterns of trait variation within and among species. In cases where specific genes are known to contribute to variation in physiological traits, analyses of intraspecific polymorphism and interspecific divergence can reveal molecular mechanisms of functional evolution and can provide insights into the possible adaptive significance of observed sequence changes. In this review, we explain how the tools and theory of quantitative genetics, population genetics, and molecular evolution can inform our understanding of mechanism and process in physiological evolution. For example, lab-based studies of polygenic inheritance can be integrated with field-based studies of trait variation and survivorship to measure selection in the wild, thereby providing direct insights into the adaptive significance of physiological variation. Analyses of quantitative genetic variation in selection experiments can be used to probe interrelationships among traits and the genetic basis of physiological trade-offs and constraints. We review approaches for characterizing the genetic architecture of physiological traits, including linkage mapping and association mapping, and systems approaches for dissecting intermediary steps in the chain of causation between genotype and phenotype. We also discuss the promise and limitations of population genomic approaches for inferring adaptation at specific loci. We end by highlighting the role of organismal physiology in the functional synthesis of evolutionary biology.

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.

Ocean acidification research in the 'post-genomic' era: Roadmaps from the purple sea urchin Strongylocentrotus purpuratus

Advances in nucleic acid sequencing technology are removing obstacles that historically prevented use of genomics within ocean change biology. As one of the first marine calcifiers to have its genome sequenced, purple sea urchins (Strongylocentrotus purpuratus) have been the subject of early research exploring genomic responses to ocean acidification, work that points to future experiments and illustrates the value of expanding genomic resources to other marine organisms in this new 'post-genomic' era. This review presents case studies of S. purpuratus demonstrating the ability of genomic experiments to address major knowledge gaps within ocean acidification. Ocean acidification research has focused largely on species vulnerability, and studies exploring mechanistic bases of tolerance toward low pH seawater are comparatively few. Transcriptomic responses to high pCO₂ seawater in a population of urchins already encountering low pH conditions have cast light on traits required for success in future oceans. Secondly, there is relatively little information on whether marine organisms possess the capacity to adapt to oceans progressively decreasing in pH. Genomics offers powerful methods to investigate evolutionary responses to ocean acidification and recent work in S. purpuratus has identified genes under selection in acidified seawater. Finally, relatively few ocean acidification experiments investigate how shifts in seawater pH combine with other environmental factors to influence organism performance. In S. purpuratus, transcriptomics has provided insight into physiological responses of urchins exposed simultaneously to warmer and more acidic seawater. Collectively, these data support that similar breakthroughs will occur as genomic resources are developed for other marine species.

The environmental genomics of metazoan thermal adaptation

Continued and accelerating change in the thermal environment places an ever-greater priority on understanding how organisms are going to respond. The paradigm of ‘move, adapt or die', regarding ways in which organisms can respond to environmental stressors, stimulates intense efforts to predict the future of biodiversity. Assuming that extinction is an unpalatable outcome, researchers have focussed attention on how organisms can shift in their distribution to stay in the same thermal conditions or can stay in the same place by adapting to a changing thermal environment. How likely these respective outcomes might be depends on the answer to a fundamental evolutionary question, namely what genetic changes underpin adaptation to the thermal environment. The increasing access to and decreasing costs of next-generation sequencing (NGS) technologies, which can be applied to both model and non-model systems, provide a much-needed tool for understanding thermal adaptation. Here we consider broadly what is already known from non-NGS studies about thermal adaptation, then discuss the benefits and challenges of different NGS methodologies to add to this knowledge base. We then review published NGS genomics and transcriptomics studies of thermal adaptation to heat stress in metazoans and compare these results with previous non-NGS patterns. We conclude by summarising emerging patterns of genetic response and discussing future directions using these increasingly common techniques.

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.