Downregulation of HOXA3 in lung adenocarcinoma and its relevant molecular mechanism analysed by RT-qPCR, TCGA and in silico analysis

Role of HOXA1-4 in the development of genetic and malignant diseases

The HOXA genes, belonging to the HOX family, encompass 11 members (HOXA1-11) and exert critical functions in early embryonic development, as well as various adult processes. Furthermore, dysregulation of HOXA genes is implicated in genetic diseases, heart disease, and various cancers. In this comprehensive overview, we primarily focused on the HOXA1-4 genes and their associated functions and diseases. Emphasis was placed on elucidating the impact of abnormal expression of these genes and highlighting their significance in maintaining optimal health and their involvement in the development of genetic and malignant diseases. Furthermore, we delved into their regulatory mechanisms, functional roles, and underlying biology and explored the therapeutic potential of targeting HOXA1-4 genes for the treatment of malignancies. Additionally, we explored the utility of HOXA1-4 genes as biomarkers for monitoring cancer recurrence and metastasis.

Identification Invasion-Related Long Non-Coding RNAs in Lung Adenocarcinoma and Analysis of Competitive Endogenous RNA Regulatory Networks

Background

Cell invasion plays a vital role in cancer development and progression. Aberrant expression of long non-coding RNAs (lncRNAs) is also critical in carcinogenesis. However, the prognostic value of invasion-related lncRNAs in lung adenocarcinoma (LUAD) remains unknown.

Methods

Differentially expressed mRNAs (DEmRNAs), lncRNAs (DElncRNAs), and microRNAs (DEmiRNAs) were between LUAD and control samples. Pearson correlation analyses were performed to screen for invasion-related DElncRNAs (DEIRLs). Univariate and multivariate Cox regression algorithms were applied to identify key genes and construct the risk score model, which was evaluated using receiver operating characteristic (ROC) curves. Gene set enrichment analysis (GSEA) was used to explore the underlying pathways of the risk model. Moreover, an invasion-related competitive endogenous RNA (ceRNA) regulatory network was constructed. Reverse transcription-quantitative polymerase chain reaction (RT-qPCR) was performed to detect the expression of prognostic lncRNAs in the LUAD and control samples.

Results

A total of 45 DElncRNAs were identified as DEIRLs. RP3-525N10.2, LINC00857, EP300-AS1, PDZRN3-AS1, and RP5-1102E8.3 were potential prognostic lncRNAs, the expression of which was verified by RT-qPCR in LUAD samples. Both the risk score model and nomogram used the prognostic lncRNAs. ROC curves showed the risk score model had moderate accuracy and the nomogram had high accuracy in predicting patient prognosis. GSEA results indicated that the risk score model was associated with many biological processes and pathways relevant to cell proliferation. A ceRNA regulatory network was constructed in which PDZRN3-miR-96-5p-CPEB1, EP300-AS1-miR-93-5p-CORO2B, and RP3-525N10.2-miR-130a-5p-GHR may be key invasion-related regulatory pathways in LUAD.

Conclusion

Our study identified five novel invasion-related prognostic lncRNAs (RP3-525N10.2, LINC00857, EP300-AS1, PDZRN3-AS1, and RP5-1102E8.3) and established an accurate model for predicting the prognosis of patients with LUAD. These findings enrich our understanding of the relationships between cell invasion, lncRNAs, and LUAD and may provide novel treatment directions.

The Profile and Clinical Significance of ITGB2 Expression in Non-Small-Cell Lung Cancer

Integrins are involved in extracellular and intracellular signaling and are often aberrantly expressed in tumors. Integrin beta 2 (ITGB2) has previously been demonstrated to be correlated with the host defense. However, the expression profile and role of ITGB2 in non-small-cell lung cancer (NSCLC) remain unclear. Here, we found that the genetic alterations in ITGB2 was predominated by gene mutation and copy number deletion using cBioPortal analysis, and its expression was downregulated in the NSCLC tissues, as validated by the UALCAN, TCGA, and GEO databases and our tissue samples. Kaplan-Meier (KM) plotter analysis revealed that patients with a lower ITGB2 expression had a shorter overall survival (OS) time (p = 0.01). Moreover, 1089 differentially expressed genes (DEGs) in the NSCLC tissues were screened using the TCGA database. The GO and KEGG enrichment analysis showed that the DEGs were closely associated with immune processes and cell adhesion. The protein-protein interaction (PPI) network revealed that 10 of 15 EMT-related genes among the DEGs might lead to the metastasis of NSCLC. Concomitantly, the expression of ITGB2 was positively correlated with the infiltration of Treg cells and Myeloid-derived suppressor cells (MDSC). Biologically, the ectopic expression of ITGB2 significantly inhibited the proliferation and metastasis of NSCLC cells. Mechanistically, we demonstrated that ITGB2 suppressed the expression of N-cadherin, Vimentin, Slug, Snail, and Twist, while it promoted E-cadherin expression, according to gain-of-function studies. In conclusion, ITGB2 can inhibit the proliferation and migration of NSCLC cells, leading to a poor prognosis, via epithelial-mesenchymal transition (EMT) signaling.

Expression Landscape and Functional Roles of HOXA4 and HOXA5 in Lung Adenocarcinoma

BACKGROUND

The role of HOXA family genes in the occurrence and progression of a variety of human cancers has been scatteredly reported. However, there is no systematic study on the differential expression, prognostic significance and potential molecular mechanism of HOXA4 and HOXA5 in LUAD.

METHODS

In-house immunohistochemistry (IHC), multi-center microarrays, RT-qPCR and RNA-seq data were incorporated for comprehensively evaluating the expression and prognostic value of HOXA4 and HOXA5 in LUAD. The mechanism of HOXA4 and HOXA5 in the formation and development of LUAD was analyzed from multiple aspects of immune correlations, upstream transcriptional regulation, functional states of single cells and co-expressed gene network. The functional roles of HOXA4 and HOXA5 in LUAD were validated by in vitro experiments.

RESULTS

As a result, in 3201 LUAD samples and 2494 non-cancer lung samples, HOXA4 and HOXA5 were significantly downexpressed (P < 0.05). The aberrant expression of HOXA5 was significantly correlated with the clinical progression of LUAD (P < 0.05). HOXA5 showed remarkable prognostic value for LUAD patients (P < 0.05). The expression of HOXA4 and HOXA5 in LUAD were negatively correlated with tumor purity and positively correlated with the infiltration of various immune cells such as B cells, T cells and macrophages. HOXA4 and HOXA5 overexpression had notable inhibitory effect on the proliferation, migration and invasion of LUAD cells.

CONCLUSIONS

In conclusion, the identified downexpressed HOXA4 and HOXA5 had significant distinguishing ability for LUAD samples and affected the cellular functions of LUAD cells. The low expression of HOXA5 indicated worse overall survival of LUAD patients. Therefore, the two HOXA family genes especially HOXA5 may serve as potential biomarkers for LUAD.

Molecular implications of HOX genes targeting multiple signaling pathways in cancer

Homeobox (HOX) genes encode highly conserved homeotic transcription factors that play a crucial role in organogenesis and tissue homeostasis. Their deregulation impacts the function of several regulatory molecules contributing to tumor initiation and progression. A functional bridge exists between altered gene expression of individual HOX genes and tumorigenesis. This review focuses on how deregulation in the HOX-associated signaling pathways contributes to the metastatic progression in cancer. We discuss their functional significance, clinical implications and ascertain their role as a diagnostic and prognostic biomarker in the various cancer types. Besides, the mechanism of understanding the theoretical underpinning that affects HOX-mediated therapy resistance in cancers has been outlined. The knowledge gained shall pave the way for newer insights into the treatment of cancer.

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.

The Identification and Validation of Two Heterogenous Subtypes and a Risk Signature Based on Ferroptosis in Hepatocellular Carcinoma

BACKGROUND

Ferroptosis is essential for tumorigenesis and progression of hepatocellular carcinoma (HCC). The heterogeneity of ferroptosis and its relationship with tumor microenvironment (TME) have still remain elusive.

METHODS

Based on 74 ferroptosis related genes (FRGs) and 3,933 HCC samples from 32 datasets, we comprehensively explored the heterogenous ferroptosis subtypes. The clinical significance, functional status, immune infiltration, immune escape mechanisms, and genomic alterations of different subtypes were further investigated.

RESULTS

We identified and validated two heterogeneous ferroptosis subtypes: C1 was metabolismlowimmunityhigh subtype and C2 was metabolismhighimmunitylow subtype. Compared to C2, C1 owned worse prognosis, and C1 tended to occur in the patients with clinical characteristics such as younger, female, advanced stage, higher grade, vascular invasion. C1 and C2 were more sensitive to immunotherapy and sorafenib, respectively. The immune escape mechanisms of C1 might be accumulating more immunosuppressive cells, inhibitory cytokines, and immune checkpoints, while C2 was mainly associated with inferior immunogenicity, defecting in antigen presentation, and lacking leukocytes. In addition, C1 was characterized by BAP1 mutation, MYC amplification, and SCD1 methylation, while C2 was characterized by the significant alterations in cell cycle and chromatin remodeling processes. We also constructed and validated a robust and promising signature termed ferroptosis related risk score (FRRS) for assessing prognosis and immunotherapy.

CONCLUSION

We identified and validated two heterogeneous ferroptosis subtypes and a reliable risk signature which used to assess prognosis and immunotherapy. Our results facilitated the understood of ferroptosis as well as clinical management and precise therapy of HCC.

The Significance of Secreted Phosphoprotein 1 in Multiple Human Cancers

Malignant tumor represents a major reason for death in the world and its incidence is growing rapidly. Developing the tools for early diagnosis is possibly a promising way to offer diverse therapeutic options and promote the survival chance. Secreted phosphoprotein 1 (SPP1), also called Osteopontin (OPN), has been demonstrated overexpressed in many cancers. However, the specific role of SPP1 in prognosis, gene mutations, and changes in gene and miRNA expression in human cancers is unclear. In this report, we found SPP1 expression was higher in most of the human cancers. Based on Kaplan-Meier plotter and the PrognoScan database, we found high SPP1 expression was significantly correlated with poor survival in various cancers. Using a large dataset of colon adenocarcinoma (COAD), head and neck cancer (HNSC), lung adenocarcinoma (LUAD), and lung squamous cell carcinoma (LUSC) patients from the Gene Expression Omnibus (GEO) and The Cancer Genome Atlas (TCGA) databases, this study identified 22 common genes and 2 common miRNAs. GO, and KEGG paths analyses suggested that SPP1 correlated genes were mainly involved in positive regulation of immune cell activation and infiltration. SPP1-associated genes and miRNAs regulatory networks suggested that their interactions may play a role in the progression of four selected cancers. SPP1 showed significant positive correlation with the immunocyte and immune marker sets infiltrating degrees. All of these data provide strong evidence that SPP1 may promote tumor progress through interacting with carcinogenic genes and facilitating immune cells’ infiltration in COAD, HNSC, LUAD, and LUSC.

Downregulation of miR‑193a‑3p via targeting cyclin D1 in thyroid cancer

Thyroid cancer (TC) is a frequently occurring malignant tumor with a rising steadily incidence. microRNA (miRNA/miR)‑193a‑3p is an miRNA that is associated with tumors, playing a crucial role in the genesis and progression of various cancers. However, the expression levels of miR‑193a‑3p and its molecular mechanisms in TC remain to be elucidated. The present study aimed to probe the expression of miR‑193a‑3p and its clinical significance in TC, including its underlying molecular mechanisms. Microarray and RNA sequencing data gathered from three major databases, specifically Gene Expression Omnibus (GEO), ArrayExpress and The Cancer Genome Atlas (TCGA) databases, and the relevant data from the literature were used to examine miR‑193a‑3p expression. Meta‑analysis was also conducted to evaluate the association between clinicopathological parameters and miR‑193a‑3p in 510 TC and 59 normal samples from the TCGA database. miRWalk 3.0, and the TCGA and GEO databases were used to predict the candidate target genes of miR‑193a‑3p. Gene Ontology, Kyoto Encyclopedia of Genes and Genomes and protein‑protein interaction network enrichment analyses were conducted by using the predicted candidate target genes to investigate the underlying carcinogenic mechanisms. A dual luciferase assay was performed to validate the targeting regulatory association between the most important hub gene cyclin D1 (CCND1) and miR‑193a‑3p. miR‑193a‑3p expression was considerably downregulated in TC compared with in the non‑cancer controls (P<0.001). The area under the curve of the summary receiver operating characteristic was 0.80. Downregulation of miR‑193a‑3p was also significantly associated with age, sex and metastasis (P=0.020, 0.044 and 0.048, respectively). Bioinformatics analysis indicated that a low miR‑193a‑3p expression may augment CCND1 expression to affect the biological processes of TC. In addition, CCND1, as a straightforward target, was validated through a dual luciferase assay. miR‑193a‑3p and CCND1 may serve as prognostic biomarkers of TC. Finally, miR‑193a‑3p may possess a crucial role in the genesis and progression of TC by altering the CCND1 expression.

Methylation in HOX Clusters and Its Applications in Cancer Therapy

HOX genes are commonly known for their role in embryonic development, defining the positional identity of most structures along the anterior–posterior axis. In postembryonic life, HOX gene aberrant expression can affect several processes involved in tumorigenesis such as proliferation, apoptosis, migration and invasion. Epigenetic modifications are implicated in gene expression deregulation, and it is accepted that methylation events affecting HOX gene expression play crucial roles in tumorigenesis. In fact, specific methylation profiles in the HOX gene sequence or in HOX-associated histones are recognized as potential biomarkers in several cancers, helping in the prediction of disease outcomes and adding information for decisions regarding the patient’s treatment. The methylation of some HOX genes can be associated with chemotherapy resistance, and its identification may suggest the use of other treatment options. The use of epigenetic drugs affecting generalized or specific DNA methylation profiles, an approach that now deserves much attention, seems likely to be a promising weapon in cancer therapy in the near future. In this review, we summarize these topics, focusing particularly on how the regulation of epigenetic processes may be used in cancer therapy.

Upregulated expression of SAC3D1 is associated with progression in gastric cancer

SAC3 domain containing 1 (SAC3D1) has been reported to be involved in numerous types of cancer. However, the role of SAC3D1 in GC has not yet been elucidated. In the present study, the mRNA expression level of SAC3D1 between GC and normal tissues were assessed with a continuous variable meta-analysis based on multiple datasets from public databases. The protein expression level of SAC3D1 in GC and normal tissues was assessed by an in-house immunohistochem-istry (IHC). The association between SAC3D1 expression and some clinical parameters was assessed based on the TCGA and IHC data. Survival analysis was performed to assess the association between SAC3D1 expression and the survival of GC patients. The co-expressed genes of SAC3D1 were determined by integrating three online tools, and the enrichment analyses were performed to determine SAC3D1-related pathways and hub co-expressed genes. SAC3D1 was significantly upregulated in GC tumor tissues in comparison to normal tissues with the SMD being 0.45 (0.12, 0.79). The IHC results also indicated that SAC3D1 protein expression in GC tissues was markedly higher than in normal tissues. The SMD following the addition of the IHC data was 0.59 (0.11, 1.07). The protein levels of SAC3D1 were positively associated with the histological grade, T stage and N stage of GC (P<0.001). The TCGA data also revealed that the SAC3D1 mRNA level was significantly associated with the N stage (P<0.001). Moreover, prognosis analysis indicated that SAC3D1 was closely associated with the prognosis of patients with GC. Moreover, 410 co-expressed genes of SAC3D1 were determined, and these genes were mainly enriched in the cell cycle. In total, 4 genes (CDK1, CCNB1, CCNB2 and CDC20) were considered key co-expressed genes. On the whole, these findings demonstrate that SAC3D1 is highly expressed in GC and may be associated with the progression of GC.

Upregulation of miR‑132‑3p in cholangiocarcinoma tissues: A study based on RT‑qPCR, The Cancer Genome Atlas miRNA sequencing, Gene Expression Omnibus microarray data and bioinformatics analyses

MicroRNAs (miRNAs/miRs) have been reported to be closely associated with numerous human diseases, including cholangiocarcinoma (CCA). However, the number of miRNAs known to be involved in CCA is limited, and the association between miR-132-3p and CCA remains unknown. In the present study, the clinical role of miR-132-3p and its potential signaling pathways were investigated by multiple approaches. Reverse transcription-quantitative PCR (RT-qPCR), CCA-associated Gene Expression Omnibus (GEO), ArrayExpress and Sequence Read Archive (SRA) miRNA-microarray or miRNA-sequencing data were screened, and meta-analyses were conducted, in order to calculate the receiver operating characteristic (ROC) curve and standardized mean difference (SMD). The predicted target genes of miR-132-3p were obtained from 12 online databases and were combined with the downregulated differentially expressed genes identified in the RNA-sequencing data of CCA. Gene Ontology annotation and pathway analysis were performed in WebGestalt. Protein-protein interaction analyses were conducted in STRING. The Cancer Genome Atlas (TCGA) mRNA expression profiles were used to validate the expression levels of hub genes at the mRNA level. The Human Protein Atlas was used to identify the protein expression levels of hub genes in CCA tissues and non-tumor biliary epithelium. The meta-analyses comprised 10 groups of RT-qPCR data, eight GEO microarray datasets and one TCGA miRNA-sequencing dataset. The SMD of miR-132-3p in CCA was 0.75 (95% CI: 0.25, 1.24), which indicated that miR-132-3p was overexpressed in CCA tissues. This finding was supported by a summary ROC value of 0.80 (95% CI: 0.76, 0.83). The pooled sensitivity and specificity were 0.81 (95% CI: 0.59, 0.93) and 0.71 (95% CI: 0.58, 0.81), respectively. The relative expression level of miR-132-3p in the early stage of CCA (stages I–II) was 6.8754±0.5279, which was markedly lower than that in the advanced stage (stages III–IVB), 7.3034±0.3267 (P=0.003). Consistently, the miR-132-3p level in low-grade CCA (grades G1-G2) was 6.7581±0.5297, whereas it was 7.1191±0.4651 in patients with high-grade CCA (grades G3-G4) (P=0.037). Furthermore, 555 potential target genes of miR-132-3p in CCA were mainly enriched in the ‘Focal Adhesion-PI3K-Akt-mTOR-signaling pathway’. In conclusion, upregulation of miR-132-3p may serve a pivotal role in the tumorigenesis and progression of CCA by targeting different pathways. Further in vitro and in vivo studies are required to support the current findings.

Hsa_circ_0008945 promoted breast cancer progression by targeting miR-338-3p

PURPOSE

To detect the expression and function of circ_0008945 in breast cancer (BC) and to explore its potential molecular mechanisms in BC tumorigenesis.

MATERIALS AND METHODS

We measured expression levels of circ_0008945, miR-338-3p and homeobox A3 (HOXA3) in BC tissue specimens and cells using quantitative reverse transcriptase polymerase chain reaction (qRT-PCR). We examined the effects of all three genes on BC cell proliferation using Cell Counting Kit-8 (CCK-8) and colony formation assays. We also performed a Transwell assay to assess the migratory and invasive ability of treated BC cells. BC cell apoptosis was assessed using flow cytometric (FCM) analysis; interaction between miR-338-3p and circ_0008945 or HOXA3 was verified by dual-luciferase reporter assay as well as by ribonucleic-acid (RNA) pulldown. Finally, we used an in vivo tumor growth assay to assess the role of circ_0008945 overexpression in BC tumor growth.

RESULTS

We found that circ_0008945 expression was significantly increased in both BC tissue specimens and cells. This increase was correlated with poor prognosis in BC patients. Knockdown of circ_0008945 inhibited BC cell proliferation, migration and invasion while promoting BC cell apoptosis in vitro. Overexpression of circ_0008945 remarkably promoted BC tumor growth in vivo. Mechanistically, circ_0008945 acted as a miRNA sponge for miR-338-3p and inhibited its expression in BC cells. Moreover, miR-338-3p targeted and inhibited HOXA3.

CONCLUSION

We found that circ_0008945 acted as a BC oncogene by physically binding miR-338-3p, which further targeted and regulated HOXA3.

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.

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.

The clinicopathological significance of decreased miR-125b-5p in hepatocellular carcinoma: evidence based on RT-qPCR, microRNA-microarray, and microRNA-sequencing

The aim of the study was to comprehensively evaluate the clinical value of miR-125b-5p in hepatocellular carcinoma (HCC) and its potential molecular mechanisms. MiR-125b-5p expression was remarkably lower as examined by real-time reverse transcription quantitative polymerase chain reaction (RT-qPCR) in 95 paired HCC and nonmalignant liver tissues in house (P<0.001), which was in accord with the results from miRNA-sequencing data with 371 cases of HCC. miRNA-chips from Gene Expression Omnibus (GEO) and ArrayExpress were screened. Among the seven included miRNA-chips, the relative expression of miR-125b-5p expression levels showed decreasing trends in HCC tissue samples compared with non-cancerous liver tissue samples. Altogether, A total of 655 cases of HCC tissues and 334 non-HCC liver tissues were included in the final meta-analysis. We observed that the expression of miR-125b-5p indeed decreased markedly in HCC tissues compared with the non-HCC tissues (SMD: -1.414, 95% CI: -1.894 to -0.935, P<0.001). The area under the SROC curve of lower expression of miR-125b-5p was 0.91 (95% CI: 0.89 to 0.94). A Kaplan-Meier survival analysis indicated that the lower expression or the absence of miR-125b-5p may be a risk factor for the poor outcome of HCC patients. Furthermore, the potential target genes of miR-125b-5p from 11 miRNA target prediction databases were intersected with 1,486 differentially expressed genes (DEGs) as calculated by RNA-sequencing data. Finally, a total of 330 GEGs were collected and enriched in the pathways of lysosome, focal adhesion, and pathways in cancer. In conclusion, this study utilizes a variety of research methods to confirm the lower level of miR-125b-5p in HCC tissues. This lower expression level of miR-125b-5p is closely related to increased disease progression in HCC patients.