- 1Yongchuan Hospital of Chongqing Medical University, Chongqing, China
- 2JiangJin Central Hosptial of Chongqing, Chongqing, China
- 3Department of Pharmacy, First People’s Hospital of Chongqing Liangjiang New District, Chongqing, China
Despite the recent progress of lung adenocarcinoma (LUAD) therapy, tumor recurrence remained to be a challenging factor that impedes the effectiveness of treatment. The objective of the present study was to predict the hub genes affecting LUAD recurrence via weighted gene co-expression network analysis (WGCNA). Microarray samples from LUAD dataset of GSE32863 were analyzed, and the modules with the highest correlation to tumor recurrence were selected. Functional enrichment analysis was conducted, followed by establishment of a protein–protein interaction (PPI) network. Subsequently, hub genes were identified by overall survival analyses and further validated by evaluation of expression in both myeloid populations and tissue samples of LUAD. Gene set enrichment analysis (GSEA) was then carried out, and construction of transcription factors (TF)–hub gene and drug–hub gene interaction network was also achieved. A total of eight hub genes (ACTR3, ARPC5, RAB13, HNRNPK, PA2G4, WDR12, SRSF1, and NOP58) were finally identified to be closely correlated with LUAD recurrence. In addition, TFs that regulate hub genes have been predicted, including MYC, PML, and YY1. Finally, drugs including arsenic trioxide, cisplatin, Jinfukang, and sunitinib were mined for the treatment of the eight hub genes. In conclusion, our study may facilitate the invention of targeted therapeutic drugs and shed light on the understanding of the mechanism for LUAD recurrence.
Introduction
As one of the most frequently diagnosed and severe tumors worldwide, lung cancer was estimated to cause almost one-quarter of all cancer deaths in 2021 (Siegel et al., 2021). Lung adenocarcinoma (LUAD) is the most common type of lung cancer, which accounts for more than 40% of lung cancer incidence (Shi et al., 2016). Although the resection of early-stage LUAD remains to be the best treatment option, over 30% of LUAD patients develop recurrence and finally lead to dismal outcomes (Carnio et al., 2013; Uramoto and Tanaka, 2014). Therefore, identifying core regulators governing LUAD recurrence and elucidating the underlying mechanisms involved in the progression of LUAD would promote the development of effective strategies for diagnosis, treatment, and prognosis of this devastating disease.
As a systematic algorithm to identify tumor-specific indicators and predict cancer-related signaling pathways, the weighted gene co-expression network analysis (WGCNA) approach provides a comprehensive strategy to elucidate the interactions of pathogenic genes and to determine the correlation between gene networks and clinical traits (Langfelder and Horvath, 2008). With the help of WGCNA, genes were classified into different modules for functional prediction, while the most central genes could be further targeted as hub genes. Other than previous bioinformatics methodologies that mainly focus on individual genes, WGCNA implements methods for both weighted and unweighted correlation networks and displays the characteristics of biological systems more precisely (Barabási et al., 2011).
With the rapid advancement in recent years, WGCNA has been widely applied in the research of various cancer types (Liang et al., 2020; Yao et al., 2020; Liu et al., 2021). In terms of LUAD, DYNLRB2 and SPTBN1 were predicted to be tumor suppressors and may serve as biomarkers for LUAD patients (Zhu et al., 2021), while KIF11 was identified to be essential for LUAD cell proliferation and metastasis (Li et al., 2021). By utilizing WGCNA, long noncoding RNA (lncRNA) SVIl-AS1 was revealed to associate with chemoresistance in LUAD by acting as competing endogenous RNA (ceRNA) (Guo et al., 2021). Moreover, high-throughput data analysis by WGCNA indicates that SPP1 may be a key regulator in LUAD, which could be directly regulated by four miRNAs and indirectly regulated by 49 lncRNAs (Luo et al., 2021). However, comprehensive identification of biomarkers correlated with LUAD recurrence using WGCNA is yet to be investigated.
As an emerging technology, single-cell RNA-sequencing (scRNA-seq) is widely accepted as a powerful approach to study the heterogeneity of gene expression in individual cells (Suvà and Tirosh, 2019). By sequencing of LUAD cells harboring EGFR mutations, elevated expression of ELF3 was observed in advanced tumor cells, which promotes tumorigenesis through PI3K/AKT/NF-κB signaling pathway (He et al., 2021). Analysis of scRNA-seq data from LUAD samples manifested as subsolid nodules suggested that malignant cells in subsolid nodules undergo a strong metabolic reprogram and immune stress (Xing et al., 2021). Moreover, a multi-region scRNA-seq study revealed that LUAD exhibits pronounced intratumor cell heterogeneity within single sites and transcriptional lineage-plasticity programs (Sinjab et al., 2021).
Here in this study, WGCNA was constructed by using the GSE32863 dataset containing gene expression profiling data from 58 LUAD tumors and their matched normal adjacent lung tissue samples. Specific modules correlated with LUAD recurrence were identified, followed by hub gene prediction and functional analyses. Noteworthy, single-cell analysis was further applied to determine the expression of hub genes in different LUAD cell types, and the potential interactions between hub genes and therapeutic drugs were also investigated.
Materials and Methods
Data Acquisition and Preprocessing
The microarray dataset GSE32863 of human LUAD with patient clinical information was extracted from Gene Expression Omnibus (GEO) database (Selamat et al., 2012). GSE32863 dataset was obtained using Illumina HumanWG-6 v3.0 expression beadchip. It includes 58 LUAD samples, which were used to construct a co-expression network followed by extraction of hub genes. R package was used to annotate raw data, generate an expression matrix, and match probes with official gene symbols. The median absolute deviations (MADs) were ranked from the largest to smallest, and the expressions of the top 25% genes with the largest differences in the samples were extracted for in-depth analysis.
Weighted Gene Co-Expression Network Analysis
We used the R package “WGCNA” to perform WGCNA on selected genes to find out the expression patterns between different genes (Langfelder and Horvath, 2008). Genes with similar expression patterns were regarded as a specific module and marked with a unique color. Subsequently, the correlation between these modules was calculated, and a heatmap was created to show the independence between each module. Then, a correlation analysis was performed to find modules related to clinical information. Modules closely related to tumor recurrence were chosen for further analysis.
Gene Ontology and Pathway Enrichment Analysis
In order to understand the biological functions and signal pathways involved in the genes from the selected modules, differentially expressed genes (DEGs) in these modules were analyzed through the R package “clusterProfiler” (Yu et al., 2012). Then, another R package “ggplot2” was used to perform the top 10 enrichment terms of Gene Ontology (GO) analysis and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis. Enriched GO terms and KEGG pathways were identified base on the cutoff criterion of p < 0.05.
Protein–Protein Interaction of the Key Module Genes
For hub gene screening, genes from the selected modules were uploaded to the STRING database to build a protein–protein interaction (PPI) network (Szklarczyk et al., 2019). The interaction score >0.4 was defined as the threshold of the key genes in the PPI network. Then, a Cytoscape plug-in cytoHubba was used to extract the top 1% targets in these modules based on the degree method (Shannon et al., 2003).
Identification and Validation of Real Hub Genes
Based on the cytoHubba analysis, the top 15 genes in the selected modules and the genes enriched in KEGG pathway were selected as candidate hub targets for further analysis. In this study, the Kaplan–Meier plotter (http://kmplot.com/analysis/) was used to draw the survival analyses of the candidate hub targets, while the Gene Expression Profiling Interactive Analysis (GEPIA) webserver (Tang et al., 2017) (http://gepia.cancer-pku.cn/) was used to verify the outcomes of survival analysis so as to screen out real hub genes. With the use of the Human Protein Atlas (http://www.proteinatlas.org) database, the real hub genes were further validated by immunohistochemistry (IHC). Moreover, Single Cell Portal (https://singlecell.broadinstitute.org/single_cell) was used to detect the expression of key targets in various cell types of lung cancer. In addition, the cBioPortal tool (http://www.cbioportal.org/; version: 2.2.0) was used to compare the genetic variations of the real hub genes in LUAD. Finally, the cBioPortal tool was also used to present the co-expression analysis of real hub genes.
Gene Set Enrichment Analysis
In the validation set GSE116959, samples of LUAD were divided into low and high groups according to the expression level of the genes. To explore the potential function of real hub genes in LUAD, gene set enrichment analysis (GSEA) was conducted and mapped into KEGG pathway enrichment database. The terms with p-value <0.05 were chosen as the cutoff criteria.
Construction of Transcription Factor–Real Hub Gene Network and Drug–Real Hub Gene Interaction
A regulating network on hub genes and LUAD was constructed by Cytoscape software. Then, the plugin iRegulon of Cytoscape was applied to forecast transcription factor (TF) regulation networks. In addition, Comparative Toxicogenomics Database (http://ctdbase.org/), a robust and publicly available database that aims to reveal how environmental exposures affect human health, was used to search for drugs for the hub genes. Chemicals supported by at least one database were selected as the potential drugs. Potential drug–key target interaction was constructed by Cytoscape. The final list only involves drugs that may interact with at least two key targets and have been reported to have an anticancer pharmacological activity.
Results
Identification of Key Modules Associated With Lung Adenocarcinoma Recurrence
The flowchart of strategy in this study is summarized in Figure 1. GSE32863, a GEO dataset containing gene expression profiling data of 58 matched LUAD and non-tumor lung samples, was analyzed using the R package WGCNA. Clinical characteristics of tumor pathological stage, gender, recurrence, egfr, kras, lkb1, and smoking status of LUAD patients were denoted. After screening by MADs arranged from large to small, the expression of the top 25% genes (6,360 genes) with the greatest differences in samples were selected for further analyses.
After classification of DEGs with similar expression patterns into modules by average linkage clustering, 14 modules were finally identified by merging similar modules when the MedissThres was set at 0.25 (Figure 2A). Among them, module eigengenes of three modules (green, purple, and brown) were found to be more related to the tumor recurrence as others. Scatter plots further indicated positive correlations between members of these three modules and gene significance for LUAD recurrence (Figures 2B–D). Herein, the green, purple, and brown modules were chosen as key modules associated with LUAD recurrence for further investigation.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/756235/fgene-12-756235-HTML/image_m/fgene-12-756235-g002.jpg)
FIGURE 2. Identification of modules associated with the clinical features of lung adenocarcinoma (LUAD). (A) Heatmap to show the correlation between modules and clinical traits with LUAD. p-Values are shown in brackets. (B–D) Scatter plot analysis to show the association between Module membership and gene significance for LUAD recurrence in green (B), purple (C), and brown modules (D).
Functional Analysis of Genes in Three Key Modules
To understand the biological functions of genes in the selected modules, GO and KEGG pathway analyses were firstly applied. As depicted in Figure 3, genes in three key modules were expected to exert their functions in terms of cell adhesion molecule binding, endosome membrane, ribosomal subunit and autophagy, etc. Meanwhile, these genes could only be enriched in one single cellular pathway called Tight Junction signaling pathway (Figure 4 and Table 1).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/756235/fgene-12-756235-HTML/image_m/fgene-12-756235-g003.jpg)
FIGURE 3. Gene Ontology (GO) enrichment analyses of gene members in three chosen modules. The colored geometry represents GO term enrichment, dots represent biological process (BP), triangles represent molecular function (CC), rectangles represent cellular component (MF), red indicates low enrichment, and green indicates high enrichment. The sizes of the geometries represent the number of genes in each GO category.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/756235/fgene-12-756235-HTML/image_m/fgene-12-756235-g004.jpg)
FIGURE 4. Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses. Genes in three modules (marked in red) could only be enriched in Tight Junction signaling pathway.
Next, genes in three key modules were extracted to establish a PPI network through STRING database and Cytoscape software. As a result, a total of 926 DEGs of the three key modules were mapped into the PPI network, including 1,426 nodes and 9,522 edges. In addition, the top 15 key targets within PPI network were selected with the help of cytoHubba plug-in in Cytoscape software based on rank of degree (Figure 5).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/756235/fgene-12-756235-HTML/image_m/fgene-12-756235-g005.jpg)
FIGURE 5. Identification of the top 15 genes from the protein–protein interaction (PPI) network of the three modules.
Identification and Validation of Hub Genes for Lung Adenocarcinoma Recurrence
According to the results above, 40 genes (25 genes could be mapped into the KEGG pathway plus 15 key targets in PPI network) remained the potentials to be the hub genes for LUAD recurrence. For the purpose of scope reduction and effect evaluation, overall survival analysis was firstly performed on these candidate genes using the survival information of Kaplan–Meier plotter. After the first round of screening, only genes with significant impact on LUAD patient survival were further validated by survival analysis with GEPIA. As a result, a total of eight candidate hub genes (ACTR3, ARPC5, RAB13, HNRNPK, PA2G4, WDR12, SRSF1, and NOP58) were identified to be significantly correlated with LUAD patient survival with both survival analysis methods (Figure 6). Remarkably, all eight candidate genes were found to be correlated with poor LUAD patient survival through GEPIA, five of which (ACTR3, ARPC5, RAB13, HNRNPK, and SRSF1) were inversely revealed to be associated with favorable outcome using Kaplan–Meier plotter.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/756235/fgene-12-756235-HTML/image_m/fgene-12-756235-g006.jpg)
FIGURE 6. Survival analyses of the eight hub gene candidates identified. (A–H) Overall survival of the eight hub gene candidates in lung adenocarcinoma (LUAD) based on Kaplan–Meier plotter. (I–P) Overall survival of the eight hub gene candidates in LUAD based on the Gene Expression Profiling Interactive Analysis (GEPIA) database. p < 0.05 was considered to indicate a statistically significant difference. ACTR3, actin-related protein 3; ARPC5, actin-related protein 2/3 complex subunit 5; HNRNPK, heterogeneous nuclear ribonucleoprotein K; NOP58, NOP58 ribonucleoprotein; PA2G4, proliferation-associated 2G4; RAB13, RAB13 member RAS oncogene family; SRSF1, serine- and arginine-rich splicing factor 1; WDR12, WD repeat domain 12.
We next sought to investigate the expression of these candidate genes in LUAD. Firstly, comprehensive single-cell RNA sequencing data mapping myeloid populations in non-small cell lung tumor and peripheral blood was utilized to detect the expression of our candidate genes in different types of immune cells (Zilionis et al., 2019). As presented in Figure 7, candidate genes ACTR3, ARPC5, and HNRNPK were highly expressed in almost all types of lung cancer immune cells, while RAB13, PA2G4, and WDR12 were found to be less abundant in the majority of myeloid populations in lung cancer. As compared with the rest of the hub genes, ACTR3, ARPC5, and HNRNPK are all characterized in regulation of actin polymerization (Yoo et al., 2006; Kopitar et al., 2019). Thus, the hyperactivation of these three genes in lung cancer cells may reflect the significance of actin polymerization during the development of cancer. Next, the distribution and expression of these eight hub gene candidates were checked using IHC staining data of both LUAD and adjacent non-tumor samples from The Human Protein Atlas database. Representative IHC images revealed elevated protein levels of HNRNPK, PA2G4, WDR12, SRSF1, and NOP58 in LUAD samples (Figure 8), suggesting the dysregulation of these genes during progression of LUAD.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/756235/fgene-12-756235-HTML/image_m/fgene-12-756235-g007.jpg)
FIGURE 7. Expression of eight hub gene candidates in various clusters of human immune cells from non-small cell lung tumor (t) and peripheral blood (b). t, tumor; b, blood; DCs, dendritic cells; pDCs, plasmacytoid DCs; RBC, red blood cell. The size of dots represents the percentage of expression; red and blue represent the level of scaled mean expression.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/756235/fgene-12-756235-HTML/image_m/fgene-12-756235-g008.jpg)
FIGURE 8. Immunohistochemical images of eight hub genes in lung adenocarcinoma (LUAD) recurrence. The protein levels of (A–E) in LUAD tissues were compared with normal lung tissues from the Human Protein Atlas database.
The genetic alterations of these candidate genes were then determined using cBioPortal. Interestingly, seven out of the eight genes were found to be amplificated in LUAD, while only half of them were reported to be partially mutated (Figures 9A,B). Gene co-expression analysis of the eight candidate genes was also performed using cBioPortal database. As shown in Figure 9C, highly positive co-expression was observed in pairs of PA2G4–NOP58, HNRNPK–SRSF1, and NOP58–SRSF1, while significantly negative co-expression was revealed in pairs of HNRNPK–ARPC5 and SRSF1–ARPC5.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/756235/fgene-12-756235-HTML/image_m/fgene-12-756235-g009.jpg)
FIGURE 9. Genetic alteration information and co-expression analysis of the eight hub genes. (A) A visual summary across a set of lung adenocarcinoma (LUAD) (data from Lung Adenocarcinoma, The Cancer Genome Atlas (TCGA), Nature 2014) showed the genetic alterations connected with the eight hub genes, which were altered in 75 (31.8%) of 230 sequenced patients (230 in total). (B) An overview of changes in the eight hub genes from the genomics datasets of LUAD in TCGA database. Summary for lung adenocarcinoma: Gene altered in 20% of 230 cases, Mutation 1.74% (4 in 230 cases), Amplification 17.39% (40 in 230 cases), and Multiple Alterations 0.87% (2 in 230 cases). (C) The co-expression analysis of the eight hub genes using the 230 samples above based on cBioPortal database.
To investigate the biological characteristics of right candidate genes associated with LUAD recurrence, the GSEA assay was further applied. The results in Figure 10 showed that upregulation of these candidate genes was enriched in multiple cellular processed such as Fc gamma R-mediated phagocytosis, regulation of actin cytoskeleton, spliceosome, and tight junction.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/756235/fgene-12-756235-HTML/image_m/fgene-12-756235-g010.jpg)
FIGURE 10. Gene set enrichment analysis (GSEA) using GSE116959 showed a positive enrichment of identified hub genes associated with four Kyoto Encyclopedia of Genes and Genomes(KEGG) pathways.
In total, all the above findings indicated that all eight candidate genes may serve as the hub genes associated with tumor recurrence in LUAD.
Constructions of Transcription Factor–Hub Gene Network and Drug–Hub Gene Interaction Network Associated With Lung Adenocarcinoma Recurrence
To enhance the significance of our study, we next sought to establish the transcriptional regulatory network of hub genes and TFs by a Cytoscape plug-in iRegulon. As revealed in Figure 11, a total of 57 TFs and eight hub genes were involved in this network, such as models of MYC–NOP58, MYC–WDR12, MYC–SRSF1, MYC–PA2G4, MYC–HNRNPK, PML–HNRNPK, TAF1–HNRNPK, YY1–HNRNPK, and TAF1–HNRNPK. To facilitate the future targeted drug screening, Comparative Toxicogenomics Database was used to search for drugs specialized in targeting the eight hub genes. Drugs that may interact with at least two hub genes and have been reported to have an anticancer pharmacological activity were selected. Finally, 11 drugs including arsenic trioxide, cisplatin, copper, ICG 001, Jinfukang, tretinoin, doxorubicin, sodium selenite, quercetin, sunitinib, and epigallocatechin gallate were discovered (Figure 12).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/756235/fgene-12-756235-HTML/image_m/fgene-12-756235-g011.jpg)
FIGURE 11. The transcriptional regulatory network of eight hub genes and TFs. TFs, transcription factors. A green hexagon node represents the TFs, a pink circular node represents hub genes, a light blue diamond node represents the lung adenocarcinoma (LUAD), and the interaction is represented by an arrow. The numbers of arrows in the networks demonstrate the contribution of 1 TF to the hub genes; and the higher the degree, the more central the nodes were within the network.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/756235/fgene-12-756235-HTML/image_m/fgene-12-756235-g012.jpg)
FIGURE 12. Interactions between drugs and the identified hub genes based on CTD database and literature validation. Light blue circular node represents drugs, red diamond node represents hub genes, their interaction is represented by line, and the degree value represents the number of targets acted on by the drug.
Discussion
Despite the novel improvements in the diagnosis and treatment of LUAD in recent years, tumor recurrence remained to be the leading cause of treatment failure and mortality after surgery (Martini et al., 1995). With the help of WGCNA approach, a total of eight hub genes were screened out as LUAD recurrence-associated hub genes in the present study. To the best of our knowledge, most of them were firstly reported as key modulators of LUAD recurrence.
Among the identified hub genes, both ACTR3 and ARPC5 encode subunits of the human Arp2/3 protein complex, which has been implicated in the control of actin polymerization in cells (Abella et al., 2016). As compared with the less studied ACTR3, the potential of ARPC5 as a prognostic biomarker has been investigated in multiple cancer types such as hepatocellular carcinoma (Huang et al., 2021) and multiple myeloma (Xiong and Luo, 2018). In this study, we further confirmed the association between both Arp2/3 components and LUAD recurrence, which may extend the clinical relevance of these biomarkers into the prognosis of lung cancer.
As a member of Rab-associated G protein family, RAB13 controls the trafficking and cellular localization of a number of key modulators in cancer and thus affects the tumorigenesis (Ioannou and McPherson, 2016). Based on The Cancer Genome Atlas (TCGA) database, elevated expression of RAB13 has been observed in the majority of cancers and is inversely correlated with patient prognosis. Moreover, by analyzing published microarray data for NCI-60 cancer cells, RAB13 was further discovered as a target gene relevant to radiosensitivity (Kim et al., 2012). Here, we observed significant amplification of RAB13 in LUAD sample, which deserves to be further investigated to unveil its role in the carcinogenesis of LUAD.
Interestingly, five out of the eight hub genes identified in the current study (HNRNPK, PA2G4, WDR12, SRSF1, and NOP58) are mainly located in the nucleus and are involved in functions of rRNA processing, RNA splicing, and ribosome biogenesis. For instance, HNRNPK is a component of the hnRNP complex and a highly conserved RNA- and DNA-binding protein (Pinol-Roma et al., 1988). Downregulation of HNRNPK in human LUAD cell line significantly reduced the formation of metastatic lung tumor nodules in mice, suggesting its role in lung metastasis (Li et al., 2019). Another example is NOP58, which forms box C/D small nucleolar ribonucleoprotein (snoRNP) with three other protein members to exert its functions in rRNA methylation and ribosome biogenesis (Shubina et al., 2016). Dysregulation of snoRNP members including NOP58 has been widely reported in cancer, which directly contributes to the overactivated ribosome biogenesis during cancer progression (Wu et al., 2020; Yi et al., 2021). Therefore, it is speculated that hyperactivation of ribosome biogenesis, which can be initiated by these hub genes, may play a crucial role in LUAD recurrence.
In addition, transcription factors (TFs) that regulate hub genes have been predicted, including MYC (Massó-Vallés et al., 2020), PML (Kuo et al., 2014), and YY1 (Lin et al., 2020). It has been reported that they were potential therapeutic targets for lung cancer and were closely related to lung cancer metastasis and recurrence. Finally, through CTD database and literature mining, it was discovered that arsenic trioxide (Huang and Zeng, 2019), Jinfukang (Que et al., 2021), sunitinib (Park and Kim, 2020), and other drugs reported to treat lung cancer could act on these eight hub genes.
In general, our current study aimed to find hub genes that might be correlated with the recurrence of LUAD using WGCNA approach. A total of eight genes including ACTR3, ARPC5, RAB13, HNRNPK, PA2G4, WDR12, SRSF1, and NOP58 were identified and further validated using multiple bioinformatics tools. It is anticipated that these hub genes could serve as biomarkers or therapeutic targets for LUAD treatment.
Data Availability Statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/Supplementary Material.
Author Contributions
CY conceived the idea. ZS and SL finished the bioinformatics analysis and drafted the manuscript. The other authors were involved in the data analysis and interpretation.
Funding
This work was supported by the Technology Innovation and Application Development Project of Chongqing Kewei United Chinese Medicine (2021ZY3591) and the Research Project of Yongchuan Hospital of Chongqing Medical University (YJJC202009).
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 thank the GEO, STRING, GEPIA, cBioPortal, TCGA, Kaplan–Meier plotter, the Human Protein Atlas, Single Cell Portal, and Comparative Toxicogenomics Database for providing their platforms and contributors for their valuable datasets.
References
Abella, J. V. G., Galloni, C., Pernier, J., Barry, D. J., Kjær, S., Carlier, M.-F., et al. (2016). Isoform Diversity in the Arp2/3 Complex Determines Actin Filament Dynamics. Nat. Cel Biol 18 (1), 76–86. doi:10.1038/ncb3286
Barabási, A.-L., Gulbahce, N., and Loscalzo, J. (2011). Network Medicine: a Network-Based Approach to Human Disease. Nat. Rev. Genet. 12 (1), 56–68. doi:10.1038/nrg2918
Carnio, S., Novello, S., Papotti, M., Loiacono, M., and Scagliotti, G. V. (2013). Prognostic and Predictive Biomarkers in Early Stage Non-small Cell Lung Cancer: Tumor Based Approaches Including Gene Signatures. Transl Lung Cancer Res. 2 (5), 372–381. doi:10.3978/j.issn.2218-6751.2013.10.05
Guo, L., Ding, L., and Tang, J. (2021). Identification of a Competing Endogenous RNA axis "SVIL‐AS1/miR‐103a/ICE1" Associated with Chemoresistance in Lung Adenocarcinoma by Comprehensive Bioinformatics Analysis. Cancer Med. 10, 6022–6034. doi:10.1002/cam4.4132
He, D., Wang, D., Lu, P., Yang, N., Xue, Z., Zhu, X., et al. (2021). Single-cell RNA Sequencing Reveals Heterogeneous Tumor and Immune Cell Populations in Early-Stage Lung Adenocarcinomas Harboring EGFR Mutations. Oncogene 40 (2), 355–368. doi:10.1038/s41388-020-01528-0
Huang, S., Li, D., Zhuang, L., Sun, L., and Wu, J. (2021). Identification of Arp2/3 Complex Subunits as Prognostic Biomarkers for Hepatocellular Carcinoma. Front. Mol. Biosci. 8, 690151. doi:10.3389/fmolb.2021.690151
Huang, W., and Zeng, Y. C. (2019). A Candidate for Lung Cancer Treatment: Arsenic Trioxide. Clin. Transl Oncol. 21 (9), 1115–1126. doi:10.1007/s12094-019-02054-6
Ioannou, M. S., and McPherson, P. S. (2016). Regulation of Cancer Cell Behavior by the Small GTPase Rab13. J. Biol. Chem. 291 (19), 9929–9937. doi:10.1074/jbc.R116.715193
Kim, H. S., Kim, S. C., Kim, S. J., Park, C. H., Jeung, H.-C., Kim, Y. B., et al. (2012). Identification of a Radiosensitivity Signature Using Integrative Metaanalysis of Published Microarray Data for NCI-60 Cancer Cells. BMC Genomics 13, 348. doi:10.1186/1471-2164-13-348
Kopitar, A. N., Markelj, G., Oražem, M., Blazina, Š., Avčin, T., Ihan, A., et al. (2019). Flow Cytometric Determination of Actin Polymerization in Peripheral Blood Leukocytes Effectively Discriminate Patients with Homozygous Mutation in ARPC1B from Asymptomatic Carriers and Normal Controls. Front. Immunol. 10, 1632. doi:10.3389/fimmu.2019.01632
Kuo, H.-Y., Huang, Y.-S., Tseng, C.-H., Chen, Y.-C., Chang, Y.-W., Shih, H.-M., et al. (2014). PML Represses Lung Cancer Metastasis by Suppressing the Nuclear EGFR-Mediated Transcriptional Activation ofMMP2. Cell Cycle 13 (19), 3132–3142. doi:10.4161/15384101.2014.949212
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
Li, M., Zhang, W., Yang, X., Liu, H., Cao, L., Li, W., et al. (2019). Downregulation of HNRNPK in Human Cancer Cells Inhibits Lung Metastasis. Anim. Models Exp. Med. 2 (4), 291–296. doi:10.1002/ame2.12090
Li, Z., Yu, B., Qi, F., and Li, F. (2021). KIF11 Serves as an Independent Prognostic Factor and Therapeutic Target for Patients with Lung Adenocarcinoma. Front. Oncol. 11, 670218. doi:10.3389/fonc.2021.670218
Liang, Q., Guan, G., Li, X., Wei, C., Wu, J., Cheng, P., et al. (2020). Profiling Pro-neural to Mesenchymal Transition Identifies a lncRNA Signature in Glioma. J. Transl Med. 18 (1), 378. doi:10.1186/s12967-020-02552-0
Lin, C.-C., Kuo, I.-Y., Wu, L.-T., Kuan, W.-H., Liao, S.-Y., Jen, J., et al. (2020). Dysregulated Kras/YY1/ZNF322A/Shh Transcriptional axis Enhances Neo-Angiogenesis to Promote Lung Cancer Progression. Theranostics 10 (22), 10001–10015. doi:10.7150/thno.47491
Liu, S., Zeng, F., Fan, G., and Dong, Q. (2021). Identification of Hub Genes and Construction of a Transcriptional Regulatory Network Associated with Tumor Recurrence in Colorectal Cancer by Weighted Gene Co-expression Network Analysis. Front. Genet. 12, 649752. doi:10.3389/fgene.2021.649752
Luo, X., Feng, L., Xu, W., Bai, X., and Wu, M. (2021). Weighted Gene Co-expression Network Analysis of Hub Genes in Lung Adenocarcinoma. Evol. Bioinform Online 17, 117693432110098. doi:10.1177/11769343211009898
Martini, N., Bains, M. S., Burt, M. E., Zakowski, M. F., McCormack, P., Rusch, V. W., et al. (1995). Incidence of Local Recurrence and Second Primary Tumors in Resected Stage I Lung Cancer. J. Thorac. Cardiovasc. Surg. 109 (1), 120–129. doi:10.1016/S0022-5223(95)70427-2
Massó-Vallés, D., Beaulieu, M.-E., and Soucek, L. (2020). MYC, MYCL, and MYCN as Therapeutic Targets in Lung Cancer. Expert Opin. Ther. Targets 24 (2), 101–114. doi:10.1080/14728222.2020.1723548
Park, K.-Y., and Kim, J. (2020). Cyclic Pentapeptide cRGDfK Enhances the Inhibitory Effect of Sunitinib on TGF-Β1-Induced Epithelial-To-Mesenchymal Transition in Human Non-small Cell Lung Cancer Cells. PLoS One 15 (8), e0232917. doi:10.1371/journal.pone.0232917
Pinol-Roma, S., Choi, Y. D., Matunis, M. J., and Dreyfuss, G. (1988). Immunopurification of Heterogeneous Nuclear Ribonucleoprotein Particles Reveals an Assortment of RNA-Binding Proteins. Genes Dev. 2 (2), 215–227. doi:10.1101/gad.2.2.215
Que, Z.-J., Yang, Y., Liu, H.-T., Shang-Guan, W.-J., Yu, P., Zhu, L.-H., et al. (2021). Jinfukang Regulates integrin/Src Pathway and Anoikis Mediating Circulating Lung Cancer Cells Migration. J. Ethnopharmacology 267, 113473. doi:10.1016/j.jep.2020.113473
Selamat, S. A., Chung, B. S., Girard, L., Zhang, W., Zhang, Y., Campan, M., et al. (2012). Genome-scale Analysis of DNA Methylation in Lung Adenocarcinoma and Integration with mRNA Expression. Genome Res. 22 (7), 1197–1211. doi:10.1101/gr.132662.111
Shannon, P., Markiel, A., Ozier, O., Baliga, N. S., Wang, J. T., Ramage, D., et al. (2003). Cytoscape: a Software Environment for Integrated Models of Biomolecular Interaction Networks. Genome Res. 13 (11), 2498–2504. doi:10.1101/gr.1239303
Shi, J., Hua, X., Zhu, B., Ravichandran, S., Wang, M., Nguyen, C., et al. (2016). Somatic Genomics and Clinical Features of Lung Adenocarcinoma: A Retrospective Study. Plos Med. 13 (12), e1002162. doi:10.1371/journal.pmed.1002162
Shubina, M. Y., Musinova, Y. R., and Sheval, E. V. (2016). Nucleolar Methyltransferase Fibrillarin: Evolution of Structure and Functions. Biochem. Mosc. 81 (9), 941–950. doi:10.1134/S0006297916090030
Siegel, R. L., Miller, K. D., Fuchs, H. E., and Jemal, A. (2021). Cancer Statistics, 2021. CA A. Cancer J. Clin. 71 (1), 7–33. doi:10.3322/caac.21654
Sinjab, A., Han, G., Treekitkarnmongkol, W., Hara, K., Brennan, P. M., Dang, M., et al. (2021). Resolving the Spatial and Cellular Architecture of Lung Adenocarcinoma by Multiregion Single-Cell Sequencing. Cancer Discov. 11, 2506–2523. doi:10.1158/2159-8290.Cd-20-1285
Suvà, M. L., and Tirosh, I. (2019). Single-Cell RNA Sequencing in Cancer: Lessons Learned and Emerging Challenges. Mol. Cel 75 (1), 7–12. doi:10.1016/j.molcel.2019.05.003
Szklarczyk, D., Gable, A. L., Lyon, D., Junge, A., Wyder, S., Huerta-Cepas, J., et al. (2019). STRING V11: Protein-Protein Association Networks with Increased Coverage, Supporting Functional Discovery in Genome-wide Experimental Datasets. Nucleic Acids Res. 47 (D1), D607–D613. doi:10.1093/nar/gky1131
Tang, Z., Li, C., Kang, B., Gao, G., Li, C., and Zhang, Z. (2017). GEPIA: a Web Server for Cancer and normal Gene Expression Profiling and Interactive Analyses. Nucleic Acids Res. 45 (W1), W98–W102. doi:10.1093/nar/gkx247
Uramoto, H., and Tanaka, F. (2014). Recurrence after Surgery in Patients with NSCLC. Transl Lung Cancer Res. 3 (4), 242–249. doi:10.3978/j.issn.2218-6751.2013.12.05
Wu, H., Qin, W., Lu, S., Wang, X., Zhang, J., Sun, T., et al. (2020). Long Noncoding RNA ZFAS1 Promoting Small Nucleolar RNA-Mediated 2′-O-Methylation via NOP58 Recruitment in Colorectal Cancer. Mol. Cancer 19 (1), 95. doi:10.1186/s12943-020-01201-w
Xing, X., Yang, F., Huang, Q., Guo, H., Li, J., Qiu, M., et al. (2021). Decoding the Multicellular Ecosystem of Lung Adenocarcinoma Manifested as Pulmonary Subsolid Nodules by Single-Cell RNA Sequencing. Sci. Adv. 7 (5). doi:10.1126/sciadv.abd9738
Xiong, T., and Luo, Z. (2018). The Expression of Actin-Related Protein 2/3 Complex Subunit 5 (ARPC5) Expression in Multiple Myeloma and its Prognostic Significance. Med. Sci. Monit. 24, 6340–6348. doi:10.12659/MSM.908944
Yao, Z., Liu, N., Zhu, X., Wang, L., Zhao, Y., Liu, Q., et al. (2020). Subanesthetic Isoflurane Abates ROS-Activated MAPK/NF-κB Signaling to Repress Ischemia-Induced Microglia Inflammation and Brain Injury. Aging 12 (24), 26121–26139. doi:10.18632/aging.202349
Yi, Y., Li, Y., Meng, Q., Li, Q., Li, F., Lu, B., et al. (2021). A PRC2-independent Function for EZH2 in Regulating rRNA 2′-O Methylation and IRES-dependent Translation. Nat. Cel Biol 23 (4), 341–354. doi:10.1038/s41556-021-00653-6
Yoo, Y., Wu, X., Egile, C., Li, R., and Guan, J.-L. (2006). Interaction of N-WASP with hnRNPK and its Role in Filopodia Formation and Cell Spreading. J. Biol. Chem. 281 (22), 15352–15360. doi:10.1074/jbc.M511825200
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
Zhu, H., Yue, H., Xie, Y., Chen, B., Zhou, Y., and Liu, W. (2021). Bioinformatics and Integrated Analyses of Prognosis-Associated Key Genes In Lung Adenocarcinoma. J. Thorac. Dis. 13, 1172–1186.
Keywords: lung adenocarcinoma, tumor recurrence, weighted gene co-expression network analysis, hub genes, transcription factor
Citation: Shen Z, Liu S, Liu J, Liu J and Yao C (2021) Weighted Gene Co-Expression Network Analysis and Treatment Strategies of Tumor Recurrence-Associated Hub Genes in Lung Adenocarcinoma. Front. Genet. 12:756235. doi: 10.3389/fgene.2021.756235
Received: 10 August 2021; Accepted: 06 October 2021;
Published: 16 November 2021.
Edited by:
Min Tang, Jiangsu University, ChinaReviewed by:
Yunpeng Xu, Rutgers, the State University of New Jersey, United StatesFeng Wang, Emory University, United States
Copyright © 2021 Shen, Liu, Liu, Liu and Yao. 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: Caoyuan Yao, ycy8210@163.com
†These authors have contributed equally to this work