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.     

Allometric relationships in morphological traits associated with foraging,swimming ability,and predator defense reveal adaptations toward brackish and freshwater environments in the threespine stickleback

Freshwater colonization by threespine stickleback has led to divergence in morphology between ancestral marine and derived freshwater populations, making them ideal for studying natural selection on phenotypes. In an open brackish–freshwater system, we previously discovered two genetically distinct stickleback populations that also differ in geometric shape: one mainly found in the brackish water lagoon and one throughout the freshwater system. As shape and size are not perfectly correlated, the aim of this study was to identify the morphological trait(s) that separated the populations in geometric shape. We measured 23 phenotypes likely to be important for foraging, swimming capacity, and defense against predation. The lateral plate morphs in freshwater displayed few significant changes in trait sizes, but the low plated expressed feeding traits more associated with benthic habitats. When comparing the completely plated genetically assigned populations, the freshwater, the hybrids, the migrants and the lagoon fish, many of the linear traits had different slopes and intercepts in trait‐size regressions, precluding our ability to directly compare all traits simultaneously, which most likely results from low variation in body length for the lagoon and migrant population. We found the lagoon stickleback population to be more specialized toward the littoral zone, displaying benthic traits such as large, deep bodies with smaller eyes compared to the freshwater completely plated morph. Further, the lagoon and migrant fish had an overall higher body coverage of lateral plates compared to freshwater fish, and the dorsal and pelvic spines were longer. Evolutionary constraints due to allometric scaling relationships could explain the observed, overall restricted, differences in morphology between the sticklebacks in this study, as most traits have diversified in common allometric trajectories. The observed differences in foraging and antipredation traits between the fish with a lagoon and freshwater genetic signature are likely a result of genetic or plastic adaptations toward brackish and freshwater environments.     

Two morphological types in the reproductive stock of three-spined stickleback,Gasterosteus aculeatus,in Lake Harutori,Hokkaido Island

Synopsis Two morphological types of the trachurus form (completely plated morph) of three-spined stickleback, Gasterosteus aculeatus, are found in Lake Harutori, Kushiro, east district of Hokkaido Island. For convenience, the two forms are referred to as ‘large type’ and ‘small type’ on the basis of body length (discrimination at 70 mm in length). The two types were examined for morphometric, meristic and reproductive characteristics. They differ in body length, the number of lateral plates, means of the 2nd dorsal spine length/body length and of the pelvic spine length/body length, the relationship between body length and head length, between body length and gonad weight, between body weight and gonad weight, and between body length and ovarian egg number; significant differences were present for each sex. These two types were compared with the anadromous stickleback migrating into the freshwater area near Lake Harutori to breed. The anadromous fish was morphologically much more similar to the large type than the small type. It is suggested that the large type is also an anadromous form and the small type is the permanent freshwater resident.     

Lateral plate number in low‐plated threespine stickleback: a study of plasticity and heritability

In the threespine stickleback Gasterosteus aculeatus model system, phenotypes are often classified into three morphs according to lateral plate number. Morph identity has been shown to be largely genetically determined, but substantial within‐morph variation in plate number exists. In this study, we test whether plate number has a plastic component in response to salinity in the low‐plated morph using a split‐clutch experiment where families were split in two, one half raised in water at 0 and the other at 30 ppt salt. We find a small salinity‐induced plastic effect on plate number in an unexpected direction, opposite to what we predicted: Fish raised in freshwater on average have slightly more plates than fish raised in saltwater. Our results confirm that heritability of plate number is high. Additionally, we find that variance in plate number at the family level can be predicted from other family level traits, which might indicate that epistatic interactions play a role in creating the observed pattern of lateral plate number variation.     

The role of calcium and predation on plate morph evolution in the three‐spined stickleback (Gasterosteus aculeatus)

While the genetic basis to plate morph evolution of the three‐spined stickleback (Gasterosteus aculeatus) is well described, the environmental variables that select for different plate and spine morphs are incompletely understood. Using replicate populations of three‐spined sticklebacks on North Uist, Scotland, we previously investigated the role of predation pressure and calcium limitation on the adaptive evolution of stickleback morphology and behavior. While dissolved calcium proved a significant predictor of plate and spine morph, predator abundance did not. Ecol. Evol., xxx, 2014 and xxx performed a comparable analysis to our own to address the same question. They failed to detect a significant effect of dissolved calcium on morphological evolution, but did establish a significant effect of predation; albeit in the opposite direction to their prediction.     

Twelve years of contemporary armor evolution in a threespine stickleback population

Abstract —Loberg Lake, Alaska was colonized by sea-run Gasterosteus aculeatus between 1983 and 1988, after the original stickleback population was exterminated. Annual samples from 1990 to 2001 reveal substantial evolution of lateral plate (armor) phenotypes. The 1990 sample was nearly monomorphic for the complete plate morph, which is monomorphic in local sea-run populations; the low plate morph, which is usually monomorphic in local freshwater populations, was absent. By 2001, the frequency of completes had declined to 11%, and lows had increased to 75%. The partial plate morph and two unusual intermediate plate phenotypes were generally rare, but occurrence of the intermediates was unexpected. These intermediate phenotypes rarely occur in other, presumably older, polymorphic populations. When low morphs first appeared, they averaged 6.8 plates per side, indicating that the ancestral plate number of low morphs is high, and their mean has subsequently declined. Contemporary evolution in this population indicates that threespine stickleback adapt to freshwater habitats within decades after invasion from the ocean, and thus phenotypes in most populations are adapted to current conditions.     

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.     

Calcium and salinity as selective factors in plate morph evolution of the three-spined stickleback (Gasterosteus aculeatus)

Identifying the causal factors underlying natural selection remains a key challenge in evolutionary biology. Although the genetic basis for the plate morph evolution of three-spined stickleback (Gasterosteus aculeatus) is well described, the environmental variables that form the basis for different morphs are not understood. We measured the effects of dissolved calcium and salinity on the growth of sticklebacks with different plate morphs from Scotland and Poland. There was a significant interaction of calcium with plate morph for fish from both regions, with complete morph sticklebacks growing more slowly at low calcium concentrations and low morph sticklebacks showing divergent responses to calcium concentration. A Scottish anadromous population showed evidence of local adaptation to high salinity, which was independent of plate morph. Polish and Scottish populations diverged in their response to salinity, suggesting a difference in osmotic regulation. The results implicate a role for calcium in selecting for plate morph evolution in sticklebacks, possibly as a limiting element in skeletal growth.     

Rapid Morphological Changes in Threespine Stickleback, Gasterosteus aculeatus, in Freshwater

Synopsis Freshwater and marine threespine stickleback, Gasterosteus aculeatus, differ remarkably in armour plate number and body shape, although differences in other morphological characters are also common. Most freshwater populations have apparently evolved after isolation of marine sticklebacks in freshwater. After colonisation of freshwater habitats, they show rapid morphological changes and associated genetic isolation within as few as eight generations. I transferred fish from marine tide pools to two isolated freshwater ponds, differing in habitat characteristics, at the beginning of the breeding season, when females had ripe ovaries and males had breeding coloration. The first generation fish that I sampled from the ponds had significantly fewer armour plates than their marine ancestors and differed in shape. I also found some significant differences between fish sampled from the larger pond and those from a smaller, adjacent pond. This extremely rapid morphological divergence suggests that either the marine sticklebacks were highly phenotypically plastic or that there was very strong natural selection acting on the first generation within freshwater habitats.     

Parallel evolution? Microsatellite variation of recently isolated marine and freshwater three-spined stickleback

Variation at nine microsatellite loci, four of which are linked to phenotypic traits (spine length and lateral plate morphology) in Canadian three-spined stickleback Gasterosteus aculeatus , are used to test for selection on marine and freshwater three-spined stickleback morphs in Iceland. There are indications of strong selection on loci linked to dorsal spine length, providing another potential example of parallel divergence at the genomic level in three-spined stickleback.     

Parasites and parallel divergence of the number of individual MHC alleles between sympatric three‐spined stickleback Gasterosteus aculeatus morphs in Iceland

Two pairs of sympatric three‐spined stickleback Gasterosteus aculeatus morphs and two single morph populations inhabiting mud and lava or rocky benthic habitats in four Icelandic lakes were screened for parasites and genotyped for MHC class IIB diversity. Parasitic infection differed consistently between G. aculeatus from different benthic habitats. Gasterosteus aculeatus from the lava or rocky habitats were more heavily infected in all lakes. A parallel pattern was also found in individual MHC allelic variation with lava G. aculeatus morphs exhibiting lower levels of variation than the mud morphs. Evidence for selective divergence in MHC allele number is ambiguous but supported by two findings in addition to the parallel pattern observed. MHC allele diversity was not consistent with diversity reported at neutral markers (microsatellites) and in Þingvallavatn the most common number of alleles in each morph was associated with lower infection levels. In the Þingvallavatn lava morph, lower infection levels by the two most common parasites, Schistocephalus solidus and Diplostomum baeri, were associated with different MHC allele numbers.     

Ontogeny of sexual dimorphism and phenotypic integration in heritable morphs

In this study we investigated the developmental basis of adult phenotypes in a non-model organism, a polymorphic damselfly (Ischnura elegans) with three female colour morphs. This polymorphic species presents an ideal opportunity to study intraspecific variation in growth trajectories, morphological variation in size and shape during the course of ontogeny, and to relate these juvenile differences to the phenotypic differences of the discrete adult phenotypes; the two sexes and the three female morphs. We raised larvae of different families in individual enclosures in the laboratory, and traced morphological changes during the course of ontogeny. We used principal components analysis to examine the effects of Sex, Maternal morph, and Own morph on body size and body shape. We also investigated the larval fitness consequences of variation in size and shape by relating these factors to emergence success. Females grew faster than males and were larger as adults, and there was sexual dimorphism in body shape in both larval and adult stages. There were also significant effects of both maternal morph and own morph on growth rate and body shape in the larval stage. There were significant differences in body shape, but not body size, between the adult female morphs, indicating phenotypic integration between colour, melanin patterning, and body shape. Individuals that emerged successfully grew faster and had different body shape in the larval stage, indicating internal (non-ecological) selection on larval morphology. Overall, morphological differences between individuals at the larval stage carried over to the adult stage. Thus, selection in the larval stage can potentially result in correlated responses in adult phenotypes and vice versa.     

Water temperature,not fish morph,determines parasite infections of sympatric Icelandic threespine sticklebacks (Gasterosteus aculeatus)

Parasite communities of fishes are known to respond directly to the abiotic environment of the host, for example, to water quality and water temperature. Biotic factors are also important as they affect the exposure profile through heterogeneities in parasite distribution in the environment. Parasites in a particular environment may pose a strong selection on fish. For example, ecological differences in selection by parasites have been hypothesized to facilitate evolutionary differentiation of freshwater fish morphs specializing on different food types. However, as parasites may also respond directly to abiotic environment the parasite risk does not depend only on biotic features of the host environment. It is possible that different morphs experience specific selection gradients by parasites but it is not clear how consistent the selection is when abiotic factors change. We examined parasite pressure in sympatric morphs of threespine stickleback (Gasterosteus aculeatus) across a temperature gradient in two large Icelandic lakes, Myvatn and Thingvallavatn. Habitat‐specific temperature gradients in these lakes are opposite. Myvatn lava rock morph lives in a warm environment, while the mud morph lives in the cold. In Thingvallavatn, the lava rock morph lives in a cold environment and the mud morph in a warm habitat. We found more parasites in fish living in higher temperature in both lakes, independent of the fish morph, and this pattern was similar for the two dominating parasite taxa, trematodes and cestodes. However, at the same time, we also found higher parasite abundance in a third morph living in deep cold–water habitat in Thingvallavatn compared to the cold‐water lava morph, indicating strong effect of habitat‐specific biotic factors. Our results suggest complex interactions between water temperature and biotic factors in determining the parasite community structure, a pattern that may have implications for differentiation of stickleback morphs.     

Body armour and lateral‐plate reduction in freshwater three‐spined stickleback Gasterosteus aculeatus: adaptations to a different buoyancy regime?

Several factors related to buoyancy were compared between one marine and two freshwater populations of three‐spined stickleback Gasterosteus aculeatus. Fish from all three populations had buoyancy near to neutral to the ambient water. This showed that neither marine nor freshwater G. aculeatus used swimming and hydrodynamic lift to prevent sinking. Comparing the swimbladder volumes showed that freshwater completely plated G. aculeatus had a significantly larger swimbladder volume than both completely plated marine and low‐plated freshwater G. aculeatus. Furthermore, body tissue density was lower in low‐plated G. aculeatus than in the completely plated marine and freshwater fish. The results show that G. aculeatus either reduce tissue density or increase swimbladder volume to adapt to lower water density. Mass measurements of lateral plates and pelvis showed that loss of body armour in low‐plated G. aculeatus could explain the tissue density difference between low‐plated and completely plated G. aculeatus. This suggests that the common occurrence of plate and armour reduction in freshwater G. aculeatus populations can be an adaptation to a lower water density.     

Genomic signatures of the plateless phenotype in the threespine stickleback

Understanding the genetic basis of traits involved in adaptive divergence and speciation is one of the most fundamental objectives in evolutionary biology. Toward that end, we look for signatures of extreme plate loss in the genome of freshwater threespine sticklebacks (Gasterosteus aculeatus). Plateless stickleback have been found in only a few lakes and streams across the world; they represent the far extreme of a phenotypic continuum (plate number) that has been studied for years, although plateless individuals have not yet been the subject of much investigation. We use a dense single nucleotide polymorphism dataset made using RADseq to study fish from three freshwater populations containing plateless and low plated individuals, as well as fish from full plated marine populations. Analyses were performed using FastStructure, sliding windows F_ST, Bayescan and latent factor mixed models to search for genomic differences between the low plated and plateless phenotypes both within and among the three lakes. At least 18 genomic regions which may contribute to within‐morph plate number variation were detected in our low plated stickleback populations. We see no evidence of a selective sweep between low and plateless fish; rather reduction of plate number within the low plated morph seems to be polygenic.     

Multivariate phenotypes and the potential for alternative phenotypic optima in wall lizard (Podarcis muralis) ventral colour morphs

A major goal in evolutionary biology is to determine how phenotypic variation arises and is maintained in natural populations. Recent studies examining the morphological, physiological and behavioural differences among discrete colour morphotypes (morphs) have revealed several mechanisms that maintain discrete variation within populations, including frequency‐dependence, density‐dependence and correlational selection. For example, trade‐offs over resource allocation to morphological, physiological and behavioural traits can drive correlational selection for morph‐specific phenotypic optima. Here, we describe a ventral colour polymorphism in the wall lizard (Podarcis muralis) and test the hypothesis that morphs differ along multivariate axes defined by trade‐offs in morphological, physiological, and immunological traits. We show that ventral colour is a discrete trait and that morphs differ in body size, prevalence of infection by parasites and infection intensity. We also find that morphs differ along multivariate phenotypic axes and experience different multivariate selection pressures. Our results suggest that multivariate selection pressures may favour alternative optimal morph‐specific phenotypes in P. muralis.     

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.     

Armor morphology and reproductive output in threespine stickleback,Gasterosteus aculeatus

Synopsis The threespine stickleback, Gasterosteus aculeatus, is an extensively armored fish inhabiting both marine and fresh waters across its holarctic distribution. Marine fish nearly always possess a full complement of bony lateral plates running from just behind the head to the tail, and a robust pelvic girdle complex. These armor features appear to constrain lateral and ventral abdominal distention, and therefore clutch volume. Freshwater populations in many areas exhibit variable reduction in lateral plate number, and in some regions the pelvic girdle is also reduced or lost. Freshwater populations also vary in the degree of abdominal distention exhibited by gravid females. We tested whether reduction in armoring might be correlated with increased clutch volume using five populations from the Cook Inlet area of Alaska. The hypothesis that populations having reduced pelvic girdle complexes would have greater size-adjusted clutch volumes was not supported. In fact, our two full-pelvic populations as a group had larger volumes. Similarly, size-adjusted clutch volumes were not related to pelvic phenotype within either of our two pelvic-reduced populations, nor to lateral plate morph within a fifth population. Other factors that may explain the interpopulation differences in clutch volume in threespine stickleback include body shape, food quantity and quality, intensity of predation, and even behavior. Except for a preliminary analysis of body shape, these possibilities remain unexplored. The concept of phenotypic integration suggests that these factors should be analyzed as a suite rather than individually.     

Colour polymorphism in relation to population dynamics of the leaf beetle, Chrysomela lapponica

We studied colour morph diversity and frequencies of light and dark morphs in non-fluctuating and fluctuating populations of willow feeding leaf beetle Chrysomela lapponica in the Kola Peninsula, NW Russia. Population-specific Shannon–Weaver diversity index positively correlated with dark morph frequencies, indicating that the larger part of colour polymorphism is related with numbers and diversity of dark morphs. Among-population variation in studied characters was not explained by pollution load or predation rates, but depended on the type of the population and the stage of density change in the fluctuating populations: both colour morph diversity and frequency of dark morphs were low in declining post-outbreak populations but equally high in non-fluctuating populations and in fluctuating populations at peak densities. In time-series, both diversity index and frequency of dark morphs decreased with post-outbreak density decline in the fluctuating population, but did not change in the non-fluctuating population. In the experiment, when adults received low quality food (plants from post-outbreak site), mortality of dark morphs during the hibernation was almost doubled relative to the mortality of light morphs, whereas on high quality food the colour morphs demonstrated similar mortality. This may indicate, that decrease in colour polymorphism extent and dark morph frequencies in the declining populations is due to selective mortality of dark morphs imposed by density dependent (induced by heavy herbivore damage during an outbreak) decrease in host-plant quality (delayed inducible resistance, DIR). DIR is known as one of the factors driving herbivore populations, but our result is the first evidence that DIR may act as a factor of natural selection. Dark morphs are not only susceptible to low food quality, but also have smaller size compared to light morphs, and therefore the dark females are presumably less fecund. Thus, decrease in frequency of low-fitness (dark) individuals in post-outbreak populations and accumulation of low-fitness phenotypes at the popu-lation peak may create feedbacks contributing to regulation of density fluctuations in Ch. lapponica.     

Non‐parallel divergence across freshwater and marine three‐spined stickleback Gasterosteus aculeatus populations

This work investigated whether multiple freshwater populations of three‐spined stickleback Gasterosteus aculeatus in different freshwater catchments in the Jutland Peninsula, Denmark, derived from the same marine populations show repeated adaptive responses. A total of 327 G. aculeatus collected at 13 sampling locations were screened for genetic variation using a combination of 70 genes putatively under selection and 26 neutral genes along with a marker linked to the ectodysplasin gene (eda), which is strongly correlated with plate armour morphs in the species. A highly significant genetic differentiation was found that was higher among different freshwater samples than between marine–freshwater samples. Tests for selection between marine and freshwater populations showed a very low degree of parallelism and no single nucleotide polymorphism was detected as outlier in all freshwater–marine pairwise comparisons, including the eda. This suggests that G. aculeatus is not necessarily the prime example of parallel local adaptation suggested in much of the literature and that important exceptions exist (i.e. the Jutland Peninsula). While marine populations in the results described here showed a high phenotype–genotype correlation at eda, a low association was found for most of the freshwater populations. The most extreme case was found in the freshwater Lake Hald where all low‐plated phenotypes were either homozygotes for the allele supposed to be associated with completely plated morphs or heterozygotes, but none were homozygotes for the putative low‐plated allele. Re‐examination of data from seven G. aculeatus studies agrees in showing a high but partial association between phenotype–genotype at eda in G. aculeatus freshwater populations and that mismatches occur everywhere in the European regions studied (higher in some areas, i.e. Denmark). This is independent of the eda marker used.