Skip to main content

Biological aging mediates the associations between urinary metals and osteoarthritis among U.S. adults

Abstract

Background

Osteoarthritis (OA) is a worldwide public health concern, mainly afflicting older adults. Although the etiology of OA remains unclear, environmental factors are increasingly considered as non-negligible risk factors. This study aims to evaluate the associations of urinary metals with OA risk and the mediated effect of biological aging.

Methods

Nine urinary metal concentrations were detected among 12,584 U.S. adults based on the National Health and Nutrition Examination Survey (NHANES), including barium (Ba), cadmium (Cd), cobalt (Co), cesium (Cs), molybdenum (Mo), lead (Pb), antimony (Sb), thallium (Tl), and uranium (Tu). Multivariable logistic regression and weighted quantile sum (WQS) regression were used to explore the associations of single metal and mixed metals with OA risk, respectively. Furthermore, biological aging was measured from different perspectives, including cell senescence (telomere length) and whole-body aging (phenotypic age and biological age). Mediation analyses were conducted to investigate the mediated effects of aging on the associations of metals with OA risk.

Results

In the single-exposure model, Cd, Co, and Cs were identified to be positively associated with OA risk, with odds ratios (OR) ranging from 1.48 to 1.64 (all P < 0.05). Mixed-exposure analyses showed consistent associations (OR 1.23, 95%CI 1.10 to 1.37) and highlighted that Cd, Co, and Cs were responsible for the outcomes. Additionally, Cd, Co, Cs, Pb, and Tl were positively associated with biological aging markers, while all biological aging markers had significant associations with OA risk. Further mediation analyses showed that the associations of single metal (mainly Cd and Cs) and mixed metals with OA risk parallelly mediated by the above biological aging markers, with the proportion of mediation ranging from 16.89 to 69.39% (all P < 0.05). Moreover, such associations were also serially mediated through telomere length-biological age path and telomere length-phenotypic age path (the proportion of mediation: 4.17–11.67%), indicating that metals accelerated cell senescence to lead to whole-body aging and finally aggravated OA progress.

Conclusions

These findings suggested that exposure to metals increased OA risk, which was possibly and partly mediated by biological aging.

Peer Review reports

Background

Osteoarthritis (OA) is a joint disorder with progressive cartilage breakdown, subchondral remodeling, and synovial inflammation. With an increasingly aging population, this syndrome is becoming a worldwide public health concern and a leading source of societal cost in older adults [1]. Epidemiologically, the global age-standardized point prevalence of OA was 3.75% in 2017 and the age-standardized prevalence from 1990 to 2017 was increased by 23.2% in the USA [2]. It is estimated that the number of U.S. people afflicted by OA will rise to 67 million by 2030 [3, 4]. Although the etiology of OA has not been fully elucidated, the non-negligible environmental factors are growingly considered to account for this disease.

Heavy metals are ubiquitously present in various environmental media, such as air, soil, drinking water, and food [5, 6]. The widespread exposure to toxic metals leads to illnesses, such as diabetes [7] and cancer [8]. Of concern, heavy metals (e.g., cadmium (Cd) and lead (Pb)) were cumulative bone-toxicants with the contribution to musculoskeletal disease [9]. Compelling toxicological evidence showed that heavy metals were linked to two major pathological changes of OA, including joint damage and inflammation [10, 11]. For example, cadmium chloride treatment caused cartilage loss and inflammation in a 3D culture model of human chondrocytes [12]. Also, a cross-sectional study showed blood Pb was significantly associated with knee OA in Korean postmenopausal women [13]. Previous epidemiological researches have presented positive associations of heavy metals with rheumatoid arthritis [14] and osteoporosis [15]. However, population study on the impact of heavy metals on OA is limited. Moreover, the role of metals in diseases usually relied on their cooperation and interactions. Individual metals cannot completely account for the occurrence and development of disease. Hence, it is worthwhile to explore the joint effects of metals on OA. Additionally, underlying mechanisms of metals in OA remain largely unrevealed.

Aging, a progressive decline in body homeostasis, is considered a target for metal action. Toxicologically, heavy metals can render cell senescence and body premature aging by mitochondrial dysfunction, oxidative stress, and DNA damage [16, 17]. Epidemiological research also reported the positive associations of Cd and Pb with telomere wear [18]. On the other hand, aging is the major pathogenesis for OA. Researchers emphasized that joint tissues are subject to aging and decay over time and the senescence of joint homeostasis-related cells contributed to the development of OA [19]. Cohort studies also showed a negative correlation between telomere length and OA incidence [20]. Importantly, elimination of senescent chondrocytes improved joint degeneration in OA mice, indicating biological aging was responsible for OA development [21]. Collectively, based on the importance of biological aging for metal exposure and OA respectively, we hypothesized that heavy metal exposure may increase OA risk by promoting biological aging.

Hence, we conducted a cross-sectional study to investigate the associations of 9 urinary metals with OA risk based on the National Health and Nutrition Examination Survey (NHANES) 1999–2016. Further, we measured biological aging status at multiple perspectives and explored the mediated effects of biological aging.

Methods

Study population

The NHANES is an ongoing national cross-sectional survey and data is available on the website of the America Centers for Disease Control and Prevention (CDC) (http://www.cdc.gov/nchs/nhanes.htm). The study protocol was approved by the National Center for Health Statistics (NCHS) Research Ethics Review Board. All participants provided written consent at the time of recruitment. In the current study, OA patients were included according to self-reported information (n = 46,736). Next, we excluded individuals whose information on 9 metals was missing (n = 34,147). Further, we excluded individuals with missing creatinine data (n = 5). Collectively, 12,584 participants were enrolled, including 1356 OA patients (Fig. 1).

Fig. 1
figure 1

Flow chart of participants selection. NHANES, National Health and Nutrition Examination Survey

OA assessment

A study showed an 81% agreement between self-report and clinical confirmation of OA [22]. In NHANES, OA was diagnosed by professionals and the information was collected by questionnaire survey. Briefly, all participants ≥ 20 years of age were asked two arthritis-related questions. First, they were inquired: “Has a doctor or other health professional ever told you that you have arthritis?” Individuals with an affirmative response were asked a follow-up question: “Which type of arthritis was it?” The participants who answered “osteoarthritis” were enrolled in the study.

Metal measurement

The detection data of 9 urinary metals were obtained from NHANES 1999–2016. Barium (Ba), Cd, cobalt (Co), cesium (Cs), molybdenum (Mo), Pb, antimony (Sb), thallium (Tl), and uranium (Tu) were primarily measured from spot urine samples by using inductively coupled plasma mass spectrometry (ICP-MS). The limit of detection (LOD) divided by the square root of two was used to replace the values below the LOD. The concentrations of metals were all corrected by urinary creatinine and were expressed as micrograms per gram creatinine.

Measurement of biological aging markers

Biological age and phenotypic age used different calculation methods and incorporated different biomarkers to measure biological aging. Biological age was proposed by Klemera and Doubal [23] based on 8 biomarkers (Ln-C-reactive protein (CRP), serum creatinine, glycosylated hemoglobin, serum albumin, serum total cholesterol, serum urea nitrogen, serum alkaline phosphatase, and systolic blood pressure) [24]. The values j and i represent the number of biomarkers and samples respectively. The values k, q, and s are the regression slope, intercept, and the root means squared error of a biomarker regressed on chronological age, respectively. The value rj2 represents the variance explained by regression of chronological age on biomarkers.

$${\text{BA}}_{E}=\frac{{\sum }_{j=1}^{m}\left({x}_{j}-{q}_{j}\right)\left(\frac{{k}_{j}}{{s}_{j}^{2}}\right)}{{\sum }_{j=1}^{m}{\left(\frac{{k}_{j}}{{s}_{j}}\right)}^{2}}$$
(1)
$${r}_{char}=\frac{{\sum }_{j=1}^{m}\frac{{r}_{j}^{2}}{\sqrt{1-{r}_{j}^{2}}}}{{\sum }_{j=1}^{m}\frac{{r}_{j}}{\sqrt{1-{r}_{j}^{2}}}}$$
(2)
$${s}_{BA}^{2}={\frac{{\sum }_{j=1}^{n}\left(\left({\text{BA}}_{Ei}-{\text{C}}{\text{A}}_{i}\right)- \frac{{\sum }_{i=1}^{n}\left({\text{BA}}_{Ei}-{\text{C}}{\text{A}}_{i}\right)}{n}\right)}{n}}^{2}-\left(\frac{1-{r}_{char}^{2}}{{r}_{char}^{2}}\right)\times \left(\frac{{\left({\text{CA}}_{max}- {CA}_{min}\right)}^{2}}{12m}\right)$$
(3)
$$\mathrm{Biological\ age}=\frac{{\sum }_{j=1}^{m}\left({x}_{j}-{q}_{j}\right)\left(\frac{{k}_{j}}{{s}_{j}^{2}}\right)+ \frac{\text{CA}}{{s}_{BA}^{2}}}{{\sum }_{j=1}^{m}{\left(\frac{{k}_{j}}{{s}_{j}}\right)}^{2}+ \frac{1}{{s}_{\text{BA}}^{2}}}$$
(4)

Phenotypic age was calculated by using 9 aging-related variables [25], including chronological age, albumin, creatinine, glucose, CRP, lymphocyte percent, mean cell volume, red blood cell distribution width, alkaline phosphatase, and white blood cell count.

$$\mathrm{Phenotypic\ age} =141.50 +\frac{\mathrm{Ln}[-0.00553 \times \mathrm{Ln}(\mathrm{exp}(\frac{-1.51714\times \mathrm{exp}(\mathrm{xb})}{0.0076927}))]}{0.09165}$$

where xb =  − 19.907 − 0.0336 × Albumin + 0.0095 × Creatinine + 0.1953 × Glucose + 0.0954 × LnCRP-0.0120 × Lymphocyte Percent + 0.0268 × Mean Cell Volume + 0.3306 × Red Cell Distribution Width + 0.00188 × Alkaline Phosphatase + 0.0554 × White Blood Cell Count + 0.0804 × Chronological Age.

The measurement method of telomere length has been reported elsewhere [26]. Briefly, blood samples were collected and determined by quantitative polymerase chain reaction (qPCR) assays to assess the telomere length relative to standard reference DNA (the T/S ratio). All sample analyses were performed in duplicate strictly following the manufacturer’s instruction and all results were checked to meet the laboratory’s standardized criteria for acceptability before being released for reporting.

Statistical analysis

All analyses were performed with SAS (version 9.4) or R (version 3.6.3) and accounted for the complex sampling design of the NHANES, and the code for analyses were in Additional file 2. Chi-square tests and t-tests were used to assess participants’ demographic characteristics by OA status. The metal concentrations were Ln-transformed to obtain approximate normal distribution (continuous variables) or categorized into four quartiles (Q1, Q2, Q3, and Q4) as categorical variables. Multivariable logistic regression was applied to estimate odds ratios (OR) and corresponding 95% confidence interval (CI) for the associations of metals and biological aging markers with OA risk. Multivariable linear regression was used to explore the associations of metals with biological aging markers. All analyses adjusted covariates for age (< 60 or ≥ 60 years), sex (male or female), race/ethnicity (Mexican American, other Hispanic, non-Hispanic White, non-Hispanic Black, or others), marital status (married/cohabiting, widowed/divorced/separated, or never married), physical activity (moderate or vigorous), drinking alcohol status (ever or never), body mass index (BMI), the ratio of family income to poverty (PIR), and serum cotinine concentration. Missing data for the covariates were coded as a missing indicator category for categorical variables and were imputed with the median for the continuous variable. We used the false-discovery rate (FDR) correction to adjust for multiple tests in the regression models. The trend test across increasing exposure groups was calculated using integer values (1, 2, 3, and 4). Pearson correlation analysis was used to assess the correlations among Ln-transformed metals. The directed acyclic graphs (R package “dagitty” and “ggdag”) showed the main relationships among exposures, outcomes, covariates, and mediators (Additional file 1: Fig. S1).

Weighted quantile sum (WQS) regression was applied to explore the overall effects of metals on OA as it performed well in characterizing environmental mixtures [27]. R package (“gWQS”) can empirically calculate the WQS index comprised of weighted sums of individual metal concentrations. The WQS index (ranged from 0 to 1) represented the mixed exposure level of metals, and the components of concern were identified by non-negligible weights. The final result was interpreted as the simultaneous effect on OA of a one-quantile increase of mixed metals.

Several sensitivity analyses were also performed. First, given the potential non-linearity and non-additive relations among metals, Bayesian kernel machine regression (BKMR) was used to assess the joint effect of all metals (“BKMR” packages) and the dose–response relationship between single metal and OA risk when fixing other metal concentrations. Second, we further adjusted for the occupation of participants (industry and agriculture, service and transportation, and others). Third, we further adjusted for age-related diseases including diabetes (yes or no), hypertension (yes or no), cardiovascular disease (yes or no) and cancer (yes or no), and OA-related medicine use (e.g., meloxicam, celecoxib) (yes or no). Fourth, we further adjusted for the survey cycle (1, 2, 3, 4, 5, 6, 7, 8, and 9). Fifth, participants with abnormal urinary creatinine were excluded. Finally, pregnant participants were excluded.

The potential mediating effects of biological aging markers on the associations of single and mixed metals (shown as WQS index) with OA risk were estimated by two mediation models, including parallel mediation (R package “mediation”) and serial mediation analyses (R package “bruceR”). Parallel mediation models used individual indicators as a mediator, while serial mediation models used a pathway as a mediator. Mediation analyses used the quasi-Bayesian Monte Carlo method with 1000 simulations based on normal approximation. The direct effect (DE) represented the effects of metal exposure on OA without a mediator. The indirect effect (IE) represented the effects of metal exposure on OA through the mediator. The proportion of mediation was calculated by using IE divided by TE (total effect).

Results

Characteristics of participants and metals distribution

Among the 12,584 adults, 1356 were diagnosed with OA. The demographic characteristics of the study participants with or without OA are listed in Table 1. Overall, age, sex, ethnicity, marital status, physical activity, drinking alcohol status, family PIR, BMI, serum cotinine, biological age, phenotypic age, and telomere length were statistically significant between OA and non-OA participants.

Table 1 Basic characteristics of participants by OA among U.S. adults, NHANES 1999–2016

The distribution of metal concentrations is shown in Additional file 1: Table S1. The detection rates of metals were greater than 75%. Pearson’s coefficient between Ln-transformed metals showed moderate correlations between Cs and Tl (r = 0.58), Ba and Co (r = 0.41), and Cd and Pb (r = 0.40), while other correlations were relatively poor (Additional file 1: Fig. S2).

Associations between metal concentration and OA risk

Figure 2 displayed the associations between creatinine-adjusted metal concentrations and OA risk by the survey-weighted logistic regression models. The highest exposure quantile of Cd (OR 1.64, 95%CI 1.20 to 2.23), Co (OR 1.59, 95%CI 1.20 to 2.10), and Cs (OR 1.48, 95%CI 1.13 to 1.93) increased the risk of OA compared to quantile 1 (all P for trend < 0.05). These associations were also presented in Ln-transformed metals and OA risk (all P < 0.05).

Fig. 2
figure 2

OR (95% CI) in OA associated with single and mixed metals levels. Models were adjusted for sex, age, race/ethnicity, education, family income-to-poverty ratio, marital status, body mass index, physical activity, drinking alcohol status, and serum cotinine. Continuous, Ln-transformed concentration of metals; Q, quartile; OA, osteoarthritis. All P for trend were FDR-adjusted

The WQS index of mixed metals was positively associated with OA risk (OR 1.23, 95%CI 1.10 to 1.37) (Fig. 2). Additionally, the highest weighted metals in WQS models were Cd (54.45%), Co (27.14%), and Cs (9.23%) (Additional file 1: Fig. S3).

In the sensitivity analyses, the BKMR model also showed a significant and positive association of mixed metals with OA risk (Additional file 1: Fig. S4). After further adjusting for the occupation of the participants, adjusting for diabetes, hypertension, cardiovascular disease, cancer, and OA-related medicine use, adjusting for the survey cycle, excluding participants with abnormal urinary creatinine, or excluding pregnant participants, the results were not materially changed (Additional file 1: Table S2-6).

Associations between metal concentration and biological aging markers

Figure 3 showed the associations of metals with biological aging markers based on the linear regression. We found the highest quartile of Cd, Co, Cs, Pb, and Tl (versus quartile 1) were associated with a growing of biological age (all P for trend < 0.001). With increasing quantiles of Cd, Co, Cs, and Pb, phenotypic age was increased (all P for trend < 0.001). Cd and Cs were negatively associated with telomere length (all P for trend < 0.05). Moreover, mixed metals had positive associations with biological age (β 4.91, 95%CI 4.52 to 5.31) and phenotypic age (β 5.90, 95%CI 5.45 to 6.34) and had a negative association with telomere length (β − 0.04, 95%CI − 0.05 to − 0.02).

Fig. 3
figure 3

Regression coefficients (95% CI) in aging markers associated with single and mixed metal levels. Models were adjusted for sex, age, race/ethnicity, education, family income-to-poverty ratio, marital status, body mass index, physical activity, drinking alcohol status, and serum cotinine. Continuous, Ln-transformed concentration of metals; Q, quartile. All P for trend were FDR-adjusted

Associations between biological aging markers and OA risk

Table 2 shows the associations of aging markers with OA risk based on the logistic regression. Each 1-year increase in biological age increased the risk of OA by 6% (95%CI 1.05 to 1.07). Similarly, each 1-year increase in phenotypic age was associated with an increased OA risk (OR 1.04, 95%CI 1.04 to 1.05). Additionally, for each unit increase in telomere length (mean T/S ratio), OR for OA decreased by 74% (95%CI 0.07 to 0.97), which was consistent with the quantile analysis (Q4 vs. Q1: OR 0.40, 95%CI 0.17 to 0.90).

Table 2 Associations of aging makers with OA risk, NHANES 1999–2016

Mediation analyses

Furthermore, parallel mediation analyses were performed to evaluate the potential mediation effects of biological aging on the associations of metals with OA risk. Biological age had significant mediated effects on the associations of Cd, Co, and Cs with OA risk, and the proportion of mediation was 69.27%, 18.43%, and 30.50% respectively (all P < 0.05). The mediated proportion of phenotypic age was 61.01%, 31.83%, and 17.00% on the associations of Cd, Co, and Cs with OA risk respectively. Telomere length also mediated the association between Cs and OA risk and the proportion of mediation was 9.81% (Additional file 1: Table S7). Moreover, biological age, phenotypic age, and telomere length parallelly mediated the associations of mixed metals with OA risk with 57.60%, 48.30%, and 9.50% proportion of mediation respectively (all P < 0.05) (Fig. 4).

Fig. 4
figure 4

Estimated proportion of the association between urinary mixed metals and OA mediated by biological age (A), phenotypic age (B), and telomere length (C) and serially mediated by telomere length-biological age path (D) and telomere length-phenotypic age path (E). Models were adjusted for sex, age, race/ethnicity, education, family income-to-poverty ratio, marital status, body mass index, physical activity, drinking alcohol status, and serum cotinine. IE, the estimate of the indirect effect; DE, the estimate of the direct effect; Proportion of mediation = IE/DE + IE

As cell senescence is a leading cause of whole-body aging, we further explored the potential paths of associations between metals and OA risk by the serial mediation models. Additional file 1: Table S8 shows the association between Cs and OA risk serially mediated by telomere length-biological age path and telomere length-phenotypic age path, with 4.55% and 4.17% proportion of mediation respectively. Besides, the serial mediated effect of telomere length-biological age path and telomere length-phenotypic age path on the association of mixed metals with OA risk were 0.70% and 0.60% respectively. And the proportion of mediation was 11.67% and 9.84%, respectively (Fig. 4).

Discussion

In the current study, we mainly provided two novel findings based on the U.S. general population. First, the single and mixed exposure models consistently emphasized Cd, Co, and Cs as strong risk factors for OA. Furthermore, biological aging was identified as a mediator for the positive associations of metals, especially mixed metals, with OA risk.

Most toxic metals are widely distributed in industrial products, soil, food, and drinking water and mainly accumulated in the liver, kidney, muscle, and bone [28, 29]. Compared to patients who suffered other bone injuries, OA patients had a higher heavy metal concentration in bone, such as Cd and Co [30]. Moreover, the toxicological study showed that Cd exposure rendered bone or joint-related injuries, such as cartilage loss [9]. Our study firstly provided epidemiological evidence on the positive and significant associations of Cd, Co, and Cs with OA risk. Similarly, our recent cross-sectional study suggested a significant association between Cd and rheumatoid arthritis risk, further indicating Cd was responsible for musculoskeletal disorders [14]. Studies of other models also provided supportive evidence of metals causing inflammation. For example, Cd, Co, or Cs exposure triggered inflammation with the upregulation of IL-6 in mice brain, heart, or macrophages [31,32,33].

Noteworthy, heavy metals usually coexist in the environment and their effects rely on their cooperation and interaction. Previous studies reported that the mixture of Cd and Pb further diminished skeletal growth and caused osteopenia compared to single metals in rats. The mixed exposure is measured by WQS with an empirically weighted index in the direction of the joint action, in which components of concern are identified by non-negligible weights. This method truly reflects the joint action of mixed exposure. In our study, Cd, Co, and Cs were identified as the strongest risk factors for OA. Also, mixed metals were positively associated with OA risk in WQS models and BKMR models, further suggesting that metal exposure may promote OA progress.

Mechanistically, our research interest focused on biological aging as a growing number of studies emphasized the role of metals in age-related diseases [34]. Our results showed that both single (Cd, Co, Cs, Pb, and Tl) and mixed metals were positively associated with biological aging. Consistent toxicological evidence presented that heavy metals can induce senescence in stromal and vascular smooth muscle cells (e.g., Co and Cd) [35, 36], as well as accelerate brain aging in animals, shown by damaged molecules and impaired energy metabolism [17]. Also, a recent population study indicated the positive associations between heavy metals (e.g., Co) and biological aging marker (DNA methylation age) [37]. Many investigations on metals and telomere length have reached similar results [38]. In addition, heavy metals may regulate senescence-related pathways, including Keap1-NRF2, nuclear translocation of NF-κB, and the high expression of HIF [11]. Collectively, the abovementioned evidence supported our finding that heavy metals may modulate aging.

On the other hand, aging is one of the crucial drivers for chronic diseases, especially OA. Aging promoted OA progress mainly through the accumulation of senescent cells, senescence-related secretory phenotype (SASP) production, and ROS-induced DNA damage [19]. Herein, we measured cellular senescence and whole-body aging using different markers and consistently found positive associations between biological aging and OA risk. Specifically, phenotypic age and biological age incorporate composite clinical biomarkers, such as inflammation (e.g., CRP), immunity (e.g., number of immune cells), and organ function (e.g., albumin and serum urea nitrogen). Both of which are useful to identify the unhealthy condition of people caused by multiple system disorders [39]. Interestingly, OA is a whole-joint disease accompanied by system inflammation, indicating the significance of whole-body aging in OA. Additionally, telomere wear is widely used to evaluate cell senescence [40]. Abundant epidemiological studies have reported the inverse associations of telomere length with OA risk, which was consistent with our findings [41]. Similarly, experimental studies reported that targeting the elimination of senescent chondrocytes or inhibiting SASP improved OA-related outcomes [21, 42]. Overall, aging is an important biological mechanism of OA onset and progress.

According to the mentioned findings, we further performed mediation analyses. We observed significant mediated effects of aging markers on the associations of metals with OA risk. Notably, three biological aging markers all mediated the positive association between Cs and OA risk, among which the mediated proportion of biological age was 30.50%, phenotypic age was 17.00%, and telomere length was 9.81%. In fact, Cs radiation has long been considered an inducer of premature aging in animals and can cause bone loss [43, 44]. What is more, the associations of mixed metals with OA risk are also mediated by these aging markers. More importantly, the mediated proportion of whole-body aging was more than 50%, indicating it was largely involved in the promoting process of metals on OA.

Interestingly, whole-body aging is mainly driven by cell senescence [45]. For example, transplantation of senescent cells into young mice caused a decline in whole-body function [46]. Inversely, anti-senescence therapeutics (e.g., senolytics and senomorphics) offered good effects on age-related diseases, such as idiopathic pulmonary fibrosis and diabetic kidney disease [47, 48]. Hence, we further explored the underlying path from metals to OA. Serial mediation analyses showed the associations of Cs and mixed metals with OA were mediated by telomere length-biological age path, and telomere length-phenotypic age path, with the mediated proportion ranging from 4.17 to 11.67%. These findings suggested that metals may accelerate cell senescence to lead to overall disorder of the body and finally aggravate OA progress.

This study has some strengths. First, we applied a variety of methods to explore the relationships between metals and OA risk in a relatively large population. What is more, we measured biological aging status from multiple perspectives. There are also limitations. First, a self-reported OA diagnosis may decrease the credibility of the results. Second, single-point measurements of metals did not fully reflect individual real exposure levels. Third, we only analyzed the correlation between 9 metals and OA risk, and further analyses were needed to identify other metals associated with OA risk. Fourth, residual and unmeasured confounders and measurement errors may bias our analyses. Fifth, although we adjusted for the survey cycle, the time span of our analysis was too wide to potentially bias the results. Finally, we conducted mediation analysis in a cross-sectional study design and it was difficult to infer a causal association. Given the limitations in the current study, these results need to be interpreted with caution and further investigations are needed to support our findings.

Conclusions

In summary, our results presented the positive associations of single and mixed metals with increased OA risk, primarily driven by Cd, Co, and Cs. Additionally, we found that metal exposure was related to aging, and aging was related to OA risk. Furthermore, mediation analyses showed that the associations of metals with OA risk may be mediated by biological aging and serially mediated by cell senescence-whole-body aging path. These findings identified risk factors for OA and indicated aging as an underlying mechanism of metal’s adverse effect on OA.

Availability of data and materials

The National Health and Nutrition Examination Survey (NHANES) data are publicly available at https://www.cdc.gov/nchs/nhanes/index.htm.

Abbreviations

Ba:

Barium

Cd:

Cadmium

CI:

Confidence interval

Co:

Cobalt

CRP:

C-reactive protein

Cs:

Cesium

ICP-MS:

Inductively coupled plasma mass spectrometry

Mo:

Molybdenum

NCHS:

National Center for Health Statistics

NHANES:

National Health and Nutrition Examination Survey

OA:

Osteoarthritis

OR:

Odds ratios

Pb:

Lead

qPCR:

Uantitative polymerase chain reaction

SASP:

Senescence-related secretory phenotype

Sb:

Antimony

SE:

Standard error

Tl:

Thallium

Tu:

Uranium

WQS:

Weighted quantile sum

References

  1. Hunter DJ, Bierma-Zeinstra S. Osteoarthritis Lancet. 2019;393(10182):1745–59.

    CAS  PubMed  Google Scholar 

  2. Safiri S, Kolahi AA, Smith E, Hill C, Bettampadi D, Mansournia MA, Hoy D, Ashrafi-Asgarabad A, Sepidarkish M, Almasi-Hashiani A, et al. Global, regional and national burden of osteoarthritis 1990–2017: a systematic analysis of the Global Burden of Disease Study 2017. Ann Rheum Dis. 2020;79(6):819–28.

    Article  PubMed  Google Scholar 

  3. Cisternas MG, Murphy L, Sacks JJ, Solomon DH, Pasta DJ, Helmick CG. Alternative methods for defining osteoarthritis and the impact on estimating prevalence in a US population-based survey. Arthritis Care Res (Hoboken). 2016;68(5):574–80.

    Article  Google Scholar 

  4. Hootman JM, Helmick CG. Projections of US prevalence of arthritis and associated activity limitations. Arthritis Rheum. 2006;54(1):226–9.

    Article  PubMed  Google Scholar 

  5. Clemens S, Ma JF. Toxic heavy metal and metalloid accumulation in crop plants and foods. Annu Rev Plant Biol. 2016;67:489–512.

    Article  CAS  PubMed  Google Scholar 

  6. Kavamura VN, Esposito E. Biotechnological strategies applied to the decontamination of soils polluted with heavy metals. Biotechnol Adv. 2010;28(1):61–9.

    Article  CAS  PubMed  Google Scholar 

  7. Soomro MH, Baiz N, Huel G, Yazbeck C, Botton J, Heude B, Bornehag CG, Annesi-Maesano I. Exposure to heavy metals during pregnancy related to gestational diabetes mellitus in diabetes-free mothers. Sci Total Environ. 2019;656:870–6.

    Article  CAS  PubMed  Google Scholar 

  8. Lim JT, Tan YQ, Valeri L, Lee J, Geok PP, Chia SE, Ong CN, Seow WJ. Association between serum heavy metals and prostate cancer risk - a multiple metal analysis. Environ Int. 2019;132:105109.

    Article  CAS  PubMed  Google Scholar 

  9. Reyes-Hinojosa D, Lozada-Pérez CA, Zamudio Cuevas Y, López-Reyes A, Martínez-Nava G, Fernández-Torres J, Olivos-Meza A, Landa-Solis C, Gutiérrez-Ruiz MC, Rojas Del Castillo E, et al. Toxicity of cadmium in musculoskeletal diseases. Environ Toxicol Pharmacol. 2019;72:103219.

    Article  CAS  PubMed  Google Scholar 

  10. Akbal A, Resorlu H, Savas Y. Toxic effects of heavy metals on bone tissue. Turk Osteoporoz Dergisi-Turkish J Osteoporosis. 2015;21(1):30–3.

    Article  Google Scholar 

  11. Paithankar JG, Saini S, Dwivedi S, Sharma A, Chowdhuri DK. Heavy metal associated health hazards: an interplay of oxidative stress and signal transduction. Chemosphere. 2021;262:128350.

    Article  CAS  PubMed  Google Scholar 

  12. YessicaEduviges ZC, Martínez-Nava G, Reyes-Hinojosa D, Mendoza-Soto L, Fernández-Torres J, López-Reyes A, Olivos-Meza A, Armienta-Hernández MA, Ruíz-Huerta EA, de Jesús G-G, et al. Impact of cadmium toxicity on cartilage loss in a 3D in vitro model. Environ Toxicol Pharmacol. 2020;74:103307.

    Article  CAS  Google Scholar 

  13. Park S, Choi NK. The relationships of blood lead level, body mass index, and osteoarthritis in postmenopausal women. Maturitas. 2019;125:85–90.

    Article  CAS  PubMed  Google Scholar 

  14. Chen L, Sun Q, Peng S, Tan T, Mei G, Chen H, Zhao Y, Yao P, Tang Y. Associations of blood and urinary heavy metals with rheumatoid arthritis risk among adults in NHANES, 1999–2018. Chemosphere. 2022;289:133147.

    Article  CAS  PubMed  Google Scholar 

  15. Li X, Li R, Yan J, Song Y, Huo J, Lan Z, Chen J, Zhang L. Co-exposure of cadmium and lead on bone health in a southwestern Chinese population aged 40–75 years. J Appl Toxicol. 2020;40(3):352–62.

    Article  CAS  PubMed  Google Scholar 

  16. Luo H, Gu R, Ouyang H, Wang L, Shi S, Ji Y, Bao B, Liao G, Xu B. Cadmium exposure induces osteoporosis through cellular senescence, associated with activation of NF-κB pathway and mitochondrial dysfunction. Environ Pollut. 2021;290:118043.

    Article  CAS  PubMed  Google Scholar 

  17. Ijomone OM, Ifenatuoha CW, Aluko OM, Ijomone OK, Aschner M. The aging brain: impact of heavy metal neurotoxicity. Crit Rev Toxicol. 2020;50(9):801–14.

    Article  CAS  PubMed  Google Scholar 

  18. Zota AR, Needham BL, Blackburn EH, Lin J, Park SK, Rehkopf DH, Epel ES. Associations of cadmium and lead exposure with leukocyte telomere length: findings from National Health and Nutrition Examination Survey, 1999–2002. Am J Epidemiol. 2015;181(2):127–36.

    Article  PubMed  Google Scholar 

  19. Xie J, Wang Y, Lu L, Liu L, Yu X, Pei F. Cellular senescence in knee osteoarthritis: molecular mechanisms and therapeutic implications. Ageing Res Rev. 2021;70:101413.

    Article  CAS  PubMed  Google Scholar 

  20. Guillén R, Otero F, Mosquera A, Vázquez-Mosquera M, Rego-Pérez I, Blanco FJ, Fernández JL. Association of accelerated dynamics of telomere sequence loss in peripheral blood leukocytes with incident knee osteoarthritis in Osteoarthritis Initiative cohort. Sci Rep. 2021.

  21. Jeon OH, Kim C, Laberge RM, Demaria M, Rathod S, Vasserot AP, Chung JW, Kim DH, Poon Y, David N, et al. Local clearance of senescent cells attenuates the development of post-traumatic osteoarthritis and creates a pro-regenerative environment. Nat Med. 2017;23(6):775–81.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  22. March LM, Schwarz JM, Carfrae BH, Bagge E. Clinical validation of self-reported osteoarthritis. Osteoarthritis Cartilage. 1998;6(2):87–93.

    Article  CAS  PubMed  Google Scholar 

  23. Klemera P, Doubal S. A new approach to the concept and computation of biological age. Mech Ageing Dev. 2006;127(3):240–8.

    Article  PubMed  Google Scholar 

  24. Levine M. Modeling the rate of senescence: can estimated biological age predict mortality more accurately than chronological age? J Gerontol A Biol Sci. 2013;68(6):667–74.

    Article  Google Scholar 

  25. Levine ME, Lu AT, Quach A, Chen BH, Assimes TL, Bandinelli S, Hou L, Baccarelli AA, Stewart JD, Li Y, et al. An epigenetic biomarker of aging for lifespan and healthspan. Aging (Albany NY). 2018;10(4):573–91.

    Article  Google Scholar 

  26. Banach M, Mazidi M, Mikhailidis DP, Toth PP, Jozwiak J, Rysz J, Watts GF. Association between phenotypic familial hypercholesterolaemia and telomere length in US adults: results from a multi-ethnic survey. Eur Heart J. 2018;39(40):3635–40.

    Article  CAS  PubMed  Google Scholar 

  27. 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 Agric Biol Environ Stat. 2015;20(1):100–20.

    Article  PubMed  Google Scholar 

  28. Satarug S, Garrett SH, Sens MA, Sens DA. Cadmium, environmental exposure, and health outcomes. Environ Health Perspect. 2010;118(2):182–90.

    Article  CAS  PubMed  Google Scholar 

  29. Rehman K, Fatima F, Waheed I, Akash MSH. Prevalence of exposure of heavy metals and their impact on health consequences. J Cell Biochem. 2018;119(1):157–84.

    Article  CAS  PubMed  Google Scholar 

  30. Scimeca M, Feola M, Romano L, Rao C, Gasbarra E, Bonanno E, Brandi ML, Tarantino U. Heavy metals accumulation affects bone microarchitecture in osteoporotic patients. Environ Toxicol. 2017;32(4):1333–42.

    Article  CAS  PubMed  Google Scholar 

  31. Guo K, Ge J, Zhang C, Lv MW, Zhang Q, Talukder M, Li JL. `Cadmium induced cardiac inflammation in chicken (Gallus gallus) via modulating cytochrome P450 systems and Nrf2 mediated antioxidant defense. Chemosphere. 2020;249:125858.

    Article  CAS  PubMed  Google Scholar 

  32. Wang J, Pan H, Lin Z, Xiong C, Wei C, Li H, Tong F, Dong X. Neuroprotective effect of fractalkine on radiation-induced brain injury through promoting the M2 polarization of microglia. Mol Neurobiol. 2021;58(3):1074–87.

    Article  CAS  PubMed  Google Scholar 

  33. Chen F, Chen R, Liu H, Sun R, Huang J, Huang Z, Jian G. BMP-7 ameliorates cobalt alloy particle-induced inflammation by suppressing Th17 responses. APMIS. 2017;125(10):880–7.

    Article  CAS  PubMed  Google Scholar 

  34. Adlard PA, Bush AI. Metals and Alzheimer’s disease: how far have we come in the clinic? J Alzheimers Dis. 2018;62(3):1369–79.

    Article  PubMed  PubMed Central  Google Scholar 

  35. Pan CH, Chen CJ, Shih CM, Wang MF, Wang JY, Wu CH. Oxidative stress-induced cellular senescence desensitizes cell growth and migration of vascular smooth muscle cells through down-regulation of platelet-derived growth factor receptor-beta. Aging (Albany NY). 2019;11(19):8085–102.

    Article  CAS  Google Scholar 

  36. Hodjat M, Rezvanfar MA, Abdollahi M. A systematic review on the role of environmental toxicants in stem cells aging. Food Chem Toxicol. 2015;86:298–308.

    Article  CAS  PubMed  Google Scholar 

  37. Xiao L, Zan G, Feng X, Bao Y, Huang S, Luo X, Xu X, Zhang Z, Yang X. The associations of multiple metals mixture with accelerated DNA methylation aging. Environ Pollut. 2021;269:116230.

    Article  CAS  PubMed  Google Scholar 

  38. Møller P, Wils RS, Jensen DM, Andersen MHG, Roursgaard M. Telomere dynamics and cellular senescence: an emerging field in environmental and occupational toxicology. Crit Rev Toxicol. 2018;48(9):761–88.

    Article  PubMed  CAS  Google Scholar 

  39. Liu Z, Kuo PL, Horvath S, Crimmins E, Ferrucci L, Levine M. A new aging measure captures morbidity and mortality risk across diverse subpopulations from NHANES IV: a cohort study. PLoS Med. 2018;15(12): e1002718.

    Article  PubMed  PubMed Central  Google Scholar 

  40. Shay JW. Telomeres and aging. Curr Opin Cell Biol. 2018;52:1–7.

    Article  CAS  PubMed  Google Scholar 

  41. Kuszel L, Trzeciak T, Richter M, Czarny-Ratajczak M. Osteoarthritis and telomere shortening. J Appl Genet. 2015;56(2):169–76.

    Article  CAS  PubMed  Google Scholar 

  42. Coryell PR, Diekman BO, Loeser RF. Mechanisms and therapeutic implications of cellular senescence in osteoarthritis. Nat Rev Rheumatol. 2021;17(1):47–57.

    Article  PubMed  Google Scholar 

  43. Nguyen HQ, To NH, Zadigue P, Kerbrat S, De La Taille A, Le Gouvello S, Belkacemi Y. Ionizing radiation-induced cellular senescence promotes tissue fibrosis after radiotherapy. A review Crit Rev Oncol Hematol. 2018;129:13–26.

    Article  PubMed  Google Scholar 

  44. Kondo H, Searby ND, Mojarrab R, Phillips J, Alwood J, Yumoto K, Almeida EA, Limoli CL, Globus RK. Total-body irradiation of postpubertal mice with (137)Cs acutely compromises the microarchitecture of cancellous bone and increases osteoclasts. Radiat Res. 2009;171(3):283–9.

    Article  CAS  PubMed  Google Scholar 

  45. Muñoz-Espín D, Serrano M. Cellular senescence: from physiology to pathology. Nat Rev Mol Cell Biol. 2014;15(7):482–96.

    Article  PubMed  CAS  Google Scholar 

  46. Xu M, Pirtskhalava T, Farr JN, Weigand BM, Palmer AK, Weivoda MM, Inman CL, Ogrodnik MB, Hachfeld CM, Fraser DG, et al. Senolytics improve physical function and increase lifespan in old age. Nat Med. 2018;24(8):1246–56.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  47. Chen L, Mei G, Jiang C, Cheng X, Li D, Zhao Y, Chen H, Wan C, Yao P, Gao C, et al. Carbon monoxide alleviates senescence in diabetic nephropathy by improving autophagy. Cell Prolif. 2021;54(6):e13052.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  48. Justice JN, Nambiar AM, Tchkonia T, LeBrasseur NK, Pascual R, Hashmi SK, Prata L, Masternak MM, Kritchevsky SB, Musi N, et al. Senolytics in idiopathic pulmonary fibrosis: results from a first-in-human, open-label, pilot study. EBioMedicine. 2019;40:554–63.

    Article  PubMed  PubMed Central  Google Scholar 

Download references

Acknowledgements

The authors thank all participants and all investigators.

Funding

This study was supported by the National Key Research and Development Project (2018YFC1604000).

Author information

Authors and Affiliations

Authors

Contributions

All authors read and approved the final manuscript. LC, FqL, and TqT conducted analyses. LC, YZ, and HmC wrote the draft of the article. PY and YhT conceived of the study design. All authors contributed to the interpretation of the results and critical revision of the manuscript for important intellectual content and approved the final version of the manuscript. YhT is the guarantor of this work and, as such, had full access to all the data in the study and takes responsibility for the integrity of the data and the accuracy of the data analysis.

Corresponding author

Correspondence to Yuhan Tang.

Ethics declarations

Ethics approval and consent to participate

The protocols of NHANES were approved by the institutional review board of the National Center for Health Statistics, CDC (Protocol #98–12 #2005–06 #2011–17). Written informed consent was obtained from each participant before participation in this study.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

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

Supplementary Information

Additional file 1:

Fig. S1. Directed acyclic graph. Fig. S2. Pearson’s correlation matrix. Fig. S3. Weighted values of urinary metals for OA in WQS models. Fig. S4. Associations of the urinary metals with OA risk estimated by Bayesian Kernel Machine Regression (BKMR). Table S1. Distributions of metals in the study population. Table S2. OR (95% CI) in OA associated with single and mixed urinary metals levels with further adjustment for occupation. Table S3. OR (95% CI) in OA associated with single and mixed urinary metals levels with further adjustment for other diseases and medicine use. Table S4. OR (95% CI) in OA associated with single and mixed urinary metals levels with further adjustment for survey cycle. Table S5. OR (95% CI) in OA associated with single and mixed urinary metals levels after excluding participants with abnormal urinary creatinine. Table S6. OR (95% CI) in OA associated with single and mixed urinary metals levels after excluding pregnant participants. Table S7. Biological aging markers as mediators in the associations of single metals with OA risk. Table S8. Telomere length and biological age/biological age as serial mediators in the associations of single metals with OA risk.

Additional file 2.

Code for analyses.

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

Chen, L., Zhao, Y., Liu, F. et al. Biological aging mediates the associations between urinary metals and osteoarthritis among U.S. adults. BMC Med 20, 207 (2022). https://doi.org/10.1186/s12916-022-02403-3

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12916-022-02403-3

Keywords