Next Article in Journal
Brevetoxin and Conotoxin Interactions with Single-Domain Voltage-Gated Sodium Channels from a Diatom and Coccolithophore
Next Article in Special Issue
Probing the Anti-Cancer Potency of Sulfated Galactans on Cholangiocarcinoma Cells Using Synchrotron FTIR Microspectroscopy, Molecular Docking, and In Vitro Studies
Previous Article in Journal
Induction of Apoptotic Cell Death in Human Leukemia U937 Cells by C18 Hydroxy Unsaturated Fatty Acid Isolated from Red Alga Tricleocarpa jejuensis
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Marine Drugs
Volume 19
Issue 3
10.3390/md19030139
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

Complementary Analytical Platforms of NMR Spectroscopy and LCMS Analysis in the Metabolite Profiling of Isochrysis galbana

by
Muhammad Safwan Ahamad Bustamam
1,
Hamza Ahmed Pantami
2,
Awanis Azizan
1,
Khozirah Shaari
1,2,
Chong Chou Min
3,
Faridah Abas
1,
Norio Nagao
3,
Maulidiani Maulidiani
4,
Sanjoy Banerjee
1,
Fadzil Sulaiman
1 and
Intan Safinar Ismail
1,2,*
1
Natural Medicine and Products Research Laboratory, Institute of Bioscience, Universiti Putra Malaysia, Serdang 43400, Selangor, Malaysia
2
Department of Chemistry, Faculty of Science, Universiti Putra Malaysia, Serdang 43400, Selangor, Malaysia
3
Department of Aquaculture, Faculty of Agriculture, Universiti Putra Malaysia, Serdang 43400, Selangor, Malaysia
4
Faculty of Science and Marine Environment, Universiti Malaysia Terengganu, Kuala Nerus 21030, Terengganu, Malaysia
*
Author to whom correspondence should be addressed.
Mar. Drugs 2021, 19(3), 139; https://0-doi-org.brum.beds.ac.uk/10.3390/md19030139
Submission received: 1 February 2021 / Revised: 25 February 2021 / Accepted: 26 February 2021 / Published: 2 March 2021
(This article belongs to the Special Issue Application of Spectroscopic Techniques in Marine Natural Products)
Browse Figures
Versions Notes

Abstract

:
This study was designed to profile the metabolites of Isochrysis galbana, an indigenous and less explored microalgae species. 1H Nuclear Magnetic Resonance (NMR) spectroscopy and Liquid Chromatography-Mass Spectrometry (LCMS) were used to establish the metabolite profiles of five different extracts of this microalga, which are hexane (Hex), ethyl acetate (EtOAc), absolute ethanol (EtOH), EtOH:water 1:1 (AqE), and 100% water (Aq). Partial least square discriminant analysis (PLS–DA) of the generated profiles revealed that EtOAc and Aq extracts contain a diverse range of metabolites as compared to the other extracts with a total of twenty-one metabolites, comprising carotenoids, polyunsaturated fatty acids, and amino acids, that were putatively identified from the NMR spectra. Meanwhile, thirty-two metabolites were successfully annotated from the LCMS/MS data, ten of which (palmitic acid, oleic acid, α-linolenic acid, arachidic acid, cholesterol, DHA, DPA, fucoxanthin, astaxanthin, and pheophytin) were similar to those present in the NMR profile. Another eleven glycerophospholipids were discovered using MS/MS-based molecular network (MN) platform. The results of this study, besides providing a better understanding of I. galbana’s chemical make-up, will be of importance in exploring this species potential as a feed ingredient in the aquaculture industry.

1. Introduction

Microalgae, as the name implies, are unicellular microscopic algae, which, depending on the species, can range from a couple of micrometres (µm) to a few hundred µm in size. They are typically found individually as well as in chains or groups in freshwater and marine systems [1]. Due to their nutritional values, microalgae have long been proposed as live feed in the aquaculture industry [2]. Among them, indigenous microalgae are widely preferred and utilized due to their high growth rate and robustness in survivability even when exposed to unpredictable and erratic tropical weathers [3]. Their lack of cell wall is also advantageous as it allows larvae, the juvenile form of animals, to easily digest and assimilate microalgae into their system [4,5,6]. Besides aquatic organisms, human beings can also benefit from supplementing the diet with microalgae as several studies have testified their health benefits in boosting the production of red blood cells along with their excellent performance as antioxidant, anti-inflammatory, anti-mutagenic, and antimicrobial agents [7].
Microalgae represent the first component of the aquatic food chain as they are a type of photosynthetic microorganism. This suggests their enormous importance and potential in generating abundant aquaculture yield by manipulating the food web causal relationship [8,9,10]. For instance, they are frequently used in the hatcheries of several aquatic animals, such as for some fish species at their primary stage of development, certain crustacean species at larval stage, various stages of bivalves (broodstock, larvae, juveniles), and post-set abalone, and as the main feeding material for zooplankton [11,12,13]. Being a natural food in the ecosystem, besides having significantly high nutritional values, makes microalgae a suitable choice of immunostimulants to be incorporated into the feed [14] as showcased by the positive effects of species such as Chlorella vulgaris and Spirulina platensis on the immunological parameters (lysozyme activity, phagocytic activity, complement activity, respiratory burst activity, and increased plasma protein) of certain fishes in multiple studies [15,16,17].
Isochrysis galbana, a golden-brown and flagellated marine microalga [18], is rich in polyunsaturated fatty acids (PUFA) such as docosahexaenoic acid (DHA) and eicosapentaenoic acid (EPA) [19] and is a valuable source of photosynthetic pigments like chlorophyll a and fucoxanthin [20], two abundantly available bioactive resources among microalgae [21,22,23,24]. Previous phytochemical studies on I. galbana selectively focused on carotenoids [25], fatty acids [26,27] and triacylglycerols [28]. Aguilera-Sáez et al. (2019) [29] recently performed an NMR-based metabolomics study on I. galbana in which amino acids, lipids, and sterols were profiled. These studies emphasised that the selection of solvents used for extraction is vital for the recovery of metabolites, a crucial stage in the workflow of a metabolomics analysis to capture the broadest number of metabolites possible [30].
Nuclear Magnetic Resonance (NMR) Spectroscopy and Liquid Chromatography-Mass Spectrometry (LCMS) are the two most common and robust analytical platforms used in metabolite profiling studies [31,32,33]. The advantages of NMR spectroscopy as an analytical tool include simple sample preparation, fast data acquisition, and high reproducibility [34,35]. On the other hand, LCMS offers other complementary features from its high sensitivity in metabolite detection (up to picomolar level), few problems with peak overlaps, and high reliability in the identification of different metabolites from a sample. Bearing in mind that there is no single analytical platform that can perform a complete identification and quantification of all molecules within a sample [36], the combination of these two powerful spectroscopic tools is certainly advantageous for the coverage, sensitivity, and reliability of the generated results. In addition, molecular networking (MN) has been an important bioinformatics platform to simulate and annotate untargeted mass spectrometry (MS) data [37,38] since its launch in 2012 [39]. This approach offers new ways of exploring the metabolome of biological samples by providing essential analog knowledge among the metabolites identified. Numerous MN applications in microalgae chemistry have been documented in recent years [40,41,42].
Hence, in this study, NMR fingerprint and LCMS profile of I. galbana extracted in five different solvents with different polarities were performed to determine the best solvent that can extract a wide range of metabolites in high concentrations. The results produced can pave the way towards developing a standardized biomass or extract from this microalgae species. To the best of our knowledge, there have been no studies done on the metabolite profiling of I. galbana using NMR in combination with LCMS as described herein.

2. Results and Discussion

2.1. Identification of Metabolites in NMR Spectra of Different Solvent Extracts

Figure 1 displays the representatives 1D 500 MHz 1H NMR spectra of I. galbana in Hex, EtOAc, EtOH, AqE and Aq extracts. The individual spectrum of each solvent extract was provided in the Supplementary File (Figure S6). Overall, the Aq spectrum showed fewer metabolite signals with low intensity as compared with the other solvents. This might be due to its low dissolution in the NMR solvent used (CD3OD:CDCl3 in the ratio of 8:2), unlike the other extracts. However, the yield of Aq extract was the highest (Table 1), while Hex extract gave the lowest yield which, suggests that I. galbana contains more polar metabolites than the less and non-polar ones. In the aromatic region, the signals were intense in EtOAc extract followed by EtOH and Hex extracts. Further interpretation of the assigned signals was completed using 1H NMR, while 2D-J-RES and HSQC experiments were utilized in order to increase the metabolites specificity and to minimize the congestion of the signals. The peaks were assigned by referring to previous studies [29,43,44,45] and by comparing them with freely available online databases such as Human Metabolome Database (HMDB) and PubChem.
A total of 21 metabolites from different chemical classes were putatively identified along with their chemical shift assignments, corresponding to multiplicities and coupling constants values as summarized in Table 2. The upfield region of the spectrum (0.8 to 3.0 ppm), which belongs to the aliphatic compounds, was mostly contributed by amino acids and fatty acids [47]. The five amino acids, which were identified as valine, isoleucine, leucine, alanine, and threonine, were clearly detected in the polar solvent, particularly in Aq extract as presented in Figure 1. These metabolites were previously reported in the same microalgae species [29] and also in another species known as Pleurochrysis carterae, which belongs to the same Haptophyta phylum as I. galbana [44].
The upfield region of all spectra, with the exception of the Aq extract, exhibited signals for both saturated and unsaturated fatty acids from the appearance of characteristic terminal methyl protons at 0.88 ppm. This signal was correlated to the last carbon attached to the methyl group at 16.4 ppm in the HSQC spectra. Meanwhile, the presence of PUFA with omega-3 FAs was confirmed with the triplet signal observed at 16.6 ppm from the terminal methyl protons of its last carbon and the HSQC cross-peaks for the carbon at 16.6 ppm [29], as shown in Figure 2a. The presence of saturated palmitic acid (C16:0) and arachidic acid (C20:0) was identified from the multiplet at 1.33 and 1.40 ppm, respectively. The unsaturated fatty acids were verified as oleic acid (18:1), α-linoleic acid (C18:3), docosahexaenoic fatty acid chain (DHA) (C22:6), and docosapentaenoic fatty acid chain (DPA) (22:5). The characteristic signal for DHA was assigned at 2.38 ppm, and it appeared as a small triplet which HSQC cross-peaks to a carbon (C-2) at 36.8 ppm in the HSQC spectra. This observation is in agreement with a study by Aguilera-Sáez et al. (2019) [29]. Another PUFA was assigned as docosapentaenoic fatty acid chain (DPA) from the multiplet peak at 1.38 ppm, which cross-peaks to C-20 at 31.2 ppm in HSQC. The J-resolved and HSQC spectra for both PUFAs (DHA and DPA) are displayed in Figure 2b. In previous studies, I. galbana was reported to contain variable levels of eicosapentaenoic fatty acid (EPA) [48,49]. However, EPA was not detected in the present study, which might be because of the high concentration of PUFAs. Culture conditions can highly affect the concentration of PUFAs [4], which may be the possible reason for the undetected EPA in this study. The level of EPA also depends on the elongation of EPA into DPA before subsequent desaturation into DHA by Δ6 desaturase occurs in the pathway of omega-3 long-chain-PUFA biosynthesis [50]. Another discovered metabolite, which was cholesterol, exhibited sterol signals in the region of 0.62–0.68 ppm, which is again similar to what was reported by Aguilera-Sáez et al. (2019) [29]. Its characteristic doublet of doublet signal was observed at 5.28 ppm after being compared with an online database.
The middle region of the spectra (3.5 and 5.5 ppm) was congested with signals of sugars, including carbohydrates, and other groups of metabolites that were discovered in Aq extract. Confirmation of the assigned metabolites was aided by information from 2D J-resolved and 1H-13C HSQC correlation, as displayed in Figure 3. These metabolites were mainly observed in the polar extracts (Aq, AqE, and EtOH). The presence of sucrose was justified by a characteristic singlet at 3.68 ppm, which cross-peaked to C-3 of cyclohexane ring at 74.5 ppm. Meanwhile, glucose was assigned based on the anomeric proton doublet of doublet peak at 3.52 ppm with its HSQC correlation to C-7 (66.5 ppm) attached to a hydroxyl group [45]. Choline, another primary metabolite found, was observed as a singlet at 3.22 ppm with its HSQC correlation to 57.0 ppm for its three methyl substituents attached to an amino skeletal [29]. The rest of the metabolites detected in this region were dimethylsulphoniopropionate (DMSP) and d-1,4/2,5-cyclohexanetetrol, which were putatively identified based on the characteristic singlet at 2.92 ppm and multiplet at 3.72 ppm with the correlation of HSQC to carbons at 28.3 ppm and 73.5 ppm (C-1,2,4,5) respectively, comparable to those reported by Aguilera-Sáez et al. (2019) [29].
The downfield region from 6.0 to 8.5 ppm displayed signals of aromatic resonances from the carotenoids and chlorophyllic constituents, specifically the olefinic protons (5.8–6.8 ppm), which had a few signals in the aliphatic region (0.8–2.5 ppm) [51]. These observations allowed the identification of several carotenoids, with much complications due to their small differences in structure; for example, the identification of fucoxanthin and violaxanthin was difficult and confusing due to the similarity in their olefinic chains. However, the presence of fucoxanthin was justified by the singlet of the terminal methyl protons at 1.37 ppm and 2.12 ppm for the cycloalkane moiety, which were not present in violaxanthin. Meanwhile, the presence of another carotenoid, astaxanthin, was confirmed via HSQC crossed-peak between proton signal at 4.15 ppm to the corresponding carbon (65.2 ppm) attached to a hydroxyl group. These carotenoids were mostly detected in the organic solvent extracts of EtOAc and EtOH, as displayed in the stacked spectra of 1H NMR (Figure 1). Another pigment identified as chlorophyll could be clearly spotted in the spectra from the isolated proton signals between 8.5 to 11.2 ppm [52]. In this study, all singlets that appeared between 8.60–9.21 ppm in the non-polar extracts were assigned as pheophytin a. Correlation with the nearest carbon between 96.2–107.2 ppm in HSQC spectra confirmed their identification, as previously reported [47].
Including I. galbana, there are four microalgae strains of the haptophyte in the order of Isochrysidales, which contain polyunsaturated long-chain alkenones (PULCAs) [53]. These unsaturated ketones are made up of a carbon chain of C37–C39 with 2 to 4 trans-type double bonds, which usually present at the intervals of seven carbon atoms [54]. Recently, Richter et al. (2019) [55] found C37–C39 alkenones in I. galbana. In the current study, the alkenones were characterized by the appearance of singlet methyl protons at 2.12 ppm. Meanwhile, the presence of trans-type double bond was confirmed by the signals of the allylic proton, which attached to carbon with a double bond at 1.95 ppm and crossed-peak to carbon at 35.4 ppm in the HSQC spectra, as reported by Iglesias et al. (2019) [46].
Overall, the 21 metabolites, including amino acids, carbohydrates, fatty acids, carotenoids, and chlorophyll, were successfully determined using NMR spectroscopy. These results verified the convenience of NMR in metabolomics profiling studies as almost all metabolites from different classes of compounds were detected in just one type of analysis [56]. Unfortunately, there is a limitation in using a single deuterated solvent: the combination of CD3OD and CDCl3 (8:2) was able to completely dissolve most of the extracts but not for Aq and Hex extracts due to the difference in polarity. Hence, this hindered the detection of the maximum number of compounds, as reported in other studies [29,46]. The selection of NMR solvent system was based on a preliminary optimization step using several different solvents combinations. The microalgae extracts were found to dissolve better in CD3OD and CDCl3 (8:2) solvent system as compared to the others. The metabolites detected in the present study cover a wide range of compound classes. Furthermore, it is mandatory to use the same solvent system in both NMR acquisition and data bucketing/binning for multivariate data analysis. All spectra were aligned and binned together before the solvent regions in each spectrum were uniformly excluded. Using a different solvent system for each solvent extract may seem wise so as to achieve high solubility; however, this will bring about an inconsistency during data processing for multivariate data analysis.

2.2. Discriminative Analysis of Solvent Extracts Based on 1H NMR Data

Partial least square-discriminant analysis (PLS-DA) was applied to classify the predefined groups of I. galbana solvent extracts and generate information on the discriminating metabolites among the groups [57]. This technique is a supervised method, requiring class label information in building the appropriate model for data interpretation. A model is interpreted as good if its cumulative R2 and Q2 values are close to 1, with R2 giving an overview of the model fitness, while Q2 describes the predictive quality of the model [58]. The validity of the models is evaluated using a permutation test (100 permutations), wherein, for a valid model, the R2 and Q2 intercepts should not exceed 0.3–0.4 and 0.05, respectively [58].
The PLS-DA score plot (Figure 4a) shows that the Aq, AqE, and EtOH extracts were discriminated from the EtOAc and Hex extracts by principal component (PC) 1. The R2 and Q2 values were 0.98 and 0.97 respectively, indicating a model of good fit, while the R2 and Q2 intercepts were between 0.21 and −0.49, respectively, further proving the robustness of this model (see Supporting Information Figure S1). Meanwhile, the metabolites that were responsible for the separation could be identified from the loading plot (Figure 4b) with their respective chemical shift binned to the width of 0.04 ppm. The metabolites such as DMSP, d-1,4/2,5-cyclohexanetetrol, choline, sucrose and glucose that were located on the positive side of PC1 are more prominent in Aq, AqE, and EtOH extracts, while those on the negative side of PC1 such as fatty acids and carotenoids were higher in EtOAc and Hex extracts. The importance and significance of the variables were determined by analysing the Variable Importance in Projection (VIP) plot with jack-knifing uncertainty bars. The variables with VIP scores greater than 1 with an error bar not crossing the baseline in the loading column plot were retained as significant and were classified as chemical markers that gave an influential contribution to the discrimination in the PLS model [59]. Based on their VIP values being greater than 1.0, thirty-two binned regions, as illustrated in the VIP plot (see Supporting Information Figure S2), significantly contributed to the class separation. Some metabolites were assigned based on several corresponding binned regions, for example, five NMR peaks (2.7, 2.9, 2.94, 3.42 and 3.46) for DMSP alone.

2.3. Relative Quantification

In order to view the distribution of the identified metabolites among different solvent extracts, the relative concentration of the metabolites was analysed and visualized using HCA. The binned regions with VIP values greater than 1 were selected from the original data set, normalized with Pareto scaling and subjected to HCA using Euclidean distance measures and Ward’s clustering algorithm. The results of the analysis were visualized as a color-coded heat map (Figure 5), which helps to reflect the distribution and relative intensities of the metabolites.
The heat map indicated that EtOAc and Hex extracts contained higher concentrations of fatty acids and carotenoids as compared to the other extracts. The Aq extract had the highest levels of DMSP, D-1,2,4,5-cyclohexanetetrol, choline, sucrose, and glucose, followed by the AqE and EtOH extracts. The heat map was in agreement with the discrimination of metabolites seen in the PLS-DA loadings plot (Figure 4b). Moreover, most of the putatively identified carotenoids and unsaturated fatty acids were discriminated on the negative side of PC1, verifying the high intensity (concentration) of these metabolites in EtOAc and Hex extracts.
The heat map also shows that the technical replicates of the same origin were grouped closely together, indicating a good reproducibility in both the sample extraction and 1H NMR measurement. The HCA dendrogram shows the samples are divided into two distinct clusters comprising EtOAc and Hex extracts in one and Aq, AqE, and EtOH extracts in another. The EtOH and AqE extracts were further grouped together and separated from the Aq extract. The only difference between EtOAc and Hex extracts was the intensity of carotenoids (fucoxanthin and astaxanthin), wherein these compounds were higher in EtOAc than Hex based on the colour intensity in the map. Hence, the HCA and PLS-DA results strongly suggest that ethyl acetate (EtOAc) is the best solvent to be used in order to extract most of the valuable fatty acids and carotenoids in high concentration from I. galbana.
To further validate the results, important metabolites were relatively quantified based on the peak area values of their characteristic chemical shifts before being analysed using two-way ANOVA test (with a threshold of p < 0.05). Out of the thirty-two metabolites identified from the VIP plot, only 10 compounds were deemed to be important in causing the group separation as their concentrations varied in a statistically significant manner. As for quantification, only peaks that did not overlap with other peaks and referred to a single metabolite were selected, such that at 3.66 ppm (s) for sucrose, 3.54 ppm (dd) for glucose, 3.22 ppm (s) for choline, 2.90 ppm (s) for DMSP, 1.82 (m) for d-1,4/2,5-cyclohexanetetrol, 1.42 ppm (m) for arachidic acid, 1.34 ppm (m) for palmitic acid, 2.82 ppm (m) for α-linolenic acid, 0.66 ppm (m) for cholesterol, and 6.62 (d) for astaxanthin. The derivatives of the carotenoids and fatty acids were excluded from the quantification process since these metabolites do not clearly display any characteristic peaks.
Relative quantification of the selected metabolites was presented in a bar chart along with their statistical difference (Figure 6). There was a similar trend among the different extracts with regard to the quantification of sucrose, glucose, choline, DMSP, and d-1,4/2,5-cyclohexanetetrol. The more polar extracts, particularly the Aq, followed by AqE and EtOH, demonstrated a high significance difference (p < 0.05) for these five metabolites. Sucrose, glucose, and d-1,4/2,5-cyclohexanetetrol were also found in EtOAc extract at 85.7%, 84.2%, and 60.6%, which were lower than that in the Aq extract. A study on Ajwa dates reported the same observations, whereby simple sugars such as sucrose and glucose, are better extracted by polar solvents like methanol at concentrationss of 50% and 100% [60]. Choline and DMSP, on the other hand, were not obvious in EtOAc and Hex extracts due to the very small signal-to-noise ratio.
Nevertheless, another 5 compounds from the class of fatty acids and carotenoids, namely arachidic acid, palmitic acid, α-linolenic acid, cholesterol, and astaxanthin, were observed the highest (p < 0.05) in EtOAc extract. Comparing between both the EtOAc and Hex extracts shows that arachidic acid and cholesterol were higher in the former with concentrations of 36.4% and 69.2%, respectively, while α-linolenic acid was higher in the latter, with the concentration of 18.9%. The effect of different solvents on the efficacy of extracting fucoxanthin was previously investigated, and ethyl acetate was shown to be more efficient than chloroform and ethanol [61]. In another study, the chloroform extract of Chaetoceros calcitrans was shown to give a higher amount of palmitic acid and cholesterol, which suggests that moderately polar solvents like ethyl acetate and chloroform are the best for fatty acids extraction [45]. Most of the fatty acids were readily detected in the Aq and AqE extracts. Thus, results from PLS-DA, HCA, and relative quantification established that EtOAc is the best solvent to extract high amounts of fatty acids and carotenoids, while amino acids and carbohydrates are more efficiently extracted using water. Hence, EtOAc and water were selected as the most suitable solvents for the recovery of high concentrations of metabolites from I. galbana, as it is crucial to obtain a wide range of compounds in the highest intensity possible [30]. This result is in agreement with the previous study on fucoxanthin, which was extracted efficiently in EtOAc solvent when compared with the other solvents [62].
In the extraction process, solvents used must have the same polarity as the solute of interest in order to effectively dissolve the targeted metabolites [63]. Different classes of solutes need certain types of extraction solvents since their chemical properties can be in different polarity and can also be thermally unstable. The extraction process is a crucial initial stage to ensure the bioactive compounds in a medicinal plant are preserved prior to further analysis [64]. The principle of solvent selection relies most significantly on the specific characteristics of the targeted compound(s) to be isolated [64]. In the present study, most of the bioactive compounds reported in I. galbana such as polyunsaturated fatty acids (PUFA), docosahexaenoic acid (DHA), and photosynthetic pigments like chlorophyll a and fucoxanthin [20] were observed to be highly present in EtOAc extract.

2.4. UHPLC–MS/MS Analysis

LCMS had been previously used to identify fatty acids, triacylglycerols, and carotenoids in I. galbana [25,28]. However, the system used was limited on the scan range and used atmospheric pressure chemical ionization (APCI) as the ionization method instead of electrospray ionization (ESI). In the present study, a UHPLC system coupled with Q-Exactive Focus Orbitrap LCMS/MS, which is a more advanced system known for its higher sensitivity and excellent mass accuracy, was utilized to achieve a more comprehensive metabolite profile of this microalgae. The total ion chromatogram (TIC) (Figure 7a,b) shows the peaks detected in the positive and negative ion modes for EtOAc extract. This extract was selected for further analysis since its NMR data showed most of the metabolite groups in high concentrations. The mass analysis was conducted in switching mode [65,66]. Using both positive and negative ionization modes to identify compounds enables a broader coverage of the metabolome than the use of a single polarity [67,68], as some analytes can only be detected in either one of these modes. The spectral interpretation software predicts and automatically generates detailed fragmentation according to the general principles of ionization, fragmentation, and rearrangement using chemical structure provided by databases such as HMDB and PubChem.

Metabolite Identification in Positive and Negative Ion Mode

The identity, retention time, UV characteristic, and fragment ion(s) for each metabolite are presented in Table 3, and the MS/MS spectra for every molecular ion detected are available in the Supporting Information (Figure S3). Thirty-two metabolites comprising eight carotenoids, five chlorophylls, four glycerophospholipids, three sphingolipids, one glycerolipid, fatty acid, and sterol were annotated in the positive mode. For the negative mode, nine main peaks were characterized as fatty acids. All metabolites detected showed high mass accuracy with their scored mass error being less than 10 ppm. However, both arachidic acid and cholesterol displayed high mass error values (>400 ppm), indicating poor accuracy as compared with their theoretical values. Nonetheless, these two metabolites’ fragments, which were ambiguously matched with the HMDB database were positively detected in NMR. Most of the carotenoids were discovered in the photodiode-array (PDA) spectrum (Figure 7c). The identification of fucoxanthin as a major carotenoid with a small mass error value was found in the positive mode at m/z 659.4283 [M + H]+ and at tR 12.07 min of the chromatogram by comparing with PubChem database, whereas the fragments of m/z 581, 411, 355, and 199 were matched with Massbank database.
The series of fragmentation pattern was proposed in Figure 8, starting with fragment 1 (via F1), which displayed the loss of H2O molecule from the parent ion, and this was followed by the inductive cleavage of the carboxylic group to yield m/z 581. The second fragmentation (F2) was suggested due to the cleavage at the carbon-oxygen bond in m/z 581 to yield m/z 411. Besides this, there were two other fragmentations (F3 and F4) that occured from the cleavage of the carbon-carbon bond of the long carbon chain and generated m/z 199 [M + H − 460]+ and m/z 355 [M + H − 304]+, respectively. These results are in agreement with the NMR data, which identified fucoxanthin as the main carotenoid in EtOAc extract.
Astaxanthin (Figure 9) was detected at tR 10.11 min with isopropanol adduct [M + IsoProp + H]+ exhibiting its ion at m/z 657.4136 [M + 61]+. Some characteristic MS fragments observed were matched with the Massbank database of astaxanthin. The first fragment (F1: Figure 9) arose from the inductive cleavage at carbon bond to produce m/z 417. This was followed by a loss of an H2O molecule, due to the hydrogen being attacked by a hydroxyl, and the cleavage of CO moiety, which resulted in the ion at m/z 371 [M + H − 180 − 18 − 28]+. A second fragment via F2 was derived after the loss of CO moiety and H2O and before the inductive cleavage at the long carbon chain yielding ion m/z 279 [M + H − 28 − 18 − 272]+. The next fragmentation (F3) was due to the loss of oxygen from the hydronium ion formation plus the cyclohexyl group [M + H − 16 − 154]+ ion cleavage, which resulted in m/z 427. Another major carotenoid, phoenicoxanthin, was detected at tR 13.63 min with its [M + H]+ ion at m/z 581.3964. This carotenoid was reported in Aurantiochytrium sp., a colorless fungus-like algae, which contains a similar composition of carotenoids and PUFAs as I. galbana [69].
Another major fatty acid was putatively identified as DHA at 12.49 min with a small mass error (−0.31 ppm) from the theoretical mass in the PubChem database. The presence of DHA, DPA, α-linolenic acid, palmitic acid, arachidic acid, and oleic acid was justified by the NMR results. Four additional acids, namely EPA, hexacosanedioic acid, 3,6,9,12,15-octadecapentaenoic acid, and stearidonic acids, which were not identified in the NMR results, were detected in the mass data. These fatty acids had been previously reported in the marine haptophytes, dinophytes, and prasinophytes [70] and are one of the main PUFAs found in Isochrysis zhangjiangensis [71].
Chlorophyll a was detected in the positive mode with an ion at m/z 893.5519 (tR 20.39 min). The fragmentation pattern of F1 (Figure 10) suggested an inductive cleavage at the long carbon chain, which yielded m/z 555, followed by the loss of oxygen at carbonyl group, which produced m/z 539 [M + H − 338 − 16]+. The second fragment (F2) was derived due to cleavage of a long carbon chain to yield m/z 615. Subsequent fragmentation occurred from m/z 615 with the loss of methoxy group (32 amu), giving an ion at m/z 583 [M + H − 278 − 32]+. This observation is in agreement with a previous large-scale microalgae production study that identified chlorophyll a as the main chlorophyll of I. galbana [72]. In another study, Crupi et al. (2013) [73] detected pheophytin a, one of the main pigments in I. galbana, by using HPLC-DAD-MS.

2.5. MS/MS-Based Molecular Networking

Molecular networking facilitates a fast comparison of mass spectrometry profiles from complicated crude extracts for successful metabolites dereplication and exploration of novel compounds, which requires a high resolution of mass spectrometry data (MS/MS) [39]. Dereplication is defined as the fast detection of defined metabolites through the comparison of experimental mass spectra with libraries. All metabolites are represented as parent ions, which are linked by the chemical fragmentation of the compound. Related compounds comprised similar parent ion fragmentation patterns, which are represented as a cosine score from 1 (extremely similar fragmentation spectra) to 0 (totally different parent ions) [39,74]. Therefore, the parent ions (nodes) are bound by edges with cosine score value, resulting in the classification of analogous or structurally related compounds in molecular clusters [39,75].
In the present study, a classical global molecular network was generated based on UPLC–MS/MS data from EtOAc extract with blank using the GNPS platform in order to focus on the lipid content of I. galbana more comprehensively as microalgae become one of the promising sources of lipid production [76]. The putative identified metabolites were achieved by manual dereplication matched with several external databases, namely HMDB, PubChem, LIPID MAPS, and Chemspider through Metabolomics Workbench platform (www.metabolomicsworkbench.org, accessed on 28 October 2020) with the lowest mass error since automated dereplication on the GNPS platform was limited and did not match any known compound. Overall, there were thirteen clusters, generated with more than two nodes per cluster as shown in the Supplementary File (Figure S5). However, most of the clusters were not fully annotated due to the limited databases search, except for one of the largest clusters, which was identified as families of glycerophospholipids and consisted of eleven putative metabolites as shown in Figure 11. Other glycerolipids such as monogalactosyldiacylglycerols (MGDG) and digalactosyldiacylglycerols (DGDS) might be identified from other clusters if other databases are explored. This will be a time-consuming task that should be thoroughly taken up in future studies.
Glycerophospholipids are made up of phosphatidic acids, compounds constructed by a glycerol molecule with fatty acid esterification in two of its hydroxyl groups and phosphoric acid esterification in the third hydroxyl. Commonly, one of the phosphate moieties’ free OH groups is esterified with another element, forming various glycerophospholipids [77]. In the current cluster, ten of them are classified as phosphatidylethanolamine (PE), since a glycerol replacement site is dominated by the phosphorylethanolamine moiety as referred to in Table 4. Unlike the others, the first metabolite, known as lysophospholipid, refers to a phospholipide that is lacking one of the two O-acyl chains. Like diacylglycerols, PE can have several combinations of fatty acids with differing lengths and saturation at positions C-1 and C-2, with the most popular are fatty acids containing 16, 18 and 20 carbons. The third structure is annotated as diglyceride (DG), consisting of two fatty acid chains, which are stearidonic acid at the C-1 and C-2 positions covalently bonded to a glycerol molecule through ester linkages without phosphorylethanolamine moiety. This finding has proven the existence of several phospholipids assigned as PE and phosphatidylcholine (PC) in the lipophilic extracts of I. galbana detected by two-dimensional 1H−31P HSQC TOCSY (heteronuclear single quantum coherence-total correlation spectroscopy) [29]. Moreover, the presence of PE can be observed significantly in other species of microalgae such as Chaetoceros gracilis, Nannochloropsis gaditana and Picochlorum atomus, especially after three days of cultivation under phosphate repletion study [78].
The aim of MN is to compare and group all MS/MS spectra in an individual or series of extracts based on their similarity expressed as cosine score value [39]. In other words, MN is capable of clustering molecules according to their structural characteristics, as its MS/MS spectrum is related to the chemical structure of the fragmented metabolites [79]. Figure 12a presents the MS/MS spectrum of LysoPE(22:4(7Z,10Z,13Z,16Z)/0:0) depicting common fragments found in MS/MS spectra of metabolites number 1 to 7 after they shared several fragments such as m/z 241 derived from glycerol and phosphate components free from its ethanolamine moiety and another fragment of the hydrocarbon chain-forming m/z 299 as proposed by Mass Frontier software. Meanwhile, metabolites numbers 8 to 11 share other common fragments, as shown in Figure 12b of the represented MS/MS spectrum of PE(16:0/18:0), showing a fragment of m/z 211 forming as glycerol and phosphate backbone with hydrogen rearrangement at carbonyl carbon and m/z 285 established as cleavage of saturated hydrocarbon chain at one of its glycerol ester oxygens. The MS/MS spectrum for all metabolites are provided as Supplementary Files (Figure S4).

2.6. Correlation between NMR and UHPLC–MS/MS Data

Advances in analytical methods of high-resolution nuclear magnetic resonance (NMR) spectroscopy and mass spectrometry (MS), with the incorporation of chemometric tools, have driven the field of metabolomics promptly for the high reproducibility of the former and high sensitivity and selectivity of the latter [80,81]. This analytical strategy in metabolomics had been applied in fields like toxicology, drug discovery, early disease detection, and food and nutrition sciences in recent years [82,83,84].
The combination of more than one analytical platform minimizes the shortcomings of using NMR or MS alone [85]. Hyphenation of NMR and MS will permit the exhibition of a diverse range of metabolites [86,87], albeit with a few duplicates, as both eventually offer a comprehensive identification with improved reliability [87,88,89,90]. Furthermore, this approach helps in optimizing the detection of unknown analytes by merging unique NMR information (chemical shifts, coupling constants) with that of MS (exact mass, molecular fragments) [91,92]. A quick approach in biological samples analysis can be achieved using NMR spectroscopy [93], as it can produce fingerprints for samples screening and classification [94]. Further profiling of the significant group of metabolites can then be done using MS.
In the current study, 21 metabolites, comprising amino acids, carbohydrates, fatty acids, and carotenoids, were successfully determined using NMR spectroscopy. Metabolites with a VIP value of more than 1 (Figure 4b) were considered significantly responsible for the separation of the different extracts. The heat map (Figure 5) demonstrates the highest intensity of fatty acids and carotenoids in EtOAc extract, while most of the amino acids were detected in the Aq extract. These results prove the superiority of NMR in identifying the different classes of compounds in fingerprinting analysis.
Low sensitivity is an inherent weakness of NMR spectroscopy, and to make up for this, the EtOAc extract was subjected to MS analysis, a more sensitive analytical tool with a detection level ranging from picomole to femtomole, to further profile the secondary metabolites [95,96]. This high sensitivity, however, may lead to a complicated forest of signals, and MS is also not a universal technique that can detect a wide range of metabolite classes in a single run; to do so, different chromatographic techniques are required [56]. The present study putatively identified thirty-two secondary metabolites by MS, which include carotenoids, chlorophylls, fatty acids, glycerophospholipids, and sphingolipids. The presence of 10 of the identified metabolites, namely palmitic acid, oleic acid, α-linolenic acid, arachidic acid, cholesterol, DHA, DPA, fucoxanthin, astaxanthin, and pheophytin a, was verified in the EtOAc extract based on the comparison of the MS data with the NMR. Other carotenoids, such as halocynthiaxanthin, diatoxanthin, and echinenone, as well as chlorophylls, such as pheophorbide a and chlorophyll a, were discovered only by MS technique. Apart from that, MS/MS-based molecular networking approach, which had never been done on I. galbana, has succeeded in the discovery of more glycerophospholipids which were classified as phosphatidylethanolamine (PE).
Both NMR and MS analytical platforms had been applied in many studies on algae, such as in the structural elucidation of carotenoids stereoisomers from Chlorococcum humicola, a green freshwater algae [97]. In this study, the structure of six closely related carotenoids was successfully characterized and quantified from a mixture. In another study, NMR and GC-MS techniques in combination were used to investigate the quality of lipids and corresponding FAMEs composition in Chlorella vulgaris, Spirulina platensis and Tetraselmis affchuii after treatment in different media composition, to explore their biodiesel potential [98].

3. Materials and Methods

3.1. Microalgae Culture and Harvest

Isolated Isochrysis galbana species (UPMC-A0009) was obtained from the Microalgal Production Laboratory, Aquatic Animal Health Unit, Faculty of Veterinary Medicine, Universiti Putra Malaysia in April, 2017. The stock culture, grown in Conway medium prepared in filtered (5 μm) and sterilized seawater [99], was maintained at 23 °C in an environmental chamber (Sanyo, Osaka, Japan) under 12–12 artificial light–dark cycle (light intensity 150 µmol/m2/s). The cultivation of I. galbana was initiated from a 25 mL aliquot of the stock culture in Erlenmeyer flasks (5 L) and was gradually scaled up to 90 L under continuous aeration in an annular photobioreactor of 120 L capacity. Biomass weight and cell count were used to monitor the microalgae growth. The former value was determined by dividing the dry weight of the filtered biomass with the filtrate volume [100]. Cell counting was performed daily, on a well-mixed sample, using a Neubauer hemocytometer (Assistant, Germany). The microalgae were harvested at the late-exponential growth phase (15 days), by centrifuging them at 10,000× g-force following Aguilera-Sáez et al. (2019) [29] with some modifications. The harvested biomass was then freeze-dried (Scanvac, Lynge, Denmark) and stored at −80 °C before use.

3.2. Preparation of Solvent Extracts

Extracts of the dried biomass were prepared in hexane (Hex), ethyl acetate (EtOAc), absolute ethanol (EtOH), EtOH:water 1:1 (AqE), and 100% water (Aq), in 6 replicates for each solvent, giving a total of thirty extracts of five different solvent polarities. Briefly, 200 mg of the dried biomass was mixed with 30 mL solvent in a 50 mL Schott test tube and sonicated in an ultrasonic bath sonicator (Kudos, Shanghai, China) for 30 min, taking care to keep the bath temperature maintained below 35 °C to avoid any unwanted decomposition or degradation of the compounds. The extract was filtered, and the residue was extracted again with fresh solvent. The extract filtrates were pooled, rotary-evaporated at 40 °C, freeze-dried, and stored under −20 °C prior to further analysis. The yield of each solvent extract was obtained from the average value of 6 sample replicates extracted in their respective solvents and freeze-dried.

3.3. Spectroscopic Measurements

Spectroscopic measurements of all extracts were separately prepared for 1H NMR spectroscopy and ultra-high-performance liquid chromatography-mass spectrometry (UHPLCMS/MS).

3.3.1. 1H NMR Analysis

Ten milligrams of each extract was vortex-mixed in microcentrifuge tubes with 700 µL of CD3OD:CDCl3 at a ratio of 8:2, with 0.05% trimethylsilylpropanoic acid (TSP) as an internal reference standard. The mixture was sonicated for 15 min at room temperature and centrifuged at 18,900× g-force for 10 min. About 600 µL of the clear supernatant was transferred to a 5 mm NMR tube for data acquisition [101]. 1H-NMR spectra were analysed using 500 MHz Varian Unity INOVA NMR spectrometer (Varian Inc., Palo Alto, CA, USA) functioning at a frequency of 499.91 MHz and maintained at 26 °C. For data acquisition, a single-pulse proton experiment with PRESAT was used with 21.0 µs pulse width, 2-s relaxation delay, 3.53 total acquisition time for 64 scans. Two-dimensional J-resolved experiment (JRES) was used to help clarify the spectral assignment. The time for the J-resolved spectrum acquisition was 50 min and 18 s, with 8 scans per 256 increments for the axis of the spin-spin coupling constant and spectral widths of 66 Hz, and 8 K data points for the chemical shift axis with spectral widths of 8012.8 Hz. The relaxation delay was set at 1.0 s. Heteronuclear single quantum coherence (HSQC) spectra were obtained using 16 scans, 1 K data points, and 256 increments at the spectral width of 13 ppm and 220 ppm for the proton and carbon dimensions, respectively. The relaxation delay was 1.0 s, giving an achievement time of 6 h, 9 min, and 9 s. The 2D NMR spectral processing for structural elucidation was carried out using MestRenova software (version 6.02-5475, Mestrelab Research, Santiago de Compostella, Spain).

3.3.2. LCMS/MS Analysis

Samples were prepared at the concentration of 2 mg/mL in LCMS-grade methanol and ultrasonicated for 10 min before being filtered through a nylon membrane (0.22 mm) into a 2 mL screw-capped sample vial to remove any precipitation. LCMS analysis was performed on a Thermo ScientificTM Q ExactiveTM Hybrid Quadrupole-Orbitrap mass spectrometer coupled to a Dionex Ultimate 3000 UHPLC system (Thermo Fisher Scientific Inc., Waltham, MA, USA), which is fitted with Acquity UPLC BEH C18 column (1.7 µm × 2.1 mm × 100 mm) (Waters, Milford, MA, USA). The mobile phase used was 0.1% formic acid in deionized water (solvent A) and 0.1% formic acid in LCMS-grade acetonitrile (solvent B). The injection volume was 5 µL, analysis time was set to 30 min, and the flow rate was 0.25 mL/min. The gradient program commenced with 10% solvent B at 0 min, 20% at 1.00 min, 30% at 2.00 min, 70% at 7.00 min, 80% at 10.00 min, 90% at 12.00 min, and 100% at 13.00–30.00 min. Molecular ion identification was obtained in switching electrospray ionization modes with a full scan range of m/z 100–1500 amu. Other MS parameters were set as follows: collision energy of 30 eV, spray voltage of 4.2 kV, capillary temperature at 350 °C, sheath gas flow rate of 50, and auxiliary nitrogen (99% pure) gas flow rate of 10. The mass resolution was set to 70,000 full widths at half maximum (FWHM). The UV detectors were set to 254, 280, 400, and 440 nm, while the PDA detector was set to 190–600 nm. Metabolite assignments were done using the retention time, UV-vis spectra, and MS data (accurate mass, negative and positive ion modes) from the compounds analysed using Thermo Xcalibur 2.0 (Thermo Fisher Scientific Inc., Waltham, MA, USA) and their comparison was performed using literature data and standard online databases (freely available), such as Metabolomics Workbench, Human Metabolome Database (HMDB), PubChem, MassBank, and Metlin. The fragmentation for all compounds discussed herein was based on the pattern derived from HighChem Mass Frontier 3.0 (Thermo Fisher Scientific Inc., Waltham, MA, USA).

3.4. Data Processing and Multivariate Data Analysis

Phasing and baseline corrections of NMR spectra were carried out using Chenomx software (version 5.1, Edmonton, AB, Canada). All NMR spectra were phase-adjusted and baseline-corrected automatically and referenced to the internal standard (TSP) at 0.00 ppm. The 1H NMR spectrum of each sample was processed and bucketed (bin width of 0.04 ppm) from the spectral region of 0.50 to 10.00 ppm. The peaks for residual water (4.80–4.90 ppm), methanol (3.30–3.32 ppm), and chloroform (7.78–7.79 ppm) were excluded from the spectral data to retain the signals from endogenous metabolites. A total of 243 integrated regions were obtained for each spectrum. The generated dataset was converted to ASCII files and imported to SIMCA-P 13.0 software (Umetrics, Umeå, Sweden) for multivariate data analysis and visualization of results. The dataset was Pareto-scaled, and correlations among the samples were then established by partial least square-discriminant analysis (PLS-DA). A score plot was constructed to visualize the separation between groups, while a loading plot was used to identify metabolites that contributed to the groups’ separation. The model was validated using the default seven-fold internal cross-validation based on the goodness-of-fit (R2X) and goodness-of-prediction (Q2) values together with the 100-permutation test. Hierarchical cluster analysis (HCA) was performed using MetaboAnalyst 3.0 (http://www.metaboanalyst.ca, accessed on 21 Novermber 2019), a public web-based platform for comprehensive analysis of metabolomics data.

3.5. Molecular Networking

The molecular networks based on MS/MS data were generated using the online workflow Global Natural Products Social Molecular Networking (GNPS) platform (http://gnps.ucsd.edu, accessed on 28 October 2020) with a registered account. Prior to uploading the data into GNPS, the raw MS data including blank were converted into mzXML format using MSConvert software downloaded from Proteowizard website (http://proteowizard.sourceforge.net/tools.shtml, accessed on 28 October 2020). Then, the converted data files were uploaded to GNPS using FileZilla 3.42.1 software (https://filezilla-project.org/, accessed on 28 October 2020). In the GNPS data analysis workflow, sample and blank data were selected as G1 and G2, respectively, with precursor ion mass tolerance set to 0.02 Da and a fragment ion mass tolerance of 0.02 Da. A network was processed with edges that were filtered to have a cosine score above 0.7 and a minimum 6 matched peaks [71]. Upon processing, the result was downloaded and the network was visualized using ChemViz 1.3 plugin (freely available at http://www.cgl.ucsf.edu/cytoscape/chemViz/, accessed on 28 October 2020) within Cytoscape 3.7.1 software (Institute of Systems Biology, Seattle, WA, USA).

3.6. Statistical Analysis

One-way analysis of variance (ANOVA) was performed using GraphPad Prism 6.0 (GraphPad Software, San Diego, CA, USA). Post-hoc analysis was carried out using Tukey’s test, wherein values with p ≤ 0.05 were considered to be statistically significant. Values were expressed as mean ± standard deviation (SD).

4. Conclusions

To the best of our knowledge, this is the first report describing the comprehensive metabolites profiling of I. galbana by NMR and supplemented with LCMS. The extraction of I. galbana in five different solvent polarities successfully exhibited the highest concentrations of fatty acids and carotenoids in EtOAc extract and amino acids and carbohydrates in Aq extract, as determined by the NMR fingerprints. Hence, both solvents are preferable for the extraction of these microalgae to obtain a wide range of compounds. The identification of important metabolites in EtOAc extract was further confirmed using UHPLC–MS/MS analysis, from which another 32 metabolites comprising fatty acids, carotenoids, glycerophospholipids, and sphingolipids were identified, with 10 of them present in both of the analyses. Further identification using MN platform on MS/MS data has discovered more glycerophospholipids that are classified as PE. The NMR fingerprinting together with MS profiling managed to fruitfully characterize a large set of metabolites with a wide range of classes while improving the accuracy of the ones identified. The understanding of metabolites content in I. galbana brought about in this study can be utilized for the preparation of a standardized feed for selected aquaculture assays.

Supplementary Materials

The following are available online at https://0-www-mdpi-com.brum.beds.ac.uk/1660-3397/19/3/139/s1, Figure S1: (a) R2 and Q2 values from summary of fit model and (b) the R2 and Q2 intercepts derived from 100 permutation tests for each of the solvent extracts; Figure S2: VIP plot with values more than 1.0 and error bars not crossing baseline; Figure S3: MS/MS spectra of molecular ion detected by LCMS/MS in positive and negative mode; Figure S4: MS/MS spectra of metabolites from a cluster of glycerophospholipids identified in EtOAc extract of I. galbana in positive mode; Figure S5: Full molecular network of EtOAc extract of I. galbana; Figure S6: Representatives 1D 500 MHz 1H NMR individual spectrum.

Author Contributions

Conceptualization of the study and experimental design, I.S.I., K.S., C.C.M., and F.A.; microalgae source and reagents S.B. and N.N.; sample preparation and spectroscopic analyses, M.S.A.B.; NMR data interpretation and data analysis, M.S.A.B., H.A.P., and A.A.; LCMS analysis, M.S.A.B. and M.M.; manuscript writing and editing, M.S.A.B., F.S., I.S.I., and K.S., financial procurement. K.S. and I.S.I. All authors have read and agreed to the published version of the manuscript.

Funding

This research has been funded by Japan Science and Technology Agency (JST)/Japan International Cooperation Agency (JICA) with Ministry of Higher Education Malaysia (MOHE) through joint research program Science and Technology Research Partnership for Sustainable Development (SATREPS) in the project for Continuous Operation System for Microalgae Production Optimized for Sustainable Tropical Aquaculture (COSMOS) under research grant scheme 6300856-12038.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Acknowledgments

We would like to acknowledge Gregory Tan Guan Yuan, Department of Aquaculture, Faculty of Agriculture, Universiti Putra Malaysia, for his guidance during microalgae cultivation and Salahudin Mohd Raof of the Spectroscopy Unit, Laboratory of Natural Products, Institute of Bioscience, Universiti Putra Malaysia, for technical assistance in NMR analysis.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Thurman, H.V. Introductory Oceanography; Prentice Hall College: Upper Saddle River, NJ, USA, 1997. [Google Scholar]
  2. Meireles, L.; Guedes, A.; Malcata, F. Lipid class composition of the microalga Pavlova lutheri: Eicosapentaenoic and docosahexaenoic acids. J. Agric. Food Chem. 2003, 51, 2237–2241. [Google Scholar] [CrossRef] [PubMed]
  3. Babuskin, S.; Radhakrishnan, K. Effect of photoperiod, light intensity and carbon sources on biomass and lipid productivities of Isochrysis galbana. Biotechnol. Lett. 2014, 36, 1653–1660. [Google Scholar] [CrossRef]
  4. Devos, M.; Poisson, L.; Ergan, F.; Pencreac’h, G. Enzymatic hydrolysis of phospholipids from Isochrysis galbana for docosahexaenoic acid enrichment. Enzym. Microb. Technol. 2006, 39, 548–554. [Google Scholar] [CrossRef]
  5. Lin, Y.; Chang, F.; Tsao, C.; Leu, J.Y. Influence of growth phase and nutrient source on fatty acid composition of Isochrysis galbana CCMP 1324 in a batch photoreactor. Biochem. Eng. J. 2007, 37, 166–176. [Google Scholar] [CrossRef]
  6. Liu, J.; Sommerfeld, M.; Hu, Q. Screening and characterization of Isochrysis strains and optimization of culture conditions for docosahexaenoic acid production. Appl. Microbiol. Biotechnol. 2013, 97, 4785–4798. [Google Scholar] [CrossRef]
  7. Nuno, K.; Vilarruel-Lopez, A.; Puebla-Perez, A.M.; Romero-Velarde, E.; Puebla-Mora, A.G.; Ascencio, F. Effects of the marine microalgae Isochrysis galbana and Nannochloropsis oculata in diabetic rats. J. Funct. Foods 2013, 5, 106–115. [Google Scholar] [CrossRef]
  8. Benemann, J. Microalgae aquaculture feeds. J. Appl. Phycol. 1992, 4, 233–245. [Google Scholar] [CrossRef]
  9. Brown, M. Nutritional value and use of microalgae in aquaculture. In Proceedings of the Avances en Nutrición Acuícola VI. Memorias del VI Simposium Internacional de Nutrición Acuícola, Cancún, Mexico, 3–6 September 2002; pp. 281–292. [Google Scholar]
  10. Pulz, O.; Gross, W. Valuable products from biotechnology of microalgae. Appl. Microbiol. Biotechnol. 2004, 65, 635–648. [Google Scholar] [CrossRef] [PubMed]
  11. Utting, S.; Millican, P. Techniques for the hatchery conditioning of bivalve broodstocks and the subsequent effect on egg quality and larval viability. Aquaculture 1997, 155, 45–54. [Google Scholar] [CrossRef]
  12. Rico-villa, B.; Woerther, P.; Mingant, C.; Pouvreau, S.; Hamon, M.; Robert, R. A flow-through rearing system for ecophysiological studies of pacific oyster Crassostrea gigas larvae. Aquaculture 2008, 282, 54–60. [Google Scholar] [CrossRef] [Green Version]
  13. Ragg, N.; King, N.; Watts, E.; Morrish, J. Optimising the delivery of the key dietary diatom Chaetoceros calcitrans to intensively cultured Greenshell (TM) mussel larvae, Perna canaliculus. Aquaculture 2010, 306, 270–280. [Google Scholar] [CrossRef]
  14. Kwak, J.H.; Baek, S.H.; Woo, Y.; Han, J.K.; Kim, B.G.; Kim, O.Y.; Lee, J.H. Beneficial immunostimulatory effect of short-term Chlorella supplementation: Enhancement of natural killer cell activity and early inflammatory response (randomized, double- blinded placebo-controlled trial). Nutr. J. 2012, 11, 53. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Adel, M.; Yeganeh, S.; Dadar, M.; Sakai, M.; Dawood, M.A.O. Effects of dietary Spirulina platensis on growth performance, humoral and mucosal immune responses and disease resistance in juvenile great sturgeon (Huso huso Linnaeus, 1754). Fish Shellfish Immunol. 2016, 56, 436–444. [Google Scholar] [CrossRef] [PubMed]
  16. Fadl, S.E.; Elgohary, M.S.; Elsadany, A.Y.; Gad, D.M.; Hanaa, F.F.; El-habashi, N.M. Contribution of microalgae-enriched fodder for the Nile tilapia to growth and resistance to infection with Aeromonas hydrophila. Algal Res. 2017, 27, 82–88. [Google Scholar] [CrossRef]
  17. Galal, A.A.A.; Reda, R.M.; Mohamed, A.A. Influences of Chlorella vulgaris dietary supplementation on growth performance, hematology, immune response and disease resistance in Oreochromis niloticus exposed to sub-lethal concentrations of penoxsulam herbicide. Fish Shellfish Immunol. 2018, 77, 445–456. [Google Scholar] [CrossRef] [PubMed]
  18. Grima, E.M.; Perez, J.A.U.; Camacho, F.G.; Fernbndez, F.G.A.; Alonso, D.L.; Castillo, C.I.S. Preservation of the marine microalga, Isochrysis galbana: Influence on the fatty acid profile. Aquaculture 1994, 123, 377–385. [Google Scholar] [CrossRef]
  19. Cho, J.Y.; Jin, H.; Lim, H.J.; Whyte, J.N.C.; Hong, Y. Growth activation of the microalga Isochrysis galbana by the aqueous extract of the seaweed Monostroma nitidum. J. Appl. Phycol. 1999, 10, 561–567. [Google Scholar] [CrossRef]
  20. Valenzuela-Espinoza, E.; Millan-Nunez, R.; Nunez-Cebrero, F. Protein, carbohydrate, lipid and chlorophyll a content in Isochrysis aff. galbana (clone T-Iso) cultured with a low cost alternative to the f/2 medium. Aquac. Eng. 2002, 25, 207–216. [Google Scholar] [CrossRef]
  21. Goh, L.P.; Loh, S.P.; Fatimah, M.Y.; Perumal, K. Bioaccessibility of carotenoids and tocopherols in marine microalgae, Nannochloropsis sp. and Chaetoceros sp. Malays. J. Nutr. 2009, 15, 77–86. [Google Scholar]
  22. Mimouni, V.; Ulmann, L.; Pasquet, V.; Mathieu, M.; Picot, L.; Bougaran, G.; Cadoret, J.-P.; Morant-Manceau, A.; Schoefs, B. The potential of microalgae for the production of bioactive molecules of pharmaceutical interest. Curr. Pharm. Biotechnol. 2012, 13, 2733–2750. [Google Scholar] [CrossRef]
  23. Zakar, T.; Laczko-Dobos, H.; Toth, T.N.; Gombos, Z. Carotenoids assist in cyanobacterial photosystem II assembly and function. Front. Plant Sci. 2016, 7, 295. [Google Scholar] [CrossRef] [Green Version]
  24. Matos, J.; Cardoso, C.; Bandarra, N.M.; Afonso, C. Microalgae as a healthy ingredient for functional food: A review. Food Funct. 2017, 8, 2672–2685. [Google Scholar] [CrossRef]
  25. Ahmed, F.; Fanning, K.; Netzel, M.; Turner, W.; Li, Y.; Schenk, P.M. Profiling of carotenoids and antioxidant capacity of microalgae from subtropical coastal and brackish waters. Food Chem. 2014, 165, 300–306. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  26. Batista, A.P.; Gouveia, L.; Bandarra, N.M.; Franco, J.M.; Raymundo, A. Comparison of microalgal biomass profiles as novel functional ingredient for food products. Algal Res. 2013, 2, 164–173. [Google Scholar] [CrossRef] [Green Version]
  27. Maadane, A.; Merghoub, N.; Ainane, T.; El Arroussi, H.; Benhima, R.; Amzazi, S.; Bakri, Y.; Wahby, I. Antioxidant activity of some Moroccan marine microalgae: Pufa profiles, carotenoids and phenolic content. J. Biotechnol. 2015, 215, 13–19. [Google Scholar] [CrossRef]
  28. Nedbalová, L.; Střížek, A.; Sigler, K.; Řezanka, T. Effect of salinity on the fatty acid and triacylglycerol composition of five haptophyte algae from the genera Coccolithophora, Isochrysis and Prymnesium determined by LC-MS/APCI. Phytochemistry 2016, 130, 64–76. [Google Scholar] [CrossRef] [PubMed]
  29. Aguilera-Sáez, L.M.; Abreu, A.C.; Camacho-Rodríguez, J.; González-López, C.V.; Del Carmen Cerón-García, M.; Fernández, I. NMR Metabolomics as an Effective Tool to Unravel the Effect of Light Intensity and Temperature on the Composition of the Marine Microalgae Isochrysis galbana. J. Agric. Food Chem. 2019, 67, 3879–3889. [Google Scholar] [CrossRef]
  30. Lu, W.; Su, X.; Klein, M.S.; Lewis, I.A.; Fiehn, O.; Rabinowitz, J.D. Metabolite Measurement: Pitfalls to Avoid and Practices to Follow. Annu. Rev. Biochem. 2017, 86, 277–304. [Google Scholar] [CrossRef]
  31. Ibanez, C.; Simo, C.; Barupal, D.; Al, E. A new metabolomic workflow for early detection of Alzheimer’s disease. J. Chromatogr. A 2013, 1302, 65–71. [Google Scholar] [CrossRef] [PubMed]
  32. Farag, M.; Wessjohann, L. Metabolome classification of commercial hypericum perforatum (St. John’s Wort) preparations via UPLC-qTOF-MS and chemometrics. Planta Med. 2012, 78, 488–496. [Google Scholar] [CrossRef] [PubMed]
  33. Ali, K.; Iqbal, M.; Yuliana, N.; Al, E. Identification of bioactive metabolites against adenosine A1 receptor using NMR-based metabolomics. Metabolomics 2013, 9, 778–785. [Google Scholar] [CrossRef]
  34. Schripsema, J. Application of NMR in plant metabolomics: Techniques, problems and prospects. Phytochem. Anal. 2010, 21, 14–21. [Google Scholar] [CrossRef] [PubMed]
  35. Kim, H.; Choi, Y.; Verpoorte, R. NMR-based plant metabolomics: Where do we stand, where do we go? Trends Biotechnol. 2011, 29, 267–275. [Google Scholar] [CrossRef]
  36. Pan, Z.; Raftery, D. Comparing and combining NMR spectroscopy and mass spectrometry in metabolomics. Anal. Bioanal. Chem. 2007, 387, 525–527. [Google Scholar] [CrossRef]
  37. Quinn, R.A.; Nothias, L.F.; Vining, O.; Meehan, M.; Esquenazi, E.; Dorrestein, P.C. Molecular networking as a drug discovery, drug metabolism and precision medicine strategy. Trends Pharmacol. Sci. 2017, 38, 143–154. [Google Scholar] [CrossRef]
  38. Traxler, M.F.; Kolter, R. A massively spectacular view of the chemical lives of microbes. Proc. Natl. Acad. Sci. USA 2012, 109, 10128–10129. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  39. Watrous, J.; Roach, P.; Alexandrov, T.; Heath, B.S.; Yang, J.Y.; Kersten, R.D.; van der Voort, M.; Pogliano, K.; Gross, H.; Raaijmakers, J. Mass spectral molecular networking of living microbial colonies. Proc. Natl. Acad. Sci. USA 2012, 109, E1743–E1752. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  40. Pantami, H.A.; Bustamam, M.S.A.; Lee, S.Y.; Ismail, I.S.; Faudzi, S.M.M.; Nakakuni, M.; Shaari, K. Comprehensive GCMS and LC-MS/MS metabolite profiling of chlorella vulgaris. Mar. Drugs 2020, 18, 367. [Google Scholar] [CrossRef]
  41. Lv, M.; Lo, C.; Hsu, C.C.; Wang, Y.; Chiang, Y.R.; Sun, Q.; Wu, Y.; Li, Y.; Chen, L.; Yu, C.P. Identification of enantiomeric byproducts during microalgae-mediated transformation of metoprolol by MS/MS spectrum based networking. Front. Microbiol. 2018, 9, 2115. [Google Scholar] [CrossRef] [PubMed]
  42. Buedenbender, L.; Astone, F.A.; Tasdemir, D. Bioactive molecular networking for mapping the antimicrobial constituents of the baltic brown alga fucus vesiculosus. Mar. Drugs 2020, 18, 311. [Google Scholar] [CrossRef]
  43. Wei, L.; Mediani, A.; Khaleeda, N.; Zolke, Z.; Wei, S.; Sa, I.; Shaari, K.; Abas, F. Phytochemistry Letters Phytochemical diversity of Clinacanthus nutans extracts and their bioactivity correlations elucidated by NMR based metabolomics. Phytochem. Lett. 2015, 14, 123–133. [Google Scholar]
  44. Chengxu, Z.; Jie, L.; Yangfang, Y.; Xiaojun, Y.; Baoning, L.; Xin, W. The metabolite profi ling of coastal coccolithophorid species Pleurochrysis carterae (Haptophyta). Chin. J. Oceanol. Limnol. 2016, 34, 749–756. [Google Scholar]
  45. Azizan, A.; Safwan, M.; Bustamam, A.; Maulidiani, M.; Shaari, K.; Ismail, I.S.; Nagao, N.; Abas, F. Metabolite Profiling of the Microalgal Diatom Chaetoceros Calcitrans and Correlation with Antioxidant and Nitric Oxide Inhibitory Activities via 1H NMR-Based Metabolomics. Mar. Drugs 2018, 16, 154. [Google Scholar] [CrossRef] [Green Version]
  46. Iglesias, M.J.; Soengas, R.; Probert, I.; Guilloud, E.; Gourvil, P.; Mehiri, M.; López, Y.; Cepas, V.; Gutiérrez-del-río, I.; Redondo-blanco, S.; et al. NMR characterization and evaluation of antibacterial and antiobiofilm activity of organic extracts from stationary phase batch cultures of five five marine microalgae (Dunaliella sp., D. salina, Chaetoceros calcitrans, C. gracilis and Tisochrysis lutea). Phytochemistry 2019, 164, 192–205. [Google Scholar] [CrossRef]
  47. Pariyani, R.; Ismail, S.; Azam, A.; Abas, F.; Shaari, K. Identification of the compositional changes in Orthosiphon stamineus leaves triggered by different drying techniques using 1H NMR metabolomics. J. Sci. Food Agric. 2017, 97, 4169–4179. [Google Scholar] [CrossRef] [PubMed]
  48. Fidalgo, J.P.; Cid, A.; Torres, E.; Sukenik, A.; Herrero, C. Effects of nitrogen source and growth phase on proximate bio- chemical composition, lipis classes and fatty acid profile of the marine microalga Isochryis galbana. Aquaculture 1998, 166, 105–116. [Google Scholar] [CrossRef] [Green Version]
  49. Tzovenis, I.; de Pauw, N.; Sorgeloos, P. Optimisation of T-ISO biomass production rich in essential fatty acids: I. Effect of different light regimes on growth and biomass production. Aquaculture 2003, 216, 203–222. [Google Scholar] [CrossRef]
  50. Adarme-Vega, T.C.; Lim, D.K.Y.; Timmins, M.; Vernen, F.; Li, Y.; Schenk, P.M. Microalgal biofactories: A promising approach towards sustainable omega-3 fatty acid production. Microb. Cell Fact. 2012, 11, 96. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  51. Englert, G.; Bjørnland, T.; Liaaen-Jensen, S. 1D and 2D NMR study of some allenic carotenoids of the fucoxanthin series. Magn. Reson. Chem. 1990, 28, 519–528. [Google Scholar] [CrossRef]
  52. Abraham, R.J.; Rowan, A.E. Nuclear magnetic resonance spectroscopy of Chlorophyll. In Chlorophylls; Hugo, S., Ed.; CRC Press: Boca Raton, FL, USA, 1991; pp. 797–834. [Google Scholar]
  53. Nakamura, H.; Sawada, K.; Araie, H.; Suzuki, I.; Shiraiwa, Y. n-Nonacosadienes from the marine haptophytes Emiliania huxleyi and Gephyrocapsa oceanica. Phytochemistry 2015, 111, 107–113. [Google Scholar] [CrossRef] [Green Version]
  54. Rechka, J.A.; Maxwell, J.R. Characterisation of Alkenone Temperature Indicators in Sediments and Organisms. Org. Geochem. 1988, 13, 727–734. [Google Scholar] [CrossRef]
  55. Richter, N.; Longo, W.M.; George, S.; Shipunova, A.; Huang, Y.; Amaral, Z.L. Phylogenetic diversity in freshwater-dwelling Isochrysidales haptophytes with implications for alkenone production. Geobiology 2019, 17, 272–280. [Google Scholar] [CrossRef] [Green Version]
  56. Emwas, A.-H.M.; Salek, R.M.; Griffin, J.L.; Merzaban, J. NMR-based metabolomics in human disease diagnosis: Applications, limitations, and recommendations. Metabolomics 2013, 9, 1048–1072. [Google Scholar] [CrossRef]
  57. Chen, C.; Krausz, K.W.; Shah, Y.M.; Idle, J.R.; Frank, J. Gonzalez Serum Metabolomics Reveals Irreversible Inhibition of Fatty Acid β-Oxidation through the Suppression of PPARα Activation as a Contributing Mechanism of Acetaminophen-Induced Hepatotoxicity. Chem. Res. Toxicol. 2009, 22, 699–707. [Google Scholar] [CrossRef] [Green Version]
  58. Eriksson, L.; Johansson, E.; Kettaneh-Wold, N.; Trygg, J.; Wikstrom, C.; Wold, S. Multi- and Megavariate Data Analysis Part 1: Basic Principles and Applications; Umetrics Academy: Umeå, Sweden, 2006. [Google Scholar]
  59. Woo, H.M.; Kim, K.M.; Choi, M.H.; Jung, B.H.; Lee, J.; Kong, G.; Nam, S.J.; Kim, S.; Bai, S.W.; Chung, B.C. Mass spectrometry based metabolomic approaches in urinary biomarker study of women’s cancers. Clin. Chim. Acta 2009, 400, 63–69. [Google Scholar] [CrossRef] [PubMed]
  60. Abdul-hamid, N.A.; Abas, F.; Ismail, I.S.; Shaari, K.; Lajis, N.H. Influence of Different Drying Treatments and Extraction Solvents on the Metabolite Profile and Nitric Oxide Inhibitory Activity of Ajwa Dates. J. Food Sci. 2015, 80, 603–611. [Google Scholar] [CrossRef] [PubMed]
  61. Kim, S.M.; Jung, Y.; Kwon, O. A Potential Commercial Source of Fucoxanthin Extracted from the Microalga Phaeodactylum tricornutum. Appl. Biochem. Biotechnol. 2012, 166, 1843–1855. [Google Scholar] [CrossRef] [PubMed]
  62. Min, S.; Kang, K.S.; Donghwa, O.K.; Pan, C.C. Fucoxanthin as a major carotenoid in Isochrysis aff. galbana: Characterization of extraction for commercial application. J. Korean Soc. Appl. Biol. Chem. 2012, 55, 477–483. [Google Scholar] [CrossRef]
  63. Altemimi, A.; Lakhssassi, N.; Baharlouei, A.; Watson, D.G.; Lightfoot, D.A. Phytochemicals: Extraction, isolation, and identification of bioactive compounds from plant extracts. Plants 2017, 6, 42. [Google Scholar] [CrossRef]
  64. Bart, H.J.; Pilz, S. Industrial Scale Natural Products Extraction, 1st ed.; Wiley-VCH Verlag GmbH & Co.: Weinheim, Germany, 2011. [Google Scholar]
  65. Tolonen, A.; Uusitalo, J. Fast screening method for the analysis of total flavonoid content in plants and foodstuffs by high-performance liquid chromatography/electrospray ionization time-of-flight mass spectrometry with polarity switching. Rapid Commun. Mass Spectrom. 2004, 18, 3113–3122. [Google Scholar] [CrossRef]
  66. Cai, F.; Xu, W.; Wei, H.; Sun, L.; Gao, S.; Yang, Q.; Feng, J.; Zhang, F.; Chen, W. Simultaneous determination of active xanthone glycosides, timosaponins and alkaloids in rat plasma after oral administration of Zi-Shen Pill extract for the pharmacokinetic study by liquid chromatography-tandem mass spectrometry. J. Chromatogr. B 2010, 878, 1845–1854. [Google Scholar] [CrossRef] [PubMed]
  67. Farag, M.; Huhman, D.; Lei, Z.; Sumner, L. Metabolic profiling and systematic identification of flavonoids and isoflavonoids in roots and cell suspension cultures of Medicago truncatula using HPLC-UV-ESI-MS and GC-MS. Phytochemistry 2007, 68, 342–354. [Google Scholar] [CrossRef] [PubMed]
  68. Nordström, A.; Want, E.; Northen, T.; Lehtiö, J.; Siuzdak, G. Multiple ionization mass spectrometry strategy used to reveal the complexity of metabolomics. Anal. Chem. 2008, 80, 421–429. [Google Scholar] [CrossRef]
  69. Yoshida, M.; Tanabe, Y.; Yonezawa, N.; Watanabe, M.M. Energy innovation potential of oleaginous microalgae. Biofuels 2012, 3, 761–781. [Google Scholar] [CrossRef]
  70. Ghioni, C.; Porter, A.E.A.; Sadler, I.H.; Tocher, D.R.; Sargent, J.R. Cultured fish cells metabolize octadecapentaenoic acid (all-cis δ3,6,9,12,15–18:5) to octadecatetraenoic acid (all-cis δ6,9,12,15–18:4) via its 2-trans intermediate (trans δ2, all-cis δ6,9,12,15–18:5). Lipids 2001, 36, 145–153. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  71. Wu, S.; Zhou, J.; Xin, Y.; Xue, S. Nutritional stress effects under different nitrogen sources on the genes in microalga Isochrysis zhangjiangensis and the assistance of Alteromonas macleodii in releasing the stress of amino acid deficiency. J. Phycol. 2015, 51, 885–895. [Google Scholar] [CrossRef]
  72. Di Lena, G.; Casini, I.; Lucarini, M.; Lombardi-Boccia, G. Carotenoid profiling of five microalgae species from large-scale production. Food Res. Int. 2018, 120, 810–818. [Google Scholar] [CrossRef]
  73. Crupi, P.; Toci, A.T.; Mangini, S.; Wrubl, F.; Rodolfi, L.; Tredici, M.R.; Coletta, A.; Antonacci, D. Determination of fucoxanthin isomers in microalgae (Isochrysis sp.) by high-performance liquid chromatography coupled with diode-array detector multistage mass spectrometry coupled with positive electrospray ionization. Rapid Commun. Mass Spectrom. 2013, 27, 1027–1035. [Google Scholar] [CrossRef]
  74. Frank, A.M.; Bandeira, N.; Shen, Z.; Tanner, S.; Briggs, S.P.; Smith, R.D.; Pevzner, P.A. Clustering millions of tandem mass spectra. J. Proteome Res. 2008, 7, 113–122. [Google Scholar] [CrossRef]
  75. 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. MS/MS networking guided analysis of molecule and gene cluster families. Proc. Natl. Acad. Sci. USA. 2013, 110, E2611–E2620. [Google Scholar] [CrossRef] [Green Version]
  76. Alishah Aratboni, H.; Rafiei, N.; Garcia-Granados, R.; Alemzadeh, A.; Morones-Ramírez, J.R. Biomass and lipid induction strategies in microalgae for biofuel production and other applications. Microb. Cell Fact. 2019, 18, 1–17. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  77. Blanco, G.; Blanco, A. Chapter 5—Lipids. In Medical Biochemistry; Academic Press: Cambridge, MA, USA, 2017; pp. 99–119. [Google Scholar]
  78. Cañavate, J.P.; Armada, I.; Hachero-Cruzado, I. Interspecific variability in phosphorus-induced lipid remodelling among marine eukaryotic phytoplankton. New Phytol. 2017, 213, 700–713. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  79. Allard, P.M.; Péresse, T.; Bisson, J.; Gindro, K.; Marcourt, L.; Pham, V.C.; Roussi, F.; Litaudon, M.; Wolfender, J.L. Integration of Molecular Networking and In-Silico MS/MS Fragmentation for Natural Products Dereplication. Anal. Chem. 2016, 88, 3317–3323. [Google Scholar] [CrossRef] [PubMed]
  80. Lindon, J.C.; Holmes, E.; Bollard, M.E.; Stanley, E.G.; Nicholson, J.K. Metabonomics technologies and their applications in physiological monitoring, drug safety assessment and disease diagnosis. Biomarkers 2004, 9, 1–31. [Google Scholar] [CrossRef]
  81. Wagner, S.; Scholz, K.; Donegan, M.; Burton, L.; Wingate, J.; Volkel, W. Metabonomics and Biomarker Discovery: LC−MS Metabolic Profiling and Constant Neutral Loss Scanning Combined with Multivariate Data Analysis for Mercapturic Acid Analysis. Anal. Chem. 2006, 78, 1296–1305. [Google Scholar] [CrossRef]
  82. Watkins, S.M.; German, J.B. Metabolomics and biochemical profiling in drug discovery and development. Curr. Opin. Mol. Ther. 2002, 4, 224–228. [Google Scholar] [PubMed]
  83. Waters, N.J.; Holmes, E.; Williams, A.; Waterfield, C.J.; Farrant, R.D.; Nicholson, J.K. NMR and pattern recognition studies on the time-related metabolic effects of alpha-naphthylisothiocyanate on liver, urine, and plasma in the rat: An integrative metabonomic approach. Chem. Res. Toxicol. 2001, 14, 1401–1412. [Google Scholar] [CrossRef]
  84. Goodacre, R.; Vaidyanathan, S.; Dunn, W.B.; Harrigan, G.G.; Kell, D.B. Metabolomics by Numbers: Acquiring and Understanding Global Metabolite Data. Trends Biotechnol. 2004, 5, 245–252. [Google Scholar] [CrossRef]
  85. Marshall, D.D.; Powers, R. Beyond the paradigm: Combining mass spectrometry and nuclear magnetic resonance for metabolomics. Prog. Nucl. Magn. Reson. Spectrosc. 2017, 100, 1–16. [Google Scholar] [CrossRef] [Green Version]
  86. Ling, Y.S.; Linag, H.J.; Chung, M.; Lin, M.; Lin, C. NMR- and MS-based metabolomics: Various organ responses following naphthalene intervention. Mol. Biosyst. 2014, 10, 1918–1931. [Google Scholar] [CrossRef]
  87. Marshall, D.; Lei, S.; Worley, B.; Huang, Y.; Garcia-Garcia, A.; Franco, R.; Dodds, E.; Powers, R. Combining DI-ESI-MS and NMR Datasets for Metabolic Profiling. Metabolomics 2015, 11, 391–402. [Google Scholar] [CrossRef] [Green Version]
  88. Bingol, K.; Brüschweiler, R. NMR/MS Translator for the Enhanced Simultaneous Analysis of Metabolomics Mixtures by NMR Spectroscopy and Mass Spectrometry: Application to Human Urine. J. Proteome Res. 2015, 14, 2642–2648. [Google Scholar] [CrossRef]
  89. Baker, J.; Ward, J.; Beale, M. Combined NMR and flow injection ESI-MS for Brassicaceae metabolomics. Methods Mol. Biol. 2012, 860, 177–191. [Google Scholar]
  90. Bingol, K.; Brüschweilerr-Li, L.; Li, D.; Zhang, B.; Xie, M.; Brüschweiler, R. Emerging new strategies for successful metabolite identification in metabolomics. Bioanalysis 2016, 8, 557–573. [Google Scholar] [CrossRef] [Green Version]
  91. Prichystal, J.; Schug, K.; Lemr, K.; Novak, J.; Havlicek, V. Structural Analysis of Natural Products. Anal. Chem. 2016, 88, 10338–10346. [Google Scholar] [CrossRef]
  92. Yang, Z. Online hyphenated liquid chromatography-nuclear magnetic resonance spectroscopy-mass spectrometry for drug metabolite and nature product analysis. J. Pharm. Biomed. Anal. 2006, 40, 516–527. [Google Scholar] [CrossRef]
  93. Van der Kooy, F.; Maltese, F.; Choi, Y. Quality control of herbal material and phytopharmaceuticals with MS and NMR based metabolic fingerprinting. Planta Med. 2009, 75, 763–775. [Google Scholar] [CrossRef] [Green Version]
  94. Windarsih, A.; Rohman, A.; Swasono, R.T. Application of H-NMR metabolite fingerprinting and chemometrics for the authentication of Curcuma longa adulterated with Curcuma manga. J. Appl. Pharm. Sci. 2018, 8, 75–81. [Google Scholar] [CrossRef]
  95. Bedair, M.; Sumner, L. Current and emerging mass-spectrometry technologies for metabolomics. Trends Anal. Chem. 2008, 27, 238–250. [Google Scholar] [CrossRef]
  96. Sumner, L.W.; Mendes, P.; Dixon, R.A. Plant metabolomics: Large-scale phytochemistry in the functional genomics era. Phytochemistry 2003, 62, 817–836. [Google Scholar] [CrossRef] [Green Version]
  97. Sivathanu, B.; Palaniswamy, S. Purification and characterization of carotenoids from green algae Chlorococcum humicola by HPLC-NMR and LC-MS-APCI. Biomed. Prev. Nutr. 2012, 2, 276–282. [Google Scholar] [CrossRef]
  98. As, S.; Costa, I.C.R.; Teixeira, C.M.L.L.; Filocomo, D.; Candido, R.; Silva, P.R.M.; Cunha, V.S.; Romeu, J. Investigation of Biodiesel Potential of Biomasses of Microalgaes Chlorella, Spirulina and Tetraselmis by NMR and GC-MS Techniques. J. Biotechnol. Biomater. 2016, 6, 1–15. [Google Scholar] [CrossRef] [Green Version]
  99. Tompkins, J.; DeVille, M.; Day, J.; Turner, M. Culture Collection of Algae and Protozoa: Catalogue of Strains; Titus Wilson and Son Ltd.: Kendal, UK, 1995. [Google Scholar]
  100. Lavens, P.; Sorgeloos, P. Manual on the production and use of live food for aquaculture. In FAO Fisheries Technical Paper No. 361; FAO: Rome, Italy, 1996. [Google Scholar]
  101. Kim, H.K.; Choi, Y.H.; Verpoorte, R. NMR-based metabolomic analysis of plants. Nat. Protoc. 2010, 5, 536–549. [Google Scholar] [CrossRef] [PubMed]
Marinedrugs 19 00139 g001 550
Figure 1. Representatives 1D 500 MHz 1H NMR spectra of hexane (Hex), ethyl acetate (EtOAc), absolute ethanol (EtOH), 50% ethanol (AqE), and aqueous (Aq) extracts of I. galbana dissolved in CD3OD:CDCl3 in the ratio of 8:2, with 0.05% trimethylsilylpropanoic acid (TSP) as an internal reference standard. The assignments of the peaks are listed in Table 2.
Figure 1. Representatives 1D 500 MHz 1H NMR spectra of hexane (Hex), ethyl acetate (EtOAc), absolute ethanol (EtOH), 50% ethanol (AqE), and aqueous (Aq) extracts of I. galbana dissolved in CD3OD:CDCl3 in the ratio of 8:2, with 0.05% trimethylsilylpropanoic acid (TSP) as an internal reference standard. The assignments of the peaks are listed in Table 2.
Marinedrugs 19 00139 g001
Marinedrugs 19 00139 g002a 550Marinedrugs 19 00139 g002b 550
Figure 2. Representative 2D-J-resolved (a) and HSQC spectra (b) of EtOAc extract. The signals were assigned as follows: 1, non-PUFA; 2, PUFA; 3, DHA; and 4, DPA.
Figure 2. Representative 2D-J-resolved (a) and HSQC spectra (b) of EtOAc extract. The signals were assigned as follows: 1, non-PUFA; 2, PUFA; 3, DHA; and 4, DPA.
Marinedrugs 19 00139 g002aMarinedrugs 19 00139 g002b
Marinedrugs 19 00139 g003 550
Figure 3. 2D-J-resolved (a) and HSQC spectra (b) of Aq extract in the region 2.9 to 3.8 ppm. The signals were assigned as follows: 1, sucrose; 2, glucose; 3, choline; 4, DMSP; and 5, d-1,4/2,5-cyclohexanetetrol.
Figure 3. 2D-J-resolved (a) and HSQC spectra (b) of Aq extract in the region 2.9 to 3.8 ppm. The signals were assigned as follows: 1, sucrose; 2, glucose; 3, choline; 4, DMSP; and 5, d-1,4/2,5-cyclohexanetetrol.
Marinedrugs 19 00139 g003
Marinedrugs 19 00139 g004 550
Figure 4. PLS-DA (a) scores and (b) loading scatter plots of the 1H NMR data representing I. galbana extracted in 5 different solvents; Aqueous (Aq), 50% ethanol (AqE), ethanol (EtOH), ethyl acetate (EtOAc) and hexane (Hex). Corresponding metabolites with VIP value more than 1:1, fatty acid derivatives; 2, DMSP; 3, d-1,4/2,5-cyclohexanetetrol; 4, amino acid derivatives; 5, choline; 6, valine; 7, threonine/DPA/fucoxanthin; 8, sucrose; 9, archidic acid; 10, leucine/α-linoleic acid; 11, palmitic acid; 12, α-linoleic acid; 13, astaxanthin; 14, cholesterol; 15, glucose; 16, carotenoid derivatives; 17, sucrose/DHA/DPA; 18, oleic acid.
Figure 4. PLS-DA (a) scores and (b) loading scatter plots of the 1H NMR data representing I. galbana extracted in 5 different solvents; Aqueous (Aq), 50% ethanol (AqE), ethanol (EtOH), ethyl acetate (EtOAc) and hexane (Hex). Corresponding metabolites with VIP value more than 1:1, fatty acid derivatives; 2, DMSP; 3, d-1,4/2,5-cyclohexanetetrol; 4, amino acid derivatives; 5, choline; 6, valine; 7, threonine/DPA/fucoxanthin; 8, sucrose; 9, archidic acid; 10, leucine/α-linoleic acid; 11, palmitic acid; 12, α-linoleic acid; 13, astaxanthin; 14, cholesterol; 15, glucose; 16, carotenoid derivatives; 17, sucrose/DHA/DPA; 18, oleic acid.
Marinedrugs 19 00139 g004
Marinedrugs 19 00139 g005 550
Figure 5. Heat map of identified metabolites of I. galbana extracted in 5 different solvents—Aqueous (Aq), 50% ethanol (AqE), ethanol (EtOH), ethyl acetate (EtOAc), and hexane (Hex). The concentration of each metabolite is colored based on a normalized scale from minimum−3 (dark green) to a maximum of 3 (dark red).
Figure 5. Heat map of identified metabolites of I. galbana extracted in 5 different solvents—Aqueous (Aq), 50% ethanol (AqE), ethanol (EtOH), ethyl acetate (EtOAc), and hexane (Hex). The concentration of each metabolite is colored based on a normalized scale from minimum−3 (dark green) to a maximum of 3 (dark red).
Marinedrugs 19 00139 g005
Marinedrugs 19 00139 g006 550
Figure 6. Relative quantification based on the mean peak area of 1H-NMR signals of identified metabolites in I. galbana extracted in 5 different solvents: Aqueous (Aq), 50% ethanol (AqE), ethanol (EtOH), ethyl acetate (EtOAc), and hexane (Hex). Data are expressed as a mean of six replicates each of the solvent systems ± standard deviation (SD). Means with different letters indicate significant differences (p < 0.05; n = 6).
Figure 6. Relative quantification based on the mean peak area of 1H-NMR signals of identified metabolites in I. galbana extracted in 5 different solvents: Aqueous (Aq), 50% ethanol (AqE), ethanol (EtOH), ethyl acetate (EtOAc), and hexane (Hex). Data are expressed as a mean of six replicates each of the solvent systems ± standard deviation (SD). Means with different letters indicate significant differences (p < 0.05; n = 6).
Marinedrugs 19 00139 g006
Marinedrugs 19 00139 g007 550
Figure 7. Total ion chromatogram (TIC) of I. galbana EtOAc extract analysed in positive mode (a), negative mode (b), and PDA (photodiode array) (c).
Figure 7. Total ion chromatogram (TIC) of I. galbana EtOAc extract analysed in positive mode (a), negative mode (b), and PDA (photodiode array) (c).
Marinedrugs 19 00139 g007
Marinedrugs 19 00139 g008 550
Figure 8. MS/MS spectrum of fucoxanthin at [M + H]+ m/z 659.4283 with proposed fragmentation patterns.
Figure 8. MS/MS spectrum of fucoxanthin at [M + H]+ m/z 659.4283 with proposed fragmentation patterns.
Marinedrugs 19 00139 g008
Marinedrugs 19 00139 g009 550
Figure 9. MS/MS spectrum of astaxanthin at m/z 657.4136, [M + IsoProp + H]+ and proposed fragmentation patterns.
Figure 9. MS/MS spectrum of astaxanthin at m/z 657.4136, [M + IsoProp + H]+ and proposed fragmentation patterns.
Marinedrugs 19 00139 g009
Marinedrugs 19 00139 g010 550
Figure 10. MS/MS spectrum of chlorophyll a with m/z 893.5519, [M + H]+ and proposed fragmentation patterns.
Figure 10. MS/MS spectrum of chlorophyll a with m/z 893.5519, [M + H]+ and proposed fragmentation patterns.
Marinedrugs 19 00139 g010
Marinedrugs 19 00139 g011 550
Figure 11. MS/MS-based molecular networking from a cluster of glycerophospholipids and identifying 11 new structures as listed in Table 4. Nodes are labeled with parent m/z values, with different size corresponding to precursor intensity. Edges are labeled with cosine scores from 0 to 1.
Figure 11. MS/MS-based molecular networking from a cluster of glycerophospholipids and identifying 11 new structures as listed in Table 4. Nodes are labeled with parent m/z values, with different size corresponding to precursor intensity. Edges are labeled with cosine scores from 0 to 1.
Marinedrugs 19 00139 g011
Marinedrugs 19 00139 g012 550
Figure 12. Represent MS/MS spectra showing common fragments shared among (a) structures 1–7 and (b) structures 8–11 acquired by UPLC–MS/MS in positive mode.
Figure 12. Represent MS/MS spectra showing common fragments shared among (a) structures 1–7 and (b) structures 8–11 acquired by UPLC–MS/MS in positive mode.
Marinedrugs 19 00139 g012
Table
Table 1. Yield for each solvent extract.
Table 1. Yield for each solvent extract.
SolventAqAqEEtOHEtOAcHex
Yield (g)0.068 ± 0.003 d0.042 ± 0.001 ab0.057 ± 0.002 c0.047 ± 0.001 b0.038 ± 0.002 a
Data are expressed as mean ± SD. Values with the different superscript letters are significantly different as determined by Turkey test (p < 0.05).
Table
Table 2. Characteristic 1H NMR signals of metabolites identified in various solvent extracts of I. galbana (s—singlet; d—doublet; t—triplet; dd—doublet of doublets; m—multiplet; J—coupling constant).
Table 2. Characteristic 1H NMR signals of metabolites identified in various solvent extracts of I. galbana (s—singlet; d—doublet; t—triplet; dd—doublet of doublets; m—multiplet; J—coupling constant).
No.Putative MetaboliteδH (ppm), Multiplicity, J (Hz)HSQC (1H-13C)HexEtOAcEtOHAqEAqReference
Amino acids
1.Valine0.99 (d, 6.2)25.1--+++[29]
1.05 (d, 2.8)21.1
2.Isoleucine1.01 (d, 3.3)19.5--+++[29]
0.96 (t, 7.6)-
3.Leucine0.97 (d, 5.0)23.8+++++[29]
4.Alanine3.74 (m)72.8--+++[44]
1.46 (d, 7.2)19.1
5.Threonine3.64 (dd, 5.7, 10.9)66.5----+[44]
1.37 (d, 6.9)-
Carbohydrates/Others
6.Sucrose 5.37 (d, 4.1)---+++[45]
4.22 (d, 2.5)76.5
3.83 (m)-
3.68 (s)-
3.44 (t, 6.7)74.5
7.Glucose3.78 (m)72.8--+++[45]
3.52 (dd, 4.6, 8.4)66.5
8.Choline3.22 (s)57.0--+++[29]
9.Dimethylsulphonio-propionate (DMSP)3.44 (t, 6.7)43.9--+++[29]
2.92 (s)28.3
2.70 (t, 6.7)33.4
10.d-1,4/2,5-cyclohexa-netetrol3.72 (m)73.2--+++[29]
1.83 (m)37.3
Saturated fatty acids
11.Palmitic acid2.35 (t, 7.6) ++++-[45]
1.61 (m)27.9
1.33 (m)33.2
0.88 (t, 6.9)16.4
12.Arachidic acid2.35 (t, 7.6) ++++-[45]
1.60 (m)27.9
1.40 (m)18.2
1.29 (m)32.2
0.88 (t, 6.9)16.4
Unsaturated fatty acids
13.Oleic acid5.34 (m)130.3++++-[45]
2.25 (t, 7.5)37.4
2.02 (m)23.3
1.27 (m)33.2
0.87 (t, 6.9)16.4
14.α-linoleic acid5.34 (m)-++++-[45]
5.32 (m)-
2.80 (m)29.1
2.25 (t, 7.5)37.4
2.02 (m)-
1.27 (m)-
0.96 (t, 7.6)16.6
15.Docosahexaenoic fatty acid chain (DHA)5.30–5.38 (m)130.3++++-HMDB
2.38 (m)36.8
2.06 (m)29.7
0.96 (t, 7.6)16.6
16.Docosapentaenoic fatty acid chain (DPA)5.30–5.38 (m)131.8++++-HMDB
2.33 (m)36.9
1.38 (m)31.3
1.30 (m)33.2
0.87 (t, 6.9)32.0
17.Cholesterol 5.28 (dd, 8.3, 17.6)-++++-HMDB
2.31 (dd, 2.8, 6.2)-
1.61 (m)-
0.82 (m)22.2
0.77 (d, 6.6)25.0
0.65 (m)-
Carotenoids
18.Fucoxanthin6.77 (dd, 4.2, 4.2)--++--[29], HMDB
6.41 (dd, 11.1, 22.2)-
2.59 (d, 18.5)-
2.12 (s)-
2.09 (s)-
1.97 (s)-
1.58 (dd, 7.2, 13.7)-
1.37 (dd, 5.7, 9.7)31.3
1.19 (s)23.1
1.08 (s)34.7
1.01 (s)26.9
19.Astaxanthin6.77 (m)-+++--[45], HMDB
6.61 (d, 3.8)-
6.63 (d, 3.2)-
4.15 (dd)65.2
2.12 (s)-
1.91 (s)13.8
1.19 (s)23.1
Chlorophylls
20.Pheophytin a9.21 (s)107.2+++--[46], PubChem
9.00 (s)99.9
8.60 (s)96.2
Alkenone
21.Polyunsaturated long-chain alkenones (PULCAs)2.12 (s)-+++--[46], PubChem
1.95 (m)35.4
Positive (+) and negative (-) signs denote present and absent, respectively; Ref. = Reference.
Table
Table 3. Putative metabolites identified in EtOAc extract of I. galbana by LCMS/MS in positive and negative mode; n.d. = not detected. (Exp. = experimental; Theo = theoretical).
Table 3. Putative metabolites identified in EtOAc extract of I. galbana by LCMS/MS in positive and negative mode; n.d. = not detected. (Exp. = experimental; Theo = theoretical).
PeaktR (min)Putative MetaboliteExp. Mass (M + H)Exp. Mass (M − H)Theo. Mass (M +/− H)Mass Error (ppm)MS Fragments (ESI+)UV (nm)MAIN CLASS
1.7.00Arachidic acid 311.1688311.3028−430293, 267, 249, 223n.dFatty acid
2.10.11Astaxanthin597.3928 597.3938−1.67597, 579, 279, 215, 109438Carotenoid
3.10.18Halocynthiaxanthin599.4077 599.4095 −3.00389, 233, 147, 109448Carotenoid
4.10.713,6,9,12,15-Octadecapentaenoic acid 273.1861273.18600.36228, 182, 133, 59n.dFatty acid
5.11.40Chlorophyll c2609.1970 609.1983−2.13591, 549, 532450Chlorophyll
6.11.21Stearidonic acid 275.2015275.2017−0.7283, 71, 59n.dFatty acid
7.11.52Chlorophyll c1611.2125 611.2139−2.25593, 551, 534446Chlorophyll
8.11.90Cholesterol387.1795 387.3548−452387, 362, 207n.dSterol
9.12.05Eicosapentaenoic acid (EPA) 301.2173301.21730187, 166, 148n.dFatty acid
10.12.07Fucoxanthin659.4281 659.4306 −3.79581, 411, 355, 199, 109448Carotenoid
11.12.08α-Linolenic acid 277.2172277.2173−0.36194, 92, 87, 59n.dFatty acid
12.12.49Docosahexaenoic acid (DHA) 327.2328327.2330−0.61213, 172, 135, 59n.dFatty acid
13.13.10Pheophorbide a593.2742 593.2758−2.69593, 533, 506, 459408, 536Chlorophyll
14.13.29Docosapentaenoic acid (DPA) 329.2487329.2486−0.30250, 226, 85, 59n.dFatty acid
15.13.37(3S,4R,3′R)-4-Hydroxyalloxanthin581.3976 581.3989−2.23563, 411, 251438Carotenoid
16.13.60Phoenicoxanthin581.3964 581.3989−4.30411, 429452Carotenoid
17.13.82Palmitic acid 255.2329255.2330−0.39246, 94, 81, 76n.dFatty acid
18.13.91Oleic Acid 281.2486281.24860101,98, 87, 64n.dFatty acid
19.15.26Hexacosanedioic acid 427.3771 427.3782 −2.57203, 139, 121, 71n.dFatty acid
20.15.83Diatoxanthin567.4182 567.4197−2.64255, 211, 119, 109460Carotenoid
21.15.92Canthaxanthin565.4023 565.4040−3.00447, 255, 119456Carotenoid
22.16.45Echinenone 551.4232 551.4247−2.72551, 502, 458, 447458Carotenoid
23.17.56PI(16:0/22:4(7Z,10Z,13Z,16Z))887.5657 887.56441.35871, 609, 591n.d Glycerophospholipids
24.17.90PC(15:1(9Z)/22:6(4Z,7Z,10Z,13Z,16Z,19Z))790.5438 790.53817.21628, 610, 356n.dGlycerophospholipids
25.18.68PI-Cer(d14:0/31:0)936.7314 936.72645.33919, 643, 591n.dSphingolipids
26.19.82PS(O-14:0/26:0)834.6580 834.6582−0.23698, 589, 543n.dGlycerophospholipids
27.20.51Chlorophyll a893.5512 893.54269.62893, 615, 539408, 536Chlorophyll
28.21.06Pheophytin a871.5711 871.5732−2.41871, 593, 533408, 536Chlorophyll
29.21.59PI-Cer(d14:0/28:0)930.6123 930.6196−7.84631, 603, 506 Sphingolipids
30.22.80NAPE(18:1(9Z)/16:1(9Z)/18:0)982.7805 982.7834−2.95921, 828, 636 Glycerophospholipids
31.24.63PI-Cer(t18:0/16:0(2OH))814.5496 814.54406.87797, 569, 543 Sphingolipids
32.28.47TG(12:0/16:0/22:5(7Z,10Z,13Z,16Z,19Z))825.6924 825.6967−5.21698, 597, 549 Glycerolipids
Table
Table 4. Putative annotation of metabolites from a cluster of glycerophospholipids identified in EtOAc extract of I. galbana.
Table 4. Putative annotation of metabolites from a cluster of glycerophospholipids identified in EtOAc extract of I. galbana.
No.Parent Mass m/z [M + H]+Mass Error (ppm)Putative AnnotationMolecular FormulaRT (min)
1.530.331313.58LysoPE(22:4(7Z,10Z,13Z,16Z)/0:0)C27H48NO7P9.24
2.788.52867.74PE(22:6(4Z,7Z,10Z,13Z,16Z,19Z)/
18:2(9Z,12Z))		C45H74NO8P14.38
3.609.4493−3.28DG(18:4(6Z,9Z,12Z,15Z)/18:4
(6Z,9Z,12Z,15Z)/0:0)		C39H60O515.00
4.794.57618.43PE(20:2(11Z,14Z)/20:3(8Z,11Z,14Z))C45H80NO8P15.47
5.840.55936.54PE(22:4(7Z,10Z,13Z,16Z)/22:6
(4Z,7Z,10Z,13Z,16Z,19Z))		C49H78NO8P14.74
6.768.56008.07PE(18:0/20:4(5Z,8Z,11Z,14Z))C43H78NO8P15.85
7.740.52878.37PE(14:0/22:4(7Z,10Z,13Z,16Z))C41H74NO8P14.93
8.720.56018.74PE(16:0/18:0)C39H78NO8P16.44
9.746.57578.44PE(18:0/18:1(9Z))C41H80NO8P16.55
10.742.54468.75PE(16:1(9Z)/20:2(11Z,14Z))C41H76NO8P15.58
11.692.5299.39PE(16:0/16:0)C37H74NO8P15.67
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Bustamam, M.S.A.; Pantami, H.A.; Azizan, A.; Shaari, K.; Min, C.C.; Abas, F.; Nagao, N.; Maulidiani, M.; Banerjee, S.; Sulaiman, F.; et al. Complementary Analytical Platforms of NMR Spectroscopy and LCMS Analysis in the Metabolite Profiling of Isochrysis galbana. Mar. Drugs 2021, 19, 139. https://0-doi-org.brum.beds.ac.uk/10.3390/md19030139

AMA Style

Bustamam MSA, Pantami HA, Azizan A, Shaari K, Min CC, Abas F, Nagao N, Maulidiani M, Banerjee S, Sulaiman F, et al. Complementary Analytical Platforms of NMR Spectroscopy and LCMS Analysis in the Metabolite Profiling of Isochrysis galbana. Marine Drugs. 2021; 19(3):139. https://0-doi-org.brum.beds.ac.uk/10.3390/md19030139

Chicago/Turabian Style

Bustamam, Muhammad Safwan Ahamad, Hamza Ahmed Pantami, Awanis Azizan, Khozirah Shaari, Chong Chou Min, Faridah Abas, Norio Nagao, Maulidiani Maulidiani, Sanjoy Banerjee, Fadzil Sulaiman, and et al. 2021. "Complementary Analytical Platforms of NMR Spectroscopy and LCMS Analysis in the Metabolite Profiling of Isochrysis galbana" Marine Drugs 19, no. 3: 139. https://0-doi-org.brum.beds.ac.uk/10.3390/md19030139

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