🫧 Lipidomic Class Profiling (Phospholipids/Sphingolipids)

Biphasic Lipid Extraction — Isolating the Membrane Lipidome from Biological Matrices

Every lipidomics experiment begins with quantitative, reproducible extraction of intact lipids from plasma, tissue, or cultured cells. Because glycerophospholipids and sphingolipids span an enormous range of polarity — from highly amphipathic phosphatidylinositols to nearly neutral ceramides — extraction chemistry must recover all classes with minimal bias while removing proteins, salts, and carbohydrates that would suppress electrospray ionization downstream. The Matyash MTBE/methanol/water partition has become the de facto standard, replacing the older chloroform-based Folch and Bligh-Dyer protocols with a safer, upper-phase organic layer that simplifies robotic automation.

• 92–97%: Organic-phase recovery (glycerophospholipids + sphingolipids)
• 13: SPLASH IS classes (deuterated/¹³C-labeled lipids)
• >99%: Protein removal (partitions to aqueous/pellet)
• 18 min: Extraction time (per 96-well plate batch)

Matyash MTBE partition and internal standard design

The extraction workflow proceeds through a defined sequence of physicochemical steps:

Step 1 — Protein precipitation and IS spiking: • 50 µL plasma (or tissue homogenate normalized to protein content by BCA assay) combined with ice-cold methanol (1:5, v/v) • SPLASH Lipidomix internal standard added before extraction: 13 stable-isotope-labeled lipids (d7-PC, d7-PE, d7-PS, d7-PI, d7-PG, d7-PA, d7-Cer, d7-SM, d7-DAG, d5-TAG, d5-Cholesteryl ester) at fixed concentrations spanning 0.1–10 µg/mL • Isotope labeling corrects for extraction efficiency, matrix effects, and ion suppression on a per-class basis

Step 2 — MTBE partition (Matyash et al. 2008): • MTBE added at 1:3 (v/v) methanol:MTBE, vortexed 10 min at room temperature • Water added (MTBE:MeOH:H2O ≈ 10:3:2.5) to induce phase separation • Upper MTBE-rich phase (less dense than water) contains >90% of glycerophospholipids and sphingolipids • Lower aqueous/methanol phase retains proteins, salts, and polar metabolites — cleanly separated for orthogonal metabolomics analysis from the same sample

Step 3 — Drying and reconstitution: • Upper phase collected, dried under gentle nitrogen stream (or SpeedVac, <30°C to prevent lipid oxidation) • Residue reconstituted in isopropanol:methanol:water (4:2:1, v/v/v) — a solvent matched to the downstream UHPLC mobile phase for optimal peak shape

Quality-control considerations: • Extraction blanks and pooled QC samples injected every 10 runs monitor carryover and drift • Antioxidant (butylated hydroxytoluene, 0.01%) added to extraction solvents suppresses lipid peroxidation artifacts, particularly for polyunsaturated species (e.g., PC 38:6, arachidonate- and DHA-containing) • Class-specific recovery bias is quantified by spike-recovery experiments: anionic lipids (PI, PA, cardiolipin) show slightly lower recovery (85–90%) than zwitterionic PC/PE (>95%) due to residual ion-pairing with divalent cations

Reverse-Phase Chromatography and High-Resolution Tandem Mass Spectrometry

Lipid class profiling can be performed by two complementary strategies: shotgun lipidomics (direct infusion, no chromatography, rapid but prone to isobaric interference) and liquid-chromatography-coupled MS (LC-MS/MS, slower but resolving lipid species by both mass and retention time). Reverse-phase UHPLC on charged-surface hybrid (CSH) C18 columns has become the workhorse approach because acyl-chain hydrophobicity elutes lipids in a predictable order within each class, adding a second dimension of specificity beyond m/z alone.

• CSH C18: Column (2.1×100mm, 1.7µm particle)
• 120,000: Orbitrap resolution (FWHM at m/z 200)
• 22 min: Gradient length (+ 3 min re-equilibration)
• ~4,200: Features detected (raw peaks before annotation)

Chromatographic separation and dual-polarity acquisition strategy

Reverse-phase separation exploits the amphipathic architecture of membrane lipids — a polar headgroup and one or two nonpolar acyl chains:

Retention behavior: • Within a class, retention time increases with acyl-chain carbon number and decreases with each additional double bond • PC 34:1 elutes before PC 36:1 (longer chain, more retained) but after PC 34:2 (more unsaturated, less retained) • This predictable elution order allows retention-time prediction models to flag misannotations and resolve isobaric/isomeric overlaps invisible to MS1 alone

Ionization polarity switching: • Positive ESI: efficiently ionizes PC, SM, and Cer as [M+H]+ or [M+HCOOH+H]+ adducts — the choline/amine headgroups readily protonate • Negative ESI: required for PE, PS, PI, PG, PA, and cardiolipin, which ionize efficiently as [M-H]- due to the phosphate headgroup • Polarity switching (fast MS scan alternation, <100 ms) within a single run doubles class coverage versus single-polarity acquisition, at the cost of ~30% duty-cycle scan time per polarity

Data-dependent (DDA) vs data-independent (DIA) MS2: • DDA: top-N most intense precursors selected per cycle for HCD fragmentation (stepped collision energy 20/30/40 eV) — provides clean, interpretable MS2 spectra but samples stochastically, missing low-abundance species run-to-run • DIA (e.g., SWATH, MSE): fragments all precursors within sequential m/z windows regardless of intensity — comprehensive but generates chimeric spectra requiring library-based deconvolution (e.g., Skyline, MS-DIAL DIA module) • Triple-quadrupole platforms (Sciex Lipidyzer, 6500+ QTRAP) instead use targeted multiple-reaction-monitoring (MRM) with class-specific neutral-loss/precursor-ion scans — e.g., PC detected by precursor ion scan of m/z 184 (phosphocholine headgroup fragment), Cer by neutral loss of the sphingoid base

Instrument comparison in routine clinical lipidomics: • Orbitrap HF-X: high resolution (120k) enables confident formula assignment and separates near-isobaric species (e.g., PC 34:1 vs PE 37:1, Δm = 0.036 Da) that unit-resolution instruments cannot distinguish • QTRAP/Lipidyzer: lower resolution but higher sensitivity and reproducibility for absolute quantification of ~1,100 pre-defined MRM transitions per 10-minute run — preferred for large clinical cohorts

From Raw Features to Lipid Species — Spectral Library Matching and FDR Control

A single UHPLC-MS/MS run generates thousands of raw chromatographic features — most representing adducts, in-source fragments, or noise rather than genuine lipid species. Converting this feature table into a confidently annotated lipidome requires isotope deconvolution, MS2 spectral library matching, and retention-time-based disambiguation, following the shorthand nomenclature standardized by the LIPID MAPS consortium and the 2020 lipidomics community update (Liebisch et al.).

• 850+: Species annotated (per plasma sample, both polarities)
• <5 ppm: Mass accuracy (MS1 precursor tolerance)
• <1%: Annotation FDR (decoy library controlled)
• 24: Lipid classes covered (LIPID MAPS categories GP + SP)

Annotation pipeline: deconvolution, spectral matching, and isobaric resolution

Software such as MS-DIAL, LipidSearch, Lipostar, or LipidXplorer converts raw .raw/.wiff files into an annotated lipid table through a standardized pipeline:

1. Peak detection and alignment: • Continuous wavelet transform or local-maximum peak picking across all samples • Retention-time alignment (dynamic time warping) corrects for chromatographic drift across a batch of 100+ injections • Isotope-pattern deconvolution groups ¹³C isotopologues and charge states into a single feature

2. MS2 spectral matching: • Fragment spectra matched against curated libraries: LIPID MAPS Structure Database (LMSD, >48,000 structures), LipidBlast (in silico predicted spectra), or instrument-specific empirical libraries • Class-diagnostic fragments confirm headgroup identity: m/z 184.07 (phosphocholine, PC/SM), m/z 241.01 (glycerophosphoinositol, PI), neutral loss of 87 Da (serine, PS) • Acyl-chain-diagnostic fragments (fatty-acyl carboxylate anions in negative mode) assign chain composition: PC 34:1 confirmed as PC 16:0_18:1 by observing m/z 255.2 (FA 16:0) and m/z 281.2 (FA 18:1) product ions

3. Shorthand nomenclature assignment (LIPID MAPS / Liebisch 2020 update): • Class abbreviation + total carbons:double bonds, e.g. "PC 34:1" (sum composition) or "PC 16:0_18:1" (molecular species) when acyl regiochemistry is resolved • Sphingolipids use the "d/t" prefix denoting dihydroxy/trihydroxy sphingoid base: "Cer d18:1/24:0" = C18 sphingosine backbone amide-linked to lignoceric acid

4. Isobaric and isomeric disambiguation: • Ether-linked plasmalogen PE (PE P-36:4) is nominally isobaric with diacyl PE (PE 36:5, Δm=0.0106 Da) — resolved by 120k-resolution MS1 plus distinct retention time and characteristic neutral-loss fragmentation • sn-1/sn-2 acyl positional isomers (e.g., PC 16:0/18:1 vs PC 18:1/16:0) require specialized ozone-induced dissociation (OzID) or UVPD — not resolved by standard HCD, and most clinical pipelines report only sum or combined molecular-species composition

5. False discovery rate control: • Decoy spectral library (mass-shifted or scrambled fragment lists) estimates the false-match rate at a given matching score threshold • Target FDR <1% typically requires dot-product spectral match score >0.7 and mass accuracy <5 ppm simultaneously

Absolute Quantification and Multivariate Modeling of the Lipidome

Raw peak areas are not directly comparable across lipid classes or batches — ionization efficiency varies by headgroup, acyl chain length, and unsaturation by more than 10-fold. Class-matched internal standard normalization converts peak area ratios into absolute concentrations (pmol per mg protein or per mL plasma), enabling multivariate statistical modeling to identify which of the ~850 annotated species differ significantly between clinical or experimental groups.

• n = 120: Cohort size (typical case-control study)
• 62: VIP > 1 species (PLS-DA variable importance)
• 38: Significant hits (BH-FDR q < 0.05)
• R²Y = 0.89: PLS-DA model fit (Q² = 0.71 cross-validated)

Normalization, scaling, and statistical testing workflow

Quantitative lipidomics statistics follow a defined chain of transformations before biological interpretation:

1. Internal standard normalization: • Each annotated species is quantified relative to its class-matched deuterated internal standard: concentration = (analyte peak area / IS peak area) × IS spiked concentration × dilution factor • Class-matching is essential — normalizing a PI species to a PC-class standard introduces systematic bias because ionization efficiency differs by headgroup chemistry (typically 2–8-fold between classes) • Where a 1:1 isotope-labeled standard is unavailable for every species, the nearest class-representative standard is used with a documented correction factor

2. Batch correction and QC filtering: • Pooled QC samples (injected every 10th run) monitor signal drift; species with QC coefficient of variation >30% are flagged or excluded • LOESS/QC-RLSC batch correction removes systematic instrument drift across long acquisition sequences (>200 injections)

3. Data transformation and scaling: • Log2 transformation stabilizes variance across the wide dynamic range of lipid concentrations (sub-nM to low-µM) • Pareto scaling (divide by square root of standard deviation) is preferred over unit-variance scaling in lipidomics because it retains proportionally more weight for high-abundance structural lipids while dampening noise inflation

4. Univariate testing: • Moderated t-tests (limma empirical Bayes) borrow variance information across all ~850 species, improving power in modest cohort sizes • Benjamini-Hochberg FDR correction controls the expected proportion of false positives among the significant calls: q<0.05 across 850 tests still permits robust discovery power

5. Multivariate modeling: • PLS-DA (partial least squares discriminant analysis) and OPLS-DA model class separation using latent variables that maximize covariance between the lipid profile and group label • Variable Importance in Projection (VIP) scores >1 flag species driving group separation; permutation testing (1,000 iterations) and 7-fold cross-validation guard against overfitting on high-dimensional, modest-n data • Random forest and elastic-net logistic regression are increasingly used as orthogonal, less overfitting-prone alternatives, often combined into consensus biomarker panels

From Lipid Class Shifts to Membrane Signaling and Disease Biomarkers

The ultimate purpose of lipidomic class profiling is mechanistic and clinical translation: connecting quantitative shifts in phospholipid and sphingolipid classes to membrane biophysics, signal-transduction pathways, and validated disease biomarkers. Ceramide and sphingomyelin ratios, phosphatidylcholine-to-phosphatidylethanolamine balance, and lysophospholipid pools are among the most reproducibly disease-associated lipid signatures across metabolic, neurodegenerative, and cardiovascular disease cohorts.

• 2.3×: Cer d18:1/16:0 fold-change (NAFLD fibrosis vs. healthy)
• ↓ 18%: PC/PE ratio (associated with membrane fragility)
• 0.86: Biomarker panel AUC (ceramide-based CV risk score)
• 3: Independent validation cohorts (total n > 900 patients)

Membrane biophysics and signaling roles of major lipid classes

Each lipid class profiled by mass spectrometry maps onto specific structural or signaling functions in the cell membrane:

• Phosphatidylcholine (PC) and phosphatidylethanolamine (PE): the two most abundant glycerophospholipid classes, together comprising 40–60% of total membrane phospholipid. PC (cylindrical, forms stable bilayers) and PE (conical, promotes negative membrane curvature) must be maintained in a defined ratio (typically PC:PE ≈ 2:1 in most mammalian membranes) — deviations are associated with hepatocyte membrane fragility and are a hallmark of nonalcoholic fatty liver disease (NAFLD) progression to steatohepatitis

• Phosphatidylserine (PS): normally confined to the inner leaflet by flippase activity; externalization to the outer leaflet is the canonical "eat-me" signal recognized by macrophages during apoptosis, making plasma PS species markers of cell turnover and clearance efficiency

• Phosphatidylinositol (PI) and phosphoinositides: precursors for PI3K/AKT and PLC/PKC signaling cascades; altered PI class abundance reflects broad perturbation of growth-factor and insulin signaling

• Cardiolipin (a tetra-acyl anionic phospholipid unique to the inner mitochondrial membrane): its acyl-chain remodeling (via the enzyme tafazzin) is essential for oxidative phosphorylation complex stability — cardiolipin abnormalities are the molecular basis of Barth syndrome and are increasingly profiled in mitochondrial disease panels

• Sphingomyelin (SM) and ceramide (Cer): SM is hydrolyzed by sphingomyelinase to generate ceramide, a potent pro-apoptotic second messenger that promotes mitochondrial outer-membrane permeabilization; elevated ceramide (particularly Cer d18:1/16:0) impairs insulin signaling in skeletal muscle and is one of the most extensively validated circulating biomarkers of cardiometabolic risk

• Sphingosine-1-phosphate (S1P): a bioactive lysophospholipid signaling through S1P receptors (targeted therapeutically by fingolimod in multiple sclerosis) that regulates lymphocyte trafficking and vascular integrity; the ceramide/S1P ratio ("sphingolipid rheostat") determines cell survival versus apoptosis fate

• Lysophosphatidylcholine (LPC): generated by phospholipase A2 hydrolysis of PC; depleted LPC pools accompany chronic inflammatory states and correlate with reduced membrane repair capacity

The plasma "Ceramide Score" — a validated four-species ratio panel (Cer d18:1/16:0, d18:1/18:0, d18:1/24:1 relative to d18:1/24:0) — is now offered as a clinical laboratory test (CERT1/CERT2, Mayo Clinic) that outperforms LDL cholesterol alone for predicting major adverse cardiovascular events, with a validated AUC of 0.86 across three independent cohorts totaling more than 900 patients. This represents one of the clearest translations of quantitative class-and-species-level lipidomic profiling into routine clinical decision-making.