Hostname: page-component-7c8c6479df-24hb2 Total loading time: 0 Render date: 2024-03-19T05:10:48.780Z Has data issue: false hasContentIssue false

Context-dependent parasite infection affects trophic niche in populations of sympatric stickleback species

Published online by Cambridge University Press:  16 May 2022

Doko-Miles J. Thorburn*
Affiliation:
School of Biological and Chemical Sciences, Queen Mary University of London, London, UK Department of Life Sciences, Imperial College London, London, UK
Thijs M. P. Bal
Affiliation:
Faculty of Biosciences and Aquaculture, Nord University, Bodø, Norway
Io S. Deflem
Affiliation:
Laboratory of Biodiversity and Evolutionary Genomics, KU Leuven, Leuven, Belgium
Filip A. M. Volckaert
Affiliation:
Laboratory of Biodiversity and Evolutionary Genomics, KU Leuven, Leuven, Belgium
Christophe Eizaguirre
Affiliation:
School of Biological and Chemical Sciences, Queen Mary University of London, London, UK
Joost A. M. Raeymaekers
Affiliation:
Faculty of Biosciences and Aquaculture, Nord University, Bodø, Norway
*
Author for correspondence: Doko-Miles J. Thorburn, E-mail: d.m.j.thorburn@qmul.ac.uk

Abstract

How parasites alter host feeding ecology remains elusive in natural populations. A powerful approach to investigate the link between infection and feeding ecology is quantifying unique and shared responses to parasite infection in related host species within a common environment. Here, 9 pairs of sympatric populations of the three-spined and nine-spined stickleback fishes were sampled across a range of freshwater and brackish habitats to investigate how parasites alter host feeding ecology: (i) biotic and abiotic determinants of parasite community composition, and (ii) to what extent parasite infection correlates with trophic niche specialization of the 2 species, using stable isotope analyses (δ15N and δ13C). It was determined that parasite community composition and host parasite load varied among sites and species and were correlated with dissolved oxygen. It was also observed that the digenean Cyathocotyle sp.'s abundance, a common directly infecting parasite with a complex life cycle, correlated with host δ13C in a fish species-specific manner. In 6 sites, correlations were found between parasite abundance and their hosts' feeding ecology. These effects were location-specific and occasionally host species or host size-specific. Overall, the results suggest a relationship between parasite infection and host trophic niche which may be an important and largely overlooked ecological factor. The population specificity and variation in parasite communities also suggest this effect is multifarious and context-dependent.

Type
Research Article
Creative Commons
Creative Common License - CCCreative Common License - BY
This is an Open Access article, distributed under the terms of the Creative Commons Attribution licence (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted re-use, distribution and reproduction, provided the original article is properly cited.
Copyright
Copyright © The Author(s), 2022. Published by Cambridge University Press

Introduction

Parasites are strong selective agents on their hosts, reducing their fitness as well as altering their behaviour, life-history traits and habitat use (Milinski, Reference Milinski1984; Miura et al., Reference Miura, Kuris, Torchin, Hechinger and Chiba2006; Pagán et al., Reference Pagán, Alonso-Blanco and García-Arenal2008; Lefèvre et al., Reference Lefèvre, Lebarbenchon, Gauthier-Clerc, Missé, Poulin and Thomas2009; Barber and Huntingford, Reference Barber and Huntingford1995). Parasites are involved in approximately 75% of food web links (Lafferty et al., Reference Lafferty, Dobson and Kuris2006), and often exceed total biomass of top predators (Kuris et al., Reference Kuris, Hechinger, Shaw, Whitney, Aguirre-Macedo, Boch, Dobson, Dunham, Fredensborg, Huspeni, Lorda, Mababa, Mancini, Mora, Pickering, Talhouk, Torchin and Lafferty2008). Despite growing evidence that parasites have a strong impact on food web ecology (Lafferty et al., Reference Lafferty, Dobson and Kuris2006; Dobson et al., Reference Dobson, Lafferty, Kuris, Hechinger and Jetz2009; Anaya-Rojas et al., Reference Anaya-Rojas, Best, Brunner, Eizaguirre, Leal, Melián, Seehausen and Matthews2019), their role in altering an individual hosts' trophic niche specialization remains largely unknown (Araújo et al., Reference Araújo, Bolnick and Layman2011; Pegg et al., Reference Pegg, Andreou, Williams and Britton2015a, Reference Pegg, Andreou, Williams and Britton2017; Britton and Andreou, Reference Britton and Andreou2016; Lockley et al., Reference Lockley, Fouda, Correia, Taxonera, Nash, Fairweather, Reischig, Durão, Dinis, Roque, Lomba, Passos, Cameron, Stiebens and Eizaguirre2020). Such parasite-mediated effects on trophic niche specialization are unresolved because (i) the effects of parasites on altering host trophic niche may be system-specific, or (ii) parasites with different life-history strategies may alter trophic interactions in different ways. For example, trophically transmitted parasites often have complex life cycles and can manipulate host behaviour, therein modifying trophic interactions (Lafferty et al., Reference Lafferty, Allesina, Arim, Briggs, De Leo, Dobson, Dunne, Johnson, Kuris, Marcogliese, Martinez, Memmott, Marquet, McLaughlin, Mordecai, Pascual, Poulin and Thieltges2008; Hammerschmidt et al., Reference Hammerschmidt, Koch, Milinski, Chubb and Parker2009). Specifically, through infection and regulation in all trophic levels, including top predators, parasites can induce or inhibit trophic cascades affecting overall ecosystem functioning (Brunner et al., Reference Brunner, Anaya-Rojas, Matthews and Eizaguirre2017; Anaya-Rojas et al., Reference Anaya-Rojas, Best, Brunner, Eizaguirre, Leal, Melián, Seehausen and Matthews2019). It should be noted that hosts are not defenceless, but mounting an immunological response can be energetically costly, and therefore an appropriate feeding strategy is necessary to compensate the costs of infection (Lee et al., Reference Lee, Cory, Wilson, Raubenheimer and Simpson2006; Brunner et al., Reference Brunner, Schmid-Hempel and Barribeau2014, Reference Brunner, Anaya-Rojas, Matthews and Eizaguirre2017).

The theoretical effects of parasite infection on trophic niche specialization have been summarized into 2 main effects: constriction and divergence (Britton and Andreou, Reference Britton and Andreou2016). Trophic niche constriction occurs when a population subgroup specializes on a restricted breadth of food items that are also consumed by conspecific generalists. Indeed, parasite infection can result in the sub-population of infected hosts being nested within the overall host's population niche (Britton and Andreou, Reference Britton and Andreou2016; Pegg et al., Reference Pegg, Andreou, Williams and Britton2017; Villalobos et al., Reference Villalobos, Nava-Bolaños, De Fuentes-Vicente, Téllez-Rendón, Huerta, Martínez-Hernández, Rocha-Ortega, Gutiérrez-Cabrera, Ibarra-Cerdeña and Córdoba-Aguilar2019). For example, when comparing the trophic niche of the generalist common roach (Rutilus rutilus) and common bream (Abramis brama) infected by Ergasilius briani, infected individuals showed a smaller niche width that is nested within that of uninfected conspecifics (Pegg et al., Reference Pegg, Andreou, Williams and Britton2017). Conversely, the overall population niche might expand due to a sub-population of infected individuals shifting diet to food items not previously exploited by conspecifics, resulting in trophic niche divergence (Britton and Andreou, Reference Britton and Andreou2016). Competitive exclusion of infected individuals would manifest as an overall population niche expansion. For instance, common carp (Cyprinus carpio) infected with the non-native Asian tapeworm (Bothriocephalus acheilognathi) can become competitively disadvantaged upon infection, and alter their diet compared to uninfected conspecifics (Britton et al., Reference Britton, Pegg and Williams2011; Pegg et al., Reference Pegg, Andreou, Williams and Britton2015a).

Identifying whether parasite-induced niche constriction or divergence are a function of environmental conditions represents a logical next step in uncovering the causal mechanisms that underlie the relationship between parasitism and host trophic ecology. Indeed, abiotic environmental conditions such as variation in physico-chemistry affect important aspects of the ecology of host–parasite interactions, including host community structure, and parasite abundance and virulence (Blanar et al., Reference Blanar, Munkittrick, Houlahan, MacLatchy and Marcogliese2009; Karvonen et al., Reference Karvonen, Kristjánsson, Skúlason, Lanki, Rellstab and Jokela2013; Mahmud et al., Reference Mahmud, Bradley and MacColl2017; Deflem et al., Reference Deflem, Bennetsen, Opedal, Calboli, Ovaskainen, Van Thuyne, Volckaert and Raeymaekers2021). For example, reductions in dissolved oxygen (DO) can induce hypoxic stress in hosts which can result in ineffective immune responses, and is associated with increased parasite diversity and infection load in fish (Lacerda et al., Reference Lacerda, Roumbedakis, Bereta Junior, Nuñer, Petrucio and Martins2018; Ojwala et al., Reference Ojwala, Otachi and Kitaka2018; Abdel-Tawwab et al., Reference Abdel-Tawwab, Monier, Hoseinifar and Faggio2019). Conversely, acidification of aquatic environments can deplete intermediate host populations as in the case in many trematode parasites and their snail hosts, and is associated with reductions in parasite richness and virulence (Lafferty and Kuris, Reference Lafferty and Kuris1999; Young and MacColl, Reference Young and MacColl2017).

To disentangle the relative contributions of host-related and environmental drivers of parasite-mediated shifts in host trophic niche, studies can focus on natural systems with environmental gradients and multiple related species persisting in a common landscape. In general, if genetically and ecologically similar sympatric species respond in a similar manner, the response is likely constrained by niche availability in the environment. However, if species living in sympatry respond differently to a shared environment, the responses are likely linked to mechanisms unique to a species (Raeymaekers et al., Reference Raeymaekers, Chaturvedi, Hablützel, Verdonck, Hellemans, Maes, De Meester and Volckaert2017; Bal et al., Reference Bal, Llanos-garrido, Chaturvedi, Verdonck, Hellemans and Raeymaekers2021).

The three-spined (Gasterosteus aculeatus Linnaeus, 1758; Gasterosteidae) and nine-spined (Pungitius pungitius Linnaeus, 1758; Gasterosteidae) sticklebacks are 2 small teleost fishes, which diverged around 26 million years ago (Varadharajan et al., Reference Varadharajan, Rastas, Löytynoja, Matschiner, Calboli, Guo, Nederbragt, Jakobsen and Merilä2019). Often found living in sympatry, both species demonstrate examples of overlapping dietary preference and habitat use, similar population genetic structure and parallel phenotypic divergence (Copp and Kováč, Reference Copp and Kováč2003; Raeymaekers et al., Reference Raeymaekers, Chaturvedi, Hablützel, Verdonck, Hellemans, Maes, De Meester and Volckaert2017; Bal et al., Reference Bal, Llanos-garrido, Chaturvedi, Verdonck, Hellemans and Raeymaekers2021). Their macroparasites are well documented, with several parasites able to infect both species, and some that are host-specific (Kalbe et al., Reference Kalbe, Wegner and Reusch2002; Zander, Reference Zander2007; Henrich et al., Reference Henrich, Benesh and Kalbe2013). Responses to infection in both species use a similar set of immune genes (Lenz et al., Reference Lenz, Eizaguirre, Kalbe and Milinski2013). The shared ecological context among the species provides a good opportunity to investigate the relative roles of host species, their environment and parasite infection on trophic niche ecology.

Using natural coexisting populations of G. aculeatus and P. pungitius across a range of habitats in Belgium and the Netherlands, we investigated the role of parasites in altering the trophic niche of host species, including (i) quantifying parasite diversity and abundance among host populations and their relationship with abiotic environmental variables, (ii) using stable isotope analysis (SIA; δ 15N and δ 13C) to assess the size of host trophic niche and investigate its relationship with parasite infection and (iii) determine whether these effects are consistent among species and sampling site, or more context-dependent. These aims test our hypotheses that infected individuals will have constricted or divergent trophic niches compared with uninfected conspecifics, and whether constriction or divergence will be linked to additional environmental variables.

Materials and methods

Sample and data collection

The coastal lowlands of the Netherlands and Belgium contain both freshwater and brackish habitats. These habitats are often dominated by 2 coexisting and phylogenetically related fish species: the 3-spined (G. aculeatus) and the 9-spined (P. pungitius) sticklebacks. Five freshwater and 4 brackish sites were sampled in March–April 2018 using dipnets (Fig. 1). Sites were selected to represent a range of salinity and were only sampled if the target species were identified living in sympatry. Sampling started with recording conductivity (μS cm−1), pH, temperature (°C), DO (mg L−1), turbidity and average water depth of the site. Fish densities were then recorded across a 100 m transect using a standardized approach of dipnetting once per meter. A total of 30 individuals per species were taken from the Belgian sites, and 25 individuals per species from the Dutch sites. In the laboratory, fish from each site and species were kept separately for up to 48 h in aerated aquaria. Lighting was controlled under a 12:12 day:night regime. The fish were euthanized with a lethal dose of tricaine methanesulphonate (MS-222) and weight (g ± 0.001), length (cm ± 0.1 cm) and sex were recorded. A detailed evaluation of parasite infection was conducted for each individual, detecting the presence of parasites using a protocol developed for sticklebacks (Kalbe et al., Reference Kalbe, Wegner and Reusch2002). First, ectoparasites were identified scanning the entire body using a stereomicroscope. Then, endoparasites were identified using a microscope and high-pressure compressorium containing the liver, gut, intestines, swim bladder, body and head, kidneys, gill arches and right eye. Finally, a muscle sample without skin was removed from the right flank of each individual, and frozen for SIA. The muscle was chosen due to the lower stable isotope turnover rate when compared with other tissues such as the liver (Guelinckx et al., Reference Guelinckx, Maes, Geysen and Ollevier2008).

Fig. 1. Map of study area. A total of 9 sites were sampled across a salinity gradient in Belgium and the Netherlands. Sites in green are considered brackish (conductivity ⩾ 1000 μS cm−1), and sites highlighted in red are freshwater. Map was produced using QGIS (QGIS.org, Reference QGIS.org2021).

Stable isotope analysis

The energy flow among trophic links can be estimated through the use of stable isotopes. In particular, nitrogen stable isotopes (δ 15N) are typically enriched by 3–4‰ between consumers and their prey (Post, Reference Post2002). Carbon stable isotopes (δ 13C) on the other hand inform on the initial source of carbon (i.e. littoral vs pelagic), enabling inferences on the feeding ecology of organisms (Post, Reference Post2002). Hence, with SIA it is possible to determine the trophic niche of individuals within and among populations and explore the link between trophic niche and parasite infection. Specifically, the study of stable isotopes and parasite infection in populations of sympatric host species can reveal whether shared parasite infections consistently alter a host's trophic niche. Here, muscle samples were dried at 60°C for 48 h, ground into a fine powder using a ball mill (Retsch UK Ltd., Hope, UK), and weighed into tin capsules (~0.8 mg; Elemental Microanalysis Ltd., Oakhampton, UK) using an ultra-microbalance (±0.001 mg; Sartorius Lab Instruments, Göttingen, Germany). Isotope ratios were analysed using continuous flow isotope mass spectrometry (Sercon, Crewe, UK). Isotope ratios are expressed in per mil (‰) relative to known reference materials for both carbon (δ 13C) and nitrogen (δ 15N). The C:N ratio is a proxy of body composition with respect to lipids, and in the data it indicates an overall low lipid content, negating the need to apply lipid correction protocols (Post et al., Reference Post, Layman, Arrington, Takimoto, Quattrochi and Montaña2007).

Parasite infection and the environment

Environmental conditions affect the distribution of parasites. Hence, the variation in parasite abundance among host populations and the relationship between parasite community composition and abiotic environmental variables was investigated. To conduct the analyses effectively, 2 datasets were created. The first dataset used square-root-transformed aggregated mean parasite abundance for each site and species (hereafter aggregated dataset), and the second was kept as an individual-based dataset which permitted site-specific inferences. The aggregated dataset was created to avoid pseudo-replication when comparing parasite community data to environmental variables, upon which a principal coordinates analysis (PCoA) was performed based on Bray–Curtis dissimilarities. Then a permutational multivariate analysis of variance (PERMANOVA) on the same Bray–Curtis dissimilarity matrix was performed to assess the effect of site, host species and the interaction between host species and site using the adonis function in the R package vegan v.2.5.5. A post hoc test was run using the pairwise.adonis function (pairwiseAdonis package v.0.4) to identify all significant pairwise combinations (Martinez Arbizu, Reference Martinez Arbizu2020). Next, a similarity percentages test based on permutations (SIMPER) was used to identify which parasites contributed most to the differences between species and site. Finally, differences were correlated among environmental conditions (i.e. DO, pH, temperature, conductivity and total stickleback density), and the parasite community that used the PCoA axes' scores as response variables in linear models.

Linking parasite infection and trophic niche

An index of parasite load, the Individual Parasite Index (IPI; Kalbe et al., Reference Kalbe, Wegner and Reusch2002), was calculated for each individual fish comparing intra-individual parasite abundance to abundances. Two datasets were created. First, all samples were pooled and IPI was calculated across all sites (IPIAll), which enables comparison of infection load across fish from all sites. However, to detect patterns of trophic niche specialization across the entire dataset, the impact of parasite infections needs to be consistent and strong. Differences in parasite community, coevolutionary histories and environmental heterogeneity all decrease the likelihood of detecting such patterns in natural systems. Hence, the second dataset was created by calculating IPI within each site separately (IPISite), permitting the investigation of site-specific patterns of parasite infection that would be obscured using the whole dataset. This is needed as IPI considers variation in parasite count in the total sample.

To link parasite infection and trophic niche, the stable isotopes δ 15N and δ 13C of both species were used to investigate the hypothesis that parasite infection affects host trophic niche specialization. Firstly, collinearity was removed from IPIAll and fish length by taking the residuals of a linear model between those 2 factors. Models were then created using δ 15N as the response variable, an interaction between IPIAll, stickleback species and fish length as the fixed factor, and sampling site as a random factor. To permit comparisons among fish sizes, we grouped fish by species and size evenly into 3 groups based on length for each site. The same model was created but changing the response variable to δ 13C. Stepwise model selections were performed based on the AIC criterion using the step function in the lmerTest package v.3.1 (Kuznetsova et al., Reference Kuznetsova, Brockhoff and Christensen2017). In addition, previous work has shown that patterns of parasite-mediated trophic niche specialization can be achieved by a single parasite species (Pegg et al., Reference Pegg, Andreou, Williams and Britton2015a, Reference Pegg, Andreou, Williams and Britton2017). To investigate whether a single parasite putatively induced trophic niche specialization here, IPIAll was swapped for the square-root-transformed abundance of the 3 most important parasites as indicated by the SIMPER analysis (Gyrodactylus sp., Neochinorhynchus sp. and Cyathocotyle sp.), and Glochidia in L05 due to its extreme abundance only in this population. A minimum of 10 individual parasites from each of the 3 most common species were required to run the model.

Analyses for each site were repeated using the same model, only swapping IPIAll for each of the following terms individually: IPISite, and the square-root-transformed abundances of Gyrodactylus sp., Neochinorhynchus sp., Cyathocotyle sp. and Glochidia in L05. Linear mixed-effects models were obtained with the R package lme4 v.1.1 (Bates et al., Reference Bates, Mächler, Bolker and Walker2015), and P values were calculated using Satterthwaite's type II degree of freedom method from the lmerTest package. As each parasite load or abundance fixed factor addressed a discrete hypothesis, we corrected P values for multiple testing for each hypothesis. For instance, does δ 15N and δ 13C differ among Cyathocotyle sp. infections groups? As a result, each hypothesis had a maximum of 18 linear models, and P values were corrected for multiple testing accordingly using false discovery rate (FDR). However, for Glochidia, only 1 site (L05) had extreme abundances, so P values were not corrected when investigating these patterns.

To test whether parasite infection can specifically mediate host trophic niche constriction or divergence, for each site, fish were grouped by species, and then by parasite infection categories (low, medium and high). The infection categories were created by grouping individuals from each species evenly into 3 groups based upon IPISite. The groupings were repeated for each of the 3 most important parasite species indicated by the SIMPER analysis and Glochidia in site L05 due to their extreme abundance only in this site. The Kolmogorov–Smirnov test was used to determine whether the low and high infection group distributions differed among the δ 15N and δ 13C axes within each site. Where appropriate, P values were corrected for multiple testing using FDR. Given the role of stickleback species denoted by previous analyses, this analysis was repeated for all sites for each host species separately. All analyses were conducted in R v4.0.2 (R Development Core Team, 2019).

Results

Parasite infection and the environment

In total, 19 species of parasites were identified, with Gyrodactylus sp. and Cyathocotyle sp. species being identified in every site (Supplementary Table S1). Notably, Proteocephalus sp., Neochinorhynchus sp., Diplostomum sp. and Anguillicoloides crassus were all identified in at least 7 of 9 sites. Raphidascaris acus, Schistocephalus solidus and Paradilepsis scolecina were the rarest parasites in G. aculeatus, with the latter being found in 1 site only (Supplementary Table S1); whereas only Schistocephalus pungitii was restricted to P. pungitius. Parasite diversity significantly differed among species and site, and in general, G. aculeatus were infected with a more diverse assemblage of parasites [analysis of variance (ANOVA) parasite diversity, F 17,480 = 16.72, P < 0.001; Table 1], with mean individual parasite diversity differing among species and site, ranging from 1.50 ± 0.19 (standard error; s.e.) for P. pungitius in site L02 to 3.61 ± 0.25 for G. aculeatus in site L05. Similarly, parasite load (IPIAll) differed among sites and species, with G. aculeatus generally having higher infection loads than P. pungitius (ANOVA, F 17,480 = 14.41, P < 0.001).

Table 1. Summary of parasite diversity and load by site

Overall, parasite community composition significantly differed among sites and species (PERMANOVA, individual-based dataset, site:species, F 8,475 = 7.36, P < 0.001; aggregated mean parasite abundance dataset, site, F 8,17 = 3.15, P < 0.001, species, F 1,17 = 2.76, P = 0.040; pairwise PERMANOVAs are reported in Supplementary Table S2; Fig. 2A). Notably, the parasite community composition among species and site, inferred from PCoA axes (PC1 and PC2), was linked to DO in both stickleback species (3S, F 1,7 = 12.03, P = 0.010; 9S, F 1,7 = 15.07, P = 0.006, Fig. 2B). However, when investigating whether parasite diversity or parasite load was correlated with DO, we found no significant correlations in either G. aculeatus (parasite diversity, F 1,7 = 1.22, P = 0.117; IPIAll, F 1,7 = 2.60, P = 0.151) or P. pungitius (parasite diversity, F 1,7 = 2.57, P = 0.153; IPIAll, F 1,7 = 2.51, P = 0.157).

Fig. 2. Linking parasite community composition and the environment. (A) PCoA of mean aggregated parasite communities for each sampling site and species based on Bray–Curtis dissimilarity. Figures (B–C) represent the relationships between parasite community composition (PC2) and the DO for each host (B) G. aculeatus, and (C) P. pungitius. Colours denote sampling site. Squares represent G. aculeatus, triangles represent P. pungitius. Grey shaded area represents the 95% confidence intervals.

Based on site by species comparisons in the individual-based dataset, SIMPER analyses revealed that the abundances of 3 parasites, Gyrodactylus sp., Neochinorhynchus sp. and Cyathocotyle sp., were significantly different in at least 25% of all comparisons (Supplementary Table S3). Next, using the aggregated dataset, Diplostomum sp. and S. solidus were significantly more abundant in G. aculeatus than in P. pungitius (Supplementary Table S3). Among the 2 stickleback species, both Gyrodactylus sp. and Neochinorhynchus sp. were generally more abundant in G. aculeatus, whereas Cyathocotyle sp. was more abundant in P. pungitius (ANOVA; Gyrodactylus; F 17,480 = 26.80, P < 0.001; Neochinorhynchus; F 17,480 = 17.17, P < 0.001; Cyathocotyle; F 17,480 = 5.84, P < 0.001; Supplementary Table S1).

Linking parasite infection and trophic niche

Across all sites and after model simplification, the best predictor for δ 15N was stickleback species, whereby P. pungitius had higher δ 15N mean than G. aculeatus (LMER, random effect variance 4.05 ± 2.01, residual variance 2.91 ± 1.70, stickleback species, F 488 = 24.09, P < 0.001, Fig. 3A). The same stickleback species effect was detected when replacing IPIAll by the most common parasites, Gyrodactylus sp., Neochinorhynchus sp. and Cyathocotyle sp. Alternatively, δ 13C was positively correlated with an interaction between fish species and size, with the slope of G. aculeatus being steeper than that of P. pungitius (LMER, random effect variance 5.73 ± 2.39, residual variance 1.24 ± 1.12, species:length, F 487 = 8.48, P = 0.004, Fig. 3B). The same minimum adequate model was reached when replacing IPIAll with Gyrodactylus sp. or Neochinorhynchus sp. abundances. However, in addition to the species:length term, δ 13C increased with Cyathocotyle abundance in G. aculeatus, but decreased with this parasite in P. pungitius (LMER, random effect variance 5.58 ± 2.36, residual variance 1.23 ± 1.11, species:length, F 485 = 10.61, P = 0.001, Cyathocotyle:species, F 485 = 6.90, P = 0.009, Fig. 3B and C).

Fig. 3. Relationship of (A) δ 15N and (B–C) δ 13C and their significant predictors across all sites. (A) δ 15N by site residuals were plotted to remove the effect of site. (B–C) The solid lines represent the significant relationship between δ 13C, stickleback species and (B) fish length or (C) Cyathocotyle sp. abundance. The colours represent the same species in all panels, squares represent G. aculeatus and triangles P. pungitius. Grey shaded area represents the 95% confidence intervals.

To reiterate, it was hypothesized that by mixing multiple distinct population together it would obscure any context-dependent patterns, which merited investigating sites individually. Using the site-specific parasite index (IPISite), in 6 of 9 sites stickleback species was the most common predictor of both δ 15N (L01, L05, L07, L14) and δ 13C variation (L07, L10, L12; all models reported in Supplementary Table S4), although fish length or a species by length interaction was also significant in some sites (δ 15N; L05* and L14+; δ 13C; L12*; * indicates a species by length interaction; + indicates length as an additive effect; Supplementary Table S4). Notably, in 2 sites δ 15N increases with either IPISite or the abundance of Glochidia (LM; L03, IPISite, F 2,46 = 4.95, FDR = 0.025; L05, Glochidia, F 1,56 = 6.18, P = 0.016; Supplementary Table S4; Fig. 4A and B). When testing δ 15N variation against Cyathocotyle sp. abundance, a host species and parasite abundance interaction was observed in 1 site, with the correlation for P. pungitius being positive and the correlation for G. aculeatus being negative (L07, F 1,54 = 8.47, FDR < 0.001; Fig. 4C). Furthermore, stable isotope variation was explained by the interaction between parasite load (IPISite) and fish length in 3 sites. Firstly, in site L02 the correlation between δ 15N and IPISite was positive for large fish and negative for small fish (L02, F 1,45 = 4.23, FDR = 0.009; Fig. 4D). Secondly, in site L14 the correlation between δ 15N and Gyrodactylus sp. was negative for large fish and positive for small fish (L14, F 4,54 = 7.55, FDR = 0.001; Fig. 4E). Thirdly, in site L10 the correlation between δ 13C and Cyathocotyle sp. abundance was positive for large fish and negative for small fish (F 4,49 = 5.40, FDR = 0.004, Fig. 4F).

Fig. 4. Site-specific parasite infection relationship with (A–E) δ 15N and (F) δ 13C. Each plot shows the significant parasite infection term in the linear mixed models (cf. Results and Supplementary Table S4). Figures (A–F) correspond to sampling sites L03, L05, L07, L02, L14 and L10, respectively. Circles are used when the depicted effect is across both species (A–B, D–F). When the species effect was significant, red squares were used for G. aculeatus and blue triangles for P. pungitius (D). When appropriate, to better visualize the parasite infection by host length interaction (D–F), hosts were grouped into 3 approximately equal categories: small (purple), medium (dark blue) and large (green). Where necessary (B, C, E), parasite abundances were corrected for collinearity (cf. Methods). Grey shaded area represents the 95% confidence intervals.

Trophic niche constriction and divergence

There were no significant differences in trophic niche distributions when comparing low and high infection groups in each site (Supplementary Table S5, Fig. S1). However, when scrutinizing sites where specific parasites were linked to trophic niche (namely, L05δ 15N-Glochidia, L07δ 15N-Cyathocotyle sp., L10δ 13C-Cyathocotyle sp. and L14δ 15N-Gyrodactylus sp.), it was found that fish in site L05 with high infection of Glochidia had a significantly reduced δ 15N breadth compared to fish with low infection (2-sample Kolmogorov–Smirnov test, D = 0.5383, P = 0.003; Fig. 5). All other trophic niche distributions were not significantly different (Supplementary Table S5). Yet, when considering trophic niche use as a character, evidence of character displacement among the stickleback species was observed in trophic niche axes in 6 of 9 sites (Supplementary Table S5). Specifically, constricted trophic niches with G. aculeatus having the broader niche were observed in 2 sites (2-sample Kolmogorov–Smirnov test; L02, Dδ 15N = 0.408, FDRδ 15N = 0.045; L07, Dδ 15N = 0.483, FDRδ 15N = 0.007, Dδ 13C = 0.690, FDRδ 13C < 0.001; Supplementary Fig. S2) and P. pungitius with the broader niche also in 2 sites (L05, Dδ 15N = 0.434, FDRδ 15N = 0.012; L07, Dδ 13C = 0.469, FDRδ 13C = 0.007; L09, Dδ 15N = 0.425, FDRδ 15N = 0.019; Supplementary Fig. S2). Finally, partial trophic niche divergence was observed in 2 sites (L10, Dδ 15N = 0.481, FDRδ 15N = 0.0100, Dδ 13C = 0.630, FDRδ 13C < 0.001; L14, Dδ 15N = 0.464, FDRδ 15N = 0.007; Supplementary Fig. S2).

Fig. 5. Differences in trophic niche among Glochidia infection groups in sampling site L05. Fish were separated into approximately equal groups based on Glochidia abundance. Density of each infection group is plotted along the side of each axis. Squares represent G. aculeatus and triangles P. pungitius.

Discussion

Focusing on 2 coexisting and phylogenetically related species and their parasites, an environmental determinant of parasite infection was identified and to what extent parasites affect trophic niche specialization was tested. Firstly, parasite community composition was correlated with DO. Secondly, parasite load and the abundance of specific parasite species were regularly negatively correlated with hosts' stable isotope values (δ 15N or δ 13C). These patterns were often site-specific, highlighting the context-dependent nature of parasite-mediated selection on niche use. Of the predicted patterns of parasite-mediated trophic niche specialization (Britton and Andreou, Reference Britton and Andreou2016), one case of trophic niche constriction was observed, where the trophic niche of the highly infected fish was entirely nested within trophic niche of the lowly infected fish. Thirdly, both species differed in trophic niche use at the majority of the study sites, albeit not consistently or in the same direction. Together, the results suggest that both host species and parasite infection simultaneously affect niche specialization and the directions of these effects are context-dependent.

Parasite community composition and the environment

Many parasite species spend the vast majority of their lifetime, and for some species their entire lifetime, inside their hosts. Therefore, whether environmental conditions directly affect parasite abundance and diversity is still debated (e.g. Sweeny et al., Reference Sweeny, Albery, Venkatesan, Fenton and Pedersen2021). Here, a correlation between parasite community composition and DO was observed in both species. Such correlations can independently or jointly shape host–parasite interactions (Johnson et al., Reference Johnson, Chase, Dosch, Hartson, Gross, Larson, Sutherland and Carpenter2007; Brunner et al., Reference Brunner, Anaya-Rojas, Matthews and Eizaguirre2017; Ojwala et al., Reference Ojwala, Otachi and Kitaka2018). Indeed, variation in assemblages of parasites has been shown to correlate with DO in other aquatic systems (Lacerda et al., Reference Lacerda, Roumbedakis, Bereta Junior, Nuñer, Petrucio and Martins2018; Ojwala et al., Reference Ojwala, Otachi and Kitaka2018; Abdel-Tawwab et al., Reference Abdel-Tawwab, Monier, Hoseinifar and Faggio2019). However, when investigating whether specific measures of parasite community composition – namely parasite diversity and load – were also correlated with DO (and the other environmental variables: pH, temperature, conductivity and total stickleback density) we found no significant correlations. Such a result highlights the complex nature of the interaction between parasite communities and their environments over a spatial scale.

Linking parasite infection and trophic niche

The effects of parasite infection on host physiology and behaviour are well documented (Barber and Dingemanse, Reference Barber and Dingemanse2010; Abdel-Tawwab et al., Reference Abdel-Tawwab, Monier, Hoseinifar and Faggio2019). However, relatively little is known about the role of parasite infection on host trophic niche use in nature (Britton and Andreou, Reference Britton and Andreou2016; Brunner et al., Reference Brunner, Anaya-Rojas, Matthews and Eizaguirre2017). There are a few empirical examples of an association between parasite infection and host trophic niche constriction and divergence (Pegg et al., Reference Pegg, Andreou, Williams and Britton2015a, Reference Pegg, Andreou, Williams and Britton2017; Villalobos et al., Reference Villalobos, Nava-Bolaños, De Fuentes-Vicente, Téllez-Rendón, Huerta, Martínez-Hernández, Rocha-Ortega, Gutiérrez-Cabrera, Ibarra-Cerdeña and Córdoba-Aguilar2019), but whether environmental conditions, such as abiotic factors, interspecific competition or host density, affect these patterns remain elusive. Across the entire dataset, a correlation between Cyathocotyle sp. abundance and δ 13C was observed, with a host-specific direction of this correlation being observed. Cyathocotyle sp. is a trematode parasite that uses snails as intermediate hosts before actively infecting fish hosts, including both stickleback species (Lenz et al., Reference Lenz, Eizaguirre, Kalbe and Milinski2013; Kvach et al., Reference Kvach, Ondračková and Jurajda2016). Trematode parasites are known to affect host locomotion and trigger host immune responses (Khan et al., Reference Khan, Spiers and Pung2003; Goodman and Johnson, Reference Goodman and Johnson2011), including in fish hosts (Barber et al., Reference Barber, Hoare and Krause2000; Binning et al., Reference Binning, Shaw and Roche2017). Such effects are consistent with species-specific shifts in δ 13C, which suggest increasing Cyathocotyle sp. infection load is associated with change in host diet and may be a proxy or driver of displacement in niche use among host species. While both species share a similar landscape, the different responses to Cyathocotyle sp. infections likely stem from species-specific immune response and may be an example of parasite-mediated trophic niche divergence.

Individual trophic niche specialization within an ecosystem can be associated with differences in adaptive and behavioural traits such as body shape, body size, gape size, gill-raker morphology and habitat utilization, differentially exposing hosts to parasites (Cucherousset et al., Reference Cucherousset, Acou, Blanchet, Britton, Beaumont and Gozlan2011; Svanbäck and Schluter, Reference Svanbäck and Schluter2012; Pegg et al., Reference Pegg, Andreou, Williams and Britton2015b; Britton and Andreou, Reference Britton and Andreou2016; Barry et al., Reference Barry, Newton, Dodd, Evans, Newton and Adams2017). Recently, parasites have also been speculated as a driver of trophic niche specialization (Pegg et al., Reference Pegg, Andreou, Williams and Britton2015a, Reference Pegg, Andreou, Williams and Britton2015b, Reference Pegg, Andreou, Williams and Britton2017; Britton and Andreou, Reference Britton and Andreou2016). Here, parasite load – specifically the abundances of Gyrodactylus sp., Cyathocotyle sp. or Glochidia – was each associated with alterations of host δ 15N or δ 13C. Parasite infections can lead to significant phenotypic changes in their hosts (Miura et al., Reference Miura, Kuris, Torchin, Hechinger and Chiba2006; Britton and Andreou, Reference Britton and Andreou2016; Anaya-Rojas et al., Reference Anaya-Rojas, Best, Brunner, Eizaguirre, Leal, Melián, Seehausen and Matthews2019). Theory also predicts infected individuals to have different optimal diets due to differences in phenotypic capacity to detect, capture and digest the available prey (Araújo et al., Reference Araújo, Bolnick and Layman2011). Hence, the occasionally host size- or host species-dependent nature of the observations may be a consequence of the variation in parasite assemblages and the resultant parasite induced phenotypes within and among sites. Additionally, it should be noted that previous studies identifying consistent patterns of parasite-mediated trophic niche specialization targeted non-native parasites on native fishes (Pegg et al., Reference Pegg, Andreou, Williams and Britton2015a, Reference Pegg, Andreou, Williams and Britton2017). In the present study, the effects of both native and non-native parasites were included, where the ecological impact of each may not be in the same magnitude or direction.

Competition within and among species in an environment is another driver of individual niche specialization (Araújo et al., Reference Araújo, Bolnick and Layman2011; Evangelista et al., Reference Evangelista, Boiche, Lecerf and Cucherousset2014; Newsome et al., Reference Newsome, Tinker, Gill, Hoyt, Doroff, Nichol and Bodkin2015). Here, host species was the strongest and most consistent determinant of trophic niche. When considering trophic niche use as a character, these results are suggestive of interspecific competition and character displacement (Schluter and McPhail, Reference Schluter and McPhail1992; Gray and Robinson, Reference Gray and Robinson2002). Given the strength of species interaction, it is difficult to evaluate the role of parasite infection. As parasites contribute to host local adaptation and even to speciation (Eizaguirre et al., Reference Eizaguirre, Lenz, Traulsen and Milinski2009; Brunner and Eizaguirre, Reference Brunner and Eizaguirre2016), it seems likely they also influence these evolutionary mechanisms through processes like parasite spillover and parasite-mediated phenotypic modifications (Britton and Andreou, Reference Britton and Andreou2016). Finally, some of the patterns observed may also be summarized by differences in adaptive potential among host species, which is underpinned by differences in demographic histories (Raeymaekers et al., Reference Raeymaekers, Chaturvedi, Hablützel, Verdonck, Hellemans, Maes, De Meester and Volckaert2017; Bal et al., Reference Bal, Llanos-garrido, Chaturvedi, Verdonck, Hellemans and Raeymaekers2021). Specifically, G. aculeatus tends to rapidly adapt to changing local environmental conditions while P. pungitius instead has the capacity to tolerate a wide range of conditions (Raeymaekers et al., Reference Raeymaekers, Chaturvedi, Hablützel, Verdonck, Hellemans, Maes, De Meester and Volckaert2017; Bal et al., Reference Bal, Llanos-garrido, Chaturvedi, Verdonck, Hellemans and Raeymaekers2021).

Trophic niche constriction and divergence

After capturing the role of the environment and host species on parasite community and host trophic niche, the main hypotheses of this study can be tested: can parasites drive trophic niche specialization in a consistent and measurable manner (Britton and Andreou, Reference Britton and Andreou2016)? Here, the majority of the results showed context-dependent patterns of parasite-mediated shifts in trophic niche, if detectable at all. Notably, in 1 site (L05), parasite-mediated trophic niche constriction was observed, whereby the trophic niche of highly Cyathocotyle sp. infected fish was nested in that of lowly infected fish. A previous study showed similar patterns of trophic niche constriction suggesting it is likely due to infected conspecifics compensating the costs of infection by consuming specific food items already within the population's dietary breadth (Pegg et al., Reference Pegg, Andreou, Williams and Britton2017). Another factor to consider is differences in dietary preferences among life stages have also been shown to differentially expose hosts to parasites (Barry et al., Reference Barry, Newton, Dodd, Evans, Newton and Adams2017). However, given the field nature of the study, we are unable to identify whether the parasitized phenotype is the result of infection or whether existing dietary preferences differently exposed subpopulations to parasites (e.g. Cucherousset et al., Reference Cucherousset, Acou, Blanchet, Britton, Beaumont and Gozlan2011; Pegg et al., Reference Pegg, Andreou, Williams and Britton2015b; Barry et al., Reference Barry, Newton, Dodd, Evans, Newton and Adams2017). Regardless, one explanation for this observation is that parasite-mediated trophic niche constriction is driven by host diet compensation, and matches the only other known study to identify such a pattern (Pegg et al., Reference Pegg, Andreou, Williams and Britton2017). Overall, such findings highlight the need for further investigations into the causal mechanisms underlying this ecological phenomenon.

Conclusions

Overall, research on the determinants of trophic niche specialization has primarily identified factors such as competition, predation and resource availability as leading causes (Araújo et al., Reference Araújo, Bolnick and Layman2011). Whilst the results principally contribute to evidence supporting the impact of interspecific competition, the role of parasitism in trophic niche specialization remains poorly resolved (Britton and Andreou, Reference Britton and Andreou2016). Hence, studies observing an instance of trophic niche constriction, and multiple observations of parasitism putatively affecting the trophic niche of hosts are important. Such results suggest that consequences of parasite infection are an important driver of niche specialization which have been largely overlooked.

Supplementary material

The supplementary material for this article can be found at https://doi.org/10.1017/S0031182022000531.

Acknowledgements

We thank Martin Kalbe, Bart Hellemans, Conny Coeckelberghs, Ria Van Houdt, Scott Dorssers and Ian Sanders for field support and technical assistance, and the Royal Netherlands Institute of Sea Research (NIOZ) for hosting us in the Netherlands. Finally, we thank Liam Dickson for helping with the map.

Author contributions

D.-M. J. T., T. M. P. B., C. E., F. V. and J. A. M. R. conceived the study. D.-M. J. T., T. M. P. B., C. E., J. A. M. R. and I. S. D. collected the data. D.-M. J. T. analysed the data. J. A. M. R. and C. E. advised the analyses. D.-M. J. T., J. A. M. R. and C. E. drafted the manuscript. All co-authors contributed to the final version of the manuscript.

Financial support

D.-M. J. T. was supported by a studentship from Queen Mary University of London. C. E. was supported by research enabling fund from Queen Mary University of London. I. D. benefited from a scholarship of the Research Foundation-Flanders (FWO-SBO, 1S48117N). Fieldwork was supported by a Young Research Talents grant by Nord University to J. A. M. R., and T. M. P. B. was supported by Nord University.

Conflict of interest

None.

Ethical standards

Sampling was done in accordance with the European directive 2010/63/EU. Explicit permission for the Belgium and Dutch sites was granted by the Flemish Agency for Nature and Forests (Agentschap voor Natuur en Bos) and by Reptile, Amphibian & Fish Conservation Netherlands (RAVON; AVD261002016604), respectively. All individuals were euthanized following directions of the KU Leuven Animal Ethics Committee (https://admin.kuleuven.be/raden/en/animal-ethics-committee).

References

Abdel-Tawwab, M, Monier, MN, Hoseinifar, SH and Faggio, C (2019) Fish response to hypoxia stress: growth, physiological, and immunological biomarkers. Fish Physiology and Biochemistry 45, 9971013.CrossRefGoogle ScholarPubMed
Anaya-Rojas, JM, Best, RJ, Brunner, FS, Eizaguirre, C, Leal, MC, Melián, CJ, Seehausen, O and Matthews, B (2019) An experimental test of how parasites of predators can influence trophic cascades and ecosystem functioning. Ecology 100, 112.CrossRefGoogle ScholarPubMed
Araújo, MS, Bolnick, DI and Layman, CA (2011) The ecological causes of individual specialisation. Ecology Letters 14, 948958.CrossRefGoogle ScholarPubMed
Bal, TMP, Llanos-garrido, A, Chaturvedi, A, Verdonck, I, Hellemans, B and Raeymaekers, JAM (2021) Adaptive divergence under gene flow along an environmental gradient in two coexisting stickleback species. Genes 12, 115.CrossRefGoogle ScholarPubMed
Barber, I and Dingemanse, NJ (2010) Parasitism and the evolutionary ecology of animal personality. Philosophical Transactions of the Royal Society B: Biological Sciences 365, 40774088.CrossRefGoogle ScholarPubMed
Barber, I and Huntingford, FA (1995) The Effect of Schistocephalus solidus (Cestoda: Pseudophyllidea) on the Foraging and Shoaling Behaviour of Three-Spined Sticklebacks, Gasterosteus aculeatus. Behaviour, 12231240. doi: 10.1163/ej.9789004170292.i-540.37CrossRefGoogle Scholar
Barber, I, Hoare, D and Krause, J (2000) Effects of parasites on fish behaviour: A review and evolutionary perspective. Reviews in Fish Biology and Fisheries 10, 131165. doi: 10.1023/A:1016658224470CrossRefGoogle Scholar
Barry, J, Newton, M, Dodd, JA, Evans, D, Newton, J and Adams, CE (2017) The effect of foraging and ontogeny on the prevalence and intensity of the invasive parasite Anguillicola crassus in the European eel Anguilla anguilla. Journal of Fish Diseases 40, 12131222.CrossRefGoogle ScholarPubMed
Bates, D, Mächler, M, Bolker, BM and Walker, SC (2015) Fitting linear mixed-effects models using lme4. Journal of Statistical Software 67, 148.CrossRefGoogle Scholar
Binning, SA, Shaw, AK and Roche, DG (2017) Parasites and host performance: Incorporating infection into our understanding of animal movement. Integrative and Comparative Biology 57, 267280. doi: 10.1093/icb/icx024CrossRefGoogle ScholarPubMed
Blanar, CA, Munkittrick, KR, Houlahan, J, MacLatchy, DL and Marcogliese, DJ (2009) Pollution and parasitism in aquatic animals: a meta-analysis of effect size. Aquatic Toxicology 93, 1828.CrossRefGoogle ScholarPubMed
Britton, JR and Andreou, D (2016) Parasitism as a driver of trophic niche specialisation. Trends in Parasitology 32, 437445.CrossRefGoogle ScholarPubMed
Britton, JR, Pegg, J and Williams, CF (2011) Pathological and ecological host consequences of infection by an introduced fish parasite. PLoS ONE 6, 18.CrossRefGoogle ScholarPubMed
Brunner, FS and Eizaguirre, C (2016) Can environmental change affect host/parasite-mediated speciation? Zoology 119, 384394.CrossRefGoogle ScholarPubMed
Brunner, FS, Schmid-Hempel, P and Barribeau, SM (2014) Protein-poor diet reduces host-specific immune gene expression in Bombus terrestris. Proceedings of the Royal Society B: Biological Sciences 281, 110.Google ScholarPubMed
Brunner, FS, Anaya-Rojas, JM, Matthews, B and Eizaguirre, C (2017) Experimental evidence that parasites drive eco-evolutionary feedbacks. Proceedings of the National Academy of Sciences of the USA 114, 36783683.CrossRefGoogle ScholarPubMed
Copp, GH and Kováč, V (2003) Sympatry between threespine Gasterosteous aculeatus and ninespine Pungitius pungitius English lowland streams. 40, 341355.Google Scholar
Cucherousset, J, Acou, A, Blanchet, S, Britton, JR, Beaumont, WRC and Gozlan, RE (2011) Fitness consequences of individual specialisation in resource use and trophic morphology in European eels. Oecologia 167, 7584.CrossRefGoogle ScholarPubMed
Deflem, IS, Bennetsen, E, Opedal, ØH, Calboli, FCF, Ovaskainen, O, Van Thuyne, G, Volckaert, FAM and Raeymaekers, JAM (2021) Predicting fish community responses to environmental policy targets. Biodiversity and Conservation 30, 14571478. doi: 10.1007/s10531-021-02154-2CrossRefGoogle Scholar
Dobson, A, Lafferty, KD, Kuris, AM, Hechinger, RF and Jetz, W (2009) Homage to Linnaeus: how many parasites? How many hosts? In the Light of Evolution 2, 6382.Google Scholar
Eizaguirre, C, Lenz, TL, Traulsen, A and Milinski, M (2009) Speciation accelerated and stabilized by pleiotropic major histocompatibility complex immunogenes. Ecology Letters 12, 512. doi: 10.1111/j.1461-0248.2008.01247.xCrossRefGoogle ScholarPubMed
Evangelista, C, Boiche, A, Lecerf, A and Cucherousset, J (2014) Ecological opportunities and intraspecific competition alter trophic niche specialization in an opportunistic stream predator. Journal of Animal Ecology 83, 10251034.CrossRefGoogle Scholar
Goodman, BA and Johnson, PTJ (2011) Disease and the extended phenotype: Parasites control host performance and survival through induced changes in body plan. PLoS ONE 6, 110. doi: 10.1371/journal.pone.0020193CrossRefGoogle ScholarPubMed
Gray, SM and Robinson, BW (2002) Experimental evidence that competition between stickleback species favours adaptive character divergence. Ecology Letters 5, 264272.CrossRefGoogle Scholar
Guelinckx, J, Maes, J, Geysen, B and Ollevier, F (2008) Estuarine recruitment of a marine goby reconstructed with an isotopic clock. Oecologia 157, 4152.CrossRefGoogle ScholarPubMed
Hammerschmidt, K, Koch, K, Milinski, M, Chubb, JC and Parker, GA (2009) When to go: optimization of host switching in parasites with complex life cycles. Evolution 63, 19761986.CrossRefGoogle ScholarPubMed
Henrich, T, Benesh, DP and Kalbe, M (2013) Hybridization between two cestode species and its consequences for intermediate host range. Parasites and Vectors 6, 19.CrossRefGoogle ScholarPubMed
Johnson, PTJ, Chase, JM, Dosch, KL, Hartson, RB, Gross, JA, Larson, DJ, Sutherland, DR and Carpenter, SR (2007) Aquatic eutrophication promotes pathogenic infection in amphibians. Proceedings of the National Academy of Sciences of the USA 104, 1578115786.CrossRefGoogle ScholarPubMed
Kalbe, M, Wegner, KM and Reusch, TBH (2002) Dispersion patterns of parasites in 0+ year three-spined sticklebacks: a cross population comparison. Journal of Fish Biology 60, 15291542.CrossRefGoogle Scholar
Karvonen, A, Kristjánsson, BK, Skúlason, S, Lanki, M, Rellstab, C and Jokela, J (2013) Water temperature, not fish morph, determines parasite infections of sympatric Icelandic threespine sticklebacks (Gasterosteus aculeatus). Ecology and Evolution 3, 15071517.CrossRefGoogle Scholar
Khan, RN, Spiers, JA and Pung, OJ (2003) Effects of the trematode Microphallus turgidus on locomotion and prey capture in the grass shrimp Palaemonetes pugio. Journal of Helminthology 77, 327330. doi: 10.1079/joh2003192CrossRefGoogle ScholarPubMed
Kuris, AM, Hechinger, RF, Shaw, JC, Whitney, KL, Aguirre-Macedo, L, Boch, CA, Dobson, AP, Dunham, EJ, Fredensborg, BL, Huspeni, TC, Lorda, J, Mababa, L, Mancini, FT, Mora, AB, Pickering, M, Talhouk, NL, Torchin, ME and Lafferty, KD (2008) Ecosystem energetic implications of parasite and free-living biomass in three estuaries. Nature 454, 515518.CrossRefGoogle ScholarPubMed
Kuznetsova, A, Brockhoff, PB and Christensen, RHB (2017) lmerTest package: tests in linear mixed effects models. Journal of Statistical Software 82, 126.CrossRefGoogle Scholar
Kvach, Y, Ondračková, M and Jurajda, P (2016) First report of metacercariae of Cyathocotyle prussica parasitising a fish host in the Czech Republic, Central Europe. Helminthologia 53, 257261.CrossRefGoogle Scholar
Lacerda, ACF, Roumbedakis, K, Bereta Junior, JGS, Nuñer, APO, Petrucio, MM and Martins, ML (2018) Fish parasites as indicators of organic pollution in southern Brazil. Journal of Helminthology 92, 322331.CrossRefGoogle ScholarPubMed
Lafferty, KD and Kuris, AM (1999) How environmental stress affects the impacts of parasites. Limnology and Oceanography 44, 925931.CrossRefGoogle Scholar
Lafferty, KD, Dobson, AP and Kuris, AM (2006) Parasites dominate food web links. Proceedings of the National Academy of Sciences of the USA 103, 1121111216.CrossRefGoogle ScholarPubMed
Lafferty, KD, Allesina, S, Arim, M, Briggs, CJ, De Leo, G, Dobson, AP, Dunne, JA, Johnson, PTJ, Kuris, AM, Marcogliese, DJ, Martinez, ND, Memmott, J, Marquet, PA, McLaughlin, JP, Mordecai, EA, Pascual, M, Poulin, R and Thieltges, DW (2008) Parasites in food webs: the ultimate missing links. Ecology Letters 11, 533546.CrossRefGoogle ScholarPubMed
Lee, KP, Cory, JS, Wilson, K, Raubenheimer, D and Simpson, SJ (2006) Flexible diet choice offsets protein costs of pathogen resistance in a caterpillar. Proceedings of the Royal Society B: Biological Sciences 273, 823829.CrossRefGoogle Scholar
Lefèvre, T, Lebarbenchon, C, Gauthier-Clerc, M, Missé, D, Poulin, R and Thomas, F (2009) The ecological significance of manipulative parasites. Trends in Ecology and Evolution 24, 4148.CrossRefGoogle ScholarPubMed
Lenz, TL, Eizaguirre, C, Kalbe, M and Milinski, M (2013) Evaluating patterns of convergent evolution and trans-species polymorphism at mhc immunogenes in two sympatric stickleback species. Evolution 67, 24002412.CrossRefGoogle ScholarPubMed
Lockley, EC, Fouda, L, Correia, SM, Taxonera, A, Nash, LN, Fairweather, K, Reischig, T, Durão, J, Dinis, H, Roque, SM, Lomba, JP, Passos, L, Cameron, SJK, Stiebens, VA and Eizaguirre, C (2020) Long-term survey of sea turtles (Caretta caretta) reveals correlations between parasite infection, feeding ecology, reproductive success and population dynamics. Scientific Reports 10, 111.CrossRefGoogle ScholarPubMed
Mahmud, MA, Bradley, JE and MacColl, ADC (2017) Abiotic environmental variation drives virulence evolution in a fish host–parasite geographic mosaic. Functional Ecology 31, 21382146.CrossRefGoogle Scholar
Martinez Arbizu, P (2020) pairwiseAdonis: pairwise multilevel comparison using adonis. R package version 0.4.Google Scholar
Milinski, M (1984) Parasites determine a predator's optimal feeding strategy. Behavioral Ecology and Sociobiology 15, 3537.CrossRefGoogle Scholar
Miura, O, Kuris, AM, Torchin, ME, Hechinger, RF and Chiba, S (2006) Parasites alter host phenotype and may create a new ecological niche for snail hosts. Proceedings of the Royal Society B: Biological Sciences 273, 13231328.CrossRefGoogle ScholarPubMed
Newsome, SD, Tinker, MT, Gill, VA, Hoyt, ZN, Doroff, A, Nichol, L and Bodkin, JL (2015) The interaction of intraspecific competition and habitat on individual diet specialization: a near range-wide examination of sea otters. Oecologia 178, 4559.CrossRefGoogle ScholarPubMed
Ojwala, RA, Otachi, EO and Kitaka, NK (2018) Effect of water quality on the parasite assemblages infecting Nile tilapia in selected fish farms in Nakuru County, Kenya. Parasitology Research 117, 34593471.CrossRefGoogle ScholarPubMed
Pagán, I, Alonso-Blanco, C and García-Arenal, F (2008) Host responses in life-history traits and tolerance to virus infection in Arabidopsis thaliana. PLoS Pathogens 4, 110.CrossRefGoogle ScholarPubMed
Pegg, J, Andreou, D, Williams, CF and Britton, JR (2015 a) Temporal changes in growth, condition and trophic niche in juvenile Cyprinus carpio infected with a non-native parasite. Parasitology 142, 15791587.CrossRefGoogle ScholarPubMed
Pegg, J, Andreou, D, Williams, CF and Britton, JR (2015 b) Head morphology and piscivory of European eels, Anguilla anguilla, predict their probability of infection by the invasive parasitic nematode Anguillicoloides crassus. Freshwater Biology 60, 19771987.CrossRefGoogle Scholar
Pegg, J, Andreou, D, Williams, CF and Britton, JR (2017) Consistent patterns of trophic niche specialization in host populations infected with a non-native copepod parasite. Parasitology 144, 945953.CrossRefGoogle ScholarPubMed
Post, DM (2002) Using stable isotopes to estimate trophic position: models, methods, and assumptions. Ecology 83, 703718.CrossRefGoogle Scholar
Post, DM, Layman, CA, Arrington, DA, Takimoto, G, Quattrochi, J and Montaña, CG (2007) Getting to the fat of the matter: models, methods and assumptions for dealing with lipids in stable isotope analyses. Oecologia 152, 179189.CrossRefGoogle Scholar
QGIS.org, 2021. QGIS Geographic Information System. QGIS Association. http://www.qgis.orgGoogle Scholar
Raeymaekers, JAM, Chaturvedi, A, Hablützel, PI, Verdonck, I, Hellemans, B, Maes, GE, De Meester, L and Volckaert, FAM (2017) Adaptive and non-adaptive divergence in a common landscape. Nature Communications 8, 18.CrossRefGoogle Scholar
R Development Core Team (2019) R: A Language and Environment for Statistical Computing. Vienna, Austria: R Foundation for Statistical Computing. Available at http://www.R-project.org/. doi:10.1007/978-3-540-74686-7.Google Scholar
Schluter, D and McPhail, JD (1992) Ecological character displacement and speciation in sticklebacks. American Naturalist 140, 85108.CrossRefGoogle ScholarPubMed
Svanbäck, R and Schluter, D (2012) Niche specialization influences adaptive phenotypic plasticity in the threespine stickleback. American Naturalist 180, 5059.CrossRefGoogle ScholarPubMed
Sweeny, AR, Albery, GF, Venkatesan, S, Fenton, A and Pedersen, AB (2021) Spatiotemporal variation in drivers of parasitism in a wild wood mouse population. Functional Ecology 35, 111.CrossRefGoogle Scholar
Varadharajan, S, Rastas, P, Löytynoja, A, Matschiner, M, Calboli, FCF, Guo, B, Nederbragt, AJ, Jakobsen, KS and Merilä, J (2019) A high-quality assembly of the nine-spined stickleback (Pungitius pungitius) genome. Genome Biology and Evolution 11, 32913308.Google ScholarPubMed
Villalobos, G, Nava-Bolaños, A, De Fuentes-Vicente, JA, Téllez-Rendón, JL, Huerta, H, Martínez-Hernández, F, Rocha-Ortega, M, Gutiérrez-Cabrera, AE, Ibarra-Cerdeña, CN and Córdoba-Aguilar, A (2019) A reduction in ecological niche for Trypanosoma cruzi-infected triatomine bugs. Parasites and Vectors 12, 19.CrossRefGoogle ScholarPubMed
Young, RE and MacColl, ADC (2017) Spatial and temporal variation in macroparasite communities of three-spined stickleback. Parasitology 144, 436449.CrossRefGoogle ScholarPubMed
Zander, CD (2007) Parasite diversity of sticklebacks from the Baltic Sea. Parasitology Research 100, 287297.CrossRefGoogle ScholarPubMed
Figure 0

Fig. 1. Map of study area. A total of 9 sites were sampled across a salinity gradient in Belgium and the Netherlands. Sites in green are considered brackish (conductivity ⩾ 1000 μS cm−1), and sites highlighted in red are freshwater. Map was produced using QGIS (QGIS.org, 2021).

Figure 1

Table 1. Summary of parasite diversity and load by site

Figure 2

Fig. 2. Linking parasite community composition and the environment. (A) PCoA of mean aggregated parasite communities for each sampling site and species based on Bray–Curtis dissimilarity. Figures (B–C) represent the relationships between parasite community composition (PC2) and the DO for each host (B) G. aculeatus, and (C) P. pungitius. Colours denote sampling site. Squares represent G. aculeatus, triangles represent P. pungitius. Grey shaded area represents the 95% confidence intervals.

Figure 3

Fig. 3. Relationship of (A) δ15N and (B–C) δ13C and their significant predictors across all sites. (A) δ15N by site residuals were plotted to remove the effect of site. (B–C) The solid lines represent the significant relationship between δ13C, stickleback species and (B) fish length or (C) Cyathocotyle sp. abundance. The colours represent the same species in all panels, squares represent G. aculeatus and triangles P. pungitius. Grey shaded area represents the 95% confidence intervals.

Figure 4

Fig. 4. Site-specific parasite infection relationship with (A–E) δ15N and (F) δ13C. Each plot shows the significant parasite infection term in the linear mixed models (cf. Results and Supplementary Table S4). Figures (A–F) correspond to sampling sites L03, L05, L07, L02, L14 and L10, respectively. Circles are used when the depicted effect is across both species (A–B, D–F). When the species effect was significant, red squares were used for G. aculeatus and blue triangles for P. pungitius (D). When appropriate, to better visualize the parasite infection by host length interaction (D–F), hosts were grouped into 3 approximately equal categories: small (purple), medium (dark blue) and large (green). Where necessary (B, C, E), parasite abundances were corrected for collinearity (cf. Methods). Grey shaded area represents the 95% confidence intervals.

Figure 5

Fig. 5. Differences in trophic niche among Glochidia infection groups in sampling site L05. Fish were separated into approximately equal groups based on Glochidia abundance. Density of each infection group is plotted along the side of each axis. Squares represent G. aculeatus and triangles P. pungitius.

Supplementary material: File

Thorburn et al. supplementary material

Tables S1-S5
Download Thorburn et al. supplementary material(File)
File 199.2 KB
Supplementary material: File

Thorburn et al. supplementary material

Figures S1-S2
Download Thorburn et al. supplementary material(File)
File 1 MB