- Open Access
Distinct type 2-high inflammation associated molecular signatures of chronic rhinosinusitis with nasal polyps with comorbid asthma
Clinical and Translational Allergy volume 10, Article number: 26 (2020)
Patients with chronic rhinosinusitis with nasal polyps (CRSwNP) and comorbid asthma have more severe disease and are difficult to treat. However, the molecular endotypes associated with CRSwNP with comorbid asthma (CRSwNP + AS) are not clear. This study aimed to investigate the characteristics of type 2 inflammation and the molecular signatures associated with CRSwNP + AS.
A total of 195 subjects; including 65 CRSwNP + AS patients, 99 CRSwNP-alone patients, and 31 healthy control subjects; were enrolled in the study. Nasal tissues from patients with CRSwNP + AS, CRSwNP-alone and control subjects were assessed for infiltration of inflammatory cells and concentrations of total IgE. Whole-transcriptome sequencing was performed and differentially expressed (DE) mRNAs and long non-coding RNAs (lncRNAs) and their associated pathways were analyzed. The correlations between type 2 cytokines and local eosinophils, tissue IgE, and transcriptome signatures were evaluated.
Significantly higher local eosinophil infiltration and higher levels of total IgE were found in nasal tissues from CRSwNP + AS patients than in nasal tissues from CRSwNP-alone patients. Furthermore, atopy and recurrence were significantly more frequent in patients with CRSwNP + AS than in patients with CRSwNP-alone (62.5% vs 28.6% and 66.7% vs 26.9%, respectively). RNA sequencing analysis identified 1988 common DE-mRNAs, and 176 common DE-lncRNAs shared by CRSwNP + AS versus control and CRSwNP-alone versus control. Weighted gene coexpression network analysis (WGCNA) identified LINC01146 as hub lncRNA dysregulated in both subtypes of CRSwNP. Overall, 968 DE-mRNAs and 312 DE-lncRNAs were identified between CRSwNP + AS and CRSwNP-alone. Both pathway enrichment analysis and WGCNA indicated that the phenotypic traits of CRSwNP + AS were mainly associated with higher activities of arachidonic acid metabolism, type 2 cytokines related pathway and fibrinolysis pathway, and lower activity of IL-17 signalling pathway. Furthermore, the expression of type 2 cytokines; IL5 and IL13, was positively correlated with local eosinophil infiltration, tissue IgE level, and the expression of DE-mRNAs that related to arachidonic acid metabolism. Moreover, WGCNA identified HK3-006 as hub lncRNA in yellow module that most positively correlated with phenotypic traits of CRSwNP + AS.
Patients with CRSwNP + AS have distinct type 2-high inflammation-associated molecular signatures in nasal tissues compared to patients with CRSwNP-alone.
Chronic rhinosinusitis (CRS), a disease characterized by chronic inflammation of the sinonasal tissue, affects 5.5–28% of the general population . CRS with nasal polyps (CRSwNP) accounts for approximately 20% of all CRS and has greater severity of clinical disease . Asthma is one of the most common chronic inflammatory disorders of the lower airway worldwide with increasing morbidity. Studies have reported that up to 60% of CRSwNP patients have comorbid asthma (CRSwNP + AS), which is one of the most challenging CRS subtypes to treat . Patients with CRSwNP + AS have greater disease severity, higher recurrence rates of nasal polyps after surgery, poorer asthma control and higher costs [4,5,6].
There is evidence that eosinophilic CRSwNP tends to have comorbid asthma more frequently . Formation of IgE, which is independent of the presence of allergy in nasal polyp tissue, is also associated with asthmatic condition in patients with CRSwNP . However, there is still no clear explanation for the association between CRSwNP and asthma, and likewise the pathogenic mechanisms leading to CRSwNP and asthma are uncertain. The united airway concept suggests that the upper and lower airway inflammation share common pathogenic mechanisms and influence each other [9, 10]. There is evidence that several features of inflammatory pattern, disrupted epithelial barrier and airway remodelling are similar in CRSwNP and asthma [11,12,13]. Thus, understanding the molecular relationship between CRSwNP and comorbid asthma may help to reveal the mechanisms that underlie airway chronic inflammation.
CRSwNP is a heterogeneous inflammatory condition with different endotypes . The majority of white patients with CRSwNP in western countries have a type 2 pattern of inflammation characterized by pronounced eosinophilia and high levels of interleukin-4 (IL-4), IL-5 and IL-13 cytokines . In contrast, Chinese patients with CRSwNP have lower type 2 inflammation and show higher degree of type 1/type 3 inflammation [16, 17]. Furthermore, evidence from previous studies on epidemiology and clinical characteristics, suggests that CRSwNP + AS may be considered a subtype of CRSwNP [18,19,20]. Thus, this study aimed to investigate the characteristics of type 2 inflammation and molecular endotypes associated with CRSwNP + AS by whole-transcriptome sequencing. Distinct type 2-high inflammation and its associated transcriptome signatures, indicated by coding mRNAs and long non-coding RNAs (lncRNAs), were found in patients with CRSwNP + AS compared to patients with CRSwNP-alone.
Materials and methods
A total of 195 subjects, including 65 CRSwNP patients with comorbid asthma (CRSwNP + AS), 99 patients with CRSwNP-alone and 31 healthy control subjects were enrolled in series from June 2017 to March 2018 in the Rhinology Department of Beijing TongRen Hospital. Patients with CRSwNP were diagnosed according to the European Position Paper on Rhinosinusitis and Nasal Polyps 2020 guidelines . The diagnosis of comorbid asthma was based on the Global Initiative for Asthma 2019 guidelines. The diagnosis of allergic rhinitis was according to Allergic Rhinitis and Its Impact on Asthma 2016 guidelines. Atopy was confirmed based on positive test for serum antigen-specific IgE (cut-off value, 0.35kUA/L), measured by Immuno-CAP 100 system (Pharmacia, Uppsala, Sweden). Patients undergoing septoplasty because of anatomic variations and without other sinonasal diseases were recruited as control subjects. All subjects were aged 18 to 70 years. Subjects with immunodeficiency, fungal sinusitis, coagulation disorder, neoplasia, pregnancy and aspirin-exacerbated respiratory disease (AERD) were excluded. None of the patients had been treated with corticosteroids, antibiotics or biologics within the 4-week period before surgery, and no patient had symptoms of infection at the time of sampling, apart from symptoms of chronic rhinosinusitis or asthma. Recurrence was defined as the presence of nasal polyps observed under nasal endoscopy, together with at least one symptom (nasal obstruction, rhinorrhea, headache/facial pain, reduction or loss of smell, sleep disturbance/fatigue) lasting at least 1 week, despite appropriate intranasal corticosteroid treatment. The postoperative follow-up period was 18 months for assessment of recurrence/non-recurrence of nasal polyps [1, 21]. All patients were followed-up after weeks 1, 2, 4, and 12, and then once every 3 months for up to 18 months by the same surgeon, who was blinded to all laboratory data . Nasal tissue samples were collected from the inferior turbinate of control subjects and the nasal polyps of CRSwNP patients during surgery. The tissue samples were processed for staining with haematoxylin and eosin, RNA sequencing and ELISA as previously described .The Ethics Committee of Beijing Tongren Hospital approved this study, and all subjects signed informed consent forms prior to enrolment in the study.
Histological evaluation of polyp tissue
Nasal polyp tissues were immediately formalin fixed after surgery, and then dehydrated and embedded in paraffin. Paraffin sections were stained with haematoxylin and eosin (H&E) and processed for histological evaluation. All sections were examined by optical microscopy at × 400 magnification. The absolute numbers and percentages of infiltrating inflammatory cells; including eosinophils, neutrophils, plasma cells, and lymphocytes; were recorded as mean of six non-overlapping regions in each section by two independent pathologists, who were blinded to the study design and clinical background of the patients.
Assessment of total IgE in nasal tissues
Concentrations of total IgE in nasal tissues were assayed using the Human IgE ELISA Kit (Arigo Biolaboratories Corporation, Taiwan). Briefly, fresh nasal tissues were placed into RIPA lysis buffer with 1% protease inhibitor cocktail (Thermo Fisher Scientific) and homogenized using a standard bench-top homogenizer (Qiagen, Valencia, CA). The tissue homogenates were centrifuged and supernatants were collected for IgE analysis, according to the manufacturer’s instructions. All samples were tested in duplicate.
RNA isolation and RNA sequencing
Nasal tissue samples of CRSwNP + AS (n = 10), CRSwNP-alone (n = 10), and control (n = 9) were randomly selected for whole-transcriptome sequencing. The nasal tissue samples collected after surgery were freshly preserved in RNAlater solution (Qiagen, Hilden, Germany) and processed for extraction and purification of total RNA using the RNeasy Kit (Qiagen, Hilden, Germany), according to the manufacturer’s instructions. The quantity and quality of the isolated RNA was determined with NanoDrop 2000 Spectrophotometer (Thermo Fischer Scientific) and 2100 TapeStation Automated Electrophoresis System (Agilent Technologies), and samples with an RNA integrity number of greater than 8.0 were chosen for sequencing. Ribosomal RNA was removed and sequencing libraries were prepared using the rRNA-depleted RNA by NEBNext UltraTM Directional RNA Library Prep Kit (New England Biolabs, USA), following the manufacturer’s instructions. RNA sequencing was performed on the Illumina Hiseq platform and 150 bp paired-end reads were generated by Novogene Bioinformatics Technology Cooperation (Beijing, China).
RNA sequencing data analysis
Adapters and low-quality tail were trimmed from reads prior to read alignment. Clean sequence reads were aligned to the human genome with Hisat2 (v2.0.5). Cufflinks (v2.2.1) was used to assemble transcripts, estimate the abundance of these transcripts, and detect differential expression among samples. For mRNA analyses, the reference genome build GRCh37 was chosen as the annotation references. For lncRNA analyses, the GENCODE v19 database was chosen as the annotation references. Fragments per kilo-base of exon per million fragments mapped (FPKM) of both lncRNAs and mRNAs in each sample was calculated based on the length of the fragments and reads count mapped to this fragment. Differential expression analysis was performed using Cuffdiff software (v2.2.1). An adjusted P < 0.05 plus fold change > 2 was used as the cut-off for significantly differentially expressed (DE) mRNAs (DE-mRNAs) and lncRNAs (DE-lncRNAs).
Pathway analysis of DE-mRNAs and function prediction of DE-lncRNAs
DE-mRNAs were loaded into Enrichr (https://amp.pharm.mssm.edu/Enrichr/)  for pathway enrichment analysis, and the top significantly enriched (P < 0.05) KEGG and Biocarta pathways with high Combined Score, provided by Enrichr, were determined.
In order to predict DE-lncRNA functions, we applied weighted gene coexpression network analysis (WGCNA)  to construct a coexpression network between DE-lncRNAs and their highly correlated DE-mRNAs based on the Pearson correlation coefficient between their normalized expression levels. Briefly, network construction and module detection were performed using the “blockwise Modules” function in the WGCNA package. A coexpression similarity matrix was calculated by computing Pearson correlations between all gene pairs, and then transformed into an adjacency matrix using a soft threshold power (β) equal to 18. A dynamic tree cut algorithm was used to detect groups of highly correlated genes. Modules were defined as the branches cut-off of the tree and each module was labelled in unique colours, of which grey colour contains probes not assigned to any module. The module eigengenes were utilized to represent each module, which was calculated via the first principal component. Using the module eigengenes, the relationship between module and tissue type (CRSwNP-alone and CRSwNP + AS) was estimated.
Enrichr was used to perform the functional enrichment analysis for each module, and the top significantly enriched pathways were determined. In interesting modules that related to the disease condition, top hub genes with high connectivity and edges with weight above a threshold of 0.1 were identified and visualized using the cytoscape network.
All data are presented as medians and interquartile range (IQR) except for age, which is presented as mean ± SD. Data analysis was performed using GraphPad Prism Version 7.0 (GraphPad Software, La Jolla, Calif). All parametric variants were analyzed using Student t tests, and nonparametric variants were analyzed by using Mann–Whitney U tests. The χ2 or Fisher exact test was used for qualitative data. Relationships between variables were evaluated using Spearman correlation analysis. Differences were considered significant at P value < 0.05.
Demographic and clinical characteristics of the subjects
Demographic and clinic characteristics of all participants enrolled in this study are presented in the Additional file 1: Table S1. There was no significant difference with regard to age, gender and smoker status between the 3 groups. Peripheral blood eosinophils and total IgE were increased in both subtypes of CRSwNP patients compared to control subjects. Assessment for the significance of differences between CRSwNP + AS group and CRSwNP-alone group indicated that atopy and recurrence were significantly more frequent in patients with CRSwNP + AS than in patients with CRSwNP-alone (62.5% vs 28.6% and 66.7% vs 26.9%, respectively). Patients with CRSwNP + AS also had a significantly decreased forced exhalation volume in one second (FEV1)/forced vital capacity (FVC) ratio and an increased fractional exhaled nitric oxide (FeNO) compared to patients with CRSwNP-alone. Similarly, CRSwNP + AS patients had a significantly higher percentage of eosinophils (6.45%) and total IgE (143.00 kU/l) in peripheral blood compared to CRSwNP-alone patients (3.00% and 53.90 kU/l, respectively; Fig. 1a, b).
Local features of nasal tissue in different patient group
To identify the local inflammatory patterns of nasal tissue in CRSwNP patients with and without comorbid asthma, the infiltrating eosinophils, neutrophils, plasma cells, and lymphocytes were assessed. Tissue sections stained with H&E demonstrated that patients with CRSwNP + AS had significantly more eosinophils in nasal polyp tissues than patients with CRSwNP-alone (Fig. 1c–e).
Given the relationship between IgE and Type 2 immune response, we further examined the total IgE levels in nasal tissues and found that IgE levels in nasal tissues from CRSwNP + AS patients were significantly higher compared to IgE levels in nasal tissues of CRSwNP-alone patients and controls (Fig. 1f). Moreover, Spearman correlation test showed that the percentage of local eosinophils was positively correlated with the concentration of tissue IgE (Additional file 1: Figure S1).
Whole transcriptome profiling of nasal tissues from CRSwNP + AS and CRSwNP-alone
To identify the gene expression profiles of CRSwNP with and without comorbid asthma, RNA sequencing was performed on nasal tissue samples from control, CRSwNP-alone and CRSwNP + AS patients. Analysis of DE-mRNAs and DE-lncRNAs demonstrated that there were 5218 DE-mRNAs and 2949 DE-lncRNAs between CRSwNP-alone and control, and 2512 DE-mRNAs and 464 DE-lncRNAs between CRSwNP + AS and control (Additional file 1: Figure S2). Additionally, 968 DE-mRNAs and 312 DE-lncRNAs were identified between CRSwNP + AS and CRSwNP-alone. Hierarchical clustering of top 500 differentially expressed genes was shown in Additional file 1: Figure S3.
Common dysregulated genes shared by CRSwNP + AS and CRSwNP-alone
A total of 1988 common DE-mRNAs were shared by CRSwNP-alone and CRSwNP + AS respectively compared to control (Fig. 2a and Additional file 1: Table S2). Assessment of these common DE-mRNAs by pathway enrichment analysis using Enrichr Demonstrated that the top enriched KEGG pathways were associated with cytokine–cytokine receptor interaction, chemokine signalling pathway, staphylococcus aureus infection, asthma, and cell adhesion molecules (Fig. 2b). The top enriched BioCarta pathways were related to eicosanoid metabolism, the co-stimulatory signal during T cell activation and IL-10 anti-inflammatory signalling pathway (Fig. 2b). Likewise, gene ontology enrichment analysis showed that common DE-mRNAs were mainly associated with regulation of immune system process, immune response, response to stimulus and cell–cell signalling (Additional file 1: Figure S4).
Similarly, a total of 176 common DE-lncRNAs were shared by CRSwNP + AS versus control and CRSwNP-alone versus control (Fig. 2c and Additional file 1: Table S2). WGCNA was applied to explore the potential functions of common DE-lncRNAs, and construction of a hierarchical clustering tree demonstrated 9 modules, each of which was labelled with a distinct colour (Fig. 2e). The size of these modules ranged from 41 to 539 genes; with the largest turquoise coloured module (Additional file 1: Figure S5) comprising 45 lncRNAs and 494 mRNAs, which were highly coexpressed. Pathway enrichment analysis showed that the genes of this turquoise module were mainly associated with cytokine–cytokine receptor interaction, asthma, T cell receptor signalling pathway, staphylococcus aureus infection, chemokine signalling pathway, and Th1, Th2, Th17 cell differentiation (Fig. 2d).
Figure 2f shows the top 50 hub genes with high connectivity and edge weigh identified and visualized by cytoscape network. LINC01146, the only one lncRNA of these top hub genes, was significantly up-regulated in both CRSwNP-alone and CRSwNP + AS. Pathway analysis of the coexpressed mRNAs of LINC01146 demonstrated that LINC01146 was mostly associated with T cell receptor signalling pathway, natural killer cell mediated cytotoxicity, Fc gamma R-mediated phagocytosis, and Th1 and Th2 cell differentiation (Additional file 1: Figure S6).
Distinct transcriptome signatures in nasal tissue of CRSwNP + AS
Overall, 212 mRNAs were down-regulated and 756 mRNAs up-regulated in nasal tissues of CRSwNP + AS patients, compared to nasal tissues of CRSwNP-alone patients (Fig. 3a). The 50 most significant DE-mRNAs are shown in Additional file 1: Table S3. Enrichr pathway analysis showed that arachidonic acid metabolism, ECM-receptor interaction, IL-17 signalling pathway, GATA3 participate in activating Th2 cytokine genes, and fibrinolysis pathway were the top significant pathways enriched by DE-mRNAs (Fig. 3b).
Detailed examination of the expression of genes involved in arachidonic acid metabolism demonstrated that 14 of 16 DE-mRNAs related to arachidonic acid metabolism; including PTGS1, TBXAS1, ALOX15B, PLA2G4D, LTC4S, GGT1, GGT5, GPX6, PLA2G4B, PLA2G6, CYP4F8, JMJD7-PLA2G4B, PTGIS and GPX3; were up-regulated in CRSwNP + AS compared to CRSwNP-alone (Fig. 3c).
Severe type 2 inflammation in nasal tissues of CRSwNP + AS
As indicated by above transcriptome data, type 1, type 2 and type 3 related signalling pathways might be differentiated between CRSwNP-alone and CRSwNP + AS. Thus, we investigated the mRNA expression of critical cytokines and their receptors, which indicate the activities of different inflammatory endotypes; namely (1) IFNG and IFNG receptor (IFNGR1 and IFNGR2) for type 1 inflammation; (2) IL5, IL4, IL13 and their receptors (IL5RA, IL5RB, IL4R, IL13RA1 and IL13RA2) for type 2 inflammation; (3) IL17A and IL17A receptor (IL17RA and IL17RC) for type 3 inflammation; as well as (4) TNF, IL1B, IL6 and their receptors (TNFRSF1A, TNFRSF1B, IL1R1, IL1R2 and IL6R) for proinflammatory activity.
We found that the expression of cytokines IL5 and IL13, and receptors IL5RA and IL5RB, indicating augmented type 2 inflammation, was significantly enhanced in both CRSwNP-alone and CRSwNP + AS compared to control (Fig. 3d and Additional file 1: Table S4). Furthermore, nasal tissue from CRSwNP + AS demonstrated significantly higher expression of IL4, IL5 and IL13; and conversely significantly lower expression of IL17A than nasal tissue of CRSwNP-alone.
Assessment of correlations between the significant cytokines and the other significant inflammatory indicators determined above, demonstrated that the expression of both IL5 and IL13 was positively correlated with the percentage of local eosinophils, concentration of tissue IgE and the expression of LTC4S, which reflected an imbalanced arachidonic acid metabolism (Fig. 4a–f). In contrast, the expression of IL17A was negatively correlated with tissue IgE and the expression of LTC4S (Fig. 4g–i).
LncRNA signatures in nasal tissue of CRSwNP + AS
We identified 229 up-regulated and 83 down-regulated DE-lncRNAs in nasal tissues of CRSwNP + AS patients compared to nasal tissues of CRSwNP-alone patients (Fig. 5a). The top 50 significant DE-lncRNAs are shown in Additional file 1: Table S5. A coexpression network constructed based on the expression of DE-lncRNAs and DE-mRNAs using WGCNA demonstrated 7 colour-coded modules underneath the cluster tree (Additional file 1: Figure S7). Assessment of the relationship between each module and tissue type (CRSwNP + AS), estimated using the module eigengenes, demonstrated that all 7 modules were significantly correlated with tissue type changes (Fig. 5b). Of these, the blue module was most negatively correlated with the phenotypic traits of CRSwNP + AS (r = − 0.73, P = 2 × 10−4), and the yellow module the most positively correlated (r = 0.66, P = 0.001).
Overall, there were 31 lncRNAs and 94 mRNAs in the blue module, and only one lncRNA (HK3-006) and 63 mRNAs in the yellow module. Pathway enrichment analysis indicated that the blue module was mainly associated with IL-17 signalling pathway and cytokine–cytokine receptor interaction (Fig. 5c), whereas the yellow module was mainly related to asthma, arachidonic acid metabolism, and signalling pathways such as calcium and JAK-STAT signalling pathway (Fig. 5d). Top 50 hub genes identified from the blue and yellow modules, respectively, are shown in Fig. 5 e, f.
Both CRSwNP and asthma are airway inflammatory disorders, which have serious effects on quality of life. A great proportion of CRSwNP patients have comorbid asthma, which makes it one of the most challenging phenotypes to treat [3, 12]. Furthermore, the molecular endotypes of CRSwNP + AS are not clear. In this regard, to our knowledge this is the first study to investigate the molecular endotypes of CRSwNP + AS compared to CRSwNP-alone, by whole-transcriptome RNA sequencing. Our study demonstrated that type 2-high inflammation patterns and their associated transcriptome features were distinct in nasal tissues of CRSwNP + AS patients compared to CRSwNP-alone patients.
Studies have indicated that CRS is a heterogeneous disease with several unclear endotypes, which are mainly characterized by type 1, 2, and 3 inflammatory patterns [15, 25,26,27]. Furthermore, some recent studies have indicated that the most prevalent endotype in CRSwNP is characterized by type 2 inflammation , and that type 2 inflammation in CRSwNP may be differentiated into moderate and severe type 2 inflammatory patterns according to the intensity of inflammation . In accordance with these studies, the findings for the differentially expressed type 1, type 2 and type 3 inflammation-related genes and pathways in the present study suggest that the CRSwNP + AS endotype is associated with a more severe type 2 inflammation, compared to CRSwNP-alone endotype. Furthermore, in accordance with the findings of Tomassen and colleagues  that severe type 2 CRSwNP increased asthma prevalence, our study has also demonstrated that in addition to more frequent atopy, more local eosinophil infiltration, and higher level of tissue IgE, the chance of recurrence of CRSwNP was also increased in patients with CRSwNP + AS compared to patients with CRSwNP-alone. Indeed, we have previously shown that high proportion of eosinophils in nasal tissue act as a reliable prognostic indicator for CRSwNP recurrence . Similarly, a more recent study has indicated that asthma in CRS patients was the only factor that increases the chance of recurrence in CRSwNP patients ; suggesting that comorbid asthma might also be a strong indicator for CRSwNP recurrence.
Airway type 2 inflammation is mainly mediated by eosinophils, mast cells, Th2 cells, ILC2s and IgE-producing B cells, in combination with increased production of cytokines IL-5, IL-4 and IL-13 . Whilst more than 80% of Western white patients with CRSwNP are characterized by type 2 inflammation, less than 50% of CRSwNP cases in East Asian countries show features of type 2 reactions [32, 33]. Like CRSwNP, asthma is also a complex and heterogeneous disease, and two major endotypes of asthma, type 2-high asthma and type 2-low asthma, have been described, based on underlying airway immune-mediated inflammation [34, 35]. In this regard approximately 50% to 60% of all patients with severe asthma in Europe and the United States account for type 2-high asthma, compared to 38.5% of the severe asthmatic patients with type 2-high asthma in China [36,37,38,39]. Thus, whilst it is possible that generally more Western patients with CRSwNP + AS are also likely to be characterised with type 2-high inflammation than Asian patients CRSwNP + AS, this nevertheless needs to be confirmed in well-designed multicenter studies in the future.
Our findings for RNA sequencing provides valuable information for exploring the general molecular mechanisms underlying the pathogenesis of CRSwNP. Similar to the findings of Peng and colleagues , our study has demonstrated that regardless of the subtypes of CRSwNP, genes and pathways that most likely contribute to the pathogenesis of CRSwNP, appear to be mainly associated with cytokine and chemokine signalling pathway, staphylococcus aureus infection, eicosanoid metabolism and cell adhesion molecules. In this respect, the well-known genes or biomarkers closely related to CRSwNP; for example CLC, POSTN, CCL18, IL13, TSLP and BPIFA1; were also identified as the top DE-mRNAs in both CRSwNP-alone and CRSwNP + AS groups.
However, the transcriptome signatures of CRSwNP + AS were characterized by distinct groups of differentially expressed genes and their enriched pathways, compared to CRSwNP-alone. In particular, the present study showed that the CRSwNP + AS endotype was mainly associated with higher activities of arachidonic acid metabolism, Th2 signalling pathway and fibrinolysis pathway, and lower activity of IL-17 signalling pathway. Beyond that, the most significant DE-mRNAs between CRSwNP + AS and CRSwNP-alone also provide important information. In line with the high concentration of tissue IgE in CRSwNP + AS, increased mRNA expression of constant region of heavy chain of IgE (IGHE) was also confirmed by RNA sequencing. Some recently identified biomarkers of CRSwNP, for example CST1 , were also found to be differentially expressed between CRSwNP + AS and CRSwNP-alone. Importantly, the present study showed that some of the most significant DE-mRNAs; including ITLN1, KCNA3 and CCR10; which were up-regulated in nasal tissue of CRSwNP + AS compared to nasal tissue of CRSwNP-alone (Additional file 1: Figure S8), are also expressed in the bronchial tissue and contribute to the pathogenesis of asthma [42,43,44]. This suggests that signature genes identified in nasal tissue of CRSwNP + AS might also play important roles in the pathogenesis of asthma. Moreover, consistent with the united airways concept, it is tempting to speculate that the pattern of gene expression in the upper airway may be influenced by the prevailing conditions in the lower airway, and as such prevalence of comorbid asthmatic conditions might exacerbate type 2 inflammation in nasal tissue of CRSwNP patients.
It has been well demonstrated that alterations in the arachidonic acid pathway play an important role in airway inflammatory conditions like rhinosinusitis, nasal polyps, allergic rhinitis, and asthma . We found a generally enhanced activity of arachidonic acid metabolism in CRSwNP + AS, indicated by the up-regulated expression of PLA2 (PLA2G4A, PLA2G4B, PLA2G4D and PLA2G6), which may promotes the release of membrane-bound arachidonic acid . The imbalanced synthesis of eicosanoids characterized by increased synthesis of cysteinyl leukotrienes (CysLTs) is correlated with the inflammatory pattern and severity of the airway inflammation [45, 47]. Consistent with this, we also showed an increased expression of LTC4S, which promotes the biosynthesis of CysLTs and indicates the presence of severe inflammation in the nasal tissues of CRSwNP + AS patients. The increased expression of PTGS1 and PTGIS may accelerate the conversion of arachidonic acid to prostacyclins (PGI2) and prostaglandins (PGD2, PGE2), which are also related to an aggravated airway inflammation and asthmatic conditions [48, 49]. Moreover, the increased expression of ALOX15B in CRSwNP + AS might have an active proinflammatory role . Collectively, these findings suggest that enhanced activity of arachidonic acid metabolism might help to augment type 2 inflammation in CRSwNP + AS [45, 51].
Th2 and Th17 signalling pathways are two major regulatory pathways, which are inversely interrelated in patients with CRSwNP . In accordanc ewith this, general observation, our RNA sequencing data also revealed higher activities of Th2 signalling pathway and lower activities of IL-17 signalling pathway in CRSwNP + AS than in CRSwNP-alone. In this respect IL-4, IL-5 and IL-13 are typical Th2 cytokines, which reflect the severity of type 2 inflammation, whereas IL-17 is a typical Th17 cytokine. Indeed, the present study further showed that the expression of IL5 and IL13 was positively correlated with other indicators of type 2 inflammation (i.e. increased local eosinophils and tissue IgE). Similarly, genes that are related to arachidonic acid metabolism were also correlated with the expression of type 2 cytokines, thus, emphasising the crucial role of type 2 inflammation in CRSwNP + AS patients.
LncRNAs play important roles in various biological processes and are emerging as reliable biomarkers and potential therapeutic targets of human chronic diseases [52, 53]. However, comparatively few studies have reported the involvement of lncRNAs in chronic nasal inflammation. Yue and colleagues  found that linc00632 was down-regulated in nasal tissues of allergic rhinitis patients and inhibited IL-13 induced inflammatory cytokine and mucus production. Wang and colleagues  showed that lncRNA XLOC_010280 might regulate the expression of CCL18 and eosinophilic inflammation in eosinophilic CRSwNP. In the present study, whole-transcriptome sequencing has revealed several dysregulated lncRNAs in both subtypes of CRSwNP, which will provide a useful pool of candidate lncRNAs in future investigations of CRSwNP endotypes.
WGCNA is an effective method of multigene analysis to construct coexpression network, and has successfully been applied for studying mRNAs and lncRNAs to distinguish dysfunctional regulatory subnetworks, select out potentially key genes, and predict lncRNA functions [56, 57]. We applied WGCNA to predict functions of DE-lncRNAs and identify hub genes in the present study and showed that common dysregulated lncRNAs in CRSwNP-alone and CRSwNP + AS have very similar predictive functions to common DE-mRNAs, as indicated by genes in the maximal module. Furthermore, we identified LINC01146 as the only top hub lncRNA, which might play a key role in the pathogenesis of CRSwNP. Although LINC01146 has also been found to be dysregulated in hepatocellular carcinoma, its precise function in this condition is not clear . Further study is needed to validate the function of common dysregulated lncRNAs in CRSwNP.
The present study also identified HK3-006 as the only top hub lncRNA that was expressed in the module that most positively correlated with phenotypic traits of CRSwNP + AS; and may be related to asthma pathway and arachidonic acid metabolism as predicted by its highly coexpressed mRNAs. Similarly, several top hub lncRNAs were identified in the module that most negatively correlated with phenotypic traits of CRSwNP + AS. Among these lncRNAs, we have identified a new lncRNA, LINC686, which is most likely to be associated with IL-17 signalling pathway. Although an increasing number of studies show that lncRNAs are involved in the regulation of cytokine signalling and inflammation [59, 60], our understanding of the functions of lncRNA is just beginning to develop.
The present study has some limitations. First, this study did not investigate any subgroups of CRSwNP + AS and CRSwNP-alone. This may be of importance as both CRSwNP and asthma are heterogeneous diseases, and one recent study has shown that patients with CRSwNP + AS can be grouped into 3 subtypes with distinct inflammatory status and disease severity . However, the present study was designed to directly compare CRSwNP + AS with CRSwNP-alone. Furthermore, as a substantial number of Asian CRSwNP patients are type 1 and/or type 3 inclined, it would be more appropriate to compare type 2 CRSwNP endotype with type 2 CRSwNP + AS endotype for a better understanding of the mechanisms and pathways operating in either endotype. Second, as the sample size used for RNA sequencing is relatively small, this might decrease the statistical power for gene signature profiling. Third, we did not evaluate lower airway inflammation by assessment of inflammatory cells and cytokines in bronchial biopsy or induced sputum, which might be helpful in understanding the association between CRSwNP and asthma. Fourth, assessment of critical cytokines and their receptors was based only on RNA expression and need to be verified according to the levels of corresponding proteins in the tissue. Fifth, whole tissue samples with multiple cell types were used for RNA sequencing and no single-cell RNA sequencing was used. As single-cell RNA sequencing allows the analysis of transcriptome from individual cells in nasal tissues , this might have provided more specific information on the role of different cell types as biomarkers of CRSwNP with or without comorbid asthma. The use of inferior turbinate tissue from controls for comparison with nasal polyp tissues from CRSwNP patients is also as a limitation of the study, particularly as nasal polyps rarely arise from inferior turbinates in healthy control subjects. However, using mucosa tissue from different locations such as middle turbinate and ethmoid tissue as controls may avoid regional variations in gene expression.
The present study demonstrated that CRSwNP patients with comorbid asthma have distinct type 2-high inflammation-associated transcriptome profiles, as indicated by differential expression profiles of key mRNAs and lncRNAs in nasal tissue, compared to patients with CRSwNP-alone. Differences in molecular mechanisms and type 2 inflammation-related molecules in CRSwNP with comorbid asthma may be useful in better understanding the mechanisms underlying the development of different endotypes of CRS, as well as development of potential biomarkers and targeted gene therapies for CRSwNP in the future.
Availability of data and materials
The complete dataset is included in this manuscript.
Chronic rhinosinusitis with nasal polyps
- CRSwNP + AS:
Chronic Rhinosinusitis with Nasal Polyps with Comorbid Asthma
Long non-coding RNA
Fragments per kilo-base of exon per million fragments mapped
Weighted gene coexpression network analysis
Kyoto Encyclopedia of Genes and Genomes
Fokkens WJ, Lund VJ, Hopkins C, Hellings PW, Kern R, Reitsma S, et al. European position paper on rhinosinusitis and nasal polyps 2020. Rhinology. 2020;58(Suppl S29):1–464.
Stevens WW, Peters AT, Hirsch AG, Nordberg CM, Schwartz BS, Mercer DG, et al. Clinical characteristics of patients with chronic rhinosinusitis with nasal polyps, asthma, and aspirin-exacerbated respiratory disease. J Allergy Clin Immunol Pract. 2017;5(4):1061–70.
Hopkins C. Chronic Rhinosinusitis with nasal polyps. N Engl J Med. 2019;381(1):55–63.
Bhattacharyya N, Villeneuve S, Joish VN, Amand C, Mannent L, Amin N, et al. Cost burden and resource utilization in patients with chronic rhinosinusitis and nasal polyps. Laryngoscope. 2019;129(9):1969–75.
Bachert C, Han JK, Desrosiers M, Hellings PW, Amin N, Lee SE, et al. Efficacy and safety of dupilumab in patients with severe chronic rhinosinusitis with nasal polyps (LIBERTY NP SINUS-24 and LIBERTY NP SINUS-52): results from two multicentre, randomised, double-blind, placebo-controlled, parallel-group phase 3 trials. Lancet. 2019;394:1638–50.
Sella GCP, Tamashiro E, Sella JA, Aragon DC, Mendonca TN, Arruda LKP, et al. Asthma Is the dominant factor for recurrence in chronic rhinosinusitis. J Allergy Clin Immunol Pract. 2020;8(1):302–9.
Cho SW, Kim DW, Kim JW, Lee CH, Rhee CS. Classification of chronic rhinosinusitis according to a nasal polyp and tissue eosinophilia: limitation of current classification system for Asian population. Asia Pac Allergy. 2017;7(3):121–30.
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.
Promsopa C, Kansara S, Citardi MJ, Fakhri S, Porter P, Luong A. Prevalence of confirmed asthma varies in chronic rhinosinusitis subtypes. Int Forum Allergy Rhinol. 2016;6(4):373–7.
Boulet LP. Influence of comorbid conditions on asthma. Eur Respir J. 2009;33(4):897–906.
Samitas K, Carter A, Kariyawasam HH, Xanthou G. Upper and lower airway remodelling mechanisms in asthma, allergic rhinitis and chronic rhinosinusitis: the one airway concept revisited. Allergy. 2018;73(5):993–1002.
Langdon C, Mullol J. Nasal polyps in patients with asthma: prevalence, impact, and management challenges. J Asthma Allergy. 2016;9:45–53.
Wang M, Tan G, Eljaszewicz A, Meng Y, Wawrzyniak P, Acharya S, et al. Laundry detergents and detergent residue after rinsing directly disrupt tight junction barrier integrity in human bronchial epithelial cells. J Allergy Clin Immunol. 2019;143(5):1892–903.
Wang C, Zhang L. Use of biologics in chronic sinusitis with nasal polyps. Curr Opin Allergy Clin Immunol. 2019;19(4):365–72.
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.
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.
Wang M, Zhang N, Zheng M, Li Y, Meng L, Ruan Y, et al. Cross-talk between TH2 and TH17 pathways in patients with chronic rhinosinusitis with nasal polyps. J Allergy Clin Immunol. 2019;144(5):1254–64.
Rix I, Hakansson K, Larsen CG, Frendo M, von Buchwald C. Management of chronic rhinosinusitis with nasal polyps and coexisting asthma: a systematic review. Am J Rhinol Allergy. 2015;29(3):193–201.
Wu D, Bleier BS, Li L, Zhan X, Zhang L, Lv Q, et al. Clinical phenotypes of nasal polyps and comorbid asthma based on cluster analysis of disease history. J Allergy Clin Immunol Pract. 2018;6(4):1297–305.
Hakansson K, Thomsen SF, Konge L, Mortensen J, Backer V, von Buchwald C. A comparative and descriptive study of asthma in chronic rhinosinusitis with nasal polyps. Am J Rhinol Allergy. 2014;28(5):383–7.
Qi S, Yan B, Liu C, Wang C, Zhang L. Predictive significance of Charcot-Leyden Crystal mRNA levels in nasal brushing for nasal polyp recurrence. Rhinology. 2020;58(2):166–74.
Li J, Jiao J, Wang M, Gao Y, Li Y, Wang Y, et al. Hypomethylation of the IL8 promoter in nasal epithelial cells of patients with chronic rhinosinusitis with nasal polyps. J Allergy Clin Immunol. 2019;144(4):993–1003.
Chen EY, Tan CM, Kou Y, Duan Q, Wang Z, Meirelles GV, et al. Enrichr: interactive and collaborative HTML5 gene list enrichment analysis tool. BMC Bioinformat. 2013;14:128.
Langfelder P, Horvath S. WGCNA: an R package for weighted correlation network analysis. BMC Bioinformat. 2008;9:559.
Stevens WW, Peters AT, Tan BK, Klingler AI, Poposki JA, Hulse KE, et al. Associations between inflammatory endotypes and clinical presentations in chronic rhinosinusitis. J Allergy Clin Immunol Pract. 2019;7(8):2812–20.
De Greve G, Hellings PW, Fokkens WJ, Pugin B, Steelant B, Seys SF. Endotype-driven treatment in chronic upper airway diseases. Clin Transl Allergy. 2017;7:22.
Tan BK, Klingler AI, Poposki JA, Stevens WW, Peters AT, Suh LA, et al. Heterogeneous inflammatory patterns in chronic rhinosinusitis without nasal polyps in Chicago. Illinois. J Allergy Clin Immunol. 2017;139(2):699–703.
Avdeeva K, Fokkens W. Precision medicine in chronic rhinosinusitis with nasal polyps. Curr Allergy Asthma Rep. 2018;18(4):25.
Bachert C, Zhang N, Hellings PW, Bousquet J. Endotype-driven care pathways in patients with chronic rhinosinusitis. J Allergy Clin Immunol. 2018;141(5):1543–51.
Lou H, Meng Y, Piao Y, Wang C, Zhang L, Bachert C. Predictive significance of tissue eosinophilia for nasal polyp recurrence in the Chinese population. Am J Rhinol Allergy. 2015;29(5):350–6.
Fahy JV. Type 2 inflammation in asthma–present in most, absent in many. Nat Rev Immunol. 2015;15(1):57–65.
Wang ET, Zheng Y, Liu PF, Guo LJ. Eosinophilic chronic rhinosinusitis in East Asians. World J Clin Cases. 2014;2(12):873–82.
Bachert C, Zhang L, Gevaert P. Current and future treatment options for adult chronic rhinosinusitis: focus on nasal polyposis. J Allergy Clin Immunol. 2015;136(6):1431–40.
Woodruff PG, Modrek B, Choy DF, Jia G, Abbas AR, Ellwanger A, et al. T-helper type 2-driven inflammation defines major subphenotypes of asthma. Am J Resp Crit Care Med. 2009;180(5):388–95.
Licari A, Castagnoli R, Brambilla I, Marseglia A, Tosca MA, Marseglia GL, et al. Asthma endotyping and biomarkers in childhood asthma. Pediatr Allergy Immunol Pulmanol. 2018;31(2):44–55.
Svenningsen S, Nair P. Asthma endotypes and an overview of targeted therapy for asthma. Front Med. 2017;4:158.
Aleman F, Lim HF, Nair P. Eosinophilic endotype of asthma. Immunol Allergy Clin North Am. 2016;36(3):559–68.
Qiu R, Xie J, Chung KF, Li N, Yang Z, He M, et al. Asthma phenotypes defined from parameters obtained during recovery from a hospital-treated exacerbation. J Allergy Clin Immunol Pract. 2018;6(6):1960–7.
Xie M, Wenzel SE. A global perspective in asthma: from phenotype to endotype. Chin Med J. 2013;126(1):166–74.
Peng Y, Zi XX, Tian TF, Lee B, Lum J, Tang SA, et al. Whole-transcriptome sequencing reveals heightened inflammation and defective host defence responses in chronic rhinosinusitis with nasal polyps. Eur Respir J. 2019;54(5):1900732.
Yan B, Lou H, Wang Y, Li Y, Meng Y, Qi S, et al. Epithelium-derived cystatin SN enhances eosinophil activation and infiltration through IL-5 in patients with chronic rhinosinusitis with nasal polyps. J Allergy Clin Immunol. 2019;144(2):455–69.
Koshy S, Huq R, Tanner MR, Atik MA, Porter PC, Khan FS, et al. Blocking KV1.3 channels inhibits Th2 lymphocyte function and treats a rat model of asthma. J Biol Chem. 2014;289(18):12623–32.
Yi L, Cheng D, Zhang K, Huo X, Mo Y, Shi H, et al. Intelectin contributes to allergen-induced IL-25, IL-33, and TSLP expression and type 2 response in asthma and atopic dermatitis. Mucosal Immunol. 2017;10(6):1491–503.
Daubeuf F, Jung F, Douglas GJ, Chevalier E, Frossard N. Protective effect of a Protein Epitope Mimetic CCR10 antagonist, POL7085, in a model of allergic eosinophilic airway inflammation. Respir Res. 2015;16:77.
Pezato R, Swierczynska-Krepa M, Nizankowska-Mogilnicka E, Derycke L, Bachert C, Perez-Novo CA. Role of imbalance of eicosanoid pathways and staphylococcal superantigens in chronic rhinosinusitis. Allergy. 2012;67(11):1347–56.
Harizi H, Corcuff JB, Gualde N. Arachidonic-acid-derived eicosanoids: roles in biology and immunopathology. Trends Mol Med. 2008;14(10):461–9.
Miyata J, Fukunaga K, Kawashima Y, Ohara O, Arita M. Cysteinyl leukotriene metabolism of human eosinophils in allergic disease. Allergol Int. 2020;69(1):28–34.
Ishiura Y, Fujimura M, Nobata K, Oribe Y, Abo M, Myou S. Prostaglandin I2 enhances cough reflex sensitivity to capsaicin in the asthmatic airway. Cough. 2007;3:2.
Harrington LS, Lucas R, McMaster SK, Moreno L, Scadding G, Warner TD, et al. COX-1, and not COX-2 activity, regulates airway function: relevance to aspirin-sensitive asthma. FASEB J. 2008;22(11):4005–10.
Magnusson LU, Lundqvist A, Karlsson MN, Skalen K, Levin M, Wiklund O, et al. Arachidonate 15-lipoxygenase type B knockdown leads to reduced lipid accumulation and inflammation in atherosclerosis. PLoS ONE. 2012;7(8):e43142.
Samuchiwal SK, Boyce JA. Role of lipid mediators and control of lymphocyte responses in type 2 immunopathology. J Allergy Clin Immunol. 2018;141(4):1182–90.
Devadoss D, Long C, Langley RJ, Manevski M, Nair M, Campos MA, et al. Long noncoding transcriptome in chronic obstructive pulmonary disease. Am J Respir Cell Mol Biol. 2019;61(6):678–88.
Schulte LN, Bertrams W, Stielow C, Schmeck B. NcRNAs in inflammatory and infectious diseases. Methods Mol Biol. 2019;1912:3–32.
Yue L, Yin X, Hao F, Dong J, Ren X, Xu O, et al. Long Noncoding RNA Linc00632 Inhibits Interleukin-13-Induced inflammatory cytokine and mucus production in nasal epithelial cells. J Innate Immun. 2020;12(1):116–28.
Wang W, Gao Z, Wang H, Li T, He W, Lv W, et al. Transcriptome analysis reveals distinct gene expression profiles in eosinophilic and noneosinophilic chronic rhinosinusitis with nasal polyps. Sci Rep. 2016;6:26604.
Liu H, Zhou T, Wang B, Li L, Ye D, Yu S. Identification and functional analysis of a potential key lncRNA involved in fat loss of cancer cachexia. J Cell Biochem. 2018;119(2):1679–88.
Duan H, Ge W, Zhang A, Xi Y, Chen Z, Luo D, et al. Transcriptome analyses reveal molecular mechanisms underlying functional recovery after spinal cord injury. Proc Natl Acad Sci USA. 2015;112(43):13360–5.
Esposti DD, Hernandez-Vargas H, Voegele C, Fernandez-Jimenez N, Forey N, Bancel B, et al. Identification of novel long non-coding RNAs deregulated in hepatocellular carcinoma using RNA-sequencing. Oncotarget. 2016;7(22):31862–77.
Carpenter S, Fitzgerald KA. Cytokines and long noncoding RNAs. Cold Spring Harb Perspect Biol. 2018;10(6):a028589.
Marques-Rocha JL, Samblas M, Milagro FI, Bressan J, Martinez JA, Marti A. Noncoding RNAs, cytokines, and inflammation-related diseases. FASEB J. 2015;29(9):3595–611.
Ordovas-Montanes J, Dwyer DF, Nyquist SK, Buchheit KM, Vukovic M, Deb C, et al. Allergic inflammatory memory in human respiratory epithelial progenitor cells. Nature. 2018;560(7720):649–54.
We thank Professor Harald Renz and Professor Diethard Gemsa for discussing the data.
This work was supported by grants from the national key R&D program of China (2018YFC0116800, 2016YFC0905200), the program for the Changjiang scholars and innovative research team (IRT13082), the national natural science foundation of China (81800882, 81630023, 81870698), Beijing municipal administration of hospitals’ mission plan (SML20150203), Beijing municipal administration of hospitals’ Dengfeng plan (DFL20190202), Beijing municipal administration of hospitals clinical medicine development of special funding support (XMLX201816) and capital health development foundation (2016-1-2052).
Ethics approval and consent to participate
The Ethics Committee of Beijing Tongren Hospital approved this study, and all subjects signed informed consent forms.
The authors declare that they have no conflicts of interest.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Demographic and clinical characteristics of study subjects. Table S2. Common DE-mRNAs and DE-lncRNAs shared by CRSwNP-alone versus control and CRSwNP+AS versus control. Table S3. Top 50 DE-mRNAs of CRSwNP+AS versus CRSwNP-alone. Table S4. Expression of key cytokines and their receptors in nasal tissues from control subjects and CRSwNP-alone and CRSwNP+AS patients. Table S5. Top 50 DE-lncRNAs of CRSwNP+AS versus CRSwNP-alone. Figure S1. Correlation of infiltrating eosinophils and total IgE in nasal tissues of patients with CRSwNP. Figure S2. Differentially expressed mRNAs and differentially expressed lncRNAs in nasal tissues of CRSwNP patients. Figure S3. Hierarchical clustering of differentially expressed genes. Figure S4. GO biological processes enriched by common dysregulated genes in CRSwNP+AS and CRSwNP-alone. FigureS5. Gene number of modules identified by weighted gene co-expression network analysis. Figure S6. The expression and potential functions of LINC01146. Figure S7. Gene modules identified by WGCNA based on expression of DE-mRNAs and DE-lncRNAs in CRSwNP+AS versus CRSwNP-alone. Figure S8. Asthma related genes expressed in nasal tissue of patients with CRSwNP+AS and patients with CRSwNP-alone.
About this article
Cite this article
Wang, M., Bu, X., Luan, G. et al. Distinct type 2-high inflammation associated molecular signatures of chronic rhinosinusitis with nasal polyps with comorbid asthma. Clin Transl Allergy 10, 26 (2020). https://doi.org/10.1186/s13601-020-00332-z
- Chronic rhinosinusitis with nasal polyps
- Type 2 inflammation
- Molecular endotype
- Transcriptome sequencing