Research Article |
Corresponding author: Kathleen D. W. Church ( kathleen.church@mail.mcgill.ca ) Academic editor: Cara Gallagher
© 2025 Kathleen D. W. Church, Steven F. Railsback, Christina A. D. Semeniuk.
This is an open access article distributed under the terms of the Creative Commons Attribution License (CC BY 4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Citation:
Church KDW, Railsback SF, Semeniuk CAD (2025) Influence of personality traits on the response of a modelled population of stream-dwelling rainbow trout (Oncorhynchus mykiss) to microplastics consumption. Individual-based Ecology 1: e137398. https://doi.org/10.3897/ibe.1.e137398
|
Microplastics in freshwater habitats are consumed by fish, including stream-dwelling salmonids, which can alter food consumption or negatively affect swimming and foraging behaviour. As population-level effects are largely unknown, a population of stream-dwelling rainbow trout (Oncorhynchus mykiss) was simulated using the agent-based model ‘inSTREAM 7’ to model population-level effects (biomass) of behavioural changes caused by microplastics consumption. Individual fish were assigned all possible combinations of two personality traits (dominance, boldness/shyness) and consumed microplastics while foraging; and their microplastics consumption, body length and abundance were tracked for three different life stages (fry, juvenile, adult) for a period of 10 years. Three scenarios were explored: a low-impact scenario with microplastics causing decreased food consumption, a medium-impact scenario which added lower swimming speed and a high-impact scenario with additional reductions of foraging efficiency. Each scenario was tested for microplastics concentrations of 0%, 1% (i.e. current levels) and 3% (i.e. future levels) of drift food. Overall, microplastics consumption did not strongly affect trout abundance. Dominant adult trout consumed disproportionally more microplastics than all other fish, especially at higher microplastic concentrations. Personality traits influenced the response of the trout to microplastics ingestion: dominant and bold adults were smaller when food consumption was reduced; shy and subordinate adults were smaller when swimming speed was lowered; and all dominant adults, regardless of boldness, were smaller when foraging efficiency was impeded, with dominant and bold fry also less abundant in this scenario. However, effects on fish body length were only found at microplastic concentrations of 3%, indicating these outcomes can be prevented, as current levels of microplastics pollution are below this concentration. Nevertheless, microplastics ingestion may become an additional stressor that interacts with the myriad of mostly anthropogenic stressors that already affect wild salmonid populations.
Agent-based models (ABMs), behaviour, microplastics, population abundance, salmonids
Microplastics, defined as plastic particles < 5 mm (i.e.
Microplastics have been found in the gut of numerous freshwater animals (
Previous research on the behavioural responses of individuals to microplastics has shown that fish that ingest microplastics often exhibit changes in their swimming behaviour. For example, zebrafish (Danio rerio) engaged in hyperactive swimming behaviour following microplastics ingestion (exposure concentration: 1 mg/l;
Mechanistic population modelling enables the scaling up of effects of microplastics consumption from the individual to the population-level. Agent Based Models (ABMs), in particular, are effective tools to evaluate the effects that emerge from behavioural changes at the individual level and upscale them to the population level (
In this study, an existing individual-based model of stream-dwelling trout population dynamics (
‘inSTREAM 7.3’: The ABM ‘inSTREAM 7.3’ (individual-based Stream Trout Environmental Assessment Model) is an individual-based model of stream trout, created for river management and research purposes to assess how changes in habitat characteristics (i.e. temperature, flow, turbidity, food production) and biological processes (i.e. competition, food intake, predation) affect trout populations (
The version of the model used in this study represents a single reach of a river (184 m long, ~ 5–10 m wide), consisting of a grid of rectangular cells that represents the spatial environment (i.e. the “mainstem” reach;
During a 10-year simulation run, daily physical inputs (stream flow, temperature and turbidity) influence trout survival and body size by affecting the habitats that individual fish select and the activities they choose to engage in, which then determines their foraging success and predation risk. The accumulated decisions of thousands of individual trout then impact population metrics like abundance, via the survival and reproductive success of numerous individual trout, as well as mean body length, through the body condition and growth of many individual trout (Fig.
Flow diagram showing the key processes of inSTREAM7 and how different factors affect trout microplastics consumption, survival, growth and abundance. Legend: Dotted lines indicate the physical environment, dashed lines show personality traits (dominance, boldness/shyness), underlined text indicates trout behavioural choices and bold text indicates model outcomes.
Simulations were run using the default, or ‘standard’ parameter values, except as noted. Standard parameter values were defined, based on a sensitivity analysis using the same simulated freshwater river habitat and were previously calibrated to obtain realistic long-term average results (as described by
To represent the effects of microplastic consumption on trout food consumption, swimming ability and foraging efficiency in our study, additional functions and modifications were made to the source code of ‘inStream 7.3’ using NetLogo (version 6.2.1;
Microplastics were added to both the water column and sediments in all three scenarios. The original “inSTREAM 7.3” included the concentrations of food available as either invertebrate drift-consumed while drift foraging or food consumed while search-feeding in the sediments, as the daily rate of available prey in grams. Microplastics were added as another type of food, represented by the parameter micro-conc and expressed as a percentage of the food available in the water column (i.e. drift food), with twice as many microplastics in the sediments (i.e. search food), in accordance with the observations in natural habitats (i.e.
In this scenario, fish are not able to swim as fast following microplastics consumption, which has been observed in laboratory studies (see below), but has not yet been evaluated under natural conditions; this is in addition to experiencing the reduced food consumption described for Scenario 1. Limited swimming ability can affect habitat choice and food intake, as fish become limited to foraging in low velocity habitats and are more likely to hide in shelter than to forage when water velocity is high. To account for these changes, a memory function, similar to the memory function used for food consumption, was added to keep track of the quantity of microplastics consumed by individual trout over the previous seven days. Consumption thresholds are based on laboratory findings from juvenile African freshwater catfish (Clarias gariepinus) exposed to microplastics and correspond to 5 mg (threshold 1) and 28 mg (threshold 2) of consumed microplastics for a mean body length of 17.28 cm (
For Scenario 3, prey capture efficiency through drift foraging is also reduced in high velocity water flow, in addition to the reduction in swimming speed and food consumption that occurs in both previous scenarios. In inSTREAM 7.3, a decreasing logistic function is used to represent foraging success as a function of the ratio of stream velocity over the maximum sustainable swimming speed. The ability to capture prey at higher velocities is represented by the parameter trout-capture-R1, the ratio of stream velocity over the maximum sustainable swimming speed at which foraging success is 10%. Reductions in foraging efficiency in high velocity habitats will further impact fitness through habitat choice and food intake, as fish become even more limited to lower velocity habitats and more likely to hide than to forage. This parameter is lowered from 1.3 to 1.17 when the first microplastics consumption threshold is reached (i.e. 1.5, 2.7 or 6.1 mg per week for fry, juveniles or adults, respectively) and is further lowered to 1.0 when the second consumption threshold is reached (i.e. 10, 15 or 33 mg per week for fry, juveniles or adults, respectively). No further decreases in prey capture efficiency occurred after the second threshold. This change in parameter values is also permanent, even though consumed microplastics are regularly excreted and do not build up in the bodies of the fish, given that microplastics are always available for consumption. The change represents a reduction of the ability of trout to capture any prey when the velocity of their habitat exceeds their maximum sustainable swimming speed (
We assigned two separate personality traits to each fish: dominance (i.e. increased access to food resources;
Values of model parameters associated with the dominance and shyness/boldness personality traits, based on default values suggested for inSTREAM 7.3 (
Model parameter | Personality trait | Levels | Parameter value |
---|---|---|---|
mort-terr-pred-hiding-factor | Shyness / Boldness | Shy | 0.90 |
Intermediate | 0.80 | ||
Bold | 0.70 | ||
mort-fish-pred-hiding-factor | Shyness / Boldness | Shy | 0.70 |
Intermediate | 0.60 | ||
Bold | 0.50 | ||
trout-react-dist-B | Dominance | Dominant | 2.25 |
Intermediate | 2.00 | ||
Subordinate | 1.75 |
Both personality traits were incorporated into all three scenarios. All trout are expected to consume greater quantities of microplastics with higher concentrations of microplastics; however, as dominant trout are defined by a higher foraging success, they are expected to consume proportionally greater quantities of microplastics, relative to intermediate and, especially, to subordinate trout. Dominant trout are expected to reach the microplastics consumption thresholds that result in permanently lowered swimming speed and reduced foraging efficiency more frequently than intermediate or subordinate trout, especially the second microplastic-consumption thresholds. Consequently, a greater difference in fitness is expected for dominant trout in the presence of microplastics relative to intermediate and subordinate trout, regardless of their boldness.
Additionally, the trade-off between risk and reward may break down for bold trout when exposed to higher concentrations of microplastics, with reductions in foraging success (driven by lower swimming speed and/or less efficient prey capture) accompanied by no changes in predation risk. As a result, trout with average boldness may experience higher fitness relative to bold or shy trout at higher microplastics’ concentrations. Shy individuals will spend less time foraging and spend more time hiding than trout with average boldness and especially bold trout, as a result of reduced swimming speed and impeded foraging efficiency. However, as trout that are both subordinate and shy tend to prioritise foraging less than the other combined personality traits, they are also expected to consume fewer microplastics and are, thus, less likely to be negatively harmed by microplastics ingestion, especially under Scenario 3. The same would apply to subordinate trout with average boldness, but to a lesser degree. Trout that are both dominant and bold, on the other hand, are more likely to be affected by microplastics as they tend to prioritise foraging more than the other combined personality traits and will, thus, be more likely to experience the negative effects of microplastic consumption relative to dominant trout that are also shy or trout with intermediate dominance, regardless of their level of boldness.
The model was evaluated to determine if it reproduced patterns observed in real-life fish regarding personality and microplastics consumption. The literature was first searched for general patterns driven by the personality traits and microplastic concentrations and their associated assumptions, that were added to the model. We then evaluated these assumptions by how well the model reproduced those patterns. Here, we present the patterns found and how these patterns were defined in a statistically testable way (detailed below).
To confirm the realism of our modelled personality traits, we used a pattern-oriented modelling approach to evaluate the recreation of four real-world behavioural patterns: i) variable survival pattern: survival for different personality traits varies over time (Dingemanse and Réale 2013), ii) behavioural syndrome under predation pattern: increased predation selects for positive correlations between boldness and dominance (i.e. behavioural syndromes; Bell and Sih (2007)), iii) risk-reward with abundant food pattern: the shy-bold behavioural axis results from trade-offs between risk and reward that tend to break down when food is abundant, as fish have access to plentiful food regardless of risk-taking (Wilson et al. 1994; Seltmann et al. 2014) and iv) dominance in a competitive environment pattern: higher competition environments select for more dominant fish (i.e. Gilmour et al. (2005)). The personality patterns were verified using Scenario 1 (see below) in the absence of microplastics (0%), using data only from the latter 8 years of the simulations to ensure the results were independent of the initial conditions.
For the variable survival pattern (i.e. personality pattern i), the coefficient of variation (CV = (SD/mean) × 100)) was used to quantify variability in survival over time or the proportion of total abundance for the nine different combinations of the two personality traits within each age group. Generally, a coefficient of variation above 20% indicates substantial variability (i.e. Reed et al. (2002)). As personality traits were assigned to trout in equal numbers, differences in abundance reflect differential survival. Welch’s t-tests were then used to assess the occurrence of the remaining three real-world behavioural patterns. For these patterns, one run of the 10-year simulation was used to represent each of the environmental conditions being compared (i.e. one simulation run for low predation vs. one simulation run for high predation), with data collected annually on the census day (i.e. Julian day: 266; only non-leap year values were used for simplicity) from the 3rd until the 10th year. For the behavioural syndrome under predation pattern (i.e. personality pattern ii), the proportional abundance of bold and dominant trout, shy and dominant trout and shy and subordinate trout were compared in a high and a low predation environment. These three combined personality traits were selected to facilitate comparison with the real-world study that reported this particular pattern (i.e. Bell and Sih (2007)). To compare high and low predation pressure, the minimum probability of surviving both terrestrial (i.e. reach-terr-pred-min) and fish predation (i.e. reach-fish-pred-min) in the riskiest combination of habitat and activity was lowered to 0.85 (i.e. high predation) and raised to 0.95 (i.e. low predation) from a standard value of 0.90 which represents a moderate quantity of predation for many stream environments (i.e.
A pattern-oriented modelling approach was also used to evaluate the occurrence of three real-life patterns of microplastic consumption that are affected by life history and optimal foraging: i) larger fish consume more microplastics than smaller fish (Horton et al. 2018; McIlwraith et al. 2021), ii) consumption of microplastics is higher when environmental concentrations are higher (
To determine whether larger fish consume more microplastics than smaller fish (i.e. microplastics consumption pattern i), a linear model was constructed with the quantity of microplastics consumed as the response variable and body size (small, large), microplastics concentration (1%, 3%) and age (0, 1, 2+) as the predictor variables. Fish of each age group were classified as small or large based on their body length relative to the mean body length for their corresponding age group. Small fish were shorter than one standard deviation below the mean body length for their age category, while large fish were longer than one standard deviation above the mean body length. Welch’s t-tests were then used to assess if higher microplastics consumption occurs when environmental concentrations are also higher (i.e. microplastics consumption pattern ii). Overall microplastics consumption of trout was compared with the lowest (i.e. 1%) and highest (i.e. 3%) microplastics concentrations, with data collected as above, with one simulation run for each of these two concentrations of microplastics. Finally, another linear model was constructed to examine if microplastics consumption is higher when food is more scarce (i.e. microplastics consumption pattern iii). Microplastics consumption, scaled to the total quantity of available food, was the response variable, while the quantity of available food (standard quantity, half the standard quantity), the concentration of microplastics and trout age were the predictor variables.
To evaluate the effects of the three different scenarios on fish biomass at the population level, data from several annual censuses conducted for each of the 10-year simulations were used, with only the last 8 years of the simulation considered to ensure results were independent from initial conditions (
Results from model simulations were analysed using linear mixed-effect models, using the ‘lme4’ package in R (
Mortality at age 1 and age 2 was determined by subtracting the abundance of fish from each age group at the end of the year (i.e. Julian Day 365), from the abundance at the beginning of the year (i.e. Julian Day 1, when all fish age by 1 year). As the simulation begins and ends in September, the first calculations for mortality used data from Julian Day 1 on the 3rd year and Julian Day 365 on the 4th year (i.e. year 4 mortality), while the final mortality calculations used data from Julian Day 1 on the 9th year and Julian Day 365 on the 10th year (i.e. year 10 mortality). The number of fish of each age class on Julian day 1 was used to determine recruitment. As personality traits are assigned to emerging fish in equal numbers in May and June, mortality for age 0 fish was not explicitly determined, but was rather inferred through fry abundance.
A sensitivity analysis (summarised in Suppl. material
Four patterns were used to verify whether our personality parameterisation reproduced known, real-world responses. For the variable survival pattern, both personality traits and all age groups exhibited substantial variability in survival over the simulation run (Suppl. material
In comparing model predictions of microplastics consumption with real-world patterns, larger fish consumed more microplastics than smaller fish within the same age category (LM, F2,288 = 26.73, p < 0.0001). As well, more microplastics were consumed by all fish when microplastics concentrations in the environment were higher (1% vs. 3% microplastics, Welch’s t-test, t = -3.78, df = 154.7, p-value = 0.0002). Relative microplastics consumption was also higher when food was more scarce, when microplastics consumption was scaled to account for food availability (LMM, F1.720 = 5.18, p = 0.023).
Microplastics consumption was also affected by significant interactions between age, personality and microplastics concentration in all three scenarios (LMMs, all p’s < 0.00001; Fig.
Mean weekly microplastics consumption (± 95% C.I.s) of age 2+ rainbow trout with different combinations of personality traits (shy, average, bold; dominant, intermediate, subordinate) at different microplastic (MP) concentrations (0%, 1%, 3%) for each scenario: A reduction in food consumption, B reduction in swimming speed and C reduction in foraging ability. Legend: sub = subordinate, avg = average-boldness, dom = dominant, int = intermediate-dominance.
As well, significant interactions were found amongst trout age, microplastics concentration and simulation year. A decrease in microplastics consumption over time was observed for age-2+ fish at concentrations of 3% under Scenario 1 (LMM, F4,5987 = 5.81, p = 0.00012; Fig.
Fish body size, as represented by body length, was significantly affected by interactions amongst age, personality and microplastics concentration, for all three scenarios (LMMs, all p’s < 0.00001; Fig.
Mean body length (± 95% C.I.s) of age 2+ rainbow trout with different combinations of personality traits (shy, average, bold; dominant, intermediate, subordinate) at different microplastics (MP) concentrations (0%, 1%, 3%), for each scenario: A reduction in food consumption, B reduction in swimming speed and C reduction in foraging ability. Legend: sub = subordinate, avg = average-boldness, dom = dominant, int = intermediate-dominance.
In Scenario 1, body length of dominant and bold age-2+ fish was negatively affected at microplastics concentrations of 3% (Tukey’s HSD, 1% vs. 3% microplastics, z = -3.36, p = 0.0022, 0% vs. 3% microplastics, z = -2.83, p = 0.013, Fig.
Interactions between age and simulation year were significant for all three scenarios (LMMs, all p’s < 0.00001), as body size of age-2+ fish decreased over time in a similar manner in each of the scenarios (Fig.
Patterns of recruitment and mortality appear to reflect general patterns of abundance, with higher mortality and higher survival experienced by fish that were more abundant (Suppl. material
Abundance was affected by interactions between age and personality for all three scenarios (LMMs, all p’s < 0.00001). Abundance differed amongst the different personality traits, but only for age-0 fish (Fig.
Abundance of age-0 rainbow trout (± 95% C.I.s) with different combinations of personality traits (shy, average, bold; dominant, intermediate, subordinate) at different microplastics (MP) concentrations (0%, 1%, 3%) for each scenario: A reduction in food consumption, B reduction in swimming speed and C reduction in foraging ability. Legend: sub = subordinate, avg = average-boldness, dom = dominant, int = intermediate-dominance.
Similarly, microplastics’ concentration significantly interacted with age in all three scenarios, but only affected abundance in age-0 fish (LMMs, all p’s < 0.01; Fig.
Interactions between trout age, microplastics concentration and year were significant in Scenario 1 (LMM, F4,6411 = 3.09, p = 0.015; Fig.
The effects of the tested parameters on population biomass are represented numerically by the magnitude of the scaled slope (summarised in Suppl. material
Microplastics consumption impacted trout abundance and fish body size through behavioural changes, namely, changes in swimming and foraging ability (i.e. prey capture efficiency), but only for certain personality traits. Overall, dominant trout were the most affected by microplastics consumption at all life stages, especially bold and dominant trout, as predicted. Only the body length of bold and dominant adults in Scenario 2 and dominant adults in Scenario 3 and the abundance of bold and dominant fry in Scenario 3, were impacted by microplastics consumption. These impacts were likely driven by the overall greater foraging success of trout with these personality traits, which results in the consumption of substantially greater quantities of microplastics, especially when higher concentrations of microplastics are present in the environment, according to predictions. However, the reduction in adult microplastics consumption in both Scenario 2 and Scenario 3 likely indicates a negative feedback loop, with the negative impacts of microplastics consumption on swimming and foraging behaviour also reducing the ability of adult fish to forage on microplastics.
However, dominance was not the only personality trait negatively affected by microplastics consumption. Contrary to predictions, the body size of adult fish that were both subordinate and shy was negatively affected by microplastics consumption, but only at the highest concentration of microplastics. This trend appears to be driven by fewer older and larger adult fish with these personality traits. These results suggest that both survival and reproductive success of shy and subordinate adults are likely to be negatively affected by high concentrations of microplastics, since larger fish are generally less vulnerable to gape-limited predators (
Our results suggest that individuals within a population of stream salmonids are not likely to suffer equally from microplastics ingestion. The much higher consumption of microplastics by dominant trout shows that the effects of microplastics consumption are not uniform across individuals. Similar to our findings, large individual differences in microplastics ingestion have been observed amongst juvenile spiny chromis damselfish (0–2000 particles per individual; Acanthochromis polyacanthus;
However, our model is limited by its lack of consideration of multiple interacting stressors, neurotoxic effects (i.e. increased susceptibility to predation, altered sensory perception, reduced food detection) and changes in personality, as personality traits were fixed in this study. Currently, the model is conservative, but can be further modified in future iterations to accommodate these factors. As well, future versions of the model would also benefit from salmonid-specific parameters on the effects of microplastics consumption on swimming and foraging ability to further improve realism, as the certainty of ABMs increases with greater data availability. Microplastics pollution represents an additional stressor for salmonid populations that are already under strain from the combined negative effects of climate change (
The results from the sensitivity analysis demonstrated that the impacts of microplastics consumption on trout populations can be determined via its effects on both swimming ability and the efficiency of prey capture, factors which largely determine foraging success for salmonids. For instance, the variable, “reaction distance for prey” represents the area over which fish capture food when drift feeding, a proxy of territory size (i.e.
More empirical research is needed to determine salmon-specific parameters of microplastics consumption and to identify the precise behavioural changes in swimming and foraging behaviour that may occur when different quantities of microplastics are consumed by stream-dwelling salmonids, such as changes in station-holding ability (i.e. resistance), reductions of swimming speed, alterations in foraging efficiency during drift foraging or increased reliance on benthic rather than drift foraging. These findings would be particularly insightful, given that other environmental factors, including sedimentation (i.e. Sweka et al. (2003); Zamor et al. (2007);
The effects of microplastics consumption can vary widely between species and were reflected in the three scenarios explored in this study. For most animals, microplastics consumption leads to a reduced intake of real food items only (
The ubiquitous presence of microplastics in aquatic habitats is generally viewed as a major threat to biodiversity (
The implications of the present study are that researchers studying and predicting the effects of microplastics should characterise fish behavioural types and incorporate these traits in their analyses to improve their understanding and predictive ability. This study found that microplastics consumption over 5–6 generations is not likely to pose a significant threat to healthy populations of stream-dwelling salmonids, like rainbow trout, via alterations to behaviour such as swimming and foraging. However, traits of dominant individuals of all ages such as body size will be disproportionally affected by the direct, negative effects of microplastics due to their much higher consumption rates. Microplastics consumption represents an additional impact for wild salmonid populations and more research is needed to assess the combined effects of multiple stressors alongside microplastic ingestion. As many studies report contrasting or neutral impacts of microplastics consumption on fish, only relatively mild impacts of microplastics ingestion were explored in this study, representative of currently available knowledge. This modelling study highlights the impacts of microplastics consumption on salmonids in freshwater stream environments and how these impacts can be mediated by personality differences. The inSTREAM model developed in the present study can be further used to explore the combined effects of multiple stressors, neurotoxic effects and changes in personality traits, on the effects of microplastics consumption on populations of stream-dwelling salmonids.
The authors are grateful to Katie Facecchia, Jenna Quinn, and Istafa Sufi for helpful discussion.
The authors have declared that no competing interests exist.
Review and/or approval by an ethics committee was not needed for this study because this is a simulation study using ABMs and no live animals were used.
This work was supported by Mitacs [IT26178]; rare Charitable Research Reserve [IT26178]; and the Natural Sciences and Engineering Research Council (NSERC) of Canada [06768].
Conceptualization: KC. Data curation: KC. Formal analysis: KC. Funding acquisition: CS, KC. Investigation: KC, SR. Methodology: CS, KC, SR. Resources: SR, CS. Software: KC, SR. Supervision: CS. Validation: KC. Visualization: KC. Writing – original draft: KC. Writing – review and editing: KC, CS, SR.
Kathleen D. W. Church https://orcid.org/0000-0002-7096-6869
Steven F. Railsback https://orcid.org/0000-0002-5923-9847
Christina A. D. Semeniuk https://orcid.org/0000-0001-5115-9853
All of the data that support the findings of this study are available in the main text or Supplementary Information.
Supplementary data
Data type: docx
Explanation note: table S1: Sensitivity analysis for stream and trout parameters altered in inSTREAM7, to evaluate the effect of each parameter on fish biomass, which accounts for abundance and body length. figure S1: Results of the sensitivity analysis for parameters altered in inSTREAM 7, and their effects on trout biomass, which accounts for abundance and body length. figure S2: The existence of real-life behavioural patterns were recreated using pattern oriented modelling. figure S3: Recruitment and mortality (± 95% C.I.s) for fish of age 1 (A, C, E) and 2 (B, D, F) nine different combinations of personality traits for the three microplastics consumption scenarios. figure S4: Recruitment and mortality over time (± 95% C.I.s) for fish of age 1 (A, C, E) and 2 (B, D, F) for the three microplastic consumption scenario. appendix S1: Methods and Results. appendix S2: Sensitivity analysis.