Comprehensive analysis of differential expression profiles reveals potential biomarkers associated with the cell cycle and regulated by p53 in human small cell lung cancer

Identifying General Tumor and Specific Lung Cancer Biomarkers by Transcriptomic Analysis

The bioinformatic pipeline previously developed in our research laboratory is used to identify potential general and specific deregulated tumor genes and transcription factors related to the establishment and progression of tumoral diseases, now comparing lung cancer with other two types of cancer. Twenty microarray datasets were selected and analyzed separately to identify hub differentiated expressed genes and compared to identify all the deregulated genes and transcription factors in common between the three types of cancer and those unique to lung cancer. The winning DEGs analysis allowed to identify an important number of TFs deregulated in the majority of microarray datasets, which can become key biomarkers of general tumors and specific to lung cancer. A coexpression network was constructed for every dataset with all deregulated genes associated with lung cancer, according to DAVID’s tool enrichment analysis, and transcription factors capable of regulating them, according to oPOSSUM´s tool. Several genes and transcription factors are coexpressed in the networks, suggesting that they could be related to the establishment or progression of the tumoral pathology in any tissue and specifically in the lung. The comparison of the coexpression networks of lung cancer and other types of cancer allowed the identification of common connectivity patterns with deregulated genes and transcription factors correlated to important tumoral processes and signaling pathways that have not been studied yet to experimentally validate their role in lung cancer. The Kaplan–Meier estimator determined the association of thirteen deregulated top winning transcription factors with the survival of lung cancer patients. The coregulatory analysis identified two top winning transcription factors networks related to the regulatory control of gene expression in lung and breast cancer. Our transcriptomic analysis suggests that cancer has an important coregulatory network of transcription factors related to the acquisition of the hallmarks of cancer. Moreover, lung cancer has a group of genes and transcription factors unique to pulmonary tissue that are coexpressed during tumorigenesis and must be studied experimentally to fully understand their role in the pathogenesis within its very complex transcriptomic scenario. Therefore, the downstream bioinformatic analysis developed was able to identify a coregulatory metafirm of cancer in general and specific to lung cancer taking into account the great heterogeneity of the tumoral process at cellular and population levels.

Disclosing Potential Key Genes, Therapeutic Targets and Agents for Non-Small Cell Lung Cancer: Evidence from Integrative Bioinformatics Analysis

Non-small-cell lung cancer (NSCLC) is considered as one of the malignant cancers that causes premature death. The present study aimed to identify a few potential novel genes highlighting their functions, pathways, and regulators for diagnosis, prognosis, and therapies of NSCLC by using the integrated bioinformatics approaches. At first, we picked out 1943 DEGs between NSCLC and control samples by using the statistical LIMMA approach. Then we selected 11 DEGs (CDK1, EGFR, FYN, UBC, MYC, CCNB1, FOS, RHOB, CDC6, CDC20, and CHEK1) as the hub-DEGs (potential key genes) by the protein–protein interaction network analysis of DEGs. The DEGs and hub-DEGs regulatory network analysis commonly revealed four transcription factors (FOXC1, GATA2, YY1, and NFIC) and five miRNAs (miR-335-5p, miR-26b-5p, miR-92a-3p, miR-155-5p, and miR-16-5p) as the key transcriptional and post-transcriptional regulators of DEGs as well as hub-DEGs. We also disclosed the pathogenetic processes of NSCLC by investigating the biological processes, molecular function, cellular components, and KEGG pathways of DEGs. The multivariate survival probability curves based on the expression of hub-DEGs in the SurvExpress web-tool and database showed the significant differences between the low- and high-risk groups, which indicates strong prognostic power of hub-DEGs. Then, we explored top-ranked 5-hub-DEGs-guided repurposable drugs based on the Connectivity Map (CMap) database. Out of the selected drugs, we validated six FDA-approved launched drugs (Dinaciclib, Afatinib, Icotinib, Bosutinib, Dasatinib, and TWS-119) by molecular docking interaction analysis with the respective target proteins for the treatment against NSCLC. The detected therapeutic targets and repurposable drugs require further attention by experimental studies to establish them as potential biomarkers for precision medicine in NSCLC treatment.

Identification of CCNB2 as A Potential Non-Invasive Breast Cancer Biomarker in Peripheral Blood Mononuclear Cells Using The Systems Biology Approach

OBJECTIVE

Breast cancer (BC) still remains an imperative clinical issue, despite advances in the diagnosis, prognosis and treatment modalities of this malignancy. Hence, progress has been made to identify non-invasive, high sensitive and specific biomarkers. Since immune system affects development of breast cancer, peripheral blood mononuclear cells (PBMCs) -a subpopulation of immune cells- can be considered as a promising tool in the field of BC biomarker research. In the current study, we initially attempted to use concept of the present shared biomarkers in solid tumors and systemic immune profile and then evaluate correlation of these biomarkers to clinical use in cancer research.

MATERIALS AND METHODS

In this experimental study, available microarray gene expression datasets of BC as well as the related PBMCs were retrieved and downloaded from the Gene Expression Omnibus (GEO) database, followed by analysis using GEO2R along with affylmGUI, a R-based package, to obtain differentially expressed genes (DEGs). Signature genes from 20 types of cancer were also applied to validate DEGs. Quantitative reverse-transcription polymerase chain reaction (qRT-PCR) was carried out to assess mRNA level of CCNB2 in PBMC of the BC patients and healthy subjects.

RESULTS

DEGs analysis for the transcription profile of BC cells and PBMCs showed two shared targets, CCNB2 and PGK1. Validation with systems biology using reweighted 20 types of cancer signature genes revealed that CCNB2 is the only common target in BC and its related PBMCs, which was further validated by qRT-PCR implying a significant increase in the level of CCNB2 in the BC patients.

CONCLUSION

Results of this study demonstrated that PBMCs are affected by BC cells and CCNB2 may be of value as a diagnostic biomarker for breast cancer. However, verification would require future detailed experimental plans.

Computational identification of biomarker genes for lung cancer considering treatment and non-treatment studies

BACKGROUND

Lung cancer is the number one cancer killer in the world with more than 142,670 deaths estimated in the United States alone in the year 2019. Consequently, there is an overreaching need to identify the key biomarkers for lung cancer. The aim of this study is to computationally identify biomarker genes for lung cancer that can aid in its diagnosis and treatment. The gene expression profiles of two different types of studies, namely non-treatment and treatment, are considered for discovering biomarker genes. In non-treatment studies healthy samples are control and cancer samples are cases. Whereas, in treatment studies, controls are cancer cell lines without treatment and cases are cancer cell lines with treatment.

RESULTS

The Differentially Expressed Genes (DEGs) for lung cancer were isolated from Gene Expression Omnibus (GEO) database using R software tool GEO2R. A total of 407 DEGs (254 upregulated and 153 downregulated) from non-treatment studies and 547 DEGs (133 upregulated and 414 downregulated) from treatment studies were isolated. Two Cytoscape apps, namely, CytoHubba and MCODE, were used for identifying biomarker genes from functional networks developed using DEG genes. This study discovered two distinct sets of biomarker genes - one from non-treatment studies and the other from treatment studies, each set containing 16 genes. Survival analysis results show that most non-treatment biomarker genes have prognostic capability by indicating low-expression groups have higher chance of survival compare to high-expression groups. Whereas, most treatment biomarkers have prognostic capability by indicating high-expression groups have higher chance of survival compare to low-expression groups.

CONCLUSION

A computational framework is developed to identify biomarker genes for lung cancer using gene expression profiles. Two different types of studies - non-treatment and treatment - are considered for experiment. Most of the biomarker genes from non-treatment studies are part of mitosis and play vital role in DNA repair and cell-cycle regulation. Whereas, most of the biomarker genes from treatment studies are associated to ubiquitination and cellular response to stress. This study discovered a list of biomarkers, which would help experimental scientists to design a lab experiment for further exploration of detail dynamics of lung cancer development.

Identification of potential target genes and crucial pathways in small cell lung cancer based on bioinformatic strategy and human samples

Small cell lung cancer (SCLC) is a carcinoma of the lungs with strong invasion, poor prognosis and resistant to multiple chemotherapeutic drugs. It has posed severe challenges for the effective treatment of lung cancer. Therefore, searching for genes related to the development and prognosis of SCLC and uncovering their underlying molecular mechanisms are urgent problems to be resolved. This study is aimed at exploring the potential pathogenic and prognostic crucial genes and key pathways of SCLC via bioinformatic analysis of public datasets. Firstly, 117 SCLC samples and 51 normal lung samples were collected and analyzed from three gene expression datasets. Then, 102 up-regulated and 106 down-regulated differentially expressed genes (DEGs) were observed. And then, functional annotation and pathway enrichment analyzes of DEGs was performed utilizing the FunRich. The protein-protein interaction (PPI) network of the DEGs was constructed through the STRING website, visualized by Cytoscape. Finally, the expression levels of eight hub genes were confirmed in Oncomine database and human samples from SCLC patients. It showed that CDC20, BUB1, TOP2A, RRM2, CCNA2, UBE2C, MAD2L1, and BUB1B were upregulated in SCLC tissues compared to paired adjacent non-cancerous tissues. These suggested that eight hub genes might be viewed as new biomarkers for prognosis of SCLC or to guide individualized medication for the therapy of SCLC.

An in-silico method leads to recognition of hub genes and crucial pathways in survival of patients with breast cancer

Breast cancer is a highly heterogeneous disorder characterized by dysregulation of expression of numerous genes and cascades. In the current study, we aim to use a system biology strategy to identify key genes and signaling pathways in breast cancer. We have retrieved data of two microarray datasets (GSE65194 and GSE45827) from the NCBI Gene Expression Omnibus database. R package was used for identification of differentially expressed genes (DEGs), assessment of gene ontology and pathway enrichment evaluation. The DEGs were integrated to construct a protein-protein interaction network. Next, hub genes were recognized using the Cytoscape software and lncRNA-mRNA co-expression analysis was performed to evaluate the potential roles of lncRNAs. Finally, the clinical importance of the obtained genes was assessed using Kaplan-Meier survival analysis. In the present study, 887 DEGs including 730 upregulated and 157 downregulated DEGs were detected between breast cancer and normal samples. By combining the results of functional analysis, MCODE, CytoNCA and CytoHubba 2 hub genes including MAD2L1 and CCNB1 were selected. We also identified 12 lncRNAs with significant correlation with MAD2L1 and CCNB1 genes. According to The Kaplan-Meier plotter database MAD2L1, CCNA2, RAD51-AS1 and LINC01089 have the most prediction potential among all candidate hub genes. Our study offers a framework for recognition of mRNA-lncRNA network in breast cancer and detection of important pathways that could be used as therapeutic targets in this kind of cancer.

Identification of APEX2 as an oncogene in liver cancer

BACKGROUND

DNA damage is one of the critical contributors to the occurrence and development of some cancers. APEX1 and APEX2 are the most important molecules in the DNA damage, and APEX1 has been identified as a diagnostic and prognostic biomarker in liver hepatocellular carcinoma (LIHC). However, the expression of APEX2 and its functional mechanisms in LIHC are still unclear.

AIM

To examine the expression of APEX2 and the potential mechanism network in LIHC.

METHODS

We conducted a pan-cancer analysis of the expression of APEX1 and APEX2 using the interactive TIMER tool. GEO datasets, including GSE14520, GSE22058, and GSE64041, were used to compare the APEX2 expression level in tumor tissues and adjacent non-tumor tissues. Then, we calculated the 5-year survival rate according to the web-based Kaplan-Meier analysis. We included the TCGA liver cancer database in GSEA analysis based on the high and low APEX2 expression, showing the potential mechanisms of APEX2 in LIHC. After that, we conducted Pearson correlation analysis using GEPIA2. Next, we performed quantitative polymerase chain reaction (qPCR) assay to examine the APEX2 levels in normal liver cell line LO2 and several liver cancer cell lines, including HepG2, Huh7, SMMC7721, and HCCLM3. APEX2 in HCCLM3 cells was knocked down using small interfering RNA. The role of APEX2 in cell viability was confirmed using CCK-8. Dual-luciferase reporter assay was performed to examine the promoter activity of CCNB1 and MYC.

RESULTS

APEX1 and APEX2 are both highly expressed in the tumor tissues of BLCA, BRCA, CHOL, COAD, ESCA, HNSC, LIHC, LUAD, LUSC, READ, and STAD. APEX2 overexpression in LIHC was validated using GSE14520, GSE22058, and GSE64041 datasets. The survival analysis showed that LIHC patients with high expression of APEX2 had a lower overall survival rate, even in the AJCC T1 patients. High level of APEX2 could indicate a lower overall survival rate in patients with or without viral hepatitis. The GSEA analysis identified that kinetochore and spindle microtubules are the two main cellular components of APEX2 in GO Ontology. APEX2 was also positively associated with molecular function regulation of chromosome segregation and DNA replication. The results of KEGG analysis indicated that APEX2 expression was positively correlated with cell cycle pathway and pro-oncogenic MYC signaling. Pearson correlation analysis showed that APEX2 had a significant positive correlation with CCNB1 and MYC. APEX2 level was higher in liver cancer cell lines than in normal liver LO2 cells. Small interfering RNA could knock down the APEX2 expression in HCCLM3 cells. Knockdown of APEX2 resulted in a decrease in the viability of HCCLM3 cells as well as the expression and promoter activity of CCNB1 and MYC.

CONCLUSION

APEX2 is overexpressed in LIHC, and the higher APEX2 level is associated with a worse prognosis in overall survival. APEX2 is closely involved in the biological processes of chromosome segregation and DNA replication. APEX2 expression is positively correlated with the pro-oncogenic pathways. Knockdown of APEX2 could inhibit the cell viability and CCNB1 and MYC pathways, suggesting that APEX2 is an oncogene in LIHC, which could be a potential pharmaceutic target in the anti-tumor therapy.

Characterization of lncRNA-Associated ceRNA Network to Reveal Potential Prognostic Biomarkers in Lung Adenocarcinoma

Lung adenocarcinoma (LUAD) is one of the most fatal malignant tumors harmful to human health. The complexity and behavior characteristics of long-non-coding RNA (lncRNA)-associated competing endogenous RNA (ceRNA) network in LUAD patients are still unclear. The purpose of this study was to elucidate the regulatory networks of dysregulated RNAs, view, and identify potential prognosis signatures involved in LUAD. The expression profiles of mRNAs, lncRNAs, and miRNAs were obtained from the TCGA database. In total, 2078 DEmRNAs, 257 DElncRNAs, and 101 DEmiRNAs were sorted out. A PPI network including 45 DEmRNAs was constructed. Ten hub genes in the PPI network associated with cell cycle-related pathways were identified and they played key roles in regulating cell proliferation. A total of three DEmiRNAs, seven DElncRNAs, and six DEmRNAs were enrolled in the ceRNA network. Except for certain genes without any published study reports, all the genes in the ceRNA network played an essential role in controlling tumor cell proliferation and were associated with prognosis in LUAD. Finally, based on step regression and Cox regression survival analysis, we identified four candidate biomarkers, including miR490, miR1293, LINC01740, and IGF2BP1, and established a risk model based on the four genes. Our study provided a global view and systematic dissection of the lncRNA-associated ceRNA network, and the identified four genes might be novel important prognostic factors involved in LUAD pathogenesis.

Molecular Signature of Subtypes of Non-Small-Cell Lung Cancer by Large-Scale Transcriptional Profiling: Identification of Key Modules and Genes by Weighted Gene Co-Expression Network Analysis (WGCNA)

Non-small-cell lung cancer (NSCLC) represents a heterogeneous group of malignancies consisting essentially of adenocarcinoma (ADC) and squamous cell carcinoma (SCC). Although the diagnosis and treatment of ADC and SCC have been greatly improved in recent decades, there is still an urgent need to identify accurate transcriptome profile associated with the histological subtypes of NSCLC. The present study aims to identify the key dysregulated pathways and genes involved in the development of lung ADC and SCC and to relate them with the clinical traits. The transcriptional changes between tumour and normal lung tissues were investigated by RNA-seq. Gene ontology (GO), canonical pathways analysis with the prediction of upstream regulators, and weighted gene co-expression network analysis (WGCNA) to identify co-expressed modules and hub genes were used to explore the biological functions of the identified dysregulated genes. It was indicated that specific gene signatures differed significantly between ADC and SCC related to the distinct pathways. Of identified modules, four and two modules were the most related to clinical features in ADC and SCC, respectively. CTLA4, MZB1, NIP7, and BUB1B in ADC, as well as GNG11 and CCNB2 in SCC, are novel top hub genes in modules associated with tumour size, SUVmax, and recurrence-free survival. Our research provides a more effective understanding of the importance of biological pathways and the relationships between major genes in NSCLC in the perspective of searching for new molecular targets.

Identification of candidate genes associated with the pathogenesis of small cell lung cancer via integrated bioinformatics analysis

The pathogenesis of small cell lung cancer (SCLC), a highly metastatic malignant tumor, remains unclear. In the present study, important genes and pathways that are involved in the pathogenesis of SCLC were identified. The following four datasets were downloaded from the Gene Expression Omnibus: GSE60052, GSE43346, GSE15240 and GSE6044. The differentially expressed genes (DEGs) between the SCLC samples and the normal samples were analyzed using R software. The limma package was used for every dataset. The RobustRankAggreg package was used to integrate the DEGs from the four datasets. Functional and pathway enrichment analyses were conducted using the Gene Ontology and Kyoto Encyclopedia of Genes and Genomes databases with FunRich software and R software, respectively. In addition, the protein-protein interaction (PPI) network of the DEGs was constructed using the STRING database and Cytoscape software. Hub genes and significant modules were identified using Molecular Complex Detection in Cytoscape software. Finally, the expression values of hub genes were determined using the Oncomine online database. In total, 412 DEGs were identified following the integration of the four datasets, with 146 upregulated genes and 266 downregulated genes. The upregulated DEGs were primarily enriched in the cell cycle, cell division and microtubule binding. The downregulated DEGs were primarily enriched in the complement and coagulation cascades, the cytokine-mediated signaling pathway and protein binding. Eight hub genes and 1 significant module correlated to the cell cycle pathway were identified based on a subset of the PPI network. Finally, five hub genes were identified as highly expressed in SCLC tissue compared with normal tissue. The cell cycle pathway may be the pathway most closely associated with the pathogenesis of SCLC. NDC80, BUB1B, PLK1, CDC20 and MAD2L1 should be the focus of follow-up studies regarding the diagnosis and treatment of SCLC.

Identification of potential agents for thymoma by integrated analyses of differentially expressed tumour-associated genes and molecular docking experiments

Thymoma, derived from the epithelial cells of the thymus, is a rare malignant tumour type. Following diagnosis with thymoma, patients generally undergo surgical treatment. However, patients with advanced-stage disease are only candidates for chemotherapy and have poor survival. Therefore, it is urgently required to explore effective chemotherapeutic agents for the treatment of thymoma. In the present study, a Bioinformatics analysis was performed to identify novel drugs for thymoma. Differentially expressed genes (DEGs) in thymoma were obtained by Gene Expression Profiling Interactive Analysis. Subsequently, these genes were processed by Connectivity Map analysis to identify suitable compounds. In addition, Metascape software was used to verify drug and target binding. Molecular docking technology was used to verify drug and target binding. Finally, absorption, distribution, metabolism and excretion parameters in the Traditional Chinese Medicine Systems Pharmacology Database and Analysis Platform database were used for drug screening and for evaluation of the potential clinical value. In total, 2,447 DEGs, including 2,204 upregulated and 243 downregulated genes, were identified from 118 thymoma patients and 339 normal samples. The top 10 drugs displaying the most significant negative correlations were fulvestrant, hesperetin, zidovudine, hydrocortisone, rolitetracycline, ellipticine, sirolimus, quinisocaine, oestradiol (estradiol) and harmine. The predicted targets of these drugs were then confirmed. The score for the association between estrogen receptor 1 (ESR1) and fulvestrant was 0.99. According to the molecular docking analysis, the total scores for the interaction between ESR1 were 10.26, and those for the interaction between tamoxifen and ESR1 were 6.60. The oral bioavailability (%), drug-likeness and drug half-life for hesperetin were 70.31, 0.27 and 15.78, respectively; those for oestradiol were 53.56, 0.32 and 3.50, respectively; and those for harmine were 56.80, 0.13 and 5.04, respectively. In conclusion, several potential therapeutic drugs for thymoma were identified in the present study. The results suggested that the compounds, including fulvestrant, estradiol, hesperetin and ellipticine, represent the most likely drugs for the treatment of thymoma. Future studies should focus on testing these novel compounds in vitro and in vivo.

Lung Cancer: One Disease or Many

OBJECTIVE

Lung cancer is classified as a single entity comprised of multiple histological subtypes. But how similar are these subtypes on a genetic level? This paper aims to address this question through a concise overview of germline and somatic differences between small cell lung cancer, lung adenocarcinoma, and lung squamous cell carcinoma.

METHODS

We reveal the weak overlap found between these 3 lung cancer subtypes using published data from one of the largest germline genetic studies on lung cancer to date and somatic mutation data from Catalogue of Somatic Mutations in Cancer (COSMIC).

RESULTS

These data indicate that these 3 subtypes share very little with each other at the genetic level. At the germline SNP level, only 24 independent SNPs from 2 chromosomes were shared across all 3 subtypes. We also demonstrate that only 30 unique cancer-specific mutations overlap the 3 subtypes from COSMIC, and that this is fewer than overlapping mutations chosen at random. Finally, we show that only 3 somatic mutational signatures are shared between these 3 subtypes.

CONCLUSION

This paper highlights that these 3 lung cancer subtypes may be distinct diseases at the genetic level. In the era of precision medicine, we feel that these genomic differences will be of utmost importance in the choice of lung cancer therapy in the future.