Joint exposure to outdoor ambient air pollutants and incident chronic kidney disease: A prospective cohort study with 90,032 older adults

Objectives: There is paucity of studies to investigate the association between combined and long-term exposure to air pollution and the risk of incident chronic kidney disease (CKD) in older adults.

Methods: A prospective cohort of 90,032 older adults who did not have CKD at baseline were followed up from January 1, 2017, to December 31, 2019. Various pollutant data, including particulate matter with diameters ≤ 2.5 mm (PM_2.5), ≤ 10 mm (PM₁₀), nitrogen dioxide (NO₂), sulfur dioxide (SO₂), Ozone (O₃), and carbon monoxide (CO), from all monitoring stations in Binhai New Area, Tianjin were considered in calculating the mean exposure concentration of each pollutant over 2 years. By summing each pollutant concentration weighted by the regression coefficients, we developed an air pollution score that assesses the combined exposure of these air pollutants. Due to the strong correlation between air pollutants, Principal Component Analysis (PCA) score was also developed. The association between air pollutants and incident CKD in the elderly was analyzed.

Results: A total of 90,032 subjects participated in this study with a median follow-up of 545 days. Among them, 22,336 (24.8%) developed CKD. The HR (95% CI) for air pollution score and incidence of CKD was 1.062 (1.060-1.063) and p <0.001 after adjusting for all confounders. The adjusted HRs for the quartile subgroups of combined air pollution score were: Q2: 1.064 (1.013–1.117); Q3: 1.141 (1.088–1.198); and Q4: 3.623 (3.482–3.770), respectively (p for trend <0.001). The adjusted HRs for the quartile subgroups of air quality index (AQI) were: Q2: 1.035 (0.985–1.086); Q3: 1.145 (1.091–1.201); and Q4: 3.603 (3.463–3.748), respectively (p for trend <0.001). When the risk score was over 86.9, it significantly rose in a steep curve. The subgroup analysis showed that male, younger or exercise were more likely to develop CKD.

Conclusion: Combined air pollution score, AQI, and PCA score were associated with an increased risk of CKD in an exposure-response relationship. Our current results might also provide evidence for developing environmental protection policies.

GRAPHICAL ABSTRACT

Graphical Abstract. The combined air pollutants score for incident chronic kidney disease.

Introduction

Air pollution is a serious health problem worldwide and has become one of the major environmental problems in China (1). Chronic kidney disease (CKD) is a long-term chronic decline in kidney function and structural kidney damage, negatively impacting the quality of life as it progresses (2, 3). The Global Burden of Disease Study has revealed that CKD was the 16^th leading cause of life expectancy loss in 2017, with a 33.7% increase compared to 2007 (4). Furthermore, 17–20% of the worldwide CKD burden could be attributed to PM_2.5 and was more likely in low- and middle-income countries (5). However, research on the association between long-term exposure to air pollution and the risk of incident CKD in developing countries is still limited. In China, air pollution has become an increasingly prominent problem that cannot be ignored and needs to be addressed as soon as possible.

Several studies have shown that air pollution particles such as particulate matter with diameters ≤ 2.5 mm (PM_2.5) (6–17), ≤ 2.5 mm (PM₁₀) (11, 12), nitrogen dioxide (NO₂) (10–12, 14, 16), sulfur dioxide (SO₂) (11, 12, 16), Ozone (O₃) (12), carbon monoxide (CO) (11, 12), and air quality complex index (AQCI) (18) are associated with an increased risk of kidney outcomes. However, these studies have focused on individual air pollutants. Some studies also lacked variables such as the Air Quality Index (AQI), O₃, and NO₂. Therefore, they could not observe combined health effects (7). As individuals are often exposed to a mixture of multiple pollutants, although there is a high correlation between the pollutants, the synergism among pollutants are greater than the sum of the effects of individual pollutants (19). For effective management of multi-pollutant air quality, epidemiological research will be crucial to evaluating pollutant synergism. However, cohort studies focusing on the relationship between combined pollutants and CKD risk are still limited. Moreover, some studies lacked data on dietary, smoking, and drinking habits, which also substantially impact the morbidity of patients (18). Some studies were cross-sectional and, therefore, presented weak causal associations (13). The relationship between variables is more easily observed in large sample cohort studies and less convincing in studies conducted with small samples (12). Additionally, the elderly are more susceptible to decreased kidney function due to air pollutants. However, some studies have not included the elderly (18).

Herein, we used a large cohort population of older people from Tianjin and combined data from air pollutant monitoring sites. Then, we evaluated the association between pollutants and incident CKD, as well as the association between PCA score and air pollution score and CKD. Our current results provided evidence for the environmental risk factors for CKD in older adults based on a large samples and might facilitate the generation of public health policies.

Methods

Study design and participants

The Binhai New Area is located in the eastern coastal region of Tianjin. Data for this study were obtained from the Tianjin Chronic Kidney Disease Cohort, including follow-up information from different primary communities in the Binhai New Area, Tianjin, China. About 300,000 people undergo physical examinations each year in this area. Data were recorded in the Tianjin Community Health Service Information System from 2013 to 2019 and were available to downloaded.

Inclusion criteria: 1. Permanent residents living in the Binhai New Area, Tianjin, for at least 6 months. 2. Able to received regular follow-up visits. Exclusion criteria: 1. Age < 60 years. 2. Missing sociodemographic and clinical data. 3. History of cancer. 4. Combined hematuria or urinary tract infection. 5. Positive urine protein and estimated glomerular filtration rate (eGFR) < 60 mL/min/1.73 m² before January 1, 2017. 6. Self-reported CKD or other kidney diseases before the study. 7. Did not undergo annual physical examinations from January 1, 2017, to December 31, 2019. A total of 89,503 subjects were excluded. The details of the study process are presented in Figure 1.

FIGURE 1

Figure 1. Flow chart of the study.

A total of 90,032 participants were finally included in the analysis. The study protocol was approved by the Ethics Committee of the Chu Hsien-I Memorial Hospital of Tianjin Medical University and was registered in the Chinese Clinical Trial Registry (ChiCTR1900023701).

Definition of outcome and stages of chronic kidney disease

Due to the lack of direct diagnostic data for CKD, we evaluated the eGFR using Chronic Kidney Disease Epidemiology Collaboration (CKD-EPI) equation. We used the eGFR as a direct indicator of renal function and proteinuria as a sensitive indicator of renal function, according to the Kidney Disease: Improving Global Outcomes (KDIGO) 2020 Guidelines (20). The outcomes were defined as follows: 1. positive urine protein: daily excretion rate of ≥ 150 mg/24 h, and other confounding factors, such as infection were excluded; 2. eGFR < 60 mL/min/1.73 m². The survival time was defined from baseline to the event occurrence or the end of the investigation. According to international guidelines (20), stage G1: eGFR ≥ 90 mL/min/1.73 m². Stage G2: 60–89 mL/min/1.73 m². Stage G3a: 45–59 mL/min/1.73 m². Stage G3b: 30–44 mL/min/1.73 m². Stage G4: 15–29 mL/min/1.73 m². Stage G5: < 15 mL/min/1.73 m².

Measurement of air pollution variables

We used data from three ground-based environmental testing sites located in the Binhai New Area, Tianjin between 2014 and 2019, where the surveyed participants had lived for at least 1 year. Seven major air quality indicators, including the AQI, PM_2.5, PM₁₀, NO₂, SO₂, O₃, and CO, were obtained from the national city air quality real-time release platform (http://113.108.142.147:20035/emcpublish/) and the homepage of Wang Xiaolei (https://quotsoft.net/air/). Urban air quality monitoring stations obtain data on the pollutants present in the air through fixed and continuous sampling. Furthermore, the AQI was used to quantitatively describe the air quality condition. The pollutants included in the AQI are SO₂, NO₂, PM₁₀, PM_2.5, O₃, CO. According to the survey design, we calculated the average daily exposure of patient for 2 years before the onset of CKD or cessation of follow-up.

Definition of joint air pollution score and PCA score

The coefficients for each pollutant were calculated using a multivariable Cox proportional risk regression model. Confounders were adjusted in the composite pollutant model, and joint pollutant scores were calculated to analyze the relationship between joint air pollution exposure and the risk of incident CKD in the elderly. Similar to previous studies (21–24), we used the β coefficients from the final COX model to create the following formula: Air pollution score = 0.44 (β1) ^* PM_2.5 + 0.281 (β2) ^* PM₁₀ + 0.865 (β3) ^* NO₂ + 1.095 (β4) ^* SO₂ + 17.33 (β5) ^* CO + (-0.521) (β6) ^* O_3. The air pollution score ranged between 80.85 and 110.77. Then, we distributed the patients into four groups based on the quartiles of the scores.

In order to address the problem of high covariance between air pollutants and their potential interaction, Principal Component Analysis (PCA) was used to analyze these pollutants. The formula of the PCA socre is: PCA score = (0.9603 ^* f1 + 0.0227 ^* f2 + 0.0099 ^* f3)/0.9929. The details are shown in Supplementary Tables 1, 2.

Measurements of covariates

Baseline data included demographic characteristics (age, gender), behavioral health habits (smoking, alcohol, exercise frequency, diet), history of diseases (diabetes, hypertension), physical examination [body mass index (BMI), waist circumference (WC), systolic and diastolic blood pressure (SBP and DBP)], and laboratory tests [white blood cell (WBC), platelet (PLT), fasting blood glucose (FBG), serum creatinine, eGFR, aspartate transaminase (AST), alanine transaminase (ALT), total cholesterol (TC), and triglyceride (TG)]. The eGFR was calculated using serum creatinine according to the modified MDRD formula. Urine protein was determined using the immunoturbidimetric method. Briefly, urine was collected after cleaning the urethral orifice and vulva before urine retention while avoiding contamination by mixing menstrual blood, leukorrhea, semen, or feces. The first 200 mL of urine was collected early in the morning, and the specimen was delivered within half an hour, with the maximum time not exceeding 2 h. All study personnel were highly trained, and strict quality control procedures were carried out.

Statistical methods

Indicators with continuous normal distribution are expressed as means (standard deviations - SDs), and a t-test was used to compare their differences. Indicators with continuous non-normal distribution are expressed as medians (25th−75th), and the Mann-Whitney U-test was used to assess their differences. Categorical indicators were expressed as counts (%), and the χ² test was used to compare them. The HR and 95% confidence intervals (CIs) were estimated using the Cox proportional risk regression model in which several potential confounders were adjusted. Potential confounders were screened by univariate COX analysis. For sensitivity analysis, the HRs (95% CI) was calculated when CKD stages 3–5 were considered as outcome, which indicates a significant decrease in eGFR and < 60 mL/min/1.73 m². Restricted cubic splines provide a powerful way to represent non-linear relationships for continuous independent variables by dividing the observation range of X variables by 3–5 knot points. Thus, we used restricted spline regression to estimate the dose-response relationships between individual pollutants and these scores (AQI, PCA score and the air pollution score) and CKD, which has the flexibility to display non-linear relationships (25). We calculated the inflection point of the cubic splines. Additionally, we conducted a stratified analysis according to gender (male, female), age (60–64, 65–69, 70–74, 75–79, ≥80 years), CKD stages (G1–G5), BMI (<25, ≥25 kg/m²), smoking (yes or no), drinking (yes or no), exercise (yes or no), diabetes (yes or no) and hypertension (yes or no) at baseline, and their effect relationships were examined based on interactions effect. The discrimination of these scores was assessed using receiver operating characteristic curve (ROC) and area under the curve (AUC). Calibration was assessed with calibration curves. All analyses were performed using R software (v 4.1.0) and results are expressed by HR and 95% CI. The p-values for all tests were two-sided, and a p <0.05 was considered statistically significant.

Results

Baseline characteristics of the study population

During a median follow-up of 558 d, a total of 22,336 patients developed CKD among the 90,032 elderly without CKD at baseline. Of these, 4,530 patients had CKD at or above stage 3. The baseline characteristics of the participants according to incident CKD are shown in Table 1. The patients who developed CKD were older (69.07 vs. 68.23 years), females (64.4 vs. 47.3%), and had a higher BMI (25.49 vs. 25.03 kg/m²) compared to patients that did not developed CKD. Patients with diabetes and hypertension at baseline were more likely to develop CKD. In the CKD group, the means (SDs) of PM_2.5, PM₁₀, AQI, NO₂, SO₂, O₃, and CO were 58.2 (4.20), 95.0 (6.29), 89.3 (3.95), 47.9 (1.95), 13.8 (1.68), 63.8 (3.25) μg/m³, and 1.20 (0.10) mg/m³, respectively. In the control group, the means (SDs) were 54.8 (3.12), 90.0 (5.15), 86.1 (3.02), 46.4 (1.50), 12.5 (1.19), 66.1 (2.09) μg/m³, and 1.12 (0.07) mg/m³, respectively. Moreover, the proportion of CKD increased with aging, and the percentage of people with stage G2 was the highest (48,836, 54.24%) (Figure 2). Additionally, most patients had only a mild decrease in kidney function.

TABLE 1

Table 1. Baseline characteristics of participants in this study.

FIGURE 2

Figure 2. Subgroups of patients based on age and CKD stages. (A) Incidence and proportion of CKD patients in different age groups. (B) The number and proportion of patients in different CKD stages. Stage G1: eGFR ≥90 ml/min/1.73 m². Stage G2: 60–89 ml/min/1.73 m². Stage G3a: 45–59 ml/min/1.73 m². Stage G3b: 30–44 ml/min/1.73 m².

The relationship between air pollutants and risk of incident CKD

The Spearman correlation analysis suggested a strong relationship between these pollutants (Figure 3). The HR and 95% CI were calculated separately for each pollutant quartile group and are presented in Table 2. As the quartiles increased, the risk of developing CKD significantly increased with the Q1 group as the control. The univariable Cox regression results showed that PM_2.5, PM₁₀, SO₂, and CO were associated with an increased risk of incident CKD. Additionally, the correlation between each pollutant and the risk of incident CKD remained significant in model 2 and 3 after adjusting for confounding factors (Table 2).

FIGURE 3

Figure 3. Spearman correlation between air pollutants.

TABLE 2

Table 2. HR and 95% CI for air pollution concentrations with the risk of incident CKD in the study.

Association between the AQI, PCA score and combined score and risk of incident CKD

Based on the results of the multivariable analysis and PCA analysis, we developed air pollution score and PCA score. Then, we evaluated the relationship between the scores and CKD risk. For each one-point increase in the air pollution score, the adjusted HR corresponding to CKD risk was 1.062 (1.060–1.063) (Table 3). The quartile grouping of air pollution scores revealed that the HRs (95% CI) for the Q2, Q3, and Q4 groups were 1.064 (1.013–1.117), 1.141 (1.088–1.198), and 3.623 (3.482–3.770), respectively (p for trend <0.001). In addition, when the outcome was CKD stage 3–5, the HRs (95% CI) for the Q2, Q3, and Q4 groups of were 1.035 (0.932–1.149), 1.117 (1.010–1.236), 3.464 (3.191–3.760), respectively (p for trend <0.001) (Table 3).

TABLE 3

Table 3. Hazard ratios and 95% confidence interval for the risk of incident CKD in the study.

Next, we investigated the relationship between the AQI and CKD (Table 3). For every one increase in the AQI, the HR (95% CI) was 1.170 (1.166–1.174) (p <0.001). The quartile grouping of AQI revealed that the HRs (95% CI) for the Q2, Q3, and Q4 groups were 1.035 (0.985–1.086), 1.145 (1.091–1.201), and 3.603 (3.463–3.748), respectively (p for trend <0.001). Then we investigated the relationship between the PCA score and CKD (Table 3). For every one increase, the HR (95% CI) was 1.264 (1.258–1.270) (p <0.001). The quartile grouping of PCA score showed that the adjusted HRs (95% CI) for the Q2, Q3, and Q4 groups were 1.064 (1.013–1.117), 1.141 (1.088–1.198), and 3.623 (3.482–3.770), respectively (p for trend <0.001). Furthermore, both the PCA score and AQI showed significant trends of HR when using CKD stage 3–5 as outcomes (Table 3).

After adjusting for potential confounders, we calculated the non-linear relationship between these scores and CKD risk using restricted spline regression. No significant change in the HR was detected when the score was ≤ 86.9 (Figure 4). Thus, the inflection point for the restricted spline regression was 86.9. The risk of CKD would rapidly increase when the score was over 86.9. Similarly, the inflection point for the PCA score and AQI were 1.5 and 85.9, respectively. Altogether, these results indicated an exposure-response relationship between these scores and the risk of incident CKD. The AUC of AQI, PCA score and combined score were 0.713, 0.715, and 0.715, respectively. Their calibration curves are also similar (Supplementary Figure 1).

FIGURE 4

Figure 4. The relationships between air pollution indicators and incident CKD by restriction spline regression. (A) PM_2.5. (B) PM₁₀. (C) SO₂. (D) O₃. (E) NO₂. (F) CO. (G) Air Quality Index (AQI). (H) combined score of all pollutants. (I). PCA score.

Stratified analyses

Furthermore, we conducted stratified analyses for age, gender, BMI, hypertension, diabetes, behavioral habits (smoking, drinking, exercise). For air pollution score and PCA score, The subgroup analysis and interaction effect analyses revealed that male, exercise or younger had a higher risk of incident CKD at the same exposure levels (p <0.05). For AQI, The subgroup analysis indicated that male and younger had a higher risk of incident CKD (p <0.05) (Figure 5).

FIGURE 5

Figure 5. Subgroup analysis and HR (95% CI) for air pollution score, air quality score and PCA score and CKD. All subgroup analyses adjusted for age, gender, BMI, smoking, drinking, exercise, history of diabetes, history of hypertension,WC, FBG, SBP, DBP, WBC, HGB, PLT, TC, TG, ALT, eGFR at baseline.

Discussion

Herein, we conducted a large cohort study in northern China and found that long-term exposure to PM_2.5, PM₁₀, SO₂, and CO were associated with an increased risk of CKD after adjusting for all potential confounding factors. We developed an air pollution score for assessing CKD risk by calculating Cox regression coefficients after adjusting for all potential confounding factors. In addition, due to the strong correlation between air pollutants, PCA score was used to analyze these pollutants. The results revealed no significant change in CKD risk when these scores was below the inflection point. On the other hand, when the risk score was higher, the CKD risk would significantly rise in a steep curve. These results indicated a dose-response relationship between AQI, PCA score and air pollution score and CKD risk, independent of traditional risk factors. Altogether, these results suggested that the prevalence of environmental pollution can significantly effect the kidney health of the elderly population. Our current results provided evidence for further developing public environmental protection policies and encouraging people to make efforts to protect the environment.

Several previous cohort studies have shown that PM_2.5 (6–10, 13, 15–17, 26), PM₁₀ (10, 26–28), SO₂ (15, 26, 27), NO₂(14, 26, 27), AQI (18), CO (26, 27) are associated with an increased risk of CKD, consistent with our current findings. Besides, a meta-analysis showed that PM_2.5, PM₁₀, NO₂, SO₂ and CO are associated with CKD (28). Some studies have shown that PM_2.5 and CO are associated with death in CKD patients (11, 17). Li et al. have shown that short- or medium-term NO₂ exposure are associated with kidney damage (12).

According to the subgroup analyses, we showed that the incidence is higher in women, however, male were more vulnerable to combined air pollution exposure than female. A meta-analysis shows increased risk of progression in men compared with women (29). The EQUAL study also demonstrate faster declines in renal function in men compared with women, even after adjustment for multiple groups of mediators, which may explain why male received more effects (30). Further researches may be needed to explain this pathogenesis. In the current study, patient who are obesity (BMI ≥ 25 kg/m²) are more susceptible to the effects of air pollution. The association between increased air pollutant concentrations and CKD prevalence is stronger in the visceral obesity group than normal group (27). According to our findings, people who do not drink alcohol or exercise are more susceptible to the effects of air pollution. As reported by Cui Guo et al., in comparison to no exercise, exposure to equivalent PM2.5 levels is associated with a greater risk of death (31), which means that people who exercise are more susceptible than those who do not exercise. Additionally, younger patients had a higher risk of developing CKD due to air pollution (p <0.001), a potential explanation for this is an imbalance in the aging population, with people of higher ages dying from other events. There are no reports on whether alcohol consumption or aging impacts the relationship between air pollution and CKD. It is only reasonable to conclude that non-drinkers and younger are more likely to be affected by air pollution than drinkers.

Although some studies based on satellite spatio-temporal models can accurately predict ground-level PM_2.5, the measurements of other pollutants, such as SO₂, NO₂, PM₁₀, AQI, and CO, are missing and cannot be directly determined (8, 9, 14, 27). In this present study, we used measurement data from surface monitoring sites, which comprehends more reliable data than satellite-based assessment. The current literature on SO₂, NO₂, PM₁₀, AQI, and CO exposure in CKD studies remains limited. Thus, we explored all variables described above. Based on these significant indicators, we established a joint risk score and PCA score, which was more reliable than a single indicator. Besides, this was a longitudinal cohort study, starting with non-CKD patients who were observed for nearly 2 years. In this period, some patients developed CKD, with stronger causal associations compared to those cross-sectional studies (13, 14).

The findings on the relationship between PM_2.5 exposure and CKD have been inconsistent, with a small number of studies not supporting the conclusion that the exposure increases CKD risk. In a retrospective cohort study (32), the results showed that long-term exposure to air pollution was unlikely to increase CKD risk. Moreover, after unit conversion, the average concentrations of PM10 and CO were 61.7 and 775 μg/m³, respectively. However, in our current study, the mean concentrations of PM₁₀ and CO was 95.03 and 1,200 μg/m³, respectively. Thus, the concentration of PM₁₀ was 1.54 times, and CO was 1.55 times higher than in the study of Hwang et al. Additionally, in the study by Hwang et al., 71.9% of the population was under 60 years old, while all patients were over 60 years old in our current research. Hence, we considered that these different results were derived from the different average age of the surveyed population and the concentrations of pollutants. Accordingly, people over 60 years old can have CKD due to greater sensitivity to air pollution. These contradictory results emphasized the necessity to increase the studies in areas with high air pollution.

To the best of our knowledge, this was the first study to assess the risk of CKD in older adults from air pollution exposure using score models based on a large cohort. Additionally, the results were still robust after adjusting for confounders. Moreover, it is possible that certain components play a major role, we are exposed to the air as a whole and cannot be exposed to one specific component alone, while our calculations showed that most component was significant, so it makes sense to assess the overall pollution score. Moreover, this study suggested an inflection score of 86.9 that can be used as a control target. Nevertheless, we believe that this inflection point only applies to Han-Chinese people over 60 years.

Our current study also has some limitations. First, it lacked measurement data for assessing indoor air pollution, older adults are exposed to indoor air for more extended period. There is a significant difference in the concentration of pollutants between indoor and outdoor air, indoor air quality is influenced by the use of air purifiers, hoods, home structure, fuel, ventilation (33). Thus, the effects of these critical factors might be overlooked, next researches should investigate the relationship between indoor air and older adults' health. Second, due to the extensive range of individual activities, we did not precisely estimate the exposure level of each individual but only averaged the pollution values of the Binhai New Area to roughly estimate the exposure level of each patient. Additionally, the follow-up period of this study was short (2 years), and the lifetime-risk still needs to be studied. In the elderly cohort, bias may exist due to unequal numbers of age groups, with the older groups have risk of death. In our study, participants were required to participate at two follow-up visits. Essentially, this means that the dead individual has been removed. Next, we will consider the effects of death and use a competitive risk model to assess potential risk factors for CKD. Moreover, More confounding factors might need to be adjusted, such as consumption of alcohol, amount and duration of smoking, and type and duration of exercise. Finally, our current study showed the incidence rate was 24.8%, which maybe overestimated. However, when we use the CKD stage 3–5 as the outcome, a similar trend has been observed. CKD is common in older people and its prevalence increases in parallel with age (34). The prevalence of CKD in China and the United States were 34.6% and 31.5–32.9% (60–89 years) (35, 36). The population in the current study was based on adults over 60 years in China, they might be more susceptible to air pollution exposure. Hence, our results do not suitable for people under 60 years and non-Chinese.

Conclusion

In summary, we demonstrated that both single and combined exposure to air pollutants was associated with an increased risk of CKD in the elderly. The air pollution score, PCA score and AQI were associated with risk of incident CKD in a dose-response relationship. This study would provide evidence for the development of environmental protection policies and emphasis the importance of persistent efforts to control air pollution.

Data availability statement

The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.

Ethics statement

The studies involving human participants were reviewed and approved by Chu Hsien-I Memorial Hospital of Tianjin Medical University. Written informed consent for participation was not required for this study in accordance with the national legislation and the institutional requirements.

Author contributions

HYL and XS: methodology, software, formal analysis, visualization, and writing-original draft. PY: supervision, conceptualization, resources, writing-review and editing, project administration, and funding acquisition. XJ: data curation, validation, and investigation. SJZ, XJL, and YL: validation and investigation. PFB, JMC, FH, YZ, CL, and HL: investigation. ZC: formal analysis and writing-review and editing. All authors read and approved the final manuscript.

Funding

This work was supported by the financial support from Tianjin Key Medical Discipline (Specialty) Construct Project (No. TJYXZDXK-032A), the Science and technology talent project of Tianjin Health Commission (No. RC20175), Scientific Research Program of Tianjin Education Commission (No. 2020KJ187), and Tianjin Science and Technology Plan Project Public Health Science and Technology Major Special Project (No. 21ZXGWSY00100). The funder was not for profit.

Acknowledgments

The authors thank the staff and participants of the study for their indispensable contributions.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher's note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpubh.2022.992353/full#supplementary-material

Abbreviations

PM_2.5, particulate matter with an aerodynamic diameter lower than 2.5 μm. NO₂, nitrogen dioxide. PM₁₀, particulate matter with an aerodynamic diameter lower than 10 μm. AQI, air quality index. NO₂, nitrogen oxide. SO₂, sulfur dioxide. O₃, Ozone. CO, carbon monoxide. FBG, Fasting blood glucose. HGB, hemoglobin. PLT, platelet. WBC, white blood cell. ALT, alanine transaminase. AST, aspartate transaminase. TC, total cholesterol. TG, triglyceride. BMI, body mass index. WC, waist circumference. SBP, systolic blood pressure. DBP, diastolic blood pressure. eGFR, estimated glomerular filtration rate. CKD, chronic kidney disease.

References

1. Liu C, Chen R, Sera F, Vicedo-Cabrera AM, Guo Y, Tong S, et al. Ambient particulate air pollution and daily mortality in 652 cities. N Engl J Med. (2019) 381:705–15. doi: 10.1056/NEJMoa1817364

PubMed Abstract | CrossRef Full Text | Google Scholar

2. Al-Aly Z, Bowe B. Air pollution and kidney disease. Clin J Am Soc Nephrol. (2020) 15:301–3. doi: 10.2215/CJN.16031219

PubMed Abstract | CrossRef Full Text | Google Scholar

3. Webster AC, Nagler EV, Morton RL, Masson P. Chronic kidney disease. Lancet. (2017) 389:1238–52. doi: 10.1016/S0140-6736(16)32064-5

PubMed Abstract | CrossRef Full Text | Google Scholar

4. GBD 2017 Causes of Death Collaborators. Global, regional, and national age-sex-specific mortality for 282 causes of death in 195 countries and territories, 1980-2017: a systematic analysis for the Global Burden of Disease Study 2017. Lancet. (2018) 392:1736–88. doi: 10.1016/S0140-6736(18)32203-7

PubMed Abstract | CrossRef Full Text | Google Scholar

5. Bowe B, Xie Y, Li T, Yan Y, Xian H, Al-Aly Z. Estimates of the 2016 global burden of kidney disease attributable to ambient fine particulate matter air pollution. BMJ Open. (2019) 9:e022450. doi: 10.1136/bmjopen-2018-022450

PubMed Abstract | CrossRef Full Text | Google Scholar

6. Blum MF, Surapaneni A, Stewart JD, Liao D, Yanosky JD, Whitsel EA, et al. Particulate matter and albuminuria, glomerular filtration rate, and incident CKD. Clin J Am Soc Nephrol. (2020) 15:311–9. doi: 10.2215/CJN.08350719

PubMed Abstract | CrossRef Full Text | Google Scholar

7. Bo Y, Brook JR, Lin C, Chang LY, Guo C, Zeng Y, et al. Reduced ambient PM(25) was associated with a decreased risk of chronic kidney disease: a longitudinal cohort study. Environ Sci Technol. (2021) 55:6876–83. doi: 10.1021/acs.est.1c00552

PubMed Abstract | CrossRef Full Text | Google Scholar

8. Bowe B, Xie Y, Li T, Yan Y, Xian H, Al-Aly Z. Particulate matter air pollution and the risk of incident CKD and progression to ESRD. J Am Soc Nephrol. (2018) 29:218–30. doi: 10.1681/ASN.2017030253

PubMed Abstract | CrossRef Full Text | Google Scholar

9. Chan TC, Zhang Z, Lin BC, Lin C, Deng HB, Chuang YC, et al. Long-term exposure to ambient fine particulate matter and chronic kidney disease: a cohort study. Environ Health Perspect. (2018) 126:107002. doi: 10.1289/EHP3304

PubMed Abstract | CrossRef Full Text | Google Scholar

10. Chen SY, Chu DC, Lee JH, Yang YR, Chan CC. Traffic-related air pollution associated with chronic kidney disease among elderly residents in Taipei City. Environ Pollut. (2018) 234:838–45. doi: 10.1016/j.envpol.2017.11.084

PubMed Abstract | CrossRef Full Text | Google Scholar

11. Jung J, Park JY, Kim YC, Lee H, Kim E, Kim YS, et al. Effects of air pollution on mortality of patients with chronic kidney disease: a large observational cohort study. Sci Total Environ. (2021) 786:147471. doi: 10.1016/j.scitotenv.2021.147471

PubMed Abstract | CrossRef Full Text | Google Scholar

12. Li A, Mei Y, Zhao M, Xu J, Li R, Zhao J, et al. Associations between air pollutant exposure and renal function: a prospective study of older adults without chronic kidney disease. Environ Pollut. (2021) 277:116750. doi: 10.1016/j.envpol.2021.116750

PubMed Abstract | CrossRef Full Text | Google Scholar

13. Li G, Huang J, Wang J, Zhao M, Liu Y, Guo X, et al. Long-term exposure to ambient PM(25) and increased risk of CKD prevalence in China. J Am Soc Nephrol. (2021) 32:448–58. doi: 10.1681/ASN.2020040517

PubMed Abstract | CrossRef Full Text | Google Scholar

14. Liang Z, Wang W, Wang Y, Ma L, Liang C, Li P, et al. Urbanization, ambient air pollution, and prevalence of chronic kidney disease: a nationwide cross-sectional study. Environ Int. (2021) 156:106752. doi: 10.1016/j.envint.2021.106752

PubMed Abstract | CrossRef Full Text | Google Scholar

15. Lin YT, Lo YC, Chiang HY, Jung CR, Wang CM, Chan TC, et al. Particulate air pollution and progression to kidney failure with replacement therapy: an advanced CKD Registry-Based Cohort Study in Taiwan. Am J Kidney Dis. (2020) 76:645–57. doi: 10.1053/j.ajkd.2020.02.447

PubMed Abstract | CrossRef Full Text | Google Scholar

16. Lin SY, Ju SW, Lin CL, Hsu WH, Lin CC, Ting IW, et al. Air pollutants and subsequent risk of chronic kidney disease and end-stage renal disease: a population-based cohort study. Environ Pollut. (2020) 261:114154. doi: 10.1016/j.envpol.2020.114154

PubMed Abstract | CrossRef Full Text | Google Scholar

17. Ran J, Sun S, Han L, Zhao S, Chen D, Guo F, et al. Fine particulate matter and cause-specific mortality in the Hong Kong elder patients with chronic kidney disease. Chemosphere. (2020) 247:125913. doi: 10.1016/j.chemosphere.2020.125913

PubMed Abstract | CrossRef Full Text | Google Scholar

18. Wang HH, Zhang SC, Wang J, Chen X, Yin H, Huang DY. Combined toxicity of outdoor air pollution on kidney function among adult women in Mianyang City, southwest China. Chemosphere. (2020) 238:124603. doi: 10.1016/j.chemosphere.2019.124603

PubMed Abstract | CrossRef Full Text | Google Scholar

19. Mauderly JL, Samet JM. Is there evidence for synergy among air pollutants in causing health effects? Environ Health Perspect. (2009) 117:1–6. doi: 10.1289/ehp.11654

PubMed Abstract | CrossRef Full Text | Google Scholar

20. Kidney Disease: Improving Global Outcomes Diabetes Work Group. KDIGO 2020 clinical practice guideline for diabetes management in chronic kidney disease. Kidney Int. (2020) 98:S1–115. doi: 10.1016/j.kint.2020.06.019

PubMed Abstract | CrossRef Full Text | Google Scholar

21. Exalto LG, Exalto LG, Biessels GJ, Karter AJ, Huang ES, Katon WJ, et al. Risk score for prediction of 10 year dementia risk in individuals with type 2 diabetes: a cohort study. Lancet Diabetes Endocrinol. (2013) 1:183–90. doi: 10.1016/S2213-8587(13)70048-2

PubMed Abstract | CrossRef Full Text | Google Scholar

22. Kivipelto M, Ngandu T, Laatikainen T, Winblad B, Soininen H, Tuomilehto J. Risk score for the prediction of dementia risk in 20 years among middle aged people: a longitudinal, population-based study. Lancet Neurol. (2006) 5:735–41. doi: 10.1016/S1474-4422(06)70537-3

PubMed Abstract | CrossRef Full Text | Google Scholar

23. Lindstrom J, Tuomilehto J. The diabetes risk score: a practical tool to predict type 2 diabetes risk. Diabetes Care. (2003) 26:725–31. doi: 10.2337/diacare.26.3.725

PubMed Abstract | CrossRef Full Text | Google Scholar

24. Wang M, Zhou T, Song Y, Li X, Ma H, Hu Y, et al. Joint exposure to various ambient air pollutants and incident heart failure: a prospective analysis in UK Biobank. Eur Heart J. (2021) 42:1582–91. doi: 10.1093/eurheartj/ehaa1031

PubMed Abstract | CrossRef Full Text | Google Scholar

25. Marrie RA, Dawson NV, Garland A. Quantile regression and restricted cubic splines are useful for exploring relationships between continuous variables. J Clin Epidemiol. (2009) 62:511–7 e1. doi: 10.1016/j.jclinepi.2008.05.015

PubMed Abstract | CrossRef Full Text | Google Scholar

26. Yang YR, Chen YM, Chen SY, Chan CC. Associations between long-term particulate matter exposure and adult renal function in the Taipei metropolis. Environ Health Perspect. (2017) 125:602–7. doi: 10.1289/EHP302

PubMed Abstract | CrossRef Full Text | Google Scholar

27. Jeong SM, Park JH, Kim HJ, Kwon H, Hwang SE. Effects of abdominal obesity on the association between air pollution and kidney function. Int J Obes. (2020) 44:1568–76. doi: 10.1038/s41366-020-0540-8

PubMed Abstract | CrossRef Full Text | Google Scholar

28. Ye JJ, Wang SS, Fang Y, Zhang XJ, Hu CY. Ambient air pollution exposure and risk of chronic kidney disease: a systematic review of the literature and meta-analysis. Environ Res. (2021) 195:110867. doi: 10.1016/j.envres.2021.110867

PubMed Abstract | CrossRef Full Text | Google Scholar

29. Minutolo R, Gabbai FB, Chiodini P, Provenzano M, Borrelli S, Garofalo C, et al. Sex differences in the progression of CKD among older patients: pooled analysis of 4 cohort studies. Am J Kidney Dis. (2020) 75:30–8. doi: 10.1053/j.ajkd.2019.05.019

PubMed Abstract | CrossRef Full Text | Google Scholar

30. Nicholas CC, Friedo WD, Marie E, Fergus JC, Claudia T, Maurizio P, et al. Renal function decline in older men and women with advanced chronic kidney disease-results from the EQUAL study. Nephrol Dial Transplant. (2021) 27:1656–63. doi: 10.1093/ndt/gfaa095

PubMed Abstract | CrossRef Full Text | Google Scholar

31. Cui G, Tsung Y, LyYun C, Changqing L, Hsiao TY, Yacong B, et al. (2021). Effects of air pollution and habitual exercise on the risk of death: a longitudinal cohort study. CMAJ. 193, E1240–9. doi: 10.1503/cmaj.211282

PubMed Abstract | CrossRef Full Text | Google Scholar

32. Hwang SY, Jeong S, Choi S, Kim DH, Kim SR, Lee G, et al. Association of air pollutants with incident chronic kidney disease in a nationally representative cohort of Korean Adults. Int J Environ Res Public Health. (2021) 18:3775. doi: 10.3390/ijerph18073775

PubMed Abstract | CrossRef Full Text | Google Scholar

33. Bi J, Wallace LA, Sarnat JA, Liu Y. Characterizing outdoor infiltration and indoor contribution of PM25 with citizen-based low-cost monitoring data. Environ Pollut. (2021) 276:116763. doi: 10.1016/j.envpol.2021.116763

PubMed Abstract | CrossRef Full Text | Google Scholar

34. Tonelli M, Riella M. Chronic kidney disease and the aging population. Am J Nephrol. (2014) 39:248–51. doi: 10.1159/000359957

PubMed Abstract | CrossRef Full Text | Google Scholar

35. Murphy D, McCulloch CE, Lin F, Banerjee T, Bragg-Gresham JL, Eberhardt MS, et al. Trends in prevalence of chronic kidney disease in the United States. Ann Intern Med. (2016) 165:473–81. doi: 10.7326/M16-0273

PubMed Abstract | CrossRef Full Text | Google Scholar

36. Wang F, He K, Wang J, Zhao MH Li Y, Zhang L, et al. Prevalence and risk factors for CKD: a comparison between the adult populations in China and the United States. Kidney Int Rep. (2018) 3:1135–43. doi: 10.1016/j.ekir.2018.05.011

PubMed Abstract | CrossRef Full Text | Google Scholar