May 2022
Volume 63, Issue 5
Open Access
Clinical and Epidemiologic Research  |   May 2022
Statin Use in Relation to Intraocular Pressure, Glaucoma, and Ocular Coherence Tomography Parameters in the UK Biobank
Author Affiliations & Notes
  • Jihye Kim
    Department of Epidemiology, Harvard T. H. Chan School of Public Health, Boston, Massachusetts, United States
  • Marianne T. Kennedy Neary
    University of Manchester & Manchester Royal Eye Hospital, Manchester, United Kingdom
  • Hugues Aschard
    Department of Epidemiology, Harvard T. H. Chan School of Public Health, Boston, Massachusetts, United States
    Institut Pasteur, Université de Paris, Department of Computational Biology, Paris, France
  • Mathew M. Palakkamanil
    Department of Ophthalmology and Visual Science, Dalhousie University, Halifax, Nova Scotia, Canada
  • Ron Do
    Charles Bronfman Institute for Personalized Medicine, Department of Genetics and Genomics, Icahn School of Medicine at Mount Sinai, New York, New York, United States
  • Janey L. Wiggs
    Department of Ophthalmology, Harvard Medical School, Massachusetts Eye and Ear Infirmary, Boston, Massachusetts, United States
  • Anthony P. Khawaja
    NIHR Biomedical Research Centre at Moorfields Eye Hospital & UCL Institute of Ophthalmology, London, United Kingdom
  • Louis R. Pasquale
    Department of Ophthalmology, Icahn School of Medicine, Mount Sinai, New York, New York, United States
  • Jae H. Kang
    Channing Division of Network Medicine, Brigham and Women's Hospital, Harvard Medical School, Boston, Massachusetts, United States
  • Correspondence: Jihye Kim, Harvard T. H. Chan School of Public Health, 655 Huntington Avenue, Boston, MA 02115, USA; jihyekim@hsph.harvard.edu 
  • Footnotes
    *  JK and MTKN contributed equally to this work as first authors.
  • Footnotes
     APK, LRP and JHK contributed equally to this work as last authors.
Investigative Ophthalmology & Visual Science May 2022, Vol.63, 31. doi:https://doi.org/10.1167/iovs.63.5.31
  • Views
  • PDF
  • Share
  • Tools
    • Alerts
      ×
      This feature is available to authenticated users only.
      Sign In or Create an Account ×
    • Get Citation

      Jihye Kim, Marianne T. Kennedy Neary, Hugues Aschard, Mathew M. Palakkamanil, Ron Do, Janey L. Wiggs, Anthony P. Khawaja, Louis R. Pasquale, Jae H. Kang, for the International Glaucoma Genetics Consortium (IGGC), UK Biobank Eye and Vision Consortium, and Modifiable Risk Factors for Glaucoma Collaboration; Statin Use in Relation to Intraocular Pressure, Glaucoma, and Ocular Coherence Tomography Parameters in the UK Biobank. Invest. Ophthalmol. Vis. Sci. 2022;63(5):31. https://doi.org/10.1167/iovs.63.5.31.

      Download citation file:


      © ARVO (1962-2015); The Authors (2016-present)

      ×
  • Supplements
Abstract

Purpose: The purpose of this study was to evaluate the relationship between statin use and glaucoma-related traits.

Methods: In a cross-sectional study, we included 118,153 UK Biobank participants with data on statin use and corneal-compensated IOP. In addition, we included 192,283 participants (8982 cases) with data on glaucoma status. After excluding participants with neurodegenerative diseases, 41,638 participants with macular retinal nerve fiber layer thickness (mRNFL) and 41,547 participants with macular ganglion cell inner plexiform layer thickness (mGCIPL) were available for analysis. We examined associations of statin use with IOP, mRNFL, mGCIPL, and glaucoma status utilizing multivariable-adjusted regression models. We assessed whether a glaucoma polygenic risk score (PRS) modified associations. We performed Mendelian randomization (MR) experiments to investigate associations with various glaucoma-related outcomes.

Results: Statin users had higher unadjusted mean IOP ± SD than nonusers, but in a multivariable-adjusted model, IOP did not differ by statin use (difference = 0.05 mm Hg, 95% confidence interval [CI] = −0.02 to 0.13, P = 0.17). Similarly, statin use was not associated with prevalent glaucoma (odds ratio [OR] = 1.05, 95% CI = 0.98 to 1.13). Statin use was weakly associated with thinner mRNFL (difference = −0.15 microns, 95% CI = −0.28 to −0.01, P = 0.03) but not with mGCIPL thickness (difference = −0.12 microns, 95% CI = −0.29 to 0.05, P = 0.17). No association was modified by the glaucoma PRS (Pinteraction 0.16). MR experiments showed no evidence for a causal association between the cholesterol-altering effect of statins and several glaucoma traits (inverse weighted variance P ≥ 0.14).

Conclusions: We found no evidence of a protective association between statin use and glaucoma or related traits after adjusting for key confounders.

Statins are a class of lipid-lowering drugs that lower 3-hydroxy-3-methylglutaryl-coenzyme A (HMG-CoA) reductase activity to decrease low density lipoprotein (LDL) cholesterol levels, and they are indicated for the primary prevention of cardiovascular disease.1,2 In the United States, approximately 39.2 million adults aged 40+ years were statin users in 2012 to 2013,1 whereas in the United Kingdom, the prescription prevalence for statins was 128 per 1000 person-years in 2013.3 Biologically, statins have pleiotropic effects,4,5 including effects on nitric oxide synthesis and thrombosis formation, and antioxidant and immunomodulatory effects. Given the widespread use of statins, there is intense interest in their impact on non-cardiovascular outcomes, including eye diseases.6,7 
Studies of statin use and glaucoma have yielded mixed results. A meta-analysis of observational studies published in 2016 suggested an inverse relation between statin use and glaucoma with short-term use.8 However, a subsequent large observational study of health professionals observed that after careful control for chronological age, the relation between statin use and incident primary open-angle glaucoma (POAG) was null,9 and a recent database study in Australia reported an adverse relation between long-term statin use and glaucoma.10 
Given the inconsistency in prior studies, we used the UK Biobank with questionnaire data, intraocular pressure (IOP) measures, ocular coherence tomography (OCT) parameters, and genetic biomarkers to evaluate the relation between statin use and various glaucoma-related traits, to test the hypothesis that statins may protect against glaucoma. We also created a genetic instrument variable that served as a proxy for HMG-CoA activity for a Mendelian Randomization (MR) study, which reflected long term propensity for having lower LDL levels and has been strongly associated with lower risk of myocardial infarction and death from cardiovascular disease.11,12 
METHODS
The UK Biobank
The UK Biobank is a population-based study that has collected health-related data on >500,000 participants aged 40 to 69 years at baseline (2006–2010). We used baseline questionnaire data (http://www.ukbiobank.ac.uk), high throughput genotyping data, ophthalmic data, and OCT imaging data. The UK Biobank was approved by the National Information Governance Board for Health and Social Care and the National Health Service North West Multicenter Research Ethics Committee (reference number 06/MRE08/65). This research was conducted using the UK Biobank Resource under application number 36741. This research study adhered to the tenets of the Declaration of Helsinki. 
Assessment of Statin and Non-Statin Hypolipidemic Drug Use
Trained nurses collected data on prescription drugs or supplements regularly used via in-person interviews at baseline. If a participant indicated on a touchscreen display that they were taking cholesterol lowering drugs, the interviewer recorded the type of medicine used. If the participant indicated they were not using any medications, the interviewer confirmed that the statement was correct. Interviewers did not collect information about dosing or duration of medication use. The statin types used by participants included simvastatin, atorvastatin, rosuvastatin, pravastatin, and fluvastatin. We also analyzed data on use of non-statin hypolipidemic drugs, which included fibric acid agents, bile acid sequestrants, niacin-derivatives, cholesterol absorption inhibitors, omega-3 fats, and other miscellaneous antihyperlipidemic agents. 
IOP and Glaucoma Status Ascertainment
In 2009 to 2013, IOP measurements in both eyes were taken using the Ocular Response Analyzer noncontact tonometer (Reichert Corp., Philadelphia, PA, USA) at 6 locations throughout the United Kingdom.13 Subjects with eye surgery or an eye infection in the prior month were excluded. As a primary outcome, we analyzed corneal-compensated IOP (IOPcc), which is least influenced by measurement artifact due to corneal biomechanics,14 and excluded extreme outliers.15 Participants with a history of glaucoma laser or surgical treatment were excluded as their untreated IOP level would not be captured. For those participants on ocular hypotensive therapy, we adopted the convention of adjusting the measured IOP upward by 30%, as we have done previously.16,17 IOPcc measurements were available for right and left eyes on 127,798 and 127,428 participants, respectively. We used the mean of all available right and left eye values. Ultimately, 118,153 eligible participants had complete data for the assessment of the relation between statin use and IOPcc. 
At baseline (2006–2010), 9198 out of 215,562 participants reported that they had glaucoma. Participants completed a touch screen questionnaire and were considered to have glaucoma if in response to the question, “Has a doctor told you that you have any of the following problems with your eyes?,” they chose glaucoma from the menu. Participants were also considered to have glaucoma if they reported a history of glaucoma surgery or laser treatment on the questionnaire or if they carried an International Classification of Diseases, Ninth Revision, (ICD-9) or ICD, 10th Revision (ICD-10) code for glaucoma (ICD-9 = 365.* or ICD-10 = H40.** [excluding H40.0* and H42.*]). Ultimately, 192,283 participants (including 8982 glaucoma cases) had sufficient data for the assessment of the relationship between statin use and glaucoma. 
OCT Data
Under the auspices of the UK Biobank Eye and Vision Consortium, from 2009 to 2010, there were 67,321 individuals who underwent macular spectral domain OCT (SD-OCT) imaging as part of the baseline examination.18 Loss of retinal ganglion cell (RGC) bodies and their axons are hallmarks of glaucomatous degeneration, and because approximately 50% of RGC bodies reside in the macula region, the macular retinal nerve fiber layer (mRNFL) and macular ganglion cell inner plexiform layer (mGCIPL) thicknesses in the macula are useful glaucoma-related biomarkers.19,20 The Topcon 3D OCT1000 Mark II was used to complete SD-OCT imaging in a dark room without pupil dilation. The 3-dimensional 6 × 6 mm2 macular volume scan mode (512 A scans per B scan and 128 horizontal B scans in a raster pattern) was used for imaging. Both eyes were imaged starting with the right eye. The OCT images were stored as downloadable electronic files from a secure portal. The mRNFL and mGCIPL were segmented using version 1.6.1.1 of the Topcon Advanced Boundary Segmentation algorithm.18 All OCT images with a quality score less than 45 were removed from the data set. Furthermore, as outlined by Khawaja et al.18 using various indicator scores that identify blinks, eye motion artifacts, and segmentation failures, we removed participants in the bottom 20th percentile of data quality. The final data set consisted of 41,638 participants for mRNFL and 41,547 participants for mGCIPL after excluding participants with Alzheimer's disease, Parkinson's disease, and multiple sclerosis who could have non-glaucomatous reductions in RNFL and GCIPL thicknesses. For mRNFL and mGCIPL measures intended for analyses, we used the mean of all available right and left eye values. 
Genotyping Data, Glaucoma Multi-Trait Analysis of Genome Wide Association Study Polygenic Risk Score and MR Experiments
Genetic data on 488,000 participants were generated using 2 closely related genotyping platforms: the Affymetrix UK BiLEVE Axiom Array, which called genotypes at 807,411 markers on approximately 50,000 individuals,21 and the Affymetrix UK Biobank Axiom Array, which generated genotypes at 825,925 markers for approximately 450,00 individuals. The smaller number of available data partially reflects post-study opt outs and other technical factors. Quality controls and imputation (the determination of genotypes at loci by inference and not by direct genotyping) were performed jointly for these platforms, as previously described.15 Imputation was based on genetic architecture ascertained in the 1000 Genomes Project, UK 10K, and the Haplotype Reference Consortium reference panels. After quality control, 92,693,895 genetic markers of 487,442 participants were available. 
To evaluate interactions with a polygenic risk score (PRS), we calculated a PRS for each participant using 2673 independent common single nucleotide polymorphisms (SNPs) associated at P ≤ 0.001 for glaucoma from a multi-trait analysis of genomewide association study (GWAS; MTAG) that included the UK Biobank.22 We applied the effect estimates from the MTAG study to generate a glaucoma 2673-gene variant PRS in the UK Biobank that also predicted earlier age at glaucoma diagnosis, glaucoma progression, and need for surgical intervention in an independent Australian data set.22 We chose the glaucoma MTAG PRS because it best captures the genetic variation for POAG, which is a complex disease with optic nerve vulnerability across the spectrum of IOP. The PRS was derived using a standard weighted sum of individual SNPs (i.e. PRS = \(\mathop \sum \nolimits_{i\ = \ 1}^{2,673} {\hat{\beta }_i} \times {\rm{SN}}{{\rm{P}}_i}\) where \({\hat{\beta }_i}\) is the estimated effect size of SNPi on glaucoma extracted from the aforementioned GWAS).22 We normalized the glaucoma PRS with mean of 0 and standard deviation (SD) of 1 for analyses. 
We also performed an MR analysis using 5 genomewide significant gene 3-HMG-CoA variants associated with lower serum LDL levels.11 Analogous to a randomized controlled trial, MR leverages the random allocation of genotypes and allows for the evaluation of associations between genetically predicted levels of a risk factor (e.g. statin use) and genetic predicted propensity for a disease outcome (e.g. glaucoma), which would assist with inferring the causal effect of the risk factor on the outcome.23 Briefly, MR leverages genotypic differences that exist at conception. The instrumental variable comprises multiple variants that capture the life-time exposure in a dose-response manner. We conducted the MR study in participants of European descent from the UK Biobank plus other cohorts to optimize the ability to find associations with glaucoma-related traits. More details regarding the MR experiments can be found in the Supplementary Methods
Statistical Analysis
We compared baseline characteristics of non-statin users and statin users by using mean (SD) difference for continuous variables and differences in frequencies for categorical variables. When variables had missing values for >5% of the study population, a missing indicator was used. If ≤5% were missing, then the median value (for continuous variables) or the largest category (for categorical variables) was used. To examine associations between statin use and IOP or OCT parameters, we used multiple linear regression models adjusted for covariates. We created five nested models to address confounding (especially confounding by indication) and detection bias by including numerous covariates, including demographic factors, conditions that are indicated for statin use (e.g. diabetes and cardiovascular disease), blood pressure, blood cholesterol level, and beta-blocker use. To evaluate the relationship between statin use and glaucoma, we conducted multiple logistic regression analyses adjusting for the same covariates. 
To assess whether the glaucoma MTAG PRS modified the relation between statin use and the various outcomes, we tested the significance of adding a multiplicative (MTAG PRS*statin) interaction term in models with the main effects of statin use and the MTAG PRS; in these models, we also included the following multiplicative interaction terms to minimize confounding: total cholesterol *MTAG PRS, cardiovascular disease *MTAG PRS, triglyceride *MTAG PRS, age *MTAG PRS, age2 *MTAG PRS, and non-statin hypolipidemic drug use *MTAG PRS.24 
We conducted additional sensitivity analyses: (1) analyses excluding those with glaucoma for analyses of IOP and of OCT parameters, and (2) analyses where the reference group were restricted to people who did not take any hypolipidemic medicines (versus those who were non-statin users). 
RESULTS
The mean age (SD) of the UK Biobank population (n = 118,153) at baseline was 56.5 (8.1) years. The population was 82.1% Caucasian and 54.4% were women. The demographics for the specific study population for each analysis described below were similar (Supplementary Table S1). 
Statin Use and IOP in the UK Biobank
For the assessment of statin use and IOP, we included 118,153 participants (Table 1). Among statin users, the majority (92.7%) reported using either simvastatin or atorvastatin. Statin users were older than non-statin users (61.4 [6.1] versus 55.9 [8.1 years), less likely to be women (38.6% vs. 56.7%), and more likely to be using beta-blockers (20.7% vs. 3.0%). Statin users had higher systolic blood pressure, higher body mass index (BMI) and were more likely to report diabetes. The average total cholesterol among statin users was lower (4.7 [0.9] mmol/L) compared to non-statin users (5.9 [1.0] mmol/L; of note, “high cholesterol” is defined as a total level >6.2 mmol/L). The crude average IOP was higher for statin users versus non-statin users (16.3 [3.9] mm Hg vs. 15.9 [3.8] mm Hg). The comparison of statin users, users of non-statin hypolipidemic drugs, and non-users of hypolipidemic medicines is provided in Supplementary Table S2
Table 1.
 
Characteristics of 118,153 UK Biobank Participants With Intraocular Pressure Measurements According to Statin Use at Baseline (2009–2013)
Table 1.
 
Characteristics of 118,153 UK Biobank Participants With Intraocular Pressure Measurements According to Statin Use at Baseline (2009–2013)
In a basic multivariable model adjusted for age, age-squared, sex, ethnicity, deprivation, spherical equivalent, and non-statin hypolipidemic medication use, statin use was associated with lower IOP (model 1: difference in IOP = −0.13 mm Hg, P = 2.9E-05; Table 2). This association was slightly attenuated with additional adjustment for cigarette smoking, alcohol use, physical activity, and caffeinated beverage consumption (model 2: difference = −0.10 mm Hg, P = 0.0027). The association remained significant after adjustment for covariates related to metabolic syndrome; namely, BMI, systolic blood pressure, diabetes, and HgA1c, as well as cardiovascular disease, which is often co-existent with serum lipid disorders (model 3: difference = −0.12 mm Hg, P = 7.4E-04). However, the relationship became nonsignificant and attenuated after additional adjustment for systemic beta-blocker use (model 4: difference = −0.06 mm Hg, P = 0.10) and remained nonsignificant in model 5 that also controlled for serum total cholesterol level and triglyceride levels (model 5: difference = 0.05 mm Hg; 95% confidence interval [CI] = −0.02 to 0.13, P = 0.17). Similarly, in model 5, specific statin types were not significantly associated with IOP (see Table 2). In sensitivity analysis excluding prevalent glaucoma cases, the relation between statin use and IOP remained nonsignificant (model 5: difference = 0.04 mm Hg, 95% CI = −0.04 to 0.11, P = 0.31; Supplementary Table S3). Results also remained unchanged in an alternative model where the reference group was restricted to participants not on any hypolipidemic agents (versus the larger group of non-statin users); furthermore, use of non-statin hypolipidemic drugs showed no association with IOP in a fully adjusted multivariable model (difference = 0.005 mm Hg, 95% CI = −0.14 to 0.15, P = 0.95; Supplementary Table S4). 
Table 2.
 
Association Between Statin Use (2006–2010) and Intraocular Pressure (Measured in 2009–2013) in the UK Biobank: Multivariable-Adjusted Difference (95% CI) in IOP (mm Hg) by Statin Use
Table 2.
 
Association Between Statin Use (2006–2010) and Intraocular Pressure (Measured in 2009–2013) in the UK Biobank: Multivariable-Adjusted Difference (95% CI) in IOP (mm Hg) by Statin Use
Statin Use and Self-Reported Glaucoma in the UK Biobank
Among 192,283 participants including 8982 cases of self-reported prevalent glaucoma, the relation between statin use and prevalent glaucoma was null in most models (model 5: OR = 1.05, 95% CI = 0.98 to 1.13; Table 3). When the reference group was restricted to participants not on any hypolipidemic treatments, the results were largely unchanged (see Supplementary Table S4). Furthermore, non-statin hypolipidemic treatment was also not associated with prevalent glaucoma in a fully adjusted multivariable model (OR = 1.04, 95% CI = 0.92 to 1.18; see Supplementary Table S4). 
Table 3.
 
Association Between Statin Use (2006–2010) and Prevalent Glaucoma in the UK Biobank (2006–2010): Multivariable-Adjusted Odds Ratio (95% CI)
Table 3.
 
Association Between Statin Use (2006–2010) and Prevalent Glaucoma in the UK Biobank (2006–2010): Multivariable-Adjusted Odds Ratio (95% CI)
Statin Use and Glaucoma-Related OCT Parameters in the UK Biobank
The average mRNL and mGCIPL thicknesses were 28.9 (3.8) microns and 75.2 (5.2) microns, respectively (see Supplementary Table S1). In the basic multivariable model adjusted for age, age-squared, sex, ethnicity, deprivation, spherical equivalent, and non-statin hypolipidemic treatments, statin use was associated with thinner mRNFL (model 1: difference = −0.32 microns, P = 8.6E-10) and mGCIPL (difference = −0.36 microns, P = 1.3E-07; Table 4). In model 5, the results for mRNFL were attenuated and only borderline significant (model 5: difference = −0.14 microns, 95% CI = −0.28 to −0.01, P = 0.03). For mGCIPL, in model 5, results were attenuated, and nonsignificant (model 5: difference = −0.12 microns, 95% CI = −0.29 to 0.05, P = 0.17). These results were nearly identical after excluding prevalent glaucoma cases (see Supplementary Table S3). In models where the reference group was restricted to participants who did not take any hypolipidemic treatments, the results were very similar (see Supplementary Table S4). 
Table 4.
 
Association Between Statin Use (2009–2010) and Macular Retinal Nerve Fiber Layer (mRNFL) and Macular Ganglion Cell Plexiform Layer (mGCIPL) Thicknesses in the Macula Region (2009–2010) Among UK Biobank Participants
Table 4.
 
Association Between Statin Use (2009–2010) and Macular Retinal Nerve Fiber Layer (mRNFL) and Macular Ganglion Cell Plexiform Layer (mGCIPL) Thicknesses in the Macula Region (2009–2010) Among UK Biobank Participants
Genetic Modification of Statin Use – Glaucoma-Related Outcomes
To determine whether the relation between statin use and glaucoma-related outcomes may differ by genetic propensity for glaucoma, we evaluated whether the MTAG PRS * statin interaction term was significant when added to model 5. The glaucoma MTAG PRS did not significantly modify the relation among statin use and IOP, prevalent glaucoma, mRNFL thickness, or mGCIPL thickness (Pinteraction ≥ 0.16; Table 5). Furthermore, the glaucoma MTAG PRS did not significantly modify the relation between non-statin hypolipidemic medicine use and the glaucoma-related outcomes (Pinteraction ≥ 0.24; Supplementary Table S5). 
Table 5.
 
Multivariable-Adjusted Interactions Between Statin Use and the Glaucoma MTAG PRS on Four Glaucoma-Related Outcomes in the UK Biobank*
Table 5.
 
Multivariable-Adjusted Interactions Between Statin Use and the Glaucoma MTAG PRS on Four Glaucoma-Related Outcomes in the UK Biobank*
Mendelian Randomization Analyses
MR analyses showed no evidence of a causal link between HMG-CoA reductase activity and any of POAG, VCDR, IOP, mRNFL, or mGCIPL (Supplementary Table S6 and Supplementary Fig. S1). Only the analysis using the POAG data set showed significant unbalanced horizontal pleiotropy as well as global heterogeneity (Cochran's Q P = 0.05, I2 = 57.2). 
Discussion
In this large study in the UK Biobank, based on the main analyses as well as the MR analysis, we did not observe any beneficial associations of statin use in relation to glaucoma prevalence or glaucoma-related traits (IOP, mRNFL, or mGCIPL); furthermore, the null associations were not modified by genetic predisposition to glaucoma. 
For IOP, our nested models revealed that the relation between statin use and IOP was subject to considerable confounding, but ultimately, in a fully adjusted model, the relationship was null (n = 118,153). Similarly, we did not observe that non-statin hypolipidemic treatment was associated with IOP. Whereas one Singaporean study (n = 10,033) observed that statin use was associated with higher IOP,25 our null findings are consistent with other studies that have observed null associations.8,26 For example, a Canadian pharmaceutical database study (n = 8548) showed no difference in adjunct topical IOP-lowering medication use between statin users and nonusers.27 In our study, one of the strongest confounders for this association was systemic beta-blocker use. This was consistent with results from another British study of 7093 participants showing that the relation between statin use and IOP was no longer significant after controlling for systemic beta-blocker use.28 Another confounder was total cholesterol; indeed, a recent meta-analysis suggested that total cholesterol is a risk factor for higher IOP.29 A strength of our study was the adjustment for systemic beta-blocker use, total cholesterol levels, and other conditions that may be indications for statin use. 
The relation between statin use and glaucoma has also been intensely studied and has yielded very mixed results ranging from inverse,8,3033 to null,9,3436 to adverse.10,37 Different ways of ascertaining statin use, dissimilar definitions of glaucoma, varied study methodologies, and uneven follow-up periods may have played a role in the inconsistent results. Furthermore, it is difficult to address residual confounding and to disentangle the effects of statins from the underlying indications for statin use, such as dyslipidemia, cardiovascular disease, and diabetes. No observational study on this subject is ideal, and our analysis lacked information on dosage or duration of statin use and was a cross-sectional study that assessed glaucoma by self-report. However, our cross-sectional study results were consistent with a meta-analysis that observed that use of statins overall was not associated with glaucoma incidence.8 In addition, our study had a large sample size, was able to control for multiple indications for statins including total serum cholesterol, evaluated multiple glaucoma outcomes, and the main analyses were augmented with an MR study.23 
Our study was among the first to assess the relation between statin use and inner retinal structure.38 Because early glaucoma can affect the structures in the macula region,39 OCT scanning of the macula region can be useful in diagnosing glaucoma19,20; indeed, a previous UK Biobank study reported that higher IOP was associated with a thinner mGCIPL (but there was no relation with mRNFL thickness).18 However, we observed no beneficial associations between statin use and macular OCT structural parameters (mRNFL and mGCIPL). An initial significant adverse association between statin use and thinner measures was attenuated after adjustment for BMI, systolic blood pressure, diabetes, and cardiovascular disease (although for mRNFL, the adverse association remained nominally significant). The findings with OCT structural parameters collectively support the notion that statin use is not associated with glaucoma. 
We also assessed whether a glaucoma MTAG PRS might interact with statin use to influence various glaucoma-related outcomes, but we did not observe any significant interactions. It is important to note that this same MTAG PRS was a powerful tool that predicted earlier age at glaucoma diagnosis, glaucoma disease progression, and the need for glaucoma surgery.22 Furthermore, in MR analyses, we used a genetic instrument to approximate HMG-CoA activity, an instrument that predicted cardiovascular disease11,40 and cancer.41 With MR analyses, we were able to preclude strong associations between statin use with various glaucoma outcomes, further supporting the lack of a strong association of statin use on glaucoma. 
Our study had several limitations. Although we had data on the type of statin use as well as data on non-statin hypolipidemic medicines, we did not have data on past statin use, which would have biased results toward the null. We also lacked data on dosage or duration of statin use, although the MR analyses provided some complementary data representing dose-response relations and life-time exposures. This was a cross-sectional study that evaluated IOP, prevalent glaucoma, and OCT measures obtained at a single timepoint; furthermore, we were unable to assess statin use in relation to peripapillary RNFL thickness measures as this data was unavailable in the UK Biobank. The definition of glaucoma was not highly specific and mainly relied on participant reporting, so our results may have been susceptible to various biases related to misclassification of the outcome. Finally, because <8% of our study population was non-White, our findings may not be generalizable to people of color. 
This study had several strengths. This was a large study with multiple glaucoma traits, covariates, and genotype information. Statin users had lower serum total cholesterol than nonusers, providing a measure of construct validity for the main exposure of interest. We also generated a genetic instrument that mimicked variation in HMG-CoA reductase activity and served as a surrogate of life-time statin use. This was also the first study to assess whether a glaucoma MTAG PRS might modify the relationship between statin use and glaucoma outcomes. 
Overall, our study provides support for the possibility that statin use is not favorably associated with glaucoma-related outcomes and that statin use may not be an effective primary glaucoma prevention strategy. 
Acknowledgments
NEI R01 EY015473; NEI R01 EY032559; NEI R01 EY022305 and unrestricted Challenge Grant from Research to Prevent Blindness (NYC) and The Glaucoma Foundation (NYC). 
Anthony P. Khawaja is Supported by a UK Research and Innovation Future Leaders Fellowship and an Alcon Research Institute Young Investigator Award. Ron Do is supported by the National Institute of General Medical Sciences of the NIH (R35 GM124836) and the National Heart, Lung, and Blood Institute of the NIH (R01 HL139865 and R01 HL155915). The sponsor or funding organization had no role in the design or conduct of this research. 
Disclosure: J. Kim, None; M.T. Kennedy Neary, None; H. Aschard, None; M.M. Palakkamanil, None; R. Do, AstraZeneca (F), Goldfinch Bio (F, S), Pensieve Health (O), Variant Bio (R); J.L. Wiggs, Allergan (R), Editas (R), Broadwing Bio (R), Maze (R), Regenxbio (R), Aerpio (F, R); A.P. Khawaja, Abbvie (R), Google Health (R), Reichert (R), Santen (R), Thea (R); L.R. Pasquale, Twenty Twenty (R), Skye Biosciences (R), Eyenovia (R); J.H. Kang, None 
References
Salami JA, Warraich H, Valero-Elizondo J, et al. National Trends in Statin Use and Expenditures in the US Adult Population From 2002 to 2013: Insights From the Medical Expenditure Panel Survey. JAMA Cardiol. 2017; 2: 56–65. [CrossRef] [PubMed]
Preventive Services Task Force US, Bibbins-Domingo K, Grossman DC, et al. Statin Use for the Primary Prevention of Cardiovascular Disease in Adults: US Preventive Services Task Force Recommendation Statement. JAMA. 2016; 316: 1997–2007. [CrossRef] [PubMed]
O'Keeffe AG, Nazareth I, Petersen I. Time trends in the prescription of statins for the primary prevention of cardiovascular disease in the United Kingdom: a cohort study using The Health Improvement Network primary care data. Clin Epidemiol. 2016; 8: 123–132. [PubMed]
Osmak M. Statins and cancer: current and future prospects. Cancer Lett. 2012; 324: 1–12. [CrossRef] [PubMed]
Blum A, Shamburek R. The pleiotropic effects of statins on endothelial function, vascular inflammation, immunomodulation and thrombogenesis. Atherosclerosis. 2009; 203: 325–330. [CrossRef] [PubMed]
Pranata R, Vania R, Victor AA. Statin reduces the incidence of diabetic retinopathy and its need for intervention: A systematic review and meta-analysis. Eur J Ophthalmol. 2021; 31: 1216–1224. [CrossRef] [PubMed]
Ludwig CA, Vail D, Rajeshuni NA, et al. Statins and the progression of age-related macular degeneration in the United States. PLoS One. 2021; 16: e0252878. [CrossRef] [PubMed]
McCann P, Hogg RE, Fallis R, Azuara-Blanco A. The Effect of Statins on Intraocular Pressure and on the Incidence and Progression of Glaucoma: A Systematic Review and Meta-Analysis. Invest Ophthalmol Vis Sci. 2016; 57: 2729–2748. [CrossRef] [PubMed]
Kang JH, Boumenna T, Stein JD, et al. Notice of Retraction and Replacement. Kang et al. Association of statin use and high serum cholesterol levels with risk of primary open-angle glaucoma. JAMA Ophthalmol. 2020; 138: 588–589. [CrossRef] [PubMed]
Yuan Y, Wang W, Shang X, et al. Association between statin use and the risks of glaucoma in Australia: a 10-year cohort study [published online ahead of print August 4, 2021]. Br J Ophthalmol, https://doi.org/10.1136/bjophthalmol-2021-318789.
Ference BA, Robinson JG, Brook RD, et al. Variation in PCSK9 and HMGCR and Risk of Cardiovascular Disease and Diabetes. N Engl J Med. 2016; 375: 2144–2153. [CrossRef] [PubMed]
Willer CJ, Schmidt EM, Sengupta S, et al. Discovery and refinement of loci associated with lipid levels. Nat Genet. 2013; 45: 1274–1283. [PubMed]
Chua SYL, Thomas D, Allen N, et al. Cohort profile: design and methods in the eye and vision consortium of UK Biobank. BMJ Open. 2019; 9: e025077. [CrossRef] [PubMed]
Luce DA. Determining in vivo biomechanical properties of the cornea with an ocular response analyzer. J Cataract Refract Surg. 2005; 31: 156–162. [CrossRef] [PubMed]
Khawaja AP, Cooke Bailey JN, Wareham NJ, et al. Genome-wide analyses identify 68 new loci associated with intraocular pressure and improve risk prediction for primary open-angle glaucoma. Nat Genet. 2018; 50: 778–782. [CrossRef] [PubMed]
MacGregor S, Ong JS, An J, et al. Genome-wide association study of intraocular pressure uncovers new pathways to glaucoma. Nat Genet. 2018; 50: 1067–1071. [CrossRef] [PubMed]
Hysi PG, Cheng CY, Springelkamp H, et al. Genome-wide analysis of multi-ancestry cohorts identifies new loci influencing intraocular pressure and susceptibility to glaucoma. Nat Genet. 2014; 46: 1126–1130. [CrossRef] [PubMed]
Khawaja AP, Chua S, Hysi PG, et al. Comparison of Associations with Different Macular Inner Retinal Thickness Parameters in a Large Cohort: The UK Biobank. Ophthalmology. 2020; 127: 62–71. [CrossRef] [PubMed]
Kim KE, Park KH. Macular imaging by optical coherence tomography in the diagnosis and management of glaucoma. Br J Ophthalmol. 2018; 102: 718–724. [CrossRef] [PubMed]
Oddone F, Lucenteforte E, Michelessi M, et al. Macular versus Retinal Nerve Fiber Layer Parameters for Diagnosing Manifest Glaucoma: A Systematic Review of Diagnostic Accuracy Studies. Ophthalmology. 2016; 123: 939–949. [CrossRef] [PubMed]
Wain LV, Shrine N, Miller S, et al. Novel insights into the genetics of smoking behaviour, lung function, and chronic obstructive pulmonary disease (UK BiLEVE): a genetic association study in UK Biobank. Lancet Respir Med. 2015; 3: 769–781. [CrossRef] [PubMed]
Craig JE, Han X, Qassim A, et al. Multitrait analysis of glaucoma identifies new risk loci and enables polygenic prediction of disease susceptibility and progression. Nat Genet. 2020; 52: 160–166. [CrossRef] [PubMed]
Davey Smith G, Hemani G. Mendelian randomization: genetic anchors for causal inference in epidemiological studies. Hum Mol Genet. 2014; 23: R89–R98. [CrossRef] [PubMed]
Keller MC. Gene x environment interaction studies have not properly controlled for potential confounders: the problem and the (simple) solution. Biol Psychiatry. 2014; 75: 18–24. [CrossRef] [PubMed]
Ho H, Shi Y, Chua J, et al. Association of Systemic Medication Use With Intraocular Pressure in a Multiethnic Asian Population: The Singapore Epidemiology of Eye Diseases Study. JAMA Ophthalmol. 2017; 135: 196–202. [CrossRef] [PubMed]
Hohn R, Mirshahi A, Nickels S, et al. Cardiovascular medication and intraocular pressure: results from the Gutenberg Health Study. Br J Ophthalmol. 2017; 101: 1633–1637. [CrossRef] [PubMed]
Iskedjian M, Walker JH, Desjardins O, et al. Effect of selected antihypertensives, antidiabetics, statins and diuretics on adjunctive medical treatment of glaucoma: a population based study. Curr Med Res Opin. 2009; 25: 1879–1888. [CrossRef] [PubMed]
Khawaja AP, Chan MP, Broadway DC, et al. Systemic medication and intraocular pressure in a British population: the EPIC-Norfolk Eye Study. Ophthalmology. 2014; 121: 1501–1507. [CrossRef] [PubMed]
Wang S, Bao X. Hyperlipidemia, Blood Lipid Level, and the Risk of Glaucoma: A Meta-Analysis. Invest Ophthalmol Vis Sci. 2019; 60: 1028–1043. [CrossRef] [PubMed]
Marcus MW, Muskens RP, Ramdas WD, et al. Cholesterol-lowering drugs and incident open-angle glaucoma: a population-based cohort study. PLoS One. 2012; 7: e29724. [CrossRef] [PubMed]
McGwin G, Jr., McNeal S, Owsley C, Girkin C, Epstein D, Lee PP. Statins and other cholesterol-lowering medications and the presence of glaucoma. Arch Ophthalmol. 2004; 122: 822–826. [CrossRef] [PubMed]
Stein JD, Newman-Casey PA, Talwar N, Nan B, Richards JE, Musch DC. The relationship between statin use and open-angle glaucoma. Ophthalmology. 2012; 119: 2074–2081. [CrossRef] [PubMed]
Talwar N, Musch DC, Stein JD. Association of Daily Dosage and Type of Statin Agent With Risk of Open-Angle Glaucoma. JAMA Ophthalmol. 2017; 135: 263–267. [CrossRef] [PubMed]
Owen CG, Carey IM, Shah S, et al. Hypotensive medication, statins, and the risk of glaucoma. Invest Ophthalmol Vis Sci. 2010; 51: 3524–3530. [CrossRef] [PubMed]
Zheng W, Dryja TP, Wei Z, et al. Systemic Medication Associations with Presumed Advanced or Uncontrolled Primary Open-Angle Glaucoma. Ophthalmology. 2018; 125: 984–993. [CrossRef] [PubMed]
Shon K, Sung KR. Dyslipidemia, Dyslipidemia Treatment, and Open-angle Glaucoma in the Korean National Health and Nutrition Examination Survey. J Glaucoma. 2019; 28: 550–556. [CrossRef] [PubMed]
Chen HY, Hsu SY, Chang YC, et al. Association Between Statin Use and Open-angle Glaucoma in Hyperlipidemia Patients: A Taiwanese Population-based Case-control Study. Medicine (Baltimore). 2015; 94: e2018. [CrossRef] [PubMed]
Marshall H, Mullany S, Qassim A, et al. Cardiovascular Disease Predicts Structural and Functional Progression in Early Glaucoma. Ophthalmology. 2021; 128: 58–69. [CrossRef] [PubMed]
Hood DC, Raza AS, de Moraes CG, et al. Glaucomatous damage of the macula. Prog Retin Eye Res. 2013; 32: 1–21. [CrossRef] [PubMed]
Schooling CM, Zhao JV, Au Yeung SL, Leung GM. Investigating pleiotropic effects of statins on ischemic heart disease in the UK Biobank using Mendelian randomisation. Elife. 2020; 9: e58567. [CrossRef] [PubMed]
Carter P, Vithayathil M, Kar S, et al. Predicting the effect of statins on cancer risk using genetic variants from a Mendelian randomization study in the UK Biobank. Elife. 2020; 9: e57191. [CrossRef] [PubMed]
Table 1.
 
Characteristics of 118,153 UK Biobank Participants With Intraocular Pressure Measurements According to Statin Use at Baseline (2009–2013)
Table 1.
 
Characteristics of 118,153 UK Biobank Participants With Intraocular Pressure Measurements According to Statin Use at Baseline (2009–2013)
Table 2.
 
Association Between Statin Use (2006–2010) and Intraocular Pressure (Measured in 2009–2013) in the UK Biobank: Multivariable-Adjusted Difference (95% CI) in IOP (mm Hg) by Statin Use
Table 2.
 
Association Between Statin Use (2006–2010) and Intraocular Pressure (Measured in 2009–2013) in the UK Biobank: Multivariable-Adjusted Difference (95% CI) in IOP (mm Hg) by Statin Use
Table 3.
 
Association Between Statin Use (2006–2010) and Prevalent Glaucoma in the UK Biobank (2006–2010): Multivariable-Adjusted Odds Ratio (95% CI)
Table 3.
 
Association Between Statin Use (2006–2010) and Prevalent Glaucoma in the UK Biobank (2006–2010): Multivariable-Adjusted Odds Ratio (95% CI)
Table 4.
 
Association Between Statin Use (2009–2010) and Macular Retinal Nerve Fiber Layer (mRNFL) and Macular Ganglion Cell Plexiform Layer (mGCIPL) Thicknesses in the Macula Region (2009–2010) Among UK Biobank Participants
Table 4.
 
Association Between Statin Use (2009–2010) and Macular Retinal Nerve Fiber Layer (mRNFL) and Macular Ganglion Cell Plexiform Layer (mGCIPL) Thicknesses in the Macula Region (2009–2010) Among UK Biobank Participants
Table 5.
 
Multivariable-Adjusted Interactions Between Statin Use and the Glaucoma MTAG PRS on Four Glaucoma-Related Outcomes in the UK Biobank*
Table 5.
 
Multivariable-Adjusted Interactions Between Statin Use and the Glaucoma MTAG PRS on Four Glaucoma-Related Outcomes in the UK Biobank*
×
×

This PDF is available to Subscribers Only

Sign in or purchase a subscription to access this content. ×

You must be signed into an individual account to use this feature.

×