Gongqi Chen1,2
Dian Chen1,2
Yuchen Feng1,2
Wenliang Wu1,2
Jiali Gao1,2Chenli Chang1,2Shengchong Chen1,2Guohua Zhen1,2*- 1Division of Respiratory and Critical Care Medicine, Department of Internal Medicine, Tongji Hospital, Tongji Medical College, Huazhong University of Science and Technology, Wuhan, China
- 2Key Laboratory of Respiratory Diseases, National Health Commission of People’s Republic of China, National Clinical Research Center for Respiratory Diseases, Wuhan, China
Background: Asthma is a heterogeneous disease with different subtypes including eosinophilic asthma (EA) and neutrophilic asthma (NA). However, the mechanisms underlying the difference between the two subtypes are not fully understood.
Methods: Microarray datasets (GSE45111 and GSE137268) were acquired from Gene Expression Omnibus (GEO) database. Differentially expressed genes (DEGs) in induced sputum between EA (n = 24) and NA (n = 15) were identified by “Limma” package. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses and Gene set enrichment analysis (GSEA) were used to explore potential signaling pathways. Weighted gene co-expression network analysis (WGCNA) were performed to identify the key genes that were strongly associated with EA and NA.
Results: A total of 282 DEGs were identified in induced sputum of NA patients compared with EA patients. In GO and KEGG pathway analyses, DEGs were enriched in positive regulation of cytokine production, and cytokine-cytokine receptor interaction. The results of GSEA showed that ribosome, Parkinson’s disease, and oxidative phosphorylation were positively correlated with EA while toll-like receptor signaling pathway, primary immunodeficiency, and NOD-like receptor signaling pathway were positively correlated with NA. Using WGCNA analysis, we identified a set of genes significantly associated NA including IRFG, IRF1, STAT1, IFIH1, IFIT3, GBP1, GBP5, IFIT2, CXCL9, and CXCL11.
Conclusion: We identified potential signaling pathways and key genes involved in the pathogenesis of the asthma subsets, especially in neutrophilic asthma.
Introduction
Asthma is a chronic airway inflammatory disease characterized by airway hyperresponsiveness, reversible airflow limitation, mucus overproduction and airway wall remodeling (King-Biggs, 2019; Gao et al., 2020). Asthma affects more than 300 million people worldwide with approximately 250,000 deaths per year (Scherzer and Grayson, 2018; Ijaz et al., 2019).
According to the differential counts of sputum inflammatory cells, asthma can be divided into four subtypes: eosinophilic asthma (EA), neutrophilic asthma (NA), mixed granulocytic asthma, and paucigranulocytic asthma (Tliba and Panettieri, 2019). EA and NA have attracted our attention on account of their differences in clinical features and therapeutic effects. EA is defined by having more than 3% of eosinophils in the sputum and with a variable severity (Erle and Sheppard, 2014; Schleich et al., 2016). Type 2 immune response plays an important role in EA. NA is defined by having more than 61% of sputum neutrophils and mostly occurs in more severe asthma (Pelaia et al., 2015; Taylor et al., 2018). Inhaled corticosteroids (ICS) is effective with airway eosinophilic inflammation, whereas neutrophilic asthma has poor response for ICS (Green et al., 2002; Jayaram et al., 2006). Although Th2 and Th17 signaling pathways are implicated in EA and NA, respectively, the underlying mechanism distinguishes these two subtypes remains unclear (Boonpiyathad et al., 2019).
Recently, bioinformatic methods have been widely applied to identify the robust differentially expressed genes (DEGs) and signaling pathways in a variety of diseases (Zeng et al., 2019). Weighted gene co-expression network analysis (WGCNA) is a widely used method in building co-expression pairwise correlation matrices (Xu et al., 2020). There are several studies on EA and NA which focus on DEGs screening, while the interconnection between the involved genes and different subtypes of asthma has not been investigated (Baines et al., 2014; Sánchez-Ovando et al., 2021).
In this study, we integrated two microarray datasets including EA and NA patients from gene expression omnibus (GEO) database. Gene Ontology (GO), Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses, and Gene set enrichment analysis (GSEA) were used to identify the potential mechanisms that distinguish between EA and NA. In order to explore the relation between gene modules and asthma subtypes, WGCNA was performed and two modules were recognized to be positively related to EA and NA groups. Finally, we sorted out a series of genes based on gene significance (GS), module membership (MM) and the protein-protein interaction (PPI) network, which might play major roles in the pathogenesis in EA and NA.
Materials and Methods
Microarray Datasets Acquisition
We obtained microarray datasets from GEO (http://www.ncbi.nlm.nih.gov/geo) utilizing the getGEO function of the GEOquery package (version 2.58.0) in R software (version 4.0.4). We adopted (eosinophilic asthma) AND (neutrophilic asthma) as the search strategy and five microarray data were obtained. The following eligibility criteria were used to include or exclude datasets and samples: 1) Expression type is expression profiling by array; 2) Sample type is induced sputum; 3) Subjects are non-severe and steroid-naïve. GSE4511 (Baines et al., 2014) and GSE137368 were extracted from the GEO database. Both platforms were GPL6104. We obtained 15 EA samples and 11 NA samples in GSE45111 and 9 EA samples and 4 NA samples in GSE137268.
Data Processing and Differentially Expressed Genes Identification
We considered consolidating two datasets since both used the same platforms. Data processing included data consolidation, batch normalization and ID conversion. The merged data were processed via ComBat normalization in sva package (version 3.38.0) in R software based on the classical Bayesian analysis to remove batch effects (Leek et al., 2012). We next identified DEGs between EA and NA group by an empirical Bayes method based on limma package (version 3.46.0) in R software (p < 0.05, |logFC|>0.5) (Ritchie et al., 2015). Volcano plot and heatmap were plotted by ggplot2 (version 3.3.3) and pheatmap (version 1.0.12) package, respectively.
Functional and Pathway Enrichment Analyses
The GO and KEGG pathway enrichment analyses of DEGs were analyzed and visualized in clusterProfiler package (version 3.18.1) in R software (Yu et al., 2012). The Metascape website (http://metascape.org), which is an online analysis tool integrated with multiple ontology sources, was implemented to conduct GO and KEGG analyses of gene modules selected by WGCNA (Zhou et al., 2019).
Gene Set Enrichment Analysis
GSEA analysis was performed using GSEA software (version 4.1.0) (Subramanian et al., 2005). KEGG was selected as the gene sets database. The gene set was deemed to significantly enriched with alpha or p values < 5% and false discovery rate (FDR) < 25% for each analysis, which was performed 1,000 times for each analysis. The parameter of GSEA are as following: The parameter “Collapse data set to gene symbols” is set to “false”. The parameter “Permutation type” is set to “phenotype”. The parameter “Enrichment statistic” is set to “weighted”, while the parameter “Metric for ranking genes” is set to “Singal2Noise”.
Weighted Gene Co-Expression Network Analysis
We implemented WGCNA by the WGCNA package in the R software (Langfelder and Horvath, 2008). We selected top 5,000 median absolute deviation (MAD) genes to construct the representation matrix and the appropriate power parameter was decided by pickSoftThreshold function. We used hierarchical clustering to identify modules of highly interconnected genes on the basis of their connectivity and covariance coefficients. The heatmap was plotted to reflect the relationships between each module and subtypes of asthma. The modules were constructed with the threshold value of the module dendrogram of 0.25, and a minimum module size of 30 genes.
Hub Genes and Key Genes Validation
Hub gene is defined as having high connectivity in a module and playing an important role in related clinical traits (Bai et al., 2021). Gene significance (GS) is defined as mediated p-value of each gene (GS = lgP) in the linear regression between gene expression and the clinical traits. Module membership (MM) is defined as the correlation between a given gene expression profile and a given model eigengene to represent intra connectivity. To identify hub genes in WGCNA, the threshold value of MM and GS was set to 0.8 and 0.2, respectively.
Hub genes in DEGs were constructed as the PPI network with the cut-off standard of interaction score >0.4 by the STRING database (http://www.string-db.org/) (Szklarczyk et al., 2021). Then, we visualized the PPI network by the Cytoscape software (version 3.6.1) (Le and Pham, 2017). In addition, we used the maximal clique centrality (MCC) computing method with the Cytohubba plugin to confirm the top 30 hub genes in the PPI network (Chin et al., 2014). Key genes were obtained by the intersection of hub genes in DEG-PPI network and selected modules eventually.
Results
Identification of Differentially Expressed Genes
To identify DEGs in induced sputum samples from EA patients compared to those from NA patients, expression profiles of GSE45111 and GSE137268 were downloaded from the GEO database. After data consolidation and removing the batch effects of the two microarray datasets, 282 DEGs (110 up-regulated and 172 down-regulated in EA patients) were screened out using the Limma package (p < 0.05, |logFC| > 0.5) (Supplementary Table S1). The result was shown with the volcano plot in Figure 1A. In addition, the logarithmic fold changes of the top 50 DEGs were showed in heatmap (Figure 1B).
FIGURE 1. Identification of differentially expressed genes (DEGs) in induced sputum from eosinophilic asthma (EA) compared to those from neutrophilic asthma (NA). (A) Volcano plots showed the DEGs by the criteria of |logFC| > 0.5 and p < 0.05. The up-regulated genes were marked in red, while the down-regulated genes were marked in blue. (B) The top 50 DEGs with the largest logFC are shown in the heatmap.
Functional and Pathway Enrichment Analyses for Differentially Expressed Genes
The enriched GO and KEGG analyses of the 282 DEGs were performed and visualized with the clusterProfiler package in R software. The biological process (BP) of GO enrichment analysis included response to interferon-γ (INFγ), response to virus and cellular response to INFγ. For the cellular component (CC), DEGs were enriched in tertiary granule, external side of plasma membrane, and specific granule. For the molecular function (MF), DEGs were significantly enriched in cytokine activity, chemokine activity, and cytokine receptor binding (Figure 2A). The result of GO enrichment of ascending and descending DEGs were shown in Supplementary Figures S1A,B, respectively. In addition, integrated DEGs were strongly involved in cytokine-cytokine receptor interaction, viral protein interaction with cytokine and cytokine receptor, and NOD-like receptor signaling pathway in KEGG pathway analysis (Figure 2B). The outcomes of KEGG pathway enrichment of ascending and descending DEGs were shown in Supplementary Figures S1C,D, respectively.
FIGURE 2. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG)enrichment annotations of the DEGs. (A) The results of GO enrichment categories included biological process (BP), cellular component (CC), and molecular function (MF). (B) The results of KEGG pathway enrichment analyses of the DEGs.
Gene set enrichment analysi
The GSEA analysis was performed to identify unique pathways involved in the pathogenesis of EA or NA. The pathways related to ribosome, Parkinson’s disease, and oxidative phosphorylation were most significantly enriched in EA group. Toll-like receptor signaling pathway, primary immunodeficiency, and NOD-like receptor signaling pathways were most strongly enriched in NA group (Figures 3A–E).
FIGURE 3. Gene set enrichment analysis (GSEA) plots of the most enriched gene sets in the EA and NA group. (A–C) The top 3 most enriched pathways in the EA group: ribosome (A); Parkinson’s disease (B); oxidative phosphorylation (C). (D–F) The top 3 most enriched pathways in the NA group: toll-like receptor significant pathway (D); primary immunodeficiency (E); NOD-like receptor signaling pathway (F).
Weighted Gene Co-Expression Network Analysis
WGCNA was performed to get a deeper insight into the association between the key modules and different asthma subtypes. As was shown in Figures 4A,B, the optimal soft-thresholding power was 12 if the correlational coefficient was >0.85. Eleven modules were generated via the average-linkage hierarchical clustering method (Figure 4C). The heatmap exhibited the correlation between different modules and subtypes of asthma. The magenta module was strongly positively correlated with EA group, so was the pink module with NA group (Figure 4D).
FIGURE 4. Identification of the gene modules related to EA and NA by weighted gene co-expression network analysis (WGCNA). (A) T scale-free fit index of soft-thresholding power. (B) Mean connectivity of various soft-thresholding powers. (C) Hierarchical clustering dendrograms of identified co-expressed genes were classified into different gene modules. Diverse colors reflected corresponding modules, and the grey module represented genes were not assigned to each network. (D) The heatmap of the relationship between each gene module and each subtype of asthma. The red represents positive correlation, while the green represents negative correlation.
Enrichment Analyses of Module Genes Identified by Weighted Gene Co-Expression Network Analysis
To further analyze the feature of the module genes, GO annotation and KEGG pathway enrichment analyses were performed using Metascape. As was shown in the bar graph and network plot (Figures 5A,C), genes in the magenta module were mainly involved in translational initiation, ribonucleoprotein complex biogenesis, mitochondrial protein complex, cytoplasmic translational initiation, and regulation of translation. The significantly enriched entries for pink module were defense response to virus, regulation of innate response, response to INFγ, and so on (Figures 5B,D).
FIGURE 5. GO and KEGG enrichment analyses of module genes recognized by the WGCNA analysis. Bar plots of the GO and KEGG enriched terms colored by p-values in magenta module (A) and pink module (B). Network of GO and KEGG enriched terms colored by genes in magenta module (C) and pink module (D).
The Identification of Hub Genes and Key Genes
In the magenta and pink modules, there was positive correlation between MM and GS (Supplementary Figures S2A,B). Under the criterial of MM ≥ 0.8 and GS ≥ 0.2, 139 and 45 hub genes were recognized in magenta and pink genes, respectively (Supplementary Table S2).
To identify the hub genes in DEGs, the PPI network, which was constructed by STRING database, was visualized by the Cytoscape software. The network of DEGs included 129 nodes and 1,054 edges (Figure 6A). The top 30 hub genes in DEGs were confirmed by the MCC method with the Cytohubba plugin. To further filtrate the significant genes, we selected the overlapped genes that were identified by hub genes in DEG-PPI network and WGCNA as key genes (Figures 6B,C). The detailed information of 10 key genes (IRFG, IRF1, STAT1, IFIH1, IFIT3, GBP1, GBP5, IFIT2, CXCL9, and CXCL11) were presented in Table 1, including 10 genes in pink module.
FIGURE 6. Identification of the hub genes and the key genes. The PPI network of the DEGs. (A) Top 30 hub genes in DEGs were confirmed with the Cytohubba plugin. The colors of the nodes reflect the degree of connectivity. The key genes are defined as the hub genes identified by both of DEG-PPI network and WGCNA method. (B) The hub genes in DEGs and in magenta module were shown using Venn diagram. No key genes were identified. (C) The hub genes in DEGs and in pink module were shown using Venn diagram. Ten key genes were identified.
TABLE 1. The information of 10 key genes.
Discussion
Airway neutrophilia is associated with asthma severity and poor responsiveness to steroid treatment (Ray and Kolls, 2017). In the present study, 282 DEGs in induced sputum samples were identified between patients with EA and NA, including 110 up-regulated genes and 172 down-regulated genes. Next, enrichment analyses, including GO, KEGG, and GSEA, were performed to explore the functions and signaling pathways related to the DEGs. After the construction of PPI network, the top 30 hub genes of DEGs were selected according to the Cytohubba algorithm. Using WGCNA analysis, pink and magenta modules were found to be strongly positively correlated with the EA and NA subsets, respectively. Finally, we defined 10 key genes according to intersection of the hub genes in DEGs and modules.
We performed enrichment analyses to explore the role of the DEGs. The significantly enriched entries for GO annotation demonstrated that DEGs were enriched in response to INFγ and cellular response to INFγ. INFγ is a Th1 cytokine that inhibits or reverses the allergic inflammation and therefore antagonizes the activating effects of Th2 cytokines including IL-4 and IL-13 in a variety of cell types (Davoodi et al., 2012). EA is orchestrated by Th2 cytokines, whereas NA is triggered by Th1 and Th17 cytokines (Erle and Sheppard, 2014; Trejo Bittar et al., 2015). This may explain the enrichment of DEGs in INFγ signaling pathway.
The GSEA data suggested that ribosome, Parkinson’s disease, and oxidative phosphorylation were mainly enriched in EA patients. Several epidemiological studies have implied that ribosome-inactivating stress are related to human mucosal epithelial illnesses (Moon, 2014). Moreover, it has been reported that intranasal neutrophilic rhinitis can be triggered by some of ribosome-inactivating trichothecenes (Islam et al., 2007; Corps et al., 2010; Carey et al., 2012). Both NA and neutrophilic rhinitis are characterized by neutrophilic inflammation. Therefore, NA subset may develop ribosomal inactivation leading to ribosome pathway enriched in EA group. Cheng’s study suggested that patients with asthma had a higher risk of developing Parkinson’s disease in their later life (Bower et al., 2006; Cheng et al., 2015). The relation between the pathogenesis of Parkinson’s disease and asthma requires further study. Toll-like receptor signaling pathway, primary immunodeficiency, and NOD-like receptor signaling pathway were mainly enriched in NA subset. Toll-like receptors (TLRs) belong to pattern recognition receptors, which play an important role in the recognition of pathogens (Xiao et al., 2012). TLRs can affect epithelial and immune cell function in asthma (Mishra et al., 2018). Of note, TLR4 is essential for Th17-driven neutrophilic airway inflammation and neutrophil recruitment (McAlees et al., 2015; Wan et al., 2020). NOD-like receptors (NLRs) are a relatively new member of the pattern recognition receptor superfamily (Song and Li, 2018). They are the key players in the innate immune responses of inflammatory lung diseases (Chaput et al., 2013). Several studies have demonstrated that NLRP3 plays an important role in asthma (Birrell and Eltom, 2011; Simpson et al., 2014; Kim et al., 2017; Chen et al., 2019). Chen and colleagues reported that blockade of the NLRP3/caspase-1 axis prevented the progression of TDI-induced NA (Chen et al., 2019).
To further investigate the relationship between co-expression genes and different asthma subtypes, we performed WGCNA. Eleven modules were defined and magenta and pink modules were most significantly positively correlated with EA and NA group, respectively. The result of the enrichment for the magenta module included translational initiation, ribonucleoprotein complex biogenesis, mitochondrial protein complex, cytoplasmic translational initiation, and regulation of translation. The result of the enrichment for the magenta module included defense response to virus, regulation of innate response, response to IFNγ, and so on. Viral infection is a common trigger for the exacerbation of asthma (Mikhail and Grayson, 2019). Neutrophils are considered to play a pivotal role in the interplay between viral infection and asthma exacerbation (Holtzman, 2012). Our data support the connection among viral infections, asthma exacerbation and NA.
After taking intersection of the hub genes in DEGs and the modules, 10 keys genes were identified, including IRFG, IRF1, STAT1, IFIH1, IFIT3, GBP1, GBP5, IFIT2, CXCL9, and CXCL11. IFNγ (IRFG) is a member of type II interferon involved in the chemotaxis of human neutrophils by up-regulating the expression of neutrophils chemokine receptors CCR1 and CCR3 (Bonecchi et al., 1999). Recent studies have revealed that IFNγ could up-regulate and down-regulate the expression of CXCL10 and SLPI, respectively, which further resulted in increased AHR and steroid resistance in severe asthma (Raundhal et al., 2015; Oriss et al., 2017). Interferon regulatory factor 1 (IRF1) is involved in a series of pathophysiological processes in viral infection, tumor immune surveillance, and proinflammatory injury (Wang et al., 2020). Increased IRF1 expression is implicated in reduced responsiveness to glucocorticoids (Tliba et al., 2008; Chapin et al., 2015). Glucocorticoids induce DUSP1 expression and downregulates MAPK activity, thereby inhibiting inflammatory response (Newton et al., 2017). However, increased DUSP1 expression also increases the activity of IRF1 and IRF1-dependent genes, including CXCL10 (Shah et al., 2016). CXCL10 can promote airway inflammation and hyperresponsiveness, and virus-induced exacerbation of asthma, resulting in poor response to glucocorticoids (Medoff et al., 2002; Wark et al., 2007). In conclusion, IRF1 could be a key gene for glucocorticoids resistance of NA. Signal transducer and activator of transcription 1 (STAT1) regulates production of Th1 cell-specific cytokine to alter inflammatory response and it is a key mediator in IFNγ signaling (Chen et al., 2017; Zhang et al., 2018). Several studies have shown that STAT1 participants in the differentiation of Th17 cells and it may become a key gene in the development of NA (Villarino et al., 2010; Kimura et al., 2008; Nielsen et al., 2015). Interferon induced with helicase C domain 1 (IFIH1) encodes MDA5, which is an intracellular sensor of viral RNA that triggers an innate immune response and is associated with the production of type I interferon and proinflammatory cytokines (Zhang et al., 2011; Mibayashi et al., 2007; Pichlmair et al., 2009). Guanylate binding protein 1 (GBP1) and guanylate binding protein 5 (GBP5) are two member of the GBP family and mediate cellular response to IFNγ in infection and inflammation (Tripal et al., 2007; Britzen-Laurent et al., 2010; Li et al., 2020). Interferon induced protein with tetratricopeptide repeats 3 (IFIT3) and interferon induced protein with tetratricopeptide repeats 2 (IFIT2) are both members of IFITs and are highly expressed in the innate immune response of cells to viral infection (Fleith et al., 2018). C-X-C motif chemokine ligand 9 (CXCL9) and C-X-C motif chemokine ligand 11 (CXCL11), ligands of chemokine receptor CXCR3, are induced by IFNγ (Proost et al., 2004; Tworek et al., 2013). These are especially involved in Th1-type response and correlates with the tissue infiltration of T cells (Koper et al., 2018). Southworth found that higher expression of CXCL11 appeared in moderate asthma after rhinovirus infection (Southworth et al., 2020). Ghebre reported that sputum CXCL9 level and serum CXCL11 level increased during asthma exacerbation (Ghebre et al., 2019).
It is intriguing to notice that IFN-stimulated genes were highly enriched in the pink module, which indicated the significant role of interferon in the pathogenesis of NA. A recent study demonstrated that IFN-stimulated genes expression is increased in severe asthma (Bhakta et al., 2018). In addition, Silva’s study revealed the overexpression of ISGs in sputum from NA (da Silva et al., 2017). The interferon family consists of three kinds of interferons, namely type I interferon, type II interferon, and type III interferon (Negishi et al., 2018). Type II interferon (IFNγ) has been discussed above. Type I interferon (interferon-alpha/beta) and Type III interferon (interferon-lambda) play crucial roles in host defense against infectious agents (Hansel et al., 2013), thus inhibiting the exacerbation of asthma. Nevertheless, the role of type I interferon and type III interferon in NA has not been reported yet (da Silva et al., 2017; Rich et al., 2020).
There are several limitations in this work. First, we have not obtained valuable clinical data, especially induced sputum cell counts. Second, cut-off standard of logFC was not high enough so that some of the expression difference of key genes between EA and NA group was subtle. Third, this work had the low number of samples and lack of the validation.
In conclusion, we have identified several enriched pathways in the EA compared to the NA. By intersection of hub genes in DEGs and modules, 10 key genes were defined. These key genes may provide new insights into the pathogenesis of NA, and become potential therapeutic targets of NA.
Data Availability Statement
The datasets analyzed during the present study are available from the Gene Expression Omnibus repository (https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE45111 and https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE137268).
Author Contributions
GC and GZ designed the research. GC, YF, WW, and CC acquired data. GC, DC, and SC analyzed the data. GC, DC, and JG visualized the data. GC and GZ interpreted the results. GC and GZ wrote the manuscript with all authors providing feedback on the revision. All authors read and approved the manuscript.
Funding
The work was supported by National Key Research and Development Program of China (2016YFC1304400); National Natural Science Foundation of China (Grants 82170036, 81670019, and 91742108).
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.
Acknowledgments
We sincerely thank Dr. Gibson PG and colleagues for sharing their data in GEO database.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmolb.2022.805570/full#supplementary-material
Supplementary Figure S1 | Function annotations analyses of DEGs. GO enrichment terms of up- and down-regulated genes are presented in A, B. KEGG pathway terms of up- and down-regulated genes are presented in C, D.
Supplementary Figure S2 | A scatter plot of the GS for the EA or NA classification versus the MM in the magenta (A) or pink (B) modules; the correlation and P-value are under the title. EA, eosinophilic asthma; NA, neutrophilic asthma; GS, gene significance; MM, module membership.
References
Bai, Z., Xu, L., Dai, Y., Yuan, Q., and Zhou, Z. (2021). ECM2 and GLT8D2 in Human Pulmonary Artery Hypertension: Fruits from Weighted Gene Co-expression Network Analysis. J. Thorac. Dis. 13 (4), 2242–2254. doi:10.21037/jtd-20-3069
Baines, K. J., Simpson, J. L., Wood, L. G., Scott, R. J., Fibbens, N. L., Powell, H., et al. (2014). Sputum Gene Expression Signature of 6 Biomarkers Discriminates Asthma Inflammatory Phenotypes. J. Allergy Clin. Immunol. 133 (4), 997–1007. doi:10.1016/j.jaci.2013.12.1091
Bhakta, N. R., Christenson, S. A., Nerella, S., Solberg, O. D., Nguyen, C. P., Choy, D. F., et al. (2018). IFN-stimulated Gene Expression, Type 2 Inflammation, and Endoplasmic Reticulum Stress in Asthma. Am. J. Respir. Crit. Care Med. 197 (3), 313–324. doi:10.1164/rccm.201706-1070oc
Birrell, M. A., and Eltom, S. (2011). The Role of the NLRP3 Inflammasome in the Pathogenesis of Airway Disease. Pharmacol. Ther. 130 (3), 364–370. doi:10.1016/j.pharmthera.2011.03.007
Bonecchi, R., Polentarutti, N., Luini, W., Borsatti, A., Bernasconi, S., Locati, M., et al. (1999)., 162. Baltimore, Md, 474–479.Up-regulation of CCR1 and CCR3 and Induction of Chemotaxis to CC Chemokines by IFN-Gamma in Human NeutrophilsJ. Immunol.1
Boonpiyathad, T., Sözener, Z. C., Satitsuksanoa, P., and Akdis, C. A. (2019). Immunologic Mechanisms in Asthma. Semin. Immunol. 46, 101333. doi:10.1016/j.smim.2019.101333
Bower, J. H., Maraganore, D. M., Peterson, B. J., Ahlskog, J. E., and Rocca, W. A. (2006). Immunologic Diseases, Anti-inflammatory Drugs, and Parkinson Disease: a Case-Control Study. Neurology 67 (3), 494–496. doi:10.1212/01.wnl.0000227906.99570.cc
Britzen-Laurent, N., Bauer, M., Berton, V., Fischer, N., Syguda, A., Reipschläger, S., et al. (2010). Intracellular Trafficking of Guanylate-Binding Proteins Is Regulated by Heterodimerization in a Hierarchical Manner. PLoS One 5 (12), e14246. doi:10.1371/journal.pone.0014246
Carey, S. A., Plopper, C. G., Hyde, D. M., Islam, Z., Pestka, J. J., and Harkema, J. R. (2012). Satratoxin-G from the Black Mold Stachybotrys chartarum Induces Rhinitis and Apoptosis of Olfactory Sensory Neurons in the Nasal Airways of Rhesus Monkeys. Toxicol. Pathol. 40 (6), 887–898. doi:10.1177/0192623312444028
Chapin, W. J., Lenkala, D., Mai, Y., Mao, Y., White, S. R., and Huang, R. S. (2015). Peripheral Blood IRF1 Expression as a Marker for Glucocorticoid Sensitivity. Pharmacogenetics and genomics 25 (3), 126–133. doi:10.1097/fpc.0000000000000116
Chaput, C., Sander, L. E., Suttorp, N., and Opitz, B. (2013). NOD-like Receptors in Lung Diseases. Front. Immunol. 4, 393. doi:10.3389/fimmu.2013.00393
Chen, S., Yao, L., Huang, P., He, Q., Guan, H., Luo, Y., et al. (2019). Blockade of the NLRP3/Caspase-1 Axis Ameliorates Airway Neutrophilic Inflammation in a Toluene Diisocyanate-Induced Murine Asthma Model. Toxicol. Sci. 170 (2), 462–475. doi:10.1093/toxsci/kfz099
Chen, Z., Pan, J., Jia, Y., Li, D., Min, Z., Su, X., et al. (2017). Effect of Memory CD4+ T Cells' Signal Transducer and Activator of Transcription (STATs) Functional Shift on Cytokine-Releasing Properties in Asthma. Cell Biol Toxicol 33 (1), 27–39. doi:10.1007/s10565-016-9357-6
Cheng, C.-M., Wu, Y.-H., Tsai, S.-J., Bai, Y.-M., Hsu, J.-W., Huang, K.-L., et al. (2015). Risk of Developing Parkinson's Disease Among Patients with Asthma: a Nationwide Longitudinal Study. Allergy 70 (12), 1605–1612. doi:10.1111/all.12758
Chin, C. H., Chen, S. H., Wu, H. H., Ho, C. W., Ko, M. T., and Lin, C. Y. (2014). cytoHubba: Identifying Hub Objects and Sub-networks from Complex Interactome. BMC Syst. Biol. 8 Suppl 4 (Suppl. 4Suppl 4), S11. doi:10.1186/1752-0509-8-S4-S11
Corps, K. N., Islam, Z., Pestka, J. J., and Harkema, J. R. (2010). Neurotoxic, Inflammatory, and Mucosecretory Responses in the Nasal Airways of Mice Repeatedly Exposed to the Macrocyclic Trichothecene Mycotoxin Roridin A. Toxicol. Pathol. 38 (3), 429–451. doi:10.1177/0192623310364026
da Silva, J., Hilzendeger, C., Moermans, C., Schleich, F., Henket, M., Kebadze, T., et al. (2017). Raised Interferon-β, Type 3 Interferon and Interferon-Stimulated Genes - Evidence of Innate Immune Activation in Neutrophilic Asthma. Clin. Exp. Allergy 47 (3), 313–323. doi:10.1111/cea.12809
Davoodi, P., Mahesh, P. A., Holla, A. D., Vijayakumar, G. S., Jayaraj, B. S., Chandrashekara, S., et al. (2012). Serum Levels of Interleukin-13 and Interferon-Gamma from Adult Patients with Asthma in Mysore. Cytokine 60 (2), 431–437. doi:10.1016/j.cyto.2012.05.012
Erle, D. J., and Sheppard, D. (2014). The Cell Biology of Asthma. J. Cel Biol 205 (5), 621–631. doi:10.1083/jcb.201401050
Fleith, R. C., Mears, H. V., Leong, X. Y., Sanford, T. J., Emmott, E., Graham, S. C., et al. (2018). IFIT3 and IFIT2/3 Promote IFIT1-Mediated Translation Inhibition by Enhancing Binding to Non-self RNA. Nucleic Acids Res. 46 (10), 5269–5285. doi:10.1093/nar/gky191
Gao, P., Tang, K., Lu, Y., Huang, Z., Wang, S., Wang, M., et al. (2020). Pentraxin 3 Promotes Airway Inflammation in Experimental Asthma. Respir. Res. 21 (1), 237. doi:10.1186/s12931-020-01499-6
Ghebre, M. A., Pang, P. H., Desai, D., Hargadon, B., Newby, C., Woods, J., et al. (2019). Severe Exacerbations in Moderate-To-Severe Asthmatics Are Associated with Increased Pro-inflammatory and Type 1 Mediators in Sputum and Serum. BMC Pulm. Med. 19 (1), 144. doi:10.1186/s12890-019-0906-7
Green, R. H., Brightling, C. E., McKenna, S., Hargadon, B., Parker, D., Bradding, P., et al. (2002). Asthma Exacerbations and Sputum Eosinophil Counts: a Randomised Controlled Trial. The Lancet 360 (9347), 1715–1721. doi:10.1016/s0140-6736(02)11679-5
Hansel, T. T., Johnston, S. L., and Openshaw, P. J. (2013). Microbes and Mucosal Immune Responses in Asthma. The Lancet 381 (9869), 861–873. doi:10.1016/s0140-6736(12)62202-8
Holtzman, M. J. (2012). Asthma as a Chronic Disease of the Innate and Adaptive Immune Systems Responding to Viruses and Allergens. J. Clin. Invest. 122 (8), 2741–2748. doi:10.1172/jci60325
Ijaz, H. M., Chowdhury, W., Lodhi, M. U., Gulzar, Q., and Rahim, M. (2019). A Case of Persistent Asthma Resistant to Available Treatment Options: Management Dilemma. Cureus 11 (3), e4194. doi:10.7759/cureus.4194
Islam, Z., Amuzie, C. J., Harkema, J. R., and Pestka, J. J. (2007). Neurotoxicity and Inflammation in the Nasal Airways of Mice Exposed to the Macrocyclic Trichothecene Mycotoxin Roridin a: Kinetics and Potentiation by Bacterial Lipopolysaccharide Coexposure. Toxicol. Sci. 98 (2), 526–541. doi:10.1093/toxsci/kfm102
Jayaram, L., Pizzichini, M. M., Cook, R. J., Boulet, L. P., Lemière, C., Pizzichini, E., et al. (2006). Determining Asthma Treatment by Monitoring Sputum Cell Counts: Effect on Exacerbations. Eur. Respir. J. 27 (3), 483–494. doi:10.1183/09031936.06.00137704
Kim, R. Y., Pinkerton, J. W., Essilfie, A. T., Robertson, A. A. B., Baines, K. J., Brown, A. C., et al. (2017). Role for NLRP3 Inflammasome-Mediated, IL-1β-Dependent Responses in Severe, Steroid-Resistant Asthma. Am. J. Respir. Crit. Care Med. 196 (3), 283–297. doi:10.1164/rccm.201609-1830oc
Kimura, A., Naka, T., Nohara, K., Fujii-Kuriyama, Y., and Kishimoto, T. (2008). Aryl Hydrocarbon Receptor Regulates Stat1 Activation and Participates in the Development of Th17 Cells. Proc. Natl. Acad. Sci. 105 (28), 9721–9726. doi:10.1073/pnas.0804231105
Koper, O., Kamińska, J., Sawicki, K., and Kemona, H. (2018). CXCL9, CXCL10, CXCL11, and Their Receptor (CXCR3) in Neuroinflammation and neurodegenerationAdvances in Clinical and Experimental Medicine. Adv. Clin. Exp. Med. 27 (6), 849–856. doi:10.17219/acem/68846
Langfelder, P., and Horvath, S. (2008). WGCNA: an R Package for Weighted Correlation Network Analysis. BMC Bioinformatics 9, 559. doi:10.1186/1471-2105-9-559
Le, D.-H., and Pham, V.-H. (2017). HGPEC: a Cytoscape App for Prediction of Novel Disease-Gene and Disease-Disease Associations and Evidence Collection Based on a Random Walk on Heterogeneous Network. BMC Syst. Biol. 11 (1), 61. doi:10.1186/s12918-017-0437-x
Leek, J. T., Johnson, W. E., Parker, H. S., Jaffe, A. E., and Storey, J. D. (2012). The Sva Package for Removing Batch Effects and Other Unwanted Variation in High-Throughput Experiments. Bioinformatics 28 (6), 882–883. doi:10.1093/bioinformatics/bts034
Li, Z., Qu, X., Liu, X., Huan, C., Wang, H., Zhao, Z., et al. (2020). GBP5 Is an Interferon-Induced Inhibitor of Respiratory Syncytial Virus. J. Virol. 94 (21). doi:10.1128/JVI.01407-20
McAlees, J. W., Whitehead, G. S., Harley, I. T., Cappelletti, M., Rewerts, C. L., Holdcroft, A. M., et al. (2015). Distinct Tlr4-Expressing Cell Compartments Control Neutrophilic and Eosinophilic Airway Inflammation. Mucosal Immunol. 8 (4), 863–873. doi:10.1038/mi.2014.117
Medoff, B. D., Sauty, A., Tager, A. M., Maclean, J. A., Smith, R. N., Mathew, A., et al. (2002). IFN-gamma-inducible Protein 10 (CXCL10) Contributes to Airway Hyperreactivity and Airway Inflammation in a Mouse Model of Asthma. J. Immunol. 168 (10), 5278–5286. doi:10.4049/jimmunol.168.10.5278
Mibayashi, M., Martínez-Sobrido, L., Loo, Y.-M., Cárdenas, W. B., Gale, M., and García-Sastre, A. (2007). Inhibition of Retinoic Acid-Inducible Gene I-Mediated Induction of Beta Interferon by the NS1 Protein of Influenza A Virus. J. Virol. 81 (2), 514–524. doi:10.1128/jvi.01265-06
Mikhail, I., and Grayson, M. H. (2019). Asthma and Viral Infections. Ann. Allergy Asthma Immunol. 123 (4), 352–358. doi:10.1016/j.anai.2019.06.020
Mishra, V., Banga, J., and Silveyra, P. (2018). Oxidative Stress and Cellular Pathways of Asthma and Inflammation: Therapeutic Strategies and Pharmacological Targets. Pharmacol. Ther. 181, 169–182. doi:10.1016/j.pharmthera.2017.08.011
Moon, Y. (2014). Ribosomal Alteration-Derived Signals for Cytokine Induction in Mucosal and Systemic Inflammation: Noncanonical Pathways by Ribosomal Inactivation. Mediators Inflamm. 2014, 708193. doi:10.1155/2014/708193
Negishi, H., Taniguchi, T., and Yanai, H. (2018). The Interferon (IFN) Class of Cytokines and the IFN Regulatory Factor (IRF) Transcription Factor Family. Cold Spring Harb Perspect. Biol. 10 (11). doi:10.1101/cshperspect.a028423
Newton, R., Shah, S., Altonsy, M. O., and Gerber, A. N. (2017). Glucocorticoid and Cytokine Crosstalk: Feedback, Feedforward, and Co-regulatory Interactions Determine Repression or Resistance. J. Biol. Chem. 292 (17), 7163–7172. doi:10.1074/jbc.r117.777318
Nielsen, J., Kofod-Olsen, E., Spaun, E., Larsen, C. S., Christiansen, M., and Mogensen, T. H. (2015). A STAT1-Gain-Of-Function Mutation Causing Th17 Deficiency with Chronic Mucocutaneous Candidiasis, Psoriasiform Hyperkeratosis and Dermatophytosis. BMJ Case Rep. 2015. doi:10.1136/bcr-2015-211372
Oriss, T. B., Raundhal, M., Morse, C., Huff, R. E., Das, S., Hannum, R., et al. (2017). IRF5 Distinguishes Severe Asthma in Humans and Drives Th1 Phenotype and Airway Hyperreactivity in Mice. JCI insight 2 (10). doi:10.1172/jci.insight.91019
Pelaia, G., Vatrella, A., Busceti, M. T., Gallelli, L., Calabrese, C., Terracciano, R., et al. (2015). Cellular Mechanisms Underlying Eosinophilic and Neutrophilic Airway Inflammation in Asthma. Mediators Inflamm. 2015, 879783. doi:10.1155/2015/879783
Pichlmair, A., Schulz, O., Tan, C.-P., Rehwinkel, J., Kato, H., Takeuchi, O., et al. (2009). Activation of MDA5 Requires Higher-Order RNA Structures Generated during Virus Infection. J. Virol. 83 (20), 10761–10769. doi:10.1128/jvi.00770-09
Proost, P., Verpoest, S., De Borne, K. V., Schutyser, E., Struyf, S., Put, W., et al. (2004). Synergistic Induction of CXCL9 and CXCL11 by Toll-like Receptor Ligands and Interferon-γ in Fibroblasts Correlates with Elevated Levels of CXCR3 Ligands in Septic Arthritis Synovial Fluids. J. Leukoc. Biol. 75 (5), 777–784. doi:10.1189/jlb.1003524
Raundhal, M., Morse, C., Khare, A., Oriss, T. B., Milosevic, J., Trudeau, J., et al. (2015). High IFN-γ and Low SLPI Mark Severe Asthma in Mice and Humans. J. Clin. Invest. 125 (8), 3037–3050. doi:10.1172/jci80911
Ray, A., and Kolls, J. K. (2017). Neutrophilic Inflammation in Asthma and Association with Disease Severity. Trends Immunology 38 (12), 942–954. doi:10.1016/j.it.2017.07.003
Rich, H. E., Antos, D., Melton, N. R., Alcorn, J. F., and Manni, M. L. (2020). Insights into Type I and III Interferons in Asthma and Exacerbations. Front. Immunol. 11, 574027. doi:10.3389/fimmu.2020.574027
Ritchie, M. E., Phipson, B., Wu, D., Hu, Y., Law, C. W., Shi, W., et al. (2015). Limma powers Differential Expression Analyses for RNA-Sequencing and Microarray Studies. Nucleic Acids Res. 43 (7), e47. doi:10.1093/nar/gkv007
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.
Scherzer, R., and Grayson, M. H. (2018). Heterogeneity and the Origins of Asthma. Ann. Allergy Asthma Immunol. 121 (4), 400–405. doi:10.1016/j.anai.2018.06.009
Schleich, F., Demarche, S., and Louis, R. (2016). Biomarkers in the Management of Difficult Asthma. Ctmc 16 (14), 1561–1573. doi:10.2174/1568026616666151015093406
Shah, S., King, E. M., Mostafa, M. M., Altonsy, M. O., and Newton, R. (2016). DUSP1 Maintains IRF1 and Leads to Increased Expression of IRF1-dependent Genes. J. Biol. Chem. 291 (41), 21802–21816. doi:10.1074/jbc.m116.728964
Simpson, J. L., Phipps, S., Baines, K. J., Oreo, K. M., Gunawardhana, L., and Gibson, P. G. (2014). Elevated Expression of the NLRP3 Inflammasome in Neutrophilic Asthma. Eur. Respir. J. 43 (4), 1067–1076. doi:10.1183/09031936.00105013
Song, N., and Li, T. (2018). Regulation of NLRP3 Inflammasome by Phosphorylation. Front. Immunol. 9, 2305. doi:10.3389/fimmu.2018.02305
Southworth, T., Pattwell, C., Khan, N., Mowbray, S. F., Strieter, R. M., Erpenbeck, V. J., et al. (2020). Increased Type 2 Inflammation post Rhinovirus Infection in Patients with Moderate Asthma. Cytokine 125, 154857. doi:10.1016/j.cyto.2019.154857
Subramanian, A., Tamayo, P., Mootha, V. K., Mukherjee, S., Ebert, B. L., Gillette, M. A., et al. (2005). Gene Set Enrichment Analysis: a Knowledge-Based Approach for Interpreting Genome-wide Expression Profiles. Proc. Natl. Acad. Sci. 102 (43), 15545–15550. doi:10.1073/pnas.0506580102
Szklarczyk, D., Gable, A. L., Nastou, K. C., Lyon, D., Kirsch, R., Pyysalo, S., et al. (2021). The STRING Database in 2021: Customizable Protein-Protein Networks, and Functional Characterization of User-Uploaded Gene/measurement Sets. Nucleic Acids Res. 49 (D1), D605–D612. doi:10.1093/nar/gkaa1074
Taylor, S. L., Leong, L. E. X., Choo, J. M., Wesselingh, S., Yang, I. A., Upham, J. W., et al. (2018). Inflammatory Phenotypes in Patients with Severe Asthma Are Associated with Distinct Airway Microbiology. J. Allergy Clin. Immunol. 141 (1), 94–103. doi:10.1016/j.jaci.2017.03.044
Tliba, O., Damera, G., Banerjee, A., Gu, S., Baidouri, H., Keslacy, S., et al. (2008). Cytokines Induce an Early Steroid Resistance in Airway Smooth Muscle Cells. Am. J. Respir. Cel Mol Biol 38 (4), 463–472. doi:10.1165/rcmb.2007-0226oc
Tliba, O., and Panettieri, R. A. (2019). Paucigranulocytic Asthma: Uncoupling of Airway Obstruction from Inflammation. J. Allergy Clin. Immunol. 143 (4), 1287–1294. doi:10.1016/j.jaci.2018.06.008
Trejo Bittar, H. E., Yousem, S. A., and Wenzel, S. E. (2015). Pathobiology of Severe Asthma. Annu. Rev. Pathol. Mech. Dis. 10, 511–545. doi:10.1146/annurev-pathol-012414-040343
Tripal, P., Bauer, M., Naschberger, E., Mörtinger, T., Hohenadl, C., Cornali, E., et al. (2007). Unique Features of Different Members of the Human Guanylate-Binding Protein Family. J. Interferon Cytokine Res. 27 (1), 44–52. doi:10.1089/jir.2007.0086
Tworek, D., Kuna, P., Młynarski, W., Górski, P., Pietras, T., and Antczak, A. (2013). MIG (CXCL9), IP-10 (CXCL10) and I-TAC (CXCL11) Concentrations after Nasal Allergen challenge in Patients with Allergic Rhinitis. aoms 5 (5), 849–853. doi:10.5114/aoms.2013.37198
Villarino, A. V., Gallo, E., and Abbas, A. K. (2010). STAT1-activating Cytokines Limit Th17 Responses through Both T-bet-dependent and -independent Mechanisms. J. Immunol., 185, 6461–6471. doi:10.4049/jimmunol.1001343
Wan, R., Jiang, J., Hu, C., Chen, X., Chen, C., Zhao, B., et al. (2020). Neutrophil Extracellular Traps Amplify Neutrophil Recruitment and Inflammation in Neutrophilic Asthma by Stimulating the Airway Epithelial Cells to Activate the TLR4/NF-Κb Pathway and Secrete Chemokines. Aging 12 (17), 16820–16836. doi:10.18632/aging.103479
Wang, J., Li, H., Xue, B., Deng, R., Huang, X., Xu, Y., et al. (2020). IRF1 Promotes the Innate Immune Response to Viral Infection by Enhancing the Activation of IRF3. J. Virol. 94 (22). doi:10.1128/JVI.01231-20
Wark, P. A. B., Bucchieri, F., Johnston, S. L., Gibson, P. G., Hamilton, L., Mimica, J., et al. (2007). IFN-γ-induced Protein 10 Is a Novel Biomarker of Rhinovirus-Induced Asthma Exacerbations. J. Allergy Clin. Immunol. 120 (3), 586–593. doi:10.1016/j.jaci.2007.04.046
Xiao, H. T., Liao, Z., Chen, L., and Tong, R. S. (2012). A Promising Approach for Asthma Treatment by Multiwayly Modulating Toll-like Receptors. Eur. Rev. Med. Pharmacol. Sci. 16 (15), 2088–2091.
Xu, M., Ouyang, T., Lv, K., and Ma, X. (2020). Integrated WGCNA and PPI Network to Screen Hub Genes Signatures for Infantile Hemangioma. Front. Genet. 11, 614195. doi:10.3389/fgene.2020.614195
Yu, G., Wang, L.-G., Han, Y., and He, Q.-Y. (2012). clusterProfiler: an R Package for Comparing Biological Themes Among Gene Clusters. OMICS: A J. Integr. Biol. 16 (5), 284–287. doi:10.1089/omi.2011.0118
Zeng, M., Liu, J., Yang, W., Zhang, S., Liu, F., Dong, Z., et al. (2019). Multiple-microarray Analysis for Identification of Hub Genes Involved in Tubulointerstial Injury in Diabetic Nephropathy. J. Cel Physiol. doi:10.1002/jcp.28313
Zhang, Y., Chen, Y., Liu, Z., and Lai, R. (2018). ERK Is a Negative Feedback Regulator for IFN-Γ/stat1 Signaling by Promoting STAT1 Ubiquitination. BMC Cancer 18 (1), 613. doi:10.1186/s12885-018-4539-7
Zhang, Z., Yuan, B., Lu, N., Facchinetti, V., and Liu, Y. J. (2011). DHX9 Pairs with IPS-1 to Sense Double-Stranded RNA in Myeloid Dendritic Cells. J. Immunol., 187, 4501–4508. doi:10.4049/jimmunol.1101307
Keywords: eosinophilic asthma, neutrophilic asthma, key genes, bioinformatics, weight gene co-expression network analysis
Citation: Chen G, Chen D, Feng Y, Wu W, Gao J, Chang C, Chen S and Zhen G (2022) Identification of Key Signaling Pathways and Genes in Eosinophilic Asthma and Neutrophilic Asthma by Weighted Gene Co-Expression Network Analysis. Front. Mol. Biosci. 9:805570. doi: 10.3389/fmolb.2022.805570
Received: 30 October 2021; Accepted: 10 January 2022;
Published: 02 February 2022.
Edited by:
Grzegorz Wegrzyn, University of Gdansk, PolandReviewed by:
Bahman Panahi, Agricultural Biotechnology Research Institute of Iran, IranMagdalena Julita Orzechowska, Medical University of Lodz, Poland
Dorota Anusewicz, Medical University of Lodz, Poland
Copyright © 2022 Chen, Chen, Feng, Wu, Gao, Chang, Chen and Zhen. 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: Guohua Zhen, ghzhen@tjh.tjmu.edu.cn