| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Original Articles |
From the Helicos BioSciences (J.F.T.), Cambridge, Mass; Molecular Medicine (J.F.T., L.S.W., S.A.P.) and Statistical Applications (C.L.H.), Pfizer Global Research and Development, Groton, Conn; Perlegen Sciences (D.A.H., D.R.C.), Mountain View, Calif; and Department of Vascular Medicine (G.K.H., J.J.P.K.), Academic Medical Center, Amsterdam, The Netherlands.
Correspondence to John J.P. Kastelein, MD, PhD, Department of Vascular Medicine, Academic Medical Center, Meibergdreef 9, Room F4-159.2, 1105 AZ Amsterdam, The Netherlands. E-mail j.j.kastelein{at}amc.uva.nl or j.s.jansen@amc.uva.nl
Received September 9, 2008; accepted January 26, 2009.
| Abstract |
|---|
|
|
|---|
Methods and Results— A total of 291 988 single-nucleotide polymorphisms (SNPs) from 1984 individuals were analyzed for association with statin response, followed by genotyping top hits in 3761 additional individuals. None was significant at the whole-genome level in either the initial or follow-up test sets for association with low-density lipoprotein cholesterol, high-density lipoprotein cholesterol, or triglyceride response. In addition to the whole-genome platform, 23 candidate genes previously associated with statin response were analyzed in these 5745 individuals. Three SNPs in apoE were most highly associated with low-density lipoprotein cholesterol response, followed by 1 in PCSK9 with a similar effect size. At the candidate gene level, SNPs in HMGCR were also significant though the effect was less than with those in apoE and PCSK9. rs7412/apoE had the most significant association (P=6x10–30), and its high significance in the whole-genome study (P=4x10–9) confirmed the suitability of this population for detecting effects. Age and gender were found to influence low-density lipoprotein cholesterol response to a similar extent as the most pronounced genetic effects.
Conclusions— Among SNPs tested with an allele frequency of at least 5%, only SNPs in apoE are found to influence statin response significantly. Less frequent variants in PCSK9 and smaller effect sizes in SNPs in HMGCR were also revealed.
Key Words: genetics hypercholesterolemia hydroxymethylglutaryl coenzyme A reductase inhibitors myocardial infarction
| Introduction |
|---|
|
|
|---|
Clinical Perspective see p 173
In our previous examination of atorvastatin response in the Atorvastatin Comparative Cholesterol Efficacy and Safety Study (ACCESS),8 SNPs in only 2 genes, apoE and ABCB1, showed a significant association with response. To extend these data and possibly reveal novel loci, we assessed a much larger cohort from the Treating to New Targets (TNT) trial in which DNA samples were available from 5745 individuals of European ancestry.9 After a 10-mg atorvastatin run-in, individuals were randomized to stay at that dose or titrated to 80 mg and response was assessed in both the run-in at 10 mg and during the titration to 80 mg.
To determine whether additional genes might be involved in statin response, we analyzed 1984 individuals who had been genotyped for
300 000 SNPs and looked for genetic associations with lipid changes. SNPs appearing most significant were then replicated in the remainder of the cohort, 3761 individuals. Independently, all individuals were genotyped for 111 SNPs in 23 candidate genes chosen based on earlier literature reports.
| Methods |
|---|
|
|
|---|
95%, and 94.2% had Hardy-Weinberg P>0.001. Principal components analysis was used to confirm that the sample set did not have substantial population structure. Consistent with the reported ancestries of the sample set, the first principal component seemed to identify European substructure and had strongest loadings on SNPs near the lactase locus. We tested the top 20 components for association with case-control status and with LDL response. None of these tests were significant after correction for the number of comparisons. Candidate gene SNPs were genotyped via SNPlex as described by the manufacturer (Applied Biosystems, Foster City, Calif). To find genetic associations, we fit linear regression models to the quantitative lipid phenotypes. Genotypes were coded as allele counts in regression models. This corresponds to fitting an additive model where each allele makes the same incremental contribution to phenotype. Models included covariates for environmental factors. For LDL-C response, covariates included age, gender, and screening LDL-C levels. No other covariates were found to be significant. The same covariates were used for the HDL-C and TG response except that the screening levels for those lipids instead of LDL-C. We selected transformations for quantitative variables to remove skew and provide more normally distributed residuals in regressions. The primary test for association consisted of an analysis of variance to assess the significance of the amount of variance explained by adding the genotype term to the model. More detailed descriptions of statistical methods and quality control criteria are described elsewhere.10
| Results |
|---|
|
|
|---|
|
A summary of demographic and metabolic phenotypes of all patients is shown in Table 1. Cases (defined as having suffered a primary event) and controls (no events) were genotyped in the whole-genome analysis and the replication cohort (comprising a cohort with no events) who were genotyped for candidate genes and for SNPs identified via the whole-genome scan. For patients remaining at 10 mg after randomization, there was no significant additional change in LDL-C. Individuals switching to 80 mg atorvastatin at 8 weeks lowered their LDL-C by a total of 55.7% and 54.8% at 20 and 60 weeks, respectively. Because only half the individuals were titrated to 80 mg, all such studies have less power than with the full cohort who could be studied at 10 mg.
|
|
For LDL-C response SNPs, 1 assay failed SNPlex design, and 2 others were omitted due to near perfect linkage disequilibrium (LD) with another SNP. No SNP identified in the initial scan replicates to the same degree in the second set. Notably, based on the effect sizes and allele frequencies, each of the top 25 SNPs had at least a 99% power to achieve a P<0.05 level of significance in the complementary TNT sample set if the initial observations were true positives with the same effect size. Only 3 SNPs achieved P<0.05, 2 of which have false discovery rates (FDRs) in excess of 25%. These SNPs fare even worse if Bonferroni corrected. The only SNP that seemed to even marginally retain significance was rs6790122, which achieves a marginal FDR of 0.057 among 22 tests. A further replication was attempted in an independently ascertained set of 1248 individuals from the ACCESS trial.8 The selection criteria for this trial were somewhat different than for TNT, but the LDL-C lowering achieved in the 2 populations was similar and the same population was used in another statin response analysis8 and yielded a significant association with LDL-C response. The power to replicate these findings in ACCESS at P<0.05 ranged from 85% to 95%, yet not a single SNP achieved this level of significance (Table 2). The only marginally promising SNP from TNT, rs6790122, had 90% power to replicate in ACCESS, but failed to do so. Its genotypic effect estimate, while in the same direction, was <30% of that observed in TNT, suggesting a real but weak effect.
Previous studies have identified numerous candidate genes associated with statin response1–4 but very few have replicated. However, some SNPs within apoE have replicated consistently and were expected to be detected. In our whole-genome scan, no SNPs in or near the apoE region were significant even at the P<0.05 level. The most strongly associated SNP in our previous study, rs7412, was not on this platform. To determine whether the absence of an association with SNPs in apoE was due to a lack of coverage in the region or to something unique to this population, 29 SNPs within 50 kb of apoE were genotyped in the cohort. In addition, we had previously tested 43 SNPs in 16 genes8 that had been reported in the literature as associated with lipid response so those SNPs with minor allele frequency >5% were retested in TNT. Earlier testing identified SNPs in additional genes (ADAMTS1, FCAR, NPC1L1, PCSK9, PON1, PPARG, and SCAP1)12–18 as potentially associated with statin response and these were also included as candidates.
The candidate SNP analysis was split into 2 stages. In the first stage, only data from the individuals who were genotyped in the whole-genome scan were analyzed so the findings would be directly comparable with the genome-wide analysis. In the second stage, both the individuals examined in the genome-wide analysis as well as all others in the TNT cohort were also analyzed. Because most candidate SNPs were not associated with statin response in the ACCESS study performed earlier, they were not expected to be associated in TNT either. SNPs that achieved an association with LDL-C response with a significance level of P<0.05, uncorrected for multiple testing, in the complete population are included in Table 3. No candidate SNPs were found to be significantly associated with HDL-C or TG response. All 111 candidate SNPs that were genotyped in the whole cohort in Hardy-Weinberg equilibrium are listed in Supplementary Table V. Thirty-eight of the 111 SNPs were also present on the whole-genome scan platform.
|
To determine whether SNPs in the apoE region could be genotyped as surrogates for rs7412, more SNPs were genotyped across the apoE region in addition to those examined in the whole-genome scan (Table 4). rs7412 was found at a frequency of only 5.6% and was not in LD with any of the whole-genome SNPs. Indeed, the only SNP tested with r2>0.2 with rs7412 was SNP17 (no dbSNP identifier, also known as ENSSNP5761165),19 with r2=0.5 (Supplementary Figure I). The entire region is noteworthy for its very low level of LD and short LD blocks. Thus, this particular region requires extensive genotyping to capture all potential functional variation. The Perlegen 322K used in this study, the Affymetrix 500K, and the Illumina 317K/550K platforms all omit rs7412 and would not detect a weak phenotypic effect. More robust phenotypes could potentially be detected on all platforms via other SNPs in weak LD with rs7412.
|
The most highly associated candidate SNPs other than those within apoE and PCSK9 include 3 within the HMGCR gene. The first 2, rs10474433 and rs17671591, are in high LD with each other (r2=0.99) and the third, rs6453131, is also linked (r2=0.68). In addition to the 3 HMGCR SNPs listed in Table 3 and 3 HMGCR SNPs in Supplementary Table V that were tested in the full set of individuals, we genotyped 12 more HMGCR SNPs in the nongenome-wide analysis subset that included black individuals. rs17238540 (also reported as SNP29) has previously been found associated with statin response in mixed cohorts or blacks.19,20 Except for the SNPs listed in Table 3, we found no association with any other HMGCR SNP in either individuals of European or black ancestry though the population size of the latter group (n=154) was not sufficient for good power. rs3846662 has also been reported to be associated with alternative splicing of HMGCR,21 but we found no association with statin response.
| Discussion |
|---|
|
|
|---|
The strongest genetic effects were observed with apoE and PCSK9, genes previously known to affect statin response. Indeed, in previous studies, only apoE has been identified as associated with LDL-C response across multiple studies in a reproducible manner.5,8 Some studies, generally with very small populations, do not show an association with apoE. The more recently identified PCSK9 is also accumulating reproducible support for an impact on statin response.15,23–25 apoE is found on chylomicrons and VLDL lipoprotein particles in the circulation and can affect their uptake in the liver. Thus, variants of apoE can affect the relative proportion of cholesterol taken up via the diet versus synthesized in the liver. Individuals synthesizing a higher proportion of cholesterol are more susceptible to its inhibition. Similarly, PCSK9 has been found to bind LDLR and affect its degradation, thus impacting cholesterol uptake in the liver and hence its synthesis.
SNPs within HMGCR, the target protein of statins, have also been sometimes found associated with LDL-C response, and we found 3 linked SNPs that achieve significance if examined as candidate SNPs. Two tightly linked SNPs, rs10474433 and rs17671591, are located
30 kb upstream of HMGCR and the third, rs6453131, is located in intron 6. None has been previously reported to be associated with statin response. The effect size of these SNPs is much smaller than that found with apoE and PCSK9 SNPs and, if the association is real, the functional SNP is not known. SNP12 (rs17244841) in intron 5 and SNP29 (rs17238540) in intron 18 are 2 highly linked SNPs that have been reported to be associated with LDL-C response in a mixed cohort with pravastatin19 and a black cohort with simvastatin20 but not associated in a European cohort with atorvastatin.8 The lack of association with atorvastatin was replicated in the present cohort though an effect specific to other statins or with specific cohorts cannot be ruled out. The combination of low frequency of SNPs and small effects makes the study of epistatic interactions unworkable. It would not be unreasonable to expect complex interactions among genes but it would likely require more extremes of response.
Although statins have variable effects among individuals (Figure), most respond if compliant with therapy. Of the >15 000 individuals who completed the 8-week run-in period in this TNT trial, a total of <5% were eliminated due to insufficient LDL-C lowering, noncompliance, or myalgia. Similarly, in the ACCESS trial in which exclusion criteria were not based on response to treatment,26 only 65 of 1427 individuals on 10 mg atorvastatin had <20% LDL-C lowering at 6 weeks. With the exception of the <1% of the population homozygous for the rare apoE or PCSK9 SNPs, a combination of age and gender is more effective at predicting statin response than any of the individual genetic factors (Table 5). These effects could be related to an underlying mechanism or simply due to variable compliance among groups. Similar results were found in the ACCESS statin trial, but other biological and genetic factors might still cause some variation.
|
The lack of strong genetic effects across a large population are likely a reflection of the complexity of lipid homeostasis and, except for individuals with rare, high-effect polymorphisms, variability in response is due to a wide range of small effects superimposed on compliance with recommended dosing. It remains to be determined whether other statins will behave similarly because sufficiently large populations have not been reported yet. With respect to predicting efficacy, genetic markers are not presently useful in setting atorvastatin dose and the driving force should continue to be treatment with an appropriate dose titration to attain the guideline-driven LDL-C goal. Moreover, these data emphasize that even comprehensive genetic analyses might miss important associations when minor allele frequencies are low and phenotypic consequences are modest and further illustrate the limits of the predictive power of genetic analyses for similar questions, especially when restricted to common variants.
| Acknowledgments |
|---|
Pfizer provided funding for the trial and subsequent analyses.
Disclosures
Drs Thompson, Hyde, Wood, and Cox and S.A. Paciga are present or former employees of Pfizer. Drs Hinds and Cox are present or former employees of Perlegen. Dr Kastelein has received grant support, lecture fees, and consulting fees from Pfizer. Dr Hovingh reports no disclosures.
| References |
|---|
|
|
|---|
2. Shitara Y, Sugiyama Y. Pharmacokinetic and pharmacodynamic alterations of 3-hydroxy-3-methylglutaryl coenzyme A (HMG-CoA) reductase inhibitors: drug–drug interactions and interindividual differences in transporter and metabolic enzyme functions. Pharmacol Ther. 2006; 112: 71–105.[CrossRef][Medline]
3. Zineh I, Johnson JA. Pharmacogenetics of chronic cardiovascular drugs: applications and implications. Expert Opin Pharmacother. 2006; 7: 1417–1427.[CrossRef][Medline]
4. Schmitz G, Schmitz-Madry A, Ugocsai P. Pharmacogenetics and pharmacogenomics of cholesterol-lowering therapy. Curr Opin Lipidol. 2007; 18: 164–173.[Medline]
5. Ordovas JM, Mooser V. The APOE locus and the pharmacogenetics of lipid response. Curr Opin Lipidol. 2002; 13: 113–117.[CrossRef][Medline]
6. Campbell H, Rudan I. Interpretation of genetic association studies in complex disease. Pharmacogenomics J. 2002; 2: 349–360.[CrossRef][Medline]
7. McGinnis B, Olson KL, Magid D, Bayliss E, Korner EJ, Brand DW, Steiner JF. Factors related to adherence to statin therapy. Ann Pharmacother. 2007; 41: 1805–1811.
8. Thompson JF, Man M, Johnson KJ, Wood LS, Lira ME, Lloyd DB, Banerjee P, Milos PM, Myrand SP, Paulauskis J, Milad MA, Sasiela WJ. An association study of 43 SNPs in 16 candidate genes with atorvastatin response. Pharmacogenomics J. 2005; 5: 352–358.[CrossRef][Medline]
9. LaRosa JC, Grundy SM, Waters DD, Shear C, Barter P, Fruchart JC, Gotto AM, Greten H, Kastelein JJ, Shepherd J, Wenger NK; Treating to New Targets (TNT) Investigators. Intensive lipid lowering with atorvastatin in patients with stable coronary disease. N Engl J Med. 2005; 352: 1425–1435.
10. Kooner JS, Chambers JC, Aguilar-Salinas CA, Hinds DA, Hyde CL, Warnes GR, Gómez Pérez FJ, Frazer KA, Elliott P, Scott J, Milos PM, Cox DR, Thompson JF. Genome-wide scan identifies variation in MLXIPL associated with plasma triglycerides. Nat Genet. 2008; 40: 149–151.[CrossRef][Medline]
11. Diabetes Genetics Initiative of Broad Institute of Harvard and MIT, Lund University, and Novartis Institutes of BioMedical Research, Saxena R, Voight BF, Lyssenko V, Burtt NP, de Bakker PI, Chen H, Roix JJ, Kathiresan S, Hirschhorn JN, Daly MJ, Hughes TE, Groop L, Altshuler D, Almgren P, Florez JC, Meyer J, Ardlie K, Bengtsson Boström K, Isomaa B, Lettre G, Lindblad U, Lyon HN, Melander O, Newton-Cheh C, Nilsson P, Orho-Melander M, Råstam L, Speliotes EK, Taskinen MR, Tuomi T, Guiducci C, Berglund A, Carlson J, Gianniny L, Hackett R, Hall L, Holmkvist J, Laurila E, Sjögren M, Sterner M, Surti A, Svensson M, Svensson M, Tewhey R, Blumenstiel B, Parkin M, Defelice M, Barry R, Brodeur W, Camarata J, Chia N, Fava M, Gibbons J, Handsaker B, Healy C, Nguyen K, Gates C, Sougnez C, Gage D, Nizzari M, Gabriel SB, Chirn GW, Ma Q, Parikh H, Richardson D, Ricke D, Purcell S. Genome-wide association analysis identifies loci for type 2 diabetes and triglyceride levels. Science. 2007; 316: 1331–1336.
12. Sabatine MS, Ploughman L, Simonsen KL, Iakoubova OA, Kirchgessner TG, Ranade K, Tsuchihashi Z, Zerba KE, Long DU, Tong CH, Packard CJ, Pfeffer MA, Devlin JJ, Shepherd J, Campos H, Sacks FM, Braunwald E. Association between ADAMTS1 matrix metalloproteinase gene variation, coronary heart disease, and benefit of statin therapy. Arterioscler Thromb Vasc Biol. 2008; 28: 562–567.
13. Iakoubova OA, Tong CH, Chokkalingam AP, Rowland CM, Kirchgessner TG, Louie JZ, Ploughman LM, Sabatine MS, Campos H, Catanese JJ, Leong DU, Young BA, Lew D, Tsuchihashi Z, Luke MM, Packard CJ, Zerba KE, Shaw PM, Shepherd J, Devlin JJ, Sacks FM. Asp92Asn polymorphism in the myeloid IgA Fc receptor is associated with myocardial infarction in two disparate populations: CARE and WOSCOPS. Arterioscler Thromb Vasc Biol. 2006; 26: 2763–2768.
14. Pisciotta L, Fasano T, Bellocchio A, Bocchi L, Sallo R, Fresa R, Colangeli I, Cantafora A, Calandra S, Bertolini S. Effect of ezetimibe coadministered with statins in genotype-confirmed heterozygous FH patients. Atherosclerosis. 2007; 194: e116–e122.[CrossRef][Medline]
15. Naoumova RP, Tosi I, Patel D, Neuwirth C, Horswell SD, Marais AD, van Heyningen C, Soutar AK. Severe hypercholesterolemia in four British families with the D374Y mutation in the PCSK9 gene: long-term follow-up and treatment Response. Arterioscler Thromb Vasc Biol. 2005; 25: 2654–2660.
16. Malin R, Laaksonen R, Knuuti J, Janatuinen T, Vesalainen R, Nuutila P, Lehtimäki T. Paraoxonase genotype modifies the effect of pravastatin on high-density lipoprotein cholesterol. Pharmacogenetics. 2001; 11: 625–633.[CrossRef][Medline]
17. Chen S, Tsybouleva N, Ballantyne CM, Gotto AM Jr, Marian AJ. Effects of PPARalpha, gamma and delta haplotypes on plasma levels of lipids, severity and progression of coronary atherosclerosis and response to statin therapy in the lipoprotein coronary atherosclerosis study. Pharmacogenetics. 2004; 14: 61–71.[CrossRef][Medline]
18. Fan YM, Laaksonen R, Janatuinen T, Vesalainen R, Nuutila P, Knuuti J, Lehtimäki T. Effects of pravastatin therapy on serum lipids and coronary reactivity are not associated with SREBP cleavage-activating protein polymorphism in healthy young men. Clin Genet. 2001; 60: 319–321.[CrossRef][Medline]
19. Chasman DI, Posada D, Subrahmanyan L, Cook NR, Stanton VP Jr, Ridker PM. Pharmacogenetic study of statin therapy and cholesterol reduction. JAMA. 2004; 291: 2821–2827.
20. Krauss RM, Mangravite LM, Smith JD, Medina MW, Wang D, Guo X, Rieder MJ, Simon JA, Hulley SB, Waters D, Saad M, Williams PT, Taylor KD, Yang H, Nickerson DA, Rotter JI. Variation in the 3-hydroxyl-3-methylglutaryl coenzyme a reductase gene is associated with racial differences in low-density lipoprotein cholesterol response to simvastatin treatment. Circulation. 2008; 117: 1537–1544.
21. Medina MW, Gao F, Ruan W, Rotter JI, Krauss RM. Alternative splicing of 3-hydroxy-3-methylglutaryl coenzyme A reductase is associated with plasma low-density lipoprotein cholesterol response to simvastatin. Circulation. 2008; 118: 355–362.
22. Thompson JF, Wood LS, Pickering EH, Dechairo B, Hyde CL. High-density genotyping and functional SNP localization in the CETP gene. J Lipid Res. 2007; 48: 434–443.
23. Dubuc G, Chamberland A, Wassef H, Davignon J, Seidah NG, Bernier L, Prat A. Statins upregulate PCSK9, the gene encoding the proprotein convertase neural apoptosis-regulated convertase-1 implicated in familial Hypercholesterolemia. Arterioscler Thromb Vasc Biol. 2004; 24: 1454–1459.
24. Rashid S, Curtis DE, Garuti R, Anderson NN, Bashmakov Y, Ho YK, Hammer RE, Moon YA, Horton JD. Decreased plasma cholesterol and hypersensitivity to statins in mice lacking Pcsk9. Proc Natl Acad Sci USA. 2005; 102: 5374–5379.
25. Berge KE, Ose L, Leren TP. Missense mutations in the PCSK9 gene are associated with hypocholesterolemia and possibly increased response to statin therapy. Arteriol Thromb Vasc Biol. 2006; 26: 1094–1100.
26. Ballantyne CM, Andrews TC, Hsia JA, Kramer JH, Shear C. Correlation of non-high-density lipoprotein cholesterol with apolipoprotein B: effect of 5 hydroxymethylglutaryl coenzyme A reductase inhibitors on non-high-density lipoprotein cholesterol levels. Am J Cardiol. 2001; 88: 265–269.[CrossRef][Medline]
| Footnotes |
|---|
Related Article
Circ Cardiovasc Genet 2009 2: 173-181.
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
|
Home | Subscriptions | Archives | Feedback | Authors | Help | Circulation Journals Home | AHA Journals Home | Search Copyright © 2009 American Heart Association, Inc. All rights reserved. Unauthorized use prohibited. |