Introduction

Diminishing global oil reserves and increased awareness of the environmental impact of burning fossil fuels have prompted renewed interest in the development of alternative fuel sources in recent years. The viability of microalgae as a carbon-neutral biofuel source has been well documented [14]. In contrast to land-based biomass, microalgae do not compete with food crops for arable land [5, 6] and can grow year-long in photobioreactors with an annual biomass productivity exceeding that of terrestrial plants [7, 8]. Additional benefits of using algae as a biofuel feedstock include possible coupling with waste-water treatment [9] and carbon dioxide sequestration [10, 11]. Despite this potential, algal biofuel production has not yet reached that of terrestrial biomass due to high production costs and energy requirements [12]. For these reasons, research has been devoted to optimizing the algae growth conditions, harvesting conditions, and lipid extraction processes [1]. As new technology is developed and oil prices continue to rise, experts believe that microalgal biofuels will become economically feasible within the next 10 to 15 years [12].

Depending on the strain of microalgae, these single-celled marine and freshwater organisms are capable of producing lipids at levels up to 20–50% of their dry cell weight under specific growth conditions [2]. These neutral lipids, mainly in the form of triacylglycerols (TAGs), can be converted to fatty acid methyl esters (FAMEs) and used as biofuel, much like lipids from terrestrial plants. In addition to TAGs, the lipid fraction of microalgae may also include chlorophylls, carotenoids, glycolipids, phospholipids, and sterols. Under optimal growth conditions, algae synthesize fatty acids for the production of membrane lipids, which include glycolipids and phospholipids. However, under stressed conditions, algae alter their biosynthetic pathways and begin to produce neutral lipids for energy storage, mainly in the form of TAGs [2, 5, 13]. Triacylglycerols are preferred over phospholipids or glycolipids for biofuel production due to their higher percentage of fatty acids and their lack of phosphorous and sulfur. Phosphate hinders the transesterification process, and it has been reported that TAGs have >99% biodiesel yield compared to <70% yield for phospholipids [14]. In addition, sulfur levels in fuels are regulated to minimize vehicle emissions. The fatty acid profile also plays an important role in determining fuel properties. Briefly, saturated fatty acids provide good oxidative stability and ignition properties, while polyunsaturated fatty acids have good cold-flow characteristics [15]. Researchers suggest that a reasonable compromise can be reached by selecting a fuel that is derived from predominantly monounsaturated fatty acids [15]. Therefore, the form of the intact lipids, in addition to the fatty acid profile, is important to establish an alga’s potential as a viable biofuel feedstock.

Conventionally, algal lipid extracts have been analyzed by gas chromatography (GC), following conversion of intact lipids to FAMEs. Detection is often carried out by flame ionization detection, where identities of FAMEs are established by comparison of retention times to known standards [1618], or by mass spectrometry, which offers fragmentation data for fatty acid identification [16, 18, 19]. Typically, GC analysis offers excellent separation efficiency and high-sensitivity, and provides a quantitative fatty acid profile of the lipid extract. However, GC analysis of FAMEs offers no information on the nature of the intact lipids. As a result, GC analysis could be misleading if relying solely on the fatty acid profile for strain selection in biofuels applications. For instance, a favorable fatty acid profile could potentially be achieved from a lipid extract composed primarily of phospholipids, which would not be ideal for biofuel development as described above. Therefore, while GC analysis of FAMEs is perfectly suited for quality control monitoring of finished oils, it is arguably not the ideal technique for many biofuel research and development activities such as strain selection.

Liquid chromatography–mass spectrometry (LC-MS) has been successfully employed for the analysis of various lipid classes such as TAGs [2024] and phospholipids [2431]. For profiling complex lipid mixtures, high-resolution mass spectrometry (HRMS) is typically required and the lower cost of most modern systems had rendered HRMS suitable for high-throughput lipid applications [3234]. In this study, lipid extracts from microalgae were directly analyzed by UHPLC-ESI-MS using a stand-alone benchtop Orbitrap mass spectrometer. Phospholipids, glycolipids, and TAGs were analyzed in the same chromatographic run, using a combination of accurate mass and diagnostic fragment ions for identification. Using this approach, comprehensive TAG profiles were obtained for six different algal strains to assess their potential for biofuel applications.

Experimental

Chemicals and reagents

Liquid chromatography solvents acetonitrile, isopropanol, and methanol were HPLC-grade or glass-distilled and purchased from Caledon (Georgetown, ON, Canada). The following triacylglycerol standards were purchased from Nu-Chek Prep (Elysian, MN, USA): trimyristin (3 × 14:0), triolein (3 × 18:1 cis-9), trilinolenin (3 × 18:3 cis-9,12,15), trieicosenoin (3 × 20:1 cis-11), and tridocosahexaenoin (3 × 22:6 cis-4,7,10,13,16,19). The following polar lipid standards were acquired from Matreya (Pleasant Gap, PA, USA): digalactosyldiacylglycerol (18:0/16:0), monogalactosyldiacylglycerol (18:0/16:0), lyso-phosphocholine (18:1), phosphoglycerol (18:1/18:1), and phosphocholine (18:1/18:1).

Algal growth conditions

Photosynthetic microalgae cultivation was conducted in 250 or 1,000 L internally illuminated Brite-Box photobioreactors [35] in batch mode at 20 °C. Seawater or lake water was sterilized by either pasteurization or tangential flow filtration. The cultures were supplemented with nutrients according to the f/2 formulation [36]. The irradiance was increased gradually from 20% at the time of inoculation to maximum intensity in mid-log phase to minimize the potential for photoinhibition. The maximum light intensity (250 μmol photons m−2 s−1) was measured in algae-free seawater with a flat quantum sensor. The following species were investigated in this study: Botryococcus braunii, Nannochloropsis gaditana, Neochloris oleoabundans, Phaeodactylum tricornutum, Porphyridium aerugineum, and Scenedesmus obliquus. The cultivated biomass was harvested as a wet paste by process centrifugation at the stationary phase of growth and lyophilized to a final moisture content of 3% (w/w). Lyophilized biomass was considered sufficiently dry for solvent assisted lipid extraction.

Lipid extraction and sample preparation

Freeze-dried biomass (homogenized with a mortar and pestle) of approximately 0.5 g was extracted using the Soxtec 2050 automated solvent extraction system (FOSS North America, Eden Prairie MN) with chloroform/methanol (2:1) at 150 °C using the following program: boiling 25 min, rinsing 40 min, solvent recovery 15 min, and pre-drying 2 min. Extracted lipids were dried at 105 °C for 1 h under nitrogen, placed in a desiccator for 1 h, and weighed to determine gravimetric lipid yield. Lipids were resolubilized in chloroform/methanol (2:1) and stored under nitrogen at −20 °C until further analysis. The extracts were then filtered using 0.22-μm Ultrafree-MC centrifugal filter devices (Amicon Bioseparations), then diluted 50-fold in 1:1 methanol/isopropanol prior to injection.

LC-MS instrumentation and conditions

LC-MS data was acquired on an Accela UHPLC system coupled to an Exactive mass spectrometer (Thermo Fisher Scientific, San Jose, CA). The Accela system consisted of a quaternary pump and autosampler, using a Thermo Hypersil Gold C8 column (100 × 2.1 mm, 1.9 μm) at 40 °C. The mobile phase consisted of A, 10 mM ammonium acetate in water (pH adjusted to 5.0 with acetic acid); B, acetonitrile; and C, isopropanol. A gradient with a flow rate of 500 μL/min was used: 30% A/70% B to 100% B from 0 to 5 min, ramped to 95% B/5% C from 5 to 6 min, 20% B/80% C from 6 to 12 min, then held at 20% B/80% C from 12 to 14.5 min. The column was then re-equilibrated back to 30% A/70% B for 5 min. Injection volume was 3.0 μL and each algal extract was analyzed in triplicate.

The Exactive mass spectrometer is a second-generation, stand-alone Orbitrap system, capable of resolution up to 100,000 (FWHM) but without the advanced scanning features offered by the linear ion trap front-end of the conventional Orbitrap mass spectrometer. The system was operated with an optional high-temperature ESI probe (HESI-II) and an “all-ion fragmentation” collision cell known as the HCD cell. In positive ion mode, the source parameters were set as follows: sheath gas, 60; auxiliary gas flow, 15; sweep gas, 0; spray voltage, 3.30 kV; capillary temperature, 380 °C; capillary voltage, 32.5 V; tube lens voltage, 85 V; skimmer voltage, 18 V; heater temperature, 350 °C. Alternating full spectrum scans and HCD scans at 50 eV were collected at 2 Hz using the “High” resolution setting (50,000 FWHM) over a mass range of 100–2,000 m/z.

Results and discussion

TAG identification strategy

Electrospray ionization (ESI) is a soft ionization technique that generates intact TAG molecular ions, typically observed as ammonium or sodium adducts depending on the chromatographic mobile phase employed. Under low-energy collision-induced dissociation (CID), TAGs fragment in a systematic fashion to produce diacylglycerol (DAG) product ions that can be used to identify TAGs through automated database searches. In our case, data were entered into a locally hosted version of the LipidMAPS database [37] that has been modified with a variety of custom fatty acids.

In order to deal with the complexity of natural lipid extracts, some form of data-dependent LC-MS/MS acquisition is typically employed, whereby the most abundant peaks in an MS survey scan are selected for MS/MS scans performed in series. However, these approaches can potentially miss low-abundance TAGs and long cycle times yield very few data points across chromatographic peaks. Therefore, our strategy involved collecting data in a completely untargeted fashion by toggling between MS scans and “all-ion fragmentation” scans offered by the Exactive mass spectrometer. The challenge of this approach is in the assignment of the correct DAG product ions to potentially co-eluting TAGs when multiple combinations are possible. However, in these rare occurrences, this can be mitigated by aligning extracted ion chromatograms of the DAG fragments with those of the intact TAGs, exploiting small shifts in retention times of nearly co-eluting TAGs.

Typical data used for TAG identification are shown in Fig. 1. The MS spectrum (Fig. 1a) shows the accurate mass of an intact lipid from the algae B. braunii at retention time 11.33 min, which can be searched against a database to yield a match within 1 ppm of an ammoniated TAG with a total of 64 carbons and four double bonds in its fatty acyl chains (denoted as TAG(64:4)). Subjecting the intact TAGs to the HCD collision cell (Fig. 1b) yields diacylglycerol fragment ions, which indicate the neutral loss of a fatty acid from the glycerol backbone. This particular TAG shows a neutral loss of fatty acids 18:1 and 28:2, indicating that the TAG(64:4) identity is 28:2/18:1/18:1. Determination of the position of fatty acyl chains within each TAG was beyond the scope of this study. It should be noted that while ammonium adducts were the predominant form of intact TAG ions in the MS scan (Fig. 1a), the relative abundance of the sodium adduct is enhanced in the all-ion fragmentation scan (Fig. 1b). However, this phenomenon is simply due to the increased stability of the sodium adducts, whereas ammonium adducts require less energy for fragmentation. Limits of detection (LODs) measured for five TAG standards ranged from 11 to 19 nmol/L, corresponding to 30–56 pg injected on-column.

Fig. 1
figure 1

Typical LC-MS spectrum (a) and “all-ion fragmentation” spectrum acquired with the HCD cell (b) for TAG 28:2/18:1/18:1 at retention time 11.33 min for B. braunii. Fragment ions in (b) corresponding to neutral losses (NL) of fatty acids 28:2 and 18:1 are indicated

Liquid chromatography optimization

Liquid chromatography conditions were optimized to separate the various lipid classes found in algal extracts as much as possible, while maintaining a reasonably short run time to accommodate moderate sample throughput. Using a tertiary mobile phase system, more polar lipids were eluted first with a water/acetonitrile gradient, followed by elution of TAGs with acetonitrile/isopropanol. As shown in Fig. 2, various classes of lipids from a B. braunii extract were separated in roughly 12 min, as outlined in the three-dimensional plot of intensity as a function of retention time and mass-to-charge ratio. Triacylglycerols elute as sharp peaks between 8 and 12 min with masses ranging from roughly m/z 750 to 1,200, while more polar lipids generally eluted in the first 8 min of the chromatogram. Although not the main objective of this study, several classes of these moderately polar lipids were identified in B. braunii and the other algal strains studied, and their position in the chromatogram is illustrated in Fig. 2. Identification was carried out by matching observed masses within 2 ppm to calculated masses in the LipidMAPS database, and lipid classes were confirmed by characteristic fragment ions. For instance, phosphocholines yield a diagnostic fragment ion at m/z 184.07, and are therefore easily revealed through extracted ion chromatograms of this mass, while phosphoglycerols exhibit a neutral loss of the sodiated headgroup (193.99 Da). A variety of other polar lipids were also identified in B. braunii, including monogalactosyldiacylglycerols, digalactosyldiacylglycerols, sulfoquinovosyldiacylglycerols, 1,2-diacylglyceryl-3-O-4′-(N,N,N-trimethyl)-homoserines, and other lipophilic compounds such as chlorophylls, as indicated in Fig. 2. With the exception of phosphocholines, which eluted as broad peaks in the 9–11-min range, all of these lipid classes were well-resolved from the TAG region of the chromatogram.

Fig. 2
figure 2

3D ion map of LC-MS chromatogram of B. braunii lipid extract. Abbreviations: PC phosphocholines, DGTS 1,2-diacylglyceryl-3-O-4′-(N,N,N-trimethyl)-homoserines, PG phosphoglycerols, MGDG monogalactosyldiacylglycerols, DGDG digalactosyldiactylglycerols, SQDG sulfoquinovosyldiacylglycerols. LC-MS conditions as described in the section “LC-MS instrumentation and conditions

Triacylglycerol profiling of algal strains

The LC-MS strategy described above was employed to profile lipids in six different algal strains that were cultured and extracted under identical conditions. As shown by the LC-MS base peak chromatograms in Fig. 3, the lipid compositions of the six strains varied significantly over the range of algal species studied. For instance, it is evident from the chromatograms that only B. braunii and S. obliquus produced significant amounts of TAGs under the growth conditions employed, as demonstrated by the relatively intense peaks in the 8–12 min range of the chromatograms. In contrast, N. oleoabundans yielded a high abundance of polar lipids in the early portion of the chromatogram, as well as broad phosphocholine peaks in the later portion of the chromatogram, but relatively small amounts of TAGs. It should be noted that due to varying degrees of ionization efficiencies for the different lipid classes, a quantitative comparison cannot be made between the different lipid classes within the same chromatogram. However, a semi-quantitative comparison can be made for the same lipid class across different algal strains.

Fig. 3
figure 3

Base peak chromatograms of lipid extracts of 6 strains of microalgae, extracted with 2:1 chloroform/methanol, in positive-mode ESI-MS. LC-MS conditions as in section “LC-MS instrumentation and conditions

As triacylglycerols offer the most efficient source of hydrocarbons for various biofuels applications, comprehensive lipid profiles were obtained for the six strains studied to evaluate their potential for biofuels feedstock. Using the TAG identification strategy described above, over 100 unique TAGs were identified in the six algal strains studied. Relative abundances of each TAG were determined by comparing peak areas of extracted ion chromatograms averaged over triplicate injections. As response factors measured with available standards for a subset of TAGs varied by up to 50%, accurate quantification of all TAGs was not possible. Therefore, a heat map was constructed using an appropriate scale spanning nearly three orders of magnitude to highlight significant differences in TAG abundances. As shown in Fig. 4, the TAG profiles varied significantly between species, with B. braunii and S. obliquus exhibiting the broadest range of TAGs. In particular, over 70 different TAGs were identified in B. braunii, with the TAGs 18:1/18:1/18:1, 18:1/18:1/18:3, 28:1/18:1/18:1 and 28:2/18:1/18:1 being the most abundant. Other TAGs containing long-chain fatty acids such as 30:1, 30:2, and 36:2, were also detected in B. braunii. Although relatively rare, these fatty acids are known to be present at significant abundance in B. braunii [38]. Further, both B. braunii and S. obliquus also have a high abundance of TAGs containing oleic acid (18:1), which is favorable for biodiesel [15]. Also listed in Fig. 4 are retention times and accurate mass measurement data for the TAGs identified. Measured masses with standard deviations were based on averages of three replicates. Mass errors were reported as the average of the absolute values of the individual mass errors [39], to eliminate positive and negative errors canceling each other.

Fig. 4
figure 4

Heat map of TAGs identified in six algae strains. BB Botryococcus braunii, NG Nannochloropsis gaditana, PT Phaeodactylum tricornutum, NO Neochloris oleoabundans, PA Porphyridium aerugineum, SO Scenedesmus obliquus

Conclusion

An ultra-high pressure liquid chromatography–mass spectrometry method was developed for profiling lipid extracts to aid in identifying algae strains suitable for biofuels applications. The 15-min chromatographic separation allowed both polar lipids and triacylglycerols to be resolved in the same run. The method requires little sample preparation and is suitable for moderate throughput applications. Profiling intact lipids offers valuable complementary data to conventional fatty acid profiling by gas chromatography for discovery-phase biofuels research. For instance, algal culturing conditions can be optimized to specifically promote TAG production, as opposed to simply overall fatty acid content. It should be noted that the algal growth conditions employed in this study were not optimized for maximum lipid yield for each strain. The generic conditions used for all strains were chosen to provide a fair comparison for lipid profiling and would be suitable for screening new or uncharacterized algal strains for which optimal growth conditions are unknown. Therefore, higher lipid content could be expected for some of the strains profiled in this study with proper optimization of growth conditions.