Mei Sun
Zhibo Yang
Boris Wawrik- 1Department of Chemistry and Biochemistry, University of Oklahoma, Norman, OK, United States
- 2Department of Botany and Microbiology, University of Oklahoma, Norman, OK, United States
Traditional approaches for the assessment of physiological responses of microbes in the environment rely on bulk filtration techniques that obscure differences among populations as well as among individual cells. Here, were report on the development on a novel micro-scale sampling device, referred to as the “Single-probe,” which allows direct extraction of metabolites from living, individual phytoplankton cells for mass spectrometry (MS) analysis. The Single-probe is composed of dual-bore quartz tubing which is pulled using a laser pipette puller and fused to a silica capillary and a nano-ESI. For this study, we applied Single-probe MS technology to the marine dinoflagellate Scrippsiella trochoidea, assaying cells grown under different illumination levels and under nitrogen (N) limiting conditions as a proof of concept for the technology. In both experiments, significant differences in the cellular metabolome of individual cells could readily be identified, though the vast majority of detected metabolites could not be assigned to KEGG pathways. Using the same approach, significant changes in cellular lipid complements were observed, with individual lipids being both up- and down-regulated under light vs. dark conditions. Conversely, lipid content increased across the board under N limitation, consistent with an adjustment of Redfield stoichiometry to reflect higher C:N and C:P ratios. Overall, these data suggest that the Single-probe MS technique has the potential to allow for near in situ metabolomic analysis of individual phytoplankton cells, opening the door to targeted analyses that minimize cell manipulation and sampling artifacts, while preserving metabolic variability at the cellular level.
Introduction
Globally, marine phytoplankton contribute ca. 45 petagrams carbon per annum to net primary production (NPP) (Falkowski et al., 1998). Phytoplankton are thus important drivers of several global biogeochemical cycles, notably of elements that are components of cellular biomass, including carbon (C), nitrogen (N), phosphorus (P), silicate (Si), and others. The abundance of phytoplankton in marine systems is shaped by the availability of nutrients as well as physical oceanic processes. Typically, phytoplankton biomass can be broadly approximated from environmental variables, such as nutrient concentrations, sea surface temperature, and solar irradiance (Antoine et al., 1996; Falkowski et al., 1998). However, at a more granular level, such as the cellular response of individual phytoplankton to dynamic oceanographic conditions, requisite adaptations are often not well-understood. In particular, whether phytoplankton are nutrient limited in the environment has long attracted the attention of oceanographers who are trying to understand the controls on NPP, given the large impact that limitation may have on the structure of marine ecosystems.
Historically, both N and P have been understood to, at times, limit productivity in marine systems (Howarth, 1988), but other nutrient, such as iron (Coale et al., 1998; Landry et al., 2000), also appear to play an important role. Additionally, phytoplankton can be co-limited by more than one nutrient due to their limited absolute abundance (i.e., kinetic limitations), or the acquisition of one nutrient may be dependent on the concentration of another (Saito et al., 2008). Understanding requisite limitations in natural systems is unfortunately not a straightforward matter, given the restricted methodological options available to researchers. For example, nutrient limitation is frequently invoked via nutrient ratios (Liebig’s law of the minimum), yet this ignores that turnover rates can be high in the face of low, but non-limiting ambient concentrations (Beardall et al., 2001; Wawrik et al., 2004). Similarly, bottle incubations (spiking nutrients) are used to infer nutrient limitation via measurements of cellular activity (e.g., carbon fixation), biomass (e.g., chl a), or photosynthetic capacity (Fv/Fm), but interpretation of requisite data can be challenging, given that measurements are made on bulk communities in which individual phytoplankton species may exhibit differential behavior. Changes in elemental uptake ratios (C: P or C: N) have been used to infer nutrient limitation (Rees et al., 1995; Beardall et al., 2001), but such analyses suffer from drawbacks similar to those of bottle incubations, in that measurements are conducted on bulk communities.
An alternate avenue has been the development of specific molecular probes for the expression of marker genes. This approach draws on a long tradition of studies that aim to develop molecular targets for specific microbial nutrient cycling activities in the environment (Scala and Kerkhof, 1998; Allen et al., 2001; Wawrik et al., 2002; and many more). For example, the expression of the global nitrogen regulator NtcA in marine cyanobacteria has been used to assess the nutritional status of natural populations of Synechococcus (Lindell and Post, 2001). Similarly, the expression of the nifH gene is widely distributed in marine systems (Turk et al., 2011), indicator of cellular alleviation of N limitation in some microbial populations. However, while this approach is powerful, it typically requires a fairly good understanding of the underlying molecular mechanisms and genetic diversity of related genes to allow for the derivation of probes or primers. Other studies have applied less targeted, transcriptomic approaches to detect nutrient limitation in marine phytoplankton (Hwang et al., 2008; Cooper et al., 2016; Harke et al., 2017), but analogous studies are more difficult in diverse natural assemblages for which genetic information is not necessarily available. Given this, our understanding of the degree and severity of nutrient limitation in the marine environment, especially at the induvial species and cellular level, remains poorly constrained.
Recently, several approaches have been developed to perform the single cell analysis. These include fluorescence, capillary electrophoresis, and mass spectrometry (MS) (Woods and Ewing, 2003; Woods et al., 2004; Cohen et al., 2008; Urban et al., 2011; Cahill et al., 2015). Here, we report the application of a novel technology that allows for the analysis of the metabolome of single phytoplankton cells to assess their physiological status. The approach utilizes the ‘Single-probe,’ a micro-scale sampling and ionization device, that is coupled to an XYZ-stage to directly insert into single phytoplankton cells to extract cellular metabolites for real-time MS analysis (Pan et al., 2014; Rao et al., 2016b). Our aims were threefold. First, we intended to establish a proof of concept, i.e., demonstrate that metabolome data could be generated for single marine phytoplankton cells. Second, we aimed to demonstrate that physiological changes at the cellular level could be detected via analysis cellular metabolites sampled. Lastly, we aimed to demonstrate that single cell metabolomics can be utilized to assess whether cells are experiencing different illumination levels and nutrient limitation.
Materials and Methods
Cultures
Non-axenic Scrippsiella trochoidea CCMP 3099 was originally obtained from the National Center for Marine Algae and Microbiota (Provasoil-Guillard NCMA, Boothbay Harbor, ME, United States). For maintenance, cultures were grown in L1 seawater media (Guillard and Ryther, 1962; Guillard et al., 1973; Guillard and Hargraves, 1993). Medium was prepared from natural seawater collected near Key West (salinity of 33), which was aged for at least 6 month in the dark and autoclaved. Maintenance and experimental cultures were grown in a light/dark incubator at 23–24°C and 30–40 μmol quanta⋅m-2 s-1 light under a 12-h light:12-h dark cycle.
Experimental Culture Setup
For the light/dark comparisons, cultures were grown under replete conditions in full L1 media containing 880 μM NaNO3 and 36 μM NaH2PO4 (N/P ratio 24:1). Experimental cultures were started as a 1:10 inoculum from exponentially growing cultures into 1 L of media in 2.5 L Pyrex Fernbach flasks without shaking and monitored daily via cell counts and chlorophyll measurements. Cell counts were conducted by addition of 1% Lugol’s iodine and direct counting of cells in 96-well microtiter plates using a dissection microscope. Dilutions were made as necessary and at least 10 wells containing 100 μL diluted culture were counted to average cell counts. Chlorophyll a was measured via fluorometry (Welschmeyer, 1994) by filtering 5 mL of culture in triplicate onto GF/F filters, over-night extraction with methanol, and quantification using a Turner Trilogy Laboratory Fluorometer. Cultures were grown into late-log phase (data not shown) and then sampled 3 h before and 3 h after the light was turned on. Cultures were sub-sampled for MS analysis, making sure to keep ‘dark’ sample exposure to light to a minimum by wrapping sampling tubes in aluminum foil.
N-deplete cultures were generated by first growing cells on L1 media with an N:P ratio of 2.4:1 (88 mM nitrate: 36 mM phosphate) analogous what has been previously described (Harke et al., 2017). This lower ratio stoichiometrically limited cultures in nitrogen and at least three transfers were performed to ensure no carryover from higher nutrient full L1 medium. Cultures were monitored daily via cell counts (see above), chlorophyll a quantification (see above), and quantification of nitrate/nitrite via a Vanadium reduction method (Miranda et al., 2001). Parallel cultures were set up in which one culture was allowed to run out of nitrogenous nutrients (N-deplete), while the culture (control) was fed additional nitrate every second day to bring total nitrate/nitrite concentrations back to starting levels (Supplementary Figure S1). Once nitrate/nitrite levels dropped below the limit of detection (∼1–2 μM) in the N-deplete culture, cultures were grown for an additional 24 h before sampling to ensure that N-depletion was complete. Both the N-deplete and replete cultures were then sampled for MS analysis of single cells.
Mass Spectrometry
Individual cells of S. trochoidea were analyzed via the ‘Single-probe’ MS techniques. Detailed fabrication protocols of the Single-probe have been previously described (Pan et al., 2014; Rao et al., 2015; Sun et al., 2017). Briefly, a Single-probe (Figure 1A) has three components: a dual-bore quartz tubing [outer diameter (OD) 500 μm; inner diameter (ID) 127 μm, Friedrich & Dimmock, Inc., Millville, NJ, United States] pulled using a laser pipette puller (P-2000 micropipette puller, Sutter Instrument, Novato, CA, United States), a fused silica capillary (OD 105 μm; ID 40 μm, Polymicro Technologies, Phoenix, AZ, United States), and a nano-ESI emitter made from the same type of fused silica capillary. A Single-probe is fabricated by embedding a fused silica capillary and a nano-ESI emitter into both of the channels of the laser-pulled dual-bore quartz needle.
FIGURE 1. Experimental setup to measure single Scrippsiella trochoidea cells using the “Single-probe” MS techniques. (A) Photograph of the Single-probe device with its different components labeled; (B) image from microscope-linked camera used to target single S. trochoidea cell with the Single-probe; (C) setup used to manipulate the Single-probe MS device with components labeled.
For the analysis, cells were deposited onto 0.2 μm polycarbonate membranes by gentle filtration, and the cells were rinsed with phosphate buffered saline (PBS; 137 mM NaCl, 2.7 mM KCl, 4.3 mM Na2HPO4, 1.47 mM KH2PO4; pH of 7.4) to remove culture medium. Filters were then placed on a home-built XYZ-translation stage system and spatial motion was controlled by a custom designed LabView software package (Lanekoff et al., 2012). The Single-probe tip (<10 μm) was then precisely insert into single S. trochoidea cells (typically ∼20–30 μm cellular diameter) using a microscope as a guide (Figure 1B). During the experiment, a syringe (250 μl; Hamilton, Co., Reno, NV, United States) was used to continuously provide the sampling solvent (acetonitrile; Sigma-Aldrich, St. Louis, MO, United States), and a liquid junction formed at the Single-probe tip to perform highly efficient extraction of cellular contents. The analytes were withdrawn by capillary action toward the nano-ESI emitter, and ionized for analysis using a Thermo LTQ Orbitrap XL mass spectrometer (Thermo Scientific, Waltham, MA, United States) (Figure 1C). Mass analyze parameters were as follows: mass resolution 60,000, +4 kV ionization voltage at positive ion mode (0.05–0.07 μA of ion current), 1 microscan, 100 ms max injection time, and automatic gain control on.
Data Analysis
The Thermo Xcalibur Qual Browser (Thermo Scientific, Waltham, MA, United States) was used to export MS data (m/z values with relative intensities) as tab-delimited data files. As a conservative approach, only relatively abundant peaks with ion intensities > 103 were exported. This approach excluded 6% of low signal peaks as background while retaining 94% of total signal intensity. The relative ion intensities were normalized to the total ion current to minimize the influence induced by fluctuations of ion signals during experiments. The Geena 2 online software tool1 was then used for peak alignment (Romano et al., 2016), and the aligned m/z values were then used for comparisons. Parameters used in Geena 2 include analysis range (from 100 to 1500 m/z), maximum number of isotopic replicas (5), maximum delta between isotopic peaks (0.05 Da), and maximum delta for aligning replicates (0.01 Da). Metaboanalyst 3.02 was used to conduct statistical data analysis, including PLS-DA (partial least squares discriminant analysis) and t-tests (Xia et al., 2015; Xia and Wishart, 2016). PLS-DA was used to visualize differences in chemical composition profiles among treatment groups, while t-tests were applied to extract molecular peaks with significant abundance changes (p < 0.05). Finally, the online database METLIN3 was used to tentatively label all ions of interest (Smith et al., 2005; Guijas et al., 2018), to perform hierarchical clustering, and to generate heat maps. Lastly, Pathos4 was used to attempt identification of significantly regulated metabolic pathways by considering all KEGG maps in all organisms (Leader et al., 2011).
Results
The main aims of this study were to develop a single-cell-based metabolomic methodology that could be applied to individual algal cells and to demonstrate that this technology can detect physiological responses to environmental stimuli. The described setup allowed us to sample individual cells by targeting them with the tip of the Single-probe (Figure 1B). Our ‘proof of concept’ experiment involved a comparison of cells in the context of diurnal illumination changes, which are known to induce significant changes in algal cellular metabolomes (Vidoudez and Pohnert, 2012).
Using the Single-probe MS technique, distinct and clear differences were observed in the S. trochoidea metabolome under different light levels. A total of 1,085 and 1,103 peaks were detected under light and dark conditions, respectively. We have tentatively labeled 581 species (i.e., 558 metabolites and 23 peptides; Supplementary Table S1). All MS data are archived using Zenodo (DOI: 10.5281/zenodo.1188486), and the descriptions of MS files were summarized in Supplementary Table S2, and a single cell mass spectrum (Supplementary Figure S2) is shown as an example. Of these detected species, 306 and 321 were differentially abundant among treatments (t-test; p < 0.05) under light and dark conditions, respectively. Results obtained from PLS-DA (Figure 2A) revealed that metabolic features formed distinct clusters (15 cells in each group), and this difference was highly significant (p = 5 × 10-4; permutation test in Metaboanalyst 3.0). Cross-validation of the PLS-DA model was conducted, and our results (Q2 > 0.5) indicate this model provided a good predictability without overfitting (Supplementary Table S3) (Worley and Powers, 2013; Triba et al., 2015). To investigate which metabolic pathways may have been significantly impacted by the difference in illumination levels, all ions with significant differences (p < 0.05) were selected and tentatively labeled by searching their m/z values in Pathos, considering all KEGG metabolites, and retaining pathways for which at least two hits were observed (Table 1). A criterion of two hits to an individual pathway was used, because the reliability of assignments for metabolites based on m/z values alone is limited. Most strikingly, only 12 (5%) and 16 (19%) of differentially abundant metabolites could be assigned to KEGG pathways for light and dark conditions respectively, indicating that the majority of the metabolic response to light changes in S. trochoidea is not captured in KEGG metabolic maps. With respect to pathways with at least two hits (Table 1), light favors biosynthesis of molecules potentially linked to the production of 12-, 14-, and 16-membered macrolides. In particular, several hits in the pathway for Avermectin were observed (Table 1A). Under dark conditions, identified pathways included those for polyketides, porphyrin, chlorophylls, terpenoids, and limonene (Table 1B). With respect to the four hits in porphyrin and chlorophyll metabolism, all were linked to the production of phycobillins.
FIGURE 2. Partial least aquares discriminant analysis (PLS-DA) of MS data. All detected metabolites were analyzed, visualizing the overall effect on the metabolome of single S. trochoidea cells by (A) light vs. dark conditions [results are reported from 15 replicates (n = 15) in each group], and (B) N-limited vs. replete conditions (n = 10 for each group).
TABLE 1. Pathways containing more than one metabolite with significantly different abundance under (A) illuminated conditions compared to cultures during (B) dark condition.
Given the successful proof of concept for light/dark conditions, we aimed to investigate whether single cell metabolomic analysis might be utilized to investigate the nutritional status of individual phytoplankton cells. As shown in Figure 2B, PLS-DA results indicated that single cellular metabolomes of S. trochoidea under N-limited condition were clearly different from those in nitrogen replete cells (n = 10). As with the light–dark treatments, N limitation induced a highly significant response (p = 1 × 10-3 and Q2> 0.5; Supplementary Table S3). KEGG pathway analysis was attempted, but yielded no metabolites in the significantly upregulated pool that could be matched to KEGG. Similarly, only 8 (<4.5%) metabolites from the down-regulated pool could be mapped with KEGG and no pathways contained more than a single hit. Little information could therefore be gleaned about requisite physiological responses via KEGG analysis, and an alternative approach to analysis was therefore taken. First, we hypothesized that N-limitation should be reflected in the C/N and N/P rations of the cellular metabolome, as cells might physiologically adjust to environmental conductions by choosing cellular metabolites with lower N content (Beardall et al., 2001; Geider and La Roche, 2002). Second, metabolites were analyzed for the putatively detected N-containing lipids (e.g., phosphoethanolamine), given that physiological responses are often manifest in the lipid pool.
C: N and C: P ratios were calculated by summing the product of the number of C, N, P atoms in each metabolite with its relative abundance across all metabolites (Figure 3). This analysis reveals that N-limited S. trochoidea has significantly higher C:N ratios in regulated metabolites as compared to replete conditions (p = 4 × 10-6) (Figure 3A). Similarly, C:P ratios were affected and were significantly higher under N-limitation (p < 0.044) (Figure 3B). Concurrently, the abundances of at least some cellular lipids were significantly affected. When all significantly regulated metabolites putatively identified as lipids are considered (Figure 4), it appears that the availability of light correlates with both up- and down-regulation of specific lipid complements (Figure 4A). Under N-limitation, however, a significant decrease of lipid abundances is observed (p < 0.05) (Figure 4B).
FIGURE 3. Elemental ratios of significantly regulated metabolites. (A) Carbon to Nitrogen ratio; (B) Carbon to Phosphorus ratio.
FIGURE 4. Heat maps generated from hierarchically clustering summarizing the cellular lipids measured from single S. trochoidea cell under different light and nutrient conditions. Only lipids for which abundances were significantly different among treatment pairs are shown. Shown are differences for (A) light vs. dark conditions and (B) replete vs. nitrogen (N) limiting conditions. Red indicates elevated and green indicates decreased signal with respect to the mean signal observed for all tested cells. PA, phosphatic acid; PE, phosphoethanolamine; PG, phosphatidylglycerol; PS, phosphatidylserine; PI, phosphatidylinositol; PC, phosphatidylcholine; MGs, monoglycerides; DGs, diglycerides; TG, triglyceride.
Discussion
Historically, the majority of oceanographic research has targeted natural populations of phytoplankton via bulk filtration techniques (e.g., filtration onto GF/F filters) to assess physiological responses to environmental factors such as light or nutrient limitation. While much has been learned using bulk filtration, requisite approaches suffer from important limitations. Most notably, individual populations or different species are not adequately resolved in this manner. It is also now well-appreciated that cellular functions, such as gene expression, proliferation, or programmed cell death, are subject to significant stochasticity, leading to high cellular chemical and phenotypic diversity at the single cell level, potentially obscuring some important patterns (Fagerer et al., 2013; Comi et al., 2017; Guillaume-Gentil et al., 2017). Single-cell analysis is therefore an attractive methodological choice when studying rare types of cells (cells available are inadequate for bulk analysis) or cells in heterogenous populations, where cell separation or sorting are impractical. Here we present the first report of a single cell-based metabolomic technology that allows for the analysis of intracellular intermediates of individual phytoplankton cells. The Single-probe directly collects cellular contents of living cells, and it does this without significant sample preparation steps (e.g., filtration or solvent extraction), thereby allowing for real-time and targeted analysis that minimizes sampling artifacts.
We note that bulk analyses [e.g., liquid chromatography (LC)/MS, gas chromatography (GC)/MS), or direct-injection MS without separation] were not conducted in this study. It is likely that the number and the types of metabolites detected here differs from those that might have been observed in bulk measurements. Traditional LC/MS and GC/MS techniques have certainly been used in the analysis of marine algae (Barofsky et al., 2009; Burriesci et al., 2012), and their value is clear in requisite studies. However, a direct comparison was not conducted, because we do not view the use of the Single-probe as a replacement for more traditional bulk biomass approaches. Rather, single-cell analysis can serve as a complementary method that allows exploration of inherent cell-to-cell variability in complex and heterogeneous systems that may not be resolved using more traditional approaches.
Given the tip size of the Single-probe, which is ca. 10 μm in diameter, analysis is currently limited to larger phytoplankton and protists. However, the currently applicable size range includes many important bloom forming algae, including toxin producing genera such as Karenia (20–40 μm) or Pseudo-nitzschia (40–175 μm). The Single-probe might therefore offer unique opportunities to help understand the biological forces that shape the success of these organisms in an ecosystem by revealing their metabolomic responses to changing environmental conditions or experimental treatments at the cellular level. We note that the expensive and bulky configuration of equipment described here still precludes easy field deployment, making the Single-probe MS technique most useful under conditions where direct access to the lab is available, or when experimental cultures are assayed. Further development in the miniaturization of high-resolution mass spectrometer or sample preservation maybe be helpful in this regard. In addition, due to the limited amount of cytoplasm found in an individual cell and limitation on the achievable sensitivity of mass spectrometers, most single cell MS studies to date are primarily focused on the analysis of relatively small molecules such as metabolites and peptides (Rubakhin et al., 2011). The detection of larger molecules, such as proteins presenting as relatively lower abundances, at single cell level remains very challenging.
Despite these limitations, single-cell MS technique holds great promise for environmental research. The Single-probe MS technique has been successfully used to study live single cancer cells (Pan et al., 2014, 2016), to map biomolecules on animal tissues with high spatial resolutions (Rao et al., 2015, 2016a,b), and to analyze the extracellular metabolites inside spheroids (Sun et al., 2017). With respect to the ability to detect a large range of cellular metabolites from single cells (i.e., from only a few pico-litters of cytoplasm sampled from a cell with a diameter of approximately10 μm), the Single-probe MS setup provides excellent detection sensitivity. We note, however, that single cell MS measurements are not strictly repeatable (cellular contents are consumed in each measurement), and that it is impractical to evaluate reproducibility in this regard. Multiple cells must therefore be measured (e.g., n = 10–15 in each group in the current study) to normalized intensities for statistical data analysis and to minimize the influence of uncertainties such as fluctuation of ion signals, minor changes of experimental tuning conditions, instrument noise, and variances associated with the matrix effect.
Here, were extend these advances by demonstrating that the physiological status of phytoplankton with respect to light and nutrients can be assessed. The availability of light was correlated with a highly significant response in the metabolite profile of S. trochoidea cells (Figure 2). As might be expected, metabolites related to porphyrin and chlorophyll pathways stood out. In particular, the detected metabolites are related to phytochrome metabolism, which is consistent with the notion that phytochromes are signal-transducing photoreceptors (Smith, 2000). Beyond phototaxis and a handful of secondary metabolite intermediates, however, few metabolites could be assigned to KEGG pathways. This, perhaps, speaks to our limited understanding of metabolism in Dinoflagellate algae. Transcriptomic analysis of S. trochoidea (same strain used here) indicated that this species makes perhaps in excess of 105 transcripts, the majority of which could not be annotated (Cooper et al., 2016).
Cellular lipids have been used to parse physiological responses of living cells under different illumination conditions (Chen et al., 2015; Wacker et al., 2016). For example, under dark conditions, mRNA for some genes involved in lipid biosynthesis have been observed to be elevated in Chlorella (Chen et al., 2015). Differential responses of lipid abundance have also been observed in the marine dinoflagellate Prorocentrum minimum under dark vs. light conditions (Manoharan et al., 1999), and light-induced significant changes in the fatty acid profiles have been reported for freshwater diatom, chrysophyte, cryptophyte, and zygnematophyte algae (Wacker et al., 2016). Differences in requisite light-dependent cellular lipid profiles are thought to be related to alterations of energy storage and the compositions of chloroplast membranes (Guihéneuf et al., 2009; Trentacoste et al., 2013; d’Ippolito et al., 2015). Consistently, our experiments indicated both up- and down-regulation of lipids in response to light in S. trochoidea. The pattern of changes in lipid composition in response to N-limitation was, however, quite different from that observed under light limitation. Most notably, all differentially regulated lipids were more abundant under the nutrient limiting conditions. This is potentially the result of limitation induced imbalances in cellular elemental composition. Under N-limitation almost all parts of central metabolism are impacted (and potentially slowed due to limiting resources), yet lipids, for the most part, do not contain N atoms and their biosynthesis might thereby proceed so long photosynthesis can proceed. Whether lipid accumulation in Scrippsiella is adaptive, however, by possibly allowing energy storage to gain advantages when nutrients are more readily available, or is simply a consequence of the onset of senescence remains to be investigated. Changes in the C:N and C:P ratios do appear to be reflected in the metabolome overall (Figure 3). This observation is consistent with the notion that elemental ratios have considerable plasticity based on nutrient availability (Geider and La Roche, 2002). Previous studies have suggested that nutrient deficiency can cause the accumulation of lipids, such as triglyceride (Takagi et al., 2000; Xin et al., 2010; El-Kassas, 2013), and that phytoplankton funnel excess NADPH to the biosynthesis of triglyceride and fatty acids to regenerate NADP+ (Thompson, 1996; Hu et al., 2008; El-Kassas, 2013). Correspondingly, we observed that a number of triglycerides significantly increased in the N-limited group (Figure 3B). We also note that C:N and C:P ratios reported here deviate from Redfield’s expectations (C:N:P ∼ 106:16:9) (Figure 3A). However, these observations are consistent with our approach. The elemental Redfield ratios are based on all cellular components, yet data analyzed here only covers a subset of all cellular constituents. For example, the Single-probe MS techniques cannot currently measure large biomolecules such as proteins and nucleic acids which are N-rich. The absence of nucleic acids in calculations likely leads to high C:P ratio estimates (Maske, 1982; Laws et al., 1983; Geider and La Roche, 2002). Values reported here are therefore not inconsistent with generally expected elemental ratios for phytoplankton biomass.
In summary, we report the development of a methodology for single-cell metabolomics of small protists such as marine dinoflagellate algae. Using the Single-probe MS technology we can efficiently monitor the cellular physiological responses of phytoplankton under different illumination and nutrient conditions. This offers the opportunity for real-time analysis of natural populations and has the potential ability to provide information about the dynamic metabolic response of individual cells to environment stimulation.
Author Contributions
BW provided funding, conducted culture work, and contributed to data analysis and writing. ZY provided funding, was involved in MS work, and contributed to data analysis and writing. MS conducted much of the work and contributed to data analysis and writing.
Funding
Support for this work came from United States National Science Foundation (Grant No. OCE-1634630) from the Biological Oceanography program.
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpls.2018.00571/full#supplementary-material
Footnotes
- ^ http://bioinformatics.hsanmartino.it/geena2/
- ^ http://www.metaboanalyst.ca/
- ^ https://metlin.scripps.edu/
- ^ http://motif.gla.ac.uk/Pathos/
References
Allen, A. E., Booth, M. G., Frischer, M. E., Verity, P. G., Zehr, J. P., and Zani, S. (2001). Diversity and detection of nitrate assimilation genes in marine bacteria. Appl. Environ. Microbiol. 67, 5343–5348. doi: 10.1128/AEM.67.11.5343-5348.2001
Antoine, D., Andre, J. M., and Morel, A. (1996). Oceanic primary production.2. Estimation at global scale from satellite (coastal zone color scanner) chlorophyll. Global Biogeochem. Cycles 10, 57–69. doi: 10.1029/95gb02832
Barofsky, A., Simonelli, P., Vidoudez, C., Troedsson, C., Nejstgaard, J. C., Jakobsen, H. H., et al. (2009). Growth phase of the diatom Skeletonema marinoi influences the metabolic profile of the cells and the selective feeding of the copepod Calanus spp. J. Plankton Res. 32, 263–272. doi: 10.1093/plankt/fbp121
Beardall, J., Young, E., and Roberts, S. (2001). Approaches for determining phytoplankton nutrient limitation. Aquat. Sci. 63, 44–69. doi: 10.1007/Pl00001344
Burriesci, M. S., Raab, T. K., and Pringle, J. R. (2012). Evidence that glucose is the major transferred metabolite in dinoflagellate–cnidarian symbiosis. J. Exp. Biol. 215, 3467–3477. doi: 10.1242/jeb.070946
Cahill, J. F., Darlington, T. K., Fitzgerald, C., Schoepp, N. G., Beld, J., Burkart, M. D., et al. (2015). Online analysis of single cyanobacteria and algae cells under nitrogen-limited conditions using aerosol time-of-flight mass spectrometry. Anal. Chem. 87, 8039–8046. doi: 10.1021/acs.analchem.5b02326
Chen, T., Liu, J., Guo, B., Ma, X., Sun, P., Liu, B., et al. (2015). Light attenuates lipid accumulation while enhancing cell proliferation and starch synthesis in the glucose-fed oleaginous microalga Chlorella zofingiensis. Sci. Rep. 5:14936. doi: 10.1038/srep14936
Coale, K. H., Johnson, K. S., Fitzwater, S. E., Blain, S. P. G., Stanton, T. P., and Coley, T. L. (1998). IronEx-I, an in situ iron-enrichment experiment: experimental design, implementation and results. Deep Sea Res. Part 2 Top. Stud. Oceanogr. 45, 919–945. doi: 10.1016/S0967-0645(98)00019-8
Cohen, D., Dickerson, J. A., Whitmore, C. D., Turner, E. H., Palcic, M. M., Hindsgaul, O., et al. (2008). Chemical cytometry: fluorescence-based single-cell analysis. Annu. Rev. Anal. Chem. 1, 165–190. doi: 10.1146/annurev.anchem.1.031207.113104
Comi, T. J., Do, T. D., Rubakhin, S. S., and Sweedler, J. V. (2017). Categorizing cells on the basis of their chemical profiles: progress in single-cell mass spectrometry. J. Am. Chem. Soc. 139, 3920–3929. doi: 10.1021/jacs.6b12822
Cooper, J. T., Sinclair, G. A., and Wawrik, B. (2016). Transcriptome analysis of Scrippsiella trochoidea CCMP 3099 reveals physiological changes related to nitrate depletion. Front. Microbiol. 7:639. doi: 10.3389/fmicb.2016.00639
d’Ippolito, G., Sardo, A., Paris, D., Vella, F. M., Adelfi, M. G., Botte, P., et al. (2015). Potential of lipid metabolism in marine diatoms for biofuel production. Biotechnol. Biofuels 8:28. doi: 10.1186/s13068-015-0212-4
El-Kassas, H. Y. (2013). Growth and fatty acid profile of the marine microalga Picochlorum sp. grown under nutrient stress conditions. Egypt. J. Aquat. Res. 39, 233–239. doi: 10.1016/j.ejar.2013.12.007
Fagerer, S. R., Schmid, T., Ibáñez, A. J., Pabst, M., Steinhoff, R., Jefimovs, K., et al. (2013). Analysis of single algal cells by combining mass spectrometry with Raman and fluorescence mapping. Analyst 138, 6732–6736. doi: 10.1039/c3an01135f
Falkowski, P. G., Barber, R. T., and Smetacek, V. V. (1998). Biogeochemical controls and feedbacks on ocean primary production. Science 281, 200–207. doi: 10.1126/science.281.5374.200
Geider, R. J., and La Roche, J. (2002). Redfield revisited: variability of C [ratio] N [ratio] P in marine microalgae and its biochemical basis. Eur. J. Phycol. 37, 1–17. doi: 10.1017/S0967026201003456
Guihéneuf, F., Mimouni, V., Ulmann, L., and Tremblin, G. (2009). Combined effects of irradiance level and carbon source on fatty acid and lipid class composition in the microalga Pavlova lutheri commonly used in mariculture. J. Exp. Mar. Biol. Ecol. 369, 136–143. doi: 10.1016/j.jembe.2008.11.009
Guijas, C., Montenegro-Burke, J. R., Domingo-Almenara, X., Palermo, A., Warth, B., Hermann, G., et al. (2018). METLIN: a technology platform for identifying knowns and unknowns. Anal. Chem. 90, 3156–3164. doi: 10.1021/acs.analchem.7b04424
Guillard, R. R. L., and Hargraves, P. E. (1993). Stichochrysis immobilis is a diatom, not a chrysophyte. Phycologia 32, 234–236.
Guillard, R. R. L., Kilham, P., and Jackson, T. A. (1973). Kinetics of silicon-limited growth in the marine diatom Thalassiosira pseudonana Hasle and Heimdal (= Cyclotella nana Hustedt). J. Phycol. 9, 233–237. doi: 10.1111/j.1529-8817.1973.tb04086.x
Guillard, R. R. L., and Ryther, J. H. (1962). Studies of marine planktonic diatoms: I. Cyclotella nana Hustedt, and Detonula confervacea (Cleve) Gran. Can. J. Microbiol. 8, 229–239. doi: 10.1139/m62-029
Guillaume-Gentil, O., Rey, T., Kiefer, P., Ibáñez, A. J., Steinhoff, R., Brönnimann, R., et al. (2017). Single-cell mass spectrometry of metabolites extracted from live cells by fluidic force microscopy. Anal. Chem. 89, 5017–5023. doi: 10.1021/acs.analchem.7b00367
Harke, M. J., Juhl, A. R., Haley, S. T., Alexander, H., and Dyhrman, S. T. (2017). Conserved transcriptional responses to nutrient stress in bloom-forming algae. Front. Microbiol. 8:1279. doi: 10.3389/fmicb.2017.01279
Howarth, R. W. (1988). Nutrient limitation of net primary production in marine ecosystems. Annu. Rev. Ecol. Syst. 19, 89–110. doi: 10.1146/annurev.es.19.110188.000513
Hu, Q., Sommerfeld, M., Jarvis, E., Ghirardi, M., Posewitz, M., Seibert, M., et al. (2008). Microalgal triacylglycerols as feedstocks for biofuel production: perspectives and advances. Plant J. 54, 621–639. doi: 10.1111/j.1365-313X.2008.03492.x
Hwang, Y. S., Jung, G., and Jin, E. (2008). Transcriptome analysis of acclimatory responses to thermal stress in Antarctic algae. Biochem. Biophys. Res. Commun. 367, 635–641. doi: 10.1016/j.bbrc.2007.12.176
Landry, M. R., Constantinou, J., Latasa, M., Brown, S. L., Bidigare, R. R., and Ondrusek, M. E. (2000). Biological response to iron fertilization in the eastern equatorial Pacific (IronEx II). III. Dynamics of phytoplankton growth and microzooplankton grazing. Mar. Ecol. Prog. Ser. 201, 57–72. doi: 10.3354/meps201057
Lanekoff, I., Heath, B. S., Liyu, A., Thomas, M., Carson, J. P., and Laskin, J. (2012). Automated platform for high-resolution tissue imaging using nanospray desorption electrospray ionization mass spectrometry. Anal. Chem. 84, 8351–8356. doi: 10.1021/ac301909a
Laws, E. A., Redalje, D. G., Karl, D. M., and Chalup, M. S. (1983). A theoretical and experimental examination of the predictions of two recent models of phytoplankton growth. J. Theor. Biol. 105, 469–491. doi: 10.1016/0022-5193(83)90188-1
Leader, D. P., Burgess, K., Creek, D., and Barrett, M. P. (2011). Pathos: a web facility that uses metabolic maps to display experimental changes in metabolites identified by mass spectrometry. Rapid Commun. Mass Spectrom. 25, 3422–3426. doi: 10.1002/rcm.5245
Lindell, D., and Post, A. F. (2001). Ecological aspects of ntcA gene expression and its use as an indicator of the nitrogen status of marine Synechococcus spp. Appl. Environ. Microbiol. 67, 3340–3349. doi: 10.1128/AEM.67.8.3340-3349.2001
Manoharan, K., Lee, T. K., Cha, J. M., Kim, J. H., Lee, W. S., Chang, M., et al. (1999). Acclimation of Prorocentrum minimum (Dinophyceae) to prolonged darkness by use of an alternative carbon source from triacylglycerides and galactolipids. J. Phycol. 35, 287–292. doi: 10.1046/j.1529-8817.1999.3520287.x
Maske, H. (1982). Ammonium-limited continuous culures of Skeletonema costatum in steady and transitional state: experimental results and model simulations. J. Mar. Biol. Assoc. U. K. 62, 919–943. doi: 10.1017/S0025315400070430
Miranda, K. M., Espey, M. G., and Wink, D. A. (2001). A rapid, simple spectrophotometric method for simultaneous detection of nitrate and nitrite. Nitric Oxide 5, 62–71. doi: 10.1006/niox.2000.0319
Pan, N., Rao, W., Kothapalli, N. R., Liu, R., Burgett, A. W., and Yang, Z. (2014). The single-probe: a miniaturized multifunctional device for single cell mass spectrometry analysis. Anal. Chem. 86, 9376–9380. doi: 10.1021/ac5029038
Pan, N., Rao, W., Standke, S. J., and Yang, Z. (2016). Using dicationic ion-pairing compounds to enhance the single cell mass spectrometry analysis using the single-probe: a microscale sampling and ionization device. Anal. Chem. 88, 6812–6819. doi: 10.1021/acs.analchem.6b01284
Rao, W., Pan, N., Tian, X., and Yang, Z. (2016a). High-resolution ambient MS imaging of negative ions in positive ion mode: using dicationic reagents with the single-probe. J. Am. Soc. Mass Spectrom. 27, 124–134.
Rao, W., Pan, N., and Yang, Z. (2015). High resolution tissue imaging using the single-probe mass spectrometry under ambient conditions. J. Am. Soc. Mass Spectrom. 26, 986–993. doi: 10.1007/s13361-015-1091-4
Rao, W., Pan, N., and Yang, Z. (2016b). Applications of the single-probe: mass spectrometry imaging and single cell analysis under ambient conditions. J. Vis. Exp. 112:e53911.
Rees, A. P., Owens, N. J. P., Heath, M. R., Plummer, D. H., and Bellerby, R. S. (1995). Seasonal nitrogen assimilation and carbon fixation in a fjordic sea loch. J. Plankton Res. 17, 1307–1324. doi: 10.1093/plankt/17.6.1307
Romano, P., Profumo, A., Rocco, M., Mangerini, R., Ferri, F., and Facchiano, A. (2016). Geena 2, improved automated analysis of MALDI/TOF mass spectra. BMC Bioinformatics 17(Suppl. 4):61. doi: 10.1186/s12859-016-0911-2
Rubakhin, S. S., Romanova, E. V., Nemes, P., and Sweedler, J. V. (2011). Profiling metabolites and peptides in single cells. Nat. Methods 8:S20–S29. doi: 10.1038/nmeth.1549
Saito, M. A., Goepfert, T. J., and Ritt, J. T. (2008). Some thoughts on the concept of colimitation: three definitions and the importance of bioavailability. Limnol. Oceanogr. 53, 276–290. doi: 10.4319/lo.2008.53.1.0276
Scala, D. J., and Kerkhof, L. J. (1998). Nitrous oxide reductase (nosZ) gene-specific PCR primers for detection of denitrifiers and three nosZ genes from marine sediments. FEMS Microbiol. Lett. 162, 61–68. doi: 10.1111/j.1574-6968.1998.tb12979.x
Smith, C. A., O’Maille, G., Want, E. J., Qin, C., Trauger, S. A., Brandon, T. R., et al. (2005). METLIN: a metabolite mass spectral database. Drug Monit. 27, 747–751. doi: 10.1097/01.ftd.0000179845.53213.39
Smith, H. (2000). Phytochromes and light signal perception by plants–an emerging synthesis. Nature 407, 585–591. doi: 10.1038/35036500
Sun, M., Tian, X., and Yang, Z. (2017). Microscale mass spectrometry analysis of extracellular metabolites in live multicellular tumor spheroids. Anal. Chem. 89, 9069–9076. doi: 10.1021/acs.analchem.7b01746
Takagi, M., Watanabe, K., Yamaberi, K., and Yoshida, T. (2000). Limited feeding of potassium nitrate for intracellular lipid and triglyceride accumulation of Nannochloris sp. UTEX LB1999. Appl. Microbiol. Biotechnol. 54, 112–117. doi: 10.1007/s002530000333
Thompson, G. A. (1996). Lipids and membrane function in green algae. Biochim. Biophys. Acta 1302, 17–45. doi: 10.1016/0005-2760(96)00045-8
Trentacoste, E. M., Shrestha, R. P., Smith, S. R., Glé, C., Hartmann, A. C., Hildebrand, M., et al. (2013). Metabolic engineering of lipid catabolism increases microalgal lipid accumulation without compromising growth. Proc. Natl. Acad. Sci. U.S.A. 110, 19748–19753. doi: 10.1073/pnas.1309299110
Triba, M. N., Le Moyec, L., Amathieu, R., Goossens, C., Bouchemal, N., Nahon, P., et al. (2015). PLS/OPLS models in metabolomics: the impact of permutation of dataset rows on the K-fold cross-validation quality parameters. Mol. Biosyst. 11, 13–19. doi: 10.1039/c4mb00414k
Turk, K. A., Rees, A. P., Zehr, J. P., Pereira, N., Swift, P., Shelley, R., et al. (2011). Nitrogen fixation and nitrogenase (nifH) expression in tropical waters of the eastern North Atlantic. ISME J. 5, 1201–1212. doi: 10.1038/ismej.2010.205
Urban, P. L., Schmid, T., Amantonico, A., and Zenobi, R. (2011). Multidimensional analysis of single algal cells by integrating microspectroscopy with mass spectrometry. Anal. Chem. 83, 1843–1849. doi: 10.1021/ac102702m
Vidoudez, C., and Pohnert, G. (2012). Comparative metabolomics of the diatom Skeletonema marinoi in different growth phases. Metabolomics 8, 654–669. doi: 10.1007/s11306-011-0356-6
Wacker, A., Piepho, M., Harwood, J. L., Guschina, I. A., and Arts, M. T. (2016). Light-induced changes in fatty acid profiles of specific lipid classes in several freshwater phytoplankton species. Front. Plant Sci. 7:264. doi: 10.3389/fpls.2016.00264
Wawrik, B., Paul, J. H., Bronk, D. A., John, D., and Gray, M. (2004). High rates of ammonium recycling drive phytoplankton productivity in the offshore Mississippi River plume. Aquat. Microb. Ecol. 35, 175–184. doi: 10.3354/ame035175
Wawrik, B., Paul, J. H., and Tabita, F. R. (2002). Real-time PCR quantification of rbcL (ribulose-1,5-bisphosphate carboxylase/oxygenase) mRNA in diatoms and pelagophytes. Appl. Environ. Microbiol. 68, 3771–3779. doi: 10.1128/AEM.68.8.3771-3779.2002
Welschmeyer, N. A. (1994). Fluorometric analysis of chlorophyll-a in the presence of chlorophyll-b and pheopigments. Limnol. Oceanogr. 39, 1985–1992. doi: 10.4319/lo.1994.39.8.1985
Woods, L. A., and Ewing, A. G. (2003). Analysis of single mammalian cells with capillary electrophoresis. Anal. Bioanal. Chem. 376, 281–283. doi: 10.1007/s00216-003-1860-y
Woods, L. A., Roddy, T. P., and Ewing, A. G. (2004). Capillary electrophoresis of single mammalian cells. Electrophoresis 25, 1181–1187. doi: 10.1002/elps.200405842
Worley, B., and Powers, R. (2013). Multivariate analysis in metabolomics. Curr. Metabolomics 1, 92–107.
Xia, J., Sinelnikov, I. V., Han, B., and Wishart, D. S. (2015). MetaboAnalyst 3.0—making metabolomics more meaningful. Nucleic Acids Res. 43, W251–W257. doi: 10.1093/nar/gkv380
Xia, J., and Wishart, D. S. (2016). Using MetaboAnalyst 3.0 for comprehensive metabolomics data analysis. Curr. Protoc. Bioinformatics 55, 14.10.1–14.10.91.
Keywords: Scrippsiella trochoidea, phytoplankton, nutrient limitation, nitrogen, marine algae, metabolomics, single-cell analysis
Citation: Sun M, Yang Z and Wawrik B (2018) Metabolomic Fingerprints of Individual Algal Cells Using the Single-Probe Mass Spectrometry Technique. Front. Plant Sci. 9:571. doi: 10.3389/fpls.2018.00571
Received: 28 December 2017; Accepted: 11 April 2018;
Published: 30 April 2018.
Edited by:
Miyako Kusano, University of Tsukuba, JapanReviewed by:
Robert Winkler, Centro de Investigación y de Estudios Avanzados del Instituto Politécnico Nacional (CINVESTAV-IPN), MexicoManoj Kumar, University of Technology Sydney, Australia
Copyright © 2018 Sun, Yang and Wawrik. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Boris Wawrik, bwawrik@ou.edu