How do you estimate the biodiversity of a whole community, or the distribution of abundances and ranges of its species, from presence/absence data in scattered samples?
It all starts with the collector's dilemma: if you double the number of samples, you will not get double the number of species, since you will find many of the same common species, and only a few new rare ones.
This non-additivity has prompted many ecologists to study the Species-Area Relationship. A common theoretical approach has been to connect this spatial pattern to the overall distribution of how common or rare a species can be. At least since Fisher's celebrated log-series [1], ecologists have been trying to, first, infer the shape of the Species Abundance Distribution, and then, use it to predict how many species should be found in a given area or a given number of samples. This has found many applications, from microbial communities to tropical forests, from estimating the number of yet-unknown species to predicting how much biodiversity may be lost if a fraction of the habitat is removed.
In this elegant work, Tovo et al. [2] propose a method that starts only from presence/absence data over a number of samples, and provides the community's diversity, as well as its abundance and range size distributions. This method is simple, analytically explicit, and accurate: the authors test it on the classic Pasoh and Barro Colorado Island tropical forest datasets, and on simulated data. They make a very laudable effort in both explaining its theoretical underpinnings, and proposing a straightforward step-by-step guide to applying it to data.
The core of Tovo et al's method is a simple property: the scale invariance of the Negative Binomial (NB) distribution. Subsampling from a NB gives another NB, where a single parameter has changed. Therefore, if the Species Abundance Distribution is close enough to some NB (which is flexible enough to accommodate all the data here), we can estimate how this parameter changes when going from (1) a single sample to (2) all the available samples, and from there, extrapolate to (3) the entire community.
This principle was first applied by the authors in a previous study [3] that required abundance data in the samples, rather than just presence/absence. Given that binary occurrence data is far more available in a variety of empirical settings, this extension is worthwhile (including its new predictions on range size distributions), and it deserves to be widely known and tested.

ADDITIONAL COMMENTS

1) To explain the novelty of the authors' contribution, it is useful to look at competing techniques.
Some ""parametric"" approaches try to infer the whole-community Species Abundance Distribution (SAD) by guessing its functional form (Gaussian, power-law, log-series...) and fitting its parameters from sampled data. The issue is that this distribution shape may not remain in the same family as we increase the sampling effort or area, so the regression problem may not be well-defined. This is where the Negative Binomial's scale invariance is useful.
Other ""non-parametric"" approaches have renounced guessing the whole SAD: they simply try to approximate of its tail of rare species, by looking at how many species are found in only one (or a few) samples. From this, they derive an estimate of biodiversity that is agnostic to the rest of the SAD. Tovo et al. [2] show the issue with these approaches: they extrapolate from the properties of individual samples to the whole community, but do not properly account for the bias introduced by the amount of sampling (the intermediate scale (2) in the summary above).

2) The main condition for all such approaches to work is well-mixedness: each sample should be sufficiently like a lot drawn from the same skewed lottery. As long as that condition applies, finding the best approach is a theoretical matter of probabilities and combinatorics that may, in time, be given a definite answer.
The authors also show that ""well-mixed"" is not as restrictive as it sounds: the method works both on real data (which is never perfectly mixed) and on simulations where species are even more spatially clustered than the empirical data. In addition, the Negative Binomial's scale invariance entails that, if it works well enough at some spatial scale, it will also work at all higher scales (until one reaches the edges of the sufficiently-well-mixed community)

3) One may ask: why the Negative Binomial as a Species Abundance Distribution?
If one wishes for some dynamical explanation, the Negative Binomial can be derived from neutral birth and death process with immigration, as shown by the authors in [3]. But to be applied to data, it should only be able to approximate the empirical distribution well enough (at all relevant scales). Depending on one's taste, this type of probabilistic approaches can be interpreted as:
- purely phenomenological, describing only the observational process of sampling from an existing state of affairs, not the ecological processes that gave rise to that state.
- a null model, from which everything in practice is expected to deviate to some extent.
- or a way to capture the statistical forces that tend to induce stable relationships between different patterns (as long as no ecological process opposes them strongly enough).

References

[1] Fisher, R. A., Corbet, A. S., & Williams, C. B. (1943). The relation between the number of species and the number of individuals in a random sample of an animal population. The Journal of Animal Ecology, 42-58. doi: 10.2307/1411
[2] Tovo, A., Formentin, M., Suweis, S., Stivanello, S., Azaele, S., & Maritan, A. (2019). Inferring macro-ecological patterns from local species' occurrences. bioRxiv, 387456, ver. 2 peer-reviewed and recommended by PCI Ecol. doi: 10.1101/387456
[3] Tovo, A., Suweis, S., Formentin, M., Favretti, M., Volkov, I., Banavar, J. R., Azaele, S., & Maritan, A. (2017). Upscaling species richness and abundances in tropical forests. Science Advances, 3(10), e1701438. doi: 10.1126/sciadv.1701438

In social species conflicts among group members typically lead to the formation of dominance hierarchies with dominant individuals outcompeting other groups members and, in some extreme cases, suppressing reproduction of subordinates. It has therefore been typically assumed that dominant individuals have a higher breeding success than subordinates. However, previous work on mammals (mostly primates) revealed high variation, with some populations showing no evidence for a link between female dominance reproductive success, and a meta-analysis on primates suggests that the strength of this relationship is stronger for species with a longer lifespan [1]. Therefore, there is now a need to understand 1) whether dominance and reproductive success are generally associated across social mammals (and beyond) and 2) which factors explains the variation in the strength (and possibly direction) of this relationship.
In their preregistration, Shivani et al. [2] plan to perform a meta-analysis on 86 social mammal species to address these two points. More specifically, they will investigate whether the relationship between female dominance and reproductive success vary according to life history traits (e.g. stronger for species with large litter size), ecological conditions (e.g. stronger when resources are limited) and the social environment (e.g. stronger for cooperative breeders than for plural breeders).
The two reviewers and I were particularly positive and enthusiastic about this preregistration and only had minor comments that were nicely addressed by the authors. We found the background well-grounded in the existing literature and that the predictions were therefore clear and well-motivated. The methods were particularly transparent with a nicely annotated R script and the authors even simulated a dataset with the same structure as the actual data in order to make sure that the coding of the data handling and statistical analyses were appropriate (without being tempted to look at model outputs from the true dataset).
Perhaps one limitation to keep in mind once we will have the chance to look at the outcome of this study if that the dataset may not be fully representative of social species with dominance hierarchies. For example, the current dataset contains only one aquatic mammal (Mirounga angustirostris) as far as I can see, which is likely due to a lack of knowledge on such systems. Furthermore, not only mammals exhibit dominance hierarchies and it will be interesting to see if the results of the proposed study hold for other social taxa (and if not, what may explain their differences).
That being said, the proposed study will already offer a much broader overview of the relationship between dominance and reproductive success in animal societies and a better understanding for its variation. The reviewers and I believe it will make an important contribution to the fields of socio-ecology and evolutionary ecology. I therefore strongly recommend this preregistration and we are particularly looking forward to seeing the outcome of this exciting study.

References

[1] Majolo, B., Lehmann, J., de Bortoli Vizioli, A., & Schino, G. (2012). Fitness‐related benefits of dominance in primates. American journal of physical anthropology, 147(4), 652-660. doi: 10.1002/ajpa.22031
[2] Shivani, Huchard, E., Lukas, D. (2020). Preregistration - The effect of dominance rank on female reproductive success in social mammals In principle acceptance by PCI Ecology of the version 1.2 on 07 July 2020. https://github.com/dieterlukas/FemaleDominanceReproductionMetaAnalysis/blob/trunk/PreregistrationMetaAnalysis_RankSuccess.Rmd

Sexual coercion can be defined as the use by a male of force, or threat of force, which increases the chances that a female will mate with him at a time when she is likely to be fertile, and/or decrease the chances that she will mate with other males, at some cost to the female (Smuts & Smuts 1993). It has been evidenced in a wide range of species and may play an important role in the evolution of sexual conflict and social systems. However, identifying sexual coercion in natural systems can be particularly challenging. Notably, while male behaviour may have immediate consequences on mating success (“harassment”), the mating benefits may be delayed in time (“intimidation”), and in such cases, evidencing coercion requires detailed temporal data at the individual level. Moreover, in some species male aggressive behaviours may be subtle or rare and hence hardly observed, yet still have important effects on female mating probability and fitness. Therefore, investigating the occurrence and consequences of sexual coercion in such species is particularly relevant but studying it in a statistically robust way is likely to require a considerable amount of time spent observing individuals.

In this paper, Smit et al. (2022) test three clear predictions of the sexual coercion hypothesis in a natural population of Mandrills, where severe male aggression towards females is rare: (1) male aggression is more likely on sexually receptive females than on females in other reproductive states, (2) receptive females are more likely to be injured and (3) male aggression directed towards females is positively related to subsequent probability of copulation between those dyads. They also tested an alternative hypothesis, the “aggressive male phenotype” under which the correlation between male aggression towards females and subsequent mating could be statistically explained by male overall aggressivity. In agreement with the three predictions of the sexual coercion hypothesis, (1) male aggression was on average 5 times more likely, and (2) injuries twice as likely, to be observed on sexually receptive females than on females in other reproductive states and (3) copulation between males and sexually receptive females was twice more likely to be observed when aggression by this male was observed on the female before sexual receptivity. There was no support for the aggressive male hypothesis.

The reviewers and I were highly positive about this study, notably regarding the way it is written and how the predictions are carefully and clearly stated, tested, interpreted, and discussed.

This study is a good illustration of a case where some behaviours may not be common or obvious yet have strong effects and likely important consequences and thus be clearly worth studying. More generally, it shows once more the importance of detailed long-term studies at the individual level for our understanding of the ecology and evolution of wild populations.

It is also a good illustration of the challenges faced, when comparing the likelihood of contrasting hypotheses means we need to alter sample sizes and/or the likelihood to observe at all some behaviours. For example, observing copulation within minutes after aggression (and therefore, showing statistical support for “harassment”) is inevitably less likely than observing copulations on the longer-term (and therefore showing statistical support for “intimidation”, when of course effort is put into recording such behavioural data on the long-term). Such challenges might partly explain some apparently intriguing results. For example, why are swollen females more aggressed by males if only aggression before the swollen period seems associated with more chances of mating? Here, the authors systematically provide effect sizes (and confidence intervals) and often describe the effects in an intuitive biological way (e.g., “Swollen females were, on average, about five times more likely to become injured”). This clearly helps the reader to not merely compare statistical significances but also the biological strengths of the estimated effects and the uncertainty around them. They also clearly acknowledge limits due to sample size when testing the harassment hypothesis, yet they provide precious information on the probability of observing mating (a rare behaviour) directly after aggression (already a rare behaviour!), that is, 3 times out of 38 aggressions observed between a male and a swollen female. Once again, this highlights how important it is to be able to pursue the enormous effort put so far into closely and continuously monitoring this wild population.

Finally, this study raises exciting new questions, notably regarding to what extent females exhibit “counter-strategies” in response to sexual coercion, notably whether there is still scope for female mate choice under such conditions, and what are the fitness consequences of these dynamic conflicting sexual interactions. No doubt these questions will sooner than later be addressed by the authors, and I am looking forward to reading their upcoming work.

References

Smit N, Baniel A, Roura-Torres B, Amblard-Rambert P, Charpentier MJE, Huchard E (2022) Sexual coercion in a natural mandrill population. bioRxiv, 2022.02.07.479393, ver. 5 peer-reviewed and recommended by Peer Community in Ecology. https://doi.org/10.1101/2022.02.07.479393

Smuts BB, Smuts R w. (1993) Male Aggression and Sexual Coercion of Females in Nonhuman Primates and Other Mammals: Evidence and Theoretical Implications. In: Advances in the Study of Behavior (eds Slater PJB, Rosenblatt JS, Snowdon CT, Milinski M), pp. 1–63. Academic Press. https://doi.org/10.1016/S0065-3454(08)60404-0

Accurately estimating survival rate and identifying the drivers of its variation is essential for our understanding of population dynamics and life history strategies (Sæther and Bakke 2000), as well as for population management and conservation (Francis et al. 1998, Doherty et al. 2014). Many studies estimate survival from capture–recapture data using the Cormack–Jolly–Seber (CJS) model (Lebreton et al. 1992). However, survival estimates are confounded with permanent emigration from the study area, which can be particularly problematic for mobile species. This is problematic, not only because CJS models under estimate true survival in populations where permanent emigration occurs (i.e. they estimate “apparent” survival), but also because some factors of interest may affect both survival and emigration (e.g., habitat quality, Paquet et al. 2020), leaving the interpretation of results challenging, for example in terms of management decisions.

Several methods have been developed to account for permanent emigration when estimating survival, for example by jointly analyzing CMR data with data on individuals’ locations at each capture/resighting site (to estimate their dispersal distances; Schaub and Royle 2013, Badia Boher et al. 2023), with telemetry data (Powel et al. 2000), mark recovery data (Burnham 1993, Fay et al. 2019), or with live-resight data (Barker 1997).

The Barker joint live-recapture/live-resight (JLRLR) model can estimate survival when resight data are continuous over a long interval and from a larger area than the capture recapture data. This model becomes particularly promising with the growing collection of data from citizen science, or remote detection tools (Dzul et al. 2023). However, as pointed out by Dzul et al., this model assumes that resight probability is homogeneous across the area where individuals can move, and this assumption is likely violated for example because of non-random movements or because of non-random location of resighting sites.

In their manuscript, Dzul et al. performed a thorough simulation study to evaluate the accuracy of survival estimates from JLRLR models under various study designs regarding the location of resight sites (global, random, fixed including the capture site, and fixed excluding the capture site). They simulated data with varying survival and movement values, varying recapture and resight probabilities, and varying sample sizes. Finally, they also developed and fitted a multi state version of the JLRLR model. They show that JLRLR models performed better than CJS models. Survival estimates were still often biased (either positively or negatively) but they were less biased when sesight sites were randomly located (rather than at fixed locations), when recapture sites were included in the resighting design, and when using the multi state JLRLR model they developed.

This study highlights (multistate) JLRLR models as an alternative to CJS models one should consider when auxiliary resight data can be collected. Moreover, it shows the importance of evaluating not only model performance, but also the efficiency of alternative sampling designs before choosing one for our studies. Hopefully, this study will help the authors and other researchers making a more informed and efficient choice of model and design to estimate survival in their study populations.

References

Jaume A. Badia-Boher, Joan Real, Joan Lluís Riera, Frederic Bartumeus, Francesc Parés, Josep Maria Bas, and Antonio Hernández-Matías. Joint estimation of survival and dispersal effectively corrects the permanent emigration bias in mark-recapture analyses. (2023) Scientific reports 13, no. 1: 6970. https://doi.org/10.1038/s41598-023-32866-0 

Richard J Barker (1997) Joint modeling of live-recapture, tag-resight, and tag-recovery data. Biometrics: 666-677. https://doi.org/10.2307/2533966 

Kenneth P. Burnham (1993) Marked Individuals in the Study of Bird Populations (ed. J.D. Lebreton), pp. 199–213. Birkhäuser, Basel

Kevin E. Doherty, David E. Naugle, Jason D. Tack, Brett L. Walker, Jon M. Graham, Jeffrey L. Beck (2014) Linking conservation actions to demography: grass height explains variation in greater sage‐grouse nest survival. Wildlife biology 20, no. 6 : 320-325. https://doi.org/10.2981/wlb.00004

Maria C. Dzul, Charles B. Yackulic, William L. Kendall (2023) The importance of sampling design for unbiased estimation of survival using joint live-recapture and live resight models. arXiv, ver.3 peer-reviewed and recommended by PCI Ecology https://doi.org/10.48550/arXiv.2312.13414

Rémi Fay, Stephanie Michler, Jacques Laesser, and Michael Schaub (2019) Integrated population model reveals that kestrels breeding in nest boxes operate as a source population. Ecography 42, no. 12: 2122-2131. https://doi.org/10.1111/ecog.04559

Charles M. Francis, John R. Sauer, Jerome R. Serie (1998) Effect of restrictive harvest regulations on survival and recovery rates of American black ducks. The Journal of Wildlife Management : 1544-1557. https://doi.org/10.2307/3802021

Jean-Dominique Lebreton, Kenneth P. Burnham, Jean Clobert, David R. Anderson (1992) Modeling survival and testing biological hypotheses using marked animals: a unified approach with case studies. Ecological monographs 62.1: 67-118. https://doi.org/10.2307/2937171

Matthieu Paquet, Debora Arlt, Jonas Knape, Matthew Low, Pär Forslund, and Tomas Pärt (2020) Why we should care about movements: Using spatially explicit integrated population models to assess habitat source–sink dynamics. Journal of Animal Ecology 89, no. 12: 2922-2933. https://doi.org/10.1111/1365-2656.13357

Larkin A. Powell, Michael J. Conroy, James E. Hines, James D. Nichols, and David G. Krementz. Simultaneous use of mark-recapture and radiotelemetry to estimate survival, movement, and capture rates. (2000) The Journal of Wildlife Management : 302-313. https://doi.org/10.2307/3803003

Bernt-Erik Sæther, Øyvind Bakke (2000) Avian life history variation and contribution of demographic traits to the population growth rate. Ecology 81.3 : 642-653. https://doi.org/10.1890/0012-9658(2000)081[0642:ALHVAC]2.0.CO;2

Michael Schaub, J. Andrew Royle. Estimating true instead of apparent survival using spatial Cormack–Jolly–Seber models (2014) Methods in Ecology and Evolution 5, no. 12: 1316-1326. https://doi.org/10.1111/2041-210X.12134

Endophytic fungi are expected to be hyperdiverse in tropical forests, and here is an article exploring their diversity, hidden in Rubiaceae leaves, in two old-growth forests of Costa Rica. Humberto Castillo-González et al. not only described their diversity, but also test for the impact of leaf development stage, tissue origin, and site location. They distinguish the different fungal lineages and do identify distinct indicators, showing that specialization of endophytic fungi could be related to other factors in tropical forests.

This article is a great example of fungal ecology in the tropics, interacting at fine and large scale with a diversity of hosts. It also invites to discuss the high specialization observed in the tropics, and the ecology of old-growth forests in Costa Rica.

References

Humberto Castillo-González, Jason C. Slot, Stephanie Yarwood, Priscila Chaverri (2025) Exploring Rubiaceae fungal endophytes across contrasting tropical forests, tree tissues, and developmental stages. bioRxiv, ver.3 peer-reviewed and recommended by PCI Ecology https://doi.org/10.1101/2024.02.13.580172