Introduction

During the last decade, translational research has highlighted intestinal microbiota as possible mediators between dietary habits and both the development and progression of coronary artery disease (CAD) [1,2,3,4,5]. The consumption of red-meat and dairy products rich in choline, betaine, carnitine, and trimethyllysine leads to the production of trimethylamine (TMA) by certain intestinal microbiota [1, 6, 7]. In a second step, TMA is oxidized to trimethylamine N-oxide (TMAO) by flavin monooxygenase in the liver [1, 8,9,10,11,12]. TMAO seems to induce systemic inflammation at least in part by the activation of the NF-κB pathway and the increased expression of pro-inflammatory cytokines including TNF-α and IL-1β [6, 13,14,15,16]. TMAO was indicated in multiple studies to accelerate atherosclerosis and enhance platelet reactivity as well as thrombosis potential [3, 17,18,19]. In addition, recent studies have documented an association between plasma TMAO concentration and the risk of death, myocardial infarction (MI), and stroke in patients with either stable CAD or acute coronary syndromes [1, 2, 4, 10, 20,21,22].

Beyond possible therapeutic opportunities, this insight suggests that TMAO and/or its precursors might have prognostic and/or diagnostic utility in the non-invasive detection of CAD, particularly the more aggressive CAD phenotype leading to myocardial ischemia during everyday activities (functionally relevant CAD, fCAD). Even more importantly, TMAO may provide prognostic utility in identifying those at incident risk for hard clinical endpoints including death and MI. For this indication, TMAO and/or its precursors may help physicians in the selection of patients for cardiac work-up including invasive or non-invasive coronary imaging [23, 24]. Given the high number of patients with very low pre-test probability for fCAD referred for sophisticated cardiac imaging including myocardial perfusion scanning, and with cardiac imaging causing annual costs of more than $500 million in the United States alone, biomarker guidance may have substantial medical and economic value [23, 25]. Therefore, the aim of this study was to prospectively assess the clinical prognostic utility performance and the diagnostic accuracy of circulating TMAO and its precursors on all-cause death, cardiovascular death, AMI and the composite endpoint of CV death and AMI, in patients with suspected fCAD.

Methods

Study design and oversight

This analysis is part of a large ongoing prospective diagnostic study (NCT01838148, clinicaltrials.gov) designed to advance the early detection of fCAD [26, 27]. The local ethics committee approved the study, which was carried out according to the principles of the Declaration of Helsinki. All patients provided written informed consent. The authors designed the study, gathered, analyzed and vouched for the data and analysis, wrote the paper, and made the decision to submit it for publication. Reported data follow STARD guidelines for studies of diagnostic accuracy [28].

Patient population

Patients were recruited from 2010 to 2016 at the University Hospital of Basel, Switzerland. Enrolled patients were suspected to have fCAD and were referred for rest/stress myocardial perfusion single-photon emission tomography/computer tomography (MPI-SPECT). MPI-SPECT/CT was the preferred cardiac imaging technique in patients with a wide range of pre-test probabilities for fCAD during that time. Patients requiring chronic dialysis were excluded.

Quantified clinical assessment

The likelihood for the presence of fCAD was quantified by the integrated clinical judgment of the treating cardiologist using a visual analogue scale (VAS) ranging from 0 to 100% twice: once before stress testing integrating all medical information available at that time, such as age, sex, cardiovascular risk factors, previous cardiac history, symptoms and baseline ECG data; second, after stress testing, integrating symptoms the patient experienced during exercise/stress, the workload achieved, and ECG changes recorded during exercise/stress. The cardiologist was blinded to both biomarker measurements and MPI-SPECT images at the time of assessment.

Blood sampling and laboratory methods

Venous non-fasting blood samples for determination of TMAO, choline, betaine, carnitine, and high-sensitivity cardiac troponin (hs-cTn) T [23, 26, 29], an established cardiovascular biomarker associated with fCAD, were obtained at rest, before stress testing. After centrifugation, samples were frozen at − 80 °C until assayed in a blinded fashion in a dedicated core laboratory. Serum TMAO, choline, betaine and carnitine were quantified using stable isotope dilution LC/MS/MS analyses as previously described using a Shimadzu Nexera Ultra High Performance Liquid Chromatograph (UHPLC) system interfaced with Shimadzu 8050 Triple Quadrupole Mass Spectrometer [1, 3, 18]. Hs-cTnT was measured with the Elecsys System on the Modular Analytics E170 or the Cobas e 602 (Roche Diagnostics, Rotkreuz, Switzerland). Limit of blank and limit of detection of this assay are 3.0 ng/L and 5.0 ng/L, respectively. The upper reference limit (URL) of a healthy reference population was 14 ng/L with an imprecision corresponding to 10% coefficient of variation (CV) at 13 ng/L [29]. Creatinine was measured by either the University Hospital of Basel central laboratory, Risch laboratories and Rothen laboratories. eGFR (estimated glomerular filtration rate) was calculated using the CKD-EPI (Chronic Kidney Disease Epidemiology Collaboration) equation provided within the transplantr-package in R. Cystatin-C was measured by SomaLogic using the SomaScan® assay, which uses single-stranded DNA-based protein affinity reagents that are modified to mimic amino acid chains, enhancing protein-nucleic acid interaction [30]. These so-called SOMAmer reagents are selected against proteins in their native folded state, which after binding and other processing steps can then be quantified by DNA quantification techniques, providing a readout in relative fluorescent units (RFU) directly proportional to the amount of the protein [30].

Adjudication of the presence of fCAD

Expert interpretation of MPI-SPECT/CT images combined with information obtained from invasive coronary angiography and whenever available fractional flow reserve measurements were used in the adjudication of fCAD.

All patients underwent a routine rest/stress dual isotope (201Tl for rest, 99mTc sestamibi for stress) or single isotope (99mTc sestamibi for stress and rest) MPI-SPECT protocol as described previously [27, 31, 32]. MPI-SPECT images were scored semi-quantitatively using a 17-segment model with a 5-point scale (0 = normal, 1 = mildly reduced tracer uptake, 2 = moderately reduced tracer uptake, 3 = severely reduced tracer uptake and 4 = no uptake). The summed rest score (SRS) and summed stress score (SSS) were calculated based on the 17 segments in the rest and stress images. The difference of SRS and SSS yielded the summed difference score (SDS), whereby an SDS of at least two or a positive transient ischemic dilation ratio (TID) was considered as inducible myocardial ischemia. Two readers derived SSS and SRS by visual assessment and compared with the software (QGS) result. Differences in the visual assessment by the two readers were resolved by consensus. In case of equivocal findings from MPI-SPECT and coronary angiography, two independent cardiologists (one interventional cardiologist, one general cardiologist) that were blinded to biomarker results reviewed the case. A positive perfusion scan was overruled when coronary angiography showed normal coronary arteries and a negative perfusion scan was overruled if coronary angiography (within three months) revealed a high-grade coronary lesion (> 75% or fractional flow reserve (FFR) lower than 0.80) [23, 26, 29].

Adjudication of major adverse cardiac events

The prognostic endpoints were all-cause death, cardiovascular death and MI during long-term follow-up. Patients were contacted by telephone or standardized follow-up letter 1 year, 2 years and 5 years after enrolment with ongoing follow-up. In case of an event, further information was obtained from hospital records, general practitioner/cardiologists records, or the national death registry blinded to biomarker concentrations.

Statistical analysis

Normality testing was done using visual assessments (histograms and QQ-plots) and as Shapiro–Wilk test. Continuous variables are presented as median and respective interquartile range (IQR) and categorical variables are presented as frequencies and respective percentages. The Agresti-Coull method was used to calculate confidence intervals of proportions. Measures of central tendency of the biomarkers were compared by Mann–Whitney U test. Baseline characteristics were compared by Kruskal test for continuous data and Fisher test for categorical data. Diagnostic accuracy of the clinical assessment, TMAO, betaine, choline, carnitine, and their combination for fCAD was quantified by the area under the receiver operating curve (AUC) and compared with the method described by DeLong et al. [33]. Based on previous findings showing an association between TMAO concentration and prevalent CAD [6, 11, 34], subgroup analysis was performed stratified according to the presence or absence of previously known CAD. We hypothesized that TMAO concentration would have highest diagnostic accuracy in patients without known CAD. Using logistic regression with fCAD as outcome, log2-transformed concentrations of TMAO or one of its precursors (betaine, choline, carnitine) were adjusted for patient characteristics, risk-factors and treatment, and combined with the quantified clinical assessment of the treating physician before and after stress testing. Multicollinearity was checked by the variance inflation factor (VIF) using the vif-function of the car package in R-Studio.

Kaplan–Meier curves for concentrations of TMAO, betaine, choline and carnitine below and above the median as well as stratified by the upper reference limit of the established biomarker hs-cTnT (URL 14 ng/L), were constructed and compared by log rank testing. Quantification of the predictive accuracy of TMAO and precursors was done by time dependent AUCs whilst accounting for censoring (timeROC package) [35]. Cox regression analysis was used to evaluate whether the analyzed biomarkers were independent predictors of patient characteristics and risk factors. Assumptions for Cox regression were tested using Schoenfeld residuals. Statistical analyses were performed with R version 4.0.2. All hypothesis testing was two-tailed, and a p value < 0.05 was considered statistically significant.

Results

Characteristics of patients

Overall, 1726 consecutive patients with suspected functionally relevant coronary artery disease (fCAD) were included in this analysis (Supplemental Fig. 1), 478 patients (28%) were adjudicated to have fCAD. A total of 421 (24%) patients underwent coronary angiography with 284 patients within the 3 months after enrolment. Anamnesis revealed 764 (44%) to have known history of CAD. During follow-up, 88 patients experienced an incident MI, 223 patients died overall, 115 of which died due to cardiovascular reasons (CV death). Characteristics of patients stratified by fCAD are shown in Table 1. Stratification of patients by experience of either an AMI, all-cause or CV death, is shown in Supplemental Table 1. Patients with fCAD tended to be older and a higher proportion were male. A significant proportion of patients had diabetes, hypertension, and a history of cardiovascular disease. To portrait renal state of the patients, baseline eGFR and cystatin-C were added to the baseline tables. As expected, eGFR was significantly lower in patients with fCAD and cystatin-C significantly higher in patients with fCAD.

Table 1 Patient baseline characteristics

Circulating TMAO, betaine, choline and carnitine concentrations

Median TMAO concentrations were significantly higher in patients with fCAD than in those without fCAD (5.33 IQR [3.55, 8.80] µM vs. 4.66 IQR [3.08, 7.24] µM, p < 0.001, Fig. 1A). Stratification by a history CAD revealed that this phenomenon was exclusively present in patients without known CAD (5.36 IQR [3.71, 8.57] µM vs 4.38 IQR [2.94, 6.83] µM, p < 0.001, Fig. 1B). Similar findings emerged for the three precursors betaine, choline and carnitine (Fig. 1C–H).

Fig. 1
figure 1

Levels of TMAO, betaine, choline and carnitine compared between patients adjudicated to have fCAD (A, C, E, G) and compared between patients adjudicated to have fCAD within subgroups of patients with and without history of CAD (B, D, F, H)

In patients with normal eGFR (n = 689), TMAO and the three precursors were significantly higher in patients adjudicated to have fCAD (Supplemental Fig. 2A–D). The baseline characteristics of this subgroup are presented in Supplemental Table 2.

TMAO concentrations were significantly, albeit weakly, correlated with hs-cTnT (Spearman’s rho 0.32, p < 0.001), NT-proBNP (Spearman’s rho 0.22, p < 0.001) and age (Spearman’s rho 0.25, p < 0.001). These correlations were weaker for betaine, choline and carnitine (Supplemental Fig. 3).

Diagnostic performance for the detection of fCAD

In the overall cohort, diagnostic accuracy for detection of fCAD was quantified and showed modest but statistically significant value for TMAO with a ROC AUC of (0.56, 95% CI 0.53–0.59, p < 0.001) and its precursors (betaine 0.55, 95% CI 0.52–0.58, p = 0.002; choline 0.54, 95% CI 0.51–0.57, p = 0.007; carnitine: 0.55, 95% CI 0.51–0.58; p = 0.001; Supplemental Fig. 4). Similar findings emerged in the subgroup of patients without known CAD (TMAO 0.59, 95% CI [0.54–0.63], p < 0.001; betaine 0.54, 95%-CI [0.50–0.58], p = 0.039; choline: 0.56, 95% CI [0.51–0.60], p = 0.006; carnitine 0.57, 95% CI [0.52–0.61], p = 0.002; Supplemental Fig. 4B). TMAO and its precursors did not significantly increase the AUC provided by the quantitative clinical assessment by the treating physician before (AUC 0.61, 95% CI [0.58–0.64]) and after (AUC 0.65, 95% CI [0.62–0.68]) stress testing (p > 0.05 for all comparisons of VAS + TMAO (or precursors) versus VAS alone, Supplemental Fig. 4C + 4D). TMAO (OR 1.19, 95% CI [1.08–1.31], p < 0.001) and its precursors were significant predictors of fCAD in the univariable model, with TMAO remaining a significant predictor even after adjusting for age, sex and CAD history (OR 1.12, 95% CI [1.01–1.24], p = 0.036, Supplemental Table 3). When adjusting for renal function (i.e., taking Cystatin-C into the model) or when adjusting for further for pre-defined patient characteristics, cardiovascular risk factors and medical history (model 2), neither TMAO nor its precursors remained significant predictors. In patients with normal eGFR, TMAO and the three precursors were only significant predictors in the univariable model but were no longer significant after adjusting further (Supplemental Table 3).

Prognostic performance for incident major adverse cardiac events

The median follow-up time was 1827.5 days (IQR [756, 1908]). The cumulative event incidence for 5-years death was 12.9% (223 events), for 5-years CV death 6.7% (115 events) and for 5-years MI 5.1% (88 events), and the composite 5-years MACE (MI and cardiovascular death) was 11.8% (203 events).

In Kaplan–Meier analyses, TMAO concentrations above the median had substantially higher all-cause, CV mortality and mace events versus patients below the median (P < 0.001, each; Supplemental Fig. 5 and 6). Similar results emerged for choline and carnitine.

Kaplan–Meier analyses revealed that stratification according to hs-TnT and TMAO concentrations below/above the URL provided incremental prognostic value with the highest all-cause and CV mortality observed in patients with both hs-cTnT and TMAO above the URL (Fig. 2A).

Fig. 2
figure 2

Kaplan Meier survival analysis of TMAO combined with low and high levels hs-Troponin (below/above URL of 14 ng/L; panel A) and TMAO combined with low and normal ranges of eGFR (> 60 mL/min/1.73m2 considered normal eGFR; panel B—in patients with available eGFR data (n = 921))

Subgroup analysis (n = 919) of marker-combination with low and normal eGFR levels, Kaplan–Meier survival analyses revealed significantly highest cumulative all-cause death events in patients with low eGFR combined with a TMAO above the median (Fig. 2B).

Time-dependent ROC curve analysis showed a consistent and moderate-to-good discriminative performance of TMAO and choline for 5-years all-cause death (e.g. TMAO at 2-years AUC: 0.67, 95% CI [0.61–0.73]; Choline at 2-years AUC: 0.64, 95%-CI [0.57–0.70]) and 5-years CV death (TMAO at 2-year, AUC: 0.73, 95% CI [0.65–0.81]; Choline at 2-year AUC: 0.67, 95% CI [0.59, 0.76]), and modest accuracy for the other markers and endpoints (Fig. 3). Furthermore, Cox-regression analyses on continuous log-transformed markers are shown in Table 2. TMAO, Carnitine and Choline remained independent predictors in the full model without adjustment of renal function (model 2). When adjusting further for renal function (model 2 + cystatin-c) both TMAO and carnitine remained independent predictors for 5-years CV death. In the cox-regression analysis performed on median stratified marker (Supplemental Table 4), only TMAO remained a significant predictor of all-cause death and CV death, even when adjusted fully and taking renal function into account (model 2 + cystatin-c).

Fig. 3
figure 3

Time-dependent AUC of ROC of the four markers and over 5-years for all-cause death (223 events), cardiovascular death (115 events), AMI (88 events) and mace (203 events). Each step represents a change in AUC due to an event or censoring. AUC in the first days (ca. 180) vary due to small number of events in early stages

Table 2 Hazard ratios of the univariate and adjusted Cox regression models with continuous log2-transformed markers for the outcomes all-cause death, cardiovascular death and acute myocardial infarction (CI—confidence interval; HR—Hazard Ratio)

Discussion

The accurate, non-invasive and inexpensive detection of fCAD is a major unmet clinical need. So too is the accurate and minimally invasive ability to identify those at increased risk for incident adverse events. The present study suggests that TMAO has clinical prognostic value for the identification of subjects at risk for incident adverse cardiac events, with modest, albeit significant, capacity to predict risk of fCAD. Recent translational insights suggest TMAO as a potential mediator of atherosclerosis [1, 2, 10, 36, 37]. Even more so, recent mechanistic studies suggest TMAO contributes to subject vulnerability for adverse events like thrombosis (MI and stroke) through both heightened platelet responsiveness, and enhanced artery wall inflammatory signaling [38, 39]. We, therefore, tested the hypothesis that TMAO, and its precursors betaine, choline and carnitine, may provide diagnostic and/or prognostic value in a large prospective study including 1726 patients with suspected fCAD referred for cardiac work-up using MPI-SPECT. We report five major findings.

First, TMAO, betaine, choline and carnitine concentrations were significantly higher in patients with fCAD versus those without fCAD. This is in line with prior pilot studies reporting that TMAO concentrations were positively associated with atherosclerotic burden in stable CAD patients [2, 40]. Interestingly, the difference observed was largely restricted to patients without known CAD history, the cohort in whom a diagnostic test is, arguably, more relevant. Also, TMAO was demonstrated to be associated with plaque instability [41]—therefore, TMAO possibly was only significantly elevated in that subgroup as the marker may play a role in soft plaques, rather than in calcified plaques. Second, TMAO and precursor marker levels were lower in patients with normal eGFR (> 60 mL/min per 1.73 m2) and patients adjudicated to have fCAD had significantly higher marker levels compared to those without. Notably, high TMAO in combination with low-eGFR were at significantly higher risk of all-cause death compared to low-TMAO combined with low-eGFR (n = 921). This is in line with previous findings demonstrating clearance of TMAO by the kidney and showing its association with reduced renal function [42,43,44,45]. Baseline eGFR additionally was significantly lower in patients with fCAD, supporting findings of lower eGFR being associated with increased risk of coronary artery disease [46,47,48]. Taking renal function, therefore, into account in the adjusted cox regression models, TMAO, choline and carnitine remained significant predictors for CV death. Third, while the AUC for fCAD for TMAO and its precursors was relatively low, it was statistically significant. Plasma levels of TMAO and its precursors (choline, betaine and carnitine) did not provide incremental diagnostic value for prediction of fCAD to clinical judgment of the cardiologist. Nonetheless, all markers were significant risk factors for fCAD in the univariate model, with TMAO remaining significant in models adjusted for age, sex and CAD history. Fourth, in contrast, concentrations of TMAO, carnitine and to a lesser extent choline were significant predictors of incident adverse events including either MI, CV-death, or all cause death. This finding is in concordance with previous reports where TMAO has been shown to be associated with all-cause and CV-death [1, 2, 4, 10, 20,21,22]. Moreover, significantly more adverse events (all-cause or CV death, or AMI) occurred in patients with high TMAO in combination with a hs-cTnT above the URL, a quantitative marker of cardiomyocyte injury strongly associated with future CV-death [23, 49]. Fifth, time-dependent ROC curve analysis showed consistent moderate-to-good prognostic accuracy of TMAO and its precursors for all-cause and CV-death during follow-up. The prognostic value for TMAO persisted also in fully adjusted models, highlighting possible implications for routine clinical practice given similar observations in recent studies [4, 20, 27, 36]. Overall, the prognostic accuracy for all-cause and CV-death was highest and most consistent for TMAO, and slightly lower and less consistent for its precursors and other outcomes. The same has been found in a stable CAD cohort where elevated TMAO likewise was associated with higher long-term mortality risk. [11]

Several therapeutic strategies could result in a reduction of TMAO concentrations: reduction in intake of red meat, rich in TMA precursors; use of targeted treatment of selected microbiota, e.g. by application of specific oral antibiotics to eradicate the intestinal microbiota responsible for TMA production; and inhibition of flavin monooxygenase in the liver [2, 3, 34, 50]. Clinical intervention studies could further elucidate whether TMAO is a true causal mediator or a bystander in atherothrombosis.

This study has important strengths, including central adjudication of fCAD, adjudication of hard outcomes, and prospective examination of the relationship of TMAO with respect to fCAD and incident MI, CVD death, and all-cause mortality over a long-term follow-up of 5 years. In addition, the study examines the relationship of multiple nutrient TMAO precursors betaine, choline and carnitine in a large prospective study of patients in need of active decision making. Also, several limitations should be considered when interpreting the findings of this study: first, these data were generated in a single-center study. While single-center studies are by definition prone to selection bias, generalizability of these findings seems high as standardized patient consenting, clinical work-up with MPI-SPECT/CT as the initial imaging modality was performed in patients with a wide range of pre-test probability for fCAD, and longer term outcome ascertainment was complete. Second, despite using a very stringent methodology for the adjudication of fCAD, we might still have misclassified a small number of patients, which would result in an underestimation of the true accuracy of TMAO and its precursors for fCAD. This is overcome by the examination of the relationship between TMAO and time dependent incident hard adverse events. Third, the cohort consisted of a predominantly Caucasian population, not allowing the analysis of other ethnicities due to their under-representation. Finally, we cannot comment on the dietary association of TMAO or precursor level, as this was not recorded in the frame of the study, nor the possible role of TMAO and its precursors in patients with terminal kidney failure on chronic dialysis, as these subjects were excluded.

In conclusion, TMAO and its precursors choline and carnitine showed potential prognostic value for short- and long-term risk-stratification for hard clinical events including CV-death, albeit with limited diagnostic utility in patients with suspected fCAD.