Ecology is the study of the distribution of organisms in space and time and their interactions. As such, there is a tradition of studies relating abiotic environmental conditions to species distribution, while another one is concerned by the effects of consumers on the abundance of their resources.  Interestingly, joining the dots appears more difficult than it would suggest: eluding the effect of species interactions on distribution remains one of the greatest challenges to elucidate nowadays (Kissling et al. 2012). Theory suggests that yes, species interactions such as predation and competition should influence range limits (Godsoe et al. 2017), but the common intuition among many biogeographers remains that over large areas such as regions and continents, environmental drivers like temperature and precipitation overwhelm their local effects. Answering this question is of primary importance in the context where species are moving around with climate warming.  Inconsistencies in food web structure may arise with asynchronized movements of consumers and their resources, leading to a major disruption in regulation and potentially ecosystem functioning. Solving this problem, however, remains very challenging because we have to rely on observational data since experiments are hard to perform at the biogeographical scale. 

The study of St-Gelais is an interesting step forward to solve this problem. Their main objective was to assess the strength of the association between phytoplankton and zooplankton communities at a large spatial scale, looking at the spatial covariation of both taxonomic and functional composition. To do so, they undertook a massive survey of more than 100 lakes across three regions of the boreal region of Québec. Species and functional composition were recorded, along with a set of abiotic variables. Classic community ecology at this point. The difficulty they faced was to disentangle the multiple causal relationships involved in the distribution of both trophic levels. Teasing apart bottom-up and top-down forces driving the assembly of plankton communities using observational data is not an easy task. On the one hand, both trophic levels could respond to variations in temperature, nutrient availability and dissolved organic carbon. The interpretation is fairly straightforward if the two levels respond to different factors, but the situation is much more complicated when they do respond similarly. There are potentially three possible underlying scenarios. First, the phyto and zooplankton communities may share the same environmental requirements, thereby generating a joint distribution over gradients such as temperature and nutrient availability. Second, the abiotic environment could drive the distribution of the phytoplankton community, which would then propagate up and influence the distribution of the zooplankton community. Alternatively, the abiotic environment could constrain the distribution of the zooplankton, which could then affect the one of phytoplankton. In addition to all of these factors, St-Gelais et al also consider that dispersal may limit the distribution, well aware of previous studies documenting stronger dispersal limitations for zooplankton communities. 

Unfortunately, there is not a single statistical approach that could be taken from the shelf and used to elucidate drivers of co-distribution. Joint species distribution was once envisioned as a major step forward in this direction (Warton et al. 2015), but there are several limits preventing the direct interpretation that co-occurrence is linked to interactions (Blanchet et al. 2020). Rather, St-Gelais used a variety of multivariate statistics to reveal the structure in their observational data. First, using a Procrustes analysis (a method testing if the spatial variation of one community is correlated to the structure of another community), they found a significant correlation between phytoplankton and zooplankton communities, indicating a taxonomic coupling between the groups. Interestingly, this observation was maintained for functional composition only when interaction-related traits were considered. At this point, these results strongly suggest that interactions are involved in the correlation, but it's hard to decipher between bottom-up and top-down perspectives. A complementary analysis performed with a constrained ordination, per trophic level, provided complementary pieces of information. First observation was that only functional variation was found to be related to the different environmental variables, not taxonomic variation. Despite that trophic levels responded to water quality variables, spatial autocorrelation was more important for zooplankton communities and the two layers appear to respond to different variables. 

It is impossible with those results to formulate a strong conclusion about whether grazing influence the co-distribution of phytoplankton and zooplankton communities. That's the mere nature of observational data. While there is a strong spatial association between them, there are also diverging responses to the different environmental variables considered. But the contrast between taxonomic and functional composition is nonetheless informative and it seems that beyond the idiosyncrasies of species composition, trait distribution may be more informative and general. Perhaps the most original contribution of this study is the hierarchical approach to analyze the data, combined with the simultaneous analysis of taxonomic and functional distributions. Having access to a vast catalog of multivariate statistical techniques, a careful selection of analyses helps revealing key features in the data, rejecting some hypotheses and accepting others. Hopefully, we will see more and more of such multi-trophic approaches to distribution because it is now clear that the factors driving distribution are much more complicated than anticipated in more traditional analyses of community data. Biodiversity is more than a species list, it is also all of the interactions between them, influencing their distribution and abundance (Jordano 2016).

References

Blanchet FG, Cazelles K, Gravel D (2020) Co-occurrence is not evidence of ecological interactions. Ecology Letters, 23, 1050–1063. https://doi.org/10.1111/ele.13525

Godsoe W, Jankowski J, Holt RD, Gravel D (2017) Integrating Biogeography with Contemporary Niche Theory. Trends in Ecology & Evolution, 32, 488–499. https://doi.org/10.1016/j.tree.2017.03.008

Jordano P (2016) Chasing Ecological Interactions. PLOS Biology, 14, e1002559. https://doi.org/10.1371/journal.pbio.1002559

Kissling WD, Dormann CF, Groeneveld J, Hickler T, Kühn I, McInerny GJ, Montoya JM, Römermann C, Schiffers K, Schurr FM, Singer A, Svenning J-C, Zimmermann NE, O’Hara RB (2012) Towards novel approaches to modelling biotic interactions in multispecies assemblages at large spatial extents. Journal of Biogeography, 39, 2163–2178. https://doi.org/10.1111/j.1365-2699.2011.02663.x

St-Gelais NF, Vogt RJ, Giorgio PA del, Beisner BE (2021) The taxonomic and functional biogeographies of phytoplankton and zooplankton communities across boreal lakes. bioRxiv, 373332, ver. 4 peer-reviewed and recommended by Peer community in Ecology. https://doi.org/10.1101/373332

Warton DI, Blanchet FG, O’Hara RB, Ovaskainen O, Taskinen S, Walker SC, Hui FKC (2015) So Many Variables: Joint Modeling in Community Ecology. Trends in Ecology & Evolution, 30, 766–779. https://doi.org/10.1016/j.tree.2015.09.007

Wisz MS, Pottier J, Kissling WD, Pellissier L, Lenoir J, Damgaard CF, Dormann CF, Forchhammer MC, Grytnes J-A, Guisan A, Heikkinen RK, Høye TT, Kühn I, Luoto M, Maiorano L, Nilsson M-C, Normand S, Öckinger E, Schmidt NM, Termansen M, Timmermann A, Wardle DA, Aastrup P, Svenning J-C (2013) The role of biotic interactions in shaping distributions and realised assemblages of species: implications for species distribution modelling. Biological Reviews, 88, 15–30. https://doi.org/10.1111/j.1469-185X.2012.00235.x