Bibliografía

Buenos Aires 01 de Febrero del 2025

Prognostic Value of Cardiovascular Biomarkers in the Population

 

 

Prognostic Value of Cardiovascular Biomarkers in the Population

 

JohannesT. Neumann, MD,PhDRaphael Twerenbold, MDJessica Weimann, MScet al   

 
JAMA. 2024;331(22):1898-1909. doi:10.1001/jama.2024.5596

              

 

Introduction

Early identification of individuals in the general population at high risk for atherosclerotic cardiovascular disease shapes primary preventive strategies to reduce the risk of developing atherosclerotic cardiovascular disease.1,2 Risk scores based on traditional risk factors for atherosclerotic cardiovascular disease (eg, the European Society of Cardiology Systematic Coronary Risk Evaluation 2 [SCORE2], the American Heart Association/American College of Cardiology Pooled Cohort Equations, and the American Heart Association Predicting Risk of Cardiovascular Disease Events [PREVENT] equations) are widely available to estimate an individual’s risk for future cardiovascular events.1-4
Cardiovascular biomarkers, such as cardiac troponin, natriuretic peptides, and C-reactive protein (CRP), are established in clinical care. Using newer, high-sensitivity cardiac troponin assays, concentrations became measurable in the general population, opening up the prospects for a broader application of this biomarker.5 Several studies have reported (1) strong associations of these biomarkers with incident atherosclerotic cardiovascular disease events in individuals with known atherosclerotic cardiovascular disease, but also, and most importantly, in apparently healthy individuals and (2) an improvement in risk stratification when these biomarkers were added to established risk prediction models.5-12
Notwithstanding the achievements of earlier studies, the actual application of routinely available cardiovascular biomarkers for risk stratification in primary prevention has not become routine clinical practice. In addition, it remains uncertain which of the established biomarkers might be best suited to predict each outcome and how such associations are influenced by age.
Therefore, this study brings together the largest multinational individual-level dataset, to date, to investigate the comparative predictive value of cardiovascular biomarkers for incident atherosclerotic cardiovascular disease events in the general population and to elucidate their differential effects according to age.

Methods

* Study Cohorts
Details on the selection of study cohorts appear in the eMethods in Supplement 1. This individual-level analysis included data from 28 multinational population-based cohorts (eTable 1 in Supplement 1). Eligible cohorts were those that included (1) individuals from the general population, in which most participants were apparently healthy (ie, had not had any major atherothrombotic cardiovascular events); (2) individuals who had at least 1 measurement of high-sensitivity cardiac troponin I, high-sensitivity cardiac troponin T, B-type natriuretic peptide (BNP), N-terminal pro-BNP (NT-proBNP), or high-sensitivity CRP; and (3) individuals with follow-up for at least 2 years. Data from all cohorts were collected and harmonized in a database. Individuals with a history of atherosclerotic cardiovascular disease events or heart failure were excluded from the analyses (Figure 1).
* Study Outcomes
The primary outcome was incident atherosclerotic cardiovascular disease, which included all fatal and nonfatal events. Incident atherosclerotic cardiovascular disease was defined by the first possible or definite coronary heart disease event, possible or definite ischemic stroke event, coronary revascularization, coronary heart disease death, ischemic stroke death, or unclassifiable death.13 The secondary outcomes were all-cause mortality, incident heart failure, incident ischemic stroke, and incident myocardial infarction. Additional information about the outcomes investigated appear in eTable 2 in Supplement 1.
* Biomarkers
For all cohorts reported by the Biomarker for Cardiovascular Risk Assessment across Europe (BiomarCaRE) consortium, serum high-sensitivity cardiac troponin I concentration was determined in the BiomarCaRE core laboratory in Hamburg, Germany, using a highly sensitive cardiac troponin I immunoassay for stored samples (Architect i2000SR; Abbott Diagnostics). The limit of detection for the immunoassay was 1.9 ng/L (range, 0-50 000 ng/L) and the assay had a coefficient of variation of 10% at a concentration of 5.2 ng/L. Measurement of NT-proBNP concentration was performed using an electrochemiluminescence immunoassay (ELECSYS 2010 and Cobas e411; Roche Diagnostics); the analytic range is 0 ng/L to 35 000 ng/L.
Measurement of high-sensitivity CRP concentration was performed using the Vario immunoassay and the Architect c8000 system (Abbott Diagnostics). For the cohorts not part of the BiomarCaRE consortium, measurements of high-sensitivity cardiac troponin I, high-sensitivity cardiac troponin T, NT-proBNP, BNP, and high-sensitivity CRP were performed as part of the local cohort–specific procedures (details on the specific assays used in each cohort appear in eTable 3 in Supplement 1).
* Statistical Analyses
Associations Between Cardiac-Specific Biomarkers and Study Outcomes
To examine the unadjusted association of the biomarkers and the primary outcome, cumulative incidence curves were computed according to biomarker quintiles. Death from causes unrelated to atherosclerotic cardiovascular disease was treated as a competing event. The curves were estimated using the Aalen-Johansen estimator. The quintiles were computed using linear quantile mixed models.14,15
Fine and Gray subdistribution hazard models were calculated. Death from causes unrelated to atherosclerotic cardiovascular disease or the secondary outcomes were treated as a competing event, respectively.16 Cox proportional hazards regression models were used for all-cause mortality. The biomarkers were used as continuous variables after applying loge transformation with hazard ratios (HRs) or subdistribution HRs computed additionally per 1-SD change to allow for comparisons of the effect size among different biomarkers. The models with high-sensitivity cardiac troponin concentration as a continuous variable were augmented with a binary variable indicating if the measured value was above or below the limit of detection.17
The regression models for all the study outcomes were adjusted for sex and cohort as stratification variables. The models were also adjusted for age, total cholesterol, high-density lipoprotein (HDL) cholesterol, current smoking, prevalent diabetes, systolic blood pressure, and self-reported use of antihypertensive drugs. For the outcomes of all-cause mortality and heart failure, the models were additionally adjusted for body mass index (calculated as weight in kilograms divided by height in meters squared). In separate analyses, high-sensitivity cardiac troponin I, NT-proBNP, and high-sensitivity CRP were combined in a multivariable model because these biomarkers represent different pathophysiological pathways and the data were readily available in the cohorts. The time-to-event models were extended by modeling the biomarkers using penalized cubic splines.
Added Predictive Value
The C statistic and net reclassification improvement (NRI) were used to quantify the added predictive value of each biomarker beyond that from the adjusted model described above. In the presence of competing risks, the calculations of the C statistic and NRI were adapted.18,19 Internal-external cross-validation was applied to control for the overoptimism of calculating performance measures on the same dataset from which the models were computed.20 Namely, each study was in turn left out and the Cox model or the Fine and Gray model was estimated in the remaining studies. Next, the models were used to estimate the event probabilities in the excluded study. The category-based NRI and the continuous NRI were calculated. The 95% CIs for the C statistic and NRI were computed by bootstrapping 500 times the internal-external cross-validation.
All statistical analyses were performed using R version 4.2.2 (R Foundation for Statistical Computing).21 Additional information about the statistical analyses appears in the eMethods in Supplement 1.

Results

* Study Population
There were 28 general population–based cohorts included (from 12 countries and 4 continents) with data on 164 054 individuals (median age, 53.1 years [IQR, 42.7-62.9 years]; 85 972 [52.4%] were women; and 9977 [6.1%] had diabetes) (Figure 1Table, and eTable 4 in Supplement 1). Of the 162 947 individuals with data for the hypertension variable, 67 719 (41.6%) had hypertension. Of the 162 139 individuals with data for the smoking variable, 40 226 (24.8%) smoked daily. The median 10-year atherosclerotic cardiovascular disease risk SCORE2 was 4.1% (IQR, 1.7%-8.6%) and the corresponding 10-year risk using the Pooled Cohort Risk Equation was 4.9% (IQR, 1.4%-13.1%).
The median biomarker concentrations were 2.5 ng/L (IQR, 1.9-4.1 ng/L) for high-sensitivity cardiac troponin I, 3.1 ng/L (IQR, 3.0-6.0 ng/L) for high-sensitivity cardiac troponin T, 43.8 ng/L (IQR, 20.6-86.2 ng/L) for NT-proBNP, 14.9 ng/L (IQR, 7.9-28.6 ng/L) for BNP, and 1.4 mg/L (IQR, 0.7-3.2 mg/L) for high-sensitivity CRP (Table and eFigure 1 in Supplement 1). Except for NT-proBNP, there was a linear relationship between age and the median biomarker concentrations (eFigure 2 in Supplement 1).
During a median follow-up of 11.8 years (IQR, 6.2-18.0 years; maximum follow-up, 28.2 years), there were 17 211 incident atherosclerotic cardiovascular disease events, 25 346 deaths from any cause, 6766 cases of heart failure, 4794 incident cases of incident ischemic stroke, and 8024 incident cases of myocardial infarction (eTable 5 in Supplement 1).
* Primary Outcome: Association of Biomarkers With Incident Atherosclerotic Cardiovascular Disease Events
After adjusting for sex and cohort and the conventional risk factors of age, total cholesterol, HDL cholesterol, smoking status, diabetes, systolic blood pressure, and self-reported use of antihypertensive drugs, the biomarker concentrations were associated with incident atherosclerotic cardiovascular disease events (subdistribution HR per 1-SD change, 1.13 [95% CI, 1.11-1.16] for high-sensitivity cardiac troponin I; 1.18 [95% CI, 1.12-1.23] for high-sensitivity cardiac troponin T; 1.21 [95% CI, 1.18-1.24] for NT-proBNP; 1.14 [95% CI, 1.08-1.22] for BNP; and 1.14 [95% CI, 1.12-1.16] for high-sensitivity CRP; Figure 2 and Figure 3). For all 5 biomarkers, there were more events per 1000 person-years in individuals with biomarker concentrations above the median compared with those with biomarker concentrations below the median (Figure 3). Additional data appear in eTable 6 in Supplement
In separate analyses, high-sensitivity cardiac troponin I, NT-proBNP, and high-sensitivity CRP were included in the same model and proved to be predictors of atherosclerotic cardiovascular disease events (adjusted subdistribution HR, 1.07 [95% CI, 1.04-1.10] for high-sensitivity cardiac troponin I; 1.19 [95% CI, 1.15-1.23] for NT-proBNP; and 1.14 [95% CI, 1.10-1.17] for high-sensitivity CRP). When stratified by biomarker quintiles, the cumulative atherosclerotic cardiovascular disease incidence gradually increased with increasing biomarker concentrations (eFigure 3 in Supplement 1).
The addition of the biomarkers to the base model, which included only conventional risk factors, was associated with an increase in the C statistics for atherosclerotic cardiovascular disease events after 1 year, 5 years, and 10 years (Figure 4 and eTable 7 in Supplement 1). The strongest increase was observed when high-sensitivity cardiac troponin I, NT-proBNP, and high-sensitivity CRP were combined in 1 model. In the reclassification analyses, the categorical NRI for the combination of high-sensitivity cardiac troponin I, NT-proBNP, and high-sensitivity CRP was 0.044 (95% CI, 0.023-0.069) (eTable 8 in Supplement 1).
The continuous NRI was 0.241 (95% CI, 0.193-0.309) for high-sensitivity cardiac troponin I, 0.201 (95% CI, 0.008-0.364) for high-sensitivity cardiac troponin T, 0.06 (95% CI, 0.016-0.093) for NT-proBNP, 0.077 (95% CI, 0.000-0.144) for BNP, 0.192 (95% CI, 0.162-0.222) for high-sensitivity CRP, and 0.23 (95% CI, 0.162-0.283) for the combination of high-sensitivity cardiac troponin I, NT-proBNP, and high-sensitivity CRP (eTable 9 in Supplement 1).
* Association of Biomarkers With Secondary Outcomes
All biomarkers were associated with all-cause mortality, incident heart failure, incident ischemic stroke, and incident myocardial infarction (Figure 3 and eTables 10-13 in Supplement 1). The associations of biomarkers with all-cause mortality and incident heart failure were larger than those for atherosclerotic cardiovascular disease. The highest subdistribution HRs were observed for incident heart failure with high-sensitivity cardiac troponin T (HR, 1.44; 95% CI, 1.38-1.51), NT-proBNP (HR, 1.62; 95% CI, 1.56-1.68), and BNP (HR, 1.59; 95% CI, 1.43-1.77). The addition of the biomarkers also improved the C statistics (Figure 4) and the appropriate classification of risk (eTables 8-9 in Supplement 1) when added to the base model for the secondary outcomes. The largest risk classification improvements were for heart failure and all-cause mortality
Sensitivity Analyses
When stratified according to the age cutoff of 65 years, older individuals (n = 34 143; aged ≥65 years) more often had diabetes and hypertension, but smoked less often than younger individuals (n = 129 456; aged <65 years) (eTable 14 in Supplement 1). Concentrations of high-sensitivity cardiac troponin I, high-sensitivity cardiac troponin T, NT-proBNP, BNP, and high-sensitivity CRP were higher, on average, in older individuals.
The association of biomarkers with atherosclerotic cardiovascular disease events remained significant in both individuals younger than 65 years of age and in those aged 65 years or older (eTable 15 in Supplement 1). The subdistribution HRs for high-sensitivity cardiac troponin I, high-sensitivity cardiac troponin T, and NT-proBNP were higher in older individuals. The subdistribution HR was lower for high-sensitivity CRP in older people. In older people, the C statistic of the base model was substantially lower compared with younger people and the addition of the biomarkers provided higher absolute increases of the C statistic in older people (eTable 16 in Supplement 1). For example, the combination of high-sensitivity cardiac troponin I, NT-proBNP, and high-sensitivity CRP increased the C statistic for 10-year atherosclerotic cardiovascular disease events from 0.812 (95% CI, 0.8021-0.8208) to 0.8194 (95% CI, 0.8089-0.8277) in younger people and from 0.6323 (95% CI, 0.5945-0.6570) to 0.6602 (95% CI, 0.6224-0.6834) in older people. Furthermore, the overall NRI with the combined biomarkers for atherosclerotic cardiovascular disease was higher in older people (NRI, 0.062; 95% CI, 0.013-0.120) compared with younger people (NRI, 0.028; 95% CI, 0.010-0.070) (eTable 17 in Supplement 1).
The HRs for all-cause mortality were also higher for NT-proBNP and BNP in older people, but slightly lower for high-sensitivity CRP and high-sensitivity cardiac troponin T (eTable 15 in Supplement 1). A similar pattern was also observed for the other outcomes of heart failure, ischemic stroke, and myocardial infarction. The overall NRI using the combination of high-sensitivity cardiac troponin I, NT-proBNP, and high-sensitivity CRP was higher in older people for the secondary outcomes of heart failure, ischemic stroke, and myocardial infarction, but lower for all-cause mortality (eTable 17 in Supplement 1).
In sensitivity analyses, the regression analyses were repeated for 10-year atherosclerotic cardiovascular disease events while removing 1 risk factor from the base model (eTable 18 in Supplement 1). The addition of high-sensitivity cardiac troponin I, NT-proBNP, and high-sensitivity CRP resulted in the highest increase in the C statistic when age was removed from the base model (C statistic difference of 0.0312; 95% CI, 0.0254-0.0373). The magnitude of change in the C statistic after adding the biomarkers to the full base model was comparable with the effect of systolic blood pressure and self-reported use of antihypertensive medication.
Sensitivity analyses also were performed that included information on cholesterol-lowering medication, which was available in 71.2% of the study population (eTable 19 in Supplement 1). These findings were consistent with the primary study findings.

Discussion

In this individual-level analysis, the value of the most commonly used biomarkers for cardiovascular risk prediction in the general population was investigated using harmonized, multinational population data (from 28 cohorts in 12 countries and 4 continents).
There were 4 salient findings. First, all investigated biomarkers were predictors not only of incident atherosclerotic cardiovascular disease events, but also of all-cause mortality, heart failure, myocardial infarction, and ischemic stroke. Even though prior studies from the general population focused on the association of biomarkers with fatal or nonfatal atherosclerotic cardiovascular disease in general, most did not consider other important outcomes.7,8,22,23 Previous analyses examining incident heart failure or ischemic stroke were limited by (1) small numbers of events, (2) the availability of aggregate data only, or (3) shorter duration of follow-up.7 Interestingly, there was a stronger association of all the investigated biomarkers with all-cause mortality, and particularly with heart failure, compared with fatal and nonfatal atherosclerotic cardiovascular disease events.
In the current dataset, all-cause mortality was the most frequently reported outcome (25 346 events), highlighting the potential competing risk of death for any regression analyses. Thus, the consideration of competing risk by using Fine and Gray regression analyses may be a possible explanation for opposite results compared with prior studies10,22,24 and strengthens the findings from the current analysis. The strong predictive value for heart failure outcomes is particularly notable given the increase in options available for preventing incident heart failure (such as intensive blood pressure control and treatment with sodium-glucose cotransporter 2 inhibitors).25 Importantly, the magnitude of change in the C statistic for the outcomes of all-cause mortality and heart failure reported in the current study is similar to other studies that investigated the addition of coronary calcium scoring to classic cardiovascular risk factors to predict atherosclerotic cardiovascular disease.26
The second significant finding was that the combination of the biomarkers high-sensitivity cardiac troponin I, NT-proBNP, and high-sensitivity CRP into 1 model provided the largest incremental predictive value and that all 3 biomarkers were independent predictors. These 3 biomarkers represent 3 different pathophysiological pathways, had the highest availability in the cohorts examined, are routinely available, and were also identified as the strongest predictors in earlier analyses of multiple biomarkers.27 Most prior studies focused on 1 specific cardiovascular biomarker and did not attempt to combine several markers into 1 model.28 The combination of high-sensitivity cardiac troponin I, NT-proBNP, and high-sensitivity CRP in the current study resulted in the biggest improvements in the C statistic for most outcomes investigated.
Interestingly, the multivariable model showed the highest HRs for NT-proBNP for all outcomes except for incident myocardial infarction for which high-sensitivity CRP showed the strongest association. This ranking of biomarkers is comparable with prior analyses from the FINRISK and Belfast PRIME cohorts,27 for which the highest HRs for atherosclerotic cardiovascular disease events were observed for NT-proBNP; however, no high-sensitivity troponin assay was available for comparison at that time.
The third important finding is a sustained association of improved risk prediction when biomarkers were added to conventional risk factors over a time horizon of more than 10 years. The long follow-up (median duration of nearly 12 years) enabled the assessment of the C statistic over a long time frame, and the incremental value of the biomarkers was apparent even beyond 10 years. This observation highlights the potential value of biomarkers for incorporation in primary prevention strategies, which ideally should address long-term effects. Prior post hoc analyses from the JUPITER, WOSCOPS, and SPRINT large clinical trials29-31 investigated the role of biomarkers, especially cardiac troponin and NT-proBNP for decision-making in preventive care. Data from the JUPITER trial29 showed that those individuals with higher concentrations of high-sensitivity cardiac troponin I or BNP were at higher risk for atherosclerotic cardiovascular disease events and may have a higher absolute risk reduction with statin treatment.
In the WOSCOPS trial,30 longitudinal measurements were available for high-sensitivity cardiac troponin I and showed an association with atherosclerotic cardiovascular disease events and also their decrease 1 year after statin treatment. Recently, post hoc analyses from the SPRINT trial revealed that individuals with elevated concentrations of high-sensitivity cardiac troponin I and NT-proBNP had a substantially increased risk of all-cause mortality and heart failure, but also had the highest absolute risk reduction with treatment compared with individuals with normal concentrations of the biomarkers.31
The fourth novel finding is the greater incremental value of biomarkers in older individuals (aged ≥65 years) compared with younger individuals (aged <65 years). Prior studies showed that with increasing age, the effect of conventional risk factors is attenuated.3,32,33 In the current study, the conventional risk factor model had a C statistic of 0.632 in older individuals vs 0.812 in younger individuals. This resulted in the development of risk prediction models specifically for older people.3,34
However, there remains substantial residual risk when predicting incident atherosclerotic cardiovascular disease events, highlighting the need for other clinically relevant risk markers. In this context, the current study findings support the relevance of cardiovascular biomarkers, especially in older individuals. This observation was primarily driven by the increasing predictive value of NT-proBNP in older people, whereas the predictive value of high-sensitivity CRP decreased. Importantly, these findings were not limited to atherosclerotic cardiovascular disease events, but were also observed for all of the secondary outcomes, especially all-cause mortality and heart failure.

Limitations

This analysis has limitations. First, this study used 5 established biomarkers that are widely available in routine clinical practice; however, the absolute measurements for high-sensitivity cardiac troponin T and BNP were limited. Second, most individuals were recruited from high-income cohorts in Europe and North America, which limits the worldwide generalizability of the findings.
Third, there were a limited number of Black participants. Non-Black participants systematically have higher concentrations of high-sensitivity CRP and higher absolute risks. Fourth, important questions remain before cardiovascular biomarkers may be considered for implementation into clinical practice. These questions include the need for cost-effectiveness analyses and the identification of a target population.

Conclusions

Cardiovascular biomarkers were strongly associated with fatal and nonfatal cardiovascular events and mortality. The addition of biomarkers to established risk factors led to only a small improvement in risk prediction metrics for atherosclerotic cardiovascular disease, but was more favorable for heart failure and mortality.

References 

1.
Visseren  FLJ, Mach  F, Smulders  YM,  et al; ESC National Cardiac Societies; ESC Scientific Document Group.  2021 ESC guidelines on cardiovascular disease prevention in clinical practice.   Eur Heart J. 2021;42(34):3227-3337. doi:10.1093/eurheartj/ehab484PubMedGoogle ScholarCrossref
2.
Arnett  DK, Blumenthal  RS, Albert  MA,  et al.  2019 ACC/AHA guideline on the primary prevention of cardiovascular disease: a report of the American College of Cardiology/American Heart Association task force on clinical practice guidelines.   Circulation. 2019;140(11):e596-e646. doi:10.1161/CIR.0000000000000678PubMedGoogle ScholarCrossref
3.
SCORE2 Working Group and ESC Cardiovascular Risk Collaboration.  SCORE2 risk prediction algorithms: new models to estimate 10-year risk of cardiovascular disease in Europe.   Eur Heart J. 2021;42(25):2439-2454. doi:10.1093/eurheartj/ehab309PubMedGoogle ScholarCrossref
4.
Khan  SS, Matsushita  K, Sang  Y,  et al; Chronic Kidney Disease Prognosis Consortium and the American Heart Association Cardiovascular-Kidney-Metabolic Science Advisory Group.  Development and validation of the American Heart Association’s PREVENT equations.   Circulation. 2024;149(6):430-449. doi:10.1161/CIRCULATIONAHA.123.067626PubMedGoogle ScholarCrossref
5.
Westermann  D, Neumann  JT, Sörensen  NA, Blankenberg  S.  High-sensitivity assays for troponin in patients with cardiac disease.   Nat Rev Cardiol. 2017;14(8):472-483. doi:10.1038/nrcardio.2017.48PubMedGoogle ScholarCrossref
6.
Neumann  JT, Twerenbold  R, Ojeda  F,  et al; COMPASS-MI Study Group.  Application of high-sensitivity troponin in suspected myocardial infarction.   N Engl J Med. 2019;380(26):2529-2540. doi:10.1056/NEJMoa1803377PubMedGoogle ScholarCrossref
7.
Willeit  P, Welsh  P, Evans  JDW,  et al.  High-sensitivity cardiac troponin concentration and risk of first-ever cardiovascular outcomes in 154,052 participants.   J Am Coll Cardiol. 2017;70(5):558-568. doi:10.1016/j.jacc.2017.05.062PubMedGoogle ScholarCrossref
8.
Zeller  T, Tunstall-Pedoe  H, Saarela  O,  et al; MORGAM Investigators.  High population prevalence of cardiac troponin I measured by a high-sensitivity assay and cardiovascular risk estimation: the MORGAM Biomarker Project Scottish Cohort.   Eur Heart J. 2014;35(5):271-281. doi:10.1093/eurheartj/eht406PubMedGoogle ScholarCrossref
9.
Everett  BM, Zeller  T, Glynn  RJ, Ridker  PM, Blankenberg  S.  High-sensitivity cardiac troponin I and B-type natriuretic peptide as predictors of vascular events in primary prevention: impact of statin therapy.   Circulation. 2015;131(21):1851-1860. doi:10.1161/CIRCULATIONAHA.114.014522PubMedGoogle ScholarCrossref
10.
Blankenberg  S, Salomaa  V, Makarova  N,  et al; BiomarCaRE Investigators.  Troponin I and cardiovascular risk prediction in the general population: the BiomarCaRE consortium.   Eur Heart J. 2016;37(30):2428-2437. doi:10.1093/eurheartj/ehw172PubMedGoogle ScholarCrossref
11.
Ridker  PM, Buring  JE, Rifai  N, Cook  NR.  Development and validation of improved algorithms for the assessment of global cardiovascular risk in women: the Reynolds Risk Score.   JAMA. 2007;297(6):611-619. doi:10.1001/jama.297.6.611
ArticlePubMed