Identification of key genes and pathways in chronic rhinosinusitis with nasal polyps using bioinformatics analysis

Transcriptomic Differentiation of Phenotypes in Chronic Rhinosinusitis and Its Implications for Understanding the Underlying Mechanisms

Chronic rhinosinusitis (CRS) is a multifaceted disease with variable clinical courses and outcomes. We aimed to determine CRS-associated nasal-tissue transcriptome in clinically well-characterized and phenotyped individuals, to gain a novel insight into the biological pathways of the disease. RNA-sequencing of tissue samples of patients with CRS with polyps (CRSwNP), without polyps (CRSsNP), and controls were performed. Characterization of differently expressed genes (DEGs) and functional and pathway analysis was undertaken. We identified 782 common CRS-associated nasal-tissue DEGs, while 375 and 328 DEGs were CRSwNP- and CRSsNP-specific, respectively. Common key DEGs were found to be involved in dendritic cell maturation, the neuroinflammation pathway, and the inhibition of the matrix metalloproteinases. Distinct CRSwNP-specific DEGs were involved in NF-kβ canonical pathways, Toll-like receptor signaling, HIF1α regulation, and the Th2 pathway. CRSsNP involved the NFAT pathway and changes in the calcium pathway. Our findings offer new insights into the common and distinct molecular mechanisms underlying CRSwNP and CRSsNP, providing further understanding of the complex pathophysiology of the CRS, with future research directions for novel treatment strategies.

Identification of key genes and pathways in chronic rhinosinusitis with nasal polyps and asthma comorbidity using bioinformatics approaches

Patients with chronic rhinosinusitis with nasal polyps (CRSwNP) and asthma comorbidity (ACRSwNP) present severe symptoms and are more likely to relapse. However, the pathogenesis of ACRSwNP is not fully understood. The aim of this study was to explore the underlying pathogenesis of ACRSwNP using bioinformatics approaches. ACRSwNP-related differentially expressed genes (DEGs) were identified by the analysis of the GSE23552 dataset. The clusterProfiler R package was used to carry out functional and pathway enrichment analysis. A protein–protein interaction (PPI) network was built using the STRING database to explore key genes in the pathogenesis of ACRSwNP. The bioinformatics analysis results were verified through qRT-PCR. The Connectivity Map (CMap) database was used to predict potential drugs for the treatment of ACRSwNP. A total of 36 DEGs were identified, which were mainly enriched in terms of regulation of immune response and detection sensory perception of taste. Thirteen hub genes including AZGP1, AQP9, GAPT, PIP, and PRR4 were identified as potential hub genes in ACRSwNP from the PPI network. Analysis of the GSE41861 dataset showed that upregulation of CST1 in nasal mucosa was associated with asthma. qRT-PCR detection confirmed the bioinformatics analysis results. Tacrolimus and spaglumic acid were identified as potential drugs for the treatment of ACRSwNP from the CMap database. The findings of this study provide insights into the pathogenesis of ACRSwNP and may provide a basis for the discovery of effective therapeutic modalities for ACRSwNP.

Analysis of competing endogenous RNA (ceRNA) crosstalk in eosinophilic chronic rhinosinusitis with nasal polyps

BACKGROUND

Chronic rhinosinusitis with nasal polyps (CRSwNP) is one of the most common chronic inflammatory diseases, and has various phenotypes. Although its pathophysiology remains obscure, evidence has shown that dysregulation of noncoding RNAs (ncRNAs) is associated with CRSwNP. ncRNAs in the cytoplasm can act as competing endogenous RNAs (ceRNAs), which are involved in many inflammatory processes. However, the ceRNA crosstalk in CRSwNP is still unclear METHODS: We investigated expression profiles of messenger RNA (mRNA), microRNAs (miRNAs), and long noncoding RNAs (lncRNAs) in eosinophilic CRSwNP and constructed a global triple ceRNA network.

RESULTS

As a result, 964 differentially expressed mRNAs (DEmRs), 207 differentially expressed miRNAs (DEmiRs), and 15 differentially expressed lncRNAs (DElncRs) were identified, and a ceRNA network containing 598 miRNA-mRNA pairs and 70 lncRNA-miRNA pairs was finally constructed. Gene set enrichment analysis (GSEA) results indicated these DEmRs were mainly enriched in "cytokine-cytokine receptor interaction," "salivary secretion," "hematopoietic cell lineage," and "chemokine signaling pathway." Moreover, we also predicted the subcellular localization of the DElncRs identified in the network via bioinformatics approaches  CONCLUSION: In summary, the present study provided the first comprehensive assessment of the ceRNA crosstalk in eosinophilic CRSwNP. These findings will be of interest to the understanding of the potential pathophysiology of this disease.

Screening of Potential Biomarkers in the Peripheral Serum for Steroid-Induced Osteonecrosis of the Femoral Head Based on WGCNA and Machine Learning Algorithms

Background. Steroid-induced osteonecrosis of the femoral head (SONFH) has produced a substantial burden of medical and social experience. However, the current diagnosis is still limited. Thus, this study is aimed at identifying potential biomarkers in the peripheral serum of patients with SONFH. Methods. The expression profile data of SONFH (number: GSE123568) was acquired from the Gene Expression Omnibus (GEO) database. Differentially expressed genes (DEGs) in SONFH were identified and used for weighted gene coexpression network analysis (WGCNA). Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses were performed to investigate the biological functions. The protein-protein interaction (PPI) network and machine learning algorithms were employed to screen for potential biomarkers. Gene set enrichment analysis (GSEA), transcription factor (TF) enrichment analysis, and competing endogenous RNA (ceRNA) network were used to determine the biological functions and regulatory mechanisms of the potential biomarkers. Results. A total of 562 DEGs, including 318 upregulated and 244 downregulated genes, were identified between SONFH and control samples, and 94 target genes were screened based on WGCNA. Moreover, biological function analysis suggested that target genes were mainly involved in erythrocyte differentiation, homeostasis and development, and myeloid cell homeostasis and development. Furthermore, GYPA, TMCC2, and BPGM were identified as potential biomarkers in the peripheral serum of patients with SONFH based on machine learning algorithms and showed good diagnostic values. GSEA revealed that GYPA, TMCC2, and BPGM were mainly involved in immune-related biological processes (BPs) and signaling pathways. Finally, we found that GYPA might be regulated by hsa-miR-3137 and that BPGM might be regulated by hsa-miR-340-3p. Conclusion. GYPA, TMCC2, and BPGM are potential biomarkers in the peripheral serum of patients with SONFH and might affect SONFH by regulating erythrocytes and immunity.

Inflammatory molecular endotypes of nasal polyps derived from Caucasian and Japanese populations

BACKGROUND

Emerging evidence suggests that chronic rhinosinusitis with nasal polyps (CRSwNP) is a highly heterogeneous disease with disparate inflammatory characteristics between different racial groups and geographies. Little is known currently about possible distinguishing factors underlying these inflammatory differences.

OBJECTIVE

To interrogate for differences between Caucasian and Japanese CRSwNP disease using whole transcriptome and single-cell RNA gene expression profiling of nasal polyps (NPs).

METHODS

We performed whole transcriptome RNA sequencing (RNA-seq) with endotype stratification of NPs from 8 Caucasian (residing in USA) and 9 Japanese (residing in Japan) patients. Reproducibility was confirmed by qPCR in an independent validation set of 46 Caucasian and 31 Japanese patients. Single-cell RNA-seq stratified key cell types for contributory transcriptional signatures.

RESULTS

Unsupervised clustering analysis identified two major endotypes present within both NP cohorts, which have previously been reported at the cytokine level: 1) type 2 endotype and 2) non-type 2 endotype. Importantly, there was a statistically significant difference in the proportion of these endotypes between these geographically distinct NP subgroups (p = 0.03). Droplet-based single-cell RNA sequencing further identified prominent type 2 inflammatory transcript expression: C-C motif chemokine ligand 13 (CCL13) and CCL18 in M2 macrophages, as well as cystatin SN (CST1) and CCL26 in basal, suprabasal, and secretory epithelial cells.

CONCLUSION

NPs from both racial groups harbor the same two major endotypes, which we determine are present in differing ratios between each cohort with CRSwNP disease. Distinct inflammatory and epithelial cells contribute to the type 2 inflammatory profiles observed.

Comprehensive analysis of pivotal biomarkers, immune cell infiltration and therapeutic drugs for steroid-induced osteonecrosis of the femoral head

Steroid-induced osteonecrosis of the femoral head (SONFH) is a progressive disease that leads to an increased disability rate. This study aimed to ascertain biomarkers, infiltrating immune cells, and therapeutic drugs for SONFH. The gene expression profile of the GSE123568 dataset was downloaded from the Gene Expression Omnibus (GEO) database. The differentially expressed genes (DEGs) were identified using the NetworkAnalyst platform. Functional enrichment, protein-protein interaction network (PPI), and module analyses were performed using Metascape tools. An immune cell abundance identifier was used to explore immune cell infiltration. Furthermore, hub genes were identified based on maximal clique centrality (MCC) evaluation using cytoHubba application and confirmed by qRT-PCR using clinical samples. Finally, the L1000 platform was used to determine potential drugs for SONFH treatment. The SONFH mouse model was used to determine the therapeutic effects of aspirin. In total, 429 DEGs were identified in SONFH samples. Functional enrichment analysis showed that these DEGs were enriched in myeloid leukocyte activation and osteoclast differentiation processes. A set of nine immune cell types was confirmed to be markedly different between the SONFH and control samples. All 10 hub genes were significantly highly expressed in the serum of SONFH patients, as shown by qRT-PCR. Finally, the therapeutic effect of aspirin on SONFH was examined in animal experiments. Taken together, our data revealed the hub genes and infiltrating immune cells in SONFH, and we also screened potential drugs for use in SONFH treatment.

Gene Expression Profiles of Circular RNAs and MicroRNAs in Chronic Rhinosinusitis With Nasal Polyps

Introduction: Chronic rhinosinusitis (CRS) is often classified primarily on the basis of the absence or presence of nasal polyps (NPs), that is, as CRS with nasal polyps (CRSwNP) or CRS without nasal polyps (CRSsNP). Additionally, according to the percentage of eosinophils, CRSwNP can be further divided into eosinophilic CRSwNP (ECRSwNP) and non-ECRSwNP. CRSwNP is a significant public health problem with a considerable socioeconomic burden. Previous research reported that the pathophysiology of CRSwNP is a complex, multifactorial disease. There have been many studies on its etiology, but its pathogenesis remains unclear. Dysregulated expression of microRNAs (miRNAs) has been shown in psoriasis, rheumatoid arthritis, pulmonary fibrosis, and allergic asthma. Circular RNAs (circRNAs) are also involved in inflammatory diseases such as rheumatoid arthritis, septic acute kidney injury, myocardial ischemia/reperfusion injury, and sepsis-induced liver damage. The function of miRNAs in various diseases, including CRSwNP, is a research hotspot. In contrast, there have been no studies on circRNAs in CRSwNP. Overall, little is known about the functions of circRNAs and miRNAs in CRSwNP. This study aimed to investigate the expression of circRNAs and miRNAs in a CRSwNP group and a control group to determine whether these molecules are related to the occurrence and development of CRSwNP.

Methods: Nine nasal mucosa samples were collected, namely, three ECRSwNP samples, three non-ECRSwNP samples, and three control samples, for genomic microarray analysis of circRNA and microRNA expression. All of the tissue samples were from patients who were undergoing functional endoscopic sinus surgery in our department. Then we selected some differentially expressed miRNAs and circRNAs for qPCR verification. Meanwhile, GO enrichment analysis and KEGG pathway analysis were applied to predict the biological functions of aberrantly expressed circRNAs and miRNAs based on the GO and KEGG databases. Receiver operating characteristic (ROC) curve analysis and principal component analysis (PCA) were performed to confirm these molecules are involved in the occurrence and development of CRSwNP.

Results: In total, 2,875 circRNAs showed significant differential expression in the CRSwNP group. Specifically, 1794 circRNAs were downregulated and 1,081 circRNAs were upregulated. In the CRSwNP group, the expression of 192 miRNAs was significantly downregulated, and none of the miRNAs were significantly upregulated. GO and KEGG analysis showed differential circRNAs and miRNAs were enriched in “amoebiasis,” “salivary secretion,” “pathways in cancer,” and “endocytosis.” Through qRT-PCR verification, the expression profiles of hsa-circ-0031593, hsa-circ-0031594, hsa-miR-132-3p, hsa-miR-145-5p, hsa-miR-146a-5p, and hsa-miR-27b-3p were shown to have statistical differences. In addition, ROC curve analysis showed that the molecules with the two highest AUCs were hsa-circ-0031593 with AUC 0.8353 and hsa-miR-145-5p with AUC 0.8690. Through PCA with the six ncRNAs, the first principal component explained variance ratio was 98.87%. The AUC of the six ncRNAs was 0.8657.

Conclusion: In our study, the expression profiles of ECRSwNP and non-ECRSwNP had no statistical differences. The differentially expressed circRNAs and miRNAs between CRSwNP and control may play important roles in the pathogenesis of CRSwNP. Altered expression of hsa-circ-0031593 and hsa-miR-145-5p have the strongest evidence for involvement in the occurrence and development of CRSwNP because their AUCs are higher than the other molecules tested in this study.

The Role of Serum Metabolomics in Distinguishing Chronic Rhinosinusitis With Nasal Polyp Phenotypes

Background: Chronic rhinosinusitis with nasal polyps (CRSwNP) is a heterogeneous disease characterized by different clinical features and treatment responsiveness. This study aimed to compare the serum metabolomics profiles between eosinophilic CRSwNP (eCRSwNP) and non-eosinophilic CRSwNP (neCRSwNP) and healthy controls (HC) and explore objective biomarkers for distinguishing eCRSwNP before surgery. Methods: Serum samples were collected from 33 neCRSwNP patients, 37 eCRSwNP patients, and 29 HC. Serum metabolomics profiles were investigated by ultra-high-performance liquid chromatography-mass spectrometry. Results: The analysis results revealed that neCRSwNP, eCRSwNP, and HC exhibited distinctive metabolite signatures. In addition, eCRSwNP could be distinguished from neCRSwNP referring to their serum metabolic profiles, and the top ten different metabolites were citrulline, choline, linoleic acid, adenosine, glycocholic acid, L-serine, triethanolamine, 4-guanidinobutyric acid, methylmalonic acid, and L-methionine, which were related to several most important pathways including arginine and proline metabolism; glycine, serine, and threonine metabolism; linoleic acid metabolism; and purine metabolism. Among these distinctive metabolites, citrulline, linoleic acid, adenosine, and 4-guanidinobutyric acid showed good predictabilities, and the serum levels of citrulline, linoleic acid, and adenosine were significantly correlated with tissue eosinophil (T-EOS) percentage and T-EOS count. Conclusion: eCRSwNP patients exhibited discriminative serum metabolic signatures in comparison with neCRSwNP patients and HC. These results suggested that metabolomics profiles contributed to understanding the pathophysiological mechanisms of CRSwNP and distinguishing its phenotypes.

Nasal Polyposis: Insights in Epithelial-Mesenchymal Transition and Differentiation of Polyp Mesenchymal Stem Cells

Chronic rhinosinusitis is a common inflammatory disease of paranasal sinuses, which causes rhinorrhea, nasal congestion, and hyposmia. The genetic predisposition or the exposure to irritants can sustain the inflammatory response and the development of nasal polyposis. Nasal polyps are benign and teardrop-shaped growths that project in the nasal cavities, and originate from the ethmoid sinuses. This inflammatory process is associated with high expression of IL-4, IL-5 and IL-13 and IgE. Antibodies targeting these cytokines or receptors represent a therapeutic strategy in the treatment of nasal polyposis in combination with corticosteroids. The molecular pathogenesis of nasal polyps in chronic rhinosinusitis (CRS) patients is associated with remodeling transition, a process in which epithelial cells lose their typical phenotype, acquiring a mesenchymal-like aspect. TGFβ/SMAD, ERK, and Wnt/β-catenin pathways are altered during the nasal tissue remodeling. miRNA and inhibitor molecules targeting these signaling pathways are able to interfere with the process; which could lead to alternative therapies. Nasal polyps are an alternative source of mesenchymal stem cells, which can be isolated from surgical biopsies. A molecular understanding of the biology of PO-MSCs will contribute to the delineating inflammatory process underlying the development of nasal polyps.

Identification of Novel Genes and Biological Pathways That Overlap in Infectious and Nonallergic Diseases of the Upper and Lower Airways Using Network Analyses

Previous genetic studies on susceptibility to otitis media and airway infections have focused on immune pathways acting within the local mucosal epithelium, and outside of allergic rhinitis and asthma, limited studies exist on the overlaps at the gene, pathway or network level between the upper and lower airways. In this report, we compared [1] pathways identified from network analysis using genes derived from published genome-wide family-based and association studies for otitis media, sinusitis, and lung phenotypes, to [2] pathways identified using differentially expressed genes from RNA-sequence data from lower airway, sinus, and middle ear tissues, in particular cholesteatoma tissue compared to middle ear mucosa. For otitis media, a large number of genes (n = 1,806) were identified as differentially expressed between cholesteatoma and middle ear mucosa, which in turn led to the identification of 68 pathways that are enriched in cholesteatoma. Two differentially expressed genes CR1 and SAA1 overlap in middle ear, sinus, and lower airway samples and are potentially novel genes for otitis media susceptibility. In addition, 56 genes were differentially expressed in both tissues from the middle ear and either sinus or lower airways. Pathways that are common in upper and lower airway diseases, whether from published DNA studies or from our RNA-sequencing analyses, include chromatin organization/remodeling, endocytosis, immune system process, protein folding, and viral process. Taken together, our findings from genetic susceptibility and differential tissue expression studies support the hypothesis that the unified airway theory wherein the upper and lower respiratory tracts act as an integrated unit also applies to infectious and nonallergic airway epithelial disease. Our results may be used as reference for identification of genes or pathways that are relevant to upper and lower airways, whether common across sites, or unique to each disease.