Abstract
Deep-sea sponge grounds are underexplored ecosystems that provide numerous goods and services to the functioning of the deep-sea. This study assessed the prokaryotic diversity in embryos, recruits, and adults of Craniella zetlandica and Craniella infrequens, common and abundant representatives of deep-sea sponge grounds in the North Atlantic. Our results reveal that symbiont transmission in the two Craniella sponge species likely occurs vertically, as highly similar microbial consortia have been identified in adults, embryos, and recruits. Moreover, transmission electron microscopy revealed high abundances of sponge-associated microorganisms, among which Chloroflexi (SAR202) were identified as common representatives by amplicon sequencing and fluorescence in situ hybridization (FISH). Equal diversity metrices, a similar overall prokaryotic community composition and a distinct dominance of the phylum Chloroflexi within all life stages are the key findings of our analyses. Information such as presented here provide understanding on the recruitment of deep-sea sponge holobionts which is needed to develop integrated management tools of such vulnerable marine ecosystems.
Introduction
A fundamental process influencing the fitness of many organisms is the establishment of host-microbial colonization during development (). Microbial consortia enable the respective animal host to expand its range of nutrition (), defense (), immunity (), and development (). Host-microbe interactions may be regulated from early developmental stages and transmitted through generations (reviewed in ; ). In general, beneficial microbial symbionts are expected to be transmitted vertically (from parents to offspring) in order to assure an establishment of a compatible partnership ( and references therein). Further, probability of vertical transmission is expected to increase with increasing dependence of the host on its microbial partners (; ; Yamamura, 1993; ; Thompson, 1994; ; Wilkinson, 2001; ). In the marine environment symbioses are common and host-microbe interactions have been assessed at different developmental stages for a variety of model organisms, including cnidarians (; ), echinoderms (; ), cephalopods (), molluscs (; Wentrup et al., 2014), and annelids (; ).
Sponges (Porifera) are evolutionarily ancient animals which first appeared between 650 and 700 million years ago (; Zumberge et al., 2018). Today, more than 9,000 species (according to the World Porifera Database) are classified into four classes (Demospongiae, Calcarea, Homoscleromorpha, and Hexactinellida) and populate virtually all marine benthic habitats from shallow tropical reefs to the polar seas and the deep ocean (van Soest et al., 2020). Sponges are prominent examples of marine holobionts, consisting of diverse prokaryotic communities with more than 63 prokaryotic phyla detected within sponge hosts (; ). The prokaryotic communities of sponges are host species-specific and distinct from seawater prokaryotic communities (; ). The transmission of microbial sponge symbionts to the next generation occurs vertically or horizontally (; ; ; Webster et al., 2010; ). Microbial communities have been detected in reproductive stages, such as oocytes, embryos, larvae, and juveniles of various shallow water species (; ; ; ; ) indicating that vertical transmission is an important mechanism for the establishment and maintenance of sponge-microbe associations over evolutionary time scales.
In comparison to sponges from shallow and often warm-water habitats, much less is known about sponges and their microbiomes from deep and cold environments. Sponges are ubiquitous and abundant components of deep-sea benthic communities and can form structurally complex and diverse ecosystems, termed “sponge grounds” (; Xavier et al., 2015; ). Sponge grounds are known to strongly support biodiversity, either directly by provision of shelter or indirectly via food webs, and further perform a variety of ecosystem functions, among others via participation in bentho-pelagic coupling and biogeochemical cycles (reviewed in ; Thompson and Fuller, 2020). Deep-sea sponge grounds (syn. “sponge aggregations” or “sponge gardens”) have been flagged as vulnerable marine ecosystems (VMEs) that are currently under pressure through fishing and other human activities (; ). The species Craniella infrequens and Craniella zetlandica are common members of such sponge grounds and have further been recorded in shelf-areas around the globe (Figure 1). Craniella sponges are spherical in shape and grow up to the size of a golf ball and larger. This study aims to describe the prokaryotic community composition and diversity of two Craniella species and their early life stages from the deep North Atlantic Ocean. In particular we sought to address whether prokaryotic community composition, diversity, and dominant microbial phyla differ between different life stages (adults, embryos, and recruits) of the two sponge species. Knowledge about deep-sea sponge holobiont recruitment and plasticity across life stages is crucial to enable, in the long term, an adequate conservation of sponge populations from vulnerable marine ecosystems.
FIGURE 1
Materials and Methods
Sampling
Sponges of different developmental stages (adults, embryos, recruits) were collected during two research cruises onboard the Norwegian research vessel RV G.O.Sars (“GS2017110” and “GS2018108”) to the North Atlantic Ocean (Schulz Bank seamount at 73.8°N, 7.5°E; and Stjernsund at 70.3°N, 22.5°E) in July 2017 and August 2018. Samples were collected by remotely operated vehicle (ROV AEgir, University of Bergen) or Agassiz trawl at an average depth of 600 m below surface. Tissue chunks of adult sponges and recruits were subsampled onboard ship with sterile scalpels from the mesohyl, rinsed, flash-frozen and stored at −80°C for molecular analyses. Additional pieces of tissue were subsampled with biopsy punches and preserved onboard ship for microscopy. Embryos were carefully picked out of the parent sponges (Figure 2) with sterile spring steel forceps using a stereomicroscope. After rinsing, 10 whole embryos were each pooled per adult sponge to account for the small biomass, fixed in RNAlaterTM Stabilizing Solution (Thermo Fisher Scientific, Waltham, MA, United States), stored at 4°C overnight, transferred to −20°C on the next day and subsequently to −80°C. In addition, intact embryos were preserved onboard ship for microscopy. In terms of sample numbers, we sampled eight adult C. infrequens sponges, of which four were brooding embryos, and nine C. infrequens recruits. During the cruise we further sampled eight adult C. zetlandica sponges of which one was brooding embryos. In addition to our cruise activities, we sampled nine C. zetlandica recruits in September 2018 from an aquarium system (Bergen, Norway). The respective parent sponges had been maintained since summer 2017 in 45 L tanks with sea-water flow through (1 L min–1) pumped from the adjacent fjord at 105 m depth below surface. Seawater was treated with a 20 μm drum filter prior to flowing through the holding tanks. Under these conditions spawning was documented frequently and the new recruits had settled in the tanks. C. zetlandica recruits from the aquarium system were smaller (and thus most likely younger) than the in situ sampled C. infrequens recruits and therefore flash-frozen as a whole. In our case the new sponge recruits had the following sizes: 1 cm mean diameter for C. infrequens; and 0.3 cm maximal diameter for C. zetlandica. Adult sponges of both sponge species had the size of a fist. In summary, for this study in total 39 sponge individuals of the two sponge species were collected: C. zetlandica (eight adults, of which one was brooding embryos, and nine recruits) and C. infrequens (eight adults, of which four were brooding embryos, and nine recruits). In addition to the sponge samples, we also sampled 18 seawater microbial communities during the two cruises (at all sponge sampling locations) and from the aquarium system as environmental reference samples. Two liters of seawater sample were filtered onto polyvinylidene fluoride (PVDF) filter membranes (with 0.22 μm pore size; Merck Millipore) and stored at −80°C.
FIGURE 2
16S Amplicon Sequencing
DNA extraction was performed on either ∼ 0.25 g of tissue or bulk of 10 pooled embryos using the DNeasy Power Soil kit (Qiagen, Venlo, Netherlands). For seawater samples DNA was extracted of one half of a PVDF filter membrane. Subsequent DNA quantity and quality controls were performed with a NanoDrop spectrophotometer and by polymerase chain reaction with universal 16S primers (27F-1429R; Weisburg et al., 1991), followed by gel electrophoresis. For amplicon sequencing, the V3-V4 variable regions of the 16S rRNA gene were amplified in a one-step PCR with the primer pair 341F-806R (dual barcoding approach; ). The primer pair has the following nucleotide sequences: 5′-CCTACGGGAGGCAGCAG-3′ () and 5′-GGACTACHVGGGTWTCTAAT-3′ (). Verification of PCR-products was performed with gel electrophoresis, followed by normalization and pooling of samples. Sequencing was conducted on a MiSeq platform (MiSeqFGx, Illumina, San Diego, CA, United States) using v3 chemistry. Zero mismatches were allowed in the barcode sequence for demultiplexing.
Raw reads were quality-checked and processed within the QIIME2 environment (version 2018.11; ). Based on forward reads (truncated to 270 nt), Amplicon Sequence Variants (ASVs) were generated using the DADA2 algorithm (python-based; ). Data were analyzed based on ASVs instead of classical operational taxonomic units (OTUs), in order to resolve fine-scale community patterns and to use a state-of-the-art method which profits from the error correction model conducted by the DADA2 algorithm. Singletons and chimeric sequences were removed during this process.
Phylogenetic trees were calculated based on the resulting ASVs via the FastTree2 plugin. Taxonomic classification of representative ASVs was achieved with the help of a primer-specific trained Naive Bayes taxonomic classifier, based on the Silva 132 99% OTUs 16S database (). Taxonomic filtering was applied on the dataset, meaning that mitochondria, chloroplasts and unassigned sequences were removed. Further, a sampling depth of 3,900 was applied to standardize the number of reads across samples. Diversity metrices were calculated, such as Shannon indices (alpha diversity) and weighted UniFrac distances (beta diversity). Non-metric multidimensional scaling (nMDS) was conducted to evaluate sample separation in ordination space (based on weighted UniFrac distances). We performed PERMANOVAs (999 permutations) on weighted UniFrac distances, and Tukey Posthoc tests on Shannon indices following ANOVA type III tests, to assess if the prokaryotic community differed significantly between sample groups. To evaluate further on the differential abundance of microbial phyla present in aquarium sponges (i.e., C. zetlandica recruits) in comparison to in situ sampled sponges (and in aquarium water vs in situ water), we ran the Linear Discriminant Analysis (LDA) Effect Size (LEfSe) algorithm (). Visualizations were created using R (version 3.0.2 with the ggplot2 package loaded; ; Wickham, 2016) and Inkscape (version 0.92.4; ). The ternary plot was created using the Ternary plot maker web tool (1).
In order to provide an independent proof-of-concept accompanying our 16S amplicon sequencing results we performed qualitative microscopic analyses. These included light microscopy, transmission electron microscopy (TEM), and fluorescence in situ hybridization of the most dominant microbial phylum detected by amplicon sequencing. Especially for the embryo-based DNA libraries we aimed at confirming by microscopy that they contained sequences which were indeed present within the embryos and not merely contaminants.
Light Microscopy and TEM
Tissue samples for light microscopy and TEM were fixed in a cocktail of glutaraldehyde (2.5%) with sodium chloride (0.34 M) and phosphate buffered saline (PBS; 0.4 M). Samples were stored at 4°C. Back on land, samples were washed 3x with buffer (at 4°C; 15 min per washing step), post-fixed in 2% osmiumtetroxide for 2 h and afterward washed three more times with buffer (at 4°C; 15 min per washing step). Subsequently partial dehydration was conducted with an ascending ethanol series, from 30% EtOH (2 × 15 min) to 50% EtOH (1 × 15 min), to 70% EtOH (2 × 15 min, including an overnight storage at 4°C), to 90% EtOH (1 × 15 min) up to 100% EtOH (2 × 15 min). After dehydration, the samples were gradually infiltrated with LR-White resin (at room temperature). The infiltration series covered the following steps: from 1 × 1 h 2:1 EtOH:LR-White, to 1 × 1 h 1:1 EtOH:LR-White, to 1 × 1 h 1:2 EtOH:LR-White, up to 2 × 2 h pure LR-White. Samples were transferred to gelatine embedding capsules filled with fresh LR-White resin and subsequently incubated for polymerization at 57°C for two days. The resulting resin blocks were manually pre-trimmed. Semithin sections (0.5 μm) were cut for light microscopy using an ultramicrotome (Reichert-Jung ULTRACUT E; Leica Camera AG, Wetzlar, Germany) which was equipped with a diamond knife (DIATOME, Biel, Switzerland). After staining with a Richardson solution (after ), semi-thin sections were mounted onto SuperFrost Ultra Plus® microscopy slides and visualized with the help of an Axio Observer.Z1 microscope (Zeiss, Göttingen, Germany).
For TEM, ultrathin sections (70 nm) were cut with the ultramicrotome, mounted on pioloform coated copper grids (75 mesh) and contrasted with uranyl acetate (20 min incubation and subsequent washing steps) as well as with Reynold’s lead citrate (3 min incubation, followed by washing steps). Visualization of the ultra-thin sections was performed on a Tecnai G2 Spirit Bio Twin transmission electron microscope (FEI Company, Hillsboro, United States) at an acceleration voltage of 80 kV.
Fluorescence in situ Hybridization
Sponge tissue samples were fixed overnight (at 4°C) onboard with 4% paraformaldehyde (PFA) in 1xPBS (pH was adjusted to approximately 6.9 with the help of NaOH and HCl). For preparation of all FISH solutions, we used filtered (0.22 μm membrane) diethylpyrocarbonate (DEPC)-treated water (Thermo Fisher Scientific, Waltham, United States) to inactivate RNase enzymes, and performed autoclaving for sterilization of buffers and solvents onboard the ship. After fixation with PFA, five washing steps (3 min) were performed with 1xPBS (4°C). Subsequently quenching of free aldehyde groups was conducted by incubation with glycine (50 mM) in 1xPBS for 15 min (room temperature). The quenching step was followed by five additional washing steps (3 min) in DEPC water (at 4°C). Dehydration was performed onboard in a ascending ethanol series (molecular grade, RNase free), covering the following steps: 30% EtOH in 1xPBS (2 × 15 min at 4°C), 50% EtOH in 1xPBS (2 × 30 min). Samples were subsequently stored at −20°C. Back on land, dehydration was continued in following steps: 50% EtOH in MilliQ water (1 × 15 min), 60% EtOH in MilliQ water (2 × 15 min), 70% EtOH in MilliQ water (2 × 15 min), 90% EtOH in MilliQ water (1 × 15 min) up to 100% EtOH (2 × 15 min). The samples were gradually infiltrated with LR-White resin, polymerized, pre-trimmed and cut into semithin sections (0.5 μm) with an ultramicrotome, in the same manner as described in section 2.3. As embryos were embedded in whole, cutting with the microtome was run for a substantial time, until the central part of the embryo was reached. At the area with widest diameter, sections were received and fixed onto two separate SuperFrost Ultra Plus® microscopy slides. While one half of the slides was directly stained with a Richardson solution and mounted using Biomount medium (Plano, Germany; see description in section “Light Microscopy and TEM”). The other half was only fixed onto microscopy slides, without any staining or covering in Biomount medium. Following the procedures described in , these latter samples were afterward labeled with Chloroflexi clade-specific probes. In particular the following probes were used: Cal825 (5′-[Cy3]-ACACCGCCCACACCTCGT-[Cy3]-3′; E. coli binding positions 825 to 843) for Caldilineae; Ana1005 (5′-[Alexa Fluor 647]-TCCGCTTTCGCTTCCGTA-[Alexa Fluor 647]-3′; E. coli binding positions 1005 to 1023) for Anaerolineae; and SAR202-104 (5′-[Alexa Fluor 488]-GTTACTCAGCCGTCTGCC-[Alexa Fluor 488]-3′; E. coli binding positions 104 to 122) for the SAR202 group. A double labeling approach at 5′ and 3′ ends was performed for all probes (Sigma-Aldrich, Germany).
Hybridization was performed in Sylgard chambers (GEOMAR in-house production) which were located inside an equilibrated humid chamber (at 46°C for 3 h). The hybridization buffer consisted of 900 mM NaCl, 20 mM Tris-HCl (pH = 7.4), 10% formamide, 0.01% sodium dodecyl sulfate and 20% dextransulfate. All Chloroflexi clade-specific probes were co-hybridized, in addition counterstaining of bacterial nucleic acids and sponge cell nuclei was performed with pre-warmed DAPI in hybridization buffer (1 ng/μL), with an incubation time of 20 min at 46°C. After hybridization, the slides were incubated in pre-warmed wash buffer (450 mM NaCl, 20 mM Tris–HCl, 0.01% sodium dodecyl sulfate) at 48°C for 25 min. The washing step was repeated twice. Subsequently, the slides were carefully rinsed with ice cold molecular grade water, dried and mounted in Mowiol medium (Mowiol® 40–88, Kuraray Europe GmbH, Tokyo, Japan). The same procedure was conducted for a negative control for non-specific binding (i.e., hybridization buffer without Chloroflexi-specific probes, but with an oligonucleotide complementary to the probe EUB338 (NON-338-probe; Wallner et al., 1993), to prevent falling into the trap of false-positive results. Fluorescence signals at the four specific wavelength ranges were detected using an Axio Observer.Z1 microscope equipped with AxioCam 506 and Zen 2 software (version 2.0.0.0; Zeiss Germany).
Results
In this study, the prokaryotic community composition of 39 Craniella sponge individuals was analyzed. A visual overview of the denoising stats is compiled in Supplementary Figure S1. The complete dataset after quality filtering consisted of 614,928 reads, which is 82% of the initial input read count. On average each sample had 15,767 reads, but as several samples contained less reads, a sampling depth of 3,900 was applied on the complete dataset for normalization. Four adult C. infrequens individuals and one adult C. zetlandica individual were found brooding embryos inside their tissue (Figure 2). For each of these brooding individuals, ten embryos had been picked and pooled. In C. infrequens, the embryos were more numerous compared to C. zetlandica and were scattered throughout the mesohyl of the parental sponge. Prokaryotic diversity (as estimated by Shannon index) was similar between the developmental stages within each sponge species (Figure 3 and Supplementary Table S1). Other alpha diversity indices (Pielou’s eveness, Faith’s PD and number of observed ASVs) led to a similar result (Supplementary Figure S2). Average Shannon indices were slighly lower in C. zetlandica than in C. infrequens, despite not statistically significantly (Figure 3 and Supplementary Table S1). The prokaryotic community compositions of both sponge species were significantly different from seawater (Figure 4A and Supplementary Table S2). Further, C. infrequens and C. zetlandica displayed species-specific microbiomes, as shown by their distinct clustering based on weighted UniFrac distances and on PERMANOVA analyses (Figure 4B, with some exceptions; Supplementary Table S2). Overall C. zetlandica samples showed a slightly higher variability than C. infrequens samples (please note in this regard that while C. infrequens specimens were all sampled at the same geographic location, C. zetlandica specimens analyzed in this study originate from different geographic locations).
FIGURE 3
FIGURE 4
With respect to the prokaryotic community composition, most bacterial phyla were shared across all life stages, and were not particularly enriched in any specific life stage (Figure 5A). Chloroflexi were the most abundant prokaryotic phylum in both sponge species and across all life stages and SAR202 was the most abundant clade, constituting on average 73% of all recorded Chloroflexi (Table 1). C. zetlandica recruits obtained from the aquarium system were a noticeable exception that were dominated by Proteobacteria and significantly enriched in this phylum compared to in situ sampled sponges (Figure 5B and Supplementary Figure S3A). Proteobacteria were however not significantly enriched in the aquarium water compared to in situ sampled seawater (Supplementary Figures S3B C). Neither were Chloroflexi significantly depleted in the aquarium water compared to in situ sampled seawater, although Chloroflexi were significantly enriched in the in situ sampled sponges compared to the aquarium sampled sponges. A total of 390 Chloroflexi ASVs were detected in the 39 Craniella samples, of which 171 were present in ≥50% of the samples per category (adult, embryo, recruit) (Figure 6A). Of these 171 common Chloroflexi ASVs, 16 were present in all life stages and both species, while being absent from seawater reference samples. These 16 ASVs are particularly abundant and account for 20–40% of the total Chloroflexi prokaryotic community in each sponge sample. Seven ASVs of the 16 ASVs belong to the SAR202 clade, while the remainder belongs to the classes JG30-KF-CM66 (4 ASVs), Anaerolineae (2 ASVs), TK17 (2 ASVs), and the order S085 (1 ASV).
FIGURE 5
TABLE 1
| Anaerolineales | Caldilineales | SAR202 clade | Other Chloroflexi | ||
| C. zetlandica | Adult | 0.1 ± 0.1 | 6.9 ± 3.4 | 77.5 ± 8.2 | 15.5 ± 7.3 |
| Embryo | 0 ± NA | 7.9 ± NA | 71.4 ± NA | 20.7 ± NA | |
| Recruit | 0 ± 0 | 15.6 ± 6.0 | 52.0 ± 14.2 | 32.4 ± 9.1 | |
| C. infrequens | Adult | 0.2 ± 0.1 | 6.0 ± 4.1 | 85.2 ± 4.4 | 8.7 ± 2.0 |
| Embryo | 0.2 ± 0.2 | 7.1 ± 0.8 | 78.9 ± 3.7 | 13.7 ± 3.5 | |
| Recruit | 0.5 ± 0.7 | 11.2 ± 6.1 | 70.0 ± 11.5 | 18.2 ± 8.1 |
Average relative abundances of Chloroflexi (Anaerolinae, Caldilineae, SAR202, others) as determined by 16S amplicon sequencing. mean, and standard deviation per group are indicated.
FIGURE 6
We next explored whether the 171 common Chloroflexi ASVs are species-specific or developmental stage-specific (Figures 6B,C). We determined 46 ASVs as being C. zetlandica specific (absent from C. infrequens), and 73 ASVs as being C. infrequens specific (absent from C. zetlandica). For C. zetlandica 26 ASVs were exclusively occurring in adults, 12 in embryos, and 10 in recruits. For C. infrequens four ASVs were occurring only in adults, 18 in embryos, and four in recruits. Interestingly, one Chloroflexi-ASV, occurred only in embryos of both sponge species. This ASV was affiliated with the Chloroflexi-class TK17, that could only be classified to class-level. Upon blasting against a sponge microbiome reference database (), this particular ASV was found in 134 sponges samples, of which 88 samples were identified as Geodia barretti. In summary, the major fraction of the detected Chloroflexi ASVs were sponge species-specific but not developmental-stage specific. Further information on the relative abundances of the ASV groups is presented in Supplementary Figure S4.
Light microscopy revealed that bacterial cells were numerous inside C. infrequens. In the embryos, the microbial cells were restricted to a narrow band at the peripheral rim (Figures 7A,B), while the interior was devoid of cellular material or bacterial morphotypes. An artifact resulting from cutting through the possibly more fluidic embryo interior can however not be excluded. Fluorescence in situ hybridization confirmed high abundances of Chloroflexi, in particular SAR202, within the “bacterial band” at the embryo rim (Figure 7C). The negative control for non-specific binding was negative (Supplementary Figures S5I–L). Bacterial cells were further found to be abundant and widely distributed in the mesohyl in adult C. infrequens (Supplementary Figure S6A). FISH confirmed high abundances of SAR202 also for adult individuals (Supplementary Figure S6B). Microbial cells within the embryos displayed diverse morphologies and were found in dividing stages (see Figure 7D for examples).
FIGURE 7
For C. zetlandica we observed an in situ (stress-induced) release of ripe propagules, visible by ROV live-footage in August 2018, during the sampling campaign “GS2018108” in Stjernsund (Supplementary Video S1: doi.pangaea.de/10.1594/PANGAEA.918218). We suspect a similar timing of spawning in summer for those sponges maintained in the aquarium system in Bergen, as high numbers of settled recruits were found in October 2018 (Supplementary Figure S7). These observations point to a synchronized reproductive cycle in aquarium and in situ thriving C. zetlandica with a spawning in late summer.
Discussion
Deep-sea sponge grounds are underexplored ecosystems that provide numerous goods and services to the functioning of the deep-sea (), such as the provision of food, habitat and refugia for commercially relevant fish species. The need to preserve these vulnerable marine ecosystems (VMEs) is thus contradictory to the high demand for socio-economic processes such as fishing in North Atlantic waters. In this study we assessed the prokaryotic diversity of two adult deep-sea sponge species of the genus Craniella and their developmental stages. These sponges are prominent although understudied members of deep-sea sponge grounds and have further been recorded in shelf-areas around the globe. While there is one report on the Proteobacteria-dominated microbiome of Craniella australiensis, which is a shallow water species in the Indo-Pacific and South China Sea (), the microbiome of C. infrequens and C. zetlandica have—to our knowledge—not been characterized to date.
Our results reveal that C. zetlandica and C. infrequens are typical high microbial abundance (HMA) sponges (; ). We base this statement on our observation that both sponge species harbor diverse and dense microbial communities within their tissues, as observed by TEM. This finding is further supported by the observation that the sponge microbiomes are dominated by Chloroflexi (SAR202), which are known indicator taxa of HMA sponges ().
The microbiomes of C. zetlandica and C. infrequens were host-species specific and distinctly different from seawater which is consistent with previous reports on sponge microbiomes (; ; ). The prokaryotic community composition and alpha diversity were very similar between adults, embryos and recruits within the same sponge species, indicating that entire microbial consortia are very likely vertically transmitted. Vertical transmission has been previously reported for a number of shallow-water HMA sponges (i.e., ; ), but is, to the best of our knowledge, novel in the context of sponge species from sponge grounds. Our results indicate that young deep-sea sponges can already contain a high diversity and density of microorganisms. Despite the observation of very similar prokaryotic community compositions between all examined life stages per species in our study, the recruits of C. zetlandica fell somewhat distant from all other samples in the nMDS plots and were dominated by Proteobacteria. It is noteworthy that these samples were collected in an aquarium system and not in the field. Indeed, previous studies have shown that the maintenance of marine organisms in captivity can affect the relative abundances of some associated microbial taxa (). In addition, microbial variability has been observed to be particularly high in recently settled recruits compared to adult individuals (; ).
The most dominant microbial phylum in all in situ sampled sponges of both sponge species and across all developmental stages were Chloroflexi. Generally, the probability of vertical symbiont-transmission is assumed to be positively related with the dependence of the host on the microbial partner partners (; ; Yamamura, 1993; ; Thompson, 1994; ; Wilkinson, 2001; ). Bacteria of the phylum Chloroflexi (in particular SAR202) are abundant and wide-spread in HMA sponges (; ) and are suggested to play an important function in the degradation of recalcitrant organic matter (; ). We propose that the 16 Craniella core ASVs detected within our study should be considered as sponge symbionts as they were found in high abundances across both species and all developmental stages, while being absent from seawater. Almost half of the 16 Craniella core ASVs belong to the SAR202 clade, making it tempting to speculate about a potential crucial functional role in respect to (recalcitrant) organic matter cycling. Our FISH analyses revealed high abundances of Chloroflexi scattered within the bacterial band at the inner periphery of sponge embryos. We conclude that the transmission of Chloroflexi through the embryonal stages reflects their tight and evolutionarily long-standing integration in the sponge holobiont. Future research efforts should be geared toward understanding the specific functional roles of Chloroflexi symbionts for deep-sea sponges.
Statements
Data availability statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found below: https://www.ncbi.nlm.nih.gov/, PRJNA612549.
Author contributions
KBu and UH designed the study, wrote the manuscript, and were responsible for the microbial pipeline. KBu, HR, and EW conducted sampling. KBa and KBu performed fluorescence in situ hybridization. HR led sponge taxonomic analysis. AF was involved in sequencing. KBu performed the data analysis (bioinformatics and visualizations). EW and KBa reviewed and edited the manuscript. All authors contributed to the article and approved the submitted version.
Funding
This study was funded by the European Union’s Horizon 2020 Research and Innovation Program under Grant Agreement No. 679849 (the SponGES project). This document reflects only the authors’ view and the Executive Agency for Small and Medium-sized Enterprises (EASME) is not responsible for any use that may be made of the information it contains.
Acknowledgments
We express our gratitude to HR, coordinator of the H2020-SponGES project on deep-sea sponge grounds of the North Atlantic, deep-sea sponge taxonomist, explorer, mentor, and friend. We thank the crews and scientific parties of RV G.O.Sars cruises “GS2017110” and “GS2018108”. On these cruises the ROV AEgir-team (University of Bergen) enabled sampling of precious samples. Ina Clefsen and Andrea Hethke provided valuable support in the lab back on land. We also thank Heikki Savolainen and Tone Ulvatn for help in aquaria facility management and technical support at the University of Bergen. Heidi Kristina Meyer’s support in the acquisition of video material is much appreciated.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest. The handling Editor declared a past co-authorship with one of the authors HR.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmars.2020.00674/full#supplementary-material
Footnotes
References
1
AnbuchezhianR.LiZ. (2014). “Marine sponge Craniella australiensis-associated bacterial,” in Encyclopedia of Metagenomics, ed.NelsonK. E. (Berlin: Springer-Verlag), 319–326. 10.1007/978-1-4899-7475-4_604
2
ApprillA.MarlowH. Q.MartindaleM. Q.RappéM. S. (2009). The onset of microbial associations in the coral Pocillopora meandrina.ISME J.3685–699. 10.1038/ismej.2009.3
3
BayerK.JahnM. T.SlabyB. M.Moitinho-SilvaL.HentschelU. (2018). Marine sponges as Chloroflexi hot spots: genomic insights and high-resolution visualization of an abundant and diverse symbiotic clade.mSystems3:e00150-18.
4
BjörkJ. R.Díez-VivesC.Astudillo-GarcíaC.ArchieE. A.MontoyaJ. M. (2019). Vertical transmission of sponge microbiota is inconsistent and unfaithful.Nat Ecol Evol.31172–1183. 10.1038/s41559-019-0935-x
5
BlomJ. (2019). Ternary Plot Maker. Available online at: https://www.ternaryplot.com/(accessed March 1, 2020).
6
BolyenE.RideoutJ. R.DillonM. R.BokulichN. A.AbnetC.Al GhalithG. A.et al (2019). Reproducible, interactive, scalable, and extensible microbiome data science using QIIME 2.Nat. Biotechnol.37852–857.
7
BullJ. J.MolineuxI. J.RiceW. R. (1991). Selection of benevolence in a host-parasite system.Evolution45875–882. 10.1111/j.1558-5646.1991.tb04356.x
8
BuschK.BeazleyL.KenchingtonE.WhoriskeyF.SlabyB.HentschelU. (2020). Microbial diversity of the glass sponge Vazella pourtalesii in response to anthropogenic activities.bioRxiv [Preprint]. 10.1101/2020.05.19.102806
9
CallahanB. J.McMurdieP. J.RosenM. J.HanA. W.JohnsonA. J. A.HolmesS. P. (2016). DADA2: high-resolution sample inference from Illumina amplicon data.Nat. Methods13581–583. 10.1038/nmeth.3869
10
CaporasoJ. G.LauberC. L.WaltersW. A.Berg-LyonsD.KnightR.LozuponeC. A.et al (2011). Global patterns of 16S rRNA diversity at a depth of millions of sequences per sample.Proc. Natl. Acad. Sci. U.S.A.1084516–4522. 10.1073/pnas.1000080107
11
CarrierT. J.DupontS.ReitzelA. M. (2019). Geographic location and food availability offer differing levels of influence on the bacterial communities associated with larval sea urchins.FEMS Microbiol. Ecol.951–9.
12
De CaraltS.UrizM. J.WijffelsR. H. (2007). Vertical transmission and successive location of symbiotic bacteria during embryo development and larva formation in Corticium candelabrum (Porifera: Demospongiae).J. Mar. Biol. Assoc. U. K.871693–1699. 10.1017/s0025315407056846
13
DouglasA. E. (1994). Symbiotic Interactions.Oxford: Oxford University Press.
14
EassonC. G.ThackerR. W. (2014). Phylogenetic signal in the community structure of host-specific microbiomes of tropical marine sponges.Front. Microbiol.5:532. 10.3389/fmicb.2014.00532
15
EberlG. (2010). A new vision of immunity: homeostasis of the superorganism.Mucosal. Immunol.3450–460. 10.1038/mi.2010.20
16
EndowK.OhtaS. (1990). Occurrence of bacteria in the primary oocytes of vesicomyid clam Calyptogena soyoae.Mar. Ecol. Prog. Ser.64309–311. 10.3354/meps064309
17
EreskovskyA. V.GonoboblevaE.VishnyakovA. (2005). Morphological evidence for vertical transmission of symbiotic bacteria in the viviparous sponge Halisarca dujardini Johnston (Porifera. Demospongiae, Halisarcida).Mar. Biol.146869–875. 10.1007/s00227-004-1489-1
18
EwaldP. W. (1987). Transmission modes and evolution of the parasitism-mutualism continuum.Ann. N. Y. Acad. Sci.503295–306. 10.1111/j.1749-6632.1987.tb40616.x
19
FlórezL. V.BiedermannP. H. W.EnglT.KaltenpothM. (2015). Defensive symbioses of animals with prokaryotic and eukaryotic microorganisms.Nat. Prod. Rep.32904–936. 10.1039/c5np00010f
20
GalacM. R.BoschI.JaniesD. A. (2016). Bacterial communities of oceanic sea star (Asteroidea: Echinodermata) larvae.Mar. Biol.163:162.
21
GalandP. E.ChapronL.MeistertzheimA.-L.PeruE.LartaudF. (2018). The effect of captivity on the dynamics of active bacterial communities differs between two deep-sea coral species.Front. Microb.9:2565. 10.3389/fmicb.2018.02565
22
GloecknerV.LindquistN.SchmittS.HentschelU. (2013). Ectyoplasia ferox, an experimentally tractable model for vertical microbial transmission in marine sponges.Microb. Ecol.65462–474. 10.1007/s00248-012-0142-7
23
GloecknerV.WehrlM.Moitinho-silvaL.SchuppP.PawlikJ. R.LindquistN. L.et al (2014). The HMA-LMA dichotomy revisited?: an electron microscopical survey of 56 sponge species.Biol. Bull.22778–88. 10.1086/bblv227n1p78
24
HarringtonB.the Inkscape Development Team (2005). Inkscape. Available online at: http://www.inkscape.org/(accessed February 18, 2019).
25
HentschelU.HopkeJ.HornM.FriedrichA. B.WagnerM.HackerJ.et al (2002). Molecular evidence for a uniform microbial community in sponges from different oceans.Appl. Environ. Microbiol. 684431–4440. 10.1128/AEM.68.9.4431-4440.2002
26
HerreE. A.KnowltonN.MuellerU. G.RehnerS. A. (1999). The evolution of mutualisms: exploring the paths between conflict and cooperation.Trends Ecol. Evol.1449–53. 10.1016/s0169-5347(98)01529-8
27
KlitgaardA. B.TendalO. S. (2004). Distribution and species composition of mass occurrences of large-sized sponges in the northeast Atlantic.Prog. Oceanogr.6157–98. 10.1016/j.pocean.2004.06.002
28
KoropatnickT. A.EngleJ. T.ApicellaM. A.StabbE. V.GoldmanW. E.Mcfall-ngaiM. J. (2004). Microbial factor-mediated development in a host-bacterial mutualism.Science3061186–1189.
29
KozichJ. J.WestcottS. L.BaxterN. T.HighlanderS. K.SchlossP. D. (2013). Development of a dual-index sequencing strategy and curation pipeline for analyzing amplicon sequence data on the Miseq Illumina sequencing platform.Appl. Environ. Microbiol.795112–5120. 10.1128/aem.01043-13
30
LandryZ.SwaB. K.HerndlG. J.StepanauskasR.GiovannoniS. J. (2017). SAR202 genomes from the dark ocean predict pathways for the oxidation of recalcitrant dissolved organic matter.mBio8e00413–17.
31
LoveG. D.GrosjeanE.StalviesC.FikeD. A.GrotzingerJ. P.BradleyA. S.et al (2009). Fossil steroids record the appearance of Demospongiae during the Cryogenian period.Nature457718–721. 10.1038/nature07673
32
MaldonadoM. (2007). Intergenerational transmission of symbiotic bacteria in oviparous and viviparous demosponges, with emphasis on intracytoplasmically-compartmented bacterial types.J. Mar. Biol. Assoc. U. K.871701–1713. 10.1017/s0025315407058080
33
MaldonadoM.AguilarR.BannisterR. J.BellJ. J.ConwayK. W.DaytonP. K.et al (2016). “Sponge grounds as key marine habitats: a synthetic review of types, structure, functional roles, and conservation concerns,” in Marine Animal Forests, edsRossiS.BramantiL.GoriA.OrejasC. (Berlin: Springer-Verlag), 1–39. 10.1007/978-3-319-17001-5_24-1
34
McFall-NgaiM.HadfieldM. G.BoschT. C. G.CareyH. V.Domazet-LošoT.DouglasA. E. (2013). Animals in a bacterial world, a new imperative for the life sciences.Proc. Natl. Acad. Sci. U.S.A.1103229–3236.
35
Moitinho-SilvaL.NielsenS.AmirA.GonzalezA.AckermannG. L.CerranoC.et al (2017a). The sponge microbiome project.Gigascience6:gix077.
36
Moitinho-SilvaL.SteinertG.NielsenS.HardoimC. C. P.WuY. C.McCormackG. P.et al (2017b). Predicting the HMA-LMA status in marine sponges by machine learning.Front. Microbiol.8:752. 10.3389/fmicb.2017.00752
37
MontgomeryM. K.McFall-NgaiM. (1994). Bacterial symbionts induce host organ morphogenesis during early postembryonic development of the squid Euprymna scolopes.Development1201719–1729.
38
MortzfeldB. M.UrbanskiS.ReitzelA. M.KünzelS.TechnauU.FrauneS. (2016). Response of bacterial colonization in Nematostella vectensis to development, environment and biogeography.Environ. Microbiol.181764–1781. 10.1111/1462-2920.12926
39
MulischM.WelschU. (2015). Romeis Mikroskopische Technik.Berlin: Springer, 1–603.
40
MurilloF. J.KenchingtonE.Koen-AlonsoM.GuijarroJ.KenchingtonT. J.SacauM.et al (2020). Mapping benthic ecological diversity and interactions with bottom-contact fishing on the Flemish Cap (northwest Atlantic).Ecol. Indic.112:106135. 10.1016/j.ecolind.2020.106135
41
MuyzerG.WaalE. C. D. E.UitierlindenA. G. (1993). Profiling of complex microbial populations by denaturing gradient gel electrophoresis analysis of polymerase chain reaction-amplified genes coding for 16S rRNA.Appl. Environ. Microbiol.591–6.
42
NicholsonJ. K.HolmesE.KinrossJ.BurcelinR.GibsonG.JiaW.et al (2012). Metabolic Interactions.Science1081262–1268.
43
NussbaumerA. D.FisherC. R.BrightM. (2006). Horizontal endosymbiont transmission in hydrothermal vent tubeworms.Nature441345–348. 10.1038/nature04793
44
PhamC. K.MurilloF. J.LiretteC.MaldonadoM.ColaçoA.OttavianiD.et al (2019). Removal of deep-sea sponges by bottom trawling in the Flemish Cap area: conservation, ecology and economic assessment.Sci. Rep.9:15843.
45
QuastC.PruesseE.YilmazP.GerkenJ.SchweerT.YarzaP.et al (2013). The SILVA ribosomal RNA gene database project: improved data processing and web-based tools.Nucleic Acids Res.41590–596.
46
R Development Core Team (2008). R: A Language and Environment for Statistical Computing.Vienna: R Foundation for Statistical Computing.
47
RosenbergE.Zilber-RosenbergI. (2011). Symbiosis and development: the hologenome concept.Birth Defects Res. Part C Embryo Today Rev.9356–66. 10.1002/bdrc.20196
48
SachsJ. L.SkophammerR. G.RegusJ. U. (2011). Evolutionary transitions in bacterial symbiosis.Proc. Natl. Acad. Sci. U.S.A.10810800–10807. 10.1073/pnas.1100304108
49
SchimakM. P.SchoffelenN. J.GruhlA.DubilierN. (2015). Transmission of Bacterial Symbionts in the Gutless Oligochaete Olavius Algarvensis.Bremen: University of Bremen, 161–199.
50
SchmittS.AngermeierH.SchillerR.LindquistN.HentschelU. (2008). Molecular microbial diversity survey of sponge reproductive stages and mechanistic insights into vertical transmission of microbial symbionts.Appl. Environ. Microbiol.747694–7708. 10.1128/aem.00878-08
51
SchmittS.WeiszJ. B.LindquistN.HentschelU. (2007). Vertical transmission of a phylogenetically complex microbial consortium in the viviparous sponge Ircinia felix.Appl. Environ. Microbiol.732067–2078. 10.1128/aem.01944-06
52
SegataN.IzardJ.WaldronL.GeversD.MiropolskyL.GarrettW. S.et al (2011). Metagenomic biomarker discovery and explanation.Genome Biol.12:R60.
53
SharpK. H.EamB.John FaulknerD.HaygoodM. G. (2007). Vertical transmission of diverse microbes in the tropical sponge Corticium sp.Appl. Environ. Microbiol.73622–629. 10.1128/aem.01493-06
54
SteinertG.RohdeS.JanussenD.BlaurockC.SchuppP. J. (2017). Host-specific assembly of sponge-associated prokaryotes at high taxonomic ranks.Sci. Rep.71–9.
55
ThomasT.Moitinho-SilvaL.LurgiM.BjörkJ. R.EassonC.Astudillo-GarcíaC.et al (2016). Diversity, structure and convergent evolution of the global sponge microbiome.Nat. Commun.7:11870.
56
ThompsonJ. N. (1994). The Coevolutionary Process.Chicago, IL: The University of Chicago Press, 1–383.
57
ThompsonT.FullerS. D. (2020). Technical Measures and Environmental Risk Assessments for Deep-Sea Sponge Conservation.Rome: FAO, 1–30.
58
van SoestR. W. M.Boury-EsnaultN.HooperJ. N. A.RützlerK.de VoogdN. J.AlvarezB.et al (2020). World Porifera Database. Accessed at: http://www.marinespecies.org/porifera(accessed February 06, 2020).
59
WallnerG.AmannR.BeiskerW. (1993). Optimizing fluorescent in situ hybridization with rRNA-targeted oligonucleotide probes for flow cytometric identification of microorganisms.Cytometry14136–143. 10.1002/cyto.990140205
60
WebsterN. S.TaylorM. W.BehnamF.LückerS.RatteiT.WhalanS.et al (2010). Deep sequencing reveals exceptional diversity and modes of transmission for bacterial sponge symbionts.Environ. Microbiol.122070–2082.
61
WeisburgW. G.BarnsS. M.PelletierD. A.LaneD. J. (1991). 16S ribosomal DNA amplification for phylogenetic study.J. Bacteriol.173697–703. 10.1128/jb.173.2.697-703.1991
62
WentrupC.WendebergA.SchimakM.BorowskiC.DubilierN. (2014). Forever competent: deep-sea bivalves are colonized by their chemosynthetic symbionts throughout their lifetime.Environ. Microbiol.163699–3713. 10.1111/1462-2920.12597
63
WickhamH. (2016). ggplot2: Elegant graphics for Data Analysis.New York, NY: Springer-Verlag, 1–213.
64
WilkinsonD. M. (2001). Horizontally acquired mutualisms, an unsolved problem in ecology??Oikos92377–384. 10.1034/j.1600-0706.2001.920222.x
65
XavierJ. R.CárdenasP.CristoboJ.Van SoestR.RappH. T. (2015). Systematics and biodiversity of deep-sea sponges of the Atlanto-Mediterranean region.J. Mar. Biol. Assoc. U. K.951285–1286.
66
YamamuraN. (1993). Vertical transmission and evolution of mutualism from parasitism.Theor. Popul. Biol.4495–109. 10.1006/tpbi.1993.1020
67
ZumbergeJ. A.LoveG. D.CárdenasP.SperlingE. A.GunasekeraS.RohrssenM.et al (2018). Demosponge steroid biomarker 26-methylstigmastane provides evidence for Neoproterozoic animals.Nat. Ecol. Evol.21709–1714. 10.1038/s41559-018-0676-2
Summary
Keywords
vulnerable marine ecosystems, sponges, Craniella, Chloroflexi, symbiosis, early life stages, amplicon sequencing, fluorescence in situ hybridization
Citation
Busch K, Wurz E, Rapp HT, Bayer K, Franke A and Hentschel U (2020) Chloroflexi Dominate the Deep-Sea Golf Ball Sponges Craniella zetlandica and Craniella infrequens Throughout Different Life Stages. Front. Mar. Sci. 7:674. doi: 10.3389/fmars.2020.00674
Received
23 March 2020
Accepted
23 July 2020
Published
21 August 2020
Volume
7 - 2020
Edited by
Joana R. Xavier, University of Porto, Portugal
Reviewed by
Cristiane Cassiolato Pires Hardoim, São Paulo State University, Brazil; Lucia Bongiorni, Institute of Marine Science (CNR), Italy
Updates
Copyright
© 2020 Busch, Wurz, Rapp, Bayer, Franke and Hentschel.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Ute Hentschel, uhentschel@geomar.de
†Deceased March 07, 2020
This article was submitted to Deep-Sea Environments and Ecology, a section of the journal Frontiers in Marine Science
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.