Resource diversity promotes among-individual diet variation, but not genomic diversity, in lake stickleback

The genetic basis of divergent melanic pigmentation in benthic and limnetic threespine stickleback

Pigmentation is an excellent trait to examine patterns of evolutionary change because it is often under natural selection. Benthic and limnetic threespine stickleback (Gasterosteus aculeatus) exhibit distinct pigmentation phenotypes, likely an adaptation to occupation of divergent niches. The genetic architecture of pigmentation in vertebrates appears to be complex. Prior QTL mapping of threespine stickleback pigmentation phenotypes has identified several candidate loci. However-relative to other morphological phenotypes (e.g., spines or lateral plates)-the genetic architecture of threespine stickleback pigmentation remains understudied. Here, we performed QTL mapping for two melanic pigmentation traits (melanophore density and lateral barring) using benthic-limnetic F2 crosses. The two traits mapped to different chromosomes, suggesting a distinct genetic basis. The resulting QTLs were additive, but explained a relatively small fraction of the total variance (~6%). QTLs maps differed by F1 family, suggesting variation in genetic architecture or ability to detect loci of small effect. Functional analysis identified enriched pathways for candidate loci. Several of the resulting candidate loci for pigmentation, including three loci in enriched pathways (bco1, sulf1, and tyms) have been previously indicated to affect pigmentation in other vertebrates. These findings add to a growing body of evidence suggesting pigmentation is often polygenic.

Destabilized host-parasite dynamics in newly founded populations

When species disperse into previously unoccupied habitats, new populations encounter unfamiliar species interactions such as altered parasite loads. Theory predicts that newly founded populations should exhibit destabilized eco-evolutionary fluctuations in infection rates and immune traits. However, to understand founder effects biologists typically rely on retrospective studies of range expansions, missing early-generation infection dynamics. To remedy this, we experimentally founded whole-lake populations of threespine stickleback. Infection rates were temporally stable in native source lakes. In contrast, newly founded populations exhibit destabilized host-parasite dynamics: high starting infection rates led to increases in a heritable immune trait (peritoneal fibrosis), suppressing infection rates. The resulting temporal auto-correlation between infection and immunity suggest that newly founded populations can exhibit rapid host-parasite eco-evolutionary dynamics.

Designing eco-evolutionary experiments for restoration projects: Opportunities and constraints revealed during stickleback introductions

Eco-evolutionary experiments are typically conducted in semi-unnatural controlled settings, such as mesocosms; yet inferences about how evolution and ecology interact in the real world would surely benefit from experiments in natural uncontrolled settings. Opportunities for such experiments are rare but do arise in the context of restoration ecology-where different "types" of a given species can be introduced into different "replicate" locations. Designing such experiments requires wrestling with consequential questions. (Q1) Which specific "types" of a focal species should be introduced to the restoration location? (Q2) How many sources of each type should be used-and should they be mixed together? (Q3) Which specific source populations should be used? (Q4) Which type(s) or population(s) should be introduced into which restoration sites? We recently grappled with these questions when designing an eco-evolutionary experiment with threespine stickleback (Gasterosteus aculeatus) introduced into nine small lakes and ponds on the Kenai Peninsula in Alaska that required restoration. After considering the options at length, we decided to use benthic versus limnetic ecotypes (Q1) to create a mixed group of colonists from four source populations of each ecotype (Q2), where ecotypes were identified based on trophic morphology (Q3), and were then introduced into nine restoration lakes scaled by lake size (Q4). We hope that outlining the alternatives and resulting choices will make the rationales clear for future studies leveraging our experiment, while also proving useful for investigators considering similar experiments in the future.

Host ecotype and rearing environment are the main drivers of threespine stickleback gut microbiota diversity in a naturalistic experiment

Host-microbiota interactions play a critical role in the hosts' biology, and thus, it is crucial to elucidate the mechanisms that shape gut microbial communities. We leveraged threespine stickleback fish (Gasterosteus aculeatus) as a model system to investigate the contribution of host and environmental factors to gut microbiota variation. These fish offer a unique opportunity for experiments in naturalistic conditions; we reared benthic and limnetic ecotypes from three different lakes in experimental ponds, allowing us to assess the relative effects of shared environment (pond), geographic origin (lake-of-origin), trophic ecology and genetics (ecotype) and biological sex on gut microbiota α- and β-diversity. Host ecotype had the strongest influence on α-diversity, with benthic fish exhibiting higher diversity than limnetic fish, followed by the rearing environment. β-diversity was primarily shaped by rearing environment, followed by host ecotype, indicating that environmental factors play a crucial role in determining gut microbiota composition. Furthermore, numerous bacterial orders were differentially abundant across ponds, underlining the substantial contribution of environmental factors to gut microbiota variation. Our study illustrates the complex interplay between environmental and host ecological or genetic factors in shaping the stickleback gut microbiota and highlights the value of experiments conducted under naturalistic conditions for understanding gut microbiota dynamics.

Viral community distribution, assembly mechanism, and associated hosts in an industrial park wastewater treatment plant

Viruses manipulate bacterial community composition and impact wastewater treatment efficiency. Some viruses pose threats to the environment and human populations through infection. Improving the efficiency of wastewater treatment and ensuring the health of the effluent and receptor pools requires an understanding of how viral communities assemble and interact with hosts in wastewater treatment plants (WWTPs). We used metagenomic analysis to study the distribution, assembly mechanism, and sensitive hosts for the viral communities in raw water, anaerobic tanks, and returned activated sludge units of a large-scale industrial park WWTP. Uroviricota (53.42% ± 0.14%) and Nucleocytoviricota (26.1% ± 0.19%) were dominant in all units. Viral community composition significantly differed between units, as measured by β diversity (P = 0.005). Compared to raw water, the relative viral abundance decreased by 29.8% in the anaerobic tank but increased by 9.9% in the activated sludge. Viral community assembly in raw water and anaerobic tanks was predominantly driven by deterministic processes (MST <0.5) versus stochastic processes (MST >0.5) in the activated sludge, indicating that differences in diffusion limits may fundamentally alter the assembly mechanisms of viral communities between the solid and liquid-phase environments. Acidobacteria was identified as the sensitive host contributing to viral abundance, exhibiting strong interactions and a mutual dependence (degree = 59). These results demonstrate the occurrence and prevalence of viruses in WWTPs, their different assembly mechanism, and sensitive hosts. These observations require further study of the mechanisms of viral community succession, ecological function, and roles in the successive wastewater treatment units.

Environmental correlates of adaptive diversification in postglacial freshwater fishes

Determining how environmental conditions contribute to divergence among populations and drive speciation is fundamental to resolving mechanisms and understanding outcomes in evolutionary biology. Postglacial freshwater fish species in the Northern Hemisphere are ideal biological systems to explore the effects of environment on diversification in morphology, ecology, and genetics (ecomorph divergences) within lakes. To date, various environmental factors have been implicated in the presence of multiple ecomorphs within particular lakes or regions. However, concerted evidence for generalizable patterns in environmental variables associated with speciation across geographical regions and across species and genera has been lacking. Here, we aimed to identify key biotic and abiotic factors associated with ecological divergence of postglacial freshwater fish species into multiple sympatric ecomorphs, focusing on species in the well-studied, widespread, and co-distributed genera Gasterosteus, Salvelinus, and Coregonus (stickleback, charr, and whitefish, respectively). We found that the presence of multiple sympatric ecomorphs tended to be associated with increasing lake surface area, maximum depth, and nutrient availability. In addition, predation, competition, and prey availability were suggested to play a role in divergence into multiple ecomorphs, but the effects of biotic factors require further study. Although we identified several environmental factors correlated with the presence of multiple ecomorphs, there were substantial data gaps across species and regions. An improved understanding of these systems may provide insight into both generalizable environmental factors involved in speciation in other systems, and potential ecological and evolutionary responses of species complexes when these variables are altered by environmental change.

Gut Microbiota Uniqueness Is Associated with Lake Size, a Proxy for Diet Diversity, in Stickleback Fish

AbstractOrganismal divergence can be driven by differential resource use and adaptation to different trophic niches. Variation in diet is a major factor shaping the gut microbiota, which is crucial for many aspects of their hosts' biology. However, it remains largely unknown how host diet diversity affects the gut microbiota, and it could be hypothesized that trophic niche width is positively associated with gut microbiota diversity. To test this idea, we sequenced the 16S ribosomal RNA gene from intestinal tissue of 14 threespine stickleback populations from lakes of varying size on Vancouver Island, Canada, that have been shown to differ in trophic niche width. Using lake size as a proxy for trophic ecology, we found evidence for higher gut microbiota uniqueness among individuals from populations with broader trophic niches. While these results suggest that diet diversity might promote gut microbiota diversity, additional work investigating diet and gut microbiota variation of the same host organisms will be necessary. Yet our results motivate the question of how host population diversity (e.g., ecological, morphological, genetic) might interact with the gut microbiota during the adaptation to ecological niches.

Opening a can of worms: a test of the co-infection facilitation hypothesis

Parasitic infections are a global occurrence and impact the health of many species. Coinfections, where two or more species of parasite are present in a host, are a common phenomenon across species. Coinfecting parasites can interact directly or indirectly via their manipulation of (and susceptibility to) the immune system of their shared host. Helminths, such as the cestode Schistocephalus solidus, are well known to suppress immunity of their host (threespine stickleback, Gasterosteus aculeatus), potentially facilitating other parasite species. Yet, hosts can evolve a more robust immune response (as seen in some stickleback populations), potentially turning facilitation into inhibition. Using wild-caught stickleback from 20 populations with non-zero S. solidus prevalence, we tested an a priori hypothesis that S. solidus infection facilitates infection by other parasites. Consistent with this hypothesis, individuals with S. solidus infections have 18.6% higher richness of other parasites compared to S. solidus-uninfected individuals from the same lakes. This facilitation-like trend is stronger in lakes where S. solidus is particularly successful but is reversed in lakes with sparse and smaller cestodes (indicative of stronger host immunity). These results suggest that a geographic mosaic of host-parasite co-evolution might lead to a mosaic of between-parasite facilitation/inhibition effects.

Temporal change in floral availability leads to periods of resource limitation and affects diet specificity in a generalist pollinator

Generalist species are core components of ecological networks and crucial for the maintenance of biodiversity. Generalist species and networks are expected to be more resilient, and therefore understanding the dynamics of specialization and generalization in ecological networks is a key focus in a time of rapid global change. Whilst diet generalization is frequently studied, our understanding of how it changes over time is limited. Here we explore temporal variation in diet specificity in the honeybee (Apis mellifera), using pollen DNA metabarcoding of honey samples, through the foraging season, over two years. We find that, overall, honeybees are generalists that visit a wide range of plants, but there is temporal variation in the degree of specialization. Temporal specialization of honeybee colonies corresponds to periods of resource limitation, identified as a lack of honey stores. Honeybees experience a lack of preferred resources in June when switching from flowering trees in spring to shrubs and herbs in summer. Investigating temporal patterns in specialization can identify periods of resource limitation that may lead to species and network vulnerability. Diet specificity must therefore be explored at different temporal scales in order to fully understand species and network stability in the face of ecological change.

Assessing the validity of fecal sampling for characterizing variation in threespine stickleback's gut microbiota

The gut microbiota is crucial for many aspects of their hosts' biology, and it has been characterized for many species across the animal kingdom. Yet, we still don't have a good understanding of whether non-lethal sampling can accurately capture the diversity of gut-associated bacterial communities, as estimated from lethal sampling of intestinal tissue. We further lack knowledge on whether non-lethal sampling methods are suitable for detecting gut microbiota shifts associated with changes in environmental factors (e.g., diet). We addressed these questions in threespine stickleback fish, a model system for evolutionary ecology, by comparing bacterial communities from intestinal tissue and feces. Despite some differences in community composition between the two sample types and considerable temporal variation among fecal samples, bacterial communities appear to largely overlap. Further, we detected consistent and significant changes of fecal bacterial communities associated with an experimental diet manipulation. This suggests that fecal sampling can represent an adequate non-lethal method to characterize the gut microbiota of threespine stickleback, but additional studies will be necessary before drawing general conclusions regarding the validity of fecal sampling for gut microbiota studies. To this end, we give recommendations to improve the characterization of the gut microbiota via fecal sampling. Fecal sampling allows studying temporal gut microbiota shifts associated with environmental change at the individual level, which increases opportunities for future experimental gut microbiota research.

Contemporary changes in phenotypic variation, and the potential consequences for eco-evolutionary dynamics

Most studies assessing rates of phenotypic change focus on population mean trait values, whereas a largely overlooked additional component is changes in population trait variation. Theoretically, eco-evolutionary dynamics mediated by such changes in trait variation could be as important as those mediated by changes in trait means. To date, however, no study has comprehensively summarised how phenotypic variation is changing in contemporary populations. Here, we explore four questions using a large database: How do changes in trait variances compare to changes in trait means? Do different human disturbances have different effects on trait variance? Do different trait types have different effects on changes in trait variance? Do studies that established a genetic basis for trait change show different patterns from those that did not? We find that changes in variation are typically small; yet we also see some very large changes associated with particular disturbances or trait types. We close by interpreting and discussing the implications of our findings in the context of eco-evolutionary studies.

Opening a can of worms: a test of the coinfection facilitation hypothesis

Parasitic infections are a global occurrence and impact the health of many species. Coinfections, where two or more species of parasite are present in a host, are a common phenomenon across species. Coinfecting parasites can interact directly, or indirectly via their manipulation of (and susceptibility to) the immune system of their shared host. Helminths, such as the cestode Schistocephalus solidus , are well known to suppress immunity of their host (threespine stickleback, Gasterosteus aculeatus ), potentially facilitating other parasite species. Yet, hosts can evolve a more robust immune response (as seen in some stickleback populations), potentially turning facilitation into inhibition. Using wild-caught stickleback from 21 populations with non-zero S. solidus prevalence, we tested an a priori hypothesis that S. solidus infection facilitates infection by other parasites. Consistent with this hypothesis, individuals with S. solidus infections have 18.6% higher richness of other parasites, compared to S. solidus -uninfected individuals from the same lakes. This facilitation-like trend is stronger in lakes where S. solidus is particularly successful but is reversed in lakes with sparse and smaller cestodes (indicative of stronger host immunity). These results suggest that a geographic mosaic of host-parasite coevolution might lead to a mosaic of between-parasite facilitation/inhibition effects.

Eukaryotes contribute more than bacteria to the recovery of freshwater ecosystem functions under different drought durations

Global climate change mostly impacts river ecosystems by affecting microbial biodiversity and ecological functions. Considering the high functional redundancy of microorganisms, the unknown relationship between biodiversity and ecosystem functions obstructs river ecological research, especially under the influence of increasing weather extremes, such as in intermittent rivers and ephemeral streams (IRES). Herein, dry-wet alternation experiments were conducted in artificial stream channels for 25 and 90 days of drought, both followed by 20 days of rewetting. The dynamic recovery of microbial biodiversity and ecosystem functions (represented by ecosystem metabolism and denitrification rate) were determined to analyse biodiversity-ecosystem-function (BEF) relationships after different drought durations. There was a significant difference between bacterial and eukaryotic biodiversity recovery after drought. Eukaryotic biodiversity was more sensitive to drought duration than bacterial, and the eukaryotic network was more stable under dry-wet alternations. Based on the establishment of partial least squares path models, we found that eukaryotic biodiversity has a stronger effect on ecosystem functions than bacteria after long-term drought. Indeed, this work represents a significant step forward for further research on the ecosystem functions of IRES, especially emphasizing the importance of eukaryotic biodiversity in the BEF relationship.

Seasonal variation of population and individual dietary niche in the avivorous bat, Ia io

The variation in niche breadth can affect how species respond to environmental and resource changes. However, there is still no clear understanding of how seasonal variability in food resources impacts the variation of individual dietary diversity, thereby affecting the dynamics of a population's dietary niche breadth. Optimal foraging theory (OFT) and the niche variation hypothesis (NVH) predict that when food resources are limited, the population niche breadth will widen or narrow due to increased within-individual dietary diversity and individual specialization or reduced within-individual dietary diversity, respectively. Here, we used DNA metabarcoding to examine the composition and seasonality of diets of the avivorous bat Ia io. Furthermore, we investigated how the dietary niches changed among seasons and how the population niche breadth changed when the availability of insect resources was reduced in autumn. We found that there was differentiation in dietary niches among seasons and a low degree of overlap, and the decrease of insect resource availability and the emergence of ecological opportunities of nocturnal migratory birds might drive dietary niche shifts toward birds in I. io. However, the population's dietary niche breadth did not broaden by increasing the within-individual dietary diversity or individual specialization, but rather became narrower by reducing dietary diversity via predation on bird resources that served as an ecological opportunity when insect resources were scarce in autumn. Our findings were consistent with the predictions of OFT, because birds as prey for bats provided extremely different resources from those of insects in size and nutritional value. Our work highlights the importance of size and quality of prey resources along with other factors (i.e., physiological, behavioral, and life-history traits) in dietary niche variation.

Integrating isotopic and nutritional niches reveals multiple dimensions of individual diet specialisation in a marine apex predator

Dietary specialisations are important determinants of ecological structure, particularly in species with high per-capita trophic influence like marine apex predators. These species are, however, among the most challenging in which to establish spatiotemporally integrated diets. We introduce a novel integration of stable isotopes with a multidimensional nutritional niche framework that addresses the challenges of establishing spatiotemporally integrated nutritional niches in wild populations, and apply the framework to explore individual diet specialisation in a marine apex predator, the white shark Carcharodon carcharias. Sequential tooth files were sampled from juvenile white sharks to establish individual isotopic (δ-space; δ¹³ C, δ¹⁵ N, δ³⁴ S) niche specialisation. Bayesian mixing models were then used to reveal individual-level prey (p-space) specialisation, and further combined with nutritional geometry models to quantify the nutritional (N-space) dimensions of individual specialisation, and their relationships to prey use. Isotopic and mixing model analyses indicated juvenile white sharks as individual specialists within a broader, generalist, population niche. Individual sharks differed in their consumption of several important mesopredator species, which suggested among-individual variance in trophic roles in either pelagic or benthic food webs. However, variation in nutrient intakes was small and not consistently correlated with differences in prey use, suggesting white sharks as nutritional specialists and that individuals could use functionally and nutritionally different prey as complementary means to achieve a common nutritional goal. We identify how degrees of individual specialisation can differ between niche spaces (δ-, p- or N-space), the physiological and ecological implications of this, and argue that integrating nutrition can provide stronger, mechanistic links between diet specialisation and its intrinsic (fitness/performance) and extrinsic (ecological) outcomes. Our time-integrated framework is adaptable for examining the nutritional consequences and drivers of food use variation at the individual, population or species level.

Genome-phenotype-environment associations identify signatures of selection in a panmictic population of threespine stickleback

Adaptive genetic divergence occurs when selection imposed by the environment causes the genomic component of the phenotype to differentiate. However, genomic signatures of natural selection are usually identified without information on which trait is responding to selection by which selective agent(s). Here, we integrate whole-genome sequencing with phenomics and measures of putative selective agents to assess the extent of adaptive divergence in threespine stickleback occupying the highly heterogeneous lake Mývatn, NE Iceland. We find negligible genome wide divergence, yet multiple traits (body size, gill raker structure and defence traits) were divergent along known ecological gradients (temperature, predatory bird densities and water depth). SNP based heritability of all measured traits was high (h²  = 0.42-0.65), indicating adaptive potential for all traits. Environment-association analyses further identified thousands of loci putatively involved in selection, related to genes linked to, for instance, neuron development and protein phosphorylation. Finally, we found that loci linked to water depth were concurrently associated with pelvic spine length variation - supporting the conclusion that divergence in pelvic spine length occurred in the face of gene flow. Our results suggest that whilst there is substantial genetic variation in the traits measured, phenotypic divergence of Mývatn stickleback is mostly weakly associated with environmental gradients, potentially as a result of substantial gene flow. Our study illustrates the value of integrative studies that combine genomic assays of multivariate trait variation with landscape genomics.

An AI-based framework for studying visual diversity of urban neighborhoods and its relationship with socio-demographic variables

This study presents a framework to study quantitatively geographical visual diversities of urban neighborhood from a large collection of street-view images using an Artificial Intelligence (AI)-based image segmentation technique. A variety of diversity indices are computed from the extracted visual semantics. They are utilized to discover the relationships between urban visual appearance and socio-demographic variables. This study also validates the reliability of the method with human evaluators. The methodology and results obtained from this study can potentially be used to study urban features, locate houses, establish services, and better operate municipalities.

Evolutionary gain and loss of a pathological immune response to parasitism

Parasites impose fitness costs on their hosts. Biologists often assume that natural selection favors infection-resistant hosts. Yet, when the immune response itself is costly, theory suggests that selection may sometimes favor loss of resistance, which may result in alternative stable states where some populations are resistant and others are tolerant. Intraspecific variation in immune costs is rarely surveyed in a manner that tests evolutionary patterns, and there are few examples of adaptive loss of resistance. Here, we show that when marine threespine stickleback colonized freshwater lakes, they gained resistance to the freshwater-associated cestode Schistocephalus solidus. Extensive peritoneal fibrosis and inflammation are a commonly observed phenotype that contributes to suppression of cestode growth and viability but also imposes a substantial cost on fecundity. Combining genetic mapping and population genomics, we find that opposing selection generates immune system differences between tolerant and resistant populations, consistent with divergent optimization.

Winter ecology of specialist and generalist morphs of European whitefish, Coregonus lavaretus, in subarctic northern Europe

European whitefish is a model species for adaptive radiation of fishes in temperate and subarctic lakes. In northern Europe the most commonly observed morphotypes are a generalist (LSR) morph and a pelagic specialist (DR) morph. The evolution of a pelagic specialist morph is something of an enigma, however, as this region is characterized by long, dark winters with pelagic primary production limited to a brief window in late summer. We conducted the first winter-based study of polymorphic whitefish populations to determine the winter ecology of both morphs, and we combined seasonal diet and stable isotope analysis with several proxies of condition in three polymorphic whitefish populations. The generalist LSR morph fed on benthic and pelagic prey in summer but was solely reliant on benthic prey in winter. This was associated with a noticeable but moderate reduction in condition, lipid content and stomach fullness in winter relative to summer. In contrast, the DR whitefish occupied a strict pelagic niche in both seasons. A significant reduction in pelagic prey during winter resulted in severe decrease in condition, lipid content and stomach fullness in DR whitefish in winter relative to summer, with the pelagic morph apparently approaching starvation in winter. We suggest that this divergent approach to seasonal foraging is associated with the divergent life-history traits of both morphs.

Population niche width is driven by within-individual niche expansion and individual specialization in introduced brook trout in mountain lakes

The width of a population's resource-use niche is determined by individual diet breadth ("within-individual component") and the degree of niche partitioning between individuals ("between-individual component"). The balance between these two factors affects ecological stability and evolutionary trajectories, and may shift as ecological opportunity permits broader population niches. Lakes in California's Sierra Nevada Mountains vary in resource diversity for introduced brook trout (Salvelinus fontinalis) due to elevation, lake morphometry, and watershed features. We compared the relative contributions of within- and between-individual niche components to two measures of the dietary niches of thirteen populations of brook trout: prey taxonomic composition and prey size distribution. For both taxonomic and size diversity of fish diets, population niche width was positively related to both the within- and between-individual components. For taxonomic diversity, the two components increased in parallel, while for size diversity, the between-individual component became more important relative to the within-individual component in populations with the greatest niche widths. Our results support the Niche Variation Hypothesis that populations with broader niches are more heterogeneous among individuals and show that individual niche width and individual specialization can operate in parallel to expand the population niche.

The evolution of trait variance creates a tension between species diversity and functional diversity

It seems intuitively obvious that species diversity promotes functional diversity: communities with more plant species imply more varied plant leaf chemistry, more species of crops provide more kinds of food, etc. Recent literature has nuanced this view, showing how the relationship between the two can be modulated along latitudinal or environmental gradients. Here we show that even without such effects, the evolution of functional trait variance can erase or even reverse the expected positive relationship between species- and functional diversity. We present theory showing that trait-based eco-evolutionary processes force species to evolve narrower trait breadths in more tightly packed, species-rich communities, in their effort to avoid competition with neighboring species. This effect is so strong that it leads to an overall reduction in trait space coverage whenever a new species establishes. Empirical data from land snail communities on the Galápagos Islands are consistent with this claim. The finding that the relationship between species- and functional diversity can be negative implies that trait data from species-poor communities may misjudge functional diversity in species-rich ones, and vice versa.

The positive relationship between species diversity and functional diversity has been shown to vary. Here, the authors use theoretical models and data from Galápagos land snail communities to show how eco-evolutionary processes can force species to evolve narrower trait breadths in more species-rich communities to avoid competition, creating a negative relationship.

A test of frequency‐dependent selection in the evolution of a generalist phenotype

A solitary population of consumers frequently evolves to the middle of a resource gradient and an intermediate mean phenotype compared to a sympatric pair of competing species that diverge to either side via character displacement. The forces governing the distribution of phenotypes in these allopatric populations, however, are little investigated. Theory predicts that the intermediate mean phenotype of the generalist should be maintained by negative frequency‐dependent selection, whereby alternate extreme phenotypes are favored because they experience reduced competition for resources when rare. However, the theory makes assumptions that are not always met, and alternative explanations for an intermediate phenotype are possible. We provide a test of this prediction in a mesocosm experiment using threespine stickleback that are ecologically and phenotypically intermediate between the more specialized stickleback species that occur in pairs. We manipulated the frequency distribution of phenotypes in two treatments and then measured effects on a focal intermediate population. We found a slight frequency‐dependent effect on survival in the predicted direction but not on individual growth rates. This result suggests that frequency‐dependent selection might be a relatively weak force across the range of phenotypes within an intermediate population and we suggest several general reasons why this might be so. We propose that allopatric populations might often be maintained at an intermediate phenotype instead by stabilizing or fluctuating directional selection.

Analysis of ancestry heterozygosity suggests that hybrid incompatibilities in threespine stickleback are environment dependent

Hybrid incompatibilities occur when interactions between opposite ancestry alleles at different loci reduce the fitness of hybrids. Most work on incompatibilities has focused on those that are "intrinsic," meaning they affect viability and sterility in the laboratory. Theory predicts that ecological selection can also underlie hybrid incompatibilities, but tests of this hypothesis using sequence data are scarce. In this article, we compiled genetic data for F2 hybrid crosses between divergent populations of threespine stickleback fish (Gasterosteus aculeatus L.) that were born and raised in either the field (seminatural experimental ponds) or the laboratory (aquaria). Because selection against incompatibilities results in elevated ancestry heterozygosity, we tested the prediction that ancestry heterozygosity will be higher in pond-raised fish compared to those raised in aquaria. We found that ancestry heterozygosity was elevated by approximately 3% in crosses raised in ponds compared to those raised in aquaria. Additional analyses support a phenotypic basis for incompatibility and suggest that environment-specific single-locus heterozygote advantage is not the cause of selection on ancestry heterozygosity. Our study provides evidence that, in stickleback, a coarse-albeit indirect-signal of environment-dependent hybrid incompatibility is reliably detectable and suggests that extrinsic incompatibilities can evolve before intrinsic incompatibilities.

Genomic and environmental influences on resilience in a cold-water fish near the edge of its range

Small, isolated populations present a challenge for conservation. The dueling effects of selection and drift in a limited pool of genetic diversity make the responses of small populations to environmental perturbations erratic and difficult to predict. This is particularly true at the edge of a species range, where populations often persist at the limits of their environmental tolerances. Populations of cisco, Coregonus artedi, in inland lakes have experienced numerous extirpations along the southern edge of their range in recent decades, which are thought to result from environmental degradation and loss of cold, well-oxygenated habitat as lakes warm. Yet, cisco extirpations do not show a clear latitudinal pattern, suggesting that local environmental factors and potentially local adaptation may influence resilience. Here, we used genomic tools to investigate the nature of this pattern of resilience. We used restriction site-associated DNA capture (Rapture) sequencing to survey genomic diversity and differentiation in southern inland lake cisco populations and compared the frequency of deleterious mutations that potentially influence fitness across lakes. We also examined haplotype diversity in a region of the major histocompatibility complex involved in stress and immune system response. We correlated these metrics to spatial and environmental factors including latitude, lake size, and measures of oxythermal habitat and found significant relationships between genetic metrics and broad and local factors. High levels of genetic differentiation among populations were punctuated by a phylogeographic break and residual patterns of isolation-by-distance. Although the prevalence of deleterious mutations and inbreeding coefficients was significantly correlated with latitude, neutral and non-neutral genetic diversity were most strongly correlated with lake surface area. Notably, differences among lakes in the availability of estimated oxythermal habitat left no clear population genomic signature. Our results shed light on the complex dynamics influencing these isolated populations and provide valuable information for their conservation.

Species interactions, environmental gradients and body size shape population niche width

Competition for shared resources is commonly assumed to restrict population-level niche width of coexisting species. However, the identity and abundance of coexisting species, the prevailing environmental conditions, and the individual body size may shape the effects of interspecific interactions on species' niche width. Here we study the effects of interspecific and intraspecific interactions, lake area and altitude, and fish body size on the trophic niche width and resource use of a generalist predator, the littoral-dwelling large, sparsely rakered morph of European whitefish (Coregonus lavaretus; hereafter LSR whitefish). We use stable isotope, diet and survey fishing data from 14 subarctic lakes along an environmental gradient in northern Norway. The isotopic niche width of LSR whitefish showed a humped-shaped relationship with increasing relative abundance of sympatric competitors, suggesting widest population niche at intermediate intensity of interspecific interactions. The isotopic niche width of LSR whitefish tended to decrease with increasing altitude, suggesting reduced niche in colder, less productive lakes. LSR whitefish typically shifted to a higher trophic position and increased reliance on littoral food resources with increasing body size, although between-lake differences in ontogenetic niche shifts were evident. In most lakes, LSR whitefish relied less on littoral food resources than coexisting fishes and the niche overlap between sympatric competitors was most evident among relatively large individuals (>250 mm). Individual niche variation was highest among >200 mm long LSR whitefish, which likely have escaped the predation window of sympatric predators. We demonstrate that intermediate intensity of interspecific interactions may broaden species' niche width, whereas strong competition for limited resources and high predation risk may suppress niche width in less productive environments. Acknowledging potential humped-shaped relationships between population niche width and interspecific interactions can help us understand species' responses to environmental disturbance (e.g. climate change and species invasions) as well as the driving forces of niche specialization.

Male and female reproductive fitness costs of an immune response in natural populations. Evolution 2021;75:2509-2523. [PMID: 33991339 PMCID: PMC8488946 DOI: 10.1111/evo.14266]

Parasites can mediate host fitness both directly, via effects on survival and reproduction, or indirectly by inducing host immune defense with costly side-effects. The evolution of immune defense is determined by a complex interplay of costs and benefits of parasite infection and immune response, all of which may differ for male and female hosts in sexual lineages. Here, we examine fitness costs associated with an inducible immune defense in a fish-cestode host-parasite system. Cestode infection induces peritoneal fibrosis in threespine stickleback (Gasterosteus aculeatus), constraining cestode growth and sometimes encasing and killing the parasite. Surveying two wild populations of stickleback, we confirm that the presence of fibrosis scar tissue is associated with reduced parasite burden in both male and female fish. However, fibrotic fish had lower foraging success and reproductive fitness (reduced female egg production and male nesting success), indicating strong costs of the lingering immunopathology. Consistent with substantial sexually concordant fitness effects of immune response, we find alignment of multivariate selection across the sexes despite sexual antagonism over morphological shape. Although both sexes experienced costs of fibrosis, the net impacts are unequal because in the two study populations females had higher cestode exposure. To evaluate whether this difference in risk should drive sex-specific immune strategies, we analyze a quantitative genetic model of host immune response to a trophically transmitted parasite. The model and empirical data illustrate how shared costs and benefits of immune response lead to shared evolutionary interests of male and female hosts, despite unequal infection risks across the sexes.

Sick of eating: Eco-evo-immuno dynamics of predators and their trophically acquired parasites

When predators consume prey, they risk becoming infected with their prey's parasites, which can then establish the predator as a secondary host. A predator population's diet therefore influences what parasites it is exposed to, as has been repeatedly shown in many species such as threespine stickleback (Gasterosteus aculeatus) (more benthic‐feeding individuals obtain nematodes from oligocheate prey, whereas limnetic‐feeding individuals catch cestodes from copepod prey). These differing parasite encounters, in turn, determine how natural selection acts on the predator's immune system. We might therefore expect that ecoevolutionary dynamics of a predator's diet (as determined by its ecomorphology) should drive correlated evolution of its immune traits. Conversely, the predator's immunity to certain parasites might alter the relative costs and benefits of different prey, driving evolution of its ecomorphology. To evaluate the potential for ecological morphology to drive evolution of immunity, and vice versa, we use a quantitative genetics framework coupled with an ecological model of a predator and two prey species (the diet options). Our analysis reveals fundamental asymmetries in the evolution of ecomorphology and immunity. When ecomorphology rapidly evolves, it determines how immunity evolves, but not vice versa. Weak trade‐offs in ecological morphology select for diet generalists despite strong immunological trade‐offs, but not vice versa. Only weak immunological trade‐offs can explain negative diet‐infection correlations across populations. The analysis also reveals that eco‐evo‐immuno feedbacks destabilize population dynamics when trade‐offs are sufficiently weak and heritability is sufficiently high. Collectively, these results highlight the delicate interplay between multivariate trait evolution and the dynamics of ecological communities.

Environmental heterogeneity promotes individual specialisation in habitat selection in a widely distributed seabird

Individual specialisations in behaviour are predicted to arise where divergence benefits fitness. Such specialisations are more likely in heterogeneous environments where there is both greater ecological opportunity and competition-driven frequency dependent selection. Such an effect could explain observed differences in rates of individual specialisation in habitat selection, as it offers individuals an opportunity to select for habitat types that maximise resource gain while minimising competition; however, this mechanism has not been tested before. Here, we use habitat selection functions to quantify individual specialisations while foraging by black-legged kittiwakes Rissa tridactyla, a marine top predator, at 15 colonies around the United Kingdom and Ireland, along a gradient of environmental heterogeneity. We find support for the hypothesis that individual specialisations in habitat selection while foraging are more prevalent in heterogeneous environments. This trend was significant across multiple dynamic habitat variables that change over short time-scales and did not arise through site fidelity, which highlights the importance of environmental processes in facilitating behavioural adaptation by predators. Individual differences may drive evolutionary processes, and therefore these results suggest that there is broad scope for the degree of environmental heterogeneity to determine current and future population, species and community dynamics.

Testing the relationship between intraspecific competition and individual specialization across both behavior and diet

Individual specialization within populations is increasingly recognized as important in both ecology and evolution, but researchers working on intraspecific variation in behavior and diet infrequently interact. This may be because individual specialization on diet and behavior was historically difficult to investigate simultaneously on the same individuals. However, approaches like stable isotope analysis that allow hindcasting past field diets for laboratory organisms may provide opportunities to unite these areas of inquiry. Here, we tested the role of intraspecific competition on individual specialization through analysis of both behavior and diet simultaneously. We focused on intraspecific competition as a mechanism that might drive individual specialization of both diet and behavior. We conducted this study in Vilas County, Wisconsin, United States (US), using rusty crayfish Faxonius rusticus from six lakes across a relative abundance gradient. We conducted six assays to measure individual specialization of behavior and used stable isotope analysis to measure individual specialization of diet. We then related both measures of individual specialization to relative abundance of F. rusticus using linear and quadratic models. We found a unimodal relationship between intraspecific competition and individual specialization of diet for F. rusticus, likely because some preferred resources are unavailable to specialize on at the highest densities of this well-studied crayfish invader. Conversely, we found greater support for a linear relationship between individual specialization of behavior and intraspecific competition, perhaps because specialization by behavior is not inherently resource-limited. Our results show that dietary and behavioral specialization may exhibit different responses to increased intraspecific competition, and demonstrate a potential technique that can be used to investigate individual specialization of diet and behavior simultaneously for the same individuals and populations.

Individual environmental niches in mobile organisms

Individual variation is increasingly recognized as a central component of ecological processes, but its role in structuring environmental niche associations remains largely unknown. Species' responses to environmental conditions are ultimately determined by the niches of single individuals, yet environmental associations are typically captured only at the level of species. Here, we develop scenarios for how individual variation may combine to define the compound environmental niche of populations, use extensive movement data to document individual environmental niche variation, test associated hypotheses of niche configuration, and examine the consistency of individual niches over time. For 45 individual white storks (Ciconia ciconia; 116 individual-year combinations), we uncover high variability in individual environmental associations, consistency of individual niches over time, and moderate to strong niche specialization. Within populations, environmental niches follow a nested pattern, with individuals arranged along a specialist-to-generalist gradient. These results reject common assumptions of individual niche equivalency among conspecifics, as well as the separation of individual niches into disparate parts of environmental space. These findings underscore the need for a more thorough consideration of individualistic environmental responses in global change research.

The interaction of resource use and gene flow on the phenotypic divergence of benthic and pelagic morphs of Icelandic Arctic charr (Salvelinus alpinus)

Conceptual models of adaptive divergence and ecological speciation in sympatry predict differential resource use, phenotype-environment correlations, and reduced gene flow among diverging phenotypes. While these predictions have been assessed in past studies, connections among them have rarely been assessed collectively. We examined relationships among phenotypic, ecological, and genetic variation in Arctic charr (Salvelinus alpinus) from six Icelandic localities that have undergone varying degrees of divergence into sympatric benthic and pelagic morphs. We characterized morphological variation with geometric morphometrics, tested for differential resource use between morphs using stable isotopes, and inferred the amount of gene flow from single nucleotide polymorphisms. Analysis of stable isotopic signatures indicated that sympatric morphs showed similar difference in resource use across populations, likely arising from the common utilization of niche space within each population. Carbon isotopic signature was also a significant predictor of individual variation in body shape and size, suggesting that variation in benthic and pelagic resource use is associated with phenotypic variation. The estimated percentage of hybrids between sympatric morphs varied across populations (from 0% to 15.6%) but the majority of fish had genotypes (ancestry coefficients) characteristic of pure morphs. Despite evidence of reduced gene flow between sympatric morphs, we did not detect the expected negative relationship between divergence in resource use and gene flow. Three lakes showed the expected pattern, but morphs in the fourth showed no detectable hybridization and had relatively low differences in resource use between them. This coupled with the finding that resource use and genetic differentiation had differential effects on body shape variation across populations suggests that reproductive isolation maintains phenotypic divergence between benthic and pelagic morphs when the effects of resource use are relatively low. Our ability to assess relationships between phenotype, ecology, and genetics deepens our understanding of the processes underlying adaptive divergence in sympatry.

Maternal effects and fitness consequences of individual variation in bottlenose dolphins' ecological niche

The niche describes the ecological and social environment that an organism lives in, as well as the behavioural tactics used to interact with its environment. A species niche is key to both ecological and evolutionary processes, including speciation, and has therefore been a central focus in ecology. Recent evidence, however, points to considerable individual variation in a species' or population's niche use, although how this variation evolves or is maintained remains unclear. We used a large longitudinal dataset to investigate the drivers and maintenance of individual variation in bottlenose dolphins' Tursiops aduncus niche. Specifically, we (a) characterised the extent of individual differences in habitat use, (b) identified whether there were maternal effects associated with this variation and (c) investigated the relationship between habitat use and calving success, a component of reproductive fitness. By examining patterns of habitat use, we provide evidence that individual dolphins vary consistently between one another in their niche. We further show that such individual variation is driven by a strong maternal effect. Finally, habitat use and calving success were not related, suggesting that use of different habitats results in similar fitness outcomes. Niche partitioning, maintained by maternal effects, likely facilitates the coexistence of multiple ecotypes within this population.

What evolutionary processes maintain MHC IIꞵ diversity within and among populations of stickleback?

Major Histocompatibility Complex (MHC) genes code for proteins that recognize foreign protein antigens to initiate T-cell-mediated adaptive immune responses. They are often the most polymorphic genes in vertebrate genomes. How evolution maintains this diversity remains of debate. Three main hypotheses seek to explain the maintenance of MHC diversity by invoking pathogen-mediated selection: heterozygote advantage, frequency-dependent selection, and fluctuating selection across landscapes or through time. Here, we use a large-scale field parasite survey in a stickleback metapopulation to test predictions derived from each of these hypotheses. We identify over 1000 MHC IIβ variants (alleles spanning paralogous genes) and find that many of them covary positively or negatively with parasite load, suggesting that these genes contribute to resistance or susceptibility. However, despite our large sample-size, we find no evidence for the widely cited stabilizing selection on MHC heterozygosity, in which individuals with an intermediate number of MHC variants have the lowest parasite burden. Nor do we observe a rare-variant advantage, or widespread fluctuating selection across populations. In contrast, we find that MHC diversity is best predicted by neutral genome-wide heterozygosity and between-population genomic divergence, suggesting neutral processes are important in shaping the pattern of metapopulation MHC diversity. Thus, although MHC IIβ is highly diverse and relevant to the type and intensity of macroparasite infection in these populations of stickleback, the main models of MHC evolution still provide little explanatory power in this system.

Among-individual diet variation within a lake trout ecotype: Lack of stability of niche use

In a polyphenic species, differences in resource use are expected among ecotypes, and homogeneity in resource use is expected within an ecotype. Yet, using a broad resource spectrum has been identified as a strategy for fishes living in unproductive northern environments, where food is patchily distributed and ephemeral. We investigated whether specialization of trophic resources by individuals occurred within the generalist piscivore ecotype of lake trout from Great Bear Lake, Canada, reflective of a form of diversity. Four distinct dietary patterns of resource use within this lake trout ecotype were detected from fatty acid composition, with some variation linked to spatial patterns within Great Bear Lake. Feeding habits of different groups within the ecotype were not associated with detectable morphological or genetic differentiation, suggesting that behavioral plasticity caused the trophic differences. A low level of genetic differentiation was detected between exceptionally large-sized individuals and other piscivore individuals. We demonstrated that individual trophic specialization can occur within an ecotype inhabiting a geologically young system (8,000-10,000 yr BP), a lake that sustains high levels of phenotypic diversity of lake trout overall. The characterization of niche use among individuals, as done in this study, is necessary to understand the role that individual variation can play at the beginning of differentiation processes.

Weighted individual-resource networks in prey-predator systems: the role of prey availability on the emergence of modular structures

Ecological networks, usually depicting interactions among species, have been recently down-scaled to the individual level, permitting description of patterns of inter-individual resource variation that are usually hindered at the species level. Optimal diet theory (ODT) models, applied to prey-predator systems, predict different patterns of nestedness and modularity in the network, depending on the available resources and intra-specific competition. The effect of resource availability on the emergence of networks structures, and ODT framework, has not yet fully been clarified. Here, we analyzed the structural patterns of individual-resource networks in 3 species of Mediterranean salamanders, in relation to changes in prey availability. We used weighted individual-resource network metrics to interpret the observed patterns, according to 3 ODT models. We found significant nestedness recurring in our study system, indicating that both selective and opportunistic individuals occur in the same population. Prey diversity, rather than abundance, was apparently related to inter-individual resource variation and promoted the emergence of significant modularity within all networks. The observed patterns of nestedness and modularity, together with the variation in resource diversity and intra-specific competition, are in agreement with the distinct preferences model of ODT. These findings suggest that in the focal prey-predator systems, individuals were able to perceive changes in prey diversity and to exploit in different ways the variations in composition of available resources, shifting their diet assembly rules accordingly. Our findings also confirm that the use of weighted individual-resource networks, in prey-predator systems, allows to disclose dynamics that are masked at the species or population level.

The genomic signature of ecological divergence along the benthic-limnetic axis in allopatric and sympatric threespine stickleback

The repeated occurrence of similar phenotypes in independent lineages (i.e., parallel evolution) in response to similar ecological conditions can provide compelling insights into the process of adaptive evolution. An intriguing question is to what extent repeated phenotypic changes are underlain by repeated changes at the genomic level and whether patterns of genomic divergence differ with the geographic context in which populations evolve. Here, we combined genomic, morphological and ecological data sets to investigate the genomic signatures of divergence across populations of threespine stickleback (Gasterosteus aculeatus) that adapted to contrasting ecological niches (benthic or limnetic) in either sympatry or allopatry. We found that genome-wide differentiation (F_ST ) was an order of magnitude higher and substantially more repeatable for sympatric benthic and limnetic specialists compared to allopatric populations with similar levels of ecological divergence. We identified genomic regions consistently differentiated between sympatric ecotypes that were also differentiated between or associated with benthic vs. limnetic niche in allopatric populations. These candidate regions were enriched on three chromosomes known to be involved in the benthic-limnetic divergence of threespine stickleback. Some candidate regions overlapped with QTL for body shape and trophic traits such as gill raker number, traits that strongly differ between benthic and limnetic ecotypes. In summary, our study shows that magnitude and repeatability of genomic signatures of ecological divergence in threespine stickleback highly depend on the geographic context. The identified candidate regions provide starting points to identify functionally important genes for the adaptation to benthic and limnetic niches.

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.

Intraspecific variation and structuring of phenotype in a lake-dwelling species are driven by lake size and elevation

The fragmented, heterogeneous and relatively depauperate ecosystems of recently glaciated lakes present contrasting ecological opportunities for resident fish. Across a species, local adaptation may induce diverse and distinct phenotypic responses to various selection pressures. We tested for intraspecific phenotypic structuring by population in a common native lake-dwelling fish species across a medium-scale geographic region with considerable variation in lake types. We investigated potential lake-characteristic drivers of trophic morphology. Using geometric morphometric techniques, we quantified the head shapes of 759 adult brown trout (Salmo trutta L.) from 28 lakes and reservoirs across Scotland. Multivariate statistical analyses showed that almost all populations differed from one another. Trout from larger and deeper lakes had deeper, but shorter heads, and smaller eyes. Higher elevation lakes were associated with fish with shorter heads and jaws. Within-population shape variation also differed by population, and was positively correlated with lake surface area and depth. Trout within reservoirs differed subtly from those in natural lakes, having larger eyes, shorter jaws and greater variability. This study documents an extraordinary morphological variation between and within populations of brown trout, and demonstrates the role of the extrinsic environment in driving phenotypic structuring over a medium-scale and varied geographic area.

Predators weaken prey intraspecific competition through phenotypic selection

Predators have a key role shaping competitor dynamics in food webs. Perhaps the most obvious way this occurs is when predators reduce competitor densities. However, consumption could also generate phenotypic selection on prey that determines the strength of competition, thus coupling consumptive and trait-based effects of predators. In a mesocosm experiment simulating fish predation on damselflies, we found that selection against high damselfly activity rates - a phenotype mediating predation and competition - weakened the strength of density dependence in damselfly growth rates. A field experiment corroborated this finding and showed that increasing damselfly densities in lakes with high fish densities had limited effects on damselfly growth rates but generated a precipitous growth rate decline where fish densities were lower - a pattern expected because of spatial variation in selection imposed by predation. These results suggest that accounting for both consumption and selection is necessary to determine how predators regulate prey competitive interactions.