Novel plasma biomarkers improve discrimination of metabolic health independent of weight

ABSTRACT We sought to determine if novel plasma biomarkers improve traditionally defined metabolic health (MH) in predicting risk of cardiovascular disease (CVD) events irrespective of

weight. Poor MH was defined in CATHGEN biorepository participants (n > 9300), a follow-up cohort (> 5600 days) comprising participants undergoing evaluation for possible ischemic heart

disease. Lipoprotein subparticles, lipoprotein-insulin resistance (LP-IR), and GlycA were measured using NMR spectroscopy (n = 8385), while acylcarnitines and amino acids were measured

using flow-injection, tandem mass spectrometry (n = 3592). Multivariable Cox proportional hazards models determined association of poor MH and plasma biomarkers with time-to-all-cause

mortality or incident myocardial infarction. Low-density lipoprotein particle size and high-density lipoprotein, small and medium particle size (HMSP), GlycA, LP-IR, short-chain

dicarboxylacylcarnitines (SCDA), and branched-chain amino acid plasma biomarkers were independently associated with CVD events after adjustment for traditionally defined MH in the overall

cohort (_p_ = 3.3 × 10−4–3.6 × 10−123), as well as within most of the individual BMI categories (_p_ = 8.1 × 10−3–1.4 × 10−49). LP-IR, GlycA, HMSP, and SCDA improved metrics of model fit

analyses beyond that of traditionally defined MH. We found that LP-IR, GlycA, HMSP, and SCDA improve traditionally defined MH models in prediction of adverse CVD events irrespective of BMI.

SIMILAR CONTENT BEING VIEWED BY OTHERS THE ROLE OF NMR-BASED CIRCULATING METABOLIC BIOMARKERS IN DEVELOPMENT AND RISK PREDICTION OF NEW ONSET TYPE 2 DIABETES Article Open access 05 September

2022 METABOLOMIC SIGNATURES OF IDEAL CARDIOVASCULAR HEALTH IN BLACK ADULTS Article Open access 20 January 2024 LATE NON-FASTING PLASMA GLUCOSE PREDICTS CARDIOVASCULAR MORTALITY INDEPENDENT

OF HEMOGLOBIN A1C Article Open access 11 May 2022 INTRODUCTION Claiming more lives in 2015 than the combined mortality of the second and third leading causes, cancer and chronic lower

respiratory disease, respectively, cardiovascular disease (CVD) remains the leading cause of morbidity and mortality within the United States1. Obesity, a known risk factor for a variety of

negative health effects including CVD and adverse CVD-related events2, remains a major public health challenge; in fact, the age-adjusted percentage of adults greater than 20 years old with

obesity in the US increased from 22.9% in 1988–1994 to 37.8% in 2013–20143. The relationship between obesity and CVD has been partially attributed to the increased prevalence of intermediate

risk factors including insulin resistance, type 2 diabetes mellitus (DM), dyslipidemia, and hypertension. However, the underpinnings remain incompletely elucidated and the heterogeneity of

prevalence of these risk factors in obesity demonstrates the complexity of the relationship4. Metabolic health (MH) has traditionally been defined by the presence or absence of a cluster of

intermediate CVD risk factors known as metabolic syndrome5, while medical literature has also classically defined being overweight and having poor metabolic health as synonymous6. Metabolic

syndrome is a known contributor to the epidemic of CVD7, however, a growing body of evidence suggests that weight is a poor metric of metabolic health4,8 even though obesity continues to be

reported as a “major determinant of metabolic syndrome9”. While intermediate CVD risk factors seen in metabolic syndrome are important, perhaps the most important metric for defining

metabolic health should be its ability to associate with hard cardiovascular endpoints. Poor metabolic health as defined by metabolic syndrome is associated with increased risk of major

adverse cardiovascular events (MACE) in the overall population7, but evidence suggests the risk of MACE is not well correlated with BMI category10, often as a result of underestimating the

risk of those patients who have components of metabolic syndrome but do not meet full criteria11,12. In order to more effectively target interventions to improve metabolic health and prevent

adverse CVD events, there is a need for better markers of metabolic health beyond traditionally defined MH that are predictive regardless of BMI. Plasma biomarkers reflect dysregulated

systemic and tissue-specific metabolism, serving as downstream read-outs of genetic, transcriptomic and proteomic variation, and thus may serve as these better markers. In fact, recent

studies in our lab using metabolomic profiling have identified a number of biomarkers that are associated with obesity, insulin resistance, DM, CVD, and response to interventions, including

low-density lipoprotein (LDL-P)13,14,15, lipoprotein-insulin resistance (LP-IR)16, GlycA17, sum of medium plus small high-density lipoprotein-particle subclass (HMSP)18, short-chain

dicarboxylacylcarnitines (SCDA)19,20, and branched-chain amino acids (BCAA)20,21,22. Thus, we hypothesized that these novel plasma biomarkers can improve predictive capabilities for incident

adverse CVD events across BMI categories and are better markers of metabolic health than traditional measures which often rely heavily on BMI. METHODS CATHGEN biorepository participants (n 

> 9300), a follow-up cohort (> 5600 days) comprising participants undergoing evaluation for possible ischemic heart disease. STUDY POPULATION Individuals undergoing cardiac

catheterization between 2001 and 2010 for the indication of possible ischemic heart disease at Duke University Medical Center provided written consent and enrolled in the CATHGEN

(CATHeterization GENetics) biorepository (n > 9300) as previously described23. For this follow-up cohort, longitudinal annual follow-up for recurrent events, hospitalizations, and vital

status has been performed and confirmed through the National Death Index and the Social Security Death Index for > 5600 days from time of CATHGEN enrollment. CATHGEN is approved by the

Duke University School of Medicine Institutional Review Board in accordance with the Declaration of Helsinki (Durham, NC, USA). Prior to administration of supplemental heparin, 50 mL of

blood was collected at the time of arterial puncture from fasting, consented participants. Samples were immediately cooled to 4 °C, processed to separate plasma, and then frozen to − 80 °C.

Clinical data including demographics, medical comorbidities, angiographic data, quantitation of left ventricular ejection fraction (LVEF), and laboratory data were collected. For this study,

participants were excluded if they were cardiac transplant recipients (n = 201), had undergone cardiac catheterization for the indication of pulmonary hypertension (n = 94), had a history

of congenital heart disease (n = 27), had missing BMI (n = 1), had no longitudinal follow-up data (n = 8), had missing biomarker data (n = 115), or had incomplete data necessary for

calculation of traditionally defined MH (n = 226). Race was classified as white or non-white. LVEF was obtained by ventriculogram at the time of cardiac catheterization or, if

ventriculography was not performed, by an echocardiogram or cardiac magnetic resonance image (MRI). Coronary artery disease was determined by the physician performing the cardiac

catheterization, and defined as present if a stenosis ≥ 50% was present in any major epicardial coronary vessel. For this study, the primary endpoint was defined as time-to-event for

all-cause mortality or MI from CATHGEN enrollment. DEFINITION OF BMI CATEGORIES AND TRADITIONALLY DEFINED POOR MH Individuals were classified into three categories based on BMI at time of

enrollment: lean (BMI < 25 kg/m2), overweight (BMI 25–30 kg/m2), or obese (BMI ≥ 30 kg/m2). Individuals were categorized as having traditionally defined poor MH if they had at least two

out of four criteria for poor MH at baseline enrollment, and considered to have traditionally defined good MH if they met 0–1 criteria (modified slightly from the Adult Treatment Panel III

(ATPIII) guidelines24 given lack of waist circumference data in CATHGEN): (1) fasting triglycerides (TG) ≥ 150 mg/dL, (2) HDL < 40 mg/dL in men or < 50 mg/dL in women, (3) history of

diabetes, insulin use, fasting glucose ≥ 100 mg/dL, or HOMA-IR > 5.13, (4) blood pressure medication use or history of patient reported hypertension. Homeostatic model assessment of

insulin resistance was calculated as follows: HOMA-IR = (fasting insulin in (µUI/mL) * glucose (mM)/22.521,25. A HOMA-IR threshold of > 5.13 (i.e., the 90th percentile) was selected and

consistent with prior work21. METABOLOMIC PROFILING Lipoprotein particle concentrations, sizes, and subclasses as well as GlycA, LP-IR, and creatinine were measured as previously

described17,20,26. For this study, we used a composite sum of medium plus small HDL-particle subclass (diameters ≤ 9.4 nm, termed HMSP) because we have previously reported HMSP as more

strongly associated with incident CVD events than individual HDL components17. GlycA is an inflammatory biomarker that consists of the NMR signal from the _N_-acetyl methyl groups of the

_N_-acetylglucosamine residues on enzymatically glycosylated acute phase proteins (primarily α1-antichymyotrypsin, α1-acid glycoprotein, haptoglobin, α1-antitrypsin, and transferrin)17. The

LP-IR score is a measure of insulin resistance derived using the weighted sum of six NMR derived lipoprotein parameters16; the multiplex validated score ranges from 0 to 100 (where 100 is

the most insulin resistant). For measurement of acylcarnitines and amino acids, targeted metabolic profiling was performed using flow-injection, tandem mass spectrometry (MS/MS) as

previously reported27. Values below the lower limits of quantitation (LLOQ) were reported and analyzed as “0,” while plasma biomarkers with > 25% of values below the LLOQ were not

analyzed. Principal components analysis (PCA) with varimax rotation was used for reduction of metabolic data into identifiable factors, with those factors having an eigenvalue of ≥ 1.0 being

retained, and plasma biomarkers with a factor load of ≥ 0.4 being reported as comprising a factor as previously described20,27,28. In our prior work, a PCA-derived factor composed of SCDA

(glutaryl carnitine (C5-DC), 3-hydroxy-cis-5-octenoyl carnitine or hexenedioyl (C6:1-DC/C8:1-OH), octenedioyl carnitine (C8:1-DC), adipoyl carnitine (C6-DC), methylmalonyl carnitine or

succinyl carnitine (Ci4-DC/C4-DC), 3-hydroxy-decanoyl carnitine or suberoyl carnitine (C10-OH/C8-DC), 3-hydroxy-dodecanoyl carnitine or sebacoyl carnitine (C12-OH/C10-DC), and citrulline) is

predictive of incident CVD events. Thus, that PCA factor was used in this study. Prior work has also shown a PCA-derived factor of BCAA and two aromatic amino acids is associated with CAD

and adverse CVD events; thus, that PCA factor was also used in this study (composed of leucine/isoleucine, valine, phenylalanine, tyrosine, and methionine)27. These PCA-factors for SCDA and

BCAA are used as measures for the composite of SCDA and BCAA in this study. STATISTICAL METHODS Baseline clinical and demographic characteristics were compared between individuals with and

without traditionally defined poor MH using _t_-tests and chi-squared tests for continuous and categorical variables, respectively. Incident events were defined as all-cause mortality or MI

(adjudicated by a central process within the Duke Databank for Cardiovascular Disease). Cox proportional hazards models were used to examine the independent and additive effects of plasma

biomarkers on time-to-event, accounting for traditionally defined poor MH. Specifically, we tested three models for time-to-event, each adjusted for BMI category (lean, overweight, obese) in

the full sample: (1) traditionally defined poor MH; (2) each individual plasma biomarker alone (LDL-P, LP-IR, GlycA, HMSP, BCAA, and SCDA), adjusted for traditionally defined poor MH; and

(3) a full model inclusive of traditionally defined poor MH and all plasma biomarkers jointly. We also stratified the data by BMI category, and tested the same models to understand the

potential heterogeneity of effects of the relationship of traditionally defined poor MH and novel plasma biomarkers across BMI categories. To evaluate medications as a confounder, we

performed sensitivity analyses on a subset of individuals with available medication data, where time-to-event models were adjusted for statin use. Statin use was determined using EHR

admission medications, i.e. an individual was considered to be on a statin if they were on a statin at admission either 365 days prior to or 10 days after the index cardiac catheterization.

Finally, while our primary goal was to identify better markers of metabolic health compared to traditionally defined MH parameters (using association with incident events as the metric to

determine “better”), in sensitivity analyses, we also created models that included traditional CVD risk factors (age, sex, race, smoking, ejection fraction, LDL-C, creatinine, family history

of CAD, and presence/absence of CAD). The assumption of proportional hazards was tested for each model; when violated, we included a term in the model with a time-varying coefficient and

compared the resulting inference to inference based on the original model. We used the Akaike Information Criterion (AIC) and area under the receiver operating characteristic curve (AUC) to

compare model fit between different models and determine the incremental discriminative capability of novel plasma markers. AIC analyses were restricted to individuals who had data needed

for all models of traditionally defined MH and for all biomarkers (n = 3222). RESULTS Table 1 shows baseline characteristics of the overall study population (n = 8671 individuals) stratified

by both metabolic health and BMI categories; 1921 (22.2%) individuals were classified as lean, 3068 individuals (35.4%) were overweight, and 3682 individuals (42.5%) were obese. Obese

individuals tended to be younger than overweight or lean individuals (59.2 ± 11.2 years vs. 62.8 ± 11.7 years vs. 64.1 ± 13.1 years, respectively) and had greater prevalence of races other

than white. Approximately 45% of obese and lean individuals were female compared with only 30% of overweight individuals. Overall, there was a high prevalence of traditionally defined poor

MH (71.3%). As expected, obese and overweight individuals had a greater prevalence of poor MH compared with lean individuals (81.5% vs. 69.1% vs. 55.3%), although the prevalence of poor MH

was high even in the lean group. Concurrently, 18.5% of obese individuals were defined as having good MH by traditional parameters, further confirming the known disconnect between metabolic

health and weight. The increasing prevalence of poor MH across increasing BMI categories was driven by all components of the metabolic health measure including higher prevalence of

hypertension and diabetes, higher TG, and lower HDL. Interestingly, LDL was lower in poor MH individuals as compared with good MH individuals across each of the three BMI categories (96.0

vs. 102.9 mg/dL lean, 101.3 vs. 106.7 mg/dL overweight, 101.4 vs. 111.9 mg/dL obese, poor vs good MH individuals respectively). In the overall cohort, individuals with poor metabolic health

had a higher prevalence of a positive family history of CAD as compared with individuals with good metabolic health (34.4% vs. 29.9%, _p_ < 0.001), which was also true in the overweight

category (35.2% vs. 31.3%, _p_ < 0.05). However, there were no significant differences in the lean and obese group. ASSOCIATION OF BMI AND TRADITIONALLY DEFINED POOR MH WITH INCIDENT CVD

EVENTS There were 3845 death or MI events (44.3% of individuals; 6.5% MI, 37.8% deaths) with a mean time to event of 2433 days (SD 1263 days). Traditionally defined poor MH was associated

with time-to-death/MI in the overall cohort (HR 1.38, 95% CI 1.28–1.49, p = 5.0 × 10−17, Fig. 1A). In the same models stratified by BMI category, traditionally defined poor MH was associated

with a stepwise increase in HR for increasing BMI category: lean: HR 1.23, 95% CI 1.08–1.39, _p_ = 1.4 × 10−3; overweight: HR 1.45, 95% CI 1.28–1.64, _p_ = 7.2 × 10−9; obese: HR 1.53, 95%

CI 1.32–1.77, _p_ = 2.0 × 10−8. 1028 events occurred in the lean population (54% of lean individuals), with 411 events occurring in lean individuals with traditionally defined good MH.

Interestingly, after accounting for smoking history, obesity was protective for events, with higher BMI categories associated with a lower risk of events (overweight: HR: 0.71, 95% CI

0.65–0.77; obese: HR 0.65, 95% CI 0.59–0.70, both compared to referent lean BMI category) (Fig. 1B). NOVEL PLASMA BIOMARKERS PREDICT INCIDENT EVENTS INDEPENDENT OF TRADITIONALLY DEFINED POOR

MH AND ACROSS ALL BMI CATEGORIES To understand the prognostic capability of novel plasma biomarkers independent of traditional metabolic health parameters, we analyzed Cox models for each

plasma biomarker with time-to-death/MI, adjusting for traditionally defined poor MH (Table 2). In these analyses, each plasma biomarker was significantly associated with time-to-death/MI in

the overall cohort; GlycA and SCDA were associated with decreased survival, while LDL-P, LP-IR, HMSP and BCAA were all associated with increased survival. In analyses stratified by BMI

category, each individual plasma biomarker remained significantly predictive of death/MI and with similar effect sizes across each BMI category except BCAA, which was associated with

time-to-event only in obese individuals. Models for plasma biomarkers LDL-P, LP-IR, GlycA, and HMSP included 8385 individuals, and those for BCAA and SCDA included 3591 individuals.

Sensitivity analyses including individuals with medication data (N = 3,174 out of 8761 [36.2%]) were conducted including statin use as a binary covariable (Supplementary Table 1). Overall,

the magnitude and significance of results remained the same in these analyses. In a full Cox model inclusive of traditionally defined poor MH and all plasma biomarkers, each individual

plasma biomarker remained independently significant with a similar magnitude of effect in the overall cohort (Supplementary Table 2). In similar full models stratified by BMI category, each

plasma biomarker other than LP-IR and LDL-P also remained significantly predictive of events in one or more obesity categories, with GlycA, HMSP, and SCDA maintaining significance across all

BMI categories (Supplementary Table 2). The proportional hazards assumption was violated in our models for plasma biomarkers LP-IR, GlycA, and HMSP in the overall population and the

BMI-stratified analyses. To examine the potential impact of non-proportional hazards related to these plasma biomarkers, we conducted sensitivity analyses allowing for a time-varying

coefficient in the model. All plasma biomarkers were still significantly associated with events in all tested models (_p_ < 10−12). Estimates of the HR at various timepoints are given in

Supplementary Table 3; these demonstrate that the magnitude of the estimated HR for each plasma biomarker increases when the time-varying coefficient is included in the model, compared to

the results in Table 2. Therefore, for simplicity in the remaining analyses and with the understanding that this may lead to an underestimation of the effect sizes for LP-IR, GlycA, and

HMSP, we utilized a Cox model without time-varying coefficients. NOVEL PLASMA BIOMARKERS INDEPENDENTLY PREDICT EVENTS IN FULL MULTIVARIABLE MODELS ADDITIONALLY ADJUSTED FOR CLINICAL RISK

FACTORS While our primary goal was to understand novel plasma biomarkers as markers of metabolic health using event prediction as a benchmark for assessing importance, there are also other

known clinical risk factors for incident CVD events that are not captured by models of traditionally defined poor MH; thus, we also created full multivariable models, including age, sex,

race, LVEF, smoking, LDL-C, creatinine, and presence/absence of CAD. In these models, which also adjusted for traditionally defined poor MH, each individual novel plasma biomarker except

BCAA remained predictive of events in the overall population (Table 2). In BMI category stratified analyses, LP-IR, GlycA, and HMSP plasma biomarkers were significantly and consistently

associated with incident CVD events across all BMI categories with a similar magnitude of effect in each category, suggesting they are important biomarkers regardless of BMI. Traditionally

defined poor MH remained predictive overall and in the overweight and obese categories, but was not significant in its prediction of events in the lean group (Table 2). In full multivariable

models inclusive of traditionally defined poor MH, CVD clinical risk factors, and all plasma biomarkers, four plasma biomarkers were significant in the overall cohort: GlycA (HR 1.23, 95%

CI 1.17–1.29, _p_ = 3.6 × 10−18); HMSP (HR 0.76, 95% CI 0.71–0.81, _p_ = 2.0 × 10−18); BCAA (HR 0.93, 95% CI 0.88–0.99, _p_ = 0.01); and SCDA (HR 1.17, 95% CI 1.08–1.28, _p_ = 2.7 × 10−4).

In analyses stratified by BMI category, only GlycA and HMSP remained significantly predictive of events across all BMI categories (Supplementary Table 4). NOVEL PLASMA BIOMARKERS

INCREMENTALLY IMPROVE CVD EVENT PREDICTION To further assess the incremental predictive capabilities of each novel plasma biomarker, we compared models using the AIC and AUC. These models

demonstrated modest improvement of model fit with addition of any individual novel plasma biomarker to traditionally defined poor MH in the overall cohort: traditionally defined poor MH (AIC

20,240, AUC 0.58), LDL-P (AIC 20,222, AUC 0.59), LP-IR (AIC 20,179, AUC 0.61), GlycA (AIC 20,162, AUC 0.63), HMSP (AIC 20,051, AUC 0.64), BCAA (AIC 20,232, AUC 0.61), SCDA (AIC 20,155, AUC

0.62) (Table 3). BMI stratified categories showed similar improvement of model fit with addition of any individual novel plasma biomarker other than LDL-P (Table 3). When all plasma

biomarkers are included, there is greater improvement in model fit than with any single variable alone in the overall cohort (AIC 19,919, AUC 0.69) and within each BMI category (Table 3). In

models inclusive of clinical covariates, we found modest improvement of model performance beyond that of traditionally defined poor MH, with further improvement in performance with addition

of any individual plasma biomarker other than LDL-P in the overall cohort: traditionally defined poor MH (AIC 19,674, AUC 0.74), LDL-P (AIC 19,673, AUC 0.74), LP-IR (AIC 19,658, AUC 0.74),

GlycA (AIC 19,598, AUC 0.75), HMSP (AIC 19,560, AUC 0.75), BCAA (AIC 19,672, AUC 0.75), SCDA (AIC 19,649, AUC 0.75) (Table 3). BMI stratified categories showed improvement of model fit with

addition of GlycA, HMSP, BCAA, and SCDA. Similar to prior, the greatest improvement in model performance occurs in the model inclusive of traditionally defined poor MH, all plasma

biomarkers, and clinical covariates (AIC 19,481, AUC 0.77) (Table 3). We also explored the potential association between the plasma biomarkers and family history of CAD (Supplementary Table

5) and found significant correlations between LP-IR (β = 0.20, _p_ = 1.18 × 10−17), LDLP (β = 0.07, _p_ = 0.002) and HMSP (β = 0.12, _p_ = 7.49 × 10−7) in the overall cohort. In the lean

group, only LP-IR (β = 0.25, _p_ = 1.82 × 10−8) and HMSP (β = 0.16, _p_ = 0.002) were associated with family history of CAD. DISCUSSION In this large cardiovascular cohort study, we

demonstrate that novel plasma biomarkers discovered through prior studies of CVD may serve as better markers of poor metabolic health as determined by their ability to predict risk of

adverse cardiovascular events across all BMI categories, independent of, and incremental to, traditionally defined poor MH parameters. Our most striking results include: (1) traditionally

defined poor MH does not predict CVD events in lean individuals in full multivariable models, but the majority of novel plasma biomarkers do; (2) GlycA and HMSP plasma biomarkers had strong

predictive capabilities in the overall cohort and across all BMI categories even after adjusting for traditionally defined poor MH and other clinical risk factors, and with the greatest

improvements in AUC on top of traditionally defined poor MH; (3) adding all novel plasma biomarkers to traditionally defined poor MH improved the AUC from 0.58 to 0.69 in the overall cohort,

and with a similar magnitude of improvement in lean and overweight BMI categories where metabolic health is perhaps the most heterogeneous. In our study, 55.3% of lean individuals met

criteria for poor MH, while 18.5% of obese individuals met criteria for having good MH. Furthermore, of the 3845 total events, 24.3% occurred in individuals with good MH and 26.7% occurred

in lean individuals. These highlight the disconnect between traditionally defined poor MH, weight, and risk of morbidity and mortality. We used hard outcomes of death or incident MI to

determine utility of traditional and novel risk factors. Overall, all of the novel plasma biomarkers assessed in this study other than BCAA were found to be predictive of incident CVD

events, even after adjusting for traditionally defined poor MH and additional clinical risk factors. This predictive capability held true in analyses stratified by BMI category for plasma

biomarkers LP-IR, GlycA, and HMSP. Furthermore, GlycA and HMSP were strongly predictive of events in the lean group, which is classically felt to be of lesser adverse cardiovascular event

risk but may actually have a higher risk. Importantly, novel plasma biomarkers also improved model fit on top of traditionally defined poor MH and additional clinical risk factors,

suggesting that they have incremental risk-predictive capabilities beyond what is easily assessed clinically. Serum concentrations of GlycA, a novel inflammatory biomarker that is an

aggregate measure of enzymatically glycosylated, acute-phase proteins, have been shown to independently predict incident CVD events in the MESA study (HR 1.33, 95% CI 1.24–1.42) and in the

JUPITER trial (HR 1.33, 95% CI 1.21–1.45)29,30, and with presence and extent of CAD17. The biologic relevance of elevations in GlycA seen in these studies may be due to effects on

glycosylation of HDL components resulting in impaired function17. In full multivariate models, a factor composed of SCDA was robust for prediction of death or MI in the overall cohort as

well as in lean and overweight individuals. These plasma biomarkers have been shown to be associated with increased risk of adverse CVD related events in the CATHGEN cohort19,20, and appear

to be reporting on dysregulated endoplasmic reticulum stress28,31. The role of HDL cholesterol in CVD is complex; studies have demonstrated an inverse relationship of HDL-C with CVD, but

pharmacologic intervention to increase HDL-C have failed to show improved outcomes32. A growing number of studies, however, are showing that NMR-derived lipoprotein particle number more

closely correlates with risk of adverse CVD outcomes than traditional LDL-C and HDL-C measures18,33. The cardioprotective effects of HMSP seen in the current study likely depend on specific

HDL subparticle characteristics that are not well accounted for by traditional HDL-C, which is preferentially affected by larger, cholesterol rich particles32; similar to that seen in a

recent study showing that the sum of medium- and small-size HDL particles was the strongest predictor of mortality18. In the current study, higher LDL-P levels were protective against death

or MI. From prior work with the CATHGEN population, we have seen evidence of a similar lipid paradox13,14,15,20. The etiology of this phenomenon remains incompletely elucidated and may be

related to medication use (e.g. patients with higher LDL levels are treated earlier and for longer with statins) but new evidence might suggest that nutritional status plays a key role34.

Unexpectedly, while HOMA-IR was higher with increasing BMI, higher LP-IR levels were protective of death or MI. This might be related to the observed obesity paradox in our secondary

prevention dataset, i.e. because obesity is protective of adverse outcomes and obese individuals have higher HOMA-IR. In addition, individuals with poor metabolic health had a higher

prevalence of a positive family history of CAD as compared with individuals with good metabolic health which was also true in the overweight category. However, there were no significant

differences in the lean and obese group, suggesting that this more simple screening tool would not be sufficient for defining CVD risk. Moreover, we found correlations between LP-IR, LDLP

and HMSP and family history of CAD in the overall cohort. In the lean group, only LP-IR and HMSP were associated with family history of CAD. This is as predicted given the known heritability

of these markers. Regarding the clinical utility of our findings and testing for these plasma biomarkers, it is relatively quick, reproducible, and of course, by their nature,

non-invasive35. Also, mass spectrometry is currently already in use for new-born screening and more severe metabolic phenotypes (e.g. phenylketonuria, maple-syrup urine disease, etc.).

Although the use of some of these plasma biomarkers is still in its infancy, there is a growing body of evidence that supports their clinical utility in cardiovascular disease35. Of note,

the NMR-derived lipoproteins are commonly used in clinical practice already. Several strengths of our study deserve mention including the relatively large sample size, availability and

length of follow-up data for hard CVD outcomes, fasting samples collected under a rigorous and similar protocol, and the breadth of clinical variables available. Importantly, many prior

studies of metabolic health measures do not have a “gold standard” to tag improvement of models to; the use of a clinically relevant outcome of death or MI is an important metric that we

were able to use to compare models. However, since this is a cross-sectional, observational study, we were unable to account for chronic dietary and environmental variations in plasma

biomarker levels, nor did we have sufficient medication data to enable adjustment for medication use (although previous work showed such adjustments had minimal influence on plasma biomarker

levels)27. Finally, although the model of metabolic health improved with the addition of clinical covariates to the plasma biomarkers, the increase was modest, which suggests that clinical

covariates contribute more to the models rather than the biomarkers. However, we note that the greatest improvement was seen in the lean group, a group that is often dismissed clinically as

having low CVD risk, and as such, biomarkers may have the greatest clinical utility for identifying higher risk individuals in this group. In conclusion, with growing evidence of the

discordance between metabolic health and BMI, it is crucial to improve clinical models for early recognition and characterization of individuals at high risk for CVD events, and thus at

greatest need for early intensive therapeutic intervention. We found plasma biomarkers had the potential to improve clinical models, to better classify overall metabolic health, and to

independently predict events, irrespective of obesity stratum. This was especially striking given that traditionally defined MH models were unable to predict events in the lean group, while

novel plasma biomarkers were successful in prediction of events; this suggests that novel plasma biomarkers may be better markers of metabolic health in this patient population than current

standard clinical models. REFERENCES * Benjamin, E. J. _et al._ Heart disease and stroke statistics-2018 update: a report from the American heart association. _Circulation_ 137, e67–e492.

https://doi.org/10.1161/cir.0000000000000558 (2018). Article  PubMed  Google Scholar  * Guh, D. P. _et al._ The incidence of co-morbidities related to obesity and overweight: a systematic

review and meta-analysis. _BMC Public Health_ 9, 88. https://doi.org/10.1186/1471-2458-9-88 (2009). Article  PubMed  PubMed Central  Google Scholar  * National Center for Health, S. in

_Health, United States, 2016: With Chartbook on Long-term Trends in Health_ (National Center for Health Statistics (US) (2017). * Wildman, R. P. _et al._ The obese without cardiometabolic

risk factor clustering and the normal weight with cardiometabolic risk factor clustering: prevalence and correlates of 2 phenotypes among the US population (NHANES 1999–2004). _Arch. Intern.

Med._ 168, 1617–1624. https://doi.org/10.1001/archinte.168.15.1617 (2008). Article  PubMed  Google Scholar  * Cornier, M. A. _et al._ The metabolic syndrome. _Endocr. Rev._ 29, 777–822.

https://doi.org/10.1210/er.2008-0024 (2008). Article  CAS  PubMed  PubMed Central  Google Scholar  * Ding, C., Chan, Z. & Magkos, F. Lean, but not healthy: the “metabolically obese,

normal-weight” phenotype. _Curr. Opin. Clin. Nutr. Metab. Care_ 19, 408–417. https://doi.org/10.1097/mco.0000000000000317 (2016). Article  CAS  PubMed  Google Scholar  * Malik, S. _et al._

Impact of the metabolic syndrome on mortality from coronary heart disease, cardiovascular disease, and all causes in United States adults. _Circulation_ 110, 1245–1250.

https://doi.org/10.1161/01.Cir.0000140677.20606.0e (2004). Article  PubMed  Google Scholar  * Samaranayake, N. R., Ong, K. L., Leung, R. Y. & Cheung, B. M. Management of obesity in the

National Health and Nutrition Examination Survey (NHANES), 2007–2008. _Ann. Epidemiol._ 22, 349–353. https://doi.org/10.1016/j.annepidem.2012.01.001 (2012). Article  PubMed  Google Scholar 

* Sherling, D. H., Perumareddi, P. & Hennekens, C. H. Metabolic syndrome. _J. Cardiovasc. Pharmacol. Ther._ 22, 365–367. https://doi.org/10.1177/1074248416686187 (2017). Article  CAS 

PubMed  Google Scholar  * Kip, K. E. _et al._ Clinical importance of obesity versus the metabolic syndrome in cardiovascular risk in women: a report from the Women’s Ischemia Syndrome

Evaluation (WISE) study. _Circulation_ 109, 706–713. https://doi.org/10.1161/01.Cir.0000115514.44135.A8 (2004). Article  PubMed  Google Scholar  * Ahmadi, A. _et al._ Is metabolic syndrome

predictive of prevalence, extent, and risk of coronary artery disease beyond its components? Results from the multinational coronary CT angiography evaluation for clinical outcome: an

international multicenter registry (CONFIRM). _PLoS ONE_ 10, e0118998. https://doi.org/10.1371/journal.pone.0118998 (2015). Article  CAS  PubMed  PubMed Central  Google Scholar  * Cavallari,

I. _et al._ Metabolic syndrome and the risk of adverse cardiovascular events after an acute coronary syndrome. _Eur. J. Prev. Cardiol._ 25, 830–838. https://doi.org/10.1177/2047487318763897

(2018). Article  PubMed  Google Scholar  * Reddy, V. S. _et al._ Relationship between serum low-density lipoprotein cholesterol and in-hospital mortality following acute myocardial

infarction (The Lipid Paradox). _Am. J. Cardiol._ 115, 557–562. https://doi.org/10.1016/j.amjcard.2014.12.006 (2015). Article  CAS  PubMed  Google Scholar  * Horwich, T. B., Hamilton, M. A.,

MacLellan, W. R. & Fonarow, G. C. Low serum total cholesterol is associated with marked increase in mortality in advanced heart failure. _J. Cardiac. Fail._ 8, 216–224.

https://doi.org/10.1054/jcaf.2002.0804216 (2002). Article  CAS  Google Scholar  * Rauchhaus, M. _et al._ The relationship between cholesterol and survival in patients with chronic heart

failure. _J. Am. Coll. Cardiol._ 42, 1933–1940. https://doi.org/10.1016/j.jacc.2003.07.016 (2003). Article  CAS  PubMed  Google Scholar  * Shalaurova, I., Connelly, M. A., Garvey, W. T.

& Otvos, J. D. Lipoprotein insulin resistance index: a lipoprotein particle-derived measure of insulin resistance. _Metab. Syndr. Relat. Disord._ 12, 422–429.

https://doi.org/10.1089/met.2014.0050 (2014). Article  CAS  PubMed  PubMed Central  Google Scholar  * McGarrah, R. W. _et al._ A novel protein glycan-derived inflammation biomarker

independently predicts cardiovascular disease and modifies the association of HDL subclasses with mortality. _Clin. Chem._ 63, 288–296. https://doi.org/10.1373/clinchem.2016.261636 (2017).

Article  CAS  PubMed  Google Scholar  * McGarrah, R. W. _et al._ High-density lipoprotein subclass measurements improve mortality risk prediction, discrimination and reclassification in a

cardiac catheterization cohort. _Atherosclerosis_ 246, 229–235. https://doi.org/10.1016/j.atherosclerosis.2016.01.012 (2016). Article  CAS  PubMed  PubMed Central  Google Scholar  * Shah, A.

A. _et al._ Metabolic profiles predict adverse events after coronary artery bypass grafting. _J. Thorac. Cardiovasc. Surg._ 143, 873–878. https://doi.org/10.1016/j.jtcvs.2011.09.070 (2012).

Article  CAS  PubMed  PubMed Central  Google Scholar  * Shah, S. H. _et al._ Baseline metabolomic profiles predict cardiovascular events in patients at risk for coronary artery disease.

_Am. Heart J._ 163, 844-850.e841. https://doi.org/10.1016/j.ahj.2012.02.005 (2012). Article  CAS  PubMed  Google Scholar  * Batch, B. C. _et al._ Branched chain amino acids are novel

biomarkers for discrimination of metabolic wellness. _Metabolism_ 62, 961–969. https://doi.org/10.1016/j.metabol.2013.01.007 (2013). Article  CAS  PubMed  PubMed Central  Google Scholar  *

Bhattacharya, S. _et al._ Validation of the association between a branched chain amino acid metabolite profile and extremes of coronary artery disease in patients referred for cardiac

catheterization. _Atherosclerosis_ 232, 191–196. https://doi.org/10.1016/j.atherosclerosis.2013.10.036 (2014). Article  CAS  PubMed  Google Scholar  * Shah, S. H. _et al._ Reclassification

of cardiovascular risk using integrated clinical and molecular biosignatures: design of and rationale for the measurement to understand the reclassification of disease of Cabarrus and

Kannapolis (MURDOCK) horizon 1 cardiovascular disease study. _Am. Heart J._ 160, 371-379.e372 (2010). Article  PubMed  Google Scholar  * Executive Summary of The Third Report of The National

Cholesterol Education Program (NCEP) Expert Panel on Detection, Evaluation, And Treatment of High Blood Cholesterol In Adults (Adult Treatment Panel III). _Jama_ 285, 2486–2497 (2001). *

Matthews, D. R. _et al._ Homeostasis model assessment: insulin resistance and beta-cell function from fasting plasma glucose and insulin concentrations in man. _Diabetologia_ 28, 412–419

(1985). Article  CAS  PubMed  Google Scholar  * Otvos, J. D. Measurement of lipoprotein subclass profiles by nuclear magnetic resonance spectroscopy. _Clin. Lab._ 48, 171–180 (2002). CAS 

PubMed  Google Scholar  * Shah, S. H. _et al._ Association of a peripheral blood metabolic profile with coronary artery disease and risk of subsequent cardiovascular events. _Circul.

Cardiovasc. Genet._ 3, 207–214. https://doi.org/10.1161/circgenetics.109.852814 (2010). Article  CAS  Google Scholar  * Kraus, W. E. _et al._ Metabolomic quantitative trait loci (mQTL)

mapping implicates the ubiquitin proteasome system in cardiovascular disease pathogenesis. _PLoS Genet_ 11, e1005553. https://doi.org/10.1371/journal.pgen.1005553 (2015). Article  CAS 

PubMed  PubMed Central  Google Scholar  * Duprez, D. A. _et al._ Comparison of the predictive value of GlycA and other biomarkers of inflammation for total death, incident cardiovascular

events, noncardiovascular and noncancer inflammatory-related events, and total cancer events. _Clin. Chem._ 62, 1020–1031. https://doi.org/10.1373/clinchem.2016.255828 (2016). Article  CAS 

PubMed  Google Scholar  * Lawler, P. R. _et al._ Circulating N-linked glycoprotein acetyls and longitudinal mortality risk. _Circ. Res._ 118, 1106–1115.

https://doi.org/10.1161/circresaha.115.308078 (2016). Article  CAS  PubMed  PubMed Central  Google Scholar  * Wang, H. _et al._ Sildenafil treatment in heart failure with preserved ejection

fraction: targeted metabolomic profiling in the RELAX trial. _JAMA Cardiol._ 2, 896–901. https://doi.org/10.1001/jamacardio.2017.1239 (2017). Article  PubMed  PubMed Central  Google Scholar

  * Rader, D. J. & Hovingh, G. K. HDL and cardiovascular disease. _Lancet (London, England)_ 384, 618–625. https://doi.org/10.1016/s0140-6736(14)61217-4 (2014). Article  CAS  Google

Scholar  * Otvos, J. D. _et al._ Low-density lipoprotein and high-density lipoprotein particle subclasses predict coronary events and are favorably changed by gemfibrozil therapy in the

Veterans Affairs High-Density Lipoprotein Intervention Trial. _Circulation_ 113, 1556–1563. https://doi.org/10.1161/circulationaha.105.565135 (2006). Article  CAS  PubMed  Google Scholar  *

Lu, Y. W. _et al._ Lipid paradox in patients with acute myocardial infarction: potential impact of malnutrition. _Clin. Nutr. (Edinburgh, Scotland)_

https://doi.org/10.1016/j.clnu.2018.10.008 (2018). Article  Google Scholar  * McGarrah, R. W., Crown, S. B., Zhang, G. F., Shah, S. H. & Newgard, C. B. Cardiovascular metabolomics.

_Circ. Res._ 122, 1238–1258. https://doi.org/10.1161/CIRCRESAHA.117.311002 (2018). Article  CAS  PubMed  PubMed Central  Google Scholar  Download references ACKNOWLEDGEMENTS We are grateful

to all participants in the CATHGEN study. FUNDING Funding for this study was supported by the National Institutes of Health Grant R01-HL127009, the American Heart Association pediatric

Strategically Focused Research Network Grant 17SFRN33670990, and National Heart, Lung and Blood Institute Grant 5P01-HL036587. AUTHOR INFORMATION Author notes * These authors contributed

equally: Stephen Ellison and Jawan W. Abdulrahim. AUTHORS AND AFFILIATIONS * Department of Anesthesiology, Duke University Medical Center, Durham, NC, USA Stephen Ellison * Duke Molecular

Physiology Institute, Duke University School of Medicine, 300 North Duke St, Durham, NC, 27701, USA Jawan W. Abdulrahim, Lydia Coulter Kwee, Nathan A. Bihlmeyer, Robert McGarrah, James R.

Bain, William E. Kraus & Svati H. Shah * Division of Cardiology, Department of Medicine, Duke University School of Medicine, Durham, NC, USA Neha Pagidipati, Robert McGarrah, William E.

Kraus & Svati H. Shah Authors * Stephen Ellison View author publications You can also search for this author inPubMed Google Scholar * Jawan W. Abdulrahim View author publications You

can also search for this author inPubMed Google Scholar * Lydia Coulter Kwee View author publications You can also search for this author inPubMed Google Scholar * Nathan A. Bihlmeyer View

author publications You can also search for this author inPubMed Google Scholar * Neha Pagidipati View author publications You can also search for this author inPubMed Google Scholar *

Robert McGarrah View author publications You can also search for this author inPubMed Google Scholar * James R. Bain View author publications You can also search for this author inPubMed 

Google Scholar * William E. Kraus View author publications You can also search for this author inPubMed Google Scholar * Svati H. Shah View author publications You can also search for this

author inPubMed Google Scholar CONTRIBUTIONS S.E. conceived the presented idea, was involved with developing the data analysis plan and interpreting the data, and was primarily responsible

for writing the manuscript. J.W.A. contributed to manuscript revision, data interpretation and wrote part of the manuscript. L.C.K. wrote part of the initial draft, edited the manuscript,

and curated and analyzed the data. N.A.B. contributed to interpretation of data, N.P. critically revised and edited the manuscript. R.M. was involved in data interpretation, writing and

editing. J.R.B. helped to plan and execute the study’s metabolomics assays and critically reviewed the manuscript. W.E.K. conceived of parent project, recruited funding, obtained samples and

data, review of analyses, editing of draft manuscript. S.H.S. secured funding, contributed to study design, supervised the project, and critically revised and edited the final manuscript.

CORRESPONDING AUTHOR Correspondence to Svati H. Shah. ETHICS DECLARATIONS COMPETING INTERESTS Dr. Pagidipati’s work has been funded by: Amarin Pharmaceutical Company; Amgen, Inc.;

AstraZeneca; Baseline Study LLC; Boehringer Ingleheim; Duke Clinical Research Institute; Eli Lilly & Company; Novo Nordisk Pharmaceutical Company; Regeneron Pharmaceuticals, Inc.;

Sanofi-S.A.; Verily Sciences Research Company. She reports consulting fees from AstraZeneca; Boehringer Ingleheim; Esperion Therapeutics. Dr. Ellison, Dr. Kwee, Dr. Bihlmeyer, Dr.

Abdulrahim, Dr. McGarrah, and Dr. Bain declare no potential competing interests. Dr. Shah and Dr. Kraus have an unlicensed patent on a related research finding and Dr. Shah receives

sponsored research funding from Baseline Study LLC. ADDITIONAL INFORMATION PUBLISHER'S NOTE Springer Nature remains neutral with regard to jurisdictional claims in published maps and

institutional affiliations. SUPPLEMENTARY INFORMATION SUPPLEMENTARY INFORMATION. RIGHTS AND PERMISSIONS OPEN ACCESS This article is licensed under a Creative Commons Attribution 4.0

International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the

source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's

Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not

permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit

http://creativecommons.org/licenses/by/4.0/. Reprints and permissions ABOUT THIS ARTICLE CITE THIS ARTICLE Ellison, S., Abdulrahim, J.W., Kwee, L.C. _et al._ Novel plasma biomarkers improve

discrimination of metabolic health independent of weight. _Sci Rep_ 10, 21365 (2020). https://doi.org/10.1038/s41598-020-78478-w Download citation * Received: 15 April 2020 * Accepted: 18

November 2020 * Published: 07 December 2020 * DOI: https://doi.org/10.1038/s41598-020-78478-w SHARE THIS ARTICLE Anyone you share the following link with will be able to read this content:

Get shareable link Sorry, a shareable link is not currently available for this article. Copy to clipboard Provided by the Springer Nature SharedIt content-sharing initiative