- 1School of Nursing and Rehabilitation, Cheeloo College of Medicine, Shandong University, Jinan, China
- 2Department of Pediatrics, The Second Hospital of Shandong University, Jinan, China
- 3Rehabilitation Center, Qilu Hospital, Cheelo College of Medicine, Shandong University, Jinan, China
Background: Currently, disease control in patients with severe eosinophilic asthma is not optimistic. Competing endogenous RNA (ceRNA) networks have been found to play a key role in asthma in recent years. However, it is unclear whether ceRNA networks play an important part in severe eosinophilic asthma.
Methods: Firstly, gene expression profiles related to severe eosinophilic asthma were downloaded from the Gene Expression Omnibus (GEO) database. Secondly, the key modules were identified by the weighted gene co-expression network analysis (WGCNA). Thirdly, genes in modules highly associated with severe eosinophilic asthma were selected for further construction of the ceRNA network. Fourthly, Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses were performed on hub genes. Finally, the results of this study were validated on the GSE143303, GSE137268, and GSE147878 datasets.
Results: 22 severe eosinophilic asthmatics and 13 healthy controls were extracted for WGCNA. We found that the genes in the black module (r = −0.75, p < 0.05) and yellow module (r = 0.65, p < 0.05) were highly associated with severe eosinophilic asthma. EP300 was discovered to serve the key connecting function in the ceRNA network. Surprisingly, lncRNAs seem to eliminate the role of EP300 in the black module and we discovered that CCT8 and miRNA-mRNA formed a circRNA-miRNA-mRNA network in the yellow module. We found that EP300 and FOXO3 in the black module were regulated by steroid hormones in the enrichment analysis, which were related to the medication used by the patient. Through validation of other datasets, we found that the hub genes in the yellow module were the key genes in the treatment of severe eosinophilic asthma. In particular, RPL17 and HNRNPK might specifically regulate severe eosinophilic asthma.
Conclusion: RPL17 and HNRNPK might particularly regulate severe eosinophilic asthma. Our results could be useful to provide potential immunotherapy targets and prognostic markers for severe eosinophilic asthma.
Introduction
Asthma is an intractable chronic inflammatory illness of the airway caused by complicated genetic and environmental factors, with a wide range of types and intensities of airway inflammation and remodeling (Papi et al., 2018). The majority of patients get standardized medication and care, and their symptoms are successfully managed; nevertheless, 5–10% of asthma patients require high-dose inhaled therapy, yet they have trouble controlling their condition, which is referred to as severe asthma (Lemière et al., 2006; Haselkorn et al., 2009). Severe asthma, which encompasses a wide range of symptoms, is becoming recognized as a highly heterogeneous illness with a wide range of molecular, biochemical, and cellular inflammatory characteristics. Eosinophilic asthma is the most prevalent subtype of severe asthma (Jatakanon et al., 2000; Lemière et al., 2006). The molecular processes underlying the incidence and progression of severe eosinophilic asthma remain unknown. As a result, we need to learn more about them to properly treat severe eosinophilic asthma.
The influence of competitive endogenous RNA (ceRNA) on asthma has recently attracted the interest of researchers as a potential new mechanism for improved asthma therapy. Salmena et al. introduced the ceRNA hypothesis as a unique regulatory mechanism between non-coding RNA (ncRNA) and coding messenger RNA (mRNA) (Thomson and Dinger, 2016). MicroRNA (miRNA)-response elements (MREs), which operate as ceRNAs and play a critical role in different clinical processes, are found in long non-coding RNA (lncRNAs), pseudogene transcripts, circular RNAs (circRNA), viral RNAs, and protein-coding transcripts (Thomson and Dinger, 2016; Fan et al., 2018). With the further development of molecular detection technology, more and more researchers have found that most lncRNAs and circRNAs are involved in the pathogenesis of asthma through the lncRNA-miRNA-mRNA axis and circRNA-miRNA-mRNA axis to regulate the Th1/Th2 balance, M2 macrophage activation, and cytokine (IL-6, IL-13, and IL-17) secretion, respectively (Lee et al., 2009; Han et al., 2019; Huang et al., 2019; Qiu et al., 2019; Shang et al., 2020). The weakness of the previous study is that it focuses mainly on the molecular mechanisms of ceRNA in asthma and not on whether ceRNA plays a role in severe eosinophilic asthma.
Weighted Gene Co-expression analysis (WGCNA) is a systems biology approach that aims to find clusters of genes that are highly correlated with external clinical features to identify potential biomarkers and provide molecular targets for the treatment of disease (Langfelder and Horvath, 2008). Previous researchers have also used the WGCNA approach to explore differential genes in asthma, but previous studies have focused on populations with mild, moderate, and severe asthma and have not explored differential genes between severe eosinophilic asthma and healthy populations (He et al., 2020; Liao et al., 2020). Second, previous studies relied just on WGCNA to identify important genes and did not go on to build ceRNA networks (Kelly et al., 2018; Zhang et al., 2021). Based on this, we used the WGCNA approach to further construct the ceRNA networks to reveal the hidden intrinsic molecular mechanism of severe eosinophilic asthma and provide evidence to support the biologically targeted therapy for severe eosinophilic asthma.
Materials and Methods
Data Acquisition and Processing
The GSE143303 (GEO Accession viewer, 2021) was a dataset linked to severe eosinophilic asthma that was retrieved from the Gene Expression Omnibus (GEO) Datasets (https://www.ncbi.nlm.nih.-gov/gds/). GPL10558 was the platform number. Endobronchial biopsies were used to compare 47 samples of distinct inflammatory phenotypes (neutrophilic, eosinophilic, and paucigranulocytic) of severe asthma to 13 healthy controls. For our study, we used endobronchial biopsies from 22 severe eosinophilic asthmatics and 13 healthy controls. 13 healthy controls had a predicted FEV1 of more than 80% and no underlying heart or lung illness. Current smokers were not allowed to participate. The global initiative for asthma (GINA) 2019 (Sánchez-Ovando et al., 2021) defined severe asthma. Patients with severe asthma had no recent history of a clinical chest or upper respiratory tract infection. They also had inhaled corticosteroid (ICS) (>500 μg fluticasone or equivalent per day) or oral corticosteroid (OCS) and were at GINA stages 4–5. The inflammatory phenotypes were defined using the 95th percentile of the differential count and total cell count, and we called it eosinophilic asthma when there were ≥3.50% eosinophils and<71.75% neutrophils (Sánchez-Ovando et al., 2020).
The R Bioconductor package affy was used to normalize raw microarray gene expression data, which was then submitted to multiple quality control techniques. Using annotation information, gene IDs were then mapped to microarray probes. The mean expression value of genes assessed by several probes was determined after probes matching more than one gene were removed from the dataset. Subsequently, the top 5000 genes from the GSE143303 dataset were then screened using the median absolute deviation (MAD) value (Zhao et al., 2021). Following that, we load clinical characteristics data such as inflammatory phenotypes, age, gender, and smoking history (Miller et al., 2010). The p value was adjusted by the false discovery rate (FDR) approach. The characteristics of the participants were listed in Supplementary Table S1.
Network Construction and Consensus Module Detection
WGCNA was able to distinguish genes into multiple clusters, and further investigate the relationship between co-expression modules and clinical phenotypes. The WGCNA package from Bioconductor was used to build co-expression networks for all genes (Langfelder and Horvath, 2008). 1) The hclust function was used to cluster samples and check for outliers; 2) The soft-thresholding power was calculated in the construction of each module using the pickSoftThreshold function of WGCNA, which calculates the scale-free topology fit index for a set of candidate powers ranging from 1 to 20 and provides a suitable power value for network construction. The suitable power was determined if the index value for the reference dataset exceeded 0.85; 3) one-step network building was performed to find co-expression modules, and the limited minimum gene number was set at 50. We assessed the relationship between modules and clinical features. Sample-specific characteristics (e.g., age and smoking) might potentially influence the relationship between gene expression and severe eosinophilic asthma, so logistic regression was used to control for confounding variables in the IBM SPSS Statistics 24.0 software.
Relating Modules to External Clinical Traits
We were able to find genes with high group significance as well as high module membership in intriguing modules using the gene significance (GS) and module membership (MM) measures. The clinically significant module for severe eosinophilic asthmatics was identified if: |GS| ≥ 0.5 and |MM| ≥ 0.6, the correlation between MM and GS in the module was statistically significant (p < 0.05).
The search tool for retrieval of interacting genes (STRING) online database (https://string-db.org/) was an online website that could build a protein-protein interaction network (PPI) based on bioinformatics predictions or biochemical experimental results (STRING, 2021). In this study, the key modules were displayed using the STRING website to create a PPI network detecting gene connections with a threshold of interaction score >0.4. Genes with |GS| > 0.5 and |MM| > 0.6 in the key module were imported to Cytoscape (version 3.8.2). The hub genes were chosen from the top 12 genes.
Enrichment Analysis
To further visualize the activities of genes in the key module, Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses were performed on hub genes in the Cytoscape plug-in ClueGo (Bindea et al., 2009). The cutoff threshold was set at a P value of less than 0.05. ClueGo classified the signal pathways discovered by enrichment analysis into groups based on functional connection; the same group was colored the same color, and the labels of each group of the most essential terms were color-coded.
Construction of Competing Endogenous RNA
The RNA Interactome Database (RNAInter) brought together experimentally validated and computationally predicted RNA interactome data from miRTarBase and starBase, including RNA–RNA, and RNA–protein interactions (Kang et al., 2021). The stronger the evidence for a link between the two genes, the closer the confidence score was to 1. The top 12 genes in the black, red, and yellow modules were entered into the RNAInter website respectively to get ceRNA network relationships, and visualization of interacting genes with a confidence score >0.55 in Cytoscape.
Validation
We validated key genes with three datasets from the GEO database. We retrieved the GSE143303 dataset from patients with severe non-eosinophilic asthma (severe neutrophilic asthma and severe paucigranulocytic asthma) (n = 25) and healthy controls (n = 13). The aim was to see if the findings were particular to severe asthma or severe eosinophilic asthma. The characteristics of the participants were listed in Supplementary Table S2. We used the GSE147878 dataset to see if the results of this study could be replicated in patients with severe asthma. The GSE147878 dataset was a cross-sectional study from endobronchial biopsies (n = 73) and induced sputum (n = 44). We extracted the transcriptomic data of bronchoscopic biopsy tissue from severe asthma (n = 42) and healthy controls (n = 13). The characteristics of the participants were listed in Supplementary Table S3. Finally, the GSE137268 dataset was induced sputum samples from asthmatics and healthy controls. To see if the findings of this study were specific to eosinophilic asthma or severe eosinophilic asthma, we took non-severe eosinophilic asthma patients (controlled and uncontrolled eosinophilic asthma) (n = 13) and healthy controls (n = 15) in GSE137268. The characteristics of the participants were listed in Supplementary Table S4. The normalized expression values of hub genes were imported into IBM SPSS Statistics 24.0 software, and the differences between the two groups were analyzed using an independent samples t-test.
Result
Co-Expression Network Construction
After clustering all samples, we discovered an outlier (sample GSE 4256708) in Supplementary Figure S1. The scale-free topology index exceeded 0.9 when the soft-power β was set to 8 (Supplementary Figure S2). There were 11 modules with sizes ranging from 55 to 1304 genes, labeled 0 through 10 in the order of decreasing size. The number 0 was set aside for genes that did not belong to any of the modules (Supplementary Figure S3).
Identification of the Clinically Significant Module and Hub Genes
The correlation between module eigengene and clinical features was used to identify module-trait associations in Figure 1. The black and red modules were shown to be adversely associated with severe eosinophilic asthma, with correlations of 0.75 and 0.52, respectively (p < 0.05). This indicated that genes in the black and red modules were predominantly downregulated in severe eosinophilic asthmatics. The yellow module was recognized as the positive module with a correlation of 0.65 (p < 0.05).
FIGURE 1. Module-trait associations. Each row corresponds to a module eigengene, column to a clinical trait. The corresponding correlations and P-values were presented.
Figures 2A,B, and Supplementary Figure S4 showed that the black, yellow, and red modules had a strong GS-MM correlation (p < 0.05), which were identified as the clinically significant module and visualized in STRING (Supplementary Figures S5A–C). In total, 94, 39, and 81 genes with |GS| > 0.5 and |MM| > 0.6 in the black, yellow, and red modules were imported into Cytoscape, respectively, the top 12 genes were filtered as the hub genes in the black, yellow, and red module respectively (Figures 3A,B, and Supplementary Figure S6). Finally, we found that the most connected hub gene was E1A-associated 300-kilodalton protein (EP300), tumor necrosis factor (TNF), and serine arginine-rich splicing factor 1 (SRSF1) in the black, red, and yellow modules.
FIGURE 2. (A) A scatterplot of Gene Significance (GS) for eosinophilic vs. Module Membership (MM) in the black module; 2 (B) A scatterplot of Gene Significance (GS) for eosinophilic vs. Module Membership (MM) in the yellow module.
FIGURE 3. (A) The top 12 genes in the black module; (B) The top 12 genes in the yellow module. The darker the color, the higher the connectivity of the gene with other genes. The lighter the color, the less connected the gene is to the other genes.
The hub genes in the red module were strongly associated with the age in Figure 1. After controlling for age, there was no statistically significant difference in the gene for the red module between the two groups of healthy controls and severe eosinophilic asthma (Supplementary Table S5). This indicated that hub genes in the red module might be closely related to the age of the patient.
CeRNA Network
EP300 was discovered to serve a key connecting function in the miRNA-mRNA network. Surprisingly, lncRNA would reduce the expression of EP300, which in turn would diminish the role of EP300 as a miRNA sponge and thus cause the down-regulation of mRNA (Figure 4). The hub genes were examined in the same way, and we discovered that CCT8 and miRNA-mRNA formed a circRNA-miRNA-mRNA network. CCT8 acted as a miRNA sponge and thus cause the up-regulation of mRNA (Figure 5). We found that miRNAs such as has-let-7t-5p and has-miR-98-5p bind to TNF, SPI1, and CCR7, limiting their expression in the red module (Supplementary Figure S7).
FIGURE 4. The lncRNA-EP300-miRNA-mRNA ceRNA network of the top 12 genes in the black module. lncRNA would reduce the expression of EP300, which in turn would diminish the role of EP300 as a miRNA sponge and thus cause the down-regulation of mRNA.
FIGURE 5. The CCT-8-miRNA-mRNA ceRNA network of the top 12 genes in the yellow module. CCT8 acted as a miRNA sponge and thus cause the up-regulation of mRNA.
Validation
Hub genes in the black module except FOXO3 were down-regulated in the GSE143303 (Supplementary Table S6), which showed hub genes other than FOXO3 might be involved in severe asthma but not particularly regulate severe eosinophilic asthma. FOXO3 may be a key gene in the regulation of eosinophilic asthma or severe eosinophilic asthma. The expression of hub genes other than FOXO3 was decreased in the GSE147878 (Supplementary Table S7), suggesting that hub genes other than FOXO3 might be key genes in the regulation of severe asthma. In addition, the expression of FOXO3 was discovered to be elevated in GSE137268 (Supplementary Table S8), which indicated that FOXO3 was not only implicated in severe eosinophilic asthma but also played a key role in eosinophilic asthma.
We discovered that the expression of hub genes in the yellow module was increased except for Ribosomal protein L17 (RPL17) and Heterogeneous nuclear ribonucleoprotein K (HNRNPK) in the GSE143303 (Supplementary Table S9), which indicated those hub genes other than RPL17 and HNRNPK play a role in severe asthma and do not particularly regulate severe eosinophilic asthma. RPL17 and HNRNPK may be key genes regulating eosinophilic asthma or severe eosinophilic asthma. The GSE147878 dataset further validated that genes other than RPL17 and HNRNPK were key genes regulating severe asthma (Supplementary Table S10). We did not find statistically significant differences in RPL17 and HNRNPK between the healthy control and non-severe eosinophilic asthma groups in the GSE137268 (Supplementary Table S11). Therefore, we speculated that RPL17 and HNRNPK might specifically regulate the specific phenotype of severe eosinophilic asthma.
The hub genes in the red module were all down-regulated in the GSE143303 dataset, but the difference was not statistically significant (Supplementary Table S12).
Functional Enrichment Analysis
The hub genes were mainly involved in the longevity regulating pathway and wnt signaling pathway in the black module. Further investigation revealed that FOXO3 was mainly enriched in pri-miRNA transcription by RNA polymerase II. FOXO3 and EP300 were also involved in chromatin binding, chromatin DNA binding, and beta-catenin binding. Interestingly, beta-catenin binding included steroid hormone mediated signaling pathway and cellular response to steroid hormone stimulus (Figure 6). The hub genes in the yellow module were mainly enriched in the spliceosome and ribosome. Further research demonstrated that CCT8 was mainly involved in unfolded protein binding. SRSF1 and RPL17 were related to the ribonucleoprotein complex. SRSF1 and HNRNPK were involved in RNA splicing, via transesterification reactions with bulged adenosine as a nucleophile. KEGG pathway analysis showed that SRSF1 and HNRNPK were mainly enriched in the Spliceosome, while RPL17 was enriched in the Ribosome (Figure 7). The genes in the red module were mainly enriched in viral protein interaction with cytokine and cytokine receptor and B cell receptor signaling pathway. The KEGG pathway of TNF was mainly enriched in Natural killer cell mediated cytotoxicity. GO enrichment analysis showed that TNF, SPI1, and CCR7 were mainly enriched in mononuclear cell migration, lymphocyte activation, regulation of leukocyte activation mononuclear cell migration, and positive regulation of leukocyte migration, negative regulation of immune effector process, and regulation of leukocyte migration (Supplementary Figure S8).
FIGURE 6. Functional enrichment of top 12 genes in the black module. The signal pathways were discovered by enrichment analysis into groups based on functional connection, the same group was colored the same color, and the labels of each group of the most essential terms were color-coded.
FIGURE 7. Functional enrichment of top 12 genes in the yellow module. The signal pathways were discovered by enrichment analysis into groups based on functional connection, the same group was colored the same color, and the labels of each group of the most essential terms were color-coded.
Discussion
Severe asthma is known to have a variety of phenotypes and endotypes, all of which have a significant impact on patients’ quality of life. The most prevalent phenotype of severe asthma is severe eosinophilic asthma. This is the first research to use WGCNA to create a ceRNA co-expression network in severe eosinophilic asthma. LncRNA-EP300 and CCT8 acting as miRNA sponges played a vital role in severe eosinophilic asthma. In addition, we found that EP300 and FOXO3 in the black module were regulated by steroid hormones, which were related to the medication used by the patient. The genes in the yellow module were the key genes in the treatment of severe eosinophilic asthma, in particular, RPL17 and HNRNPK might specifically regulate severe eosinophilic asthma. This suggests that genes regulating severe eosinophilic asthma and genes regulating severe asthma in the yellow module act in concert in the ceRNA network. We hope that our discoveries will help us better understand and treat severe eosinophilic asthma in the future.
In this study, the median dose of ICS in patients with severe eosinophilic asthma was 1600ug/day, and 5/23 patients received OCS. We found that FOXO3 was involved in cellular response to steroid hormone stimulus. However, the result of validation in GSE137268 showed that FOXO3 was up-regulated and statistically significant in non-severe eosinophilic asthma. FOXO3 is located on chromosome 6q21 and is protein-encoding gene-regulating aging, apoptosis, and tumor. Lützner et al. (2012) found that FOXO3 is a glucocorticoid receptor target that has two functional glucocorticoid responsive regions in its promoter. In the presence of glucocorticoids, FOXO3 stimulates its expression via a positive autoregulatory feedback loop. Yuan et al. (2020) showed that FOXO3 increased in mild asthma patients and decreased in severe patients. This indicated that FOXO3 might be a predictor of sensitivity to glucocorticoids in asthmatics. Non-severe eosinophilic asthma was still sensitive to glucocorticoid medication and had a good treatment outcome. In contrast, patients who were not sensitive to glucocorticoids and had poor treatment outcomes developed severe eosinophilic asthma. Therefore, FOXO3 may be a very significant biomarker for the progression of common eosinophilic asthma to severe eosinophilic asthma.
As we all know, histone modification is a crucial mechanism of epigenetic transcriptional control, histone acetylation (HAT) and histone deacetylases (HDACs) are enzymes that control the acetylation and deacetylation of histones (Cheng et al., 2019). EP300 is a well-known example of endogenous HAT. The levels of EP300 were substantially higher while the levels of histone deacetylases (HDACs) were much lower in asthma. In our study, we found that EP300 was involved in steroid hormone mediated signaling pathways and cellular response to steroid hormone stimulus. Ito et al. (2002) have previously shown that the increased expression of several inflammatory genes in asthma may be due to an increase in HAT activity. In addition, when asthmatic patients were given inhaled steroids, HAT activity was lowered to control levels.
SRSF1 is the archetype member of the SR protein family of splicing regulators (Das and Krainer, 2014). Recent studies demonstrated that SRSF1 is involved in inflammation. Fu et al. (2021) showed that SRSF1 expression was elevated in LPS-induced acute lung injury. However, it is unknown the role of SRSF1 that played in the development and progression of severe eosinophilic asthma. Only a few studies available so far showed that SRSF1 might be a potential biomarker for asthma (Maghsoudloo et al., 2020). The specific mechanism was to be further discovered.
CCT8 is a circRNA and is a member of the TCP-1 chaperone protein (CCT) family of genes. The chaperonin CCT8 controls proteostasis essential for T cell maturation, selection, and function (Oftedal et al., 2021). Current studies have found that CCT8 was overexpressed in cancer, but no study so far has reported any association of CCT8 with asthma. CircRNAs have attracted extensive attention in the pathogenesis of asthma in recent years (Huang et al., 2021). However, the specific roles of these circRNAs in severe asthma were not fully clear.
RPL17 is a member of the L22 family of ribosomal proteins and is the only ribosomal protein that interacts with all six structural domains of 23S rRNA and plays an important role in guiding the proper folding and conformation of 23S rRNA (Pool, 2009). RPL17 was mainly involved in viral gene expression, viral transcription, and cotranslational protein targeting the membrane. There is no evidence that RPL17 plays a role in asthma. However, some studies have explored the relationship between ribosomal proteins and asthma. Dong et al. (2017) showed that ribosomal protein S3 (RPS3) siRNA lessened HDM-induced airway mucus hypersecretion, cytokine production, and serum IgE elevation. Future studies are yet to be required to further explore the connection between RPL17 and severe eosinophilic asthma.
HNRNPK is a DNA/RNA-binding protein and regulates a wide range of biological processes and disease pathogenesis (Wang et al., 2020). The mechanism of HNRNPK involvement in severe eosinophilic asthma is currently unclear. Meng et al. (2021) found that HNRNPK acted downstream of TNFα-TNFR2 signaling, and was involved in the inflammatory response. Additional studies will be required to further explore the relationship between HNRNPK and severe eosinophilic asthma.
After adjusting for age, the expression of genes in the red module was not significantly different between the two groups. The genes in the red module were found to be mainly involved in inflammatory and immune response processes in the enrichment analysis. Increased levels of circulating cytokines and proinflammatory indicators including IL-6, IL-1, and TNF are linked to aging (Michaud et al., 2013). Previous studies have found that both TNF, SPI1, and CCR7 were up-regulated in severe asthma (Wittwer et al., 2006; Niessen et al., 2021; Wang et al., 2021). In addition, Xing et al. (2002) found dexamethasone resulted in significant inhibition of pro-inflammatory cytokines, T-cell stimulation, chemokines, and chemokine receptors. We speculated that the downregulation of hub genes in the red module may be due to differences in the age of the participants and the use of ICS and OCS.
Although the findings of our study have important clinical implications, the limitations must also be noted. First, we did not find a dataset that could validate the relationship between lncRNAs like MALAT1 and NEAT1 and severe eosinophilic asthma. Future studies are therefore expected to deeply explore the relationship between lncRNAs and severe eosinophilic asthma. Second, we could not obtain the dose of ICS and OCS taken by each patient with severe eosinophilic asthma, so we could not control the effect of drug use on the results. However, we also found the effect of drug use on the results by other methods such as dataset validation and enrichment analysis. At the same time, we also identified some genes in the black and red modules that were regulated when conventional drugs were used to treat severe eosinophilic asthma, and potential regulatory genes in the yellow modules of severe eosinophilic asthma that remain uncontrolled by drug therapy, which provided evidence for subsequent targeted and precise treatment of severe eosinophilic asthma.
Conclusion
In conclusion, our study reported that the ceRNA network played an essential and significant role in severe eosinophilic asthma. The hub genes in the yellow module were the key genes in the treatment of severe eosinophilic asthma. In particular, RPL17 and HNRNPK might specifically regulate severe eosinophilic asthma. Our results could be useful to provide potential immunotherapy targets and prognostic markers for severe eosinophilic asthmatics patients. Further mechanistic studies were required to validate and elucidate this result.
Data Availability Statement
Publicly available datasets were analyzed in this study. The names of the repository/repositories and accession number(s) can be found in the article/Supplementary Material.
Ethics Statement
Ethical review and approval was not required for the study on human participants in accordance with the local legislation and institutional requirements. The patients/participants provided their written informed consent to participate in this study.
Author Contributions
HW, ZZ, SY, and OC designed the study; YM, YJ, BM, and JG analyzed the data, HW wrote this article. All the authors read and approved the final manuscript.
Funding
This work was supported by the National Natural Science Foundation of China (8217091118); the Natural Science Foundation of Shandong Province (ZR2020MH006), and the Key Research and Development Program of Shandong Province (2019GSF108198).
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/fphar.2022.852536/full#supplementary-material
Abbreviations
ceRNA, competing endogenous RNA; circRNA, circular RNAs; EP300, E1A-associated 300-kilodalton protein; FDR, false discovery rate; GEO, gene expression omnibus; GINA, global initiative for asthma; GO, gene ontology; GS, gene significance; HAT, histone acetylation; HDACs, histone deacetylases; KEGG, kyoto encyclopedia of genes and genomes; lncRNAs, long non-coding RNA; MAD, median absolute deviation; miRNA, viral RNAs, microRNA; MM, module membership; MREs, microRNA response elements; mRNA, messenger RNA; ncRNA, non-coding RNA; PBMC, peripheral blood mononuclear cells; PPI, protein-protein interaction network; RNAInter, RNA interactome database; RPS3, ribosomal protein S3; SRSF1, serine arginine-rich splicing factor 1; STRING, search tool for the retrieval of interacting genes; TFs, transcription factors; WGCNA, weighted gene co-expression analysis.
References
Bindea, G., Mlecnik, B., Hackl, H., Charoentong, P., Tosolini, M., Kirilovsky, A., et al. (2009). ClueGO: a Cytoscape Plug-In to Decipher Functionally Grouped Gene Ontology and Pathway Annotation Networks. Bioinformatics 25 (8), 1091–1093. doi:10.1093/bioinformatics/btp101
Cheng, Q., Shang, Y., Huang, W., Zhang, Q., Li, X., and Zhou, Q. (2019). p300 Mediates the Histone Acetylation of ORMDL3 to Affect Airway Inflammation and Remodeling in Asthma. Int. Immunopharmacol. 76, 105885. doi:10.1016/j.intimp.2019.105885
Das, S., and Krainer, A. R. (2014). Emerging Functions of SRSF1, Splicing Factor and Oncoprotein, in RNA Metabolism and Cancer. Mol. Cancer Res. 12 (9), 1195–1204. doi:10.1158/1541-7786.MCR-14-0131
Dong, J., Liao, W., Peh, H. Y., Chan, T. K., Tan, W. S., Li, L., et al. (2017). Ribosomal Protein S3 Gene Silencing Protects against Experimental Allergic Asthma. Br. J. Pharmacol. 174 (7), 540–552. doi:10.1111/bph.13717
Fan, C. N., Ma, L., and Liu, N. (2018). Systematic Analysis of lncRNA-miRNA-mRNA Competing Endogenous RNA Network Identifies Four-lncRNA Signature as a Prognostic Biomarker for Breast Cancer. J. Transl. Med. 16, 264. doi:10.1186/s12967-018-1640-2
Fu, C., Hu, Y., Liu, J., and Yang, Y. (2021). The Splicing Factor SRSF1 Stabilizes the mRNA of TSLP to Enhance Acute Lung Injury. CELL Immunol. 362, 104299. doi:10.1016/j.cellimm.2021.104299
GEO Accession viewer (2021). “GEO Accession Viewer,” in 2021: NCBI's Gene Expression Omnibus (GEO) Is a Public Archive and Resource for Gene Expression Data.
Han, X., Huang, S., Xue, P., Fu, J., Liu, L., Zhang, C., et al. (2019). LncRNA PTPRE-AS1 Modulates M2 Macrophage Activation and Inflammatory Diseases by Epigenetic Promotion of PTPRE. Sci. Adv. 5 (12), eaax9230. doi:10.1126/sciadv.aax9230
Haselkorn, T., Fish, J. E., Zeiger, R. S., Szefler, S. J., Miller, D. P., Chipps, B. E., et al. (2009). Consistently Very Poorly Controlled Asthma, as Defined by the Impairment Domain of the Expert Panel Report 3 Guidelines, Increases Risk for Future Severe Asthma Exacerbations in the Epidemiology and Natural History of Asthma: Outcomes and Treatment Regimens (TENOR) Study. J. Allergy Clin. Immunol. 124 (5), 895–902.e1–4. doi:10.1016/j.jaci.2009.07.035
He, L. L., Xu, F., Zhan, X. Q., Chen, Z. H., and Shen, H. H. (2020). Identification of Critical Genes Associated with the Development of Asthma by Co-expression Modules Construction. Mol. Immunol. 123, 18–25. doi:10.1016/j.molimm.2020.01.015
Huang, Z., Cao, Y., Zhou, M., Qi, X., Fu, B., Mou, Y., et al. (2019). Hsa_circ_0005519 Increases IL-13/IL-6 by Regulating Hsa-Let-7a-5p in CD4+ T Cells to Affect Asthma. Clin. Exp. Allergy 49 (8), 1116–1127. doi:10.1111/cea.13445
Huang, Z., Fu, B., Qi, X., Xu, Y., Mou, Y., Zhou, M., et al. (2021). Diagnostic and Therapeutic Value of Hsa_circ_0002594 for T Helper 2-Mediated Allergic Asthma. Int. Arch. Allergy Immunol. 182 (5), 388–398. doi:10.1159/000511612
Ito, K., Caramori, G., Lim, S., Oates, T., Chung, K. F., Barnes, P. J., et al. (2002). Expression and Activity of Histone Deacetylases in Human Asthmatic Airways. Am. J. Respir. Crit. Care Med. 166 (3), 392–396. doi:10.1164/rccm.2110060
Jatakanon, A., Lim, S., and Barnes, P. J. (2000). Changes in Sputum Eosinophils Predict Loss of Asthma Control. Am. J. Respir. Crit. Care Med. 161, 64–72. doi:10.1164/ajrccm.161.1.9809100
Kang, J., Tang, Q., He, J., Li, L., Yang, N., Yu, S., et al. (2021). RNAInter v4.0: RNA Interactome Repository with Redefined Confidence Scoring System and Improved Accessibility. NUCLEIC ACIDS Res. 50 (D1), D326–D332. doi:10.1093/nar/gkab997
Kelly, R. S., Chawes, B. L., Blighe, K., Virkud, Y. V., Croteau-Chonka, D. C., McGeachie, M. J., et al. (2018). An Integrative Transcriptomic and Metabolomic Study of Lung Function in Children with Asthma. CHEST 154 (2), 335–348. doi:10.1016/j.chest.2018.05.038
Langfelder, P., and Horvath, S. (2008). WGCNA: an R Package for Weighted Correlation Network Analysis. BMC Bioinforma. 9, 559. doi:10.1186/1471-2105-9-559
Lee, S. H., Jeong, H. M., Choi, J. M., Cho, Y. C., Kim, T. S., Lee, K. Y., et al. (2009). Runx3 Inhibits IL-4 Production in T Cells via Physical Interaction with NFAT. Biochem. Biophys. Res. Commun. 381 (2), 214–217. doi:10.1016/j.bbrc.2009.02.026
Lemière, C., Ernst, P., Olivenstein, R., Yamauchi, Y., Govindaraju, K., Ludwig, M. S., et al. (2006). Airway Inflammation Assessed by Invasive and Noninvasive Means in Severe Asthma: Eosinophilic and Noneosinophilic Phenotypes. J. Allergy Clin. Immunol. 118, 1033–1039. doi:10.1016/j.jaci.2006.08.003
Liao, Y., Li, P., Wang, Y., Chen, H., Ning, S., and Su, D. (2020). Construction of Asthma Related Competing Endogenous RNA Network Revealed Novel Long Non-coding RNAs and Potential New Drugs. Respir. Res. 21 (1), 14. doi:10.1186/s12931-019-1257-x
Lützner, N., Kalbacher, H., Krones-Herzig, A., and Rösl, F. (2012). FOXO3 Is a Glucocorticoid Receptor Target and Regulates LKB1 and its Own Expression Based on Cellular AMP Levels via a Positive Autoregulatory Loop. PLOS ONE 7 (7), e42166. doi:10.1371/journal.pone.0042166
Maghsoudloo, M., Azimzadeh Jamalkandi, S., Najafi, A., and Masoudi-Nejad, A. (2020). An Efficient Hybrid Feature Selection Method to Identify Potential Biomarkers in Common Chronic Lung Inflammatory Diseases. GENOMICS 112 (5), 3284–3293. doi:10.1016/j.ygeno.2020.06.010
Meng, Y., Zhao, Q., An, L., Jiao, S., Li, R., Sang, Y., et al. (2021). A TNFR2-hnRNPK Axis Promotes Primary Liver Cancer Development via Activation of YAP Signaling in Hepatic Progenitor Cells. CANCER Res. 81 (11), 3036–3050. doi:10.1158/0008-5472.CAN-20-3175
Michaud, M., Balardy, L., Moulis, G., Gaudin, C., Peyrot, C., Vellas, B., et al. (2013). Proinflammatory Cytokines, Aging, and Age-Related Diseases. J. Am. Med. Dir. Assoc. 14 (12), 877–882. doi:10.1016/j.jamda.2013.05.009
Miller, J. A., Horvath, S., and Geschwind, D. H. (2010). Divergence of Human and Mouse Brain Transcriptome Highlights Alzheimer Disease Pathways. Proc. Natl. Acad. Sci. U. S. A. 107, 12698–12703. doi:10.1073/pnas.0914257107
Niessen, N. M., Gibson, P. G., Baines, K. J., Barker, D., Yang, I. A., Upham, J. W., et al. (2021). Sputum TNF Markers Are Increased in Neutrophilic and Severe Asthma and Are Reduced by Azithromycin Treatment. ALLERGY 76 (7), 2090–2101. doi:10.1111/all.14768
Oftedal, B. E., Maio, S., Handel, A. E., White, M. P. J., Howie, D., Davis, S., et al. (2021). The Chaperonin CCT8 Controls Proteostasis Essential for T Cell Maturation, Selection, and Function. Commun. Biol. 4 (1), 681. doi:10.1038/s42003-021-02203-0
Papi, A., Brightling, C., Pedersen, S. E., and Reddel, H. K. (2018). Asthma. LANCET 391 (10122), 783–800. doi:10.1016/S0140-6736(17)33311-1
Pool, M. R. (2009). A Trans-membrane Segment inside the Ribosome Exit Tunnel Triggers RAMP4 Recruitment to the Sec61p Translocase. J. Cell Biol. 185 (5), 889–902. doi:10.1083/jcb.200807066
Qiu, Y. Y., Wu, Y., Lin, M. J., Bian, T., Xiao, Y. L., and Qin, C. (2019). LncRNA-MEG3 Functions as a Competing Endogenous RNA to Regulate Treg/Th17 Balance in Patients with Asthma by Targeting microRNA-17/RORγt. Biomed. Pharmacother. 111, 386–394. doi:10.1016/j.biopha.2018.12.080
Sánchez-Ovando, S., Baines, K. J., Barker, D., Wark, P. A., and Simpson, J. L. (2020). Six Gene and TH2 Signature Expression in Endobronchial Biopsies of Participants with Asthma. Immun. Inflamm. Dis. 8 (1), 40–49. doi:10.1002/iid3.282
Sánchez-Ovando, S., Simpson, J. L., Barker, D., Baines, K. J., and Wark, P. A. B. (2021). Transcriptomics of Biopsies Identifies Novel Genes and Pathways Linked to Neutrophilic Inflammation in Severe Asthma. Clin. Exp. allergy J. Br. Soc. Allergy Clin. Immunol. 51 (10), 1279–1294.
Shang, Y., Sun, Y., Xu, J., Ge, X., Hu, Z., Xiao, J., et al. (2020). Exosomes from Mmu_circ_0001359-Modified ADSCs Attenuate Airway Remodeling by Enhancing FoxO1 Signaling-Mediated M2-like Macrophage Activation. Mol. Ther. Nucleic Acids 19, 951–960. doi:10.1016/j.omtn.2019.10.049
Thomson, D. W., and Dinger, M. E. (2016). Endogenous microRNA Sponges: Evidence and Controversy. Nat. Rev. Genet. 17 (5), 272–283. doi:10.1038/nrg.2016.20
Wang, C. H., Weng, C. M., Huang, T. T., Lee, M. J., Lo, C. Y., Chen, M. C., et al. (2021). Anti‐IgE Therapy Inhibits Chemotaxis, Proliferation and Transformation of Circulating Fibrocytes in Patients with Severe Allergic Asthma. Respirology 26 (9), 842–850. doi:10.1111/resp.14096
Wang, Z., Qiu, H., He, J., Liu, L., Xue, W., Fox, A., et al. (2020). The Emerging Roles of hnRNPK. J. CELL Physiol. 235 (3), 1995–2008. doi:10.1002/jcp.29186
Wittwer, J., Marti-Jaun, J., and Hersberger, M. (2006). Functional Polymorphism in ALOX15 Results in Increased Allele-specific Transcription in Macrophages through Binding of the Transcription Factor SPI1. Hum. Mutat. 27 (1), 78–87. doi:10.1002/humu.20273
Xing, N., L Maldonado, M. L. M. L., Bachman, L. A., McKean, D. J., Kumar, R., and Griffin, M. D. (2002). Distinctive Dendritic Cell Modulation by Vitamin D(3) and Glucocorticoid Pathways. Biochem. Biophys. Res. Commun. 297 (3), 645–652. doi:10.1016/s0006-291x(02)02262-3
Yuan, L., Wang, L., Du, X., Qin, L., Yang, M., Zhou, K., et al. (2020). The DNA Methylation of FOXO3 and TP53 as a Blood Biomarker of Late-Onset Asthma. J. Transl. Med. 18 (1), 467. doi:10.1186/s12967-020-02643-y
Zhang, Z., Wang, J., and Chen, O. (2021). Identification of Biomarkers and Pathogenesis in Severe Asthma by Coexpression Network Analysis. BMC Med. Genomics 14 (1), 51. doi:10.1186/s12920-021-00892-4
Keywords: severe eosinophilic asthma, competing endogenous RNA, co-expression network analysis, WGCNA (weighted gene co-expression network analyses), circRNA
Citation: Wang H, Zhang Z, Ma Y, Jia Y, Ma B, Gu J, Chen O and Yue S (2022) Construction of Severe Eosinophilic Asthma Related Competing Endogenous RNA Network by Weighted Gene Co-Expression Network Analysis. Front. Pharmacol. 13:852536. doi: 10.3389/fphar.2022.852536
Received: 11 January 2022; Accepted: 22 April 2022;
Published: 11 May 2022.
Edited by:
José Antonio Cornejo-García, Universidad de Málaga, SpainReviewed by:
Catalina Sanz, University of Salamanca, SpainAkul Singhania, La Jolla Institute for Immunology (LJI), United States
Copyright © 2022 Wang, Zhang, Ma, Jia, Ma, Gu, Chen and Yue. 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: Shouwei Yue, c2hvdXdlaXlAc2R1LmVkdS5jbg==; Ou Chen, Y2hlbm91QHNkdS5lZHUuY24=