Baseline characteristics
The final cohort included 360 children, with 180 records in each treatment group. No statistically significant differences were observed in sex distribution, age, disease duration, family history of allergy, food allergy, AR severity, allergen sensitization profiles, stable asthma, or atopic dermatitis (all p > 0.05, Table 1). Age and disease duration were comparable between the two treatment groups (Figure 1). All patients met the predefined criteria for persistent AR. Baseline pulmonary function was within the predefined normal range in both groups, with no significant between-group differences in FEV1 (96.8 ± 8.7% vs 97.4 ± 8.5% predicted, p = 0.508), FVC (98.1 ± 8.9% vs 98.7 ± 9.1% predicted, p = 0.524), or FEV1/FVC (88.6 ± 4.7% vs 88.9 ± 4.5%, p = 0.532). Treatment adherence was also comparable (95.6 ± 4.8% vs 95.1 ± 5.0%, p = 0.331).
| Variable | Group A (n = 180) | Group B (n = 180) | χ²/t | p-value |
| Gender, male/female | 86/94 | 93/87 | 0.55 | 0.46 |
| Age (years) | 8.84 ± 1.57 | 8.97 ± 1.58 | 0.78 | 0.434 |
| Disease duration (months) | 7.56 ± 1.43 | 7.48 ± 1.38 | 0.54 | 0.59 |
| Family history of allergy | 78 (43.3) | 82 (45.6) | 0.18 | 0.67 |
| Food allergy | 44 (24.4) | 41 (22.8) | 0.14 | 0.71 |
| AR severity, moderate/severe | 35/145 | 38/142 | 0.15 | 0.694 |
| House dust mite sensitization | 121 (67.2) | 124 (68.9) | 0.11 | 0.735 |
| Pollen sensitization | 69 (38.3) | 73 (40.6) | 0.19 | 0.666 |
| Animal dander/mold sensitization | 31 (17.2) | 29 (16.1) | 0.08 | 0.777 |
| Stable asthma | 28 (15.6) | 31 (17.2) | 0.18 | 0.669 |
| Atopic dermatitis | 21 (11.7) | 24 (13.3) | 0.23 | 0.633 |
Table 1: Baseline demographic and clinical characteristics of the two groups.
Values are mean ± SD or n (%), as appropriate. Allergen-sensitization categories are not mutually exclusive.

Figure 1: Comparison of baseline age and disease duration between the two treatment groups. Panel A shows age, expressed in years, and Panel B shows disease duration, expressed in months. Bars represent mean ± SD. Please click here to view a larger version of this figure.
Nasal symptom scores
Baseline TNSS values were comparable between Group A and Group B (10.82 ± 1.38 vs 10.96 ± 1.30 points, p = 0.327). Repeated-measures ANOVA showed significant effects of group (F (1,358) = 4.30, p = 0.039, partial η2 = 0.012), time (F (1,358) = 1919.55, p < 0.001, partial η2 = 0.843), and group × time interaction (F (1,358) = 9.05, p = 0.003, partial η2 = 0.025). Post-treatment TNSS was lower in Group B than in Group A (4.65 ± 2.12 vs 5.32 ± 2.09; mean difference, 0.67 points; 95% CI, 0.24–1.11; Cohen’s d = 0.32; p = 0.003; Figure 2).

Figure 2. Comparison of Total Nasal Symptom Score (TNSS) between groups before and after treatment. Error bars represent SD. * = p < 0.001 versus pre-treatment within the same group; # = p = 0.003 versus Group A at the same time point. Pairwise comparisons were Bonferroni-adjusted. Please click here to view a larger version of this figure.
Non-nasal symptom scores
Baseline TNNSS values were comparable between Group A and Group B (7.38 ± 1.46 vs 7.31 ± 1.49 points, p = 0.668). Repeated-measures ANOVA showed significant effects of group (F(1,358) = 7.88, p = 0.005, partial η2 = 0.022), time (F (1,358) = 2684.15, p < 0.001, partial η2 = 0.882), and group × time interaction (F (1,358) = 5.17, p = 0.024, partial η2 = 0.014). Post-treatment TNNSS was lower in Group B than Group A (1.96 ± 1.43 vs 2.48 ± 1.01; mean difference, 0.52 points; 95% CI, 0.26–0.77; Cohen’s d = 0.42; p < 0.001; Figure 3).

Figure 3: Comparison of Total Non-Nasal Symptom Score (TNNSS) between groups before and after treatment. Error bars represent SD. * = p < 0.001 versus pre-treatment within the same group; # = p < 0.001 versus Group A at the same time point. Pairwise comparisons were Bonferroni-adjusted. Please click here to view a larger version of this figure.
Time to symptom improvement
Based on daily symptom records, all patients in both groups met the predefined ≥50% improvement threshold for each of the four nasal symptoms during the 28-day treatment period; therefore, no observations were censored. Group B had a shorter time to symptom improvement than Group A. The mean between-group difference (Group A minus Group B) was 3.10 days for rhinorrhea (95% CI, 2.77–3.43; Cohen’s d = 1.97; p < 0.001), 0.47 days for nasal pruritus (95% CI, 0.10–0.84; Cohen’s d = 0.26; p = 0.014), 0.64 days for nasal congestion (95% CI, 0.35–0.92; Cohen’s d = 0.47; p < 0.001), and 0.48 days for sneezing (95% CI, 0.07–0.89; Cohen’s d = 0.24; p = 0.023; Table 2).
| Rhinorrhea | Nasal Pruritus | Nasal Congestion | Sneezing |
| Group A | 13.48 ± 1.20 | 12.43 ± 1.68 | 10.28 ± 1.26 | 8.93 ± 2.39 |
| Group B | 10.38 ± 1.87 | 11.96 ± 1.89 | 9.64 ± 1.47 | 8.45 ± 1.48 |
| t-value | 18.7 | 2.48 | 4.43 | 2.28 |
| p-value | p < 0.001 | p = 0.014 | p < 0.001 | p = 0.023 |
Table 2: Time to symptom improvement in the two treatment groups, expressed in days.
The overall event rates were compared using Fisher’s exact test (p = 0.623).
Inflammatory biomarkers
Baseline IgE, EOS, and IL-13 levels were comparable between groups (IgE: 224.80 ± 28.60 vs 229.60 ± 27.90 U/mL, p = 0.108; EOS: 1.34 ± 0.35 vs 1.31 ± 0.39 ×109/L, p = 0.441; IL-13: 154.72 ± 24.86 vs 157.63 ± 25.14 ng/L, p = 0.270). Significant time and group × time interaction effects were observed for all three biomarkers. The group main effect was significant for EOS and IL-13 but not for IgE. For IgE, the group, time, and interaction effects were F (1,358) = 1.57 (p = 0.211, partial η2 = 0.004), F (1,358) = 9416.28 (p < 0.001, partial η2 = 0.963), and F (1,358) = 74.28 (p < 0.001, partial η2 = 0.172), respectively. For EOS, the corresponding effects were F (1,358) = 5.38 (p = 0.021, partial η2 = 0.015), F (1,358) = 3744.79 (p < 0.001, partial η2 = 0.913), and F(1,358) = 35.78 (p < 0.001, partial η2 = 0.091). For IL-13, they were F (1,358) = 7.68 (p = 0.006, partial η2 = 0.021), F(1,358) = 1194.55 (p < 0.001, partial η2 = 0.769), and F (1,358) = 19.76 (p < 0.001, partial η2 = 0.052). Post-treatment IgE, EOS, and IL-13 were lower in Group B; for IL-13, the mean between-group difference was 12.42 ng/L (95% CI, 8.07–16.77; Cohen’s d = 0.59; p < 0.001; Figure 4).

Figure 4: Longitudinal comparison of inflammatory biomarkers between the two treatment groups.
(A) Serum immunoglobulin E (IgE) levels before and after treatment. (B) Peripheral blood eosinophil (EOS) counts before and after treatment. (C) Serum interleukin-13 (IL-13) levels before and after treatment. Error bars represent the standard deviation (SD). * = p < 0.001 versus before treatment within the same group; # = p < 0.001 versus Group A at the same time point. Please click here to view a larger version of this figure.
Adverse events
Group A had 3 documented adverse events (1 nausea/vomiting and 2 somnolence), whereas Group B had 1 documented adverse event (somnolence). The overall adverse-event rates were 1.67% and 0.56%, respectively, with no statistically significant difference (Fisher’s exact test, p = 0.623; Table 3). None of the four documented adverse events resulted in treatment discontinuation, and no discontinuations because of intolerance or lack of efficacy were documented in either group. In the targeted retrospective review, apart from somnolence, no documented cases of insomnia, nightmares or abnormal dreams, irritability, agitation, anxiety, mood or behavioral changes, suicidal ideation, or self-harm were identified. Because these symptoms were ascertained retrospectively rather than with a validated prospective assessment instrument, under-reporting of mild or transient symptoms cannot be excluded.
| Adverse event | Group A (n = 180) | Group B (n = 180) |
| Nausea/vomiting | 1 | 0 |
| Dizziness/headache | 0 | 0 |
| Somnolence | 2 | 1 |
| Diarrhea | 0 | 0 |
| Total | 3 (1.67%) | 1 (0.56%) |
Table 3: Adverse events were recorded in the two treatment groups.
The de-identified participant-level analytical dataset used for the analyses reported in this manuscript is provided as Supplementary File 1, and the corresponding variable dictionary is provided as Supplementary File 2. Direct patient identifiers and calendar dates are excluded from the shared dataset to protect participant privacy. Patient_ID represents an arbitrary de-identified study identifier used solely to link variables within the supplementary dataset and has no analytical meaning. The clinical and outcome variables provided in Supplementary File 1 correspond to those used in the analyses reported in this manuscript.
Supplementary File 1: Participant-level analytical dataset. The dataset contains participant-level data for the 360 children included in the retrospective cohort, including treatment group and regimen, demographic and baseline clinical characteristics, treatment adherence, pre- and post-treatment TNSS and TNNSS values, time to improvement of individual nasal symptoms, inflammatory biomarker measurements, adverse events, neuropsychiatric-event documentation, and follow-up status. Please click here to download this file.
Supplementary File 2: Variable dictionary for the participant-level analytical dataset. The file defines the variables included in Supplementary File 1 and provides the corresponding descriptions, variable types, coding or ranges, units, and assessment time points. Please click here to download this file.