Introduction

Octocorals are common benthic suspension feeders in tropical, subtropical, temperate and polar regions1,2,3. The majority of octocoral species are found in waters deeper than 50 m4 where they create dense single- or multi-species aggregations, structuring three-dimensional and highly heterogenous habitats known as coral gardens3,5. These communities provide essential habitat for a variety of associated fauna3,6.

The Azores Archipelago, located in the central North Atlantic, harbors an extremely rich biodiversity of cold-water octocorals, reaching a total of 101 species which represent the highest octocoral species richness known so far in North Atlantic7,8. Coral gardens constitute the most prominent cold-water coral (CWC) habitat in the Azores, with monospecific or multispecific octocoral communities frequently colonizing seamounts and island slopes7,9. Because of their life-history traits, including slow growth and high longevity, recovery of octocoral communities from fisheries and other disturbances can be very slow10,11 and thus coral gardens have been classified as vulnerable marine ecosystems (VMEs) in need of protection12,13. However, effective conservation of VMEs requires knowledge on the species biology and ecology which in the case of deep-sea octocorals is scarce1.

Resource acquisition is a key factor in the biology of suspension feeders14,15, ultimately determining population dynamics and species distributions16,17,18. Thus, knowledge on feeding biology of key habitat formers such as octocorals is pivotal to understand local ecosystems. Octocorals can feed on a variety of prey including microplankton, nanoeukaryots, as well as detritus19,20,21,22. In some cases, their diet varies seasonally following the cycles of local phytoplankton and zooplankton communities15,23. Although most octocoral species seem to be able to ingest and utilize phytoplankton20,23,24, small zooplankton with low mobility is the main component of the natural diet of many temperate species22,25,26. While considerable knowledge exists on the feeding biology and ecophysiology of shallow octocorals, such information is scarcer for deep-sea octocorals with a few studies so far focusing mainly on Antarctic ecosystems20,27.

In this study, we examined the feeding biology of two common habitat-forming deep octocoral species in the Azores Archipelago: Dentomuricea aff. meteor and Viminella flagellum. The two species form dense coral gardens (Fig. 1) on seamounts between 200–600 m and very frequently occur in mixed populations7,9. The objectives of the study were (1) to determine the ability of the two species to exploit different food sources and (2) to evaluate if assimilation of different food sources affects their metabolic activity. We employed an experimental approach, with the use of aquaria flumes with steady flow velocity and four different food treatments including provision of live phytoplankton, dissolved organic matter (DOM) and live zooplankton, as well as fasting (deprivation of particulate food). Food sources were artificially enriched with 13C and 15N and food utilization was quantified (a) as the appearance of 13C15N in the coral tissue, indicating tracer carbon (C) and nitrogen (N) incorporation from the provided food, (b) as the production of 13C-enriched dissolved inorganic carbon (DIC) by the coral, indicating tracer C respiration, and (c) as the production of 13C15N-enriched particulate organic carbon and nitrogen (POC, PON) by the coral, indicating tracer POC and PON release.

Figure 1
figure 1

(a) Mixed coral garden of the octocorals Viminella flagellum and Dentomuricea aff. meteor (Gavin Newman, Greenpeace); coral fragment of D. aff. meteor (b) and its polyps (c); coral fragment of V. flagellum (d) and its polyps (e). Scale bar 1 mm.

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).

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 (F13, F15) and carbon to nitrogen (C/N) ratio. CC: Chaetoceros calcitrans, NG: Nannochloropsis gaditana, BP: Branchionus plicatilis.

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.

Figure 2
figure 2

Tracer C (a) and N (b) incorporation (average ± SD) of the octocoral species Dentomuricea aff. meteor and Viminella flagellum upon provision of different food sources enriched with 13C and 15N. Axis breaks are used to highlight the large differences of tracer among treatments. PHYTO: phytoplankton; Chaetoceros calcitrans; DOM: dissolved organic matter; ZOO: zooplankton Branchionus plicatilis.

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.

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.

Figure 3
figure 3

Oxygen consumption (average ± SD) of the octocoral species Dentomuricea aff. meteor and Viminella flagellum upon provision of different food sources. Coral fragments were fed with the respective food source for four days and oxygen was measured in closed-cell incubations that took place immediately after feeding on day four and lasted for approximately 12-14 h. FAST: no food provision; PHYTO: phytoplankton Chaetoceros calcitrans; ZOO: zooplankton Branchionus plicatilis; DOM: dissolved organic matter.

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).

Figure 4
figure 4

Tracer fluxes (average ± SD), including tracer C respiration (a), tracer C release (b) and tracer N release (c) of the octocoral species Dentomuricea aff. meteor and Viminella flagellum upon provision of different food sources. Numbers below bars represent the number of coral fragments for which positive estimates were obtained (max 7). Axis breaks are used to highlight large differences in scale among some treatments. PHYTO: phytoplankton Chaetoceros calcitrans; DOM: dissolved organic matter; ZOO: zooplankton Branchionus plicatilis.

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 the provided C was traced back, reaching 23% under the ZOO treatment, 0.48% under the DOM treatment and 0.17% under the PHYTO treatment.

Figure 5
figure 5

Tracer utilized by fragments of two octocoral species Dentomuricea aff. meteor and Viminella flagellum, expressed as a percentage of the provided carbon of different food sources: PHYTO: phytoplankton Chaetoceros calcitrans; DOM: dissolved organic matter; ZOO: zooplankton 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 scleractinians28,29,30,31. However, only a few deep-sea octocorals, mainly of the family Primnoidae, have been shown to base their diet on microzooplankton32,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 waves35,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 herbivory24,37,38. Moreover, in previous studies, both target octocorals have been shown to occupy low trophic levels, placed between primary and secondary consumers39. Thus, the lower incorporation 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 pertusum40) which displayed a rather unselective feeding behavior under feeding with similar food sources and flow velocities41. 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 octocorals42,43 and scleractinians44,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 cycles46, 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 scarce47. 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 change48. A number of previous studies have highlighted the variable and seasonally-dependent diet of octocoral species in temperate ecosystems19,49,50 while similar seasonality has been also reported in benthic Antarctic ecosystems27,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 sources31. 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 rises16,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 availability46, 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 pathogens54. 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 productivity54,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.

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 process59,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 pattern61. 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 costs14,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. pertusa63. 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) species18,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 assemblages64. 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 species65,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 communities67. 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 webs68,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 level70. Species distribution modelling has shed light to the distribution of different deep-sea coral groups, including octocorals71, but comparatively little effort has been made to delve further into the biological and physiological characteristics which shape these distributions72,73. Taking into account these traits will not only improve predictions on species occurrences and help to identify priority areas for conservation and management74, 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 seafloor48. 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 depth7. It can reach heights of up to 1.5 m and can create dense monospecific or mixed populations with other species (Fig. 1a).

Viminella flagellum75 is a whip coral of the family Ellisellidae. It creates monopodial colonies without branches which can grow up to 3 m height76. Its distribution includes the eastern North Atlantic coast, islands of the Macaronesia and the Mediterranean Sea77,78,79. It is usually encountered between 120 and 500 m depth77,78 and can form monospecific or mixed coral aggregations, often with branching octocorals such as D. aff. meteor, Callogorgia verticillata and Acanthogorgia armata80.

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, by cooling systems connected to temperature controllers. Corals were fed daily with a 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 Seamount81. 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 species22,26. Due to the known capacity of cnidarians to utilize DOM47,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).

Figure 6
figure 6

Experimental design of the two feeding experiments with the octocoral species Dentomuricea aff. meteor and Viminella flagellum. Exp cycle: Experimental cycles; Aq: Aquaria; PHYTO: phytoplankton Chaetoceros calcitrans; ZOO: zooplankton Branchionus plicatilis; DOM: dissolved organic matter; FAST: no food provision. Rectangles represent experimental aquaria and black dots represent coral fragments.

Before each experimental cycle, zooplankton and phytoplankton food sources were prepared by adding enhanced levels of the stable isotope tracers 13C and 15N 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% 15N-sodium nitrate (NaNO3, Cambridge Isotopes) and 100% 13C-bicarbonate (NaHCO4, 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 13C15N-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 did not affect cell concentration significantly. An algal-derived product of dissolved amino-acids (Cambridge Isotopes, U 13C 97–99%, U 15N 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 O2 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 13C-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, 13C-POC, PON and 15N-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. 13C: 12C) of each coral fragment (Rsample) was calculated as Rsample = ([δtracerCsample/1000] + 1) × Rref, where Rref = 0.0111802 for organic C (OC) and Rref = RN2 = 0.0036782 for organic N (ON). Fractional abundance of 13C and 15N (e.g. F13 = 13C/[12C + 13C]) was expressed as Ftracer = Rsample/(Rsample + 1). Experimental 13C and 15N 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 13C incorporation was calculated by multiplying 13C enrichment with tissue OC content (μmol 13C fragment−1) and tracer 15N incorporation was obtained by multiplying 15N enrichment with tissue ON content (μmol 15N 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 (F13 or F15) 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 (H3PO4, 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 δ13C. 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 δ13C of and δ15N 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 13C and 15N 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. Model diagnostics 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 treatments85. LME models were build using the packages LME486 and nlme87.

Data availability

The datasets generated and analyzed in the current study are available in the Pangaea repository, under the following link: https://doi.pangaea.de/10.1594/PANGAEA.913184.