# How Environmental Factors Affect the Production of Guanidine Alkaloids by the Mediterranean Sponge Crambe crambe

Eva Ternon, Erica D Perino, Renata D Manconi, Roberto D Pronzato, Olivier P Thomas

<figure><img src="image_1.png" /></figure>

## To cite this version:

Eva Ternon, Erica D Perino, Renata D Manconi, Roberto D Pronzato, Olivier P Thomas. How Environmental Factors Affect the Production of Guanidine Alkaloids by the Mediterranean Sponge Crambe crambe. Marine drugs, 2017, 15 (6), 10.3390/md15060181. [DOI: 10.3390/md15060181](https://doi.org/10.3390/md15060181)

HAL Id: hal-01854419

https://hal.science/hal-01854419

Submitted on 18 Dec 2020

HAL is a multi-disciplinary open access archive for the deposit and dissemination of scientific research documents, whether they are published or not. The documents may come from teaching and research institutions in France or abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est destinée au dépôt et à la diffusion de documents scientifiques de niveau recherche, publiés ou non, émanant des établissements d'enseignement et de recherche français ou étrangers, des laboratoires publics ou privés.

<figure><img src="image_2.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_3.png" /></figure>

Article

# How Environmental Factors Affect the Production of Guanidine Alkaloids by the Mediterranean Sponge *Crambe crambe*

Eva Ternon <fnref n="1" />,*, Erica Perino <fnref n="2" />, Renata Manconi <fnref n="3" />, Roberto Pronzato <fnref n="2" /> and Olivier P. Thomas <fnref n="1" />,<fnref n="4" />,*

Academic Editor: Paul Long

Received: 6 December 2016; Accepted: 12 June 2017; Published: 16 June 2017

**Abstract:** Most marine sponges are known to produce a large array of low molecular-weight metabolites which have applications in the pharmaceutical industry. The production of so-called specialized metabolites may be closely related to environmental factors. In this context, assessing the contribution of factors like temperature, nutrients or light to the metabolomes of sponges provides relevant insights into their chemical ecology as well as the supply issue of natural sponge products. The sponge *Crambe crambe* was chosen as a model due to its high content of specialized metabolites belonging to polycyclic guanidine alkaloids (PGA). First results were obtained with field data of both wild and farmed specimens collected in two seasons and geographic areas of the North-Western Mediterranean. Then, further insights into factors responsible for changes in the metabolism were gained with sponges cultivated under controlled conditions in an aquarium. Comparative metabolomics showed a clear influence of the seasons and to a lesser extent of the geography while no effect of depth or farming was observed. Interestingly, sponge farming did not limit the production of PGA, while ex situ experiments did not show significant effects of several abiotic factors on the specialized metabolome at a one-month time scale. Some hypotheses were finally proposed to explain the very limited variations of PGA in *C. crambe* placed under different environmental conditions.

**Keywords:** sponge; abiotic factors; specialized metabolome; metabolomics; sponge farming

## 1. Introduction

Marine sponges are known to produce a large array of small molecules also called specialized metabolites that are likely to be involved in allelopathic interactions [1–6]. From a biotechnological point of view, such metabolites are frequently investigated, mostly for pharmaceutical applications [7–9]. Several biotechnological approaches have been developed for the large-scale production of these high value-added marine products which include in situ and ex situ farming, primmorph or cell culture [10,11]. Farming explants in situ has been considered as one of the most cost-effective ways to guarantee sufficient sustainable bioactive metabolites [12]. Since spatial and temporal variations in concentrations of major specialized metabolites have been documented previously [13–15], the ideal farming location that promotes their biosynthesis as well as the harvest season have to be assessed.

<footnote n="1">¹ Université Côte d'Azur, CNRS, OCA, IRD, Géoazur, 250 rue Albert Einstein, 06560 Valbonne, France</footnote>

<footnote n="2">² Dipartimento di Scienze della Terra, dell'Ambiente e della Vita, Università di Genova, Corso Europa 26, 16132 Genoa, Italy; ericaperino@yahoo.it (E.P.); pronzato@dipteris.unige.it (R.P.)</footnote>

<footnote n="3">³ Dipartimento di Scienze della Natura e del Territorio, Università di Sassari, Via Muroni 25, 07100 Sassari, Italy; r.manconi@uniss.it</footnote>

<footnote n="4">⁴ Marine Biodiscovery, School of Chemistry, National University of Ireland Galway, University Road, H91 TK33 Galway, Ireland</footnote>

<footnote>* Correspondence: eva.ternon@unice.fr (E.T.); olivier.thomas@nuigalway.ie (O.P.T.); Tel.: +353-(0)91493563 (O.P.T.)</footnote>
The natural variations in the production of the specialized metabolites have been mostly explained by environmental and ecological factors like biotic interactions [16–19], abiotic factors [13,20–22] or more rarely, the organism's life cycle [15,23] as well as physiological state [24,25]. Overall, the effects of abiotic factors on the production of specialized metabolites have produced contrasting results so far [20,26–29]. As highly effective filter feeders, marine sponges are able to take up a large range (<0.2–50 µm) of inorganic and organic particles from seawater [30]. Therefore, the environmental conditions, such as the seawater temperature or the availability of nutrients, are expected to influence the sponge metabolism, and assessment of their impact on the biosynthesis and concentrations of specialized metabolites is of high relevance in the context of sponge farming.

In contrast to previous studies where a targeted analysis of a small number of identified compounds was performed, comparative metabolomics allows the monitoring of environmental changes using a broader range of metabolites [31]. Metabolomics refers to the set of low molecular-weight molecules (<1.5 kDa) biosynthesized by the cells in association with its metabolism, defined as the metabolome [32]. Providing a snapshot of the metabolic phenotype at a given time, metabolomics has the potential to determine biomarkers indicative of the effects of abiotic stressors such as natural and anthropogenic factors [33]. Therefore, the comparison of metabolomes for the investigation of organism-environmental interactions using high resolution analytical technology constitutes a powerful tool. In this study, we used a metabolomics approach using UHPLC-HRMS to assess the influence of seasonality and geographical distribution on the specialized metabolome of the Mediterranean encrusting sponge *Crambe crambe* (Porifera, Demospongiae) for wild populations and in situ farms. In a second set of experiments performed ex situ, we investigated the environmental factors likely to be responsible for the detected variations in situ (temperature, light, and nutrients availability).

The specialized metabolome of *C. crambe* has been well characterized during the past decades. It is composed of two families of polycyclic guanidine alkaloids (PGA, Supplemental information Figure 1); the crambescins (one or two cycles) and crambescidins (five cycles) that both exhibit substantial biological activities [34–39]. From *C. crambe*, several crambescidins (named from their molecular weight i.e., 800, 816 and 830) have been isolated and characterized [36,40,41] together with three sub-families of crambescins A–C differing in the presence of a pyrrolidine ring (A), a spiroaminal (B) and a linear 3-hydroxypropyl chain (C) [37,40,42]. We recently showed that crambescidins might form a chemical shield around the sponge that contributes to the sponge's protection and ecological success [6]. Alteration of this chemical shield due to changes in environmental conditions will affect the allelopathic interactions that the sponge maintains with its habitat.

## 2. Materials and Methods

* Collection of wild and farmed sponges

Triplicates of wild specimens of *C. crambe* were collected at different depths (5, 10 and 25 m) by SCUBA diving in the Ligurian Sea (Tigullio Gulf Punta del Faro 44°29′84.83" N, 9°21′86.02" E) in April (spring) and September 2013 (autumn) and in the Sardinia Sea (Porto Conte Bay, 40°35′43.69" N, 8°12′31.44" E) (Figure 1(A-a,B-c)) in April 2013 (spring). Autumn and spring correspond respectively to the assumed maximum and minimum production of guanidine alkaloids by *C. crambe* [17]. Farmed specimens of *C. crambe* were collected by SCUBA diving in the Ligurian Sea (Tigullio Gulf, Punta Pedale 44°19′15.35" N, 9°12′57.00" E) and in the Sardinia Sea (Tramariglio Bay, 40°35′32.60" N, 8°10′11.72" E) (Figure 1(A-b,B-c)) in April and September 2013. After collection, the sponges were immediately frozen at –20 °C, freeze-dried and maintained at –20 °C until extraction.

Mariculture of *C. crambe* was performed using the methodology USAMA® [43] using underwater modular apparatus composed of different structures: (1) a stainless steel plant, composed of 1.5 × 1.5 m² modules that can be combined in different ways (Supplementary Materials Figure S2A,B). Several pivots are used to settle, at equal distance, the ropes with the sponge explants (necklaces); (2) a PVC modular plant, made by rectangular or square forms (Supplementary Materials Figure S2C,D) anchored by ropes to structures. Pivots are inserted to anchor the sponge necklaces. The sponge
explants were assembled in the 'sponge necklaces' [44]: (i) travertine tile method, composed of 10 cm × 10 cm travertine tiles anchored to the ropes and the rope drilled into the corner of rocky tablets. The sponge explants were pasted on each tile using a non-toxic underwater bi-component epoxy putty (SubCoat, Veneziani Yacht Systems, Italy) (Supplementary Materials Figure S2E); (ii) nylon mesh method, made of a posh-like net with 2 cm meshes fixed to ropes. The sponge explants were enclosed in the meshes (Supplementary Materials Figure S2F).

<figure><img src="image_4.png" /><figcaption>**Figure 1.** Collection sites of *Crambe crambe* in the Ligurian Sea (A) (a Punta del Faro, b Punta Pedale); and Sardinia Sea (B) (c Capo Caccia).</figcaption></figure>

The abiotic sea parameters were obtained by the multiparameter probe 'Idronaut' in the Portofino Promontory area from 0 to 35 m (Povero pers. comm.). At both the Liguria and Sardinia sites, water temperatures were recorded by means of an underwater HOBO® Data Logger (Onset, MA, USA) installed in the plants and in the wild. The temperature was recorded every 6 h. Monthly average temperature values were then considered.

* **Aquarium experiments**

*Crambe crambe* specimens were sampled in the bay of Villefranche-sur-mer (20 m depth along the rocky coast) in October 2013 (Temperature & Light experiment) and in January 2014 (Nutrient experiment). Specimens were carefully transferred to the Observatory of Villefranche-sur-mer and acclimated for one month in aquaria supplied by continuous seawater flow (2 L·min⁻¹) pumped at a 5 m depth in the bay of Villefranche-sur-mer with a natural light/dark cycle as described previously [45]. Each sponge specimen was cut into small replicates of similar size (~5 cm²) in order to perform the two following sets of experiments: (1) "Temperature and Light experiment" and (2) "Nutrient experiment" that investigated the effects of variation in those environmental parameters on the specialized metabolism of *C. crambe*. Being collected from an encrusted rocky substrate, the specimens were carefully cut into different pieces with a hammer and a chisel to produce genetically identical biological replicates.

## 2.1. *Temperature and Light Experiment*

This first experiment was performed in November 2013, using a running seawater flow rate of 2 L·min⁻¹. The five following treatments were implemented: (i) control (T = 16 °C; natural light cycle), (ii) light 24/24 h (300 lux), (iii) dark 24/24 h (0 lux), (iv) 14 °C, and (v) 20 °C. The collection of sponge specimens was performed after 0, 7 and 30 days, in each aquarium and in triplicate. Water temperature
and light intensity were recorded in all aquaria by means of an underwater HOBO® Data Logger (Onset, MA, USA). The probes were set to record temperature 2 times per day (every 12 h).

## 2.2. *Nutrient Experiment*

The second experiment was conducted in March 2014, using running seawater at a flow rate of 0.7 L·min⁻¹. Five aquaria each containing 4 sponge replicates were set up: 2 control aquaria with only running seawater and 3 treatment aquaria with running seawater plus the addition of a nutrient mixture. The nutrient mixture was made up of MQ water and inorganic salts, i.e., phosphate (KH₂PO₄) and nitrate (NaNO₃), and was brought to the three treated aquaria by a peristaltic pump at a flow rate of 1 mL·min⁻¹ in order to reach a constant concentration of phosphate (PO₄³⁻ = 2.5 µmol·L⁻¹) and nitrate (NO₃⁻ = 18 µmol·L⁻¹). These concentrations were chosen as they are consistent with the anthropogenized site of Cortiou, close to the city of Marseille in the Northwestern Mediterranean Sea [46]. At T0 and after 5, 11, 18 and 24 days, sponge replicates were collected from each aquarium, and seawater was sampled from one control and one treated aquaria. Seawater samples were spiked with mercuric chloride and kept at +4 °C. Dissolved inorganic salts (PO₄³⁻, NO₃⁻, NO₂⁻ and SiOH₄) in seawater were then analyzed by a colorimetric method, using a Continuous Flow Automated Analyzer Technicon AutoAnalyzer II following a method developed previously [47]. For each set of experiments, sponge samples were flash-frozen using liquid nitrogen and maintained at -26 °C until chemical treatment.

* **Metabolomic analysis**

Methanol, acetonitrile and dichloromethane (Chromasolv®, analytical grade), formic acid (puriss. P.a. ~98%) and ammonium formate (LC-MS Ultra for UHPLC-MS) were purchased from Sigma Aldrich (Saint-Quentin Fallavier, France). Ultrapure water was prepared using a Milli-Q water system (Millipore Ltd., Billerica, MA, USA).

All sponge samples were freeze-dried and ground to obtain a dry powder which was extracted three times with a 15 mL mixture of MeOH/CH₂Cl₂ (1:1, *v*/*v*) in an ultrasonic bath (35 kHz) at room temperature. Due to the encrusting growth form of the sponge *C. crambe*, the rocky substrate was burnt at 550 °C for 24 h to obtain the sponge organic weight (Total weight - rock weight = sponge weight). Crude extracts were filtered through 8 µm paper filter (Whatman Grade 2V), dried using a rotary evaporator and further fractionated by solid phase extraction over a 500 mg C18 silica SPE cartridge (Phenomenex, Torrance, CA, USA) with a step gradient of H₂O (20 mL), MeOH (20 mL), CH₂Cl₂ (20 mL). The methanol fraction was dried using an SPD111 SpeedVac (Thermosavant, Model RH12-28), concentrated to 10 mg·mL⁻¹, and stored at -26 °C until further dilution prior to LC-MS analysis.

*UHPLC-HRMS analysis*

Prior to the UHPLC-HRMS analysis, all samples were solubilized in 1 mL of MeOH and were diluted 1/1000. In addition, five Quality Control (QC) samples were prepared by combining 10 µL of each sample in a 2 mL vial and these pools were injected every seven samples to allow chromatogram alignment during data treatment. In-line UHPLC-UV-HRMS analysis was performed using a Dionex system Ultimate 3000 equipped with an autosampler and a Dionex Ultimate 3000 diode array detector (210 and 280 nm detection), connected to a qToF mass spectrometer with an electrospray ionization interface (Bruker Impact II). Mass spectra were recorded in the positive mode. UHPLC separation was achieved on an analytical Nucleodur PolarTec column (100 × 2 mm, 1.8 µm, Macherey Nagel) using a linear elution gradient of H₂O/CH₃CN/formic acid to which was added 10 mM of ammonium formate from 80:20:0.1 (*v*/*v*/*v*, isocratic from 0 to 2 min) to 40:60:0.1 (*v*/*v*/*v*, isocratic from 8 to 10 min), at a flow rate 0.45 mL·min⁻¹ for a total of 14 min. The injected volume was set at 10 µL and detection at 280 nm. The mass spectrometer analyzer parameters were set as follows: nebulizer sheath gas, N₂ (2.1 bar); dry gas, N₂ (8 L/min); capillary temperature, 200 °C; capillary voltage, 2500 V; end plate offset, 500 V; collision gas, He; collision energy, 4 eV. Data were acquired in the 50 to 1200 m/z range.
UHPLC-HRMS *data processing for untargeted analyses*

Following UHPLC-ESI-HRMS data acquisition, MS chromatograms (Base Peak chromatogram—BPC) were exported as line spectra and converted into the netCDF file format to process the data in centroid mode with XCMS [48] using R software (version 3.2.2.). The XCMS approach involved several steps necessary to generate a final matrix: (1) Peak picking (peakwidth = c(2, 20), ppm = 10) without threshold prefilter [49], (2) retention time correction (method = "obiwarp"), (3) grouping (bw = 10, minfrac = 0.3, minsamp = 1), (4) Fillpeaks using a parameter of "noise level" that sets the minimum intensity for a centroid data point to be considered as part of a peak (noise level = 10³), and finally (5) report and data matrix generation (ions/Retention time × sample) that was exported using Microsoft Excel. Two successive filtering steps were applied to the matrix in R in order to suppress data of high analytical variability: (i) signal/noise ratio > 10 between variables matching in both extracts and blank samples (considers that ions are part of the noise if their mean intensity is less than 10-fold the blank value) and (ii) coefficient of variation in the intensity of the variables for the QC samples <20%. Each peak area was further normalized (according to the drift of equivalent ion of pooled samples), mean-centered (the mean ion intensity of all samples was subtracted from each sample ion intensity) and finally log-transformed to make individual features more comparable, using on-line MetaboAnalyst 3.0 [50]. Multivariate analysis was conducted (Principal components analysis (PCA) and Partial Least Squares-Discriminant Analysis (PLS-DA)) and PCA was further chosen since it gave the best discrimination. Very important peaks (VIPs) were determined according to the PCA loading plots.

*Metabolomic targeted analysis*

Chemical identification was performed only for compounds of interest (crambescins and crambescidins, Supplementary Materials Figure S2) that have been extensively described. The chromatogram obtained between 5.5 and 8 min of retention time was considered our metabolomic window, and the major compounds of this area were thoroughly identified by comparing their MS pattern with a standard (Supplementary Materials Table S1) [5,37]. Because several ion adducts can be detected for a single compound with electrospray ionization, the following mono- and di-charged adducts were also taken into account in the metabolite annotation: [M + H]⁺, [M + Na]⁺, [M + NH₄]⁺, [M + 2H]²⁺/2, [M + 2NH₄]²⁺/2. They are the most commonly observed in positive electrospray ionization mode with ammonium formate and formic acid as eluent modifiers. The theoretical *m/z* values for typical adduct species were compared with the experimental values to ensure the identification of the compounds (once per treatment and fraction and once for the QC samples). For spectra presenting multiple adducts, we summed the area of all adducts to obtain the total intensity corresponding to one compound. Analytical blanks confirmed that no memory effects or sample contamination skewed results. The concentration of the metabolites was assessed using a calibration with standards of crambescidin-816 and crambescin-A2-462, assuming all compounds of the same family had similar ionization potential. One-way ANOVA was performed on the data to investigate the statistical differences between conditions.

**3. Results**

* In situ study with wild and farmed sponges

To investigate depth, seasonal and geographic effects on the specialized metabolome of *C. crambe*, wild specimens were collected in the Ligurian Sea during both spring and autumn and in the Sardinia Sea during spring. A first metabolomics approach showed a significant effect of the season on the metabolome of wild specimens of *C. crambe* as revealed by the Principal Component Analysis (PCA) plot (Figure 2). Variations of the metabolic content in relation with geography or depth were not statistically significant. *Crambe crambe* was also farmed in the Ligurian Sea using both mesh and tiles protocols, and samples were collected in autumn and spring. The metabolomics approach shows that
the specialized metabolome of specimens farmed are not significantly different from wild specimens, and the highest variability was attributed to seasonal changes (Figure 2). On the other hand, the farming protocols did not show significant differences in the specialized metabolome of *C. crambe*. Therefore, growing conditions and location have small impact on the content of the specialized metabolites of the sponge.

<figure><img src="image_5.png" /></figure>

<table><caption>**Figure 2.** PCA plot obtained from the metabolomics study on wild (Wi) and farmed (Fa, Tile or Mesh protocols) samples harvested at different locations (Sar = Sardinia, Lig = Ligurian) in both autumn (Aut) and spring (Sp). The table lists the five most important ions contributing to the distribution of samples in the PCA plot.</caption><thead><tr><th>m/z</th><th>RT</th><th>[M + H]<sup>+</sup></th></tr></thead><tbody><tr><td>490.3373</td><td>6.90</td><td>979.6746</td></tr><tr><td>483.3403</td><td>6.98</td><td>963.6806</td></tr><tr><td>246.2087</td><td>7.05</td><td>489.4174</td></tr><tr><td>223.1689</td><td>4.86</td><td>445.3378</td></tr><tr><td>263.2011</td><td>6.18</td><td>525.4022</td></tr></tbody></table>

The five most contributive ions to this PCA were identified (Table S1) using full scan chromatograms for validation. Interestingly, most of the ions eluted between 4.8 and 7 min, corresponding to compounds slightly more polar than crambescins and crambescidins. None of these compounds were assigned to known metabolites. Therefore, we concluded that the variability of the metabolic content is not driven by the major metabolites (i.e., crambescins and crambescidins) but rather by minor metabolites.

To better investigate the variations in the concentrations of the two major families of specialized metabolites, the crambescins and crambescidins, a targeted metabolomics analysis was undertaken. From this study, we outlined a lower concentration of both families during spring compared to autumn (39 against 391 µmol·g⁻¹ sp, *p* = 1.55 × 10⁻⁶ ***, Ligurian samples; Figure 3A). The concentration of specialized metabolites in autumn was supported by both the crambescin and crambescidin families (up to 178 and 215 µmol·g⁻¹ sp, respectively). A similar seasonal trend was shared by
farmed specimens for both the crambescins and crambescidins families (Figure 3A). Wild and farmed
samples showed almost significant (p = 0.0657) differences for both families of compounds; these
differences being mostly supported by samples harvested in autumn (93 against 204 µmol·g⁻¹ sp for
the crambescins and 169 against 284 µmol·g⁻¹ sp for the crambescidins).

<figure><img src="image_6.png" /><figcaption>**Figure 3.** Concentrations of crambescidins and crambescins in *C. crambe* in µmol·g⁻¹ sp obtained from the targeted metabolomics study for (A) wild and farmed specimens, (B) Temperature and Light experiment, and (C) Nutrient experiment (*** p < 5 × 10⁻⁴; ** p < 0.005; * p < 0.05, one-way ANOVA).</figcaption></figure>

The spatial location of *C. crambe* slightly impacted the content of crambescidins since their concentrations in Sardinia were statistically higher than those in the Ligurian Sea (45 µmol·g⁻¹ sp, *p* = 0.0338 *), a result mostly supported by the concentrations of crambescidins. It is noteworthy that the crambescins did not show any significant difference between the two regions (23 and 22 µmol·g⁻¹ sp in Liguria and Sardinia, respectively).

* Ex situ experiments

Some experiments were also conducted in aquaria to study the variability of the specialized
metabolome of *C. crambe* when submitted to diverse environmental factors.
### 3.1. *Temperature and Light Experiments*

Although we attempted to keep fixed conditions of temperature in all aquaria, this was made difficult by an open running seawater system that is affected by natural variations. The temperature of the control aquaria naturally fluctuated between 15 and 17 °C, and a maximum of one-degree variation was recorded for both the controlled "14 °C" (between 13.5 °C and 14.5 °C) and the "20 °C" (between 19.5 °C and 20.5 °C) experiments. Light intensity was set at 300 Lux (~1200 W/m²) during the 12 h of a light-dark cycle in the temperature experiment. For the light experiments, the same cycle was maintained for the control experiment, and the results of the metabolomic analyses were compared with a permanent light experiment and permanent dark experiment.

A PCA plot revealed no clear difference between control, temperature and light experiments (Figure 4). It is noteworthy that the inter-sample variability within a treatment increased with the experimental time as shown by the distribution of the replicates. To be noted, the five most contributive ions to the PCA presented mid-polarity, very close to the PGA (5.5 < RT < 6.5), but were not assigned to any known metabolites.

<figure><img src="image_7.png" /><figcaption>**Figure 4.** PCA plot obtained from an untargeted metabolomics study of the temperature/light experiment. C = control, D = Dark, L = Light, 14 °C = Temperature 14 °C and 20 °C = Temperature 20 °C. The table shows the five most important peaks that contribute to the distribution of samples on the PCA plot, along with their retention time.</figcaption></figure>
The targeted metabolomics study confirmed the absence of statistically significant differences in the concentrations of the crambescins and crambescidins of *C. crambe* specimens across the treatments. The summed concentration of the major specialized metabolites was very close to the one measured in samples harvested in the field, although the proportion of crambescidin represented up to 50–70% of the total content (Figure 3B). Differences between treatments, for instance, the concentration of crambescidins in control at T2 (increase by a factor 2.7), were mainly driven by inter-individual variability.

## 3.2. *Nutrient Experiment*

### 3.2.1. Nutrient Concentrations

Although only nitrate and phosphate were added to the aquaria, four inorganic salts were measured in the seawater during the experiment: nitrate, nitrite, phosphate and silicate (SI, Figure S3). No difference in the concentration of non-added salts was noticeable between control and treated aquaria; those salts being most probably controlled by the sponge. The nutrient mixture was efficiently added to the treated aquarium as shown by the nitrate and phosphate concentrations (23 and 3 µmol·L⁻¹ at T4, respectively). However, attempted homogeneous concentrations in treated aquaria were only reached at day 18 after important variations (PO₄³⁻: 0.7–1.5 µmol·L⁻¹; NO₃⁻: 6–13 µmol·L⁻¹). This instability may be explained by the technical issues encountered during the delivery of the nutrient mixture by the peristaltic pump.

### 3.2.2. Metabolomics Study

The distribution of the samples on the PCA showed no clear difference in the specialized metabolome of control and nutrient-treated samples (Figure 5). However, a change in the metabolome at T4 seemed to be initiated and should be considered. It is worth noticing that the most contributive ions differed in this experiment from the previous one. In this case, we were able to identify crambescin C1 with a lower concentration at T4 in nutrient-treated samples.

<figure><img src="image_8.png" /><figcaption>Figure 5. *Cont.*</figcaption></figure>
<table><thead><tr><th>m/z</th><th>RT</th><th>[M+H]+</th><th>Known</th></tr></thead><tbody><tr><td>280.1541</td><td>7.03</td><td>559.3082</td><td>no</td></tr><tr><td>408.1482</td><td>8.13</td><td>408.1482</td><td>no</td></tr><tr><td>326.2444</td><td>5.8</td><td>326.2444</td><td>no</td></tr><tr><td>358.2473</td><td>1.15</td><td>358.2473</td><td>no</td></tr><tr><td>234.1909</td><td>6.5</td><td>467.3818</td><td>yes</td></tr></tbody></table>

**Figure 5.** PCA plot obtained from an untargeted metabolomics study of the nutrients. C = control and N = nutrient-treated samples and triplicates of biological samples are represented by a colored point. The table shows the five most important ions contributing to explain the sample distribution on the PCA plot along with their retention time and their expected [M + H]$^+$.

The targeted metabolomics analysis showed that the differences observed between control and nutrient-treated samples were not statistically different for the major metabolites (Figure 3C). The concentrations measured in this experiment were very close to those measured in in situ samples.

## 4. Discussion

The Mediterranean Sea is considered an oligotrophic region where spatial distribution of biogeochemical elements [51] as well as plankton species [52] is dependent on meso-scale circulation of water masses. Variations of such biotic and abiotic parameters might contribute to variations in the production of the specialized metabolites in common organisms like sponges. This aspect is of interest as sponges are recognized as important producers of metabolites with pharmaceutical interest. We selected the Mediterranean sponge *C. crambe* as a model due to the presence of a high quantity of bioactive alkaloids named Polycyclic Guanidinic Alkaloids from the crambescin and crambescidin families.

In order to assess the variability, specimens of both wild and farmed *Crambe crambe* sponges were harvested from two distinct Mediterranean sub-basins (Ligurian Sea and Sardinia Sea), located from ~500 kilometers along a latitudinal cline (44°18′ N vs. 40°35′ N, respectively), thus differing by the intensity of the solar radiation they receive as well as by their hydrological characteristics (temperature, salinity) and primary productivity [53]. The production of specialized metabolites in wild specimens was shown to be very similar between the two sampled locations (Figure 2). Therefore, the role played by environmental parameters in the concentration of specialized metabolites is limited in this case. This limited intraspecific metabolic variability observed between the two locations might be related to the genetic diversity as previously suggested [1,33,54–57]. In accordance with previous studies [13–15,20,23,58], a marked seasonal trend of the metabolic content showed a maximum in autumn and a minimum in spring (specimens of the Ligurian Sea, Figure 2). Interestingly, the same seasonal pattern was observed in both wild and farmed specimens cultured in this area, suggesting a role of internal physiological factors as main drivers of the variations of the specialized metabolism. Such factors could either be growth [21,24,25] or reproduction [57,58]. It has been often hypothesized that the release of larvae in late summer signals the end of the energy investment in reproduction, allowing a switch towards an enhanced production of specialized metabolites during autumn and throughout the winter [17]. This trade-off between the reproduction and the production of specialized metabolites was hypothesized by [24] and [17] for *C. crambe*. Our data set also reveals a significant difference in the metabolic content between wild and farmed specimens in autumn (Figure 2), with higher concentration of the toxic crambescidin family in wild sponges (Figure 3A). Thus, a combination of internal factors acting together with biotic factors such as competition for space [1,16,18,19] may also explain the seasonal variability of the production of the specialized metabolites. Abiotic factors such as temperature or nutrient availability might trigger other processes like sponge reproduction or other constraints of primary metabolism [25,58,59], rather than directly affecting the specialized metabolome.

The metabolomics approach allowed encompassing of the variability of the whole set of specialized metabolites produced by *C. crambe* with both the minor and the major metabolites. It is
interesting to note that none of the major specialized metabolites belong to the five most contributive
ions. Minor specialized metabolites supported the variability of the metabolic content. However, since
the seasonal and geographical trends of the major specialized metabolites revealed by the targeted
analysis are very close to those observed for the minor metabolites, the variation of both contents is
most likely closely related. This observation supports the idea of a combined contribution of both
some associated bacteria and their host sponge cells. Indeed, these results are in accordance with
relatively fast changes in the microbial production of key intermediate minor compounds under
changing environmental conditions. However, the production of the major compounds PGA would be
slower as they would require some key enzymes potentially present in sponge cells.

Our assumptions of the minor role played by abiotic factors for the analysis of the variability of
the metabolic content was confirmed by the two sets of experiments conducted in aquaria, where no
relevant effects of light, temperature or nutrient availability on the whole specialized metabolome were
observed (Figure 3B,C, Figures 4 and 6). The absence of obvious links between the light intensity, the
seawater temperature and the specialized metabolome is consistent with previous studies conducted
either in situ [15,20,26,29] or in aquaria [28]. Although no relationship between the nutrient availability
and the content of the specialized metabolism was demonstrated, DIN (dissolved inorganic nitrogen)
consumption is known to be a key mechanism of nutrient acquisition in temperate sponges [60].
Since most, if not all, specialized metabolites of the sponge *C. crambe* belong to the PGA, an increase in
DIN concentrations could have been expected to enhance their biosynthesis. The relationship between
nutrient availability and specialized metabolites is believed to be more subtle than single uptake
and use, especially since there is increasing evidence that sponge micro-symbionts are the first to
consume the nutrients [61–63]. Indeed, some sponges were shown to obtain a significant portion
of their nutrients from the bacterial symbionts themselves [64]. Likewise, symbionts may somehow
be involved in the biosynthesis of specialized metabolites in several sponge species. If we recently
observed that crambescins and crambescidins are stored in sponge cells [6], their biosynthesis could
be at least be partly inferred to be related to some microbial symbionts. Betaproteobacteria, recently
identified as major symbiont in *C. crambe*, contain a long-chain fatty acid ligase [65] that is likely to be
involved in the biosynthesis of the upper alkyl chain of the crambescins [5]. A biosynthetic pathway
for both families of compounds, implying several biosynthetic intermediates, e.g., the crambescin
A1 carboxylic acid intermediate, was proposed (Figure 6). The corresponding [M + H]⁺ = 308.2347
was found after extraction of the base peak chromatogram within the same RT (6.14 min) of the
crambescins A1 (462). Although variations in the concentration of this assumed intermediate was
found across the seasons, following the trend observed for crambescin A (data not shown), the lack of
clear identification [66] does not allow thorough conclusions. Additions of DIP and DIN could have
eventually stimulated the production of first biosynthetic links through microbial use.

<figure><img src="image_9.png" /><figcaption>**Figure 6.** Proposition of a biosynthetic intermediate for crambescin A1.</figcaption></figure>

Previous observations allowed us to propose a hypothesis for the low variations in the PGA
concentrations under different experimental conditions. Indeed, the major PGA are likely to be stored
at high concentration in spherulous cells [6], causing changes in the absolute concentrations of these
compounds that are relatively non-significant over a short-term period. Small inputs of compounds in
a large pool of stored ones will not strongly affect their final concentrations. We therefore expect that longer term experiments (months) would enable a clearer view of changes in the overall concentrations of PGA, but their stability at a shorter time (weeks) is remarkable. This storage buffering effect would also comprise a long-term defense strategy, especially against pathogenic microbes.

**5. Conclusions**

The combination of the untargeted and targeted metabolomics approaches revealed that the specialized metabolome of *C. crambe* varies in two periods of the year, and those variations are not related to abiotic factors (temperature, light and nutrients). Biotic factors, like the reproduction cycle, but also additional abiotic factors other than light, temperature or nutrients, may contribute to the seasonal variability in the metabolome.

Ex situ experiments confirmed the limited impact of abiotic factors on the metabolic content and revealed stability in the concentrations of the major specialized metabolites at a short time scale. We assume that this stability could be related to a buffering effect due to the storage of compounds in spherulous cells, which would mitigate the metabolic changes.

In the context of sponge-based biotechnologies, farming or ex situ cultures of *C. crambe* did not improve the biosynthesis of its specialized metabolites but did not significantly lower it. Since in situ cultivation of *C. crambe* by standardized methodology in sea-based farming structures enhances its mean growth (reaches around ten times the initial volume per year [44]; Manconi, pers. obs.), an increase in the specialized metabolite content over the years could be expected, promoting this cultivation protocol as a sustainable way to harvest specialized metabolites with low environmental impacts. Combining this process with our sustainable extraction of PGA from *C. crambe* [45] may represent a highly efficient process for the medium-scale production of PGA.

**Supplementary Materials:** The following are available online at www.mdpi.com/1660-3397/15/6/181/s1.

**Acknowledgments:** This study was funded by European funding within the BAMMBO project (H2020, www.bammbo.org sustainable production of Biologically Active Marine Metabolite Based Origin). We warmly thank Ornella Passafiume (OOV) for her help with the nutrient analysis as well as Fabio Damiano Ledda and Andrea Padiglia (University of Sassari) for their participation with the mariculture. We also acknowledge David Luquet and EMBRC for the collection of sponge specimens and their maintenance in aquaria at the OOV. The metabolomics study was supported by the MALLABAR platform (Marseille, France).

**Author Contributions:** All authors participated to the experimental design and conception; E.T. and E.P. performed the experiments; E.T. and E.P. analyzed the data; E.T. and O.P.T. wrote the paper.

**Conflicts of Interest:** The authors declare no conflict of interest.

**References**

1. Becerro, M.A.; Turon, X.; Uriz, M.J. Multiple Functions for Secondary Metabolites in Encrusting Marine Invertebrates. *J. Chem. Ecol.* **1997**, *23*, 1527–1547. [CrossRef]

2. Assmann, M.; Lichte, E.; Pawlik, J.R.; Kck, M. Chemical defenses of the Caribbean sponges *Agelas wiedenmayeri* and *Agelas conifera*. *Mar. Ecol. Prog. Ser.* **2000**, *207*, 255–262. [CrossRef]

3. Engel, S.; Pawlik, J.R. Allelopathic activities of sponge extracts. *Mar. Ecol. Prog. Ser.* **2000**, *207*, 273–281. [CrossRef]

4. Loh, T.-L.; Pawlik, J.R. Chemical defenses and resource trade-offs structure sponge communities on Caribbean coral reefs. *Proc. Natl. Acad. Sci. USA* **2014**, *111*, 4151–4156. [CrossRef] [PubMed]

5. Genta-Jouve, G.; Croué, J.; Weinberg, L.; Cocandeau, V.; Holderith, S.; Bontemps, N.; Suzuki, M.; Thomas, O.P. Two-dimensional ultra high pressure liquid chromatography quadrupole/time-of-flight mass spectrometry for semi-targeted natural compounds identification. *Phytochem. Lett.* **2014**, *10*, 318–323. [CrossRef]

6. Ternon, E.; Zarate, L.; Chenesseau, S.; Croué, J.; Dumollard, R.; Suzuki, M.T.; Thomas, O.P. Spherulization as a process for the exudation of chemical cues by the encrusting sponge *C. crambe*. *Sci. Rep.* **2016**, *6*. [CrossRef] [PubMed]

7. Braekman, J.-C.; Daloze, D. Chemical and biological aspects of sponge secondary metabolites. *Phytochem. Rev.* **2004**, *3*, 275–283. [CrossRef]
8. Wijffels, R.H. Potential of sponges and microalgae for marine biotechnology. *Trends Biotechnol. 2008*, **26**, 26–31. [CrossRef] [PubMed]

9. Sipkema, D.; Franssen, M.C.R.; Osinga, R.; Tramper, J.; Wijffels, R.H. Marine Sponges as Pharmacy. *Mar. Biotechnol. 2005*, **7**, 142. [CrossRef] [PubMed]

10. Koopmans, M.; Martens, D.; Wijffels, R.H. Towards Commercial Production of Sponge Medicines. *Mar. Drugs* **2009**, *7*, 787–802. [CrossRef] [PubMed]

11. Sipkema, D.; Yosef, N.; Adamczewski, M.; Osinga, R.; Mendola, R.; Tramper, J.; Wijffels, R.H. Large-scale production of pharmaceuticals by marine sponges: Sea, cell, or synthesis? *Biotechnol. Bioeng.* **2005**, *90*, 201–222. [CrossRef] [PubMed]

12. Duckworth, A.; Battershill, C. Sponge aquaculture for the production of biologically active metabolites: The influence of farming protocols and environment. *Aquaculture* **2003**, *221*, 311–329. [CrossRef]

13. Page, M.; West, L.; Northcote, P.; Battershill, C.; Kelly, M. Spatial and Temporal Variability of Cytotoxic Metabolites in Populations of the New Zealand Sponge *Mycale hentscheli*. *J. Chem. Ecol.* **2005**, *31*, 1161–1174. [CrossRef] [PubMed]

14. Ivanisevic, J.; Thomas, O.P.; Pedel, L.; Pénez, N.; Ereskovsky, A.; Culioli, G.; Pérez, T. Biochemical Trade-Offs: Evidence for Ecologically Linked Secondary Metabolism of the Sponge Oscarella balibaloi. *PLoS ONE* **2011**, *6*, e28059. [CrossRef] [PubMed]

15. De Caralt, S.; Agell, G.; Uriz, M.-J. Sources of Secondary Metabolite Variation in *Dysidea avara* (Porifera: Demospongiae): The Importance of Having Good Neighbors. *Mar. Drugs* **2013**, *11*, 489–503. [CrossRef] [PubMed]

16. Becerro, M.A.; Turon, X.; Uriz, M.J. Natural variation of toxicity in encrusting sponge *Crambe crambe* (Schmidt) in relation to size and environment. *J. Chem. Ecol.* **1995**, *21*, 1931–1946. [CrossRef] [PubMed]

17. Becerro, M.A.; Uriz, M.J.; Turon, X. Chemically-mediated interactions in benthic organisms: The chemical ecology of *Crambe crambe* (Porifera, Poecilosclerida). In *Interactions and Adaptation Strategies of Marine Organisms*; Naumov, A.D., Hummel, H., Sukhotin, A.A., Ryland, J.S., Eds.; Springer: Dordrecht, The Netherlands, 1997; pp. 77–89. [CrossRef]

18. Luter, H.M.; Duckworth, A.R. Influence of size and spatial competition on the bioactivity of coral reef sponges. *Biochem. Syst. Ecol.* **2010**, *38*, 146–153. [CrossRef]

19. Luter, H.M.; Duckworth, A.R.; Syms, C. Cytotoxic and anti-microbial activity of the sponge *Iotrochota* sp. as a function of size and spatial competitors. *Mar. Biol. Res.* **2007**, *3*, 312–318. [CrossRef]

20. Abdo, D.A.; Motti, C.A.; Battershill, C.N.; Harvey, E.S. Temperature and Spatiotemporal Variability of Salicylihalamide A in the Sponge *Haliclona* sp. *J. Chem. Ecol.* **2007**, *33*, 1635–1645. [CrossRef] [PubMed]

21. Ferretti, C.; Vacca, S.; De Ciucis, C.; Marengo, B.; Duckworth, A.R.; Manconi, R.; Pronzato, R.; Domenicotti, C. Growth dynamics and bioactivity variation of the Mediterranean demosponges *Agelas oroides* (Agelasida, Agelasidae) and *Petrosia ficiformis* (Haplosclerida, Petrosiidae). *Mar. Ecol.* **2009**, *30*, 327–336. [CrossRef]

22. Turon, X.; Martí, R.; Uriz, M.J. Chemical bioactivity of sponges along an environmental gradient in a Mediterranean cave. *Sci. Mar.* **2009**, *73*, 387–397. [CrossRef]

23. Ivanisevic, J.; Pérez, T.; Ereskovsky, A.V.; Barnathan, G.; Thomas, O.P. Lysophospholipids in the Mediterranean Sponge *Oscarella tuberculata*: Seasonal Variability and Putative Biological Role. *J. Chem. Ecol.* **2011**, *37*, 537. [CrossRef] [PubMed]

24. Turon, X.; Tarjuelo, I.; Uriz, M.J. Growth dynamics and mortality of the encrusting sponge *Crambe crambe* (Poecilosclerida) in contrasting habitats: Correlation with population structure and investment in defence. *Funct. Ecol.* **1998**, *12*, 631–639. [CrossRef]

25. Leong, W.; Pawlik, J. Evidence of a resource trade-off between growth and chemical defenses among Caribbean coral reef sponges. *Mar. Ecol. Prog. Ser.* **2010**, *406*, 71–78. [CrossRef]

26. Becerro, M.A.; Paul, V.J. Effects of depth and light on secondary metabolites and cyanobacterial symbionts of the sponge *Dysidea granulosa*. *Mar. Ecol. Prog. Ser.* **2004**, *280*, 115–128. [CrossRef]

27. Sacristán-Soriano, O.; Banaigs, B.; Becerro, M.A. Temporal Trends in the Secondary Metabolite Production of the Sponge *Aplysina aerophoba*. *Mar. Drugs* **2012**, *10*, 677–693. [CrossRef] [PubMed]

28. Duckworth, A.; West, L.; Vansach, T.; Stubler, A.; Hardt, M. Effects of water temperature and pH on growth and metabolite biosynthesis of coral reef sponges. *Mar. Ecol. Prog. Ser.* **2012**, *462*, 67–77. [CrossRef]
29. Thoms, C.; Horn, M.; Wagner, M.; Hentschel, U.; Proksch, P. Monitoring microbial diversity and natural product profiles of the sponge *Aplysina cavernicola* following transplantation. *Mar. Biol.* **2003**, *142*, 685-692.
[CrossRef]

30. Maldonado, M.; Ribes, M.; van Duyl, F.C. *Nutrient Fluxes through Sponges: Biology, Budgets, and Ecological Implications*; Academic Press: Cambridge, MA, USA, 2012.

31. Kuhlisch, C.; Pohnert, G. Metabolomics in chemical ecology. *Nat. Prod. Rep.* **2015**, *32*, 937-955.
[CrossRef]
[PubMed]

32. Griffiths, W.; Karu, K.; Hornshaw, M.; Woffendin, G.; Wang, Y. Metabolomics and metabolite profiling: Past heroes and future developments. *Eur. J. Mass Spectrom.* **2007**, *13*, 45-50.
[CrossRef]
[PubMed]

33. Bundy, J.G.; Davey, M.P.; Viant, M.R. Environmental metabolomics: A critical review and future perspectives. *Metabolomics* **2009**, *5*, 3.
[CrossRef]

34. Aoki, S.; Kong, D.; Matsui, K.; Kobayashi, M. Erythroid Differentiation in K562 Chronic Myelogenous Cells Induced by Crambescidin 800, a Pentacyclic Guanidine Alkaloid. *Anticancer Res.* **2004**, *24*, 2325-2330.
[PubMed]

35. Aron, Z.D.; Pietraszkiewicz, H.; Overman, L.E.; Valeriote, F.; Cuevas, C. Synthesis and anticancer activity of side chain analogs of the crambescidin alkaloids. *Bioorg. Med. Chem. Lett.* **2004**, *14*, 3445-3449.
[CrossRef]
[PubMed]

36. Berlinck, R.G.; Braekman, J.C.; Daloze, D.; Bruno, I.; Riccio, R.; Ferri, S.; Spampinato, S.; Speroni, E. Polycyclic guanidine alkaloids from the marine sponge *Crambe crambe* and Ca++ channel blocker activity of crambescidin 816. *J. Nat. Prod.* **1993**, *56*, 1007-1015.
[CrossRef]
[PubMed]

37. Bondu, S.; Genta-Jouve, G.; Leiros, M.; Vale, C.; Guigonis, J.-M.; Botana, L.M.; Thomas, O.P. Additional bioactive guanidine alkaloids from the Mediterranean sponge *Crambe crambe*. *RSC Adv.* **2012**, *2*, 2828-2835.
[CrossRef]

38. Martín, V.; Vale, C.; Bondu, S.; Thomas, O.P.; Vieytes, M.R.; Botana, L.M. Differential Effects of Crambescins and Crambescidin 816 in Voltage-Gated Sodium, Potassium and Calcium Channels in Neurons. *Chem. Res. Toxicol.* **2013**, *26*, 169-178.
[CrossRef]
[PubMed]

39. Rubiolo, J.A.; Ternon, E.; Lopez-Alonso, H.; Thomas, O.P.; Vega, F.V.; Vieytes, M.R.; Botana, L.M. Crambescidin-816 Acts as a Fungicidal with More Potency than Crambescidin-800 and -830, Inducing Cell Cycle Arrest, Increased Cell Size and Apoptosis in Saccharomyces cerevisiae. *Mar. Drugs* **2013**, *11*, 4419-4434.
[CrossRef]
[PubMed]

40. Jares-Erijman, E.A.; Sakai, R.; Rinehart, K.L. Crambescidins: New antiviral and cytotoxic compounds from the sponge *Crambe crambe*. *J. Org. Chem.* **1991**, *56*, 5712-5715.
[CrossRef]

41. Jares-Erijman, E.A.; Ingrum, A.L.; Carney, J.R.; Rinehart, K.L.; Sakai, R. Polycyclic guanidine-containing compounds from the mediterranean sponge *Crambe crambe*: The structure of 13,14,15-isocrambescidin 800 and the absolute stereochemistry of the pentacyclic guanidine moieties of the crambescidins. *J. Org. Chem.* **1993**, *58*, 4805-4808.
[CrossRef]

42. Berlinck, R.G.; Braekman, J.C.; Daloze, D.; Bruno, I.; Ruccio, R.; Rogeau, D.; Amade, P. Crambines C1 and C2: Two further ichthyotoxic guanidine alkaloids from the sponge *Crambe crambe*. *J. Nat. Prod.* **1992**, *55*, 528-532.
[CrossRef]
[PubMed]

43. Tipologie di Impianto Modulare per la Spongicoltura Subacquea Anche in Policoltura U.S.A.M.A. (Underwater Sponge Aquaculture Modular System). Available online: https://iris.unige.it/handle/11567/243308#.WUNwWmdLeUk (accessed on 16 June 2017).

44. Perino, E. Sustainable Production of Bioactive Compounds from Sponges. Ph.D. Thesis, Universita' degli studi di Genova, Genova, Italy, 2014; p. 187.

45. Pérez-López, P.; Ternon, E.; Gonzalez-Garcia, S.; Genta-Jouve, G.; Feijoo, G.; Thomas, O.P.; Moreira, M.T.L. Environmental solutions for the sustainable production of bioactive natural products from the marine sponge *Crambe crambe*. *Sci. Total Environ.* **2014**, *475*, 71-82.
[CrossRef]
[PubMed]

46. Tedetti, M.; Longhitano, R.; Garcia, N.; Guigue, C.; Ferretto, N.; Goutx, M. Fluorescence properties of dissolved organic matter in coastal Mediterranean waters influenced by a municipal sewage effluent (Bay of Marseilles, France). *Environ. Chem.* **2012**, *9*, 438-449.
[CrossRef]

47. Kerouel, R.; Aminot, A. Procédure Optimisée Hors-contaminations pour l'Analyse des éléments nutritifs dissous dans l'eau de mer. *Mar. Environ. Res.* **1987**, *22*, 19-32.
[CrossRef]
48. Smith, C.A.; Want, E.J.; O'Maille, G.; Abagyan, R.; Siuzdak, G. XCMS: Processing Mass Spectrometry Data for Metabolite Profiling Using Nonlinear Peak Alignment, Matching, and Identification. *Anal. Chem.* **2006**, *78*, 779-787. [CrossRef] [PubMed]

49. Patti, G.J.; Tautenhahn, R.; Siuzdak, G. Meta-analysis of untargeted metabolomic data from multiple profiling experiments. *Nat. Protoc.* **2012**, *7*, 508–516. [CrossRef] [PubMed]

50. Xia, J.; Sinelnikov, I.V.; Han, B.; Wishart, D.S. MetaboAnalyst 3.0—Making metabolomics more meaningful. *Nucleic Acids Res.* **2015**, *43*, W251–W257. [CrossRef] [PubMed]

51. Lévy, M. Mesoscale variability of phytoplankton and of new production: Impact of the large-scale nutrient distribution. *J. Geophys. Res. Oceans* **2003**, *108*, 3358. [CrossRef]

52. Lafuente, J.G.; Garcia, A.; Mazzola, S.; Quintanilla, L.; Delgado, J.; Cuttita, A.; Patti, B. Hydrographic phenomena influencing early life stages of the Sicilian Channel anchovy. *Fish. Oceanogr.* **2002**, *11*, 31–44. [CrossRef]

53. Bosc, E.; Bricaud, A.; Antoine, D. Seasonal and interannual variability in algal biomass and primary production in the Mediterranean Sea, as derived from 4 years of SeaWiFS observations. *Glob. Biogeochem. Cycles* **2004**, *18*, GB1005. [CrossRef]

54. Ivanišević, J.; Thomas, O.P.; Lejeusne, C.; Chevaldonné, P.; Pérez, T. Metabolic fingerprinting as an indicator of biodiversity: Towards understanding inter-specific relationships among Homoscleromorpha sponges. *Metabolomics* **2011**, *7*, 289–304. [CrossRef]

55. López-Legentil, S.; Turon, X. Population genetics, phylogeography and speciation of Cystodytes (Ascidiaceae) in the western Mediterranean Sea. *Biol. J. Linn. Soc.* **2006**, *88*, 203–214. [CrossRef]

56. Paul, V.J.; Arthur, K.E.; Ritson-Williams, R.; Ross, C.; Sharp, K. Chemical Defenses: From Compounds to Communities. *Biol. Bull.* **2007**, *213*, 226–251. [CrossRef] [PubMed]

57. López-Legentil, S.; Turon, X. How do morphotypes and chemotypes relate to genotypes? The colonial ascidian Cystodytes (Polycitoridae). *Zool. Scr.* **2005**, *34*, 3–14. [CrossRef]

58. Turon, X.; Becerro, M.A.; Uriz, M.J. Seasonal Patterns of Toxicity in Benthic Invertebrates: The Encrusting Sponge *Crambe crambe* (Poecilosclerida). *Oikos* **1996**, *75*, 33–40. [CrossRef]

59. Coma, R.; Ribes, M. Seasonal energetic constraints in Mediterranean benthic suspension feeders: Effects at different levels of ecological organization. *Oikos* **2003**, *101*, 205–215. [CrossRef]

60. Jiménez, E.; Ribes, M. Sponges as a source of dissolved inorganic nitrogen: Nitrification mediated by temperate sponges. *Limnol. Oceanogr.* **2007**, *52*, 948–958. [CrossRef]

61. Wilkinson, C.R. Net Primary Productivity in Coral Reef Sponges. *Science* **1983**, *219*, 410–412. [CrossRef] [PubMed]

62. Ribes, M.; Jimenez, E.; Yahel, G.; Lopez-Sendino, P.; Diez, B.; Massana, R.; Sharp, J.H.; Coma, R. Functional convergence of microbes associated with temperate marine sponges. *Environ. Microbiol.* **2012**, *14*, 1224–1239. [CrossRef] [PubMed]

63. Zhang, F.; Blasiak, L.C.; Karolin, J.O.; Powell, R.J.; Geddes, C.D.; Hill, R.T. Phosphorus sequestration in the form of polyphosphate by microbial symbionts in marine sponges. *Proc. Natl. Acad. Sci. USA* **2015**, *112*, 4381–4386. [CrossRef] [PubMed]

64. Wang, G. Diversity and biotechnological potential of the sponge-associated microbial consortia. *J. Ind. Microbiol. Biotechnol.* **2006**, *33*, 545. [CrossRef] [PubMed]

65. Croué, J.; West, N.J.; Escande, M.-L.; Intertaglia, L.; Lebaron, P.; Suzuki, M. A single betaproteobacterium dominates the microbial community of the crambescidine-containing sponge *Crambe crambe*. *Sci. Rep.* **2013**, *3*, 2583. [CrossRef] [PubMed]

66. Sumner, L.W.; Amberg, A.; Barrett, D.; Beale, M.H.; Beger, R.; Daykin, C.A.; Fan, T.; Fiehn, O.; Goodacre, R.; Griffin, J.L.; et al. Proposed minimum reporting standards for chemical analysis. *Metabolomics* **2007**, *3*, 211–221. [CrossRef] [PubMed]

<figure><img src="image_10.png" /></figure>