BackgroundColorectal cancer (CRC) cell lines are widely used pre-clinical model systems. Comprehensive insights into their molecular characteristics may improve model selection for biomedical studies.MethodsWe have performed DNA, RNA and protein profiling of 34 cell lines, including (i) targeted deep sequencing (n = 612 genes) to detect single nucleotide variants and insertions/deletions; (ii) high resolution DNA copy number profiling; (iii) gene expression profiling at exon resolution; (iv) small RNA expression profiling by deep sequencing; and (v) protein expression analysis (n = 297 proteins) by reverse phase protein microarrays.ResultsThe cell lines were stratified according to the key molecular subtypes of CRC and data were integrated at two or more levels by computational analyses. We confirm that the frequencies and patterns of DNA aberrations are associated with genomic instability phenotypes and that the cell lines recapitulate the genomic profiles of primary carcinomas. Intrinsic expression subgroups are distinct from genomic subtypes, but consistent at the gene-, microRNA- and protein-level and dominated by two distinct clusters; colon-like cell lines characterized by expression of gastro-intestinal differentiation markers and undifferentiated cell lines showing upregulation of epithelial-mesenchymal transition and TGFβ signatures. This sample split was concordant with the gene expression-based consensus molecular subtypes of primary tumors. Approximately ¼ of the genes had consistent regulation at the DNA copy number and gene expression level, while expression of gene-protein pairs in general was strongly correlated. Consistent high-level DNA copy number amplification and outlier gene- and protein- expression was found for several oncogenes in individual cell lines, including MYC and ERBB2.ConclusionsThis study expands the view of CRC cell lines as accurate molecular models of primary carcinomas, and we present integrated multi-level molecular data of 34 widely used cell lines in easily accessible formats, providing a resource for preclinical studies in CRC. Electronic supplementary materialThe online version of this article (doi:10.1186/s12943-017-0691-y) contains supplementary material, which is available to authorized users.

No abstract

Patterns of glycosylation are important in cancer, but the molecular mechanisms that drive changes are often poorly understood. The androgen receptor drives prostate cancer (PCa) development and progression to lethal metastatic castration-resistant disease. Here we used RNA-Seq coupled with bioinformatic analyses of androgen-receptor (AR) binding sites and clinical PCa expression array data to identify ST6GalNAc1 as a direct and rapidly activated target gene of the AR in PCa cells. ST6GalNAc1 encodes a sialytransferase that catalyses formation of the cancer-associated sialyl-Tn antigen (sTn), which we find is also induced by androgen exposure. Androgens induce expression of a novel splice variant of the ST6GalNAc1 protein in PCa cells. This splice variant encodes a shorter protein isoform that is still fully functional as a sialyltransferase and able to induce expression of the sTn-antigen. Surprisingly, given its high expression in tumours, stable expression of ST6GalNAc1 in PCa cells reduced formation of stable tumours in mice, reduced cell adhesion and induced a switch towards a more mesenchymal-like cell phenotype in vitro. ST6GalNAc1 has a dynamic expression pattern in clinical datasets, being significantly up-regulated in primary prostate carcinoma but relatively down-regulated in established metastatic tissue. ST6GalNAc1 is frequently upregulated concurrently with another important glycosylation enzyme GCNT1 previously associated with prostate cancer progression and implicated in Sialyl Lewis X antigen synthesis. Together our data establishes an androgen-dependent mechanism for sTn antigen expression in PCa, and are consistent with a general role for the androgen receptor in driving important coordinate changes to the glycoproteome during PCa progression.

BackgroundApproximately 15% of primary colorectal cancers have DNA mismatch repair deficiency, causing a complex genome with thousands of small mutations—the microsatellite instability (MSI) phenotype. We investigated molecular heterogeneity and tumor immunogenicity in relation to clinical endpoints within this distinct subtype of colorectal cancers.MethodsA total of 333 primary MSI+ colorectal tumors from multiple cohorts were analyzed by multilevel genomics and computational modeling—including mutation profiling, clonality modeling, and neoantigen prediction in a subset of the tumors, as well as gene expression profiling for consensus molecular subtypes (CMS) and immune cell infiltration.ResultsNovel, frequent frameshift mutations in four cancer-critical genes were identified by deep exome sequencing, including in CRTC1, BCL9, JAK1, and PTCH1. JAK1 loss-of-function mutations were validated with an overall frequency of 20% in Norwegian and British patients, and mutated tumors had up-regulation of transcriptional signatures associated with resistance to anti-PD-1 treatment. Clonality analyses revealed a high level of intra-tumor heterogeneity; however, this was not associated with disease progression. Among the MSI+ tumors, the total mutation load correlated with the number of predicted neoantigens (P = 4 × 10−5), but not with immune cell infiltration—this was dependent on the CMS class; MSI+ tumors in CMS1 were highly immunogenic compared to MSI+ tumors in CMS2-4. Both JAK1 mutations and CMS1 were favorable prognostic factors (hazard ratios 0.2 [0.05–0.9] and 0.4 [0.2–0.9], respectively, P = 0.03 and 0.02).ConclusionsMultilevel genomic analyses of MSI+ colorectal cancer revealed molecular heterogeneity with clinical relevance, including tumor immunogenicity and a favorable patient outcome associated with JAK1 mutations and the transcriptomic subgroup CMS1, emphasizing the potential for prognostic stratification of this clinically important subtype.See related research highlight by Samstein and Chan 10.1186/s13073-017-0438-9 Electronic supplementary materialThe online version of this article (doi:10.1186/s13073-017-0434-0) contains supplementary material, which is available to authorized users.

Background: In breast cancer, activation of bone morphogenetic protein (BMP) signaling and elevated levels of BMP-antagonists have been linked to tumor progression and metastasis. However, the simultaneous upregulation of BMPs and their antagonist, and the fact that both promote tumor aggressiveness seems contradictory and is not fully understood. Methods: We analyzed the transcriptomes of the metastatic 66cl4 and the non-metastatic 67NR cell lines of the 4T1 mouse mammary tumor model to search for factors that promote metastasis. CRISPR/Cas9 gene editing was used for mechanistic studies in the same cell lines. Furthermore, we analyzed gene expression patterns in human breast cancer biopsies obtained from public datasets to evaluate co-expression and possible relations to clinical outcome. Results: We found that mRNA levels of the BMP-antagonist Grem1, encoding gremlin1, and the ligand Bmp4 were both significantly upregulated in cells and primary tumors of 66cl4 compared to 67NR. Depletion of gremlin1 in 66cl4 could impair metastasis to the lungs in this model. Furthermore, we found that expression of Grem1 correlated with upregulation of several stem cell markers in 66cl4 cells compared to 67NR cells. Both in the mouse model and in patients, expression of GREM1 associated with extracellular matrix organization, and formation, biosynthesis and modification of collagen. Importantly, high expression of GREM1 predicted poor prognosis in estrogen receptor negative breast cancer patients. Analyses of large patient cohorts revealed that amplification of genes encoding BMP-antagonists and elevation of the corresponding transcripts is evident in biopsies from more than half of the patients and much more frequent for the secreted BMP-antagonists than the intracellular inhibitors of SMAD signaling. Conclusion: In conclusion, our results show that GREM1 is associated with metastasis and predicts poor prognosis in ER-negative breast cancer patients. Gremlin1 could represent a novel target for therapy.

Given a graph with nonnegative edge weights and a set D of node pairs, the 2-path network problem requires a minimum weight set of edges such that the induced subgraph contains a path with one or two edges connecting each pair in D. The problem is NP-hard. We present two integer programming models for the problem and study properties of associated polytopes, including cutting planes. Two approximation algorithms are suggested and analyzed. Some computational experience is reported.

We investigated the intra-patient heterogeneity of driver gene mutations among colorectal liver metastases by sequencing 479 tumor samples from 106 patients. A near-perfect intra-patient concordance was found in the mutation status of the primary tumor and multiple metastatic lesions of KRAS/NRAS/BRAF and PIK3CA when high-sensitivity methods were applied. Mutations in KRAS alone and KRAS/NRAS/BRAF combined had a negative prognostic impact after liver resection. Background: The prevalence and clinical implications of genetic heterogeneity in patients with multiple colorectal liver metastases remain largely unknown. In a prospective series of patients undergoing resection of colorectal liver metastases, the aim was to investigate the inter-metastatic and primary-to-metastatic heterogeneity of mutations in KRAS, NRAS, BRAF, and PIK3CA and their prognostic impact. Patients and Methods: We analyzed the mutation status among 372 liver metastases and 78 primary tumors from 106 patients by methods used in clinical routine testing, by Sanger sequencing, by next-generation sequencing (NGS), and/or by droplet digital polymerase chain reaction. The 3-year cancer-specific survival (CSS) was analyzed using the Kaplan-Meier method. Results: Although Sanger sequencing indicated inter-metastatic mutation heterogeneity in 14 of 97 patients (14%), almost all cases were refuted by high-sensitive NGS. Also, heterogeneity among metastatic deposits was concluded only for PIK3CA in 2 patients. Similarly, primary-to-metastatic heterogeneity was indicated in 8 of 78 patients (10%) using Sanger sequencing but for only 2 patients after NGS, showing the emergence of 1 KRAS and 1 PIK3CA mutation in the metastatic lesions. KRAS mutations were present in 53 of 106 patients (50%) and were associated with poorer 3-year CSS after liver resection (37% vs. 61% for KRAS wild-type; P ¼ .004). Poor prognostic associations were found also for the combination of KRAS/NRAS/BRAF mutations compared with triple wild-type (P ¼ .002). Conclusion: Intrapatient mutation heterogeneity was virtually undetected, both between the primary tumor and the liver metastases and among the metastatic deposits. KRAS mutations separately, and KRAS/NRAS/BRAF mutations combined, were associated with poor patient survival after partial liver resection.

BackgroundA fusion gene is a hybrid gene consisting of parts from two previously independent genes. Chromosomal rearrangements leading to gene breakage are frequent in high-grade serous ovarian carcinomas and have been reported as a common mechanism for inactivating tumor suppressor genes. However, no fusion genes have been repeatedly reported to be recurrent driver events in ovarian carcinogenesis. We combined genomic and transcriptomic information to identify novel fusion gene candidates and aberrantly expressed genes in ovarian carcinomas.MethodsExamined were 19 previously karyotyped ovarian carcinomas (18 of the serous histotype and one undifferentiated). First, karyotypic aberrations were compared to fusion gene candidates identified by RNA sequencing (RNA-seq). In addition, we used exon-level gene expression microarrays as a screening tool to identify aberrantly expressed genes possibly involved in gene fusion events, and compared the findings to the RNA-seq data.ResultsWe found a DPP9-PPP6R3 fusion transcript in one tumor showing a matching genomic 11;19-translocation. Another tumor had a rearrangement of DPP9 with PLIN3. Both rearrangements were associated with diminished expression of the 3′ end of DPP9 corresponding to the breakpoints identified by RNA-seq. For the exon-level expression analysis, candidate fusion partner genes were ranked according to deviating expression compared to the median of the sample set. The results were collated with data obtained from the RNA-seq analysis. Several fusion candidates were identified, among them TMEM123-MMP27, ZBTB46-WFDC13, and PLXNB1-PRKAR2A, all of which led to stronger expression of the 3′ genes. In view of our previous findings of nonrandom rearrangements of chromosome 19 in this cancer type, particular emphasis was given to changes of this chromosome and a DDA1-FAM129C fusion event was identified.ConclusionsWe have identified novel fusion gene candidates in high-grade serous ovarian carcinoma. DPP9 was involved in two different fusion transcripts that both resulted in deregulated expression of the 3′ end of the transcript and thus possible loss of the active domains in the DPP9 protein. The identified rearrangements might play a role in tumorigenesis or tumor progression.Electronic supplementary materialThe online version of this article (10.1186/s12885-017-3625-6) contains supplementary material, which is available to authorized users.