Skip to main content

ORIGINAL RESEARCH article

Front. Immunol., 21 December 2022
Sec. Mucosal Immunity

A diagnostic model for predicting type 2 nasal polyps using biomarkers in nasal secretion

Zaichuan Wang,&#x;Zaichuan Wang1,2†Qiqi Wang&#x;Qiqi Wang1†Su Duan,&#x;Su Duan1,3†Yuling Zhang,Yuling Zhang1,2Limin Zhao,Limin Zhao1,2Shujian Zhang,Shujian Zhang1,2Liusiqi Hao,Liusiqi Hao1,2Yan LiYan Li1Xiangdong Wang,Xiangdong Wang1,2Chenshuo Wang,Chenshuo Wang1,2Nan ZhangNan Zhang4Claus BachertClaus Bachert4Luo Zhang,,*Luo Zhang1,2,3*Feng Lan*Feng Lan1*
  • 1Beijing Key Laboratory of Nasal Disease, Beijing Institute of Otolaryngology, Beijing, China
  • 2Department of Otolaryngology Head and Neck Surgery, Beijing TongRen Hospital, Capital Medical University, Beijing, China
  • 3Department of Allergy, Beijing TongRen Hospital, Capital Medical University, Beijing, China
  • 4Upper Airways Research Laboratory, Department of Otorhinolaryngology, Ghent University, Ghent, Belgium

Background: Predicting type 2 chronic rhinosinusitis with nasal polyps (CRSwNP) may help for selection of appropriate surgical procedures or pharmacotherapies in advance. However, an accurate non-invasive method for diagnosis of type 2 CRSwNP is presently unavailable.

Methods: To optimize the technique for collecting nasal secretion (NasSec), 89 CRSwNP patients were tested using nasal packs made with four types of materials. Further, Th2low and Th2highCRSwNP defined by clustering analysis in another 142 CRSwNP patients using tissue biomarkers, in the meanwhile, inflammatory biomarkers were detected in NasSec of the same patients collected by the selected nasal pack. A diagnostic model was established by machine learning algorithms to predict Th2highCRSwNP using NasSecs biomarkers.

Results: Considering the area under receiver operating characteristic curve (AUC) for IL-5 in NasSec, nasal pack in polyvinyl alcohol (PVA) was superior to other materials for NasSec collection. When Th2low and Th2highCRSwNP clusters were defined, logistic regression and decision tree model for prediction of Th2highCRSwNP demonstrated high AUCs values of 0.92 and 0.90 respectively using biomarkers of NasSecs. Consequently, the pre-pruned decision tree model; based on the levels of IL-5 in NasSec (≤ 15.04 pg/mL), blood eosinophil count (≤ 0.475*109/L) and absence of comorbid asthma; was chosen to define Th2lowCRSwNP from Th2highCRSwNP for routine clinical use.

Conclusions: Taken together, a decision tree model based on a combination of NasSec biomarkers and clinical features can accurately define type 2 CRSwNP patients and therefore may be of benefit to patients in receiving appropriate therapies in daily clinical practice.

Introduction

Chronic rhinosinusitis with nasal polyps (CRSwNP) is an inflammatory disease of the nasal and paranasal cavities, which can be distinguished from chronic rhinosinusitis without nasal polyps (CRSsNP), the other phenotype of CRS, according to patients’ symptoms and by endoscopy and computed tomography (CT) scan. However, following the “one-size-fits-all” treatment strategy in clinic, some CRSwNP patients with similar symptoms may have different therapeutic outcomes, suggesting that there may be different underling pathomechanisms (1). In this regard, endotyping helps in the recognition of the complexity and heterogeneity of disease and selection of appropriate surgical procedures and pharmacotherapy for patients (2).

Type 2 CRSwNP patients are characterized by the presence of staphylococcal enterotoxin immunoglobulin E (SE-IgE) and type 2 cytokines such as IL-4, IL-5 and IL-13, as well as high asthma comorbidity and recurrence (35). Indeed in general, non-type 2 CRSwNP and type 2 CRSwNP require different pharmacological, surgical approaches and biologic treatments (6). “Reboot” approach or extensive surgery, which is not classical functional endoscopic sinus surgery, can significantly reduce the recurrence of type 2 CRSwNP by removing the dysfunctional and inflamed mucosa without interfering with re-epithelialization of sinus wall (2, 7). Biologics targeting antigen-specific IgE, IL-5 or IL-4 receptor are able to control symptoms and complications of severe type 2 CRSwNP (1, 810). Thus, distinguishing type 2 CRSwNP from non-type 2 CRSwNP is likely to benefit specialists in selecting appropriate effective treatments to improve the outcome in patients.

Whilst tissue inflammatory molecular and cellular biomarkers are generally used to distinguish type 2 CRSwNP from non-type 2 CRSwNP patients, a precise technique for diagnosing type 2 CRSwNP prior to treatment is still not available (4, 1113). In this regard the European Forum for Research and Education in Allergy and Airway Diseases (EUFOREA) (14) has described an algorithm based on comorbid asthma and blood eosinophil count, which can define type 2 CRSwNP in Europe with 86.4% sensitivity and 72.7% specificity. However, inflammatory signatures of CRSwNP vary around the world, and a lower ratio of comorbid asthma is noticed in China, Korea and Thailand compared to Europe (15). Thus, whether the EUFOREA (14) recommended algorithm can be used for type 2 CRSwNP patients in Asian countries awaits further investigations. Assessment of local inflammation in nasal secretion (NasSec) or nasal mucus may be a useful method in the diagnosis of type 2 CRSwNP. Reportedly, assessment of local inflammation in NasSec obtained from common meatus or nasal mucus was shown to be non-invasive and reproducible for local immunological analysis (16). Furthermore, inflammatory markers in NasSecs have been used to evaluate the efficacy of biologics in CRSwNP patients receiving Dupilumab treatment (17). However, current sampling techniques to obtain NasSec or nasal mucus vary widely with regard to the material and preparation of the sample for assessment of inflammatory markers; with reference ranges for inflammatory mediators also not been elucidated for routine clinical use (1821).

In this study, we therefore attempted to standardize the process for NasSec collection from CRSwNP patients, with particular reference to the optimal method for collection. Furthermore, we attempted to establish a diagnostic model for routine clinical use, based on the patient’s clinical features and biomarkers of local inflammation in NasSec, to accurately define type 2 CRSwNP patients prior to surgery.

Materials and methods

Patients and study design

Patients with CRSwNP scheduled to undergo nasal surgery were recruited into the study from the Department of Rhinology and Allergy, Beijing TongRen Hospital from June 2019 to May 2021. The diagnosis of CRSwNP was made according to the European Position Paper on Rhinosinusitis and Nasal polyps 2012 guidelines (22). Patients with cystic fibrosis, fungal sinusitis, or primary ciliary dyskinesia were excluded. The diagnosis of asthma was confirmed by pulmonologists according to symptoms and respiratory function tests. Patients were considered to be allergic based on the presence of serum specific IgE to inhalant allergens, including Dermatophagoides pteronyssinus (Der p), Dermatophagoides farinae (Der f), Mugwort, Penicillium notatum, Candida albicans, Alternaria, Cladosporium, and Aspergillus. None of the patients had taken oral or intranasal glucocorticoids or antibiotics within 4 weeks before surgery. Twice-daily nasal irrigation and mucolytics were used for patients to alleviate clinical symptoms if necessary, and leukotriene modifiers and antihistamines were also prescribed to CRSwNP patients with allergy. Prior to surgery, NasSec and blood samples were collected from each CRSwNP patient without acute infection, and nasal polyp tissues were collected during surgery; and all samples were stored by freezing at -80°C until further analysis. This study was approved by the Ethics Committee of Beijing TongRen Hospital (TRECKY2019-027), and written informed consent was obtained from all patients before enrolment in the study.

Subjects were assigned to two cross-sectional groups. The clinical characteristics are presented in Table 1. In brief, 89 patients were included in Group 1 and were randomized to receive one of the four nasal packs for NasSec collection to evaluate the appropriate nasal pack for NasSec collection. Based on that, nasal secretion of 142 patients in Group 2 was collected using the appropriate material (PVA, as mentioned in results section), in the meanwhile, histologic samples were obtained for further analysis. To estimate appropriate nasal pack for NasSec collection, the sample size was calculated using linear correlation analysis (correlation coefficient = 0.7), and in developing a diagnostic model (machine learning), area under the curve (AUC) was regarded as the main index when calculating the sample size.

TABLE 1
www.frontiersin.org

Table 1 Patient characterization.

Clinical features of CRSwNP patients

The clinical characteristics documented for these patients included nasal symptoms, body mass index (BMI), endoscopic and CT assessment, peripheral blood cell counts, histological evaluation and grading eosinophilia in nasal smears as follows: i. The nasal symptoms were determined using visual analogue scale (VAS) score. Each symptom was marked by patients on a 10-cm scale. Zero means no discomfort and 10 means the most unbearable symptoms. ii. BMI was calculated as the weight in kilograms divided by the square of the height in meters for each patient. iii. Bilateral nasal polyp score was evaluated according to the 5-point scale (23). The overall CT score was determined according to the Lund-Mackay staging system (24); and the ratio of total ethmoid sinus score to total maxillary score (E/M ratio) was also calculated. iv. A complete blood cell count in CRSwNP patients without infection was performed before surgery; and nasal smears collected from the middle portion of the inferior turbinate in both nostrils were stained with May-Grünwald-Giemsa stain for assessment of eosinophils by light microscopy. v. Polyp tissue sections were stained with hematoxylineosin (H&E), eosinophils, neutrophils, plasma cells and lymphocytes were assessed by bright-field light microscopy (Olympus, BX51) at ×400 magnification. Immune cells were counted in 3 non-overlapping fields in the lamina propria by two pathologists. The numbers of inflammatory cells were expressed as the percentage of total inflammatory cells present. vi. Eosinophilia in nasal smears were graded according to the literature (25), resulting in scores of 0 to 4: 0, none eosinophils; 0.5, occasional eosinophils, 1, 1-5 eosinophils/high power field (HPF) or small eosinophils clumps; 2, 5-15 eosinophils/HPF or larger eosinophils clumps; 3, 5-15 eosinophils/HPF or clumps of eosinophils that do not cover the entire field; 4, >20 eosinophils/HPF with clumps covering the entire field.

NasSec collection

NasSec were collected using four types of commercially available nasal packs of the same size (2.5 cm *0.5 cm *0.67 cm) and made from different materials; including polyvinyl acetate (PVAc, IVALON, Fabco®, New London, CT, USA), cross-linked polyvinyl alcohol (cl-PVA, MEROCEL®, Medtronic, Florida, USA), polyvinyl alcohol (PVA, JCMED, Tianjin, China), and polyvinyl fluoride (PVF, Kuaijie, Beijing YINGJIA Co., Beijing, China). Generally, all nasal packs were made from organic polymers. PVA is a water-soluble synthetic polymer; PVAc is the precursor for PVA with poor water absorption; PVF is made from PVA by reaction with formaldehyde with the feature of water resistance and solubility of non-polar solvents; cl-PVA is a PVA sponge with polyethylene film. Prior to collection of NasSec, the nasal cavity was cleaned by wiping with cotton wool soaked in 0.9% sodium chloride (NaCl) to clear any pre-secreted secretions, and then the specified nasal pack was visually placed between the inferior turbinate and the septum in each nostril for 5 min without local or general anaesthesia. No adverse effects or intolerance were reported. Afterwards, the nasal pack was carefully retrieved using forceps and transferred into 15 mL tube containing 3 mL 0.9% NaCl, for 2 hours at 4°C to dissolve the NasSec in the saline as indicated by Berings and colleagues (19). At the end of this period, the contents of the tube were transferred into a 5 mL syringe and the secretion was recovered from the nasal pack using the plunger, followed by centrifugation at 1500 g at 4°C for 15 min. The supernatant from each sample was collected and aliquoted before storage at -80°C until further analysis. To evaluate the standard IgE recovery rate of different nasal packs, different doses (10, 20, 40, 80, 160, 320, 640 ng/mL) of standard IgE (human IGHE, Sino Biological, Beijing, China) were diluted into 3mL 0.9% NaCl. And then nasal packs were merged into standard IgE solution for 2h, IgE was measured in the solution after centrifuge.

Measurement of inflammatory markers

Inflammatory cytokines were measured in nasal tissue and NasSec samples. Prior to assessment, nasal polyp tissues were suspended in 0.9% NaCl solution with a complete protease inhibitor cocktail (Roche, Mannheim, Germany) based on tissue weight, and homogenized mechanically using a Tissue Lyser LT (Qiagen, Hilden, Germany), as described previously (26). The presence of total IgE and eosinophil cationic protein (ECP) was measured in tissue homogenates, NasSecs or standard IgE in 0.9%NaCl using the ImmunoCAP system (Immunodiagnostics; Thermo Fisher Scientific, Uppsala, Sweden), according to manufacturer’s instructions. The inflammatory markers IL-5, CCL5, CCL24, CCL26, IFN-γ, IL-17 and periostin were assessed by Luminex xMAP suspension array technology (R&D, Minneapolis, Minnesota, USA) in a Bio-Plex 200 system (Bio-Rad, MI). The cut off value of cytokines was defined according to half of the detection limit of commercial kit as described in the literature (11) (0.27 pg/mL for IL-5). Concentrations of the biomarkers in tissue homogenates were expressed as mass versus volume after multiplication by a homogenization dilution factor of 11. Values for total IgE ≥ 112.2 KU/L in tissue were regarded as positive (4).

CRSwNP patients clustering by Partitioning Around Medoids algorithm

Tissue inflammatory biomarkers such as IL-5, ECP and total IgE in polyp tissue, comorbid-asthma, peripheral blood eosinophil count, and revision surgery were used for type 2 CRSwNP clustering analysis as documented in the literature (2729). Continuous variables were transformed to Z-score, while categorical variables were expressed as 0 or 1, and a ‘gower’ dissimilarity matrix, which allows both continuous and categorical data, was calculated. The subjects were sorted into groups using the PAM algorithm, with cluster numbers ranging from 2 to 15 generated. The optical number of clusters was determined by the “elbow” (maximum change) in the wss plot or by the maximized average silhouette (Supplementary Figure 1). To visualize distinction among clusters, principal component analysis (PCA) was performed by plotting subjects in first two dimensions after multidimensional scaling. A modified heat-map was used for descriptive characterization of individual clusters.

Machine learning

To establish a prediction model for type 2 CRSwNP model, two common machine learning algorithms (decision tree and logistic regression) were used as described before (30). Eighteen clinical features and biomarkers, as described in Table 2, were used for machine learning. Missing values were imputed using the mean value from 5 nearest neighbouring values (31). The naive random over sampling method was applied to adjust for imbalance in two groups. For logistic regression model, the least absolute shrinkage and selection operator (LASSO) regulation was performed for feature selection. 5-folds cross-validation was implemented to avoid overfitting. We used the AUC to evaluate the performance of the model and the model in which the cross-validation AUC was within one standard error of the maximum was selected. When λ = 0.052, 5 factors were screened out by the logistic regression model; including IL-5 and CCL26 in nasal secretion, comorbid asthma, revision surgery and blood eosinophil count; exhibiting a high goodness-of-fit in the validation dataset (cross-validation AUC = 0.90, Supplementary Figures 2A, B). Each decision node of a decision tree corresponds to a single predictor variable and a split cut-off on that variable. The cut-off of variables was determined according to the minimized Gini index when building the decision tree. Gini index formula was giniindex=1P12P22 , which P1 and P2 represent the proportion of Th2high CRSwNP in two subcategories. Pre-pruning was used to modify the decision tree model by limiting the leaf nodes. Optimal leaf nodes were generated based on the AUC value in 5-folds cross-validation. Model with 4 leaf nodes was selected, showing a high cross-validation AUC (0.89, Supplementary Figure 2C). Machine learning was performed using the Scikit-learn, imbalanced-learn and glmnet package.

TABLE 2
www.frontiersin.org

Table 2 Demographic and clinical characteristics of CRSwNP patients investigated for development of the predicting model.

Statistical analysis

The Chi-square or Fisher’s exact tests were used for comparing categorical variables. Continuous variables were analyzed by the Mann-Whitney U test or Kruskal-Wallis H test. A Benjamini-Hochberg tail-area-based correction for multiple comparisons was applied to strictly control for false positive. Correlations for cytokine concentrations between NasSec and tissue homogenate were analyzed by Pearson correlation and p value less than 0.05 was considered statistically significant. The diagnostic ability of each biomarker or predicting model was calculated based on the AUC. Statistical analysis was performed using Python programming language (version 3.8.8, Python Software Foundation) and R software (version 4.0.3, R Development Core Team). Data were plotted using the ggplot2 and pROC package for R.

Results

NasSec collected by PVA nasal pack can be used for immune assessments in CRSwNP patients

NasSec was collected from 89 CRSwNP patients (Group 1) using four different nasal packs. The clinical characteristics of patients were shown in Table 1. Association assessments indicated that total IgE in tissue homogenate is positively correlated to that in NasSec collected by PVA and PVAc (Figure 1A). Nevertheless, there was a close relationship between IL-5 in tissue homogenate and IL-5 in NasSec only collected by PVA (r=0.376; p=0.0636) (Figure 1A). In terms of total IgE measurement, the AUCs of NasSec collected by PVF (AUC=0.773) were superior to PVA (AUC=0.651) (Figure 1B). When comparing the sensitivity and specificity of IL-5 measurement in NasSec, PVA was found to be superior to other materials for NasSec collection (AUC=0.870) (Figure 1B). Due to the low sensitivity and specificity of total IgE measurement in NasSec, we also took different doses of standard IgE in 0.9% NaCl as positive controls to clarify the recovery rate of IgE using four nasal packs. As shown in Figure 1C, both PVA and PVF had higher recovery rates of standard IgE than PVAc and cl-PVA. Taking together, PVA was selected for the appropriate nasal packs for NasSec collection to evaluate type 2 inflammation.

FIGURE 1
www.frontiersin.org

Figure 1 Comparison of IL-5 and IgE levels in nasal secretions collected by 4 types of nasal packs. (A) Correlation between total IgE and IL-5 (Group1, n = 89) in nasal secretion and tissue collected by 4 types of nasal packs. (B) Receiver operating characteristic curves of IL-5, total IgE in nasal secretion (Group1, n = 89), using concentrations in tissue homogenate as gold standards. (C) Correlation and recovery rate of standard IgE in nasal secretion collected by four types of nasal packs. PVA, polyvinyl alcohol; PVAc, polyvinyl acetate; cl-PVA, cross-linked PVA; PVF, polyvinyl fluoride.

Type 2 CRSwNP patients defined by clustering analysis using tissue inflammatory biomarkers

The clinical characteristics of patients in Group 2 are shown in Table 1. Clustering 142 CRSwNP patients based on their tissue inflammatory biomarkers resulted in 8 separate clusters (Figure 2A). Clusters 1 and 2 were characterized by a low or undetectable concentration of IL-5, ECP and total IgE, and were thus designated as Th2lowCRSwNP (Figure 2B). In clusters 3 to 8, CRSwNP patients with highest levels of type 2 biomarkers were regarded as Th2highCRSwNP (Figure 2B). Among the 8 clusters, CRSwNP patients in cluster 2 expressed an elevated level of tissue IL-17 and a high proportion of patients had revision surgery. Cluster 6 demonstrated the highest level of local type 2 biomarkers and type 2 systemic features such as high blood eosinophils count, serum IgE, and proportion of comorbid asthma and allergy (Figure 2B). Generally, the scores of nasal blockage and smell loss were higher in Th2highCRSwNP patients than that in Th2lowCRSwNP patients. As expected, polyp tissue of Th2highCRSwNP patients had more eosinophils and fewer lymphocytes in comparison to the tissue of Th2lowCRSwNP patients (p<0.001, Table 2).

FIGURE 2
www.frontiersin.org

Figure 2 Clustering of 142 CRSwNP patients (Group 2). (A) Principal component analysis (PCA) plot of patients plotted in a 2-dimensional space after multidimensional scaling. (B) Modified heat map of data relevant to each cluster. Rows define clusters of patients with CRSwNP. Categorical variables (asthma, allergy and revision surgery) are shown as positive ratio. Continuous variables (tissue biomarkers, blood eosinophils count and blood total IgE) are shown as geometric mean. Each q-value represents Benjamini-Hochberg adjusted p-value.

NasSec IL-5, NasSec total IgE, comorbid asthma and blood eosinophils counts are useful markers for diagnosis of Th2highCRSwNP

Th2highCRSwNP patients in Beijing were predicted following the same algorithm recommended by EUFOREA (14), which confirmed the sensitivity and specificity to be 85.1% and 69.1% respectively (Supplementary Figure 3). Considering that, biomarkers in NasSec were included in to improve the accuracy of type 2 CRSwNP patient definition. Assessment of the inflammatory biomarkers analyzed in NasSec in different CRSwNP clusters (Group 2) further demonstrated the concentrations of IL-5 to be consistently significantly elevated in all Th2high clusters (cluster 3 to 8), but not in Th2low clusters (cluster 1 to 2) (Figure 3A). ROC analysis further demonstrated that NasSec IL-5, NasSec total IgE, comorbid asthma and blood eosinophils counts were good predictors of Th2highCRSwNP (Figures 3B-D).

FIGURE 3
www.frontiersin.org

Figure 3 NasSec and clinic markers for diagnosis of Th2highCRSwNP. (A). Modified heat map of biomarkers in nasal secretions from each cluster (Group 2, n = 142). (B-D). Receiver operating characteristic curve (ROC) analysis for prediction of Th2highCRSwNP (Group 2, n = 142). AUC: area under ROC curve; E/M: ratio of Lund-Mackay scores for the ethmoid sinus and maxillary sinus.

A decision tree prediction model using NasSec IL-5 in combination with comorbid-asthma and blood eosinophil count can accurately predict Th2highCRSwNP

In order to better predict Th2highCRSwNP using biomarkers in NasSec, two machine learning models were established based on clinical features and biomarkers in Group 2. A logistic regression model was expressed as y = 0.0070 × IL-5 in NasSec (pg/mL) + 0.00138 × CCL26 in NasSec (pg/mL) + 2.73 × asthma + 2.44 × blood eosinophil count (109/L) – 0.3 × previous sinus surgery (with presence of asthma as 1 and absence of asthma as 0; with previous sinus surgery as 1 and without previous sinus surgery as 0); and a decision tree model was established as shown in Figure 4A. Considering the relatively low sensitivity of the logistic regression model (73.6%), the pre-pruned decision tree was finally selected as the preferred prediction model; with an absence of asthma comorbidity plus blood eosinophil count ≤ 0.475*109/L and IL-5 concentration ≤ 15.04 pg/mL in NasSec shown to be suitable for distinguishing Th2lowCRSwNP patients from Th2highCRSwNP patients. Furthermore, the AUC value for this prediction model was 0.90, with sensitivity and specificity of 83.9% and 89.1%, respectively (Figure 4B).

FIGURE 4
www.frontiersin.org

Figure 4 Machine learning models for non-invasive diagnosis of Th2highCRSwNP. (A) The decision tree model for predicting type 2 CRSwNP (Group 2, n = 142); IL-5 = concentration of the cytokine in nasal secretion collected by polyvinyl alcohol (PVA) nasal packs (pg/mL). EOS# = peripheral blood eosinophil absolute count (109/L). (B) Receiver operating characteristic curve (ROC) analysis of logistic regression and decision tree model for prediction of Th2highCRSwNP (Group 2, n = 142).

Discussion

This study has demonstrated that Th2highCRSwNP patients can be identified by cluster analysis using reported clinical variables and biomarker cytokines. Using these parameters, we demonstrated for the first time that a decision tree based on presence or absence of comorbid asthma, the number of blood eosinophils and IL-5 concentration in NasSec collected by PVA nasal packs provides an accurate way of distinguishing Th2highCRSwNP patients from Th2lowCRSwNP patients prior to treatment.

Indeed, NasSec is easier to obtain and a good non-invasive alternative to nasal tissue in assessing local inflammation. However, the biomarkers in nasal secretions and in nasal tissue are not always consistent (32). The techniques for collection and processing of NasSec vary widely and may influence the final analysis. In this regard, this study is the first to validate the common meatus NasSec collection procedure by comparing the suitability of four types of commonly employed nasal packs in Group 1 patients and to demonstrate that nasal packs using PVA were the most efficient for collecting NasSec for analysis of the inflammatory cytokine IL-5. In terms of total IgE, the NasSec collected by PVA and PVAc was high correlated with that in nasal tissue, while PVA and PVF had higher recovery rates of standard IgE. Thus, in light of this evidence, PVA was eventually chosen as the most appropriate material for collection of NasSecs to evaluate local inflammation.

The AUC of the algorithm recommended by EUFOREA (14) to predict type 2 CRSwNP in China was found to be 0.85, similar to that in Europe. However, the specificity (69.1%) of this algorithm is relatively low. Although we cannot parallelly compare our 8 clusters to already known clusters described in the literature (11) due to the exclusion of CRSsNP patients during the subject enrollments. Still, we noticed that a Th2lowCRSwNP cluster (cluster 2) in our study has a high recurrence rate, which is similar to the described Th2lowCRSwNP clusters in Tomassen et al.’s study (11). Among all Th2highCRSwNP clusters (cluster 3-8), patients in cluster 7 and 8 are characterized by a higher level of IL-5, ECP and total IgE, while 92% of patients in Th2highCRSwNP clusters are less type 2 inflammation. It is further evident that the inflammatory signatures of CRSwNP varied in the whole world. Thus, to establish an easy-to-use diagnostic model with higher specificity, we included selected a variety of inflammatory cytokines, CCL26, CCL5 and periostin in NasSec, CT score, and useful clinical features as variables for type 2 CRSwNP diagnosis in our model. Our data demonstrated that the use of NasSec IL-5 alone or peripheral blood eosinophil count could assess nasal type 2 inflammation (AUC for IL-5 alone = 0.748; AUC for blood eosinophil count = 0.776). Eventually, a pruned decision tree model, which provided a higher AUC (0.90), was therefore established to predict type 2 CRSwNP accurately.

This study is, however, somewhat limited in that the selection of the specific inflammatory cytokines for identification of type 2 CRSwNP was generally based on the data available in the literature, and that these cytokines were mainly evaluated in the nasal secretion for type 2 identification. In particular, some typical type 2 related biomarkers such as IL-4, IL-13 and alarmins were not included in the present study as IL-4 and IL-13 are not easily detectable in NasSec or tissue of same type 2 CRSwNP patients. Furthermore, as only CRSwNP patients from Beijing were enrolled in this study, the applicability of our findings to other regions in the world remains to be elucidated.

Conclusion

This is the first study to compare the suitability of different materials in nasal packs for collection of NasSec, and to demonstrate that PVA appears to be the most appropriate material in this regard. This study has further demonstrated that with comorbid asthma, the number of blood eosinophils > 0.475*109/L or IL-5 in NasSec > 15.04 ng/mL can be used to define Th2highCRSwNP. This may be of particular clinic relevance because the ability to accurately diagnose type 2 CRSwNP or non-type 2 CRSwNP prior to commencement of treatment may aid in the selection of appropriate surgical procedures and pharmacotherapy accordingly.

Data availability statement

The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.

Ethics statement

The studies involving human participants were reviewed and approved by the Ethics Committee of Beijing TongRen Hospital (TRECKY2019-027), and written informed consent was obtained from all patients before enrollment in the study.

Author contributions

FL and LuZ conceived and designed the experiments. SD, CW and XW collected and handled the samples. QW, YZ, LiZ, and SZ performed the experiments. SD, QW, FL, LH and ZW analyzed the data. SD, FL, QW, ZW, NZ, CB, LiZ and LuZ drafted and revised the manuscript. All authors contributed to the article and approved the submitted version.

Funding

This work was supported by grants from the National Nature Science Foundation of China (81970851, 82271140), the Beijing Nova Program of Science and Technology (Z191100001119117); the Young top-notch talent (2018000021223ZK12); the Beijing Hospitals Authority Youth Programme (QML20180201); CAMS Innovation Fund for Medical Sciences (2019-I2M-5-022); Public Welfare Development and Reform Pilot Project (2019-10) and the National Key Research and Development Program of China (2022YFC2504100).

Acknowledgments

We gratefully acknowledge Prof. Guorong Wu, Southwest University, China, for guidance on establishing the machine learning model and appropriate statistical analysis.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fimmu.2022.1054201/full#supplementary-material

References

1. Bachert C, Mannent L, Naclerio RM, Mullol J, Ferguson BJ, Gevaert P, et al. Effect of subcutaneous dupilumab on nasal polyp burden in patients with chronic sinusitis and nasal polyposis: A randomized clinical trial. JAMA (2016) 315(5):469–79. doi: 10.1001/jama.2015.19330

PubMed Abstract | CrossRef Full Text | Google Scholar

2. Alsharif S, Jonstam K, van Zele T, Gevaert P, Holtappels G, Bachert C. Endoscopic sinus surgery for type-2 crs wnp: An endotype-based retrospective study. Laryngoscope (2019) 129(6):1286–92. doi: 10.1002/lary.27815

PubMed Abstract | CrossRef Full Text | Google Scholar

3. Van Zele T, Gevaert P, Watelet JB, Claeys G, Holtappels G, Claeys C, et al. Staphylococcus aureus colonization and ige antibody formation to enterotoxins is increased in nasal polyposis. J Allergy Clin Immunol (2004) 114(4):981–3. doi: 10.1016/j.jaci.2004.07.013

PubMed Abstract | CrossRef Full Text | Google Scholar

4. Wang X, Zhang N, Bo M, Holtappels G, Zheng M, Lou H, et al. Diversity of Th cytokine profiles in patients with chronic rhinosinusitis: A multicenter study in Europe, Asia, and Oceania. J Allergy Clin Immunol (2016) 138(5):1344–53. doi: 10.1016/j.jaci.2016.05.041

PubMed Abstract | CrossRef Full Text | Google Scholar

5. DeConde AS, Mace JC, Levy JM, Rudmik L, Alt JA, Smith TL. Prevalence of polyp recurrence after endoscopic sinus surgery for chronic rhinosinusitis with nasal polyposis. Laryngoscope (2017) 127(3):550–5. doi: 10.1002/lary.26391

PubMed Abstract | CrossRef Full Text | Google Scholar

6. Bachert C, Marple B, Schlosser RJ, Hopkins C, Schleimer RP, Lambrecht BN, et al. Adult chronic rhinosinusitis. Nat Rev Dis Primers (2020) 6(1):86. doi: 10.1038/s41572-020-00218-1

PubMed Abstract | CrossRef Full Text | Google Scholar

7. Jonstam K, Alsharif S, Bogaert S, Suchonos N, Holtappels G, Jae-Hyun Park J, et al. Extent of inflammation in severe nasal polyposis and effect of sinus surgery on inflammation. Allergy (2021) 76(3):933–6. doi: 10.1111/all.14550

PubMed Abstract | CrossRef Full Text | Google Scholar

8. Legrand F, Klion AD. Biologic therapies targeting eosinophils: Current status and future prospects. J Allergy Clin Immunol Pract (2015) 3(2):167–74. doi: 10.1016/j.jaip.2015.01.013

PubMed Abstract | CrossRef Full Text | Google Scholar

9. Bachert C, Zinreich SJ, Hellings PW, Mullol J, Hamilos DL, Gevaert P, et al. Dupilumab reduces opacification across all sinuses and related symptoms in patients with crswnp. Rhinology (2020) 58(1):10–7. doi: 10.4193/Rhin18.282

PubMed Abstract | CrossRef Full Text | Google Scholar

10. Bachert C, Sousa AR, Lund VJ, Scadding GK, Gevaert P, Nasser S, et al. Reduced need for surgery in severe nasal polyposis with mepolizumab: Randomized trial. J Allergy Clin Immunol (2017) 140(4):1024–31.e14. doi: 10.1016/j.jaci.2017.05.044

PubMed Abstract | CrossRef Full Text | Google Scholar

11. Tomassen P, Vandeplas G, Van Zele T, Cardell LO, Arebro J, Olze H, et al. Inflammatory endotypes of chronic rhinosinusitis based on cluster analysis of biomarkers. J Allergy Clin Immunol (2016) 137(5):1449–56.e4. doi: 10.1016/j.jaci.2015.12.1324

PubMed Abstract | CrossRef Full Text | Google Scholar

12. Hoggard M, Waldvogel-Thurlow S, Zoing M, Chang K, Radcliff FJ, Wagner Mackenzie B, et al. Inflammatory endotypes and microbial associations in chronic rhinosinusitis. Front Immunol (2018) 9:2065. doi: 10.3389/fimmu.2018.02065

PubMed Abstract | CrossRef Full Text | Google Scholar

13. Liao B, Liu JX, Li ZY, Zhen Z, Cao PP, Yao Y, et al. Multidimensional endotypes of chronic rhinosinusitis and their association with treatment outcomes. Allergy (2018) 73(7):1459–69. doi: 10.1111/all.13411

PubMed Abstract | CrossRef Full Text | Google Scholar

14. Bachert C, Han JK, Wagenmann M, Hosemann W, Lee SE, Backer V, et al. Euforea expert board meeting on uncontrolled severe chronic rhinosinusitis with nasal polyps (Crswnp) and biologics: Definitions and management. J Allergy Clin Immunol (2021) 147(1):29–36. doi: 10.1016/j.jaci.2020.11.013

PubMed Abstract | CrossRef Full Text | Google Scholar

15. Zhang Y, Gevaert E, Lou H, Wang X, Zhang L, Bachert C, et al. Chronic rhinosinusitis in Asia. J Allergy Clin Immunol (2017) 140(5):1230–9. doi: 10.1016/j.jaci.2017.09.009

PubMed Abstract | CrossRef Full Text | Google Scholar

16. Watelet JB, Gevaert P, Holtappels G, Van Cauwenberge P, Bachert C. Collection of nasal secretions for immunological analysis. Eur Arch Otorhinolaryngol (2004) 261(5):242–6. doi: 10.1007/s00405-003-0691-y

PubMed Abstract | CrossRef Full Text | Google Scholar

17. Jonstam K, Swanson BN, Mannent LP, Cardell LO, Tian N, Wang Y, et al. Dupilumab reduces local type 2 pro-inflammatory biomarkers in chronic rhinosinusitis with nasal polyposis. Allergy (2019) 74(4):743–52. doi: 10.1111/all.13685

PubMed Abstract | CrossRef Full Text | Google Scholar

18. Gilles-Stein S, Beck I, Chaker A, Bas M, McIntyre M, Cifuentes L, et al. Pollen derived low molecular compounds enhance the human allergen specific immune response in vivo. Clin Exp Allergy (2016) 46(10):1355–65. doi: 10.1111/cea.12739

PubMed Abstract | CrossRef Full Text | Google Scholar

19. Berings M, Arasi S, De Ruyck N, Perna S, Resch Y, Lupinek C, et al. Reliable mite-specific ige testing in nasal secretions by means of allergen microarray. J Allergy Clin Immunol (2017) 140(1):301–3.e8. doi: 10.1016/j.jaci.2016.11.047

PubMed Abstract | CrossRef Full Text | Google Scholar

20. Kramer MF, Burow G, Pfrogner E, Rasp G. In vitro diagnosis of chronic nasal inflammation. Clin Exp Allergy (2004) 34(7):1086–92. doi: 10.1111/j.1365-2222.2004.01989.x

PubMed Abstract | CrossRef Full Text | Google Scholar

21. Riechelmann H, Deutschle T, Friemel E, Gross HJ, Bachem M. Biological markers in nasal secretions. Eur Respir J (2003) 21(4):600–5. doi: 10.1183/09031936.03.00072003

PubMed Abstract | CrossRef Full Text | Google Scholar

22. Fokkens WJ, Lund VJ, Mullol J, Bachert C, Alobid I, Baroody F, et al. Epos 2012: European position paper on rhinosinusitis and nasal polyps 2012. a summary for otorhinolaryngologists. Rhinology (2012) 50(1):1–12. doi: 10.4193/Rhino12.000

PubMed Abstract | CrossRef Full Text | Google Scholar

23. Gevaert P, Calus L, Van Zele T, Blomme K, De Ruyck N, Bauters W, et al. Omalizumab is effective in allergic and nonallergic patients with nasal polyps and asthma. J Allergy Clin Immunol (2013) 131(1):110–6 e1. doi: 10.1016/j.jaci.2012.07.047

PubMed Abstract | CrossRef Full Text | Google Scholar

24. Lund VJ, Mackay IS. Staging in rhinosinusitus. Rhinology (1993) 31(4):183–4.

PubMed Abstract | Google Scholar

25. Howarth PH, Persson CG, Meltzer EO, Jacobson MR, Durham SR, Silkoff PE. Objective monitoring of nasal airway inflammation in rhinitis. J Allergy Clin Immunol (2005) 115(3 Suppl 1):S414–41. doi: 10.1016/j.jaci.2004.12.1134

PubMed Abstract | CrossRef Full Text | Google Scholar

26. Zhang N, Holtappels G, Claeys C, Huang G, van Cauwenberge P, Bachert C. Pattern of inflammation and impact of staphylococcus aureus enterotoxins in nasal polyps from southern China. Am J Rhinol (2006) 20(4):445–50. doi: 10.2500/ajr.2006.20.2887

PubMed Abstract | CrossRef Full Text | Google Scholar

27. Chung YW, Cha J, Han S, Chen Y, Gucek M, Cho HJ, et al. Apolipoprotein e and periostin are potential biomarkers of nasal mucosal inflammation. a parallel approach of in vitro and in vivo secretomes. Am J Respir Cell Mol Biol (2020) 62(1):23–34. doi: 10.1165/rcmb.2018-0248OC

PubMed Abstract | CrossRef Full Text | Google Scholar

28. Yao Y, Yang C, Yi X, Xie S, Sun H. Comparative analysis of inflammatory signature profiles in eosinophilic and noneosinophilic chronic rhinosinusitis with nasal polyposis. Biosci Rep (2020) 40(2):BSR20193101. doi: 10.1042/BSR20193101

PubMed Abstract | CrossRef Full Text | Google Scholar

29. Soler ZM, Schlosser RJ, Bodner TE, Alt JA, Ramakrishnan VR, Mattos JL, et al. Endotyping chronic rhinosinusitis based on olfactory cleft mucus biomarkers. J Allergy Clin Immunol (2021) 147(5):1732–41.e1. doi: 10.1016/j.jaci.2021.01.021

PubMed Abstract | CrossRef Full Text | Google Scholar

30. Christodoulou E, Ma J, Collins GS, Steyerberg EW, Verbakel JY, Van Calster B. A systematic review shows no performance benefit of machine learning over logistic regression for clinical prediction models. J Clin Epidemiol (2019) 110:12–22. doi: 10.1016/j.jclinepi.2019.02.004

PubMed Abstract | CrossRef Full Text | Google Scholar

31. Takeda K, Sakakibara S, Yamashita K, Motooka D, Nakamura S, El Hussien MA, et al. Allergic conversion of protective mucosal immunity against nasal bacteria in patients with chronic rhinosinusitis with nasal polyposis. J Allergy Clin Immunol (2019) 143(3):1163–75.e15. doi: 10.1016/j.jaci.2018.07.006

PubMed Abstract | CrossRef Full Text | Google Scholar

32. Workman AD, Nocera AL, Mueller SK, Otu HH, Libermann TA, Bleier BS. Translating transcription: Proteomics in chronic rhinosinusitis with nasal polyps reveals significant discordance with messenger rna expression. Int Forum Allergy Rhinol (2019) 9(7):776–86. doi: 10.1002/alr.22315

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: chronic rhinosinusitis with nasal polyps, diagnostic model, nasal secretion, IL-5, blood eosinophils

Citation: Wang Z, Wang Q, Duan S, Zhang Y, Zhao L, Zhang S, Hao L, Li Y, Wang X, Wang C, Zhang N, Bachert C, Zhang L and Lan F (2022) A diagnostic model for predicting type 2 nasal polyps using biomarkers in nasal secretion. Front. Immunol. 13:1054201. doi: 10.3389/fimmu.2022.1054201

Received: 26 September 2022; Accepted: 08 December 2022;
Published: 21 December 2022.

Edited by:

Glenis Kathleen Scadding, University College London, United Kingdom

Reviewed by:

Huabin Li, Fudan University, China
Marek Sanak, Jagiellonian University Medical College, Poland

Copyright © 2022 Wang, Wang, Duan, Zhang, Zhao, Zhang, Hao, Li, Wang, Wang, Zhang, Bachert, Zhang and Lan. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Luo Zhang, dr.luozhang@139.com; Feng Lan, fenglanent@126.com

These authors have contributed equally to this work

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.