Abstract:Background: Lung cancer is the most commonly diagnosed cancer and the leading cause of cancer-related death worldwide. Moreover, it is highly susceptible to distant metastasis, which is the main cause of pain in advanced lung cancer, and frequently occurs in the bone. This study aimed to identify the differentially expressed genes (DEGs) related to metastatic bone disease in lung cancer using bioinformatics methods and to analyze the risk factors influencing the incidence of secondary bone metastasis in lung cancer. Methods: Gene expression profiles from the GSE175601 and GSE10799 datasets in the Gene Expression Omnibus (GEO) database were analyzed to screen for the DEGs associated with lung cancer bone metastasis. The STRING database was used to construct a protein-protein interaction (PPI) network, and the MCODE plugin was used to identify the key genes. The expression of these important genes in lung tumor tissues and their correlation with prognosis were validated in The Cancer Genome Atlas (TCGA) database. An examination of clinical data from patients diagnosed with stage IV lung adenocarcinoma treated at the Anhui No. 2 Provincial People's Hospital was conducted. Immunohistochemistry was used to examine the expression of key genes in lung cancer tumor tissues. A binary logistic regression analysis was conducted to examine the interactions in the expression of critical genes associated with bone metastasis in lung carcinoma patients. Results: In total, 59 DEGs were identified in the GSE175601 and GSE10799 datasets through Venn diagram construction. The PPI network analysis revealed two significant modules and eight candidate genes (LAPTM5, LCP2, CD53, ARHGAP25, C1QA, DES, MYH11, and VIM). According to TCGA database analysis, in carcinogenic tissues of the lung, the expression of these eight critical genes is downregulated. Further, only the lung cancer patients who had high expressions of ARHGAP25 had an improved progress-free interval (PFI) (P<0.05), disease-specific survival (DSS), and overall survival (OS). Of the 49 with stage IV lung adenocarcinoma patients included in the study, 27 (55.10%) developed bone metastasis. The immunohistochemical (IHC) results indicated that the expression score of ARHGAP25 was significantly lower in the group with bone metastasis (3.93±2.95) than the group without bone metastasis (6.64±3.62) (P=0.006). The proportion of patients with low ARHGAP25 expression was significantly higher in the group with bone metastasis (70.37%, 19/27) than the group without bone metastasis (31.82%, 7/22) (P=0.007). The binary logistic regression analysis identified serum alkaline phosphatase (ALP) and ARHGAP25 expression levels as independent risk factors for the occurrence of secondary bone metastatic disease in lung carcinoma patients. Conclusions: The key gene ARHGAP25 identified through bioinformatics for lung cancer bone metastasis was significantly downregulated. Its low expression constitutes an independent risk factor for secondary bone metastatic disease in patients with lung carcinoma.

Bioinformatics Analysis for the Identification of Key Genes and Long Non-Coding RNAs Related to Bone Metastasis in Breast Cancer.

The oncogenic potentials and diagnostic significance of long non-coding RNA LINC00310 in breast cancer.

Exploring the molecular mechanism associated with breast cancer bone metastasis using bioinformatic analysis and microarray genetic interaction network

New insight into long non-coding RNAs associated with bone metastasis of breast cancer based on an integrated analysis

Identification of Candidate Biomarkers Correlated With the Pathogenesis and Prognosis of Non-small Cell Lung Cancer via Integrated Bioinformatics Analysis

Identification of Key Genes and Pathways in Ewing's Sarcoma Patients Associated with Metastasis and Poor Prognosis

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

Identification of Key Genes Associated with Brain Metastasis from Breast Cancer: A Bioinformatics Analysis

Identification of Genes Associated with Breast Cancer Metastasis to Bone on A Protein-Protein Interaction Network with A Shortest Path Algorithm

Bioinformatics Analysis of Prognosis-Related Long Non-Coding RNAs in Invasive Breast Carcinoma

Identification of Hub Genes to Regulate Breast Cancer Spinal Metastases by Bioinformatics Analyses

Identification of Key Gene Signatures Associated With Bone Metastasis in Castration-Resistant Prostate Cancer Using Co-Expression Analysis

Identification of hub genes for early detection of bone metastasis in breast cancer

Bioinformatic analysis of differentially expressed genes in lung cancer bone metastasis and their implications for disease progression in lung cancer patients

A signal-based method for finding driver modules of breast cancer metastasis to the lung

Gene Expression Profile Analysis of the Bone Microenvironment in Patients with Spinal Metastases

Comprehensive Analysis of lncRNA and miRNA Regulatory Network Reveals Potential Prognostic Non-coding RNA Involved in Breast Cancer Progression

A gene expression signature-based nomogram model in prediction of breast cancer bone metastases

Development and Analysis of Long Non-Coding RNA-associated Competing Endogenous RNA Network for Osteosarcoma Metastasis

Identification of Potential Crucial Genes and Key Pathways in Breast Cancer Using Bioinformatic Analysis

Construction of Bone Metastasis-Specific Regulation Network Based on Prognostic Stemness-Related Signatures in Breast Invasive Carcinoma