Constraints on utilization of the EDA-signaling pathway in threespine stickleback evolution

Many traits evolve in parallel in widely separated populations. The evolutionary radiation of threespine sticklebacks provides a powerful model for testing the molecular basis of parallel evolution in vertebrates. Although marine sticklebacks are completely covered with bony armor plates, most freshwater populations have dramatic reductions in plates. Recent genetic studies have shown that major changes in armor patterning are likely due to regulatory alterations in the gene encoding the secreted signaling molecule ectodysplasin (EDA). In mammals, mutations in many different components of the EDA-signaling pathway produce similar changes in hair, teeth, sweat glands, and dermal bones. To test whether other genes in the EDA pathway also control natural variation in armor plates, we identified and mapped stickleback EDA Receptor (EDAR), the EDAR-Associated Death Domain adaptor, Tumor Necrosis Factor Receptor (TNFR) SuperFamily member 19, its adaptor TNFR-Associated Factor 6, and the downstream regulator nuclear factor kappa B Essential Modulator (NEMO). In contrast to the diversity of genes underlying ectodermal dysplasia disease phenotypes in humans, none of these EDA pathway components map to chromosomes previously shown to modify armor plates in natural populations, though EDAR showed a small but significant effect on plate number. We further investigated whether these genes exhibit differences in copy number, target size, or genomic organization that might make them less suitable targets for evolutionary change. In comparison with EDA, all these genes have smaller surrounding noncoding (putative regulatory) regions, with fewer evolutionarily conserved regions. We suggest that the presence of highly modular cis-acting control sequences may be a key factor influencing the likelihood that particular genes will serve as the basis of major phenotypic changes in nature.     

Adaptive evolution of lateral plates in three‐spined stickleback Gasterosteus aculeatus: a case study in functional analysis of natural variation

The three‐spined stickleback Gasterosteus aculeatus is a model species for studying questions in ecology and evolution. The rapid diversification of G. aculeatus in post‐glacial freshwater environments, combined with recently developed molecular tools, provides a unique opportunity to study the functional basis of fitness variation in natural populations. In derived freshwater populations, a number of morphological traits have diverged in parallel from the marine ancestral state, including the number of lateral armour plates. Evolution of reduced armour in freshwater populations is due to positive selection from both abiotic and biotic mechanisms. The major effect gene (ectodysplasin‐A or Eda), along with several minor effect genetic regions, has recently been shown to control lateral plate variation. Field experiments have further determined the fitness consequences of allelic variation at the major effect locus. This work helps elucidate the mechanisms connecting genetic variation with phenotypic variation and fitness in the wild, a synthesis that should be applicable to many other phenotypic traits and species of fishes.     

Parallel evolution of Pitx1 underlies pelvic reduction in Scottish threespine stickleback (Gasterosteus aculeatus)

Little is known about the genetic and molecular mechanisms that underlie adaptive phenotypic variation in natural populations or whether similar genetic and molecular mechanisms are utilized when similar adaptive phenotypes arise in independent populations. The threespine stickleback (Gasterosteus aculeatus) is a good model system to investigate these questions because these fish display a large amount of adaptive phenotypic variation, and similar adaptive phenotypes have arisen in multiple, independent stickleback populations. A particularly striking pattern of parallel evolution in sticklebacks is reduction of skeletal armor, which has occurred in numerous freshwater locations around the world. New genetic and genomic tools for the threespine stickleback have made it possible to identify genes that underlie loss of different elements of the skeletal armor. Previous work has shown that regulatory mutations at the Pitx1 locus are likely responsible for loss of the pelvic structures in independent stickleback populations from North America and Iceland. Here we show that the Pitx1 locus is also likely to underlie pelvic reduction in a Scottish population of threespine stickleback, which has apparently evolved pelvic reduction under a different selection regime than the North American populations.     

Heterozygosity and asymmetry: Ectodysplasin as a form of genetic stress in marine threespine stickleback

Genome‐wide heterozygosity has long been hypothesized to play a role in buffering organisms against developmental perturbations, potentially resulting in increased symmetry. If true, this could in part explain the maintenance of standing genetic variation in wild populations. Marine threespine sticklebacks (Gasterosteus aculeatus) were sampled across their eastern Pacific coastal distribution from Alaska to California and variations in asymmetry for both structural and nonstructural armor traits (lateral plates) were assessed. Structural plates consistently showed less asymmetry than nonstructural plates, but standardized measures of heterozygosity were not correlated with the extent of asymmetry expressed by a fish. Fish that were heterozygous for the major‐effect gene controlling lateral plate variation (Ectodysplasin) had higher occurrences of asymmetry, even when the individuals were phenotypically fully plated. Collectively, this suggests that heterozygosity at a major‐effect locus can have a greater impact on asymmetry than heterozygosity sampled across the genome.     

Parallel evolution and inheritance of quantitative traits

Parallel phenotypic evolution, the independent evolution of the same trait in closely related lineages, is interesting because it tells us about the contribution of natural selection to phenotypic evolution. Haldane and others have proposed that parallel evolution also results from a second process, the similarly biased production of genetic variation in close relatives, an idea that has received few tests. We suggest that influence of shared genetic biases should be detectable by the disproportionate use of the same genes in independent instances of parallel phenotypic evolution. We show how progress in testing this prediction can be made through simple tests of parallel inheritance of genetic differences: similar additive, dominance, and epistasis components in analysis of line means and similar effective numbers of loci. We demonstrate parallel inheritance in two traits, lateral plate number and body shape, in two lineages of threespine stickleback that have adapted independently to freshwater streams on opposite sides of the Pacific Ocean. Notably, reduction of plate number in freshwater involves a substitution at the same major locus in both lineages. Our results represent only a first step in the study of the genetics of parallel phenotypic evolution in sticklebacks. Nevertheless, we have shown how such studies can be employed to test the genetic hypothesis of parallel evolution and how study of parallel evolution might yield insights into the roles of both selection and genetic constraint in phenotypic evolution.     

Parallel evolution by correlated response: lateral plate reduction in threespine stickleback

Recent work has revealed the molecular mechanisms governing one of the most dramatic examples of parallel evolution in nature: the repeated loss of lateral plate armor in freshwater populations of threespine stickleback. Yet, the ecological mechanisms responsible for armor loss remain unclear. Using a balanced experimental design, we examined Heuts' (1947) hypothesis that selection due to differences in salinity indirectly drive the reduction of lateral plate armor in fresh water while maintaining armor in the sea. We measured two fitness-related traits, hatching success and juvenile growth rate, in offspring of reduced (low and partial) and complete lateral plate morphs from two polymorphic populations when raised in either fresh water or salt water. In contrast to Heuts' results, there was little difference among morphs in hatching success. However, salinity strongly influenced juvenile growth: offspring of reduced lateral plate morphs grew substantially faster (up to 65%) than offspring of completely plated morphs in fresh water, but there was little difference in salt water. We suggest that the parallel loss of lateral plates in fresh water has arisen through a correlated response to selection for faster growth during lateral plate development, but the effect of salinity on hatching success and juvenile growth rate cannot explain the predominance of completely plated morphs in marine populations.     

Modular Skeletal Evolution in Sticklebacks Is Controlled by Additive and Clustered Quantitative Trait Loci

Understanding the genetic architecture of evolutionary change remains a long-standing goal in biology. In vertebrates, skeletal evolution has contributed greatly to adaptation in body form and function in response to changing ecological variables like diet and predation. Here we use genome-wide linkage mapping in threespine stickleback fish to investigate the genetic architecture of evolved changes in many armor and trophic traits. We identify >100 quantitative trait loci (QTL) controlling the pattern of serially repeating skeletal elements, including gill rakers, teeth, branchial bones, jaws, median fin spines, and vertebrae. We use this large collection of QTL to address long-standing questions about the anatomical specificity, genetic dominance, and genomic clustering of loci controlling skeletal differences in evolving populations. We find that most QTL (76%) that influence serially repeating skeletal elements have anatomically regional effects. In addition, most QTL (71%) have at least partially additive effects, regardless of whether the QTL controls evolved loss or gain of skeletal elements. Finally, many QTL with high LOD scores cluster on chromosomes 4, 20, and 21. These results identify a modular system that can control highly specific aspects of skeletal form. Because of the general additivity and genomic clustering of major QTL, concerted changes in both protective armor and trophic traits may occur when sticklebacks inherit either marine or freshwater alleles at linked or possible “supergene” regions of the stickleback genome. Further study of these regions will help identify the molecular basis of both modular and coordinated changes in the vertebrate skeleton.     

Adapting to winter in wheat: a long-term study follows parallel phenotypic and genetic changes in three experimental wheat populations

Drawing a direct connection between adaptive evolution at the phenotypic level and underlying genetic factors has long been a major goal of evolutionary biologists, but the genetic characterization of adaptive traits in natural populations is notoriously difficult. The study of evolution in experimental populations offers some help — initial conditions are known and changes can be tracked for extended periods under conditions more controlled than wild populations and more realistic than laboratory or greenhouse experiments. In this issue of Molecular Ecology , researchers studying experimental wheat populations over a 12-year period have demonstrated evolution in a major adaptive trait, flowering time, and parallel changes in underlying genetic variation ( Rhoné et al . 2008 ). Their work suggests that cis -regulatory mutations at a single gene may explain most of the flowering time variation in these populations.     

ENVIRONMENT SPECIFIC PLEIOTROPY FACILITATES DIVERGENCE AT THE ECTODYSPLASIN LOCUS IN THREESPINE STICKLEBACK

Adaptive radiation occurs when divergent natural selection in different environments leads to phenotypic differentiation. The pleiotropic effects of underlying genes can either promote or constrain this diversification. Identifying the pleiotropic effects of genes responsible for divergent traits, and testing how the environment influences these effects, can therefore help to provide an understanding of how ecology drives evolutionary change between populations. Positive selection on low-armor alleles at the Ectodysplasin ( Eda ) locus in threespine stickleback has led to the repeated evolution of reduced armor in populations following freshwater colonization by fully armored marine sticklebacks. Here, we demonstrate that Eda has environmentally determined pleiotropic effects on armor and growth. When raised in freshwater, reduced armor sticklebacks carrying "low" alleles at Eda had increased growth rate relative to fully armored sticklebacks carrying "complete" alleles. In saltwater treatments this growth advantage was present during juvenile growth but lost during adult growth, suggesting that in this environment stickleback are able to develop full armor plates without sacrificing overall growth rate. The environment specific pleiotropic effects of Eda demonstrate that ecological factors can mediate the influence of genetic architecture in driving phenotypic evolution. Furthermore, because size is important for mate choice in stickleback, the growth rate differences influenced by Eda may have effects on reproductive isolation between marine and freshwater populations.     

MULTIPLE EVOLUTIONARY PATHWAYS TO DECREASED LATERAL PLATE COVERAGE IN FRESHWATER THREESPINE STICKLEBACKS

Adaptation to novel environments can be based either on standing genetic variation or variation attributable to new mutations. When standing genetic variation for a functional adaptation is lacking, and variation due to new mutations is not yet available, adaptation is possible only through alternative functional solutions. Reduction in the number of bony lateral plates as an adaptation to freshwater colonization by marine threespine sticklebacks (Gasterosteus aculeatus) has occurred in numerous independent cases through allelic substitution in the ectodysplasin‐a (Eda) gene. Studying the phenotypic and genetic variation in plate number and size in five marine and six freshwater threespine stickleback populations, we found that when variation in Eda was limiting (i.e., alleles associated with the low‐plate morph were missing or in extremely low frequency), plate number reduction did not take place in freshwater populations, but reduced lateral plate coverage was achieved by a reduction in the size of lateral plates. Our results suggest that this phenotypically and genetically discrete "small‐plated" threespine stickleback—which is the dominant form in three northern European freshwater populations—may be functionally equivalent to the low‐plated morph and hence, serve as an example of convergent evolution toward functional similarity in the face of genetic constraints.     

Independent axes of genetic variation and parallel evolutionary divergence of opercle bone shape in threespine stickleback

Evolution of similar phenotypes in independent populations is often taken as evidence of adaptation to the same fitness optimum. However, the genetic architecture of traits might cause evolution to proceed more often toward particular phenotypes, and less often toward others, independently of the adaptive value of the traits. Freshwater populations of Alaskan threespine stickleback have repeatedly evolved the same distinctive opercle shape after divergence from an oceanic ancestor. Here we demonstrate that this pattern of parallel evolution is widespread, distinguishing oceanic and freshwater populations across the Pacific Coast of North America and Iceland. We test whether this parallel evolution reflects genetic bias by estimating the additive genetic variance-covariance matrix (G) of opercle shape in an Alaskan oceanic (putative ancestral) population. We find significant additive genetic variance for opercle shape and that G has the potential to be biasing, because of the existence of regions of phenotypic space with low additive genetic variation. However, evolution did not occur along major eigenvectors of G, rather it occurred repeatedly in the same directions of high evolvability. We conclude that the parallel opercle evolution is most likely due to selection during adaptation to freshwater habitats, rather than due to biasing effects of opercle genetic architecture.     

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.     

Parallel and nonparallel aspects of ecological, phenotypic, and genetic divergence across replicate population pairs of lake and stream stickleback

Parallel (or convergent) evolution provides strong evidence for a deterministic role of natural selection: similar phenotypes evolve when independent populations colonize similar environments. In reality, however, independent populations in similar environments always show some differences: some nonparallel evolution is present. It is therefore important to explicitly quantify the parallel and nonparallel aspects of trait variation, and to investigate the ecological and genetic explanations for each. We performed such an analysis for threespine stickleback (Gasterosteus aculeatus) populations inhabiting lake and stream habitats in six independent watersheds. Morphological traits differed in the degree to which lake-stream divergence was parallel across watersheds. Some aspects of this variation were correlated with ecological variables related to diet, presumably reflecting the strength and specifics of divergent selection. Furthermore, a genetic scan revealed some markers that diverged between lakes and streams in many of the watersheds and some that diverged in only a few watersheds. Moreover, some of the lake-stream divergence in genetic markers was associated within some of the lake-stream divergence in morphological traits. Our results suggest that parallel evolution, and deviations from it, are primarily the result of natural selection, which corresponds in only some respects to the dichotomous habitat classifications frequently used in such studies.     

A gene with major phenotypic effects as a target for selection vs. homogenizing gene flow

Genes with major phenotypic effects facilitate quantifying the contribution of genetic vs. plastic effects to adaptive divergence. A classical example is Ectodysplasin (Eda), the major gene controlling lateral plate phenotype in three‐spined stickleback. Completely plated marine stickleback populations evolved repeatedly towards low‐plated freshwater populations, representing a prime example of parallel evolution by natural selection. However, many populations remain polymorphic for lateral plate number. Possible explanations for this polymorphism include relaxation of selection, disruptive selection or a balance between divergent selection and gene flow. We investigated 15 polymorphic stickleback populations from brackish and freshwater habitats in coastal North‐western Europe. At each site, we tracked changes in allele frequency at the Eda gene between subadults in fall, adults in spring and juveniles in summer. Eda genotypes were also compared for body size and reproductive investment. We observed a fitness advantage for the Eda allele for the low morph in freshwater and for the allele for the complete morph in brackish water. Despite these results, the differentiation at the Eda gene was poorly correlated with habitat characteristics. Neutral population structure was the best predictor of spatial variation in lateral plate number, suggestive of a substantial effect of gene flow. A meta‐analysis revealed that the signature of selection at Eda was weak compared to similar studies in stickleback. We conclude that a balance between divergent selection and gene flow can maintain stickleback populations polymorphic for lateral plate number and that ecologically relevant genes may not always contribute much to local adaptation, even when targeted by selection.     

Asymmetry in structural defenses: insights into selective predation in the wild

Assessment of geographical patterns in fluctuating asymmetry (small, random differences between sides of bilateral characters) among populations shows promise as a tool to resolve the relative biomechanical importance of traits, in addition to being a possible indicator of habitat quality. We used 115 endemic freshwater populations of threespine stickleback (Gasterosteus aculeatus) from Haida Gwaii (Queen Charlotte Islands), British Columbia, Canada, to explore the degree of concordance between geographical variation of asymmetry in a predator defense structure (bony lateral plates) and geographical variation in several indirect measures of predation regime as well as several abiotic habitat variables. We found a geographical cline in the population frequency of lateral plate asymmetries, with reduced asymmetry in the southern clear-water regions of the archipelago characterized by long reaction distances and greater chance of capture by predators, and elevated asymmetry in the northern stained-water regions with poor visibility and low chances of capture. Lateral plate asymmetry was strongly correlated with expression of several defensive armor traits, including total plate numbers among populations, mean cross-sectional diameter of stickleback with the dorsal and pelvic spines erect, and mean degree of overlap between the plates and spine supports. There were no significant correlations between frequency of asymmetric fish and any of our abiotic habitat variables. Stickleback with structural plate asymmetries had fewer trout-induced scars than symmetric fish in the significant majority of populations, and there was a decrease in structural plate asymmetry with age in stained-water habitats, suggesting that trout predators may be selectively removing asymmetric fish in some lakes. This study provides evidence that geographical variation in developmental stability of threespine stickleback, as seen in the frequencies of asymmetry, reflects differences among populations in the importance of structural defenses to fitness rather than differences in habitat quality, and that asymmetry may be a target of selection by predators in wild populations.     

The utility of QTL-Linked markers to detect selective sweeps in natural populations--a case study of the EDA gene and a linked marker in threespine stickleback

Sequence polymorphisms in coding genes and variability in quantitative trait loci (QTL)-linked markers can be used to uncover the evolutionary mechanisms of traits involved in adaptive processes. We studied sequence variation in the EDA gene and allelic variation in 18 microsatellites - one of which (Gac4174) is linked with the EDA QTL - in low, partially and completely plated morphs from eight threespine stickleback European populations. The results agree with previous studies in that EDA polymorphism is closely related to plate number variation: EDA sequences grouped populations into low and completely plated morphs, whereas microsatellites failed to do so. Furthermore, partially plated fish were heterozygous with respect to the distinctive EDA alleles for completely and low plated morphs, indicating that completely plated morph alleles are not entirely dominant in controlling the expression of lateral plate number. An examination of population differentiation in plate number with quantitative genetic methods revealed that the degree of differentiation exceeded that expected from genetic drift alone (Q(ST) > F(ST)). Our results support the adaptive genetic differentiation of plate morphs and the view that distinctive EDA gene polymorphism occurs in similar sites across the distribution range of this species. Yet, allele frequency differentiation in the Gac4174 microsatellite locus, informative in experimental crosses for plate number variation, did not differ from that of neutral markers and, was therefore unable to detect the signature of natural selection responsible for population divergence.     

GENETIC ARCHITECTURE OF SKELETAL EVOLUTION IN EUROPEAN LAKE AND STREAM STICKLEBACK

Advances in genomic techniques are greatly facilitating the study of molecular signatures of selection in diverging natural populations. Connecting these signatures to phenotypes under selection remains challenging, but benefits from dissections of the genetic architecture of adaptive divergence. We here perform quantitative trait locus (QTL) mapping using 488 F₂ individuals and 2011 single nucleotide polymorphisms (SNPs) to explore the genetic architecture of skeletal divergence in a lake‐stream stickleback system from Central Europe. We find QTLs for gill raker, snout, and head length, vertebral number, and the extent of lateral plating (plate number and height). Although two large‐effect loci emerge, QTL effect sizes are generally small. Examining the neighborhood of the QTL‐linked SNPs identifies several genes involved in bone formation, which emerge as strong candidate genes for skeletal evolution. Finally, we use SNP data from the natural source populations to demonstrate that some SNPs linked to QTLs in our cross also exhibit striking allele frequency differences in the wild, suggesting a causal role of these QTLs in adaptive population divergence. Our study paves the way for comparative analyses across other (lake‐stream) stickleback populations, and for functional investigations of the candidate genes.     

Five anthocyanin polymorphisms are associated with an R2R3-MYB cluster in Mimulus guttatus (Phrymaceae)

? Premise of study: Botanists have long been interested in the reasons for genetic variation among individuals, populations, and species of plants. The anthocyanin pathway is ideal for studying the evolution of such phenotypic variation. ? Methods: We used a combination of quantitative trait loci mapping and association studies to understand the genetic basis of variation in five anthocyanin phenotypes including calyx, corolla, and leaf coloration patterns that vary within and among populations of Mimulus guttatus. We then examined what genes might be responsible for this phenotypic variation and whether one of the traits, calyx spotting, is randomly distributed across the geographic range of the species. ? Key results: All five phenotypes in M. guttatus were primarily controlled by the same major locus (PLA1), which contains a tandem array of three R2R3-MYB genes known to be involved in the evolution of flower color in a related species of Mimulus. Calyx spotting was nonrandomly distributed across the range of M. guttatus and correlated with multiple climate variables. ? Conclusions: The results of this study suggest that variation in R2R3-MYB genes is the primary cause of potentially important anthocyanin phenotypic variation within and among populations of M. guttatus, a finding consistent with recent theoretical and empirical research on flower color evolution.     

Parallel adaptation to climate above the 35th parallel

Independent or parallel evolution of similar traits is key to understanding the genetics and limitations of adaptation. Adaptation from the same genetic changes in different populations defines parallel evolution. Such genetic changes can derive from standing ancestral variation or de novo mutations and excludes instances of adaptive introgression. In this issue of Molecular Ecology, Walden et al.(2020) investigate the scale of parallel climate adaptation from standing genetic variation between two North American Arabidopsis lyrata lineages, each formed by a distinct evolutionary history during the last glacial cycle. By identifying adaptive variants correlated with three ecologically significant climatic gradients, they show that instead of the same genetic variants or even genes, parallel evolution is only observed at the level of biological processes. The evolution of independent adaptive variants to climate in two genetically close lineages is explained by their different post‐glacial demographic histories. Separate glacial refugia and strong population bottlenecks were probably sufficient to change the landscape of shared allele frequencies, hindering the possibility of parallel evolution.