A medication-wide association study to identify medications associated with incident clinically significant diabetic retinopathy

Abstract

Background:

Diabetic retinopathy, a common microvascular complication of diabetes mellitus, is one of the leading causes of vision loss worldwide. Although some oral drugs have been suggested to affect the risk of diabetic retinopathy, systematic evaluation about the associations between medications and diabetic retinopathy is still absent.

Objective:

To comprehensively investigate associations of systemic medications with incident clinically significant diabetic retinopathy (CSDR).

Design:

Population-based cohort study.

Methods:

From 2006 to 2009, more than 26 000 participants residing in New South Wales were enrolled in the 45 and Up study. Diabetic participants with self-reported physician diagnosis or records of anti-diabetic medication prescriptions were finally included in the current analysis. CSDR was defined as diabetic retinopathy cases requiring retinal photocoagulation recorded in the Medicare Benefits Schedule database from 2006 to 2016. Prescriptions of systemic medication from 5 years to 30 days prior to CSDR were retrieved from the Pharmaceutical Benefits Scheme. The study participants were equally split into training and testing datasets. Logistic regression analyses were performed for the association between each of systemic medication and CSDR in the training dataset. After controlling the false discovery rate (FDR), significant associations were further validated in the testing dataset.

Results:

The 10-year incidence of CSDR was 3.9% (n = 404). A total of 26 systemic medications were found to be positively associated with CSDR, among which 15 were validated by the testing dataset. Additional adjustments for pertinent comorbidities suggested that isosorbide mononitrate (ISMN) (OR: 1.87, 95%CI: 1.00–3.48), calcitriol (OR: 4.08, 95% CI: 2.02–8.24), three insulins and analogues (e.g., intermediate-acting human insulin, OR: 4.28, 95% CI: 1.69–10.8), five anti-hypertensive medications (e.g., furosemide, OR: 2.53, 95% CI: 1.77–3.61), fenofibrate (OR: 1.96, 95% CI: 1.36–2.82) and clopidogrel (OR: 1.72, 95% CI: 1.15–2.58) were independently associated with CSDR.

Conclusion:

This study investigated the association of a full spectrum of systemic medications with incident CSDR. ISMN, calcitriol, clopidogrel, a few subtypes of insulin, anti-hypertensive and cholesterol-lowering medications were found to be associated with incident CSDR.

Keywords

45 and Up Study diabetic retinopathy medication-wide association study

Introduction

Diabetic retinopathy (DR) is one of the most common microvascular complications of diabetes mellitus (DM) and the leading cause of acquired vision loss worldwide among the working-age population.¹ With the increasing global burden of DM and population aging, an estimated 245 million people will have DR globally by 2045.² Visual impairment and blindness due to DR have a significant impact on the quality of life and impose a tremendous socioeconomic burden for both individuals and society.³

With the rapid development of the pharmaceutical industry, an increasing number of medications are available for use and polypharmacy has become more common. This leads to growing interest in the associations between systemic medications and DR. Recently, the effects of anti-hypertensive,^4,5 cholesterol control,^6–9 anti-diabetic,¹⁰ and anti-depressant¹¹ drugs on DR have been intensively investigated. However, previous studies mostly focused on only one or a few systemic medications and estimated medication exposures based on questionnaires, which might be subject to selective reporting bias and recall bias.

Genome-wide associated study (GWAS) is an unbiased and statistically powerful method and has been widely used for identifying genotype–phenotype associations through testing millions of genetic variations. We have proposed a medication-wide association study (MWAS) which is analogous to GWAS in basic principles. This approach allows us to comprehensively search for and validate the associations between all available medications and a specific disease, thus avoiding selective reporting bias.¹² A comprehensive understanding of the association between systemic medications and DR may provide insights into the underlying mechanisms of DR and facilitate the development of new treatments.

Therefore, we aimed to comprehensively investigate associations of systemic medications with clinically significant diabetic retinopathy (CSDR) in working-age Australian populations.

Methods

Data sources

This study was conducted based on the data obtained from the Sax Institution’s 45 and Up Study, which is the largest population-based cohort study in the Southern Hemisphere. The methodology of this study has been described in detail elsewhere.¹³ In brief, participants aged 45 years and above residing in New South Wales were randomly sampled from the Services Australia (formerly the Australian Government Department of Human Services) Medicare enrollment database and invited to join the study by mail. Baseline recruitment commenced in February 2006 and was completed by the end of 2009. Participants joined the study by completing questionnaires and written consent forms for linkage of their medical information to health databases and long-term follow-up. A total of 267,153 participants, about 10% of the New South Wales general population aged 45 years and above, were enrolled in the study at the baseline, and the overall response rate was estimated to be 18%.

The demographic and social characteristics and general health-related data were collected from the self-report questionnaire at baseline, which is available at https://www.saxinstitute.org.au/solutions/45-and-up-study/use-the-45-and-up-study/data-and-technical-information/. Medical information was extracted from two Australian healthcare databases, the Medicare Benefits Schedule (MBS) and Pharmaceutical Benefits Scheme (PBS) databases, which were linked to the 45 and Up Study participants by the use of an identification number assigned by the Services Australia. The MBS database records all the information on claims for medical and diagnostic services provided through Medicare, while the PBS database records the claims for prescribed medications. Both the MBS and PBS data were available from January 2004 to December 2016.

Study population

Diabetic participants aged 45–64 years were included in the present study. The methods used to identify diabetes status were consistent with previous studies, which had been proven to be of high sensitivity and specificity.^14,15 Specifically, participants who reported physician-diagnosed diabetes or had records of anti-diabetic medication prescriptions before the baseline recruitment were identified as diabetic participants. Exclusion criteria of this study included (1) gestational diabetes, defined as no anti-diabetic medication records after the last childbirth but report a history of diabetes in the questionnaire;¹⁴ (2) a history of vitrectomy prior to 1 January 2006 (MBS code: 42725); and (3) a history of retinal photocoagulation (RPC) treatment prior to 1 January 2006 (MBS code: 42809). The process of participant selection is summarized in Figure 1.

Figure 1.

Methodology to identify eligible participants and scan for systemic medications associated with clinically significant diabetic retinopathy.

Definition of CSDR

CSDR, representing proliferative diabetic retinopathy (PDR) or diabetic macular edema (DME) requiring laser therapy, was defined as incident RPC treatment from 1 January 2006 to 31 December 2016 (MBS code, 42809). The index date was defined as the date of the first RPC treatment.

Definition of medication exposure

We obtained medication information from the PBS database for each participant during their identification period, that is, from 5 years to 30 days prior to the index date for participants with CSDR, and from study entry to 31 December 2016 for those without CSDR. Medications prescribed for ophthalmic indications or used by less than 1% of the included participants were excluded from the analysis. Systemic medications were analyzed by the fifth level code of the Anatomical Therapeutic Chemical (ATC) classification system. Participants with at least three fillings of a particular systemic medication during the identification period were classified as exposed to that medication; participants without such records were classified as non-exposed. Those with one or two fillings were classified as having uncertain exposure and were excluded from the analysis for that particular medication analysis.

Potential confounders

Potential confounders included sociodemographic characteristics (age, gender, education level, household income per year, health insurance status, and remoteness index) and comorbidities (history of hypertension, dyslipidemia, cardiovascular diseases [CVD], and osteoporosis). Age was classified into four groups: 45–49, 50–54, 55–59, and 60–64 years. Household income per year was categorized into four groups: <AU$20,000, AU$20,000–AU$40,000, AU$40,000–AU$70,000, and >AU$70,000. Education level was divided into three classes: no qualifications, completion of a certificate or diploma or trade but lower than a university degree, and a university degree or above. Health insurance types included private with extras, private with no extras, veteran’s card, health care concession card, and others. The remoteness index was classified into major cities, inner regions, outer regions, and remote and very remote areas. A history of CVD was defined as a self-reported physician diagnosis of any type of heart disease or stroke, or any treatment received in the last month for heart attack, angina, or other heart diseases. Histories of hypertension, dyslipidemia, and osteoporosis were based on information from questionnaires.

Statistical analysis

The software used for statistical analysis in the study included SAS version 9.4 and R version 3.3.0 (R Foundation for Statistical Computing, www.R-project.org). To validate our results, a random 50:50 split was performed to divide participants into the training dataset and the testing dataset (Figure 1). Baseline characteristics of the participants were reported using means and standard deviations (SDs) for continuous variables and numbers and proportions for categorical variables. Student’s t test for continuous data was used for comparison between participants in the training and testing dataset, while Pearson’s chi-square test was used for the comparison of categorical data. A two-sided p value < 0.05 was considered statistically significant.

Logistic regression models (model 1) adjusting for age, gender, education level, household income, health insurance status, and remoteness index were used separately to assess the association of each systemic medication with incident CSDR. To reduce the risk of overadjustment, pertinent comorbidities including the history of hypertension, dyslipidemia, CVD, and osteoporosis were adjusted for in further analyses (model 2) of those medications which were significantly associated with CSDR in the first analyses (model 1). The odds ratios (ORs) and 95% confidence intervals (CIs) were calculated. In the training dataset, the Benjamin–Hochberg procedure was used to control the false discovery rate (FDR) at the level of 5%. If a systemic medication achieved FDR <5% significant in the training dataset and nominal statistical significance in the testing dataset (p < 0.05), it was deemed tentatively validated.

To improve the generalizability of our findings, a sensitivity analysis was performed to investigate the association between tentatively validated systemic medications and incident CSDR among participants stratified by the remoteness index. Participants who lived in major cities were selected as the training dataset, while the remaining participants were classified as the testing dataset. Subgroup analyses were conducted to explore whether pertinent comorbidities could modify the associations between systemic medications and incident CSDR. The interaction between medication and pertinent comorbidity was considered statistically significant when p value for their interaction was less than 0.05 and then a further subgroup analysis would be conducted.

Results

Study population

Among 267,153 participants at baseline, a total of 10,339 participants met the inclusion and exclusion criteria. Figure 1 provides more details on the participant selection process. The mean age of included participants was 57.2 years with a standard deviation of 5.2 years, and 45.2% of participants were female. Comparisons of baseline participants’ characteristics between the training and testing datasets are shown in Table 1. Participants in the training and testing datasets had similar baseline characteristics.

Table 1.

Baseline characteristics of participants in the training and testing datasets.

	Total	Training dataset	Testing dataset	p value
No.	10339	5170	5169	–
Age (years), %
45–49	1206 (11.7)	610 (11.8)	596 (11.5)	0.858
50–54	2138 (20.7)	1058 (20.5)	1080 (20.9)
55–59	3115 (30.1)	1569 (30.3)	1546 (29.9)
60–64	3880 (37.5)	1933 (37.4)	1947 (37.7)
Female, %	4675 (45.2)	2340 (45.3)	2335 (45.2)	0.928
Household income (AU$/year), %
<20,000	2383 (23.0)	1185 (22.9)	1198 (23.2)	0.768
20,000–40,000	1693 (16.4)	832 (16.1)	861 (16.7)
40,000–70,000	1955 (18.9)	994 (19.2)	961 (18.6)
>70,000	2301 (22.3)	1164 (22.5)	1137 (22.0)
Education, %
No qualifications	1581 (15.3)	798 (15.4)	783 (15.1)	0.925
Certificate/diploma/trade	6635 (64.2)	3307 (64.0)	3328 (64.4)
University degree	1988 (19.2)	994 (19.2)	994 (19.2)
Health insurance type, %
Private with extras	4175 (40.4)	2118 (41.0)	2057 (39.8)	0.243
Private no extras	843 (8.2)	408 (7.9)	435 (8.4)
Veterans card	133 (1.3)	70 (1.4)	63 (1.2)
Concession card	3226 (31.2)	1613 (31.2)	1613 (31.2)
None of above	1769 (17.1)	864 (16.7)	905 (17.5)
Remoteness index, %
Major cities	5479 (53.0)	2767 (53.5)	2712 (52.5)	0.235
Inner regional	3436 (33.2)	1711 (33.1)	1725 (33.4)
Outer regional	1109 (10.7)	530 (10.3)	579 (11.2)
Remote or very remote	123 (1.2)	56 (1.1)	67 (1.3)
Incidence of CSDR, n	3.9 (404)	3.8 (197)	4.0 (270)	0.610

CSDR, clinically significant diabetic retinopathy.

Associations between medications and incident CSDR

During the study period, the number of participants who developed CSDR was 404 (3.9%). There was no significant difference in the incidence of CSDR between training and testing datasets (3.8% versus 4.0%, p = 0.610). Participants in the present study were exposed to 209 prescribed systemic medications in total. The three most commonly prescribed systemic medications in the training dataset were metformin, atorvastatin, and gliclazide, which were used by 73%, 49%, and 38% of participants, respectively. All systemic medications used by the participants in the training dataset are presented in Supplementary Table 1.

A total of 26 systemic medications were significantly associated with incident CSDR (FDR <5%) in the training dataset, and 15 of them were tentatively validated (p < 0.05) in the testing dataset. Figure 2 and Supplementary Table 1 show the associations between each of the analyzed systemic medications and incident CSDR in the training dataset. The results for associations between systemic medications and incident CSDR in the testing dataset are presented in Supplementary Table 2.

Figure 2.

Manhattan plot shows results for the association between systemic medications and incident clinically significant diabetic retinopathy in the training dataset.

Table 2 presents the 15 validated systemic medications positively associated with incident CSDR. After further adjustments for demographic characteristics and pertinent comorbidities mentioned above (model 2), insulins and its analogs were associated with the largest risk of developing CSDR. The risk of CSDR in this medication subcategory ranged from the insulin glargine (OR: 2.22, 95% CI: 1.08–4.54) to intermediate-acting human insulin (OR: 4.28, 95% CI: 1.69–10.8). Five anti-hypertensive medications, including two diuretics (e.g., furosemide, OR: 2.53, 95% CI: 1.77–3.61), one calcium channel blocker (CCB) (amlodipine, OR: 1.47, 95% CI: 1.02–2.11), and two α-antiadrenergic agents (e.g. moxonidine, OR: 2.41, 95% CI: 1.51–3.85), were significantly associated with an increased risk of CSDR. One cholesterol-lowering agent (fenofibrate, OR: 1.96, 95% CI: 1.36–2.82), one vasodilator (isosorbide mononitrate [ISMN], OR: 1.87, 95% CI: 1.00–3.48), one vitamin D analog (calcitriol, OR: 4.08, 95% CI: 2.02–8.24), and one platelet aggregation inhibitor (clopidogrel, OR: 1.72, 95% CI: 1.15–2.58) also had significant positive associations with incident CSDR.

Table 2.

Systemic medications with positive associations with clinically significant diabetic retinopathy progression in the training and testing datasets.

Drug name	Numbers of participants^a		Multiple-adjusted model 1^b				Multiple-adjusted model 2^c
	Training dataset	Testing dataset	Training dataset		Testing dataset		Training dataset		Testing dataset
	N = 5170	N = 5169	OR (95%CI)	FDR	OR (95%CI)	p value	OR (95%CI)	p value	OR (95%CI)	p value
Insulin human, intermediate-acting	287	645	4.62 (3.16–6.75)	2.21E−13	3.58 (2.63–4.88)	7.05E−16	11.12 (4.01–30.80)	<1.00E−20	4.28 (1.69–10.83)	2.16E−03
Insulin lispro	139	1171	4.40 (2.60–7.45)	1.38E−06	3.19 (2.40–4.24)	1.88E−15	4.24 (0.86–20.95)	7.67E−02	3.76 (1.14–12.36)	2.92E−02
Insulin aspart	582	306	4.30 (3.13–5.89)	2.99E−17	3.28 (2.20–4.89)	5.20E−09	9.92 (3.87–25.40)	<1.00E−20	3.15 (1.45–6.87)	3.86E−03
Insulin human, fast-acting	93	655	4.12 (2.19–7.74)	2.57E−04	2.48 (1.76–3.50)	2.20E−07	21.10 (4.69–95.00)	7.00E−05	3.23 (0.72–14.54)	1.27E−01
Insulin glargine	1094	130	3.77 (2.82–5.05)	5.51E−17	2.90 (1.59–5.28)	4.96E−04	5.66 (2.38–13.46)	9.00E−05	2.22 (1.08−4.54)	2.91E−02
Insulin human, intermediate- or long-acting	204	126	2.75 (1.66–4.55)	1.04E−03	2.34 (1.23–4.45)	9.53E−03	3.58 (0.70–18.25)	1.25E−01	2.97 (0.82–10.74)	9.62E−02
Fenofibrate	622	194	2.21 (1.56–3.14)	2.45E−04	3.88 (2.46–6.13)	5.88E−09	2.15 (1.50–3.07)	2.60E−05	1.96 (1.36–2.82)	3.22E–04
Furosemide	680	569	3.24 (2.32–4.54)	3.82E–10	1.98 (1.38–2.84)	1.90E–04	3.28 (2.32–4.63)	<1.00E−20	2.53 (1.77–3.61)	<1.00E−20
Spironolactone	324	759	2.50 (1.63–3.84)	4.42E−04	1.51 (1.06–2.15)	2.34E−02	2.51 (1.62–3.90)	4.00E−05	2.18 (1.40–3.40)	5.72E−04
Amlodipine	747	70	2.05 (1.46–2.86)	4.42E−04	4.16 (2.07–8.36)	6.16E−05	1.97 (1.40–2.78)	1.06E−04	1.47 (1.02–2.11)	3.83E−02
Moxonidine	245	321	2.59 (1.62–4.13)	9.93E−04	2.23 (1.44–3.44)	3.03E−04	2.45 (1.52–3.95)	2.27E−04	2.41 (1.51–3.85)	2.45E−04
Prazosin	280	594	2.01 (1.23–3.28)	4.36E−02	1.63 (1.12–2.38)	1.14E−02	1.88 (1.14–3.09)	1.32E−02	2.16 (1.35–3.46)	1.42E−03
ISMN	183	180	3.25 (1.97–5.38)	1.29E−04	1.85 (1.02–3.36)	4.17E−02	3.35 (1.97–5.72)	9.00E−06	1.87 (1.00–3.48)	4.90E−02
Clopidogrel	594	294	1.91 (1.32–2.77)	6.50E−03	2.21 (1.39–3.52)	8.61E−04	1.98 (1.34–2.93)	6.54E−04	1.72 (1.15–2.58)	8.62E−03
Calcitriol	72	554	3.59 (1.74–7.42)	6.50E−03	2.63 (1.88–3.69)	1.81E−08	3.49 (1.68–7.28)	8.41E−04	4.08 (2.02–8.24)	8.60E−05

CI, confidence interval; FDR, false discovery rate; ISMN, isosorbide mononitrate; OR, odds ratio. Boldface indicates statistical significance.

Number of participants: the number of participants exposed to each drug in the dataset.

Model 1: Adjusting for demographic covariates, including age, gender, educational level, household income, health insurance status, and remoteness index.

Model 2: Adjusting for covariates in model 1 plus hypertension, dyslipidemia, cardiovascular disease, and osteoporosis.

Sensitivity and subgroup analysis

In the sensitivity analysis (Supplementary Table 3), the results remained consistent with those of the main analyses apart from fast-acting human insulin (urban dataset: OR: 2.02, 95% CI: 0.96–4.23, p = 0.063) and prazosin (rural dataset: OR: 1.34, 95% CI: 0.77–2.34, p = 0.299). In the subgroup analyses, no interaction effects between medications and their pertinent comorbidities were statistically significant, except for fenofibrate and dyslipidemia (p-value for interaction = 0.007). Fenofibrate was consistently associated with an increased risk of CSDR for participants without dyslipidemia (OR: 2.79, 95% CI: 2.04–3.81), but no significant association was found for those with dyslipidemia (OR: 1.29, 95% CI: 0.84–1.98) (Supplementary Table 4).

Discussion

The present study was the first to comprehensively explore associations of systemic medications with incident CSDR using an approach analogous to GWAS. Exposures to subtypes of anti-diabetic (human insulin, insulin aspart, and insulin glargine), anti-hypertensive (diuretics, CCBs, and α-antiadrenergic agents), and cholesterol-lowering drugs (fenofibrate) were found to be associated with increased risks of developing CSDR. Novel associations of incident CSDR with ISMN and calcitriol were also identified, which had not been reported before.

One of the challenges in interpreting the association between medication use and disease outcome is to differentiate whether this association is caused by the medication or the disease for which the medication was indicated. Our study found that the use of diuretics (furosemide and spironolactone), CCBs (amlodipine), or α-antiadrenergic agents (moxonidine) was associated with an increased risk of developing CSDR. However, no significant association was found between the risk of developing CSDR and the most prescribed anti-hypertensive drug (perindopril: number of participants with perindopril exposure: 2416; training dataset: OR: 1.63, 95% CI: 1.20–2.21, FDR = 0.017; testing dataset: OR: 1.23, 95% CI: 0.89–1.69, p = 0.214). In addition, after adjusting for comorbidities including hypertension, the association between these anti-hypertensive drugs and incident CSDR remained significant. No significant interaction effects between hypertension and anti-hypertensive drugs were found. These findings support the hypothesis that these subtypes of anti-hypertensive drugs, rather than hypertension, were associated with increased risks of CSDR.¹⁶ Similar findings were also observed for cholesterol-lowering drugs. After adjusting for relevant comorbidities such as CVD and hypertension, the associations remained significant.

We found that exposure to ISMN was associated with an increased risk of CSDR, which had not been reported before. Previous studies reported that the long-term use of oral nitrate might result in adverse cardiovascular events,^17,18 possibly through nitrate tolerance, endothelial dysfunction, and oxidative stress.^17,19,20 These pathophysiological processes may partly explain the positive association between ISMN and CSDR, as systemic and retinal vessels have many similarities in structure and physiology,²¹ although further studies are needed to explore the exact mechanisms.

Calcitriol, the active form of vitamin D, was shown to be positively associated with CSDR in our study. Some indirect evidence has suggested that calcitriol may have an inhibitory role in the pathogenesis of DR. For example, a few epidemiological studies have found an increased risk of DR among participants with vitamin D deficiency.^22–24 Vitamin D receptor and the enzymes involved in vitamin D metabolisms, such as 1a-hydroxylase and vitamin D3–inactivating enzyme 24-hydroxylase, were detected in the retina, suggesting it may function as a paracrine/autocrine regulator.^25–27 In addition, its effects in reducing retinal neovascularization^28,29 and inflammation,^30,31 and improving glucose homeostasis and insulin sensitivity³² provide further support for this association. However, the role of calcitriol supplementation in the formation or progression of DR still lacks direct evidence. We might speculate that calcitriol or additional vitamin D may not be beneficial to DR as patients taking calcitriol may have had a deficiency of vitamin D levels. Nevertheless, we could not obtain information on the serum vitamin D level and could not adjust for conditions requiring calcitriol treatment other than osteoporosis (i.e. renal disease).

Consistent with previous studies, insulins and its analogs were significantly associated with CSDR.³³ One possible explanation is the ‘early worsening’ phenomenon, which is the deterioration in DR precipitated by rapid reductions in serum glucose level and upregulation of vascular endothelial growth factor (VEGF) levels.^33,34 Other potential explanations for this association include the insulin-associated vascular leakage mediated by betacellulin and signaling via the epidermal growth factor receptor³⁵ or the fact that insulin users may be more likely to have suboptimal glycemic control than those who were not on insulin. Regarding anti-hypertensive drugs, their associations with CSDR remain inconclusive.^4,36,37 Positive associations of diuretics and CCBs with CSDR observed in our study were partially consistent with findings from Lin and Lai,⁴ who were able to show angiotensin-converting enzyme inhibitors and/or angiotensin receptor blockers, CCBs, and diuretics were associated with higher risks of developing sight-threatening diabetic retinopathy in comparison with β-blockers. The discrepancy in associations between various anti-hypertensive drugs and CSDR may be partly explained by their differential effects on glycemic and metabolic control.³⁷ Spironolactone has been reported to have a glycemic worsening potential of 0.2% on patients’ hemoglobin A1c (HbA1c),³⁸ and this poor glycemic control may have negative effects on CSDR. As for clopidogrel, there is no direct evidence from human research on associations between clopidogrel and CSDR, although aspirin, another kind of antiplatelet agent, had no impact on DR progression or the rates of vitreous hemorrhage.³⁹ Findings on the association between fibrates and DR remain controversial.^6,8,9,40,41 Two large clinical trials showed that fenofibrate slowed the progression of DR and reduced the need for laser treatment in type 2 diabetic patients with pre-existing retinopathy,^9,40 whereas two real-world observational studies did not find such associations.^6,8 The discrepancy between our findings and those of previous studies may be explained by different patient populations, different confounding factors controlled for, and the lack of data on lipid profiles in our study.

To the best of our knowledge, we were the first to investigate the association of a full spectrum of medications with incident CSDR using the MWAS approach. The MWAS approach used in the present study has several advantages. First, compared with previous studies that explored only one or a few medications at a time, MWAS enables investigations of associations between a wide range of medications and the target disease, thus reducing the risk of selective reporting bias. Second, the internal cross-validation by splitting randomly or based on the remoteness index in the included participants and FDR controlling add to the robustness of the results. Third, MWAS is a systematic and thorough method that may be applied to a wide range of other diseases in future studies. Last but not the least, novel associations identified by MWAS can provide insights into the biochemical and pathogenesis pathways behind diseases, as well as help to identify new therapeutic targets and potential side effects of a specific drug. A comprehensive understanding of all available medications’ association with a specific disease can guide clinicians and pharmacists in making management decisions. In addition to the advantages of the MWAS methodology, other strengths of the present analysis include the large sample size, the wide range of covariates adjusted for in the final models, and the use of the national claim data to estimate the exposure to systemic medications.

Several limitations of our study should also be noted. First, due to the data-driven and hypothesis-independent nature of MWAS, it is hard to determine whether these identified associations are casual or due to relevant comorbidities, which is also a limitation of other observational studies. The list of CSDR-associated systemic medications needs to be replicated using an independent sample and validated by further prospective interventional trials. Second, fundus imaging data were not available in this study; hence, the records of RPC were used to define CSDR, leading to potential biases. Some patients with CSDR may be treated with anti-VEGF or vitrectomy. However, according to the PBS schedule and Australian guidelines for the management of diabetic retinopathy in 2008, aflibercept or ranibizumab must be used in combination with RPC when treating DME, and RPC is also recommended before vitrectomy. Thus, almost all the patients receiving intravitreal injections or vitrectomy should be covered within the cases we tracked by RPC, although the actual number might be relatively underestimated. Third is the fact that information on other baseline systematic diseases (i.e. renal disease), types of diabetes, blood tests (i.e. HbA1c), and reasons for the prescription was not available in the 45 and Up study and was hence not taken into account in this analysis. The possibility of residual confounding could not be completely ruled out in this research, despite the fact that as many confounders as possible were adjusted for in this analysis. Further well-designed, hypothesis-dependent studies are required to validate the current findings. Finally, records from the PBS database could only provide information on the prescription of medications rather than the use of medications by participants. In future studies, electronic medication monitoring devices may be able to more accurately estimate the exact medication exposure.

In summary, the present analysis is the first to systematically explore associations between a full spectrum of systemic medications and CSDR to date. The novel positive associations observed between ISMN or calcitriol and incident CSDR may inspire future research into the relevant biochemical and pathogenic pathways underlying DR. The associations of insulin, anti-hypertensive, and cholesterol-lowering medications with incident CSDR may guide future clinical practice. Additional observational and experimental studies are needed to further validate these associations.

Supplemental Material

sj-docx-1-oed-10.1177_25158414221139002 – Supplemental material for A medication-wide association study to identify medications associated with incident clinically significant diabetic retinopathy

Supplemental material, sj-docx-1-oed-10.1177_25158414221139002 for A medication-wide association study to identify medications associated with incident clinically significant diabetic retinopathy by Ruilin Xiong, Wei Wang, Xianwen Shang, Yixiong Yuan, Yifan Chen, Lei Zhang, Katerina V Kiburg, Zhuoting Zhu and Mingguang He in Therapeutic Advances in Ophthalmology

Footnotes

Acknowledgements

This research was completed using data collected through the 45 and Up Study (). The 45 and Up Study is managed by the Sax Institute in collaboration with major partner Cancer Council NSW and the following partners: the Heart Foundation; NSW Ministry of Health; NSW Department of Communities and Justice; and Australian Red Cross Lifeblood. We acknowledge that Services Australia supplied the Medicare Benefits Schedule (MBS) and Pharmaceutical Benefits Scheme (PBS) data to the Sax Institute. We thank the many thousands of people participating in the 45 and Up Study.

Declarations

ORCID iD

Katerina V Kiburg

Supplemental material

Supplemental material for this article is available online.

References

Klein

BE.

Overview of epidemiologic studies of diabetic retinopathy. Ophthalmic Epidemiol 2007; 14: 179–183.

IDF. IDF diabetes atlas 8th edition, 2017, https://diabetesatlas.org/atlas/eighth-edition/

DAEP Ltd. The economic impact of diabetic macular oedema in Australia, 2015, https://www2.deloitte.com/au/en/pages/economics/articles/economic-impact-diabetic-eye-disease.html

Lin

Lai

MS.

Antihypertensive drugs and diabetic retinopathy in patients with type 2 diabetes. Ophthalmologica 2016; 235: 87–96.

Chaturvedi

Porta

Klein

, et al. Effect of candesartan on prevention (DIRECT-Prevent 1) and progression (DIRECT-Protect 1) of retinopathy in type 1 diabetes: randomised, placebo-controlled trials. Lancet 2008; 372: 1394–1402.

Kawasaki

Konta

Nishida

Lipid-lowering medication is associated with decreased risk of diabetic retinopathy and the need for treatment in patients with type 2 diabetes: a real-world observational analysis of a health claims database. Diabetes Obes Metab 2018; 20: 2351–2360.

Kang

Chen

Garg

, et al. Association of statin therapy with prevention of vision-threatening diabetic retinopathy. JAMA Ophthalmol 2019; 137: 363–371.

Vail

Callaway

Ludwig

, et al. Lipid-lowering medications are associated with lower risk of retinopathy and ophthalmic interventions among United States patients with diabetes. Am J Ophthalmol 2019; 207: 378–384.

Keech

Mitchell

Summanen

, et al. Effect of fenofibrate on the need for laser treatment for diabetic retinopathy (FIELD study): a randomised controlled trial. Lancet 2007; 370: 1687–1697.

10.

Tang

Zhao

, et al. Comparisons of diabetic retinopathy events associated with glucose-lowering drugs in patients with type 2 diabetes mellitus: a network meta-analysis. Diabetes Obes Metab 2018; 20: 1262–1279.

11.

Yekta

Xie

Bogner

, et al. The association of antidepressant medications and diabetic retinopathy among people with diabetes. J Diabetes Complications 2015; 29: 1077–1084.

12.

Tzoulaki

Patel

Okamura

, et al. A nutrient-wide association study on blood pressure. Circulation 2012; 126: 2456–2464.

13.

Banks

Redman

Jorm

, et al. Cohort profile: the 45 and Up Study. Int J Epidemiol 2008; 37: 941–947.

14.

Comino

Tran

Haas

, et al. Validating self-report of diabetes use by participants in the 45 and Up Study: a record linkage study. BMC Health Serv Res 2013; 13: 481.

15.

Yan

Han

, et al. Does daily dietary intake affect diabetic retinopathy progression? 10-year results from the 45 and Up Study. Br J Ophthalmol 2019; 104: 1774–1780.

16.

Zheng

Dryja

Wei

, et al. Systemic medication associations with presumed advanced or uncontrolled primary open-angle glaucoma. Ophthalmology 2018; 125: 984–993.

17.

Yiu

Pong

Siu

, et al. Long-term oral nitrate therapy is associated with adverse outcome in diabetic patients following elective percutaneous coronary intervention. Cardiovasc Diabetol 2011; 10: 52.

18.

Nakamura

Moss

Brown

, et al. Long-term nitrate use may be deleterious in ischemic heart disease: a study using the databases from two large-scale postinfarction studies. Am Heart J 1999; 138: 577–585.

19.

Oelze

Knorr

Kroller-Schon

, et al. Chronic therapy with isosorbide-5-mononitrate causes endothelial dysfunction, oxidative stress, and a marked increase in vascular endothelin-1 expression. Eur Heart J 2013; 34: 3206–3216.

20.

Thomas

DiFabio

Gori

, et al. Once daily therapy with isosorbide-5-mononitrate causes endothelial dysfunction in humans: evidence of a free-radical-mediated mechanism. J Am Coll Cardiol 2007; 49: 1289–1295.

21.

Sun

Wang

Mackey

, et al. Retinal vascular caliber: systemic, environmental, and genetic associations. Surv Ophthalmol 2009; 54: 74–95.

22.

Herrmann

Sullivan

Veillard

, et al. Serum 25-hydroxyvitamin D: a predictor of macrovascular and microvascular complications in patients with type 2 diabetes. Diabetes Care 2015; 38: 521–528.

23.

Kaur

Donaghue

Chan

, et al. Vitamin D deficiency is associated with retinopathy in children and adolescents with type 1 diabetes. Diabetes Care 2011; 34: 1400–1402.

24.

Luo

Gao

Qin

LL.

The association between vitamin D deficiency and diabetic retinopathy in type 2 diabetes: a meta-analysis of observational studies. Nutrients 2017; 9: 307.

25.

Jamali

Wang

Darjatmoko

, et al. Vitamin D receptor expression is essential during retinal vascular development and attenuation of neovascularization by 1, 25(OH)2D3. PLoS ONE 2017; 12: e0190131.

26.

Alsalem

Patel

Susarla

, et al. Characterization of vitamin D production by human ocular barrier cells. Invest Ophthalmol Vis Sci 2014; 55: 2140–2147.

27.

Skowron

Pawlicka

Gil

The role of vitamin D in the pathogenesis of ocular diseases. Folia Med Cracov 2018; 58: 103–118.

28.

Ren

Zhao

, et al. The impact of 1,25-dihydroxy vitamin D3 on the expressions of vascular endothelial growth factor and transforming growth factor-β in the retinas of rats with diabetes. Diabetes Res Clin Pract 2012; 98: 474–480.

29.

Merrigan

Kennedy

BN.

Vitamin D receptor agonists regulate ocular developmental angiogenesis and modulate expression of dre-miR-21 and VEGF. Br J Pharmacol 2017; 174: 2636–2651.

30.

Lee

Rekhi

Hoh Kam

, et al. Vitamin D rejuvenates aging eyes by reducing inflammation, clearing amyloid beta and improving visual function. Neurobiol Aging 2012; 33: 2382–2389.

31.

Tang

Zhou

Luger

, et al. Calcitriol suppresses antiretinal autoimmunity through inhibitory effects on the Th17 effector response. J Immunol 2009; 182: 4624–4632.

32.

Yaribeygi

Maleki

Sathyapalan

, et al. The molecular mechanisms by which vitamin D improve glucose homeostasis: a mechanistic review. Life Sci 2020; 244: 117305.

33.

Chantelau

Kimmerle

Meyer-Schwickerath

Insulin, insulin analogues and diabetic retinopathy. Arch Physiol Biochem 2008; 114: 54–62.

34.

Early worsening of diabetic retinopathy in the Diabetes Control and Complications Trial. Arch Ophthalmol 1998; 116: 874–886.

35.

Sugimoto

Cutler

Shen

, et al. Inhibition of EGF signaling protects the diabetic retina from insulin-induced vascular leakage. Am J Pathol 2013; 183: 987–995.

36.

Wang

Zhang

, et al. Effects of RAS inhibitors on diabetic retinopathy: a systematic review and meta-analysis. Lancet Diabetes Endocrinol 2015; 3: 263–274.

37.

Standl

Erbach

Schnell

What should be the antihypertensive drug of choice in diabetic patients and should we avoid drugs that increase glucose levels? Pro and cons. Diabetes Metab Res Rev 2012; 28: 60–66.

38.

Yamaji

Tsutamoto

Kawahara

, et al. Effect of eplerenone versus spironolactone on cortisol and hemoglobin A(1)(c) levels in patients with chronic heart failure. Am Heart J 2010; 160: 915–921.

39.

Effects of aspirin treatment on diabetic retinopathy: ETDRS report number 8 – Early Treatment Diabetic Retinopathy Study Research Group. Ophthalmology 1991; 98: 757–765.

40.

Chew

Davis

Danis

, et al. The effects of medical management on the progression of diabetic retinopathy in persons with type 2 diabetes: the Action to Control Cardiovascular Risk in Diabetes (ACCORD) Eye Study. Ophthalmology 2014; 121: 2443–2451.

41.

Mozetic

Pacheco

Latorraca

COC

, et al. Statins and/or fibrates for diabetic retinopathy: a systematic review and meta-analysis. Diabetol Metab Syndr 2019; 11: 92.

Supplementary Material

Please find the following supplemental material available below.

For Open Access articles published under a Creative Commons License, all supplemental material carries the same license as the article it is associated with.

For non-Open Access articles published, all supplemental material carries a non-exclusive license, and permission requests for re-use of supplemental material or any part of supplemental material shall be sent directly to the copyright owner as specified in the copyright notice associated with the article.

0.00 MB

0.06 MB