J Cancer 2020; 11(2):388-402. doi:10.7150/jca.31636

Research Paper

TBX2 Identified as a Potential Predictor of Bone Metastasis in Lung Adenocarcinoma via Integrated Bioinformatics Analyses and Verification of Functional Assay

Huajian Yu1*, Fangyu Zhao1*, Jing Li1, Kechao Zhu2, Hechun Lin1, Zhen Pan2, Miaoxin Zhu1, Ming Yao1✉, Mingxia Yan1,3✉

1. State Key Laboratory of Oncogenes and Related Genes, Shanghai Cancer Institute, Renji Hospital, Shanghai Jiao Tong University School of Medicine, Shanghai, China;
2. Department of Orthopedics, Shanghai Jiao Tong University Affiliated Sixth People's Hospital, Shanghai, China.
3. Fudan University Shanghai Cancer Center, Shanghai Medical College, Fudan University, Shanghai, China.
*These authors contributed equally to this work.

This is an open access article distributed under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0/). See http://ivyspring.com/terms for full terms and conditions.
Citation:
Yu H, Zhao F, Li J, Zhu K, Lin H, Pan Z, Zhu M, Yao M, Yan M. TBX2 Identified as a Potential Predictor of Bone Metastasis in Lung Adenocarcinoma via Integrated Bioinformatics Analyses and Verification of Functional Assay. J Cancer 2020; 11(2):388-402. doi:10.7150/jca.31636. Available from https://www.jcancer.org/v11p0388.htm

File import instruction

Abstract

Objective: Bone metastasis from patients with advanced lung adenocarcinoma (LAC) is a very serious complication. To better understand the molecular mechanism, our current study sheds light on identification of hub genes mediating bone metastatic spread by combining bioinformatic analysis with functional verification.

Methods: First, we downloaded a lung adenocarcinoma dataset (GSE76194) from Gene Expression Omnibus, analyzed differentially expressed genes (DEGs) through Limma package in R software and constructed a protein-protein interaction network. Based on that preliminary data, we further performed modular and topological analysis using Cystoscope to obtain biological connected genes. Through literature searching and performing mRNA expression analysis on the other independent public dataset (GSE10799), we finally focused on TBX2. Functional effects of TBX2 were performed in tumorigenicity assays including migration and invasion assays, cell proliferation assay, and cell cycle assay. In addition, mechanically, we found enriched pathways related to bone metastasis using Gene Set Enrichment Analysis (GSEA) and validated our results by western blot.

Result: A total of 1132 significant genes were sorted initially. We selected common significant genes (log FC>2; p<0.01) from both the biological network data and microarray data. In total, 44 such genes were identified. we found TBX2, along with 10 other genes, to be reported with relevance to bone metastasis in other cancer types. Moreover, TBX2 showed significantly higher expression levels in patients that were found positive for metastasis to bone marrow compared to patients that did not exhibit this type of metastasis in the other separated cohort (GSE10799). Thus, we finally focused on TBX2. We found that TBX2 had detectable expression in LAC cell lines and silencing endogenous TBX2 expression in A549 and H1299 cell lines markedly suppressed migration and invasion, cell proliferation and arrested cell-cycle. Pathway enrichment analyses suggested that TBX2 drove LAC oncogenesis and metastasis through various pathways with epithelial mesenchymal transition (EMT) figuring prominently in the bone metastatic group, which was evidenced by western blot.

Conclusion: Collectively, TBX2 plays as a potential predictor of bone metastasis from LAC, yielding a better promise view towards “driver” gene responsible for bone metastasis.

Keywords: Integrated bioinformatics analysis, TBX2, Bone metastasis, Lung adenocarcinoma.