The life history of 21 breast cancers.

Cancer evolves dynamically as clonal expansions supersede one another driven by shifting selective pressures, mutational processes, and disrupted cancer genes. These processes mark the genome, such that a cancer's life history is encrypted in the somatic mutations present. We developed algorithms to decipher this narrative and applied them to 21 breast cancers. Mutational processes evolve across a cancer's lifespan, with many emerging late but contributing extensive genetic variation. Subclonal diversification is prominent, and most mutations are found in just a fraction of tumor cells. Every tumor has a dominant subclonal lineage, representing more than 50% of tumor cells. Minimal expansion of these subclones occurs until many hundreds to thousands of mutations have accumulated, implying the existence of long-lived, quiescent cell lineages capable of substantial proliferation upon acquisition of enabling genomic changes. Expansion of the dominant subclone to an appreciable mass may therefore represent the final rate-limiting step in a breast cancer's development, triggering diagnosis.

The circular RNome of primary breast cancer.

Circular RNAs (circRNAs) are a class of RNAs that is under increasing scrutiny, although their functional roles are debated. We analyzed RNA-seq data of 348 primary breast cancers and developed a method to identify circRNAs that does not rely on unmapped reads or known splice junctions. We identified 95,843 circRNAs, of which 20,441 were found recurrently. Of the circRNAs that match exon boundaries of the same gene, 668 showed a poor or even negative (R < 0.2) correlation with the expression level of the linear gene. In silico analysis showed only a minority (8.5%) of circRNAs could be explained by known splicing events. Both these observations suggest that specific regulatory processes for circRNAs exist. We confirmed the presence of circRNAs of CNOT2, CREBBP, and RERE in an independent pool of primary breast cancers. We identified circRNA profiles associated with subgroups of breast cancers and with biological and clinical features, such as amount of tumor lymphocytic infiltrate and proliferation index. siRNA-mediated knockdown of circCNOT2 was shown to significantly reduce viability of the breast cancer cell lines MCF-7 and BT-474, further underlining the biological relevance of circRNAs. Furthermore, we found that circular, and not linear, CNOT2 levels are predictive for progression-free survival time to aromatase inhibitor (AI) therapy in advanced breast cancer patients, and found that circCNOT2 is detectable in cell-free RNA from plasma. We showed that circRNAs are abundantly present, show characteristics of being specifically regulated, are associated with clinical and biological properties, and thus are relevant in breast cancer.

The prognostic and predictive value of ESR1 fusion gene transcripts in primary breast cancer.

BACKGROUND: In breast cancer (BC), recurrent fusion genes of estrogen receptor alpha (ESR1) and AKAP12, ARMT1 and CCDC170 have been reported. In these gene fusions the ligand binding domain of ESR1 has been replaced by the transactivation domain of the fusion partner constitutively activating the receptor. As a result, these gene fusions can drive tumor growth hormone independently as been shown in preclinical models, but the clinical value of these fusions have not been reported. Here, we studied the prognostic and predictive value of different frequently reported ESR1 fusion transcripts in primary BC. METHODS: We evaluated 732 patients with primary BC (131 ESR1-negative and 601 ESR1-positive cases), including two ER-positive BC patient cohorts: one cohort of 322 patients with advanced disease who received first-line endocrine therapy (ET) (predictive cohort), and a second cohort of 279 patients with lymph node negative disease (LNN) who received no adjuvant systemic treatment (prognostic cohort). Fusion gene transcript levels were measured by reverse transcriptase quantitative PCR. The presence of the different fusion transcripts was associated, in uni- and multivariable Cox regression analysis taking along current clinico-pathological characteristics, to progression free survival (PFS) during first-line endocrine therapy in the predictive cohort, and disease- free survival (DFS) and overall survival (OS) in the prognostic cohort. RESULTS: The ESR1-CCDC170 fusion transcript was present in 27.6% of the ESR1-positive BC subjects and in 2.3% of the ESR1-negative cases. In the predictive cohort, none of the fusion transcripts were associated with response to first-line ET. In the prognostic cohort, the median DFS and OS were respectively 37 and 93 months for patients with an ESR1-CCDC170 exon 8 gene fusion transcript and respectively 91 and 212 months for patients without this fusion transcript. In a multivariable analysis, this ESR1-CCDC170 fusion transcript was an independent prognostic factor for DFS (HR) (95% confidence interval (CI): 1.8 (1.2-2.8), P = 0.005) and OS (HR (95% CI: 1.7 (1.1-2.7), P = 0.023). CONCLUSIONS: Our study shows that in primary BC only ESR1-CCDC170 exon 8 gene fusion transcript carries prognostic value. None of the ESR1 fusion transcripts, which are considered to have constitutive ER activity, was predictive for outcome in BC with advanced disease treated with endocrine treatment.

Detection of tumor-derived extracellular vesicles in plasma from patients with solid cancer.

BACKGROUND: Extracellular vesicles (EVs) are actively secreted by cells into body fluids and contain nucleic acids of the cells they originate from. The goal of this study was to detect circulating tumor-derived EVs (ctEVs) by mutant mRNA transcripts (EV-RNA) in plasma of patients with solid cancers and compare the occurrence of ctEVs with circulating tumor DNA (ctDNA) in cell-free DNA (cfDNA). METHODS: For this purpose, blood from 20 patients and 15 healthy blood donors (HBDs) was collected in different preservation tubes (EDTA, BCT, CellSave) and processed into plasma within 24 h from venipuncture. EVs were isolated with the ExoEasy protocol from this plasma and from conditioned medium of 6 cancer cell lines and characterized according to MISEV2018-guidelines. RNA from EVs was isolated with the ExoRNeasy protocol and evaluated for transcript expression levels of 96 genes by RT-qPCR and genotyped by digital PCR. RESULTS: Our workflow applied on cell lines revealed a high concordance between cellular mRNA and EV-RNA in expression levels as well as variant allele frequencies for PIK3CA, KRAS and BRAF. Plasma CD9-positive EV and GAPDH EV-RNA levels were significantly different between the preservation tubes. The workflow detected only ctEVs with mutant transcripts in plasma of patients with high amounts (> 20%) of circulating tumor DNA (ctDNA). Expression profiling showed that the EVs from patients resemble healthy donors more than tumor cell lines supporting that most EVs are derived from healthy tissue. CONCLUSIONS: We provide a workflow for ctEV detection by spin column-based generic isolation of EVs and PCR-based measurement of gene expression and mutant transcripts in EV-RNA derived from cancer patients' blood plasma. This workflow, however, detected tumor-specific mutations in blood less often in EV-RNA than in cfDNA.

Optimization of Pancreatic Juice Collection: A First Step Toward Biomarker Discovery and Early Detection of Pancreatic Cancer.

INTRODUCTION: Imaging-based surveillance programs fail to detect pancreatic ductal adenocarcinoma at a curable stage, creating an urgent need for diagnostic biomarkers. METHODS: Secretin-stimulated pancreatic juice (PJ) was collected from the duodenal lumen during endoscopic ultrasound. The yield of biomarkers and organoids was compared for 2 collection techniques (endoscope suction channel vs catheter-based) and 3 periods (0-4 vs 4-8 vs 8-15 minutes). RESULTS: Collection through the endoscope suction channel was superior to collection with a catheter. Collection beyond 8 minutes reduced biomarker yield. PJ-derived organoid culture was feasible. DISCUSSION: The optimal protocol for secretin-stimulated PJ collection is through the endoscope suction channel for 8 minutes allowing biomarker detection and organoid culture.

Androgen receptor expression in circulating tumor cells of patients with metastatic breast cancer.

The androgen receptor (AR) has potential clinical relevance in metastatic breast cancer (mBC) since it might be a treatment target and has been associated with endocrine resistance. A minimal-invasive way to determine AR expression on metastatic tumor cells is by characterization of circulating tumor cells (CTCs). Here, we assessed AR mRNA expression in CTCs (CTC-AR) and in matched primary tumor samples from mBC patients representing different breast cancer subtypes. In addition, we explored CTC-AR-status in relation to outcome on endocrine therapy. AR, and 92 AR or estrogen receptor (ER) related genes, were measured in CellSearch-enriched CTCs from 124 mBC patients and in 52 matched FFPE primary tissues using quantitative reverse-transcriptase PCR. AR in CTCs was considered positive if the expression was 1 standard deviation higher than the expression measured in 11 healthy blood donors. A total of 31% of the mBC patients had AR-positive (AR+) CTCs. 58% of the matched CTC and primary tumor samples were discordant with respect to AR status, observing both switches from AR+ to AR-negative (AR-) and vice versa. There was no statistically significant difference in progression-free survival for patients treated with ER-targeting drugs and CTC-AR-status (13 AR+/ 37 AR- cases, p = 0.28). Thus, AR can be determined in RNA isolated from CTCs, with in our set 31% AR-positive samples. Given the discordance between AR status in CTC samples and corresponding primary tumors, determination of AR expression in CTCs might be a promising tool to select mBC patients for AR inhibiting agents.

Analysis of clonal expansions through the normal and premalignant human breast epithelium reveals the presence of luminal stem cells.

It is widely accepted that the cell of origin of breast cancer is the adult mammary epithelial stem cell; however, demonstrating the presence and location of tissue stem cells in the human breast has proved difficult. Furthermore, we do not know the clonal architecture of the normal and premalignant mammary epithelium or its cellular hierarchy. Here, we use deficiency in the mitochondrial enzyme cytochrome c oxidase (CCO), typically caused by somatic mutations in the mitochondrial genome, as a means to perform lineage tracing in the human mammary epithelium. PCR sequencing of laser-capture microdissected cells in combination with immunohistochemistry for markers of lineage differentiation was performed to determine the clonal nature of the mammary epithelium. We have shown that in the normal human breast, clonal expansions (defined here by areas of CCO deficiency) are typically uncommon and of limited size, but can occur at any site within the adult mammary epithelium. The presence of a stem cell population was shown by demonstrating multi-lineage differentiation within CCO-deficient areas. Interestingly, we observed infrequent CCO deficiency that was restricted to luminal cells, suggesting that niche succession, and by inference stem cell location, is located within the luminal layer. CCO-deficient areas appeared large within areas of ductal carcinoma in situ, suggesting that the rate of clonal expansion was altered in the premalignant lesion. © 2017 The Authors. The Journal of Pathology published by John Wiley & Sons Ltd on behalf of Pathological Society of Great Britain and Ireland.

Gene length corrected trimmed mean of M-values (GeTMM) processing of RNA-seq data performs similarly in intersample analyses while improving intrasample comparisons.

BACKGROUND: Current normalization methods for RNA-sequencing data allow either for intersample comparison to identify differentially expressed (DE) genes or for intrasample comparison for the discovery and validation of gene signatures. Most studies on optimization of normalization methods typically use simulated data to validate methodologies. We describe a new method, GeTMM, which allows for both inter- and intrasample analyses with the same normalized data set. We used actual (i.e. not simulated) RNA-seq data from 263 colon cancers (no biological replicates) and used the same read count data to compare GeTMM with the most commonly used normalization methods (i.e. TMM (used by edgeR), RLE (used by DESeq2) and TPM) with respect to distributions, effect of RNA quality, subtype-classification, recurrence score, recall of DE genes and correlation to RT-qPCR data. RESULTS: We observed a clear benefit for GeTMM and TPM with regard to intrasample comparison while GeTMM performed similar to TMM and RLE normalized data in intersample comparisons. Regarding DE genes, recall was found comparable among the normalization methods, while GeTMM showed the lowest number of false-positive DE genes. Remarkably, we observed limited detrimental effects in samples with low RNA quality. CONCLUSIONS: We show that GeTMM outperforms established methods with regard to intrasample comparison while performing equivalent with regard to intersample normalization using the same normalized data. These combined properties enhance the general usefulness of RNA-seq but also the comparability to the many array-based gene expression data in the public domain.

Gene Expression Differences between Ductal Carcinoma in Situ with and without Progression to Invasive Breast Cancer.

To understand the molecular alterations driving the progression of ductal carcinoma in situ (DCIS), we compared patients with pure DCIS and patients with DCIS and synchronous invasive breast cancer (IBC). Twelve patients with extensive pure DCIS were included as a representation of indolent lesions with limited invasive capacity. These cases were matched with 12 patients with a limited DCIS component and IBC, representing lesions with a high invasive potential. Matching included age and surrogate DCIS subtypes. Gene expression profiling was performed on DCIS cells to identify transcriptional differences between these two groups. The identified genes were validated by immunohistochemistry. Nine genes showed significantly different expression. Most of these genes were highly expressed in DCIS samples with IBC, including PLAU (P = 0.002), COL1A1 (P = 0.006), KRT81 (P = 0.009), S100A7 (P = 0.015), SCGB1D2 (P = 0.023), KRT18 (P = 0.029), and NOTCH3 (P = 0.044), whereas EGFR and CXCL14 showed a higher expression in cases with pure DCIS (P = 0.015 and P = 0.028, respectively). This difference was only significant for SCGB1D2 (P = 0.009). Hierarchical clustering revealed distinct clustering of patients with and without invasion. Patients with pure DCIS have a different gene expression pattern as compared to patients with DCIS and synchronous IBC. These genes may pinpoint to driver pathway(s) that play an important role in DCIS progression.

T lymphocytes facilitate brain metastasis of breast cancer by inducing Guanylate-Binding Protein 1 expression.

The discovery of genes and molecular pathways involved in the formation of brain metastasis would direct the development of therapeutic strategies to prevent this deadly complication of cancer. By comparing gene expression profiles of Estrogen Receptor negative (ER-) primary breast tumors between patients who developed metastasis to brain and to organs other than brain, we found that T lymphocytes promote the formation of brain metastases. To functionally test the ability of T cells to promote brain metastasis, we used an in vitro blood-brain barrier (BBB) model. By co-culturing T lymphocytes with breast cancer cells, we confirmed that T cells increase the ability of breast cancer cells to cross the BBB. Proteomics analysis of the tumor cells revealed Guanylate-Binding Protein 1 (GBP1) as a key T lymphocyte-induced protein that enables breast cancer cells to cross the BBB. The GBP1 gene appeared to be up-regulated in breast cancer of patients who developed brain metastasis. Silencing of GBP1 reduced the ability of breast cancer cells to cross the in vitro BBB model. In addition, the findings were confirmed in vivo in an immunocompetent syngeneic mouse model. Co-culturing of ErbB2 tumor cells with activated T cells induced a significant increase in Gbp1 expression by the cancer cells. Intracardial inoculation of the co-cultured tumor cells resulted in preferential seeding to brain. Moreover, intracerebral outgrowth of the tumor cells was demonstrated. The findings point to a role of T cells in the formation of brain metastases in ER- breast cancers, and provide potential targets for intervention to prevent the development of cerebral metastases.

CCAT2, a novel noncoding RNA mapping to 8q24, underlies metastatic progression and chromosomal instability in colon cancer.

The functional roles of SNPs within the 8q24 gene desert in the cancer phenotype are not yet well understood. Here, we report that CCAT2, a novel long noncoding RNA transcript (lncRNA) encompassing the rs6983267 SNP, is highly overexpressed in microsatellite-stable colorectal cancer and promotes tumor growth, metastasis, and chromosomal instability. We demonstrate that MYC, miR-17-5p, and miR-20a are up-regulated by CCAT2 through TCF7L2-mediated transcriptional regulation. We further identify the physical interaction between CCAT2 and TCF7L2 resulting in an enhancement of WNT signaling activity. We show that CCAT2 is itself a WNT downstream target, which suggests the existence of a feedback loop. Finally, we demonstrate that the SNP status affects CCAT2 expression and the risk allele G produces more CCAT2 transcript. Our results support a new mechanism of MYC and WNT regulation by the novel lncRNA CCAT2 in colorectal cancer pathogenesis, and provide an alternative explanation of the SNP-conferred cancer risk.

An 8-gene mRNA expression profile in circulating tumor cells predicts response to aromatase inhibitors in metastatic breast cancer patients.

BACKGROUND: Molecular characterization of circulating tumor cells (CTC) is promising for personalized medicine. We aimed to identify a CTC gene expression profile predicting outcome to first-line aromatase inhibitors in metastatic breast cancer (MBC) patients. METHODS: CTCs were isolated from 78 MBC patients before treatment start. mRNA expression levels of 96 genes were measured by quantitative reverse transcriptase polymerase chain reaction. After applying predefined exclusion criteria based on lack of sufficient RNA quality and/or quantity, the data from 45 patients were used to construct a gene expression profile to predict poor responding patients, defined as disease progression or death <9 months, by a leave-one-out cross validation. RESULTS: Of the 45 patients, 19 were clinically classified as poor responders. To identify them, the 75% most variable genes were used to select genes differentially expressed between good and poor responders. An 8-gene CTC predictor was significantly associated with outcome (Hazard Ratio [HR] 4.40, 95% Confidence Interval [CI]: 2.17-8.92, P < 0.001). This predictor identified poor responding patients with a sensitivity of 63% and a positive predictive value of 75%, while good responding patients were correctly predicted in 85% of the cases. In multivariate Cox regression analysis, including CTC count at baseline, the 8-gene CTC predictor was the only factor independently associated with outcome (HR 4.59 [95% CI: 2.11-9.56], P < 0.001). This 8-gene signature was not associated with outcome in a group of 71 MBC patients treated with systemic treatments other than AI. CONCLUSIONS: An 8-gene CTC predictor was identified which discriminates good and poor outcome to first-line aromatase inhibitors in MBC patients. Although results need to be validated, this study underscores the potential of molecular characterization of CTCs.

Integrative analysis of genomics and proteomics data on clinical breast cancer tissue specimens extracted with acid guanidinium thiocyanate-phenol-chloroform.

Acid guanidinium thiocyanate, phenol, and chloroform extraction (AGPC) is a commonly used procedure to extract RNA from fresh frozen tissues and cell lines. In addition, DNA and proteins can be recovered, which makes AGPC an attractive source for integrative analysis on tissues of which little material is available, such as clinical specimens. Despite this potential, AGPC has only scarcely been used for proteomic analysis, mainly due to difficulties in extracting proteins. We have used a quantitative mass spectrometry method to show that proteins can readily be recovered from AGPC extracted tissues with high recovery and repeatability, which allows this method to be used for global proteomic profiling. Protein expression data for a selected number of clinically relevant markers, of which transcript and protein levels are known to be correlated, were in agreement with genomic and transcriptomic data obtained from the same AGPC lysate. Furthermore, global proteomic profiling successfully discriminated breast tumor tissues according to their clinical subtype. Lastly, a reference gene set of differentially expressed transcripts was strongly enriched in the differentially abundant proteins in our cohort. AGPC lysates are therefore well suited for comparative protein and integrative analyses.

DC-SCRIPT is a novel regulator of the tumor suppressor gene CDKN2B and induces cell cycle arrest in ERα-positive breast cancer cells.

Breast cancer is one of the most common causes of cancer-related deaths in women. The estrogen receptor (ERα) is well known for having growth promoting effects in breast cancer. Recently, we have identified DC-SCRIPT (ZNF366) as a co-suppressor of ERα and as a strong and independent prognostic marker in ESR1 (ERα gene)-positive breast cancer patients. In this study, we further investigated the molecular mechanism on how DC-SCRIPT inhibits breast cancer cell growth. DC-SCRIPT mRNA levels from 190 primary ESR1-positive breast tumors were related to global gene expression, followed by gene ontology and pathway analysis. The effect of DC-SCRIPT on breast cancer cell growth and cell cycle arrest was investigated using novel DC-SCRIPT-inducible MCF7 breast cancer cell lines. Genome-wide expression profiling of DC-SCRIPT-expressing MCF7 cells was performed to investigate the effect of DC-SCRIPT on cell cycle-related gene expression. Findings were validated by real-time PCR in a cohort of 1,132 ESR1-positive breast cancer patients. In the primary ESR1-positive breast tumors, DC-SCRIPT expression negatively correlated with several cell cycle gene ontologies and pathways. DC-SCRIPT expression strongly reduced breast cancer cell growth in vitro, breast tumor growth in vivo, and induced cell cycle arrest. In addition, in the presence of DC-SCRIPT, multiple cell cycles related genes were differentially expressed including the tumor suppressor gene CDKN2B. Moreover, in 1,132 primary ESR1-positive breast tumors, DC-SCRIPT expression also correlated with CDKN2B expression. Collectively, these data show that DC-SCRIPT acts as a novel regulator of CDKN2B and induces cell cycle arrest in ESR1-positive breast cancer cells.

Complex landscapes of somatic rearrangement in human breast cancer genomes.

Multiple somatic rearrangements are often found in cancer genomes; however, the underlying processes of rearrangement and their contribution to cancer development are poorly characterized. Here we use a paired-end sequencing strategy to identify somatic rearrangements in breast cancer genomes. There are more rearrangements in some breast cancers than previously appreciated. Rearrangements are more frequent over gene footprints and most are intrachromosomal. Multiple rearrangement architectures are present, but tandem duplications are particularly common in some cancers, perhaps reflecting a specific defect in DNA maintenance. Short overlapping sequences at most rearrangement junctions indicate that these have been mediated by non-homologous end-joining DNA repair, although varying sequence patterns indicate that multiple processes of this type are operative. Several expressed in-frame fusion genes were identified but none was recurrent. The study provides a new perspective on cancer genomes, highlighting the diversity of somatic rearrangements and their potential contribution to cancer development.

The challenge of gene expression profiling in heterogeneous clinical samples.

Almost all samples used in tumor biology, such as tissues and bodily fluids, are heterogeneous, i.e., consist of different cell types. Evaluating the degree of heterogeneity in samples can increase our knowledge on processes such as clonal selection and metastasis. In addition, generating expression profiles from specific sub populations of cells can reveal their distinct functions. Tissue heterogeneity also poses a challenge, as it can confound the interpretation of gene expression data. This chapter will (1) give a brief overview on how heterogeneity may influence gene expression profiling data and (2) describe the methods that are currently available to assess transcriptional biomarkers in a heterogeneous cell population.

KRAS and BRAF mutation status in circulating colorectal tumor cells and their correlation with primary and metastatic tumor tissue.

Although anti-EGFR therapy has established efficacy in metastatic colorectal cancer, only 10-20% of unselected patients respond. This is partly due to KRAS and BRAF mutations, which are currently assessed in the primary tumor. To improve patient selection, assessing mutation status in circulating tumor cells (CTCs), which possibly better represent metastases than the primary tumor, could be advantageous. We investigated the feasibility of KRAS and BRAF mutation detection in colorectal CTCs by comparing three sensitive methods and compared mutation status in matching primary tumor, liver metastasis and CTCs. CTCs were isolated from blood drawn from 49 patients before liver resection using CellSearch™. DNA and RNA was isolated from primary tumors, metastases and CTCs. Mutations were assessed by co-amplification at lower denaturation temperature-PCR (Transgenomic™), real-time PCR (EntroGen™) and nested Allele-Specific Blocker (ASB-)PCR and confirmed by Sanger sequencing. In 43 of the 49 patients, tissue RNA and DNA was of sufficient quantity and quality. In these 43 patients, discordance between primary and metastatic tumor was 23% for KRAS and 7% for BRAF mutations. RNA and DNA from CTCs was available from 42 of the 43 patients, in which ASB-PCR was able to detect the most mutations. Inconclusive results in patients with low CTC counts limited the interpretation of discrepancies between tissue and CTCs. Determination of KRAS and BRAF mutations in CTCs is challenging but feasible. Of the tested methods, nested ASB-PCR, enabling detection of KRAS and BRAF mutations in patients with as little as two CTCs, seems to be superior.

Loss of E-cadherin is not a necessity for epithelial to mesenchymal transition in human breast cancer.

Epithelial to mesenchymal transition (EMT) is typically defined by the acquisition of a spindle cell morphology in combination with loss of E-cadherin and upregulation of mesenchymal markers. However, by studying E-cadherin inactivation in 38 human breast cancer cell lines, we noted that not all cell lines that had undergone EMT had concomitantly lost E-cadherin expression. We further investigated this discrepancy functionally and in clinical breast cancer specimens. Interestingly, reconstitution of wild-type E-cadherin cDNA in a E-cadherin negative cell line that had undergone EMT (MDA-MB-231) did not revert the spindle morphology back to an epithelial morphology. Neither were changes observed in the expression of several markers known to be involved in the EMT process. Similarly, upregulation of E-cadherin via global DNA demethylation in eleven cell lines that had undergone EMT did not induce a change in cell morphology, nor did it alter the expression of EMT markers in these cells. Next, we extracted genes differentially expressed between cell lines that had undergone EMT versus cell lines that had not undergone EMT. Caveolin-1 was identified to be an excellent marker for EMT, irrespective of E-cadherin status (specificity and sensitivity of 100 %). Consistent with our observations in the breast cancer cell lines, expression of Caveolin-1 identified a subset of basal breast cancers, particularly of metaplastic pathology, and only 50 % of these lacked E-cadherin expression. The discrepancy between E-cadherin loss and EMT was thus reproduced in clinical samples. Together, these results indicate that in human breast cancer loss of E-cadherin is not causal for EMT and even not a necessity.

Ancestry-shift refinement mapping of the C6orf97-ESR1 breast cancer susceptibility locus.

We used an approach that we term ancestry-shift refinement mapping to investigate an association, originally discovered in a GWAS of a Chinese population, between rs2046210[T] and breast cancer susceptibility. The locus is on 6q25.1 in proximity to the C6orf97 and estrogen receptor alpha (ESR1) genes. We identified a panel of SNPs that are correlated with rs2046210 in Chinese, but not necessarily so in other ancestral populations, and genotyped them in breast cancer case:control samples of Asian, European, and African origin, a total of 10,176 cases and 13,286 controls. We found that rs2046210[T] does not confer substantial risk of breast cancer in Europeans and Africans (OR = 1.04, P = 0.099, and OR = 0.98, P = 0.77, respectively). Rather, in those ancestries, an association signal arises from a group of less common SNPs typified by rs9397435. The rs9397435[G] allele was found to confer risk of breast cancer in European (OR = 1.15, P = 1.2 x 10(-3)), African (OR = 1.35, P = 0.014), and Asian (OR = 1.23, P = 2.9 x 10(-4)) population samples. Combined over all ancestries, the OR was 1.19 (P = 3.9 x 10(-7)), was without significant heterogeneity between ancestries (P(het) = 0.36) and the SNP fully accounted for the association signal in each ancestry. Haplotypes bearing rs9397435[G] are well tagged by rs2046210[T] only in Asians. The rs9397435[G] allele showed associations with both estrogen receptor positive and estrogen receptor negative breast cancer. Using early-draft data from the 1,000 Genomes project, we found that the risk allele of a novel SNP (rs77275268), which is closely correlated with rs9397435, disrupts a partially methylated CpG sequence within a known CTCF binding site. These studies demonstrate that shifting the analysis among ancestral populations can provide valuable resolution in association mapping.

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.

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.