- 1Changing Oceans Group, School of Geosciences, The University of Edinburgh, Edinburgh, United Kingdom
- 2Singapore Centre for Environmental Life Sciences Engineering, Nanyang Technological University, Singapore, Singapore
- 3School of Energy, Geoscience, Infrastructure and Society, Heriot-Watt University, Riccarton, United Kingdom
Holobionts formed by a host organism and associated symbionts are key biological units in marine ecosystems where they are responsible for fundamental ecosystem services. Therefore, understanding anthropogenic impacts on holobionts is essential. Sponges (Phylum Porifera) are ideal holobiont models. They host a complex microbial community and provide ecosystem services including nutrient cycling. At bathyal depths, sponges can accumulate forming dense sponge ground habitats supporting biodiverse associated communities. However, the impacts of spilled oil and dispersants on sponge grounds cannot be understood without considering exposures mediated through sponge filtration of marine snow particles. To examine this, we exposed the model sponge Halichondria panicea to oil, dispersant and “marine oil snow” contaminated seawater and elucidate the complex molecular response of the holobiont through metatranscriptomics. While the host response included detoxification and immune response pathways, the bacterial symbiotic response differed and was at least partially the result of a change in the host environment rather than a direct response to hydrocarbon exposure. As the sponge host reduced its pumping activity and internal tissue oxygen levels declined, the symbionts changed their metabolism from aerobic to anaerobic pathways possibly via quorum sensing. Furthermore, we found evidence of hydrocarbon degradation by sponge symbionts, but sponge mortality (even when exposed to low concentrations of hydrocarbons) implied this may not provide the holobiont with sufficient resilience against contaminants. Given the continued proposed expansion of hydrocarbon production into deep continental shelf and slope settings where sponge grounds form significant habitats it is important that dispersant use is minimised and that environmental impact assessments carefully consider the vulnerability of sponge holobionts.
Introduction
Symbiosis as defined by De Bary (1879) refers to close, long-term and stable associations between organisms of different species. In marine ecosystems, symbiosis between microbes and multicellular organisms is ubiquitous and constitutes a key process driving evolution (Björk et al., 2019; Apprill, 2020). From the concept of symbiosis, the term holobiont emerged in the 1990s and is defined as a host and its associated symbionts forming a single biological unit (Dittami et al., 2021). The term holobiont is now widely used in the field of marine ecology (Dittami et al., 2021) and it is understood that holobionts fulfil important ecosystem services in marine environments (Apprill, 2020). Despite their importance, there is minimal understanding of how marine holobionts respond to environmental stress and how stressors affect the ecosystem services they deliver. In tropical corals (the most studied marine holobiont), environmental stress can lead to coral “bleaching,” the dissociation between the coral host and its microalgal symbionts, ultimately degrading tropical coral reef habitats and highlighting the key role of the symbionts in the holobiont health (e.g., Apprill, 2020). Understanding the impact of environmental stress and pollutants on marine holobionts in poorly understood systems is therefore crucial.
Sponges are ideal holobiont models to study the importance of symbiosis in marine organisms (Pita et al., 2016, 2018; Webster and Thomas, 2016). Sponges are key filter-feeding organisms found in many benthic environments (Van Soest et al., 2012) and host complex microbial communities that actively contribute to holobiont metabolism and health (Pita et al., 2018). Sponge associated microbial communities can account for more than 40% of holobiont biomass (Webster and Taylor, 2012) and contribute to core biological functions including carbon and nitrogen cycling (Hoffmann et al., 2009; Maldonado et al., 2012; Li et al., 2014). Mediated through such close metabolic links, sponge holobionts are integral in benthic ecosystem functioning including bentho-pelagic coupling, i.e., the transfer of energy and nutrients from the benthos to the water column (De Goeij et al., 2013). When present at high densities, sponges form important habitats known as sponge grounds. These habitats are found at all depths and regions, including the polar oceans and support high levels of biodiversity, meeting the criteria as both vulnerable marine ecosystems (VMEs) by the UN Food and Agriculture Organisation and as ecologically and biologically significant marine areas (EBSAs) by the Convention on Biological Diversity (Maldonado et al., 2015). But despite their ecological significance and vulnerability, the impacts of pollutants on sponge holobionts remain poorly understood.
Offshore hydrocarbon exploration and production activities are taking place in areas where sponge grounds are present (Williams et al., 2018; Kazanidis et al., 2019; Vad et al., 2019), yet the potential impacts of an oil spill on sponges are not fully understood. During an oil spill, dispersants are often applied to reduce the formation of surface slicks by increasing the dissolution of hydrocarbons into the water column (National Research Council [NRC], 2005). For example, during the Deepwater Horizon oil spill which took place in the Gulf of Mexico in 2010, a record seven million litres of dispersants were applied both at the surface and near the well (Girard and Fisher, 2018). So far, only two studies have investigated the impacts of hydrocarbons and dispersant exposures on adult sponges and sponge larvae (Luter et al., 2019; Vad et al., 2020). Both studies produced water accommodated oil fractions (WAFs – mixtures of oil in seawater) and chemically enhanced WAFs (CEWAFs – mixtures of dispersant and oil in seawater) to expose sponges to hydrocarbons and dispersant. Both studies showed significant changes in sponge gene expression profiles (Luter et al., 2019; Vad et al., 2020). There is also anecdotal evidence that sponges host oil-degrading bacteria (Rubin-Blum et al., 2017). However, no study yet has investigated whether these symbionts provide enhanced resilience to oil contamination. Even more importantly, there are no studies that have examined the effect of oil and dispersant on “marine snow,” a highly significant nutritional source actively filtered and ingested from seawater by sponges.
Marine snow (MS) is composed of organic and inorganic material, as well as prokaryotic and phytoplankton cells (Simon et al., 2002). It has a pivotal role in the transport of organic material from the upper reaches of the water column as it settles to the seafloor (Simon et al., 2002). At the seafloor, marine snow forms a critical food source for many benthic filter and suspension-feeders (Newell et al., 2005). Marine oil snow (MOS), obtained when small oil droplets get enclosed in marine snow aggregates (Duran Suja et al., 2017), was first observed during the Deepwater Horizon oil spill (Passow et al., 2012) during which MOS was detected within 2 weeks of the spill (Fu et al., 2014). While MOS has also been successfully produced in laboratory experiments by mixing crude oil, dispersant and seawater (Fu et al., 2014; Duran Suja et al., 2017), the impact of MOS on benthic organisms is not well characterised. One laboratory study found that MOS significantly reduced the survival of macro-invertebrates through exposure to oil and hypoxic conditions (Van Eenennaam et al., 2018). Cold-water corals were significantly affected by MOS in the field following the Deepwater Horizon oil spill, displaying a strong stress response at molecular level and visual signs of injury with limited recovery even 7 years after the spill (DeLeo et al., 2018; Girard and Fisher, 2018). To date no study has examined the impact of MOS on the sponge holobiont.
We addressed this knowledge gap by conducting an experiment during which the widely distributed intertidal sponge model Halichondria panicea (commonly known as the breadcrumb sponge) was exposed to MS, MOS, and CEWAF. Metatranscriptomic sequencing was undertaken after 3 days of exposure to determine the holobiont response to hydrocarbon contamination. H. panicea symbiotic community composition is dominated by alphaproteobacterium “Candidatus Halichondribacter symbioticus” (Althoff et al., 1998; Wichels et al., 2006; Knobloch et al., 2019), of which the reference genome has been recently investigated (Knobloch et al., 2020). To inform the metatranscriptomic results, treatment properties (pH, dissolved oxygen concentration and hydrocarbon concentrations) were monitored throughout the exposure. In addition, sponge pumping activity was assessed by measuring tissue oxygen concentration depth profiles. We investigated three hypotheses: (1) the sponge holobiont can be exposed to hydrocarbons through consumption of MOS; (2) the response to hydrocarbons is shared between the sponge host and its associated microbial symbionts; and (3) the sponge associated microbial community is capable of hydrocarbon degradation.
Materials and Methods
Sampling and Experimental Design
Sponge and seawater samples were collected from the Berwickshire coast of the North Sea as detailed in Supplementary Information (SI). Sponge samples were identified by observation of their spicules under light microscopy and the organisation of their skeleton following nitric acid dissolution of the sponge tissue. A flow-through experimental apparatus was built in which each sponge sample was kept in an individual incubation chamber (as described in Vad et al., 2020 and in the Supplementary Information). Sponges were left in the incubation chambers to acclimatise for 48 h in seawater before the beginning of the experiment. The experimental apparatus was kept in the dark to avoid photo oxidation and in a temperature-controlled room (10°C). At the start of the experiment, header flasks connected to the incubation chambers, were filled with relevant treatment solutions (see next section). The experiment lasted for 7 days, exposure was continuous, and the achieved hydrocarbon concentrations in each treatment were measured on day 1, 3, and 5 of the exposure (see section entitled “Physico-Chemical Characterisation of Treatment Solutions”).
CEWAF, MS and MOS Preparation
CEWAF solutions were produced following the Chemical Response to Oil Spills: Ecological Research Forum methodology (Singer et al., 2000) with Schiehallion crude oil (BP) and dispersant Slickgone NS (Dasic International). Further details regarding Schiehallion oil and dispersant Slickgone NS are given in the Supplementary Information. In this study, a CEWAF solution (30 L in total) at a nominal oil loading of 1 g of crude oil per L was prepared and Slickgone was applied at a volume ratio of 1:10 as advised by the manufacturers. The mixtures were then mixed with a magnetic stirrer at a speed of 300 rpm (small vortex visible) for 18 h at 10°C in the dark. The mixtures were allowed to stand for 3 h and the aqueous phases (avoiding non-dispersed/solubilised oil or dispersant) were sub-sampled into clean (autoclaved and acid-washed with 5% nitric acid) screw-capped glass bottles with Teflon caps. The CEWAF solutions were stored at 4°C until used in the experiment. MS aggregates were obtained by rotating three bottles of seawater (to allow for enough aggregates to form; as established in preliminary trials) at 10°C throughout the length of the experiment. MOS aggregates were similarly obtained by rotation three bottles of CEWAF solutions at 10°C throughout the duration of the experiment. During preliminary trials, it was established that the MS and MOS aggregates took 2 to 3 days to synthesise, so the preparation of these treatments were started at the same time as the acclimatisation period. Each day during the exposure, three to five MS and MOS aggregates of 2–3 mm across in diameter were gently pipetted out of the bottles and added to the MS and MOS treated incubation chambers.
Physico-Chemical Characterisation of Treatment Solutions
Seawater pH and dissolved oxygen concentrations were measured in each incubation chamber throughout the experiment. Changes in dissolved oxygen concentrations and pH were detected by performing an ANOVA at each time point after checking for normality and homoscedasticity. Tukey’s all-pair comparisons tests were then used to identify statistical differences (Hothorn et al., 2008; R Core Team, 2021).
Seawater samples (90 mL) were also collected for liquid-liquid hydrocarbon extraction using analytical grade dichloromethane. Following extraction, samples were sent to Terra Tek for speciated USEPA PAHs GC-MS analysis targetting naphthalene, acenaphtylene, acenaphthene, fluorene, phenanthrene, anthracene, fluoranthene, pyrene, benzo[a]anthracene, chrysene, benzo[b]fluoranthene, benzo[k]fluoranthene, benzo[a]pyrene, indeno[1,2,3-cd]pyrene, dibenzo[ah]anthracene and benzo[ghi]perylene. One millilitre of extracts was spiked with internal standard in a 2 mL vial, capped and analysed. The GC-MS analysis was performed on an Agilent 6890 GC with 5973 inert MSD detector in SIM mode, on a 30 m × 0.25 mm ID × 0.25 μm df column with 5 m guard column (Restek p/n 13623-124). PERMANOVA analysis (Oksanen et al., 2020) were applied to test for statistical differences in hydrocarbon concentrations between treatments and time points.
Sponge Tissue Dissolved Oxygen Concentration Profile
Dissolved oxygen concentration was measured within sponge tissues as a proxy for pumping activity (Hoffmann et al., 2008) using a Presens needle spot mounted on a Presens micromanipulator. Measurements were taken at the sponge’s surface and then every 100 μm within its tissue. Three profiles were measured for each sponge at every time point and each profile location was chosen randomly on each sponge. At every time point, statistical differences in dissolved oxygen concentration through the sponge tissue depth between treatments were determined by constructing repeated-measures linear mixed effect models where oxygen concentration was the response variable and depth and treatment were the explanatory variables (Pinheiro et al., 2021).
RNA Extraction and Sequencing
Total RNA was extracted from sponge tissue samples collected after 3 days of exposure and preserved in RNA later at −20°C. Extractions were performed using Qiagen Total RNA Blood and Tissue extraction kits following manufacturer’s instructions. DNA was removed from the extractions using Qiagen DNase kits and RNA extracts were eluded into 30 μL of DNA/RNA free sterile water. RNA quality and quantity were then, respectively, assessed by spectrophotometer using a NanoDrop and by fluorometer using a Qiagen Qubit kit. Only RNA samples with 260/230 and 260/280 ratios between 1.8 and 2.2 were submitted for sequencing. One sample (MS3) failed this quality check and was therefore not sent to sequencing. RNA extractions were sent to Edinburgh Genomics for library preparation and NovaSeq 100PE sequencing. Library preparation was performed by completing a RiboZero treatment with prokaryotic and eukaryotic rRNA depletions. Paired-end reads were run to 100 bp to yield at least 750 M + 750 M reads.
Metatranscriptomic Analysis
When sequencing data were received, read quality was assessed using FastQC (v0.11.9; Andrews, 2015). Short (<75 bp) and low-quality sequences (Phred score < 30) were removed using Cutadapt (v3.1; Martin, 2011). Reads from all samples were then assembled de novo with Megahit (v1.2.9 with default settings; Li et al., 2015) on the University of Edinburgh Linux Compute Cluster Eddie. Short contigs (<300 bp) were removed from the meta-assembly and the meta-assembly was then annotated with Diamond Blast (v2.0.9 with default settings; Buchfink et al., 2021) against the nr database (accessed in November 2020). An E value cut-off of 1e-5 was used to select the best blast hits. Bowtie2 (v2.4.2) was used to align the reads back to the meta-assembly with setting “sensitive local” (Langmead and Salzberg, 2012). Raw and normalised transcript abundance was calculated using RSEM with default settings (v1.3.2; Li and Dewey, 2011). Taxonomic identification of the meta-assembly contigs was assessed using MEGAN (v6.21.2; Huson et al., 2016) to separate sponge (host) and bacterial (symbionts) assemblies. In addition, presence of key alphaproteobacterium symbiont “Candidatus Halichondribacter symbioticus” was assessed by mapping all reads to the bacterial reference genome (Knobloch et al., 2019) with Bowtie2 with setting “sensitive local” (v2.4.2). Differential expression analysis was then performed on the sponge and symbiotic dataset separately in R using DESeq2 (v3.13; Love et al., 2014). Differentially expressed genes compared to the control group were identified at a FDR adjusted value of 0.01. KEGG annotations were obtained with GhostKoala (Kanehisa et al., 2016; accessed in January 2021). To support the differential expression analysis, bacterial metagenomes were compared between all treatments. Prokka (v1.14.5 with default metagenome settings) was first used to convert bacterial transcripts obtained in each treatment to protein sequences (Seemann, 2014). MicrobeAnnotator (v2.0.5 with diamond method; Ruiz-Perez et al., 2021) was then used to annotate the proteins through searches in the KoFamscan, UniProt Swissprot, RefSeq and Trembl databases (last accessed in May 2022). This enhanced annotation enables an in-depth comparisons of metagenome completeness. Further details on how figures were generated for this manuscript is provided in the Supplementary Information.
Results
Gross Observations
Sponges in the control and MS treatment retained their bright yellow colour and appeared healthy throughout the experiment (7-day exposure). In the MOS, CEWAF and CEWAF + MOS treatments, sponge samples slowly turned grey and dark and eventually died (Figure 1A). Samples in the CEWAF + MOS and CEWAF treatment all died after 5 days of exposure while samples exposed to MOS alone died after 7 days.
Figure 1. (A) Change in tissue colouration in the CEWAFMOS treatment compared to control over time. (B) Tissue oxygen concentration (μmol/L) across sponge tissue depth (μm) in each treatment conditions after 3 days of exposure. Note that similar profile was obtained with measurements taken after 5 days of exposure (no statistically significant differences were found between time points within each treatment). Different lowercase letters indicate a significant difference between treatments.
Physicochemistry of Treatments
Overall, seawater dissolved oxygen concentration and pH remained constant in the control, MS and MOS treatment throughout the experiment (Supplementary Figure 1). Dissolved oxygen concentration remained above 320 μmol/L and pH varied between 7.7 and 8.1 in the control, MS and MOS treatments. In contrast, dissolved oxygen concentration and pH were significantly lower in the CEWAF and CEWAF + MOS at day 1, 3, and 6 of exposure (Supplementary Figure 1). Dissolved oxygen reached a minimum of 256.2 ± 3.8 μmol/L (final day of exposure) and 259.4 ± 14.9 μmol/L (first day of exposure) in the CEWAF and CEWAF + MOS treatment, respectively, while pH decreased to 7.3 ± 0.1 on the final day of the exposure in both treatments.
None of the 16 polycyclic aromatic hydrocarbons (PAHs) measured in this study were detected in the control and MS treatments throughout the experiment. In the MOS, CEWAF, and CEWAF + MOS treatments, only naphthalene, fluorene, and phenanthrene were measured above limits of detection (Supplementary Figure 2). Concentrations of all three hydrocarbons were highest in the CEWAF + MOS treatment where naphthalene concentrations rose to 8.8 ± 5.3 μg/L after 1 day of exposure while fluorene and phenanthrene concentrations, respectively, reached 0.6 ± 0.3 μg/L and 3.8 ± 0.8 μg/L after 1 day of exposure. As in the CEWAF + MOS treatment, naphthalene was the most common hydrocarbon detected in the MOS and CEWAF treatments where the concentration reached a maximum of 0.6 ± 0.1 μg/L and 7.0 ± 2.4 μg/L after 3 days of exposure, respectively, (Supplementary Figure 2). Accordingly, Σ16PAHs were highest in the CEWAF + MOS treatment where it rose to 13.3 ± 5.5 μg/L after 1 day of exposure while Σ16PAHs peaked at 11.5 ± 2.9 μg/L and 0.7 ± 0.0 μg/L in the CEWAF and MOS treatments after 3 days of exposure, respectively. The PERMANOVA analysis revealed that hydrocarbon concentrations varied significantly between treatment and time points (Supplementary Table 1).
Sponge Tissue Oxygen Concentration Profiles
Oxygen concentration within the sponge tissue dropped in all treatments with tissue depth and were consistent across all time points considered (Figure 1B). However, significant differences were identified in the sponge oxygen depth profile across treatments. In control conditions, the dissolved oxygen concentration in the sponges stabilised at ∼250 μmol/L at ∼500 μm tissue depth. In the presence of MS, the sponge oxygen depth profile reached a plateau at ∼400 μmol/L at ∼500 μm depth. When exposed to hydrocarbons (MOS, CEWAF or CEWAF + MOS treatments), the sponge oxygen concentrations decreased rapidly with tissue depth and reached 0 from ∼600 μm tissue depth (Figure 1B). Statistically significant differences between oxygen depth profiles were found between the control-MS conditions and the MOS-CEWAF-CEWAF + MOS treatments after both 3 and 5 days of exposure (Figure 1B).
Metatranscriptomic Analysis
The de novo metatranscriptomic assembly was composed of 557,483 contigs, built from a total of 789,274,613 paired sequence reads. The assembly GC content was 40.42% and the assembly N50 reached 1,615 base pairs. Overall, 68.9% of the contigs (383,995 contigs) were identified against the nr database. The annotated sequences included 48,221 eukaryotic (12.6%), 325,636 bacterial (84.8%), 1,331 archaeal (0.34%) and 1,209 viral (0.31%) sequences. In total, 14,791 eukaryotic sequences were identified from poriferan sequences (30.7% of eukaryotic sequences) while 10,747 were derived from cnidarian sequences (22.3% of eukaryotic sequences). Bacterial sequences were dominated by sequences identified in Proteobacteria (175,877 sequences equivalent to 54.0% of bacterial sequences) and in the Fibrobacteres-Chlorobi-Bacteroidetes group (57,251 sequences equivalent to 17.6% of bacterial sequences). Presence of key symbiont “Ca. H. symbioticus” was assessed in all samples and overall read alignment scores with the reference genome were consistent across treatments varying between 16.4% ± 0.8 and 20.4% ± 3.6 of all reads (sponge and bacterial; Table 1).
Table 1. Overall alignment scores of all reads against “Ca. H. symbioticus” reference genome in each sample.
In this study, the holobiont expression profiles changed significantly depending on treatment condition. The first component of the PCA analysis, which explained 91% of the variation in expression profiles, separated the CEWAF and CEWAF + MOS exposed sponges from samples exposed to other conditions (Figure 2). The second component of the PCA analysis explained 4% of the variation in expression and separated control samples from sponges exposed to MS and MOS (Figure 2). The largest number of differentially expressed genes (DEGs) were identified in the CEWAF + MOS treatment in the sponge host and bacterial associated community (Figure 3). Overall, the majority of DEGs were found across the CEWAF and CEWAF + MOS treatments.
Figure 2. Principal component plot of metatranscriptome-wide expression profiles in H. panicea samples exposed to control, MS, MOS, CEWAF, and CEWAF + MOS conditions.
Figure 3. Upset plots showing the number of differentially expressed genes in (A) the sponge host transcriptome and (B) the symbiotic bacterial metatranscriptome across treatments relative to control conditions. The abbreviation C.MOS stands for CEWAF + MOS. The nodes below each bar plot illustrates which treatments are compared. The bar plot therefore show the number of differentially expressed genes relative to control conditions common to the treatments being compared.
Several Kyoto Encyclopaedia of Genes and Genomes (KEGG) pathways were significantly enriched amongst sponge and bacterial DEGs. In the sponge host, several core signalling pathways (NF-kappa B, IL-17, NOD-like receptor, TNF and RIG-I-like receptor) were found significantly enriched in samples exposed to all treatments tested in this study (Supplementary Figure 3). Furthermore, the KEGG pathways “fatty acid biosynthesis” and “biotin metabolism” were significantly enriched in sponges exposed to MS and MOS (Supplementary Figure 3). In the CEWAF and CEWAF + MOS treatment, KEGG pathways involved in the MAPK signalling pathway, apoptosis, and tissue organisation (cell adhesion molecules and Hedgehog signalling pathways) as well as immune system (antigen processing and presentation) were found significantly enriched (Supplementary Figure 3). Amongst the bacterial DEGs, KEGG pathways involved in flagellar assembly and NOD-like receptor signalling pathways were found significantly enriched in all four treatments (Supplementary Figures 4, 5). In the MS treatment, many DEGs were involved in nitrogen metabolism while in the MOS treatment, DEGs contributed to the bacterial secretion system (Supplementary Figure 4). The KEGG pathways enrichment networks produced from the bacterial genes in the CEWAF and CEWAF + MOS treatment were almost identical and bacterial DEGs found in these two treatments contributed predominantly to amino acid synthesis, the carbon cycle and quorum sensing (Figure 4 and Supplementary Figure 5). Although KEGG functions associated with hydrocarbon degradation were not significantly enriched, evidence for hydrocarbon degradation was nevertheless identified amongst the bacterial DEGs pool. In total, 134 DEGs linked to hydrocarbon degradation functions were detected and these genes were significantly up-regulated in samples exposed to the CEWAF and CEWAF + MOS treatments (Supplementary Figure 6).
Figure 4. Networks of significantly enriched KEGG pathways amongst the differentially expressed genes in the bacterial symbionts exposed to CEWAF + MOS. Note that KEGG enrichment results for the CEWAF treatment were similar.
To further explore the changes in metabolic activities found through differential expression analysis, metagenome annotation and comparisons were performed through Prokka and MicrobeAnnotator (Supplementary Figure 7). Overall, module completeness significantly increased in all treatments compared to controls, including for core metabolic activities such as amino acid synthesis, aromatics degradation, carbon fixation and nitrogen metabolism amongst others (Supplementary Figure 7).
Discussion
Scale of Holobiont Response Depends on Hydrocarbon Concentrations
Exposure to MOS alone led to a muted but detectable response in the sponge holobiont compared to the significant changes found in gene expression profiles of sponges in the CEWAF and CEWAF + MOS treatments. This is likely due to hydrocarbon concentrations in the MOS treatment being low throughout the exposure, as only three to five marine oil snow aggregates of 2–3 mm diameter were added to the incubation chambers during each day of the exposure. In fact, Σ16PAHs concentrations measured in the MOS treatment (ΣPAH16 0.131 to 0.712 μg/L) were below the range observed (ΣPAH16 0.88 to 6.28 μg/L) in the field during the Deepwater Horizon oil spill (Diercks et al., 2010). Furthermore, the metatranscriptomic profiles of the sponge holobionts exposed to MOS alone were similar to the profiles of sponges exposed to MS. Evidence of sponge feeding was seen in both treatments as KEGG functions such as fatty acid biosynthesis were found significantly enriched amongst the host DEGs pool. Sponge feeding in the MOS and MS treatment is also indicated by the bacterial metatranscriptomic profile. Indeed, the KEGG function flagellar assembly identified as significantly enriched in the bacterial DEGs pool, could, in fact, be indicative of the regeneration of pumping cells (choanocytes) in the sponge tissue. Choanocytes, which are equipped with flagella and constitute the pumping engines of the sponge, are organised in an epithelium known as the choanoderm along a sponge’s internal pumping chambers (De Goeij et al., 2009). These choanocytes are replaced approximately every 6 h to support pumping activities – a rapid cell cycle that matches rates seen in unicellular rather than multicellular organisms (De Goeij et al., 2009). Evidence of feeding in our experiment means that intracellular exposure to hydrocarbons through MOS is therefore highly likely. While only small quantities of MOS aggregates were added to the MOS and CEWAF + MOS treatments, large quantities of MOS were observed in situ following the Deepwater Horizon oil spill (Passow et al., 2012; Halanych et al., 2021). Following the well blowout, benthic species including cold-water coral colonies were observed completely covered in MOS aggregates (DeLeo et al., 2018). Despite the MOS treatment in this study prompting only a relatively small metatranscriptomic response in the sponge holobiont, exposures to higher quantities of MOS (and so higher hydrocarbon concentrations) could lead to stronger changes in gene expression profiles, akin to those seen here in the CEWAF and CEWAF + MOS treatments.
Host VS Symbiont Responses to Hydrocarbon Exposure
This study is the first to apply metatranscriptomics to reveal how sponge symbionts alter their metabolic activities in response to changes in the host tissue as an indirect consequence of environmental pollution. In the CEWAF and CEWAF + MOS treatments, the response to the hydrocarbon exposure between the sponge host and its associated microbial community differed. The KEGG functions found enriched in the host DEG pool highlighted the reliance of the sponge on detoxification, oxidative stress response and immune response pathways. Specifically, the MAPK signalling pathway was found to be a core component of the sponge molecular response to hydrocarbon contamination. This is in line with what has already been described in other marine invertebrates (e.g., Jenny et al., 2016; Park et al., 2020), other sponge species (Châtel et al., 2011) and in H. panicea (Vad et al., 2020). In addition, we detected evidence of sponge host tissue degradation as cell adhesion function and the hedgehog signalling pathway were both significantly increased in the CEWAF and CEWAF + MOS treatments. Indeed, the hedgehog signalling pathway has been identified in the genome of the tropical sponge Amphimedon queenslandica and is known to be involved in tissue development (Adamska et al., 2007). However, in our study, we also found that the bacterial response outweighed the host response by a factor of 8 to 10 in terms of number of DEGs detected. As part of the bacterial response, the DEGs identified were involved in key metabolic activities such as nitrogen and carbon fixation as illustrated by the relatively small number of significantly enriched KEGG functions detected. Furthermore, the list of enriched KEGG functions identified revealed a change from aerobic to anaerobic activities. Metagenome comparisons performed through MicrobeAnnotator confirmed the changes in metabolic activities identified by differential gene expression analysis. It is therefore possible that the bacterial community responded strongly to the severe drop in dissolved oxygen concentration inside sponge tissues, a consequence of the host halting its pumping activity when exposed to hydrocarbons. A change in bacterial community composition is likely to have contributed to the shifts in metabolic activities found in our study; however, stable alignment scores against “Ca. H. symbioticus” reference genome does not indicate a severe shift in microbial community composition. Additional metagenomics analysis would be required to fully address shifts in microbial community composition.
The Role of Quorum Sensing in the Sponge Holobiont
In the CEWAF and CEWAF + MOS treatments, quorum sensing played a key role in symbionts’ response. Quorum sensing is a form of cell-cell communication defined as the regulation of gene expression in response to the release of chemical signal molecules (autoinducers) beyond a minimal stimulatory threshold concentration (e.g., Miller and Bassler, 2001). In this study, it is possible that quorum sensing enabled the symbiotic community response leading to dramatic changes in metabolic processes and a rapid switch from aerobic to anaerobic metabolism. Although quorum sensing functions were significantly enriched within the bacterial DEG pool, quorum sensing could facilitate symbiont-symbiont as well as host-symbiont interactions (Shiner et al., 2005). Quorum sensing compounds have already been isolated from sponge associated bacteria (Zan et al., 2012) or from sponge extracts (Batista et al., 2018) and identified in sponge metabolic profiles (Papale et al., 2020). However, to the best of our knowledge, this is the first study to identify the specific role of quorum sensing in the response of a marine holobiont to environmental stress.
Hydrocarbon Degradation Does Not Lead to Hydrocarbon Resilience in Sponges
Hydrocarbon degradation was identified in this study in sponges exposed to higher concentrations of hydrocarbons (CEWAF and CEWAF + MOS treatment). Oil degrading bacteria are present at low concentration in the marine environment, both in seawater and in sediments (Joye et al., 2016). Hydrocarbon degrading symbionts have also been found in sponges collected from deep-sea hydrocarbon seeps (Rubin-Blum et al., 2017). However, we present the first evidence of hydrocarbon degradation in the intertidal sponge H. panicea. Despite this, the sponges in the hydrocarbon treatments (CEWAF, CEWAF + MOS and MOS) all died within 7 days of exposure. This demonstrates that the degradation of hydrocarbon by symbionts may not provide the holobiont with any significant increase in resilience to hydrocarbons pollution. This confirms observations made during the Deepwater Horizon oil spill, reporting high mortality in sponge assemblages following the well blowout (Beyer et al., 2016; McClain et al., 2019).
Conclusion
Sponge grounds cover up to hundreds of km2 and support high level of biodiversity by increasing the three-dimensionality of the local benthos, altering local hydrodynamics and contributing to nutrient cycling (Maldonado et al., 2015). Nevertheless, hydrocarbon production activities are already taking place, or are being proposed, within sponge grounds and remediation strategies in place at those fields include the use of dispersants (BP, 2014; BP and Genesis Oil and Gas Consultants Ltd, 2018). For example, oil production activities in the NE Atlantic Faroe-Shetland Channel are being carried within the Sponge Belt Nature Conservation Marine Protect Area (MPA) designated to protect sponge grounds (BP, 2014; BP and Genesis Oil and Gas Consultants Ltd, 2018; Kazanidis et al., 2019; Vad et al., 2019). Further developments such as the CAMBO oil field, located close to the MPA, are currently being considered by the United Kingdom government and dispersants are still listed in the development environmental impact assessment (Siccar Point Energy E&P Ltd, 2020). In addition, sponge grounds form a crucial habitat in benthic Artic environments (Maldonado et al., 2015), which also hold significant untapped oil and gas reservoirs (USGS, 2008) and new major shipping routes due to the furthering losses of Arctic sea ice (Bambulyak and Ehlers, 2020). Sampling at sponge grounds is difficult as they are often found in the deep-sea and require the use of remotely operated vehicles. In this experimental study, we used intertidal sponge model H. panicea as a proxy for deeper species. As the molecular pathways involved in the response of H. panicea rely on core molecular function and are well preserved across metazoans (Böhm et al., 2001), we can hypothesise that results on deeper species found at sponge grounds such as Geodia spp would be similar. Here, sponges exposed to very low concentrations of hydrocarbons (MOS only treatment) died after 7 days of exposure. Our study therefore demonstrates that by contributing to the production of MOS and increasing the concentrations of hydrocarbons in seawater, the use of dispersants increased the risk posed by hydrocarbon contamination to sponges and should therefore be limited within areas that contain sponge grounds. Furthermore, our study also highlights the importance of considering holobiont in ecotoxicology investigations. Future experiments should therefore aim at understanding how changes in the functioning of key holobionts can lead to shifts in the vulnerable habitats they provide, such as sponge grounds.
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.ebi.ac.uk/ena, ERP134402.
Author Contributions
JV and LD ran the experiment, analysed the data, and wrote the first manuscript draft. SS contributed to the metatranscriptomics analysis. TH and JR secured funding for the study. All authors contributed to the conceptualisation of the study and the methodology, reviewed, and edited the manuscript.
Funding
This study received support from the Marine Alliance for Science and Technology for Scotland (grant number OGSG14) and from the NERC Natural Environment Research Council Centre for Doctoral Training in Oil & Gas. This study was part of the European Union’s Horizon 2020 Research and Innovation Programme under grant agreement No. 678760 (ATLAS) and No. 818123 (iAtlantic). This paper reflects the authors’ view alone and the European Union cannot be held responsible for any use that may be made of the information contained.
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.
Publisher’s Note
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.
Acknowledgments
We would like to thank S. Hennige for comments on the manuscript prior to submission and to our reviewers whose feedback improved this study. We made use of the resources provided by the Edinburgh Compute and Data Facility (ECDF; http://www.ecdf.ed.ac.uk/).
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2022.909853/full#supplementary-material
References
Adamska, M., Matus, D. Q., Adamski, M., Green, K., Rokhsar, D. S., Martindale, M. Q., et al. (2007). The evolutionary origin of hedgehog proteins. Curr. Biol. 2007, R836–R837. doi: 10.1016/j.cub.2007.08.010
Althoff, K., Schütt, C., Steffen, R., Batel, R., and Müller, W. E. G. (1998). ‘Evidence for a symbiosis between bacteria of the genus Rhodobacter and the marine sponge Halichondria panicea: Harbor also for putatively toxic bacteria?’ Mar. Biol. 130, 529–536. doi: 10.1007/s002270050273
Andrews, S. (2015). FastQC: a quality control tool for high throughput sequence data. Cambridge, MA: Babraham Institute.
Apprill, A. (2020). ‘The Role of Symbioses in the Adaptation and Stress Responses of Marine Organisms’. Annu. Rev. Mar. Sci. 12, 291–314. doi: 10.1146/annurev-marine-010419-010641
Bambulyak, A., and Ehlers, S. (2020). Oil spill damage: a collision scenario and financial liability estimations for the Northern Sea Route area. Ship Technol. Res. 148–164. doi: 10.1080/09377255.2020.1786932
Batista, D., Costa, R., Carvalho, A. P., Batista, W. R., Rua, C. P. J., de Oliveira, L., et al. (2018). ‘Environmental conditions affect activity and associated microorganisms of marine sponges’. Mar. Env. Res. 142, 59–68. doi: 10.1016/j.marenvres.2018.09.020
Beyer, J., Trannum, H. C., Bakke, T., Hodson, P. V., and Collier, T. K. (2016). ‘Environmental effects of the Deepwater Horizon oil spill: a review’. Mar. Pollut. Bull. 110, 28–51. doi: 10.1016/j.marpolbul.2016.06.027
Björk, J. R., Díez-Vives, C., Astudillo-García, C., Archie, E. A., and Montoya, J. M. (2019). Vertical transmission of sponge microbiota is inconsistent and unfaithful. Nat. Ecol. Evol. 2019:935. doi: 10.1038/s41559-019-0935-x
Böhm, M., Hentschel, U., Friedrich, A., Fieseler, L., Steffen, R., Gamulin, V., et al. (2001). ‘Molecular response of the sponge Suberites domuncula to bacterial infection’. Mar. Biol. 139, 1037–1045. doi: 10.1007/s002270100656
BP (2014). Schiehallion Offshore Oil Pollution Emergency Plan and Justification Document. London: BP.
BP and Genesis Oil and Gas Consultants Ltd (2018). Alligin Field Development Environmental Statement. London: BP.
Buchfink, B., Reuter, K., and Drost, H. G. (2021). ‘Sensitive protein alignments at tree-of-life scale using DIAMOND’. Nat. Methods 18, 366–368. doi: 10.1038/s41592-021-01101-x
Châtel, A., Talarmin, H., Hamer, B., Schröder, H. C., Müller, W. E. G., and Dorange, G. (2011). ‘MAP kinase cell signaling pathway as biomarker of environmental pollution in the sponge Suberites domuncula’. Ecotoxicology 20, 1727–1740. doi: 10.1007/s10646-011-0706-1
De Goeij, J. M., De Kluijver, A., Van Duyl, F. C., Vacelet, J., Wijffels, R. H., De Goeij, F. P. M., et al. (2009). Cell kinetics of the marine sponge Halisarca caerulea reveal rapid cell turnover and shedding. J. Exp. Biol. 212, 3892–3900. doi: 10.1242/jeb.034561
De Goeij, J. M., van Oevelen, D., Vermeij, M. J. A., Osinga, R., Middelburg, J. J., de Goeij, A. F. P. M., et al. (2013). ‘Surviving in a Marine Desert: the Sponge Loop Retains Resources Within Coral Reefs’. Science 342, 108–110. doi: 10.1126/science.1241981
DeLeo, D. M., Herrera, S., Lengyel, S. D., Quattrini, A. M., Kulathinal, R. J., and Cordes, E. E. (2018). ‘Gene expression profiling reveals deep-sea coral response to the Deepwater Horizon oil spill’. Mole. Ecol. 27, 4066–4077. doi: 10.1111/mec.14847
Diercks, A.-R., Highsmith, R. C., Asper, V. L., Joung, D., Zhou, Z., Guo, L., et al. (2010). ‘Characterization of subsurface polycyclic aromatic hydrocarbons at the Deepwater Horizon site’. Geophys. Res. Lett. 37:L20602. doi: 10.1029/2010GL045046
Dittami, S. M., Leblanc, C., Arboleda, E., Auguet, J.-C., Bigalke, A., Briand, E., et al. (2021). A community perspective on the concept of marine holobionts: current status, challenges, and future directions. PeerJ 17:10911. doi: 10.7717/peerj.10911
Duran Suja, L., Summers, S., and Gutierrez, T. (2017). ‘Role of EPS, Dispersant and Nutrients on the Microbial Response and MOS Formation in the Subarctic Northeast Atlantic. Front. Microb. 8, 1–15. doi: 10.3389/fmicb.2017.00676
Fu, J., Gong, Y., Zhao, X., O’Reilly, S. E., and Zhao, D. (2014). ‘Effects of oil and dispersant on formation of marine oil snow and transport of oil hydrocarbons’. Env. Sci. Tech. 48, 14392–14399. doi: 10.1021/es5042157
Girard, F., and Fisher, C. R. (2018). ‘Long-term impact of the Deepwater Horizon oil spill on deep-sea corals detected after seven years of monitoring’. Biol. Conserv. 225, 117–127. doi: 10.1016/J.BIOCON.2018.06.028
Halanych, K., Ainsworth, C., Cordes, E., Dodge, R., Huettel, M., Mendelssohn, I., et al. (2021). ‘Effects of Petroleum By-Products and Dispersants on Ecosystems’. Oceanography 34, 152–163. doi: 10.5670/oceanog.2021.123
Hoffmann, F., Radax, R., Woebken, D., Holtappels, M., Lavik, G., Rapp, H. T., et al. (2009). ‘Complex nitrogen cycling in the sponge Geodia barretti’. Env. Microb. 11, 2228–2243. doi: 10.1111/j.1462-2920.2009.01944.x
Hoffmann, F., Røy, H., Bayer, K., Hentschel, U., Pfannkuchen, M., Brümmer, F., et al. (2008). ‘Oxygen dynamics and transport in the Mediterranean sponge Aplysina aerophoba’. Mar. Biol. 153, 1257–1264. doi: 10.1007/s00227-008-0905-3
Hothorn, T., Bretz, F., and Westfall, P. (2008). ‘Simultaneous Inference in General Parametric Models’. Biometr. J. 50, 346–363. doi: 10.1002/bimj.200810425
Huson, D. H., Beier, S., Flade, I., Górska, A., El-Hadidi, M., Mitra, S., et al. (2016). ‘MEGAN Community Edition - Interactive Exploration and Analysis of Large-Scale Microbiome Sequencing Data’. PLoS Comput. Biol. 12:e1004957. doi: 10.1371/journal.pcbi.1004957
Jenny, M. J., Walton, W. C., Payton, S. L., Powers, J. M., Findlay, R. H., Shields, B. O., et al. (2016). ‘Transcriptomic evaluation of the American oyster, Crassostrea virginica, deployed during the Deepwater Horizon oil spill: evidence of an active hydrocarbon response pathway. Mar. Env. Res. 120, 166–181. doi: 10.7266/N7CR5R80
Joye, S. B., Kleindienst, S., Gilbert, J. A., Handley, K. M., Weisenhorn, P., Overholt, W. A., et al. (2016). ‘Responses of Microbial Communities to Hydrocarbon Exposures Oceanography’. Oceanography 29, 136–149. doi: 10.5670/oceanog.2016.78
Kanehisa, M., Sato, Y., and Morishima, K. (2016). ‘BlastKOALA and GhostKOALA: KEGG tools for functional characterization of genome and metagenome sequences’. J. Mole. Biol. 428, 726–731. doi: 10.1016/j.jmb.2015.11.006
Kazanidis, G., Vad, J., Henry, L.-A., Neat, F., Berx, B., Georgoulas, K., et al. (2019). Distribution of Deep-Sea Sponge Aggregations in an Area of Multisectoral Activities and Changing Oceanic Conditions. Front. Mar. Sci. 6:163. doi: 10.3389/fmars.2019.00163
Knobloch, S., Jóhannsson, R., and Marteinsson, V. (2019). ‘Bacterial diversity in the marine sponge Halichondria panicea from Icelandic waters and host-specificity of its dominant symbiont “Candidatus Halichondribacter symbioticus. FEMS Microb. Ecol. 95:220. doi: 10.1093/FEMSEC/FIY220
Knobloch, S., Jóhannsson, R., and Marteinsson, V. P. (2020). Genome analysis of sponge symbiont ‘Candidatus Halichondribacter symbioticus’ shows genomic adaptation to a host-dependent lifestyle. Environ. Microb. 22, 483–498. doi: 10.1111/1462-2920.14869
Langmead, B., and Salzberg, S. (2012). ‘Fast gapped-read alignment with Bowtie2’. Nat. Methods 9, 357–359.
Li, B., and Dewey, C. N. (2011). ‘RSEM: accurate transcript quantification from RNA-Seq data with or without a reference genome’. BMC Bioinform. 12:323. doi: 10.1186/1471-2105-12-323
Li, D., Liu, C.-M., Luo, R., Sadakane, K., and Lam, T.-W. (2015). MEGAHIT: an ultra-fast single-node solution for large and complex metagenomics assembly via succinct de Bruijn graph. Bioinformatics 2015:33. doi: 10.1093/bioinformatics/btv033
Li, Z.-Y., Wang, Y.-Z., He, L.-M., and Zheng, H.-J. (2014). ‘Metabolic profiles of prokaryotic and eukaryotic communities in deep-sea sponge Lamellomorpha sp. indicated by metagenomics’. Sci. Rep. 4, 1–11. doi: 10.1038/srep03895
Love, M. I., Huber, W., and Anders, S. (2014). ‘Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2’. Gen. Biol. 15:550. doi: 10.1186/s13059-014-0550-8
Luter, H. M., Whalan, S., Andreakis, N., Abdul Wahab, M., Botté, E. S., Negri, A. P., et al. (2019). ‘The Effects of Crude Oil and Dispersant on the Larval Sponge Holobiont. mSystems 4, 1–17. doi: 10.1128/msystems.00743-19
Maldonado, M., Aguilar, R., Bannister, R. J., Bell, J. J., Conway, K. W., Dayton, P. K., et al. (2015). “Sponge Grounds as Key Marine Habitats: A Synthetic Review of Types, Structure, Functional Roles, and Conservation Concerns,” in Marine Animal Forests, eds S. Rossi, L. Bramanti, A. Gori, and C. Orejas (Cham: Springer International Publishing), 1–39. doi: 10.1007/978-3-319-17001-5_24-1
Maldonado, M., Ribes, M., and van Duyl, F. C. (2012). Nutrient Fluxes through Sponges. Biology, Budgets, and Ecological Implications. Adv. Mar. Biol. 2012:5. doi: 10.1016/B978-0-12-394283-8.00003-5
Martin, M. (2011). ‘Cutadapt removes adapter sequences from high-throughput sequencing reads’. EMBnet. J. 17:10. doi: 10.14806/ej.17.1.200
McClain, C. R., Nunnally, C., and Benfield, M. C. (2019). ‘Persistent and substantial impacts of the Deepwater Horizon oil spill on deep-sea megafauna’. R. Soc. Open Sci. 6:191164. doi: 10.1098/rsos.191164
Miller, M. B., and Bassler, B. L. (2001). Quorum sensing in bacteria. Annu. Rev. Microbiol. 2001, 165–199. doi: 10.1146/annurev.micro.55.1.165.National
National Research Council [NRC] (2005). Oil Spill Dispersants: Efficacy and Effects. Washington DC: National Academy Press.
Newell, C. R., Pilskaln, C. H., Robinson, S. M., and MacDonald, B. A. (2005). ‘The contribution of marine snow to the particle food supply of the benthic suspension feeder, Mytilus edulis. J. Exp. Mar. Biol. Ecol. 321, 109–124. doi: 10.1016/j.jembe.2005.01.006
Oksanen, J., Blanchet, G., Friendly, M., Kindt, R., Legendre, P., McGlinn, D., et al. (2020). vegan: Community Ecology Package’. R package version 2.5-7. Available online at: https://CRAN.R-project.org/package=vegan (accessed date 2022-04-17).
Papale, M., Rizzo, C., Fani, R., Bertolino, M., Costa, G., Paytuví-Gallart, A., et al. (2020). Exploring the Diversity and Metabolic Profiles of Bacterial Communities Associated With Antarctic Sponges (Terra Nova Bay, Ross Sea). Front. Ecol. Evol. 8:268. doi: 10.3389/fevo.2020.00268
Park, J. C., Hagiwara, A., Park, H. G., and Lee, J. S. (2020). ‘The glutathione S-transferase genes in marine rotifers and copepods: identification of GSTs and applications for ecotoxicological studies’. Mar. Pollut. Bull. 156:111080. doi: 10.1016/j.marpolbul.2020.111080
Passow, U., Ziervogel, K., Asper, V., and Diercks, A. (2012). Marine snow formation in the aftermath of the Deepwater Horizon oil spill in the Gulf of Mexico. Env. Res. Lett. 7:035301. doi: 10.1088/1748-9326/7/3/035301
Pinheiro, J., Bates, D., Debroy, S., and Sarkar, D., and R Core Team. (2021). nlme: Linear and Nonlinear Mixed Effects Models_’. R package version 3.1-152. Available online at: https://CRAN.R-project.org/package=nlme (accessed date: 2022-03-25).
Pita, L., Fraune, S., and Hentschel, U. (2016). Emerging sponge models of animal-microbe symbioses. Front. Microbiol. 2016:2102. doi: 10.3389/fmicb.2016.02102
Pita, L., Rix, L., Slaby, B. M., Franke, A., and Hentschel, U. (2018). ‘The sponge holobiont in a changing ocean: from microbes to ecosystems’. Microbiome 6:46. doi: 10.1186/s40168-018-0428-1
R Core Team (2021). R: A language and environment for statistical computing. Vienna: R Foundation for statistical Computing.
Rubin-Blum, M., Antony, C. P., Borowski, C., Sayavedra, L., Pape, T., Sahling, H., et al. (2017). Short-chain alkanes fuel mussel and sponge Cycloclasticus symbionts from deep-sea gas and oil seeps. Nat. Microbiol. 2:93. doi: 10.1038/nmicrobiol.2017.93
Ruiz-Perez, C. A., Conrad, R. E., and Konstantinidis, K. T. (2021). ‘MicrobeAnnotator: a user-friendly, comprehensive functional annotation pipeline for microbial genomes’. BMC Bioinform. 22, 1–16. doi: 10.1186/S12859-020-03940-5/FIGURES/4
Seemann, T. (2014). ‘Prokka: rapid prokaryotic genome annotation’. Bioinformatics 30, 2068–2069. doi: 10.1093/BIOINFORMATICS/BTU153
Shiner, E. K., Rumbaugh, K. P., and Williams, S. C. (2005). Interkingdom signalling: deciphering the language of acyl homoserine lactones. FEMS Microbiol. Rev. 2005, 935–947. doi: 10.1016/j.femsre.2005.03.001
Siccar Point Energy E&P Ltd (2020). Cambo Oil Field, UKCS Blocks 204/4a, 204/5a, 204/9a and 204/10a Environmental Impact Assessment (EIA). ABERDEEN: Siccar Point Energy E&P Ltd
Simon, M., Grossart, H. P., Schweitzer, B., and Ploug, H. (2002). ‘Microbial ecology of organic aggregates in aquatic ecosystems’. Aquatic Microb. Ecol. 28, 175–211. doi: 10.3354/ame028175
Singer, M. M., Aurand, D., Bragin, G. E., Clark, J. R., Coelho, G. M., Sowby, M. L., et al. (2000). ‘Standardization of the preparation and quantitation of water-accommodated fractions of petroleum for toxicity testing’. Mar. Pollut. Bull. 40, 1007–1016. doi: 10.1016/S0025-326X(00)00045-X
USGS (2008). Circum-Arctic Resource Appraisal: Estimates of Undiscovered Oil and Gas North of the Arctic Circle. USGS Fact Sheet. Reston: U.S. Department of the Interior, U.S. Geological Survey, 2008–3049.
Vad, J., Dunnett, F., Liu, F., Montagner, C. C., Roberts, J. M., and Henry, T. B. (2020). Soaking up the oil: Biological impacts of dispersants and crude oil on the sponge Halichondria panicea. Chemosphere 257:127109. doi: 10.1016/j.chemosphere.2020.127109
Vad, J., Kazanidis, G., Henry, L.-A., Jones, D. O. B., Gates, A. R., and Roberts, J. M. (2019). ‘Environmental controls and anthropogenic impacts on deep-sea sponge grounds in the Faroe-Shetland Channel, NE Atlantic: the importance of considering spatial scale to distinguish drivers of change. ICES J. Mar. Sci. 77, 451–461. doi: 10.1093/icesjms/fsz185
Van Eenennaam, J. S., Rahsepar, S., Radović, J. R., Oldenburg, T. B. P., Wonink, J., Langenhoff, A. A. M., et al. (2018). ‘Marine snow increases the adverse effects of oil on benthic invertebrates’. Mar. Pollut. Bull. 126, 339–348. doi: 10.1016/j.marpolbul.2017.11.028
Van Soest, R. W. M., Boury-Esnault, N., Vacelet, J., Dohrmann, M., Erpenbeck, D., De Voogd, N. J., et al. (2012). ‘Global Diversity of Sponges (Porifera)’. PLoS One 7:e35105. doi: 10.1371/journal.pone.0035105
Webster, N. S., and Taylor, M. W. (2012). ‘Marine sponges and their microbial symbionts: love and other relationships’. Env. Microbiol. 14, 335–346. doi: 10.1111/j.1462-2920.2011.02460.x
Webster, N. S., and Thomas, T. (2016). ‘Defining the Sponge Hologenome. mBio 7, 1–14. doi: 10.1128/mBio.00135-16.Invited
Wichels, A., Würtz, S., Döpke, H., Schütt, C., and Gerdts, G. (2006). ‘Bacterial diversity in the breadcrumb sponge Halichondria panicea (Pallas)’. FEMS Microb. Ecol. 56, 102–118. doi: 10.1111/j.1574-6941.2006.00067.x
Williams, A., Althaus, F., MacIntosh, H., Loo, M., Gowlett-Holmes, K., Tanner, J. E., et al. (2018). ‘Characterising the invertebrate megafaunal assemblages of a deep-sea (200–3000 m) frontier region for oil and gas exploration: the Great Australian Bight, Australia. Deep Sea Res. Part II 15, 78–91. doi: 10.1016/j.dsr2.2018.07.015
Zan, J., Cicirelli, E. M., Mohamed, N. M., Sibhatu, H., Kroll, S., Choi, O., et al. (2012). ‘A complex LuxR-LuxI type quorum sensing network in a roseobacterial marine sponge symbiont activates flagellar motility and inhibits biofilm formation’. Mole. Microbiol. 85, 916–933. doi: 10.1111/j.1365-2958.2012.08149.x
Keywords: holobiont, marine snow, marine oil snow, metatranscriptomics, marine sponges, Halichondria panicea
Citation: Vad J, Duran Suja L, Summers S, Henry TB and Roberts JM (2022) Marine Sponges in a Snowstorm – Extreme Sensitivity of a Sponge Holobiont to Marine Oil Snow and Chemically Dispersed Oil Pollution. Front. Microbiol. 13:909853. doi: 10.3389/fmicb.2022.909853
Received: 31 March 2022; Accepted: 06 June 2022;
Published: 15 July 2022.
Edited by:
Cristiane Cassiolato Pires Hardoim, São Paulo State University, BrazilReviewed by:
Roger C. Prince, Stonybrook Apiary, United StatesAnoop Alex, University of Porto, Portugal
Stephen Knobloch, Matis Ltd., Iceland
Copyright © 2022 Vad, Duran Suja, Summers, Henry and Roberts. 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: Johanne Vad, Johanne.vad@ed.ac.uk
†ORCID: Johanne Vad, orcid.org/0000-0003-2215-0491; Stephen Summers, orcid.org/0000-0001-7680-701X; Theodore B. Henry, orcid.org/0000-0002-9675-9454; J. Murray Roberts, orcid.org/0000-0003-1688-5133