Submit a preprint

Direct submissions to PCI Ecology from bioRxiv.org are possible using the B2J service

Latest recommendations

IdTitle * Authors * Abstract * Picture * Thematic fields * RecommenderReviewersSubmission date
01 Mar 2022
article picture

Dissimilarity of species interaction networks: quantifying the effect of turnover and rewiring

How to evaluate and interpret the contribution of species turnover and interaction rewiring when comparing ecological networks?

Recommended by ORCID_LOGO based on reviews by Ignasi Bartomeus and 1 anonymous reviewer

A network includes a set of vertices or nodes (e.g., species in an interaction network), and a set of edges or links (e.g., interactions between species). Whether and how networks vary in space and/or time are questions often addressed in ecological research. 

Two ecological networks can differ in several extents: in that species are different in the two networks and establish new interactions (species turnover), or in that species that are present in both networks establish different interactions in the two networks (rewiring). The ecological meaning of changes in network structure is quite different according to whether species turnover or interaction rewiring plays a greater role. Therefore, much attention has been devoted in recent years on quantifying and interpreting the relative changes in network structure due to species turnover and/or rewiring.

Poisot et al. (2012) proposed to partition the global variation in structure between networks, \( \beta_{WN} \) (WN = Whole Network) into two terms: \( \beta_{OS} \) (OS = Only Shared species) and \( \beta_{ST} \) (ST = Species Turnover), such as \( \beta_{WN} = \beta_{OS} + \beta_{ST} \).

The calculation lays on enumerating the interactions between species that are common or not to two networks, as illustrated on Figure 1 for a simple case. Specifically, Poisot et al. (2012) proposed to use a Sorensen type measure of network dissimilarity, i.e., \( \beta_{WN} = \frac{a+b+c}{(2a+b+c)/2} -1=\frac{b+c}{2a+b+c} \) , where \( a \) is the number of interactions shared between the networks, while \( b \) and \( c \) are interaction numbers unique to one and the other network, respectively. \( \beta_{OS} \) is calculated based on the same formula, but only for the subnetworks including the species common to the two networks, in the form \( \beta_{OS} = \frac{b_{OS}+c_{OS}}{2a_{OS}+b_{OS}+c_{OS}} \) (e.g., Fig. 1). \( \beta_{ST} \) is deduced by subtracting \( \beta_{OS} \) from \( \beta_{WN} \) and represents in essence a "dissimilarity in interaction structure introduced by dissimilarity in species composition" (Poisot et al. 2012).

Figure 1. Ecological networks exemplified in Fründ (2021) and discussed in Poisot (2022). a is the number of shared links (continuous lines in right figures), while b+c is the number of edges unique to one or the other network (dashed lines in right figures).

Alternatively, Fründ (2021) proposed to define \( \beta_{OS} = \frac{b_{OS}+c_{OS}}{2a+b+c} \) and \( \beta_{ST} = \frac{b_{ST}+c_{ST}}{2a+b+c} \), where \( b_{ST}=b-b_{OS} \)  and \( c_{ST}=c-c_{OS} \) , so that the components \( \beta_{OS} \) and \( \beta_{ST} \) have the same denominator. In this way, Fründ (2021) partitioned the count of unique \( b+c=b_{OS}+b_{ST}+c_{ST} \) interactions, so that \( \beta_{OS} \) and \( \beta_{ST} \) sums to \( \frac{b_{OS}+c_{OS}+b_{ST}+c_{ST}}{2a+b+c} = \frac{b+c}{2a+b+c} = \beta_{WN} \). Fründ (2021) advocated that this partition allows a more sensible comparison of \( \beta_{OS} \) and \( \beta_{ST} \), in terms of the number of links that contribute to each component.

For instance, let us consider the networks 1 and 2 in Figure 1 (left panel) such as \( a_{OS}=2 \) (continuous lines in right panel), \( b_{ST} + c_{ST} = 1 \) and \( b_{OS} + c_{OS} = 1 \) (dashed lines in right panel), and thereby \( a = 2 \), \( b+c=2 \), \( \beta_{WN} = 1/3 \). Fründ (2021) measured \( \beta_{OS}=\beta_{ST}=1/6 \) and argued that it is appropriate insofar as it reflects that the number of unique links in the OS and ST components contributing to network dissimilarity (dashed lines) are actually equal. Conversely, the formula of Poisot et al. (2012) yields \( \beta_{OS}=1/5 \), hence \( \beta_{ST} = \frac{1}{3}-\frac{1}{5}=\frac{2}{15}<\beta_{OS} \). Fründ (2021) thus argued that the method of Poisot tends to underestimate the contribution of species turnover.

To clarify and avoid misinterpretation of the calculation of \( \beta_{OS} \) and \( \beta_{ST} \) in Poisot et al. (2012), Poisot (2022) provides a new, in-depth mathematical analysis of the decomposition of \( \beta_{WN} \). Poisot et al. (2012) quantify in \( \beta_{OS} \) the actual contribution of rewiring in network structure for the subweb of common species. Poisot (2022) thus argues that \( \beta_{OS} \) relates only to the probability of rewiring in the subweb, while the definition of \( \beta_{OS} \) by Fründ (2021) is relative to the count of interactions in the global network (considered in denominator), and is thereby dependent on both rewiring probability and species turnover. Poisot (2022) further clarifies the interpretation of \( \beta_{ST} \). \( \beta_{ST} \) is obtained by subtracting \( \beta_{OS} \) from \( \beta_{WN} \) and thus represents the influence of species turnover in terms of the relative architectures of the global networks and of the subwebs of shared species. Coming back to the example of Fig.1., the Poisot et al. (2012) formula posits that \( \frac{\beta_{ST}}{\beta_{WN}}=\frac{2/15}{1/3}=2/5 \), meaning that species turnover contributes two-fifths of change in network structure, while rewiring in the subweb of common species contributed three fifths.  Conversely, the approach of Fründ (2021) does not compare the architectures of global networks and of the subwebs of shared species, but considers the relative contribution of unique links to network dissimilarity in terms of species turnover and rewiring. 

Poisot (2022) concludes that the partition proposed in Fründ (2021) does not allow unambiguous ecological interpretation of rewiring. He provides guidelines for proper interpretation of the decomposition proposed in Poisot et al. (2012).

References

Fründ J (2021) Dissimilarity of species interaction networks: how to partition rewiring and species turnover components. Ecosphere, 12, e03653. https://doi.org/10.1002/ecs2.3653

Poisot T, Canard E, Mouillot D, Mouquet N, Gravel D (2012) The dissimilarity of species interaction networks. Ecology Letters, 15, 1353–1361. https://doi.org/10.1111/ele.12002

Poisot T (2022) Dissimilarity of species interaction networks: quantifying the effect of turnover and rewiring. EcoEvoRxiv Preprints, ver. 4 peer-reviewed and recommended by Peer Community in Ecology. https://doi.org/10.32942/osf.io/gxhu2

Dissimilarity of species interaction networks: quantifying the effect of turnover and rewiringTimothée Poisot<p style="text-align: justify;">Despite having established its usefulness in the last ten years, the decomposition of ecological networks in components allowing to measure their β-diversity retains some methodological ambiguities. Notably, how to ...Biodiversity, Interaction networks, Theoretical ecologyFrançois Munoz2021-07-31 00:18:41 View
01 Jun 2018
article picture

Data-based, synthesis-driven: setting the agenda for computational ecology

Some thoughts on computational ecology from people who I’m sure use different passwords for each of their accounts

Recommended by based on reviews by Matthieu Barbier and 1 anonymous reviewer

Are you an ecologist who uses a computer or know someone that does? Even if your research doesn’t rely heavily on advanced computational techniques, it likely hasn’t escaped your attention that computers are increasingly being used to analyse field data and make predictions about the consequences of environmental change. So before artificial intelligence and robots take over from scientists, now is great time to read about how experts think computers could make your life easier and lead to innovations in ecological research. In “Data-based, synthesis-driven: setting the agenda for computational ecology”, Poisot and colleagues [1] provide a brief history of computational ecology and offer their thoughts on how computational thinking can help to bridge different types of ecological knowledge. In this wide-ranging article, the authors share practical strategies for realising three main goals: (i) tighter integration of data and models to make predictions that motivate action by practitioners and policy-makers; (ii) closer interaction between data-collectors and data-users; and (iii) enthusiasm and aptitude for computational techniques in future generations of ecologists. The key, Poisot and colleagues argue, is for ecologists to “engage in meaningful dialogue across disciplines, and recognize the currencies of their collaborations.” Yes, this is easier said than done. However, the journey is much easier with a guide and when everyone involved serves to benefit not only from the eventual outcome, but also the process.

References

[1] Poisot, T., Labrie, R., Larson, E., & Rahlin, A. (2018). Data-based, synthesis-driven: setting the agenda for computational ecology. BioRxiv, 150128, ver. 4 recommended and peer-reviewed by PCI Ecology. doi: 10.1101/150128

Data-based, synthesis-driven: setting the agenda for computational ecologyTimothée Poisot, Richard Labrie, Erin Larson, Anastasia Rahlin<p>Computational ecology, defined as the application of computational thinking to ecological problems, has the potential to transform the way ecologists think about the integration of data and models. As the practice is gaining prominence as a way...Meta-analyses, Statistical ecology, Theoretical ecologyPhillip P.A. Staniczenko2018-02-05 20:51:41 View
22 May 2019
article picture

Sex makes them sleepy: host reproductive status induces diapause in a parasitoid population experiencing harsh winters

The response of interacting species to biotic seasonal cues

Recommended by and based on reviews by Anne Duplouy and 1 anonymous reviewer

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

Sex makes them sleepy: host reproductive status induces diapause in a parasitoid population experiencing harsh wintersTougeron K., Brodeur J., van Baaren J., Renault D. and Le Lann C.<p>When organisms coevolve, any change in one species can induce phenotypic changes in traits and ecology of the other species. The role such interactions play in ecosystems is central, but their mechanistic bases remain underexplored. Upper troph...Coexistence, Evolutionary ecology, Experimental ecology, Host-parasite interactions, PhysiologyAdele Mennerat2018-07-18 18:51:03 View
29 May 2023
article picture

Using integrated multispecies occupancy models to map co-occurrence between bottlenose dolphins and fisheries in the Gulf of Lion, French Mediterranean Sea

Mapping co-occurence of human activities and wildlife from multiple data sources

Recommended by based on reviews by Mason Fidino and 1 anonymous reviewer

Two fields of research have grown considerably over the past twenty years: the investigation of human-wildlife conflicts (e.g. see Treves & Santiago-Ávila 2020), and multispecies occupancy modelling (Devarajan et al. 2020). In their recent study, Lauret et al. (2023) combined both in an elegant methodological framework, applied to the study of the co-occurrence of fishing activities and bottlenose dolphins in the French Mediterranean.

A common issue with human-wildlife conflicts (and, in particular, fishery by-catch) is that data is often only available from those conflicts or interactions, limiting the validity of the predictions (Kuiper et al. 2022). Lauret et al. use independent data sources informing the occurrence of fishing vessels and dolphins, combined in a Bayesian multispecies occupancy model where vessels are "the other species". I particularly enjoyed that approach, as integration of human activities in ecological models can be extremely complex, but can also translate in phenomena that can be captured as one would of individuals of a species, as long as the assumptions are made clearly. Here, the model is made more interesting by accounting for environmental factors (seabed depth) borrowing an approach from Generalized Additive Models in the Bayesian framework. While not pretending to provide (yet) practical recommendations to help conserve bottlenose dolphins (and other wildlife conflicts), this study and the associated code are a promising step in that direction.

REFERENCES

Devarajan, K., Morelli, T.L. & Tenan, S. (2020), Multi-species occupancy models: review, roadmap, and recommendations. Ecography, 43: 1612-1624. https://doi.org/10.1111/ecog.04957

Kuiper, T., Loveridge, A.J. and Macdonald, D.W. (2022), Robust mapping of human–wildlife conflict: controlling for livestock distribution in carnivore depredation models. Anim. Conserv., 25: 195-207. https://doi.org/10.1111/acv.12730

Lauret V, Labach H, David L, Authier M, & Gimenez O (2023) Using integrated multispecies occupancy models to map co-occurrence between bottlenose dolphins and fisheries in the Gulf of Lion, French Mediterranean Sea. Ecoevoarxiv, ver. 2 peer-reviewed and recommended by PCI Ecology. https://doi.org/10.32942/osf.io/npd6u

Treves, A. & Santiago-Ávila, F.J. (2020). Myths and assumptions about human-wildlife conflict and coexistence. Conserv. Biol. 34, 811–818.  https://doi.org/10.1111/cobi.13472

Using integrated multispecies occupancy models to map co-occurrence between bottlenose dolphins and fisheries in the Gulf of Lion, French Mediterranean SeaValentin Lauret, Hélène Labach, Léa David, Matthieu Authier, Olivier Gimenez<p style="text-align: justify;">In the Mediterranean Sea, interactions between marine species and human activities are prevalent. The coastal distribution of bottlenose dolphins (<em>Tursiops truncatus</em>) and the predation pressure they put on ...Marine ecology, Population ecology, Species distributionsPaul Caplat2022-10-21 11:13:36 View
22 Apr 2021
article picture

The hidden side of the Allee effect: correlated demographic traits and extinction risk in experimental populations

Allee effects under the magnifying glass

Recommended by ORCID_LOGO based on reviews by Tom Van Dooren, Dani Oro and 1 anonymous reviewer

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

The hidden side of the Allee effect: correlated demographic traits and extinction risk in experimental populationsVercken Elodie, Groussier Géraldine, Lamy Laurent, Mailleret Ludovic<p style="text-align: justify;">Because Allee effects (i.e., the presence of positive density-dependence at low population size or density) have major impacts on the dynamics of small populations, they are routinely included in demographic models ...Demography, Experimental ecology, Population ecologyDavid Alonso2020-09-30 16:38:29 View
05 Apr 2019
article picture

Using a large-scale biodiversity monitoring dataset to test the effectiveness of protected areas at conserving North-American breeding birds

Protected Areas effects on biodiversity: a test using bird data that hopefully will give ideas for much more studies to come

Recommended by based on reviews by Willson Gaul and 1 anonymous reviewer

In the face of worldwide declines in biodiversity, evaluating the effectiveness of conservation practices is an absolute necessity. Protected Areas (PA) are a key tool for conservation, and the question “Are PA effective” has been on many a research agenda, as the introduction to this preprint will no doubt convince you. A challenge we face is that, until now, few studies have been explicitly designed to evaluate PA, and despite the rise of meta-analyses on the topic, our capacity to quantify their effect on biodiversity remains limited.
This study by Cazalis et al. [1] uses the rich dataset of the North-American Breeding Bird Survey and a sound paired design to investigate how PA change bird assemblages. The methodological care brought to the study in itself is worth the read, and the results are insightful. I will not spoil too much by revealing here that things are “complicated”, and that effects – or lack thereof – depend on the type of ecosystem, and the type of species considered.
If you are interested in conservation, bird communities, species life-history, or like beautiful plots: go and read it.

References

[1] Cazalis, V., Belghali, S., & Rodrigues, A. S. (2019). Using a large-scale biodiversity monitoring dataset to test the effectiveness of protected areas at conserving North-American breeding birds. bioRxiv, 433037, ver. 4 peer-reviewed and recommended by PCI Ecology. doi: 10.1101/433037

Using a large-scale biodiversity monitoring dataset to test the effectiveness of protected areas at conserving North-American breeding birdsVictor Cazalis, Soumaya Belghali, Ana S.L. Rodrigues<p>Protected areas currently cover about 15% of the global land area, and constitute one of the main tools in biodiversity conservation. Quantifying their effectiveness at protecting species from local decline or extinction involves comparing prot...Biodiversity, Conservation biology, Human impact, Landscape ecology, MacroecologyPaul Caplat2018-10-04 08:43:34 View
02 Oct 2018
article picture

How optimal foragers should respond to habitat changes? On the consequences of habitat conversion.

Optimal foraging in a changing world: old questions, new perspectives

Recommended by ORCID_LOGO based on reviews by Frederick Adler, Andrew Higginson and 1 anonymous reviewer

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

How optimal foragers should respond to habitat changes? On the consequences of habitat conversion.Vincent Calcagno, Frederic Hamelin, Ludovic Mailleret, Frederic GrognardThe Marginal Value Theorem (MVT) provides a framework to predict how habitat modifications related to the distribution of resources over patches should impact the realized fitness of individuals and their optimal rate of movement (or patch residen...Behaviour & Ethology, Dispersal & Migration, Foraging, Landscape ecology, Spatial ecology, Metacommunities & Metapopulations, Theoretical ecologyFrancois-Xavier Dechaume-Moncharmont2018-03-05 10:42:11 View
14 Dec 2018
article picture

Recommendations to address uncertainties in environmental risk assessment using toxicokinetics-toxicodynamics models

Addressing uncertainty in Environmental Risk Assessment using mechanistic toxicological models coupled with Bayesian inference

Recommended by based on reviews by Andreas Focks and 2 anonymous reviewers

Environmental Risk Assessment (ERA) is a strategic conceptual framework to characterize the nature and magnitude of risks, to humans and biodiversity, of the release of chemical contaminants in the environment. Several measures have been suggested to enhance the science and application of ERA, including the identification and acknowledgment of uncertainties that potentially influence the outcome of risk assessments, and the appropriate consideration of temporal scale and its linkage to assessment endpoints [1].
Baudrot & Charles [2] proposed to approach these questions by coupling toxicokinetics-toxicodynamics models, which describe the time-course of processes leading to the adverse effects of a toxicant, with Bayesian inference. TKTD models separate processes influencing an organismal internal exposure (´toxicokinetics´, i.e., the uptake, bioaccumulation, distribution, biotransformation and elimination of a toxicant) from processes leading to adverse effects and ultimately its death (´toxicodynamics´) [3]. Although species and substance specific, the mechanistic nature of TKTD models facilitates the comparison of different toxicants, species, life stages, environmental conditions and endpoints [4].
Baudrot & Charles [2] investigated the use of a Bayesian framework to assess the uncertainties surrounding the calibration of General Unified Threshold Models of Survival (a category of TKTD) with data from standard toxicity tests, and their propagation to predictions of regulatory toxicity endpoints such as LC(x,t) [the lethal concentration affecting any x% of the population at any given exposure duration of time t] and MF(x,t) [an exposure multiplication factor leading to any x% effect reduction due to the contaminant at any time t].
Once calibrated with empirical data, GUTS models were used to explore individual survival over time, and under untested exposure conditions. Lethal concentrations displayed a strong curvilinear decline with time of exposure. For a given total amount of contaminant, pulses separated by short time intervals yielded higher mortality than pulses separated by long time intervals, as did few pulses of high amplitude when compared to multiple pulses of low amplitude. The response to a pulsed contaminant exposure was strongly influenced by contaminant depuration times. These findings highlight one important contribution of TKTD modelling in ecotoxicology: they represent just a few of the hundreds of exposure scenarios that could be mathematically explored, and that would be unfeasible or even unethical to conduct experimentally.
GUTS models were also used for interpolations or extrapolations of assessment endpoints, and their marginal distributions. A case in point is the incipient lethal concentration. The responses of model organisms to contaminants in standard toxicity tests are typically assessed at fixed times of exposure (e.g. 24h or 48h in the Daphnia magna acute toxicity test). However, because lethal concentrations are strongly time-dependent, it has been suggested that a more meaningful endpoint would be the incipient (i.e. asymptotic) lethal concentration when time of exposure increases to infinity. The authors present a mathematical solution for calculating the marginal distribution of such incipient lethal concentration, thereby providing both more relevant information and a way of comparing experiments, compounds or species tested for different periods of time.
Uncertainties were found to change drastically with time of exposure, being maximal at extreme values of x for both LC(x,t) and MF(x,t). In practice this means that assessment endpoints estimated when the effects of the contaminant are weak (such as LC10, the contaminant concentration resulting in the mortality of 10% of the experimental population), a commonly used assessment value in ERA, are prone to be highly variable.
The authors end with recommendations for improved experimental design, including (i) using assessment endpoints at intermediate values of x (e.g., LC50 instead of LC10) (ii) prolonging exposure and recording mortality over the course of the experiment (iii) experimenting one or few peaks of high amplitude close to each other when assessing pulsed exposure. Whereas these recommendations are not that different from current practices, they are based on a more coherent mechanistic grounding.
Overall, this and other contributions from Charles, Baudrot and their research group contribute to turn TKTD models into a real tool for Environmental Risk Assessment. Further enhancement of ERA´s science and application could be achieved by extending the use of TKTD models to sublethal rather than lethal effects, and to chronic rather than acute exposure, as these are more controversial issues in decision-making regarding contaminated sites.

References

[1] Dale, V. H., Biddinger, G. R., Newman, M. C., Oris, J. T., Suter, G. W., Thompson, T., ... & Chapman, P. M. (2008). Enhancing the ecological risk assessment process. Integrated environmental assessment and management, 4(3), 306-313. doi: 10.1897/IEAM_2007-066.1
[2] Baudrot, V., & Charles, S. (2018). Recommendations to address uncertainties in environmental risk assessment using toxicokinetics-toxicodynamics models. bioRxiv, 356469, ver. 3 peer-reviewed and recommended by PCI Ecol. doi: 10.1101/356469
[3] EFSA Panel on Plant Protection Products and their Residues (PPR), Ockleford, C., Adriaanse, P., Berny, P., Brock, T., Duquesne, S., ... & Kuhl, T. (2018). Scientific Opinion on the state of the art of Toxicokinetic/Toxicodynamic (TKTD) effect models for regulatory risk assessment of pesticides for aquatic organisms. EFSA Journal, 16(8), e05377. doi: 10.2903/j.efsa.2018.5377
[4] Jager, T., Albert, C., Preuss, T. G., & Ashauer, R. (2011). General unified threshold model of survival-a toxicokinetic-toxicodynamic framework for ecotoxicology. Environmental science & technology, 45(7), 2529-2540. doi: 10.1021/es103092a

Recommendations to address uncertainties in environmental risk assessment using toxicokinetics-toxicodynamics modelsVirgile Baudrot and Sandrine Charles<p>Providing reliable environmental quality standards (EQS) is a challenging issue for environmental risk assessment (ERA). These EQS are derived from toxicity endpoints estimated from dose-response models to identify and characterize the environm...Chemical ecology, Ecotoxicology, Experimental ecology, Statistical ecologyLuis Schiesari2018-06-27 21:33:30 View
27 Nov 2023
article picture

Modeling Tick Populations: An Ecological Test Case for Gradient Boosted Trees

Gradient Boosted Trees can deliver more than accurate ecological predictions

Recommended by ORCID_LOGO based on reviews by 2 anonymous reviewers

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.

References

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

Modeling Tick Populations: An Ecological Test Case for Gradient Boosted TreesWilliam Manley, Tam Tran, Melissa Prusinski, Dustin Brisson<p style="text-align: justify;">General linear models have been the foundational statistical framework used to discover the ecological processes that explain the distribution and abundance of natural populations. Analyses of the rapidly expanding ...Parasitology, Species distributions, Statistical ecologyTimothée PoisotAnonymous, Anonymous2023-03-23 23:41:17 View
12 May 2022
article picture

Riparian forest restoration as sources of biodiversity and ecosystem functions in anthropogenic landscapes

Complex but positive diversity - ecosystem functioning relationships in Riparian tropical forests

Recommended by ORCID_LOGO based on reviews by 2 anonymous reviewers

Many ecological drivers can impact ecosystem functionality and multifunctionality, with the latter describing the joint impact of different functions on ecosystem performance and services. It is now generally accepted that taxonomically richer ecosystems are better able to sustain high aggregate functionality measures, like energy transfer, productivity or carbon storage (Buzhdygan 2020, Naeem et al. 2009), and different ecosystem services (Marselle et al. 2021) than those that are less rich. Antonini et al. (2022) analysed an impressive dataset on animal and plant richness of tropical riparian forests and abundances, together with data on key soil parameters. Their work highlights the importance of biodiversity on functioning, while accounting for a manifold of potentially covarying drivers. Although the key result might not come as a surprise, it is a useful contribution to the diversity - ecosystem functioning topic, because it is underpinned with data from tropical habitats. To date, most analyses have focused on temperate habitats, using data often obtained from controlled experiments. 

The paper also highlights that diversity–functioning relationships are complicated. Drivers of functionality vary from site to site and each measure of functioning, including parameters as demonstrated here, can be influenced by very different sets of predictors, often associated with taxonomic and trait diversity. Single correlative comparisons of certain aspects of diversity and functionality might therefore return very different results. Antonini et al. (2022) show that, in general, using 22 predictors of functional diversity, varying predictor subsets were positively associated with soil functioning. Correlational analyses alone cannot resolve the question of causal link. Future studies should therefore focus on inferring precise mechanisms behind the observed relationships, and the environmental constraints on predictor subset composition and strength.

References

Antonini Y, Beirão MV, Costa FV, Azevedo CS, Wojakowski MM, Kozovits AR, Pires MRS, Sousa HC de, Messias MCTB, Fujaco MA, Leite MGP, Vidigal JP, Monteiro GF, Dirzo R (2022) Riparian forest restoration as sources of biodiversity and ecosystem functions in anthropogenic landscapes. bioRxiv, 2021.09.08.459375, ver. 3 peer-reviewed and recommended by Peer Community in Ecology. https://doi.org/10.1101/2021.09.08.459375

Buzhdygan OY, Meyer ST, Weisser WW, Eisenhauer N, Ebeling A, Borrett SR, Buchmann N, Cortois R, De Deyn GB, de Kroon H, Gleixner G, Hertzog LR, Hines J, Lange M, Mommer L, Ravenek J, Scherber C, Scherer-Lorenzen M, Scheu S, Schmid B, Steinauer K, Strecker T, Tietjen B, Vogel A, Weigelt A, Petermann JS (2020) Biodiversity increases multitrophic energy use efficiency, flow and storage in grasslands. Nature Ecology & Evolution, 4, 393–405. https://doi.org/10.1038/s41559-020-1123-8

Marselle MR, Hartig T, Cox DTC, de Bell S, Knapp S, Lindley S, Triguero-Mas M, Böhning-Gaese K, Braubach M, Cook PA, de Vries S, Heintz-Buschart A, Hofmann M, Irvine KN, Kabisch N, Kolek F, Kraemer R, Markevych I, Martens D, Müller R, Nieuwenhuijsen M, Potts JM, Stadler J, Walton S, Warber SL, Bonn A (2021) Pathways linking biodiversity to human health: A conceptual framework. Environment International, 150, 106420. https://doi.org/10.1016/j.envint.2021.106420

Naeem S, Bunker DE, Hector A, Loreau M, Perrings C (Eds.) (2009) Biodiversity, Ecosystem Functioning, and Human Wellbeing: An Ecological and Economic Perspective. Oxford University Press, Oxford. https://doi.org/10.1093/acprof:oso/9780199547951.001.0001

Riparian forest restoration as sources of biodiversity and ecosystem functions in anthropogenic landscapesYasmine Antonini, Marina Vale Beirao, Fernanda Vieira Costa, Cristiano Schetini Azevedo, Maria Wojakowski, Alessandra Kozovits, Maria Rita Silverio Pires, Hildeberto Caldas Sousa, Maria Cristina Teixeira Braga Messias, Maria Augusta Goncalves Fuja...<ol> <li style="text-align: justify;">Restoration of tropical riparian forests is challenging, since these ecosystems are the most diverse, dynamic, and complex physical and biological terrestrial habitats. This study tested whether biodiversity ...Biodiversity, Community ecology, Ecological successions, Ecosystem functioning, Terrestrial ecologyWerner Ulrich2021-09-10 10:51:23 View