Introduction

Together with population aging, frailty has garnered considerable interest among clinical researchers and public health professionals, because its reversible and dynamic nature offer promise for disability prevention. Frailty is a geriatric syndrome characterized by declining physiological reserves and increasing vulnerability, in which disrupted homeostasis due to accumulated inter-systemic deficits results in disproportionate health consequences in response to minor stressors1. This complex age-related pathophysiology is considerably influenced by combined environmental and genetic factors2.

Air pollution is the foremost reversible environmental factor associated with premature death or disability3. Accruing epidemiological evidence supports associations between exposure to air pollution and adverse outcomes such as cardiovascular events4,5,6, cancer7, Alzheimer’s disease8 and mortality9. Fine particulate matter < 2.5 μm in diameter (PM2.5) is of greatest health concern, not only such small particles can penetrate deeper into the lungs than larger ones do, but also they are likely to be made up of a more toxic mix4.

Possible biological mechanisms explaining the health consequences of PM2.5 exposure, include oxidative stress, endothelial dysfunction, and pro-inflammatory effects such as elevated interleukin 6 levels4,10; frailty shares similar pathoetiology11,12. Previous studies have investigated potential associations between fine airborne particulates and frailty13,14,15. A study of 2059 non-smokers from the United States National Health and Nutrition Examination Survey found a positive association between secondhand tobacco smoke and frailty13. Others confirmed an association between PM2.5 and frailty, measured using a frailty index, among post-myocardial infarction patients who were inherently vulnerable to frailty14,15. However, it remains unknown whether or not PM2.5 exposure is associated with frailty in the general population. Hence, we investigated associations between annual PM2.5 exposure and frailty among community-dwelling older people, and in potentially vulnerable subgroups.

Results

The prevalence of frailty in this study sample was 5.2%; compared to the robust people, 1,080 frail participants were significantly older, predominantly female, had higher disease burden, lived rurally, and had lower annual PM2.5 exposure (Table 1). Adjusted for age, sex, smoking, CCI, and urbanicity, frail and pre-frail statuses were significantly associated with high annual PM2.5 exposure. Compared to robust older adults, Odds ratio (OR) for those with prefrail and frailty were 1.4 (95% confidence interval (CI) 1.3–1.5) and 1.5 (95% CI 1.2–0.9). Weakness (OR 1.4, 95%CI 1.3–1.6) and low activity (OR 1.6, 95%CI 1.4–1.7) were associated with high PM2.5 exposure, but other phenotypic components of frailty were not (Table 2).

Table 1 Demographic and health-related characteristics by frailty status.
Table 2 Association between frailty status and individual phenotypes, and PM2.5 in logistic regression analyses.

Figure 1 and Supplementary Table S1 in the supplement summarize the results of multivariate multinomial logistic analysis of associations between frailty and PM2.5 across various subgroups. The likelihood of people with high PM2.5 exposure being frail was generally higher among vulnerable groups such as those who were older, had higher disease burden, or with low household income. Since the lowest PM2.5 exposure level in rural areas was zero, we used subsamples from urban areas to test associations between frailty and PM2.5, with results similar to those for the entire sample.

Figure 1
figure 1

Association between high PM2.5 exposure and frailty status by age, sex, smoking, disease burden, and household income.

Frailty reflects aging and the progression of disease, including comorbid chronic conditions. The association between frailty and PM2.5 seems to be driven by weakness and low activity—which presumably were also associated with chronic disease (respiratory disease, cardiovascular disease, cancer, etc.) and were also likely associated with PM2.5. Hence, we conducted a sensitive analysis that adjusted for chronic diseases of hypertension, diabetes, heart disease, stroke and chronic kidney disease. Results were similar (Table supplement S2).

Discussion

This study discovered a significant association between frailty and fine airborne particulates in a large population of community-dwelling elders. The associations were stronger among men, people older than 75 years, with higher disease burden, non-smokers, and those with lower household income.

Most studies on air pollution have focused on conventional disease-oriented outcomes, such as cardiometabolic disease, cancer, or mortality5,7,9, rather than on functional ones; few studied frailty outcomes. However, frailty per se is emerging as a prime candidate for disability prevention in the modern era of function-centered medicine. The role of frailty in modifying pulmonary function deterioration16 and post-myocardial infarction mortality15 has been reported. Ten-year follow-up of 848 non-frail myocardial infarction patients showed increased likelihood of becoming frail14. Hong Kong investigators found differences in frailty prevalence between administrative districts but failed to prove an association with air pollution, due to limited exposure data17. Our findings affirm an association between PM2.5 and frailty, and extend this from heart attack survivors to the general population.

Phenotypic manifestations of frailty are better suited to investigating pathological mechanisms1 than are operational definitions based on a frailty index18. As far as we know, this is the first evidence to of an association between PM2.5 exposure and phenotypically-defined frailty in the general population. Among five established frailty phenotypes, weakness was significantly associated with PM2.5 exposure levels. Weakness, which is an essential component of mobility type frailty and sarcopenia19,20, increased accuracy in predicting cardiovascular mortality in previous studies21. Given the strong association between handgrip strength and mortality, weakness may be an intermediate state on the way to eventual mortality due to air pollution.

Frailty, disease burden, and socioeconomic status have synergistic effects on disability22,23. We found the likelihood of people with high PM2.5 exposure levels being frail to be around four–five-fold higher among vulnerable groups, such as those with lower household income and higher disease burden, consistent with previous studies14,15. Likewise, the association was stronger in people aged ≥ 75 years compared to younger ones, corroborating other evidence that PM2.5 exposure has a greater impact in older people15. Taken together, the results of this study suggest that a stronger association between PM2.5 exposure and frailty in vulnerable groups may be particularly inimical to health outcomes.

Although frailty is more prevalent in women, the association between frailty and PM2.5 exposure was only evident in men, even after adjusting for smoking. A time-series study in Paris found that men, but not women, had significantly increased risk of hospitalization for respiratory complaints associated with air pollution24; this was consistent with the sex difference we observed, but the mechanism is unknown and warrants further investigation.

Non-smokers with high PM2.5 exposure in our study population were more likely than smokers to be frail. This may reflect an unmeasured residual competing risk of active smoking that masked the health hazard of PM2.5. Misclassification bias might be another plausibility. Similarly, the European Study of Cohorts for Air Pollution Effects found a stronger association between air pollution and stroke in non-smokers than in smokers6.

Putative mechanisms explaining the health consequences of PM2.5 exposure include oxidative stress, inflammatory responses, and gene/epigenetic modifications—factors also implicated in aging and frailty pathogenesis1. Airborne pollutants may disrupt homeostasis and accelerate age-related declines in functional performance and capacity at cellular, organ and system levels25; multi-systemic homeostatic disruption would conduce to frailty and hinder successful aging. Previous studies that investigated the relationship between PM2.5 exposure and frailty, used an operational definition based on a frailty index which measures cumulative deficits18 rather than making direct phenotypic assessments11. Our study affirmed an association with directly-measured phenotypic frailty.

Although the mean annual PM2.5 level in New Taipei city was much lower than in many other parts of Asia, which routinely exceed 35 μg/m3,4 it is considerably higher than the optimal target of < 10 μg/m3 set by World Health Organization air quality guidelines26 and routine values of < 12 μg/m3 in the United States and Canada4. Findings from a Hong Kong study of 3,240 community-dwelling adults ≥ 65 years old, suggests that green neighborhood space may mitigate exposure to air pollution and reverse frailty27, which could be the basis of a national level strategy. Personal level monitoring may contribute to increased awareness and protection4.

This study had limitations. First, the cross-sectional design precluded exploration of a reverse causality between frailty and PM2.5 exposure. Second, air pollution data were obtained from ambient air monitoring stations rather than personal exposure levels, which prevented distinguishing between indoor and outdoor pollution. Third, urbanicity was determined according to an urbanization index developed a decade ago, which may not accurately represent shifted demographic characteristics.

Conclusion

Our study detected a significant association between PM2.5 exposure and frailty in among the general population, which was stronger in vulnerable groups. These results are consistent with the hypothesis that exposure to PM2.5 disrupts multi-systemic homeostasis, suggesting that frailty prevention and intervention strategies should incorporate reducing air pollution.

Methods

Data sources and participants

Study data were excerpted from the New Taipei City Elderly Health Examination Database (NTCHD) and the Taiwan Air Quality-Monitoring Database (TAQMD). The NTCHD was established for early detection of physical conditions and promote senior health. All older adults ≥ 65 years old residing in New Taipei City could receive voluntary government-funded annual examinations face-to-face by physicians, and these examinations include anthropometry, health-related behaviors, physical and mental performance and biochemistry results.

Upon enrollment, participants gave written informed consent authorizing the New Taipei City Government to process health examination data for research and policy purposes. New Taipei City Department of Health removed all potentially identifying information to protect privacy and generate the anonymized NTCHD. Full details of recruitment and data collections procedures are described elsewhere22.

The present study included 26,026 potentially eligible NTCHD participants in 2016. Having excluded 5,420 with incomplete data, the residential areas of the final analytic sample of 20,606 were linked to the locations of TAQMD ambient air-quality monitoring stations to estimate daily PM2.5 exposures.

This study was designed and conducted in accordance with the principles of the Declaration of Helsinki. New Taipei City Department of Health approved the use of this anonymized dataset for research purposes, and waived the requirement for Institutional Review Board approval. The design and reporting format follow Strengthening the Reporting of Observational Studies in Epidemiology (STROBE) guidelines.

Frailty

The Cardiovascular Health Study frailty criteria comprise unintentional weight loss, exhaustion, weakness, slowness, and low activity; people with three or more of these frailty phenotypes are classed as frail, those with one or two as prefrail, and those with none as robust11. Unintentional weight loss is defined as losing > 5% of body weight over the previous 12 months. Exhaustion was determined by affirmation of two questions from the Center for Epidemiologic Studies Depression Scale questionnaire28—‘‘I felt everything I did was an effort’’ and ‘‘I could not get going’’—on ≥ 3 days per week. Physical activity was calculated from the International Physical Activity Questionnaire short-form score, based on self-reported exercise and leisure time physical activities, and expressed as weekly energy expenditure29. Energy expenditure was calculated as the metabolic equivalent relevant to designated type of exercise times body weight and days of exercise. In this study, energy expenditure below 383 kcal/week in men or 270 kcal/week in women was defined as low physical activity according to Cardiovascular Health Study11. Walking speed was measured by a six-metre walk test at a normal pace from a moving start without deceleration. Weakness was defined as maximum dominant handgrip strength of < 26 kg for men or < 18 kg for women, and slowness as 6-m walk speed < 0.8 m/s30.

Urbanization

An urbanization index developed by the Taiwan National Health Research Institute was used to determine the urbanization level of each participant’s residential area. This index, which was based on national census data including population density and aging, education levels, medical resources, and agricultural employment, classified 359 regions throughout Taiwan into seven strata from most urbanized (level 1) to the least (level 7)31. Urbanicity was dichotomized as urban (level ≤ 3) or rural (level ≥ 4), based on a previous study32.

Other variables

Smoking was defined as using tobacco during the previous 6 months. Charlson Comorbidity Index (CCI) ≥ 2 indicated severe disease burden33. Household income data from the Taiwan Ministry of Finance was used as a proxy for socioeconomic status; median annual income equivalent to less than the population median of 19,677 United States Dollars (USD) was considered low.

PM2.5 exposure

Daily PM2.5 level was excerpted from TAQMD. Average annual fine particulate matter concentration was calculated by accumulating daily PM2.5 levels from the 12 months prior to the index interview date; exposure in the lowest quartile (< 16.3 μg/m3) was classed as low.

Statistical analysis

All analyses were performed using SAS statistics software, Version 9.4 for Windows (SAS Institute, Inc., Cary, NC, USA). A two-sided p-value of < 0.05 or 95% confidence interval (CI) that did not include the null hypothesis value were considered statistically significant.

The Kolmogorov–Smirnov D test was used to check whether numerical variables were normally distributed. Continuous variables with normal distributions were expressed as means plus/minus standard deviation, and those with non-normal distributions as median (first quartile, third quartile). Categorical variables were expressed as frequency/proportions. One-way ANOVA, Kruskal–Wallis Test, chi-square analysis, or Fisher exact test were used as appropriate to compare descriptive characteristics. Cochran-Armitage trend test was used to test for trends. PM2.5 concentrations were transformed logarithmically to approximately linearity assumption for logistic regression. Urbanicity exhibited collinearity with household income, hence income was not included in the statistical model. Univariable and multivariable multinomial logistic regression analyses, adjusted for age, sex, smoking, CCI, and urbanicity, were used to explore associations between frailty status and PM2.5. Multivariable multinomial logistic analysis employed statistical analysis on association between frailty and PM2.5. Frailty as a categorical variable severed as a dependent variable; PM2.5 and other confounders referred as independent variable. The statistical results showed ORs and CI comparing prefrail vs. robust and frail vs. robust, respectively. Log PM2.5 was input to the model after categorizing low/high PM2.5 for sensitivity analysis. Prespecified subgroup analyses included age (< 75 vs. ≥ 75 years), sex (male vs. female), disease burden (CCI < 2 vs. ≥ 2), smoking (no vs. yes), household income (≥ 19,677 USD vs. < 19,677 USD) and urbanicity (urbanization index ≤ 3 vs. ≥ 4).