- 1Department of Respiratory and Critical Care Medicine, The First Affiliated Hospital of Guangxi Medical University, Nanning, China
- 2Department of Respiratory and Critical Care Medicine, The First Affiliated Hospital of Henan University of Science and Technology, Luoyang, China
Background and Objective: Corticosteroid resistance is a major barrier to chronic obstructive pulmonary disease (COPD), but the exact mechanism of corticosteroid resistance in COPD has been less well studied.
Methods: The microarray dataset GSE11906, which includes genomic and clinical data on COPD, was downloaded from the Gene Expression Omnibus (GEO) database, and the differentially expressed genes (DEGs) were identified using R software. Gene set enrichment analysis (GSEA) and Kyoto Encyclopedia of Genes (KEGG) were utilized to enrich and analyze the gene cohort related to the response to steroid hormones, respectively. The Connectivity Map (CMap) database was used to screen corticosteroid resistance-related drugs that might exert a potential therapeutic effect. STRING was used to construct a protein-protein interaction (PPI) network of the gene cohort, and the CytoHubba plug-in of Cytoscape was used to screen the hub genes in the PPI network. The expression levels of hub genes in cigarette smoke extract (CSE)-stimulated bronchial epithelial cells were assayed by quantitative real-time PCR and western blotting.
Results: Twenty-one genes were found to be correlated with the response to steroid hormones. In the CMap database, 32 small-molecule compounds that might exert a therapeutic effect on corticosteroid resistance in COPD were identified. Nine hub genes were extracted from the PPI network. The expression levels of the BMP4, FOS, FN1, EGFR, and SPP1 proteins were consistent with the microarray data obtained from molecular biology experiments. Scopoletin significantly restrained the increases in the levels of AKR1C3, ALDH3A1, FN1 and reversed the decreases of phosphorylated GR and HDAC2 caused by CSE exposure.
Conclusion: The BMP4, FOS, FN1, EGFR, and SPP1 genes are closely correlated with CSE-induced glucocorticoid resistance in airway epithelial cells. Scopoletin may be a potential drug for the treatment of glucocorticoid resistance caused by CSE.
Introduction
Chronic obstructive pulmonary disease (COPD) is estimated to become the third leading cause of death worldwide by 2020. Environmental exposure can cause COPD, and the primary factor is cigarette smoking (López-Campos et al., 2016). Chronic inflammation in the airways and lung parenchyma is a characteristic of COPD, and various inflammatory cytokines are upregulated in the lungs (López-Campos et al., 2016; Barnes, 2019). Although glucocorticoids are the front-line anti-inflammatory therapies for chronic inflammation, most patients with COPD have a poor response (Rossios et al., 2012). The molecular mechanisms of steroid resistance are due to the decreased expression of the glucocorticoid receptor (GR), which is the receptor for corticosteroids (Barnes, 2017). Several molecular pathways are involved in steroid resistance and lead to the downregulation of GR, and these pathways include the PIK3CD/Akt, AP-1, P38 MAPK, and JNK pathways (Barnes, 2010), whose activity and expression are upregulated in the lungs of smokers or patients with COPD (Kirkham and Barnes, 2013; Barnes et al., 2015; Fischer et al., 2015). Histone deacetylase-2 (HDAC2) can be recruited by corticosteroids, which results in the effective repression of activated inflammatory genes. Cigarette smoking markedly suppresses the activity and expression of HDAC2 induced by oxidative/nitrative stress and results in inflammation showing resistance to glucocorticoids (Barnes, 2009; Li et al., 2012). However, the effects of the current treatments for corticosteroid resistance are poor. Therefore, finding new targets for treating corticosteroid resistance has potential clinical value.
With the recent advancement of gene chips, genome-wide expression profiles can be obtained using microarray platforms and analyzed using R language and unbiased bioinformatic tools to identify novel hub genes (Sun et al., 2019). Exporting microarray data to modern pathway profiling software could reveal novel targets for treatment (Wei et al., 2019; Xie et al., 2019). Gene set enrichment analysis (GSEA) can help analyze and interpret microarray platform data (Powers et al., 2018). The power of GSEA, which summarizes microarray data of gene expression changes into gene sets, provides a list of genes that are significantly enriched based on their correlation with biological functions and biological processes. This method can enrich genes that are significantly associated with a pathway and phenotype by enabling the detection of gene sets, and its practicality for many applications has been demonstrated (Subramanian et al., 2007; Powers et al., 2018). This method has been successfully used to identify key pathways involved in pancreatic carcinoma (Zhou et al., 2018), reveal seven pathways associated with methyladenosine-related genes in hepatocellular carcinoma (Zhou et al., 2019), and elucidate the role of the TGF-β signaling pathway in preeclampsia (Wu et al., 2019). These results show the potential use of GSEA for studying COPD-related genes and extracting genes involved in corticosteroid resistance from the whole genome.
Here, we applied a bioinformatics method to reveal novel genes associated with corticosteroid resistance in patients with COPD. The “limma” package was used to identify the differentially expressed genes (DEGs), and a GSEA was performed to reveal corticosteroid resistance-related genes. A protein-protein interaction (PPI) network was established using String (https://string-db.org/), and the hub genes were extracted using the Cytoscape software. Furthermore, the corticosteroid resistance-related genes were verified by quantitative real-time polymerase chain reaction (qPCR) and western blot (WB) analyses.
Materials and Methods
Microarray Data Collection
The raw microarray dataset of COPD (GSE11906), which was obtained using the Affymetrix human genome U133 plus 2.0 array, was downloaded from the Gene Expression Omnibus (GEO) database (https://www.ncbi.nlm.nih.gov/geo/). The tissue source of these samples was the bronchial airway epithelium. The mRNA expression data were normalized by log2 transformation to obtain the COPD/normal ratio in the GSE11906 dataset.
GSEA
The “limma” package of R software (version 3.5.1) was used to screen GSE11906 to identify the DEGs, and mRNAs with a P-value < 0.05 were considered to be genes with significant differential expression. For the identification of genes associated with the response to steroid hormones, the DEGs were analyzed by a GSEA using Gene Ontology (GO) gene sets accessed from MSigDB. Java GSEA software was used for GSEA, the gene sets from COPD vs normal samples were compared, and the gene sets with an FDR <0.1 and a nominal p < 0.05 were regarded as significantly enriched. The distinct expression profiles of genes associated with the response to steroid hormones are shown as a heatmap.
Functional Enrichment Analysis
To further research the functional implications of the DEGs associated with the response to steroid hormones, we performed Kyoto Encyclopedia of Genes and Genomes (KEGG) functional enrichment analyses. The results from the KEGG enrichment analyses were based on the criterion P-value < 0.05. The significant enrichment results were visualized using R software.
PPI of Corticosteroid Resistance-Related Genes and Identification of Hub Genes
After screening for genes associated with the response to steroid hormones based on module connectivity, the genes that were found to be enriched using the GSEA software were uploaded to the Search Tool for the Retrieval of Interacting Genes (STRING) database to study the interactions among these genes, as well as NR3C1/GR and HDAC2. The interaction data for the above-mentioned genes were downloaded from STRING and analyzed using Cytoscape (version: 3.5.1) to construct a PPI network. Furthermore, Cytoscape software was used to analyze the PPI network to identify the hub genes (Yi et al., 2018; Zhou et al., 2018). CytoHubba in Cytoscape was used to select the hub genes in the PPIs, and the analysis was performed using the default parameters.
Screening of Potential Drugs for Corticosteroid Resistance Treatment Based on a Connectivity Map
The CMap database was used to screen small molecules with therapeutic potential in corticosteroid resistance. The DEGs associated with corticosteroid resistance were uploaded to the CMap database and analyzed, and a negative connectivity score suggested that the drug might reverse the signature biology of human cells and show its potential therapeutic value (Liu et al., 2020). The small molecules were identified based on the following criteria: P-value < 0.05 and a negative connectivity score.
Validation of Key Genes Through Biological Experiments
BEAS-2B Cell Culture and Treatment
The human bronchial epithelial cell line BEAS-2B (ATCC® CRL-9609™, Shanghai Cell Bank, Chinese Academy of Sciences) was maintained in RPMI 1640 medium (Gibco, New York, NY, United States) supplemented with 10% fetal bovine serum (Gibco) at 37°C under 5% CO2 conditions. BEAS-2B cells were treated with cigarette smoke extract (CSE) for 24 h and collected for the extraction of RNA and proteins for qPCR and WB assays, respectively.
CSE Preparation
CSE was prepared according to the protocol from our previous study (Li et al., 2012). The smoke of 10 unfiltered cigarettes was slowly dissolved in 20 ml of RPMI-1640 medium. The pH of the CSE solution was adjusted to 7.4, and the CSE solution was filtered through a 0.22-μm leach to remove bacteria. The optical density was measured at a wavelength of 320 nm. The CSE solution was used in the experiments within 1 h of extraction. A CSE concentration of 0.1% was used for the study based on the results obtained using the Cell Counting Kit-8 (CCK-8) assay.
Validation of the Results by qPCR and WB Analyses
Total RNA from the CSE-treated and control cells was extracted using the TRIzol reagent (Invitrogen, Carlsbad, CA, United States). The primer sequences are listed in Table 1, and the total RNA quality and quantity were analyzed with a spectrophotometer (Thermo, United States). RNA was reverse transcribed into cDNA using the PrimeScriptTM RT reagent kit with gDNA Eraser (TaKaRa, Dalian, China), and mRNA was prepared using kits. qPCR was performed using SYBR® Premix Ex TaqTM II (Takara, Dalian, China). The procedures were performed according to the manufacturers’ instructions. The protein concentration in the lysate of BEAS-2B cells stimulated with 0.1% CSE for 24 h was determined by the bicinchoninic acid method. Twenty micrograms of protein was electrophoresed on a sodium dodecyl sulfate polyacrylamide gel and then electroblotted onto polyvinylidene difluoride membranes. After blocking with 5% skim milk, the membranes were incubated with the corresponding primary antibodies overnight in a shaker at 4°C. The following primary antibodies were used in this study: anti-BMP4 (1:1000; A11405, ABclonal), SPP1 (1:1000; A11405, ABclonal), ALDH3A1 (1:1000; A13275, ABclonal), AKR1C3 (1:1000; A1781, ABclonal), c-Fos (1:1000; #2250, Cell Signaling Technology), FN1 (1:1000; #26836, Cell Signaling Technology), EGFR (1:1000; #4267, Cell Signaling Technology), GR (1:1000; #12041, Cell Signaling Technology), phospho (P)-GR (1:1000; #4161, Cell Signaling Technology), HDAC2 (1:1000; #2540, Cell Signaling Technology), and GAPDH (1:4,000; AF7021, Affinity). The expression of proteins was detected using the Odyssey system after incubation with a secondary antibody (1:20,000; #5151, Cell Signaling Technology). Scopoletin (S2500, Sigma-Aldrich) was considered a compound that could potentially alleviate steroid resistance according to the results of CMap database analysis. WB analysis was used to confirm the influence of scopoletin on the hub genes.
Statistical Analysis
The qPCR and WB data are presented as the means ±SDs. Three repeated independent experiments were performed from each group. Statistical comparisons were performed using the SPSS 18.0 software (Chicago, IL, United States). A two-tailed Student’s t-test was used for comparison between two groups. One-way analysis of variance was conducted to evaluate significant differences among multiple groups. The results with p < 0.05 were considered statistically significant.
Results
Identification of Differentially Expressed Genes
Thirty-three COPD samples from long-term smokers and 72 normal samples from nonsmokers were included in this study. The main clinical features of the participants are listed in Table 2. A total of 335 DEGs, including 187 upregulated and 148 downregulated genes, were obtained from GSE11906 using the limma package. Volcano plots of the DEGs are shown in Figure 1 (fold change ≥1.5 and P-value ≤ 0.05).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/672065/fphar-12-672065-HTML/image_m/fphar-12-672065-t002.jpg)
TABLE 2. Clinical information of the study participants. GOLD: The Global Initiative for Chronic Obstructive Lung Disease.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/672065/fphar-12-672065-HTML/image_m/fphar-12-672065-g001.jpg)
FIGURE 1. Volcano plot of DEG expression in the microarray. The red points in the volcano plot are the mRNAs showing significantly higher expression, and the green points show the mRNAs with significantly lower expression. The differences are set as a fold change ≥1.5 and a P-value ≤ 0.05. FC, fold change.
GSEA of DEGs
GSEA of the DEGs was conducted using GSEA software based on the following filter conditions: P-value < 0.05 and FDR q-value < 0.05. The GSEA plots of the genes enriched in the response to steroid hormones are shown in Figure 2. The 21 genes related to the response to steroid hormones that were found to be enriched by the GSEA are shown in Table 3, and the expression of these 21 genes is shown in the heatmap diagram presented in Figure 3.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/672065/fphar-12-672065-HTML/image_m/fphar-12-672065-g003.jpg)
FIGURE 3. Heatmap of DEGs related to the response to steroid hormones. The blue coded bar above the heat map represents the normal samples and the red coded bar represents the copd samples. Red dots represent upregulated DEGs, and blue dots represent downregulated DEGs.
KEGG Enrichment Analysis
To study the functional implications of the 21 above-mentioned genes, we performed a KEGG functional enrichment analysis. The results showed that these genes were significantly enriched in the immune system, signal transduction, cytokine signaling in the immune system, fluid shear stress, atherosclerosis and pathways in cancer (Figure 4).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/672065/fphar-12-672065-HTML/image_m/fphar-12-672065-g004.jpg)
FIGURE 4. KEGG enrichment analysis of the 21 genes related to the response to steroid hormones. The outer circle represents the number of genes enriched in the corresponding pathway. Red dots indicate upregulated DEGs and blue dots indicate downregulated DEGs. The inner circle represents the Z-score. Red indicates more enrichment of upregulated DEGs, and purple represents more enrichment of downregulated DEGs.
Construction of the PPI Network and Identification of Hub Genes
The 21 above-mentioned genes, as well as NR3C1 and HDAC2, were uploaded to STRING to study their interactions, and PPIs were plotted using Cytoscape based on the interaction data obtained from STRING (Figure 5A). In the PPI network of the 21 genes related to the response to steroid hormones in COPD and NR3C1, the top nine hub genes with the highest criticality were EGFR, FOS, SPP1, FN1, CCL2, IL1B, NQ O 1, BMP4, and C3 (Figure 5B). Therefore, we speculated that these nine hub genes might play an important role in corticosteroid resistance in smoke-induced COPD.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/672065/fphar-12-672065-HTML/image_m/fphar-12-672065-g005.jpg)
FIGURE 5. PPI network of genes associated with glucocorticoid resistance. (A): PPI network of 21 genes associated with glucocorticoid resistance. The node size indicates the level of the gene in the PPI network, and the node color depth indicates the level of the clustering coefficient. The thickness of the edge lines indicates the size of the coexpression relationship between the genes, and the color depth of the edge lines indicates the levels of the interaction coefficients between the genes. (B): PPI of nine hub genes in the 21-gene cohort associated with glucocorticoid resistance. A darker node color indicates that the gene plays a more critical role in the network.
Verification of Hub Genes
To verify the hub genes extracted from the PPI network, we determined their mRNA and protein expression levels in CSE-treated human bronchial epithelial cells using qPCR and western blot analyses. In addition, we verified the expression of ALDH3A1 and AKR1C3, which were significantly upregulated in patients with COPD and belonged to the core enrichment genes according to the microarray analysis. The qPCR results indicated that the expression of all hub genes, as well as ALDH3A1 and AKR1C3, was consistent with the data from the selected microarray datasets (Figure 6). The WB data showed that the expression levels of the C-FOS, EGFR, GR, P-GR, and HDAC2 proteins were significantly downregulated in the CSE group. The expression levels of the BMP4, FN1, SPP1, ALDH3A1, and AKR1C3 proteins in BEAS-2B cells increased after stimulation with 0.1% CSE (p < 0.05; Figure 7). No significant differences in the expression of the CCL2, IL1B, C3, and NQO1 proteins were found between the CSE and control groups.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/672065/fphar-12-672065-HTML/image_m/fphar-12-672065-g006.jpg)
FIGURE 6. Relative mRNA expression levels of the hub genes in 0.1% CSE-stimulated Beas-2B cells (*p < 0.05, **p < 0.01, ***p < 0.001).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/672065/fphar-12-672065-HTML/image_m/fphar-12-672065-g007.jpg)
FIGURE 7. Relative protein expression levels determined by western blot analysis. (A) BMP4, (B) AKR1C3, (C) ALDH3A1, (D) C-FOS, (E) FN1, (F) EGFR, (G) SPP1, (H) GR, (I) P-GR and (J) HDAC2 (*p < 0.05, **p < 0.01).
Prediction and Verification of Potential Agents for the Treatment of CSE-Induced Corticosteroid Resistance
Small-molecule compounds that might alleviate corticosteroid resistance were predicted through CMap analysis. The results suggest that 32 compounds might exert a therapeutic effect on corticosteroid resistance induced by oxidative stress in COPD (Table 4). The most promising small-molecule compounds are scopoletin, harmol, AG-012559, AH-23848, perhexiline, and fluorometholone. We chose scopoletin as a potential agent for its ability to inhibit the activation of nuclear factor-κB (NF-κB) against COPD-associated oxidative stress (Leema and Tamizhselvi, 2018). A scopoletin concentration of 100 μM was chosen as the final intervention concentration according to the results of our CCK-8 assay, and the compound was added to the medium 2 h before CSE treatment. The results indicated that scopoletin significantly restrained the increases in the levels of AKR1C3, ALDH3A1 and FN1 caused by CSE exposure. In addition, scopoletin obviously reversed the decreases in P-GR and HDAC2 levels (Figure 8).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/672065/fphar-12-672065-HTML/image_m/fphar-12-672065-t004.jpg)
TABLE 4. Small-molecule compounds identified as potential drugs for corticosteroid resistance treatment by CMap analysis.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/672065/fphar-12-672065-HTML/image_m/fphar-12-672065-g008.jpg)
FIGURE 8. Scopoletin contributes to the recovery of the indicators of CSE-induced glucocorticoid resistance. (A) AKR1C3, (B) ALDH3A1, (C) FN1, (D) HDAC2, (F) P-GR (*p < 0.05, **p < 0.01, ***p < 0.001).
Discussion
The imbalance of pulmonary oxidants and antioxidants caused by long-term exposure to cigarette smoke is thought to be the main driving mechanism of COPD (Kirkham and Barnes, 2013). Several studies have demonstrated that oxidative stress and inflammation in COPD are not only local but also interdependent and closely related processes in systemic inflammation in the lungs (Sundar et al., 2013). Cigarette smoke-mediated oxidative stress and inflammation induce a number of kinase signaling pathways that are considered to be associated with steroid resistance (Liao et al., 2020). Bronchial/lung epithelial cells are the first barrier against inhaled cigarette smoke. Oxidative stress induced by cigarette smoke promotes the release of various cytokines and chemokines in epithelial cells. Epithelial damage and repair play a central role in the pathogenesis of COPD (Heijink et al., 2014). Previous studies have shown that the stimulation of lung epithelial cells with cigarette smoke can result in corticosteroid insensitivity. An important molecular mechanism leading to glucocorticoid resistance has been shown to be a result of decreases in the expression and activity of GR and HDAC2 caused by cigarette smoke exposure (Barnes, 2017).
Our study extracted DEGs from the GEO database of COPD at the whole gene level, and 21 enriched genes correlated with the response to steroid hormones were identified from a GSEA of these DEGs. This gene cohort included NQ O 1, ALDH3A1, AKR1C3, NR2F1, ABHD2, NR0B1, SPP1, BMP4, CCL2, DEFB1, FN1, CA2, NR4A3, IL1B, NR2F6, FOSB, FOS, EGFR, GHR, C3, and NTS. The nine core genes in the gene cohort, which were considered to make a major contribution to the enrichment score of the gene set, were NQO1, ALDH3A1, AKR1C3, NR2F1, ABHD2, NR0B1, SPP1, BMP4, and CCL2. We also analyzed the 21-gene cohort through the KEGG database. The results showed that the genes in the cohort were enriched in the immune system, signal transduction, cytokine signaling in the immune system, fluid shear stress, atherosclerosis, and pathways in cancer. These results indicate that the functions of the 21 genes are closely correlated with the immune response. Furthermore, we constructed a PPI network to assess the relationships among the 21 genes identified in the cohort, as well as HDAC2 and GR. The following hub genes of the 21-gene PPI network were extracted using the Cytoscape plug-in of CytoHubba: FN1, EGFR, CCL2, IL1B, C3, SPP1, NQO1, FOS, and BMP4.
A previous study demonstrated that miR-183 regulates LPS-induced oxidative stress in rat hippocampal neurons by regulating the expression of FN1 to reduce oxidative damage (Xie et al., 2020). Nicotine and cigarette smoke obviously stimulate the expression of FN1 in lung fibroblasts and mouse lung tissues (Roman et al., 2004; Zhao et al., 2019). EGFR plays an important role in cancer progression and regulates inflammation and oxidative stress (Fang et al., 2016). EGFR was shown to be activated by the meprin alpha metalloproteinase to mediate oxidative stress in macrophages (Wang et al., 2017). Additionally, EGFR suppresses the activation of transcription factors, such as NF-κB, and induces proinflammatory gene transcription in macrophages. CCL2 and IL1B are common cytokines involved in chronic inflammation in COPD (Sethi et al., 2012; Barnes, 2016). Cigarette smoke can induce macrophages to release inflammatory mediators, including CXCL1, CXCL8, and CCL2. The increased release of IL1 leads to an increase in the number of macrophages in the lungs of smokers and patients with COPD (Barnes, 2016). Cigarette smoke-induced reactive oxygen species (ROS) can induce IL1B generation by activating Toll-like receptors and inducing immune responses (Zuo et al., 2014). C3 activation is closely related to oxidative stress in mice. Mice with transient focal cerebral ischemia with lower ROS levels show a better neurological outcome, and the inhibition of C3 activation in mice leads to an anti-inflammatory effect (Yang et al., 2013). SPP1 is a molecule with a complex function that decreases the oxidative stress effect and has been implicated in inflammation in the kidney (Trostel et al., 2018). Additionally, SPP1 might exert proinflammatory effects in vascular endothelial cells induced by cigarette smoke and might contribute to inflammation in smokers (Bishop et al., 2012).
NQO1 protein is highly expressed in various cells, such as vascular endothelial cells, epithelial cells, and adipocytes. A number of studies have shown that NQO1 exerts antioxidant effects by reducing the level of ROS (Jo et al., 2016). However, no significant upregulation was observed at the protein level in CSE-treated BEAS-2B cells. FOS is a subunit of AP-1, which is a proinflammatory transcription factor sensitive to oxidative stress that interacts with GR proteins (Frenkel et al., 2015). The overexpression of AP-1 can prevents the binding of GR to glucocorticoid response elements and other proinflammatory DNA-binding transcription factors via PPIs with GR (Pujolsa et al., 2009; Jacques et al., 2010). Previous studies have shown that glucocorticoid resistance in asthma is correlated with the overexpression of AP-1 (Jacques et al., 2010). BMP4 is a new proinflammatory marker that stimulates the production of ROS and induces monocyte adhesion (Tian et al., 2012). BMP4 was found to be upregulated in the airway epithelium of asymptomatic smokers and patients with COPD. The accumulation of BMP4 aggravates airway remodeling by changing the phenotypes of basal stem/progenitor cells (Zuo et al., 2019). AKR1C3 is believed to be closely related to the metabolism of steroids. It can be activated by oxidative stress factors and leads to radiotherapy resistance in lung cancer tissues (Xie et al., 2013). ALDH3A1 was found to be upregulated in lung tissue after exposure to carcinogenic aldehydes (Patel et al., 2008).
To study whether the nine hub genes are altered by CSE stimulation in BEAS-2B cells, we determined the expression of these genes using qPCR and WB assays. Our results demonstrated that the changes in the mRNA expression of the nine genes were consistent with the results from the analysis of the microarray data. The expression levels of the BMP4, FOS, FN1, EGFR, SPP1, ALDH3A1, and AKR1C3 proteins were consistent with the mRNA results. The results demonstrate that CSE could induce glucocorticoid resistance by regulating these genes. The inconsistency between the mRNA and protein levels of CCL2, IL1B, C3, and NQO1 may be related to post-transcriptional modification or post-translational protein processing, specific mechanism needs further study.
CMap is a database covering diseases, gene expression patterns, and interactions between small-molecule compounds (Gao et al., 2019; Wei et al., 2020). This database is an effective genome-based tool for screening new chemopreventive drugs and is commonly used to find small-molecule compounds that exert therapeutic effects on certain diseases (Wei et al., 2020; Yu et al., 2020). We used CMap for drug discovery and identified small-molecule compounds as potential drugs for corticosteroid resistance treatment. Our results showed that 32 small-molecule compounds might exert potential therapeutic effects. Some compounds have been found to exert antioxidative stress effects in other fields. For example, in yeast cells, harmol and harmalol have a significant protective effect against oxidative stress induced by H2O2 and paraquat. The antigenotoxic and antimutagenic effects of harmol result from its ability to scavenge hydroxyl radicals (Moura et al., 2007). Scopoletin efficiently quenches oxidative stress by stabilizing the nuclear factor erythroid 2-related factor 2 (Nrf2)/antioxidant responsive element (ARE) pathway through enhancement of the nuclear translocation and phosphorylation of Nrf2 (Narasimhan et al., 2019). We chose scopoletin as a candidate compound that might alleviate steroid resistance based on the results of the CMap database analysis. The results of treatment indicated that scopoletin significantly restrained the increases in the AKR1C3, ALDH3A1, and FN1 levels caused by CSE exposure and reversed the decreases in the P-GR and HDAC2 levels. Thus, scopoletin may contribute to the reduction of CSE-induced glucocorticoid resistance in airway epithelial cells. However, the specific molecular regulatory mechanism and how these genes interact to regulate corticosteroid resistance still need to be clarified in future experiments.
In summary, we extracted and validated seven genes that were correlated with corticosteroid resistance and found that they interacted with GR and HDAC2. However, the study has some limitations. First, we have no direct evidence showing the mutual regulatory relationships among the hub genes, and further studies are needed to validate these relationships. Second, animal model validation is needed for future mechanistic research.
Conclusion
In this study, we extracted 21 DEGs by integrating the microarray data of COPD and identified a set of enriched genes related to the response to steroid hormones. It was found that the BMP4, FOS, FN1, EGFR, and SPP1 genes were closely correlated with CSE-induced glucocorticoid resistance in airway epithelial cells. In addition, scopoletin was found to be a potential drug for the treatment of glucocorticoid resistance caused by CSE. We hope this research can contribute to the development of corticosteroid resistance therapy for COPD in the future.
Data Availability Statement
The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding author.
Author Contributions
ZH conceived the ideas. JD designed and supervised the research. GP and NM performed the research and wrote the manuscript. LG and FC collected and analyzed the data. JB revised the manuscript. All authors read and approved the final version of the manuscript.
Funding
The study was supported by the National Natural Science Foundation of China (Grant No. 81860010) and Guangxi Medical High-level Backbone Personnel Training “139” Plan Training Project (Grant No. G201903008, G202002002).
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.
References
Barnes, P. J., Burney, P. G., Silverman, E. K., Celli, B. R., Vestbo, J., Wedzicha, J. A., et al. (2015). Chronic Obstructive Pulmonary Disease. Nat. Rev. Dis. Primers 1, 15076. doi:10.1038/nrdp.2015.76
Barnes, P. J. (2017). Glucocorticosteroids. Handb Exp. Pharmacol. 237, 93–115. doi:10.1007/164_2016_62
Barnes, P. J. (2009). Histone Deacetylase-2 and Airway Disease. Ther. Adv. Respir. Dis. 3, 235–243. doi:10.1177/1753465809348648
Barnes, P. J. (2016). Inflammatory Mechanisms in Patients with Chronic Obstructive Pulmonary Disease. J. Allergy Clin. Immunol. 138, 16–27. doi:10.1016/j.jaci.2016.05.011
Barnes, P. J. (2010). Mechanisms and Resistance in Glucocorticoid Control of Inflammation. J. Steroid Biochem. Mol. Biol. 120, 76–85. doi:10.1016/j.jsbmb.2010.02.018
Bishop, E., Theophilus, E. H., and Fearon, I. M. (2012). In Vitro and Clinical Studies Examining the Expression of Osteopontin in Cigarette Smoke-Exposed Endothelial Cells and Cigarette Smokers. BMC Cardiovasc. Disord. 12, 75. doi:10.1186/1471-2261-12-75
Fang, Q., Zou, C., Zhong, P., Lin, F., Li, W., Wang, L., et al. (2016). EGFR Mediates Hyperlipidemia-Induced Renal Injury via Regulating Inflammation and Oxidative Stress: the Detrimental Role and Mechanism of EGFR Activation. Oncotarget 7, 24361–24373. doi:10.18632/oncotarget.8222
Fischer, B. M., Voynow, J. A., and Ghio, A. J. (2015). COPD: Balancing Oxidants and Antioxidants. Int. J. Chron. Obstruct Pulmon Dis. 10, 261–276. doi:10.2147/COPD.S42414
Frenkel, B., White, W., and Tuckermann, J. (2015). Glucocorticoid-Induced Osteoporosis. Adv. Exp. Med. Biol. 872, 179–215. doi:10.1007/978-1-4939-2895-8_8
Gao, Y., Kim, S., Lee, Y. I., and Lee, J. (2019). Cellular Stress-Modulating Drugs Can Potentially Be Identified by In Silico Screening with Connectivity Map (CMap). Int. J. Mol. Sci. 20, 5601. doi:10.3390/ijms20225601
Heijink, I., van Oosterhout, A., Kliphuis, N., Jonker, M., Hoffmann, R., Telenga, E., et al. (2014). Oxidant-induced Corticosteroid Unresponsiveness in Human Bronchial Epithelial Cells. Thorax 69, 5–13. doi:10.1136/thoraxjnl-2013-203520
Jacques, E., Semlali, A., Boulet, L. P., and Chakir, J. (2010). AP-1 Overexpression Impairs Corticosteroid Inhibition of Collagen Production by Fibroblasts Isolated from Asthmatic Subjects. Am. J. Physiol. Lung Cel Mol PhysiolLung Cell. Mol. Physiol. 299, L281–L287. doi:10.1152/ajplung.00360.2009
Jo, H. S., Kim, D. S., Ahn, E. H., Kim, D. W., Shin, M. J., Cho, S. B., et al. (2016). Protective Effects of Tat-NQO1 against Oxidative Stress-Induced HT-22 Cell Damage, and Ischemic Injury in Animals. BMB Rep. 49, 617–622. doi:10.5483/bmbrep.2016.49.11.117
Kirkham, P. A., and Barnes, P. J. (2013). Oxidative Stress in COPD. Chest 144, 266–273. doi:10.1378/chest.12-2664
Leema, G., and Tamizhselvi, R. (2018). Protective Effect of Scopoletin against Cerulein-Induced Acute Pancreatitis and Associated Lung Injury in Mice. Pancreas 47, 577–585. doi:10.1097/MPA.0000000000001034
Li, M., Zhong, X., He, Z., Wen, M., Li, J., Peng, X., et al. (2012). Effect of Erythromycin on Cigarette-Induced Histone Deacetylase Protein Expression and Nuclear Factor-Κb Activity in Human Macrophages In Vitro. Int. Immunopharmacol 12, 643–650. doi:10.1016/j.intimp.2011.12.022
Liao, W., Lim, A. Y. H., Tan, W. S. D., Abisheganaden, J., and Wong, W. S. F. (2020). Restoration of HDAC2 and Nrf2 by Andrographolide Overcomes Corticosteroid Resistance in Chronic Obstructive Pulmonary Disease. Br. J. Pharmacol. 177, 3662–3673. doi:10.1111/bph.15080
Liu, J., Feng, M., Li, S., Nie, S., Wang, H., Wu, S., et al. (2020). Identification of Molecular Markers Associated with the Progression and Prognosis of Endometrial Cancer: a Bioinformatic Study. Cancer Cel Int 20, 59. doi:10.1186/s12935-020-1140-3
López-Campos, J. L., Tan, W., and Soriano, J. B. (2016). Global burden of COPD. Respirology 21, 14–23. doi:10.1111/resp.12660
Moura, D. J., Richter, M. F., Boeira, J. M., Pêgas Henriques, J. A., and Saffi, J. (2007). Antioxidant Properties of Beta-Carboline Alkaloids Are Related to Their Antimutagenic and Antigenotoxic Activities. Mutagenesis 22, 293–302. doi:10.1093/mutage/gem016
Narasimhan, K. K. S., Jayakumar, D., Velusamy, P., Srinivasan, A., Mohan, T., Ravi, D. B., et al. (2019). Morinda citrifolia and its Active Principle Scopoletin Mitigate Protein Aggregation and Neuronal Apoptosis through Augmenting the DJ-1/Nrf2/ARE Signaling Pathway. Oxid Med. Cel Longev 2019, 2761041. doi:10.1155/2019/2761041
Patel, M., Lu, L., Zander, D. S., Sreerama, L., Coco, D., and Moreb, J. S. (2008). ALDH1A1 and ALDH3A1 Expression in Lung Cancers: Correlation with Histologic Type and Potential Precursors. Lung Cancer 59, 340–349. doi:10.1016/j.lungcan.2007.08.033
Powers, R. K., Goodspeed, A., Pielke-Lombardo, H., Tan, A. C., and Costello, J. C. (2018). GSEA-InContext: Identifying Novel and Common Patterns in Expression Experiments. Bioinformatics 34, i555–i564. doi:10.1093/bioinformatics/bty271
Pujolsa, L., Mullol, J., and Picado, C. (2009). Glucocorticoid Receptor in Human Respiratory Epithelial Cells. Neuroimmunomodulation 16, 290–299. doi:10.1159/000216187
Roman, J., Ritzenthaler, J. D., Gil-Acosta, A., Rivera, H. N., and Roser-Page, S. (2004). Nicotine and Fibronectin Expression in Lung Fibroblasts: Implications for Tobacco-Related Lung Tissue Remodeling. FASEB J. 18, 1436–1438. doi:10.1096/fj.03-0826fje
Rossios, C., To, Y., Osoata, G., Ito, M., Barnes, P. J., and Ito, K. (2012). Corticosteroid Insensitivity Is Reversed by Formoterol via Phosphoinositide-3-Kinase Inhibition. Br. J. Pharmacol. 167, 775–786. doi:10.1111/j.1476-5381.2012.01864.x
Sethi, S., Mahler, D. A., Marcus, P., Owen, C. A., Yawn, B., and Rennard, S. (2012). Inflammation in COPD: Implications for Management. Am. J. Med. 125, 1162–1170. doi:10.1016/j.amjmed.2012.06.024
Subramanian, A., Kuehn, H., Gould, J., Tamayo, P., and Mesirov, J. P. (2007). GSEA-P: a Desktop Application for Gene Set Enrichment Analysis. Bioinformatics 23, 3251–3253. doi:10.1093/bioinformatics/btm369
Sun, G., Li, Y., Peng, Y., Lu, D., Zhang, F., Cui, X., et al. (2019). Identification of a Five-Gene Signature with Prognostic Value in Colorectal Cancer. J. Cel Physiol. 234, 3829–3836. doi:10.1002/jcp.27154
Sundar, I. K., Yao, H., and Rahman, I. (2013). Oxidative Stress and Chromatin Remodeling in Chronic Obstructive Pulmonary Disease and Smoking-Related Diseases. Antioxid. Redox Signal. 18, 1956–1971. doi:10.1089/ars.2012.4863
Tian, X. Y., Yung, L. H., Wong, W. T., Liu, J., Leung, F. P., Liu, L., et al. (2012). Bone Morphogenic Protein-4 Induces Endothelial Cell Apoptosis through Oxidative Stress-dependent p38MAPK and JNK Pathway. J. Mol. Cel Cardiol 52, 237–244. doi:10.1016/j.yjmcc.2011.10.013
Trostel, J., Truong, L. D., Roncal-Jimenez, C., Miyazaki, M., Miyazaki-Anzai, S., Kuwabara, M., et al. (2018). Different Effects of Global Osteopontin and Macrophage Osteopontin in Glomerular Injury. Am. J. Physiol. Ren. Physiol 315, F759–f768. doi:10.1152/ajprenal.00458.2017
Wang, L., Huang, Z., Huang, W., Chen, X., Shan, P., Zhong, P., et al. (2017). Inhibition of Epidermal Growth Factor Receptor Attenuates Atherosclerosis via Decreasing Inflammation and Oxidative Stress. Sci. Rep. 8, 45917. doi:10.1038/srep45917
Wei, C., Liang, Q., Li, X., Li, H., Liu, Y., Huang, X., et al. (2019). Bioinformatics Profiling Utilized a Nine Immune-Related Long Noncoding RNA Signature as a Prognostic Target for Pancreatic Cancer. J. Cel Biochem 120, 14916–14927. doi:10.1002/jcb.28754
Wei, G., Sun, H., Liu, J. L., Dong, K., Liu, J., and Zhang, M. (2020). Indirubin, a Small Molecular Deriving from Connectivity Map (CMAP) Screening, Ameliorates Obesity-Induced Metabolic Dysfunction by Enhancing Brown Adipose Thermogenesis and white Adipose browning. Nutr. Metab. (Lond) 17, 21. doi:10.1186/s12986-020-00440-4
Wu, K., Liu, F., Wu, W., Chen, Y., and Zhang, W. (2019). Bioinformatics Approach Reveals the Critical Role of TGF-β Signaling Pathway in Pre-eclampsia Development. Eur. J. Obstet. Gynecol. Reprod. Biol. 240, 130–138. doi:10.1016/j.ejogrb.2019.06.034
Xie, L., Yu, J., Guo, W., Wei, L., Liu, Y., Wang, X., et al. (2013). Aldo-keto Reductase 1C3 May Be a New Radioresistance Marker in Non-small-cell Lung Cancer. Cancer Gene Ther. 20, 260–266. doi:10.1038/cgt.2013.15
Xie, Q., Feng, L., Zhang, C., and Gan, L. (2020). MicroRNA-183 Regulates Lipopolysaccharide-Induced Oxidative Stress of Hippocampal Neurons by Targeting the Fibronectin 1 Gene. Neurosci. Lett. 725, 134902. doi:10.1016/j.neulet.2020.134902
Xie, Y., Lu, W., Wang, S., Tang, X., Tang, H., Zhou, Y., et al. (2019). Validation of the 12-gene Predictive Signature for Adjuvant Chemotherapy Response in Lung Cancer. Clin. Cancer Res. 25, 150–157. doi:10.1158/1078-0432.CCR-17-2543
Yang, J., Ahn, H. N., Chang, M., Narasimhan, P., Chan, P. H., and Song, Y. S. (2013). Complement Component 3 Inhibition by an Antioxidant Is Neuroprotective after Cerebral Ischemia and Reperfusion in Mice. J. Neurochem. 124, 523–535. doi:10.1111/jnc.12111
Yi, G., Liang, M., Li, M., Fang, X., Liu, J., Lai, Y., et al. (2018). A Large Lung Gene Expression Study Identifying IL1B as a Novel Player in Airway Inflammation in COPD Airway Epithelial Cells. Inflamm. Res. 67, 539–551. doi:10.1007/s00011-018-1145-8
Yu, L., Xu, F., and Gao, L. (2020). Predict New Therapeutic Drugs for Hepatocellular Carcinoma Based on Gene Mutation and Expression. Front. Bioeng. Biotechnol. 8, 8. doi:10.3389/fbioe.2020.00008
Zhao, J., Li, M., Wang, Z., Chen, J., Zhao, J., Xu, Y., et al. (2019). Role of PM2.5 in the Development and Progression of COPD and its Mechanisms. Respir. Res. 20, 120. doi:10.1186/s12931-019-1081-3
Zhou, Y., Yin, Z., Hou, B., Yu, M., Chen, R., Jin, H., et al. (2019). Expression Profiles and Prognostic Significance of RNA N6-Methyladenosine-Related Genes in Patients with Hepatocellular Carcinoma: Evidence from Independent Datasets. Cancer Manag. Res. 11, 3921–3931. doi:10.2147/CMAR.S191565
Zhou, Z., Cheng, Y., Jiang, Y., Liu, S., Zhang, M., Liu, J., et al. (2018). Ten Hub Genes Associated with Progression and Prognosis of Pancreatic Carcinoma Identified by Co-expression Analysis. Int. J. Biol. Sci. 14, 124–136. doi:10.7150/ijbs.22619
Zuo, L., He, F., Sergakis, G. G., Koozehchian, M. S., Stimpfl, J. N., Rong, Y., et al. (2014). Interrelated Role of Cigarette Smoking, Oxidative Stress, and Immune Response in COPD and Corresponding Treatments. Am. J. Physiol. Lung Cel Mol PhysiolLung Cell. Mol. Physiol. 307, L205–L218. doi:10.1152/ajplung.00330.2013
Keywords: microarray analysis, bioinformatics analysis, corticosteroid resistance, hub genes, human airway epithelial cells, scopoletin
Citation: Pei G, Ma N, Chen F, Guo L, Bai J, Deng J and He Z (2021) Screening and Identification of Hub Genes in the Corticosteroid Resistance Network in Human Airway Epithelial Cells via Microarray Analysis. Front. Pharmacol. 12:672065. doi: 10.3389/fphar.2021.672065
Received: 28 February 2021; Accepted: 27 September 2021;
Published: 11 October 2021.
Edited by:
Ali Önder Yildirim, Helmholtz-Gemeinschaft Deutscher Forschungszentren (HZ), GermanyReviewed by:
Henriette Uhlenhaut, Helmholtz-Gemeinschaft Deutscher Forschungszentren (HZ), GermanyGizem Günes, Helmholtz Center Munich, Germany
Copyright © 2021 Pei, Ma, Chen, Guo, Bai, Deng and He. 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: Zhiyi He, zhiyi-river@163.com
†These authors have contributed equally to this work