Genetic signature of adaptive peak shift in threespine stickleback

Genetic Architecture of Flowering Time Differs Between Populations With Contrasting Demographic and Selective Histories

Understanding the evolutionary factors that impact the genetic architecture of traits is a central goal of evolutionary genetics. Here, we investigate how quantitative trait variation accumulated over time in populations that colonized a novel environment. We compare the genetic architecture of flowering time in Arabidopsis populations from the drought-prone Cape Verde Islands and their closest outgroup population from North Africa. We find that trait polygenicity is severely reduced in the island populations compared to the continental North African population. Further, trait architectures and reconstructed allelic histories best fit a model of strong directional selection in the islands in accord with a Fisher-Orr adaptive walk. Consistent with this, we find that large-effect variants that disrupt major flowering time genes (FRI and FLC) arose first, followed by smaller effect variants, including ATX2 L125F, which is associated with a 4-day reduction in flowering time. The most recently arising flowering time-associated loci are not known to be directly involved in flowering time, consistent with an omnigenic signature developing as the population approaches its trait optimum. Surprisingly, we find no effect in the natural population of EDI-Cvi-0 (CRY2 V367M), an allele for which an effect was previously validated by introgression into a Eurasian line. Instead, our results suggest the previously observed effect of the EDI-Cvi-0 allele on flowering time likely depends on genetic background, due to an epistatic interaction. Altogether, our results provide an empirical example of the effects demographic history and selection has on trait architecture.

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.

Species divergence under competition and shared predation

Species competing for resources also commonly share predators. While competition often drives divergence between species, the effects of shared predation are less understood. Theoretically, competing prey species could either diverge or evolve in the same direction under shared predation depending on the strength and symmetry of their interactions. We took an empirical approach to this question, comparing antipredator and trophic phenotypes between sympatric and allopatric populations of threespine stickleback and prickly sculpin fish that all live in the presence of a trout predator. We found divergence in antipredator traits between the species: in sympatry, antipredator adaptations were relatively increased in stickleback but decreased in sculpin. Shifts in feeding morphology, diet and habitat use were also divergent but driven primarily by stickleback evolution. Our results suggest that asymmetric ecological character displacement indirectly made stickleback more and sculpin less vulnerable to shared predation, driving divergence of antipredator traits between sympatric species.

Characterizing phenotypic diversity in marine populations of the threespine stickleback

The threespine stickleback (Gasterosteus aculeatus) is an important model for studying the evolution of vertebrate morphology. Sticklebacks inhabit freshwater, brackish, and marine northern hemisphere waters. Anadromous and marine populations (hereafter marine) are assumed to have remained unchanged morphologically from ancestral marine sticklebacks, despite marine environments varying on regional and local scales. Recent studies suggest that genetic and phenotypic structure exists in marine populations, yet the scale of this variation, and its ecological causes remain unclear. Our goal was to assess morphological trait variation in marine stickleback populations around Southern British Columbia (BC) and determine if oceanographic and habitat characteristics were associated with this variation. Between May-July 2019, we sampled 534 sticklebacks from 15 sites around Vancouver Island, a region characterized by a large diversity of oceanographic and habitat features. We characterized trait variation using two-dimensional (2D) geometric morphometric analysis, comparing individuals between oceanographic regions and habitats. We focused on head and body shape. We found that marine sticklebacks varied morphologically among and between regions and habitats, but the variation did not appear to be related to environmental variation. Sexual dimorphism was the largest source of variation, but oceanographic and habitat variables influenced differences between sexes. We concluded that marine sticklebacks offer abundant opportunities for expanding our knowledge of drivers of morphology.

Genetics of adaptation

The rediscovery of Mendel’s work showing that the heredity of phenotypes is controlled by discrete genes was followed by the reconciliation of Mendelian genetics with evolution by natural selection in the middle of the last century with the Modern Synthesis. In the past two decades, dramatic advances in genomic methods have facilitated the identification of the loci, genes, and even individual mutations that underlie phenotypic variants that are the putative targets of natural selection. Moreover, these methods have also changed how we can study adaptation by flipping the problem around, allowing us to first examine what loci show evidence of having been under selection, and then connecting these genetic variants to phenotypic variation. As a result, we now have an expanding list of actual genetic changes that underlie potentially adaptive phenotypic variation. Here, we synthesize how considering the effects of these adaptive loci in the context of cellular environments, genomes, organisms, and populations has provided new insights to the genetic architecture of adaptation.

Genomic changes underlying repeated niche shifts in an adaptive radiation

In adaptive radiations, single lineages rapidly diversify by adapting to many new niches. Little is known yet about the genomic mechanisms involved, that is, the source of genetic variation or genomic architecture facilitating or constraining adaptive radiation. Here, we investigate genomic changes associated with repeated invasion of many different freshwater niches by threespine stickleback in the Haida Gwaii archipelago, Canada, by resequencing single genomes from one marine and 28 freshwater populations. We find 89 likely targets of parallel selection in the genome that are enriched for old standing genetic variation. In contrast to theoretical expectations, their genomic architecture is highly dispersed with little clustering. Candidate genes and genotype-environment correlations match the three major environmental axes predation regime, light environment, and ecosystem size. In a niche space with these three dimensions, we find that the more divergent a new niche from the ancestral marine habitat, the more loci show signatures of parallel selection. Our findings suggest that the genomic architecture of parallel adaptation in adaptive radiation depends on the steepness of ecological gradients and the dimensionality of the niche space.

The evolution of larger size in high-altitude Drosophila melanogaster has a variable genetic architecture

Important uncertainties persist regarding the genetic architecture of adaptive trait evolution in natural populations, including the number of genetic variants involved, whether they are drawn from standing genetic variation, and whether directional selection drives them to complete fixation. Here, we take advantage of a unique natural population of Drosophila melanogaster from the Ethiopian highlands, which has evolved larger body size than any other known population of this species. We apply a bulk segregant quantitative trait locus mapping approach to 4 unique crosses between highland Ethiopian and lowland Zambian populations for both thorax length and wing length. Results indicated a persistently variable genetic basis for these evolved traits (with largely distinct sets of quantitative trait loci for each cross), and at least a moderately polygenic architecture with relatively strong effects present. We complemented these mapping experiments with population genetic analyses of quantitative trait locus regions and gene ontology enrichment analysis, generating strong hypotheses for specific genes and functional processes that may have contributed to these adaptive trait changes. Finally, we find that the genetic architectures indicated by our quantitative trait locus mapping results for size traits mirror those from similar experiments on other recently evolved traits in this species. Collectively, these studies suggest a recurring pattern of polygenic adaptation in this species, in which causative variants do not approach fixation and moderately strong effect loci are present.

Adaptive divergence and the evolution of hybrid trait mismatch in threespine stickleback

Selection against mismatched traits in hybrids is the phenotypic analogue of intrinsic hybrid incompatibilities. Mismatch occurs when hybrids resemble one parent population for some phenotypic traits and the other parent population for other traits, and is caused by dominance in opposing directions or from segregation of alleles in recombinant hybrids. In this study, we used threespine stickleback fish (Gasterosteus aculeatus L.) to test the theoretical prediction that trait mismatch in hybrids should increase with the magnitude of phenotypic divergence between parent populations. We measured morphological traits in parents and hybrids in crosses between a marine population representing the ancestral form and twelve freshwater populations that have diverged from this ancestral state to varying degrees according to their environments. We found that trait mismatch was greater in more divergent crosses for both F₁ and F₂ hybrids. In the F₁, the divergence–mismatch relationship was caused by traits having dominance in different directions, whereas it was caused by increasing segregating phenotypic variation in the F₂. Our results imply that extrinsic hybrid incompatibilities accumulate as phenotypic divergence proceeds.

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.

Allen Orr and the genetics of adaptation

Over most of the 20th century, evolutionary biologists predominantly subscribed to a strong form of "micro-mutationism," in which adaptive phenotypic divergence arises from allele frequency changes at many loci, each with a small effect on the phenotype. To be sure, there were well-known examples of large-effect alleles contributing to adaptation, yet such cases were generally regarded as atypical and unrepresentative of evolutionary change in general. In 1998, Allen Orr published a landmark theoretical paper in Evolution, which showed that both small- and large-effect mutations are likely to contribute to "adaptive walks" of a population to an optimum. Coupled with a growing set of empirical examples of large-effect alleles contributing to divergence (e.g., from QTL studies), Orr's paper provided a mathematical formalism that converted many evolutionary biologists from micro-mutationism to a more pluralistic perspective on the genetic basis of evolutionary change. We revisit the theoretical insights emerging from Orr's paper within the historical context leading up to 1998, and track the influence of this paper on the field of evolutionary biology through an examination of its citations over the last two decades and an analysis of the extensive body of theoretical and empirical research that Orr's pioneering paper inspired.

Longer or shorter spines: Reciprocal trait evolution in stickleback via triallelic regulatory changes in Stanniocalcin2a

Vertebrates have repeatedly modified skeletal structures to adapt to their environments. The threespine stickleback is an excellent system for studying skeletal modifications, as different wild populations have either increased or decreased the lengths of their prominent dorsal and pelvic spines in different freshwater environments. Here we identify a regulatory locus that has a major morphological effect on the length of stickleback dorsal and pelvic spines, which we term Maser (major spine enhancer). Maser maps in a closely linked supergene complex that controls multiple armor, feeding, and behavioral traits on chromosome IV. Natural alleles in Maser are differentiated between marine and freshwater sticklebacks; however, alleles found among freshwater populations are also differentiated, with distinct alleles found in short- and long-spined freshwater populations. The distinct freshwater alleles either increase or decrease expression of the bone growth inhibitor gene Stanniocalcin2a in developing spines, providing a simple genetic mechanism for either increasing or decreasing spine lengths in natural populations. Genomic surveys suggest many recurrently differentiated loci in sticklebacks are similarly specialized into three or more distinct alleles, providing multiple ancient standing variants in particular genes that may contribute to a range of phenotypes in different environments.

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

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

Genome-Wide Association Study of Body Shape-Related Traits in Large Yellow Croaker (Larimichthys crocea)

Large yellow croaker (Larimichthys crocea) is one of the most important cultured marine fish on the southeast coast of China. Its body shape is important for the aquaculture industry since it affects the behavior such as swimming, ingesting, and evading, as well as customer preference. Due to the greater consumer demand of small head, slender body large yellow croaker, selecting and breeding of slender individuals with the assistance of genetic markers will benefit the industry quickly. In this study, several traits were employed to represent body shape, including body depth/body length (BD/BL), body thickness/body length (BT/BL), caudal peduncle depth/caudal peduncle length (CPDLR), tail length/body length (TL/BL), and body area/head area (BA/HA). Genome-wide association study was conducted with a panmictic population of 280 individuals to identify SNP and genes potentially associated with body shape. A set of 20 SNPs on 12 chromosomes were identified to be significantly associated with body shape-related traits. Besides, 5 SNPs were identified to be suggestive associated with CPDLR and BT/BL. Surrounding these SNPs, we found some body shape-related candidate genes, including fabp1, acrv1, bcor, mstn, bambi, and neo1, which involved in lipid metabolism, TGF-β signaling, and BMP pathway and other important regulatory pathways. These results will be useful for the understanding of the genetic basis of body shape formation and helpful for body shape controlling of large yellow croaker by using marker-assisted selection.

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

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

A Single Interacting Species Leads to Widespread Parallel Evolution of the Stickleback Genome

Biotic interactions are potent, widespread causes of natural selection and divergent phenotypic evolution and can lead to genetic differentiation with gene flow among wild populations ("isolation by ecology") [1-4]. Biotic selection has been predicted to act on more genes than abiotic selection thereby driving greater adaptation [5]. However, difficulties in isolating the genome-wide effect of single biotic agents of selection have limited our ability to identify and quantify the number and type of genetic regions responding to biotic selection [6-9]. We identified geographically interspersed lakes in which threespine stickleback fish (Gasterosteus aculeatus) have repeatedly adapted to the presence or absence of a single member of the ecological community, prickly sculpin (Cottus asper), a fish that is both a competitor and a predator of the stickleback [10]. Whole-genome sequencing revealed that sculpin presence or absence accounted for the majority of genetic divergence among stickleback populations, more so than geography. The major axis of genomic variation within and between the two lake types was correlated with multiple traits, indicating parallel natural selection across a gradient of biotic environments. A large proportion of the genome-about 1.8%, encompassing more than 600 genes-differentiated stickleback from the two biotic environments. Divergence occurred in 141 discrete genomic clumps located mainly in regions of low recombination, suggesting that genes brought to lakes by the colonizing ancestral population often evolved together in linked blocks. Strong selection and a wealth of standing genetic variation explain how a single member of the biotic community can have such a rapid and profound evolutionary impact.

Complexity Theory: An Overview with Potential Applications for the Social Sciences

Systems theory has been challenged in the recent literature due to its perceived disconnection from today’s research and practice demands. Moving away from the reductionist frameworks and the complicated domain predominated by known unknowns and order, a call is being made to the social sciences to begin adopting complexity theory and newer connectionist methods that better address complexity and open social systems. Scholars and scholar-practitioners will continue to find the need to apply complexity theory as wicked problems become more prevalent in the social sciences. This paper differentiates between general systems theory (GST) and complexity theory, as well as identifies advantages for the social sciences in incorporating complexity theory as a formal theory. Complexity theory is expanded upon and identified as providing a new perspective and a new method of theorizing that can be practiced by disciplines within the social sciences. These additions could better position the social sciences to address the complexity associated with advancing technology, globalization, intricate markets, cultural change, and the myriad of challenges and opportunities to come.

Efficient CRISPR-Cas9 editing of major evolutionary loci in sticklebacks

BACKGROUND

Stickleback fish are widely used to study the genetic and ecological basis of phenotypic evolution. Although several major loci have now been identified that contribute to evolutionary differences between wild populations, further study of the phenotypes associated with particular genes and mutations has been limited by the difficulty of generating targeted mutations at precise locations in the stickleback genome.

APPROACH AND AIMS

We compared different methods of expressing single-guide RNAs (sgRNAs) and Cas9 activity in fertilized stickleback eggs. We used an easily scored pigmentation gene (SLC24A5) to screen for molecular lesions, phenotypic effects, and possible germline transmission of newly induced alleles. We then used the optimized CRISPR methods to target two major evolutionary loci in sticklebacks, KITLG and EDA. We hypothesized that coding region mutations in the KITLG gene would alter body pigmentation and possibly sex determination, and that mutations in the EDA gene would disrupt the formation of most armor plates, fin rays, spines, teeth, and gill rakers.

RESULTS

Targeted deletions were successfully induced at each target locus by co-injecting one-cell stage stickleback embryos with either Cas9 mRNA or Cas9 protein, together with sgRNAs designed to protein-coding exons. Founder animals were typically mosaic for multiple mutations, which they transmitted through the germline at overall rates of 21 to 100%. We found that the copy of KITLG on the X chromosome (KITLGX) has diverged from the KITLG on the Y chromosome (KITLGY). Predicted loss-of-function mutations in the KITLGX gene dramatically altered pigmentation in both external skin and internal organ, but the same was not true for KITLGY mutations. Predicted loss-of-function mutations in either the KITLGX or KITLGY genes did not lead to sex reversal or prevent fertility. Homozygous loss-of-function mutations in the EDA gene led to complete loss of armor plates, severe reduction or loss of most soft rays in the dorsal, anal, and caudal fins, and severe reductions in tooth and gill raker number. In contrast, long dorsal and pelvic spines remained intact in EDA mutant animals, suggesting that common co-segregation of plate loss and spine reduction in wild populations is unlikely to be due to pleiotropic effects of EDA mutations.

CONCLUSION

CRISPR-Cas9 approaches can be used to induce germline mutations in key evolutionary loci in sticklebacks. Targeted coding region mutations confirm an important role for KITLG and EDA in skin pigmentation and armor plate reduction, respectively. They also provide new information about the functions of these genes in other body structures.

RADseq and mate choice assays reveal unidirectional gene flow among three lamprey ecotypes despite weak assortative mating: Insights into the formation and stability of multiple ecotypes in sympatry

Adaptive divergence with gene flow often results in complex patterns of variation within taxa exhibiting substantial ecological differences among populations. One example where this may have occurred is the parallel evolution of freshwater-resident nonparasitic lampreys from anadromous-parasitic ancestors. Previous studies have focused on transitions between these two phenotypic extremes, but here, we considered more complex evolutionary scenarios where an intermediate freshwater form that remains parasitic is found sympatrically with the other two ecotypes. Using population genomic analysis (restriction-associated DNA sequencing), we found that a freshwater-parasitic ecotype was highly distinct from an anadromous-parasitic form (Q_lake-P  = 96.8%, F_st  = 0.154), but that a freshwater-nonparasitic form was almost completely admixed in Loch Lomond, Scotland. Demographic reconstructions indicated that both freshwater populations likely derived from a common freshwater ancestor. However, while the nonparasitic ecotype has experienced high levels of introgression from the anadromous-parasitic ecotype (Q_anad-P  = 37.7%), there is no evidence of introgression into the freshwater-parasitic ecotype. Paradoxically, mate choice experiments predicted high potential for gene flow: Males from all ecotypes were stimulated to spawn with freshwater-parasitic females, which released gametes in response to all ecotypes. Differentially fixed single nucleotide polymorphisms identified genes associated with growth and development, which could possibly influence the timing of metamorphosis, resulting in significant ecological differences between forms. This suggests that multiple lamprey ecotypes can persist in sympatry following shifts in adaptive peaks, due to environmental change during their repeated colonization of post-glacial regions, followed by periods of extensive gene flow among such diverging populations.

Keystone Genes

The keystone species concept is used in ecology to describe individual species with disproportionately large effects on their communities. We extend this idea to the level of genes with disproportionately large effects on ecological processes. Such 'keystone genes' (KGs) would underlie traits involved in species interactions or causing critical biotic and/or abiotic changes that influence emergent community and ecosystem properties. We propose a general framework for how KGs could be identified, while keeping KGs under the umbrella of 'ecologically important genes' (EIGs) that also include categories such as 'foundation genes', 'ecosystem engineering genes', and more. Although likely rare, KGs and other EIGs could dominate certain ecological processes; thus, their discovery and study are relevant for understanding eco-evolutionary dynamics.

Contemporary ancestor? Adaptive divergence from standing genetic variation in Pacific marine threespine stickleback

Background

Populations that have repeatedly colonized novel environments are useful for studying the role of ecology in adaptive divergence – particularly if some individuals persist in the ancestral habitat. Such “contemporary ancestors” can be used to demonstrate the effects of selection by comparing phenotypic and genetic divergence between the derived population and their extant ancestors. However, evolution and demography in these “contemporary ancestors” can complicate inferences about the source (standing genetic variation, de novo mutation) and pace of adaptive divergence. Marine threespine stickleback (Gasterosteus aculeatus) have colonized freshwater environments along the Pacific coast of North America, but have also persisted in the marine environment. To what extent are marine stickleback good proxies of the ancestral condition?

Results

We sequenced > 5800 variant loci in over 250 marine stickleback from eight locations extending from Alaska to California, and phenotyped them for platedness and body shape. Pairwise F_ST varied from 0.02 to 0.18. Stickleback were divided into five genetic clusters, with a single cluster comprising stickleback from Washington to Alaska. Plate number, Eda, body shape, and candidate loci showed evidence of being under selection in the marine environment. Comparisons to a freshwater population demonstrated that candidate loci for freshwater adaptation varied depending on the choice of marine populations.

Conclusions

Marine stickleback are structured into phenotypically and genetically distinct populations that have been evolving as freshwater stickleback evolved. This variation complicates their usefulness as proxies of the ancestors of freshwater populations. Lessons from stickleback may be applied to other “contemporary ancestor”-derived population studies.

Electronic supplementary material

The online version of this article (10.1186/s12862-018-1228-8) contains supplementary material, which is available to authorized users.

Natural selection and the predictability of evolution in Timema stick insects

Predicting evolution remains difficult. We studied the evolution of cryptic body coloration and pattern in a stick insect using 25 years of field data, experiments, and genomics. We found that evolution is more difficult to predict when it involves a balance between multiple selective factors and uncertainty in environmental conditions than when it involves feedback loops that cause consistent back-and-forth fluctuations. Specifically, changes in color-morph frequencies are modestly predictable through time (r² = 0.14) and driven by complex selective regimes and yearly fluctuations in climate. In contrast, temporal changes in pattern-morph frequencies are highly predictable due to negative frequency-dependent selection (r² = 0.86). For both traits, however, natural selection drives evolution around a dynamic equilibrium, providing some predictability to the process.

Speciation through the lens of biomechanics: locomotion, prey capture and reproductive isolation

Speciation is a multifaceted process that involves numerous aspects of the biological sciences and occurs for multiple reasons. Ecology plays a major role, including both abiotic and biotic factors. Whether populations experience similar or divergent ecological environments, they often adapt to local conditions through divergence in biomechanical traits. We investigate the role of biomechanics in speciation using fish predator-prey interactions, a primary driver of fitness for both predators and prey. We highlight specific groups of fishes, or specific species, that have been particularly valuable for understanding these dynamic interactions and offer the best opportunities for future studies that link genetic architecture to biomechanics and reproductive isolation (RI). In addition to emphasizing the key biomechanical techniques that will be instrumental, we also propose that the movement towards linking biomechanics and speciation will include (i) establishing the genetic basis of biomechanical traits, (ii) testing whether similar and divergent selection lead to biomechanical divergence, and (iii) testing whether/how biomechanical traits affect RI. Future investigations that examine speciation through the lens of biomechanics will propel our understanding of this key process.

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'.

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.

Genetic dissection of adaptive form and function in rapidly speciating cichlid fishes

Genes of major phenotypic effects and strong genetic correlations can facilitate adaptation, direct selective responses, and potentially lead to phenotypic convergence. However, the preponderance of this type of genetic architecture in repeatedly evolved adaptations remains unknown. Using hybrids between Haplochromis chilotes (thick-lipped) and Pundamilia nyererei (thin-lipped) we investigated the genetics underlying hypertrophied lips and elongated heads, traits that evolved repeatedly in cichlids. At least 25 loci of small-to-moderate and mainly additive effects were detected. Phenotypic variation in lip and head morphology was largely independent. Although several QTL overlapped for lip and head morphology traits, they were often of opposite effects. The distribution of effect signs suggests strong selection on lips. The fitness implications of several detected loci were demonstrated using a laboratory assay testing for the association between genotype and variation in foraging performance. The persistence of low fitness alleles in head morphology appears to be maintained through antagonistic pleiotropy/close linkage with positive-effect lip morphology alleles. Rather than being based on few major loci with strong positive genetic correlations, our results indicate that the evolution of the Lake Victoria thick-lipped ecomorph is the result of selection on numerous loci distributed throughout the genome.

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.

Characterizing phenotypic divergence using three-dimensional geometric morphometrics in four populations of threespine stickleback (Gasterosteus aculeatus; Pisces: Gasterosteidae) in Katmai National Park and Preserve, Alaska

The threespine stickleback (Gasterosteus aculeatus L., 1758) is a vertebrate model for the study of the relationship between phenotype and environment in facilitating rapid evolutionary change. Using four populations from a system of lakes in Katmai National Park and Preserve, Alaska, and microcomputed tomography and three-dimensional geometric morphometrics, we test the hypothesis that stickleback populations inhabiting freshwater environments display cranial phenotypes that are intermediate between the putative ancestral form and the low-plated freshwater populations that demonstrate substantial divergence toward new phenotypic optima. We further test the hypothesis that phenotypic covariance structure is disrupted in the context of such putatively recent adaptive events. We report significant phenotypic differences among all four populations that includes a component of sexual dimorphism. Furthermore, we show evidence of disrupted phenotypic covariance structure among these populations. Taken together, these findings indicate the importance of phenotypic quantification as a key step in elucidating both the ecological processes responsible for rapid adaptive radiations and the role of developmental mechanisms in biasing evolutionary change.

Quantitative trait locus analysis of body shape divergence in nine-spined sticklebacks based on high-density SNP-panel

Heritable phenotypic differences between populations, caused by the selective effects of distinct environmental conditions, are of commonplace occurrence in nature. However, the actual genomic targets of this kind of selection are still poorly understood. We conducted a quantitative trait locus (QTL) mapping study to identify genomic regions responsible for morphometric differentiation between genetically and phenotypically divergent marine and freshwater nine-spined stickleback (Pungitius pungitius) populations. Using a dense panel of SNP-markers obtained by restriction site associated DNA sequencing of an F2 recombinant cross, we found 22 QTL that explained 3.5-12.9% of phenotypic variance in the traits under investigation. We detected one fairly large-effect (PVE = 9.6%) QTL for caudal peduncle length-a trait with a well-established adaptive function showing clear differentiation among marine and freshwater populations. We also identified two large-effect QTL for lateral plate numbers, which are different from the lateral plate QTL reported in earlier studies of this and related species. Hence, apart from identifying several large-effect QTL in shape traits showing adaptive differentiation in response to different environmental conditions, the results suggest intra- and interspecific heterogeneity in the genomic basis of lateral plate number variation.

Factors influencing the effect size distribution of adaptive substitutions

The distribution of effect sizes of adaptive substitutions has been central to evolutionary biology since the modern synthesis. Early theory proposed that because large-effect mutations have negative pleiotropic consequences, only small-effect mutations contribute to adaptation. More recent theory suggested instead that large-effect mutations could be favoured when populations are far from their adaptive peak. Here we suggest that the distributions of effect sizes are expected to differ among study systems, reflecting the wide variation in evolutionary forces and ecological conditions experienced in nature. These include selection, mutation, genetic drift, gene flow, and other factors such as the degree of pleiotropy, the distance to the phenotypic optimum, whether the optimum is stable or moving, and whether new mutation or standing genetic variation provides the source of adaptive alleles. Our goal is to review how these factors might affect the distribution of effect sizes and to identify new research directions. Until more theory and empirical work is available, we feel that it is premature to make broad generalizations about the effect size distribution of adaptive substitutions important in nature.

Genomics of Rapid Incipient Speciation in Sympatric Threespine Stickleback

Ecological speciation is the process by which reproductively isolated populations emerge as a consequence of divergent natural or ecologically-mediated sexual selection. Most genomic studies of ecological speciation have investigated allopatric populations, making it difficult to infer reproductive isolation. The few studies on sympatric ecotypes have focused on advanced stages of the speciation process after thousands of generations of divergence. As a consequence, we still do not know what genomic signatures of the early onset of ecological speciation look like. Here, we examined genomic differentiation among migratory lake and resident stream ecotypes of threespine stickleback reproducing in sympatry in one stream, and in parapatry in another stream. Importantly, these ecotypes started diverging less than 150 years ago. We obtained 34,756 SNPs with restriction-site associated DNA sequencing and identified genomic islands of differentiation using a Hidden Markov Model approach. Consistent with incipient ecological speciation, we found significant genomic differentiation between ecotypes both in sympatry and parapatry. Of 19 islands of differentiation resisting gene flow in sympatry, all were also differentiated in parapatry and were thus likely driven by divergent selection among habitats. These islands clustered in quantitative trait loci controlling divergent traits among the ecotypes, many of them concentrated in one region with low to intermediate recombination. Our findings suggest that adaptive genomic differentiation at many genetic loci can arise and persist in sympatry at the very early stage of ecotype divergence, and that the genomic architecture of adaptation may facilitate this.

Intraguild predation leads to genetically based character shifts in the threespine stickleback

Intraguild predation is a common ecological interaction that occurs when a species preys upon another species with which it competes. The interaction is potentially a mechanism of divergence between intraguild prey (IG-prey) populations, but it is unknown if cases of character shifts in IG-prey are an environmental or evolutionary response. We investigated the genetic basis and inducibility of character shifts in threespine stickleback from lakes with and without prickly sculpin, a benthic intraguild predator (IG-predator). Wild populations of stickleback sympatric with sculpin repeatedly show greater defensive armor and water column height preference. We laboratory-raised stickleback from lakes with and without sculpin, as well as marine stickleback, and found that differences between populations in armor, body shape, and behavior persisted in a common garden. Within the common garden, we raised stickleback half-families from multiple populations in the presence and absence of sculpin. Although the presence of sculpin induced trait changes in the marine stickleback, we did not observe an induced response in the freshwater stickleback. Behavioral and morphological trait differences between freshwater populations thus have a genetic basis and suggest an evolutionary response to intraguild predation.

Fundamental insights into ontogenetic growth from theory and fish

The fundamental features of growth may be universal, because growth trajectories of most animals are very similar, but a unified mechanistic theory of growth remains elusive. Still needed is a synthetic explanation for how and why growth rates vary as body size changes, both within individuals over their ontogeny and between populations and species over their evolution. Here, we use Bertalanffy growth equations to characterize growth of ray-finned fishes in terms of two parameters, the growth rate coefficient, K, and final body mass, m∞. We derive two alternative empirically testable hypotheses and test them by analyzing data from FishBase. Across 576 species, which vary in size at maturity by almost nine orders of magnitude, K scaled as [Formula: see text]. This supports our first hypothesis that growth rate scales as [Formula: see text] as predicted by metabolic scaling theory; it implies that species that grow to larger mature sizes grow faster as juveniles. Within fish species, however, K scaled as [Formula: see text]. This supports our second hypothesis, which predicts that growth rate scales as [Formula: see text] when all juveniles grow at the same rate. The unexpected disparity between across- and within-species scaling challenges existing theoretical interpretations. We suggest that the similar ontogenetic programs of closely related populations constrain growth to [Formula: see text] scaling, but as species diverge over evolutionary time they evolve the near-optimal [Formula: see text] scaling predicted by metabolic scaling theory. Our findings have important practical implications because fish supply essential protein in human diets, and sustainable yields from wild harvests and aquaculture depend on growth rates.

Genomic architecture of habitat-related divergence and signature of directional selection in the body shapes of Gnathopogon fishes

Evolution of ecomorphologically relevant traits such as body shapes is important to colonize and persist in a novel environment. Habitat-related adaptive divergence of these traits is therefore common among animals. We studied the genomic architecture of habitat-related divergence in the body shape of Gnathopogon fishes, a novel example of lake-stream ecomorphological divergence, and tested for the action of directional selection on body shape differentiation. Compared to stream-dwelling Gnathopogon elongatus, the sister species Gnathopogon caerulescens, exclusively inhabiting a large ancient lake, had an elongated body, increased proportion of the caudal region and small head, which would be advantageous in the limnetic environment. Using an F2 interspecific cross between the two Gnathopogon species (195 individuals), quantitative trait locus (QTL) analysis with geometric morphometric quantification of body shape and restriction-site associated DNA sequencing-derived markers (1622 loci) identified 26 significant QTLs associated with the interspecific differences of body shape-related traits. These QTLs had small to moderate effects, supporting polygenic inheritance of the body shape-related traits. Each QTL was mostly located on different genomic regions, while colocalized QTLs were detected for some ecomorphologically relevant traits that are proxy of body and caudal peduncle depths, suggesting different degree of modularity among traits. The directions of the body shape QTLs were mostly consistent with the interspecific difference, and QTL sign test suggested a genetic signature of directional selection in the body shape divergence. Thus, we successfully elucidated the genomic architecture underlying the adaptive changes of the quantitative and complex morphological trait in a novel system.

Genetic architecture of pollination syndrome transition between hummingbird-specialist and generalist species in the genus Rhytidophyllum (Gesneriaceae)

Adaptation to pollinators is a key factor of diversification in angiosperms. The Caribbean sister genera Rhytidophyllum and Gesneria present an important diversification of floral characters. Most of their species can be divided in two major pollination syndromes. Large-open flowers with pale colours and great amount of nectar represent the generalist syndrome, while the hummingbird-specialist syndrome corresponds to red tubular flowers with a less important nectar volume. Repeated convergent evolution toward the generalist syndrome in this group suggests that such transitions rely on few genes of moderate to large effect. To test this hypothesis, we built a linkage map and performed a QTL detection for divergent pollination syndrome traits by crossing one specimen of the generalist species Rhytidophyllum auriculatum with one specimen of the hummingbird pollinated R. rupincola. Using geometric morphometrics and univariate traits measurements, we found that floral shape among the second-generation hybrids is correlated with morphological variation observed between generalist and hummingbird-specialist species at the genus level. The QTL analysis showed that colour and nectar volume variation between syndromes involve each one major QTL while floral shape has a more complex genetic basis and rely on few genes of moderate effect. Finally, we did not detect any genetic linkage between the QTLs underlying those traits. This genetic independence of traits could have facilitated evolution toward optimal syndromes.

RAD-QTL Mapping Reveals Both Genome-Level Parallelism and Different Genetic Architecture Underlying the Evolution of Body Shape in Lake Whitefish (Coregonus clupeaformis) Species Pairs

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

Comparative analysis of Japanese three-spined stickleback clades reveals the Pacific Ocean lineage has adapted to freshwater environments while the Japan Sea has not

Divergent selection and adaptive divergence can increase phenotypic diversification amongst populations and lineages. Yet adaptive divergence between different environments, habitats or niches does not occur in all lineages. For example, the colonization of freshwater environments by ancestral marine species has triggered adaptive radiation and phenotypic diversification in some taxa but not in others. Studying closely related lineages differing in their ability to diversify is an excellent means of understanding the factors promoting and constraining adaptive evolution. A well-known example of the evolution of increased phenotypic diversification following freshwater colonization is the three-spined stickleback. Two closely related stickleback lineages, the Pacific Ocean and the Japan Sea occur in Japan. However, Japanese freshwater stickleback populations are derived from the Pacific Ocean lineage only, suggesting the Japan Sea lineage is unable to colonize freshwater. Using stable isotope data and trophic morphology, we first show higher rates of phenotypic and ecological diversification between marine and freshwater populations within the Pacific Ocean lineage, confirming adaptive divergence has occurred between the two lineages and within the Pacific Ocean lineage but not in the Japan Sea lineage. We further identified consistent divergence in diet and foraging behaviour between marine forms from each lineage, confirming Pacific Ocean marine sticklebacks, from which all Japanese freshwater populations are derived, are better adapted to freshwater environments than Japan Sea sticklebacks. We suggest adaptive divergence between ancestral marine populations may have played a role in constraining phenotypic diversification and adaptive evolution in Japanese sticklebacks.

Evolved tooth gain in sticklebacks is associated with a cis-regulatory allele of Bmp6

Developmental genetic studies of evolved differences in morphology have led to the hypothesis that cis-regulatory changes often underlie morphological evolution. However, because most of these studies focus on evolved loss of traits, the genetic architecture and possible association with cis-regulatory changes of gain traits are less understood. Here we show that a derived benthic freshwater stickleback population has evolved an approximate twofold gain in ventral pharyngeal tooth number compared with their ancestral marine counterparts. Comparing laboratory-reared developmental time courses of a low-toothed marine population and this high-toothed benthic population reveals that increases in tooth number and tooth plate area and decreases in tooth spacing arise at late juvenile stages. Genome-wide linkage mapping identifies largely separate sets of quantitative trait loci affecting different aspects of dental patterning. One large-effect quantitative trait locus controlling tooth number fine-maps to a genomic region containing an excellent candidate gene, Bone morphogenetic protein 6 (Bmp6). Stickleback Bmp6 is expressed in developing teeth, and no coding changes are found between the high- and low-toothed populations. However, quantitative allele-specific expression assays of Bmp6 in developing teeth in F1 hybrids show that cis-regulatory changes have elevated the relative expression level of the freshwater benthic Bmp6 allele at late, but not early, stages of stickleback development. Collectively, our data support a model where a late-acting cis-regulatory up-regulation of Bmp6 expression underlies a significant increase in tooth number in derived benthic sticklebacks.

Fisher's geometric model with a moving optimum

Fisher's geometric model has been widely used to study the effects of pleiotropy and organismic complexity on phenotypic adaptation. Here, we study a version of Fisher's model in which a population adapts to a gradually moving optimum. Key parameters are the rate of environmental change, the dimensionality of phenotype space, and the patterns of mutational and selectional correlations. We focus on the distribution of adaptive substitutions, that is, the multivariate distribution of the phenotypic effects of fixed beneficial mutations. Our main results are based on an “adaptive-walk approximation,” which is checked against individual-based simulations. We find that (1) the distribution of adaptive substitutions is strongly affected by the ecological dynamics and largely depends on a single composite parameter γ, which scales the rate of environmental change by the “adaptive potential” of the population; (2) the distribution of adaptive substitution reflects the shape of the fitness landscape if the environment changes slowly, whereas it mirrors the distribution of new mutations if the environment changes fast; (3) in contrast to classical models of adaptation assuming a constant optimum, with a moving optimum, more complex organisms evolve via larger adaptive steps.

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

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