Mutualistic interactions are the weird uncles of population and community ecology. They are everywhere, from the microbes aiding digestion in animals’ guts to animal-pollination services in ecosystems; They increase productivity through facilitation; They fascinate us when small birds pick the teeth of a big-mouthed crocodile. Yet, mutualistic interactions are far less studied and understood than competition or predation. Possibly because we are naively convinced that there is no mystery here: isn’t it obvious that mutualistic interactions necessarily facilitate species coexistence? Since mutualistic species benefit from one another, if one species evolves, the other should just follow, isn’t that so?

It is not as simple as that, for several reasons. First, because simple mutualistic Lotka-Volterra models showed that most of the time mutualistic systems should drift to infinity and be unstable (e.g. Goh 1979). This is not what happens in natural populations, so something is missing in simple models. At a larger scale, that of communities, this is even worse, since we are still far from understanding the link between the topology of mutualistic networks and the stability of a community. Second, interactions are context-dependent: mutualistic species exchange resources, and thus from the point of view of one species the interaction is either beneficial or not, depending on the net gain of energy (e.g. Holland and DeAngelis 2010). In other words, considering interactions as mutualistic per se is too caricatural. Third, since evolution is blind, the evolutionary response of a species to an environmental change can have any effect on its mutualistic partner, and not necessarily a neutral or positive effect. This latter reason is particularly highlighted by the paper by A. Weinbach et al. (2021).

Weinbach et al. considered a simple two-species mutualistic Lotka-Volterra model and analyzed the evolutionary dynamics of a trait controlling for the rate of interaction between the two species by using the classical Adaptive Dynamics framework. They showed that, depending on the form of the trade-off between this interaction trait and its effect on the intrinsic growth rate, several situations can occur at evolutionary equilibrium: species can stably coexist and maintain their interaction, or the interaction traits can evolve to zero where species can coexist without any interactions.

Weinbach et al. then investigated the fate of the two-species system if a partner species is strongly affected by environmental change, for instance, a large decrease of its growth rate. Because of the supposed trade-off between the interaction trait and the growth rate, the interaction trait in the focal species tends to decrease as an evolutionary response to the decline of the partner species. If environmental change is too large, the interaction trait can evolve to zero and can lead the partner species to extinction. An “evolutionary murder”.

Even though Weinbach et al. interpreted the results of their model through the lens of plant-pollinators systems, their model is not specific to this case. On the contrary, it is very general, which has advantages and caveats. By its generality, the model is informative because it is a proof of concept that the evolution of mutualistic interactions can have unexpected effects on any category of mutualistic systems. Yet, since the model lacks many specificities of plant-pollinator interactions, it is hard to evaluate how their result would apply to plant-pollinators communities.

I wanted to recommend this paper as a reminder that it is certainly worth studying the evolution of mutualistic interactions, because i) some unexpected phenomenons can occur, ii) we are certainly too naive about the evolution and ecology of mutualistic interactions, and iii) one can wonder to what extent we will be able to explain the stability of mutualistic communities without accounting for the co-evolutionary dynamics of mutualistic species.

References

Goh BS (1979) Stability in Models of Mutualism. The American Naturalist, 113, 261–275. http://www.jstor.org/stable/2460204.

Holland JN, DeAngelis DL (2010) A consumer–resource approach to the density-dependent population dynamics of mutualism. Ecology, 91, 1286–1295. https://doi.org/10.1890/09-1163.1

Weinbach A, Loeuille N, Rohr RP (2021) Eco-evolutionary dynamics further weakens mutualistic interaction and coexistence under population decline. bioRxiv, 570580, ver. 5 peer-reviewed and recommended by Peer Community in Ecology. https://doi.org/10.1101/570580

Climate change is accelerating (IPCC 2022), and so applies ever stronger selective pressures on biodiversity (Segan et al. 2016). Possible responses include range shifts or adaptations to new climatic conditions (Bellard et al. 2012), but there is still much uncertainty about the extent of most species' adaptive capacities and the impact of extreme climatic events.
 
The pine processionary is a major pest of pine trees in the Mediterranean area. It is notably one of the few species for which a clear link between recent climate change and its northward expansion has been established (Battisti et al. 2005), and as such is often considered as globally benefitting from climate change. However, recent results show a retraction of its range at the southern limit (Bourougaaoui et al. 2021), exposed to high warming (+1.4°C in Tunisia since 1901 as opposed to +1.12°C on average in the Northern hemisphere) and extreme summer temperature events (Verner et al. 2013). Thus, it is possible that the species' adaptive abilities are being challenged at the southern limit of its native range by the magnitude of observed climate change.
 
In this work, Bourougaaoui et al. (2022) investigate how climate change over the last 30 years has impacted the reproductive success of the pine processionary moth in Tunisia. A major methodological interest of this study is that they used data both from historical collections and from recent samplings, which raised a challenge for running a longitudinal analysis as sampling locations differed between the two periods. By applying a grouping method to local climatic data, the authors were able to define several large climatic clusters within the country, and analyze long-term data from different sites within the same clusters. They find that both fecundity and hatching rate decreased over the period, while at the same time both the average temperature increased and climate variability increased. One of the main conclusions is that recurrent episodes of extreme heat during summer might have a larger impact than the long-term increase of average temperature, which strongly echoes how the intensification of weather extremes is currently proving one of the most important dimensions of climate change.
 
However, a most interesting hypothesis also arises from the analysis of the differences between climatic clusters: preexisting adaptations to heat, for instance, phenological shifts that allow the most sensitive stages to develop earlier in the season before the extreme heat events are most likely to occur, might actually reduce impacts in the historically warmest areas. Thus the greatest climate vulnerability might not always stand where one expects it.
 
References

Battisti A, Stastny M, Netherer S, Robinet C, Schopf A, Roques A, Larsson S (2005) Expansion of Geographic Range in the Pine Processionary Moth Caused by Increased Winter Temperatures. Ecological Applications, 15, 2084–2096. https://doi.org/10.1890/04-1903

Bellard C, Bertelsmeier C, Leadley P, Thuiller W, Courchamp F (2012) Impacts of climate change on the future of biodiversity. Ecology Letters, 15, 365–377. https://doi.org/10.1111/j.1461-0248.2011.01736.x

Bourougaaoui A, Ben Jamâa ML, Robinet C (2021) Has North Africa turned too warm for a Mediterranean forest pest because of climate change? Climatic Change, 165, 46. https://doi.org/10.1007/s10584-021-03077-1

Bourougaaoui A, Robinet C, Jamaa MLB, Laparie M (2022) Effects of climate warming on the pine processionary moth at the southern edge of its range: a retrospective analysis on egg survival in Tunisia. bioRxiv, 2021.08.17.456665, ver. 5 peer-reviewed and recommended by Peer Community in Ecology. https://doi.org/10.1101/2021.08.17.456665

IPCC. 2022. Climate Change 2022: Impacts, Adaptation, and Vulnerability. Contribution of Working Group II to the Sixth Assessment Report of the Intergovernmental Panel on Climate Change [H.-O. Pörtner, D.C. Roberts, M. Tignor, E.S. Poloczanska, K. Mintenbeck, A. Alegría, M. Craig, S. Langsdorf, S. Löschke, V. Möller, A. Okem, B. Rama (eds.)]. Cambridge University Press. In Press.

Segan DB, Murray KA, Watson JEM (2016) A global assessment of current and future biodiversity vulnerability to habitat loss–climate change interactions. Global Ecology and Conservation, 5, 12–21. https://doi.org/10.1016/j.gecco.2015.11.002

Verner D (2013) Tunisia in a Changing Climate : Assessment and Actions for Increased Resilience and Development. World Bank, Washington, DC. https://doi.org/10.1596/978-0-8213-9857-9  

“Pattern, like beauty, is to some extent in the eye of the beholder” (Grant 1977 in Wiens, 1989)

The preprint "Environmental perturbations and transitions between ecological and evolutionary equilibria: an eco-evolutionary feedback framework" by Coulson (2020) presents a general framework for evolutionary ecology, useful to interpret patterns of selection and evolutionary responses to environmental transitions. The paper is written in an accessible and intuitive manner. It reviews important concepts which are at the heart of evolutionary ecology. Together, they serve as a worldview which you can carry with you to interpret patterns in data or observations in nature. I very much appreciate it that Coulson (2020) presents his personal intuition laid bare, the framework he uses for his research and how several strong concepts from theoretical ecology fit in there. Overviews as presented in this paper are important to understand how we as researchers put the pieces together.
A main message of the paper is that resource detection and acquisition traits, broadly called "resource accrual traits" are at the core of evolutionary dynamics. These traits and the processes they are involved in often urge some degree of individual specialization. Not all traits are resource accrual traits all the time. Guppies are cited as an example, which have traits in high predation environments that make foraging easier for them, such as being less conspicuous to predators. In the absence of predators, these same traits might be neutral. Their colour pattern might then contribute much less to the odds of obtaining resources.
"Resource accrual" reminds me of discussions of resource holding potential (Parker 1974), which can be for example the capacity to remain on a bird feeder without being dislodged. However, the idea is much broader and aggression does not need to be important for the acquisition of resources. Evolutionary success is reserved for those steadily obtaining resources. This recalls the pessimization principle of Metz et al. (2008), which applies in a restricted set of situations and where the strategy which persists at the lowest resource levels systematically wins evolutionary contests. If this principle would apply universally, the world then inherently become the worst possible. Resources determine energy budgets and different life history strategies allocate these differently to maximize fitness. The fine grain of environments and the filtration by individual histories generate a lot of variation in outcomes. However, constraint-centered approaches (Kempes et al. 2019, Kooijman 2010) are mentioned but are not at the core of this preprint. Evolution is rather seen as dynamic programming optimization with interactions within and between species. Coulson thus extends life history studies such as for example Tonnabel et al. (2012) with eco-evolutionary feedbacks. Examples used are guppies, algae-rotifer interactions and others. Altogether, this makes for an optimistic paper pushing back the pessimization principle.
Populations are expected to spend most of the time in quasi-equilibrium states where the long run stochastic growth rate is close to zero for all genotypes, alleles or other chosen classes. In the preprint, attention is given to reproductive value calculus, another strong tool in evolutionary dynamics (Grafen 2006, Engen et al. 2009), which tells us how classes within a population contribute to population composition in the distant future. The expected asymptotic fitness of an individual is equated to its expected reproductive value, but this might require particular ways of calculating reproductive values (Coulson 2020). Life history strategies can also be described by per generation measures such as R0 (currently on everyone's radar due to the coronavirus pandemic), generation time etc. Here I might disagree because I believe that this focus in per generation measures can lead to an incomplete characterization of plastic and other strategies involved in strategies such as bet-hedging. A property at quasi-equilibrium states is precise enough to serve as a null hypothesis which can be falsified: all types must in the long run leave equal numbers of descendants. If there is any property in evolutionary ecology which is useful it is this one and it rightfully merits attention.
However, at quasi-equilibrium states, directional selection has been observed, often without the expected evolutionary response. The preprint aims to explain this and puts forward the presence of non-additive gene action as a mechanism. I don't believe that it is the absence of clonal inheritance which matters very much in itself (Van Dooren 2006) unless genes with major effect are present in protected polymorphisms. The preprint remains a bit unclear on how additive gene action is broken, and here I add from the sphere in which I operate. Non-additive gene action can be linked to non-linear genotype-phenotype maps (Van Dooren 2000, Gilchrist and Nijhout 2001) and if these maps are non-linear enough to create constraints on phenotype determination, by means of maximum or minimum phenotypes which cannot be surpassed for any combination of the underlying traits, then they create additional evolutionary quasi-equilibrium states, with directional selection on a phenotype such as body size. I believe Coulson hints at this option (Coulson et al. 2006), but also at a different one: if body size is mostly determined by variation in resource accrual traits, then the resource accrual traits can be under stabilizing selection while body size is not. This requires that all resource accrual traits affect other phenotypic or demographic properties next to body size. In both cases, microevolutionary outcomes cannot be inferred from inspecting body sizes alone, either resource accrual traits need to be included explicitly, or non-linearities, or both when the map between resource accrual and body size is non-linear (Van Dooren 2000).
The discussion of the phenotypic gambit (Grafen 1984) leads to another long-standing issue in evolutionary biology. Can predictions of adaptation be made by inspecting and modelling individual phenotypes alone? I agree that with strongly non-linear genotype-phenotype maps they cannot and for multivariate sets of traits, genetic and phenotypic correlations can be very different (Hadfield et al. 2007). However, has the phenotypic gambit ever claimed to be valid globally or should it rather be used locally for relatively small amounts of variation? Grafen (1984) already contained caveats which are repeated here. As a first approximation, additivity might produce quite correct predictions and thus make the gambit operational in many instances. When important individual traits are omitted, it may just be misspecified. I am interested to see cases where the framework Coulson (2020) proposes is used for very large numbers of phenotypic and genotypic attributes. In the end, these highly dimensional trait distributions might basically collapse to a few major axes of variation due to constraints on resource accrual.
I highly recommend reading this preprint and I am looking forward to the discussion it will generate.

References

[1] Coulson, T. (2020) Environmental perturbations and transitions between ecological and evolutionary equilibria: an eco-evolutionary feedback framework. bioRxiv, 509067, ver. 4 peer-reviewed and recommended by PCI Ecology. doi: 10.1101/509067
[2] Coulson, T., Benton, T. G., Lundberg, P., Dall, S. R. X., and Kendall, B. E. (2006). Putting evolutionary biology back in the ecological theatre: a demographic framework mapping genes to communities. Evolutionary Ecology Research, 8(7), 1155-1171.
[3] Engen, S., Lande, R., Sæther, B. E. and Dobson, F. S. (2009) Reproductive value and the stochastic demography of age-structured populations. The American Naturalist 174: 795-804. doi: 10.1086/647930
[4] Gilchrist, M. A. and Nijhout, H. F. (2001). Nonlinear developmental processes as sources of dominance. Genetics, 159(1), 423-432.
[5] Grafen, A. (1984) Natural selection, kin selection and group selection. In: Behavioural Ecology: An Evolutionary Approach,2nd edn (JR Krebs & NB Davies eds), pp. 62–84. Blackwell Scientific, Oxford.
[6] Grafen, A. (2006). A theory of Fisher's reproductive value. Journal of mathematical biology, 53(1), 15-60. doi: 10.1007/s00285-006-0376-4
[7] Hadfield, J. D., Nutall, A., Osorio, D. and Owens, I. P. F. (2007). Testing the phenotypic gambit: phenotypic, genetic and environmental correlations of colour. Journal of evolutionary biology, 20(2), 549-557. doi: 10.1111/j.1420-9101.2006.01262.x
[8] Kempes, C. P., West, G. B., and Koehl, M. (2019). The scales that limit: the physical boundaries of evolution. Frontiers in Ecology and Evolution, 7, 242. doi: 10.3389/fevo.2019.00242
[9] Kooijman, S. A. L. M. (2010) Dynamic Energy Budget theory for metabolic organisation. University Press, third edition.
[10] Metz, J. A. J., Mylius, S.D. and Diekman, O. (2008) When does evolution optimize?. Evolutionary Ecology Research 10: 629-654.
[11] Parker, G. A. (1974). Assessment strategy and the evolution of fighting behaviour. Journal of theoretical Biology, 47(1), 223-243. doi: 10.1016/0022-5193(74)90111-8
[12] Tonnabel, J., Van Dooren, T. J. M., Midgley, J., Haccou, P., Mignot, A., Ronce, O., and Olivieri, I. (2012). Optimal resource allocation in a serotinous non‐resprouting plant species under different fire regimes. Journal of Ecology, 100(6), 1464-1474. doi: 10.1111/j.1365-2745.2012.02023.x
[13] Van Dooren, T. J. M. (2000). The evolutionary dynamics of direct phenotypic overdominance: emergence possible, loss probable. Evolution, 54(6), 1899-1914. doi: 10.1111/j.0014-3820.2000.tb01236.x
[14] Van Dooren, T. J. M. (2006). Protected polymorphism and evolutionary stability in pleiotropic models with trait‐specific dominance. Evolution, 60(10), 1991-2003. doi: 10.1111/j.0014-3820.2006.tb01837.x

In this study, the authors report a new method for estimating the abundance of the Pyrenean brown bear population. Precisely, the methodology involved aims to apply Pollock's closed robust design (PCRD) capture-recapture models to estimate population abundance and trends over time. Overall, the results encourage the use of PCRD to study populations' demographic rates, while minimizing biases due to inter-individual heterogeneity in detection probabilities.

Estimating the size and trends of animal population over time is essential for informing conservation status and management decision-making (Nichols & Williams 2006). This is particularly the case when the population is small, geographically scattered, and threatened. Although several methods can be used to estimate population abundance, they may be difficult to implement when individuals are rare, elusive, solitary, largely nocturnal, highly mobile, and/or occupy large home ranges in remote and/or rugged habitats. Moreover, in such standard methods,

• the population is assumed to be closed both geographically (no immigration nor emigration) and demographically (no births nor deaths) and
• all individuals are assumed to have identical detection probabilities regardless of their individual attributes (e.g., age, body mass, social status) and habitat features (home-range location and composition) (Otis et al. 1978).

However, these conditions are rarely met in real populations, such as wild mammals (e.g., Bellemain et al. 2005; Solbert et al. 2006), and therefore the risk of underestimating population size can rapidly increase because the assumption of perfect detection of all individuals in the population is violated.

Focusing on the critically endangered Pyrenean brown bear that was close to extinction in the mid-1990s, the study by Vanpe et al. (2022), uses protocol and opportunistic data to describe a statistical modeling exercise to construct mark-recapture histories from 2008 to 2020. Among the data, the authors collected non-invasive samples such as a mixture of hair and scat samples used for genetic identification, as well as photographic trap data of recognized individuals. These data are then analyzed in RMark to provide detection and survival estimates. The final model (i.e. PCRD capture-recapture) is then used to provide Bayesian population estimates. Results show a five-fold increase in population size between 2008 and 2020, from 13 to 66 individuals. Thus, this study represents the first published annual abundance and temporal trend estimates of the Pyrenean brown bear population since 2008.

Then, although the results emphasize that the PCRD estimates were broadly close to the MRS counts and had reasonably narrow associated 95% Credibility Intervals, they also highlight that the sampling effort is different according to individuals. Indeed, as expected, the detection of an individual depends on

• the intraspecific home range size variation that results in individuals that move the most being most likely to be detected and
• the mortality rate which is higher on cubs than on adults and subadults (due to infanticide by males, predation, death of the mother, or abandonment).

Overall, the PCRD capture-recapture modelling approach, involved in this study, provides robust estimates of abundance and demographic rates of the Pyrenean brown bear population (with associated uncertainty) while minimizing and considering bias due to inter-individual heterogeneity in detection probabilities.

The authors conclude that mark-recapture provides useful population estimates and urge wildlife ecologists and managers to use robust approaches, such as the RDPC capture-recapture model, when studying large mammal populations. This information is essential to inform management decisions and assess the conservation status of populations.

References

Bellemain, E.V.A., Swenson, J.E., Tallmon, D., Brunberg, S. and Taberlet, P. (2005). Estimating population size of elusive animals with DNA from hunter-collected feces: four methods for brown bears. Cons. Biol. 19(1), 150-161. https://doi.org/10.1111/j.1523-1739.2005.00549.x

Nichols, J.D. and Williams, B.K. (2006). Monitoring for conservation. Trends Ecol. Evol. 21(12), 668-673. https://doi.org/10.1016/j.tree.2006.08.007

Otis, D.L., Burnham, K.P., White, G.C. and Anderson, D.R. (1978). Statistical inference from capture data on closed animal populations. Wildlife Monographs (62), 3-135.

Solberg, K.H., Bellemain, E., Drageset, O.M., Taberlet, P. and Swenson, J.E. (2006). An evaluation of field and non-invasive genetic methods to estimate brown bear (Ursus arctos) population size. Biol. Conserv. 128(2), 158-168. https://doi.org/10.1016/j.biocon.2005.09.025

Vanpé C, Piédallu B, Quenette P-Y, Sentilles J, Queney G, Palazón S, Jordana IA, Jato R, Elósegui Irurtia MM, de la Torre JS, and Gimenez O (2022) Estimating abundance of a recovering transboundary brown bear population with capture-recapture models. bioRxiv, 2021.12.08.471719, ver. 4 recommended and peer-reviewed by PCI Ecology. https://doi.org/10.1101/2021.12.08.471719