Polycyclic Aromatic Hydrocarbons and the Risk of Kidney Stones in US Adults: An Exposure-Response Analysis of NHANES 2007-2012.

BACKGROUND: Polycyclic aromatic hydrocarbons (PAHs) exposure may cause various diseases. However, the association between PAHs exposure and kidney stones remains unclear. The purpose of this study was to examine the relationship between PAHs and the risk of kidney stones in the US population.
METHODS: The study included a total of 30,442 individuals (≥20 years) from the 2007-2012 National Health and Nutrition Examination Survey (NHANES). Nine urinary PAHs were included in this study. Logistic regression and dose-response curves were used to evaluate the association between PAHs and the risk of kidney stones.
RESULTS: We selected 4385 participants. The dose-response curves showed a significant positive association between total PAHs, 2-hydroxynaphthalene, 1-hydroxyphenanthrene, 2-hydroxyphenanthrene, 9-hydroxyfluorene and the risk of kidney stones after adjusting for confounding factors. Compared with the low group, an increased risk of kidney stones was observed in the high group of total PAHs [OR (95% CI), 1.32 (1.06-1.64), P=0.013], 2-hydroxynaphthalene [OR (95% CI), 1.37 (1.10-1.71), P=0.005], 1-hydroxyphenanthrene [OR (95% CI), 1.24 (1.00-1.54), P=0.046] and 9-hydroxyfluorene [OR (95% CI), 1.36 (1.09-1.70), P=0.007].
CONCLUSION: High levels of PAHs were positively associated with the risk of kidney stones in the US population.

Introduction

Kidney stones are a common and frequently occurring disease of the urinary system, which are caused by the abnormal accumulation of crystal substances in the kidney.1 In recent years, the prevalence of kidney stones is 12% and 5% in adult men and women in western countries, respectively, and the incidence rate has increased.2 Kidney stones can cause serious complications, pain, hematuria, infection, decreased kidney function, and even kidney failure.3,4 Dietary and lifestyle factors play an important role in the risk of kidney stones, especially high animal protein intake increases the risk of kidney stone formation.5

Polycyclic aromatic hydrocarbons (PAHs) are hydrocarbon compounds containing more than two fused aromatic rings. To date, it has been found that PAHs consist of more than 200 compounds, such as naphthalene, anthracene, phenanthrene and pyrene.6,7 PAHs are generated from incomplete combustion of organic polymer compounds, including coal, petrochemicals, rubber, plastics, wood and tobacco, etc.8 PAHs have been detected in the human production and living environment.9 They are a typical class of persistent pollutants that can be rapidly enriched in organisms and cause various hazards.10 PAHs can be absorbed into the body by lungs, gut and skin, and be metabolized and eliminated in the urine within a short period.11 Therefore, urinary PAHs concentrations are considered biomarkers of exposure to PAHs.12

As a typical persistent organic pollutant, previous studies have shown that the main target organ of PAHs is the lung, but the liver, with its detoxification function, and the kidney, with its excretion function, may also be affected by PAHs exposure.13 Studies on the mechanism of damage caused by PAHs have shown that PAHs can enhance the body’s oxidative stress and increase the body’s reactive oxygen species.14–16 However, the association between PAHs and the risk of kidney stones remains unclear. We aim to investigate whether the level of PAHs in a representative population sample from the National Health and Nutrition Survey (NHANES) is associated with the risk of kidney stones.

Patients and Methods

Patients Selection

NHANES is a nationally representative research program that uses a complex stratified multistage sampling design to assess children and adults’ health and nutritional status in the United States. The project has been conducted every two years since 1999, and data are collected through household interviews and health screenings. The NHANES Institutional Review Board approved the study protocol for the project, and all participants signed informed consent forms during the survey. The study included 30,442 individuals (≥20 years) from the NHANES 2007–2008, 2009–2010, and 2011–2012 cycles. More details are available on the web ().

The exclusion criteria were as follows: (a) Patients who had not completed the survey (kidney stones KIQ026) (n=12,784); (b) unknown PAHs (n=12,439); (c) unknown family poverty ratio (n=478); (d) unknown BMI (n=51); (e) unknown creatinine (n=259); (f) unknown uric acid (n=2); (g) abnormal PAHs (n=44). After excluding all factors, 4385 participants were selected.

Study Variables and Outcome

Total PAHs are the main predictors in this study, including 1-hydroxynaphthalene, 2-hydroxynaphthalene, 3-hydroxyfluorene, 2-hydroxyfluorene, 3-hydroxyphenanthrene, 1-hydroxyphenanthrene, 2-hydroxyphenanthrene, 1-hydroxypyrene and 9-hydroxyfluorene. According to the website (wwwn.cdc.gov/nchs/data/nhanes/2011-2012/labmethods/pah_g_met.pdf), urinary PAHs were measured by standard procedures.

Besides, we included other variables, such as gender, age, race, education level, marital status, family poverty ratio, body mass index (BMI), hypertension, diabetes, recreational activities. Based on the Chronic Kidney Disease Epidemiology Collaboration equation, we calculate the estimated glomerular filtration rate (eGFR) of all participants.17

Male: GFR = 141 × min (Scr/0.9, 1)−0.411 × max (Scr/0.9, 1)−1.209 × 0.993Age × 1.159 (if black).

Female: GFR = 141 × min (Scr/0.7, 1)−0.329 × max (Scr/0.7, 1)−1.209 × 0.993Age × 1.018×1.159 (if black).

The endpoints of the study were the history of kidney stones. We extracted the relevant information from the questionnaire data file. In the questionnaire, participants responded to “Have you ever had kidney stones?”. Those who answered “No” were considered to have no history of kidney stones.

Statistical Analysis

Continuous variables were expressed as mean ± standard deviation, and t-tests for slope were used in the generalized linear model. Categorical variables were defined as n (%) and analyzed by Chi-square analysis. Nine kinds of PAHs and total PAHs were divided into two groups. Binary logistic regression models were used to evaluate the corrected odds ratios (ORs) and 95% confidence interval (CI) intervals for factors associated with kidney stones. In the extended model, based on the correlation of demographic characteristics with PAHs, we adjusted gender, age, race, education levels, marital status, family poverty, body mass index, hypertension, diabetes, recreational activities, blood urea nitrogen, uric acid, creatinine and eGFR of the participants.

The restricted cubic spline function is a powerful tool for describing dose-response relationships between continuous variables and outcomes. The restricted cubic spline function was applied to describe the dose-response relationship between PAHs and the risk of kidney stones and adjusted for the model variables.18 The study was analyzed using SPSS software (version 24.0) and RStudio software (version 1.2.5033), GraphPad software (version 8.0.2) and Origin 2018 software. P<0.05 was considered statistically significant.

Results

Participant Characteristics

A total of 30,442 NHANES participants (2007–2012) were included. Of these participants, 5165 had complete data on urinary PAHs and stone investigations. According to the exclusion criteria, 4385 participants were finally included in the study (Figure 1). The characteristics and demographics of participants were shown in Table 1. Among them, 376 participants (8.6%) had the stone formation, and 4009 participants (91.4%) had no stone formation. The chi-square test indicated that several variables were significantly different between non-stone formers and stone formers, including gender, age, race, BMI, hypertension, diabetes, vigorous recreational activities, blood urea nitrogen, creatinine, Uric acid, eGFR and total PAHs. Participants with a history of kidney stones were mostly concentrated in men (60.6%), ≥50 years (58.8%), non-Hispanic whites (59.0%), BMI≥25 kg/m2 (75.0%), hypertension-positive (43.9%), diabetes-positive (17.3%) and less vigorous recreational activities (84.8%). Moreover, stone formers had higher blood urea nitrogen, creatinine, uric acid, total PAHs and lower eGFR than none-stone formers.

Figure 1

Schematic flow diagram of inclusion and exclusion criteria for our study cohort.

Table 1

Baseline Characteristics of NHANES Participants Between 2007 and 2012 a

Characteristic	Total	None-Stone Formers	Stone Formers	P value
	No. (%)	No. (%)	No. (%)
Total participants	4385	4009 (91.4)	376 (8.6)
Gender	<0.001
Male	2196 (50.1)	1968 (49.1)	228 (60.6)
Female	2189 (49.9)	2041 (50.9)	148 (39.4)
Age	<0.001
Mean (SD)	48.57 (17.70)	48.03 (17.70)	54.32 (16.62)	<0.001
<50 years	2314 (52.8)	2159 (53.9)	155 (41.2)
≥50 years	2071 (47.2)	1850 (46.1)	221 (58.8)
Race	<0.001
Non-Hispanic white	2009 (45.8)	1787 (44.6)	222 (59.0)
Non-Hispanic black	887 (20.2)	844 (21.1)	43 (11.4)
Mexican American	649 (14.8)	610 (15.2)	39 (10.4)
Other Hispanic	449 (10.2)	399 (10.0)	50 (13.3)
Other	391 (8.0)	369 (9.2)	22 (5.9)
Education	0.516
Less than high school	1130 (25.8)	1040 (25.9)	90 (23.9)
High school or equivalent	1010 (23.0)	915 (22.8)	95 (25.3)
College or above	2240 (51.1)	2050 (51.1)	190 (50.5)
Other	5 (0.1)	4 (0.1)	1 (0.3)
Marital status	0.133
Married	2275 (51.9)	2066 (51.5)	209 (55.6)
Unmarried	2110 (48.1)	1943 (48.5)	167 (44.4)
Family poverty ratio (%)	0.689
<1.3%	1461 (33.3)	1343 (33.5)	118 (31.4)
1.3–3.5%	1572 (35.8)	1435 (35.8)	137 (36.4)
≥3.5%	1352 (30.8)	1231 (30.7)	121 (32.2)
BMI (kg/m2)	0.025
Mean (SD)	28.98 (6.79)	28.92 (6.82)	29.67 (6.46)	0.040
<25 kg/m2	1319 (30.1)	1225 (30.6)	94 (25.0)
≥25 kg/m2	3066 (69.9)	2784 (69.4)	282 (75.0)
Hypertension	<0.001
Yes	1491 (34.0)	1326 (33.1)	165 (43.9)
No/Unknown	2894 (66.0)	2683 (66.9)	211 (56.1)
Diabetes	<0.001
Yes	503 (11.5)	438 (10.9)	65 (17.3)
No/Unknown	3882 (88.5)	3571 (89.1)	311 (82.7)
Vigorous recreational activities	0.002
Yes	941 (21.5)	884 (22.1)	57 (15.2)
No	3444 (78.5)	3125 (77.9)	319 (84.8)
Moderate recreational activities	0.146
Yes	1740 (39.7)	1604 (40.0)	136 (36.2)
No	2645 (60.3)	2405 (60.0)	240 (63.8)
Blood urea nitrogen, mg/dL	13.10 (5.78)	12.95 (5.63)	14.64 (7.06)	<0.001
Creatinine, mg/dL	0.89 (0.34)	0.89 (0.33)	0.94 (0.35)	0.003
Uric acid, mg/dL	5.50 (1.45)	5.48 (1.44)	5.73 (1.57)	0.001
eGFR	94.99 (23.6)	95.60 (23.57)	88.41 (23.62)	<0.001
Total PAHs, ng/L	147.41 (188.43)	144.96 (180.40)	173.50 (257.87)	0.005

Notes: For categorical variables, P values were analyzed by chi-square tests. For continuous variables, the t-test for slope was used in generalized linear models.

Abbreviations: BMI, body mass index; PAHs, polycyclic aromatic hydrocarbons; eGFR, estimated glomerular filtration rate.

Based on the total PAHs levels, we divided all participants into two groups: <79.66 ng/L (Low group) and ≥79.66 ng/L (High group). The clinical characteristics of participants in each group were shown in Table 2. Participants with higher total PAHs were concentrated in female, college or above, unmarried, family poverty ratio <1.3%, less recreational activities, kidney stones, higher creatinine, higher blood urea nitrogen, higher uric acid and lower eGFR. Moreover, similar results were observed in non-stone formers and stone formers ( and ).

Table 2

Characteristics of the Study Population by Categories of Total Polycyclic Aromatic Hydrocarbons (PAHs) Levels in NHANES 2007–2012 a

Characteristic	Total	Total PAHs Level		P valuea
	No. (%)	Low Group	High Group
Total participants	4385	2193 (50.0)	2192 (50.0)
Gender	<0.001
Male	2196 (50.1)	1194 (54.4)	1002 (45.7)
Female	2189 (49.9)	999 (45.6)	1190 (54.3)
Age	0.269
Mean (SD)	48.57 (17.70)	48.24 (17.17)	48.89 (18.21)	0.227
<50 years	2314 (52.8)	1139 (51.9)	1175 (53.6)
≥50 years	2071 (47.2)	1054 (48.1)	1017 (46.4)
Race	0.581
Non-Hispanic white	2009 (45.8)	1014 (46.2)	995 (45.4)
Non-Hispanic black	887 (20.2)	438 (20.0)	449 (20.5)
Mexican American	649 (14.8)	314 (14.3)	335 (15.3)
Other Hispanic	449 (10.2)	219 (10.0)	230 (10.5)
Other	391 (8.0)	208 (9.5)	183 (8.3)
Education	<0.001
Less than high school	1130 (25.8)	449 (20.5)	681 (31.1)
High school or equivalent	1010 (23.0)	466 (21.2)	544 (24.8)
College or above	2240 (51.1)	1275 (58.1)	965 (44.0)
Other	5 (0.1)	3 (0.1)	2 (0.1)
Marital status	<0.001
Married	2275 (51.9)	1225 (55.9)	1050 (47.9)
Unmarried	2110 (48.1)	968 (44.1)	1142 (52.1)
Family poverty ratio (%)	<0.001
<1.3%	1461 (33.3)	619 (28.2)	842 (38.4)
1.3–3.5%	1572 (35.8)	778 (35.5)	794 (36.2)
≥3.5%	1352 (30.8)	796 (36.3)	556 (25.4)
BMI (kg/m2)	0.091
Mean (SD)	28.98 (6.79)	28.86 (6.75)	29.09 (6.83)	0.288
<25 kg/m2	1319 (30.1)	634 (28.9)	685 (31.2)
≥25 kg/m2	3066 (69.9)	1559 (71.1)	1507 (68.8)
Hypertension	0.983
Yes	1491 (34.0)	746 (34.0)	745 (34.0)
No/Unknown	2894 (66.0)	1447 (66.0)	1447 (66.0)
Diabetes	0.117
Yes	503 (11.5)	235 (10.7)	268 (12.2)
No/Unknown	3882 (88.5)	1958 (89.3)	1924 (87.8)
Vigorous recreational activities	<0.001
Yes	941 (21.5)	563 (25.7)	378 (17.2)
No	3444 (78.5)	1630 (74.3)	1814 (82.8)
Moderate recreational activities	<0.001
Yes	1740 (39.7)	964 (44.0)	776 (35.4)
No	2645 (60.3)	1229 (56.0)	1416 (64.6)
Kidney stone	0.023
Yes	376 (8.6)	167 (7.6)	209 (9.5)
No	4009 (91.4)	2026 (92.4)	1983 (90.5)
Blood urea nitrogen, mg/dL	13.10 (5.78)	12.69 (5.70)	13.50 (5.84)	<0.001
Creatinine, mg/dL	0.89 (0.34)	5.42 (1.47)	5.58 (1.44)	<0.001
Uric acid, mg/dL	5.50 (1.45)	0.87 (0.36)	0.91 (0.31)	<0.001
eGFR	94.99 (23.6)	96.59 (23.34)	93.38 (23.87)	<0.001

Notes: aFor categorical variables, P values were analyzed by chi-square tests. For continuous variables, the t-test for slope was used in generalized linear models.

Abbreviations: BMI, body mass index; PAHs, polycyclic aromatic hydrocarbons; eGFR, estimated glomerular filtration rate.

Profiles of Urinary PAHs

Spearman’s rank correlation coefficients of PAHs, blood urea nitrogen, creatinine, uric acid and eGFR were shown in . The distribution levels of PAHs in this study were shown in . We found that the concentrations of 1-hydroxynaphthalene and 2-hydroxynaphthalene were higher among the nine PAHs ().

Urinary PAHs and Kidney Stone

Logistic regression analysis was performed to further identify risk factors associated with the prevalence of kidney stones. In the extended model, after adjusting for gender, age, race, education levels, marital status, family poverty, body mass index, hypertension, diabetes, recreational activities, blood urea nitrogen, uric acid, creatinine and eGFR, we found a significant positive correlation between total PAHs, 2-hydroxyphenanthrene, 1-hydroxyphenanthrene, 9-hydroxyfluorene and the risk of kidney stones, and the adjusted ORs for the prevalence of kidney stones increased significantly with the increase of total PAHS [OR (95% CI), 1.32 (1.06–1.64), p = 0.013], 2-hydroxynaphthalene [OR (95% CI), 1.37 (1.10–1.71), P = 0.005], 1-hydroxyphenanthrene [OR (95% CI), 1.24 (1.00–1.54), P=0.046] and 9-hydroxyfluorene [OR (95% CI), 1.36 (1.09–1.70), P=0.007] (Table 3, Figure 2). The nonlinear dose-response risk curve showed that the risk of kidney stones increased with the increase of total PAHs, 2-hydroxyphenanthrene, 1-hydroxyphenanthrene and 9-hydroxyfluorene. (Table 4, Figure 3).

Table 3

Adjusted Odds Ratios for Associations Between the Polycyclic Aromatic Hydrocarbons and the Presence of Kidney Stone in NHANES 2007–2012 a

Polycyclic Aromatic Hydrocarbons	Basic Model		Core Model		Extended Model
	aOR (95% CI)	P	aOR (95% CI)	P	aOR (95% CI)	P
Total PAHs
Low group	1.00	1.00	1.00
High group	1.35 (1.09–1.68)	0.006	1.38 (1.11–1.72)	0.004	1.32 (1.06–1.64)	0.013
1-hydroxynaphthalene
Low group	1.00	1.00	1.00
High group	1.13 (0.91–1.40)	0.272	1.12 (0.90–1.39)	0.317	1.11 (0.89–1.38)	0.353
2-hydroxynaphthalene
Low group	1.00	1.00	1.00
High group	1.33 (1.07–1.65)	0.009	1.41 (1.13–1.75)	0.003	1.37 (1.10–1.71)	0.005
3-hydroxyfluorene
Low group	1.00	1.00	1.00
High group	1.18 (0.95–1.46)	0.131	1.19 (0.95–1.48)	0.123	1.18 (0.94–1.47)	0.147
2-hydroxyfluorene
Low group	1.00	1.00	1.00
High group	1.25 (1.01–1.55)	0.039	1.27 (1.02–1.58)	0.033	1.25 (1.00–1.55)	0.051
3-hydroxyphenanthrene
Low group	1.00	1.00	1.00
High group	1.14 (0.92–1.41)	0.224	1.09 (0.88–1.36)	0.426	1.09 (0.87–1.35)	0.463
1-hydroxyphenanthrene
Low group	1.00	1.00	1.00
High group	1.33 (1.07–1.65)	0.009	1.27 (1.02–1.58)	0.034	1.24 (1.00–1.54)	0.046
2-hydroxyphenanthrene
Low group	1.00	1.00	1.00
High group	1.32 (1.07–1.64)	0.011	1.17 (0.94–1.46)	0.172	1.16 (0.93–1.44)	0.200
1-hydroxypyrene
Low group	1.00	1.00	1.00
High group	1.21 (0.97–1.50)	0.089	1.16 (0.93–1.45)	0.190	1.15 (0.92–1.44)	0.212
9-hydroxyfluorene
Low group	1.00	1.00	1.00
High group	1.41 (1.14–1.75)	0.002	1.38 (1.11–1.72)	0.004	1.36 (1.09–1.70)	0.007

Notes: p value <0.05 are shown in bold. aAdjusted covariates: Basic model: gender, age, and body mass index; Core model: basic model plus race, education levels, marital status and family poverty; Extended model: core model plus hypertension, diabetes, vigorous recreational activities, moderate recreational activities, blood urea nitrogen, uric acid, creatinine and estimated glomerular filtration rate (eGFR).

Abbreviations: CI, confidence interval; aOR, adjusted odds ratio.

Figure 2

Adjusted odds ratios for associations between the PAHs and the presence of kidney stone in NHANES 2007–2012. (A) Total PAHs; (B) 1-hydroxynaphthalene; (C) 2-hydroxynaphthalene; (D) 3-hydroxyfluorene; (E) 2-hydroxyfluorene; (F) 3-hydroxyphenanthrene; (G) 1-hydroxyphenanthrene; (H) 2-hydroxyphenanthrene; (I) 1-hydroxypyrene; (J) 9-hydroxyfluorene.

Table 4

Weighted Odds Ratios and 95% Confidence Intervals of Kidney Stones by Levels of PAHs, NHANES 2007–2012 a

Total PAHs, ng/L	200	400	600	800	1000	1200
Overall	1.10 (0.87–1.40)	1.26 (0.90–1.76)	1.43 (1.03–2.00)	1.63 (1.11–2.39)	1.86 (1.16–2.98)	2.12 (1.18–3.80)
1-hydroxynaphthalene, ng/L	100	200	400	600	800	1000
Overall	1.31 (0.99–1.73)	1.49 (1.05–2.11)	1.68 (1.18–2.38)	1.89 (1.23–2.89)	2.12 (1.24–3.64)	2.38 (1.22–4.65)
2-hydroxynaphthalene, ng/L	50	100	200	300	400	450
Overall	1.04 (1.00–1.08)	1.07 (0.83–1.37)	1.19 (0.86–1.65)	1.37 (0.94–2.01)	1.58 (0.91–2.74)	1.70 (0.89–3.25)
3-hydroxyfluorene, ng/L	2.5	5	10	15	20	25
Overall	1.26 (1.07–1.49)	1.55 (1.15–2.10)	1.67 (1.19–2.35)	1.53 (1.04–2.24)	1.39 (0.81–2.37)	1.27 (0.61–2.63)
2-hydroxyfluorene, ng/L	5	10	20	30	40	50
Overall	1.20 (1.04–1.37)	1.45 (1.08–1.96)	1.50 (1.07–2.09)	1.38 (0.92–2.06)	1.26 (0.71–2.24)	1.16 (0.53–2.53)
3-hydroxyphenanthrene, ng/L	2	4	6	8	10	12
Overall	1.23 (0.92–1.64)	1.33 (0.96–1.83)	1.39 (0.98–1.96)	1.45 (0.96–2.19)	1.52 (0.93–2.49)	1.59 (0.88–2.86)
1-hydroxyphenanthrene, ng/L	3	6	9	12	15	18
Overall	1.20 (0.91–1.56)	1.38 (1.02–1.88)	1.57 (1.12–2.19)	1.78 (1.20–2.63)	2.02 (1.26–3.24)	2.28 (1.30–4.01)
2-hydroxyphenanthrene, ng/L	1	2	4	6	8	10
Overall	1.14 (0.99–1.32)	1.34 (0.99–1.83)	1.41 (1.02–1.95)	1.47 (0.96–2.24)	1.53 (0.88–2.67)	1.60 (0.78–3.27)
1-hydroxypyrene, ng/L	2.5	5	10	15	20	25
Overall	1.53 (1.20–1.97)	1.74 (1.25–2.41)	1.51 (1.01–2.28)	1.32 (0.71–2.44)	1.15 (0.49–2.72)	1.00 (0.32–3.10)
9-hydroxyfluorene, ng/L	10	20	30	40	50	60
Overall	1.31 (0.97–1.77)	1.37 (0.99–1.89)	1.39 (1.00–1.95)	1.42 (0.98–2.06)	1.45 (0.95–2.21)	1.48 (0.92–2.38)

Notes: Adjusted covariates: gender, age, race, education levels, marital status, family poverty, body mass index, hypertension, diabetes, vigorous recreational activities, moderate recreational activities, blood urea nitrogen, uric acid, creatinine and estimated glomerular filtration rate (eGFR).

Figure 3

Relative risk for kidney stones based on PAHs level. The solid black lines represent aORs according to restricted cubic splines for PAHs level. The shaded areas represent upper and lower 95% CIs. Adjustment factors are as same as which presented in extended model of Table 3. (A) Total PAHs; (B) 1-hydroxynaphthalene; (C) 2-hydroxynaphthalene; (D) 3-hydroxyfluorene; (E) 2-hydroxyfluorene; (F) 3-hydroxyphenanthrene; (G) 1-hydroxyphenanthrene; (H) 2-hydroxyphenanthrene; (I) 1-hydroxypyrene; (J) 9-hydroxyfluorene.

Discussion

Kidney stones are a very common and recurring disease. It is considered to be caused by the abnormal accumulation of crystals in the kidneys. Among the various types of renal stones, those composed of Calcium oxalate (CaOx) are the commonest.1 Although there are many treatment methods for stones, the incidence and recurrence rate of kidney stones are still high, which exerts serious effects on patients’ lives.19 Therefore, exploring the risk factors of kidney stones is of great significance for reducing the social and economic burden and improving patients’ postoperative recovery.

PAHs are hydrocarbons containing more than two benzene rings. They are widely distributed in our lives through the incomplete combustion of fossil fuels such as coal, petroleum and the release of tar, coal, coke, pitch, and cigarettes.20 Many studies have reported that PAHs may cause damage to the respiratory system, digestive system, circulatory system, and urinary system.21–23 Mallin et al found that PAHs exposure would increase the occurrence risk of bladder cancer.24 Polycyclic aromatic hydrocarbon metabolites may lead to insulin resistance and b-cell dysfunction.25 Participants with higher total PAHs were concentrated in female, college or above, unmarried, family poverty ratio <1.3%, less recreational activities, kidney stones, higher creatinine, higher blood urea nitrogen, higher uric acid and lower eGFR. Li et al reported that PAHs exposure may be associated with albuminuria.26 PAHs may cause kidney damage by exacerbating kidney dysfunction. Polycyclic aromatic hydrocarbon metabolites can induce oxidative stress, even cause DNA damage in cells, and promote cardiovascular and cerebrovascular diseases.27 PAHs are absorbed by humans mainly through the respiratory tract, digestive tract, skin, and some will be excreted in urine or feces after being metabolized.28,29 Some of them further covalently bind with DNA or proteins, thus destroying the structure and function of proteins and nucleic acids, breaking the oxidation-antioxidant balance in the body.30 Only a few unmetabolized PAH prototypes will be eliminated from the body through urine or bile and feces.31 Studying the relationship between PAH metabolites and related diseases can more effectively explain the pathogenic mechanism.

At present, more and more studies on the occurrence of kidney stones are related to oxidative stress and inflammation.32 Oxidative stress is defined as the imbalance between oxidants and antioxidants. Reactive oxygen species (ROS) can mediate oxidative stress in mitochondria. At the same time, reactive oxygen species participate in CaOx crystals to promote renal tubular epithelial cell damage and apoptosis. Apoptosis and cell membrane structure changes can cause crystal adhesion and finally initiate a cascade of kidney stones.33,34 Huang’s research showed that lipid peroxidation was related to renal tubular epithelial cell damage, which was related to free radicals produced by patients with CaOx stones.35 PAHs can produce a large amount of ROS in the metabolic process of the human body. When it exceeds the body’s compensation, the excessive ROS can attack the biological macromolecules, causing oxidative stress, inflammatory response and immune response. Vijaya et al found that PAHs activate the oxidative stress mechanism after exposure, and the active molecules generated in the reaction can activate multiple signaling pathways in cells closely related to the immune cell function.36 Our study provides potential results for studying the relationship between PAHs exposure and the prevalence of kidney stones. The metabolism of PAHs may increase the incidence of kidney stones by increasing oxidative stress.

Previous studies have also examined the association between PAHs and celiac disease and kidney stones in adults. By collecting information from 5,560 adult participants in NHANES 2011–2012, Shiue found that after adjusting for urinary creatinine, age, sex, BMI, household income to poverty ratio, serum cotinine, alcohol consumption habits, education level, and physical activity level, urinary 2-hydroxyfluorene, 3-hydroxyfluorene, 1-hydroxyphenanthrene, 1-hydroxypyrene, and 2-hydroxynaphthalene were associated with kidney stones.37 In our study, we expanded the years of data collection and adjusted for sex, age, race, education level, marital status, household income to poverty ratio, BMI, hypertension, diabetes, physical activity level (vigorous and moderate), blood urea nitrogen, uric acid, creatinine, and eGFR. We found that total PAHs, 2-hydroxynaphthalene, 1-hydroxyphenanthrene, 2-hydroxyphenanthrene, and 9-hydroxy fluorene were significantly and positively correlated with the risk of kidney stones. Compared with the low group, the high group showed significant associations in total PAHs [OR (95% CI), 1.32 (1.06–1.64)], 2-hydroxynaphthalene [OR (95% CI), 1.37 (1.10–1.71)], 1-hydroxyphenanthrene [OR (95% CI), 1.24 (1.00–1.54)] and 9-hydroxyfluorene [OR (95% CI), 1.36 (1.09–1.70)]. In addition, we constructed dose-response curves between each PAHs and kidney stones to more visually demonstrate the risk relationship between each PAHs and kidney stones. Furthermore, we found that stone formers had higher blood urea nitrogen, creatinine, uric acid, total PAHs and lower eGFR than non-stone formers, possibly because PAHs may cause kidney damage by exacerbating kidney dysfunction, which in turn affects the excretion of PAHs.37,38

Using data from a large sample of NHANES in the United States, we can explore the relationship between PAHs and the prevalence of kidney stones while adjusting potential variables. Later, we can further discover its toxicity, pathogenic mechanism, and metabolic pathways, which will provide new methods for preventing and treating related diseases. There were several limitations of this study. First, due to the cross-sectional nature of the study, we could not fully determine the relationship between PAH exposure and kidney stones, which requires further prospective studies to verify its accuracy. Second, there is no information about the type of kidney stones. In addition, data on kidney stones were based on self-reporting by participants rather than professional diagnosis by physicians.

Conclusions

Our study shows that High levels of total PAHs were positively associated with the risk of kidney stones in the US population. Longitudinal studies are still needed to determine the exact relationship between PAHs and kidney stones.