Perfluoroalkyl and polyfluoroalkyl substances (PFASs) are widely used in consumer products. However, the role of PFAS in infertility is still poorly understood. A total of 788 women from the 2013-2016 nationally representative NHANES were included to explore the association between PFAS exposure and self-reported infertility. Six PFAS, including PFDE, PFNA, PFHxS, n-PFOA, n-PFOS, and Sm-PFOS, were detected by online SPE-HPLC-TIS-MS/MS. We used the generalized linear regression model (GLM), generalized additive models (GAM), and Bayesian kernel machine regression (BKMR) to assess the single effects, non-linear relationships, and mixed effects on women's infertility, respectively. The prevalence of self-reported infertility was 15.54% in this study. In GLM, n-PFOA showed a negative association with self-reported infertility in women for the Q3 (OR: 0.396, 95% CI: 0.119, 0.788) and Q4 (OR: 0.380, 95% CI: 0.172-0.842) compared with Q1 (p for trend = 0.013). A negative trend was also observed in n-PFOS and ∑PFOS (p for trend < 0.05). In GAM, a non-linear relationship was revealed in Sm-PFOS, which exhibits a U-shaped relationship. The BKMR model indicated that there might be a joint effect between PFAS and women's infertility, to which PFNA contributed the highest effect (PIP = 0.435). Moreover, age stratification analysis showed a different dose-response curve in under and above 35 years old. Women under the age of 35 have a more noticeable U-shaped relationship with infertility. Therefore, the relatively low level of mixed PFAS exposure was negatively associated with self-reported infertility in women in general, and the impact of PFAS on infertility may vary among women of different age groups. Further studies are needed to determine the etiological relationship.
Keywords: Bayesian kernel machine regression (BKMR); PFAS; generalized additive models (GAM); generalized linear model (GLM); infertility; mixed effect.