RecurrenceOnline: an online analysis tool to determine breast cancer recurrence and hormone receptor status using microarray data

Tripartite-motif 3 represses ovarian cancer progression by downregulating lactate dehydrogenase A and inhibiting AKT signaling

The E3 ubiquitin ligase Tripartite-motif 3 (TRIM3) is known to play a crucial role in tumor suppression in various tumors through different mechanisms. However, its function and mechanism in ovarian cancer have yet to be elucidated. Our study aims to investigate the expression of TRIM3 in ovarian cancer and evaluate its role in the development of the disease. Our findings revealed a significant decrease in TRIM3 mRNA and protein levels in ovarian cancer tissues and cells when compared to normal ovarian epithelial tissues and cells. Furthermore, we observed a negative correlation between the protein level of TRIM3 and the FIGO stage, as well as a positive correlation with the survival of ovarian cancer patients. Using gain and loss of function experiments, we demonstrated that TRIM3 can inhibit cell proliferation, migration and invasion of the ovarian cancer cells in vitro, as well as suppress tumor growth in vivo. Mechanistic studies showed that TRIM3 interacts with lactate dehydrogenase A, a key enzyme in the glycolytic pathway, through its B-box and coiled-coil domains and induces its ubiquitination and proteasomal degradation, leading to the inhibition of glycolytic ability in ovarian cancer cells. RNA-sequencing analysis revealed significant alterations in the phosphatidylinositol signaling pathways upon TRIM3 overexpression. Additionally, overexpression of TRIM3 inhibited the phosphorylation of AKT. In conclusion, our study demonstrated that TRIM3 exerts a tumor-suppressive effect in ovarian cancer, at least partially, by downregulating LDHA and inhibiting the AKT signaling pathway, and thus leading to the inhibition of glycolysis and limiting the growth of ovarian cancer cells.

Suppression of Bcl3 Disrupts Viability of Breast Cancer Cells through Both p53-Dependent and p53-Independent Mechanisms via Loss of NF-κB Signalling

The NF-κB co-factor Bcl3 is a proto-oncogene that promotes breast cancer proliferation, metastasis and therapeutic resistance, yet its role in breast cancer cell survival is unclear. Here, we sought to determine the effect of Bcl3 suppression alone on breast cancer cell viability, with a view to informing future studies that aim to target Bcl3 therapeutically. Bcl3 was suppressed by siRNA in breast cancer cell lines before changes in viability, proliferation, apoptosis and senescence were examined. Bcl3 suppression significantly reduced viability and was shown to induce apoptosis in all cell lines tested, while an additional p53-dependent senescence and senescence-associated secretory phenotype was also observed in those cells with functional p53. The role of the Bcl3/NF-κB axis in this senescence response was confirmed via siRNA of the non-canonical NF-κB subunit NFKB2/p52, which resulted in increased cellular senescence and the canonical subunit NFKB1/p50, which induced the senescence-associated secretory phenotype. An analysis of clinical data showed a correlation between reduced relapse-free survival in patients that expressed high levels of Bcl3 and carried a p53 mutation. Together, these data demonstrate a dual role for Bcl3/NF-κB in the maintenance of breast cancer cell viability and suggests that targeting Bcl3 may be more beneficial to patients with tumours that lack functional p53.

The Breast Cancer Single-Cell Atlas: Defining cellular heterogeneity within model cell lines and primary tumors to inform disease subtype, stemness, and treatment options

PURPOSE

Breast Cancer (BC) is the most diagnosed cancer in women; however, through significant research, relative survival rates have significantly improved. Despite progress, there remains a gap in our understanding of BC subtypes and personalized treatments. This manuscript characterized cellular heterogeneity in BC cell lines through scRNAseq to resolve variability in subtyping, disease modeling potential, and therapeutic targeting predictions.

METHODS

We generated a Breast Cancer Single-Cell Cell Line Atlas (BSCLA) to help inform future BC research. We sequenced over 36,195 cells composed of 13 cell lines spanning the spectrum of clinical BC subtypes and leveraged publicly available data comprising 39,214 cells from 26 primary tumors.

RESULTS

Unsupervised clustering identified 49 subpopulations within the cell line dataset. We resolve ambiguity in subtype annotation comparing expression of Estrogen Receptor, Progesterone Receptor, and Human Epidermal Growth Factor Receptor 2 genes. Gene correlations with disease subtype highlighted S100A7 and MUCL1 overexpression in HER2 + cells as possible cell motility and localization drivers. We also present genes driving populational drifts to generate novel gene vectors characterizing each subpopulation. A global Cancer Stem Cell (CSC) scoring vector was used to identify stemness potential for subpopulations and model multi-potency. Finally, we overlay the BSCLA dataset with FDA-approved targets to identify to predict the efficacy of subpopulation-specific therapies.

CONCLUSION

The BSCLA defines the heterogeneity within BC cell lines, enhancing our overall understanding of BC cellular diversity to guide future BC research, including model cell line selection, unintended sample source effects, stemness factors between cell lines, and cell type-specific treatment response.

Bioinformatics Analysis and Validation of the Role of Lnc-RAB11B-AS1 in the Development and Prognosis of Hepatocellular Carcinoma

Lnc-RAB11B-AS1 is reported to be dysregulated in several types of cancers and can function as both an oncogene and tumor suppressor gene. To evaluate the potential role of lnc-RAB11B-AS1 in hepatocellular carcinoma (HCC), we investigated and evaluated its expression in HCC based on the data mining of a series of public databases, including TCGA, GEO, ICGC, HPA, DAVID, cBioPortal, GeneMIANA, TIMER, and ENCORI. The data showed downregulation of lnc-RAB11B-AS1 in HCC and was accompanied by the synchronous downregulation of the targeted RAB11B mRNA and its protein. Low expression of lnc-RAB11B-AS1 was associated with shorter overall survival (OS) and disease-free survival (DFS) of HCC patients, PD1/PD-L1 was correlated with low expression of RAB11B. Furthermore, Gene Ontology (GO) functional annotation and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis showed a correlation between immune cell change and non-alcoholic fatty liver disease. The above findings revealed that lnc-RAB11B-AS1 was down-regulated in HCC and closely associated with the clinical stage of the HCC patients, suggesting that lnc-RAB11B-AS1 could be a possible predictor for HCC and a potential new therapeutic target for the treatment of HCC.

GIT1 protects against breast cancer growth through negative regulation of Notch

Hyperactive Notch signalling is frequently observed in breast cancer and correlates with poor prognosis. However, relatively few mutations in the core Notch signalling pathway have been identified in breast cancer, suggesting that as yet unknown mechanisms increase Notch activity. Here we show that increased expression levels of GIT1 correlate with high relapse-free survival in oestrogen receptor-negative (ER(-)) breast cancer patients and that GIT1 mediates negative regulation of Notch. GIT1 knockdown in ER(-) breast tumour cells increased signalling downstream of Notch and activity of aldehyde dehydrogenase, a predictor of poor clinical outcome. GIT1 interacts with the Notch intracellular domain (ICD) and influences signalling by inhibiting the cytoplasm-to-nucleus transport of the Notch ICD. In xenograft experiments, overexpression of GIT1 in ER(-) cells prevented or reduced Notch-driven tumour formation. These results identify GIT1 as a modulator of Notch signalling and a guardian against breast cancer growth.

Notch signalling is reported to be hyperactivated in oestrogen receptor-negative (ER-) breast cancer. Here the authors show that G protein-coupled receptor kinase-interacting protein 1 (GIT1) negatively regulates Notch signalling and tumour growth in ER- breast cancer by blocking Notch ICD nuclear translocation.

Identification of epigenetic modifications mediating the antagonistic effect of selenium against cadmium-induced breast carcinogenesis

The antagonistic effect of selenium (Se) against cadmium (Cd)-induced breast carcinogenesis was reported, but underlying mechanisms were unclear. The aim of this study was to identify the epigenetically regulated genes and biological pathways mediating the antagonistic effect. We exposed MCF-7 cells to Cd and Se alone or simultaneously. Cell proliferation was assessed by MTT assay, and differential epigenome (DNA methylation, microRNA, and long non-coding RNA) was obtained by microarrays. We cross-verified the epigenetic markers with differential transcriptome, and the ones modulated by Cd and Se in opposite directions were regarded to mediate the antagonistic effect. The epigenetically regulated genes were validated by using gene expression data in human breast tissues. We further assessed the biological functions of these validated genes. Our results showed that Se alleviated the proliferative effect of Cd on MCF-7 cell. A total of 10 epigenetically regulated genes were regarded to mediate the antagonistic effect, including APBA2, KIAA0895, DHX35, CPEB3, SVIL, MYLK, ZFYVE28, ABLIM2, GRB10, and PCDH9. Biological function analyses suggested that these epigenetically regulated genes were involved in multiple cancer-related pathways, such as focal adhesion and PI3K/Akt pathway. In conclusion, we provided evidence that Se antagonized the Cd-induced breast carcinogenesis via epigenetic modification and revealed the critical pathways.

gp130 Cytokines Activate Novel Signaling Pathways and Alter Bone Dissemination in ER+ Breast Cancer Cells

Breast cancer cells frequently home to the bone marrow, where they encounter signals that promote survival and quiescence or stimulate their proliferation. The interleukin-6 (IL-6) cytokines signal through the co-receptor glycoprotein130 (gp130) and are abundantly secreted within the bone microenvironment. Breast cancer cell expression of leukemia inhibitory factor (LIF) receptor (LIFR)/STAT3 signaling promotes tumor dormancy in the bone, but it is unclear which, if any of the cytokines that signal through LIFR, including LIF, oncostatin M (OSM), and ciliary neurotrophic factor (CNTF), promote tumor dormancy and which signaling pathways are induced. We first confirmed that LIF, OSM, and CNTF and their receptor components were expressed across a panel of breast cancer cell lines, although expression was lower in estrogen receptor-negative (ER^- ) bone metastatic clones compared with parental cell lines. In estrogen receptor-positive (ER⁺ ) cells, OSM robustly stimulated phosphorylation of known gp130 signaling targets STAT3, ERK, and AKT, while CNTF activated STAT3 signaling. In ER^- breast cancer cells, OSM alone stimulated AKT and ERK signaling. Overexpression of OSM, but not CNTF, reduced dormancy gene expression and increased ER⁺ breast cancer bone dissemination. Reverse-phase protein array revealed distinct and overlapping pathways stimulated by OSM, LIF, and CNTF with known roles in breast cancer progression and metastasis. In breast cancer patients, downregulation of the cytokines or receptors was associated with reduced relapse-free survival, but OSM was significantly elevated in patients with invasive disease and distant metastasis. Together these data indicate that the gp130 cytokines induce multiple signaling cascades in breast cancer cells, with a potential pro-tumorigenic role for OSM and pro-dormancy role for CNTF. © 2021 American Society for Bone and Mineral Research (ASBMR).

The cell cycle-related genes RHAMM, AURKA, TPX2, PLK1, and PLK4 are associated with the poor prognosis of breast cancer patients

Breast cancer is the third most common type of cancer diagnosed. Cell cycle is a complex but highly organized and controlled process, in which normal cells sense mitogenic growth signals that instruct them to enter and progress through their cell cycle. This process culminates in cell division generating two daughter cells with identical amounts of genetic material. Uncontrolled proliferation is one of the hallmarks of cancer. In this study, we analyzed the expression of the cell cycle-related genes receptor for hyaluronan (HA)-mediated motility (RHAMM), AURKA, TPX2, PLK1, and PLK4 and correlated them with the prognosis in a collective of 3952 breast cancer patients. A high messenger RNA expression of all studied genes correlated with a poor prognosis. Stratifying the patients according to the expression of hormonal receptors, we found that in patients with estrogen and progesterone receptor-positive and human epithelial growth factor receptor 2-negative tumors, and Luminal A and Luminal B tumors, the expression of the five analyzed genes correlates with worse survival. qPCR analysis of a panel of breast cancer cell lines representative of major molecular subtypes indicated a predominant expression in the luminal subtype. In vitro experiments showed that radiation influences the expression of the five analyzed genes both in luminal and triple-negative model cell lines. Functional analysis of MDA-MB-231 cells showed that small interfering RNA knockdown of PLK4 and TPX2 and pharmacological inhibition of PLK1 had an impact on the cell cycle and colony formation. Looking for a potential upstream regulation by microRNAs, we observed a differential expression of RHAMM, AURKA, TPX2, PLK1, and PLK4 after transfecting the MDA-MB-231 cells with three different microRNAs. Survival analysis of miR-34c-5p, miR-375, and miR-142-3p showed a different impact on the prognosis of breast cancer patients. Our study suggests that RHAMM, AURKA, TPX2, PLK1, and PLK4 can be used as potential targets for treatment or as a prognostic value in breast cancer patients.

Role of Rho guanine nucleotide exchange factors in non-small cell lung cancer

Conventionally, Rho guanine nucleotide exchange factors (GEFs) are known activators of Rho guanosine triphosphatases (GTPases) that promote tumorigenesis. However, the role of Rho GEFs in non-small cell lung cancer (NSCLC) remains largely unknown. Through the screening of 81 Rho GEFs for their expression profiles and correlations with survival, four of them were identified with strong significance for predicting the prognosis of NSCLC patients. The four Rho GEFs, namely ABR, PREX1, DOCK2 and DOCK4, were downregulated in NSCLC tissues compared to normal tissues. The downregulation of ABR, PREX1, DOCK2 and DOCK4, which can be attributfed to promoter methylation, is correlated with poor prognosis. The underexpression of the four key Rho GEFs might be related to the upregulation of MYC signaling and DNA repair pathways, leading to carcinogenesis and poor prognosis. Moreover, overexpression of ABR was shown to have a tumor-suppressive effect in PC9 and H1703 cells. In conclusion, the data reveal the unprecedented role of ABR as tumor suppressor in NSCLC. The previously unnoticed functions of Rho GEFs in NSCLC will inspire researchers to investigate the distinct roles of Rho GEFs in cancers, in order to provide critical strategies in clinical practice.

microRNA-140-3p modulates invasiveness, motility, and extracellular matrix adhesion of breast cancer cells by targeting syndecan-4

Syndecan-4, a predicted target of the microRNA miR-140-3p, plays an important role in multiple steps of tumor progression and is the second most abundant heparan sulfate proteoglycan produced by breast carcinoma cell lines. To investigate the potential functional relationship of miR-140-3p and syndecan-4, MDA-MB-231, SKBR3, and MCF-7 breast cancer (BC) cells were transiently transfected with pre-miR-140-3p, syndecan-4 small interfering RNAJ, or control reagents, respectively. Altered cell behavior was monitored by adhesion, migration, and invasion chamber assays. Moreover, the prognostic value of syndecan-4 was assessed by Kaplan-Maier Plotter analysis of gene expression data from tumor samples of 4929 patients. High expression of syndecan-4 was associated with better relapse-free survival in the whole collective of BC patients, but correlated with a worse survival in the subgroup of estrogen receptor negative and estrogen/progesterone-receptor negative patients. miR-140-3p expression was associated with improved survival irrespective of hormone receptor status. miR-140-3p overexpression induced posttranscriptional downregulation of syndecan-4, as demonstrated by quantitative real-time PCR (qPCR), flow cytometry, and luciferase assays, resulting in decreased BC cell migration and matrigel invasiveness. Furthermore, miR-140-3p overexpression and syndecan-4 silencing increased the adhesion of BC to fibronectin and laminin. qPCR analysis demonstrated that syndecan-4 silencing leads to altered gene expression of adhesion-related molecules, such as fibronectin and focal adhesion kinase, as well as in the gene expression of the proinvasive factors matrix metalloproteinase 2 and heparanase (also known as HPSE). We conclude that syndecan-4 is a novel target of miR-140-3p that regulates BC cell invasiveness and cell-matrix interactions in the tumor microenvironment.

A 95-gene signature stratifies recurrence risk of invasive disease in ER-positive, HER2-negative, node-negative breast cancer with intermediate 21-gene signature recurrence scores

PURPOSE

A subset of patients with intermediate 21-gene signature assay recurrence score may benefit from adjuvant chemoendocrine therapy, but a predictive strategy is needed to identify such patients. The 95-gene signature assay was tested to stratify patients with intermediate RS into high (95GC-H) and low (95GC-L) groups that were associated with invasive recurrence risk.

METHODS

Patients with ER-positive, HER2-negative, node-negative breast cancer and RS 11-25 who underwent definitive surgery and adjuvant endocrine therapy without any cytotoxic agents were included. RNA was extracted from archived formalin-fixed, paraffin-embedded samples, and 95-gene signature was calculated.

RESULTS

206 patients had RS of 11-25 (95GC-L, N = 163; 95GC-H, N = 43). In Cox proportional hazards model, 95GC-H was significantly associated with shorter time to recurrence than was 95GC-L (HR 5.94; 95%CI 1.81-19.53; P = 0.005). The correlation between 95-gene signature and 21-gene signature assay scores was not strong (correlation coefficient r = 0.27), which might suggest that 95-gene signature reflects biological characteristics differing from what 21-gene signature shows.

CONCLUSIONS

The 95-gene signature stratifies patients with ER-positive, HER2-negative, node-negative invasive breast cancer and intermediate RS of 11-25 into high and low groups that are associated with recurrence risk of invasive disease. Further retrospective analysis in the prospectively accrued TAILORx population is warranted to confirm that 95-gene signature can identify patients who would benefit from adjuvant chemoendocrine therapy.

Clinical implications of TPO and AOX1 in pediatric papillary thyroid carcinoma

BACKGROUND

Thyroid carcinoma is a common pediatric head and neck cancer, of which papillary thyroid cancer (PTC) is the most common type. Previously, we found that thyroid peroxidase (TPO) and aldehyde oxidase 1 (AOX1) were differentially expressed in PTC. This study explored the clinical importance of TPO and AOX1 in the diagnosis and prognosis of PTC in children.

METHODS

Both TPO and AOX1 expression in PTC were analyzed using datasets from Gene Expression Omnibus (GEO). TPO and AOX1 protein levels in plasma from patients with PTC and non-tumor controls were detected via enzyme-linked immunosorbent assay (ELISA). The diagnostic accuracy of TPO and AOX1 was assessed using receiver operating characteristic (ROC) curve analysis. The association between gene expression levels and patient survival was explored using the Kaplan-Meier plotter online database.

RESULTS

The results revealed that TPO and AOX1 expression was significantly downregulated in four independent datasets (GSE33630, GSE27155, GSE3678, and GSE3467). TPO and AOX1 protein levels in blood plasma were significantly decreased in patients with PTC. Quantitative analysis demonstrated that TPO and AOX1 levels in plasma had satisfactory predictive performance and the ability to discriminate PTC from healthy samples. Prognostic analysis demonstrated that low levels of TPO and AOX1 were markedly associated with poor survival in patients with PTC.

CONCLUSIONS

In summary, these results implied that TPO and AOX1 could serve as novel biomarkers for the diagnosis and prognosis of pediatric PTC.

Integrated bioinformatic analysis identifies UBE2Q1 as a potential prognostic marker for high grade serous ovarian cancer

Background

High grade serous ovarian cancer (HGSOC) accounts for nearly 60% of total cases of epithelial ovarian cancer. It is the most aggressive subtype, which shows poor prognosis and low patient survival. For better management of HGSOC patients, new prognostic biomarkers are required to facilitate improved treatment strategies and ensure suitable healthcare decisions.

Methods

We performed genome wide expression analysis of HGSOC patient samples to identify differentially expressed genes (DEGs) using R based Limma package, Clust and other statistical tools. The identified DEGs were subjected to weighted gene co-expression network analysis (WGCNA) to identify co-expression patterns of relevant genes. Module trait and gene ontology analyses were performed to establish important gene co-expression networks and their biological functions. Overlapping the most relevant DEG cluster 4 with prominent WGCNA cyan module identified strongest correlation of UBE2Q1 with ovarian cancer and its prognostic significance on survival probability of ovarian cancer patients was investigated. The predictive value of UBE2Q1 as a potential biomarker was analysed by correlating its expression with 12-months relapse free survival of patients in response to platin/taxane, the standard first-line chemotherapy for ovarian cancer, and analysing area under the ROC curve.

Results

An integrated gene expression analysis and WGCNA, identified UBE2Q1 as a potential prognostic marker associated with poor relapse-free survival and response outcome to platin/taxane treatment of patients with high grade serous ovarian cancer.

Conclusions

Our study identifies a potential UBE2Q1 – B4GALT3 functional axis in ovarian cancer, where only the E2 conjugating enzyme showed a poor prognostic impact on the disease.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12885-021-07928-z.

Differential epigenetic profiles induced by sodium selenite in breast cancer cells

OBJECTIVES

Selenium (Se) was a potential anticancer micronutrient with proposed epigenetic effect. However, the Se-induced epigenome in breast cancer cells was yet to be studied.

METHODS

The profiles of DNA methylation, microRNA (miRNA), long non-coding RNA (lncRNA), and message RNA (mRNA) in breast cancer cells treated with sodium selenite were examined by microarrays. We verified the epigenetic modifications by integrating their predicted target genes and differentially expressed mRNAs. The epigenetically regulated genes were further validated in a breast cancer cohort by associating with tumor progression. We conducted a series of bioinformatics analyses to assess the biological function of these validated genes and identified the critical genes.

RESULTS

The Se-induced epigenome regulated the expression of 959 genes, and 349 of them were further validated in the breast cancer cohort. Biological function analyses suggested that these validated genes were enriched in several cancer-related pathways, such as PI3K/Akt and metabolic pathways. Based on the degrees of expression change, hazard ratio difference, and connectivity, NEDD4L and FMO5 were identified as the critical genes.

CONCLUSIONS

These results confirmed the epigenetic effects of sodium selenite and revealed the epigenetic profiles in breast cancer cells, which would help understand the mechanisms of Se against breast cancer.

Initial Identification of UDP-Glucose Dehydrogenase as a Prognostic Marker in Breast Cancer Patients, Which Facilitates Epirubicin Resistance and Regulates Hyaluronan Synthesis in MDA-MB-231 Cells

UDP-glucose-dehydrogenase (UGDH) synthesizes UDP-glucuronic acid. It is involved in epirubicin detoxification and hyaluronan synthesis. This work aimed to evaluate the effect of UGDH knockdown on epirubicin response and hyaluronan metabolism in MDA-MB-231 breast cancer cells. Additionally, the aim was to determine UGDH as a possible prognosis marker in breast cancer. We studied UGDH expression in tumors and adjacent tissue from breast cancer patients. The prognostic value of UGDH was studied using a public Kaplan–Meier plotter. MDA-MB-231 cells were knocked-down for UGDH and treated with epirubicin. Epirubicin-accumulation and apoptosis were analyzed by flow cytometry. Hyaluronan-coated matrix and metabolism were determined. Autophagic-LC3-II was studied by Western blot and confocal microscopy. Epirubicin accumulation increased and apoptosis decreased during UGDH knockdown. Hyaluronan-coated matrix increased and a positive modulation of autophagy was detected. Higher levels of UGDH were correlated with worse prognosis in triple-negative breast cancer patients that received chemotherapy. High expression of UGDH was found in tumoral tissue from HER2^--patients. However, UGDH knockdown contributes to epirubicin resistance, which might be associated with increases in the expression, deposition and catabolism of hyaluronan. The results obtained allowed us to propose UGDH as a new prognostic marker in breast cancer, positively associated with development of epirubicin resistance and modulation of extracellular matrix.

Intratumoral heterogeneity of second-harmonic generation scattering from tumor collagen and its effects on metastatic risk prediction

Background

Metastases are the leading cause of breast cancer-related deaths. The tumor microenvironment impacts cancer progression and metastatic ability. Fibrillar collagen, a major extracellular matrix component, can be studied using the light scattering phenomenon known as second-harmonic generation (SHG). The ratio of forward- to backward-scattered SHG photons (F/B) is sensitive to collagen fiber internal structure and has been shown to be an independent prognostic indicator of metastasis-free survival time (MFS). Here we assess the effects of heterogeneity in the tumor matrix on the possible use of F/B as a prognostic tool.

Methods

SHG imaging was performed on sectioned primary tumor excisions from 95 untreated, estrogen receptor-positive, lymph node negative invasive ductal carcinoma patients. We identified two distinct regions whose collagen displayed different average F/B values, indicative of spatial heterogeneity: the cellular tumor bulk and surrounding tumor-stroma interface. To evaluate the impact of heterogeneity on F/B’s prognostic ability, we performed SHG imaging in the tumor bulk and tumor-stroma interface, calculated a 21-gene recurrence score (surrogate for OncotypeDX®, or S-ODX) for each patient and evaluated their combined prognostic ability.

Results

We found that F/B measured in tumor-stroma interface, but not tumor bulk, is prognostic of MFS using three methods to select pixels for analysis: an intensity threshold selected by a blinded observer, a histogram-based thresholding method, and an adaptive thresholding method. Using both regression trees and Random Survival Forests for MFS outcome, we obtained data-driven prediction rules that show F/B from tumor-stroma interface, but not tumor bulk, and S-ODX both contribute to predicting MFS in this patient cohort. We also separated patients into low-intermediate (S-ODX < 26) and high risk (S-ODX ≥26) groups. In the low-intermediate risk group, comprised of patients not typically recommended for adjuvant chemotherapy, we find that F/B from the tumor-stroma interface is prognostic of MFS and can identify a patient cohort with poor outcomes.

Conclusions

These data demonstrate that intratumoral heterogeneity in F/B values can play an important role in its possible use as a prognostic marker, and that F/B from tumor-stroma interface of primary tumor excisions may provide useful information to stratify patients by metastatic risk.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12885-020-07713-4.

Differential epigenetic and transcriptional profile in MCF-7 breast cancer cells exposed to cadmium

Cadmium (Cd) has been confirmed to be associated with breast carcinogenesis, but the mechanism was not clarified yet. Given that epigenetic modification was speculated as underlying mechanism, we examined the differential epigenome caused by Cd in breast cancer cells. Profiles of DNA methylation, microRNA (miRNA), long non-coding RNA (lncRNA), and message RNA (mRNA) were derived from Cd-treated and untreated MCF-7 breast cancer cells by microarray. We identified 997 target genes epigenetically regulated by Cd through cross-verification with the differential epigenome and transcriptome, and 400 of them were further validated in a breast cancer cohort. Biological function analyses suggested that several pathways were involved in Cd-induced breast carcinogenesis, such as Wnt signaling, metabolism, and human papilloma virus (HPV) infection. TXNRD1 and CCT3 were further identified as the critical genes based on the degree of expression change, hazard ratio difference, and connectivity. The present study revealed that Cd epigenetically regulated several pathways involving in breast carcinogenesis, particularly the Wnt signaling and metabolic pathways, among which TXNRD1 and CCT3 might play critical roles. It was also suggested that Cd and HPV infection might jointly participate in breast tumorigenesis.

A Combined Systemic Strategy for Overcoming Cisplatin Resistance in Head and Neck Cancer: From Target Identification to Drug Discovery

Simple Summary

The efficiency of cisplatin is limited by drug resistance in head–neck cancer (HNC) patients. In this study, we established a cisplatin resistance (CR) cell model, generated CR related transcriptome profiling, and combined application of bioinformatics methodology to discover a possible way to overcome CR. Analysis of the functional pathway revealed that mitotic division is a novel mechanism significantly contributing to CR. Spindle pole body component 25 (SPC25), a kinetochore protein, was overexpressed in CR cells and significantly correlated with worse HNC patient survival. The silencing of SPC25 increased cisplatin sensitivity and reduced cancer stemness property. Integration of CR transcriptome profiling and drug database discovered a natural extract compound, celastrol, possessing a potent cytotoxic effect in CR cells to reverse CR. Thus, we combined systemic strategies to demonstrated that a novel biological process (mitotic cell division), a hub gene (SPC25), and a natural compound (celastrol) as novel strategies for the treatment of refractory HNC.

Abstract

Cisplatin is the first-line chemotherapy agent for head and neck cancer (HNC), but its therapeutic effects are hampered by its resistance. In this study, we employed systemic strategies to overcome cisplatin resistance (CR) in HNC. CR cells derived from isogenic HNC cell lines were generated. The CR related hub genes, functional mechanisms, and the sensitizing candidates were globally investigated by transcriptomic and bioinformatic analyses. Clinically, the prognostic significance was assessed by the Kaplan–Meier method. Cellular and molecular techniques, including cell viability assay, tumorsphere formation assay, RT-qPCR, and immunoblot, were used. Results showed that these CR cells possessed highly invasive and stem-like properties. A total of 647 molecules was identified, and the mitotic division exhibited a novel functional mechanism significantly related to CR. A panel of signature molecules, MSRB3, RHEB, ULBP1, and spindle pole body component 25 (SPC25), was found to correlate with poor prognosis in HNC patients. SPC25 was further shown as a prominent molecule, which markedly suppressed cancer stemness and attenuated CR after silencing. Celastrol, a nature extract compound, was demonstrated to effectively inhibit SPC25 expression and reverse CR phenotype. In conclusion, the development of SPC25 inhibitors, such as the application of celastrol, maybe a novel strategy to sensitize cisplatin for the treatment of refractory HNC.

Identification of NUF2 and FAM83D as potential biomarkers in triple-negative breast cancer

Background

Breast cancer is a heterogeneous disease. Compared with other subtypes of breast cancer, triple-negative breast cancer (TNBC) is easy to metastasize and has a short survival time, less choice of treatment options. Here, we aimed to identify the potential biomarkers to TNBC diagnosis and prognosis.

Material/Methods

Three independent data sets (GSE45827, GSE38959, GSE65194) were downloaded from the Gene Expression Omnibus (GEO). The R software packages were used to integrate the gene profiles and identify differentially expressed genes (DEGs). A variety of bioinformatics tools were used to explore the hub genes, including the DAVID database, STRING database and Cytoscape software. Reverse transcription quantitative PCR (RT-qPCR) was used to verify the hub genes in 14 pairs of TNBC paired tissues.

Results

In this study, we screened out 161 DEGs between 222 non-TNBC and 126 TNBC samples, of which 105 genes were up-regulated and 56 were down-regulated. These DEGs were enriched for 27 GO terms and two pathways. GO analysis enriched mainly in “cell division”, “chromosome, centromeric region” and “microtubule motor activity”. KEGG pathway analysis enriched mostly in “Cell cycle” and “Oocyte meiosis”. PPI network was constructed and then 10 top hub genes were screened. According to the analysis results of the Kaplan-Meier survival curve, the expression levels of only NUF2, FAM83D and CENPH were associated with the recurrence-free survival in TNBC samples (P < 0.05). RT-qPCR confirmed that the expression levels of NUF2 and FAM83D in TNBC tissues were indeed up-regulated significantly.

Conclusions

The comprehensive analysis showed that NUF2 and FAM83D could be used as potential biomarkers for diagnosis and prognosis of TNBC.

Long-Term Exposure of Early-Transformed Human Mammary Cells to Low Doses of Benzo[a]pyrene and/or Bisphenol A Enhances Their Cancerous Phenotype via an AhR/GPR30 Interplay

It is of utmost importance to decipher the role of chronic exposure to low doses of environmental carcinogens on breast cancer progression. The early-transformed triple-negative human mammary MCF10AT1 cells were chronically (60 days) exposed to low doses (10⁻¹⁰ M) of Benzo[a]pyrene (B[a]P), a genotoxic agent, and/or Bisphenol A (BPA), an endocrine disruptor. Our study revealed that exposed MCF10AT1 cells developed, in a time-dependent manner, an acquired phenotype characterized by an increase in cancerous properties (anchorage independent growth and stem-like phenotype). Co-exposure of MCF10AT1 cells to B[a]P and BPA led to a significantly greater aggressive phenotype compared to B[a]P or BPA alone. This study provided new insights into the existence of a functional interplay between the aryl hydrocarbon receptor (AhR) and the G protein-coupled receptor 30 (GPR30) by which chronic and low-dose exposure of B[a]P and/or BPA fosters the progression of MCF10AT1 cells into a more aggressive substage. Experiments using AhR or GPR30 antagonists, siRNA strategies, and RNAseq analysis led us to propose a model in which AhR signaling plays a “driver role” in the AhR/GPR30 cross-talk in mediating long-term and low-dose exposure of B[a]P and/or BPA. Retrospective analysis of two independent breast cancer cohorts revealed that the AhR/GPR30 mRNA expression signature resulted in poor breast cancer prognosis, in particular in the ER-negative and the triple-negative subtypes. Finally, the study identified targeting AhR and/or GPR30 with specific antagonists as a strategy capable of inhibiting carcinogenesis associated with chronic exposure to low doses of B[a]P and BPA in MCF10AT1 cells. Altogether, our results indicate that the engagement of both AhR and GPR30 functions, in particular in an ER-negative/triple-negative context of breast cells, favors tumor progression and leads to poor prognosis.

SETD3 acts as a prognostic marker in breast cancer patients and modulates the viability and invasion of breast cancer cells

In several carcinomas, the SET Domain Containing 3, Actin Histidine Methyltransferase (SETD3) is associated with oncogenesis. However, there is little knowledge about the role of SETD3 in the progression and prognosis of breast cancer. In this study, we first analyzed the prognostic value of SETD3 in breast cancer patients using the database of the public Kaplan-Meier plotter. Moreover, in vitro assays were performed to assess the role of SETD3 in the viability and capacity of invasion of human breast cancer cell lines. We observed that the high expression of SETD3 was associated with better relapse-free survival (RFS) of the whole collective of 3,951 patients, of Estrogen Receptor-positive, and of Luminal A-type breast cancer patients. However, in patients lacking expression of estrogen-, progesterone- and HER2-receptor, and those affected by a p53-mutation, SETD3 was associated with poor RFS. In vitro analysis showed that SETD3 siRNA depletion affects the viability of triple-negative cells as well as the cytoskeletal function and capacity of invasion of highly invasive MDA-MB-231 cells. Interestingly, SETD3 regulates the expression of other genes associated with cancer such as β-actin, FOXM1, FBXW7, Fascin, eNOS, and MMP-2. Our study suggests that SETD3 expression can act as a subtype-specific biomarker for breast cancer progression and prognosis.

Prognostic value of DLGAP5 in colorectal cancer

PURPOSE

DLG7 (disc large homolog 7) is a microtubule-associated protein encoded by DLGAP5 (DLG associated protein 5) gene and has an important role during spindle assembly. Spindle assembly deregulation is a well-known cause of genomic instability. The aim of this study was to investigate the influence of DLGAP5 expression on survival and to evaluate its potential use as a biomarker in colorectal cancer (CRC).

METHODS

DLGAP5 expression was measured in the primary tumor and corresponding normal mucosa samples from 109 patients with CRC and correlated to clinical and pathological data. The results were validated in a second, publically available patient cohort. Molecular effects of DLG7/DLGAP5 in CRC were analyzed via functional assays in knockdown cell lines.

RESULTS

DLGAP5 downregulation led to a significant reduction of the invasion and migration potential in CRC. In addition, DLGAP5 expression correlates with nodal status and advanced UICC stage (III-IV).Subgroup analyses revealed a correlation between DLGAP5 overexpression and poor survival in patients with non-metastatic disease (M0). Furthermore, overexpression of DLGAP5 is associated with worse overall survival in distinct molecular CRC subtypes.

CONCLUSIONS

The results of this study suggest the importance of DLGAP5 in defining a more aggressive CRC phenotype. DLG7/DLGAP5 represents a potential biomarker for CRC in molecular subgroups of CRC.

ROCplot.org: Validating predictive biomarkers of chemotherapy/hormonal therapy/anti-HER2 therapy using transcriptomic data of 3,104 breast cancer patients

Systemic therapy of breast cancer can include chemotherapy, hormonal therapy and targeted therapy. Prognostic biomarkers are able to predict survival and predictive biomarkers are able to predict therapy response. In this report, we describe the initial release of the first available online tool able to identify gene expression-based predictive biomarkers using transcriptomic data of a large set of breast cancer patients. Published gene expression data of 36 publicly available datasets were integrated with treatment data into a unified database. Response to therapy was determined using either author-reported pathological complete response data (n = 1,775) or relapse-free survival status at 5 years (n = 1,329). Treatment data includes chemotherapy (n = 2,108), endocrine therapy (n = 971) and anti-human epidermal growth factor receptor 2 (HER2) therapy (n = 267). The transcriptomic database includes 20,089 unique genes and 54,675 probe sets. Gene expression and therapy response are compared using receiver operating characteristics and Mann-Whitney tests. We demonstrate the utility of the pipeline by cross-validating 23 paclitaxel resistance-associated genes in different molecular subtypes of breast cancer. An additional set of established biomarkers including TP53 for chemotherapy in Luminal breast cancer (p = 1.01E-19, AUC = 0.769), HER2 for trastuzumab therapy (p = 8.4E-04, AUC = 0.629) and PGR for hormonal therapy (p = 8.6E-05, AUC = 0.7), are also endorsed. The tool is designed to validate and rank new predictive biomarker candidates in real time. By analyzing the selected genes in a large set of independent patients, one can select the most robust candidates and quickly eliminate those that are most likely to fail in a clinical setting. The analysis tool is accessible at www.rocplot.org.

Examination of the Biomark assay as an alternative to Oncotype DX for defining chemotherapy benefit

Currently the 21-gene recurrence score (RS) assay called Oncotype DX is recommended by the National Comprehensive Cancer Network guideline for defining the benefit of chemotherapy. To overcome the cost disadvantages of the Oncotype DX assay and the turnaround time, a multigene assay was examined to compare the correlation of the RS and the predicted score (PS) of the present study. Paraffin-embedded tissues of 50 cases with early-stage estrogen receptor (ER)-positive breast cancer, who underwent the Oncotype DX test were used. A total of 149 candidate genes with high correlation to the RS were identified, in another project (Lee et al, unpublished data). Reverse transcription-quantitative polymerase chain reaction biomark assays were conducted using the dynamic array integrated fluidic circuit and the correlation analysis was performed with BRB ArrayTools. A predictive model was developed by the coefficient and gene expression, and 41 genes were identified. If the cut-off was ≥18, the predicted model was 18/50 cases, and the RS was 19, indicating that the differential rate of predicted response against RS was 2%. If the cutoff was ≥11, the predicted model was 38/50 cases and the RS was 34, indicating a difference of 8%. Genes common to the Oncotype DX and the Biomark assay include marker of proliferation Ki-67, aurora kinase A, Erb-B2 receptor tyrosine kinase 2, glutathione S-transferase Mu 1, estrogen receptor 1, progesterone receptor, B-cell lymphoma 2, signal peptide CUB domain EGF-like 2 and 5 reference genes. The remaining 28 genes are involved in various pathways and functions. This result indicates that there is a significant correlation between PS and RS scores, although validation of results is required to accurately determine the risk of distant recurrence. The Biomark assay is an easy and inexpensive way to measure mRNA expression. The present study demonstrates the possibility of the Biomark assay as an alternative for defining chemotherapy benefit in individual patients with ER-positive early-stage breast cancer.

Independent validation of induced overexpression efficiency across 242 experiments shows a success rate of 39

Although numerous studies containing induced gene expression have already been published, independent authentication of their results has not yet been performed. Here, we utilized available transcriptomic data to validate the achieved efficiency in overexpression studies. Microarray data of experiments containing cell lines with induced overexpression in one or more genes were analyzed. All together 342 studies were processed, these include 242 different genes overexpressed in 184 cell lines. The final database includes 4,755 treatment-control sample pairs. Successful gene induction (fold change induction over 1.44) was validated in 39.3% of all genes at p < 0.05. Number of repetitions within a study (p < 0.0001) and type of used vector (p = 0.023) had significant impact on successful overexpression efficacy. In summary, over 60% of studies failed to deliver a reproducible overexpression. To achieve higher efficiency, robust and strict study design with multi-level quality control will be necessary.

Exploratory Analysis of Single-Gene Predictive Biomarkers in HERA DASL Cohort Reveals That C8A mRNA Expression Is Prognostic of Outcome and Predictive of Benefit of Trastuzumab

Purpose

The Herceptin Adjuvant study is an international multicenter randomized trial that compared 1 or 2 years of trastuzumab given every 3 weeks with observation in women with human epidermal growth factor 2-positive (HER2+) breast cancer after chemotherapy. Identification of biomarkers predictive of a benefit from trastuzumab will minimize overtreatment and lower health care costs.

Methods

To identify possible single-gene biomarkers, an exploratory analysis of 3,669 gene probes not expected to be expressed in normal breast tissue was conducted. Disease-free survival (DFS) was used as the end point in a Cox regression model, with the interaction term between C8A mRNA and treatment as a categorical variable split on the cohort mean.

Results

A significant interaction between C8A mRNA and treatment was detected (P < .001), indicating a predictive response to trastuzumab treatment. For the C8A-low subgroup (mRNA expression lower than the cohort mean), no significant treatment benefit was observed (P = .73). In the C8A-high subgroup, patients receiving trastuzumab experienced a lower hazard of a DFS event by approximately 75% compared with those in the observation arm (hazard ratio [HR], 0.25; P < .001). A significant prognostic effect of C8A mRNA also was seen (P < .001) in the observation arm, where the C8A-high group hazard of a DFS event was three times the respective hazard of the C8A-low group (HR, 3.27; P < .001). C8A mRNA is highly prognostic in the Hungarian Academy of Science HER2+ gastric cancer cohort (HR, 1.72; P < .001).

Conclusion

C8A as a single-gene biomarker prognostic of DFS and predictive of a benefit from trastuzumab has the potential to improve the standard of care in HER2+ breast cancer if validated by additional studies. Understanding the advantage of overexpression of C8A related to the innate immune response can give insight into the mechanisms that drive cancer.

Distinct prognostic values of mRNA expression of glutathione peroxidases in non-small cell lung cancer

Introduction

Glutathione peroxidases (GPxs) constitutes an enzyme family which has the ability to reduce free hydrogen peroxide to water and lipid hydroperoxides to their corresponding alcohols, and its main biological roles are to protect organisms from oxidative stress-induced damage. GPxs include eight members in different tissues of the body, and they play essential roles in carcinogenesis. However, the prognostic value of individual GPx in non-small cell lung cancer (NSCLC) remains elusive.

Materials and methods

In the current study, we investigated the prognostic value of GPxs in NSCLC patients through the “Kaplan–Meier plotter” database, wherein updated gene expression data and survival information from a total of 1,926 NSCLC patients are included.

Results

High expression of GPx1 mRNA was correlated with worse overall survival (OS) in adenocarcinoma patients. High expression of GPx2 mRNA was correlated with worse OS for all NSCLC patients. In contrast, high expression of GPx3 mRNA was associated with better OS for all NSCLC patients. High expression of GPx4 mRNA was significantly correlated with worsening adenocarcinoma in these patients. GPx5 mRNA high expression correlated with worsening OS for all NSCLC patients.

Discussion

The current findings of prognostic values of individual mRNA expression of GPxs in NSCLC patients indicate some GPxs may have prognostic value in NSCLC patients, and this needs further study.

Clinical and molecular aspects of breast cancer: Targets and therapies

Breast Cancer is a complex disease characterized by the occurrence of multiple molecular alterations. Currently, some molecular markers are in use for breast cancer diagnostic, prognostic, and predictive purposes. Thus, genetic signatures are available for improving the decision-making. The biomarkers are also essential as therapeutic approaches, but many questions remain due to the lack of efficacy on breast cancer treatment, mainly for triple-negative breast cancer subtype. Since the genetic profile of breast cancer can also be related to different ethnic groups and geographic areas, the reference populations of the genetic assays and clinical trials need to include a broader population beyond the European and North American patients. In this review, we analyzed the current and potential molecular markers that could help to improve the strategies for breast cancer therapy.

Identification of endothelial selectin as a potential prognostic marker in breast cancer

Endothelial selectin (ELAM1 or CD62E) has been previously reported as being associated with the prognosis of multiple types of cancer. However, its prognostic value in breast cancer (BC) remains unclear. The aim of the present study was to investigate the prognostic value of ELAM1 mRNA expression in BC tissue. The prognostic value of ELAM1 mRNA was assessed in patients with BC using the Kaplan-Meier plotter (KM-plot) database. The KM-plot generated updated ELAM1 mRNA expression data and survival analysis from a total of 3,951 patients with BC, gathered from 35 datasets. Low expression of ELAM1 mRNA was correlated with a poorer overall survival in 1,402 patients with BC followed for 20 years [hazard ratio (HR), 0.71; 95% confidence interval (CI), 0.57–0.88; log-rank P=0.0016]. Low expression of ELAM1 was also correlated with poorer relapse-free survival (HR, 0.69; 95% CI, 0.62–0.77; log-rank P=2.2e-11) in 3,951 patients and poorer distant metastasis-free survival (HR, 0.79; 95% CI, 0.65–0.96; log-rank P=0.02) in 1,746 patients with BC followed for 20 years. Results from the Metabolic gEne RApid visualizer database indicated that ELAM1 mRNA expression was elevated in normal tissue. The results of the present study suggest that ELAM1 mRNA is a potential prognostic and metastatic marker in patients with BC.

Construction of a novel multi-gene assay (42-gene classifier) for prediction of late recurrence in ER-positive breast cancer patients

PURPOSE

Prediction models for late (> 5 years) recurrence in ER-positive breast cancer need to be developed for the accurate selection of patients for extended hormonal therapy. We attempted to develop such a prediction model focusing on the differences in gene expression between breast cancers with early and late recurrence.

METHODS

For the training set, 779 ER-positive breast cancers treated with tamoxifen alone for 5 years were selected from the databases (GSE6532, GSE12093, GSE17705, and GSE26971). For the validation set, 221 ER-positive breast cancers treated with adjuvant hormonal therapy for 5 years with or without chemotherapy at our hospital were included. Gene expression was assayed by DNA microarray analysis (Affymetrix U133 plus 2.0).

RESULTS

With the 42 genes differentially expressed in early and late recurrence breast cancers in the training set, a prediction model (42GC) for late recurrence was constructed. The patients classified by 42GC into the late recurrence-like group showed a significantly (P = 0.006) higher late recurrence rate as expected but a significantly (P = 1.62 × E-13) lower rate for early recurrence than non-late recurrence-like group. These observations were confirmed for the validation set, i.e., P = 0.020 for late recurrence and P = 5.70 × E-5 for early recurrence.

CONCLUSION

We developed a unique prediction model (42GC) for late recurrence by focusing on the biological differences between breast cancers with early and late recurrence. Interestingly, patients in the late recurrence-like group by 42GC were at low risk for early recurrence.

Immunohistochemical Ki67 after short-term hormone therapy identifies low-risk breast cancers as reliably as genomic markers

Background

The purpose of this study was to test whether immunohistochemical (IHC) Ki67 levels after short-term preoperative hormone therapy (post-Ki67) predict similar numbers of patients with favorable prognoses as genomic markers.

Results

Thirty paired cases (60 samples) were enrolled in this study. Post-Ki67 levels were significantly lower than pre-treatment Ki67 levels (P < 0.001). Post-Ki67 predicted more low-risk cases (83.3%, 25/30) than pre-genomic surrogate signature(GSS) (66.7%: 20/30), but the difference in predictive power was not significant (P = 0.233). Proliferation (MKI67, STK15, Survivin, CCNB1, and MYBL2) and estrogen (ER, PGR, BCL2, and SCUBE2) related signatures were significantly downregulated after therapy (P < 0.001 and 0.041, respectively).

Materials and Methods

Core needle biopsy specimens of primary breast cancer were collected at Okayama University Hospital from hormone receptor-positive and human epidermal growth factor 2-negative patients that subsequently received two weeks of neoadjuvant hormone therapy. Paired post-treatment specimens from surgical samples were also collected. IHC Ki67 levels and GSS were compared between pre- and post-hormone treatment samples. Changes of gene expression pattern in short-term hormone therapy were also assessed.

Conclusions

IHC based post-Ki67 levels may have distinct predictive power compared with the naïve IHC Ki67. Future studies with larger cohorts and longer follow-up periods may be needed to validate our results.

The non-receptor tyrosine kinase TNK2/ACK1 is a novel therapeutic target in triple negative breast cancer

Breast cancer is the most prevalent cancer in women worldwide. About 15-20% of all breast cancers do not express estrogen receptor, progesterone receptor or HER2 receptor and hence are collectively classified as triple negative breast cancer (TNBC). These tumors are often relatively aggressive when compared to other types of breast cancer, and this issue is compounded by the lack of effective targeted therapy. In our previous phosphoproteomic profiling effort, we identified the non-receptor tyrosine kinase TNK2 as activated in a majority of aggressive TNBC cell lines. In the current study, we show that high expression of TNK2 in breast cancer cell lines correlates with high proliferation, invasion and colony forming ability. We demonstrate that knockdown of TNK2 expression can substantially suppress the invasiveness and proliferation advantage of TNBC cells in vitro and tumor formation in xenograft mouse models. Moreover, inhibition of TNK2 with small molecule inhibitor (R)-9bMS significantly compromised TNBC proliferation.

Finally, we find that high levels of TNK2 expression in high-grade basal-like breast cancers correlates significantly with poorer patient outcome. Taken together, our study suggests that TNK2 is a novel potential therapeutic target for the treatment of TNBC.

Cross-validation of survival associated biomarkers in gastric cancer using transcriptomic data of 1,065 patients

Introduction

Multiple gene expression based prognostic biomarkers have been repeatedly identified in gastric carcinoma. However, without confirmation in an independent validation study, their clinical utility is limited. Our goal was to establish a robust database enabling the swift validation of previous and future gastric cancer survival biomarker candidates.

Results

The entire database incorporates 1,065 gastric carcinoma samples, gene expression data. Out of 29 established markers, higher expression of BECN1 (HR = 0.68, p = 1.5E-05), CASP3 (HR = 0.5, p = 6E-14), COX2 (HR = 0.72, p = 0.0013), CTGF (HR = 0.72, p = 0.00051), CTNNB1 (HR = 0.47, p = 4.3E-15), MET (HR = 0.63, p = 1.3E-05), and SIRT1 (HR = 0.64, p = 2.2E-07) correlated to longer OS. Higher expression of BIRC5 (HR = 1.45, p = 1E-04), CNTN1 (HR = 1.44, p = 3.5E- 05), EGFR (HR = 1.86, p = 8.5E-11), ERCC1 (HR = 1.36, p = 0.0012), HER2 (HR = 1.41, p = 0.00011), MMP2 (HR = 1.78, p = 2.6E-09), PFKB4 (HR = 1.56, p = 3.2E-07), SPHK1 (HR = 1.61, p = 3.1E-06), SP1 (HR = 1.45, p = 1.6E-05), TIMP1 (HR = 1.92, p = 2.2E- 10) and VEGF (HR = 1.53, p = 5.7E-06) were predictive for poor OS.

MATERIALS AND METHODS

We integrated samples of three major cancer research centers (Berlin, Bethesda and Melbourne datasets) and publicly available datasets with available follow-up data to form a single integrated database. Subsequently, we performed a literature search for prognostic markers in gastric carcinomas (PubMed, 2012–2015) and re-validated their findings predicting first progression (FP) and overall survival (OS) using uni- and multivariate Cox proportional hazards regression analysis.

Conclusions

The major advantage of our analysis is that we evaluated all genes in the same set of patients thereby making direct comparison of the markers feasible. The best performing genes include BIRC5, CASP3, CTNNB1, TIMP-1, MMP-2, SIRT, and VEGF.

Down-regulation of tensin2 enhances tumorigenicity and is associated with a variety of cancers

Tensin family members, including tensin2 (TNS2), are present as major components of the focal adhesions. The N-terminal end of TNS2 contains a C1 region (protein kinase C conserved region 1) that is not found in other tensin members. Three isoforms of TNS2 have been identified with previous reports describing the shortest V3 isoform as lacking the C1 region. Although TNS2 is known to regulate cell proliferation and migration, its role in tumorigenicity is controversial. By gain-of-function overexpression approaches, results supporting either promotion or reduction of cancer cell tumorigenicity were reported. Here we report that the complete V3 isoform also contains the C1 region and describe the expression patterns of the three human TNS2 isoforms. By loss-of-function approaches, we show that silencing of TNS2 up-regulates the activities of Akt, Mek, and IRS1, and increases tumorigenicities in A549 and Hela cells. Using public database analyses we found that TNS2 is down-regulated in head and neck, esophageal, breast, lung, liver, and colon cancer. In addition, patients with low TNS2 expression showed poor relapse-free survival rates for breast and lung cancers. These results strongly suggest a role of tensin2 in suppressing cell transformation and reduction of tumorigenicity.

VEGFA activates an epigenetic pathway upregulating ovarian cancer-initiating cells

The angiogenic factor, VEGFA, is a therapeutic target in ovarian cancer (OVCA). VEGFA can also stimulate stem‐like cells in certain cancers, but mechanisms thereof are poorly understood. Here, we show that VEGFA mediates stem cell actions in primary human OVCA culture and OVCA lines via VEGFR2‐dependent Src activation to upregulate Bmi1, tumor spheres, and ALDH1 activity. The VEGFA‐mediated increase in spheres was abrogated by Src inhibition or SRC knockdown. VEGFA stimulated sphere formation only in the ALDH1⁺ subpopulation and increased OVCA‐initiating cells and tumor formation in vivo through Bmi1. In contrast to its action in hemopoietic malignancies, DNA methyl transferase 3A (DNMT3A) appears to play a pro‐oncogenic role in ovarian cancer. VEGFA‐driven Src increased DNMT3A leading to miR‐128‐2 methylation and upregulation of Bmi1 to increase stem‐like cells. SRC knockdown was rescued by antagomir to miR‐128. DNMT3A knockdown prevented VEGFA‐driven miR‐128‐2 loss, and the increase in Bmi1 and tumor spheres. Analysis of over 1,300 primary human OVCAs revealed an aggressive subset in which high VEGFA is associated with miR‐128‐2 loss. Thus, VEGFA stimulates OVCA stem‐like cells through Src‐DNMT3A‐driven miR‐128‐2 methylation and Bmi1 upregulation.

linc00673 (ERRLR01) is a prognostic indicator of overall survival in breast cancer

LncRNAs are novel noncoding RNAs involved in the epigenetic regulation of gene expression by recruiting ribonucleoprotein complexes to specific genomic loci to initiate histone methylation and/or other chromatin modifications. LncRNAs themselves function as tumor suppressors or oncogenes, depending on the gene regulatory networks they govern. We identified lnc00673 (ERRLR01) as a marker of overall survival (OS) in breast cancer patients. Specifically, ERRLR01 levels were elevated in triple-negative breast cancer (TNBC) as compared with Luminal-A, Luminal-B, and HER2 breast cancer subtypes. ERRLR01 levels were also inversely correlated with breast cancer survival across all breast cancer patients. Upon stratification, OS in ERα^- tumors correlated with negative overall survival, while in ERα⁺ tumors, ERRLR01 correlated with positive outcomes. This suggests ERRLR01 is modulated by hormone signaling in breast cancer. Gene-network analysis revealed ERRLR01 correlated with distinct pathways including "epithelial development" and "cellular differentiation." These data suggest ERRLR01 operates as an oncogene in TNBC, as well as a biomarker in breast cancer patients.

Association of sirtuins with clinicopathological parameters and overall survival in gastric cancer

To evaluate the associations of sirtuins (SIRT1-7) with clinicopathological parameters in gastric cancer, sirtuins expression profile in NCBI GEO datasets, GSE62254 and GSE15459, was integrated and analyzed. The results suggested that SIRT4, SIRT6, and SIRT7 were associated with Lauren classification and SIRT3-5 were associated with pStage in gastric cancer. Then an online database derived from 1,065 gastric cancer cases, Kaplan-Meier plotter, was used to explore the associations of the mRNA levels of sirtuins with overall survival in gastric cancer. Survival curves generated from Kaplan-Meier plotter suggested that high expression of SIRT1 mRNA was favorable for overall survival in gastric cancer (SIRT1: HR = 0.64, 95% CI = 0.54–0.76, P = 2.2E-07), high expressions of SIRT2-4 and SIRT6-7 were poor for overall survival (SIRT2: HR = 2.31, 95% CI = 1.87–2.87, P = 3.6E-15; SIRT3: HR = 1.99, 95% CI = 1.62–2.45, P = 2.6E-11; SIRT4: HR = 1.41, 95% CI = 1.19–1.68, P = 6.6E-05; SIRT6: HR = 2.02, 95% CI = 1.66–2.47, P = 1.7E-12; SIRT7: HR = 1.96, 95% CI = 1.63–2.35, P = 2.7E-13), whereas no significant association existed between SIRT5 mRNA expression and overall survival. Further analyses stratified by gender, stages, Lauren classification, differentiation, treatment, and HER2 status were also performed. In summary, high SIRT1 mRNA level was associated with better overall survival, SIRT2-4 and 6–7 were associated with poor overall survival, whereas SIRT5 did not show significant association with overall survival in gastric cancer.

DUSP4 is associated with increased resistance against anti-HER2 therapy in breast cancer

The majority of patients develop resistance against suppression of HER2-signaling mediated by trastuzumab in HER2 positive breast cancer (BC). HER2 overexpression activates multiple signaling pathways, including the mitogen-activated protein kinase (MAPK) cascade. MAPK phosphatases (MKPs) are essential regulators of MAPKs and participate in many facets of cellular regulation, including proliferation and apoptosis. We aimed to identify whether differential MKPs are associated with resistance to targeted therapy in patients previously treated with trastuzumab. Using gene chip data of 88 HER2-positive, trastuzumab treated BC patients, candidate MKPs were identified by Receiver Operator Characteristics analysis performed in R. Genes were ranked using their achieved area under the curve (AUC) values and were further restricted to markers significantly associated with worse survival. Functional significance of the two strongest predictive markers was evaluated in vitro by gene silencing in HER2 overexpressing, trastuzumab resistant BC cell lines SKTR and JIMT-1. The strongest predictive MKPs were DUSP4/MKP-2 (AUC=0.75, p=0.0096) and DUSP6/MKP-3 (AUC=0.77, p=5.29E-05). Higher expression for these correlated to worse survival (DUSP4: HR=2.05, p=0.009 and DUSP6: HR=2, p=0.0015). Silencing of DUSP4 had significant sensitization effects – viability of DUSP4 siRNA transfected, trastuzumab treated cells decreased significantly compared to scramble-siRNA transfected controls (SKTR: p=0.016; JIMT-1: p=0.016). In contrast, simultaneous treatment with DUSP6 siRNA and trastuzumab did not alter cell proliferation. Our findings suggest that DUSP4 may represent a new potential target to overcome trastuzumab resistance.

Prognostic roles of mRNA expression of notch receptors in non-small cell lung cancer

Notch signalling is aberrantly activated in human non-small cell lung cancer (NSCLC). Nevertheless, the prognostic roles of mRNA expression of four Notch receptors in NSCLC patients remain elusive. In this report, we reported the prognostic roles of Notch receptors in a total of 1,926 NSCLC patients through “The Kaplan-Meier plotter” (KM plotter) database which is capable to assess the effect of 22,277 genes on survival of NSCLC patients. We found that mRNA high expression level of Notch1 was associated with better overall survival (OS) for all NSCLC patients, hazard ratio (HR) 0.78 (0.69-0.89), p=0.00019, better OS in adenocarcinoma (Ade) patients, HR 0.59 (0.46-0.75), p=1.5e-05, as well as in squamous cell carcinoma (SCC) patients, HR 0.78 (0.62-0.99), p=0.044. mRNA high expression levels of Notch2 and Notch3 were associated with worsen OS for all NSCLC patients, as well as in Ade, but not in SCC patients. mRNA high expression level of Notch4 was not found to be associated with to OS for all NSCLC patients. In addition, mRNA high expression levels of Notch2, Notch3, but Notch4 are significantly associated with the NSCLC patients who have different smoking status. Our results indicate that mRNA expression of Notch receptors may have distinct prognostic values in NSCLC patients. These results will benefit for developing tools to accurately predict the prognosis of NSCLC patients.

Cyclin D1 overexpression correlates with poor tumor differentiation and prognosis in gastric cancer

Overexpression of cyclin D is associated with the molecular tumorigenesis of gastric cancer. The purpose of the present study was to investigate the expression of cyclin D in human gastric cancer and to determine the potential correlations between cyclin D expression and clinicopathological characteristics of specific histological types, as well as its prognostic significance. In the present study, the expression of the cyclin D1 (CCND1), cyclin D2 (CCND2) and cyclin D3 (CCND3) genes in gastric cancer patients was explored using the Oncomine database, and their correlation with overall survival (OS) and progression-free survival (PFS) was evaluated using Kaplan-Meier analysis. The prognostic significance of CCND1 protein expression was evaluated by western blot analysis of 32 matched specimens of gastric adenocarcinomas and normal tissues obtained from patients treated at the National Cheng Kung University Hospital (Tainan, Taiwan). Analysis of the Oncomine cancer microarray database revealed that CCND1 gene expression was significantly increased in gastric intestinal-type adenocarcinoma, while CCND2 was significantly increased in diffuse gastric adenocarcinoma, gastric intestinal-type adenocarcinoma and gastric mixed adenocarcinoma. Kaplan-Meier analysis indicated that overexpression of CCND1 was associated with reduced OS and PFS. In addition, overexpression of CCND1 and downregulation of CCND2 were significantly correlated with receptor tyrosine-protein kinase erb-2-negative tumors and poor differentiation. The ratio of relative CCND1 expression (expressed as the CCND1/β-actin ratio) in tumor tissues compared with that in normal tissues was correlated with poor differentiation (P=0.0018). In summary, CCND1 overexpression is associated with shorter survival in patients with gastric cancer and with poorly differentiated tumors.

VEGFA links self-renewal and metastasis by inducing Sox2 to repress miR-452, driving Slug

Cancer stem cells (CSC) appear to have increased metastatic potential, but mechanisms underlying this are poorly defined. Here we show that VEGFA induction of Sox2 promotes EMT and tumor metastasis. In breast lines and primary cancer culture, VEGFA rapidly upregulates SOX2 expression, leading to SNAI2 induction, EMT, increased invasion and metastasis. We show Sox2 downregulates miR-452, which acts as a novel metastasis suppressor to directly target the SNAI2 3′-untranslated region (3′-UTR). VEGFA stimulates Sox2- and Slug-dependent cell invasion. VEGFA increases lung metastasis in vivo, and this is abrogated by miR-452 overexpression. Furthermore, SNAI2 transduction rescues metastasis suppression by miR-452. Thus, in addition to its angiogenic action, VEGFA upregulates Sox2 to drive stem cell expansion, together with miR-452 loss and Slug upregulation, providing a novel mechanism whereby cancer stem cells acquire metastatic potential. Prior work showed EMT transcription factor overexpression upregulates CSC. Present work indicates that stemness and metastasis are a two-way street: Sox2, a major mediator of CSC self-renewal, also governs the metastatic process.

Benzyl Isothiocyanate potentiates p53 signaling and antitumor effects against breast cancer through activation of p53-LKB1 and p73-LKB1 axes

Functional reactivation of p53 pathway, although arduous, can potentially provide a broad-based strategy for cancer therapy owing to frequent p53 inactivation in human cancer. Using a phosphoprotein-screening array, we found that Benzyl Isothiocynate, (BITC) increases p53 phosphorylation in breast cancer cells and reveal an important role of ERK and PRAS40/MDM2 in BITC-mediated p53 activation. We show that BITC rescues and activates p53-signaling network and inhibits growth of p53-mutant cells. Mechanistically, BITC induces p73 expression in p53-mutant cells, disrupts the interaction of p73 and mutant-p53, thereby releasing p73 from sequestration and allowing it to be transcriptionally active. Furthermore, BITC-induced p53 and p73 axes converge on tumor-suppressor LKB1 which is transcriptionally upregulated by p53 and p73 in p53-wild-type and p53-mutant cells respectively; and in a feed-forward mechanism, LKB1 tethers with p53 and p73 to get recruited to p53-responsive promoters. Analyses of BITC-treated xenografts using LKB1-null cells corroborate in vitro mechanistic findings and establish LKB1 as the key node whereby BITC potentiates as well as rescues p53-pathway in p53-wild-type as well as p53-mutant cells. These data provide first in vitro and in vivo evidence of the integral role of previously unrecognized crosstalk between BITC, p53/LKB1 and p73/LKB1 axes in breast tumor growth-inhibition.

KRAS driven expression signature has prognostic power superior to mutation status in non-small cell lung cancer

KRAS is the most frequently mutated oncogene in non‐small cell lung cancer (NSCLC). However, the prognostic role of KRAS mutation status in NSCLC still remains controversial. We hypothesize that the expression changes of genes affected by KRAS mutation status will have the most prominent effect and could be used as a prognostic signature in lung cancer.

We divided NSCLC patients with mutation and RNA‐seq data into KRAS mutated and wild type groups. Mann‐Whitney test was used to identify genes showing altered expression between these cohorts. Mean expression of the top five genes was designated as a “transcriptomic fingerprint” of the mutation. We evaluated the effect of this signature on clinical outcome in 2,437 NSCLC patients using univariate and multivariate Cox regression analysis.

Mutation of KRAS was most common in adenocarcinoma. Mutation status and KRAS expression were not correlated to prognosis. The transcriptomic fingerprint of KRAS include FOXRED2, KRAS, TOP1, PEX3 and ABL2. The KRAS signature had a high prognostic power. Similar results were achieved when using the second and third set of strongest genes. Moreover, all cutoff values delivered significant prognostic power (p < 0.01). The KRAS signature also remained significant (p < 0.01) in a multivariate analysis including age, gender, smoking history and tumor stage.

We generated a “surrogate signature” of KRAS mutation status in NSCLC patients by computationally linking genotype and gene expression. We show that secondary effects of a mutation can have a higher prognostic relevance than the primary genetic alteration itself.

What's new?

As many as one‐quarter of patients with lung adenocarcinoma (AC), a form of non‐small cell lung cancer (NSCLC), exhibit tumor‐associated mutations in KRAS. Whether KRAS mutation status and expression are correlated to prognosis, however, remains unclear. In this study, a surrogate signature of KRAS mutation status was generated for NSCLC by relating genotype to gene‐expression signature. The approach led to the identification of a significant correlation between overall survival in lung AC and the transcriptomic fingerprint of somatic KRAS mutations. Three genes strongly influenced by KRAS mutation may be relevant to the search for novel NSCLC drug targets.

Colon cancer subtypes: concordance, effect on survival and selection of the most representative preclinical models

Multiple gene-expression-based subtypes have been proposed for the molecular subdivision of colon cancer in the last decade. We aimed to cross-validate these classifiers to explore their concordance and their power to predict survival. A gene-chip-based database comprising 2,166 samples from 12 independent datasets was set up. A total of 22 different molecular subtypes were re-trained including the CCHS, CIN25, CMS, ColoGuideEx, ColoGuidePro, CRCassigner, MDA114, Meta163, ODXcolon, Oncodefender, TCA19, and V7RHS classifiers as well as subtypes established by Budinska, Chang, DeSousa, Marisa, Merlos, Popovici, Schetter, Yuen, and Watanabe (first authors). Correlation with survival was assessed by Cox proportional hazards regression for each classifier using relapse-free survival data. The highest efficacy at predicting survival in stage 2-3 patients was achieved by Yuen (p = 3.9e-05, HR = 2.9), Marisa (p = 2.6e-05, HR = 2.6) and Chang (p = 9e-09, HR = 2.35). Finally, 61 colon cancer cell lines from four independent studies were assigned to the closest molecular subtype.

Validation of RNAi Silencing Efficiency Using Gene Array Data shows 18.5% Failure Rate across 429 Independent Experiments

No independent cross-validation of success rate for studies utilizing small interfering RNA (siRNA) for gene silencing has been completed before. To assess the influence of experimental parameters like cell line, transfection technique, validation method, and type of control, we have to validate these in a large set of studies. We utilized gene chip data published for siRNA experiments to assess success rate and to compare methods used in these experiments. We searched NCBI GEO for samples with whole transcriptome analysis before and after gene silencing and evaluated the efficiency for the target and off-target genes using the array-based expression data. Wilcoxon signed-rank test was used to assess silencing efficacy and Kruskal–Wallis tests and Spearman rank correlation were used to evaluate study parameters. All together 1,643 samples representing 429 experiments published in 207 studies were evaluated. The fold change (FC) of down-regulation of the target gene was above 0.7 in 18.5% and was above 0.5 in 38.7% of experiments. Silencing efficiency was lowest in MCF7 and highest in SW480 cells (FC = 0.59 and FC = 0.30, respectively, P = 9.3E−06). Studies utilizing Western blot for validation performed better than those with quantitative polymerase chain reaction (qPCR) or microarray (FC = 0.43, FC = 0.47, and FC = 0.55, respectively, P = 2.8E−04). There was no correlation between type of control, transfection method, publication year, and silencing efficiency. Although gene silencing is a robust feature successfully cross-validated in the majority of experiments, efficiency remained insufficient in a significant proportion of studies. Selection of cell line model and validation method had the highest influence on silencing proficiency.

A gene expression signature of retinoblastoma loss-of-function is a predictive biomarker of resistance to palbociclib in breast cancer cell lines and is prognostic in patients with ER positive early breast cancer

Palbociclib is a CDK4/6 inhibitor that received FDA approval for treatment of hormone receptor positive (HR+) HER2 negative (HER2neg) advanced breast cancer. To better personalize patients treatment it is critical to identify subgroups that would mostly benefit from it. We hypothesize that complex alterations of the Retinoblastoma (Rb) pathway might be implicated in resistance to CDK4/6 inhibitors and aim to investigate whether signatures of Rb loss-of-function would identify breast cancer cell lines resistant to palbociclib. We established a gene expression signature of Rb loss-of-function (RBsig) by identifying genes correlated with E2F1 and E2F2 expression in breast cancers within The Cancer Genome Atlas. We assessed the RBsig prognostic role in the METABRIC and in a comprehensive breast cancer meta-dataset. Finally, we analyzed whether RBsig would discriminate palbociclib-sensitive and -resistant breast cancer cells in a large RNA sequencing-based dataset. The RBsig was associated with RB1 genetic status in all tumors (p <7e-32) and in luminal or basal subtypes (p < 7e-11 and p < 0.002, respectively). The RBsig was prognostic in the METABRIC dataset (discovery: HR = 1.93 [1.5-2.4] p = 1.4e-08; validation: HR = 2.01 [1.6-2.5] p = 1.3e-09). Untreated and endocrine treated patients with estrogen receptor positive breast cancer expressing high RBsig had significantly worse recurrence free survival compared to those with low RBsig (HR = 2.37 [1.8 - 3.2] p = 1.87e-08 and HR = 2.62 [1.9- 3.5] p = 8.6e-11, respectively). The RBsig was able to identify palbociclib resistant and sensitive breast cancer cells (ROC AUC = 0,7778). Signatures of RB loss might be helpful in personalizing treatment of patients with HR+/HER2neg breast cancer. Further validation in patients receiving palbociclib is warranted.

ExSurv: A Web Resource for Prognostic Analyses of Exons Across Human Cancers Using Clinical Transcriptomes

Survival analysis in biomedical sciences is generally performed by correlating the levels of cellular components with patients' clinical features as a common practice in prognostic biomarker discovery. While the common and primary focus of such analysis in cancer genomics so far has been to identify the potential prognostic genes, alternative splicing - a posttranscriptional regulatory mechanism that affects the functional form of a protein due to inclusion or exclusion of individual exons giving rise to alternative protein products, has increasingly gained attention due to the prevalence of splicing aberrations in cancer transcriptomes. Hence, uncovering the potential prognostic exons can not only help in rationally designing exon-specific therapeutics but also increase specificity toward more personalized treatment options. To address this gap and to provide a platform for rational identification of prognostic exons from cancer transcriptomes, we developed ExSurv (https://exsurv.soic.iupui.edu), a web-based platform for predicting the survival contribution of all annotated exons in the human genome using RNA sequencing-based expression profiles for cancer samples from four cancer types available from The Cancer Genome Atlas. ExSurv enables users to search for a gene of interest and shows survival probabilities for all the exons associated with a gene and found to be significant at the chosen threshold. ExSurv also includes raw expression values across the cancer cohort as well as the survival plots for prognostic exons. Our analysis of the resulting prognostic exons across four cancer types revealed that most of the survival-associated exons are unique to a cancer type with few processes such as cell adhesion, carboxylic, fatty acid metabolism, and regulation of T-cell signaling common across cancer types, possibly suggesting significant differences in the posttranscriptional regulatory pathways contributing to prognosis.

The effects of lymph node status on predicting outcome in ER+ /HER2- tamoxifen treated breast cancer patients using gene signatures

BACKGROUND

Lymph node (LN) status is the most important prognostic variable used to guide ER positive (+) breast cancer treatment. While a positive nodal status is traditionally associated with a poor prognosis, a subset of these patients respond well to treatment and achieve long-term survival. Several gene signatures have been established as a means of predicting outcome of breast cancer patients, but the development and indication for use of these assays varies. Here we compare the capacity of two approved gene signatures and a third novel signature to predict outcome in distinct LN negative (-) and LN+ populations. We also examine biological differences between tumours associated with LN- and LN+ disease.

METHODS

Gene expression data from publically available data sets was used to compare the ability of Oncotype DX and Prosigna to predict Distant Metastasis Free Survival (DMFS) using an in silico platform. A novel gene signature (Ellen) was developed by including patients with both LN- and LN+ disease and using Prediction Analysis of Microarrays (PAM) software. Gene Set Enrichment Analysis (GSEA) was used to determine biological pathways associated with patient outcome in both LN- and LN+ tumors.

RESULTS

The Oncotype DX gene signature, which only used LN- patients during development, significantly predicted outcome in LN- patients, but not LN+ patients. The Prosigna gene signature, which included both LN- and LN+ patients during development, predicted outcome in both LN- and LN+ patient groups. Ellen was also able to predict outcome in both LN- and LN+ patient groups. GSEA suggested that epigenetic modification may be related to poor outcome in LN- disease, whereas immune response may be related to good outcome in LN+ disease.

CONCLUSIONS

We demonstrate the importance of incorporating lymph node status during the development of prognostic gene signatures. Ellen may be a useful tool to predict outcome of patients regardless of lymph node status, or for those with unknown lymph node status. Finally we present candidate biological processes, unique to LN- and LN+ disease, that may indicate risk of relapse.

The prognostic roles of ALDH1 isoenzymes in gastric cancer

Increased aldehyde dehydrogenase 1 (ALDH1) activity has been determined to be present in the stem cells of several kinds of cancers including gastric cancer (GC). Nevertheless, which ones of ALDH1's isoenzymes are leading to ALDH1 activity remains elusive. In this study, we examined the prognostic value and hazard ratio (HR) of individual ALDH1 isoenzymes in patients with GC using "The Kaplan-Meier plotter" database. mRNA high expression level of ALDH1A1 was not found to be significantly correlated with the overall survival (OS) of all patients with GC followed for 20 years, HR =0.86 (95% confidence interval [CI]: 0.7-1.05), P=0.13. mRNA high expression level of ALDH1A2 was also not significantly correlated with OS for all patients with GC, HR =1.13 (95% CI: 0.91-1.41), P=0.25. mRNA high expression level of ALDH1A3 was found to be significantly correlated with worsened OS in either intestinal-type patients, HR =2.24 (95% CI: 1.44-3.49), P=0.00026, or diffuse-type patients, HR =1.91 (95% CI: 1.02-3.59), P=0.04. Interestingly, mRNA high expression level of ALDH1B1 was found to be significantly correlated with better OS for all patients with GC, HR =0.66 (95% CI: 0.53-0.81), P=7.8e-05, and mRNA high expression level of ALDH1L1 was found to be significantly correlated with worsened OS for all patients with GC, HR =1.23 (95% CI: 1-1.51), P=0.048. Furthermore, our results also indicate that ALDH1A3 and ALDH1L1 are potential major contributors to the ALDH1 activity in GC, since mRNA high expression levels of ALDH1A3 and ALDH1L1 were found to be significantly correlated with worsened OS for all patients with GC. Based on our study, ALDH1A3 and ALDH1L1 are potential prognostic markers and therapeutic targets for patients with GC.

Mining distinct aldehyde dehydrogenase 1 (ALDH1) isoenzymes in gastric cancer

Aldehyde dehydrogenase 1 (ALDH1) consists of a family of intracellular enzymes, highly expressed in stem cells populations of leukemia and some solid tumors. Up to now, 6 isoforms of ALDH1 have been reported. However, the expression patterns and the identity of ALDH1 isoenzymes contributing to ALDH1 activity, as well as the prognostic values of ALDH1 isoenzymes in cancers all remain to be elucidated. Here, we studied the expressions of ALDH1 transcripts in gastric cancer (GC) compared with the normal controls using the ONCOMINE database. Through the Kaplan-Meier plotter database, which contains updated gene expression data and survival information of 876 GC patients, we also investigated the prognostic values of ALDH1 isoenzymes in GC patients. It was found that when compared with normal tissues, ALDH1A1 mRNA expression was downregulated, whereas ALDH1A3 and ALDH1B1 were upregulated in GC patients. In survival analyses, high ALDH1A1 and ALDH1B1 expressions were associated with better overall survival (OS) in all GC patients. In addition, high transcription activity of ALDH1A1 predicted better OS in gastric intestinal type adenocarcinoma, but not in diffuse gastric adenocarcinoma. GC patients with high mRNA level of ALDH1B1 showed better OS in gastric intestinal type, and worse OS in diffuse type. Oppositely, high transcription activities of ALDH1A2, ALDH1A3 and ALDH1L1 predicted worsen overall survival in GC patients, suggesting that these isoenzymes might be responsible mainly for the ALDH1 activities in GC. These data provides ALDH1A2, ALDH1A3 and ALDH1L1 as excellent potential targets for individualized treatment of GC patients.