Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

PCSK9 loss of function is protective against extra-coronary atherosclerotic cardiovascular disease in a large multi-ethnic cohort

  • Aeron M. Small ,

    Contributed equally to this work with: Aeron M. Small, Jennifer E. Huffman

    Roles Conceptualization, Data curation, Formal analysis, Methodology, Writing – original draft, Writing – review & editing

    Affiliations Corporal Michael J. Crescenz VA Medical Center, Philadelphia, PA, United States of America, Department of Medicine, Yale University School of Medicine, New Haven, CT, United States of America

  • Jennifer E. Huffman ,

    Contributed equally to this work with: Aeron M. Small, Jennifer E. Huffman

    Roles Conceptualization, Formal analysis, Methodology, Writing – review & editing

    Affiliation Massachusetts Veterans Epidemiology Research and Information Center (MAVERIC), VA Boston Healthcare System, Boston, MA, United States of America

  • Derek Klarin,

    Roles Conceptualization, Formal analysis

    Affiliations Boston VA Healthcare System, Boston, MA, United States of America, Center for Genomic Medicine, Massachusetts General Hospital, Harvard Medical School, Boston, MA, United States of America, Program in Medical and Population Genetics, Broad Institute of MIT and Harvard, Cambridge, MA, United States of America, Department of Surgery, Massachusetts General Hospital, Boston, MA, United States of America

  • Julie A. Lynch,

    Roles Writing – review & editing

    Affiliations Department of Veterans Affairs, Salt Lake City Health Care System, Salt Lake City, UT, United States of America, University of Massachusetts College of Nursing & Health Sciences, Boston, MA, United States of America, Center for Healthcare Organization and Implementation Research, Edith Nourse Rogers Memorial VA Hospital, Bedford, MA, United States of America

  • Themistocles Assimes,

    Roles Conceptualization, Writing – review & editing

    Affiliations Department of Medicine, VA Palo Alto Health Care System, Palo Alto, CA, United States of America, Stanford University School of Medicine, Stanford, CA, United States of America

  • Scott DuVall,

    Roles Conceptualization, Writing – review & editing

    Affiliations Department of Veterans Affairs, Salt Lake City Health Care System, Salt Lake City, UT, United States of America, Veterans Affairs Informatics and Computing Infrastructure

  • Yan V. Sun,

    Roles Conceptualization, Writing – review & editing

    Affiliations Atlanta VA Health Care System, Decatur, GA, United States of America, Department of Epidemiology, Emory University Rollins School of Public Health, Atlanta, GA, United States of America, Department of Biomedical Informatics, Emory University School of Medicine, Atlanta, GA, Unites States of America

  • Labiba Shere,

    Roles Writing – review & editing

    Affiliation Palo Alto Epidemiology Research and Information Center for Genomics, VA Palo Alto Health Care System

  • Pradeep Natarajan,

    Roles Conceptualization, Writing – review & editing

    Affiliations Boston VA Healthcare System, Boston, MA, United States of America, Center for Genomic Medicine, Massachusetts General Hospital, Harvard Medical School, Boston, MA, United States of America, Cardiovascular Research Center, Massachusetts General Hospital, Boston, MA, United States of America

  • Michael Gaziano,

    Roles Data curation, Writing – review & editing

    Affiliations Boston VA Healthcare System, Boston, MA, United States of America, Department of Medicine, Brigham and Women’s Hospital, Harvard Medical School, Boston, MA, United States of America

  • Daniel J. Rader,

    Roles Writing – review & editing

    Affiliation Department of Medicine, Perlman School of Medicine, University of Pennsylvania, Philadelphia, PA, United States of America

  • Peter W. F. Wilson,

    Roles Conceptualization, Writing – review & editing

    Affiliations Atlanta VA Health Care System, Decatur, GA, United States of America, Emory Clinical Cardiovascular Research Institute, Atlanta, GA, United States of America

  • Philip S. Tsao,

    Roles Conceptualization, Writing – review & editing

    Affiliations Department of Medicine, VA Palo Alto Health Care System, Palo Alto, CA, United States of America, Stanford University School of Medicine, Stanford, CA, United States of America

  • Kyong-Mi Chang,

    Roles Conceptualization, Writing – review & editing

    Affiliations Corporal Michael J. Crescenz VA Medical Center, Philadelphia, PA, United States of America, Department of Medicine, Perlman School of Medicine, University of Pennsylvania, Philadelphia, PA, United States of America

  • Kelly Cho,

    Roles Conceptualization, Writing – review & editing

    Affiliation Boston VA Healthcare System, Boston, MA, United States of America

  • Christopher J. O’Donnell,

    Roles Conceptualization, Writing – review & editing

    Affiliations Boston VA Healthcare System, Boston, MA, United States of America, Cardiovascular Medicine Division, Department of Medicine, Brigham and Women’s Hospital, Harvard Medical School, Boston, MA, United States of America

  • Juan P. Casas,

    Roles Conceptualization, Writing – review & editing

    Affiliation Boston VA Healthcare System, Boston, MA, United States of America

  • Scott M. Damrauer ,

    Roles Conceptualization, Formal analysis, Investigation, Methodology, Writing – review & editing

    scott.damrauer@va.gov

    Affiliations Corporal Michael J. Crescenz VA Medical Center, Philadelphia, PA, United States of America, Department of Surgery, Perelman School of Medicine, University of Pennsylvania, Philadelphia, PA, United States of America

  •  [ ... ],
  • on behalf of the VA Million Veteran Program

    Membership of the VA Million Veteran Program is listed in the S1 Data.

  • [ view all ]
  • [ view less ]

Abstract

Background

Therapeutic inhibition of PCSK9 protects against coronary artery disease (CAD) and ischemic stroke (IS). The impact on other diseases remains less well characterized.

Methods

We created a genetic risk score (GRS) for PCSK9 using four single nucleotide polymorphisms (SNPs) at or near the PCSK9 locus known to impact lower LDL-Cholesterol (LDL-C): rs11583680, rs11591147, rs2479409, and rs11206510. We then used our GRS to calculate weighted odds ratios reflecting the impact of a genetically determined 10 mg/dL decrease in LDL-C on several pre-specified phenotypes including CAD, IS, peripheral artery disease (PAD), abdominal aortic aneurysm (AAA), type 2 diabetes, dementia, chronic obstructive pulmonary disease, and cancer. Finally, we used our weighted GRS to perform a phenome-wide association study.

Results

Genetic and electronic health record data that passed quality control was available in 312,097 individuals, (227,490 White participants, 58,907 Black participants, and 25,700 Hispanic participants). PCSK9 mediated reduction in LDL-C was associated with a reduced risk of CAD and AAA in trans-ethnic meta-analysis (CAD OR 0.83 [95% CI 0.80–0.87], p = 6.0 x 10−21; AAA OR 0.76 [95% CI 0.68–0.86], p = 2.9 x 10−06). Significant protective effects were noted for PAD in White individuals (OR 0.83 [95% CI 0.71–0.97], p = 2.3 x 10−04) but not in other genetic ancestries. Genetically reduced PCSK9 function associated with a reduced risk of dementia in trans-ethnic meta-analysis (OR 0.86 [95% CI 0.78–0.93], p = 5.0 x 10−04).

Conclusions

Genetically reduced PCSK9 function results in a reduction in risk of several important extra-coronary atherosclerotic phenotypes in addition to known effects on CAD and IS, including PAD and AAA. We also highlight a novel reduction in risk of dementia, supporting a well-recognized vascular component to cognitive impairment and an opportunity for therapeutic repositioning.

Introduction

The discovery of loss-of-function (LoF) variants in PCSK9 provided strong evidence that therapeutic manipulation of PCSK9 can prevent coronary artery disease (CAD) through reduction of low-density lipoprotein cholesterol (LDL-C) [1]. This functional hypothesis has been confirmed by the reduction in CAD and ischemic stroke (IS) reported with PCSK9 inhibition in randomized trials [2, 3]. Evaluation of the association between genetic variation in PCSK9 and diverse outcomes provides the opportunity to identify unexpected effects of PCSK9 inhibition that are not likely to be evident during the short follow up time of clinical trials. Notably, there are several such reports demonstrating an increased risk of type 2 diabetes (T2D) for individuals with PCSK9 LoF [46].

In this study we leveraged the large size and diversity of the VA Million Veteran Program (MVP) to identify the phenotypic consequences of genetic variation in PCSK9 function. We specifically considered CAD, peripheral artery disease (PAD), ischemic stroke (IS), and abdominal aortic aneurysm (AAA), as well as pre-specified non-atherosclerotic diseases, including T2D, dementia, chronic obstructive pulmonary disease (COPD), and cancer, which are common in the clinical population likely to be treated with PCSK9 inhibitors. We subsequently performed a phenome wide association study (PheWAS) in order to identify additional associated phenotypes.

Materials and methods

This study was approved by the Department of Veteran’s Affairs Central Institutional Review Board.

All phenotypes were determined as of enrollment in the MVP. Atherosclerotic phenotypes were developed from a combination of diagnosis and procedure codes present in the VA electronic health record (EHR) of MVP participants (S1 Table in S1 File). Pre-selected non-atherosclerotic phenotypes were developed from a combination of appropriate diagnostic billing codes (S2 Table in S1 File). Phenotyping for PheWAS was performed using the phecode method, described elsewhere [7]. Lipid phenotypes represent maximum (LDL-C, total cholesterol, triglycerides) or minimum (HDL-C) values present in the EHR at any time prior to enrollment. Statin use was adjudicated for all MVP participants as of enrollment.

We created a genetic risk score (GRS) for PCSK9 function using four single nucleotide polymorphisms (SNPs) at or near the PCSK9 locus that have been previously demonstrated in a genetic risk score to strongly associate with lower LDL-C in individuals of European genetic ancestry [4]: rs11583680, rs11591147, rs2479409, and rs11206510.

Using participant level data, we combined these SNPs in a gene-centric score weighted by each variant’s effect on LDL-C within MVP [8] and calculated the association of a one-standard deviation this score with lipid traits (total cholesterol, LDL-C, HDL-C, and triglycerides), controlling for age, sex, and 5 ancestry-specific principal components using linear regression. We then tested the association of a one standard deviation change in the gene-centric score with pre-specified atherosclerotic and non-atherosclerotic traits. Using these results we then calculated adjusted odds ratios to reflect the odds that a phenotype (outcome) occurs given a genetic risk score burden equivalent to a 10 mg/dL difference in LDL-C from the population mean (exposure). We considered a Bonferroni corrected p-value of 6.3x10-3 (p = 0.05/n = 8) significant. Statin use was added as a covariate for sensitivity analysis for all atherosclerotic phenotypes and for select other phenotypes meeting our significance threshold in the primary analysis.

PheWAS was performed using the PheWAS package in R [9]. For PheWAS, we limited our analysis to traits which had greater than 200 cases in all three genetic ancestries (n = 507) and used a Bonferroni corrected p-value of 9.x10-5 (p = 0.05/n = 507) as a significance threshold. All PheWAS analyses were adjusted for age, sex, and 5 ancestry-specific principal components. Ancestry specific PheWAS were meta-analyzed using inverse variance weighted meta-analysis.

Results

Genetic and EHR data passing quality control were available in 312,097 individuals [8] including 227,490 (72.9%) White participants, 58,907 (18.9%) Black participants, and 25,700 (8.2%) Hispanic participants. Minor allele frequencies by genetic ancestry for each SNP in the GRS are described in S3 Table in S1 File (S3 Table in S1 File). For every standard deviation change in GRS, we observed a reduction in LDL-C of 2.65 mg/dL, 1.45 mg/dL, and 2.34 mg/dL in White, Black, and Hispanic individuals, respectively (Fig 1).

thumbnail
Fig 1. Effect of 1 standard deviation difference in PSCK9 genetic risk score on lipid traits in mg/dL.

Beta and 95% confidence interval is displayed for the effect of our PCSK9 genetic risk score on lipid traits for White individuals, Black individuals, Hispanic individuals, and inverse variance weighted meta-analysis of White, Black, and Hispanic populations (IVWMA). All analyses were performed controlling for age, sex, and 5 ancestry-specific principal components.

https://doi.org/10.1371/journal.pone.0239752.g001

Fig 2 summarizes the effect of the genetic variation in PSCK9 on atherosclerotic diseases and Fig 3 summarizes the same for pre-selected non-atherosclerotic endpoints. PCSK9 mediated reductions in LDL-C were associated with an approximately 20% reduced risk of both CAD and AAA in trans-ethnic meta-analysis. Significant protective effects were noted for PAD in White individuals, but not in other genetic ancestries or in trans-ethnic meta-analysis. A nominal, but not experiment-wide, significant association was seen between PCSK9 mediated LDL-C reduction and IS in White individuals and in meta-analysis. Of our non-atherosclerotic endpoints, genetically determined reduced PCSK9 function was associated with a reduced risk of dementia in trans-ethnic meta-analysis. In contrast to previous reports, there was no evidence for increased risk of T2D associated with genetically diminished PCSK9 function.

thumbnail
Fig 2. Weighted odds ratios reflecting the impact of a genetically determined 10 mg/dL decrease in LDL-C on pre-specified atherosclerotic traits.

Odds Ratio and 95% confidence interval is displayed for the effect of a genetically determined 10 mg/dL decrease in LDL-C on primary atherosclerotic traits for White individuals, Black individuals, Hispanic individuals, and inverse variance weighted meta-analysis of White, Black, and Hispanic populations (IVWMA). All analyses were performed controlling for age, sex, and 5 ancestry-specific principal components.

https://doi.org/10.1371/journal.pone.0239752.g002

thumbnail
Fig 3. Weighted odds ratios reflecting the impact of a genetically determined 10 mg/dL decrease in LDL-C on safety and efficacy endpoints.

Odds ratio and 95% confidence interval is displayed for the effect of a genetically determined 10 mg/dL decrease in LDL-C on safety and efficacy endpoints for White individuals, Black individuals, Hispanic individuals, and inverse variance weighted meta-analysis of White, Black, and Hispanic populations (IVWMA). All analyses were performed controlling for age, sex, and 5 ancestry-specific principal components.

https://doi.org/10.1371/journal.pone.0239752.g003

The addition of statin use as a covariate attenuated the strength of association between our GRS and all selected phenotypes. Of the three phenotypes meeting our pre-specified Bonferroni correction in primary analysis (CAD, AAA, and dementia), only CAD and AAA maintained a significant association with our GRS after adjusting for statin use, though both the effect estimates and significance were diminished (S4 Table in S1 File).

A total of 11 traits reached our significance threshold in trans-ancestry PheWAS including traits related to dyslipidemia (hyperlipidemia, disorders of lipid metabolism, hypercholesterolemia, mixed hyperlipidemia), coronary disease (coronary atherosclerosis, ischemic heart disease, other chronic ischemic heart disease, angina pectoris, and unstable angina), and arterial aneurysms (other aneurysm, aortic aneurysm) (S5 Table in S1 File).

Discussion

We generated a GRS predicting PCSK9 function and calculated its impact on risk for a variety of atherosclerotic and non-atherosclerotic phenotypes within the MVP. Our GRS confirm and extend evidence that reduction in LDL-C reduces risk for CAD and IS, observed in previous genetic studies and in clinical trials of PCSK9 inhibitors. We also demonstrated that genetic reduction in PCSK9 function is associated with a reduced risk of other atherosclerotic diseases, including PAD, and AAA. Lastly, we demonstrate a novel association between genetically determined reduction in PCSK9 function and protection against dementia.

AAA is prevalent in 4–8% of the population overall, and is more common in the elderly and in males [10, 11]. We demonstrated that genetic reduction in PCSK9 function significantly reduces the risk of AAA, a finding that was previously described in PheWAS of the rs11591147 variant and Mendelian randomization experiments in MVP [8]. This result is further supported by a meta-analysis of independent experiments using GRSs of lipid trait-associated SNPs to evaluate AAA risk [12]. Taken together, these data highlight a potential role for LDL-C modulating therapies such as PCSK9 inhibition in preventing AAA and identifies avenues for future clinical trials.

Our data further demonstrated that within White individuals, PCSK9 inhibition reduced the risk of atherosclerotic PAD. This finding is consistent with results from the Further Cardiovascular Outcomes Research With PCSK9 Inhibition in Subjects with Elevated Risk (FOURIER) trial, which demonstrated that lowering of LDL-C with Evolocumab reduced the risk of major adverse limb events [13]. Our result extends the FOURIER finding by further suggesting that LDL-C modulating therapies may decrease lifetime burden of PAD. In fact, there are few studies which evaluate the impact of genetic LDL-C modulation on PAD. A 2009 study of 13,634 EUR individuals demonstrated a significant reduction in risk for prevalent PAD, which did not extend to incident disease [3]. We again confirm this finding in White individuals. Results were not significant for Black or Hispanic individuals, perhaps due to lower power to detect associations in these populations.

We found a significant reduction in the risk of all-cause dementia amongst individuals with genetically reduced PCSK9 expression. Both vascular dementia and non-vascular dementias, including Alzheimer Disease, have a recognized vascular component underlying their pathophysiology [14]. In the recent genetic meta-analysis of Alzheimer’s disease, several genes were implicated in lipid processing, including APOM, APOA5, and ABCA1 [15]. Elevated LDL-C increases the risk of both vascular dementia and Alzheimer disease, although treatment with statins has not been shown to consistently reduce the risk of dementia [16]. Our findings add to the literature by providing evidence of a protective effect of LDL-C modulation on the risk of dementia. This hypothesis could be confirmed in additional observational studies and tested in future clinical trials.

We did not replicate the increased risk of diabetes seen in other studies examining pleiotropy of PCSK9 LoF. One possible explanation is heterogeneity in the definitions used for T2D. Whereas several previous studies used a diabetes definition involving laboratory data, ours utilized only ICD codes, allowing for the possibility of case-control misclassification. Further, our study had less power to detect a difference in T2D risk when compared to other larger published meta-analyses [4]. These considerations aside, our study is consistent with results from FOURIER trial subgroup analysis, which demonstrated no difference in risk for new onset diabetes in patients treated with Evolocumab [17].

The addition of statin use as a covariate diminished the effect estimates and strength of association for all atherosclerotic phenotypes. For dementia, the effect estimate was reversed, though the strength of association no longer met our significance threshold. These results highlight the benefits of pharmacologic control of lipids in individuals with higher baseline LDL-C, with HMG CoA reductase inhibition flattening the risk profile between those with and without genetic PCSK9 LoF. Notably, in spite of controlling for statin use, individuals with genetic PCSK9 LoF had further reduction in risk from both CAD and AAA, arguing for an additional benefit of PCSK9 inhibition in modifying risk of atherosclerotic vascular disease in individuals with suboptimal lipid control who are already taking a statin.

Limitations of our work include that the MVP cohort is predominantly male, with reduced power to examine for differences in females. Further, it is possible that pleiotropy could impact our results due to the high phenotypic correlations among atherosclerotic vascular diseases [12]. Additional observational and intervention studies are warranted to more definitely establish causality. Finally, our GRS was established in White population. SNPs included in the score had varying MAFs by genetic ancestry and while results are scaled by the ancestry specific GRS effect on LDL-C, care should be taken in interpreting results in Black and Hispanic populations.

In summary, we generated a GRS predicting PCSK9 function and demonstrated a reduction in risk of several important extra-coronary atherosclerotic phenotypes in addition to known effects on CAD, including PAD, AAA, and ischemic stroke. We also highlight a novel reduction in risk of dementia, supporting a well-recognized vascular component to cognitive impairment. Our work emphasizes the benefits of biobank-driven genetics research and highlights the power of the MVP, a unique and diverse biobank of US veterans.

References

  1. 1. JC C., Sequence Variations in PCSK9, Low LDL, and Protection against Coronary Heart Disease. N Engl J Med, 2006. 23.
  2. 2. Hopewell J.C., et al., Differential effects of PCSK9 variants on risk of coronary disease and ischaemic stroke. Eur Heart J, 2018. 39(5): p. 354–359. pmid:29020353
  3. 3. Folsom A.R., et al., Variation in PCSK9, low LDL cholesterol, and risk of peripheral arterial disease. Atherosclerosis, 2009. 202(1): p. 211–5. pmid:18436227
  4. 4. Schmidt A.F., et al., PCSK9 genetic variants and risk of type 2 diabetes: a mendelian randomisation study. The Lancet Diabetes & Endocrinology, 2017. 5(2): p. 97–105. pmid:27908689
  5. 5. Ference B.A., et al., Variation in PCSK9 and HMGCR and Risk of Cardiovascular Disease and Diabetes. N Engl J Med, 2016. 375(22): p. 2144–2153. pmid:27959767
  6. 6. Ference B.A., Mendelian Randomization Study of ACLY and Cardiovascular Disease. N Engl J Med, 2019. 380. pmid:30865797
  7. 7. Denny J.C., et al., PheWAS: demonstrating the feasibility of a phenome-wide scan to discover gene-disease associations. Bioinformatics, 2010. 26(9): p. 1205–10. pmid:20335276
  8. 8. Klarin D., et al., Genetics of blood lipids among ~300,000 multi-ethnic participants of the Million Veteran Program. Nat Genet, 2018. 50(11): p. 1514–1523. pmid:30275531
  9. 9. Carroll R.J., Bastarache L., and Denny J.C., R PheWAS: data analysis and plotting tools for phenome-wide association studies in the R environment. Bioinformatics, 2014. 30(16): p. 2375–6. pmid:24733291
  10. 10. Scott R.A.P., The Multicentre Aneurysm Screening Study (MASS) into the effect of abdominal aortic aneurysm screening on mortality in men: a randomised controlled trial. The Lancet, 2002. 360(9345): p. 1531–1539. pmid:12443589
  11. 11. Sakalihasan N., Limet R., and Defawe O.D., Abdominal aortic aneurysm. The Lancet, 2005. 365(9470): p. 1577–1589.
  12. 12. Harrison S.C., et al., Genetic Association of Lipids and Lipid Drug Targets With Abdominal Aortic Aneurysm: A Meta-analysis. JAMA Cardiol, 2018. 3(1): p. 26–33. pmid:29188294
  13. 13. Bonaca M.P., et al., Low-Density Lipoprotein Cholesterol Lowering With Evolocumab and Outcomes in Patients With Peripheral Artery Disease: Insights From the FOURIER Trial (Further Cardiovascular Outcomes Research With PCSK9 Inhibition in Subjects With Elevated Risk). Circulation, 2018. 137(4): p. 338–350. pmid:29133605
  14. 14. Arvanitakis Z., et al., Relation of cerebral vessel disease to Alzheimer's disease dementia and cognitive function in elderly people: a cross-sectional study. The Lancet Neurology, 2016. 15(9): p. 934–943. pmid:27312738
  15. 15. Kunkle B.W., et al., Genetic meta-analysis of diagnosed Alzheimer's disease identifies new risk loci and implicates Abeta, tau, immunity and lipid processing. Nat Genet, 2019. 51(3): p. 414–430. pmid:30820047
  16. 16. McGuinness B., Statins for Prevention of Dementia. Cochrane Database of Systematic Reviews, 2016(1). pmid:26727124
  17. 17. Sabatine M.S., et al., Cardiovascular safety and efficacy of the PCSK9 inhibitor evolocumab in patients with and without diabetes and the effect of evolocumab on glycaemia and risk of new-onset diabetes: a prespecified analysis of the FOURIER randomised controlled trial. The Lancet Diabetes & Endocrinology, 2017. 5(12): p. 941–950. pmid:28927706