Philippa Kaur

Balvanera, P. et al. Bioscience 64, 49–57 (2014).Article 

Google Scholar 
Hooper, D. U. et al. Ecol. Monogr. 75, 3–35 (2005).Article 

Google Scholar 
Balvanera, P. et al. Ecol. Lett. 9, 1146–1156 (2006).Article 
PubMed 

Google Scholar 
Cardinale, B. J. et al. Am. J. Bot. 98, 572–592 (2011).Article 
PubMed 

Google Scholar 
Cardinale, B. J. B. J. et al. Nature 486, 59–67 (2012).Article 
CAS 
PubMed 

Google Scholar 
Manning, P. et al. in Advances in Ecological Research (eds Eisenhauer N. et al.) 323–356 (Academic, 2019).Le Provost, G. et al. Nat. Ecol. Evol. https://doi.org/10.1038/s41559-022-01918-5 (2022).Felipe-Lucia, M. R. et al. Proc. Natl Acad. Sci. USA 117, 28140–28149 (2020).Article 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Foley, J. A. et al. Science 309, 570–574 (2005).Article 
CAS 
PubMed 

Google Scholar 
Cardinale, B. J. et al. Ecology 94, 1697–1707 (2013).Article 
PubMed 

Google Scholar 
Teles da Mota, V. & Pickering, C. J. Outdoor Recreat. Tour. 30, 100295 (2020).Article 

Google Scholar 
Mitchell, M. G. E. et al. Trends Ecol. Evol. 30, 190–198 (2015).Article 
PubMed 

Google Scholar 
Raudsepp-Hearne, C. & Peterson, G. D. Ecol. Soc. 21, 16 (2016).Article 

Google Scholar 
Chaplin-Kramer, R. & Kremen, C. Ecol. Appl. 22, 1936–1948 (2012).Article 
PubMed 

Google Scholar  More

Chauhan, B. S. Grand challenges in weed management. Front. Agron. https://doi.org/10.3389/fagro.2019.00003 (2020).Article 

Google Scholar 
Oerke, E. C. Crop losses to pests. J. Agric. Sci. 144, 31–43 (2006).
Google Scholar 
Samejima, H. & Sugimoto, Y. Recent research progress in combatting root parasitic weeds. Biotechnol. Biotechnol. Equip. 32(2), 221–240 (2018).CAS 

Google Scholar 
Aly, R. Conventional and biotechnological approaches for control of parasitic weeds. In Vitro Cell. Dev. Biol. Plant 43(4), 304–317 (2007).
Google Scholar 
Fernández-Aparicio, M., Delavault, P. & Timko, M. P. Management of infection by parasitic weeds: A review. Plants 9(9), 1184 (2020).PubMed Central 

Google Scholar 
Rodenburg, J., Demont, M., Zwart, S. J. & Bastiaans, L. Parasitic weed incidence and related economic losses in rice in Africa. Agric. Ecosyst. Environ. 235, 306–317 (2016).
Google Scholar 
Weisberger, D., Nichols, V. & Liebman, M. Does diversifying crop rotations suppress weeds? A meta-analysis. PLoS One 14(7), e0219847 (2019).CAS 
PubMed 
PubMed Central 

Google Scholar 
Ejeta, G. The Striga scourge in Africa: A growing pandemic. In Integrating New Technologies for Striga Control: Towards Ending the Witch-hunt 3–16 (World Scientific, 2007). https://doi.org/10.1142/9789812771506_0001.Chapter 

Google Scholar 
Netting, R. M. & Stone, M. P. Agro-diversity on a farming frontier: Kofyar smallholders on the Benue plains of central Nigeria. Africa 66(1), 52–70 (1996).
Google Scholar 
Pimentel, D. et al. Conserving biological diversity in agricultural and forestry systems. Bioscience 42, 354–362 (1992).
Google Scholar 
Khoshbakht, K. & Hammer, K. How many plant species are cultivated?. Genet. Resour. Crop Evol. 55(7), 925–928. https://doi.org/10.1007/s10722-008 (2008).Article 

Google Scholar 
Hajjar, R., Jarvis, D. I. & Gemmill-Herren, B. The utility of crop genetic diversity in maintaining ecosystem services. Agric. Ecosyst. Environ. 123(4), 261–270 (2008).
Google Scholar 
He, H. M. et al. Crop diversity and pest management in sustainable agriculture. J. Integr. Agric. 18(9), 1945–1952 (2019).
Google Scholar 
Ofori, F. & Stern, W. R. Cereal–legume intercropping systems. Adv. Agron. 41, 41–90 (1987).
Google Scholar 
Tanveer, M., Anjum, S. A., Hussain, S., Cerdà, A. & Ashraf, U. Relay cropping as a sustainable approach: Problems and opportunities for sustainable crop production. Environ. Sci. Pollut. Res. 24(8), 6973–6988 (2017).
Google Scholar 
Hartwig, N. L. & Ammon, H. U. Cover crops and living mulches. Weed Sci. 50(6), 688–699 (2002).CAS 

Google Scholar 
Raseduzzaman, M. D. & Jensen, E. S. Does intercropping enhance yield stability in arable crop production? A meta-analysis. Eur. J. Agron. 91, 25–33 (2017).
Google Scholar 
Davis, A. S., Hill, J. D., Chase, C. A., Johanns, A. M. & Liebman, M. Increasing cropping system diversity balances productivity, profitability and environmental health. PLoS One 7(10), e47149 (2012).ADS 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Himmelstein, J., Ares, A., Gallagher, D. & Myers, J. A meta-analysis of intercropping in Africa: Impacts on crop yield, farmer income, and integrated pest management effects. Int. J. Agric. Sustain. 15(1), 1–10 (2017).
Google Scholar 
Abson, D. J., Fraser, E. D. & Benton, T. G. Landscape diversity and the resilience of agricultural returns: A portfolio analysis of land-use patterns and economic returns from lowland agriculture. Agric. Food Secur. 2(1), 1–15 (2013).
Google Scholar 
Renard, D. & Tilman, D. National food production stabilized by crop diversity. Nature 571(7764), 257–260 (2019).ADS 
CAS 
PubMed 

Google Scholar 
Gaudin, A. C. et al. Increasing crop diversity mitigates weather variations and improves yield stability. PLoS One 10(2), e0113261 (2015).PubMed 
PubMed Central 

Google Scholar 
Bowles, T. M. et al. Long-term evidence shows that crop-rotation diversification increases agricultural resilience to adverse growing conditions in North America. One Earth 2(3), 284–293 (2020).ADS 

Google Scholar 
Chauhan, B. S., Singh, R. G. & Mahajan, G. Ecology and management of weeds under conservation agriculture: A review. Crop Prot. 38, 57–65 (2012).
Google Scholar 
Nichols, V., Verhulst, N., Cox, R. & Govaerts, B. Weed dynamics and conservation agriculture principles: A review. Field Crop Res. 183, 56–68 (2015).
Google Scholar 
Banik, P., Midya, A., Sarkar, B. K. & Ghose, S. S. Wheat and chickpea intercropping systems in an additive series experiment: Advantages and weed smothering. Eur. J. Agron. 24(4), 325–332 (2006).
Google Scholar 
Workayehu, T. & Wortmann, C. S. Maize–bean intercrop weed suppression and profitability in Southern Ethiopia. Agron. J. 103(4), 1058–1063 (2011).
Google Scholar 
Haugaard-Nielsen, H., Ambus, P. & Jensen, E. S. Interspecific competition, N use and interference with weeds in pea barley intercropping. Field Crop Res. 70, 101–109 (2001).
Google Scholar 
Jensen, E. S. Intercropping of Cereals and Grain Legumes for Increased Production, Weed Control, Improved Product Quality and Prevention of N-losses in European Organic Farming Systems, Final Report on Intercrop Project (QLK5-CT-2002-02352) (Risø National Laboratory, 2006).Arlauskienė, A., Šarūnaitė, L., Kadžiulienė, Ž, Deveikytė, I. & Maikštėnienė, S. Suppression of annual weeds in pea and cereal intercrops. Agron. J. 106(5), 1765–1774 (2014).
Google Scholar 
Szumigalski, A. & van Acker, R. Weed suppression and crop production in annual intercrops. Weed Sci. 53(6), 813–825 (2005).CAS 

Google Scholar 
Stoltz, E. & Nadeau, E. Effects of intercropping on yield, weed incidence, forage quality and soil residual N in organically grown forage maize (Zea mays L.) and faba bean (Vicia faba L.). Field Crop Res. 169, 21–29 (2014).
Google Scholar 
Sauerborn, J., Müller-Stöver, D. & Hershenhorn, J. The role of biological control in managing parasitic weeds. Crop Prot. 26(3), 246–254 (2007).
Google Scholar 
Jamil, M., Rodenburg, J., Charnikhova, T. & Bouwmeester, H. J. Pre-attachment Striga hermonthica resistance of New Rice for Africa (NERICA) cultivars based on low strigolactone production. New Phytol. 192(4), 964–975. https://doi.org/10.1111/j.1469-8137.2011.03850.x (2011).Article 
CAS 
PubMed 

Google Scholar 
Yoneyama, K. et al. Nitrogen deficiency as well as phosphorus deficiency in sorghum promotes the production and exudation of 5-deoxystrigol, the host recognition signal for arbuscular mycorrhizal fungi and root parasites. Planta 227(1), 125–132. https://doi.org/10.1007/s00425-007-0600-5 (2007).Article 
CAS 
PubMed 

Google Scholar 
Sauerborn, J. Legumes used for weed control in agroecosystems in the tropics. Plant Res. Dev. 50, 74–82 (1999).
Google Scholar 
Ejeta, G. & Butler, L. G. Host-parasite interactions throughout the Striga life cycle, and their contributions to Striga resistance. Afr. Crop Sci. J. 1(2), 75–80. https://doi.org/10.4314/acsj.v1i2.69889 (1993).Article 

Google Scholar 
Carsky, R. J., Singh, L. & Ndikawa, R. Suppression of Striga hermonthica on sorghum using a cowpea intercrop. Exp. Agric. 30(3), 349–358. https://doi.org/10.1017/s0014479700024467 (1994).Article 

Google Scholar 
Hsiao, A. I., Worsham, A. D. & Moreland, D. E. Effects of temperature and dl-strigol on seed conditioning and germination of witchweed (Striga asiatica). Ann. Bot. 61(1), 65–72. https://doi.org/10.1093/oxfordjournals.aob.a087528 (1988).Article 
CAS 

Google Scholar 
Patterson, D. T. Effects of Environment on Growth and Reproduction of Witchweed and the Ecological Range of Witchweed (Monograph Series of the Weed Science Society of America, 1990).Stewart, G. R. & Press, M. C. The physiology and biochemistry of parasitic angiosperms. Annu. Rev. Plant Biol. 41(1), 127–151. https://doi.org/10.1146/annurev.pp.41.060190.001015 (1990).Article 
CAS 

Google Scholar 
Anil, L., Park, R. H. P. & Miller, F. A. Temperate intercropping of cereals for forage: A review of the potential for growth and utilization with particular reference to the UK. Grass Forage Sci. 53, 301–317 (1998).
Google Scholar 
Mamolos, A. & Kalburtji, K. Significance of allelopathy in crop rotation. J. Crop Prod. 4, 197–218 (2001).
Google Scholar 
Khan, T. D., Chung, M. I., Xuan, T. D. & Tawata, S. The exploitation of crop allelopathy in sustainable agricultural production. J. Agron. Crop Sci. 191(3), 172–184 (2005).
Google Scholar 
Cissoko, M., Boisnard, A., Rodenburg, J., Press, M. C. & Scholes, J. D. New Rice for Africa (NERICA) cultivars exhibit different levels of post-attachment resistance against the parasitic weeds Striga hermonthica and Striga asiatica. New Phytol. 192(4), 952–963 (2011).CAS 
PubMed 

Google Scholar 
Rodenburg, J. et al. Do NERICA rice cultivars express resistance to Striga hermonthica (Del.) Benth. and Striga asiatica (L.) Kuntze under field conditions?. Field Crop Res. 170, 83–94 (2015).
Google Scholar 
Randrianjafizanaka, M. T., Autfray, P., Andrianaivo, A. P., Ramonta, I. R. & Rodenburg, J. Combined effects of cover crops, mulch, zero-tillage and resistant varieties on Striga asiatica (L.) Kuntze in rice-maize rotation systems. Agric. Ecosyst. Environ. 256, 23–33 (2018).
Google Scholar 
Rodenburg, J. et al. Genetic variation and host–parasite specificity of Striga resistance and tolerance in rice: The need for predictive breeding. New Phytol. 214(3), 1267–1280. https://doi.org/10.1111/nph.14451 (2017).Article 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Nickrent, D. L. & Musselman, L. J. Introduction to parasitic flowering plants. Plant Health Instr. 13(6), 300–315 (2004).
Google Scholar 
Parker, C. Parasitic weeds: A world challenge. Weed Sci. 60(2), 269–276 (2012).CAS 

Google Scholar 
Shai Vaingast 2014. im2graph. Retrieved from: https://www.im2graph.co.il/free-downloads/windows-3264bit/ (2014).Google Maps 2021. https://maps.google.com [Accessed February 2021–December 2022].Kambach, S. et al. Consequences of multiple imputation of missing standard deviations and sample sizes in meta-analysis. Ecol. Evol. 10(20), 11699–11712 (2020).PubMed 
PubMed Central 

Google Scholar 
Nakagawa, S. & Freckleton, R. P. Missing inaction: The dangers of ignoring missing data. Trends Ecol. Evol. 23(11), 592–596 (2008).PubMed 

Google Scholar 
Idris, N. R. N. & Robertson, C. The effects of imputing the missing standard deviations on the standard error of meta analysis estimates. Commun. Stat. Simul. Comput. 38(3), 513–526. https://doi.org/10.1080/03610910802556106 (2009).Article 
MathSciNet 
MATH 

Google Scholar 
van Buuren, S. & Groothuis-Oudshoorn, K. mice: Multivariate imputation by chained equations in R. J. Stat. Softw. 45, 1–67 (2011).
Google Scholar 
van Buuren, S. Flexible Imputation of Missing Data (CRC Press, 2018).MATH 

Google Scholar 
Fick, S. E. & Hijmans, R. J. WorldClim 2: New 1-km spatial resolution climate surfaces for global land areas. Int. J. Climatol. 37, 4302–4315. https://doi.org/10.1002/joc.5086 (2017).Article 

Google Scholar 
O’Donnell, M. S. & Ignizio, D. A. Bioclimatic predictors for supporting ecological applications in the conterminous United States. US Geol. Surv. Data Ser. 691(10), 4–9 (2012).
Google Scholar 
Reuter, H. I., Nelson, A. & Jarvis, A. An evaluation of void filling interpolation methods for SRTM data. Int. J. Geogr. Inf. Sci. 21(9), 983–1008 (2007).
Google Scholar 
CGIAR—Consortium for Spatial Information. http://srtm.csi.cgiar.org © 2004–2021. Accessed September 19, 2021, via: http://srtm.csi.cgiar.org/srtmdata/.QGIS Development Team. QGIS Geographic Information System http://qgis.osgeo.org (Open Source Geospatial Foundation Project, 2020).Kuznetsova, A., Brockhoff, P. B. & Christensen, R. H. B. lmerTest Package: Tests in linear mixed effects models. J. Stat. Softw. 82(13), 26. https://doi.org/10.18637/jss.v082.i13 (2017).Article 

Google Scholar 
Song, C., Peacor, S. D., Osenberg, C. W. & Bence, J. R. An assessment of statistical methods for non-independent data in ecological meta-analyses. Ecology 101(12), e03184. https://doi.org/10.1002/ecy.3184 (2020).Article 
PubMed 

Google Scholar 
Del Rey, A. C. compute.es: Compute Effect Sizes. R package version 0.2-2. https://cran.r-project.org/package=compute.es (2013).R Core Team. R: A language and environment for statistical computing. http://www.R-project.org/ (R Foundation for Statistical Computing, 2020).Wickham, H., Francois, R., Henry, L. & Müller, K. dplyr: A grammar of data manipulation. R package version 0.4. 3 (2015)Bates, D., Mächler, M., Bolker, B. & Walker, S. Fitting linear mixed-effects models using lme4. J. Stat. Softw. 67(1), 48. https://doi.org/10.18637/jss.v067.i01 (2015).Article 

Google Scholar 
Liebman, M. & Dyck, E. Crop rotation and intercropping strategies for weed management. Ecol. Appl. 3(1), 92–122 (1993).PubMed 

Google Scholar 
Pumariño, L. et al. Effects of agroforestry on pest, disease and weed control: A meta-analysis. Basic Appl. Ecol. 16(7), 573–582 (2015).
Google Scholar 
Kuyah, S., Whitney, C. W., Jonsson, M., Sileshi, G. W., Öborn, I., Muthuri, C. W. & Luedeling, E. Agroforestry delivers a win-win solution for ecosystem services in sub-Saharan Africa. A meta-analysis (2019).Evidente, A., Fernandez-Aparicio, M., Andolfi, A., Rubiales, D. & Motta, A. Trigoxazonane, a mono-substituted trioxazonane from Trigonella foenum-graecum root exudates, inhibits Orobanche crenata seed germination. Phytochemistry 68, 2487–2492 (2007).CAS 
PubMed 

Google Scholar 
Khan, Z. R. et al. Control of witchweed Striga hermonthica by intercropping with Desmodium spp., and the mechanism defined as allelopathic. J. Chem. Ecol. 28(9), 1871–1885 (2002).CAS 
PubMed 

Google Scholar 
Nakagawa, S. et al. Methods for testing publication bias in ecological and evolutionary meta-analyses. Methods Ecol. Evol. 13(1), 4–21 (2022).
Google Scholar 
Bakker, A. et al. Beyond small, medium, or large: Points of consideration when interpreting effect sizes. Educ. Stud. Math. 102(1), 1–8 (2019).
Google Scholar 
Scott, D. et al. Mapping the drivers of parasitic weed abundance at a national scale: A new approach applied to Striga asiatica in the mid-west of Madagascar. Weed Res. 60(5), 323–333 (2020).
Google Scholar 
Scott, D. et al. Identifying existing management practices in the control of Striga asiatica within rice–maize systems in mid-west Madagascar. Ecol. Evol. 11(19), 13579–13592 (2021).PubMed 
PubMed Central 

Google Scholar 
Rubiales, D. & Fernández-Aparicio, M. Innovations in parasitic weeds management in legume crops. A review. Agron. Sustain. Dev. 32(2), 433–449 (2012).CAS 

Google Scholar 
Bir, M. S. H. et al. Weed population dynamics under climatic change. Weed Turfgrass Sci. 3(3), 174–182 (2014).
Google Scholar 
Mohamed, K. I., Bolin, J. F., Musselman, L. J. & Townsend, P. A. Genetic diversity of Striga and implications for control and modelling future distributions. In Integrating New Technologies for Striga Control—Towards Ending the Witch-Hunt (eds Ejeta, G. & Gressel, J.) 71–84 (World Scientific, 2007).
Google Scholar 
Mandumbu, R., Mutengwa, C. S., Mabasa, S. & Mwenje, E. Predictions of the Striga scourge under new climate in southern Africa. J. Biol. Sci. 17, 192–201. https://doi.org/10.3923/jbs.2017.194.201 (2017).Article 

Google Scholar 
Mudereri, B. T. et al. Multi-source spatial data-based invasion risk modelling of Striga (Striga asiatica) in Zimbabwe. GIScience Remote Sens. 57(4), 553–571. https://doi.org/10.1080/15481603.2020.1744250 (2020).Article 

Google Scholar  More

Study siteWe conducted laboratory studies at the field site in Finca Irlanda, which is a 300-hectare organic shaded coffee farm located at 1100-m altitude, in the municipality of Tapachula, the state of Chiapas in Southern Mexico (92° 20′ 29″ W and 15° 10′ 65″ N). For the laboratory experiments, all organisms were freshly collected from Finca Irlanda or reared in the lab from insects collected from the field close by. The lab and field work was performed with a permit from the farm owner the Peters family.Ant aggression testTo examine the effect of phorid flies (P. lascinosus) on the aggressivity of ants (A. sericeasur) towards the parasitoids of the beetle larvae (A. orbigera), we conducted an ant aggression test with two treatments: with and without phorids. In the first treatment, a small coffee branch containing two leaves with scale insects (C. viridis) and 20 ant workers were both introduced into a one-liter plastic container. This was done to mimic as much as possible field conditions where the ants are tending scale insects. After waiting for at least 15 min for the ants to calm down and start tending the scale insects, one third- or fourth-instar larva of the beetle was introduced. In the second treatment, all settings were the same except for the addition of 3–4 phorid flies. Once the two treatments were established, one female parasitoid wasp (H. shuvakhinae) was released into each container. During a forty-minute trial, each time that a parasitoid wasp encountered an ant worker, the response of the ant individual was recorded. Ant responses to parasitoids were classified into two categories: (1) the ant ignores the wasp; (2) the ant attacks the wasp. All insects were used for a single replicate and then discarded. A total of four replicates were completed for both the presence and absence of phorids. For each trial, we calculated the proportion of actions (either aggressive or none) by ants when encountering the parasitoid wasp in the treatments with and without phorid flies. We used R36 to conduct a two-sample Mann–Whitney U test on the proportion of ant actions.Parasitism experiments and analysesTo examine the parasitoid wasp’s host preference and the effect of the 1st degree and the 2nd degree HOIs on the beetle’s parasitism and sex ratio, we conducted a laboratory experiment in insect tents (60 cm × 60 cm × 60 cm) with three treatments: (1) no ants (no HOIs but only the wasp and the beetle larvae), (2) ants (1st degree HOI), and (3) ants and phorids (1st and 2nd degree HOIs) (Fig. 1-B). We randomly assigned insect tents to each treatment in each trial, and the tents for each treatment were also shuffled in each trial. All beetle larvae used for these experiments were reared in the laboratory for at least two generations from freshly collected beetle adults. In each tent we placed a coffee branch with 4–6 leaves infested with approximately 100 scale insects inside a plastic container at the center of an insect tent. The set up for the three treatments of species combinations were as follows: (1) 4–5 third or fourth instar beetle larvae and a parasitoid wasp; (2) 4–5 third or fourth instar beetle larvae, a parasitoid wasp, and about 60–80 ant workers; (3) 4–5 third or fourth instar beetle larvae, a parasitoid wasp, about 60–80 ant workers and 3–4 phorid flies. Organism densities in these treatments were close to those observed in the field. To allow for acclimation, we introduced organisms into the tents in the following order: first, we introduced the coffee branch containing scales, immediately followed by the ants (in treatments 2 and 3). After the ants settled down and started tending the scale insects, we introduced the beetle larvae. Once the larvae began moving on the coffee leaves, we introduced the phorids (in treatment 3). When the three treatments were established, and the organisms exhibit normal behavior, we released one lab-reared female parasitoid wasp (H. shuvakhinae) in each tent (treatments 1, 2, and 3). We allowed the organisms to interact for 24 h. After 24 h, we collected all beetle larvae in each treatment and reared them with sufficient scale insects as food, until beetle adults emerged or parasitism symptoms appeared (parasitized larvae turned into hardened black mummies). The treatments of no HOI and 1st + 2nd degree HOI were repeated for 10 consecutive times, and the treatment of 1st degree HOI was repeated for 11 consecutive times, with new individuals of each organism. We recorded parasitism instances and beetle sexes upon emergence. To estimate the sex ratio without parasitoid influence, 78 randomly selected beetle individuals were reared on coffee leaves with scale insects without any interaction with other organisms.To analyze the effect of the parasitoid, the ant and the phorid fly on the parasitism rate and the sex ratio of the beetle, we developed a nested model, starting from$$logitleft(widehat{P}(S)right)=a+bA$$where (widehat{P}(S)) is the probability of an individual being parasitized, A is a binary variable, standing for the absence (0) and presence (1) of ants, a is the baseline probability of parasitism, and b is the magnitude of parasitism altered by ants in the logistic function. We further hypothesized that phorid attacks modify the strength of the interaction modification that ants exert upon the host-parasitoid interaction. Therefore,$$b=g+hP$$where P is another binary variable, standing for the presence (1) and absence (0) of phorids. Substituting b, we obtain the following function,$$logitleft(widehat{P}(S)right)=a+gA+hAP$$where g represents the effect of ants on the parasitism rate of A. orbigera larvae, and h represents the effect of the fly’s facilitation, via interfering with the ant’s interference on the parasitism rate of A. orbigera larvae. We used binary responses (1: survival; 0: parasitized) of all available beetle individuals across the three treatments. We performed model selection based on the Akaike Information Criterion (AIC) and likelihood ratio tests. For the latter, we started model selection by fitting the full model and preceding each step by eliminating the term that had the least significance (the greatest p-value) on the explanation of the dependent variable. The analysis was performed with the application of the bbmle package in R. By doing this, we determined the maximum likelihood estimates of survival probability of the beetle, (widehat{P}(S)), in the three treatments: (1) A = 0, AP = 0 (no HOI); (2) A = 1, AP = 0 (one HOI: ant interference) and (3) A = 1, AP = 1 (interacting HOIs: phorid interference with ant interference), and errors associated with these estimates.The same idea applies to the sex ratio of the beetle under the influence of various organisms. We developed the following equation,$$logitleft(widehat{P}(F|S)right)= r+mA+nAP$$where (widehat{P}(F|S)) is the probability of a parasitism survivor being female. A and P are both binary variables. Respectively, they represent the ant and the phorid fly, and the numeric attributes, 0 and 1, denote their absence and presence. As before, model selection and parameter estimates were conducted with AIC. By doing this, we determined (widehat{P}(F|S)), the estimate of being a female beetle given survival, for the three treatments: (1) A = 0, AP = 0 (no HOI); (2) A = 1, AP = 0 (one HOI: ant interference) and (3) A = 1, AP = 1 (interacting HOIs: phorid interference with ant interference), and errors associated with these estimates. We employed the mle2 function in the bbmle package in R to estimate the female probability (1) in the absence of HOI (the beetle and the parasitoid alone), (2) in the presence of the 1st degree HOI (the beetle, the parasitoid and the ant), and (3) in the presence of the 1st and the 2nd degree HOIs (the beetle, the parasitoid, the ant and the phorid fly).Probabilities of per capita female and per capita male survival from parasitism under the influence of ant and the phorid flyTo test whether the sex ratio of beetle survivors’ population is due to sex-differential survival probability, Bayes’ theorem was employed. Per capita female survival probability from parasitism in each treatment of the parasitism experiment was derived based on (widehat{P}(F)), (widehat{P}left(F|Sright),) and (widehat{P}(S)), and per capita male survival probability was derived based on (widehat{P}(M)), (widehat{P}left(M|Sright),) and (widehat{P}(S)). According to the Central Limit Theorem, the estimates of proportions, (widehat{P}left(S|Fright)) and (widehat{P}left(S|Mright)), are approximately normally distributed,$$widehat{P}left(S|Fright)sim Nleft(widehat{P}left(S|Fright), sqrt{frac{widehat{P}(S|F)times left(1-widehat{P}left(S|Fright)right)}{{n}^{*}}}right)$$$$widehat{P}left(S|Mright)sim Nleft(widehat{P}left(S|Mright), sqrt{frac{widehat{P}(S|M)times left(1-widehat{P}left(S|Mright)right)}{{n}^{*}}}right)$$with means (widehat{P}left(S|Fright)) and (widehat{P}(S|M)), and standard deviations (sqrt{frac{widehat{P}left(S|Fright)times (1-widehat{P}left(S|Fright))}{{n}^{*}}}) and (sqrt{frac{widehat{P}left(S|Mright)times (1-widehat{P}left(S|Mright))}{{n}^{*}}}), where (widehat{P}(S|F)) and (widehat{P}(S|M)), respectively, are the population proportions of females and males. Here we employ n*, the smallest sample size among those of the three variables in the Bayesian formulas for males and females. Since the three variables have different sample sizes, n* guarantees a conservative estimate of standard error, and thus confidence interval, of each derived probability. More

The plots of Igeo, PERI, and PLI of HMs in the topsoil of the tourist area of Sayram Lake (Fig. 5) reveal the degree of HM pollution and eco-risk in this study area on the one hand and, on the other hand, indicate the direction for the relevant agencies to target soil environmental protection and HM pollution prevention and control measures. In this study, the Igeo results showed that Cd was the most highly enriched HM, and Pb, Zn, Cd, and Ni were slightly enriched in a few sample sites. The unnatural accumulation of these elements is usually closely associated with human activities in the area34. Tourism is the main economic activity in the district, and published studies have reported that tourism infrastructure construction (e.g., roads, buildings, etc.) and tourism wastes (e.g., plastic bags, batteries, hotel wastewater) release Cd into the soil35. Additionally, the accumulation of Pb, Zn, Cu and Ni in soils is usually associated with traffic emissions36. The PERI showed that the study area was at low risk overall, with only point ss04 exhibiting medium risk; however, this result was caused by the abnormally high Cd concentration value (Fig. 4) at point ss04 (Cd (concentration): 1.08 mg/kg, Cd (background): 0.34 mg/kg). This anomalous concentration value has a large influence on the PERI calculated based on the measured concentration, the background value and the toxicity coefficient. Therefore, references to this point can be appropriately removed when considering eco-risk. The PLI of each sampling point was greater than 1 and less than 2, which means that the area was in a moderately contaminated state. In general, the degree of soil HM contamination in this area was low; however, due to HM toxicity, bioaccumulation, and persistence37, the HM contamination of this area still requires sustained attention.Figure 5Contamination and ecological risk indices: (a) geoaccumulation index (Igeo) of HMs; (b) ecological risk of individual HMs; (c) potential ecological risk index (PERI) of HMs; (d) pollution load index (PLI) of HMs.Full size imageCorrelation analysis is an efficient way to reveal correlations among HMs through Pearson correlation coefficients, and HMs with significant correlations may originate from the same source38. As shown in Table S5, the elemental pairs Cd-Cu (p  More

Réale, D., Reader, S. M., Sol, D., McDougall, P. T. & Dingemanse, N. J. Integrating animal temperament within ecology and evolution. Biol. Rev. 82, 291–318 (2007).PubMed 

Google Scholar 
Kaiser, M. I. & Müller, C. What is an animal personality?. Biol. Philos. 36, 1 (2021).
Google Scholar 
Sih, A., Bell, A. & Johnson, J. C. Behavioral syndromes: An ecological and evolutionary overview. Trends Ecol. Evol. 19, 372–378 (2004).PubMed 

Google Scholar 
Gosling, S. D. From mice to men: What can we learn about personality from animal research?. Psychol. Bull. 127, 45–86 (2001).PubMed 

Google Scholar 
Biro, P. A. & Stamps, J. A. Are animal personality traits linked to life-history productivity?. Trends Ecol. Evol. 23, 361–368 (2008).PubMed 

Google Scholar 
Koolhaas, J. M. Coping style and immunity in animals: Making sense of individual variation. Brain Behav. Immun. 22, 662–667 (2008).PubMed 

Google Scholar 
Réale, D. et al. Personality and the emergence of the pace-of-life syndrome concept at the population level. Philos. Trans. R. Soc. B Biol. Sci. 365, 4051–4063 (2010).
Google Scholar 
Stamps, J. A. Growth-mortality tradeoffs and ‘personality traits’ in animals. Ecol. Lett. 10, 355–363 (2007).PubMed 

Google Scholar 
Finkemeier, M. A., Langbein, J. & Puppe, B. Personality research in mammalian farm animals: Concepts, measures, and relationship to welfare. Front Vet. Sci. 10(5), 355–363 (2018).
Google Scholar 
Murphy, E., Nordquist, R. E. & van der Staay, F. J. A review of behavioural methods to study emotion and mood in pigs. Sus. Scrofa. Appl. Anim. Behav. Sci 159, 9–28 (2014).
Google Scholar 
Lauber, M. C. Y., Hemsworth, P. H. & Barnett, J. L. The effects of age and experience on behavioural development in dairy calves. Appl. Anim. Behav. Sci. 99, 41–52 (2006).
Google Scholar 
Neave, H. W., Costa, J. H. C., Weary, D. M. & von Keyserlingk, M. A. G. Personality is associated with feeding behavior and performance in dairy calves. J. Dairy Sci. 101, 7437–7449 (2018).PubMed 

Google Scholar 
Foris, B., Zebunke, M., Langbein, J. & Melzer, N. Evaluating the temporal and situational consistency of personality traits in adult dairy cattle. Plos One 13, e0204619 (2018).PubMed 
PubMed Central 

Google Scholar 
Dingemanse, N. J. & Dochtermann, N. A. Quantifying individual variation in behaviour: Mixed-effect modelling approaches. J. Anim. Ecol. 82, 39–54 (2013).PubMed 

Google Scholar 
Dingemanse, N. J., Kazem, A. J. N., Réale, D. & Wright, J. Behavioural reaction norms: Animal personality meets individual plasticity. Trends Ecol. Evol. 25, 81–89 (2010).PubMed 

Google Scholar 
Nakagawa, S. & Schielzeth, H. Repeatability for Gaussian and non-Gaussian data: A practical guide for biologists. Biol. Rev. https://doi.org/10.1111/j.1469-185X.2010.00141.x (2010).Article 
PubMed 

Google Scholar 
Bell, A. M., Hankison, S. J. & Laskowski, K. L. The repeatability of behaviour: A meta-analysis. Anim. Behav. 77, 771–783 (2009).PubMed 
PubMed Central 

Google Scholar 
Neave, H. W., Costa, J. H. C., Benetton, J. B., Weary, D. M. & von Keyserlingk, M. A. G. Individual characteristics in early life relate to variability in weaning age, feeding behavior, and weight gain of dairy calves automatically weaned based on solid feed intake. J. Dairy Sci. 102, 10250–10265 (2019).PubMed 

Google Scholar 
Berckmans, D. Precision livestock farming technologies for welfare management in intensive livestock systems. Rev. Sci. Tech. OIE 33, 189–196 (2014).
Google Scholar 
Carslake, C., Vázquez-Diosdado, J. A. & Kaler, J. Machine learning algorithms to classify and quantify multiple behaviours in dairy calves using a sensor: Moving beyond classification in precision livestock. Sensors 21, 88 (2020).ADS 
PubMed Central 

Google Scholar 
Hertel, A. G., Niemelä, P. T., Dingemanse, N. J. & Mueller, T. A guide for studying among-individual behavioral variation from movement data in the wild. Mov. Ecol. 8(1), 1–18 (2020).
Google Scholar 
Occhiuto, F., Vázquez-Diosdado, J. A., Carslake, C. & Kaler, J. Personality and predictability in farmed calves using movement and space-use behaviours quantified by ultra-wideband sensors. R. Soc. Open Sci. 9, 212019 (2022).ADS 
PubMed 
PubMed Central 

Google Scholar 
Carslake, C., Occhiuto, F., Vázquez-Diosdado, J. A. & Kaler, J. Repeatability and predictability of calf feeding behaviors—quantifying between- and within-individual variation for precision livestock farming. Front. Vet. Sci. https://doi.org/10.3389/fvets.2022.827124 (2022).Article 
PubMed 
PubMed Central 

Google Scholar 
Tolkamp, B. J. & Kyriazakis, I. To split behaviour into bouts, log-transform the intervals. Anim. Behav. 57, 807–817 (1999).PubMed 

Google Scholar 
Houslay, T. M. & Wilson, A. J. Avoiding the misuse of BLUP in behavioural ecology. Behav. Ecol. 28, 948 (2017).PubMed 
PubMed Central 

Google Scholar 
R Core Team. R. Preprint at (2021).Bürkner, P.-C. Advanced bayesian multilevel modeling with the R package brms. R. J. 10, 395 (2018).
Google Scholar 
Dancey, C. P. & Reidy, J. Statistics without maths for psychology (Pearson education, 2007).
Google Scholar 
von Keyserlingk, M. A. G., Brusius, L. & Weary, D. M. Competition for teats and feeding behavior by group-housed dairy calves. J. Dairy Sci. 87, 4190–4194 (2004).
Google Scholar 
Fraley, R. C. & Roberts, B. W. Patterns of continuity: A dynamic model for conceptualizing the stability of individual differences in psychological constructs across the life course. Psychol. Rev. 112, 60–74 (2005).PubMed 

Google Scholar 
Ashcroft, J., Semmler, C., Carnell, S., van Jaarsveld, C. H. M. & Wardle, J. Continuity and stability of eating behaviour traits in children. Eur. J. Clin. Nutr. 62, 985–990 (2008).PubMed 

Google Scholar 
Neave, H. W., Costa, J. H. C., Weary, D. M. & von Keyserlingk, M. A. G. Long-term consistency of personality traits of cattle. R. Soc. Open Sci. 7, 191849 (2020).ADS 
PubMed 
PubMed Central 

Google Scholar 
Müller, R. & von Keyserlingk, M. A. G. Consistency of flight speed and its correlation to productivity and to personality in Bos taurus beef cattle. Appl. Anim. Behav. Sci. 99, 193–204 (2006).
Google Scholar 
Neja, W., Sawa, A., Jankowska, M., Bogucki, M. & Krężel-Czopek, S. Effect of the temperament of dairy cows on lifetime production efficiency. Arch. Anim. Breed 58, 193–197 (2015).
Google Scholar 
Haskell, M. J., Simm, G. & Turner, S. P. Genetic selection for temperament traits in dairy and beef cattle. Front Genet. https://doi.org/10.3389/fgene.2014.00368 (2014).Article 
PubMed 
PubMed Central 

Google Scholar 
Whalin, L., Neave, H. W., Føske Johnsen, J., Mejdell, C. M. & Ellingsen-Dalskau, K. The influence of personality and weaning method on early feeding behavior and growth of Norwegian red calves. J. Dairy Sci. 105, 1369–1386 (2022).PubMed 

Google Scholar 
Dammhahn, M., Dingemanse, N. J., Niemelä, P. T. & Réale, D. Pace-of-life syndromes: A framework for the adaptive integration of behaviour, physiology and life history. Behav. Ecol. Sociobiol. 72(3), 1–8 (2018).
Google Scholar 
Kelly, D. N. et al. Large variability in feeding behavior among crossbred growing cattle. J. Anim. Sci. 98, 1–10 (2020).
Google Scholar 
Neave, H. W., Weary, D. M. & von Keyserlingk, M. A. G. Review: Individual variability in feeding behaviour of domesticated ruminants. Animal 12, S419–S430 (2018).PubMed 

Google Scholar 
DeVries, T. J., von Keyserlingk, M. A. G., Weary, D. M. & Beauchemin, K. A. Measuring the feeding behavior of lactating dairy cows in early to peak lactation. J. Dairy Sci. 86, 3354–3361 (2003).PubMed 

Google Scholar 
Kelly, D. N., Sleator, R. D., Murphy, C. P., Conroy, S. B. & Berry, D. P. Phenotypic and genetic associations between feeding behavior and carcass merit in crossbred growing cattle. J. Anim. Sci. 99, skab285 (2021).PubMed 

Google Scholar 
Weary, D. M., Huzzey, J. M. & von Keyserlingk, M. A. G. Board-invited review: Using behavior to predict and identify ill health in animals. J. Anim. Sci. 87, 770–777 (2009).PubMed 

Google Scholar 
Carter, A. J., Feeney, W. E., Marshall, H. H., Cowlishaw, G. & Heinsohn, R. Animal personality: What are behavioural ecologists measuring?. Biol. Rev. 88, 465–475 (2013).PubMed 

Google Scholar 
Biro, P. A. Do rapid assays predict repeatability in labile (behavioural) traits?. Anim Behav 83, 1295–1300 (2012).
Google Scholar 
Percie du Sert, N. et al. Reporting animal research: Explanation and elaboration for the ARRIVE guidelines 20. Plos Biol. 18, e3000411 (2020).PubMed 
PubMed Central 

Google Scholar  More

Paerl, H. W. Mitigating toxic planktonic cyanobacterial blooms in aquatic ecosystems facing increasing anthropogenic and climatic pressures. Toxins. 10, 1–16 (2018).
Google Scholar 
Harke, M. J. et al. A review of the global ecology, genomics, and biogeography of the toxic cyanobacterium Microcystis spp. Harmful Algae 54, 4–20. https://doi.org/10.1016/j.hal.2015.12.007 (2016).Article 
PubMed 

Google Scholar 
Paerl, H. W. & Barnard, M. A. Mitigating the global expansion of harmful cyanobacterial blooms: Moving targets in a human- and climatically-altered world. Harmful Algae 96, 101845. https://doi.org/10.1016/j.hal.2020.101845 (2020).Article 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Paerl, H. W. Mitigating harmful cyanobacterial blooms in a human- and climatically-impacted world. Life. 4, 988–1012 (2014).PubMed 
PubMed Central 

Google Scholar 
Burford, M. A. et al. Perspective: Advancing the research agenda for improving understanding of cyanobacteria in a future of global change. Harmful Algae 91, 101601. https://doi.org/10.1016/j.hal.2019.04.004 (2020).Article 
CAS 
PubMed 

Google Scholar 
Havens, K. E., James, R. T., East, T. L. & Smith, V. H. N: P ratios, light limitation, and cyanobacterial dominance in a subtropical lake impacted by non-point source nutrient pollution. Environ. Pollut. 122, 379–390 (2003).CAS 
PubMed 

Google Scholar 
Bernard, C. Cyanobacteria and cyanotoxins. Rev. Franç. Lab. 2014, 53–68 (2014).
Google Scholar 
Paerl, H. W. & Otten, T. G. Harmful cyanobacterial blooms: Causes, consequences, and controls. Microb. Ecol. 65, 995–1010 (2013).CAS 
PubMed 

Google Scholar 
Dolman, A. M. et al. Cyanobacteria and cyanotoxins: The influence of nitrogen versus phosphorus. PLoS ONE 7, 38575 (2012).
Google Scholar 
Svirčev, Z. et al. Global geographical and historical overview of cyanotoxin distribution and cyanobacterial poisonings. Arch. Toxicol. https://doi.org/10.1007/s00204-019-02524-4 (2019).Article 
PubMed 

Google Scholar 
Massey, I. Y. & Yang, F. A mini review on microcystins and bacterial degradation. Toxins 12, 268 (2020).CAS 
PubMed Central 

Google Scholar 
Paerl, H. W. et al. Mitigating eutrophication and toxic cyanobacterial blooms in large lakes: The evolution of a dual nutrient (N and P) reduction paradigm. Hydrobiologia 847, 4359–4375. https://doi.org/10.1007/s10750-019-04087-y (2020).Article 
CAS 

Google Scholar 
Sapp, M. et al. Species-specific bacterial communities in the phycosphere of microalgae?. Microb. Ecol. 53, 683–699 (2007).PubMed 

Google Scholar 
Cai, H., Jiang, H., Krumholz, L. R. & Yang, Z. Bacterial community composition of size-fractioned aggregates within the phycosphere of cyanobacterial blooms in a eutrophic freshwater lake. PLoS ONE 9, 102879 (2014).ADS 

Google Scholar 
Grant, M. A. A., Kazamia, E., Cicuta, P. & Smith, A. G. Direct exchange of vitamin B 12 is demonstrated by modelling the growth dynamics of algal-bacterial cocultures. ISME J. Nat. Publ. Group 8, 1418–1427 (2014).CAS 

Google Scholar 
Shi, L., Cai, Y., Kong, F. & Yu, Y. Specific association between bacteria and buoyant Microcystis colonies compared with other bulk bacterial communities in the eutrophic Lake Taihu, China. Environ. Microbiol. Rep. 4, 669–678 (2012).CAS 
PubMed 

Google Scholar 
Brunberg, A. K. Contribution of bacteria in the mucilage of Microcystis spp (Cyanobacteria) to benthic and pelagic bacterial production in a hypereutrophic lake. FEMS Microbiol. Ecol. 29, 13–22 (1999).CAS 

Google Scholar 
Shao, K. et al. The responses of the taxa composition of particle-attached bacterial community to the decomposition of Microcystis blooms. Sci. Total. Environ. 488–489, 236–242. https://doi.org/10.1016/j.scitotenv.2014.04.101 (2014).Article 
ADS 
CAS 
PubMed 

Google Scholar 
Jankowiak, J. G. & Gobler, C. J. The composition and function of microbiomes within microcystis colonies are significantly different than native bacterial assemblages in two North American lakes. Front. Microbiol. 11, 1–26 (2020).
Google Scholar 
Bauer, A. & Forchhammer, K. Bacterial predation on cyanobacteria. Microb. Physiol. 99, 108 (2021).
Google Scholar 
Ndlela, L. L., Oberholster, P. J., Van Wyk, J. H. & Cheng, P. H. Bacteria as biological control agents of freshwater cyanobacteria: Is it feasible beyond the laboratory?. Appl. Microbiol. Biotechnol. 102, 9911–9923 (2018).CAS 
PubMed 

Google Scholar 
Yang, C. et al. Distinct network interactions in particle-associated and free-living bacterial communities during a Microcystis aeruginosa bloom in a plateau lake. Front. Microbiol. 8, 1–15 (2017).
Google Scholar 
Xu, H. et al. Contrasting network features between free-living and particle-attached bacterial communities in Taihu Lake. Microb. Ecol. 76, 303–313 (2018).PubMed 

Google Scholar 
Liu, M. et al. Community dynamics of free-living and particle-attached bacteria following a reservoir Microcystis bloom. Sci. Total Environ. 660, 501–511. https://doi.org/10.1016/j.scitotenv.2018.12.414 (2019).Article 
ADS 
CAS 
PubMed 

Google Scholar 
Parveen, B. et al. Bacterial communities associated with Microcystis colonies differ from free-living communities living in the same ecosystem. Environ. Microbiol. Rep. 5, 716–724 (2013).CAS 
PubMed 

Google Scholar 
Louati, I. et al. Structural diversity of bacterial communities associated with bloom-forming freshwater cyanobacteria differs according to the cyanobacterial genus. PLoS ONE 10, 0140614 (2015).
Google Scholar 
Zwirglmaier, K., Keiz, K., Engel, M., Geist, J. & Raeder, U. Seasonal and spatial patterns of microbial diversity along a trophic gradient in the interconnected lakes of the Osterseen Lake District, Bavaria. Front. Microbiol. 6, 1–18 (2015).
Google Scholar 
Scherer, P. I. et al. Temporal dynamics of the microbial community composition with a focus on toxic cyanobacteria and toxin presence during harmful algal blooms in two South German lakes. Front. Microbiol. 8, 1–17 (2017).
Google Scholar 
Kokocinski, M., Dziga, D., Antosiak, A. & Soininen, J. Are bacterio- and phytoplankton community compositions related in lakes differing in their cyanobacteria contribution and physico-chemical properties?. Genes 12, 855 (2021).CAS 
PubMed 
PubMed Central 

Google Scholar 
Dziga, D. et al. Correlation between specific groups of heterotrophic bacteria and microcystin biodegradation in freshwater bodies of central Europe. FEMS Microbiol. Ecol. https://doi.org/10.1111/j.1574-6941.1999.tb00594.x (2019).Article 
PubMed 

Google Scholar 
Jurczak, T. et al. Elimination of microcystins by water treatment processes: Examples from Sulejow Reservoir, Poland. Water Res. 39, 2394–2406 (2005).CAS 
PubMed 

Google Scholar 
Mankiewicz-Boczek, J. et al. Detection and monitoring toxigenicity of cyanobacteria by application of molecular methods. Environ Toxicol. 21, 380–387 (2006).ADS 
CAS 
PubMed 

Google Scholar 
Rajaniemi-Wacklin, P. et al. Correspondence between phylogeny and morphology of Snowella spp. and Woronichinia naegeliana, cyanobacteria commonly occurring in lakes. J. Phycol. 42, 226–232 (2006).
Google Scholar 
DrobacBacković, D. et al. Cyanobacteria, cyanotoxins, and their histopathological effects on fish tissues in Fehérvárcsurgó reservoir Hungary. Environ. Monit. Assess. https://doi.org/10.1007/s10661-021-09324-3 (2021).Article 

Google Scholar 
Kallscheuer, N. et al. Analysis of bacterial communities in a municipal duck pond during a phytoplankton bloom and isolation of Anatilimnocola aggregata gen. nov., sp. Nov., Lacipirellula limnantheis sp. Nov. and Urbifossiella limnaea gen. nov. sp. nov. belonging to the phylum. Environ. Microbiol. 23, 1379–1396 (2021).CAS 
PubMed 

Google Scholar 
Davis, T. W. et al. Effects of nitrogenous compounds and phosphorus on the growth of toxic and non-toxic strains of Microcystis during cyanobacterial blooms. Aquat. Microb. Ecol. 61, 149–162 (2010).
Google Scholar 
Gobler, C. J., Davis, T. W., Coyne, K. J. & Boyer, G. L. Interactive influences of nutrient loading, zooplankton grazing, and microcystin synthetase gene expression on cyanobacterial bloom dynamics in a eutrophic New York lake. Harmful Algae 6, 119–133 (2007).CAS 

Google Scholar 
Mankiewicz-Boczek, J. et al. Cyanophages infection of microcystis bloom in lowland dam reservoir of Sulejów, Poland. Microb. Ecol. 71, 315–325 (2016).CAS 
PubMed 

Google Scholar 
Davis, T. W., Berry, D. L., Boyer, G. L. & Gobler, C. J. The effects of temperature and nutrients on the growth and dynamics of toxic and non-toxic strains of Microcystis during cyanobacteria blooms. Harmful Algae 8, 715–725 (2009).CAS 

Google Scholar 
Yoshida, M., Yoshida, T., Takashima, Y., Hosoda, N. & Hiroishi, S. Dynamics of microcystin-producing and non-microcystin-producing Microcystis populations is correlated with nitrate concentration in a Japanese lake. FEMS Microbiol. Lett. 266, 49–53 (2007).CAS 
PubMed 

Google Scholar 
Sezenna, M. L. Proteobacteria: Phylogeny, Metabolic Diversity and Ecological Effects (Nova Science Publishers, Inc., 2011).
Google Scholar 
Rilling, J. I., Acuña, J. J., Sadowsky, M. J. & Jorquera, M. A. Putative nitrogen-fixing bacteria associated with the rhizosphere and root endosphere of wheat plants grown in an andisol from southern Chile. Front. Microbiol. 9, 1–13 (2018).
Google Scholar 
Lukumbuzya, M. et al. A refined set of rRNA-targeted oligonucleotide probes for in situ detection and quantification of ammonia-oxidizing bacteria. Water Res. 186, 116375 (2020).
Google Scholar 
Prosser, J. I., Head, I. M. & Stein, L. Y. The family Nitrosomonadaceae. In The Prokaryotes: Alphaproteobacteria and Betaproteobacteria (eds Rosenberg, E. et al.) 901–918 (Springer, 2014). https://doi.org/10.1007/978-3-642-30197-1_372.Chapter 

Google Scholar 
Jia, L., Jiang, B., Huang, F. & Hu, X. Nitrogen removal mechanism and microbial community changes of bioaugmentation subsurface wastewater infiltration system. Bioresour. Technol. 294, 122140. https://doi.org/10.1016/j.biortech.2019.122140 (2019).Article 
CAS 
PubMed 

Google Scholar 
Daft, M. J. & Stewart, W. D. P. Bacterial pathogens of freshwater blue-green algae. New Phytol. 70, 819–829 (1971).
Google Scholar 
Chun, S. J. et al. Network analysis reveals succession of Microcystis genotypes accompanying distinctive microbial modules with recurrent patterns. Water Res. 170, 115326. https://doi.org/10.1016/j.watres.2019.115326 (2020).Article 
CAS 
PubMed 

Google Scholar 
Parulekar, N. N. et al. Characterization of bacterial community associated with phytoplankton bloom in a eutrophic lake in South Norway using 16S rRNA gene amplicon sequence analysis. PLoS ONE 12, 1–22 (2017).
Google Scholar 
Guedes, I. A. et al. Close link between harmful cyanobacterial dominance and associated bacterioplankton in a tropical eutrophic reservoir. Front. Microbiol. 9, 424 (2018).PubMed 
PubMed Central 

Google Scholar 
Allgaier, M. & Grossart, H. P. Seasonal dynamics and phylogenetic diversity of free-living and particle-associated bacterial communities in four lakes in northeastern Germany. Aquat. Microb. Ecol. 45, 115–128 (2006).
Google Scholar 
Chen, S. et al. Disentangling the drivers of Microcystis decomposition: Metabolic profile and co-occurrence of bacterial community. Sci. Total Environ. 739, 140062. https://doi.org/10.1016/j.scitotenv.2020.140062 (2020).Article 
ADS 
CAS 
PubMed 

Google Scholar 
Leflaive, J. & Ten-Hage, L. Algal and cyanobacterial secondary metabolites in freshwaters: A comparison of allelopathic compounds and toxins. Freshw. Biol. 52, 199–214 (2007).CAS 

Google Scholar 
Song, H. et al. Biological and chemical factors driving the temporal distribution of cyanobacteria and heterotrophic bacteria in a eutrophic lake (West Lake, China). Appl. Microbiol. Biotechnol. 101, 1685–1696. https://doi.org/10.1007/s00253-016-7968-8 (2017).Article 
CAS 
PubMed 

Google Scholar 
Bagatini, I. L. et al. Host-specificity and dynamics in bacterial communities associated with bloom-forming freshwater phytoplankton. PLoS ONE 9, 85957 (2014).ADS 

Google Scholar 
Kohler, E. et al. Biodegradation of microcystins during gravity-driven membrane (GDM) ultrafiltration. PLoS ONE 9, 111794 (2014).ADS 

Google Scholar 
Wu, X. et al. Culturing of “unculturable” subsurface microbes: Natural organic carbon source fuels the growth of diverse and distinct bacteria from groundwater. Front. Microbiol. 11, 1–10 (2020).CAS 

Google Scholar 
Morotomi, M., Nagai, F. & Watanabe, Y. Parasutterella secunda sp. no., isolated from human faeces and proposal of Sutterellaceae fam. nov. in the order Burkholderiales. Int. J. Syst. Evol. Microbiol. 61, 637–643 (2011).CAS 
PubMed 

Google Scholar 
Kiedrzyńska, E. et al. Point sources of nutrient pollution in the lowland river catchment in the context of the baltic Sea eutrophication. Ecol. Eng. 70, 337–348 (2014).
Google Scholar 
Hwang, W. M., Ko, Y., Kim, J. H. & Kang, K. Ahniella affigens gen Nov, sp. nov., a gammaproteobacterium isolated from sandy soil near a stream. Int. J. Syst. Evol. Microbiol. 68, 2478–2484 (2018).CAS 
PubMed 

Google Scholar 
Qian, H. et al. Spatial variability of cyanobacteria and heterotrophic bacteria in Lake Taihu (China). Bull. Environ. Contam. Toxicol. 99, 380–384 (2017).CAS 
PubMed 

Google Scholar 
Humbert, J. F. et al. Comparison of the structure and composition of bacterial communities from temperate and tropical freshwater ecosystems. Environ. Microbiol. 11, 2339–2350 (2009).CAS 
PubMed 

Google Scholar 
Newton, R. J., Jones, S. E., Eiler, A., McMahon, K. D. & Bertilsson, S. A guide to the natural history of freshwater lake Bacteria. Microbiol. Mol. Biol. Rev. 1, 1–10 (2011).
Google Scholar 
Parveen, B., Mary, I., Vellet, A., Ravet, V. & Debroas, D. Temporal dynamics and phylogenetic diversity of free-living and particle-associated Verrucomicrobia communities in relation to environmental variables in a mesotrophic lake. FEMS Microbiol. Ecol. 83, 189–201 (2013).CAS 
PubMed 

Google Scholar 
Henson, M. W., Lanclos, V. C., Faircloth, B. C. & Thrash, J. C. Cultivation and genomics of the first freshwater SAR11 (LD12) isolate. ISME J. 12, 1846–1860 (2018).CAS 
PubMed 
PubMed Central 

Google Scholar 
Yang, C. et al. The characteristics and algicidal mechanisms of cyanobactericidal bacteria, a review. World J. Microbiol. Biotechnol. 36, 1–10. https://doi.org/10.1007/s11274-020-02965-5 (2020).Article 

Google Scholar 
Izydorczyk, K. et al. Influence of abiotic and biotic factors on microcystin content in Microcystis aeruginosa cells in a eutrophic temperate reservoir. J. Plankton Res. 30, 393–400 (2008).CAS 

Google Scholar 
Mankiewicz-Boczek, J. et al. Bacteria homologus to Aeromonas capable of microcystin degradation. Open Life Sci. 10, 106–116 (2015).CAS 

Google Scholar 
Jaskulska, A., Font Nájera, A., Czarny, P., Serwecińska, L. & Mankiewicz-boczek, J. Daily dynamic of transcripts abundance of Ma-LMM01-like cyanophages in two lowland European reservoirs. Ecohydrol. Hydrobiol. 21, 543–548 (2021).
Google Scholar 
Gągała, I. et al. Role of environmental factors and toxic genotypes in the regulation of microcystins-producing cyanobacterial blooms. Microb. Ecol. 67, 465–479 (2014).PubMed 

Google Scholar 
Klindworth, A. et al. Evaluation of general 16S ribosomal RNA gene PCR primers for classical and next-generation sequencing-based diversity studies. Nucleic Acids Res. 41, 1–11 (2013).
Google Scholar 
Illumina. 16S Metagenomic Sequencing Library Preparation. (2013). http://support.illumina.com/content/dam/illumina-support/documents/documentation/chemistry_documentation/16s/16s-metagenomic-library-prep-guide-15044223-b.pdf.Frangeul, L. et al. Highly plastic genome of Microcystis aeruginosa PCC 7806, a ubiquitous toxic freshwater cyanobacterium. BMC Genomics 9, 1–20 (2008).
Google Scholar 
Hammer, Ø., Harper, D. A. T. & Ryan, P. D. Past: Paleontological statistics software package for education and data analysis even a cursory glance at the recent paleontological literature should convince anyone tha. Palaeontol. Electron. 4, 1–9 (2001).
Google Scholar 
Suzuki, M. T., Taylor, L. T. & DeLong, E. F. Quantitative analysis of small-subunit rRNA genes in mixed microbial populations via 5’-nuclease assays. Appl. Environ. Microbiol. 66, 4605–4614. https://doi.org/10.1128/AEM.66.11.4605-4614.2000 (2000).ADS 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Neilan B. A et al. rRNA sequences and evolutionary relationships among toxic and nontoxic cyanobacteria of the genus Microcystis Int J Syst Bacteriol 47(3), 693–697, https://doi.org/10.1099/00207713-47-3-693 (1997).
Google Scholar  More

Vega Thurber, R. et al. Deciphering coral disease dynamics: Integrating host, microbiome, and the changing environment. Front. Ecol. Evol. 2020, 8 (2020).
Google Scholar 
Groner, M. L. et al. Managing marine disease emergencies in an era of rapid change. Philos. Trans. R. Soc. Lond. B Biol. Sci. 371, 1689 (2016).
Google Scholar 
Richardson, L. L. Coral diseases: What is really known?. Trends Ecol. Evol. 13, 438–443 (1998).CAS 
PubMed 

Google Scholar 
Miller, M. W., Lohr, K. E., Cameron, C. M., Williams, D. E. & Peters, E. C. Disease dynamics and potential mitigation among restored and wild staghorn coral, Acropora cervicornis. PeerJ https://doi.org/10.7287/peerj.preprints.328 (2014).Article 
PubMed 
PubMed Central 

Google Scholar 
Teplitski, M. & Ritchie, K. How feasible is the biological control of coral diseases?. Trends Ecol. Evol. 24, 378–385 (2009).PubMed 

Google Scholar 
Zhu, Y. et al. Genetic diversity and disease control in rice. Nature 406, 718–722 (2000).ADS 
CAS 
PubMed 

Google Scholar 
Altermatt, F. & Ebert, D. Genetic diversity of Daphnia magna populations enhances resistance to parasites. Ecol. Lett. 11, 918–928 (2008).PubMed 

Google Scholar 
Aronson, R. B. & Precht, W. F. White-band disease and the changing face of Caribbean coral reefs. In (ed Porter, J. W.) The Ecology and Etiology of Newly Emerging Marine Diseases 25–38 (Springer Netherlands, 2001).Ruiz-Moreno, D. et al. Global coral disease prevalence associated with sea temperature anomalies and local factors. Dis. Aquat. Organ. 100, 249–261 (2012).PubMed 

Google Scholar 
Precht, W. F., Gintert, B. E., Robbart, M. L., Fura, R. & van Woesik, R. Unprecedented disease-related coral mortality in Southeastern Florida. Sci. Rep. 6, 31374 (2016).ADS 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Gignoux-Wolfsohn, S. A., Marks, C. J. & Vollmer, S. V. White Band Disease transmission in the threatened coral, Acropora cervicornis. Sci. Rep. 2, 804 (2012).ADS 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Aronson, R., Bruckner, A., Moore, J., Precht, B. & Weil, E. Acropora cervicornis. IUCN Red List of Threatened Species https://doi.org/10.2305/iucn.uk.2008.rlts.t133381a3716457.en (2008).Alvarez-Filip, L., González-Barrios, F. J., Pérez-Cervantes, E., Molina-Hernández, A. & Estrada-Saldívar, N. Stony coral tissue loss disease decimated Caribbean coral populations and reshaped reef functionality. Commun. Biol. 5, 440 (2022).PubMed 
PubMed Central 

Google Scholar 
Heres, M. M., Farmer, B. H., Elmer, F. & Hertler, H. Ecological consequences of Stony Coral Tissue Loss Disease in the Turks and Caicos Islands. Coral Reefs 40, 609–624 (2021).
Google Scholar 
Neely, K. L., Shea, C. P., Macaulay, K. A., Hower, E. K. & Dobler, M. A. Short- and long-term effectiveness of coral disease treatments. Front. Mar. Sci. 2021, 8 (2021).
Google Scholar 
Neely, K. L., Macaulay, K. A., Hower, E. K. & Dobler, M. A. Effectiveness of topical antibiotics in treating corals affected by stony coral tissue loss disease. PeerJ 8, e9289 (2020).PubMed 
PubMed Central 

Google Scholar 
Shilling, E. N., Combs, I. R. & Voss, J. D. Assessing the effectiveness of two intervention methods for stony coral tissue loss disease on Montastraea cavernosa. Sci. Rep. 11, 8566 (2021).CAS 
PubMed 
PubMed Central 

Google Scholar 
Walker, B. K., Turner, N. R., Noren, H. K. G., Buckley, S. F. & Pitts, K. A. Optimizing stony coral tissue loss disease (SCTLD) intervention treatments on Montastraea cavernosa in an Endemic Zone. Front. Mar. Sci. 8, 666224 (2021).
Google Scholar 
Forrester, G. E., Arton, L., Horton, A., Nickles, K. & Forrester, L. M. Antibiotic treatment ameliorates the impact of stony coral tissue loss disease (SCTLD) on coral communities. Front. Mar. Sci. 2022, 9 (2022).
Google Scholar 
Lee-Hing, C. et al. Management responses in Belize and Honduras, as stony coral tissue loss disease expands its prevalence in the Mesoamerican reef. Front. Mar. Sci. 9, 1 (2022).ADS 

Google Scholar 
Young, C. N., Schopmeyer, S. A. & Lirman, D. A review of reef restoration and coral propagation using the threatened genus Acropora in the Caribbean and Western Atlantic. Bull. Mar. Sci. 88, 1075–1098 (2012).
Google Scholar 
Lirman, D. & Schopmeyer, S. Ecological solutions to reef degradation: Optimizing coral reef restoration in the Caribbean and Western Atlantic. PeerJ 4, e2597 (2016).PubMed 
PubMed Central 

Google Scholar 
Baums, I. B. et al. Considerations for maximizing the adaptive potential of restored coral populations in the western Atlantic. Ecol. Appl. 29, e01978 (2019).PubMed 
PubMed Central 

Google Scholar 
Rosales, S. M. et al. Microbiome differences in disease-resistant vs susceptible Acropora corals subjected to disease challenge assays. Sci. Rep. 9, 18279 (2019).ADS 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Pinzón, C. J. H., Beach-Letendre, J., Weil, E. & Mydlarz, L. D. Relationship between phylogeny and immunity suggests older caribbean coral lineages are more resistant to disease. PLoS ONE 9, e104787. https://doi.org/10.1371/journal.pone.0104787 (2014).Article 
ADS 
CAS 

Google Scholar 
Drury, C. et al. Genomic patterns in Acropora cervicornis show extensive population structure and variable genetic diversity. Ecol. Evol. 7, 6188–6200 (2017).PubMed 
PubMed Central 

Google Scholar 
Maneval, P., Jacoby, C. A., Harris, H. E. & Frazer, T. K. Genotype, nursery design, and depth influence the growth of Acropora cervicornis fragments. Front. Mar. Sci. 8, 1 (2021).
Google Scholar 
Wright, R. M. et al. Intraspecific differences in molecular stress responses and coral pathobiome contribute to mortality under bacterial challenge in Acropora millepora. Sci. Rep. 7, 2609 (2017).ADS 
PubMed 
PubMed Central 

Google Scholar 
Vollmer, S. V. & Kline, D. I. Natural disease resistance in threatened staghorn corals. PLoS ONE 3, e3718 (2008).ADS 
PubMed 
PubMed Central 

Google Scholar 
Miller, N., Maneval, P., Manfrino, C., Frazer, T. K. & Meyer, J. L. Spatial distribution of microbial communities among colonies and genotypes in nursery-reared Acropora cervicornis. PeerJ 8, e9635 (2020).PubMed 
PubMed Central 

Google Scholar 
Klinges, G., Maher, R. L., Vega-Thurber, R. L. & Muller, E. M. Parasitic, “Candidatus Aquarickettsia rohweri” is a marker of disease susceptibility in Acropora cervicornis but is lost during thermal stress. Environ. Microbiol. 22, 5341–5355 (2020).CAS 
PubMed 
PubMed Central 

Google Scholar 
Miller, M. W. et al. Genotypic variation in disease susceptibility among cultured stocks of elkhorn and staghorn corals. PeerJ 7, e6751 (2019).PubMed 
PubMed Central 

Google Scholar 
Rohr, J. R. et al. Towards common ground in the biodiversity-disease debate. Nat. Ecol. Evol. 4, 24–33 (2020).PubMed 

Google Scholar 
Shearer, T. L., Porto, I. & Zubillaga, A. L. Restoration of coral populations in light of genetic diversity estimates. Coral Reefs 28, 727–733 (2009).ADS 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Ostfeld, R. S. & Keesing, F. Biodiversity and disease risk: The case of lyme disease. Conserv. Biol. 14, 722–728 (2000).
Google Scholar 
Lively, C. M. The effect of host genetic diversity on disease spread. Am. Nat. 175, E149–E152 (2010).PubMed 

Google Scholar 
Ostfeld, R. S. & Keesing, F. Effects of host diversity on infectious disease. Annu. Rev. Ecol. Evol. Syst. 43, 157–182 (2012).
Google Scholar 
King, K. C. & Lively, C. M. Does genetic diversity limit disease spread in natural host populations?. Heredity 109, 199–203 (2012).CAS 
PubMed 
PubMed Central 

Google Scholar 
Acevedo-Whitehouse, K., Gulland, F., Greig, D. & Amos, W. Inbreeding: Disease susceptibility in California sea lions. Nature 422, 35 (2003).ADS 
CAS 
PubMed 

Google Scholar 
O’Brien, S. J. et al. Genetic basis for species vulnerability in the cheetah. Science 227, 1428–1434 (1985).ADS 
PubMed 

Google Scholar 
Pearman, P. B. & Garner, T. W. J. Susceptibility of Italian agile frog populations to an emerging strain of Ranavirus parallels population genetic diversity. Ecol. Lett. 8, 401–408 (2005).
Google Scholar 
Reber, A., Castella, G., Christe, P. & Chapuisat, M. Experimentally increased group diversity improves disease resistance in an ant species. Ecol. Lett. 11, 682–689 (2008).PubMed 

Google Scholar 
Mundt, C. C. Use of multiline cultivars and cultivar mixtures for disease management. Annu. Rev. Phytopathol. 40, 381–410 (2002).CAS 
PubMed 

Google Scholar 
Elton, C. S. The Ecology of Invasions by Animals and Plants (University of Chicago Press, 2000).
Google Scholar 
Schopmeyer, S. A. et al. Regional restoration benchmarks for Acropora cervicornis. Coral Reefs 36, 1047–1057 (2017).ADS 

Google Scholar 
Baums, I. B., Miller, M. W. & Hellberg, M. E. Geographic variation in clonal structure in a reef-building Caribbean coral, Acropora palmata. Ecol. Monogr. 76, 503–519 (2006).
Google Scholar 
Gignoux-Wolfsohn, S. A., Precht, W. F., Peters, E. C., Gintert, B. E. & Kaufman, L. S. Ecology, histopathology, and microbial ecology of a white-band disease outbreak in the threatened staghorn coral Acropora cervicornis. Dis. Aquat. Organ. 137, 217–237 (2020).PubMed 

Google Scholar 
Gignoux-Wolfsohn, S. A. & Vollmer, S. V. Identification of candidate coral pathogens on white band disease-infected staghorn coral. PLoS ONE 10, e0134416 (2015).PubMed 
PubMed Central 

Google Scholar 
Brooks, M. et al. GlmmTMB balances speed and flexibility among packages for zero-inflated generalized linear mixed modeling. R J. 9, 378–400 (2017).
Google Scholar 
Fox, J. & Weisburg, S. An R Companion to Applied Regression 3rd edn. (Sage, 2019).
Google Scholar  More

To our knowledge, this is the only work done on the terrestrial biodiversity status in the direct vicinity of a limited duration volcanic eruption. In this contribution, we document and assess the impact on the main plant and animal groups within the ecosystems during a volcanic eruption (Table 1). While some groups were clearly disadvantaged: ferns and herbaceous plants as well as invertebrates and saurians (lizards and geckos); other groups such as conifers and woody shrubs showed better resilience, as did the birds.This study is particularly important because of its location in a Mediterranean biodiversity hotspot13,14, harbouring a unique ecosystem of oceanic island organisms (38% of the Canary archipelago endemicity). Islands indeed exhibit a disproportionate amount of the world’s biodiversity but unfortunately a high number of extinctions have also occurred there14. The biodiversity in the south of the island is poorer than in the north. This is probably explained in part by the relatively frequent volcanic activity featuring seven major eruptions since 1585, including this one in 2021 (see15), which led to alternating destruction and neo-colonization processes.Concerning the flora, the Canary pine forest was the most affected ecosystem and vegetation type, as it is dominant in the vicinity of the new volcanic vents. The southern slopes of this forest were the most disturbed area due to the location of the volcano, combined with the prevailing northeasterly trade winds (Fig. 1). Tephra fallout and sulphurous gases were the main factors that affected the pine forest, over a vast surface area. Furthermore, the local xerophytic and thermophilous habitats also lost much of their surface area. In contrast to the pine forest, this drastic reduction was caused by the progressive downslope expansion of the lava flows.The Canary Island pine was thus notably affected by tephra fall, sulphuric acid aerosol12, and short episodes of acid rain. However, this conifer shows high resistance to temperature, confirming its great adaptation to volcanic events16, which is probably also one of the keys to its resistance to the more frequent present-day wildfires17. This pine species has evolved among volcanoes for the last 13 My16 and has adapted successfully to high temperatures. Moreover, thunderstorms with lightning occur in the Canaries together with abundant rainfall; consequently, wild forest fires should presumably not have been so frequent in the island’s past, before human colonization. In this habitat it is also remarkable that epiphytic lichens (U. articulata) apparently resisted on the pines until the 12th week, considering their high sensitivity to anthropogenic pollution18.The life cycle of flowering plants was drastically disrupted due to all the above factors, with great impact on foliage, photosynthesis, and growth. However, soil changes due to the deposition of tephra and its lixiviation by rain is one of the most dramatic factors affecting plants and a long-term impact of volcanic eruptions19. The nearest individuals to the crater were most directly affected by intense tephra falls and concentrated volcanic gases (SO2, HCl, HF, CO2). However, plants located in the nearest 200 m to the lava flows but at more than 2 km from the crater were presumably more disturbed by the high temperature of the slow-cooling lava and its lesser gas emissions.Large woody plants exhibited a better frequency of survival than smaller ones in the face of this extreme stress (Table S1 and19). In the Hekla area (Iceland), most trees have thickened trunks, indicating that those trees that survive have had a long life subjected to frequent volcanic damage19. Secondary woodiness of island plants (sensu20) has been traditionally related to drought20,21, ecological shift22 or a counter-selection of inbreeding depression in founding island populations23. However, this adaptation also favours the resistance of many shrubby plants to high temperatures close to craters and lava flows but primarily their resistance to the intense tephra falls that affect a much larger area. In addition, plant and stem height plays a fundamental role in overcoming the deep layers of deposits. This latter effect was particularly important up to 2.5 km from the crater (tephra thickness  > 30 cm) (Figs. 1 and 2), as the herbaceous plants were completely buried, sometimes to more than 1.5 m depth. Therefore, the seed bank has also probably been rendered largely non-functional. However, deposits were recorded over almost the whole island, indicating that longer lasting or more intense eruptions would severely affect an even larger area. Such events have been hitherto ignored in the intensely discussed “island woodiness” debate21,23,24,25,26,27. We found surviving populations of endemic woody taxa heavily impacted by tephra deposits close to lava flows, across a wide range of genera such as Rumex (R. lunaria), Echium (E. brevirame), Euphorbia (E. lamarckii, E. canariensis and E. balsamifera), Aeonium (A. davidbramwellii), Rubia (R. fruticosa), Schizogyne (S. sericea), Carlina (C. falcata) or Sonchus (S. hierrensis) (Table S2), which coincide with the general list of woody Canary plants20. Most members of these genera in other ecosystems on continents are mainly herbaceous. As such eruptions and their impacts due to ash depositions are frequent events on volcanic islands, e.g. several times within a century on La Palma, this is a “frequent” selective process at evolutionary time scales.With regard to the fauna, the invertebrate community collapsed during the first two weeks (Table S2), probably due to rapid deterioration of the growth state of plants. These changes in the invertebrates were caused by the tephra contacting the cuticular lipid layer28 and water loss due to tegument abrasion29. In this period, many insect pests (especially whitefly pupae) in banana plantations (farmers’ observations) were drastically reduced. This sudden decrease in insect populations affected the whole food web and probably caused part of the ecological collapse of saurian and some passerine communities30. In the case of lizards, smaller individuals seem to resist the adverse conditions better than large ones, as observed in other eruptions3. This could be linked to their lower food requirements and greater ease in finding refuges. Loss of body condition in lizards post-eruption has been recorded and negatively affects reproduction quality31. However, some lizards have shown a good ability to find food in the tephra substrate32. We found abundant tephra particles in some vertebrate droppings (lizards, birds, and mammals) during the eruption, probably involuntarily ingested. At least in bats, ingestion during feeding produces physiological stress that is likely related to baldness, high ectoparasite loads or possible mineral deficiencies33.As described in the Canary Islands, some passerines show high fidelity to their territories (see34). During the eruption, Sardinian warblers (Curruca melanocephala) maintained their territories until the imminent arrival of lava flows. Larger birds (kestrels F. tinnunculus, ravens C. corax and buzzards B. buteo) were well able to continue flying in the areas surrounding the crater. Furthermore, some cases like F. tinnunculus showed great feeding plasticity in the first couple of weeks. At least six times, kestrels tried to catch birds (especially small passerines and doves), contrary to their usual diet based on abundant lizards and insects35. Widening of trophic niches in island organisms has traditionally been interpreted as linked to disharmony in island ecosystems36,37,38. However, this plasticity is tremendously beneficial in ecological catastrophes, where food becomes exceptionally scarce. In the case of bats, their flight is limited by the delicate structure of their patagium, which can be damaged by the frequent pyroclastic tephra fall. Furthermore, scarcity of insects in the first few kilometres from the crater probably led to their displacement to other more distant and richer food resource zones.As we learned from the movement capacity of the vertebrate animals that still inhabited the affected area, those with greater mobility, birds and bats, resisted the eruptive process much better than those with less mobility, e.g. saurians.Lastly, during this destructive event on La Palma, we had the opportunity to increase our knowledge of how ecological-evolutionary adaptations have favoured the survival of insular organisms. Such responses are traditionally mentioned in the context of island biology. As already mentioned, one of the most interesting findings verifies the remarkable adaptation of Canary Island pine trees (P. canariensis) to volcanism (see16), including extremely harsh ecological conditions. Other insular trends related to the prevalence of woodiness in insular flowering plants20,21, or the high trophic plasticity of some vertebrates on oceanic islands36, have not previously been associated with their potential evolution along with volcanic processes. However, such evolutionary adaptations most likely played an important role in the survival of plants and animals affected by the volcano. For this reason, it is worth considering and debating whether these previously mentioned evolutionary processes are in fact also linked to repeated volcanic episodes on oceanic islands. More