3.1. Population characteristics
Table 1 presented an overview of the general characteristics of the study population. The analysis included a total of 331 participants, with 40 individuals exhibiting high NfL levels. Comparing the low NfL and high NfL groups, no significant differences were found in terms of gender or ethnicity. However, there was a significant difference in age between these two groups. Furthermore, the serum concentration of hexanaldehyde was notably higher in participants with high NfL levels compared to those in the low NfL group. The Spearman correlation analysis revealed a partially dense correlation pattern among certain aldehydes within the 331 subjects, as shown in Supplementary Figure S2. The Spearman correlation coefficients ranged from −0.06 to 0.54, indicating various degrees of correlation. Notably, a strong correlation was observed between exposure to butyraldehyde and propanaldehyde (ρ = 0.54), as well as between isopentanaldehyde and propanaldehyde (ρ = 0.54), as depicted in the heatmap.
Table 1 Characteristics of the overall target population according to NfL concentrations.
|
Overall
(n=331)
|
Low NfL
(n=291)
|
High NfL
(n=40)
|
p
|
Age
|
43.00 (33.00, 56.00)
|
42.00 (32.00, 55.00)
|
58.00 (40.14, 67.76)
|
0.002
|
Sex(%)
|
|
|
|
0.417
|
Male
|
159 (51.6)
|
137(50.4)
|
22(63.3)
|
|
Female
|
172(48.4)
|
154(49.6)
|
18(36.7)
|
|
Ethnicity(%)
|
|
|
|
0.789
|
White
|
130(60.3)
|
116(60.6)
|
14(57.3)
|
|
Asian
|
40 (5.2)
|
35(5.1)
|
5 (6.3)
|
|
Black
|
46(9.6)
|
40(9.7)
|
6(9.0)
|
|
Mexican
|
69 (15.0)
|
58(14.5)
|
11(19.6)
|
|
Other
|
46(9.9)
|
42(10.1)
|
4(7.8)
|
|
Benzaldehyde
|
1.28 (0.82, 1.76)
|
1.28 (0.80, 1.77)
|
1.32 (0.90, 1.69)
|
0.877
|
Butyraldehyde
|
0.54 (0.39, 0.70)
|
0.53 (0.39, 0.69)
|
0.60 (0.43, 0.78)
|
0.161
|
Heptanaldehyde
|
0.49 (0.43, 0.58)
|
0.49 (0.43, 0.58)
|
0.51 (0.44, 0.55)
|
0.647
|
Hexanaldehyde
|
2.10 (1.71, 2.62)
|
2.08 (1.70, 2.53)
|
2.59 (2.00, 2.87)
|
0.002
|
Isopentanaldehyde
|
0.43 (0.30, 0.68)
|
0.43 (0.29, 0.67)
|
0.49 (0.32, 1.15)
|
0.14
|
Propanaldehyde
|
1.93 (1.44, 2.48)
|
1.90 (1.40, 2.44)
|
2.30 (1.75, 2.80)
|
0.064
|
Notes: Data are expressed as median (interquartile range), or n (percentage),n was unweighted sample size, percentage was weighted.
Wilcoxon rank sum test and Chi-square tests were used for comparisons among groups. Sampling design complexity is taken into account in all analyses.
Statistical significance (p<0.05).
Abbreviations: NfL ,serum neurofilament light chain.
3.2. Multivariable logistic regression model to assess the association between aldehydes and NfL concentrations
Multivariable logistic regression was utilized to evaluate the impact of individual aldehydes on NfL concentrations. Hexanaldehyde and isopentanaldehyde displayed significant associations with NfL concentrations in the top two quartiles in the unadjusted model. After adjusting for all covariates, hexanaldehyde and isopentanaldehyde continued to display significant associations with NfL concentrations in the upper two quartiles. The other aldehydes did not show any significant association with NfL concentrations (Table 2).
Table 2
Multivariable associations of selected aldehydes with NfL concentrations in US population 2013–2014.
|
Q1
|
Q2
|
Q3
|
Q4
|
p for trend
|
Benzaldehyde
|
|
|
|
|
|
Model 1
|
Reference
|
1.22(0.32,4.69)
|
2.17(0.58,8.14)
|
0.75(0.18,3.19)
|
0.805
|
Model 2
|
Reference
|
1.39(0.24,7.94)
|
2.45(0.35,17.15)
|
0.85(0.13,5.80)
|
0.987
|
Butyraldehyde
|
|
|
|
|
|
Model 1
|
Reference
|
0.38(0.05,2.60)
|
1.03(0.15,6.89)
|
2.41(0.57,10.12)
|
0.229
|
Model 2
|
Reference
|
0.30(0.02,4.19)
|
1.07(0.09,13.27)
|
2.00(0.37,10.89)
|
0.164
|
Heptanaldehyde
|
|
|
|
|
|
Model 1
|
Reference
|
0.73(0.34,1.55)
|
1.68(0.45,6.19)
|
0.97(0.25,3.82)
|
0.774
|
Model 2
|
Reference
|
0.64(0.22,1.80)
|
1.50(0.34,6.64)
|
1.15(0.31,4.24)
|
0.569
|
Hexanaldehyde
|
|
|
|
|
|
Model 1
|
Reference
|
2.93(0.46,18.60)
|
3.60(1.12,11.57)
|
7.54(2.96,19.21)
|
< 0.001
|
Model 2
|
Reference
|
2.98(0.28,32.20)
|
3.45(0.81,14.80)
|
8.31(3.09,22.34)
|
< 0.001
|
Isopentanaldehyde
|
|
|
|
|
|
Model 1
|
Reference
|
4.41(1.48,13.14)
|
2.35(0.51,10.79)
|
4.15(1.66,10.35)
|
0.109
|
Model 2
|
Reference
|
3.47(1.01,11.95)
|
1.88(0.39,9.13)
|
5.04(1.63,15.64)
|
0.061
|
Propanaldehyde
|
|
|
|
|
|
Model 1
|
Reference
|
3.42(0.73,16.05)
|
4.27(0.89,20.38)
|
3.44(0.94,12.65)
|
0.080
|
Model 2
|
Reference
|
3.12(0.49,19.71)
|
4.65(0.68,31.52)
|
3.53(0.75,16.60)
|
0.060
|
Model 1: unadjusted |
Model 2: adjusted for sex, age, ethnicity |
Notes: Data were expressed as median (interquartile range).Sampling design complexity was taken into account in all analyses. |
Statistical significance (p < 0.05). |
Abbreviations: NfL ,serum neurofilament light chain;Q, quartile. |
3.3. BKMR model to assess the association between aldehydes and NfL concentrations
The ln-transformed concentration of each aldehyde was considered as a continuous variable in our study. We applied the BKMR model to examine the combined impact of aldehyde exposures on binary outcomes. Figure 1.a illustrated the overall association between aldehyde mixtures and NfL concentrations.When the percentile of all aldehyde mixtures reached 55% or higher, it demonstrated a notable elevation in NfL concentrations, thereby underscoring a positive correlation between mixed aldehydes and NfL levels. Figure 1.b revealed that isopentanaldehyde and hexanaldehyde exhibited an increase in NfL concentrations when other aldehydes were maintained at median values. The posterior inclusion probabilities (PIP) derived from the BKMR model revealed that isopentanaldehyde had a PIP of 0.855, and hexanaldehyde had a PIP of 0.600,suggesting that isopentanaldehyde and hexanaldehyde made the most substantial contributions to the mixed exposure model (Supplementary Table S1). To further investigate the potential relationship between serum aldehydes, we generated bivariate cross-sections of the exposure-response function. Supplementary Figure S3 confirms that there was no evidence of interaction between any two types of aldehydes based on parallel exposure-response relationships. Consequently, We proceeded to assess the individual effect of each aldehyde by estimating univariate summaries of the change in NfL concentrations associated with a shift from the 25th percentile to the 75th percentile for a single aldehyde, while keeping all other aldehydes fixed at specific thresholds (25th, 50th, or 75th percentile). Our analysis revealed a significant positive effect for isopentanaldehyde and hexanaldehyde (Supplementary Figure S4).
3.4. WQS regression model to assess the association between aldehydes and NfL concentrations
The WQS index was statistically associated with NfL concentrations in the positive direction (OR = 2.46, 95%CI: 1.27–4.73) ,while no statistically significant association was observed in the negative direction (OR = 1.05, 95%CI: 0.70–2.04). Detailed findings from the adjusted models can be found in Table 3. Figure 2 presents the estimated weights of aldehydes for each WQS index. In the positive model, the top three aldehydes with the highest weights were butyraldehyde, isopentanaldehyde, and hexanaldehyde (weighted 0.47, 0.29, and 0.19, respectively). It is worth noting that butyraldehyde differed from the BKMR results, whereas isopentanaldehyde and hexanaldehyde were consistent with the BKMR findings..
Table 3
Association between WQS regression index and NfL concentrations in US population 2013–2014.
Variables
|
OR
|
95% CI of OR
|
p
|
Positive
|
2.46
|
1.27,4.73
|
0.008
|
Negative
|
1.05
|
0.70 ,2.04
|
0.514
|
Notes:Analysis was adjusted for sex, age, ethnicity.OR estimates represent the odds ratios of NfL concentrations when the WQS index was increased by one quartile. |
Statistical significance (p < 0.05). |
Abbreviations: NfL ,serum neurofilament light chain.OR: odds ratio; CI: confidence interval; |
3.5. The results of functional enrichment analysis and the identification of hub genes in aldehyde-induced axonal damage
Figure 3.a illustrated an overlap between genes associated with aldehydes and those linked to axonal damage, comprising a total of 99 genes.The top five enrichment results were inflammatory response,cell activation,regulation of cell-cell adhesion,positive regulation of cell death,regulation of mononuclear cell migration and regulation of neuron death(Figure 3.b &c) . Moreover, Figure 3.d showed the top ten hub genes in aldehyde-induced axonal damage, including TNF, IL6, CCL2, CXCL10, CCL3, CCL5, CXCL8, CXCL1, IL13, and FOS.A considerable number of aldehyde-related diseases documented in the CTD database were found to be associated with axonal damage (Supplementary Figure S5, Figure S6, and Figure S7)
3.6. Mediation analyses results
Furthermore, we conducted mediation analyses to investigate whether inflammation parameters, including WBC, LYMPH, MONO, and NEUT, played a mediating role in the associations between isopentanaldehyde, hexanaldehyde, and NfL concentrations. The result indicated that WBC, MONO, and NEUT had significant mediation effects on the relationships between isopentanaldehyde and NfL concentrations, accounting for 20.12%, 16.79%, and 4.23% of the total effect, respectively (all P<0.05). However, no significant mediation effects were observed for WBC, LYMPH, MONO, and NEUT in the relationships between hexanaldehyde and NfL concentrations (all P>0.05). For more detailed information, please refer to Figure 4 and Supplementary Figure S8.