Search

In contrast to primary colorectal cancer (CRC) little is known about the genomic landscape of metastasized CRC. Here we present whole genome sequencing data of metastases of 429 CRC patients participating in the pan-cancer CPCT-02 study (NCT01855477). Unsupervised clustering using mutational signature patterns highlights three major patient groups characterized by signatures known from primary CRC, signatures associated with received prior treatments, and metastasis-specific signatures. Compared to primary CRC, we identify additional putative (non-coding) driver genes and increased frequencies in driver gene mutations. In addition, we identify specific genes preferentially affected by microsatellite instability. CRC-specific 1kb-10Mb deletions, enriched for common fragile sites, and LINC00672 mutations are associated with response to treatment in general, whereas FBXW7 mutations predict poor response specifically to EGFR-targeted treatment. In conclusion, the genomic landscape of mCRC shows defined changes compared to primary CRC, is affected by prior treatments and contains features with potential clinical relevance.

New therapies are an urgent medical need in all breast cancer subgroups. Metabotropic glutamate receptor 1 (mGluR1) is suggested as a potential new molecular target. We e

INTRODUCTION: Local T-cell immunity is recognized for its contribution to the evolution and therapy response of various carcinomas. Here, we investigated characteristics of tumor-infiltrating lymphocytes (TILs), as well as T-cell evasive mechanisms in different soft tissue sarcoma (STS) subtypes. METHODS: Liposarcoma, gastrointestinal stromal tumor (GIST), leiomyosarcoma, myxofibrosarcoma and pleomorphic sarcomas were assessed for T-cell numbers and phenotypes using flow cytometry. Next-generation sequencing was used to analyze T-cell receptor repertoire, mutational load, immune cell frequencies, and expression of immune-related genes. RESULTS: GIST, myxofibrosarcoma and pleomorphic sarcoma showed high numbers of CD8+ TILs, with GIST having the lowest fraction of effector memory T cells. These TILs coexpress the immune checkpoints PD1, TIM3, and LAG3 in myxofibrosarcoma and pleomorphic sarcoma, yet TILs coexpressing these checkpoints were near negligible in GIST. Fractions of dominant T-cell clones among STS subtypes were lowest in GIST and liposarcoma, whereas mutational load was relatively low in all STS subtypes. Furthermore, myeloid-derived cells and expression of the costimulatory ligands CD86, ICOS-L and 41BB-L were lowest in GIST when compared with other STS subtypes. CONCLUSION: STS subtypes differ with respect to number and phenotypical signs of antitumor responsiveness of CD8+ TILs. Notably, GIST, myxofibrosarcoma and pleomorphic sarcoma harbor high numbers of CD8+ T cells, yet in the GIST microenvironment, these T cells are less differentiated and non-exhausted, which is accompanied with a relatively low expression of costimulatory ligands.

BACKGROUND: There are profound individual differences in clinical outcomes between colorectal cancers (CRCs) presenting with identical stage of disease. Molecular stratification, in conjunction with the traditional TNM staging, is a promising way to predict patient outcomes. We investigated the interconnectivity between tumor stage and tumor biology reflected by the Consensus Molecular Subtypes (CMSs) in CRC, and explored the possible value of these insights in patients with stage II colon cancer. METHODS: We performed a retrospective analysis using clinical records and gene expression profiling in a meta-cohort of 1040 CRC patients. The interconnectivity of tumor biology and disease stage was assessed by investigating the association between CMSs and TNM classification. In order to validate the clinical applicability of our findings we employed a meta-cohort of 197 stage II colon cancers. RESULTS: CMS4 was significantly more prevalent in advanced stages of disease (stage I 9.8% versus stage IV 38.5%, p < 0.001). The observed differential gene expression between cancer stages is at least partly explained by the biological differences as reflected by CMS subtypes. Gene signatures for stage III-IV and CMS4 were highly correlated (r = 0.77, p < 0.001). CMS4 cancers showed an increased progression rate to more advanced stages (CMS4 compared to CMS2: 1.25, 95% CI: 1.08-1.46). Patients with a CMS4 cancer had worse survival in the high-risk stage II tumors compared to the total stage II cohort (5-year DFS 41.7% versus 100.0%, p = 0.008). CONCLUSIONS: Considerable interconnectivity between tumor biology and tumor stage in CRC exists. This implies that the TNM stage, in addition to the stage of progression, might also reflect distinct biological disease entities. These insights can potentially be utilized to optimize identification of high-risk stage II colon

The genetically heterogeneous triple-negative breast cancer (TNBC) continues to be an intractable disease, due to lack of effective targeted therapies. Gene amplification is a major event in tumorigenesis. Genes with amplification-dependent expression are being explored as therapeutic targets for cancer treatment. In this study, we have applied Analytical Multi-scale Identification of Recurring Events analysis and transcript quantification in the TNBC genome across 222 TNBC tumors and identified 138 candidate genes with positive correlation in copy number gain (CNG) and gene expression. siRNA-based loss-of-function screen of the candidate genes has validated EGFR, MYC, ASAP1, IRF2BP2, and CCT5 genes as drivers promoting proliferation in different TNBC cells. MYC, ASAP1, IRF2BP2, and CCT5 display frequent CNG and concurrent expression over 2173 breast cancer tumors (cBioPortal dataset). More frequently are MYC and ASAP1 amplified in TNBC tumors (>30%, n = 320). In particular, high expression of ASAP1, the ADP-ribosylation factor GTPase-activating protein, is significantly related to poor metastatic relapse-free survival of TNBC patients (n = 257, bc-GenExMiner). Furthermore, we have revealed that silencing of ASAP1 modulates numerous cytokine and apoptosis signaling components, such as IL1B, TRAF1, AIFM2, and MAP3K11 that are clinically relevant to survival outcomes of TNBC patients. ASAP1 has been reported to promote invasion and metastasis in various cancer cells. Our findings that ASAP1 is an amplification-dependent TNBC driver gene promoting TNBC cell proliferation, functioning upstream apoptosis components, and correlating to clinical outcomes of TNBC patients, support ASAP1 as a potential actionable target for TNBC treatment.

Background: Breast cancer (BC) immune infiltrates play a critical role in tumor progression and response to treatment. Besides stromal tumor infiltrating lymphocytes (sTILs) which have recently reached level 1B evidence as a prognostic marker in triple negative BC, a plethora of methods to assess immune infiltration exists, and it is unclear how these compare to each other and if they can be used interchangeably. Methods: Two experienced pathologists scored sTIL, intra-tumoral TIL (itTIL), and 6 immune cell types (CD3+, CD4+, CD8+, CD20+, CD68+, FOXP3+) in the International Cancer Genomics Consortium breast cancer cohort using hematoxylin and eosin-stained (n = 243) and immunohistochemistry-stained tissue microarrays (n = 254) and whole slides (n = 82). The same traits were evaluated using transcriptomic- and methylomic-based deconvolution methods or signatures. Results: The concordance correlation coefficient (CCC) between pathologists for sTIL was very good (0.84) and for cell-specific immune infiltrates slightly lower (0.63-0.66). Comparison between tissue microarray and whole slide pathology scores revealed systematically higher values in whole slides (ratio 2.60-5.98). The Spearman correlations between microscopic sTIL and transcriptomic- or methylomic-based assessment of immune infiltr

Despite the success of imatinib in advanced gastrointestinal stromal tumor (GIST) patients, 50% of the patients experience resistance within two years of treatment underscoring the need to get better insight into the mechanisms conferring imatinib resistance. Here the microRNA and mRNA expression profiles in primary (imatinib-naïve) and imatinib-resistant GIST were examined. Fifty-three GIST samples harboring primary KIT mutations (exon 9; n = 11/exon 11; n = 41/exon 17; n = 1) and comprising imatinib-naïve (IM-n) (n = 33) and imatinib-resistant (IM-r) (n = 20) tumors, were analyzed. The microRNA expression profiles were determined and from a subset (IM-n, n = 14; IM-r, n = 15) the mRNA expression profile was established. Ingenuity pathway analyses were used to unravel biochemical pathways and gene networks in IM-r GIST. Thirty-five differentially expressed miRNAs between IM-n and IM-r GIST samples were identified. Additionally, miRNAs distinguished IM-r samples with and without secondary KIT mutations. Furthermore 352 aberrantly expressed genes were found in IM-r samples. Pathway and network analyses revealed an association of differentially expressed genes with cell cycle progression and cellular proliferation, thereby implicating genes and pathways involved in imatinib resistance in GIST. Differentially expressed miRNAs and mRNAs between IM-n and IM-r GIST were identified. Bioinformatic analyses provided insight into the genes and biochemical pathways involved in imatinib-resistance and highlighted key genes that may be putative treatment targets.

(1) Background: Markers identifying which patients with metastatic, castration-resistant prostate cancer (mCRPC) will benefit from cabazitaxel therapy are currently lacking. Therefore, the aim of this study was to identify markers associated with outcome to cabazitaxel therapy based on counts and gene expression profiles of circulating tumor cells (CTCs). (2) Methods: From 120 mCRPC patients, CellSearch enriched CTCs were obtained at baseline and after 6 weeks of cabazitaxel therapy. Furthermore, 91 genes associated with prostate cancer were measured in mRNA of these CTCs. (3) Results: In 114 mCRPC patients with an evaluable CTC count, the CTC count was independently associated with poor progression-free survival (PFS) and overall survival (OS) in multivariable analysis with other commonly used variables associated with outcome in mCRPC (age, prostate specific antigen (PSA), alkaline phosphatase, lactate dehydrogenase (LDH), albumin, hemoglobin), together with alkaline phosphatase and hemoglobin. A five-gene expression profile was generated to predict for outcome to cabazitaxel therapy. However, even though this signature was associated with OS in univariate analysis, this was not the case in the multivariate analysis for OS nor for PFS. (4) Conclusion: The established five-gene expression profile in CTCs was not independently associated with PFS nor OS. However, along with alkaline phosphatase and hemoglobin, CTC-count is independently associated with PFS and OS in mCRPC patients who are treated with cabazitaxel.

Triple-negative breast cancer (TNBC) is an aggressive subtype of breast cancer with poor clinical prognosis and limited targeted treatment strategies. Kinase inhibitor screening of a panel of 20 TNBC cell lines uncovered three critical TNBC subgroups: 1) sensitive to only MEK inhibitors; 2) sensitive to only Akt inhibitors; 3) resistant to both MEK/Akt inhibitors. Using genomic, transcriptomic and proteomic datasets of these TNBC cell lines we unravelled molecular features associated with the MEK and Akt drug resistance. MEK inhibitor-resistant TNBC cell lines were discriminated from Akt inhibitor-resistant lines by the presence of PIK3CA/PIK3R1/PTEN mutations, high p-Akt and low p-MEK levels, yet these features could not distinguish double-resistant cells. Gene set enrichment analyses of transcriptomic and proteomic data of the MEK and Akt inhibitor response groups revealed a set of cell cycle-related genes associated with the double-resistant phenotype; these genes were overexpressed in a subset of breast cancer patients. CDK inhibitors targeting the cell cycle programme could overcome the Akt and MEK inhibitor double-resistance. In conclusion, we uncovered molecular features and alternative treatment strategies for TNBC that are double-resistant to Akt and MEK inhibitors.

Background: The effective treatment of triple-negative breast cancer (TNBC) remains a profound clinical challenge. Despite frequent epidermal growth factor receptor (EGFR) overexpression and reliance on downstream signalling pathways in TNBC, resistance to EGFR-tyrosine kinase inhibitors (TKIs) remains endemic. Therefore, the identification of targeted agents, which synergise with current therapeutic options, is paramount. Methods: Compound-based, high-throughput, proliferation screening was used to profile the response of TNBC cell lines to EGFR-TKIs, western blotting and siRNA transfection being used to examine the effect of inhibitors on EGFR-mediated signal transduction and cellular dependence o

Ttriple-negative breast cancer (TNBC) is an aggressive and highly metastatic breast cancer subtype. Enhanced TNBC cell motility is a prerequisite of TNBC cell dissemination. Here, we apply an imaging-based RNAi phenotypic cell migration screen using two highly motile TNBC cell lines (Hs578T and MDA-MB-231) to provide a repository of signaling determinants that functionally drive TNBC cell motility. We have screened ~4,200 target genes individually and discovered 133 and 113 migratory modulators of Hs578T and MDA-MB-231, respectively, which are linked to signaling networks predictive for breast cancer progression. The splicing factors PRPF4B and BUD31 and the transcription factor BPTF are essential for cancer cell migration, amplified in human primary breast tumors and associated with metastasis-free survival. Depletion of PRPF4B, BUD31 and BPTF causes primarily down regulation of genes involved in focal adhesion and ECM-interaction pathways. PRPF4B is essential for TNBC metastasis formation in vivo, making PRPF4B a candidate for further drug development.

MicroRNAs (miRNAs) regulate gene expression post-transcriptionally. In this way they might influence whether a cell is sensitive or resistant to a certain drug. So far, only a limited number of relatively small scale studies comprising few cell lines and/or drugs have been performed. To obtain a broader view on miRNAs and their association with drug response, we investigated the expression levels of 411 miRNAs in relation to drug sensitivity in 36 breast cancer cell lines. For this purpose IC50 values of a drug screen involving 34 drugs were associated with miRNA expression data of the same breast cancer cell lines. Since molecular subtype of the breast cancer cell lines is considered a confounding factor in drug association studies, multivariate analysis taking subtype into account was performed on significant miRNA-drug associations which retained 13 associations. These associations consisted of 11 different miRNAs and eight different drugs (among which Paclitaxel, Docetaxel and Veliparib). The taxanes, Paclitaxel and Docetaxel, were the only drugs having miRNAs in common: hsa-miR-187-5p and hsa-miR-106a-3p indicative of drug resistance while Paclitaxel sensitivity alone associated with hsa-miR-556-5p. Tivantinib was associated with hsalet-7d-5p and hsa-miR-18a-5p for sensitivity and hsa-miR-637 for resistance. Drug sensitivity was associated with hsa-let-7a-5p for Bortezomib, hsa-miR-135a-3p for JNJ-707 and hsa-miR-185-3p for Panobinostat. Drug resistance was associated with hsa-miR-182-5p for Veliparib and hsa-miR-629-5p for Tipifarnib. Pathway analysis for significant miRNAs was performed to reveal biological roles, aiding to find a potential mechanistic link for the observed associations with drug response. By doing so hsa-miR-187-5p was linked to the cell cycle G2-M checkpoint in line with this checkpoint being the target of taxanes. In conclusion, our study shows that miRNAs could potentially serve as biomarkers for intrinsic drug resistance and that pathway

Splicing factors (SFs) act in dynamic macromolecular complexes to modulate RNA processing. To understand the complex role of SFs in cancer progression, we performed a systemic analysis of the co-regulation of SFs using primary tumor RNA sequencing data. Co-regulated SFs were associated with aggressive breast cancer phenotypes and enhanced metastasis formation, resulting in the classification of Enhancer- (21 genes) and Suppressor-SFs (64 genes). High Enhancer-SF levels were related to distinct splicing patterns and expression of known oncogenic pathways such as respiratory electron transport, DNA damage and cell cycle regulation. Importantly, largely identical SF co-regulation was observed in almost all major cancer types, including lung, pancreas and prostate cancer. In conclusion, we identified cancer-associated co-regulated expression of SFs that are associated with aggressive phenotypes. This study increases the global understanding of the role of the spliceosome in cancer progression and also contributes to the development of strategies to cure cancer patients.

Global loss of DNA methylation and CpG island (CGI) hypermethylation are key epigenomic aberrations in cancer. Global loss manifests itself in partially methylated domains (PMDs) which extend up to megabases. However, the distribution of PMDs within and between tumor types, and their effects on key functional genomic elements including CGIs are poorly defined. We comprehensively show that loss of methylation in PMDs occurs in a large fraction of the genome and represents the prime source of DNA methylation variation. PMDs are hypervariable in methylation level, size and distribution, and display elevated mutation rates. They impose intermediate DNA methylation levels incognizant of functional genomic elements including CGIs, underpinning a CGI methylator phenotype (CIMP). Repression effects on tumor suppressor genes are negligible as they are generally excluded from PMDs. The genomic distribution of PMDs reports tissue-of-origin and may represent tissue-specific silent regions which tolerate instability at the epigenetic, transcriptomic and genetic level.

There is an increasing awareness of the importance of tumor – immune cell interactions to the evolution and therapy responses of breast cancer (BC). Not surprisingly, numerous studies are currently assessing the clinical value of immune modulation for BC patients. However, till now durable clinical responses are only rarely observed. It is important to realiz

There is an increasing awareness of the importance of tumor – immune cell interactions to the evolution and therapy responses of breast cancer (BC). Not surprisingly, numerous studies are currently assessing the clinical value of immune modulation for BC patients. However, till now durable clinical responses are only rarely observed. It is important to realiz

The human mitochondrial DNA (mtDNA) encodes 37 genes, including thirteen proteins essential for the respiratory chain, and RNAs functioning in the mitochondrial translation apparatus. The total number of mtDNA molecules per cell (mtDNA content) is variable between tissue types and also between tumors and their normal counterparts. For breast cancer, tumors tend to be depleted in their mtDNA content compared to adjacent normal mammary tissue. Various studies have shown that primary breast tumors harbor somatic mtDNA variants. A decrease in mtDNA content or the presence of somatic variants could indicate a reduced mitochondrial function within breast cancer. In this explorative study we aimed to further understand genomic changes and expression of the mitochondrial genome within breast cancer, by analyzing RNA sequencing data of primary breast tumor specimens of 344 cases. We demonstrate that somatic variants detected at the mtRNA level are representative for somatic variants in the mtDNA. Also, the number of somatic variants within the mitochondrial transcriptome is not associated with mutational processes impacting the nuclear genome, but is positi

The large number of non-coding RNAs (ncRNAs) and their breadth of functionalities has fuelled many studies on their roles in cancer. We previously linked four microRNAs to breast cancer prognosis. One of these microRNAs, hsa-miR-7, was found to be regulated by another type of ncRNA, the circular non-coding RNA (circRNA) CDR1-AS, which contains multiple hsa-miR-7 binding sites. Based on this finding, we studied the potential clinical value of this circRNA on breast cancer prognosis in a cohort based on a cohort that was previously analysed for hsa-miR-7 and in an adjuvant hormone-naïve cohort for 1st-line tamoxifen treatment outcomes, in which we also analysed hsa-miR-7. A negative correlation was observed between hsa-miR-7 and CDR1-AS in both cohorts. Despite associations with various clinical metrics (e.g., tumour grade, tumour size, and relapse location), CDR1-AS was neither prognostic nor predictive of relevant outcomes in our cohorts. However, we did observe stromal CDR1-AS expression, suggesting a possible cell-type specific interaction. Next to the known association of hsa-miR-7 expression with poor prognosis in primary breast cancer, we found that high hsa-miR-7 expression was predictive of an adverse response to tamoxifen therapy and poor progression-free and post-relapse overall survival in patients with recurrent disease.

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.

The identification of patients with high-risk stage II colon cancer who may benefit from adjuvant therapy may allow the clinical approach to be tailored for these patients based on an understanding of tumour biology. MicroRNAs have been proposed as markers of the prognosis or treatment response in colorectal cancer. Recently, a 2-microRNA signature (l et-7i and miR-10b) was proposed to identify colorectal cancer patients at risk of developing distant metastasis. We assessed the prognostic value of this signature and additional candidate microRNAs in an independent, clinically well-defined, prospectively collected cohort of primary colon cancer patients including stage I-II colon cancer without and stage III colon cancer with adjuvant treatment. The 2-microRNA signature specifically predicted hepatic recurrence in the stage I-II group, but not the overall ability to develop distant metastasis. The addition of miR-30b to the 2-microRNA signature allowed the prediction of both distant metastasis and hepatic recurrence in patients with stage I-II colon cancer who did not receive adjuvant chemotherapy. Available gene expression data allowed us to associate m iR-30b expression with axon guidance and l et-7i expression with cell adhesion, migration, and motility.

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

Genomic rearrangements that give rise to oncogenic gene fusions can offer actionable targets for cancer therapy. Here we present a systematic analysis of oncogenic gene fusions among a clinically well-characterized, prospectively collected set of 278 primary colon cancers spanning diverse tumor stages and clinical outcomes. Gene fusions and somatic genetic variations were identified in fresh frozen clinical specimens by Illumina RNA-sequencing, the STAR fusion gene detection pipeline, and GATK RNA-seq variant calling. We considered gene fusions to be pathogenically relevant when recurrent, producing divergent gene expression (outlier analysis), or as functionally important (e.g., kinase fusions). Overall, 2.5% of all specimens were defined as harboring a relevant gene fusion (kinase fusions 1.8%). Novel configurations of BRAF, NTRK3, and RET gene fusions resulting from chromosomal translocations were identified. An R-spondin fusion was found in only one tumor (0.35%), much less than an earlier reported frequency of 10% in colorectal cancers. We also found a novel fusion involving USP9X-ERAS formed by chromothripsis and leading to high expression of ERAS, a constitutively active RAS protein normally expressed only in embryonic stem cells. This USP9X–ERAS fusion appeared highly oncogenic on the basis of its ability to activate AKT signaling. Oncogenic fusions were identified only in lymph node–negative tumors that lacked BRAF or KRAS mutations. In summary, we identified several novel oncogenic gene fusions in colorectal cancer that may drive malignant development and offer new targets for personalized therapy.

Somatic rearrangements contribute to the mutagenized landscape of cancer genomes. Here, we systematically interrogated rearrangements in 560 breast cancers by using a piecewise constant fitting approach. We identified 33 hotspots of large (>100 kb) tandem duplications, a mutational signature associated with homologous-recombination-repair deficiency. Notably, these tandem-duplication hotspots were enriched in breast cancer germline susceptibility loci (odds ratio (OR) = 4.28) and breast-specific 'super-enhancer' regulatory elements (OR = 3.54). These hotspots may be

In a previous study, we detected a significant association between phosphoserine aminotransferase 1 (PSAT1) hyper-methylation and mRNA levels to outcome to tamoxifen treatment in recurrent disease. We here aimed to study the association of PSAT1 protein levels to outcome upon tamoxifen treatment and to obtain more insight in its role in tamoxifen resistance. A cohort of ER positive, hormonal therapy naïve primary breast carcinomas was immunohistochemically (IHC) stained for PSAT1. Staining was analyzed for association with patient's time to progression (TTP) and overall response on first-line tamoxifen for recurrent disease. PSAT1 mRNA levels were also assessed by reverse transcriptase quantitative polymerase chain reaction (RT-qPCR; n = 161) and Affymetrix GeneChip (n = 155). Association of PSAT1 to biological pathways on tamoxifen outcome were assessed by global test. PSAT1 protein and mRNA levels were significantly associated to poor outcome to tamoxifen treatment. When comparing PSAT1 p

Both healthy and cancerous breast tissue is heterogeneous, which is a bottleneck for proteomics-based biomarker analysis, as it obscures the cellular origin of a measured protein. We therefore aimed at obtaining a protein-level interpretation of malignant transformation through global proteome analysis of a variety of laser capture microdissected cells originating from benign and malignant breast tissues. We compared proteomic differences between these tissues, both from cells of epithelial origin and the stromal environment, and performed string analysis. Differences in protein abundances corresponded with several hallmarks of cancer, including loss of cell adhesion, transformation to a migratory phenotype, and enhanced energy metabolism. Furthermore, despite enriching for (tumor) epithelial cells, many changes to the extracellular matrix were detected in microdissected cells of epithelial origin. The stromal compartment was heterogeneous and richer in the number of fibroblast and immune cells in malignant sections, compared to benign tissue sections. Furthermore, stroma could be clearly divided into reactive and nonreactive based on extracellular matrix disassembly proteins. We conclude that proteomics analysis of both microdissected epithelium and stroma gives an additional layer of information and more detailed insight into malignant transformation.

OBJECTIVE: Overall and splice specific expression of Spleen Tyrosine Kinase (SYK) has been posed as a marker predicting both poor and favorable outcome in various epithelial malignancies. However, its role in colorectal cancer is largely unknown. The aim of this study was to explore the prognostic role of SYK in three cohorts of colon cancer patients.METHODS: Total messenger RNA (mRNA) expression of SYK, SYK(T), and mRNA expression of its two splice variants SYK short (S) and SYK long (L) were measured using quantitative reverse transcriptase (RT-qPCR) in 240 primary colon cancer patients (n = 160 patients with chemonaive lymph node negative [LNN] and n = 80 patients with adjuvant treated lymph node positive [LNP] colon cancer) and related to microsatellite instability (MSI), known colorectal cancer mutations, and disease-free (DFS), hepatic metastasis-free (HFS) and overall survival (OS). Two independent cohorts of patients with respectively 48 and 118 chemonaive LNN colon cancer were used for validation.RESULTS: Expression of SYK and its splice variants was significantly lower in tumors with MSI, and in KRAS wild type, BRAF mutant and PTEN mutant tumors. In a multivariate Cox regression analysis, as a continuous variable, increasing SYK(S) mRNA expression was associated with worse HFS (Hazard Ratio[HR] = 1.83; 95% Confidence Interval[

Background: Identification of specific risk groups for recurrence after surgery for isolated colorectal liver metastases (CRLM) remains challenging due to the heterogeneity of the disease. Classical clinicopathologic parameters have limited prognostic value. The aim of this study was to identify a gene expression signature measured in CRLM discriminating early from late recurrence after partial hepatectomy. Methods: CRLM from two patient groups were collected: I) with recurrent disease ≤12 months after surgery (N = 33), and II) without recurrences and disease free for ≥36 months (N = 30). The patients were clinically homogeneous; all had a low clinical risk score (0-2) and did not receive (neo-) adjuvant chemotherapy. Total RNA was hybridised to Illumina arrays, and processed for analysis. A leave-one-out cross validation (LOOCV) analysis was performed to identify a prognostic gene expression signature. Results: LOOCV yielded an 11-gene profile with prognostic value in relation to recurrent disease ≤12 months after partial hepatectomy. This signature had a sensitivity of 81.8%, with a specificity of 66.7% for predicting recurrences (≤12 months) versus no recurrences for at least 36 months after surgery (X2 P < 0.0001). Conclusion: The current study yielded an 11-gene signature at mRNA level in CRLM discriminating early from late or no relapse after partial hepatectomy.

Estrogen receptor (ER) positive tumors represent the majority of breast malignancies, and are effectively treated with hormonal therapies, such as tamoxifen. However, in the recurrent disease resistance to tamoxifen therapy is common and a major cause of death. In recent years, in-depth proteome analyses have enabled identification of clinically useful biomarkers, particularly, when heterogeneity in complex tumor tissue was reduced using laser capture microdissection (LCM). In the current study, we performed high resolution proteomic analysis on two cohorts of ER positive breast tumors derived from patients who either manifested good or poor outcome to tamoxifen treatment upon recurrence. A total of 112 fresh frozen tumors were collected from multiple medical centers and divided into two sets: an in-house training and a multi-center test set. Epithelial tumor cells were enriched with LCM and analyzed by nano-LC Orbitrap mass spectrometry (MS), which yielded >3000 and >4000 quantified proteins in the training and test sets, respectively. Raw data are available via ProteomeXchange with identifiers PXD000484 and PXD000485. Statistical analysis showed differential abundance of 99 proteins, of which a subset of 4 proteins was selected through a multivariate step-down to develop a predictor for tamoxifen treatment outcome. The 4-protein signature significantly predicted poor outcome patients in the test set, independent of predictive histopathological characteristics (hazard ratio [HR] = 2.17; 95% confidence interval [CI] = 1.15 to 4.17; multivariate Cox regression p value = 0.017). Immunohistochemical (IHC) staining of PDCD4, one of the signature proteins, on an independent set of formalin-fixed paraffin-embedded tumor tissues provided and independent technical validation (HR = 0.72; 95% CI = 0.57 to 0.92; multivariate Cox regression p value = 0.009). We hereby report the first validated protein predictor for tamoxifen treatment outcome in recurrent ER-positive breast ca

We analysed whole-genome sequences of 560 breast cancers to advance understanding of the driver mutations conferring clonal advantage and the mutational processes generating somatic mutations. We found that 93 protein-coding cancer genes carried probable driver mutations. Some non-coding regions exhibited high mutation frequencies, but most have distinctive structural features probably causing elevated mutation rates and do not contain driver mutations. Mutational signature analysis was extended to genome rearrangements and revealed twelve base substitution and six rearrangement signatures. Three rearrangement signatures, characterized by tandem duplications or deletions, appear associated with defective homologous-recombination-based DNA repair: one with deficient BRCA1 function, another with deficient BRCA1 or BRCA2 function, the cause of the third is unknown. This analysis of all classes of somatic mutation across exons, introns and intergenic regions highlights the repertoire of cancer genes and mutational processes operating, and progresses towards a comprehensive account of the somatic genetic basis of breast cancer.

We have previously identified UMP-CMP kinase (CMPK1) as a prognostic marker for triple negative breast cancer (TNBC) by mass spectrometry (MS). In this study we evaluated CMPK1 association to prognosis in an independent set of samples by immunohistochemistry (IHC) and assessed biological pathways associated to its expression through gene set enrichment analysis (GSEA). A total of 461 TNBC paraffin-embedded tissues were collected from different academic hospitals in Europe, incorporated into tissue micro-arrays (TMA), and stained for CMPK1 expression. We also collected gene expression data of 60 samples, which were also present in the TMA, for GSEA correlation analysis. CMPK1 IHC staining showed both cytoplasmic and nuclear components. While cytoplasmic CMPK1 did not show any association to metastasis free survival (MFS), nuclear CMPK1 was associated to poor prognosis independently from other prognostic factors in stratified Cox regression analyses. GSEA correlation analysis of the nuclear CMPK1-stratified gene expression dataset showed a significant enrichment of extracellular matrix (ECM; positive correlation) and cell cycle (negative correlation) associated genes. We have shown here that nuclear CMPK1 is indicative of poor prognosis in TNBCs and that its expression may be related to dysregulation of ECM and cell cycle molecules.

In breast cancer, GATA3 mutations have been associated with a favorable prognosis and the response to neoadjuvant aromatase inhibitor treatment. Therefore, we investigated whether GATA3 mutations predict the outcome of tamoxifen treatment in the advanced setting. In a retrospective study consisting of 235 hormone-naive patients with ER-positive breast cancer who received tamoxifen as first-line treatment for recurrent disease, GATA3 mutations (in 14.0% of patients) did not significantly associate with either the overall response rate (ORR) or with the length of progression-free survival (PFS) after the start of tamoxifen therapy. Interestingly, among 148 patients for whom both mutation and mRNA expression data were available, GATA3 mutations associated with an increased expression of GATA3. However, only 23.7% of GATA3 high tumors had a mutation. Evaluation of the clinical significance of GATA3 mRNA revealed that it was associated with prolonged PFS, but not with the ORR, also in multivariate analysis. Thus, GATA3 mRNA expression, but not GATA3 mutation, is an independent predictor of prolonged PFS in ER-positive breast cancer patients who received first-line tamoxifen for recurrent disease. Besides GATA3 mutation, other mechanisms must exist that underlie increased GATA3 levels.

Background: CTCs are a promising alternative for metastatic tissue biopsies for use in precision medicine approaches. We investigated to what extent the molecular characteristics of circulat

Background: Studies in taxane and/or anthracycline refractory metastatic breast cancer (mBC) patients have shown approximately 30% response rates to irinotecan. Hence, a significant number of patients will experience irinotecan-induced side effects without obtaining any benefit. The aim of this study was to lay the groundwork for development of predictive biomarkers for irinotecan treatment in BC. Methods: We established BC cell lines with acquired or de novo resistance to SN-38, by exposing the human BC cell lines MCF-7 and MDA-MB-231 to either stepwise increasing concentrations over 6months or an initial high dose of SN-38 (the active metabolite of irinotecan), respectively. The resistant cell lines were analyzed for cross-resistance to other anti-cancer drugs, global gene expression, growth rates, TOP1 and TOP2A gene copy numbers and protein expression, and inhibition of the breast cancer resistance protein (ABCG2/BCRP) drug efflux pump. Results: We found that the resistant cell lines showed 7-100 fold increased resistance to SN-38 but remained sensitive to docetaxel and the non-camptothecin Top1 inhibitor LMP400. The resistant cell lines were characterized by Top1 down-regulation, changed isoelectric points of Top1 and reduced growth rates. The gene and protein expression of ABCG2/BCRP was up-regulated in the resistant sub-lines and functional assays revealed BCRP as a key mediator of SN-38 resistance. Conclusions: Based on our preclinical results, we suggest analyzing the predictive value of the BCRP in breast cancer patients scheduled for irinotecan treatment. Moreover, LMP400 should be tested in a clinical setting in breast cancer patients with resistance to irinotecan.

Tamoxifen therapy resistance constitutes a major cause of death in patients with recurrent estrogen receptor (ER) positive breast cancer. Through high resolution mass spectrometry (MS), we previously generated a 4-protein predictive signature for tamoxifen therapy outcome in recurrent breast cancer. ANXA1 and CALD1, which were not included in the classifier, were however the most differentially expressed proteins. We first evaluated the clinical relevance of these markers in our MS cohort, followed by immunohistochemical (IHC) staining on an independent set of tumors incorporated in a tissue microarray (TMA) and regression analysis in relation to time to progression (TTP), clinical benefit and objective response. In order to assess which mechanisms ANXA1 and CALD1 might been involved in, we performed Ingenuity pathway analysis (IPA) on ANXA1 and CALD1 correlated proteins in our MS cohort. ANXA1 (Hazard ratio [HR] = 1.83; 95% confidence interval [CI]: 1.22-2.75; P = 0.003) and CALD1 (HR = 1.57; 95% CI: 1.04-2.36; P = 0.039) based patient stratification showed significant association to TTP, while IHC staining on TMA showed that both ANXA1 (HR = 1.82; 95% CI: 1.12-3.00; P = 0.016) and CALD1 (HR = 2.29; 95% CI: 1.40-3.75; P = 0.001) expression was associated with shorter TTP independently of traditional predictive fac

INTRODUCTION: Recent advances in molecular diagnostics allow diffuse gliomas to be classified based on their genetic changes into distinct prognostic subtypes. However, a systematic analysis of all molecular markers has thus far not been performed; most classification schemes use a predefined and select set of genes/molecular markers. Here, we have analysed the TCGA dataset (combined glioblastoma (GBM) and lower grade glioma (LGG) datasets) to identify all prognostic genetic markers in diffuse gliomas in order to generate a comprehensive classification scheme. RESULTS: Of the molecular markers investigated (all genes mutated at a population frequency >1.7 % and frequent chromosomal imbalances) in the entire glioma dataset, 57 were significantly associated with overall survival. Of these, IDH1 or IDH2 mutations are associated with lowest hazard ratio, which confirms IDH as the most important prognostic marker in diffuse gliomas. Subsequent subgroup analysis largely confirms many of the currently used molecular classification schemes for diffuse gliomas (ATRX or TP53 mutations, 1p19q codeletion). Our analysis also identified PI3-kinase mutations as markers of poor prognosis in IDH-mutated + ATRX/TP53 mutated diffuse gliomas, median survival 3.7 v. 6.3 years (P = 0.02, Hazard rate (HR) 2.93, 95 % confidence interval (CI) 1.16 - 7.38). PI3-kinase mutations were also prognostic in two independent datasets. In our analysis, no additional molecular markers were identified that further refine the molecular classification of diffuse gliomas. Interestingly, these molecular classifiers do not fully explain the variability in survival observed for diffuse glioma patients. We demonstrate that tumor grade remains an important prognostic factor for overall survival in diffuse gliomas, even within molecular glioma subtypes. Tumor grade was correlated with the mutational load (the number of no

BACKGROUND: Acquired uniparental disomy (aUPD) can lead to homozygosity for tumor suppressor genes or oncogenes. Our purpose is to determine the frequency and profile aUPD regions in serous ovarian cancer (SOC) and investigated the association of aUPD with clinical features and patient outcomes.METHODS: We analyzed single nucleotide polymorphism (SNP) array-based genotyping data on 532 SOC specimens from The Cancer Genome Atlas database to identify aUPD regions. Cox univariate regression and Cox multivariate proportional hazards analyses were performed for survival analysis.RESULTS: We found that 94.7% of SOC samples harbored aUPD; the most common aUPD regions were in chromosomes 17q (76.7%), 17p (39.7%), and 13q (38.3%). In Cox univariate regression analysis, two independent regions of aUPD on chromosome 17q (A and C), and whole-chromosome aUPD were associated with shorter overall survival (OS), and five regions on chromosome 17q (A, D-G) and BRCA1 were associated with recurrence-free survival time. In Cox multivariable proportional hazards analysis, whole-chromosome aUPD was associated with shorter OS. One region of aUPD on chromosome 22q (B) was associated with unilateral disease. A statistically significant association was found between aUPD at TP53 loci and homozygous mutation of TP53 (p < 0.0001).CONCLUSIONS: aUPD is a common event and some recurrent loci are associated with a poor outcome for patients with serous ovarian cancer.

Background: Acquired uniparental disomy (aUPD) can lead to homozygosity for tumor suppressor genes or oncogenes. Our purpose is to determine the frequency and profile aUPD regions in serous ovarian cancer (SOC) and investigated the association of aUPD with clinical features and patient outcomes. Methods: We analyzed single nucleotide polymorphism (SNP) array-based genotyping data on 532 SOC specimens from The Cancer Genome Atlas database to identify aUPD regions. Cox univariate regression and Cox multivariate proportional hazards analyses were performed for survival analysis. Results: We found that 94.7% of SOC samples harbored aUPD; the most common aUPD regions were in chromosomes 17q (76.7%), 17p (39.7%), and 13q (38.3%). In Cox univariate regression analysis, two independent regions of aUPD on chromosome 17q (A and C), and whole-chromosome aUPD were associated with shorter overall survival (OS), and five regions on chromosome 17q (A, D-G) and BRCA1 were associated with recurrence-free survival time. In Cox multivariable proportional hazards analysis, whole-chromosome aUPD was associated with shorter OS. One region of aUPD on chromosome 22q (B) was associated with unilateral disease. A statistically significant association was found between aUPD at TP53 loci and homozygous mutation of TP53 (p < 0.0001). Conclusions: aUPD is a common event and some recurrent loci are associated with a poor outcome for patients with serous ovarian cancer.

The molecular characterization of circulating tumor cells (CTCs) is a promising tool for the repeated and non-invasive evaluation of predictive and prognostic factors. Challenges associated with CTC characterization using the only FDA approved method for CTC enumeration, the CellSearch technique, include the presence of an excess of leukocytes in CTC-enriched blood fractions. Here we aimed to identify colorectal tumor-specific gene expression levels in the blood of patients with and without detectable CTCs according to CellSearch criteria. Materials and methods: Blood of 30 healthy donors (HDs) and 142 metastatic colorectal cancer (mCRC) patients was subjected to CellSearch CTC enumeration and isolation. In all samples, 95 mRNAs were measured by reverse transcriptase quantitative PCR (RT-qPCR). HD blood samples and patient samples with three or more CTCs were compared to identify CTC-specific mRNAs. Patient samples without detectable CTCs were separately analyzed. Results: Thirty-four CTC-specific mRNAs were higher expressed in patients with ≥3 CTCs compared with HDs (Mann-Whitney U-test P<0.05). Among patients without detectable CTCs, a HD-unlike subgroup was identified which could be distinguished from HDs by the expression of epithelial genes such as KRT19, KRT20 and AGR2. Also, in an independent patient set, a similar HD-unlike group could be identified among the patients without detectable CTCs according to the CellSearch system. Conclusion: Extensive molecular characterization of colorectal CTCs is feasible and a subgroup of patients without detectable CTCs according to CellSearch criteria bears circulating tumor load, which may have clinical consequences. This CTC-specific gene panel for mCRC patients may enable the exploration of CTC characterization as a novel means to further individualize cancer treatment.

We here describe two proteomic datasets deposited in ProteomeXchange via PRIDE partner repository [1] with dataset identifiers PXD000484 (defined as "training") and PXD000485 (defined as "test") that have been used for the development of a tamoxifen outcome predictive signature [2]. Both datasets comprised 56 fresh frozen estrogen receptor (ER) positive primary breast tumor specimens derived from patients who received tamoxifen as first line therapy for recurrent disease. Patient groups were defined based on time to progression (TTP) after start of tamoxifen therapy (6 months cutoff): 32 good and 24 poor treatment outcome patients were comprised in the training set, respectively. The test set included 41 good and 15 poor treatment outcome patients. All specimens were subjected to laser capture microdissection (LCM) to enrich for epithelial tumor cells prior to high resolution mass spectrometric (MS) analysis. Protein identificati

Background: Gene expression profiling of tumors is a successful tool for the discovery of new cancer biomarkers and potential targets for the development of new therapeutic strategies. Reliable profiling is preferably performed on fresh frozen (FF) tissues in which the quality of nucleic acids is better preserved than in formalin-fixed paraffin-embedded (FFPE) material. However, since snap-freezing of biopsy materials is often not part of daily routine in pathology laboratories, one may have to rely on archival FFPE material. Procedures to retrieve the RNAs from FFPE materials have been developed and therefore, datasets obtained from FFPE and FF materials need to be made compatible to ensure reliable comparisons are possible. Aim: To develop an efficient method to compare gene expression profiles obtained from FFPE and FF samples using the same platform. Methods: Twenty-six FFPE-FF sample pairs of the same tumors representing various cancer types, and two FFPE-FF sample pairs of breast cancer cell lines, were included. Total RNA was extracted and gene expression profiling was carried out using Illumina's Whole-Genome cDNA-mediated Annealing, Selection, extension and Ligation (WG-DASL) V3 arrays, enabling the simultaneous detection of 24,526 mRNA transcripts. A sample exclusion criterion was created based on the expression of 11 stably expressed reference genes. Pearson correlation at the probe level was calculated for paired FFPE-FF, and three cut-off values were chosen. Spearman correlation coefficients between the matched FFPE and FF samples were calculated for three probe lists with varying levels of significance and compared to the correlation based on all measured probes. Unsupervised hierarchical cluster analysis was performed to verify performance of the included probe lists to compare matched FPPE-FF samples. Results: Twenty-seven FFPE-FF pairs passed the sample exclusion criterion. From the profiles of 27 FFPE and FF matched sampl

BRCA1-mutated breast carcinomas may have distinct biological features, suggesting the involvement of specific oncogenic pathways in tumor development. The identification of genomic aberrations characteristic for BRCA1-mutated breast carcinomas could lead to a better understanding of BRCA1-associated oncogenic events and could prove valuable in clinical testing for BRCA1-involvement in patients. Methods: For this purpose, genomic and gene expression profiles of basal-like BRCA1-mutated breast tumors (n=27) were compared with basal-like familial BRCAX (non-. BRCA1/. 2/. CHEK2*1100delC) tumors (n=14) in a familial cohort of 120 breast carcinomas. Results: Genome wide copy number profiles of the BRCA1-mutated breast carcinomas in our data appeared heterogeneous. Gene expression analyses identifie

Background: Over-treatment of estrogen receptor positive (ER+), lymph node-negative (LNN) breast cancer patients with chemotherapy is a pressing clinical problem that can be addressed by improving techniques to predict tumor metastatic potential. Here we demonstrate that analysis of second harmonic generation (SHG) emission direction in primary tumor biopsies can provide prognostic information about the metastatic outcome of ER+, LNN breast cancer, as well as stage 1 colorectal adenocarcinoma. Methods: SHG is an optical signal produced by fibrillar collagen. The ratio of the forward-to-backward emitted SHG signals (F/B) is sensitive to changes in structure of individual collagen fibers. F/B from excised primary tumor tissue was measured in a retrospective study of LNN breast cancer patients who had received no adjuvant systemic therapy and related to metastasis-free survival (MFS) and overall survival (OS) rates. In addition, F/B was studied for its ass

Epithelial ovarian cancer is a highly heterogeneous disease and remains the most lethal gynaecological malignancy in the Western world. Therapeutic approaches need to account for inter-patient and intra-tumoural heterogeneity and detailed characterization of in vitro models representing the different histological an

Background: Clinical outcome of patients with triple-negative breast cancer (TNBC) is highly variable. This study aims to identify and validate a prognostic protein signature for TNBC patients to reduce unnecessary adjuv

Introduction: Breast cancer is a genetically and phenotypically complex disease. To understand the role of miRNAs in this molecular complexity, we performed miRNA expression analysis in a cohort of molecularly well-characterized human breast cancer cell lines to identify miRNAs associated with the most common molecular subtypes and the most frequent genetic aberrations. Methods: Using a microarray carrying LNA™ modified oligonucleotide capture probes), expression levels of 725 human miRNAs were measured in 51 breast cancer cell lines. Differential miRNA expression was explored by unsupervised cluster analysis and was then associated with the molecular subtypes and genetic aberrations commonly present in breast cancer. Results: Unsupervised cluster analysis using the most variably expressed miRNAs divided the 51 breast cancer cell lines into a major and a minor cluster predominantly mirroring the luminal and basal intrinsic subdivision of breast cancer cell lines. One hundred and thirteen miRNAs were differentially expressed between these two main clusters. Forty miRNAs were differentially expressed between basal-like and normal-like/claudin-low cell lines. Within the luminal-group, 39 miRNAs were associated with ERBB2 overexpression and 24 with E-cadherin gene mutations, which are frequent in this subtype of breast cancer cell lines. In contrast, 31 miRNAs were associated with E-cadherin promoter hypermethylation, which, contrary to E-cadherin mutation, is exclusively observed in breast cancer cell lines that are not of luminal origin. Thirty miRNAs were associated with p16INK4 status while only a few

Purpose: Previously, we developed a radiosensitivity molecular signature [radiosensitivity index (RSI)] that was clinically validated in 3 independent datasets (rectal, esophageal, and head and neck) in 118 patients. Here, we test RSI in radiotherapy (RT)-treated breast cancer patients. Experimental Design: RSI was tested in 2 previously published breast cancer datasets. Patients were treated at the Karolinska University Hospital (n = 159) and Erasmus Medical Center (n = 344). RSI was applied as previously described. Results: We tested RSI in RT-treated patients (Karolinska). Patients predicted to be radiosensitive (RS) had an improved 5-year relapse-free survival when compared with radioresistant (RR) patients (95% vs. 75%, P = 0.0212), but there was no difference between RS/RR patients treated without RT (71% vs. 77%, P = 0.6744), consistent with RSI being RT-specific (interaction term RSI - RT, P = 0.05). Similarly, in the Erasmus dataset, RT-treated RS patients had an improved 5-year distant metastasis-free survival over RR patients (77% vs. 64%, P = 0.0409), but no difference was observed in patients treated without RT (RS vs. RR, 80% vs. 81%, P = 0.9425). Multivariable analysis showed RSI is the strongest variable in RT-treated patients (Karolinska, HR

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 associ

Background: Genetic alterations in cellular signaling networks are a hallmark of cancer, however, effective methods to discover them are lacking. A novel form of abnormality called acquired uniparental disomy (aUPD) was recently found to pinpoint the region of mutated genes in various cancers, thereby identifying the region for next-generation sequencing. Methods/Principal Findings: We retrieved large genomic data sets from the Gene Expression Omnibus database to perform genome-wide analysis of aUPD in breast tumor samples and cell lines using approaches that can reliably detect aUPD. Aupd was identified in 52.29% of the tumor samples. The most frequent aUPD regions were located at chromosomes 2q, 3p, 5q, 9p, 9q, 10q, 11q, 13q, 14q and 17q. We evaluated the data for any correlation between the most frequent aUPD regions and HER2/neu, ER, and PR status, and found a statistically significant correlation between the recurrent regions of aUPD and triple negative (TN) breast cancers. aUPD at chromosome 17q (VEZF1, WNT3), 3p (SUMF1, GRM7), 9p (MTAP, NFIB) and 11q (CASP1, CASP4, CASP5) are predictors for TN. The frequency of aUPD was found to be significantly higher in TN breast cancer cases compared to HER2/neu-positive and/or ER or PR-positive cases. Furthermore, using previously published mutation data, we found TP53 homozygously mutated in cell lines having aUPD in that locus. Conclusions/Significance: We conclude that aUPD is a common and non-random molecular feature of breast cancer that is most prominent in triple negative cases. As aUPD regions are different among the main pathological subtypes, specific aUPD regions may aid the sub-classification of breast cancer. In addition, we provide statistical support using TP53 as an example that identifying aUPD regions can be an effective approach in finding aberrant genes. We thus conclud