Figure. Kaplan-Meier time to diabetic macular edema (DME) curves according to thiazolidinedione use with or without insulin. The log-rank test gives an χ2 statistic of 373 (P < .001), which shows a clear difference in DME incidence according to thiazolidinedione use. In a comparison of thiazolidinedione use with nonuse, the hazard ratio was 5.19 (95% CI, 4.31-6.25).
Idris I, Warren G, Donnelly R. Association Between Thiazolidinedione Treatment and Risk of Macular Edema Among Patients With Type 2 Diabetes. Arch Intern Med. 2012;172(13):1005-1011. doi:10.1001/archinternmed.2012.1938
Author Affiliations: Sherwood Forest Hospitals Foundation Trust, Nottinghamshire, England (Dr Idris); Trent Research Design Service for the East Midlands, Nottingham, England (Dr Warren); and School of Graduate Entry Medicine and Health, University of Nottingham, Nottingham (Drs Idris and Donnelly).
Background Findings of prior studies have been inconclusive about the ocular effects of thiazolidinediones on diabetic macular edema (DME). We evaluated, in patients with type 2 diabetes (T2D), the short-term and long-term risks of developing DME among users vs nonusers of thiazolidinediones.
Methods A retrospective cohort study of 103 368 patients with T2D and no DME at baseline using The Health Improvement Network (THIN) database. Clinical, biochemical, and demographic information was obtained for the period January 1, 2000, through November 30, 2009.
Results At 1 year, the incidence of DME was 1.3% (n = 41) and 0.2% (n = 227) among thiazolidinedione users (n = 3227) and nonusers (n = 100 141), respectively (odds ratio [OR], 5.7 [95% CI, 4.1-7.9]). After Cox multiple regression analysis (adjusted for age; systolic blood pressure; levels of lipids and hemoglobin A1c; and use of aspirin, fibrates, insulin, oral antidiabetic drugs, or renin-angiotensin system blockers), multiple imputation analysis to adjust for missing values, and propensity score analysis to exclude for any selection bias, thiazolidinedione use was associated with an increased risk of DME at 1-year follow-up (OR, 2.3 [95% CI, 1.5-3.6]) and 10-year follow-up (hazard ratio [HR], 2.3; [95% CI, 1.7-3.0]). The effect was similar for pioglitazone and rosiglitazone. Combination therapy with insulin plus a thiazolidinedione was associated with a higher risk of DME after propensity score adjustment (HR, 3.0 [95% CI, 1.5-5.9]), while aspirin use (HR, 0.6 [95% CI, 0.4-0.9]) and angiotensin-converting enzyme inhibitor use (HR, 0.4 [95% CI, 0.2-0.7]) were associated with a reduced risk of DME.
Conclusion Among patients with T2D, treatment with a thiazolidinedione was associated with an increased risk of DME at 1-year and 10-year follow-up evaluations.
The thiazolidinediones, pioglitazone and rosiglitazone, are peroxisome proliferator–activated receptor γ (PPAR-γ) agonists that ameliorate peripheral and hepatic insulin resistance and are effective glucose-lowering treatments for patients with type 2 diabetes mellitus.1 The American Diabetes Association (ADA) and the European Association for the Study of Diabetes (EASD) recommend thiazolidinediones as second- or third-line therapy in combination with other oral agents or insulin to achieve target levels of glycemic control.2 The risk-benefit ratio for thiazolidinedione use has been the subject of intense discussion following a series of metabolic and cardiovascular outcome studies.3,4 Several analyses have highlighted the major adverse effects of thiazolidinediones, including an increased incidence of bone fractures,5 fluid retention and edema,6 increased risk of heart failure,7 and, more recently, potentially increased risk of bladder cancer.8
The PPAR-γ receptors are expressed abundantly in endothelial, renal, and retinal vascular tissues.9,10 Thiazolidinedione-induced peripheral and pulmonary edema is well recognized, and the underlying mechanisms are likely to include vasodilation, renal sodium reabsorption, and direct effects on vascular endothelial permeability.11- 13 More recently, however, small clinical studies have suggested an association between thiazolidinediones and diabetic macular edema (DME),14 one of the major sight-threatening complications that affects up to 20% of patients with type 2 diabetes.15
A cohort study from the Kaiser Permanente database reported an increased incidence of DME following short-term (1-year) exposure to pioglitazone,16 but no association between thiazolidinediones and DME was observed in a cross-sectional analysis of baseline data in the ACCORD eye substudy,17 and neither of these studies investigated the potential longer-term effects of thiazolidinedione use. Accordingly, the aim of this study was to evaluate the short- and long-term effects of thiazolidinedione treatment (pioglitazone and rosiglitazone) on the risks of developing DME in a large population cohort and to identify any risk factors that may influence visual outcome in patients with type 2 diabetes treated with a thiazolidinedione.
The study was approved by the National Research Ethics Committee East of England–Cambridge South (reference No. 10/H0305/46). This was a retrospective cohort study using The Health Improvement Network (THIN) database. THIN is an established database that collects anonymized electronic data from a volunteer sample of United Kingdom general practices using the Vision primary care computer system (In Practice Systems, London, England). THIN includes medical records on 9.1 million patients (3.4 million of whom are alive and registered with a practice currently contributing to THIN at the date of last collection) attending approximately 479 general practices in England and Wales. Data recorded in THIN include demographic information, medical diagnoses, prescribed medications, laboratory results, lifestyle characteristics, and referrals to specialists. Past and current medical diagnoses are recorded using READ codes (a thesaurus of coded medical terms maintained and distributed by the United Kingdom Terminology Center). THIN has been shown to be highly representative of the wider United Kingdom population. It has been well validated at both the practice and data set level by comparison with other national data sources for demographics and morbidity, mortality, prevalence, and geographic rates.18,19
Patients selected for analysis had a diagnosis of type 2 diabetes (defined as diabetes not requiring insulin therapy within 6 months of diagnosis), were 18 years or older, and were registered with the medical practices for at least 1 year from the index date of January 2000. A diagnosis of diabetes was based on READ codes for “diabetes” diagnosis and/or use of glucose-lowering therapies. The index date was chosen based on the licensing of both rosiglitazone and pioglitazone in the United Kingdom. The last date of follow-up for this study was November 30, 2009. According to figures from the THIN database, 95% of repeat prescriptions in general practice have a periodicity (frequency) of 6 months or less. Patients with less than 6 months of exposure to a thiazolidinedione were excluded from the study. Individuals were included in the analysis if they had no history of DME before baseline. The patient cohort was identified using relational database tools that involve strategies to design, code, test, and run programs to identify patients' information and outcome parameters.
Medical and demographic codes for the parameters of interest included age; body mass index (BMI; calculated as weight in kilograms divided by height in meters squared); body weight; levels of hemoglobin A1c (HbA1c) and serum lipids (total cholesterol, low-density lipoprotein cholesterol [LDL-C], high-density lipoprotein cholesterol [HDL-C], and triglycerides); use of specific antihypertensive drugs, lipid-lowering therapies, insulin, and aspirin; and systolic and diastolic blood pressure (BP).
The primary outcome measure of the study was a new diagnosis of DME as indicated by READ codes (including all codes relating to a diagnosis of DME and/or focal laser photocoagulation for DME) at 1 year and 10 years from the patient's entry into the study. A diagnosis of DME recorded on the THIN database would typically originate from a specialist ophthalmologic assessment of the patient following routine annual eye screening. Specialist-derived diagnoses of DME on THIN have been validated for accuracy,20 but some degree of measurement error is still possible.
Comprehensive diabetic retinopathy screening in the United Kingdom is free to patients and is based on a standardized protocol, independent of drug prescription, comorbidities, and diabetes consultations. National screening services aim to routinely screen all patients with diabetes yearly for the presence of diabetic retinopathy, including DME. Beginning in 2005, all retinopathy screening performed in England and Wales has been based on digital photographic images, in line with recommendations from the National Diabetic Retinopathy Screening program. Patients with suspected DME at screening are routinely referred for specialist ophthalmologic assessment at the Hospital Eye Service. The coded diagnosis of DME, or a coded diagnosis of focal laser photocoagulation (a specific therapy for DME), is then recorded in the THIN database. Information on visual acuity is unfortunately not available in THIN.
The THIN database was queried to determine the use of prescribed thiazolidinediones prior to the diagnosis of DME or prior to censuring of patients with no diagnosis of DME. Information on drug prescription was extracted from the database.
Baseline clinical and demographic variables were compared between users and nonusers of thiazolidinediones using the t test or the Wilcoxon rank-sum test for continuous variables, and the χ2 test or the 2-tailed Fisher exact test for categorical variables.
A 2 × 2 contingency table was used to calculate the unadjusted odds ratio (OR) of DME incidence at 1 year with rosiglitazone and pioglitazone exposure. Multiple logistic regression analysis was used as the primary analysis to calculate the adjusted OR controlling for the following covariates: systolic and diastolic BP, HbA1c and lipid levels; body weight; BMI; and the use of insulin, oral antidiabetic therapies, ACE inhibitors, angiotensin II receptor blockers, lipid-lowering therapies, and/or aspirin. All values were measured at baseline. The possible confounding variables were prechosen on clinical grounds, including the thiazolidinedione-insulin interaction, and variables to be included in the final model were chosen by comparing the log-likelihoods of competing models. A diagnostic plot (log[time] vs log[−log(survival function)]) confirmed that the constant proportionality assumptions were met. In addition, the Hosmer and Lemeshow parameter did not show evidence for lack of fit (P = .44).
Kaplan-Meier curves were constructed using 10-year outcome data to illustrate graphically the difference in DME incidence according to exposure to a thiazolidinedione. A log-rank test was used to check for statistically significant differences between users and nonusers of thiazolidinediones for the 10-year follow-up and subgroup analysis for 6-month follow-up.
Cox regression analysis was used to calculate DME incidence hazard ratios (HRs) adjusted for confounding variables. All possible confounders, including the thiazolidinedione-insulin interaction and aspirin and/or ACE-inhibitor use, were preplanned, and the model was selected by comparing model log-likelihoods. Threshold significance was given at the conventional level of P = .05, and 95% CIs are shown for individual HRs.
To address missing values in the Cox regression analysis, multiple imputation was performed using a linear regression method to impute 10 HbA1c values per missing case. Multiple imputation was not used in the final logistic regression model because HbA1c level, the variable with a significant number of missing values, was not included in the final model.
In the United Kingdom Prospective Diabetes Study, HbA1c level was shown to increase with diabetes duration, regardless of intensive or conventional glucose control.21 We therefore performed a sensitivity analysis for both the multiple logistic and Cox regression results by performing propensity score analyses,22 as determined by HbA1c level and other covariates to adjust for treatment selection bias and to ensure that both cohorts were comparable. A nonparsimonious logistic regression analysis was used to generate a propensity score, which is a predictor for thiazolidinedione prescription based on the observed covariates (sex, Townsend index, age, BMI, systolic and diastolic BP, HbA1c level, and aspirin and/or insulin use). Previous studies have shown that stratification by propensity score into 5 strata is sufficient to remove 90% to 95% of the selection bias.23 An alternative method of kernel matching in which the outcome for each exposed individual was compared with a weighted mean of the outcomes of unexposed individuals was used to check for the sensitivity of the method of using the propensity scores to reduce bias.24
STATA 11 software (StataCorp LP) was used for the statistical analysis and IBM Modeler 13 for the database work. Results are expressed as ORs for logistic regression and HRs for Cox regression analysis.
The analysis identified 109 295 patients in 464 practices who fulfilled the initial selection criteria. Of these, 5927 patients were not included in the analysis because the baseline date for these records was earlier than the acceptable mortality reporting date for the practice, which made them potentially less reliable. Thus, the cohort for further analysis included 103 368 patients with type 2 diabetes. The baseline characteristics for users and nonusers of thiazolidinediones are listed in Table 1.
The incidence of DME at 1 year of follow-up according to thiazolidinedione use is listed in Table 2. Patients exposed to a thiazolidinedione had a significantly higher risk of developing DME (OR, 5.7 [95% CI, 4.1-7.9]) (Fisher exact test P < .001).
Following multiple logistic regression to correct for confounding factors and multiple imputation analysis to correct for missing values, we found that exposure to a thiazolidinedione was associated with a significantly increased risk of DME at 1 year (OR, 3.3 [95% CI, 2.2-5.0]). The analysis for pioglitazone and rosiglitazone individually showed that each drug was separately associated with a significantly increased risk of DME at 1 year, and there was no significant difference between the 2 thiazolidinediones (pioglitazone OR, 3.6 [95% CI, 2.0-6.6]; rosiglitazone OR, 3.1 [95% CI, 1.9-5.1]).
The increased risk of DME at 1 year was still present when propensity score methods used to remove selection bias (stratification OR, 2.3 [95% CI, 1.5-3.6]; kernel-matching OR, 2.7 [95% CI, 1.6-4.4]). The potential confounders that did not change the magnitude of associations of interest and were therefore excluded were age, HbA1c, BMI, lipid measurements, and use of aspirin, fibrates or angiotensin receptor blockers (Table 3).
Ten years of follow-up data for the patient cohort (n = 103 359) were used to construct Kaplan-Meier curves for time to DME incidence and were further analyzed by Cox regression analysis. The log-rank test yielded an χ2 statistic of 373 (P < .001) and showed a clear difference in DME incidence between users and nonusers of thiazolidinediones (HR, 5.2 [95% CI, 4.3-6.3]) (P < .001) (Figure). A log-rank test was also performed for time to DME up to 6 months with P < .001, which shows more clearly that the increased risk was also present at short-term follow-up.
Following adjustment for confounding factors using Cox regression and multiple imputation analysis to adjust for missing values, we found that there was a significantly increased risk of developing DME among patients exposed to a thiazolidinedione (HR, 2.8 [95% CI, 2.1-3.5]) (P < .001); after adjusting for selection bias using propensity score stratification, we found a similar result (HR, 2.3 [95% CI, 1.7-3.0]). In this multiple regression model, HbA1c, systolic BP, aspirin use, ACE-inhibitor use, and serum triglyceride levels were variables that produced substantial changes in the magnitudes of the associations of interest (Table 4). A plot to test for the proportional hazard assumption showed that the observed and predicted curves agreed and only deviated when the number of DME events in the thiazolidinedione group small at longer-term follow-up.
Further Cox regression and multiple imputation analysis using an interaction model showed that combination therapy with insulin and a thiazolidinedione increased the HR of DME even further (HR, 4.4 [95% CI, 2.5-7.8]; or HR with propensity score adjustment, 3.0 [95% CI, 1.5-5.9]). Conversely, concurrent use of aspirin (HR, 0.7 [95% CI, 0.5-0.9]) and an ACE inhibitor (HR, 0.5 [95% CI, 0.3-0.9]) was associated with a reduction in the HR of developing DME independent of other covariates and after propensity score stratification: HR, 0.6 (95% CI, 0.4-0.9) and HR, 0.4 (95% CI, 0.2-0.7), respectively (Table 4).
To further strengthen the analysis, propensity scores were derived to exclude indication bias and produce comparable subgroups for a stratified analysis. The potential confounding variables for users and nonusers of thiazolidinediones were generally balanced across the propensity score quintiles (Table 5).
This large retrospective cohort study analyzed the primary care electronic medical records of more than 100 000 patients with type 2 diabetes and showed that, even after adjustment for various confounding factors known to influence diabetic retinopathy, exposure to a thiazolidinedione is associated with an increased risk of developing DME. The association was evident with both pioglitazone and rosiglitazone. A previous short-term study from Kaiser Permanente reported a similar increase in risk of DME after 1 year (OR, 2.6),16 but the analysis did not adjust for several key variables, in particular the use of BP-lowering drugs such as angiotensin receptor blockers and ACE inhibitors, which may affect the progression of diabetic retinopathy.25,26
Diabetic macular edema is a sight-threatening chronic condition,14,15 so an important new observation in the present study is that the increased risk of DME continued to accrue during the 10-year of follow-up. The patients at greatest risk of developing DME were those taking thiazolidinediones in combination with insulin. Diabetic macular edema has been associated in clinical practice with those patients who develop thiazolidinedione-induced pedal edema,14,27 and in large clinical trials the highest incidence of fluid retention and peripheral edema (16%) was observed in the subgroup of patients undergoing combination therapy with a thiazolidinedione and insulin.28
Although glycemic control did not modulate the short-term (1-year) risk of DME among users of thiazolidinediones, HbA1c was significant in the longer-term analysis. In addition, the use of aspirin and/or ACE inhibitor drugs was associated with a modest reduction in the HR of DME, which is consistent with previous clinical and experimental studies.25
Diabetic macular edema is the leading cause of blindness among patients with type 2 diabetes.15,29 The absolute rate of new-onset DME was much lower in the present study (about 1%) than has been reported in previous clinical trials of diabetic retinopathy.29 This probably reflects the low detection rates in routine clinical practice and highlights the importance of a national systematic screening program for diabetic retinopathy. The cause of DME is complex, but disruption to the blood-retinal barrier (BRB), mediated in part by local release of vascular endothelial growth factor (VEGF),30 results in chronic leakage of fluid and reduced visual acuity.15 Important regulators of VEGF include hypoxia, hyperglycemia, and angiotensin II, but PPAR-γ agonists also increase VEGF expression,31 and our research group12 has previously shown that thiazolidinediones increase vascular endothelial permeability in vitro. A number of systemic changes associated with thiazolidinediones may also contribute (eg, sodium and fluid retention, vasodilation, changes in BP, microvascular perfusion).32 Thus, there are plausible mechanisms, both local and systemic, by which thiazolidinediones might adversely affect the risk of DME, particularly in insulin-treated patients.
The strengths of this cohort analysis include a very large sample size and a long duration of follow-up. To increase the reliability of the study, we included only patients who were apparently compliant with thiazolidinedione therapy (>6 months and in receipt of repeat prescriptions), and we measured new-onset DME during a 10-year observation period. Diabetic macular edema is a specialist diagnosis that would normally be made only by an ophthalmologist when the patient is referred following routine annual eye screening.
It is an important limitation of this study that information about duration of diabetes and duration of individual patient exposure to a thiazolidinedione was not available. However, the baseline characteristics of patients in the 2 cohorts were comparable in terms of age, sex, and BMI. The percentage of patients receiving insulin and mean HbA1c levels were slightly higher in the thiazolidinedione user group, but this was balanced by slightly higher BP and more patients receiving aspirin, ACE inhibitors, and statins in the nonuser group. Furthermore, the propensity score analysis mitigates against indication bias and strengthens the interpretation of a true association between DME and use of a thiazolidinedione.
During the study period (2000-2009), thiazolidinediones were routinely and widely used in clinical practice in patients with diabetes of varying duration (eg, as dual or initial combination therapy with metformin in those with short-duration disease and as third-line agents added to metformin/sulfonylurea combination therapy).
Despite propensity score analysis and the similar baseline profiles, it is possible, but in our view unlikely, that the thiazolidinedione users were significantly different than the nonusers either in diabetes duration or in their cardiovascular risk profile. The present study confirms and extends other observations implicating thiazolidinedione use with an increased risk of DME, but the data require cautious interpretation.
In conclusion, this large-population study has shown that, even after adjustment for a range of confounding variables, thiazolidinediones are associated with an increased risk of DME in patients with type 2 diabetes, especially those undergoing insulin therapy. The risk increased continuously over the 10-year duration of this study. The explanation for this association is unclear, but plausible mechanisms described previously exist by which these drugs might adversely influence this important ocular complication of type 2 diabetes. A larger and more detailed meta-analysis of randomized controlled trials (ideally in high-risk patients) will be needed to clearly establish the risk-benefit profile of thiazolidinediones in patients with, or at risk of, DME. In addition, future clinical trials to evaluate new PPAR-γ (or dual α,γ) agonists should prospectively evaluate the ocular safety of these drugs. Clinicians should be vigilant in the clinical screening for DME among those patients taking thiazolidinediones.
Correspondence: Iskandar Idris, MD, FRCP, FRCP(Edin), Department of Diabetes and Endocrinology, Kings Treatment Centre, Clinic 5, Sherwood Forest Hospitals Foundation Trust, Nottinghamshire, England (email@example.com).
Accepted for Publication: March 13, 2012.
Published Online: June 11, 2012. doi:10.1001/archinternmed.2012.1938
Author Contributions: Dr Idris had full access to all the data in the study and takes responsibility for the integrity of the data and the accuracy of the data analysis. Study concept and design: Idris and Donnelly. Acquisition of data: Idris. Analysis and interpretation of data: Idris and Warren. Drafting of the manuscript: Idris and Donnelly. Critical revision of the manuscript for important intellectual content: Idris, Warren, and Donnelly. Statistical analysis: Idris and Warren. Obtained funding: Idris. Administrative, technical, and material support: Donnelly. Study supervision: Donnelly.
Financial Disclosure: None reported.
Funding/Support: Sherwood Forest Hospitals Charitable Trust Fund (no external funding).
Role of the Sponsor: The sponsor's role was only to provide indemnity arrangement.
Additional Contributions: We thank Mary Thomson, MD, research manager at THIN, for providing access to the database.