Extent of QTL Reuse During Repeated Phenotypic Divergence of Sympatric Threespine Stickleback

Pleiotropy increases parallel selection signatures during adaptation from standing genetic variation

The phenomenon of parallel evolution, whereby similar genomic and phenotypic changes occur across replicated pairs of populations or species, is widely studied. Nevertheless, the determining factors of parallel evolution remain poorly understood. Theoretical studies have proposed that pleiotropy, the influence of a single gene on multiple traits, is an important factor. In order to gain a deeper insight into the role of pleiotropy for parallel evolution from standing genetic variation, we characterized the interplay between parallelism, polymorphism, and pleiotropy. The present study examined the parallel gene expression evolution in 10 replicated populations of Drosophila simulans, which were adapted from standing variation to the same new temperature regime. The data demonstrate that the parallel evolution of gene expression from standing genetic variation is positively correlated with the strength of pleiotropic effects. The ancestral variation in gene expression is, however, negatively correlated with parallelism. Given that pleiotropy is also negatively correlated with gene expression variation, we conducted a causal analysis to distinguish cause and correlation and evaluate the role of pleiotropy. The causal analysis indicated that both direct (causative) and indirect (correlational) effects of pleiotropy contribute to parallel evolution. The indirect effect is mediated by historic selective constraint in response to pleiotropy. This results in parallel selection responses due to the reduced standing variation of pleiotropic genes. The direct effect of pleiotropy is likely to reflect a genetic correlation among adaptive traits, which in turn gives rise to synergistic effects and higher parallelism.

Inversions contribute disproportionately to parallel genomic divergence in dune sunflowers

The probability of parallel genetic evolution is a function of the strength of selection and constraints imposed by genetic architecture. Inversions capture locally adapted alleles and suppress recombination between them, which limits the range of adaptive responses. In addition, the combined phenotypic effect of alleles within inversions is likely to be greater than that of individual alleles; this should further increase the contributions of inversions to parallel evolution. We tested the hypothesis that inversions contribute disproportionately to parallel genetic evolution in independent dune ecotypes of Helianthus petiolaris. We analysed habitat data and identified variables underlying parallel habitat shifts. Genotype-environment association analyses of these variables indicated parallel responses of inversions to shared selective pressures. We also confirmed larger seed size across the dunes and performed quantitative trait locus mapping with multiple crosses. Quantitative trait loci shared between locations fell into inversions more than expected by chance. We used whole-genome sequencing data to identify selective sweeps in the dune ecotypes and found that the majority of shared swept regions were found within inversions. Phylogenetic analyses of shared regions indicated that within inversions, the same allele typically was found in the dune habitat at both sites. These results confirm predictions that inversions drive parallel divergence in the dune ecotypes.

Non-random mating behaviour between diverging littoral and pelagic three-spined sticklebacks in an invasive population from Upper Lake Constance

Adaptive divergence and increased genetic differentiation among populations can lead to reproductive isolation. In Lake Constance, Germany, a population of invasive three-spined stickleback (Gasterosteus aculeatus) is currently diverging into littoral and pelagic ecotypes, which both nest in the littoral zone. We hypothesized that assortative mating behaviour contributes to reproductive isolation between these ecotypes and performed a behavioural experiment in which females could choose between two nest-guarding males. Behaviour was recorded, and data on traits relevant to mate choice were collected. Both females of the same and different ecotypes were courted with equal vigour. However, there was a significant interaction effect of male and female ecotypes on the level of aggression in females. Littoral females were more aggressive towards pelagic males, and pelagic females were more aggressive towards littoral males. This indicates rejection of males of different ecotypes in spite of the fact that littoral males were larger, more intensely red-coloured and more aggressive than the pelagic males-all mating traits female sticklebacks generally select for. This study documents the emergence of behavioural barriers during early divergence in an invasive and rapidly diversifying stickleback population and discusses their putative role in facilitating reproductive isolation and adaptive radiation within this species.

Heterogeneous genomic architecture of skeletal armour traits in sticklebacks

Whether populations adapt to similar selection pressures using the same underlying genetic variants depends on population history and the distribution of standing genetic variation at the metapopulation level. Studies of sticklebacks provide a case in point: when colonizing and adapting to freshwater habitats, three-spined sticklebacks (Gasterosteus aculeatus) with high gene flow tend to fix the same adaptive alleles in the same major loci, whereas nine-spined sticklebacks (Pungitius pungitius) with limited gene flow tend to utilize a more heterogeneous set of loci. In accordance with this, we report results of quantitative trait locus (QTL) analyses using a backcross design showing that lateral plate number variation in the western European nine-spined sticklebacks mapped to 3 moderate-effect QTL, contrary to the major-effect QTL in three-spined sticklebacks and different from the 4 QTL previously identified in the eastern European nine-spined sticklebacks. Furthermore, several QTL were identified associated with variation in lateral plate size, and 3 moderate-effect QTL with body size. Together, these findings indicate more heterogenous and polygenic genetic underpinnings of skeletal armour variation in nine-spined than three-spined sticklebacks, indicating limited genetic parallelism underlying armour trait evolution in the family Gasterostidae.

Mate choice in the brain: species differ in how male traits 'turn on' gene expression in female brains

Mate choice plays a fundamental role in speciation, yet we know little about the molecular mechanisms that underpin this crucial decision-making process. Stickleback fish differentially adapted to limnetic and benthic habitats are reproductively isolated and females of each species use different male traits to evaluate prospective partners and reject heterospecific males. Here, we integrate behavioural data from a mate choice experiment with gene expression profiles from the brains of females actively deciding whether to mate. We find substantial gene expression variation between limnetic and benthic females, regardless of behavioural context, suggesting general divergence in constitutive gene expression patterns, corresponding to their genetic differentiation. Intriguingly, female gene co-expression modules covary with male display traits but in opposing directions for sympatric populations of the two species, suggesting male displays elicit a dynamic neurogenomic response that reflects known differences in female preferences. Furthermore, we confirm the role of numerous candidate genes previously implicated in female mate choice in other species, suggesting evolutionary tinkering with these conserved molecular processes to generate divergent mate preferences. Taken together, our study adds important new insights to our understanding of the molecular processes underlying female decision-making critical for generating sexual isolation and speciation.

A theory of oligogenic adaptation of a quantitative trait

Rapid phenotypic adaptation is widespread in nature, but the underlying genetic dynamics remain controversial. Whereas population genetics envisages sequential beneficial substitutions, quantitative genetics assumes a collective response through subtle shifts in allele frequencies. This dichotomy of a monogenic and a highly polygenic view of adaptation raises the question of a middle ground, as well as the factors controlling the transition. Here, we consider an additive quantitative trait with equal locus effects under Gaussian stabilizing selection that adapts to a new trait optimum after an environmental change. We present an analytical framework based on Yule branching processes to describe how phenotypic adaptation is achieved by collective changes in allele frequencies at the underlying loci. In particular, we derive an approximation for the joint allele-frequency distribution conditioned on the trait mean as a comprehensive descriptor of the adaptive architecture. Depending on the model parameters, this architecture reproduces the well-known patterns of sequential, monogenic sweeps, or of subtle, polygenic frequency shifts. Between these endpoints, we observe oligogenic architecture types that exhibit characteristic patterns of partial sweeps. We find that a single compound parameter, the population-scaled background mutation rate Θbg, is the most important predictor of the type of adaptation, while selection strength, the number of loci in the genetic basis, and linkage only play a minor role.

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.

Genetic basis of ecologically relevant body shape variation among four genera of cichlid fishes

Divergence in body shape is one of the most widespread and repeated patterns of morphological variation in fishes and is associated with habitat specification and swimming mechanics. Such ecological diversification is the first stage of the explosive adaptive radiation of cichlid fishes in the East African Rift Lakes. We use two hybrid crosses of cichlids (Metriaclima sp. × Aulonocara sp. and Labidochromis sp. × Labeotropheus sp., >975 animals total) to determine the genetic basis of body shape diversification that is similar to benthic-pelagic divergence across fishes. Using a series of both linear and geometric shape measurements, we identified 34 quantitative trait loci (QTL) that underlie various aspects of body shape variation. These QTL are spread throughout the genome, each explaining 3.2-8.6% of phenotypic variation, and are largely modular. Further, QTL are distinct both between these two crosses of Lake Malawi cichlids and compared to previously identified QTL for body shape in fishes such as sticklebacks. We find that body shape is controlled by many genes of small effect. In all, we find that convergent body shape phenotypes commonly observed across fish clades are most likely due to distinct genetic and molecular mechanisms.

How predictable is adaptation from standing genetic variation? Experimental evolution in Drosophila highlights the central role of redundancy and linkage disequilibrium

Experimental evolution is well-suited to test the predictability of evolution without the confounding effects of inaccurate forecasts about future environments. Most of the literature about parallel (and thus predictable) evolution has been carried out in asexual microorganisms, which adapt by de novo mutations. Nevertheless, parallel evolution has also been studied in sexual species at the genomic level. Here, I review the evidence for parallel evolution in Drosophila, the best-studied obligatory outcrossing model for adaptation from standing genetic variation in the laboratory. Similar to asexual microorganisms, evidence for parallel evolution varies between the focal hierarchical levels. Selected phenotypes consistently respond in a very predicable way, but the underlying allele frequency changes are much less predictable. The most important insight is that the predictability of the genomic selection response for polygenic traits depends highly on the founder population and to a much lesser extent on the selection regime. This implies that predicting adaptive genomic response is challenging and requires a good understanding of the adaptive architecture (including linkage disequilibrium) in the ancestral populations. This article is part of the theme issue 'Interdisciplinary approaches to predicting evolutionary biology'.

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.

Towards evolutionary predictions: Current promises and challenges

Evolution has traditionally been a historical and descriptive science, and predicting future evolutionary processes has long been considered impossible. However, evolutionary predictions are increasingly being developed and used in medicine, agriculture, biotechnology and conservation biology. Evolutionary predictions may be used for different purposes, such as to prepare for the future, to try and change the course of evolution or to determine how well we understand evolutionary processes. Similarly, the exact aspect of the evolved population that we want to predict may also differ. For example, we could try to predict which genotype will dominate, the fitness of the population or the extinction probability of a population. In addition, there are many uses of evolutionary predictions that may not always be recognized as such. The main goal of this review is to increase awareness of methods and data in different research fields by showing the breadth of situations in which evolutionary predictions are made. We describe how diverse evolutionary predictions share a common structure described by the predictive scope, time scale and precision. Then, by using examples ranging from SARS-CoV2 and influenza to CRISPR-based gene drives and sustainable product formation in biotechnology, we discuss the methods for predicting evolution, the factors that affect predictability and how predictions can be used to prevent evolution in undesirable directions or to promote beneficial evolution (i.e. evolutionary control). We hope that this review will stimulate collaboration between fields by establishing a common language for evolutionary predictions.

Variation in the genomic basis of parallel phenotypic and ecological divergence in benthic and pelagic morphs of Icelandic Arctic charr (Salvelinus alpinus)

Sympatric adaptive phenotypic divergence should be underlain by genomic differentiation between sub-populations. When divergence drives similar patterns of phenotypic and ecological variation within species we expect evolution to draw on common allelic variation. We investigated divergence histories and genomic signatures of adaptive divergence between benthic and pelagic morphs of Icelandic Arctic charr. Divergence histories for each of four populations were reconstructed using coalescent modelling and 14,187 single nucleotide polymorphisms. Sympatric divergence with continuous gene flow was supported in two populations while allopatric divergence with secondary contact was supported in one population; we could not differentiate between demographic models in the fourth population. We detected parallel patterns of phenotypic divergence along benthic-pelagic evolutionary trajectories among populations. Patterns of genomic differentiation between benthic and pelagic morphs were characterized by outlier loci in many narrow peaks of differentiation throughout the genome, which may reflect the eroding effects of gene flow on nearby neutral loci. We then used genome-wide association analyses to relate both phenotypic (body shape and size) and ecological (carbon and nitrogen stable isotopes) variation to patterns of genomic differentiation. Many peaks of genomic differentiation were associated with phenotypic and ecological variation in the three highly divergent populations, suggesting a genomic basis for adaptive divergence. We detected little evidence for a parallel genomic basis of differentiation as most regions and outlier loci were not shared among populations. Our results show that adaptive divergence can have varied genomic consequences in populations with relatively recent common origins, similar divergence histories, and parallel phenotypic divergence.

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.

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.

Pool-GWAS on reproductive dormancy in Drosophila simulans suggests a polygenic architecture

The genetic basis of adaptation to different environments has been of long-standing interest to evolutionary biologists. Dormancy is a well-studied adaptation to facilitate overwintering. In Drosophila melanogaster, a moderate number of genes with large effects have been described, which suggests a simple genetic basis of dormancy. On the other hand, genome-wide scans for dormancy suggest a polygenic architecture in insects. In D. melanogaster, the analysis of the genetic architecture of dormancy is complicated by the presence of cosmopolitan inversions. Here, we performed a genome-wide scan to characterize the genetic basis of this ecologically extremely important trait in the sibling species of D. melanogaster, D. simulans that lacks cosmopolitan inversions. We performed Pool-GWAS in a South African D. simulans population for dormancy incidence at 2 temperature regimes (10 and 12°C, LD 10:14). We identified several genes with SNPs that showed a significant association with dormancy (P-value < 1e-13), but the overall modest response suggests that dormancy is a polygenic trait with many loci of small effect. Our results shed light on controversies on reproductive dormancy in Drosophila and have important implications for the characterization of the genetic basis of this trait.

Pleiotropy facilitates parallel adaptation in sticklebacks

Highly pleiotropic genes are predicted to be used less often during adaptation, as mutations in these loci are more likely to have negative fitness consequences. Following this logic, we tested whether pleiotropy impacts the probability that a locus will be used repeatedly in adaptation. We used two proxies to estimate pleiotropy: number of phenotypic traits affected by a given genomic region and gene connectivity. We first surveyed 16 independent stream‐lake and three independent benthic‐limnetic ecotype pairs of threespine stickleback to estimate genome‐wide patterns in parallel genomic differentiation. Our analysis revealed parallel divergence across the genome; 30%–37% of outlier regions were shared between at least two independent pairs in either the stream‐lake or benthic‐limnetic comparisons. We then tested whether parallel genomic regions are less pleiotropic than nonparallel regions. Counter to our a priori prediction, parallel genomic regions contained genes with significantly more pleiotropy; that is, influencing a greater number of traits and more highly connected. The increased pleiotropy of parallel regions could not be explained by other genomic factors, as there was no significant difference in mean gene count, mutation or recombination rates between parallel and nonparallel regions. Interestingly, although nonparallel regions contained genes that were less connected and influenced fewer mapped traits on average than parallel regions, they also tended to contain the genes that were predicted to be the most pleiotropic. Taken together, our findings are consistent with the idea that pleiotropy only becomes constraining at high levels and that low or intermediate levels of pleiotropy may be beneficial for adaptation.

Analysis of ancestry heterozygosity suggests that hybrid incompatibilities in threespine stickleback are environment dependent

Hybrid incompatibilities occur when interactions between opposite ancestry alleles at different loci reduce the fitness of hybrids. Most work on incompatibilities has focused on those that are "intrinsic," meaning they affect viability and sterility in the laboratory. Theory predicts that ecological selection can also underlie hybrid incompatibilities, but tests of this hypothesis using sequence data are scarce. In this article, we compiled genetic data for F2 hybrid crosses between divergent populations of threespine stickleback fish (Gasterosteus aculeatus L.) that were born and raised in either the field (seminatural experimental ponds) or the laboratory (aquaria). Because selection against incompatibilities results in elevated ancestry heterozygosity, we tested the prediction that ancestry heterozygosity will be higher in pond-raised fish compared to those raised in aquaria. We found that ancestry heterozygosity was elevated by approximately 3% in crosses raised in ponds compared to those raised in aquaria. Additional analyses support a phenotypic basis for incompatibility and suggest that environment-specific single-locus heterozygote advantage is not the cause of selection on ancestry heterozygosity. Our study provides evidence that, in stickleback, a coarse-albeit indirect-signal of environment-dependent hybrid incompatibility is reliably detectable and suggests that extrinsic incompatibilities can evolve before intrinsic incompatibilities.

Parallel genomic responses to historical climate change and high elevation in East Asian songbirds

Parallel evolution can be expected among closely related taxa exposed to similar selective pressures. However, parallelism is typically stronger at the phenotypic level, while genetic solutions to achieve these phenotypic similarities may differ. For polygenic traits, the availability of standing genetic variation (i.e., heterozygosity) may influence such genetic nonparallelism. Here, we examine the extent to which high-elevation adaptation is parallel-and whether the level of parallelism is affected by heterozygosity-by analyzing genomes of 19 Paridae species distributed across East Asia with a dramatic east-west elevation gradient. We find that western highlands endemic parids have consistently lower levels of heterozygosity-likely the result of late-Pleistocene demographic contraction-than do parids found exclusively in eastern lowlands, which remained unglaciated during the late Pleistocene. Three widespread species (east to west) have high levels of heterozygosity similar to that observed in eastern species, although their western populations are less variable than eastern ones. Comparing genomic responses to extreme environments of the Qinghai-Tibet Plateau, we find that the most differentiated genomic regions between each high-elevation taxon and its low-elevation relative are significantly enriched for genes potentially related to the oxygen transport cascade and/or thermogenesis. Despite no parallelism at particular genes, high similarity in gene function is found among comparisons. Furthermore, parallelism is not higher in more heterozygous widespread parids than in highland endemics. Thus, in East Asian parids, parallel functional response to extreme elevation appears to rely on different genes, with differences in heterozygosity having no effect on the degree of genetic parallelism.

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.

The genetic architecture of temperature adaptation is shaped by population ancestry and not by selection regime

BACKGROUND

Understanding the genetic architecture of temperature adaptation is key for characterizing and predicting the effect of climate change on natural populations. One particularly promising approach is Evolve and Resequence, which combines advantages of experimental evolution such as time series, replicate populations, and controlled environmental conditions, with whole genome sequencing. Recent analysis of replicate populations from two different Drosophila simulans founder populations, which were adapting to the same novel hot environment, uncovered very different architectures-either many selection targets with large heterogeneity among replicates or fewer selection targets with a consistent response among replicates.

RESULTS

Here, we expose the founder population from Portugal to a cold temperature regime. Although almost no selection targets are shared between the hot and cold selection regime, the adaptive architecture was similar. We identify a moderate number of targets under strong selection (19 selection targets, mean selection coefficient = 0.072) and parallel responses in the cold evolved replicates. This similarity across different environments indicates that the adaptive architecture depends more on the ancestry of the founder population than the specific selection regime.

CONCLUSIONS

These observations will have broad implications for the correct interpretation of the genomic responses to a changing climate in natural populations.

Genetic population structure constrains local adaptation in sticklebacks

Repeated and independent adaptation to specific environmental conditions from standing genetic variation is common. However, if genetic variation is limited, the evolution of similar locally adapted traits may be restricted to genetically different and potentially less optimal solutions or prevented from happening altogether. Using a quantitative trait locus (QTL) mapping approach, we identified the genomic regions responsible for the repeated pelvic reduction (PR) in three crosses between nine-spined stickleback populations expressing full and reduced pelvic structures. In one cross, PR mapped to linkage group 7 (LG7) containing the gene Pitx1, known to control pelvic reduction also in the three-spined stickleback. In the two other crosses, PR was polygenic and attributed to 10 novel QTL, of which 90% were unique to specific crosses. When screening the genomes from 27 different populations for deletions in the Pitx1 regulatory element, these were only found in the population in which PR mapped to LG7, even though the morphological data indicated large-effect QTL for PR in several other populations as well. Consistent with the available theory and simulations parameterized on empirical data, we hypothesize that the observed variability in genetic architecture of PR is due to heterogeneity in the spatial distribution of standing genetic variation caused by >2× stronger population structuring among freshwater populations and >10× stronger genetic isolation by distance in the sea in nine-spined sticklebacks as compared to three-spined sticklebacks.

Independent evolution toward larger body size in the distinctive Faroe Island mice

Most phenotypic traits in nature involve the collective action of many genes. Traits that evolve repeatedly are particularly useful for understanding how selection may act on changing trait values. In mice, large body size has evolved repeatedly on islands and under artificial selection in the laboratory. Identifying the loci and genes involved in this process may shed light on the evolution of complex, polygenic traits. Here, we have mapped the genetic basis of body size variation by making a genetic cross between mice from the Faroe Islands, which are among the largest and most distinctive natural populations of mice in the world, and a laboratory mouse strain selected for small body size, SM/J. Using this F2 intercross of 841 animals, we have identified 111 loci controlling various aspects of body size, weight and growth hormone levels. By comparing against other studies, including the use of a joint meta-analysis, we found that the loci involved in the evolution of large size in the Faroese mice were largely independent from those of a different island population or other laboratory strains. We hypothesize that colonization bottleneck, historical hybridization, or the redundancy between multiple loci have resulted in the Faroese mice achieving an outwardly similar phenotype through a distinct evolutionary path.

The genomic signature of ecological divergence along the benthic-limnetic axis in allopatric and sympatric threespine stickleback

The repeated occurrence of similar phenotypes in independent lineages (i.e., parallel evolution) in response to similar ecological conditions can provide compelling insights into the process of adaptive evolution. An intriguing question is to what extent repeated phenotypic changes are underlain by repeated changes at the genomic level and whether patterns of genomic divergence differ with the geographic context in which populations evolve. Here, we combined genomic, morphological and ecological data sets to investigate the genomic signatures of divergence across populations of threespine stickleback (Gasterosteus aculeatus) that adapted to contrasting ecological niches (benthic or limnetic) in either sympatry or allopatry. We found that genome-wide differentiation (F_ST ) was an order of magnitude higher and substantially more repeatable for sympatric benthic and limnetic specialists compared to allopatric populations with similar levels of ecological divergence. We identified genomic regions consistently differentiated between sympatric ecotypes that were also differentiated between or associated with benthic vs. limnetic niche in allopatric populations. These candidate regions were enriched on three chromosomes known to be involved in the benthic-limnetic divergence of threespine stickleback. Some candidate regions overlapped with QTL for body shape and trophic traits such as gill raker number, traits that strongly differ between benthic and limnetic ecotypes. In summary, our study shows that magnitude and repeatability of genomic signatures of ecological divergence in threespine stickleback highly depend on the geographic context. The identified candidate regions provide starting points to identify functionally important genes for the adaptation to benthic and limnetic niches.

Polygenic adaptation: a unifying framework to understand positive selection

Most adaption processes have a polygenic genetic basis, but even with the recent explosive growth of genomic data we are still lacking a unified framework describing the dynamics of selected alleles. Building on recent theoretical and empirical work we introduce the concept of adaptive architecture, which extends the genetic architecture of an adaptive trait by factors influencing its adaptive potential and population genetic principles. Because adaptation can be typically achieved by many different combinations of adaptive alleles (redundancy), we describe how two characteristics - heterogeneity among loci and non-parallelism between replicated populations - are hallmarks for the characterization of polygenic adaptation in evolving populations. We discuss how this unified framework can be applied to natural and experimental populations.

Distinct Patterns of Selective Sweep and Polygenic Adaptation in Evolve and Resequence Studies

In molecular population genetics, adaptation is typically thought to occur via selective sweeps, where targets of selection have independent effects on the phenotype and rise to fixation, whereas in quantitative genetics, many loci contribute to the phenotype and subtle frequency changes occur at many loci during polygenic adaptation. The sweep model makes specific predictions about frequency changes of beneficial alleles and many test statistics have been developed to detect such selection signatures. Despite polygenic adaptation is probably the prevalent mode of adaptation, because of the traditional focus on the phenotype, we are lacking a solid understanding of the similarities and differences of selection signatures under the two models. Recent theoretical and empirical studies have shown that both selective sweep and polygenic adaptation models could result in a sweep-like genomic signature; therefore, additional criteria are needed to distinguish the two models. With replicated populations and time series data, experimental evolution studies have the potential to identify the underlying model of adaptation. Using the framework of experimental evolution, we performed computer simulations to study the pattern of selected alleles for two models: 1) adaptation of a trait via independent beneficial mutations that are conditioned for fixation, that is, selective sweep model and 2) trait optimum model (polygenic adaptation), that is adaptation of a quantitative trait under stabilizing selection after a sudden shift in trait optimum. We identify several distinct patterns of selective sweep and trait optimum models in populations of different sizes. These features could provide the foundation for development of quantitative approaches to differentiate the two models.

Rapid repeatable phenotypic and genomic adaptation following multiple introductions

Uncovering the genomic basis of repeated adaption can provide important insights into the constraints and biases that limit the diversity of genetic responses. Demographic processes such as admixture or bottlenecks affect genetic variation underlying traits experiencing selection. The impact of these processes on the genetic basis of adaptation remains, however, largely unexamined empirically. We here test repeatability in phenotypes and genotypes along parallel climatic clines within the native North American and introduced European and Australian Ambrosia artemisiifolia ranges. To do this, we combined multiple lines of evidence from phenotype-environment associations, F_ST -like outlier tests, genotype-environment associations and genotype-phenotype associations. We used 853 individuals grown in common garden from 84 sampling locations, targeting 19 phenotypes, >83 k SNPs and 22 environmental variables. We found that 17%-26% of loci with adaptive signatures were repeated among ranges, despite alternative demographic histories shaping genetic variation and genetic associations. Our results suggest major adaptive changes can occur on short timescales, with seemingly minimum impacts due to demographic changes linked to introduction. These patterns reveal some predictability of evolutionary change during range expansion, key in a world facing ongoing climate change, and rapid invasive spread.

Parallel genomic architecture underlies repeated sexual signal divergence in Hawaiian Laupala crickets

When the same phenotype evolves repeatedly, we can explore the predictability of genetic changes underlying phenotypic evolution. Theory suggests that genetic parallelism is less likely when phenotypic changes are governed by many small-effect loci compared to few of major effect, because different combinations of genetic changes can result in the same quantitative outcome. However, some genetic trajectories might be favoured over others, making a shared genetic basis to repeated polygenic evolution more likely. To examine this, we studied the genetics of parallel male mating song evolution in the Hawaiian cricket Laupala. We compared quantitative trait loci (QTL) underlying song divergence in three species pairs varying in phenotypic distance. We tested whether replicated song divergence between species involves the same QTL and whether the likelihood of QTL sharing is related to QTL effect size. Contrary to theoretical predictions, we find substantial parallelism in polygenic genetic architectures underlying repeated song divergence. QTL overlapped more frequently than expected based on simulated QTL analyses. Interestingly, QTL effect size did not predict QTL sharing, but did correlate with magnitude of phenotypic divergence. We highlight potential mechanisms driving these constraints on cricket song evolution and discuss a scenario that consolidates empirical quantitative genetic observations with micro-mutational theory.

QTG-Finder: A Machine-Learning Based Algorithm To Prioritize Causal Genes of Quantitative Trait Loci in Arabidopsis and Rice

Linkage mapping is one of the most commonly used methods to identify genetic loci that determine a trait. However, the loci identified by linkage mapping may contain hundreds of candidate genes and require a time-consuming and labor-intensive fine mapping process to find the causal gene controlling the trait. With the availability of a rich assortment of genomic and functional genomic data, it is possible to develop a computational method to facilitate faster identification of causal genes. We developed QTG-Finder, a machine learning based algorithm to prioritize causal genes by ranking genes within a quantitative trait locus (QTL). Two predictive models were trained separately based on known causal genes in Arabidopsis and rice. An independent validation analysis showed that the models could recall about 64% of Arabidopsis and 79% of rice causal genes when the top 20% ranked genes were considered. The top 20% ranked genes can range from 10 to 100 genes, depending on the size of a QTL. The models can prioritize different types of traits though at different efficiency. We also identified several important features of causal genes including paralog copy number, being a transporter, being a transcription factor, and containing SNPs that cause premature stop codon. This work lays the foundation for systematically understanding characteristics of causal genes and establishes a pipeline to predict causal genes based on public data.

Genetic divergence outpaces phenotypic evolution among threespine stickleback populations in old freshwater habitats

Species such as threespine stickleback (Gasterosteus aculeatus) that inhabit divergent selective environments and that have diversified on different time scales can be of value for understanding evolutionary processes. Here we synthesize high-resolution genotypic and phenotypic data to explore a largely unstudied distribution of threespine stickleback populations living in oceanic and freshwater habitats along coastal and inland regions of Oregon. Many inland aquatic habitats of Oregon remained unglaciated during the last ice age, meaning that some extant Oregon lake and river stickleback may have descended from freshwater populations established long before more well-studied, post-glacial freshwater populations. To address the degree of congruence between genetic and phenotypic divergence, we directly compared Oregon stickleback to much younger (post-glacial) Alaskan populations. We found phenotypic variation in Oregon stickleback to be primarily partitioned between oceanic and freshwater habitats, as has been documented in other stickleback systems. However, the main axis of genetic divergence was between coastal and inland regions regardless of habitat type. Furthermore, when comparing patterns between Oregon and Alaska we found similar levels of phenotypic divergence, but much greater genetic divergence among Oregon's populations. The Oregon stickleback system therefore appears well suited for future studies linking genotypic and phenotypic change, further extending the utility of this small fish to provide general insights into evolutionary processes.

A way forward with eco evo devo: an extended theory of resource polymorphism with postglacial fishes as model systems

A major goal of evolutionary science is to understand how biological diversity is generated and altered. Despite considerable advances, we still have limited insight into how phenotypic variation arises and is sorted by natural selection. Here we argue that an integrated view, which merges ecology, evolution and developmental biology (eco evo devo) on an equal footing, is needed to understand the multifaceted role of the environment in simultaneously determining the development of the phenotype and the nature of the selective environment, and how organisms in turn affect the environment through eco evo and eco devo feedbacks. To illustrate the usefulness of an integrated eco evo devo perspective, we connect it with the theory of resource polymorphism (i.e. the phenotypic and genetic diversification that occurs in response to variation in available resources). In so doing, we highlight fishes from recently glaciated freshwater systems as exceptionally well-suited model systems for testing predictions of an eco evo devo framework in studies of diversification. Studies on these fishes show that intraspecific diversity can evolve rapidly, and that this process is jointly facilitated by (i) the availability of diverse environments promoting divergent natural selection; (ii) dynamic developmental processes sensitive to environmental and genetic signals; and (iii) eco evo and eco devo feedbacks influencing the selective and developmental environments of the phenotype. We highlight empirical examples and present a conceptual model for the generation of resource polymorphism - emphasizing eco evo devo, and identify current gaps in knowledge.

John ME. How to Investigate the Origins of Novelty: Insights Gained from Genetic, Behavioral, and Fitness Perspectives

Biologists are drawn to the most extraordinary adaptations in the natural world, often referred to as evolutionary novelties, yet rarely do we understand the microevolutionary context underlying the origins of novel traits, behaviors, or ecological niches. Here we discuss insights gained into the origins of novelty from a research program spanning biological levels of organization from genotype to fitness in Caribbean pupfishes. We focus on a case study of the origins of novel trophic specialists on San Salvador Island, Bahamas and place this radiation in the context of other rapid radiations. We highlight questions that can be addressed about the origins of novelty at different biological levels, such as measuring the isolation of novel phenotypes on the fitness landscape, locating the spatial and temporal origins of adaptive variation contributing to novelty, detecting dysfunctional gene regulation due to adaptive divergence, and connecting behaviors with novel traits. Evolutionary novelties are rare, almost by definition, and we conclude that integrative case studies can provide insights into this rarity relative to the dynamics of adaptation to more common ecological niches and repeated parallel speciation, such as the relative isolation of novel phenotypes on fitness landscapes and the transient availability of ecological, genetic, and behavioral opportunities.

Genomic divergence and adaptive convergence in Drosophila simulans from Evolution Canyon, Israel

Biodiversity refugia formed by unique features of the Mediterranean arid landscape, such as the dramatic ecological contrast of "Evolution Canyon," provide a natural laboratory in which local adaptations to divergent microclimate conditions can be investigated. Significant insights have been provided by studies of Drosophila melanogaster diversifying along the thermal gradient in Evolution Canyon, but a comparative framework to survey adaptive convergence across sister species at the site has been lacking. To fill this void, we present an analysis of genomic polymorphism and evolutionary divergence of Drosophila simulans, a close relative of Drosophila melanogaster with which it co-occurs on both slopes of the canyon. Our results show even deeper interslope divergence in D. simulans than in D. melanogaster, with extensive signatures of selective sweeps present in flies from both slopes but enhanced in the population from the hotter and drier south-facing slope. Interslope divergence was enriched for genes related to electrochemical balance and transmembrane transport, likely in response to increased selection for dehydration resistance on the hotter slope. Both species shared genomic regions that underwent major selective sweeps, but the overall level of adaptive convergence was low, demonstrating no shortage of alternative genomic solutions to cope with the challenges of the microclimate contrast. Mobile elements were a major source of genetic polymorphism and divergence, affecting all parts of the genome, including coding sequences of mating behavior-related genes.

(Non)Parallel Evolution

Parallel evolution across replicate populations has provided evolutionary biologists with iconic examples of adaptation. When multiple populations colonize seemingly similar habitats, they may evolve similar genes, traits, or functions. Yet, replicated evolution in nature or in the laboratory often yields inconsistent outcomes: Some replicate populations evolve along highly similar trajectories, whereas other replicate populations evolve to different extents or in distinct directions. To understand these heterogeneous outcomes, biologists are increasingly treating parallel evolution not as a binary phenomenon but rather as a quantitative continuum ranging from parallel to nonparallel. By measuring replicate populations’ positions along this (non)parallel continuum, we can test hypotheses about evolutionary and ecological factors that influence the extent of repeatable evolution. We review evidence regarding the manifestation of (non)parallel evolution in the laboratory, in natural populations, and in applied contexts such as cancer. We enumerate the many genetic, ecological, and evolutionary processes that contribute to variation in the extent of parallel evolution.

Differential predation alters pigmentation in threespine stickleback (Gasterosteus aculeatus)

Animal pigmentation plays a key role in many biological interactions, including courtship and predator avoidance. Sympatric benthic and limnetic ecotypes of threespine stickleback (Gasterosteus aculeatus) exhibit divergent pigment patterns. To test whether differential predation by cutthroat trout contributes to the differences in pigmentation seen between the ecotypes, we used a within-generation selection experiment on F₂ benthic-limnetic hybrids. After 10 months of differential selection, we compared the pigmentation of fish under trout predation to control fish not exposed to trout predation. We found that stickleback exhibited more lateral barring in ponds with trout predation. Ponds with trout were also less turbid, and a greater degree of barring was negatively correlated with the magnitude of turbidity across pond replicates. A more benthic diet, a proxy for habitat use, was also correlated with greater lateral barring and green dorsal pigmentation. These patterns suggest that differential exposure to cutthroat trout predation may explain the divergence in body pigmentation between benthic and limnetic ecotypes.

Evolution as an ecosystem process: insights from genomics

Evolution is a fundamental ecosystem process. The study of genomic variation of organisms can not only improve our understanding of evolutionary processes, but also of contemporary and future ecosystem dynamics. We argue that integrative research between the fields of genomics and ecosystem ecology could generate new insights. Specifically, studies of biodiversity and ecosystem functioning, evolutionary rescue, and eco-evolutionary dynamics could all benefit from information about variation in genome structure and the genetic architecture of traits, whereas genomic studies could benefit from information about the ecological context of evolutionary dynamics. We propose new ways to help link research on functional genomic diversity with (reciprocal) interactions between phenotypic evolution and ecosystem change. Despite numerous challenges, we anticipate that the wealth of genomic data being collected on natural populations will improve our understanding of ecosystems.

The genetic and molecular architecture of phenotypic diversity in sticklebacks

A major goal of evolutionary biology is to identify the genotypes and phenotypes that underlie adaptation to divergent environments. Stickleback fish, including the threespine stickleback (Gasterosteus aculeatus) and the ninespine stickleback (Pungitius pungitius), have been at the forefront of research to uncover the genetic and molecular architecture that underlies phenotypic diversity and adaptation. A wealth of quantitative trait locus (QTL) mapping studies in sticklebacks have provided insight into long-standing questions about the distribution of effect sizes during adaptation as well as the role of genetic linkage in facilitating adaptation. These QTL mapping studies have also provided a basis for the identification of the genes that underlie phenotypic diversity. These data have revealed that mutations in regulatory elements play an important role in the evolution of phenotypic diversity in sticklebacks. Genetic and molecular studies in sticklebacks have also led to new insights on the genetic basis of repeated evolution and suggest that the same loci are involved about half of the time when the same phenotypes evolve independently. When the same locus is involved, selection on standing variation and repeated mutation of the same genes have both contributed to the evolution of similar phenotypes in independent populations.This article is part of the themed issue 'Evo-devo in the genomics era, and the origins of morphological diversity'.

Genetic Coupling of Female Mate Choice with Polygenic Ecological Divergence Facilitates Stickleback Speciation

Ecological speciation with gene flow is widespread in nature [1], but it presents a conundrum: how are associations between traits under divergent natural selection and traits that contribute to assortative mating maintained? Theoretical models suggest that genetic mechanisms inhibiting free recombination between loci underlying these two types of traits (hereafter, "genetic coupling") can facilitate speciation [2-4]. Here, we perform a direct test for genetic coupling by mapping both divergent traits and female mate choice in a classic model of ecological speciation: sympatric benthic and limnetic threespine stickleback (Gasterosteus aculeatus). By measuring mate choice in F2 hybrid females, we allowed for recombination between loci underlying assortative mating and those under divergent ecological selection. In semi-natural mating arenas in which females had access to both benthic and limnetic males, we found that F2 females mated with males similar to themselves in body size and shape. In addition, we found two quantitative trait loci (QTLs) associated with female mate choice that also predicted female morphology along the benthic-limnetic trait axis. Furthermore, a polygenic genetic model that explains adaptation to contrasting benthic and limnetic feeding niches [5] also predicted F2 female mate choice. Together, these results provide empirical evidence that genetic coupling of assortative mating with traits under divergent ecological selection helps maintain species in the face of gene flow, despite a polygenic basis for adaptation to divergent environments.

Significant Synteny and Colocalization of Ecologically Relevant Quantitative Trait Loci Within and Across Species of Salmonid Fishes

The organization of functional regions within genomes has important implications for evolutionary potential. Considerable research effort has gone toward identifying the genomic basis of phenotypic traits of interest through quantitative trait loci (QTL) analyses. Less research has assessed the arrangement of QTL in the genome within and across species. To investigate the distribution, extent of colocalization, and the synteny of QTL for ecologically relevant traits, we used a comparative genomic mapping approach within and across a range of salmonid species. We compiled 943 QTL from all available species [lake whitefish (Coregonus clupeaformis), coho salmon (Oncorhynchus kisutch), rainbow trout (O. mykiss), Chinook salmon (O. tshawytscha), Atlantic salmon (Salmo salar), and Arctic charr (Salvelinus alpinus)]. We developed a novel analytical framework for mapping and testing the distribution of these QTL. We found no correlation between QTL density and gene density at the chromosome level but did at the fine-scale. Two chromosomes were significantly enriched for QTL. We found multiple synteny blocks for morphological, life history, and physiological traits across species, but only morphology and physiology had significantly more than expected. Two or three pairs of traits were significantly colocalized in three species (lake whitefish, coho salmon, and rainbow trout). Colocalization and fine-scale synteny suggest genetic linkage between traits within species and a conserved genetic basis across species. However, this pattern was weak overall, with colocalization and synteny being relatively rare. These findings advance our understanding of the role of genomic organization in the renowned ecological and phenotypic variability of salmonid fishes.

Adaptive and non-adaptive divergence in a common landscape

Species in a common landscape often face similar selective environments. The capacity of organisms to adapt to these environments may be largely species specific. Quantifying shared and unique adaptive responses across species within landscapes may thus improve our understanding of landscape-moderated biodiversity patterns. Here we test to what extent populations of two coexisting and phylogenetically related fishes—three-spined and nine-spined stickleback—differ in the strength and nature of neutral and adaptive divergence along a salinity gradient. Phenotypic differentiation, neutral genetic differentiation and genomic signatures of adaptation are stronger in the three-spined stickleback. Yet, both species show substantial phenotypic parallelism. In contrast, genomic signatures of adaptation involve different genomic regions, and are thus non-parallel. The relative contribution of spatial and environmental drivers of population divergence in each species reflects different strategies for persistence in the same landscape. These results provide insight in the mechanisms underlying variation in evolutionary versatility and ecological success among species within landscapes.

The three-spined stickleback is a model species for the study of adaptive divergence. Here, Raeymaekers et al. compare how the three-spined stickleback and its relative the nine-spined stickleback vary at the phenotypic and genomic levels in response to the same spatial and environmental drivers.

Novel Candidate Genes Underlying Extreme Trophic Specialization in Caribbean Pupfishes

The genetic changes responsible for evolutionary transitions from generalist to specialist phenotypes are poorly understood. Here we examine the genetic basis of craniofacial traits enabling novel trophic specialization in a sympatric radiation of Cyprinodon pupfishes endemic to San Salvador Island, Bahamas. This recent radiation consists of a generalist species and two novel specialists: a small-jawed "snail-eater" and a large-jawed "scale-eater." We genotyped 12 million single nucleotide polymorphisms (SNPs) by whole-genome resequencing of 37 individuals of all three species from nine populations and integrated genome-wide divergence scans with association mapping to identify divergent regions containing putatively causal SNPs affecting jaw size-the most rapidly diversifying trait in this radiation. A mere 22 fixed variants accompanied extreme ecological divergence between generalist and scale-eater species. We identified 31 regions (20 kb) containing variants fixed between specialists that were significantly associated with variation in jaw size which contained 11 genes annotated for skeletal system effects and 18 novel candidate genes never previously associated with craniofacial phenotypes. Six of these 31 regions showed robust signs of hard selective sweeps after accounting for demographic history. Our data are consistent with predictions based on quantitative genetic models of adaptation, suggesting that the effect sizes of regions influencing jaw phenotypes are positively correlated with distance between fitness peaks on a complex adaptive landscape.

Regulatory Architecture of Gene Expression Variation in the Threespine Stickleback Gasterosteus aculeatus

Much adaptive evolutionary change is underlain by mutational variation in regions of the genome that regulate gene expression rather than in the coding regions of the genes themselves. An understanding of the role of gene expression variation in facilitating local adaptation will be aided by an understanding of underlying regulatory networks. Here, we characterize the genetic architecture of gene expression variation in the threespine stickleback (Gasterosteus aculeatus), an important model in the study of adaptive evolution. We collected transcriptomic and genomic data from 60 half-sib families using an expression microarray and genotyping-by-sequencing, and located expression quantitative trait loci (eQTL) underlying the variation in gene expression in liver tissue using an interval mapping approach. We identified eQTL for several thousand expression traits. Expression was influenced by polymorphism in both cis- and trans-regulatory regions. Trans-eQTL clustered into hotspots. We did not identify master transcriptional regulators in hotspot locations: rather, the presence of hotspots may be driven by complex interactions between multiple transcription factors. One observed hotspot colocated with a QTL recently found to underlie salinity tolerance in the threespine stickleback. However, most other observed hotspots did not colocate with regions of the genome known to be involved in adaptive divergence between marine and freshwater habitats.

The genetic architecture of novel trophic specialists: larger effect sizes are associated with exceptional oral jaw diversification in a pupfish adaptive radiation

The genetic architecture of adaptation is fundamental to understanding the mechanisms and constraints governing diversification. However, most case studies focus on loss of complex traits or parallel speciation in similar environments. It is still unclear how the genetic architecture of these local adaptive processes compares to the architecture of evolutionary transitions contributing to morphological and ecological novelty. Here, we identify quantitative trait loci (QTL) between two trophic specialists in an excellent case study for examining the origins of ecological novelty: a sympatric radiation of pupfishes endemic to San Salvador Island, Bahamas, containing a large-jawed scale-eater and a short-jawed molluscivore with a skeletal nasal protrusion. These specialized niches and trophic traits are unique among over 2000 related species. Measurements of the fitness landscape on San Salvador demonstrate multiple fitness peaks and a larger fitness valley isolating the scale-eater from the putative ancestral intermediate phenotype of the generalist, suggesting that more large-effect QTL should contribute to its unique phenotype. We evaluated this prediction using an F2 intercross between these specialists. We present the first linkage map for pupfishes and detect significant QTL for sex and eight skeletal traits. Large-effect QTL contributed more to enlarged scale-eater jaws than the molluscivore nasal protrusion, consistent with predictions from the adaptive landscape. The microevolutionary genetic architecture of large-effect QTL for oral jaws parallels the exceptional diversification rates of oral jaws within the San Salvador radiation observed over macroevolutionary timescales and may have facilitated exceptional trophic novelty in this system.

Genetics of Skeletal Evolution in Unusually Large Mice from Gough Island

Organisms on islands often undergo rapid morphological evolution, providing a platform for understanding mechanisms of phenotypic change. Many examples of evolution on islands involve the vertebrate skeleton. Although the genetic basis of skeletal variation has been studied in laboratory strains, especially in the house mouse Mus musculus domesticus, the genetic determinants of skeletal evolution in natural populations remain poorly understood. We used house mice living on the remote Gough Island-the largest wild house mice on record-to understand the genetics of rapid skeletal evolution in nature. Compared to a mainland reference strain from the same subspecies (WSB/EiJ), the skeleton of Gough Island mice is considerably larger, with notable expansions of the pelvis and limbs. The Gough Island mouse skeleton also displays changes in shape, including elongations of the skull and the proximal vs. distal elements in the limbs. Quantitative trait locus (QTL) mapping in a large F₂ intercross between Gough Island mice and WSB/EiJ reveals hundreds of QTL that control skeletal dimensions measured at 5, 10, and/or 16 weeks of age. QTL exhibit modest, mostly additive effects, and Gough Island alleles are associated with larger skeletal size at most QTL. The QTL with the largest effects are found on a few chromosomes and affect suites of skeletal traits. Many of these loci also colocalize with QTL for body weight. The high degree of QTL colocalization is consistent with an important contribution of pleiotropy to skeletal evolution. Our results provide a rare portrait of the genetic basis of skeletal evolution in an island population and position the Gough Island mouse as a model system for understanding mechanisms of rapid evolution in nature.

Partially repeatable genetic basis of benthic adaptation in threespine sticklebacks

The extent to which convergent adaptation to similar ecological niches occurs by a predictable genetic basis remains a fundamental question in biology. Threespine stickleback fish have undergone an adaptive radiation in which ancestral oceanic populations repeatedly colonized and adapted to freshwater habitats. In multiple lakes in British Columbia, two different freshwater ecotypes have evolved: a deep-bodied benthic form adapted to forage near the lake substrate, and a narrow-bodied limnetic form adapted to forage in open water. Here, we use genome-wide linkage mapping in marine × benthic F2 genetic crosses to test the extent of shared genomic regions underlying benthic adaptation in three benthic populations. We identify at least 100 Quantitative Trait Loci (QTL) harboring genes influencing skeletal morphology. The majority of QTL (57%) are unique to one cross. However, four genomic regions affecting eight craniofacial and armor phenotypes are found in all three benthic populations. We find that QTL are clustered in the genome and overlapping QTL regions are enriched for genomic signatures of natural selection. These findings suggest that benthic adaptation has occurred via both parallel and nonparallel genetic changes.

Genetic Architecture of Conspicuous Red Ornaments in Female Threespine Stickleback

Explaining the presence of conspicuous female ornaments that take the form of male-typical traits has been a longstanding challenge in evolutionary biology. Such female ornaments have been proposed to evolve via both adaptive and nonadaptive evolutionary processes. Determining the genetic underpinnings of female ornaments is important for elucidating the mechanisms by which such female traits arise and persist in natural populations, but detailed information about their genetic basis is still scarce. In this study, we investigated the genetic architecture of two ornaments, the orange-red throat and pelvic spine, in the threespine stickleback (Gasterosteus aculeatus). Throat coloration is male-specific in ancestral marine populations but has evolved in females in some derived stream populations, whereas sexual dimorphism in pelvic spine coloration is variable among populations. We find that ornaments share a common genetic architecture between the sexes. At least three independent genomic regions contribute to red throat coloration, and harbor candidate genes related to pigment production and pigment cell differentiation. One of these regions is also associated with spine coloration, indicating that both ornaments might be mediated partly via pleiotropic genetic mechanisms.