Association between pyrethroid exposure and osteoarthritis: a national population-based cross-sectional study in the US

Background With the restriction of organophosphorus and other insecticides, pyrethroids are currently the second most-used group of insecticides worldwide due to their advantages such as effectiveness and low toxicity for mammalian. Animal studies and clinical case reports have documented associations between adverse health outcomesand exposure to pyrethroids. At present, the association between chronic pyrethroid exposure and osteoarthritis (OA) remains elusive. Methods Cross-sectional data from the National Health and Nutrition Examination Survey 1999–2002 and 2007–2014 were used to explore the associations of pyrethroid exposure and OA. Urinary level of 3-phenoxybenzoic acid (3-PBA) in urine samples was used to evaluate the exposure of pyrethroid, and OA was determined on the basis of self-reported physician diagnoses. Multivariable logistic regression models were used to investigate the association between pyrethroid exposure and OA. Results Among the 6528 participants, 650 had OA. The weighted geometric mean of urinary volume-based 3-PBA concentration were 0.45 µg/L. With adjustments for major confounders, compared to participants in the lowest quartile of urinary volume-based 3-PBA, those in the highest quartilehad higher odds of OA (odds ratio, 1.39; 95% confidence interval: 1.01, 1.92). There was no nonlinear relationship between urinary volume-based 3-PBA and OA (P for non-linearity = 0.89). Conclusion High urinary 3-PBA concentration was associated with increased OA odds in the US adults. Pyrethroid exposure in the population should be monitored regularly. Supplementary Information The online version contains supplementary material available at 10.1186/s12889-023-16225-2.


Background
Osteoarthritis (OA), mainly characterized by loss of cartilage and changes of the articular structures, is the most commonly occurring form of arthritis [1,2].OA is becoming a worldwide public health concern.Because of the increasing aging population and the increasing rates of obesity, the age-standardized prevalence of OA in the US was increased by 23.2% from 1990 to 2017 [3].It is estimated that the number of adults with OA in U.S. will increase to 67 million by 2030.Although the complex etiology of OA has not been fully crystallized, the environmental factors are increasingly considered as the causes of the disease [4][5][6].
Because of advantages such as effectiveness and low toxicity for mammalian, pyrethroids are the second most used pesticides in the world and are widely used for pest control in fabric manufacturing, indoor environments, and agricultural environments, accounting for more than 30% of the global pesticide use [7,8].Both dietary and non-dietary factors are important sources for human exposure to pyrethroids [9,10].The latest data from the U.S. Food and Drug Administration's Pesticide Residue Monitoring Program show that pesticide residues were detected in about 47% of domestic food and 49% of imported food samples during examination in 2016 [11].Dietary determinants of pyrethroid urine metabolite concentrations include fruits and vegetables, poultry and seafood [12].In addition, potential non-dietary determinants include pyrethroid residues from residential and agricultural applications, which can adhere to dust and soil.Notably, indoor pyrethroid residue concentrations are usually higher than those detected in soil, water, and sediment [13,14].
Adverse effects of pyrethroids on health outcomes are increasingly documented in recent years [15,16].The hormone-like activity of pyrethroid metabolites has been reported to be greater than that of their parent structures, which implies a potential endocrine disrupting effect [17,18].Moreover, thyroid hormone and endocrine hormones, such as total testosterone (TT), estrogen, sex hormone-binding globulin (SHBG), are involved in the pathogenesis of OA [19,20].Therefore, we hypothesis that pyrethroid exposure is associated with increased OA odds.
To our knowledge, there remains elusive about the association between chronic exposure to pyrethroids and OA.To fill this knowledge gap, we investigated the association of urinary 3-phenoxybenzoic acid (3-PBA), a biomarker for measuring pyrethroids exposure, with OA in US adults: first, we investigated the association of urinary 3-PBA with OA ; second, we examined the nonlinear association of urinary 3-PBA with OA; finally, we explored the association of urinary 3-PBA with OA in subgroups.

Study design and population
The National Health and Nutrition Examination Survey (NHANES) is a cross-sectional survey organized by the National Center for Health Statistics of the Centers for Disease Control and Prevention (CDC).Data of NHANES was obtained using a complex, stratified, clustered multi-staged probability sampling design [21].

Definition of OA
Data of OA status was obtained from questionnaire investigation.Briefly, participants aged 20 years and older were asked "Has a doctor or other health professional ever told you that you have arthritis?",participants who answered "yes" in this question were then asked "Which type of arthritis was it?".Participants who answered "osteoarthritis" were classified as having OA, and others were classified as normal [25].A study showed that agreement between self-reported and clinically confirmed OA has reached 81% [26].As healthcare becomes increasingly accessible and people pay more attention to health, the consistency between self-reported and clinical diagnosis is likely to improve [27,28].

Covariates
Age was regarded as a continuous variable; race/ethnicity was grouped as non-Hispanic Black, Hispanic, non-Hispanic White, and others; education attainment was divided into three categories (less than high school graduate, high school graduate or GED, and some college or above); the family income-to-poverty ratio was defined as the total household income divided by the poverty line and then divided as three level (< 1.30, 1.30-3.49,and ≥ 3.50) [29].
Diet quality was evaluated using the Healthy Eating Index-2015 (HEI-2015), with a total score from 0 to 100, and higher scores indicate higher diet quality.Leisure time physical activity (LTPA) was self-reported and was estimated as minutes of moderate recreational physical activity plus twice the minutes of vigorous recreational physical activity.Based on the 2018 Physical Activity Guidelines for Americans [30], participants with LTPA less than 150 min/wk and with LTPA of 150 min/ wk or more were classified as inactive and active, respectively [31].The body mass index (BMI) was calculated as weight in kilograms divided by height in meters squared.Drinking status was classified as "never", "former" (a history of daily binge drinking), "current heavy" (≥ 3 drinks per day for females, ≥ 4 drinks per day for males, or binge drinking [≥ 4 drinks on same occasion for females, ≥ 5 drinks on same occasion for males] on 5 or more days per month), "current moderate" (≥ 2 drinks per day for females, ≥ 3 drinks per day for males, or binge drinking ≥ 2 days per month), or "current mild" (current drinking but not meeting the standard of moderate and heavy drinking) [32].Smoking status was classified as "never" (smoked less than 100 cigarettes in life), "former" (smoked more than 100 cigarettes in life, but they were not currently smoking cigarettes.), or "current smokers" (smoked more than 100 cigarettes in life and smoke some days or every day) [33].

Statistical analysis
We used survey-specific sample weights to account for unequal selection probabilities, oversampling of certain subgroups, and nonresponse adjustment to ensure nationally representative estimates.According to NHANES Analytic and Reporting Guidelines, we incorporated the sample weights for 1999-2002 and 2007-2014 [35].
Geometric mean and its 95% confidence interval (CI) were shown for volume-based and creatinine-corrected urinary 3-PBA concentrations, respectively.Continuous variables were expressed as weighted mean with CIs, and categorical variables were expressed as weighted frequency distribution.Baseline characteristics were compared using the survey-weighted linear regression for continuous variables and the survey-weighted Chisquare test for categorical variables.
Notably, concentrations of creatinine may vary with age, sex, and race/ethnicity.For a diversified population from NHANES, it may not be optimal to use the creatinine-correction value of urinary 3-PBA concentration for analysis [15].Therefore, we used volume-based 3-PBA concentration (unadjusted for creatinine) in the analysis with urinary creatinine added as a separate independent variable [15,36].In addition, we used creatinine-corrected 3-PBA as the primary exposure in sensitivity analyses to test the robustness of identified associations.Urinary 3-PBA was modeled both as a continuous variable with log-transformation, and as quartiles with the lowest quartile regarded as the reference group.Three multivariate logistic regression models were constructed to estimate the odds ratios (ORs) and 95% CIs for associations of 3-PBA and OA.Model 1 was adjusted for urinary creatinine only; model 2 was additionally adjusted for age, gender, race/ethnicity, family poverty-to-income ratio, and education attainment; model 3 was additionally adjusted for diet quality (HEI-2015 score), LTPA, smoking status, drinking status, BMI, DMs, and hypertension.
Linear trends (dose-response relationships) were estimated by assigning the median 3-PBA value for each 3-PBA quartiles as a continuous variable.Moreover, logistic regression based on restricted cubic splines with three knots was conducted to examine the nonlinear association between 3-PBA and OA.

Basic characteristics of study participants
Among the 6528 participants, 650 had OA.The weighted geometric mean (95% CI) of volume-based and creatinine-corrected urinary 3-PBA concentrations were 0.45 µg/L (95% CI: 0.43, 0.48) and 0.47 µg/g (95% CI: 0.44, 0.50), respectively.Volume-based urinary 3-PBA concentrations were higher in non-Hispanic blacks, current smokers than their counterparts, and creatininecorrected 3-PBA concentrations were higher in non-Hispanic whites, current or former smokers, and those over 50 years than their counterparts (Table S1).Compared with participants without OA, those with OA are prone to female, non-Hispanic white, having higher family income-to-poverty ratio, former smoking, former drinking, having hypertension, and having DM.Notably, the levels both of urinary volume-based 3-PBA and urinary creatinine-corrected 3-PBA in participants with OA were significantly higher than those in participants without OA (Table 1).The details of quartiles of the 3-PBA level are shown in Table 2.

Associations between 3-PBA and OA
The OR of OA was positively correlated with urinary volume-based 3-PBA concentrations in all models (Table 3).Each unit increase in ln (urinary volume-based 3-PBA) was associated with 9% increased OA odds.Furthermore, compared with participants in the lowest quartile, those with the highest quartile of urinary volume-based 3-PBA had higher odds of OA (OR, 1.39; 95% CI: 1.00, 1.92) (P for trend = 0.07).Regressions based on restricted cubic splines showed that there was no nonlinear relationship between urinary volume-based 3-PBA and OA (P for non-linearity = 0.89) (Fig. 1).

Sensitivity analysis
To examine the robustness of association between urinary 3-PBA concentration and OA, we regarded urinary creatinine-corrected 3-PBA as the primary exposure and conducted same analyses as volume-based 3-PBA.As expected, we got consistent results (Table S2, Figure S2, and Table S3).

Discussion
In this nationally representative sample of US adults, we documented that increased levels of urinary 3-PBA concentrations were associated with increased OA odds.
Moreover, these association remained consistent in sensitivity analyses and in most subgroups.
Different availability of pyrethroid-contained products or environments may confer the different pyrethroids exposure for population in different country.In this study, we found that the geometric mean of creatinine-corrected 3-PBA concentration was 0.47 µg/g creatine in US adults.It was lower than reported in China (0.93 µg/g) [37] and South Korea (1.95 µg/g) [38], but was higher than that in Japan (0.40 µg/g) [39].Notably, an increasing trend of urinary 3-PBA levels in US children and adults was documented [8], suggesting that exposure of pyrethroid has been an urgent public health issue.
3-PBA has anti-thyroid hormone activity and exhibits strong thyroid hormone receptor antagonistic capacity [40].Long-term pyrethroid exposure in mice affects levels of thyroxine and triiodothyronine in the serum and brain tissue [41].Furthermore, population study has documented that high urinary 3-PBA levels are associated with the decrease of thyroid hormones level [38].Notably, thyroid hormones play a key role in the development of OA through stimulating terminal chondrocyte differentiation and transglutaminase activity in articular cartilage [42,43].Therefore, pyrethroids exposure with anti-thyroid hormone activity may exert a prominent impact on development of OA.
Pyrethroids can also affect sex hormone.Pyrethroids and their metabolites showed anti-estrogenic activity in human and rat estrogen receptor α [18,44], which may lead to type II collagen degradation and related structural alterations [45].Moreover, high urinary 3-PBA is associated with increased SHBG and TT, and levels of TT and SHBG are positively associated with OA risk [46,47].Therefore, we reasonably speculate that 3-PBA increases the odds of OA through affecting thyroid hormone and sex hormone levels in the body.
At present, the risk factors confirmed for OA including increase of age, genetics, obesity, females, and exercise [5,6].Interestingly, we found that 3-PBA was significantly associated with the odds of OA in these corresponding subgroups, indicating that pyrethroids or their metabolites may interact with these risk factors for OA.As shown in our subgroup results, we observed a significant interaction between sex and 3-PBA for OA, indicating a stronger endocrine disruptor activity of pyrethroids in female than in male.Indeed, female are at higher risk of OA than male, which may be conferred by difference of genetics and anatomic.Heritability of OA is stronger in female than in male [48], moreover, the thickness of articular cartilage of distal femur and patella in female is lower than that in male, and femurs in female are narrower than those in male [48].In addition, changes in hormone levels may mediate age-sex interaction in OA risk.For example, hormone levels in the body changes The bold values mean statistical significance Abbreviations: OA, osteoarthritis; 3-PBA, 3-phenoxybenzoic acid; LTPA, leisure-time physical activity; Q, quartile; BMI, body mass index; DMs, diabetes mellitus status; 3-PBA, 3-phenoxybenzoic acid greatly after menopause, leading to an increased OA risk [47]; causal relationships of TT and SHBG with OA were more likely to be observed in female than in male [46,47].Similarly, association between 3-PBA and TT varies by sex, that a positive association among female but no association among male [49].Moreover, TT is mainly secreted by testes in male, while TT is directly secreted by adrenal gland or ovary, or is indirectly transformed from testosterone precursor in female [49,50].In this case, pyrethroids and their metabolites may target different action pathways for female and male, leading to heterogeneity of association between 3-PBA and OA.However, the potential relationship between pyrethroids and OA remains to be explored.There existed heterogeneity of association among others subgroups, including age, race, BMI, LTPA, and DMs.In addition to the lower power of test caused by small samples, several biological and behavioral factors can partly explain this heterogeneity.We observed a significant positive association between urinary 3-PBA and OA odds in overweight or obese participants rather than in normal weight participants.One reasonable explanation may be that extra mass increases the joint load and causes additional stress on the articular cartilage, thus leading to degenerative changes in the weight bearing joints and OA susceptibility [51].Another explanation may be that effect of endocrine disrupting chemicals may be modified by metabolic changes caused by fat distribution, such as modulations in insulin resistance or leptin production [14].Interestingly, we found that participants with actively LTPA seem to be more susceptible to the effect of pyrethroid exposure.Indeed, the relationship between physical activity and OA is complex.Although the role of physical activity in the treatment of OA has been highlighted [52], high intensity physical activity seems to increase the odds of OA, and it is unclear whether this association is due to physical activity or injury [6].The potential modifying effect of physical exercise on association between pyrethroid exposure and OA needs to be further clarified.
Some strengths in this study must be underscored.We provided novel insights into the association between pyrethroid exposure and health outcomes.Moreover, we considered a wide range of potential confounders including socioeconomic status, dietary, lifestyle, and comorbidities for analysis.Limitations should be acknowledged for interpreting our findings.Firstly, the identified association between pyrethroid exposure and OA may not be causal because of a cross-sectional design; secondly, the data of 3-PBA in the NHANES comes from one single spot urine samples, and the short half-life of pyrethroids in humans may damage the credibility of 3-PBA as a reflection of long-term exposure status [53]; finally, the assessment of OA was collected by questionnaire survey, thereby increasing the risk of measurement bias.

Conclusion
Urinary 3-PBA concentration was positively associated with OA odds in the US adults.Our study highlights for the first time the association between chronic exposure to pyrethroids and OA.Pyrethroid exposure in the population should be monitored regularly.

Table 1
Baseline characteristics of study participants (N = 6528)

Table 2
Details of quartiles of the 3-PBA level

Table 3
OR (95% CIs for association between urinary volume-based 3-PBA and OA Model 2: Adjusted for age (continuous, years), sex (male, female), race/ethnicity (Hispanic, non-Hispanic white, non-Hispanic black, and others), education (less than high school, high school and college or higher), and family poverty-to-income ratio (< 1.3, 1.3-3.5, and > 3.5) plus variables in Model 1 Model 3: Adjusted for diet quality (continuous, HEI-2015 score), LTPA (inactive and active), smoking status (never, former, and current), drinking status (never, former, mild, moderate, and heavy), BMI, DMs (diabetes, prediabetes, and normal) and hypertension (yes and no) plus variables in Model 2 Abbreviation: OR, odds ratio; CIs, confidence intervals; OA, osteoarthritis; 3-PBA, 3-phenoxybenzoic acid; Q, quartile This study utilized publicly available data from the National Health and Nutrition Examination Survey.This data can be accessed using the following link: https://www.cdc.gov/nchs/nhanes/index.htm.
DeclarationsEthics approval and consent to participateThe survey was approved by the National center for health statistics (NCHS) Research Ethics Review Board.All participants provided written informed consent.NCHS Research Ethics Review Board approval was waived for this analysis because of the publicly available and de-identified data.All methods

Table 4
OR (95% CIs) for association between urinary volume-based 3-PBA and OA in subgroups All models were adjusted for urinary creatinine, age, sex, race/ethnicity, education, family poverty-to-income ratio, HEI-2015, LTPA, smoking status, drink status, BMI,DMs and hypertensionThe bold values mean statistical significance Abbreviation: OR, odds ratio; CIs, confidence intervals; OA, osteoarthritis; 3-PBA, 3-phenoxybenzoic acid; Q, quartile; BMI, body mass index; LTPA, leisure-time physical activity; DMs, diabetes mellitus status were carried out in accordance with relevant guidelines and regulations (Declaration of Helsinki).