Testing the relationship between microbiome composition and flux of carbon and nutrients in Caribbean coral reef sponges
Microbiome volume 7, Article number: 124 (2019)
Sponges are important suspension-feeding members of reef communities, with the collective capacity to overturn the entire water column on shallow Caribbean reefs every day. The sponge-loop hypothesis suggests that sponges take up dissolved organic carbon (DOC) and, via assimilation and shedding of cells, return carbon to the reef ecosystem as particulate organic carbon (POC). Sponges host complex microbial communities within their tissues that may play a role in carbon and nutrient cycling within the sponge holobiont. To investigate this relationship, we paired microbial community characterization (16S rRNA analysis, Illumina Mi-Seq platform) with carbon (DOC, POC) and nutrient (PO4, NOx, NH4) flux data (specific filtration rate) for 10 common Caribbean sponge species at two distant sites (Florida Keys vs. Belize, ~ 1203 km apart).
Distance-based linear modeling revealed weak relationships overall between symbiont structure and carbon and nutrient flux, suggesting that the observed differences in POC, DOC, PO4, and NOx flux among sponges are not caused by variations in the composition of symbiont communities. In contrast, significant correlations between symbiont structure and NH4 flux occurred consistently across the dataset. Further, several individual symbiont taxa (OTUs) exhibited relative abundances that correlated with NH4 flux, including one OTU affiliated with the ammonia-oxidizing genus Cenarchaeum.
Combined, these results indicate that microbiome structure is uncoupled from sponge carbon cycling and does not explain variation in DOC uptake among Caribbean coral reef sponges. Accordingly, differential DOC assimilation by sponge cells or stable microbiome components may ultimately drive carbon flux in the sponge holobiont.
Sponges (phylum Porifera) are sessile invertebrates that have long been integral members of benthic communities, appearing in the late Precambrian fossil record (580 million years ago) , and contribute to important present-day ecological functions in marine ecosystems. For example, sponges contribute to benthic-pelagic coupling and the mediation of primary production and nitrification via prokaryotic symbionts . Sponges have an extraordinary ability to pump large volumes of water [3,4,5,6] and often retain greater than 80% of filtered particulates , thus serving the important function of nutrient transfer from pelagic to benthic communities. In coral reef habitats, sponge population densities may be increasing , with sponges as abundant as reef-building corals on many Caribbean reefs  and expected to play an increasingly important role in coral reef ecology during a period of broad declines in reef habitats [2, 10].
Sponges support rich and intricate microbial communities [11,12,13,14,15] that differ from seawater and sediments [12, 16,17,18,19]. The microbial communities of sponges can comprise up to 38% of the total tissue biomass (108–1010 bacteria per gram sponge, 2–4 orders of magnitude greater than seawater) in high microbial abundance (HMA) sponges [20, 21]. In low microbial abundance (LMA) sponges, microbial symbionts occur at similar concentrations as seawater (106–108 bacteria cells per gram of sponge) , exhibit lower diversity , and are generally dominated by Proteobacteria [23,24,25]. The microbial communities of both HMA and LMA sponges are generally specific to each sponge species (even across great distances [26, 27]), with host species explaining 64% of observed beta diversity variation in microbial communities  and representing the dominant factor structuring sponge microbial communities for both abundant and rare bacterial taxa . Given the high abundance and host-specificity of microbial communities in sponges, it is suspected that microbes may contribute to carbon and nutrient flux and processing in host sponges.
Recent work suggests that sponges recycle carbon within coral reef ecosystems, consuming dissolved organic matter and excreting particulate matter, a process potentially involving microbial symbionts. This hypothesis, termed the sponge-loop hypothesis , is similar to the microbial loop theory proposed by Azam and colleagues , in which free-living marine microbes use dissolved organic carbon (DOC), one of the largest organic carbon reservoirs on earth , and convert DOC into particulate organic carbon (POC) for consumption by higher trophic levels in pelagic food webs. Similarly, the sponge-loop hypothesis proposes the cycling of DOC into POC, advancing a two-step process by which sponges mediate this conversion: (1) DOC uptake and (2) detritus excretion . Using isotopic tracers, a recent study demonstrated that DOC from coral mucus can be directly transferred into the bulk tissue of the warm-water sponge Mycale fistulifera and the cold-water sponge Hymedesmia coriacea . Other studies have suggested that DOC exuded by corals , algae , or seagrasses  can be consumed by sponges, and these data provide support for the first component of the sponge-loop hypothesis in different hosts (HMA and LMA) and environments (shallow-water and deep-water reefs). The second component, detritus excretion via rapid cell turnover (cell proliferation and shedding), may not be as widespread a phenomenon as DOC uptake, occurring in cryptic, encrusting sponge species [29, 32, 33], but absent in massive, emergent sponge species .
While the role of microbial symbionts in the sponge-loop is unknown, the abundance and composition of microbial symbionts have been shown to affect some aspects of sponge physiology and nutrient cycles. Such patterns most notably occur across the HMA-LMA spectrum, where microbiome distinctions correlate with differences in feeding behaviors. For example, LMA and HMA sponges exhibit different pumping rates, with LMA sponges pumping faster and exhibiting greater choanocyte density than HMA sponges [35,36,37]. Further, differences in pumping rates, aquiferous system density, and microbial communities between HMA and LMA sponges may also affect host carbon and nutrient flux. In general, HMA sponges exhibit greater uptake rates of DOC than LMA sponges [35, 38, 39] and different inorganic nitrogen fluxes . Indeed, recent work suggests that differential nitrogen cycling between sympatric LMA and HMA sponges results in trophic niche separation, thereby facilitating co-existence and efficient nutrient utilization in oligotrophic environments . The relationships between microbial symbiont communities, DOC flux, and nutrient cycling in coral reef sponges may be important for understanding recent changes to coral reef ecosystems, particularly in the Caribbean .
Comparisons of carbon and nutrient flux across the HMA-LMA spectrum have yielded important insights into sponge-mediated nutrient cycling, yet most studies lack comprehensive microbiome characterizations within these broad host categories. Sponge microbiomes differ markedly within the HMA and LMA categories, as most exhibit a species-specific signature . Accordingly, a direct test of the relationship between microbiome composition and sponge nutrient flux has not been conducted. This study investigated whether sponges with different rates of nutrient flux hosted different microbial symbionts, at the community and individual taxa levels, by assessing correlations between DOC, POC, and nutrient flux and the structure and composition of microbial communities in sponges. We characterized the microbiomes of 10 common, emergent (i.e., non-cryptic and non-encrusting) Caribbean coral reef sponge species from two locations to investigate the drivers of inter- and intra-specific variation in sponge microbiomes and determine the relationship between symbiont community variation and carbon and nutrient fluxes.
Ambient seawater and sponge tissue were collected from two geographically distant locations (~ 1203 km apart): Conch Reef, Florida (24° 56.9′ N, 80° 27.2′ W), and Carrie Bow Cay, Belize (16° 48.14′ N, 88° 4.79′ W). Ten of the most common Caribbean coral reef sponges were sampled in Key Largo, Florida (Conch Reef), and eight of these species were sampled in Belize (Carrie Bow Cay) from 13 to 23 m depths in June and July 2016, respectively (Table 1). Only apparently healthy sponge individuals (i.e., no evidence of disease, tissue damage, algal colonization, or epibionts) with a single osculum were sampled (except in the case of Agelas tubulata which had multiple oscula) . All sponge tissue samples were collected in separate bags, preserved in 100% ethanol, and stored at − 20 °C until processing. Seawater samples (1 L) were collected at each sampling site and day of sponge tissue sample collection, concentrated onto 0.2 μm filters, preserved in 100% ethanol, and stored at − 20 °C until processing.
All ten sponge species represent common Caribbean coral reef species and were identified morphologically following Zea et al. . To confirm identifications made in the field, PCR amplification of the partial mitochondrial cytochrome c oxidase subunit I (COI) gene was amplified using the forward primer LCO1490 and reverse primer HCO2198 for species barcoding . PCR amplification reactions contained 0.5 μl of each primer (10 μM), 12.5 μl (0.5 units) of MyTaq™ Red Mix DNA polymerase (Bioline), 1 μl of DNA template, and PCR water for a total reaction volume of 25 μl. The thermocycler conditions included initial denaturation step (95 °C, 1 min) followed by 35 cycles of denaturation (95 °C, 15 s), annealing (45 °C, 15 s), and extension (72 °C, 10 s), with a final extension step (72 °C, 1 min) and 6 °C hold. The COI amplicons were used in a sequencing PCR with BigDye version 3.1 (Applied Biosystems) and a thermocycler program consisting of an initial denaturation step (96 °C, 1 min), 25 cycles of annealing (50 °C, 5 s), extension (60 °C, 4 min), and denaturation (96 °C, 10 s), followed by a final annealing (50 °C, 5 s), extension (60 °C, 4 min), and 10 °C hold. Amplicons were cleaned using BigDye® XTerminator™ Purification Kit (Thermo Fischer), following the manufacturer’s protocol, and sequenced on an AB 3500 Gene Analyzer (Applied Biosystems) at the UNCW Center for Marine Science. Forward and reverse sequences were aligned in Geneious version 8.1.9  to create consensus sequences and compared to the GenBank database using the nucleotide-nucleotide Basic Local Alignment Search Tool (BLASTn). Sequence data were deposited in GenBank under the accession numbers MH297440 to MH297461.
For each sponge tissue sample that was collected, metadata on sponge pumping rates, sponge volumes, and carbon fluxes were collected on in situ colonies (prior to tissue sampling) and processed as reported previously . Briefly, sponge pumping rates were measured using an acoustic Doppler velocimeter (SonTek)  and sponge tissue volume was estimated using measurements of the dimensions of each sponge . Paired incurrent (ambient) and excurrent seawater samples (1.5 L) were collected via syringe and subsequently filtered (Whatman GF/F). POC on filters was quantified using a CE Elantech NC2100 elemental analyzer [35, 46], and DOC in the filtrate of each sample was quantified using a Shimadzu TOC 5050 analyzer . In Belize, an additional 40 mL of the filtrate from each incurrent and excurrent seawater sample was collected and stored frozen until quantification of NOx, NH4, and PO4 using a Bran + Luebbe AutoAnalyzer III following standard protocols . The specific filtration rate (SFR, μmol C, or nutrient/s/L sponge) for each carbon and nutrient species (i.e., POC, DOC, NOx, NH4, and PO4) was calculated as:
where Cin and Cex are the incurrent and excurrent concentrations of each carbon pool or nutrient type (C/mL), Vsponge is the sponge tissue volume (L), and Q is the pumping rate for each sponge (mL/s); thus, positive values indicate net consumption and negative values indicate net production of a particular carbon pool or nutrient type. Carbon data is available in McMurray et al.  and nutrient data in Additional file 1.
DNA extraction and sequence processing
Ethanol-preserved tissue samples were dissected into 2 mm3 cubes that included interior and exterior sponge tissue and were extracted using the DNeasy® Blood & Tissue Kit (Qiagen) following the manufacturer protocols. Partial (V4) 16S rRNA gene sequences were amplified using the 515f forward primer and 806r reverse primer  and sequenced on an Illumina MiSeq platform at Molecular Research LP (Shallowater, TX). Illumina sequence reads were processed in mothur v1.38.0  using a modified version of the bioinformatics pipeline described in Weigel and Erwin . Briefly, raw sequences (n = 13.9 million) were demultiplexed, quality-filtered, aligned, classified, and clustered into operational taxonomic units (OTUs) at 97% sequence identity (nOTU = 25,712). Sequence libraries for each sample were subsampled to the lowest read count (n = 15,825), and all data analyses were based on the subsampled dataset. Sequence data were deposited in the Sequence Read Archive of the National Center for Biotechnology under accession number SRP142647.
Microbial community diversity
Diversity statistics (Shannon-Weaver, OTU richness, Simpson) were calculated in mothur using OTU relative abundance data. Two-way nested analyses of variance (ANOVA) were used to test for significant differences in diversity indices for two factors: “source” (sponge species or seawater) and “location” (Key Largo or Belize) nested within source in JMP (version 12.0), followed by Tukey’s honest significant difference (HSD) tests to assess multiple post hoc comparisons of means.
Microbial community structure
Bray-Curtis similarity matrices were constructed using square root transformed OTU relative abundance data to give more even representation of rare and abundant taxa in community comparisons using Primer-e (version 6.1.11) and visualized in a cluster dendrogram and a two-dimensional non-metric multidimensional scaling (nMDS) plot. A permutational multivariate analysis of variance (PERMANOVA, version 1.0.1) was conducted to test for significant differences in microbial community structure across two factors: “source” and “location” nested within source, with significance determined by Monte Carlo asymptotic P values.
Similarity percentage (SIMPER) analysis was used to identify individual OTUs driving the overall dissimilarity between microbial communities within each sponge species, using OTU relative abundance matrices and a cumulative dissimilarity cutoff percentage of 0.70. Significant differences in OTU relative abundances were assessed using Metastats  in mothur with 1000 permutations. OTUs of interest were then taxonomically identified and compared to OTUs identified in nutrient correlations (see below).
Correlations between microbial community diversity and carbon/nutrient flux
Spearman rank-order correlations were conducted in SigmaPlot (version 11) to assess relationships between microbial diversity indices (Shannon-Weaver, OTU richness, Simpson) and carbon/nutrient flux data (DOC, POC, NOx, NH4, and PO4).
Correlations between microbial community structure and carbon/nutrient flux
Distance-based linear models (DistLM) were conducted in Primer-e to assess correlations between microbial community (Bray-Curtis) similarity and carbon/nutrient flux data (DOC, POC, NOx, NH4, and PO4) and visualized with distance-based redundancy analysis (dbRDA) plots. Specifically, carbon and nutrient flux measurements estimated from specific filtration rates (SFR, μmol C, or nutrient/s/L sponge) were tested as marginal predictor variables for microbial community similarity within and among sponge species. Analyses were repeated using Bray-Curtis similarity matrices constructed from untransformed (raw) data, showing identical significance patterns and thus minimal impact of data transformation on statistical results (Additional files 2 and 3). Analyses were also repeated using carbon and nutrient uptake (Cin–Cex) instead of specific filtration rates, showing similar statistical patterns when considering a metric not influenced by pumping rates (Additional files 4 and 5).
Inter-specific comparisons were conducted at three levels: all sponge species, within each sponge category (HMA, LMA), and between each pairwise species comparison. Pearson correlations were run in JMP to compare nutrient flux data (POC and NH4) in the form of SFR and OTU abundance counts for the first 1000 OTUs within each sponge species.
Functional gene PCR screening
To test for the presence of signature functional genes involved in nitrogen cycling processes (nitrification: ammonia monooxygenase, amoA/amoB; nitrogen fixation: nitrogenase, nifH; and denitrification: nitrite reductase, nirS), the following primer pairs and cited procedures were used: Arch-amoAF and Arch-amoAR , AmoBMF and AmoBMR , nif1/2 and nif3/4 , nirS1F and nirS6R . For all sponge species, at least two replicates were tested for the presence of each functional gene, except for Aplysina archeri (only one sample available). If positive, all remaining replicates of that species were tested at both sites. If negative, the PCR was repeated twice for verification.
Host species effects (inter-specific variation)
Molecular barcoding confirmed the species identifications from the field (Additional file 6). Microbial communities differed significantly in richness and diversity among sponge species and seawater sources (Shannon-Weaver, P < 0.0001; OTU Richness, P < 0.0001; Simpson, P < 0.0001, Table 1). Similarly, significant differences in microbial community structure were detected among sponge species and seawater sources (P = 0.001). A dendrogram based on Bray-Curtis similarity of microbial communities revealed two main branches, one consisting of all HMA sponge species and a second branch consisting of all LMA sponge species plus seawater (Fig. 1). Within each main branch, samples clustered by source (sponge species), except the closely related species Verongula gigantea and V. reiswigi (Fig. 1). Similar clustering patterns were revealed by nMDS ordination and highlighted a distinct cluster for Mycale laxissima within the LMA sponge species (Additional file 7). Differences in microbial community composition were also observed between HMA and LMA sponges, with LMA microbiomes composed primarily of Alpha- and Gammaproteobacteria and dominated by a small number of OTUs within these taxonomic lineages (Fig. 2). For example, a single OTU (001), an Alphaproteobacterium in the genus Roseivivax, comprised 60% of all Niphates digitalis sequence reads and 77% of sequence reads in M. laxissima was affiliated with the orders Rhodospirillales (50.4%, class Alphaproteobacteria) and Oceanospirillales (26.4%, class Gammaproteobacteria). In contrast, HMA sponges exhibited a distribution of phyla that was more even, including Actinobacteria, Acidobacteria, Chloroflexi, Cyanobacteria, Poribacteria, and Proteobacteria (Fig. 2).
Biogeographic effects (intra-specific variation)
Microbial communities differed significantly in richness and diversity among the same sponge species collected from different locations (Shannon-Weaver, P = 0.004; richness, P < 0.001; Simpson P = 0.010). Pairwise tests revealed that these trends were driven by a single species for each metric: only the microbial communities in M. laxissima showed significant differences in diversity across locations (higher in Florida; Shannon, P = 0.041; Simpson, P = 0.005), while OTU richness was only significantly different within microbial communities of Callyspongia plicifera between sites (higher in Florida; P = 0.007). Accordingly, most species exhibited similar levels of alpha-diversity between locations (Table 1). In contrast, the community structure of sponge microbiomes differed significantly within host species across locations (P = 0.001), with three LMA sponge species, one HMA species, and seawater all exhibiting significant differences in microbial community structure between locations: C. plicifera (P = 0.006), N. digitalis (P = 0.019), M. laxissima (P = 0.027), Xestospongia muta (P = 0.012), and seawater (P = 0.024). Further, the Bray-Curtis similarity-based dendrogram revealed that all samples clustered by location (Florida or Belize) within each sponge species or seawater source (Fig. 1).
Individual symbiont OTUs driving these community-level biogeographic trends were determined using similarity percentage (SIMPER) and differential relative abundance (Metastats) analyses. A total of 58 OTUs exhibited significantly different abundances between sampling locations across all sponge species, with nearly all (98%) representing rare members of the microbiome with low individual contributions to community dissimilarity (< 3%, Additional file 8). However, a single OTU (003, Alphaproteobacteria) that dominated symbiont communities in M. laxissima exhibited significantly higher relative abundance in sponges from Belize (55.2%) than Florida (39.6%) and accounted for > 17% of community dissimilarity between locations. Of these 58 OTUs, 6 were present in 2 or more sponge species, with 1 OTU (006, Synechococcus) significantly differing between locations for all LMA sponges (Additional file 8). Additionally, OTU 006 was also detected at significantly higher relative abundances in seawater samples from Florida compared to Belize, comprising 7.2% of sequence reads within Florida seawater samples. Twenty-nine OTUs were found to significantly differ between locations in HMA sponges (7 in A. tubulata, 4 in V. gigantea, 4 in V. reiswigi, and 15 in X. muta), and 31 OTUs were found to significantly differ between locations in LMA sponges (15 in C. plicifera, 4 in Callyspongia vaginalis, 10 in M. laxissima, and 8 in N. digitalis; Additional file 8). Notably, these OTUs contributed to twice the community dissimilarity in LMA sponges (15%) compared to HMA hosts (7.2%). Nearly one third of these OTUs (n = 18, 31%) were affiliated with the phylum Proteobacteria, with the majority (n = 10, 55%) of proteobacterial OTUs belonging specifically to the class Gammaproteobacterium.
Carbon, nutrient, and microbial abundance correlations
The comparative DistLM analyses showed that microbiome structure varied significantly with POC flux among, but not within, sponge species (Table 2). Significant correlations were detected at inter-specific levels when considering the microbial communities of all host species (P = 0.001) and only LMA hosts (P = 0.002), while variation in POC flux among HMA hosts did not correlate with variability in microbiome structure (P = 0.403, Table 2). Accordingly, pairwise species comparisons were significant for half of the LMA pairs (notably, all involving M. laxissima) and most of the LMA-HMA pairs (except those involving C. vaginalis and Verongula spp.), but none of the HMA pairs (Table 2). In contrast, no significant DistLM correlations were observed between microbiome structure and host DOC flux at any level (inter- and intraspecific, Table 2). No significant correlations were detected with host PO4 and NOx flux, but host NH4 flux correlated significantly with microbiome structure among sponge species (P = 0.002, Table 3). Inter-specific comparisons within LMA and HMA categories were not significant (P = 0.073 and 0.549, respectively), with half of the individually paired LMA-HMA sponge species comparisons significant (P > 0.05, Table 3). Comparisons of carbon/nutrient flux and microbiome diversity revealed similar patterns: no significant correlations with host DOC, PO4, and NOx flux for any alpha-diversity metric, but POC correlated significantly with observed richness (P < 0.001) and NH4 with Shannon-Weaver and Simpson indices (P < 0.001, Additional file 9). Differential carbon and nutrient consumption between LMA and HMA sponges was visualized in dbRDA plots, highlighting the strong correlations between microbiome structure and POC and NH4 flux, and the weak correlations with DOC, PO4, and NOx (Fig. 3).
In addition to differences at the community level, individual symbiont OTUs that exhibited significant correlations with NH4 (n = 132, Additional file 10) and POC (n = 288, Additional file 11) flux were identified within each host species. Most of these OTUs (> 94%) were rare microbiome members (< 1% relative abundance) with only 11 OTUs also exhibiting significant variation within host species (and where available) across locations (Metastats, P < 0.05; Table 4). Of these 11 OTUs, 1 OTU each belonged to A. tubulata, V. gigantea, and V. reiswigi; 3 belonged to N. digitalis; and 5 belonged to X. muta (Table 4). Only 8 of these 11 OTUs comprised ≥ 1% of average relative abundance in the microbial community of their respective sponge: OTU 060 (Piscirickettsiaceae, V. gigantea), 070 (Acidomicrobiales, V. reiswigi), 035 (Nitrospiraceae, X. muta), 067 (AncK6, X. muta), 081 (Chloroflexi, X. muta), 090 (Cenarchaeum, A. tubulata), and 011 (Alphaproteobacteria, N. digitalis; Table 4).
Functional gene PCR screening
Archaeal ammonia monooxygenase (amoA) genes were successfully amplified in nearly all HMA sponge species (A. tubulata, V. gigantea, V. reiswigi, and X. muta), while nitrite reductase (nirS) genes were only amplified in A. tubulata (Florida and Belize) and V. reiswigi (Belize, Table 5). Nitrogenase (nifH) and bacterial ammonia monooxygenase (amoB) genes were not successfully amplified in any samples, despite triplicate reactions attempted.
Matching emerging trends in the field of sponge microbiology, the microbial communities in 10 common Caribbean sponge species were significantly different from seawater [16, 17] and exhibited a high degree of host specificity [16, 26], with greater intra- than interspecific similarity across locations [28, 37, 55]. Significant differences in diversity and composition between microbial communities of HMA and LMA sponges were also observed, with LMA sponge microbial communities exhibiting lower diversity and higher relative abundances of Proteobacteria, consistent with previous work . One LMA sponge species, M. laxissima, clustered separately from other LMA hosts and was recently reported to exhibit slower pumping rates and differential nutrient utilization compared to other LMA sponge species from the same Caribbean sites . In addition, a significant effect of location on microbial diversity and composition was detected within each host sponge species, indicating small intra-specific shifts in microbial communities across sites. These biogeographic differences were primarily manifested in beta-diversity metrics (i.e., compositional differences), although significant differences in alpha-diversity were observed for two LMA hosts (richness in C. plicifera and Shannon-Weaver index in M. laxissima). Accordingly, our microbiome dataset encompassed variability among and within sponge species and allowed for correlations at multiple levels between variability in microbiome structure and variation in carbon and nutrient fluxes.
Surprisingly, no significant correlations between symbiont structure and DOC flux were detected among or within any sponge species. POC flux did track with symbiont structure across the HMA-LMA spectrum but explained little of the variation in POC flux among sponge species within these categories and none of the variation within sponge species. Matching these findings, individual OTUs found to correlate with POC flux were mostly rare (< 1% relative abundance), with only four OTUs comprising > 1% relative abundance. Together, these results indicate that microbial communities do not exhibit clear and consistent structural shifts with variable rates of holobiont carbon flux and thus do not explain variability in DOC uptake among sponges as hypothesized previously [35, 39]. These findings may indicate that sponge cells, not microbial cells, are the primary site of DOC uptake, as suggested by recent experiments tracing DOC incorporation into bacteria-specific and sponge-specific phospholipid-derived fatty acids . Alternatively, DOC uptake may only occur in a specific portion of the sponge microbiome (e.g., stable or functionally redundant taxa), thereby dissociating overall symbiont structure and carbon cycling. Indeed, recent genomic evidence from sponge symbiont phyla enriched in HMA sponges revealed a complex suite of carbohydrate-degrading genes [56, 57]. In either case, our results suggest that divergent physiologies and symbiont abundances across the HMA-LMA spectrum play a greater role in DOC uptake than the fine-scale (OTU level) composition of the microbial consortia in any given species. Compositional insights from amplicon sequence data represent an important first step in characterizing microbiomes and assessing the impacts of environmental factors on microbiome structure; however, there are also technical limitations associated with the nature of relative abundance data . Thus, correlations between absolute symbiont abundance and DOC uptake may not be detected by the methods employed herein and remain a target for future study.
In contrast to the observed decoupling of symbiont structure and carbon flux, microbial community structure correlated strongly with NH4 flux across the HMA-LMA spectrum. Previous work has shown greater NH4 consumption by HMA than LMA sponges [34, 40, 41] and different nitrogen trophic levels for Caribbean HMA versus LMA sponges , indicating differential nitrogen cycling pathways in HMA and LMA sponge microbiomes. Our results show that variations in NH4 flux between HMA and LMA sponge hosts specifically track with microbiome composition and suggest that nitrification is a key energy generation process in the sponge microbiome with cascading effects on the entire prokaryotic symbiont community. Supporting these conclusions, symbiont taxa affiliated with ammonia-oxidizing (Thaumarchaeota ) and nitrite-oxidizing (Nitrospirae) lineages were more common (an order of magnitude greater relative abundance) in HMA compared to LMA microbiomes. Further, we identified individual symbiont OTUs that correlated with NH4 flux, including the archaeal OTU-090 (Thaumarchaeota) in the HMA host A. tubulata whose relative abundance increased with greater levels of NH4 uptake. In general, prokaryotic functional guilds involved in nitrogen cycling are restricted to specific taxonomic groups, whereas DOC uptake is a general process common to nearly all prokaryotic taxa. Thus, the strong link between symbiont structure and nitrogen cycling observed among the holobionts of the sponges investigated herein, and the weak link with carbon cycling, may result from the narrow phylogenetic distribution of functional guilds in the nitrogen cycle.
The presence and expression of functional genes encoding for key metabolic enzymes offers additional insight into nutrient cycling within the sponge microbiome. Herein, we detected ammonia monooxygenase (amoA) genes in most HMA sponges (four of six species), supporting the presence of nitrifying symbionts that use ammonia as a substrate for energy generation. Functional gene screening also revealed the presence of nitrite reductase (nirS) genes in two HMA sponge species, indicating the potential for denitrification pathways in some sponge microbiomes. Indeed, sponge microbial communities have been shown to contain members capable of numerous nitrogen transformation pathways, including nitrogen fixation, denitrification, and nitrification . Notably, our study did not detect nitrogen-fixing bacterial taxa or nifH genes in Caribbean sponges, despite their previous detection in the same sponge species (Ircinia strobilina and M. laxissima) from the same site (Key Largo, Florida) , and primers targeting amoB gene subunits were negative, despite amplification of amoA gene subunits. Technical reasons may account for the lack of gene detection herein, and further metagenomic studies are required to confirm these results. Similar PCR screening for genes involved in carbon flux (e.g., DOC transporters) is complicated by the vast diversity of these genes, but is now approachable using metagenomic (or metatranscriptomic) techniques (e.g., Poretsky et al. ). Such approaches to symbiont characterization will aid in clarifying the presence and activity of symbiont functional guilds in the sponge microbiome, as well as the relationship between symbiont activity and carbon and nutrient flux.
Our results reveal novel insights into the relationship between symbiont structure and holobiont carbon and nutrient flux, while also confirming previous findings regarding the drivers of sponge microbial community structure and differences across the HMA-LMA spectrum. Our results show sponge microbial communities are structured by host species and, to a lesser degree, biogeography, yet inter- and intra-specific variation in sponge microbiomes are uncoupled from sponge carbon flux. As such, these results do not support previous theories that variations in microbial community structure are specifically related to differential DOC flux in host sponges [39, 62, 63]. Rather, the relationship between ambient DOC, seawater flux, and the thresholds in DOC flux previously observed [38, 41, 46] is likely due to differential assimilation by sponge cells [33, 64] or is mediated by stable members of sponge microbiomes. Future research assessing the differential regulation of genes involved in DOC uptake (e.g., DOC transporters) may provide additional insight into the physiology of stable microbiome members and cascading effects on holobiont carbon flux. In contrast to carbon flux, the correlational and functional gene data support previous findings of sponge microbial communities having the ability to participate in nitrogen cycling [65,66,67] and implicate nitrification as a key metabolic process affecting overall microbiome structure in Caribbean sponges . Together, these results indicate that the sponge microbiome may play important roles in nutrient cycling within coral reef ecosystems and should be considered when assessing the ecological impacts of sponges on reef communities, including nutrient feedback loops causing ecosystem phase shifts on coral reefs.
Availability of data and materials
The 16S rRNA sequence datasets generated and analyzed in the current paper are available in the Sequence Read Archives (SRP142647) May 31, 2019, at the following link https://www.ncbi.nlm.nih.gov/sra/SRP142647. The host sequences generated in the current study have been submitted to GenBank under the title “Microbial Symbionts, Carbon and Nutrient Cycling in Caribbean Coral Reef Sponges”.
Li C-W, Chen J-Y, Hua T-E. Precambrian sponges with cellular structures. Science. 1998;279(5352):879–82.
Diaz MC, Rutzler K. Sponges: an essential component of Caribbean coral reefs. Bull Mar Sci. 2001;69(2):535–46.
Reiswig HM. Particle feeding in natural populations of three marine Demosponges. Biol Bull. 1971;141(3):568–91.
Turon X, Galera J, Uriz MJ. Clearance rates and aquiferous systems in two sponges with contrasting life-history strategies. J Exp Zool. 1997;278(1):22–36.
Ribes M, Coma R, Gili J-M. Natural diet and grazing rate of the temperate sponge Dysidea avara (Demospongiae, Dendroceratida) throughout an annual cycle. Mar Ecol Prog Ser. 1999;176:179–90.
McMurray SE, Pawlik JR, Finelli CM. Trait-mediated ecosystem impacts: how morphology and size affect pumping rates of the Caribbean giant barrel sponge. Aquat Biol. 2014;23(1):1–13.
Maldonado M, Ribes M, van Duyl FC. Nutrient fluxes through sponges: biology, budgets, and ecological implications. In: Adv Mar Biol Elsevier; 2012. p. 113–82.
McMurray SE, Henkel TP, Pawlik JR. Demographics of increasing populations of the giant barrel sponge Xestospongia muta in the Florida Keys. Ecology. 2010;91(2):560–70.
Loh T-L, McMurray SE, Henkel TP, Vicente J, Pawlik JR. Indirect effects of overfishing on Caribbean reefs: sponges overgrow reef-building corals. Peer J. 2015;3:e901.
Bell JJ. The functional roles of marine sponges. Estuar Coast Shelf Sci. 2008;79(3):341–53.
Hentschel U, Piel J, Degnan SM, Taylor MW. Genomic insights into the marine sponge microbiome. Nat Rev Microbiol. 2012;10(9):641–54.
Taylor MW, Radax R, Steger D, Wagner M. Sponge-associated microorganisms: evolution, ecology, and biotechnological potential. Microbiol Mol Biol Rev. 2007;71(2):295–347.
Taylor MW, Tsai P, Simister RL, Deines P, Botte E, Ericson G, Schmitt S, Webster NS. ‘Sponge-specific’ bacteria are widespread (but rare) in diverse marine environments. ISME J. 2013;7(2):438–43.
Thomas T, Moitinho-Silva L, Lurgi M, Bjork JR, Easson C, Astudillo-Garcia C, Olson JB, Erwin PM, López-Legentil S, Luter H, et al. Diversity, structure and convergent evolution of the global sponge microbiome. Nat Commun. 2016;7:11870.
Thomas T, Rusch D, DeMaere MZ, Yung PY, Lewis M, Halpern A, Heidelberg KB, Egan S, Steinberg PD, Kjelleberg S. Functional genomic signatures of sponge bacteria reveal unique and shared features of symbiosis. ISME J. 2010;4(12):1557–67.
Hentschel U, Hopke J, Horn M, Friedrich AB, Wagner M, Hacker J, Moore BS. Molecular evidence for a uniform microbial community in sponges from different oceans. Appl Environ Microbiol. 2002;68(9):4431–40.
Weigel BL, Erwin PM. Intraspecific variation in microbial symbiont communities of the sun sponge, Hymeniacidon heliophila, from intertidal and subtidal habitats. Appl Environ Microbiol. 2016;82(2):650–8.
Montalvo NF, Hill RT. Sponge-associated bacteria are strictly maintained in two closely related but geographically distant sponge hosts. Appl Environ Microbiol. 2011;77(20):7207–16.
Fieseler L, Horn M, Wagner M, Hentschel U. Discovery of the novel candidate phylum “Poribacteria” in marine sponges. Appl Environ Microbiol. 2004;70(6):3724–32.
Vacelet J. Etude en microscopie electronique de l’association entre bacteries et spongiaires du genre Verongia (Dictyoceratida). J Microsc Biol Cell. 1975;23(3):271–88.
Hentschel U, Usher KM, Taylor MW. Marine sponges as microbial fermenters. FEMS Microbiol Ecol. 2006;55(2):167–77.
Gloeckner V, Wehrl M, Moitinho-Silva L, Gernert C, Schupp P, Pawlik JR, Lindquist NL, Erpenbeck D, Wörheide G, Hentschel U. The HMA-LMA dichotomy revisited: an electron microscopical survey of 56 sponge species. Biol Bull. 2014;227(1):78–88.
Webster SN, Hill TR. The culturable microbial community of the Great Barrier Reef sponge Rhopaloeides odorabile is dominated by an α-Proteobacterium. Mar Biol. 2001;138(4):843–51.
Croué J, West NJ, Escande M-L, Intertaglia L, Lebaron P, Suzuki MT. A single betaproteobacterium dominates the microbial community of the crambescidine-containing sponge Crambe crambe. Sci Rep. 2013;3:2583.
Giles EC, Kamke J, Moitinho-Silva L, Taylor MW, Hentschel U, Ravasi T, Schmitt S. Bacterial community profiles in low microbial abundance sponges. FEMS Microbiol Ecol. 2013;83(1):232–41.
Lee OO, Wang Y, Yang J, Lafi FF, Al-Suwailem A, Qian P-Y. Pyrosequencing reveals highly diverse and species-specific microbial communities in sponges from the Red Sea. ISME J. 2011;5(4):650–64.
Pita L, Turon X, López-Legentil S, Erwin PM. Host rules: spatial stability of bacterial communities associated with marine sponges (Ircinia spp.) in the Western Mediterranean Sea. FEMS Microbiol Ecol. 2013;86(2):268–76.
Reveillaud J, Maignien L, Eren MA, Huber JA, Apprill A, Sogin ML, Vanreusel A. Host-specificity among abundant and rare taxa in the sponge microbiome. ISME J. 2014;8(6):1198–209.
de Goeij JM, van Oevelen D, Vermeij MJA, Osinga R, Middelburg JJ, de Goeij AFPM, Admiraal W. Surviving in a marine desert: the sponge loop retains resources within coral reefs. Science. 2013;342(6154):108–10.
Azam F, Fenchel T, Field J, Gray J, Meyer-Reil L, Thingstad F. The ecological role of water-column microbes in the sea. Mar Ecol Prog Ser. 1983;10(3):257–63.
Hansell DA, Carlson CA, Repeta DJ, Schlitzer R. Dissolved organic matter in the ocean: a controversy stimulates new insights. Oceanography. 2009;22(4):202–11.
Rix L, de Goeij JM, Mueller CE, Struck U, Middelburg JJ, van Duyl FC, Al-Horani FA, Wild C, Naumann MS, van Oevelen D. Coral mucus fuels the sponge loop in warm- and cold-water coral reef ecosystems. Sci Rep. 2016;6:18715.
Rix L, de Goeij JM, van Oevelen D, Struck U, Al-Horani FA, Wild C, Naumann MS. Differential recycling of coral and algal dissolved organic matter via the sponge loop. Funct Ecol. 2017;31(3):778–89.
Weisz JB, Hentschel U, Lindquist N, Martens CS. Linking abundance and diversity of sponge-associated microbial communities to metabolic differences in host sponges. Mar Biol. 2007;152(2):475–83.
McMurray SE, Stubler AD, Erwin PM, Finelli CM, Pawlik JR. A test of the sponge-loop hypothesis for emergent Caribbean reef sponges. Mar Ecol Prog Ser. 2018;588:1–14.
Weisz JB, Lindquist N, Martens CS. Do associated microbial abundances impact marine demosponge pumping rates and tissue densities? Oecologia. 2008;155(2):367–76.
Poppell E, Weisz J, Spicer L, Massaro A, Hill A, Hill M. Sponge heterotrophic capacity and bacterial community structure in high- and low-microbial abundance sponges. Mar Ecol. 2014;35(4):414–24.
McMurray SE, Pawlik JR, Finelli CM. Demography alters carbon flux for a dominant benthic suspension feeder, the giant barrel sponge, on Conch Reef, Florida Keys. Funct Ecol. 2017;31(11):2188–98.
Hoer DR, Gibson PJ, Tommerdahl JP, Lindquist NL, Martens CS. Consumption of dissolved organic carbon by Caribbean reef sponges. Limnol Oceanogr. 2017;63(1):337–51.
Southwell MW, Weisz JB, Martens CS, Lindquist N. In situ fluxes of dissolved inorganic nitrogen from the sponge community on Conch Reef, Key Largo, Florida. Limnol Oceanogr. 2008;53(3):986–96.
Morganti T, Coma R, Yahel G, Ribes M. Trophic niche separation that facilitates co-existence of high and low microbial abundance sponges is revealed by in situ study of carbon and nitrogen fluxes. Limnol Oceanogr. 2017;62(5):1963–83.
Pawlik JR, Burkepile DE, Thurber RV. A vicious circle? Altered carbon and nutrient cycling may explain the low resilience of Caribbean coral reefs. BioScience. 2016;66(6):470–6.
Zea S, Henkel TP, Pawlik JR. The Sponge Guide: a picture guide to Caribbean sponges, vol. vol.3; 2014.
Folmer O, Black M, Hoeh W, Lutz R, Vrijenhoek R. DNA primers for amplification of mitochondrial cytochrome c oxidase subunit I from diverse metazoan invertebrates. Mol Mar Biol Biotechnol. 1994;3(5):294.
Kearse M, Moir R, Wilson A, Stones-Havas S, Cheung M, Sturrock S, Buxton S, Cooper A, Markowitz S, Duran C. Geneious Basic: an integrated and extendable desktop software platform for the organization and analysis of sequence data. Bioinformatics. 2012;28(12):1647–9.
McMurray SE, Johnson ZI, Hunt DE, Pawlik JR, Finelli CM. Selective feeding by the giant barrel sponge enhances foraging efficiency. Limnol Oceanogr. 2016;61(4):1271–86.
EPA US. Methods for the determination of chemical substances in marine and estuarine environmental matrices. 2nd ed. Washington, DC: Edited by Agency USEP; 1997.
Caporaso JG, Lauber CL, Walters WA, Berg-Lyons D, Lozupone CA, Turnbaugh PJ, Fierer N, Knight R. Global patterns of 16S rRNA diversity at a depth of millions of sequences per sample. PNAS. 2011;108(Supplement 1):4516–22.
Schloss PD, Westcott SL, Ryabin T, Hall JR, Hartmann M, Hollister EB, Lesniewski RA, Oakley BB, Parks DH, Robinson CJ, et al. Introducing mothur: open-source, platform-independent, community-supported software for describing and comparing microbial communities. Appl Environ Microbiol. 2009;75(23):7537–41.
White JR, Nagarajan N, Pop M. Statistical methods for detecting differentially abundant features in clinical metagenomic samples. PLoS Comput Biol. 2009;5(4):e1000352.
Francis CA, Roberts KJ, Beman JM, Santoro AE, Oakley BB. Ubiquity and diversity of ammonia-oxidizing archaea in water columns and sediments of the ocean. PNAS. 2005;102(41):14683–8.
Calvó L, Garcia-Gil LJ. Use of amoB as a new molecular marker for ammonia-oxidizing bacteria. J Microbiol Methods. 2004;57(1):69–78.
Mohamed NM, Colman AS, Tal Y, Hill RT. Diversity and expression of nitrogen fixation genes in bacterial symbionts of marine sponges. Environ Microbiol. 2008;10(11):2910–21.
Yang Z, Li Z. Spatial distribution of prokaryotic symbionts and ammoxidation, denitrifier bacteria in marine sponge Astrosclera willeyana. Sci Rep. 2012;2:528.
Gantt SE, López-Legentil S, Erwin PM. Stable microbial communities in the sponge Crambe crambe from inside and outside a polluted Mediterranean harbor. FEMS Microbiol Lett. 2017;364(11):fnx105.
Kamke J, Sczyrba A, Ivanova N, Schwientek P, Rinke C, Mavromatis K, Woyke T, Hentschel U. Single-cell genomics reveals complex carbohydrate degradation patterns in poribacterial symbionts of marine sponges. ISME J. 2013;7:2287.
Astudillo-García C, Slaby BM, Waite DW, Bayer K, Hentschel U, Taylor MW. Phylogeny and genomics of SAUL, an enigmatic bacterial lineage frequently associated with marine sponges. Environ Microbiol. 2018;20(2):561–76.
Tkacz A, Hortala M, Poole PS. Absolute quantitation of microbiota abundance in environmental samples. Microbiome. 2018;6(1):110.
Weisz JB. Measuring impacts of associated microbial communities on Caribbean reef sponges : searching for symbiosis. Chapel Hill: University of North Carolina; 2006.
Park S-J, Ghai R, Martín-Cuadrado A-B, Rodríguez-Valera F, Chung W-H, Kwon K, Lee J-H, Madsen EL, Rhee S-K. Genomes of two new ammonia-oxidizing archaea enriched from deep marine sediments. PLoS One. 2014;9(5):e96449.
Poretsky RS, Sun S, Mou X, Moran MA. Transporter genes expressed by coastal bacterioplankton in response to dissolved organic carbon. Environ Microbiol. 2010;12(3):616–27.
Arrieta JM, Mayol E, Hansman RL, Herndl GJ, Dittmar T, Duarte CM. Dilution limits dissolved organic carbon utilization in the deep ocean. Science. 2015;348(6232):331–3.
Fiore CL, Jarett JK, Lesser MP. Symbiotic prokaryotic communities from different populations of the giant barrel sponge, Xestospongia muta. Microbiol Open. 2013;2(6):938–52.
de Goeij JM, Moodley L, Houtekamer M, Carballeira NM, Van Duyl FC. Tracing 13C-enriched dissolved and particulate organic carbon in the bacteria-containing coral reef sponge Halisarca caerulea: evidence for DOM-feeding. Limnol Oceanogr. 2008;53(4):1376–86.
Hoffmann F, Radax R, Woebken D, Holtappels M, Lavik G, Rapp HT, Schläppy M-L, Schleper C, Kuypers MM. Complex nitrogen cycling in the sponge Geodia barretti. Environ Microbiol. 2009;11(9):2228–43.
Schläppy M-L, Schöttner SI, Lavik G, Kuypers MM, de Beer D, Hoffmann F. Evidence of nitrification and denitrification in high and low microbial abundance sponges. Mar Biol. 2010;157(3):593–602.
Fiore CL, Baker DM, Lesser MP. Nitrogen biogeochemistry in the Caribbean sponge, Xestospongia muta: a source or sink of dissolved inorganic nitrogen? PLoS One. 2013;8(8):e72961.
We thank the staff of FIU’s Aquarius Reef Base in Key Largo, FL, and the Smithsonian Institution’s Carrie Bow Cay Field Station in Belize for logistical support. We thank R. Whitehead for the help with seawater sample analysis and current and previous members of the Pawlik, Erwin, and Finelli labs for the help with fieldwork and sample processing. Research in the Florida Keys National Marine Sanctuary was performed under permit FKNMS-2015-165; research conducted in Belize was performed under permit 000042-16 from the Belize Fisheries Department.
The research was funded by a National Science Foundation, Biological Oceanography Program grant (OCE-1558580).
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.
Sponge volumes, pumping rates, nutrient uptake, and nutrient flux data from Belize sponge holobionts used in DistLM analyses. (DOCX 19 kb)
Comparison of DistLM results using Bray-Curtis similarity based on square-root transformed vs. untransformed (raw) data, showing analysis for all species and by category (HMA/LMA). (DOCX 15 kb)
Comparison of DistLM results using Bray-Curtis similarity based on square-root transformed vs. untransformed (raw) data, showing analysis for all species and by category (HMA/LMA). (DOCX 15 kb)
Comparison of DistLM results using POC and DOC flux (specific filtration rates, SFR) vs. POC and DOC uptake (Cin – Cex, In-Ex), showing analysis for all species and by category (HMA/LMA). (DOCX 14 kb)
Comparison of DistLM results using nutrient flux (specific filtration rates, SFR) vs. nutrient uptake (Cin – Cex, In-Ex), showing analysis for all species and by category (HMA/LMA). (DOCX 15 kb)
Field (morphology-based) species identifications and top BLASTn species identifications from amplified barcoding sequences. (DOCX 14 kb)
Bray-Curtis dissimilarity matrices based 2D NMDs cluster plot of microbial communities for all samples (sponge and seawater) across locations (Belize and Florida). (DOCX 78 kb)
Taxonomic classifications and top BLASTn matches (sources in parenthesis, sponge species names in italics) of common microbial OTUs between locations (Belize and Florida). (DOCX 38 kb)
Spearman (rank order) correlations between carbon (DOC, POC) and nutrient flux (NH4, NOx, and PO4 specific filtration rates, μmol/s/L sponge) and microbial community (alpha) diversity metrics, showing correlation coefficients (R) and significance (P) values. (DOCX 20 kb)
Highly correlated microbial abundances and NH4 specific filtration rates for all sponges with more than 2 individuals sampled in Belize with a significant P-value (< 0.05). (DOCX 33 kb)
Significant differences in correlated microbial abundance and POC specific filtration rates for all sponges sampled at both sites that had more than 2 individuals sampled at each site and a significant P-value (< 0.05). (DOCX 56 kb)
About this article
Cite this article
Gantt, S.E., McMurray, S.E., Stubler, A.D. et al. Testing the relationship between microbiome composition and flux of carbon and nutrients in Caribbean coral reef sponges. Microbiome 7, 124 (2019). https://doi.org/10.1186/s40168-019-0739-x