Participant Characteristics and Genotype Distribution
The clinical baseline data of the two groups of subjects were statistically analyzed. The results indicated that the serum levels of IgE, IL-4, IL-6, and IL-10 showed significant differences. Compared with the non-AR volunteer group, the concentration of serum IgE, IL-4, and IL-6 in the AR patient group was increased, while the IL-10 level was decreased (Table 1). The genotype distribution of miR-146a rs2910164 was detected using the TaqMan-qPCR method. The analysis results showed that the GG genotype (χ2=7.018, OR (95% CI) = 0.471 (0.269-0.825), p=0.008) and the G allele (χ2=7.474, OR (95% CI) = 0.637 (0.506-0.894), p=0.006) of rs2910164 were more prevalent in the AR patient group compared to the non-AR group. The genotype distribution at rs2910164 locus followed the Hardy-Weinberg equilibrium (pHWE = 0.191) (Table 2).
Plasma miR-146a Expression Analysis
The expression level of plasma miR-146a across the two study cohorts was measured. A marked reduction in miR-146a level was observed in AR patients (Figure 1A). The ROC curve analysis demonstrated that miR-146a has a good diagnostic efficacy for AR (Figure 1B). The AUC was 0.891. At the optimal cutoff value of 0.785, the sensitivity was 83.89%, the specificity was 79.44%, and the Youden index was 0.633. Based on the miR-146a rs2910164 genotype, the two groups of subjects were respectively divided into three subgroups: GG, GC, and CC. The miR-146a expression differences in each subgroup were compared. The statistical results showed that in the non-AR group, miR-146a expression showed no significant statistical difference among different genotype subgroups (Figure 1C); In the AR group, the expression level of plasma miR-146a in the GG genotype subgroup was significantly lower than that in the GC and CC genotype subgroups (Figure 1D). The above data indicate that the expression of plasma miR-146a in AR patients carrying the GG genotype of rs2910164 is reduced.
Serum Cytokine Analysis and Correlation with miR-146a Expression
Next, the serum levels of IgE, IL-4, IL-6, and IL-10 in the different genotype subgroups of the subjects were statistically analyzed. The statistical results revealed no significant differences in serum levels of IgE, IL-4, IL-6, and IL-10 among different genotype subgroups (Figure 2A-D); In the AR group, the levels of serum IgE, IL-4, and IL-6 in patients carrying the GG genotype were significantly higher than those in patients with the GC and CC genotypes, while the level of serum IL-10 was lower (Figure 2A-D). Pearson correlation analysis was used to investigate the linear correlations between the expression of plasma miR-146a and the TNSS score, RQLQ score, as well as the serum levels of IgE, IL-4, IL-6, and IL-10. The results indicated that the expression of miR-146a was significantly negatively correlated with the TNSS score, RQLQ score, and serum levels of IgE, IL-4, and IL-6, while it was positively correlated with the serum level of IL-10 (Table 3). These findings suggest that the expression level of miR-146a is closely related to the severity of AR as well as the levels of allergic reactions and inflammation.
Risk Factor Analysis for Allergic Rhinitis
Univariate logistic regression was used to analyze the risk factors of AR, and the results indicated that IgE (OR = 1.719, 95% CI = 1.151-2.568, p = 0.008), IL-4 (OR = 1.602, 95% CI = 1.073-2.392, p = 0.021), IL-6 (OR = 1.566, 95% CI = 1.048-2.341, p = 0.029), IL-10 (OR = 0.649, 95% CI = 0.435-0.968, p = 0.034), and miR-146a expression (OR = 0.149, 95% CI = 0.096-0.233, p < 0.001) had impact on the occurrence of AR. After controlling for confounding factors through multivariate logistic regression analysis, it was statistically found that IgE (OR=1.724, 95% CI = 1.092-2.722, p = 0.019), IL-4 (OR = 1.943, 95% CI = 1.206-3.131, p = 0.006), IL-10 (OR = 0.479, 95% CI = 0.296-0.773, p = 0.003), and miR-146a expression (OR = 0.148, 95% CI = 0.093-0.235, p < 0.001) were independent risk factors for AR (Table 4).

Figure 1: Schematic diagram of the workflow. Please click here to view a larger version of this figure.

Figure 2. Plasma miR-146a expression and diagnostic performance in AR.
(A) Relative plasma miR-146a expression in non-AR volunteers and patients with AR. (B) ROC curve evaluating the diagnostic performance of plasma miR-146a for AR. (C) Plasma miR-146a expression among GG, GC, and CC genotype subgroups in the non-AR group. (D) Plasma miR-146a expression among GG, GC, and CC genotype subgroups in the AR group. Data are presented as mean ± SD. ***p < 0.001. Please click here to view a larger version of this figure.

Figure 3. Serum cytokine concentration of subjects with different genotypes of miR-146a rs2910164 in non-AR and AR groups.
Serum levels of (A) IgE, (B) IL-4, (C) IL-6, and (D) IL-10 in GG, GC, and CC genotype subgroups of the non-AR and AR groups. Data are presented as mean ± SD. *p < 0.05, **p < 0.01, ***p < 0.001; ns, not significant. Please click here to view a larger version of this figure.
| Information | non-AR | AR | p-value |
| Age (years) | 42.37±12.43 | 41.81±11.86 | 0.316 |
| BMI (kg/m2) | 23.72±2.36 | 23.58±2.71 | 0.302 |
| Sex (male, n) | 87 | 97 | 0.641 |
| Smoker | 81 | 94 | 0.929 |
| IgE (IU/ml) | 36.17±10.44 | 116.52±24.38 | <0.001 |
| IL-4 (ng/L) | 1.76±0.61 | 18.74±4.62 | <0.001 |
| IL-6 (ng/L) | 1.29±0.57 | 9.63±3.79 | <0.001 |
| IL-10 (ng/L) | 3.26±0.68 | 1.74±0.73 | <0.001 |
| TNSS | / | 5.92±1.82 | |
| RQLQ | / | 2.43±1.16 | |
| Abbreviation: AR, Allergic rhinitis; BMI, Body Mass Index; IgE, Immunoglobulin E; IL, Interleukin; TNSS, Total Nasal Symptom Score; RQLQ, Rhinoconjunctivitis Quality of Life Questionnaire; p<0.05 indicates statistical significance. |
Table 1: Clinical characteristics of study participants.
Comparison of demographic characteristics, serum cytokine levels, and clinical assessment scores between the non-AR and AR groups. Abbreviations are defined in the table footnote.
| rs2910164 | non-AR (n=180) (%) | AR (n=211) (%) | Χ2 | p-value | OR (95% CI) |
| Codominant model | | | | | |
| GG | 31 (17.22%) | 57 (27.01%) | | | 1 |
| GC | 82 (45.56%) | 96 (45.50%) | 2.832 | 0.092 | 0.637 (0.376-1.079) |
| CC | 67 (37.22%) | 58 (27.49%) | 7.018 | 0.008 | 0.471 (0.269-0.825) |
| Dominant model | | | | | |
| GG+ GC | 113 (62.78%) | 153 (72.51%) | | | 1 |
| CC | 67 (38.22%) | 58 (27.49%) | 4.232 | 0.04 | 0.639 (0.417-0.980) |
| Recessive model | | | | | |
| GG | 31 (17.22%) | 57 (27.01%) | | | 1 |
| GC + CC | 149 (82.78%) | 154 (72.99%) | 5.34 | 0.021 | 0.562 (0.344-0.919) |
| Alleles | | | | | |
| G | 144 (40.00%) | 210 (49.76%) | | | 1 |
| C | 216 (60.00%) | 212 (50.24%) | 7.474 | 0.006 | 0.637 (0.506-0.894) |
| pHWE | 0.494 | 0.191 | | | |
Abbreviation: AR, Allergic rhinitis. OR, Odds Ratio; pHWE>0.05 suggests that the group is consistent with the Hardy-Weinberg equilibrium. p < 0.05 indicates statistical significance.
Note: The GG genotype was set as the reference group (OR = 1.000) for all genetic models. |
Table 2: Association between miR-146a rs2910164 polymorphism and allergic rhinitis.
Distribution of rs2910164 genotypes and alleles in the non-AR and AR groups analyzed using codominant, dominant, recessive, and allele genetic models.
| Items | miR-146a | TNSS | RQLQ | IgE | IL-4 | IL-6 | IL-10 |
| Pearson correlation | 1 | -0.742 | -0.736 | -0.806 | -0.756 | -0.727 | 0.757 |
| p-value | / | <0.001 | <0.001 | <0.001 | <0.001 | <0.001 | <0.001 |
| Abbreviation: TNSS, Total Nasal Symptom Score; RQLQ, Rhinoconjunctivitis Quality of Life Questionnaire; IgE, Immunoglobulin E; IL, Interleukin; p<0.05 indicates statistical significance. |
Table 3: Correlation of plasma miR-146a expression with clinical and inflammatory parameters.
Pearson correlation analysis of plasma miR-146a expression with TNSS score, RQLQ score, serum IgE level, and cytokine levels in study participants.
| Items | Univariate regression analysis | Multivariate regression analysis |
| OR | 95% CI | p-value | OR | 95% CI | p-value |
| Age | 0.984 | 0.661-1.465 | 0.937 | | | |
| BMI | 1.301 | 0.799-2.119 | 0.29 | | | |
| Sex | 1.099 | 0.738-1.638 | 0.641 | | | |
| Smoker | 1.018 | 0.683-1.519 | 0.929 | | | |
| IgE | 1.719 | 1.151-2.568 | 0.008 | 1.724 | 1.092-2.722 | 0.019 |
| IL-4 | 1.602 | 1.073-2.392 | 0.021 | 1.943 | 1.206-3.131 | 0.006 |
| IL-6 | 1.566 | 1.048-2.341 | 0.029 | 1.489 | 0.924-2.353 | 0.088 |
| IL-10 | 0.649 | 0.435-0.968 | 0.034 | 0.479 | 0.296-0.773 | 0.003 |
| miR-146a expression | 0.149 | 0.096-0.233 | <0.001 | 0.148 | 0.093-0.235 | <0.001 |
| Abbreviation: AR, Allergic rhinitis; BMI, Body Mass Index; IgE, Immunoglobulin E; IL, Interleukin. |
Table 4: Logistic regression analysis of factors associated with allergic rhinitis.
Univariate and multivariate logistic regression analyses evaluating demographic variables, serum cytokines, and plasma miR-146a expression as factors associated with AR.
Supplementary File 1: Raw Data Please click here to download this file.