What is needed for next-generation ecological and evolutionary genomics?

The Genomics Revolution in Nonmodel Species: Predictions vs. Reality for Salmonids

The increasing feasibility of whole-genome sequencing has been highly anticipated, promising to transform our understanding of the biology of nonmodel species. Notably, dramatic cost reductions beginning around 2007 with the advent of high-throughput sequencing inspired publications heralding the 'genomics revolution', with predictions about its future impacts. Although such predictions served as useful guideposts, value is added when statements are evaluated with the benefit of hindsight. Here, we review 10 key predictions made early in the genomics revolution, highlighting those realised while identifying challenges limiting others. We focus on predictions concerning applied aspects of genomics and examples involving salmonid species which, due to their socioeconomic and ecological significance, have been frontrunners in applications of genomics in nonmodel species. Predicted outcomes included enhanced analytical power, deeper insights into the genetic basis of phenotype and fitness variation, disease management and breeding program advancements. Although many predictions have materialised, several expectations remain unmet due to technological, analytical and knowledge barriers. Additionally, largely unforeseen advancements, including the identification and management applicability of large-effect loci, close-kin mark-recapture, environmental DNA and gene editing have added under-anticipated value. Finally, emerging innovations in artificial intelligence and bioinformatics offer promising new directions. This retrospective evaluation of the impacts of the genomic revolution offers insights into the future of genomics for nonmodel species.

Comparative transcriptomics reveals potential regulators of climate adaptation in Santalum album L. (Indian Sandalwood)

Santalum album L. (Indian Sandalwood), a valued tree species known for its fragrant heartwood and essential oil is facing increasing threat due to severe anthropogenic pressures compounded by climate change which has resulted in depletion of its adaptive gene pool. The present study investigates the transcriptome-level responses of nine sandalwood genotypes sourced from diverse climatic zones to identify adaptive genes in the species. Comparative transcriptomics predicted 727, 1141 and 479 differentially expressed transcripts (DETs) across wet vs. dry; monsoon vs. dry and wet vs. monsoon conditions, respectively, and majority of DETs were up-regulated in samples sourced from high rainfall areas. Transcripts including heat shock proteins, Zinc finger binding protein, ribosomal proteins, transcription factors and protein kinase were identified as probable regulators of climate adaptation in S. album. The expression changes of eight selected transcripts were further validated by real-time quantitative PCR. Protein-protein interaction analysis revealed key hub transcripts involved in climate response, while alternative splicing events in transcripts such as SURP and G-patch domain-containing protein 1-like protein, G-type lectin S-receptor-like serine/threonine protein kinase B120, Tetraspanin-3 and ARM repeat superfamily protein indicated the probable role of alternate splicing in increasing the transcript diversity during adaptation. This study presents the first insight into the molecular mechanisms of climate adaptation in the species and can form the basis for specific interventions such as selective breeding, genetic manipulation, and habitat management for conservation and long-term survival of sandalwood.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13205-025-04218-4.

Genomic Basis of Freshwater Adaptation in the Palaemonid Prawn Genus Macrobrachium: Convergent Evolution Following Multiple Independent Colonization Events

Adaptation to different salinity environments can enhance morphological and genomic divergence between related aquatic taxa. Species of prawns in the genus Macrobrachium naturally inhabit different osmotic niches and possess distinctive lifecycle traits associated with salinity tolerance. This study was conducted to investigate the patterns of adaptive genomic divergence during freshwater colonization in 34 Macrobrachium species collected from four continents; Australia, Asia, North and South America. Genotyping-by-sequencing (GBS) technique identified 5018 loci containing 82,636 single nucleotide polymorphisms (SNPs) that were used to reconstruct a phylogenomic tree. An additional phylogeny was reconstructed based on 43 candidate genes, previously identified as being potentially associated with freshwater adaptation. Comparison of the two phylogenetic trees revealed contrasting topologies. The GBS tree indicated multiple independent continent-specific invasions into freshwater by Macrobrachium lineages following common marine ancestry, as species with abbreviated larval development (ALD), i.e., species having a full freshwater life history, appeared reciprocally monophyletic within each continent. In contrast, the candidate gene tree showed convergent evolution for all ALD species worldwide, forming a single, well-supported clade. This latter pattern is likely the result of common evolutionary pressures selecting key mutations favored in continental freshwater habitats Results suggest that following multiple independent invasions into continental freshwaters at different evolutionary timescales, Macrobrachium taxa experienced adaptive genomic divergence, and in particular, convergence in the same genomic regions with parallel shifts in specific conserved phenotypic traits, such as evolution of larger eggs with abbreviated larval developmental.

Characterization, Selection, and Trans-Species Polymorphism in the MHC Class II of Heermann’s Gull (Charadriiformes)

The major histocompatibility complex (MHC) enables vertebrates to cope with pathogens and maintain healthy populations, thus making it a unique set of loci for addressing ecology and evolutionary biology questions. The aim of our study was to examine the variability of Heermann’s Gull MHC class II (MHCIIB) and compare these loci with other Charadriiformes. Fifty-nine MHCIIB haplotypes were recovered from sixty-eight Heermann’s Gulls by cloning, of them, twelve were identified as putative true alleles, forty-five as unique alleles, and two as pseudogenes. Intra and interspecific relationships indicated at least two loci in Heermann’s Gull MHCIIB and trans-species polymorphism among Charadriiformes (coinciding with the documented evidence of two ancient avian MHCIIB lineages, except in the Charadriidae family). Additionally, sites under diversifying selection revealed a better match with peptide-binding sites inferred in birds than those described in humans. Despite the negative anthropogenic activity reported on Isla Rasa, Heermann’s Gull showed MHCIIB variability consistent with population expansion, possibly due to a sudden growth following conservation efforts. Duplication must play an essential role in shaping Charadriiformes MHCIIB variability, buffering selective pressures through balancing selection. These findings suggest that MHC copy number and protected islands can contribute to seabird conservation.

Riverscape Genomics Clarifies Neutral and Adaptive Evolution in an Amazonian Characin Fish (Triportheus albus)

Understanding the role of natural selection in the evolution of wild populations is challenging due to the spatial complexity of natural systems. The richest diversity of freshwater fishes in the world is found in the Amazon Basin, a system where marked hydrochemical differences exist at the interface of major rivers with distinct “water colors” (i.e., black, white, and clear water). We hypothesize that divergent natural selection associated with these “aquatic ecotones” influences population-level adaptive divergence in the non-migratory Amazonian fish fauna. This hypothesis was tested using a landscape genomics framework to compare the relative contribution of environmental and spatial factors to the evolutionary divergence of the Amazonian characin fish Triportheus albus. The framework was based on spatial data, in situ hydrochemical measurements, and 15,251 filtered SNPs (single nucleotide polymorphisms) for T. albus sampled from three major Amazonian rivers. Gradient Forest, redundancy analysis (RDA) and BayPass analyses were used to test for signals of natural selection, and model-based and model-free approaches were used to evaluate neutral population differentiation. After controlling for a signal of neutral hierarchical structure which was consistent with the expectations for a dendritic system, variation in turbidity and pH were key factors contributing to adaptive divergence. Variation in genes involved in acid-sensitive ion transport pathways and light-sensitive photoreceptor pathways was strongly associated with pH and turbidity variability. This study improves our understanding of how natural selection and neutral evolution impact on the distribution of aquatic biodiversity from the understudied and ecologically complex Amazonia.

Genetic divergence and local adaptation of Liriodendron driven by heterogeneous environments

Ecological adaptive differentiation alters both the species diversity and intraspecific genetic diversity in forests, thus affecting the stability of forest ecosystems. Therefore, knowledge of the genetic underpinnings of the ecological adaptive differentiation of forest species is critical for effective species conservation. In this study, single-nucleotide polymorphisms (SNPs) from population transcriptomes were used to investigate the spatial distribution of genetic variation in Liriodendron to assess whether environmental variables can explain genetic divergence. We examined the contributions of environmental variables to population divergence and explored the genetic underpinnings of local adaptation using a landscape genomic approach. Niche models and statistical analyses showed significant niche divergence between L. chinense and L. tulipifera, suggesting that ecological adaptation may play a crucial role in driving interspecific divergence. We detected a new fine-scale genetic structure in L. chinense, and divergence of the six groups occurred during the late Pliocene to early Pleistocene. Redundancy analysis (RDA) revealed significant associations between genetic variation and multiple environmental variables. Environmental association analyses identified 67 environmental association loci (EALs; nonsynonymous SNPs) that underwent interspecific or intraspecific differentiation, 28 of which were associated with adaptive genes. These 28 candidate adaptive loci provide substantial evidence for local adaptation in Liriodendron. Our findings reveal ecological adaptive divergence pattern between Liriodendron species and provide novel insight into the role of heterogeneous environments in shaping genetic structure and driving local adaptation among populations, informing future L. chinense conservation efforts.

The parasite Schistocephalus solidus secretes proteins with putative host manipulation functions

BACKGROUND

Manipulative parasites are thought to liberate molecules in their external environment, acting as manipulation factors with biological functions implicated in their host's physiological and behavioural alterations. These manipulation factors are part of a complex mixture called the secretome. While the secretomes of various parasites have been described, there is very little data for a putative manipulative parasite. It is necessary to study the molecular interaction between a manipulative parasite and its host to better understand how such alterations evolve.

METHODS

Here, we used proteomics to characterize the secretome of a model cestode with a complex life cycle based on trophic transmission. We studied Schistocephalus solidus during the life stage in which behavioural changes take place in its obligatory intermediate fish host, the threespine stickleback (Gasterosteus aculeatus). We produced a novel genome sequence and assembly of S. solidus to improve protein coding gene prediction and annotation for this parasite. We then described the whole worm's proteome and its secretome during fish host infection using LC-MS/MS.

RESULTS

A total of 2290 proteins were detected in the proteome of S. solidus, and 30 additional proteins were detected specifically in the secretome. We found that the secretome contains proteases, proteins with neural and immune functions, as well as proteins involved in cell communication. We detected receptor-type tyrosine-protein phosphatases, which were reported in other parasitic systems to be manipulation factors. We also detected 12 S. solidus-specific proteins in the secretome that may play important roles in host-parasite interactions.

CONCLUSIONS

Our results suggest that S. solidus liberates molecules with putative host manipulation functions in the host and that many of them are species-specific.

Comparative transcriptomics of spotted seatrout (Cynoscion nebulosus) populations to cold and heat stress

Resilience to climate change depends on a species' adaptive potential and phenotypic plasticity. The latter can enhance survival of individual organisms during short periods of extreme environmental perturbations, allowing genetic adaptation to take place over generations. Along the U.S. East Coast, estuarine-dependent spotted seatrout (Cynoscion nebulosus) populations span a steep temperature gradient that provides an ideal opportunity to explore the molecular basis of phenotypic plasticity. Genetically distinct spotted seatrout sampled from a northern and a southern population were exposed to acute cold and heat stress (5 biological replicates in each treatment and control group), and their transcriptomic responses were compared using RNA-sequencing (RNA-seq). The southern population showed a larger transcriptomic response to acute cold stress, whereas the northern population showed a larger transcriptomic response to acute heat stress compared with their respective population controls. Shared transcripts showing significant differences in expression levels were predominantly enriched in pathways that included metabolism, transcriptional regulation, and immune response. In response to heat stress, only the northern population significantly upregulated genes in the apoptosis pathway, which could suggest greater vulnerability to future heat waves in this population as compared to the southern population. Genes showing population-specific patterns of expression, including hpt, acot, hspa5, and hsc71, are candidates for future studies aiming to monitor intraspecific differences in temperature stress responses in spotted seatrout. Our findings contribute to the current understanding of phenotypic plasticity and provide a basis for predicting the response of a eurythermal fish species to future extreme temperatures.

The power and limitations of gene expression pathway analyses toward predicting population response to environmental stressors

Rapid environmental changes impact the global distribution and abundance of species, highlighting the urgency to understand and predict how populations will respond. The analysis of differentially expressed genes has elucidated areas of the genome involved in adaptive divergence to past and present environmental change. Such studies however have been hampered by large numbers of differentially expressed genes and limited knowledge of how these genes work in conjunction with each other. Recent methods (broadly termed "pathway analyses") have emerged that aim to group genes that behave in a coordinated fashion to a factor of interest. These methods aid in functional annotation and uncovering biological pathways, thereby collapsing complex datasets into more manageable units, providing more nuanced understandings of both the organism-level effects of modified gene expression, and the targets of adaptive divergence. Here, we reanalyze a dataset that investigated temperature-induced changes in gene expression in marine-adapted and freshwater-adapted threespine stickleback (Gasterosteus aculeatus), using Weighted Gene Co-expression Network Analysis (WGCNA) with PANTHER Gene Ontology (GO)-Slim overrepresentation and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis. Six modules exhibited a conserved response and six a divergent response between marine and freshwater stickleback when acclimated to 7°C or 22°C. One divergent module showed freshwater-specific response to temperature, and the remaining divergent modules showed differences in height of reaction norms. PPARAa, a transcription factor that regulates fatty acid metabolism and has been implicated in adaptive divergence, was located in a module that had higher expression at 7°C and in freshwater stickleback. This updated methodology revealed patterns that were not found in the original publication. Although such methods hold promise toward predicting population response to environmental stressors, many limitations remain, particularly with regard to module expression representation, database resources, and cross-database integration.

Evaluating Population Genomic Candidate Genes Underlying Flowering Time in Arabidopsis thaliana Using T-DNA Insertion Lines

Population genomic scans have emerged as a powerful tool to detect regions of the genome that are potential targets of selection. Despite the success of genomic scans in identifying novel lists of loci potentially underlying adaptation, few studies proceed to validate the function of these candidate genes. In this study, we used transfer-DNA (T-DNA) insertion lines to evaluate the effects of 27 candidate genes on flowering time in North American accessions of Arabidopsis thaliana. We compared the flowering time of T-DNA insertion lines that knock out the function of a candidate gene obtained from population genomic studies to a wild type under long- and short-day conditions. We also did the same for a collection of randomly chosen genes that had not been identified as candidates. We validated the well-known effect of long-day conditions in accelerating flowering time and found that gene disruption caused by insertional mutagenesis tends to delay flowering. Surprisingly, we found that knockouts in random genes were just as likely to produce significant phenotypic effects as knockouts in candidate genes. T-DNA insertions at a handful of candidate genes that had previously been identified as outlier loci showed significant delays in flowering time under both long and short days, suggesting that they are promising candidates for future investigation.

Comparative transcriptomics of an arctic foundation species, tussock cottongrass (Eriophorum vaginatum), during an extreme heat event

Tussock cottongrass (Eriophorum vaginatum) is a foundation species for much of the arctic moist acidic tundra, which is currently experiencing extreme effects of climate change. The Arctic is facing higher summer temperatures and extreme weather events are becoming more common. We used Illumina RNA-Seq to analyse cDNA libraries for differential expression of genes from leaves of ecologically well-characterized ecotypes of tussock cottongrass found along a latitudinal gradient in the Alaskan Arctic and transplanted into a common garden. Plant sampling was performed on a typical summer day and during an extreme heat event. We obtained a de novo assembly that contained 423,353 unigenes. There were 363 unigenes up-regulated and 1,117 down-regulated among all ecotypes examined during the extreme heat event. Of these, 26 HSP unigenes had >log2-fold up-regulation. Several TFs associated with heat stress in previous studies were identified that had >log2-fold up- or down-regulation during the extreme heat event (e.g., DREB, NAC). There was consistent variation in DEGs among ecotypes, but not specifically related to whether plants originated from taiga or tundra ecosystems. As the climate changes it is essential to determine ecotypic diversity at the genomic level, especially for widespread species that impact ecosystem function.

A comparative synthesis of transcriptomic analyses reveals major differences between WSSV-susceptible Litopenaeus vannamei and WSSV-refractory Macrobrachium rosenbergii

Since the 1990s White Spot Syndrome Virus (WSSV) has severely affected shrimp aquaculture worldwide causing a global pandemic of White Spot Disease (WSD) in penaeid culture. However, not all decapod species that can be infected by WSSV show the same susceptibility to the virus, thus raising interesting questions regarding the potential genetic traits that might confer resistance to WSSV. In order to shed light into the genetic markers of WSSV resistance, we employed a dual approach: i) we initially analysed the transcriptomes derived from the hepatopancreas of two species, the susceptible white shrimp Litopenaeus vannamei and the refractory fresh water prawn Macrobrachium rosenbergii, both infected with WSSV. We found a large number of differentially expressed genes (DEGs) belonging to the immune system (mostly anti-microbial peptides and haemolymph clotting components) that were generally up-regulated in M. rosenbergii and down-regulated in L. vannamei. Further, in both species we identified many up-regulated DEGs that were related to metabolism (suggesting a metabolic shift during the infection) and, interestingly, in L. vannamei only, we found several DEGs that were related to moult and suggested an inhibition of the moult cycle in this species following WSSV infection. ii) we then identified a limited number of genetic markers putatively linked with WSD tolerance by employing an ecological genomics approach in which we compared published reports with our own RNA-seq datasets for different decapod species infected with WSSV. Using this second comparative approach, we found nine candidate genes which are consistently down-regulated in susceptible species and up-regulated in refractory species and which have a role in immune response. Together our data offer novel insights into gene expression differences that can be found in susceptible and refractory decapod species infected with WSSV and provide a valuable resource towards our understanding of the potential genetic basis of tolerance to WSSV.

Gene expression response to a nematode parasite in novel and native eel hosts

Invasive parasites are involved in population declines of new host species worldwide. The high susceptibilities observed in many novel hosts have been attributed to the lack of protective immunity to the parasites which native hosts acquired during their shared evolution. We experimentally infected Japanese eels (Anguilla japonica) and European eels (Anguilla anguilla) with Anguillicola crassus, a nematode parasite that is native to the Japanese eel and invasive in the European eel. We inferred gene expression changes in head kidney tissue from both species, using RNA-seq data to determine the responses at two time points during the early stages of infection (3 and 23 days postinfection). At both time points, the novel host modified the expression of a larger and functionally more diverse set of genes than the native host. Strikingly, the native host regulated immune gene expression only at the earlier time point and to a small extent while the novel host regulated these genes at both time points. A low number of differentially expressed immune genes, especially in the native host, suggest that a systemic immune response was of minor importance during the early stages of infection. Transcript abundance of genes involved in cell respiration was reduced in the novel host which may affect its ability to cope with harsh conditions and energetically demanding activities. The observed gene expression changes in response to a novel parasite that we observed in a fish follow a general pattern observed in amphibians and mammals, and suggest that the disruption of physiological processes, rather than the absence of an immediate immune response, is responsible for the higher susceptibility of the novel host.

Effects of marker type and filtering criteria on Q ST-F ST comparisons

Comparative studies of quantitative and neutral genetic differentiation (Q ST-F ST tests) provide means to detect adaptive population differentiation. However, Q ST-F ST tests can be overly liberal if the markers used deflate F ST below its expectation, or overly conservative if methodological biases lead to inflated F ST estimates. We investigated how marker type and filtering criteria for marker selection influence Q ST-F ST comparisons through their effects on F ST using simulations and empirical data on over 18 000 in silico genotyped microsatellites and 3.8 million single-locus polymorphism (SNP) loci from four populations of nine-spined sticklebacks (Pungitius pungitius). Empirical and simulated data revealed that F ST decreased with increasing marker variability, and was generally higher with SNPs than with microsatellites. The estimated baseline F ST levels were also sensitive to filtering criteria for SNPs: both minor alleles and linkage disequilibrium (LD) pruning influenced F ST estimation, as did marker ascertainment. However, in the case of stickleback data used here where Q ST is high, the choice of marker type, their genomic location, ascertainment and filtering made little difference to outcomes of Q ST-F ST tests. Nevertheless, we recommend that Q ST-F ST tests using microsatellites should discard the most variable loci, and those using SNPs should pay attention to marker ascertainment and properly account for LD before filtering SNPs. This may be especially important when level of quantitative trait differentiation is low and levels of neutral differentiation high.

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.

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.

Opportunities and limitations of reduced representation bisulfite sequencing in plant ecological epigenomics

Contents Summary 738 I. Introduction 738 II. RRBS loci as genome-wide epigenetic markers 739 III. Exploiting functional annotation of RRBS loci 739 IV. Limitations of RRBS methods for nonmodel species 740 V. Maximising the impact of RRBS in plants 741 VI. Conclusions 741 Acknowledgements 741 SUMMARY: Investigating the features and implications of epigenetic mechanisms across the breadth of organisms and ecosystems is important for understanding the ecological relevance of epigenetics. Several cost-effective reduced representation bisulfite sequencing approaches (RRBS) have been recently developed and applied to different organisms that lack a well annotated reference genome. These new approaches improve the assessment of epigenetic diversity in ecological settings and may provide functional insights. We assess here the opportunities and limitations of RRBS in nonmodel plant species. Well thought out experimental designs that include complementary gene expression studies, and the improvement of genomics resources for the target group, promise to maximize the effect of future RRBS studies.

R(NA)-tistic expression: The art of matching unknown mRNA and proteins to environmental response in ecological genomics

A challenge of modern ecological genomics is reducing uncertainty surrounding the biological inferences from gene expression. For example, approximately 40% of proteins in eukaryotic model organisms do not contain characterized domains (Gollery et al., 2006). Even proteins of "known function" are typically only characterized in the sense that they have a domain function, but provide no information on their biological role within the cell (e.g., activation, pathways or targets). Yet, as molecular ecologists, a common objective is to elucidate how organisms respond to environmental variation through changes in gene expression, including homoeostatic, acclimatory, and adaptive responses to environmental stressors, a challenge increased by poor protein ecological annotation. Now, in this issue of Molecular Ecology, Orsini et al. (2017) use the quintessential Daphnia system to characterize the differences in stress response in three genotypic backgrounds to common biotic and abiotic stressors found in nature. Using an optimized weighted gene co-expression network analysis, they link genes of unknown function to genes that they co-activate with and enrich for gene ontology. Determining the functional networks of genes that behave in genotype- and treatment-specific responses gives insight into possible pathways and respective ecological roles, helping pave the way for the next generation of transcriptomic studies in molecular ecology.

Temperature gradient affects differentiation of gene expression and SNP allele frequencies in the dominant Lake Baikal zooplankton species

Local adaptation and phenotypic plasticity are main mechanisms of organisms' resilience in changing environments. Both are affected by gene flow and are expected to be weak in zooplankton populations inhabiting large continuous water bodies and strongly affected by currents. Lake Baikal, the deepest and one of the coldest lakes on Earth, experienced epilimnion temperature increase during the last 100 years, exposing Baikal's zooplankton to novel selective pressures. We obtained a partial transcriptome of Epischura baikalensis (Copepoda: Calanoida), the dominant component of Baikal's zooplankton, and estimated SNP allele frequencies and transcript abundances in samples from regions of Baikal that differ in multiyear average surface temperatures. The strongest signal in both SNP and transcript abundance differentiation is the SW-NE gradient along the 600+ km long axis of the lake, suggesting isolation by distance. SNP differentiation is stronger for nonsynonymous than synonymous SNPs and is paralleled by differential survival during a laboratory exposure to increased temperature, indicating directional selection operating on the temperature gradient. Transcript abundance, generally collinear with the SNP differentiation, shows samples from the warmest, less deep location clustering together with the southernmost samples. Differential expression is more frequent among transcripts orthologous to candidate thermal response genes previously identified in model arthropods, including genes encoding cytoskeleton proteins, heat-shock proteins, proteases, enzymes of central energy metabolism, lipid and antioxidant pathways. We conclude that the pivotal endemic zooplankton species in Lake Baikal exists under temperature-mediated selection and possesses both genetic variation and plasticity to respond to novel temperature-related environmental pressures.

Comparative transcriptomics of cyprinid minnows and carp in a common wild setting: a resource for ecological genomics in freshwater communities

Comparative transcriptomics can now be conducted on organisms in natural settings, which has greatly enhanced understanding of genome–environment interactions. Here, we demonstrate the utility and potential pitfalls of comparative transcriptomics of wild organisms, with an example from three cyprinid fish species (Teleostei:Cypriniformes). We present extensively filtered and annotated transcriptome assemblies that provide a valuable resource for studies of genome evolution (e.g. polyploidy), ecological and morphological diversification, speciation, and shared and unique responses to environmental variation in cyprinid fishes. Our results and analyses address the following points: (i) ‘essential developmental genes’ are shown to be ubiquitously expressed in a diverse suite of tissues across later ontogenetic stages (i.e. juveniles and adults), making these genes are useful for assessing the quality of transcriptome assemblies, (ii) the influence of microbiomes and other exogenous DNA, (iii) potentially novel, species-specific genes, and (iv) genomic rearrangements (e.g. whole genome duplication). The data we present provide a resource for future comparative work in cypriniform fishes and other taxa across a variety of sub-disciplines, including stress response, morphological diversification, community ecology, ecotoxicology, and climate change.

The utility of transcriptomics in fish conservation

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

The Rapid Evolution of an Ohnolog Contributes to the Ecological Specialization of Incipient Yeast Species

Identifying the molecular changes that lead to ecological specialization during speciation is one of the major goals of molecular evolution. One question that remains to be thoroughly investigated is whether ecological specialization derives strictly from adaptive changes and their associated trade-offs, or from conditionally neutral mutations that accumulate under relaxed selection. We used whole-genome sequencing, genome annotation and computational analyses to identify genes that have rapidly diverged between two incipient species of Saccharomyces paradoxus that occupy different climatic regions along a south-west to north-east gradient. As candidate loci for ecological specialization, we identified genes that show signatures of adaptation and accelerated rates of amino acid substitutions, causing asymmetric evolution between lineages. This set of genes includes a glycyl-tRNA-synthetase, GRS2, which is known to be transcriptionally induced under heat stress in the model and sister species S. cerevisiae. Molecular modelling, expression analysis and fitness assays suggest that the accelerated evolution of this gene in the Northern lineage may be caused by relaxed selection. GRS2 arose during the whole-genome duplication (WGD) that occurred 100 million years ago in the yeast lineage. While its ohnolog GRS1 has been preserved in all post-WGD species, GRS2 has frequently been lost and is evolving rapidly, suggesting that the fate of this ohnolog is still to be resolved. Our results suggest that the asymmetric evolution of GRS2 between the two incipient S. paradoxus species contributes to their restricted climatic distributions and thus that ecological specialization derives at least partly from relaxed selection rather than a molecular trade-off resulting from adaptive evolution.

Ecological plant epigenetics: Evidence from model and non-model species, and the way forward

Growing evidence shows that epigenetic mechanisms contribute to complex traits, with implications across many fields of biology. In plant ecology, recent studies have attempted to merge ecological experiments with epigenetic analyses to elucidate the contribution of epigenetics to plant phenotypes, stress responses, adaptation to habitat, and range distributions. While there has been some progress in revealing the role of epigenetics in ecological processes, studies with non-model species have so far been limited to describing broad patterns based on anonymous markers of DNA methylation. In contrast, studies with model species have benefited from powerful genomic resources, which contribute to a more mechanistic understanding but have limited ecological realism. Understanding the significance of epigenetics for plant ecology requires increased transfer of knowledge and methods from model species research to genomes of evolutionarily divergent species, and examination of responses to complex natural environments at a more mechanistic level. This requires transforming genomics tools specifically for studying non-model species, which is challenging given the large and often polyploid genomes of plants. Collaboration among molecular geneticists, ecologists and bioinformaticians promises to enhance our understanding of the mutual links between genome function and ecological processes.

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

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

Conserved Transcription Factors Steer Growth-Related Genomic Programs in Daphnia

Ecological genomics aims to understand the functional association between environmental gradients and the genes underlying adaptive traits. Many genes that are identified by genome-wide screening in ecologically relevant species lack functional annotations. Although gene functions can be inferred from sequence homology, such approaches have limited power. Here, we introduce ecological regulatory genomics by presenting an ontology-free gene prioritization method. Specifically, our method combines transcriptome profiling with high-throughput cis-regulatory sequence analysis in the water fleas Daphnia pulex and Daphnia magna. It screens coexpressed genes for overrepresented DNA motifs that serve as transcription factor binding sites, thereby providing insight into conserved transcription factors and gene regulatory networks shaping the expression profile. We first validated our method, called Daphnia-cisTarget, on a D. pulex heat shock data set, which revealed a network driven by the heat shock factor. Next, we performed RNA-Seq in D. magna exposed to the cyanobacterium Microcystis aeruginosa. Daphnia-cisTarget identified coregulated gene networks that associate with the moulting cycle and potentially regulate life history changes in growth rate and age at maturity. These networks are predicted to be regulated by evolutionary conserved transcription factors such as the homologues of Drosophila Shavenbaby and Grainyhead, nuclear receptors, and a GATA family member. In conclusion, our approach allows prioritising candidate genes in Daphnia without bias towards prior knowledge about functional gene annotation and represents an important step towards exploring the molecular mechanisms of ecological responses in organisms with poorly annotated genomes.

Major host transitions are modulated through transcriptome-wide reprogramming events in Schistocephalus solidus, a threespine stickleback parasite

Parasites with complex life cycles have developed numerous phenotypic strategies, closely associated with developmental events, to enable the exploitation of different ecological niches and facilitate transmission between hosts. How these environmental shifts are regulated from a metabolic and physiological standpoint, however, still remain to be fully elucidated. We examined the transcriptomic response of Schistocephalus solidus, a trophically transmitted parasite with a complex life cycle, over the course of its development in an intermediate host, the threespine stickleback, and the final avian host. Results from our differential gene expression analysis show major reprogramming events among developmental stages. The final host stage is characterized by a strong activation of reproductive pathways and redox homoeostasis. The attainment of infectivity in the fish intermediate host-which precedes sexual maturation in the final host and is associated with host behaviour changes-is marked by transcription of genes involved in neural pathways and sensory perception. Our results suggest that un-annotated and S. solidus-specific genes could play a determinant role in host-parasite molecular interactions required to complete the parasite's life cycle. Our results permit future comparative analyses to help disentangle species-specific patterns of infection from conserved mechanisms, ultimately leading to a better understanding of the molecular control and evolution of complex life cycles.

Eigenvector metabolite analysis reveals dietary effects on the association among metabolite correlation patterns, gene expression, and phenotypes

INTRODUCTION

'Multi-omics' datasets obtained from an organism of interest reared under different environmental treatments are increasingly common. Identifying the links among metabolites and transcripts can help to elucidate our understanding of the impact of environment at different levels within the organism. However, many methods for characterizing physiological connections cannot address unidentified metabolites.

OBJECTIVES

Here, we use Eigenvector Metabolite Analysis (EvMA) to examine links between metabolomic, transcriptomic, and phenotypic variation data and to assess the impact of environmental factors on these associations. Unlike other methods, EvMA can be used to analyze datasets that include unidentified metabolites and unannotated transcripts.

METHODS

To demonstrate the utility of EvMA, we analyzed metabolomic, transcriptomic, and phenotypic datasets produced from 20 Drosophila melanogaster genotypes reared on four dietary treatments. We used a hierarchical distance-based method to cluster the metabolites. The links between metabolite clusters, gene expression, and overt phenotypes were characterized using the eigenmetabolite (first principal component) of each cluster.

RESULTS

EvMA recovered chemically related groups of metabolites within the clusters. Using the eigenmetabolite, we identified genes and phenotypes that significantly correlated with each cluster. EvMA identifies new connections between the phenotypes, metabolites, and gene transcripts. Conclusion EvMA provides a simple method to identify correlations between metabolites, gene expression, and phenotypes, which can allow us to partition multivariate datasets into meaningful biological modules and identify under-studied metabolites and unannotated gene transcripts that may be central to important biological processes. This can be used to inform our understanding of the effect of environmental mechanisms underlying physiological states of interest.

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.

Using an integrative approach to investigate the evolution of behaviour

Behaviour is a central focus of interest in biology because it has an impact on several aspects of an organism's life. Evolutionary biologists have realised the advantage of an integrative approach that jointly studies the molecular, cellular and physiological levels of an individual to link them with the organismal behavioural phenotype. First, this mechanistic information helps in understanding physiological and evolutionary constraints acting on the behavioural response to the environment and its evolution. Second, it furthers our understanding of the process of molecular convergent evolution. Finally, we learn about natural variation in molecular, cellular and physiological traits present in wild populations and their underlying genetic basis, which can be a substrate for selection to act on. I illustrate these points using our work on behaviour variation in fishes. The information on the mechanistic bases of behaviour variation in various species and behaviours will contribute to an ecological annotation of genes and to uncover new mechanisms implicated in how this astonishing behavioural diversity arose, is maintained and will evolve.

Connecting genotypes, phenotypes and fitness: harnessing the power of CRISPR/Cas9 genome editing

One of the fundamental goals in evolution and ecology is to identify the genetic basis of adaptive phenotypes. Unfortunately, progress towards this goal has been hampered by a lack of genetic tools available for nonmodel organisms. The exciting new development of the CRISPR (clustered regularly interspaced short palindromic repeat)/Cas9 (CRISPR-associated nuclease 9) genome-editing system now promises to transform the field of molecular ecology by providing a versatile toolkit for manipulating the genome of a wide variety of organisms. Here, we review the numerous applications of this groundbreaking technology and provide a practical guide to the creation of genetic knockouts, transgenics and other related forms of gene manipulation in nonmodel organisms. We also specifically discuss the potential uses of the CRISPR/Cas9 system in ecological and evolutionary studies, which will further advance the field towards the long-standing goal of connecting genotypes, phenotypes and fitness.

Towards the identification of the loci of adaptive evolution

1. Establishing the genetic and molecular basis underlying adaptive traits is one of the major goals of evolutionary geneticists in order to understand the connection between genotype and phenotype and elucidate the mechanisms of evolutionary change. Despite considerable effort to address this question, there remain relatively few systems in which the genes shaping adaptations have been identified. 2. Here, we review the experimental tools that have been applied to document the molecular basis underlying evolution in several natural systems, in order to highlight their benefits, limitations and suitability. In most cases, a combination of DNA, RNA and functional methodologies with field experiments will be needed to uncover the genes and mechanisms shaping adaptation in nature.

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.

Transcriptomic responses to emamectin benzoate in Pacific and Atlantic Canada salmon lice Lepeophtheirus salmonis with differing levels of drug resistance

Salmon lice Lepeophtheirus salmonis are an ecologically and economically important parasite of wild and farmed salmon. In Scotland, Norway, and Eastern Canada, L. salmonis have developed resistance to emamectin benzoate (EMB), one of the few parasiticides available for salmon lice. Drug resistance mechanisms can be complex, potentially differing among populations and involving multiple genes with additive effects (i.e., polygenic resistance). Indicators of resistance development may enable early detection and countermeasures to avoid the spread of resistance. Here, we collect sensitive Pacific L. salmonis and sensitive and resistant Atlantic L. salmonis from salmon farms, propagate in laboratory (F1), expose to EMB in bioassays, and evaluate either baseline (Atlantic only) or induced transcriptomic differences between populations. In all populations, induced responses were minor and a cellular stress response was not identified. Pacific lice did not upregulate any genes in response to EMB, but downregulated degradative enzymes and transport proteins at 50 ppb EMB. Baseline differences between sensitive and now resistant Atlantic lice were much greater than responses to exposures. All resistant lice overexpressed degradative enzymes, and resistant males, the most resistant group, overexpressed collagenases to the greatest extent. These results indicate an accumulation of baseline expression differences related to resistance.

Functional genomics of acclimation and adaptation in response to thermal stress in Daphnia

Background

Gene expression regulation is one of the fundamental mechanisms of phenotypic plasticity and is expected to respond to selection in conditions favoring phenotypic response. The observation that many organisms increase their stress tolerance after acclimation to moderate levels of stress is an example of plasticity which has been long hypothesized to be based on adaptive changes in gene expression. We report genome-wide patterns of gene expression in two heat-tolerant and two heat-sensitive parthenogenetic clones of the zooplankton crustacean Daphnia pulex exposed for three generations to either optimal (18°C) or substressful (28°C) temperature.

Results

A large number of genes responded to temperature and many demonstrated a significant genotype-by-environment (GxE) interaction. Among genes with a significant GxE there were approximately equally frequent instances of canalization, i.e. stronger plasticity in heat-sensitive than in heat-tolerant clones, and of enhancement of plasticity along the evolutionary vector toward heat tolerance. The strongest response observed is the across-the-board down-regulation of a variety of genes occurring in heat-tolerant, but not in heat-sensitive clones. This response is particularly obvious among genes involved in core metabolic pathways and those responsible for transcription, translation and DNA repair.

Conclusions

The observed down-regulation of metabolism, consistent with previous findings in yeast and Drosophila, may reflect a general compensatory stress response. The associated down-regulation of DNA repair pathways potentially creates a trade-off between short-term benefits of survival at high temperature and long-term costs of accelerated mutation accumulation.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-859) contains supplementary material, which is available to authorized users.

Parallel signatures of selection in temporally isolated lineages of pink salmon

Studying the effect of similar environments on diverse genetic backgrounds has long been a goal of evolutionary biologists with studies typically relying on experimental approaches. Pink salmon, a highly abundant and widely ranging salmonid, provide a naturally occurring opportunity to study the effects of similar environments on divergent genetic backgrounds due to a strict two-year semelparous life history. The species is composed of two reproductively isolated lineages with overlapping ranges that share the same spawning and rearing environments in alternate years. We used restriction-site-associated DNA (RAD) sequencing to discover and genotype approximately 8000 SNP loci in three population pairs of even- and odd-year pink salmon along a latitudinal gradient in North America. We found greater differentiation within the odd-year than within the even-year lineage and greater differentiation in the southern pair from Puget Sound than in the northern Alaskan population pairs. We identified 15 SNPs reflecting signatures of parallel selection using both a differentiation-based method (BAYESCAN) and an environmental correlation method (BAYENV). These SNPs represent genomic regions that may be particularly informative in understanding adaptive evolution in pink salmon and exploring how differing genetic backgrounds within a species respond to selection from the same natural environment.

Applying evolutionary concepts to wildlife disease ecology and management

Existing and emerging infectious diseases are among the most pressing global threats to biodiversity, food safety and human health. The complex interplay between host, pathogen and environment creates a challenge for conserving species, communities and ecosystem functions, while mediating the many known ecological and socio-economic negative effects of disease. Despite the clear ecological and evolutionary contexts of host-pathogen dynamics, approaches to managing wildlife disease remain predominantly reactionary, focusing on surveillance and some attempts at eradication. A few exceptional studies have heeded recent calls for better integration of ecological concepts in the study and management of wildlife disease; however, evolutionary concepts remain underused. Applied evolution consists of four principles: evolutionary history, genetic and phenotypic variation, selection and eco-evolutionary dynamics. In this article, we first update a classical framework for understanding wildlife disease to integrate better these principles. Within this framework, we explore the evolutionary implications of environment-disease interactions. Subsequently, we synthesize areas where applied evolution can be employed in wildlife disease management. Finally, we discuss some future directions and challenges. Here, we underscore that despite some evolutionary principles currently playing an important role in our understanding of disease in wild animals, considerable opportunities remain for fostering the practice of evolutionarily enlightened wildlife disease management.

Transcriptome resources for the white-footed mouse (Peromyscus leucopus): new genomic tools for investigating ecologically divergent urban and rural populations

Genomic resources are important and attainable for examining evolutionary change in divergent natural populations of nonmodel species. We utilized two next-generation sequencing (NGS) platforms, 454 and SOLiD 5500XL, to assemble low-coverage transcriptomes of the white-footed mouse (Peromyscus leucopus), a widespread and abundant native rodent in eastern North America. We sequenced liver mRNA transcripts from multiple individuals collected from urban populations in New York City and rural populations in undisturbed protected areas nearby and assembled a reference transcriptome using 1 080 065 954 SOLiD 5500XL (75 bp) reads and 3 052 640 454 GS FLX + reads. The reference contained 40 908 contigs with a N50 = 1044 bp and a total content of 30.06 Megabases (Mb). Contigs were annotated from Mus musculus (39.96% annotated) Uniprot databases. We identified 104 655 high-quality single nucleotide polymorphisms (SNPs) and 65 single sequence repeats (SSRs) with flanking primers. We also used normalized read counts to identify putative gene expression differences in 10 genes between populations. There were 19 contigs significantly differentially expressed in urban populations compared to rural populations, with gene function annotations generally related to the translation and modification of proteins and those involved in immune responses. The individual transcriptomes generated in this study will be used to investigate evolutionary responses to urbanization. The reference transcriptome provides a valuable resource for the scientific community using North American Peromyscus species as emerging model systems for ecological genetics and adaptation.

The shape of things to come in the study of the origin of species?

Perhaps Darwin would agree that speciation is no longer the mystery of mysteries that it used to be. It is now generally accepted that evolution by natural selection can contribute to ecological adaptation, resulting in the evolution of reproductive barriers and, hence, to the evolution of new species (Schluter & Conte 2009; Meyer 2011; Nosil 2012). From genes that encode silencing proteins that cause infertility in hybrid mice (Mihola et al. 2009), to segregation distorters linked to speciation in fruit flies (Phadnis & Orr 2009), or pollinator-mediated selection on flower colour alleles driving reinforcement in Texan wildflowers (Hopkins & Rausher 2012), characterization of the genes that drive speciation is providing clues to the origin of species (Nosil & Schluter 2011). It is becoming apparent that, while recent work continues to overturn historical ideas about sympatric speciation (e.g. Barluenga et al. 2006), ecological circumstances strongly influence patterns of genomic divergence, and ultimately the establishment of reproductive isolation when gene flow is present (Elmer & Meyer 2011). Less clear, however, are the genetic mechanisms that cause speciation, particularly when ongoing gene flow is occurring. Now, in this issue, Franchini et al. (2014) employ a classic genetic mapping approach augmented with new genomic tools to elucidate the genomic architecture of ecologically divergent body shapes in a pair of sympatric crater lake cichlid fishes. From over 450 segregating SNPs in an F2 cross, 72 SNPs were linked to 11 QTL associated with external morphology measured by means of traditional and geometric morphometrics. Annotation of two highly supported QTL further pointed to genes that might contribute to ecological divergence in body shape in Midas cichlids, overall supporting the hypothesis that genomic regions of large phenotypic effect may be contributing to early-stage divergence in Midas cichlids.

Divergent immunity and energetic programs in the gills of migratory and resident Oncorhynchus mykiss

Divergent life history strategies occur in steelhead or rainbow trout Oncorhynchus mykiss, and many populations produce both migrant (anadromous fish that move to the ocean after rearing) and resident (do not migrate and remain in fresh water) individuals. Mechanisms leading to each type are only partially understood; while the general tendency of a population is heritable, individual tendency may be plastic, influenced by local environment. Steelhead hatchery programmes aim to mitigate losses in wild stocks by producing trout that will migrate to the ocean and not compete with wild trout for limited freshwater resources. To increase our understanding of gill function in these migratory or resident phenotypes, here we compare gill transcriptome profiles of hatchery-released fish either at the release site (residents) or five river kilometres downstream while still in full fresh water (migrants). To test whether any of these genes can be used as predictive markers for smoltification, we compared these genes between migrant-like and undifferentiated trout while still in the hatchery in a common environment (prerelease). Results confirmed the gradual process of smoltification, and the importance of energetics, gill remodelling and ion transport capacity for migrants. Additionally, residents overexpressed transcripts involved in antiviral defences, potentially for immune surveillance via dendritic cells in the gills. The best smoltification marker candidate was protein s100a4, expression of which was highly correlated with Na(+) , K(+) ATPase (NKA) activity and smolt-like morphology in pre- and postrelease trout gills.

Challenges and prospects in genome-wide quantitative trait loci mapping of standing genetic variation in natural populations

A considerable challenge in evolutionary genetics is to understand the genetic mechanisms that facilitate or impede evolutionary adaptation in natural populations. For this, we must understand the genetic loci contributing to trait variation and the selective forces acting on them. The decreased costs and increased feasibility of obtaining genotypic data on a large number of individuals have greatly facilitated gene mapping in natural populations, particularly because organisms whose genetics have been historically difficult to study are now within reach. Here we review the methods available to evolutionary ecologists interested in dissecting the genetic basis of traits in natural populations. Our focus lies on standing genetic variation in outbred populations. We present an overview of the current state of research in the field, covering studies on both plants and animals. We also draw attention to particular challenges associated with the discovery of quantitative trait loci and discuss parallels to studies on crops, livestock, and humans. Finally, we point to some likely future developments in genetic mapping studies.

Integrating phenotypic plasticity within an Ecological Genomics framework: recent insights from the genomics, evolution, ecology, and fitness of plasticity

E.B. Ford's 1964 book Ecological Genetics was a call for biologists to engage in multidisciplinary work in order to elucidate the link between genotype, phenotype, and fitness for ecologically relevant traits. In this review, we argue that the integration of an ecological genomics framework in studies of phenotypic plasticity is a promising approach to elucidate the causal links between genes and the environment, particularly during colonization of novel environments, environmental change, and speciation. This review highlights some of the questions and hypotheses generated from a mechanistic, evolutionary, and ecological perspective, in order to direct the continued and future use of genomic tools in the study of phenotypic plasticity.