Skip to main content

Identifying key genes and functionally enriched pathways in Th2-high asthma by weighted gene co-expression network analysis

Abstract

Background

Asthma is a chronic lung disease characterized by reversible inflammation of the airways. The imbalance of Th1/Th2 cells plays a significant role in the mechanisms of asthma. The aim of this study was to identify asthma-related key genes and functionally enriched pathways in a Th2-high group by using weighted gene coexpression network analysis (WGCNA).

Methods

The gene expression profiles of GSE4302, which included 42 asthma patients and 28 controls, were selected from the Gene Expression Omnibus (GEO). A gene network was constructed, and genes were classified into different modules using WGCNA. Gene ontology (GO) was performed to further explore the potential function of the genes in the most related module. In addition, the expression profile and diagnostic capacity (ROC curve) of hub genes of interest were verified by dataset GSE67472.

Results

In dataset GSE4302, subjects with asthma were divided into Th2-high and Th2-low groups according to the expression of the SERPINB2, POSTN and CLCA1 genes. A weighted gene coexpression network was constructed, and genes were classified into 7 modules. Among them, the red module was most closely associated with Th2-high asthma, which contained 60 genes. These genes were significantly enriched in different biological processes and molecular functions. A total of 8 hub genes (TPSB2, CPA3, ITLN1, CST1, SERPINB10, CEACAM5, CHD26 and P2RY14) were identified, and the expression levels of these genes (except TPSB2) were confirmed in dataset GSE67472. ROC curve analysis validated that the expression of these 8 genes exhibited excellent diagnostic efficiency for Th2-high asthma and Th2-low asthma.

Conclusions

The study provides a novel perspective on Th2-high asthma by WGCNA, and the hub genes and potential pathways involved may be beneficial for the diagnosis and management of Th2-high asthma.

Peer Review reports

Background

Asthma is a common chronic inflammatory disease of the lower airway characterized by airflow obstruction, bronchial hyperresponsiveness, and airway inflammation [1]. The CD4 T-helper cell type 2 (Th2)-mediated pathway in the airway epithelium is critical in allergic asthma [2]. Drugs targeting Th2 cytokines (IL-4, IL-5, and IL-13) have been developed and shown efficacy in recent studies [3, 4]. However, the etiology and progression of asthma are still under investigation.

Although multiple studies have shown differences in cytokine expression between Th2-high and Th2-low asthma [4, 5], it is not known whether the gene signature of bronchial airway epithelium can identify Th2-high and Th2-low endotypes. Weighted gene coexpression network analysis (WGCNA) is a useful tool to identify functional pathways and candidate biomarkers and has recently been used to explore the complex relationships between gene network signatures and complicated phenotypes [6]. Meanwhile, WGCNA has been successfully adopted to investigate the gene-network signature, coexpression modules, and hub genes in respiratory diseases, such as severe asthma [7,8,9] and chronic obstructive pulmonary disease (COPD) [10,11,12].

In this study, we analyzed gene expression in asthmatic patients and healthy control subjects and measured the expression of genes related to type 2 inflammation. We aimed to identify hub genes through WGCNA to further explore the potential mechanism underlying Th2-high asthma.

Materials and methods

Microarray data collection

We downloaded mRNA expression profiles GSE4302 and GSE67472 from the Gene Expression Omnibus (GEO, https://www.ncbi.nlm.nih.gov/geo). Specifically, airway epithelial brushings for microarray analysis of 42 subjects with asthma (not on inhaled steroids) and 28 healthy controls were obtained from GSE4302. Airway epithelial brushings for microarray analysis in GSE67472 were obtained from 62 subjects with mild-to-moderate asthma (not on inhaled steroids) and 43 healthy controls. Asthma subjects in the GSE4302 microarray dataset were divided into Th2-high and Th2-low asthma based on their expression of a three-gene signature of type 2 inflammation: SERPINB2, POSTN and CLCA1 [2, 12].

WGCNA

The coefficient of variation (CV) of each gene in the GSE4302 microarray dataset was calculated, and genes with CV > 5% (4133 genes) were used to construct a gene coexpression network using the WGCNA package (version 1.69) [6] in R. Firstly, the outlier samples were removed by hierarchical cluster analysis. Secondly, we constructed a scale-free network (R2 = 0.9) based on the criteria that soft-thresholding power β was set as 10 using the pickSoftThreshold function. Subsequently, the adjacency gene network was transformed into a topological overlap matrix (TOM), and the corresponding dissimilarity (1-TOM) was calculated. Lastly, dynamicTreeCut algorithm was used to identify modules, and 30 was chosen as the minimum number of genes in each module. A value of 0.25 was defined as the threshold for cut height to merge possible similar modules.

Identification of clinical significant modules

To further analyze the module, module eigengenes (ME) was calculated, which is the principal component of each gene module and could be considered as a representative of all genes in a given module. The correlation between genes and ME was defined as module membership (MM), and the correlation between the genes and clinical trait was defined as the gene significance (GS). The ME values were correlated with control, Th2-high and Th2-low groups by Pearson's correlation. Finally, the module highly correlated with clinical traits of Th2-high asthma was selected for further analysis.

Functional enrichment analysis

To obtain further insights into the function of the genes in the module most related to Th2-high asthma, Gene Ontology (GO) analysis was performed and visualized with the R package “Clusterprofiler” (version 3.18.0) [13]. The ontology contains three categories: molecular function (MF), biological process (BP), and cellular component (CC). A Q value < 0.05 was set as the cutoff.

Identification and validation of hub genes

In the module-trait correlation analysis, genes with absolute GS greater than 0.4 and high within-module connectivity of the modules (absolute MM > 0.8) were considered hub genes. The hub genes were further validated in datasets GSE67472. We screened the differentially expressed genes (DEGs) from Th2-high asthma vs. Th2-low asthma and Th2-high asthma vs. healthy controls using the “limma” R package. DEGs were considered genes with a false discovery rate (FDR) < 0.05 and absolute log2 fold change (FC) ≥ 1. Among the hub genes, eight genes (TPSB2, CPA3, ITLN1, CST1, SERPINB10, CEACAM5, CHD26, P2RY14) were further validated in another dataset GSE67472. The area under the receiver operating characteristics (ROC) curve (AUC) was calculated, and an ROC curve was plotted with the “ROCR” package to evaluate the capability of hub genes to distinguish Th2-high asthma and Th2-low asthma.

Results

Asthma phenotyping based on epithelial markers

Seventy subjects (28 healthy control subjects and 42 asthma subjects) in datasets GSE4302 underwent unsupervised hierarchical clustering based on the microarray expression levels of SERPINB2, POSTN and CLCA1 (Fig. 1). Approximately half of the asthma subjects (n = 22) with consistently high expression of a signature of Type 2 inflammation (SERPINB2, POSTN and CLCA1) were identified as Th2-high asthma, and 20 asthma subjects were identified as Th2-low asthma.

Fig. 1
figure 1

Heatmap showing unsupervised hierarchical clustering of POSTN, CLCA1, and SERPINB2 expression levels in bronchial epithelium. Red represents a high expression level, and green represents a low expression level

Weighted coexpression network construction

The GSE4302 samples were clustered using hierarchical agglomerative clustering with average linkage, and one outlier sample (GSM98201) was removed (Fig. 2a). In this study, GSE4302 was used to construct coexpression networks. The power of β = 10 (scale-free R2 = 0.9) was selected as the correlation coefficient threshold to ensure a scale-free network (Fig. 2b). A total of seven modules were identified through WGCNA (Fig. 2c).

Fig. 2
figure 2

a Clustering dendrogram of 69 samples. b Determination of soft-thresholding power in the WGCNA. Analysis of the scale-free fit index for various soft-thresholding powers (β) and analysis of the mean connectivity for various soft-thresholding powers. c Clustering dendrogram of differentially expressed genes based on a dissimilarity measure (1-TOM)

Key modules and hub genes identification

Module-trait relationship analyses showed that multiple modules were related to Th2-high asthma, and the red module, which included 60 genes, was most significantly associated with Th2-high asthma (Fig. 3a). This module was identified as the clinically significant module for further analysis. Figure 3b shows the significance of these genes in the red module, and 15 genes with high connectivity in the clinically significant module were identified as hub genes based on the cutoff criteria (|GS|> 0.4 and |MM|> 0.8). The detailed gene information of each module was listed in Additional file 1.

Fig. 3
figure 3

a Heatmap of the correlation between module and clinical traits (each cell contained the correlation coefficient and corresponding P value). b The gene significance for Th2-high asthma in the red module (one dot represents one gene in the red module)

Gene ontology of the key co-expression module

GO functional enrichment analysis of the genes in the clinically significant module showed that the genes were mainly enriched in BP and were involved in negative regulation of hydrolase activity, T cell activation and negative regulation of peptidase and endopeptidase activity (Fig. 4). The genes in the MF category were mainly enriched in endopeptidase/peptidase inhibitor activity and endopeptidase/peptidase regulator activity (Fig. 4). However, the genes in the CC group showed no enrichment results. These results indicated that these genes could be related to multiple biological pathways orchestrating Th2-high asthma development.

Fig. 4
figure 4

Gene ontology of red module genes including biological process (BP) analysis and molecular function (MF) analysis

Validation of hub genes

In dataset GSE4302, 27 DEGs were identified in the group of Th2-high asthma vs. Th2-low asthma (Fig. 5a), and 38 DEGs were identified in the Th2-high asthma vs. healthy group (Fig. 5a). The expression of eight hub genes (CST1, P2RY14, SERPINB10, CEACAM5, CDH26, CPA3, TPSB2, and ITLN1) was significantly upregulated in Th2-high asthmatic subjects compared to Th2-low asthmatic subjects and healthy subjects (Fig. 5a). Moreover, the expression levels of the seven hub genes were investigated in another dataset, GSE67472 (the gene TPSB2 expression data were not found in GSE67472). The seven genes were also significantly upregulated in the airway epithelium of Th2-high asthma patients compared to Th2-low asthma patients and controls in dataset GSE67472 (Fig. 5b). In addition, the ROC curve indicated that these hub genes exhibited excellent diagnostic efficiency for Th2-high asthma and Th2-low asthma in datasets GSE4302 and GSE67472 (Fig. 5c).

Fig. 5
figure 5

Validation of hub genes in datasets GSE4302 and GSE67472. a Volcano plot for differential gene expression in dataset GSE4302. b The expression levels of CST1, P2RY14, SERPINB10, CEACAM5, CDH26, CPA3 and ITLN1 were significantly upregulated in Th2-high asthma patients in dataset GSE67472. The Kruskal–Wallis test was used to evaluate the statistical significance of differences. c ROC curve of hub genes in dataset GSE4302 and dataset GSE67472

Discussion

Bronchial asthma a highly heterogenous disease. Based on the characteristics of inflammation, bronchial asthma can be identified into two endotypes, type 2 asthma and nontype 2 asthma [14,15,16]. Th2 cells and cytokines (an imbalance of Th1/Th2) play a major role in type 2 inflammation, which is considered a major pathophysiological mechanism of asthma [2, 4]. Several reports have identified a typical Th2 feature characterized by elevated Th2-type inflammation in the lower airway [7, 17]. In another study, gene expression in patients with mild asthma and healthy control subjects found that Th2 signature genes in lower airway epithelial cells were associated with atopy and eosinophilic airway inflammation [18, 19].

There is growing interest in seeking biomarkers to distinguish Th2-high and Th2-low phenotypes and to predict responsiveness to treatment [20]. At present, several biomarkers, such as fractional exhaled nitric oxide (FeNO), serum IgE, blood or sputum eosinophils and serum periostin, have been adopted [21]. Bhakta et al. reported that the three-gene mean of periostin, CLCA1 and serpinB2 in airway epithelial brushings identifies Th2-high and Th2-low populations and predicts FEV1 improvement with inhaled corticosteroid (ICS) [22].

WGCNA can provide a deep understanding of gene networks, hub genes and pathogenesis associated with asthma. In this study, we speculated that assessing gene expression in airway epithelial cells using WGCNA may identify the gene-network signature, provide potential biomarkers of Th2-high asthma and further predict treatment responses.

Seven modules were identified as significantly correlated with asthma status. The red module observed in this study was positively correlated with asthma traits and regarded as the most important module in Th2-high asthma pathogenesis. It was further analyzed using a GO enrichment method and the results showed it to be mainly enriched in endopeptidase and peptidase inhibitor/regulator activity. Then, key genes of TPSB2, CPA3, ITLN1, CST1, SERPINB10, CEACAM5, CHD26, and P2RY14 were identified.

The above eight hub genes were validated in another dataset GSE67472. We found that their expression levels (except TPSB2, which had not been studied in dataset GSE67472) were also significantly upregulated in Th2-high asthma patients compared to controls in dataset GSE67472. The diagnostic efficiencies of hub genes were evaluated by ROC curve analysis. Notably, the AUC of ROC curves were all more than 0.750 for the eight hub genes, indicating a high diagnostic value. This suggested that these hub genes probably participate in the pathogenesis and immune imbalance and could be used as potential biomarkers for the diagnosis and therapeutic targets of Th2-high asthma.

Although the mechanisms underlying asthma are unclear, some of these hub genes are under investigation [23,24,25,26,27,28,29]. For example, CST1 is a member of the type 2 cystatin (CST) superfamily, which is known to inhibit the proteolytic activities of cysteine proteases. Increased CPA3 expression in intraepithelial mast cells among Th2-high asthma patients was observed [23]. Yan et al. discussed the coexpressed DEGs of CST1 and CPA3 and found that there were significant correlations of the CST1 and CPA3 genes with novel biomarkers involved in the comorbidity of rhinitis and asthma [24]. In addition, Natasha A Winter et al. found that the gene expression of TPSAB1/TPSB2 and CPA3 correlated with sputum mast cells/basophils, and in severe asthma, TPSB2 and CPA3 were associated with eosinophilic airway inflammation and blood eosinophils [25].

The hub gene ITLN1, which is expressed in airway epithelial cells and involved in inflammatory pathways downstream of IL-13, was found to be significantly upregulated in the bronchial brushings of patients with asthma compared with controls [26, 27]. It has been reported that ITLN1 is a biomarker associated with disease susceptibility in asthma [28]. A study on ERPINB10-knockdown mice found that these groups of mice had diminished numbers of Th2 cells and were more susceptible to apoptosis, indicating that SERPINB10 may contribute to allergic inflammation and the Th2 response of asthma by inhibiting the apoptosis of Th2 cells [29].

In general, our study identified some hub genes through WGCNA to further explore the potential mechanism underlying Th2-high asthma. This may provide a better molecular characterization of asthma and consequent establishment of new pharmacological targets. Furthermore, the Th2-high subgroup of asthma is corticosteroid sensitive and benefits from biologic agents [4, 30]. In this respect, future research on the function of these genes not only helps promote our understanding of Th2-high and Th2-low endotypes of asthma but may also help to guide treatment strategies and predict response to treatment.

However, there were some limitations to our study. On one hand, the study only presented data mining and data analysis, and did not perform experiments to further study the exact mechanism of the identified hub genes in asthma. On the other hand, the analysis was based on limited genetic data with small sample sizes. Thus, the expression of these identified hub genes still needs to be validated in more datasets, and the exact mechanisms also need more exploration.

Conclusion

In summary, our study performed WGCNA to provide a framework of coexpression gene modules of Th2-high asthma and led to the identification of some hub genes that may help us to better understand the underlying pathogenesis and guide treatment decisions for the Th2-high subgroup of asthma. The exact molecular mechanisms of the hub genes and functional pathways in Th2 asthma still need to be further explored.

Availability of data and materials

Our data can be found in the Gene Expression Omnibus (GEO, https://www.ncbi.nlm.nih.gov/geo/, GSE4302 and GSE67472) database.

Abbreviations

WGCNA:

Weighted gene coexpression network analysis

GEO:

Gene Expression Omnibus

GO:

Gene ontology

Th2:

T-helper cell type 2

COPD:

Chronic obstructive pulmonary disease

CV:

Coefficient of variation

MF:

Molecular function

BP:

Biological process

CC:

Cellular component

DEGs:

Differentially expressed genes

FDR:

False discovery rate

FC:

Fold change

FeNO:

Fractional exhaled nitric oxide

ICS:

Inhaled corticosteroid

References

  1. Liu Z, Li M, Fang X, Shen L, Yao W, Fang Z, et al. Identification of surrogate prognostic biomarkers for allergic asthma in nasal epithelial brushing samples by WGCNA. J Cell Biochem. 2018;120(4):5137–50.

    Article  Google Scholar 

  2. Woodruff PG, Modrek B, Choy DF, Jia G, Abbas AR, Ellwanger A, et al. T-helper type 2–driven inflammation defines major subphenotypes of asthma. Am J Respir Crit Care Med. 2009;180(5):388–95.

    CAS  Article  Google Scholar 

  3. Busse WW, Kraft M, Rabe KF, Deniz Y, Rowe PJ, Ruddy M, et al. Understanding the key issues in the treatment of uncontrolled persistent asthma with type 2 inflammation. Eur Respir J. 2021;58(2):2003393. https://doi.org/10.1183/13993003.03393-2020.

  4. Matucci A, Vivarelli E, Nencini F, Maggi E, Vultaggio A. Strategies targeting type 2 inflammation: from monoclonal antibodies to JAK-inhibitors. Biomedicines. 2021;9(10):1497.

    CAS  Article  Google Scholar 

  5. Peters MC, Mekonnen ZK, Yuan S, Bhakta NR, Woodruff PG, Fahy JV. Measures of gene expression in sputum cells can identify TH2-high and TH2-low subtypes of asthma. J Allergy Clin Immunol. 2014;133(2):388–94.

    CAS  Article  Google Scholar 

  6. Langfelder P, Horvath S. WGCNA: an R package for weighted correlation network analysis. BMC Bioinform. 2008;9(1):1–13.

    Article  Google Scholar 

  7. Modena BD, Bleecker ER, Busse WW, Erzurum SC, Gaston BM, Jarjour NN, et al. Gene expression correlated with severe asthma characteristics reveals heterogeneous mechanisms of severe disease. Am J Respir Crit Care Med. 2017;195(11):1449–63.

    CAS  Article  Google Scholar 

  8. Zhang Z, Wang J, Chen O. Identification of biomarkers and pathogenesis in severe asthma by coexpression network analysis. BMC Med Genomics. 2021;14(1):1–9.

    CAS  Article  Google Scholar 

  9. He L-L, Xu F, Zhan X-Q, Chen Z-H, Shen H-H. Identification of critical genes associated with the development of asthma by co-expression modules construction. Mol Immunol. 2020;123:18–25.

    CAS  Article  Google Scholar 

  10. Morrow JD, Qiu W, Chhabra D, Rennard SI, Belloni P, Belousov A, et al. Identifying a gene expression signature of frequent COPD exacerbations in peripheral blood using network methods. BMC Med Genomics. 2015;8(1):1–11.

    Article  Google Scholar 

  11. Yi G, Liang M, Li M, Fang X, Liu J, Lai Y, et al. A large lung gene expression study identifying IL1B as a novel player in airway inflammation in COPD airway epithelial cells. Inflamm Res. 2018;67(6):539–51.

    CAS  Article  Google Scholar 

  12. Christenson SA, Steiling K, van den Berge M, Hijazi K, Hiemstra PS, Postma DS, et al. Asthma–COPD overlap. Clinical relevance of genomic signatures of type 2 inflammation in chronic obstructive pulmonary disease. Am J Respir Crit Care Med. 2015;191(7):758–66.

    CAS  Article  Google Scholar 

  13. Yu G, Wang L-G, Han Y, He Q-Y. clusterProfiler: an R package for comparing biological themes among gene clusters. OMICS J Integr Biol. 2012;16(5):284–7.

    CAS  Article  Google Scholar 

  14. Ricciardolo FLM, Carriero V, Bertolini F. Which therapy for non-type(T)2/T2-low asthma. J Pers Med. 2021;12(1):10.

    Article  Google Scholar 

  15. Santus P, Saad M, Damiani G, Patella V, Radovanovic D. Current and future targeted therapies for severe asthma: managing treatment with biologics based on phenotypes and biomarkers. Pharmacol Res. 2019;146:104296.

    CAS  Article  Google Scholar 

  16. Lötvall J, Akdis CA, Bacharier LB, Bjermer L, Casale TB, Custovic A, et al. Asthma endotypes: a new approach to classification of disease entities within the asthma syndrome. J Allergy Clin Immunol. 2011;127(2):355–60.

    Article  Google Scholar 

  17. Bhakta NR, Christenson SA, Nerella S, Solberg OD, Nguyen CP, Choy DF, et al. IFN-stimulated gene expression, type 2 inflammation, and endoplasmic reticulum stress in asthma. Am J Respir Crit Care Med. 2018;197(3):313–24.

    CAS  Article  Google Scholar 

  18. Woodruff G, Subtypes P. of asthma defined by epithelial cell expression of messenger RNA and MicroRNA. Ann Am Thorac Soc. 2013;10(Supplement):S186–9.

    Article  Google Scholar 

  19. Woodruff PG, Boushey HA, Dolganov GM, Barker CS, Yang YH, Donnelly S, et al. Genome-wide profiling identifies epithelial cell genes associated with asthma and with treatment response to corticosteroids. Proc Natl Acad Sci. 2007;104(40):15858–63.

    CAS  Article  Google Scholar 

  20. Robinson D, Humbert M, Buhl R, Cruz AA, Inoue H, Korom S, et al. Revisiting Type 2-high and Type 2-low airway inflammation in asthma: current knowledge and therapeutic implications. Clin Exp Allergy. 2017;47(2):161–75.

    CAS  Article  Google Scholar 

  21. Diamant Z, Tufvesson E, Bjermer L. Which biomarkers are effective for identifying Th2-driven inflammation in asthma? Curr Allergy Asthma Rep. 2013;13(5):477–86.

    CAS  Article  Google Scholar 

  22. Bhakta NR, Solberg OD, Nguyen CP, Nguyen CN, Arron JR, Fahy JV, et al. A qPCR-based metric of Th2 airway inflammation in asthma. Clin Transl Allergy. 2013;3(1):24.

    CAS  Article  Google Scholar 

  23. Dougherty RH, Sidhu SS, Raman K, Solon M, Solberg OD, Caughey GH, et al. Accumulation of intraepithelial mast cells with a unique protease phenotype in TH2-high asthma. J Allergy Clin Immunol. 2010;125(5):1046–53.

    CAS  Article  Google Scholar 

  24. Yan Z, Liu L, Jiao L, Wen X, Liu J, Wang N. Bioinformatics analysis and identification of underlying biomarkers potentially linking allergic rhinitis and asthma. Med Sci Monit. 2020;26:e924934.

    CAS  PubMed  PubMed Central  Google Scholar 

  25. Winter NA, Qin L, Gibson PG, McDonald VM, Baines KJ, Faulkner J, et al. Sputum mast cell/basophil gene expression relates to inflammatory and clinical features of severe asthma. J Allergy Clin Immunol. 2021;148(2):428–38.

    CAS  Article  Google Scholar 

  26. Kuperman DA, Lewis CC, Woodruff PG, Rodriguez MW, Yang YH, Dolganov GM, et al. Dissecting asthma using focused transgenic modeling and functional genomics. J Allergy Clin Immunol. 2005;116(2):305–11.

    CAS  Article  Google Scholar 

  27. Watanabe T, Chibana K, Shiobara T, Tei R, Koike R, Nakamura Y, et al. Expression of intelectin-1 in bronchial epithelial cells of asthma is correlated with T-helper 2 (Type-2) related parameters and its function. Allergy Asthma Clin Immunol. 2017;13(1):1–11.

    Article  Google Scholar 

  28. Pemberton AD, Rose-Zerilli MJ, Holloway JW, Gray RD, Holgate ST. A single-nucleotide polymorphism in intelectin 1 is associated with increased asthma risk. J Allergy Clin Immunol. 2008;122(5):1033–4.

    CAS  Article  Google Scholar 

  29. Mo Y, Ye L, Cai H, Zhu G, Wang J, Zhu M, et al. SERPINB10 contributes to asthma by inhibiting the apoptosis of allergenic Th2 cells. Respir Res. 2021;22(1):1–12.

    Article  Google Scholar 

  30. Fahy JV. Type 2 inflammation in asthma—present in most, absent in many. Nat Rev Immunol. 2014;15(1):57–65.

    Article  Google Scholar 

Download references

Acknowledgements

Not applicable.

Funding

No funding was received for this article.

Author information

Affiliations

Authors

Contributions

YC, LC and YW are responsible for the design of the study and prepared for the manuscript. LL, LY and XP collected the data and revised the manuscript. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Lina Chen.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1

. The detailed gene information of each module.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visithttp://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Cao, Y., Wu, Y., Lin, L. et al. Identifying key genes and functionally enriched pathways in Th2-high asthma by weighted gene co-expression network analysis. BMC Med Genomics 15, 110 (2022). https://doi.org/10.1186/s12920-022-01241-9

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12920-022-01241-9

Keywords

  • Asthma
  • Weighted gene coexpression network analysis (WGCNA)
  • Hub gene
  • Gene ontology