Skip to main content

Air pollution, residential greenness, and metabolic dysfunction biomarkers: analyses in the Chinese Longitudinal Healthy Longevity Survey

Abstract

Background

We hypothesize higher air pollution and fewer greenness exposures jointly contribute to metabolic syndrome (MetS), as mechanisms on cardiometabolic mortality.

Methods

We studied the samples in the Chinese Longitudinal Healthy Longevity Survey. We included 1755 participants in 2012, among which 1073 were followed up in 2014 and 561 in 2017. We used cross-sectional analysis for baseline data and the generalized estimating equations (GEE) model in a longitudinal analysis. We examined the independent and interactive effects of fine particulate matter (PM2.5) and Normalized Difference Vegetation Index (NDVI) on MetS. Adjustment covariates included biomarker measurement year, baseline age, sex, ethnicity, education, marriage, residence, exercise, smoking, alcohol drinking, and GDP per capita.

Results

At baseline, the average age of participants was 85.6 (SD: 12.2; range: 65–112). Greenness was slightly higher in rural areas than urban areas (NDVI mean: 0.496 vs. 0.444; range: 0.151–0.698 vs. 0.133–0.644). Ambient air pollution was similar between rural and urban areas (PM2.5 mean: 49.0 vs. 49.1; range: 16.2–65.3 vs. 18.3–64.2). Both the cross-sectional and longitudinal analysis showed positive associations of PM2.5 with prevalent abdominal obesity (AO) and MetS, and a negative association of NDVI with prevalent AO. In the longitudinal data, the odds ratio (OR, 95% confidence interval-CI) of PM2.5 (per 10 μg/m3 increase) were 1.19 (1.12, 1.27), 1.16 (1.08, 1.24), and 1.14 (1.07, 1.21) for AO, MetS and reduced high-density lipoprotein cholesterol (HDL-C), respectively. NDVI (per 0.1 unit increase) was associated with lower AO prevalence [OR (95% CI): 0.79 (0.71, 0.88)], but not significantly associated with MetS [OR (95% CI): 0.93 (0.84, 1.04)]. PM2.5 and NDVI had a statistically significant interaction on AO prevalence (pinteraction: 0.025). The association between PM2.5 and MetS, AO, elevated fasting glucose and reduced HDL-C were only significant in rural areas, not in urban areas. The association between NDVI and AO was only significant in areas with low PM2.5, not under high PM2.5.

Conclusions

We found air pollution and greenness had independent and interactive effect on MetS components, which may ultimately manifest in pre-mature mortality. These study findings call for green space planning in urban areas and air pollution mitigation in rural areas.

Peer Review reports

Background

Metabolic syndrome (MetS) is a risk factor for morbidity and mortality. Specifically, it is a group of pathologic conditions that precede non-communicable diseases, including cardiovascular disease (CVD) and diabetes [1]. It has become a global problem with the increasing prevalence in both developed and developing countries [2]. There are plenty of amenable causes of MetS. An increasing number of studies have been focusing on environmental determinants.

Fine particulate matter (PM2.5) is an independent risk factor for mortality in many locations and exposure levels [3]. PM2.5 has been implicated in causing systemic inflammation and altered metabolism of lipids and glucose [4,5,6]. At the same time, living in areas with higher greenness is associated with a reduced risk of mortality and cardiovascular disease [7]. However, there was no established evidence on the association between PM2.5 and MetS according to current controversial findings in various countries [8, 9]. A limited number of research findings in China were inconsistent [10, 11]. Compared to air pollution, much less attention has been paid to greenness and MetS worldwide, especially for the older adults aged 80 or older, and there was also little agreement [12,13,14]. Some prior findings showed combined or synergistic effects of PM2.5 and greenness on mortality [15, 16]. No studies looked at their interaction on MetS based on our knowledge.

The relationship between air pollution and residential greenness can be complex and need additional analyses for generalizability in different climates, income levels, and places with varying population density. A recent study based on a Canadian cohort of 2.4 million individuals found adjustment of greenness attenuated the effect of PM2.5. The effect of air pollution on cardiovascular mortality was the largest in places with the least greenness. Studies that do not account for greenness may overstate the harmful effect of air pollution on mortality [15]. In a seven metropolitan cities study in South Korea, the effect of PM10 was higher in areas of lower greenness for cardiovascular-related mortality, but not for non-accidental mortality and respiratory-related mortality [17]. A cohort study spanning 22 provinces in China of elderly individuals found that people living in urban areas experienced higher health benefits of greenness. People living in rural regions were more likely to be harmed by air pollution [16]. Not all studies found a significant interaction between greenness and air pollution. An Israel-based study found the incorporation of greenness into the PM2.5 model did not improve the cardiovascular disease predictions for stroke and myocardial infarction, although air pollution and greenness had strong independent effects on these outcomes [18]. As for MetS, KORA F4/FF4 cohort in Germany and Whitehall II study in the UK found the association between greenness and MetS was reversed and became positive after adjusting for PM2.5 in the model. In contrast, 33 Communities Chinese Health Study (33CCHS) in China found this association was only partly attenuated after adjusting for air pollution [12,13,14].

Large uncertainty still exists about the pattern and mechanisms of greenness and air pollution impact on MetS. With the rapid urbanization and population aging in developing countries, including China, the role of these environmental determinants is yet to be determined. Using a cohort of older adults in eight regions in China, we aim to (1) estimate the prevalence of MetS and its components based on measured biomarkers, (2) determine the independent effects of PM2.5 and greenness on metabolic syndrome biomarkers, (3) assess the interactive effect of PM2.5 and greenness, and (4) to assess effect modification by age, gender, and urban versus rural regions. These analyses are anticipated to generate insights that can improve our limited understanding of whether and how the two important environmental factors related to urbanization affect metabolic syndrome, a health problem with increasing prevalence in rapidly developing parts of the world.

Methods

Study population

We used data from the sub-cohort of the Chinese Longitudinal Healthy Longevity Survey: Healthy Ageing and Biomarkers Cohort Study (HABCS). The study collected blood samples for biomarker examinations during 2008 to 2017 in eight places designated as longevity areas (Laizhou City of Shandong Province, Xiayi County of Henan Province, Zhongxiang City of Hubei Province, Mayang County of Hunan Province, Yongfu County of Guangxi Autonomous Area, Sanshui District of Guangdong Province, Chengmai County of Hainan Province and Rudong County of Jiangsu Province). The published cohort profile described the study design and sample method [19]. The waist circumference was measured since 2012. We set the study baseline at 2012 and excluded 85 participants aged younger than 65, 286 participants with missing biomarker value, 91 participants with missing NDVI or PM2.5 value, and 222 participants with missing covariates value (Fig. S1). We finally included 1755 participants at baseline. During 2012–2017, 1115 participants were followed up at least twice, and 519 participants were followed up three times.

Air pollution and residential greenness measurements

Ground-level PM2.5 concentrations were estimated by the Atmospheric Composition Analysis Group. They combined aerosol optical depth retrievals from the National Aeronautics and Space Administration’s Moderate Resolution Imaging Spectroradiometer, Multi-angle Imaging SpectroRadiometer, and Sea-viewing Wide field-of-view Sensor satellite instruments; vertical profiles derived from the GEOS-Chem chemical transport model; and calibration to ground-based observations of PM2.5 using geographically weighted regression [20]. The resultant PM2.5 concentration estimates were highly consistent (R2 = 0.81) with out-of-sample cross-validated PM2.5 concentrations from monitors. We matched the annual average PM2.5 concentrations in a 1 km × 1 km grid to each participant’s residence [21].

We calculated Normalized Difference Vegetation Index (NDVI) with a 500-m radius around each participant’s residence to quantify greenness exposure. We used satellite images from the Moderate-Resolution Imaging Spectro-Radiometer (MODIS) in the National Aeronautics and Space Administration’s Terra Satellite. The NDVI calculation formula is near-infrared radiation minus visible radiation divided by near-infrared radiation plus visible radiation, ranging from − 1.0 to 1.0, with larger values indicating higher vegetative density levels. There are two NDVI values for January, April, July, and October between 2008 and 2014 in our database to reflect the seasonal variation of greenness. We linked NDVI imagery to the longitude and latitude of each residential address and calculated greenness in 500 m radii.

We matched time-varying annual PM2.5 and NDVI of 2008–2014 to the data. We calculated the average value of one-year, three-year, and five-year exposure time windows as long-term cumulative exposures measurements. We used the same exposure results as the 2014 wave for the 2017 wave since we lacked the environmental exposure data from 2014 to 2017.

Biomarker measurements

The participants provided the blood sample at the same time as the interview time in 2012, 2014, and 2017. The medical technician tested blood plasma biomarkers included fasting glucose, glycated serum protein (GSP), total cholesterol (TC), triglyceride (TG), and high-density lipoprotein cholesterol (HDL-C) using an Automatic Biochemistry Analyzer (Hitachi 7180, Japan) with commercially available diagnostic kits (Roche Diagnostic, Mannheim, Germany) at Capital Medical University in Beijing. Low-density lipoprotein cholesterol (LDL-C) was calculated using the formula of Friedewald et al.: LDL-C = TC-(HDL-C)-TG/5 [22].

Trained medical staff performed anthropometric measurements for the participants, including waist circumference, and two blood pressure measurements with at least a one-minute interval between them. We used the mean value of the two blood pressure measurements.

Definition of metabolic syndrome (MetS) and components

We defined the MetS using the Adult Treatment Panel III of the National Cholesterol Education Program (ATP III) guidelines, modified in accordance with the waist circumference cutoff points proposed by World Health Organization (WHO) for Asian populations (modified ATP III). It was defined as the presence of at least three of the following criteria: elevated fasting glucose (fasting glucose≥100 mg/dL), abdominal obesity (AO: Waist circumference ≥ 90 cm for males and ≥ 80 cm for females), hypertension (SBP ≥ 130/DBP ≥ 85 mmHg), hypertriglyceridemia (TG ≥ 150 mg/dL), and reduced HDL-C (HDLC< 40 mg/dL for males and < 50 mg/dL for females) [23, 24]. We also did sensitivity analysis for the MetS defined by the Joint Interim Societies [25].

Baseline covariates

We categorized the ethnicity as Han Chinese or ethnic minorities. We used years in schools as a measure of literacy level. We classified marital status into two categories: currently married and living with the spouse, or not married (widowed/separated/divorced/never married/married but not living with the spouse). We classified city and town as “Urban”, and village as “Rural.” We firstly divided the regular exercise, smoking, and alcohol drinking status into three categories: “Current,” “Former,” and “Never”. For example, participants were asked, “do you do exercise regularly at present (planned exercise like walking, playing balls, running and so on)?” and/or “did you do exercise regularly in the past?”. We defined the regular exercise status as “Current” for participants who answered “Yes” to the first question, “Former” for who answered “No” to the first question and “Yes” to the second question, and “Never” for who answered “No” to both two questions. Then we further quantified the current smoker based on the number of times smoke (or smoked) per day: < 20 times/day and ≥ 20 times/day. We also quantified the current alcohol drinker based on the kind of alcohol and how much they drank per day. The unit of alcohol was a Chinese unit of weight called ‘Liang’ [50 g (g)]. The level of alcohol consumption was calculated as drinks of alcohol per day, based on the beverage type and amount, assuming the following alcohol content by volume (v/v) typically seen in China: strong liquor 53%, weak liquor 38%, grape wine 12%, rice wine 15%, and beer 4% [26]. A standard drink was equal to 14.0 g of pure alcohol according to the criterion of the Center for Disease Control and Prevention in the USA, and moderate drinking is up to 1 drink per day for women and up to 2 drinks per day for men according to Dietary Guidelines for Americans 2015–2020. Therefore, we defined those who drank equal or less than 14 g pure alcohol per day for the female or 28 g per day for the male as light drinkers, otherwise heavy drinkers. We collected Gross Domestic Product (GDP) per capita by county/district from the local statistical yearbook.

Statistical analysis

We described univariate statistics of our exposure, outcome variables, and covariates in eight areas. We built the multivariate logistic regression model in the cross-sectional analysis to analyze the association between residential environment (residential greenness and ambient air pollution) and baseline MetS and each component. For the longitudinal analysis, we used generalized estimating equations (GEE) to assess the association between the repeatedly measured residential environment and the repeatedly measured metabolic biomarkers. For each biomarker: firstly, we built the single exposure model to regress only one environment factor on the biomarker; Second, we built the two-exposure model to regress both greenness and air pollution on the biomarker; Third, we added the product term of centered greenness and air pollution (NDVI×PM2.5) in the model to assess their interaction and one exposure’s association with the outcome under another exposure’s mean level. We adjusted for biomarker measurement year, baseline age, sex, ethnicity, education, marriage, residence, exercise, smoking, alcohol drinking, and GDP per capita in these models. Considering gender difference plays a vital role in the health of the old population, we further examined the greenness, air pollution, and gender three-way interaction by adding the term “NDVI×PM2.5 × Sex” in the model. We performed sensitivity analyses using environment exposure of different time windows (1 year or five-year average NDVI or PM2.5). Given the selection bias due to lost to follow-up, we also built models for those with at least one follow-up. We conducted stratified analyses based on age, sex, and residence to test the possible modification. We set the nominal significance level at 0.05. We used R 4.0.0 to run all the analyses.

Results

Population characteristics and environmental exposure level

We studied 1755 participants aged 65 to 112 years old, with a mean age of 85 (SD:12.2); 53.8% were female. Most were Han participants (92.3%), lived in rural areas (83.1%), never had regular exercise (81.9%), never smoked (75.4%), and never drank alcohol (77.9%). There were 370 (21.1%) participants who fit the criteria for MetS, 583 (33.2%) for abdominal obesity (AO), 307 (17.5%) for elevated fasting glucose, 1285 (73.2%) for hypertension, 157 (8.9%) for hypertriglyceridemia, and 679 (38.7%) for reduced HDL-C (Table 1). Those who were lost of follow-up were older, more likely to be female, living in areas with higher GDP, not currently married, and without formal education (Table S1).

Table 1 Baseline population characteristics

PM2.5 was not associated with NDVI (Pearson correlation coefficient: 0.0004; p > 0.05). The three-year NDVI (0.1 unit) of the rural area was slightly higher than the urban area (mean: 4.96 vs. 4.44; range: 1.51–6.98 vs. 1.33–6.44), and the mean of three-year PM2.5 (10 μg/m3) were almost the same in the rural and urban areas (mean: 4.90 vs. 4.91; range: 1.62–6.53 vs. 1.83–6.42) of our sample (Table 1). The mean of the three-year NDVI (0.1 unit) of the eight counties was 4.88 (SD: 0.94), ranging from 3.36 (0.81) in Sanshui to 5.37 (0.59) in Rudong. The mean of three-year PM2.5 (10 μg/m3) of the eight areas was 4.90 (SD: 1.53), ranging from 1.83 μg/m3 (SD: 0.03) in Chengmai to 6.42 μg/m3 (SD: 0.02) in Xiayi (Fig. 1, Table S2).

Fig. 1
figure 1

The NDVI and PM2.5 level in the eight sample districts. Note: We used “ggplot2” and “sf” packages in R 4.0.0 (URL https://www.R-project.org/) to draw the map

Environmental exposure and MetS

In both the cross-sectional and longitudinal analyses, higher PM2.5 was associated with higher odds of MetS [OR (95%CI): 1.17 (1.07, 1.28) and 1.16 (1.08, 1.24) respectively], and the association between NDVI and MetS tended to be negative but was not statistically significant [OR (95%CI): 0.94 (0.81, 1.09) and 0.93 (0.84, 1.04) respectively]. These associations did not change when adding both PM2.5 and NDVI in the model, and there was no significant interaction between PM2.5 and NDVI on MetS (Table 2 & Table S3).

Table 2 The association between the greenness and air pollution with the metabolic syndrome and the components (Binary outcome) in the longitudinal analysisa

Environmental exposure and MetS components

In both the cross-sectional and longitudinal analyses, higher PM2.5 was associated with higher odds of AO [OR (95%CI): 1.25 (1.16, 1.36) and 1.19 (1.12, 1.27) respectively], while higher NDVI was associated with lower odds of AO [OR (95% CI): 0.81 (0.71, 0.92) and 0.79 (0.71, 0.88) respectively] (Table 2 & Table S3). In addition, higher PM2.5 was associated with higher waist circumference [mean difference (95% CI): 1.12 (0.83, 1.40)] while higher NDVI was associated with lower waist circumference [mean difference (95% CI): − 1.21 (− 1.76, − 0.66)] (Table S4).

For the lipids, higher PM2.5 was only associated with higher odds of reduced HDL-C [OR (95%CI): 1.14 (1.07, 1.21)] in the longitudinal analyses. There were no significant association between PM2.5 and TG or hypertriglyceridemia, or between NDVI and TG, HDL-C, hypertriglyceridemia or reduced HDL-C. Besides, PM2.5 and NDVI were both negatively associated with TC and LDL-C (Table 2, Table S4). The association between PM2.5 and elevated fasting glucose were not statistically significant in either cross-sectional or longitudinal analyses [OR (95%CI): 1.08 (0.99, 1.19) and 1.06 (0.99, 1.13) respectively]. NDVI showed a negative association with the odds of elevated fasting glucose only in the cross-sectional analyses [OR (95%CI): 0.84 (0.72, 0.99)] (Table 2, Table S3). Both PM2.5 and NDVI were not associated with hypertension in either cross-sectional or longitudinal analyses. These results also persisted in the two-exposure model (Table 2, Table S3 and S4).

Sensitivity analyses

Using the one-year and five-year average exposure window, the above associations persisted except for that the positive association between one-year PM2.5 and odds of elevated fasting glucose became statistically significant (Table S5). Among those with at least one follow-up, the results did not change significantly either (Table S6). The findings based on the Joint Interim Societies definition of MetS were also similar (Table S7).

Possible effect modification

We found a significant interaction of PM2.5 and NDVI on AO (beta estimate of interaction term = − 0.088, P = 0.025) and waist circumference (beta estimate of interaction term = − 0.396, P = 0.031) (Table 2, Table S4). Higher PM2.5 was associated with a higher probability of AO, and the association for exposure beyond 30 μg/m3 became stronger with the increase of the greenness level. Higher NDVI was associated with a lower probability of AO and the association was stronger under relatively higher PM2.5 exposure (Fig. 2). For the three-way interaction of air pollution, greenness, and gender on metabolic biomarkers, we only found a significant three-way interaction on GSP. In areas with low NDVI, the association strength and direction of PM2.5 with GSP in the females were different from males, and applies in areas with high NDVI (Fig. S2).

Fig. 2
figure 2

The interaction model of PM2.5 and NDVI on abdominal obesity in the longitudinal analysis. Note: The figure was based on the logistic regression for abdominal obesity including the interaction term of PM2.5 and NDVI adjusting for biomarker measurement year, baseline age, sex, ethnicity, education, marriage, residence, exercise, smoking, alcohol drinking, and GDP per capita. Higher PM2.5 was associated with higher probability of AO, and the effect size decreased with the increase of the greenness level for exposure beyond 30 μg/m3. Higher NDVI was associated with lower probability of AO and the effect size was stronger under relatively higher PM2.5 exposure. We used R package "interactions" to draw the figure.

In the stratified analysis, the association between PM2.5 and AO was weaker in areas with high NDVI exposure than areas with low NDVI [OR (95%CI): 1.17 (1.08, 1.28) vs. 1.25 (1.13, 1.39)]. The association between NDVI and AO was only significant in areas with low PM2.5 [OR (95%CI): 0.61 (0.52, 0.73)]. PM2.5 shown a harmful association with MetS, AO, elevated fasting glucose, and reduced HDL-C only in rural areas [OR (95%CI): 1.18 (1.09, 1.28) for MetS, 1.22 (1.14, 1.30) for AO, 1.08 (1.01, 1.16) for elevated fasting glucose, and 1.15 (1.07, 1.23) for reduced HDL-C], not in urban areas. NDVI’s protective association with AO was a little stronger in urban areas than rural areas. The association between PM2.5 with MetS, AO, reduced HDL-C were stronger in the male than female, and the association between NDVI with AO were similar for males and females. The association between PM2.5 and MetS as well as its components were all more significant in the population aged younger than 80 compared to those aged 80 or older. NDVI was still not associated with MetS in the two different age groups, but had a stronger association with AO in those younger than 80 (Table 3).

Table 3 The association between the greenness and air pollution with the metabolic syndrome and the components (Binary outcome) in the longitudinal analysis stratified by PM2.5, NDVI, age, sex, and residencea

Discussion

We found air pollution could increase the risk of MetS, AO, and reduced HDL-C while residential greenness could decrease the risk of AO. We further identified an environment-environment interaction of air pollution-greenness on AO. The association strength for air pollution decreased along with the increase of greenness. The association for greenness was stronger under high-level air pollution exposure than that under low-level air pollution.

Two recent meta-analysis studies on air pollution and MetS showed inconsistent findings. One found PM2.5 (per 10 μg/m3 increase) was not significantly associated with MetS prevalence [OR (95% CI): 1.34 (0.96, 1.89), P = 0.09] or MetS incidence [Hazard ratio (HR): 2.78 (95% CI: 0.70, 11.02), P = 0.15] [8], while another one found annual PM2.5 (per 5 μg/m3 increase) was associated with 14% of MetS risk increase [Risk Ratio (RR): 1.14 (95% CI: 1.03, 1.25)] [9]. The included studies reported associations of different sizes in varied areas. Some studies were conducted in areas with a mean PM2.5 higher than 50 μg/m3. A study in northern rural China reported the adjusted OR of MetS for per 5 μg/m3 increment in PM2.5 was 1.42 (95% CI: 1.36, 1.49) [11], while another study only found borderline associations and reported the adjusted odds ratio of MetS per 10 μg/m3 increment in PM2.5 was 1.09 (95% CI: 1.00, 1.18) in northern urban China [10]. A Korean national cohort found PM2.5 level was significantly associated with a higher risk for developing MetS [HR (95% CI): 1.07 (1.03, 1.11)] [27]. Some studies were conducted in areas with a mean PM2.5 lower than 50 μg/m3. The study in Saudi Arabian population in Jeddah observed a significant association between a 10 μg/m3 increase in PM2.5 and increased risks for MetS [RR (95% CI): 1.12 (1.06, 1.19)] [28]. Another study in the highly urbanized German Ruhr Area reported the OR of per interquartile range (IQR = 1.5 μg/m3) PM2.5 was 1.04 (95% CI: 0.92, 1.17) for MetS prevalence and 1.21 (95% CI: 0.99, 1.48) for MetS incidence [29]. A 1-μg/m3 increase of PM2.5 was associated with a higher risk of developing MetS [HR (95% CI): 1.27 (1.06, 1.52)] in an US older men cohort [27]. We found PM2.5 was only significantly associated with MetS in rural areas [OR (95%CI) for 10 μg/m3 increment in PM2.5: 1.18 (1.09, 1.28)], and not in urban populations. More studies on air pollution-MetS risk association, especially in low−/middle-income countries, are warranted.

There are a few meta-analyses demonstrated the association between PM2.5 and MetS composition biomarker: long-term exposure of PM2.5 was associated with a higher level of BMI with the pooled β (95% CI) of 0.34 (0.30, 0.38) per 10 mg/m3 increment [30], higher type 2 diabetes incidence [HR (95% CI): 1.10 (1.04, 1.17) per 10 μg/m3 increment] [6], and higher hypertension prevalence [OR (95% CI):1.05 (1.01, 1.09)] [31]. A few studies found air pollutants only significantly associated with TC, not with HDL-C or TG [5]. A previous CLHLS study reported higher 3-year average exposure to PM2.5 was associated with higher fasting blood glucose [32]. In our research, we also found higher PM2.5 associated with AO, reduced HDL-C and elevated fasting glucose, which was robust among different age and sex groups. However, we only saw PM2.5 increased the risk for elevated fasting glucose in rural areas, and risk for hypertriglyceridemia in the population aged younger than 80. We found no significant association between PM2.5 and hypertension.

The negative association between greenness and MetS tended to be insignificant in the elderly based on previous studies, which congruent to our observation. KORA F4/FF4 cohort in German found a negative association between greenness and MetS in both cross-sectional and longitudinal analysis in German but both were insignificant [14]. The 33CCHS conducted in northern urban China found the adjusted OR of MetS per IQR increase in 500 m buffer NDVI of August was 0.81 (95% CI: 0.70, 0.93) for the total population aged 18–74 years, but the association disappeared in subgroup participants aged ≥65 [13]. Whitehall II study in the UK (aged 45–69 years at baseline) found a significant negative association [12]. We did not find a significant association of NDVI on MetS in any subgroup in urban or rural areas, for female or male, aged from 65 to 80 or older than 80.

For MetS composition biomarker, a recent meta-analysis showed higher NDVI was associated with lower odds of overweight/obesity [OR (95% CI): 0.88 (0.84, 0.91)], and most studies were from developed nations (88%) [33]. We also found NDVI associated with lower odds of AO. The possible pathway can be that green spaces could decrease sympathetic nervous system activation [34]. A study in urban northeastern China found higher greenness was consistently associated with lower TC, TG, LDL-C levels, higher HDL-C levels [35], and lower fasting glucose levels [36]. We also found greenness negatively associated with TC, LDL-C, but not associated with TG, HDL-C, or fasting glucose.

We found PM2.5 and NDVI were both associated with the metabolic biomarkers. The association varied in different age, sex, and residence categories. PM2.5 inhalation could cause pulmonary and systemic inflammation. According to the animal findings, rats that were exposed to Beijing’s highly polluted air experienced the following changes: perivascular and peribronchial inflammation in the lungs, increased tissue and systemic oxidative stress, dyslipidemia, and enhanced proinflammatory status of epididymal fat. TLR2/4-dependent inflammatory activation and lipid oxidation in the lung can spill over systemically, leading to metabolic dysfunction and weight gain [37]. The pathways linking greenness to health include physical activity (50 studies), air pollution (43 studies), social interaction/cohesion (27 studies), mental health/stress/well-being (17 studies), perceived greenness/use (16 studies), and physical health/biomarker (14 studies) and other factors according to the latest review of previous empirical studies [38]. Greenness may decrease the risk for obesity by promoting exercise. Greenness and air pollution may act in separate pathways since our two exposure models showed no major mediation effect according to the similar estimates of the single exposure and two-exposure models.

For the relationship between air pollution and greenness, a longitudinal study in China found a significant interaction between PM2.5 and NDVI on all-cause mortality, and individuals living in areas with more greenness appear to be affected more by air pollution, but it showed no monotonic trend [16]. An ecological study in Greece found a significant inverse interaction between PM2.5 and NDVI on cardiovascular mortality with the PM2.5 effects decreasing in areas with higher greenery, and they found no interaction on natural-cause mortality [39]. Previous studies have related both greenness and PM2.5 with metabolic syndrome and biomarkers. However, most studies only considered PM2.5 as a mediator of greenness. There has been no study reported on the interaction of air pollution and greenness on metabolic biomarkers. We reported NDVI had a significant interaction with PM2.5 on AO, but no interaction on metabolic syndrome.

Our study has several strengths. First, our cohort has a relatively older mean age than previous studies, and it has a large sample of centenarians which is rare in the world. Secondly, a limited number of studies focused on greenness and the multiple exposures of both air pollution and greenness. While individual studies on environmental predictors exist, ours is a novel approach to assessing the interaction of air pollution and greenness on metabolic syndrome biomarkers. Third, many previous studies were conducted in specific regions like rural or urban areas. We identified high-risk vulnerable older adults from different geographic regions of China. Fourth, we had repeat measurements of a variety of individual metabolic biomarkers. Fifth, we calculated the greenness and air pollution level at the individual residence level, and we tested different exposure time windows before the health outcome. We also surveyed a wide range of lifestyle and district factors to adjust for possible confounding.

There are several limitations to our study. The specific oldest-old population also limited the generalizability of our findings. Those who were lost to follow-up were older, with a possible selection bias. Thus, we did sensitivity analysis only for those with at least one follow-up, and the results persisted. We lacked the exposure data from 2015 to 2017 and used the same exposure as the 2014 wave for the 2017 wave. We found this should not affect our results much since the trend of PM2.5 across 2008–2014 was steady within each area. The sensitivity analysis showed no significant difference among one-year, three-year, and five-year exposure windows. There is also no extensive heterogeneity of PM2.5 measurement among participants within each area. This possible misclassification usually attenuates the association to null, which means the exposure of higher resolution may show a stronger association with the health outcomes. In addition, we have no indoor air pollution measurements or greenness accessibility data to account for the dynamic personal exposure, which limited the accuracy of the exposure measurement. For the outcome, we lack the metabolism medication information to better define the metabolic syndrome, which may cause underestimating MetS prevalence. We presented the real-world observational evidence, and there may be residual confounding like the diet. We conducted multiple comparisons without correction, for which we exercised caution by presenting confidence intervals and exact p-value.

Conclusions

Our findings contributed to the evidence of harmful association of PM2.5 and protective association of NDVI with specific MetS components in an oldest-old population, newly identified a significant interaction between PM2.5 and NDVI on AO, and demonstrated the difference between urban and rural areas. Other than the personal actionable lifestyle risk factors, it is also necessary to incorporate environmental determinants into metabolic diseases prevention. This study emphasized the importance of green space planning in urban areas and air pollution mitigation in rural areas to decrease the CVD burden contributed by MetS biomarkers for the policymakers. Further studies can examine if PM2.5 and NDVI only interact or if their effect can counteract each other and explore the underlying biology pathway.

Availability of data and materials

The CLHLS datasets are available upon request to the public from the Peking University Open Research Data on CLHLS (http://opendata.pku.edu.cn/dataverse/CHADS).

Abbreviations

PM2.5 :

Fine particulate matter

MetS:

Metabolic syndrome

CLHLS:

Chinese Longitudinal Healthy Longevity Survey

NDVI:

Normalized Difference Vegetation Index

GSP:

Glycated serum protein

TC:

Total cholesterol

TG:

Triglyceride

HDL-C:

High-density lipoprotein cholesterol

LDL-C:

Low-density lipoprotein cholesterol

AO:

Abdominal obesity

SBP:

Systolic blood pressure

DBP:

Diastolic blood pressure

GEE:

Generalized estimating equations

OR:

Odds ratio

CI:

Confidence interval

IQR:

Interquartile range

References

  1. Cornier MA, Dabelea D, Hernandez TL, Lindstrom RC, Steig AJ, Stob NR, et al. The metabolic syndrome. Endocr Rev. 2008;29:777–822.

    Article  CAS  Google Scholar 

  2. Saklayen MG. The global epidemic of the metabolic syndrome. Curr Hypertens Rep. 2018;20:12. https://doi.org/10.1007/s11906-018-0812-z.

    Article  PubMed  PubMed Central  Google Scholar 

  3. Orellano P, Reynoso J, Quaranta N, Bardach A, Ciapponi A. Short-term exposure to particulate matter (PM10 and PM2.5), nitrogen dioxide (NO2), and ozone (O3) and all-cause and cause-specific mortality: systematic review and meta-analysis. Environ Int. 2020;142:105876. https://doi.org/10.1016/j.envint.2020.105876.

    Article  CAS  PubMed  Google Scholar 

  4. Xing YF, Xu YH, Shi MH, Lian YX. The impact of PM2.5 on the human respiratory system. J Thorac Dis. 2016;8:E69–74.

    Article  Google Scholar 

  5. Gaio V, Roquette R, Dias CM, Nunes B. Ambient air pollution and lipid profile: Systematic review and meta-analysis. Environ Pollut. 2019;254:113036.

    Article  CAS  Google Scholar 

  6. Yang BY, Fan S, Thiering E, Seissler J, Nowak D, Dong GH, et al. Ambient air pollution and diabetes: a systematic review and meta-analysis. Environ Res. 2020;180:108817. https://doi.org/10.1016/j.envres.2019.108817.

    Article  CAS  PubMed  Google Scholar 

  7. Yuan Y, Huang F, Lin F, Zhu P, Zhu P. Green space exposure on mortality and cardiovascular outcomes in older adults: a systematic review and meta-analysis of observational studies. Aging Clin Exp Res. 2021;33:1783–97. https://doi.org/10.1007/s40520-020-01710-0.

    Article  PubMed  Google Scholar 

  8. Zang ST, Luan J, Li L, Wu QJ, Chang Q, Dai HX, et al. Air pollution and metabolic syndrome risk: evidence from nine observational studies. Environ Res. 2021;202:111546. https://doi.org/10.1016/j.envres.2021.111546.

    Article  CAS  PubMed  Google Scholar 

  9. Ning J, Zhang Y, Hu H, Hu W, Li L, Pang Y, et al. Association between ambient particulate matter exposure and metabolic syndrome risk: a systematic review and meta-analysis. Sci Total Environ. 2021;782:146855. https://doi.org/10.1016/j.scitotenv.2021.146855.

    Article  CAS  PubMed  Google Scholar 

  10. Yang BY, Qian ZM, Li S, Fan S, Chen G, Syberg KM, et al. Long-term exposure to ambient air pollution (including PM 1 ) and metabolic syndrome: The 33 Communities Chinese Health Study (33CCHS). Environ Res. 2018;164:204–11. https://doi.org/10.1016/j.envres.2018.02.029.

    Article  CAS  PubMed  Google Scholar 

  11. Hou J, Liu X, Tu R, Dong X, Zhai Z, Mao Z, et al. Long-term exposure to ambient air pollution attenuated the association of physical activity with metabolic syndrome in rural Chinese adults: a cross-sectional study. Environ Int. 2020;136:105459. https://doi.org/10.1016/j.envint.2020.105459.

    Article  CAS  PubMed  Google Scholar 

  12. de Keijzer C, Basagaña X, Tonne C, Valentín A, Alonso J, Antó JM, et al. Long-term exposure to greenspace and metabolic syndrome: a Whitehall II study. Environ Pollut. 2019;255:113231. https://doi.org/10.1016/j.envpol.2019.113231.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  13. Yang BY, Liu KK, Markevych I, Knibbs LD, Bloom MS, Dharmage SC, et al. Association between residential greenness and metabolic syndrome in Chinese adults. Environ Int. 2020;135:105388. https://doi.org/10.1016/j.envint.2019.105388.

    Article  PubMed  Google Scholar 

  14. Voss S, Schneider A, Huth C, Wolf K, Markevych I, Schwettmann L, et al. Long-term exposure to air pollution, road traffic noise, residential greenness, and prevalent and incident metabolic syndrome: results from the population-based KORA F4/FF4 cohort in Augsburg, Germany. Environ Int. 2021;147:106364. https://doi.org/10.1016/j.envint.2020.106364.

    Article  CAS  PubMed  Google Scholar 

  15. Crouse DL, Pinault L, Balram A, Brauer M, Burnett RT, Martin RV, et al. Complex relationships between greenness, air pollution, and mortality in a population-based Canadian cohort. Environ Int. 2019;128:292–300. https://doi.org/10.1016/j.envint.2019.04.047.

    Article  CAS  PubMed  Google Scholar 

  16. Ji JS, Zhu A, Lv Y, Shi X. Interaction between residential greenness and air pollution mortality: analysis of the Chinese longitudinal healthy longevity survey. Lancet Planet Health. 2020;4:e107–15. https://doi.org/10.1016/S2542-5196(20)30027-9.

    Article  PubMed  PubMed Central  Google Scholar 

  17. Kim S, Kim H, Lee JT. Interactions between ambient air particles and greenness on cause-specific mortality in seven Korean metropolitan cities, 2008-2016. Int J Environ Res Public Health. 2019;16:1866.

    Article  Google Scholar 

  18. Yitshak-Sade M, Kloog I, Novack V. Do air pollution and neighborhood greenness exposures improve the predicted cardiovascular risk? Environ Int. 2017;107:147–53. https://doi.org/10.1016/j.envint.2017.07.011.

    Article  CAS  PubMed  Google Scholar 

  19. Lv Y, Mao C, Yin Z, Li F, Wu X, Shi X. Healthy Ageing and Biomarkers Cohort Study ( HABCS ): a cohort profile; 2019. p. 1–8.

    Google Scholar 

  20. Van Donkelaar A, Martin RV, Brauer M, Hsu NC, Kahn RA, Levy RC, et al. Global estimates of fine particulate matter using a combined geophysical-statistical method with information from satellites, models, and monitors. Environ Sci Technol. 2016;50:3762–72.

    Article  Google Scholar 

  21. Li T, Zhang Y, Wang J, Xu D, Yin Z, Chen H, et al. All-cause mortality risk associated with long-term exposure to ambient PM2·5 in China: a cohort study. Lancet Public Health. 2018;3:e470–7. https://doi.org/10.1016/S2468-2667(18)30144-0.

    Article  PubMed  Google Scholar 

  22. Friedewald WT, Levy RI, Fredrickson DS. Estimation of the concentration of low-density lipoprotein cholesterol in plasma, without use of the preparative ultracentrifuge. Clin Chem. 1972;18:499–502. https://doi.org/10.1093/clinchem/18.6.499.

    Article  CAS  PubMed  Google Scholar 

  23. Cleeman JI. Executive summary of the third report of the National Cholesterol Education Program (NCEP) expert panel on detection, evaluation, and treatment of high blood cholesterol in adults (adult treatment panel III). J Am Med Assoc. 2001;285:2486–97.

    Article  Google Scholar 

  24. Liu L, Cao Z, Lu F, Liu Y, Lv Y, Qu Y, et al. Vitamin D deficiency and metabolic syndrome in elderly Chinese individuals: evidence from CLHLS. Nutr Metab. 2020;17:1–11.

    Article  Google Scholar 

  25. Alberti KGMM, Eckel RH, Grundy SM, Zimmet PZ, Cleeman JI, Donato KA, et al. Harmonizing the metabolic syndrome: a joint interim statement of the international diabetes federation task force on epidemiology and prevention; national heart, lung, and blood institute; American heart association; World heart federation; International. Circulation. 2009;120:1640–5.

    Article  CAS  Google Scholar 

  26. Millwood IY, Li L, Smith M, Guo Y, Yang L, Bian Z, et al. Alcohol consumption in 0.5 million people from 10 diverse regions of China: prevalence, patterns and socio-demographic and health-related correlates. Int J Epidemiol. 2013;42:816–27.

    Article  Google Scholar 

  27. Wallwork RS, Colicino E, Zhong J, Kloog I, Coull BA, Vokonas P, et al. Ambient fine particulate matter, outdoor temperature, and risk of metabolic syndrome. Am J Epidemiol. 2017;185:30–9. https://doi.org/10.1093/aje/kww157.

    Article  PubMed  Google Scholar 

  28. Shamy M, Alghamdi M, Khoder MI, Mohorjy AM, Alkhatim AA, Alkhalaf AK, et al. Association between exposure to ambient air particulates and metabolic syndrome components in a Saudi Arabian population. Int J Environ Res Public Health. 2018;15:27.

    Article  Google Scholar 

  29. Matthiessen C, Lucht S, Hennig F, Ohlwein S, Jakobs H, Jöckel KH, et al. Long-term exposure to airborne particulate matter and NO2 and prevalent and incident metabolic syndrome – results from the Heinz Nixdorf recall study. Environ Int. 2018;116:74–82.

    Article  CAS  Google Scholar 

  30. Huang S, Zhang X, Huang J, Lu X, Liu F, Gu D. Ambient air pollution and body weight status in adults: a systematic review and meta-analysis. Environ Pollut. 2020;265:114999. https://doi.org/10.1016/j.envpol.2020.114999.

    Article  CAS  PubMed  Google Scholar 

  31. Yang BY, Qian Z, Howard SW, Vaughn MG, Fan SJ, Liu KK, et al. Global association between ambient air pollution and blood pressure: a systematic review and meta-analysis. Environ Pollut. 2018;235:576–88. https://doi.org/10.1016/j.envpol.2018.01.001.

    Article  CAS  PubMed  Google Scholar 

  32. Zhang Y, Li T, Ma R, Yin Z, Wang J, He MZ, et al. Long-term exposure to ambient fine particulate matter and fasting blood glucose level in a Chinese elderly cohort. Sci Total Environ. 2020;717:137191. https://doi.org/10.1016/j.scitotenv.2020.137191.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  33. Luo YN, Huang WZ, Liu XX, Markevych I, Bloom MS, Zhao T, et al. Greenspace with overweight and obesity: a systematic review and meta-analysis of epidemiological studies up to 2020. Obes Rev. 2020;21:1–28.

    Article  Google Scholar 

  34. Yeager R, Riggs DW, DeJarnett N, Tollerud DJ, Wilson J, Conklin DJ, et al. Association between residential greenness and cardiovascular disease risk. J Am Heart Assoc. 2018;7:e009117.

    Article  CAS  Google Scholar 

  35. Yang BY, Markevych I, Heinrich J, Bloom MS, Qian Z, Geiger SD, et al. Residential greenness and blood lipids in urban-dwelling adults: the 33 communities Chinese health study. Environ Pollut. 2019;250:14–22. https://doi.org/10.1016/j.envpol.2019.03.128.

    Article  CAS  PubMed  Google Scholar 

  36. Yang BY, Markevych I, Heinrich J, Bowatte G, Bloom MS, Guo Y, et al. Associations of greenness with diabetes mellitus and glucose-homeostasis markers: the 33 communities Chinese health study. Int J Hyg Environ Health. 2019;222:283–90. https://doi.org/10.1016/j.ijheh.2018.12.001.

    Article  CAS  PubMed  Google Scholar 

  37. Wei Y, Zhang J, Li Z, Gow A, Chung KF, Hu M, et al. Chronic exposure to air pollution particles increases the risk of obesity and metabolic syndrome: findings from a natural experiment in Beijing. FASEB J. 2016;30:2115–22.

    Article  CAS  Google Scholar 

  38. Dzhambov AM, Browning MHEM, Markevych I, Hartig T, Lercher P. Analytical approaches to testing pathways linking greenspace to health: a scoping review of the empirical literature. Environ Res. 2020;186:109613. https://doi.org/10.1016/j.envres.2020.109613.

    Article  CAS  PubMed  Google Scholar 

  39. Kasdagli MI, Katsouyanni K, de Hoogh K, Lagiou P, Samoli E. Associations of air pollution and greenness with mortality in Greece: an ecological study. Environ Res. 2020;196:110348. https://doi.org/10.1016/j.envres.2020.110348.

    Article  CAS  PubMed  Google Scholar 

Download references

Acknowledgments

The authors thank all the participants and workers of the CLHLS study.

Funding

The Chinese Longitudinal Healthy Longevity Study (CLHLS) datasets analyzed in this paper are jointly supported by the National Key R&D Program of China (2018YFC2000400), National Natural Sciences Foundation of China (72061137004,71490732), Duke/Duke-NUS/RECA(Pilot)/2019/0051, and the U.S. National Institute of Aging of National Institute of Health (P01AG031719). The funders had no role in this study analysis, interpretation of data, or writing the manuscript.

Author information

Authors and Affiliations

Authors

Contributions

J.S.J. and LX.L. conceptualized the study, conducted statistical analysis, drafted and edited the article; Y.ZENG and XM.S. acquired the data; LJ.Y., YB.L., Y.ZHANG., TT.L., CR.H, HD.K., JF.Z., Y. ZENG, XM.S. interpreted the results and revised the article. All authors provided critical insights and reviewed the article. The author(s) read and approved the final manuscript.

Corresponding author

Correspondence to John S. Ji.

Ethics declarations

Ethics approval and consent to participate

The research ethics committees of Duke University and Peking University approved the study (IRB00001052–13074). All participants in the study have given informed consents.

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: Table S1.

Population characteristics between those followed up and lost follow-up.

Additional file 2: Table S2.

Baseline population characteristics across different counties.

Additional file 3: Table S3.

The association between the greenness and air pollution with the 2012 baseline metabolic biomarkers (binary outcome).

Additional file 4: Table S4.

The association between the greenness and air pollution with the metabolic biomarkers (continuous outcome) in the longitudinal analysis.

Additional file 5: Table S5.

The association between air pollution with the metabolic biomarkers (One-year and five-year exposure) in the longitudinal analysis.

Additional file 6: Table S6.

The association between greenness, air pollution with the metabolic biomarkers among the participants with at least one follow-up.

Additional file 7: Table S7.

The association between the greenness and air pollution with the metabolic syndrome and the components (binary outcome) in the longitudinal analysis using the Joint Interim Societies’ definition of MetS for Chinese populations.

Additional file 8: Figure S1.

Study population.

Additional file 9: Figure S2.

The three-way interaction model of PM2.5, NDVI, and gender on glycated serum protein (GSP).

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

Liu, L., Yan, L.L., Lv, Y. et al. Air pollution, residential greenness, and metabolic dysfunction biomarkers: analyses in the Chinese Longitudinal Healthy Longevity Survey. BMC Public Health 22, 885 (2022). https://doi.org/10.1186/s12889-022-13126-8

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12889-022-13126-8

Keywords