- 1Department of Dermatology, The First Affiliated Hospital of Hunan Normal University/Hunan Provincial People’s Hospital, Changsha, China
- 2Department of Dermatology, Chinese Traditional Hospital of Changsha, Changsha, China
- 3Nursing Department, Hunan Provincial People’s Hospital/The First Affiliated Hospital of Hunan Normal University, Changsha, China
Background: The etiopathogenesis of chronic spontaneous urticaria (CSU) has not been fully understood, and there has been extensive interest in the interaction between inflammatory dermatosis and pyroptosis. This study intends to investigate the molecular mechanism of pyroptosis-related genes in CSU via bioinformatic ways, aiming at identifying the potential key biomarker.
Methods: GSE72540, the RNA expression profile dataset of CSU, was utilized as the training set, and GSE57178 as the validation set. Differently expressed pyroptosis-related genes (DEPRGs), GO, KEGG, and DO analyses were performed. The hub genes were explored by the protein–protein interaction analysis. Moreover, CIBERSORT was employed for estimating immune cell types and proportions. Then, we constructed a DEmRNA–miRNA–DElncRNA ceRNA network and a drug–gene interaction network. Finally, ELISA was used for gene expression analysis.
Results: We recognized 17 DEPRGs, whose enrichment analyses showed that they were mostly enriched in inflammatory response and immunomodulation. Moreover, 5 hub genes (IL1B, TNF, and IRF1 are upregulated, HMGB1 and P2RX7 are downregulated) were identified via the PPI network and verified by a validation set. Then immune infiltration analysis displayed that compared with normal tissue, CSU owned a significantly higher proportion of mast cells activated, but a lower proportion of T cells CD4 naive and so on. Furthermore, IL1B was statistically and positively associated with mast cells activated in CSU, and SNHG3, the upstream factor of IL1B in the ceRNA we constructed, also related with mast cells in CSU. Further analysis exhibited that the protein subcellular localization of IL1B was extracellular, according with its intercellular regulation role; IL1B was significantly correlated with key immune checkpoints; and the NOD-like receptor signaling pathway was the mainly involved pathway of IL1B based on the couple databases. What is more, the result of ELISA of CSU patients was the same as the above analyses about IL1B. In addition, the drug–gene interaction network contained 15 potential therapeutic drugs targeting IL1B, and molecular docking might make this relationship viable.
Conclusion: IL1B and its related molecules might play a key role in the development of CSU and could be potential biomarkers in CSU.
Introduction
Chronic spontaneous urticaria (CSU) is one of chronic inflammatory dermatosis, which is delineated as, for recognized or unrecognized causes, angioedema, wheal, or both occurring spontaneously for more than 6 weeks (1). The CSU’s prevalence is approximately 1% of the population (lifetime prevalence = 1.4%; point prevalence = 0.7%) (2). Moreover, CSU will get the increase of risk for comorbid autoimmune diseases like autoimmune thyroid disease (3). The frequently recurring symptoms, pruritus, urticaria, and angioedema, severely affect patients’ performance at school and work and impair their quality of life, which brings much encumbrance to both their households and society (4). Unlike acute urticaria, which is usually caused by an identifiable agent like an allergic reaction to a drug or other, the cause and pathogenesis of CSU are complex and remain largely unclear (5). Consequently, it is of great significance for individualized and effective treatment to reveal the pathogenesis and recognize key biomarkers of CSU.
The CSU’s etiopathogenesis has not been totally uncovered, but the existing studies suggest that maladjustment of inflammatory cells (such as mast cells and basophils) is the potentially core contributor (6). It is well known that a series of intracellular signaling cascades result in mast cell activation, after IgE bonds to the high-affinity IgE receptor. The activated mast cell releases proteases, histamines, and cytokines with the generation of platelet-activating mediators and other arachidonic acid metabolites (leukotrienes C4, D4, and E4 and prostaglandin D2). These cytokines give rise to vascular permeability and increased vasodilation, ensuing interstitial edema and sensory nerve stimulation and causing the obvious itchiness, redness, and swelling (7, 8). Besides, some CSU patients could show signs of activation of the coagulation/fibrinolytic system, such as significant elevation of serum factors like D-dimer, sICAM-1, and sVCAM-1 (9, 10). In this setting, some treatment strategies were developed, such as antihistamine, biological agent, and immunosuppressant. The first-line symptomatic treatment for CSU is largely depending on modern 2nd-generation H1 antihistamines, but standard-dosed antihistamines are ineffective in about 40% of the patients (11). Omalizumab, a humanized anti-IgE antibody, is the first licensed biological treatment by the Food and Drug Administration (FDA) for patients with CSU refractory to H1 antihistamines. However, relapse rates following the withdrawal of omalizumab are high (12). Ciclosporin is recommended for combination treatment in patients with severe disease refractory. Nevertheless, it is a problem which cannot be ignored that long-term use of ciclosporin leads to serious side effects (13). It can thus be seen that complete control of symptoms in the majority of patients remains a worldwide challenge. Consequently, further and fuller exploring the inflammatory reaction pathogenesis of CSU is scientifically significant to the clinical therapy.
There has been extensive interest in the interaction between inflammatory reaction and pyroptosis. The pyroptosis is denoted as inflammasome-dependent cell death (14). It was found that it makes a critical difference in the development of numerous inflammatory skin diseases. A research reported that Mdivi-1 significantly suppressed the pyroptotic cell death of keratinocytes and inhibited NLRP3 inflammasome activation to play a protective role in atopic dermatitis (15). Deng et al. indicated that cycloastragenol could inhibit the liberation of inflammatory mediators and macrophage infiltration in psoriasis by inhibiting the pyroptosis which NLRP3 mediated (16). In addition, previous research has established that aberrant NLRP3 inflammasome activation in masts cell contributes to histamine-independent urticaria by production of IL-1β in cryopyrin-associated periodic syndromes (CAPSs) (17). However, our understanding of CSU with pyroptosis is still pretty limited.
Microarray technology has been widely utilized into biological studies in recent years, and the data generated by it like the mRNA dataset could be an advantageous instrument for discovering critical factors of etiopathogenesis of diseases, which offers valuable insulation and foundation for further novel studies (18, 19). In this study, via applying the bioinformatic method, we analyzed the data of the Gene Expression Omnibus (GEO) (20), which has its origin in microarray technology, to explore immune cell infiltration and ceRNA network, and further reveal the molecular mechanism of pyroptosis-related genes in CSU and identify key biomarkers.
Material and Method
Microarray Data Source
The analysis process of this study is shown in Figure 1. We downloaded the datasets (GSE72540, GSE57178) from the GEO database (Table 1). GSE72540 contained 31 samples’ RNA expression profiling, selecting 10 CSU samples and 8 control samples, and GSE57178 contained 18 samples, selecting 6 CSU samples and 7 control samples.
Identifying Differently Expressed Pyroptosis-Related Genes
We normalized and preprocessed data and identified the different expression genes (DEGs) among the CSU sample and control via the GEO2R tool (21). |log2 FC| >1 and p < 0.05 as the cutoff. The 161 pyroptosis-related genes (PRGs) were downloaded from the GeneCards database (Supplementary Table S1) (22). Altogether consistent genes between DEGs and PRGs were identified as differently expressed pyroptosis-related genes (DEPRGs).
GO, KEGG, and DO Enrichment Analyses of DEPRGs
GO enrichment analysis [included MF (molecular function), BP (biological process), and CC (cellular component)] and KEGG pathway analysis were executed via the Metascape database (23). Min Enrichment ≥1.5, Min Overlap ≥3, and p < 0.01 were considered as the threshold. The WebGestalt tool (24) was used for DO enrichment analysis, and FDR ≤ 0.05 as the significance level.
Protein–Protein Interaction Network and Module Analyses
To investigate the protein–protein interaction (PPI) network, we used the STRING tool (25) and visualized it and analyzed the interactions of DEPRGs by the Cytoscape software (26). The Molecular Complex Detection (MCODE) plug-in was utilized for the module analysis of the PPI network. The cytoHubba tool was used for identifying the hub genes. The hub genes’ GO enrichment analysis was performed through the ClueGO plug-in.
Data Verification
The RNA expressed dataset GSE57178, containing 6 CSU lesion samples and 7 healthy control samples, was utilized as the validation set to verify the reliability of hub genes.
Immune Infiltration Analysis
The immune infiltration was calculated by the web tool CIBERSORT (27), which is a deconvolution algorithm that can evaluate the proportion of 22 infiltrating lymphocyte subsets in a large number of tissue samples. The GraphPad Prism 8.0.2 (San Diego, CA, USA) tool (28) was utilized for the correlation analysis between different immune cells, and between immune cells and hub genes, calculating the ratio of every kind of immune cell in CSU tissue and control.
Exploration ceRNA Network of the Hub Genes
To explore the miRNA–mRNA interaction of the ceRNA network, the potential miRNAs targeting the hub gene were identified via the TargetScan (29), miRNet (30), and DIANA TOOLS TarBase v.8 databases (31). If this was concurrently recognized in each database, the result was considered as true. Next, the possible lncRNAs targeting the miRNA were predicted through the miRNet database, which was cross-checked with the differently expressed lncRNA (DElncRNA) of CSU. LncRNA subcellular localization was predicted using lncLocator (32). The web-based tools, Wei Sheng Xin (http://www.bioinformatics.com.cn) and Draw Venn Diagram (http://bioinformatics.psb.ugent.be/webtools/Venn/), were used for data visualization.
Gene Set Enrichment Analysis
The Gene Set Enrichment Analysis (GSEA) tool (33) was used for exploring the molecular signaling pathway in which IL1B might be involved in CSU. The pathway enrichment analysis utilized the c2.cp.kegg.v7.3.symbols.gmt gene sets of the official website. False discovery rate q-value <0.01 was regarded as difference.
Analysis of Protein Subcellular Localization and Correlation With Immune Checkpoints
The protein subcellular localization of IL1B was predicted using the Cell-PLoc 2.0 tool (34), which is a package of web servers for predicting subcellular localization of proteins in different organisms. The correlation between IL1B and key immune checkpoints (35) such as HAVCR2(TIM3), LAG3, CTLA4, CD274(PDL1), PDCD1(PD1), and TIGIT were analyzed via Pearson’s correlation coefficient in GraphPad Prism 8.0.2.
Drug–Gene Interaction and Molecular Docking Analysis
To explore the drug–gene interaction, the DrugBank database (36) was utilized for identifying existing or/and potentially associated drug substances. Moreover, the Cytoscape software was utilized for data visualization. The molecular structure of the ligand and the target protein were obtained from the PubChem database (37) and PDB database (38). Docking simulations were conducted through AutoDock Vina (39) to generate the docking energy. The PyMOL software (40) was performed to visualize docked complexes.
Enzyme-Linked Immunosorbent Assay
To examine the protein levels of IL1B, NLRP3, and mast cell tryptase (MCT), serums from 10 CSU patients and 10 healthy controls were harvested for enzyme-linked immunosorbent assay (ELISA) (the CSU patients without medical treatment within 2 weeks and concomitant autoimmune diseases, whose detailed information is in Supplementary Table S2). Specific ELISA kits for IL1B (Neobioscience, Shenzhen, China), NLRP3 (uscnk, Wuhan, China), and MCT (Fufeng, Shanghai, China) were used according to the instructions of the manufacturer. Briefly, the standard samples, which were offered through the kit, of known concentration and the samples from two experimental groups were added to the kit plate and then incubated with the kit reagents (41). The OD450 values were detected using the microplate reader (Huisong, Shenzhen, China).
Statistical Analysis
The unpaired Student’s t-test was performed for data analysis of two groups. The potential correlation between the two variables was detected by Pearson’s correlation coefficient. p < 0.05 was considered as significance level. GraphPad Prism 8.0.2 was performed as the statistical software.
Results
Recognition of DEPRGs of CSU
The CSU RNA expression profile dataset (GSE72540) was normalized as shown in Figure 2A. 1,297 DEGs (containing 1,033 DEmRNAs and 173 DElncRNAs) were identified in the GSE72540 dataset (Supplementary Table S3), and their volcano plot is shown in Figure 2B. As shown in Figure 2C and Table 2, we identified 17 congruous DEPRGs via integrated bioinformatics analysis, including 13 congruously upregulated and 4 congruously downregulated. The heat map of DEPRGs is shown in Figure 2D.
Figure 2 Identification of DEPRGs of CSU. (A) Normalization of selected samples of GSE72540. (B) The differentially expressed genes of GSE72540. (C) The DEPRGs of CSU. (D) The heat map of the DEPRGs.
Function Enrichment Analyses of the DEPRGs
The GO analysis of DEPRGs was performed to reveal their biology functions. As shown, in the GO BP category, most of DEPRGs were mostly involved into regulation of cytokine production, interleukin-1 beta production, interleukin-1 production, etc. (Figure 3A). In the GO CC category, most of the DEPRGs were enriched into the membrane microdomain and membrane raft (Figure 3B). In the GO MF category, the main DEPRGs were enriched in cytokine activities, cytokine receptor binding, and signaling receptor activator activity, etc. (Figure 3C). The results of KEGG pathway enrichment exhibited that the mostly involved pathways were the NF-kappa B signaling pathway, TNF signaling pathway, and NOD-like receptor signaling pathway (Figure 3D). Utilizing the WebGestalt online database to further explore the function of DEPRGs, the result of DO enrichment showed that dermatomyositis, leishmaniasis, cutaneous, ulcerative colitis, etc., were the major diseases that DEPRGs participated in (Figure 3E). These suggested that inflammation and immune response were the major function of DEPRGs.
Figure 3 The top ten lists of function enrichment analyses of DEPRGs. (A) GO BP; (B) GO CC; (C) GO MF; (D) KEGG signaling pathway; (E) DO enrichment. **p < 0.01; ***p < 0.001.
PPI Network and Hub Gene Analyses
To reveal the interaction of each protein, the PPI network of the DEPRGs was built according to the STRING database, including 15 nodes and 34 edges. In the protein network graph, each node represented a protein, and the edge represented a connection between two proteins. Moreover, among the 15 nodes, 3 nodes were downregulated, and 12 were upregulated (Figure 4A). The targets were sorted by target connectivity from small to large in the PPI network, as shown in Figure 4B. The most important module was selected, including 10 edges and 8 nodes (Figure 4C). Hub genes were detected consistently via four algorithms (degree, MNC, stress, and MCC) of cytoHubba (Figure 4D). The top five gene scores were considered to be hub genes of CSU: IL1B, TNF, IRF1, HMGB1, and P2RX7 (Figure 4E and Table 3). Because the more closely knitted gene in the network is more fundamental to regulation, we further investigated the functions of the hub genes through the ClueGO plug-in. As Figure 4F shows, they were still primarily enriched in immunomodulation including regulation of adaptive immune response, lymphocyte proliferation, and regulation of phagocytosis. According to GSE57178, the mRNA expression of each hub gene manifested that, compared with the control, IL1B, IRF1, and TNF were significantly overexpressed while P2RX7 had a significantly lower expression in CSU, which was the same as the above results and indicated that IL1B, IRF1, P2RX7, and TNF were the key genes of CSU (Figure 4G).
Figure 4 The PPI network and hub gene analyses. (A) The PPI network of the DEPRGs, the bigger sizes of the edge and node mean the higher degree. The red means upregulated, and green means downregulated. (B) The connectivity rank of genes. (C) The first module of the PPI network. (D) Four algorithms were utilized to identified hub genes. (E) 5 hub genes of CSU. (F) The biological process of hub genes via the ClueGO tool. (G) Data validation of hub genes by GSE57178. *P < 0.05; **P < 0.01; ***P < 0.001.
Immune Infiltration Analysis
We investigated the difference among CSU tissues and control to explore the panorama of immune infiltration of CSU via the CIBERSORT algorithm. The ratio of 22 immune cells of samples is shown in Figure 5A. The correlation between each of immune cells is shown in Figure 5B, among which T cells CD4 memory activated were significantly correlated with macrophages M1 and macrophages M2, and eosinophils were statistically correlated with dendritic cells resting. At the side of control tissue, CSU owned a higher ratio of mast cells activated, and T cells CD4 naive, plasma cells, and B cells memory were significantly lower (Figure 5C). Next, we revealed the relation among the abundance of the immune cells and hub gene expression through the Pearson’s correlation coefficient. The results displayed that mast cells activated were statistically positively related to the levels of IL1B and TNF, but negatively to HMGB1’s; B cells memory and plasma cells were positively correlated with HMGB1 and P2RX7, but negatively with IL1B, IRF1, and TNF; T cells CD4 naive were positively correlated with HMGB1 and P2RX7, but negatively with IL1B (Figure 6).
Figure 5 Immune infiltration analysis of CSU. (A) The ratio of 22 immune cells of each sample of CSU. (B) The correlation between each of immune cells. (C) The proportion of immune cells in CSU and control.
Figure 6 The correlation between the hub gene and the immune cell. (A) HMGB1; (B) IL1B; (C) IRF1; (D) P2RX7; and (E) TNF.
The mRNA–miRNA–lncRNA ceRNA Network of CSU
The non-coding RNA (ncRNA) never participates in encoding proteins but was discovered to be involved in many biological functions, and perturbation of mRNA–miRNA–lncRNA ceRNA networks may affect diseases. The miRNA-targeting hub gene was concurrently recognized by all object databases as true, and their Venn diagrams are shown in Figure 7A. In addition, the possible lncRNA targeting the miRNA was predicted via miRNet online databases and got the intersection with 173 DElncRNAs. The lncRNA and mRNA of ceRNA must have a consistent expression trend according to the ceRNA mechanism. Then, we got 9 unique DElncRNAs based on the above. The lncRNAs, which compete with miRNAs by acting as ceRNAs to regulate the expression of mRNA targets, should be in the cytoplasm. As a result, only 4 DElncRNAs (HOTAIR, SCARNA9, SNHG3, and TUG1) were predicted in the cytoplasm by lncLocator (Figure 7B). Finally, an 18-axis ceRNA network (containing HMGB1/hsa-mir-17-5p/HOTAIR, IL1B/hsa-mir-21-5p/SNHG3, P2RX7/hsa-mir-20a-5p/HOTAIR and so on) was identified (Figure 7C). In addition, HOTAIR, SCARNA9, SNHG3, and TUG1 were significantly related to the major infiltration cell of CSU (Figure 8).
Figure 7 The construction of the lncRNA–miRNA–mRNA ceRNA network of CSU. (A) Venn diagram of miRNAs targeting each hub gene. (B) The subcellular localization of lncRNA of ceRNA. (C) The alluvial diagram of the ceRNA network.
Figure 8 The correlation between the lncRNA of ceRNA and immune cells. (A) HOTAIR; (B) TUG1; (C) SNG3; (D) SCARNA9.
GSEA of IL1B
Due to the fact that IL1B had been verified and that it played a role in immune infiltration and the ceRNA network of CSU, and simultaneously log2FC of IL1B was maximal in the hub genes, we chose IL1B for further analysis. The result of GSEA further verified the above results. As Figure 9 shows, besides ubiquitin-mediated proteolysis, arachidonic acid metabolism, cytosolic DNA sensing pathway, galactose metabolism, and apoptosis, IL1B was still mainly enriched in the NOD-like receptor signaling pathway.
Figure 9 The GSEA of IL1B. (A) NOD-like receptor signaling pathway. (B) Apoptosis. (C) Arachidonic acid metabolism. (D) Cytosolic DNA sensing pathway. (E) Galactose metabolism. (F) Ubiquitin-mediated proteolysis.
Protein Subcellular Localization and Correlation With Immune Checkpoint Analyses of IL1B
Different subcellular localizations of protein decide different biological functions. The protein subcellular localization of IL1B predicted by Cell-PLoc 2.0 was extracellular (Figure 10A). As displayed in Figure 10B, IL1B was significantly correlated with familiar immune checkpoints such as CD274(PDL1), CTLA4, HAVCR2(TIM3), and TIGIT, which indicated the important effect of IL1B in immune response further.
Figure 10 Integrated analyses of IL1B. (A) protein subcellular localization of IL1B. (B) The correlation between immune checkpoints and IL1B. (C) Drug–gene interaction network of IL1B. (D) Molecular docking between IL1B and minocycline.
Drug–Gene Interaction and Molecular Docking Analyses of IL1B
Developing potential therapeutic drugs for targeting IL1B provides a specific treatment strategy. The drug–gene interaction network of IL1B is exhibited in Figure 10C, in which there were 15 potential therapeutic drugs identified and 7 of them were approved (Table 4). Then, we worked out the molecular binding site of IL1B and minocycline, one of approved small-molecule drugs (Figure 10D).
IL1B Might Participate in Activation of Mast Cells via the NLRP3 in CSU
Compared with the healthy control, IL1B showed a significant overexpression in serum of CSU patients via ELISA, which is the same as our bioinformatic prediction (Figure 11A). The result of the ROC curve analysis showed that the area under the curve was 0.87 (p 0.01), which suggested the role of IL1B in diagnosis of CSU and further that it may be a potential biomarker in CSU (Figure 11B). Moreover, the result exhibited that MCT was overexpressed in CSU, and expression of MCT was statistically correlated with IL1B (Figures 11C, D). Due to the fact that MCT is the key marker of mast cell activation (42), it indicated that IL1B may participate in mast cell activation. NLRP3 is a subtype of NOD-like receptors and is famous for one of the key pyroptosis cytokines (43). The further ELISA result of CSU patient serum displayed that NLRP3 was significantly overexpressed and correlated with IL1B and MCT (Figures 11E–G). It testified the above bioinformatic prediction again which IL1B participated in, in the NOD-like receptor signaling pathway. Moreover, it also advised that it could be via the NLRP3 that IL1B participates in activation of mast cells.
Figure 11 IL1B might participate in activation of mast cells via the NLRP3 in CSU. (A) The expression of IL1B in CSU and control. (B) The ROC curve of IL1B. (C) The expression of MCT in CSU and control. (D) The correlation between IL1B and MCT in CSU. (E) The expression of NLRP3 in CSU and control. (F) The correlation between IL1B and NLRP3 in CSU. (G) The correlation between NLRP3 and MCT in CSU. **p < 0.01, ***p < 0.001, ****p < 0.0001.
Discussion
CSU is a common chronic inflammatory dermatosis, which has significantly negative impacts on the quality of people’s life owing to its repeated outbreaks and protracted course (44). Although current treatments of CSU get little effectiveness, how to more effectively mitigate and avert recurrence is still a global challenge as there are still many unknowns in its genesis. In addition, it has reached a consensus in the last guidelines that further research in some areas of CSU is needed, such as identification of mast cell/basophil-activating factors, identification of serum biomarkers of urticarial activity/mast cell activation, and identification of new histological markers (1). Remarkably, pyroptosis was one of deaths associated with cell membrane rupture. The increased number of cell membranes in mast cells might lead to the liberation of intracellular β-hexosaminidase and histamine (45). Since pyroptosis showed a great research prospect in inflammatory skin diseases, this work tries to identify and verify the potential key biomarkers of CSU from the standpoint of pyroptosis-related genes through bioinformatics ways, especially in inflammatory response, to provide a new perspective for the etiopathogenesis and therapeutic approaches of CSU.
In the present research, we recognized 1,297 DEGs from the CSU RNA expression profile. Then crossing the DEGS with pyroptosis-related genes, 17 DEPRGs (containing 4 downregulated genes and 13 upregulated genes) were recognized and then performed into gene function analysis. As shown, the DEPRGs were mostly involved in inflammatory response, as well as in pro-inflammatory effects (such as positive regulation of interleukin-8 production, positive regulation of interleukin-1 beta production, positive regulation of interleukin-1 production) and biological regulation (including signaling receptor activator activity, receptor ligand activity, cytokine activity), the majority of which are generally accredited to constituent parts of the development of CSU. A research containing 153 CSU patients suggested that the IL1 gene had a significant role in the susceptibility to CSU (46). Kasperska-Zajac et al. indicated that severity of systemic inflammation of CSU was related to elevated il-8 (47). The DEPRGs mainly participated in inflammatory pathways according to KEGG, likely the NF-kappa B signaling pathway, TNF signaling pathway, and NOD-like receptor signaling pathway. In addition, the result of DO further confirms the above. The DEPRGs were majorly enriched in inflammatory diseases like ulcerative colitis and fever. This advises that, to some extent, the DEPRGs could have a function to participate in the systemic inflammation of CSU.
Through the PPI network and module analyses, we identified five hub genes, namely, IL1B, TNF, IRF1 (all upregulated genes) and HMGB1, P2RX72 (both downregulated genes). They are the common inflammatory cytokines, but most of them have not been reported to be implicated in the development of CSU, so this would be a new finding. To fully explore the maladjustment of inflammatory cells of CSU, we executed immune infiltration analysis. It was found that CSU tissue owned a higher proportion of mast cells activated, but relatively lower ones of T cells CD4 naive, plasma cells, and B cells memory. Previous studies indicated that CSU was considered to be principally a mast cell-driven disease (48). Moreover, it has also been reported that the etiopathogenesis of CSU was closely associated with the dynamical unbalance of Th1/Th2 cells of CD4+T cells (49). However, since there have been few research reports, the relationship between CSU and plasma cells or/and B cells memory might be an interesting finding. Moreover, our research further showed that each of the hub genes (IL1B, TNF, IRF1, HMGB1, and P2RX7) was statistically related to major infiltration cells. Especially, IL1B and TNF were statistically and positively associated with mast cells activated, which suggests that they are related to maladjustment of inflammatory cells of CSU and might be its possible immunomodulation pivots. In addition, to reveal a systemically interactive modulation in CSU, we constructed a ceRNA network, in which there were 8 axes, such as IL1B/miR-21-5p/SNHG3, HMGB1/miR-34a-5p/TUG1, and P2RX7/miR-588/TUG1. It is worth noting that we also found that SNHG3 was significantly and negatively correlated with mast cells resting and plasma cells, and TUG1 was negatively related to mast cells activated, which further verified that the DEmRNA–miRNA–DElncRNA ceRNA network did have a critical role in maladjustment of inflammatory cells of CSU.
We chose IL1B to do further analysis for three reasons. First, it had the biggest fold change in hub genes. Second, it was verified by the validation set GSE57178. Third, both it and its upstream factor SNHG3 were related to the activation of mast cells. These indicated that IL1B could be in a more critical position in the development of CSU. IL1B (IL-1β) is a potent pro-inflammatory cytokine and plays a role in the innate and adaptive immunity of humans (50). Under the stimulation of immune response, inflammation, and infection, IL1B is released from monocytes, macrophages, and dendritic cells, affects local cells by paracrine, and targets distant cells via endocrine, ultimately leading to a series of inflammatory cascade responses like activation of immune cells and pyroptosis (51). Abnormal IL1B-related signaling pathways have been shown to be connected with some immune inflammatory diseases like SLE and UC (52, 53). These were in agreement with the result of our research. We predicted that the protein subcellular localization of IL1B was extracellular, which is in accord with its intercellular regulation role. Moreover, it was found that IL1B was significantly correlated with familiar immune checkpoints in CSU, such as PDL1, CTLA4, TIM3, and TIGIT, which someway showed its role of immune regulation in CSU. Furthermore, IL1B was an overexpression examined in clinical CSU patients by ELISA, and the ROC curve analysis confirmed the dependability of its diagnostic value; the point was that mast cells were significantly activated in CSU and IL1B did correlate with it, which verified our bioinformatic analyses and suggests that IL1B could be a potential prognostic and diagnostic biomarker in CSU.
The NOD-like receptor signaling pathway was the major involved pathway according to the enrichment analysis of IL1B in a couple of databases. Moreover, NLRP3 is a subtype of NOD-like receptors and is famous for one of key pyroptosis cytokines; moreover, it is a well-known activator of IL1B (54). The study of Guo et al. showed that the increased expression of NLRP3 in mast cells leads to the activation of caspase-1 and ultimately to production and secretion of IL-1β in endometriosis (55). These are further supported in our work. In clinical CSU patients, NLRP3 was statistically overexpressed and related to the activation of mast cells. NLRP3 was also significantly correlated with IL1B, which might advise that the pyroptosis-related signaling pathway is activated in CSU and it might be related to the activation of mast cells. What is more, it might be via the NOD-like receptor signaling pathway, NLRP3, that IL1B participates in activation of mast cells. Moreover, we further identified 15 potential therapeutic drugs targeting IL1B, which provides a possible therapeutic strategy for CSU. Molecular docking revealed that the exact molecular binding makes this relationship more reliable.
Our study also had some limitations. We measured gene expression levels using sera from clinical CSU patients rather than tissues, which is not good enough but still clinically representative. Besides, we will perform the experiments in vivo and in vitro to further confirm our results in the future.
In sum, we identified 5 hub genes, IL1B, TNF, IRF1, HMGB1, and P2RX7, from pyroptosis-related genes, which are mainly involved in the inflammatory response and maladjustment of inflammatory cells of CSU. Particularly, ILB and its ceRNA axis might play a role in the activation of mast cells of CUS, and this might be achieved via the NOD-like receptor signaling pathway (NLRP3). Therefore, ILB and its related molecules might be potential key biomarkers in the development of CSU, and our study would provide a new perspective for the etiopathogenesis and therapeutic programs of CSU.
Data Availability Statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/Supplementary Material.
Ethics Statement
Written informed consent was obtained from the individual(s)’ and minor(s)’ legal guardian/next of kin, for the publication of any potentially identifiable images or data included in this article.
Author Contributions
SP and HF conceptualized the study design. SP and TZ drafted the manuscript. HF revised the manuscript. SP, SZ, QT, and YY collected data and performed the analysis. All authors contributed to the article and approved the submitted version.
Funding
This work was financed by the Hunan Provincial Health Commission scientific research projects (B2016229; B2019055; 202114051780) and Hunan Provincial Innovation Foundation for Postgraduate (CX20200547).
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fimmu.2022.850993/full#supplementary-material
Supplementary Table S1 | 161 pyroptosis-related genes
Supplementary Table S2 | Population characteristics
Supplementary Table S3 | 1033 DEmRNAs and 173 DElncRNAs of CSU
Supplementary Table S4 | The DEGs of GSE72540
Supplementary Table S5 | The GO and KEGG enrichment of DEPRGs
Supplementary Table S6 | PPI topology table
Supplementary Table S7 | Hub genes-mirRNAs
Supplementary Table S8 | MiRNAs-lncRNAs
References
1. Zuberbier T, Abdul Latiff AH, Abuzakouk M, Aquilina S, Asero R, Baker D, et al. The International EAACI/GA²LEN/EuroGuiDerm/APAAACI Guideline for the Definition, Classification, Diagnosis, and Management of Urticaria. Allergy (2022) 77:734–66. doi: 10.1111/all.15090
2. Saini SS. Chronic Spontaneous Urticaria: Etiology and Pathogenesis. Immunol Allergy Clinics North America (2014) 34(1):33–52. doi: 10.1016/j.iac.2013.09.012
3. Kolkhir P, Altrichter S, Asero R, Daschner A, Ferrer M, Giménez-Arnau A, et al. Autoimmune Diseases Are Linked to Type IIb Autoimmune Chronic Spontaneous Urticaria. Allergy Asthma Immunol Res (2021) 13(4):545–59. doi: 10.4168/aair.2021.13.4.545
4. Gonçalo M, Gimenéz-Arnau A, Al-Ahmad M, Ben-Shoshan M, Bernstein JA, Ensina LF, et al. The Global Burden of Chronic Urticaria for the Patient and Society. Br J Dermatol (2021) 184(2):226–36. doi: 10.1111/bjd.19561
5. Zuberbier T, Aberer W, Asero R, Abdul Latiff AH, Baker D, Ballmer-Weber B, et al. The EAACI/GA²LEN/EDF/WAO Guideline for the Definition, Classification, Diagnosis and Management of Urticaria. Allergy (2018) 73(7):1393–414. doi: 10.1111/all.13397
6. Giménez-Arnau AM, DeMontojoye L, Asero R, Cugno M, Kulthanan K, Yanase Y, et al. The Pathogenesis of Chronic Spontaneous Urticaria: The Role of Infiltrating Cells. J Allergy Clin Immunol Pract (2021) 9(6):2195–208. doi: 10.1016/j.jaip.2021.03.033
7. Maurer M, Eyerich K, Eyerich S, Ferrer M, Gutermuth J, Hartmann K, et al. Urticaria: Collegium Internationale Allergologicum (CIA) Update 2020. Int Arch Allergy Immunol (2020) 181(5):321–33. doi: 10.1159/000507218
8. He L, Yi W, Huang X, Long H, Lu Q. Chronic Urticaria: Advances in Understanding of the Disease and Clinical Management. Clin Rev Allergy Immunol (2021) 61(3):424–48. doi: 10.1007/s12016-021-08886-x
9. Asero R. Severe CSU and Activation of the Coagulation/Fibrinolysis System: Clinical Aspects. Eur Ann Allergy Clin Immunol (2019) 52:15–7. doi: 10.23822/EurAnnACI.1764-1489.109
10. Puxeddu L, Panza F, Pratesi F, Bartaloni D, Casigliani R, Rocchi V, et al. CCL5/RANTES, sVCAM-1, and sICAM-1 in Chronic Spontaneous Urticaria. Int Arch Allergy Immunol (2013) 162(4):330–4. doi: 10.1159/000354922
11. Altrichter S, Staubach P, Pasha M, Singh B, Chang AT, Bernstein JA, et al. An Open-Label, Proof-Of-Concept Study of Lirentelimab for Antihistamine-Resistant Chronic Spontaneous and Inducible Urticaria. J Allergy Clin Immunol (2018) 98(7):641–7. doi: 10.2340/00015555-2941
12. Zhao Z, Cai T, Chen H, Chen L, Chen Y, Gao X, et al. Expert Consensus on the Use of Omalizumab in Chronic Urticaria in China. World Allergy Organ J (2021) 14(11):100610. doi: 10.1016/j.waojou.2021.100610
13. Hon KL, Leung AKC, Ng WGG, Loo SK. Chronic Urticaria: An Overview of Treatment and Recent Patents. Recent Patents Inflamm Allergy Drug Discovery (2019) 13(1):27–37. doi: 10.2174/1872213X13666190328164931
14. McKenzie BA, Dixit VM, Power C. Fiery Cell Death: Pyroptosis in the Central Nervous System. Trends Neurosci (2020) 43(1):55–73. doi: 10.1016/j.tins.2019.11.005
15. Li L, Mu Z, Liu P, Wang Y, Yang F, Han X. Mdivi-1 Alleviates Atopic Dermatitis Through the Inhibition of NLRP3 Inflammasome. Exp Dermatol (2021) 30(12):1734–44. doi: 10.1111/exd.14412
16. Deng G, Chen W, Wang P, Zhan T, Zheng W, Gu Z, et al. Inhibition of NLRP3 Inflammasome-Mediated Pyroptosis in Macrophage by Cycloastragenol Contributes to Amelioration of Imiquimod-Induced Psoriasis-Like Skin Inflammation in Mice. Int Immunopharmacol (2019) 74:105682. doi: 10.1016/j.intimp.2019.105682
17. Kambe N, Nakamura Y, Saito M, Nishikomori R. The Inflammasome, an Innate Immunity Guardian, Participates in Skin Urticarial Reactions and Contact Hypersensitivity. Allergol Int: Off J Jpn Soc Allergol (2010) 59(2):105–13. doi: 10.2332/allergolint.09-RAI-0160
18. Cortés-Ciriano I, Gulhan DC, Lee JJ, Melloni GEM, Park PJ. Computational Analysis of Cancer Genome Sequencing Data. Nat Rev Genet (2022) 6(2):e2100785. doi: 10.1002/smtd.202100785
19. Abbott TR, Dhamdhere G, Liu Y, Lin X, Goudy L, Zeng L, et al. Development of CRISPR as an Antiviral Strategy to Combat SARS-CoV-2 and Influenza. Cell (2020) 181(4):865–76.e12. doi: 10.1016/j.cell.2020.04.020
20. Clough E, Barrett T. The Gene Expression Omnibus Database. Methods Mol Biol (Clifton NJ) (2016) 1418:93–110. doi: 10.1007/978-1-4939-3578-9_5
21. Barrett T, Wilhite SE, Ledoux P, Evangelista C, Kim IF, Tomashevsky M, et al. NCBI GEO: Archive for Functional Genomics Data Sets–Update. Nucleic Acids Res (2013) 41(Database issue):D991–5. doi: 10.1093/nar/gks1193
22. Shao Y, Jia H, Huang L, Li S, Wang C, Aikemu B, et al. An Original Ferroptosis-Related Gene Signature Effectively Predicts the Prognosis and Clinical Status for Colorectal Cancer Patients. Front Oncol (2021) 11:711776. doi: 10.3389/fonc.2021.711776
23. Zhou Y, Zhou B, Pache L, Chang M, Khodabakhshi AH, Tanaseichuk O, et al. Metascape Provides a Biologist-Oriented Resource for the Analysis of Systems-Level Datasets. Nat Commun (2019) 10(1):1523. doi: 10.1038/s41467-019-09234-6
24. Liao Y, Wang J, Jaehnig EJ, Shi Z, Zhang B. WebGestalt 2019: Gene Set Analysis Toolkit With Revamped UIs and APIs. Nucleic Acids Res (2019) 47(W1):W199–w205. doi: 10.1093/nar/gkz401
25. Szklarczyk D, Gable AL, Nastou KC, Lyon D, Kirsch R, Pyysalo S, et al. The STRING Database in 2021: Customizable Protein-Protein Networks, and Functional Characterization of User-Uploaded Gene/Measurement Sets. Nucleic Acids Res (2021) 49(D1):D605–d12. doi: 10.1093/nar/gkaa1074
26. Doncheva NT, Morris JH, Gorodkin J, Jensen LJ. Cytoscape StringApp: Network Analysis and Visualization of Proteomics Data. J Proteome Res (2019) 18(2):623–32. doi: 10.1021/acs.jproteome.8b00702
27. Chen B, Khodadoust MS, Liu CL, Newman AM, Alizadeh AA. Profiling Tumor Infiltrating Immune Cells With CIBERSORT. Methods Mol Biol (Clifton NJ) (2018) 1711:243–59. doi: 10.1007/978-1-4939-7493-1_12
28. Mitteer DR, Greer BD, Randall KR, Briggs AM. Further Evaluation of Teaching Behavior Technicians to Input Data and Graph Using GraphPad Prism. Behav Anal (Washington DC) (2020) 20(2):81–93. doi: 10.1037/bar0000172
29. Shi Y, Yang F, Wei S, Xu G. Identification of Key Genes Affecting Results of Hyperthermia in Osteosarcoma Based on Integrative ChIP-Seq/TargetScan Analysis. Med Sci Monit: Int Med J Exp Clin Res (2017) 23:2042–8. doi: 10.12659/MSM.901191
30. Chang L, Zhou G, Soufan O, Xia J. Mirnet 2.0: Network-Based Visual Analytics for miRNA Functional Analysis and Systems Biology. Nucleic Acids Res (2020) 48(W1):W244–51. doi: 10.1093/nar/gkaa467
31. Bozgeyik I. Therapeutic Potential of miRNAs Targeting SARS-CoV-2 Host Cell Receptor ACE2. Meta Gene (2021) 27:100831. doi: 10.1016/j.mgene.2020.100831
32. Lin Y, Pan X, Shen HB. Lnclocator 2.0: A Cell-Line-Specific Subcellular Localization Predictor for Long non-Coding RNAs With Interpretable Deep Learning. Bioinf (Oxford England) (2021) 1–9. doi: 10.1093/bioinformatics/btab127
33. Subramanian A, Tamayo P, Mootha VK, Mukherjee S, Ebert BL, Gillette MA, et al. Gene Set Enrichment Analysis: A Knowledge-Based Approach for Interpreting Genome-Wide Expression Profiles. Proc Natl Acad Sci United States America (2005) 102(43):15545–50. doi: 10.1073/pnas.0506580102
34. Chou KC, Shen HB. Cell-PLoc: A Package of Web Servers for Predicting Subcellular Localization of Proteins in Various Organisms. Nat Protoc (2008) 3(2):153–62. doi: 10.1038/nprot.2007.494
35. Liu D. Cancer Biomarkers for Targeted Therapy. biomark Res (2019) 7:25. doi: 10.1186/s40364-019-0178-7
36. Wishart DS, Feunang YD, Guo AC, Lo EJ, Marcu A, Grant JR, et al. DrugBank 5.0: A Major Update to the DrugBank Database for 2018. Nucleic Acids Res (2018) 46(D1):D1074–d82. doi: 10.1093/nar/gkx1037
37. Kim S, Chen J, Cheng T, Gindulyte A, He J, He S, et al. PubChem in 2021: New Data Content and Improved Web Interfaces. Nucleic Acids Res (2021) 49(D1):D1388–d95. doi: 10.1093/nar/gkaa971
38. Karuppasamy MP, Venkateswaran S, Subbiah P. PDB-2-PBv3.0: An Updated Protein Block Database. J Bioinf Comput Biol (2020) 18(2):2050009. doi: 10.1142/S0219720020500092
39. Nguyen NT, Nguyen TH, Pham TNH, Huy NT, Bay MV, Pham MQ, et al. Autodock Vina Adopts More Accurate Binding Poses But Autodock4 Forms Better Binding Affinity. J Chem Inf Modeling (2020) 60(1):204–11. doi: 10.1021/acs.jcim.9b00778
40. Lam WW, Siu SW. PyMOL Mcontrol: Manipulating Molecular Visualization With Mobile Devices. Biochem Mol Biol Educ: Bimonthly Publ Int Union Biochem Mol Biol (2017) 45(1):76–83. doi: 10.1002/bmb.20987
41. Aydin S. A Short History, Principles, and Types of ELISA, and Our Laboratory Experience With Peptide/Protein Analyses Using ELISA. Peptides (2015) 72:4–15. doi: 10.1016/j.peptides.2015.04.012
42. Payne V, Kam PC. Mast Cell Tryptase: A Review of its Physiology and Clinical Significance. Anaesthesia (2004) 59(7):695–703. doi: 10.1111/j.1365-2044.2004.03757.x
43. Swanson KV, Deng M, Ting JP. The NLRP3 Inflammasome: Molecular Activation and Regulation to Therapeutics. Nat Rev Immunol (2019) 19(8):477–89. doi: 10.1038/s41577-019-0165-0
44. Hoskin B, Ortiz B, Paknis B, Kavati A. Humanistic Burden of Refractory and Nonrefractory Chronic Idiopathic Urticaria: A Real-world Studyin the United States. Clin Ther (2019) 41(2) 205–20. doi: 10.1016/j.clinthera.2018.12.004
45. Liu M, Lu J, Chen Y, Shi X, Li Y, Yang S, et al. Sodium Sulfite-Induced Mast Cell Pyroptosis and Degranulation. J Agric Food Chem (2021) 69(27):7755–64. doi: 10.1021/acs.jafc.1c02436
46. Brzoza Z, Rymarczyk B, Grzeszczak W, Trautsolt W, Pluta-Kubicz M, Moczulski D. Interleukin 1 Gene Polymorphisms Presumably Participate in the Pathogenesis of Chronic Spontaneous Autoreactive Urticaria. J Interferon Cytokine Res: Off J Int Soc Interferon Cytokine Res (2020) 40(10):497–500. doi: 10.1089/jir.2020.0103
47. Kasperska-Zając A, Damasiewicz-Bodzek A, Grzanka R, Skrzypulec-Frankel A, Bieniek K, Sikora-Żydek A, et al. Circulating Soluble LIGHT/TNFSF14 is Increased and Associated With IL-8 Concentration in Chronic Spontaneous Urticaria. Int J Immunopathol Pharmacol (2018) 32:2058738418784431. doi: 10.1177/2058738418784431
48. Ulambayar B, Chen YH, Ban GY, Lee JH, Jung CG, Yang EM, et al. Detection of Circulating IgG Autoantibody to FcεRIα in Sera from Chronic Spontaneous Urticaria Patients. J Microbiol Immunol Infect (2020) 53(1):141–7. doi: 10.1016/j.jmii.2017.10.003
49. Chen Q, Zhong H, Chen WC, Zhai Z, Zhou Z, Song Z, et al. Different Expression Patterns of Plasma Th1-, Th2-, Th17- and Th22-Related Cytokines Correlate With Serum Autoreactivity and Allergen Sensitivity in Chronic Spontaneous Urticaria. J Eur Acad Dermatol Venereol: JEADV (2018) 32(3):441–8. doi: 10.1111/jdv.14541
50. Dinarello CA. A Clinical Perspective of IL-1β as the Gatekeeper of Inflammation. Eur J Immunol (2011) 41(5):1203–17. doi: 10.1002/eji.201141550
51. Wang Y, Che M, Xin J, Zheng Z, Li J, Zhang S. The Role of IL-1β and TNF-α in Intervertebral Disc Degeneration. Biomed Pharmacother = Biomed Pharmacother (2020) 131:110660. doi: 10.1016/j.biopha.2020.110660
52. Mohammadoo-Khorasani M, Salimi S, Tabatabai E, Sandoughi M, Zakeri Z, Farajian-Mashhadi F. Interleukin-1β (IL-1β) & IL-4 Gene Polymorphisms in Patients With Systemic Lupus Erythematosus (SLE) & Their Association With Susceptibility to SLE. Indian J Med Res (2016) 143(5):591–6. doi: 10.4103/0971-5916.187107
53. Dinallo V, Marafini I, Di Fusco D, Laudisi F, Franzè E, Di Grazia A, et al. Neutrophil Extracellular Traps Sustain Inflammatory Signals in Ulcerative Colitis. J Crohn’s Colitis (2019) 13(6):772–84. doi: 10.1093/ecco-jcc/jjy215
54. Mao L, Kitani A, Hiejima E, Montgomery-Recht K, Zhou W, Fuss I, et al. Bruton Tyrosine Kinase Deficiency Augments NLRP3 Inflammasome Activation and Causes IL-1β-Mediated Colitis. J Clin Invest (2020) 130(4):1793–807. doi: 10.1172/JCI128322
Keywords: chronic spontaneous urticaria (CSU), pyroptosis-related genes, IL1B, bioinformatics, inflammation, immunology
Citation: Peng S, Zhang T, Zhang S, Tang Q, Yan Y and Feng H (2022) Integrated Bioinformatics and Validation Reveal IL1B and Its Related Molecules as Potential Biomarkers in Chronic Spontaneous Urticaria. Front. Immunol. 13:850993. doi: 10.3389/fimmu.2022.850993
Received: 08 January 2022; Accepted: 22 February 2022;
Published: 18 March 2022.
Edited by:
Richard Williams, University of Oxford, United KingdomReviewed by:
Basavaraj Mallikarjunayya Vastrad, KLE Society’s College of Pharmacy, IndiaIlaria Puxeddu, University of Pisa, Italy
Copyright © 2022 Peng, Zhang, Zhang, Tang, Yan and Feng. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Hao Feng, ZG9jdG9yZmVuZ2hhb0AxMjYuY29t