Skip to main content
Download PDF

Association between organophosphorus insecticides exposure and osteoarthritis in patients with arteriosclerotic cardiovascular disease

BMC Public Health volume 24, Article number: 1873 (2024) Cite this article

Abstract

Background

Organic phosphorus insecticides (OPPs) are a class of environmental pollutants widely used worldwide with potential human health risks. We aimed to assess the association between exposure to OPPs and osteoarthritis (OA) particularly in participants with atherosclerotic cardiovascular disease (ASCVD).

Methods

Participants’ information was obtained from data in the National Health and Nutrition Examination (NHANES). Weighted logistic regression models were utilized to detect associations between OPPs metabolites and OA. Restricted cubic spline plots (RCS) were drawn to visualize the dose-response relationship between each metabolite and OA prevalence. Weighted quantile sum (WQS) regression and Bayesian kernel-machine regression (BKMR), were applied to investigate the joint effect of mixtures of OPPs on OA.

Results

A total of 6871 samples were included in our study, no significant associations between OPPs exposure and OA incidence were found in whole population. However, in a subset of 475 individuals with ASCVD, significant associations between DMP (odds ratio [OR] as a continuous variable = 1.22, 95% confidence interval [CI]: 1.07,1.28), DEP ((odds ratio [OR] of the highest tertile compared to the lowest = 2.43, 95% confidence interval [CI]: 1.21,4.86), and OA were observed. DMP and DEP showed an increasing dose-response relationship to the prevalence of OA, while DMTP, DETP, DMDTP and DEDTP showed a nonlinear relationship. Multi-contamination modeling revealed a 1.34-fold (95% confidence intervals:0.80, 2.26) higher prevalence of OA in participants with high co-exposure to OPPs compared to those with low co-exposure, with a preponderant weighting (0.87) for the dimethyl dialkyl phosphate metabolites (DMAPs). The BKMR also showed that co-exposure of mixed OPPs was associated with an increased prevalence of OA, with DMP showing a significant dose-response relationship.

Conclusion

High levels of urine dialkyl phosphate metabolites (DAP) of multiple OPPs are associated with an increased prevalence of OA in patients with ASCVD, suggesting the need to prevent exposure to OPPs in ASCVD patients to avoid triggering OA and further avoid the occurrence of cardiovascular events caused by OA.

Peer Review reports

Introduction

Pesticides are widely used in environments such as crop cultivation, horticulture, and household life because of their ability to prevent pathogenic bacterial infections and kill pests [1]. Organophosphorus insecticides (OPPs) composed about 33% of all pesticide applications in the United States due to their low cost and high efficiency [2]. Especially due to the need to control pests that can spread diseases, these pesticides have been applied in the field of public health, increasing exposure in the population [3]. The large amount of OPPs remaining in food, drinking water, soil, and air can enter the human body through oral, skin, inhalation, and eye contact, causing health hazards [4]. Most OPPs are composed of phosphates (or thiophosphates or thiophosphates) and an organic group. In most cases, these phosphates (or thiophosphates or thiophosphates) are replaced by O, O-dialkyl groups, where the alkyl group is usually dimethyl or diethyl. Once inside the body, OPPs can be converted into oxygenated forms by enzymes and then react with available cholinesterase enzymes. The oxon may undergo either enzymatic or spontaneous hydrolysis to produce of dimethyl dialkyl phosphate metabolites (DMAPs) or diethyl dialkyl phosphate metabolites (DEAPs) and organic group moieties. These metabolites and/or their glucuronate or sulfate conjugates are excreted in the urine [5]. Quantification of exposure to OPPs is critical as exposure to OPPs typically occurs through multiple pathways. Measurement of six common DAP metabolites in urine can provide information on exposure to OPPs or exposure to DAP itself, which may be present in the environment as breakdown products of OPPs (environmental DAP), which is the most commonly used quantification method for assessing exposure to OPPs [6]. Mechanistically speaking, OPPs can cause acute and chronic toxicity by inhibiting acetylcholinesterase activity or cholinesterase independent pathways (inducing protein cross-linking into high molecular weight aggregates) [7]. Since the last century, evidence on the neurotoxicity of OPPs has gradually been confirmed and increasingly reported. Growth and neurodevelopmental limitations in children, as well as cognitive impairment and neurodegenerative diseases in adults, such as Parkinson’s disease and Alzheimer’s disease, are all associated with chronic exposure to OPPs [8, 9]. Recent studies have shown that OPPs can also induce diseases beyond the nervous system, including metabolic diseases (obesity and type 2 diabetes) [10], asthma [11], hearing impairment [12], non-Hodgkin’s lymphoma [13], immune response suppression [14] and acute nephrotoxicity [15].

Osteoarthritis (OA), a degenerative disease characterized by loss of cartilage and changes in joint structure, is estimated to afflict 7.6% of the world’s population (595 million people) in 2020, and is expected to torment 1 billion individuals by 2050, posing a significant burden on public health [16]. As the most common type of arthritis, the incidence of OA has nearly doubled in the industrial era and continues to increase. This suggests that in addition to classic risk factors such as obesity and aging, environmental factors play a significant role in the occurrence of osteoarthritis. Epidemiologic studies suggest that polychlorinated biphenyls (PCBs), perfluoroalkyls (PFAS), lead, and other heavy metals are possibly associated with OA [17]. However, there are no large population-based studies on the association of OPPs metabolites with OA, except for a single-center descriptive study that included only 25 patients with OA, which reported that OPPs were detectable in more than 80% of patient serum samples and were the most prominent of all pesticides [18].

Arteriosclerotic cardiovascular disease (ASCVD) most commonly referred to either coronary artery disease, peripheral artery disease, cerebrovascular disease, or aortic atherosclerosis, with plaques forming in any number of vascular beds [19]. ASCVD has been the leading cause of mortality in developed countries since the mid-20th century, imposing a huge burden on modern societies and becoming a growing public health problem [20]. The comorbidities (e.g., chronic kidney disease, hypertension, diabetes mellitus, dyslipidemia, etc.) of ASCVD are also important risk factors for cardiovascular disease, and these disorders create a mutually reinforcing positive feedback loop of endothelial inflammation, injury, and dysfunction, ultimately delivering multiple hits to the patient [21]. Therefore, prevention and treatment of comorbidities in patients with ASCVD is essential to reduce cardiovascular-caused mortality. Patients with cardiovascular disease often have some genetic changes, which are not only crucial to the development of atherosclerosis, but also participate in the detoxification process of OPPs [22,23,24,25,26]. Therefore, cardiovascular disease patients may be more sensitive to exposure to OPPs, and studying the relationship between pesticide exposure and arthritis occurrence in this population can help evaluate comorbidity risk and provide patients with more comprehensive health management. Reducing pesticide exposure may help improve the quality of life of cardiovascular disease patients and reduce pain and discomfort caused by arthritis. However, no current studies have focused on the relationship between OPPs exposure and OA in a population of patients with cardiovascular disease. Our study aimed firstly to address the issue of whether exposure to OPPs is associated with suffering OA, and secondly to investigate whether OPPs are associated with concomitant OA in ASCVD, a group of patients prone to comorbidities, through publicly available large-scale survey data, and thus provide epidemiologic evidence for public health decision-making.

Materials and methods

Data source and study population

NHANES is a continuous, cross-sectional survey conducted by the Centers for Disease Control and Prevention that combines interviews and physical examinations to study the health and nutritional status of US population. Using a complex, multistage, stratified, clustered probability sampling design, NHANES collects sample data every 2 years to represent the entire population of the United States. A random one-third subsample of all survey participants in eight 2-year cycles from 1999 to 2018 (except 2009 to 2010) had urinary OPPs metabolite measurements and were included in this study.

Measurements of urinary organophosphorus insecticides metabolites

Six measurable dialkyl phosphate (DAP) metabolites were considered non-specific biomarkers of exposure to at least 28 OPPs, including dimethylphosphate (DMP); dimethylthiophosphate (DMTP); dimethyldithiophosphate (DMDTP); diethylphosphate (DEP); diethylthiophosphate (DETP); and diethyldithiophosphate (DEDTP). These metabolites were quantitatively measured by freeze-drying and chemical derivatization of urine specimens followed by isotope-dilution gas chromatography-tandem mass spectrometry (GC-MS/MS) analysis. Urine samples were processed and analyzed according to the method of Bravo et al. [6]. Specifically, 4 milliliters of urine were added to an isotopic labeled internal standard mixture and then concentrated to dryness using azeotropic distillation with acetonitrile. The residue is dissolved in acetonitrile, and dialkyl phosphates (DAPs) are derived into their respective chloropropyl esters using 1-chloro-3-iodopropane and potassium carbonate. The solution containing chloropropyl ester was concentrated and analyzed using gas chromatography-positive chemical ionization-tandem mass spectrometry. Quantification of DAP metabolites using isotope-dilution calibration method. The urine creatinine value of each sample was extracted and used to correct for the concentration of DAPs. The details of laboratory measurements can be found in a previous literature.

Assessment of ASCVD

ASCVD assessment was obtained from the Medical Conditions questionnaire, which was self- and proxy-reported by participants on a broad range of health conditions and medical history for both children and adults. The NHANES MCQ section is generally modeled on the “Medical Conditions” questionnaire section of the U.S. National Health Interview Survey.

The following 4 questions answered by participants are the basis for our assessment:

  1. 1.

    “Ever told you had coronary heart disease?”

  2. 2.

    “Ever told you had angina/angina pectoris?”

  3. 3.

    “Ever told you had heart attack?”.

  4. 4.

    “Ever told you had a stroke?”.

Participants who answered “yes” to one of the questions were considered as ASCVD group. The remaining participants were classified as nonregular users who answered “No” to the second question and those who had never ASCVD.”

Assessment of arthritis

The NHANES medical conditions (variable name prefix MCQ) section provides self- and proxy-reported personal interview data on the conditions of arthritis, including Osteoarthritis or degenerative arthritis (OA), Rheumatoid arthritis (RA), Psoriatic arthritis and others.

The following 2 questions were used to obtain the morbidity of arthritis in participants :

  1. 1.

    “Doctor ever said you had arthritis”.

  2. 2.

    “Which type of arthritis was it?”

Participants who answered “yes” to questions 1 were considered as who had Arthritis. The remaining participants were classified as no Arthritis group.” Furthermore, the Arthritis group was classified as OA, RA, and others.

Other covariates

Demographic variables

Individuals’ self-reported age, sex, race, annual household income, family income-to-poverty ratio, education level, and marital status at interview were obtained from the NHANES demographic files.

Body mass index

Body mass index data of study subjects were recorded in NHANES body measurements, which were collected by trained health technicians in the MEC.

High-density lipoprotein (HDL)& triglycerides& total cholesterol

The HDL-C concentration、Triglycerides concentration and Total cholesterol concentration were calculated from the blood lipids measurements in NHANES survey include total cholesterol, high-density lipoprotein cholesterol (HDL-C), low-density lipoproteins cholesterol (LDL-C), and triglycerides. Data on total cholesterol are provided in Cholesterol - Total (TCHOL_J) file, and HDL-C data are provided in Cholesterol - High - Density Lipoprotein (HDL_J).

Tobacco use

Smoking status was constructed from responses to 2 items from the questionnaire section : 1) “Have you smoked at least 100 cigarettes in your entire life?” and 2) “Do you now smoke cigarettes?” Respondents who reported smoking every day or some days and had smoked > 100 cigarettes were categorized as “current smokers”; respondents who reported currently not smoking but having smoked > 100 cigarettes in the past were categorized as former smokers; and respondents who reported having smoked < 100 cigarettes ever were categorized as nonsmokers.

Alcohol use

Alcohol use data were collected during the MEC interview. A drink was defined as a “12 oz. beer, 5 oz. glass of wine, or one and a half ounces of liquor”. Participants were asked if they had consumed at least 12 drinks in their entire life and in any one year. They were also asked, “In the past 12 months,”: (1) “how often did you drink any type of alcoholic beverage?” (drinking frequency – days per week, month, or year); and (2) “on those days that you drank alcoholic beverages, on the average, how many drinks did you have?” .In this paper, the specific classification criteria of alcohol use refer to the methods in previous literature [27, 28].

Estimated glomerular filtration rate

We calculated estimated glomerular filtration (eGFR) rate using the Chronic Kidney Disease Epidemiology Collaboration equation, as reported in a previous study [29].

Diabetes mellitus

A subject meeting any of the following criteria will be diagnosed with diabetes [18]: (1) Doctor told you that you have diabetes; (2) glycohemoglobin HbA1c(%) > 6.5; (3) fasting glucose (mmol/L) ≥ 7.0; (4) random blood glucose (mmol/L) ≥ 11.1; (5) 2-hour Oral Glucose Tolerance Test (OGTT) blood glucose (mmol/L) ≥ 11.1; 6)use of antidiabetic agents.

Hypertension

Hypertension in NHANES was defined as participant self-reported hypertension (answer “yes” to the question “Have you ever been told by a doctor or other health professional that you had hypertension, also called high blood pressure?”) or elevated blood pressure during physical examination (mean systolic blood pressure ≥ 140 mm Hg, or mean diastolic blood pressure ≥ 90 mmHg).

Physical activity

Physical activity was collected from the Physical Activity questionnaire (PAQ) in NHANES. PAQ was asked in the home, by trained interviewers, using the Computer Assisted Personal Interview (CAPI) system. Different types of physical activity have different MET (Metabolic equivalent) values, and NHANES provides recommended MET values for different types of exercise. The PAQ survey including Vigorous work-related activity (MET = 8), Moderate work-related activity (MET = 4), Walking or bicycling for Transportation (MET = 4), Vigorous leisure-time physical activity (MET = 8) and Moderate leisure-time physical activity (MET = 4). PA can be calculated based on the MET value, activity type, weekly frequency, and duration. We calculated the value of PA according to the following formula: PA (MET-h/wk) = MET × weekly frequency × duration of each physical activity [30].

All adjustment variables and their details are summarized presented in Table S1.

Data Analysis

We utilized the R software (version 4.1.2, Bell Laboratories) for all statistical procedures. Due to the complex sampling design of NHANES, we considered strata, primary sampling unit, and sampling weights in our data analysis lined with the National Center for Health Statistics directions to estimate statistics for the entire US population. Creatinine-corrected concentrations of urinary OPPs metabolites were used for further analysis after logarithmic transformation to a normally distributed distribution.

According to the self-reported medical condition, the included subjects were considered to have one or more diseases (ASCVD, OA, hypertension, diabetes). The presence of OA in the subjects and the concentration of OPPs metabolites were used as outcome variables and exposure variables, respectively, in all analysis processes. All subjects were divided into two subgroups (non ASCVD and ASCVD) based on whether they had ASCVD. Statistical differences between cases with OA and non-OA were evaluated for categorical and continuous variables using chi-square and t tests, respectively. Multivariate survey-weighted logistic regressions adjusted for those covariates that were significantly different between the two groups (OA patients and non-OA) were used to assess the association of OPPs exposure with prevalence of OA in whole populations or subsets of included samples. Concentrations of OPPs metabolites were used as continuous and categorical variables (divided into tertiles) to calculate the odds ratio of OA individually.

In the subsets with ASCVD, subgroup analyses were performed to ascertain whether the association varied across gender, age group, and ethnicity. Nonlinear relationships between levels of various OPPs metabolites and overall OA events were tested by drawing restricted cubic spline (RCS) curves with three knots. Furthermore, we used a weighted quantile sum (WQS) regression, which is widely used in the field of environmental epidemiology [31], to determine the health effects of mixed exposures to these six metabolites of OPPs and the relative weights of single chemical, based on methodologies described in previous literature [32]. Another novel statistical method, Bayesian kernel-machine regression (BKMR), was also used to assess the exposure-response relationship for a mixture of chemicals [33, 34], which addressed two issues: (1) the effect of individual OPPs metabolites on OA when the concentration of other OPPs metabolites was fixed at the median; and (2) the overall effect of OPPs metabolites mixture at different percentile on OA, compared to when all OPPs metabolites were kept at the 50th percentile.

Results

Of the 111,797 participants, 64,165 individuals were included in our primary selection after removing participants with missing OA and ASCVD prevalence information. From these, we further screened a subsample of 12,793 participants registered for urine DAP metabolite testing, and after excluding records with missing values for covariates (n = 5651) and missing sample weights (n = 271), 6871 subjects were included in our final analysis, of whom 475 had ASCVD (Figure S1 shows the flowchart of sample selection).

Patients with ASCVD were older, had a higher BMI, were less physically active, had a predominant history of smoking and alcohol consumption, had elevated ALP and GGT values, had elevated triglycerides and decreased total cholesterol and HDL, and were more often associated with hypertension and diabetes compared to the non-ASCVD population (Table S2). In the non-ASCVD patient population as well as in the overall population, multivariate logistic regression results did not reveal a significant association of these six DAP metabolites with OA (Table S3).

These ASCVD patients have an average age of 61.5 years, with 305 males and 170 females, of which over 80% are overweight (30 kg/m2 > BMI ≥ 25 kg/m2) or obese (BMI ≥ 30 kg/m2). Proportionality of coronary heart disease, angina, myocardial infarction, and stroke., are 47.09%, 33.74%, 44.47%, and 36.00%, respectively. Table 1 and Table S4 summarizes the weighted and unweighted baseline characteristics of the included subjects, respectively. The majority of the selected subjects were non-Hispanic white populations (56.63%), male (64.21%), with a high school education or above (74.31%), and had a history of smoking and alcohol use. Among these ASCVD patients, 76.42% were accompanied by hypertension and 33.68% were accompanied by diabetes. Among these ASCVD patients, 76.42% were accompanied by hypertension and 33.68% were accompanied by diabetes. The highest median concentrations of these OPPs metabolites in the urine samples of the subjects were DMTP, DMP, and DEP, with lower exposure levels for the remaining three metabolites. Of particular note, we can see from Table 1 that 23.01% of ASCVD patients also have OA. Compared to non OA patients, OA patients tend to be elderly, female, with low bilirubin levels, and low eGFR. There were no significant differences in the concentrations of various OPPs metabolites and other covariates between the two groups.

Table 1 Weighted basic characteristic of the participants with ASCVD
Full size table

Each OPPs metabolite concentration was treated as a continuous and categorical variable, respectively, and its relationship with OA occurrence was presented in Table 2. In the crude model, the incidence of OA increased continuously with higher concentrations of DEP and DMP, both as continuous and tertiary categorical variables. In the fully adjusted model, only the concentration of DMP was independently associated with prevalence of OA. An increase in the natural logarithmic units of DMP was associated with an increased prevalence of OA (odds ratio [OR] = 1.22, 95% confidence interval [CI]: 1.07,1.28). When categorized by tertiles, intermediate to high levels of DMP exposure were associated with an elevated prevalence of OA, although this was only statistically significant in the intermediate group (middle tertile compared with lowest tertile, odds ratio [OR] = 2.94, 95% confidence interval [CI]: 1.42,6.09 ). Model 2, adjusted only for age and sex, led to the same conclusion. Unexpectedly, in either model, the prevalence of OA was lower in the moderately DETP-exposed population compared with the lowly exposed population.

Table 2 Association between OPPs exposure and OA among ASCVD patients
Full size table

The results of the subgroup analyses after stratification by sex, age, and race are presented in Table 3. Compared with the lowest tertile, the odds ratio [OR] of DEP in the highest tertile and DMP in the middle tertile for OA among participants older than 60 years was 2.76 (95% confidence interval [CI]: 1.18,6.47) and 3.62 (95% confidence interval [CI]: 1.49,8.82), respectively. In the female population, the odds ratios [OR] of continuous ln-converted DMP as well as intermediate tertile DMP for OA were 1.36 (95% confidence interval [CI]: 1.01,1.82) and 3.39 (95% confidence interval [CI]: 1.50, 7.62), respectively. Among non-Hispanic white populations, the odds ratio [OR] of moderate compared with lowest levels of DMP for OA was 3.61 (95% confidence interval [CI]: 1.55,8.39).

Table 3 Subgroup analyses on the weighted relationship between DEP and DMP with the prevalence of OA among ASCVD patients
Full size table

The RCS curves (Fig. 1) resolved the linearly increasing relationship between ln-transformed urinary DEP or DMP concentrations and OA, which was broadly consistent with the results of multivariate logistic regression. In addition, the remaining four OPPs metabolites showed a U-shaped nonlinear relationship, except for DMTP, which was not significant.

Fig. 1
figure 1

Restricted cubic spline models representing the associations between ln-transformed urinary OPPs metabolite concentrations and osteoarthritis. The model was adjusted for sex, age, bilirubin and eGFR (completely-adjusted model). A: DEP; B: DMP; C: DMTP; D: DETP; E: DMDTP; F: DEDTP

Full size image

Furthermore, two multi-contaminant statistical strategies (WQS regression and BKMR) were used to assess the effect of the OPPs mixture on OA. In the fully adjusted model, the WQS index of the mixture of OPPs was positively associated with an increased prevalence of OA in the overall participants (OR 1.34; 95% CI:0.80, 2.26). DMDTP (35.92%), DMTP (31.83%), and DMP (19.03%) contributed to the WQS index in the top three, followed by DEDTP (5.57%), DEP (4.08%), and DETP (3.56%) (Fig. 2A). The BKMR model illustrated the effect of individual OPP exposure and mixed exposure on OA. We found that exposure to DMP alone was associated with an increased prevalence of OA. The univariate exposure-response curves for DEP, DMTP, DETP, DMDTP, and DEDTP exposure and OA prevalence were relatively flat (Fig. 2B).The cumulative effect of exposure to OPPs was positive, indicating that mixtures of OPPs were positively associated with OA prevalence (Fig. 2C). Although the overall effect of the metabolism of the six OPPs on OA prevalence was not statistically significant when all OPPs metabolites were above their 50th percentile, there was still a clear trend toward an increase.

Fig. 2
figure 2

(A)The weights of each OPPs metabolite in WQS model regression index for osteoarthritis. (B) Dose-response function (95%CI) between each chemical of the OPPs metabolite and osteoarthritis while fixing other chemicals at median values. All OPPs metabolite concentrations were ln-transformed prior to BKMR analysis. (C) Overall effect of the OPE mixture on osteoarthritis by BKMR analysis. All models were adjusted for sex, age, bilirubin and eGFR (completely-adjusted model)

Full size image

Discussion

Our study is the first large-scale cross-sectional investigation of the relationship between OPPs exposure and OA. Our preliminary analyses showed no significant relationship between OPPs exposure and OA prevalence in the whole population and in a subsample of patients without ASCVD (Table S3), and we therefore refocused more on the effect of OPPs on the incidence of OA in patients with ASCVD. DMP and DEP were significantly related to higher prevalence of OA in multivariate regression models, and a linear dose-response relationship was observed. In subsequent subgroup analyses, older adults, females, and non-Hispanic white populations appeared to be more susceptible to the arthrotoxicity of OPPs. Although the correlation between other metabolites is not significant, the RCS curve reveals an increasing trend in the incidence of OA after the concentration of metabolites exceeds a certain threshold. Synergistic exposure to multiple OPPs increased the prevalence of OA in ASCVD patients in multi-pollutant models (WQS and BKMR). Of these, three dimethyl dialkyl phosphate metabolites (DMAPs) contributed the most to the WQS index, implying that the toxicity of these three DMAPs may be more significant, with the toxicity of DMP appearing to be dose-dependent.

We found that the predominant metabolites in the U.S. ASCVD population were DMP, DMTP, DEP, whose urinary concentrations were an order of magnitude higher than those of other DAPs. Despite numerical differences due to heterogeneity of assay methods and biological samples, several population-wide epidemiologic surveys from North America [5, 35], Asia [36], and Europe [37] showed general agreement with the trend of our results. Ye, M. et al. found that DMP, DMTP, and DEP were detected in more than 65% of urine samples from 3,466 Canadian participants aged 20–79 years, with median creatinine-corrected concentrations of 27.5, 14.5, and 17.9 nmol/g, respectively [35]. The main metabolites of OPPs in the Israeli population are DMP and DMTP, whose exposure levels are dozens of times higher than those of adults in the U.S. and Canada attributable to the widespread agricultural application of organic phosphorus pesticides in Israel and the high consumption of fresh and self-produced fruits and vegetables. Notably, analyte concentrations of total DAP (median, 17.6–24.9 µg/g creatinine) were significantly higher in the occupationally exposed population than in the general population [37, 38]. In addition, the concentrations total DAP in the urine of pregnant women in agricultural community and specific urban residential area were much higher than the exposure level among women of childbearing age in general regions [39,40,41]. The participants in the CHAMACOS birth cohort were pregnant women from the agricultural region of Salinas Valley, California, whose urine levels of DAP metabolites were 2.5 times higher than those of NHANES women, with a median concentration of 112.7 nmol/g creatinine [39]. The Generation R study from Rotterdam, Netherlands reported a higher median total DAP exposure in pregnant women compared to CHAMACOS and NHANES, reaching 316.0 nmol/g creatinine [40]. Further analysis revealed that increased intake of fruits and vegetables in the CHAMCOS and Generation R cohorts were the determining factors for elevated levels of DAP metabolites in pregnant women [42, 43]. Noting that certain populations are more susceptible to the damaging effects of environmental exposure to non-high levels of OPPs due to specific genetic characteristics [44]. [39,40,41]. Therefore the next studies should emphasize the toxicological effects in individuals susceptible to OPPs. Previous studies have demonstrated that mutations in the paraoxonase (PON) gene family polymorphisms act as protective factors against coronary heart disease, with a 7% decrease in prevalence in patients with coronary heart disease compared to control subjects [45]. PONs are not only involved in preventing atherosclerosis prophylaxis, but also have a wide range of physiological hydrolytic activities, including drug metabolism and detoxification of neurotoxic agents, and play a role in regulating sensitivity to specific pesticides or neurotoxic agents [22, 24]. Other cardiovascular disease-related gene signatures, such as ABCB1 and GST, have also been proposed to be associated with susceptibility to pesticide toxicity [23, 25, 26]. Our study found that compared to non ASCVD patients, only ASCVD patients observed a significant positive correlation between OA occurrence and OPPs exposure at almost the same dose of OPPs exposure levels (Table 2, Table S2, Table S3). The genetic changes in ASCVD patients seem to provide a new research perspective to explain the susceptibility differences of OPPs to joint toxicity, but further research is needed to demonstrate the correlation differences between the two populations.

The stratified analysis results showed that the susceptibility of joints to OPPs exposure is gender dependent. There was no significant difference in exposure levels of OPPs between female ASCVD patients and female ASCVD patients (Table S5), but in female ASCVD patients, OPPs exposure was significantly positively correlated with concurrent OA (Table 3). Similar findings have been made in other studies regarding the gender differences in the toxicity of organophosphorus compounds. Shaffo, Frances C et al. reported that a widely used organophosphorus pesticide, chlorpyrifos, induces airway hyperresponsiveness in female rats at lower doses and rapid response times [46]. Another study of organophosphorus pollutants found that female zebrafish showed more pronounced phenotypic alterations such as intestinal microbiota compared to males despite lower concentrations of 2-ethylhexyl diphenyl phosphate (EHDPHP) [47]. This phenotypic gender difference that occurs at extremely low exposure levels implies that different target distributions of the action of OPPs between genders, as well as different mediating mechanisms, may be a topic worthy of deeper investigation in the future.

Osteoarthritis is a whole joint disease characterized by alterations in the composition, structure, and material properties of single or multiple joint tissues, including periarticular bone, articular cartilage, and synovial bursae, caused by factors of cellular or matrix origin [48]. Little research exists on the association of OPPs exposure with OA, except for one study that reported detectable serum concentrations of several pesticides, including OPPs, in a population with osteoarthritis [18]. Previous studies have reported that OPPs trigger immunosuppression/dysregulation through different mechanisms, thereby increasing host susceptibility to disease and a range of stressful conditions [49]. OPPs, such as chlorpyrifos and malathion, induce the release of pro-inflammatory cytokines such as IL-1β, IL-6, TNF-α and INF-γ. These cytokines are responsible for the activation and reassembly of nuclear transcription factors such as NF-κB and are involved in inflammation and apoptosis of damaged cells, ultimately lead to inflammatory changes in target organs related to OPPs (such as lungs, nerves, liver, tongue, etc.) [50,51,52,53]. Considering that pro-inflammatory cytokines, one of the key mediators of the pathophysiological disordered process of OA, especially IL-1β and TNF-α, control the degeneration of the articular cartilage matrix [54], it is reasonable to hypothesize that OPPs-mediated secretion of inflammatory mediators may be crucial for the development of OA. Although there is no direct evidence from previous studies indicating that OPPs have an impact on joint tissue, the mechanisms of OPPs in other diseases can interpret the link between OA and OPPs exposure observed in this study and may be a direction for further investigate.

Notably, one of the identified important risk factors for osteoarthritis is obesity [55, 56]. Epidemiologic investigations have shown that chronic exposure to organophosphorus pesticides is associated with increased waist circumference, body mass index, and prevalence of diabetes [57, 58]. Several ex vivo and in vivo experimental studies have also confirmed that OPPs affect glucose and lipid metabolic homeostasis, mainly in the form of hypertriglycerides, hyperglycemia, and hyperinsulinemia [59]. Whereas abnormalities in glucose and lipid metabolism are widely present in ASCVD patients [60], our study also reported that ASCVD patients had higher body mass index, triglyceride, and glycated hemoglobin levels compared to the non-ASCVD population. Therefore, we hypothesize that ASCVD patients with pre-existing glycolipid metabolism problems, such as obesity, are further affected by exposure to OPPs, increasing the frequency of OA. This metabolic factor may be one of the reasons why ASCVD patients are more susceptible to the arthrotoxicity of OPPs compared to the general population, and more in-depth studies are expected to provide more supporting evidence in the future.

The main strength of this study is the combination of multiple statistical methods to investigate the potential impact of exposure to OPPs (alone and mixed) on OA, which could enhance the robustness of our results. As far as we know, this is the first attempt to conduct a large-scale epidemiologic study to examine the relationship between OPPs exposure and OA (both in the general population and in the ASCVD population). Our findings are instructive for follow-up studies and emphasize the need for in vivo and in vitro experiments to elucidate underlying biological mechanisms. However, several limitations of this study should be noted. First, the nature of the cross-sectional design did not allow us to infer a causal relationship between open exposure and OA. In addition, the dialkyl phosphate (DAP) metabolites of OPPs are nonspecific, and we cannot know exactly which OPPs are the most harmful and need to be restricted for daily life use. Subsequent studies should therefore focus on developing measurements of specific biomarkers. Finally, despite adjusting for previously identified confounders in our statistical model, unmeasured or unknown factors may still play a non-negligible role in the association of OPPs exposure with OA.

Conclusion

The current study suggests a possible link between exposure to OPPs (especially DMP and DEP) and an increased prevalence of OA in ASCVD patients. We also provide suggestive evidence for the cumulative effect of OPPs mixtures on OA, which predicts that ASCVD patients should be particularly careful in preventing the arthrotoxicity of OPPs. Further large-scale cohort studies and researches into the biological mechanisms are needed to confirm and clarify our findings.

Data availability

The datasets generated and analysed during the current study are available online in the NHANES repository [https://wwwn.cdc.gov/nchs/nhanes].

References

  1. Tudi M, Ruan HD, Wang L et al. Agriculture Development, Pesticide Application and its impact on the Environment. Int J Env Res Pub He 2021;18(3).

  2. EPA. Environmental protection agency- pesticides. 2017; https://www.epa.gov/sites/default/files/2017-01/documents/pesticides-industry-sales-usage-2016_0.pdf. Accessed 13th Oct, 2023.

  3. El-Nahhal I, El-Nahhal Y. Pesticide residues in drinking water, their potential risk to human health and removal options. J Environ Manage. 2021;299:113611.

    Article  CAS  PubMed  Google Scholar 

  4. Kim KH, Kabir E, Jahan SA. Exposure to pesticides and the associated human health effects. Sci Total Environ. 2017;575:525–35.

    Article  CAS  PubMed  Google Scholar 

  5. Barr DB, Bravo R, Weerasekera G, et al. Concentrations of dialkyl phosphate metabolites of organophosphorus pesticides in the U.S. population. Environ Health Persp. 2004;112(2):186–200.

    Article  CAS  Google Scholar 

  6. Bravo R, Driskell WJ, Whitehead RD, Needham LL, Barr DB. Quantitation of dialkyl phosphate metabolites of organophosphate pesticides in human urine using GC-MS-MS with isotopic internal standards. J Anal Toxicol. 2002;26(5):245–52.

    Article  CAS  PubMed  Google Scholar 

  7. Lockridge O, Schopfer LM, Review. Organophosphorus toxicants, in addition to inhibiting acetylcholinesterase activity, make covalent adducts on multiple proteins and promote protein crosslinking into high molecular weight aggregates. Chem-Biol Interact 2023;376.

  8. Eskenazi B, Bradman A, Castorina R. Exposures of children to organophosphate pesticides and their potential adverse health effects. Environ Health Persp. 1999;107:409–19.

    Article  CAS  Google Scholar 

  9. Sarailoo M, Afshari S, Asghariazar V, Safarzadeh E, Dadkhah M. Cognitive impairment and Neurodegenerative Diseases Development Associated with Organophosphate pesticides exposure: a review study. Neurotox Res. 2022;40(5):1624–43.

    Article  PubMed  Google Scholar 

  10. Czajka M, Matysiak-Kucharek M, Jodlowska-Jedrych B et al. Organophosphorus pesticides can influence the development of obesity and type 2 diabetes with concomitant metabolic changes. Environ Res 2019;178.

  11. Shaffo FC, Grodzki AC, Fryer AD, Lein PJ. Mechanisms of organophosphorus pesticide toxicity in the context of airway hyperreactivity and asthma. Am J Physiol-Lung C. 2018;315(4):L485–501.

    Article  CAS  Google Scholar 

  12. Zou MS, Huang M, Zhang JY, Chen R. Exploring the effects and mechanisms of organophosphorus pesticide exposure and hearing loss. Front Public Health 2022;10.

  13. Hu LQ, Luo D, Zhou TT, Tao Y, Feng JW, Mei SR. The association between non-hodgkin lymphoma and organophosphate pesticides exposure: a meta-analysis. Environ Pollut. 2017;231:319–28.

    Article  CAS  PubMed  Google Scholar 

  14. Li Q, Kawada T. The mechanism of Organophosphorus Pesticide-Induced inhibition of cytolytic activity of Killer cells. Cell Mol Immunol. 2006;3(3):171–8.

    CAS  PubMed  Google Scholar 

  15. Sobolev VE, Sokolova MO, Jenkins RO, Goncharov NV. Molecular mechanisms of Acute Organophosphate Nephrotoxicity. Int J Mol Sci 2022;23(16).

  16. Collaborators GBDO. Global, regional, and national burden of osteoarthritis, 1990–2020 and projections to 2050: a systematic analysis for the global burden of Disease Study 2021. Lancet Rheumatol. 2023;5(9):e508–22.

    Article  Google Scholar 

  17. Deprouw C, Courties A, Fini JB et al. Pollutants: a candidate as a new risk factor for osteoarthritis-results from a systematic literature review. Rmd Open 2022;8(2).

  18. Li AJ, Banjabi AA, Takazawa M, Kumosani TA, Yousef JM, Kannan K. Serum concentrations of pesticides including organophosphates, pyrethroids and neonicotinoids in a population with osteoarthritis in Saudi Arabia. Sci Total Environ 2020;737.

  19. Barrett TJ. Macrophages in atherosclerosis regression. Arterioscler Thromb Vasc Biol. 2020;40(1):20–33.

    Article  CAS  PubMed  Google Scholar 

  20. Rosamond W, Flegal K, Furie K, et al. Heart disease and stroke statistics–2008 update: a report from the American Heart Association Statistics Committee and Stroke Statistics Subcommittee. Circulation. 2008;117(4):e25–146.

    PubMed  Google Scholar 

  21. Feinstein MJ. Multihit Interactions of Antigens, Immune responses, and comorbidities in Cardiovascular Disease Pathogenesis: methods and potential mechanisms. Arterioscl Throm Vas. 2021;41(1):523–5.

    Article  CAS  Google Scholar 

  22. Harel M, Aharoni A, Gaidukov L, et al. Structure and evolution of the serum paraoxonase family of detoxifying and anti-atherosclerotic enzymes. Nat Struct Mol Biol. 2004;11(5):412–9.

    Article  CAS  PubMed  Google Scholar 

  23. Mostafalou S, Abdollahi M. The susceptibility of humans to neurodegenerative and neurodevelopmental toxicities caused by organophosphorus pesticides. Arch Toxicol. 2023;97(12):3037–60.

    Article  CAS  PubMed  Google Scholar 

  24. Costa LG, Cole TB, Jarvik GP, Furlong CE. Functional genomic of the paraoxonase (PON1) polymorphisms: effects on pesticide sensitivity, cardiovascular disease, and drug metabolism. Annu Rev Med. 2003;54:371–92.

    Article  CAS  PubMed  Google Scholar 

  25. Schumacher T, Benndorf RA. ABC Transport Proteins in Cardiovascular Disease-A brief Summary. Molecules 2017;22(4).

  26. Li Y, Li L, Fan D, Wang Z, Cui Y. Effects of GST null genotypes on individual susceptibility to atherosclerotic cardiovascular diseases: a meta-analysis. Free Radic Res. 2020;54(8–9):567–73.

    Article  CAS  PubMed  Google Scholar 

  27. Rattan P, Penrice DD, Ahn JC, et al. Inverse Association of Telomere length with Liver Disease and Mortality in the US Population. Hepatol Commun. 2022;6(2):399–410.

    Article  CAS  PubMed  Google Scholar 

  28. Hicks CW, Wang D, Matsushita K, Windham BG, Selvin E. Peripheral neuropathy and all-cause and Cardiovascular Mortality in U.S. adults: a prospective cohort study. Ann Intern Med. 2021;174(2):167–74.

    Article  PubMed  Google Scholar 

  29. Levey AS, Stevens LA, Schmid CH, et al. A new equation to estimate glomerular filtration rate. Ann Intern Med. 2009;150(9):604–12.

    Article  PubMed  PubMed Central  Google Scholar 

  30. Tian X, Xue B, Wang B, et al. Physical activity reduces the role of blood cadmium on depression: a cross-sectional analysis with NHANES data. Environ Pollut. 2022;304:119211.

    Article  CAS  PubMed  Google Scholar 

  31. Keil AP, Buckley JP, O’Brien KM, Ferguson KK, Zhao S, White AJ. A quantile-based g-Computation Approach to addressing the effects of exposure mixtures. Environ Health Persp 2020;128(4).

  32. Carrico C, Gennings C, Wheeler DC, Factor-Litvak P. Characterization of Weighted Quantile Sum Regression for highly correlated data in a risk analysis setting. J Agr Biol Envir St. 2015;20(1):100–20.

    Article  Google Scholar 

  33. Bobb JF, Valeri L, Claus Henn B, et al. Bayesian kernel machine regression for estimating the health effects of multi-pollutant mixtures. Biostatistics. 2015;16(3):493–508.

    Article  PubMed  Google Scholar 

  34. Bobb JF, Henn BC, Valeri L, Coull BA. Statistical software for analyzing the health effects of multiple concurrent exposures via bayesian kernel machine regression. Environ Health-Glob 2018;17.

  35. Ye M, Beach J, Martin JW, Senthilselvan A. Urinary dialkyl phosphate concentrations and lung function parameters in adolescents and adults: results from the Canadian Health measures Survey. Environ Health Persp. 2016;124(4):491–7.

    Article  CAS  Google Scholar 

  36. Berman T, Goldsmith R, Göen T, et al. Urinary concentrations of organophosphate pesticide metabolites in adults in Israel: demographic and dietary predictors. Environ Int. 2013;60:183–9.

    Article  CAS  PubMed  Google Scholar 

  37. Koureas M, Tsakalof A, Tzatzarakis M, Vakonaki E, Tsatsakis A, Hadjichristodoulou C. Biomonitoring of organophosphate exposure of pesticide sprayers and comparison of exposure levels with other population groups in Thessaly (Greece). Occup Environ Med. 2014;71(2):126–33.

    Article  PubMed  Google Scholar 

  38. Auid LS, Pandiyan A, Vanka J et al. Residual dialkyl phosphate metabolites concentrations of organophosphate pesticides among Indian farm-workers: implication of exposure and hazard assessment. J Occup Environ Med 2023;September 4(1536–5948 (Electronic)).

  39. Bradman A, Eskenazi B, Barr DB, et al. Organophosphate urinary metabolite levels during pregnancy and after delivery in women living in an agricultural community. Environ Health Persp. 2005;113(12):1802–7.

    Article  CAS  Google Scholar 

  40. Ye XB, Pierik FH, Hauser R, et al. Urinary metabolite concentrations of organophosphorous pesticides, bisphenol A, and phthalates among pregnant women in Rotterdam, the Netherlands: the Generation R study. Environ Res. 2008;108(2):260–7.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  41. Sagiv SK, Harris MH, Gunier RB et al. Prenatal Organophosphate Pesticide exposure and traits related to Autism Spectrum disorders in a Population Living in Proximity to Agriculture. Environ Health Persp 2018;126(4).

  42. van den Dries MA, Pronk A, Guxens M, et al. Determinants of organophosphate pesticide exposure in pregnant women: a population-based cohort study in the Netherlands. Int J Hyg Environ Health. 2018;221(3):489–501.

    Article  PubMed  PubMed Central  Google Scholar 

  43. Kalantzi OI, Castorina R, Gunier RB, et al. Determinants of organophosphorus pesticide urinary metabolite levels in pregnant women from the CHAMACOS cohort. Sci Total Environ. 2023;854:158551.

    Article  CAS  PubMed  Google Scholar 

  44. Stykel MG, Humphries K, Kirby MP, et al. Nitration of microtubules blocks axonal mitochondrial transport in a human pluripotent stem cell model of Parkinson’s disease. FASEB J. 2018;32(10):5350–64.

    Article  CAS  PubMed  Google Scholar 

  45. Oliveira SA, Mansur AP, Ribeiro CC, Ramires JA, Annichino-Bizzacchi JM. PON1 M/L55 mutation protects high-risk patients against coronary artery disease. Int J Cardiol. 2004;94(1):73–7.

    Article  PubMed  Google Scholar 

  46. Shaffo FC, Grodzki AC, Schelegle ES, Lein PJ. The Organophosphorus Pesticide Chlorpyrifos induces sex-specific Airway Hyperreactivity in adult rats. Toxicol Sci. 2018;165(1):244–53.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  47. Yang R, Yang Y, Yang L, Chen H, Zhong W, Zhu L. First insight into the sex-dependent accumulation, tissue distribution and potential toxicities of 2-ethylhexyl diphenyl phosphate and its metabolites in adult zebrafish. J Hazard Mater. 2023;452:131299.

    Article  CAS  PubMed  Google Scholar 

  48. Martel-Pelletier J, Barr AJ, Cicuttini FM et al. Osteoarthritis. Nat Reviews Disease Primers 2016;2(1).

  49. Camacho-Pérez MR, Covantes-Rosales CE, Toledo-Ibarra GA et al. Organophosphorus pesticides as modulating substances of inflammation through the Cholinergic pathway. Int J Mol Sci 2022;23(9).

  50. Lasram MM, Lamine AJ, Dhouib IB, et al. Antioxidant and anti-inflammatory effects of N-acetylcysteine against malathion-induced liver damages and immunotoxicity in rats. Life Sci. 2014;107(1–2):50–8.

    Article  CAS  PubMed  Google Scholar 

  51. El-Sayed NM, Ahmed AAM, Selim MAA. Cytotoxic effect of chlorpyrifos is associated with activation of Nrf-2/HO-1 system and inflammatory response in tongue of male Wistar rats. Environ Sci Pollut R. 2018;25(12):12072–82.

    Article  CAS  Google Scholar 

  52. Tigges J, Worek F, Thiermann H, Wille T. Organophosphorus pesticides exhibit compound specific effects in rat precision-cut lung slices (PCLS): mechanisms involved in airway response, cytotoxicity, inflammatory activation and antioxidative defense. Arch Toxicol. 2022;96(1):321–34.

    Article  CAS  PubMed  Google Scholar 

  53. Banks CN, Lein PJ. A review of experimental evidence linking neurotoxic organophosphorus compounds and inflammation. Neurotoxicology. 2012;33(3):575–84.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  54. Kapoor M, Martel-Pelletier J, Lajeunesse D, Pelletier JP, Fahmi H. Role of proinflammatory cytokines in the pathophysiology of osteoarthritis. Nat Rev Rheumatol. 2011;7(1):33–42.

    Article  CAS  PubMed  Google Scholar 

  55. Wang Y, Simpson JA, Wluka AE, et al. Relationship between body adiposity measures and risk of primary knee and hip replacement for osteoarthritis: a prospective cohort study. Arthritis Res Ther. 2009;11(2):R31.

    Article  PubMed  PubMed Central  Google Scholar 

  56. Yusuf E, Nelissen RG, Ioan-Facsinay A, et al. Association between weight or body mass index and hand osteoarthritis: a systematic review. Ann Rheum Dis. 2010;69(4):761–5.

    Article  PubMed  Google Scholar 

  57. Raafat N, Abass MA, Salem HM. Malathion exposure and insulin resistance among a group of farmers in Al-Sharkia governorate. Clin Biochem. 2012;45(18):1591–5.

    Article  CAS  PubMed  Google Scholar 

  58. Starling AP, Umbach DM, Kamel F, Long S, Sandler DP, Hoppin JA. Pesticide use and incident diabetes among wives of farmers in the Agricultural Health Study. Occup Environ Med. 2014;71(9):629–35.

    Article  PubMed  Google Scholar 

  59. Yang F, Li J, Pang G, Ren F, Fang B. Effects of Diethyl Phosphate, a non-specific metabolite of Organophosphorus pesticides, on serum lipid, hormones, inflammation, and gut microbiota. Molecules 2019;24(10).

  60. Skrha J. Jr. Diabetes, lipids, and CV Risk. Curr Atheroscler Rep. 2021;23(3):8.

    Article  PubMed  Google Scholar 

Download references

Acknowledgements

Not applicable.

Funding

This work was funded by the National Natural Science Foundation of China (82300305), Natural Science Foundation of Jiangsu Province (BK20220707), Youth Innovation Research Fund of Jiangning Hospital Affiliated to Nanjing Medical University (JNYYZXKY202203), Research and Development Fund of Kangda College of Nanjing Medical University (KD2022KYJJZD074), and Science and Technology Development Found of Nanjing Medical University (NMUB20220166).

Author information

Author notes

  1. Shenhao Zhu, Yang Zhou and Menglin Chao contributed equally to this work.

Authors and Affiliations

  1. Department of Urology, The First Affiliated Hospital of Nanjing Medical University, No 300 Guangzhou Road, Nanjing, 210029, China

    Shenhao Zhu

  2. Department of Sports Medicine and Joint Surgery, Nanjing First Hospital, Nanjing Medical University, Nanjing, 210006, China

    Yang Zhou & Hongyao Xu

  3. Department of Endocrinology, The First Affiliated Hospital of Nanjing Medical University, No 300 Guangzhou Road, Nanjing, 210029, China

    Menglin Chao

  4. Department of Cardiology, The Affiliated Jiangning Hospital with Nanjing Medical University, Nanjing, 211100, China

    Yuqing Zhang, Weili Cheng, Lai Zhang & Qin Tao

  5. Department of Pharmacy, The First Affiliated Hospital of Nanjing Medical University, No 300 Guangzhou Road, Nanjing, 210029, China

    Qiang Da

Authors
  1. Shenhao Zhu
    View author publications

    You can also search for this author in PubMed Google Scholar

  2. Yang Zhou
    View author publications

    You can also search for this author in PubMed Google Scholar

  3. Menglin Chao
    View author publications

    You can also search for this author in PubMed Google Scholar

  4. Yuqing Zhang
    View author publications

    You can also search for this author in PubMed Google Scholar

  5. Weili Cheng
    View author publications

    You can also search for this author in PubMed Google Scholar

  6. Hongyao Xu
    View author publications

    You can also search for this author in PubMed Google Scholar

  7. Lai Zhang
    View author publications

    You can also search for this author in PubMed Google Scholar

  8. Qin Tao
    View author publications

    You can also search for this author in PubMed Google Scholar

  9. Qiang Da
    View author publications

    You can also search for this author in PubMed Google Scholar

Contributions

Conceptualization, L.Z., Q.T. and Q.D.; methodology, S.Z.; software, Y.Z. and H.X.; formal analysis, S.Z., H.X., and M.C.; data curation, Y.Q.Z. and L.Z.; writing—original draft preparation, S.Z.; writing—review and editing, S.Z., Y.Z., and W.C.; supervision, L.Z. and Q.D. ;All authors reviewed and approved the final manuscript as submitted.

Corresponding authors

Correspondence to Hongyao Xu, Lai Zhang, Qin Tao or Qiang Da.

Ethics declarations

Ethics approval and consent to participate

Data survey conducted by NHANES has been approved by the NCHS Research Ethics Review Board (ERB). All information from the NHANES program is available and free for public, so an individual investigator was not necessary to obtain approval from the institution’s internal ethics review board.

Consent for publication

Not Applicable.

Competing interests

The authors declare no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Electronic supplementary material

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Zhu, S., Zhou, Y., Chao, M. et al. Association between organophosphorus insecticides exposure and osteoarthritis in patients with arteriosclerotic cardiovascular disease. BMC Public Health 24, 1873 (2024). https://doi.org/10.1186/s12889-024-19414-9

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12889-024-19414-9

Keywords

BMC Public Health

ISSN: 1471-2458

Contact us