Abstract
Introduction The pathways underlying chronic rhinosinusitis with nasal polyps (CRSwNP) are unclear. We conducted genome-wide gene expression analysis to determine pathways and candidate gene sets associated with CRSwNP.
Methods We performed whole-transcriptome RNA sequencing on 42 polyp (CRSwNP-NP) and 33 paired nonpolyp inferior turbinate (CRSwNP-IT) tissues from patients with CRSwNP and 28 inferior turbinate samples from non-CRS controls (CS-IT). We analysed the differentially expressed genes (DEGs) and the gene sets that were enriched in functional pathways.
Results Principal component-informed analysis revealed cilium function and immune regulation as the two main Gene Ontology (GO) categories differentiating CRSwNP patients from controls. We detected 6182 and 1592 DEGs between CRSwNP-NP versus CS-IT and between CRSwNP-NP versus CRSwNP-IT tissues, respectively. Atopy status did not have a major impact on gene expression in various tissues. GO analysis on these DEGs implicated extracellular matrix (ECM) disassembly, O-glycan processing, angiogenesis and host viral response in CRSwNP pathogenesis. Ingenuity Pathway Analysis identified significant enrichment of type 1 interferon signalling and axonal guidance canonical pathways, angiogenesis, and collagen and fibrotic changes in CRSwNP (CRSwNP-NP and CRSwNP-IT) tissues compared with CS-IT. Finally, gene set enrichment analysis implicated sets of genes co-regulated in processes associated with inflammatory response and aberrant cell differentiation in polyp formation.
Conclusions Gene signatures involved in defective host defences (including cilia dysfunction and immune dysregulation), inflammation and abnormal metabolism of ECM are implicated in CRSwNP. Functional validation of these gene expression patterns will open opportunities for CRSwNP therapeutic interventions such as biologics and immunomodulators.
Abstract
Whole-transcriptome sequencing shows gene signatures associated with interferon signalling response and host defence are implicated in chronic rhinosinusitis with nasal polyps, opening up new opportunities for therapeutic interventions http://bit.ly/2KyF3dT
Introduction
Chronic rhinosinusitis (CRS) is a debilitating inflammatory disease of the nasal mucosa [1]. Nasal polyps (NP) are the most common comorbidity, affecting ∼30% of patients with CRS [2], and cause a considerable socioeconomic burden and impaired quality of life [3]. Clinical management of CRS with NP (CRSwNP) is largely ineffective, partly due to the limited understanding of the underlying pathogenic factors and a paucity of effective therapeutic interventions [4, 5]. Understanding the mechanisms that underlie CRSwNP pathogenesis may help identify targets in critical candidate pathways for therapeutic interventions.
The known biological processes implicated in CRSwNP are multifaceted, including airway inflammation [2], tight junction impairment [6], pathogen infection [7, 8] and a defective host defence [9–13]. The complexity of gene expression patterns in CRSwNP has not been adequately addressed, partly because the reference tissues often used have been suboptimal for biomarker profiling. For instance, differences in gene signatures between NP and control tissues, i.e. inferior turbinate (IT), from healthy subjects might stem from genetic or environmental factors, or both [14–17]. Additionally, Gene Ontology (GO) analyses have mostly been based on differentially expressed genes (DEGs) identified by microarray assays, but not on gene sets or pathways [18]. Finally, the limited sample sizes have markedly decreased the statistical power for gene signature profiling [14, 15]. Ordovas-Montanes et al. [19], however, reported strong type 2 immune-mediated signatures associated with barrier tissue dysfunction in NP, despite a relatively small sample size that did not compromise the analysis because of deep single-cell profiling. Whole-genome approaches that use well-powered sample sizes are necessary to thoroughly elucidate gene signatures in CRSwNP.
RNA sequencing is a powerful means to explore candidate transcripts in an unbiased manner [20–22], and to identify cellular and molecular pathways. We hypothesised that the gene signatures and pathways identified by RNA sequencing would differ considerably between pairwise comparisons of the IT of patients with CRSwNP and control subjects. In this study, we thus performed RNA sequencing to identify the transcriptomic signatures associated with CRSwNP, by incorporating a three-way comparison of the nasal tissues from CRSwNP patients and controls as well as a comprehensive analytical pipeline to determine the DEGs and gene sets that enriched in functional pathways. The gene signatures and candidate pathways would help identify disease targets (i.e. abnormal matrix metabolism, immune dysfunction) and novel treatment options such as biologics for CRSwNP [14–16, 18, 20].
Methods
Study design and participants
Subjects were recruited from individual sites and provided written informed consent. The Institutional Review Boards of Guangdong General Hospital (Guangzhou, China), The Second Hospital of Shandong University (Jinan, China), The First Affiliated Hospital of Nanchang University (Nanchang, China) and the National Healthcare Group Domain Specific Review Board of Singapore (Singapore) approved the study protocol.
CRSwNP diagnosis was based on the European Position Paper on Rhinosinusitis and Nasal Polyps 2012 [1]. All CRSwNP patients had bilateral NP (scored 1–4 points bilaterally based on the NP size score, evaluated via nasal endoscopy) [23, 24]. Patients with antrochoanal polyps, fungal sinusitis and/or recurrent lower airway infections were excluded. None of the participants had CRS without NP or aspirin-exacerbated respiratory diseases. We recruited control subjects who underwent septoplasty for anatomical variations without otolaryngologist-diagnosed CRS (free from clinical symptoms; confirmed via inspection of nasal endoscopy). Neither atopy nor allergic rhinitis were mandatory for the exclusion criteria for both CRSwNP patients and control subjects. Furthermore, none of the study participants had an upper respiratory tract infection or had used intranasal or systemic corticosteroids or antibiotics within 4 weeks before the surgery. In total, we recruited 44 CRSwNP patients and 41 control subjects. All but three study participants (Singaporeans) were of Chinese ethnicity. Data from 42 CRSwNP-NP samples, 33 CRSwNP-IT samples and 28 CS-IT samples passed the initial quality control of RNA sequencing and were therefore included in further analysis.
In this cross-sectional study (figure 1), paired IT and NP was biopsied from patients with CRSwNP, and IT was biopsied from control subjects. Following biopsy, tissues were immediately split for RNA sequencing and histological assessment. Haematoxylin/eosin staining and immunohistochemistry were performed as previously described for eosinophil and neutrophil enumeration of paraffin-embedded slides, respectively [25]. Eosinophilic or neutrophilic NP denoted the percentage of eosinophils or neutrophils >10%, respectively [25, 26].
Study workflow from subject recruitment and sample collection to data analysis and interpretation. CRSwNP-NP: nasal polyps tissues from patients with chronic rhinosinusitis with nasal polyps; CS-IT: inferior turbinate in control subjects; CRSwNP-IT: paired inferior turbinate tissues from patients with chronic rhinosinusitis with nasal polyps. More information on the detailed methods is provided in the Methods section.
RNA sequencing and data processing
Tissues were pre-treated with RNAlater (Thermo Fisher Scientific, Waltham, MA, USA) for overnight incubation at 4°C and then stored in −80°C freezers before RNA extraction. Total RNA was extracted by a double extraction protocol: RNA isolation by acid guanidinium thiocyanate–phenol–chloroform extraction (TRIzol; Thermo Fisher Scientific), followed by RNeasy Micro clean-up (Qiagen, Hilden, Germany). All human RNAs were analysed on an Agilent Bioanalyzer (Agilent, Santa Clara, CA, USA) for quality assessment by RNA integrity number (range 5.8–9.1; median 7.3). cDNA libraries were prepared with 2 ng total RNA input and 1 μL of a 1:50 000 dilution of ERCC RNA Spike-in Controls (Thermo Fisher Scientific) using the SMARTSeq version 2 protocol with the following modifications: 1) 20 µM template switch oligonucleotides and 2) 250 pg cDNA with 1/5 reaction of a Nextera XT kit (Illumina, San Diego, CA, USA). The length distribution of the cDNA libraries was monitored using a DNA High Sensitivity Reagent Kit on a Labchip GX system (Perkin Elmer, Waltham, MA, USA). All samples were subjected to an indexed paired-end sequencing run of 2×151 cycles on an HiSeq 4000 platform (Illumina) (30 samples per lane).
Statistical analyses
We expected the gene signatures identified through whole-transcriptome sequencing to be independent of the CRSwNP phenotypes. To eliminate interindividual differences arising from genetic or environmental factors, we conducted a three-way comparison between the IT in control subjects (CS-IT) versus the paired NP (CRSwNP-NP) and IT (CRSwNP-IT) in patients with CRSwNP. To minimise the heterogeneity of the study population, we recruited participants without physician-diagnosed or self-reported asthma. We stratified patients with CRSwNP by different inflammatory phenotypes (i.e. eosinophilic or neutrophilic NP phenotypes) to determine whether the inflammatory phenotypes markedly confound gene expression patterns. We recruited an adequate sample size of participants (n>25 for each group) to ensure sufficient statistical power to investigate the interactions of different gene sets and pathways. Importantly, we designed an analytical pipeline to identify different gene signatures between CRSwNP-IT and CS-IT within the same individuals with CRSwNP. These novel study designs have collectively underpinned the validity of our findings.
HiSeq FASTQ files were mapped to the human genome build hg38 using STAR (splice junction annotations based on GENCODE version 26; www.gencodegenes.org). Gene counts were obtained with featureCounts using GENCODE version 26 annotations. Principal component analysis (PCA) was carried out in R version 3.3.3 (www.r-project.org) using log2 reads per kilobase million (RPKM) to determine the differences between the samples in each group. DEG analysis was performed using edgeR as pairwise comparisons on genes filtered for log2 RPKM interquartile range >0.5. We presented the adjusted p-values to derive more conservative significance thresholds given the large number of tests performed. We corrected for multiple testing with the Benjamini–Hochberg algorithm (adjusted p-value <0.05 was considered statistically significant). Venn diagrams were generated using R software. GO enrichment analysis using the Bioconductor (www.bioconductor.org) package TopGO was performed to identify biological processes and protein functional groups enriched in groups of DEGs (background was the filtered list of genes used for DEG analysis). Ingenuity Pathway Analysis (www.nihlibrary.nih.gov/resources/tools/ingenuity-pathways-analysis-ipa) was used to define gene sets-determined pathways associated with DEGs. Significance threshold was set to default. Gene Set Enrichment Analysis (GSEA) was conducted using the Bioconductor package fgsea to identify gene sets of enrichment patterns using the Molecular Signatures Database (MSigDB; http://software.broadinstitute.org/gsea/msigdb). All software was run using default parameters, unless otherwise stated.
Results
Subject recruitment and clinical characteristics
We recruited 44 CRSwNP patients and 41 control subjects whose surgically resected samples were subject to RNA sequencing. Overall, CRSwNP patients were older than control subjects and had a lower proportion of atopy (both p<0.05) (supplementary table E1). 15 patients with CRSwNP (34.1%) had tissue eosinophilia. The mean±sd NP size score was 4.7±1.1. None of the study participants had physician-diagnosed or self-reported asthma. Subanalysis of the 42 CRSwNP patients and 28 control subjects whose RNA sequencing data passed the initial quality control is also shown in supplementary table E1, which was not materially altered compared with that of the overall study population.
Global transcriptomic profiling reveals that ciliary function and immunomodulation segregate the nasal tissues
We first evaluated the overall distribution of gene expression across the three tissues types by PCA on the 103 tissue samples that passed quality control standards (figure 2). No technical effects in PCA were identified. Confounders such as batch effects, RNA integrity number and sample collection processes did not confer any major influence on the analysis. However, the global expression patterns revealed a spatial separation on the PCA matrix based on tissue category, with the tissues from the CRSwNP patients separating from the control subject tissues. The strongest variance was observed between CRSwNP-NP and CS-IT, and a comparatively minor separation was observed between the nonpolyp CRSwNP-IT and CS-IT tissues (figure 2a–c).
Principal component (PC)-informed analysis of the global gene expression patterns. CRSwNP-NP: nasal polyps tissues from patients with chronic rhinosinusitis with nasal polyps; CRSwNP-IT: paired inferior turbinate tissues from patients with chronic rhinosinusitis with nasal polyps; CS-IT: inferior turbinate in control subjects; PLC: phospholipase C; GPCR: G-protein coupled receptor; IFN: interferon; ECM: extracellular matrix. a) PC plot with sample clustering using the variance of axis 1 (PC1) and 2 (PC2). Box plots of b) PC1 and c) PC2 stratified by tissues: CRSwNP-NP, CRSwNP-IT and CS-IT. Gene Ontology of the biological processes associated with the top loading genes from d) PC1 and e) PC2.
We then analysed the top loading genes accounting for this observed variance (supplementary table E2). The first principal component (PC1) comprised genes regulating ciliogenesis and ciliary function (i.e. cilium assembly, ultrastructural protein assembly, cilium motility, intraciliary transport), whereas PC2 mainly consisted of genes modulating inflammatory and immune responses (i.e. chemokine-mediated signalling, innate and adaptive immune response, T-cell proliferation, co-stimulation), and collagen and extracellular matrix (ECM) metabolism (i.e. complex of collagen trimmers, ECM structural constituent and organisation) (figure 2d and e). We then performed a GO analysis to further evaluate these identified genes using the known evidence from previously reported associations.
Next, we performed clustering analysis on the global gene expression data to further examine the heterogeneity of gene expression among individual participants independent of their disease status. Here, we observed tighter clustering between CRSwNP-NP and CRSwNP-IT subjects than with CS-IT (supplementary figure E1). This analysis validated the PCA findings, despite that CRSwNP-IT shared certain gene expression patterns with CRSwNP-NP.
Collectively, these findings highlight a more pronounced difference in gene signatures (particularly cilia-associated and inflammatory genes) in CRSwNP patients independent of the tissue type. This is critical since the inflammation in CRSwNP-IT could also contribute to disease.
Pairwise comparisons identify more DEGs in NP than in control IT tissue
Having delineated the global gene expression patterns among each tissue type, we next evaluated the DEGs for pairwise comparison among the tissues. For this purpose, we performed edgeR analysis on the quality controlled sequencing data using transcript-level counts. We first sought to evaluate the differences among nonpolyp IT tissues. Here, we only identified 915 DEGs for the CRSwNP-IT and CS-IT comparison that survived multiple testing (supplementary table E3), which corroborated our expectations from the PC-informed analysis. The top GO sets enriched by these DEGs were involved in the interferon (IFN) signalling pathway and viral responses (table 1 and supplementary table E4). By contrast, we identified a notably larger number of DEGs (n=6182) when comparing CRSwNP-NP with CS-IT. These DEGs were involved in collagen processing and organisation, inflammatory responses, and O-glycan processing (table 1 and supplementary table E4). We then made a comparison between tissues from CRSwNP patients only, which revealed 1592 DEGs. Interestingly, the GO analysis identified a broader set of significantly enriched pathways from these DEGs (table 1 and supplementary table E4) that ranged from collagen organisation, angiogenesis, ECM disassembly to inflammatory responses.
Gene Ontology (GO) analysis for the biological processes associated with differentially expressed genes for pairwise analysis of nasal epithelial tissues
We next thoroughly examined the top DEGs (based on log fold change) for the pairwise comparisons (figure 3a–c). Interestingly, the DEGs in this analysis identified biological processes that were also identified in the previous analyses, such as collagen organisation and inflammatory responses (particularly IFN signalling). We additionally observed other significantly enriched biological processes, such as ECM disassembly, angiogenesis, ossification and neutrophil chemotaxis, that are reportedly involved in CRSwNP pathogenesis. Furthermore, we identified certain tissue-specific gene signatures such as CCL18 levels that were enriched in CRSwNP-NP compared with CRSwNP-IT and CS-IT (log fold change 6.48 and 6.85, respectively; both false discovery rate <0.05) (figure 3d–f). The top DEGs for pairwise comparisons are detailed in supplementary tables E5–E7 and all DEGs for pairwise comparisons are listed in supplementary table E8.
Differentially expressed genes (DEGs) stratified by tissue type. CRSwNP-NP: nasal polyps from patients with chronic rhinosinusitis with nasal polyps; CS-IT: inferior turbinate in control subjects; CRSwNP-IT: paired inferior turbinate tissues from patients with chronic rhinosinusitis with nasal polyps. Volcano plots showing the DEGs for a) CRSwNP-NP versus CS-IT, b) CRSwNP-NP versus CRSwNP-IT and c) CRSwNP-IT versus CS-IT. Positive log fold changes are shown in red and negative log fold changes are shown in green. d–f) Heat maps of the top 10 positive log fold changes and the top 10 negative log fold changes for each of the above comparisons.
Surprisingly, further analysis including the covariates indicative of the inflammatory phenotypes (eosinophilia and neutrophilia) did not reveal strong DEG signatures (supplementary table E9). After adjusting for multiple testing, only 16 and two DEGs were detected for eosinophilic and neutrophilic phenotypes, respectively. Analysis accounting for coexisting allergic rhinitis revealed only 35 DEGs, highlighting that allergic rhinitis was not a major confounder for our transcriptomic analysis. Moreover, atopic status conferred limited impact on the gene expression in CRSwNP patients and control subjects (supplementary table E10). However, caution should be exercised when interpreting this covariate analysis given the limited samples size for the stratified samples.
Finally, we compared the DEGs in each pairwise analysis and those that were shared among CRSwNP-IT, CRSwNP-NP and CS-IT (figure 4). The most unique DEGs were observed in the CRSwNP-NP versus CS-IT analysis; here, 4114 (66.54%) DEGs were not detected in other analyses. Only 124 DEGs were shared among the three pairwise analyses and, as expected, the lowest count of unique DEGs (75 (8.19%)) was found between CRSwNP-IT with CS-IT, suggesting strong differential expression in nonpolyp tissues.
Venn diagram representing the overlap between the differentially expressed genes (DEGs) among the three pairwise comparisons: paired inferior turbinate tissues from patients with chronic rhinosinusitis with nasal polyps (CRSwNP-IT) versus inferior turbinate in control subjects (CS-IT), nasal polyps from patients with chronic rhinosinusitis with nasal polyps (CRSwNP-NP) versus CRSwNP-IT and CRSwNP-NP versus CS-IT.
Taken together, these findings revealed strong DEG signals related to inflammation, viral responses and airway remodelling in CRSwNP that were not confounded by the inflammatory phenotypes.
IPA highlights axonal guidance signalling as the common pathway in CRSwNP-NP, CRSwNP-IT and CS-IT
To determine the relevance of the identified DEGs to CRSwNP, we performed IPA on two sets of genes. For the first set, we used the 1387 DEGs shared (124+1263) between the CRSwNP-NP versus CS-IT comparison and the CRSwNP-NP versus CRSwNP-IT comparison to identify enriched pathways. The top canonical pathways (determined by the percentage overlap for between-group comparisons) enriched for these sets of DEGs consisted of agranulocyte adhesion and diapedesis, axonal guidance and interleukin (IL)-8 signalling, and phagosome formation (table 2). For the second set, we analysed the 4795 (681+4114) DEGs detected in the CRSwNP-NP versus CS-IT comparison, but not the CRSwNP-NP versus CRSwNP-IT comparison (table 2). These genes represent canonical pathways including molecular mechanisms of cancer, N-formylmethionyl-leucyl-phenylalanine signalling in neutrophils, and phosphate and tension homology deleted on chromosome 10 and eukaryotic initiation factor-2 signalling. Axonal guidance signalling was a common canonical pathway to both analyses.
The top canonical pathways for the differentially expressed gene (DEGs) identified between nasal polyps from patients with chronic rhinosinusitis with nasal polyps (CRSwNP-NP) versus inferior turbinate in control subjects (CS-IT) and CRSwNP-NP versus paired inferior turbinate tissues from patients with chronic rhinosinusitis with nasal polyps (CRSwNP-IT)
These findings imply that exaggerated inflammatory responses and epithelial growth may predispose NP formation and that axonal guidance signalling represents a critical target for NP management.
Enrichment analysis revealed gene sets associated with the inflammatory response and aberrant cell differentiation potentially involved in NP formation
Finally, we performed an enrichment analysis using GSEA to identify concordant gene expression patterns in CRSwNP-NP compared with CRSwNP-IT and CS-IT. Representative candidate enriched gene sets for CRSwNP-NP versus CS-IT revealed myosin filament, microtubule movement, and other epithelium integrity and function categories (table 3). In addition, we detected strong inflammatory response genes and those involved in inflammasome-associated IL-1 secretion. Interestingly, immune responses, such as immunological memory and monocyte or macrophage immune responses, were highly enriched in the CRSwNP-NP versus CRSwNP-IT comparisons, whereas the candidate gene sets of IFN-associated inflammation and keratinisation of epithelial cells were enriched in the CRSwNP-IT versus CS-IT comparisons (supplementary tables E11 and E12). We observed coherent and consistent gene sets involved in the adaptive immune response that were upregulated in CRSwNP-NP versus CRSwNP-IT (figure 5a), including regulatory T-cell (Treg) molecules such as FOXP3, complement genes such as C4BPA and C4BPB, ILs such as IL4, IL4R and IL1RL1, and protein kinase PTPRC. Another interesting candidate pathway between CRSwNP-NP and CRSwNP-IT involved IFN signatures that were usually observed in cancer (figure 5b). Importantly, genes such as STAT1, IFIT1, IFIT3 and IRF9 were all consistently downregulated in CRSwNP-NP compared with paired IT tissue (CRSwNP-IT), as observed in many cancer tissues. These findings indicate that the inflammation in CRSwNP-NP due to infections might be mediated by IFN-stimulated genes.
Candidate enriched gene sets associated with the comparison between nasal polyps from patients with chronic rhinosinusitis with nasal polyps (CRSwNP-NP) and paired inferior turbinate tissues from patients with chronic rhinosinusitis with nasal polyps (CRSwNP-IT)
Enrichment plots of key pathways identified from the gene set analysis. CRSwNP-NP: nasal polyps from patients with chronic rhinosinusitis with nasal polyps; CRSwNP-IT: paired inferior turbinate tissues from patients with chronic rhinosinusitis with nasal polyps; CS-IT: inferior turbinate in control subjects. a, b) Selected enrichment plots for the pairwise CRSwNP-NP versus CRSwNP-IT comparisons: a) negative regulation of adaptive immune response and b) Einav interferon signature in cancer. Selected pathways for the c) CRSwNP-NP versus CS-IT and d) CRSwNP-NP versus CS-IT comparisons: c) molecular mediators in inflammatory response and d) regulation of tolerance induction.
We found a more generic inflammatory response network identified through GSEA when comparing CRSwNP-NP with CS-IT (figure 5c). However, this inflammatory response could be associated with CRS pathogenesis, NP formation or both. Lastly, when comparing CRSwNP-IT and CS-IT, we found evidence of persistent inflammation mediated through a lack of tolerance induction (figure 5d). Interestingly, key genes involved in immunoregulation and Treg function included FOXP3, IL2RA, IDO1, IRAK3 and HLAG, suggesting strong inflammation in the nonpolyp tissues in patients with CRSwNP, and the lack of tolerance might potentially lead to NP formation. Details of all three comparisons are provided in supplementary table E13.
Taken together, these data indicate that gene sets associated with the inflammatory and immune responses are enriched in NP, and support that inflammation and aberrant cell differentiation underlie early-stage NP formation.
Discussion
We performed whole-transcriptomic sequencing followed by pairwise comparisons of CRSwNP-NP, CRSwNP-IT and CS-IT, which allowed us to identify candidate genes and canonical pathways. The global gene expression patterns differed most dramatically between CRSwNP-NP and CS-IT, whereas CRSwNP-IT shared a certain degree of overlapping gene signatures with CS-IT. Apart from the identification of DEGs associated with ECM disassembly, O-glycan processing and host viral response, we revealed a novel association for NP to axonal guidance signalling as a common pathway and the potential induction of immune intolerance associated with persistent airway inflammation in NP.
CRS pathogenesis with NP formation follows sequential events, including environmental triggers (i.e. allergen sensitisation, Staphylococcus aureus infections), heightened inflammation, angiogenesis and tissue remodelling [1, 27]. Intact epithelial barrier function is essential to protect the airways from allergens or pathogenic insults. Indeed, our PC-informed analysis revealed genes associated with ciliary assembly and function and immune responses, which echoed previous studies regarding aberrant ciliary ultrastructural marker expression, mislocalisation of ciliogenesis markers and aberrant immune function in CRSwNP [28, 29]. Importantly, our GO analysis highlighted the interactions between inflammation (and possibly chronic latent viral infection), defective host defences (cellular immunity and impaired mucociliary clearance) and airway remodelling in driving CRSwNP. Because these signatures could have stemmed from the interaction between impaired mucociliary clearance and tissue-specific inflammation, we made further comparisons between CRSwNP-NP and CRSwNP-IT [2]. These comparisons revealed DEGs involved in collagen and ECM metabolism, inflammatory responses, ossification, and angiogenesis, reaffirming the roles of a heightened IFN signalling response [30], angiogenesis [31], ossification [32], increased matrix metalloproteinase expression [33], and collagen synthesis and deposition [34] in modulating CRSwNP formation.
Pathway analyses revealed canonical pathways that mediate NP formation. The DEGs shared between CRSwNP-NP versus CS-IT and CRSwNP-NP versus CRSwNP-IT comparisons revealed the axonal guidance signalling pathway [35], which reportedly correlates with chronic airway inflammation [36]. Interestingly, the unique DEGs for the aforementioned comparisons revealed that pathways associated with the molecular mechanisms of cancer and neutrophil N-formylmethionyl-leucyl-phenylalanine signalling drive polyp formation. Suppressing exaggerated cell growth might, therefore, be an effective management approach for CRSwNP in the future, given that the rapid recurrence of CRSwNP mimics benign tumour growth.
We further explored the gene sets enriched among the the various nasal tissues. Gene sets associated with epithelial differentiation and IFN responses were implicated in CRSwNP-NP versus CS-IT, which partially echoed our findings by IPA. Importantly, comparing CRSwNP-NP and CRSwNP-IT revealed the gene sets associated with inflammatory (particularly IL-1) responses, O-glycan biosynthesis and NP formation, which are likely the crucial risk factors that modulate early-stage CRSwNP. Significantly altered IFN-stimulated gene expression implicates the biological significance of abnormal antiviral responses in inducing NP formation.
Although eosinophilia reportedly predisposes NP formation, our analysis based on eosinophil enumeration surprisingly revealed finite DEGs (n=35 at false discovery rate <0.05), suggesting that eosinophilia was not a major confounder. Despite a previous study identifying the top DEGs such as CST1 and CCL18 in eosinophilic NP, no three-way pairwise comparisons such as ours have been made to control for false positives [37]. Similarly, coexisting allergic rhinitis did not yield strong DEGs.
Our study has important clinical implications. Our analytical pipeline has robustly identified the gene signatures and their associated signalling pathways and gene sets based on three pairwise comparisons with CS-IT as nondisease controls and CRSwNP-IT as disease controls. This approach minimised false-positive findings arising from the interactions between environmental and genetic factors. These DEGs represent a valuable database for future functional validation of gene expression patterns (i.e. cilia and mucin genes) in vivo and in vitro. The immune intolerance has shed light on a novel target for the management of CRSwNP, indicating the potential roles of immunomodulators. Furthermore, our findings may guide the selection of gene signatures to determine the efficacy of therapeutic interventions (including biologics such as mepolizumab, benralizumab and omalizumab) for future clinical trials.
Some limitations should be addressed. First, we sequenced heterogeneous tissues with multiple cell types and cell deconvolution was not conducted. Future studies that account for various cell types are needed. Second, we did not repeatedly sample the tissues, which might help address the differential gene expression before and after surgery or medication therapy (i.e. biologic exposure), nor did we identify the cardinal culprit DEGs responsible for recurrence of CRSwNP. Third, we only evaluated upper airway inflammation, but not systemic inflammation (i.e. eosinophilic versus noneosinophilic) via blood sampling. Fourth, the IT is not an ideal anatomical site to compare gene expression patterns with NP; however, the middle turbinate is frequently occupied by large-sized NP (NP size score >1), which renders surgical sampling challenging. Nonetheless, we found strong evidence of similar inflammatory responses in CRSwNP-NP and CRSwNP-IT, and the airway epithelium of CRSwNP-NP was similar overall to CRSwNP-IT. The prominent gene signatures associated with inflammatory responses in CRSwNP-IT indicate that the IT (control group) remains a valid source for sequencing investigations. Finally, we did not adjust RNA expression levels with age because it explained only a minority of the total variation.
In summary, we performed whole-transcriptomic sequencing with a novel three-way pairwise comparison, and identified in the CRSwNP-IT versus CS-IT comparison the important roles of gene sets orchestrating immune intolerance and a dysregulated immune response (i.e. defective Treg response), which may provide targets for therapeutic interventions. Our study provides valuable information for future investigations to functionally validate the gene expression patterns of defective host defences (including cilia dysfunction and immune dysregulation), inflammation and abnormal metabolism of the ECM, and to explore the biological responses to therapeutic interventions in CRSwNP.
Supplementary material
Supplementary Material
Please note: supplementary material is not edited by the Editorial Office, and is uploaded as it has been supplied by the author.
Supplementary material ERJ-00732-2019.SUPPLEMENT
Supplementary material figure E1 erj-00732-2019.supplement_figureE1
Supplementary material table E1 erj-00732-2019.supplement_tableE1
Supplementary material table E2 erj-00732-2019.supplement_tableE2
Supplementary material table E3 erj-00732-2019.supplement_tableE3
Supplementary material table E4 erj-00732-2019.supplement_tableE4
Supplementary material table E5 erj-00732-2019.supplement_tableE5
Supplementary material table E6 erj-00732-2019.supplement_tableE6
Supplementary material table E7 erj-00732-2019.supplement_tableE7
Supplementary material table E8 erj-00732-2019.supplement_tableE8
Supplementary material table E9 erj-00732-2019.supplement_tableE9
Supplementary material table E10 erj-00732-2019.supplement_tableE10
Supplementary material table E11 erj-00732-2019.supplement_tableE11
Supplementary material table E12 erj-00732-2019.supplement_tableE12
Supplementary material table E13 erj-00732-2019.supplement_tableE13
Shareable PDF
Supplementary Material
This one-page PDF can be shared freely online.
Shareable PDF ERJ-00732-2019.Shareable
Acknowledgements
The authors would like to thank the study participants who donated their samples for the advancement of scientific knowledge. The authors would also like to thank Jessica Tamanini from Insight Editing London (London, UK) for editing assistance for the manuscript.
Footnotes
This article has supplementary material available from erj.ersjournals.com
Author contributions: Conceived and designed the experiments: D.Y. Wang, A.K. Andiappan and W-J. Guan. Collection of samples: Y. Peng, X-X. Zi, T-F. Tian, L. Shi, J. Ye and Q-H. Qiu. Performed the experiments: Y. Peng, X-X. Zi, T-F. Tian and A.K. Andiappan. Data analysis: A.K. Andiappan, Y. Peng, X-X. Zi, T-F. Tian, W-J. Guan, B. Lee, J. Lum, S.A. Tang and K.S. Tan. Contributed reagents/materials/analysis tools: A.K. Andiappan and D.Y. Wang. Wrote the paper: W-J. Guan, A.K. Andiappan, K.S. Tan, Y. Peng, X-X. Zi, T-F. Tian. Critical review and approval: A.K. Andiappan, D.Y. Wang and W-J. Guan.
Conflict of interest: Y. Peng has nothing to disclose.
Conflict of interest: X-X. Zi has nothing to disclose.
Conflict of interest: T-F. Tian has nothing to disclose.
Conflict of interest: B. Lee has nothing to disclose.
Conflict of interest: J. Lum has nothing to disclose.
Conflict of interest: S.A. Tang has nothing to disclose.
Conflict of interest: K.S. Tan has nothing to disclose.
Conflict of interest: Q-H. Qiu has nothing to disclose.
Conflict of interest: J. Ye has nothing to disclose.
Conflict of interest: L. Shi reports grants from The Key Research Development Program of Shandong Province (2018CXGC1214), outside the submitted work.
Conflict of interest: W-J. Guan reports grants from Guangdong Province Universities and Colleges Pearl River Scholar Funded Scheme 2017, during the conduct of the study.
Conflict of interest: A.K. Andiappan reports grants from National Medical Research Council of Singapore (OFYIRG17nov065), during the conduct of the study.
Conflict of interest: D.Y. Wang reports grants from National Medical Research Council of Singapore (NMRC/CIRG/1458/2016), outside the submitted work.
Support statement: Supported by National Medical Research Council of Singapore (NMRC/CIRG/1458/2016; D.Y. Wang), National Medical Research Council of Singapore (OFYIRG17nov065; A.K. Andiappan), The Key Research Development Program of Shandong Province (2018CXGC1214; L. Shi), Guangdong Province Universities and Colleges Pearl River Scholar Funded Scheme 2017 (201710010097; W-J. Guan), and the National Nature Science Foundation of China (81873690; Q-H. Qiu). Funding information for this article has been deposited with the Crossref Funder Registry.
- Received April 11, 2019.
- Accepted August 8, 2019.
- Copyright ©ERS 2019