Search

Species distribution models (SDMs) are widely used to predict range shifts but could be unreliable under climate change scenarios because they do not account for evolution. The thermal physiology of a species is a key determinant of its range and thus incorporating thermal trait evolution into SDMs might be expected to alter projected ranges. We identified a genetic basis for physiological and behavioural traits that evolve in response to temperature change in natural populations of threespine stickleback (Gasterosteus aculeatus). Using these data, we created geographical range projections using a mechanistic niche area approach under two climate change scenarios. Under both scenarios, trait data were either static ("no evolution" models), allowed to evolve at observed evolutionary rates ("evolution" models) or allowed to evolve at a rate of evolution scaled by the trait variance that is explained by quantitative trait loci (QTL; "scaled evolution" models). We show that incorporating these traits and their evolution substantially altered the projected ranges for a widespread panmictic marine population, with over 7-fold increases in area under climate change projections when traits are allowed to evolve. Evolution-informed SDMs should improve the precision of forecasting range dynamics under climate change, and aid in their application to management and the protection of biodiversity.

Epigenetic mechanisms underlying phenotypic change are hypothesized to contribute to population persistence and adaptation in the face of environmental change. To date, few studies have explored the heritability of intergenerationally stable methylation levels in natural populations, and little is known about the relative contribution of cis- and trans-regulatory changes to methylation variation. Here, we explore the heritability of DNA methylation, and conduct methylation quantitative trait loci (meQTLs) analysis to investigate the genetic architecture underlying methylation variation between marine and freshwater ecotypes of threespine stickleback (Gasterosteus aculeatus). We quantitatively measured genome-wide DNA methylation in fin tissue using reduced representation bisulfite sequencing of F1 and F2 crosses, and their marine and freshwater source populations. We identified cytosines (CpG sites) that exhibited stable methylation levels across generations. We found that additive genetic variance explained an average of 24–35% of the methylation variance, with a number of CpG sites possibly autonomous from genetic control. We also detected both cis- and trans-meQTLs, with only trans-meQTLs overlapping with previously identified genomic regions of high differentiation between marine and freshwater ecotypes. Finally, we identified the genetic architecture underlying two key CpG sites that were differentially methylated between ecotypes. These findings demonstrate a potential role for DNA methylation in facilitating adaptation to divergent environments and improve our understanding of the heritable basis of population epigenomic variation.

The most remarkable feature of our planet is the diversity of its life forms, ranging from viruses and nanobacteria to blue whales and giant sequoias to satanic leaf‐tailed geckos and leafy seadragons (look them up!). Life is found in essentially all environments on earth, and the number of species living on our planet is many times greater than we could have imagined a century ago. A well‐regarded estimate pegs the number of eukaryotic species on earth at 8.7 million (±1.3 million), of which fewer than 15% are currently described (Mora et al., 2011). The diversity of prokaryotes is less clear (and highly controversial), but an analysis of 1.6 billion 16S ribosomal RNA sequences estimated that 0.8–1.6 million prokaryotic operational taxonomic units exist globally (Louca et al., 2019). While we do not know how many species are currently extant, or have existed in the past, we do know that this biodiversity is valuable, providing food, fibre and medicine, furnishing ecosystem services such as water and air purification, nutrient cycling, pollination and carbon uptake, and contributing to technological innovations ranging from biotechnology to robotics to material science. Moreover, biodiversity underlies the cultural identity of human populations and is important to human health and well‐being. Geographically, species richness increases from the Polar Regions to the tropics in terrestrial and surface marine ecosystems. Thus, some countries, especially those in tropical and subtropical regions, are endowed with much greater biodiversity than others. Unfortunately, benefits arising from the access and utilization of this biodiversity have been unequally shared, with (paradoxically) biodiversity‐poor countries often accruing the lion's share of economic gains. There can be imbalances within countries as well, wherein some segments of the population obtain greater economic benefits from biodiversity and associated traditional knowledge than indigenous peoples. The “Nagoya Protocol on Access to Genetic Resources and the Fair and Equitable Sharing of Benefits Arising from their Utilization,” which came into force in 2014, is an international agreement designed to ensure that the benefits arising from biodiversity are shared equitably (https://www.cbd.int/abs/). However, few scientific journals require compliance with the Nagoya Protocol or the reporting of benefits from biodiversity research. In this editorial, we (the editors of Molecular Ecology and Molecular Ecology Resources) express support for the Nagoya Protocol and the principle of benefit sharing. We believe that scientific journals publishing research on biodiversity can play an important role in implementing the Nagoya Protocol and in reporting on benefits generated from such research. Below, we provide background on the Nagoya Protocol, discuss the kinds of benefits that may arise from biodiversity research, describe the rationale for reporting on these benefits and introduce changes to the journals’ Data Accessibility Statements to incorporate the requirements and goals of the Nagoya Protocol.

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.

The molecular mechanisms underlying behavioural evolution following colonization of novel environments are largely unknown. Molecules that interact to control equilibrium within an organism form physiological regulatory networks. It is essential to determine whether particular components of physiological regulatory networks evolve or if the network as a whole is affected in populations diverging in behavioural responses, as this may affect the nature, amplitude and number of impacted traits. We studied the regulation of four physiological regulatory networks in freshwater and marine populations of threespine stickleback raised in a common environment, which were previously characterized as showing evolutionary divergence in behaviour and stress reactivity. We measured nineteen components of these networks (ligands and receptors) using mRNA and monoamine levels in the brain, pituitary and interrenal gland, as well as hormone levels. Freshwater fish showed higher expression in the brain of adrenergic (adrb2a), serotonergic (htr2a) and dopaminergic (DRD2) receptors, but lower expression of the htr2b receptor. Freshwater fish also showed higher expression of the mc2r receptor of the glucocorticoid axis in the interrenals. Collectively, our results suggest that the inheritance of the regulation of these networks may be implicated in the evolution of behaviour and stress reactivity in association with population divergence. Our results also suggest that evolutionary change in freshwater threespine stickleback may be more associated with the expression of specific receptors rather than with global changes of all the measured constituents of the physiological regulatory networks.

Genome-wide heterozygosity has long been hypothesized to play a role in buffering organisms against developmental perturbations, potentially resulting in increased symmetry. If true, this could in part explain the maintenance of standing genetic variation in wild populations. Marine threespine sticklebacks (Gasterosteus aculeatus) were sampled across their eastern Pacific coastal distribution from Alaska to California and variations in asymmetry for both structural and nonstructural armor traits (lateral plates) were assessed. Structural plates consistently showed less asymmetry than nonstructural plates, but standardized measures of heterozygosity were not correlated with the extent of asymmetry expressed by a fish. Fish that were heterozygous for the major-effect gene controlling lateral plate variation (Ectodysplasin) had higher occurrences of asymmetry, even when the individuals were phenotypically fully plated. Collectively, this suggests that heterozygosity at a major-effect locus can have a greater impact on asymmetry than heterozygosity sampled across the genome.

Conservation and management of endemic species may increasingly involve efforts to prevent hybridization with other species. Native westslope cutthroat trout (Oncorhynchus clarkii lewisi) management in western North America is based largely on admixture estimates with introduced rainbow trout (O. mykiss), with the highest conservation priority given to cutthroat populations that do not exhibit admixture. This study examined the hypothesis that such ancestry quotients are dependent upon the genetic background of reference rainbow trout populations. We used 10 microsatellite loci to estimate admixture within westslope cutthroat trout collected from 39 sites from Alberta, Canada, using three genetically distinct (pairwise FST = 0.100–0.281) rainbow trout genetic backgrounds: a wild (introduced) population from Alberta, two wild (native) populations from British Columbia, and a present-day hatchery broodstock line. Ancestry quotients were significantly impacted by genetic background, whereby the extent of admixture was highest with locally introduced (wild, naturalized) rainbow trout lines and lowest with the hatchery lines. Our results suggest that future studies ought to explore the possibility that local adaptation or drift in introduced rainbow trout populations may contribute to decreased reproductive isolation with geographically proximal cutthroat trout populations.

Size-selective harvest may lead to over-exploitation of commercial fisheries, but the population genetic and evolutionary consequences of such practices remain poorly understood. We investigated the role of within-generation selection in a historically over-exploited Lake Whitefish (Coregonus clupeaformis) population associated with fisheries-induced evolution in Lesser Slave Lake, Alberta, Canada. DNA from archived scales of Lake Whitefish collected between 1986 and 1999 were genotyped at 20 microsatellites and 51 gene-coding SNPs associated with growth and reproduction. We found that the Lake Whitefish in Lesser Slave Lake consisted of a single genetic stock, with microsatellites revealing more temporal than spatial variation in allele frequencies. A comparative genome scan among replicate cohorts from commercially harvested versus random survey samples identified one candidate SNP under divergent selection. This SNP localized within a gene encoding nucleoside diphosphate kinase A, a protein associated with differential growth. Collectively, the results highlight the utility of within-generation genome scans towards investigating the evolutionary consequences of harvest in the wild.

Degraded DNA from suboptimal field sampling is common in molecular ecology. However, its impact on techniques that use restriction site associated next-generation DNA sequencing (RADSeq, GBS) is unknown. We experimentally examined the effects of in situDNA degradation on data generation for a modified double-digest RADSeq approach (3RAD). We generated libraries using genomic DNA serially extracted from the muscle tissue of 8 individual lake whitefish (Coregonus clupeaformis) following 0-, 12-, 48- and 96-h incubation at room temperature posteuthanasia. This treatment of the tissue resulted in input DNA that ranged in quality from nearly intact to highly sheared. All samples were sequenced as a multiplexed pool on an Illumina MiSeq. Libraries created from low to moderately degraded DNA (12-48 h) performed well. In contrast, the number of RADtags per individual, number of variable sites, and percentage of identical RADtags retained were all dramatically reduced when libraries were made using highly degraded DNA (96-h group). This reduction in performance was largely due to a significant and unexpected loss of raw reads as a result of poor quality scores. Our findings remained consistent after changes in restriction enzymes, modified fold coverage values (2- to 16-fold), and additional read-length trimming. We conclude that starting DNA quality is an important consideration for RADSeq; however, the approach remains robust until genomic DNA is extensively degraded.

Disentangling the origin of putatively introduced/invasive species is of increasing priority in conservation biology. The presence of a previously undocumented species may be due to an undetected recent population or range expansion associated with environmental change, or due to an introduction by humans. We used molecular tools to address the origin of a population of threespine stickleback (Gasterosteus aculeatus) first identified in central Alberta, Canada in 1980 from a single lake. We characterized this inland, high elevation Alberta population in comparison to samples from five representative geographic regions worldwide, using mtDNA and nine microsatellite loci to elucidate genetic structure and estimate divergence times. We found significantly lower levels of genetic variation in the Alberta population, which could reflect either a recent colonization or periodic bottlenecks associated with winterkills. While we did find that the Alberta samples were most closely related to the North American West Coast populations, we did not uncover a putative source population. Alberta samples formed a clade in phylogenetic analyses, with divergence time estimates between the Alberta and British Columbia samples ca. 250–750 kya. The hypothesis that the Alberta population represents a natural colonization during North America’s last glacial recession could not be rejected. Collectively, these data suggest that the genetic signature of colonizing populations following introductions may be similar to populations at their range limit, contributing to difficulties establishing population origins.

Phenotypic plasticity is predicted to facilitate individual survival and/or evolve in response to novel environments. Plasticity that facilitates survival should both permit colonization and act as a buffer against further evolution, with contemporary and derived forms predicted to be similarly plastic for a suite of traits. On the other hand, given the importance of plasticity in maintaining internal homeostasis, derived populations that encounter greater environmental heterogeneity should evolve greater plasticity. We tested the evolutionary significance of phenotypic plasticity in coastal British Columbian postglacial populations of threespine stickleback (Gasterosteus aculeatus) that evolved under greater seasonal extremes in temperature after invading freshwater lakes from the sea. Two ancestral (contemporary marine) and two derived (contemporary freshwater) populations of stickleback were raised near their thermal tolerance extremes, 7 and 22 °C. Gene expression plasticity was estimated for more than 14,000 genes. Over five thousand genes were similarly plastic in marine and freshwater stickleback, but freshwater populations exhibited significantly more genes with plastic expression than marine populations. Furthermore, several of the loci shown to exhibit gene expression plasticity have been previously implicated in the adaptive evolution of freshwater populations, including a gene involved in mitochondrial regulation (PPARAa). Collectively, these data provide molecular evidence that highlights the importance of plasticity in colonization and adaptation to new environments.

SummaryThe role of natural selection in the maintenance of genetic variation in wild populations remains a major problem in evolution. The influence of disruptive natural selection on genetic variation is especially interesting because it might lead to the evolution of assortative mating or dominance [1, 2]. In theory, variation can persist at a gene under disruptive natural selection, but the process is little studied and there are few examples [3, 4]. We report a stable polymorphism in the bony armor of threespine stickleback maintained with a deficit of heterozygotes at the major underlying gene, Ectodysplasin (Eda) [5]. The deficit vanishes at the embryo life stage only to re-emerge in adults, indicating that disruptive natural selection, rather than nonrandom mating, is the cause. The mechanism enabling long-term persistence of the polymorphism is unknown, but disruptive selection is predicted to be frequency dependent, favoring homozygous genotypes when they become rare. Further research on the ecological and evolutionary processes affecting individual genes will ultimately lead to a better understanding of the causes of genetic variation in populations.

The discipline of molecular ecology has undergone enormous changes since the journal bearing its name was launched approximately two decades ago. The field has seen great strides in analytical methods development, made groundbreaking discoveries and experienced a revolution in genotyping technology. Here, we provide brief perspectives on the main subdisciplines of molecular ecology, describe key questions and goals, discuss common challenges, predict future research directions and suggest research priorities for the next 20 years., We thank the participants in the Molecular Ecology Symposium and Online Forum for many of the ideas put forward in this article.

Populations at the edge of their range often invoke taxonomic confusion and are increasingly considered to harbour cryptic genetic diversity of significant adaptive potential. In the Peace River region of northwestern Canada, three sculpin species have been reported: spoonhead (Cottus ricei), slimy (Cottus cognatus) and prickly (Cottus asper) sculpin. Prickly sculpin occurrence in this region represents the most eastern edge of its distribution, but its status has remained uncertain following its initial discovery in 1989. These populations may represent an independently evolving lineage of special conservation concern, or be the consequence of an ongoing range expansion, possibly accompanied by interspecific hybridization with local species. Using a combination of mtDNA sequencing and microsatellite analyses, we did not find peripheral population differentiation or interspecific hybridization, suggesting that the Albertan Peace River population belongs to the same genetic group as its western counterparts. Future studies will benefit from a greater understanding of whether demographically independent prickly sculpin populations established in Alberta without the typical genetic signatures of expansion at the periphery of their range.

Understanding the genomic basis of adaptive divergence in the presence of gene flow remains a major challenge in evolutionary biology. In prickly sculpin (Cottus asper), an abundant euryhaline fish in northwestern North America, high genetic connectivity among brackish-water (estuarine) and freshwater (tributary) habitats of coastal rivers does not preclude the build-up of neutral genetic differentiation and emergence of different life history strategies. Because these two habitats present different osmotic niches, we predicted high genetic differentiation at known teleost candidate genes underlying salinity tolerance and osmoregulation. We applied whole genome sequencing of pooled DNA samples (Pool-Seq) to explore adaptive divergence between two estuarine and two tributary habitats. Paired-end sequence reads were mapped against genomic contigs of European Cottus and the gene content of candidate regions was explored based on comparisons with the threespine stickleback genome. Genes showing signals of repeated differentiation among brackish-water and freshwater habitats included functions such as ion-transport and structural permeability in freshwater gills, which suggests that local adaptation to different osmotic niches might contribute to genomic divergence among habitats. Overall, the presence of both repeated and unique signatures of differentiation across many loci scattered throughout the genome is consistent with polygenic adaptation from standing genetic variation and locally variable selection pressures in the early stages of life history divergence. This article is protected by copyright. All rights reserved.

Thermal pollution from industrial processes can have negative impacts on the spawning and development of cold-water fish. Point sources of thermal effluent may need to be managed to avoid affecting discrete populations. Correspondingly, we examined fine-scale ecological and genetic population structure of two whitefish species (Coregonus clupeaformis and Prosopium cylindraceum) on Lake Huron, Canada, in the immediate vicinity of thermal effluent from nuclear power generation. Niche metrics using δ13C and δ15N stable isotopes showed high levels of overlap (48.6 to 94.5%) in resource use by adult fish captured in areas affected by thermal effluent compared to nearby reference locations. Isotopic niche size, a metric of resource use diversity, was 1.3- to 2.8-fold higher than reference values in some thermally affected areas, indicative of fish mixing. Microsatellite analyses of genetic population structure (Fst, STRUCTURE and DAPC) indicated that fish captured at all locations in the vicinity of the power plant were part of a larger population extending beyond the study area. In concert, ecological and genetic markers do not support the presence of an evolutionarily significant unit in the vicinity of the power plant. Thus, future research should focus on the potential impacts of thermal emissions on development and recruitment.

Despite recent progress, we still know relatively little about the genetic architecture that underlies adaptation to divergent environments. Determining whether the genetic architecture of phenotypic adaptation follows any predictable patterns requires data from a wide variety of species. However, in many organisms, genetic studies are hindered by the inability to perform genetic crosses in the laboratory or by long generation times. Admixture mapping is an approach that circumvents these issues by taking advantage of hybridization that occurs between populations or species in the wild. Here, we demonstrate the utility of admixture mapping in a naturally occurring hybrid population of threespine sticklebacks (Gasterosteus aculeatus) from Enos Lake, British Columbia. Until recently, this lake contained two species of sticklebacks adapted to divergent habitats within the lake. This benthic-limnetic species pair diverged in a number of phenotypes, including male nuptial coloration and body shape, which were previously shown to contribute to reproductive isolation between them. However, recent ecological disturbance has contributed to extensive hybridization between the species, and there is now a single, admixed population within Enos Lake. We collected over 500 males from Enos Lake and found that most had intermediate nuptial color and body shape. By genotyping males with nuptial color at the two extremes of the phenotypic distribution, we identified seven genomic regions on three chromosomes associated with divergence in male nuptial color. These genomic regions are also associated with variation in body shape, suggesting that tight linkage and/or pleiotropy facilitated adaptation to divergent environments in this benthic-limnetic species pairs.

Transition of an evolving population to a new adaptive optimum is predicted to leave a signature in the distribution of effect sizes of fixed mutations. If they affect many traits (are pleiotropic), large effect mutations should contribute more when a population evolves to a farther adaptive peak than to a nearer peak. We tested this prediction in wild threespine stickleback fish (Gasterosteus aculeatus) by comparing the estimated frequency of large effect genetic changes underlying evolution as the same ancestor adapted to two lake types since the end of the ice age. A higher frequency of large effect genetic changes (quantitative trait loci) contributed to adaptive evolution in populations that adapted to lakes representing a more distant optimum than to lakes in which the optimum phenotype was nearer to the ancestral state. Our results also indicate that pleiotropy, not just optimum overshoot, contributes to this difference. These results suggest that a series of adaptive improvements to a new environment leaves a detectable mark in the genome of wild populations. Although not all assumptions of the theory are likely met in natural systems, the prediction may be robust enough to the complexities of natural environments to be useful when forecasting adaptive responses to large environmental changes.

As populations adapt to novel environments, divergent selection will promote heterogeneous genomic differentiation via reductions in gene flow for loci underlying adaptive traits. Using a data set of over 100 SNP markers, genome scans were performed to investigate the effect of natural selection maintaining differentiation in five lakes harbouring sympatric pairs of normal and dwarf lake whitefish (Coregonus clupeaformis). A variable proportion of SNPs (between 0% and 12%) was identified as outliers, which corroborated the predicted intensity of competitive interactions unique to each lake. Moreover, strong reduction in heterozygosity was typically observed for outlier loci in dwarf but not in normal whitefish, indicating that directional selection has been acting on standing genetic variation more intensively in dwarf whitefish. SNP associations in backcross hybrid progeny identified 16 genes exhibiting genotype-phenotype associations for four adaptive traits (growth, swimming activity, gill rakers and condition factor). However, neither simple relationship between elevated levels of genetic differentiation with adaptive phenotype nor conspicuous genetic signatures for parallelism at outlier loci were detected, which underscores the importance of independent evolution among lakes. The integration of phenotypic, transcriptomic and functional genomic information identified two candidate genes (sodium potassium ATPase and triosephosphate isomerase) involved in the recent ecological divergence of lake whitefish. Finally, the identification of several markers under divergent selection suggests that many genes, in an environment-specific manner, are recruited by selection and ultimately contributed to the repeated ecological speciation of a dwarf phenotype.

All-female 'species' of fish have been shown to be great models in ecological and evolutionary studies because of the insights they can provide into the origin and evolution of asexuality, the ecology of hybrids, associations between genotype and environment, and the maintenance of sex. Gynogenetic organisms that evolved from sexual ancestors, and combine the disadvantageous traits from sexuality and asexuality, have long baffled evolutionary biologists trying to understand their origin and persistence with their sympatric sexual counterparts. In this issue, a new study using an integrated molecular phylogenetic and classical genetic approach has uncovered compelling evidence regarding the obscure asexual origin of the Amazon molly, Poecilia formosa. By performing an extensive phylogeographic analysis, Stock et al. (2010) provide evidence that the Amazon molly arose only once within its history, with monophyly being strongly supported by mitochondrial DNA and microsatellite analyses. This result, combined with an elaborate failed attempt to resynthesize the lineage, suggests that vertebrate gynogens such as the Amazon molly are not rare because they are at a disadvantage to their sexual counterparts, but because the genomic conditions under which they arise are rare. Organisms that apparently combine the disadvantages of both sexuality and asexuality remain difficult to understand from both an ecological and an evolutionary perspective, and Stock et al. (2010) highlight several outstanding important questions. Nonetheless, given that we now have a better knowledge of the origin and history of this unique 'species', this should allow researchers to better understand how these frozen F1's can persist amidst the masterpiece of nature.

Ecological speciation is the process by which barriers to gene flow between populations evolve due to adaptive divergence via natural selection. A relatively unexplored area in ecological speciation is the role of gene expression. Gene expression may be associated with ecologically important phenotypes not evident from morphology and play a role during colonization of new environments. Here we review two potential roles of gene expression in ecological speciation: (1) its indirect role in facilitating population persistence and (2) its direct role in contributing to genetically based reproductive isolation. We find indirect evidence that gene expression facilitates population persistence, but direct tests are lacking. We also find clear examples of gene expression having effects on phenotypic traits and adaptive genetic divergence, but links to the evolution of reproductive isolation itself remain indirect. Gene expression during adaptive divergence seems to often involve complex genetic architectures controlled by gene networks, regulatory regions, and “eQTL hotspots.” Nonetheless, we review how approaches for isolating the functional mutations contributing to adaptive divergence are proving to be successful. The study of gene expression has promise for increasing our understanding ecological speciation, particularly when integrative approaches are applied.

Growing knowledge of the molecular basis of adaptation in wild populations is expanding the study of natural selection. We summarize ongoing efforts to infer three aspects of natural selection—mechanism, form and history—from the genetics of adaptive evolution in threespine stickleback that colonized freshwater after the last ice age. We tested a mechanism of selection for reduced bony armour in freshwater by tracking genotype and allele frequency changes at an underlying major locus (Ectodysplasin) in transplanted stickleback populations. We inferred disruptive selection on genotypes at the same locus in a population polymorphic for bony armour. Finally, we compared the distribution of phenotypic effect sizes of genes underlying changes in body shape with that predicted by models of adaptive peak shifts following colonization of freshwater. Studies of the effects of selection on genes complement efforts to identify the molecular basis of adaptive differences, and improve our understanding of phenotypic evolution.

The accessibility of new genomic resources, high-throughput molecular technologies and analytical approaches such as genome scans have made finding genes contributing to fitness variation in natural populations an increasingly feasible task. Once candidate genes are identified, we argue that it is necessary to take a mechanistic approach and work up through the levels of biological organization to fully understand the impacts of genetic variation at these candidate genes. We demonstrate how this approach provides testable hypotheses about the causal links among levels of biological organization, and assists in designing relevant experiments to test the effects of genetic variation on phenotype, whole-organism performance capabilities and fitness. We review some of the research programs that have incorporated mechanistic approaches when examining naturally occurring genetic and phenotypic variation and use these examples to highlight the value of developing a comprehensive understanding of the relationship between genotype and fitness. We give suggestions to guide future research aimed at uncovering and understanding the genetic basis of adaptation and argue that further integration of mechanistic approaches will help molecular ecologists better understand the evolution of natural populations.

We used microarrays and a previously established linkage map to localize the genetic determinants of brain gene expression for a backcross family of lake whitefish species pairs (Coregonus sp.). Our goals were to elucidate the genomic distribution and sex specificity of brain expression QTL (eQTL) and to determine the extent to which genes controlling transcriptional variation may underlie adaptive divergence in the recently evolved dwarf (limnetic) and normal (benthic) whitefish. We observed a sex bias in transcriptional genetic architecture, with more eQTL observed in males, as well as divergence in genome location of eQTL between the sexes. Hotspots of nonrandom aggregations of up to 32 eQTL in one location were observed. We identified candidate genes for species pair divergence involved with energetic metabolism, protein synthesis, and neural development on the basis of colocalization of eQTL for these genes with eight previously identified adaptive phenotypic QTL and four previously identified outlier loci from a genome scan in natural populations. Eighty-eight percent of eQTL-phenotypic QTL colocalization involved growth rate and condition factor QTL, two traits central to adaptive divergence between whitefish species pairs. Hotspots colocalized with phenotypic QTL in several cases, revealing possible locations where master regulatory genes, such as a zinc-finger protein in one case, control gene expression directly related to adaptive phenotypic divergence. We observed little evidence of colocalization of brain eQTL with behavioral QTL, which provides insight into the genes identified by behavioral QTL studies. These results extend to the transcriptome level previous work illustrating that selection has shaped recent parallel divergence between dwarf and normal lake whitefish species pairs and that metabolic, more than morphological, differences appear to play a key role in this divergence.

Mapping of expression quantitative trait loci (eQTL) is a powerful means for elucidating the genetic architecture of gene regulation. Yet, eQTL mapping has not been applied toward investigating the regulation architecture of genes involved in the process of population divergence, ultimately leading to speciation events. Here, we conducted an eQTL mapping experiment to compare the genetic architecture of transcript regulation in adaptive traits, differentiating the recently evolved limnetic (dwarf) and benthic (normal) species pairs of lake whitefish. The eQTL were mapped in three data sets derived from an F1 hybrid-dwarf backcrossed family: the entire set of 66 genotyped individuals and the two sexes treated separately. We identified strikingly more eQTL in the female data set (174), compared to both male (54) and combined (33) data sets. The majority of these genes were not differentially expressed between male and female progeny of the backcross family, thus providing evidence for a strong pleiotropic sex-linked effect in transcriptomic regulation. The subtelomeric region of a linkage group segregating in females encompassed >50% of all eQTL, which exhibited the most pronounced additive effects. We also conducted a direct comparison of transcriptomic profiles between pure dwarf and normal progeny reared in controlled conditions. We detected 34 differentially expressed transcripts associated with eQTL segregating only in sex-specific data sets and mostly belonging to functional groups that differentiate dwarf and normal whitefish in natural populations. Therefore, these eQTL are not related to interindividual variation, but instead to the adaptive and historical genetic divergence between dwarf and normal whitefish. This study exemplifies how the integration of genetic and transcriptomic data offers a strong means for dissecting the functional genomic response to selection by separating mapping family-specific effects from genetic factors under selection, potentially involved in the phenotypic divergence of natural populations.

Adaptive evolutionary change is contingent on variation and selection; thus, understanding adaptive divergence and ultimately speciation requires information on both the genetic basis of adaptive traits as well as an understanding of the role of divergent natural selection on those traits. The lake whitefish (Coregonus clupeaformis) consists of several sympatric "dwarf" (limnetic) and normal (benthic) species pairs that co-inhabit northern postglacial lakes. These young species pairs have evolved independently and display parallelism in life history, behavioral, and morphological divergence associated with the use of distinct trophic resources. We identified phenotype-environment associations and determined the genetic architecture and the role of selection modulating population genetic divergence in sympatric dwarf and normal lake whitefish. The genetic architecture of 9 adaptive traits was analyzed in 2 hybrid backcrosses individually phenotyped throughout their life history. Significant quantitative trait loci (QTL) were associated with swimming behavior (habitat selection and predator avoidance), growth rate, morphology (condition factor and gill rakers), and life history (onset of maturity and fecundity). Genome scans among 4 natural sympatric pairs, using loci segregating in the map, revealed a signature of selection for 24 loci. Loci exhibiting a signature of selection were associated with QTL relative to other regions of the genome more often than expected by chance alone. Two parallel QTL outliers for growth and condition factor exhibited segregation distortion in both mapping families, supporting the hypothesis that adaptive divergence contributing to parallel reductions of gene flow among natural populations may cause genetic incompatibilities. Overall, these findings offer evidence that the genetic architecture of ecological speciation is associated with signatures of selection in nature, providing strong support for the hypothesis that divergent natural selection is currently maintaining adaptive differentiation and promoting ecological speciation in lake whitefish species pairs.

Elucidating the genetic architecture of population divergence may reveal the evolution of reproductive barriers and the genomic regions implicated in the process. We assembled genetic linkage maps for the dwarf and Normal lake whitefish species complex and their hybrids. A total of 877 AFLP loci and 30 microsatellites were positioned. The homology of mapped loci between families supported the existence of 34 linkage groups (of 40n expected) exhibiting 83% colinearity among linked loci between these two families. Classes of AFLP markers were not randomly distributed among linkage groups. Both AFLP and microsatellites exhibited deviations from Mendelian expectations, with 30.4% exhibiting significant segregation distortion across 28 linkage groups of the four linkage maps in both families (P < 0.00001). Eight loci distributed over seven homologous linkage groups were significantly distorted in both families and the level of distortion, when comparing homologous loci of the same phase between families, was correlated (Spearman R = 0.378, P = 0.0021). These results suggest that substantial divergence incurred during allopatric glacial separation and subsequent sympatric ecological specialization has resulted in several genomic regions that are no longer complementary between dwarf and Normal populations issued from different evolutionary glacial lineages.

Parallel changes in body shape may evolve in response to similar environmental conditions, but whether such parallel phenotypic changes share a common genetic basis is still debated. The goal of this study was to assess whether parallel phenotypic changes could be explained by genetic parallelism, multiple genetic routes, or both. We first provide evidence for parallelism in fish shape by using geometric morphometrics among 300 fish representing five species pairs of Lake Whitefish. Using a genetic map comprising 3438 restriction site−associated DNA sequencing single-nucleotide polymorphisms, we then identified quantitative trait loci underlying body shape traits in a backcross family reared in the laboratory. A total of 138 body shape quantitative trait loci were identified in this cross, thus revealing a highly polygenic architecture of body shape in Lake Whitefish. Third, we tested for evidence of genetic parallelism among independent wild populations using both a single-locus method (outlier analysis) and a polygenic approach (analysis of covariation among markers). The single-locus approach provided limited evidence for genetic parallelism. However, the polygenic analysis revealed genetic parallelism for three of the five lakes, which differed from the two other lakes. These results provide evidence for both genetic parallelism and multiple genetic routes underlying parallel phenotypic evolution in fish shape among populations occupying similar ecological niches.

As natural selection must act on underlying genetic variation, discovering the number and location of loci under the influence of selection is imperative towards understanding adaptive divergence in evolving populations. Studies employing genome scans have hypothesized that the action of divergent selection should reduce gene flow at the genomic locations implicated in adaptation and speciation among natural populations, yet once ‘outlier’ patterns of variation have been identified the function and role of such loci needs to be confirmed. We integrated adaptive QTL mapping and genomic scans among diverging sympatric pairs of the lake whitefish (Coregonus clupeaformis) species complex in order to test the hypothesis that differentiation between dwarf and normal ecotypes at growth-associated QTL was maintained by directional selection. We found evidence of significantly high levels of molecular divergence among eight growth QTL where two of the strongest candidate loci under the influence of directional selection exhibited parallel reductions of gene flow over multiple populations.

Coregonine fish represent the most successful evolutionary lineage of salmonids with Coregonus as the most speciose salmonid genus inhabiting numerous postglacial lakes across the northern hemisphere. We isolated and characterized 31 polymorphic microsatellite loci in Coregonus clupeaformis with an average number of 5.3 alleles per locus (range three to eight) and an overall expected heterozygosity of 0.74 ± 0.11. Two loci revealed significant linkage associations through analyses of mapping families. Six additional salmonid taxa assessed for cross-species amplification revealed between 18 and 26 positive amplifications and between two and 12 polymorphic loci per species.

The phylogeography history and contemporary agents of selection for many marine fisheries, characterized by widespread species distributions in the face of significant harvest, remains poorly understood. Chub mackerel (Scomber japonicus) are a widespread species in the Indo-Pacific and represent one of the top five commercially fished species in the world, yet their phylogeographic history remains unknown. We characterized the genetic diversity, structure and demographic history of S. japonicus throughout adjacent Chinese seas (from the Yellow Sea to the South China Sea). Using 220 individuals from 11 sites, we inferred 55 distinct haplotypes from complete mitochondrial cytochrome b gene sequences. Haplotype diversity ranged from 0.505 to 0.967 and nucleotide diversity ranged from 0.00056 to 0.01042. Genetic differentiation (F

The concept of the designatable unit (DU) affords a practical approach to identifying diversity below the species level for conservation prioritization. However, its suitability for defining conservation units in ecologically diverse, geographically widespread and taxonomically challenging species complexes has not been broadly evaluated. The lake whitefish species complex (Coregonus spp.) is geographically widespread in the Northern Hemisphere, and it contains a great deal of variability in ecology and evolutionary legacy within and among populations, as well as a great deal of taxonomic ambiguity. Here, we employ a set of hierarchical criteria to identify DUs within the Canadian distribution of the lake whitefish species complex. We identified 36 DUs based on (i) reproductive isolation, (ii) phylogeographic groupings, (iii) local adaptation and (iv) biogeographic regions. The identification of DUs is required for clear discussion regarding the conservation prioritization of lake whitefish populations. We suggest conservation priorities among lake whitefish DUs based on biological consequences of extinction, risk of extinction and distinctiveness. Our results exemplify the need for extensive genetic and biogeographic analyses for any species with broad geographic distributions and the need for detailed evaluation of evolutionary history and adaptive ecological divergence when defining intraspecific conservation units.

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.

Adaptation and reproductive isolation, the engines of biological diversity, are still elusive when discussing the genetic bases of speciation. Namely, the number of genes and magnitude of selection acting positively or negatively on genomic traits implicated in speciation is contentious. Here, we describe the first steps of an ongoing research program aimed at understanding the genetic bases of population divergence and reproductive isolation in the lake whitefish (Coregonus clupeaformis). A preliminary linkage map originating from a hybrid cross between dwarf and normal ecotypes is presented, whereby some of the segregating AFLP markers were found to be conserved among natural populations. Maximum-likelihood was used to estimate hybrid indices from non-diagnostic markers at 998 AFLP loci. This allowed identification of the most likely candidate loci that have been under the influence of selection during the natural hybridisation of whitefish originating from different glacial races. As some of these loci could be identified on the linkage map, the possibility that selection of traits in natural populations may eventually be correlated to specific chromosomal regions was demonstrated. The future prospects and potential of these approaches to elucidate the genetic bases of adaptation and reproductive isolation among sympatric ecotypes of lake whitefish is discussed.

Experimental estimates of the effects of selection on genes determining adaptive traits add to our understanding of the mechanisms of evolution. We measured selection on genotypes of the Ectodysplasin locus, which underlie differences in lateral plates in threespine stickleback fish. A derived allele (low) causing reduced plate number has been fixed repeatedly after marine stickleback colonized freshwater from the sea, where the ancestral allele (complete) predominates. We transplanted marine sticklebacks carrying both alleles to freshwater ponds and tracked genotype frequencies over a generation. The low allele increased in frequency once lateral plates developed, most likely via a growth advantage. Opposing selection at the larval stage and changing dominance for fitness throughout life suggest either that the gene affects additional traits undergoing selection or that linked loci also are affecting fitness.

• Premise of the study: Microsatellite markers were isolated and characterized in Plectritis congesta for studying the evolution of this highly variable species. • METHODS: and Results: We used 454 sequencing of DNA enriched for microsatellite repeats to develop microsatellite markers. This produced 262079 reads with an average length of 324 bp, representing approximately 800 microsatellite regions from which 48 primers were tested. Eleven markers reliably amplified without optimization. These primer pairs showed a high degree of heterozygosity and allelic diversity. Unexpectedly, half of the markers contained multiple peaks, with up to four alleles per individual, which suggests that either polyploidy or isolated gene duplication has occurred within this clade. These primers successfully cross-amplified in P. macrocera, indicating the utility of these markers for the genus. • CONCLUSIONS: With variation in mating system and habitat, a mix of duplicated and nonduplicated markers, and high genetic variance, Plectritis is an ideal candidate model genus for studying the ecological and evolutionary consequences of gene duplication.

Adaptive radiation occurs when divergent natural selection in different environments leads to phenotypic differentiation. The pleiotropic effects of underlying genes can either promote or constrain this diversification. Identifying the pleiotropic effects of genes responsible for divergent traits, and testing how the environment influences these effects, can therefore help to provide an understanding of how ecology drives evolutionary change between populations. Positive selection on low-armor alleles at the Ectodysplasin (Eda) locus in threespine stickleback has led to the repeated evolution of reduced armor in populations following freshwater colonization by fully armored marine sticklebacks. Here, we demonstrate that Eda has environmentally determined pleiotropic effects on armor and growth. When raised in freshwater, reduced armor sticklebacks carrying "low" alleles at Eda had increased growth rate relative to fully armored sticklebacks carrying "complete" alleles. In saltwater treatments this growth advantage was present during juvenile growth but lost during adult growth, suggesting that in this environment stickleback are able to develop full armor plates without sacrificing overall growth rate. The environment specific pleiotropic effects of Eda demonstrate that ecological factors can mediate the influence of genetic architecture in driving phenotypic evolution. Furthermore, because size is important for mate choice in stickleback, the growth rate differences influenced by Eda may have effects on reproductive isolation between marine and freshwater populations.

Studies of phenotype-environment associations in adaptive radiation have focused largely on morphological traits related to resource-based phenotypic differences. The genetic basis of adaptive behaviors implicated in population divergence remains poorly understood, as few studies have tested the hypothesis of behavioral phenotype-environment associations. We provide evidence of a phenotype-environment association for differential adaptive swimming behaviors through experiments conducted on dwarf, normal, and hybrid lake whitefish (Coregonus clupeaformis). Highly significant differences were observed for depth selection, directional changes, and burst swimming, implicating a genetic basis for these behaviors. Hybrid crosses revealed that depth selection is under additive genetic control, while dominance effects were suggested for directional changes and burst swimming. Estimates for the genetic basis of behavioral differentiation from an animal model were consistent with these observations. Comparative estimates of behavioral differentiation (Q(ST)) against neutral expectations (F(ST)) revealed pronounced departures from neutral expectations in all three behavioral phenotypes, consistent with the hypothesis that directional selection has driven the divergence of behavior in dwarf and normal lake whitefish ecotypes.