WU Jun Hui , WU Yao , WANG Zi Jing , TIAN Yao Hua , WU Yi Qun , WU Tao , WANG Meng Ying,WANG Xiao Wen, WANG Jia Ting, and HU Yong Hua,2,#
Cardiac arrhythmia is a serious public health problem in many countries[1].Previous studies estimated that 33.5 million people are affected by cardiac arrhythmia worldwide,and this number will continue to grow as society ages[2]. Despite improvements in diagnostic and therapeutic interventions in electrophysiology,the disease burden,incidence and prevalence of cardiac arrhythmia continue to increase and have aroused public health concern.Increasing evidence has indicated that air pollution may be associated with cardiac autonomic nervous system[3].Among air pollutants,particulate matter with an aerodynamic diameter of ≤ 2.5 μm (PM2.5) is considered exert more toxicity than other air pollutants,as it provides a larger surface area and absorbs or condenses more toxic substances per unit mass[4].
However,most of previous studies were performed in western countries where PM2.5concentrations are much lower than those in developing countries[5].Furthermore,these studies focused primarily on hospital admissions,which may affected by many factors like scheduled appointments and availability of hospital beds.Thus,we use a time-series design to evaluate whether PM2.5pollution levels were related to the risk of cardiac arrhythmia between 2010 and 2012 in Beijing,China.To our knowledge,this is the first time that the short-term effects of PM2.5levels on arrhythmia-associated outpatient visits,hospital admissions,and emergency admissions have been studied simultaneously at the city level.
We obtained the daily numbers of hospital visits for cardiac arrhythmia from the citywide database,Beijing Medical Claim Data for Employees (BMCDE),which records medical claim data for all working or retired employees (aged ≥ 18 years) who are covered by basic medical insurance in Beijing.By the end of 2017,nearly 80% of Beijing’s residents (17.8 million)were included in the database.Basic characteristics(e.g.,sex and age) of individuals diagnosed with cardiac arrhythmia were obtained from the BMCDE between January 1,2010 and June 30,2012 (912 days).As in previous studies,cardiac arrhythmia was identified by ICD-10 codes (Supplementary Material available in www.besjournal.com)[3].Patients younger than 18 years old were excluded.
Daily PM2.5concentration data were retrieved from an air-monitoring station of the US Embassy(http://www.stateair.net) in Chaoyang District,Beijing.China did not include PM2.5in the national first-report monitoring data until 2013;thus,the monitoring data from the US Embassy were the only open source of PM2.5concentration data during the study period.Previous studies have reported that the PM2.5levels detected by the US Embassy were approximately comparable with citywide PM2.5levels in Beijing[6].In order to reduce exposure misclassification,it could be considered that the maximum distance from the US Embassy to the hospital is about 40 km[7].Within that radius,the monitoring data covering all the region with high population density in Beijing (> 5,000 residents/km2) and 79.2% of Beijing's residents[6].In the present study,we used the daily average PM2.5concentrations (24 h) as a proxy for population exposure levels.Moreover,daily meteorological data including temperature (°C) and relative humidity (%)were collected from the Chinese Meteorological Bureau.
The association between PM2.5concentrations and cardiac arrhythmia-associated hospital visits was assessed using a quasi-Poisson regression model,which has been widely adopted in time-series studies[8].The following covariates were included in the main model:a penalized spline function of calendar time with 7 degrees of freedom (df) per year to control for underlying time trends,public holiday,and day of the week as categorical variables to account for short-term variations,and penalized spline functions with 3dffor temperature and 3dffor relative humidity.The formula applied in the present study is as follows:

whereE(Yt) refers to the expected number of cardiac arrhythmia-associated hospital visits on dayt,α represents the model intercept,β denotes the log(relative risk) of morbidity relative to unit increase in PM2.5level;andsrepresents a smoother based on the penalized splines.Following the approaches applied in previous studies,we selected the degrees of freedom (df) for calendar time,temperature,and relative humidity[3].
We applied a penalized cubic regression spline of the PM2.5concentration with 3dfto assess the concentration-response association.We also explored the relationship between PM2.5and cardiac arrhythmia-associated hospital visits by building models with a single-day lag from the current day (lag 0) up to the previous 3 days (lag 1,lag 2,and lag 3),and with 2-day (lag 0-1),3-day (lag 0-2),and 4-day (lag 0-3) moving average concentrations.To examine potential effects per subgroup,analyses were stratified by age,sex,and season.The warm season was defined as April to September,and the cool season was defined as October to March.We used the Z-test to assess the statistical significance of the subgroup differences.
All results were presented as the percentage change and 95%CIof the daily cardiac arrhythmiaassociated hospital visits for each 10-μg/m3increase in ambient PM2.5.We applied the“mgcv”and“nlme”packages in R 3.2.2 for all the analyses[9].
In total,1,435,139 outpatient visits,29,837 hospital admissions,and 43,347 emergency visits between January 1,2010 to June 30,2012 were identified from the BMCDE database.Table 1 summarizes the descriptive statistics for the cardiac arrhythmia-associated hospital visits,PM2.5concentrations,and weather conditions.The daily mean counts for the outpatient visits,hospital admissions and emergency visits were 1,573 (1,132),33 (22),and 48 (27) respectively. For PM2.5concentration,the annual average was 99.5(75.3) μg/m3,and the maximum was 493.0 μg/m3.
Table 2 presents the acute elevations in PM2.5concentrations were related to increased hospitalvisits for cardiac arrhythmia.A 10-μg/m3increase in PM2.5level on lag days 0-3 corresponded with 0.71%(95%CI: 0.43%-0.47%) and 0.48% (95%CI:0.59%-0.82%) increases in outpatient and hospital admissions for cardiac arrhythmia,respectively.A 10 μg/m3increase in PM2.5level on current day corresponded with 0.17% (95%CI:0.10%-0.23%)increases in emergency visits for cardiac arrhythmia.Our findings provide new evidence which may help improve targeted intervention strategies for cardiac arrhythmia in China.Previous western studies have explored associations between PM2.5and hospitalizations for cardiac arrhythmia,which support the results of our study.A recent metaanalysis of PM2.5and daily hospitalizations for cardiac arrhythmia reported that arrhythmiaassociated hospitalizations were related to increases in PM2.5(relative risk=1.02 per 10 μg/m3)[5].However,most of previous studies have been conducted in the US and Europe,while Beijing has higher PM2.5levels than those areas,which might lead to differences in effect estimates.

Table 1.Summary statistics for daily cardiac arrhythmia-associated hospital visits,daily PM2.5 concentrations and meteorological data
Table 3 lists the estimates for subgroup analyses.The association between PM2.5and arrhythmiaassociated outpatient visits was greater during the warm season (0.94%,95%CI:0.91%-0.97%),and the estimated PM2.5effect was stronger in men (0.82%,95%CI:0.66%-0.99%).For hospital admissions,the estimates were higher for men (0.48%,95%CI:0.45%-0.51%) and during the warm season (0.52%,95%CI:0.50%-0.54%).For emergency visits,the estimates were greater for men (0.19%,95%CI:0.10%-0.29%) and during the warm season (0.22%,95%CI:0.14%-0.30%).Similarly,a recent study reported that the short-term effect of PM2.5on arrhythmias was relatively stronger for men than for women in all lag models[10],but the underlying mechanism of the sex difference remains elusive and requires more study. Differences in seasonal estimates may be due to the behavioral differences.Beijing residents are more likely to participate in outdoor activities and ventilate their homes duringthe warm season,bringing monitoring data closer to their individual exposure levels.Another explanation for the higher correlation is that high temperatures may affect the chemical conversion and deposition processes of air pollutants.

Table 2.Percentage change with 95% confidence interval in cardiac arrhythmia-associated hospital visits for each 10 μg/m3 increase in PM2.5 levels for different lag days

Table 3.Increases in cardiac arrhythmia-associated hospital visits for each 10 μg/m3 increase in PM2.5 for lag days 0-3 (for outpatient visits and hospital admissions)and current-day (for emergency visits) with stratification
In the present citywide time-series analysis,significant and positive associations between PM2.5levels and hospital visits for cardiac arrhythmia were found in Beijing.To our knowledge,this is the first citywide research in Beijing,or even in other highly air polluted region,to comprehensively reveal the acute effect of PM2.5on cardiac arrhythmiaassociated outpatient,emergency,and hospital admissions in a same study.
Some limitations should be acknowledged.First,data from one fixed monitoring station could notexactly represent the actual situation of personal exposure to PM2.5.Second,because the information was limited,we did not explore other potential modifiers,including health-related and nutritional factors,which are reported to be associated with arrhythmias[3].Third,the US Embassy only provided PM2.5monitoring data during the study period.The lack of authoritative records of other air pollutants has limited our exploration of the independent effects of PM2.5. Finally,due to information limitations,we are unable to distinguish all subtypes of arrhythmias.Future research is warranted to investigate the modification effect between different cardiac arrhythmia subtypes.
#Correspondence should be addressed to Professor HU Yong Hua,PHD,MD,Tel:86-10-82801189,E-mail:yhhu@bjmu.edu.cn
Biographical note of the first author:WU Jun Hui,female,born in 1992,PHD Candidate,majoring in genetic epidemiology.
Supplementary Material
We used ICD-10 codes to identify cardiac arrhythmia and these codes included I44 (atrioventricular and left bundle-branch block),I45 (other conduction disorders),I46 (cardiac arrest),I47 (paroxysmal tachycardia),I48(atrial fibrillation and flutter),and I49 (other cardiac arrhythmia)(Dominici et al.2006[1],Zheng et al.2018[2]).
Reference
1.Dominici F,Peng RD,Bell ML,et al.Fine particulate air pollution and hospital admission for cardiovascular and respiratory diseases.JAMA,2006;295,1127-34.
2.Zheng Q,Liu H,Zhang J,et al.The effect of ambient particle matters on hospital admissions for cardiac arrhythmia:a multi-city casecrossover study in China.Environ Health,2018;17,60.
Biomedical and Environmental Sciences2021年7期