Baseline characteristics of study populations
Weighted distributions of the baseline characteristics were displayed in Table 1, comprising demographic, examination, laboratory, and questionnaire information from NHANES 2011-2018. The mean age of retained participants was 12.9 years old, and non-Hispanic White people represented the majority of the population. Then FENO were quartiled for further analysis. The distributions of gender, ethnicity, BMI, waist, TC, TG and caffeine intake were statistically significant (p < 0.05), whereas the distributions of age, educational background, family income, family smoke status, vitamin B12, vitamin C, vitamin D, vitamin K, folate, iron, copper, zinc selenium, serum iron, LDL and HDL were not statistically distinct (p > 0.05). Relative to groups with higher FENO, groups with lower FENO exhibited a higher caffeine intake.
Table 1. Weighted characteristics of study populations in disaggregated by quartiles of FENO
|
Q1
|
Q2
|
Q3
|
Q4
|
P value
|
n
|
192
|
248
|
242
|
246
|
|
Gender (%)
|
|
|
|
|
0.003
|
Male
|
88 (45.8)
|
139 (56.0)
|
144 (59.5)
|
155 (63.0)
|
|
Female
|
104 (54.2)
|
109 (44.0)
|
98 (40.5)
|
91 (37.0)
|
|
Age (mean (SD))
|
11.84 (3.75)
|
12.23 (3.51)
|
12.34 (3.57)
|
12.57 (3.34)
|
0.193
|
Ethnicity (%)
|
|
|
|
|
<0.001
|
White
|
73 (38.0)
|
87 (35.1)
|
58 (24.0)
|
49 (19.9)
|
|
Black
|
41 (21.4)
|
69 (27.8)
|
85 (35.1)
|
100 (40.7)
|
|
Mexican American
|
43 (22.4)
|
39 (15.7)
|
32 (13.2)
|
36 (14.6)
|
|
Others
|
35 (18.2)
|
53 (21.4)
|
67 (27.7)
|
61 (24.8)
|
|
Education (%)
|
|
|
|
|
0.14
|
Less than 7th grade
|
116 (60.4)
|
135 (54.4)
|
124 (51.2)
|
123 (50.0)
|
|
More than 7th grade
|
76 (39.6)
|
113 (45.6)
|
118 (48.8)
|
123 (50.0)
|
|
Family income (%)
|
|
|
|
|
0.348
|
$0-24,999
|
71 (37.0)
|
71 (28.6)
|
73 (30.2)
|
77 (31.3)
|
|
$25,000-64,999
|
68 (35.4)
|
105 (42.3)
|
85 (35.1)
|
82 (33.3)
|
|
$65,000-99,999
|
29 (15.1)
|
44 (17.7)
|
45 (18.6)
|
48 (19.5)
|
|
$100,000 and more
|
24 (12.5)
|
28 (11.3)
|
39 (16.1)
|
39 (15.9)
|
|
Family smoke (%)
|
|
|
|
|
0.207
|
Yes
|
34 (17.7)
|
44 (17.7)
|
45 (18.6)
|
30 (12.2)
|
|
No
|
158 (82.3)
|
204 (82.3)
|
197 (81.4)
|
216 (87.8)
|
|
BMI (mean (SD))
|
22.23 (6.26)
|
22.84 (6.51)
|
23.24 (7.08)
|
21.68 (5.43)
|
0.039
|
Waist (mean (SD))
|
75.58 (16.75)
|
78.09 (17.60)
|
77.79 (17.70)
|
74.25 (14.19)
|
0.036
|
Vitamin B12 (mean (SD))
|
9.15 (5.12)
|
9.34 (6.32)
|
9.17 (5.75)
|
9.00 (5.97)
|
0.937
|
Vitamin C (mean (SD))
|
144.38 (115.92)
|
141.72 (127.15)
|
147.43 (136.83)
|
161.94 (133.08)
|
0.316
|
Vitamin D (mean (SD))
|
9.97 (7.58)
|
9.94 (7.27)
|
9.39 (6.90)
|
9.69 (7.61)
|
0.82
|
Vitamin K (mean (SD))
|
127.23 (130.91)
|
115.81 (116.15)
|
123.59 (135.40)
|
117.70 (117.57)
|
0.759
|
Folate (mean (SD))
|
980.36 (613.24)
|
1003.10 (565.60)
|
1014.51 (744.01)
|
1056.30 (636.52)
|
0.648
|
Iron (mean (SD))
|
26.19 (12.70)
|
27.52 (14.80)
|
27.79 (15.46)
|
27.10 (13.57)
|
0.677
|
Copper (mean (SD))
|
1.81 (0.77)
|
1.86 (0.97)
|
1.89 (0.92)
|
1.81 (0.84)
|
0.704
|
Zinc (mean (SD))
|
19.33 (9.29)
|
19.76 (14.38)
|
19.62 (10.98)
|
19.39 (9.37)
|
0.976
|
Selenium (mean (SD))
|
183.83 (87.98)
|
180.21 (78.24)
|
186.80 (96.67)
|
180.89 (79.81)
|
0.826
|
Caffeine (mean (SD))
|
72.08 (99.95)
|
61.54 (90.11)
|
57.50 (96.60)
|
43.82 (74.53)
|
0.011
|
Serum iron (mean (SD))
|
14.81 (7.20)
|
15.33 (7.41)
|
14.54 (6.43)
|
14.98 (5.91)
|
0.827
|
TC (mean (SD))
|
161.37 (29.51)
|
158.50 (28.19)
|
160.50 (30.53)
|
153.68 (27.39)
|
0.034
|
TG (mean (SD))
|
103.24 (52.50)
|
97.45 (57.14)
|
97.38 (83.75)
|
80.33 (54.67)
|
0.037
|
LDL (mean (SD))
|
93.11 (26.19)
|
87.52 (28.25)
|
88.60 (30.33)
|
79.63 (21.79)
|
0.092
|
HDL (mean (SD))
|
51.32 (13.41)
|
51.06 (12.82)
|
52.60 (12.12)
|
52.85 (11.91)
|
0.364
|
The associations of caffeine intake and FENO
We utilized three weighted linear regression models to examine the relationship between caffeine intake and FENO in school-aged asthmatic children (Table 2). From the results, we noticed a statistically significant negative relationship between caffeine intake and FENO in all models. For each additional unit of caffeine intake (mg), FENO decreased by a certain amount, with 0.03 (-0.04,-0.02) ppb for model 1 which adjusted no covariates, 0.03 (-0.04,-0.01) ppb for model 2 which adjusted for sex, age, ethnicity, educational background, family income, and 0.03 (-0.05,-0.01) ppb for model 3 adjusted for BMI, waist, TC, TG, HDL based on model 2. The trend test was statistically significant in all models, with p-value 0.0039 for model 2, 0.0017 for model 3, and less than 0.001 for model 1, which indicated caffeine intake was linearly associated with FENO in the models.
Table 2. Three weighted linear regression models explicate the relationship of caffeine intake with FENO
|
Model 1
|
Model 2
|
Model 3
|
β (95% CI) P value
|
β (95% CI) P value
|
β (95% CI) P value
|
ALL
|
-0.03(-0.04,-0.02) <0.001
|
-0.03(-0.04,-0.01) <0.001
|
-0.03(-0.05,-0.01) 0.009
|
Caffeine intake
|
|
|
|
Q1
|
Reference
|
Reference
|
Reference
|
Q2
|
-2.71( -8.84, 3.42) 0.378
|
-1.38( -7.65,4.90) 0.659
|
-5.73(-13.99,2.53) 0.166
|
Q3
|
-6.26(-12.12,-0.41) 0.037
|
-4.97(-11.01,1.07) 0.103
|
-7.24(-15.67,1.20) 0.089
|
Q4
|
-9.97(-14.52,-5.43) <0.001
|
-8.86(-13.9,-3.82) 0.001
|
-12.66(-19.22,-6.1) <0.001
|
P for trend
|
<0.001
|
<0.001
|
0.0013
|
Stratification relationship of caffeine intake with FENO
We further analyzed stratification relationship of caffeine intake and FENO in different subgroups by sex, age, ethnicity, education background, income, BMI, waist, TC and HDL to ensure that the results of the linear regression analysis were reliable (Figure 2). According to stratified analysis results, we discovered that a negative connection of caffeine intake with FENO in the specific populations, who were aged 6-12 years and 13-18 years, Female, less than 7th grade, more than 7th grade, BMI less than 20 kg/m2, waist less than 70 cm, TC 146-174 mg and HDL less than 46 mg. Moreover, no significant difference of the interaction test was found in all subgroup analyses (p for interaction > 0.05).
Assessment of the relative importance of selected variables by the XGBoost model
In the stages of model development and verification, the XGBoost algorithmic model of machine learning was used to assess the relative importance of selected variables associated with FENO. The selected variables consisted of age, education, income, BMI, vitamin A, vitamin B12, vitamin C, vitamin D, vitamin K, folate, iron, copper, zinc, calcium, selenium and caffeine intake. According to the assessment results of the contribution of the XGBoost model to each variable, we discovered that vitamin C, iron, caffeine intake, vitamin B12 and vitamin D were the five most influential factors in FENO (Figure 3, Appendix Table 1). As a relatively critical variable, caffeine intake was subsequently included in the general additive model and segmented regression model to further assess the reliability of the linear regression analysis results in our survey.
Evaluating the linear or nonlinear associations of caffeine intake and FENO
The general additive model (GAM), which has been proven to be extremely sensitive to determining whether the correlation is linear or non-linear, was used to explore linear or nonlinear associations between caffeine intake and FENO to verify the trustworthiness of the results of the regression analysis. As shown in Figure 4, a smooth fit curve was generated from model 3 to illustrate the potential connections. After controlling for other variables, a linear relationship between caffeine intake and FENO was observed in asthmatic children.
In addition, a segmented regression model was implemented to confirm the linearity or non-linearity of the association between caffeine intake and FENO (Table 3, Figure 5). Our study revealed a lack of statistical significance for the inflection point (K=3.38) because the log-likelihood ratio was greater than 0.05. Furthermore, there were no statistically significant distinctions between the one-line model and the segmented regression model. It follows that the one-line model was a more appropriate method to clarify the correlation between caffeine intake and FENO. All the above results suggested a linear and inverse relationship between caffeine intake and FENO.
Table 3. Analysis of the threshold effect of caffeine intake and FENO implementing the two-piecewise linear regression model
|
β (95% CI) P value
|
Model 1 (one-line model)
|
linear effect
|
-0.048(-0.087,-0.009) 0.016
|
Model 2 (two-piecewise linear regression model)
|
Inflection point (K)
|
3.38
|
Caffeine intake < K
|
-1.16(-5.02,2.69) 0.55
|
Caffeine intake > K
|
-0.033(-0.079,0.014) 0.084
|
Log likelihood ratio
|
0.267
|
The model 1 and 2 all adjusted for sex, age, ethnicity, educational background, and income.
|