ORIGINAL RESEARCH article

Front. Pharmacol., 02 June 2022

Sec. Pharmacology of Infectious Diseases

Volume 13 - 2022 | https://doi.org/10.3389/fphar.2022.857730

Bioinformatics and System Biology Approach to Reveal the Interaction Network and the Therapeutic Implications for Non-Small Cell Lung Cancer Patients With COVID-19

  • 1. Guangzhou University of Chinese Medicine, Guangzhou, China

  • 2. Department of Oncology, The First Affiliated Hospital of Guangzhou University of Chinese Medicine, Guangzhou, China

Article metrics

View details

7

Citations

5,5k

Views

1,8k

Downloads

Abstract

Background: Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), the leading cause of coronavirus disease-2019 (COVID-19), is an emerging global health crisis. Lung cancer patients are at a higher risk of COVID-19 infection. With the increasing number of non-small-cell lung cancer (NSCLC) patients with COVID-19, there is an urgent need of efficacious drugs for the treatment of COVID-19/NSCLC.

Methods: Based on a comprehensive bioinformatic and systemic biological analysis, this study investigated COVID-19/NSCLC interactional hub genes, detected common pathways and molecular biomarkers, and predicted potential agents for COVID-19 and NSCLC.

Results: A total of 122 COVID-19/NSCLC interactional genes and 21 interactional hub genes were identified. The enrichment analysis indicated that COVID-19 and NSCLC shared common signaling pathways, including cell cycle, viral carcinogenesis, and p53 signaling pathway. In total, 10 important transcription factors (TFs) and 44 microRNAs (miRNAs) participated in regulations of 21 interactional hub genes. In addition, 23 potential candidates were predicted for the treatment of COVID-19 and NSCLC.

Conclusion: This study increased our understanding of pathophysiology and screened potential drugs for COVID-19 and NSCLC.

Introduction

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), responsible for pandemic coronavirus disease-2019 (COVID-19), is a novel beta-coronavirus belonging to the subgenus Sarbecovirus (Ciotti et al., 2020; Umakanthan et al., 2020). As an emerging global health crisis, SARS-CoV-2 shares a similar transmission mode with other respiratory viruses, mainly through air droplets and close contact (Shang et al., 2021). Due to the rapidly evolving nature of SARS-CoV-2, as on 21 November 2021, there had been 257,788,585 confirmed cases with mortality calculated at 2% (Worldometer, 2021). Recently, the development of Delta and Omicron variants has further complicated the control of the pandemic (Lopez Bernal et al., 2021; Wang and Powell, 2021). About 80% of COVID-19 patients exhibit mild to moderate clinical manifestations, including fever, dyspnea, dry cough, and acute pneumonia (Liu et al., 2020; Mallah et al., 2021; Shang et al., 2021). The case fatality rate for COVID-19 shows a close connection between age, underlying disease status, and immune state (Olloquequi, 2020). A growing body of evidence shows that cancer patients harbor a higher risk of COVID-19 infection, along with severe events and unfavorable outcomes (Passaro et al., 2021; Sinha and Kundu, 2021). A cohort study of 1,590 COVID-19 cases reported that lung cancer was the most frequent type of cancer [5 (28%) of 18 cancer patients] due to the inherent associated pulmonary fragility (Liang et al., 2020). Moreover, lung cancer patients with smoking-related lung damage, significant cardiovascular or respiratory comorbidities, and older age are more likely to develop COVID-19 severity (Berlin et al., 2020; Liang et al., 2020; Passaro et al., 2020). According to the Global Cancer Statistics in 2020, lung cancer remains to be the leading cause of cancer incidence and mortality, representing 11.4% and 18.0% of all cases, respectively (Sung et al., 2021). Accounting for approximately 85% of lung cancers, non-small-cell lung cancer (NSCLC) is comprised of several histological subtypes such as lung adenocarcinoma, squamous-cell carcinoma, and large-cell carcinoma (Zappa and Mousa, 2016). Despite chemotherapy and targeted therapies being widely applied for the treatment of NSCLC, the 5-year survival rate has remained abysmally low (16%) for the last four decades (Suresh et al., 2019).

Cell proliferation is a vital and fundamental mechanism for growth, development, and regeneration of eukaryotic organisms (Diaz-Moralli et al., 2013). Dysregulation of the cell cycle leads to aberrant cell proliferation, which is found in various malignancies (Williams and Stoeber, 2012). Most NSCLCs have detectable cell cycle abnormalities, and the more defective the cell cycle becomes, the more severe the consequences would be (Sterlacci et al., 2012). Driven by complex interactions between host factors, tumorigenesis creates an ideal tumor microenvironment and promotes tumor formation (Dzobo, 2021). Prior studies have pointed out that the tumor microenvironment not only played an important part in tumor development at primary and metastatic sites but also deteriorated viral infection (Tian et al., 2020; Malkani and Rashid, 2021). In lung cancer patients, the tumor microenvironment supports SARS-CoV-2 proteins by activating cytokine storm- and cellular metabolic variation-related pathways, which further accelerate infection and weaken the immune system. Angiotensin-converting enzyme 2 (ACE2), an entry receptor for SARS-CoV-2, almost ubiquitously present in human organs, but primarily expressed in alveolar epithelial type II cells, secretes surfactant and plays a crucial part in pulmonary gas exchange (Zhou P. et al., 2020; Walls et al., 2020). In lung cancer tissues, significantly upregulated ACE2 caused lung parenchyma to become vulnerable to SARS-CoV-2 attack (Zhang L. et al., 2020; Mason, 2020). By triggering associated cascades, SARS-CoV-2 infections increased inflammatory mediators, which could induce paracrine senescence through prolonging cytokine signaling. Accumulative evidence suggests that cellular senescence damages vascular functions persistently by impairing endothelium and decreasing angiogenesis. Vascular dysfunction strikes the balance between anti- and procoagulant pathways and further increases the risk of abnormal coagulation (thrombosis) (Bochenek et al., 2016; Escher et al., 2020). To conclude, persistent inflammation and cellular senescence are associated with pulmonary parenchyma injury, which could potentially restrict blood flow to the lungs, activate coagulation, induce capillary damage, and eventually contribute to hypoxemia and acute respiratory distress in COVID-19 and NSCLC patients.

In the context of the highly contagious COVID-19 pandemic, the number of NSCLC patients with COVID-19 continues to rise (Worldometer, 2021). Although remdesivir (Eastman et al., 2020; Hung et al., 2020), ribavirin, and hydroxychloroquine (Yao et al., 2020) have been approved for COVID-19, none of them have been suggested to be specific. With the widespread application of high-throughput technologies, a large amount of biological data has been generated. We performed a comprehensive bioinformatic and systemic biological analysis to further understand mechanisms and seek potential efficacious drugs to combat COVID-19/NSCLC (Figure 1).

FIGURE 1

FIGURE 1

Overview of the workflow of the present study.

Materials and Methods

Identification of Genetic Interrelations Between COVID-19 and Non-Small-Cell Lung Cancer

To determine genetic interrelations shared by COVID-19 and NSCLC, we searched a variety of databases to collect available data for each disease. We searched six online databases for COVID-19–related genes, including Online Mendelian Inheritance in Man (OMIM, https://www.omim.org/) (Hamosh et al., 2005), Therapeutic Target Database (TTD, http://db.idrblab.net/ttd/) (Wang et al., 2020), PubChem (https://pubchem.ncbi.nlm.nih.gov/) (Kim S. et al., 2019), DisGeNET (https://www.disgenet.org/covid/diseases/summary/) (Piñero et al., 2020), GeneCards (https://www.genecards.org/) (Stelzer et al., 2016), and Comparative Toxicogenomics Database (CTD, http://ctdbase.org/) (Davis et al., 2021). Transcriptomic RNA-sequencing (RNA-seq) datasets of COVID-19 [GSE147507 (Daamen et al., 2021), GSE157103 (Overmyer et al., 2021), and GSE166190 (Vono et al., 2021)] were downloaded from Gene Expression Omnibus (GEO, https://www.ncbi.nlm.nih.gov/geo/) (Clough and Barrett, 2016). Then, we used the R package Deseq2 (Love et al., 2014) to extract differentially expressed genes (DEGs) at the threshold |log2FoldChange| >1 and adjusted the p-value <0.05. COVID-19–related genes were identified by intersecting related genes from online databases and GEO datasets.

In addition, five online platforms were searched for NSCLC-related genes, including OMIM (Hamosh et al., 2005), CTD (Davis et al., 2021), TTD (Wang et al., 2020), DisGeNET (Piñero et al., 2020), and GeneCards (Stelzer et al., 2016). RNA-seq datasets of NSCLC were obtained from The Cancer Genome Atlas (TCGA) data portal on the UCSC Xena database (https://xenabrowser.net/datapages/) (Goldman et al., 2020). Data of 1,135 tissues (1,027 cancer tissues and 108 para-cancerous tissues) from 1,016 NSCLC patients were obtained.

To obtain the robust and biologically significant DEG list, we extracted DEGs from TCGA-NSCLC–related genes by combining data from lung cancers with different pathological classifications, according to the previous studies (Han et al., 2019; Zhang J. et al., 2020; Zhang et al., 2021). We employed the R package Deseq2 (Love et al., 2014), set the threshold |log2FoldChange| at >1, and adjusted the p-value at <0.05. NSCLC-related genes were collected by intersecting related genes from public platforms and TCGA data portal. After intersecting the COVID-19–related genes and NSCLC-related genes, we identified genes shared by COVID-19 and NSCLC as interactional genes. These genes might play important roles in pathophysiological processes of COVID-19 and NSCLC and serve as important clues for screening candidate drugs for NSCLC patients with COVID-19.

Protein–Protein Interaction Network Analysis and Sub-Network Analysis

To interpret associated cellular machinery operations and explore protein mechanisms, PPI network analysis was conducted using STRING (https://string-db.org/) based on proteins derived from COVID-19/NSCLC interactional genes. Active interaction sources of the PPI mainly included text-mining, experiments, databases, co-expression, neighborhood, gene fusion, and co-occurrence. Moreover, to assure the highest confidence of the network, the minimum confidence score was set at 0.90 for network construction. Disconnected nodes were removed from the network. Then Cytoscape 3.9.0 software (Shannon et al., 2003) was applied for network visualization. In the network, the significance of each node was evaluated by betweenness centrality, degree and closeness centrality. In addition, a network analysis module named cytoHubba (Chin et al., 2014) was employed for detecting interactional hub genes. Analytical methods including betweenness, bottleneck, closeness, degree, density of maximum neighborhood component (DMNC), eccentricity, edge-percoalated component (EPC), maximal clique centrality (MCC), maximum neighborhood component (MNC), radiality, and stress were applied to obtain top 10 genes of each method. Based on previous studies, hub genes were defined as genes with a degree value twice the median or more in the whole network (Guo et al., 2015; Yu et al., 2018; Zhuang et al., 2020). Hub genes and genes obtained by the cytoHubba module were intersected, and duplication was removed. By employing these 12 approaches, COVID-19/NSCLC interactional hub genes were identified. Furthermore, sub-network analyses were conducted by the Molecular Complex Detection (MCODE) (Bader and Hogue, 2003) module of Cytoscape 3.9.0 software (Shannon et al., 2003) to identify important gene clusters in the whole PPI network. MCODE is a useful tool for detecting densely connected regions in large protein–protein interaction networks, which helps to identify gene clusters and understand the connectivity and proximity of genes. In the present study, the parameters of MCODE were set as follows: degree cutoff: 2; cluster finding method: haircut; node score cutoff: 0.2; K-Core: 2; and maximal depth: 100.

Gene Ontology and Kyoto Encyclopedia of Genes and Genomes Pathway Analysis

Based on COVID-19/NSCLC interactional genes, GO analysis and KEGG pathway analyses were carried out by the R package clusterProfiler (Wu et al., 2021) to explore potential pathogenic mechanisms of COVID-19 and NSCLC. GO analysis classified associated mechanisms into three categories, including biological processes (BP), cellular components (CC), and molecular functions (MF), and then KEGG pathway analysis further specified related mechanisms. p-value <0.05 was set as a standard metric to quantify the most closely related GO and KEGG terms. Furthermore, to explore the function of interactional genes in the biggest sub-network, GO terms and KEGG pathway terms (p-value < 0.05) relevant to interactional genes were also identified by the R package clusterProfiler (Wu et al., 2021).

Identification of Transcription Factors and MiRNAs Interacting With Interactional Hub Genes

To study underlying regulatory mechanisms at the transcriptional level and identify hub protein’s regulatory molecules, a comprehensive network-based method was employed to decipher regulatory TFs and miRNAs. ChIP-X Enrichment Analysis 3 (ChEA3, https://amp.pharm.mssm.edu/chea3/) (Keenan et al., 2019) is an open-access tool for identifying TFs that are in control of observed alters in gene expression. ChEA3 assembles TFs data from ENCODE (Davis et al., 2018), ReMap (Chèneby et al., 2020), Genotype-Tissue Expression (GTEx) (Carithers and Moore, 2015), ARCHS4 (Lachmann et al., 2018), and Enrichr (Kuleshov et al., 2016). Evidence mainly from ChIP-seq experiments assures the interaction between TFs and customized input genes. The integration method MeanRank performed the best in the ChEA3 benchmark and was recommended as the preferable method to present results from different TFs libraries. Results from each library were sorted in ascending order by score using ChEA3. The lower the score, the closer the connection between a gene and a TF. In the study, TFs of COVID-19/NSCLC interactional hub genes were retrieved from ChEA3. Moreover, top 10 TFs with the lowest scores were selected for each gene after removing duplication. Experimentally supported data about the relation of miRNAs to COVID-19/NSCLC interactional hub genes were obtained from Tarbase (Karagkouni et al., 2018) and miRTarbase databases (Chou et al., 2018). To assure the accuracy and robustness of results, we searched the Tarbase database and selected miRNAs that were identified by low-throughput experimental techniques. In the miRTarbase database, the source of miRNAs included a reporter assay, western blot, and qt-PCR. Species was set to Homo sapiens in these two databases. Subsequently, regulatory networks of TFs, miRNAs, and COVID-19/NSCLC interactional hub genes were also constructed by Cytoscape 3.9.0 software (Shannon et al., 2003).

Evaluation of Applicant Drugs for COVID-19 and Non-Small-Cell Lung Cancer

Evaluation of protein–drug interactions is a crucial strategy to detect structural features and respond to receptor sensitivity of proteins, which can pave the way for drug development. To explore promising drugs for COVID-19/NSCLC, drug–proteins interactions were retrieved from the Enrichr database (https://maayanlab.cloud/Enrichr/) (Kuleshov et al., 2016) and the ShinyGO v0.75 web tool (http://bioinformatics.sdstate.edu/go/) (Ge et al., 2020).

Enrichr, containing 192 gene-set libraries and counting, is a comprehensive web portal for gene-set enrichment analysis (Kuleshov et al., 2016). Two series of Enrichr library were downloaded and analyzed, including drug perturbations from GEO (DPFG) and a collection of Drug Signatures Database (DSigDB). DPFG was built based on experimentally supported data from GEO in which gene expression levels were assessed before and after drug administration. Two types of data were retrieved from DPFG for subsequent analyses, including data on the relation of a drug to upregulated genes and data on the relation of a drug to downregulated genes. Species was set to Homo sapiens for filtering data. DSigDB is a collection of 22,527 drug-related gene sets with 17,389 compounds covering genes (Yoo et al., 2015). ShinyGO v0.75 is a graphical tool for gene-set enrichment analysis, which is accessible to KEGG and STRING. It applies false discovery rate (FDR) as an adjustment method for p-values of enrichment terms. To obtain robust protein–drug interactions, we merely selected drug molecules which targeted more than half of COVID-19/NSCLC interactional hub genes. Enrichment terms from ShinyGO v0.75 with an adjusted p-value <0.05 were included. Finally, drug molecules obtained from Enrichr libraries and ShinyGO v0.75 were carefully searched on PubChem and PubMed databases. Drug molecules with anticancer and/or antivirus effect were finally selected by browsing relevant articles carefully.

Result

Identification of Interactional Genes Between COVID-19 and Non-Small-Cell Lung Cancer

COVID-19–related genes retrieved from OMIM, TTD, PubChem, DisGeNET, GeneCards, and CTD were 2, 78, 628, 1,832, 2,572, and 9,860, respectively (Supplementary File S1). After removing duplication, a total of 11,392 related genes of COVID-19 were obtained from these online databases. COVID-19–related DEGs obtained from the RNA-seq dataset GSE147507, GSE157103, and GSE166190 were 4,825 (3,064 upregulated and 1,761 downregulated), 3,089 (1,422 upregulated and 1,667 downregulated), and 494 (278 upregulated and 216 downregulated), respectively (Table 1; Supplementary File S2). By intersecting COVID-19–related genes from online databases and DEGs from RNA-seq datasets, we obtained 3,795 COVID-19–related genes (Figure 2A).

TABLE 1

SourcePlatformGroup and sample countDEGs upregulatedDEGs downregulatedDEGs allTotal DEG upregulatedTotal DEG downregulated
GSE147507GPL18573A549-ACE2 SARS-CoV-2 (6) vs. A549-ACE2 mock (6)18425622,40430641761
A549 SARS-CoV-2 (6) vs. A549 mock (13)225236461
Calu3 SARS-CoV-2 (3) vs. Calu3 mock (3)1,2936801973
COVID-19 lung biopsy (2) vs. Healthy lung biopsy (2)385456841
NHBE SARS-CoV-2 (3) vs. NHBE mock (7)6949118
GSE157103GPL24676Female-COVID-ICU (17) vs. female-non COVID-ICU (7)4542286821,4221,667
Female-COVID-non-ICU (21) vs. female-non-COVID-non-ICU (6)5285391,067
Male-COVID-ICU (33) vs. male-non-COVID-ICU (8)587104691
Male-COVID-non-ICU (29) vs. male-non-COVID-non-ICU (4)8231,3542,177
GSE166190GPL20301Adult-COVID-19-positive vs. Adult-COVID-negative242202444278216
Child-COVID-19-Positive vs.child-COVID-19-negative401454
TCGA-NSCLCIllumina HiSeqCancer (1,027) vs. normal (108)8,2292,13910,3688,2292,139

DEGs in different RNA-seq datasets.

FIGURE 2

FIGURE 2

Identification of COVID-19/NSCLC interactional genes by intersecting COVID-19–related genes and NSCLC-related genes from public databases and DEGs from RNA-seq datasets. (A) Identification of COVID-19-related genes. (B) Identification of NSCLC-related genes. (C) Identification of COVID-19/NSCLC interactional genes.

NSCLC-related genes collected from OMIM, CTD, TTD, DisGeNET, and GeneCards were 15, 144, 167, 447, and 4,536, respectively (Supplementary File S3). After removing duplication, a total of 4,725 related genes of NSCLC were obtained. A total of 10,368 DEGs (8,229 upregulated and 2,139 downregulated) of NSCLC were identified from the TCGA (Supplementary File S4). By intersecting NSCLC-related genes from online platforms and DEGs from the TCGA, we obtained 618 NSCLC-related target genes (Figure 2B). Finally, a total of 122 COVID-19/NSCLC interactional genes were identified by intersecting COVID-19–related genes and NSCLC-related genes (Figure 2C).

Systemic Biological Significance of Interactional COVID-19/Non-Small-Cell Lung Cancer Genes

Systemic biology provides multivariate approaches to analyze the larger interactive network of biological pathways holistically and identify important players in disease onset and progression (Starchenko and Lauffenburger, 2018). To reveal characteristics shared by COVID-19 and NSCLC at the genetic level based on systemic biology, GO and KEGG pathway enrichment analyses were performed based on 122 COVID-19/NSCLC interactional genes. As a result, 392 GO terms (BP: 352; CC: 26; and MF: 14) and 23 KEGG pathways were highlighted. The top 10 GO terms of each ontology and 23 KEGG pathways are shown in Figure 3. Representative BP terms included cell aging, regulation of mitotic cell cycle phase transition, and regulation of cell cycle phase transition; representative CC terms included condensed chromosome kinetochore, kinetochore, and condensed chromosome; representative MF terms included cytokine activity, receptor ligand activity, and DNA binding and bending. In addition, representative pathways included the Janus kinase/signal transducers and activators of transcription (JAK/STAT) signaling pathway, interleukin (IL)-17 signaling pathway, chemical carcinogenesis receptor activation, cell cycle, and cellular senescence. These GO terms and KEGG pathways may exert a synergistic effect in the morbidity of COVID-19/NSCLC, which could be clues of therapeutic strategies for these two diseases.

FIGURE 3

FIGURE 3

Functional annotation of COVID-19/NSCLC interactional genes. (A) GO analysis of COVID-19/NSCLC interactional genes. (B) KEGG pathway analysis of COVID-19/NSCLC interactional genes. Note: rich factor is defined as the ratio of input genes that are annotated in a term to all genes that are annotated in this term. The computational formula of rich factor is as follows: Rich factor = number of input genes under this pathway term/number of all annotated genes under this pathway term. The greater the rich factor, the greater the degree of pathway enrichment.

Identification of Hub Genes and Exploration of Their Interactions

To explore interrelationships of 122 COVID-19/NSCLC interactional genes, a PPI network containing 59 nodes and 347 edges with the highest confidence scores has been visualized (Figures 4A, 5A). Key parameters of the network were as follows: mean betweenness: 0.05, mean degree: 11.7, and mean closeness: 0.59. A total of 25 DEGs with a degree twice the median or more were identified by the PPI network (Figure 4B). In addition, the other 11 analytical network approaches were employed to identify the top 10 genes for each approach (Table 2). After intersecting genes identified by different approaches, a total of 21 COVID-19/NSCLC interactional hub genes were finally obtained for subsequent analyses. The detail information of these 21 genes is shown in Table 3.

FIGURE 4

FIGURE 4

PPI network analysis based on COVID-19/NSCLC interactional genes. (A) PPI network containing 59 nodes and 347 edges. (B) Bubble chart of the genes with degree value more than the two-fold median degree value in the whole network. Note: nodes represent interactional genes, and edges represent interaction relationships in panel (A).

FIGURE 5

FIGURE 5

PPI network analysis and sub-network analysis based on identified COVID-19/NSCLC interactional genes. (A) PPI network containing 59 nodes and 347 edges. (B) Biggest sub-network of the PPI network. (C) Small sub-network of the PPI network. Note: nodes represent interactional genes, and edges represent interaction relationships. The depth of the color of the node is positively correlated with the degree value in panel (A).

TABLE 2

MethodGeneGene count
BetweennessAURKB, BIRC5, CCNA2, CCNB2, CDC20, FOXM1, MYBL2, PLK1, TOP2A, and UBE2C10
BottleneckBIRC5, CCNA2, CDC20, FOXM1, MYBL2, PLK1, TOP2A, and UBE2C8
ClosenessAURKB, BIRC5, BUB1B, CCNA2, CCNB2, CDC20, DLGAP5, PLK1, TOP2A, and TPX210
DegreeAURKB, BIRC5, BUB1B, CCNA2, CCNB2, CDC20, DLGAP5, PLK1, TOP2A, and TPX210
DMNCBIRC5, HJURP, KIF2C, KIF4A, NCAPH, NDC80, NEK2, SPAG5, TTK, and UBE2C10
EccentricityBIRC5, CCNA2, CCNB2, CDKN3, E2F8, FOXM1, NEK2, PLK1, RAD51, and UHRF110
EPCAURKB, BUB1B, CCNA2, CCNB2, CDC20, DLGAP5, TOP2A, TPX2, TTK, and UBE2C10
MCCAURKB, BUB1B, CCNA2, CCNB2, CDC20, DLGAP5, NDC80, TOP2A, TPX2, and TTK10
MNCAURKB, BUB1B, CCNA2, CCNB2, CDC20, DLGAP5, RRM2, TOP2A, TPX2, and TTK10
RadialityAURKB, BIRC5, BUB1B, CCNA2, CCNB2, CDC20, DLGAP5, PLK1, TOP2A, and TPX210
StressAURKB, BIRC5, BUB1B, CCNA2, CCNB2, CDC20, FOXM1, PLK1, TOP2A, and UBE2C10

Top 10 important genes identified by different analytical methods.

TABLE 3

Hub geneDescriptionEnsembl gene IDEntrezGene typeChrPosition (Mbp)
CDC20Cell division cycle 20ENSG00000117399991Protein coding143.358981
KIF2CKinesin family member 2CENSG0000014294511,004Protein coding144.739818
NEK2NIMA-related kinase 2ENSG000001176504,751Protein coding1211.658657
RRM2Ribonucleotide reductase regulatory subunit M2ENSG000001718486,241Protein coding210.120698
NCAPHNon-SMC condensin I complex subunitENSG0000012115223,397Protein coding296.335766
HJURPHolliday junction recognition proteinENSG0000012348555,355Protein coding2233.833416
CCNA2Cyclin A2ENSG00000145386890Protein coding4121.816444
TTKTTK protein kinaseENSG000001127427,272Protein coding680.003887
FOXM1Forkhead box M1ENSG000001112062,305Protein coding122.85768
DLGAP5DLG-associated protein 5ENSG000001267879,787Protein coding1455.148112
BUB1BBUB1 mitotic checkpoint serine/threonine kinase BENSG00000156970701Protein coding1540.161023
CCNB2Cyclin B2ENSG000001574569,133Protein coding1559.105126
PLK1Polo-like kinase 1ENSG000001668515,347Protein coding1623.677656
AURKBAurora kinase BENSG000001789999,212Protein coding178.204733
SPAG5Sperm-associated antigen 5ENSG0000007638210,615Protein coding1728.577565
TOP2ADNA topoisomerase II alphaENSG000001317477,153Protein coding1740.388525
BIRC5Baculoviral IAP repeat-containing 5ENSG00000089685332Protein coding1778.214186
NDC80NDC80 kinetochore complex componentENSG0000008098610,403Protein coding182.571557
TPX2TPX2 microtubule nucleation factorENSG0000008832522,974Protein coding2031.739271
UBE2CUbiquitin-conjugating enzyme E2 CENSG0000017506311,065Protein coding2045.812576
KIF4AKinesin family member 4AENSG0000009088924,137Protein codingX70.290104

COVID-19/NSCLC interactional hub genes.

To further understand the main systemic biological significance of the network, a sub-network analysis based on nodes of the whole network was conducted. With the recommended parameters in the molecular complex detection (MCODE) module (degree cutoff: 2, node score cutoff: 0.2, K-Core: 2, and maximum depth = 100), two sub-networks were extracted from the primary one (Figures 5B,C). Based on the biggest sub-network containing 23 nodes and 232 edges (Figure 5B), 283 GO terms (BP: 235; CC: 35; and MF: 13) and 12 KEGG pathways were highlighted (Supplementary File S5). The top 10 terms of each GO and all the KEGG pathways are shown in Figures 6A–D. Typical BP terms mainly involved mitotic spindle organization, organelle fission, and chromosome segregation; typical CC terms mainly involved chromosomal region, midbody, and spindle; typical MF terms mainly involved protein threonine kinase activity, microtubule binding, and kinase regulator activity. In addition, the main KEGG pathways involved cell cycle, viral carcinogenesis, p53 signaling pathway, and cellular senescence. By comparing the KEGG pathways from the whole PPI network and the main sub-network, functions of COVID-19/NSCLC interactional genes in cell cycle, cellular senescence, and p53 signaling pathway were underlined.

FIGURE 6

FIGURE 6

GO and KEGG pathway analysis based on the interactional genes in the biggest sub-network. (A) Top 10 BP terms of GO analysis. (B) Top 10 CC terms of GO analysis. (C)Top 10 MF terms of GO analysis. (D) Twelve pathway terms of KEGG analysis. Note: nodes represent genes or pathways, and edges represent enrichment relationships. The size of term nodes is positively correlated with the number of enriched genes. The size of gene nodes is positively correlated with the number of enriched terms.

Determination of Regulatory Signatures

A network-based approach was employed to identify regulatory TFs and miRNAs and to have a glimpse of changes happening at the transcriptional level. A total of 1,632 TFs were filtered from the ChEA3 database and sorted by the MeanRank method for the most robust result (Supplementary File S6). The data source of TFs of each COVID-19/NSCLC interactional hub gene is shown in Figure 7A. The top 10 TFs of the interactional hub gene recommended by the ChEA3 database were retrieved for subsequent analyses (Figure 7B). Next, a TF–gene interaction network containing 30 nodes and 186 edges was constructed (Figure 7C). Centromere protein A (CENPA), DNA methyltransferase 1 (DNMT1), MYB proto-oncogene-like 2 (MYBL2), transcription factor Dp-1 (TFDP1), and zinc finger protein 367 (ZNF367) were identified as the most influential regulatory factors since they targeted all interactional hub genes.

FIGURE 7

FIGURE 7

Determination of regulatory signatures (TFs). (A) Data source of TFs of the COVID-19/NSCLC interactional hub gene. (B) Top 10 TFs identified by the ChEA3 database. (C) TFs–interactional hub genes interaction network. Note: red nodes represent TFs, green nodes represent genes, and edges represent interaction relationships in panel (C).

Furthermore, we obtained 17 and 40 miRNAs from Tarbase and mirTarbase, respectively (Supplementary File S7). After removing duplicated outcomes, a total of 44 miRNAs engaged with interactional hub genes of COVID-19 and NSCLC were identified. A miRNA–gene interaction network comprising 53 nodes and 51 edges was built (Figure 8). Representative miRNAs included hsa-miR-24-3p (targeting at four interactional hub genes), hsa-miR-16-5p (targeting at two interactional hub genes), hsa-let-7a-5p (targeting at two interactional hub genes), hsa-miR-34a-5p (targeting at two interactional hub genes), and hsa-miR-10b-3p (targeting at two interactional hub genes). Taken together, these two networks indicated that 10 important TFs and 44 post-transcriptional regulatory signatures (miRNAs) provided evidence for exploring regulatory mechanisms of COVID-19 and NSCLC by participating in the regulation of interactional hub genes.

FIGURE 8

FIGURE 8

MiRNAs–interactional hub genes interaction network. Note: purple nodes represent genes, orange nodes represent MiRNAs, and edges represent the interactions between nodes.

Identification of Candidate Drugs for COVID-19 and Non-Small-Cell Lung Cancer

To identify prospective agents for COVID-19 and NSCLC, we selected potential candidates from the Enrichr database, including DSigDB, DPFG, and ShinyGO v0.75. Studies of drug molecules from these three sources were retrieved from the PubMed database for screening drugs with anticancer and/or antivirus effect (Supplementary Table S1). According to the result, the number of drug molecules obtained from DSigDB, DPFG, and ShinyGO v0.75 was 15, 10, and 3, respectively. A total of 23 candidates were identified after removing duplication. Among them, nine drugs have already been used in NSCLC, while six drugs have been registered to clinical trials for COVID-19. Resveratrol was identified both in DSigDB and DPFG, lucanthone was identified both in DSigDB and ShinyGO v0.75, and vemurafenib was identified both in DPFG and ShinyGO v0.75. Notably, natural small molecules included resveratrol and quercetin worth focus because of their binding affinity with molecular ACE-2.

Discussion

The present outbreak of SARS-CoV-2 has spread to 224 countries, with more than 250 million people infected (Worldometer, 2021). The elderly and those with pre-existing complications or smoking-related lung damage are susceptible to COVID-19, while they are precisely the profile of lung cancer patients (Berlin et al., 2020; Liao et al., 2020). The emergence of COVID-19 has resulted in detrimental alterations including (epi-)genetic alterations and changes at the transcriptional level and related signaling pathways, which have changed lung cancer treatment provision (Alnajeebi et al., 2021). Various candidate drugs and vaccines are being studied for the treatment of COVID-19; nevertheless, no therapy has shown a specific effect (Majumder and Minko, 2021). Therefore, in this study, we conducted a comprehensive systemic biological and bioinformatic analysis to identify potential interactional biomarkers that might help understand co-pathogenic factors of these two diseases and provide clues for the treatment of COVID-19 and NSCLC.

Identified Related Genes Indicated Genetic Interrelationships of COVID-19 and Non-Small-Cell Lung Cancer

Bioinformatic and systemic biological analysis can integrate data of gene expressions and protein interactions, analyze gene-regulatory pathways, and construct protein–protein networks, which is an effective tool for molecular mechanisms (Durmuş et al., 2015; Ahmed, 2020) and drug development (van Driel and Brunner, 2006), especially in the field of cancer treatment (Lee et al., 2017; Wooller et al., 2017) and infectious diseases (Josset et al., 2010; Li J. et al., 2016). NSCLC is a genomic disease which involves a loss of control over vital cellular functions and could be deteriorated by COVID-19 infection due to differential expressions of multiple host genes (Ahmed, 2020). Biomarker identification is helpful in managing NSCLC or COVID-19 (Samprathi and Jayashree, 2020; Wang et al., 2021). Inherent cancer-related changes in NSCLC genome and the aberrant expression of host factors caused by COVID-19 infection could interact with each other and impact the prognosis of patients (Passaro et al., 2021). In this study, we identified 122 COVID-19/NSCLC interactional genes through bioinformatic and systemic biological approach, which would pave the way to reveal common pathophysiological mechanisms shared by COVID-19 and NSCLC (Ahmed, 2020; Samprathi and Jayashree, 2020; Wang et al., 2021).

Gene Ontology and Kyoto Encyclopedia of Genes and Genomes Analyses Highlighted Systemic Biological Significance

To elucidate the roles of the identified 122 COVID-19/NSCLC interactional genes, GO and KEGG pathway analyses were conducted. In GO analysis, biological processes were mostly enriched in sister chromatid segregation, regulation of cell aging, spindle organization, and mitotic nuclear division. Dysregulated cell–cycle machinery, uncontrolled proliferation, and resistance to programmed cell death feature cancers, including NSCLC (Sherr, 1996; Hanahan and Weinberg, 2000; Smolle et al., 2019). Due to the paucity of well-proven therapies or vaccines, inhibiting cell division may be a good way to control COVID-19, which shares a similar goal with NSCLC treatment (Borcherding et al., 2020). Some of interactional genes participated in the BP have been found to be closely associated with the development and prognosis of NSCLC and COVID-19. For example, aurora kinase B (AURKB) has been found among the DEGs of SARS-CoV-2 in Caco-2 cells (Bock and Ortea, 2020). In addition, an in vitro experiment suggested that AURKB inhibition caused cell cycle arrest and polyploidy, followed by extensive cell death in the NSCLC model (Bertran-Alamillo et al., 2019). During SARS-CoV-2 infection, the PLK1 gene and its translational products performed significant tasks in the normal cell division cycle, promoting the spread of virus (Chen et al., 2021). The inhibition of PLK1 could alter the immune microenvironment of NSCLC by increasing dendritic cell maturation and enriching T-cell infiltration (Zhou et al., 2021).

KEGG links a group of genes in the genome with the higher-order functional information via computerizing current knowledge about cellular processes and standardizing gene annotations (Kanehisa and Goto, 2000). Based on KEGG pathway analysis, we found pathways shared by COVID-19 and NSCLC, including cell cycle, celluar senescence, IL-17 signaling pathway, chemical carcinogenesis-receptor activation, p53 signaling pathway, and Janus kinase/signal transducer and activator of the transcription (JAK-STAT) signaling pathway. Previous studies suggested that JAK-STAT, IL-17, and p53 signaling pathways were not only closely associated with the tumorigenesis of NSCLC via regulating cell proliferation, differentiation, and apoptosis but also involved in the promotion of cytokine storm which deteriorated conditions in COVID-19 patients (Shibabaw, 2020; Hu H. et al., 2021; Satarker et al., 2021). Regulated by interferon (IFN) signaling, the endocytosis of SARS-CoV-2 in lung epithelial cells is initiated by the ACE2 receptor. A cross-sectional study showed that in severe COVID-19 patients, dysregulated type I/III IFNs and JAK/STAT signaling caused impaired antiviral responses (Hadjadj et al., 2020). During the development of COVID-19, SARS-CoV-2 infection induced an excessive immune response and released a variety of pro-inflammatory cytokines through the JAK/STAT pathway, such as IL-2, IL-6, and granulocyte colony-stimulating factor signaling (Solimani et al., 2021). In addition, high levels of pro-inflammatory cytokines and chemokines accelerated disease progression, thus becoming the main risk factor of cellular senescence and age-related diseases (Bruunsgaard et al., 2003; Hu X. et al., 2021). Previous studies showed that the severity of COVID-19 and NSCLC was closely associated with hyperinflammation that could drive lung and multiorgan injury and increase mortality via cytokine storm and sepsis (Gomes et al., 2014; Tay et al., 2020). Thus, targeting inflammation-related signaling pathways may be a good option for the treatment of COVID-19 and NSCLC.

Among the interactional genes involved in the aforementioned pathways, cyclin A2 (CCNA2) and cyclin B2 (CCNB2) are members of the cyclin protein family, which play critical roles in controlling cell cycle, cell senescence, and viral infection (Nam and van Deursen, 2014). Overexpressed CCNA2 and CCNB2 are correlated with a poor prognosis in NSCLC (Cooper et al., 2009; Takashima et al., 2014; Qian et al., 2015). Meanwhile, during SARS-CoV-2 infection, the increasing level of CCNA2 not only provides favorable conditions for virus transmission but also hints at chromosomal abnormalities and other genetic material damage in the host cells (Chen et al., 2021). Cellular senescence, an irreversible state of cell cycle arrest in response to damaging stimuli, secretes factors known as senescence-associated secretory phenotype or SASP (González-Gualda et al., 2021; Kumari and Jat, 2021). In this case, cells maintain active metabolism without responses of mitogenic and apoptotic signals (González-Gualda et al., 2021). In the elderly or patients upon consistent and chronic damages, cellular senescence may impair regenerative ability and accelerate the progression of inflammation and lung cancer (Hernandez-Segura et al., 2018; Kuźnar-Kamińska et al., 2018). In lung cancer tumorigenesis, senescent cells produced soluble signaling factors [including interleukins (Rabinovich et al., 2007), chemokines (Coppé et al., 2010), and growth factors (Yang et al., 2006)], proteases (Finkel et al., 2007), insoluble proteins, and extracellular matrix components (Acosta et al., 2008a; Acosta et al., 2008b), which mediated cell proliferation, invasion, and migration (Wald et al., 2011; Han et al., 2015; Kuźnar-Kamińska et al., 2016). Furthermore, a strong connection between cellular senescence and SARS-CoV-2 has been found, as cellular senescence increased the risk of developing severe COVID-19 (Nehme et al., 2020). SARS-CoV-2 infection could induce paracrine senescence by increasing IFN secretion in infected cells, and danger-associated molecular patterns released in cells undergoing pyroptosis or necroptosis further accelerated senescence in the environment via SASP factors (Kim et al., 2009; Acosta et al., 2013; Kandhaya-Pillai et al., 2017; Nehme et al., 2020). Cellular senescence, in turn, 1) led to a weak adaptive immunity by increasing senescent-like T cells and B lymphocytes (Frasca et al., 2017; De Biasi et al., 2020); 2) enhanced aberrant healing response and tissue fibrosis in the respiratory system (Mason, 2020); 3) resulted in vascular dysfunction via decreasing angiogenesis and increasing thrombosis and inflammatory responses in COVID-19 patients (Ungvari et al., 2018; Escher et al., 2020). Dysregulation of the cell cycle and cellular senescence has been observed in both COVID-19 and lung cancer. Thus, the importance of understanding the homeostasis maintenance and pathological alteration of these biological processes should be emphasized.

In summary, the aforementioned genes and pathways are commonly and significantly altered in both COVID-19 and NSCLC, preparing the ground for pathophysiological studies and drug development.

Protein–Protein Interaction and Sub-Network Analysis Explored In-Depth Interrelationships of Interactional Hub Genes

On the basis of interactional genes, we built a PPI network presenting in-depth biological characteristics. In the biggest sub-network, CCNA2, CCNB2, AURKB, DNA topoisomerase II alpha (TOP2A), and baculoviral IAP repeat containing 5 (BIRC5) were closely connected with other genes. TOP2A, a target for cytotoxic drugs (etoposide, anthracyclines) and a key regulator of chromosome condensation and chromatid separation, prevents DNA replication and transcription (Nicoś et al., 2021). It is generally believed that TOP2A is a prognostic indicator of NSCLC, and its level is negatively correlated with the prognosis (Hou et al., 2017; Ma et al., 2019). In consistence with the result of our study, TOP2A has been identified as a hub gene that could govern many cellular processes by protein–protein interactions in COVID-19 patients (Hasan et al., 2022).

Comparing to the result of KEGG analysis based on 122 interactional genes, the biggest sub-network identified additional pathways, such as forkhead box, sub-group O (FoxO) signaling pathway, ubiquitin-mediated proteolysis, viral carcinogenesis, and apoptosis. As tumor suppressors, FoxO proteins increased the expression level of death receptor ligands such as Fas ligand and tumor necrosis factor (TNF) apoptosis ligand, engaged with pro-apoptotic pathways, and blocked cell cycle progression (Farhan et al., 2017). According to the result, viral carcinogenesis was one of the most important KEGG pathways, indicating that co-pathogenic factors shared by COVID-19 and NSCLC could be the basis of finding potential drugs for synergistic treatment.

Interaction of Transcription Factors and MicroRNAs With Interactional Hub Genes

Genetic alterations have long been blamed for malignancy but cannot fully explain tumor development. Epigenetic dysregulations give new insights into how heritable changes in gene expression happen without involving changes in nucleotide sequence, and finally promote carcinogenesis (Yuan et al., 2016). MiRNAs and TFs are the largest families of trans-acting gene regulatory species and pivotal players in a complex regulatory network (Sharma et al., 2020). Herein, we investigated TF–gene and miRNA–genes interaction that might help in learning more about disease development.

Transcription factors are key cellular components that regulate multiple genes over a long distance by maintaining proliferation status, driving cellular differentiation, and determining cell fate (Essebier et al., 2017). The identified top 10 TFs included CENPA, DNMT1, MYBL2, TFDP1, ZNF367, high mobility group AT-hook 1 (HMGA1), proliferation-associated 2G4 (PA2G4), high mobility group AT-hook 2 (HMGA2), zinc finger protein 878 (ZNF878), and FOXM1. Among them, CENPA (Mullen et al., 2020), DNMT1 (Wu et al., 2020), MYBL2 (Mullen et al., 2020), TFDP1 (Zhan et al., 2017), ZNF367 (Liu Z. et al., 2018), HMGA1 (Zhang Z. et al., 2015), HMGA2 (Gao et al., 2018), and FOXM1 (Mullen et al., 2020) were previously shown to be highly associated with NSCLC progression. It has been demonstrated that CENPA, MYBL2, and FOXM1 were linked to numerous cancer-specific enhancers, and their elevated expression levels were associated with a poor survival rate of NSCLC patients (Mullen et al., 2020). Meanwhile, the level of DNMT1 was found significantly downregulated in SARS-COV-2–infected epithelial cells (Muhammad et al., 2021). SARS-COV-2 infection mainly affects the molecular mechanisms of aging centered on HMGA1 and HMGA2 proteins, and their interactions may impair or trigger inflammatory pathways, leading to various responses in different age groups (Mercatelli et al., 2021). To conclude, DNMT1, HMGA1, and HMGA2 might be central TFs in the TF–gene regulatory network in COVID-19/NSCLC.

MiRNAs are endogenous small non-coding RNA molecules that can regulate the expression of non-coding sequences and genes involved in oncogenesis (Liu et al., 2014). Various studies have identified miRNAs as key players in the pathogenesis and therapeutics of viral diseases. Moreover, as part of host–pathogen interactions, miRNA can scan target SARS-CoV-2 genes as well as host inflammatory machinery to counter-act the impairing effects of infection (Ghosh et al., 2009; Khokhar et al., 2022). The top 10 most significant miRNAs mainly involved in respiratory diseases (Table 4), including NSCLC [hsa-miR-24-3p (Wei et al., 2021), hsa-miR-34a-5p (Zhou X. et al., 2020), hsa-miR-10b-3p (Liu et al., 2017), hsa-miR-20a-5p, and hsa-miR-17-5p (Leidinger et al., 2016)], idiopathic pulmonary fibrosis [hsa-miR-524-5p (Li Q. et al., 2020), COVID-19 (hsa-miR-20a-5p, hsa-miR-17-5p (Li C. et al., 2020)), and hsa-miR-16-5p (Kim et al., 2020)], pneumoconiosis [hsa-let-7a-5p (Zhang et al., 2018)], and asthma [(hsa-let-7a-5p (Huang et al., 2019)]. Interestingly, compared with the healthy controls, the levels of hsa-miR-20a-5p and hsa-miR-17-5p were significantly downregulated in patients with COVID-19; therefore, they were considered as essential modulators of viral replication (Li C. et al., 2020). In addition, according to a study based on high-throughput qRT-PCR validation, hsa-miR-20a-5p and hsa-miR-17-5p were the top markers that could distinguish NSCLC patients from unaffected controls with 94.5% accuracy (Leidinger et al., 2016).

TABLE 4

MiRNA nameRespiratory disease typeReferences
hsa-miR-24-3pNSCLCWei et al. (2021)
hsa-miR-16-5pCOVID-19Kim et al. (2020)
hsa-let-7a-5pPneumoconiosis and asthmaZhang et al. (2018), Huang et al. (2019)
hsa-miR-34a-5pNSCLCXingni Zhou et al. (2020)
hsa-miR-10b-3pPneumoconiosis and asthmaLiu et al. (2017)
hsa-miR-20a-5pNSCLC and COVID-19Leidinger et al. (2016), Caixia Li et al. (2020)
hsa-miR-17-5pNSCLC and COVID-19Leidinger et al. (2016), Caixia Li et al. (2020)
hsa-miR-524-5pIdiopathic pulmonary fibrosisQilong Li et al. (2020)
hsa-miR-376a-5pHas not been reportedNot available
hsa-miR-483-3pHas not been reportedNot available

Top 10 miRNAs involved in various respiratory diseases.

Notably, hsa-miR-24-3p, hsa-miR-16-5p, and hsa-let-7a-5p targeted AURKB, central regulator of cell division (Bertran-Alamillo et al., 2019), indicating that these miRNAs may be potential targets to control cell cycle and cellular senescence. In addition, miR-34a-5p regulated the G1/S checkpoint in NSCLC cells (Gupta et al., 2020). Forced expression of miR-34a-5p enhanced p21 expression and promoted cellular senescence, whereas downregulated miR-34a-5p decreased senescence and increased apoptosis by targeting B-cell lymphoma-2 (BCL2), myelocytomatosis oncogene (MYC), mesenchymal-epithelial transition (MET), and p53 (Gallardo et al., 2009; Kasinski and Slack, 2012). According to the result, BIRC5, FOXM1, CCNA2, and PLK1 showed a higher degree of interaction with miRNAs, which may play crucial roles in the interaction networks.

To sum up, TF–gene intereactions are reactors that regulate gene expression by binding with target genes and miRNAs (Zhang HM. et al., 2015). The constructed TF–gene and miRNA–gene interaction networks help further understand the direct regulatory relationship of miRNA and TFs in NSCLC and COVID-19, while interactional hub genes may be crucial biomarkers and therapeutic targets.

Potential Drugs Provided Possible Treatments for COVID-19 and Non-Small-Cell Lung Cancer

In this study, a total of 23 candidates containing 18 drugs and 5 natural compounds were identified. Among them, nine drugs have already been administered to NSCLC patients, including irinotecan (Murakami et al., 2018), dasatinib (Xu et al., 2015), 5-fluorouracil (Nokihara et al., 2017), etoposide (Liang et al., 2017), carboplatin (Hasegawa et al., 2020), ascorbic acid (Pathak et al., 2005), azacitidine (Cheng et al., 2021), decitabine (Chu et al., 2013), and palbociclib (Nie et al., 2019). According to the ClinicalTrials database (https://clinicaltrials.gov/), cyclosporin A (NCT04540926), dasatinib (NCT04830735), quercetin (NCT04851821), etoposide (NCT04356690), ascorbic acid (NCT04401150), and decitabine (NCT04482621) have been registered to clinical trials for evaluating their efficacy in COVID-19 patients. Therefore, it can be inferred that dasatinib, etoposide, ascorbic acid, and decitabine are likely to be the promising agents for NSCLC patients with COVID-19, which provides support for our predictions.

Quercetin and resveratrol are two natural compounds which have both antiviral and anti-NSCLC effects (Zhao et al., 2010; Wu et al., 2015; Lin et al., 2017; Chai et al., 2021). Quercetin, a well-known natural polyphenol with anti-inflammatory, antioxidant, and immunomodulatory properties, is involved in a variety of diseases, such as viral infections, respiratory diseases, allergies, asthma, and cancer (Li Y. et al., 2016; Xu et al., 2019). Previous network pharmacology and molecular docking studies anticipated that quercetin could interfere with SARS-CoV-2 replication by interacting with 3-chymotrypsin-like protease (3CLpro), papain-like protease (PLpro), and S proteins (Derosa et al., 2021; DI Pierro et al., 2021b). An open-label, randomized controlled trial included 152 outpatients with confirmed SARS-CoV-2 infection but without severe COVID-19 symptoms found that formulated quercetin treatments reduced the frequency and length of hospitalization, the need of non-invasive oxygen therapy, and the number of deaths (DI Pierro et al., 2021a). In addition, available experimental studies suggested that quercetin could directly modulate multiple lung cancer-relevant miRNAs and DNA methylation (Kim DH. et al., 2019; Kedhari Sundaram et al., 2019). Another ingredient resveratrol is a potent antioxidant that could inhibit platelet aggregation and vasodilation, reduce blood viscosity, and maintain blood flow (Olas and Wachowicz, 2005; de la Lastra and Villegas, 2007). In vitro and in vivo studies indicated that resveratrol induced cell apoptosis and inhibited proliferation, growth, and metastasis in NSCLC (Yousef et al., 2017; Tang et al., 2020). Given the thrombin-inhibitory and anti-inflammatory effects, resveratrol deserves further study for the treatment of COVID-19 and NSCLC.

It is noteworthy that a computational structure-based study identified trichostatin A as a potential SARS-CoV-2 Mpro inhibitor. The result has been further validated by an essay which has suggested that trichostatin A could reduce the viral RNA load, viral antigen expression, and infectious virus particle formation (Wen et al., 2021). Other identified ingredients also had an anti-respiratory virus effect. For instance, cyclosporine at non-cytotoxic concentrations could induce a strong inhibition of the replication of specific coronaviruses in vitro, including SARS-CoV, MERS-CoV, and human coronavirus 229E (HCoV-229E) (Poulsen et al., 2020). A retrospective observational study recently found that in COVID-19 hospitalized patients, cyclosporine was significantly associated with a decrease in mortality, probably due to its combined activity of immunosuppression and antiviral activity (Guisado-Vasco et al., 2020). Dyall et al. (2014) found that dasatinib was active against both MERS-CoV and SARS-CoV in vitro, and it might minimize immunotoxicity as it blocked viral replication. Vemurafenib interfered the cellular Raf/MEK/ERK signaling cascade by binding to the ATP-binding site of BRAF (V600E) kinase and inhibiting its function (Pleschka et al., 2001). Interestingly, since Raf/MEK/ERK signaling pathways mediated the increasing SARS-CoV-1 replication (Dyall et al., 2014), it may be a therapeutic target for host-directed SARS-CoV-2 antivirals. Thus, vemurafenib may also be a potential anti-COVID-19 drug. To conclude, according to available literature and the result of our study, identified molecules may be prospective agents in the treatment of COVID-19 and NSCLC.

Based on the comprehensive bioinformatic and systemic biological analysis, the present study carried out a study framework to reveal interaction networks and therapeutic implications for NSCLC patients with COVID-19, the findings of which would provide evidence and shed light for the further research on COVID-19/NSCLC. However, there were still some limitations in the present study. First, the sample size of certain disease studies might be insufficient to capture all of the critical disease-related genes for identifying the common DEGs. Second, incompleteness of available interactome data and limitation of computational methods might make the conclusions less dependable and accurate. Finally, the results of the present study were derived from multiple computational approaches, and future in vivo and in vitro experiments are required to fully assess the biological relevance of candidates.

Conclusion

In summary, based on bioinformatic analyses, we predicted promising therapeutic ingredients and drugs for NSCLC patients with COVID-19. We detected interactional hub genes enriched in regulating biological processes and signaling pathways, which were mainly relevant to cell division, cell aging, cell cycle, and cellular senescence. This study has identified 1) potential therapeutic targets, including CCNA2, CCNB2, AURKB, TOP2A, and BIRC5; 2) signaling pathways primarily related to cell cycle, cell aging, viral carcinogenesis, and p53 signaling pathway; 3) potential agents for the treatment of COVID-19 and NSCLC, including quercetin, resveratrol, cyclosporine, dasatinib, etoposide, ascorbic acid, and decitabine. Future experimental and clinical studies should be carried out with predicted agents to explore pharmacological mechanisms and to inform possible interventions for COVID-19 and NSCLC.

Statements

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

ZL and ZZ conceived and designed the study. ZZ and QC performed data analysis and data interpretation. ZZ, QC, and XZ conducted the bioinformatic and statistical analyses. ZZ and QC plotted the figures. QC, XZ, and HC prepared the manuscript.

Funding

This work was supported by the “Double First-Class” and High-level University Discipline Collaborative Innovation Team Project of Guangzhou University of Chinese Medicine (Grant No. 2021xk23).

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors, and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fphar.2022.857730/full#supplementary-material

Supplementary File S1

COVID-19–related genes obtained from public databases.

Supplementary File S2

COVID-19–related DEGs obtained from the RNA-seq dataset.

Supplementary File S3

NSCLC-related genes obtained from public databases.

Supplementary File S4

DEGs of NSCLC obtained from TCGA database.

Supplementary File S5

GO and KEGG terms based on the biggest sub-network with the associated statistics.

Supplementary File S6

TFs filtered from the ChEA3 database.

Supplementary File S7

MiRNAs obtained from Tarbase and mirTarbase.

Supplementary Table S1

Antiviral and/or anticancer activities of candidate drugs for COVID-19/NSCLC from Enrichr database and ShinyGO v0.75 database.

Glossary

  • AURKB

    Aurora kinase B

  • BIRC5

    Baculoviral IAP repeat-containing 5

  • BP

    Biological processes

  • BUB1B

    BUB1 mitotic checkpoint serine/threonine kinase B

  • CC

    Cellular components

  • CCNA2

    Cyclin A2

  • CCNB2

    Cyclin B2

  • CDC20

    Cell division cycle 20

  • CENPA

    Centromere protein A

  • ChEA3

    ChIP-X Enrichment Analysis 3

  • COVID-19

    Coronavirus disease-2019

  • CTD

    Comparative Toxicogenomics Database

  • DEGs

    Differentially expressed genes

  • DLGAP5

    DLG-associated protein 5

  • DNMT1

    DNA methyltransferase 1

  • DPFG

    Drug Perturbations from GEO

  • DSigDB

    Drug Signatures Database

  • FoxO

    Forkhead box, sub-group O

  • GEO

    Gene Expression Omnibus

  • GO

    Gene Ontology

  • IFN

    Interferon

  • IL

    Interleukin

  • JAK

    STAT Janus kinase/signal transducer and activator of transcription

  • KEGG

    Kyoto Encyclopedia of Genes and Genomes

  • MERS-CoV

    Middle East respiratory syndrome coronavirus

  • MF

    Molecular functions

  • MYBL2

    MYB proto-oncogene-like 2

  • NSCLC

    Non-small-cell lung cancer ACE2 angiotensin-converting enzyme 2

  • OMIM

    Online Mendelian Inheritance in Man

  • PLK1

    Polo-like kinase 1

  • PPI

    Protein–protein interaction

  • RNA-seq

    Transcriptomic RNA-sequencing

  • SARS-CoV-2

    Severe acute respiratory syndrome coronavirus 2

  • SASP

    Senescence-associated secretory phenotype

  • TCGA

    The Cancer Genome Atlas

  • TFDP1

    Transcription factor Dp-1

  • TFs

    Transcription factors

  • TOP2A

    DNA topoisomerase II alpha

  • TTD

    Therapeutic Target Database

  • ZNF367

    Zinc finger protein 367

References

  • 1

    AcostaJ. C.BanitoA.WuestefeldT.GeorgilisA.JanichP.MortonJ. P.et al (2013). A Complex Secretory Program Orchestrated by the Inflammasome Controls Paracrine Senescence. Nat. Cel. Biol.15 (8), 978990. 10.1038/ncb2784

  • 2

    AcostaJ. C.O'LoghlenA.BanitoA.GuijarroM. V.AugertA.RaguzS.et al (2008a). Chemokine Signaling via the CXCR2 Receptor Reinforces Senescence. Cell133 (6), 10061018. 10.1016/j.cell.2008.03.038

  • 3

    AcostaJ. C.O'LoghlenA.BanitoA.RaguzS.GilJ. (2008b). Control of Senescence by CXCR2 and its Ligands. Cell Cycle7 (19), 29562959. 10.4161/cc.7.19.6780

  • 4

    AfrasiabiE.AhlgrenJ.BergelinN.TörnquistK. (2008). Phorbol 12-Myristate 13-Acetate Inhibits FRO Anaplastic Human Thyroid Cancer Cell Proliferation by Inducing Cell Cycle Arrest in G1/S Phase: Evidence for an Effect Mediated by PKCdelta. Mol. Cel. Endocrinol.292 (1-2), 2635. 10.1016/j.mce.2008.04.018

  • 5

    AhmedF. (2020). A Network-Based Analysis Reveals the Mechanism Underlying Vitamin D in Suppressing Cytokine Storm and Virus in SARS-CoV-2 Infection. Front. Immunol.11, 590459. 10.3389/fimmu.2020.590459

  • 6

    AlnajeebiA. M.AlharbiH. F. H.AlelwaniW.BabteenN. A.AlansariW.ShamlanG.et al (2021). COVID-19 Candidate Genes and Pathways Potentially Share the Association with Lung Cancer. Cchts25, 110. 10.2174/1386207324666210712092649

  • 7

    ArgentaD. F.SilvaI. T.BassaniV. L.KoesterL. S.TeixeiraH. F.SimõesC. M. (2015). Antiherpes Evaluation of Soybean Isoflavonoids. Arch. Virol.160 (9), 23352342. 10.1007/s00705-015-2514-z

  • 8

    BaderG. D.HogueC. W. (2003). An Automated Method for Finding Molecular Complexes in Large Protein Interaction Networks. BMC Bioinformatics4, 2. 10.1186/1471-2105-4-2

  • 9

    BadiaR.AnguloG.Riveira-MuñozE.PujantellM.PuigT.RamirezC.et al (2016). Inhibition of Herpes Simplex Virus Type 1 by the CDK6 Inhibitor PD-0332991 (Palbociclib) through the Control of SAMHD1. J. Antimicrob. Chemother.71 (2), 387394. 10.1093/jac/dkv363

  • 10

    BeachL. B.RawsonJ. M.KimB.PattersonS. E.ManskyL. M. (2014). Novel Inhibitors of Human Immunodeficiency Virus Type 2 Infectivity. J. Gen. Virol.95 (Pt 12), 27782783. 10.1099/vir.0.069864-0

  • 11

    BergmannS.ElbaheshH. (2019). Targeting the Proviral Host Kinase, FAK, Limits Influenza a Virus Pathogenesis and NFkB-Regulated Pro-Inflammatory Responses. Virology534, 5463. 10.1016/j.virol.2019.05.020

  • 12

    BerlinI.ThomasD.Le FaouA. L.CornuzJ. (2020). COVID-19 and Smoking. Nicotine Tob. Res.22 (9), 16501652. 10.1093/ntr/ntaa059

  • 13

    Bertran-AlamilloJ.CattanV.SchoumacherM.Codony-ServatJ.Giménez-CapitánA.CanteroF.et al (2019). AURKB as a Target in Non-Small Cell Lung Cancer with Acquired Resistance to Anti-EGFR Therapy. Nat. Commun.10 (1), 1812. 10.1038/s41467-019-09734-5

  • 14

    BochenekM. L.SchützE.SchäferK. (2016). Endothelial Cell Senescence and Thrombosis: Ageing Clots. Thromb. Res.147, 3645. 10.1016/j.thromres.2016.09.019

  • 15

    BockJ. O.OrteaI. (2020). Re-analysis of SARS-CoV-2-Infected Host Cell Proteomics Time-Course Data by Impact Pathway Analysis and Network Analysis: A Potential Link with Inflammatory Response. Aging (Albany NY)12 (12), 1127711286. 10.18632/aging.103524

  • 16

    BorcherdingN.JethavaY.VikasP. (2020). Repurposing Anti-Cancer Drugs for COVID-19 Treatment. Drug Des. Devel. Ther.14, 50455058. 10.2147/dddt.S282252

  • 17

    BruunsgaardH.Andersen-RanbergK.HjelmborgJv.PedersenB. K.JeuneB. (2003). Elevated Levels of Tumor Necrosis Factor Alpha and Mortality in Centenarians. Am. J. Med.115 (4), 278283. 10.1016/s0002-9343(03)00329-2

  • 18

    CarewJ. S.EspitiaC. M.EsquivelJ. A.2ndMahalingamD.KellyK. R.ReddyG.et al (2011). Lucanthone Is a Novel Inhibitor of Autophagy that Induces Cathepsin D-Mediated Apoptosis. J. Biol. Chem.286 (8), 66026613. 10.1074/jbc.M110.151324

  • 19

    CarithersL. J.MooreH. M. (2015). The Genotype-Tissue Expression (GTEx) Project. Biopreserv. Biobank13 (5), 307308. 10.1089/bio.2015.29031.hmm

  • 20

    ChaiR.XuC.LuL.LiuX.MaZ. (2021). Quercetin Inhibits Proliferation of and Induces Apoptosis in Non-Small-Cell Lung Carcinoma via the lncRNA SNHG7/miR-34a-5p Pathway. Immunopharmacol. Immunotoxicol.43 (6), 693703. 10.1080/08923973.2021.1966032

  • 21

    ChanK. K. L.SiuM. K. Y.JiangY. X.WangJ. J.LeungT. H. Y.NganH. Y. S. (2018). Estrogen Receptor Modulators Genistein, Daidzein and ERB-041 Inhibit Cell Migration, Invasion, Proliferation and Sphere Formation via Modulation of FAK and PI3K/AKT Signaling in Ovarian Cancer. Cancer Cel. Int.18, 65. 10.1186/s12935-018-0559-2

  • 22

    ChenJ.-C.XieT.-A.LinZ.-Z.LiY.-Q.XieY.-F.LiZ.-W.et al (2021). Identification of Key Pathways and Genes in SARS-CoV-2 Infecting Human Intestines by Bioinformatics Analysis. Biochem. Genet.119. 10.1007/s10528-021-10144-w

  • 23

    ChènebyJ.MénétrierZ.MestdaghM.RosnetT.DouidaA.RhalloussiW.et al (2020). ReMap 2020: A Database of Regulatory Regions from an Integrative Analysis of Human and Arabidopsis DNA-Binding Sequencing Experiments. Nucleic Acids Res.48 (D1), D180d188. 10.1093/nar/gkz945

  • 24

    ChengH.ZouY.ShahC. D.FanN.BhagatT. D.GucalpR.et al (2021). First-in-Human Study of Inhaled Azacitidine in Patients with Advanced Non-Small Cell Lung Cancer. Lung Cancer154, 99104. 10.1016/j.lungcan.2021.02.015

  • 25

    ChinC. H.ChenS. H.WuH. H.HoC. W.KoM. T.LinC. Y. (2014). cytoHubba: Identifying Hub Objects and Sub-Networks from Complex Interactome. BMC Syst. Biol.8 (Suppl. 4), S11. 10.1186/1752-0509-8-s4-s11

  • 26

    ChouC. H.ShresthaS.YangC. D.ChangN. W.LinY. L.LiaoK. W.et al (2018). miRTarBase Update 2018: A Resource for Experimentally Validated microRNA-Target Interactions. Nucleic Acids Res.46 (D1), D296d302. 10.1093/nar/gkx1067

  • 27

    ChuB. F.KarpenkoM. J.LiuZ.AimiuwuJ.Villalona-CaleroM. A.ChanK. K.et al (2013). Phase I Study of 5-Aza-2'-Deoxycytidine in Combination with Valproic Acid in Non-Small-Cell Lung Cancer. Cancer Chemother. Pharmacol.71 (1), 115121. 10.1007/s00280-012-1986-8

  • 28

    CiottiM.CiccozziM.TerrinoniA.JiangW. C.WangC. B.BernardiniS. (2020). The COVID-19 Pandemic. Crit. Rev. Clin. Lab. Sci.57 (6), 365388. 10.1080/10408363.2020.1783198

  • 29

    CloughE.BarrettT. (2016). “The Gene Expression Omnibus Database,” in Statistical Genomics (New York, NY: Springer), 93110. 10.1007/978-1-4939-3578-9_5

  • 30

    Colunga BiancatelliR. M. L.BerrillM.MarikP. E. (2020). The Antiviral Properties of Vitamin C. Expert Rev. Anti Infect. Ther.18 (2), 99101. 10.1080/14787210.2020.1706483

  • 31

    CooperW. A.Kohonen-CorishM. R.McCaughanB.KennedyC.SutherlandR. L.LeeC. S. (2009). Expression and Prognostic Significance of Cyclin B1 and Cyclin A in Non-Small Cell Lung Cancer. Histopathology55 (1), 2836. 10.1111/j.1365-2559.2009.03331.x

  • 32

    CoppéJ. P.DesprezP. Y.KrtolicaA.CampisiJ. (2010). The Senescence-Associated Secretory Phenotype: The Dark Side of Tumor Suppression. Annu. Rev. Pathol.5, 99118. 10.1146/annurev-pathol-121808-102144

  • 33

    DaamenA. R.BachaliP.OwenK. A.KingsmoreK. M.HubbardE. L.LabonteA. C.et al (2021). Comprehensive Transcriptomic Analysis of COVID-19 Blood, Lung, and Airway. Sci. Rep.11 (1), 7052. 10.1038/s41598-021-86002-x

  • 34

    DavisA. P.GrondinC. J.JohnsonR. J.SciakyD.WiegersJ.WiegersT. C.et al (2021). Comparative Toxicogenomics Database (CTD): Update 2021. Nucleic Acids Res.49 (D1), D1138d1143. 10.1093/nar/gkaa891

  • 35

    DavisC. A.HitzB. C.SloanC. A.ChanE. T.DavidsonJ. M.GabdankI.et al (2018). The Encyclopedia of DNA Elements (ENCODE): Data portal Update. Nucleic Acids Res.46 (D1), D794d801. 10.1093/nar/gkx1081

  • 36

    De BiasiS.MeschiariM.GibelliniL.BellinazziC.BorellaR.FidanzaL.et al (2020). Marked T Cell Activation, Senescence, Exhaustion and Skewing towards TH17 in Patients with COVID-19 Pneumonia. Nat. Commun.11 (1), 3434. 10.1038/s41467-020-17292-4

  • 37

    de Carvalho Scharf SantanaN.LimaN. A.DesotiV. C.BidóiaD. L.de Souza Bonfim MendonçaP.RattiB. A.et al (2016). Vitamin K3 Induces Antiproliferative Effect in Cervical Epithelial Cells Transformed by HPV 16 (SiHa Cells) through the Increase in Reactive Oxygen Species Production. Arch. Gynecol. Obstet.294 (4), 797804. 10.1007/s00404-016-4097-7

  • 38

    de la LastraC. A.VillegasI. (2007). Resveratrol as an Antioxidant and Pro-Oxidant Agent: Mechanisms and Clinical Implications. Biochem. Soc. Trans.35 (Pt 5), 11561160. 10.1042/bst0351156

  • 39

    DerosaG.MaffioliP.D'AngeloA.Di PierroF. (2021). A Role for Quercetin in Coronavirus Disease 2019 (COVID-19). Phytother. Res.35 (3), 12301236. 10.1002/ptr.6887

  • 40

    DI PierroF.DerosaG.MaffioliP.BertuccioliA.TogniS.RivaA.et al (2021a). Possible Therapeutic Effects of Adjuvant Quercetin Supplementation Against Early-Stage COVID-19 Infection: A Prospective, Randomized, Controlled, and Open-Label Study. Int. J. Gen. Med.14, 23592366. 10.2147/ijgm.S318720

  • 41

    DI PierroF.KhanA.BertuccioliA.MaffioliP.DerosaG.KhanB. A.et al (2021b). Quercetin Phytosome® as a Potential Candidate for Managing COVID-19. Minerva Gastroenterol. (Torino)67 (2), 190195. 10.23736/s2724-5985.20.02771-3

  • 42

    Diaz-MoralliS.Tarrado-CastellarnauM.MirandaA.CascanteM. (2013). Targeting Cell Cycle Regulation in Cancer Therapy. Pharmacol. Ther.138 (2), 255271. 10.1016/j.pharmthera.2013.01.011

  • 43

    DidelotC.MirjoletJ. F.Barberi-HeyobM.RamacciC.TeitenM. H.MerlinJ. L. (2003). Oncoprotein Expression of E6 and E7 Does Not Prevent 5-Fluorouracil (5FU) Mediated G1/S Arrest and Apoptosis in 5FU Resistant Carcinoma Cell Lines. Int. J. Oncol.23 (1), 8187. 10.3892/ijo.23.1.81

  • 44

    DuJ.CullenJ. J.BuettnerG. R. (2012). Ascorbic Acid: Chemistry, Biology and the Treatment of Cancer. Biochim. Biophys. Acta1826 (2), 443457. 10.1016/j.bbcan.2012.06.003

  • 45

    DurmuşS.ÇakırT.ÖzgürA.GuthkeR. (2015). A Review on Computational Systems Biology of Pathogen-Host Interactions. Front. Microbiol.6, 235. 10.3389/fmicb.2015.00235

  • 46

    DyallJ.ColemanC. M.HartB. J.VenkataramanT.HolbrookM. R.KindrachukJ.et al (2014). Repurposing of Clinically Developed Drugs for Treatment of Middle East Respiratory Syndrome Coronavirus Infection. Antimicrob. Agents Chemother.58 (8), 48854893. 10.1128/AAC.03036-14

  • 47

    DzoboK. (2021). The Role of Viruses in Carcinogenesis and Molecular Targeting: From Infection to Being a Component of the Tumor Microenvironment. OMICS: A J. Integr. Biol.25 (6), 358371. 10.1089/omi.2021.0052

  • 48

    EastmanR. T.RothJ. S.BrimacombeK. R.SimeonovA.ShenM.PatnaikS.et al (2020). Remdesivir: A Review of its Discovery and Development Leading to Emergency Use Authorization for Treatment of COVID-19. ACS Cent. Sci.6 (5), 672683. 10.1021/acscentsci.0c00489

  • 49

    EscherR.BreakeyN.LämmleB. (2020). Severe COVID-19 Infection Associated with Endothelial Activation. Thromb. Res.190, 62. 10.1016/j.thromres.2020.04.014

  • 50

    EssebierA.LamprechtM.PiperM.BodénM. (2017). Bioinformatics Approaches to Predict Target Genes from Transcription Factor Binding Data. Methods131, 111119. 10.1016/j.ymeth.2017.09.001

  • 51

    FarhanM.WangH.GaurU.LittleP. J.XuJ.ZhengW. (2017). FOXO Signaling Pathways as Therapeutic Targets in Cancer. Int. J. Biol. Sci.13 (7), 815827. 10.7150/ijbs.20052

  • 52

    FinkelT.SerranoM.BlascoM. A. (2007). The Common Biology of Cancer and Ageing. Nature448 (7155), 767774. 10.1038/nature05985

  • 53

    FrascaD.DiazA.RomeroM.BlombergB. B. (2017). Human Peripheral Late/Exhausted Memory B Cells Express a Senescent-Associated Secretory Phenotype and Preferentially Utilize Metabolic Signaling Pathways. Exp. Gerontol.87 (Pt A), 113120. 10.1016/j.exger.2016.12.001

  • 54

    FujitaM.HasegawaA.YamamoriM.OkamuraN. (2017). In Vitro and In Vivo Cytotoxicity of Troglitazone in Pancreatic Cancer. J. Exp. Clin. Cancer Res.36 (1), 91. 10.1186/s13046-017-0557-6

  • 55

    FukanoK.TsukudaS.OshimaM.SuzukiR.AizakiH.OhkiM.et al (2018). Troglitazone Impedes the Oligomerization of Sodium Taurocholate Cotransporting Polypeptide and Entry of Hepatitis B Virus into Hepatocytes. Front. Microbiol.9, 3257. 10.3389/fmicb.2018.03257

  • 56

    GallardoE.NavarroA.ViñolasN.MarradesR. M.DiazT.GelB.et al (2009). miR-34a as a Prognostic Marker of Relapse in Surgically Resected Non-Small-Cell Lung Cancer. Carcinogenesis30 (11), 19031909. 10.1093/carcin/bgp219

  • 57

    GaoN.WangF. X.WangG.ZhaoQ. S. (2018). Targeting the HMGA2 Oncogene by miR-498 Inhibits Non-Small Cell Lung Cancer Biological Behaviors. Eur. Rev. Med. Pharmacol. Sci.22 (6), 16931699. 10.26355/eurrev_201803_14582

  • 58

    GarbeC.EigentlerT. K. (2018). Vemurafenib. Recent Results Cancer Res.211, 7789. 10.1007/978-3-319-91442-8_6

  • 59

    GeS. X.JungD.YaoR. (2020). ShinyGO: A Graphical Gene-Set Enrichment Tool for Animals and Plants. Bioinformatics36 (8), 26282629. 10.1093/bioinformatics/btz931

  • 60

    GencerE. B.UralA. U.AvcuF.BaranY. (2011). A Novel Mechanism of Dasatinib-Induced Apoptosis in Chronic Myeloid Leukemia; Ceramide Synthase and Ceramide Clearance Genes. Ann. Hematol.90 (11), 12651275. 10.1007/s00277-011-1212-5

  • 61

    GhoshZ.MallickB.ChakrabartiJ. (2009). Cellular versus Viral microRNAs in Host-Virus Interaction. Nucleic Acids Res.37 (4), 10351048. 10.1093/nar/gkn1004

  • 62

    GoldmanM. J.CraftB.HastieM.RepečkaK.McDadeF.KamathA.et al (2020). Visualizing and Interpreting Cancer Genomics Data via the Xena Platform. Nat. Biotechnol.38 (6), 675678. 10.1038/s41587-020-0546-8

  • 63

    GomesM.TeixeiraA. L.CoelhoA.AraújoA.MedeirosR. (2014). The Role of Inflammation in Lung Cancer. Adv. Exp. Med. Biol.816, 123. 10.1007/978-3-0348-0837-8_1

  • 64

    González‐GualdaE.BakerA. G.FrukL.Muñoz‐EspínD. (2021). A Guide to Assessing Cellular Senescence In Vitro and In Vivo. Febs j288 (1), 5680. 10.1111/febs.15570

  • 65

    Guisado-VascoP.Valderas-OrtegaS.Carralón-GonzálezM. M.Roda-SantacruzA.González-CortijoL.Sotres-FernándezG.et al (2020). Clinical Characteristics and Outcomes Among Hospitalized Adults with Severe COVID-19 Admitted to a Tertiary Medical Center and Receiving Antiviral, Antimalarials, Glucocorticoids, or Immunomodulation with Tocilizumab or Cyclosporine: A Retrospective Observational Study (COQUIMA Cohort). EClinicalMedicine28, 100591. 10.1016/j.eclinm.2020.100591

  • 66

    GuoQ.ZhongM.XuH.MaoX.ZhangY.LinN. (2015). A Systems Biology Perspective on the Molecular Mechanisms Underlying the Therapeutic Effects of Buyang Huanwu Decoction on Ischemic Stroke. Rejuvenation Res.18 (4), 313325. 10.1089/rej.2014.1635

  • 67

    GuptaS.SilveiraD. A.MombachJ. C. M. (2020). ATM/miR-34a-5p Axis Regulates a P21-dependent Senescence-Apoptosis Switch in Non-small Cell Lung Cancer: A Boolean Model of G1/S Checkpoint Regulation. FEBS Lett.594 (2), 227239. 10.1002/1873-3468.13615

  • 68

    HadjadjJ.YatimN.BarnabeiL.CorneauA.BoussierJ.SmithN.et al (2020). Impaired Type I Interferon Activity and Inflammatory Responses in Severe COVID-19 Patients. Science369 (6504), 718724. 10.1126/science.abc6027

  • 69

    HamamotoI.HarazakiK.InaseN.TakakuH.TashiroM.YamamotoN. (2013). Cyclosporin A Inhibits the Propagation of Influenza Virus by Interfering with a Late Event in the Virus Life Cycle. Jpn. J. Infect. Dis.66 (4), 276283. 10.7883/yoken.66.276

  • 70

    HamoshA.ScottA. F.AmbergerJ. S.BocchiniC. A.McKusickV. A. (2005). Online Mendelian Inheritance in Man (OMIM), a Knowledgebase of Human Genes and Genetic Disorders. Nucleic Acids Res.33 (Database Issue), D514D517. 10.1093/nar/gki033

  • 71

    HanN.YuanX.WuH.XuH.ChuQ.GuoM.et al (2015). DACH1 Inhibits Lung Adenocarcinoma Invasion and Tumor Growth by Repressing CXCL5 Signaling. Oncotarget6 (8), 58775888. 10.18632/oncotarget.3463

  • 72

    HanX.HuangC.QuX.LiuS.YangX.WangY.et al (2019). Tripartite Motif-Containing 15 Overexpression in Non-Small Cell Lung Cancer Is Associated with Poor Patient Prognoses. J. Cancer10 (4), 843852. 10.7150/jca.27856

  • 73

    HanahanD.WeinbergR. A. (2000). The Hallmarks of Cancer. Cell100 (1), 5770. 10.1016/s0092-8674(00)81683-9

  • 74

    HasanM. I.RahmanM. H.IslamM. B.IslamM. Z.HossainM. A.MoniM. A. (2022). Systems Biology and Bioinformatics Approach to Identify Blood Based Signatures Molecules and Drug Targets of Patient with COVID-19. Inform. Med. Unlocked28, 100840. 10.1016/j.imu.2021.100840

  • 75

    HasegawaT.SuzukiH.AbeJ.SakuradaA.EndoC.SatoN.et al (2020). A Phase 2 Study of Adjuvant Carboplatin Plus S-1 Followed by Maintenance S-1 Therapy for Patients with Completely Resected Stage II/IIIA Non-small Cell Lung Cancer-Japanese Northern East Area Thoracic Surgery Study Group JNETS1302 Study. J. Thorac. Dis.12 (7), 35913601. 10.21037/jtd-20-715

  • 76

    Hernandez-SeguraA.NehmeJ.DemariaM. (2018). Hallmarks of Cellular Senescence. Trends Cel. Biol.28 (6), 436453. 10.1016/j.tcb.2018.02.001

  • 77

    HoG. Y.WoodwardN.CowardJ. I. (2016). Cisplatin versus Carboplatin: Comparative Review of Therapeutic Management in Solid Malignancies. Crit. Rev. Oncol. Hematol.102, 3746. 10.1016/j.critrevonc.2016.03.014

  • 78

    HolzbergM.BoergelingY.SchräderT.LudwigS.EhrhardtC. (2017). Vemurafenib Limits Influenza A Virus Propagation by Targeting Multiple Signaling Pathways. Front. Microbiol.8, 2426. 10.3389/fmicb.2017.02426

  • 79

    HouG. X.LiuP.YangJ.WenS. (2017). Mining Expression and Prognosis of Topoisomerase Isoforms in Non-Small-Cell Lung Cancer by Using Oncomine and Kaplan-Meier Plotter. PLoS One12 (3), e0174515. 10.1371/journal.pone.0174515

  • 80

    HuH.WangK.WangL.DuY.ChenJ.LiY.et al (2021). He-Jie-Shen-Shi Decoction as an Adjuvant Therapy on Severe Coronavirus Disease 2019: A Retrospective Cohort and Potential Mechanistic Study. Front. Pharmacol.12, 700498. 10.3389/fphar.2021.700498

  • 81

    HuX.LiJ.FuM.ZhaoX.WangW. (2021). The JAK/STAT Signaling Pathway: From Bench to Clinic. Signal. Transduct. Target. Ther.6 (1), 402. 10.1038/s41392-021-00791-1

  • 82

    HuangZ.CaoY.ZhouM.QiX.FuB.MouY.et al (2019). Hsa_circ_0005519 Increases IL-13/IL-6 by Regulating Hsa-Let-7a-5p in CD4+ T Cells to Affect Asthma. Clin. Exp. Allergy49 (8), 11161127. 10.1111/cea.13445

  • 83

    HungI. F.LungK. C.TsoE. Y.LiuR.ChungT. W.ChuM. Y.et al (2020). Triple Combination of Interferon Beta-1b, Lopinavir-Ritonavir, and Ribavirin in the Treatment of Patients Admitted to Hospital with COVID-19: An Open-Label, Randomised, Phase 2 Trial. Lancet395 (10238), 16951704. 10.1016/s0140-6736(20)31042-4

  • 84

    JossetL.TextorisJ.LoriodB.FerrarisO.MoulesV.LinaB.et al (2010). Gene Expression Signature-Based Screening Identifies New Broadly Effective Influenza a Antivirals. PLoS One5 (10), e13169. 10.1371/journal.pone.0013169

  • 85

    Kandhaya-PillaiR.Miro-MurF.Alijotas-ReigJ.TchkoniaT.KirklandJ. L.SchwartzS. (2017). Tnfα-Senescence Initiates a STAT-Dependent Positive Feedback Loop, Leading to a Sustained Interferon Signature, DNA Damage, and Cytokine Secretion. Aging (Albany NY)9 (11), 24112435. 10.18632/aging.101328

  • 86

    KanehisaM.GotoS. (2000). KEGG: Kyoto Encyclopedia of Genes and Genomes. Nucleic Acids Res.28 (1), 2730. 10.1093/nar/28.1.27

  • 87

    KaragkouniD.ParaskevopoulouM. D.ChatzopoulosS.VlachosI. S.TastsoglouS.KanellosI.et al (2018). DIANA-TarBase V8: A Decade-Long Collection of Experimentally Supported miRNA-Gene Interactions. Nucleic Acids Res.46 (D1), D239d245. 10.1093/nar/gkx1141

  • 88

    KasinskiA. L.SlackF. J. (2012). miRNA-34 Prevents Cancer Initiation and Progression in a Therapeutically Resistant K-Ras and P53-Induced Mouse Model of Lung Adenocarcinoma. Cancer Res.72 (21), 55765587. 10.1158/0008-5472.Can-12-2001

  • 89

    Kedhari SundaramM.HussainA.HaqueS.RainaR.AfrozeN. (2019). Quercetin Modifies 5'CpG Promoter Methylation and Reactivates Various Tumor Suppressor Genes by Modulating Epigenetic Marks in Human Cervical Cancer Cells. J. Cel. Biochem.120 (10), 1835718369. 10.1002/jcb.29147

  • 90

    KeenanA. B.TorreD.LachmannA.LeongA. K.WojciechowiczM. L.UttiV.et al (2019). ChEA3: Transcription Factor Enrichment Analysis by Orthogonal Omics Integration. Nucleic Acids Res.47 (W1), W212w224. 10.1093/nar/gkz446

  • 91

    KhokharM.TomoS.PurohitP. (2022). MicroRNAs Based Regulation of Cytokine Regulating Immune Expressed Genes and Their Transcription Factors in COVID-19. Meta Gene31, 100990. 10.1016/j.mgene.2021.100990

  • 92

    KimD. H.KhanH.UllahH.HassanS. T. S.ŠmejkalK.EfferthT.et al (2019). MicroRNA Targeting by Quercetin in Cancer Treatment and Chemoprotection. Pharmacol. Res.147, 104346. 10.1016/j.phrs.2019.104346

  • 93

    KimK. S.KangK. W.SeuY. B.BaekS. H.KimJ. R. (2009). Interferon-Gamma Induces Cellular Senescence Through P53-Dependent DNA Damage Signaling in Human Endothelial Cells. Mech. Ageing Dev.130 (3), 179188. 10.1016/j.mad.2008.11.004

  • 94

    KimS.ChenJ.ChengT.GindulyteA.HeJ.HeS.et al (2019). PubChem 2019 Update: Improved Access to Chemical Data. Nucleic Acids Res.47 (D1), D1102d1109. 10.1093/nar/gky1033

  • 95

    KimW. R.ParkE. G.KangK. W.LeeS. M.KimB.KimH. S. (2020). Expression Analyses of MicroRNAs in Hamster Lung Tissues Infected by SARS-CoV-2. Mol. Cell43 (11), 953963. 10.14348/molcells.2020.0177

  • 96

    KimuraS.KuramotoK.HomanJ.NaruokaH.EgoT.NogawaM.et al (2012). Antiproliferative and Antitumor Effects of Azacitidine against the Human Myelodysplastic Syndrome Cell Line SKM-1. Anticancer Res.32 (3), 795798.

  • 97

    KishoreC.SundaramS.KarunagaranD. (2019). Vitamin K3 (Menadione) Suppresses Epithelial-Mesenchymal-Transition and Wnt Signaling Pathway in Human Colorectal Cancer Cells. Chem. Biol. Interact309, 108725. 10.1016/j.cbi.2019.108725

  • 98

    KuleshovM. V.JonesM. R.RouillardA. D.FernandezN. F.DuanQ.WangZ.et al (2016). Enrichr: a Comprehensive Gene Set Enrichment Analysis Web Server 2016 Update. Nucleic Acids Res.44 (W1), W90W97. 10.1093/nar/gkw377

  • 99

    KumariR.JatP. (2021). Mechanisms of Cellular Senescence: Cell Cycle Arrest and Senescence Associated Secretory Phenotype. Front. Cel. Dev. Biol.9, 645593. 10.3389/fcell.2021.645593

  • 100

    Kuźnar-KamińskaB.Mikuła-PietrasikJ.KsiążekK.TykarskiA.Batura-GabryelH. (2018). Lung Cancer in Chronic Obstructive Pulmonary Disease Patients: Importance of Cellular Senescence. Pol. Arch. Intern. Med.128 (7-8), 462468. 10.20452/pamw.4297

  • 101

    Kuźnar-KamińskaB.Mikuła-PietrasikJ.SosińskaP.KsiążekK.Batura-GabryelH. (2016). COPD Promotes Migration of A549 Lung Cancer Cells: the Role of Chemokine CCL21. Int. J. Chron. Obstruct. Pulmon. Dis.11, 10611066. 10.2147/copd.S96490

  • 102

    LachmannA.TorreD.KeenanA. B.JagodnikK. M.LeeH. J.WangL.et al (2018). Massive Mining of Publicly Available RNA-Seq Data from Human and Mouse. Nat. Commun.9 (1), 1366. 10.1038/s41467-018-03751-6

  • 103

    LeeB.MinJ. A.NashedA.LeeS. O.YooJ. C.ChiS. W.et al (2019). A Novel Mechanism of Irinotecan Targeting MDM2 and Bcl-xL. Biochem. Biophys. Res. Commun.514 (2), 518523. 10.1016/j.bbrc.2019.04.009

  • 104

    LeeH.YoonN. E.JungB. H. (2017). Metabolomics Study of Cancer Targeting Small Molecule Kinase Inhibitors in Cell Culture. Drug Metab. Pharmacokinet.32 (1), S77. 10.1016/j.dmpk.2016.10.304

  • 105

    LeidingerP.BrefortT.BackesC.KrappM.GalataV.BeierM.et al (2016). High-throughput qRT-PCR Validation of Blood microRNAs in Non-small Cell Lung Cancer. Oncotarget7 (4), 46114623. 10.18632/oncotarget.6566

  • 106

    LiC.HuX.LiL.LiJ. H. (2020). Differential microRNA Expression in the Peripheral Blood from Human Patients with COVID-19. J. Clin. Lab. Anal.34 (10), e23590. 10.1002/jcla.23590

  • 107

    LiJ.ZhengS.ChenB.ButteA. J.SwamidassS. J.LuZ. (2016). A Survey of Current Trends in Computational Drug Repositioning. Brief Bioinform17 (1), 212. 10.1093/bib/bbv020

  • 108

    LiQ.LiM.ZhengK.LiH.YangH.MaS.et al (2020). Detection of microRNA Expression Levels Based on Microarray Analysis for Classification of Idiopathic Pulmonary Fibrosis. Exp. Ther. Med.20 (4), 30963103. 10.3892/etm.2020.9068

  • 109

    LiX.MeiQ.NieJ.FuX.HanW. (2015). Decitabine: a Promising Epi-Immunotherapeutic Agent in Solid Tumors. Expert Rev. Clin. Immunol.11 (3), 363375. 10.1586/1744666x.2015.1002397

  • 110

    LiY.YaoJ.HanC.YangJ.ChaudhryM. T.WangS.et al (2016). Quercetin, Inflammation and Immunity. Nutrients8 (3), 167. 10.3390/nu8030167

  • 111

    LiangJ.BiN.WuS.ChenM.LvC.ZhaoL.et al (2017). Etoposide and Cisplatin versus Paclitaxel and Carboplatin with Concurrent Thoracic Radiotherapy in Unresectable Stage III Non-small Cell Lung Cancer: a Multicenter Randomized Phase III Trial. Ann. Oncol.28 (4), 777783. 10.1093/annonc/mdx009

  • 112

    LiangW.GuanW.ChenR.WangW.LiJ.XuK.et al (2020). Cancer Patients in SARS-CoV-2 Infection: a Nationwide Analysis in China. Lancet Oncol.21 (3), 335337. 10.1016/s1470-2045(20)30096-6

  • 113

    LiaoZ.Rivin Del CampoE.SalemA.PangQ.LiuH.Lopez GuerraJ. L. (2020). Optimizing Lung Cancer Radiation Treatment Worldwide in COVID-19 Outbreak. Lung Cancer146, 230235. 10.1016/j.lungcan.2020.05.029

  • 114

    LimW.JeongM.BazerF. W.SongG. (2017). Coumestrol Inhibits Proliferation and Migration of Prostate Cancer Cells by Regulating AKT, ERK1/2, and JNK MAPK Cell Signaling Cascades. J. Cel. Physiol.232 (4), 862871. 10.1002/jcp.25494

  • 115

    LinS. C.HoC. T.ChuoW. H.LiS.WangT. T.LinC. C. (2017). Effective Inhibition of MERS-CoV Infection by Resveratrol. BMC Infect. Dis.17 (1), 144. 10.1186/s12879-017-2253-8

  • 116

    LiuB.LiJ.CairnsM. J. (2014). Identifying miRNAs, Targets and Functions. Brief Bioinform.15 (1), 119. 10.1093/bib/bbs075

  • 117

    LiuM.LiuH.ChenJ. (2018). Mechanisms of the CDK4/6 Inhibitor Palbociclib (PD 0332991) and its Future Application in Cancer Treatment (Review). Oncol. Rep.39 (3), 901911. 10.3892/or.2018.6221

  • 118

    LiuQ.YuZ.YuanS.XieW.LiC.HuZ.et al (2017). Circulating Exosomal microRNAs as Prognostic Biomarkers for Non-small-cell Lung Cancer. Oncotarget8 (8), 1304813058. 10.18632/oncotarget.14369

  • 119

    LiuY.WangK.MassoudT. F.PaulmuruganR. (2020). SARS-CoV-2 Vaccine Development: An Overview and Perspectives. ACS Pharmacol. Transl. Sci.3 (5), 844858. 10.1021/acsptsci.0c00109

  • 120

    LiuZ.JiangL.ZhangG.LiS.JiangX. (2018). MiR-24 Promotes Migration and Invasion of Non-small Cell Lung Cancer by Targeting ZNF367. J. Buon.23 (5), 14131419.

  • 121

    LongleyD. B.HarkinD. P.JohnstonP. G. (2003). 5-fluorouracil: Mechanisms of Action and Clinical Strategies. Nat. Rev. Cancer3 (5), 330338. 10.1038/nrc1074

  • 122

    Lopez BernalJ.AndrewsN.GowerC.GallagherE.SimmonsR.ThelwallS.et al (2021). Effectiveness of Covid-19 Vaccines against the B.1.617.2 (Delta) Variant. N. Engl. J. Med.385 (7), 585594. 10.1056/NEJMoa2108891

  • 123

    LoveM. I.HuberW.AndersS. (2014). Moderated Estimation of Fold Change and Dispersion for RNA-Seq Data with DESeq2. Genome Biol.15 (12), 550. 10.1186/s13059-014-0550-8

  • 124

    MaW.WangB.ZhangY.WangZ.NiuD.ChenS.et al (2019). Prognostic Significance of TOP2A in Non-small Cell Lung Cancer Revealed by Bioinformatic Analysis. Cancer Cel. Int.19, 239. 10.1186/s12935-019-0956-1

  • 125

    MajumderJ.MinkoT. (2021). Recent Developments on Therapeutic and Diagnostic Approaches for COVID-19. Aaps j23 (1), 14. 10.1208/s12248-020-00532-2

  • 126

    MalkaniN.RashidM. U. (2021). SARS-COV-2 Infection and Lung Tumor Microenvironment. Mol. Biol. Rep.48 (2), 19251934. 10.1007/s11033-021-06149-8

  • 127

    MallahS. I.GhorabO. K.Al-SalmiS.AbdellatifO. S.TharmaratnamT.IskandarM. A.et al (2021). COVID-19: Breaking Down a Global Health Crisis. Ann. Clin. Microbiol. Antimicrob.20 (1), 35. 10.1186/s12941-021-00438-7

  • 128

    MasonR. J. (2020). Pathogenesis of COVID-19 from a Cell Biology Perspective. Eur. Respir. J.55 (4), 2000607. 10.1183/13993003.00607-2020

  • 129

    MercatelliD.PedaceE.VeltriP.GiorgiF. M.GuzziP. H. (2021). Exploiting the Molecular Basis of Age and Gender Differences in Outcomes of SARS-CoV-2 Infections. Comput. Struct. Biotechnol. J.19, 40924100. 10.1016/j.csbj.2021.07.002

  • 130

    MontecuccoA.ZanettaF.BiamontiG. (2015). Molecular Mechanisms of Etoposide. Excli j14, 95108. 10.17179/excli2015-561

  • 131

    MuhammadJ. S.Saheb Sharif-AskariN.CuiZ. G.HamadM.HalwaniR. (2021). SARS-CoV-2 Infection-Induced Promoter Hypomethylation as an Epigenetic Modulator of Heat Shock Protein A1L (HSPA1L) Gene. Front. Genet.12, 622271. 10.3389/fgene.2021.622271

  • 132

    MullenD. J.YanC.KangD. S.ZhouB.BorokZ.MarconettC. N.et al (2020). TENET 2.0: Identification of Key Transcriptional Regulators and Enhancers in Lung Adenocarcinoma. Plos Genet.16 (9), e1009023. 10.1371/journal.pgen.1009023

  • 133

    MurakamiS.SaitoH.KondoT.ItoH.OshitaF.YamadaK. (2018). Phase II Study of Nedaplatin and Irinotecan as Adjuvant Chemotherapy for Completely Resected Non-small Cell Lung Cancer. Cancer Chemother. Pharmacol.81 (1), 8187. 10.1007/s00280-017-3460-0

  • 134

    NamH. J.van DeursenJ. M. (2014). Cyclin B2 and P53 Control Proper Timing of Centrosome Separation. Nat. Cel. Biol.16 (6), 538549. 10.1038/ncb2952

  • 135

    NehmeJ.BorghesanM.MackedenskiS.BirdT. G.DemariaM. (2020). Cellular Senescence as a Potential Mediator of COVID-19 Severity in the Elderly. Aging Cell19 (10), e13237. 10.1111/acel.13237

  • 136

    NicośM.Rolska-KopińskaA.KrawczykP.GrendaA.BożykA.SzczyrekM.et al (2021). Effect of TOP2A and ERCC1 Gene Polymorphisms on the Efficacy and Toxicity of Cisplatin and Etoposide-Based Chemotherapy in Small Cell Lung Cancer Patients. Arch. Med. Sci.17 (2), 474480. 10.5114/aoms.2020.92572

  • 137

    NieH.ZhouX.ShuzhangD.NieC.ZhangX.HuangJ. (2019). Palbociclib Overcomes Afatinib Resistance in Non-small Cell Lung Cancer. Biomed. Pharmacother.109, 17501757. 10.1016/j.biopha.2018.10.170

  • 138

    NokiharaH.LuS.MokT. S. K.NakagawaK.YamamotoN.ShiY. K.et al (2017). Randomized Controlled Trial of S-1 versus Docetaxel in Patients with Non-small-cell Lung Cancer Previously Treated with Platinum-Based Chemotherapy (East Asia S-1 Trial in Lung Cancer). Ann. Oncol.28 (11), 26982706. 10.1093/annonc/mdx419

  • 139

    OlasB.WachowiczB. (2005). Resveratrol, a Phenolic Antioxidant with Effects on Blood Platelet Functions. Platelets16 (5), 251260. 10.1080/09537100400020591

  • 140

    OlloquequiJ. (2020). COVID-19 Susceptibility in Chronic Obstructive Pulmonary Disease. Eur. J. Clin. Invest.50 (10), e13382. 10.1111/eci.13382

  • 141

    OvermyerK. A.ShishkovaE.MillerI. J.BalnisJ.BernsteinM. N.Peters-ClarkeT. M.et al (2021). Large-Scale Multi-Omic Analysis of COVID-19 Severity. Cell Syst.12 (1), 23e7. e27. 10.1016/j.cels.2020.10.003

  • 142

    ParsonsP. G.LeanJ.KhooS. K.LarkJ. (1989). Effects of Adriamycin and Etoposide on the Replication of Adenovirus 5 in Sensitive and Resistant Human Tumour Cells. Biochem. Pharmacol.38 (1), 3137. 10.1016/0006-2952(89)90145-7

  • 143

    PassaroA.BestvinaC.Velez VelezM.GarassinoM. C.GaronE.PetersS. (2021). Severity of COVID-19 in Patients with Lung Cancer: Evidence and Challenges. J. Immunother. Cancer9 (3), e002266. 10.1136/jitc-2020-002266

  • 144

    PassaroA.PetersS.MokT. S. K.AttiliI.MitsudomiT.de MarinisF. (2020). Testing for COVID-19 in Lung Cancer Patients. Ann. Oncol.31 (7), 832834. 10.1016/j.annonc.2020.04.002

  • 145

    PathakA. K.BhutaniM.GuleriaR.BalS.MohanA.MohantiB. K.et al (2005). Chemotherapy Alone vs. Chemotherapy Plus High Dose Multiple Antioxidants in Patients with Advanced Non Small Cell Lung Cancer. J. Am. Coll. Nutr.24 (1), 1621. 10.1080/07315724.2005.10719438

  • 146

    PiñeroJ.Ramírez-AnguitaJ. M.Saüch-PitarchJ.RonzanoF.CentenoE.SanzF.et al (2020). The DisGeNET Knowledge Platform for Disease Genomics: 2019 Update. Nucleic Acids Res.48 (D1), D845d855. 10.1093/nar/gkz1021

  • 147

    PleschkaS.WolffT.EhrhardtC.HobomG.PlanzO.RappU. R.et al (2001). Influenza Virus Propagation Is Impaired by Inhibition of the Raf/MEK/ERK Signalling cascade. Nat. Cel. Biol.3 (3), 301305. 10.1038/35060098

  • 148

    PoulsenN. N.von BrunnA.HornumM.Blomberg JensenM. (2020). Cyclosporine and COVID-19: Risk or Favorable?Am. J. Transpl.20 (11), 29752982. 10.1111/ajt.16250

  • 149

    QianX.SongX.HeY.YangZ.SunT.WangJ.et al (2015). CCNB2 Overexpression Is a Poor Prognostic Biomarker in Chinese NSCLC Patients. Biomed. Pharmacother.74, 222227. 10.1016/j.biopha.2015.08.004

  • 150

    QiaoC.QianH.WangJ.ZhaoT.MaP.WangS.et al (2020). PD173074 Blocks G1/S Transition via CUL3-Mediated Ubiquitin Protease in HepG2 and Hep3B Cells. PLoS One15 (6), e0234708. 10.1371/journal.pone.0234708

  • 151

    RabinovichA.MedinaL.PiuraB.SegalS.HuleihelM. (2007). Regulation of Ovarian Carcinoma SKOV-3 Cell Proliferation and Secretion of MMPs by Autocrine IL-6. Anticancer Res.27 (1a), 267272.

  • 152

    RawsonJ. M.LandmanS. R.ReillyC. S.BonnacL.PattersonS. E.ManskyL. M. (2015). Lack of Mutational Hot Spots during Decitabine-Mediated HIV-1 Mutagenesis. Antimicrob. Agents Chemother.59 (11), 68346843. 10.1128/aac.01644-15

  • 153

    RawsonJ. M.RothM. E.XieJ.DalyM. B.ClouserC. L.LandmanS. R.et al (2016). Synergistic Reduction of HIV-1 Infectivity by 5-azacytidine and Inhibitors of Ribonucleotide Reductase. Bioorg. Med. Chem.24 (11), 24102422. 10.1016/j.bmc.2016.03.052

  • 154

    SamprathiM.JayashreeM. (2020). Biomarkers in COVID-19: An Up-To-Date Review. Front. Pediatr.8, 607647. 10.3389/fped.2020.607647

  • 155

    SatarkerS.TomA. A.ShajiR. A.AlosiousA.LuvisM.NampoothiriM. (2021). JAK-STAT Pathway Inhibition and Their Implications in COVID-19 Therapy. Postgrad. Med.133 (5), 489507. 10.1080/00325481.2020.1855921

  • 156

    ShangZ.ChanS. Y.LiuW. J.LiP.HuangW. (2021). Recent Insights into Emerging Coronavirus: SARS-CoV-2. ACS Infect. Dis.7 (6), 13691388. 10.1021/acsinfecdis.0c00646

  • 157

    ShannonP.MarkielA.OzierO.BaligaN. S.WangJ. T.RamageD.et al (2003). Cytoscape: a Software Environment for Integrated Models of Biomolecular Interaction Networks. Genome Res.13 (11), 24982504. 10.1101/gr.1239303

  • 158

    SharmaR.UpadhyayS.BhatB.SinghG.BhattacharyaS.SinghA. (2020). Abiotic Stress Induced miRNA-TF-Gene Regulatory Network: A Structural Perspective. Genomics112 (1), 412422. 10.1016/j.ygeno.2019.03.004

  • 159

    SherrC. J. (1996). Cancer Cell Cycles. Science274 (5293), 16721677. 10.1126/science.274.5293.1672

  • 160

    ShibabawT. (2020). Inflammatory Cytokine: IL-17A Signaling Pathway in Patients Present with COVID-19 and Current Treatment Strategy. J. Inflamm. Res.13, 673680. 10.2147/jir.S278335

  • 161

    SinhaS.KunduC. N. (2021). Cancer and COVID-19: Why Are Cancer Patients More Susceptible to COVID-19?Med. Oncol.38 (9), 101. 10.1007/s12032-021-01553-3

  • 162

    SmolleE.Fink-NeuboeckN.LindenmannJ.Smolle-JuettnerF.PichlerM. (2019). The Biological and Clinical Relevance of Inhibitor of Growth (ING) Genes in Non-small Cell Lung Cancer. Cancers (Basel)11 (8), 1118. 10.3390/cancers11081118

  • 163

    SolimaniF.MeierK.GhoreschiK. (2021). Janus Kinase Signaling as Risk Factor and Therapeutic Target for Severe SARS-CoV-2 Infection. Eur. J. Immunol.51 (5), 10711075. 10.1002/eji.202149173

  • 164

    SongX.WuJ. Q.YuX. F.YangX. S.YangY. (2018). Trichostatin A Inhibits Proliferation of Triple Negative Breast Cancer Cells by Inducing Cell Cycle Arrest and Apoptosis. Neoplasma65 (6), 898906. 10.4149/neo_2018_181212N476

  • 165

    StantchevT. S.MarkovicI.TelfordW. G.ClouseK. A.BroderC. C. (2007). The Tyrosine Kinase Inhibitor Genistein Blocks HIV-1 Infection in Primary Human Macrophages. Virus. Res.123 (2), 178189. 10.1016/j.virusres.2006.09.004

  • 166

    StarchenkoA.LauffenburgerD. A. (2018). In Vivo systems Biology Approaches to Chronic Immune/inflammatory Pathophysiology. Curr. Opin. Biotechnol.52, 916. 10.1016/j.copbio.2018.02.006

  • 167

    StelzerG.RosenN.PlaschkesI.ZimmermanS.TwikM.FishilevichS.et al (2016). The GeneCards Suite: From Gene Data Mining to Disease Genome Sequence Analyses. Curr. Protoc. Bioinformatics54, 1.30.11.30.33. 10.1002/cpbi.5

  • 168

    SterlacciW.FieglM.TzankovA. (2012). Prognostic and Predictive Value of Cell Cycle Deregulation in Non-small-cell Lung Cancer. Pathobiology79 (4), 175194. 10.1159/000336462

  • 169

    SungH.FerlayJ.SiegelR. L.LaversanneM.SoerjomataramI.JemalA.et al (2021). Global Cancer Statistics 2020: GLOBOCAN Estimates of Incidence and Mortality Worldwide for 36 Cancers in 185 Countries. CA Cancer J. Clin.71 (3), 209249. 10.3322/caac.21660

  • 170

    SureshD.ZambreA.MukherjeeS.GhoshdastidarS.JiangY.JoshiT.et al (2019). Silencing AXL by Covalent siRNA-Gelatin-Antibody Nanoconjugate Inactivates mTOR/EMT Pathway and Stimulates P53 for TKI Sensitization in NSCLC. Nanomedicine20, 102007. 10.1016/j.nano.2019.04.010

  • 171

    TakashimaS.SaitoH.TakahashiN.ImaiK.KudoS.AtariM.et al (2014). Strong Expression of Cyclin B2 mRNA Correlates with a Poor Prognosis in Patients with Non-small Cell Lung Cancer. Tumour Biol.35 (5), 42574265. 10.1007/s13277-013-1556-7

  • 172

    TangS. M.DengX. T.ZhouJ.LiQ. P.GeX. X.MiaoL. (2020). Pharmacological Basis and New Insights of Quercetin Action in Respect to its Anti-cancer Effects. Biomed. Pharmacother.121, 109604. 10.1016/j.biopha.2019.109604

  • 173

    TayM. Z.PohC. M.RéniaL.MacAryP. A.NgL. F. P. (2020). The trinity of COVID-19: Immunity, Inflammation and Intervention. Nat. Rev. Immunol.20 (6), 363374. 10.1038/s41577-020-0311-8

  • 174

    TianS.HuW.NiuL.LiuH.XuH.XiaoS. Y. (2020). Pulmonary Pathology of Early-phase 2019 Novel Coronavirus (COVID-19) Pneumonia in Two Patients with Lung Cancer. J. Thorac. Oncol.15 (5), 700704. 10.1016/j.jtho.2020.02.010

  • 175

    TianqinG.ChunleiC.JingjingW. (2016). Synergistic Anti-glioma Effects In Vitro and In Vivo of Enediyne Antibiotic Neocarzinostatin and Paclitaxel via Enhanced Growth Delay and Apoptosis-Induction. Biol. Pharm. Bull.39 (10), 16231630. 10.1248/bpb.b16-00285

  • 176

    UmakanthanS.SahuP.RanadeA. V.BukeloM. M.RaoJ. S.Abrahao-MachadoL. F.et al (2020). Origin, Transmission, Diagnosis and Management of Coronavirus Disease 2019 (COVID-19). Postgrad. Med. J.96 (1142), 753758. 10.1136/postgradmedj-2020-138234

  • 177

    UngvariZ.TarantiniS.KissT.WrenJ. D.GilesC. B.GriffinC. T.et al (2018). Endothelial Dysfunction and Angiogenesis Impairment in the Ageing Vasculature. Nat. Rev. Cardiol.15 (9), 555565. 10.1038/s41569-018-0030-z

  • 178

    van DrielM. A.BrunnerH. G. (2006). Bioinformatics Methods for Identifying Candidate Disease Genes. Hum. Genomics2 (6), 429432. 10.1186/1479-7364-2-6-429

  • 179

    VonoM.HuttnerA.LemeilleS.Martinez-MurilloP.MeyerB.BaggioS.et al (2021). Robust Innate Responses to SARS-CoV-2 in Children Resolve Faster Than in Adults without Compromising Adaptive Immunity. Cell Rep.37 (1), 109773. 10.1016/j.celrep.2021.109773

  • 180

    WaldO.IzharU.AmirG.KirshbergS.ShlomaiZ.ZamirG.et al (2011). Interaction between Neoplastic Cells and Cancer-Associated Fibroblasts through the CXCL12/CXCR4 axis: Role in Non-small Cell Lung Cancer Tumor Proliferation. J. Thorac. Cardiovasc. Surg.141 (6), 15031512. 10.1016/j.jtcvs.2010.11.056

  • 181

    WallsA. C.ParkY. J.TortoriciM. A.WallA.McGuireA. T.VeeslerD. (2020). Structure, Function, and Antigenicity of the SARS-CoV-2 Spike Glycoprotein. Cell181 (2), 281e6. e286. 10.1016/j.cell.2020.02.058

  • 182

    WangM.HerbstR. S.BoshoffC. (2021). Toward Personalized Treatment Approaches for Non-small-cell Lung Cancer. Nat. Med.27 (8), 13451356. 10.1038/s41591-021-01450-2

  • 183

    WangX.PowellC. A. (2021). How to Translate the Knowledge of COVID‐19 into the Prevention of Omicron Variants. Clin. Translational Med11 (12), e680. 10.1002/ctm2.680

  • 184

    WangY.ZhangS.LiF.ZhouY.ZhangY.WangZ.et al (2020). Therapeutic Target Database 2020: Enriched Resource for Facilitating Research and Early Development of Targeted Therapeutics. Nucleic Acids Res.48 (D1), D1031d1041. 10.1093/nar/gkz981

  • 185

    WeiD.SunL.FengW. (2021). hsa_circ_0058357 Acts as a ceRNA to Promote Non-small C-ell L-ung C-ancer P-rogression via the hsa-miR-24-3p/AVL9 axis. Mol. Med. Rep.23 (6), 470. 10.3892/mmr.2021.12109

  • 186

    WenL.TangK.ChikK. K.ChanC. C.TsangJ. O.LiangR.et al (2021). In Silico structure-based Discovery of a SARS-CoV-2 Main Protease Inhibitor. Int. J. Biol. Sci.17 (6), 15551564. 10.7150/ijbs.59191

  • 187

    WilliamsG. H.StoeberK. (2012). The Cell Cycle and Cancer. J. Pathol.226 (2), 352364. 10.1002/path.3022

  • 188

    WoollerS. K.Benstead-HumeG.ChenX.AliY.PearlF. M. G. (2017). Bioinformatics in Translational Drug Discovery. Biosci. Rep.37 (4), BSR20160180. 10.1042/bsr20160180

  • 189

    Worldometer (2021). COVID 19 Coronavirus Pandemic. [Online]. Available: https://www.worldometers.info/coronavirus/ (Accessed November 21, 2021).

  • 190

    WuT.HuE.XuS.ChenM.GuoP.DaiZ.et al (2021). clusterProfiler 4.0: A Universal Enrichment Tool for Interpreting Omics Data. The Innovation2 (3), 100141. 10.1016/j.xinn.2021.100141

  • 191

    WuW.LiR.LiX.HeJ.JiangS.LiuS.et al (2015). Quercetin as an Antiviral Agent Inhibits Influenza A Virus (IAV) Entry. Viruses8 (1), 6. 10.3390/v8010006

  • 192

    WuX. Y.ChenH. C.LiW. W.YanJ. D.LvR. Y. (2020). DNMT1 Promotes Cell Proliferation via Methylating hMLH1 and hMSH2 Promoters in EGFR-Mutated Non-small Cell Lung Cancer. J. Biochem.168 (2), 151157. 10.1093/jb/mvaa034

  • 193

    XingX. L.LuY.QiuH. L. (2016). Effect of Cyclosporin A Particles of Varying Diameters on Gastric Cancer Cell Apoptosis. Genet. Mol. Res.15 (2), gmr8085. 10.4238/gmr.15028085

  • 194

    XuC.BuczkowskiK. A.ZhangY.AsahinaH.BeauchampE. M.TeraiH.et al (2015). NSCLC Driven by DDR2 Mutation Is Sensitive to Dasatinib and JQ1 Combination Therapy. Mol. Cancer Ther.14 (10), 23822389. 10.1158/1535-7163.Mct-15-0077

  • 195

    XuD.HuM. J.WangY. Q.CuiY. L. (2019). Antioxidant Activities of Quercetin and its Complexes for Medicinal Application. Molecules24 (6), 1123. 10.3390/molecules24061123

  • 196

    YangG.RosenD. G.ZhangZ.BastR. C.Jr.MillsG. B.ColacinoJ. A.et al (2006). The Chemokine Growth-Regulated Oncogene 1 (Gro-1) Links RAS Signaling to the Senescence of Stromal Fibroblasts and Ovarian Tumorigenesis. Proc. Natl. Acad. Sci. U S A.103 (44), 1647216477. 10.1073/pnas.0605752103

  • 197

    YaoX.YeF.ZhangM.CuiC.HuangB.NiuP.et al (2020). In Vitro Antiviral Activity and Projection of Optimized Dosing Design of Hydroxychloroquine for the Treatment of Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2). Clin. Infect. Dis.71 (15), 732739. 10.1093/cid/ciaa237

  • 198

    YooM.ShinJ.KimJ.RyallK. A.LeeK.LeeS.et al (2015). DSigDB: Drug Signatures Database for Gene Set Analysis. Bioinformatics31 (18), 30693071. 10.1093/bioinformatics/btv313

  • 199

    YousefM.VlachogiannisI. A.TsianiE. (2017). Effects of Resveratrol against Lung Cancer: In Vitro and In Vivo Studies. Nutrients9 (11), 1231. 10.3390/nu9111231

  • 200

    YuG.WangW.WangX.XuM.ZhangL.DingL.et al (2018). Network Pharmacology-Based Strategy to Investigate Pharmacological Mechanisms of Zuojinwan for Treatment of Gastritis. BMC Complement. Altern. Med.18 (1), 292. 10.1186/s12906-018-2356-9

  • 201

    YuanS.YuZ.LiuQ.ZhangM.XiangY.WuN.et al (2016). GPC5, a Novel Epigenetically Silenced Tumor Suppressor, Inhibits Tumor Growth by Suppressing Wnt/β-Catenin Signaling in Lung Adenocarcinoma. Oncogene35 (47), 61206131. 10.1038/onc.2016.149

  • 202

    ZappaC.MousaS. A. (2016). Non-small Cell Lung Cancer: Current Treatment and Future Advances. Transl. Lung Cancer Res.5 (3), 288300. 10.21037/tlcr.2016.06.07

  • 203

    ZhanW.WangW.HanT.XieC.ZhangT.GanM.et al (2017). COMMD9 Promotes TFDP1/E2F1 Transcriptional Activity via Interaction with TFDP1 in Non-small Cell Lung Cancer. Cell Signal30, 5966. 10.1016/j.cellsig.2016.11.016

  • 204

    ZhangH.ShaoH.GolubovskayaV. M.ChenH.CanceW.AdjeiA. A.et al (2016). Efficacy of Focal Adhesion Kinase Inhibition in Non-small Cell Lung Cancer with Oncogenically Activated MAPK Pathways. Br. J. Cancer115 (2), 203211. 10.1038/bjc.2016.190

  • 205

    ZhangH. M.KuangS.XiongX.GaoT.LiuC.GuoA. Y. (2015). Transcription Factor and microRNA Co-regulatory Loops: Important Regulatory Motifs in Biological Processes and Diseases. Brief Bioinform16 (1), 4558. 10.1093/bib/bbt085

  • 206

    ZhangJ.LiD.ZhangY.DingZ.ZhengY.ChenS.et al (2020). Integrative Analysis of mRNA and miRNA Expression Profiles Reveals Seven Potential Diagnostic Biomarkers for Non-small C-ell L-ung C-ancer. Oncol. Rep.43 (1), 99112. 10.3892/or.2019.7407

  • 207

    ZhangJ. J.HongJ.MaY. S.ShiY.ZhangD. D.YangX. L.et al (2021). Identified GNGT1 and NMU as Combined Diagnosis Biomarker of Non-small-cell Lung Cancer Utilizing Bioinformatics and Logistic Regression. Dis. Markers2021, 6696198. 10.1155/2021/6696198

  • 208

    ZhangL.HanX.ShiY. (2020). Comparative Analysis of SARS-CoV-2 Receptor ACE2 Expression in Multiple Solid Tumors and Matched Non-diseased Tissues. Infect. Genet. Evol.85, 104428. 10.1016/j.meegid.2020.104428

  • 209

    ZhangL.HaoC.YaoS.TangR.GuoW.CongH.et al (2018). Exosomal miRNA Profiling to Identify Nanoparticle Phagocytic Mechanisms. Small14 (15), e1704008. 10.1002/smll.201704008

  • 210

    ZhangZ.WangQ.ChenF.LiuJ. (2015b). Elevated Expression of HMGA1 Correlates with the Malignant Status and Prognosis of Non-small Cell Lung Cancer. Tumour Biol.36 (2), 12131219. 10.1007/s13277-014-2749-4

  • 211

    ZhaoW.BaoP.QiH.YouH. (2010). Resveratrol Down-Regulates Survivin and Induces Apoptosis in Human Multidrug-Resistant SPC-A-1/CDDP Cells. Oncol. Rep.23 (1), 279286.

  • 212

    ZhongW.YanT.WebberM. M.OberleyT. D. (2004). Alteration of Cellular Phenotype and Responses to Oxidative Stress by Manganese Superoxide Dismutase and a Superoxide Dismutase Mimic in RWPE-2 Human Prostate Adenocarcinoma Cells. Antioxid. Redox Signal.6 (3), 513522. 10.1089/152308604773934279

  • 213

    ZhouJ.YangQ.LuL.TuoZ.ShouZ.ChengJ. (2021). PLK1 Inhibition Induces Immunogenic Cell Death and Enhances Immunity against NSCLC. Int. J. Med. Sci.18 (15), 35163525. 10.7150/ijms.60135

  • 214

    ZhouP.YangX. L.WangX. G.HuB.ZhangL.ZhangW.et al (2020). A Pneumonia Outbreak Associated with a New Coronavirus of Probable Bat Origin. Nature579 (7798), 270273. 10.1038/s41586-020-2012-7

  • 215

    ZhouX.ZhangZ.LiangX. (2020). Regulatory Network Analysis to Reveal Important miRNAs and Genes in Non-small Cell Lung Cancer. Cell J21 (4), 459466. 10.22074/cellj.2020.6281

  • 216

    ZhuangZ.ChenQ.HuangC.WenJ.HuangH.LiuZ. (2020). A Comprehensive Network Pharmacology-Based Strategy to Investigate Multiple Mechanisms of HeChan Tablet on Lung Cancer. Evid. Based Complement. Alternat. Med.2020, 7658342. 10.1155/2020/7658342

Summary

Keywords

COVID-19, non-small-cell lung cancer, interaction network, bioinformatic analysis, systemic biological analysis

Citation

Zhuang Z, Zhong X, Chen Q, Chen H and Liu Z (2022) Bioinformatics and System Biology Approach to Reveal the Interaction Network and the Therapeutic Implications for Non-Small Cell Lung Cancer Patients With COVID-19. Front. Pharmacol. 13:857730. doi: 10.3389/fphar.2022.857730

Received

19 January 2022

Accepted

28 March 2022

Published

02 June 2022

Volume

13 - 2022

Edited by

Francesca Caccuri, University of Brescia, Italy

Reviewed by

Ettore Mosca, National Research Council, Italy

Sucheendra K. Palaniappan, The Systems Biology Institute, Japan

Updates

Copyright

*Correspondence: Zhanhua Liu,

This article was submitted to Pharmacology of Infectious Diseases, a section of the journal Frontiers in Pharmacology

Disclaimer

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.

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics