Search

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.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 byAs a result, in 3201 LUAD samples and 2494 non-cancer lung samples, HOXA4 and HOXA5 were significantly downexpressed (P0.05). The aberrant expression of HOXA5 was significantly correlated with the clinical progression of LUAD (P0.05). HOXA5 showed remarkable prognostic value for LUAD patients (P0.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.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.

Background. Currently, the benefits of nasopharyngeal carcinoma (NPC) therapy are limited, and it is necessary to further explore possible therapeutic targets. Aryl hydrocarbon receptor nuclear translocator 2 (ARNT2) has been extensively studied in other cancer species, but little has been explored in NPC. The aim of this study was to verify the expression level of ARNT2 and its underlying mechanism in NPC. Methods. Datasets containing ARNT2 mRNA expression levels were retrieved and collected from various databases to explore the expression status of ARNT2 in NPC. ARNT2-related coexpressed genes, differential expressed genes, and target genes were obtained for functional enrichment analysis. The potential target gene of ARNT2 and their regulatory relationship were studied through ChIP-seq data. CIBERSORTx was used to assess the immune infiltration of NPC, and the association with ARNT2 expression was calculated through correlation analysis. Results. ARNT2 was upregulated and possessed an excellent discriminatory capability in NPC samples. ARNT2 positively correlated target genes were clustered in pathways in cancer, while negatively correlated target genes were enriched in immune-related pathway. The ChIP-seq information of ARNT2 and histone showed that prostaglandin-endoperoxide synthase 2 (PTGS2) was a potential target gene of ARNT2. CIBERSORTx revealed the immunity status in NPC, and ARNT2 expression was correlated with infiltration of five immune cells. Conclusions. ARNT2 is overexpressed in NPC and may regulate PTGS2 to participate in the cancer process. ARNT2 serves as a key oncogenic target in NPC patients.

BackgroundEsophageal squamous cell carcinoma (ESCC) ranks the sixth in mortality rates in cancers due to a lack of a specific target of diagnosis and treatment in the early stages. Although Forkhead box M1 (FOXM1) has been reported to be differentially expressed in ESCC, its clinical role and function in ESCC remained unclarified.MethodsData from our hospital and public databases (n = 1906) were combined to estimate how FOXM1 overexpression showed its discriminatory ability between ESCC and non-ESCC esophageal tissues. Downstream targets of FOXM1 were predicted by using Cistrome database. Functional enrichment analyses were performed to explore the potential signaling pathways related to FOXM1 in ESCC. Based on the available clinical parameters, we investigated the prognosis potential of FOXM1 and its targets.ResultsThe pooled standard mean difference (SMD) for FOXM1 is 2.62 (95% CI: 2.08–3.16), indicating that FOXM1 is upregulated in ESCC. FOXM1 has an extremely high discrimination potential in ESCC because the area under the curve (AUC) of the summary receiver operating characteristic curve (sROC) is 0.99 (95% CI: 0.97–0.99). A total of 168 downstream targets were identified, and nine hub genes were screened from them. We found that FOXM1 and its targets were significantly enriched in the cell cycle. Additionally, the correlation between FOXM1 and clinical parameters had not been observed, except for age.Conclusions FOXM1 is upregulated in ESCC and has an extremely high discrimination potential in ESCC.

Jing-Xiao Li,1 Jin-Shu Pang,2 Bin-Tong Yin,1 Gang Chen,1 Jun-Hong Chen,1 Jia-Yuan Luo,1 Xia Yang,1 Li-Ting Qin,1 Jiang-Hui Zeng,3 Peng Chen,4 Jia-Bo Chen,4 Deng Tang1 1Department of Pathology, The First Affiliated Hospital of Guangxi Medical University, Nanning, Guangxi Zhuang Autonomous Region, People’s Republic of China; 2Department of Medical Ultrasonics, The First Affiliated Hospital of Guangxi Medical University, Nanning, Guangxi Zhuang Autonomous Region, People’s Republic of China; 3Department of Clinical Laboratory, The Third Affiliated Hospital of Guangxi Medical University/Nanning Second People’s Hospital, Nanning, Guangxi Zhuang Autonomous Region, 530031, People’s Republic of China; 4Department of Pediatric Surgery, The First Affiliated Hospital of Guangxi Medical University, Nanning, Guangxi Zhuang Autonomous Region, 530031, People’s Republic of ChinaCorrespondence: Deng TangDepartment of Pathology, The First Affiliated Hospital of Guangxi Medical University, Nanning, Guangxi Zhuang Autonomous Region, People’s Republic of ChinaTel +86-7715356534Email tdnever@163.comPurpose: The molecular mechanisms and signal pathways of ferroptosis in hepatoblastoma (HB) have not yet been clarified. In previous studies, activating transcription factor 3 (ATF3) was reported to be correlated with several tumors, but the clinical significance of ATF3 has never been determined. Herein, we investigated the clinicopathological value and mechanisms of ATF3 in regulating ferroptosis in HB.Methods: The mRNA microarray and RNA-sequencing data of 402 samples from our hospital and public databases were used to estimate ATF3 expression and assess its clinical role in HB. The standard mean difference (SMD) and summary receiver operating characteristic curves were utilized to judge the discrimination ability of ATF3 between HB and non-HB liver tissues. We examined the expression variation of ATF3 in HB cells after the treatment with erastin. We also predicted the target genes of ATF3 as a transcriptional factor from public Chromatin Immunoprecipitation-sequencing data and selected the ferroptosis-related genes for a signaling pathway analysis.Results: In ten series, the pooled SMD for ATF3 was − 0.91, demonstrating that ATF3 expression was predominantly lower in HB than in non-HB liver tissues. ATF3 down-regulation showed moderate potential to distinguish HB from non-HB liver tissues (area under curves = 0.83, 95% confidence interval = 0.79– 0.86). Altogether, 4855 putative targets of ATF3 as a transcriptional factor were collected, among which, 60 genes were ferroptosis-related.Conclusion: The down-regulated ATF3 expression may play a vital role in the occurrence of HB possible partially by regulating ferroptosis.Keywords: transcription factor, ATF3, hepatoblastoma, ferroptosis

Background Acid phosphatase type 6 (ACP6) is a mitochondrial lipid phosphate phosphatase that played a role in regulating lipid metabolism and there is still blank in the clinico-pathological significance and functional roles of ACP6 in human cancers. No investigations have been conducted on ACP6 in hepatocellular carcinoma (HCC) up to date. Methods Herein, we appraised the clinico-pathological significance of ACP6 in HCC via organizing expression profiles from globally multi-center microarrays and RNA-seq datasets. The molecular basis of ACP6 in HCC was explored through multidimensional analysis. We also carried out in vitro and in vivo experiment on nude mice to investigate the effect of knocking down ACP6 expression on biological functions of HCC cells, and to evaluate the expression variance of ACP6 in xenograft of HCC tissues before and after the treatment of NC. Results ACP6 displayed significant overexpression in HCC samples (standard mean difference (SMD) = 0.69, 95% confidence interval (CI) = 0.56–0.83) and up-regulated ACP6 performed well in screening HCC samples from non-cancer liver samples. ACP6 expression was also remarkably correlated with clinical progression and worse overall survival of HCC patients. There were close links between ACP6 expression and immune cells including B cells, CD8 + T cells and naive CD4 + T cells. Co-expressed genes of ACP6 mainly participated in pathways including cytokine-cytokine receptor interaction, glucocorticoid receptor pathway and NABA proteoglycans. The proliferation and migration rate of HCC cells transfected with ACP6 siRNA was significantly suppressed compared with those transfected with negative control siRNA. ACP6 expression was significantly inhibited by nitidine chloride (NC) in xenograft HCC tissues. Conclusions ACP6 expression may serve as novel clinical biomarker indicating the clinical development of HCC and ACP6 might be potential target of anti-cancer effect by NC in HCC.

Background: The clinical value of pyroptosis-related genes (PRGs) in lung adenocarcinoma (LUAD) remains obscure. Objective: The study attempts to explore PRGs in LUAD, which will enable an understanding of LUAD from the perspective of PRGs. Method: Lung adenocarcinoma patients were diagnosed using pathology, and their clinical information was collected from several public databases. A PRGs prognostic signature (PPS) for LUAD patients was established based on a multivariate Cox regression analysis. The differential expression of PRGs was identified using standardized mean differences in 6,958 samples. The area under the curve (AUC) was used to evaluate the predictive effects of the PPS to determine the survival rate of LUAD patients. Decision curve analysis was utilized to assess the clinical significance of the PPS in LUAD. Results: The PPS consists of five PRGs, namely CASP3, CASP9, GSDMB, NLRP1, and TNF. The prognostic effect of the PPS is evident in all the predicted one-, three-, and five-year survival rates (AUCs ≥ 0.58). The PPS represents an independent risk factor for the prognosis of LUAD patients (hazard ratio > 1; 95% confidence interval excluding 1). The PPS risk score can predict the prognosis of LUAD patients more accurately than PRGs of the PPS and multiple clinical parameters, such as age, tumor stage, and clinical stage. The decision curve analysis revealed that the nomogram based on the PPS and clinical parameters might result in better clinical decisions. Conclusion: The PPS makes it feasible to distinguish LUAD from non-LUAD. Thus, the underlying significance of the PPS in distinguishing LUAD from non-LUAD is promising.

miR-15b-5p has frequently been reported to function as a biomarker in some malignancies; however, the function of miR-15b-5p in hepatocellular carcinoma (HCC) and its molecular mechanism are still not well understood. The present study was designed to confirm the clinical value of miR-15b-5p and further explore its underlying molecular mechanism. A comprehensive investigation of the clinical value of miR-15b-5p in HCC was investigated by data mining The Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO) datasets as well as literature. In addition, intersected target genes of miR-15b-5p were predicted using the miRWalk database and differentially expressed genes of HCC from TCGA. Furthermore, gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses were carried out. Then, a protein-protein interaction network (PPI) was constructed to reveal the interactions between some hub target genes of miR-15b-5p. The miR-15b-5p expression level in HCC was predominantly overexpressed compared with non-HCC tissues samples (SMD=0.618, 95% CI: 0.207, 1.029; P

Background Lung adenocarcinoma (LUAD) is the main pathological type of pulmonary malignant tumors. Synaptojanin 2 (SYNJ2), a member of the synaptojanin family, is considered as an attractive underlying therapeutic target in several types of cancer. The purpose of the current research was to comprehensively investigate the molecular function and underlying mechanism of SYNJ2 in LUAD with the goal of providing a promising therapeutic target for LUAD. Methods We applied tissue microarrays, immunohistochemistry (IHC), high-throughput RNA sequencing (RNA-seq) data, and gene microarrays of public databases to comprehensively examine the protein level, mRNA expression, clinical significance, and prognosis value of SYNJ2 in LUAD. The standard mean difference (SMD) and summary receiver operating characteristic curves (SROC) were computed to evaluate the comprehensive expression value of SYNJ2 in LUAD. Kaplan-Meier survival curves were plotted to determine the overall survival of SYNJ2 with different expression levels in LUAD. Univariate and multivariate cox regression analyses were implemented to investigate the dependent prognostic ability of SYNJ2 expression and the clinically-relevant traits of LUAD. Moreover, SYNJ2 co-expressed genes were screened for Gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) function analysis and protein–protein interactions (PPI) network construction. Results SYNJ2 overexpression in LUAD was confirmed by 2,299 LUAD samples and 1,387 non-tumor tissues obtained from tissue microarrays along with TCGA and public gene microarrays. The pooled SMD was 1.29 (95% CI: 1.09-1.48, I-square=78.6%, P

Objective. miRNA has gained attention as a therapeutic target in various malignancies. The proposal of this study was to investigate the biological functions of key miRNAs and target genes in cancers of the digestive tract which include esophageal carcinoma (ESCA), gastric adenocarcinoma (GAC), colon adenocarcinoma (COAD), and rectal adenocarcinoma (READ). Materials and Methods. After screening differentially expressed miRNAs (DEMIs) and differentially expressed mRNAs (DEMs) in four digestive cancers from The Cancer Genome Atlas (TCGA) database, the diagnostic value of above DEMIs was evaluated by receiver-operating characteristic (ROC) curve analysis. Then, corresponding DEMIs’ target genes were predicted by miRWalk 2.0. Intersection of predicted target genes and DEMs was taken as the target genes of DEMIs, and miRNA-mRNA regulatory networks between DEMIs and target genes were constructed. Meanwhile, the univariate Cox risk regression model was used to screen miRNAs with distinct prognostic value, and Kaplan–Meier analysis was used to determine their significance of prognosis. Furthermore, we performed bioinformatics methods including protein-protein interaction (PPI) networks, gene ontology (GO) annotation, Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis, and gene group RIDA analysis by Gene-Cloud of Biotechnology Information (GCBI) to explore the function and molecular mechanisms of DEMIs and predicted target genes in tumor development. Results. Eventually, 3 DEMIs (miR-7-3, miR-328, and miR-323a) with significant prognostic value were obtained. In addition, 3 DEMIs (miR-490-3p, miR-133a-3p, and miR-552-3p) and 281 target genes were identified, and the 3 DEMIs showed high diagnostic value in READ and moderate diagnostic value in ESCA, GAC, and COAD. Also, the miRNA-mRNA regulatory network with 3 DEMIs and 281 overlapping genes was successfully established. Functional enrichment analysis showed that 281 overlapping genes were mainly related to regulation of cell proliferation, cell migration, and PI3K-Akt signaling pathway. Conclusion. The diagnostic value and prognostic value of significant DEMIs in cancers of the digestive tract were identified, which may provide a novel direction for treatment and prognosis improvement of cancers of the digestive tract.

MicroRNAs (miRNAs) play key roles in the regulation of gene expression during multiple physiological processes, including early development, differentiation, and ageing. However, their involvement in age-related thymic involution is not clear. In this study, we profiled the global transcriptome and miRNAome of thymic epithelial cells in 1- and 3-month-old male and female mice, and predicted the possible transcription factors and target genes of the four most significantly differentially expressed miRNAs (DEMs) (miR-183-5p, miR-199b-5p, miR-205-5p, and miR-200b-3p) by performing bioinformatics analyses. We also evaluated the relationships between the significantly DEMs and mRNAs. We performed quantitative polymerase chain reaction to confirm the changes in the expression of the miRNAs and their predicted target genes. We found that miR-183-5p, miR-199b-5p, miR-205-5p, and miR-200b-3p can be used as a biomarker group for mouse thymus development and involution. In addition, the predicted target genes (Ptpn4, Slc2a9, Pkib, Pecam1, and Prkdc), which were identified by mRNA sequencing analysis, were mainly involved in growth, development, and accelerated senescence. In conclusion, miRNAs and their predicted target genes likely play important roles in thymus development and involution. To the best of our knowledge, this is the first study to systematically analyze the relevance of miRNAs and their targets by mRNA sequencing in mouse thymic epithelial cells. © 2018 IUBMB Life, 70(7):678-690, 2018.

To practically verify the clinical value of miR-26a-5p and thoroughly explore its target genes as well as its potential functions in hepatocellular carcinoma (HCC).HCC and adjacent non-cancerous hepatic tissues of 95 HCC patients were collected for analysis using reverse transcription quantitative real-time PCR (qRT-PCR). For the bioinformatics analysis, we identified potential target genes for miR-26a-5p from differentially expressed genes (DEGs) in Gene Expression Omnibus (GEO) data sets and miRWalk predicted database. Gene ontology (GO), Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway and protein-protein interaction (PPI) analyses were applied to analyze the prospective mechanisms of the predicted target genes.MiR-26a-5p showed a significantly lower expression level in HCC tissues (1.56±1.07) than adjacent benign liver tissues (2.28±1.06, P0.001). The area under the curve (AUC) of the receiver operating characteristic (ROC) was 0.665 (95% CI: 0.588-0.743, P0.001). Significant correlations between miR-26a-5p expression and clinicopathological features such as gender (r=0.275, P0.01), clinical TNM stage (r=-0.306, P0.01), and metastasis (r=-0.321, P0.01) were observed. To examine potential target genes, we obtained 175 genes for further function analysis, by attaining the intersection of 2062 up-regulated DEGs and 1390 online-predicted target genes. The GO and KEGG pathway annotation indicated focal adhesion, regulation of actin cytoskeleton and the PI3K-Akt signaling pathway as significant prospective mechanisms. The PPI network indicated that NRAS was the most essential hub gene in the whole network.Down-regulated miR-26a-5p was closely correlated with the status of metastasis and the progression of HCC. MiR-26a-5p might play protective roles by targeting diverse genes and pathways.

Kidney renal papillary cell carcinoma (KIRP) accounts for 10%-15% of renal cell carcinoma (RCC), patients with KIRP tend to have a poor prognosis, and there was a lack of effective prognostic indicators for this type of cancer. Currently, owing to the availability of The Cancer Genome Atlas (TCGA), long non-coding RNAs (LncRNAs) have been discovered to indicate a prognostic value in some tumors. In that regard, we analyzed lncRNA-sequencing data of KIRP in TCGA, and among 780 differentially-expressed lncRNAs, we selected 37 lncRNAs which were able to assist the prognosis. In addition, by using the multivariate cox regression analysis, the prognosis index (PI) that consisted of 7 lncRNAs (including AFAP1-AS1, GAS6-AS1, RP11-1C8.7, RP11-21L19.1, RP11-503C24.1, RP11-536I6.2, and RP11-63A11.1) could predict the progression and outcomes of KIRP with accuracy. More importantly, the PI was considered an independent indicator for prognostication of KIRP. Moreover, having categorized patients with KIRP into cohorts of high risk and low risk, according to the PI, we found that the key genes and pathways varied in these two groups. Overall, these LncRNAs, especially the PI, may be conceived as biomarkers and helpful for determining the different pathological stages for KIRP patients. However, their biological functions need to be further confirmed.

Lung squamous cell carcinoma (LUSC) accounts for a significant proportion of lung cancer and there have been few therapeutic alternatives for recurrent LUSC due to the lack of specific driver molecules. To investigate the prospective role of lncRNAs in the tumorigenesis and progression of LUSC, the aberrantly expressed lncRNAs were calculated based on The Cancer Genome Atlas RNA-seq data. Of 7589 lncRNAs with 504 LUSC cases, 884 lncRNAs were identified as being aberrantly expressed (|log2 fold change| >2 and adjusted P

Long non-coding RNAs (lncRNAs) expression profile signature for survival assessment in lung squamous cell carcinoma (LUSC) are largely inconsistent due to distinct detecting approaches and small sample size. Systematic and integrative investigation of RNA-Seq based data from The Cancer Genome Atlas (TCGA) herein was performed to determine candidate lncRNAs for prognosis evaluation of LUSC. A total of 60483 genes, including 7589 lncRNAs were assessed in a cohort including 478 LUSC cases with follow-up data. Firstly, 4225 differentially expressed lncRNAs were obtained via R packages. Next, univariate and multivariate Cox proportional hazards regression revealed that 41 lncRNAs were closely related to the survival of LUSC. Finally, lncRNA based prognosis index (PI) could predict overall survival of LUSC with high accuracy (AUC = 0.652, CI: 0.598, 0.705), PI = expCYP4F26P*βCYP4F26P+expRP11-108M12.3*βRP11-108M12.3+expRP11-38M8.1*βRP11-38M8.1+expRP11-54H7.4*βRP11-54H7.4+expZNF503-AS1*βZNF503-AS1. Furthermore, it was confirmed that the five-lncRNA signature could act as an independent prognostic indicator for LUSC (HR = 2.068, p < 0.001 with univariate analysis, HR = 1.928, p = 0.038 with multivariate). Besides, we constructed a weighted gene co-expression network analysis (WGCNA) of key lncRNA RP11-54H7.4 according to the p-value of related genes' weight. This study provides a RNA-Seq based prognostic signature with five lncRNAs for further clinical application to LUSC patients.

// Jiang-hui Zeng 1, * , Dan-dan Xiong 1, * , Yu-yan Pang 1 , Yu Zhang 1 , Rui-xue Tang 1 , Dian-zhong Luo 1 , Gang Chen 1 1 Department of Pathology, First Affiliated Hospital of Guangxi Medical University, Nanning, Guangxi Zhuang Autonomous Region 530021, China * These authors contributed equally to this work Correspondence to: Gang Chen, email: chen_gang_triones@163.com Dian-zhong Luo, email: 13878802796@163.com Keywords: esophageal carcinoma, DEMs, DEGs, TCGA, diagnosis Received: November 08, 2016 Accepted: February 28, 2017 Published: March 09, 2017 ABSTRACT Esophageal carcinoma (ESCA) is one of the most common malignancies worldwide, and its pathogenesis is complex. In this study, we identified differentially expressed miRNAs (DEMs) and genes (DEGs) of ESCA from The Cancer Genome Atlas (TCGA) database. The diagnostic values of DEMs were determined by receiver operating characteristic (ROC) analyses and validated based on data from Gene Expression Omnibus (GEO). The top five DEMs with the best diagnostic values were selected, and their potential targets were predicted by various in silico methods. These target genes were then identified among the DEGs from TCGA. Furthermore, the overlapping genes were subjected to protein-protein interaction (PPI) analysis, Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses. The miRNA-transcription factor (TF) regulatory relations were determined using CircuitsDB and TransmiR. Finally, the regulatory networks of miRNA-TF and miRNA-gene were constructed and analyzed. A total of 136 DEMs and 3541 DEGs were identified in ESCA. The top five DEMs with the highest area under the receiver operating characteristic curve (AUC) values were miRNA-93 (0.953), miRNA-21 (0.928), miRNA-4746 (0.915), miRNA-196a-1 (0.906) and miRNA-196a-2 (0.906). The combined AUC of these five DEMs was 0.985. The KEGG analysis with 349 overlapping genes showed that the calcium signaling pathway and the neuroactive ligand-receptor interaction were the most relevant pathways. The regulatory networks of miRNA-TF and miRNA-gene, including 38 miRNA-TF and 560 miRNA-gene pairs, were successfully established. Our findings may provide new insights into the molecular mechanisms of ESCA pathogenesis. Future research will aim to explore the role of novel miRNAs in the pathogenesis and improve the early diagnosis of ESCA.

The cancer susceptibility candidate 9 (CASC9) gene has been reported to exert an oncogenic effect in several types of cancer. However, its role in lung squamous cell carcinoma (LUSC) is unknown. Therefore, the present study examined the expression of CASC9 in LUSC and non-cancer tissues by reverse transcription-quantitative polymerase chain reaction assays and by data mining of high-throughput public databases, including The Cancer Genome Atlas, the Gene Expression Omnibus, ArrayExpress and the Cancer Cell Line Encyclopedia. In vitro experiments were conducted to investigate the effects of CASC9 on the viability and the proliferation of LUSC cells. Furthermore, consulting the alteration status of CASC9 in LUSC from cBioPortal, functional enrichment analysis of co-expressed genes, prediction of potential transcription factors, and inspection of adjacent protein-coding genes were conducted to explore the potential molecular mechanism of CASC9 in LUSC. The results revealed that CASC9 was overexpressed in LUSC tissue, and significantly associated with the malignant progression of LUSC. In vitro experiments demonstrated that CASC9 knockdown by RNA interference attenuated the viability and proliferation of LUSC cells. Multiple copies of CASC9 gene were detected in 4 of 179 available sequenced LUSC cases. A functional enrichment analysis of 200 co-expressed genes indicated that these genes were significantly associated with terms, including 'cell-cell junction organization', 'desmosome organization', 'epidermis development', 'Hippo signaling pathway', 'pathogenic Escherichia coli infection' and 'PID HIF1 TF pathway'. Three genes, Fos-related antigen 2 (FOSL2), SWI/SNF complex subunit SMARCC2, and transcription factor COE1 (EBF1), were predicted by lncRNAMap to be associated with CASC9. Among these, the expression of FOSL2 and EBF1 was positively and negatively correlated with the expression of CASC9, respectively. Two adjacent protein-coding genes, cysteine-rich secretory protein LCCL domain-containing 1 and hepatocyte nuclear factor 4-γ, were also positively correlated with CASC9 expression. In conclusion, the present data suggest that CASC9 serves as an oncogene in LUSC and may be a promising target for alternative therapeutic options for patients with this condition.

Background The scientific understanding of long non-coding RNAs (lncRNAs) has improved in recent decades. Nevertheless, there has been little research into the role that lncRNAs play in clear cell renal cell carcinoma (ccRCC). More lncRNAs are assumed to influence the progression of ccRCC via their own molecular mechanisms. Methods This study investigated the prognostic significance of differentially expressed lncRNAs by mining high-throughput lncRNA-sequencing data from The Cancer Genome Atlas (TCGA) containing 13,198 lncRNAs from 539 patients. Differentially expressed lncRNAs were assessed using the R packages edgeR and DESeq. The prognostic significance of lncRNAs was measured using univariate Cox proportional hazards regression. ccRCC patients were then categorized into high- and low-score cohorts based on the cumulative distribution curve inflection point the of risk score, which was generated by the multivariate Cox regression model. Samples from the TCGA dataset were divided into training and validation subsets to verify the prognostic risk model. Bioinformatics methods, gene set enrichment analysis, and protein–protein interaction networks, Gene Ontology, and Kyoto Encyclopedia of Genes and Genomes analyses were subsequently used. Results It was found that the risk score based on 6 novel lncRNAs (CTA-384D8.35, CTD-2263F21.1, LINC01510, RP11-352G9.1, RP11-395B7.2, RP11-426C22.4) exhibited superior prognostic value for ccRCC. Moreover, we categorized the cases into two groups (high-risk and low-risk), and also examined related pathways and genetic differences between them. Kaplan–Meier curves indicated that the median survival time of patients in the high-risk group was 73.5 months, much shorter than that of the low-risk group (112.6 months; P

Although the molecular studies of single gastrointestinal tumors have been widely reported by media, it is not clear about the function of small nucleolar RNA (snoRNA) in the progression, development and prognostic significance in colon adenocarcinoma, and its certain molecular mechanisms and functions remain to be studied. This study aims to dig out the gene expression data profile of colon adenocarcinoma and construct the prognostic molecular pathology prediction-evaluation, ultimately revealing the clinical prognostic value of snoRNA in colon adenocarcinoma. 932 differentially expressed snoRNAs of the colon adenocarcinoma were obtained by edgeR R package. Only 4 prognostically-significant snoRNAs (SNORD14E, SNORD67, SNORD12C, and SNORD17) (P 0.05) were discovered after univariate COX regression mode analysis. Moreover, through multivariate COX regression mode analysis, 2 prognostically-significant snoRNAs (SNORD14E and SNORD67) (P 0.05) were obtained. Using the above 473 COAD samples, a prognostic model of risk score was constructed. The inflection point of the prognostic risk score acted as a boundary to divide the patients into high-risk and low-risk groups. The K-M survival curve of the prognostic model of risk score revealed that high risk group has a lower survival rate (P 0.05). The research has successfully provided valuable prognostic factors and prognostic models for patients with malignant colon tumor.

In order to determine the diagnostic efficacy of microRNA (miR)-122-5p and to identify the potential molecular signaling pathways underlying the function of miR-122-5p in hepatocellular carcinoma (HCC), the expression profiles of data collected from The Cancer Genome Atlas (TCGA), Gene Expression Omnibus (GEO) and literature databases were analyzed, along with any associations between clinicopathological characteristics and the diagnostic value of miR-122-5p in HCC. The intersection of 12 online prediction databases and differentially expressed genes from TCGA and GEO were utilized in order to select the prospective target genes of miR-122-5p in HCC. Gene Ontology (GO), Kyoto Encyclopedia of Genes and Genomes (KEGG) and protein-protein interaction network (PPI) analyses were subsequently performed based on the selected target genes. The average expression level of miR-122-5p was decreased in HCC patients compared with controls from TCGA database (P

Malignant tumors of the digestive tract include esophageal, gastric, and colorectal carcinomas, which all have high global mortality rates. A clinical role for small nuclear RNA (snRNA), a type of small non-coding RNA, has not yet been documented for digestive tract pan-adenocarcinomas. Therefore, the aim of the study was to identify differentially expressed snRNAs and to explore their prognostic implications in pan-adenocarcinomas from the esophagus, stomach, colon, and rectum. The pan-carcinoma RNA-sequencing data of four types of digestive tract cancers with 1, 102 cases obtained from The Cancer Genome Atlas (TCGA) project were analyzed and the differentially expressed snRNAs were evaluated using the edgeR package. The prognostic value of each of the selected snRNAs was determined by univariate and multivariate Cox regression analyses. All the digestive tract pan-adenocarcinomas showed differential expression of three snRNAs: the up-regulated RNU1-106 P and RNU6-850 P and the down-regulated RNU6-529 P. Interestingly, RNU6-101 P appeared to be a risk factor for esophageal adenocarcinoma (ESAD) and RNVU1-4 was potentially a protective factor for stomach adenocarcinoma (STAD) survival. This consistent finding of differential expression of all three snRNAs in all four types of digestive system cancers suggests potential roles for these snRNAs in the tumorigenesis of digestive system cancers. RNU6-101 P could play a pivotal role in the progression of ESAD and RNVU1-4 could perform a protective role in STAD. However, since the current findings were based on RNA-sequencing data mining, more studies are needed for verification.

Background/Aims: Left- and right-sided colon cancers are considered to be two different diseases and have altered outcomes. However, specific molecules to predict the prognosis of left- and right-sided colon cancers are currently lacking. Methods: Expression profiling of colon cancer were downloaded from The Cancer Genome Atlas (TCGA). Differentially expressed genes (DEGs) of left- and right-sided colon cancers were compared by DESeq analysis. The prognostic values of DEGs were assessed by univariate and multivariate Cox regression. Prognostic index models of two side colon cancers were conducted with prognostic values genes, respectively. Interaction of DEGs was then analyzed by the protein-protein interaction (PPI). Different biology function of two sides of colon cancer was assessed by Gene Set Enrichment Analysis (GSEA). Results: A total of 167 DEGs were identified between left- and right-sided colon cancers based on TCGA data. Using univariate COX regression analysis, five genes (PHACTR3, CKMT2, CYP2W1, ERFE, HOXC4) were related to overall survival in left-sided, and eight distinguishable genes (EREG, ERFE, HOXC6, SLC22A31, TFF1, GFI1, ZG16, RASL10B) in right-sided. Further, left-sided prognostic model was established with PHACTR3 and CKMT2 (HR=2.040; 95%CI=1.004-4.145; P=0.049). Distinguishable prognostic signature for right-sided colon cancer was established based on EREG, ERFE, GFI1, and RASL10B (HR=3.530; 95%CI: 1.934-6.444; P< 0.001) in multivariate analysis. PPI analysis of 167 DEGs showed that CCL5, GNG4, GNLY, GZMH, DRD2, and FASLG genes were at the core of interaction network. In GSEA function analysis, four pathways, including antigen processing and presentation, natural killer cell mediated cytotoxicity, intestinal immune network for Iga production, and type I diabetes mellitus, were significantly enriched in the DEGs of the right-sided colon cancer. Conclusions: This study constructs a panel of potential prognostic model of left- and right-sided colon cancers, respectively. We also provide molecular biological alterations between left- and right-sided colon cancers.