Contrasting metabolic strategies of two co-occurring deep-sea octocorals

The feeding biology of deep-sea octocorals remains poorly understood, as attention is more often directed to reef building corals. The present study focused on two common deep-water octocoral species in the Azores Archipelago, Dentomuricea aff. meteor and Viminella flagellum, aiming at determining their ability to exploit different food sources. We adopted an experimental approach, with three different food sources, including live phytoplankton, live zooplankton and dissolved organic matter (DOM), that were artificially enriched with 13C and 15N (C and N tracers). The presence of tracers was subsequently followed in the coral tissue, C respiration and particulate organic C and N (POC and PON) release. In both species, feeding with zooplankton resulted in significantly higher incorporation of tracers in all measured variables, compared to the other food sources, highlighting the importance of zooplankton for major physiological processes. Our results revealed contrasting metabolic strategies between the two species, with D. aff. meteor acquiring higher amounts of prey and allocating higher percentage to respiration and release of POC and PON than V. flagellum. Such metabolic differences can shape species fitness and distributions and have further ecological implications on the ecosystem function of communities formed by different octocoral species.


Results
All utilized food sources were significantly enriched above background, i.e. above non-labelled food (Table 1). Target C concentrations within the aquaria were successfully achieved in the case of the DOM and ZOO treatment, however they were below the target for the PHYTO treatment, by 34% for D. aff. meteor and 44% for V. flagellum respectively (Table 1).
Coral fragments incorporated tracer C and N from all food sources in their tissue (Fig. 2). All results on tracers are presented as average ± standard deviation. Both species incorporated significantly higher tracer C and N under the ZOO treatment followed by the DOC treatment (Fig. 2, Table 2). Coral fragments of D. aff. meteor and V. flagellum under the ZOO treatment incorporated 422% and 453% more tracer C than under the PHYTO treatment, respectively. Viminella flagellum displayed lower tracer incorporation compared to D. aff. meteor in all treatments, reaching on average 64% lower C incorporation and 70% lower N incorporation than D. aff. meteor (Fig. 2). Total tracer C incorporation varied between 0.19 ± 0.24 μmol tracer C mmol tissue C −1 under the PHYTO treatment to 80.2 ± 24.5 μmol tracer C mmol tissue C −1 under the ZOO treatment for D. aff. meteor and from 0.05 ± 0.02 μmol tracer C mmol tissue C −1 under the PHYTO treatment to 21.8 ± 13.9 μmol tracer C mmol tissue C −1 under the ZOO treatment for V. flagellum. www.nature.com/scientificreports/ In both species, oxygen consumption was significantly higher in the ZOO treatment compared to the other food treatments (Table 2), reaching on average 0.308 ± 0.042 μmol Ο 2 mmol tissue C −1 h −1 for D. aff. meteor and 0.151 ± 0.036 μmol Ο 2 mmol tissue C −1 h −1 for V. flagellum (Fig. 3). Overall, oxygen consumption was almost two times higher in fragments of D. aff. meteor compared to V. flagellum.
Both species utilized tracer C and N derived from the provided food for C respiration, POC and PON release (Fig. 4). Tracer C respiration was significantly higher under the ZOO treatment for both species (Fig. 4, Table 2). Similarly, for both species tracer POC and PON release were higher under the ZOO treatment while fragments under the DOM and PHYTO treatments displayed very low average values of POC and PON release, which did not differ significantly from zero ( Table 2).
To take into account differences in the provided C quantity, the amount of utilized C tracer is provided as a percentage of the total provided C in Fig. 5. For D. aff. meteor 45% of the provided zooplankton-derived tracer C could be traced back in tissue incorporation, DIC respiration and POC release, with phytoplankton-derived and DOM-derived tracer C reaching 1.08% and 0.26%, respectively (Fig. 5). Under the PHYTO treatment, corals utilized most of the C tracer for POC release, while in the two other treatments they utilized most of the C tracer for tissue incorporation, followed by tracer C respiration (Fig. 5). In V. flagellum, a smaller percentage of Table 1. Characteristics of enriched food sources used in feeding experiments of Dentomuricea aff. meteor and Viminella flagellum, including average ± SD of the created C concentrations in the aquaria, total provided C per coral fragment and per coral organic carbon, fractional abundance (F 13 , F 15 ) and carbon to nitrogen (C/N) ratio. CC: Chaetoceros calcitrans, NG: Nannochloropsis gaditana, BP: Branchionus plicatilis.

Discussion
To our knowledge this is the first study addressing the feeding biology of deep octocoral species by employing an experimental approach with the use of stable isotopes. The strong differences in tracer C and N incorporation among the food treatments suggest a higher efficiency in capture and ingestion of zooplankton prey when this resource is available, compared to the other food treatments. Previous studies have highlighted zooplankton as an important dietary component of deep-sea scleractinians [28][29][30][31] . However, only a few deep-sea octocorals, mainly of the family Primnoidae, have been shown to base their diet on microzooplankton 32,33 while most studied species rely mostly on phytodetritus and particulate organic matter (POM) 28,34 . Our results provide a direct demonstration of the importance of zooplankton for some deep-sea octocorals and highlight that they might be more selective than previously thought. In areas with strong hydrodynamics, such as seamounts, fresh phytoplankton can be directly transported to great depths by rapid downwelling and tidal waves 35,36 , and therefore it might be often available for the two target octocorals. The closely positioned pinnules of octocoral tentacles and their relatively weak nematocysts generally suggest herbivory 24,37,38 . Moreover, in previous studies, both target octocorals have been shown to occupy low trophic levels, placed between primary and secondary consumers 39 . Thus, the lower incorporation Table 2. Coefficients of constructed models to explore the effect of four food treatments (FAST: fasting; PHYTO: diatom Chaetoceros calcitrans; DOM: dissolved organic matter; ZOO: rotifer Branchionus plicatilis) on each dependent variable after analysis of collected data on tissue, respiration and excretion of two octocoral species Dentomuricea aff. meteor and Viminella flagellum. If food source was excluded from the respective model during model construction, it was assumed it had no significant effect on the response variable in question and therefore no coefficients are provided. SE: Standard error. www.nature.com/scientificreports/ of phytoplankton was unexpected. It is possible that during the experiment, coral fragments fed additionally on small particles (< 1 μm) that passed the filtration system. This could explain the lower DOM and phytoplankton utilization, but it further supports the hypothesis that both species can display selective feeding. The feeding behavior of the two target octocorals is in strong contrast to the results of similar feeding experiments with the deep-sea scleractinian Lophelia pertusa (recently synonymized to Desmophyllum pertusum 40 )   www.nature.com/scientificreports/ which displayed a rather unselective feeding behavior under feeding with similar food sources and flow velocities 41 . These differences might be due to the use of processed (freeze-dried) versus live prey. The use of isotopically enriched live prey has been used before to study the feeding preferences of octocorals 42,43 and scleractinians 44,45 . It allows more realistic comparisons compared to dry food, as it takes into consideration both the capture and ingestion ability of the study species. On the other hand, it includes a considerable error in determining and standardizing provided C quantities since C content can vary among culture batches. In the present study, the available C in the aquaria of the DOM and ZOO treatments was 30-40% higher than in the PHYTO treatment, thus a proportionately higher utilization of the DOM and ZOO food sources was expected. While this can explain the small differences in tracer C utilization between the PHYTO and DOM treatments, it cannot explain the disproportionally larger tracer utilization under the ZOO treatment, strongly indicating more efficient feeding on zooplankton.
Since zooplankton dynamics in the study area follow seasonal phytoplankton productivity cycles 46 , it is likely that this important food source displays strong seasonal fluctuations. Both target species displayed the ability to utilize food sources variable in size and composition, therefore during the rest of the year it is very likely that the two species sustain their metabolism through feeding on other sources such as phytoplankton and DOM. Dissolved organic matter has proven to be an important food source for cold-water scleractinians when particulate food sources are scarce 47 . The ability to utilize DOM can be very important in oligotrophic deep-sea environments where food availability is seasonal and is likely to be further affected by climate change 48 . A number of previous studies have highlighted the variable and seasonally-dependent diet of octocoral species in temperate ecosystems 19,49,50 while similar seasonality has been also reported in benthic Antarctic ecosystems 27,51 .
The two target octocorals displayed higher oxygen consumption and tracer C respiration upon feeding with zooplankton, highlighting the importance of zooplankton to meet their respiratory and metabolic demands. Similar results have been reported for Desmophyllum dianthus which displayed lower oxygen respiration, calcification and TOC release after exclusion of zooplankton from available food sources 31 . Because of the seasonal availability of zooplankton and its importance for tissue incorporation and metabolism, it appears likely that physiological processes that require the development of C and N rich tissues, such as growth and reproduction, may also undergo strong seasonality. This is a common phenomenon for octocorals in temperate areas, which display seasonal cycles in their biochemical levels e.g. 52 and often pass through periods of metabolic dormancy in summer months when available seston is scarce and temperature rises 16,53 . Dentomuricea aff. meteor and V. flagellum display gamete presence all year round with frequent seasonal peaks in spring and autumn (M. Rakka, unpubl.), which could be related to higher zooplankton availability 46 , but more studies on abiotic conditions and physiological cycles are essential to unravel their ecophysiology.
Fragments fed with zooplankton also displayed higher tracer POC and PON release in both species, which in corals is associated with mucus production, essential in processes such as feeding, cleansing and protection from epibionts and pathogens 54 . Mucus production can be extremely important to protect the corals against mechanical and chemical disturbance due to bottom trawling, oil extraction and mining (e.g. mine tailing and drill cutting resuspension) [55][56][57] . Moreover, coral mucus has been identified as an ecologically important element for CWC communities, since it enhances microbial activity and therefore mineralization, recycling and overall ecosystem productivity 54,58 . The increased POC release under the ZOO treatment showcases how feeding on zooplankton can enhance the contribution of octocoral species to C recycling and highlights their importance for bentho-pelagic coupling. www.nature.com/scientificreports/ Overall, D. aff. meteor appeared to acquire higher percentage of the provided C compared to V. flagellum. Moreover, when fed with the most effectively utilized food source (zooplankton), D. aff. meteor allocated a higher proportion of captured C to respiration and POC release compared to V. flagellum. This highlights the different strategies adopted by the two species. While V. flagellum appeared more conservative in resource allocation, storing most of the captured C in tissue, and minimizing losses, the pattern displayed by D. aff. meteor is indicative of "sloppy feeding", in which high amounts of the captured C are lost during the feeding process 59,60 . Feeding behaviour and metabolic rates can be influenced by an array of factors, such as environmental variables, physiological status and life stage. The two target species were haphazardly collected, maintained in similar conditions and were expected to be in similar stages in their reproductive cycle (M. Rakka, unpubl.). Thus, variations in colony health, age and maturity are more likely to explain the observed variance within each species, while the marked difference between the two octocorals may be attributed to species-specific characteristics, such as morphology and growth pattern 61 . The species D. aff. meteor has a branching pattern with high surface to volume ratio and possesses a large number of small polyps which can increase both capture rates and metabolic costs 14,62 . On the other hand, V. flagellum displays an erect growing pattern with bigger polyps and a lower surface to volume ratio that may have lower maintenance costs.
Both target species incorporated a lower amount of tracer from phytoplankton compared to L. pertusa 63 . In contrast, fragments fed with zooplankton displayed higher tracer incorporation than L. pertusa, i.e. 10 times higher for D. aff. meteor and two times higher for V. flagellum. Although these differences might be attributed to the use of dry versus live prey, they highlight that resource acquisition strategies are species-specific, as also demonstrated for several tropical and cold-water coral (CWC) species 18,28,33 . Metabolism is tightly connected to the ecological niche of a species, and different responses to food supply can explain distributions of species and species assemblages 64 . For example, D. aff. meteor, due to its metabolic strategy presented herein, is expected to have an advantage under high food concentration, however, it is unlikely to outperform V. flagellum under low food conditions. Similarly to reefs formed by cold-water scleractinian species 65,66 , the rich organic excretion of octocorals promotes organic cycling and plays an important role supporting a diverse community of associated fauna in coral gardens and adjacent deep-sea communities 67 . Taking into account that the two species displayed different strategies in respect to the respired and released C, we hypothesize that the communities dominated by D. aff. meteor are likely to be characterized by higher C and N recycling whereas communities of V. flagellum will likely have higher residence time of C and N in the coral tissue, with further consequences for the local C cycle. Sloppy feeding is known as an important behavior for the support of C cycles, fueling the microbial loop and supporting local food webs 68,69 . Thus, the role of these species to local and global marine biogeochemical cycles should be further investigated.
In conclusion, the present study provides important knowledge on the resource utilization and metabolic strategies of two important habitat forming octocorals, that can help understand patterns at the species, population and community level 70 . Species distribution modelling has shed light to the distribution of different deepsea coral groups, including octocorals 71 , but comparatively little effort has been made to delve further into the biological and physiological characteristics which shape these distributions 72,73 . Taking into account these traits will not only improve predictions on species occurrences and help to identify priority areas for conservation and management 74 , but will also provide a more robust understanding of the ecology of deep-sea corals. Coral resource use and metabolism is likely to change under future conditions of increased seawater temperature, stratification of water masses and consequent reduction in the quantity and quality of POM flux to the seafloor 48 . Further studies on the ecophysiology of octocoral species under present and future scenarios of climate change are therefore essential to improve our understanding of the distribution and ecological function of deep-sea communities.

Methods
Target species. The species Dentomuricea aff. meteor is a fan-shaped alcyonacean of the family Plexauridae.
Its known distribution is limited to seamounts close to the Mid-Atlantic ridge where it is typically encountered between 200 and 400 m depth 7 . It can reach heights of up to 1.5 m and can create dense monospecific or mixed populations with other species (Fig. 1a).
Viminella flagellum 75 is a whip coral of the family Ellisellidae. It creates monopodial colonies without branches which can grow up to 3 m height 76 . Its distribution includes the eastern North Atlantic coast, islands of the Macaronesia and the Mediterranean Sea [77][78][79] . It is usually encountered between 120 and 500 m depth 77,78 and can form monospecific or mixed coral aggregations, often with branching octocorals such as D. aff. meteor, Callogorgia verticillata and Acanthogorgia armata 80 . Colony collection and maintenance. Live colonies of both V. flagellum and D. aff. meteor were collected as by-catch from long-line fisheries on R/V Archipelago (ARQDAÇO monitoring program, University of the Azores) and on commercial fishing vessels through a fisheries observer program, during September-November 2017. Collection was performed in various seamounts within the Azorean EEZ (Supplementary Table S1). Colonies were transferred to the aquaria facilities of IMAR (DeepSeaLab) in coolers and distributed in three 170 L aquaria in a thermo-regulated room, in darkness. Corals were inspected for potential tissue injuries from the collection process and colonies with unhealthy tissue were discarded. Aquaria were supplied with seawater pumped from 5 m depth in continuous flow-through open systems. Before entering the aquaria, water was treated with UV-light (P10 UVsystem and Vecton 600, TMC) and was repeatedly filtered (mesh size: 50 μm and 1 μm). Temperature was maintained at 14 ± 0.7 °C, which is similar to the temperature recorded at coral gardens of the two species 81  www.nature.com/scientificreports/ frozen mixture of microalgae, microzooplankton and frozen thawed macrozooplankton (Artemia nauplii, Mysis shrimps), which was enriched with live microalgae and rotifers 2-3 times per week. The collected colonies were left to acclimatize for approximately three months in the aforementioned conditions. Subsequently, colonies were divided in 8-10 cm fragments and mounted to bases made of epoxy one month before the experiments (Fig. 1b,d). During this period, fragments were closely monitored to ensure that they had vibrant colour, intact tissue, and displayed polyp activity.
Feeding experiment. Four different food treatments were created, based on: a phytoplankton derived source (PHYTO), a zooplankton derived source (ZOO), a dissolved organic carbon (DOM) derived source and a fourth treatment where no particulate food source (size > 1 μm) was provided (fasting, FAST). Food treatments were created based on current knowledge of the species biology and on available food sources in their natural environment. The diatom species Chaetoceros calcitrans was selected as a phytoplankton-derived food source. Species of the genus Chaetoceros are common components of spring blooms in some of the sampling sites, e.g. Condor Seamount 81 . The rotifer Branchionus plicatilis was selected as zooplankton-derived food source due to its small size (140-330 μm) and slow swimming capacity, which correspond to the characteristics of zooplankton prey usually captured by octocoral species 22,26 . Due to the known capacity of cnidarians to utilize DOM 47,82 this food source was also used as a food treatment. Lastly, the absence of additional food (FAST) aimed at measuring the basal metabolic activity of the corals, in the absence of particulate food. The feeding experiment was run in four 33 L flumes designed to keep live prey in continuous circulation (Supplementary Fig. S1), which allowed the use of a multilevel experimental design (Fig. 6). The different number of available specimens for the two species led to a slightly different experimental design for each (Fig. 6), however the same rationale was followed for both species. One month before the experiment, each mother colony (n = 5 for D. aff. meteor and n = 12 for V. flagellum) was divided in smaller fragments (n = 5 for each colony of D. aff meteor and n = 4 for each colony of V. flagellum) and these were randomly distributed to the four food treatments. The characteristics of fragments in the different treatments are presented in Supplementary Table S2. Because of the limited number of experimental flumes, we repeated the experimental work several times, in order to have more than one aquaria replicates for each food treatment (n = 3 for D. aff. meteor and n = 4 for V. flagellum). Each repetition is referred to as experimental cycle (Fig. 6). At the beginning of each experimental cycle, coral fragments were randomly positioned in the available aquaria for each treatment. This led to a total of 15 coral fragments for each food treatment for the case of D. aff. meteor and 12 coral fragments for each food treatment for V. flagellum. Dependence among fragments on the colony and aquaria level was treated statistically, by the use of mixed effects models (see statistical analysis).
Before each experimental cycle, zooplankton and phytoplankton food sources were prepared by adding enhanced levels of the stable isotope tracers 13 C and 15 N to the respective culture media. Two microalgae species, C. calcitrans and Nannochloropsis gaditana were cultured using artificial seawater and an F/2 culture medium containing 50% 15 N-sodium nitrate (NaNO 3 , Cambridge Isotopes) and 100% 13 C-bicarbonate (NaHCO 4 , Cambridge Isotopes) for three weeks. Subsequently, cultures were harvested by filtering with membrane filters (0.2 μm), rinsed with filtered SW (0.2 μm) and re-suspended in artificial SW. Rotifer starter cultures (concentration: 45 rotifers ml −1 ) were inoculated in filtered seawater (1 μm), and continuously fed with 13 C 15 N-enriched N. gaditana, cultured as described above, for 6 days. Rotifer cultures were harvested by filtering (nylon filters, 40 μm), rinsed and re-suspended in artificial seawater. Preliminary analysis was performed to ensure that harvesting procedures www.nature.com/scientificreports/ did not affect cell concentration significantly. An algal-derived product of dissolved amino-acids (Cambridge Isotopes, U 13C 97-99%, U 15 N 97-99%, CNLM-452-0.5) was used as DOM food source. Prey was provided to a target carbon (C) concentration of 10 μmol L −1 , similar to POM concentrations previously recorded in mixed gardens of the two species (A. Colaço, pers. comm.). Cultures were scheduled to reach the desired prey concentration, corresponding to the desired target C concentration, on the day of delivery and were harvested a few hours before provision. To monitor the experimental food concentrations, aliquots of the provided food were taken before provision and analyzed a posteriori for DW and carbon content.
For both species, each experimental cycle consisted of five days. At the start of each cycle, fragments were positioned in the aquaria one next to the other, perpendicular to the current, avoiding branch overlapping. Once per day a predefined quantity of food was provided to reach a concentration of 10 μmol C L −1 . Flow of 4 cm s −1 was established for one hour and water renewal was paused just before supplying aquaria with food. This flow speed was selected as it allowed both species to capture live prey (Rakka et al., unpublished data) and kept prey in suspension for 12 h without affecting its concentration. After 12 h, water renewal was reestablished and all remaining food was cleaned by siphoning.
In the last day of each experimental cycle and immediately after the end of feeding, closed cell incubations were performed to measure oxygen consumption, DIC respiration and POC/PON release. Seven coral fragments from each food treatment were transferred to 450 ml glass chambers with pre-filtered seawater (0.7 μm) and glass-coated magnetic stirrers. All chambers were placed in a water bath keeping temperature at 14 ± 0.5 °C. Another two chambers were left without coral fragments and served as controls. Respiration rates were derived from depletion of dissolved O 2 during the incubation, measured by a single channel oxygen meter (Fibox4) with a PSt3 sensor (PreSens, Germany). The chamber with the larger fragment was continuously connected to the sensor to monitor oxygen saturation levels during the incubations. Each incubation lasted for approximately 14 h in which oxygen saturation never dropped below 80%. Oxygen consumption was standardized to the tissue organic carbon content (OC), i.e. without taking into account the main skeletal axis. These values were adjusted for rates recorded in chambers without coral fragments to account for microbial respiration. Water samples were taken before and after each incubation to determine the concentration of DIC and 13 C-DIC. Samples were kept in 10 mL headspace vials with 10 μL of a saturated mercury chloride solution and stored at 4 °C until analysis. Lastly, the remaining water (300 ml) from each chamber after the end of the incubation was filtered through precombusted, preweighted GF/F (0.7 μm, WHATMAN) filters to estimate POC, 13 C-POC, PON and 15 N-PON release. Filters were freeze-dried and kept at room temperature until analysis.
Upon completion of each experimental cycle, fragments were freeze-dried and stored at -80 ºC. Fragments were dissected to separate the tissue from the skeleton. The tissue was ground by mortar and pestle, and a subsample was analyzed for total C and N content and isotopic ratios using an elemental analyzer (Thermo Electron Flash 1112) coupled to an isotope ratio mass spectrometer (EA-IRMS, DELTA-V, THERMO Electron Corporation). A second subsample was acidified stepwise with drops of HCl to remove the inorganic C fraction, and all remaining material was analyzed on the elemental analyzer for organic C content and isotopic ratio.
Calculation of tracer C and N incorporation, i.e. calculation of the amount of C and N which came from the provided food source and was incorporated in coral tissue, was performed as described in Maier et al. 63 . Tissue C and N content of each fragment (i.e. without coral skeleton) was standardized to DW, and expressed as mmol C or N (g DW) −1 . The heavy/light isotope ratio (e.g. 13 C: 12 C) of each coral fragment (Rsample) was calculated as Rsample = ([δtracerCsample/1000] + 1) × Rref, where Rref = 0.0111802 for organic C (OC) and Rref = RN 2 = 0.0036782 for organic N (ON). Fractional abundance of 13 C and 15 N (e.g. F 13 = 13 C/[ 12 C + 13 C]) was expressed as Ftracer = Rsample/(Rsample + 1). Experimental 13 C and 15 N enrichment of each coral fragment tissue was expressed in relation to the fractional abundance of the respective fasting (non-enriched) fragment or average of fasting fragments of the same colony. Finally, tracer 13 C incorporation was calculated by multiplying 13 C enrichment with tissue OC content (μmol 13 C fragment −1 ) and tracer 15 N incorporation was obtained by multiplying 15 N enrichment with tissue ON content (μmol 15 N fragment −1 ). The total amount of C or N incorporated into coral tissue from the provided labelled food source (tracer C and N incorporation) was calculated by dividing the tracer C or N incorporation of each fragment with the fractional abundance (F 13 or F 15 ) of the respective food source. Final tracer C and N incorporation rates were normalized to the OC (mmol) of each coral fragment.
Concentration of DIC was determined on an Apollo SciTech AS-C3 analyzer, after transforming DIC to gaseous carbon dioxide by addition of concentrated phosphoric acid (H 3 PO 4 , volume: 10 μL mL −1 ) in each headspace vial. Subsequently, a 10 μL subsample of headspace gas was obtained from each vial and analyzed in the isotope ratio mass spectrometer, as described above, to obtain measurements of δ 13 C. Filters collected for POC and PON measurements were weighted (accuracy 0.1 mg) and analyzed with the isotope ratio mass spectrometer to obtain concentrations of POC and PON, as well as δ 13 C of and δ 15 N respectively. Determination of fluxes, i.e. DIC respiration, POC and PON release was estimated in two steps. Firstly, the bulk fluxes, i.e. the amount of total C and N released during the incubation period were calculated as the respective concentration difference between start and end water sample. Subsequently, tracer fluxes, i.e. the amount of C and N derived from the provided food were estimated, by multiplying bulk fluxes by their relative enrichment in 13 C and 15 N during the incubation, and dividing by the food enrichment following Maier et al. 83 . A final tracer C budget was compiled by estimating tracer C incorporation, tracer C respiration and tracer C release for the duration of the whole experiment for each treatment and is reported as percentage of the provided C. Statistical analysis. Data exploration was done following Zuur et al. 84 to select the most appropriate statistical modelling method. To test if independent factors had a significant effect on the dependent variables in question, the former were added progressively to the models and the Akaike Information Criterion (AIC) along with maximum likelihood ratio (MLR) tests were used to select the most appropriate model. www.nature.com/scientificreports/ were inspected to detect potential violation of model assumptions. Statistical analysis was performed in R 3.5.0 (R Core Team, 2018). We provide detailed results of the MLR tests in the supplementary information (Supplementary Table S3) and coefficients of the best models in Table 2. Linear Mixed Effects Models (LMEs) and GLSs were used to analyze all response variables. Colony and experimental aquaria were incorporated as crossed random factors to deal with dependence related with: (1) the existence of multiple coral fragments that originated from the same colony and (2) the fact that multiple coral fragments were positioned in the same aquaria. Whenever the assumption of homogeneity of variance was not fulfilled, variance structure components were added to the models to allow the variance to differ among tested treatments 85 . LME models were build using the packages LME4 86 and nlme 87 .

Data availability
The datasets generated and analyzed in the current study are available in the Pangaea repository, under the following link: https:// doi. panga ea. de/ 10. 1594/ PANGA EA. 913184.