Evaluation of the HOXA11 level in patients with lung squamous cancer and insights into potential molecular pathways via bioinformatics analysis

Role of homeobox d10 gene targeted signaling pathways in cancers

Homeobox D10 (HOXD10) is a transcription factor from the homeobox gene family that controls cell differentiation and morphogenesis throughout development.Due to their functional interaction, changes in HOXD10 gene expression might induce tumors. This narrative review focuses on how and why the dysregulation in the signaling pathways linked with HOXD10 contributes to the metastatic development of cancer. Organ development and tissue homeostasis need highly conserved homeotic transcription factors from homeobox (HOX) genes. Their dysregulation disrupts regulatory molecule action, causing tumors. The HOXD10 gene is upregulated in breast, gastric, hepatocellular, colorectal, bladder, cholangiocellular carcinoma and prostate cancer. Tumor signaling pathways are affected by HOXD10 gene expression changes. This study examines HOXD10-associated signaling pathway dysregulation, which may alter metastatic cancer signaling. In addition, the theoretical foundations that alter HOXD10-mediated therapeutic resistance in malignancies has been presented. New cancer therapy methods will be simpler to develop with the newly discovered knowledge. This review showed that HOXD10 may be a tumor suppressor gene and a new cancer treatment target signaling pathway.

Systematic Analysis of the Prognostic Significance and Roles of the Integrin Alpha Family in Non-Small Cell Lung Cancers

INTRODUCTION

Lung cancer is one of the most common cancer malignancies and the principal cause of cancer-associated deaths worldwide. Non-small cell lung cancers (NSCLCs) account for more than 80% of all lung cancer cases. Recent studies showed that the genes of the integrin alpha (α) (ITGA) subfamily play a fundamental role in various cancers. However, little is known about the expression and roles of distinct ITGA proteins in NSCLCs.

METHODS

Gene Expression Profiling Interactive Analysis and UALCAN (University of ALabama at Birmingham CANcer) web resources and The Cancer Genome Atlas (TCGA), ONCOMINE, cBioPortal, GeneMANIA, and Tumor Immune Estimation Resource databases were used to evaluate differential expression, correlations between the expression levels of individual genes, the prognostic value of overall survival (OS) and stage, genetic alterations, protein-protein interactions, and the immune cell infiltration of ITGAs in NSCLCs. We used R (v. 4.0.3) software to conduct gene correlation, gene enrichment, and clinical correlation of RNA sequencing data of 1016 NSCLCs from TCGA. To evaluate the expression of ITGA5/8/9/L at the expression and protein levels, qRT-PCR, immunohistochemistry (IHC), and hematoxylin and eosin (H&E) were performed, respectively.

RESULTS

Upregulated levels of ITGA11 messenger RNA and downregulated levels of ITGA1/3/5/7/8/9/L/M/X were observed in the NSCLC tissues. Lower expression of ITGA5/6/8/9/10/D/L was discovered to be expressively associated with advanced tumor stage or poor patient prognosis in patients with NSCLC. A high mutation rate (44%) of the ITGA family was observed in the NSCLCs. Gene Ontology functional enrichment analyses results revealed that the differentially expressed ITGAs could be involved in roles related to extracellular matrix (ECM) organization, collagen-containing ECM cellular components, and ECM structural constituent molecular functions. The results of the Kyoto Encyclopedia of Genes and Genomes analysis revealed that ITGAs may be involved in focal adhesion, ECM-receptor interaction, and amoebiasis; the expression of ITGAs was significantly correlated with the infiltration of diverse immune cells in NSCLCs. ITGA5/8/9/L was also highly correlated with PD-L1 expression. The validation results for marker gene expression in NSCLC tissues by qRT-PCR, IHC, and H&E staining indicated that the expression of ITGA5/8/9/L decreased compared with that in normal tissues.

CONCLUSION

As potential prognostic biomarkers in NSCLCs, ITGA5/8/9/L may fulfill important roles in regulating tumor progression and immune cell infiltration.

An insight into the diagnostic and prognostic value of HOX A13 ’s expression in non‐muscle invasive bladder cancer

Background

Several studies have interrogated the molecular pathways and their interacting genes underlying bladder cancer (BCa) tumorigenesis, yet, the role of homeobox genes is still poorly understood. Specifically, HOXA13, which plays an important role as a major actor in the urogenital tract's development.

Methods

Immunohistochemical (IHC) staining was performed to inspect the differential expression of HOXA13 protein in non‐muscle‐invasive bladder cancer (NMIBC) and non‐tumoral tissues. A semiquantitative scoring system was adopted to evaluate the IHC labeling. Correlation to clinical parameters was performed by descriptive statistics. Overall survival was estimated by the Kaplan–Meier method and Cox regression model. The functional HOX A13 protein association networks (PPI) were obtained using String 11.0 database.

Results

HOX A13 exhibited cytoplasmic and nuclear staining. Its expression levels were lower in high‐grade NMIBC (HG NMIBC) compared to low‐grade ones (LG NMIBC). The expression of HOX A13 was correlated to tumor grade (LG/HG) (p = 0.036) and stage (TA/T1) (p = 0.036). Nevertheless, its expression was not correlated to clinical parameters and was not able to predict the overall survival of patients with HG NMIBC. Finally, PPI analysis revealed that HOX A13 seems to be a part of a molecular network holding mainly PBX1, MEIS, ALDH1A2, HOX A10, and HOX A11.

Conclusion

The deregulation of HOX A13 is not associated with the prognosis of BCa. It seems to be rather implicated in the early initiation of urothelial tumorigenesis and thus may serve as a diagnostic marker in patients with NMIBC. Further experimentations on larger validation sets are mandatory.

One sample fits all: a microfluidic-assisted methodology for label-free isolation of CTCs with downstream methylation analysis of cfDNA in lung cancer

Lung cancer (LC) is a major cause of mortality. Late diagnosis, associated with limitations in tissue biopsies for adequate tumor characterization contribute to limited survival of lung cancer patients. Liquid biopsies have been introduced to improve tumor characetrization through the analysis of biomarkers, including circulating tumour cells (CTCs) and cell-free DNA (cfDNA). Considering their availability in blood, several enrichment strategies have been developed to augment circulating biomarkers for improving diagnostic, prognostic and treament efficacy assessment; often, however, only one biomarker is tested. In this work we developed and implemented a microfluidic chip for label-free enrichment of CTCs with a methodology for subsequent cfDNA analysis from the same cryopreserved sample. CTCs were successfully isolated in 38 of 42 LC patients with the microfluidic chip. CTCs frequency was significantly higher in LC patients with advanced disease. A cut-off of 1 CTC per mL was established for diagnosis (sensitivity = 76.19%, specificity = 100%) and in patients with late stage lung cancer, the presence of ≥5 CTCs per mL was significantly associated with shorter overall survival. MIR129-2me and ADCY4me panel of cfDNA methylation performed well for LC detection, whereas MIR129-2me combined with HOXA11me allowed for patient risk stratification. Analysis of combinations of biomarkers enabled the definition of panels for LC diagnosis and prognosis. Overall, this study demonstrates that multimodal analysis of tumour biomarkers via microfluidic devices may significantly improve LC characterization in cryopreserved samples, constituting a reliable source for continuous disease monitoring.

UCHL1 acts as a potential oncogene and affects sensitivity of common anti-tumor drugs in lung adenocarcinoma

Background

Lung adenocarcinoma is the leading cause of cancer death worldwide. Recently, ubiquitin C-terminal hydrolase L1 (UCHL1) has been demonstrated to be highly expressed in many tumors and plays the role of an oncogene. However, the functional mechanism of UCHL1 is unclear in lung adenocarcinoma progression.

Methods

We analyzed the differential expression of the UCHL1 gene in lung adenocarcinoma and normal lung tissues, and the correlation between the UCHL1 gene and prognosis was also analyzed by the bioinformatics database TCGA. Meanwhile, we detected and analyzed the expression of UCHL1 and Ki-67 protein in a tissue microarray (TMA) containing 150 patients with lung adenocarcinoma by immunohistochemistry (IHC) and clinicopathological characteristics by TCGA database. In vitro experiments, we knocked down the UCHL1 gene of A549 cells and detected the changes in cell migration, invasion, and apoptosis. At the same time, we analyzed the effect of UCHL1 on anti-tumor drug sensitivity of lung adenocarcinoma by a bioinformatics database. In terms of the detection rate of lung adenocarcinoma indicators, we analyzed the impact of UCHL1 combined with common clinical indicators on the detection rate of lung adenocarcinoma through a bioinformatics database.

Results

In this study, the analysis of UCHL1 protein expression in lung adenocarcinoma proved that obviously higher UCHL1 protein level was discovered in lung adenocarcinoma tissues. The expression of UCHL1 was closely related to poor clinical outcomes. Interestingly, a significantly positive correlation between the expression of UCHL1 and Ki-67-indicated UCHL1 was associated with tumor migration and invasion. Through executing loss of function tests, we affirmed that silencing of UCHL1 expression significantly inhibited migration and invasion of lung adenocarcinoma cells in vitro. Furthermore, lung adenocarcinoma cells with silenced UCHL1 showed a higher probability of apoptosis. In terms of the detection rate of lung adenocarcinoma indicators, we discovered UCHL1 could improve the detection rate of clinical lung adenocarcinoma and affect drug sensitivity.

Conclusion

In lung adenocarcinoma, UCHL1 promotes tumor migration, invasion, and metastasis by inhibiting apoptosis and has an important impact on the clinical drug treatment of lung adenocarcinoma. In addition, UCHL1 can improve the detection rate of clinical lung adenocarcinoma. Above all, UCHL1 may be a new marker for the diagnosis of lung adenocarcinoma and provide a new target for the treatment of clinical diseases.

Comprehensive analysis of angiogenesis subtype of squamous cell carcinoma

Background

Squamous cell carcinoma (SCC) is a disease with distinct management complexities as it displays a remarkably heterogeneous molecular subtype. However, the landscape of angiogenesis for SCC is not fully investigated.

Method and materials

The angiogenesis-related subtypes of SCC were established by using the ConsensusClusterPlus package based on angiogenesis-related genes and TCGA data. We analyzed the alteration of genes and miRNAs as well as pathways associated with angiogenesis subtypes. Next, the regulation network, the correlation with genomic characteristics, immune microenvironment, and clinical features of the angiogenesis subtypes were further investigated. Finally, the prognostic impact of the angiogenesis-related subtypes for SCC was also analyzed.

Results

A total of 1368 SCC samples were included in this study. Two angiogenesis subtypes were then identified based on the one hundred and sixty-three angiogenesis-related genes with subtype1 (angiogenesis subtype) of 951 SCC patients and subtype2 (non-angiogenesis subtype) of 417 SCC. GSEA revealed that angiogenesis and epithelial-mesenchymal transition, inflammatory response, and hypoxia were enriched in the angiogenesis subtype. Eight of the 15 immune checkpoints (ADORA2A, BTLA, CD276, CYBB, HAVCR2, SIGLEC7, SIGLEC9, and VTCN1) were significantly upregulated while C10orf54 were significantly downregulated in the angiogenesis subtype. The survival analysis revealed that the patients in the angiogenesis subtype have poorer survival outcomes than those in the non-angiogenesis subtype (P = 0.017 for disease-free interval and P = 0.00013 for overall survival).

Conclusion

Our analysis revealed a novel angiogenesis subtype classification in SCC and provides new insights into a hallmark of SCC progression.

Immune and Metabolic Dysregulated Coding and Non-coding RNAs Reveal Survival Association in Uterine Corpus Endometrial Carcinoma

Uterine corpus endometrial carcinoma (UCEC) is one of the most common gynecologic malignancies, but only a few biomarkers have been proven to be effective in clinical practice. Previous studies have demonstrated the important roles of non-coding RNAs (ncRNAs) in diagnosis, prognosis, and therapy selection in UCEC and suggested the significance of integrating molecules at different levels for interpreting the underlying molecular mechanism. In this study, we collected transcriptome data, including long non-coding RNAs (lncRNAs), microRNAs (miRNAs), and messenger RNAs (mRNAs), of 570 samples, which were comprised of 537 UCEC samples and 33 normal samples. First, differentially expressed lncRNAs, miRNAs, and mRNAs, which distinguished invasive carcinoma samples from normal samples, were identified, and further analysis showed that cancer- and metabolism-related functions were enriched by these RNAs. Next, an integrated, dysregulated, and scale-free biological network consisting of differentially expressed lncRNAs, miRNAs, and mRNAs was constructed. Protein-coding and ncRNA genes in this network showed potential immune and metabolic functions. A further analysis revealed two clinic-related modules that showed a close correlation with metabolic and immune functions. RNAs in the two modules were functionally validated to be associated with UCEC. The findings of this study demonstrate an important clinical application for improving outcome prediction for UCEC.

MiR-30b-5p Influences Chronic Exercise Arthritic Injury by Targeting Hoxa1

We identified the role of miR-30b-5p in chronic exercise arthritic injury. Rats with chronic exercise arthritic injury received treatment with miR-30b-5p antagomiR. H&E and Safranin O-fast green staining were performed. The levels of tumor necrosis factor-α (TNF-α) and interleukin-1β (IL-1β) were detected. The binding relationship between homeobox A1 (Hoxa1) and miR-30b-5p was revealed. After manipulating the expressions of miR-30b-5p and/or Hoxa1 in chondrocytes, the viability, apoptosis and migration of chondrocytes were assessed. The levels of molecules were determined by qRT-PCR or Western blot. MiR-30b-5p antagomiR ameliorated articular cartilage lesion and destruction, reduced Mankin's score and the levels of TNF-α, IL-1β, miR-30b-5p, matrix metallopeptidase 13 (MMP-13), and cleaved caspase-3, and increased relative thickness and the levels of Hoxa1, Aggrecan and type II collagen (COLII) in model rats. MiR-30b-5p up-regulation decreased Hoxa1 level, viability, migration and induced apoptosis, whereas miR-30b-5p down-regulation produced the opposite effects. MiR-30b-5p up-regulation increased the levels of MMP-13 and cleaved caspase-3, but decreased those of Aggrecan and COLII in chondrocytes. However, the action of miR-30b-5p up-regulation on chondrocytes was reversed by Hoxa1 overexpression. In conclusion, miR-30b-5p is involved in cartilage degradation in rats with chronic exercise arthritic injury and regulates chondrocyte apoptosis and migration by targeting Hoxa1.

Pan-Cancer Genome-Wide DNA Methylation Analyses Revealed That Hypermethylation Influences 3D Architecture and Gene Expression Dysregulation in HOXA Locus During Carcinogenesis of Cancers

DNA methylation dysregulation during carcinogenesis has been widely discussed in recent years. However, the pan-cancer DNA methylation biomarkers and corresponding biological mechanisms were seldom investigated. We identified differentially methylated sites and regions from 5,056 The Cancer Genome Atlas (TCGA) samples across 10 cancer types and then validated the findings using 48 manually annotated datasets consisting of 3,394 samples across nine cancer types from Gene Expression Omnibus (GEO). All samples’ DNA methylation profile was evaluated with Illumina 450K microarray to narrow down the batch effect. Nine regions were identified as commonly differentially methylated regions across cancers in TCGA and GEO cohorts. Among these regions, a DNA fragment consisting of ∼1,400 bp detected inside the HOXA locus instead of the boundary may relate to the co-expression attenuation of genes inside the locus during carcinogenesis. We further analyzed the 3D DNA interaction profile by the publicly accessible Hi-C database. Consistently, the HOXA locus in normal cell lines compromised isolated topological domains while merging to the domain nearby in cancer cell lines. In conclusion, the dysregulation of the HOXA locus provides a novel insight into pan-cancer carcinogenesis.

Expression Profile and Prognostic Values of HOXA Family Members in Laryngeal Squamous Cell Cancer

The homeobox A cluster (HOXA) gene family, comprising 11 members, is involved in a wide spectrum of biological functions in human cancers. However, there is little research on the expression profile and prognostic values of HOXA genes in laryngeal squamous cell cancer (LSCC). Based on updated public resources and integrative bioinformatics analysis, we assessed the expression profile and prognostic values of the HOXA family members. Expression and methylation data on HOXA family members were obtained from The Cancer Genome Atlas (TCGA). The prognostic values of HOXA members and clinical features were identified. A gene set enrichment analysis (GSEA) was conducted to explore the mechanism underlying the involvement of HOXA members in LSCC. The associations between tumor immune infiltrating cells (TIICs) and the HOXA family members were evaluated using the Tumor Immune Estimation Resource (TIMER) database. HOXA2 and HOXA4 were downregulated and HOXA7 and HOXA9–13 were upregulated in LSCC. Upregulation of HOXA10, HOXA11, and HOXA13, along with two clinical characteristics (M stage and gender), were associated with a poor LSCC prognosis based on the results of univariate and multivariate Cox proportional hazards regression analyses. Although there were no significant correlations between TIICs and HOXA members, the GSEA results indicated that HOXA members participate in multiple biological processes underlying tumorigenesis. This study comprehensively analyzed the HOXA members, providing insights for further investigation of the HOXA family members as potential targets in LSCC.

A self-enforcing HOXA11/Stat3 feedback loop promotes stemness properties and peritoneal metastasis in gastric cancer cells

Rationale: Peritoneal metastasis is one of the most common and life-threatening metastases in gastric cancer patients. The disseminated gastric cancer cells forming peritoneal metastasis exhibit a variety of characteristics that contrast with those of adjacent epithelial cell of gastric mucosa and even primary gastric cancer cells. We hypothesized that the gene expression profiles of peritoneal foci could reveal the identities of genes that might function as metastatic activator.

Methods: In this study, we show, using in vitro, in vivo, in silico and gastric cancer tissues studies in humans and mice, that Homoebox A11 (HOXA11) potently promote peritoneal metastasis of gastric cancer cells.

Results: Its mechanism of action involves alternation of cancer stemness and subsequently enhancement of the adhesion, migration and invasion and anti-apoptosis. This is achieved, mainly, through formation of a positive feedback loop between HOXA11 and Stat3, which is involved in the stimulation of Stat3 signaling pathway.

Conclusions: These observations uncover a novel peritoneal metastatic activator and demonstrate the association between HOXA11, Stat3 and cancer stemness of gastric cancer cells, thereby revealing a previously undescribed mechanism of peritoneal metastasis.

Diagnostic and prognostic values of integrin α subfamily mRNA expression in colon adenocarcinoma

The integrin α (ITGA) subfamily genes play a fundamental role in various cancers. However, the potential mechanism and application values of ITGA genes in colon adenocarcinoma (COAD) remain elusive. The present study investigated the significance of the expression of ITGA genes in COAD from the perspective of diagnosis and prognosis. A COAD RNA-sequencing dataset was obtained from The Cancer Genome Atlas. The present study investigated the biological function of the ITGA subfamily genes through bioinformatics analysis. Reverse transcription-quantitative polymerase chain reaction was applied to investigate the distribution of integrin α8 (ITGA8) expression in COAD tumors and adjacent normal tissues. Bioinformatics analysis indicated that ITGA genes were noticeably enriched in cell adhesion and the integrin-mediated signaling pathway, and co-expressed with each other. It was also revealed through observation that the majority of gene expression was significantly low in tumor tissues (P<0.05), and diagnostic receiver operating characteristic curves revealed that most of the genes could serve as significant diagnostic markers in COAD (P<0.05), especially ITGA8 which had a high diagnostic value with an area under curve (AUC) of 0.989 [95% confidence interval (CI) 0.980–0.997] in COAD (P<0.0001). In addition, ITGA8 expression was verified in clinical samples and it was revealed that it was higher in adjacent normal tissues (P=0.041) compared to COAD tissues, and the AUC was 0.704 (95% CI, 0.577–0.831; P<0.0085). Multivariate survival analysis indicated that integrin α (ITGA5) may be an independent prognostic indicator for COAD overall survival. Gene set enrichment analysis indicated that ITGA5 may participate in multiple biological processes and pathways. The present study revealed that ITGA genes were associated with the diagnosis and prognosis of COAD. The mRNA expression of ITGA8 may be a potential diagnosis biomarker and ITGA5 may serve as an independent prognosis indicator for COAD.

Protective potential of miR-146a-5p and its underlying molecular mechanism in diverse cancers: a comprehensive meta-analysis and bioinformatics analysis

Background/aims

Studies have shown that miR-146a-5p was differentially expressed in diverse cancers, but the associations between miR-146a-5p expression and prognosis across multiple types of cancer as well its potential targets and downstream pathways have not been comprehensively analyzed. In this study, we performed the first meta-analysis of the prognostic value of miR-146a-5p expression in diverse malignancies and explored prospective targets of miR-146a-5p and related signaling pathways.

Methods

A thorough search for articles related to miR-146a-5p was performed, and RNA-seq data from The Cancer Genome Atlas (TCGA) and microarray data from gene expression omnibus profiles were used to collect information about the prognostic value of miR-146a-5p. A comprehensive meta-analysis was conducted. Twelve platforms in miRWalk 2.0 were applied to predict targets of miR-146a-5p. TCGA RNA-seq data were used to validate the inverse relationships between miR-146a-5p and its likely targets. Subsequently, gene ontology and pathway analyses were conducted using Funrich version 3.1.3. Potential protein–protein interaction (PPI) networks were constructed. Potential target genes of miR-146a-5p in lung cancer were validated by RT-qPCR.

Results

We included 10 articles in the meta-analysis. In a pooled analysis, the high miR-146a-5p expression group showed a better overall survival in solid cancers, particularly in reproductive system cancers and digestive system cancers. A total of 120 predicted target genes were included in a bioinformatics analysis. Five pathways involving phospholipase C (PLC) and aquaporins (AQPs) were the most significantly enriched Kyoto Encyclopedia of Genes and Genomes pathways. Moreover, the PPI network displayed the related signaling pathways and interactions among proteins. AQP1 and FYN were validated by RT-qPCR to be potential targets of miR-146a-5p in lung cancer.

Conclusion

There is a close link between high miR-146a-5p expression and better overall survival in 21 types of solid cancer, especially in reproductive system and digestive system cancers. Furthermore, miR-146a-5p could inhibit diverse malignancies by modulating pathways linked to PLC or AQPs. In summary, miR-146a-5p is a potential prognostic biomarker and therapeutic target for various cancers.

An inverse interaction between HOXA11 and HOXA11-AS is associated with cisplatin resistance in lung adenocarcinoma

HOXA11, which is a member of the homeobox (HOX) gene family, and its natural antisense transcript (NAT) HOXA11-AS have been reported to be closely related to the development of lung cancer. We aimed to investigate their specific roles in cisplatin (DDP) resistance in lung adenocarcinoma (LUAD). First, we found that HOXA11 is hypermethylated and significantly downregulated in a DDP-resistant A549 cell line (A549/DDP) and LUAD tissues, while the HOXA11-AS expression level is elevated. Although HOXA11 and HOXA11-AS mRNA overlap in the 5'-untranslated region (5' UTR) and share two CpG islands, DNA methylation only regulates the expression of HOXA11. Then, we found that HOXA11 and HOXA11-AS have an inverse interaction by transfecting their siRNAs and overexpression vectors into A549 and A549/DDP cells. A dual-luciferase reporter assay further confirmed that the overlapping 5'UTR is essential for the bidirectional regulation between HOXA11 and HOXA11-AS. Functional analysis showed that knockdown of HOXA11 expression in A549 cells induced DDP resistance and activated Akt/β-catenin signaling, while overexpression of HOXA11 in A549/DDP cells increased DDP sensitivity and inhibited Akt/β-catenin signaling. Moreover, HOXA11-AS knockdown in A549 cells increased DDP sensitivity and inhibited Akt/β-catenin signaling, while the overexpression of HOXA11-AS in A549/DDP cells induced DDP resistance and activated Akt/β-catenin signaling. In conclusion, our study demonstrates that the inverse interaction between HOXA11 and HOXA11-AS promotes DDP resistance in LUAD.

Deciphering The Potential Role of Hox Genes in Pancreatic Cancer

The Hox gene family plays an important role in organogenesis and animal development. Currently, 39 Hox genes that are clustered in four chromosome regions have been identified in humans. Emerging evidence suggests that Hox genes are involved in the development of the pancreas. However, the expression of Hox genes in pancreatic tumor tissues has been investigated in only a few studies. In addition, whether specific Hox genes can promote or suppress cancer metastasis is not clear. In this article, we first review the recent progress in studies on the role of Hox genes in pancreatic cancer. By comparing the expression profiles of pancreatic cancer cells isolated from genetically engineered mice established in our laboratory with three different proliferative and metastatic abilities, we identified novel Hox genes that exhibited tumor-promoting activity in pancreatic cancer. Finally, a potential oncogenic mechanism of the Hox genes was hypothesized.

Comprehensive clinical implications of homeobox A10 in 3,199 cases of non-small cell lung cancer tissue samples combining qRT-PCR, RNA sequencing and microarray data

In the current study, we proposed to explore the potential role and related signaling pathways of Homobox A10 (HOXA10) in non-small cell lung cancer (NSCLC). HOXA10 levels in lung adenocarcinoma (LUAD) and lung squamous cell carcinoma (LUSC) were detected by qRT-PCR and the expression of HOXA10 was significantly up-regulated in the NSCLC tissue of all 55 pairs (P = 0.037). Overexpression of HOXA10 was closely correlated with the clinical stage of LUSC (P = 0.011). HOXA10 expression in RNA sequencing data based on 1, 077 cases exhibited concordant significant up-regulation in NSCLC, LUAD and LUSC (P < 0.001). In NSCLC, HOXA10 expression was closely correlated to patient T stage (P = 0.006). In LUAD, HOXA10 expression was compactly correlated to patient N stage (P = 0.02). Some of the microarrays from Gene Expression Omnibus (GEO) and ArrayExpress showed consistent over-expression of HOXA10 levels in NSCLCs. More importantly, the combined SMD value was 0.052 (95% CI: 0.29-0.75, P < 0.001) generated by meta-analysis from 47 datasets based on 4, 616 cases of NSCLC. The area under the curve (AUC) of SROC supported the over-expression of HOXA10 in NSCLC as being 0.88 (95% CI: 0.81-0.93), with sensitivity and specificity of 0.88 (95% CI: 0.81-0.93) and 0.56 (95% CI: 0.44-0.66), respectively. In addition, 111 co-expressed genes were collected from cBioPortal and enriched in "cell cycle", "cell adhesion molecules", "p53 signaling", and "adherens junction". Interestingly, an up-regulation trend of HOXA10 protein expression was also observed in NSCLC through tissue chips and immunohistochemistry. In conclusion, the overexpression of HOXA10 may play a pivotal role in the tumorigenesis of NSCLC, and this effect is observed more obviously in LUSC than in LUAD.