Across the temporal expanse of history, the impact of human activities on global landscapes has manifested as a complex interplay of ecological alterations. From the advent of early agricultural practices to the successive waves of industrialization characterizing the 18th and 19th centuries, anthropogenic forces have exerted profound and enduring transformations upon Earth's ecosystems. Indeed, by 2017, more than 80% of the terrestrial biosphere was transformed by human populations and land use, and just 19% remains as wildlands (Ellis et al. 2021).
 
Urbanization engenders profound alterations in environmental conditions, exerting substantial impacts on biological communities. The expansion of built infrastructure, modification of land use patterns, and the introduction of impervious surfaces and habitat fragmentation are key facets of urbanization (Faeth et al. 2011). These alterations generate biodiversity loss, changes in the composition of biological communities, disruptions in access and availability of food and nutrients, and a loss of efficiency in the immune system's control of infections, etc. (Reyes et al. 2013).
 
In this study, Caizergues et al. (2023) investigated the prevalence and diversity of avian malaria parasites (Plasmodium/Haemoproteus sp. and Leucocytozoon sp.) in great tits (Parus major) living across an urbanization gradient. The study reveals nuanced patterns of avian malaria prevalence and lineage diversity in great tits across urban and non-urban environments. While overall parasite diversity remains consistent, there are marked differences in prevalence between life stages and habitats. They observed a high prevalence in adult birds (from 95% to 100%), yet lower prevalence in fledglings (from 0% to 38%). Notably, urban nestlings exhibit higher parasite prevalence than their non-urban counterparts, suggesting a potential link between early malaria infection and the urban heat island effect. This finding underscores the importance of considering both spatial and temporal aspects of urbanization in understanding disease dynamics. Parasite lineages were not habitat-specific. The results suggest a potential parasitic burden in more urbanized areas, with a marginal but notable effect of nest-level urbanization on Plasmodium prevalence. This challenges the common perception of lower parasitic prevalence in urban environments and highlights the need for further investigation into the factors influencing parasite prevalence at finer spatial scales.
 
The discussion emphasizes the significance of examining vector distributions, abundance, and diversity in urban areas, which may be influenced by ecological niches and the presence of suitable habitats such as marshes. The identification of habitat-specific Haemosporidian lineages, particularly those occurring more frequently in urban areas, raises intriguing questions about the factors influencing parasite diversity. The presence of rare lineages in urban environments, such as AFR065, DELURB4, and YWT4, suggests a potential connection between urban bird communities and specific parasite strains.
 
Future research should empirically demonstrate these relationships to enhance our understanding of urban parasitology. This finding has broader implications for wildlife epidemiology, especially when introducing or keeping exotic wildlife in contact with native species. The study highlights the importance of considering not only the prevalence but also the specific lineages of parasites in understanding the dynamics of avian malaria in urban and non-urban habitats. This preprint contributes valuable insights to the ongoing discourse on the intricate interplay between ecological repercussions of human-induced changes (urbanization), biological communities, and the prevalence of vector-borne diseases.
 
References

Caizergues AE, Robira B, Perrier C, Jeanneau M, Berthomieu A, Perret S, Gandon S, Charmantier A (2023) Cities as parasitic amplifiers? Malaria prevalence and diversity in great tits along an urbanization gradient. bioRxiv, 2023.05.03.539263, ver. 3 peer-reviewed and recommended by Peer Community in Ecology. https://doi.org/10.1101/2023.05.03.539263

Ellis EC, Gauthier N, Klein Goldewijk K, Bliege Bird R, Boivin N, Díaz S, Fuller DQ, Gill JL, Kaplan JO, Kingston N, Locke H, McMichael CNH, Ranco D, Rick TC, Shaw MR, Stephens L, Svenning JC, Watson JEM. People have shaped most of terrestrial nature for at least 12,000 years. Proc Natl Acad Sci U S A. 2021 Apr 27;118(17):e2023483118. https://doi.org/10.1073/pnas.2023483118. 

Faeth  SH, Bang  C, Saari  S (2011) Urban biodiversity: Patterns and mechanisms. Ann N Y Acad Sci 1223:69–81. https://doi.org/10.1111/j.1749-6632.2010.05925.x

Faeth  SH, Bang  C, Saari  S (2011) Urban biodiversity: Patterns and mechanisms. Ann N Y Acad Sci 1223:69–81. https://doi.org/10.1111/j.1749-6632.2010.05925.x

Reyes  R, Ahn  R, Thurber  K, Burke  TF (2013) Urbanization and Infectious Diseases: General Principles, Historical Perspectives, and Contemporary Challenges. Challenges Infect Dis 123. https://doi.org/10.1007/978-1-4614-4496-1_4

Dispersal is a key element of the natural dynamics of meta-communities, and plays a central role in the success of populations colonizing new landscapes. Understanding how demographic processes may affect the speed at which alien species spread through environmentally-heterogeneous habitat fragments is therefore of key importance to manage biological invasions. This requires studying together the complex interplay of dispersal and population processes, two inextricably related phenomena that can produce many possible outcomes. Stochastic models offer an opportunity to describe this kind of process in a meaningful way, but to ensure that they are realistic (sensu Levins 1966) it is also necessary to combine model simulations with empirical data (Snäll et al. 2007).

Morel-Journel et al. (2023) put together stochastic models and experimental data to study how population density may affect the speed at which alien species spread through a heterogeneous landscape. They do it by focusing on what they call ‘colonisation debt’, which is merely the impact that population density at the invasion front may have on the speed at which the species colonizes patches of different carrying capacities. They investigate this issue through two largely independent approaches. First, a stochastic model of dispersal throughout the patches of a linear, 1-dimensional landscape, which accounts for different degrees of density-dependent growth. And second, a microcosm experiment of a parasitoid wasp colonizing patches with different numbers of host eggs. In both cases, they compare the velocity of colonization of patches with lower or higher carrying capacity than the previous one (i.e. what they call upward or downward gradients).

Their results show that density-dependent processes influence the speed at which new fragments are colonized is significantly reduced by positive density dependence. When either population growth or dispersal rate depend on density, colonisation debt limits the speed of invasion, which turns out to be dependent on the strength and direction of the gradient between the conditions of the invasion front, and the newly colonized patches. Although this result may be quite important to understand the meta-population dynamics of dispersing species, it is important to note that in their study the environmental differences between patches do not take into account eventual shifts in the scenopoetic conditions (i.e. the values of the environmental parameters to which species niches’ respond to; Hutchinson 1978, see also Soberón 2007). Rather, differences arise from variations in the carrying capacity of the patches that are consecutively invaded, both in the in silico and microcosm experiments. That is, they account for potential differences in the size or quality of the invaded fragments, but not on the costs of colonizing fragments with different environmental conditions, which may also determine invasion speed through niche-driven processes. This aspect can be of particular importance in biological invasions or under climate change-driven range shifts, when adaptation to new environments is often required (Sakai et al. 2001; Whitney & Gabler 2008; Hill et al. 2011).

The expansion of geographical distribution ranges is the result of complex eco-evolutionary processes where meta-community dynamics and niche shifts interact in a novel physical space and/or environment (see, e.g., Mestre et al. 2020). Here, the invasibility of native communities is determined by niche variations and how similar are the traits of alien and native species (Hui et al. 2023). Within this context, density-dependent processes will build upon and heterogeneous matrix of native communities and environments (Tischendorf et al. 2005), to eventually determine invasion success. What the results of Morel-Journel et al. (2023) show is that, when the invader shows density dependence, the invasion process can be slowed down by variations in the carrying capacity of patches along the dispersal front. This can be particularly useful to manage biological invasions; ongoing invasions can be at least partially controlled by manipulating the size or quality of the patches that are most adequate to the invader, controlling host populations to reduce carrying capacity. But further, landscape manipulation of such kind could be used in a preventive way, to account in advance for the effects of the introduction of alien species for agricultural exploitation or biological control, thereby providing an additional safeguard to practices such as the introduction of parasitoids to control plagues. These practical aspects are certainly worth exploring further, together with a more explicit account of the influence of the abiotic conditions and the characteristics of the invaded communities on the success and speed of biological invasions.

REFERENCES

Hill, J.K., Griffiths, H.M. & Thomas, C.D. (2011) Climate change and evolutionary adaptations at species' range margins. Annual Review of Entomology, 56, 143-159. https://doi.org/10.1146/annurev-ento-120709-144746

Hui, C., Pyšek, P. & Richardson, D.M. (2023) Disentangling the relationships among abundance, invasiveness and invasibility in trait space. npj Biodiversity, 2, 13. https://doi.org/10.1038/s44185-023-00019-1

Hutchinson, G.E. (1978) An introduction to population biology. Yale University Press, New Haven, CT.

Levins, R. (1966) The strategy of model building in population biology. American Scientist, 54, 421-431. 

Mestre, A., Poulin, R. & Hortal, J. (2020) A niche perspective on the range expansion of symbionts. Biological Reviews, 95, 491-516. https://doi.org/10.1111/brv.12574

Morel-Journel, T., Haond, M., Duan, L., Mailleret, L. & Vercken, E. (2023) Colonisation debt: when invasion history impacts current range expansion. bioRxiv, 2022.11.13.516255, ver. 3 peer-reviewed and recommended by Peer Community in Ecology. https://doi.org/10.1101/2022.11.13.516255

Snäll, T., B. O'Hara, R. & Arjas, E. (2007) A mathematical and statistical framework for modelling dispersal. Oikos, 116, 1037-1050. https://doi.org/10.1111/j.0030-1299.2007.15604.x

Sakai, A.K., Allendorf, F.W., Holt, J.S., Lodge, D.M., Molofsky, J., With, K.A., Baughman, S., Cabin, R.J., Cohen, J.E., Ellstrand, N.C., McCauley, D.E., O'Neil, P., Parker, I.M., Thompson, J.N. & Weller, S.G. (2001) The population biology of invasive species. Annual Review of Ecology and Systematics, 32, 305-332. https://doi.org/10.1146/annurev.ecolsys.32.081501.114037

Soberón, J. (2007) Grinnellian and Eltonian niches and geographic distributions of species. Ecology Letters, 10, 1115-1123. https://doi.org/10.1111/j.1461-0248.2007.01107.x

Tischendorf, L., Grez, A., Zaviezo, T. & Fahrig, L. (2005) Mechanisms affecting population density in fragmented habitat. Ecology and Society, 10, 7. https://doi.org/10.5751/ES-01265-100107

Whitney, K.D. & Gabler, C.A. (2008) Rapid evolution in introduced species, 'invasive traits' and recipient communities: challenges for predicting invasive potential. Diversity and Distributions, 14, 569-580. https://doi.org/10.1111/j.1472-4642.2008.00473.x

Different sources of uncertainty are known to affect our ability to predict ecological dynamics (Petchey et al. 2015). However, the consequences of uncertainty on prediction biases have been less investigated, especially when predictions are scaled up to higher levels of organisation as is commonly done in ecology for instance. The study of Orr et al. (2020) addresses this issue. It shows that, in complex systems, the uncertainty of unbiased predictions at a lower level of organisation (e.g. species level) leads to a bias towards underestimation of change at higher level of organisation (e.g. ecosystem level). This bias is strengthened by larger uncertainty and by higher dimensionality of the system.
This general result has large implications for many fields of science, from economics to energy supply or demography. In ecology, as discussed in this study, these results imply that the uncertainty of predictions of species’ change increases the probability of underestimation of changes of diversity and stability at community and ecosystem levels, especially when species richness is high. The uncertainty of predictions of species’ change also increases the probability of underestimation of change when multiple ecosystem functions are considered at once, or when the combined effect of multiple stressors is considered.
The consequences of species diversity on ecosystem functions and stability have received considerable attention during the last decades (e.g. Cardinale et al. 2012, Kéfi et al. 2019). However, since the bias towards underestimation of change increases with species diversity, future studies will need to investigate how the general statistical effect outlined by Orr et al. might affect our understanding of the well-known relationships between species diversity and ecosystem functioning and stability in response to perturbations.

References

Cardinale BJ, Duffy JE, Gonzalez A, Hooper DU, Perrings C, Venail P, Narwani A, Mace GM, Tilman D, Wardle DA, Kinzig AP, Daily GC, Loreau M, Grace JB, Larigauderie A, Srivastava DS, Naeem S (2012) Biodiversity loss and its impact on humanity. Nature, 486, 59–67. https://doi.org/10.1038/nature11148
Kéfi S, Domínguez‐García V, Donohue I, Fontaine C, Thébault E, Dakos V (2019) Advancing our understanding of ecological stability. Ecology Letters, 22, 1349–1356. https://doi.org/10.1111/ele.13340
Orr JA, Piggott JJ, Jackson A, Arnoldi J-F (2020) Why scaling up uncertain predictions to higher levels of organisation will underestimate change. bioRxiv, 2020.05.26.117200. https://doi.org/10.1101/2020.05.26.117200
Petchey OL, Pontarp M, Massie TM, Kéfi S, Ozgul A, Weilenmann M, Palamara GM, Altermatt F, Matthews B, Levine JM, Childs DZ, McGill BJ, Schaepman ME, Schmid B, Spaak P, Beckerman AP, Pennekamp F, Pearse IS (2015) The ecological forecast horizon, and examples of its uses and determinants. Ecology Letters, 18, 597–611. https://doi.org/10.1111/ele.12443

Biodiversity loss can have important consequences for ecosystem functions, as exemplified by a large body of literature spanning at least three decades [1–3]. While connections between species diversity and ecosystem functions are now well-defined and understood, the importance of diversity within species is more elusive. Despite a surge in theoretical work on how intraspecific diversity can affect coexistence in simple community types [4,5], not much is known about how intraspecific diversity drives ecosystem processes in more complex community types. One particular challenge is that intraspecific diversity can be expressed as observable variation of functional traits, or instead subsist as genetic variation of which the consequences for ecosystem processes are harder to grasp.
Raffard et al. [6] examined how intraspecific biodiversity loss in a consumer fish changes species diversity at lower trophic levels and ecosystem processes in pond mesocosms. An interesting feature of this experiment is that it crosses functional and genetic intraspecific diversity. To do so, Raffard and colleagues measured and genotyped European minnow (P. phoxinus) individuals sampled from streams across southern France. Combining these collected specimens into experimental ponds allowed them to control functional (population variance of body size) and genetic intraspecific richness (number of genotypes).
Effects on minnow biomass production were mostly small; biomass was significantly reduced only when lowering both functional and genetic richness. However, the consequences for lower trophic levels (zooplankton and macroinvertebrates) were more pronounced and – importantly – not intuitive. For instance, the macroinvertebrate community was less species-diverse at higher minnow functional richness. If minnows with different body sizes would be the main regulator factors [7] explaining macroinvertebrate interactions, one would expect a more diverse set of minnow body sizes (i.e. higher functional minnow richness) to permit higher instead of lower macroinvertebrate richness. At the same time, the macroinvertebrate community was more species-diverse at higher minnow genotype richness, which could indicate unobserved minnow traits determining macroinvertebrate diversity more than the usual suspects (functional consumer richness). Such unobserved traits could be behavioral traits, allowing for resource partitioning among fish.
The consequences of functional minnow diversity loss on zooplankton diversity were negative, as expected in case body size differences among fish would facilitate coexistence of their zooplankton prey, as explained above. However, this was only the case when genetic diversity was high, suggesting nonstraightforward interactive effects of observed and non-observed traits on prey diversity.
The effects of functional and genetic minnow diversity loss on invertebrate (macroinvertebrates and zooplankton) abundance were more consistent than for invertebrate diversity. This suggests again nonstraightforward relationships in this experimental ecosystem, but now between invertebrate diversity and abundance. When using abundance as a proxy for an ecosystem process (which the authors did not), this result illustrates that biodiversity loss in multitrophic communities can have consequences that are challenging to interpret, let alone predict [8,9]. Path analyses showed how the observed changes of invertebrate diversity and abundance co-determined decomposition, a key ecosystem function. These path analyses had highest explanatory power show when including both kinds of intraspecific diversity.
Taken together, the results by Raffard and colleagues suggest that genetic consumer richness can drive species diversity of connected trophic levels and ecosystem processes with similar magnitude as functional diversity. Indeed, the effects of genetic consumer richness were shown to be so strong as to compensate or exacerbate the loss of observed functional richness. The exact mechanisms explaining these effects remain to be identified, however. The possibility that fish grazing by fish with different (observed or not observed) traits regulates coexistence among invertebrate prey, for instance, would depend on how strong fish consumption feeds back on prey growth during a 30-week experiment. As the authors indicate, detailed studies on resource partitioning among consumers (e.g. using stable isotope labelling) can shed light on these matters. Doing so may address a more fundamental question, which is if the mechanisms linking intraspecific diversity to function are different from those linking interspecific diversity to function, and at what time scales.

References

[1] Tilman D, Downing JA (1994) Biodiversity and stability in grasslands. Nature, 367, 363–365. https://doi.org/10.1038/367363a0
[2] Cardinale BJ, Duffy JE, Gonzalez A, Hooper DU, Perrings C, Venail P, Narwani A, Mace GM, Tilman D, Wardle DA, Kinzig AP, Daily GC, Loreau M, Grace JB, Larigauderie A, Srivastava DS, Naeem S (2012) Biodiversity loss and its impact on humanity. Nature, 486, 59–67. https://doi.org/10.1038/nature11148
[3] De Laender F, Rohr JR, Ashauer R, Baird DJ, Berger U, Eisenhauer N, Grimm V, Hommen U, Maltby L, Meliàn CJ, Pomati F, Roessink I, Radchuk V, Brink PJV den (2016) Reintroducing Environmental Change Drivers in Biodiversity–Ecosystem Functioning Research. Trends in Ecology & Evolution, 31, 905–915. https://doi.org/10.1016/j.tree.2016.09.007
[4] Hart SP, Schreiber SJ, Levine JM (2016) How variation between individuals affects species coexistence. Ecology Letters, 19, 825–838. https://doi.org/10.1111/ele.12618
[5] Barabás G, D’Andrea R (2016) The effect of intraspecific variation and heritability on community pattern and robustness. Ecology Letters, 19, 977–986. https://doi.org/10.1111/ele.12636
[6] Raffard A, Cucherousset J, Montoya JM, Richard M, Acoca-Pidolle S, Poésy C, Garreau A, Santoul F, Blanchet S (2020) Intraspecific diversity loss in a predator species alters prey community structure and ecosystem functions. bioRxiv, 2020.06.10.144337, ver. 3 peer-reviewed and recommended by PCI Ecology. https://doi.org/10.1101/2020.06.10.144337
[7] Pásztor L, Botta-Dukát Z, Magyar G, Czárán T, Meszéna G. Theory-Based Ecology: A Darwinian approach. Oxford University Press. https://doi.org/10.1093/acprof:oso/9780199577859.001.0001
[8] Binzer A, Guill C, Rall BC, Brose U (2016) Interactive effects of warming, eutrophication and size structure: impacts on biodiversity and food-web structure. Global Change Biology, 22, 220–227. https://doi.org/10.1111/gcb.13086
[9] Schwarz B, Barnes AD, Thakur MP, Brose U, Ciobanu M, Reich PB, Rich RL, Rosenbaum B, Stefanski A, Eisenhauer N (2017) Warming alters energetic structure and function but not resilience of soil food webs. Nature Climate Change, 7, 895–900. https://doi.org/10.1038/s41558-017-0002-z

Marginal value theorem (MVT) is an archetypal model discussed in every behavioural ecology textbook. Its popularity is largely explained but the fact that it is possible to solve it graphically (at least in its simplest form) with the minimal amount of equations, which is a sensible strategy for an introductory course in behavioural ecology [1]. Apart from this heuristic value, one may be tempted to disregard it as a naive toy model. After a burst of interest in the 70's and the 80's, the once vivid literature about optimal foraging theory (OFT) has lost its momentum [2]. Yet, OFT and MVT have remained an active field of research in the parasitoidologists community, mostly because the sampling strategy of a parasitoid in patches of hosts and its resulting fitness gain are straightforward to evaluate, which eases both experimental and theoretical investigations [3].
This preprint [4] is in line with the long-established literature on OFT. It follows two theoretical articles [5,6] in which Vincent Calcagno and co-authors assessed the effect of changes in the environmental conditions on optimal foraging strategy. This time, they did not modify the shape of the gain function (describing the diminishing return of the cumulative intake as a function of the residency time in a patch) but the relative frequencies of good and bad patches. At first sight, that sounds like a minor modification of their earlier models. Actually, even the authors initially were fooled by the similarities before spotting the pitfalls. Here, they genuinely point out the erroneous verbal prediction in their previous paper in which some non-trivial effects of the change in patch frequencies have been overlooked. The present study indeed provides a striking example of ecological fallacy, and more specifically of Simpson's paradox which occurs when the aggregation of subgroups modifies the apparent pattern at the scale of the entire population [7,8]. In the case of MVT under constraints of habitat conversion, the increase of the residency times in both bad and good patches can result in a decrease of the average residency time at the level of the population. This apparently counter-intuitive property can be observed, for instance, when the proportion of bad quality patches strongly increases, which increases the probability that the individual forages on theses quickly exploited patches, and thus decreases its average residency time on the long run.
The authors thus put the model on the drawing board again. Proper assessment of the effect of change in the frequency of patch quality is more mathematically challenging than when one considers only changes in the shape of the gain function. The expected gain must be evaluated at the scale of the entire habitat instead of single patch. Overall, this study, which is based on a rigorous formalism, stands out as a warning against too rapid interpretations of theoretical outputs. It is not straightforward to generalize the predictions of previous models without careful evaluating their underlying hypotheses. The devil is in the details: some slight, seemingly minor, adjustments of the assumptions may have some major consequences.
The authors discussed the general conditions leading to changes in residency times or movement rates. Yet, it is worth pointing out again that it would be a mistake to blindly consider these theoretical results as forecasts for the foragers' behaviour in natura. OFT models has for a long time been criticized for sweeping under the carpet the key questions of the evolutionary dynamics and the maintenance of the optimal strategy in a population [9,10]. The distribution of available options is susceptible to change rapidly due to modifications of the environmental conditions or, even more simply, the presence of competitors which continuously remove the best options from the pool of available options [11]. The key point here is that the constant monitoring of available options implies cognitive (neural tissue is one of the most metabolically expensive tissues) and ecological costs: assessment and adjustment to the environmental conditions requires time, energy, and occasional mistakes (cost of naiveté, [12]). While rarely considered in optimal analyses, these costs should severely constraint the evolution of the subtle decision rules. Under rapidly fluctuating conditions, it could be more profitable to maintain a sub-optimal strategy (but performing reasonably well on the long run) than paying the far from negligible costs implied by the pursuit of optimal strategies [13,14]. For instance, in the analysis presented in this preprint, it is striking how close the fitness gains of the plastic and the non-plastic forager are, particularly if one remembers that the last-mentioned cognitive and ecological costs have been neglected in these calculations.
Yet, even if one can arguably question its descriptive value, such models are worth more than a cursory glance. They still have normative value insofar that they provide upper bounds for the response to modifications of the environmental conditions. Such insights are precious to design future experiments on the question. Being able to compare experimentally measured behaviours with the extremes of the null model (stubborn non-plastic forager) and the optimal strategy (only achievable by an omniscient daemon) informs about the cognitive bias or ecological costs experienced by real life foragers. I thus consider that this model, and more generally most OFT models, are still a valuable framework which deserves further examination.

References

[1] Fawcett, T. W. & Higginson, A. D. 2012 Heavy use of equations impedes communication among biologists. Proc. Natl. Acad. Sci. 109, 11735–11739. doi: 10.1073/pnas.1205259109
[2] Owens, I. P. F. 2006 Where is behavioural ecology going? Trends Ecol. Evol. 21, 356–361. doi: 10.1016/j.tree.2006.03.014
[3] Louâpre, P., Fauvergue, X., van Baaren, J. & Martel, V. 2015 The male mate search: an optimal foraging issue? Curr. Opin. Insect Sci. 9, 91–95. doi: 10.1016/j.cois.2015.02.012
[4] Calcagno, V., Hamelin, F., Mailleret, L., & Grognard, F. (2018). How optimal foragers should respond to habitat changes? On the consequences of habitat conversion. bioRxiv, 273557, ver. 4 peer-reviewed and recommended by PCI Ecol. doi: 10.1101/273557
[5] Calcagno, V., Grognard, F., Hamelin, F. M., Wajnberg, É. & Mailleret, L. 2014 The functional response predicts the effect of resource distribution on the optimal movement rate of consumers. Ecol. Lett. 17, 1570–1579. doi: 10.1111/ele.12379
[6] Calcagno, V., Mailleret, L., Wajnberg, É. & Grognard, F. 2013 How optimal foragers should respond to habitat changes: a reanalysis of the Marginal Value Theorem. J. Math. Biol. 69, 1237–1265. doi: 10.1007/s00285-013-0734-y
[7] Galipaud, M., Bollache, L., Wattier, R., Dechaume-Moncharmont, F.-X. & Lagrue, C. 2015 Overestimation of the strength of size-assortative pairing in taxa with cryptic diversity: a case of Simpson's paradox. Anim. Behav. 102, 217–221. doi: 10.1016/j.anbehav.2015.01.032
[8] Kievit, R. A., Frankenhuis, W. E., Waldorp, L. J. & Borsboom, D. 2013 Simpson's paradox in psychological science: a practical guide. Front. Psychol. 4, 513. doi: 10.3389/fpsyg.2013.00513
[9] Bolduc, J.-S. & Cézilly, F. 2012 Optimality modelling in the real world. Biol. Philos. 27, 851–869. doi: 10.1007/s10539-012-9333-3
[10] Pierce, G. J. & Ollason, J. G. 1987 Eight reasons why optimal foraging theory is a complete waste of time. Oikos 49, 111–118. doi: 10.2307/3565560
[11] Dechaume-Moncharmont, F.-X., Brom, T. & Cézilly, F. 2016 Opportunity costs resulting from scramble competition within the choosy sex severely impair mate choosiness. Anim. Behav. 114, 249–260. doi: 10.1016/j.anbehav.2016.02.019
[12] Snell-Rood, E. C. 2013 An overview of the evolutionary causes and consequences of behavioural plasticity. Anim. Behav. 85, 1004–1011. doi: 10.1016/j.anbehav.2012.12.031
[13] Fawcett, T. W., Fallenstein, B., Higginson, A. D., Houston, A. I., Mallpress, D. E. W., Trimmer, P. C. & McNamara, J. M. 2014 The evolution of decision rules in complex environments. Trends Cogn. Sci. 18, 153–161. doi: 10.1016/j.tics.2013.12.012
[14] Marshall, J. A. R., Trimmer, P. C., Houston, A. I. & McNamara, J. M. 2013 On evolutionary explanations of cognitive biases. Trends Ecol. Evol. 28, 469-473. doi: 10.1016/j.tree.2013.05.013