Reproductive isolation in a three-way contact zone.

Contact zones between divergent forms within a species provide insight into the role of gene flow in adaptation and speciation. Previous work has focused on contact zones involving only two divergent forms, but in nature, many more than two populations may overlap simultaneously and experience gene flow. Patterns of introgression in wild populations are, therefore, likely much more complicated than is often assumed. We begin to address this gap in current knowledge by investigating patterns of divergence and introgression across a complex natural contact zone. We use phenotypic and genomic data to confirm the existence of a three-way contact zone among divergent freshwater resident, saltwater resident and saltwater migratory three-spined stickleback (Gasterosteus aculeatus) on the island of North Uist, Scottish Western Isles. We find evidence for hybridization, mostly between saltwater resident and saltwater migratory forms. Despite hybridization, genomic analyses reveal pairwise islands of divergence between all forms that are maintained across the contact zone. Genomic cline analyses also provide evidence for selection and/or hybrid incompatibilities in divergent regions. Divergent genomic regions occur across multiple chromosomes and involve many known adaptive loci and several chromosomal inversions. We also identify distinct immune gene expression profiles between forms, but no evidence for transgressive expression in hybrids. Our results suggest that reproductive isolation is maintained in this three-way contact zone, despite some hybridization, and that reduced recombination in chromosomal inversions may play an important role in maintaining this isolation.

On the Origins of Phenotypic Parallelism in Benthic and Limnetic Stickleback.

Rapid evolution of similar phenotypes in similar environments, giving rise to in situ parallel adaptation, is an important hallmark of ecological speciation. However, what appears to be in situ adaptation can also arise by dispersal of divergent lineages from elsewhere. We test whether two contrasting phenotypes repeatedly evolved in parallel, or have a single origin, in an archetypal example of ecological adaptive radiation: benthic-limnetic three-spined stickleback (Gasterosteus aculeatus) across species pair and solitary lakes in British Columbia. We identify two genomic clusters across freshwater populations, which differ in benthic-limnetic divergent phenotypic traits and separate benthic from limnetic individuals in species pair lakes. Phylogenetic reconstruction and niche evolution modeling both suggest a single evolutionary origin for each of these clusters. We detected strong phylogenetic signal in benthic-limnetic divergent traits, suggesting that they are ancestrally retained. Accounting for ancestral state retention, we identify local adaptation of body armor due to the presence of an intraguild predator, the sculpin (Cottus asper), and environmental effects of lake depth and pH on body size. Taken together, our results imply a predominant role for retention of ancestral characteristics in driving trait distribution, with further selection imposed on some traits by environmental factors.

Ecological constraint, rather than opportunity, promotes adaptive radiation in three-spined stickleback (Gasterosteus aculeatus) on North Uist.

The context and cause of adaptive radiations have been widely described and explored but why rapid evolutionary diversification does not occur in related evolutionary lineages has yet to be understood. The standard answer is that evolutionary diversification is provoked by ecological opportunity and that some lineages do not encounter the opportunity. Three-spined sticklebacks on the Scottish island of North Uist show enormous diversification, which seems to be associated with the diversity of aquatic habitats. Sticklebacks on the neighboring island of South Uist have not been reported to show the same level of evolutionary diversity, despite levels of environmental variation that we might expect to be similar to North Uist. In this study, we compared patterns of morphological and environmental diversity on North and South Uist. Ancestral anadromous sticklebacks from both islands exhibited similar morphology including size and bony "armor." Resident sticklebacks showed significant variation in armor traits in relation to pH of water. However, North Uist sticklebacks exhibited greater diversity of morphological traits than South Uist and this was associated with greater diversity in pH of the waters of lochs on North Uist. Highly acidic and highly alkaline freshwater habitats are missing, or uncommon, on South Uist. Thus, pH appears to act as a causal factor driving the evolutionary diversification of stickleback in local adaptation in North and South Uist. This is consistent with diversification being more associated with ecological constraint than ecological opportunity.

The maintenance of standing genetic variation: Gene flow vs. selective neutrality in Atlantic stickleback fish.

Adaptation to derived habitats often occurs from standing genetic variation. The maintenance within ancestral populations of genetic variants favourable in derived habitats is commonly ascribed to long-term antagonism between purifying selection and gene flow resulting from hybridization across habitats. A largely unexplored alternative idea based on quantitative genetic models of polygenic adaptation is that variants favoured in derived habitats are neutral in ancestral populations when their frequency is relatively low. To explore the latter, we first identify genetic variants important to the adaptation of threespine stickleback fish (Gasterosteus aculeatus) to a rare derived habitat-nutrient-depleted acidic lakes-based on whole-genome sequence data. Sequencing marine stickleback from six locations across the Atlantic Ocean then allows us to infer that the frequency of these derived variants in the ancestral habitat is unrelated to the likely opportunity for gene flow of these variants from acidic-adapted populations. This result is consistent with the selective neutrality of derived variants within the ancestor. Our study thus supports an underappreciated explanation for the maintenance of standing genetic variation, and calls for a better understanding of the fitness consequences of adaptive variation across habitats and genomic backgrounds.

Courtship behavior, nesting microhabitat, and assortative mating in sympatric stickleback species pairs.

The maintenance of reproductive isolation in the face of gene flow is a particularly contentious topic, but differences in reproductive behavior may provide the key to explaining this phenomenon. However, we do not yet fully understand how behavior contributes to maintaining species boundaries. How important are behavioral differences during reproduction? To what extent does assortative mating maintain reproductive isolation in recently diverged populations and how important are "magic traits"? Assortative mating can arise as a by-product of accumulated differences between divergent populations as well as an adaptive response to contact between those populations, but this is often overlooked. Here we address these questions using recently described species pairs of three-spined stickleback (Gasterosteus aculeatus), from two separate locations and a phenotypically intermediate allopatric population on the island of North Uist, Scottish Western Isles. We identified stark differences in the preferred nesting substrate and courtship behavior of species pair males. We showed that all males selectively court females of their own ecotype and all females prefer males of the same ecotype, regardless of whether they are from species pairs or allopatric populations. We also showed that mate choice does not appear to be driven by body size differences (a potential "magic trait"). By explicitly comparing the strength of these mating preferences between species pairs and single-ecotype locations, we were able to show that present levels of assortative mating due to direct mate choice are likely a by-product of other adaptations between ecotypes, and not subject to obvious selection in species pairs. Our results suggest that ecological divergence in mating characteristics, particularly nesting microhabitat may be more important than direct mate choice in maintaining reproductive isolation in stickleback species pairs.

Intercontinental genomic parallelism in multiple three-spined stickleback adaptive radiations.

Parallelism, the evolution of similar traits in populations diversifying in similar conditions, provides strong evidence of adaptation by natural selection. Many studies of parallelism focus on comparisons of different ecotypes or contrasting environments, defined a priori, which could upwardly bias the apparent prevalence of parallelism. Here, we estimated genomic parallelism associated with components of environmental and phenotypic variation at an intercontinental scale across four freshwater adaptive radiations (Alaska, British Columbia, Iceland and Scotland) of the three-spined stickleback (Gasterosteus aculeatus). We combined large-scale biological sampling and phenotyping with restriction site associated DNA sequencing (RAD-Seq) data from 73 freshwater lake populations and four marine ones (1,380 fish) to associate genome-wide allele frequencies with continuous distributions of environmental and phenotypic variation. Our three main findings demonstrate that (1) quantitative variation in phenotypes and environments can predict genomic parallelism; (2) genomic parallelism at the early stages of adaptive radiations, even at large geographic scales, is founded on standing variation; and (3) similar environments are a better predictor of genome-wide parallelism than similar phenotypes. Overall, this study validates the importance and predictive power of major phenotypic and environmental factors likely to influence the emergence of common patterns of genomic divergence, providing a clearer picture than analyses of dichotomous phenotypes and environments.

Prior exposure to long-day photoperiods alters immune responses and increases susceptibility to parasitic infection in stickleback.

Seasonal disease and parasitic infection are common across organisms, including humans, and there is increasing evidence for intrinsic seasonal variation in immune systems. Changes are orchestrated through organisms' physiological clocks using cues such as day length. Ample research in diverse taxa has demonstrated multiple immune responses are modulated by photoperiod, but to date, there have been few experimental demonstrations that photoperiod cues alter susceptibility to infection. We investigated the interactions among photoperiod history, immunity and susceptibility in laboratory-bred three-spined stickleback (a long-day breeding fish) and its external, directly reproducing monogenean parasite Gyrodactylus gasterostei. We demonstrate that previous exposure to long-day photoperiods (PLD) increases susceptibility to infection relative to previous exposure to short days (PSD), and modifies the response to infection for the mucin gene muc2 and Treg cytokine foxp3a in skin tissues in an intermediate 12 L : 12 D photoperiod experimental trial. Expression of skin muc2 is reduced in PLD fish, and negatively associated with parasite abundance. We also observe inflammatory gene expression variation associated with natural inter-population variation in resistance, but find that photoperiod modulation of susceptibility is consistent across host populations. Thus, photoperiod modulation of the response to infection is important for host susceptibility, highlighting new mechanisms affecting seasonality of host-parasite interactions.

Relationships between immune gene expression and circulating cytokine levels in wild house mice.

Quantitative PCR (qPCR) has been commonly used to measure gene expression in a number of research contexts, but the measured RNA concentrations do not always represent the concentrations of active proteins which they encode. This can be due to transcriptional regulation or post-translational modifications, or localization of immune environments, as can occur during infection. However, in studies using free-living non-model species, such as in ecoimmunological research, qPCR may be the only available option to measure a parameter of interest, and so understanding the quantitative link between gene expression and associated effector protein levels is vital.Here, we use qPCR to measure concentrations of RNA from mesenteric lymph node (MLN) and spleen tissue, and multiplex ELISA of blood serum to measure circulating cytokine concentrations in a wild population of a model species, Mus musculus domesticus.Few significant correlations were found between gene expression levels and circulating cytokines of the same immune genes or proteins, or related functional groups. Where significant correlations were observed, these were most frequently within the measured tissue (i.e., the expression levels of genes measured from spleen tissue were more likely to correlate with each other rather than with genes measured from MLN tissue, or with cytokine concentrations measured from blood).Potential reasons for discrepancies between measures including differences in decay rates and transcriptional regulation networks are discussed. We highlight the relative usefulness of different measures under different research questions and consider what might be inferred from immune assays.

Immune state is associated with natural dietary variation in wild mice Mus musculus domesticus.

The ability, propensity and need to mount an immune response vary both among individuals and within a single individual over time.A wide array of parameters has been found to influence immune state in carefully controlled experiments, but we understand much less about which of these parameters are important in determining immune state in wild populations.Diet can influence immune responses, for example when nutrient availability is limited. We therefore predict that natural dietary variation will play a role in modulating immune state, but this has never been tested.We measured carbon and nitrogen stable isotope ratios in an island population of house mice Mus musculus domesticus as an indication of dietary variation, and the expression of a range of immune-related genes to represent immune state.After accounting for potential confounding influences such as age, sex and helminth load, we found a significant association between carbon isotope ratio and levels of immune activity in the mesenteric lymph nodes, particularly in relation to the inflammatory response.This association demonstrates the important interplay between diet and an animal's response to immune challenges, and therefore potentially its susceptibility to disease. A plain language summary is available for this article.

Geographical location influences the composition of the gut microbiota in wild house mice (Mus musculus domesticus) at a fine spatial scale.

The composition of the mammalian gut microbiota can be influenced by a multitude of environmental variables such as diet and infections. Studies investigating the effect of these variables on gut microbiota composition often sample across multiple separate populations and habitat types. In this study we explore how variation in the gut microbiota of the house mouse (Mus musculus domesticus) on the Isle of May, a small island off the east coast of Scotland, is associated with environmental and biological factors. Our study focuses on the effects of environmental variables, specifically trapping location and surrounding vegetation, as well as the host variables sex, age, body weight and endoparasite infection, on the gut microbiota composition across a fine spatial scale in a freely interbreeding population. We found that differences in gut microbiota composition were significantly associated with the trapping location of the host, even across this small spatial scale. Sex of the host showed a weak association with microbiota composition. Whilst sex and location could be identified as playing an important role in the compositional variation of the gut microbiota, 75% of the variation remains unexplained. Whereas other rodent studies have found associations between gut microbiota composition and age of the host or parasite infections, the present study could not clearly establish these associations. We conclude that fine spatial scales are important when considering gut microbiota composition and investigating differences among individuals.

Admixture between Ancient Lineages, Selection, and the Formation of Sympatric Stickleback Species-Pairs.

Ecological speciation has become a popular model for the development and maintenance of reproductive isolation in closely related sympatric pairs of species or ecotypes. An implicit assumption has been that such pairs originate (possibly with gene flow) from a recent, genetically homogeneous ancestor. However, recent genomic data have revealed that currently sympatric taxa are often a result of secondary contact between ancestrally allopatric lineages. This has sparked an interest in the importance of initial hybridization upon secondary contact, with genomic reanalysis of classic examples of ecological speciation often implicating admixture in speciation. We describe a novel occurrence of unusually well-developed reproductive isolation in a model system for ecological speciation: the three-spined stickleback (Gasterosteus aculeatus), breeding sympatrically in multiple lagoons on the Scottish island of North Uist. Using morphological data, targeted genotyping, and genome-wide single-nucleotide polymorphism data, we show that lagoon resident and anadromous ecotypes are strongly reproductively isolated with an estimated hybridization rate of only ∼1%. We use palaeoecological and genetic data to test three hypotheses to explain the existence of these species-pairs. Our results suggest that recent, purely ecological speciation from a genetically homogeneous ancestor is probably not solely responsible for the evolution of species-pairs. Instead, we reveal a complex colonization history with multiple ancestral lineages contributing to the genetic composition of species-pairs, alongside strong disruptive selection. Our results imply a role for admixture upon secondary contact and are consistent with the recent suggestion that the genomic underpinning of ecological speciation often has an older, allopatric origin.

Internal embryonic development in a non-copulatory, egg-laying teleost, the three-spined stickleback, Gasterosteus aculeatus.

The switch from egg-laying to retaining and giving birth to live young is a major transition in the history of life. Despite its repeated evolution across the fishes, records of intermediate phenotypes are vanishingly rare, with only two known cases in existence of normally egg-laying fish species retaining embryos within the ovaries. We report the discovery of a third occurrence, in which well-developed embryos were found in the ovaries of a three-spined stickleback (Gasterosteus aculeatus), a non-copulatory, normally oviparous species. Extracted from the parent fish, these embryos hatched and grew to adulthood. Genetic and physiological examination of the parent fish and offspring ruled out development by parthenogenesis and hermaphroditism, therefore implicating internal fertilisation by a male stickleback. This extremely rare phenomenon may have been facilitated in this population by an unusual tendency for females to become egg-bound, and suggests that some major transitions may arise almost spontaneously.

Predictable genome-wide sorting of standing genetic variation during parallel adaptation to basic versus acidic environments in stickleback fish.

Genomic studies of parallel (or convergent) evolution often compare multiple populations diverged into two ecologically different habitats to search for loci repeatedly involved in adaptation. Because the shared ancestor of these populations is generally unavailable, the source of the alleles at adaptation loci, and the direction in which their frequencies were shifted during evolution, remain elusive. To shed light on these issues, we here use multiple populations of threespine stickleback fish adapted to two different types of derived freshwater habitats-basic and acidic lakes on the island of North Uist, Outer Hebrides, Scotland-and the present-day proxy of their marine ancestor. In a first step, we combine genome-wide pooled sequencing and targeted individual-level sequencing to demonstrate that ecological and phenotypic parallelism in basic-acidic divergence is reflected by genomic parallelism in dozens of genome regions. Exploiting data from the ancestor, we next show that the acidic populations, residing in ecologically more extreme derived habitats, have adapted by accumulating alleles rare in the ancestor, whereas the basic populations have retained alleles common in the ancestor. Genomic responses to selection are thus predictable from the ecological difference of each derived habitat type from the ancestral one. This asymmetric sorting of standing genetic variation at loci important to basic-acidic divergence has further resulted in more numerous selective sweeps in the acidic populations. Finally, our data suggest that the maintenance in marine fish of standing variation important to adaptive basic-acidic differentiation does not require extensive hybridization between the marine and freshwater populations. Overall, our study reveals striking genome-wide determinism in both the loci involved in parallel divergence, and in the direction in which alleles at these loci have been selected.

DNA fragility in the parallel evolution of pelvic reduction in stickleback fish.

Evolution generates a remarkable breadth of living forms, but many traits evolve repeatedly, by mechanisms that are still poorly understood. A classic example of repeated evolution is the loss of pelvic hindfins in stickleback fish (Gasterosteus aculeatus). Repeated pelvic loss maps to recurrent deletions of a pelvic enhancer of the Pitx1 gene. Here, we identify molecular features contributing to these recurrent deletions. Pitx1 enhancer sequences form alternative DNA structures in vitro and increase double-strand breaks and deletions in vivo. Enhancer mutability depends on DNA replication direction and is caused by TG-dinucleotide repeats. Modeling shows that elevated mutation rates can influence evolution under demographic conditions relevant for sticklebacks and humans. DNA fragility may thus help explain why the same loci are often used repeatedly during parallel adaptive evolution.

Otolith development in wild populations of stickleback: Jones & Hynes method does not apply to most populations.

This paper critiques Jones & Hynes (1950) findings by analysing sequential samples of otoliths from three wild populations of Gasterosteus aculeatus from North Uist, Scotland and Nottingham, England. Contrary to Jones & Hynes (1950), but coincident with the finding of later researchers, our results showed that no true translucent ring formed in the otolith of G. aculeatus during their first summer. The first translucent ring was probably starting to be formed by the end of summer and was completed by the end of their first winter. There was no second opaque ring in the otoliths of G. aculeatus before they passed their first winter. The second opaque ring was just starting to appear by early April of the second year in the southern population i.e. Nottingham, but later, by May, in the northern populations i.e. North Uist. Formation of the opaque ring in G. aculeatus mostly occurs in spring and summer, with younger fish starting earlier. In contrast, the formation of translucent rings is mostly during autumn and winter, but can be more widespread through the year, possibly as a result of reproductive investment.

A genetics-based approach confirms immune associations with life history across multiple populations of an aquatic vertebrate (Gasterosteus aculeatus).

Understanding how wild immune variation covaries with other traits can reveal how costs and trade-offs shape immune evolution in the wild. Divergent life history strategies may increase or alleviate immune costs, helping shape immune variation in a consistent, testable way. Contrasting hypotheses suggest that shorter life histories may alleviate costs by offsetting them against increased mortality, or increase the effect of costs if immune responses are traded off against development or reproduction. We investigated the evolutionary relationship between life history and immune responses within an island radiation of three-spined stickleback, with discrete populations of varying life histories and parasitism. We sampled two short-lived, two long-lived and an anadromous population using qPCR to quantify current immune profile and RAD-seq data to study the distribution of immune variants within our assay genes and across the genome. Short-lived populations exhibited significantly increased expression of all assay genes, which was accompanied by a strong association with population-level variation in local alleles and divergence in a gene that may be involved in complement pathways. In addition, divergence around the eda gene in anadromous fish is likely associated with increased inflammation. A wider analysis of 15 populations across the island revealed that immune genes across the genome show evidence of having diverged alongside life history strategies. Parasitism and reproductive investment were also important sources of variation for expression, highlighting the caution required when assaying immune responses in the wild. These results provide strong, gene-based support for current hypotheses linking life history and immune variation across multiple populations of a vertebrate model.

Response to Comment on "Precipitation drives global variation in natural selection".

The comment by Myers-Smith and Myers focuses on three main points: (i) the lack of a mechanistic explanation for climate-selection relationships, (ii) the appropriateness of the climate data used in our analysis, and (iii) our focus on estimating climate-selection relationships across (rather than within) taxonomic groups. We address these critiques in our response.

What Are the Environmental Determinants of Phenotypic Selection? A Meta-analysis of Experimental Studies.

Although many selection estimates have been published, the environmental factors that cause selection to vary in space and time have rarely been identified. One way to identify these factors is by experimentally manipulating the environment and measuring selection in each treatment. We compiled and analyzed selection estimates from experimental studies. First, we tested whether the effect of manipulating the environment on selection gradients depends on taxon, trait type, or fitness component. We found that the effect of manipulating the environment was larger when selection was measured on life-history traits or via survival. Second, we tested two predictions about the environmental factors that cause variation in selection. We found support for the prediction that variation in selection is more likely to be caused by environmental factors that have a large effect on mean fitness but not for the prediction that variation is more likely to be caused by biotic factors. Third, we compared selection gradients from experimental and observational studies. We found that selection varied more among treatments in experimental studies than among spatial and temporal replicates in observational studies, suggesting that experimental studies can detect relationships between environmental factors and selection that would not be apparent in observational studies.

Precipitation drives global variation in natural selection.

Climate change has the potential to affect the ecology and evolution of every species on Earth. Although the ecological consequences of climate change are increasingly well documented, the effects of climate on the key evolutionary process driving adaptation-natural selection-are largely unknown. We report that aspects of precipitation and potential evapotranspiration, along with the North Atlantic Oscillation, predicted variation in selection across plant and animal populations throughout many terrestrial biomes, whereas temperature explained little variation. By showing that selection was influenced by climate variation, our results indicate that climate change may cause widespread alterations in selection regimes, potentially shifting evolutionary trajectories at a global scale.

Eda haplotypes in three-spined stickleback are associated with variation in immune gene expression.

Haplotypes underlying local adaptation and speciation are predicted to have numerous phenotypic effects, but few genes involved have been identified, with much work to date concentrating on visible, morphological, phenotypes. The link between genes controlling these adaptive morphological phenotypes and the immune system has seldom been investigated, even though changes in the immune system could have profound adaptive consequences. The Eda gene in three-spined stickleback is one of the best studied major adaptation genes; it directly controls bony plate architecture and has been associated with additional aspects of adaptation to freshwater. Here, we exposed F2 hybrids, used to separate Eda genotype from genetic background, to contrasting conditions in semi-natural enclosures. We demonstrate an association between the Eda haplotype block and the expression pattern of key immune system genes. Furthermore, low plated fish grew less and experienced higher burdens of a common ectoparasite with fitness consequences. Little is currently known about the role of the immune system in facilitating adaptation to novel environments, but this study provides an indication of its potential importance.