T cell immune regulator 1 is a prognostic marker associated with immune infiltration in glioblastoma multiforme

Glioma is the most common primary brain tumor and glioblastoma multiforme (GBM) is the most malignant brain glioma with the worst prognosis. T cell immune regulator 1 (TCIRG1) constitutes the V₀a₃ subunit of vacuolar ATPase (V-ATPase), and the function of V-ATPase in malignant tumors, such as breast cancer, melanoma and hepatocellular carcinoma, has been reported. However, the effect of the TCIRG1 subunit on GBM remains to be fully elucidated. mRNA levels of TCIRG1 in different cancer types and the corresponding normal tissues were extracted from the Oncomine and Tumor Immune Estimation Resource (TIMER) databases. The Gene Expression Omnibus (access number: {"type":"entrez-geo","attrs":{"text":"GSE16011","term_id":"16011"}}GSE16011), the Chinese Glioma Genome Atlas and The Cancer Genome Atlas were used to investigate the mRNA level of TCIRG1 in glioma. Protein level validation in glioma was performed using western blotting. The Database for Annotation, Visualization and Integrated Discovery was used to analyze Gene Ontology (GO) categories for genes correlated with TCIRG1 in GBM. Protein-protein interaction (PPI) networks and module analyses were performed using Cytoscape software and the MCODE plugin. The correlation between tumor immune cell infiltration and TCIRG1 expression was explored using the TIMER database. Additionally, the correlation between TCIRG1 and the gene signature of immune infiltration was explored through TIMER and Gene Expression Profiling Interactive Analysis. External validation of TCIRG1 expression according to immune signatures in GBM was performed using the {"type":"entrez-geo","attrs":{"text":"GSE16011","term_id":"16011"}}GSE16011 dataset with the GlioVis online tool. It was found that TCIRG1 expression was increased in GBM and numerous malignant tumors and may serve as a biomarker of the mesenchymal subtype of GBM. GO category analysis of positively correlated genes revealed that TCIRG1 was correlated with the immune response in GBM. PPI network and module analyses also supported the potential function of TCIRG1 in the local immune response. The expression of TCIRG1 was associated with various immune markers. It was therefore speculated that TCIRG1 is associated with glioma malignancy and may be a marker of unfavorable prognosis in patients with GBM, and it could be regarded as a prognostic biomarker and an indicator of immune infiltration in GBM.     

Circular RNA hsa_circ_0005114-miR-142-3p/miR-590-5p-adenomatous polyposis coli protein axis as a potential target for treatment of glioma

Glioma is the most common type of brain tumor and is associated with a high mortality rate. Despite recent advances in treatment options, the overall prognosis in patients with glioma remains poor. Studies have suggested that circular (circ)RNAs serve important roles in the development and progression of glioma and may have potential as therapeutic targets. However, the expression profiles of circRNAs and their functions in glioma have rarely been studied. The present study aimed to screen differentially expressed circRNAs (DECs) between glioma and normal brain tissues using sequencing data collected from the Gene Expression Omnibus database ({"type":"entrez-geo","attrs":{"text":"GSE86202","term_id":"86202"}}GSE86202 and {"type":"entrez-geo","attrs":{"text":"GSE92322","term_id":"92322"}}GSE92322 datasets) and explain their mechanisms based on the competing endogenous (ce)RNA regulatory hypothesis. In total, 424 commonly downregulated DECs (with the Gene_symbol annotated in the circBase database) in these two datasets were identified. Using the CircInteractome and Starbase databases, 18 micro (mi)RNAs (miRs) were predicted to interact with DECs, while 22 glioma-related genes obtained from the Comparative Toxicogenomics Database were predicted to be regulated by 15 miRNAs via the miRwalk 2.0 database. A ceRNA network was established based on 115 DECs, 15 miRNAs and 22 mRNAs. LinkedOmics online analysis using The Cancer Genome Atlas (TCGA) data showed that hsa-miR-142-3p/hsa-miR-590-5p and their target gene adenomatous polyposis coli protein (APC) were all significantly associated with overall survival rate and their prognosis trend was opposite, revealing that high expression levels of hsa-miR-142-3p/hsa-miR-590-5 were associated with a poor overall survival rate, while high APC expression with a good overall survival rate. UALCAN analysis using TCGA data of glioblastoma multiforme and the {"type":"entrez-geo","attrs":{"text":"GSE25632","term_id":"25632"}}GSE25632 and {"type":"entrez-geo","attrs":{"text":"GSE103229","term_id":"103229"}}GSE103229 microarray datasets showed that hsa-miR-142-3p/hsa-miR-590-5p was upregulated and APC was downregulated. Thus, hsa-miR-142-3p/hsa-miR-590-5p-APC-related circ/ceRNA axes may be important in glioma, and hsa_circ_0005114 interacted with both of these miRNAs. Functional analysis showed that hsa_circ_0005114 was involved in insulin secretion, while APC was associated with the Wnt signaling pathway. In conclusion, hsa_circ_0005114-miR-142-3p/miR-590-5p-APC ceRNA axes may be potential targets for the treatment of glioma.     

Bioinformatics and immunohistochemistry analyses of expression levels and clinical significance of CXCL2 and TANs in an oral squamous cell carcinoma tumor microenvironment of Prophyromonas gingivalis infection

The present study aimed to detect the immunoexpression and clinical significance of Porphyromonas gingivalis (P. gingivalis) in the tumor microenvironment (TME) of oral squamous cell carcinoma (OSCC). The immunoexpression of P. gingivalis in OSCC tissues was detected via immunohistochemistry (IHC) after P. gingivalis was infected into the TME of OSCC. To identify the differentially expressed genes in the carcinogenesis and progression of OSCC with P. gingivalis infection, microarray datasets ({"type":"entrez-geo","attrs":{"text":"GSE87539","term_id":"87539"}}GSE87539 and {"type":"entrez-geo","attrs":{"text":"GSE138206","term_id":"138206"}}GSE138206) were downloaded from the Gene Expression Omnibus database. The immunoexpression levels of C-X-C motif chemokine ligand 2 (CXCL2) and tumor-associated neutrophils (TANs) were also evaluated via IHC, and the immunoexpression levels of all three clinical variables were analyzed using χ² or Fisher''s exact tests. The survival rates were calculated using the Kaplan-Meier method and the survival curves were compared using log-rank tests. Predominantly strong immunoexpression of P. gingivalis was identified in OSCC samples. CXCL2 was considered to be a differential gene in the two datasets. Immunoexpression of P. gingivalis was positively associated with CXCL2 and TANs expression. Furthermore, P. gingivalis was associated with survival status (P<0.001) and differentiation (P<0.001). CXCL2 was associated with age (P=0.038) and survival status (P=0.003), while TANs were associated with T stage (P=0.015) and clinical stage (P=0.002). These clinical variables were considered to be independent risk factors for the poor prognosis of patients with OSCC. Collectively, the results suggested that the immunoexpression of P. gingivalis may be positively associated with CXCL2 and TANs. In addition, the strong immunoexpression levels of P. gingivalis, CXCL2 and TANs may be associated with a poor prognosis in patients with OSCC.     

A novel five-gene score to predict complete pathological response to neoadjuvant chemotherapy in ER-positive/HER2-negative breast cancer

Neoadjuvant Chemotherapy (NAC) is not frequently used in ER-positive/HER2-negative breast cancer (BC) because around 10% patients achieve pathological complete response (pCR). Since NAC can result in cancer downstaging both in the breast and axilla and prevent a morbid surgery, thus a score to predict pCR in this population will be crucial to identify patients who can benefit from this approach. A total of 4038 patients from cohorts; {"type":"entrez-geo","attrs":{"text":"GSE25066","term_id":"25066"}}GSE25066, {"type":"entrez-geo","attrs":{"text":"GSE20194","term_id":"20194"}}GSE20194, Hess, {"type":"entrez-geo","attrs":{"text":"GSE20181","term_id":"20181"}}GSE20181, TCGA-BRCA and METBRIC were analyzed. The score was generated by the 5 most highly expressed genes in the Hallmark E2F targets gene set amongst patients in the {"type":"entrez-geo","attrs":{"text":"GSE25066","term_id":"25066"}}GSE25066 cohort with ER-positive/HER2-negative BC who achieved pCR. The area under the curve was significantly higher in the score than that for the E2F targets score. High score ER-positive/HER2-negative BCs were significantly associated with higher Nottingham pathological grade, AJCC cancer stage, MKI67 expression levels, intratumor heterogeneity, homologous recombination defects, mutation burden, neoantigen load, and infiltration of anti-cancer immune cells (CD4⁺, T helper type1, plasmacytoid dendritic cells, M1 macrophages). They also expressed lower abundance of stromal cells including fibroblasts, lymphatic endothelial cells, pericytes and adipocytes consistently in {"type":"entrez-geo","attrs":{"text":"GSE25066","term_id":"25066"}}GSE25066, TCGA and METABRIC cohorts. All cell proliferation-related gene sets, G2M checkpoint, E2F targets, MYC targets v1 and v2, Mitotic Spindle, were strongly enriched in high score BCs consistently in 3 cohorts. The gene score was significantly associated with high pCR rate consistently in the {"type":"entrez-geo","attrs":{"text":"GSE25066","term_id":"25066"}}GSE25066 (38%, P < 0.001), {"type":"entrez-geo","attrs":{"text":"GSE20194","term_id":"20194"}}GSE20194 (16%, P = 0.006), and Hess cohort (23%, P = 0.037). In conclusion, the 5-gene score reflects cancer cell proliferation and immune cell infiltration, and predicts pCR after NAC in ER-positive/HER2-negative breast cancer.     

Identification of gene signatures and potential therapeutic targets for acquired chemotherapy resistance in gastric cancer patients

BackgroundGastric cancer (GC) is the most common type of gastrointestinal cancer, and has been studied extensively. However, resistance to chemotherapeutic agents has become a major problem, leading to treatment failure. This study aimed to investigate the molecular mechanisms mediating acquired resistance to cisplatin and fluorouracil (CF) combination-based chemotherapy in GC patients.MethodsThe microarray datasets ({"type":"entrez-geo","attrs":{"text":"GSE14209","term_id":"14209"}}GSE14209, {"type":"entrez-geo","attrs":{"text":"GSE30070","term_id":"30070"}}GSE30070) were downloaded from the Gene Expression Omnibus (GEO) database to identify differentially expressed genes (DEGs) and differentially expressed miRNAs (DEMs) using the limma package in R/Bioconductor. Possible targets of the DEMs were predicted using miRWalk, and the putative miRNA-mRNA regulatory network was constructed using Cytoscape software. Gene Ontology (GO), Kyoto Encyclopedia of Genes and Genomes (KEGG), and protein-protein interaction (PPI) analyses were then conducted and visualized using the Search Tool for Retrieval of Interacting Genes (STRING) and Cytoscape. The prognostic value of hub genes was revealed by Kaplan-Meier Plotter. The causal relationships and interactions between proteins were displayed using DisNor. Finally, similarity analysis was conducted using the Connectivity Map (CMap) profiles to predict a group of small molecules in GC treatment.ResultsA total of 394 DEGs and 31 DEMs were identified after analysis of pre- and post-treatment samples of clinical responders to CF therapy. TM9SF4, hsa-miR-185-5p, and hsa-miR-145-5p were found to be critical in the miRNA-mRNA regulatory network. The DEGs were found to be mainly enriched in the processes of ribonucleoprotein complex assembly, catalytic activity acting on RNA, mitochondrial matrix, and thermogenesis. The DEMs were predominantly found to be involved in single-stranded RNA binding and endoplasmic reticulum lumen. HDAC5, DDX17, ILF3, and SDHC were identified as hub genes in the PPI network. Of these, HDAC5, DDX17, and ILF3 were found to be closely related to the overall survival of GC patients. DisNor identified the first neighbors of the key genes. Furthermore, CMap profiles predicted a group of small molecules, including several histone deacetylase inhibitors (HDACIs), menadione, and mibefradil, which could serve as promising therapeutic agents to reverse acquired resistance to CF therapy.ConclusionsOur findings reveal new targets and alternative therapies to overcome the acquired resistance of GC patients to CF treatment.     

Molecular subtypes based on CNVs related gene signatures identify candidate prognostic biomarkers in lung adenocarcinoma

The classical factors for predicting prognosis currently cannot meet the developing requirements of individualized and accurate prognostic evaluation in lung adenocarcinoma (LUAD). With the rapid development of high-throughput DNA sequencing technologies, genomic changes have been discovered. These sequencing data provide unprecedented opportunities for identifying cancer molecular subtypes. In this article, we classified LUAD into two distinct molecular subtypes (Cluster 1 and Cluster 2) based on Copy Number Variations (CNVs) and mRNA expression data from the Cancer Genome Atlas (TCGA) based on non-negative matrix factorization. Patients in Cluster 1 had worse outcomes than that in Cluster 2. Molecular features in subtypes were assessed to explain this phenomenon by analyzing differential expression genes expression pattern, which involved in cellular processes and environmental information processing. Analysis of immune cell populations suggested different distributions of CD4+ T cells, CD8+ T cells, and dendritic cells in the two subtypes. Subsequently, two novel genes, TROAP and RASGRF1, were discovered to be prognostic biomarkers in TCGA, which were confirmed in {"type":"entrez-geo","attrs":{"text":"GSE31210","term_id":"31210"}}GSE31210 and Tianjin Medical University Cancer Institute and Hospital LUAD cohorts. We further proved their crucial roles in cancers by vitro experiments. TROAP mediates tumor cell proliferation, cycle, invasion, and migration, not apoptosis. RASGRF1 has a significant effect on tumor microenvironment. In conclusion, our study provides a novel insight into molecular classification based on CNVs related genes in LUAD, which may contribute to identify new molecular subtypes and target genes.     

Four potential microRNAs affect the progression of pancreatic ductal adenocarcinoma by targeting MET via the PI3K/AKT signaling pathway

Pancreatic ductal adenocarcinoma (PDAC) is the most common tumor subtype of pancreatic cancer, which exhibits poor patient prognosis due to the lack of effective biomarkers in the diagnosis and treatment. The present study aimed to identify the potential biomarkers of PDAC carcinogenesis and progression using three microarray datasets, {"type":"entrez-geo","attrs":{"text":"GSE15471","term_id":"15471"}}GSE15471, {"type":"entrez-geo","attrs":{"text":"GSE16515","term_id":"16515"}}GSE16515 and {"type":"entrez-geo","attrs":{"text":"GSE28735","term_id":"28735"}}GSE28735, which were downloaded from the Gene Expression Omnibus database. The datasets were analyzed to screen out differentially expressed genes (DEGs) in PDAC tissues and adjacent normal tissues. A total of 143 DEGs were identified, including 132 upregulated genes and 11 downregulated genes. Gene Ontology and Kyoto Encyclopedia of Genes and Genomes functional and signaling pathway enrichment analyses were performed on the DEGs, and the Search Tool for the Retrieval of Interacting Genes/Proteins database was used to construct a protein-protein interaction network. The main functions of DEGs include extracellular matrix degradation, and regulation of matrix metalloproteinase activity and the PI3K-Akt signaling pathway. The five hub genes were subsequently screened using Cytoscape software, and survival analysis demonstrated that abnormal expression levels of the hub genes was associated with poor disease-free survival and overall survival. Biological experiments were performed to confirm whether mesenchymal-to-epithelial transition (MET) factors promote the proliferation, migration and invasion of PDAC cells via the PI3K/AKT signaling pathway. In addition, six MET-targeted microRNAs (miRNAs) were identified, four of which had conserved binding sites with MET. Based on the signaling pathway enrichment analysis of these miRNAs, it is suggested that they can affect the progression of PDAC by targeting MET via the PI3K/AKT signaling pathway. In conclusion, the hub genes and miRNAs that were identified in the present study contribute to the molecular mechanisms of PDAC carcinogenesis and progression. They also provide candidate biomarkers for early diagnosis and treatment of patients with PDAC.     

IL-17 signaling pathway plays a key role in laryngeal squamous cell carcinoma with ethnic specificity

Laryngeal squamous cell carcinoma (LSCC) is a common aggressive head and neck squamous cell carcinoma (HNSCC) and racial disparities have been reported to exist in it. However, its molecular mechanism and associated ethnic specificity are still unclear. Here, we leveraged mRNA expression data from 2 gene expression omnibus datasets ({"type":"entrez-geo","attrs":{"text":"GSE142083","term_id":"142083"}}GSE142083 & {"type":"entrez-geo","attrs":{"text":"GSE117005","term_id":"117005"}}GSE117005) of Chinese samples and the cancer genome atlas (TCGA) datasets of Caucasian samples to demonstrate the expression signature of LSCC. The {"type":"entrez-geo","attrs":{"text":"GSE142083","term_id":"142083"}}GSE142083 dataset was used as the discovery set since it had 53 pairs of LSCC tissues and matched adjacent normal tissues, and the {"type":"entrez-geo","attrs":{"text":"GSE117005","term_id":"117005"}}GSE117005 dataset was treated as the validation set with 5 pairs of tissues. Differential gene expression analysis and enrichment pathway analysis were performed. Besides, we employed weighted gene co-expression network analysis to identify hub genes in validated pathways. The TCGA datasets were used to evaluate ethnic specificity. Immunohistochemistry (IHC) was employed to further validate the hub gene. Overall, the IL-17 signaling pathway was significantly enriched for upregulated genes in two Chinese datasets while not in TCGA datasets; and IL17RC, MAPK3, S100A8, MMP3, CXCL8, and TNFA1P3 were hub genes regulating such pathway. Therein, IL17RC might be the most important one and the IHC results displayed that the IL17RC gene upregulated in the LSCC tissue. IL-17 signaling pathway has an ethnicity-specific effect in LSCC where it is upregulated in the Chinese while not in the Caucasians and IL17RC might play a key role. Targeting genes located in the IL-17 signaling pathway may be beneficial for Chinese LSCC patients.     

Identification of a 3-gene signature for predicting the prognosis of stage II colon cancer based on microsatellite status

BackgroundThough colon cancer (CC) is one of the most malignant tumors across the world, CC patients with microsatellite instability-high (MSI-H) in stage II seem to have a better prognosis. However, the molecular mechanisms underlying the phenomena haven’t been elucidated yet.MethodsThis study enrolled 322 CCs with known microsatellite status from {"type":"entrez-geo","attrs":{"text":"GSE143985","term_id":"143985"}}GSE143985, {"type":"entrez-geo","attrs":{"text":"GSE39582","term_id":"39582"}}GSE39582 and {"type":"entrez-geo","attrs":{"text":"GSE92921","term_id":"92921"}}GSE92921 in the Gene Expression Omnibus (GEO) database. Robust rank aggregation (RRA) analysis, univariate Cox regression analysis and multivariate Cox stepwise regression analysis were performed to identify genes and construct risk score signature. Kaplan-Meier and receiver operating characteristic (ROC) curves analyses were used to evaluate the prognostic value of the signature. The potential mechanisms underlying this signature were assessed in the Metascape database, gene set enrichment analysis (GSEA) and immune infiltration analysis.ResultsRRA analysis identified 40 differently expressed genes (DEGs). A 3-gene risk score signature (MKQ signature) associated with disease-free survival (DFS) was generated. DFS was significantly longer in CC patients with lower than higher scores (P=0.0046). The areas under curves (AUCs) of the time-dependent ROC curves of MKQ signature at 1-, 3- and 5-year DFS were 1, 0.963 and 0.961 respectively. Recurrence-free survival (RFS) was significantly longer in patients in {"type":"entrez-geo","attrs":{"text":"GSE39582","term_id":"39582"}}GSE39582 with lower than higher risk scores (P=0.032). The AUCs for 1-, 3- and 5-year RFS in {"type":"entrez-geo","attrs":{"text":"GSE39582","term_id":"39582"}}GSE39582 were 0.63, 0.618 and 0.583, respectively, validating the value of the MKQ signature. Functional annotation and GSEA revealed that the MKQ signature was associated with multiple immune-related pathways. Immune cell infiltration was found to differ in patients differing in the MKQ signature.ConclusionsGene expression and microsatellite status identified a 3-gene signature (MKQ signature) that could facilitate risk-stratified management in patients with stage II CC. Dysregulation of MSMB, KRT23, and QPRT can serve as prognostic markers in stage II CC.     

Low intratumoral genetic neutrophil-to-lymphocyte ratio (NLR) is associated with favorable tumor immune microenvironment and with survival in triple negative breast cancer (TNBC)

Patients with triple negative breast cancer (TNBC) have a poor prognosis. A novel prognostic biomarker may guide management by appropriately selecting patients for particular treatments. Peripheral blood neutrophil-to-lymphocyte ratio (NLR) was reported to associate with cancer progression, thus we hypothesized that intratumor genetic NLR will reflect tumor immune microenvironment (TIME) and breast cancer biology. The intratumoral genetic NLR previously defined as the ratio of CD66b (CEACAM8) and CD8 (CD8A) gene expressions was utilized to analyze total of 2,994 patients from METABRIC, TCGA, {"type":"entrez-geo","attrs":{"text":"GSE21094","term_id":"21094"}}GSE21094, {"type":"entrez-geo","attrs":{"text":"GSE22358","term_id":"22358"}}GSE22358, {"type":"entrez-geo","attrs":{"text":"GSE25088","term_id":"25088"}}GSE25088, {"type":"entrez-geo","attrs":{"text":"GSE32646","term_id":"32646"}}GSE32646, and {"type":"entrez-geo","attrs":{"text":"GSE2603","term_id":"2603"}}GSE2603 cohorts. Intratumoral genetic NLR did not correlate with cancer stage nor clinical parameters of cancer cell proliferation such as Nottingham histological grade or MKI67 expression levels in neither the METABRIC or TCGA cohorts. Intratumoral genetic NLR-high breast cancer was not associated with pathologic complete response (pCR) after neoadjuvant chemotherapy in 5 independent cohorts with different regimens. Despite these results, intratumoral genetic NLR-high TNBC demonstrated worse disease-free, disease-specific, and overall survival. Intratumoral genetic NLR-low TNBC enriched multiple immune-related gene sets, was associated with higher favorable immune-related scores and with a favorable TIME, whereas no gene sets enriched to NLR-high TNBC. In conclusion, intratumoral genetic NLR-low TNBC was associated with favorable TIME and with better survival.     

Anti-malarial drug dihydroartemisinin downregulates the expression levels of CDK1 and CCNB1 in liver cancer

Liver cancer is the third leading cause of cancer-associated mortality worldwide. By the time liver cancer is diagnosed, it is already in the advanced stage. Therefore, novel therapeutic strategies need to be identified to improve the prognosis of patients with liver cancer. In the present study, the profiles of {"type":"entrez-geo","attrs":{"text":"GSE84402","term_id":"84402"}}GSE84402, {"type":"entrez-geo","attrs":{"text":"GSE19665","term_id":"19665"}}GSE19665 and {"type":"entrez-geo","attrs":{"text":"GSE121248","term_id":"121248"}}GSE121248 were used to screen differentially expressed genes (DEGs). Subsequently, Gene Ontology and Kyoto Encyclopedia of Genes and Genomes pathway enrichment analyses for DEGs were conducted using the Database for Annotation, Visualization and Integrated Discovery. The protein-protein interaction network was established to screen the hub genes associated with liver cancer. Additionally, the expression levels of hub genes were validated using the Gene Expression Profiling Interactive Analysis and Oncomine databases. In addition, the prognostic value of hub genes in patients with liver cancer was analyzed using Kaplan-Meier Plotter. It was demonstrated that 132 and 246 genes were upregulated and downregulated, respectively, in patients with liver cancer. Among these DEGs, 10 hub genes with high connected node values were identified, which were AURKA, BIRC5, BUB1B, CCNA2, CCNB1, CCNB2, CDC20, CDK1, DLGAP5 and MAD2L1. CDK1 and CCNB1 had the most connection nodes and the highest score and were therefore, the most significantly expressed. In addition, it was demonstrated that high expression levels of CDK1 and CCNB1 were associated with poor overall survival time of patients with liver cancer. Dihydroartemisinin (DHA) is a Food and Drug Administration-approved drug, which is derived from the traditional Chinese medicine Artemisia annua Linn. DHA inhibits cell proliferation in numerous cancer types, including liver cancer. In our previous study, it was revealed that DHA inhibited the proliferation of HepG2215 cells. In the present study, it was further demonstrated that DHA reduced the expression levels of CDK1 and CCNB1 in liver cancer. Overall, CDK1 and CCNB1 were the potential therapeutic targets of liver cancer, and DHA reduced the expression levels of CDK1 and CCNB1, and inhibited the proliferation of liver cancer cells.     

Establishment and validation of an autophagy-related prognostic signature for survival predicting in cutaneous melanoma

Existing staging system for prognosis evaluating for Skin Cutaneous Melanoma (SKCM) patients had defects of subjective, inaccuracy and inconsistently, therefore, to identify specific and applicable prognostic markers and promote personalized therapeutic interventions is urgently required. This study aims to build a robust autophagy-related genes (ARGs) signature for prognosis monitoring of SKCM patients. We determined 26 ARGs as differentially expressed autophagy-related genes (DEARGs) from 103 SKCM and 23 normal skin samples in {"type":"entrez-geo","attrs":{"text":"GSE15605","term_id":"15605"}}GSE15605 and {"type":"entrez-geo","attrs":{"text":"GSE3189","term_id":"3189"}}GSE3189 datasets. Optimal prognostic DEARGs composed the risk model were screened and verified in 458 SKCM patients in TCGA cohort as the training cohort and 209 patients in {"type":"entrez-geo","attrs":{"text":"GSE65904","term_id":"65904"}}GSE65904 as the test cohort. Finally, 4 optimal independent prognostic DEARGs (CAPNS1, DAPK2, PARP1 and PTK6) were filtered out in the training cohort to establish the risk model. A prognostic nomogram was established for quantitative survival prediction. The risk model grouped high-risk SKCM cancer patients exhibited significantly shorter survival times in both training and test cohorts. The area under the ROC curve for risk score model was 0.788 and 0.627 in the training and test cohorts indicated the risk model was relatively accurate for prognosis monitoring. Clinical correlation analysis exhibited that risk score was an independent predictor for prognosis significantly associated with T/N classification. The prognostic value of the 4 risk genes formed the risk model was also validated respectively. We identified a novel autophagy-related signature for prognosis monitoring. It has the potential to be an independent prognostic indicator and can benefit targeted therapy.     

A 5-gene classifier from the carcinoma-associated fibroblast transcriptomic profile and clinical outcome in colorectal cancer

Based on 108 differentially expressed genes between carcinoma-associated fibroblasts (CAFs) and paired normal colonic fibroblasts we recently reported, a 5-gene classifier for relapse prediction in Stage II/III colorectal cancer (CRC) was developed. Its predictive value was validated in datasets {"type":"entrez-geo","attrs":{"text":"GSE17538","term_id":"17538"}}GSE17538, {"type":"entrez-geo","attrs":{"text":"GSE33113","term_id":"33113"}}GSE33113 and {"type":"entrez-geo","attrs":{"text":"GSE14095","term_id":"14095"}}GSE14095. An additional validation was performed in a metacohort (n=317) and 142 CRC patients by means of RT-PCR. The 5-gene classifier was significantly associated with increased relapse risk and death from CRC across all validation series of Stage II/III patients used. Multivariate Cox regression analyses confirmed the independent prognostic value of the stromal classifier (HR=2.67; P=0.002). Post-test probabilities provided evidence of the suitability of the 5-gene classifier in clinical practice, identifying a subgroup of Stage-II patients who were at high risk of relapse. Moreover, the a priory worst prognosis mesenchymal subtype of tumours can be stratified according to the physiological status of their carcinoma-associated fibroblasts. In conclusion the CAFs-derived 5-gene classifier provides more accurate information about outcome than conventional clinicopathological criteria and it could be useful to take clinical decisions, especially in Stage II. Additionally, the classifier put into relevance the CAF''s intratumoral heterogeneity and might contribute to find relevant targets for depleting adequate CAFS subtypes.     

Prostate cancer genes associated with TMPRSS2�CERG gene fusion and prognostic of biochemical recurrence in multiple cohorts

Background:

Recent studies have indicated that prostate cancer patients with the TMPRSS2–ERG gene fusion have a higher risk of recurrence. To identify markers associated with TMPRSS2–ERG fusion and prognostic of biochemical recurrence, we analysed a cohort of 139 men with prostate cancer for 502 molecular markers.

Methods:

RNA from radical prostatectomy tumour specimens was analysed using cDNA-mediated, annealing, selection, extension and ligation (DASL) to determine mRNAs associated with TMPRSS2–ERG T1/E4 fusion and prognostic of biochemical recurrence. Differentially expressed mRNAs in T1/E4-positive tumours were determined using significance analysis of microarrays (false discovery rate (FDR) <5%). Univariate and multivariate Cox regression determined genes, gene signatures and clinical factors prognostic of recurrence (P-value <0.05, log–rank test). Analysis of two prostate microarray studies ({"type":"entrez-geo","attrs":{"text":"GSE1065","term_id":"1065"}}GSE1065 and {"type":"entrez-geo","attrs":{"text":"GSE8402","term_id":"8402"}}GSE8402) validated the findings.

Results:

In the 139 patients from this study and from a 455-patient Swedish cohort, 15 genes in common were differentially regulated in T1/E4 fusion-positive tumours (FDR <0.05). The most significant mRNAs in both cohorts coded ERG. Nine genes were found prognostic of recurrence in this study and in a 596-patient Minnesota cohort. A molecular recurrence score was significant in prognosticating recurrence (P-value 0.000167) and remained significant in multivariate analysis of a mixed clinical model considering Gleason score and TMPRSS2–ERG fusion status.

Conclusions:

TMPRSS2–ERG T1/E4 fusion-positive tumours had differentially regulated mRNAs observed in multiple studies, the most significant one coded for ERG. Several mRNAs were consistently associated with biochemical recurrence and have potential clinical utility but will require further validation for successful translation.