In temperate regions, food abundance and quality vary greatly throughout the year, and the ability of organisms to synchronise their phenology to these changes is a key determinant of their reproductive success. Successful synchronisation requires that cues are perceived prior to change, leaving time for physiological adjustments.
But what are the cues used to anticipate seasonal changes? Abiotic factors like temperature and photoperiod are known for their driving role in the phenology of a wide range of plant an animal species [1,2] . Arguably though, biotic cues directly linked to upcoming changes in food abundance could be as important as abiotic factors, but the response of organisms to these cues remains relatively unexplored.
Biotic cues may be particularly important for higher trophic levels because of their tight interaction with the hosts or preys they depend on. In this study Tougeron and colleagues [3] address this topic using interacting insects, namely herbivorous aphids and the parasitic wasps (or parasitoids) that feed on them. The key finding of the study by Tougeron et al. [3] is that the host morph in which parasitic wasp larvae develop is a major driver of diapause induction. More importantly, the aphid morph that triggers diapause in the wasp is the one that will lay overwintering eggs in autumn at the onset of harsh winter conditions. Its neatly designed experimental setup also provides evidence that this response may vary across populations as host-dependent diapause induction was only observed in a wasp population that originated from a cold area. As the authors suggests, this may be caused by local adaptation to environmental conditions because, relative to warmer regions, missing the time window to enter diapause in colder regions may have more dramatic consequences. The study also shows that different aphid morphs differ greatly in their chemical composition, and points to particular types of metabolites like sugars and polyols as specific cues for diapause induction.
This study provides a nice example of the complexity of biological interactions, and of the importance of phenological synchrony between parasites and their hosts. The authors provide evidence that phenological synchrony is likely to be achieved via chemical cues derived from the host. A similar approach was used to demonstrate that the herbivorous beetle Leptinotarsa decemlineata uses plant chemical cues to enter diapause [4]. Beetles fed on plants exposed to pre-wintering conditions entered diapause in higher proportions than those fed on control plants grown at normal conditions. As done by Tougeron et al. [3], in [4] the authors associated diapause induction to changes in the composition of metabolites in the plant. In both studies, however, the missing piece is to unveil the particular chemical involved, an answer that may be provided by future experiments.
Latitudinal clines in diapause induction have been described in a number of insect species [5]. Correlative studies, in which the phenology of different trophic levels has been monitored, suggest that these clines may in part be governed by lower trophic levels. For example, Phillimore et al. [6] explored the relative contribution of temperature and of host plant phenology on adult flight periods of the butterfly Anthocharis cardamines. Tougeron et al. [3], by using aphids and their associated parasitoids, take the field further by moving from observational studies to experiments. Besides, aphids are not only a tractable host-parasite system in the laboratory, they are important agricultural pests. Improving our basic knowledge of their ecological interactions may ultimately contribute to improving pest control techniques. The study by Tougeron et al. [3] exemplifies the multiple benefits that can be gained from addressing fundamental questions in species that are also directly relevant to society.

References

[1] Tauber, M. J., Tauber, C. A., and Masaki, S. (1986). Seasonal Adaptations of Insects. Oxford, New York: Oxford University Press.
[2] Bradshaw, W. E., and Holzapfel, C. M. (2007). Evolution of Animal Photoperiodism. Annual Review of Ecology, Evolution, and Systematics, 38(1), 1–25. doi: 10.1146/annurev.ecolsys.37.091305.110115
[3] Tougeron, K., Brodeur, J., Baaren, J. van, Renault, D., and Lann, C. L. (2019b). Sex makes them sleepy: host reproductive status induces diapause in a parasitoid population experiencing harsh winters. bioRxiv, 371385, ver. 6 peer-reviewed and recommended by PCI Ecology. doi: 10.1101/371385
[4] Izzo, V. M., Armstrong, J., Hawthorne, D., and Chen, Y. (2014). Time of the season: the effect of host photoperiodism on diapause induction in an insect herbivore, Leptinotarsa decemlineata. Ecological Entomology, 39(1), 75–82. doi: 10.1111/een.12066
[5] Hut Roelof A., Paolucci Silvia, Dor Roi, Kyriacou Charalambos P., and Daan Serge. (2013). Latitudinal clines: an evolutionary view on biological rhythms. Proceedings of the Royal Society B: Biological Sciences, 280(1765), 20130433. doi: 10.1098/rspb.2013.0433
[6] Phillimore, A. B., Stålhandske, S., Smithers, R. J., and Bernard, R. (2012). Dissecting the Contributions of Plasticity and Local Adaptation to the Phenology of a Butterfly and Its Host Plants. The American Naturalist, 180(5), 655–670. doi: 10.1086/667893

For decades, the effect of population density on individual performance has been studied by ecologists using both theoretical, observational, and experimental approaches. The generally accepted definition of the Allee effect is a positive correlation between population density and average individual fitness that occurs at low population densities, while individual fitness is typically decreased through intraspecific competition for resources at high population densities.  Allee effects are very relevant in conservation biology because species at low population densities would then be subjected to much higher extinction risks. 

However, due to all kinds of stochasticity, low population numbers are always more vulnerable to extinction than larger population sizes. This effect by itself cannot be necessarily ascribed to lower individual performance at low densities, i.e, Allee effects. Vercken and colleagues (2021) address this challenging question and measure the extent to which average individual fitness is affected by population density analyzing 30 experimental populations. As a model system, they use populations of parasitoid wasps of the genus Trichogramma. They report Allee effect in 8 out 30 experimental populations. Vercken and colleagues's work has several strengths. 

First of all, it is nice to see that they put theory at work. This is a very productive way of using theory in ecology. As a starting point, they look at what simple theoretical population models say about Allee effects (Lewis and Kareiva 1993; Amarasekare 1998; Boukal and Berec 2002). These models invariably predict a one-humped relation between population-density and per-capita growth rate. It is important to remark that pure logistic growth, the paradigm of density-dependence, would never predict such qualitative behavior. It is only when there is a depression of per-capita growth rates at low densities that true Allee effects arise. Second, these authors manage to not only experimentally test this main prediction but also report additional demographic traits that are consistently affected by population density. 

In these wasps, individual performance can be measured in terms of the average number of individuals every adult is able to put into the next generation ---the lambda parameter in their analysis. The first panel in figure 3 shows that the per-capita growth rates are lower in populations presenting Allee effects, the ones showing a one-humped behavior in the relation between per-capita growth rates and population densities (see figure 2). Also other population traits, such maximum population size and exitinction probability, change in a correlated and consistent manner. 

In sum, Vercken and colleagues's results are experimentally solid and based on theory expectations. However, they are very intriguing. They find the signature of Allee effects in only 8 out 30 populations, all from the same genus Trichogramma, and some populations belonging to the same species (from different sampling sites) do not show consistently Allee effects. Where does this population variability comes from? What are the reasons underlying this within- and between-species variability? What are the individual mechanisms driving Allee effects in these populations? Good enough, this piece of work generates more intriguing questions than the question is able to clearly answer. Science is not a collection of final answers but instead good questions are the ones that make science progress. 

References

Amarasekare P (1998) Allee Effects in Metapopulation Dynamics. The American Naturalist, 152, 298–302. https://doi.org/10.1086/286169

Boukal DS, Berec L (2002) Single-species Models of the Allee Effect: Extinction Boundaries, Sex Ratios and Mate Encounters. Journal of Theoretical Biology, 218, 375–394. https://doi.org/10.1006/jtbi.2002.3084

Lewis MA, Kareiva P (1993) Allee Dynamics and the Spread of Invading Organisms. Theoretical Population Biology, 43, 141–158. https://doi.org/10.1006/tpbi.1993.1007

Vercken E, Groussier G, Lamy L, Mailleret L (2021) The hidden side of the Allee effect: correlated demographic traits and extinction risk in experimental populations. HAL, hal-02570868, ver. 4 peer-reviewed and recommended by Peer community in Ecology. https://hal.archives-ouvertes.fr/hal-02570868

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

Tick-borne diseases are an important burden on public health all over the globe, making accurate forecasts of tick population a key ingredient in a successful public health strategy. Over long time scales, tick populations can undergo complex dynamics, as they are sensitive to many non-linear effects due to the complex relationships between ticks and the relevant (numerical) features of their environment.

But luckily, capturing complex non-linear responses is a task that machine learning thrives on. In this contribution, Manley et al. (2023) explore the use of Gradient Boosted Trees to predict the distribution (presence/absence) and abundance of ticks across New York state.

This is an interesting modelling challenge in and of itself, as it looks at the same ecological question as an instance of a classification problem (presence/absence) or of a regression problem (abundance). In using the same family of algorithm for both, Manley et al. (2023) provide an interesting showcase of the versatility of these techniques. But their article goes one step further, by setting up a multi-class categorical model that estimates jointly the presence and abundance of a population. I found this part of the article particularly elegant, as it provides an intermediate modelling strategy, in between having two disconnected models for distribution and abundance, and having nested models where abundance is only predicted for the present class (see e.g. Boulangeat et al., 2012, for a great description of the later).

One thing that Manley et al. (2023) should be commended for is their focus on opening up the black box of machine learning techniques. I have never believed that ML models are more inherently opaque than other families of models, but the focus in this article on explainable machine learning shows how these models might, in fact, bring us closer to a phenomenological understanding of the mechanisms underpinning our observations.

There is also an interesting discussion in this article, on the rate of false negatives in the different models that are being benchmarked. Although model selection often comes down to optimizing the overall quality of the confusion matrix (for distribution models, anyway), depending on the type of information we seek to extract from the model, not all types of errors are created equal. If the purpose of the model is to guide actions to control vectors of human pathogens, a false negative (predicting that the vector is absent at a site where it is actually present) is a potentially more damaging outcome, as it can lead to the vector population (and therefore, potentially, transmission) increasing unchecked.

Boulangeat I, Gravel D, Thuiller W. Accounting for dispersal and biotic interactions to disentangle the drivers of species distributions and their abundances: The role of dispersal and biotic interactions in explaining species distributions and abundances. Ecol Lett. 2012;15: 584-593.
https://doi.org/10.1111/j.1461-0248.2012.01772.x

Manley W, Tran T, Prusinski M, Brisson D. (2023) Modeling tick populations: An ecological test case for gradient boosted trees. bioRxiv, 2023.03.13.532443, ver. 3 peer-reviewed and recommended by Peer Community in Ecology. https://doi.org/10.1101/2023.03.13.532443

Scale is a big topic in ecology [1]. Environmental variation happens at particular scales. The typical scale at which organisms disperse is species-specific, but, as a first approximation, an ensemble of similar species, for instance, trees, could be considered to share a typical dispersal scale. Finally, characteristic spatial scales of species interactions are, in general, different from the typical scales of dispersal and environmental variation. Therefore, conceptually, we can distinguish these three characteristic spatial scales associated with three different processes: species selection for a given environment (E), dispersal (D), and species interactions (I), respectively.  

From the famous species-area relation to the spatial distribution of biomass and species richness, the different macro-ecological patterns we usually study emerge from an interplay between dispersal and local interactions in a physical environment that constrains species establishment and persistence in every location. To make things even more complicated, local environments are often modified by the species that thrive in them, which establishes feedback loops.  It is usually assumed that local interactions are short-range in comparison with species dispersal, and dispersal scales are typically smaller than the scales at which the environment varies (I < D < E, see [2]), but this should not always be the case. 

The authors of this paper [2] relax this typical assumption and develop a theoretical framework to study how diversity and ecosystem functioning are affected by different relations between the typical scales governing interactions, dispersal, and environmental variation. This is a huge challenge. First, diversity and ecosystem functioning across space and time have been empirically characterized through a wide variety of macro-ecological patterns. Second, accommodating local interactions, dispersal and environmental variation and species environmental preferences to model spatiotemporal dynamics of full ecological communities can be done also in a lot of different ways. One can ask if the particular approach suggested by the authors is the best choice in the sense of producing robust results, this is, results that would be predicted by alternative modeling approaches and mathematical analyses [3]. The recommendation here is to read through and judge by yourself.  

The main unusual assumption underlying the model suggested by the authors is non-local species interactions. They introduce interaction kernels to weigh the strength of the ecological interaction with distance, which gives rise to a system of coupled integro-differential equations. This kernel is the key component that allows for control and varies the scale of ecological interactions. Although this is not new in ecology [4], and certainly has a long tradition in physics ---think about the electric or the gravity field, this approach has been widely overlooked in the development of the set of theoretical frameworks we have been using over and over again in community ecology, such as the Lotka-Volterra equations or, more recently, the metacommunity concept [5].

In Physics, classic fields have been revised to account for the fact that information cannot travel faster than light. In an analogous way, a focal individual cannot feel the presence of distant neighbors instantaneously. Therefore, non-local interactions do not exist in ecological communities. As the authors of this paper point out, they emerge in an effective way as a result of non-random movements, for instance, when individuals go regularly back and forth between environments (see [6], for an application to infectious diseases), or even migrate between regions. And, on top of this type of movement, species also tend to disperse and colonize close (or far) environments. Individual mobility and dispersal are then two types of movements, characterized by different spatial-temporal scales in general. Species dispersal, on the one hand, and individual directed movements underlying species interactions, on the other, are themselves diverse across species, but it is clear that they exist and belong to two distinct categories. 

In spite of the long and rich exchange between the authors' team and the reviewers, it was not finally clear (at least, to me and to one of the reviewers) whether the model for the spatio-temporal dynamics of the ecological community (see Eq (1) in [2]) is only presented as a coupled system of integro-differential equations on a continuous landscape for pedagogical reasons, but then modeled on a discrete regular grid for computational convenience. In the latter case, the system represents a regular network of local communities,  becomes a system of coupled ODEs, and can be numerically integrated through the use of standard algorithms. By contrast,  in the former case, the system is meant to truly represent a community that develops on continuous time and space, as in reaction-diffusion systems. In that case, one should keep in mind that numerical instabilities can arise as an artifact when integrating both local and non-local spatio-temporal systems. Spatial patterns could be then transient or simply result from these instabilities. Therefore, when analyzing spatiotemporal integro-differential equations, special attention should be paid to the use of the right numerical algorithms. The authors share all their code at https://zenodo.org/record/5543191, and all this can be checked out. In any case, the whole discussion between the authors and the reviewers has inherent value in itself, because it touches on several limitations and/or strengths of the author's approach,  and I highly recommend checking it out and reading it through.

Beyond these methodological issues, extensive model explorations for the different parameter combinations are presented. Several results are reported, but, in practice, what is then the main conclusion we could highlight here among all of them?  The authors suggest that "it will be difficult to manage landscapes to preserve biodiversity and ecosystem functioning simultaneously, despite their causative relationship", because, first, "increasing dispersal and interaction scales had opposing
effects" on these two patterns, and, second, unexpectedly, "ecosystems attained the highest biomass in scenarios which also led to the lowest levels of biodiversity". If these results come to be fully robust, this is, they pass all checks by other research teams trying to reproduce them using alternative approaches, we will have to accept that we should preserve biodiversity on its own rights and not because it enhances ecosystem functioning or provides particular beneficial services to humans. 

References

[1] Levin, S. A. 1992. The problem of pattern and scale in ecology. Ecology 73:1943–1967. https://doi.org/10.2307/1941447

[2] Yuval R. Zelnik, Matthieu Barbier, David W. Shanafelt, Michel Loreau, Rachel M. Germain. 2023. Linking intrinsic scales of ecological processes to characteristic scales of biodiversity and functioning patterns. bioRxiv, ver. 2 peer-reviewed and recommended by Peer Community in Ecology.  https://doi.org/10.1101/2021.10.11.463913

[3] Baron, J. W. and Galla, T. 2020. Dispersal-induced instability in complex ecosystems. Nature Communications  11, 6032. https://doi.org/10.1038/s41467-020-19824-4

[4] Cushing, J. M. 1977. Integrodifferential equations and delay models in population dynamics 
 Springer-Verlag, Berlin. https://doi.org/10.1007/978-3-642-93073-7

[5] M. A. Leibold, M. Holyoak, N. Mouquet, P. Amarasekare, J. M. Chase, M. F. Hoopes, R. D. Holt, J. B. Shurin, R. Law, D. Tilman, M. Loreau, A. Gonzalez. 2004. The metacommunity concept: a framework for multi-scale community ecology. Ecology Letters, 7(7): 601-613. https://doi.org/10.1111/j.1461-0248.2004.00608.x

[6] M. Pardo-Araujo, D. García-García, D. Alonso, and F. Bartumeus. 2023. Epidemic thresholds and human mobility. Scientific reports 13 (1), 11409. https://doi.org/10.1038/s41598-023-38395-0