If you don't remember your password, you can reset it by entering your email address and clicking the Reset Password button. You will then receive an email that contains a secure link for resetting your password
If the address matches a valid account an email will be sent to __email__ with instructions for resetting your password
1 Present address of S. E. Hancock: School of Medical Sciences, University of New South Wales, Sydney, Australia.
Sarah E. Hancock
Footnotes
1 Present address of S. E. Hancock: School of Medical Sciences, University of New South Wales, Sydney, Australia.
Affiliations
School of Medicine and Molecular Horizons, University of Wollongong, Wollongong, AustraliaIllawarra Health and Medical Research Institute, Wollongong, Australia
School of Medicine and Molecular Horizons, University of Wollongong, Wollongong, AustraliaIllawarra Health and Medical Research Institute, Wollongong, Australia
Secretions from meibomian glands located within the eyelid (commonly known as meibum) are rich in nonpolar lipid classes incorporating very-long (22–30 carbons) and ultra-long (>30 carbons) acyl chains. The complex nature of the meibum lipidome and its preponderance of neutral, nonpolar lipid classes presents an analytical challenge, with typically poor chromatographic resolution, even between different lipid classes. To address this challenge, we have deployed differential mobility spectrometry (DMS)-MS to interrogate the human meibum lipidome and demonstrate near-baseline resolution of the two major nonpolar classes contained therein, namely wax esters and cholesteryl esters. Within these two lipid classes, we describe ion mobility behavior that is associated with the length of their acyl chains and location of unsaturation. This capability was exploited to profile the molecular speciation within each class and thus extend meibum lipidome coverage. Intriguingly, structure-mobility relationships in these nonpolar lipids show similar trends and inflections to those previously reported for other physicochemical properties of lipids (e.g., melting point and phase-transition temperatures). Taken together, these data demonstrate that differential ion mobility provides a powerful orthoganol separation technology for the analysis of neutral lipids in complex matrices, such as meibum, and may further provide a means to predict physicochemical properties of lipids that could assist in inferring their biological function(s).
Meibomian glands, located within the eyelids, produce a waxy, lipid-rich substance known as meibum. Meibum is thought to supply the outermost tear film lipid layer, which helps prevent the evaporation of the aqueous layer beneath it. The composition of the tear film lipid layer influences its stability (
The International Workshop on Meibomian Gland Dysfunction: Report of the Subcommittee on Tear Film Lipids and Lipid–Protein Interactions in Health and Disease.
). The meibum lipidome is quite complex and is unusual in that it contains mostly neutral, nonpolar lipid classes that contain very-long (22–30 carbons) to ultra-long (>30 carbons) acyl chains. The two most abundant of these, the wax esters (WEs) and cholesteryl esters (CEs), make up around 75 mol% of total meibum lipids (
Two MS strategies have previously been used to identify and quantify meibomian lipids: direct infusion with tandem MS (also known as shotgun lipidomics) and LC/MS. Both approaches use soft ionization techniques, such as atmospheric pressure chemical ionization or ESI. Shotgun lipidomics has been successful in identifying and quantifying a range of lipid species within meibum (
), but the presence of isobaric species (i.e., lipids with the same nominal mass but differing chemical composition) can complicate spectra and impede the discovery of new compounds (e.g., CEs of the form CE x:2 are isobaric with WEs of the form [x+ 26]:0). LC/MS-based workflows have also been successful in mapping much of the meibum lipidome (
). Nevertheless, LC cannot completely resolve the neutral meibum lipid classes. Meibum separation by normal-phase LC typically displays a single broad peak under which all of the nonpolar lipid classes elute (
Ion mobility spectrometry (IMS) coupled to MS has recently been applied to the separation of a range of isomeric and isobaric compounds in the gas phase (
Quantitative structural multi-class lipidomics using differential mobility: electron impact excitation of ions from organics (EIEIO) mass spectrometry.
Comprehensive lipidomic analysis of human plasma using multidimensional liquid- and gas-phase separations: two-dimensional liquid chromatography–mass spectrometry vs. liquid chromatography–trapped-ion-mobility–mass spectrometry.
) in both standard mixtures and biological extracts. Conventional ion mobility such as drift-time ion mobility spectrometry (DTIMS) measures the migration time of ions through an inert buffer gas under low electric field conditions, with the mobility of a given ion being inversely related to its collisional cross-section (CCS), a function of the mass-to-charge ratio (m/z) of the ion, its shape and size, and its interactions with drift gas (
). Traveling-wave ion mobility spectrometry (TWIMS) is a variation on DTIMS in which ions are subjected to an electric field that is applied to one section and set to sequentially sweep through the length of the ion mobility cell in the direction of ion migration. In both approaches, libraries of drift time and CCS are generated from standards to aid in the identification of isomeric or isobaric lipids. Separation of isomeric and/or isobaric lipids by DTIMS or TWIMS is limited by the resolving power (i.e., a function of drift-tube distance and pressure), with even high-resolution DTIMS being unable to separate several isomeric phospholipids (
). More recently, an ultra-high resolution TWIMS platform known as structures for lossless ion manipulations achieved baseline separation of phosphatidylcholine double bond positional isomers and near-baseline separation of stereoisomers, a feat that was enabled by two passes through the serpentine ion path of the structures for lossless ion manipulations platform, equivalent to a drift-tube distance of 30.6 m (
In contrast to ion mobility measured with DTIMS or TWIMS, differential mobility spectrometry (DMS), also known as field asymmetric ion mobility, separates ions in a given population by exploiting their differential mobility under low- and high-field conditions when exposed to an alternating asymmetric waveform known as a separation voltage (SV). Under the influence of this SV, ions will oscillate and drift toward one of the electrodes unless their trajectory is corrected by a superimposed DC compensation voltage (CV). By scanning through CV one can derive an ionogram based on the mass spectrometric detection of the ion populations being transmitted through the DMS cell for a given SV. The separations observed in DMS are often found to have a lower m/z dependence and thus greater structural sensitivity than either DTIMS or TWIMS. DMS has been used as a filter to separate lipid classes and/or isomers/isobars in both synthetic mixes and biological extracts prior to structural determination by tandem MS (
Quantitative structural multi-class lipidomics using differential mobility: electron impact excitation of ions from organics (EIEIO) mass spectrometry.
). These combination technologies have been applied to separate lipids based upon within-class structural characteristics such as chain length and unsaturation (
Herein, we show the ability of DMS-MS to separate the two main lipid classes present in meibum, WEs and CEs, with near-baseline resolution, a feat that is not achievable using common chromatographic strategies. The DMS-MS workflow applied here reveals several characteristics driving this separation by DMS-MS, including the influence of carbon chain length and double bond position on mobility behavior. By combining DMS with OzID, 123 unique WE and CE species were characterized within meibum, including four double bond isomers within each of the seven CE sum compositions. In addition, the relationship between ion mobility behavior and carbon chain length in both WEs and CEs and the double bond position in WEs follows a similar trend to that observed for their phase-transition temperature (TM). Therefore, DMS may not only provide an orthogonal approach to the separation of nonpolar lipids that can be readily incorporated into shotgun or LC/MS lipidomic workflows but may also be a way of predicting the physiochemical properties of newly discovered lipid molecular species within complex lipidomes.
MATERIALS AND METHODS
Materials
LC/MS-grade methanol and HPLC-grade chloroform were obtained from VWR International (Brisbane, Australia). Ammonium acetate and lithium acetate of the highest purity available were purchased from Sigma-Aldrich (Sydney, Australia). All experiments were conducted in compliance with the tenets of the Declaration of Helsinki and were approved by the University of New South Wales Human Research Ethics Advisory Panel. The meibum sample was a pooled sample from 170 human donors collected and extracted by approved protocols as described previously (
). Briefly, meibum samples were collected from patients by meibomian gland evaluator, with no topical anesthesia used. The meibomian gland evaluator was placed below the eyelash line of the lower eyelid and was held in this position for 10 s. A stainless-steel ethanol-cleaned and heat-sterilized spatula was used to collect the expressed meibum by gently sliding across the orifices. Meibum collected on the spatula was dissolved in a glass vial containing ∼1 ml chloroform. This sample was then transferred to a sleeved glass vial, dried under N2, and resuspended in 200 μl chloroform. Samples were pooled prior to analysis.
Nomenclature
Lipid nomenclature used throughout is based on general literature recommendations (
). WEs are denoted as X1:Y1/X2:Y2, where X1:Y1 define the number of carbons and double bonds present on fatty acid and X2:Y2 the fatty alcohol. CEs are described as X:Y, where X is the total number of carbons and Y is the number of double bonds in the fatty acyl chain. The double bond position is described by either n-x or Δy nomenclature, where x and y are the number of carbon-carbon bonds from the methyl end of the fatty acid (i.e., omega position) or carboxylate moiety, respectively.
DMS-MS
DMS-MS was conducted using a QTRAP 5500 equipped with a SelexION DMS interface and TurboV ESI source controlled by Analyst 1.5.2 (SCIEX, Concord, Canada). ESI voltage was set at +5.5 kV, with the source equilibrated to room temperature. Nitrogen was used for nebulizing (15 psi), curtain (20 psi), and collision (9 arb units) gas. The pooled meibum extract was diluted 100-fold in methanol spiked with 7.5 mM ammonium acetate and directly infused into the mass spectrometer at a rate of 5 µl/min. In a typical DMS experiment SV (+4500 V), resolving gas (nitrogen; 25 psi), and DMS cell temperature (225°C) were held at an optimized value while CV was scanned and MS data were acquired. No gas-phase modifiers were used. DMS-MS data were acquired from m/z 50–1200 with a declustering and exit potential of +100 V and +10 V, respectively. A wide CV range (−20 to +30 V) and large step size (0.5 V) were used initially to determine the optimal transmission range for meibum lipids ionized as [M+NH4]+ cations, after which the CV range and step size were narrowed (+6 to +22; step size = 0.1 V). Precursor ion scans for CEs and WEs were obtained using settings reported previously (
). An ozone generator (Titan 30; Absolute Ozone, Alberta, Canada) was used to produce ozone (∼14.7% w/w O2) from oxygen supplied at a constant flow rate (0.2 l/min), pressure (20 psi), and power setting 30% that was then delivered to the collision cell of the mass spectrometer via a variable leak valve (Nenion, Lustenau, Austria) through the collision gas line. For OzID experiments, meibum extract was diluted 50-fold in methanol spiked with 2 mM lithium acetate and infused into the mass spectrometer at a rate of 5 µl/min. DMS settings were as described above, while experimental details for DMS-MRM of [M+Li]+ monounsaturated CEs are described in supplemental Table S5. For DMS-OzID a modified time-delayed fragmentation scan was used to mass-select monounsaturated CE [M+Li]+ ions in Q1 that were then transferred into Q2 with minimal collision energy (e.g., 5 eV). Target ions were trapped in Q2 with ozone and allowed to react for 1 s before being ejected into Q3. OzID fragment and residual precursor ions were then scanned out for mass analysis at 1,000 Th/s. The observed neutral losses in the resulting OzID mass spectra were used to assign the site(s) of unsaturation in the corresponding lipids using the previously predicted transitions listed in Table 1 (
). DMS-OzID of lithiated CE species were acquired using the SV, temperature, and pressure settings described above while ramping the CV from +6 to +12 V (step size = 0.01 V; n = 4)
Fig. 1A diagram of the QTRAP5500 mass spectrometer illustrating the location and configuration of the DMS device used in these experiments and highlighting the hardware modifications made to the instrument to facilitate OzID.
R Core Team. 2013. R: A language and environment for statistical computing. R Foundation for Statistical Computing. Available from https://www.r-project.org.
). Gaussian filtering was applied with a window size determined by a combination of the normalized SD and number of data points across the full width at half maximum (i.e., CV step size) and were visualized using the ggplot2 package (
). Data presented were averaged across four to five measurements from the same sample on differing days and are presented as mean ± RSD. Piecewise regression between carbon chain length and CV was conducted in R using the segmented package with statistical significance set at P < 0.05 (
). Given the challenge of separating the neutral lipid-dominated lipidome of meibum by LC, we sought to determine if DMS could achieve greater separation of the major two meibum lipid classes, CEs and WEs. To that end, we subjected a methanolic solution of meibum spiked with ammonium acetate to ESI to generate abundant [M+NH4]+ ions from WE and CE species. In initial tests, DMS-MS data were acquired while ramping the CV from +6 to +22 V (step size = 0.1 V). The resulting plot of CV against total ion current detected by the mass spectrometer (hereafter referred to as an ionogram) is shown in Fig. 2A. This ionogram displayed a bimodal distribution, suggesting the separation of two distinct ion populations within the gas phase by DMS. We probed the identity of each ion population by holding the CV at a given voltage and acquiring the resulting mass spectrum (Fig. 2B, D). The mass spectrum obtained at CV +11.2 V (Fig. 2B) displayed a number of ions with mass-to-charge ratios (m/z) that closely matched the pattern of [M+NH4]+ ions obtained from a precursor ion scan targeting CEs (m/z 369; Fig. 2C). Likewise, the masses observed at CV +16.1 V (Fig. 2D) were similar to those of [M+NH4]+ WE obtained by a precursor ion scan for a protonated 18:1 fatty acid (m/z 283; Fig. 2E). Taken together, these findings indicate that DMS-MS can separate the CE and WE lipid populations present within meibum as gas-phase ions.
Fig. 2DMS separation of human meibomian gland secretions results in the bimodal ionogram (A). When the CV is held at +11.2 V, the resulting mass spectrum (B) closely resembles a precursor ion scan for CEs (m/z 369; C). Similarly, holding the CV at +16:1 V produces a mass spectrum (D) that resembles a precursor ion scan for WE species containing an 18:1 fatty acid (m/z 283; E).
To confirm the separation of CE and WE classes in meibum by DMS-MS, we ran a series of DMS-MRM experiments targeting known meibum CE and WE species as [M+NH4]+ ions using transitions described by previous studies (supplemental Tables S1–S3), and WE species present within meibum at >0.5% abundance of total WEs are shown in Fig. 4. Meibum WEs generally had maximum transmission through the DMS cell at CVs that were inversely proportional to their carbon chain length as might be expected (
), with each carbon increase resulting in a ∼0.3 V decrease in CV (Fig. 4). Unsaturated WEs with a double bond in either the fatty acid or fatty alcohol domain were typically transmitted through the DMS cell at a lower CV compared with their saturated counterparts. For example, saturated WE 16:0/24:0 had a maximum CV of +17.7 ± 0.1 V, while monounsaturated WEs 16:1/24:0 and 16:1/24:1 had maximal CVs of +17.2 ± 0.5 V and +15.7 ± 0.4 V, respectively. Interestingly, for unsaturated WEs, the location of the double bond within the fatty acid or fatty alcohol also influenced its mobility. Unsaturated WEs characteristically displayed lower maximum CVs than their saturated equivalents, with a double bond in the fatty alcohol portion of the WEs decreasing maximum CV by ∼1.3 V on average compared with a reduction of only 0.6 V when the double bond was located in the fatty acid (Fig. 3).
Fig. 4Ionograms of WEs detected by DMS-MRM present at an abundance of >0.5% of total WEs in meibum. WEs are shown grouped by fatty acyl chain (boxes) and fatty alcohols (x axis). WE species are presented as mean maximal CV transmission (closed circles) ± RSD (n = 4). The shaded areas indicate the CV range and relative intensity at which the WE species migrated through the DMS cell.
Fig. 3DMS-MRM of known CE and WE species present within meibum shows near-baseline separation. DMS-MRM data were acquired using a separation voltage of +4500 V, a resolving gas pressure of 25 psi, and a DMS cell temperature of 225°C. CV was scanned between +6 and +22 V with a step size of 0.1 V. Settings optimized for DMS-MRM experiments of both WEs and CEs are available in supplemental Tables S1–S4.
For CEs, close to 60 discrete [M+NH4]+ species present in meibum were surveyed by DMS-MRM (supplemental Table S4), with the ionograms for CEs present at >0.5% abundance of total CEs shown in Fig. 5. The ion mobility behavior seen with CEs in many ways paralleled that observed for WEs (see above), where transmission by DMS was strongly associated with the mass of the molecule and its degree of unsaturation. Like WEs, the CV for the maximum transmission of a given saturated CE species was inversely related to its mass (Fig. 5A). Saturated CEs migrated through the DMS cell at CVs ranging from +8 V to +16 V, with a noticeable plateau in this relationship for very-long chain CEs (i.e., 22–30 carbons). For long-chain saturated CEs 16 to 21 carbons in length, the maximum CV decreased on average by 0.5 V for each single carbon increase in chain length. In contrast, for the very-long chain saturated CEs the relationship between CV and chain length attenuates to less than 0.1 V per carbon number increase. Piecewise regression confirmed this breakpoint seen in the relationship between saturated CE mass and maximum CV after 21 carbons (R2 = 0.992).
Fig. 5Ionograms of (A) saturated, (B) monounsaturated, and (C) polyunsaturated meibum CE species (CE 20:2 and CE 24:2) present in meibum at an abundance of >0.5% of total CEs. CE species are shown as mean maximal CV transmission (closed circles) ± RSD (n = 4). The shaded areas indicate the CV range and relative intensity at which the CE species migrated through the DMS cell. Solid lines show fitted line(s) from piecewise regression.
Monounsaturated CEs (i.e., CEs with one double bond in their acyl chain) were transmitted through the DMS cell at a slightly lower CV range than their saturated counterparts (+7.5 to +13.7 V; Fig. 5B). However, an unusual U-shaped relationship between carbon chain length and mean maximum CV was observed for monounsaturated CEs, with piecewise regression identifying two breakpoints: the first between monounsaturated CEs 22 and 24 carbons in length and the second for CEs with acyl chains 28 carbons and longer (R2 = 0.990). For long-chain CEs with monounsaturated fatty acyl chains between 18 and 22 carbons in length, the average maximum CV decreased by 0.9 V for each two-carbon addition. This trend reversed entirely for monounsaturated CEs with 24 to 28 carbons in their acyl chain, with the mean maximal CV increasing by 0.7 V for every two-carbon elongation. For very-long to ultra-long chain monounsaturated CEs with fatty acids longer than 28 carbons there was little observable change in maximal CV with each two-carbon increase, with CEs 28:1, 30:1, and 32:1 averaging a similar maximum CV as CE 18:1. As monounsaturated CEs are known to exist as a number of double bond positional isomers in meibum (
), it follows the that the site of unsaturation in the fatty acid of the CE could be driving the discontinuous relationship seen between mass and maximum CV for monounsaturated CEs.
Influence of double bond position on CE separation by DMS
To test our hypothesis that the position of the double bond on the acyl chain of monounsaturated CEs was affecting its ion mobility, we characterized monounsaturated CE double bond isomers by OzID (
). During OzID, mass-selected CE ions were trapped within the collision cell of the mass spectrometer in the presence of ozone, which resulted in characteristic pairs of product ions indicative of double bond position (i.e., Criegee and aldehyde ions spaced by 16 Da; Table 1). The [M+NH4]+ ions formed from monounsaturated CEs produced few neutral loss fragments indicative of double bond position, even after up to 5 s trapping with ozone (supplemental Fig. S1A). However, cationization of monounsaturated CEs with lithium resulted in abundant diagnostic Criegee and aldehyde ions within 1 s of trapping time with ozone (supplemental Fig. S1B). Accordingly, DMS-MRM transitions of [M+Li]+ monounsaturated CEs were acquired (n = 5; supplemental Table S5) to confirm similar ion mobility behavior, with the resulting ionogram displaying an analogous U-shaped relationship between mass and CV to that observed for [M+NH4]+ monounsaturated CEs acquired under the same DMS-MS conditions (supplemental Fig. S2). Some small differences in maximal CV transmission for specific CE species were seen between the two ion adducts, including a lower CV range for [M+Li]+ adducts (+6.7 to +11.9 V) and a higher maximum CV for CE 30:1 versus CE 18:1 (Fig. 5B).
Subsequently, DMS-OzID spectra were acquired on mass-selected monounsaturated CE [M+Li]+ ions while ramping CV from +6 to +12 V (n = 4). Monitoring the extracted ion chromatograms for diagnostic ions indicative of double bond position revealed the presence of four different double bond isomers within each monounsaturated CE, including n-5, n-7, n-9, and n-11 (see supplemental Fig. S3 for an example using CE 20:1). These extracted ion chromatograms for each CE double bond isomer were then used to produce an ionogram that compares maximal CV between monounsaturated CE species with the same double bond position (supplemental Fig. S4). Surprisingly, monounsaturated CEs with the same n double bond position had distinct mean maximal CVs, suggesting that the double bond position from the methyl end was not a major factor driving the separation of monounsaturated CEs by DMS (supplemental Fig. S4). The double bond position in lipids can also be referenced from the carboxylate moiety of the acyl chain (i.e., Δy). When Δ in double bond position was plotted against the mean maximal CV of isomeric monounsaturated CE species (Fig. 6), lipids with a common position of the double bond with respect to the carboxylate moiety were found to group together, suggesting this to be a better predictor of CE differential ion mobility behavior.
Fig. 6Ionograms extracted from DMS-OzID spectra of monounsaturated CE meibum species grouped by their Δ double bond position (i.e., number of carbons present between the double bond and the carboxylate moiety). CE species are shown as mean maximal CV transmission (closed circles) ± RSD (n = 4). The shaded areas indicate the CV range and relative intensity at which the CE species migrated through the DMS cell.
Herein we have shown the application of DMS-MS to the separation of lipids present within meibum (i.e., secretions from human meibomian glands), a substance that is composed of a principally nonpolar lipidome. Using this approach, we detected 123 unique meibum WE and CE species in total (Table 2). Analyzing lipidomes that contain predominately nonpolar lipid species typically presents an analytical challenge, as the lipid species contained therein are inadequately separated by commonly used LC/MS methods. While DMS-MS separation achieved near-baseline resolution of the two major lipid classes present in meibum, WEs and CEs, individual lipids within these two classes were also dispersed based on carbon chain length and the position of carbon-carbon double bonds.
TABLE 2Summary of all WE and CE species detected in meibum extracts by DMS-MS
The relationship between the size of a lipid (as a function of the number of carbons in its acyl chain) and its time of arrival during low-field ion mobility is well documented (
Comprehensive lipidomic analysis of human plasma using multidimensional liquid- and gas-phase separations: two-dimensional liquid chromatography–mass spectrometry vs. liquid chromatography–trapped-ion-mobility–mass spectrometry.
). In conventional IMS such as DTIMS and TWIMS, larger or more extended molecules have a greater CCS and therefore interact more with the drift gas, resulting in longer arrival-time distributions. When DMS is conducted without any gas-phase modifiers, lipids typically exhibit type-C differential ion mobility behavior, where the CV at which gas-phase ions successfully traverse the DMS cell rises with increasing SV (
). Type-C behavior is dominated by hard-sphere interactions, with minimal ion clustering during the low-field portion of the waveform. Under these conditions, physical characteristics of the molecule such as CCS largely drive the mobility of the ion (
). Consequently, a general trend between ion mobility behavior and mass has been reported for DMS of lipids with no gas-phase modifiers, where larger lipids typically have lower CVs (
). Thus far, linear relationships between CV and mass for several lipid classes have been reported in the literature with the exception of TGs, which show a slight attenuation of the relationship between m/z and CV with increasing size (
). In the present study, we observed an inverse but primarily linear relationship between WE mass and CV (Fig. 4), whereas both saturated and unsaturated CEs had a discontinuous relationship between size and CV (Fig. 5). For saturated CEs, we observed a distinct plateau in the relationship between CV and size for CEs with very-long acyl chains (Fig. 5A). This change in differential ion mobility behavior cannot be solely attributed to the mass of the lipid given that meibum CEs are, on average, up to about two-thirds the mass of the largest TG. However, the fatty acyls found in meibum CEs are substantially longer than those usually found in other animal tissues, which typically have a maximum of 24 carbons in any given fatty acid. Interestingly, the relationship between saturated CE carbon chain length and CV bears a resemblance to the relationship between acyl chain length and melting point for long-chain and ultra-long-chain free fatty acids (
) and shows reasonable agreeance when the two values are correlated (supplemental Fig. S5). This suggests that the same physiochemical properties that determine the temperature at which saturated lipids transition from the solid to liquid phase may also affect their ion mobility behavior under an asymmetric waveform. Of note, increased rotamer disorder within the lipid hydrocarbons as the acyl chain lengthens would lead to a larger cross-sectional area, thus influencing its ion mobility. We hypothesize that at a certain hydrocarbon length the difference in CCS between the ordered and disordered state would become comparatively smaller, resulting in the attenuated differential ion mobility behavior seen for very-long chain CEs (Fig. 5A).
Unsaturated lipids also display a reduction in ion mobility transit time during both DTIMS (
), which is a product of the smaller CCS imparted by the double bond. Similar behavior has also been reported for a number of other lipid classes during DMS, in which unsaturated lipids have lower CVs compared with their saturated counterparts, including phosphatidylcholines (but not lysophosphatidylcholines), phosphatidylethanolamines, TGs, and diacylglycerols (
). In the present study, we observed that the location of the carbon-carbon double bond within a WE had a significant impact on differential ion mobility, with a double bond being present in the fatty alcohol of the WE, resulting in a greater reduction in CV than a double bond in the fatty acid (Fig. 4). This behavior might be partially explained by the extreme asymmetrical structure of meibum WEs, where the carbon chain length of fatty alcohol can be anywhere from 2 to 16 carbons longer than the fatty acid. Indeed, variability in ion mobility behavior for asymmetric versus symmetric phospholipids with saturated acyl chains has been well documented with DMS, with greater differences in CV seen when the discrepancy between the two acyl chains is larger (
). Greater differences in CV can be seen between isomers when the longer acyl chain contains a double bond, allowing the resolution of sn-positional isomers as [M+Ag]+ ions by DMS (
). This long carbon chain might impart greater flexibility to the molecule, thus affecting the CCS of the ion and its ability to interact with any background gases as the ion moves through the high- and low-field conditions of the DMS cell. A similar relationship has also been described between transition temperature (TM) and double bond position in WEs, where a WE with a double bond present in the fatty alcohol has a lower TM than a WE with a double bond in the fatty acid (
). This further supports our hypothesis that the same physicochemical properties that influence lipid physical properties such as TM also direct their mobility as type-C ions under an asymmetric waveform.
The most striking ion mobility behavior was observed for the DMS-MS of monounsaturated CE species. Monounsaturated CE displayed a U-shaped relationship between their carbon chain length and the CV at which they migrated through the DMS cell (Fig. 5B). Moreover, this U-shaped relationship was best correlated with the position of the double bond from the carboxylate moiety (Fig. 6). Atomistic modeling has shown that liquid-phase cholesteryl oleate exists on a spectrum of folded to unfolded ensembles, with the greatest flexibility being present between carbons 3 and 5 on the acyl chain (
). In previous work, we determined that an interaction exists between the site of charge localization on the carbonyl and both the endocyclic and acyl carbon-carbon double bonds in cholesteryl oleate that potentially aids folding of the structure in the gas phase, resulting in a smaller CCS and enhanced ion mobility (
). Therefore, it stands to reason that the distance between the sterol head group and double bond may also affect the folding of the steryl ester in the gas phase, with elongated distances between the sterol and double bond resulting in compact structures when the double bond is located >15 carbons from the carboxylate moiety. To our knowledge no measurement of melting point or TM for lipids with ultra-long monounsaturated fatty acids like those typically found in meibum have been reported in the literature, preventing more definitive conclusions from being drawn on the relationship between differential ion mobilty behavior and this property of ultra-long-chain monounsaturated CEs.
CONCLUSIONS
We have described the ability of DMS to achieve near-baseline separation of the two dominant neutral lipid classes present in human meibum, WEs and CEs. This strategy is shown to be effective in removing isobaric overlaps in direct infusion analysis of nonpolar lipidomes. Combined with collision and ozone-induced dissociation DMS-MS/MS revealed 123 individual molecular contributors to WE and CE classes, including four double bond isomers within seven monounsauturated CEs. DMS separation of discrete molecular species within these two neutral classes revealed the influence of carbon chain length and double bond position on differential ion mobility behavior, where an inverse relationship between lipid acyl chain length and CV were observed. We also observed a distinct effect of unsaturation position on ion mobility behavior, where the inclusion of a double bond in either the fatty acid or alcohol in WEs influenced the CV at which it migrated through the DMS cell. Unsaturated CEs showed the most striking DMS behavior, with the position of the double bond relative to the carboxylate moiety best correlating with ion mobility. Many of the observed behaviors under the asymmetric ion mobility field bore similar relationships between lipid physical properties and their TM, indicating that the same physicochemical properties influence the differential ion mobility behavior. Taken together, this suggests that DMS-MS may capable of predicting the physicochemical properties of discrete lipid molecular species in complex lipidomes.
Acknowledgments
S.E.H. and T.W.M. acknowledge the use of the Mass Spectrometry User Resource and Research Facility within the School of Chemistry (University of Wollongong). B.L.J.P. and S.J.B. acknowledge the support of the Central Analytical Research Facility operated by the Institute for Future Environments (Queensland University of Technology).
The International Workshop on Meibomian Gland Dysfunction: Report of the Subcommittee on Tear Film Lipids and Lipid–Protein Interactions in Health and Disease.
Quantitative structural multi-class lipidomics using differential mobility: electron impact excitation of ions from organics (EIEIO) mass spectrometry.
Comprehensive lipidomic analysis of human plasma using multidimensional liquid- and gas-phase separations: two-dimensional liquid chromatography–mass spectrometry vs. liquid chromatography–trapped-ion-mobility–mass spectrometry.
R Core Team. 2013. R: A language and environment for statistical computing. R Foundation for Statistical Computing. Available from https://www.r-project.org.
This work was supported by Australian Research Council Linkage Project GrantLP140100711 (with industry support from Allergan), Discovery Project Grants DP150101715 and DP190101486 and Future Fellowship Grant FT110100249 (T.W.M.).
The online version of this article (available at http://www.jlr.org) contains a supplement.
32295-1&id=gr1.jpg){kind=link}
32295-1&id=gr2.jpg){kind=link}
32295-1&id=gr4.jpg){kind=link}
32295-1&id=gr3.jpg){kind=link}
32295-1&id=gr5.jpg){kind=link}
32295-1&id=gr6.jpg){kind=link}