Smid, M.

KLF6-SV1 Drives Breast Cancer Metastasis and Is Associated with Poor Survival (Article)

2013-01-23T00:00:00Z

Metastasis is the major cause of cancer mortality. A more thorough understanding of the mechanisms driving this complex multistep process will aid in the identification and characterization of therapeutically targetable genetic drivers of disease progression. We demonstrate that KLF6-SV1, an oncogenic splice variant of the KLF6 tumor suppressor gene, is associated with increased metastatic potential and poor survival in a cohort of 671 lymph node-negative breast cancer patients. KLF6-SV1 overexpression in mammary epithelial cell lines resulted in an epithelial-to-mesenchymal-like transition and drove aggressive multiorgan metastatic disease in multiple in vivo models. Additionally, KLF6-SV1 loss-of-function studies demonstrated reversion to an epithelial and less invasive phenotype. Combined, these findings implicate KLF6-SV1 as a key driver of breast cancer metastasis that distinguishes between indolent and lethal early-stage disease and provides a potential therapeutic target for invasive breast cancer.

High TWIST1 mRNA expression is associated with poor prognosis in lymph node-negative and estrogen receptor-positive human breast cancer and is co-expressed with stromal as well as ECM related genes (Article)

2012-09-11T00:00:00Z

Introduction: The TWIST homolog 1 (TWIST1) is a transcription factor that induces epithelial to mesenchymal transition (EMT), a key process in metastasis. The purpose of this study was to investigate whether TWIST1 expression predicts disease progression in a large breast cancer cohort with long-term clinical follow-up, and to reveal the biology related to TWIST1 mediated disease progression.Methods: TWIST1 mRNA expression level was analyzed by quantitative real-time reverse polymerase chain reaction (RT-PCR) in 1,427 primary breast cancers. In uni- and multivariate analysis using Cox regression, TWIST1 mRNA expression level was associated with metastasis-free survival (MFS), disease-free survival (DFS) and overall survival (OS). Separate analyses in lymph node-negative patients (LNN, n = 778) who did not receive adjuvant systemic therapy, before and after stratification into estrogen receptor (ER)-positive (n = 552) and ER-negative (n = 226) disease, were also performed. The association of TWIST1 mRNA with survival endpoints was assessed using Kaplan-Meier analysis. Using gene expression arrays, genes showing a significant Spearman rank correlation with TWIST1 were used to identify overrepresented Gene Ontology (GO) terms and Kyoto Encyclopedia of Genes and Genomes (KEGG)-annotated biological pathways.Results: Increased mRNA expression level of TWIST1 analyzed as a continuous variable in both uni- and multivariate analysis was associated with shorter MFS in all patients (hazard ratio (HR): 1.17, 95% confidence interval, (95% CI):1.09 to 1.26; and HR: 1.17, 95% CI: 1.08 to 1.26; respectively), in LNN patients (HR: 1.22, 95% CI: 1.09 to 1.36; and HR: 1.21, 95% CI: 1.07 to 1.36; respectively) and in the ER-positive subgroup of LNN patients (HR: 1.34, 95% CI: 1.17 to 1.53; and HR: 1.32, 95% CI: 1.14 to 1.53; respectively). Similarly, high TWIST1 expression was associated with shorter DFS and OS in all patients and in the LNN/ER-positive subgroup. In contrast, no association of TWIST1 mRNA expression with MFS, DFS or OS was observed in ER-negative patients. Genes highly correlated with TWIST1 were significantly enriched for cell adhesion and ECM-related signaling pathways. Furthermore, TWIST1 mRNA was highly expressed in tumor stroma and positively related to tumor stromal content (P <0.001).Conclusions: TWIST1 mRNA expression is an independent prognostic factor for poor prognosis in LNN/ER-positive breast cancer. The biological associations suggest an involvement of the tumor microenvironment in TWIST1's adverse role in breast cancer.

mRNA and microRNA expression profiles in circulating tumor cells and primary tumors of metastatic breast cancer patients (Article)

2011-06-01T00:00:00Z

Purpose: Molecular characterization of circulating tumor cells (CTC) holds great promise. Unfortunately, routinely isolated CTC fractions currently still contain contaminating leukocytes, which makes CTC-specific molecular characterization extremely challenging. In this study, we determined mRNA and microRNA (miRNA) expression of potentially CTC-specific genes that are considered to be clinically relevant in breast cancer. Experimental Design: CTCs were isolated with the epithelial cell adhesion molecule-based CellSearch Profile Kit. Selected genes were measured by real-time reverse transcriptase PCR in CTCs of 50 metastatic breast cancer patients collected before starting first-line systemic therapy in blood from 53 healthy blood donors (HBD) and in primary tumors of 8 of the patients. The molecular profiles were associated with CTC counts and clinical parameters and compared with the profiles generated from the corresponding primary tumors. Results: We identified 55 mRNAs and 10 miRNAs more abundantly expressed in samples from 32 patients with at least 5 CTCs in 7.5 mL of blood compared with samples from 9 patients without detectable CTCs and HBDs. Clustering analysis resulted in 4 different patient clusters characterized by 5 distinct gene clusters. Twice the number of patients from cluster 2 to 4 had developed both visceral and nonvisceral metastases. Comparing transcript levels in CTCs with those measured in corresponding primary tumors showed clinically relevant discrepancies in estrogen receptor and HER2 levels. Conclusions: Our study shows that molecular profiling of low numbers of CTCs in a high background of leukocytes is feasible and shows promise for further studies on the clinical relevance of molecular characterization of CTCs.

Gene expression profiling assigns CHEK2 1100delC breast cancers to the luminal intrinsic subtypes (Article)

2011-05-26T00:00:00Z

CHEK2 1100delC is a moderate-risk cancer susceptibility allele that confers a high breast cancer risk in a polygenic setting. Gene expression profiling of CHEK2 1100delC breast cancers may reveal clues to the nature of the polygenic CHEK2 model and its genes involved. Here, we report global gene expression profiles of a cohort of 155 familial breast cancers, including 26 CHEK2 1100delC mutant tumors. In line with previous work, all CHEK2 1100delC mutant tumors clustered among the hormone receptor-positive breast cancers. In the hormone receptor-positive subset, a 40-gene CHEK2 signature was subsequently defined that significantly associated with CHEK2 1100delC breast cancers. The identification of a CHEK2 gene signature implies an unexpected biological homogeneity among the CHEK2 1100delC breast cancers. In addition, all 26 CHEK2 1100delC tumors classified as luminal intrinsic subtype breast cancers, with 8 luminal A and 18 luminal B tumors. This biological make-up of CHEK2 1100delC breast cancers suggests that a relatively limited number of additional susceptibility alleles are involved in the polygenic CHEK2 model. Identification of these as-yet-unknown susceptibility alleles should be aided by clues from the 40-gene CHEK2 signature.

Pathway analysis of gene lists associated with platinum-based chemotherapy resistance in ovarian cancer: the big picture (Article)

2010-05-01T00:00:00Z

Abstract OBJECTIVE: Ovarian cancer is the leading cause of death from gynecological cancers in the Western world (Parkin et al., 2005). The overall 5-year survival is only 30% (Moss and Kaye, 2002), which is for a significant part due to platinum-based chemotherapy resistance. In this study, we performed a pathway analysis on nine published gene sets associated with platinum resistance in ovarian cancer, including a study by us. With this exploratory study, we aim to identify overlapping pathways associated with platinum-based chemotherapy resistance mechanisms in ovarian cancer. METHODS: Gene Ontology (GO) analysis and Ingenuity Pathway Analysis (IPA) were performed to determine which functional processes were differentially represented in the combined gene lists of nine studies (457 genes) compared to all Unigene identifiers or the Ingenuity knowledge base. RESULTS: The GO and IPA analysis resulted in the generation of 23 gene networks, and showed that 13 GO processes (>or=2 times enriched), 71 canonical pathways (p<0.05,), eight toxicity pathways (p<0.05) and 74 biological functions (p<0.005) are significantly associated with the 9-study gene set. CONCLUSION AND RECOMMENDATIONS: Several pathways identified have previously been shown to be associated with therapy resistance: these include 'oxidative stress response mediated by Nrf2,' 'TP53 signaling' and 'TGFbeta signaling.' The role of TGFbeta signaling and related miRNAs identified in the network analysis in epithelial-to-mesenchymal transition (EMT) and stemness as well as the possible relation with platin-based chemotherapy resistance are further discussed in detail. We propose that future international cooperation should aim at a uniform pooled analysis of the wealth of ovarian cancer array data already available. This will enhance the power of each separate ovarian cancer study and can lead to promising results.

Distinct gene mutation profiles among luminal-type and basal-type breast cancer cell lines (Article)

2010-05-01T00:00:00Z

Breast cancer has for long been recognized as a highly diverse tumor group, but the underlying genetic basis has been elusive. Here, we report an extensive molecular characterization of a collection of 41 human breast cancer cell lines. Protein and gene expression analyses indicated that the collection of breast cancer cell lines has retained most, if not all, molecular characteristics that are typical for clinical breast cancers. Gene mutation analyses identified 146 oncogenic mutations among 27 well-known cancer genes, amounting to an average of 3.6 mutations per cell line. Mutations in genes from the p53, RB and PI3K tumor suppressor pathways were widespread among all breast cancer cell lines. Most important, we have identified two gene mutation profiles that are specifically associated with luminal-type and basal-type breast cancer cell lines. The luminal mutation profile involved E-cadherin and MAP2K4 gene mutations and amplifications of Cyclin D1, ERBB2 and HDM2, whereas the basal mutation profile involved BRCA1, RB1, RAS and BRAF gene mutations and deletions of p16 and p14ARF. These subtype-specific gene mutation profiles constitute a genetic basis for the heterogeneity observed among human breast cancers, providing clues for their underlying biology and providing guidance for targeted pharmacogenetic intervention in breast cancer patients.

Patterns and incidence of chromosomal instability and their prognostic relevance in breast cancer subtypes (Article)

2010-01-01T00:00:00Z

One of the hallmarks of human solid tumors is chromosomal instability (CIN). We studied global patterns as well as individual levels of CIN and determined the prognostic relevance among breast cancer subtypes. For this, we used single nucleotide polymorphism copy number data of 313 primary lymph-node negative breast cancers. The level of CIN for individual samples was determined by counting the total number of chromosomal segments showing a gain or loss per specimen. Hierarchical clustering resulted in four groups showing distinct patterns of abnormalities, predominantly characterized by 1q gain, 8q gain, 1q&8q gain, or no gain of these loci. Estrogen receptor (ER)-positive and ER-negative samples showed an uneven distribution (statistically significant) across the cluster-groups, as did the molecular subtypes and triple-negative tumors (negative for estrogen-, progesterone-, and her2/neu-receptor). The CIN-score was significantly higher in ER-negative and triple-negative samples. Among luminal cancers, luminal B had a higher CIN-score than luminal A. The CIN-score was significantly associated with prognosis, measured by the time to distant metastasis, in ER-positive, luminal B, and her2/neu subtypes, but not in ER-negative patients. Our study points to a multifaceted role for CIN in breast cancer. Within ER-negative samples, CIN is likely related to the onset but other factors govern the progression of the disease. In contrast, CIN is clearly associated with progression in ER-positive, luminal B, and her2/neu subtypes; thus, assessing CIN in these subtypes may contribute to personalized patient management.

CITED2 and NCOR2 in anti-oestrogen resistance and progression of breast cancer (Article)

2009-12-01T00:00:00Z

Background:Endocrine therapies of breast cancer are effective but ultimately fail because of the development of treatment resistance. We have previously revealed several genes leading to tamoxifen resistance in vitro by retroviral insertion mutagenesis. To understand the manner in which these genes yield tamoxifen resistance, their effects on global gene expression were studied and those genes resulting in a distinct gene expression profile were further investigated for their clinical relevance.Methods:Gene expression profiles of 69 human breast cancer cell lines that were made tamoxifen resistant through retroviral insertion mutagenesis were obtained using oligonucleotide arrays and analysed with bioinformatic tools. mRNA levels of NCOR2 and CITED2 in oestrogen receptor-positive breast tumours were determined by quantitative RT-PCR. mRNA levels were evaluated for association with metastasis-free survival (MFS) in 620 patients with lymph node-negative primary breast cancer who did not receive systemic adjuvant therapy, and with clinical benefit in 296 patients receiving tamoxifen therapy for recurrent breast cancer.Results:mRNA expression profiles of most tamoxifen-resistant cell lines were strikingly similar, except for the subgroups of cell lines in which NCOR2 or CITED2 were targeted by the retrovirus. Both NCOR2 and CITED2 mRNA levels were associated with MFS, that is, tumour aggressiveness, independently of traditional prognostic factors. In addition, high CITED2 mRNA levels were predictive for a clinical benefit from first-line tamoxifen treatment in patients with advanced disease.Conclusions: Most retrovirally targeted genes yielding tamoxifen resistance in our cell lines do not impose a distinctive expression profile, suggesting that their causative role in cell growth may be accomplished by post-transcriptional processes. The associations of NCOR2 and CITED2 with outcome in oestrogen receptor-positive breast cancer patients underscore the clinical relevance of functional genetic screens to better understand disease progression, which may ultimately lead to the development of improved treatment options.

Functional identification of genes causing estrogen independence of human breast cancer cells (Article)

2009-03-01T00:00:00Z

Endocrine treatment of breast cancer is widely applied and effective. However, in advanced disease cases, the tumors will eventually progress into an estrogen-independent and therapy-resistant phenotype. To elucidate the molecular mechanisms underlying this endocrine therapy failure, we applied retroviral insertion mutagenesis to identify the main genes conferring estrogen independence to human breast cancer cells. Estrogen-dependent ZR-75-1 cells were infected with replication-defective retroviruses followed by selection with the anti-estrogen 4-hydroxy-tamoxifen. In the resulting panel of 79 tamoxifen-resistant cell lines, the viral integrations were mapped within the human genome. Genes located in the immediate proximity of the retroviral integration sites were characterized for altered expression and their capacity to confer anti-estrogen resistance when transfected into breast cancer cells. Out of 15 candidate BCAR (breast cancer anti-estrogen resistance) genes, seven (AKT1, AKT2, BCAR1, BCAR3, EGFR, GRB7, and TRERF1/BCAR2) were shown to directly underlie estrogen independence. Our results show that insertion mutagenesis is a powerful tool to identify BCAR loci, which may provide insights into the molecular and cellular mechanisms of breast tumor progression and therapy resistance thereby offering novel targets for the development of tailor-made therapeutical and prevention strategies.

Urokinase receptor splice variant uPAR-del4/5-associated gene expression in breast cancer: Identification of rab31 as an independent prognostic factor (Article)

2008-09-01T00:00:00Z

Purpose: To evaluate the pure prognostic impact of the uPA-receptor splice variant uPAR-del4/5 for lymph node-negative breast cancer patients, and to identify differentially expressed genes associated with high or low uPAR-del4/5 mRNA levels. Patients and methods: mRNA transcript levels were measured by real-time PCR in tumor samples from 280 node-negative breast cancer patients who had not received adjuvant systemic therapy. Endpoints were distant metastasis-free survival (DMFS) and overall survival (OS). Gene expression analysis was performed with RNA isolated from breast cancer tissue and breast cancer cell lines using Affymetrix U133a GeneChips. Results: In multivariate analysis, uPAR-del4/5 significantly contributed to the base model of traditional prognostic factors for DMFS (HR = 3.29, P < 0.001) and OS (HR = 2.87, P = 0.002). Using microarrays, seven genes were found to be up-regulated in tumor samples and cancer cell lines with high uPAR-del4/5 mRNA expression. The gene encoding rab31, a member of the Ras oncogene family, was selected for quantitative analysis of mRNA expression in the set of 280 patients. High rab31 values were significantly associated with worse outcome of patients for DMFS (HR = 2.27, P < 0.001) and OS (HR = 2.01, P = 0.008) in multivariate analysis, independent from uPAR-del4/5. The patient subgroup with high uPAR-del4/5 and rab31 levels showed the worst DMFS and OS (P < 0.001, both) compared with tumors with low values of both factors. Conclusions: Our results suggest that uPAR-del4/5 and rab31 mRNA represent independent prognostic markers in breast cancer and may be components of different, but possibly associated, tumor-relevant signaling pathways.

Gene expression profiles associated with the presence of a fibrotic focus and the growth pattern in lymph node-negative breast cancer (Article)

2008-05-15T00:00:00Z

Purpose: A fibrotic focus, the scar-like area found in the center of an invasive breast tumor, is a prognostic parameter associated with an expansive growth pattern, hypoxia, and (lymph) angiogenesis. Little is known about the molecular pathways involved. Experimental Design: Sixty-five patients were selected of whom microarray data of the tumor and H&E slides for histologic analysis were available. The growth pattern and the presence and size of a fibrotic focus were assessed. Differences in biological pathways were identified with global testing. The correlations of growth pattern and fibrotic focus with common breast cancer signatures and with clinicopathologic variables and survival were investigated. Results: Tumors with a large fibrotic focus showed activation of Ras signaling and of the hypoxia-inducible factor-1 α pathway. Furthermore, unsupervised hierarchical cluster analysis with hypoxia- and (lymph) angiogenesis-related genes showed that hypoxia-inducible factor-1 α vascular endothelial growth factor A, and carbonic anhydrase 9 were overexpressed. The presence of a fibrotic focus, especially a large fibrotic focus, was associated with the basal-like subtype (P = 0.009), an activated wound-healing signature (P = 0.06), and a poor-prognosis 76-gene signature (P = 0.004). The presence of a fibrotic focus (P = 0.02) and especially of a large fibrotic focus (P = 0.004) was also associated with early development of distant metastasis. Conclusions: Our results sustain the hypothesis that hypoxia-driven angiogenesis is essential in the biology of a fibrotic focus. Ras and Akt might play a role as downstream modulators. Our data furthermore suggest that vascular endothelial growth factor A does not only drive angiogenesis but also lymphangiogenesis in tumors with a fibrotic focus. Our data also show an association between the presence of a fibrotic focus and infaust molecular signatures.

Storing, linking, and mining microarray databases using SRS. (Article)

2005-07-27T00:00:00Z

BACKGROUND: SRS (Sequence Retrieval System) has proven to be a valuable platform for storing, linking, and querying biological databases. Due to the availability of a broad range of different scientific databases in SRS, it has become a useful platform to incorporate and mine microarray data to facilitate the analyses of biological questions and non-hypothesis driven quests. Here we report various solutions and tools for integrating and mining annotated expression data in SRS. RESULTS: We devised an Auto-Upload Tool by which microarray data can be automatically imported into SRS. The dataset can be linked to other databases and user access can be set. The linkage comprehensiveness of microarray platforms to other platforms and biological databases was examined in a network of scientific databases. The stored microarray data can also be made accessible to external programs for further processing. For example, we built an interface to a program called Venn Mapper, which collects its microarray data from SRS, processes the data by creating Venn diagrams, and saves the data for interpretation. CONCLUSION: SRS is a useful database system to store, link and query various scientific datasets, including microarray data. The user-friendly Auto-Upload Tool makes SRS accessible to biologists for linking and mining user-owned databases.

Molecular portrait of the progestagenic and estrogenic actions of tibolone: behavior of cellular networks in response to tibolone. (Article)

2005-02-01T00:00:00Z

Tibolone is a synthetic steroid with estrogenic effects on brain, vagina, and bone without stimulating the endometrium. During tibolone treatment, it is thought that the progestagenic properties of tibolone stimulate cell differentiation, which effectively counterbalances the growth-stimulating effects of the estrogenic properties of tibolone. The objective of this study was to characterize the expression profile that reflects the endometrial responses to the separated estrogenic (growth-inducing) and progestagenic (growth-inhibiting) actions of tibolone, thus gaining insight into the counteracting effect of these properties of tibolone on the endometrium. The estrogenic action of tibolone was studied in the estrogen-responsive ECC1 cell line (expressing estrogen receptor alpha), and the progestagenic action was studied in the progesterone-responsive cell line Ishikawa PRAB-36 (expressing PRA and PRB). The data showed that the progestagenic and estrogenic effects of tibolone produce different expression profiles with a narrow overlap in genes; however, both properties modulate the same biological processes. The final genetic network analysis indicated that the estrogenic effect of tibolone is potentially counterbalanced by the progestagenic metabolite of tibolone via differential regulation of similar cellular processes. For example, both progestagenic and estrogenic properties stimulate proliferation, but they exert the opposite effect on apoptosis. The apoptosis network was stimulated by the progestagenic properties of tibolone; in contrast, the estrogenic effect of tibolone suppressed the apoptosis network. The current results indicate that this differential regulation is realized through modulation of a different group of genes and rarely via contraregulation of the same set of genes.

Breast cancer oestrogen independence mediated by BCAR1 or BCAR3 genes is transmitted through mechanisms distinct from the oestrogen receptor signalling pathway or the epidermal growth factor receptor signalling pathway. (Article)

2005-01-01T00:00:00Z

INTRODUCTION: Tamoxifen is effective for endocrine treatment of oestrogen receptor-positive breast cancers but ultimately fails due to the development of resistance. A functional screen in human breast cancer cells identified two BCAR genes causing oestrogen-independent proliferation. The BCAR1 and BCAR3 genes both encode components of intracellular signal transduction, but their direct effect on breast cancer cell proliferation is not known. The aim of this study was to investigate the growth control mediated by these BCAR genes by gene expression profiling. METHODS: We have measured the expression changes induced by overexpression of the BCAR1 or BCAR3 gene in ZR-75-1 cells and have made direct comparisons with the expression changes after cell stimulation with oestrogen or epidermal growth factor (EGF). A comparison with published gene expression data of cell models and breast tumours is made. RESULTS: Relatively few changes in gene expression were detected in the BCAR-transfected cells, in comparison with the extensive and distinct differences in gene expression induced by oestrogen or EGF. Both BCAR1 and BCAR3 regulate discrete sets of genes in these ZR-75-1-derived cells, indicating that the proliferation signalling proceeds along distinct pathways. Oestrogen-regulated genes in our cell model showed general concordance with reported data of cell models and gene expression association with oestrogen receptor status of breast tumours. CONCLUSIONS: The direct comparison of the expression profiles of BCAR transfectants and oestrogen or EGF-stimulated cells strongly suggests that anti-oestrogen-resistant cell proliferation is not caused by alternative activation of the oestrogen receptor or by the epidermal growth factor receptor signalling pathway.

GO-Mapper: functional analysis of gene expression data using the expression level as a score to evaluate Gene Ontology terms. (Article)

2004-11-01T00:00:00Z

MOTIVATION: Retrieval of information on biological processes from large-scale expression data is still a time-consuming task. An automated analysis utilizing all expression information would greatly increase our understanding of the samples under study. RESULTS: We describe here a novel method to obtain a functional analysis of complex gene expression data. Instead of applying a predefined expression threshold, Gene Ontology (GO) terms are weighted using the actual measured levels of expression of all associated genes. Based on this concept, the application GO-Mapper was developed to quantitatively link gene expression levels to GO-terms for multiple experiments in an automated way. The applicability of GO-Mapper was developed and validated on in house and public human microarray data and mouse SAGE data. We demonstrate that the GO-Mapper allows for interrelating relevant biological functions with the experiments under study. AVAILABILITY: The GO-Mapper application is free of charge available from our website.