Int J Med Sci 2026; 23(5):1761-1767. doi:10.7150/ijms.126734 This issue Cite

Research Paper

Synergistic Effect of Ambient PM2.5 Exposure and Waist-Hip Ratio on Non-Alcoholic Fatty Liver Disease Risk in a Taiwanese Population

Shih-Tsung Chang1,2, Ying-Hsiang Chou2,3,4, Chih-Da Wu5,6, Oswald Ndi Nfor7, Wen Yu Lu7, Chien-Ning Huang1,8 Corresponding address, Yung-Po Liaw1,7,9,10,11 Corresponding address

1. Institute of Medicine, Chung Shan Medical University, Taichung City, 40201, Taiwan.
2. Department of Radiation Oncology, Chung Shan Medical University Hospital, Taichung City, 40201, Taiwan.
3. School of Medicine, Chung Shan Medical University, Taichung City,40201, Taiwan.
4. School of Medical Imaging and Radiological Sciences, Chung Shan Medical University, Taichung City, 40201, Taiwan.
5. National Institute of Environmental Health Sciences, National Health Research Institutes, Miaoli, Taiwan.
6. Department of Geomatics, National Cheng Kung University, Tainan, Taiwan.
7. Department of Public Health, Institute of Public Health, Chung Shan Medical University, Taichung City 40201, Taiwan.
8. Department of Internal Medicine, Chung Shan Medical University Hospital, Taichung City, 40201, Taiwan.
9. Department of Medical Imaging, Chung Shan Medical University Hospital, Taichung City 40201, Taiwan.
10. Department of Public Health, College of Public Health, China Medical University, 406040, No. 100, Sec. 1, Jingmao Rd., Beitun Dist., Taichung, Taiwan.
11. Department of Medical Research, China Medical University Hospital, China Medical University, Taichung City, Taiwan.

Received 2025-10-14; Accepted 2026-3-24; Published 2026-3-30

Citation:
Chang ST, Chou YH, Wu CD, Nfor ON, Lu WY, Huang CN, Liaw YP. Synergistic Effect of Ambient PM2.5 Exposure and Waist-Hip Ratio on Non-Alcoholic Fatty Liver Disease Risk in a Taiwanese Population. Int J Med Sci 2026; 23(5):1761-1767. doi:10.7150/ijms.126734. https://www.medsci.org/v23p1761.htm
Other styles

File import instruction

Abstract

Graphic abstract

Epidemiological studies have established a connection between environmental factors and the prevalence of non-alcoholic fatty liver disease (NAFLD). Fine particulate matter with an aerodynamic diameter ≤2.5 μm (PM2.5), a major component of ambient air pollution, has been implicated in systemic inflammation and metabolic dysfunction. The waist-hip ratio (WHR), an important measure of body composition, is a key risk factor for metabolic syndrome and NAFLD. We investigated the association between ambient PM2.5 exposure, WHR, and NAFLD. We used data from the Taiwan Biobank collected between 2016 and 2020, involving 15,049 individuals aged 30 to 70 years. Ambient PM2.5 exposure was estimated using a land-use regression model and categorized into quartiles. NAFLD was identified using liver function tests and imaging, and WHR was classified as normal or abnormal based on sex-specific cutoffs. Participants with NAFLD (n = 7179) had higher body mass index, a greater prevalence of abnormal WHR, and less favorable lipid profiles compared with those without NAFLD (n = 7870). Higher PM2.5 exposure was associated with increased odds of NAFLD, with participants in the highest quartile (PM2.5 > 39.85 μg/m³) showing an odds ratio (OR) of 1.251 (95% CI: 1.114-1.404). An abnormal WHR was also associated with higher odds of NAFLD (OR = 1.641, 95% CI: 1.507-1.787). A statistically significant interaction between PM2.5 exposure and WHR was observed, with individuals with abnormal WHR exhibiting consistently higher odds of NAFLD across all PM2.5 quartiles. In the highest PM2.5 quartile, participants with abnormal WHR had an OR of 2.020 (95% CI: 1.708-2.389), compared with an OR of 1.442 (95% CI: 1.215-1.712) among those with normal WHR. These findings suggest that central obesity may amplify the association between PM2.5 exposure and NAFLD, highlighting the importance of integrated public health strategies targeting both air pollution and metabolic health.

Keywords: PM2.5 exposure, waist-hip ratio, non-alcoholic fatty liver disease, land-use regression model

Introduction

Air pollution and obesity are significant public health issues that have a profound impact on metabolic disorders worldwide. PM2.5 is particularly harmful because its small size allows it to penetrate deeply into the lungs and enter the systemic circulation, thereby contributing to respiratory and cardiovascular diseases. Increasing evidence also suggests that PM2.5 exposure may play an essential role in metabolic dysregulation, including the development of NAFLD [1, 2].

Epidemiological studies have increasingly linked environmental factors to the rising global prevalence of NAFLD. Among these factors, PM2.5 has gained attention due to its adverse effects on liver health. VoPham et al. (2022) [3] and Guo et al. (2022) [4] reported that PM2.5 exposure is associated with an increased risk of metabolic dysfunction-associated fatty liver disease (MAFLD), potentially mediated through systemic inflammation and oxidative stress that promote hepatic steatosis and fibrosis. Additional epidemiological evidence further supports the association between PM2.5 exposure and NAFLD risk [5].

Given that PM2.5 exposure contributes to metabolic dysfunction, individual susceptibility factors—particularly those related to body fat distribution—may modify its impact on liver disease risk. WHR, an important indicator of central adiposity, is a well-established risk factor for metabolic syndrome and NAFLD [6-8]. Emerging evidence suggests that central obesity may exacerbate the deleterious effects of PM2.5 on liver health. Excess visceral fat promotes oxidative stress, chronic inflammation, and adipocytokine imbalance, mechanisms that may amplify PM2.5-induced hepatic injury and accelerate NAFLD development [9, 10].

In Taiwan, where the prevalence of NAFLD has been steadily increasing [11], the Taiwan Biobank (TWB) provides a valuable population-based resource for examining the interplay between environmental exposures, anthropometric measures, and metabolic diseases. Previous studies, including those by Yang et al. (2021) [2] have leveraged TWB data to identify associations among environmental factors, genetic predispositions, lifestyle characteristics, and liver disease risk.

Building on this growing body of evidence, the present study aimed to investigate the synergistic effect of ambient PM2.5 exposure and WHR on NAFLD prevalence in a Taiwanese population. By integrating air quality data with biometric and health information from the TWB, we sought to provide a comprehensive assessment of how environmental and metabolic factors jointly influence NAFLD risk. Understanding these interactions is particularly important in densely populated regions such as Taiwan, where industrial pollution remains a critical public health concern. Ultimately, our findings may inform targeted prevention strategies and public health policies addressing both environmental and metabolic determinants of NAFLD. Accordingly, we analyzed TWB data collected between 2016 and 2020 from 15,049 participants aged 30-70 years to examine the relationship between ambient PM2.5 exposure, WHR, and NAFLD prevalence.

Materials and Methods

Study design and participants

This study employed a cross-sectional design using baseline data from the TWB collected between 2016 and 2020. Participants were community-dwelling adults aged 30-70 years who completed baseline health examinations, including waist and hip circumference measurements, and had available PM2.5 exposure estimates. Individuals assessed in different enrollment years represent independent participants, and no longitudinal follow-up or repeated measurements were conducted. Participants with incomplete data were excluded, resulting in a final analytical sample of 15,049 subjects.

The study protocol was approved by the institutional review committee of Chung Shan Medical University (CS1-22209, CS1-20009, and CS1-23101). All TWB participants provided written informed consent at enrollment.

Exposure assessment and outcome measures

The TWB does not directly collect air pollution data. Air pollution metrics, including PM2.5 concentrations, were obtained from the Taiwan Environmental Protection Administration (EPA), which operates approximately 71 automated monitoring stations across Taiwan and records daily average concentrations of various pollutants. For this study, PM2.5 concentrations from 2000 to 2016 were used to characterize long-term ambient exposure based on participants' residential locations, reflecting cumulative environmental exposure before health assessment.

Annual average PM2.5 concentrations, expressed in micrograms per cubic meter (μg/m³), were derived from daily measurements. To account for spatial and temporal variability across Taiwan, PM2.5 exposure levels were estimated for 349 geographic areas using a machine learning-coupled land-use regression (LUR) model, as described previously [12]. Participants were assigned PM2.5 exposure estimates based on their residential locations.

We acknowledge that the temporal gap between the end of the exposure period (2016) and health assessments for some participants (up to 2020) may raise concerns. Accordingly, the PM2.5 exposure metric was defined to reflect historical long-term exposure rather than contemporaneous PM2.5 levels, consistent with prior epidemiological studies examining the chronic health effects of air pollution. This approach assumes that cumulative exposure to ambient PM2.5 contributes to long-term metabolic disturbances, including the development of NAFLD.

PM2.5 exposure levels were categorized into quartiles as follows: PM2.5 ≤ Q1 (PM2.5 ≤ 28.90 μg/m3), Q1 < PM2.5 ≤ Q2 (28.90 < PM2.5 ≤ 35.99 μg/m3), Q2 < PM2.5 ≤ Q3 (35.99 < PM2.5 ≤ 39.85 μg/m3), and PM2.5 > Q3 (PM2.5 > 39.85 μg/m3).

NAFLD status was determined using liver enzyme measurements and imaging data (e.g., ultrasound) available in the TWB. Potential NAFLD cases were identified based on abnormal liver enzyme values or imaging evidence of hepatic steatosis. Self-reported disease information in the biobank reflects physician-diagnosed conditions.

Waist circumference was measured at the narrowest point between the ribs and iliac crest, and hip circumference at the widest point over the buttocks. WHR was calculated by dividing waist circumference by hip circumference. WHR was classified as normal or abnormal according to sex-specific cutoffs: ≤0.90 for men and ≤0.85 for women were considered normal, while values above these thresholds were considered abnormal.

Covariates

Covariate data were collected using standardized questionnaires and clinical measurements and included age, sex, smoking status, exercise habits, body mass index (BMI), vegetarian diet, coffee consumption, alanine aminotransferase/aspartate aminotransferase (ALT/AST) ratio, total cholesterol, triglycerides, high-density lipoprotein cholesterol (HDL-C), and low-density lipoprotein cholesterol (LDL-C). These variables were included as potential confounders based on prior literature and biological plausibility.

Statistical analysis

All statistical analyses were conducted using SAS software version 9.4. Participant characteristics were summarized using chi-square tests for categorical variables and Student's t-tests for continuous variables. Logistic regression models were used to estimate crude and multivariable-adjusted ORs and 95% CIs for the associations between PM2.5 exposure, WHR, and NAFLD.

Crude models estimated unadjusted associations by including each exposure variable of interest (e.g., PM2.5 quartiles or WHR) separately without additional covariates. Multivariable-adjusted logistic regression models were subsequently constructed to control for potential confounding factors selected a priori based on previous literature and biological plausibility. These covariates (detailed above) are predefined risk factors for NAFLD and may also be associated with exposure to air pollution or body fat distribution.

BMI and WHR were included simultaneously to distinguish overall adiposity from central fat distribution. PM2.5 exposure was categorized into quartiles, with the lowest quartile serving as the reference group. Effect modification by WHR was evaluated by including a multiplicative interaction term between PM2.5 exposure and WHR. Stratified analyses by PM2.5 quartiles were also performed to examine the association between WHR and NAFLD across different exposure levels.

Results

Demographic and baseline characteristics

The study included 15,049 participants, divided into two groups based on the presence of NAFLD: 7870 participants without NAFLD and 7179 with NAFLD (Table 1). Significant differences were observed between the two groups across various demographic and physiological measures. Participants with NAFLD exhibited a higher body mass index (BMI) (26.089 ± 0.044 vs. 22.844 ± 0.034, p < 0.0001), a higher WHR (70.71% with an abnormal ratio vs. 40.00% in the non-NAFLD group, p < 0.0001), and altered lipid profiles, including higher triglycerides (152.600 ± 1.397 vs. 93.982 ± 0.630, p < 0.0001) and lower high-density lipoprotein (49.634 ± 0.136 vs. 59.030 ± 0.154, p < 0.0001). The prevalence of NAFLD was also higher among older adults (73.38% aged ≥ 50 vs. 66.25% in the non-NAFLD group, p < 0.0001) and men (36.51% vs. 25.24%, p < 0.0001).

Association of PM2.5 exposure and waist-hip ratio with NAFLD

In univariate analyses (Table 2), higher PM2.5 exposure was associated with an increased likelihood of NAFLD. Participants in the highest PM2.5 quartile (PM2.5 > Q3) had a significantly higher crude odds of NAFLD compared with those in the lowest quartile (OR = 1.251, 95% CI: 1.114-1.404, p < 0.0001), with a significant trend across quartiles (p-trend < 0.0001). Similarly, an abnormal WHR was associated with higher odds of NAFLD (OR = 1.641, 95% CI: 1.507-1.787, p < 0.0001) compared with a normal ratio. Crude analyses also indicated that men had lower odds of NAFLD than women (OR = 0.848, 95% CI: 0.765-0.940, p = 0.0017), and participants aged ≥50 years had higher odds compared with those <50 years (OR = 1.438, 95% CI: 1.306-1.584, p < 0.0001). Lifestyle factors, including smoking and exercise, were not significantly associated with NAFLD in these univariate analyses.

 Table 1 

Demographic, lifestyle, and clinical characteristics of participants with and without non-alcoholic fatty liver disease

VariablesNo NAFLDNAFLDp-value
(n = 7870)(n = 7179)
Mean PM2.5 (μg/m3), %0.3603
PM2.5 ≤ Q11714 (21.78)1591 (22.16)
Q1 < PM2.5 ≤ Q22085 (26.49)1871 (26.06)
Q2 < PM2.5 ≤ Q32184 (27.75)1925 (26.81)
PM2.5 > Q31887 (23.98)1792 (24.96)
Waist-hip ratio, n (%)<0.0001
Normal4722 (60.00)2103 (29.29)
Abnormal3148 (40.00)5076 (70.71)
Gender, n (%)<0.0001
Women5884 (74.76)4558 (63.49)
Men1986 (25.24)2621 (36.51)
Age, n (%)<0.0001
Age < 502656 (33.75)1911 (26.62)
Age ≥ 505214 (66.25)5268 (73.38)
Smoking status, n (%)<0.0001
Non-smoker6998 (88.92)5962 (83.05)
Current smoker872 (11.08)1217 (16.95)
Exercise, n (%)<0.0001
No4044 (51.39)3936 (54.83)
Yes3826 (48.61)3243 (45.17)
BMI, (kg/m2)22.844 ± 0.03426.089 ± 0.044<0.0001
Vegetarian, n (%)0.3179
No7137 (90.69)6544 (91.15)
Yes733 (9.31)635 (8.85)
Coffee consumption, n (%)0.0727
No4603 (58.49)4095 (57.04)
Yes3267 (41.51)3084 (42.96)
ALT/AST ratio0.763 ± 0.0031.022 ± 0.004<0.0001
Total Cholesterol196.200 ± 0.397199.900 ± 0.442<0.0001
Triglyceride93.982 ± 0.630152.600 ± 1.397<0.0001
High-density lipoprotein59.030 ± 0.15449.634 ± 0.136<0.0001
Low-density lipoprotein118.600 ± 0.347125.300 ± 0.394<0.0001

Abbreviations: NAFLD, nonalcoholic fatty liver disease; PM, particulate matter; n, sample size; %, percent; BMI, body mass index; kg, kilogram; m2, square meters; ALT, alanine aminotransferase; AST, aspartate aminotransferase.

A statistically significant interaction between PM2.5 exposure and WHR on the odds of NAFLD was observed (p for interaction = 0.0388). In stratified multivariable-adjusted analyses by PM2.5 quartiles (Table 3), an abnormal WHR was consistently associated with higher odds of NAFLD across all levels of PM2.5 exposure. The adjusted odds ratios for abnormal versus normal WHR were 1.737 (95% CI: 1.448-2.084) for PM2.5 ≤ Q1, 1.636 (95% CI: 1.383-1.935) for Q1 < PM2.5 ≤ Q2, 1.730 (95% CI: 1.465-2.044) for Q2 < PM2.5 ≤ Q3, and 1.508 (95% CI: 1.275-1.783) for PM2.5 > Q3.

Joint exposure analyses further supported this interaction (Table 4). After adjustment for potential confounders, participants with abnormal WHR and PM2.5 > Q3 had the highest odds of NAFLD (OR = 2.020, 95% CI: 1.708-2.389) compared with the reference group (normal WHR and PM2.5 ≤ Q1). In contrast, among participants with normal WHR, higher PM2.5 exposure was associated with a more modest increase in NAFLD odds (OR = 1.442, 95% CI: 1.215-1.712).

 Table 2 

Univariate (crude) logistic regression analyses examining the associations between PM2.5 exposure, WHR, and other individual factors with odds of NAFLD.

VariablesOR95% CIp-value
Mean PM2.5
PM2.5 ≤ Q1 (ref)---
Q1 < PM2.5 ≤ Q20.9640.860-1.0800.5223
Q2 < PM2.5 ≤ Q31.0340.923-1.1580.5625
PM2.5 > Q31.2511.114-1.4040.0001
p-trendP-value <0.0001
Waist-hip ratio
Normal (ref)---
Abnormal1.6411.507-1.787<0.0001
Sex
Women (ref)---
Men0.8480.765-0.9400.0017
Age
Age < 50 (ref)---
Age ≥ 501.4381.306-1.584<0.0001
Smoking status
Non-smoker (ref)---
Current smoker0.9170.806-1.0440.1913
Exercise
No (ref)---
Yes1.0010.922-1.0870.9804
BMI1.1891.172-1.205<0.0001
Vegetarian diet
No (ref)---
Yes0.9410.820-1.0800.3860
Coffee consumption
No (ref)---
Yes1.0690.897-1.1580.0997
ALT/AST ratio6.8665.898-7.994<0.0001
Total Cholesterol0.9800.975-0.986<0.0001
Triglyceride1.0101.009-1.011<0.0001
High-density lipoprotein1.0010.994-1.0080.7423
Low-density lipoprotein1.0231.017-1.028<0.0001
interaction: PM2.5*waist-hip ratiop-value = 0.0388

OR, odds ratio; CI, confidence interval; ref, reference; BMI, body mass index; ALT, alanine aminotransferase; AST, aspartate aminotransferase; WHR, waist-hip ratio.

Discussion

Our study found a significant interaction between ambient PM2.5 levels and WHR in the risk of developing NAFLD among adults in Taiwan. Individuals with an abnormal WHR coupled with higher quartiles of PM2.5 exposure displayed an increased propensity for NAFLD. Notably, the highest quartile of PM2.5 exposure was associated with a marked increase in NAFLD risk, particularly among those with an abnormal WHR. These findings suggest that increased central adiposity, as indicated by WHR, may exacerbate the hepatotoxic effects of air pollution, aligning with broader evidence on the health impacts of particulate matter [13].

 Table 3 

Associations between WHR and odds of NAFLD stratified by PM2.5 quartiles, estimated using multivariable-adjusted logistic regression models.

VariablesPM2.5 ≤ Q1Q1 < PM2.5 ≤ Q2Q2 < PM2.5 ≤ Q3PM2.5 > Q3
OR95% CIp-valueOR95% CIp-valueOR95% CIp-valueOR95% CIp-value
Waist-hip ratio
Normal (ref)------------
Abnormal1.7371.448-2.084<0.00011.6361.383-1.935<0.00011.7301.465-2.044<0.00011.5081.275-1.783<0.0001
Sex
Women (ref)------------
Men0.9280.744-1.1570.50650.8030.656-0.9840.03410.9060.744-1.1040.32750.7780.633-0.9570.0176
Age
Age < 50 (ref)------------
Age ≥ 501.4851.190-1.8540.00051.4531.206-1.750<0.00011.4411.199-1.732<0.00011.3761.136-1.6670.0011

OR, odds ratio; CI, confidence interval; ref, reference; BMI, body mass index; ALT, alanine aminotransferase; AST, aspartate aminotransferase; WHR, waist-hip ratio.

Adjusted for smoking status, exercise, BMI, diet, coffee intake, ALT, Total cholesterol, HDL-C, LDL-C, and Triglycerides.

 Table 4 

Joint associations of WHR and PM2.5 exposure quartiles with odds of NAFLD, estimated using multivariable-adjusted logistic regression models

VariablesOR95% CIP-value
Waist-hip ratio and mean PM2.5
Normal waist-hip ratio, PM2.5 ≤ Q1 (ref)---
Normal waist-hip ratio, Q1 < PM2.5 ≤ Q21.0100.848-1.2040.9098
Normal waist-hip ratio, Q2 < PM2.5 ≤ Q31.0640.894-1.2650.4843
Normal waist-hip ratio, PM2.5 > Q31.4421.215-1.712<0.0001
Abnormal waist-hip ratio, PM2.5 ≤ Q11.8191.531-2.160<0.0001
Abnormal waist-hip ratio, Q1 < PM2.5 ≤ Q21.6851.432-1.984<0.0001
Abnormal waist-hip ratio, Q2 < PM2.5 ≤ Q31.8331.558-2.157<0.0001
Abnormal waist-hip ratio, PM2.5 > Q32.0201.708-2.389<0.0001

Abbreviation: OR, odds ratio; CI, confidence interval; ref, reference; BMI, body mass index; ALT, alanine aminotransferase; AST, aspartate aminotransferase; WHR, waist-hip ratio.

Adjusted for age, sex, smoking, exercise, BMI, diet, coffee intake, ALT, Total cholesterol, HDL-C, LDL-C, and Triglycerides.

The higher risk for NAFLD among individuals with abnormal WHR and higher levels of PM2.5 may be attributed to several potential mechanisms. Firstly, abnormal WHR, indicative of central obesity, is often associated with adipose tissue dysfunction, which can disrupt lipid metabolism and result in increased release of free fatty acids (FFAs) into the bloodstream. These FFAs can be taken up by the liver, promoting hepatic fat accumulation. This aligns with previous studies highlighting the role of adipose tissue dysfunction in NAFLD development [14, 15]. Secondly, exposure to higher levels of PM2.5 has been linked to systemic inflammation and oxidative stress [16]. Inflammatory mediators and oxidative stress can induce liver injury and inflammation, contributing to the development and progression of NAFLD [17, 18]. Furthermore, abnormal WHR is often associated with insulin resistance, a key factor in the pathogenesis of NAFLD. Exposure to PM2.5 has been shown to exacerbate insulin resistance and impair glucose homeostasis, further contributing to liver fat accumulation and NAFLD development [19, 20].

Regarding the specific association between PM2.5, WHR, and NAFLD, there is limited prior evidence focusing on WHR as a modifier of PM2.5-related NAFLD risk. While several studies have examined the impact of general obesity on the health effects of air pollution, few have specifically investigated how central obesity, as measured by WHR, influences these effects. Our study adds a novel dimension to the existing literature by highlighting the synergistic contribution of central adiposity and air pollution to NAFLD risk.

Beyond its hepatic consequences, NAFLD is increasingly recognized as a multisystem disease, with cardiovascular disease representing its most common and clinically significant comorbidity. Central obesity and metabolic dysfunction are key contributors to both NAFLD and cardiovascular risk. Environmental exposures such as ambient PM2.5 may exacerbate these risks by promoting systemic inflammation and atherosclerosis, thereby linking liver and vascular pathology. Indeed, chronic PM2.5 exposure has been strongly associated with cardiovascular disease incidence and mortality, highlighting its systemic impact [21]. Evidence suggests that NAFLD severity correlates with early markers of atherosclerosis, such as carotid intima-media thickness (IMT). For example, a study by Tarantino et al. [22] demonstrated that carotid IMT is predicted by combined eotaxin levels and hepatic steatosis severity in obese patients with NAFLD, suggesting a direct link between liver fat accumulation and vascular changes. Our findings of synergistic effects between PM2.5 and central obesity on NAFLD risk may therefore have broader implications, potentially contributing to the increased cardiovascular risk observed in individuals with NAFLD, particularly in populations exposed to high levels of air pollution.

Taken together, our findings suggest that the combination of abnormal WHR and higher levels of PM2.5 may synergistically contribute to an increased risk for NAFLD. Adipose tissue dysfunction, inflammation, oxidative stress, and insulin resistance may serve as underlying mechanisms linking these factors. These findings have significant clinical and public health implications. Clinically, individuals with abnormal WHR and higher levels of PM2.5 may benefit from specific interventions. Lifestyle modifications, such as regular exercise programs and weight management strategies, could be recommended to mitigate NAFLD risk in these individuals. Additionally, counseling on measures to reduce exposure to air pollution, such as avoiding high-pollution areas and using air purifiers, may also be beneficial.

A multidisciplinary approach is crucial to address the complex relationship between abnormal WHR, PM2.5, and NAFLD. Collaboration between healthcare providers, public health agencies, and environmental organizations is essential to develop comprehensive strategies that address both individual and environmental risk factors, ultimately contributing to the prevention of NAFLD and improvement of population health.

One of the strengths of our study is the use of a large, well-characterized cohort from the TWB, which allowed for a detailed assessment of both PM2.5 exposure and health outcomes. The application of LUR models for estimating PM2.5 exposure also adds precision to our exposure assessment (Huang et al., 2019). However, our study has limitations. The cross-sectional design limits our ability to establish causality between PM2.5 exposure, WHR, and NAFLD. Furthermore, reliance on self-reported data for some variables may introduce reporting bias. Lastly, despite adjustment for several potential confounders, residual confounding by unmeasured factors, such as dietary intake and genetic predispositions, cannot be ruled out. Nevertheless, our findings provide valuable insights into the clinical and public health implications of the association between abnormal WHR, PM2.5, and NAFLD, and lay the groundwork for future investigations in this area.

Conclusion

In summary, this study underscores the significant correlation between higher levels of PM2.5 exposure and the increased risk of NAFLD, with the greatest risk observed in individuals in the highest quartile of exposure. The consistent elevation in NAFLD risk among individuals with an abnormal WHR across all levels of PM2.5 exposure highlights the critical interplay between environmental and metabolic factors. These findings support integrated public health strategies that address both air pollution and metabolic health to mitigate the burden of NAFLD and its associated cardiovascular complications. Future research should continue to explore these interactions to develop more targeted and effective prevention and management strategies.

Abbreviations

OR: odds ratio; CI: confidence interval; ref: reference; BMI: body mass index; ALT: alanine aminotransferase; AST: aspartate aminotransferase; WHR: waist-hip ratio; NAFLD: non-alcoholic fatty liver disease; PM: particulate matter; WHR: waist-hip ratio.

Acknowledgements

The authors thank Chung Shan Medical University Hospital and the Ministry of Science and Technology for supporting this study financially.

Funding

This work was supported by Chung Shan Medical University Hospital (grant number CSH-2025-D-001) and the Ministry of Science and Technology (grant numbers MOST 111-2121-M-040-002, NSTC 112-2121-M-040-002, NSTC 113-2121-M-040-001, NSTC 114-2121-M-040-002).

Ethics approval and consent to participate

The institutional review committee of Chung Shan Medical University approved the current study (CS1-22209, CS1-20009, and CS1-23101). Taiwan Biobank participants had provided written informed consent during enrollment.

Availability of data and materials

Data are available from the authors upon reasonable request and with permission of the Taiwan Biobank.

Competing Interests

The authors have declared that no competing interest exists.

References

1. Zhang C, Ma T, Liu C, Ma D, Wang J, Liu M. et al. PM2. 5 induced liver lipid metabolic disorders in C57BL/6J mice. Frontiers in Endocrinology. 2023;14:1212291

2. Chen J, Wu L, Yang G, Zhang C, Liu X, Sun X. et al. The influence of PM2. 5 exposure on non-alcoholic fatty liver disease. Life Sciences. 2021;270:119135

3. VoPham T, Kim NJ, Berry K, Mendoza JA, Kaufman JD, Ioannou GN. PM2. 5 air pollution exposure and nonalcoholic fatty liver disease in the Nationwide Inpatient Sample. Environmental research. 2022;213:113611

4. Guo B, Guo Y, Nima Q, Feng Y, Wang Z, Lu R. et al. Exposure to air pollution is associated with an increased risk of metabolic dysfunction-associated fatty liver disease. Journal of hepatology. 2022;76:518-25

5. Li FR, Liao J, Zhu B, Li X, Cheng Z, Jin C. et al. Long-term exposure to air pollution and incident non-alcoholic fatty liver disease and cirrhosis: A cohort study. Liver Int. 2023;43:299-307

6. Atri A, Jiwanmall SA, Nandyal MB, Kattula D, Paravathareddy S, Paul TV. et al. The Prevalence and Predictors of Non-alcoholic Fatty Liver Disease in Morbidly Obese Women - A Cross-sectional Study from Southern India. Eur Endocrinol. 2020;16:152-5

7. Andreasson A, Carlsson AC, Önnerhag K, Hagström H. Waist/hip ratio better predicts development of severe liver disease within 20 years than body mass index: a population-based cohort study. Clinical Gastroenterology and Hepatology. 2017;15:1294-301 e2

8. Özhan B, Ersoy B, Özkol M, Kiremitci S, Ergin A. Waist to height ratio: a simple screening tool for nonalcoholic fatty liver disease in obese children. The Turkish journal of pediatrics. 2016;58:518-23

9. Tarantino G, Capone D, Finelli C. Exposure to ambient air particulate matter and non-alcoholic fatty liver disease. World journal of gastroenterology: WJG. 2013;19:3951

10. Bo Y, Lin C, Guo C, Wong M, Huang B, Lau A. et al. Chronic exposure to ambient air pollution and the risk of non-alcoholic fatty liver disease: A cross-sectional study in Taiwan and Hong Kong. Ecotoxicology and Environmental Safety. 2024;275:116245

11. Hsu C-S, Kao J-H. Non-alcoholic fatty liver disease: an emerging liver disease in Taiwan. Journal of the Formosan Medical Association. 2012;111:527-35

12. Wong P-Y, Su H-J, Lee H-Y, Chen Y-C, Hsiao Y-P, Huang J-W. et al. Using land-use machine learning models to estimate daily NO2 concentration variations in Taiwan. Journal of Cleaner Production. 2021;317:128411

13. Brook RD, Rajagopalan S, Pope III CA, Brook JR, Bhatnagar A, Diez-Roux AV. et al. Particulate matter air pollution and cardiovascular disease: an update to the scientific statement from the American Heart Association. Circulation. 2010;121:2331-78

14. Trayhurn P. Adipose tissue in obesity—an inflammatory issue. Endocrinology. 2005;146:1003-5

15. Jensen MD. Role of body fat distribution and the metabolic complications of obesity. The Journal of Clinical Endocrinology & Metabolism. 2008;93:s57-s63

16. Patel V, Sanyal AJ. Nonalcoholic fatty liver disease. Liver Immunology: Principles and Practice: Springer. 2013 p. 345-59

17. Li N, Sioutas C, Cho A, Schmitz D, Misra C, Sempf J. et al. Ultrafine particulate pollutants induce oxidative stress and mitochondrial damage. Environmental health perspectives. 2003;111:455-60

18. Tilg H, Moschen AR. Evolution of inflammation in nonalcoholic fatty liver disease: the multiple parallel hits hypothesis. Hepatology. 2010;52:1836-46

19. Gastaldelli A, Cusi K, Pettiti M, Hardies J, Miyazaki Y, Berria R. et al. Relationship between hepatic/visceral fat and hepatic insulin resistance in nondiabetic and type 2 diabetic subjects. Gastroenterology. 2007;133:496-506

20. Xu X, Liu C, Xu Z, Tzan K, Zhong M, Wang A. et al. Long-term exposure to ambient fine particulate pollution induces insulin resistance and mitochondrial alteration in adipose tissue. Toxicological Sciences. 2011;124:88-98

21. Feng Y, Fang G, Ouyang J, Guo Y, Yang X, Ji W, et al. Risk of PM<sub>2.5</sub> and its constituents for cardiovascular diseases: A cohort study of over 5.8 million participants. iScience. 2026; 29

22. Tarantino G, Costantini S, Finelli C, Capone F, Guerriero E, La Sala N. et al. Carotid intima-media thickness is predicted by combined eotaxin levels and severity of hepatic steatosis at ultrasonography in obese patients with nonalcoholic fatty liver disease. PLoS One. 2014;9:e105610

Author contact

Corresponding address Corresponding authors: Yung-Po Liaw, Department of Public Health, Institute of Public Health, Chung Shan Medical University, No. 110, Sec. 1 Jianguo N. Rd., Taichung City 40201, Taiwan, Email: Liawypedu.tw, Tel.: +886-4-36097501, Orcid: https://orcid.org/0000-0003-2046-4964. Chien-Ning Huang, Chung Shan Medical University Hospital, No. 110, Sec. 1 Jianguo N. Rd., Taichung City 40201, Taiwan. Email: cshy049com, Tel.: +886-4-24730022#11010, Orcid: https://orcid.org/0009-0009-5936-7900.


Citation styles

APA
Chang, S.T., Chou, Y.H., Wu, C.D., Nfor, O.N., Lu, W.Y., Huang, C.N., Liaw, Y.P. (2026). Synergistic Effect of Ambient PM2.5 Exposure and Waist-Hip Ratio on Non-Alcoholic Fatty Liver Disease Risk in a Taiwanese Population. International Journal of Medical Sciences, 23(5), 1761-1767. https://doi.org/10.7150/ijms.126734.

ACS
Chang, S.T.; Chou, Y.H.; Wu, C.D.; Nfor, O.N.; Lu, W.Y.; Huang, C.N.; Liaw, Y.P. Synergistic Effect of Ambient PM2.5 Exposure and Waist-Hip Ratio on Non-Alcoholic Fatty Liver Disease Risk in a Taiwanese Population. Int. J. Med. Sci. 2026, 23 (5), 1761-1767. DOI: 10.7150/ijms.126734.

NLM
Chang ST, Chou YH, Wu CD, Nfor ON, Lu WY, Huang CN, Liaw YP. Synergistic Effect of Ambient PM2.5 Exposure and Waist-Hip Ratio on Non-Alcoholic Fatty Liver Disease Risk in a Taiwanese Population. Int J Med Sci 2026; 23(5):1761-1767. doi:10.7150/ijms.126734. https://www.medsci.org/v23p1761.htm

CSE
Chang ST, Chou YH, Wu CD, Nfor ON, Lu WY, Huang CN, Liaw YP. 2026. Synergistic Effect of Ambient PM2.5 Exposure and Waist-Hip Ratio on Non-Alcoholic Fatty Liver Disease Risk in a Taiwanese Population. Int J Med Sci. 23(5):1761-1767.

This is an open access article distributed under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0/). See https://ivyspring.com/terms for full terms and conditions.
Popup Image