Association Between Visit-to-Visit Mean Arterial Pressure Variability and the Risk of Ischemic Heart Disease and Ischemic Stroke Among Patients With Hypertension in Thailand

Abstract

Background:

Blood pressure (BP) variability has been recognized as a significant risk factor for cardiovascular diseases (CVD). We aim to evaluate the association between mean arterial pressure (MAP) variability and the increased risk of ischemic heart disease (IHD) and ischemic stroke (IS) among hypertensive patients in Thailand.

Methods:

We analyzed data from the Thailand DM/HT study, which included hypertensive patients nationwide in 2014 to 2015 and 2018. MAP variability was computed based on the MAP values across 3 visits within 1 year and expressed as standard deviation (SD). We used multivariable log-binomial regression models to evaluate the associations between MAP variability and the risk of IHD and IS.

Results:

Among 92 854 individuals, 594 new-onset IHD events (0.64%) and 187 IS incidents among 95 486 individuals (0.20%). Compared to the lowest quartile (Q1), higher quartiles of SD were associated with increased risk of IHD, with adjusted risk ratios (aRRs) of 1.06 (95% confidence interval [CI]: 0.82-1.38) for Q2, 1.35 (95% CI: 1.06-1.72) for Q3, and 1.50 (95% CI: 1.18-1.90) for Q4. Similarly, higher SD quartiles raised the risk of IS, with aRRs of 1.35 (95% CI: 0.83-2.20) for Q2, 1.56 (95% CI: 0.98-2.48) for Q3, and 1.97 (95% CI: 1.26-3.07) for Q4, when compared to Q1.

Conclusion:

Our study demonstrated that higher visit-to-visit MAP variability in hypertensive patients was strongly associated with an increased risk of CVD. We emphasize the importance of incorporating BP variability into management strategies to help reduce the risk of CVD in these patients.

Keywords

blood pressure variability mean arterial pressure ischemic heart disease ischemic stroke hypertension

Background

Prior studies consistently demonstrated the strong association between high blood pressure (BP) and the risk of cardiovascular disease (CVD), including ischemic heart disease (IHD) and ischemic stroke (IS).^1
-6 Additionally, high BP variability was shown to significantly raise the risk of CVD and mortality, independent of mean BP levels.^7
-11 Consequently, BP variability has emerged as a crucial risk factor for CVD and can provide more accurate predictions for clinical outcomes in adults.¹²

The cyclic BP curve is determined by 2 distinct phenotypes: pulse pressure, which corresponds to large artery stiffness and wave reflections, and mean arterial pressure (MAP), which reflects the resistance of the microvascular network.¹³ MAP provides a more comprehensive view of the average pressure in the arteries throughout the cardiac cycle. In contrast, systolic blood pressure (SBP) and diastolic blood pressure (DBP) only capture the peak pressures during systole and diastole, respectively.¹⁴ Previous studies revealed that MAP serves as a more effective predictor of CVD risk than SBP and DBP.¹⁵ Furthermore, prior research reported that MAP and overall CVD events were positively associated, including IHD, IS, and CVD deaths.^15,16

Few studies assessed the potential impact of MAP variability on CVD outcomes. For example, recent evidence from the Atherosclerosis Risk in Communities (ARIC) study indicated that increased visit-to-visit MAP variability is a significant risk factor for heart failure (HF) and all-cause mortality.¹⁷ Furthermore, although machine learning has been implemented using various variables in the existing dataset to predict CVD, such as coronary heart disease, MAP variability has not yet been incorporated into these models.¹⁸ This omission reflects insufficient evidence regarding the relationship between MAP variability and CVD.

Hypertension (HTN) is a common chronic medical condition affecting approximately one-fourth of Thai adults.¹⁹ Among Thai patients with HTN, IHD and IS are the most prevalent CVD complications, and they remain the leading causes of mortality in the Thai population.^5,6,20
-22 Despite the high prevalence of HTN, limited information is available on the BP variability, particularly in MAP, and its association with CVD outcomes in patients with HTN. Therefore, our study seeks to investigate the association between visit-to-visit MAP variability and the risk of IHD and IS among patients with HTN who receive continuous care in Thailand.

Methods

Study Design and Subjects

We employed data from the 2014, 2015, and 2018 Thailand DM/HT studies. The study was supported by the National Health Security Office (NHSO) to conduct a series of annual surveys designed to evaluate the clinical outcomes of individuals with HTN and/or type 2 diabetes (T2D) aged 20 years or older who were receiving care at clinics in hospital nationwide in Thailand.²³ The study sample comprised 101 661 individuals with HTN who were enrolled in 2014, 2015, and 2018.

Figure 1 presents the eligible participants and the total number of participants included in the analysis. Regarding the availability of data based on Thailand DM/HT study data, BP information for each patient’s last 3 visits within 1 year was extracted from the medical records. To examine the association between visit-to-visit MAP variability and IHD incidence, we excluded 4011 individuals with a history of IHD at baseline (before 12 months from the date of data collection) and those without complete 3 BP measures within 1 year. This resulted in a final analysis sample of 92 854 participants without IHD at baseline and having 3 BP measurement records. Moreover, we assessed the association between visit-to-visit MAP variability and IS incidence by excluding 1255 individuals with a prior history of IS and those lacking complete 3 BP measurements within 1 year. This led to a final analysis sample of 95 486 participants without IS at baseline and having 3 BP measurement records.

Figure 1.

Study flowchart with detailed study exclusion information.

Data Collection

Sakboonyarat et al published the details of data collection for the Thailand DM/HT study.^23,24 The study involved a proficient registered nurse meticulously reviewing and extracting the data using a case report form (CRF) based on a standardized protocol. Then, the data were sent to the data management unit. The study participants’ characteristics encompassed several factors, such as sex, age, health insurance scheme, geographical region, hospital levels where patients receive HTN care, comorbidities (such as T2D, dyslipidemia, renal insufficiency, and atrial fibrillation), smoking status, weight, height, and antihypertensive medications, which were duly documented for analysis.

Typically, the attending physician documents information regarding patients with HTN diagnosed with CVD complications using the International Classification of Diseases, Tenth Revision (ICD-10) codes. Additionally, well-trained hospital staff conduct audits of the medical record reviews to ensure that the ICD-10 codes and diagnoses are used accurately. In the Thailand DM/HT study, the skilled registered nurse abstracted information on CVD complications, including IHD and IS, based on the ICD-10 codes and completed the CRF.

Outcomes

According to data from the Thailand DM/HT study, information on IHD and IS was extracted from medical records using these ICD-10 codes. The incidence of IHD was defined as the occurrence of IHD as a new onset in patients with hypertension within a span of 12 months. The definition was based on the ICD-10 I20-I25 or a recorded history of coronary revascularization.²⁵ Furthermore, the incidence of IS was defined as a new case of IS among hypertensive individuals within a period of 12 months, determined in accordance with the ICD-10: I63.²⁵

Exposures

The study participants’ BP information, which included SBP and DBP, was collected in the last 3 visits within 1 year and utilized for analysis purposes. The MAP was calculated through the following formula: MAP = DBP + [1/3(SBP – DBP)].²⁶ Then, the average MAP levels were computed across 3 visits for each participant. The visit-to-visit MAP variability measurements were calculated as follows: the SD_MAP = $\sqrt[2]{\frac{\sum_{i = 1}^{3} {(M A P_{i} - m e a n M A P)}^{2}}{2}}$ . Furthermore, we considered the visit-to-visit MAP var-iability, calculated as the coefficient of variation (CV) = $\frac{S D_{M A P}}{M e a n_{M A P}}$ in the sensitivity analysis.

Covariates

In this study, the covariates considered comprised demographic variables and health-related information. Age was categorized into 2 categories: <65 years and ⩾65 years. Regions included North, Central, Northeast, and South. The health insurance scheme encompassed universal health coverage (UHC), civil servant medical benefits (CSMB), social security (SS), and others. The clinic’s location was classified as a regional hospital, provincial hospital, and community hospital. T2D, dyslipidemia (DLP), renal insufficiency, and a history of atrial fibrillation were defined based on medical records. Moreover, tobacco use was determined based on the smoking status recorded in the medical records as never smoked and ever smoked. Furthermore, height and weight data were collected and utilized to calculate the body mass index (BMI) and categorized as <25 kg/m² and ⩾25 kg/m². Duration of HTN was designated as <10 years, 10–19 years, and 20 years and over. Antihypertensive medication included angiotensin-converting enzyme inhibitors/angiotensin receptor blockers (ACEIs/ARBs), beta-blockers use, calcium channel blockers (CCBs) use, and diuretics use.

Statistical Analysis

Descriptive statistics were employed to determine the characteristics of the study participants. Categorical variables were presented as percentages, whereas continuous variables were reported as mean, standard deviation (SD), median, and interquartile range (Q1-Q3).

According to the Thailand DM/HT study data, which included a short follow-up of a closed cohort (1-year follow-up) with a binary outcome, we used log-binomial regression to estimate adjusted risk ratios (aRRs) and 95% confidence intervals (CIs) for the incidence of IHD and IS associated with higher MAP variability. Adjusted risk difference (RD) and 95% CI were estimated using logistic regression analysis with margins command (adjrr). The multivariable analysis models were adjusted for potential confounders as follows: model 1: age, sex; model 2: variables in model 1 plus regions, hospital level, health insurance scheme, T2D, DLP, renal insufficiency, history of atrial fibrillation, HTN duration, ACEIs/ARBs use, beta-blockers use, CCBs use, diuretics use, smoking status, BMI; and model 3: variables in model 2 plus the mean of MAP visits 1 to 3. Furthermore, a trend analysis in each model was conducted to evaluate the dose-response relationship between MAP variability and the outcomes, presenting it as p for trend. SD_MAP spline curves for the average adjusted prediction of IHD and IS incidence were created using the margins command after the multivariable log-binomial regression model 3.

Subgroup analyses of key variables (sex, age, T2D, dyslipidemia, renal insufficiency, BMI, and smoking status) were performed to compare the incidence of IHD and IS between the highest (SD Q4) and the other 3 quartiles (SD Q1 + Q2 + Q3) of MAP variability. An interaction term between key variable and MAP variability was individually added to the adjusted log-binomial regression model 3, and the aRRs and CIs for these associations were estimated.

Sensitivity analyses were performed using the CV as an alternative method to evaluate MAP variability and its association with the risk of IHD and IS. In addition to dividing SD_MAP and CV_MAP into quartiles, we categorized the exposures into tertiles to evaluate the association between MAP variability and the risk of IHD and IS. To address the issue of missing BP records, we employed marginal structural models (MSMs) to analyze the association between MAP variability and the risk of IHD and IS. The methods for MSMs are detailed in the Supplementary Data . All statistical analyses were conducted using StataCorp. 2021. Stata Statistical Software: Release 17, College Station, TX: StataCorp LLC. All tests were 2-sided, and a P-value < .05 was considered statistically significant.

Ethics Consideration

The study was reviewed and approved by the Institutional Review Board, the Royal Thai Army Medical Department, in compliance with international guidelines such as the Declaration of Helsinki, the Belmont Report, CIOMS Guidelines, and ICH-GCP (approval number S055h/65_Exp). A waiver of documentation of informed consent was obtained for the use of secondary data and was granted by the Institutional Review Board, the Royal Thai Army Medical Department.

Results

Characteristics of Study Participants

Table 1 lists the baseline characteristics of the study participants. Around two-thirds of the participants were women. The mean age of the participants was 63.9 ± 11.7 years and 64.1 ± 11.8 years in IHD and stroke outcome cohorts, respectively. The mean values of the average SBP, DBP, and MAP across 3 visits were 133.8 ± 11.5, 76.4 ± 8.4, and 95.5 ± 8.0 mmHg, respectively.

Table 1.

Baseline characteristics of study participants.

Variables	Baseline for ischemic heart disease	Baseline for ischemic stroke
Variables	n (%)	n (%)
Total	92 854	95 486
Sex
Men	35 019 (37.7)	36 197 (37.9)
Women	57 835 (62.3)	59 289 (62.1)
Age, years
Mean (SD)	63.9 (11.7)	64.1 (11.8)
Median (Q1-Q3)	64.0 (56.0-72.0)	64.0 (56.0-73.0)
<65	48 337 (52.1)	49 049 (51.4)
⩾65	44 462 (47.9)	46 380 (48.6)
Geographical regions
North	23 713 (25.5)	24 409 (25.6)
Central	31 762 (34.2)	32 807 (34.4)
Northeast	23 052 (24.8)	23 519 (24.6)
South	14 327 (15.4)	14 751 (15.5)
Health scheme
Universal health coverage	68 538 (73.8)	70 570 (73.9)
Civil servant medical benefits	19 311 (20.8)	19 843 (20.8)
Social security	4250 (4.6)	4300 (4.5)
Others	755 (0.8)	773 (0.8)
Hospital level
Regional hospital	6784 (7.3)	6959 (7.3)
Provincial hospital	18 658 (20.1)	19 235 (20.1)
Community hospital	67 412 (72.6)	69 292 (72.6)
Diabetes
No	79 816 (86.0)	82 036 (85.9)
Yes	13 038 (14.0)	13 450 (14.1)
Dyslipidemia
No	10 322 (11.1)	10 572 (11.1)
Yes	82 532 (88.9)	84 914 (88.9)
Renal insufficiency
No	81 474 (87.7)	83 545 (87.5)
Yes	11 380 (12.3)	11 941 (12.5)
History of atrial fibrillation
No	91 927 (99.0)	94 449 (98.9)
Yes	927 (1.0)	1037 (1.1)
Hypertension duration, years
<10	75 590 (81.4)	77 496 (81.2)
10-19	16 448 (17.7)	17 134 (17.9)
20 and over	816 (0.9)	856 (0.9)
Body mass index, kg/m ²
Mean (SD)	25.0 (4.7)	25.0 (4.7)
Smoking status
Never	77 887 (83.9)	79 952 (83.7)
Ever smoke	14 967 (16.1)	15 534 (16.3)
ACEIs/ARBs use
No	38 959 (42.0)	39 964 (41.9)
Yes	53 895 (58.0)	55 522 (58.2)
Beta-blockers use
No	69 152 (74.5)	70 326 (73.7)
Yes	23 702 (25.5)	25 160 (26.4)
Calcium channel blockers use
No	38 708 (41.7)	40 360 (42.3)
Yes	54 146 (58.3)	55 126 (57.7)
Diuretics use
No	69 167 (74.5)	71 075 (74.4)
Yes	23 687 (25.5)	24 411 (25.6)
Average systolic BP, mmHg
Mean (SD)	133.8 (11.5)	133.8 (11.5)
Median (Q1-Q3)	133.3 (126.7-140.0)	133.3 (126.7-140.0)
Average diastolic BP, mmHg
Mean (SD)	76.4 (8.4)	76.4 (8.4)
Median (Q1-Q3)	76.7 (70.7-82.0)	76.7 (70.7-82.0)
Average MAP, mmHg
Mean (SD)	95.5 (8.0)	95.5 (8.0)
Median (Q1-Q3)	95.4 (90.2-100.4)	95.3 (90.1-100.3)
SD of MAP
Mean (SD)	5.7 (3.6)	5.7 (3.6)
Median (Q1-Q3)	5.1 (3.1-7.5)	5.1 (3.1-7.6)
CV of MAP
Mean (SD)	0.06 (0.04)	0.06 (0.04)
Median (Q1-Q3)	0.05 (0.03-0.08)	0.05 (0.03-0.08)

Abbreviations: ACEIs, angiotensin-converting enzyme inhibitors; ARBs, angiotensin receptor blockers; BP, blood pressure; MAP, mean arterial pressure; SD, standard deviation; CV, coefficient of variation

Association Between MAP Variability and the Risk of IHD Among Patients With HTN

Table 2 presents the association between MAP variability and the risk of IHD. Among the 92 854 patients with HTN without a history of IHD, 594 new-onset IHD events were observed (0.64%). After adjusting for potential confounders, the higher quartiles of SD Q3 and Q4 were associated with a 35% and 50% increased risk of incident IHD, respectively, compared to the lowest quartile of SD (Q1; aRR: 1.35 and 95% CI: 1.06-1.72; aRR: 1.50 and 95% CI: 1.18-1.90) Also, a dose-response relationship was noted between MAP variability measured by SD and incident IHD, with P for trend <.001. Figure 2A reveals that the average adjusted prediction of IHD increased with MAP variability measured by SD when it reached a higher level. Figure 3A and Supplemental Table 1 demonstrate the association between MAP variability and the risk of IHD on an additive scale. Subgroup analyses of key variables showed that the risk of IHD in the highest quartile of SD (Q4) was greater than that in the lower quartiles of SD (Q1 + Q2 + Q3; Table 3).

Table 2.

Association of visit-to-visit MAP variability measured by SD with the risk of ischemic heart disease and ischemic stroke.

MAP Variability	Total	Incidence	Univariable analysis		Model 1*		Model 2*		Model 3*
MAP Variability	N	N (%)	crude RR (95% CI)	P-value	adjusted RR (95% CI)	P-value	adjusted RR (95% CI)	P-value	adjusted RR (95% CI)	P-value
Ischemic heart disease
Total	92 854	594 (0.64)
SD Q1	24 146	117 (0.48)	Ref.		Ref.		Ref.		Ref.
SD Q2	22 540	118 (0.52)	1.08 (0.84-1.39)	.552	1.07 (0.83-1.39)	.580	1.06 (0.82-1.38)	.655	1.06 (0.82-1.38)	.656
SD Q3	23 164	162 (0.70)	1.44 (1.14-1.83)	.002	1.41 (1.11-1.79)	.004	1.35 (1.06-1.72)	.017	1.35 (1.06-1.72)	.016
SD Q4	23 004	197 (0.86)	1.77 (1.41-2.22)	<.001	1.66 (1.32-2.09)	<.001	1.48 (1.18-1.88)	.001	1.50 (1.18-1.90)	.001
P for trend		<0.001	<0.001		<0.001		<0.001		<0.001
Ischemic stroke
Total	95 486	187 (0.20)
SD Q1	24 696	30 (0.12)	Ref.		Ref.		Ref.		Ref.
SD Q2	23 127	38 (0.16)	1.35 (0.84-2.18)	.216	1.35 (0.84-2.17)	.222	1.35 (0.83-2.19)	.229	1.35 (0.83-2.20)	.224
SD Q3	23 874	50 (0.21)	1.72 (1.10-2.71)	.018	1.69 (1.08-2.66)	.023	1.58 (0.99-2.51)	.055	1.56 (0.98-2.48)	.061
SD Q4	23 789	69 (0.29)	2.39 (1.56-3.66)	<.001	2.28 (1.48-3.50)	<.001	2.10 (1.35-3.28)	.001	1.97 (1.26-3.07)	.003
P for trend		<0.001	<0.001		<0.001		0.001		0.002

Multivariable analysis; Model 1: adjusted for age and sex; Model 2: adjusted for Model 1 + regions, hospital level, health scheme, diabetes, dyslipidemia, renal insufficiency, history of atrial fibrillation, hypertension duration, ACEIs/ARBs use, beta-blockers use, calcium channel blockers use, diuretics use, smoking status, body mass index; Model 3: adjusted for Model 2+ mean of MAP visit 1 to 3.

Abbreviations: SD, standard deviation; MAP, mean arterial pressure; RR, risk ratio; CI, confidence interval

Figure 2.

Spline curve of SD_MAP for the average adjusted prediction of (A) ischemic heart disease incidence and (B) ischemic stroke incidence and 95% CI: (A) Average adjusted prediction of ischemic heart disease incidence among hypertensive patients, adjusting for age and sex, regions, hospital level, health scheme, diabetes, dyslipidemia, renal insufficiency, history of atrial fibrillation, hypertension duration, ACEIs/ARBs use, beta-blockers use, calcium channel blockers use, diuretics use, smoking status, body mass index, and mean of MAP visit 1 to 3, (B) Average adjusted prediction of ischemic stroke incidence among hypertensive patients, adjusting for age and sex, regions, hospital level, health scheme, diabetes, dyslipidemia, renal insufficiency, history of atrial fibrillation, hypertension duration, ACEIs/ARBs use, beta-blockers use, calcium channel blockers use, diuretics use, smoking status, body mass index, and mean of MAP visit 1 to 3. SD: standard deviation, MAP: mean arterial pressure, CI: confidence interval, aRR: Adjusted risk ratio for 1 unit of SD increase.

Figure 3.

Association between MAP variability measured by standard deviation and (A) ischemic heart disease incidence and (B) ischemic stroke incidence on an additive scale.

Table 3.

Association of the highest quartile of MAP variability group (SD Q4) compared with the lower MAP variability group (SD Q1 + Q2 + Q3) for ischemic heart disease in subgroups.

Variables	adjusted RR (95% CI)*	P-value	P for interaction
Overall	1.31 (1.10-1.57)	.003
Sex			.300
Men	1.21 (0.92-1.58)	.172
Women	1.42 (1.12-1.80)	.004
Age, years			.173
<65	1.09 (0.80-1.47)	.586
⩾65	1.50 (1.20-1.87)	<.001
Type 2 diabetes			.294
No	1.38 (1.14-1.68)	.001
Yes	0.98 (0.60-1.58)	.922
Dyslipidemia			.033
No	2.07 (1.28-3.34)	.003
Yes	1.22 (1.01-1.48)	.042
Renal insufficiency			.213
No	1.22 (1.00-1.50)	.056
Yes	1.68 (1.15-2.45)	.007
Body mass index, kg/m ²			.083
<25	1.53 (1.22-1.91)	<.001
⩾25	1.01 (0.74-137)	.950
Smoking status			.151
Never	1.39 (1.14-1.70)	.001
Ever smoke	1.05 (0.70-1.57)	.824

Adjusted for age and sex, regions, hospital level, health scheme, diabetes, dyslipidemia, renal insufficiency, history of atrial fibrillation, hypertension duration, ACEIs/ARBs use, beta-blockers use, calcium channel blockers use, diuretics use, smoking status, body mass index, and mean of MAP visit 1 to 3.

Abbreviations: SD, standard deviation; MAP, mean arterial pressure; RR, risk ratio; CI, confidence interval

Association Between MAP Variability and the Risk of IS Among Patients With HTN

Out of the 95 486 patients with HTN without a history of IS, 187 incidents of new-onset IS were observed (0.20%). After controlling for potential confounding factors, the results demonstrated that patients in the higher quartile of MAP variability (SD Q4) had a considerably higher risk of developing incident IS than those in the lowest quartile of SD (Q1), with aRR of 1.97 (95% CI: 1.26-3.07; Table 2). A dose-response relationship between MAP variability measured by SD and incident IS was observed (P for trend = .002). Figure 2B illustrates that the average adjusted prediction of IS increased with MAP variability measured by SD when it reached a higher level. Figure 3B and Supplemental Table 1 reveal the association between MAP variability and the risk of IS on an additive scale. Subgroup analyses of key variables presented differential risks of IS in subgroups of sex, age, T2D, DLP, renal insufficiency, BMI, and smoking status. However, the interaction testing revealed no heterogeneity (Table 4).

Table 4.

Association of the highest quartile of MAP variability group (SD Q4) compared with the lower MAP variability group (SD Q1 + Q2 + Q3) for ischemic stroke in subgroups.

Variables	adjusted RR (95% CI)*	P-value	P for interaction
Overall	1.51 (1.10-2.06)	.010
Sex			.547
Men	1.71 (1.08-2.71)	.023
Women	1.35 (0.88-2.07)	.167
Age, years			.165
<65	1.11 (0.66-1.87)	.693
⩾65	1.86 (1.25-2.77)	.002
Type 2 diabetes			.978
No	1.50 (1.06-2.11)	.022
Yes	1.55 (0.73-3.29)	.252
Dyslipidemia			.103
No	0.72 (0.24-2.19)	.562
Yes	1.63 (1.17-2.26)	.004
Renal insufficiency			.097
No	1.34 (0.95-1.90)	.100
Yes	2.52 (1.19-5.34)	.016
Body mass index, kg/m ²			.723
<25	1.51 (1.02-2.22)	.038
⩾25	1.57 (0.92-2.67)	.097
Smoking status			.275
Never	1.39 (0.97-1.98)	.072
Ever smoke	2.00 (1.03-3.90)	.041

Abbreviations: SD, standard deviation; MAP, mean arterial pressure; RR, risk ratio; CI, confidence interval

Sensitivity Analyses

The sensitivity analyses assessed the association between MAP variability, measured by the CV, and the risk of IHD and IS. As detailed in Supplemental Table 2, these findings consistently supported the primary analysis results, revealing a substantial association between the highest quartile of MAP variability (CV Q4) and an increased risk of IHD and IS. In an alternative approach that categorized MAP variability into tertiles, we observed associations between MAP variability—both measured by SD and CV—and the risk of IHD and IS, as presented in Supplemental Tables 3 and 4. Furthermore, MSMs were employed to account for patients with HTN who did not have 3 BP readings documented in their medical records. The results were consistent with the primary analysis, demonstrating that higher MAP variability is associated with an increased risk of IHD and IS, as shown in Supplemental Table 5.

Discussion

In the present study, a positive association was observed between visit-to-visit MAP variability and the risk of both IHD and IS among patients with HTN receiving continuous care in Thailand. This relationship remained significant even after adjusting for potential confounding factors, including mean MAP, and was robust across a range of sensitivity analyses. Our findings suggest that visit-to-visit MAP variability may represent a potentially important risk factor for IHD and IS in patients with HTN. Furthermore, they may offer a more accurate estimate of risk than the mean MAP alone.

To the best of our knowledge, this study represents the first comprehensive investigation to date into the relationship between visit-to-visit MAP variability and IHD and IS. Our findings reveal that individuals in the highest quartile of SD_MAP have a 50% increased risk of developing IHD—equating to an additional 3 cases per 1000 individuals—and a 97% increased risk of experiencing IS, which corresponds to an additional case per 1000 individuals, compared to those in the lowest quartile of SD_MAP.

MAP is a fundamental driving force for the perfusion of vital organs and is also associated with end-organ damage, such as ischemia and infarction, which can occur if MAP drops for an extended period.²⁶ Prior studies reported that MAP and overall CVD events are positively associated, including IHD, stroke, and CVD mortality.^15,16 However, few studies evaluated the potential impact of MAP variability on CVD outcomes. Recent findings from the ARIC study stated that long-term MAP variability and HF and all-cause mortality are related.¹⁷ Our study demonstrates the association of visit-to-visit MAP variability with IHD and IS in a large population with HTN in Thailand. This finding aligns with the evidence that MAP variability may be a potential risk factor for IHD and IS. Therefore, minimizing MAP variability to achieve optimal BP management may be beneficial.

Previous studies demonstrated that SBP trajectories over time are superior predictors of stroke, coronary heart disease, and cardiovascular mortality compared to a single SBP assessment.^10,11,27 For instance, a large-scale study conducted on 3 285 684 US veterans revealed that elevated SBP variability in individuals with and without HTN was linked to higher risks of IHD, stroke, and all-cause mortality.¹¹ The current study builds upon these earlier findings. It demonstrates that 1-year visit-to-visit MAP variability, similar to SBP variability, is associated with the risks of IHD and IS. The findings support and emphasize the importance of BP variability for health-related outcomes, particularly concerning MAP.

Regarding the subgroup analysis of key variables, our findings demonstrated that the highest quartile of MAP variability was associated with a higher risk for IHD and IS. Moreover, the effect of measure modification of DLP on the association between MAP variability and IHD on a multiplicative scale was observed. Our analysis indicated a weaker association between MAP variability and IHD among those with DLP. This finding may be explained by DLP being a significant risk factor for IHD.²⁸ Individuals with DLP are more likely to have a higher risk for IHD than those without DLP; therefore, it dilutes the effect of higher MAP variability on the risk of IHD.

Consistent with prior research, elevated BP variability was linked to an enhanced likelihood of cardiovascular complications. This phenomenon may be elucidated by a number of proposed mechanisms, including alterations in the elastic characteristics of blood vessels and aortic distensibility,²⁹ as well as lifestyle factors such as noncompliance with antihypertensive medications³⁰ and exposure to emotional stimuli and psychological stressors.^31,32 Moreover, a prior investigation on patients with T2D revealed that BP variability exhibited a positive relationship with pulse wave velocity and a reverse association with ankle-brachial index, thereby revealing a potential association between BP variability and vascular dysfunction.^33,34 Notably, increased BP variability exerts stress on blood vessels, leading to endothelial dysfunction, a possible risk factor for the atherosclerotic process and a plausible contributor to developing IHD and IS.^35,36 Furthermore, the MAP variability can result in a direct decrease in myocardial perfusion and impact blood perfusion in cerebral tissues.²⁶ As MAP acts as a primary driving force for vital organ perfusion, repeated transient reductions in MAP over time could potentially raise the risk of relative hypoperfusion in both cardiac and cerebral tissues.^26,37

Another possible mechanism that could explain our findings is the role of autonomic dysfunction related to CVD. Autonomic dysfunction pertains to issues within the autonomic nervous system, including both the sympathetic and parasympathetic branches. This dysfunction may manifest as increased BP variability.³⁸ Consequently, BP variability could serve as a proxy indicator, allowing for noninvasive measurements to monitor changes in autonomic function.³⁹

Our findings display that MAP variability, including SD and CV, can be applied in clinical practice to reduce CVD risk among patients with HTN. Healthcare providers may incorporate MAP variability into HTN treatment by adjusting antihypertensive medications based on MAP variability as a target. This approach focuses on reducing variability and lowering mean BP to achieve the optimal target BP.⁴⁰ Additionally, BP variability may help identify individuals at high risk for CVD among patients with HTN. In our subgroup analysis, we found a stronger relationship between MAP variability and IHD and IS in individuals aged 65 and older (aRR:1.50 for IHD; aRR: 1.86 for IS) compared to those under 65 (aRR of 1.09 for IHD and 1.11 for IS). Therefore, older patients with HTN with higher BP variability should be closely monitored for CVD complications, including IHD and IS. Therefore, we suggest conducting further detailed studies that clarify the role of MAP variability and the underlying mechanisms involved in developing IHD and IS to better our understanding of optimal BP management and provide better guidance for clinical practice.

Several strengths characterized our study. First, we were able to analyze a large sample size of almost 100 thousand individuals, providing a robust dataset for our analyses. Additionally, the study population was highly representative of Thai patients with HTN receiving continuous care nationwide, thereby enhancing the generalizability of our findings. Furthermore, this study, while focused on Thai individuals, may yield findings relevant to other populations regarding HTN management, mainly through the lens of BP variability rather than BP levels.

However, our study also had several limitations. One major limitation was the relatively short duration of visit-to-visit MAP measurements, which were taken over just 1 year. Therefore, the findings indicated a cross-sectional association. Moreover, due to the potential confounding factors inherent in an observational study, causality between MAP variability, IHD, and IS cannot be established. Typically, BP measurements at outpatient clinics follow the Thai guidelines for treating HTN, which suggest the standard procedure for BP measurement. Nonetheless, the data used in this study was collected from the Thailand DM/HT study, which may have introduced measurement errors in BP readings taken by hospital staff nationwide. Therefore, misclassification of the exposure may occur and may bias the results. Despite our efforts to account for potential confounding factors using multivariable analyses, residual confounding remains, particularly concerning dietary behavior, physical activity, and adherence to antihypertensive medication, which were not available in the Thailand DM/HT study data. This may introduce bias into our findings. This study included Thai patients with HTN receiving continuous care nationwide; therefore, the generalizability may be limited to this population. However, subgroup analyses of key variables and several sensitivity analyses yielded consistent results, thereby bolstering the robustness of our findings.

Conclusion

Our study involving Thai patients with HTN receiving continuous care revealed that higher visit-to-visit MAP variability was strongly associated with an increased risk of IHD and IS. This finding suggests that MAP variability could be a potential risk factor for both IHD and IS. In addition to monitoring BP levels, it is crucial to account for BP variability in managing HTN. This approach can guide medication adjustments to achieve optimal BP targets and help identify patients at high cardiovascular risk.

Supplemental Material

sj-docx-1-cic-10.1177_11795468251342338 – Supplemental material for Association Between Visit-to-Visit Mean Arterial Pressure Variability and the Risk of Ischemic Heart Disease and Ischemic Stroke Among Patients With Hypertension in Thailand

Supplemental material, sj-docx-1-cic-10.1177_11795468251342338 for Association Between Visit-to-Visit Mean Arterial Pressure Variability and the Risk of Ischemic Heart Disease and Ischemic Stroke Among Patients With Hypertension in Thailand by Boonsub Sakboonyarat, Jaturon Poovieng and Ram Rangsin in Clinical Medicine Insights: Cardiology

Footnotes

Acknowledgements

The authors thank the National Health Security Office in Thailand for supporting the Thailand DM/HT study dataset.

Abbreviations

HTN: Hypertension

CVD: Cardiovascular disease

IHD: Ischemic heart disease

IS: Ischemic stroke

BP: Blood pressure

SBP: Systolic blood pressure

DBP: Diastolic blood pressure

MAP: Mean arterial pressure

SD: Standard deviation

CV: Coefficient of variation

T2D: Type 2 diabetes

DLP: Dyslipidemia

ORCID iD

Boonsub Sakboonyarat

Ethical Considerations

This study was reviewed and approved by the Institutional Review Board, Royal Thai Army Medical Department, following international guidelines including the Declaration of Helsinki, the Belmont Report, CIOMS Guidelines, and the International Conference on Harmonization of Technical Requirements for Registration of Pharmaceuticals for Human Use–Good Clinical Practice (ICH–GCP) (approval number S055h/65_Exp).

Consent to Participate

Due to the use of secondary data, a waiver of documentation of informed consent was utilized. The Institutional Review Board, Royal Thai Army Medical Department, approved an informed consent waiver.

Consent for Publication

Consent for publication is not applicable.

Author Contributions

BS conducted the study concept, analyzed and interpreted the data, and was a major contributor to writing the manuscript. JP conducted the study concept and interpretation of the data. RR performed the study concept/design, data collection, and interpretation of the data. All authors read and approved the final manuscript.

Funding

The author(s) disclosed receipt of the following financial support for the research, authorship, and/or publication of this article: This study was funded by the Phramongkutklao College of Medicine, Bangkok, Thailand.

Declaration of Conflicting Interests

The author(s) declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.

Data Availability Statement

Data cannot be shared publicly because the data set contains identifying information; additionally, the data belong to the Thailand DM/HT study of the Medical Research Network of the Consortium of Thai Medical Schools (MedResNet). Thus, ethical restrictions exist on the data set. Data are available from the Thai National Health Security Office (NHSO), Bangkok, Thailand (contact via ) for researchers who meet the criteria for access to confidential data.

Supplemental Material

Supplemental material for this article is available online.

References

Razo

Welgan

Johnson

, et al. Effects of elevated systolic blood pressure on ischemic heart disease: a burden of proof study. Nat Med. 2022;28(10):2056-2065.

Malik

Georgakis

Vujkovic

, et al. Relationship between blood pressure and incident cardiovascular disease: linear and nonlinear mendelian randomization analyses. Hypertension. 2021;77(6):2004-2013.

Luo

Cheng

Zhang

, et al. Association between high blood pressure and long term cardiovascular events in young adults: systematic review and meta-analysis. BMJ. 2020;370:m3222.

Brunström

Carlberg

Association of blood pressure lowering with mortality and cardiovascular disease across blood pressure levels: a systematic review and meta-analysis. JAMA Intern Med. 2018;178(1):28-36.

Sakboonyarat

Rangsin

Prevalence and associated factors of ischemic heart disease (IHD) among patients with diabetes mellitus: a nation-wide, cross-sectional survey. BMC Cardiovasc Disord. 2018;18(1):151.

Chantkran

Chaisakul

Rangsin

Mungthin

Sakboonyarat

Prevalence of and factors associated with stroke in hypertensive patients in Thailand from 2014 to 2018: A nationwide cross-sectional study. Sci Rep. 2021;11:17614.

Hata

Arima

Rothwell

, et al. Effects of visit-to-visit variability in systolic blood pressure on macrovascular and microvascular complications in patients with type 2 diabetes mellitus: the ADVANCE trial. Circulation. 2013;128(12):1325-1334.

Muntner

Whittle

Lynch

, et al. Visit-to-visit variability of blood pressure and coronary heart disease, stroke, heart failure, and mortality: a cohort study. Ann Intern Med. 2015;163(5):329-338.

Hastie

Jeemon

Coleman

, et al. Long-term and ultra long–term blood pressure variability during follow-up and mortality in 14 522 patients with hypertension. Hypertension. 2013;62(4):698-705.

10.

Rothwell

Howard

Dolan

, et al. Prognostic significance of visit-to-visit variability, maximum systolic blood pressure, and episodic hypertension. Lancet. 2010;375(9718):895-905.

11.

Gosmanova

Mikkelsen

Molnar

, et al. Association of systolic blood pressure variability with mortality, coronary heart disease, stroke, and renal disease. J Am Coll Cardiol. 2016;68(13):1375-1386.

12.

Messerli

Hofstetter

Rimoldi

Rexhaj

Bangalore

Risk factor variability and cardiovascular outcome: JACC review topic of the week. J Am Coll Cardiol. 2019;73(20):2596-2603.

13.

Safar

Boudier

HS.

Vascular development, pulse pressure, and the mechanisms of hypertension. Hypertension. 2005;46(1):205-209.

14.

Kandil

Soliman

Alghamdi

Jennings

El-Baz

Using mean arterial pressure in hypertension diagnosis versus using either systolic or diastolic blood pressure measurements. Biomedicines. 2023;11(3):849.

15.

Sesso

Stampfer

Rosner

, et al. Systolic and diastolic blood pressure, pulse pressure, and mean arterial pressure as predictors of cardiovascular disease risk in men. Hypertension. 2000;36(5):801-807.

16.

Hadaegh

Shafiee

Hatami

Azizi

Systolic and diastolic blood pressure, mean arterial pressure and pulse pressure for prediction of cardiovascular events and mortality in a Middle Eastern population. Blood Press. 2012;21(1):12-18.

17.

Liu

Chen

Zhang

, et al. Long-term visit-to-visit mean arterial pressure variability and the risk of heart failure and all-cause mortality. Front Cardiovasc Med. 2021;8:665117.

18.

Sapra

Sandhu

Goyal

Intelligent method for detection of coronary artery disease with ensemble approach. In: Lecture Notes in Electrical Engineering, vol. 668, 2021.

19.

Aekplakorn

Thai National Health Examination Survey VI (2019-2020). 2019. Accessed October 3, 2023. https://online.fliphtml5.com/bcbgj/znee/#p=187

20.

Feigin

Stark

Johnson

CO.

Global, regional, and national burden of stroke and its risk factors, 1990–2019: a systematic analysis for the Global Burden of Disease Study 2019. Lancet Neurol. 2021;20(10):795-820.

21.

Naghavi

Abajobir

Abbafati

, et al. Global, regional, and national age-sex specific mortality for 264 causes of death, 1980–2016: a systematic analysis for the Global Burden of Disease Study 2016. Lancet. 2017;390(10100):1151-1210.

22.

Wang

Mubarik

, et al. Global burden of ischemic heart disease and attributable risk factors, 1990–2017: a secondary analysis based on the global burden of disease study 2017. Clin Epidemiol. 2021;13:859-870.

23.

The National Health Security Office T. DAMUS. Data archival for maximum utilization system. 2022. Accessed October 7, 2022. https://dmht.thaimedresnet.org/

24.

Sakboonyarat

Pima

Chokbumrungsuk

, et al. National trends in the prevalence of glycemic control among patients with type 2 diabetes receiving continuous care in Thailand from 2011 to 2018. Sci Rep. 2021;11(1):14260.

25.

World Health Organization. International Statistical Classification of Diseases and Related Health Problems, 10th Revision. 5th ed. 2016. Accessed November 18, 2022. https://apps.who.int/iris/handle/10665/246208

26.

DeMers

Wachs

Physiology, mean arterial pressure. StatPearls [Internet]. Published online 2020.

27.

Tielemans

Geleijnse

Menotti

, et al. Ten-year blood pressure trajectories, cardiovascular mortality, and life years lost in 2 extinction cohorts: the Minnesota business and professional men study and the Zutphen study. J Am Heart Assoc. 2015;4(3):e001378.

28.

Lorenzatti

Toth

PP.

New perspectives on atherogenic dyslipidaemia and cardiovascular disease. Eur Cardiol Rev. 2020;15:1-9.

29.

Shimbo

Shea

McClelland

, et al. Associations of aortic distensibility and arterial elasticity with long-term visit-to-visit blood pressure variability: the multi-ethnic study of atherosclerosis (MESA). Am J Hypertens. 2013;26(7):896-902.

30.

Muntner

Levitan

Joyce

, et al. Association between antihypertensive medication adherence and visit-to-visit variability of blood pressure. J Clin Hypertens. 2013;15(2):112-117.

31.

Tomitani

Kanegae

Kario

The effect of psychological stress and physical activity on ambulatory blood pressure variability detected by a multisensor ambulatory blood pressure monitoring device. Hypertens Res. 2023;46(4):916-921.

32.

Eguchi

Hoshide

Schwartz

Shimada

Kario

Visit-to-visit and ambulatory blood pressure variability as predictors of incident cardiovascular events in patients with hypertension. Am J Hypertens. 2012;25(9):962-968.

33.

Okada

Fukui

Tanaka

, et al. Visit-to-visit variability in systolic blood pressure is correlated with diabetic nephropathy and atherosclerosis in patients with type 2 diabetes. Atherosclerosis. 2012;220(1):155-159.

34.

McClary

Massey

Ankle brachial index. Published online 2019.

35.

Eto

Toba

Akishita

, et al. Reduced endothelial vasomotor function and enhanced neointimal formation after vascular injury in a rat model of blood pressure lability. Hypertens Res. 2003;26(12):991-998.

36.

Davignon

Ganz

Role of endothelial dysfunction in atherosclerosis. Circulation. 2004;109(23_suppl_1):III27-32.

37.

Vedel

Holmgaard

Rasmussen

, et al. Perfusion Pressure cerebral infarct (PPCI) trial-the importance of mean arterial pressure during cardiopulmonary bypass to prevent cerebral complications after cardiac surgery: study protocol for a randomised controlled trial. Trials. 2016;17(1):247-311.

38.

Hadaya

Ardell

JL.

Autonomic modulation for cardiovascular disease. Front Physiol. 2020;11:617459.

39.

Ooi

Lim

Seng

, et al. Non-invasive parameters of autonomic function using beat-to-beat cardiovascular variations and arterial stiffness in hypertensive individuals: a systematic review. Biomed Eng Online. 2024;23(1):23.

40.

Sheikh

Sobotka

Garg

, et al. Blood pressure variability in clinical practice: past, present and the future. J Am Heart Assoc. 2023;12(9):e029297.

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.04 MB