Suitability of Phytosterols Alongside Fatty Acids as Chemotaxonomic Biomarkers for Phytoplankton

The composition and abundance of phytoplankton is an important factor defining ecological status of marine and freshwater ecosystems. Chemotaxonomic markers (e.g., pigments and fatty acids) are needed for monitoring changes in a phytoplankton community and to know the nutritional quality of seston for herbivorous zooplankton. Here we investigated the suitability of sterols along with fatty acids as chemotaxonomic markers using multivariate statistics, by analyzing the sterol and fatty acid composition of 10 different phytoplankton classes including altogether 37 strains isolated from freshwater lakes. We were able to detect a total of 47 fatty acids and 29 sterols in our phytoplankton samples, which both differed statistically significantly between phytoplankton classes. Due to the high variation of fatty acid composition among Cyanophyceae, taxonomical differentiation increased when Cyanophyceae were excluded from statistical analysis. Sterol composition was more heterogeneous within class than fatty acids and did not improve separation of phytoplankton classes when used alongside fatty acids. However, we conclude that sterols can provide additional information on the abundance of specific genera within a class which can be generated by using fatty acids. For example, whereas high C16 ω-3 PUFA (polyunsaturated fatty acid) indicates the presence of Chlorophyceae, a simultaneous high amount of ergosterol could specify the presence of Chlamydomonas spp. (Chlorophyceae). Additionally, we found specific 4α-methyl sterols for distinct Dinophyceae genera, suggesting that 4α-methyl sterols can potentially separate freshwater dinoflagellates from each other.


INTRODUCTION
The composition and abundance of phytoplankton is an important factor defining secondary and primary production, since phytoplankton synthesize many essential biomolecules (e.g., fatty acids, sterols, amino acids), which consumers cannot synthesize de novo (Harrison, 1990;Sargent et al., 1999;Lafont, 2000;Brett et al., 2009;Martin-Creuzburg and Von Elert, 2009). However, in aquatic ecology phytoplankton are often assumed to be a homogenous group and a uniform food source for higher trophic levels, primarily in studies using stable isotopes (Karlsson et al., 2003;Taipale et al., 2008;Solomon et al., 2011). The dietary quality of distinct phytoplankton classes for zooplankton is nowadays known to vary primarily due to the different availability of ω-3 and ω-6 fatty acids (Jonasdottir, 1994;Brett et al., 2009), but also because of lack of sterols (Von Elert et al., 2003). Compared to the laborious work of microscopic identification of phytoplankton, chemotaxonomic markers (e.g., pigments and fatty acids) provide a faster method for monitoring changes in a phytoplankton community (Mackey et al., 1996;Strandberg et al., 2015), but also provide information on the nutritional quality of seston for the whole food web. The use of sterols as chemotaxonomic biomarkers alongside fatty acids could potentially enhance separation of phytoplankton classes, but also provide information of the potential limitation of sterols in zooplankton diet. However, the similarity of sterol profiles among distinct phytoplankton classes is not yet analyzed using modern multivariate statistics, which have been previously used for fatty acids in freshwater phytoplankton, marine macrophytes, and macroalgae (Galloway et al., 2012;Kelly and Scheibling, 2012;Taipale et al., 2013).
Due to the great nutritional value and conservative transfer of fatty acids on higher trophic levels, fatty acids have been used as trophic markers (FATM, Fatty Acid Trophic Markers; Dalsgaard et al., 2003) for providing insight into consumer diets in pelagic and benthic food webs (Stott et al., 1997;Kelly and Scheibling, 2012). Fatty acid profiles of freshwater microalgae and marine macroalgae are mainly phylogenydependent, and environmental conditions influence the fatty acid composition only slightly (Galloway et al., 2012;Taipale et al., 2013;Galloway and Winder, 2015). Thus, a Bayesian mixing model based on fatty acids has been used successfully for describing phytoplankton community structure (Strandberg et al., 2015) as well as for taxonomic primary production (Dijkman et al., 2009). Non-metric multidimensional scaling analysis (NMS) of freshwater phytoplankton can separate green algae (Chlorophyceae, Trebouxiophyceae, Conjugatophyceae), diatoms (Diatomophyceae), and euglenoids (Euglenophcyeae) into their own groups, but clusters Cryptophyceae, Chrysophyceae, and Raphidophyceae in one group (Taipale et al., 2013). Additionally, fatty acid profiles of Cyanophyceae are highly variable (Los and Mironov, 2015), and thus chemotaxonomic markers for better phytoplankton phylogenetic separation are needed.
Sterols could give additional information on phytoplankton community structure, but also together with the essential ω-3 and ω-6 fatty acids on the biochemical quality of phytoplankton for aquatic food webs since both are essential for herbivorous consumers (Brett et al., 2009;Martin-Creuzburg and Von Elert, 2009). Phytosterols have been intensively studied since they were found in the 1930s (Carter et al., 1939;Patterson, 1991). But, perhaps due to the long research history of sterols, the same sterol can have several names in the literature, thus making sterol nomenclature confusing and complex. Phytosterols are sterols that are structurally similar to cholesterol but are synthesized in plants and phytoplankton. Usually, they are characterized by alkyl groups (methyl or ethyl group) at C-24, by side-chain double bonds at C-22 or C-24 or/and by nuclear double bonds at positions other than C-5 (Morreau et al., 2002). Among phytoplankton, the Dinophyceae are an exception because they can have a methyl group in the C-4 or C-23 position (Withers et al., 1978;Piretti et al., 1997). Some Chlorophyceae contain unusual 7 -sterols such as chondrillasterol (Patterson, 1991;Thompson, 1996), whereas some Chlorophyceae in the order Chlamydomonadales have ergosterol (Gealt et al., 1981;Brumfield et al., 2010), which is also found in Euglenophyceae (Patterson, 1991). Full chemical names for sterols are included at Table 1. There is a great diversity of sterols in Diatomophyceae, but none of them are found exclusively in diatoms (Rampen et al., 2010). For example brassicasterol is also a major sterol in Cryptophyceae (Rampen et al., 2010) and stigmasterol, and β-sitosterol are typical for Chrysophyceae and Raphidophyceae in addition to Diatomophyceae (Cranwell et al., 1988;Rampen et al., 2010;Leblond et al., 2013).
Here we have studied the fatty acid and sterol profiles of 35 major freshwater phytoplankton genera (Tables 1, 2). Our 37 phytoplankton strains belong to 10 phytoplankton classes: Chlorophyceae, Conjugatophyceae, Cryptophyceae, Cyanophyceae, Diatomophyceae, Dinophyceae, Euglenophyceae, Raphidophyceae, Synurophyceae, and Trebouxiophyceae. We used multivariate statistics to describe the similarities and differences in fatty acid and sterol composition among freshwater phytoplankton classes and explored if the addition of sterols improves the differentiation between phytoplankton classes.

Phytoplankton Culturing
The algal classification follows the taxonomy and trivial names of AlgaeBase (Guiry and Guiry, 2014). Class, order, and species of the cultured strains are presented in Table 3. The cultured phytoplankton genera were chosen based on the most abundant taxa in the phytoplankton database of the Finnish Environment Institute (SYKE). The phytoplankton data represent on average 81% of the total phytoplankton biomass in Finnish lakes during open water season (June-September, in 1975, and they also represent the most abundant phytoplankton genera in lakes across the boreal zone, including Northern Great Britain, Norway, Sweden, and Finland (Järvinen et al., 2013;Maileht et al., 2013). Phytoplankton composition varies in boreal lakes by total phosphorus (µg l −1 ), water color (mg l −1 Pt), and Secchi-disc transparency (Lepistö and Rosenström, 1998).
Most of the phytoplankton strains were grown at 18 • C under a 14:10 h light:dark cycle with light intensity of 30-70 µmol m −2 s −1 . Each strain was cultured in a medium specific to that strain (Table 3). Additionally, Tabellaria (strain 19 and 20) was cultured in two different media (Chu10 and Z8) and Euglena gracilis (strain 27 and 28) with (EG) and without (AF6) organic substrates ( Table 3). The strains of Cryptomonas Relative retention time (RRT) was calculated in comparison to the retention time (RT) of cholestane (RRT = 1). The systematic name follows the IUPAC nomenclature and the trivial name is the most commonly used trivial name of the sterol. Position of the double bond refers to the carbon number which has a covalent double bond and C-24 alkylation shows what kind of alkyl group is bound to carbon number 24 (LIPID MAPS). Sterols were run as trimethylsilyl ethers (TMS) and were identified using specific ions, including the mass ion (TMS M + ), and the most abundant mass ion (base peak; Jones et al., 1994).
Frontiers in Plant Science | www.frontiersin.org  Fatty acids are presented in the order they eluted (RT, retention time in minutes) from the DB-23 column. Reference ion was used for fatty acid identification and major peak for quantification. Names for fatty acids are presented as C:D:ω; C, the number of carbon in the chain; D, the number of double bonds; ω, position of first double bond. Additionally, c and t cites refer to cis or trans transfiguration. Abbreviation i-(iso) and a-(anteiso) cites to fatty acids with methyl branched position at second or third carbon from the end. Cycites to the cyclopropane fatty acid.
(strain 17), Gonyostomum semen (strain 18), E. gracilis (strain 27 and 28), and Selenastrum sp. (strain 35) were cultured at 20 • C with light instensity of 40 µmol m −2 s −1 and using a 16:8 h light:dark cycle. We used plastic or glass flasks, volume > 200 ml. Depending on the cell density, 1-6 ml of the phytoplankton stock was inoculated per 200 ml of fresh culture media every 2 weeks. The samples for phytoplankton analyses were harvested in the late phase of exponential growth, i.e., 2-3 weeks after the inoculation.

Fatty Acid Analyses
Lipids from freeze-dried, homogenized phytoplankton (1-3 mg) were extracted using chloroform: methanol: NaCl mixture with volumes of 2, 1, 0.75 mL of chloroform, methanol, and 2% NaCl as reported in Parrish (1999). Phytoplankton samples were sonicated for 10 min and vortexed (2-3X) and the lower organic phases were removed after which 2 mL of chloroform was added. After centrifuging the lower organic phase was removed and pooled, and the organic solvent was evaporated to dryness. For the formation of fatty acid methyl esters (FAME), 1 mL of toluene, and a 2 mL of sulphuric acid-methanol solution were used with dried lipids (all of each sample), and incubated at 80 • C in a water bath for 2 h. Samples were neutralized by adding 2 mL of 2% KHCO 3 and FAMEs diluted in 5 mL of hexane. The upper layer was transferred into a clean centrifuge tube after which this step was repeated, and finally hexane was evaporated to dryness under nitrogen flow and the sample was transferred to autosampler vials with inserts (Agilent, Santa Clara, United States) with 200 µL hexane. FAMEs were analyzed with a gas chromatograph (Shimadzu Ultra, Kyoto, Japan) equipped with mass detector (GC-MS) and using helium as a carrier gas. An Agilent R (Santa Clara, California, U.S.A.) DB-23 column (30 m × 0.25 mm × 0.15 µm) was used with the following temperature program: 60 • C was maintained for 1.5 min, then the temperature was increased at 10 • C min −1 to 100 • C, followed by 2 • C min −1 at 140 • C, and 1 • C min −1 at 180 • C and finally heated at 2 • C min −1 to 210 • C and held for 6 min. Fatty acids were identified by the retention times (RT) and using specific ions ( Table 2). Fatty acid concentrations were calculated using calibration curves based on known standard solutions of a FAME standard mixture (GLC standard mixture 566c, Nu-Chek Prep, Elysian, Minnesota, U.S.A.). The Pearson correlation coefficient was >0.99 for each individual fatty acid calibration curve (details of the method in Taipale et al., 2013).

Sterol Analyses
Replicates of dried total lipid samples (subsamples of the ones used for fatty acids) of phytoplankton were used for sterol analysis. To separate sterols from fatty acids, lipid samples of cultured strains were saponified by adding 2 mL of 1.2 N methanolic KOH and samples was kept at 70 • C for 2 h (in water bath). After this saponification, 1 mL of 2% sodium chloride and 3 mL of hexane were added to the samples, and they were vortexed and centrifuged. The upper organic layer was transferred to a new acetone-washed tube (precombusted at 440 • C for 4 h) and 3 mL of hexane were added, and then this step was repeated to maximize the yield of the non-polar phase. After evaporating hexane to dryness, samples were dried with acetone for removing any residues of water or moisture and to maximize the silylation reaction in the next stage.
For silylation, 0.2 mL of N,O-bis[trimethylsilyltrifluoroacetamide] (BSTFA) with 1% (w) trimethylchlorosilane (TMCS) (Fluka Sigma-Aldrich, St. Louis, Missouri, U.S.A.) and 0.2 mL of pyridine (Sigma-Aldrich, St. Louis, Missouri, U.S.A.) were added to the samples and incubated for 2 h at 70 • C, after which the samples were evaporated and transferred to glass vials with inserts. The samples were dissolved in 50-200 µl of hexane. Trimethylsilyl (TMS) derivatives of sterols were analyzed with a gas chromatograph (Shimadzu) equipped with a mass detector. A Phenomenex R (Torrance, California, U.S.A.) ZB-5 Guardian column (30 m × 0.25 mm × 0.25 µm) was used with the following temperature program: 150 • C was maintained for 1 min, then the temperature was increased at 15 • C min −1 to 280 • C, finally heated by 2 • C min −1 to 320 • C, and held for 10.3 min. The total length of the temperature program was 40 min. Retention times and relative retention times (RRT) for detected sterols are presented in Table 1. Helium was used as a carrier gas with average velocity of 34 cm sec −1 . Sterols were quantified using authentic standard solutions of plant sterol mixture from Larodan (Solna, Sweden; including 53% β-sitosterol, 7% stigmasterol, 26% of campesterol, 13% of brassicasterol), and cholesterol and fucosterol from Sigma-Aldrich. Sterol concentrations were calculated using four point calibration curves based on known standard concentrations (between 1.3 and 159 µg/µL). The Pearson correlation coefficient was >0.99 for each individual sterol calibration curve. Recovery percentage of sterol samples was calculated using 5-α-cholestane (Sigma-Aldrich) as an internal standard.
The identification of TMS ethers of sterols was based on the standard mixes described above, retention times together with the NIST Mass Spectral Database 11 (http://chemdata. nist.gov) and literature (Rahier and Benveniste, 1989;Jones et al., 1994;Goad and Akhisa, 1997). The molecular ion peak, molecular ion −15 (M+-15), base peak and other used ions for identification are presented in Table 1.

Statistical Analyses
We used non-metric multidimensional scaling (NMS) to visualize the variation in fatty acid and sterol composition and in a combined dataset of sterols and fatty acids (contribution of each sterol and fatty acids was calculated from the sum concentration of all sterols and fatty acids) among phytoplankton classes. Permutational multivariate analysis of variance (PERMANOVA, Anderson et al., 2008) was used to test whether these differences in the sterol and/or fatty acid composition were statistically significant. PERMANOVA was run with unrestricted permutation of raw data and type III sums of squares. Similarity percentages (SIMPER) were used to identify the characteristic sterols and fatty acids of each phytoplankton class. We used PERMDISP (Anderson, 2006) to investigate the within-class variation in fatty acid and/or sterol composition. All the multivariate analyses were operated on Euclidean distances of untransformed data with the program PRIMER-E (v.6; Ivybridge, United Kingdom) and the PERMANOVA+ add-on.
The within-class multivariate dispersion varied between different phytoplankton classes in sterol, fatty acid, and combined signature datasets (Figure 3). Generally, dispersion varied among classes and the variation within phytoplankton class was higher in sterol signatures [PERMDISP; F (9, 29) =5.862] than in fatty acids [PERMDISP; F (9, 29) =9.6037] or in combined signatures [PERMDISP; F (9, 29) =9.131] (Figure 4). Sterol signatures were most similar in the Synurophyceae and Euglenophyceae, otherwise the mean distance-to-centroid in sterol signatures was two or even three times higher than in fatty acid signatures. The similarity of class was equal when only fatty acid data was used and also in the combined dataset of fatty acids and sterols, but in some cases sterols increased the integrity of the class (e.g., in Diatomophyceae).
Five sterols and fatty acids that contributed most to within-group (phytoplankton class) similarity and their mean proportions were identified with the SIMPER analysis ( Table 5). The "Contributions" are the percent contributions of sterols or fatty acids to the similarities within the taxa in that class. When SIMPER analysis was performed with the combined datasets of sterols and fatty acids, only the classes Synurophyceae, Euglenophyceae, and Conjugatophyceae had sterols among the five most important classifying contributors. Sterols and fatty acids were most responsible for within-group similarities and also played an important role in separating the phytoplankton groups in the NMS (Figure 3).
NMS of the combined data did not differ from the NMS ordination of the fatty acid data, and axes 1, 2, and 3 correlated with the same fatty acids as in the fatty acid data alone ( Table 6). The only difference was that NMS 3 correlated positively (r > 0.44, p < 0.05) also with cholesterol.

DISCUSSION
Previous studies of the fatty acids of freshwater phytoplankton has shown separation at the class level (Taipale et al., 2013), but the sterol composition of different freshwater phytoplankton taxa has not been previously studied systematically using multivariate statistical analysis. Our analysis of 37 dominant freshwater phytoplankton species from 10 classes showed that even though both fatty acids and sterols varied significantly among phytoplankton classes, the integrity of sterol signatures within classes was lower than in fatty acids. The combined dataset FIGURE 1 | Fatty acid profiles of cultured freshwater phytoplankton strains. Strain numbers are given in Table 2. Fatty acids are presented as major groups: saturated fatty acids (SAFA), C 16 and C 18 monounsaturated fatty acids (C 16 MUFA, C 18 MUFA), branched fatty acids, and C 16 , C 18 , C 20 , and C 22 polyunsaturated fatty acids (PUFA).
of sterols and fatty acids showed similar multivariate dispersion to the fatty acid data (Figure 3), but in some cases mean distance-to-centroids among class decreased with the addition of sterol data (Figure 4). Phytoplankton class explained ∼70% of the variation in fatty acid signatures and combined fatty acid and sterol signatures when Cyanophyceae were excluded. The FIGURE 2 | Sterol profiles of the studied phytoplankton strains. Orange-brown represents 7 unsaturation, yellow represents 7, 22 diunsaturation, blue represents 5 unsaturation, green represents 5, 22 diunsaturation, black represents 5, 24(28) diunsaturation, red represents 5, 7, 22 triunsaturation, and violet represents 4α-methyl and 23,24 -trimethyl sterols.
Frontiers in Plant Science | www.frontiersin.org Statistically significant differences between two classes (Cya, Cyanophyceae; Cry, Cryptophyceae; Dia, Diatomophyceae; Chlo, Chlorophyceae; Con, Conjugatophyceae) are marked with stars. The factor Class explained most of variation when only fatty acids were included into analyses.  Table 3. Detailed list of correlations is presented in Table 6.
percentage decreased 20% when Cyanophyceae were included in the fatty acid or the combined data ( Table 4). While fatty acid profiles of phytoplankton are characteristic at class level, sterols seem to give information of the genus level for some phytoplankton classes. The average numbers of detected fatty acids and sterols from individual freshwater phytoplankton strain were 19 and 4, respectively. Rampen et al. (2010) found 44 sterols in 106 marine diatoms, of which 15 were abundant and the rest were only found in less than 3% of their strains. In addition to phytoplankton, in many other organisms only one or a few sterols dominated the profile (Goad and Akhisa, 1997;Martin-Creuzburg and Von Elert, 2009). While the absence or presence of one fatty acid has a minor role in dissimilarities among taxa, a difference of one sterol could result in high dissimilarity in multivariate statistics. Therefore, we also investigated the combined data of fatty acids and sterols. We found that the contribution of sterols to total biomass of fatty acids and sterols was generally low and thus sterols had a minor role in separating phytoplankton classes. In NMS analysis (Figure 3) cholesterol was the only sterol correlating with any of the three NMS dimensions, and it correlated positively together with EPA and DHA with the NMS 3 that was responsible for separating Dinophyceae from other phytoplankton classes.
The multivariate analysis and Bayesian mixing modeling (Strandberg et al., 2015) separate different phytoplankton phylogenetical groups based on similarities and dissimilarities in biomolecule profiles. Therefore, models do not emphasize specific single biomarkers, and thus e.g. dinosterol, was not a significant component for separating Dinophyceae from other phytoplankton groups. It should also be remembered that even if phytoplankton are lacking some of the fatty acids and sterols characteristic to other phytoplankton groups, NMS clusters samples using similarities based on the whole profiles. For example Aphanothece and Synechococcus did not contain any sterols, or ω-3 or ω-6 fatty acids, but NMS clustered them together with Diatomophyceae due to their high 16:1ω7 content in these cyanophyceael strains. Therefore, for better separation of Cyanophyceae from algae more Cyanophyceae and other phytoplankton taxa should be analyzed for sterols and fatty acids. Also, a higher number of studied strains in all classes would give better separation and more detailed understanding on the relationship between sterols and phytoplankton phylogeny.

Suitability of Lipid Biomolecules of Metabolic Pathways for Biomarkers
Specific chemotaxonomic fatty acids or sterols for phytoplankton taxa are often by-products of metabolic pathways. For example, the polyunsaturated fatty acids 16:2ω6 and 16:3ω3 are precursors of 16:4ω3 during the ω3 desaturation pathway from 16:1ω9 (Erwin, 1973), and all of these fatty acids are found in Chlorophyceae, Conjugatophyceae, and Trebouxiophyceae. Furthermore, Diatomophyceae contain 6 desaturase activity that can desaturate 16:1ω9 to16:2ω4 and 16:3ω4, but 16ω3 PUFA are not found in Diatomophyceae (Perkins and Witting, 1975). Both ω3 and ω6 C 18 PUFA are synthesized from oleic acid via desaturation and elongation. SDA (stearidonic acid, 18:4ω3) is synthesized from ALA (α-linolenic acid, 18:3ω3) by desaturation, but Dinophyceae prefer to retroconvert EPA to 18:5ω3 than desaturate SDA to 18:5ω3 (Perkins and Witting, 1975). However, whereas marine and brackish water dinoflagellates have been reported to have 18:5ω3 and DHA as their major FA (Leblond et al., 2006;Mooney et al., 2007), we found only trace amounts of 18:5ω3 in Dinophyceae, and more SDA and ALA. We also found higher amounts of EPA especially in Peridinium cintum than in marine or brackish water Dinophyceae strains, and Peridinium did not contain any 18:5ω3, indicating that the higher EPA content can be the result of lower retroconversion activity. Among phytoplankton, the Cyanophyceae is the most diverse phytoplankton class according their FA profiles (Figures 3, 4), and especially their ability to desaturate oleic acid (18:1ω9) to C 18 PUFA varies greatly between distinct genera, which have been used in previous studies in classifying Cyanophyceae into four groups (Kenyon, 1972, Murata et al., 1992, Los and Mironov, 2015. The cultured Cyanophyceae, Aphanothece, and Synechococcus, lacked C 18 PUFA, but contained high amounts of C 16 MUFA (monounsaturated fatty acid), and thus they belong FIGURE 4 | PERMIDISP was used to evaluate dispersion within each phytoplankton class. Highest mean distance-to-centroid was measured when sterols (STE) were used alone, but fatty acids (FA) and combined data of sterols and fatty acids (STE + FA) have similar dispersion.
to the first group. Additionally, even though Phormidium was able to synthesize C 18 PUFA, we found only limited amount of 18 PUFA in this genus. The second group of Cyanophyceae has the first double bond of triunsaturated C 18 PUFA at the position ω3 and group 3 at the position ω6. The only genus of cultured Cyanophyceae in our study that belonged to group 2 was Microcystis, whereas the rest of the cultured strains (Snowella, Anabaena, Planktothrix, Limnothrix, Pseudanabaena) belonged to group 3. The fourth group had four acyl-lipid fatty acid desaturases and can synthesize SDA (Los and Mironov, 2015). However, none of our cyanobacterial strains belonged to this group.
Algal sterol biosynthesis is much more complicated than that of fatty acids, having more steps toward the final product and resulting in many subproducts and a high number of different phytosterols (Nes and McKean, 1977). Algal sterols are synthesized via the mevalonate (MVA) or methyl-D-erythritol 4-phosphate (MEP) pathways of isoprenoid biosynthesis, and the synthesis requires approximately 30 different enzymes. Our goal was to analyze the major phytosterols in phytoplankton but not to analyze all sterols existing in trace amounts that are the precursors of some other sterols. Chemotaxonomical differences in synthesized sterols were found among classes even though some classes did not differ statistically from each other. Chlamydomonas has the ability to synthesize egrosterol and corbisterol from cycloartenol that differs from the fungal acetate-mevalonate pathway (Miller et al., 2012). In addition to Chlamydomonas we found high amounts of ergosterol and corbisterol in Euglena. Ergosterol and corbisterol (double bond positions at 5, 7, 22 ) have been previously found in some other marine Chlorophyceae genera (Dunaliella, Polytoma, Haematococcus, Chlorella) (Patterson, 1991). Unusual 7 and 7, 22 -sterols were only synthesized among Chlorophyceae in Eudorina, Monoraphidium, Pediastrum, Acutodesmus, and Selenastrum, which are the most common green algae in addition to Chlamydomonas in boreal lakes. Additionally, 7 and 7, 22 -sterols were the major sterols in Closterium and trace amounts of these sterols were found in Cosmarium, both members of Conjugatophyceae. Furthermore, other studies have found 7 and 7, 22 -sterols in two other marine Chlorophyceae: Ankistrodesmus and Oocystis (Patterson, 1991), that are also found in freshwater systems. Most dinoflagellates contain high amounts of 4α-methyl sterols, which are rare among other phytoplankton taxa (Mansour et al., 1999;Volkman, 2003). Among the two cultured Dinophyceae, Peridinium had more 4α-methyl sterols than Ceratium, which contained <5% of 4α-methyl sterols. Altogether, it could be concluded that by including only the fatty acids or sterols synthesized through pathways common only in one phytoplankton group, like here in Dinophyceae, the discrimination among classes can be enhanced in multivariate analysis.

Sterols and Fatty Acids as Chemotaxonomic Biomarkers
Biomolecules can be useful biomarkers in food web studies for separating phytoplankton, terrestrial matter and bacteria, and to quantify their importance in the diets of zooplankton. Fatty acids have already been used for defining phytoplankton composition (Strandberg et al., 2015), but also for dietary analysis of zooplankton, fish and seals (Sargent et al., 1987;Iverson et al., 1997;Galloway et al., 2014). This is possible since the fatty acid analysis of freshwater phytoplankton and marine macrophytes has shown that the variation in fatty acid profiles is largely explained by taxonomic identity, and thus it seems that environmental conditions such as nutrients and light intensity have relatively lower impact on fatty acid composition (Galloway et al., 2012;Taipale et al., 2013;Bi et al., 2014;Galloway and Winder, 2015). Sterols can be used only for phytoplankton community analysis and not for trophic transfer since cholesterol is a major sterol in zooplankton and fish, is a precursor for steroid hormones (Grieneisen, 1994), and is a required for developmental patterning of animal embryonic structures (Porter et al., 1996).
Our PERMANOVA analysis showed that 50% of the variation in phytoplankton sterol composition was explained by class, and thus the major part of sterol variation was explained by phylogeny. Environmental conditions may affect the concentrations of individual sterols, and for example light, phosphorus, and temperature influence the abundance of 7sterols and 7, 22 -sterols (fungisterol, chondrillasterol and 22dihydrochondrillassterol) in Scenedesmus (Veron et al., 1996;Piepho et al., 2010Piepho et al., , 2012. In these studies Scenedesmus did not synthesize 5 -sterols under any conditions, however, the growth phase influenced dinoflagellates to synthesize different sterols (Amo et al., 2010), thus showing some flexibility in sterol synthesis in freshwater algae.
The traditional biomarker approach is looking for specific biomolecules that are absent or rare in other organisms, but as mentioned above present multivariate models do not require specific biomarkers. Nevertheless, these specific biomarkers can be used to track biomass of specific phytoplankton class or genus. Previous studies have identified dinosterol and related sterols as unique biomarkers for dinoflagellates (Volkman et al., 1993;Piretti et al., 1997). Additionally, our previous study (Taipale et al., 2013) of freshwater phytoplankton showed C 15 to C 17 PUFA to be specific and useful biomarkers for Chlorophyceae and Trebouxiophyceae (16:2ω6, 16:3ω3 and 16:4ω3), Diatomophyceae (16:3ω4), and Euglenophyceae (C 15 and C 17 PUFA). However, simultaneous monitoring of fatty acids and sterols can provide more detailed information of phylogeny but also of nutritional adequacy. The simultaneous presence of 16ω3 PUFAs and ergosterol would indicate Chlamydomonas whereas 5, 22 -sterols would point to Sphaerocystis, and 7sterols to other genera of Chlorophyceae and Conjugatophyceae. Another useful approach could be to monitor 4α-desmethyl sterols simultaneously with DHA for tracking the abundance of Dinophyceae; this could also provide information on nutritional quality since DHA is physiologically the most important FA for copepod zooplankton and fish.
Previous studies have shown that sterol-coding genes (Villanueva et al., 2014) and 18S rRNA sequences together with sterol composition of marine diatoms (Rampen et al., 2010) can be used for forming algal groups. However, presently it's not possible to use same primers for all phytoplankton classes (including Cyanophyceae and microalgae), and thus we were forced to rely on present taxonomy. Phytoplankton taxonomy is based more on morphology than on gene sequencing, and therefore taxonomy relies on subjective opinions and present consensus (Reynolds et al., 2002). Thus, dissimilarities in sterol composition within present taxa do not necessarily mean an inconsistent pattern of sterol synthesis, but may indicate problems in species identification and classification. Of the strains we used, Botryococcus is placed under the Chlorophyceae in some classifications (Christensen, 1980), but gene sequencing separated it into the Trebouxiophyceae (Senoysu et al., 2004). Our sterol analysis showed that Botryococcus together with Sphaerocystis has 5 -sterols, and thus these genera together with Chlamydomonas differ chemotaxonomically from the other green algae (Pröschold et al., 2001). However, the somatic growth and reproduction of herbivorous zooplankton was similar with Acutodesmus and Chlamydomonas diets, indicating nutritional similarity in spite of difference in sterol composition (Taipale et al., 2014). More studies of biochemical analysis of different freshwater strains together with molecular tools would be of utmost importance for more precise phytoplankton chemotaxonomy and the true phylogeny.
Multivariate statistics do not only reveal chemotaxonomic relationship of freshwater phytoplankton, but also can reveal larger groups which are able to synthesize essential fatty acids or sterols. The abundance of ALA, EPA, and DHA are important diet components defining food quality for herbivorous zooplankton and the whole food web. In our analysis NMS axis 3 (Figure 3) correlated strongly with DHA thus revealing taxa which can synthesize DHA: Dinophyceae, Diotomophyceae, Cryptophyceae, and Chrysophyceae. For sterols the dietary difference is not so clear, but in laboratory studies Martin-Creuzburg et al. (2014) showed that fucosterol, brassicasterol, ergosterol, and stigmasterol supported the somatic growth of Daphnia better than cholesterol, lathosterol, cholestanol, or 7-dehydrocholesterol, thus green algae (Chlorophyceae, Conjugatophyceae, Figure 3) are a potentially weaker quality diet for herbivorous zoonplankton. This is mostly the case in eutrophic lakes, since green algae are not so abundant in other lake types of the boreal zone (Lepistö and Rosenström, 1998). However, more studies with zooplankton are required.
Altogether, our conclusion is that fatty acid and sterol composition of freshwater phytoplankton differs at class level with some dissimilarity in sterols profiles. Actually, whereas fatty acid profiles can be used for defining phytoplankton composition at class level, additional sterol analysis can tell genus level composition in some cases. For example, ergosterol might be an excellent biomarker for Chlamydomonas in small boreal lakes to discriminate between different Chlorophyceae, where Euglena is usually absent. Additionally, 4α-methyl sterols can be used to track Dinophyceae, but also potentially to distinguish between Peridinium and Ceratium whose 4α-methyl sterols differed from each other. The usability of the fatty acid-based Bayesian mixing modeling (e.g., FASTAR, Fatty Acid Source Tracking Algorithm in R; Strandberg et al., 2015) could be improved by adding sterols, but for that more sterols profiles of freshwater strains should be analyzed.

AUTHOR CONTRIBUTIONS
ST, EP, and KV designed the study. ST, EP, and KV cultured different phytoplankton strains and ST extracted and analyzed sterols and fatty acids from phytoplankton. KV calculated major phytoplankton composition using phytoplankton database of Finnish Environment Institute. ST and MH did statistical analysis. ST wrote the paper. All authors discussed the results and commented on the manuscript.