Next Article in Journal
Towards Standardization of Data Normalization Strategies to Improve Urinary Metabolomics Studies by GC×GC-TOFMS
Previous Article in Journal
Metabolomic Profiling in Neuromyelitis Optica Spectrum Disorder Biomarker Discovery
Previous Article in Special Issue
Lipidomics Issues on Human Positive ssRNA Virus Infection: An Update
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Chasing the Major Sphingolipids on Earth: Automated Annotation of Plant Glycosyl Inositol Phospho Ceramides by Glycolipidomics

1
Department of Analytical Chemistry, Faculty of Chemistry, University of Vienna, Waehringer Str. 38, 1090 Vienna, Austria
2
Vienna Metabolomics Center (VIME), University of Vienna, Althanstraße 14, 1090 Vienna, Austria
3
Chemistry Meets Microbiology, University of Vienna, Althanstraße 14, 1090 Vienna, Austria
4
Institute of Pharmaceutical Sciences, University of Graz, Universitätsplatz 1/I, 8010 Graz, Austria
*
Author to whom correspondence should be addressed.
Metabolites 2020, 10(9), 375; https://doi.org/10.3390/metabo10090375
Submission received: 20 August 2020 / Revised: 16 September 2020 / Accepted: 16 September 2020 / Published: 19 September 2020
(This article belongs to the Special Issue Advances in Lipidomics: Biomedicine, Nutrients and Methodology)

Abstract

:
Glycosyl inositol phospho ceramides (GIPCs) are the major sphingolipids on earth, as they account for a considerable fraction of the total lipids in plants and fungi, which in turn represent a large portion of the biomass on earth. Despite their obvious importance, GIPC analysis remains challenging due to the lack of commercial standards and automated annotation software. In this work, we introduce a novel GIPC glycolipidomics workflow based on reversed-phase ultra-high pressure liquid chromatography coupled to high-resolution mass spectrometry. For the first time, automated GIPC assignment was performed using the open-source software Lipid Data Analyzer (LDA), based on platform-independent decision rules. Four different plant samples (salad, spinach, raspberry, and strawberry) were analyzed and the results revealed 64 GIPCs based on accurate mass, characteristic MS2 fragments and matching retention times. Relative quantification using lactosyl ceramide for internal standardization revealed GIPC t18:1/h24:0 as the most abundant species in all plants. Depending on the plant sample, GIPCs contained mainly amine, N-acetylamine or hydroxyl residues. Most GIPCs revealed a Hex-HexA-IPC core and contained a ceramide part with a trihydroxylated t18:0 or a t18:1 long chain base and hydroxylated fatty acid chains ranging from 16 to 26 carbon atoms in length (h16:0–h26:0). Interestingly, four GIPCs containing t18:2 were observed in the raspberry sample, which was not reported so far. The presented workflow supports the characterization of different plant samples by automatic GIPC assignment, potentially leading to the identification of new GIPCs. For the first time, automated high-throughput profiling of these complex glycolipids is possible by liquid chromatography-high-resolution tandem mass spectrometry and subsequent automated glycolipid annotation based on decision rules.

Graphical Abstract

1. Introduction

The sphingolipidome of plants contains glycosyl inositol phospho ceramides (GIPCs), glycosylceramides and ceramides, whereas sphingomyelin, globosides, sulfatides or gangliosides are absent. GIPCs were characterized as the major sphingolipid on earth due to their high abundance in plants and fungi, which comprise a large portion of the biomass of the biosphere [1]. GIPCs were first described more than 60 years ago as “phytoglycolipids” [2]. The total plant lipid content can consist of up to 40% GIPCs [3]. The structure of these plant sphingolipids has three major subunits: (1) a polar inositol containing part, (2) the sphingoid backbone with a long-chain base (amino-alcohol) linked by an amide bond to a (3) fatty acyl chain moiety [2,4]. The terms d, t and q refer to the hydroxylation state of the whole ceramide or long-chain base (LCB) moiety, ranging from two (d) to four (q) hydroxy groups. The term h denotes a hydroxylation of the fatty acyl group (i.e., the ceramide moiety q40:1 can correspond to a t18:1 LCB connected to a h22:0 fatty acyl). Di- and trihydroxylation of LCBs with t18:0, t18:1(8Z and 8E) (the main sphingoid base in some species), and d18:0, d18:1(8Z and 8E), d18:2 (4E/8Z and 4E/8E) and fatty acid components varying in chain-length, saturation and hydroxylation state (h16:0–h26:1, 20:0 to 28:0) have been reported in plant GIPCs [5,6]. Different GIPC core structures were determined from higher plants ranging from simple high-abundant A-series species with Hex-HexA-IPC and HexN(Ac)-HexA-IPC (Hex = hexose, HexA = hexuronic acid, IPC = inositol phospho ceramide, HexN = hexosamine, and HexNAc = N-acetyl hexosamine) to low abundant F-series species containing several arabinoses and hexoses [3,7]. Despite the fact that GIPCs are an integral part of the plant plasma membrane, there is still little knowledge concerning its molecular organization and the way this organization is involved in signaling processes necessary for cellular adaptation [1]. To understand the interplay of GIPCs with different enzymes and their detailed function in the plasma membrane in plants, comprehensive structural information provided by observation tools such as NMR or MS are necessary.
Even though GIPCs were discovered 60 years ago, their analysis remains challenging due to the lack of available standards, automated annotation software and reference databases. For example, CHEBI [8] does not provide any GIPCs and the comprehensive LIPID MAPS Structure Database (LMSD) contains only one GIPC (A-NH2-t18:1/h24:0) [9]. As GIPCs consist of a sugar head group linked to a lipid subunit causing amphiphilic properties, they are neither well covered by common glycomics nor lipidomics workflows. Consequently, specialized glycolipidomics analysis strategies are required, e.g., applying a mixture of 2-propanol (IPA), hexane and water [10]. The combination of liquid chromatography and mass spectrometry (LC-MS) has been used due to its unprecedented potential to annotate GIPCs by m/z, retention time and fragmentation pattern [7,11]. Unambiguous GIPC identification requires both retention time evaluation and detection of structural subunits by tandem mass spectrometry (MS2), due to the absence of commercial standards. Most GIPC LC-MS-based analysis workflows were performed almost a decade ago by electrospray ionization followed by analysis with low resolution mass spectrometers (QQQ, QTRAP) [7,11]. Meanwhile, high-resolution mass spectrometers (such as TOF, orbitrap, FTICR) have been established with up to 1,000,000 resolution enabling GIPC analysis by accurate mass [12]. Additionally, ultra-high pressure liquid chromatography (up to 1500 bar) with sub 2-µm particles provides high chromatographic resolution and excellent sensitivity. Up to now, GIPC analysis has been performed by tedious manual annotation and curation [1,7,12,13] and expert knowledge was necessary to interpret glycosphingolipid tandem mass spectrometry fragmentation patterns [14,15,16]. The instrumentation advancements of the recent years paved the way for automated high-throughput GIPC analysis. In this work, a variety of plants, i.e., iceberg lettuce (Lactuca sativa var. capitata nidus tenerimma), deep frozen spinach (Spinacia oleracea), raspberries (Rubus idaeus), and strawberries (Fragaria) were analyzed by the combination of reversed-phase (RP) ultra-high pressure liquid chromatography (UHPLC) and high-resolution mass spectrometry (HRMS). For the first time, automated GIPC annotation will be performed using the Lipid Data Analyzer (LDA) and platform-independent decision rules [17].

2. Results

Here we describe a novel workflow by RP-HRMS/MS using the open-source program LDA [17] for automated GIPC assignment. Method development considerations and guidelines for the automated structural analysis of GIPCs are provided. Finally, we test the developed glycolipidomics workflow for different plant samples, leading to a reference database of GIPCs, including fragmentation and retention time information.

2.1. Method Development for Automated GIPC Assignment

GIPCs were extracted by a mixture of IPA, n-hexane and water [18]. So far, most LC-MS-based GIPC chromatographic separations relied on the use of tetrahydrofuran (THF) containing solvents [7,11,12,13]. However, the usage of THF in the eluent system has some drawbacks: (1) it is aprotic and cannot donate a proton; thus, for ionization, pairing with a protic solvent (usually water) is necessary; (2) it can attack tubing (especially PEEK tubing); (3) it tends to polymerize (usually in APCI mode); and (4) it is highly flammable. In order to avoid the use of THF, we developed a novel GIPC method based on RP-HRMS/MS, facilitating a 30 min isopropanol gradient (detailed information can be found in the Materials and Methods Section 4.3). GIPC detection was performed using both negative and positive electrospray ionization and high-resolution Orbitrap MS (see Materials and Methods Section 4.4). Importantly, GIPC analysis requires relatively high RF voltages (S-lens RF level of 45) to ensure efficient transport of medium size glycolipids in the mass spectrometer. Figure 1 shows the extracted ion chromatogram of GIPCs in salad samples analyzed by RP-HRMS, based on data-dependent MS2 (ddMS2) in positive and negative ion modes. The GIPCs displayed in Figure 1 belong to the A-series (Hex(R1)-HexA-IPC) with R1 being a hydroxyl group and the ceramide portion consisting of a hydroxylated saturated fatty acyl chain attached to a t18:1 long chain base.
As no commercial standards are available, GIPC assignment has to be conducted with caution. In such a situation, the use of the equivalent carbon number model (ECN) is required [19,20]. The ECN model originates from state of the art lipidomics workflows and is based on elution orders observed in RP columns: (1) longer fatty acid chains will increase the retention time (see Figure 1) and (2) more double bonds will decrease the retention time [21] (see Table S1). To increase the level of confidence in GIPC annotation, we accepted only GIPCs that: (1) were detectable by accurate mass (±5 ppm) in MS1 at the same retention time in both positive and negative ion modes (Figure 1); (2) showed MS2 spectra with characteristic fragments for the ceramide and sugar part in at least one ion mode and; (3) fulfilled the ECN model.

2.2. Structural Elucidation and GIPC Annotation Based on MS2 Information

In this work, we introduce the first automated GIPC annotation workflow based on structural information provided by acquired MS2 spectra. Structural analysis and automated GIPC annotation was performed based on a set of in-house developed decision rules for the freely available software LDA [17,22]. As no standards were available, blank extractions (no GIPC annotations found) and GIPC annotations in salad [13] and spinach [12] reported in the literature were used to validate GIPC assignments (Figure 1, Table A1 and Table S2). Various LCBs (d18:0, d18:1, d18:2, t18:0, and t18:1) and fatty acids (FAs) (16–26) with or without hydroxylation have been reported [5,13]. Moreover, R1 in Figure 2A can either be a hydroxyl (OH), an amine (NH2) or an N-acetylamine (NAc) group, increasing the number of putative GIPCs even within a single series.
The final decision rule set was based on well-defined fragments (fragment rules) and their intensity relationships (intensity rules) (Folder S1). The characteristic fragments [IP] (m/z 259) and [IP-H2O] (m/z 241) are mandatory in negative ion mode (e.g.: Figure 2B). However, these fragments are not specific, since they are produced by other phosphoinositol-containing lipids too. Thus, for a confident identification, negative or positive ion mode fragments indicating the sugar or ceramide part have to be detected.
In the majority of cases (see level 2 annotations, Table A1 and Table S2), MS2 spectra with GIPC fragmentation patterns were detected in both negative and positive mode. Depending on the fragmentation pattern and the level of confidence [24] of the structural elucidation, GIPCs are assigned as either: (1) series-R1-hydroxylation stage-carbon number (LCB + FA)-number of double bonds (LCB + FA) if the exact ceramide composition is not known or (2) series-R1-LCB/FA. Figure 2B displays an exemplary ddMS2 spectrum of A-OH-q42:1 with m/z 1260.7237, in salad recorded in negative ion mode. The positive ion mode fragmentation pattern of the [M + H]+ precursor (m/z 1262.7389, Figure 2C) revealed further structural details, based on the identification of [W]+, [W-H2O]+ and [W-2H2O]+ fragments, indicating an A-OH-t18:1/h24:0 GIPC. Additional GIPC confirmation is possible by Z0 fragments ([Z0]+, [Z0-H2O]+) of the [M + H]+ precursor and by the sodium adduct [M + Na]+ (Figure A1), where sugar fragments are readily observable. GIPCs were annotated based on single ionization information only if (1) in negative ion mode in addition to the apparent [IP]/[IP-H2O]/[H2PO4] fragments at m/z 259, 241 and 97, other characteristic fragments were detectable e.g., [C3PO3] (m/z 596 − R1 = NH2, m/z 597 − R1 = OH, m/z 638 − R1 = NAc), [C3PO3-C1-CO2] (m/z 373) or [C3PO3-C1-CO2-H2O] (m/z 355) or (2) in positive ion mode the [IP]+ (m/z 261)/[IP + Na]+ (m/z 283) and fragments indicating the ceramide moiety (e.g., Z0) were identified by LDA. The detailed fragment information used for GIPC annotation can be found in Table S2.
GIPC annotation can be hampered by the presence of isobaric masses for qX:Y NH2 and t(X − 2):(Y − 1) NAc (where X refers to the carbon number (LCB + FA) and Y refers to the number of double bonds (LCB + FA), respectively). This may result in false positive GIPC identifications, because these classes share the same characteristic fragments m/z 241, 259, 355, 373 and 417. The correct structural elucidation is possible if additional fragments such as [C3PO3] (R1 = OH − m/z 597, R1 = NH2m/z 596, R1 = NAc − m/z 638) in negative ion mode or if LCBs in positive ion mode can be identified based on [W]+, [W-H2O]+ and [W-2H2O]+ fragments. In the ddMS2 spectra of the [M + H]+-precursor, trihydroxylated LCBs are characterized by the presence of three W fragments ([W]+, [W-H2O]+ and [W-2H2O]+), such as t18:0 (m/z 300, 282 and 264) and t18:1 (m/z 298, 280, 262), while dihydroxylated species miss the [W-2H2O]+ fragment, e.g., d18:0 (m/z 284, 266), d18:1 (m/z 282, 264) and d18:2 (m/z 280, 262). As such, both LCB hydroxylation levels can be clearly distinguished.

2.3. Analysis of Different Plant GIPCs by UHPLC-HRMS Suggesting t18:2 LCB

The novel RP-HRMS/MS and GIPC annotation workflow was used to analyze different plant samples, namely salad (Lactuca sativa var. capitata nidus tenerimma), deep frozen spinach (Spinacia oleracea), raspberries (Rubus idaeus) and strawberries (Fragaria). As glycosphingolipid analysis is not negatively impacted by alkaline hydrolysis [10], alkaline hydrolysis was performed to simplify lipid profiles by removing the phospholipid background in the unknown plant samples (strawberry and raspberry, detailed information can be found in the Materials and Methods Section 4.2.2). Figure A2 shows the RP-HRMS/MS GIPC profile for the five most abundant GIPCs determined in spinach, strawberry and raspberry samples. For the sake of clarity, the five most abundant GIPCs in salad (A-NAc-t18:1/h24:0, A-NH2-t18:1/h24:0, A-OH-t18:1 h22:0 and h24:0, A-OH-t18:0/h24:0) are not displayed in Figure A2. Irrespective of the plant sample, the species group A-R1-t18:1/h24:0 was always the most abundant one. While in spinach R1 was always N-acetylamine (A-NAc-t18:1 h22:0 to h26:0) for the five dominating GIPCs, in strawberries the major GIPCs contained a hydroxyl group as R1 (A-OH-t18:1 h23:0 to h26:0 and A-OH-t18:0/h24:0). In contrast to that, raspberries had an amine group as R1 for four out of five shown GIPCs (A-NH2-t18:0/h24:0, A-NH2-t18:1 h22:0 and h24:0, A-NH2-t18:2/h24:0 and A-OH-t18:1/h24:0), emphasizing the structural diversity of GIPCs in different plants.
By analyzing different GIPCs, the NAc, NH2 and OH-species from the A series could be detected (Figure 3A–C) with high confidence by (1) accurate determination of mass, (2) matching retention times of ion modes, (3) characteristic fragments and (4) the ECN model. We recommend checking isotopic patterns to avoid false positive hits. For a comprehensive overview of the annotated GIPCs see Table A1.
Due to the absence of commercially available GIPC standards, relative quantification of the individual species was performed using C16 lactosyl(ß) ceramide (d18:1/16:0) as the internal standard. This compound is similar in structure (sugar and ceramide moiety) and retention time (14 min). Even though lactosyl ceramide (d18:1/16:0) may be present in plants, we could not detect it in our samples, thus, making it suitable as the internal standard in our workflow. Normalization by the internal standard (area ratio) and dry weight was performed for MS1-based relative quantification by Skyline [25] (Figure 3A–C). Estimated concentrations in the nmol to µmol range per gram dry weight were observed, which is consistent with the literature [12,18].
In summary, 64 GIPCs in salad (19), spinach (8), strawberry (10) and raspberry (27) were annotated (Table A1). Ranking of the GIPC annotations was performed according to the guidelines of the metabolomics society [24,26], leading to 48 level 2 (matching accurate masses and MS2 in negative and positive mode) GIPCs, 13 level 3 (MS2 in one ion mode with matching accurate masses in both ion modes) GIPCs and 3 level 3** (matching accurate masses in both ion modes, MS2 in one ion mode but lacking information on IP fragments in positive ion mode or lacking sugar information in negative ion mode) GIPCs. The annotations found in spinach and salad are in accordance with literature [12,13]. To the best of our knowledge, this is the first report on GIPCs in strawberries and raspberries.
Within the annotated GIPCs with structural information on LCB and fatty acyl composition, t18:1 followed by t18:0 and t18:2 were the most prominent LCBs in the analyzed plants (Table A1). For GIPCs containing N-acetylamine residues t18:1 was the most abundant LCB with regard to normalized ratios per gram dry weight (Figure 3A). The same holds true for the amine or hydroxyl group containing GIPCs, with additional high abundance of t18:0 LCBs (Figure 3B,C). While in spinach solely t18:1 LCBs were detected, salad, strawberries and raspberries show more variation in terms of LCB composition with presence of both t18:1 and t18:0 (Table A1).
Interestingly, besides the expected t18:0 and t18:1 LCBs (R1 = NAc, NH2, OH), we additionally annotated four t18:2 (R1 = NAc, NH2, OH) species in raspberries (Figure 3A–C). These annotations are verified by coinciding retention times in positive ion mode (Figure A3A), detection of characteristic fragments in MS2 spectra (Figure A3B,C) and conformity with the ECN model (Table S1). However, we could not find any report in the literature of t18:2 species, which can be explained as up to now no automated GIPC annotation was possible and t18:2 GIPC species were only detected in raspberries. As no standards are available, it is difficult to prove the presence of this species and further investigation is needed. A confirmed t18:2 LCB would indicate a much higher diversity in sphingolipids than anticipated in the past. Another hint for the complex nature of GIPCs in raspberries is the additional annotation of GIPCs with di- and trihydroxylated variants compared to all other analyzed plants.
Concerning LCB and fatty acyl combinations, t18:1/h22:0 and t18:1/h24:0 showed equal annotation numbers for N-acetylamine or amine containing GIPCs (Figure 3A,B). Independent of the NH2, OH or NAc functional group, only two odd chain fatty acids (h23:0, h25:0) were detected and no fatty acids with a length from 17 to 21 carbon atoms were found. For the hydroxyl group containing GIPC variants, the combinations t18:1/h22:0, t18:1/h23:0, t18:0/h24:0, t18:1/h24:0 and t18:1/h25:0 were found in equal annotation numbers. Overall, plant GIPCs with a combination of t18:1 LCB and a h24:0 fatty acyl moiety were the most abundant ones in terms of normalized ratios per gram dry weight (Figure 3A–C, Table A1).

3. Discussion

GIPCs are the major sphingolipids on earth [1]. Hence, it is important to understand their function and distribution in plants and fungi. However, GIPC analysis remains extremely challenging, as tailored extraction strategies for this glycolipid class are necessary. GIPC analysis is in its infancy due to the lack of standards and databases. In this work, we present the first automated high-throughput GIPC annotation workflow which is based on RP-HRMS/MS. By using a novel 30 min gradient based on isopropanol with a reversed-phase column, packed with sub 2-µm particles, fast GIPC analysis was possible at the same time avoiding standard eluent use of tetrahydrofuran. Four different plant samples were analyzed. For salad and spinach, literature information has been available [12,13], while for raspberry and strawberry, GIPC profiles were completely uncharacterized. Using strict filtering by (1) accurate mass determination (±5 ppm) with matching retention times for both ion modes in MS1, (2) MS2 spectra with characteristic fragments and (3) expected retention time series, we produced a database of 64 GIPCs (Table A1). As no GIPC standards are available, only GIPC annotation hits with level 2 and 3 confidence [24] were possible. The most prominent MS2 fragments for GIPCs are [IP] fragments in both ion modes ([H]: m/z 241, 259; [H]+: m/z 261; [Na]+: m/z 283). However, additional sugar or ceramide fragments are essential for correct GIPC annotation. The high MS2 mass range coverage (m/z 65 to 2500) provided by the Orbitrap was beneficial to determine GIPC low mass fragments such as m/z 79 [PO3] or 97 [H2PO4], besides high mass precursors such as 1261 [M − H] (Figure 2).
Relative quantification with the internal standard lactosyl ceramide revealed GIPC t18:1/h24:0 as the most abundant species, independent of the plant sample. Depending on the plant sample, GIPCs contained mainly amine, N-acetyl or hydroxyl residues. Most GIPCs showed a Hex-HexA-IPC core with a trihydroxylated t18:0 or t18:1 long-chain base ceramide part a and hydroxylated fatty acid chains ranging from h16:0 to h26:0. Interestingly, in raspberry, four GIPCs contained t18:2, which was not reported so far. This finding would suggest the existence of more complex sphingolipid species in nature than previously anticipated. Further analysis by orthogonal methods such as NMR, GC-MS or IMS and available GIPC standards would be necessary to confirm the presence of the t18:2 GIPC group. Different analytical strategies could also resolve potential isomeric species and provide comprehensive details on the sugar moiety present in GIPCs. Nevertheless, this example shows the power of this workflow to detect promising novel GIPC candidates in an automated fashion. In order to support LC-MS-based GIPC analysis in general, we provide the mass lists for GIPCs in positive and negative ion modes (Tables S3 and S4), as well as the fragmentation rules (Folder S1) for setting up the automated GIPC analysis by Lipid Data Analyzer. Even though we confirmed the GIPCs exclusively from the A-series, the presented strategy is also suitable to determine less or more complex GIPC series, such as 0, B, C, D, E and F. However, extended analytical workflows (e.g., multi-stage fragmentation/MSn) and additional software method development might be necessary. Precursor mass lists for positive ([M + H]+) and negative ([M − H]) ion modes comprising series 0–F, LCBs d18:0, d18:1, d18:2, t18:0 and t18:1 and fatty acyls h15:0–h26:0, h15:1–h26:1 and n20:0–n28:0 (n = non-hydroxylated), as reported in the literature [5,13], can be found in Tables S5 and S6. In general, we believe that LC-HRMS/MSn combined with automated annotation based on decision rules will pave the way for more complex glycolipidomics profiling.

4. Materials and Methods

4.1. Material

The plant material used was derived from salad (Lactuca sativa var. capitata nidus tenerimma), deep frozen spinach (Spinacia oleracea), raspberries (Rubus idaeus) and strawberries (Fragaria). (A more detailed description of plant samples can be found in Table A2.)
All chemicals were of LC-MS grade. Acetonitrile (ACN), methanol (MeOH), IPA and water were bought from Honeywell (Offenbach, Germany) and n-hexane was bought from VWR (Vienna, Austria). Butylated hydroxytoluene (BHT) was purchased from Sigma-Aldrich (Vienna, Austria), ammonium formate (AF) from Sigma-Aldrich (Vienna, Austria) and formic acid from VWR (Vienna, Austria). C16 Lactosyl(ß) Ceramide (d18:1/16:0) (D-lactosyl-ß-1,1’ N-palmitoyl-D-erythro-sphingosine) was purchased from Avanti Polar Lipids, Inc. (Alabaster, Albama, USA), was used as internal standard (IS) and dissolved in an appropriate amount of IPA to achieve a concentration of 100 µM.

4.2. Sample Preparation

Salad was manually cut into small pieces before being weighed into falcon tubes (50 mL, VWR, Vienna, Austria) using a CPA225D balance (Sartorius, Vienna, Austria). Raspberries and strawberries (whole fruits) were homogenized with a hand blender (Tefal/SEB, Ecully, France). Raspberries, strawberries and deep-frozen homogenized spinach were directly weighed into 10 mL glass vials (more details can be found in Table A2). In order to prevent potential oxidation of lipids, 3 mL of an approximately 0.01% BHT solution in IPA were added and samples were mixed. Subsequently 30 µL IS were spiked into all samples except for one replicate (to test for potential IS presence in plants). Salad samples were homogenized using an ultra-turax (miccra d-1, Heitersheim, Germany) which was cleaned with 70% IPA and dried between the samples. In order to inhibit lipase activity, all samples were incubated at 75 °C for 30 min under constant shaking [27]. The warm salad samples were subsequently transferred into glass vials. The following sections provide a detailed overview of the extraction strategies that were applied.

4.2.1. One-Phase Extraction

The extraction of GIPCs from salad and spinach was performed as previously reported [18] using a mixture of IPA, n-hexane and water. Amounts of 3.47 mL IPA, 0.6 mL n-hexane and 1.93 mL water were added to the salad and spinach samples. In order to ensure sufficient accessibility of the plant material, samples were vortexed and manually shaken prior to incubation at 60 °C for 15 min under constant shaking.

4.2.2. One-Phase Extraction Combined with Alkaline Hydrolysis

To avoid the occurrence of glycerophospholipids, which might reduce GIPC ionization efficiency and lead to potential false identifications, alkaline hydrolysis was applied for the raspberry and strawberry samples, using an adapted workflow [28]. After incubating the plant material with the BHT solution for 30 min at 75 °C under constant shaking, 3.47 mL IPA and 0.6 mL n-hexane were added. Samples were vortexed and put on a shaker for 15 min at 60 °C. As soon as the samples had reached room temperature 707 µL 1 M KOH in MeOH was added and the solution was vortexed. After shaking the samples for 2 h at 37 °C, they were left at room temperature. Subsequently 100% formic acid was added until a pH of ~6–7 was reached and 1.93 mL water was added before repeating the incubation step.

4.2.3. Centrifugation, Drying and Reconstitution

Irrespective of the extraction strategy, the warm samples were centrifuged at 1000 rpm for 10 min at 4 °C and the supernatant was transferred into a separate glass vial. The solvent was evaporated to dryness overnight in a Genevac EZ-2 Series Personal Evaporator (SP Scientific, Ipswich, UK) and the dried residue was reconstituted in 2 mL IPA:H2O (65:35) [13]. Samples were vortexed prior and after ultrasonication at 30 °C for 15 min. Subsequently, 500 µL of this solution was filtered directly into HPLC vials through a ClariStep filter (Sartorius, Vienna, Austria). Pools were prepared separately for each plant by pipetting 50 µL of each biological replicate into a separate HPLC vial. A quality control pool was prepared by combining 30 µL of the pooled samples.

4.3. Reversed-Phase Chromatography

Liquid chromatography was performed using a C18 Acquity UHPLC HSS T3 reversed phase column (2.1 × 150 mm, 100 Å, 1.8 µm, Waters, Vienna, Austria) equipped with a VanGuard Pre-column (2.1 × 5 mm, 100 Å, 1.8 µm, Waters, Vienna, Austria) at a column temperature of 40 °C. The flow rate was 0.25 mL/min and the backpressure was 460 bar at the starting conditions. Gradient elution with a total runtime of 30 min was performed using the solvent A: ACN:H2O (3:2, v/v) and the solvent B: IPA:ACN (9:1, v/v), both of which contained 0.1% formic acid and 10 mM ammonium formate.
The gradient can be described as follows: 0–2 min 30% B, 2–3 min ramp to 55% B, 3–17 min ramp to 67% B, 17–22 min ramp to 100% B, 22–26 min 100% B, followed by an equilibration step from 26 to 30 min using 30% B. A Vanquish Duo UHPLC system (Thermo Fisher Scientific, Germering, Germany) was used and injections were performed with an autosampler. An injection volume of 10 µL was chosen and the injector needle was flushed with 75% IPA and 1% formic acid in between the injections.

4.4. High-Resolution Mass Spectrometry

The LC system was coupled to a Q Exactive HF (Thermo Fisher Scientific, Bremen, Germany) high resolution mass spectrometer, applying a HESI ion source with an S-lens RF level of 45. Measurements were carried out in positive and negative modes using different parameters. The following settings were applied in positive mode: spray voltage: 3.5 kV, capillary temperature 220 °C, sheath gas flow rate: 30, and auxiliary flow rate: 5. In negative mode parameters were adapted as follows: spray voltage: 2.8 kV, capillary temperature 250 °C, sheath gas flow rate: 35 (a.u.), and auxiliary flow rate: 10 (a.u.). The top 10 data-dependent MS2 spectra were obtained at a scan range of 500 to 3000 m/z with HCD using normalized collision energies of 35 (+35 in positive ion mode, −35 in negative mode), an MS1 resolution of 15,000 or 30,000 with an AGC target of 1e6 and MS2 resolution of 15,000 with an AGC target of 1e5. MS2 spectra were acquired based on an inclusion list (“do not pick others” option) containing the GIPC series 0–F (m/z values were calculated using enviPat Web 2.4 [29]). A more comprehensive picture of the GIPC composition of the analyzed plant material was obtained using several rounds of automatically generated exclusions lists for the sample pools [30].

4.5. Data Analysis

The GIPC assignment was performed using LDA (version 2.8.0) [17]; corresponding settings (Table A3), mass lists (Tables S3 and S4) and decision rule sets for series A (Folder S1) can be found in the Appendix A and Supplementary Materials. The correct GIPC annotation was ensured by a manual inspection of the results. MS1-based relative quantification of annotated GIPCs was performed with Skyline [25]. Total areas were divided by the corresponding calculated dry weights and areas of the IS, resulting in normalized ratios per g dry weight, of which the average was taken based on the number of replicates (3 for salad and spinach, 4 for strawberries and raspberries). More information can be found in Appendix B.

Supplementary Materials

The following are available online at https://www.mdpi.com/2218-1989/10/9/375/s1, Table S1. Application of the ECN model to ensure correct GIPC annotation, Table S2. List of annotated GIPCs including all detected MS2 fragments, Table S3. LDA mass list used for automated annotation of the series A GIPCs in positive mode, Table S4. LDA mass list used for automated annotation of the series A GIPCs in negative mode, Table S5. List of [M + H]+ precursors comprising the GIPC series 0–F, Table S6. List of [M − H] precursors comprising the GIPC series 0–F, Folder S1. Fragmentation rules for GIPC analysis by LDA.

Author Contributions

Conceptualization, E.R.; methodology, L.P. and E.R.; software, L.P., S.S., N.T. and J.H.; validation, L.P., E.R. and J.H.; formal analysis, L.P., S.S., N.T. and E.R.; investigation, L.P. and E.R.; resources, E.R. and G.K.; data curation, L.P. and E.R.; writing—original draft preparation, E.R. and L.P.; writing—review and editing, E.R., L.P., J.H. and G.K.; visualization, L.P.; supervision, E.R.; project administration, E.R. and G.K. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Acknowledgments

This work was supported by the University of Vienna, the Faculty of Chemistry, the Vienna Metabolomics Center (VIME; http://metabolomics.univie.ac.at/), the research platform Chemistry Meets Microbiology and the Mass Spectrometry Centre of the University of Vienna. The authors thank (1) the Department of Food Chemistry and Toxicology (University of Vienna) for sharing their sample preparation equipment, (2) Sophia Mundigler for the preparation of a literature database for glycolipid fragmentation (3) Martin Schaier for graphical abstract support, as well as (4) all members of the Koellensperger lab (University of Vienna) for continuous support.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Table A1. Overview of GIPCs annotated in salad, spinach, strawberries and raspberries. The precursor ion m/z and retention times are listed as provided by the LDA display results function. In cases where only matching retention times but no m/z were explicitly shown by the LDA, because only one MS2 spectrum was annotated (level 3, level 3**), corresponding values (marked with an asterisk*) were manually assigned at the peak maximum using Thermo Scientific FreeStyle. In the Level column the levels of identification are listed. For all annotated GIPCs, accurate mass and retention times were observed. Their level of identification depends on the ddMS2 spectra. Level 2 annotation is based on the ddMS2 spectra in both ion modes. For level 3, ddMS2 spectra with characteristic fragments could only be detected in one ion mode. Putative hits, which cannot be annotated with such high confidence, because they were only observed with a sugar fragment in positive mode, but did not show the [IP]+/[IP+Na]+ fragment, are listed as level 3** at the end of the table and were not included in Figure 3.
Table A1. Overview of GIPCs annotated in salad, spinach, strawberries and raspberries. The precursor ion m/z and retention times are listed as provided by the LDA display results function. In cases where only matching retention times but no m/z were explicitly shown by the LDA, because only one MS2 spectrum was annotated (level 3, level 3**), corresponding values (marked with an asterisk*) were manually assigned at the peak maximum using Thermo Scientific FreeStyle. In the Level column the levels of identification are listed. For all annotated GIPCs, accurate mass and retention times were observed. Their level of identification depends on the ddMS2 spectra. Level 2 annotation is based on the ddMS2 spectra in both ion modes. For level 3, ddMS2 spectra with characteristic fragments could only be detected in one ion mode. Putative hits, which cannot be annotated with such high confidence, because they were only observed with a sugar fragment in positive mode, but did not show the [IP]+/[IP+Na]+ fragment, are listed as level 3** at the end of the table and were not included in Figure 3.
CompositionPlantm/z
[M − H]
m/z
[M + H]+
Rt_neg [min]Rt_pos [min]LevelNormalized Ratio/g dw [g−1]CV [%]
A-NAc-t18:1/h16:0Spinach1189.62461191.63888.008.0023.1525
A-NAc-t18:1/h16:0Salad1189.62461191.63978.018.0022.1028
A-NAc-t18:1/h22:0Spinach1273.71791275.735914.2514.25227.4319
A-NAc-t18:1/h22:0Salad1273.71931275.734414.1914.2623.5744
A-NAc-t18:1/h22:0Raspberry1273.71721275.734114.2114.1920.946
A-NAc-t18:1/h23:0Spinach1287.73251289.750415.6415.65211.8321
A-NAc-t18:0/h24:0Raspberry1303.76481305.781317.9418.1722.134
A-NAc-t18:1/h24:0Spinach1301.75021303.765716.9816.98270.8921
A-NAc-t18:1/h24:0Salad1301.74981303.763716.9816.94217.2332
A-NAc-t18:1/h24:0Raspberry1301.74911303.764916.9616.9423.405
A-NAc-t18:2/h24:0Raspberry1299.73221301.74914.9814.9721.273
A-NAc-t18:1/h25:0Spinach1315.76521317.781118.4118.41211.0820
A-NAc-t18:1/h26:0Spinach1329.78111331.79719.7519.75220.5416
A-NH2-t34:0Raspberry1133.63461135.650610.9110.9221.096
A-NH2-d34:1Raspberry1115.6248 *1117.641510.8910.9031.857
A-NH2-t34:1Raspberry1131.6196 *1133.636510.2710.2730.758
A-NH2-t18:1/h16:0Salad1147.61481149.63028.648.5226.465
A-NH2-q40:0Salad1233.7225 *1235.737216.8216.7633.8727
A-NH2-t18:0/h22:0Raspberry1233.72421235.738916.6416.5523.337
A-NH2-t18:1/h22:0Spinach1231.7058 *1233.722915.5615.5231.8718
A-NH2-t18:1/h22:0Salad1231.70831233.722615.5115.49210.4511
A-NH2-t18:1/h22:0Raspberry1231.70821233.723415.4115.2523.877
A-NH2-t18:2/h22:0Raspberry1229.6945 *1231.707613.4213.4231.014
A-NH2-t18:1/h23:0Salad1245.7242 *1247.740016.7816.8734.0415
A-NH2-t18:1/h23:0Raspberry1245.72481247.738816.8916.8522.797
A-NH2-t18:0/h24:0Raspberry1261.75511263.770419.6219.6426.929
A-NH2-t18:0/h24:0Salad1261.7531*1263.770119.8219.5333.8027
A-NH2-t18:1/h24:0Spinach1259.73911261.756518.3218.3025.0618
A-NH2-t18:1/h24:0Salad1259.74001261.755418.5218.44232.006
A-NH2-t18:1/h24:0Raspberry1259.74001261.754318.4618.42212.288
A-NH2-t18:2/h24:0Raspberry1257.72311259.739715.9816.2524.228
A-NH2-t18:1/h25:0Raspberry1273.7539 *1275.769219.5519.7831.5511
A-NH2-q43:2Raspberry1271.7419 *1273.757117.7117.7330.879
A-OH-t34:0Raspberry1134.61921136.636910.2110.2120.588
A-OH-t18:0/h16:0Salad1150.61481152.63138.848.8321.930
A-OH-t18:1/h16:0Salad1148.59841150.61388.148.13210.2716
A-OH-q34:1Strawberry1148.59821150.61448.118.1321.516
A-OH-q40:0Salad1234.70951236.7239 *15.6615.6033.002
A-OH-t18:0/h22:0Raspberry1234.70841236.724715.3615.6321.376
A-OH-t18:1/h22:0Salad1232.69211234.707114.4714.44215.1619
A-OH-t18:1/h22:0Strawberry1232.69211234.708314.2514.4323.692
A-OH-t18:1/h22:0Raspberry1232.69241234.705514.4314.4321.446
A-OH-q41:0Strawberry1248.72321250.7395 *17.0116.9830.832
A-OH-t18:1/h23:0Raspberry1246.70761248.7822115.8215.7721.467
A-OH-t18:1/h23:0Salad1246.70791248.720915.8515.8424.4518
A-OH-t18:1/h23:0Strawberry1246.70761248.723815.5315.7628.201
A-OH-t18:0/h24:0Salad1262.73971264.756418.5518.50212.304
A-OH-t18:0/h24:0Raspberry1262.73851264.75418.5118.5322.586
A-OH-t18:0/h24:0Strawberry1262.73871264.755318.4418.3224.531
A-OH-t18:1/h24:0Salad1260.72371262.738917.3017.26250.7714
A-OH-t18:1/h24:0Strawberry1260.72271262.739117.2617.00222.152
A-OH-t18:1/h24:0Raspberry1260.72401262.738817.2517.2323.896
A-OH-q42:2Strawberry1258.70731260.7238 *15.2115.2231.222
A-OH-t18:1/h24:1Salad1258.70671260.724314.3614.3322.0410
A-OH-t18:2/h24:0Raspberry1258.70821260.723415.2315.2321.475
A-OH-t18:0/h25:0Strawberry1276.75421278.768719.8019.8021.335
A-OH-t18:1/h25:0Salad1274.73871276.753318.7618.7523.7211
A-OH-t18:1/h25:0Strawberry1274.73711276.755718.4218.4825.832
A-OH-t18:1/h25:0Raspberry1274.73811276.754818.7118.7320.849
A-OH-t18:1/h26:0Salad1288.75531290.770620.0120.0128.6518
A-OH-t18:1/h26:0Strawberry1288.7540 *1290.769719.7619.9435.492
A-NH2-d18:2/16:0Raspberry1113.6100 *1115.62559.899.893 **0.5710
A-NH2-q41:0Raspberry1247.7383 *1249.754418.1618.173 **2.037
A-NH2-q41:2Raspberry1243.7076 *1245.723214.7714.803 **0.877
Figure A1. The ddMS2 spectrum of the [M + Na]+ of GIPC A-OH-t18:1/h24:0 (m/z 1284.7193, Rt 17.30 min), measured in positive ion mode, showing the characteristic [IP + Na]+ and additional sugar fragments.
Figure A1. The ddMS2 spectrum of the [M + Na]+ of GIPC A-OH-t18:1/h24:0 (m/z 1284.7193, Rt 17.30 min), measured in positive ion mode, showing the characteristic [IP + Na]+ and additional sugar fragments.
Metabolites 10 00375 g0a1
Figure A2. Comparison of the RP-HRMS/MS GIPC profiles in spinach (green), strawberry (rose) and raspberry (dark-red), showing the five most abundant GIPCs found in each plant sample measured in positive ion mode (detailed information can be found in Table A1).
Figure A2. Comparison of the RP-HRMS/MS GIPC profiles in spinach (green), strawberry (rose) and raspberry (dark-red), showing the five most abundant GIPCs found in each plant sample measured in positive ion mode (detailed information can be found in Table A1).
Metabolites 10 00375 g0a2
Figure A3. (A) Extracted ion chromatograms of GIPCs in positive ion mode, having a t18:2 LCB annotation on levels 2 and 3 in raspberries (A-NAc-t18:2/h24:0, A-NH2-t18:2 h22:0 and h24:0, as well as A-OH-t18:2/h24:0). The ddMS2 spectra of GIPCs A-NH2-t18:2/h24:0 at the retenion time of 16.1 min (B) A-OH-t18:2/h24:0 at the retention time of 15.2 min and (C) in negative and positive mode, showing characteristic fragments.
Figure A3. (A) Extracted ion chromatograms of GIPCs in positive ion mode, having a t18:2 LCB annotation on levels 2 and 3 in raspberries (A-NAc-t18:2/h24:0, A-NH2-t18:2 h22:0 and h24:0, as well as A-OH-t18:2/h24:0). The ddMS2 spectra of GIPCs A-NH2-t18:2/h24:0 at the retenion time of 16.1 min (B) A-OH-t18:2/h24:0 at the retention time of 15.2 min and (C) in negative and positive mode, showing characteristic fragments.
Metabolites 10 00375 g0a3aMetabolites 10 00375 g0a3bMetabolites 10 00375 g0a3c
Table A2. The description of the plant samples, including plant species, origin, number of biological replicates, and the average fresh- and dry weights [g]. The extraction of GIPCs from strawberries and raspberries was performed one day after the collection was performed (28 June 2020).
Table A2. The description of the plant samples, including plant species, origin, number of biological replicates, and the average fresh- and dry weights [g]. The extraction of GIPCs from strawberries and raspberries was performed one day after the collection was performed (28 June 2020).
Plant SpeciesOriginReplicatesFresh Weight [g]Dry Weight [g]
SaladLocal supermarket4~1~0.04
SpinachLocal supermarket4~1.2~0.08
Strawberries47°58′ N, 16°6′ O5~1.1~0.10
Raspberries47°58′ N, 16°6′ O5~1.1~0.16
Table A3. The exemplary LDA parameters and settings used for automated GIPC annotation in negative ion mode.
Table A3. The exemplary LDA parameters and settings used for automated GIPC annotation in negative ion mode.
ParameterSetting
Time before tol.1 min
Time after tol.1 min
Rel. Base-peak cutoff0.1‰
Rt-shift0.0 min
Isotopic quantitation of _ isotopes where _ isotopic peak(s) have to match2, 1
Find molecules where retention time is unknownyes
LDA-version2.8.0
machineNameOrbiTrap_exactive
neutronMass1.005
coarseChromMzTolerance0.015
MS2true
basePeakCutoff0.1
massShift0.0
threeDViewerDefaultTimeResolution2
threeDViewerDefaultMZResolution0.005
ms2PrecursorTolerance0.013
ms2MzTolerance0.02
ms2MinIntsForNoiseRemoval100
ms2IsobarSCExclusionRatio0.01
ms2IsobarSCFarExclusionRatio0.1
ms2IsobaricOtherRtDifference2.0
chainCutoffValue0.01
ms2ChromMultiplicationFactorForInt10
threeDViewerMs2DefaultTimeResolution1
threeDViewerMs2DefaultMZResolution1
maxFileSizeForChromTranslationAtOnce500
chromMultiplicationFactorForInt1000
chromLowestResolution1
chromSmoothRange8.0
chromSmoothRepeats4
use3Dtrue
isotopeCorrectionfalse
removeFromOtherIsotopestrue
respectIsotopicDistributiontrue
checkChainLabelCombinationFromSpeciesNamefalse
useNoiseCutofftrue
noiseCutoffDeviationValue2.0
scanStep2
profileMzRangeExtraction0.05
profileTimeTolerance5.0
profileIntThreshold5.0
broaderProfileTimeTolerance3.0
profileSmoothRange0.0025
profileSmoothRepeats1
profileMeanSmoothRepeats2
profileMzMinRange0.002
profileSteepnessChange11.5
profileSteepnessChange21.8
profileIntensityCutoff10.15
profileIntensityCutoff20.2
profileGeneralIntCutoff0.03
profilePeakAcceptanceRange0.012
profileSmoothingCorrection0.0
profileMaxRange0.03
smallChromMzRange0.004
smallChromSmoothRepeats3
smallChromMeanSmoothRepeats0
smallChromSmoothRange2.0
smallChromIntensityCutoff0.03
broadChromSmoothRepeats5
broadChromMeanSmoothRepeats0
broadChromSmoothRange2
broadChromIntensityCutoff0.0
broadChromSteepnessChangeNoSmall1.33
broadChromIntensityCutoffNoSmall0.05
finalProbeTimeCompTolerance0.1
finalProbeMzCompTolerance5.0E-4
overlapDistanceDeviationFactor1.5
overlapPossibleIntensityThreshold0.15
overlapSureIntensityThreshold0.7
overlapPeakDistanceDivisor3.0
overlapFullDistanceDivisor6.0
peakDiscardingAreaFactor1000
isotopeInBetweenTime30
isoInBetweenAreaFactor3.0
isoNearNormalProbeTime30
relativeAreaCutoff0.05
relativeFarAreaCutoff0.05
relativeFarAreaTimeSpace30
relativeIsoInBetweenCutoff0.5
isoInBetweenMaxTimeDistance300
twinPeakMzTolerance0.01
closePeakTimeTolerance10
twinInBetweenCutoff0.95
unionInBetweenCutoff0.8
sparseDatafalse

Appendix B

Automated GIPC annotation was performed using LDA (version 2.8.0) [17] with the settings provided in Table A3. The mass-to-charge ratios included in the mass lists (see Tables S3 and S4) were calculated separately for negative and positive ion modes, with enviPat Web 2.4 [29] and decision rules (see Folder S1) were created based on fragments reported in the literature [12,13]. Please note that the raw data acquired in negative ion mode has to be analyzed using the mass list of Table S4 and the fragmentation rules ending with ‘-H.frag’, while for positive mode the mass list of Table S3 and corresponding fragmentation rules (‘H.frag’ and ‘Na.frag’) should be used. Further information on working with the LDA can be found in [31].

References

  1. Gronnier, J.; Germain, V.; Gouguet, P.; Cacas, J.-L.; Mongrand, S. GIPC: Glycosyl Inositol Phospho Ceramides, the major sphingolipids on earth. Plant Signal. Behav. 2016, 11, e1152438. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Carter, H.E.; Gigg, R.H.; Laws, J.H. Structure of Phytoglyolipide. J. Biol. Chem. 1958, 233, 1309–1314. [Google Scholar] [PubMed]
  3. Cacas, J.; Buré, C.; Grosjean, K.; Gerbeau-Pissot, P.; Lherminier, J.; Rombouts, Y.; Maes, E.; Bossard, C.; Gronnier, J.; Furt, F.; et al. Revisiting Plant Plasma Membrane Lipids in Tobacco: A Focus on Sphingolipids. Plant Physiol. 2016, 170, 367–384. [Google Scholar] [CrossRef] [PubMed]
  4. Carter, H.E.; Strobach, D.R.; Hawthorne, J.N. Biochemistry of the Sphingolipids. XVIII. Complete Structure of Tetrasaccharide Phytoglycolipid. Biochemistry 1969, 8, 383–388. [Google Scholar] [CrossRef] [PubMed]
  5. Buré, C.; Cacas, J.; Mongrand, S.; Schmitter, J.-M. Characterization of glycosyl inositol phosphoryl ceramides from plants and fungi by mass spectrometry. Anal. Bioanal. Chem. 2014, 406, 995–1010. [Google Scholar] [CrossRef]
  6. LIPID MAPS® Lipidomics Gateway. Available online: https://lipidmaps.org/resources/lipidweb/index.php?page=lipids/sphingo/glyP_ino/index.htm (accessed on 4 August 2020).
  7. Buré, C.; Cacas, J.-L.; Wang, F.; Gaudin, K.; Domergue, F.; Mongrand, S.; Schmitter, J.-M. Fast screening of highly glycosylated plant sphingolipids by tandem mass spectrometry. Rapid Commun. Mass Spectrom. 2011, 25, 3131–3145. [Google Scholar] [CrossRef]
  8. Hastings, J.; Owen, G.; Dekker, A.; Ennis, M.; Kale, N.; Muthukrishnan, V.; Turner, S.; Swainston, N.; Mendes, P.; Steinbeck, C. ChEBI in 2016: Improved services and an expanding collection of metabolites. Nucleic Acids Res. 2016, 44, D1214–D1219. [Google Scholar] [CrossRef]
  9. Sud, M.; Fahy, E.; Cotter, D.; Brown, A.; Dennis, E.A.; Glass, C.K.; Merrill, A.H.; Murphy, R.C.; Raetz, C.R.H.; Russell, D.W.; et al. LMSD: LIPID MAPS structure database. Nucleic Acids Res. 2007, 35, 527–532. [Google Scholar] [CrossRef] [Green Version]
  10. Barrientos, R.C.; Zhang, Q. Recent advances in the mass spectrometric analysis of glycosphingolipidome—A review. Anal. Chim. Acta 2020. [Google Scholar] [CrossRef]
  11. Markham, J.E.; Jaworski, J.G. Rapid measurement of sphingolipids from Arabidopsis thaliana by reversed-phase high-performance liquid chromatography coupled to electrospray ionization tandem mass spectrometry. Rapid Commun. Mass Spectrom. 2007, 21, 1304–1314. [Google Scholar] [CrossRef]
  12. Blaas, N.; Humpf, H.U. Structural profiling and quantitation of glycosyl inositol phosphoceramides in plants with fourier transform mass spectrometry. J. Agric. Food Chem. 2013, 61, 4257–4269. [Google Scholar] [CrossRef] [PubMed]
  13. Cacas, J.; Buré, C.; Furt, F.; Maalouf, J.; Badoc, A.; Cluzet, S.; Schmitter, J.; Antajan, E.; Mongrand, S. Biochemical survey of the polar head of plant glycosylinositolphosphoceramides unravels broad diversity. Phytochemistry 2013, 96, 191–200. [Google Scholar] [CrossRef] [PubMed]
  14. Domon, B.; Costello, C.E. Structure Elucidation of Glycosphingolipids and Gangliosides Using High-Performance Tandem Mass Spectrometry. Biochemistry 1988, 27, 1534–1543. [Google Scholar] [CrossRef] [PubMed]
  15. Domon, B.; Costello, C.E. A systematic nomenclature for carbohydrate fragmentations in FAB-MS/MS spectra of glycoconjugates. Glycoconj. J. 1988, 5, 397–409. [Google Scholar] [CrossRef]
  16. Hsu, F.; Turk, J. Glycosphingolipids as Lithiated Adducts by Dissociation on a Triple Stage Quadrupole Instrument. J. Am. Soc. Mass Spectrom. 2001, 12, 61–79. [Google Scholar] [CrossRef] [Green Version]
  17. Hartler, J.; Triebl, A.; Ziegl, A.; Trötzmüller, M.; Rechberger, G.N.; Zeleznik, O.A.; Zierler, K.A.; Torta, F.; Cazenave-Gassiot, A.; Wenk, M.R.; et al. Deciphering lipid structures based on platform-independent decision rules. Nat. Methods 2017, 14, 1171–1174. [Google Scholar] [CrossRef]
  18. Markham, J.E.; Li, J.; Cahoon, E.B.; Jaworski, J.G. Separation and Identification of Major Plant Sphingolipid Classes from Leaves. J. Biol. Chem. 2006, 281, 22684–22694. [Google Scholar] [CrossRef] [Green Version]
  19. Dugo, P.; Cacciola, F.; Kumm, T.; Dugo, G.; Mondello, L. Comprehensive multidimensional liquid chromatography: Theory and applications. J. Chromatogr. A 2008, 1184, 353–368. [Google Scholar] [CrossRef]
  20. Lísa, M.; Holčapek, M. Triacylglycerols profiling in plant oils important in food industry, dietetics and cosmetics using high-performance liquid chromatography–atmospheric pressure chemical ionization mass spectrometry. J. Chromatogr. A 2008, 1198–1199, 115–130. [Google Scholar] [CrossRef]
  21. Ovčačíková, M.; Lísa, M.; Cífková, E.; Holčapek, M. Retention behavior of lipids in reversed-phase ultrahigh-performance liquid chromatography–electrospray ionization mass spectrometry. J. Chromatogr. A 2016, 1450, 76–85. [Google Scholar] [CrossRef]
  22. Hartler, J.; Trötzmüller, M.; Chitraju, C.; Spener, F.; Köfeler, H.C.; Thallinger, G.G. Lipid Data Analyzer: Unattended identification and quantitation of lipids in LC-MS data. Bioinformatics 2011, 27, 572–577. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. LIPID MAPS® Lipidomics Gateway. Available online: https://www.lipidmaps.org/data/LMSDRecord.php?LMID=LMSP03030005 (accessed on 15 March 2020).
  24. Salek, R.M.; Steinbeck, C.; Viant, M.R.; Goodacre, R.; Dunn, W.B. The role of reporting standards for metabolite annotation and identification in metabolomic studies. Gigascience 2013, 2, 2047-217X-2-13. [Google Scholar] [CrossRef] [PubMed]
  25. Adams, K.J.; Pratt, B.; Bose, N.; Dubois, L.G.; St. John-Williams, L.; Perrott, K.M.; Ky, K.; Kapahi, P.; Sharma, V.; Maccoss, M.J.; et al. Skyline for Small Molecules: A Unifying Software Package for Quantitative Metabolomics. J. Proteome Res. 2020, 19, 1447–1458. [Google Scholar] [CrossRef] [PubMed]
  26. Metabolomics Society. Available online: http://metabolomicssociety.org/ (accessed on 14 July 2020).
  27. Shiva, S.; Enninful, R.; Roth, M.R.; Tamura, P.; Jagadish, K.; Welti, R. An efficient modified method for plant leaf lipid extraction results in improved recovery of phosphatidic acid. Plant Methods 2018, 14, 1–8. [Google Scholar] [CrossRef] [Green Version]
  28. Peng, B.; Weintraub, S.T.; Coman, C.; Ponnaiyan, S.; Sharma, R.; Tews, B.; Winter, D.; Ahrends, R. A Comprehensive High-Resolution Targeted Workflow for the Deep Profiling of Sphingolipids. Anal. Chem. 2017, 89, 12480–12487. [Google Scholar] [CrossRef] [PubMed]
  29. Loos, M.; Gerber, C.; Corona, F.; Hollender, J.; Singer, H. Accelerated Isotope Fine Structure Calculation Using Pruned Transition Trees. Anal. Chem. 2015, 87, 5738–5744. [Google Scholar] [CrossRef] [PubMed]
  30. Koelmel, J.P.; Kroeger, N.M.; Gill, E.L.; Ulmer, C.Z.; Bowden, J.A.; Patterson, R.E.; Yost, R.A.; Garrett, T.J. Expanding Lipidome Coverage Using LC-MS/MS Data-Dependent Acquisition with Automated Exclusion List Generation. J. Am. Soc. Mass Spectrom. 2017, 28, 908–917. [Google Scholar] [CrossRef] [Green Version]
  31. LDA User Manual. Available online: http://genome.tugraz.at/lda2/2.6/LDA_2.6.pdf (accessed on 17 September 2020).
Figure 1. Extracted ion chromatogram of glycosyl inositol phospho ceramides (GIPCs) in salad samples analyzed by RP-HRMS/MS analysis using ddMS2 in positive (red) and negative (blue) ion modes. Assigned GIPCs belong to the Hex-HexA-IPC series with a t18:1 long-chain base (LCB) and varying chain length of the hydroxylated saturated fatty acids. Retention times coincided in positive and negative ion modes. Increasing carbon numbers result in belated elution.
Figure 1. Extracted ion chromatogram of glycosyl inositol phospho ceramides (GIPCs) in salad samples analyzed by RP-HRMS/MS analysis using ddMS2 in positive (red) and negative (blue) ion modes. Assigned GIPCs belong to the Hex-HexA-IPC series with a t18:1 long-chain base (LCB) and varying chain length of the hydroxylated saturated fatty acids. Retention times coincided in positive and negative ion modes. Increasing carbon numbers result in belated elution.
Metabolites 10 00375 g001
Figure 2. Overview of the GIPC fragmentation for the example of GIPC A-OH-t18:1/h24:0 in salad: (A) The fragment assignment of GIPC A-OH-t18:1/h24:0 (adapted from [23]). The W fragment is shown in a light blue color. Please note that a full structural characterization is not possible by RP-HRMS/MS, (B) The product ion spectrum in negative ion mode at m/z 1260.7237, showing characteristic fragments m/z 241 and 259, 355, 373 and 417. The sugar head group was confirmed by the [C3PO3] fragment (m/z 597, R1 = OH). [Z0PO3] and [Y1-H] fragments prove the ceramide moiety. (C) The positive ion mode ddMS2 spectrum of the [M + H]+ precursor, exhibiting the [W]+, [W-H2O]+ and [W-2H2O]+ fragments at m/z 298, 280 and 262, which are characteristic for the t18:1 LCB.
Figure 2. Overview of the GIPC fragmentation for the example of GIPC A-OH-t18:1/h24:0 in salad: (A) The fragment assignment of GIPC A-OH-t18:1/h24:0 (adapted from [23]). The W fragment is shown in a light blue color. Please note that a full structural characterization is not possible by RP-HRMS/MS, (B) The product ion spectrum in negative ion mode at m/z 1260.7237, showing characteristic fragments m/z 241 and 259, 355, 373 and 417. The sugar head group was confirmed by the [C3PO3] fragment (m/z 597, R1 = OH). [Z0PO3] and [Y1-H] fragments prove the ceramide moiety. (C) The positive ion mode ddMS2 spectrum of the [M + H]+ precursor, exhibiting the [W]+, [W-H2O]+ and [W-2H2O]+ fragments at m/z 298, 280 and 262, which are characteristic for the t18:1 LCB.
Metabolites 10 00375 g002aMetabolites 10 00375 g002bMetabolites 10 00375 g002c
Figure 3. The normalized ratio per gram dry weight for annotated GIPCs in salad (light-green), spinach (green), strawberries (rose) and raspberries (dark-red), by using different substituents for the functional group R1: (A) NAc, (B) NH2, and (C) OH (more detailed information can be found in Table A1).
Figure 3. The normalized ratio per gram dry weight for annotated GIPCs in salad (light-green), spinach (green), strawberries (rose) and raspberries (dark-red), by using different substituents for the functional group R1: (A) NAc, (B) NH2, and (C) OH (more detailed information can be found in Table A1).
Metabolites 10 00375 g003aMetabolites 10 00375 g003bMetabolites 10 00375 g003c

Share and Cite

MDPI and ACS Style

Panzenboeck, L.; Troppmair, N.; Schlachter, S.; Koellensperger, G.; Hartler, J.; Rampler, E. Chasing the Major Sphingolipids on Earth: Automated Annotation of Plant Glycosyl Inositol Phospho Ceramides by Glycolipidomics. Metabolites 2020, 10, 375. https://doi.org/10.3390/metabo10090375

AMA Style

Panzenboeck L, Troppmair N, Schlachter S, Koellensperger G, Hartler J, Rampler E. Chasing the Major Sphingolipids on Earth: Automated Annotation of Plant Glycosyl Inositol Phospho Ceramides by Glycolipidomics. Metabolites. 2020; 10(9):375. https://doi.org/10.3390/metabo10090375

Chicago/Turabian Style

Panzenboeck, Lisa, Nina Troppmair, Sara Schlachter, Gunda Koellensperger, Jürgen Hartler, and Evelyn Rampler. 2020. "Chasing the Major Sphingolipids on Earth: Automated Annotation of Plant Glycosyl Inositol Phospho Ceramides by Glycolipidomics" Metabolites 10, no. 9: 375. https://doi.org/10.3390/metabo10090375

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

Article Metrics

Back to TopTop