Household Air Pollution and Blood Pressure, Vascular Damage, and Subclinical Indicators of Cardiovascular Disease in Older Chinese Adults

Abstract Background Limited data suggest that household air pollution from cooking and heating with solid fuel (i.e., coal and biomass) stoves may contribute to the development of hypertension and vascular damage. Methods Using mixed-effects regression models, we investigated the associations of household air pollution with blood pressure (BP) and vascular function in 753 adults (ages 40–79 years) from 3 diverse provinces in China. We conducted repeated measures of participants’ household fuel use, personal exposure to fine particulate air pollution (PM2.5), BP, brachial–femoral pulse wave velocity (bfPWV), and augmentation index. Ultrasound images of the carotid arteries were obtained to assess intima–media thickness (CIMT) and plaques. Covariate information on sociodemographics, health behaviors, 24-h urinary sodium, and blood lipids was also obtained. Results Average estimated yearly personal exposure to PM2.5 was 97.5 µg/m3 (SD: 79.2; range: 3.5–1241), and 65% of participants cooked with solid fuel. In multivariable models, current solid fuel use was associated with higher systolic (2.4 mm Hg, 95% CI: −0.4, 4.9) and diastolic BP (1.4 mm Hg, 95% CI: −0.1, 3.0) and greater total area of plaques (1.7 mm2, 95% CI: −6.5, 9.8) compared with exclusive use of electricity or gas stoves. A 1 − ln(µg/m3) increase in PM2.5 exposure was associated with higher systolic (1.5 mm Hg, 95% CI: 0.2, 2.7) and diastolic BP (1.0 mm Hg, 95% CI: 0.4, 1.7) and with greater CIMT (0.02 mm, 95% CI: 0.00, 0.04) and total area of plaques (4.7 mm2, 95% CI: −2.0, 11.5). We did not find associations with arterial stiffness, except for a lower bfPWV (−1.5 m/s, 95% CI: −3.0, −0.0) among users of solid fuel heaters. Conclusions These findings add to limited evidence that household air pollution is associated with higher BP and with greater CIMT and total plaque area.


BACKGROUND
Limited data suggest that household air pollution from cooking and heating with solid fuel (i.e., coal and biomass) stoves may contribute to the development of hypertension and vascular damage.

METHODS
Using mixed-effects regression models, we investigated the associations of household air pollution with blood pressure (BP) and vascular function in 753 adults (ages 40-79 years) from 3 diverse provinces in China.We conducted repeated measures of participants' household fuel use, personal exposure to fine particulate air pollution (PM 2.5 ), BP, brachial-femoral pulse wave velocity (bfPWV), and augmentation index.Ultrasound images of the carotid arteries were obtained to assess intima-media thickness (CIMT) and plaques.Covariate information on sociodemographics, health behaviors, 24-h urinary sodium, and blood lipids was also obtained.

CONCLUSIONS
These findings add to limited evidence that household air pollution is associated with higher BP and with greater CIMT and total plaque area.
High blood pressure (BP) is the leading cause of death, responsible for 10.8 million premature deaths in 2019 worldwide. 1In China, 250 million people live with hypertension, which is responsible for 24% of China's deaths and 14% of disability-adjusted life years. 2Urban air pollution shows strong and consistent associations with high BP and the development of cardiovascular disease in studies mostly conducted in high-income cities. 3,4 Less understood is whether household air pollution from biomass and coal (i.e., solid fuel) stoves also contributes to hypertension and vascular damage, their underlying preclinical conditions, and the exposure-response relationships between them. 5ver 2.5 billion people globally and 600 million Chinese 3 use solid fuel cookstoves and heaters that emit air pollution into homes and communities.Several cohort studies indicate that the primary use of a solid fuel stove increases the risk of mortality and nonfatal cardio-respiratory events by 7%-19%. 5These studies are supported by analyses showing higher BP and greater carotid intima-media thickness (CIMT) in women using wood-burning cookstoves, [6][7][8] and that switching to less-polluting stoves lowers BP. 9 However, despite widespread exposure to household air pollution, there remain very few exposure-response studies 10 and even fewer that included men, 5 which form a significant knowledge gap for policy.
Leveraging repeated measurements from the INTERMAP China Prospective (ICP) study, 11 a geographically diverse cohort in rural and peri-urban villages, we investigated associations of estimated yearly exposure to fine particulate matter <2.5 microns (PM 2.5 ) with the severity of arteriosclerosis and vascular disease, as measured by BP, arterial stiffness, CIMT, and carotid plaque features.We also assessed whether adults using solid fuel stoves at present and over the previous 2 decades had greater vascular damage than adults using only gas or electric (i.e., clean fuel) stoves.

Study population
The ICP Study is an ongoing longitudinal study to investigate the nutritional and environmental risk factors for cardiovascular disease, including BP and vascular function. 11It is an extension of the INTERMAP study, a 1996 investigation of nutrient intake and BP in 839 adults aged 40-59 years that lived in 17 villages in northern (Beijing and Shanxi) and southern (Guangxi) China. 12In 2015, we recontacted 680 of the surviving INTERMAP study participants and enrolled 574 into the ICP Study (participation rate: 84%). 11An additional 208 adults ages 40-59 years (i.e., same ages as INTERMAP participants in 1996) were randomly selected from village lists and enrolled into the study (participation rate: 86%) to ensure sufficient sample size as some INTERMAP participants had died or moved away.Participants lived in rural and peri-urban villages that were within 20-30 km of urban areas, and provided written informed consent.Study protocols received ethics committee approvals from all investigator institutions.

Data collection
Participants attended village clinics on up to 7 occasions during a total of 6 data collection campaigns across our study sites that took place between August 2015 and November 2016.We conducted at least 1 winter and summer campaign in the northern sites to capture seasonal differences, which resulted 1 full campaign in Guangxi and 2 full campaigns in Beijing and Shanxi.A full campaign consisted of 2 days of data collection, including repeated measurements of BP.During the first clinic visit in each campaign, staff verbally administered individual and household questionnaires to participants.During the second and third clinic visits in each campaign, we repeatedly measured participants' anthropometrics, exposure to air pollution, BP, and arterial stiffness, and collected 24 hours urine for biochemistry.Ultrasound images of the carotid arteries and whole blood samples were obtained once in June (Shanxi), September (Beijing), and November (Guangxi) of 2016 (Supplementary Figure S1 online).The images were collected at the same time as other study measurements in Beijing and Guangxi, and in a separate third reduced data collection campaign in Shanxi where BP and surveys were also conducted.Details about measurements are provided below and published in a study protocol. 11

Current and historical household fuel use and intensity of use
We developed visual libraries of stoves and fuels for each site that were used to collect information on both current and historical fuel use during the first clinic visit.For each stove pictured, participants indicated whether they had used it during the past 20 years and, if so, answered questions about the years (in 5-year intervals) and frequency of use, stove location in the home, and fuel.Stove-fuel categories were combined with reported use data to construct fuel-based exposure metrics (Table 1).Cookstoves refer to devices used for household cooking or water boiling activities, whereas heating stoves are those primarily used for space heating during winter. 13

Personal exposure to air pollution
We measured personal exposure to PM 2.5 for up to 96 hours in 2 seasons to estimate seasonal and yearly exposures (Table 1).Participants wore waist packs with air samplers that collected PM 2.5 on PTFE filters.In a random 80% subsample, we monitored compliance in wearing the air samplers by placing pedometers (HJ-321 Tri-Axis, Omron, Japan) inside the waist packs.Filters were analyzed for mass at the Wisconsin State Laboratory of Hygiene.Details on PM 2.5 measurement are published elsewhere. 11,14

Cardiovascular outcomes
][17][18] Following at least 5-min of quiet rest, BP was measured on the participant's right arm at least 3 times using an automated oscillometric device (Omron HEM-907).Up to 5 measurements were taken if the difference between the last 2 was >5 mm Hg.BP was taken on 2 consecutive days during each measurement campaign, except for Shanxi's third campaign which had only 1 day of BP measurement due to study logistics.The average of the last 2 measurements from consecutive days of each campaign was used for statistical analysis, preserving season-specific BP and arterial stiffness variations.Hypertension was defined as use of antihypertension medication, or systolic (≥140 mm Hg) and/or diastolic (≥90 mm Hg) BP.Up to 6 measurements, 3 per campaign, of bfPWV and AIx (Vicorder, Smart Medical, UK) were taken in the supine position with cuffs around the participant's right upper arm and right upper thigh.
Measurements of CIMT and carotid artery plaques were conducted by a single technician in B-mode ultrasonography using a portable Mindray Z6 Ultrasound and 7-10 MHz linear array transducer.At least 3 images of the far wall of 3 10-mm segments of the bilateral carotid artery were obtained for both arteries, i.e., 18 images total. 18A clinician analyzed each image offline using an automated software (Carotid Analyzer Vascular Research Tools 6, Medical Imaging Applications, Coralville) to assess CIMT, total area of plaques, and GSM, and the average of all measurements was used for statistical analysis.A plaque was defined as a localized thickness of CIMT ≥1.5 mm or a focal raised lesion of ≥0.5 mm, and at least 3 images were obtained.Plaques were measured by manually tracing the perimeter and computing their cross-sectional area.Blood pressure, ultrasound measurements, and image-based analyses were conducted by staff who were blinded to the participants' exposure status.

Covariates
Staff administered individual and household questionnaires in each campaign to obtain information on household income and sociodemographics and assessed individual cardiovascular risk factors.Serum specimens were collected once and processed in a centralized testing facility for lipids, and 24-hour urine samples were collected twice in each campaign and analyzed for sodium.Anthropometrics were measured in each campaign using a height ruler attached to a mechanical scale and SECA measuring tapes.

Statistical analysis
Mixed-effects regression models with participant-and village-specific random intercepts were used to estimate associations between household air pollution, measured by fuel use and exposure to PM 2.5 , and outcomes with repeated measures (BP, bfPWV, AIx).These models account for the lack of independence between repeated measurements in the same participant and between participants in the same village.Crosssectional models with village-specific random intercepts were conducted for outcomes with single measurements (CIMT, total plaque area, GSM).Known and suspected confounders a Clean fuel stoves included those powered by gas or electricity and solid fuel stoves included those powered by coal, wood, crop residues, or other forms of biomass.
b This group was excluded from the analysis given the small sample size.
included in the final models were age, age, sex, household income, waist circumference, height (AIx, bfPWV), smoking status, secondhand smoke exposure, alcohol consumption, physical activity, 24-h urinary sodium (BP, bfPWV, AIx), time of BP measurement (PM 2.5 models), and total/HDL cholesterol ratio. 19Waist circumference and body mass index were highly correlated (r = 0.8); thus, we adjusted only for waist circumference which was a stronger confounder.We did not adjust for outdoor temperature because it was strongly correlated with PM 2.5 .We also did not adjust for previous diagnosis of diabetes or hypertension in the main analysis given that they could be confounders or along the causal pathway.Seasonal measurements of PM 2.5 were time-weighted (7 nonheating and 5 heating months) to estimate average yearly exposure for models predicting CIMT and plaque features.Homogeneity of variance and normality assumptions were evaluated by residual vs. fitted plots and QQ plots.Data were inspected to determine if participants with missing variables were systematically different from rest of the study population based on other information collected, and found that the missing variables appeared to be missing at random.We handled missing variables with multiple imputation for income (n = 87), lipids (n = 107), urinary sodium (n = 26), waist circumference (n = 16), and body mass index (n = 17) (Supplementary Table S1 online).Median village-and campaign-specific values were used for participants missing the time of BP measurement (n = 29).
We assessed effect modification for variables previously shown to modify air pollution-vascular associations including by sex, age, antihypertensive medication use, and geographic region for all outcomes and by season for BP and arterial stiffness.Sensitivity analyses included (i) limiting fuel use analyses to indoor stoves; (ii) excluding participants with self-reported diabetes, which could be a confounder or along the causal pathway; (iii) adjusting for heating fuel in models with cooking fuel as the exposure and vice versa; and, (iv) adjusting for hypertension status in models of arterial stiffness, CIMT, and plaque features since high BP could be a confounder or along the causal pathway.Statistical analyses were conducted in R (www.r-project.org,v. 1.2.1335).

Characteristics of the study population
At the time of the first survey, participants were 40-79 years old (median = 63 years), 55% female, and 23% were current tobacco smokers (Table 2).Exclusive users of clean fuel, on average, had higher education and income, were more likely to smoke, less likely to regularly engage in physical activity, and more likely to report taking antihypertensive medication and have diabetes compared with solid fuel users.Current exclusive users of clean fuel had lower personal exposures to PM 2.5 than current users of solid fuel cookstoves [mean (SD): 91.7 (63.0) vs. 100.6µg/m 3 (86.4)].
Over half (52%) of the participants had hypertension, 72% of whom reported taking antihypertensive medication.Mean BP and brachial-femoral PWV were higher in winter (Table 3).Blood pressure (P = <0.01)and AIx (P = <0.01),but not bfPWV (P = 0.16), differed across provinces.CIMT was greater in northern participants compared with those in Guangxi (P = 0.02).Nearly 75% of participants with CIMT measurement had at least 1 carotid plaque (P < 0.01).Sex and province-specific demographics are provided in Supplementary Table S2 online.During the study period, 7 participants reported diagnoses of incident coronary heart disease and 5 reported having a stroke.
Exposure to solid fuel stove use and PM 2. 5   Most participants cooked with solid fuel (65%), and 82% of northern participants heated with solid fuel (Table 3).On average, participants reported 249 (SD: 214) indoor solid fuel "stove-use days" in the previous year, with greater intensity of use in Shanxi and Guangxi.Long-term (previous 20 years) intensity of indoor solid fuel stove use was greatest in Shanxi [mean (SD): 18.9 stove-use years (6.5) vs. 15.5 (8.5) in Beijing and 14.4 (8.5) in Guangxi].Average estimated yearly exposures to PM 2.5 were 97.5 µg/m 3 , which is ~10 times the WHO's PM 2.5 guideline (10 µg/m 3 ).We observed a moderately high correlation between seasonal and estimated yearly average exposure to PM 2.5 among participants with air pollution exposure assessment in more than 1 campaign (r = 0.78).

Associations of solid fuel use with blood pressure and vascular function
Current use of solid fuel cookstoves was associated with higher BP (systolic: 2.4 mm Hg, 95% CI: −0.4,4.9; diastolic: 1.4 mm Hg, 95% CI: −0.1, 3.0) compared with exclusive clean fuel cooking (Figure 1), though both confidence intervals included zero.We did not find consistent associations between current use of solid fuel heaters and intensity of stove use with BP or between fuel use and arterial stiffness, except for an inverse association between current use of solid fuel heaters and lower bfPWV (−1.5 m/s, 95% CI: −3.0, 0.0).Current and long-term solid fuel use and greater intensity of use were all associated with greater CIMT and features of plaques, but the coefficients were small and confidence intervals included zero (Figure 2).

Sensitivity analyses
Limiting our analysis of current stove use to only indoor stoves attenuated the associations between stove use and BP, likely due to greater exposure misclassification (Supplementary Table S3 online).Additional sensitivity analyses did not change our overall findings, including adjustment for hypertension status or excluding those with diabetes, both which had little effect on our point estimates.
Additional adjusting for heating fuel in models with cooking fuel as the exposure and vice versa also did not change our overall findings (Supplementary Table S4 online).

DISCUSSION
In this multi-provincial study of Chinese adults, exposure to household air pollution was associated with higher BP and with greater CIMT and total area of carotid plaques.The associations between household air pollution and arterial stiffness were mostly null, with the exception of lower bfPWV in users of solid fuel heaters.
Our study has a number of important strengths.Our relatively large sample included men, who comprise 35% of household air pollution-attributable cardiovascular diseaserelated burden 20 but are understudied. 21We controlled for a comprehensive set of confounders including socioeconomic status and dietary indicators that were not considered in Abbreviations: CCA, common carotid artery; ECA, external carotid artery; ICA, internal carotid artery; PM 2.5 , fine particulate matter <2.5 µm.Survey data indicated no major change in fuel or stove use between campaigns, such that the same fuel use exposure variables were used for each season in the repeated measures analysis.previous studies. 11,12,14We measured up to 96 hours of exposure to PM 2.5 across seasons, allowing us to estimate yearly averages, which are relevant for long-term outcomes like CIMT and plaques.Finally, our outcomes reflect different domains of functional and structural vascular damage, and are recognized risk factors for cardiovascular and renal disease.S3 online.S3 online.
][27][28] Household air pollution was not associated with arterial stiffness in our study, with the exception of an unexpected finding of lower bfPWV among users of solid fuel heaters.Previously, our group found a null association between household air pollution and arterial stiffness, but outdoor air pollution has been linked with vascular damage through ventricular repolarization-induced myocardial vulnerability and arrhythmias. 25,29A possible explanation is residual confounding by home (indoor) temperature, which is associated with BP and arterial stiffness. 30he PM 2.5 -CIMT coefficient [0.02 mm per 1 ln(µg/m 3 ) increase in PM 2.5 ] in our study is similar to CIMT differences between biomass and clean fuel users in Peru (0.03 mm) 31 and Nigeria (0.04 mm). 6Our CIMT results are supported by our finding that PM 2.5 and solid fuel stove use were associated with a greater total plaque area, which are consistent with a study of biomass use and prevalence of atherosclerotic plaque in Peru. 31Outdoor PM 2.5 has also been linked with greater carotid plaque (HR: 1.78, 95% CI: 1.55, 2.03) in a cohort of previously carotid plaque-free individuals living in Shijiazhuang, China. 32Our study provides a new dimension to the existing literature on air pollution and atherosclerosis, which is mostly comprised studies in European and North American populations exposed to much lower air pollution (e.g., 0.02 mm per 10 µg/m 3 , with PM 2.5 ranging from 13 to 23 µg/m 3 ). 33ur finding linking lower GSM, an indicator of plaque instability and rupture risk, 17 with higher PM 2.5 in older participants is novel but should be interpreted with caution.Greater than 90% of the women in our study sample were postmenopausal and are at greater risk for vascular remodeling. 34As such, older participants (≥63 years) were more likely to have plaques than younger participants (80% vs. 69%), though no significant interactions between longer exposures and age were observed.Plaques can become more vulnerable to rupture after exposure to PM 2.5 via blood thrombogenicity and the promotion of foam cell formation. 35However, this finding was not supported by our fuel use analysis, and our ultrasound images were not designed for advanced plaque tissue characterization.
Differences in underlying health and cardiovascular risk profiles of participants from northern vs. southern China could explain the larger associations in our northern participants.Northern China has higher overall cardiovascular risk including stroke. 36Differences in fuel types and PM 2.5 composition may also differentially affect hypertension and vascular function of participants in northern vs. southern China. 37Biomass stoves were used in all sites, whereas coal was more common in the north.However, we could not differentiate between fuel types in our analysis as most northern households using biomass also used coal.Notably, restricting our fuel use exposure measure to only indoor stoves attenuated the associations between stove use and BP.Air exchanges between the outdoor and indoor environments are high in rural China 38 such that pollution from outdoor stoves could still impact personal exposures.It is possible that excluding outdoor stoves introduced exposure misclassification which led to bias toward the null.
Based on our multivariable models, the difference in BP and vascular outcomes for participants in the 90th vs. 10th percentile of PM 2.5 exposure are 2.5 mm Hg systolic BP, 1.8 mm Hg diastolic BP, 0.03 mm CIMT, and 6.6 mm 2 total area of plaques.Although associations of this magnitude are relatively small, they translate into substantial clinically meaningful reduction in BP at the population level.For example, 10 mm Hg higher systolic BP is associated with a 21% increase in cardiovascular mortality in Chinese adults. 39Even a 2 mm Hg reduction in systolic BP can reduce coronary heart disease by 4%, and stroke mortality by 6%, 40 and a 2 mm Hg reduction in diastolic BP is associated with a 14% reduction in stroke risk, and 6% reduction in coronary heart disease in 35-64-year US men. 41ur study did not capture the change in various indicators of vascular damage due to household air pollution given the relatively short follow-up (<9 months) between each campaign, but a 0.1 mm decrease in CIMT translates to a 15% decrease in risk myocardial infraction and a 17% decrease in stroke. 42ur study has several limitations to consider for future studies.Stove use was self-reported and retrospectively assessed, and thus subject to recall bias.The mostly null findings for long-term solid fuel use, in contrast to more consistent associations with current practices, may be due to exposure misclassification.Measurement error in personal exposure to PM 2.5 is also likely given the large withinindividual variability in exposure among our participants (σ 2 = 0.80), 14 though this error most likely nondifferential and would bias our results toward the null.Residual confounding is possible in this observational study.Our repeated-measures design accounted for nonvarying factors, and we measured most known BP and vascular risk factors, including behavioral and socioeconomic factors not previously measured.Though we were unable to measure home temperature due to study logistics, which may have confounded some wintertime associations.Further, the findings from our repeated measures analysis should be interpreted with caution as we obtained only one measurement of some time-varying covariates like lipids and did not collect information on indoor temperature or BP medication compliance.Finally, given that our study is an observational and cross-sectional study, it is limited in its ability to identify causal effects.Future studies in this cohort can assess if household air pollution affects the development of hypertension or the progression of atherosclerosis and plaque instability.
These findings add to a small but rapidly growing body of evidence that household air pollution is associated with subclinical cardiovascular outcomes including higher blood pressure, and our study contributes to the relatively limited exposure-response data on personal exposure to household air pollution and the severity of atherosclerosis and vascular disease, particularly for men.China's cardiovascular burden exceeds the global average, and is expected to increase by 3.4 million excess deaths by 2030 as its population ages. 43eduction of elevated BP can effectively slow the progression of end-stage renal and cardiovascular disease. 44Population interventions that reduce BP and the prevalence of hypertension are urgently needed, and our study contributes to evidence that clean household energy could be part of a comprehensive hypertension prevention strategy.

Figure 1 .
Figure1.Associations of household fuel use and personal exposure to PM 2.5 with blood pressure and arterial stiffness.Reference groups: 1 = current exclusive use of clean fuel cookstoves; 2 = current exclusive use of clean fuel heating stoves.Mixed effects regression models with participant-and village-specific random intercepts adjusting for age, sex, income, waist circumference, height (AIx, bfPWV), smoking status, secondhand smoke exposure, alcohol consumption, physical activity, 24-h urinary sodium, time of day of BP measurement (PM 2.5 models only), and total/HDL cholesterol ratio were used.Results not shown for no heating stoves group (n = 19).Personal exposure to PM 2.5 refers to the average season-specific exposures to PM 2.5 over 2 consecutive days during each campaign.Point estimates are available on Supplementary TableS3online.

Figure 2 .
Figure 2. Associations of household fuel use and personal exposure to PM 2.5 with atherosclerosis and plaques.Reference groups: 1 = current exclusive use of clean fuel cookstoves; 2 = current exclusive use of clean fuel heating stoves.Mixed effects regressions models with village-specific random intercepts adjusting for age, sex, yearly household income, waist circumference, alcohol consumption, smoking status, secondhand smoke exposure, physical activity, and total/HDL cholesterol ratio were used.Personal exposure to PM 2.5 refers to estimated yearly exposure based on time-weighted averages of season-specific measurements of personal exposure to PM 2.5 (7 nonheating and 5 heating months).Point estimates are available on Supplementary TableS3online.

Table 1 .
Household air pollution exposure variables and their descriptions a

Table 2 .
Characteristics of study participants by current cooking fuel use (n (%) or mean (SD)) ≥90 mm Hg) blood pressure; Physical activity = exercise and farm-based physical activity domains.Characteristics of our study population are provided by province and sex in SupplementaryTable S2 online, and variables with missing data are provided in Supplementary Table S1 online.

Table 3 .
Descriptive statistics for blood pressure and vascular outcomes and exposure to household air pollution, by study site and campaign (mean (SD), median [min, max], or n (%))