Evolution of stickleback in 50 years on earthquake-uplifted islands

Parallel evolution despite low genetic diversity in three-spined sticklebacks

When populations repeatedly adapt to similar environments they can evolve similar phenotypes based on shared genetic mechanisms (parallel evolution). The likelihood of parallel evolution is affected by demographic history, as it depends on the standing genetic variation of the source population. The three-spined stickleback (Gasterosteus aculeatus) repeatedly colonized and adapted to brackish and freshwater. Most parallel evolution studies in G. aculeatus were conducted at high latitudes, where freshwater populations maintain connectivity to the source marine populations. Here, we analysed southern and northern European marine and freshwater populations to test two hypotheses. First, that southern European freshwater populations (which currently lack connection to marine populations) lost genetic diversity due to bottlenecks and inbreeding compared to their northern counterparts. Second, that the degree of genetic parallelism is higher among northern than southern European freshwater populations, as the latter have been subjected to strong drift due to isolation. The results show that southern populations exhibit lower genetic diversity but a higher degree of genetic parallelism than northern populations. Hence, they confirm the hypothesis that southern populations have lost genetic diversity, but this loss probably happened after they had already adapted to freshwater conditions, explaining the high degree of genetic parallelism in the south.

Morphological description of spontaneous pelvic fin loss in a neotropical cichlid fish

Pelvic fins are a characteristic structure of the vertebrate Bauplan. Yet, pelvic fin loss has occurred repeatedly across a wide diversity of other lineages of tetrapods and at least 48 times in teleost fishes. This pelvic finless condition is often associated with other morphological features such as body elongation, loss of additional structures, and bilateral asymmetry. However, despite the remarkable diversity in the several thousand cichlid fish species, none of them are characterized by the complete absence of pelvic fins. Here, we examined the musculoskeletal structure and associated bilateral asymmetry in Midas cichlids (Amphilophus cf. citrinellus) that lost their pelvic fins spontaneously in the laboratory. Due to this apparent mutational loss of the pelvic girdle and fins, the external and internal anatomy are described in a series of "normal" Midas individuals and their pelvic finless sibling tankmates. First, other traits associated with teleost pelvic fin loss, the genetic basis of pelvic fin loss, and the potential for pleiotropic effects of these genes on other traits in teleosts were all reviewed. Using these traits as a guide, we investigated whether other morphological differences were associated with the pelvic girdle/fin loss. The mean values of the masses of muscle of the pectoral fin, fin ray numbers in the unpaired fins, and oral jaw tooth numbers did not differ between the two pelvic fin morphotypes. However, significant differences in meristic values of the paired traits assessed were observed for the same side of the body between morphotypes. Notably, bilateral asymmetry was found exclusively for the posterior lateral line scales. Finally, we found limited evidence of pleiotropic effects, such as lateral line scale numbers and fluctuating asymmetry between the Midas pelvic fin morphotypes. The fast and relatively isolated changes in the Midas cichlids suggest minor but interesting pleiotropic effects could accompany loss of cichlid pelvic fins.

Extreme events, trophic chain reactions, and shifts in phenotypic selection

Demographic consequences of rapid environmental change and extreme climatic events (ECEs) can cascade across trophic levels with evolutionary implications that have rarely been explored. Here, we show how an ECE in high Arctic Svalbard triggered a trophic chain reaction, directly or indirectly affecting the demography of both overwintering and migratory vertebrates, ultimately inducing a shift in density-dependent phenotypic selection in migratory geese. A record-breaking rain-on-snow event and ice-locked pastures led to reindeer mass starvation and a population crash, followed by a period of low mortality and population recovery. This caused lagged, long-lasting reductions in reindeer carrion numbers and resultant low abundances of Arctic foxes, a scavenger on reindeer and predator of migratory birds. The associated decrease in Arctic fox predation of goose offspring allowed for a rapid increase in barnacle goose densities. As expected according to r- and K-selection theory, the goose body condition (affecting reproduction and post-fledging survival) maximising Malthusian fitness increased with this shift in population density. Thus, the winter ECE acting on reindeer and their scavenger, the Arctic fox, indirectly selected for higher body condition in migratory geese. This high Arctic study provides rare empirical evidence of links between ECEs, community dynamics and evolution, with implications for our understanding of indirect eco-evolutionary impacts of global change.

On the Origins of Phenotypic Parallelism in Benthic and Limnetic Stickleback

Rapid evolution of similar phenotypes in similar environments, giving rise to in situ parallel adaptation, is an important hallmark of ecological speciation. However, what appears to be in situ adaptation can also arise by dispersal of divergent lineages from elsewhere. We test whether two contrasting phenotypes repeatedly evolved in parallel, or have a single origin, in an archetypal example of ecological adaptive radiation: benthic-limnetic three-spined stickleback (Gasterosteus aculeatus) across species pair and solitary lakes in British Columbia. We identify two genomic clusters across freshwater populations, which differ in benthic-limnetic divergent phenotypic traits and separate benthic from limnetic individuals in species pair lakes. Phylogenetic reconstruction and niche evolution modeling both suggest a single evolutionary origin for each of these clusters. We detected strong phylogenetic signal in benthic-limnetic divergent traits, suggesting that they are ancestrally retained. Accounting for ancestral state retention, we identify local adaptation of body armor due to the presence of an intraguild predator, the sculpin (Cottus asper), and environmental effects of lake depth and pH on body size. Taken together, our results imply a predominant role for retention of ancestral characteristics in driving trait distribution, with further selection imposed on some traits by environmental factors.

Sex-specific transgenerational plasticity: developmental temperatures of mothers and fathers have different effects on sons and daughters

Each parent can influence offspring phenotype via provisioning of the zygote or sex-specific DNA methylation. Transgenerational plasticity may therefore depend on the environmental conditions experienced by each parent. We tested this hypothesis by conducting a fully factorial experiment across three generations of guppies (Poecilia reticulata), determining the effects of warm (28°C) and cold (21°C) thermal backgrounds of mothers and fathers on mass and length, and thermal performance (sustained and sprint swimming speeds, citrate synthase and lactate dehydrogenase activities; 18, 24, 28, 32 and 36°C test temperatures) of sons and daughters. Offspring sex was significant for all traits except for sprint speed. Warmer mothers produced sons and daughters with reduced mass and length, and warmer fathers produced shorter sons. Sustained swimming speed (Ucrit) of male offspring was greatest when both parents were raised at 28°C, and warmer fathers produced daughters with greater Ucrit. Similarly, warmer fathers produced sons and daughters with greater metabolic capacity. We show that the thermal variation experienced by parents can modify offspring phenotype, and that predicting the impacts of environmental change on populations would require knowledge of the thermal background of each mother and father, particularly where sexes are spatially segregated.

How can physiology best contribute to wildlife conservation in a warming world?

Global warming is now predicted to exceed 1.5°C by 2033 and 2°C by the end of the 21st century. This level of warming and the associated environmental variability are already increasing pressure on natural and human systems. Here we emphasize the role of physiology in the light of the latest assessment of climate warming by the Intergovernmental Panel on Climate Change. We describe how physiology can contribute to contemporary conservation programmes. We focus on thermal responses of animals, but we acknowledge that the impacts of climate change are much broader phylogenetically and environmentally. A physiological contribution would encompass environmental monitoring, coupled with measuring individual sensitivities to temperature change and upscaling these to ecosystem level. The latest version of the widely accepted Conservation Standards designed by the Conservation Measures Partnership includes several explicit climate change considerations. We argue that physiology has a unique role to play in addressing these considerations. Moreover, physiology can be incorporated by institutions and organizations that range from international bodies to national governments and to local communities, and in doing so, it brings a mechanistic approach to conservation and the management of biological resources.

Phylogenetic relationships of the North American catfishes (Ictaluridae, Siluriformes): Investigating the origins and parallel evolution of the troglobitic species

Insular habitats have played an important role in developing evolutionary theory, including natural selection and island biogeography. Caves are insular habitats that place extreme selective pressures on organisms due to the absence of light and food scarcity. Therefore, cave organisms present an excellent opportunity for studying colonization and speciation in response to the unique abiotic conditions that require extreme adaptations. One vertebrate family, the North American catfishes (Ictaluridae), includes four troglobitic species that inhabit the karst region bordering the western Gulf of Mexico. The phylogenetic relationships of these species have been contentious, and conflicting hypotheses have been proposed to explain their origins. The purpose of our study was to construct a time-calibrated phylogeny of Ictaluridae using first-occurrence fossil data and the largest molecular dataset on the group to date. We test the hypothesis that troglobitic ictalurids have evolved in parallel, thus resulting from repeated cave colonization events. We found that Prietella lundbergi is sister to surface-dwelling Ictalurus and that Prietella phreatophila + Trogloglanis pattersoni are sister to surface-dwelling Ameiurus, suggesting that ictalurids colonized subterranean habitats at least twice in evolutionary history. The sister relationship between Prietella phreatophila and Trogloglanis pattersoni may indicate that these two species diverged from a common ancestor following a subterranean dispersal event between Texas and Coahuila aquifers. We recovered Prietella as a polyphyletic genus and recommend P. lundbergi be removed from this genus. With respect to Ameiurus, we found evidence for a potentially undescribed species sister to A. platycephalus, which warrants further investigation of Atlantic and Gulf slope Ameiurus species. In Ictalurus, we identified shallow divergence between I. dugesii and I. ochoterenai, I. australis and I. mexicanus, and I. furcatus and I. meridionalis, indicating a need to reexamine the validity of each species. Lastly, we propose minor revisions to the intrageneric classification of Noturus including the restriction of subgenus Schilbeodes to N. gyrinus (type species), N. lachneri, N. leptacanthus, and N. nocturnus.

Dimensionality and Modularity of Adaptive Variation: Divergence in Threespine Stickleback from Diverse Environments

AbstractPopulations are subjected to diverse environmental conditions that affect fitness and induce evolutionary or plastic responses, resulting in phenotypic divergence. Some authors contend that such divergence is concentrated along a single major axis of trait covariance even if that axis does not lead populations directly toward a fitness optimum. Other authors argue that divergence can occur readily along many phenotype axes at the same time. We use populations of threespine stickleback (Gasterosteus aculeatus) from 14 lakes with contrasting ecological conditions to find some resolution along the continuum between these two extremes. Unlike many previous studies, we included several functional suites of traits (defensive, swimming, trophic) potentially subject to different sources of selection. We find that populations exhibit dimensionality of divergence that is high enough to preclude a history of constraint along a single axis-both for divergence in multivariate mean trait values and for the structure of trait covariances. Dimensionality varied among trait suites and were strongly influenced by the inclusion of specific traits, and integration of trait suites varied between populations. We leverage this variation into new insights about the process of divergence and suggest that similar analyses could increase understanding of other adaptive radiations.

Selection on early survival does not explain germination rate clines in Mimulus cardinalis

PREMISE

Many traits covary with environmental gradients to form phenotypic clines. While local adaptation to the environment can generate phenotypic clines, other nonadaptive processes may also. If local adaptation causes phenotypic clines, then the direction of genotypic selection on traits should shift from one end of the cline to the other. Traditionally, genotypic selection on non-Gaussian traits like germination rate have been hampered because it is challenging to measure their genetic variance.

METHODS

Here we used quantitative genetics and reciprocal transplants to test whether a previously discovered cline in germination rate showed additional signatures of adaptation in the scarlet monkeyflower (Mimulus cardinalis). We measured genotypic and population level covariation between germination rate and early survival, a component of fitness. We developed a novel discrete log-normal model to estimate genetic variance in germination rate.

RESULTS

Contrary to our adaptive hypothesis, we found no evidence that genetic variation in germination rate contributed to variation in early survival. Across populations, southern populations in both gardens germinated earlier and survived more.

CONCLUSIONS

Southern populations have higher early survival but it is not caused by faster germination. This pattern is consistent with nonadaptive forces driving the phenotypic cline in germination rate, but future work will need to assess whether there is selection at other life stages. This statistical framework should help expand quantitative genetic analyses for other waiting-time traits.

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

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

Selection in the city: Rapid and fine-scale evolution of urban eastern water dragons

Oceanic archipelagos have long been treated as a Petri dish for studies of evolutionary and ecological processes. Like archipelagos, cities exhibit similar patterns and processes, such as the rapid phenotypic divergence of a species between urban and nonurban environments. However, on a local scale, cities can be highly heterogenous, where geographically close populations can experience dramatically different environmental conditions. Nevertheless, we are yet to understand the evolutionary and ecological implications for populations spread across a heterogenous cityscape. To address this, we compared neutral genetic divergence to quantitative trait divergence within three native riparian and four city park populations of an iconic urban adapter, the eastern water dragon. We demonstrated that selection is likely acting to drive divergence of snout-vent length and jaw width across native riparian populations that are geographically isolated and across city park populations that are geographically close yet isolated by urbanization. City park populations as close as 0.9 km exhibited signs of selection-driven divergence to the same extent as native riparian populations isolated by up to 114.5 km. These findings suggest that local adaptation may be occurring over exceptionally small geographic and temporal scales within a single metropolis, demonstrating that city parks can act as archipelagos for the study of rapid evolution.

Parallel recolonizations generate distinct genomic sectors in kelp following high-magnitude earthquake disturbance

Large-scale disturbance events have the potential to drastically reshape biodiversity patterns. Notably, newly vacant habitat space cleared by disturbance can be colonized by multiple lineages, which can lead to the evolution of distinct spatial "sectors" of genetic diversity within a species. We test for disturbance-driven sectoring of genetic diversity in intertidal southern bull kelp, Durvillaea antarctica (Chamisso) Hariot, following the high-magnitude 1855 Wairarapa earthquake in New Zealand. Specifically, we use genotyping-by-sequencing (GBS) to analyse fine-scale population structure across the uplift zone and apply machine learning to assess the fit of alternative recolonizaton models. Our analysis reveals that specimens from the uplift zone carry distinctive genomic signatures potentially linked to post-earthquake recolonization processes. Specifically, our analysis identifies two parapatric spatial-genomic sectors of D. antarctica at Turakirae Head, which experienced the most dramatic uplift. Based on phylogeographical modelling, we infer that bull kelp in the Wellington region was probably a source for recolonization of the heavily uplifted Turakirae Head coastline, via two parallel, eastward recolonization events. By identifying multiple parapatric genotypic sectors within a recently recolonized coastal region, the current study provides support for the hypothesis that competing lineage expansions can generate striking spatial structuring of genetic diversity, even in highly dispersive taxa.

Allopatric origin of sympatric whitefish morphs with insights on the genetic basis of their reproductive isolation

The European whitefish (Coregonus lavaretus) species complex is a classic example of recent adaptive radiation. Here, we examine a whitefish population introduced to northern Finnish Lake Tsahkal in the late 1960s, where three divergent morphs (viz. littoral, pelagic, and profundal feeders) were found 10 generations after. Using demographic modeling based on genomic data, we show that whitefish morphs evolved during a phase of strict isolation, refuting a rapid sympatric divergence scenario. The lake is now an artificial hybrid zone between morphs originated in allopatry. Despite their current syntopy, clear genetic differentiation remains between two of the three morphs. Using admixture mapping, we identify five SNPs associated with gonad weight variation, a proxy for sexual maturity and spawning time. We suggest that ecological adaptations in spawning time evolved in allopatry are currently maintaining partial reproductive isolation in the absence of other barriers to gene flow.

Chromosomal inversions can limit adaptation to new environments

Chromosomal inversions are often thought to facilitate local adaptation and population divergence because they can link multiple adaptive alleles into non-recombining genomic blocks. Selection should thus be more efficient in driving inversion-linked adaptive alleles to high frequency in a population, particularly in the face of maladaptive gene flow. But what if ecological conditions and hence selection on inversion-linked alleles change? Reduced recombination within inversions could then constrain the formation of optimal combinations of pre-existing alleles under these new ecological conditions. Here, we outline this idea of inversions limiting adaptation and divergence when ecological conditions change across time or space. We reason and use simulations to illustrate that the benefit of inversions for local adaptation and divergence under one set of ecological conditions can come with a concomitant constraint for adaptation to novel sets of ecological conditions. This limitation of inversions to adaptation may contribute to the maintenance of polymorphism within species.

Freshwater Colonization, Adaptation, and Genomic Divergence in Threespine Stickleback

The Threespine Stickleback is ancestrally a marine fish, but many marine populations breed in fresh water (i.e., are anadromous), facilitating their colonization of isolated freshwater habitats a few years after they form. Repeated adaptation to fresh water during at least 10 My and continuing today has led to Threespine Stickleback becoming a premier system to study rapid adaptation. Anadromous and freshwater stickleback breed in sympatry and may hybridize, resulting in introgression of freshwater-adaptive alleles into anadromous populations, where they are maintained at low frequencies as ancient standing genetic variation. Anadromous stickleback have accumulated hundreds of freshwater-adaptive alleles that are disbursed as few loci per marine individual and provide the basis for adaptation when they colonize fresh water. Recent whole-lake experiments in lakes around Cook Inlet, Alaska have revealed how astonishingly rapid and repeatable this process is, with the frequency of 40% of the identified freshwater-adaptive alleles increasing from negligible (∼1%) in the marine founder to ≥50% within ten generations in fresh water, and freshwater phenotypes evolving accordingly. These high rates of genomic and phenotypic evolution imply very intense directional selection on phenotypes of heterozygotes. Sexual recombination rapidly assembles freshwater-adaptive alleles that originated in different founders into multilocus freshwater haplotypes, and regions important for adaptation to freshwater have suppressed recombination that keeps advantageous alleles linked within large haploblocks. These large haploblocks are also older and appear to have accumulated linked advantageous mutations. The contemporary evolution of Threespine Stickleback has provided broadly applicable insights into the mechanisms that facilitate rapid adaptation.

The ectodysplasin-A receptor is a candidate gene for lateral plate number variation in stickleback fish

Variation in lateral plating in stickleback fish represents a classical example of rapid and parallel adaptation in morphology. The underlying genetic architecture involves polymorphism at the ectodysplasin-A gene (EDA). However, lateral plate number is influenced by additional loci that remain poorly characterized. Here, we search for such loci by performing genome-wide differentiation mapping based on pooled whole-genome sequence data from a European stickleback population variable in the extent of lateral plating, while tightly controlling for the phenotypic effect of EDA. This suggests a new candidate locus, the EDA receptor gene (EDAR), for which additional support is obtained by individual-level targeted Sanger sequencing and by comparing allele frequencies among natural populations. Overall, our study illustrates the power of pooled whole-genome sequencing for searching phenotypically relevant loci and opens opportunities for exploring the population genetics and ecological significance of a new candidate locus for stickleback armor evolution.

Advancing human disease research with fish evolutionary mutant models

Model organism research is essential to understand disease mechanisms. However, laboratory-induced genetic models can lack genetic variation and often fail to mimic the spectrum of disease severity. Evolutionary mutant models (EMMs) are species with evolved phenotypes that mimic human disease. EMMs complement traditional laboratory models by providing unique avenues to study gene-by-environment interactions, modular mutations in noncoding regions, and their evolved compensations. EMMs have improved our understanding of complex diseases, including cancer, diabetes, and aging, and illuminated mechanisms in many organs. Rapid advancements of sequencing and genome-editing technologies have catapulted the utility of EMMs, particularly in fish. Fish are the most diverse group of vertebrates, exhibiting a kaleidoscope of specialized phenotypes, many that would be pathogenic in humans but are adaptive in the species' specialized habitat. Importantly, evolved compensations can suggest avenues for novel disease therapies. This review summarizes current research using fish EMMs to advance our understanding of human disease.

Habitat fragmentation induces rapid divergence of migratory and isolated sticklebacks

The adaptive capacity of many organisms is seriously challenged by human-imposed environmental change, which currently happens at unprecedented rates and magnitudes. For migratory fish, habitat fragmentation is a major challenge that can compromise their survival and reproduction. Therefore, it is important to study if fish populations can adapt to such modifications of their habitat. Here, we study whether originally anadromous three-spined stickleback populations (Gasterosteus aculeatus; “migrants”) changed in behavior and morphology in response to human-induced isolation. We made use of a natural field-experiment, where the construction of pumping stations and sluices in the 1970s unintendedly created replicates of land-locked stickleback populations (“resident”) in the Netherlands. For two years, we systematically tested populations of residents and migrants for differences in morphology and behavioral traits (activity, aggressiveness, exploration, boldness, and shoaling) in lab-based assays. We detected differences between migrant and resident populations in virtually all phenotypic traits studied: compared with the ancestral migrants, residents were smaller in size, had fewer and smaller plates and were significantly more active, aggressive, exploratory and bolder, and shoaled less. Despite large ecological differences between 2018 and 2019, results were largely consistent across the two years. Our study shows that human-induced environmental change has led to the rapid and consistent morphological and behavioral divergence of stickleback populations in about 50 generations. Such changes may be adaptive but this remains to be tested.

The maintenance of standing genetic variation: Gene flow vs. selective neutrality in Atlantic stickleback fish

Adaptation to derived habitats often occurs from standing genetic variation. The maintenance within ancestral populations of genetic variants favourable in derived habitats is commonly ascribed to long‐term antagonism between purifying selection and gene flow resulting from hybridization across habitats. A largely unexplored alternative idea based on quantitative genetic models of polygenic adaptation is that variants favoured in derived habitats are neutral in ancestral populations when their frequency is relatively low. To explore the latter, we first identify genetic variants important to the adaptation of threespine stickleback fish (Gasterosteus aculeatus) to a rare derived habitat—nutrient‐depleted acidic lakes—based on whole‐genome sequence data. Sequencing marine stickleback from six locations across the Atlantic Ocean then allows us to infer that the frequency of these derived variants in the ancestral habitat is unrelated to the likely opportunity for gene flow of these variants from acidic‐adapted populations. This result is consistent with the selective neutrality of derived variants within the ancestor. Our study thus supports an underappreciated explanation for the maintenance of standing genetic variation, and calls for a better understanding of the fitness consequences of adaptive variation across habitats and genomic backgrounds.

Contrasted patterns in climate change risk for Mediterranean fisheries

Climate change is rapidly becoming one of the biggest threats to marine life, and its impacts have the potential to strongly affect fisheries upon which millions of people rely. This is particularly crucial for the Mediterranean Sea, which is one of the world's biodiversity hotspots, one of the world's most overfished regions, and where temperatures are rising 25% faster than in the rest of the ocean on average. In this study, we calculated a vulnerability index for 100 species that compose 95% of the Mediterranean catches, through a trait-based approach. The Climate Risk Assessment (CRA) methodology was subsequently used to assess the risks due to climate change of Mediterranean fisheries. We found that the northern Mediterranean fisheries target more vulnerable species than their southern counterparts. However, when combining this catch-based vulnerability with a suite of socio-economic parameters, north African countries stand out as the most vulnerable to climate change impacts. Indeed, considering countries' exposure of the fisheries sector and their vulnerability to climate change, a sharp contrast between northern and southern Mediterranean appears, with Egypt and Tunisia scoring the highest risk. By integrating a trait-based approach on targeted marine species with socio-economic features, our analysis helps to better understand the ramifications of climate change consequences on Mediterranean fisheries and highlights the regions that could potentially be particularly affected.

Threespine Stickleback: A Model System For Evolutionary Genomics

The repeated adaptation of oceanic threespine sticklebacks to fresh water has made it a premier organism to study parallel evolution. These small fish have multiple distinct ecotypes that display a wide range of diverse phenotypic traits. Ecotypes are easily crossed in the laboratory, and families are large and develop quickly enough for quantitative trait locus analyses, positioning the threespine stickleback as a versatile model organism to address a wide range of biological questions. Extensive genomic resources, including linkage maps, a high-quality reference genome, and developmental genetics tools have led to insights into the genomic basis of adaptation and the identification of genomic changes controlling traits in vertebrates. Recently, threespine sticklebacks have been used as a model system to identify the genomic basis of highly complex traits, such as behavior and host-microbiome and host-parasite interactions. We review the latest findings and new avenues of research that have led the threespine stickleback to be considered a supermodel of evolutionary genomics.

Rapid Parallel Adaptation to Anthropogenic Heavy Metal Pollution

The impact of human-mediated environmental change on the evolutionary trajectories of wild organisms is poorly understood. In particular, capacity of species to adapt rapidly (in hundreds of generations or less), reproducibly and predictably to extreme environmental change is unclear. Silene uniflora is predominantly a coastal species, but it has also colonized isolated, disused mines with phytotoxic, zinc-contaminated soils. To test whether rapid, parallel adaptation to anthropogenic pollution has taken place, we used reduced representation sequencing (ddRAD) to reconstruct the evolutionary history of geographically proximate mine and coastal population pairs and found largely independent colonization of mines from different coastal sites. Furthermore, our results show that parallel evolution of zinc tolerance has occurred without gene flow spreading adaptive alleles between mine populations. In genomic regions where signatures of selection were detected across multiple mine-coast pairs, we identified genes with functions linked to physiological differences between the putative ecotypes, although genetic differentiation at specific loci is only partially shared between mine populations. Our results are consistent with a complex, polygenic genetic architecture underpinning rapid adaptation. This shows that even under a scenario of strong selection and rapid adaptation, evolutionary responses to human activities (and other environmental challenges) may be idiosyncratic at the genetic level and, therefore, difficult to predict from genomic data.

Physiological insight into the evolution of complex phenotypes: aerobic performance and the O2 transport pathway of vertebrates

Evolutionary physiology strives to understand how the function and integration of physiological systems influence the way in which organisms evolve. Studies of the O2 transport pathway - the integrated physiological system that transports O2 from the environment to mitochondria - are well suited to this endeavour. We consider the mechanistic underpinnings across the O2 pathway for the evolution of aerobic capacity, focusing on studies of artificial selection and naturally selected divergence among wild populations of mammals and fish. We show that evolved changes in aerobic capacity do not require concerted changes across the O2 pathway and can arise quickly from changes in one or a subset of pathway steps. Population divergence in aerobic capacity can be associated with the evolution of plasticity in response to environmental variation or activity. In some cases, initial evolutionary divergence of aerobic capacity arose exclusively from increased capacities for O2 diffusion and/or utilization in active O2-consuming tissues (muscle), which may often constitute first steps in adaptation. However, continued selection leading to greater divergence in aerobic capacity is often associated with increased capacities for circulatory and pulmonary O2 transport. Increases in tissue O2 diffusing capacity may augment the adaptive benefit of increasing circulatory O2 transport owing to their interactive influence on tissue O2 extraction. Theoretical modelling of the O2 pathway suggests that O2 pathway steps with a disproportionately large influence over aerobic capacity have been more likely to evolve, but more work is needed to appreciate the extent to which such physiological principles can predict evolutionary outcomes.

Genomic Architecture of Rapid Parallel Adaptation to Fresh Water in a Wild Fish

Rapid adaptation to novel environments may drive changes in genomic regions through natural selection. However, the genetic architecture underlying these adaptive changes is still poorly understood. Using population genomic approaches, we investigated the genomic architecture that underlies rapid parallel adaptation of Coilia nasus to fresh water by comparing four freshwater-resident populations with their ancestral anadromous population. Linkage disequilibrium network analysis and population genetic analyses revealed two putative large chromosome inversions on LG6 and LG22, which were enriched for outlier loci and exhibited parallel association with freshwater adaptation. Drastic frequency shifts and elevated genetic differentiation were observed for the two chromosome inversions among populations, suggesting that both inversions would undergo divergent selection between anadromous and resident ecotypes. Enrichment analysis of genes within chromosome inversions showed significant enrichment of genes involved in metabolic process, immunoregulation, growth, maturation, osmoregulation, and so forth, which probably underlay differences in morphology, physiology and behavior between the anadromous and freshwater-resident forms. The availability of beneficial standing genetic variation, large optimum shift between marine and freshwater habitats, and high efficiency of selection with large population size could lead to the observed rapid parallel adaptive genomic change. We propose that chromosomal inversions might have played an important role during the evolution of rapid parallel ecological divergence in the face of environmental heterogeneity in C. nasus. Our study provides insights into the genomic basis of rapid adaptation of complex traits in novel habitats and highlights the importance of structural genomic variants in analyses of ecological adaptation.

Reinventing the wheel? Reassessing the roles of gene flow, sorting and convergence in repeated evolution

Biologists have long been intrigued by apparently predictable and repetitive evolutionary trajectories inferred across a variety of lineages and systems. In recent years, high-throughput sequencing analyses have started to transform our understanding of such repetitive shifts. While researchers have traditionally categorized such shifts as either "convergent" or "parallel," based on relatedness of the lineages involved, emerging genomic insights provide an opportunity to better describe the actual evolutionary mechanisms at play. A synthesis of recent genomic analyses confirms that convergence is the predominant driver of repetitive evolution among species, whereas repeated sorting of standing variation is the major driver of repeated shifts within species. However, emerging data reveal numerous notable exceptions to these expectations, with recent examples of de novo mutations underpinning convergent shifts among even very closely related lineages, while repetitive sorting processes have occurred among even deeply divergent taxa, sometimes via introgression. A number of very recent analyses have found evidence for both processes occurring on different scales within taxa. We suggest that the relative importance of convergent versus sorting processes depends on the interplay between gene flow among populations, and phylogenetic relatedness of the lineages involved.

Predicting future from past: The genomic basis of recurrent and rapid stickleback evolution

Similar forms often evolve repeatedly in nature, raising long-standing questions about the underlying mechanisms. Here, we use repeated evolution in stickleback to identify a large set of genomic loci that change recurrently during colonization of freshwater habitats by marine fish. The same loci used repeatedly in extant populations also show rapid allele frequency changes when new freshwater populations are experimentally established from marine ancestors. Marked genotypic and phenotypic changes arise within 5 years, facilitated by standing genetic variation and linkage between adaptive regions. Both the speed and location of changes can be predicted using empirical observations of recurrence in natural populations or fundamental genomic features like allelic age, recombination rates, density of divergent loci, and overlap with mapped traits. A composite model trained on these stickleback features can also predict the location of key evolutionary loci in Darwin's finches, suggesting that similar features are important for evolution across diverse taxa.

Repeatability of Adaptive Radiation Depends on Spatial Scale: Regional Versus Global Replicates of Stickleback in Lake Versus Stream Habitats

The repeatability of adaptive radiation is expected to be scale-dependent, with determinism decreasing as greater spatial separation among "replicates" leads to their increased genetic and ecological independence. Threespine stickleback (Gasterosteus aculeatus) provide an opportunity to test whether this expectation holds for the early stages of adaptive radiation-their diversification in freshwater ecosystems has been replicated many times. To better understand the repeatability of that adaptive radiation, we examined the influence of geographic scale on levels of parallel evolution by quantifying phenotypic and genetic divergence between lake and stream stickleback pairs sampled at regional (Vancouver Island) and global (North America and Europe) scales. We measured phenotypes known to show lake-stream divergence and used reduced representation genome-wide sequencing to estimate genetic divergence. We assessed the scale dependence of parallel evolution by comparing effect sizes from multivariate models and also the direction and magnitude of lake-stream divergence vectors. At the phenotypic level, parallelism was greater at the regional than the global scale. At the genetic level, putative selected loci showed greater lake-stream parallelism at the regional than the global scale. Generally, the level of parallel evolution was low at both scales, except for some key univariate traits. Divergence vectors were often orthogonal, highlighting possible ecological and genetic constraints on parallel evolution at both scales. Overall, our results confirm that the repeatability of adaptive radiation decreases at increasing spatial scales. We suggest that greater environmental heterogeneity at larger scales imposes different selection regimes, thus generating lower repeatability of adaptive radiation at larger spatial scales.

Ancient and modern stickleback genomes reveal the demographic constraints on adaptation

Adaptation is typically studied by comparing modern populations with contrasting environments. Individuals persisting in the ancestral habitat are typically used to represent the ancestral founding population; however, it has been questioned whether these individuals are good proxies for the actual ancestors.¹ To address this, we applied a paleogenomics approach² to directly access the ancestral genepool: partially sequencing the genomes of two 11- to 13,000-year-old stickleback recovered from the transitionary layer between marine and freshwater sediments of two Norwegian isolation lakes³ and comparing them with 30 modern stickleback genomes from the same lakes and adjacent marine fjord, in addition to a global dataset of 20 genomes.⁴ The ancient stickleback shared genome-wide ancestry with the modern fjord population, whereas modern lake populations have lost substantial ancestral variation following founder effects, and subsequent drift and selection. Freshwater-adaptive alleles found in one ancient stickleback genome have not risen to high frequency in the present-day population from the same lake. Comparison to the global dataset suggested incomplete adaptation to freshwater in our modern lake populations. Our findings reveal the impact of population bottlenecks in constraining adaptation due to reduced efficacy of selection on standing variation present in founder populations.

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

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

Heritability of DNA methylation in threespine stickleback (Gasterosteus aculeatus)

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.

A framework for understanding gene expression plasticity and its influence on stress tolerance

Phenotypic plasticity can serve as a stepping stone towards adaptation. Recently, studies have shown that gene expression contributes to emergent stress responses such as thermal tolerance, with tolerant and susceptible populations showing distinct transcriptional profiles. However, given the dynamic nature of gene expression, interpreting transcriptomic results in a way that elucidates the functional connection between gene expression and the observed stress response is challenging. Here, we present a conceptual framework to guide interpretation of gene expression reaction norms in the context of stress tolerance. We consider the evolutionary and adaptive potential of gene expression reaction norms and discuss the influence of sampling timing, transcriptomic resilience, as well as complexities related to life history when interpreting gene expression dynamics and how these patterns relate to host tolerance. We highlight corals as a case study to demonstrate the value of this framework for non-model systems. As species face rapidly changing environmental conditions, modulating gene expression can serve as a mechanistic link from genetic and cellular processes to the physiological responses that allow organisms to thrive under novel conditions. Interpreting how or whether a species can employ gene expression plasticity to ensure short-term survival will be critical for understanding the global impacts of climate change across diverse taxa.

Environmental correlates of genetic variation in the invasive European starling in North America

Populations of invasive species that colonize and spread in novel environments may differentiate both through demographic processes and local selection. European starlings (Sturnus vulgaris) were introduced to New York in 1890 and subsequently spread throughout North America, becoming one of the most widespread and numerous bird species on the continent. Genome-wide comparisons across starling individuals and populations can identify demographic and/or selective factors that facilitated this rapid and successful expansion. We investigated patterns of genomic diversity and differentiation using reduced-representation genome sequencing of 17 winter-season sampling sites. Consistent with this species' high dispersal rate and rapid expansion history, we found low geographical differentiation and few F_ST outliers even at a continental scale. Despite starting from a founding population of ~180 individuals, North American starlings show only a moderate genetic bottleneck, and models suggest a dramatic increase in effective population size since introduction. In genotype-environment associations we found that ~200 single-nucleotide polymorphisms are correlated with temperature and/or precipitation against a background of negligible genome- and range-wide divergence. Given this evidence, we suggest that local adaptation in North American starlings may have evolved rapidly even in this wide-ranging and evolutionarily young system. This survey of genomic signatures of expansion in North American starlings is the most comprehensive to date and complements ongoing studies of world-wide local adaptation in these highly dispersive and invasive birds.

Genetic response to human-induced habitat changes in the marine environment: A century of evolution of European sprat in Landvikvannet, Norway

Habitat changes represent one of the five most pervasive threats to biodiversity. However, anthropogenic activities also have the capacity to create novel niche spaces to which species respond differently. In 1880, one such habitat alterations occurred in Landvikvannet, a freshwater lake on the Norwegian coast of Skagerrak, which became brackish after being artificially connected to the sea. This lake is now home to the European sprat, a pelagic marine fish that managed to develop a self-recruiting population in barely few decades. Landvikvannet sprat proved to be genetically isolated from the three main populations described for this species; that is, Norwegian fjords, Baltic Sea, and the combination of North Sea, Kattegat, and Skagerrak. This distinctness was depicted by an accuracy self-assignment of 89% and a highly significant F ST between the lake sprat and each of the remaining samples (average of ≈0.105). The correlation between genetic and environmental variation indicated that salinity could be an important environmental driver of selection (3.3% of the 91 SNPs showed strong associations). Likewise, Isolation by Environment was detected for salinity, although not for temperature, in samples not adhering to an Isolation by Distance pattern. Neighbor-joining tree analysis suggested that the source of the lake sprat is in the Norwegian fjords, rather than in the Baltic Sea despite a similar salinity profile. Strongly drifted allele frequencies and lower genetic diversity in Landvikvannet compared with the Norwegian fjords concur with a founder effect potentially associated with local adaptation to low salinity. Genetic differentiation (F ST) between marine and brackish sprat is larger in the comparison Norway-Landvikvannet than in Norway-Baltic, which suggests that the observed divergence was achieved in Landvikvannet in some 65 generations, that is, 132 years, rather than gradually over thousands of years (the age of the Baltic Sea), thus highlighting the pace at which human-driven evolution can happen.

Fitness maps to a large-effect locus in introduced stickleback populations

Mutations of small effect underlie most adaptation to new environments, but beneficial variants with large fitness effects are expected to contribute under certain conditions. Genes and genomic regions having large effects on phenotypic differences between populations are known from numerous taxa, but fitness effect sizes have rarely been estimated. We mapped fitness over a generation in an F2 intercross between a marine and a lake stickleback population introduced to a freshwater pond. A quantitative trait locus map of the number of surviving offspring per F2 female detected a single, large-effect locus near Ectodysplasin (Eda), a gene having an ancient freshwater allele causing reduced bony armor and other changes. F2 females homozygous for the freshwater allele had twice the number of surviving offspring as homozygotes for the marine allele, producing a large selection coefficient, s = 0.50 ± 0.09 SE. Correspondingly, the frequency of the freshwater allele increased from 0.50 in F2 mothers to 0.58 in surviving offspring. We compare these results to allele frequency changes at the Eda gene in an Alaskan lake population colonized by marine stickleback in the 1980s. The frequency of the freshwater Eda allele rose steadily over multiple generations and reached 95% within 20 y, yielding a similar estimate of selection, s = 0.49 ± 0.05, but a different degree of dominance. These findings are consistent with other studies suggesting strong selection on this gene (and/or linked genes) in fresh water. Selection on ancient genetic variants carried by colonizing ancestors is likely to increase the prevalence of large-effect fitness variants in adaptive evolution.

Insects in high-elevation streams: Life in extreme environments imperiled by climate change

Climate change is altering conditions in high-elevation streams worldwide, with largely unknown effects on resident communities of aquatic insects. Here, we review the challenges of climate change for high-elevation aquatic insects and how they may respond, focusing on current gaps in knowledge. Understanding current effects and predicting future impacts will depend on progress in three areas. First, we need better descriptions of the multivariate physical challenges and interactions among challenges in high-elevation streams, which include low but rising temperatures, low oxygen supply and increasing oxygen demand, high and rising exposure to ultraviolet radiation, low ionic strength, and variable but shifting flow regimes. These factors are often studied in isolation even though they covary in nature and interact in space and time. Second, we need a better mechanistic understanding of how physical conditions in streams drive the performance of individual insects. Environment-performance links are mediated by physiology and behavior, which are poorly known in high-elevation taxa. Third, we need to define the scope and importance of potential responses across levels of biological organization. Short-term responses are defined by the tolerances of individuals, their capacities to perform adequately across a range of conditions, and behaviors used to exploit local, fine-scale variation in abiotic factors. Longer term responses to climate change, however, may include individual plasticity and evolution of populations. Whether high-elevation aquatic insects can mitigate climatic risks via these pathways is largely unknown.

The relationship between morphology and behavior in mixed-species flocks of island birds

Understanding how co‐occurring species divide ecological space is a central issue in ecology. Functional traits have the potential to serve as a means for quantitatively assessing niche partitioning by different species based on their ecological attributes, such as morphology, behavior, or trophic habit. This enables testing ecological and evolutionary questions using functional traits at spatio‐temporal scales that are not feasible using traditional field methods. Both rapid evolutionary change and inter‐ and intraspecific competition, however, may limit the utility of morphological functional traits as indicators of how niches are partitioned. To address how behavior and morphology interact, we quantified foraging behavior of mixed‐species flocks of birds in the Solomon Islands to test whether behavior and morphology are correlated in these flocks. We find that foraging behavior is significantly correlated with morphological traits (p = .05), but this correlation breaks down after correcting for phylogenetic relatedness (p = .66). These results suggest that there are consistent correlations between aspects of behavior and morphology at large taxonomic scales (e.g., across genera), but the relationship between behavior and morphology depends largely on among‐clade differences and may be idiosyncratic at shallower scales (e.g., within genera). As a result, general relationships between behaviors and morphology may not be applicable when comparing close relatives.

The tree and the table: Darwin, Mendeleev and the meaning of 'theory'

Darwin and Mendeleev revolutionized their respective disciplines by organizing diverse facts into simple, pictorial representations-a tree and a table. Each representation provides a foundation for a scientific theory for two reasons. First, a successful theory unites diverse phenomena under a single explanatory framework. Second, it does so in a way that defines paths for future inquiry that extends its reach and tests its limits. For Mendeleev, this meant creating a table that accommodated the current understanding of the elements but also contained blanks that predicted the discovery of previously unknown elements. More importantly, the structure of the table helped shape future research to define the structure of matter. For Darwin, envisioning life as a tree meant defining the rules that govern the origin of adaptations, species and shape the constantly shifting diversity of life. At the same time, his theory inspired research into the laws of inheritance and created diverse new areas of research, like behaviour, sexual selection and biogeography. The shared property of Darwin and Mendeleev's contributions was to provide a unifying rational explanation for natural phenomena. This article is part of the theme issue 'Mendeleev and the periodic table'.

Post-glacial establishment of locally adapted fish populations over a steep salinity gradient

Studies of colonization of new habitats that appear from rapidly changing environments are interesting and highly relevant to our understanding of divergence and speciation. Here, we analyse phenotypic and genetic variation involved in the successful establishment of a marine fish (sand goby, Pomatoschistus minutus) over a steep salinity drop from 35 PSU in the North Sea (NE Atlantic) to two PSU in the inner parts of the post-glacial Baltic Sea. We first show that populations are adapted to local salinity in a key reproductive trait, the proportion of motile sperm. Thereafter, we show that genome variation at 22,190 single nucleotide polymorphisms (SNPs) shows strong differentiation among populations along the gradient. Sequences containing outlier SNPs and transcriptome sequences, mapped to a draft genome, reveal associations with genes with relevant functions for adaptation in this environment but without overall evidence of functional enrichment. The many contigs involved suggest polygenic differentiation. We trace the origin of this differentiation using demographic modelling and find the most likely scenario is that at least part of the genetic differentiation is older than the Baltic Sea and is a result of isolation of two lineages prior to the current contact over the North Sea-Baltic Sea transition zone.

Estimating the potential for coral adaptation to global warming across the Indo-West Pacific

The potential of reef-building corals to adapt to increasing sea-surface temperatures is often debated but has rarely been comprehensively modeled on a region-wide scale. We used individual-based simulations to model adaptation to warming in a coral metapopulation comprising 680 reefs and representing the whole of the Central Indo-West Pacific. Encouragingly, some reefs-most notably Vietnam, Japan, Taiwan, New Caledonia and the southern half of the Great Barrier Reef-exhibited high capacity for adaptation and, in our model, maintained coral cover even under a rapid "business-as-usual" warming scenario throughout the modeled period (200 years). Higher resilience of these reefs was observed under all tested parameter settings except the models prohibiting selection and/or migration during warming. At the same time, the majority of reefs in the region tended to collapse within the first 100 years of warming. The adaptive potential (odds of maintaining high coral cover) of a given reef could be predicted based on two metrics: the reef's present-day temperature, and the proportion of recruits immigrating from warmer locations. The latter metric explains the most variation in adaptive potential, and significantly correlates with actual coral cover changes observed throughout the region between the 1970s and the early 2000s. These findings will help prioritize coral conservation efforts and plan assisted gene flow interventions to boost the adaptive potential of specific coral populations.

Genomic release-recapture experiment in the wild reveals within-generation polygenic selection in stickleback fish

How rapidly natural selection sorts genome-wide standing genetic variation during adaptation remains largely unstudied experimentally. Here, we present a genomic release-recapture experiment using paired threespine stickleback fish populations adapted to selectively different lake and stream habitats. First, we use pooled whole-genome sequence data from the original populations to identify hundreds of candidate genome regions likely under divergent selection between these habitats. Next, we generate F2 hybrids from the same lake-stream population pair in the laboratory and release thousands of juveniles into a natural stream habitat. Comparing the individuals surviving one year of stream selection to a reference sample of F2 hybrids allows us to detect frequency shifts across the candidate regions toward the genetic variants typical of the stream population-an experimental outcome consistent with polygenic directional selection. Our study reveals that adaptation in nature can be detected as a genome-wide signal over just a single generation.

Architecture of Parallel Adaptation in Ten Lacustrine Threespine Stickleback Populations from the White Sea Area

Adaptation of threespine stickleback to freshwater involves parallel recruitment of freshwater alleles in clusters of closely linked sites, or divergence islands (DIs). However, it remains unclear to what extent the DIs and the alleles that constitute them coincide between populations that underwent adaptation to freshwater independently. We examine threespine sticklebacks from ten freshwater lakes that emerged 500–1500 years ago in the White Sea basin, with the emphasis on repeatability of genomic patterns of adaptation among the lake populations and the role of local recombination rate in the distribution and structure of DIs. The 65 detected DIs are clustered in the genome, forming 12 aggregations, and this clustering cannot be explained by the variation of the recombination rate. Only 21 of the DIs are present in all the freshwater populations, likely being indispensable for successful colonization of freshwater environment by the ancestral marine population. Within most DIs, the same set of single nucleotide polymorphisms (SNPs) distinguish marine and freshwater haplotypes in all the lake populations; however, in some DIs, freshwater alleles differ between populations, suggesting that they could have been established by recruitment of different haplotypes in different populations.

Asymmetric Isolation and the Evolution of Behaviors Influencing Dispersal: Rheotaxis of Guppies above Waterfalls

Populations that are asymmetrically isolated, such as above waterfalls, can sometimes export emigrants in a direction from which they do not receive immigrants, and thus provide an excellent opportunity to study the evolution of dispersal traits. We investigated the rheotaxis of guppies above barrier waterfalls in the Aripo and Turure rivers in Trinidad-the later having been introduced in 1957 from a below-waterfall population in another drainage. We predicted that, as a result of strong selection against downstream emigration, both of these above-waterfall populations should show strong positive rheotaxis. Matching these expectations, both populations expressed high levels of positive rheotaxis, possibly reflecting contemporary (rapid) evolution in the introduced Turure population. However, the two populations used different behaviors to achieve the same performance of strong positive rheotaxis, as has been predicted in the case of multiple potential evolutionary solutions to the same functional challenge (i.e., "many-to-one mapping"). By contrast, we did not find any difference in rheotactic behavior above versus below waterfalls on a small scale within either river, suggesting constraints on adaptive divergence on such scales.

A Few Stickleback Suffice for the Transport of Alleles to New Lakes

Threespine stickleback populations provide a striking example of local adaptation to divergent habitats in populations that are connected by recurrent gene flow. These small fish occur in marine and freshwater habitats throughout the Northern Hemisphere, and in numerous cases the smaller freshwater populations have been established “de novo” from marine colonists. Independently evolved freshwater populations exhibit similar phenotypes that have been shown to derive largely from the same standing genetic variants. Geographic isolation prevents direct migration between the freshwater populations, strongly suggesting that these shared locally adaptive alleles are transported through the marine population. However it is still largely unknown how gene flow, recombination, and selection jointly impact the standing variation that might fuel this adaptation. Here we use individual-based, spatially explicit simulations to determine the levels of gene flow that best match observed patterns of allele sharing among habitats in stickleback. We aim to better understand how gene flow and local adaptation in large metapopulations determine the speed of adaptation and re-use of standing genetic variation. In our simulations we find that repeated adaptation uses a shared set of alleles that are maintained at low frequency by migration-selection balance in oceanic populations. This process occurs over a realistic range of intermediate levels of gene flow that match previous empirical population genomic studies in stickleback. Examining these simulations more deeply reveals how lower levels of gene flow leads to slow, independent adaptation to different habitats, whereas higher levels of gene flow leads to significant mutation load – but an increased probability of successful population genomic scans for locally adapted alleles. Surprisingly, we find that the genealogical origins of most freshwater adapted alleles can be traced back to the original generation of marine individuals that colonized the lakes, as opposed to subsequent migrants. These simulations provide deeper context for existing studies of stickleback evolutionary genomics, and guidance for future empirical studies in this model. More broadly, our results support existing theory of local adaptation but extend it by more completely documenting the genealogical history of adaptive alleles in a metapopulation.

Evidence of a landlocked reproducing population of the marine pejerrey Odontesthes argentinensis (Actinopterygii; Atherinopsidae)

In South America, the order Atheriniformes includes the monophyletic genus Odontesthes with 20 species that inhabit freshwater, estuarine and coastal environments. Pejerrey Odontesthes argentinensis is widely distributed in coastal and estuarine areas of the Atlantic Ocean and is known to foray into estuaries of river systems, particularly in conditions of elevated salinity. However, to our knowledge, a landlocked self-sustaining population has never been recorded. In this study, we examined the pejerrey population of Salada de Pedro Luro Lake (south-east of Buenos Aires Province, Argentina) to clarify its taxonomic identity. An integrative taxonomic analysis based on traditional meristic, landmark-based morphometrics and genetic techniques suggests that the Salada de Pedro Luro pejerrey population represents a novel case of physiological and morphological adaptation of a marine pejerrey species to a landlocked environment and emphasises the environmental plasticity of this group of fishes.