Next Article in Journal
Nile Tilapia Derived TP4 Shows Broad Cytotoxicity toward to Non-Small-Cell Lung Cancer Cells
Previous Article in Journal
Pro-Apoptotic Activity of French Polynesian Padina pavonica Extract on Human Osteosarcoma Cells
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Marine Drugs
Volume 16
Issue 12
10.3390/md16120505
marinedrugs-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

MS/MS-Based Molecular Networking Approach for the Detection of Aplysiatoxin-Related Compounds in Environmental Marine Cyanobacteria

by
Chi Ying Gary Ding
1,†,
Li Mei Pang
2,†,
Zhao-Xun Liang
2,
Kau Kiat Kelvin Goh
3,
Evgenia Glukhov
4,
William H. Gerwick
4 and
Lik Tong Tan
1,*
1
Natural Sciences and Science Education, National Institute of Education, Nanyang Technological University, Singapore 637616, Singapore
2
School of Biological Sciences, Nanyang Technological University, Singapore 637551, Singapore
3
Singapore Phenome Centre, Lee Kong Chian School of Medicine, Nanyang Technological University, Singapore 636921, Singapore
4
Center for Marine Biotechnology and Biomedicine, Scripps Institution of Oceanography, University of California at San Diego, La Jolla, CA 92093, USA
*
Author to whom correspondence should be addressed.
†
These two authors contributed equally to this article.
Mar. Drugs 2018, 16(12), 505; https://0-doi-org.brum.beds.ac.uk/10.3390/md16120505
Submission received: 12 November 2018 / Revised: 5 December 2018 / Accepted: 10 December 2018 / Published: 13 December 2018
(This article belongs to the Special Issue “Harmful Algae: Natural Products Research and Medical Application Exploration”―Selected Papers from the 18th International Conference on Harmful Algae (ICHA 2018))
Browse Figures
Versions Notes

Abstract

:
Certain strains of cyanobacteria produce a wide array of cyanotoxins, such as microcystins, lyngbyatoxins and aplysiatoxins, that are associated with public health issues. In this pilot study, an approach combining LC-MS/MS and molecular networking was employed as a rapid analytical method to detect aplysiatoxins present in four environmental marine cyanobacterial samples collected from intertidal areas in Singapore. Based on 16S-ITS rRNA gene sequences, these filamentous cyanobacterial samples collected from Pulau Hantu were determined as Trichodesmium erythraeum, Oscillatoria sp. PAB-2 and Okeania sp. PNG05-4. Organic extracts were prepared and analyzed on LC-HRMS/MS and Global Natural Product Social Molecular Networking (GNPS) for the presence of aplysiatoxin-related molecules. From the molecular networking, six known compounds, debromoaplysiatoxin (1), anhydrodebromoaplysiatoxin (2), 3-methoxydebromoaplysiatoxin (3), aplysiatoxin (4), oscillatoxin A (5) and 31-noroscillatoxin B (6), as well as potential new analogues, were detected in these samples. In addition, differences and similarities in molecular networking clusters related to the aplysiatoxin molecular family were observed in extracts of Trichodesmium erythraeum collected from two different locations and from different cyanobacterial species found at Pulau Hantu, respectively.

1. Introduction

Marine filamentous cyanobacteria are prolific sources of novel secondary metabolites, and a number of these compounds possess significant biological activity, having agricultural and pharmaceutical applications [1,2]. A notable example is a series of highly cytotoxic lead compounds belonging to the dolastatin class of molecules [3]. For instance, dolastatin 10 and synthetic derivatives, soblidotin (TZT-1027) and tasidotin, entered clinical trials in the 1990s as anticancer agents but dropped out in Phase II due to toxicity and/or lack of efficacy in trials [4]. A synthetic analogue, monomethyl auristatin E, was subsequently formulated as an antibody-drug conjugate, brentuximab vendotin (Adecetris®), for the treatment of Hodgkin’s lymphoma and systemic anaplastic large cell lymphoma in 2012 [5].
In recent years, various studies revealed that certain marine cyanobacterial strains produce compounds that concurrently provide beneficial and detrimental effects to human health. One such compound is the polyketide-based cyanotoxin, aplysiatoxin. Cyanotoxins collectively cause a wide spectrum of health conditions, ranging from dermatitis to severe carcinoma [6,7,8]. Aplysiatoxin and its related compounds, which includes oscillatoxin A and nhatrangins, have been reported from Lyngbya majuscula, Oscillatoria nigro-viridis and Trichodesmium erythraeum [6]. The tumor promoting activity of aplysiatoxin and debromoaplysiatoxin is attributed to their specific molecular roles as protein kinase C activators. Interestingly, a simplified synthetic analogue of aplysiatoxin, aplog-1, was recently reported to possess anti-cancer properties with mode of action similar to bryostatin 1 [9,10,11,12,13]. In addition, studies revealed several aplysiatoxin-related compounds to display therapeutic activities, such as anti-Chikungunya and modulation of potassium ion channels [14,15]. A recent study by Richard et al. demonstrated that aplysiatoxin and its debrominated analogue induced expression of latent HIV-1 provirus in both cell line and primary cell models [16].
Due to the varied biological activities of the aplysiatoxins, they make good candidates for biomedical research as well as in their detection and monitoring of toxins during algal bloom occurrence. Given the correlation between global warming and occurrence of cyanobacterial blooms, the incidence of illness caused by contact with cyanotoxins can potentially be increased [8]. In contrast to research on hepatotoxins (e.g., microcystins and nodularin) as well as other neurotoxins (e.g., β-N-methylamino-l-alanine), very little work has been conducted on the development of detection methods of other cyanotoxins, such as the aplysiatoxins [17]. In this report, we employed MS/MS-based molecular networking approach for the rapid detection of aplysiatoxin class of molecules from environmental marine cyanobacterial samples obtained from four intertidal areas in Singapore. The detection method, based on either HPLC-HRMS/MS or UPLC-HRMS/MS combined with data analysis using Global Natural Product Social Molecular Networking (GNPS), led to the detection of both known and new analogues related to the aplysiatoxin class of cyanotoxins.

2. Results

2.1. Phylogenetic Analysis of Environmental Marine Cyanobacterial Samples

Of the four marine cyanobacterial samples analyzed in this study, three samples were collected between June and August 2016 at three intertidal areas located at Pulau Hantu Besar (TLT/PHB/001), Pulau Hantu Central Lagoon area 1 (TLT/PHC/001) and Pulau Hantu Central Lagoon area 2 (TLT/PHC/002). A fourth cyanobacterial sample, Trichodesmium erythraeum, was previously collected from Pulau Seringat Kias (TLT/PSK/001) in August 2012 [8]. A phylogenetic tree comprising of the three cyanobacterial samples from Pulau Hantu was constructed against seven cyanobacterial strains and is presented in Figure 1. Gene sequence data, obtained after amplification of the 16S rRNA gene, for TLT/PHB/001, TLT/PHC/001 and TLT/PHC/002 was assigned to Trichodesmium erythraeum with a 97.8% sequence identity (1465 bp, 100% gene sequence coverage), Okeania sp. PNG05-4 with a 98.8% sequence identity (1399 bp, 95.8% gene sequence coverage), and Oscillatoria sp. PAB-21 with a 99.1% sequence identity (1363 bp, 93.4% gene sequence coverage), respectively. Individual score and coverage/identity giving corresponding closest homologs based on 16S-ITS rRNA sequences from TLT/PHB, TLT/PHC001 and TLT/PHC002 were tabulated in Table S1 (Supplementary Data), and complete 16S rRNA sequence in Table S2 (Supplementary Data).

2.2. MS/MS-Based Molecular Networks of Marine Cyanobacterial Extracts

Organic crude extracts from each of the cyanobacterial samples were prepared and analyzed initially by UPLC–HRMS/MS Waters system (Milford, MA, USA). In addition, organic crude extract from TLT/PSK/001 were analyzed on both HPLC-HRMS/MS Thermo–Finnigan system (Waltham, MA, USA) (TLT/PSK/001a) and UPLC-HRMS/MS Waters system (TLT/PSK/001b). Mass spectrometric data obtained from these analyses were then used to generate molecular networking clusters using the online platform Global Natural Products Social Molecular Networking (GNPS). Debromoaplysiatoxin (1) (Figure 2), isolated from the extract of TLT/PHC/001, was used as seed molecule in the analysis (HRESIMS of sodiated molecular ion m/z 615.3005, calculated for C32H48O10Na+ 615.2987, difference = 2.9 ppm). Molecular networks of the extracts (an example is shown for TLT/PHB/001 in Figure 3), visualized using the software Cytoscape v3.5.1 (The Cytoscape Consortium, New York, NY, USA), revealed the presence of clusters related to the molecular family of the aplysiatoxin class in all four environmental cyanobacterial samples (clusters A to E in Figure 4). In addition to the presence of the seeded debromoaplysiatoxin, five other related compounds (Figure 2), namely, anhydrodebromoaplysiatoxin (2), 3-methoxydebromoaplysiatoxin (3), aplysiatoxin (4), oscillatoxin A (5) and 31-noroscillatoxin B (6), were detected in clusters A, B, C and E based on their HRMS (High-Resolution Mass Spectrometry) data (Table 1 and Figure 4). Furthermore, potential new aplysiatoxin analogues, based on their relatedness to 1, were detected from the clusters. These potential new compounds are indicated as 15 nodes from cluster A, eight nodes from both cluster B and C, seven nodes from cluster D and 17 nodes from cluster E (Figure 4). Some putative assignments for the new aplysiatoxin analogues were included, based on the HRMS m/z values.

3. Discussion

Certain cyanobacteria strains are known to produce a wide array of toxic metabolites or cyanotoxins, such as aplysiatoxins, microcystins and nodularins [7]. Many of these cyanotoxins affect human health through their actions as tumor promoters, diarrhetic, and pro-inflammatory activities [18]. In addition, these metabolites may also limit various aquatic recreational activities, thereby affecting related economies. Such impacts are particularly acute in certain geographical regions, such as Singapore, where its location in the tropics provides a conducive environment for the occurrence, growth, and more importantly blooms of marine cyanobacteria [19]. For instance, regular occurrence of marine cyanobacterial blooms has been sighted at several offshore islands in Singapore, including Pulau Semakau, P. Seringat-Kias, P. Hantu, P. Ubin and P. Sakijang Bendera. The cyanobacterial samples used in this study were obtained from islands designated as tourist destinations with aquatic recreational activities, including SCUBA diving and swimming, as the main attraction. The persistent occurrence of marine cyanobacteria at these islands (personal observation by Lik Tong Tan) warrants setting up of a rapid and relatively straightforward analytical method for the detection of cyanotoxins from environmental samples. A rapid detection method is integral and equally imperative to prevent the general public’s exposure to these toxic compounds.
In this study, molecular networking clusters of various environmental marine cyanobacterial extracts were analyzed using a MS/MS-based molecular networking approach [20]. Initial chemical investigation of the organic extract obtained from the cyanobacterial sample TLT/PHC/001 collected at Pulau Hantu resulted in the isolation of the known cyanotoxin debromoaplysiatoxin (1). This molecule was subsequently used as a seed molecule for the creation of the aplysiatoxin molecular family in the molecular network. Each node within the clusters represent an ionized compound and dereplication was performed by comparison with the m/z values of known aplysiatoxin-related analogues (Table 2).
In addition to debromoaplysiatoxin, a number of known analogues, anhydrodebromoaplysiatoxin (2), 3-methoxydebromoaplysiatoxin (3), aplysiatoxin (4), oscillatoxin A (5) and 31-noroscillatoxin B (6), were found to be present in the molecular network clusters A, B, C and E. The identity of these known compounds was confirmed by comparing the calculated m/z values with that obtained from HRMS/MS (differences of less than 5 ppm) as shown in Table 1 and Table 2. Compounds 1–3 were present in clusters A, B, C and E. Additional compounds 4 and 5 were also present in cluster A while 6 was found in cluster E.
The phylogenetic results of the cyanobacterial samples, coupled with careful analysis of the molecular networking clusters of their extracts, revealed a number of interesting observations. Though both TLT/PHB/001 and TLT/PSK/001 cyanobacterial samples have been identified as Trichodesmium erythraeum, their clusters (A and E), related to the aplysiatoxin molecular family based on the seeded 1, were different (Figure 4). Clusters A and E gave 15 and 17 nodes, respectively, both revealing the presence of all three compounds 1–3. The difference in the complexity of these two Trichodesmium erythraeum-derived molecular network clusters and the number of nodes present may be attributed to the versatility of the biosynthetic gene clusters present in the same cyanobacteria, and different environmental triggers present in the two locations where the cyanobacteria were collected. On the other hand, aplysiatoxin-related clusters B and C produced from two different cyanobacterial species, TLT/PHC/001 (identified as Okeania sp. PNG05-4) and TLT/PHC/002 (identified as Oscillatoria sp. PAB-21), respectively, with only slight difference in the m/z (within 5 ppm) observed between corresponding nodes. Interestingly, these two samples were collected from sites that were just 100 m apart at Pulau Hantu. In addition, when comparing the full molecular networks of these two different cyanobacterial species (Figure 5), minimal differences in the compounds, aside from the aplysiatoxin-related cluster, were observed. These observations suggest that similar compounds are largely produced as a result of the same environmental triggers, though by different cyanobacteria, and this can be validated with the help of results obtained through the investigations of related biosynthetic gene clusters of TLT/PHC/001 and TLT/PHC/002 in producing the aplysiatoxin-related compounds.
It has been shown that the MS/MS-based molecular networking is useful for dereplication of known compounds as well as the discovery of potentially new natural products [21]. In molecular networking, MS/MS spectra are compared based on the similarity and differences in fragmentation patterns. Structurally related compounds can then be clustered as nodes due to certain degrees of similarity between the compounds expressed as edges in the molecular network. In addition, the availability of HRMS/MS spectral data of known cyanotoxins in GNPS database could facilitate the detection process [22]. Recently, Teta and co-workers used a similar MS-based molecular networking approach to detect new microcystin variants from the 2014 freshwater cyanobacterial harmful bloom at Green Lake, Seattle [23]. In another study by Briand et al., molecular networking was employed to track changes in secondary metabolic profiles, including microcystins and other peptides, of Microcystis aeruginosa strains due to intraspecific interactions [24]. In this present study, two nodes with m/z 629.219 and 629.220 observed in clusters B, C and E are likely to be compound 3. However, it is evident that only the node with observed value of 629.22 fits within the accepted ppm difference for compound 3 (Figure 4). Considering the presence of both nodes found within the same cluster, the node of m/z 629.219 could be a new analogue related to 2. A similar relationship was observed from the nodes of m/z 615.218 and 615.298 present in the same set of clusters (Figure 4). Some putative assignments were provided to unknown compounds based on observed HRMS of the nodes within each cluster. An example is the compound with m/z value of 631.247 in clusters B and C that could correspond to an analogue of 3 having a reduced carbonyl moiety or cleavage to the oxide functionality of either tetrahydropyran rings.
The advantages of use and advancements in instrumentation provided a variety of MS/MS systems used by laboratories today. An interesting technical approach adopted in this study was to employ two different systems in the analysis, namely a low-resolution Thermo-Finnigan HPLC-MS/MS and a high-resolution Waters UPLC-MS/MS. The main reason was to make comparisons of the molecular networks created based on each system. The organic extract from TLT/PSK/001 was randomly chosen and molecular networks were then constructed individually using the MS/MS data acquired from both MS systems using similar parameters. Analysis of the molecular network clusters TLT/PSK/001a (D) and TLT/PSK/001b (E) using Thermo–Finnigan and Waters systems, respectively, revealed several differences. The first difference was the size of each molecular network due to a difference in quantity of nodes present, with 7 and 17 nodes in clusters D and E, respectively. Another difference was observed comparing compounds 1, 2 and 6 present in both clusters, where each corresponding node in both clusters gave different m/z values. A case in point is the node that corresponds to 2, with m/z of 597.267 and 597.261 observed in clusters D and E, respectively. The chemical identity of the corresponding nodes was achieved and confirmed to be 2 in both clusters, where a mass difference of within 5 ppm was determined based on the observed and calculated m/z obtained from its respective system. The last difference was the detection of different compounds from the same cyanobacterial strain in the two molecular networks. Compounds 4 and 5 were observed exclusively in cluster D obtained from the Thermo–Finnigan system; compound 3 observed in cluster E created using the Waters system. While the exact explanation could not be provided, these differences between the clusters may be attributed largely to the different MS instruments’ operating parameters (i.e., resolution, type of ionization and detector), hence care should be taken while comparing molecular networks from different systems.
Current aquatic biomonitoring and ecotoxicology programs in Singapore are confined to the analysis of freshwater cyanobacterial strains for hepatoxins, such as microcystins and nodularin. Methods, such as mass spectrometry, capillary electrophoresis and qPCR, have been used for the detection of microcystins as well as for assessing cyanotoxin risk in freshwater cyanobacterial strains obtained from reservoirs in Singapore [25,26,27]. The unprecedented nature of this pilot study underscores the importance of the MS/MS-based molecular networking approach as a rapid analytical technique to detect non-peptidic cyanotoxins, such as the aplysiatoxins, from environmental samples. Moreover, recent reports showed that several aplysiatoxin analogues possess therapeutic activities, such as antiviral and potassium channel inhibition activities [14,15]. In light of these discoveries, the MS/MS-based molecular networking approach employed in this study could aid in the detection of new aplysiatoxin analogues that harbors varied biomedical potential, contributing to the growing list of marine natural products as therapeutic agents.

4. Materials and Methods

4.1. General Experimental Procedures

Proton NMR spectra were recorded in CDCl3 on a 400 MHz Bruker (Billerica, MA, USA) NMR spectrometer using the residual solvent signals (δH at 7.28 ppm) as internal standards. High-performance liquid chromatography isolation and purification of compound 1 was conducted on a Shimadzu (Nakagyo-ku, Kyoto, Japan) LC-8A preparative LC coupled to a Shimadzu SPD-M10A VP diode array detector. The UPLC-HRMS/MS analysis was conducted on a Waters (Milford, MA, USA) Xevo G2-XS qTOF with an ESI positive ion mode and data-dependent acquisition mode to obtain MS and MS/MS data. The HPLC-MS/MS analysis was conducted on a Thermo–Finnigan (San Jose, CA, USA) LCQ Advantage ion trap mass spectrometer with an ESI positive ion mode and data-dependent manner to obtain MS and MS/MS data. High-performance liquid chromatography preparative and analytical columns was obtained from Phenomenex (Torrance, CA, USA), and a UPLC analytical column was obtained from Waters. All solvents were HPLC-grade purchased from Thermo-Fisher Scientific (Waltham, MA, USA) and HPLC-grade water obtained by filtration using a Milli-Q Direct water purification system from Millipore (Billerica, MA, USA) unless otherwise stated. Deuterated NMR solvents were purchased from Cambridge Isotope Laboratories (Tewksbury, MA, USA).

4.2. Marine Cyanobacterial Collections

Samples of the benthic filamentous marine cyanobacteria were collected from offshore islands of Pulau Seringat Kias on August 7, 2012 (TLT/PSK/001), Pulau Hantu Besar, on 9 June 2016 (TLT/PHB/001), and Pulau Hantu Central Lagoon area 1 on 6 August 2016 (TLT/PHC/001) and Pulau Hantu Central Lagoon area 2 on 9 June 2016 (TLT/PHC/002). Voucher samples are stored at Natural Products Chemistry Laboratory, National Institute of Education, Nanyang Technological University, Singapore as TLT/PSK/07AUG2012/001, TLT/PHB/09JUNE2016/001, TLT/PHC/06AUG2016/001 and TLT/PHC/09JUNE2016/002, respectively. A small amount of the cyanobacterial samples collected in 2016 (TLT/PHB/001, TLT/PHC/001 and TLT/PHC/002) was preserved for genetic analysis in RNAlater® (Ambion, Austin, TX, USA). Chemistry samples were preserved in 70% aqueous EtOH at −20 °C before chemical extraction.

4.3. DNA Isolation and PCR Amplification

Genetically preserved biomass (~ 150 mg) of each of the three filamentous marine cyanobacteria samples were washed 3 times with 3M NaCl solution. Genomic DNA was extracted using CTAB procedure as described previously, with increased lysozyme concentration of 5 mg/mL. The 16S rRNA genes and the 16S-23S ITS regions were PCR-amplified using primers 27F (5′-AGAGTTTGATCMTGGCTCAG-3′) and 340 (5′-CTCTGTGTGCCTAGGTATCC-3′) [28,29]. The PCR reaction volumes were 25 µL containing: 1 µL (~100 ng) of DNA, 5 µL of 5 × Q5 reaction buffer, 0.5 µL of (10 mM) dNTP mix, 1.25 µL of each primer (10 µM), 0.25 µL of Q5 High-Fidelity DNA polymerase and 15.75 µL of dH2O. The PCR reactions were performed in a T100TM Thermal Cycler (Bio-Rad, Hercules, CA, USA) using the following program: Initial denaturation for 30 s at 98 °C, amplification by 30 cycles of 10 s at 98 °C, 15 s at 63 °C and 2 min at 72 °C, and final elongation for 10 min at 72 °C. The PCR products were analyzed on 1% agarose gel in TAE (Tris-acetate-EDTA) buffer and visualized by gel red. The PCR products were purified using NucleoSpin Gel and PCR Clean-up kit (Macherey-Nagel, Duren, Germany). All purified PCR products of 27F and 340 primers were sequenced with 883R primer (5′-ATTAAACCACATACTCCACC-3′), FR primer (5′-ACGAGCTGACGACAGCCATG-3′), HR primer (5′-AAGGAGGTGATCCAGCCGCA-3′) and 1052F (5′-CATCATTAAGTTGGGCACTC-3′) [22]. Sequencing was performed at 1stBASE Laboratories Sdn BhD (Singapore).

4.4. Phylogenetic Analysis

Gene sequences were aligned using MUSCLE (European Bioinformatics Institute, EMBL-EBI, Cambridge, United Kingdom) [30]. The phylogenetic tree was plotted using maximum likelihood (ML) in MEGA, v6 (Pennsylvania State University, State College, PA, USA) [31]. The 16S rRNA gene sequence of each cyanobacterial sample was multiply aligned with 7 cyanobacterial strains obtained from the GenBank/EMBL/DDBJ databases. Bootstrap re-sampling was performed with 1000 replicates.

4.5. Extraction of Crude Extracts and Isolation of Debromoaplysiatoxin (1)

Each of the four environmental cyanobacterial samples (about 10 g each) were exhaustively extracted with 2:1 CH2Cl2/MeOH. The crude extract from TLT/PHC/001 (1.421 g) was chosen and dried under vacuum and applied to a silica gel VLC column. The column was eluted to produce nine sub-fractions: C1 (elution by 100% hexanes), C2 (elution by 10% EtOAc/hexanes), C3 (elution by 30% EtOAc/hexanes), C4 (elution by 50% EtOAc/hexanes), C5 (elution by 60% EtOAc/hexanes), C6 (elution by 80% EtOAC/hexanes), C7 (elution by 100% EtOAc), C8 (elution by 10% MeOH/EtOAc) and C9 (20% MeOH/EtOAc). Sub-fractions were subjected to Brine Shrimp Toxicity Assay and sub-fraction C6 displayed mortality at 93.2% and 81.9% when tested at 100 ppm and 10 ppm, respectively. The remaining sub-fraction C6 (121.4 mg) was filtered over a C18 SPE cartridge by elution with 50% MeOH/H2O. The filtered sample (29.0 mg) was subjected to preparative reversed-phase HPLC separation using a Phenomenex Luna phenyl-hexyl column (250 mm × 100 mm, 5 μm). Elution at 3.0 mL/min using eluent A: H2O; eluent B: CH3CN; gradient: 38% B for 1.5 min, 38→80% B in 5.5 min, 80→100% B in 5 min, then held at 100% B for 18 min to yield 1 (8 mg, tR = 20.024 min). Structure and identity of 1 was confirmed based on proton NMR in CDCl3 (Figure S1, Supplementary Data) and HRESIMS m/z 615.3005 (calculated for C32H48O10Na+, 615.2987, difference = 2.9 ppm). Mass spectrometry parameters included a spray voltage of 3 kV, cone voltage of 40 V, source temperature of 120 °C, desolvation temperature of 450 °C, cone gas flow of 50 L/h and desolvation gas flow of 950 L/h.

4.6. UPLC-HRMS/MS Analysis

Crude extracts were re-dissolved in 1 mL CH3CN with vortex mixing over 5 min before transferring into separate Eppendorf tubes. Tubes were then centrifuged at 10,000 rpm at 4 °C over 10 min and the supernatant was aliquoted and finally diluted with CH3CN to 10,000× dilution. Then, 1.5 μL of each 10,000× dilution sample was subjected to UPLC-HRMS/MS performed using the same high-resolution spectrometer and parameters stated above, together with a Waters ACQUITY BEH C1850 mm × 2.1 mm, 1.7 μm column and maintained at a column temperature of 40 °C and sample temperature at 4 °C using a step elution program of ‘1’ based on the Waters system. Step elution program was as follows: Mobile phase of 98% CH3CN in 0.1% aq. HCOOH/0.1% aq. HCOOH, run time of 14 min and flow rate of 0.5 mL/min. All the mass spectra were recorded in the positive-ion mode. Data were collected in the data-dependent acquisition mode, in which the first ten most intense ions of a full-scan mass spectrum were subjected to MS/MS analysis. Tandem mass spectrometry scans were obtained for selected ions with a mass range of 100–2000 Da, MS scan time of 0.33 s over 12 min, MS/MS scan time of 0.10 s and a collision energy ramp of 10–50 V.

4.7. HPLC-MS/MS Analysis

TLT/PSK/001 was re-dissolved in MeOH and 20 μL of the solution was subjected (no-waste mode) to HPLC-MS/MS using a Thermos–Finnigan system with a Phenomenex Kinetex C 18 100 mm × 4.6 mm, 5 μm column. Elution at 0.5 mL/min using eluent A: 0.1% aq. HCOOH in H2O; eluent B: CH3CN; gradient: 30→100% B over 20 min, then held at 100% B for 10 min. The HPLC eluate was electrospray ionized at 35 eV and analyzed for positive ions, and MS/MS spectra were obtained in a data-dependent manner using collision induced dissociation (CID) at 35 eV.

4.8. Molecular Networking Creation

The mass spectral data were converted to .mzXml format using the program msconvert (ProteoWizard Software Foundation, Los Angeles, CA, USA). A molecular network was created using the online workflow at GNPS website [32]. The data were filtered by removing all MS/MS peaks within ± 17 Da of the precursor m/z. Tandem mass spectrometry spectra were window filtered by choosing only the top 6 peaks in the ± 50 Da window throughout the spectrum. The data were 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.65 and more than 10 matched peak. 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. The spectra in the network were then searched against GNPS’ spectral libraries. The library spectra were filtered in the same manner as the input data. The data were then imported into Cytoscape v3.5.1 (The Cytoscape Consortium, New York, NY, USA) and displayed as a network of nodes and edges [26]. All matches kept between network spectra and library spectra were required to have a score above 0.65 and at least 1 matched peak. Analogue search was enabled against the library with a maximum mass shift of 100.0 Da.

Supplementary Materials

The following are available online at https://0-www-mdpi-com.brum.beds.ac.uk/1660-3397/16/12/505/s1, 1H NMR of 1, and phylogenetic data of TLT/PHB/001, TLT/PHC/001 and TLT/PHC/002.

Author Contributions

C.Y.G.D. performed the extraction of the cyanobacterial samples, the isolation, purification and structure elucidation of 1, and constructed the molecular networks. L.M.P. performed the phylogenetic analysis of the cyanobacterial samples. C.Y.G.D. and K.K.K.G. performed the UPLC-MS/MS analysis on Waters system. E.G. performed the HPLC-MS/MS analysis on Thermos–Finnigan system. C.Y.G.D. and L.M.P. wrote the paper. Z.-X.L., W.H.G. and L.T.T. edited and vetted the paper. All authors read and corrected the paper.

Funding

This research was funded by the National Research Foundation, Prime Minister’s Office, Singapore under its Marine Science Research and Development Programme (MSRDP-P15) awarded to L.T.T.

Acknowledgments

This work was supported by the National Research Foundation, Prime Minister’s Office, Singapore under its Marine Science Research and Development Programme (Award No. MSRDP-P15). The authors would also like to acknowledge Angela Ah Buan Lim-Ang, Anna Say Guek Yee-Lee and Yuen Ming Gary Chan (NSSE, NIE) for providing valuable technical assistance.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Tan, L.T. Bioactive natural products from marine cyanobacteria for drug discovery. Phytochemistry 2007, 7, 954–979. [Google Scholar] [CrossRef]
  2. Abed, R.M.; Dobretsov, S.; Sudesh, K. Applications of cyanobacteria in biotechnology. J. Appl. Microbiol. 2009, 1, 1–12. [Google Scholar] [CrossRef]
  3. Ducat, D.C.; Way, J.C.; Silver, P.A. Engineering cyanobacteria to generate high-value products. Trends Biotechnol. 2011, 2, 95–103. [Google Scholar] [CrossRef]
  4. Pettit, G.R.; Kamano, Y.; Herald, C.L.; Tuinman, A.A.; Boettner, F.E.; Kizu, H.; Schmidt, J.M.; Baczynskyj, L.; Tomer, K.B.; Bontems, R.J. The isolation and structure of a remarkable marine animal antineoplastic constituent: dolastatin 10. J. Am. Chem. Soc. 1987, 22, 6883–6885. [Google Scholar] [CrossRef]
  5. Ju, W.; Zhang, M.; Wilson, K.M.; Petrus, M.N.; Bamford, R.N.; Zhang, X.; Guha, R.; Ferrer, M.; Thomas, C.J.; Waldmann, T.A. Augmented efficacy of Brentuximab Vedotin combined with Ruxolitinib and/or Navitoclax in a murine model of human Hodgkin’s lymphoma. Proc. Natl. Acad. Sci. USA 2016, 6, 1624–1629. [Google Scholar] [CrossRef]
  6. Pearson, L.A.; Dittmann, E.; Mazmouz, R.; Ongley, S.E.; D’Agostino, P.M.; Neilan, B.A. The genetics, biosynthesis and regulation of toxic specialized metabolites of cyanobacteria. Harmful Algae 2016, 54, 98–111. [Google Scholar] [CrossRef]
  7. Zanchett, G.; Oliveira-Filho, E.C. Cyanobacteria and cyanotoxins: from impacts on aquatic ecosystems and human health to anticarcinogenic effects. Toxins 2013, 10, 1896–1917. [Google Scholar] [CrossRef]
  8. Nagarajan, M.; Maruthanayagam, V.; Sundararaman, M. A review of pharmacological and toxicological potentials of marine cyanobacterial metabolites. J. Appl. Toxicol. 2012, 3, 153–185. [Google Scholar] [CrossRef]
  9. Yanagita, R.C.; Kamachi, H.; Tanaka, K.; Murakami, A.; Nakagawa, Y.; Tokuda, H.; Nagai, H.; Irie, K. Role of the phenolic hydroxyl group in the biological activities of simplified analogue of aplysiatoxin with antiproliferative activity. Bioorg. Med. Chem. Lett. 2010, 20, 6064–6066. [Google Scholar] [CrossRef] [Green Version]
  10. Nakagawa, Y.; Kikumori, M.; Yanagita, R.; Murakami, A.; Tokuda, H.; Nagai, H.; Irie, K. Synthesis and biological evaluation of the 12,12-dimethyl derivative of aplog-1, an anti-proliferative analog of tumor-promoting aplysiatoxin. Biosci. Biotechnol. Biochem. 2011, 75, 1167–1173. [Google Scholar] [CrossRef]
  11. Kamachi, H.; Tanaka, K.; Yanagita, R.C.; Murakami, A.; Murakami, K.; Tokuda, H.; Suzuki, N.; Nakagawa, Y.; Irie, K. Structure-activity studies on the side chain of a simplified analog of aplysiatoxin (aplog-1) with anti-proliferative activity. Bioorg. Med. Chem. 2013, 10, 2695–2702. [Google Scholar] [CrossRef]
  12. Kikumori, M.; Yanagita, R.C.; Tokuda, H.; Suzuki, N.; Nagai, H.; Suenaga, K.; Irie, K. Structure-activity studies on the spiroketal moiety of a simplified analogue of debromoaplysiatoxin with antiproliferative activity. J. Med. Chem. 2012, 11, 5614–5626. [Google Scholar] [CrossRef]
  13. Nakagawa, Y.; Yanagita, R.C.; Hamada, N.; Murakami, A.; Takahashi, H.; Saito, N.; Nagai, H.; Irie, K. A simple analogue of tumor-promoting aplysiatoxin is an antineoplastic agent rather than a tumor promoter: Development of a synthetically accessible protein kinase C activator with bryostatin-like activity. J. Am. Chem. Soc. 2009, 22, 7573–7579. [Google Scholar] [CrossRef]
  14. Gupta, D.K.; Kaur, P.; Leong, S.T.; Tan, L.T.; Prinsep, M.R.; Chu, J.J. Anti-Chikungunya viral activities of aplysiatoxin-related compounds from the marine cyanobacterium Trichodesmium erythraeum. Mar. Drugs 2014, 1, 115–127. [Google Scholar] [CrossRef]
  15. Han, B.-N.; Liang, T.-T.; Keen, L.J.; Fan, T.-T.; Zhang, X.-D.; Xu, L.; Zhao, Q.; Wang, S.-P.; Lin, H.-W. Two marine cyanobacterial aplysiatoxin polyketides, neo-debromoaplysiatoxin A and B, with K+ channel inhibition activity. Org. Lett. 2018, 3, 578–581. [Google Scholar] [CrossRef]
  16. Richard, K.; Williams, D.; de Silva, E.; Brockman, M.; Brumme, Z.; Andersen, R.; Tietjen, I. Identification of novel Hiv-1 latency-reversing agents from a library of marine natural products. Viruses 2018, 10, 348. [Google Scholar] [CrossRef]
  17. Kaushik, R.; Balasubramanian, R. Methods and approaches used for detection of cyanotoxins in environmental samples: a review. Crit. Rev. Environ. Sci. Technol. 2013, 13, 1349–1383. [Google Scholar] [CrossRef]
  18. Nagai, H.; Kan, Y.; Fujita, T.; Sakamoto, B.; Hokama, Y. Manauealide C and anhydrodebromoaplysiatoxin, toxic constituents of the Hawaiian red alga, Gracilaria coronopifolia. Biosci. Biotechnol. Biochem. 1998, 5, 1011–1013. [Google Scholar] [CrossRef]
  19. Paerl, H.W.; Huisman, J. Climate change: A catalyst for global expansion of harmful cyanobacterial blooms. Environ. Microbiol. Rep. 2009, 1, 27–37. [Google Scholar] [CrossRef]
  20. Nguyen, D.D.; Wu, C.-H.; Moree, W.J.; Lamsa, A.; Medema, M.H.; Zhao, X.; Gavilan, R.G.; Aparicio, M.; Atencio, L.; Jackson, C.; et al. MS/MS networking guided analysis of molecule and gene cluster families. Proc. Nat. Acad. Sci. U.S.A. 2013, 28, E2611–E2620. [Google Scholar] [CrossRef]
  21. Gerwick, W.H. The face of a molecule. J. Nat. Prod. 2017, 9, 2583–2588. [Google Scholar] [CrossRef]
  22. Bertin, M.J.; Schwartz, S.L.; Lee, J.; Korobeynikov, A.; Dorrestein, P.C.; Gerwick, L.; Gerwick, W.H. Spongosine production by a Vibrio harveyi strain associated with the sponge Tectitethya crypta. J. Nat. Prod. 2015, 3, 493–499. [Google Scholar] [CrossRef]
  23. Teta, R.; Della Sala, G.; Glukhov, E.; Gerwick, L.; Gerwick, W.H.; Mangoni, A.; Costantino, V. Combined LC-MS/MS and molecular networking approach reveals new cyanotoxins from the 2014 cyanobacterial bloom in Green Lake, Seattle. Environ. Sci. Technol. 2015, 24, 14301–14310. [Google Scholar] [CrossRef]
  24. Briand, E.; Bormans, M.; Gugger, M.; Dorrestein, P.C.; Gerwick, W.H. Changes in secondary metabolic profiles of Microcystis aeruginosa strains in response to intraspecific interactions. Environ. Microbiol. 2015, 2, 384–400. [Google Scholar] [CrossRef]
  25. Benke, P.I.; Vinay Kumar, M.C.S.; Pan, D.; Swarup, S. A mass spectrometry-based unique fragment approach for the identification of microcystins. Analyst 2015, 4, 1198–1206. [Google Scholar] [CrossRef]
  26. Te, S.H.; Chen, E.Y.; Gin, K.Y.-H. Comparison of quantitative PCR and droplet digital PCR multiplex assays for two genera of bloom-forming cyanobacteria, Cylindrospermopsis and Microcystis. Appl. Environ. Microl. 2015, 15, 5203–5211. [Google Scholar] [CrossRef]
  27. Birungi, G.; Yau Li, S.F. Determination of cyanobacterial cyclic peptide hepatotoxins in drinking water using CE. Electrophoresis 2009, 15, 2737–2742. [Google Scholar] [CrossRef]
  28. Graca, A.P.; Bondoso, J.; Gaspar, H.; Xavier, J.R.; Monteiro, M.C.; de la Cruz, M.; Oves-Costales, D.; Vicente, F.; Lage, O.M. Antimicrobial activity of heterotrophic bacterial communities from the marine sponge Erylus discophorus (Astrophorida, Geodiidae). PLoS ONE 2013, 11, e78992. [Google Scholar]
  29. Mincer, T.J.; Jensen, P.R.; Kauffman, C.A.; Fenical, W. Widespread and persistent populations of a major new marine actinomycete taxon in ccean sediments. Appl. Environ. Microbiol. 2002, 10, 5005–5011. [Google Scholar] [CrossRef]
  30. Edgar, R.C. MUSCLE: Multiple sequence alignment with high accuracy and high throughput. Nucleic Acids. Res. 2004, 5, 1792–1797. [Google Scholar] [CrossRef]
  31. Tamura, K.; Stecher, G.; Peterson, D.; Filipski, A.; Kumar, S. MEGA 6: Molecular Evolutionary Genetics Analysis version 6.0. Mol. Biol. Evol. 2013, 12, 2725–2729. [Google Scholar] [CrossRef]
  32. Wang, M.; Carver, J.J.; Phelan, V.V.; Sanchez, L.M.; Garg, N.; Peng, Y.; Nguyen, D.D.; Watrous, J.; Kapono, C.A.; Luzzatto-Knaan, T.; et al. Sharing and community curation of mass spectrometry data with Global Natural Products Social Molecular Networking. Nat. Biotechnol. 2016, 34, 828–837. [Google Scholar] [CrossRef] [Green Version]
Marinedrugs 16 00505 g001 550
Figure 1. Phylogenetic tree (maximum likelihood) for three environmental marine cyanobacterial strains, TLT/PHB/001, TLT/PHC/001 and TLT/PHC/002. The scale bar indicates 0.005 expected nucleotide substitutions per site.
Figure 1. Phylogenetic tree (maximum likelihood) for three environmental marine cyanobacterial strains, TLT/PHB/001, TLT/PHC/001 and TLT/PHC/002. The scale bar indicates 0.005 expected nucleotide substitutions per site.
Marinedrugs 16 00505 g001
Marinedrugs 16 00505 g002 550
Figure 2. Chemical structures of aplysiatoxin (4), and related analogues debromoaplysiatoxin (1), anhydrodebromoaplysiatoxin (2), 3-methoxydebromoaplysiatoxin (3), oscillatoxin A (5) and 31-noroscillatoxin B (6).
Figure 2. Chemical structures of aplysiatoxin (4), and related analogues debromoaplysiatoxin (1), anhydrodebromoaplysiatoxin (2), 3-methoxydebromoaplysiatoxin (3), oscillatoxin A (5) and 31-noroscillatoxin B (6).
Marinedrugs 16 00505 g002
Marinedrugs 16 00505 g003 550
Figure 3. Molecular network of the extract obtained from TLT/PHB/001 generated using Cytoscape v3.5.1. (The Cytoscape Consortium, New York, NY, USA).
Figure 3. Molecular network of the extract obtained from TLT/PHB/001 generated using Cytoscape v3.5.1. (The Cytoscape Consortium, New York, NY, USA).
Marinedrugs 16 00505 g003
Marinedrugs 16 00505 g004 550
Figure 4. Cyanobacterial molecular network clusters based on the seeded 1 (nodes outlined) with a cosine similarity cutoff of 0.65 obtained from extracts. (A): TLT/PHB/001 (square nodes), (B): TLT/PHC/001 (diamond nodes), (C): TLT/PHC/002 (V-shape nodes), (D): TLT/PSK/001a (ellipse nodes; based on Thermo–Finnigan system) and (E): TLT/PSK/001b (ellipse nodes; based on Waters system). Italicized labels indicate putative assignments; numbers in bold indicate compounds 1 to 6 (refer to Figure 2). Edge thickness corresponds to relative cosine similarity between nodes.
Figure 4. Cyanobacterial molecular network clusters based on the seeded 1 (nodes outlined) with a cosine similarity cutoff of 0.65 obtained from extracts. (A): TLT/PHB/001 (square nodes), (B): TLT/PHC/001 (diamond nodes), (C): TLT/PHC/002 (V-shape nodes), (D): TLT/PSK/001a (ellipse nodes; based on Thermo–Finnigan system) and (E): TLT/PSK/001b (ellipse nodes; based on Waters system). Italicized labels indicate putative assignments; numbers in bold indicate compounds 1 to 6 (refer to Figure 2). Edge thickness corresponds to relative cosine similarity between nodes.
Marinedrugs 16 00505 g004
Marinedrugs 16 00505 g005 550
Figure 5. Overlaying of two molecular networks based on the extracts obtained from TLT/PHC/001 and TLT/PHC/002. Red and green nodes represent compounds produced by TLT/PHC/001 and TLT/PHC/002, respectively. Relative proportion of the compound belonging to each extract was represented as a pie chart in each node.
Figure 5. Overlaying of two molecular networks based on the extracts obtained from TLT/PHC/001 and TLT/PHC/002. Red and green nodes represent compounds produced by TLT/PHC/001 and TLT/PHC/002, respectively. Relative proportion of the compound belonging to each extract was represented as a pie chart in each node.
Marinedrugs 16 00505 g005
Table
Table 1. Mass differences (Δ) of observed and calculated m/z values of known compounds found in clusters, based on HRMS data obtained from Waters system.
Table 1. Mass differences (Δ) of observed and calculated m/z values of known compounds found in clusters, based on HRMS data obtained from Waters system.
ClusterCompoundMolecular Formula of Observed [M + Na]+Calculated m/zObserved m/zaΔ (ppm)
A1C32H48O10Na615.2987615.2988+0.1
2C32H46O9Na597.3040597.3047+1.1
3C33H50O10Na629.3302629.3308+0.9
4C32H47O10BrNa695.2236695.2227−0.4
5C31H46O10Na601.2989601.2982−1.2
B1C32H48O10Na615.2987615.3005+2.9
2C32H46O9Na597.3040597.3024−2.7
3C33H50O10Na629.3302629.3307+0.8
C1C32H48O10Na615.2987615.2986−0.1
2C32H46O9Na597.3040597.3058+3.0
3C33H50O10Na629.3302629.3290−1.9
E1C32H48O10Na615.2987615.2985−0.3
2C32H46O9Na597.3040597.3027−2.1
3C33H50O10Na629.3302629.3308+0.9
6C31H44O10Na599.2832599.2843+1.8
aΔ = Observed m/z − Calculated m/z.
Table
Table 2. A list of known aplysiatoxin and related analogues isolated from marine cyanobacteria with the m/z values of commonly occurring proton and sodium adducts [1,7,14,15,16].
Table 2. A list of known aplysiatoxin and related analogues isolated from marine cyanobacteria with the m/z values of commonly occurring proton and sodium adducts [1,7,14,15,16].
CompoundMolecular FormulaMolecular MassCalculated m/z [M + H]+Calculated m/z [M + Na]+
AplysiatoxinC32H47O10Br671.621672.622694.604
DebromoaplysiatoxinC32H48O10592.725593.727615.709
Oscillatoxin AC31H46O10578.698579.700601.682
Anhydro-19-bromoaplysiatoxinC32H44O9Br2732.502733.502755.484
Anhydro-19,21-dibromoaplysiatoxinC32H43O9Br3811.398812.398834.380
17-Bromooscillatoxin AC31H45O10Br657.594658.595680.577
17,19-Dibromooscillatoxin AC31H44O10Br2736.491737.491759.473
19-BromoaplysiatoxinC32H46O10Br2750.517751.518773.500
Oscillatoxin B1C32H46O10590.709591.711613.692
Oscillatoxin B2C32H46O10590.709591.711613.692
31-Noroscillatoxin BC31H44O10576.683577.684599.666
Oscillatoxin DC31H42O8542.668543.670565.651
3-MethoxyaplysiatoxinC33H49O10Br685.641686.648709.639
3-MethoxydebromoaplysiatoxinC32H46O10606.745607.753629.735
Nhatrangin AC21H32O8412.474413.482435.465
Nhatrangin BC21H31O8Br491.370492.378514.360
AnhydrodebromoaplysiatoxinC32H46O9574.702575.711597.693
30-Methyloscillatoxin DC32H44O8556.695557.696579.678
Manauealide AC32H47O10Cl627.168628.171650.153
Manauealide BC32H47O10Br671.621672.622694.604
Manauealide CC34H50O11634.762635.763657.745
Anhydrooscillatoxin AC31H44O9560.683561.685583.667
AnhydroaplysiatoxinC32H45O9Br653.606654.607676.589
Neo-debromoaplysiatoxin AC32H46O10590.703591.712613.693
Neo-debromoaplysiatoxin BC27H38O6458.588459.596481.578

Share and Cite

MDPI and ACS Style

Ding, C.Y.G.; Pang, L.M.; Liang, Z.-X.; Goh, K.K.K.; Glukhov, E.; Gerwick, W.H.; Tan, L.T. MS/MS-Based Molecular Networking Approach for the Detection of Aplysiatoxin-Related Compounds in Environmental Marine Cyanobacteria. Mar. Drugs 2018, 16, 505. https://0-doi-org.brum.beds.ac.uk/10.3390/md16120505

AMA Style

Ding CYG, Pang LM, Liang Z-X, Goh KKK, Glukhov E, Gerwick WH, Tan LT. MS/MS-Based Molecular Networking Approach for the Detection of Aplysiatoxin-Related Compounds in Environmental Marine Cyanobacteria. Marine Drugs. 2018; 16(12):505. https://0-doi-org.brum.beds.ac.uk/10.3390/md16120505

Chicago/Turabian Style

Ding, Chi Ying Gary, Li Mei Pang, Zhao-Xun Liang, Kau Kiat Kelvin Goh, Evgenia Glukhov, William H. Gerwick, and Lik Tong Tan. 2018. "MS/MS-Based Molecular Networking Approach for the Detection of Aplysiatoxin-Related Compounds in Environmental Marine Cyanobacteria" Marine Drugs 16, no. 12: 505. https://0-doi-org.brum.beds.ac.uk/10.3390/md16120505

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop