Association of Visit-to-Visit Variability of Blood Pressure with Cardiovascular Disease among Type 2 Diabetes Mellitus Patients: A Cohort Study.

BACKGROUND: Increasing evidence has shown that visit-to-visit variability (VVV) of blood pressure (BP) is associated with an increased risk of cardiovascular disease (CVD). The objective of this study was to evaluate the impact of VVV of systolic blood pressure (SBP) and diastolic blood pressure (DBP) on the risk of CVD among patients with type 2 diabetes mellitus (T2DM) in China.
METHODS: We conducted a retrospective cohort study of 10,163 T2DM patients who were not previously diagnosed with CVD from January 2008 to December 2012 in Ningbo, China. The VVV of BP was calculated using five metrics, including standard deviation (SD), coefficient of variation (CV), variation independent of mean, average real variability, and successive variability (SV) of measurements, obtained over a 24-month measurement period. Hazard ratios and 95% confidence intervals (CIs) were estimated by Cox proportional hazards regression models for the associations of variability in BP with risk of CVD.
RESULTS: A total of 894 CVD events were observed during a median follow-up of 49.5 months. The hazard ratio in the highest quintile of SD of SBP was 1.24 (95% CI, 1.01 to 1.52) compared with patients in the lowest quintile. The association between higher VVV of DBP and risk of CVD was not consistent across different metrics and sensitivity analyses.
CONCLUSION: Higher VVV of SBP was associated with an increased risk of CVD, irrespective of the mean SBP level. Future studies are needed to confirm these findings.

INTRODUCTION

Cardiovascular disease (CVD) is a major cause of death and disability in patients with diabetes and consumes a significant portion of diabetes resources at a national level [1]. The approach recommended by the diabetic management guidelines is to maintain optimal blood pressure (BP) [23]. Recently, increasing attention has been paid to the prognostic value of visit-to-visit variability (VVV) for CVD events [45].

VVV of BP was first investigated by Grove et al. [6]. Rothwell et al. [7] found that VVV in systolic blood pressure (SBP) and maximum SBP were strong predictors for stroke, independent of mean SBP. Several prospective as well as retrospective studies of VVV of SBP were conducted in various patient populations [7891011]. A recent meta-analysis that included 23 studies concluded that VVV of SBP was associated with a risk of CVD incidence [12]. A limited number of studies have reported the association between VVV of SBP and CVD among patients with diabetes [1213141516], but only one study reported that VVV of diastolic blood pressure (DBP) was a significant predictor of CVD incidence among patients with type 2 diabetes mellitus (T2DM) [17]. Moreover, the effect of VVV of SBP on cardiovascular events may demonstrate heterogeneity between Chinese and non-Chinese populations. A study conducted in the UK revealed a lower risk of CVD among diabetic patients of Chinese ethnicity compared to the risk of patients of other ethnicities, such as Caucasian and Indian [18]. Evidence on the association between VVV of SBP and CVD among Chinese diabetic patients remains limited.

Thus, population-based studies among patients with diabetes are still needed to increase the body of evidence on the association of VVV of BP with risk of CVD events. We conducted a large retrospective cohort study aimed to evaluate the association between VVV of BP and CVD risk among T2DM patients in China.

METHODS

Study design

This population-based retrospective cohort study was conducted between January 1, 2008 and December 31, 2012 in Yinzhou District, Ningbo, Zhejiang province, China. Yinzhou is the largest district in Ningbo and has an adult population of approximately 961,000 people. Yinzhou has established a Health Information System (Yinzhou HIS) covering all health services, including public hospitals and primary-care clinics [1920]. In Yinzhou, CVDs and diabetes were required to be reported for disease surveillance and management by local general practitioners (GPs). GPs were also required to follow up with these patients at an interval of approximately 3 months, and outcome records were uploaded into the Chronic Disease Registry System. For CVD cases treated outside the Yinzhou region, the events were recorded and updated during the follow-up with the GP.

Data on T2DM patients were retrieved from the Chronic Disease Registry System in the Yinzhou HIS, which were coded as E11.X according to the 10th version of the International Classification of Diseases (ICD-10). A unique and encoded identifier was created for each patient in the Yinzhou HIS. We used this identifier to link all related databases to obtain the relevant covariates. The date of the first BP recording was defined as the baseline. The BP measurement period was defined as 24 months after baseline; thus, a total of nine BP measurements were used to calculate the VVV of BP. In the current study, T2DM patients with the first BP measurement record between January 1, 2008 and December 31, 2012 were included. The exclusion criteria were as follows: (1) Patients who had CVD events before baseline or during the BP measurement period and (2) Patients who were not followed regularly at an interval of approximately 3 months or had more than one missing BP value during the measurement period.

CVD events consisted of myocardial infarction and cerebrovascular disease, which were coded as I20 to I25 and I60 to I69 according to the ICD-10 in the Yinzhou HIS. Each patient was followed until the date of a cardiovascular event, death, the last known follow-up date or December 31, 2016, whichever occurred first.

BP variability measurements and visit-to-visit variability

The procedures for measuring BP during the BP measurement period followed standardized guidelines that included quality control measures [21]. The BPs of all the subjects were measured using a mercury sphygmomanometer while the patients were in a seated position after resting for 5 minutes. The first and fifth Korotkoff sounds were used to measure SBP and DBP, respectively. The quality control measures implemented for BP measurement were as follows. All hospitals and clinics in the region followed the same standardized protocol for BP measurement. At least two BP measurements were taken at each visit with a break between the measurements of at least 1 minute, and the average value was recorded by community doctors. The mercury sphygmomanometer was calibrated every six months. Talking and a full bladder were avoided during the measurement of BP. The average of the BP measurements was recorded by the doctors and uploaded into the Yinzhou HIS.

The standard deviation (SD) was used as the main metric of visit-to-visit blood pressure variability (BPV), and the other four metrics included the coefficient of variation (CV), successive variability (SV), average real variability (ARV), and variance independent of the mean (VIM). SD and CV account for the total variability around the mean BP across visits, while ARV and SV measure the differences in BP from one visit to the next. VIM was first used by Rothwell mainly to reduce the correlation between the measurement of BPV (take SD as an example) and the mean BP value. VIM was defined as SD/(meanλ), while λ was calculated through curve fitting. The mean BP value for each patient, averaged over all BP values recorded during the BP measurement period, was considered as a covariate. Detailed information for all metrics of VVV BP is listed in Supplementary Table 1.

Baseline covariates included the demographic characteristics, clinical parameters, and treatment modalities of the patients. Demographic characteristics included age, sex, body mass index (BMI), and smoking and alcohol drinking status. Current smoking was defined as having smoked at least one cigarette per day for 1 year or more. Current alcohol drinkers were defined as patients who consumed at least 100 g of any type of alcohol per week for the last year. Clinical parameters consisted of duration of diabetes and lipid profiles (including fasting plasma glucose, total cholesterol, triglycerides, high density lipoprotein cholesterol, and low density lipoprotein cholesterol). Treatment modalities included the use of antihypertensive drugs, oral antidiabetic drugs or insulin, and the use of lipid lowering treatments.

Statistical analysis

Descriptive analyses were used to analyze the baseline characteristics. Univariate analyses were used to test the differences among the quintiles of the SDs of SBP of patients and between patients included and excluded from the present study. Pearson's correlation coefficients were estimated between measures of VVV of BP and mean BP during the measurement period. Multivariate Cox proportional hazards regression models were performed to evaluate the effect of BP variability on cardiovascular outcomes after adjusting for age (continuous variable), sex (man or woman), duration of diabetes mellitus (continuous variable), smoking status (current or not), alcohol drinking status (current or not), BMI (continuous variable), use of glucose treatment (yes or no), use of BP lowering treatment (yes or no), and use of lipid lowering treatment (yes or no) (model I), or variables in model I plus mean BP values (continuous variable) (model II). Each BP variability measure was categorized by quintiles with a reference level of the lowest quintile. Linear trends across the quintiles were calculated by including the quintile as a continuous variable (ranging from 1 to 5) in the models.

Subgroup analyses were conducted through stratification by age (>60 and ≤60 years), sex (male or female), use of anti-hypertensive drugs (yes or no), and use of lipid lowering treatments to minimize potential bias. We also performed sensitivity analyses to examine the stability of our results. Patients with eight (n=5,642) and nine (n=4,521) BP readings were evaluated separately. Sensitivity analyses were also performed by restricting the analysis to patients with BP ≤140/90 and ≤130/80 mm Hg during the measurement period. We also excluded patients who had no information regarding their lipid profile at baseline, who were diagnosed with diabetes mellitus less than 1 year prior or who were followed for less than 1 year after the measurement period. In addition, we classified the SDs of SBP into five groups (≤5, 5.1 to 7.5, 7.6 to 10, 10.1 to 15, and >15 mm Hg) to determine whether maintaining an SBP within 10 mm Hg was an optimal target for BP variability control [13]. To explore the nonlinear relationships between CVD and the SDs of SBP as a continuous variable, we used restricted cubic splines with 4 knots placed at 5, 7.5, 10, and 15 mm Hg. All analyses were conducted in R software 3.3.2 (https://www.r-project.org/). A two-sided P<0.05 was considered statistically significant.

Ethics approval and consent to participate

This study conformed to the guidelines of the Declaration of Helsinki, and the study procedures were reviewed and approved by the ethics committee of the School of Medicine, Zhejiang University (No. 2018-028). Each patient agreed to participate and signed the informed consent form.

RESULTS

A total of 10,163 T2DM patients were included in the final analysis. Fig. 1 shows a flow chart of this study. We excluded 450 patients who had a history of CVD at baseline, 387 patients who developed CVD during the BP measurement period and 3,167 patients who were not regularly followed at an interval of 3 months or who had more than one missing BP value. The mean age at baseline was 61.3 years, 60.9% were female, and the median duration of DM was 40 months. The mean BP was 130.1/80.2 mm Hg at baseline. Table 1 summarizes the baseline characteristics of all patients included. The characteristics of the patients excluded from the current study are listed in Supplementary Table 2. During a median of 49.5 months of follow-up, a total of 894 cardiovascular events were recorded. Pearson's correlation coefficients between mean BP level and VVV BP metrics are shown in Supplementary Table 3.

Fig. 1

Study design for the investigation of the association of visit-to-visit variability (VVV) in blood pressure (BP) and cardiovascular disease (CVD). The BP measurements during the measurement period (0, 3, 6, 9, 12, 15, 18, 21, 24 months after baseline) were used to calculate the mean and VVV. The median follow-up period was 49.5 months.

Table 1

Baseline characteristics of included patients in the present study

Characteristic	Total	SD of SBP Quintile 1 (<4.51 mm Hg)	SD of SBP Quintile 2 (≥4.51 and <5.98 mm Hg)	SD of SBP Quintile 3 (≥5.98 and <7.62 mm Hg)	SD of SBP Quintile 4 (≥7.62 and <10.28 mm Hg)	SD of SBP Quintile 5 (≥10.28 mm Hg)	P value for difference
No. of patients	10,163	2,035	2,023	2,040	2,032	2,033
No. of cases	894	162	127	159	168	278
Cumulative incidence rate, <0.001	8.8	8.0	6.3	7.8	8.3	13.7	<0.001
Median follow-up time, mo, median (IQR)	49.5 (13.1)	49.2 (13.6)	49.3 (11.1)	49.4 (13.0)	49.4 (12.1)	50.4 (14.0)	<0.001
Person-years	60,420	12,067	12,040	12,145	12,036	12,133
Incidence density, case/1,000 person-yr	14.8	13.4	10.5	13.1	14.0	22.9	<0.001
Patients characteristics
Age	61.3±10.5	60.7±10.4	60.5±10.9	61.4±10.2	62.1±10.6a	64.1±10.0a	<0.001
Female sex	6,185 (60.9)	1,206 (59.3)	1,179 (58.3)	1,219 (59.8)	1,262 (62.1)a	1,319 (64.9)a	<0.001
Body mass index, kg/m2	24.1±2.6	24.2±2.4	24.1±2.6	24.0±2.7	24.1±2.6	24.2±2.6a	0.967
Duration of diabetes mel- litus, mo, median (IQR)	40.1 (9.3)	36.3 (7.9)	38.7 (8.2)	39.4 (9.1)	40.6 (10.8)a	45.4 (11.7)a	<0.001
Smoking	2,373 (23.3)	510 (25.1)	491 (24.3)	483 (23.7)	467 (23.0)	422 (20.8)a	0.016
Alcohol drinking	1,317 (13.0)	267 (13.1)	251 (12.4)	265 (13.0)	267 (13.1)	267 (13.1)	0.950
Education							<0.001
Illiterate	485 (4.8)	111 (5.5)	106 (5.2)	103 (5.0)	100 (4.9)	65 (3.2)
Literate	2,146 (21.1)	462 (22.7)	443 (21.9)	468 (22.9)	424 (20.9)	349 (17.2)
Elementary school	5,417 (53.3)	1,131 (55.6)	1,088 (53.8)	1,064 (52.2)	1,043 (51.3)	1,091 (53.7)
High school	1,945 (19.1)	302 (14.8)	346 (17.1)	379 (18.6)	424 (20.9)	494 (24.3)
College or above	169 (1.7)	29 (1.4)	40 (2.0)	26 (1.3)	40 (2.0)	34 (1.7)
Fasting plasma glucose, mmol/L	8.1±4.0	8.0±3.8	7.7±3.7a	8.0±4.4	8.4±4.7	8.6±3.3a	<0.001
Lipid profile	6.7±1.4	7.6±1.7	6.9±1.5	6.3±1.4a	6.0±1.3a	6.6±1.4a	<0.001
Total cholesterol, mmol/L Triglycerides, mmol/L	2.6±1.3	2.5±1.2	2.6±1.4	2.4±1.2	2.2±1.1a	2.2±1.3a	<0.001
High-density lipoprotein, mmol/L	2.1±0.6	2.3±0.7	2.2±0.6	2.1±0.6a	1.9±0.5a	2.2±0.6a	<0.001
Low-density lipoprotein, mmol/L	4.6±0.8	4.8±1.6	4.6±1.4	4.5±1.4a	4.2±1.3a	4.7±1.5a	<0.001
Any use of blood glucose treatment	9,625 (94.7)	1,929 (94.8)	1,926 (95.2)	1,929 (94.6)	1,915 (94.2)	1,926 (94.7)	0.736
Any use of BP lowering treatment	4,981 (49.0)	876 (43.0)	892 (44.1)	963 (47.2)a	1,048 (51.6)a	1,202 (59.1)a	<0.001
Any use of lipid lowering treatment	925 (9.1)	160 (7.9)	148 (7.3)	183 (9.0)a	206 (10.1)a	228 (11.1)a	<0.001
Blood pressure at baseline
SBP, mm Hg	130.1±11.7	129.8±6.4	129.6±7.4	129.9±9.2	131.4±11.4a	134.7±18.8a	<0.001
DBP, mm Hg	80.2±13.0	80.3±6.0	80.4±15.0	80.5±12.3	80.9±7.1	82.2±19.6a	<0.001
Blood pressure during measurement period
SBP, mm Hg	131.1±7.1	129.4±5.0	128.9±5.0	129.1±5.3	130.2±6.6a	132.9±10.9a	<0.001
DBP, mm Hg	80.9±6.5	80.2±6.7	79.8±6.0	79.9±5.7	80.2±6.5	80.9±7.4a	<0.001

Values are presented as mean±standard deviation or number (%) unless otherwise indicated. Continuous variables were tested by ANOVA or Kruskal-Wallis test, and categorical variables were tested by chi-square test.

SD, standard deviation; SBP, systolic blood pressure; IQR, interquartile range; DBP, diastolic blood pressure.

aP<0.05 for post hoc analysis with SD of SBP quintile 1 as a reference.

Table 2 presents the association of SBP VVV with the risk of cardiovascular events after adjusting for age, sex, duration of diabetes mellitus, BMI, smoking status, alcohol drinking status, fasting plasma glucose, use of any BP lowering treatments, use of any blood glucose lowering treatments and use of any lipid lowering treatments (model I). Patients in the highest quintile of SBP SD had a 29% higher risk of CVD compared with the lowest quintile (hazard ratio [HR], 1.29; 95% confidence interval [CI], 1.06 to 1.57). After additional adjustment for mean SBP during the measurement period (model 2), the corresponding HR for the highest quintile versus lowest was 1.24 (95% CI, 1.01 to 1.52), and the trend remained significant (P for trend <0.001). The results did not significantly change when we modeled the other four VVV SBP parameters. A similar but statistically nonsignificant association was observed between SD DBP and the risk of CVD (Q5 vs. Q1: HR, 1.22; 95% CI, 0.99 to 1.50) (Table 3). The results of BPV modeled as a continuous variable in different ranges of SBP and DBP are listed in Supplementary Table 4.

Table 2

HRs and 95% CIs for cardiovascular disease, associated with visit-to-visit SBP variability

Variable	Model Ia		Model IIb
	Adjusted HR (95% CI)	P for trend	Adjusted HR (95% CI)	P for trend
SD SBP, mm Hg		<0.001		<0.001
Q1 (<4.51)	1.00 (reference)		1.00 (reference)
Q2 (≥4.51 and <5.98)	0.76 (0.60–1.02)		0.77 (0.61–1.03)
Q3 (≥5.98 and <7.62)	0.92 (0.74–1.15)		0.95 (0.76–1.18)
Q4 (≥7.62 and <10.28)	0.90 (0.72–1.12)		0.90 (0.72–1.12)
Q5 (≥10.28)	1.29 (1.06–1.57)c		1.24 (1.01–1.52)c
CV SBP, mm Hg		<0.001		0.026
Q1 (<3.49)	1.00 (reference)		1.00 (reference)
Q2 (≥3.49 and <4.62)	0.77 (0.61–1.02)		0.79 (0.62–1.01)
Q3 (≥4.62 and <5.89)	0.96 (0.77–1.19)		0.98 (0.79–1.22)
Q4 (≥5.89 and <7.84)	0.89 (0.71–1.10)		0.90 (0.72–1.11)
Q5 (≥7.84)	1.24 (1.02–1.51)c		1.23 (1.01–1.50)c
VIM SBP, mm Hg		0.007		0.027
Q1 (<1.72)	1.00 (reference)		1.00 (reference)
Q2 (≥1.72 and <2.27)	0.74 (0.58–1.05)		0.79 (0.60–1.04)
Q3 (≥2.27 and <2.90)	0.93 (0.75–1.16)		0.95 (0.77–1.18)
Q4 (≥2.90 and <3.85)	0.89 (0.72–1.11)		0.91 (0.73–1.12)
Q5 (≥3.85)	1.20 (0.99–1.47)		1.20 (1.04–1.46)c
ARV SBP, mm Hg		0.031		0.025
Q1 (<4.00)	1.00 (reference)		1.00 (reference)
Q2 (≥4.00 and <5.71)	1.05 (0.84–1.30)		1.06 (0.85–1.32)
Q3 (≥5.71 and <7.63)	0.95 (0.76–1.19)		0.97 (0.77–1.21)
Q4 (≥7.63 and <10.57)	0.95 (0.77–1.19)		0.96 (0.77–1.20)
Q5 (≥10.57)	1.29 (1.06–1.58)c		1.24 (1.01–1.52)c
SV SBP, mm Hg		0.046		0.037
Q1 (<5.45)	1.00 (reference)		1.00 (reference)
Q2 (≥5.45 and <7.46)	0.94 (0.75–1.17)		0.96 (0.77–1.20)
Q3 (≥7.46 and <9.74)	0.86 (0.69–1.07)		0.89 (0.71–1.11)
Q4 (≥9.74 and <13.44)	0.88 (0.70–1.09)		0.89 (0.72–1.11)
Q5 (≥13.44)	1.26 (1.03–1.54)c		1.24 (1.01–1.52)c

HR, hazard ratio; CI, confidence interval; SBP, systolic blood pressure; SD, standard deviation; CV, coefficient of variation; VIM, variation independent of mean; ARV, absolute real variability; SV, successive variation.

aModel I included adjustment for age, sex, duration of diabetes mellitus, smoking status, alcohol drinking status, body mass index, baseline fasting plasma glucose, use of glucose treatment, use of antihypertensive treatment and use of lipid lowering treatment, bModel II included the same adjustment as Model I plus mean SBP during measurement period, cP<0.05 with Quintile 1 as reference.

Table 3

HRs and 95% CIs for cardiovascular disease, associated with visit-to-visit DBP variability

Variable	Model Ia		Model IIb
	Adjusted HR (95% CI)	P for trend	Adjusted HR (95% CI)	P for trend
SD DBP, mm Hg		0.006		0.007
Q1 (<3.50)	1.00 (reference)		1.00 (reference)
Q2 (≥3.50 and <4.46)	1.01 (0.80–1.26)		1.01 (0.81–1.26)
Q3 (≥4.46 and <5.39)	0.92 (0.73–1.15)		0.92 (0.74–1.16)
Q4 (≥5.39 and <6.78)	1.08 (0.87–1.35)		1.09 (0.88–1.35)
Q5 (≥6.78)	1.23 (1.00–1.51)c		1.22 (0.99–1.50)
CV DBP, mm Hg		0.020		0.018
Q1 (<4.35)	1.00 (reference)		1.00 (reference)
Q2 (≥4.35 and <5.59)	1.11 (0.89–1.38)		1.11 (0.89–1.38)
Q3 (≥5.59 and <6.80)	1.01 (0.81–1.27)		1.02 (0.81–1.28)
Q4 (≥6.80 and <8.47)	1.03 (0.83–1.29)		1.04 (0.84–1.30)
Q5 (≥8.47)	1.24 (1.00–1.52)c		1.24 (1.01–1.53)c
VIM DBP, mm Hg		0.099		0.070
Q1 (<0.74)	1.00 (reference)		1.00 (reference)
Q2 (≥0.74 and <0.97)	1.01 (0.81–1.25)		1.01 (0.81–1.26)
Q3 (≥0.97 and <1.20)	1.02 (0.82–1.27)		1.03 (0.83–1.28)
Q4 (≥1.20 and <1.49)	0.99 (0.80–1.23)		1.01 (0.81–1.25)
Q5 (≥1.49)	1.11 (0.90–1.37)		1.12 (0.91–1.38)
ARV DBP		0.684		0.719
Q1 (<3.13)	1.00 (reference)		1.00 (reference)
Q2 (≥3.13 and <4.29)	1.11 (0.90–1.36)		1.11 (0.90–1.36)
Q3 (≥4.29 and <5.57)	0.92 (0.74–1.15)		0.92 (0.74–1.15)
Q4 (≥5.57 and <7.14)	0.92 (0.74–1.14)		0.92 (0.74–1.14)
Q5 (≥7.14)	1.12 (0.92–1.38)		1.12 (0.91–1.37)
SV DBP		0.183		0.194
Q1 (<4.25)	1.00 (reference)		1.00 (reference)
Q2 (≥4.25 and <5.67)	1.24 (0.99–1.54)		1.25 (0.99–1.54)
Q3 (≥5.67 and <7.07)	0.99 (0.79–1.24)		1.00 (0.80–1.25)
Q4 (≥7.07 and <9.06)	0.97 (0.78–1.21)		0.98 (0.78–1.22)
Q5 (≥9.06)	1.29 (1.05–1.58)c		1.28 (1.04–1.58)c

HR, hazard ratio; CI, confidence interval; DBP, diastolic blood pressure; SD, standard deviation; CV, coefficient of variation; VIM, variation independent of mean; ARV, absolute real variability; SV, successive variation.

aModel I included adjustment for age, sex, duration of diabetes mellitus, smoking status, alcohol drinking status, body mass index, baseline fasting plasma glucose, use of glucose treatment, use of antihypertensive treatment and use of lipid lowering treatment, bModel II included the same adjustment as Model I plus mean DBP during measurement period, cP<0.05 with Quintile 1 as reference.

Fig. 2 shows that the VVV of SBP is an independent predictor of CVD across several sensitivity analyses, except when analyses were restricted to patients aged ≤60 years or patients with a maximum BP ≤130/80 mm Hg during the measurement period. When patients were classified into five groups (≤5 mm Hg, n=2,769; 5.1 to 7.5 mm Hg, n=3,196; 7.6 to 10 mm Hg, n= 2,018; 10.1 to 15 mm Hg, n=1,668; and >15 mm Hg, n=512 for SD of SBP), the results showed that HRs for the 10.1 to 15 and >15 mm Hg groups were statistically significant (1.29 [95% CI, 1.03 to 1.61] and 1.56 [95% CI, 1.16 to 2.09], respectively), compared with the lowest group (≤5 mm Hg). The results of the restricted cubic splines in the Cox regression models showed that the adverse effect reached significance when the SD of SBP was higher than approximately 9 mm Hg (Supplementary Fig. 1).

Fig. 2

Hazard ratio (95% confidence interval [CI]) for cardiovascular disease of Q5 compared to Q1 of standard deviation of systolic blood pressure (SBP) in different subgroups. Data were adjusted for sex, age, smoking status, alcohol drinking status, duration of diabetes mellitus, body mass index, baseline fasting plasma glucose, use of blood glucose treatment, use of antihypertensive drugs, use of lipid lowering drugs and mean SBP values. Baseline low density lipoprotein cholesterol, high density lipoprotein cholesterol, total cholesterol, and triglycerides were further adjusted when patients without baseline lipid profiles were excluded. BP, blood pressure; DM, diabetes mellitus.

DISCUSSION

In the current study, we found positive associations between VVV of SBP and risk of CVD among T2DM patients. The associations were consistent after multivariable adjustment and across several sensitivity analyses, demonstrating that VVV of SBP was an independent predictor of CVD incidence irrespective of mean SBP values. The results of VVV of DBP were not consistent across sensitivity analyses. Our findings could provide additional evidence for the prognostic value of VVV of BP for CVD events in diabetic patients.

A number of studies have evaluated the associations between VVV of BP and risk of CVD events. The majority of these studies were conducted in the general population or patients with hypertension [222324252627], and several studies also demonstrated the adverse effect of high VVV in diabetes [28293031]. Ohkuma et al. [14] used data from the Action in Diabetes and Vascular Diseases: Preterax and Diamicron Modified Release Controlled Evaluation-Observational study (ADVANCE-ON) and found that VVV in SBP improved risk prediction beyond predictions using traditional risk factors, including mean SBP. A territory-wide retrospective cohort study among Chinese primary care T2DM patients has identified that greater SBP variability may increase the risk of all-cause CVD mortality [13], which is consistent with our results. Moreover, to our knowledge, this is the first study that presents the association between VVV of DBP and CVD events among Chinese T2DM patients. Significant associations were observed in the highest quintile of DBP compared with the lowest quintile when CV and SV were used. Additional studies are needed to confirm whether the VVV of DBP is a prognostic predictor of CVD and whether it can provide additional information beyond that provided by the VVV of SBP.

Our study found that CVD risk increased only in the highest quintile of VVV of SBP, and CVD risk was insignificantly lower in the second and third quintiles of VVV of SBP than in the first quintile (The survival curve stratified by quintiles of SD SBP was shown in Supplementary Fig. 2). Although the nonlinearity test for variability metrics of SBP did not reach statistical significance, this finding suggests that the relationship between BP variability and CVD may not be perfectly linear, which has been reported in several studies [2232]. There may be a threshold effect for the impact of BPV, and further studies are needed to address this issue. When we classified patients into five groups by arbitrary SBP values, group 4 (10.1 to 15 mm Hg) and group 5 (>15 mm Hg) both had a higher risk of CVD compared with that of the lowest group (SD SBP <5 mm Hg), and group 5 had the highest risk. When we modeled SD SBP as a continuous variable in Supplementary Table 4, the results showed that every 1 SD of SBP increase (6.5 mm Hg) was associated with a 4% increase in the risk of CVD. The other four metrics (CV, VIM, ARV, SV) were also statistically significant when modeled as continuous variables. We also found that SBP variability was more closely associated with cerebrovascular disease after stratification by disease subtypes (data not shown), but these findings merit further investigation because the number of ischemic heart disease cases was relatively small in this study.

Evidence for the association between diastolic BPV and CVD was limited and varied compared to the association of CVD and systolic BPV [3334]. However, it has been reported that SBP is a more important prognostic factor than DBP in the middle-aged population [3536]. In this study, we found that mean SBP during the measurement period had prognostic value, whereas mean DBP did not (data not shown). On the other hand, SBP was easier to accurately measure than DBP by using a mercury manometer [37]. Thus, we believe that VVV of SBP is a more robust prognostic indicator than VVV of DBP.

Although evidence for the detrimental effect of long-term SBP variability has been consistent, it has still not been applied in clinical practice [38]. Possible reasons may include the following: the variability needs months or years to be assessed, and the metrics of VVV BP cannot be calculated easily. Our study showed that the adverse effect of VVV of SBP reached significance when the SD of SBP was used, which was consistent with the previous finding that the optimal target for SBP SD should be <10 mm Hg [13]. This finding might be used in clinical settings for BP control among T2DM patients, as patients with diabetes usually have BP measurements taken in the hospital.

One possible mechanism of the association between SBP variability and CVD was reported as suboptimal BP control [3940]. However, our results remained consistent after we restricted our analyses to patients with BP ≤140/90 mm Hg (n=6,555). This finding indicated that T2DM patients with a higher VVV of BP still had a higher risk of CVD even when BP values were under control. When the included patients were further restricted to those with SBP ≤130/80 mm Hg during the measurement period (n=1,480), no associations were observed between SD of SBP and CVD. Future studies with a large sample size are needed to validate this finding.

Arterial stiffness, endothelial dysfunction and coronary artery calcification have also been reported as potential explanations for the adverse effects of high VVV [4142434445]. Faramawi et al. [43] reported that increased BP variability was associated with endothelial dysfunction, which accelerates atherosclerosis, thus increasing the risk of CVD. Shimbo et al. [42] found that SDs of SBP were inversely associated with aortic distensibility. Unfortunately, information on vascular wall status was not available in the current study. Further studies are required to better underlie the mechanistic relationships.

Our study has several strengths, including the large sample size, BP variability measurements and the relatively long-term follow-up time. We used a rigorous inclusion criterion of at least eight visits to calculate BP variability because a recent meta-analysis reported that the number of visits used to quantify VVV of BP affects the magnitude of the VVV itself [46]. We also used five metrics to measure VVV of BP, which added to the robustness of our work. However, limitations of our study should be noted. First, causality cannot be concluded due to the natural limitation of the retrospective observational cohort design. Patient characteristics and clinical information collected in this study were uploaded into the HIS in the past when patients were undergoing physical examinations or other health tests, which reduced the occurrence of recall bias. Further prospective studies are required to confirm our results. Second, selection bias cannot be ruled out because we excluded 3,167 patients due to irregular follow-up or missing BP values. However, the baseline characteristics of patients who were excluded and of those who were included were similar; thus, the risk of selection bias of the current study is acceptable. Third, there may be ascertainment bias of CVD because the accuracy of the Yinzhou HIS was not validated. This ascertainment bias could influence CVD events at baseline and during the follow-up period, especially for the patients who left Yinzhou district and did not return again, who could only be followed via telephone by local GPs. These patients were easily lost to follow-up. Fourth, information on patients' adherence to treatment was not available in this study. However, because the frequency of visits to clinics somehow reflects the adherence to treatment, our inclusion of only patients with 8 or 9 BP readings reduced the heterogeneity of adherence to treatment to some extent. Fifth, data on potential cofounders, such as certain kind of anti-diabetic drugs, anti-hypertensive drugs and hemoglobin, were not available in the current study, which may limit the generalizability of our results. Finally, although the standard guideline for the BP measurement process was consolidated between all clinics and hospitals as mentioned above, manual measurement using a mercury sphygmomanometer may cause more human error than measurements conducted by an automated device. Thus, the method of BP measurement may be a confounding factor of BPV and merits further investigation.

In summary, our study provides evidence to support the association between VVV of SBP and risk of CVD among T2DM patients, regardless of mean BP levels. Moreover, our findings confirmed that an optimal target for a VVV SD of SBP control might be <10 mm Hg. Future studies are required to explore the potential mechanism underlying the increased risk of CVD associated with increased VVV of SBP.