ORIGINAL RESEARCH article

Front. Public Health, 10 July 2025

Sec. Environmental Health and Exposome

Volume 13 - 2025 | https://doi.org/10.3389/fpubh.2025.1602566

Developing an interpretable machine learning predictive model of chronic obstructive pulmonary disease by serum PFAS concentration

  • XS

    Xiaomei Shao 1† *

  • LZ

    Ling Zhang 2

  • YW

    Yuting Wang 1

  • YY

    Youmei Ying 1

  • XC

    Xueqin Chen 3*

  • 1. Nanjing Jiangbei Hospital, Affiliated Nanjing Jiangbei Hospital of Xinglin College, Nantong University, Jiangsu, China

  • 2. Huai'an No. 3 People's Hospital, Huaian Second Clinical College of Xuzhou Medical University, Jiangsu, China

  • 3. The Affiliated Taizhou People's Hospital of Nanjing Medical University, Taizhou School of Clinical Medicine, Nanjing Medical University, Taizhou, Jiangsu, China

Abstract

Background:

Chronic obstructive pulmonary disease (COPD) is a leading cause of morbidity and mortality worldwide, with limited early detection strategies. While previous studies have examined the relationship between per- and polyfluoroalkyl substances (PFAS) and COPD, limited research has applied interpretable machine learning (ML) techniques to this association.

Methods:

We investigated the association between PFAS exposure and COPD risk in 4,450 National Health and Nutrition Examination Survey (NHANES) participants from 2013 to 2018. After excluding missing covariates and extreme PFAS values and applying K-nearest neighbors (KNN) imputation, nine ML models, including CatBoost, were built and evaluated using metrics like accuracy, area under the curve (AUC), sensitivity, and specificity. The best-performing model was further analyzed using partial dependence plots (PDP) and SHapley additive exPlanations (SHAP) analysis. To enhance clinical applicability, the final model was deployed as a publicly accessible web-based risk calculator.

Results:

CatBoost emerged as the best model, achieving an accuracy of 84%, AUC of 0.89, sensitivity of 81%, and specificity of 84%. PDP revealed that higher perfluorooctane sulfonic acid (PFOS) and perfluoroundecanoic acid (PFUA) levels were associated with reduced COPD risk, whereas perfluorooctanoic acid (PFOA) and 2-(N-Methyl-perfluorooctane sulfonamido) acetic acid (MPAH) showed positive associations with COPD. perfluorononanoic acid (PFNA), perfluorodecanoic acid (PFDE), and perfluorohexane sulfonic acid (PFHxS) demonstrated mixed or non-linear effects. SHAP analysis provided insights into individual predictions and overall variable contributions, clarifying the complex PFAS-COPD relationship. The deployed web-based calculator enables interactive prediction and risk interpretation, supporting potential public health applications.

Conclusion:

CatBoost identified PFOS and PFUA as protective factors against COPD, while PFOA and MPAH increased risk of COPD. These findings emphasize the need for stricter PFAS regulation and highlight the potential of machine learning in guiding prevention strategies.

Introduction

Global burden and trends of COPD

Chronic obstructive pulmonary disease (COPD) is a major global health issue, affecting an estimated 328 million people worldwide (). While smoking is the leading cause, other factors such as biomass fuel exposure, occupational hazards, and air pollution also contribute significantly, especially in low- and middle-income countries (, ). Despite its high prevalence, 70%−80% of COPD cases remain undiagnosed due to the challenges in early detection ().

Machine learning in disease prediction

Machine learning (ML) has emerged as a transformative tool for COPD screening and risk assessment by analyzing complex, multi-dimensional healthcare data (). For instance, Lin et al. () developed a machine learning-based decision system using gradient boosting classifiers (CatBoost, LightGBM, and XGBoost), which achieved an area under the curve (AUC) of 99.85% in identifying high-risk COPD groups. Similarly, Wang et al. () created a COPD risk screening model using logistic regression and generalized additive models, with an AUC exceeding 0.8, showing strong predictive performance. Zeng et al. () developed a ML model using data from over 43,000 COPD patients, achieving an AUC of 0.866 for predicting severe exacerbations within 1 year, outperforming previous models. These studies highlight the potential of ML to improve COPD screening, enhance diagnostic accuracy, and support more effective interventions.

Environmental exposures and respiratory health

Global environmental pollution exposure is widespread, with 91% of the world's population living in areas exceeding WHO safety guidelines for pollutants like PM2.5 and ozone (, ). Environmental conditions are linked to 24% of all deaths globally, with air pollution alone causing 400,000 premature deaths annually in Europe and reducing average life expectancy by 1 year (, ). Niu et al. () found that particulate matter exposure increased COPD exacerbation risk, particularly in younger and severe COPD patients. Yan et al. () demonstrated that higher blood cadmium and lead levels were associated with increased COPD risk, while anthocyanidin intake above 11.56 mg/day reduced cadmium-related COPD risk by 27%. Madani et al. () showed that volatile organic compounds from local sources significantly increased respiratory disease-related emergency room visits, with ethylbenzene having the greatest impact on asthma and COPD. Environmental pollutants pose significant respiratory health risks globally, with effects varying by pollutant type and population vulnerability.

PFAS exposure: background and health impacts

Per- and polyfluoroalkyl substances (PFAS) are a group of synthetic chemicals widely used in industrial and consumer products due to their exceptional chemical stability, water resistance, and heat resistance (). However, their persistence in the environment and bioaccumulation in human tissues have raised significant public health concerns (). PFAS exposure has been linked to various adverse health outcomes, including metabolic disorders (), liver damage (), immune dysfunction (), and respiratory diseases, such as asthma () and reduced lung function (, ). Recent studies have also explored the relationship between PFAS and COPD. For instance, Wang et al. analyzed data from the National Health and Nutrition Examination Survey (NHANES) 2007–2018 and found that perfluorooctanoic acid (PFOA) and PFNA exposure significantly increased COPD risk, particularly in males, with a J-shaped dose-response relationship (, ). Their study further identified serum albumin as a mediator in the association between PFOA and COPD, with a mediation proportion of 17.94%, suggesting potential pathways involving oxidative stress and chronic inflammation (). Despite these advancements, research on PFAS and COPD remains limited, and limited studies have applied ML approaches to investigate this relationship or develop predictive models.

Rationale for model interpretability in public health

Despite emerging evidence linking PFAS exposure to COPD, current research remains limited in both scope and methodology (, ). Most existing studies rely on conventional statistical models, which may not fully capture the complex, non-linear relationships between PFAS and COPD risk, nor do they provide individualized risk estimation (, ). Moreover, few have explored the use of machine learning to enhance predictive performance or model interpretability in this domain. To address these gaps, our study aims to systematically evaluate the relationship between PFAS exposure and COPD risk using advanced ML approaches. By leveraging nationally representative data from the 2013–2018 NHANES, we developed interpretable ML models to predict individual COPD risk, focusing on performance metrics such as AUC, sensitivity, and specificity. We further applied SHAP and partial dependence analyses to uncover both global and personalized insights into how specific PFAS contribute to COPD risk. Finally, to support real-world application, we translated our findings into an accessible online risk calculator, facilitating early screening and informing prevention strategies in public health practice.

Method

Study population

The National Health and Nutrition Examination Survey (NHANES) is a program conducted by the CDC to study the health and nutrition of people living in the United States (). For this study, we used data from three NHANES cycles (2013–2018), which included 29,400 participants. After excluding individuals with missing covariates or serum PFAS concentration data, 4,844 participants remained. Missing values, present in < 20% of the data, were addressed using the K-nearest neighbors (KNN) imputation method. To ensure robust results, we further excluded extreme PFAS values below the 1st percentile and above the 99th percentile (), leaving a final sample of 4,450 participants, as shown in Figure 1. All participants provided written informed consent, and the study was approved by the National Center for Health Statistics Research Ethics Review Board.

Figure 1

Serum PFAS

The PFAS analyzed were 2-(N-Methyl-perfluorooctane sulfonamido) acetic acid (MPAH), perfluorodecanoic acid (PFDE), perfluorohexane sulfonic acid (PFHxS), perfluorononanoic acid (PFNA), perfluorooctane sulfonic acid (PFOS), perfluorooctanoic acid (PFOA), and perfluoroundecanoic acid (PFUA). Total concentrations of PFOS and PFOA were calculated by combining their isomers: linear (n-PFOA) and branched (Sb-PFOA) for PFOA, and linear (n-PFOS) and monomethyl branched (Sm-PFOS) for PFOS. Pearson correlation coefficients were used to evaluate relationships among the seven PFAS.

Covariates

This study included age, gender, race, education level, marital status, body mass index (BMI), family income, and smoking status as covariates. Race was divided into five categories: Mexican American, other Hispanic, non-Hispanic White, non-Hispanic Black, and other. Education was grouped into two levels: high school or less, and more than high school. Marital status options were married, widowed, divorced, separated, never married, and living with a partner. Family income was calculated as a ratio of family income to poverty guidelines, with any value above 5 recorded as 5. Smoking status was defined as having smoked at least 100 cigarettes over a lifetime. To assess multicollinearity among covariates, we calculated the variance inflation factor (VIF). Variables with a VIF < 10 were retained for model construction, consistent with prior methodological recommendations for avoiding instability in multivariate models ().

ML model construction and evaluation

The ML models were built using 15 variables, comprising 10 continuous variables (age, family income, BMI, and seven PFAS biomarkers: MPAH, PFDE, PFHxS, PFNA, PFOA, PFOS, and PFUA) and five categorical variables (gender, race, education level, marital status, and smoking status). Continuous variables were standardized using StandardScaler from scikit-learn to ensure zero mean and unit variance. Categorical variables were encoded as integers without additional transformation. The dataset was randomly split into training (80%, n = 3,560) and testing (20%, n = 890) sets using stratified sampling to maintain the proportion of COPD cases in both sets.

Nine machine learning algorithms were implemented using Python 3.9.19 and scikit-learn 1.3.0: random forest (RF), support vector machine (SVM), decision tree (DT), K-nearest neighbors (KNN), multilayer perceptron (MLP), voting classifier (VC), light gradient boosting machine (LightGBM), CatBoost, and Extreme Gradient Boosting (XGBoost). These models were chosen based on their demonstrated performance in prior studies involving clinical or environmental health prediction tasks (, ).

Hyperparameter tuning was performed using grid search, with the optimized parameters provided in Supplementary Table S1. The workflow of the study is shown in Figure 1. Model performance was evaluated using metrics such as the receiver operating characteristic (ROC) curve, area under the curve (AUC), accuracy, sensitivity (recall), specificity, false-positive rate (FPR), false-negative rate (FNR), positive predictive value (PPV), negative predictive value (NPV), and F1 score. These metrics are widely used in medical machine learning studies to assess both discriminatory power and classification balance, especially under imbalanced conditions (, ).

ML model interpretation

To analyze the impact of individual PFAS on COPD risk, partial dependence plots (PDPs) were created using the sklearn.inspection module with a grid resolution of 50 points. These plots demonstrate how a specific feature influences the model's predictions while holding other variables constant. Using the trained CatBoost model, the relationship between selected features and COPD risk was calculated and visualized. The trends were smoothed using B-spline interpolation (scipy.interpolate.splrep with smoothing parameter s = 30) to enhance readability, and individual variability was highlighted through sample-specific curves. Additionally, rug plots were included to show the distribution of feature values, providing a deeper understanding of their range within the dataset.

SHapley Additive exPlanations (SHAP) analysis was applied to understand how individual features influenced the predictions made by the trained CatBoost model (). The SHAP values, calculated using “TreeExplainer,” provided a breakdown of each feature's contribution to the model output. A combined visualization was created, consisting of a dot plot to display the distribution and direction of feature impacts and a bar plot to rank features by their average contribution. This dual representation provided a clear view of the importance and variability of each feature, offering valuable insights into the factors driving COPD risk predictions. All analysis code and data are made publicly available at https://huggingface.co/spaces/MLML202512/COPD/tree/main for reproducibility.

Web-based risk calculator development

To translate the trained machine learning model into a user-friendly application, we developed an interactive web-based COPD risk calculator using the Gradio framework (https://www.gradio.app/). The calculator was built based on the final CatBoost model, which was trained using selected demographic, socioeconomic, lifestyle, and PFAS biomarker variables. Only the numeric features were standardized using StandardScaler, consistent with the model training pipeline, while categorical variables were kept in their original format as encoded integers. The interface allows users to input raw values for 15 features, including five categorical (gender, race, education level, marital status, and smoking) and 10 numeric variables (age, family income, BMI, and seven PFAS biomarkers: MPAH, PFDE, PFHxS, PFNA, PFOA, PFOS, and PFUA). Upon input, the backend applies the same preprocessing pipeline and uses the trained CatBoost model to generate a binary prediction (COPD or Healthy), a probability score, and a qualitative risk level categorized as low, medium, or high.

Statistical analysis

Continuous variables were reported as means with standard deviations (SD), and categorical variables as counts with percentages. T-tests and chi-square tests were used to compare PFAS levels and demographics between COPD and non-COPD groups. Analyses were performed using Python (3.9.19) and R (4.4.0), with p-value < 0.05 considered significant ().

Result

Baseline characteristics

Among 4,450 participants, as shown in Table 1, 180 (4.0%) had COPD. Participants with COPD were older (64.6 ± 11.5 vs. 49.0 ± 17.6 years, p-value < 0.001) and more likely to be non-Hispanic White (61.7 vs. 36.5%, p-value < 0.001) or have a lower education level (60.0 vs. 43.7%, p-value < 0.001). Marital status also differed, with more widowed individuals in the COPD group (15.6% vs. 6.8%, p-value < 0.001). While smoking prevalence was lower in the COPD group (13.9 vs. 59.7%, p-value < 0.001), this may reflect smoking cessation after diagnosis or survivor bias. PFAS analysis showed higher levels of MPAH (p-value < 0.001), lower PFDE (p-value = 0.004), and lower PFUA (p-value = 0.006) in the COPD group, with no significant differences for PFHxS, PFNA, PFOA, or PFOS.

Table 1

CharacteristicsNon-COPD (N = 4,270)COPD (N = 180)Total (N = 4,450)p value
Age, year49.0 (17.6)64.6 (11.5)49.6 (17.6)< 0.001
Gender, n(%)
 Male2,012 (47.1%)96 (53.3%)2,108 (47.4%)0.263
 Female2,258 (52.9%)84 (46.7%)2,342 (52.6%)
Race, n(%)
 Mexican American690 (16.2%)8 (4.4%)698 (15.7%)< 0.001
 Other Hispanic481 (11.3%)15 (8.3%)496 (11.1%)
 Non-Hispanic White1,558 (36.5%)111 (61.7%)1,669 (37.5%)
 Non-Hispanic Black909 (21.3%)26 (14.4%)935 (21.0%)
 Other632 (14.8%)20 (11.1%)652 (14.7%)
Education, n(%)
High school or lower1,864 (43.7%)108 (60.0%)1,972 (44.3%)< 0.001
More than high school2,406 (56.3%)72 (40.0%)2,478 (55.7%)
Marital status, n(%)
 Married2,191 (51.3%)79 (43.9%)2,270 (51.0%)< 0.001
 Windowed291 (6.8%)28 (15.6%)319 (7.2%)
 Divorced448 (10.5%)37 (20.6%)485 (10.9%)
 Separated158 (3.7%)5 (2.8%)163 (3.7%)
 Never married808 (18.9%)27 (15.0%)835 (18.8%)
 Living with partner374 (8.8%)4 (2.2%)378 (8.5%)
Weight, kg82.7 (22.1)86.2 (27.8)82.8 (22.4)0.558
Height, cm167 (10.2)166 (10.4)167 (10.2)0.985
BMI, kg/m229.7 (7.13)31.1 (9.73)29.8 (7.26)0.480
Family income2.52 (1.61)1.80 (1.22)2.49 (1.60)< 0.001
Smoke, n (%)2,550 (59.7%)25 (13.9%)2,575 (57.9%)< 0.001
MPAH, ng/ml0.161 (0.184)0.231 (0.237)0.164 (0.187)< 0.001
PFDE, ng/ml0.236 (0.202)0.194 (0.182)0.234 (0.202)0.004
PFHxS, ng/ml1.63 (1.45)1.72 (1.42)1.63 (1.44)0.403
PFNA, ng/ml0.695 (0.470)0.636 (0.437)0.692 (0.469)0.188
PFOA, ng/ml1.88 (1.16)2.00 (1.18)1.89 (1.16)0.323
PFOS, ng/ml6.63 (5.32)6.75 (5.07)6.64 (5.31)0.782
PFUA, ng/ml0.153 (0.138)0.120 (0.0970)0.152 (0.137)0.006

Demographic and clinical features of the participants.

BMI, body mass index; MPAH, 2-(N-Methyl-perfluorooctane sulfonamido) acetic acid; PFDE, perfluorodecanoic acid; PFHxS, perfluorohexane sulfonic acid; PFNA, perfluorononanoic acid; PFOS, perfluorooctane sulfonic acid; PFOA, perfluorooctanoic acid; PFUA, perfluoroundecanoic acid.

Serum PFAS concentrations showed significant changes from 2013 to 2018 (p-value < 0.001), as shown in Table 2. PFHxS, PFNA, PFOA, and PFOS levels declined over time, with PFOS dropping from 6.91 ng/ml in 2013–2014 to 6.22 ng/ml in 2017–2018, and PFOA from 2.23 to 1.62 ng/ml. MPAH, PFDE, and PFUA levels remained relatively stable. These trends suggest reduced PFAS exposure, likely due to regulatory measures and shifts in industrial practices. The Pearson correlation analysis showed strong relationships between PFUA and PFDE (r = 0.74) and PFOS with PFNA (r = 0.62), while MPAH exhibited weak correlations with other PFAS (Supplementary Figure S1). These results suggest shared sources or pathways for certain PFAS.

Table 2

PFASNHANES cyclesp value
2013–2014 (N = 1,444)2015–2016 (N = 1,506)2017–2018 (N = 1,500)Total (N = 4,450)
MPAH, ng/ml0.169 (0.194)0.153 (0.183)0.171 (0.183)0.164 (0.187)< 0.001
PFDE, ng/ml0.250 (0.210)0.209 (0.199)0.245 (0.194)0.234 (0.202)< 0.001
PFHxS, ng/ml1.82 (1.54)1.62 (1.46)1.46 (1.31)1.63 (1.44)< 0.001
PFNA, ng/ml0.817 (0.494)0.730 (0.473)0.535 (0.390)0.692 (0.469)< 0.001
PFOA, ng/ml2.23 (1.29)1.82 (1.09)1.62 (1.01)1.89 (1.16)< 0.001
PFOS, ng/ml6.91 (5.21)6.79 (5.48)6.22 (5.21)6.64 (5.31)< 0.001
PFUA, ng/ml0.157 (0.152)0.133 (0.121)0.165 (0.135)0.152 (0.137)< 0.001

Serum concentration of PFAS from 2013 to 2018.

MPAH, 2-(N-Methyl-perfluorooctane sulfonamido) acetic acid; PFAS, Per- and polyfluoroalkyl substances; PFDE, perfluorodecanoic acid; PFHxS, perfluorohexane sulfonic acid; PFNA, perfluorononanoic acid; PFOS, perfluorooctane sulfonic acid; PFOA, perfluorooctanoic acid; PFUA, perfluoroundecanoic acid.

ML models construction and evaluation

Nine ML models, including RF, SVM, DT, KNN, MLP, VC, LGB, CB, and XGB, were constructed and evaluated to predict COPD risk. Performance metrics such as AUC, accuracy, sensitivity, and specificity were used to assess the models, as shown in Table 3. Among these, CatBoost emerged as the best-performing model, achieving the highest accuracy (84%), AUC (0.89), sensitivity (81%), and specificity (84%). The ROC curves in Figure 2 further confirmed the robust performance of CatBoost, showing minimal overfitting and consistent AUC values between training and testing datasets. In contrast, other models like KNN exhibited significant overfitting, with a large performance gap between training (AUC = 0.92) and testing (AUC = 0.69). Given its superior performance, CatBoost was selected as the final model for further analysis.

Table 3

MetricsRFSVMDTKNNMLPVCLGBCBXGB
AUC0.870.820.810.690.870.870.890.890.88
Accuracy (%)846062747171738473
Sensitivity/Recall0.690.940.920.610.920.890.890.810.89
Specificity0.850.590.610.750.700.700.730.840.73
FPR0.150.410.390.250.300.300.270.160.27
FNR0.310.060.080.390.080.110.110.190.11
PPV0.160.090.090.090.110.110.120.180.12
NPV0.991.000.990.980.990.990.990.990.99
F1 score0.260.160.170.160.200.200.210.290.21

Discrimination characteristics among nine ML models.

RF, random forest; SVM, support vector machine; DT, decision tree; KNN, K-nearest neighbors; MLP, multilayer perceptron; VC, voting classifier; LGB, light gradient boosting machine; CB, CatBoost; XGB, extreme gradient boosting; AUC, area under the receiver operating characteristic curve; Accuracy (%), Percentage of correctly classified samples; Sensitivity/Recall, true positive rate; Specificity, true negative rate; FPR, false positive rate; FNR, false negative rate; PPV, positive predictive value; NPV, negative predictive value; F1 score, Harmonic mean of precision and recall.

Figure 2

ML models interpretation

To investigate the relationship between specific PFAS exposure and COPD risk, we performed partial dependence analysis in the trained CatBoost model (Figure 3). The results revealed varying, non-linear associations for different PFAS. COPD risk decreased with higher levels of PFOS and PFUA, suggesting a potential protective effect, while PFOA and MPAH showed a positive association, with risk increasing at higher concentrations. PFNA exhibited a U-shaped relationship, indicating increased risk at both low and high levels, while moderate levels were associated with lower risk. PFDE demonstrated a decreasing trend in risk at moderate levels, followed by an increase at higher concentrations. PFHxS showed a fluctuating pattern without a clear monotonic trend. These findings highlighted the complex influence of PFAS on COPD risk, suggesting that different PFAS may affect the disease through distinct mechanisms.

Figure 3

To further interpret the contributions of individual features to COPD risk, SHAP analysis was performed. Figure 4A illustrated a waterfall plot, which highlighted the impact of key features on an individual prediction. Smoking status had the largest positive contribution to COPD risk, followed by PFNA and MPAH. Conversely, family income and PFUA were associated with reduced risk. The plot clearly showed how individual features influenced the model's prediction for a specific instance. Figure 4B presented a summary plot of SHAP values across the entire dataset, ranking features by their overall importance. Age was the most significant contributor to COPD risk, with older age associated with higher risk. Among PFAS, PFUA, PFHxS, and PFOS demonstrated negative contributions, indicating that lower levels of these PFAS were linked to higher COPD risk. Conversely, MPAH and PFOA showed positive contributions, meaning that higher levels were associated with increased risk. PFNA and PFDE exhibited a mixed effect, with both low and high levels contributing differently to the risk. The SHAP summary plot illustrated these trends, with red indicating feature values that increase COPD risk and blue indicating values that decrease COPD risk, providing a clear and detailed understanding of the directionality of each PFAS's impact on COPD risk.

Figure 4

Web-based risk calculator

To enhance accessibility and clinical applicability, we implemented a web-based COPD risk calculator using the Gradio framework. This interactive tool integrates the trained CatBoost model and allows users to input raw demographic, lifestyle, and PFAS biomarker data through a browser interface (Figure 5). The calculator automatically standardizes numeric features in the backend and provides real-time predictions, including binary classification (COPD or Healthy), probability of risk, and a qualitative risk level (low, medium, or high). The web-based calculator serves as a user-friendly prototype for personalized risk assessment and may assist clinicians or public health professionals in early identification and stratification of COPD risk, particularly in PFAS-exposed populations (https://huggingface.co/spaces/MLML202512/COPD).

Figure 5

Discussion

Summary of main findings and model performance

This study is the first to use interpretable ML techniques to investigate the association between PFAS exposure and COPD risk, utilizing data from the US NHANES (2013–2018). Among the nine ML models tested, CatBoost emerged as the best performer, achieving an accuracy of 84%, AUC of 0.89, sensitivity of 81%, and specificity of 84%, making it the optimal choice for predicting COPD risk. To provide deeper insights, feature importance analysis, partial dependence plots and SHAP analysis were conducted to evaluate how individual PFAS and other factors influence COPD risk. These findings underscored the importance of regulating PFAS exposure to mitigate health risks and demonstrated the potential of interpretable ML methods to identify high-risk populations, guiding targeted interventions and improving public health outcomes.

PFAS as key predictors of COPD risk: consistency with prior studies

Previous research highlights that PFOA and PFNA are strongly associated with increased COPD risk, particularly among males, exhibiting a characteristic nonlinear and J-shaped dose-response relationship for PFOA exposure (). Similarly, Pan et al. demonstrated significant associations between serum levels of PFOS and PFOA and increased COPD risk, noting differential impacts based on sex, age, and smoking status, and indicating protective roles of moderate-intensity physical activity in mitigating PFAS-related COPD risk (). Our study aligned with these findings, as the CatBoost model identified PFAS, particularly MPAH and PFOA, as significant predictors of COPD risk. Notably, our study uniquely identified PFOS and PFUA as potentially protective against COPD risk, differing from findings reported by Wang et al. () and Pan et al. (), who found positive associations for PFOS. These discrepancies may result from variations in demographic characteristics, exposure measurement methodologies, or different adjustments for confounding variables across studies. While previous literature suggests that PFAS may influence COPD development through inflammation and oxidative stress pathways (, ), the specific biological roles of individual PFAS compounds like PFOS and PFUA remain complex and heterogeneous. Thus, further longitudinal and mechanistic studies are needed to clarify these differences and establish causality. Moreover, SHAP analysis in our study highlighted the notable contribution of PFAS to COPD risk, alongside demographic and socioeconomic factors. These results reinforced the hypothesis that PFAS may influence COPD development through mechanisms such as inflammation and oxidative stress (, ), further emphasizing the need for stricter PFAS regulation and further exploration of their impact on respiratory health.

Biological mechanisms underlying PFAS–COPD associations

The observed relationships between PFAS levels and COPD risk in our study can be explained by underlying biological mechanisms, including inflammation (), oxidative stress (), and PFAS interactions with albumin and lung tissues (). For PFOS and PFUA, the protective association at higher concentrations may reflect their ability to stabilize pulmonary surfactants and reduce oxidative stress (). Albumin, known to bind PFOS and PFUA, could facilitate their targeted delivery to lung tissues (), while moderate and lower levels might help maintain epithelial integrity () and mitigate inflammation (), key drivers of COPD progression. In contrast, PFOA and MPAH were positively associated with COPD risk at higher concentrations, which aligned with their known pro-inflammatory and oxidative effects (). PFOA has been shown to activate the NLRP3 inflammasome and increase cytokine production, including IL-6 and TNF-α, leading to sustained inflammation in lung tissues (). MPAH may exert similar effects by disrupting epithelial barriers and exacerbating oxidative stress (), contributing to airway damage and disease progression. These findings highlighted the role of chronic inflammation and oxidative damage as central mechanisms linking higher PFOA and MPAH levels to increased COPD risk.

Non-linear effects of PFNA, PFDE, and PFHxS

The U-shaped relationship observed with PFNA and the mixed pattern with PFDE reflected their dual roles in COPD risk. At moderate concentrations, PFNA and PFDE may exhibit stabilizing effects on lung tissues, potentially reducing inflammation and oxidative stress. However, at very low or high concentrations, these PFAS may disrupt immune homeostasis and amplify inflammatory responses, leading to increased COPD risk (). The fluctuating trend for PFHxS likely stems from its complex interplay with inflammatory and antioxidant pathways, which may vary depending on individual susceptibility and exposure levels (). These findings emphasized the nuanced and concentration-dependent effects of PFAS on COPD risk, highlighting the importance of further mechanistic studies to better understand their roles in respiratory health. These findings emphasized the need for further toxicological studies to elucidate the specific mechanisms by which different PFAS contribute to COPD risk. Experimental research is also needed to determine whether certain PFAS exhibit synergistic or antagonistic effects, particularly in cases of mixed exposure. Understanding these interactions will be critical for developing targeted strategies to mitigate the health impacts of PFAS exposure and for informing regulatory policies aimed at reducing risks associated with these persistent environmental pollutants.

Study limitations

This study has several limitations. First, as NHANES used a multi-stage stratified sampling design, the findings may not fully represent the entire U.S. population. Second, while our machine learning models demonstrated strong predictive performance, they lack external validation on independent datasets, which is essential to assess model stability and generalizability. Third, COPD status in NHANES was based on self-reported questionnaire data rather than spirometry or clinical diagnosis, which may lead to recall bias or disease misclassification. Additionally, smoking status was also self-reported and may be subject to underreporting, particularly among certain demographic groups. Fourth, although we adjusted for several known covariates, potential unmeasured confounders such as physical activity, dietary factors, occupational exposures, and access to healthcare services were not available in our dataset. These variables could influence both PFAS exposure and COPD risk and may have biased the observed associations. Fifth, PFAS concentrations were measured at a single time point, which may not accurately reflect long-term or cumulative exposure levels. Given the chronic nature of COPD, longer-term exposure assessments would provide a more accurate understanding of causal relationships. Furthermore, the exclusion of participants with missing data may have introduced sampling bias, and the lack of access to detailed healthcare records—such as medication history, comorbidities, or imaging findings—limited our ability to fully characterize disease severity or differentiate COPD subtypes. Moreover, this study did not formally compare models with and without PFAS variables, which may limit the assessment of their specific contribution to COPD risk prediction. Finally, cultural and regional differences in environmental exposure, healthcare access, and disease awareness may limit the generalizability of these findings to other populations or countries. These limitations underscore the need for further longitudinal studies incorporating detailed clinical records, long-term exposure measurements, and more comprehensive confounding adjustment to validate and expand upon our findings.

Conclusion

This study explored the relationship between PFAS exposure and COPD risk using NHANES (2013–2018) data, applying interpretable machine learning techniques for the first time. Among the nine models, CatBoost performed best, achieving an accuracy of 84%, an AUC of 0.89, a sensitivity of 81%, and a specificity of 84%, making it the optimal model. PDP analysis revealed that higher PFOS and PFUA levels were associated with reduced COPD risk, while higher PFOA and MPAH increased risk. PFNA, PFHxS, and PFDE showed complex, non-linear associations. SHAP analysis provided individual risk predictions and overall variable contributions, while an interactive web-based calculator was deployed for real-time risk assessment. This is the first study to integrate interpretable ML algorithms with large-scale epidemiological data to examine concentration-dependent effects of individual PFAS compounds on COPD risk. By combining advanced modeling with user-friendly tools, our approach bridges data science and clinical application. These results emphasize the need for PFAS regulatory actions and demonstrate how transparent ML can enhance precision risk stratification in chronic respiratory diseases, providing a scalable framework adaptable to other environmental exposures and health outcomes.

Statements

Data availability statement

The original contributions presented in the study are included in the article/Supplementary material, further inquiries can be directed to the corresponding authors.

Ethics statement

The studies involving humans were approved by all procedures and protocols for the National Health and Nutrition Examination Survey (NHANES) have been reviewed and approved by the National Center for Health Statistics (NCHS) Research Ethics Review Board. The studies were conducted in accordance with the local legislation and institutional requirements. The participants provided their written informed consent to participate in this study.

Author contributions

XS: Methodology, Conceptualization, Formal analysis, Project administration, Software, Writing – original draft, Writing – review & editing. LZ: Conceptualization, Data curation, Methodology, Software, Writing – original draft, Writing – review & editing. YW: Supervision, Validation, Visualization, Writing – original draft, Writing – review & editing. YY: Conceptualization, Methodology, Software, Writing – original draft, Writing – review & editing. XC: Supervision, Writing – original draft, Writing – review & editing.

Funding

The author(s) declare that no financial support was received for the research and/or publication of this article.

Acknowledgments

We appreciate the contribution of all staffs and participants in the U.S. National Health and Nutrition Examination Survey (NHANES).

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Generative AI statement

The author(s) declare that no Gen AI was used in the creation of this manuscript.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpubh.2025.1602566/full#supplementary-material

References

Summary

Keywords

chronic obstructive pulmonary disease, machine learning, partial dependence plot, SHapley additive exPlanations, environment pollution

Citation

Shao X, Zhang L, Wang Y, Ying Y and Chen X (2025) Developing an interpretable machine learning predictive model of chronic obstructive pulmonary disease by serum PFAS concentration. Front. Public Health 13:1602566. doi: 10.3389/fpubh.2025.1602566

Received

30 March 2025

Accepted

18 June 2025

Published

10 July 2025

Volume

13 - 2025

Edited by

Ling Zhang, Wuhan University of Science and Technology, China

Reviewed by

Chengyong Jia, Albert Einstein College of Medicine, United States

Najm Alsadat Madani, State University of New York, United States

Updates

Copyright

*Correspondence: Xiaomei Shao Xueqin Chen

†These authors have contributed equally to this work

Disclaimer

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics