Understanding why and how species coexist in local communities is one of the central questions in ecology. There is general agreement that species distribution and coexistence are determined by a number of key mechanisms, including the environmental requirements of species, dispersal, evolutionary constraints, resource availability and selection, metapopulation dynamics, and biotic interactions (e.g. Soberón & Nakamura 2009; Colwell & Rangel 2009; Ricklefs 2015). These factors are however intricately intertwined in a scale-structured fashion (Hortal et al. 2010; D’Amen et al. 2017), making it particularly difficult to tease apart the effects of each one of them. This could be addressed by the novel field of Joint Species Distribution Modelling (JSDM; Okasvainen & Abrego 2020), as it allows assessing the effects of several sets of factors and the co-occurrence and/or covariation in abundances of potentially interacting species at the same time (Pollock et al. 2014; Ovaskainen et al. 2016; Dormann et al. 2018). However, the development of JSDM has been hampered by the general lack of good-quality detailed data on species co-occurrences and abundances (see Hortal et al. 2015).
Facon et al. (2021) use a particularly large compilation of field surveys to study the abundance and co-occurrence of Tephritidae fruit flies in c. 400 orchards, gardens and natural areas throughout the island of Réunion. Further, they combine such information with lab data on their host-selection fundamental niche (i.e. in the absence of competitors), codifying traits of female choice and larval performances in 21 host species. They use Poisson Log-Normal models, a type of mixed model that allows one to jointly model the random effects associated with all species, and retrieve the covariations in abundance that are not explained by environmental conditions or differences in sampling effort. Then, they use a series of models to evaluate the effects on these matrices of ecological covariates (date, elevation, habitat, climate and host plant), species interactions (by comparing with a constrained residual variance-covariance matrix) and the species’ host-selection fundamental niches (through separate models for each fly species).
The eight Tephritidae species inhabiting Réunion include both generalists and specialists in Solanaceae and Cucurbitaceae with a known history of interspecific competition. Facon et al. (2021) use a comprehensive JSDM approach to assess the effects of different factors separately and altogether. This allows them to identify large effects of plant hosts and the fundamental host-selection niche on species co-occurrence, but also to show that ecological covariates and weak –though not negligible– species interactions are necessary to account for all residual variance in the matrix of joint species abundances per site. Further, they also find evidence that the fitness per host measured in the lab has a strong influence on the abundances in each host plant in the field for specialist species, but not for generalists. Indeed, the stronger effects of competitive exclusion were found in pairs of Cucurbitaceae specialist species. However, these analyses fail to provide solid grounds to assess why generalists are rarely found in Cucurbitaceae and Solanaceae. Although they argue that this may be due to Connell’s (1980) ghost of competition past (past competition that led to current niche differentiation), further data on the evolutionary history of these fruit flies is needed to assess this hypothesis.
Finding evidence for the effects of competitive interactions on species’ occurrences and spatial distributions is often difficult, perhaps because these effects occur over longer time scales than the ones usually studied by ecologists (Yackulic 2017). The work by Facon and colleagues shows that weak effects of competition can be detected also at the short ecological timescales that determine coexistence in local communities, under the virtuous combination of good-quality data and sound analytical designs that account for several aspects of species’ niches, their biotopes and their joint population responses. This adds a new dimension to the application of Hutchinson’s (1978) niche framework to understand the spatial dynamics of species and communities (see also Colwell & Rangel 2009), although further advances to incorporate dispersal-driven metacommunity dynamics (see, e.g., Ovaskainen et al. 2016; Leibold et al. 2017) are certainly needed. Nonetheless, this work shows the potential value of in-depth analyses of species coexistence based on combining good-quality field data with well-thought out JSDM applications. If many studies like this are conducted, it is likely that the uprising field of Joint Species Distribution Modelling will improve our understanding of the hierarchical relationships between the different factors affecting species coexistence in ecological communities in the near future.
Colwell RK, Rangel TF (2009) Hutchinson’s duality: The once and future niche. Proceedings of the National Academy of Sciences, 106, 19651–19658. https://doi.org/10.1073/pnas.0901650106
Connell JH (1980) Diversity and the Coevolution of Competitors, or the Ghost of Competition Past. Oikos, 35, 131–138. https://doi.org/10.2307/3544421
D’Amen M, Rahbek C, Zimmermann NE, Guisan A (2017) Spatial predictions at the community level: from current approaches to future frameworks. Biological Reviews, 92, 169–187. https://doi.org/10.1111/brv.12222
Dormann CF, Bobrowski M, Dehling DM, Harris DJ, Hartig F, Lischke H, Moretti MD, Pagel J, Pinkert S, Schleuning M, Schmidt SI, Sheppard CS, Steinbauer MJ, Zeuss D, Kraan C (2018) Biotic interactions in species distribution modelling: 10 questions to guide interpretation and avoid false conclusions. Global Ecology and Biogeography, 27, 1004–1016. https://doi.org/10.1111/geb.12759
Facon B, Hafsi A, Masselière MC de la, Robin S, Massol F, Dubart M, Chiquet J, Frago E, Chiroleu F, Duyck P-F, Ravigné V (2021) Joint species distributions reveal the combined effects of host plants, abiotic factors and species competition as drivers of community structure in fruit flies. bioRxiv, 2020.12.07.414326. ver. 4 peer-reviewed and recommended by Peer community in Ecology. https://doi.org/10.1101/2020.12.07.414326
Hortal J, de Bello F, Diniz-Filho JAF, Lewinsohn TM, Lobo JM, Ladle RJ (2015) Seven Shortfalls that Beset Large-Scale Knowledge of Biodiversity. Annual Review of Ecology, Evolution, and Systematics, 46, 523–549. https://doi.org/10.1146/annurev-ecolsys-112414-054400
Hortal J, Roura‐Pascual N, Sanders NJ, Rahbek C (2010) Understanding (insect) species distributions across spatial scales. Ecography, 33, 51–53. https://doi.org/10.1111/j.1600-0587.2009.06428.x
Hutchinson, G.E. (1978) An introduction to population biology. Yale University Press, New Haven, CT.
Leibold MA, Chase JM, Ernest SKM (2017) Community assembly and the functioning of ecosystems: how metacommunity processes alter ecosystems attributes. Ecology, 98, 909–919. https://doi.org/10.1002/ecy.1697
Ovaskainen O, Abrego N (2020) Joint Species Distribution Modelling: With Applications in R. Cambridge University Press, Cambridge. https://doi.org/10.1017/9781108591720
Ovaskainen O, Roy DB, Fox R, Anderson BJ (2016) Uncovering hidden spatial structure in species communities with spatially explicit joint species distribution models. Methods in Ecology and Evolution, 7, 428–436. https://doi.org/10.1111/2041-210X.12502
Pollock LJ, Tingley R, Morris WK, Golding N, O’Hara RB, Parris KM, Vesk PA, McCarthy MA (2014) Understanding co-occurrence by modelling species simultaneously with a Joint Species Distribution Model (JSDM). Methods in Ecology and Evolution, 5, 397–406. https://doi.org/10.1111/2041-210X.12180
Ricklefs RE (2015) Intrinsic dynamics of the regional community. Ecology Letters, 18, 497–503. https://doi.org/10.1111/ele.12431
Soberón J, Nakamura M (2009) Niches and distributional areas: Concepts, methods, and assumptions. Proceedings of the National Academy of Sciences, 106, 19644–19650. https://doi.org/10.1073/pnas.0901637106
Yackulic CB (2017) Competitive exclusion over broad spatial extents is a slow process: evidence and implications for species distribution modeling. Ecography, 40, 305–313. https://doi.org/10.1111/ecog.02836
This is really a superb piece of science; I join both reviewers in their congratulations for one of the best JSDM studies I've seen so far. Really good data, good design, and top-notch interpretation of the results. We all three liked it a lot.
All that said, both reviewers have a number of concerns about the current version of your work. Most of them are minor comments and/or related to clarity, but Carsten Dormann raises a good point about the potential bias in the selection of model parameters caused by the focus on the best models that you may want to consider. As he admits, this is mainly a phylosophical question, but if your study is going to set up a higher standard in JSM, a better account of the uncertainty in model selection would help building a stronger discipline.
Also, both reviewers point to the possibility that there may be some spatial and temporal structure in your data that could be due to unaccounted-for processes or factors. I believe that an additional assessment of whether there is some structure remaining in the residuals of the association matrix could be informative about the existence or not of other processes. Indeed, Joaquín Calatayud highlights that these structured effects may be negligible once you account for environment and co-existence processes, but you need to at least discuss that - and having supplementary analyses to support such discussion would round up your work.
In any case, this is really a great work, and I am looking forward for the resubmission of a new version that accounts and/or discuss the concerns raised by both reviewers. If their comments are properly addressed or discussed it is most likely that I can recommend your preprint in PCI Ecology.