The genetics of adaptive shape shift in stickleback: pleiotropy and effect size

Evolved Bmp6 enhancer alleles drive spatial shifts in gene expression during tooth development in sticklebacks

Mutations in enhancers have been shown to often underlie natural variation but the evolved differences in enhancer activity can be difficult to identify in vivo. Threespine sticklebacks (Gasterosteus aculeatus) are a robust system for studying enhancer evolution due to abundant natural genetic variation, a diversity of evolved phenotypes between ancestral marine and derived freshwater forms, and the tractability of transgenic techniques. Previous work identified a series of polymorphisms within an intronic enhancer of the Bone morphogenetic protein 6 (Bmp6) gene that are associated with evolved tooth gain, a derived increase in freshwater tooth number that arises late in development. Here, we use a bicistronic reporter construct containing a genetic insulator and a pair of reciprocal two-color transgenic reporter lines to compare enhancer activity of marine and freshwater alleles of this enhancer. In older fish, the two alleles drive partially overlapping expression in both mesenchyme and epithelium of developing teeth, but the freshwater enhancer drives a reduced mesenchymal domain and a larger epithelial domain relative to the marine enhancer. In younger fish, these spatial shifts in enhancer activity are less pronounced. Comparing Bmp6 expression by in situ hybridization in developing teeth of marine and freshwater fish reveals similar evolved spatial shifts in gene expression. Together, these data support a model in which the polymorphisms within this enhancer underlie evolved tooth gain by shifting the spatial expression of Bmp6 during tooth development, and provide a general strategy to identify spatial differences in enhancer activity in vivo.

Sexual dimorphism in an adaptive radiation: Does intersexual niche differentiation result in ecological character displacement?

Evolutionary radiations are one plausible explanation for the rich biodiversity on Earth. Adaptive radiations are the most studied form of evolutionary radiations, and ecological opportunity has been identified as one factor permitting them. Competition among individuals is supposedly highest in populations of conspecifics. Divergent modes of resource use might minimize trophic overlap, and thus intersexual competition, resulting in ecological character displacement between sexes. However, the role of intersexual differentiation in speciation processes is insufficiently studied. The few studies available suggest that intersexual niche differentiation exists in adaptive radiations, but their role within the radiation, and the extent of differentiation within the organism itself, remains largely unexplored. Here, we test the hypothesis that multiple morphological structures are affected by intersexual niche differentiation in "roundfin" Telmatherina, the first case where intersexual niche differentiation was demonstrated in an adaptive fish radiation. We show that sexes of two of the three morphospecies differ in several structural components of the head, all of these are likely adaptive. Sexual dimorphism is linked to the respective morphospecies-specific ecology and affects several axes of variation. Trait variation translates into different feeding modes, processing types, and habitat usages that add to interspecific variation in all three morphospecies. Intrasexual selection, that is, male-male competition, may contribute to variation in some of the traits, but appears unlikely in internal structures, which are invisible to other individuals. We conclude that intersexual variation adds to the adaptive diversity of roundfins and might play a key role in minimizing intersexual competition in emerging radiations.

Pleiotropy data resource as a primer for investigating co-morbidities/multi-morbidities and their role in disease

Most current biomedical and protein research focuses only on a small proportion of genes, which results in a lost opportunity to identify new gene-disease associations and explore new opportunities for therapeutic intervention. The International Mouse Phenotyping Consortium (IMPC) focuses on elucidating gene function at scale for poorly characterized and/or under-studied genes. A key component of the IMPC initiative is the implementation of a broad phenotyping pipeline, which is facilitating the discovery of pleiotropy. Characterizing pleiotropy is essential to identify gene-disease associations, and it is of particular importance when elucidating the genetic causes of syndromic disorders. Here we show how the IMPC is effectively uncovering pleiotropy and how the new mouse models and gene function hypotheses generated by the IMPC are increasing our understanding of the mammalian genome, forming the basis of new research and identifying new gene-disease associations.

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.

Pervasive admixture and the spread of a large-lipped form in a cichlid fish radiation

Adaptive radiations have proven important for understanding the mechanisms and processes underlying biological diversity. The convergence of form and function, as well as admixture and adaptive introgression, are common in adaptive radiations. However, distinguishing between these two scenarios remains a challenge for evolutionary research. The Midas cichlid species complex (Amphilophus spp.) is a prime example of adaptive radiation, with phenotypic diversification occurring at various stages of genetic differentiation. One species, A. labiatus, has large fleshy lips, is associated with rocky lake substrates, and occurs patchily within Lakes Nicaragua and Managua. By contrast, the similar, but thin-lipped, congener, A. citrinellus, is more common and widespread. We investigated the evolutionary history of the large-lipped form, specifically regarding whether the trait has evolved independently in both lakes from ancestral thin-lipped populations, or via dispersal and/or admixture events. We collected samples from distinct locations in both lakes, and assessed differences in morphology and ecology. Using RAD-seq, we genotyped thousands of SNPs to measure population structure and divergence, demographic history, and admixture. We found significant between-species differences in ecology and morphology, local intraspecific differences in body shape and trophic traits, but only limited intraspecific variation in lip shape. Despite clear ecological differences, our genomic approach uncovered pervasive admixture between the species and low genomic differentiation, with species within lakes being genetically more similar than species between lakes. Taken together, our results suggest a single origin of large-lips, followed by pervasive admixture and adaptive introgression, with morphology being driven by local ecological opportunities, despite ongoing gene-flow.

The genetic basis of coordinated plasticity across functional units in a Lake Malawi cichlid mapping population

Adaptive radiations are often stereotypical, as populations repeatedly specialize along conserved environmental axes. Phenotypic plasticity may be similarly stereotypical, as individuals respond to environmental cues. These parallel patterns of variation, which are often consistent across traits, have led researchers to propose that plasticity can facilitate predictable patterns of evolution along environmental gradients. This "flexible stem" model of evolution raises questions about the genetic nature of plasticity, including how complex is the genetic basis for plasticity? Is plasticity across traits mediated by many distinct loci, or few "global" regulators? To address these questions, we reared a hybrid cichlid mapping population on alternate diet regimes mimicking an important environmental axis. We show that plasticity across an array of ecologically relevant traits is generally morphologically integrated, such that traits respond in a coordinated manner, especially those with overlapping function. Our genetic data are more ambiguous. While our mapping experiment provides little evidence for global genetic regulators of plasticity, these data do contain a genetic signal for the integration of plasticity across traits. Overall, our data suggest a compromise between genetic modularity, whereby plasticity may evolve independently across traits, and low level but widespread genetic integration, establishing the potential for plasticity to experience coordinated evolution.

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.

Adaptation via pleiotropy and linkage: Association mapping reveals a complex genetic architecture within the stickleback Eda locus

Genomic mapping of the loci associated with phenotypic evolution has revealed genomic "hotspots," or regions of the genome that control multiple phenotypic traits. This clustering of loci has important implications for the speed and maintenance of adaptation and could be due to pleiotropic effects of a single mutation or tight genetic linkage of multiple causative mutations affecting different traits. The threespine stickleback (Gasterosteus aculeatus) is a powerful model for the study of adaptive evolution because the marine ecotype has repeatedly adapted to freshwater environments across the northern hemisphere in the last 12,000 years. Freshwater ecotypes have repeatedly fixed a 16 kilobase haplotype on chromosome IV that contains Ectodysplasin (Eda), a gene known to affect multiple traits, including defensive armor plates, lateral line sensory hair cells, and schooling behavior. Many additional traits have previously been mapped to a larger region of chromosome IV that encompasses the Eda freshwater haplotype. To identify which of these traits specifically map to this adaptive haplotype, we made crosses of rare marine fish heterozygous for the freshwater haplotype in an otherwise marine genetic background. Further, we performed fine-scale association mapping in a fully interbreeding, polymorphic population of freshwater stickleback to disentangle the effects of pleiotropy and linkage on the phenotypes affected by this haplotype. Although we find evidence that linked mutations have small effects on a few phenotypes, a small 1.4-kb region within the first intron of Eda has large effects on three phenotypic traits: lateral plate count, and both the number and patterning of the posterior lateral line neuromasts. Thus, the Eda haplotype is a hotspot of adaptation in stickleback due to both a small, pleiotropic region affecting multiple traits as well as multiple linked mutations affecting additional traits.

Genetic, morphological and ecological variation across a sharp hybrid zone between two alpine butterfly species

Identifying the mechanisms involved in the formation and maintenance of species is a central question in evolutionary biology, and distinguishing the selective drivers of populations' divergence from demographic processes is of particular interest to better understand the speciation process. Hybrid zones are recognized to provide ideal places to investigate the genetic architecture of speciation and to identify the mechanisms allowing diverging species to maintain their integrity in the face of gene flow. Here, we studied two alpine butterfly species, Coenonympha macromma and C. gardetta, which can be found flying together and hybridizing in narrow contact zones in the southern French Alps. We characterized the genomic composition of individuals, their morphology and their local habitat requirements, within and around a hybrid zone. Genetic diversity analysis at 794 SNPs revealed that all individuals within the hybrid zone were highly admixed, which was not the case outside the hybrid zone. Cline analysis showed that, despite ongoing hybridization, 56 out of 122 loci differentially fixed or nearly so between the two species were impermeable to introgression across the sharp hybrid zone (9 km wide). We also found concordance in cline position and width among genetic, morphological and environmental variation, suggesting a coupling of different reproductive barriers. Habitat characteristics such as the presence of trees and shrubs and the start of the growing season were strongly associated with the genetic variation, and we found evidence of divergence at genetic markers associated with morphology and physiology, putatively involved in visual or environmental reproductive isolation. We discuss the various behavioural and ecological factors that might interplay to maintain current levels of divergence and gene flow between this species pair.

Why does allometry evolve so slowly?

Morphological allometry is striking due to its evolutionary conservatism, making it an example of a certain sort of evolutionary stasis. Organisms that vary in size, whether for developmental, environmental, or evolutionary reasons, adopt shapes that are predictable from that size alone. There are two major hypotheses to explain this. It may be that natural selection strongly favors each allometric pattern, or that organisms lack the development and genetic capacity to produce variant shapes for selection to act on. Using a high-throughput system for measuring the size and shape of Drosophila wings, we documented an allometric pattern that has been virtually unchanged for 40 million years. We performed an artificial selection experiment on the static allometric slope within one species. In just 26 generations, we were able to increase the slope from 1.1 to 1.4, and decrease it to 0.8. Once artificial selection was suspended, the slope rapidly evolved back to a value near the initial static slope. This result decisively rules out the hypothesis that allometry is preserved due to a lack of genetic variation, and provides evidence that natural selection acts to maintain allometric relationships. On the other hand, it seems implausible that selection on allometry in the wing alone could be sufficiently strong to maintain static allometries over millions of years. This suggests that a potential explanation for stasis is selection on a potentially large number of pleiotropic effects. This seems likely in the case of allometry, as the sizes of all parts of the body may be altered when the allometric slope of one body part is changed. Unfortunately, hypotheses about pleiotropy have been very difficult to test. We lay out an approach to begin the systematic study of pleiotropic effects using genetic manipulations and high-throughput phenotyping.

Phenotypic divergence among threespine stickleback that differ in nuptial coloration

By studying systems in their earliest stages of differentiation, we can learn about the evolutionary forces acting within and among populations and how those forces could contribute to reproductive isolation. Such an understanding would help us to better discern and predict how selection leads to the maintenance of multiple morphs within a species, rather than speciation. The postglacial adaptive radiation of the threespine stickleback (Gasterosteus aculeatus) is one of the best-studied cases of evolutionary diversification and rapid, repeated speciation. Following deglaciation, marine stickleback have continually invaded freshwater habitats across the northern hemisphere and established resident populations that diverged innumerable times from their oceanic ancestors. Independent freshwater colonization events have yielded broadly parallel patterns of morphological differences in freshwater and marine stickleback. However, there is also much phenotypic diversity within and among freshwater populations. We studied a lesser-known freshwater "species pair" found in southwest Washington, where male stickleback in numerous locations have lost the ancestral red sexual signal and instead develop black nuptial coloration. We measured phenotypic variation in a suite of traits across sites where red and black stickleback do not overlap in distribution and at one site where they historically co-occurred. We found substantial phenotypic divergence between red and black morphs in noncolor traits including shape and lateral plating, and additionally find evidence that supports the hypothesis of sensory drive as the mechanism responsible for the evolutionary switch in color from red to black. A newly described third "mixed" morph in Connor Creek, Washington, differs in head shape and size from the red and black morphs, and we suggest that their characteristics are most consistent with hybridization between anadromous and freshwater stickleback. These results lay the foundation for future investigation of the underlying genetic basis of this phenotypic divergence as well as the evolutionary processes that may drive, maintain, or limit divergence among morphs.

Genome-wide association analyses reveal polygenic genomic architecture underlying divergent shell morphology in Spanish Littorina saxatilis ecotypes

Gene flow between diverging populations experiencing dissimilar ecological conditions can theoretically constrain adaptive evolution. To minimize the effect of gene flow, alleles underlying traits essential for local adaptation are predicted to be located in linked genome regions with reduced recombination. Local reduction in gene flow caused by selection is expected to produce elevated divergence in these regions. The highly divergent crab-adapted and wave-adapted ecotypes of the marine snail Littorina saxatilis present a model system to test these predictions. We used genome-wide association (GWA) analysis of geometric morphometric shell traits associated with microgeographic divergence between the two L. saxatilis ecotypes within three separate sampling sites. A total of 477 snails that had individual geometric morphometric data and individual genotypes at 4,066 single nucleotide polymorphisms (SNPs) were analyzed using GWA methods that corrected for population structure among the three sites. This approach allowed dissection of the genomic architecture of shell shape divergence between ecotypes across a wide geographic range, spanning two glacial lineages. GWA revealed 216 quantitative trait loci (QTL) with shell size or shape differences between ecotypes, with most loci explaining a small proportion of phenotypic variation. We found that QTL were evenly distributed across 17 linkage groups, and exhibited elevated interchromosomal linkage, suggesting a genome-wide response to divergent selection on shell shape between the two ecotypes. Shell shape trait-associated loci showed partial overlap with previously identified outlier loci under divergent selection between the two ecotypes, supporting the hypothesis of diversifying selection on these genomic regions. These results suggest that divergence in shell shape between the crab-adapted and wave-adapted ecotypes is produced predominantly by a polygenic genomic architecture with positive linkage disequilibrium among loci of small effect.

Genetics of adaptation: Experimental test of a biotic mechanism driving divergence in traits and genes

The genes underlying adaptations are becoming known, yet the causes of selection on genes-a key step in the study of the genetics of adaptation-remains uncertain. We address this issue experimentally in a threespine stickleback species pair showing exaggerated divergence in bony defensive armor in association with competition-driven character displacement. We used semi-natural ponds to test the role of a native predator in causing divergent evolution of armor and two known underlying genes. Predator presence/absence altered selection on dorsal spines and allele frequencies at the Msx2a gene across a generation. Evolutionary trajectories of alleles at a second gene, Pitx1, and the pelvic spine trait it controls, were more variable. Our experiment demonstrates how manipulation of putative selective agents helps to identify causes of evolutionary divergence at key genes, rule out phenotypic plasticity as a sole determinant of phenotypic differences, and eliminate reliance on fitness surrogates. Divergence of predation regimes in sympatric stickleback is associated with coevolution in response to resource competition, implying a cascade of biotic interactions driving species divergence. We suggest that as divergence proceeds, an increasing number of biotic interactions generate divergent selection, causing more evolution in turn. In this way, biotic adaptation perpetuates species divergence through time during adaptive radiation in an expanding number of traits and genes.

Marker-trait associations and genomic predictions of interspecific pear (Pyrus) fruit characteristics

Interspecific pear (Pyrus spp.) hybrid populations are often used to develop novel cultivars. Pear cultivar breeding is a lengthy process because of long juvenility and the subsequent time required for reliable fruit phenotyping. Molecular techniques such as genome-wide association (GWA) and genomic selection (GS) provide an opportunity to fast-forward the development of high-value cultivars. We evaluated the genetic architecture of 10 pear fruit phenotypes (including sensory traits) and the potential of GS using genotyping-by-sequencing of 550 hybrid seedlings from nine interrelated full-sib families. Results from GWA suggested a complex polygenic nature of all 10 traits as the maximum variance explained by each marker was less than 4% of the phenotypic variance. The effect-size of SNPs for each trait suggested many genes of small effect and few of moderate effect. Some genomic regions associated with pear sensory traits were similar to those reported for apple - possibly a result of high synteny between the apple and pear genomes. The average (across nine families) GS accuracy varied from 0.32 (for crispness) to 0.62 (for sweetness), with an across-trait average of 0.42. Further efforts are needed to develop larger genotype-phenotype datasets in order to predict fruit phenotypes of untested seedlings with sufficient efficiency.

Eco-evolutionary feedbacks between prey densities and linkage disequilibrium in the predator maintain diversity

Diversity occurs at multiple scales. Within a single population, there is diversity in genotypes and phenotypes. At a larger scale, within ecological communities, there is diversity in species. A number of studies have investigated how diversity at these two scales influence each other through what has been termed eco-evolutionary feedbacks. Here we study a three-species ecological module called apparent competition, in which the predator is evolving in a trait that determines its interaction with two prey species. Unlike previous studies on apparent competition, which employed evolutionary frameworks with very simple genetics, we study an eco-evolutionary model in which the predator's trait is determined by two recombining diallelic loci, so that its mean and variance can evolve, as well as associations (linkage disequilibrium) between the loci. We ask how eco-evolutionary feedbacks with these two loci affect the coexistence of the prey species and the maintenance of polymorphisms within the predator species. We uncover a novel eco-evolutionary feedback between the prey densities and the linkage disequilibrium between the predator's loci. Through a stability analysis, we demonstrate how these feedbacks affect polymorphisms at both loci and, among others, may generate stable cycling.

Comparative transcriptome analysis reveals potential evolutionary differences in adaptation of temperature and body shape among four Percidae species

Considering the divergent temperature habitats and morphological traits of four Percidae species: yellow perch (Perca flavescens), Eurasian perch (Perca fluviatilis), pike perch (Sander lucioperca), and ruffe (Gymnocephalus cernua), we stepped into the transcriptome level to discover genes and mechanisms that drive adaptation to different temperature environments and evolution in body shape. Based on 93,566 to 181,246 annotated unigenes of the four species, we identified 1,117 one-to-one orthologous genes and subsequently constructed the phylogenetic trees that are consistent with previous studies. Together with the tree, the ratios of nonsynonymous to synonymous substitutions presented decreased evolutionary rates from the D. rerio branch to the sub-branch clustered by P. flavescens and P. fluviatilis. The specific 93 fast-evolving genes and 57 positively selected genes in P. flavescens, compared with 22 shared fast-evolving genes among P. fluviatilis, G. cernua, and S. lucioperca, showed an intrinsic foundation that ensure its adaptation to the warmer Great Lakes and farther south, especially in functional terms like “Cul4-RING E3 ubiquitin ligase complex.” Meanwhile, the specific 78 fast-evolving genes and 41 positively selected genes in S. lucioperca drew a clear picture of how it evolved to a large and elongated body with camera-type eyes and muscle strength so that it could occupy the highest position in the food web. Overall, our results uncover genetic basis that support evolutionary adaptation of temperature and body shape in four Percid species, and could furthermore assist studies on environmental adaptation in fishes.

Assessing the Architecture of Drosophila mojavensis Locomotor Evolution with Bulk Segregant Analysis

Behavior is frequently predicted to be especially important for evolution in novel environments. If these predictions are accurate, there might be particular patterns of genetic architecture associated with recently diverged behaviors. Specifically, it has been predicted that behaviors linked to population divergence should be underpinned by a few genes of relatively large effect, compared to architectures of intrapopulation behavioral variation, which is considered to be highly polygenic. More mapping studies of behavioral variation between recently diverged populations are needed to continue assessing the generality of these predictions. Here, we used a bulk segregant mapping approach to dissect the genetic architecture of a locomotor trait that has evolved between two populations of the cactophilic fly Drosophila mojavensis We created an F8 mapping population of 1,500 individuals from advanced intercross lines and sequenced the 10% of individuals with the highest and lowest levels of locomotor activity. Using three alternative statistical approaches, we found strong evidence for two relatively large-effect QTL that is localized in a region homologous to a region of densely packed behavior loci in Drosophila melanogaster, suggesting that clustering of behavior genes may display relatively deep evolutionary conservation. Broadly, our data are most consistent with a polygenic architecture, though with several loci explaining a high proportion of variation in comparison to similar behavioral traits. We further note the presence of several antagonistic QTL linked to locomotion and discuss these results in light of theories regarding behavioral evolution and the effect size and direction of QTL for diverging traits in general.

Pleiotropic Effects of ebony and tan on Pigmentation and Cuticular Hydrocarbon Composition in Drosophila melanogaster

Pleiotropic genes are genes that affect more than one trait. For example, many genes required for pigmentation in the fruit fly Drosophila melanogaster also affect traits such as circadian rhythms, vision, and mating behavior. Here, we present evidence that two pigmentation genes, ebony and tan, which encode enzymes catalyzing reciprocal reactions in the melanin biosynthesis pathway, also affect cuticular hydrocarbon (CHC) composition in D. melanogaster females. More specifically, we report that ebony loss-of-function mutants have a CHC profile that is biased toward long (>25C) chain CHCs, whereas tan loss-of-function mutants have a CHC profile that is biased toward short (<25C) chain CHCs. Moreover, pharmacological inhibition of dopamine synthesis, a key step in the melanin synthesis pathway, reversed the changes in CHC composition seen in ebony mutants, making the CHC profiles similar to those seen in tan mutants. These observations suggest that genetic variation affecting ebony and/or tan activity might cause correlated changes in pigmentation and CHC composition in natural populations. We tested this possibility using the Drosophila Genetic Reference Panel (DGRP) and found that CHC composition covaried with pigmentation as well as levels of ebony and tan expression in newly eclosed adults in a manner consistent with the ebony and tan mutant phenotypes. These data suggest that the pleiotropic effects of ebony and tan might contribute to covariation of pigmentation and CHC profiles in Drosophila.

Genetic dissection of assortative mating behavior

The evolution of new species is made easier when traits under divergent ecological selection are also mating cues. Such ecological mating cues are now considered more common than previously thought, but we still know little about the genetic changes underlying their evolution or more generally about the genetic basis for assortative mating behaviors. Both tight physical linkage and the existence of large-effect preference loci will strengthen genetic associations between behavioral and ecological barriers, promoting the evolution of assortative mating. The warning patterns of Heliconius melpomene and H. cydno are under disruptive selection due to increased predation of nonmimetic hybrids and are used during mate recognition. We carried out a genome-wide quantitative trait locus (QTL) analysis of preference behaviors between these species and showed that divergent male preference has a simple genetic basis. We identify three QTLs that together explain a large proportion (approximately 60%) of the difference in preference behavior observed between the parental species. One of these QTLs is just 1.2 (0-4.8) centiMorgans (cM) from the major color pattern gene optix, and, individually, all three have a large effect on the preference phenotype. Genomic divergence between H. cydno and H. melpomene is high but broadly heterogenous, and admixture is reduced at the preference-optix color pattern locus but not the other preference QTLs. The simple genetic architecture we reveal will facilitate the evolution and maintenance of new species despite ongoing gene flow by coupling behavioral and ecological aspects of reproductive isolation.

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.

Repeated Selection of Alternatively Adapted Haplotypes Creates Sweeping Genomic Remodeling in Stickleback

Heterogeneous genetic divergence can accumulate across the genome when populations adapt to different habitats while still exchanging alleles. How long does diversification take and how much of the genome is affected? When divergence occurs in parallel from standing genetic variation, how often are the same haplotypes involved? We explore these questions using restriction site-associated DNA sequencing genotyping data and show that broad-scale genomic repatterning, fueled by copious standing variation, can emerge in just dozens of generations in replicate natural populations of threespine stickleback fish (Gasterosteus aculeatus). After the catastrophic 1964 Alaskan earthquake, marine stickleback colonized newly created ponds on seismically uplifted islands. We find that freshwater fish in these young ponds differ from their marine ancestors across the same genomic segments previously shown to have diverged in much older lake populations. Outside of these core divergent regions the genome shows no population structure across the ocean-freshwater divide, consistent with strong local selection acting in alternative environments on stickleback populations still connected by significant gene flow. Reinforcing this inference, a majority of divergent haplotypes that are at high frequency in ponds are detectable in the sea, even across great geographic distances. Building upon previous population genomics work in this model species, our data suggest that a long history of divergent selection and gene flow among stickleback populations in oceanic and freshwater habitats has maintained polymorphisms of alternatively adapted DNA sequences that facilitate parallel evolution.

Mapping morphological shape as a high-dimensional functional curve

Detecting how genes regulate biological shape has become a multidisciplinary research interest because of its wide application in many disciplines. Despite its fundamental importance, the challenges of accurately extracting information from an image, statistically modeling the high-dimensional shape and meticulously locating shape quantitative trait loci (QTL) affect the progress of this research. In this article, we propose a novel integrated framework that incorporates shape analysis, statistical curve modeling and genetic mapping to detect significant QTLs regulating variation of biological shape traits. After quantifying morphological shape via a radius centroid contour approach, each shape, as a phenotype, was characterized as a high-dimensional curve, varying as angle θ runs clockwise with the first point starting from angle zero. We then modeled the dynamic trajectories of three mean curves and variation patterns as functions of θ. Our framework led to the detection of a few significant QTLs regulating the variation of leaf shape collected from a natural population of poplar, Populus szechuanica var tibetica. This population, distributed at altitudes 2000-4500 m above sea level, is an evolutionarily important plant species. This is the first work in the quantitative genetic shape mapping area that emphasizes a sense of 'function' instead of decomposing the shape into a few discrete principal components, as the majority of shape studies do.

Dorsal spine evolution in threespine sticklebacks via a splicing change in MSX2A

BACKGROUND

Dorsal spine reduction in threespine sticklebacks (Gasterosteus aculeatus) is a classic example of recurrent skeletal evolution in nature. Sticklebacks in marine environments typically have long spines that form part of their skeletal armor. Many derived freshwater populations have evolved shorter spines. Changes in spine length are controlled in part by a quantitative trait locus (QTL) previously mapped to chromosome 4, but the causative gene and mutations underlying the repeated evolution of this interesting skeletal trait have not been identified.

RESULTS

Refined mapping of the spine length QTL shows that it lies near the MSX2A transcription factor gene. MSX2A is expressed in developing spines. In F1 marine × freshwater fish, the marine allele is preferentially expressed. Differences in expression can be attributed to splicing regulation. Due to the use of an alternative 5 ^' splice site within the first exon, the freshwater allele produces greater amounts of a shortened, non-functional transcript and makes less of the full-length transcript. Sequence changes in the MSX2A region are shared by many freshwater fish, suggesting that repeated evolution occurs by reuse of a spine-reduction variant. To demonstrate the effect of full-length MSX2A on spine length, we produced transgenic freshwater fish expressing a copy of marine MSX2A. The spines of the transgenic fish were significantly longer on average than those of their non-transgenic siblings, partially reversing the reduced spine lengths that have evolved in freshwater populations.

CONCLUSIONS

MSX2A is a major gene underlying dorsal spine reduction in freshwater sticklebacks. The gene is linked to a separate gene controlling bony plate loss, helping explain the concerted effects of chromosome 4 on multiple armor-reduction traits. The nature of the molecular changes provides an interesting example of morphological evolution occurring not through a simple amino acid change, nor through a change only in gene expression levels, but through a change in the ratio of splice products encoding both normal and truncated proteins.

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

Developmental timing differences underlie armor loss across threespine stickleback populations

Comparing ontogenetic patterns within a well-described evolutionary context aids in inferring mechanisms of change, including heterochronies or deletion of developmental pathways. Because selection acts on phenotypes throughout ontogeny, any within-taxon developmental variation has implications for evolvability. We compare ontogenetic order and timing of locomotion and defensive traits in three populations of threespine stickleback that have evolutionarily divergent adult forms. This analysis adds to the growing understanding of developmental genetic mechanisms of adaptive change in this evolutionary model species by delineating when chondrogenesis and osteogenesis in two derived populations begin to deviate from the developmental pattern in their immediate ancestors. We found that differences in adult defensive morphologies arise through abolished or delayed initiation of these traits rather than via an overall heterochronic shift, that intra-population ontogenetic variation is increased for some derived traits, and that altered armor developmental timing differentiates the derived populations from each other despite parallels in adult lateral plate armor phenotypes. We found that changes in ossified elements of the pelvic armor are linked to delayed and incomplete development of an early-forming pelvic cartilage, and that this disruption likely presages the variable pelvic vestiges documented in many derived populations.

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.

Body shape differences in a pair of closely related Malawi cichlids and their hybrids: Effects of genetic variation, phenotypic plasticity, and transgressive segregation

Phenotypic differences may have genetic and plastic components. Here, we investigated the contributions of both for differences in body shape in two species of Lake Malawi cichlids using wild-caught specimens and a common garden experiment. We further hybridized the two species to investigate the mode of gene action influencing body shape differences and to examine the potential for transgressive segregation. We found that body shape differences between the two species observed in the field are maintained after more than 10 generations in a standardized environment. Nonetheless, both species experienced similar changes in the laboratory environment. Our hybrid cross experiment confirmed that substantial variation in body shape appears to be genetically determined. The data further suggest that the underlying mode of gene action is complex and cannot be explained by simple additive or additive-dominance models. Transgressive phenotypes were found in the hybrid generations, as hybrids occupied significantly more morphospace than both parentals combined. Further, the body shapes of transgressive individuals resemble the body shapes observed in other Lake Malawi rock-dwelling genera. Our findings indicate that body shape can respond to selection immediately, through plasticity, and over longer timescales through adaptation. In addition, our results suggest that hybridization may have played an important role in the diversification of Lake Malawi cichlids through creating new phenotypic variation.

What can genome-wide association studies tell us about the evolutionary forces maintaining genetic variation for quantitative traits?

Contents 21 I. 21 II. 22 III. 24 IV. 25 V. 29 30 References 30 SUMMARY: Understanding the evolutionary forces that shape genetic variation within species has long been a goal of evolutionary biology. Integrating data for the genetic architecture of traits from genome-wide association mapping studies (GWAS) along with the development of new population genetic methods for identifying selection in sequence data may allow us to evaluate the roles of mutation-selection balance and balancing selection in shaping genetic variation at various scales. Here, we review the theoretical predictions for genetic architecture and additional signals of selection on genomic sequence for the loci that affect traits. Next, we review how plant GWAS have tested for the signatures of various selective scenarios. Limited evidence to date suggests that within-population variation is maintained primarily by mutation-selection balance while variation across the landscape is the result of local adaptation. However, there are a number of inherent biases in these interpretations. We highlight these challenges and suggest ways forward to further understanding of the maintenance of variation.

Quantitative gene-gene and gene-environment mapping for leaf shape variation using tree-based models

Leaf shape traits have long been a focus of many disciplines, but the complex genetic and environmental interactive mechanisms regulating leaf shape variation have not yet been investigated in detail. The question of the respective roles of genes and environment and how they interact to modulate leaf shape is a thorny evolutionary problem, and sophisticated methodology is needed to address it. In this study, we investigated a framework-level approach that inputs shape image photographs and genetic and environmental data, and then outputs the relative importance ranks of all variables after integrating shape feature extraction, dimension reduction, and tree-based statistical models. The power of the proposed framework was confirmed by simulation and a Populus szechuanica var. tibetica data set. This new methodology resulted in the detection of novel shape characteristics, and also confirmed some previous findings. The quantitative modeling of a combination of polygenetic, plastic, epistatic, and gene-environment interactive effects, as investigated in this study, will improve the discernment of quantitative leaf shape characteristics, and the methods are ready to be applied to other leaf morphology data sets. Unlike the majority of approaches in the quantitative leaf shape literature, this framework-level approach is data-driven, without assuming any pre-known shape attributes, landmarks, or model structures.

Non-reef environments impact the diversification of extant jacks, remoras and allies (Carangoidei, Percomorpha)

Various factors may impact the processes of diversification of a clade. In the marine realm, it has been shown that coral reef environments have promoted diversification in various fish groups. With the exception of requiem sharks, all the groups showing a higher level of diversity in reefs than in non-reef habitats have diets based predominantly on plankton, algae or benthic invertebrates. Here we explore the pattern of diversification of carangoid fishes, a clade that includes numerous piscivorous species (e.g. trevallies, jacks and dolphinfishes), using time-calibrated phylogenies as well as ecological and morphological data from both extant and fossil species. The study of carangoid morphospace suggests that reef environments played a role in their early radiation during the Eocene. However, contrary to the hypothesis of a reef-association-promoting effect, we show that habitat shifts to non-reef environments have increased the rates of morphological diversification (i.e. size and body shape) in extant carangoids. Piscivory did not have a major impact on the tempo of diversification of this group. Through the ecological radiation of carangoid fishes, we demonstrate that non-reef environments may sustain and promote processes of diversification of different marine fish groups, at least those including a large proportion of piscivorous species.

The genetic architecture of local adaptation and reproductive isolation in sympatry within the Mimulus guttatus species complex

The genetic architecture of local adaptation has been of central interest to evolutionary biologists since the modern synthesis. In addition to classic theory on the effect size of adaptive mutations by Fisher, Kimura and Orr, recent theory addresses the genetic architecture of local adaptation in the face of ongoing gene flow. This theory predicts that with substantial gene flow between populations local adaptation should proceed primarily through mutations of large effect or tightly linked clusters of smaller effect loci. In this study, we investigate the genetic architecture of divergence in flowering time, mating system-related traits, and leaf shape between Mimulus laciniatus and a sympatric population of its close relative M. guttatus. These three traits are probably involved in M. laciniatus' adaptation to a dry, exposed granite outcrop environment. Flowering time and mating system differences are also reproductive isolating barriers making them 'magic traits'. Phenotypic hybrids in this population provide evidence of recent gene flow. Using next-generation sequencing, we generate dense SNP markers across the genome and map quantitative trait loci (QTLs) involved in flowering time, flower size and leaf shape. We find that interspecific divergence in all three traits is due to few QTL of large effect including a highly pleiotropic QTL on chromosome 8. This QTL region contains the pleiotropic candidate gene TCP4 and is involved in ecologically important phenotypes in other Mimulus species. Our results are consistent with theory, indicating that local adaptation and reproductive isolation with gene flow should be due to few loci with large and pleiotropic effects.

Adaptation and acclimation of traits associated with swimming capacity in Lake Whitefish (coregonus clupeaformis) ecotypes

Background

Improved performance in a given ecological niche can occur through local adaptation, phenotypic plasticity, or a combination of these mechanisms. Evaluating the relative importance of these two mechanisms is needed to better understand the cause of intra specific polymorphism. In this study, we reared populations of Lake Whitefish (Coregonus clupeaformis) representing the’normal’ (benthic form) and the ‘dwarf’ (derived limnetic form) ecotypes in two different conditions (control and swim-training) to test the relative importance of adaptation and acclimation in the differentiation of traits related to swimming capacity. The dwarf whitefish is a more active swimmer than the normal ecotype, and also has a higher capacity for aerobic energy production in the swimming musculature. We hypothesized that dwarf fish would show changes in morphological and physiological traits consistent with reductions in the energetic costs of swimming and maintenance metabolism.

Results

We found differences in traits predicted to decrease the costs of prolonged swimming and standard metabolic rate and allow for a more active lifestyle in dwarf whitefish. Dwarf whitefish evolved a more streamlined body shape, predicted to lead to a decreased drag, and a smaller brain, which may decrease their standard metabolic rate. Contrary to predictions, we also found evidence of acclimation in liver size and metabolic enzyme activities.

Conclusion

Results support the view that local adaptation has contributed to the genetically-based divergence of traits associated with swimming activity. Presence of post-zygotic barriers limiting gene flow between these ecotype pairs may have favoured repeated local adaptation to the limnetic niches.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-016-0732-y) contains supplementary material, which is available to authorized users.

Ectodysplasin signalling genes and phenotypic evolution in sculpins (Cottus)

Despite their deeply conserved function among vertebrates, ectodysplasin (Eda) signalling genes are involved in microevolutionary change in humans and sticklebacks. If such a dual role is common, Eda signalling genes constitute hotspots for morphological evolution. Variation in sculpin (Cottus) skin prickling and body shape resembles patterns caused by variation in Eda signalling in sticklebacks. We mapped Eda signalling genes and performed quantitative trait locus mapping in crosses between Cottus rhenanus and Cottus perifretum. A genomic region containing the Eda receptor (Edar) was strongly associated with prickling and contributed to shape. The expression of Edar in developing prickles and skeletal elements in Cottus was confirmed by in situ hybridization. Coding sequence changes between Edar alleles in C. rhenanus and C. perifretum exceeded sequence differentiation in other vertebrates. However, it is likely that additional genetic elements besides coding changes affect the phenotypic variation. Although the phenotype in a natural hybrid lineage between C. rhenanus and C. perifretum resembles C. perifretum, the respective coding Edar alleles are not fully fixed (88.6%). Hence, our results support an involvement of Eda signalling in microevolutionary changes, but imply that the Edar gene is affected by multiple evolutionary processes that vary among freshwater sculpins.

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.

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.

Innate immune responses to gut microbiota differ between oceanic and freshwater threespine stickleback populations

Animal hosts must co-exist with beneficial microbes while simultaneously being able to mount rapid, non-specific, innate immune responses to pathogenic microbes. How this balance is achieved is not fully understood, and disruption of this relationship can lead to disease. Excessive inflammatory responses to resident microbes are characteristic of certain gastrointestinal pathologies such as inflammatory bowel disease (IBD). The immune dysregulation of IBD has complex genetic underpinnings that cannot be fully recapitulated with single-gene-knockout models. A deeper understanding of the genetic regulation of innate immune responses to resident microbes requires the ability to measure immune responses in the presence and absence of the microbiota using vertebrate models with complex genetic variation. Here, we describe a new gnotobiotic vertebrate model to explore the natural genetic variation that contributes to differences in innate immune responses to microbiota. Threespine stickleback, Gasterosteus aculeatus, has been used to study the developmental genetics of complex traits during the repeated evolution from ancestral oceanic to derived freshwater forms. We established methods to rear germ-free stickleback larvae and gnotobiotic animals monoassociated with single bacterial isolates. We characterized the innate immune response of these fish to resident gut microbes by quantifying the neutrophil cells in conventionally reared monoassociated or germ-free stickleback from both oceanic and freshwater populations grown in a common intermediate salinity environment. We found that oceanic and freshwater fish in the wild and in the laboratory share many intestinal microbial community members. However, oceanic fish mount a strong immune response to residential microbiota, whereas freshwater fish frequently do not. A strong innate immune response was uniformly observed across oceanic families, but this response varied among families of freshwater fish. The gnotobiotic stickleback model that we have developed therefore provides a platform for future studies mapping the natural genetic basis of the variation in immune response to microbes.

Summary: We developed a gnotobiotic threespine stickleback system to identify the cellular and genetic bases for variation in inflammatory responses to microbiota.

How predation shaped fish: the impact of fin spines on body form evolution across teleosts

It is well known that predators can induce morphological changes in some fish: individuals exposed to predation cues increase body depth and the length of spines. We hypothesize that these structures may evolve synergistically, as together, these traits will further enlarge the body dimensions of the fish that gape-limited predators must overcome. We therefore expect that the orientation of the spines will predict which body dimension increases in the presence of predators. Using phylogenetic comparative methods, we tested this prediction on the macroevolutionary scale across 347 teleost families, which display considerable variation in fin spines, body depth and width. Consistent with our predictions, we demonstrate that fin spines on the vertical plane (dorsal and anal fins) are associated with a deeper-bodied optimum. Lineages with spines on the horizontal plane (pectoral fins) are associated with a wider-bodied optimum. Optimal body dimensions across lineages without spines paralleling the body dimension match the allometric expectation. Additionally, lineages with longer spines have deeper and wider body dimensions. This evolutionary relationship between fin spines and body dimensions across teleosts reveals functional synergy between these two traits and a potential macroevolutionary signature of predation on the evolutionary dynamics of body shape.

The Genetic Architecture of Domestication in Animals

Domestication has been essential to the progress of human civilization, and the process itself has fascinated biologists for hundreds of years. Domestication has led to a series of remarkable changes in a variety of plants and animals, in what is termed the “domestication phenotype.” In domesticated animals, this general phenotype typically consists of similar changes in tameness, behavior, size/morphology, color, brain composition, and adrenal gland size. This domestication phenotype is seen in a range of different animals. However, the genetic basis of these associated changes is still puzzling. The genes for these different traits tend to be grouped together in clusters in the genome, though it is still not clear whether these clusters represent pleiotropic effects, or are in fact linked clusters. This review focuses on what is currently known about the genetic architecture of domesticated animal species, if genes of large effect (often referred to as major genes) are prevalent in driving the domestication phenotype, and whether pleiotropy can explain the loci underpinning these diverse traits being colocated.

Iterative development and the scope for plasticity: contrasts among trait categories in an adaptive radiation

Phenotypic plasticity can influence evolutionary change in a lineage, ranging from facilitation of population persistence in a novel environment to directing the patterns of evolutionary change. As the specific nature of plasticity can impact evolutionary consequences, it is essential to consider how plasticity is manifested if we are to understand the contribution of plasticity to phenotypic evolution. Most morphological traits are developmentally plastic, irreversible, and generally considered to be costly, at least when the resultant phenotype is mis-matched to the environment. At the other extreme, behavioral phenotypes are typically activational (modifiable on very short time scales), and not immediately costly as they are produced by constitutive neural networks. Although patterns of morphological and behavioral plasticity are often compared, patterns of plasticity of life history phenotypes are rarely considered. Here we review patterns of plasticity in these trait categories within and among populations, comprising the adaptive radiation of the threespine stickleback fish Gasterosteus aculeatus. We immediately found it necessary to consider the possibility of iterated development, the concept that behavioral and life history trajectories can be repeatedly reset on activational (usually behavior) or developmental (usually life history) time frames, offering fine tuning of the response to environmental context. Morphology in stickleback is primarily reset only in that developmental trajectories can be altered as environments change over the course of development. As anticipated, the boundaries between the trait categories are not clear and are likely to be linked by shared, underlying physiological and genetic systems.

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

How predictable is the genetic basis of phenotypic adaptation? Answering this question begins by estimating the repeatability of adaptation at the genetic level. Here, we provide a comprehensive estimate of the repeatability of the genetic basis of adaptive phenotypic evolution in a natural system. We used quantitative trait locus (QTL) mapping to discover genomic regions controlling a large number of morphological traits that have diverged in parallel between pairs of threespine stickleback (Gasterosteus aculeatus species complex) in Paxton and Priest lakes, British Columbia. We found that nearly half of QTL affected the same traits in the same direction in both species pairs. Another 40% influenced a parallel phenotypic trait in one lake but not the other. The remaining 10% of QTL had phenotypic effects in opposite directions in the two species pairs. Similarity in the proportional contributions of all QTL to parallel trait differences was about 0.4. Surprisingly, QTL reuse was unrelated to phenotypic effect size. Our results indicate that repeated use of the same genomic regions is a pervasive feature of parallel phenotypic adaptation, at least in sticklebacks. Identifying the causes of this pattern would aid prediction of the genetic basis of phenotypic evolution.

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.

Where the lake meets the sea: strong reproductive isolation is associated with adaptive divergence between lake resident and anadromous three-spined sticklebacks

Contact zones between divergent forms of the same species are often characterised by high levels of phenotypic diversity over small geographic distances. What processes are involved in generating such high phenotypic diversity? One possibility is that introgression and recombination between divergent forms in contact zones results in greater phenotypic and genetic polymorphism. Alternatively, strong reproductive isolation between forms may maintain distinct phenotypes, preventing homogenisation by gene flow. Contact zones between divergent freshwater-resident and anadromous stickleback (Gasterosteus aculeatus L.) forms are numerous and common throughout the species distribution, offering an opportunity to examine these contrasting hypotheses in greater detail. This study reports on an interesting new contact zone located in a tidally influenced lake catchment in western Ireland, characterised by high polymorphism for lateral plate phenotypes. Using neutral and QTL-linked microsatellite markers, we tested whether the high diversity observed in this contact zone arose as a result of introgression or reproductive isolation between divergent forms: we found strong support for the latter hypothesis. Three phenotypic and genetic clusters were identified, consistent with two divergent resident forms and a distinct anadromous completely plated population that migrates in and out of the system. Given the strong neutral differentiation detected between all three morphotypes (mean F_ST = 0.12), we hypothesised that divergent selection between forms maintains reproductive isolation. We found a correlation between neutral genetic and adaptive genetic differentiation that support this. While strong associations between QTL linked markers and phenotypes were also observed in this wild population, our results support the suggestion that such associations may be more complex in some Atlantic populations compared to those in the Pacific. These findings provide an important foundation for future work investigating the dynamics of gene flow and adaptive divergence in this newly discovered stickleback contact zone.