Individuals' diet diversity influences gut microbial diversity in two freshwater fish (threespine stickleback and Eurasian perch)

Vertebrates' diets profoundly influence the composition of symbiotic gut microbial communities. Studies documenting diet‐microbiota associations typically focus on univariate or categorical diet variables. However, in nature individuals often consume diverse combinations of foods. If diet components act independently, each providing distinct microbial colonists or nutrients, we expect a positive relationship between diet diversity and microbial diversity. We tested this prediction within each of two fish species (stickleback and perch), in which individuals vary in their propensity to eat littoral or pelagic invertebrates or mixtures of both prey. Unexpectedly, in most cases individuals with more generalised diets had less diverse microbiota than dietary specialists, in both natural and laboratory populations. This negative association between diet diversity and microbial diversity was small but significant, and most apparent after accounting for complex interactions between sex, size and diet. Our results suggest that multiple diet components can interact non‐additively to influence gut microbial diversity.     

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.     

Admixture mapping of male nuptial colour and body shape in a recently formed hybrid population of threespine stickleback

Despite recent progress, we still know relatively little about the genetic architecture that underlies adaptation to divergent environments. Determining whether the genetic architecture of phenotypic adaptation follows any predictable patterns requires data from a wide variety of species. However, in many organisms, genetic studies are hindered by the inability to perform genetic crosses in the laboratory or by long generation times. Admixture mapping is an approach that circumvents these issues by taking advantage of hybridization that occurs between populations or species in the wild. Here, we demonstrate the utility of admixture mapping in a naturally occurring hybrid population of threespine sticklebacks (Gasterosteus aculeatus) from Enos Lake, British Columbia. Until recently, this lake contained two species of sticklebacks adapted to divergent habitats within the lake. This benthic–limnetic species pair diverged in a number of phenotypes, including male nuptial coloration and body shape, which were previously shown to contribute to reproductive isolation between them. However, recent ecological disturbance has contributed to extensive hybridization between the species, and there is now a single, admixed population within Enos Lake. We collected over 500 males from Enos Lake and found that most had intermediate nuptial colour and body shape. By genotyping males with nuptial colour at the two extremes of the phenotypic distribution, we identified seven genomic regions on three chromosomes associated with divergence in male nuptial colour. These genomic regions are also associated with variation in body shape, suggesting that tight linkage and/or pleiotropy facilitated adaptation to divergent environments in benthic–limnetic species pairs.     

Connectivity in a pond system influences migration and genetic structure in threespine stickleback

Neutral genetic structure of natural populations is primarily influenced by migration (the movement of individuals and, subsequently, their genes) and drift (the statistical chance of losing genetic diversity over time). Migration between populations is influenced by several factors, including individual behavior, physical barriers, and environmental heterogeneity among populations. However, drift is expected to be stronger in populations with low immigration rate and small effective population size. With the technological advancement in geological information systems and spatial analysis tools, landscape genetics now allows the development of realistic migration models and increased insight to important processes influencing diversity of natural populations. In this study, we investigated the relationship between landscape connectivity and genetic distance of threespine stickleback (Gasterosteus aculeatus) inhabiting a pond complex in Belgjarskógur, Northeast Iceland. We used two landscape genetic approaches (i.e., least-cost-path and isolation-by-resistance) and asked whether gene flow, as measured by genetic distance, was more strongly associated with Euclidean distance (isolation-by-distance) or with landscape connectivity provided by areas prone to flooding (as indicated by Carex sp. cover)? We found substantial genetic structure across the study area, with pairwise genetic distances among populations (D_PS) ranging from 0.118 to 0.488. Genetic distances among populations were more strongly correlated with least-cost-path and isolation-by-resistance than with Euclidean distance, whereas the relative contribution of isolation-by-resistance and Euclidian distance could not be disentangled. These results indicate that migration among stickleback populations occurs via periodically flooded areas. Overall, this study highlights the importance of transient landscape elements influencing migration and genetic structure of populations at small spatial scales.     

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.     

Among‐lake reciprocal transplants induce convergent expression of immune genes in threespine stickleback

Geographic variation in parasite communities can drive evolutionary divergence in host immune genes. However, biotic and abiotic environmental variation can also induce plastic differences in immune function among populations. At present, there is little information concerning the relative magnitudes of heritable vs. induced immune divergence in natural populations. We examined immune gene expression profiles of threespine stickleback (Gasterosteus aculeatus) from six lakes on Vancouver Island, British Columbia. Parasite community composition differs between lake types (large or small, containing limnetic‐ or benthic‐like stickleback) and between watersheds. We observed corresponding differences in immune gene expression profiles among wild‐caught stickleback, using a set of seven immune genes representing distinct branches of the immune system. To evaluate the role of environmental effects on this differentiation, we experimentally transplanted wild‐caught fish into cages in their native lake, or into a nearby foreign lake. Transplanted individuals' immune gene expression converged on patterns typical of their destination lake, deviating from their native expression profile. Transplant individuals' source population had a much smaller effect, suggesting relatively weak genetic underpinning of population differences in immunity, as viewed through gene expression. This strong environmental regulation of immune gene expression provides a counterpoint to the large emerging literature documenting microevolution and genetic diversification of immune function. Our findings illustrate the value of studying immunity in natural environmental settings where the immune system has evolved and actively functions.     

Quantifying the constraining influence of gene flow on adaptive divergence in the lake-stream threespine stickleback system

The constraining effect of gene flow on adaptive divergence is often inferred but rarely quantified. We illustrate ways of doing so using stream populations of threespine stickleback (Gasterosteus aculeatus) that experience different levels of gene flow from a parapatric lake population. In the Misty Lake watershed (British Columbia, Canada), the inlet stream population is morphologically divergent from the lake population, and presumably experiences little gene flow from the lake. The outlet stream population, however, shows an intermediate phenotype and may experience more gene flow from the lake. We first used microsatellite data to demonstrate that gene flow from the lake is low into the inlet but high into the outlet, and that gene flow from the lake remains relatively constant with distance along the outlet. We next combined gene flow data with morphological and habitat data to quantify the effect of gene flow on morphological divergence. In one approach, we assumed that inlet stickleback manifest well-adapted phenotypic trait values not constrained by gene flow. We then calculated the deviation between the observed and expected phenotypes for a given habitat in the outlet. In a second approach, we parameterized a quantitative genetic model of adaptive divergence. Both approaches suggest a large impact of gene flow, constraining adaptation by 80-86% in the outlet (i.e., only 14-20% of the expected morphological divergence in the absence of gene flow was observed). Such approaches may be useful in other taxa to estimate how important gene flow is in constraining adaptive divergence in nature.     

Evolution and assessment of colour patterns in stream-resident and anadromous male threespine stickleback Gasterosteus aculeatus from three regions

We compared the colour patterns of free swimming, reproductively active male threespine stickleback Gasterosteus aculeatus of the anadromous and stream ecotypes from three geographically distinct regions. Consistent with the hypothesis of environmentally mediated selection, our results indicate ecologically replicated differences in G. aculeatus coloration between anadromous and stream-resident populations, and that G. aculeatus probably have the visual acuity to discriminate colour pattern differences between anadromous and stream-resident fish.     

Strong and parallel salinity‐induced phenotypic plasticity in one generation of threespine stickleback

Phenotypic plasticity is a major factor contributing to variation of organisms in nature, yet its evolutionary significance is insufficiently understood. One example system where plasticity might have played an important role in an adaptive radiation is the threespine stickleback (Gasterosteus aculeatus), a fish that has diversified after invading freshwater lakes repeatedly from the marine habitat. The parallel phenotypic changes that occurred in this radiation were extremely rapid. This study evaluates phenotypic plasticity in stickleback body shape in response to salinity in fish stemming from a wild freshwater population. Using a split‐clutch design, we detected surprisingly large phenotypically plastic changes in body shape after one generation. Fish raised in salt water developed shallower bodies and longer jaws, and these changes were consistent and parallel across families. Although this work highlights the effect of phenotypic plasticity, we also find indications that constraints may play a role in biasing the direction of possible phenotypic change. The slopes of the allometric relationship of individual linear traits did not change across treatments, indicating that plastic change does not affect the covariation of traits with overall size. We conclude that stickleback have a large capacity for plastic phenotypic change in response to salinity and that plasticity and evolutionary constraints have likely contributed to the phenotypic diversification of these fish.     

An experimental test of the effect of predation upon behaviour and trait correlations in the threespine stickleback

Benthic and limnetic threespine stickleback (Gasterosteus aculeatus) are a classic example of ecological speciation. Behavioural and armour divergence between these species has been predicted to be the result of divergent selection driven in part by differential predation from cutthroat trout (Oncorhynchus clarki). To experimentally test this prediction, we reared split families of benthic–limnetic hybrids in the presence or absence of trout predation. Our results show that the presence of trout had little effect upon stickleback behaviour. We then compared performance in behavioural assays among stickleback that varied in armour to test if armour morphology correlates with behaviour. Our measurements revealed trait correlations between several behaviours and components of armour morphology. Trout predation did not result in an increased correlation between traits, therefore differential trout predation between benthics and limnetics is unlikely to be the cause of these correlations. The presence of trait correlations in advanced generation hybrids suggests that pleiotropy or linkage between genes underlying behaviour and armour morphology may be greater than previously appreciated.     

Nonlinear effects of temperature on body form and developmental canalization in the threespine stickleback

Theoretical models predict that nonlinear environmental effects on the phenotype also affect developmental canalization, which in turn can influence the tempo and course of organismal evolution. Here, we used an oceanic population of threespine stickleback (Gasterosteus aculeatus) to investigate temperature‐induced phenotypic plasticity of body size and shape using a paternal half‐sibling, split‐clutch experimental design and rearing offspring under three different temperature regimes (13, 17 and 21 °C). Body size and shape of 466 stickleback individuals were assessed by a set of 53 landmarks and analysed using geometric morphometric methods. At approximately 100 days, individuals differed significantly in both size and shape across the temperature groups. However, the temperature‐induced differences between 13 and 17 °C (mainly comprising relative head and eye size) deviated considerably from those between 17 and 21 °C (involving the relative size of the ectocoracoid, the operculum and the ventral process of the pelvic girdle). Body size was largest at 17 °C. For both size and shape, phenotypic variance was significantly smaller at 17 °C than at 13 and 21 °C, indicating that development is most stable at the intermediate temperature matching the conditions encountered in the wild. Higher additive genetic variance at 13 and 21 °C indicates that the plastic response to temperature had a heritable basis. Understanding nonlinear effects of temperature on development and the underlying genetics are important for modelling evolution and for predicting outcomes of global warming, which can lead not only to shifts in average morphology but also to destabilization of development.     

Computational identification of miRNAs,their targets and functions in three‐spined stickleback (Gasterosteus aculeatus)

An intriguing question in biology is how the evolution of gene regulation is shaped by natural selection in natural populations. Among the many known regulatory mechanisms, regulation of gene expression by microRNAs (miRNAs) is of critical importance. However, our understanding of their evolution in natural populations is limited. Studying the role of miRNAs in three‐spined stickleback, an important natural model for speciation research, may provide new insights into adaptive polymorphisms. However, lack of annotation of miRNA genes in its genome is a bottleneck. To fill this research gap, we used the genome of three‐spined stickleback to predict miRNAs and their targets. We predicted 1486 mature miRNAs using the homology‐based miRNA prediction approach. We then performed functional annotation and enrichment analysis of these targets, which identified over‐represented motifs. Further, a database resource (GAmiRdb) has been developed for dynamically searching miRNAs and their targets exclusively in three‐spined stickleback. Finally, the database was used in two case studies focusing on freshwater adaptation in natural populations. In the first study, we found 44 genomic regions overlapping with predicted miRNA targets. In the second study, we identified two SNPs altering the MRE seed site of sperm‐specific glyceraldehyde‐3‐phosphate gene. These findings highlight the importance of the GAmiRdb knowledge base in understanding adaptive evolution.     

Major Histocompatibility Complex class IIb polymorphism influences gut microbiota composition and diversity

Animals harbour diverse communities of symbiotic bacteria, which differ dramatically among host individuals. This heterogeneity poses an immunological challenge: distinguishing between mutualistic and pathogenic members of diverse and host‐specific microbial communities. We propose that Major Histocompatibility class II (MHC) genotypes contribute to recognition and regulation of gut microbes, and thus, MHC polymorphism contributes to microbial variation among hosts. Here, we show that MHC IIb polymorphism is associated with among‐individual variation in gut microbiota within a single wild vertebrate population of a small fish, the threespine stickleback. We sampled stickleback from Cedar Lake, on Vancouver Island, and used next‐generation sequencing to genotype the sticklebacks’ gut microbiota (16S sequencing) and their MHC class IIb exon 2 sequences. The presence of certain MHC motifs was associated with altered relative abundance (increase or decrease) of some microbial Families. The effect sizes are modest and entail a minority of microbial taxa, but these results represent the first indication that MHC genotype may affect gut microbiota composition in natural populations (MHC‐microbe associations have also been found in a few studies of lab mice). Surprisingly, these MHC effects were frequently sex‐dependent. Finally, hosts with more diverse MHC motifs had less diverse gut microbiota. One implication is that MHC might influence the efficacy of therapeutic strategies to treat dysbiosis‐associated disease, including the outcome of microbial transplants between healthy and diseased patients. We also speculate that macroparasite‐driven selection on MHC has the potential to indirectly alter the host gut microbiota, and vice versa.     

Microhabitat contributes to microgeographic divergence in threespine stickleback

Since the New Synthesis, most migration-selection balance theory has predicted that there should be negligible differentiation over small spatial scales (relative to dispersal), because gene flow should erode any effect of divergent selection. Nevertheless, there are classic examples of microgeographic divergence, which theory suggests can arise under specific conditions: exceptionally strong selection, phenotypic plasticity in philopatric individuals, or nonrandom dispersal. Here, we present evidence of microgeographic morphological variation within lake and stream populations of threespine stickleback (Gasterosteus aculeatus). It seems reasonable to assume that a given lake or stream population of fish is well-mixed. However, we found this assumption to be untenable. We examined trap-to-trap variation in 34 morphological traits measured on stickleback from 16 lakes and 16 streams. Most traits varied appreciably among traps within populations. Both between-trap distance and microhabitat characteristics such as depth and substrate explained some of the within-population morphological variance. Microhabitat was also associated with genotype at particular loci but there was no genetic isolation by distance, implying that heritable habitat preferences may contribute to microgeographic variation. Our study adds to growing evidence that microgeographic divergence can occur across small spatial scales within individuals’ daily dispersal neighborhood where gene flow is expected to be strong.     

Lateral plate evolution in the threespine stickleback: getting nowhere fast

Gasterosteus aculeatus is a small Holarctic fish with marine, anadromous, and freshwater populations. Marine and anadromous populations apparently have changed little in the past 10 million years and exhibit limited geographical variation. In contrast, freshwater isolates have been founded repeatedly by marine and anadromous populations, and post-glacial isolates have undergone extraordinary adaptive radiation. Stickleback traits that have diversified during post-glacial radiation, including the lateral plates (LP), can evolve substantially within decades after colonization of fresh water or when the environment (particularly predation regime) changes. Although highly divergent freshwater isolates of G. aculeatus have existed for at least 10 million years, they have rarely experienced sustained evolutionary divergence leading to formation of widespread, phenotypically distinct species. The paradox of rapid LP evolution without sustained divergence has resulted from selective extinction of highly divergent populations, because they are specialized for conditions in small, isolated habitats that tend to dry up within limited periods. Biological species of G. aculeatus may also evolve within decades, and are also prone to extinction because they are endemic to and specialized for small, ephemeral habitats. The high rate of evolution observed in contemporary threespine stickleback populations may not be unique to this species complex and has important implications for use of post-glacial populations in comparative studies, speciation rate, and discrimination of sympatric and allopatric speciation.     

Red coloration and male parental behaviour in the threespine stickleback

In an anadromous population of the threespine stickleback, male red coloration was correlated with agonistic and egg care behaviour of parental males, but not with hatching success. When one male took over a neighbour's territory, the dominant male had more intense red coloration.     

Ecological causes of morphological evolution in the three‐spined stickleback

The central assumption of evolutionary theory is that natural selection drives the adaptation of populations to local environmental conditions, resulting in the evolution of adaptive phenotypes. The three‐spined stickleback (Gasterosteus aculeatus) displays remarkable phenotypic variation, offering an unusually tractable model for understanding the ecological mechanisms underpinning adaptive evolutionary change. Using populations on North Uist, Scotland we investigated the role of predation pressure and calcium limitation on the adaptive evolution of stickleback morphology and behavior. Dissolved calcium was a significant predictor of plate and spine morph, while predator abundance was not. Stickleback latency to emerge from a refuge varied with morph, with populations with highly reduced plates and spines and high predation risk less bold. Our findings support strong directional selection in three‐spined stickleback evolution, driven by multiple selective agents.