Statins use and risk of new-onset diabetes in hypertensive patients: a population-based retrospective cohort study in Yinzhou district, Ningbo city, People's Republic of China.

BACKGROUND: Reports have suggested that statin use is associated with an increased incidence of type 2 diabetes mellitus (T2DM). Guidelines suggested that statins should be prescribed in hypertensive patients for primary prevention. However, there were very few studies on the risk of T2DM associated with statin use among patients with hypertension in mainland People's Republic of China.
PURPOSE: To determine the association between statin use and new-onset diabetes mellitus among patients with hypertension in mainland People's Republic of China. PATIENTS AND METHODS: We performed a retrospective cohort study of hypertensive patients using the Yinzhou regional health care database from January 1, 2010, to August 31, 2016. Patients aged 30-90 years old without T2DM were eligible for inclusion. We identified new statin initiators and nonusers by using prescription records of inpatients and outpatients. Multivariate Cox model and propensity score methods were used to adjust potential confounders, including age, sex, body mass index, comorbidities, lifestyle characteristics, and baseline antihypertensive drug use. The risk of incident T2DM among statin initiators compared to nonusers was estimated by the Cox proportional hazards model. Propensity scores for statin use were then developed using logistic regression, statin initiators were matched 1:1 with nonusers according to propensity scores with the nearest neighbor matching method within 0.2 caliper width, and Cox regression was again conducted.
RESULTS: Among 67,993 patients (21,551 statin initiators; 46,442 nonusers), the unadjusted incidence rate of incident T2DM was higher in statin initiators than nonusers (25.68 versus 14.19 events/1,000 person-years; adjusted hazard ratio: 1.55; 95% confidence interval: 1.44-1.66). After propensity score 1:1 matching (19,818 statin initiators; 19,818 nonusers), baseline characteristics between 2 groups were balanced except that the nonusers group was 0.53 years older on average (P<0.001). Then statin use was still associated with a significant increased risk for T2DM in the matched cohort (adjusted hazard ratio: 1.54; 95% confidence interval: 1.41-1.67). Subgroup analyses also demonstrated similar findings.
CONCLUSION: Our study indicated an association between statin use and an increased risk of new-onset diabetes mellitus. It provides better understanding of statin and new-onset diabetes mellitus association among hypertensive patients in real-word setting. As an observational study, our findings were prone to unmeasured confounding and bias.

Introduction

There is strong clinical evidence that statins can effectively reduce all-cause mortality and cardiovascular events in patients at risk of cardiovascular diseases (CVDs).1–3 Thus, statins are widely used for primary or secondary prevention of CVD. However, studies of statin therapy have provided conflicting results about the association between statins use and risk for incident new-onset diabetes mellitus (NOD). The first study (West of Scotland Coronary Prevention Study) that assessed this association reported a protective effect (hazard ratio [HR]: 0.70; 95% CI: 0.50–0.98) but used nonstandardized criteria for diabetes diagnosis. However, results from subsequent hypothesis-testing trials have been inconsistent. Thus, several meta-analyses have been performed to evaluate this association. Rajpathak et al4 have conducted a meta-analysis of statins compared with placebo and included 57,593 patients from 6 trials and found an increased risk for incident diabetes mellitus (relative risk: 1.13; 95% CI: 1.03–1.23). Sattar et al5 reported a meta-analysis which showed statins were associated with a 9% increased risk of diabetes in 91,140 patients from 13 trials.

Based on the evidence from randomized controlled trials,6–9 systematic reviews,4,5 and a few observational studies,10,11 the United States Food and Drug Administration mandated labeling changes for all statins in February 2012. Labels now include information concerning diabetogenic effects, including inducing diabetes and increases in hemoglobin A1c or fasting glucose level. In November 2012, the China Food and Drug Administration warned that statins users may have an increased risk of raised blood sugar levels and the development of type 2 diabetes mellitus (T2DM). Although it is based mainly on evidence from clinical trials and meta-analyses, there are limited data regarding the long-term effects of these associations in mainland People’s Republic of China. Besides, several reports suggested that Chinese are more sensitive to statins at a certain dose.12–14 The guidelines on treatment and health care system are different from other countries, even Taiwan, in mainland People’s Republic of China.15,16 On the other hand, hypertensive patients are at increased risk of developing diabetes mellitus.17,18 And both hypertension and diabetes are the major causes of CVD.19,20 Evidence, expert consensus, and guidelines supported that statins should be prescribed in hypertensive patients for primary prevention.21–23 However, few studies have examined whether statin use was associated with an increased incidence of T2DM in hypertensive patients. Accordingly, the objective of this study was to assess the association between statin use and NOD among patients with hypertension in mainland china.

Patients and methods

Study design

We performed a population-based retrospective cohort study among patients aged 30–90 years old with hypertension in Yinzhou district, Ningbo city, People’s Republic of China. Hypertension was defined based on the diagnosis records of outpatients or inpatients (10th revision of the International Classification of Diseases [ICD-10]: I10). Patients who did not receive statins at any time during the study were attributed to the nonusers group. Statin use was defined as statin therapy for at least 90 days during the study period (January 1, 2011, to August 31, 2016), and the date of first statin prescription was termed the index date. To avoid immortal time bias, the index date for nonusers was assigned according to the distribution of users’ index date.24

We defined our diabetes endpoint as the first diagnosis of T2DM (ICD-10: E11) or a prescription for insulin or an oral antidiabetic medication as based on the data in the outpatients or inpatients records. The primary endpoint of the study was defined as the diabetes endpoint after 180 days of statin treatment during the follow-up period. We assumed that NOD within the 180 days after statin use could not be attributed to treatment. Thus, patients who had a diagnosis of diabetes within the first 180 days were censored in our time-to-event analysis. A sensitivity analysis was conducted in which all NOD cases after initiation of statin therapy are evaluated. In order to better identify the new statins users and NOD, patients were excluded from the study if they were diagnosed as having T2DM or prescribed statins before the index date.

Our study setup a look-back period (1 year preceding their index date) to define some important confounders, namely lifestyle, comorbidities, and antihypertensive drug use.

Data resources

Yinzhou is the largest district in Ningbo, an eastern coastal city of People’s Republic of China, with a population of nearly 1.2 million residents. In 2005, Yinzhou Health and Family Planning Commission intended to imitate the primary care model seen in UK. They setup 287 village- or community-based health service stations. Each station has more than 2 doctors who are responsible for the primary care and chronic diseases follow-up such as hypertension, diabetes, CVD, and cancers. In the meantime, they initiated a Health Information System, which collects the health care information of the residents. By 2010, this system covered all the hospitals (3 general hospitals, 24 community health centers, and 287 community health service stations) in Yinzhou. Moreover, it integrated hospital data, public health data, and community health data which were stored in the administrative health care database and included outpatient and inpatient records, prescribing information, as well as demographic characteristics. The outpatient and inpatient records comprised patient ID, date of diagnosis, diagnosis name, and diagnosis code (ICD-10). Prescribing information consisted of product name, brand, generic names, dosage, and prescription date. The lifestyle characteristics were recorded using a standard questionnaire by trained staff. Before obtaining the data set, all the data from patients was anonymized and no personal details were shared with health researchers for privacy protection. This database has been used and described in a previous study.25–27 This study was approved by the Peking University Health Science Center Ethics Committee.

Statistical analysis

In baseline analysis, age and BMI are presented as mean (standard deviation) and median (interquartile range). Differences in age and BMI were analyzed by using the Student’s t-test. Other variables are set as categorical or discrete, which were expressed as counts and percentages and evaluated by χ2 or Fisher’s exact test between the groups as appropriate.

In the primary analysis, Cox proportional hazard regression model was used to estimate the relation between statin use and risk of NOD. We developed a multivariable model to adjust the potential confounders, including age, sex, body mass index (BMI), comorbidities (dyslipidemia, myocardial infarction, coronary artery disease, cerebrovascular accidents, heart failure, and angina pectoris), lifestyle characteristics (drinking, smoking, physical activity), and baseline antihypertensive drug use (β-blockers, diuretics, angiotensin-converting enzyme inhibitors, angiotensin receptor blockers, calcium channel blockers).

Additionally, propensity score matching (PSM) approach was applied to balance the systematic difference between 2 groups as secondary analyses by using R package “MatchIt.” First, a propensity score, representing the probability of being in the statin-user or non-statin user group, was calculated for each patient via a logistic regression model. Variables were selected for inclusion in the model if they were associated with both the exposure and the outcome, or associated with the outcome. Age, sex, BMI, comorbidities, lifestyles, and baseline antihypertensive drug use were finally included. Second, statin users were then 1:1 matched to nonusers according to propensity scores with the nearest neighbor matching method within 0.2 caliper width. Third, Cox model was applied again and those factors remaining imbalanced after matching would be included in the model (eg, age).

To assess the possible effect modification of variables related to both the exposure and the outcome, we carried out subgroup analyses by stratifying the study cohort according to sex (female or male), age groups (≤49, 50–59, 60–69, and ≥70), BMI categories (<24, 24–28 and ≥28 kg/m2), smoking status (never, past, current), alcohol use (never, occasional, often), physical activity (<1 d/wk, 1–4 d/wk, >4 d/wk), individual statins (fluvastatin, rosuvastatin, simvastatin, and atorvastatin), and potency of statins (low, moderate, and high). We defined the potency based on recommendations from the 2013 American College of Cardiology/American Heart Association guidelines (high dose: atorvastatin ≥40 mg and rosuvastatin ≥20 mg; moderate dose: atorvastatin 10–20 mg, rosuvastatin 5–10 mg, simvastatin 20–40 mg, and fluvastatin 40–80 mg; and low dose: simvastatin 10 mg and fluvastatin 20–40 mg). The HRs were estimated from Cox models, adjusted for age, sex, BMI, comorbidities, lifestyle, and baseline antihypertensive drug use. Age was excluded in models within age subgroups. Similarly, corresponding variables were individually excluded in models within corresponding subgroups. All tests were 2 sided, and an α level of 0.05 was considered statistically significant. All time-to-event analyses were performed with SAS, version 9.4 (SAS Institute Inc., Cary, NC, USA).

Results

Overall, our study identified 95,327 patients who were aged 30–90 and diagnosed as hypertension in 2010 by using the Yinzhou database (Figure 1). A total of 67,993 patients were finally enrolled after restricting the analysis to patients without missing data, newly exposed to statins, and without history of diabetes (21,551 statin users and 46,442 nonusers).

Figure 1

Flow chart of the study population.

Table 1 describes the baseline characteristics of statin users and nonusers before and after PSM. In the original cohort, the statins users group and nonusers group had a significant difference in all the baseline characteristics. After PSM, baseline characteristics of two matched groups (19,819 statin users and 19,819 nonusers) were balanced, except that the nonusers group was a bit older. The mean age of nonusers group was 61.93 years, significantly older than statin group (61.40 years). Figure 2 shows the histograms of the density of propensity scores for original cohort and matched cohort. During the follow-up period 4,172 (6.14%) patients developed NOD in the total cohort. According to our definition of outcome, there were 3,421 events including 1,884 events in statin group and 1,537 events among nonusers. Patients who developed NOD were older, were more overweight, and were more likely to have dyslipidemia (Table 2).

Table 1

Baseline characteristics of the original and matched cohorts of hypertensive patients without history of diabetes

Characteristics	Overall patients			Matched patients
	Non-statin	Statin	P-value	Non-statin	Statin	P-value
N (%)	46,442 (100)	21,551 (100)		19,819 (100)	19,819 (100)
Age
Mean ± SD	58.87 (13.0)	61.48 (11.2)	<0.001	61.93 (12.8)	61.40 (11.2)	<0.001
Median (Q1–Q3)	58.00 (49–67)	61.00 (54–69)		61.00 (53–72)	61.00 (54–69)
Sex (female)	24,487 (52.7)	11,989 (55.6)	<0.001	11,065 (55.8)	10,990 (55.5)	0.992
Dyslipidemia	1,085 (2.3)	2,734 (12.7)	<0.001	1,084 (5.5)	1,099 (5.5)	0.982
MI	13 (0.0)	28 (0.1)	<0.001	11 (0.1)	8 (0.0)	0.522
CAD	1,353 (2.9)	1,210 (5.6)	<0.001	582 (2.9)	630 (3.2)	0.385
CVA	428 (0.9)	518 (2.4)	<0.001	247 (1.2)	253 (1.3)	0.818
Heart failure	437 (0.9)	418 (1.9)	<0.001	215 (1.1)	227 (1.1)	0.995
Angina pectoris	407 (0.9)	489 (2.3)	<0.001	227 (1.1)	219 (1.1)	0.257
BMI, kg/m2
Mean ± SD	23.34 (4.26)	23.85 (5.92)	<0.001	23.71 (5.28)	23.79 (4.00)	0.090
Median (Q1–Q3)	23.10 (21.50–24.90)	23.50 (21.80–25.50)		23.50 (21.60–25.40)	23.50 (21.80–25.50)
BMI, kg/m2			<0.001			0.811
<24	29,305 (63.1)	12,029 (55.8)		11,071 (55.9)	11,150 (56.3)
24–28	14,345 (30.9)	7,642 (35.5)		7,086 (35.8)	6,989 (35.3)
≥28	2,792 (6.0)	1,880 (8.7)		1,662 (8.4)	1,680 (8.5)
Cigarette smoking			<0.001			0.421
Never	34,918 (75.2)	16,302 (75.6)		15,026 (75.8)	15,003 (75.7)
Past	4,156 (8.9)	2,268 (10.5)		2,129 (10.7)	2,068 (10.4)
Current	7,368 (15.9)	2,981 (13.8)		2,664 (13.4)	2,748 (13.9)
Alcohol consumption			0.024			0.952
Never	36,328 (78.2)	16,945 (78.6)		15,576 (78.6)	15,586 (78.6)
Occasional	2,723 (5.9)	1,183 (5.5)		1,064 (5.4)	1,096 (5.5)
Often	7,391 (15.9)	3,423 (15.9)		3,179 (16.0)	3,137 (15.8)
Physical activity			<0.001			0.682
<1 d/wk	20,359 (43.8)	9,796 (45.5)		9,012 (45.5)	8,973 (45.3)
1–4 d/wk	15,475 (33.3)	7,279 (33.8)		6,674 (33.7)	6,690 (33.8)
>4 d/wk	10,608 (22.8)	4,476 (20.8)		4,133 (20.9)	4,156 (21.0)
β blockers	5,995 (12.9)	3,684 (17.1)	<0.001	3,295 (16.6)	3,332 (16.8)	0.989
Diuretics	8,654 (18.6)	4,982 (23.1)	<0.001	4,577 (23.1)	4,523 (22.8)	0.659
ACEI	10,277 (22.1)	5,780 (26.8)	<0.001	5,310 (26.8)	5,244 (26.5)	0.407
ARBs	6,863 (14.8)	4,300 (20.0)	<0.001	3,900 (19.7)	3,934 (19.8)	0.881
CCBs	5,405 (11.6)	3,385 (15.7)	<0.001	3,093 (15.6)	3,098 (15.6)	0.978
Type of statin
Non-statin	46,442 (100.0)	0 (0.0)		19,819 (100.0)	0 (0.0)
Atorvastatin	0 (0.0)	9,263 (43.0)		0 (0.0)	8,516 (43.0)
Fluvastatin	0 (0.0)	170 (0.8)		0 (0.0)	157 (0.8)
Rosuvastatin	0 (0.0)	1,403 (6.5)		0 (0.0)	1,326 (6.7)
Simvastatin	0 (0.0)	10,715 (49.7)		0 (0.0)	9,820 (49.5)

Note: Unless indicated all values in this table are presented as n (%).

Abbreviations: MI, myocardial infarction; CAD, coronary artery disease; CVA, cerebrovascular accident; BMI, body mass index; d, day; wk, week; ACEI, angiotensin-converting enzyme inhibitors; ARBs, angiotensin receptor blockers; CCBs, calcium channel blockers.

Figure 2

Histograms of the density of propensity scores for cohort before and after matching.

Table 2

Baseline characteristics of patients who developed T2DM

Characteristics	Non-T2DM	T2DM	P-value
N (%)	63,821 (100)	4,172 (100)
Age
Mean ± SD	59.63 (12.63)	60.81 (10.76)	<0.001
Median (Q1–Q3)	59.00 (51–68)	61.00 (54–68)	<0.001
Sex (female)	34,124 (53.5)	2,352 (56.4)	<0.001
Dyslipidemia	3,346 (5.2)	473 (11.3)	<0.001
MI	36 (0.1)	5 (0.1)	0.196
CAD	2,274 (3.6)	289 (6.9)	<0.001
CVAs	834 (1.3)	112 (2.7)	<0.001
Heart failure	758 (1.2)	97 (2.3)	<0.001
Angina pectoris	811 (1.3)	85 (2.0)	<0.001
BMI, kg/m2
Mean ± SD	23.45 (4.93)	24.30 (3.27)	<0.001
Median (Q1–Q3)	23.20 (21.50–25.10)	24.00 (22.20–26.00)
BMI, kg/m2			<0.001
<24	39,283 (61.6)	2,051 (49.2)
24–28	4,173 (6.5)	499 (12.0)
≥28	20,365 (31.9)	1,622 (38.9)
Cigarette smoking			<0.001
Never	48,064 (75.3)	3,156 (75.6)
Past	5,946 (9.3)	478 (11.5)
Current	9,811 (15.4)	538 (12.9)
Alcohol consumption			0.359
Never	50,039 (78.4)	3,234 (77.5)
Occasional	3,651 (5.7)	255 (6.1)
Often	10,131 (15.9)	683 (16.4)
Physical activity			<0.001
<1 d/wk	28,137 (44.1)	2,018 (48.4)
1–4 d/wk	21,338 (33.4)	1,416 (33.9)
>4 d/wk	14,346 (22.5)	738 (17.7)
β blockers	9,101 (14.3)	578 (13.9)	0.481
Diuretics	12,727 (19.9)	909 (21.8)	0.004
ACEI	14,959 (23.4)	1,098 (26.3)	<0.001
ARBs	10,344 (16.2)	819 (19.6)	<0.001
CCBs	8,163 (12.8)	627 (15.0)	<0.001
Type of statin			<0.001
Non-statin	44,226 (69.3)	2,216 (53.1)
Atorvastatin	8,553 (13.4)	710 (17.0)
Fluvastatin	153 (0.2)	17 (0.4)
Rosuvastatin	1,321 (2.1)	82 (2.0)
Simvastatin	9,568 (15.0)	1,147 (27.5)

Note: Unless indicated all values in this table are presented as n (%).

Abbreviations: T2DM, type 2 diabetes mellitus; MI, myocardial infarction; CAD, coronary artery disease; CVAs, cerebrovascular accidents; BMI, body mass index; d, day; wk, week; ACEI, angiotensin-converting enzyme inhibitors; ARBs, angiotensin receptor blockers; CCBs, calcium channel blockers.

Table 3 shows the results of primary analysis for 2 different definitions of endpoints. For NOD after 180 days of statin initiation, the crude incidence of NOD per 1,000 person-years by statin exposure status was 25.68 for statin group and 14.19 for nonusers. Compared to nonusers, the crude NOD HR was 1.82 (95% CI: 1.70–1.95). Adjustments were made in Cox model for age, sex, BMI, comorbidities, lifestyle, and antihypertensive drug. The adjusted NOD HR was 1.55 (95% CI: 1.44–1.66). After matching for propensity score and adjusting for imbalanced age, statin users showed a 54% increase in the risk for NOD with an adjusted HR of 1.54 (95% CI: 1.41–1.67). For all NOD after statin initiation, the crude incidence of NOD in both groups was found to be higher (16.69 for statin group, 32.68 for nonusers group). The crude NOD HR was 1.96 (95% CI: 1.84–2.08), and the adjusted NOD HR was 1.75 (95% CI: 1.64–1.87). In the matched cohort, the adjusted NOD HR was 1.65 (95% CI: 1.53–1.78). As can be seen from Figure 3, the cumulative incidence of NOD was higher among new users of statins compared with non-statin users during the follow-up period.

Table 3

Incidence rates of NOD by statin exposure status

Group	Number of patients	Number of outcome	Total person-years	Incidence outcome per 1,000 PYs	Adjusted HR (95%)
NOD after 180 days
Non	46,442	1,884	132,767.5	14.19	1.55 (1.44–1.66)
Statin	21,551	1,537	59,848.48	25.68
Non (matched)	19,818	917	56,445.33	16.25	1.54 (1.41–1.67)
Statin (matched)	19,818	1,334	55,237.63	24.15
All NOD
Non	46,442	2,216	132,767.5	16.69	1.75 (1.64–1.87)
Statin	21,551	1,956	59,848.48	32.68
Non (matched)	19,818	1,079	56,445.33	19.06	1.65 (1.53–1.78)
Statin (matched)	19,818	1,729	55,237.63	31.31

Abbreviations: NOD, now onset of diabetes mellitus; HR, hazard ratio.

Figure 3

Cumulative incidence of NOD between new users of statin and non-statin users during the follow-up period.

Abbreviations: NOD, new-onset of diabetes mellitus; grp, group.

We found consistent results in subgroup analysis assessing the association between statins use and increase incident of NOD (Table 4). Sex-specific analysis showed that the adjusted NOD HR was 1.58 (95% CI: 1.42–1.77) among male patients, which was similar to that among female patients (adjusted HR: 1.52, 95% CI: 1.38–1.67). In addition, age-specific HRs decreased in older people, even though the association in each stratum was statistically significant after adjustment for confounding factors. We observed significantly increased risk of NOD across subgroups by lifestyle (smoking, alcohol use, and physical activity). The statin-user group showed a significantly increased risk for NOD in patients with lower BMI level. Among different types of statin agents, fluvastatin had the lowest, and most statistically insignificant, risk of NOD (adjusted HR: 1.44, 95% CI: 0.83–2.48), while other three types of statins were associated with a significantly increased risk of NOD (rosuvastatin: adjusted HR: 1.47, 95% CI: 1.19–1.83; simvastatin: adjusted HR: 1.64, 95% CI: 1.51–1.78; atorvastatin: adjusted HR: 1.52, 95% CI: 1.37–1.68). Compared to nonusers, the adjusted HRs were 1.44 (95% CI: 1.29–1.61), 1.65 (95% CI: 1.53–1.79), and 1.81 (95% CI: 0.86–2.82) for the low-, moderate-, and high-intensity statin users, respectively.

Table 4

Risk for incident NOD by statin use, within sex, age, BMI, smoke status, alcohol use, exercise frequency, statin types, and potency subgroups

Subgroups	N (%)	Unadjusted HR (95% CI)	Adjusted HR (95% CI)		P-value
Sex
Male	31,517 (46.35)	1.86 (1.68–2.06)	1.58(1.42–1.77)		<0.0001
Female	36,476 (53.65)	1.78 (1.62–1.94)	1.52 (1.38–1.67)		<0.0001
Age
≤49	15,404 (22.66)	2.39 (2.00–2.84)	1.90 (1.58–2.29)		<0.0001
50–59	19,618 (28.85)	1.88 (1.66–2.12)	1.64 (1.45–1.87)		<0.0001
60–69	17,400 (25.59)	1.43 (1.27–1.60)	1.28 (1.13–1.45)		<0.0001
≥70	15,571 (22.90)	1.70 (1.47–1.96)	1.39 (1.19–1.62)		<0.0001
BMI, kg/m2
<24	42,371 (62.32)	1.92 (1.75–2.12)	1.65 (1.50–1.83)		<0.0001
24–28	21,245 (31.25)	1.59 (1.42–1.77)	1.43 (1.27–1.60)		<0.0001
≥28	4,377 (6.44)	1.53(1.26–1.86)	1.44 (1.17–1.77)		<0.0001
Cigarette smoking
Never	51,220 (75.33)	1.77 (1.64–1.92)	1.57 (1.44–1.70)		<0.0001
Past	6,424 (9.45)	1.67 (1.37–2.04)	1.56 (1.27–1.93)		<0.0001
Current	10,349 (15.22)	2.16 (1.79–2.60)	1.79 (1.46–2.18)		<0.0001
Alcohol consumption
Never	53,273 (78.35)	1.82 (1.68–1.96)	1.60 (1.48–1.73)		<0.0001
Occasional	3,906 (5.74)	2.05 (1.55–2.71)	1.70 (1.26–2.29)		<0.0001
Often	10,814 (15.90)	1.76 (1.49–2.09)	1.57 (1.31–1.87)		<0.0001
Physical activity
<1 d/wk	30,155 (44.35)	1.99 (1.70–2.35)	1.71 (1.44–2.03)		<0.0001
1–4 d/wk	22,754 (33.47)	1.83 (1.63–2.05)	1.66 (1.47–1.87)		<0.0001
>4 d/wk	15,084 (22.18)	1.73 (1.57–1.91)	1.54 (1.39–1.71)		<0.0001
Type of statin
Fluvastatin	170 (0.79)	1.57 (0.91–2.70)	1.44 (0.83–2.48)		0.4903
Rosuvastatin	1,403 (6.51)	1.50 (1.15–1.97)	1.47 (1.19–1.83)		0.0004
Simvastatin	10,715 (49.72)	1.85 (1.71–2.00)	1.64 (1.51–1.78)		<0.0001
Atorvastatin	9,263 (42.98)	1.80 (1.64–1.98)	1.52 (1.37–1.68)		<0.0001
Potency of statin
Low	7,450 (34.57)	1.74 (1.57–1.92)	1.44 (1.29–1.61)		0.0042
Moderate	13,976 (64.85)	1.86 (1.72–2.00)	1.65 (1.53–1.79)		<0.0001
High	125 (0.58)	1.00 (0.32–3.10)	1.80 (0.86–2.82)		0.1432

Abbreviations: NOD, new-onset diabetes mellitus; HR, hazard ratio; CI, confidence interval; BMI, body mass index; d, day; wk, week.

Discussion

In this population-based retrospective cohort study, we found that statin therapy was associated with an increased risk of NOD compared to nonusers in patients with hypertension. These results are consistent with the findings of previously performed studies. For the Chinese population, 2 Taiwanese studies have shown that statin use was significantly associated with 15% and 27% increase in risk of incident diabetes, respectively.28,29 Our results are similar to the one reported in a British study (HR: 1.57, 95% CI: 1.54–1.59), which was based on the Clinical Practice Research Datalink.30 However, the risk seems even higher in two Korean studies (HR: 1.84, 95% CI: 1.63–2.09; HR: 1.99, 95% CI: 1.36–2.92).31,32 Compared to the results of these observational studies, we observed a moderately increased risk of NOD in statin users. On the other hand, the incidence of NOD was higher in our study than that in previous studies. One potential reason for this higher incidence could be that our study subjects were patients with hypertension who were at higher risk for incident diabetes than the general population.28,30

Since this was an observational study, the way we defined the cohort may induce bias. When researchers attempt to use a simple cohort approach to emulate the randomized controlled trial design, a form of bias called immortal time bias may be produced.33 Immortal time bias has been previously described, and several solutions have been proposed to prevent such bias including prescription time-distribution matching and time-dependent analysis.24

In our study, we used prescription time-distribution matching method, which means the overall distribution of index date of the nonusers is matched to the users’ index date (date of first statin prescription).33 The advantage of this design method is that it avoids the imbalance of the index date distribution between the 2 groups without apparent loss of study efficiency. With the balanced index date between 2 groups, we could better define the look-back period to identify the comorbidities and concomitant medicine intake.

We also conducted the subgroup analysis by sex, age, BMI, statin types, and statin potency. In multivariate Cox model, statin use was associated with increased risk of NOD in males (adjusted HR: 1.58, 95% CI: 1.42–1.77) and females (adjusted HR: 1.52, 95% CI: 1.38–1.67). A previous study suggested that incident NOD associated with statin use may be more common in women.10 We observed consistently increased risk of NOD by statin use across subgroups by age. This association was decreased but remained significant in older patients. As far as we know, there were 2 reports that demonstrated results similar to our finding.28,30 Macedo et al30 found that incidence rates of T2DM increased with age in nonusers, whereas in statin users incidence rates remained high. One of the potential mechanisms is that younger people at high risk of developing diabetes are accurately targeted for statin treatment, which would moderate the expected age gradient. Besides, our finding showed that patients with a BMI lower than 24 were at greater risk (adjusted HR: 1.65, 95% CI: 1.50–1.83) for NOD than those with BMI of 28 or higher (adjusted HR: 1.44, 95% CI: 1.17–1.77), even though a significant association was observed in all strata. Similar results had been reported by a study done in postmenopausal women.10 These findings have not been well explained. A further study to verify the weight effect on statin-related NOD is warranted. However, these finding suggested that we should pay more attention to the prescription of statins to younger patients with normal BMI.

Many studies have investigated the impact of specific statins on the incidence of NOD and have suggested that atorvastatin and simvastatin are more prone to be associated with increased risk of NOD.10,31,32 Our findings seemed to be consistent with these findings. The results of the multivariate Cox models demonstrated that atorvastatin users and simvastatin users had the highest risk of NOD. Due to the fact that only a few patients were treated with fluvastatin and rosuvastatin, the 95% CI of HR for these two subgroups were wider. The multivariate adjusted HRs were 1.44 (95% CI: 0.83–2.48) and 1.47 (95% CI: 1.19–1.83), respectively. Some studies have reported that statin-related NOD risk is likely potency dependent.34,35 We found that risk for NOD increased with increasing intensity of statins in our study. Compared to nonusers, the adjusted HRs were 1.44 (95% CI: 1.29–1.61), 1.65 (95% CI: 1.53–1.79), and 1.80 (95% CI: 0.86–2.82) for low-, moderate-, and high-potency statin users, respectively.

Limitations and strengths

Our results should be interpreted cautiously because of some limitations. First, we identified the patients with hypertension, diabetes mellitus, or other comorbidities by using the inpatient and outpatient diagnosis on ICD-10 codes. Our T2DM outcome definition, consisting of a diagnosis of T2DM or a prescription for a diabetes medication in outpatient and inpatient records, was probably unable to capture the actual diabetic patients. However, most of the observational studies based on the existing database suffered the same problems.34,36 Second, there were some confounding factors that we could not control. The information about the levels of dyslipidemia and hypertension was lacking, which were potentially important risk factors for the initiation of T2DM.17,18,37 Besides, dietary patterns are crucial determinants of T2DM,38 and data regarding this were not collected in the Yinzhou database. This might be the reason that we found relatively high HRs in our study. Third, we could only identify if the doctors prescribed drugs for the patients from the database; we have no idea whether the patients actually took the medications.

Our study also had some strengths. First, our study was the first to investigate the risk of statins in hypertensive population of Chinese mainland by using a well-developed health care database which reflects the physician behavior in a real-life setting. Second, we conducted a new users and new-onset T2DM design by excluding the patients who had took statins or had already been being diagnosed as having T2DM. Third, we use prescription time-distribution matching method to avoid immortal time bias and better identify the baseline characteristics. Forth, compared to other Asian studies,28,29,31,32 some crucial confounders (smoking status, alcohol use, and exercise frequency) were included in our analysis.

Conclusion

Our findings indicated that statin therapy was associated with an increased risk for incident NOD compared to nonusers of statins among hypertensive patients. Regardless of patients’ characteristics, statin types, or dosage, the association remained. The risk seems to be smaller among old and overweight patients who are at higher risk for CVD. These data provide better evidence about the association between statins and NOD in a real-world context in mainland People’s Republic of China. As the benefits of statins overwhelm the risk, especially for patients who are in a higher risk for CVD or with existing CVD, clinicians should followed the guideline to consider the statin therapy. As for those younger than 50 years old, the tradeoff between statin’s benefits and adverse effects should be well discussed. These findings should be interpreted cautiously as observational studies are subject to bias and residual confounding that cannot be ruled out.