- Open Access
Nutrient dynamics in coral symbiosis depend on both the relative and absolute abundance of Symbiodiniaceae species
Microbiome volume 10, Article number: 192 (2022)
Symbionts provide a variety of reproductive, nutritional, and defensive resources to their hosts, but those resources can vary depending on symbiont community composition. As genetic techniques open our eyes to the breadth of symbiont diversity within myriad microbiomes, symbiosis research has begun to consider what ecological mechanisms affect the identity and relative abundance of symbiont species and how this community structure impacts resource exchange among partners. Here, we manipulated the in hospite density and relative ratio of two species of coral endosymbionts (Symbiodinium microadriaticum and Breviolum minutum) and used stable isotope enrichment to trace nutrient exchange with the host, Briareum asbestinum.
The patterns of uptake and translocation of carbon and nitrogen varied with both density and ratio of symbionts. Once a density threshold was reached, carbon acquisition decreased with increasing proportions of S. microadriaticum. In hosts dominated by B. minutum, nitrogen uptake was density independent and intermediate. Conversely, for those corals dominated by S. microadriaticum, nitrogen uptake decreased as densities increased, and as a result, these hosts had the overall highest (at low density) and lowest (at high density) nitrogen enrichment.
Our findings show that the uptake and sharing of nutrients was strongly dependent on both the density of symbionts within the host, as well as which symbiont species was dominant. Together, these complex interactive effects suggest that host regulation and the repression of in hospite symbiont competition can ultimately lead to a more productive mutualism.
Symbiosis is a source of evolutionary innovation in which an entire ecological/biological repertoire can be acquired through close association with other species. These partnerships yield trait combinations that can expand access to nutrients, provide protection from enemies, and/or aid in mobility and dispersal. In corals, endosymbiotic dinoflagellate algae of the family Symbiodiniaceae supply the majority of the host’s energetic needs via photosynthesis. Symbiodiniaceae benefit from access to a light-rich habitat and inorganic nitrogen sources in the form of the host’s metabolic waste. This combination of trophic strategies (i.e., algal photosynthesis and coral heterotrophy) and efficient recycling of nutrients between the partners underpins the success of corals in nutrient-limited tropical oceans. However, anthropogenic stressors have exposed the vulnerability of this critical symbiosis to environmental change resulting in global coral bleaching and mortality events [1, 2], but flexibility within this symbiosis may prove to be its salvation .
The identity and relative abundance of algal symbionts define many aspects of coral physiology, including growth [4,5,6], disease , and bleaching . Comparative studies of corals and other cnidarians hosting one Symbiodiniaceae species or another have demonstrated differences among symbiont species in their contributions to host metabolism, with variations in both the amount of carbon [9,10,11,12] and assimilated inorganic nitrogen  transferred to the coral. As genetic techniques open our eyes to the breadth of symbiont diversity within myriad microbiomes, our understanding of the ecological mechanisms that restrict and maintain diversity in symbiosis has not kept pace . Most corals maintain some capacity for symbiont flexibility , and diversity of symbionts at background abundances is common . Yet it remains poorly understood how ecological interactions among co-occurring symbionts, and between each symbiont and their host, will influence the balance of costs and benefits of the association.
Natural selection on symbionts alone is not expected to select for beneficial partners, instead, as in all ecological systems, selection should favor symbionts best able to exploit their host and outcompete other symbionts . As such, unregulated symbiont diversity not only exposes the host to sub-optimal symbionts but also to the consequences of antagonism among symbionts competing for space within the host . Recent work demonstrated that symbionts alter their metabolism in vast ways when a competing symbiont species is introduced in culture, shifting investment of newly assimilated nutrients between compounds that support growth vs. storage . If those competitive behaviors persist for symbionts within coral tissues, they may negatively affect the quality and persistence of the symbiosis. The result of competition may prevent the coral from associating with the most beneficial symbionts, divert nutrients from the coral-algal exchange, and/or potentially expose the host to allelopathic compounds. Although multiple symbiont cells can reside within a single host cell, individual symbionts are bound by a multilayered symbiosome membrane that serves as the bridge for transport of nutrients, gases, and photosynthetically fixed carbon. This membrane helps to maintain metabolic stability within symbiosis  and may also inhibit the ability of symbionts to detect the presence of potential competitors.
It is critical for the stability of the mutualism that corals and other hosts develop mechanisms to both avoid conflict among symbionts and to avoid associations with less beneficial symbionts. To ensure the highest return on their investment, many hosts have evolved mechanisms to allocate more resources to highly cooperative symbionts and/or to inhibit associations with more parasitic symbionts. These host sanctions and rewards are two mechanisms by which symbiotic interactions can be stabilized on ecological timescales . For example, soybeans can reduce growth of intracellular rhizobia that fail to fix nitrogen , while other legumes can reward specific fungal partners with more carbohydrates . Corals have some capacity to regulate nutrient availability to their symbionts to inhibit over-proliferation of less beneficial symbionts in hospite [23, 24] and ensure the release of sufficient photosynthetic carbon to the host . But the ability of corals to use these mechanisms to regulate community composition remains unknown.
In most corals, symbiosis is initiated via exogeneous uptake of symbionts in the early stages of ontogeny. Over time, what can initially be a diverse community winnows towards a relatively predictable, stable association [26, 27]. The octocoral Briareum asbestinum provides a well-studied model for symbiont dynamics in an obligate mutualism, in part because of its ability to associate with several symbiont species, including but not limited to Symbiodinium microadriaticum and Breviolum minutum. Furthermore, B. asbestinum juveniles have documented patterns of symbiont uptake and succession in both the field and the laboratory [26, 28], with outcomes that indicate that host control , symbiont competition , and environmental context  play a role in shaping the symbiont community. Here, we use stable isotope tracers of carbon and nitrogen to track the nutritional currency exchange between host and symbionts as a function of symbiont community diversity to examine how nutrient sharing in the symbiosis is affected by symbiont diversity and better understand the ecological processes that shape this symbiosis.
Briareum asbestinum is a common Caribbean octocoral that releases aposymbiotic larvae that acquire symbionts from the environment in the following weeks-years . In May of 2013, larvae were collected from the surfaces of approximately 30 B. asbestinum colonies on a shallow reef (~2 m depth) along Long Key, Florida, and transported back to the Keys Marine Laboratories (KML; Long Key, FL, USA). Once in the laboratory, larvae were rinsed and maintained for the duration of the experiment in artificial seawater (ASW; Instant Ocean®; salinity: 36 ppt) to prohibit exposure to environmental sources of Symbiodiniaceae. For the first 48 h, larvae were maintained in plastic containers with noncirculating but regularly refreshed ASW. The larvae were transferred to 1.5-L sterile, plastic water bottles (~300 larvae L−1) and transported to the University of Buffalo culture facilities where they settled onto cleaned gorgonian axial branches, a common natural settlement substrate, which had been collected in the field and sterilized. Once larvae had settled and metamorphosed into polyps (~2 weeks following collection), hereafter referred to as “recruits,” the branches were distributed among 1-L containers with 600 mL ASW (n = 32 containers with 30–40 B. asbestinum recruits each).
Two cultured Symbiodiniaceae species were selected from the BURR Culture Collection at the University at Buffalo: Symbiodinium microadriaticum (Culture ID: 04-503; cpType: A194, ; ITStype: A1 ;) and Breviolum minutum (Culture ID: Mf01.05b; cpType: B184; ITStype: B1). Each of these species has the ability to establish sustained symbioses with B. asbestinum in previous laboratory studies [28, 30]. They have also been found in naturally established B. asbestinum juveniles in the field but are generally displaced by other Symbiodiniaceae species within months following settlement [26, 33]. However, those homologous symbionts have yet to be successfully brought into culture.
Previous labwork demonstrated that S. microadriaticum is displaced from the symbiosis with B. asbestinum following exposure to Breviolum species including B. minutum [28, 30]. Therefore, we varied the timing of initial exposure to B. minutum (months to weeks before the isotope tracer experiment) as well as the relative ratio of symbiont species inocula (S. microadriaticum:B. minutum target ratios — 100:0, 90:10, 75:25, 50:50, 25:75, 10:90, 0:100) to generate recruits with an array of single and mixed symbiont species associations. The actual abundance and ratio of symbiont species were determined for each individual recruit using qPCR (see below). For 3 months, recruits were maintained in artificial seawater (ASW; Instant Ocean®, Blacksburg, VA, USA) at 28 °C, beneath four fluorescent T5 grow lights (measured at 165–190 photons m−2 s−1 with a LI-COR optical sensor) on a 14/10 h light:dark cycle. Throughout this time, regular (2–3 times per week) maintenance included full ASW changes and symbiont inoculations at total densities of 200 cells/ml. Recruits were also fed weekly with Zeigler® Larval AP100 Diet (Zeigler Bros., Inc.) during the first 2 months of maintenance, but feeding was discontinued 1 month prior to the isotope tracer experiment.
Pulse chase with dual isotope tracers
A pulse-chase experimental design was used to trace the assimilation and transfer of nitrate (15NO3) and bicarbonate (H13CO3) in B. asbestinum recruits (Fig. 1A). First, B. asbestinum were acclimated for 2 days in ASW supplemented with non-isotopically enriched nitrate and bicarbonate to stimulate the induction of nitrate-reducing enzymes . They were then transferred to isotopically enriched ASW for a period of 12 h of light (pulse phase). Finally, they were rinsed under flowing ASW to remove unincorporated compounds from their surfaces and transferred to ASW for 24 h (chase period). Pilot studies were used to determine the appropriate duration of the pulse and chase period and to determine the amount of tissue necessary for accurate nutrient readings.
Each experimental chamber consisted of a 1-l glass beaker with 0.7 L of stirring seawater. Within this, we added branch(es) with 30–40 B. asbestinum recruits. A 17-l stock solution of media for each stage of the experiment (acclimation, pulse, chase) was prepared at the start and used for all replicates. Nutrient supplementation at [NO3−] = 10 μM and [HCO3−] = 0.4 mM was included in the acclimation media, using non-isotopically enriched compounds, and in the pulse media using isotopically enriched compounds. No supplemental nutrients were added in chase media (Table S1). Laboratory lighting was maintained at (165–190 photons m−2 s−1; Licor) on a 14:10 h light:dark cycle. Four chambers were run in parallel, and the entire acclimation/pulse/chase routine was repeated five times over 2 weeks for a total of 20 experimental replicates. As a control, we included a set of dark incubations which were covered by a black-out box for the entirety of the pulse period to inhibit photosynthesis and confirm the symbiont as primary source of traced nutrient incorporation (Fig. S1). Water samples were taken from each of the stocks throughout the experiment for analysis of NO3-N, NH3-N, TDN, and DIC at the Penn State Institutes of Energy and the Environment (Table S1).
At the end of the experiment, B. asbestinum were removed at their base from branches, rinsed twice in 0.2-μm filtered DI water, transferred into 600 μl of 0.2-μm filtered DI water, and macerated with a Tissue Master 125 (Omni International). One-hundred microliters of homogenized tissues was subsampled, preserved in 90% EtOH, and stored at −20 °C for subsequent genetic analysis. The remaining homogenate was separated into host and symbiont fractions by a series of centrifugation and washing steps (Fig. 1B). First symbionts were pelleted at 1000 rcf, and the supernatant with host tissue was removed and frozen. The pellet was then rinsed twice in DI water and frozen. A total of 630 B. asbestinum recruits were processed.
Genomic DNA was extracted from homogenate tissue subsamples using a CTAB protocol optimized for Symbiodiniaceae extraction . Quantitative PCR (qPCR) using genera-specific primers and probe sets were used to quantify the total copies of each of the Symbiodinium and Breviolum actin genes . To convert actin gene abundance to absolute cell numbers, and account for interspecies differences in actin copy number, standards made from genomic DNA extracted from cultured S. microadriaticum and B. minutum were also prepared from 2000, 8000, 16,000, 64,000, 19,2000, and 384,000 cells.
Reaction volumes were 10 μl with 5 μl TaqMan Genotyping Master Mix and 1 μl genomic DNA template. Assays were optimized for each target including the following: Symbiodinium — 150 nM forward primer, 100 nM reverse primer, and 150 nM probe and Breviolum — 200 nM forward primer, 300 nM of reverse primer, and 100 nM probe. All qPCR reactions were performed using a CFX Real-Time PCR Machine (Biorad) with 384-well module with an initial incubation (2 min at 50 °C, 10 min at 90 °C) followed by 40 cycles of 10 s at 95 °C and 1 min at 60 °C. Cycle threshold (CT) values for each assay were calculated using an automatic baseline interval and relative fluorescence threshold of 0.01. Each sample was run in duplicate, and each plate included a full set of replicated standards for each symbiont species. Spearman’s rank correlation was used to test for a correlation between symbiont density (the number of cells per recruit) and symbiont ratio (the relative ratio of S. microadriaticum to B. minutum).
Stable isotope analysis (SIA)
We determined that 8–10 B. asbestinum recruits were required to obtain enough host and symbiont biomass for accurate determination of C and N isotopic ratios. We used the symbiont ratio as a basis to combine 8–10 paired host or symbiont tissues from individual recruits onto a filter for SIA analyses (Fig. 1B). Samples were prepared for SIA by filtering onto precombusted 0.7 μm glass fiber filters (Whatman GF/F), rinsed once with 1% HCl to remove inorganic carbonate from residual media, rinsed three times with DI water, and then folded, wrapped in foil, frozen, and transferred to SUNY-ESF.
At the Environmental Science Stable Isotope Laboratory (EaSSIL) at SUNY-ESF, samples were lyophilized (−60 °C, 48 h, 150 torr) and packed into tin capsules. Using a Costech elemental analyzer linked via a ThermoFinnigan ConFlo III interface to a Finnigan MAT Delta XL Plus stable isotope mass spectrometer (EA-IRMS), we determined total organic carbon content, total organic nitrogen content, δ13C, and δ15N values of host, symbionts, and nutrient solutions. The precision of the stable isotope measurements (expressed in the standard per mil notation) was verified using the National Institutes of Standards and Technology RM8573 (δ13C = −26.4 ± 0.1‰, δ15N = −4.5 ± 0.3‰, n = 38), and RM8574 (δ13C = +37.3 ± 0.3‰, δ15N = +47.6 ± 0.3‰, n = 38). Daily precision of the instrument was verified by repeated analyses of internal laboratory standards including acetanilide (δ13C = −33.4 ± 0.2‰, δ15N = −1.0 ± 0.1‰, n = 32), valine (δ13C = −10.7 ± 0.2‰, δ15N = −6.6 ± 0.1‰, n = 9), and daphnia (δ13C = −24.4 ± 0.2‰, δ15N = +17.4 ± 0.3‰, n = 8), during the sample runs.
Enriched values are reported as atom percent of the heavy isotope (AP15N & AP13C), which is calculated by the count of the heavy isotope relative to the total number of atoms of N or C in the sample. We also calculated the total milligrams of newly assimilated nitrogen per recruit (Nnew), by multiplying AP15N (per mille) by the total organic nitrogen (mg) for each of the tissue fractions and summing them with the assumption that N is highly retained in the symbiosis through the duration of the experiment, whereas carbon may be respired over time .
Multiple linear regression was used to test if symbiont density (number of symbionts per recruit), symbiont community composition (symbiont ratio = the relative ratio of S. microadriaticum and B. minutum cells [S. m.:B. m.] determined by qPCR), or their interaction were significant predictors of the response variables AP13Chost, AP13Csym AP15Nhost, AP15Nsym, and Total Nnew/recruit. When a significant interaction term was identified, we used the interflex package in R (Hainmuller et al. 2021) to generate marginal effect plots with 95% confidence intervals and visualize how the effect of one predictor variable changed with each level of the other predictor variable. Using the confidence intervals, we approximated where the slope was significantly different from 0, partitioned the data accordingly, and plotted the data with the significant and nonsignificant linear relationships of partitioned groups. All data were confirmed to meet the assumptions of homogeneity of variance and linearity, with the residuals approximately normally distributed.
Symbiont density and relative abundance
Single-strain inoculations yielded ~50 recruits each hosting either S. microadriaticum or B. minutum. Of these, 20 recruits with only S. microadriaticum and 20 with only B. minutum were used as dark controls and the rest in the main experiment. The regime of varying the timing of exposure and ratio of S. microadriaticum and B. minutum used for mixed inoculations was successful in establishing ~350 recruits of mixed symbiont species assemblages at various symbiont ratios (S. m.:B. m.) (Fig. 2A). However, an artifact of delaying exposure to B. minutum was an overrepresentation of recruits with higher ratios of S. microadriaticum, including an additional 25 recruits in which only S. microadriaticum was detected. This occurred most often in treatments in which B. minutum exposure occurred only in the final weeks preceding the pulse-chase experiment. Mean symbiont densities across all recruits ranged from 268,738 to 1,387,584 cells/recruit, and the density of symbionts was not correlated with symbiont ratio (Spearman’s rank correlation rho(584) = 0.01, p = 0.83; Fig. 2B). Therefore, we grouped 8–10 recruits onto filters based on their similarity in S. m.:B. m. This resulted in 59 pairs of host tissue and symbiont tissue filters on which we performed stable isotope analysis. The maximum range of S. m.:B. m. values among individual recruits combined onto a single filter was 4.7% with a mean range of S. m.:B. m. of 2% (Fig. 2A).
A high proportion and abundance of B. minutum lead to higher carbon assimilation rates in both host and symbiont tissues. In host tissues, both an increasing proportion of B. minutum (p < 0.001) and, to a lesser extent, an increasing density of symbionts overall (p = 0.021) were positively correlated with AP13Chost (F(3.44) = 28.6, p < 0.001) with the model explaining 54% of the variance (Fig. 3A and B). However, in symbiont tissues, there was a significant interaction between symbiont ratio and symbiont density on AP13Csym (p = 0.035). In this case, the model explained 17% of the variance (F(3.44) = 4.26, p < 0.01) increasing to 37% of the variance explained where two identified outliers were removed (F(3.42) = 9.86, p < 0.01). Marginal effect plots indicated that as B. minutum increased in proportion, so did carbon AP13Csym, but that this only occurred when symbiont densities were higher than ~700,000 cells per recruit (Fig. 3C and D, Fig. S2).
Symbiont ratio and density were shown to be interactive in predicting the amount of newly assimilated nitrogen in both host and symbiont tissues. In host tissues, the full interactive model significantly (F(3.43) = 3.25, p = 0.03) explained 17% of the variance in measured AP15Nhost with a significant effect of symbiont density (p < 0.01) and a significant interaction term (p = 0.04; Fig. 4A and B). Marginal effect plots indicate that the overall trend in nitrogen enrichment varied depending on the dominant symbiont species (Fig. S3). In recruits with > 50% B. minutum, AP15Nhost did not scale with symbiont densities. However, in recruits with > 50% S. microadriaticum, increasing densities of symbionts per host led to higher values of AP15Nhost.
For symbiont tissues, significant model terms included symbiont density (p < 0.001) symbiont ratio (p < 0.001) and their interaction (p = 0.002) explaining 31% of the variance (F(3.44) = 8.16, p < 0.001; Fig. 4C and D). Marginal effect plots indicate that within less densely populated hosts, i.e., those harboring less than ~800,000 symbionts per host, AP15Nsym was lower in hosts dominated by B. minutum and increased with increasing proportions of S. microadriaticum (Fig. S4). Concomitantly, within those hosts dominated by S. microadriaticum, increasing symbiont densities resulted in significant decreases in AP15Nsym (Fig. S5).
We examined the relationship between host and symbiont tissue enrichment by plotting AP13Csym by AP13Chost and modeling that relationship for S. microadriaticum-dominated and B. minutum-dominated hosts (Fig. 5A). Both showed a positive correlation between symbiont and host enrichment, but the slope of the relationship was higher in B. minutum-dominated hosts than those dominated by S. microadriaticum (0.83 and 0.28, respectively; Fig. 5A). Similar trends were seen for the relationship between AP15Nsym and AP15Nhost dominated by either B. minutum or S. microadriaticum (1.07 and 0.62, respectively; Fig. 5B).
The sum of the newly assimilated nitrogen will significantly be affected by symbiont density (p = 0.03) with a significant interaction between symbiont density and ratio (p = 0.02); this full interaction model explained 20% of the variance (F(3.41) = 4.61, p < 0.01; Fig. 6). Marginal effect plots indicate that in hosts dominated by S. microadriaticum, increasing symbiont densities was negatively correlated with Total Nnew (Fig. S6).
Recent research has highlighted the complex and principal role of nutrient exchange among coral hosts and their algal symbionts in maintaining a functional symbiosis [6, 23, 38]. In this study, we found that nutrients were largely affected by interactive effects of total symbiont density and the relative ratio of symbiont species. For example, only after a symbiont density threshold was achieved was carbon shown to scale with the relative ratio of symbionts, increasing with the relative ratio of B. minutum. Conversely, nitrogen dynamics were characterized by the dominant symbiont species. Those hosts dominated by B. minutum received consistent nitrogen resources independent of symbiont density, whereas increasing symbiont densities negatively impacted nitrogen assimilation in hosts dominated by S. microadriaticum. But rather than evidence for interspecific competition among symbionts, these species-specific patterns in nutrient cycling provide a mechanism by which nutrient dynamics may be involved in host regulation of symbiont community structure.
Recent study of Symbiodiniaceae species in co-culture demonstrated that interspecific competition altered nutrient assimilation rates and subsequent compound production of Cladocopium goreaui and Durusdinium trenchii . If these interactions are sustained within the host, where nitrogen is limited , they have the potential to destabilize the symbiosis through alterations in both nutrient acquisition and sharing. Instead, in hospite cohabitation of S. microadriaticum and B. minutum within the tissues of host B. asbestinum did not show a competitive shift in metabolism. After a threshold density of symbionts was reached, the production and sharing of photosynthetically derived carbon were well predicted by the ratio of each symbiont type with mixed communities providing intermediate levels of carbon resources (Fig. 3A and C). An increased number of co-dominated (S. m.:B. m. ~1:1) recruits may provide a more clear and potentially nonlinear fit to this relationship; however, the observed effect remains far smaller than that expected from co-cultured Symbiodiniaceae experiments (McIlroy et al. 2020). While multiple symbiont cells can exist within a single host cell , the symbiosome may serve to physically isolate individual symbionts . The additive relationship observed here suggests that interactions among in hospite symbionts are indeed restricted. A similar pattern was seen in radiotracer experiments in the giant sea anemone where significant differences in carbon translocation (14C) occurred between anemones harboring “type A” or “type B” Symbiodiniaceae, while mixed symbiont assemblages (A + B) translocated intermediate levels of carbon relative to the two monophyletic groups . Both our study and Loram et al. (2007) assessed isotopic values of bulk tissues of mixed symbiont communities. It is therefore possible that antipodal responses of symbiont species would have gone undetected . Ultimately however, interactions among codominant symbionts, if any, did not interfere with the balance of carbon in the symbiosis more generally. Carbon translocation (AP13Chost) mirrored that of symbiont tissues (AP13Csym) with carbon enrichment values increasing linearly with increasing proportions of in hospite B. minutum to S. microadriaticum symbionts (Fig. 3A and C).
We found that S. microadriaticum was not only less productive but also more selfish with carbon resources (Fig. 5A), as has been reported in other hosts [42, 43]. In Briareum asbestinum, hosts that harbor B. minutum have been shown to have higher survival rates than those that harbour other symbiont species . This link between variation in symbiont productivity and host survival provides a basis for natural selection in favour of hosts that can limit long-term associations with S. microadriaticum and/or promote associations with B. minutum. Indeed, in both the lab and field, B. asbestinum transition to Breviolum-dominated symbiont communities [26, 28]. This occurs in spite of the fact that S. microadriaticum can infect corals sooner and at higher densities than other symbiont species . While other studies have suggested that increasing symbiont densities can cause self-shading, and limit photosynthetically driven carbon assimilation , this was not seen in our study. Instead, symbiont densities did not affect AP13Csym and led to slight but significant increases in AP13Chost (Fig. 3). Symbiodiniaceae are capable of maintaining high rates of photosynthesis across light conditions by increasing the efficiency of light harvesting and utilization [44, 45]. Hosts with higher symbiont densities may also represent later stages in the onset of the symbiosis wherein symbiont cells drastically reduce cell replication and are able to generate more photosynthates in excess [6, 36].
Although carbon fixation via photosynthesis is essential to the coral-algal symbiosis, in oligotrophic tropical waters, nitrogen availability often limits reef productivity (reviewed in [46, 47]) and has a primary role in regulating symbiont abundance [38, 39]. In our study, the dominant (i.e., the strain present at > 50% ratio) defined patterns of nitrogen assimilation which affected both AP15Nsym and AP15Nhost. Where S. microadriaticum was the dominant strain, increasing cell densities resulted in decreasing AP15Nhost and AP15Nsym (Fig. 4). This shows that S. microadriaticum are sensitive to nutrient dynamics with the ability to ramp up assimilation where nitrogen is in excess, a functional trait that may allow them to proliferate quickly in newly settled recruits . Given a limited pool of nitrogen within the host habitat , density-dependent uptake may indicate that competition among symbionts (i.e., intraspecific competition among S. microadriaticum) drives nitrogen assimilation rates. However, a closer look at the cumulative amount of newly assimilated nitrogen for an individual holobiont, accumulated across both host and symbiont tissues, shows that total amount of new nitrogen assimilated per recruit drops precipitously with increasing abundances of cells in S. microadriaticum-dominated hosts (Fig. 6). This suggests that less nitrogen overall was available to be assimilated as S. microadriaticum increased in abundance. Critically, this nitrogen limitation was not evident in hosts dominated by B. minutum (Fig. 4B and D) where AP15Nsym and AP15Nhost remained stable across densities.
The cumulative effect of these limitations on nutrient dynamics of the symbiosis was ecologically complex. While AP15N in symbiont tissues was higher or equal in S. microadriaticum-dominated communities, relative to B. minutum-dominated communities, that functional difference did not consistently benefit the host. Instead, at low densities, S. microadriaticum communities were most beneficial to the host in terms of nitrogen (highest AP15Nhost), but at higher symbiont densities, hosts dominated by S. microadriaticum received less nitrogen than those same densities of B. minutum-dominated communities (Fig. 4D). This is not the first evidence of selfishness or parasitism of Symbiodinium spp. (see [30, 43, 48]).
For long-lived hosts, sequential partnering with multiple symbiont species over time, including sub-optimal species, may have advantages. Acacia trees, for example, benefit from associations with a sterilizing ant species which increases survivorship in early ontogeny and is later replaced by ants that increase fecundity . These accumulated effects on lifetime fitness, however, were dependent on the timing and efficiency of turnover between different partners. While poorly understood, a more drastic example of symbiont turnover occurs in B. asbestinum where, by 4 years of age, juveniles have switched between two species of Breviolum ; a similar phenomenon occurs in Acropora juveniles at 3 or more years old . Initial uptake and winnowing may also be an important part of this sequence. Mortality is extremely high for newly settled coral recruits with rates that decrease as coral size increases . Therefore, recruits may initially benefit from highly infectious, and quickly proliferating symbiont species that are ultimately sub-optimal, but only if they can later be replaced by more optimal symbionts. In this case, S. microadriaticum at low densities provided the highest nitrogen benefits to their hosts but became less beneficial relative to B. minutum as symbiont densities increased. Symbiont types within Breviolum ultimately provide the greatest growth and survivorship benefits to B. asbestinum hosts . Ontogenetic flexibility in the regulation of symbiont densities  and identity  has been demonstrated. While a mechanistic understanding of symbiont turnover and the function of diversity in symbiosis is lacking, our stable isotope tracing experiments have provided a snapshot of nutrient transactions between host and symbionts across an array of symbiont community profiles. A host’s ability to regulate symbiont densities in hospite is critical for maintaining a stable symbiosis (Cunning & Baker 2014); here, we suggest that there is a species-specific mechanism for this. Technological advances in stable isotope tracing (e.g., NanoSIMS, compound specific stable isotope analyses) and further combinations with genetic techniques (e.g., qPCR, FISH probing, and flow cytometry sorting [18, 53] can provide further insights into the maintenance and restriction of diversity in symbioses across systems.
This study provides insight into the complex nutrient dynamics that occur between corals and the Symbiodiniaceae that sustain them. While intraspecific competition seemed to be repressed in the host environment, rates of nitrogen assimilation showed dynamic and species-specific regulation. While theory has long suggested that hosts evolve strategies to stabilize mutualisms, these dynamics have been difficult to track. On the basis of our observations, we suggest that the regulation of endosymbiotic communities by coral hosts can lead to more productive mutualisms. As genetic techniques open our eyes to the breadth of symbiont diversity within myriad microbiomes, our understanding of the ecological mechanisms that restrict and maintain diversity in symbiosis must advance. Methodological approaches which provide a more holistic picture of potential conflict and synergy will continue to reveal critical aspects of the ecological and evolutionary dynamics of symbiosis.
Availability of data and materials
All data and code are available at https://github.com/shelby26/NutrientDynamics.
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.
Brown B. Coral bleaching: causes and consequences. Coral Reefs. 1997;16:S129–38.
Oppen MJHV, Oliver JK, Putnam HM, Gates RD. Building coral reef resilience through assisted evolution. Proc Natl Acad Sci. 2015;112:2307–13.
Little AF, Oppen MJH v, Willis BL. Flexibility in algal endosymbioses shapes growth in reef corals. Science. 2004;304:1492–4.
McIlroy SE, Gillette P, Cunning R, Klueter A, Capo T, Baker AC, et al. The effects of Symbiodinium (Pyrrhophyta) identity on growth, survivorship, and thermal tolerance of newly settled coral recruits. Lin S, editor. J Phycol. 2016;52:1114–24.
Tivey TR, Parkinson JE, Weis VM. Host and symbiont cell cycle coordination is mediated by symbiotic state, nutrition, and partner identity in a model cnidarian-dinoflagellate symbiosis. Medina M, Ruby EG, editors. mBio. 2020;11:e02626–19.
Rouzé H, Lecellier G, Saulnier D, Lecellier VB. Symbiodinium clades A and D differentially predispose Acropora cytherea to disease and Vibrio spp. colonization. Ecol Evol. 2016;6:560–72.
Sampayo EM, Ridgway T, Bongaerts P, Hoegh-Guldberg O. Bleaching susceptibility and mortality of corals are determined by fine-scale differences in symbiont type. Proc Natl Acad Sci. 2008;105:10444–9.
Loram JE, Trapido-Rosenthal HG, Douglas AE. Functional significance of genetically different symbiotic algae Symbiodinium in a coral reef symbiosis. Mol Ecol. 2007;16:4849–57 Available from: http://onlinelibrary.wiley.com/doi/10.1111/j.1365-294X.2007.03491.x/abstract.
Cantin NE, Oppen MJH v, Willis BL, Mieog JC, Negri AP. Juvenile corals can acquire more carbon from high-performance algal symbionts. Coral Reefs. 2009;28:405–14.
Baker DM, Andras JP, Jordán-Garza AG, Fogel ML. Nitrate competition in a coral symbiosis varies with temperature among Symbiodinium clades. ISME J. 2013;7:1248–51.
Wall CB, Kaluhiokalani M, Popp BN, Donahue MJ, Gates RD. Divergent symbiont communities determine the physiology and nutrition of a reef coral across a light-availability gradient. ISME J. 2020;14:945–58.
Palmer T, Stanton M, Young T. Competition and coexistence: exploring mechanisms that restrict and maintain diversity within mutualist guilds. Am Nat. 2003;162:S63–79 Available from: http://www.jstor.org/stable/3473453.
Hartmann AC, Baird AH, Knowlton N, Huang D. The paradox of environmental symbiont acquisition in obligate mutualisms. Curr Biol. 2017;27:3711–3716.e3.
Silverstein RN, Correa AMS, Baker AC. Specificity is rarely absolute in coral-algal symbiosis: implications for coral response to climate change. Proc R Soc London Ser B-Biol Sci. 2012;279:2609–18.
Frank S. Host-symbiont conflict over the mixing of symbiotic lineages. Proc R Soc London Ser B-Biol Sci. 1996;263:339–44.
Douglas AE. The symbiotic habit; 2009. p. 1–23.
McIlroy SE, Wong JCY, Baker DM. Competitive traits of coral symbionts may alter the structure and function of the microbiome. ISME J. 2020;14:2424–32.
Barott KL, Venn AA, Perez SO, Tambutté S, Tresguerres M. Coral host cells acidify symbiotic algal microenvironment to promote photosynthesis. Proc Natl Acad Sci. 2015;112:607–12.
Kiers ET, Heijden MGA v d. Mutualistic stability in the arbuscular mycorrhizal symbiosis: exploring hypotheses of evolutionary cooperation. Ecology. 2006;87:1627–36.
Kiers ET, Rousseau RA, West SA, Denison RF. Host sanctions and the legume–rhizobium mutualism. Nature. 2003;425:78–81.
Kiers ET, Duhamel M, Beesetty Y, Mensah JA, Franken O, Verbruggen E, et al. Reciprocal rewards stabilize cooperation in the mycorrhizal symbiosis. Science. 2011;333:880–2.
Cui G, Liew YJ, Li Y, Kharbatia N, Zahran NI, Emwas A-H, et al. Host-dependent nitrogen recycling as a mechanism of symbiont control in Aiptasia. PLoS Genet. 2019;15:e1008189.
Ezzat L, Towle E, Irisson J-O, Langdon C, Ferrier-Pagès C. The relationship between heterotrophic feeding and inorganic nutrient availability in the scleractinian coral T. reniformisunder a short-term temperature increase. Limnol Oceanogr. 2015;61:89–102.
Dubinsky Z, Berman-Frank I. Uncoupling primary production from population growth in photosynthesizing organisms in aquatic ecosystems. Aquat Sci. 2001;63:4–17.
Poland DM, Mansfield JM, Hannes AR, Lewis C, Shearer TL, Connelly SJ, et al. Variation in Symbiodinium communities in juvenile Briareum asbestinum (Cnidaria: Octocorallia) over temporal and spatial scales. Mar Ecol Prog Ser. 2013;476:23–37.
Abrego D, Oppen MJHV, Willis BL. Onset of algal endosymbiont specificity varies among closely related species of Acropora corals during early ontogeny. Mol Ecol. 2009;18:3532–43.
Poland DM, Coffroth MA. Host growth and survivorship varies with endosymbiotic algal partner in developing cnidarians. Mar Ecol Prog Ser Oldendorf. 2019;612:87–100.
Poland DM, Coffroth MA. Trans-generational specificity within a cnidarian–algal symbiosis. Coral Reefs. 2017;36:119–29.
McIlroy SE, Cunning R, Baker AC, Coffroth MA. Competition and succession among coral endosymbionts. Ecol Evol. 2019;9:12767–78.
Santos S, Kinzie R, Sakai K, Coffroth M. Molecular characterization of nuclear small subunit (18S)-rDNA pseudogenes in a symbiotic dinoflagellate (Symbiodinium, Dinophyta). J Eukaryot Microbiol. 2003;50:417–21.
LaJeunesse T. Investigating the biodiversity, ecology, and phylogeny of endosymbiotic dinoflagellates in the genus Symbiodinium using the its region: in search of a “species” level marker. J Phycol. 2001;37:866–80.
Coffroth MA, Lewis CF, Santos SR, Weaver JL. Environmental populations of symbiotic dinoflagellates in the genus Symbiodinium can initiate symbioses with reef cnidarians. Curr Biol. 2006;16:R985–7.
Grover R, Maguer JF, Allemand D, Ferrier-Pages C. Nitrate uptake in the scleractinian coral Stylophora pistillata. Limnol Oceanogr. 2003;48:2266–74.
Coffroth MA, Lasker HR, Diamond ME, Bruenn JA, Bermingham E. DNA fingerprints of a gorgonian coral: a method for detecting clonal structure in a vegetative species. Mar Biol. 1992;114:317–25.
Cunning R, Baker AC. Excess algal symbionts increase the susceptibility of reef corals to bleaching. Nat Climate Change. 2013;3:259–62.
Tanaka Y, Suzuki A, Sakai K. The stoichiometry of coral-dinoflagellate symbiosis: carbon and nitrogen cycles are balanced in the recycling and double translocation system. ISME J. 2018;12:860–8.
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 National Acad Sci. 2021;118:e2022653118.
Krueger T, Horwitz N, Bodin J, Giovani M-E, Escrig S, Fine M, et al. Intracellular competition for nitrogen controls dinoflagellate population density in corals. Proc R Soc B-Biol Sci. 2020;287:20200049.
Muscatine L, Ferrier-Pages C, Blackburn A, Gates RD, Baghdasarian G, Allemand D. Cell-specific density of symbiotic dinoflagellates in tropical anthozoans. Coral Reefs. 1998;17:329–37.
Venn AA, Tambutté E, Lotto S, Zoccola D, Allemand D, Tambutté S. Imaging intracellular pH in a reef coral and symbiotic anemone. Proc Natl Acad Sci U S A. 2009;106:16574–9.
Ezzat L, Fine M, Maguer JF, Grover R, Ferrier-Pagès C. Carbon and nitrogen acquisition in shallow and deep holobionts of the scleractinian coral S. pistillata. Frontiers in Marine. Science. 2017;4:2273.
Baker DM, Freeman CJ, Wong JCY, Fogel ML, Knowlton N. Climate change promotes parasitism in a coral symbiosis. ISME J. 2018;12:921–30.
Cooper TF, Berkelmans R, Ulstrup KE, Weeks S, Radford B, Jones AM, et al. Environmental factors controlling the distribution of Symbiodinium harboured by the coral Acropora millepora on the Great Barrier Reef. PLoS One. 2011;6:e25536.
Ziegler M, Roder C, Büchel C, Voolstra CR. Niche acclimatization in Red Sea corals is dependent on flexibility of host-symbiont association. Mar Ecol Prog Ser. 2015;533:149–61.
Muscatine L, Porter J. Reef corals: mutualistic symbioses adapted to nutrient-poor environments. Bioscience. 1977;27:454–60 Available from: http://www.jstor.org/stable/1297526.
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.
Sachs JL, Wilcox TP. A shift to parasitism in the jellyfish symbiont Symbiodinium microadriaticum. Proc R Soc B. 2006;273:425–9.
Palmer TM, Doak DF, Stanton ML, Bronstein JL, Kiers ET, Young TP, et al. Synergy of multiple partners, including freeloaders, increases host fitness in a multispecies mutualism. Proc Natl Acad Sci U S A. 2010;107:17234–9.
Edmunds PJ, Gates RD. Size-dependent differences in the photophysiology of the reef coral Porites astreoides. Biol Bull. 2004;206:61–4.
Lecointe A, Domart-Coulon I, Paris A, Meibom A. Cell proliferation and migration during early development of a symbiotic scleractinian coral. Proc R Soc B-Biol Sci. 2016;283:20160206.
McIlroy SE, Coffroth MA. Coral ontogeny affects early symbiont acquisition in laboratory-reared recruits. Coral Reefs. 2017;36:927–32.
McIlroy SE, Smith GJ, Geller JB. FISH-flow: a quantitative molecular approach for describing mixed clade communities of Symbiodinium. Coral Reefs. 2014;33:157–67.
We would like to thank Lindsey Mattick, DJ Valint, Jesse Crandall, Anke Kleuter, and the staff at the Long Key Marine Laboratories for field and technical support. David Baker provided feedback on early versions of the manuscript.
Funding for study design and collection of data was provided by NSF OCE-09-26822 (MAC) and NSF DDIG 1311562 (SEM). Data analysis and writing were supported by RGC GRF 17108620 and 17117221 (SEM) and NSF OCE-1559105 (CPt).
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Additional file 1: Table S1.
Nutrient and isotopic compositions of artificial seawater media used for each stage of the experiment. Concentrations were based on Tanaka et al. (2006). All compounds were manufactured and purchased from Sigma-Aldrich ® (St. Louis, MO, USA). Fig. S1. Isotope baseline and dark controls. A set of recruits was tested to provide a baseline isotope value and to confirm that symbionts were the primer drivers of inorganic carbon and nitrogen into the symbiosis. Baseline recruits were sampled prior to the start of the experiment, Dark Enriched samples were exposed to isotopically enriched enriched nutriens but kept in the dark to inhibit photosynthesis; Enriched samples show the cumulative data from the final experiment. Fig. S2. The significant interaction effect of symbiont ratio and symbiont density on atom percent 13C values of symbiont tissues was examined by plotting the marginal effect of symbiont ratio (S.m.:B.m.) as moderated by symbiont density (cells/recruit). Green bar indicates the range of symbiont densities in which there is a significant positive effect of symbiont ratio on AP13Csym. Gray area shows the 95% CI. Fig. S3. The significant interaction effect of symbiont ratio and symbiont density on atom percent 15N values of host tissues was examined by plotting the marginal effect of symbiont density as moderated by symbiont ratio (S.m.:B.m.). Red bar indicates the range of symbiont ratios in which there is a significant negative effect of symbiont density on AP15Nhost. Gray area shows 95% CI. Fig. S4. The significant interaction effect of symbiont ratio and symbiont density on atom percent 15N values of symbiont tissues was examined by plotting the marginal effect of symbiont ratio (S.m.:B.m.) as moderated by symbiont density (cells/recruit). Red bar indicates the range of symbiont densities in which there is a significant negative effect of symbiont ratio on AP15Nsym. Gray area shows the 95% CI. Fig. S5. The significant interaction effect of symbiont ratio and symbiont density on atom percent 15N values of symbiont tissues was examined by plotting the marginal effect of symbiont density as moderated by symbiont ratio (S.m.:B.m.). Red bar indicates the range of symbiont ratios in which there is a significant negative effect of symbiont density on AP15Nsym. Gray area shows 95% CI. Fig. S6. The significant interaction effect of symbiont ratio and symbiont density on total assimilated nitrogen (mg/recruit) was examined by plotting the marginal effect of symbiont density as moderated by symbiont ratio (S.m.:B.m.). Red bar indicates the range of symbiont ratios in which there is a significant negative effect of symbiont density on total assimialted nitrogen. Gray area shows 95% CI.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.
About this article
Cite this article
McIlroy, S.E., terHorst, C.P., Teece, M. et al. Nutrient dynamics in coral symbiosis depend on both the relative and absolute abundance of Symbiodiniaceae species. Microbiome 10, 192 (2022). https://doi.org/10.1186/s40168-022-01382-0
- Symbiosis ecology
- Nutrient exchange
- Stable isotopes