A 23-gene prognostic classifier for prediction of recurrence and survival for Asian breast cancer patients

We report a 23- gene-classifier profiled from Asian women, with the primary purpose of assessing its clinical utility towards improved risk stratification for relapse for breast cancer patients from Asian cohorts within 10 years’ following mastectomy. Four hundred and twenty-two breast cancer patients underwent mastectomy and were used to train the classifier on a logistic regression model. A subset of 197 patients were chosen to be entered into the follow-up studies post mastectomy who were examined to determine the patterns of recurrence and survival analysis based on gene expression of the gene classifier, age at diagnosis, tumor stage and lymph node status, over a 5 and 10 years follow-up period. Metastasis to lymph node (N2-N3) with N0 as the reference (N2 vs. N0 hazard ratio: 2.02 (1.05–8.70), N3 vs. N0 hazard ratio: 4.32 (1.41–13.22) for 5 years) and gene expression of the 23-gene panel (P=0.06, 5 years and 0.02, 10 years, log-rank test) were found to have significant discriminatory effects on the risk of relapse (HR (95%CI):2.50 (0.95–6.50)). Furthermore, survival curves for subgroup analysis with N0-N1 and T1-T2 predicted patients with higher risk scores. The study provides robust evidence of the effectiveness of the 23-gene-classifier and could be used to determine the risk of relapse event (locoregional and distant recurrence) in Asian patients, leading to a meaningful reduction in chemotherapy recommendations.     

A prognostic gene expression profile that predicts circulating tumor cell presence in breast cancer patients

The detection of circulating tumor cells (CTCs) in the peripheral blood and microarray gene expression profiling of the primary tumor are two promising new technologies able to provide valuable prognostic data for patients with breast cancer. Meta-analyses of several established prognostic breast cancer gene expression profiles in large patient cohorts have demonstrated that despite sharing few genes, their delineation of patients into "good prognosis" or "poor prognosis" are frequently very highly correlated, and combining prognostic profiles does not increase prognostic power. In the current study, we aimed to develop a novel profile which provided independent prognostic data by building a signature predictive of CTC status rather than outcome. Microarray gene expression data from an initial training cohort of 72 breast cancer patients for which CTC status had been determined in a previous study using a multimarker QPCR-based assay was used to develop a CTC-predictive profile. The generated profile was validated in two independent datasets of 49 and 123 patients and confirmed to be both predictive of CTC status, and independently prognostic. Importantly, the "CTC profile" also provided prognostic information independent of the well-established and powerful '70-gene' prognostic breast cancer signature. This profile therefore has the potential to not only add prognostic information to currently-available microarray tests but in some circumstances even replace blood-based prognostic CTC tests at time of diagnosis for those patients already undergoing testing by multigene assays.     

Predicting breast cancer using an expression values weighted clinical classifier

Background

Clinical data, such as patient history, laboratory analysis, ultrasound parameters-which are the basis of day-to-day clinical decision support-are often used to guide the clinical management of cancer in the presence of microarray data. Several data fusion techniques are available to integrate genomics or proteomics data, but only a few studies have created a single prediction model using both gene expression and clinical data. These studies often remain inconclusive regarding an obtained improvement in prediction performance. To improve clinical management, these data should be fully exploited. This requires efficient algorithms to integrate these data sets and design a final classifier.LS-SVM classifiers and generalized eigenvalue/singular value decompositions are successfully used in many bioinformatics applications for prediction tasks. While bringing up the benefits of these two techniques, we propose a machine learning approach, a weighted LS-SVM classifier to integrate two data sources: microarray and clinical parameters.

Results

We compared and evaluated the proposed methods on five breast cancer case studies. Compared to LS-SVM classifier on individual data sets, generalized eigenvalue decomposition (GEVD) and kernel GEVD, the proposed weighted LS-SVM classifier offers good prediction performance, in terms of test area under ROC Curve (AUC), on all breast cancer case studies.

Conclusions

Thus a clinical classifier weighted with microarray data set results in significantly improved diagnosis, prognosis and prediction responses to therapy. The proposed model has been shown as a promising mathematical framework in both data fusion and non-linear classification problems.     

Network based consensus gene signatures for biomarker discovery in breast cancer

Diagnostic and prognostic biomarkers for cancer based on gene expression profiles are viewed as a major step towards a better personalized medicine. Many studies using various computational approaches have been published in this direction during the last decade. However, when comparing different gene signatures for related clinical questions often only a small overlap is observed. This can have various reasons, such as technical differences of platforms, differences in biological samples or their treatment in lab, or statistical reasons because of the high dimensionality of the data combined with small sample size, leading to unstable selection of genes. In conclusion retrieved gene signatures are often hard to interpret from a biological point of view. We here demonstrate that it is possible to construct a consensus signature from a set of seemingly different gene signatures by mapping them on a protein interaction network. Common upstream proteins of close gene products, which we identified via our developed algorithm, show a very clear and significant functional interpretation in terms of overrepresented KEGG pathways, disease associated genes and known drug targets. Moreover, we show that such a consensus signature can serve as prior knowledge for predictive biomarker discovery in breast cancer. Evaluation on different datasets shows that signatures derived from the consensus signature reveal a much higher stability than signatures learned from all probesets on a microarray, while at the same time being at least as predictive. Furthermore, they are clearly interpretable in terms of enriched pathways, disease associated genes and known drug targets. In summary we thus believe that network based consensus signatures are not only a way to relate seemingly different gene signatures to each other in a functional manner, but also to establish prior knowledge for highly stable and interpretable predictive biomarkers.     

A prognostic eight-gene expression signature for patients with breast cancer receiving adjuvant chemotherapy

Breast cancer is a popularly diagnosed malignant tumor. Genomic profiling studies suggest that breast cancer is a disease with heterogeneity. Chemotherapy is one of the chief means to treat breast cancer, while its responses and clinical outcomes vary largely due to the conventional clinicopathological factors and inherent chemosensitivity of breast cancer. Using the least absolute shrinkage and selection operator (LASSO) Cox regression model, our study established a multi-mRNA-based signature model and constructed a relative nomogram in predicting distant-recurrence-free survival for patients receiving surgery and following chemotherapy. We constructed a signature of eight mRNAs (IPCEF1, SYNDIG1, TIGIT, SPESP1, C2CD4A, CLCA2, RLN2, and CCL19) with the LASSO model, which was employed to separate subjects into groups with high- and low-risk scores. Obvious differences of distant-recurrence-free survival were found between these two groups. This eight-mRNA-based signature was independently associated with the prognosis and had better prognostic value than classical clinicopathologic factors according to multivariate Cox regression results. Receiver operating characteristic results demonstrated excellent performance in diagnosing 3-year distant-recurrence by the eight-mRNA signature. A nomogram that combined both the eight-mRNA-based signature and clinicopathological risk factors was constructed. Comparing with an ideal model, the nomograms worked well both in the training and validation sets. Through the results that the eight-mRNA signature effectively classified patients into low- and high-risk of distant recurrence, we concluded that this eight-mRNA-based signature played a promising predictive role in prognosis and could be clinically applied in breast cancer patients receiving adjuvant chemotherapy.     

Distinct genes related to drug response identified in ER positive and ER negative breast cancer cell lines

Breast cancer patients have different responses to chemotherapeutic treatments. Genes associated with drug response can provide insight to understand the mechanisms of drug resistance, identify promising therapeutic opportunities, and facilitate personalized treatment. Estrogen receptor (ER) positive and ER negative breast cancer have distinct clinical behavior and molecular properties. However, to date, few studies have rigorously assessed drug response genes in them. In this study, our goal was to systematically identify genes associated with multidrug response in ER positive and ER negative breast cancer cell lines. We tested 27 human breast cell lines for response to seven chemotherapeutic agents (cyclophosphamide, docetaxel, doxorubicin, epirubicin, fluorouracil, gemcitabine, and paclitaxel). We integrated publicly available gene expression profiles of these cell lines with their in vitro drug response patterns, then applied meta-analysis to identify genes related to multidrug response in ER positive and ER negative cells separately. One hundred eighty-eight genes were identified as related to multidrug response in ER positive and 32 genes in ER negative breast cell lines. Of these, only three genes (DBI, TOP2A, and PMVK) were common to both cell types. TOP2A was positively associated with drug response, and DBI was negatively associated with drug response. Interestingly, PMVK was positively associated with drug response in ER positive cells and negatively in ER negative cells. Functional analysis showed that while cell cycle affects drug response in both ER positive and negative cells, most biological processes that are involved in drug response are distinct. A number of signaling pathways that are uniquely enriched in ER positive cells have complex cross talk with ER signaling, while in ER negative cells, enriched pathways are related to metabolic functions. Taken together, our analysis indicates that distinct mechanisms are involved in multidrug response in ER positive and ER negative breast cells.     

Network analysis of gene lists for finding reproducible prognostic breast cancer gene signatures

Many genome-scale studies in molecular biology deliver results in the form of a ranked list of gene names, accordingly to some scoring method. There is always the question how many top-ranked genes to consider for further analysis, for example, in order creating a diagnostic or predictive gene signature for a disease. This question is usually approached from a statistical point of view, without considering any biological properties of top-ranked genes or how they are related to each other functionally. Here we suggest a new method for selecting a number of genes in a ranked gene list such that this set forms the Optimally Functionally Enriched Network (OFTEN), formed by known physical interactions between genes or their products. The method allows associating a network with the gene list, providing easier interpretation of the results and classifying the genes or proteins accordingly to their position in the resulting network. We demonstrate the method on four breast cancer datasets and show that 1) the resulting gene signatures are more reproducible from one dataset to another compared to standard statistical procedures and 2) the overlap of these signatures has significant prognostic potential. The method is implemented in BiNoM Cytoscape plugin (http://binom.curie.fr).     

Joint classifier and feature optimization for comprehensive cancer diagnosis using gene expression data.

Recent research has demonstrated quite convincingly that accurate cancer diagnosis can be achieved by constructing classifiers that are designed to compare the gene expression profile of a tissue of unknown cancer status to a database of stored expression profiles from tissues of known cancer status. This paper introduces the JCFO, a novel algorithm that uses a sparse Bayesian approach to jointly identify both the optimal nonlinear classifier for diagnosis and the optimal set of genes on which to base that diagnosis. We show that the diagnostic classification accuracy of the proposed algorithm is superior to a number of current state-of-the-art methods in a full leave-one-out cross-validation study of five widely used benchmark datasets. In addition to its superior classification accuracy, the algorithm is designed to automatically identify a small subset of genes (typically around twenty in our experiments) that are capable of providing complete discriminatory information for diagnosis. Focusing attention on a small subset of genes is useful not only because it produces a classifier with good generalization capacity, but also because this set of genes may provide insights into the mechanisms responsible for the disease itself. A number of the genes identified by the JCFO in our experiments are already in use as clinical markers for cancer diagnosis; some of the remaining genes may be excellent candidates for further clinical investigation. If it is possible to identify a small set of genes that is indeed capable of providing complete discrimination, inexpensive diagnostic assays might be widely deployable in clinical settings.     

A novel autophagy-related lncRNA prognostic risk model for breast cancer

Long non-coding RNAs (lncRNAs) are well known as crucial regulators to breast cancer development and are implicated in controlling autophagy. LncRNAs are also emerging as valuable prognostic factors for breast cancer patients. It is critical to identify autophagy-related lncRNAs with prognostic value in breast cancer. In this study, we identified autophagy-related lncRNAs in breast cancer by constructing a co-expression network of autophagy-related mRNAs-lncRNAs from The Cancer Genome Atlas (TCGA). We evaluated the prognostic value of these autophagy-related lncRNAs by univariate and multivariate Cox proportional hazards analyses and eventually obtained a prognostic risk model consisting of 11 autophagy-related lncRNAs (U62317.4, LINC01016, LINC02166, C6orf99, LINC00992, BAIAP2-DT, AC245297.3, AC090912.1, Z68871.1, LINC00578 and LINC01871). The risk model was further validated as a novel independent prognostic factor for breast cancer patients based on the calculated risk score by Kaplan-Meier analysis, univariate and multivariate Cox regression analyses and time-dependent receiver operating characteristic (ROC) curve analysis. Moreover, based on the risk model, the low-risk and high-risk groups displayed different autophagy and oncogenic statues by principal component analysis (PCA) and Gene Set Enrichment Analysis (GSEA) functional annotation. Taken together, these findings suggested that the risk model of the 11 autophagy-related lncRNAs has significant prognostic value for breast cancer and might be autophagy-related therapeutic targets in clinical practice.     

NEFL mRNA expression level is a prognostic factor for early-stage breast cancer patients

Neurofilament, light polypeptide (NEFL) was demonstrated to be ectopically expressed in breast cancer tissues and decreased in lymph node metastases compared to the paired primary breast cancers in our previous study. Moreover, in several studies, NEFL was regarded as a tumor suppressor gene, and its loss of heterozygosity (LOH) was related to carcinogenesis and metastasis in several types of cancer. To explore the role of NEFL in the progression of breast cancer and to evaluate its clinical significance, we detected the NEFL mRNA level in normal breast tissues, primary breast cancer samples and lymph node metastases, and then analyzed the association between the NEFL expression level and several clinicopathological parameters and disease-free survival (DFS). NEFL mRNA was found to be expressed in 92.3% of breast malignancies and down-regulated in lymph node metastases compared to the paired primary tumors. NEFL mRNA level was lower in primary breast cancers with positive lymph nodes than in cancers with negative lymph nodes. Moreover, a low expression level of NEFL mRNA indicated a poor five-year DFS for early-stage breast cancer patients. Thus, NEFL mRNA is ectopically expressed in breast malignancies and could be a potential prognostic factor for early-stage breast cancer patients.     

A deep learning classifier for digital breast tomosynthesis

PurposeTo develop a computerized detection system for the automatic classification of the presence/absence of mass lesions in digital breast tomosynthesis (DBT) annotated exams, based on a deep convolutional neural network (DCNN).Materials and MethodsThree DCNN architectures working at image-level (DBT slice) were compared: two state-of-the-art pre-trained DCNN architectures (AlexNet and VGG19) customized through transfer learning, and one developed from scratch (DBT-DCNN). To evaluate these DCNN-based architectures we analysed their classification performance on two different datasets provided by two hospital radiology departments. DBT slice images were processed following normalization, background correction and data augmentation procedures. The accuracy, sensitivity, and area-under-the-curve (AUC) values were evaluated on both datasets, using receiver operating characteristic curves. A Grad-CAM technique was also implemented providing an indication of the lesion position in the DBT slice.Results Accuracy, sensitivity and AUC for the investigated DCNN are in-line with the best performance reported in the field. The DBT-DCNN network developed in this work showed an accuracy and a sensitivity of (90% ± 4%) and (96% ± 3%), respectively, with an AUC as good as 0.89 ± 0.04. A k-fold cross validation test (with k = 4) showed an accuracy of 94.0% ± 0.2%, and a F1-score test provided a value as good as 0.93 ± 0.03. Grad-CAM maps show high activation in correspondence of pixels within the tumour regions.Conclusions We developed a deep learning-based framework (DBT-DCNN) to classify DBT images from clinical exams. We investigated also a possible application of the Grad-CAM technique to identify the lesion position.     

Comprehensive analysis of gene expression profiles to identify differential prognostic factors of primary and metastatic breast cancer

Breast cancer accounts for nearly half of all cancer-related deaths in women worldwide. However, the molecular mechanisms that lead to tumour development and progression remain poorly understood and there is a need to identify candidate genes associated with primary and metastatic breast cancer progression and prognosis. In this study, candidate genes associated with prognosis of primary and metastatic breast cancer were explored through a novel bioinformatics approach. Primary and metastatic breast cancer tissues and adjacent normal breast tissues were evaluated to identify biomarkers characteristic of primary and metastatic breast cancer. The Cancer Genome Atlas-breast invasive carcinoma (TCGA-BRCA) dataset (ID: HS-01619) was downloaded using the mRNASeq platform. Genevestigator 8.3.2 was used to analyse TCGA-BRCA gene expression profiles between the sample groups and identify the differentially-expressed genes (DEGs) in each group. For each group, Gene Ontology and Kyoto Encyclopedia of Genes and Genomes pathway enrichment analyses were used to determine the function of DEGs. Networks of protein–protein interactions were constructed to identify the top hub genes with the highest degree of interaction. Additionally, the top hub genes were validated based on overall survival and immunohistochemistry using The Human Protein Atlas. Of the top 20 hub genes identified, four (KRT14, KIT, RAD51, and TTK) were considered as prognostic risk factors based on overall survival. KRT14 and KIT expression levels were upregulated while those of RAD51 and TTK were downregulated in patients with breast cancer. The four proposed candidate hub genes might aid in further understanding the molecular changes that distinguish primary breast tumours from metastatic tumours as well as help in developing novel therapeutics. Furthermore, they may serve as effective prognostic risk markers based on the strong correlation between their expression and patient overall survival.     

Nuclear NHERF1 expression as a prognostic marker in breast cancer

Our purpose was to investigate whether Na⁺/H⁺ exchanger regulatory factor 1 (NHERF1) expression could be linked to prognosis in invasive breast carcinomas. NHERF1, an ezrin-radixin-moesin (ERM) binding phosphoprotein 50, is involved in the linkage of integral membrane proteins to the cytoskeleton. It is therefore believed to have an important role in cell signaling associated with changes in cell cytoarchitecture. NHERF1 expression is observed in various types of cancer and is related to tumor aggressiveness. To date the most extensive analyses of the influence of NHERF1 in cancer development have been performed on breast cancer. However, the underlying mechanism and its prognostic significance are still undefined. NHERF1 expression was studied by immunohistochemistry (IHC) in a cohort of 222 breast carcinoma patients. Association of cytoplasmic and nuclear NHERF1 expression with survival was analyzed. Disease-free survival (DFS) and overall survival (OS) were determined based on the Kaplan–Meier method. Cytoplasmic NHERF1 expression was associated with negative progesterone receptor (PgR) (P=0.017) and positive HER2 expression (P=0.023). NHERF1 also showed a nuclear localization and this correlated with small tumor size (P=0.026) and positive estrogen receptor (ER) expression (P=0.010). Multivariate analysis identified large tumor size (P=0.011) and nuclear NHERF1 expression (P=0.049) to be independent prognostic variables for DFS. Moreover, the nuclear NHERF1(−)/ER(−) immunophenotype (27%) was statistically associated with large tumor size (P=0.0276), high histological grade (P=0.0411), PgR-negative tumors (P<0.0001) and high proliferative activity (P=0.0027). These patients had worse DFS compared with patients with nuclear NHERF1(+)/ER(+) tumors (75.4% versus 92.6% P=0.010). These results show that the loss of nuclear NHERF1 expression is associated with reduced survival, and the link between nuclear NHERF1 and ER expression may serve as a prognostic marker for the routine clinical management of breast cancer patients.     

A prognostic eight-lncRNA expression signature in predicting recurrence of ER-positive breast cancer receiving endocrine therapy

Long noncoding RNAs (lncRNAs) have the main role in the tumorigenesis of breast cancer. In the present study, lncRNA expression profiling was collected to identify a lncRNA expression signature from the Gene Expression Omnibus database. An eight-lncRNA signature was established to predict the survival of patients with estrogen receptor (ER)-positive breast cancer receiving endocrine therapy. Patients were separated into a low-risk group and a high-risk group based on this signature. Patients in high-risk group have worse survival compared to those in low-risk group using Kaplan–Meier curve analysis with log-rank test. Receiver operating characteristic analysis suggested good diagnostic efficiency of the eight-lncRNA signature. When adjusting the clinical features, including age, grade, lymph node status, and tumor size, this signature was independently associated with the relapse-free survival. The prognostic value of the lncRNA prognostic model was then validated in validation sets. When validated in a cohort of patients treated with neoadjuvant chemotherapy and endocrine therapy, this signature demonstrated good performance as well. Besides, we have built a nomogram that integrated the conventional clinicopathological features and the eight-lncRNA-based signature. To sum up, our results indicated that the eight-lncRNA prognostic model was a reliable tool to group patients at high and low risk of disease relapse. This signature may have possible implication in prognostic evaluations of patients with ER-positive breast cancer receiving endocrine therapy.     

A prognostic 11 long noncoding RNA expression signature for breast invasive carcinoma

Breast cancer, the most common cancer in women worldwide, is associated with high mortality. The long non-coding RNAs (lncRNAs) with a little capacity of coding proteins is playing an increasingly important role in the cancer paradigm. Accumulating evidences demonstrate that lncRNAs have crucial connections with breast cancer prognosis while the studies of lncRNAs in breast cancer are still in its primary stage. In this study, we collected 1052 clinical patient samples, a comparatively large sample size, including 13 159 lncRNA expression profiles of breast invasive carcinoma (BRCA) from The Cancer Genome Atlas database to identify prognosis-related lncRNAs. We randomly separated all of these clinical patient samples into training and testing sets. In the training set, we performed univariable Cox regression analysis for primary screening and played the model for Robust likelihood-based survival for 1000 times. Then 11 lncRNAs with a frequency more than 600 were selected for prediction of the prognosis of BRCA. Using the analysis of multivariate Cox regression, we established a signature risk-score formula for 11 lncRNA to identify the relationship between lncRNA signatures and overall survival. The 11 lncRNA signature was validated both in the testing and the complete set and could effectively classify the high-/low-risk group with different OS. We also verified our results in different stages. Moreover, we analyzed the connection between the 11 lncRNAs and the genes of ESR1, PGR, and Her2, of which protein products (ESR, PGR, and HER2) were used to classify the breast cancer subtypes widely. The results indicated correlations between 11 lncRNAs and the gene of PGR and ESR1. Thus, a prognostic model for 11 lncRNA expression was developed to classify the BRAC clinical patient samples, providing new avenues in understanding the potential therapeutic methods of breast cancer.     

Regulatory mechanisms for abnormal expression of the human breast cancer specific gene 1 in breast cancer cells

Breast cancer specific gene 1, also referred as synu-clein γ, was originally isolated from a human breasttumor cDNA library[1]. It reveals extensive sequencehomology to a family of neuronal cytosolic proteins,synuclein α and synuclein β[2,3]. Synuclein…     

Isolation and characterization of a metastatic hybrid cell line generated by ER negative and ER positive breast cancer cells in mouse bone marrow

Background

The origin and the contribution of breast tumor heterogeneity to its progression are not clear. We investigated the effect of a growing orthotopic tumor formed by an aggressive estrogen receptor (ER)-negative breast cancer cell line on the metastatic potential of a less aggressive ER-positive breast cancer cell line for the elucidation of how the presence of heterogeneous cancer cells might affect each other''s metastatic behavior.

Methods

ER positive ZR-75-1/GFP/puro cells, resistant to puromycin and non-tumorigenic/non-metastatic without exogenous estrogen supplementation, were injected intracardiacally into mice bearing growing orthotopic tumors, formed by ER negative MDA-MB-231/GFP/Neo cells resistant to G418. A variant cell line B6, containing both estrogen-dependent and -independent cells, were isolated from GFP expressing cells in the bone marrow and re-inoculated in nude mice to generate an estrogen-independent cell line B6TC.

Results

The presence of ER negative orthotopic tumors resulted in bone metastasis of ZR-75-1 without estrogen supplementation. The newly established B6TC cell line was tumorigenic without estrogen supplementation and resistant to both puromycin and G418 suggesting its origin from the fusion of MDA-MB-231/GFP/Neo and ZR-75-1/GFP/puro in the mouse bone marrow. Compared to parental cells, B6TC cells were more metastatic to lung and bone after intracardiac inoculation. More significantly, B6TC mice also developed brain metastasis, which was not observed in the MDA-MB-231/GFP/Neo cell-inoculated mice. Low expression of ERα and CD24, and high expression of EMT-related markers such as Vimentin, CXCR4, and Integrin-β1 along with high CD44 and ALDH expression indicated stem cell-like characteristics of B6TC. Gene microarray analysis demonstrated a significantly different gene expression profile of B6TC in comparison to those of parental cell lines.

Conclusions

Spontaneous generation of the novel hybrid cell line B6TC, in a metastatic site with stem cell-like properties and propensity to metastasize to brain, suggest that cell fusion can contribute to tumor heterogeneity.     

Prognostic classification of breast cancer and gene expression profiling

Clinical and pathological heterogeneity of breast cancer, partly responsible of therapeutic failures, reflects complex and combinatory molecular alterations until now poorly documented by classical investigation tools. Thorough molecular typing is crucial. The advent of DNA microarray-based gene expression profiling allowed consistent progresses in this direction. A novel molecular taxonomy of breast cancer has been defined, signatures that predict clinical outcome or therapeutic response have been identified, some of them being tested in ongoing prospective clinical trials. In this review, we present the main results and their potential clinical applications. We also discuss their current limits and future hopes in the therapeutic management of patients.