Feb 29, 2024
LC/MS analysis of plasma samples from PPMI
- John Chen1,
- Romeo Maciuca1,
- Paul Benton1,
- Jung Suh1,
- Sonnet Davis1,
- Sarah Huntwork-Rodriguez1,
- Ning Xia2,3,
- Michael A. Schwarzschild2,3
- 1Denali Therapeutics Inc, San Francisco, CA, USA;
- 2Molecular Neurobiology Laboratory, Department of Neurology, Massachusetts General Hospital, Boston, MA, USA;
- 3Aligning Science Across Parkinson’s (ASAP) Collaborative Research Network
- Metabolomics Protocols & WorkflowsTech. support email: bbmisraccb@gmail.com
Protocol Citation: John Chen, Romeo Maciuca, Paul Benton, Jung Suh, Sonnet Davis, Sarah Huntwork-Rodriguez, Ning Xia, Michael A. Schwarzschild 2024. LC/MS analysis of plasma samples from PPMI. protocols.io https://dx.doi.org/10.17504/protocols.io.81wgbxe2ylpk/v1
Manuscript citation:
License: This is an open access protocol distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited
Protocol status: Working
We use this protocol and it's working
Created: January 26, 2024
Last Modified: February 29, 2024
Protocol Integer ID: 94181
Keywords: LC/MS analysis, plasma samples, PPMI
Funders Acknowledgements:
Aligning Science Across Parkinson’s (ASAP)
Grant ID: Grant ID: ASAP-000312
Abstract
Plasma samples from PPMI were analyzed by liquid chromatography with mass spectrometry (LC/MS) for a variety of metabolites (including piperine) and lipids as interrelated markers of Parkinson's disease and its pathophysiology.
Plasma sample preparation for LC/MS analysis
Plasma sample preparation for LC/MS analysis
Thawed plasma samples were spun down at 3.5k*g for 10 min at 4°C.
Plasma (10 µL) samples were transferred to an Agilent Captiva collection plate and processed via the Agilent Bravo Metabolomics Sample Prep Platform (Agilent Technologies, Santa Clara, CA).
Bravos Low Volume Plasma Metabolite/Lipid protocols were applied. 1:1 methanol/ethanol (v/v) (112.5 µL) containing internal standards were added into each.
Plates were shaken for 1 min at 1000 rpm and incubated for 10 min at room temperature.
Water (82.5 µL) was added to each well, shaken for 1 min at 500 rpm, and incubated for 10 min at room temperature.
The sample (200 µL) was transferred to a Captiva EMR-Lipid plate to collect metabolite filtrate and remove protein and lipids.
The Captiva EMR-Lipid plate was washed twice with 250 µL of 1:1:1 water/ethanol/methanol (v/v/v).
The metabolite filtrate was dried under N2 gas for 4 h then reconstituted with 200 µL methanol.
To collect the lipids retained on the Captiva EMR-Lipid plate, 1:2 dichloromethane/methanol (v/v) (1.8 mL) was added to the Captiva EMR-Lipid plate for elution.
Lipid filtrate was collected in a 2mL glass coated microplate, dried under N2 gas for 4 h, and reconstituted with 200 µL methanol.
For analysis of the GlcCer/GalCer panel, a 25 µL aliquot from the methanol lipid fraction was dried under N2 gas for 4 h and resuspended in 50 μL of 92.5/5/2.5 acetonitrile/isopropanol/water) with 5mM ammonium formate and 0.5% formic Acid.
Data reporting: Targeted LC-MS/MS data analysis
Data reporting: Targeted LC-MS/MS data analysis
Both targeted and untargeted metabolomics and lipidomics analysis were analyzed in 8 different batches.
Samples included in each batch were selected such that independent factors and covariates of interests were randomized and evenly distributed between batches.
Additionally, 3 pooled control and 3 pooled PD case plasma samples all samples were created by PPMI and these were distributed evenly within and between batches.
Each of the pooled samples were injected 3-6 times within each batch and values from replicate extractions and injections were used to calculate assay CV and inter-batch deviations.
Missing peak areas results were identified and visualized using the naniar r package (version 0.6.1).
Missing values were imputed using k- nearest neighbors imputation method using aggregated information generated from 5 closest donor values.
Imputations were performed using the VIM r package (version 6.1.1).
Following data imputation, areas for endogenous metabolites and lipids were divided by areas
of analyte specific spiked surrogate stable isotope internal standards (See method supplement table) and converted into area ratios.
Subsequently, inter-batch variances of area ratios were corrected by using the Combat function provided in the sva r package (version 3.38).
Relative log expression plots of area ratios were plotted to visualize variations within and across
batches. Following correction for between batch effects, imputed values were removed.
For each non-imputed quantified analyte, un-normalized peak area (UNITS = area), ratios of endogenous metabolite area to surrogate internal standards (UNITS = area_ratio) and batch adjusted area ratios (UNITS = adjusted_area_ratio) were reported.
Batch adjusted area ratios are recommended for downstream analysis. However, unnormalized areas and area ratios provided should allow for different normalization schemes to be adopted by end data consumers.
Note that analytes with missing values in >70% of the samples were not reported.
Data analysis and reports were generated using R version 4.0.2.
Untargeted LC-MS data analysis
Untargeted LC-MS data analysis
Data files were extracted from Thermo ‘raw’ format to mzML using proteowizard msconvert
(version-3.0.21334).
A table of files and batch number was produced using the file name nomenclature.
This was feed into the apLCMS algorithm using the latest update (version 6.6.8) for multi-batch
processing (two.step.hybrid), as described in Liu et al.(1)
Briefly, this method performs peak detection on each file within a batch, while also
gain accuracy and power by using knowledge from other samples within the batch.
Once peak detection is finished retention time alignment is performed.
The algorithm finds a medium retention time deviation sample within the batch and
uses this sample to align the other samples to this medium sample.
When all the batches have been processed a pairwise retention time alignment (similar to
the above method) and weak-signal recovery is performed across the batches.
The weak-signal recovery adds accuracy to downstream data analysis while also
adding alignment across batches.
Protocol references
Liu, Q., Walker, D., Uppal, K.et al. Addressing the batch effect issue for LC/MS metabolomics data in data preprocessing. Sci Rep10,13856 (2020). https://doi.org/10.1038/s41598-020-70850-0