Mimicry between species has long attracted the attention of scientists. Over a century ago, Bates first proposed that palatable species should gain a benefit by resembling unpalatable species (Bates 1862). Not long after, Müller suggested that there could also be a mutual advantage for two unpalatable species to mimic one another to reduce predator error (Müller 1879). These forms of mimicry, Batesian and Müllerian, are now widely studied, providing broad insights into behaviour, ecology and evolution.

Numerous taxa, including both invertebrates and vertebrates, show examples of Batesian or Müllerian mimicry. Bees and wasps provide a particularly interesting case due to the differences in defence between females and males of the same species. While both males and females may display warning colours, only females can sting and inject venom to cause pain and allow escape from predators. Therefore, males are palatable mimics and can resemble females of their own species or females of another species (dual sex-limited mimicry). This asymmetry in defence could have impacts on both population structure and community assembly, yet research into mimicry largely focuses on systems without sex differences.

Here, Boutin and colleagues (2023) use a differential equations model to explore the effect of mimicry on population structure and community assembly for sex-limited defended species. Specifically, they address three questions, 1) how do female noxiousness and sex-ratio influence the extinction risk of a single species?; 2) what is the effect of mimicry on species co-existence? and 3) how does dual sex-limited mimicry influence species co-existence? Their results reveal contexts in which populations with undefended males can persist, the benefit of Müllerian mimicry for species coexistence and that dual sex-limited mimicry can have a destabilising impact on species coexistence.

The results not only contribute to our understanding of how mimicry is maintained in natural systems but also demonstrate how changes in relative abundance or population structure of one species could impact another species. Further insight into the population and community dynamics of insects is particularly important given the current population declines (Goulson 2019; Seibold et al 2019).

References

Bates, H. W. 1862. Contributions to the insect fauna of the Amazon Valley, Lepidoptera: Heliconidae. Trans. Linn. Soc. Lond. 23:495- 566. https://doi.org/10.1111/j.1096-3642.1860.tb00146.x

Boutin, M., Costa, M., Fontaine, C., Perrard, A., Llaurens, V. 2022 Influence of sex-limited mimicry on extinction risk in Aculeata: a theoretical approach. bioRxiv, ver. 2 peer-reviewed and recommended by Peer Community in Ecology. https://doi.org/10.1101/2022.10.21.513153

Goulson, D. 2019. The insect apocalypse, and why it matters. Curr. Biol. 29: R967-R971. https://doi.org/10.1016/j.cub.2019.06.069

Müller, F. 1879. Ituna and Thyridia; a remarkable case of mimicry in butterflies. Trans. Roy. Entom. Roc. 1879:20-29.

Seibold, S., Gossner, M. M., Simons, N. K., Blüthgen, N., Müller, J., Ambarlı, D., ... & Weisser, W. W. 2019. Arthropod decline in grasslands and forests is associated with landscape-level drivers. Nature, 574: 671-674. https://doi.org/10.1038/s41586-019-1684-3

The secretion of extracellular polymeric substances (EPS) enables microorganisms to shape and interact with their environment [1]. EPS support cell adhesion and motility, offer protection from unfavorable conditions, and facilitate nutrient acquisition and transfer between microorganisms [2]. EPS production and consumption thus control the formation and structural organization of biofilms [3]. However, in marine environments, our understanding of the sources and composition of EPS is limited.
 
In this study, Hubas et al. [4] compare the carbon and nitrogen isotope ratios in EPS with the carbon isotope ratios of fatty acid biomarkers to identify the main EPS producers in intertidal sediments. The authors find pronounced differences in the diversity, composition, isotope signatures, and production/consumption dynamics of EPS between muddy and sandy environments. While the contribution of diatoms was highest in the bound fraction of EPS in muddy environments, diatom contribution was highest in the colloidal fraction of EPS in sandy environments. These differences between sites likely reflect the functional differences in EPS dynamics of epipelic and episammic sediment communities.
 
Taken together, the innovative approach of the authors provides insights into the diversity and origin of EPS in microphytobenthic communities and highlights the importance of different microbial groups in EPS production. These findings are vital for understanding EPS dynamics in microbial interactions and their role in the functioning of coastal ecosystems.

References

1. Flemming, H.-C. (2016) EPS-then and now. Microorganisms 4, 41 https://doi.org/10.3390/microorganisms4040041
2. Wolfaardt, G.M. et al. (1999) Function of EPS. In Microbial Extracellular Polymeric Substances, pp. 171–200, Springer Berlin Heidelberg https://doi.org/10.1007/978-3-642-60147-7
3. Flemming, H.-C. et al. (2007) The EPS matrix: the “house of biofilm cells.” J. Bacteriol. 189, 7945–7947 https://doi.org/10.1128/jb.00858-07
4. Hubas, C. et al. (2022) Identification of microbial exopolymer producers in sandy and muddy intertidal sediments by compound-specific isotope analysis. bioRxiv, ver. 2 peer-reviewed and recommended by Peer Community in Ecology. https://doi.org/10.1101/2022.12.02.516908

Tick-borne pathogens are considered as one of the major threats to public health – Lyme borreliosis being a well-known example of such disease. Global change – from climate change to changes in land use or invasive species – is playing a role in the increased risk associated with these pathogens. An important aspect of our knowledge of ticks and their associated pathogens is seasonality – one component being the phenology of within-year tick occurrences. This is important both in terms of health risk – e.g., when is the risk of encountering ticks high – and ecological understanding, as tick dynamics may depend on the weather as well as different hosts with their own dynamics and habitat use.

Hoch et al. (2023) provide a detailed data set on the phenology of one species of tick, Ixodes ricinus, in 6 different locations of France. Whereas relatively cool sites showed a clear peak in spring-summer, warmer sites showed in addition relatively high occurrences in fall-winter, with a minimum in late summer-early fall. Such results add robust data to the existing evidence of the importance of local climatic patterns for explaining tick phenology.

Recent analyses have shown that the phenology of Lyme borreliosis has been changing in northern Europe in the last 25 years, with seasonal peaks in cases occurring now 6 weeks earlier (Goren et al. 2023). The study by Hoch et al. (2023) is of too short duration to establish temporal changes in phenology (“only” 8 years, 2014-2021, see also Wongnak et al 2021 for some additional analyses; given the high year-to-year variability in weather, establishing phenological changes often require longer time series), and further work is needed to get better estimates of these changes and relate them to climate, land-use, and host density changes. Moreover, the phenology of ticks may also be related to species-specific tick phenology, and different tick species do not respond to current changes in identical ways (see for example differences between the two Ixodes species in Finland; Uusitalo et al. 2022). An efficient surveillance system requires therefore an adaptive monitoring design with regard to the tick species as well as the evolving causes of changes.

References

Goren, A., Viljugrein, H., Rivrud, I. M., Jore, S., Bakka, H., Vindenes, Y., & Mysterud, A. (2023). The emergence and shift in seasonality of Lyme borreliosis in Northern Europe. Proceedings of the Royal Society B: Biological Sciences, 290(1993), 20222420. https://doi.org/10.1098/rspb.2022.2420

Hoch, T., Madouasse, A., Jacquot, M., Wongnak, P., Beugnet, F., Bournez, L., . . . Agoulon, A. (2023). Seasonality of host-seeking Ixodes ricinus nymph abundance in relation to climate. bioRxiv, ver.4 peer-reviewed and recommended by Peer Community In Ecology. https://doi.org/10.1101/2022.07.25.501416

Uusitalo, R., Siljander, M., Lindén, A., Sormunen, J. J., Aalto, J., Hendrickx, G., . . . Vapalahti, O. (2022). Predicting habitat suitability for Ixodes ricinus and Ixodes persulcatus ticks in Finland. Parasites & Vectors, 15(1), 310. https://doi.org/10.1186/s13071-022-05410-8

Wongnak, P., Bord, S., Jacquot, M., Agoulon, A., Beugnet, F., Bournez, L., . . . Chalvet-Monfray, K. (2022). Meteorological and climatic variables predict the phenology of Ixodes ricinus nymph activity in France, accounting for habitat heterogeneity. Scientific Reports, 12(1), 7833. https://doi.org/10.1038/s41598-022-11479-z

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

Michael Raatz (2023) provides here a relevant exploration of higher-order interactions, i.e. interactions involving more than two related species (Terry et al. 2019), in the case of food web and competition interactions. More precisely, he shows by modeling that essential resources may significantly mediate focal species' persistence. Simultaneously, the provision of essential resources may strongly affect the resulting community structure, by driving to extinction first the predator and then, depending on the higher-order interaction, potentially also the associated competitor. 

Today, all ecologists should be aware of the potential effects of high-order interactions on species' (and likely on ecosystem's) fate (Golubski et al. 2016, Grilli et al. 2017). Yet, we should soon be prepared to include any high-order interaction into any interaction network (i.e. not only between species, but also between species and abiotic components, and between biotic, anthropogenic and abiotic components too). For this purpose, we will need innovative approaches such as hypergraphs (Golubski et al. 2016) and discrete-event models (Gaucherel and Pommereau 2019, Thomas et al. 2022) able to manage highly complex interactions, with numerous interacting components and variables. Such a rigorous study is a necessary and preliminary step in taking into account such a higher complexity. 

References

Gaucherel, C. and F. Pommereau. 2019. Using discrete systems to exhaustively characterize the dynamics of an integrated ecosystem. Methods in Ecology and Evolution 00:1–13. https://doi.org/10.1111/2041-210X.13242

Golubski, A. J., E. E. Westlund, J. Vandermeer, and M. Pascual. 2016. Ecological Networks over the Edge: Hypergraph Trait-Mediated Indirect Interaction (TMII) Structure trends in Ecology & Evolution 31:344-354. https://doi.org/10.1016/j.tree.2016.02.006

Grilli, J., G. Barabas, M. J. Michalska-Smith, and S. Allesina. 2017. Higher-order interactions stabilize dynamics in competitive network models. Nature 548:210-213. https://doi.org/10.1038/nature23273

Raatz, M. 2023. Provision of essential resources as a persistence strategy in food webs. bioRxiv, ver. 3 peer-reviewed and recommended by Peer Community in Ecology. https://doi.org/10.1101/2023.01.27.525839

Terry, J. C. D., R. J. Morris, and M. B. Bonsall. 2019. Interaction modifications lead to greater robustness than pairwise non-trophic effects in food webs. Journal of Animal Ecology 88:1732-1742. https://doi.org/10.1111/1365-2656.13057

Thomas, C., M. Cosme, C. Gaucherel, and F. Pommereau. 2022. Model-checking ecological state-transition graphs. PLoS Computational Biology 18:e1009657. https://doi.org/10.1371/journal.pcbi.1009657