Ecology

Abiotic parametersThe temperature regime experienced by the embryos was purposefully natural and therefore varied between the three air exposure treatments. The subtidal treatment, where embryos were continuously submerged in water, remained around 9.5 °C for the duration of the experiment, while the intertidal treatments experienced dips in temperature during outside air exposure down to 2.5 °C and 0.8 °C, for the low and high intertidal respectively (Fig. 3A). Accumulated thermal units (ATU; days × temperature post collection until hatch) for each air exposure treatment were 79.6, 75.4 and 65.3 for subtidal, low intertidal and high intertidal, respectively. Despite differences in thermal regime, peak hatch was on the same day (March 14, 2021) for all air exposure and CO2 treatments, estimated at 11 dpf.Figure 3Temperature and pH experienced by the herring embryos and larvae throughout the experiment. Hourly measurements of (A) air/water temperature experienced by herring embryos in each of the tidal treatments (subtidal: continuous immersion in 9.5 °C water; low intertidal: 2 × 2 h air exposure; high intertidal: 5 + 9 h air exposure) and (B) pH levels in the tanks for each of the CO2 treatments (greens = control, 400 µatm CO2, yellows = medium, 1600 µatm CO2, reds = high, 3000 µatm CO2); dots are pH levels measured in the individual jars during larval incubation.Full size imageThe pH levels in the tanks were measured hourly and were stable over the course of the embryonic incubation period, with no overlap between treatments, although there was some overlap between individual jars. Control treatment was consistently around a pH of 8, the medium treatment had a pH of 7.4 and the high CO2 treatment had a mean pH of 7.1 (Table 1). After hatch, when the larvae were transferred to the jars, circulation and gas exchange between jars and tank were not as high and CO2 accumulated in the jars over time, leading to pH levels deviating from tank pH levels (Fig. 3B). Although oxygen levels remained high (7–9 mg/L), the pH dropped from a mean 8–7.6 in the control on two occasions, and was brought back up with a partial water exchange from the incubation tank water. The pH in the medium and high CO2 treatments were not as affected (Fig. 3B), however, final water chemistry measurements after completion of the experiment (2 days post water exchange) revealed much higher CO2 levels in all treatments (Table 1: day 15).Table 1 Mean water parameters for each treatment (mean of 3 tanks ± S.D.) at the beginning (day 1, 2021-03-06) and end (day 6, 2021-03-12) of embryonic incubation and mean parameters in the jars (N = 9) at the end of larval incubation (day 15, 2021-03-19); Temperature (T), salinity, pCO2, total CO2 (TCO2) measured at distinct sampling intervals with the BoL; total alkalinity (TA) and pH (on the total scale) calculated with CO2SYS.Full size tableEffect of air exposure and CO2 treatment during embryonic developmentNeither embryonic survival nor growth were significantly affected by treatment in our experiment. Percent daily embryonic mortality was low and not significantly affected by CO2 treatment or air exposure (CO2: p = 0.088, F2 = 2.45; Tide: p = 0.11, F2 = 2.19; CO2*Tide: p = 0.18, F2 = 1.59) . Egg diameter at 6 dpf was also not significantly affected by treatment (CO2: p = 0.38, X2 (2, N = 30) = 1.92; Tide: p = 0.83, X2 (2, N = 30) = 0.33; CO2*Tide: p = 0.08, X2 (2, N = 30) = 8.25). Metabolic rate, as indicated by embryonic heart rate, was significantly affected by air exposure at 6 dpf (p  *; 0.1  >).Full size image More

Pfennig, D. The adaptive significance of an environmentally-cued developmental switch in an anuran tadpole. Oecologia 85, 101–107 (1990).ADS 
PubMed 
Article 

Google Scholar 
Brönmark, C. & Miner, J. G. Predator-induced phenotypical change in body morphology in Crucian carp. Science 258, 1348–1350 (1992).ADS 
PubMed 
Article 

Google Scholar 
Wikelski, M. & Thom, C. Marine iguanas shrink to survive El Niño. Nature 403, 37–38 (2000).ADS 
CAS 
PubMed 
Article 

Google Scholar 
Agrawal, A. A. Phenotypic plasticity in the interactions and evolution of species. Science 294, 321–326 (2001).ADS 
CAS 
PubMed 
Article 

Google Scholar 
Huchard, E., English, S., Bell, M. B. V., Thavarajah, N. & Clutton-Brock, T. Competitive growth in a cooperative mammal. Nature 533, 532–534 (2016).ADS 
CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Travis, J. Evaluating the adaptive role of morphological plasticity. In: Ecological morphology (pp. 99–122) (The University of Chicago Press, Chicago, 1994).Lázaro, J., Dechmann, D. K. N., LaPoint, S., Wikelski, M. & Hertel, M. Profound reversible seasonal changes of individual skull size in a mammal. Curr. Biol. 27, R1106–R1107 (2017).PubMed 
Article 
CAS 

Google Scholar 
Lázaro, J. & Dechmann, D. K. Dehnel’s phenomenon. Ecol. Evol. 31, R463–R465 (2021).
Google Scholar 
Bronstein, J. L. The evolution of facilitation and mutualism. J. Ecol. 97, 1160–1170 (2009).Article 

Google Scholar 
Leigh, J. The evolution of mutualism. J. Environ. Biol. 23, 2507–2528 (2010).
Google Scholar 
Liu, C., Yang, D. R. & Peng, Y. Q. Body size in a pollinating fig wasp and implications for stability in a fig-pollinator mutualism. Entomol. Exper. Appl. 138, 249–255 (2011).Article 

Google Scholar 
Pires, M. M., Guimarães, P. R., Galetti, M. & Jordano, P. Pleistocene megafaunal extinctions and the functional loss of long-distance seed-dispersal services. Ecography 41, 153–163 (2018).Article 

Google Scholar 
Boucher, D., James, S. & Keeler, K. The ecology of mutualism. Annu. Rev. Ecol. Syst. 13, 315–347 (1982).Article 

Google Scholar 
Irwin, R. E. & Brody, A. K. Nectar robbing in Ipomopsis aggregata: effects on pollinator behavior and plant fitness. Oecologia 116, 519–527 (1998).ADS 
PubMed 
Article 

Google Scholar 
Allen, G. The Anemonefishes: their classification and biology (T.F.H. Publications, 1972).
Google Scholar 
Fautin, D.G. & Allen, G.R. Field guide to anemonefishes and their host sea anemones. (Western Australian Museum, Perth, 1992).Ollerton, J., McCollin, D., Fautin, D. G. & Allen, G. R. Finding NEMO: nestedness engendered by mutualistic organization in anemonefish and their hosts. Proc. R. Soc. B Biol. Sci. 274, 591–598 (2006).Article 

Google Scholar 
Fricke, H. & Fricke, S. Monogamy and sex change by aggressive dominance in coral reef fish. Nature 266, 830–832 (1977).ADS 
CAS 
PubMed 
Article 

Google Scholar 
Buston, P. M. Size and growth modification in clownfish. Nature 424, 145–146 (2003).ADS 
CAS 
PubMed 
Article 

Google Scholar 
Mariscal, R. N. The nature of the symbiosis between Indo-Pacific anemone fishes and sea anemones. Mar. Biol. 6, 58–65 (1970).Article 

Google Scholar 
Elliott, J. K., Elliott, J. M. & Mariscal, R. N. Host selection, location, and association behaviors of anemonefishes in field settlement experiments. Mar. Biol. 122, 377–389 (1995).Article 

Google Scholar 
Verde, A. E., Cleveland, A. & Lee, R. W. Nutritional exchange in a tropical tripartite symbiosis II: direct evidence for the transfer of nutrients from host anemone and zooxanthellae to anemonefish. Mar. Biol. 162, 2409–2429 (2015).Article 
CAS 

Google Scholar 
Cleveland, A., Verde, E. A. & Lee, R. W. Nutritional exchange in a tropical tripartite symbiosis: direct evidence for the transfer of nutrients from anemonefish to host anemone and zooxanthellae. Mar. Biol. 158, 589–602 (2011).Article 

Google Scholar 
Sale, P. F. Effect of cover on agonistic behavior of a reef fish: a possible spacing mechanism. Ecology 53, 753–758 (1972).Article 

Google Scholar 
Fricke, H. W. & Holzberg, S. Social units and hermaphroditism in a pomacentrid fish. Naturwissenschaften 61, 367–368 (1974).ADS 
Article 

Google Scholar 
Fricke, H. W. Control of different mating systems in a coral reef fish by one environmental factor. Anim. Behav. 28, 561–569 (1980).Article 

Google Scholar 
Mitchell, J. S. & Dill, L. M. Why is group size correlated with the size of the host sea anemone in the false clown anemonefish?. Canad. J. Zool. 83, 372–376 (2005).Article 

Google Scholar 
Chausson, J., Srinivasan, M. & Jones, G. P. Host anemone size as a determinant of social group size and structure in the orange clownfish (Amphiprion percula). PeerJ 6, e5841 (2018).PubMed 
PubMed Central 
Article 

Google Scholar 
Reed, C., Branconi, R., Majoris, J., Johnson, C. & Buston, P. Competitive growth in a social fish. Biol. Lett. 15, 20180737 (2019).PubMed 
PubMed Central 
Article 

Google Scholar 
Buston, P. M. Mortality is associated with social rank in the clown anemonefish (Amphiprion percula). Mar. Biol. 143, 811–815 (2003).Article 

Google Scholar 
Branconi, R. et al. Ecological and social constraints combine to promote evolution of non-breeding strategies in clownfish. Comm. Biol. 3, 1–7 (2020).Article 
CAS 

Google Scholar 
Schmiege, P. F., D’Aloia, C. C. & Buston, P. M. Anemonefish personalities influence the strength of mutualistic interactions with host sea anemones. Mar. Biol. 164, 24 (2017).Article 

Google Scholar 
Barbasch, T. A. & Buston, P. M. Plasticity and personality of parental care in the clown anemonefish. Anim. Behav. 136, 65–73 (2018).Article 

Google Scholar 
Abramoff, M. D., Magalhaes, P. J. & Ram, S. J. Image PROcessing with ImageJ. Biophoto. Int. 11, 36–42 (2004).
Google Scholar 
R Core Team. R: A language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria. https://www.R-project.org/ (2020).Goodrich, B., Gabry, J., Ali I. & Brilleman, S. Rstanarm: Bayesian applied regression modeling via Stan. R package version 2.21.1 https://mc-stan.org/rstanarm (2020).Weatherley, A. H. Approaches to understanding fish growth. Trans. Am. Fish. Soc. 119, 662–672 (1990).Article 

Google Scholar 
Gabry, J. Shinystan: interactive visual and numerical diagnostics and posterior analysis for Bayesian models. R package version 2.5.0. https://CRAN.R-project.org/package=shinystan (2018).Gelman, A., Goodrich, B., Gabry, J. & Vehtari, A. R-squared for Bayesian regression models. Am. Stat. 3, 307–309 (2018).MathSciNet 

Google Scholar 
Lüdecke, D., Ben-Shachar, M. S., Patil, I., Waggoner, P. & Makowski, D. Performance: an R package for assessment, comparison and testing of statistical models. J. Open Sour. Softw. 6, 60 (2021).
Google Scholar 
Gabry, J., Simpson, D., Vehtari, A., Betancourt, M. & Gelman, A. Visualization in Bayesian workflow. J. R. Stat. Soc. Ser. A Stat. Soc. 182, 389–402 (2019).MathSciNet 
Article 

Google Scholar 
Gabry, J. & Mahr, T. Bayesplot: plotting for bayesian models. R package version 1.8.0. https://mc-stan.org/bayesplot/ (2021).Elliott, J. K. & Mariscal, R. N. Coexistence of nine anemonefish species: differential host and habitat utilization, size and recruitment. Mar. Biol. 138, 23–36 (2001).Article 

Google Scholar 
Buston, P. M. Forcible eviction and prevention of recruitment in the clown anemonefish. Behav. Ecol. 14, 576–582 (2003).Article 

Google Scholar 
Fautin, D. G. & Allen, G. R. Anemone fishes and their host sea anemones: a guide for aquarists and divers. Sea Challengers (1997).Beldade, R., Blandin, A., O’Donnell, R. & Mills, S. C. Cascading effects of thermally-induced anemone bleaching on associated anemonefish hormonal stress response and reproduction. Nat. Commun. 8, 1–9 (2017).ADS 
CAS 
Article 

Google Scholar 
Cortese, D. et al. Physiological and behavioural effects of anemone bleaching on symbiont anemonefish in the wild. Funct. Ecol. 35, 663–674 (2021).Article 

Google Scholar 
Scherbatskoy, E. C. et al. Characterization of a novel picornavirus isolated from moribund aquacultured clownfish. J. Gen. Virol. 101, 735–745 (2020).CAS 
PubMed 
Article 

Google Scholar 
Saenz-Agudelo, P., Jones, G. P., Thorrold, S. R. & Planes, S. Mothers matter: contribution to local replenishment is linked to female size, mate replacement and fecundity in a fish metapopulation. Mar. Biol. 162, 3–14 (2014).Article 

Google Scholar 
Barbasch, T. A. et al. Substantial plasticity of reproduction and parental care in response to local resource availability in a wild clownfish population. Oikos 129, 1844–1855 (2020).Article 

Google Scholar 
Sebens, K. P. The ecology of indeterminate growth in animals. A. Rev. Ecol. Syst. 18, 371–407 (1987).Article 

Google Scholar 
Buston, P. M. & García, M. B. An extraordinary life span estimate for the clown anemonefish Amphiprion percula. J. Fish Biol. 70, 1710–1719 (2007).Article 

Google Scholar 
Chamberlain, S. A., Kilpatrick, J. R. & Holland, J. N. Do extrafloral nectar resources, species abundances, and body sizes contribute to the structure of ant–plant mutualistic networks?. Oecologia 164, 741–750 (2010).ADS 
PubMed 
Article 

Google Scholar 
Marting, P. R., Kallman, N. M., Wcislo, W. T. & Pratt, S. C. Ant-plant sociometry in the Azteca-Cecropia mutualism. Sci. Rep. 8, 1–15 (2018).Article 
CAS 

Google Scholar 
Fordyce, J. A. The evolutionary consequences of ecological interactions mediated through phenotypic plasticity. J. Exp. Biol. 209, 2377–2383 (2006).PubMed 
Article 

Google Scholar 
West-Eberhard, M. J. Developmental plasticity and evolution (Oxford University Press, 2003).Book 

Google Scholar 
West-Eberhard, M. J. Phenotypic accommodation: adaptive innovation due to developmental plasticity. J. Exp. Zool. B Mol. Develop. Evol. 304, 610–618 (2005).Article 

Google Scholar 
Moczek, A. P. et al. The role of developmental plasticity in evolutionary innovation. Proc. R. Soc. B Biol. Sci. 278, 2705–2713 (2011).Article 

Google Scholar  More

This section is devoted to show results and to highlight eventual effects of the interplay between the nonlinearity characterizing the system dynamics and the presence of noisy fluctuations for the irradiance variable.Analysis of experimental dataThe need of taking into account noisy fluctuations of such an environmental variable is well demonstrated in Fig. 1. In the first panel (a) the experimental time behaviour of the irradiance is shown. This noisy curve is based on the experimental data (purple points) of the Boussole buoy located in the Gulf of Lion, collected over a period of nine years, precisely from 2004 to 2013. The time series of the experimental data presents quite a few gaps in time due to the malfunction of the buoy. This aspect has been remedied by merging the experimental data with those of the OASIM model validated for the Boussole site61 (yellow points). The latter is a multispectral atmospheric radiative transfer model that is in turn forced by experimental-model data based on ECMWF ERAINTERIM reanalyses which provide, for example, cloud cover data. The radiative model is partly stochastic since it considers the effects stemming from the presence of clouds, averaged along a single day (this explains why the yellow points are slightly less scattered). We see that the OASIM model accurately reproduces the profile which emerges from the experimental data. Further, we stress that the experimental data are only used in this initial analysis. In the biogeochemical simulations the irradiance signal is fully reconstructed starting from a realistic seasonal cycle combined with a range of different random fluctuations, and the information from OASIM is not used. In the second panel (b) the daily (black points) as well as the three-month (red points) running mean of the experimental series are plotted. Figure 1c shows the irradiance noisy fluctuations (INF) which have been obtained by subtracting the three-month running mean curve (3MRM, red curve in Fig. 1b) from the daily running mean one (DRM, black curve in Fig. 1b) and normalizing with respect to the mean of the 3MRM ((overline{3MRM})), namely (INF = (DRM – 3MRM) / overline{3MRM}). We see that a seasonal overall trend with higher oscillations during the winter time can be seen, implying that the characteristics of the noise may change over the year. Moreover, a slight imbalance between positive and negative values of the noisy fluctuations (that is, different values of the maximum fluctuation intensity) is present. The physical reason for the occurrence of such an aspect can be ascribed to the fact that the maximum value of solar irradiance corresponds to that measured during a sunny day. Conversely, the minimum level tends to zero corresponding to a dense darkness. While the former is close to the mean value of the solar irradiance (most of all in summer), the latter is much further away and then a natural asymmetry arises in the random fluctuations. However, it should be noted that, apart from the intense spikes, the asymmetry is not so pronounced, as proved by the mean value (red line in Fig. 1c) which is practically zero, namely (0.4%) of the (overline{3MRM}). Therefore, basing on this last observation, to model the noise affecting the irradiance dynamics, as a first approximation we consider a symmetric Gaussian autocorrelated noise as described in the next subsection.On the basis of such experimental results, we postulate the hypothesis that random fluctuations of light cannot be neglected, most of all in the study of ecological systems where light profoundly determines the system dynamics, governing fundamental processes at the basis of of the food web.Figure 1(a) Experimental data (purple points) of the stochastic solar irradiance collected by the Boussole buoy in a time-window of 9 years (2004-2013); the yellow points are the data generated by the OASIM model used to fill the gaps present in the experimental time-series due to malfunctioning of the buoy. (b) Daily (black points) and three-month (red points) running mean of the light curve in panel (a). (c) Irradiance noisy fluctuations (INF), obtained by subtracting the three-month running mean curve (3MRM) from the daily running mean one (DRM) and normalizing with respect to the mean value of 3MRM ((overline{3MRM})), namely (INF = (DRM – 3MRM) / overline{3MRM}); the red line represents the mean value of such fluctuations. Data already presented and validated in61.Full size imageSolar irradianceThe solar irradiance forcing is derived considering a deterministic seasonal oscillation combined with an Ornstein-Uhlenbeck process. The coefficient of variation (CV) of simulated light forcing, Fig. 2, (CV=sigma / mu) ((mu) and (sigma) being mean value and standard deviation calculated over both time and numerical realizations), is shown for 231 (D-tau) pairs. D and (tau) represent the intensity of a Gaussian noise source and the auto-correlation time of the fluctuations, respectively (see Eqs. (2) and (3)).Each pixel represents the mean value on time of CV calculated with respect to 1000 different stochastic realizations. Figure 2Coefficient of variation ((CV=sigma / mu)) of irradiance resulting from numerical integration of model equations for 231 (D-tau) different scenarios.Full size imageIt is easy to see the agreement between the results obtained from the numerical integration and the theoretical ones derivable from Eq. (5) by putting (text {var}{F_L(0)}=0) and (t gg 1), getting (sigma ^2_L=D / 2tau). In Fig. 2, indeed, the maximum values of (sigma) lie in the upper left part of the plot corresponding to small (high) values of (tau) (D). As it is clear the values of D have been chosen in order to obtain a relative standard deviation ranging from (5%mu) to (60%mu). We underline that, in this case, it is possible to interchangeably consider (sigma) and CV since the dependence of CV on D and (tau) does not differ from that of (sigma) (meaning that the dependence of (sigma) is not altered by dividing by (mu)) (results not shown).Effects on population dynamicsIn this section the noise-induced effects on the population dynamics are examined. The nine planktonic populations present a different qualitative behaviour of the CV, compared to that of the irradiance. In this case, the CV is characterized by a strong non-monotonic dependence on the parameter (tau). This aspect can be appreciated in Fig. 3 where different curves of CV versus the time correlation parameter are shown for different fixed values of D.Figure 3Coefficient of variation ((CV=sigma / mu)) of the nine planktonic populations resulting from numerical integration of model equations plotted versus the considered values of (tau); the different curves are related to different values of the noise intensity D.Full size imageThe existence of a maximum value for CV can be appreciated for each species. Although the qualitative behaviour is the same for all strains, particular attention has to be payed on diatoms and nanoflagellates. All the other species, indeed, present a percent variation of standard deviation between (2%) and (15%). In the case of nanoflagellates, instead, the D-dependent range is (20-90%), while diatoms reach values over the (100%) for the highest values of D. Therefore, these two species, in particular, and the whole system, in general, are extremely sensitive to the auto-correlation time which characterizes the noise.We note that the different curves related to the different selected values of D approach the horizontal axis, tending asymptotically to vanish as (tau) increases. Such a behaviour can be explained by the fact that high values of (tau) give rise to a more correlated dynamics, so that (tau rightarrow infty) implies fully correlated time-behaviours corresponding to the deterministic case. In this instance, then, all the different realizations give the same results, making the standard deviation vanish. The same happens, independently of the value of (tau), for low values of noise intensity for which the corresponding curves approach the same almost vanishing value (see orange, gray and yellow lines). Differently from the previous case, when (tau rightarrow 0) the noise tends to a delta-correlated noise, that is a white noise; for (tau ne 0), instead, the noise spectrum is not flat, being characterized by a Cauchy-Lorentz distribution. The strong nonmonotonicity of CV with respect to (tau), emerging when there are relatively high values of CV, implies a greater variability of the system biomass. Lower values of CV indicate that the system dynamics is less influenced by the presence of noise where very little or no differences with respect to the deterministic case are present. Conversely, high values of CV clearly demonstrate the remarkable signature of the presence of an impacting noise source. It is interesting to note that the noise influence on the ecosystem strongly depends on both (tau) and D, that is, just an intense noise is not enough to generate a greater response of the ecosystem. In particular, experimental data are characterized by a CV approximately equal to 0.361, which corresponds to values of D and (tau) lying on the diagonal strip in Fig. 2 ranging from ((tau ,D)=(0.5,10^4)) to ((tau ,D)=(365,10^7)). Finally we note the presence of a noise suppression effect. High values of D, indeed, can generate slight effects when the correlation time (tau) does not take on suitable values.The results shown here are an extension of the previous work by Benincà et al.56. There, the authors analyse a simpler, less realistic model of two interacting populations, whose dynamics is affected by a randomly fluctuating temperature. In that case, moreover, the deterministic oscillations of the temperature are suppressed, and the system exhibits intrinsic Lotka-Volterra oscillations whose frequency match with the characteristic one(s) of the noise. On the contrary, here, the observed maximum response (see Fig. 3) cannot be interpreted as a synchronization effect, since our model does not present intrinsic Lotka-Volterra-like oscillations and the periodic population variability is only due to the deterministic forcing(s).The nonmonotonic behaviour of the CV can be then interpreted as the signature of the intimate interplay between the ecological system and the noise. This interplay, indeed, has a pivotal role in both determining the dynamics of the populations and defining the characteristics of the ecosystem.In Fig. 3 it can be observed that the value of (tau) for which CV is maximum strongly depends on the noise intensity D. In particular, it is possible to note that the peaks in Fig. 3 move towards higher values of (tau) as the noise intensity increases. Thus, Fig. 3 demonstrates that the maximum-response effect to the random fluctuations is sensitive to the noise intensity D.However, it is important to underline that the response of the system to the noisy signal does not depend on the yearly oscillations induced by the deterministic forcings. Indeed, by considering constant the deterministic part of all external forcings (temperature, irradiance, wind and salinity), the non monotonic behaviour of CV with respect to both (tau) and D is still present, provided that the populations are not extinct (plot not shown). In this scenario indeed, besides dinoflagellates, diatoms and nanoflagellates are practically extinct as well, exhibiting thus a constant vanishing variance. All the other strains, instead, present qualitatively the same nonmonotonicity with only slight differences (shift of the peaks and different mean values of the CV curves), probably due to the extinction of diatoms and nanoflagellates which causes relevant differences in the system dynamics. More specifically, the system’s response seems to depend on both the noise intensity and the correlation time (see Fig. 3).In this scenario (absence of seasonal driving) we have studied the dependence on both parameters D and (tau) of the probability density functions (PDFs) of the non-vanishing populations. In Fig. 4, the PDFs of bacteria (B1), picophytoplankton (P3), microzooplankton (Z5) and etherotrophic nanoflagellates (Z6) are plotted for (tau =0.5) and eight different values of the parameter D.Figure 4Dependence of the probability density functions of non-vanishing populations on the parameter D for (tau =0.5). The curves are normalized within the interval taken into account. For this reason the relative peaks of the curves in the bottom panels have different values compared to those of the top panels. However, the figure aims at showing the existence of the value of the noise intensity for which the system is more sensitive as well as the generation of a stationary out-of-equilibrium state induced by the noise.Full size imageWe see that the mean value and the variance of these populations are strongly affected by the presence of random fluctuations in the irradiance. Specifically, as the noise intensity increases the mean values of picophytoplankton and bacteria concentrations exhibit a shift. In particular, the results indicate that picophytoplankton is disavantaged by the presence of a noisy component in the irradiance, which indeed tends to inhibit its ability to absorbe the solar light, slowing down its growth. As a consequence, since phytoplankton and bacteria compete for the same resources, as the former declines the latter are favoured, with a compensation mechanism which allows their predators (zooplankton populations) to be almost not affected by the noisy behaviour of the irradiance. Further, we note that for intermediate values of the noise intensity ((D = 10^4 – 10^5)) a maximum of the variance occurs (the PDFs are clearly spread on a wider range of values). Such an effect indicates that the noisy behaviour of irradiance strongly influences the whole ecosystem dynamics. Moreover, the nonmonotonic behaviour of the variance (its PDFs become larger and then tighter again as the noise intensity increases) indicates that the noise pushes the ecosystem away from equilibrium, driving it towards a non-equilibrium steady state. Finally, we note that the nonmonotonic behaviour of CV as a function of the noise intensity remains also in the presence of seasonal driving.Figure 5Coefficient of variation ((CV=sigma / mu)) of nine planktonic populations resulting from numerical integration of model equations plotted versus the considered values of D; different curves correspond to different values of the correlation time (tau).Full size imageFigure 5 shows indeed the nonmonotonic response of the ecosystem to the change of D when the deterministic seasonal cycling of the four environmental parameters (temperature, irradiance, wind and salinity) is present. It is easy to observe that also in this instance the major noise-induced effect appears in nanoflagellates and diatoms with a percent standard deviation of 50(%) and 100(%), respectively. The coalescence of different curves (related to different values of (tau)), as D decreases, is due to the fact that for (D rightarrow 0) the impact of the noise is negligible and the evolution of the system practically resembles the deterministic one. On the contrary, for higher values of D remarkable differences arise and clear peaks of CV appear in the considered range of variation.These plots show that, for a fixed value of (tau), there exists a value of the noise intensity for which the planktonic concentrations are maximally spread around their mean values (corresponding to the maximum value of CV and then of the variance). Moreover, such a nonmonotonic behaviour suggests the presence of a resonance, which can be interpreted as the effect of the interplay between the nonlinearity of the system and the environmental random fluctuations.Also in this case, the interplay between the two parameters D and (tau) in determining and characterizing the dynamics of the ecosystem transparently emerges. The value of D corresponding to the maximum value of CV, indeed, basically depends on the specific value of (tau).Finally, we point out that the different dynamic scenarios identified by the D-(tau) couples can be experienced by the system during the year, since the two parameters may seasonally vary depending on the different weather conditions. In other words, a seasonally varying noise (see Fig. 1c) may cause the nine populations explore different regions of the D-(tau) space during the year. Therefore, the results reported in this paper can highlight the detectable yearly variability of a marine ecosystem which does not stem from the deterministic seasonal variation of environmental parameters.Effects on the organic carbonIn this subsection the effects of the irradiance noise on the biogechemistry are analysed. In Fig. 6 the dependence on (tau) of both the CV [panel (a)] and the mean value concentration [panel (b)] of detritus, labile dissolved organic carbon (L-DOC), semi-labile dissolved organic carbon (SL-DOC) and gross primary production (GPP) are shown. All these biogeochemical properties are correlated with carbon cycling. Gross primary production is related to the amount of carbon entering in the ecosystem, and is related to the maximum energy available in the ecosystem progressively dissipated in the trophic web. Gross primary production is directly affected by light fluctuation and its CV shape is very similar to that of the irradiance, Fig. 2. We selected also detritus and DOC because they are important indicators for the carbon cycling dynamics and are related to the cycling of chemicals like heavy metals62. The different curves, related to different values of D, approach the same (vanishing) value for large (tau). As previously discussed for the CV [Fig. 6(a)] of biomass concentrations, this circumstance is due to the fact that, in this case, the system dynamics tends to the deterministic case, characterized by a unique possible realization implying a vanishing standard deviation. For high correlation times thus the system is insensitive to the noise intensity. On the contrary, for small values of (tau), different values of D lead to significant differences of the variance. In particular, detritus, L-DOC and SL-DOC exhibit a clear non-monotonic behaviour whose maximum value depends on the combined values of D-(tau). Only the GPP presents a decreasing monotonic behaviour.The dependence of the mean value concentration on (tau), instead, is qualitatively the same for all the four parameters. Also in this case we can note a diversification with respect to D occurring at small (tau) and a (deterministic) constant value arising for low (high) values of D ((tau)).These results manifest that not only the population dynamics, but also all the biogeochemical processes are profoundly affected by the presence of stochastic environmental variables. The values and the behaviour of the examined quantities are indeed determined by the intimate interplay between the intensity and the time correlation of the noise fluctuations.Figure 6(a) Coefficient of variation ((CV=sigma / mu)) and (b) mean value concentration ((mu)) of detritus, labile dissolved organic carbon (L-DOC), semi-labile dissolved organic carbon (SL-DOC) and gross primary production (GPP) resulting from numerical integration of model equations plotted versus the considered values of (tau); the different curves are related to different values of the correlation time D.Full size image More

Armbrust, E. V. The life of diatoms in the world’s oceans. Nature 459, 185–192 (2009).CAS 
Article 

Google Scholar 
Mann, D. G. The species concept in diatoms. Phycologia 38, 437–495 (1999).Article 

Google Scholar 
Smetacek, V. Diatoms and the ocean carbon cycle. Protist 150, 25–32 (1999).CAS 
Article 

Google Scholar 
Rynearson, T. A. et al. Major contribution of diatom resting spores to vertical flux in the sub-polar North Atlantic. Deep. Res. Part I Oceanogr. Res. Pap. 82, 60–71 (2013).CAS 
Article 

Google Scholar 
Allen, J. T. et al. Diatom carbon export enhanced by silicate upwelling in the northeast Atlantic. Nature 437, 728–732 (2005).CAS 
Article 

Google Scholar 
Boyd, P. W., Strzepek, R., Fu, F. & Hutchins, D. A. Environmental control of open-ocean phytoplankton groups: Now and in the future. Limnol. Oceanogr. 55, 1353–1376 (2010).CAS 
Article 

Google Scholar 
Hátún, H., Somavilla, R., Rey, F., Johnson, C. & Mathis, M. The subpolar gyre regulates silicate concentrations in the North Atlantic. Sci. Rep. 1–9 https://doi.org/10.1038/s41598-017-14837-4 (2017).Bopp, L. Response of diatoms distribution to global warming and potential implications: a global model study. Geophys. Res. Lett. 32, 2–5 (2005).Article 
CAS 

Google Scholar 
Warner, A. J. & Hays, G. C. Sampling by the Continuous Plankton Recorder survey. Prog. Oceanogr. 6611, 237–256 (1994).Article 

Google Scholar 
Edwards, M., Beaugrand, G., Reid, P. C., Rowden, A. A. & Jones, M. B. Ocean climate anomalies and the ecology of the North Sea. Mar. Ecol. Prog. Ser. 239, 1–10 (2002).Article 

Google Scholar 
Allen, S. et al. Interannual stability of phytoplankton community composition in the North-East atlantic. Mar. Ecol. Prog. Ser. 655, 43–57 (2020).Article 

Google Scholar 
Edwards, M. & Richardson, A. J. Impact of climate change on marine pelagic phenology and trophic mismatch. Nature 430, 881–884 (2004).CAS 
Article 

Google Scholar 
Wihsgott, J. U. et al. Observations of vertical mixing in autumn and its effect on the autumn phytoplankton bloom. Prog. Oceanogr. 177, 1157–1165 (2019).Article 

Google Scholar 
Kamykowski, D. & Zentara, S. J. Predicting plant nutrient concentrations from temperature and sigma-T in the upper kilometer of the world ocean. Deep. Res. Part A-Oceanogr. Res. Pap. 33, 89–105 (1986).CAS 
Article 

Google Scholar 
Kamykowski, D. & Zentara, S. J. Changes in world ocean nitrate availability through the 20th century. Deep. Res. Part I-Oceanogr. Res. Pap. 52, 1719–1744 (2005).Article 

Google Scholar 
Behrenfeld, M. J. et al. Climate-driven trends in contemporary ocean productivity. Nature 444, 752–755 (2006).CAS 
Article 

Google Scholar 
López-Urrutia, A., San Martin, E., Harris, R. P. & Irigoien, X. Scaling the metabolic balance of the oceans. Proc. Natl. Acad. Sci. USA 103, 8739–8744 (2006).Article 
CAS 

Google Scholar 
Richardson, A. J. & Schoeman, D. S. Climate impact on plankton ecosystems in the Northeast Atlantic. Science 305, 1609–1612 (2004).CAS 
Article 

Google Scholar 
Edwards, M., Johns, D., Leterme, S., Svendsen, E. & Richardson, A. Regional climate change and harmful algal blooms in the northeast Atlantic. Limnol. Oceanogr. 51, 820–829 (2006).Article 

Google Scholar 
Hinder, S. L. et al. Changes in marine dinoflagellate and diatom abundance under climate change. Nat. Clim. Chang. 2, 271–275 (2012).Article 

Google Scholar 
Batten, S. et al. CPR sampling: the technical background, materials and methods, consistency and comparability. Prog. Oceanogr. 58, 193–215 (2003).Article 

Google Scholar 
Reid, P. C. et al. The Continuous Plankton Recorder: concepts and history, from plankton indicator to undulating recorders. Prog. Oceanogr. 58, 117–173 (2003).Article 

Google Scholar 
Edwards, M., Beaugrand, G., Hays, G. C., Koslow, J. A. & Richardson, A. J. Multi-decadal oceanic ecological datasets and their application in marine policy and management. Trends Ecol. Evol. 25, 602–610 (2010).Article 

Google Scholar 
Edwards, M. et al. North Atlantic warming over six decades drives decreases in krill abundance with no associated range shift. Commun. Biol. 4, 1–10 (2021).Article 

Google Scholar 
Richardson, A. J. et al. Using continuous plankton recorder data. Prog. Oceanogr. 68, 27–74 (2006).Article 

Google Scholar 
Hélaouët, P., Beaugrand, G. & Reygondeau, G. Reliability of spatial and temporal patterns of C. finmarchicus inferred from the CPR survey. J. Mar. Syst. 153, 18–24 (2016).Article 

Google Scholar 
Owens, N. J. P. et al. All plankton sampling systems underestimate abundance: response to “Continuous plankton recorder underestimates zooplankton abundance” by J.W. Dippner and M. Krause. J. Mar. Syst. 128, 240–242 (2013).Article 

Google Scholar 
Jonas, T. D., Walne, A., Beaugrand, G., Gregory, L. & Hays, G. C. The volume of water filtered by a Continuous Plankton Recorder sample: the effect of ship speed. J. Plankton Res. 26, 1499–1506 (2004).Article 

Google Scholar 
O’Reilly, C. H., Zanna, L. & Woollings, T. Assessing external and internal sources of Atlantic multidecadal variability using models, proxy data, and early instrumental indices. J. Clim 32, 7727–7745 (2019).Article 

Google Scholar 
Qin, M., Dai, A. & Hua, W. Quantifying contributions of internal variability and external forcing to atlantic multidecadal variability since 1870. Geophys. Res. Lett. 47, 1–11 (2020).
Google Scholar 
Mann, M. E., Steinman, B. A. & Miller, S. K. Absence of internal multidecadal and interdecadal oscillations in climate model simulations. Nat. Commun. 11, 1–9 (2020).Article 
CAS 

Google Scholar 
Enfield, D. B., Mestas-Nuñez, A. M. & Trimble, P. J. The Atlantic Multidecadal Oscillation and its relation to rainfall and river flows in the continental U.S. Geophys. Res. Lett. 28, 2077–2080 (2001).Article 

Google Scholar 
Gray, S. T., Graumlich, L. J., Betancourt, J. L. & Pederson, G. T. A tree-ring based reconstruction of the Atlantic Multidecadal Oscillation since 1567 A.D. Geophys. Res. Lett. 31, 2–5 (2004).Article 

Google Scholar 
Edwards, M., Beaugrand, G., Helaouët, P., Alheit, J. & Coombs, S. Marine ecosystem response to the Atlantic Multidecadal Oscillation. PLoS One 8, e57212 (2013).CAS 
Article 

Google Scholar 
Jones, P. D., New, M., Parker, D. E. & Martin, S. & Rigor, I. G. Surface air temperature and its changes over the past 150 years. Rev. Geophys. 37, 173–199 (1999).Article 

Google Scholar 
Beaugrand, G. et al. Reorganization of North Atlantic marine copepod biodiversity and climate. Science296, 1692–1694 (2002).CAS 
Article 

Google Scholar 
Alvain, S., Moulin, C., Dandonneau, Y. & Bréon, F. M. Remote sensing of phytoplankton groups in case 1 waters from global SeaWiFS imagery. Deep Sea Res. Part I Oceanogr. Res. Pap. 52, 1989–2004 (2005).Article 

Google Scholar 
Kalnay, E. et al. The NCEP/NCAR 40-year reanalysis project. Bull. Am. Meteorol. Soc. 77, 437–471 (1996).Article 

Google Scholar 
Huang, B. et al. Extended reconstructed sea surface temperature, version 5 (ERSSTv5): upgrades, validations, and intercomparisons. J. Clim. 30, 8179–8205 (2017).Article 

Google Scholar 
Lam, N. S. N. Spatial interpolation methods: a review. Am. Cartogr. 10, 129–150 (1983).Article 

Google Scholar 
Beaugrand, G., McQuatters-Gollop, A., Edwards, M. & Goberville, E. Long-term responses of North Atlantic calcifying plankton to climate change. Nat. Clim. Chang. 3, 263–267 (2012).Article 
CAS 

Google Scholar 
Beaugrand, G. et al. Prediction of unprecedented biological shifts in the global ocean. Nat. Clim. Chang. 9, 237–243 (2019).Article 

Google Scholar  More

Chapron, G. et al. Recovery of large carnivores in Europe’s modern human-dominated landscapes. Science 346, 1517–1519 (2014).ADS 
CAS 
PubMed 
Article 

Google Scholar 
Lute, M. L., Carter, N. H., López-Bao, J. V. & Linnell, J. D. C. Conservation professionals’ views on governing for coexistence with large carnivores. Biol. Cons. 248, 108668 (2020).Article 

Google Scholar 
Gantchoff, M. G. & Belant, J. L. Regional connectivity for recolonizing American black bears (Ursus americanus) in southcentral USA. Biol. Cons. 214, 66–75 (2017).Article 

Google Scholar 
Ripple, W. J. et al. Status and ecological effects of the world’s largest carnivores. Science 343, 25 (2014).Article 
CAS 

Google Scholar 
Kays, R. et al. Does hunting or hiking affect wildlife communities in protected areas?. J. Appl. Ecol. 54, 242–252 (2017).Article 

Google Scholar 
Schipper, J. et al. The status of the world’s land and marine mammals: diversity, threat, and knowledge. Science 322, 225–230 (2008).ADS 
CAS 
PubMed 
Article 

Google Scholar 
Smith, J. A., Wang, Y. & Wilmers, C. C. Top carnivores increase their kill rates on prey as a response to human-induced fear. Proc. R. Soc. B Biol. Sci. 282, 20142711 (2015).Article 

Google Scholar 
Stillfried, M., Belant, J. L., Svoboda, N. J., Beyer, D. E. & Kramer-Schadt, S. When top predators become prey: Black bears alter movement behaviour in response to hunting pressure. Behav. Proc. 120, 30–39 (2015).Article 

Google Scholar 
Støen, O.-G. et al. Physiological evidence for a human-induced landscape of fear in brown bears (Ursus arctos). Physiol. Behav. 152, 244–248 (2015).PubMed 
Article 
CAS 

Google Scholar 
Evans, M. J., Rittenhouse, T. A. G., Hawley, J. E. & Rego, P. W. Black bear recolonization patterns in a human-dominated landscape vary based on housing: New insights from spatially explicit density models. Landsc. Urban Plan. 162, 13–24 (2017).Article 

Google Scholar 
LaRue, M. A. et al. Cougars are recolonizing the midwest: Analysis of cougar confirmations during 1990–2008. J. Wildl. Manag. 76, 1364–1369 (2012).Article 

Google Scholar 
Cove, M. V., Fergus, C., Lacher, I., Akre, T. & McShea, W. J. Projecting mammal distributions in response to future alternative landscapes in a rapidly transitioning region. Remote Sens. 11, 2482 (2019).ADS 
Article 

Google Scholar 
Frid, A. & Dill, L. Human-caused disturbance stimuli as a form of predation risk. Conserv. Ecol. 6, 25 (2002).
Google Scholar 
Clinchy, M. et al. Fear of the human “super predator” far exceeds the fear of large carnivores in a model mesocarnivore. Behav. Ecol. 27, 1826–1832 (2016).
Google Scholar 
Suraci, J. P., Clinchy, M., Zanette, L. Y. & Wilmers, C. C. Fear of humans as apex predators has landscape-scale impacts from mountain lions to mice. Ecol. Lett. 22, 1578–1586 (2019).PubMed 
Article 

Google Scholar 
Gaynor, K. M., Hojnowski, C. E., Carter, N. H. & Brashares, J. S. The influence of human disturbance on wildlife nocturnality. Science 360, 1232–1235 (2018).ADS 
CAS 
PubMed 
Article 

Google Scholar 
Smith, J. A., Thomas, A. C., Levi, T., Wang, Y. & Wilmers, C. C. Human activity reduces niche partitioning among three widespread mesocarnivores. Oikos 127, 890–901 (2018).Article 

Google Scholar 
Tucker, M. A. et al. Moving in the Anthropocene: Global reductions in terrestrial mammalian movements. Science 359, 466–469 (2018).ADS 
CAS 
PubMed 
Article 

Google Scholar 
Carter, N. H., Brown, D. G., Etter, D. R. & Visser, L. G. American black bear habitat selection in northern Lower Peninsula, Michigan, USA, using discrete-choice modeling. Ursus 21, 57–71 (2010).Article 

Google Scholar 
Naidoo, R. & Burton, A. C. Relative effects of recreational activities on a temperate terrestrial wildlife assemblage. Conserv. Sci. Pract. 2, e271 (2020).
Google Scholar 
Geffroy, B., Samia, D. S. M., Bessa, E. & Blumstein, D. T. How nature-based tourism might increase prey vulnerability to predators. Trends Ecol. Evol. 30, 755–765 (2015).PubMed 
Article 

Google Scholar 
Geffroy, B. et al. Evolutionary dynamics in the Anthropocene: Life history and intensity of human contact shape antipredator responses. PLoS Biol. 18, e3000818 (2020).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Beeco, J. A., Hallo, J. C. & Brownlee, M. T. J. GPS visitor tracking and recreation suitability mapping: tools for understanding and managing visitor use. Landsc. Urban Plan. 127, 136–145 (2014).Article 

Google Scholar 
Thorsen, N. H. et al. Smartphone app reveals that lynx avoid human recreationists on local scale, but not home range scale. Sci. Rep. 12, 1–13 (2022).Article 
CAS 

Google Scholar 
Evans, M. J., Hawley, J. E., Rego, P. W. & Rittenhouse, T. A. G. Hourly movement decisions indicate how a large carnivore inhabits developed landscapes. Oecologia 190, 11–23 (2019).ADS 
PubMed 
Article 

Google Scholar 
Carlos, A. W. D., Bright, A. D., Teel, T. L. & Vaske, J. J. Human-black bear conflict in urban areas: an integrated approach to management response. Hum. Dimens. Wildl. 14, 174–184 (2009).Article 

Google Scholar 
Johnson, H. E. et al. Human development and climate affect hibernation in a large carnivore with implications for human–carnivore conflicts. J. Appl. Ecol. 55, 663–672 (2018).Article 

Google Scholar 
Gould, N. P., Powell, R., Olfenbuttel, C. & DePerno, C. S. Growth and reproduction by young urban and rural black bears. J. Mammal. 102, 1165–1173 (2021).Article 

Google Scholar 
Ditmer, M. A., Noyce, K. V., Fieberg, J. R. & Garshelis, D. L. Delineating the ecological and geographic edge of an opportunist: The American black bear exploiting an agricultural landscape. Ecol. Model. 387, 205–219 (2018).Article 

Google Scholar 
McFadden-Hiller, J. E. Jr. & Belant, J. L. Spatial distribution of black bear incident reports in michigan. PLoS One 11, e0154474 (2016).PubMed 
PubMed Central 
Article 
CAS 

Google Scholar 
Ladle, A., Steenweg, R., Shepherd, B. & Boyce, M. S. The role of human outdoor recreation in shaping patterns of grizzly bear-black bear co-occurrence. PLoS One 13, e0191730 (2018).PubMed 
PubMed Central 
Article 

Google Scholar 
Wilbur, R. C., Lischka, S. A., Young, J. R. & Johnson, H. E. Experience, attitudes, and demographic factors influence the probability of reporting human–black bear interactions. Wildl. Soc. Bull. 42, 22–31 (2018).Article 

Google Scholar 
Lustig, E. J., Lyda, S. B., Leslie, D. M., Luttbeg, B. & Fairbanks, W. S. Resource selection by recolonizing American Black Bears. J. Wildl. Manage. 85, 531–542 (2021).Article 

Google Scholar 
Sun, C. C., Fuller, A. K., Hare, M. P. & Hurst, J. E. Evaluating population expansion of black bears using spatial capture-recapture. J. Wildl. Manage. 81, 814–823 (2017).Article 

Google Scholar 
Kautz, T. M. et al. Large carnivore response to human road use suggests a landscape of coexistence. Glob. Ecol. Conserv. 30, e01772 (2021).Article 

Google Scholar 
Michigan Department of Natural Resources (MIDNR) (2021).Blount, J. D., Chynoweth, M. W., Green, A. M. & Şekercioğlu, Ç. H. Review: COVID-19 highlights the importance of camera traps for wildlife conservation research and management. Biol. Cons. 256, 108984 (2021).Article 

Google Scholar 
Weather Atlas. https://www.weather-atlas.com/enEvans, J. S. Spatial Analysis and Modelling Utilities. Package ‘spatialEco’. https://cran.r-project.org/web/packages/spatialEco/spatialEco.pdf (2021).Díaz-Ruiz, F., Caro, J., Delibes-Mateos, M., Arroyo, B. & Ferreras, P. Drivers of red fox (Vulpes vulpes) daily activity: prey availability, human disturbance or habitat structure?. J. Zool. 298, 128–138 (2016).Article 

Google Scholar 
Moore, J. F. et al. Comparison of species richness and detection between line transects, ground camera traps, and arboreal camera traps. Anim. Conserv. 23, 561–572 (2020).Article 

Google Scholar 
Parsons, A. W. et al. Urbanization focuses carnivore activity in remaining natural habitats, increasing species interactions. J. Appl. Ecol. 56, 1894–1904 (2019).Article 

Google Scholar 
Allen, M. L., Sibarani, M. C., Utoyo, L. & Krofel, M. Terrestrial mammal community richness and temporal overlap between tigers and other carnivores in Bukit Barisan Selatan National Park, Sumatra. Anim. Biodiv. Conserv. 43(1), 97–107 (2020).Article 

Google Scholar 
Tian, C. et al. Temporal niche patterns of large mammals in Wanglang National Nature Reserve, China. Glob. Ecol. Conserv. 22, e01015 (2020).Article 

Google Scholar 
Meredith, M. & Ridout, M. Estimates of coefficient of overlapping for animal activity patterns. Package ‘overlap’. https://cran.r-project.org/web/packages/overlap/overlap.pdf (2020).RStudio Team. RStudio: Integrated Development for R. RStudio, PBC, Boston, MA. http://www.rstudio.com/ (2021).Ridout, M. S. & Linkie, M. Estimating overlap of daily activity patterns from camera trap data. JABES 14, 322–337 (2009).MathSciNet 
MATH 
Article 

Google Scholar 
Lashley, M. A. et al. Estimating wildlife activity curves: comparison of methods and sample size. Sci. Rep. 8, 1–11 (2018).CAS 
Article 

Google Scholar 
Rowcliffe, M. Animal Activity Statistics. Package ‘activity’. https://cran.r-project.org/web/packages/activity/activity.pdf (2021).MacKenzie, D. I., Nichols, J. D., Hines, J. E., Knutson, M. G. & Franklin, A. B. Estimating site occupancy, colonization, and local extinction when a species is detected imperfectly. Ecology 84, 2200–2207 (2003).Article 

Google Scholar 
Wei, T., & Simko, V. Visualization of a Correlation Matrix. Package ‘corrplot’. https://cran.r-project.org/web/packages/corrplot/corrplot.pdf (2017).Norton, D. C. et al. Female American black bears do not alter space use or movements to reduce infanticide risk. PLoS One 13, e0203651 (2018).PubMed 
PubMed Central 
Article 
CAS 

Google Scholar 
Ditmer, M. A. et al. Behavioral and physiological responses of American black bears to landscape features within an agricultural region. Ecosphere 6, 1–21 (2015).Article 

Google Scholar 
Clark, D. et al. Using machine learning methods to predict the movement trajectories of the Louisiana black bear. SMU Data Sci. Rev. 5, 25 (2021).
Google Scholar  More

Growth conditions and plant materialsTwo experiments were conducted concurrently (sites A and B) in the same screen house in 2019 between the months of May and July at the Landmark University Greenhouse and Hydroponic Technology Center, a section of the Teaching and Research Farm of the University in Omu-Aran, Kwara State Nigeria. Experiment at site B was conducted simultaneously as A so as to validate the results of experiment A. Landmark University lies within Latitude 8° 7′ 26.21388″ and 5° 5′ 0.1788″. Both experiments (A & B) involved tomato (Solanum lycopersicum L. variety cherry) and cucumber (Cucumis sativus L. variety marketer) crops. For each crop, seeds were sown into a separate seed tray filled with coco peat (Coco peat, SRIMATHI EXPORT, INDIA). Cocopeat is the mesocarp tissue or husk after the grinding of coconut fruit. It has a lightweight and high water and nutrient holding capacities, it has an acceptable pH, electrical conductivity, and other chemical attributes27. Rice husk is the by-product of rice after milling. The rice husk used was collected from the rice processing mill of Landmark University. Rice husk is a highly porous and light weighted material with a very high specific area28.Two sets of seed trays (one for organic and another for inorganic fertilizers) were used each for tomato and cucumber crops in the nursery. Both were raised in the nursery for two weeks before transplanting. Black grow bags (30 × 17 cm) filled with a coco peat/rice husk (1:4 ratio by volume) mixture with a weight of about 10 kg were arranged in a screen house. Both the nursery and establishment of crop proper take place in a screen house. The screen house has a galvanized iron as the frame, a UV covering on top, side net for screening insect pests the floor fairly covered with granite. Temperature and relative humidity within the screen house during the period of the experiment was monitored using a Thermograph and a Barograph, and they were at an average of 31 °C and 75%, respectively.The grow bags were randomly placed in the screen house for the unbiased application of amendments. For both tomato and cucumber crops, the treatment comprised of six (6) levels of liquid organic fertilizer (5, 15, 25, 35, 45, 55 mL), in-organic fertilizer, and a control (ordinary borehole water). Levels of organic fertilizers were selected based on the recommendation of 20 mL of liquid organic fertilizer by29. The eight (8) treatments both for tomato and cucumber were arranged in a Completely Randomized Design replicated three times. One healthy plant was maintained per grow bag and four grow bags represent a treatment and there were 32 plants per block each for tomato and cucumber. For both crops, the experiment lasted for 90 days.Organic and in-organic nutrient solutionsThe liquid organic fertilizer used was obtained from the biomass of Mexican sunflower (Tithonia diversifolia). Fresh biomass (mainly leaves and stems) of the plant was collected from the Teaching and Research Farms of Landmark University, Nigeria. After rinsing, they were cut with a sterile knife into pieces of ≤ 1 cm size. A sample was taken for initial physicochemical analyses by grinding in a sterile mortal, diluted with sterile water and analyzed. The biomass was then soaked in sterile water inside a clean container, and allowed to ferment spontaneously for a period of 14 days. During the fermentation, samples were taken every 4 days for microbial analyses of the major players during the fermentation. At the end of fermentation, the mixture was separated using a sieve of mesh size ≤ 2 mm. The liquid portion was then refrigerated prior to the planting regime while another sample was taken to ascertain the physicochemical and microbial qualities of the produced liquid fertilizer. The chemical analysis is presented in Table 4. For inorganic fertilizer, Water soluble fertilizers employed in hydroponics were used (Hydroponics fertilizer, Anmol chemicals, India); calcium nitrate 650 mg L−1, potassium nitrate 450 mg L−1, magnesium 400 mg L−1, chelate 20 mg L−1, mono-ammonium phosphate 400 mg L−1. The electrical conductivity (EC) of the solution was 1.9 dS m-1.Irrigation and fertigationThe tomato and cucumber plants were fertigated morning and evening daily for one hour on each occasion according to the treatments. Preparation of the nutrient solution was with borehole water and was supplied to plants by an online pressure drip irrigation system set at 2.0 L h-1 using an arrowhead on each tomato and cucumber plant. Different tanks (250 L) were installed according to the various treatments making a total of 8 tanks. The organic fertilizer was diluted according to the various treatments equivalent to 1.25, 3.75, 6.25, 8.75, 11.25, and 13.75 L per 250 L of water respectively for 5, 15, 25, 35, 45, and 55 mL treatments. The nutrient solutions were refilled when the consumption is less than 20% of the initial volume (250 L) in the tank. One day per week, crops were irrigated with ordinary water to wash out pipes and prevent deposits of salts. The same concentration of nutrient was used from transplanting to the termination of the study for both tomato and cucumber crops, however, at the flowering of the crops, the volume of fertigation was increased to 3.0 L h-1 to be able to cope with the size of the plants.Trellising, pest and diseases controlFor both tomato and cucumber crops, plant vines were supported by twisting them around a wire that is- attached to the roof of the screen house and 2 m from the ground. Lateral outgrowths were cut off every week to ensure a sturdy single stem. Pests and diseases were scouted every day. Whiteflies, aphids, and other insects were controlled with orizon (Producer, location of producer) (active ingredient, acetamiprid, and abamectin) using 0.133% v/v. Fungi were controlled using ridomil gold (Producer, Location of producer) at 2% w/v.Determination of growth and yield of tomato and cucumberThree tomato and cucumber plants were randomly selected for each treatment for the determination of growth parameters (plant height, leaf area, number of leaves per plant, and stem diameter) at mid the flowering stage of tomato and cucumber plants.The leaf area of tomato was calculated using the model (A = KL2) developed by Lyon30, where L = Length of tomato leaf, K = constant which is 0.1551, and A = leaf area of tomato. Similarly, the leaf area of cucumber was calculated using A = 0.88LW – 4.27, where L = cucumber leaf length and W = cucumber leaf width, A = leaf area of cucumber31.Tomato fruits were ready for harvest from 65 days after transplanting, harvestings were done twice every week (Mondays and Fridays) for up to 85 days after transplanting. Similarly, harvesting of cucumber fruits started 35 days after transplanting and harvestings were also done twice a week (Mondays and Fridays), harvesting was carried out till 60 days after transplanting. Tomato and cucumber fruit yields were counted and weighed at each harvest.Analysis of tomato and cucumber leaves and fruitsAt the 50% flowering stage of tomato and cucumber plants, ten leaf samples were collected from each treatment. The leaf samples were oven-dried at 75 °C for 24 h and thereafter grounded. The grounded samples were later analyzed for nitrogen (N), phosphorous (P), potassium (K), calcium (Ca), and magnesium (Mg) content using the method of described by32. At harvest, four matured tomato and cucumber fruits of uniform size were selected per treatment, and their nutrient compositions were determined using the method of33.Statistical analysisAll data collected on the growth, yield, leaf, and fruit nutrient contents of tomato and cucumber were subjected to analysis of variance (ANOVA). The SPSS V 21.0 (New York, USA) software was used to perform ANOVA and Duncan’s multiple range test (DMRT) was used to compare means at a 5% probability level.
Ethical approvalI confirm that all the research meets ethical guidelines and adheres to the legal requirements of the study country.Compliance with international, national and/or institutional guidelinesExperimental research (either cultivated or wild), comply with relevant institutional, national, and international guidelines and legislation. Experimental studies were carried out in accordance with relevant institutional, national or international guidelines or regulation. More

Gifford-Gonzalez, D. Bones are not enough: Analogues, knowledge, and interpretive strategies in zooarchaeology. J. Anthropol. Archaeol. 10, 215–254 (1991).Article 

Google Scholar 
Pobiner, B. L. The zooarchaeology and paleoecology of early hominin scavenging. Evol. Anthropol. 2, 68–82 (2020).Article 

Google Scholar 
Rodriguez, A. et al. Right or left? Determining the hand holding the tool from use traces. J. Archaeol. Sci. Rep. 31, 102316 (2020).
Google Scholar 
Feix, T., Kivell, T. L., Pouydebat, E. & Dollar, A. M. Estimating thumb-index finger precision grip and manipulation potential in extant and fossil primates. J. R. Soc. Interface. https://doi.org/10.1098/rsif.2015.0176 (2015).Article 
PubMed 
PubMed Central 

Google Scholar 
Bardo, A. et al. The implications of thumb movements for Neanderthal and modern human manipulation. Sci. Rep. 10, 19323 (2020).ADS 
CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Stout, D., Semaw, S., Rogers, M. J. & Cauche, D. Technological variation in the earliest Oldowan from Gona, Afar, Ethiopia. J. Hum. Evol. 58, 474–491 (2010).PubMed 
Article 

Google Scholar 
Tennie, C., Premo, L. S., Braun, D. R. & McPherron, S. P. Resetting the null hypothesis: Early stone tools and cultural transmission. Curr. Anthrop. 58, 652–672 (2017).Article 

Google Scholar 
Tennie, C. The zone of latent solution (ZLS) account remains the most parsimonious explanation for early stone tools. Curr. Anthrop. 60, 331–332 (2019).
Google Scholar 
Tennie, C., Braun, D. R., Premo, L. S. & McPherron, S. P. The Island Test for Cumulative Culture in Paleolithic Cultures. In The Nature of Culture. Series: Vertebrate Paleobiology and Paleoanthropology (eds Haidle, M. N. et al.) (Springer, 2016).
Google Scholar 
Perreault, C. The Quality of the Archaeological Record (University of Chicago Press, 2019).Book 

Google Scholar 
Proffitt, T. et al. Wild monkeys flake stone tools. Nature 539, 85–88 (2016).ADS 
CAS 
PubMed 
Article 

Google Scholar 
Carvalho, S., Cunha, E., Sousa, C. & Matsuzawa, T. Chaînes opératoires and resource-exploitation strategies in chimpanzee (Pan troglodytes) nut cracking. J. Hum. Evol. 55, 148–163 (2008).PubMed 
Article 

Google Scholar 
Westergaard, G. C. & Suomi, S. J. The stone tools of capuchins (Cebus apella). Int. J. Primatol. 16, 1017–1024 (1995).Article 

Google Scholar 
De la Torre, I. & Mora, R. Technological Strategies in the Lower Pleistocene at Olduvai Beds I and II (Service de Prehistoire, Universite de Liege, 2005).
Google Scholar 
M. D. O. M. Í. Dominguez-Rodrigo, 3.3-Million-Year-Old Stone Tools and Butchery Traces? More Evidence Needed. PaleoAnthropology. 9 (2016).Harmand, S. et al. 3.3-million-year-old stone tools from Lomekwi 3, West Turkana, Kenya. Nature 521, 310–315 (2015).ADS 
CAS 
PubMed 
Article 

Google Scholar 
Andrefsky, W. Lithics: Macroscopic Approaches to Analysis (Cambridge University Press, 2005).Book 

Google Scholar 
Malaivijitnond, S. et al. Stone-tool usage by Thai long-tailed macaques (Macaca fascicularis). Am. J. Primatol. 69, 227–233 (2007).PubMed 
Article 

Google Scholar 
Luncz, L. V. et al. Resource depletion through primate stone technology. eLife 6, e23647 (2017).PubMed 
PubMed Central 
Article 

Google Scholar 
Leca, J. B., Gunst, N. & Huffman, M. Complexity in object manipulation by Japanese macaques (Macaca fuscata): A cross-sectional analysis of manual coordination in stone handling patterns. J. Comp. Psychol. 125, 61 (2011).PubMed 
Article 

Google Scholar 
Toth, N., Schick, K. D., Savage-Rumbaugh, E. S., Sevcik, R. A. & Rumbaugh, D. M. Pan the tool-maker: Investigations into the stone tool-making and tool-using capabilities of a bonobo (Pan paniscus). J. Archaeol. Sci. 20, 81–91 (1993).Article 

Google Scholar 
Wright, R. V. S. Imitative learning of a flaked stone technology-The case of an orangutan. Mankind 8, 296–306 (2009).
Google Scholar 
Bandini, E. et al. Naïve, unenculturated chimpanzees fail to make and use flaked stone tools. Open Res. Eur. https://doi.org/10.12688/openreseurope.13186.2 (2021).Article 
PubMed 
PubMed Central 

Google Scholar 
J. Henrich, C. Tennie, in Chimpanzees and Human Evolution, M. Muller, R. Wrangham, D. Pilbeam, Eds. (Harvard University Press, Cambridge, MA, (2017), 645–702.Breuer, T., Ndoundou-Hockemba, M. & Fishlock, V. First observation of tool use in wild gorillas. PLoS Biol. 3, 2041–2043 (2005).CAS 
Article 

Google Scholar 
Wittiger, L., Society, W. C., River, C. & Project, G. Tool use during display behavior in wild cross river gorillas. Am. J. Primat. 5, 1–5 (2007).
Google Scholar 
Kinani, J. F. & Zimmerman, D. Tool use for food acquisition in a wild mountain gorilla (Gorilla beringei beringei). Am. J. Primat. 77, 353–357 (2015).Article 

Google Scholar 
Grueter, C. C., Robbins, M. M., Ndagijimana, F. & Stoinski, T. S. Possible tool use in a mountain gorilla. Behav. Processes. 100, 160–162 (2013).PubMed 
Article 

Google Scholar 
Parker, S. T., Kerr, M., Markowitz, H. & Gould, J. A survey of tool use in zoo gorillas. In The Mentalities of Gorillas and Orangutans: Comparative Perspectives (eds Parker, S. T. et al.) (Cambridge University Press, 1999).Chapter 

Google Scholar 
Shumaker, R. W., Walkup, K. R. & Beck, B. B. Animal Tool Behavior: The Use and Manufacture of Tools by Animals (Johns Hopkins University Press, 2011).
Google Scholar 
Pouydebat, E., Berge, C., Gorce, P. & Coppens, Y. Use and manufacture of tools to extract food by captive Gorilla gorilla gorilla: Experimental approach. Folia Primat. 76, 180–183. https://doi.org/10.1159/000084381 (2005).Article 

Google Scholar 
Haslam, M. ‘Captivity bias’ in animal tool use and its implications for the evolution of hominin technology. PTRBAE 368, 20120421 (2013).PubMed 
PubMed Central 

Google Scholar 
Van Schaik, C. P., Deaner, R. O. & Merrill, M. Y. The conditions for tool use in primates: Implications for the evolution of material culture. J. Hum. Evol. 36, 719–741 (1999).PubMed 
Article 

Google Scholar 
Pouydebat, E., Gorce, P., Coppens, Y. & Bels, V. Biomechanical study of grasping according to the volume of the object: Human versus non-human primates. J. Biomech. 42, 266–272 (2009).PubMed 
Article 

Google Scholar 
Pouydebat, E., Laurin, M., Gorce, P. & Bels, V. Evolution of grasping among anthropoids. J. Evol. Bio. 21, 1732–1743 (2008).CAS 
Article 

Google Scholar 
Bardo, A., Cornette, R., Borel, A. & Pouydebat, E. Manual function and performance in humans, gorillas and orangutans during the same tool use task. Am. J. Phys. Anthropol. https://doi.org/10.1002/ajpa.2332 (2017).Article 
PubMed 

Google Scholar 
A. Bardo, A. Borel, H. Meunier, J. P. Guéry, E. Pouydebat, Manual abilities in great apes during a tool use task. Am. J. Phys. Anthropol. doi: 10.1002 (2016).W. C. McGrew, Why is ape tool use so confusing. Comparative socioecology: the behavioural ecology of humans and other mammals. 457–472 (1989).Cipolletta, C. et al. Termite feeding by Gorilla gorilla gorilla at Bai Hokou, Central African Republic. Int. J. Primatol. 28, 457–476 (2007).Article 

Google Scholar 
Salmi, R., Rahman, U. & Doran-Sheehy, D. M. Hand preference for a novel bimanual coordinated task during termite feeding in wild western gorillas (Gorilla gorilla gorilla). Int. J. Primatol. 37, 200–212 (2016).Article 

Google Scholar 
Masi, S. et al. The influence of seasonal frugivory on nutrient and energy intake in wild western gorillas. PLoS ONE 10, e0129254 (2015).PubMed 
PubMed Central 
Article 
CAS 

Google Scholar 
Redford, K. H. & Dorea, J. G. The nutritional value of invertebrates with emphasis on ants and termites as food for mammals. J. Zool. 203, 385–395 (1984).CAS 
Article 

Google Scholar 
McGrew, W. C. The ‘other faunivory’revisited: Insectivory in human and non-human primates and the evolution of human diet. J. Hum. Evol. 71, 4–11 (2014).PubMed 
Article 

Google Scholar 
Tennie, C., O’Malley, R. C. & Gilby, I. C. Why do chimpanzees hunt? Considering the benefits and costs of acquiring and consuming vertebrate versus invertebrate prey. J. Hum. Evol. 71, 38–45 (2014).PubMed 
Article 

Google Scholar 
McBrearty, S. Consider the humble termite: Termites as agents of post-depositional disturbance at African archaeological sites. J. Archaeol. Sci. 17, 111–143 (1990).Article 

Google Scholar 
Okwakol, M. J. N. Effects of Cubitermes testaceus (Williams) on some physical and chemical properties of soil in a grassland area of Uganda. Afr. J. Ecol. 25, 147–153 (1987).Article 

Google Scholar 
Altmann, J. Observational study of behavior: Sampling methods. Behavior 49, 227–267 (1974).CAS 
Article 

Google Scholar 
Robira, B. et al. Handedness in gestural and manipulative actions in male hunter-gatherer Aka pygmies from Central African Republic. Am. J. Phys. Anthropol. 166(481–491), 19 (2018).
Google Scholar 
Meguerditchian, A., Calcutt, S. E., Lonsdorf, E. V., Ross, S. R. & Hopkins, W. D. Brief communication: Captive gorillas are right-handed for bimanual feeding. Am. J. Phys. Anthropol. 141, 638–645 (2010).PubMed 
PubMed Central 

Google Scholar 
Dapena, J. E. S. Ú. S., William, J., Anderst, N. P. & Toth, The biomechanics of the arm swing in Oldowan stone flaking. In The Oldowan: Case Studies into the Earliest Stone Age (No. 1). Gosport (eds Toth, N. P. & Schick, K. D.) (Stone Age Institute Press, 2006).
Google Scholar 
Nowell, A. A. & Fletcher, A. W. The development of feeding behaviour in wild western lowland gorillas (Gorilla gorilla gorilla). Behaviour 145, 171–193 (2008).Article 

Google Scholar 
Pouydebat, E., Gorce, P., Coppens, Y. & Bels, V. Substrate optimization in nuts cracking by capuchin monkeys. Am. J. Primatol. 68, 1017–1024 (2006).CAS 
PubMed 
Article 

Google Scholar 
Boinski, S., Quatrone, R. P. & Swarttz, H. Substrate and tool use by brown capuchins in Suriname: Ecological contexts and cognitive bases. Am. Anthropol. 102, 741–761 (2000).Article 

Google Scholar 
Panger, M. A. Object-use in free-ranging white-faced capuchins (Cebus capucinus) in Costa Rica. Am. J. Phys. Anthropol. 106, 311–321 (1998).CAS 
PubMed 
Article 

Google Scholar 
Parker, S. T. & Gibson, K. R. Object manipulation, tool use and sensorimotor intelligence as feeding adaptations in cebus monkeys and great apes. J. Hum. Evol. 6, 623–641 (1977).Article 

Google Scholar 
Heldstab, S. A., Isler, K., Schuppli, C. & van Schaik, C. P. When ontogeny recapitulates phylogeny: Fixed neurodevelopmental sequence of manipulative skills among primates. Sci. Adv. 6, eabb4685 (2020).ADS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Clutton-Brock, T. H. Some aspects of intraspecific variation in feeding and ranging behaviour in primates. In Primate Ecology Studies of Feeding And Ranging Behavior in Lemurs, Monkeys and Apes (ed. Clutton-Brock, T. H.) (Academic Press, 1977).
Google Scholar 
Key, C. & Ross, C. Sex differences in energy expenditure in non-human primates. Proc. R. Soc. Lond. B. 266, 2479–2485 (1999).CAS 
Article 

Google Scholar 
Lockman, J. J. A perception–action perspective on tool use development. Child Dev. 71, 137–144 (2000).CAS 
PubMed 
Article 

Google Scholar 
Masi, S. et al. Unusual feeding behavior in wild great apes, a window to understand origins of self-medication in humans: Role of sociality and physiology on learning process. Physiol. Behav. 105, 337–349 (2012).CAS 
PubMed 
Article 

Google Scholar 
Falótico, T. & Ottoni, E. B. The manifold use of pounding stone tools by wild capuchin monkeys of Serra da Capivara National Park, Brazil. Behaviour 153, 421–442 (2016).Article 

Google Scholar 
Falótico, T. & Ottoni, E. B. Stone throwing as a sexual display in wild female bearded capuchin monkeys, Sapajus libidinosus. PLoS ONE 8, e79535 (2013).ADS 
PubMed 
PubMed Central 
Article 
CAS 

Google Scholar 
Mannu, M. & Ottoni, E. B. The enhanced tool-kit of two groups of wild bearded capuchin monkeys in the Caatinga: Tool making, associative use, and secondary tools. Am. J. Primatol. 71, 242–251 (2009).PubMed 
Article 

Google Scholar 
Gumert, M. D., Kluck, M. & Malaivijitnond, S. Thephysical characteristics and usage patterns of stoneaxe and pounding hammers used by long-tailedmacaques in the Andaman Sea region of Thailand. Am. J. Primatol. 71, 594–608. https://doi.org/10.1002/ajp.20694 (2009).Article 
PubMed 

Google Scholar 
Marzke, M. W. Precision grips, hand morphology, and tools. Am. J. Phys. Anthropol. 102, 91–110 (1997).CAS 
PubMed 
Article 

Google Scholar 
Matsuzawa, T. Chimpanzee Intelligence in Nature and in Captivity Isomorphism of Symbol Use and Tool Use (Cambridge University Press, 1996).Book 

Google Scholar 
Westergaard, G. C. & Suomi, S. J. A simple stone-tool technology in monkeys. J. Hum. Evol. 27, 399–404 (1994).Article 

Google Scholar 
Liu, Q. et al. Kinematics and energetics of nut-cracking in wild capuchin monkeys (Cebus libidinosus) in Piauí, Brazil. Am. J. Phys. Anthropol. 138, 210–220 (2009).CAS 
PubMed 
Article 

Google Scholar 
Günther, M. M. & Boesch, C. Energetic Cost of Nut-cracking Behaviour in Wild Chimpanzees. In Hands of Primates 109–129 (Springer, 1993).
Google Scholar 
Roach, N. T., Venkadesan, M., Rainbow, M. J. & Lieberman, D. E. Elastic energy storage in the shoulder and the evolution of high-speed throwing in Homo. Nature 498, 483–486 (2013).ADS 
CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Young, N. M., Capellini, T. D., Roach, N. T. & Alemseged, Z. Fossil hominin shoulders support an African ape-like last common ancestor of humans and chimpanzees. PNAS 112, 11829–11834 (2015).ADS 
CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Doran-Sheehy, D., Mongo, P., Lodwick, J. & Conklin-Brittain, N. L. Male and female western gorilla diet: Preferred foods, use of fallback resources, and implications for ape versus old world monkey foraging strategies. Am. J. Phys. Anthropol. 140, 727–738 (2009).CAS 
PubMed 
Article 

Google Scholar 
Breuer, T., Hockemba, M. B. N., Olejniczak, C., Parnell, R. J. & Stokes, E. J. Physical maturation, life-history classes and age estimates of free-ranging western gorillas – Insights from Mbeli Bai, Republic of Congo. Am. J. Primatol. 71, 106–119 (2009).PubMed 
Article 

Google Scholar 
Hopkins, W. D. et al. The use of bouts and frequencies in the evaluation of hand preferences for a coordinated bimanual task in chimpanzees (Pan troglodytes): An empirical study comparing two different indices of laterality. J. Comp. Psychol. 115, 294–299 (2001).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Byrne, R. W. & Byrne, J. M. Manual dexterity in the gorilla: bimanual and digit role differentiation in a natural task. Anim. Cogn. 4, 347–361 (2001).CAS 
PubMed 
Article 

Google Scholar  More

Pogoreutz C, Voolstra CR, Rädecker N, Weis V, Cardenas A, Raina J-B. The coral holobiont highlights the dependence of cnidarian animal hosts on their associated microbes. In Bosch TCG, Hadfield MG, editors. Cellular Dialogues in the Holobiont. CRC Press; 2020. pp. 91–118. https://doi.org/10.1201/9780429277375-7Rohwer F, Seguritan V, Azam F, Knowlton N. Diversity and distribution of coral-associated bacteria. Mar Ecol Prog Ser. 2002;243:1–10.
Google Scholar 
LaJeunesse TC, Parkinson JE, Gabrielson PW, Jeong HJ, Reimer JD, Voolstra CR, et al. Systematic revision of Symbiodiniaceae highlights the antiquity and diversity of coral endosymbionts. Curr Biol. 2018;28:2570–80.e6CAS 
PubMed 

Google Scholar 
Muscatine L, Porter JW. Reef corals: mutualistic symbioses adapted to nutrient-poor environments. Bioscience 1977;27:454–60.
Google Scholar 
Christian R, Voolstra DJ, Suggett RS, Peixoto JE, Parkinson KM, Quigley CB, et al. Extending the natural adaptive capacity of coral holobionts. Nature Reviews Earth & Environment. 2021;2:747–762. https://doi.org/10.1038/s43017-021-00214-3Article 

Google Scholar 
Bourne DG, Morrow KM, Webster NS. Insights into the coral microbiome: underpinning the health and resilience of reef ecosystems. Annu Rev Microbiol. 2016;70:317–40.CAS 
PubMed 

Google Scholar 
Rädecker N, Pogoreutz C, Voolstra CR, Wiedenmann J, Wild C. Nitrogen cycling in corals: the key to understanding holobiont functioning? Trends Microbiol. 2015;23:490–7.PubMed 

Google Scholar 
Matthews JL, Raina JB, Kahlke T, Seymour JR, van Oppen MJ, Suggett DJ. Symbiodiniaceae‐bacteria interactions: rethinking metabolite exchange in reef‐building corals as multi‐partner metabolic networks. Environ Microbiol 2020;22:1675–87.PubMed 

Google Scholar 
Kimes NE, Van Nostrand JD, Weil E, Zhou J, Morris PJ. Microbial functional structure of Montastraea faveolata, an important Caribbean reef‐building coral, differs between healthy and yellow‐band diseased colonies. Environ Microbiol. 2010;12:541–56.CAS 
PubMed 

Google Scholar 
Neave MJ, Apprill A, Ferrier-Pagès C, Voolstra CR. Diversity and function of prevalent symbiotic marine bacteria in the genus Endozoicomonas. Appl Environ Micro. 2016;100:8315–24.CAS 

Google Scholar 
Neave MJ, Michell CT, Apprill A, Voolstra CR. Endozoicomonas genomes reveal functional adaptation and plasticity in bacterial strains symbiotically associated with diverse marine hosts. Sci Rep. 2017;7:1–12.
Google Scholar 
Krediet CJ, Ritchie KB, Alagely A, Teplitski M. Members of native coral microbiota inhibit glycosidases and thwart colonization of coral mucus by an opportunistic pathogen. ISME J. 2013;7:980–90.CAS 
PubMed 

Google Scholar 
Raina J-B, Tapiolas D, Motti CA, Foret S, Seemann T, Tebben J, et al. Isolation of an antimicrobial compound produced by bacteria associated with reef-building corals. PeerJ 2016;4:e2275.PubMed 
PubMed Central 

Google Scholar 
Diaz JM, Hansel CM, Apprill A, Brighi C, Zhang T, Weber L, et al. Species-specific control of external superoxide levels by the coral holobiont during a natural bleaching event. Nat Commun. 2016;7:1–10.
Google Scholar 
Dunlap WC, Shick JM. Ultraviolet radiation‐absorbing mycosporine‐like amino acids in coral reef organisms: a biochemical and environmental perspective. J Phycol. 1998;34:418–30.
Google Scholar 
Webster NS, Smith LD, Heyward AJ, Watts JE, Webb RI, Blackall LL, et al. Metamorphosis of a scleractinian coral in response to microbial biofilms. Appl Environ Microbiol. 2004;70:1213–21.CAS 
PubMed 
PubMed Central 

Google Scholar 
Gómez-Lemos LA, Doropoulos C, Bayraktarov E, Diaz-Pulido G. Coralline algal metabolites induce settlement and mediate the inductive effect of epiphytic microbes on coral larvae. Sci Rep. 2018;8:1–11.
Google Scholar 
Pernice M, Raina J-B, Rädecker N, Cárdenas A, Pogoreutz C, Voolstra CR. Down to the bone: the role of overlooked endolithic microbiomes in reef coral health. ISME J. 2020;14:325–34.PubMed 

Google Scholar 
Ricci F, Marcelino VR, Blackall LL, Kühl M, Medina M, Verbruggen H. Beneath the surface: community assembly and functions of the coral skeleton microbiome. Microbiome 2019;7:1–10.
Google Scholar 
Marcelino VR, Verbruggen H. Multi-marker metabarcoding of coral skeletons reveals a rich microbiome and diverse evolutionary origins of endolithic algae. Sci Rep. 2016;6:1–9.
Google Scholar 
Verbruggen H, Marcelino VR, Guiry MD, Cremen MCM, Jackson CJ. Phylogenetic position of the coral symbiont Ostreobium (Ulvophyceae) inferred from chloroplast genome data. J Phycol. 2017;53:790–803.CAS 
PubMed 

Google Scholar 
Del Campo J, Pombert J-F, Šlapeta J, Larkum A, Keeling PJ. The ‘other’coral symbiont: Ostreobium diversity and distribution. ISME J 2017;11:296–9.PubMed 

Google Scholar 
Massé A, Domart-Coulon I, Golubic S, Duché D, Tribollet A. Early skeletal colonization of the coral holobiont by the microboring Ulvophyceae Ostreobium sp. Sci Rep. 2018;8:1–11.
Google Scholar 
Halldal P. Photosynthetic capacities and photosynthetic action spectra of endozoic algae of the massive coral Favia. Biol Bull. 1968;134:411–24.CAS 

Google Scholar 
Fork D, Larkum A. Light harvesting in the green alga Ostreobium sp., a coral symbiont adapted to extreme shade. Mar Biol. 1989;103:381–5.
Google Scholar 
Fine M, Steindler L, Loya Y. Endolithic algae photoacclimate to increased irradiance during coral bleaching. Mar Freshw Res. 2004;55:115–21.CAS 

Google Scholar 
Fine M, Roff G, Ainsworth T, Hoegh-Guldberg O. Phototrophic microendoliths bloom during coral “white syndrome”. Coral Reefs. 2006;25:577–81.
Google Scholar 
Galindo-Martínez CT, Weber M, Avila-Magaña V, Enríquez S, Kitano H, Medina M, et al. The role of the endolithic alga Ostreobium spp. during coral bleaching recovery. Sci Rep. 2022;12:1–12.
Google Scholar 
Fine M, Loya Y. Endolithic algae: an alternative source of photoassimilates during coral bleaching. Proc R Soc B Biol Sci. 2002;269:1205–10.
Google Scholar 
Schlichter D, Zscharnack B, Krisch H. Transfer of photoassimilates from endolithic algae to coral tissue. Naturwissenschaften 1995;82:561–4.CAS 

Google Scholar 
Sangsawang L, Casareto BE, Ohba H, Vu HM, Meekaew A, Suzuki T, et al. 13C and 15N assimilation and organic matter translocation by the endolithic community in the massive coral Porites lutea. R Soc Open Sci. 2017;4:171201.PubMed 
PubMed Central 

Google Scholar 
Marcelino VR, Morrow KM, van Oppen MJ, Bourne DG, Verbruggen H. Diversity and stability of coral endolithic microbial communities at a naturally high pCO2 reef. Mol Ecol. 2017;26:5344–57.CAS 
PubMed 

Google Scholar 
Marcelino VR, Van Oppen MJ, Verbruggen H. Highly structured prokaryote communities exist within the skeleton of coral colonies. ISME J. 2018;12:300–3.PubMed 

Google Scholar 
Yang S-H, Tandon K, Lu C-Y, Wada N, Shih C-J, Hsiao SS-Y, et al. Metagenomic, phylogenetic, and functional characterization of predominant endolithic green sulfur bacteria in the coral Isopora palifera. Microbiome 2019;7:1–13.
Google Scholar 
Ferrer L, Szmant A, editors. Nutrient regeneration by the endolithic community in coral skeletons. Proceedings of the 6th International Coral Reef Symposium; 1988: AIMS Townsville, Australia.Eakin CM, Devotta D, Heron S, Connolly S, Liu G, Geiger E, et al. The 2014-17 global coral bleaching event: The most severe and widespread coral reef destruction. Research Square. 2022. https://doi.org/10.21203/rs.3.rs-1555992/v1Article 

Google Scholar 
Hughes TP, Kerry JT, Álvarez-Noriega M, Álvarez-Romero JG, Anderson KD, Baird AH, et al. Global warming and recurrent mass bleaching of corals. Nature 2017;543:373–7.CAS 
PubMed 

Google Scholar 
Hughes TP, Anderson KD, Connolly SR, Heron SF, Kerry JT, Lough JM, et al. Spatial and temporal patterns of mass bleaching of corals in the Anthropocene. Science 2018;359:80–3.CAS 
PubMed 

Google Scholar 
Hughes TP, Kerry JT, Baird AH, Connolly SR, Dietzel A, Eakin CM, et al. Global warming transforms coral reef assemblages. Nature 2018;556:492–6.CAS 
PubMed 

Google Scholar 
Veron J, Stafford-Smith M, Corals of the World, Volumes 1-3. Australian Institute of Marine Science. Odyssey Publishing; 2000.Brown B, Dunne R, Phongsuwan N, Patchim L, Hawkridge J. The reef coral Goniastrea aspera: a ‘winner’becomes a ‘loser’during a severe bleaching event in Thailand. Coral Reefs. 2014;33:395–401.
Google Scholar 
Klepac C, Barshis D. Reduced thermal tolerance of massive coral species in a highly variable environment. Proc R Soc B Biol Sci. 2020;287:20201379.CAS 

Google Scholar 
Nicolas R, Evensen CR, Voolstra M, Fine G, Perna C, Buitrago-López A, et al. Empirically derived thermal thresholds of four coral species along the Red Sea using a portable and standardized experimental approach. Coral Reefs. 2022;41:239–52. https://doi.org/10.1007/s00338-022-02233-yArticle 

Google Scholar 
Madin JS, Anderson KD, Andreasen MH, Bridge TC, Cairns SD, Connolly SR, et al. The Coral Trait Database, a curated database of trait information for coral species from the global oceans. Sci Data. 2016;3:1–22.
Google Scholar 
Roth F, Karcher DB, Rädecker N, Hohn S, Carvalho S, Thomson T, et al. High rates of carbon and dinitrogen fixation suggest a critical role of benthic pioneer communities in the energy and nutrient dynamics of coral reefs. Funct Ecol. 2020;34:1991–2004.
Google Scholar 
Harrison PJ, Waters RE, Taylor F. A broad spectrum artificial sea water medium for coastal and open ocean phytoplankton. J Phycol. 1980;16:28–35.
Google Scholar 
Andersson AF, Lindberg M, Jakobsson H, Bäckhed F, Nyrén P, Engstrand L. Comparative analysis of human gut microbiota by barcoded pyrosequencing. PLoS One. 2008;3:e2836.PubMed 
PubMed Central 

Google Scholar 
Bayer T, Neave MJ, Alsheikh-Hussain A, Aranda M, Yum LK, Mincer T, et al. The microbiome of the Red Sea coral Stylophora pistillata is dominated by tissue-associated Endozoicomonas bacteria. Appl Environ Microbiol. 2013;79:4759–62.CAS 
PubMed 
PubMed Central 

Google Scholar 
Callahan BJ, McMurdie PJ, Rosen MJ, Han AW, Johnson AJA, Holmes SP. DADA2: high-resolution sample inference from Illumina amplicon data. Nat Methods. 2016;13:581–3.CAS 
PubMed 
PubMed Central 

Google Scholar 
Quast C, Pruesse E, Yilmaz P, Gerken J, Schweer T, Yarza P, et al. The SILVA ribosomal RNA gene database project: improved data processing and web-based tools. Nucleic Acids Res. 2012;41:D590–D6.PubMed 
PubMed Central 

Google Scholar 
Wickham H. ggplot2. Wiley Interdiscip Rev Comput Stat. 2011;3:180–5.
Google Scholar 
McMurdie PJ, Holmes S. phyloseq: an R package for reproducible interactive analysis and graphics of microbiome census data. PLoS One. 2013;8:e61217.CAS 
PubMed 
PubMed Central 

Google Scholar 
Dixon P. The vegan package. J Veg Sci. 2003;14:927–30.
Google Scholar 
Lin H, Peddada SD. Analysis of compositions of microbiomes with bias correction. Nat Commun. 2020;11:1–11.CAS 

Google Scholar 
Bolger AM, Lohse M, Usadel B. Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics 2014;30:2114–20.CAS 
PubMed 
PubMed Central 

Google Scholar 
Li D, Luo R, Liu C-M, Leung C-M, Ting H-F, Sadakane K, et al. MEGAHIT v1. 0: a fast and scalable metagenome assembler driven by advanced methodologies and community practices. Methods 2016;102:3–11.CAS 
PubMed 

Google Scholar 
Nurk S, Meleshko D, Korobeynikov A, Pevzner PA. metaSPAdes: a new versatile metagenomic assembler. Genome Res. 2017;27:824–34.CAS 
PubMed 
PubMed Central 

Google Scholar 
Peng Y, Leung HC, Yiu S-M, Chin FY. Meta-IDBA: a de Novo assembler for metagenomic data. Bioinformatics 2011;27:i94–i101.CAS 
PubMed 
PubMed Central 

Google Scholar 
Gurevich A, Saveliev V, Vyahhi N, Tesler G. QUAST: quality assessment tool for genome assemblies. Bioinformatics 2013;29:1072–5.CAS 
PubMed 
PubMed Central 

Google Scholar 
Hyatt D, Chen G-L, LoCascio PF, Land ML, Larimer FW, Hauser LJ. Prodigal: prokaryotic gene recognition and translation initiation site identification. BMC Bioinform. 2010;11:1–11.
Google Scholar 
Patro R, Duggal G, Love MI, Irizarry RA, Kingsford C. Salmon provides fast and bias-aware quantification of transcript expression. Nat Methods. 2017;14:417–9.CAS 
PubMed 
PubMed Central 

Google Scholar 
Menzel P, Ng KL, Krogh A. Fast and sensitive taxonomic classification for metagenomics with Kaiju. Nat Commun. 2016;7:1–9.
Google Scholar 
Aramaki T, Blanc-Mathieu R, Endo H, Ohkubo K, Kanehisa M, Goto S, et al. KofamKOALA: KEGG ortholog assignment based on profile HMM and adaptive score threshold. Bioinformatics 2020;36:2251–2.CAS 
PubMed 

Google Scholar 
Hill MO. Diversity and evenness: a unifying notation and its consequences. Ecology 1973;54:427–32.
Google Scholar 
Bates D, Sarkar D, Bates MD, Matrix L. The lme4 package. R Package Version. 2007;2:74.
Google Scholar 
Mandal S, Van Treuren W, White RA, Eggesbø M, Knight R, Peddada SD. Analysis of composition of microbiomes: a novel method for studying microbial composition. Micro Ecol Health Dis. 2015;26:27663.
Google Scholar 
Rivera-Pinto J, Egozcue JJ, Pawlowsky-Glahn V, Paredes R, Noguera-Julian M, Calle ML. Balances: a new perspective for microbiome analysis. mSystems. 2018;3:e00053–18.CAS 
PubMed 
PubMed Central 

Google Scholar 
Langmead B, Salzberg SL. Fast gapped-read alignment with Bowtie 2. Nat Methods. 2012;9:357.CAS 
PubMed 
PubMed Central 

Google Scholar 
Kang DD, Li F, Kirton E, Thomas A, Egan R, An H, et al. MetaBAT 2: an adaptive binning algorithm for robust and efficient genome reconstruction from metagenome assemblies. PeerJ 2019;7:e7359.PubMed 
PubMed Central 

Google Scholar 
Alneberg J, Bjarnason BS, De Bruijn I, Schirmer M, Quick J, Ijaz UZ, et al. Binning metagenomic contigs by coverage and composition. Nat Methods. 2014;11:1144–6.CAS 
PubMed 

Google Scholar 
Wu Y-W, Simmons BA, Singer SW. MaxBin 2.0: an automated binning algorithm to recover genomes from multiple metagenomic datasets. Bioinformatics 2016;32:605–7.CAS 
PubMed 

Google Scholar 
Parks DH, Imelfort M, Skennerton CT, Hugenholtz P, Tyson GW. CheckM: assessing the quality of microbial genomes recovered from isolates, single cells, and metagenomes. Genome Res. 2015;25:1043–55.CAS 
PubMed 
PubMed Central 

Google Scholar 
Uritskiy GV, DiRuggiero J, Taylor J. MetaWRAP—a flexible pipeline for genome-resolved metagenomic data analysis. Microbiome 2018;6:1–13.
Google Scholar 
Olm MR, Brown CT, Brooks B, Banfield JF. dRep: a tool for fast and accurate genomic comparisons that enables improved genome recovery from metagenomes through de-replication. ISME J. 2017;11:2864–8.CAS 
PubMed 
PubMed Central 

Google Scholar 
Chaumeil P-A, Mussig AJ, Hugenholtz P, Parks DH. GTDB-Tk: a toolkit to classify genomes with the Genome Taxonomy Database. Bioinformatics 2020;36:1925–7.CAS 

Google Scholar 
Parks DH, Chuvochina M, Waite DW, Rinke C, Skarshewski A, Chaumeil P-A, et al. A standardized bacterial taxonomy based on genome phylogeny substantially revises the tree of life. Nat Biotechnol. 2018;36:996–1004.CAS 
PubMed 

Google Scholar 
Seemann T. Prokka: rapid prokaryotic genome annotation. Bioinformatics 2014;30:2068–9.CAS 
PubMed 

Google Scholar 
Zhou Z, Tran P, Briester AM, Liu Y, Kieft K, Cowley ES, et al. METABOLIC: high-throughput profiling of microbial genomes for functional traits, metabolism, biogeochemistry, and community-scale functional networks. Microbiome. 2022;10:33 https://doi.org/10.1186/s40168-021-01213-8CAS 
Article 
PubMed 
PubMed Central 

Google Scholar 
Na S-I, Kim YO, Yoon S-H, Ha S-M, Baek I, Chun J. UBCG: up-to-date bacterial core gene set and pipeline for phylogenomic tree reconstruction. J Microbiol. 2018;56:280–5.CAS 
PubMed 

Google Scholar 
Morel J, Jay S, Féret J-B, Bakache A, Bendoula R, Carreel F, et al. Exploring the potential of PROCOSINE and close-range hyperspectral imaging to study the effects of fungal diseases on leaf physiology. Sci Rep. 2018;8:1–13.CAS 

Google Scholar 
Calamita F, Imran HA, Vescovo L, Mekhalfi ML, La, Porta N. Early identification of root rot disease by using hyperspectral reflectance: the case of pathosystem Grapevine/Armillaria. Remote Sens. 2021;13:2436.
Google Scholar 
Brumfield KD, Huq A, Colwell RR, Olds JL, Leddy MB. Microbial resolution of whole-genome shotgun and 16S amplicon metagenomic sequencing using publicly available NEON data. PLoS One. 2020;15:e0228899.PubMed 
PubMed Central 

Google Scholar 
Khachatryan L, de Leeuw RH, Kraakman ME, Pappas N, Te Raa M, Mei H, et al. Taxonomic classification and abundance estimation using 16S and WGS—A comparison using controlled reference samples. Forensic Sci Int Genet. 2020;46:102257.CAS 
PubMed 

Google Scholar 
Cardénas A, Voolstra C. 75 Coral Endolith Bacterial Genomes (MAGs) from Red Sea corals Goniastrea edwardsi and Porites lutea (Version 1) [Data set]. Zenodo. 2021. https://doi.org/10.5281/zenodo.5606932Article 

Google Scholar 
Branson O, Bonnin EA, Perea DE, Spero HJ, Zhu Z, Winters M, et al. Nanometer-scale chemistry of a calcite biomineralization template: Implications for skeletal composition and nucleation. Proc Natl Acad Sci. 2016;113:12934–9.CAS 
PubMed 
PubMed Central 

Google Scholar 
Sauvage T, Schmidt WE, Suda S, Fredericq S. A metabarcoding framework for facilitated survey of endolithic phototrophs with tufA. BMC Ecol. 2016;16:1–21.
Google Scholar 
Wegley L, Edwards R, Rodriguez‐Brito B, Liu H, Rohwer F. Metagenomic analysis of the microbial community associated with the coral Porites astreoides. Environ Microbiol. 2007;9:2707–19.CAS 
PubMed 

Google Scholar 
Robbins S, Song W, Engelberts J, Glasl B, Slaby BM, Boyd J, et al. A genomic view of the microbiome of coral reef demosponges. ISME J 2021;15:1641–54.CAS 
PubMed 
PubMed Central 

Google Scholar 
Yang S-Y, Lu C-Y, Tang S-L, Das RR, Sakai K, Yamashiro H, et al. Effects of ocean acidification on coral endolithic bacterial communities in Isopora palifera and Porites lobata. Front Mar Sci. 2020;7:603293.
Google Scholar 
Yang SH, Lee ST, Huang CR, Tseng CH, Chiang PW, Chen CP, et al. Prevalence of potential nitrogen‐fixing, green sulfur bacteria in the skeleton of reef‐building coral Isopora palifera. Limnol Oceanogr. 2016;61:1078–86.
Google Scholar 
Cai L, Zhou G, Tian R-M, Tong H, Zhang W, Sun J, et al. Metagenomic analysis reveals a green sulfur bacterium as a potential coral symbiont. Sci Rep. 2017;7:1–11.
Google Scholar 
Kühl M, Holst G, Larkum AW, Ralph PJ. Imaging of oxygen dynamics within the endolithic algal community of the massive coral Porites Lobata. J Phycol. 2008;44:541–50.PubMed 

Google Scholar 
Roberty S, Bailleul B, Berne N, Franck F, Cardol P. PSI Mehler reaction is the main alternative photosynthetic electron pathway in Symbiodinium sp., symbiotic dinoflagellates of cnidarians. N. Phytol. 2014;204:81–91.CAS 

Google Scholar 
Shigeoka S, Ishikawa T, Tamoi M, Miyagawa Y, Takeda T, Yabuta Y, et al. Regulation and function of ascorbate peroxidase isoenzymes. J Exp Bot. 2002;53:1305–19.CAS 
PubMed 

Google Scholar 
Roberty S, Fransolet D, Cardol P, Plumier J-C, Franck F. Imbalance between oxygen photoreduction and antioxidant capacities in Symbiodinium cells exposed to combined heat and high light stress. Coral Reefs. 2015;34:1063–73.
Google Scholar 
Petersen JM, Zielinski FU, Pape T, Seifert R, Moraru C, Amann R, et al. Hydrogen is an energy source for hydrothermal vent symbioses. Nature 2011;476:176–80.CAS 
PubMed 

Google Scholar 
McCollom T, Amend J. A thermodynamic assessment of energy requirements for biomass synthesis by chemolithoautotrophic micro‐organisms in oxic and anoxic environments. Geobiology 2005;3:135–44.CAS 

Google Scholar 
Heijnen J, Van Dijken J. In search of a thermodynamic description of biomass yields for the chemotrophic growth of microorganisms. Biotechnol Bioeng. 1992;39:833–58.CAS 
PubMed 

Google Scholar 
Bar-Even A, Noor E, Milo R. A survey of carbon fixation pathways through a quantitative lens. J Exp Bot. 2012;63:2325–42.CAS 
PubMed 

Google Scholar 
Schulze E-D, Mooney HA, Biodiversity and ecosystem function: Springer Science & Business Media; 2012.Lawton JH, Brown VK, Redundancy in ecosystems. Biodiversity and Ecosystem Function: Springer; 1994. p. 255–70.Mori AS, Furukawa T, Sasaki T. Response diversity determines the resilience of ecosystems to environmental change. Biol Rev. 2013;88:349–64.PubMed 

Google Scholar 
Nyström M. Redundancy and response diversity of functional groups: implications for the resilience of coral reefs. Ambio 2006;35:30–5.PubMed 

Google Scholar 
Rädecker N, Pogoreutz C, Gegner HM, Cárdenas A, Roth F, Bougoure J, et al. Heat stress destabilizes symbiotic nutrient cycling in corals. Proc Natl Acad Sci. 2021;118:e2022653118.PubMed 
PubMed Central 

Google Scholar 
Ziegler M, Grupstra CG, Barreto MM, Eaton M, BaOmar J, Zubier K, et al. Coral bacterial community structure responds to environmental change in a host-specific manner. Nat Commun. 2019;10:1–11.CAS 

Google Scholar 
Dikou A, Van, Woesik R. Survival under chronic stress from sediment load: spatial patterns of hard coral communities in the southern islands of Singapore. Mar Pollut Bull. 2006;52:1340–54.CAS 
PubMed 

Google Scholar 
Hennige SJ, Smith DJ, Walsh S-J, McGinley MP, Warner ME, Suggett DJ. Acclimation and adaptation of scleractinian coral communities along environmental gradients within an Indonesian reef system. J Exp Mar Biol Ecol. 2010;391:143–52.
Google Scholar 
Cárdenas A, Neave MJ, Haroon MF, Pogoreutz C, Rädecker N, Wild C, et al. Excess labile carbon promotes the expression of virulence factors in coral reef bacterioplankton. ISME J. 2018;12:59–76.PubMed 

Google Scholar 
Cárdenas A, Ye J, Ziegler M, Payet JP, McMinds R, Thurber RV, et al. Coral-associated viral assemblages from the Central Red Sea align with host species and contribute to holobiont genetic diversity. Front Microbiol. 2020;11:572534.PubMed 
PubMed Central 

Google Scholar 
McCook GD-PLJ. The fate of bleached corals: patterns and dynamics of algal recruitment. Mar Ecol Prog Ser. 2002;232:115–28.
Google Scholar 
Reshef L, Koren O, Loya Y, Zilber‐Rosenberg I, Rosenberg E. The coral probiotic hypothesis. Environ Microbiol. 2006;8:2068–73.CAS 
PubMed 

Google Scholar 
Voolstra CR, Ziegler M. Adapting with microbial help: microbiome flexibility facilitates rapid responses to environmental change. BioEssays 2020;42:2000004.
Google Scholar 
Rosenberg E, Zilber-Rosenberg I. The hologenome concept of evolution after 10 years. Microbiome 2018;6:1–14.
Google Scholar 
Wiedenmann J, D’Angelo C, Smith EG, Hunt AN, Legiret F-E, Postle AD, et al. Nutrient enrichment can increase the susceptibility of reef corals to bleaching. Nat Clim Change. 2013;3:160–4.CAS 

Google Scholar 
DeCarlo TM, Gajdzik L, Ellis J, Coker DJ, Roberts MB, Hammerman NM, et al. Nutrient-supplying ocean currents modulate coral bleaching susceptibility. Sci Adv. 2020;6:eabc5493.PubMed 
PubMed Central 

Google Scholar 
Pogoreutz C, Rädecker N, Cardenas A, Gärdes A, Voolstra CR, Wild C. Sugar enrichment provides evidence for a role of nitrogen fixation in coral bleaching. Glob Change Biol 2017;23:3838–48.
Google Scholar  More