Study population
The Centers for Disease Control and Prevention (CDC) conducts the National Health and Nutrition Examination Survey (NHANES), a nationally representative survey of the noninstitutionalized U.S. civilian population. NHANES includes interviews, physical examinations, and laboratory tests. Extensive data were collected through interviews, including demographic, socioeconomic, dietary and health-related questionnaire data (37). Laboratory tests and medical, dental and physiological examinations are performed by trained medical personnel. The NHANES dataset is publicly available at (https://www.cdc.gov/nchs/nhanes/index.htm).
Our study analyzed NHANES data from 2007 to 2016, during which serum α-Klotho levels were measured. We combined five NHANES cycles (2007–2016) with 50588 subjects (37). Among 50588 subjects, serum α-Klotho levels were analyzed in participants aged 40–79 (n = 13764). Of these participants, only 11036 had simultaneously measured blood lead and mercury levels. The final study population included 11032 participants (Fig. 1).
The NHANES research protocol was approved by the National Center for Health Statistics Institutional Review Board (https://www.cdc.gov/nchs/nhanes/irba98.htm). Additionally, all participants provided verbal and written consent for future research.
Determination of serum α-Klotho levels
Serum α-Klotho levels were analyzed in frozen serum samples aged 40–79 collected from the NHANES database between 2007 and 2016. Fresh-frozen serum samples, stored at − 80°C, were analyzed by the Centers for Disease Control and Prevention to the Northwest Lipid Metabolism and Diabetes Research Laboratories, Division of Metabolism, Endocrinology, and Nutrition, University of Washington. Serum α-Klotho levels were detected using an ELISA kit produced by Japan’s IBL International Company. The average of two repeated analyses was used for quality assurance, and samples with duplicate results that differed by more than 10% were reanalyzed. The resulting detection sensitivity was 4.33 pg/mL. Detection was performed in 114 samples from healthy volunteers, with a reference range of 285.8–1638.6 pg/ml and an average value of 698.0 pg/mL (38).
Determination of blood lead and blood mercury levels
All determination results were completed at the National Center for Environmental Health, and Centers for Disease Control and Prevention for analysis in accordance with the Laboratory Procedures Manual. Mass spectrometry was used to determine lead and mercury, and inductively coupled plasma mass spectrometry (ICP-MS) was used to determine lead and mercury levels in blood. The lowest detection limits for lead and mercury are 0.07ug/dL and 0.28ug/L respectively. The results field for analytes with values below the lower limits of detection were filled with an imputed fill value equal to the lower limit of detection divided by the square root of 2 (LLOD/sqrt [2]) (0.05 µg/dL for blood lead, 0.14 µg/L for blood mercury). From 2013 onwards, for participants 12 years and older, special sample weights were created for the subsample. These special weights accounted for the additional probability of selection into the subsample, as well as the additional nonresponse to these lab tests. Therefore, if a participant 12 years and older was selected as part of the one-half subsample, but did not provide a blood specimen, he/she would have the sample weight value assigned as “0” in his/her record. The accuracy and precision of all detection methods were controlled to adhere to the Scientific Quality Control and Quality Assurance Performance Standards for Environmental Health Laboratories (39).
Renal function assessment
Renal function assessment in this study was mainly performed using two parameters. Namely, eGFR and urinary albumin-to-creatinine ratio (UACR), the main focus of this study was eGFR, which was determined by the 2021 Chronic Kidney Disease Epidemic Collaboration equation (CKD-EPI) (40), which relies on serum creatinine levels. UACR was used as the outcome measure for sensitivity analysis, and its calculation formula (41) is:
$${UACR}_{\left(\frac{mg}{g}\right)}=\frac{{urinary albumin level}_{\left(\frac{mg}{dL}\right)}}{{urinary creatinine level}_{\left(\frac{g}{dL}\right)}}$$
Urinary albumin and creatinine levels were obtained from urine samples.
Covariates
Covariate data included demographic questionnaire data, socioeconomic status, lifestyle, and habit variables. Age is divided into two groups (40–59 years, 60–79 years). Body mass index is calculated by dividing weight (kilograms) by the square of height (meters) (18.5–25 is normal, 25–30 is overweight, and > 30 is obese). Race is divided into White, Black, Mexican, and other races. Education Degree is divided into Less than high school (less than 9th, High school (9–11th grade [including 12th with no diploma], high school or equivalent), and College or higher (some college or associate’s degree or college graduate or above). Sports activities are divided into Vigorous, Moderate, and low).
Health status included hypertension, defined as mean systolic blood pressure (SBP) ≥ 140 mmHg and/or mean diastolic blood pressure (DBP) ≥ 90 mmHg, or self-reported diagnosis of hypertension and taking hypertension medication. Diabetes is defined as (1) Doctor diagnosed diabetes, (2) Glycated hemoglobin > 6.5%, (3) Fasting blood glucose ≥ 7.0mmol/L, (4) Random blood glucose ≥ 11.1mmol/L, (5) Glucose tolerance test 2-hour blood glucose ≥ 11.1mmol /L, (6) Use diabetes drugs or insulin. Atherosclerotic heart disease refers to a history of coronary heart disease, angina, stroke, and heart attack.
Statistical analysis
First, α-Klotho, blood lead, blood mercury, and eGFR were log-transformed to reduce skewness. Categorical variables were assessed by calculation of (percent), while continuous variables were determined by (M ± SD), and nonparametric continuous variables were analyzed by calculating (IQR). Pearson correlation coefficients were calculated for eGFR, blood lead, blood mercury, and Klotho levels and their corresponding log-transformed values. Linear regression analysis was used to evaluate the correlation between renal function as the dependent variable and blood lead and blood mercury as independent variables. Beta coefficients and 95% CI were calculated. Generalized additive models (GAM) were used for nonlinear relationships. Models were adjusted for age, sex, race, education, income, BMI and physical activity, hypertension, diabetes, and cardiovascular disease.
Second, we explored the potential mediating role of Klotho in the association between metals and renal function. Causal mediation analysis was used to conduct mediation analysis linear regression models for metals (blood lead, blood mercury), mediators, and outcomes. Use the bootstrap method to estimate confidence intervals. The logarithm of metals (lead, mercury) is the exposure variable, the logarithm of serum anti-aging protein α-Klotho is the mediating factor, and the logarithm of eGFR is the outcome variable. The following results were measured:
-
The effect of metals on α-Klotho;
-
The effect of α-Klotho on renal function;
-
The effect of metals on overall renal function;
-
Considering the overall impact of metals on renal function in the case of α-Klotho protein, the proportion of intermediary factors is calculated using the following formula:
$$\left(\beta total effect - \beta direct effect\right)/\beta total effect \times 100$$
Third, subgroup analysis was performed to explore whether age or gender modulates the mediating effect of serum α-Klotho on the correlation between metals and eGFR. In mediation analyses, stratification was performed by age of analysis (40–59 years, 60–79 years) or gender (male, female). Since MICE can handle variable types (continuous variables, binary, unordered categorical, ordinal categorical), we performed multiple imputations on covariate data with < 5% missing data. The imputation procedure was performed using a linear regression method for continuous variables, and an ordinal or binary logistic regression model for categorical variables, chained equations with m = 5. Use Rubin’s rules to summarize the effects of the five data sets to obtain a data set after the combined effects with the most minor error (42, 43). To verify the robustness of the results after imputation, we used sensitivity analysis to explore the role of serum α-Klotho in the relationship between metals and UACR. Based on NHANES guidelines, we used weighted estimation with MEC weights. All analytical calculation processes were performed in R software (Version 4.3.1), and two-sided p < 0.05 indicated statistical significance.