Evaluating the Effects of an Organic Extract from the Mediterranean Sponge Geodia cydonium on Human Breast Cancer Cell Lines

Marine sponges are an excellent source of bioactive secondary metabolites for pharmacological applications. In the present study, we evaluated the chemistry, cytotoxicity and metabolomics of an organic extract from the Mediterranean marine sponge Geodia cydonium, collected in coastal waters of the Gulf of Naples. We identified an active fraction able to block proliferation of breast cancer cell lines MCF-7, MDA-MB231, and MDA-MB468 and to induce cellular apoptosis, whereas it was inactive on normal breast cells (MCF-10A). Metabolomic studies showed that this active fraction was able to interfere with amino acid metabolism, as well as to modulate glycolysis and glycosphingolipid metabolic pathways. In addition, the evaluation of the cytokinome profile on the polar fractions of three treated breast cancer cell lines (compared to untreated cells) demonstrated that this fraction induced a slight anti-inflammatory effect. Finally, the chemical entities present in this fraction were analyzed by liquid chromatography high resolution mass spectrometry combined with molecular networking.


Introduction
Marine sponges are conspicuous members of the marine benthos, occurring worldwide from polar and temperate to tropical seas and from the intertidal to deep-sea environments [1,2].They have been considered a gold mine for the discovery of marine natural products during the past 50 years, with about 4851 compounds identified to date, contributing to nearly 30% of all marine natural products discovered so far [3].Many of these products are of considerable biotechnological interest as pharmaceutical agents [4,5].Marine sponges have given rise to two commercial anticancer compounds: (i) cytarabine or Ara-C, from the Caribbean sponge Tethya crypta, approved in 1969 and still used to treat acute myelocytic leukemia and non-Hodgkin's lymphoma; and (ii) Eribulin or Halaven, a synthetic derivative based on the structure of Halichondrin B from the sponge Halichondria okadai, approved in 2010 for the treatment of drug refractory breast cancer.Another commercial drug obtained from sponges is vidarabine or Ara-A, a synthetic analog of spongouridine originally isolated in 1950 from the sponge Tethya crypta, commercialized as Vira-A ® and aciclovir (Zovirax ® ) for the treatment of herpes simplex and herpes zoster viruses infections [6].Other sponge compounds are in phase I clinical trials (Hemiasterlin E7974) against cancer, or in preclinical trials against a number of other pathologies: Geodisterol sulfates (fungal infections), Plakortin (malaria), Homogenistic acid (malaria), Hymenidin (tuberculosis), Ggyrosanols (viral infections), Dysidine (diabetes), Capnellene (inflammation), Callyspongidiol (immunity), Calyculin A (nervous system), and Dysideamine (nervous system) [7].The search for new drug candidates from sponges is therefore on the rise [8,9].Here, we explore the biological activity of the Mediterranean sponge Geodia cydonium (Porifera, Demospongiae, Astrophorida, Geodidae).In a previous study, we demonstrated that a methanol extract of G. cydonium had an anti-inflammatory effect on the human breast cancer MCF-7 cell line by decreasing the levels of several pro-inflammatory cytokines without inducing cytotoxic effects thus indicating the potential of G. cydonium against human breast cancer [10].
The present study aimed to study the effects of the organic extract of G. cydonium on human breast cancer cell lines MCF-7, MDA-MB231 and MDA-MB468, and on normal breast cell line, MCF-10A, used as a control.The extract was fractionated by chromatography, and of the three most abundant fractions obtained, only one (fraction 3) was able to block cell proliferation of all breast cancer cell lines, with no effects on the normal breast cell line.Cellular studies were then performed to verify if this fraction induced apoptosis and/or blockage of the cell cycle.Moreover, metabolomic profiling on cells treated with fraction 3 has allowed for the identification of the metabolite pathways modulated by this fraction.Anti-or pro-inflammatory effects on three breast cancer cell lines have also been tested by cytokinome evaluation on the related cellular polar fractions.The active fraction was analyzed by liquid chromatography (LC)-high resolution mass spectrometry (HRMS) and tandem mass spectrometry (HRMS/MS) together with the molecular networking technique [11].

Cell Proliferation
Sulforhodamine B (SRB) assay was used to identify the concentrations at which cell growth was inhibited by 50% in the breast cancer cell lines MCF-7, MDA-MB231 and MDA-MB468.None of the three fractions blocked cell proliferation in the normal cell line MCF-10A after 24 and 48 h of incubation (Figure S1).Moreover, two of the three fractions had no effects on all three human cancer cells after 24 and 48 h of treatment (Figures S2-S4).Interestingly, only one fraction blocked cell proliferation in all three human cancer cell lines.
In particular, MCF-7 cells showed a growth inhibition with the half minimal (50%) Inhibitory Concentration (IC 50 ) of 72 and 67 µg/mL after 24 and 48 h of treatment, respectively (Figure 1a); MDA-MB231 cells showed a growth inhibition with an IC 50 at 73 and 44 µg/mL after 24 and 48 h of treatment, respectively (Figure 1b); MDA-MB468 cells showed a growth inhibition with an IC 50 at 80 and 70 µg/mL after 24 and 48 h of treatment, respectively (Figure 1c).IC 50 values after 48 h treatment were lower compared to those obtained after 24 h.

Characterization by Liquid Chromatography-High Resolution Mass Spectrometry and Tandem Mass Spectrometry (LC-HRMS and LC-HRMS/MS) Combined with Bioinformatic Analyses (GNPS)
The active fraction was analyzed by LC-HRMS and LC-HRMS/MS (Figure 2).Data obtained were used to generate a molecular network using the Global Natural Product Social Molecular Network (GNPS).The mass spectral molecular networking resulted in a fast identification of known metabolites from natural extracts (dereplication) as well as new analogs [12].More specifically, the spectra from one or more LC-HRMS/MS runs have been compared pairwise, and each spectrum also compared with MS/MS spectra of known natural products in GNPS libraries.The two-dimensional network obtained from the active fraction from G. cydoniumis shown in Figure 2. A node represents a single chemical entity and its relatedness with other compounds present in the mixture is represented by an edge.The network contains six clusters ranging from two to thirteen nodes.Dereplication and search for analogs resulted in the identification of several known metabolites (Figure 2).

Characterization by Liquid Chromatography-High Resolution Mass Spectrometry and Tandem Mass Spectrometry (LC-HRMS and LC-HRMS/MS) Combined with Bioinformatic Analyses (GNPS)
The active fraction was analyzed by LC-HRMS and LC-HRMS/MS (Figure 2).Data obtained were used to generate a molecular network using the Global Natural Product Social Molecular Network (GNPS).The mass spectral molecular networking resulted in a fast identification of known metabolites from natural extracts (dereplication) as well as new analogs [12].More specifically, the spectra from one or more LC-HRMS/MS runs have been compared pairwise, and each spectrum also compared with MS/MS spectra of known natural products in GNPS libraries.The two-dimensional network obtained from the active fraction from G. cydonium is shown in Figure 2. A node represents a single chemical entity and its relatedness with other compounds present in the mixture is represented by an edge.The network contains six clusters ranging from two to thirteen nodes.Dereplication and search for analogs resulted in the identification of several known metabolites (Figure 2).
A more traditional approach based on the search against databases of natural and specific marine compounds (Metlin at https://metlin.scripps.edu/and MarinLit at http://pubs.rsc.org/marinlit/) was also performed when no correspondence between molecular formula and compound name was found.Molecular networking indeed correlated compounds that share similar structural features and the dereplication workflow followed the same principle; therefore, a match within GNPS libraries does not always imply identity, but gives a significant clue on the structure of the molecule.In addition, even if libraries contain a large number of MS/MS spectra of known natural products, this number is still limited.A more traditional approach based on the search against databases of natural and specific marine compounds (Metlin at https://metlin.scripps.edu/and MarinLit at http://pubs.rsc.org/marinlit/) was also performed when no correspondence between molecular formula and compound name was found.Molecular networking indeed correlated compounds that share similar structural features and the dereplication workflow followed the same principle; therefore, a match within GNPS libraries does not always imply identity, but gives a significant clue on the structure of the molecule.In addition, even if libraries contain a large number of MS/MS spectra of known natural products, this number is still limited.
More specifically, cluster (a) containing thirteen nodes was the cluster containing nucleosides and nucleobases.Dereplication indicated the presence in the cluster of methyladenosine (m/z 282) and methyladenine (m/z 150).The two nodes directly connecting to these (m/z 268 and m/z 136 respectively) were also recognized by GNPS as methyladenosine and methyladenine but with a mass difference of 14 amu (atomic mass unit) with the library compounds, corresponding to a missing methyl group.Analyses of the HRMS spectra indicated the molecular formulas as C10H13O4N5 and C5H5N5, therefore confirming the identity of the compounds as adenosine and adenine, respectively.HRMS/MS spectra of the nucleosides revealed their characteristic fragmentation pattern: cleavage of the glycosidic bonds gave protonated bases (adenine and methyladenine) with a sugar moiety as the neutral fragment [13].According to LC-HRMS, the relative abundance of methyladenosine is 20 times higher than adenosine, while methyladenine is double that of adenine.
Cluster (b) contained two nodes; the one with m/z 166 was dereplicated as phenylalanine (Phe), as confirmed also by the corresponding molecular formula.The second node (m/z 132) was identified as 2 amu different from pipecolinic acid (C6H11O2N); its molecular formula in the HRMS spectrum of C6H14O2N for the pseudomolecular ion peak [M+H] + indicated the presence of two hydrogen atoms more than pipecolinic acid, and consequently suggested that the molecule had an open ring, and was therefore leucine (Leu).The two amino acids were found in comparable amounts in the active fraction, with Phe more abundant than Leu.
In cluster (c) the node with m/z 190.0498 has been indicated as 5-aminosalycilic acid (5-asa) in GNPS.The difference of 36 amu of the molecule with 5-asa corresponds to a C3 unit.Therefore, the molecule contains three more carbon atoms than 5-asa as well as three additional unsaturations.Metlin recognized the molecule as 3-hydroxyquinaldic acid (3-HQA) that perfectly matched the More specifically, cluster (a) containing thirteen nodes was the cluster containing nucleosides and nucleobases.Dereplication indicated the presence in the cluster of methyladenosine (m/z 282) and methyladenine (m/z 150).The two nodes directly connecting to these (m/z 268 and m/z 136 respectively) were also recognized by GNPS as methyladenosine and methyladenine but with a mass difference of 14 amu (atomic mass unit) with the library compounds, corresponding to a missing methyl group.Analyses of the HRMS spectra indicated the molecular formulas as C 10 H 13 O 4 N 5 and C 5 H 5 N 5 , therefore confirming the identity of the compounds as adenosine and adenine, respectively.HRMS/MS spectra of the nucleosides revealed their characteristic fragmentation pattern: cleavage of the glycosidic bonds gave protonated bases (adenine and methyladenine) with a sugar moiety as the neutral fragment [13].According to LC-HRMS, the relative abundance of methyladenosine is 20 times higher than adenosine, while methyladenine is double that of adenine.
Cluster (b) contained two nodes; the one with m/z 166 was dereplicated as phenylalanine (Phe), as confirmed also by the corresponding molecular formula.The second node (m/z 132) was identified as 2 amu different from pipecolinic acid (C 6 H 11 O 2 N); its molecular formula in the HRMS spectrum of C 6 H 14 O 2 N for the pseudomolecular ion peak [M+H] + indicated the presence of two hydrogen atoms more than pipecolinic acid, and consequently suggested that the molecule had an open ring, and was therefore leucine (Leu).The two amino acids were found in comparable amounts in the active fraction, with Phe more abundant than Leu.
In cluster (c) the node with m/z 190.0498 has been indicated as 5-aminosalycilic acid (5-asa) in GNPS.The difference of 36 amu of the molecule with 5-asa corresponds to a C 3 unit.Therefore, the molecule contains three more carbon atoms than 5-asa as well as three additional unsaturations.Metlin recognized the molecule as 3-hydroxyquinaldic acid (3-HQA) that perfectly matched the requests.In addition, analysis of the HRMS/MS spectrum of the molecule showed the loss of water and formic acid (HCOOH) to form the ions at m/z 172.0390 and 144.0440, but not the loss of ammonia as expected by 5-asa.
Finally, in cluster (d) there were phosphatidylethanolamine (PE) and its analogs, whereas cluster (e) and cluster (f) both comprised unknown compounds.
Due to the small amount of the networked compounds, it was not possible to establish the structure of the unknown compounds.

Apoptosis Increase in MCF-7, MDA-MB231 and MDA-MB468 Cells after Treatment
Based on the results obtained above, the activity of the active fraction as an apoptosis-inducer was evaluated after treatment with IC 50 concentrations obtained after 48 h.As shown in Table 1, an increase in the number of apoptotic cells (51.2% for MCF-7, 63.1% for MDA-MB231 and 56.6% for MDA-MB468 cells) was observed compared to untreated cells used as controls.

RT-qPCR Analysis on MCF-7, MDA-MB231 and MDA-MB468 Cells after Treatment
To further elucidate the molecular mechanism through which the active fraction was able to induce apoptosis in breast cancer cells, we have examined mRNA expressions of some genes involved in the intrinsic mitochondrial pathway such as p53, Bax, p38, and caspase-3; extrinsic death receptor pathways were analyzed though the expression levels of caspase-3 and caspase-8 (Table 2), as already performed in our recent paper [14].RT-qPCR was used to detect the mRNA expression after treatment with IC 50 concentrations of the active fraction in MCF-7, MDA-MB231 and MDA-MB468 cells.Expression changes were normalized on β-actin mRNA expression (Figure 3).Results showed that the mRNA expression of p53, Bax, p38, caspase-3 and caspase 8 genes increased significantly after treatment for 48 h with the active fraction in MCF-7, MDA-MB231 and MDA-MB468 cells.Taken together, these observations indicate that the active fraction induced apoptosis through both intrinsic and extrinsic pathways.

Analyses of the Active Fraction on the Cell Cycle of MCF-7, MDA-MB231 and MDA-MB468 Cells
Evaluation of cell cycle phases was performed on all three breast cancer lines after treatment with IC50 concentrations obtained after 48 h.As shown in Table 3, our results indicated that there were no visible changes for all three human breast cancer cell lines, suggesting that the active fraction had cytotoxic effects and was able to induce apoptosis but not blockage of cell cycle progression.

Metabolomic Profiling on Breast Cancer Cells
1 H-NMR spectra were acquired on polar extracts from MCF-7, MDA-MB231 and MDA-MB468 cells before and after treatment with IC50 concentrations of the active fraction.Table 4 shows the proton resonances related to the metabolites identified in three breast cancer cells.In detail, the spectral region from 0.5 to 3 ppm contains signals from alanine, arginine, aspartate, glutamate, glutamine, isoleucine, lactate, leucine, lipids, lysine, proline, threonine and valine.The spectral region from 3 to 5.5 ppm contains mainly signals from choline, α-glucose, β-glucose, glycine, Results showed that the mRNA expression of p53, Bax, p38, caspase-3 and caspase 8 genes increased significantly after treatment for 48 h with the active fraction in MCF-7, MDA-MB231 and MDA-MB468 cells.Taken together, these observations indicate that the active fraction induced apoptosis through both intrinsic and extrinsic pathways.

Analyses of the Active Fraction on the Cell Cycle of MCF-7, MDA-MB231 and MDA-MB468 Cells
Evaluation of cell cycle phases was performed on all three breast cancer lines after treatment with IC 50 concentrations obtained after 48 h.As shown in Table 3, our results indicated that there were no visible changes for all three human breast cancer cell lines, suggesting that the active fraction had cytotoxic effects and was able to induce apoptosis but not blockage of cell cycle progression.cells before and after treatment with IC 50 concentrations of the active fraction.Table 4 shows the proton resonances related to the metabolites identified in three breast cancer cells.In detail, the spectral region from 0.5 to 3 ppm contains signals from alanine, arginine, aspartate, glutamate, glutamine, isoleucine, lactate, leucine, lipids, lysine, proline, threonine and valine.The spectral region from 3 to 5.5 ppm contains mainly signals from choline, α-glucose, β-glucose, glycine, glycero-phosphocholine and phosphocholine.The 5.5-8.5 ppm region contains the resonances of histidine, phenylalanine and tyrosine.Orthogonal Projections to Latent Structures discriminant analysis (OPLS-DA) plot indicated that the spectra obtained for the three cell lines clustered in different groups (Figure 4A); each cell line showed statistically different proton signals and metabolites after treatment compared to controls (see variable importance in projection (VIP) score plot) (Figure 4B-D).
Data showed that the level of lactate increased after treatment in all three cell lines whereas α-glucose, β-glucose, choline, glycerophosphocholine, glutamine, glutamate and lipids decreased.The level of other metabolites also decreased: proline in MCF-7 cells, threonine in MDA-MB231, asparagine and lysine in MDA-MB468 cells, whereas glycine in both MDA-MB231 and MDA-MB468 cells.
Pathway analysis showed that the metabolites modulated by the active fraction were mainly involved in glycolysis, and lipid and amino acid metabolism (Table S1).

Evaluation of Cytokine Levels in Breast Cancer Cells
In a previous study [10], we demonstrated that a methanol extract of G. cydonium had an anti-inflammatory effect on the human breast cancer MCF-7 cell line by decreasing the levels of several pro-inflammatory cytokines thereby indicating the potential of G. cydonium against human breast cancer.In that paper [10], we evaluated the levels of the following 27 cytokines (IL-1β, IL-1ra, IL-2, IL-4, IL-5, IL-6, IL-7, CCL2, CCL11, CXCL10, CXCL8, IFN-γ, IL-9, IL-10, IL-12 (p70), IL-13, IL-15, IL-17, basic FGF, G-CSF, GM-CSF, MIP-1α, MIP-1β, PDGF-ββ, RANTES, TNF-α, and VEGF) by multiplex biometric Enzyme-Linked Immunosorbent Assay (ELISA)-based immunoassay.Hence, since the active fraction in this study derives from the same methanol extract of G. cydonium [10], we evaluated the same panel of 27 cytokines that comprises thirteen pro-inflammatory and anti-inflammatory interleukins, seven chemokines, five growth factors, one interferon and one tumor necrosis factor.More specifically, cytokine levels were evaluated in MCF-7, MDA-MB231 and MDA-MB468 cellular supernatants after treatment with the half minimal (50%) Inhibitory Concentration (IC 50 ) concentrations of the active fraction after 48 h.Untreated cells were used as control.As shown in Figure 5, the levels of VEGF (vascular endothelial growth factor) and of two pro-inflammatory cytokines/chemokines, IL-8

Discussion
Whereas several studies have focused on marine sponges with the aim of identifying natural compounds with pharmacological activity, few studies have been performed on the marine sponge G. cydonium.Very old studies identified steroidal ketones from this species [15].
Our data revealed that an active fraction obtained from G. cydonium induced apoptosis, differentially increasing the number of apoptotic cells in the three different breast cancer cell lines.These differences in the percentage of apoptotic cells for the three breast cancer cell lines may be attributed to their different nature.In fact, MCF-7 cells are estrogen-receptor-positive, whereas MDA-MB231 and MDA-MB468 cells are estrogen-receptor-negative.In addition, the triple-negative human breast cancer cells, MDA-MB231 and MDA-MB468, are characterized by a more malignant phenotype [16,17] because they constitutively express mutations of the tumor protein p53, whereas MCF-7 cells have the wild-type p53.
An important outcome from our studies was that metabolomic profiling revealed different metabolites in the three breast cancer cells MCF-7, MDA-MB231 and MDA-MB468.Pathway analysis showed that the metabolites modulated by the active fraction were mainly involved in glycolysis, and lipid and amino acid metabolism.It is well known that most tumors are highly dependent on glucose to support bioenergetic and macromolecular synthesis.In fact, tumor cells have a higher rate of glucose consumption through a different glycolysis pathway compared to normal cells, in which there is the conversion of pyruvate to lactate [18].Thus, glucose consumption is an important step for cancer cells because the deprivation of glucose can induce their death [19,20].Our data show that the active sponge fraction produced a decrease in glucose levels and an increase in lactate deriving from the unbalanced conversion of glucose to lactate in all three breast cancer cells.Considering that glutamine, glycine and glutamate originate from glucose metabolism via glycolysis intermediates, a decrease in their levels may be correlated to a decrease in glucose levels.
We also observed a decrease in choline levels after treatment with sponge extract.Choline is necessary for normal membrane formation, and its abnormal metabolism is considered a metabolic hallmark associated to oncogenesis and tumor progression [21].Its decrease in the three breast cancer cell lines after treatment suggests that our active sponge fraction may be able to block cancer progression.
Different cytokines were also decreased by the sponge extract, including VEGF, CXCL10, and IL-8.VEGF is a powerful angiogenic factor and has been shown to have a role in tumor angiogenesis.In particular, it has a key role in the development, progression and metastasis of different types of cancers [22].CXCL10 is a pro-inflammatory chemokine which modulates innate and adaptive

Discussion
Whereas several studies have focused on marine sponges with the aim of identifying natural compounds with pharmacological activity, few studies have been performed on the marine sponge G. cydonium.Very old studies identified steroidal ketones from this species [15].
Our data revealed that an active fraction obtained from G. cydonium induced apoptosis, differentially increasing the number of apoptotic cells in the three different breast cancer cell lines.These differences in the percentage of apoptotic cells for the three breast cancer cell lines may be attributed to their different nature.In fact, MCF-7 cells are estrogen-receptor-positive, whereas MDA-MB231 and MDA-MB468 cells are estrogen-receptor-negative.In addition, the triple-negative human breast cancer cells, MDA-MB231 and MDA-MB468, are characterized by a more malignant phenotype [16,17] because they constitutively express mutations of the tumor protein p53, whereas MCF-7 cells have the wild-type p53.
An important outcome from our studies was that metabolomic profiling revealed different metabolites in the three breast cancer cells MCF-7, MDA-MB231 and MDA-MB468.Pathway analysis showed that the metabolites modulated by the active fraction were mainly involved in glycolysis, and lipid and amino acid metabolism.It is well known that most tumors are highly dependent on glucose to support bioenergetic and macromolecular synthesis.In fact, tumor cells have a higher rate of glucose consumption through a different glycolysis pathway compared to normal cells, in which there is the conversion of pyruvate to lactate [18].Thus, glucose consumption is an important step for cancer cells because the deprivation of glucose can induce their death [19,20].Our data show that the active sponge fraction produced a decrease in glucose levels and an increase in lactate deriving from the unbalanced conversion of glucose to lactate in all three breast cancer cells.Considering that glutamine, glycine and glutamate originate from glucose metabolism via glycolysis intermediates, a decrease in their levels may be correlated to a decrease in glucose levels.
We also observed a decrease in choline levels after treatment with sponge extract.Choline is necessary for normal membrane formation, and its abnormal metabolism is considered a metabolic hallmark associated to oncogenesis and tumor progression [21].Its decrease in the three breast cancer cell lines after treatment suggests that our active sponge fraction may be able to block cancer progression.
Different cytokines were also decreased by the sponge extract, including VEGF, CXCL10, and IL-8.VEGF is a powerful angiogenic factor and has been shown to have a role in tumor angiogenesis.In particular, it has a key role in the development, progression and metastasis of different types of cancers [22].CXCL10 is a pro-inflammatory chemokine which modulates innate and adaptive responses of the immune system and/or regulation of cell growth and angiogenesis.Since its induction is associated with several human diseases including cancer, it is currently in use for cancer target therapy [23].IL-8 is a pro-inflammatory chemokine, and its increased expression in different cell types such as endothelial cells, infiltrating neutrophils, cancer cells, and tumor-associated macrophages suggests that it may be an important regulatory factor within the tumor microenvironment which promotes proliferation, survival, and migration of endothelial, cancer cells, and infiltrating leukocites [24].
We also observed a slight increase in levels of IL-4 and IL-10.These are cytokines with anti-inflammatory activity and properties involved in antitumor response [25].Our data therefore corroborate previous findings showing that a methanolic extract from the same sponge decreased pro-inflammatory cytokine levels in the human breast cancer cell line MCF-7 [10].
The main constituents of the active sponge fraction were identified as: (i) methyladenine, methyladenosine, adenine and adenosine in cluster (a); (ii) leucine and phenylalanine in cluster (b); (iii) 3-HQA in cluster (c); and (iv) phosphatidylethanolamine in cluster (d).The metabolites identified in cluster (a) are nucleosides (methyladenine and methyladenosine) and nucleotides (adenine and adenosine) that are vital components of all living cells involved in several key biological processes (nucleic acid synthesis).Marine nucleosides have already been shown to have antiviral, anticancer, vasodilator, muscle relaxant, and hypertensive activities and are currently providing new lead compounds for drug design [26].The amino acids (leucine and phenylalanine) in cluster (b) are known to occur both in the free-state and as basic units of proteins and other metabolites (peptides) in marine organisms such as sponges and associated bacteria.Marine amino acid derivatives and peptides have been shown to possess interesting biological properties; in particular, they have attracted much attention due to their high specificity against cancer cell lines [27].3-HQA in cluster (c) is found in several bisintercalator natural products of marine origin (such as thiocoraline, triostin, SW-163 and echinomycin/quinomycin) [28], as one of the two key chromophores (together with quinoxalinic acid) for binding to duplex DNA by insertion between the bases to allow proper placement of the peptidic core into the DNA minor groove [29].Phosphatidylethanolamine (PE) in cluster (d) belongs to the family of phospholipids and glycolipids that are the main constituents of sponge cell membranes [30,31], and have been shown to possess immune-modulating and antitumor activity [32,33].
In conclusion, the presence of these compounds in our extract suggests that they could be responsible for the observed changes in the metabolic profiling and cytokine secretion observed in the three cancer cell lines.Our results represent an interesting finding considering that triple negative breast cancer cells (MDA-MB231 and MDA-MB468) represent 15% of all breast cancer cases and are associated with a high malignancy and low chances of survival compared to estrogen-positive cells.Further studies are necessary to understand the mechanism of action of the single compounds present in the active fraction from G. cydonium and to evaluate their efficacy compared to other drugs used for the treatment of breast cancer cells.It would also be interesting to evaluate the possibility of a synergistic effect of a pool of these compounds to improve the effect of chemotherapy treatments against the proliferation of breast cancer cells, as already shown for hepatocellular carcinoma [34].

Collection, Extraction, and Separation
Several samples of the sponge Geodia cydonium (Order Tetractinellida, Family Geodiidae) were collected on July 2015 at 20 m depth by scuba diving in the "Parco Sommerso di Baia" (Gulf of Naples, Italy).After collection, the samples were immediately frozen and stored at −20 • C until extraction.The sponges (509 g wet weight) were homogenized and extracted with MeOH (2 × 2 L), MeOH and CHCl 3 in different ratios (2:1, 1:1, 1:2) and then with CHCl 3 (2 × 2 L) [35].The MeOH extracts were partitioned between H 2 O and n-BuOH [36]; the BuOH layer was combined with the CHCl 3 extracts and concentrated in vacuo.The resulting organic extract (3.9 g) was chromatographed by Droplet CounterCurrent Chromatography (DCCC) using CHCl 3 /CH 3 OH/H 2 O (7:13:8) in the ascending mode; 6 mL fractions were collected and combined in ten fractions based on their similar Thin Layer Chromatography (TLC) retention times.The most abundant fractions were tested to evaluate their anti-proliferative activity on three human breast cancer cell lines, MDA-MB231, MDA-MB468 and MCF-7, and the normal human breast epithelial cell line MCF-10A.

LC-HRMS and LC-HRMS/MS and Molecular Networking Analyses
Experiments were performed using a Thermo LTQ Orbitrap XL high-resolution ElectroSpray Ionization (ESI) mass spectrometer coupled to an Agilent model 1100 LC system, which included a solvent reservoir, in-line degasser, binary pump, and refrigerated autosampler.A 5 µm Kinetex C18 column (50 × 2.1 mm), maintained at 25 • C, was operated using a gradient elution of H 2 O and MeOH running at 200 µL/min.The gradient program was as follows: 10% MeOH for 5 min, 10-100% MeOH over 25 min, 100% MeOH for 13 min.All the mass spectra were recorded in the positive-ion mode.MS parameters were a spray voltage of 4.8 kV, a capillary temperature of 285 • C, a sheath gas rate of 32 units N 2 (ca.320 mL/min), and an auxiliary gas rate of 15 units N 2 (ca.150 mL/min).Data were collected in the data-dependent acquisition (DDA) mode, in which the first and second most intense ions of a full-scan mass spectrum were subjected to tandem mass spectrometry (MS/MS) analysis.MS/MS scans were obtained for selected ions with collision induced dissociation (CID) fragmentation, isolation width 2.0, normalized collision energy 36, Activation Q 0.250, and activation time 30 ms.Mass data were analyzed using the Thermo Xcalibur software (ThermoFisher Scientific, Waltham, MA, USA).
A molecular network was created using the online workflow at GNPS.The data was then clustered with MS-Cluster with a parent mass tolerance of 2.0 Da and a MS/MS fragment ion tolerance of 0.5 Da to create consensus spectra.Further, consensus spectra that contained less than 2 spectra were discarded.A network was then created where edges were filtered to have a cosine score above 0.6 and more than 6 matched peaks.Further edges between two nodes were kept in the network if and only if each of the nodes appeared in each other's respective top 10 most similar nodes.For dereplication purposes the spectra in the network were then searched against GNPS' spectral libraries.All matches kept between network spectra and library spectra were required to have a score above 0.6 and at least 6 matched peaks.Analog search was enabled against the library with a maximum mass shift of 100.0 Da.The data were then imported into Cytoscape 3.2.1 (Available online: http://www.cytoscape.org/) and displayed as a network of nodes and edges.The network was organized with the preferred layout plug-in.

Cell Treatment and Cell Proliferation Assay
Cell proliferation of cancer cells was assessed in presence and absence of the methanol extract from G. cydonium by colorimetric assay with sulforhodamine B (SRB, Sigma-Aldrich).This extract was first dissolved in dimethyl sulfoxide (DMSO 100 mM, Sigma-Aldrich) at concentrations <0.1%, so as not to induce toxic effects on cells.Thus, a stock solution (100 mg/mL) and its serial dilutions had a final concentration of DMSO of 0.05%.
Cancer cells were plated in 96 well tissue culture plates at a concentration of 15 × 10 3 cells per well and allowed to attach for 24 h.Cells were then treated with different concentrations of the methanol extract (2.5, 5, 10, 25, 50, 100, and 200 µg/mL) and incubated for 24 and 48 h.These concentrations were selected on the basis of our recent paper [10].After 48 h of treatment, cells were fixed with trichloroacetic acid (Sigma-Aldrich) for 1 h at 4 • C. Subsequently they were stained for 30 min with 0.4% (wt/vol) sulforhodamine B (SRB, Sigma-Aldrich) dissolved in 1% acetic acid.The number of viable cells was directly proportional to the amount of protein bound-dye, which was then solubilized with 10 mM Tris base solution (pH 10.5) and measured at 540 nm using the ELISA fluorometric assay (Bio-Rad, Hercules, CA, USA; Microplate Reader).All experiments were performed in duplicate and repeated three times.The IC 50 was assessed from the dose-response curves.

Apoptosis Evaluation
After counting, 3 × 10 5 cells were harvested and washed twice with ice-cold PBS.Cells were labeled with an Annexin V and Dead Cell Assay kit according to the manufacturer's instructions (Merck Millipore, Darmstadt, Germany).The kit detects the externalization of phosphatidylserine (PS) in apoptotic cells using fluorescently-labeled Annexin V in combination with the dead cell marker 7-aminoactinomycin D (7-AAD).We identified four populations of cells: (1) viable cells that did not undergo detectable apoptosis: Annexin V (−) and dead cell marker (−); (2) early apoptotic cells: Annexin V (+) and dead cell marker (−); (3) late apoptotic cells: Annexin V (+) and dead cell marker (+); and (4) cells that died via non-apoptotic pathways: Annexin V (−) and dead cell marker (+).Cells were counted using the Muse™ Cell Analyzer (Merck Millipore) and analyzed using software provided by Merck Millipore.

Cell Cycle Assay
In total, 1 × 10 6 cells were counted for the Muse™ Cell Cycle Assay that consisted in the use of the nuclear DNA intercalating stain RNAse A and propidium iodide (PI) in a proprietary formulation.The latter was used to discriminate cells in different phases of the cell cycle, based on differential DNA content in the presence of RNAse to increase the specificity of DNA staining.After treatment with the active G. cydonium fraction, cells were washed with phosphate buffered saline (PBS) and centrifuged.The supernatant was removed and 1 mL of ice cold 70% ethanol was added to the re-suspended cell pellet.Samples were capped and frozen at −20 • C for at least 3 h prior to staining.Ethanol-fixed cells were washed with PBS and incubated with 200 µL of Muse™ Cell Cycle Reagent for 30 min at room temperature, in the dark.After staining, cells were processed for cell cycle analysis.

RNA Extraction and Real Time qPCR (RT-qPCR)
RNA isolation and cDNA preparation were performed as previously described [10].The reverse-transcribed products were used to perform a qPCR in order to evaluate the expression level of transcripts of selected genes.Sequences for mRNAs from the nucleotide data bank (National Center for Biotechnology Information, Bethesda, MD, USA) were used to design primer pairs for RT-qPCR (Primer Express, Applied Biosystems, Foster City, CA, USA).Oligonucleotides were obtained from Sigma-Aldrich.The primer sequences are provided in Table 2.
Starting with 2 µg of total RNA, we prepared 20-fold dilution of the resulting cDNA to achieve the concentration equivalent to 100 ng of RNA (Life Technologies-Invitrogen), according to the manufacturer's instructions.Ten nanograms of cDNA were amplified in a total volume of 25 µL containing 1X SYBR Green PCR Master Mix (Applied Biosystems) and 300 nM of forward and reverse primers.The thermal profile conditions were as follows: 5 min of denaturation at 95 • C followed by 44 cycles at 95 • C for 30 s and 60 • C for 1 min.We have added one cycle for melting curve analysis at 95 • C for 15 s, 60 • C for 15 s and 95 • C for 15 s to verify the presence of a single product.Melting-curve analysis was carried out after amplification to verify the validity of the amplicon.Each assay included a no-template control for each primer pair.To capture intra-assay variability, all RT-qPCR reactions were carried out in triplicate.For all RT-qPCR experiments, the data from each cDNA sample were normalized using β-actin mRNA as endogenous level [10].Sample ∆C t values were calculated as the difference between the means of gene markers C t and housekeeping assay C t from the same sample.The 2x-fold expression level was chosen as the threshold for significance of target genes.Statistical analyses (paired Student's t) were performed using Prism software (Graphpad Software, La Jolla, CA, USA).

Extraction of the Polar Fractions in Untreated and Treated Cancer Cells
All cancer cell lines were plated in cell culture flasks (~2 × 10 6 cells/flask) and treated with the active fraction at the IC 50 concentration obtained after 48 h treatment.After incubation time (48 h), cellular supernatants were collected and stored at −80 • C for further investigation.Cell pellets obtained by trypsin digestion were washed twice in Phosphate buffered saline and deuterated water (PBS-D 2 O) and refrigerated at −80 • C. Subsequently they were re-suspended in 170 µL of H 2 O and 700 µL of methanol and were sonicated for 30 s.Then, 350 µL of chloroform was added and cell pellets were shaken on ice in an orbital shaker for 10 min.H 2 O/chloroform (350 µL, 1:1, v/v) was added to each cell suspension and centrifuged at 10,000 rpm for 10 min at 4 • C. Thereafter, the aqueous (polar) and lipophilic (apolar) phases were collected separately and evaporated by SpeedVac system.

1 H-NMR Metabolomic Analysis of the Cellular Polar Fractions
A 600-MHz Bruker Avance DRX spectrometer with a cryoprobe was used to acquire 1 H spectra on the cellular polar fractions.They were dissolved in 630 µL of PBS-D 2 O with the pH adjusted to 7.20, and 70 µL of sodium salt of 3-(trimethylsilyl)-1-propanesulfonic acid (1% in D 2 O) used as the internal standard.
All 1 H-NMR spectra were acquired at 300 K with the excitation sculpting pulse sequence to suppress water resonance.A double-pulsed field gradient echo was used, with a soft square pulse of 4 ms at the water resonance frequency and with gradient pulses of 1 ms duration, adding 128 transients of 64,000 complex points, with an acquisition time of 4 s/transient.Time domain data were all zero-filled to 256,000 complex points and an exponential amplification of 0.6 Hz was applied prior to Fourier transformation.

Statistical and Pathway Analysis
The spectral 0.50-8.60ppm region of 1 H-NMR spectra was integrated in buckets of 0.04 ppm by the AMIX package (Bruker, Biospin GmbH, Rheinstetten, Germany).The water resonance region (4.5-5.2 ppm) was excluded during the analysis and the bucketed region was normalized to the total spectrum area using Pareto scaling.Orthogonal Projections to Latent Structures discriminant analysis (OPLS-DA) was used to compare the spectra obtained on the polar phases from three breast cancer cell lines before and after treatment because OPLS-DA can more effectively cope with chemical shift variation in full-resolution 1 H-NMR datasets [37] without requiring binning or alignment steps.Pathway analysis on the metabolites that were modulated after treatment was performed using Metabo Analyst tool [38].

Bio-Plex Assay
Several cytokines, chemokines, and growth factor levels were evaluated concurrently with the Bio-Plex assay that containing dyed microspheres conjugated with a monoclonal antibody highly specific for a target protein.The method was carried out according to the manufacturer's instructions (Bio-Plex Bio-Rad) to assess the cytokines levels.The Bio-Plex Pro Human Cytokine 27-Plex Immunoassay has been used on supernatants of the three lines of human breast cancer after treatment with sponge extract concentrations.This panel consists of: IL-1β, IL-1ra, IL-2, IL-4, IL-5, IL-6, IL-7, CCL2, CCL11, CXCL10, CXCL8, IFN-γ, IL-9, IL-10, IL-12 (p70), IL-13, IL-15, IL-17, basic FGF, G-CSF, GM-CSF, MIP-1α, MIP-1β, PDGF-ββ, RANTES, TNF-α, and VEGF.Protein levels were determined using a Bio-Plex array reader (Luminex, Austin, TX, USA) that quantifies multiplex immunoassays in a 96-well format with very small fluid volumes.The analyte level was calculated using a standard curve, with software provided by the manufacturer (Bio-Plex Manager Software, Austin, TX, USA).A Bio-Plex array reader (Luminex, Austin, TX, USA) that quantifies multiplex immunoassays in a 96-well format with very small fluid volumes, has been used for protein level determination.The levels of the analytes were calculated using a standard curve, with the Bio-Plex Manager Software provided by the manufacturer.

Figure 2 .
Figure 2. (Left) Two-dimensional molecular network of the active fraction from the sponge Geodia cydonium.In the clusters indicated with (a-f), nodes are labeled with parent m/z ratio (M+H) + ions; edge thickness is related to cosine similarity score; (Right) Liquid chromatography-high resolution mass spectrometry (LC-HRMS) profile of the active fraction; the most abundant ions are represented by colored bold contour squares.

Figure 2 .
Figure 2. (Left) Two-dimensional molecular network of the active fraction from the sponge Geodia cydonium.In the clusters indicated with (a-f), nodes are labeled with parent m/z ratio (M+H) + ions; edge thickness is related to cosine similarity score; (Right) Liquid chromatography-high resolution mass spectrometry (LC-HRMS) profile of the active fraction; the most abundant ions are represented by colored bold contour squares.

Figure 3 .
Figure 3. RT-qPCR analysis: mRNA fold changes were evaluated as ratios between the expression levels of five genes in three breast cancer cell lines, MCF-7, MDA-MB231 and MDA-MB468, after treatment with the active fraction compared to untreated cells.

Figure 3 .
Figure 3. RT-qPCR analysis: mRNA fold changes were evaluated as ratios between the expression levels of five genes in three breast cancer cell lines, MCF-7, MDA-MB231 and MDA-MB468, after treatment with the active fraction compared to untreated cells.

Figure 4 .
Figure 4. Orthogonal Projections to Latent Structures discriminant analysis (OPLS-DA) plots: (A) OPLS-DA and variable importance in projection (VIP) analysis where the metabolites increased or decreased in the endo-metabolome of: (B) MCF-7; (C) MDA-MB231; and (D) MDA-MB468 cells after treatment with the active fraction from the sponge Geodia cydonium, compared to untreated cells.

Figure 4 .
Figure 4. Orthogonal Projections to Latent Structures discriminant analysis (OPLS-DA) plots: (A) OPLS-DA and variable importance in projection (VIP) analysis where the metabolites increased or decreased in the endo-metabolome of: (B) MCF-7; (C) MDA-MB231; and (D) MDA-MB468 cells after treatment with the active fraction from the sponge Geodia cydonium, compared to untreated cells.

Figure 5 .
Figure 5. Cytokine levels in breast cancer cells.Scheme reporting the cytokines modulated by the active fraction obtained from the sponge Geodia cydonium.In particular, anti-inflammatory cytokines that are increased after treatment are reported in red and pro-inflammatory and pro-angiogenic cytokines that are decreased after treatment are reported in green.Blank cells indicate that the cytokines are not modulated by active fraction.

Figure 5 .
Figure 5. Cytokine levels in breast cancer cells.Scheme reporting the cytokines modulated by the active fraction obtained from the sponge Geodia cydonium.In particular, anti-inflammatory cytokines that are increased after treatment are reported in red and pro-inflammatory and pro-angiogenic cytokines that are decreased after treatment are reported in green.Blank cells indicate that the cytokines are not modulated by active fraction.

Table 1 .
Apoptosis.Percentage of live, apoptotic and dead cells expressed as mean ± standard deviation by the Muse Annexin V and Dead Cell assay in MCF-7, MDA-MB231 and MDA-MB468 cells after treatment with IC 50 concentrations obtained after 48 h.Untreated cells were used as control.

Table 2 .
Primer sequences of the genes used in this study.

Table 4 .
List of 1 H chemical shift (ppm) of metabolites found in three breast cancer cells.