Scale-dependent effects of host patch traits on species composition in a stickleback parasite metacommunity.

A core goal of ecology is to understand the abiotic and biotic variables that regulate species distributions and community composition. A major obstacle is that the rules governing species distributions can change with spatial scale. Here, we illustrate this point using data from a spatially nested metacommunity of parasites infecting a metapopulation of threespine stickleback fish from 34 lakes on Vancouver Island, British Columbia. Like most parasite metacommunities, the composition of stickleback parasites differs among host individuals within each host population, and differs between host populations. The distribution of each parasite taxon depends, to varying degrees, on individual host traits (e.g., mass, diet) and on host-population characteristics (e.g., lake size, mean host mass, mean diet). However, in most cases in this data set, a given parasite was regulated by different factors at the host-individual and host-population scales, leading to scale-dependent patterns of parasite-species co-occurrence.

Resistance and tolerance: A hierarchical framework to compare individual versus family-level host contributions in an experimental amphibian-trematode system.

Hosts have two general strategies for mitigating the fitness costs of parasite exposure and infection: resistance and tolerance. The resistance-tolerance framework has been well developed in plant systems, but only recently has it been applied to animal-parasite interactions. However, difficulties associated with estimating fitness, controlling parasite exposure, and quantifying parasite burden have limited application of this framework to animal systems. Here, we used an experimental approach to quantify the relative influence of variation among host individuals and genetic families in determining resistance and tolerance within an amphibian-trematode system. Importantly, we used multiple, alternative metrics to assess each strategy, and employed a Bayesian analytical framework to compare among responses while incorporating uncertainty. Relative to unexposed hosts, exposure to the pathogenic trematode (Ribeiroia ondatrae) reduced the survival and growth of California newts (Taricha torosa) (survival: 93% vs. 74%; growth: 0.29 vs. -0.5 vs mm day -1). Similarly, parasite infection success (the inverse of resistance) ranged from 8% to 100%. Yet despite this broad variation in host resistance and tolerance among individual newts, we found no evidence for transmissable, among-family variation in any of the resistance or tolerance metrics. This suggests that opportunities for evolution of these traits is limited, likely requiring significant increases in mutation, gene flow, or environmental heterogeneity. Our study provides a quantitative framework for evaluating the importance of alternative metrics of resistance and tolerance across multiple time points in the study of host-parasite interactions in animal systems.

Using multi-response models to investigate pathogen coinfections across scales: insights from emerging diseases of amphibians.

Associations among parasites affect many aspects of host-parasite dynamics, but a lack of analytical tools has limited investigations of parasite correlations in observational data that are often nested across spatial and biological scales.Here we illustrate how hierarchical, multiresponse modeling can characterize parasite associations by allowing for hierarchical structuring, offering estimates of uncertainty, and incorporating correlational model structures. After introducing the general approach, we apply this framework to investigate coinfections among four amphibian parasites (the trematodes Ribeiroia ondatrae and Echinostoma spp., the chytrid fungus Batrachochytrium dendrobatidis, and ranaviruses) and among >2000 individual hosts, 90 study sites, and five amphibian host species.Ninety-two percent of sites and 80% of hosts supported two or more pathogen species. Our results revealed strong correlations between parasite pairs that varied by scale (from among hosts to among sites) and classification (microparasite versus macroparasite), but were broadly consistent across taxonomically diverse host species. At the host-scale, infection by the trematode R. ondatrae correlated positively with the microparasites, B. dendrobatidis and ranavirus, which were themselves positively associated. However, infection by a second trematode (Echinostoma spp.) correlated negatively with B. dendrobatidis and ranavirus, both at the host- and site-level scales, highlighting the importance of differential relationships between micro- and macroparasites.Given the extensive number of coinfecting symbiont combinations inherent to natural systems, particularly across multiple host species, multiresponse modeling of cross-sectional field data offers a valuable tool to identify a tractable number of hypothesized interactions for experimental testing while accounting for uncertainty and potential sources of co-exposure. For amphibians specifically, the high frequency of co-occurrence and coinfection among these pathogens - each of which is known to impair host fitness or survival - highlights the urgency of understanding parasite associations for conservation and disease management.

Gene expression stasis and plasticity following migration into a foreign environment.

Selection against migrants is key to maintaining genetic differences between populations linked by dispersal. However, migrants may mitigate fitness costs by proactively choosing among available habitats, or by phenotypic plasticity. We previously reported that a reciprocal transplant of lake and stream stickleback (Gasterosteus aculeatus) found little support for divergent selection. Here, we revisit that experiment to test whether phenotypic plasticity in gene expression may have helped migrants adjust to unfamiliar habitats. We measured gene expression profiles in stickleback via TagSeq and tested whether migrants between lake and stream habitats exhibited a plastic response to their new environment that allowed them to converge on the expression profile of adapted natives. We report extensive gene expression differences between genetically divergent lake and stream stickleback, despite gene flow. But for many genes, expression was highly plastic. Fish transplanted into the adjoining habitat partially converged on the expression profile typical of natives from their new habitat. This suggests that expression plasticity may soften the impact of migration. Nonetheless, lake and stream fish differed in survival rates and parasite infection rates in our study, implying that expression plasticity is not fast or extensive enough to fully homogenize fish performance.

Frequency dependence limits divergent evolution by favouring rare immigrants over residents.

Two distinct forms of natural selection promote adaptive biological diversity. Divergent selection occurs when different environments favour different phenotypes, leading to increased differences between populations. Negative frequency-dependent selection occurs when rare variants within a population are favoured over common ones, increasing diversity within populations. These two diversifying forces promote genetic variation at different spatial scales, and may act in opposition, but their relative effects remain unclear because they are rarely measured concurrently. Here we show that negative frequency-dependent selection within populations can favor rare immigrants over locally adapted residents. We reciprocally transplanted lake and stream ecotypes of threespine stickleback into lake and stream habitats, while manipulating the relative abundance of residents versus immigrants. We found negative frequency-dependence: survival was highest for the locally rare ecotype, rather than natives. Also, individuals with locally rare major histocompatibility complex (MHC) class IIb genotypes were infected by fewer parasites. This negative frequency-dependent selection will tend to favour rare immigrants over common residents, amplifying the effect of migration and undermining the efficacy of divergent natural selection to drive population differences. The only signal of divergent selection was a tendency for foreign fish to have higher parasite loads than residents, after controlling for MHC genotype rarity. Frequency-dependent ecological interactions have long been thought to promote speciation. Our results suggest a more nuanced view in which negative frequency dependence alters the fate of migrants to promote or constrain evolutionary divergence between populations.

Natural selection on MHC IIß in parapatric lake and stream stickleback: Balancing, divergent, both or neither?

Major histocompatibility complex (MHC) genes encode proteins that play a central role in vertebrates' adaptive immunity to parasites. MHC loci are among the most polymorphic in vertebrates' genomes, inspiring many studies to identify evolutionary processes driving MHC polymorphism within populations and divergence between populations. Leading hypotheses include balancing selection favouring rare alleles within populations, and spatially divergent selection. These hypotheses do not always produce diagnosably distinct predictions, causing many studies of MHC to yield inconsistent or ambiguous results. We suggest a novel strategy to distinguish balancing vs. divergent selection on MHC, taking advantage of natural admixture between parapatric populations. With divergent selection, individuals with immigrant alleles will be more infected and less fit because they are susceptible to novel parasites in their new habitat. With balancing selection, individuals with locally rare immigrant alleles will be more fit (less infected). We tested these contrasting predictions using three-spine stickleback from three replicate pairs of parapatric lake and stream habitats. We found numerous positive and negative associations between particular MHC IIß alleles and particular parasite taxa. A few allele-parasite comparisons supported balancing selection, and others supported divergent selection between habitats. But, there was no overall tendency for fish with immigrant MHC alleles to be more or less heavily infected. Instead, locally rare MHC alleles (not necessarily immigrants) were associated with heavier infections. Our results illustrate the complex relationship between MHC IIß allelic variation and spatially varying multispecies parasite communities: different hypotheses may be concurrently true for different allele-parasite combinations.

Partitioning the effects of isolation by distance, environment, and physical barriers on genomic divergence between parapatric threespine stickleback.

Genetic divergence between populations is shaped by a combination of drift, migration, and selection, yielding patterns of isolation-by-distance (IBD) and isolation-by-environment (IBE). Unfortunately, IBD and IBE may be confounded when comparing divergence across habitat boundaries. For instance, parapatric lake and stream threespine stickleback (Gasterosteus aculeatus) may have diverged due to selection against migrants (IBE), or mere spatial separation (IBD). To quantitatively partition the strength of IBE and IBD, we used recently developed population genetic software (BEDASSLE) to analyze partial genomic data from three lake-stream clines on Vancouver Island. We find support for IBD within each of three outlet streams (unlike prior studies of lake-stream stickleback). In addition, we find evidence for IBE (controlling for geographic distance): the genetic effect of habitat is equivalent to geographic separation of ∼1.9 km of IBD. Remarkably, of our three lake-stream pairs, IBE is strongest where migration between habitats is easiest. Such microgeographic genetic divergence would require exceptionally strong divergent selection, which multiple experiments have failed to detect. Instead, we propose that nonrandom dispersal (e.g., habitat choice) contributes to IBE. Supporting this conclusion, we show that the few migrants between habitats are a nonrandom subset of the phenotype distribution of the source population.

Multilevel Models for the Distribution of Hosts and Symbionts.

Symbiont occurrence is influenced by host occurrence and vice versa, which leads to correlations in host-symbiont distributions at multiple levels. Interactions between co-infecting symbionts within host individuals can cause correlations in the abundance of two symbiont species across individual hosts. Similarly, interactions between symbiont transmission and host population dynamics can drive correlations between symbiont and host abundance across habitat patches. If ignored, these interactions can confound estimated responses of hosts and symbionts to other factors. Here, we present a general hierarchical modeling framework for distributions of hosts and symbionts, estimating correlations in host-symbiont distributions at the among-site, within-site, among-species, and among-individual levels. We present an empirical example from a multi-host multi-parasite system involving amphibians and their micro- and macroparasites. Amphibian hosts and their parasites were correlated at multiple levels of organization. Macroparasites often co-infected individual hosts, but rarely co-infected with the amphibian chytrid fungus. Such correlations may result from interactions among parasites and hosts, joint responses to environmental factors, or sampling bias. Joint host-symbiont models account for environmental constraints and species interactions while partitioning variance and dependence in abundance at multiple levels. This framework can be adapted to a wide variety of study systems and sampling designs.

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.

Stepwise threshold clustering: a new method for genotyping MHC loci using next-generation sequencing technology.

Genes of the vertebrate major histocompatibility complex (MHC) are of great interest to biologists because of their important role in immunity and disease, and their extremely high levels of genetic diversity. Next generation sequencing (NGS) technologies are quickly becoming the method of choice for high-throughput genotyping of multi-locus templates like MHC in non-model organisms. Previous approaches to genotyping MHC genes using NGS technologies suffer from two problems:1) a "gray zone" where low frequency alleles and high frequency artifacts can be difficult to disentangle and 2) a similar sequence problem, where very similar alleles can be difficult to distinguish as two distinct alleles. Here were present a new method for genotyping MHC loci--Stepwise Threshold Clustering (STC)--that addresses these problems by taking full advantage of the increase in sequence data provided by NGS technologies. Unlike previous approaches for genotyping MHC with NGS data that attempt to classify individual sequences as alleles or artifacts, STC uses a quasi-Dirichlet clustering algorithm to cluster similar sequences at increasing levels of sequence similarity. By applying frequency and similarity based criteria to clusters rather than individual sequences, STC is able to successfully identify clusters of sequences that correspond to individual or similar alleles present in the genomes of individual samples. Furthermore, STC does not require duplicate runs of all samples, increasing the number of samples that can be genotyped in a given project. We show how the STC method works using a single sample library. We then apply STC to 295 threespine stickleback (Gasterosteus aculeatus) samples from four populations and show that neighboring populations differ significantly in MHC allele pools. We show that STC is a reliable, accurate, efficient, and flexible method for genotyping MHC that will be of use to biologists interested in a variety of downstream applications.

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.

Contrasting patterns of phenotype-dependent parasitism within and among populations of threespine stickleback.

Variation in infection rate arises from variation in host exposure and resistance to parasites both within and among populations. All things being equal, phenotypes that increase exposure risk should covary positively with infection among individuals. It might therefore be expected that populations with mean phenotypes that increase exposure might also have higher rates of infection. However, such positive covariance between exposure and infection at the population level might be undermined by other factors such as geographic variation in parasite abundance or host resistance, negating or reversing in between-population comparisons. We studied rates of infection of two parasites among 18 populations of threespine stickleback (Gasterosteus aculeatus). As predicted, within populations, trophic morphology covaries with infection of two trophically transmitted parasites: individuals with benthic (or limnetic) phenotypes were more likely to be infected with a benthic (or limnetic) parasite. However, across populations, the relationship between morphology and infection rate was absent (limnetic parasite) or reversed (benthic parasite). Our results confirm the importance of phenotype-dependent exposure, but stress different factors or processes, such as the evolution of reduced susceptibility, might shape variation in infection at larger spatial scales.

Does intraspecific size variation in a predator affect its diet diversity and top-down control of prey?

It has long been known that intraspecific variation impacts evolutionary processes, but only recently have its potential ecological effects received much attention. Theoretical models predict that genetic or phenotypic variance within species can alter interspecific interactions, and experiments have shown that genotypic diversity in clonal species can impact a wide range of ecological processes. To extend these studies to quantitative trait variation within populations, we experimentally manipulated the variance in body size of threespine stickleback in enclosures in a natural lake environment. We found that body size of stickleback in the lake is correlated with prey size and (to a lesser extent) composition, and that stickleback can exert top-down control on their benthic prey in enclosures. However, a six-fold contrast in body size variance had no effect on the degree of diet variation among individuals, or on the abundance or composition of benthic or pelagic prey. Interestingly, post-hoc analyses revealed suggestive correlations between the degree of diet variation and the strength of top-down control by stickleback. Our negative results indicate that, unless the correlation between morphology and diet is very strong, ecological variation among individuals may be largely decoupled from morphological variance. Consequently we should be cautious in our interpretation both of theoretical models that assume perfect correlations between morphology and diet, and of empirical studies that use morphological variation as a proxy for resource use diversity.

Foraging trait (co)variances in stickleback evolve deterministically and do not predict trajectories of adaptive diversification.

How does natural selection shape the structure of variance and covariance among multiple traits, and how do (co)variances influence trajectories of adaptive diversification? We investigate these pivotal but open questions by comparing phenotypic (co)variances among multiple morphological traits across 18 derived lake-dwelling populations of threespine stickleback, and their marine ancestor. Divergence in (co)variance structure among populations is striking and primarily attributable to shifts in the variance of a single key foraging trait (gill raker length). We then relate this divergence to an ecological selection proxy, to population divergence in trait means, and to the magnitude of sexual dimorphism within populations. This allows us to infer that evolution in (co)variances is linked to variation among habitats in the strength of resource-mediated disruptive selection. We further find that adaptive diversification in trait means among populations has primarily involved shifts in gill raker length. The direction of evolutionary trajectories is unrelated to the major axes of ancestral trait (co)variance. Our study demonstrates that natural selection drives both means and (co)variances deterministically in stickleback, and strongly challenges the view that the (co)variance structure biases the direction of adaptive diversification predictably even over moderate time spans.

Ecological release from interspecific competition leads to decoupled changes in population and individual niche width.

A species's niche width reflects a balance between the diversifying effects of intraspecific competition and the constraining effects of interspecific competition. This balance shifts when a species from a competitive environment invades a depauperate habitat where interspecific competition is reduced. The resulting ecological release permits population niche expansion, via increased individual niche widths and/or increased among-individual variation. We report an experimental test of the theory of ecological release in three-spine stickleback (Gasterosteus aculeatus). We factorially manipulated the presence or absence of two interspecific competitors: juvenile cut-throat trout (Oncorhynchus clarki) and prickly sculpin (Cottus asper). Consistent with the classic niche variation hypothesis, release from trout competition increased stickleback population niche width via increased among-individual variation, while individual niche widths remained unchanged. In contrast, release from sculpin competition had no effect on population niche width, because increased individual niche widths were offset by decreased between-individual variation. Our results confirm that ecological release from interspecific competition can lead to increases in niche width, and that these changes can occur on behavioural time scales. Importantly, we find that changes in population niche width are decoupled from changes in the niche widths of individuals within the population.

Phenotype-dependent native habitat preference facilitates divergence between parapatric lake and stream stickleback.

Adaptive divergence between adjoining populations reflects a balance between the diversifying effect of divergent selection and the potentially homogenizing effect of gene flow. In most models of migration-selection balance, gene flow is assumed to reflect individuals' inherent capacity to disperse, without regard to the match between individuals' phenotypes and the available habitats. However, habitat preferences can reduce dispersal between contrasting habitats, thereby alleviating migration load and facilitating adaptive divergence. We tested whether habitat preferences contribute to adaptive divergence in a classic example of migration-selection balance: parapatric lake and stream populations of three-spine stickleback (Gasterosteus aculeatus). Using a mark-transplant-recapture experiment on morphologically divergent parapatric populations, we showed that 90% of lake and stream stickleback returned to their native habitat, reducing migration between habitats by 76%. Furthermore, we found that dispersal into a nonnative habitat was phenotype dependent. Stream fish moving into the lake were morphologically more lake-like than those returning to the stream (and the converse for lake fish entering the stream). The strong native habitat preference documented here increases the extent of adaptive divergence between populations two- to fivefold relative to expectations with random movement. These results illustrate the potential importance of adaptive habitat choice in driving parapatric divergence.