Fine scale diversity in the lava: genetic and phenotypic diversity in small populations of Arctic charr Salvelinus alpinus

BACKGROUND

A major goal in evolutionary biology is to understand the processes underlying phenotypic variation in nature. Commonly, studies have focused on large interconnected populations or populations found along strong environmental gradients. However, studies on small fragmented populations can give strong insight into evolutionary processes in relation to discrete ecological factors. Evolution in small populations is believed to be dominated by stochastic processes, but recent work shows that small populations can also display adaptive phenotypic variation, through for example plasticity and rapid adaptive evolution. Such evolution takes place even though there are strong signs of historical bottlenecks and genetic drift. Here we studied 24 small populations of the freshwater fish Arctic charr (Salvelinus alpinus) found in groundwater filled lava caves. Those populations were found within a few km2-area with no apparent water connections between them. We studied the relative contribution of neutral versus non-neutral evolutionary processes in shaping phenotypic divergence, by contrasting patterns of phenotypic and neutral genetic divergence across populations in relation to environmental measurements. This allowed us to model the proportion of phenotypic variance explained by the environment, taking in to account the observed neutral genetic structure.

RESULTS

These populations originated from the nearby Lake Mývatn, and showed small population sizes with low genetic diversity. Phenotypic variation was mostly correlated with neutral genetic diversity with only a small environmental effect.

CONCLUSIONS

Phenotypic diversity in these cave populations appears to be largely the product of neutral processes, fitting the classical evolutionary expectations. However, the fact that neutral processes did not explain fully the phenotypic patterns suggests that further studies can increase our understanding on how neutral evolutionary processes can interact with other forces of selection at early stages of divergence. The accessibility of these populations has provided the opportunity for long-term monitoring of individual fish, allowing tracking how the environment can influence phenotypic and genetic divergence for shaping and maintaining diversity in small populations. Such studies are important, especially in freshwater, as habitat alteration is commonly breaking populations into smaller units, which may or may not be viable.

Genomic underpinnings of head and body shape in Arctic charr ecomorph pairs

Across its Holarctic range, Arctic charr (Salvelinus alpinus) populations have diverged into distinct trophic specialists across independent replicate lakes. The major aspect of divergence between ecomorphs is in head shape and body shape, which are ecomorphological traits reflecting niche use. However, whether the genomic underpinnings of these parallel divergences are consistent across replicates was unknown but key for resolving the substrate of parallel evolution. We investigated the genomic basis of head shape and body shape morphology across four benthivore-planktivore ecomorph pairs of Arctic charr in Scotland. Through genome-wide association analyses, we found genomic regions associated with head shape (89 SNPs) or body shape (180 SNPs) separately and 50 of these SNPs were strongly associated with both body and head shape morphology. For each trait separately, only a small number of SNPs were shared across all ecomorph pairs (3 SNPs for head shape and 10 SNPs for body shape). Signs of selection on the associated genomic regions varied across pairs, consistent with evolutionary demography differing considerably across lakes. Using a comprehensive database of salmonid QTLs newly augmented and mapped to a charr genome, we found several of the head- and body-shape-associated SNPs were within or near morphology QTLs from other salmonid species, reflecting a shared genetic basis for these phenotypes across species. Overall, our results demonstrate how parallel ecotype divergences can have both population-specific and deeply shared genomic underpinnings across replicates, influenced by differences in their environments and demographic histories.

The effects of host ecology and phylogeny on gut microbiota (non)parallelism across birds and mammals

IMPORTANCE

What are the roles of determinism and contingency in evolution? The paleontologist and evolutionary biologist Stephen J. Gould raised this question in his famous thought experiment of "replaying life's tape." Settings where independent lineages have repeatedly adapted to similar ecological niches (i.e., parallel evolution) are well suited to address this question. Here, we quantified whether repeated ecological shifts across 53 mammalian and 50 avian host species are associated with parallel gut microbiota changes. Our results indicate that parallel shifts in host diet are associated with greater gut microbiota parallelism (i.e., more deterministic). While further research will be necessary to obtain a comprehensive picture of the circumstances under which deterministic gut microbiota changes might be expected, our study can be instrumental in motivating the use of more quantitative methods in microbiota research. This, in turn, can help us better understand microbiota dynamics during adaptive evolution of their hosts.

Patterns in antipredator armature reduction and maintenance in isolated spring populations of an amphipod crustacean

Organisms colonizing new habitats can undergo adaptive change due to novel selective landscapes encountered in the new environment. Examples in nature where the development of the same traits has repeatedly occurred on multiple independent occasions upon colonizing a novel habitat represent instances of parallel evolution. Here we test whether the colonization of spring habitat by the principally lacustrine amphipod crustacean Pallaseopsis quadrispinosa has resulted in parallel evolution in armature traits using empirical data on morphology and mitochondrial DNA and through a breeding experiment. Analysis of mtDNA CO1 sequences shows that the spring populations share no common history and have evolved in isolation from each other and from their neighbouring lake populations since deglaciation approximately 12,000 years ago and are now fixed for different haplogroups. Dorsal spines and lateral projections were absent or less developed in all spring populations than in lake populations. Variation in armature development also could be explained by predator presence as populations with fish predators exhibited more developed spines than those without fish. In a laboratory breeding experiment, hybrid Spring × Lake F1 offspring had intermediate development of armature compared to offspring of Lake × Lake and Spring × Spring matings. The results support the hypothesis that armature reduction has independently evolved on multiple occasions in P. quadrispinosa. Recent research has questioned the degree to which parallel evolution actually explains variance in traits. Taking into account the predation regime, sexual dimorphism and mineral composition of the trait, a more precise understanding of the factors influencing parallel evolution emerges.

Differential genotype response to increased resource abundance helps explain parallel evolution of Daphnia populations in the wild

Under controlled laboratory conditions, previous studies have shown that selection can produce repeatable evolutionary trajectories. Yet, the question remains for many of these studies if, given identical starting populations, evolution in the wild proceeds in a non-random direction. In the present study, we investigated the extent to which rapid evolution in the wild is parallel by monitoring the genetic composition of replicate populations of Daphnia in field mesocosms containing two clonal genotypes. We found parallel changes across all nine mesocosms, in which the same genotype increased in frequency. To probe whether genotype-specific response to resource abundance could have led to this frequency change, we conducted a life-history assay under high-resource abundance and low-resource abundance. We found that resource exploitation differed by genotype, in that, while one genotype (the winner in the field mesocosms) was more fit than the other genotype at high resources, the other genotype performed slightly better at low resources. We suspect that levels of resource abundance found in the summer field mesocosms had values in which the genotype better with abundant resources had the advantage. These findings suggest that variation in certain traits associated with resource acquisition can drive genotype frequency change.

Genomic responses to parallel temperature gradients in the eelgrass Zostera marina in adjacent bays

The extent of parallel genomic responses to similar selective pressures depends on a complex array of environmental, demographic, and evolutionary forces. Laboratory experiments with replicated selective pressures yield mixed outcomes under controlled conditions and our understanding of genomic parallelism in the wild is limited to a few well-established systems. Here, we examine genomic signals of selection in the eelgrass Zostera marina across temperature gradients in adjacent embayments. Although we find many genomic regions with signals of selection within each bay there is very little overlap in signals of selection at the SNP level, despite most polymorphisms being shared across bays. We do find overlap at the gene level, potentially suggesting multiple mutational pathways to the same phenotype. Using polygenic models we find that some sets of candidate SNPs are able to predict temperature across both bays, suggesting that small but parallel shifts in allele frequencies may be missed by independent genome scans. Together, these results highlight the continuous rather than binary nature of parallel evolution in polygenic traits and the complexity of evolutionary predictability.

Repeated genetic divergence plays a minor role in repeated phenotypic divergence of lake-stream stickleback

Recent studies have shown that the repeated evolution of similar phenotypes in response to similar ecological conditions (here "parallel evolution") often occurs through mutations in the same genes. However, many previous studies have focused on known candidate genes in a limited number of systems. Thus, the question of how often parallel phenotypic evolution is due to parallel genetic changes remains open. Here, we used quantitative trait locus (QTL) mapping in F2 intercrosses between lake and stream threespine stickleback (Gasterosteus aculeatus) from four independent watersheds on Vancouver Island, Canada to determine whether the same QTL underlie divergence in the same phenotypes across, between, and within watersheds. We find few parallel QTL, even in independent crosses from the same watershed or for phenotypes that have diverged in parallel. These findings suggest that different mutations can lead to similar phenotypes. The low genetic repeatability observed in these lake-stream systems contrasts with the higher genetic repeatability observed in other stickleback systems. We speculate that differences in evolutionary history, gene flow, and/or the strength and direction of selection might explain these differences in genetic parallelism and emphasize that more work is needed to move beyond documenting genetic parallelism to identifying the underlying causes.

The role of plastic and evolved DNA methylation in parallel adaptation of threespine stickleback (Gasterosteus aculeatus)

Repeated phenotypic patterns among populations undergoing parallel evolution in similar environments provide support for the deterministic role of natural selection. Epigenetic modifications can mediate plastic and evolved phenotypic responses to environmental change and might make important contributions to parallel adaptation. While many studies have explored the genetic basis of repeated phenotypic divergence, the role of epigenetic processes during parallel adaptation remains unclear. The parallel evolution of freshwater ecotypes of threespine stickleback fish (Gasterosteus aculeatus) following colonization of thousands of lakes and streams from the ocean is a classic example of parallel phenotypic and genotypic adaptation. To investigate epigenetic modifications during parallel adaptation of threespine stickleback, we reanalysed three independent data sets that investigated DNA methylation variation between marine and freshwater ecotypes. Although we found widespread methylation differentiation between ecotypes, there was no significant tendency for CpG sites associated with repeated methylation differentiation across studies to be parallel versus nonparallel. To next investigate the role of plastic versus evolved changes in methylation during freshwater adaptation, we explored if CpG sites exhibiting methylation plasticity during salinity change were more likely to also show evolutionary divergence in methylation between ecotypes. The directions of divergence between ecotypes were generally in the opposite direction to those observed for plasticity when ecotypes were challenged with non-native salinity conditions, suggesting that most plastic responses are likely to be maladaptive during colonization of new environments. Finally, we found a greater number of CpG sites showing evolved changes when ancestral marine ecotypes are acclimated to freshwater environments, whereas plastic changes predominate when derived freshwater ecotypes transition back to their ancestral marine environments. These findings provide evidence for an epigenetic contribution to parallel adaptation and demonstrate the contrasting roles of plastic and evolved methylation differences during adaptation to new environments.

Compiling forty years of guppy research to investigate the factors contributing to (non)parallel evolution

Examples of parallel evolution have been crucial for our understanding of adaptation via natural selection. However, strong parallelism is not always observed even in seemingly similar environments where natural selection is expected to favour similar phenotypes. Leveraging this variation in parallelism within well-researched study systems can provide insight into the factors that contribute to variation in adaptive responses. Here we analyse the results of 36 studies reporting 446 average trait values in Trinidadian guppies, Poecilia reticulata, from different predation regimes. We examine how the extent of predator-driven phenotypic parallelism is influenced by six factors: sex, trait type, rearing environment, ecological complexity, evolutionary history, and time since colonization. Analyses show that parallel evolution in guppies is highly variable and weak on average, with only 24.7% of the variation among populations being explained by predation regime. Levels of parallelism appeared to be especially weak for colour traits, and parallelism decreased with increasing complexity of evolutionary history (i.e., when estimates of parallelism from populations within a single drainage were compared to estimates of parallelism from populations pooled between two major drainages). Suggestive - but not significant - trends that warrant further research include interactions between the sexes and different trait categories. Quantifying and accounting for these and other sources of variation among evolutionary 'replicates' can be leveraged to better understand the extent to which seemingly similar environments drive parallel and nonparallel aspects of phenotypic divergence.

Hybrid evolution repeats itself across environmental contexts in Texas sunflowers (Helianthus)

To what extent is evolution repeatable? Little is known about whether the evolution of hybrids is more (or less) repeatable than that of nonhybrids. We used field experimental evolution in annual sunflowers (Helianthus) in Texas to ask the extent to which hybrid evolution is repeatable across environments compared to nonhybrid controls. We created hybrids between Helianthus annuus (L.) and H. debilis (Nutt.) and grew plots of both hybrids and nonhybrid controls through eight generations at three sites in Texas. We collected seeds from each generation and grew each generation × treatment × home site combination at two final common gardens. We estimated the strength and direction of evolution in terms of fitness and 24 traits, tested for repeated versus nonrepeated evolution, and assessed overall phenotypic evolution across lineages and in relation to a locally adapted phenotype. Hybrids consistently evolved higher fitness over time, while controls did not, although trait evolution varied in strength across home sites. Repeated evolution was more evident in hybrids versus nonhybrid controls, and hybrid evolution was often in the direction of the locally adapted phenotype. Our findings have implications for both the nature of repeatability in evolution and the contribution of hybridization to evolution across environmental contexts.

Non-parallel morphological divergence following colonization of a new host plant

Adaptation to new ecological niches is known to spur population diversification and may lead to speciation if gene flow is ceased. While adaptation to the same ecological niche is expected to be parallel, it is more difficult to predict whether selection against maladaptive hybridization in secondary sympatry results in parallel divergence also in traits that are not directly related to the ecological niches. Such parallelisms in response to selection for reproductive isolation can be identified through estimating parallelism in reproductive character displacement across different zones of secondary contact. Here, we use a host shift in the phytophagous peacock fly Tephritis conura, with both host races represented in two geographically separate areas East and West of the Baltic Sea to investigate convergence in morphological adaptations. We asked (i) if there are consistent morphological adaptations to a host plant shift and (ii) if the response to secondary sympatry with the alternate host race is parallel across contact zones. We found surprisingly low and variable, albeit significant, divergence between host races. Only one trait, the length of the female ovipositor, which serves an important function in the interaction with the hosts, was consistently different between host races. Instead, co-existence with the other host race significantly affected the degree of morphological divergence, but the divergence was largely driven by different traits in different contact zones. Thus, local stochastic fixation or reinforcement could generate trait divergence, and additional evidence is needed to conclude whether divergence is locally adaptive.

Heterosis counteracts hybrid breakdown to forestall speciation by parallel natural selection

In contrast to ecological speciation, where reproductive isolation evolves as a consequence of divergent natural selection, speciation by parallel natural selection has been less thoroughly studied. To test whether parallel evolution drives speciation, we leveraged the repeated evolution of benthic and limnetic ecotypes of threespine stickleback fish and estimated fitness for pure crosses and within-ecotype hybrids in semi-natural ponds and in laboratory aquaria. In ponds, we detected hybrid breakdown in both ecotypes but this was counterbalanced by heterosis and the strength of post-zygotic isolation was nil. In aquaria, we detected heterosis in limnetic crosses and breakdown in benthic crosses, which is suggestive of process- and ecotype-specific environment-dependence. In ponds, heterosis and breakdown were three times greater in limnetic crosses than in benthic crosses, contrasting the prediction that the fitness consequences of hybridization should be greater in crosses among more derived ecotypes. Consistent with a primary role for stochastic processes, patterns differed among crosses between populations from different lakes. Yet, the observation of qualitatively similar patterns of heterosis and hybrid breakdown for both ecotypes when averaging the lake pairs indicates that the outcome of hybridization is repeatable in a general sense.

The genetics of adaptation in freshwater Eurasian shad (Alosa)

Studying the genetics of phenotypic convergence can yield important insights into adaptive evolution. Here, we conducted a comparative genomic study of four lineages (species and subspecies) of anadromous shad (Alosa) that have independently evolved life cycles entirely completed in freshwater. Three naturally diverged (A. fallax lacustris, A. f. killarnensis, and A. macedonica), and the fourth (A. alosa) was artificially landlocked during the last century. To conduct this analysis, we assembled and annotated a draft of the A. alosa genome and generated whole‐genome sequencing for 16 anadromous and freshwater populations of shad. Widespread evidence for parallel genetic changes in freshwater populations within lineages was found. In freshwater A. alosa, which have only been diverging for tens of generations, this shows that parallel adaptive evolution can rapidly occur. However, parallel genetic changes across lineages were comparatively rare. The degree of genetic parallelism was not strongly related to the number of shared polymorphisms between lineages, thus suggesting that other factors such as divergence among ancestral populations or environmental variation may influence genetic parallelism across these lineages. These overall patterns were exemplified by genetic differentiation involving a paralog of ATPase‐α1 that appears to be under selection in just two of the more distantly related lineages studied, A. f. lacustris and A. alosa. Our findings provide insights into the genetic architecture of adaptation and parallel evolution along a continuum of population divergence.

Inconsistent evolution and growth-survival tradeoffs in Gambusia affinis

Growth-survival tradeoffs may be a generalizable mechanism influencing trajectories of prey evolution. Here, we investigate evolutionary contributions to growth and survival in western mosquitofish (Gambusia affinis) from 10 populations from high- and low-predation ancestral environments. We assess (i) the degree to which evolutionary components of growth and survival are consistent or inconsistent across populations within ancestral predation environments, and (ii) whether growth and survival trade off at the population level. We measure growth and survival on groups of common-reared mosquitofish in pond mesocosms. We find that evolution of growth is consistent, with fish from low-predation ancestral environments showing higher growth, while the evolution of survival is inconsistent, with significant population-level divergence unrelated to ancestral predation environment. Such inconsistency prevents a growth-survival tradeoff across populations. Thus, the generalizability of contemporary evolution probably depends on local context of evolutionary tradeoffs, and a continued focus on singular selective agents (e.g. predators) without such local context will impede insights into generalizable evolutionary patterns.

Detecting (non)parallel evolution in multidimensional spaces: angles, correlations and eigenanalysis

Parallelism between evolutionary trajectories in a trait space is often seen as evidence for repeatability of phenotypic evolution, and angles between trajectories play a pivotal role in the analysis of parallelism. However, properties of angles in multidimensional spaces have not been widely appreciated by biologists. To remedy this situation, this study provides a brief overview on geometric and statistical aspects of angles in multidimensional spaces. Under the null hypothesis that trajectory vectors have no preferred directions (i.e. uniform distribution on hypersphere), the angle between two independent vectors is concentrated around the right angle, with a more pronounced peak in a higher-dimensional space. This probability distribution is closely related to t- and beta distributions, which can be used for testing the null hypothesis concerning a pair of trajectories. A recently proposed method with eigenanalysis of a vector correlation matrix can be connected to the test of no correlation or concentration of multiple vectors, for which simple test procedures are available in the statistical literature. Concentration of vectors can also be examined by tools of directional statistics such as the Rayleigh test. These frameworks provide biologists with baselines to make statistically justified inferences for (non)parallel evolution.

A phylogeny of Antirrhinum reveals parallel evolution of alpine morphology

Parallel evolution of similar morphologies in closely related lineages provides insight into the repeatability and predictability of evolution. In the genus Antirrhinum (snapdragons), as in other plants, a suite of morphological characters are associated with adaptation to alpine environments. We tested for parallel trait evolution in Antirrhinum by investigating phylogenetic relationships using restriction-site associated DNA (RAD) sequencing. We then associated phenotypic information to our phylogeny to reconstruct the patterns of morphological evolution and related this to evidence for hybridisation between emergent lineages. Phylogenetic analyses showed that the alpine character syndrome is present in multiple groups, suggesting that Antirrhinum has repeatedly colonised alpine habitats. Dispersal to novel environments happened in the presence of intraspecific and interspecific gene flow. We found support for a model of parallel evolution in Antirrhinum. Hybridisation in natural populations, and a complex genetic architecture underlying the alpine morphology syndrome, support an important role of natural selection in maintaining species divergence in the face of gene flow.

Population genetics and independently replicated evolution of predator-associated burst speed ecophenotypy in mosquitofish

Many species show replicated ecophenotypy due to recurring patterns of natural selection. Based on the presence or absence of pursuit predators, at least 17 species of fish repeatedly differentiated in body shape in a manner that increases burst swimming speed and the likelihood of predator escape. The predator-associated burst speed (PABS) ecophenotype is characterized by a small head and trunk and enlarged caudal region. Mechanisms promoting replicated phenotype-environment association include selection (without evolution), a single instance of adaptive evolution followed by biased habitat occupation, repeated instances of local adaptation, or adaptive phenotypic plasticity. Common garden rearing of mosquitofish, Gambusia affinis, demonstrated a likely heritable basis for PABS phenotypy, but it is unknown whether populations are otherwise genetically distinct or whether replicated ecophenotypy represents a single or replicated instances of adaptation. To genetically characterize the populations and test hypotheses of single or multiple adaptations, we characterized variation in 12 polymorphic DNA microsatellites in the previously studied G. affinis populations. Populations were genetically distinct by multilocus analysis, exhibited high allelic diversity, and were heterozygote deficient, which effects were attributed to G. affinis's shoaling nature and habitat patchiness. Genetic and phenotypic distances among populations were correlated for non-PABS but not PABS morphology. Multilocus analysis demonstrated ecophenotype polyphyly and scattered multivariate genetic structure which support only the replicated-adaptation model. As all of the diverse tests performed demonstrated lack of congruence between patterns of molecular genetic and PABS differentiation, it is likely that divergent natural selection drove multiple instances of adaptive evolution.

Rapid, parallel evolution of field mustard (Brassica rapa) under experimental drought

Climate change is driving evolutionary and plastic responses in populations, but predicting these responses remains challenging. Studies that combine experimental evolution with ancestor-descendant comparisons allow assessment of the causes, parallelism, and adaptive nature of evolutionary responses, although such studies remain rare, particularly in a climate change context. Here, we created experimental populations of Brassica rapa derived from the same natural population and exposed these replicated populations to experimental drought or watered conditions for four generations. We then grew ancestors and descendants concurrently, following the resurrection approach. Experimental populations under drought showed rapid evolution of earlier flowering time and increased specific leaf area, consistent with a drought escape strategy and observations in natural populations. Evolutionary shifts followed the direction of selection and increased fitness under drought, indicative of adaptive evolution. Evolution to drought also occurred largely in parallel among replicate populations. Further, traits showed phenotypic plasticity to drought, but the direction and effect size of plasticity varied. Our results demonstrate parallel evolution to experimental drought, suggesting that evolution to strong, consistent selection may be predictable. Broadly, our study demonstrates the utility of combining experimental evolution with the resurrection approach to investigate responses to climate change.

Parallel shape divergence between ecotypes of the limpet Nacella concinna along the Antarctic Peninsula: a new model species for parallel evolution?

Parallel phenotypic divergence is the independent differentiation between phenotypes of the same lineage or species occupying ecologically similar environments in different populations. We tested in the Antarctic limpet Nacella concinna the extent of parallel morphological divergence in littoral and sublittoral ecotypes throughout its distribution range. These ecotypes differ in morphological, behavioural and physiological characteristics. We studied the lateral and dorsal outlines of shells and the genetic variation of the mitochondrial gene Cytochrome Oxidase subunit I from both ecotypes in 17 sample sites along more than 2,000 km. The genetic data indicate that both ecotypes belong to a single evolutionary lineage. The magnitude and direction of phenotypic variation differ between ecotypes across sample sites; completely parallel ecotype-pairs (i.e., they diverge in the same magnitude and in the same direction) were detected in 84.85% of lateral and 65.15% in dorsal view comparisons. Besides, specific traits (relative shell height, position of shell apex, and elliptical/pear-shape outline variation) showed high parallelism. We observed weak morphological covariation between the two shape shell views, indicating that distinct evolutionary forces and environmental pressures could be acting on this limpet shell shape. Our results demonstrate there is a strong parallel morphological divergence pattern in N. concinna along its distribution, making this Antarctic species a suitable model for the study of different evolutionary forces shaping the shell evolution of this limpet.

Speciation dynamics and extent of parallel evolution along a lake-stream environmental contrast in African cichlid fishes

Understanding the dynamics of speciation is a central topic in evolutionary biology. Here, we investigated how morphological and genomic differentiation accumulated along the speciation continuum in the African cichlid fish Astatotilapia burtoni. While morphological differentiation was continuously distributed across different lake-stream population pairs, we found that there were two categories with respect to genomic differentiation, suggesting a “gray zone” of speciation at ~0.1% net nucleotide divergence. Genomic differentiation was increased in the presence of divergent selection and drift compared to drift alone. The quantification of phenotypic and genetic parallelism in four cichlid species occurring along a lake-stream environmental contrast revealed parallel and antiparallel components in rapid adaptive divergence, and morphological convergence in species replicates inhabiting the same environments. Furthermore, we show that the extent of parallelism was higher when ancestral populations were more similar. Our study highlights the complementary roles of divergent selection and drift on speciation and parallel evolution.

Phenotypic and genotypic parallel evolution in parapatric ecotypes of Senecio

The independent and repeated adaptation of populations to similar environments often results in the evolution of similar forms. This phenomenon creates a strong correlation between phenotype and environment and is referred to as parallel evolution. However, we are still largely unaware of the dynamics of parallel evolution, as well as the interplay between phenotype and genotype within natural systems. Here, we examined phenotypic and genotypic parallel evolution in multiple parapatric Dune‐Headland coastal ecotypes of an Australian wildflower, Senecio lautus. We observed a clear trait‐environment association in the system, with all replicate populations having evolved along the same phenotypic evolutionary trajectory. Similar phenotypes have arisen via mutational changes occurring in different genes, although many share the same biological functions. Our results shed light on how replicated adaptation manifests at the phenotypic and genotypic levels within populations, and highlight S. lautus as one of the most striking cases of phenotypic parallel evolution in nature.

Parallel evolution of morphological traits and body shape in littoral and pelagic brook charr, Salvelinus fontinalis, along a gradient of interspecific competition

The parallel evolution of similar ecotypes in response to comparable environmental conditions is believed to reveal the importance of divergent selection in phenotypic diversifying processes. Systems characterized by the presence of multiple replicate populations expressing resource polymorphism thus provide an ideal opportunity to address the occurrence and factors affecting the parallel evolution of ecotypes. Previous studies have shown that brook charr (Salvelinus fontinalis) exhibit resource polymorphism in some Canadian Shield lakes, where a littoral ecotype feeds mainly on zoobenthos and a pelagic ecotype feeds mostly on zooplankton. Using morphological traits and geometric morphometric analyses on 18 native brook charr populations, we explicitly tested (i) whether brook charr ecotypes show parallel evolution across populations (i.e. the same morphological traits discriminate ecotypes among lakes) and (ii) whether interspecific competition decreases the amplitude of morphological differentiation between ecotypes, if any, because brook charr experience some level of competitive exclusion from the littoral habitat in the presence of creek chub or white sucker. We observed a low level of parallel evolution, where the littoral ecotype was overall stouter with longer fins and smaller eyes than the pelagic ecotype. Interspecific competition had no clear impacts on the amplitude of morphological differentiation. We also observed that inter-lake morphological differences are greater than between ecotypes within lakes, suggesting an important effect of local environmental factors on population morphology. Early-stage of diversification as well as phenotypic plasticity and morphological integration could explain why resource polymorphism is still subtle in brook charr populations.

Natural and anthropogenic sources of habitat variation influence exploration behaviour, stress response, and brain morphology in a coastal fish

Evolutionary ecology aims to better understand how ecologically important traits respond to environmental heterogeneity. Environments vary both naturally and as a result of human activities, and investigations that simultaneously consider how natural and human-induced environmental variation affect diverse trait types grow increasingly important as human activities drive species endangerment. Here, we examined how habitat fragmentation and structural habitat complexity affect disparate trait types in Bahamas mosquitofish Gambusia hubbsi inhabiting tidal creeks. We tested a priori predictions for how these factors might influence exploratory behaviour, stress reactivity and brain anatomy. We examined approximately 350 adult Bahamas mosquitofish from seven tidal-creek populations across Andros Island, The Bahamas that varied in both human-caused fragmentation (three fragmented and four unfragmented) and natural habitat complexity (e.g. fivefold variation in rock habitat). Populations that had experienced severe human-induced fragmentation, and thus restriction of tidal exchange from the ocean, exhibited greater exploration of a novel environment, stronger physiological stress responses to a mildly stressful event and smaller telencephala (relative to body size). These changes matched adaptive predictions based mostly on (a) reduced chronic predation risk and (b) decreased demands for navigating tidally dynamic habitats. Populations from sites with greater structural habitat complexity showed a higher propensity for exploration and a relatively larger optic tectum and cerebellum. These patterns matched adaptive predictions related to increased demands for navigating complex environments. Our findings demonstrate environmental variation, including recent anthropogenic impacts (<50 years), can significantly affect complex, ecologically important traits. Yet trait-specific patterns may not be easily predicted, as we found strong support for only six of 12 predictions. Our results further highlight the utility of simultaneously quantifying multiple environmental factors-for example had we failed to account for habitat complexity, we would not have detected the effects of fragmentation on exploratory behaviours. These responses, and their ecological consequences, may be complex: rapid and adaptive phenotypic responses to anthropogenic impacts can facilitate persistence in human-altered environments, but may come at a cost of population vulnerability if ecological restoration was to occur without consideration of the altered traits.

Capturing the rapidly evolving study of adaptation

Research on the genomics of adaptation is rapidly changing. In the last few decades, progress in this area has been driven by methodological advances, not only in the way increasingly large amounts of molecular data are generated (e.g. with high-throughput sequencing), but also in the way these data are analysed. This includes a growing appreciation and quantitative treatment of covariation among units within the same data type (e.g. genes) or across data types (e.g. genes and phenotypes). The development and adoption of more and more integrative tools have resulted in richer and more interesting empirical work. This special issue - comprising methodological, empirical, and review papers - aims to capture a 'snapshot' of this rapidly evolving field. We discuss in particular three important themes in the study of adaptation: the genetic architecture of adaptive variation, protein-coding and regulatory changes, and parallel evolution. We highlight how more traditional key themes in the study of genetic architecture (e.g. the number of loci underlying adaptive traits and the distribution of their effects) are now being complemented by other factors (e.g. how patterns of linkage and number of loci interact to affect the ability to adapt). Similarly, apart from addressing the relative importance of protein-coding and regulatory changes, we now have the tools to look in-depth at specific types of regulatory variation to gain a clearer picture of regulatory networks. Finally, parallel evolution has always been central to the study of adaptation, but now we are often able to address the question of whether - and to what extent - parallelism at the organismal or phenotypic level is matched by parallelism at the genetic level. Perhaps most importantly, we can now determine what mechanisms are driving parallelism (or lack thereof) across levels of biological organization. All these recent methodological developments open up new directions for future studies of adaptive changes across traits, levels of biological organization, demographic contexts and time scales.

The complex ecology of genitalia: Gonopodium length and allometry in the Trinidadian guppy

Male genitalia present an extraordinary pattern of rapid divergence in animals with internal fertilization, which is usually attributed to sexual selection. However, the effect of ecological factors on genitalia divergence could also be important, especially so in animals with nonretractable genitalia because of their stronger interaction with the surrounding environment in comparison with animals with retractable genitalia. Here, we examine the potential of a pervasive ecological factor (predation) to influence the length and allometry of the male genitalia in guppies. We sampled guppies from pairs of low-predation (LP) and high-predation (HP) populations in seven rivers in Trinidad, and measured their body and gonopodium length. A key finding was that HP adult males do not have consistently longer gonopodia than do LP adult males, as had been described in previous work. However, we did find such divergence for juvenile males: HP juveniles have longer gonopodia than do LP juveniles. We therefore suggest that an evolutionary trend toward the development of longer gonopodia in HP males (as seen in the juveniles) is erased after maturity owing to the higher mortality of mature males with longer gonopodia. Beyond these generalities, gonopodium length and gonopodium allometry were remarkably variable among populations even within a predation regime, thus indicating strong context dependence to their development/evolution. Our findings highlight the complex dynamics of genitalia evolution in Trinidadian guppies.

Adaptive Divergence under Gene Flow along an Environmental Gradient in Two Coexisting Stickleback Species

There is a general and solid theoretical framework to explain how the interplay between natural selection and gene flow affects local adaptation. Yet, to what extent coexisting closely related species evolve collectively or show distinctive evolutionary responses remains a fundamental question. To address this, we studied the population genetic structure and morphological differentiation of sympatric three-spined and nine-spined stickleback. We conducted genotyping-by-sequencing and morphological trait characterisation using 24 individuals of each species from four lowland brackish water (LBW), four lowland freshwater (LFW) and three upland freshwater (UFW) sites in Belgium and the Netherlands. This combination of sites allowed us to contrast populations from isolated but environmentally similar locations (LFW vs. UFW), isolated but environmentally heterogeneous locations (LBW vs. UFW), and well-connected but environmentally heterogenous locations (LBW vs. LFW). Overall, both species showed comparable levels of genetic diversity and neutral genetic differentiation. However, for all three spatial scales, signatures of morphological and genomic adaptive divergence were substantially stronger among populations of the three-spined stickleback than among populations of the nine-spined stickleback. Furthermore, most outlier SNPs in the two species were associated with local freshwater sites. The few outlier SNPs that were associated with the split between brackish water and freshwater populations were located on one linkage group in three-spined stickleback and two linkage groups in nine-spined stickleback. We conclude that while both species show congruent evolutionary and genomic patterns of divergent selection, both species differ in the magnitude of their response to selection regardless of the geographical and environmental context.

Intercontinental genomic parallelism in multiple three-spined stickleback adaptive radiations

Parallelism, the evolution of similar traits in populations diversifying in similar conditions, provides strong evidence of adaptation by natural selection. Many studies of parallelism focus on comparisons of different ecotypes or contrasting environments, defined a priori, which could upwardly bias the apparent prevalence of parallelism. Here, we estimated genomic parallelism associated with components of environmental and phenotypic variation at an intercontinental scale across four freshwater adaptive radiations (Alaska, British Columbia, Iceland and Scotland) of the three-spined stickleback (Gasterosteus aculeatus). We combined large-scale biological sampling and phenotyping with restriction site associated DNA sequencing (RAD-Seq) data from 73 freshwater lake populations and four marine ones (1,380 fish) to associate genome-wide allele frequencies with continuous distributions of environmental and phenotypic variation. Our three main findings demonstrate that (1) quantitative variation in phenotypes and environments can predict genomic parallelism; (2) genomic parallelism at the early stages of adaptive radiations, even at large geographic scales, is founded on standing variation; and (3) similar environments are a better predictor of genome-wide parallelism than similar phenotypes. Overall, this study validates the importance and predictive power of major phenotypic and environmental factors likely to influence the emergence of common patterns of genomic divergence, providing a clearer picture than analyses of dichotomous phenotypes and environments.

Putative chromosomal rearrangements are associated primarily with ecotype divergence rather than geographic separation in an intertidal, poorly dispersing snail

Littorina saxatilis is becoming a model system for understanding the genomic basis of ecological speciation. The parallel formation of crab-adapted ecotypes that exhibit partial reproductive isolation from wave-adapted ecotypes has enabled genomic investigation of conspicuous shell traits. Recent genomic studies suggest that chromosomal rearrangements may enable ecotype divergence by reducing gene flow. However, the genomic architecture of traits that are divergent between ecotypes remains poorly understood. Here, we use 11,504 single nucleotide polymorphism (SNP) markers called using the recently released L. saxatilis genome to genotype 462 crab ecotype, wave ecotype and phenotypically intermediate Spanish L. saxatilis individuals with scored phenotypes. We used redundancy analysis to study the genetic architecture of loci associated with shell shape, shape corrected for size, shell size and shell ornamentation, and to compare levels of co-association among different traits. We discovered 341 SNPs associated with shell traits. Loci associated with trait divergence between ecotypes were often located inside putative chromosomal rearrangements recently characterized in Swedish L. saxatilis. In contrast, we found that shell shape corrected for size varied primarily by geographic site rather than by ecotype and showed little association with these putative rearrangements. We conclude that genomic regions of elevated divergence inside putative rearrangements were associated with divergence of L. saxatilis ecotypes along steep environmental axes-consistent with models of adaptation with gene flow-but were not associated with divergence among the three geographical sites. Our findings support predictions from models indicating the importance of genomic regions of reduced recombination allowing co-association of loci during ecological speciation with ongoing gene flow.

The role of introgression and ecotypic parallelism in delineating intraspecific conservation units

Parallel evolution can occur through selection on novel mutations, standing genetic variation or adaptive introgression. Uncovering parallelism and introgressed populations can complicate management of threatened species as parallelism may have influenced conservation unit designations and admixed populations are not generally considered under legislations. We examined high coverage whole-genome sequences of 30 caribou (Rangifer tarandus) from across North America and Greenland, representing divergent intraspecific lineages, to investigate parallelism and levels of introgression contributing to the formation of ecotypes. Caribou are split into four subspecies and 11 extant conservation units, known as designatable units (DUs), in Canada. Using genomes from all four subspecies and six DUs, we undertake demographic reconstruction and confirm two previously inferred instances of parallel evolution in the woodland subspecies and uncover an additional instance of parallelism of the eastern migratory ecotype. Detailed investigations reveal introgression in the woodland subspecies, with introgressed regions found spread throughout the genomes encompassing both neutral and functional sites. Our investigations using whole genomes highlight the difficulties in unequivocally demonstrating parallelism through adaptive introgression in nonmodel species with complex demographic histories, with standing variation and introgression both potentially involved. Additionally, the impact of parallelism and introgression on conservation policy for management units needs to be considered in general, and the caribou designations will need amending in light of our results. Uncovering and decoupling parallelism and differential patterns of introgression will become prevalent with the availability of comprehensive genomic data from nonmodel species, and we highlight the need to incorporate this into conservation unit designations.

Parallelism in eco-morphology and gene expression despite variable evolutionary and genomic backgrounds in a Holarctic fish

Understanding the extent to which ecological divergence is repeatable is essential for predicting responses of biodiversity to environmental change. Here we test the predictability of evolution, from genotype to phenotype, by studying parallel evolution in a salmonid fish, Arctic charr (Salvelinus alpinus), across eleven replicate sympatric ecotype pairs (benthivorous-planktivorous and planktivorous-piscivorous) and two evolutionary lineages. We found considerable variability in eco-morphological divergence, with several traits related to foraging (eye diameter, pectoral fin length) being highly parallel even across lineages. This suggests repeated and predictable adaptation to environment. Consistent with ancestral genetic variation, hundreds of loci were associated with ecotype divergence within lineages of which eight were shared across lineages. This shared genetic variation was maintained despite variation in evolutionary histories, ranging from postglacial divergence in sympatry (ca. 10-15kya) to pre-glacial divergence (ca. 20-40kya) with postglacial secondary contact. Transcriptome-wide gene expression (44,102 genes) was highly parallel across replicates, involved biological processes characteristic of ecotype morphology and physiology, and revealed parallelism at the level of regulatory networks. This expression divergence was not only plastic but in part genetically controlled by parallel cis-eQTL. Lastly, we found that the magnitude of phenotypic divergence was largely correlated with the genetic differentiation and gene expression divergence. In contrast, the direction of phenotypic change was mostly determined by the interplay of adaptive genetic variation, gene expression, and ecosystem size. Ecosystem size further explained variation in putatively adaptive, ecotype-associated genomic patterns within and across lineages, highlighting the role of environmental variation and stochasticity in parallel evolution. Together, our findings demonstrate the parallel evolution of eco-morphology and gene expression within and across evolutionary lineages, which is controlled by the interplay of environmental stochasticity and evolutionary contingencies, largely overcoming variable evolutionary histories and genomic backgrounds.

Expanding the mutualistic niche: parallel symbiont turnover along climatic gradients

Keystone mutualisms, such as corals, lichens or mycorrhizae, sustain fundamental ecosystem functions. Range dynamics of these symbioses are, however, inherently difficult to predict because host species may switch between different symbiont partners in different environments, thereby altering the range of the mutualism as a functional unit. Biogeographic models of mutualisms thus have to consider both the ecological amplitudes of various symbiont partners and the abiotic conditions that trigger symbiont replacement. To address this challenge, we here investigate 'symbiont turnover zones'--defined as demarcated regions where symbiont replacement is most likely to occur, as indicated by overlapping abundances of symbiont ecotypes. Mapping the distribution of algal symbionts from two species of lichen-forming fungi along four independent altitudinal gradients, we detected an abrupt and consistent β-diversity turnover suggesting parallel niche partitioning. Modelling contrasting environmental response functions obtained from latitudinal distributions of algal ecotypes consistently predicted a confined altitudinal turnover zone. In all gradients this symbiont turnover zone is characterized by approximately 12°C average annual temperature and approximately 5°C mean temperature of the coldest quarter, marking the transition from Mediterranean to cool temperate bioregions. Integrating the conditions of symbiont turnover into biogeographic models of mutualisms is an important step towards a comprehensive understanding of biodiversity dynamics under ongoing environmental change.

Effects of predation risk on egg steroid profiles across multiple populations of threespine stickleback

Predation often has consistent effects on prey behavior and morphology, but whether the physiological mechanisms underlying these effects show similarly consistent patterns across different populations remains an open question. In vertebrates, predation risk activates the hypothalamic-pituitary-adrenal (HPA) axis, and there is growing evidence that activation of the maternal HPA axis can have intergenerational consequences via, for example, maternally-derived steroids in eggs. Here, we investigated how predation risk affects a suite of maternally-derived steroids in threespine stickleback eggs across nine Alaskan lakes that vary in whether predatory trout are absent, native, or have been stocked within the last 25 years. Using liquid chromatography coupled with mass spectroscopy (LC-MS/MS), we detected 20 steroids within unfertilized eggs. Factor analysis suggests that steroids covary within and across steroid classes (i.e. glucocorticoids, progestogens, sex steroids), emphasizing the modularity and interconnectedness of the endocrine response. Surprisingly, egg steroid profiles were not significantly associated with predator regime, although they were more variable when predators were absent compared to when predators were present, with either native or stocked trout. Despite being the most abundant steroid, cortisol was not consistently associated with predation regime. Thus, while predators can affect steroids in adults, including mothers, the link between maternal stress and embryonic development is more complex than a simple one-to-one relationship between the population-level predation risk experienced by mothers and the steroids mothers transfer to their eggs.

Genomic architecture of parallel ecological divergence: Beyond a single environmental contrast

The study of parallel ecological divergence provides important clues to the operation of natural selection. Parallel divergence often occurs in heterogeneous environments with different kinds of environmental gradients in different locations, but the genomic basis underlying this process is unknown. We investigated the genomics of rapid parallel adaptation in the marine snail Littorina saxatilis in response to two independent environmental axes (crab-predation versus wave-action and low-shore versus high-shore). Using pooled whole-genome resequencing, we show that sharing of genomic regions of high differentiation between environments is generally low but increases at smaller spatial scales. We identify different shared genomic regions of divergence for each environmental axis and show that most of these regions overlap with candidate chromosomal inversions. Several inversion regions are divergent and polymorphic across many localities. We argue that chromosomal inversions could store shared variation that fuels rapid parallel adaptation to heterogeneous environments, possibly as balanced polymorphism shared by adaptive gene flow.

Evidence for host effect on the intestinal microbiota of whitefish (Coregonus sp.) species pairs and their hybrids

Investigating relationships between microbiota and their host is essential toward a full understanding of how animal adapt to their environment. Lake Whitefish offers a powerful system to investigate processes of adaptive divergence where the dwarf, limnetic species evolved repeatedly from the normal, benthic species. We compared the transient intestinal microbiota between both species from the wild and in controlled conditions, including their reciprocal hybrids. We sequenced the 16s rRNA gene V3-V4 regions to (a) test for parallelism in the transient intestinal microbiota among sympatric pairs, (b) test for transient intestinal microbiota differences among dwarf, normal, and hybrids reared under identical conditions, and (c) compare intestinal microbiota between wild and captive whitefish. A significant host effect on microbiota taxonomic composition was observed when all lakes were analyzed together and in three of the five species pairs. In captive whitefish, host effect was also significant. Microbiota of both reciprocal hybrids fell outside of that observed in the parental forms. Six genera formed a bacterial core which was present in captive and wild whitefish, suggesting a horizontal microbiota transmission. Altogether, our results complex interactions among the host, the microbiota, and the environment, and we propose that these interactions define three distinct evolutionary paths of the intestinal microbiota.

100-year time series reveal little morphological change following impoundment and predator invasion in two Neotropical characids

Human activities are dramatically altering ecosystems worldwide, often resulting in shifts in selection regimes. In response, natural populations sometimes undergo rapid phenotypic changes, which, if adaptive, can increase their probability of persistence. However, in many instances, populations fail to undergo any phenotypic change, which might indicate a variety of possibilities, including maladaptation. In freshwater ecosystems, the impoundment of rivers and the introduction of exotic species are among the leading threats to native fishes. We examined how the construction of the Panama Canal, which formed Lake Gatun, and the subsequent invasion of the predatory Cichla monoculus influenced the morphology of two native fishes: Astyanax ruberrimus and Roeboides spp. Using a 100-year time series, we studied variation in overall body shape over time (before vs. after impoundment and invasion) as well as across space (between an invaded and an uninvaded reservoir). In addition, we examined variation in linear morphological traits associated with swim performance and predator detection/avoidance. Notwithstanding a few significant changes in particular traits in particular comparisons, we found only limited evidence for morphological change associated with these two stressors. Most observed changes were subtle, and tended to be site- and species-specific. The lack of a strong morphological response to these stressors, coupled with dramatic population declines in both species, suggests they may be maladapted to the anthropogenically perturbed environment of Lake Gatun, but direct measures of fitness would be needed to test this. In general, our results suggest that morphological responses to anthropogenic disturbances can be very limited and, when they do occur, are often complex and context-dependent.

Independent lineages in a common environment: the roles of determinism and contingency in shaping the migration timing of even- versus odd-year pink salmon over broad spatial and temporal scales

Studies of parallel evolution are seldom able to disentangle the influence of cryptic environmental variation from that of evolutionary history; whereas the unique life history of pink salmon (Oncorhynchus gorbuscha) presents an opportunity to do so. All pink salmon mature at age two and die after breeding. Hence, pink salmon bred in even years are completely reproductively isolated from those bred in odd years, even if the two lineages bred in same location. We used time series (mean = 7 years, maximum = 74 years) of paired even- and odd-year populations from 36 rivers spanning over 2000 km to explore parallelism in migration timing, a trait with a strong genetic basis. Migration timing was highly parallel, being determined almost entirely by local environmental differences among rivers. Interestingly, interannual changes in migration timing different somewhat between lineages. Overall, our findings indicate very strong determinism, with only a minor contribution of contingency.

Molecular and morphological convergence to sulfide-tolerant fishes in a new species of Jenynsia (Cyprinodontiformes: Anablepidae), the first extremophile member of the family

Freshwater sulfide springs have extreme environmental conditions that only few vertebrate species can tolerate. These species often develop a series of morphological and molecular adaptations to cope with the challenges of life under the toxic and hypoxic conditions of sulfide springs. In this paper, we described a new fish species of the genus Jenynsia, Anablepidae, from a sulfide spring in Northwestern Argentina, the first in the family known from such extreme environment. Jenynsia sulfurica n. sp. is diagnosable by the lack of scales on the pre-pelvic area or the presence of a single row of scales, continuous or not, from the isthmus to the bases of the pelvic fins. Additionally, it presents a series of morphological and molecular characteristics that appear convergent with those seen in other fish species (e.g., Poeciliids) inhabiting sulfide springs. Most notably, J. sulfurica has an enlarged head and postorbital area compared to other fish of the genus and a prognathous lower jaw with a hypertrophied lip, thought to facilitate respiration at the air-water interface. Analyses of cox1 sequence showed that J. sulfurica has two unique mutations resulting in amino acid substitutions convergent to those seen in Poeciliids from sulfide springs and known to provide a physiological mechanism related to living in sulfide environments. A phylogenetic analysis, including molecular and morphological characters, placed J. sulfurica as sister taxa to J. alternimaculata, a species found in nearby, non-sulfide habitats directly connected to the sulfide springs. Thus, it can be inferred that the selection imposed by the presence of H2S has resulted in the divergence between these two species and has potentially served as a barrier to gene flow.

A roadmap for urban evolutionary ecology

Urban ecosystems are rapidly expanding throughout the world, but how urban growth affects the evolutionary ecology of species living in urban areas remains largely unknown. Urban ecology has advanced our understanding of how the development of cities and towns change environmental conditions and alter ecological processes and patterns. However, despite decades of research in urban ecology, the extent to which urbanization influences evolutionary and eco-evolutionary change has received little attention. The nascent field of urban evolutionary ecology seeks to understand how urbanization affects the evolution of populations, and how those evolutionary changes in turn influence the ecological dynamics of populations, communities, and ecosystems. Following a brief history of this emerging field, this Perspective article provides a research agenda and roadmap for future research aimed at advancing our understanding of the interplay between ecology and evolution of urban-dwelling organisms. We identify six key questions that, if addressed, would significantly increase our understanding of how urbanization influences evolutionary processes. These questions consider how urbanization affects nonadaptive evolution, natural selection, and convergent evolution, in addition to the role of urban environmental heterogeneity on species evolution, and the roles of phenotypic plasticity versus adaptation on species' abundance in cities. Our final question examines the impact of urbanization on evolutionary diversification. For each of these six questions, we suggest avenues for future research that will help advance the field of urban evolutionary ecology. Lastly, we highlight the importance of integrating urban evolutionary ecology into urban planning, conservation practice, pest management, and public engagement.

Parallel genetic evolution and speciation from standing variation

Adaptation often proceeds from standing variation, and natural selection acting on pairs of populations is a quantitative continuum ranging from parallel to divergent. Yet, it is unclear how the extent of parallel genetic evolution during adaptation from standing variation is affected by the difference in the direction of selection between populations. Nor is it clear whether the availability of standing variation for adaptation affects progress toward speciation in a manner that depends on the difference in the direction of selection. We conducted a theoretical study investigating these questions and have two primary findings. First, the extent of parallel genetic evolution between two populations rapidly declines as selection changes from fully parallel toward divergent, and this decline is steeper in organisms with more traits (i.e., greater dimensionality). This rapid decline happens because small differences in the direction of selection greatly reduce the fraction of alleles that are beneficial in both populations. For example, populations adapting to optima separated by an angle of 33° might have only 50% of potentially beneficial alleles in common. Second, relative to when adaptation is from only new mutation, adaptation from standing variation improves hybrid fitness under parallel selection and reduces hybrid fitness under divergent selection. Under parallel selection, genetic parallelism from standing variation reduces the phenotypic segregation variance in hybrids, thereby increasing mean fitness in the parental environment. Under divergent selection, larger pleiotropic effects of alleles fixed from standing variation cause maladaptive transgressive phenotypes when combined in hybrids. Adaptation from standing genetic variation therefore slows progress toward speciation under parallel selection and facilitates progress toward speciation under divergent selection.

Individual variation in reproductive behaviour is linked to temporal heterogeneity in predation risk

Variation in predation risk is a major driver of ecological and evolutionary change, and, in turn, of geographical variation in behaviour. While predation risk is rarely constant in natural populations, the extent to which variation in predation risk shapes individual behaviour in wild populations remains unclear. Here, we investigated individual differences in reproductive behaviour in 16 Trinidadian guppy populations and related it to the observed variation in predator biomass each population experienced. Our results show that high heterogeneity in predator biomass is linked to individual behavioural diversification. Increased within-population heterogeneity in predator biomass is also associated with behavioural polymorphism. Some individuals adjust the frequency of consensual mating behaviour in response to differences in sex ratio context, while others display constantly at elevated frequencies. This pattern is analogous to a 'live fast, die young' pace-of-life syndrome. Notably, both high and low mean differences in predator biomass led to a homogenization of individual frequency of consensual mating displays. Overall, our results demonstrate that individual behavioural variation is associated with heterogeneity in predator biomass, but not necessarily with changes in mean values of predator biomass. We suggest that heterogeneity in predator biomass is an informative predictor of adaptive responses to changes in biotic conditions.

Do replicates of independent guppy lineages evolve similarly in a predator-free laboratory environment?

The Trinidadian guppy is emblematic of parallel and convergent evolution, with repeated demonstrations that predation regime is a driver of adaptive trait evolution. A classic and foundational experiment in this system was conducted by John Endler 40 years ago, where male guppies placed into low-predation environments in the laboratory evolved increased color in a few generations. However, Endler's experiment did not employ the now typical design for a parallel/convergent evolution study, which would employ replicates of different ancestral lineages. We therefore implemented an experiment that seeded replicate mesocosms with small founding populations of guppies originating from high-predation populations of two very different lineages. The different mesocosms were maintained identically, and male guppy color was quantified every four months. After one year, we tested whether male color had increased, whether replicates within a lineage had parallel phenotypic trajectories, and whether the different lineages converged on a common phenotype. Results showed that male guppy color generally increased through time, primarily due to changes in melanic color, whereas the other colors showed inconsistent and highly variable trajectories. Most of the nonparallelism in phenotypic trajectories was among mesocosms containing different lineages. In addition to this mixture of parallelism and nonparallelism, convergence was not evident in that the variance in color among the mesocosms actually increased through time. We suggest that our results reflect the potential importance of high variation in female preference and stochastic processes such as drift and founder effects, both of which could be important in nature.