In most vertebrate species, the period of life spanning from departure from the growing site until reaching a more advanced life stage (immature or adult) is critical. During this period, juveniles are often highly vulnerable because they have not reached the morphological, physiological and behavioural maturity levels of adults yet and are therefore at high risk of mortality, e.g. through starvation, depredation or competition (e.g. Marchetti & Price 1989, Wunderle 1991, Naef-Daenzer & Grüebler 2016). In line with this, juvenile survival is most often far lower than adult survival (e.g. Wooller et al. 1992). In species with parental care, juveniles have to acquire behavioural independence from their parents and possibly establish their own territory during this period of life. Very often, this is also the period that is least well-known in the life cycle (Cox et al. 2014, Naef-Daenzer & Grüebler 2016) because of reduced accessibility to individuals and/or adoption of low conspicuous behaviours. Therefore, our understanding of how juveniles acquire typical adult behaviours and how this progressively increases their survival prospects is still very limited (Naef-Daenzer & Grüebler 2016), and questions such as the length of this transition period or the cognitive (e.g. learning, memorization) mechanisms involved remain largely unresolved. This is particularly true regarding the acquisition of independent foraging behaviour (Marchetti & Price 1989).

Because direct observations of juvenile behaviours are usually very difficult except in specific situations or at the cost of an enormous effort, the use of remote tracking devices can be particularly appealing in this context (e.g. Ponchon et al. 2013, Kays et al. 2015). Over the past decades, technical advances have allowed the monitoring of not only individuals’ movements at both large and small spatial scales but also their activities and behaviours based on different parameters recording e.g. speed of movement or diving depth (Whitford & Klimley 2019). Device miniaturization has in particular allowed smaller species to be equipped and/or longer periods of time to be monitored (e.g. Naef-Daenzer et al. 2005). This has opened up whole fields of research, and has been particularly used on marine seabirds. In these species, individuals are most often inaccessible when at sea, representing most of the time outside (and even within) the breeding season, and the life cycle of these long-lived species can include an extended immature period (up to many years) during which most of them will remain unseen, until they come back as breeders or pre-breeders (e.g. Wooller et al. 1992, Oro & Martínez-Abraín 2009). Survival has been found to increase gradually with age in these species before reaching high values characteristic of the adult stage. However, the mechanisms underlying this increase are still to be deciphered.

The study by Delord et al. (2023) builds upon the hypothesis that juveniles gradually learn foraging techniques and movement strategies, improving their foraging efficiency, as previous data on flight parameters seemed to show in different long-lived bird species. Yet, these previous studies obtained data over a limited period of time, i.e. a few months at best. Whether these data could capture the whole dynamics of the progressive acquisition of foraging and movement skills can only be assessed by measuring behaviour over a longer time period and comparing it to similar data in adults, to account for seasonal variation in relation to both resource availability and energetic demands, e.g. due to molt.

The present study (Delord et al. 2023) addresses these questions by taking advantage of longer-lasting recordings of the location and activity of juvenile, immature and adult birds obtained simultaneously to investigate changes over time in juvenile behaviour and thereby provide hints about how young progressively acquire foraging skills. This study is performed on Amsterdam albatrosses, a highly endangered long-lived sea bird, with obvious conservation issues (Thiebot et al. 2015). The results show progressive changes in foraging effort over the first two months after departure from the birth colony, but large differences remain between life stages over a much longer time frame. They also reveal strong variations between sexes and over time in the year. Overall, this study, therefore, confirms the need for very long-term data to be collected in order to address the question of progressive behavioural maturation and associated survival consequences in such species with strongly deferred maturity. Ideally, the same individuals should be monitored over different life stages, from the juvenile period up to adulthood, but this would require further technical development to release the issue of powering duration limitation.

As reviewers emphasized in the first review round, one main challenge now remains to ascertain the outcome of the observed behavioural changes in foraging behaviour: we expect them to reflect improvement in foraging skills and thus performance of juveniles over time, but this would need to be tested. Collecting data on foraging efficiency is yet another challenge, that future technical developments may also help overcome. Importantly also, data were available only for individuals that could be caught again because the tracking device had to be retrieved from the bird. Here, a substantial fraction of the loggers (one-fifth) could not be found again (Delord et al. 2023). To what extent the birds for which no data could be obtained are a random sample of the equipped birds would also need to be assessed. The further development of remote tracking techniques allowing data to be downloaded from a long distance should help further exploration of behavioural ontogeny of juveniles while maturing and its survival consequences. Because the maturation process explored here is likely to show very different characteristics (e.g. timing and speed) in smaller / shorter-lived species (see Cox et al. 2014, Naef-Daenzer & Grüebler 2016), the development of miniaturization is also expected to allow further investigation of post-fledging behavioural maturation in a wider range of bird species. Our understanding of this crucial life phase in different types of species should thus continue to progress in the coming years.

References

Cox W. A., Thompson F. R. III, Cox A. S. & Faaborg J. 2014. Post-fledging survival in passerine birds and the value of post-fledging studies to conservation. Journal of Wildlife Management, 78: 183-193. https://doi.org/10.1002/jwmg.670

Delord K., Weimerskirch H. & Barbraud C. 2023. The challenges of independence: ontogeny of at-sea behaviour in a long-lived seabird. bioRxiv, ver. 6 peer-reviewed and recommended by Peer Community in Ecology. https://doi.org/10.1101/2021.10.23.465439

Kays R., Crofoot M. C., Jetz W. & Wikelski M. 2015. Terrestrial animal tracking as an eye on life and planet. Science, 348 (6240). https://doi.org/10.1126/science.aaa2478

Marchetti K: & Price T. 1989. Differences in the foraging of juvenile and adult birds: the importance of developmental constraints. Biological Reviews, 64: 51-70. https://doi.org/10.1111/j.1469-185X.1989.tb00638.x

Naef-Daenzer B., Fruh D., Stalder M., Wetli P. & Weise E. 2005. Miniaturization (0.2 g) and evaluation of attachment techniques of telemetry transmitters. The Journal of Experimental Biology, 208: 4063–4068. https://doi.org/10.1242/jeb.01870

Naef-Daenzer B. & Grüebler M. U. 2016. Post-fledging survival of altricial birds: ecological determinants and adaptation. Journal of Field Ornithology, 87: 227-250. https://doi.org/10.1111/jofo.12157

Oro D. & Martínez-Abraín A. 2009. Ecology and behavior of seabirds. Marine Ecology, pp.364-389.

Ponchon A., Grémillet D., Doligez B., Chambert T., Tveera T., Gonzàles-Solìs J & Boulinier T. 2013. Tracking prospecting movements involved in breeding habitat selection: insights, pitfalls and perspectives. Methods in Ecology and Evolution, 4: 143-150. https://doi.org/10.1111/j.2041-210x.2012.00259.x

Thiebot J.-B., Delord K., Barbraud C., Marteau C. & Weimerskirch H. 2015. 167 individuals versus millions of hooks: bycatch mitigation in longline fisheries underlies conservation of Amsterdam albatrosses. Aquatic Conservation 26: 674-688. https://doi.org/10.1002/aqc.2578

Whitford M & Klimley A. P. An overview of behavioral, physiological, and environmental sensors used in animal biotelemetry and biologging studies. Animal Biotelemetry, 7: 26. https://doi.org/10.1186/s40317-019-0189-z

Wooller R.D., Bradley J. S. & Croxall J. P. 1992. Long-term population studies of seabirds. Trends in Ecology and Evolution, 7: 111-114. https://doi.org/10.1016/0169-5347(92)90143-y

Wunderle J. M. 1991. Age-specific foraging proficiency in birds. Current Ornithology, 8: 273-324.

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

Humans are dramatically modifying many aspects of our planet via increasing concentrations of carbon dioxide in the atmosphere, patterns of land-use change, and unsustainable exploitation of the planet’s resources. These changes impact the abundance of species of wild organisms, with winners and losers. Identifying how different species and groups of species are influenced by anthropogenic activity in different biomes, continents, and habitats, has become a pressing scientific question with many publications reporting analyses of disparate data on species population sizes. Many conclusions are based on the linear analysis of rather short time series of organismal abundances.
 
There has been particular interest in how arthropods are impacted by environmental change, with several recent papers reporting contradictory results. To investigate why these contradictions might arise, Duchenne et al. (2022) conducted an analysis of four published data sets along with a series of experimental analyses of simulated time series to examine the power of widely used statistical analyses to gain inference on temporal trends. Their important paper reveals that accurate inference on dynamics, particularly of species that exhibit large temporal fluctuations in abundance, requires time series that are substantially longer than are typically collected, as well as careful thought as to whether linear models are appropriate. Linear analyses of short time series are susceptible to providing unreliable inference as trends can be strongly influenced by points at either end of the time series. 
 
Duchenne et al.’s paper provides important insight on the conditions when strong inference on temporal trends of arthropod (and other species) abundances can be made, and when they should be treated with caution. They do not doubt that many insect and arachnid species are changing their abundances, and that patterns in these changes may vary spatially. What their results do say is that we should treat grand claims of population recovery or rapid declines apparently to extinction with caution when they are based on short time series, particularly of species that show significant boom and bust dynamics. In many ways, these results are not unexpected, but it is nice to see such careful and thoughtful analyses and interpretation. More data are required for most arthropod species before clear assessments of abundance trends can be made. Given our reliance on many arthropods for food, pollination, and numerous ecosystem services, and the ability of other species to spread devastating human diseases such as dengue and malaria, it is advisable that we slow our modification of their habitats while additional data are collected to allow us to better characterise the trajectory of arthropod populations to understand what the consequences of our actions on the natural world are likely to be.  
 
References

Duchenne F, Porcher E, Mihoub J-B, Loïs G, Fontaine C (2022) Controversy over the decline of arthropods: a matter of temporal baseline? bioRxiv, 2022.02.09.479422, ver. 3 peer-reviewed and recommended by Peer Community in Ecology. https://doi.org/10.1101/2022.02.09.479422

"Good and evil are so close as to be chained together in the soul"
Robert Louis Stevenson, Dr. Jekyll and Mr. Hyde

Maternally inherited symbionts—microorganisms that pass from a female host to her progeny—have two main ways of increasing their own fitness. First, they can increase the fecundity or viability of infected females. This “positive fitness effects” strategy is the one commonly used by mutualistic symbionts, such as Buchnera aphidicola—the bacterial endosymbiont of the pea aphid, Acyrthosiphon pisum [4]. Second, maternally inherited symbionts can manipulate the reproduction of infected females in a way that enhances symbiont transmission at the expense of host fitness. A famous example of this “reproductive parasitism” strategy is the cytoplasmic incompatibility (CI) [3] induced by bacteria of the genus Wolbachia in their arthropod and nematode hosts. CI works as a toxin-antidote system, whereby the sperm of infected males is modified in a lethal way (toxin) that can only be reverted if the egg is also infected (antidote) [1]. As a result, CI imposes a kind of conditional sterility on their hosts: while infected females are compatible with both infected and uninfected males, uninfected females experience high offspring mortality if (and only if) they mate with infected males [7].

These two symbiont strategies (positive fitness effects versus reproductive parasitism) have been traditionally studied separately, both empirically and theoretically. However, it has become clear that the two strategies are not mutually exclusive, and that a reproductive parasite can simultaneously act as a mutualist—an infection type that has been dubbed “Jekyll and Hyde” [6], after the famous novella by Robert Louis Stevenson about kind scientist Dr. Jekyll and his evil alter ego, Mr. Hyde. In important previous work, Zug and Hammerstein [7] analyzed the consequences of positive fitness effects on the dynamics of different kind of infections, including “Jekyll and Hyde” infections characterized by CI and other reproductive parasitism strategies. Building on this and related modeling framework, Karisto et al. [2] re-investigate and expand on the interplay between positive fitness effects and reproductive parasitism in Wolbachia infections by focusing on CI in both diplodiploid and haplodiploid populations, and by paying particular attention to the mathematical assumption structure underlying their results.

Karisto et al. begin by reviewing classic models of Wolbachia infections in diplodiploid populations that assume a “negative fitness effect” (modeled as a fertility penalty on infected females), characteristic of a pure strategy of reproductive parasitism. Together with the positive frequency-dependent effects due to CI (whereby the fitness benefits to symbionts infecting females increase with the proportion of infected males in the population) this results in population dynamics characterized by two stable equilibria (the Wolbachia-free state and an interior equilibrium with a high frequency of Wolbachia-carrying hosts) separated by an unstable interior equilibrium. Wolbachia can then spread once the initial frequency is above a threshold or an invasion barrier, but is prevented from fixing by a proportion of infections failing to be passed on to offspring. Karisto et al. show that, given the assumption of negative fitness effects, the stable interior equilibrium can never feature a Wolbachia prevalence below one-half. Moreover, they convincingly argue that a prevalence greater than but close to one-half is difficult to maintain in the presence of stochastic fluctuations, as in these cases the high-prevalence stable equilibrium would be too close to the unstable equilibrium signposting the invasion barrier.

Karisto et al. then relax the assumption of negative fitness effects and allow for positive fitness effects (modeled as a fertility premium on infected females) in a diplodiploid population. They show that positive fitness effects may result in situations where the original invasion threshold is now absent, the bistable coexistence dynamics are transformed into purely co-existence dynamics, and Wolbachia symbionts can now invade when rare. Karisto et al. conclude that positive fitness effects provide a plausible and potentially testable explanation for the low frequencies of symbiont-carrying hosts that are sometimes observed in nature, which are difficult to reconcile with the assumption of negative fitness effects. 

Finally, Karisto et al. extend their analysis to haplodiploid host populations (where all fertilized eggs develop as females). Here, they investigate two types of cytoplasmic incompatibility: a female-killing effect, similar to the CI effect studied in diplodiploid populations (the “Leptopilina type” of Vavre et al. [5]) and a masculinization effect, where CI leads to the loss of paternal chromosomes and to the development of the offspring as a male (the “Nasonia type” of Vavre et al. [5]). The models are now two-sex, which precludes a complete analytical treatment, in particular regarding the stability of fixed points. Karisto et al. compensate by conducting large numerical analyses that support their claims. Importantly, all main conclusions regarding the interplay between positive fitness effects and reproductive parasitism continue to hold under haplodiploidy. 

All in all, the analysis and results by Karisto et al. suggest that it is not necessary to resort to classical (but depending on the situation, unlikely) mechanisms, such as ongoing invasion or source-sink dynamics, to explain arthropod populations featuring low-prevalent Wolbachia infections. Instead, low-frequency equilibria might be simply due to reproductive parasites conferring beneficial fitness effects, or Wolbachia symbionts playing Dr. Jekyll (positive fitness effects) and Mr. Hyde (cytoplasmatic incompatibility). 

References

[1] Beckmann JF, Bonneau M, Chen H, Hochstrasser M, Poinsot D, Merçot H, Weill M, Sicard M, Charlat S (2019) The Toxin–Antidote Model of Cytoplasmic Incompatibility: Genetics and Evolutionary Implications. Trends in Genetics, 35, 175–185. https://doi.org/10.1016/j.tig.2018.12.004

[2] Karisto P, Duplouy A, Vries C de, Kokko H (2022) Positive fitness effects help explain the broad range of Wolbachia prevalences in natural populations. bioRxiv, 2022.04.11.487824, ver. 5 peer-reviewed and recommended by Peer Community in Ecology. https://doi.org/10.1101/2022.04.11.487824

[3] Laven H (1956) Cytoplasmic Inheritance in Culex. Nature, 177, 141–142. https://doi.org/10.1038/177141a0

[4] Perreau J, Zhang B, Maeda GP, Kirkpatrick M, Moran NA (2021) Strong within-host selection in a maternally inherited obligate symbiont: Buchnera and aphids. Proceedings of the National Academy of Sciences, 118, e2102467118. https://doi.org/10.1073/pnas.2102467118

[5] Vavre F, Fleury F, Varaldi J, Fouillet P, Bouletreau M (2000) Evidence for Female Mortality in Wolbachia-Mediated Cytoplasmic Incompatibility in Haplodiploid Insects: Epidemiologic and Evolutionary Consequences. Evolution, 54, 191–200. https://doi.org/10.1111/j.0014-3820.2000.tb00019.x

[6] Zug R, Hammerstein P (2015) Bad guys turned nice? A critical assessment of Wolbachia mutualisms in arthropod hosts. Biological Reviews, 90, 89–111. https://doi.org/10.1111/brv.12098

[7] Zug R, Hammerstein P (2018) Evolution of reproductive parasites with direct fitness benefits. Heredity, 120, 266–281. https://doi.org/10.1038/s41437-017-0022-5

The mean trophic level of catches in world fisheries has generally declined over the 20th century, a phenomenon called "fishing down the food web" (Pauly et al. 1998). Several mechanisms have been proposed to explain this decline including the collapse of, or decline in, higher trophic level stocks leading to the inclusion of lower trophic level stocks in the fishery. Fishing down the food web may lead to a reduction in the resilience, i.e., the capacity to rebound from change, of the fished community, which is concerning given the necessity of resilience in the face of climate change. 

The practice of adaptive harvesting, which involves fishing stocks based on their availability, can also result in a reduction in the average trophic level of a fishery (Branch et al. 2010). Adaptive harvesting, similar to adaptive foraging, can affect the resilience of fisheries. Generally, adaptive foraging acts as a stabilizing force in communities (Valdovinos et al. 2010), however it is not clear how including harvesters as the adaptive foragers will affect the resilience of the system.

Tromeur and Loeuille (2023) analyze the effects of adaptively harvesting a trophic community. Using a system of ordinary differential equations representing a predator-prey model where both species are harvested, the researchers mathematically analyze the impact of increasing fishing effort and adaptive harvesting on the mean trophic level and resilience of the fished community. This is achieved by computing the equilibrium densities and equilibrium allocation of harvest effort.  In addition, the researchers numerically evaluate adaptive harvesting in a tri-trophic system (predator, prey, and resource). The study focuses on the effect of adaptively distributing harvest across trophic levels on the mean trophic level of catches, the propensity for regime shifts to occur, the ability to return to equilibrium after a disturbance, and the speed of this return. 

The results indicate that adaptive harvesting leads to a decline in the mean trophic level of catches, resulting in “fishing down the food web”. Furthermore, the study shows that adaptive harvesting may harm the overall resilience of the system. Similar results were observed numerically in a tri-trophic community.

While adaptive foraging is generally a stabilizing force on communities, the researchers found that adaptive harvesting can destabilize the harvested community. One of the key differences between adaptive foraging models and the model presented here, is that the harvesters do not exhibit population dynamics. This lack of a numerical response by the harvesters to decreasing population sizes of their stocks leads to regime shifts. The realism of a fishery that does not respond numerically to declining stock is debatable, however it is very likely that there will a least be significant delays due to social and economic barriers to leaving the fishery, that will lead to similar results.

This study is not unique in demonstrating the ability of adaptive harvesting to result in “fishing down the food web”. As pointed out by the researchers, the same results have been shown with several different model formulations (e.g., age and size structured models). Similarly, this study is not unique to showing that increasing adaptation speeds decreases the resilience of non-linear predator-prey systems by inducing oscillatory behaviours. Much of this can be explained by the destabilising effect of increasing interaction strengths on food webs (McCann et al. 1998). 

By employing a straightforward model, the researchers were able to demonstrate that adaptive harvesting, a common strategy employed by fishermen, can result in a decline in the average trophic level of catches, regime shifts, and reduced resilience in the fished community. While previous studies have observed some of these effects, the fact that the current study was able to capture them all with a simple model is notable. This modeling approach can offer insight into the role of human behavior on the complex dynamics observed in fisheries worldwide.

References

Branch, T. A., R. Watson, E. A. Fulton, S. Jennings, C. R. McGilliard, G. T. Pablico, D. Ricard, et al. 2010. The trophic fingerprint of marine fisheries. Nature 468:431–435. https://doi.org/10.1038/nature09528

Tromeur, E., and N. Loeuille. 2023. Effects of adaptive harvesting on fishing down processes and resilience changes in predator-prey and tritrophic systems. bioRxiv 290460, ver 5 peer-reviewed and recommended by PCI Ecology. https://doi.org/10.1101/290460

McCann, K., A. Hastings, and G.R. Huxel. 1998. Weak trophic interactions and the balance of nature. Nature 395: 794-798. https://doi.org/10.1038/27427

Pauly, D., V. Christensen, J. Dalsgaard, R. Froese, and F. Torres Jr. 1998. Fishing down marine food webs. Science 279:860–86. https://doi.org/10.1126/science.279.5352.860

Valdovinos, F.S., R. Ramos-Jiliberto, L. Garay-Naravez, P. Urbani, and J.A. Dunne. 2010. Consequences of adaptive behaviour for the structure and dynamics of food webs. Ecology Letters 13: 1546-1559. https://doi.org/10.1111/j.1461-0248.2010.01535.x