|
| |||||||||||||||||||||||||||||||||||||||||||||||
Background Common single-nucleotide polymorphisms (SNPs) that are associated with blood low-density lipoprotein (LDL) or high-density lipoprotein (HDL) cholesterol modestly affect lipid levels. We tested the hypothesis that a combination of such SNPs contributes to the risk of cardiovascular disease.
Methods We studied SNPs at nine loci in 5414 subjects from the cardiovascular cohort of the Malmö Diet and Cancer Study. We first validated the association between SNPs and either LDL or HDL cholesterol and subsequently created a genotype score on the basis of the number of unfavorable alleles. We used Cox proportional-hazards models to determine the time to the first cardiovascular event in relation to the genotype score.
Results All nine SNPs showed replication of an association with levels of either LDL or HDL cholesterol. With increasing genotype scores, the level of LDL cholesterol increased from 152 mg to 171 mg per deciliter (3.9 to 4.4 mmol per liter), whereas HDL cholesterol decreased from 60 mg to 51 mg per deciliter (1.6 to 1.3 mmol per liter). During follow-up (median, 10.6 years), 238 subjects had a first cardiovascular event. The genotype score was associated with incident cardiovascular disease in models adjusted for covariates including baseline lipid levels (P<0.001). The use of the genotype score did not improve the clinical risk prediction, as assessed by the C statistic. However, there was a significant improvement in risk classification with the use of models that included the genotype score, as compared with those that did not include the genotype score.
Conclusions A genotype score of nine validated SNPs that are associated with modulation in levels of LDL or HDL cholesterol was an independent risk factor for incident cardiovascular disease. The score did not improve risk discrimination but did modestly improve clinical risk reclassification for individual subjects beyond standard clinical factors.
These observations suggest two hypotheses. First, a DNA sequence variant that is related to blood lipoprotein levels may influence the risk of cardiovascular disease. Second, since DNA sequence variants represent an index of lifelong exposure to altered lipoprotein levels,12 whereas plasma measurements vary, DNA sequence variants may add predictive information beyond a single measurement of blood lipids.
The potential of such genetic testing is limited by the fact that each sequence variant that has been discovered to date explains only a modest fraction of the variance (2% or less) in lipid levels.13 However, the combination of several lipid-related polymorphisms could improve risk prediction. Accordingly, we sought first to validate the association between single-nucleotide polymorphisms (SNPs) at nine loci with the level of LDL or HDL cholesterol and then to evaluate the ability of a panel of validated lipid-modulating SNPs to predict a first cardiovascular event in a large, community-based cohort.
Methods
Study Sample and Risk Factors
The Malmö Diet and Cancer Study (MDCS) is a community-based prospective epidemiologic cohort of 28,449 subjects who were recruited for baseline examination between 1991 and 1996.14 From this cohort, 6103 subjects were randomly selected to participate in a cardiovascular cohort (MDCS-CC), which seeks to investigate risk factors for cardiovascular disease. The MDCS-CC study protocols were approved by the ethics committee at Lund University. All subjects provided written informed consent.
All subjects underwent a review of their medical history, a physical examination, and a laboratory assessment of cardiovascular risk factors. The status of cigarette smoking was elicited by a self-administered questionnaire and was coded as never, former, or current. Current cigarette smoking was defined as any smoking within the past year. A parental or sibling history of myocardial infarction (regardless of the age at onset) was also elicited by questionnaire. Blood pressure was measured once after 10 minutes of rest in the supine position. Diabetes mellitus was defined by self-report on the basis of a physician's diagnosis, a fasting blood glucose level of more than 6.0 mmol per liter (108 mg per deciliter), or the use of antidiabetes medication. C-reactive protein was measured by a high-sensitivity assay (Tina-quant CRP, Roche Diagnostics). In fasting venous blood samples, we measured total cholesterol, HDL cholesterol, and triglycerides according to standard procedures at the Department of Clinical Chemistry, University Hospital Malmö. Levels of LDL cholesterol were calculated according to Friedewald's formula, with the assignment of a missing value to subjects with a triglyceride level of more than 400 mg per deciliter (4.5 mmol per liter). The use of lipid-lowering and antihypertensive therapy was recorded.
SNP Selection and Genotyping
We selected for study 11 SNPs in 9 genes on the basis of literature reports of individual associations and from a genomewide association study. We selected six SNPs from six loci on the basis of their association with levels of LDL or HDL cholesterol in at least one previous study. These six SNPs were, for association with LDL cholesterol, APOB (apolipoprotein B, rs693),4 PCSK9 (proprotein convertase subtilisin/kexin type 9, rs11591147),6 and LDLR (low-density lipoprotein receptor, rs688)5; and for association with HDL cholesterol, CETP (cholesteryl ester transfer protein, rs1800775),8 LIPC (hepatic lipase, rs1800588),9 and LPL (lipoprotein lipase, rs328).15 We selected five SNPs from a genomewide association study of levels of blood LDL and HDL cholesterol in 2931 subjects from Sweden and Finland (a sample distinct from the subjects in the MDCS-CC).13 In the genomewide association study, we found SNPs in APOB (rs7575840), APOE cluster (apolipoprotein E, rs4420638), HMGCR (3-hydroxy-3-methylglutaryl-coenzyme A reductase, rs12654264), and LDLR (rs1529729) to be nominally associated with LDL cholesterol (P=1x10–4, P=3.4x10–13, P=4x10–4, and P=0.002, respectively) and ABCA1 (ATP-binding cassette subfamily A member 1, rs3890182) with HDL cholesterol (P=3.3x10–5).16
Genotyping was performed with the use of matrix-assisted laser desorption–ionization time-of-flight mass spectrometry on a MassARRAY platform (Sequenom), as described previously.17 All SNPs were in Hardy–Weinberg equilibrium (P>0.001). The genotyping success rate was 96%. Using 15 samples analyzed in quadruplicate, we found the genotyping error rate to be <0.7%.
Outcomes
Cardiovascular events were ascertained through linkage of the 10-digit personal identification number of each Swedish citizen with three registries: the Swedish Hospital Discharge Register, the Swedish Cause of Death Register, and the Stroke Register of Malmö.18 The prespecified composite end point of cardiovascular events was defined as myocardial infarction, ischemic stroke, and death from coronary heart disease. Myocardial infarction was defined on the basis of codes 410 and I21 in the International Classification of Diseases, 9th Revision and 10th Revision (ICD-9 and ICD-10), respectively. Ischemic stroke was defined on the basis of codes 434 or 436 (ICD-9) and I63 or I64 (ICD-10). Before inclusion, ischemic stroke events in the Swedish Hospital Discharge Register were validated with the use of the Stroke Register of Malmö, in which original medical records (including imaging studies, when available) were examined.18 We considered only ischemic strokes and excluded cases of intracerebral or subarachnoid hemorrhage. Death from coronary heart disease was defined on the basis of codes 412 and 414 (ICD-9) or I22–I23 and I25 (ICD-10) in the Swedish Cause of Death Register. Follow-up extended to December 31, 2003.
Statistical Analysis
For each subject, levels of LDL or HDL cholesterol were adjusted for age, sex, and diabetes status to create a multivariable-adjusted residual level of LDL or HDL cholesterol. The multivariable-adjusted residual lipid level for each subject served as the phenotype in genotype–phenotype association analyses. We conducted multivariable linear regression analyses to test the null hypothesis that LDL or HDL cholesterol residuals did not differ according to an increasing copy number in the minor allele. Study subjects who were receiving lipid-lowering medication at the time of the baseline examination were excluded from these analyses.
On the basis of the results of the genotype–phenotype association analyses, we identified nine SNPs that were independently associated with lipid levels. Two SNPs that we initially tested were correlated with other SNPs in the panel and thus excluded.
We assumed that an individual SNP with a modest effect on lipid levels would have limited power to show an effect on incident cardiovascular events. We therefore constructed a genotype score on the basis of the number of unfavorable alleles (those associated with higher LDL cholesterol levels or lower HDL cholesterol levels) that were carried by each subject for each of the nine SNPs. We calculated crude incidence rates of first cardiovascular events according to strata of genotype scores.
We constructed multivariable proportional-hazards models to examine the association between the genotype score and the time to the first cardiovascular event, excluding subjects who had had a previous myocardial infarction or ischemic stroke. We first confirmed that the proportional-hazards assumption was met.19 The hazard ratio for the genotype score as a continuous measure was estimated in a model adjusting for all 14 available baseline covariates. Cumulative incidence curves were constructed according to the genotype score with the use of Cox regression analysis.
To evaluate the ability of the genotype score to classify risk, we plotted receiver-operating-characteristic (ROC) curves for the 14 baseline covariates with or without the genotype score.20 We used the incidence of cardiovascular events at 10 years as the outcome. The C statistic, a measure of the area under the ROC curve, was calculated with and without the genotype score. We also evaluated the ability of the genotype score to reclassify risk across categories outlined by the Adult Treatment Panel III of the National Cholesterol Education Program21 with the use of two separate approaches suggested by Cook22 and Pencina et al.23
All genetic association analyses were performed with the use of SPSS software, version 14.0. The tests for the proportional-hazards assumption were performed with the use of the survival package for R; the ROC curves and C statistics were generated with the use of the ROCR package for R (www.r-project.org).
Results
Baseline Characteristics
Of the 6103 MDCS-CC subjects, 563 did not provide a baseline plasma sample, 56 did not have data regarding levels of either LDL or HDL cholesterol, 18 did not have DNA available, and 52 did not have genotypes for at least one SNP.24 This left 5414 subjects who were eligible for the present investigation.
Baseline characteristics of the subjects are shown in Table 1. The subjects were all of self-reported European ancestry. The mean (±SD) age was 58±6 years, and 59% were women. A small percentage of the subjects were receiving lipid-lowering therapy (2.3%) or had had a prevalent myocardial infarction or stroke (2.3%).
|
For validation of the association between SNPs and lipid levels, 5287 subjects who were not receiving lipid-lowering therapy at baseline were available. We validated the association between each of 11 SNPs with levels of either LDL cholesterol or HDL cholesterol (Table 2). The proportion of residual trait variance that was explained by each SNP was small, ranging from 0.1 to 2.5%. Among the SNPs that were associated with LDL cholesterol levels, the difference in LDL cholesterol level between homozygote classes ranged from 3 mg to 72 mg per deciliter (0.1 to 1.9 mmol per liter). Among the SNPs related to HDL cholesterol, the difference ranged from 3 mg to 5 mg per deciliter.
|
We then constructed a genotype score for each subject by counting the number of unfavorable alleles carried by that subject for each SNP. With zero, one, or two unfavorable alleles possible for each of nine SNPs, the potential genotype score for each subject ranged from 0 to 18. Since only a small number of subjects had a score of 6 or less, these groups were a priori collapsed into one group for the purposes of data display and analysis. In a similar way, subjects with a score of 13 or more were collapsed into one group.
SNPs and Cardiovascular Events
Complete data for all covariates in the multivariable Cox regression analyses were available for 4232 subjects who had not had prevalent myocardial infarction or stroke. During a median of 10.6 years of follow-up, 137 men and 101 women in this cohort had an incident cardiovascular event, prespecified as myocardial infarction (131 subjects), ischemic stroke (96), or death from coronary heart disease (11).
Table 3 shows the lipid levels and crude incidence rates for the first cardiovascular event stratified according to genotype score. LDL cholesterol levels increased from 152 mg per deciliter (3.9 mmol per liter) for those with a score of 6 or less to 171 mg per deciliter (4.4 mmol per liter) for those with a score of 13 or more (P for trend, 2x10–18). HDL cholesterol levels decreased from 60 mg per deciliter (1.6 mmol per liter) for those with a score of 6 or less to 51 mg per deciliter (1.3 mmol per liter) for those with a score of 13 or more (P for trend, 3x10–24). There was an increase in the crude incidence rate of first cardiovascular events according to genotype score from 3.1 per 1000 person-years (among those with a score
6) to 11.0 per 1000 person-years (among those with a score
13).
|
After accounting for baseline lipid levels and a full set of available risk factors, we found that the genotype score remained associated with the first cardiovascular event (adjusted hazard ratio, 1.15 per copy of an unfavorable allele; 95% confidence interval [CI], 1.07 to 1.24; P=3x10–4 by multivariable analysis) (Table 4). As compared with subjects with a genotype score of 9 or less (43.2% of the sample), subjects with a score of 11 or more (33.6%) had an increased risk of a cardiovascular event by a factor of 1.63 (95% CI, 1.21 to 2.19; P=0.001). The predicted cumulative freedom from cardiovascular events according to the genotype score, with adjustment for baseline covariates, is shown in Figure 1.
|
|
|
To further characterize the relationship between genotype and cardiovascular events, we performed several additional exploratory analyses. These included analyses comparing the incremental predictive value of the SNPs combined in the genotype score with that of the SNPs modeled as a set, analyses of the incremental predictive value of the genotype score for individual components of the cardiovascular event end point, analyses of whether specific SNPs that comprised the score appeared to drive the association with cardiovascular disease, and analyses of the separate predictive value of the SNPs related to LDL cholesterol levels and those related to HDL cholesterol levels. All the additional analyses are discussed in the Supplementary Appendix.
Discussion
In a sample of approximately 5400 subjects from a community-based cohort, we observed that SNPs from nine loci were associated with the baseline level of LDL or HDL cholesterol; four of these associations represented a first replication. During more than 10 years of follow-up, we found that a panel of nine lipid-associated polymorphisms also predicted risk for a first cardiovascular event and provided independent information beyond that obtained from baseline lipid levels and other established risk factors. Nevertheless, the genotype score did not improve the accuracy of the clinical risk prediction, as assessed by the C statistic. Among subjects in the intermediate Adult Treatment Panel III risk category, 26% were reclassified into a higher or lower risk category when the genotype score was added to the other covariates.
These results highlight several important principles. First, for continuous, normally distributed quantitative traits like LDL and HDL cholesterol levels, multiple common alleles influence trait variation, with each allele conferring a modest effect.25 Second, although each SNP exerts a modest effect, a combination of SNPs, in aggregate, can have a substantial influence on lipid levels. Third, lipid-associated SNPs may provide incremental information about cardiovascular risk beyond lipid levels.
Why might a panel of SNPs associated with levels of LDL or HDL cholesterol help predict cardiovascular outcomes even after accounting for baseline lipid levels? It has been proposed that a single measurement of blood lipids may inadequately reflect the lifetime exposure to altered lipid levels.12,26 Lipid measurements may be affected by dietary, diurnal, and random factors.27 By contrast, genotypes are invariant over time, can be precisely measured on a single assay, and can exert an effect on lipid levels over a lifetime. An alternative possibility is that SNPs may influence atherogenesis through mechanisms unrelated to either LDL or HDL cholesterol.
When added to conventional risk factors, the current SNP panel did not improve the ability to determine which patients will have cardiovascular disease and which will remain disease-free, although it did modestly improve risk prediction, as measured by clinical risk reclassification. It remains possible that with the addition of more SNPs related to LDL or HDL cholesterol or SNPs directly related to the risk of myocardial infarction (e.g., the chromosome 9p21 SNPs recently associated with coronary artery disease and myocardial infarction28,29,30), a panel of SNPs could prove useful in risk discrimination or targeting therapies.
Several limitations of our study deserve mention. Since this analysis is among the first in the genetic-association literature to report on a panel of lipid-modulating SNPs,31 we may have overestimated the effect size per unfavorable allele. Additional studies should clarify the true effect size per unfavorable allele.
The specific SNP set in this panel may need refinement. We selected a set of loci and SNPs on the basis of previous studies and a strong biologic rationale. However, we did not test other loci meeting these criteria. At each of the selected loci, the particular SNPs that we tested may not actually be causal but rather correlated to the causal SNP.
In addition, due to the rapid progress of sequencing and genomewide association approaches, a flood of new association discoveries is expected in the near future.32,33,34 It will be important to systematically test the contribution of such newly identified loci to future risk-prediction models.
Our study sample was entirely of European ancestry, and some of the SNPs we studied are not polymorphic in other ethnic groups (e.g., rs11591147 at PCSK9 is monomorphic in persons of African ancestry).12 This factor may limit the generalizability of our findings. However, several of the SNPs (including the LPL and LIPC variants) have been shown to relate to levels of LDL or HDL cholesterol in a similar fashion in both whites and nonwhites.35
Moreover, in the interest of simplicity, we assigned an equal weight in our genotype score to each unfavorable allele, regardless of the strength of the effect on cholesterol. However, the effect on risk may be more accurately estimated by weighting the contribution of an allele to the score on the basis of its contribution to levels of LDL or HDL cholesterol and on the relationship of such levels to cardiovascular events. In addition, the effect of each SNP on the risk of cardiovascular disease may not relate to its effect on cholesterol in a linear fashion.
Finally, reclassification analysis may be affected in important ways by the choice of risk categories. For example, end points predicted in Adult Treatment Panel III differed from those in our study, so alternative risk-categorization schemes may be appropriate. To overcome this limitation, we also assessed the integrated discrimination index, a metric that is not dependent on the choice of categories.
In summary, in our analysis, nine SNPs from nine genes were found to be associated with levels of either LDL or HDL cholesterol. This panel of polymorphisms was independently associated with a risk of first myocardial infarction, ischemic stroke, or death from coronary heart disease. The SNP panel did not improve risk discrimination but did modestly improve clinical risk reclassification for individual subjects beyond standard clinical factors.
Supported by a Doris Duke Charitable Foundation Clinical Scientist Development Award, the Fannie E. Rippel Foundation, the Donovan Family Foundation, and a grant (K23-HL083102) from the National Institutes of Health (NIH), all to Dr. Kathiresan; by the Swedish Medical Society, the Swedish Heart and Lung Foundation, the Ernhold Lundströms Research Foundation, the Mossfelt Foundation, and the Region Skane, all to Dr. Melander; by the Swedish Medical Research Council, the Medical Faculty of Lund University, Malmö University Hospital, the Albert Påhlsson Research Foundation, and the Crafoord Foundation, all to Drs. Melander and Orho-Melander; by the Sigrid Juselius Foundation, the Finnish Diabetes Research Foundation, the Folkhalsan Research Foundation, the Clinical Research Institute of Helsinki University Central Hospital, and a Linné grant from the Swedish Research Council, all to Dr. Groop; by a Burroughs Wellcome Fund Clinical Scholar in Translational Research award and a Distinguished Clinical Scholar award from the Doris Duke Charitable Foundation, both to Dr. Altshuler; and by a Doris Duke Charitable Foundation Clinical Scientist Development Award, a Burroughs Wellcome Fund Career Award for Medical Scientists, and a grant (K23-HL080025) from the NIH, all to Dr. Newton-Cheh.
Dr. Kathiresan reports receiving consulting fees from Novartis and Pfizer; Dr. Hirschhorn, having an equity interest in and receiving consulting fees from Correlagen Diagnostics; Dr. Groop, receiving consulting fees from Tethys, lecture fees from Lilly and GlaxoSmithKline, and grant support from Novartis; and Dr. Altshuler, receiving consulting fees from Merck and MPM Capital and grant support from Novartis and having an equity interest in MPM Capital. No other potential conflict of interest relevant to this article was reported.
Source Information
From the Cardiology Division (S.K., C.N.-C.) and the Center for Human Genetic Research and Department of Molecular Biology (D.M.A.), Massachusetts General Hospital; the Departments of Medicine (S.K., D.M.A., C.N.-C.), Pediatrics (J.N.H.), and Genetics (J.N.H., D.M.A.), Harvard Medical School; and the Division of Endocrinology and Genetics, Children's Hospital (J.N.H.) — all in Boston; the Broad Institute of the Massachusetts Institute of Technology and Harvard University, Cambridge, MA (S.K., C.G., N.P.B., J.N.H., D.M.A., C.N.-C.); the Departments of Clinical Sciences, Hypertension, and Cardiovascular Diseases (O.M.), Diabetes and Endocrinology (D.A., C.R., L.G., M.O.-M.), Internal Medicine (G.B.), and Epidemiological Research (B.H.), University Hospital Malmö, Lund University, Malmö, and Chalmers University of Technology, Göteborg (D.A.) — both in Sweden; and Helsinki University Hospital, Helsinki (L.G.).
Drs. Kathiresan and Orho-Melander contributed equally to this article.
Address reprint requests to Dr. Kathiresan at the Cardiovascular Disease Prevention Center, Cardiology Division, Massachusetts General Hospital, 25 New Chardon St., Suite 301, Boston, MA 02114, or at skathiresan{at}partners.org.
References
| |||||||||||||||||||||||||||||||||||||||||||||||
Related Letters:
Cholesterol Gene Polymorphisms and Cardiovascular Events
Vergeer M., El-Harchaoui K., Stroes E. S.G., Tiret L., Tregouet D.-A., Cambien F., Kathiresan S., Newton-Cheh C., Orho-Melander M.
Extract |
Full Text |
PDF
N Engl J Med 2008;
359:92-93, Jul 3, 2008.
Correspondence
This article has been cited by other articles:
HOME | SUBSCRIBE | SEARCH | CURRENT ISSUE | PAST ISSUES | COLLECTIONS | PRIVACY | TERMS OF USE | HELP | beta.nejm.org Comments and questions? Please contact us. The New England Journal of Medicine is owned, published, and copyrighted © 2009 Massachusetts Medical Society. All rights reserved. |