Introduction

Obesity-related illness is an increasingly important global health issue that places a tremendous economic burden on society1. The negative health effects of prolonged obesity are partly fuelled by chronic low-grade inflammation, which contributes to cardiometabolic and kidney pathophysiology2,3,4. However, the exact mechanisms that link obesity with cardiometabolic and kidney diseases are unclear and remain a subject of intensive research. The search for biomarkers that assist in the identification of novel disease-related pathways is critical to develop new therapies that are tailored to subpopulations particularly prone to obesity-related pathophysiology.

Disease-related biomarkers are often identified and quantified in blood-derived plasma or serum5,6. Preparation of plasma and serum requires the removal of cellular components by centrifugation. Generation of plasma is preceded by the addition of an anti-coagulant (e.g. EDTA, heparin or citrate) to the whole blood. By contrast, the blood used for serum is allowed to clot before centrifugation, resulting in lower concentrations of clotting factors (such as fibrinogen and coagulation cascade proteins) in serum than in plasma. The World Health Organization generally recommends using plasma as this more accurately reflects the physiological and/or pathophysiological state of the patient7. However, biomarkers are often reported to have better detectability in serum8 despite the fact that serum has a slightly lower total protein concentration than plasma9. Indeed, some intracellularly stored proteins and lipids are only detectable upon coagulation-induced release from leukocytes and platelets, and serum is preferred in assays detecting, for example, cardiac troponins10,11,12. Importantly, the choice of biofluid is not merely a question of detectability, but it may also affect the conclusions drawn from a study. For example, Alsaif et al. showed that of 16 proteins (identified in either plasma or serum) that were differentially expressed between healthy controls and subjects with bipolar disorder, only two showed differential expression in both serum and plasma13.

The aim of our study was to determine whether the use of plasma or serum would yield different results when screening for obesity-related biomarkers. We analyzed proteins and lipids that have previously been suggested to play a role in obesity-related cardiometabolic diseases in plasma and serum from age- and sex-matched groups of lean and obese humans. Our results show that the use of plasma or serum may have an effect on study outcome when screening for obesity-related biomarkers and we identify key markers that highlight this issue.

Results and Discussion

Detectability of proteins in plasma versus serum

We used four Olink multiplex protein panels (inflammation, cardiometabolic, cardiovascular II, cardiovascular III) selected on the basis of their relevance to obesity-related diseases to measure protein concentrations in plasma and serum from 11 obese subjects and 11 age- and sex-matched lean controls. The characteristics of the human cohort are presented in Table 1. Of the 368 proteins analyzed (10 of which were measured in duplicate panels, see Supplementary Table S1 for the full list), one protein (BDNF) was excluded due to technical issues, nine proteins (IL-1 alpha, IL-2, TSLP, IL-22 RA1, IL-13, TNF, IL-20, IL-33, IFN-gamma) were excluded because they were undetectable in both plasma and serum, and 23 additional proteins were excluded because values were missing in >30% of the samples in all of the four groups (lean plasma, lean serum, obese plasma, obese serum; Supplementary Table S1). Detectability issues with one of the excluded proteins, NT-proBNP, have previously been reported14. In total, 335 proteins were included in the comparative analyses (Supplementary Fig. S1).

Table 1 Summary of cohort demographics.

For the majority of proteins, their concentrations were similar between plasma and serum (Supplementary Fig. S2a,b). After adjusting for multiple comparisons using the stringent Holm-Bonferroni test, we found significantly different concentrations between plasma and serum for 23.5% and 33.4% of proteins in the lean and obese cohorts, respectively [adjusted (adj.) p < 0.05, Fig. 1]. Most of these proteins were present at higher concentrations in serum, which may partly be explained by the clotting-induced volume displacement effect15,16 and by the fact that coagulation elicits release of platelet granules and intracellularly stored cytokines17,18,19. The intracellularly stored protein MCP-1, for example, exhibited significantly higher concentrations in serum compared with plasma (in both the inflammation and the cardiovascular III panels) in the lean and obese groups. Of note, we did not record female menstruation cycle and/or menopausal state, which may affect platelet activation, although conflicting results have been shown20,21,22,23,24,25.

Figure 1
figure 1

Detectability of proteins in plasma versus serum. Heatmaps showing protein biomarkers that exhibited significantly different concentrations in plasma versus serum in (a) lean subjects (n = 11) and (b) obese subjects (n = 11) after adjustment for multiple comparisons using the method of Holm-Bonferroni at adj. p < 0.05. Proteins that are significantly different in only one of the groups (lean or obese) are marked in bold. For proteins that are present in duplicate protein panels, the panel is indicated in parentheses: I, inflammation; CVII, cardiovascular II; and CVIII, cardiovascular III. Relative protein concentrations are reported as z-scores.

A subset of proteins with significantly different concentrations in plasma and serum (including HSP-27, PAR-1, 4E-BP1 and SRC) exhibited lower concentrations in serum (Fig. 1). HSP-27 has been proposed as a biomarker for both cardiometabolic disease and cancer26, although controversial results have been reported27. Of note, a recent study showed that the concentration of HSP-27 increased by about three-fold with just one freeze-thaw cycle in plasma but was more stable in serum28. All of our samples underwent two freeze-thaw cycles, which could explain the higher detectability of HSP-27 in plasma. During coagulation, PAR-1 is cleaved29 and SRC30 and 4-EBP131,32 become prone to degradation through proteolytic pathways, which likely explains their lower concentrations in serum. Furthermore, two proteins, AXIN1 and STK4, passed our cut-off criteria for detection in plasma but not in serum (Supplementary Table S1). Enrichment analysis of all the proteins that were significantly altered between plasma and serum confirmed the enrichment of pathways involved in neutrophil chemotaxis and platelet activation (Supplementary Fig. S2c, Supplementary Table S3).

Sensitivity of plasma versus serum when screening for obesity-related protein biomarkers

Most of the biomarkers that showed significantly different concentrations between the obese and lean groups were present at higher levels in the obese group; however, a small number (including IGFBP-1 and GH) showed lower concentrations in the obese group, with significant differences observed in both plasma and serum (Table 2). The number of proteins with significantly different concentrations between the lean and obese groups was greater in serum (Table 2), in agreement with (1) an earlier study that reported higher sensitivity of serum to detect diabetes-associated differences in metabolite concentration33 and (2) the fact that obesity is associated with higher leukocyte and platelet counts and increased platelet activation34,35. MCP-3 was present at higher concentrations in the obese versus lean group in serum but not plasma (Table 2), and showed low detectability in plasma in both groups (Supplementary Table S1). However, concentrations of PAI-1 were only significantly higher in the obese versus lean group in plasma despite showing higher detectability in serum (Table 2). This difference in sensitivity versus detectability for PAI-1 was confirmed by ELISA (Supplementary Fig. S3). PAI-1 inhibits fibrinolysis and has been proposed to be an important biomarker in cardiometabolic and diabetes research, although, as recently reviewed, conflicting results have been reported36. A possible explanation for this discrepancy, at least in part, may be due to the interchangeable use of plasma versus serum; indeed, studies comparing lean versus obese and/or diabetic groups have reported differences in PAI-1 levels when using plasma37 but not serum38. Coagulation-induced secretion of intracellular PAI-1 is likely responsible for the high serum levels of PAI-1, which may mask the differences between the lean and obese groups.

Table 2 Proteins that exhibited significant differences in concentrations between obese and lean groups in plasma and/or serum.

Protein correlations in plasma versus serum

For the correlation analysis, 316 proteins survived the cut-off criteria (Supplementary Fig. S1). We observed significant correlations between plasma and serum samples for most (68.8%) of the proteins analyzed in the lean and obese groups combined (Table 3), although fewer significant correlations were seen when dividing the cohort into obese and lean (Supplementary Table S4). Of the 10 proteins that were measured in duplicate panels, eight displayed similar correlations between plasma and serum. However, MCP-1 and uPA only showed a significant correlation between plasma and serum in one of the duplicate panels.

Table 3 Correlations of protein concentrations in plasma versus serum in all subjects.

We observed good correlations between plasma and serum samples for leptin (r = 1.00, adj. p < 0.001) and IGFBP-1 (r = 0.98, adj. p < 0.001), which are proteins that exhibited obesity-associated differences in concentration (Fig. 2a,b). Some proteins showed poor correlations, such as PCSK9 (r = 0.27, ns) and FGF-23 (r = 0.43 and 0.64 in the inflammation and cardiovascular II panels, respectively, both ns) (Fig. 2c,d). PCSK9 binds to the receptor for low-density lipoprotein and PCSK9 inhibitors are therefore of intense interest to pharmaceutical companies39,40. Studies interchangeably measure PCSK9 in plasma41,42 and serum43,44, but our result indicates that the choice of biofluid could potentially have a significant impact on the conclusions drawn. Our panels also included the FDA-approved biomarkers KIM-1 and osteopontin, which are used to monitor kidney disease45,46. KIM-1 was well correlated between plasma and serum (r = 0.99, adj. p < 0.001) but osteopontin displayed a poor correlation (r = 0.46, ns) (Fig. 2e,f).

Figure 2
figure 2

Protein correlations in plasma versus serum. Pearson correlations (r) between normalized protein expression (NPX) values for proteins in plasma and serum samples. Each data point is from one individual (open triangles: obese; closed triangles: lean). p values were adjusted by the Holm-Bonferroni multiple comparison test.

Lipids in plasma versus serum, and in lean versus obese groups

We also performed targeted lipidomics of inflammation-related lipids in plasma and serum from the lean and obese groups. Of the 76 lipids analyzed (see Supplementary Table S2), two were excluded as they did not survive the cut-off criteria for the comparative analysis (Supplementary Fig. S4). For most of the lipids, there were no major differences in concentration between plasma and serum (Supplementary Fig. S5). We observed that concentrations of 21.6% of the lipids in the lean cohort and 18.9% of the lipids in the obese cohort were significantly higher in serum than in plasma (after FDR adjustment, adj. p < 0.05); none of the lipids showed lower concentrations in serum (Fig. 3a,b). In total, 73 lipids survived the cut-off for the correlation analyses; we observed significant correlations between plasma and serum for 64% of the analyzed lipids when analyzed in the lean and obese groups combined (Table 4), and fewer significant correlations were seen when dividing the cohort into obese and lean (Supplementary Table S5).

Figure 3
figure 3

Oxylipins in plasma versus serum, and in lean versus obese groups. Heatmaps showing lipids that exhibited significantly different concentrations in plasma versus serum in (a) lean subjects (n = 11) and (b) obese subjects (n = 11) after adjustment for multiple comparisons using the false discovery rate (FDR) test at adj. p < 0.05. Lipids that are significantly different in only one of the groups (lean or obese) are marked in bold. Relative lipid concentrations are reported as z-scores. (c) Lipids that showed significantly different concentrations between the obese and lean groups in plasma and/or serum after FDR adjustment.

Table 4 Correlations of lipid concentrations in plasma versus serum in all subjects.

Four lipids showed significantly different concentrations between the obese and lean groups in plasma and/or serum (Fig. 3c). Concentrations of AEA and 19,20-DiHDoPA were significantly different (higher for AEA and lower for 19,20-DiHDoPA in the obese group) in both plasma and serum, but concentrations of 15-HETE and 11-HETE were significantly different (both higher in the obese group) only in plasma (Fig. 3c).

Concluding remarks

In this study, we investigated whether the use of plasma or serum would yield different results when screening for obesity-related biomarkers. For most of the proteins and lipids, their concentrations showed good correlations between plasma and serum. However, it is important to note that PCSK9 concentrations did not correlate between plasma and serum, indicating that caution must be taken when comparing studies that use different biofluids. Although most of the protein and lipids had similar concentrations in plasma and serum, those that did differ were generally present at higher concentrations in serum. Importantly, we observed significantly higher concentrations of the key disease-associated biomarker PAI-1 in the obese group only in plasma and not in serum, despite the protein showing higher detectability in serum. This result highlights that sensitivity does not necessarily parallel detectability. Furthermore, some obesity-induced changes, for example of MCP-3 concentrations, were only detected in serum. Collectively, these findings show that care should be taken when choosing biofluids for the study of biomarkers, particularly those for which we report differences in sensitivity/detectability between plasma and serum.

Methods

Study participants

We recruited obese subjects [body mass index (BMI) 35–55 kg/m2, aged 18–65 years] from a cohort scheduled to undergo gastric bypass surgery, as well as age- and sex-matched lean subjects (BMI 18.5–24.9 kg/m2). Subjects were excluded if they were taking anti-inflammatory and/or immunosuppressive drugs, currently smoked, or had been diagnosed with significant gastrointestinal disease or inflammatory bowel disease. Study participants were enrolled in accordance with the Helsinki Declaration and provided written informed consent. The study was approved by the Gothenburg Ethical Review Board #682-14 (ClinicalTrials.gov NCT02322073).

Blood collection

Venous blood samples obtained from study participants after an overnight fast were collected in eitherplasma tubes spray coated with K2EDTA (Greiner Bio One) or serum tubes containing inert separator gel and silica particles as clot activator (Greiner Bio One). Plasma samples were centrifuged immediately whereas serum samples were allowed to clot for 30 min at room temperature before centrifugation (10 min at room temperature, 3,000 rpm Hettich EBA200). Samples were snap-frozen in liquid nitrogen and stored at −80 °C until analysis.

Multiplex protein assay

Protein biomarkers were analyzed using the proximity extension assay, using four protein panels (inflammation, cardiometabolic, cardiovascular II and cardiovascular III) (Olink Proteomics, Uppsala, Sweden) at the Clinical Biomarkers Facility at Science for Life Laboratory (Uppsala University, Sweden) according to the manufacturer’s instructions. Briefly, 1 µl plasma or serum was incubated with a mixture of 92 proximity antibody pairs tagged with oligonucleotides in a 96-well plate. In this assay, once a pair of antibodies binds to their corresponding antigens in close proximity, linked oligonucleotides hybridize into double stranded DNA, which is further extended and amplified, and ultimately quantified by high-throughput real-time PCR (BioMark™ HD System, Fluidigm Corporation). To avoid intra-assay variability, plasma and serum samples were analyzed on the same plate.

ELISA

Plasma and serum PAI-1 levels were measured using a commercially available ELISA for Human Total Serpin E1/PAI-1 (#DY9387-05, R&D), according to the manufacturer’s instructions. To ensure that the protein was quantified within the linear range of the standard curve, plasma and serum were diluted 1:100 and 1:500, respectively.

Measurements of oxylipins, endocannabinoids and ceramides

Oxylipins, endocannabinoids, and ceramides in plasma and serum were isolated and quantified using modifications of published protocols47,48,49. Briefly, plasma or serum aliquots (40 µl) were spiked with deuterated oxylipin, endocannabinoid and ceramide surrogates, mixed with butylated hydroxyl toluene and ethylene diamine tetraacetic acid, and extracted with 200 µl isopropanol containing the internal standards 1-cyclohexyl ureido, 3-dodecanoic acid and 1-phenyl ureido 3-hexanoic acid in isopropanol. The homogenate was then centrifuged (10 min, 4 °C, 15,000 g) and the isopropanol supernatant was collected and stored at −20 °C until analysis.

Analytes were separated using a Waters Acquity ultra-performance liquid chromatography (UPLC; Waters, Milford, MA) on a 2.1 mm × 150 mm, 1.7 µm BEH C18 column (Waters) for analysis of oxylipins and endocannabinoids, and 2.1 mm × 150 mm, 1.7 µm BEH C8 column (Waters) for analysis of ceramides. Separated analytes were detected by tandem mass-spectrometry, using electrospray ionization with multi reaction monitoring on an API 6500 QTRAP (Sciex, Redwood City, CA) for oxylipins and endocannabinoids, and an API 4000 QTRAP (Sciex) for ceramides. Analytes were quantified using internal standard methods and 7–9 point calibration curves of authentic standards.

Measurement of non-esterified fatty acids

Non-esterified fatty acids in plasma and serum were isolated and converted to fatty acid methyl esters (FAMEs) as previously reported47. Briefly, plasma or serum aliquots (50 µl) were spiked with lipid class surrogates, mixed with 410 µl isopropanol, followed by 520 µl cyclohexane and 570 µl 0.1 M ammonium acetate. Samples were then centrifuged (5 min, 4 °C, 15,000 g), the upper organic phase was collected, and the remainder was re-extracted with a second 520 µl cyclohexane aliquot. The samples were then dried by vacuum centrifugation and reconstituted in 100 µl toluene and 180 µl methanol. To prepare FAMEs, 280 µl of toluene/methanol extracts were enriched with 20 µl methanol containing 60 µM C15:1n5 and incubated with 45 µl 2 M TMS-diazomethane in hexane (Sigma-Aldrich, St. Louis MO) for 30 min at room temperature. Samples were dried under vacuum and the residue was dissolved in 100 µl hexane containing 4 µM C23:0, which acted as an internal standard. Samples were then stored at −20 °C until analysis.

FAMEs were separated on a 30 m × 0.25 mm × 0.25 µm DB-225 ms column in a 6890 gas chromatogram interfaced with a 5973A mass selective detector (Agilent Technologies, Santa Clara, CA). All fatty acids were quantified against a 7-point calibration curves of authentic standards. Peak identifications were based on retention times and m/z ratios, with peak confirmation by inspection of simultaneously acquired full scan spectra collected from 50–400 m/z. Calibrants and internal standards were purchased from NuchekPrep (Elysian, MN), Sigma-Aldrich, or Avanti Polar Lipids. Data were quantified using Chemstation vE.02.14 (Agilent Technologies) against 6–8 point calibration curves.

Statistical analysis

Data are reported for proteins and lipids that had <30% missing values in: (1) at least one of the four groups (lean plasma, lean serum, obese plasma, obese serum) for the comparative analyses or (2) all of the four groups for the plasma-serum correlations.

Statistical analysis of the protein multiplex data was done in the R environment (version 3.5.1) using packages gplots (3.0.1) and gdata (2.18.0)50. For the proteins reported, missing values were replaced with limit of detection (LOD) values. Hierarchical clustering with Pearson correlation distance and complete linkage confirmed that the dataset did not include outliers. Concentrations of proteins are reported as normalized protein expression (NPX) values, an arbitrary unit on a log2 scale. Heatmaps were generated using hierarchical clustering based on correlation distance and Ward’s (ward.D2) clustering. Comparisons of protein levels using Student’s t-test were paired when comparing individual donor plasma versus serum values and unpaired when comparing the lean versus obese groups; p values were adjusted for multiple comparisons using either the stringent Holm-Bonferroni test or the commonly used false discovery rate (FDR) test as indicated (adjusted p values < 0.05 were considered significant). Pearson coefficient of correlation (r) values were calculated and p values were adjusted by the Holm-Bonferroni multiple comparison test.

The pathway enrichment analysis for proteins was done using Metascape51. Briefly, Gene IDs corresponding to significantly altered proteins were analysed, using the 325 unique proteins that survived the cut-off criteria as the background list. A Gene Ontology category was deemed significantly enriched if the p value was lower than 0.01 and displayed a minimum enrichment of 1.5.

Statistical analysis of the lipidomics data was done in MetaboAnalyst52. For the lipids reported, missing values were replaced with half of the lowest reported value. Fatty acid data normalization was optimized in Jmp Pro v 12.0 and confirmed using the Shapiro-Wilk normality test. For statistical analysis, data points underwent log transformation and pareto scaling. Heatmaps were generated using hierarchical clustering based on Euclidean’s method of distance calculation and Ward’s clustering. Unadjusted p values were adjusted using FDR (adjusted p values < 0.05 were considered significant).