Multi-omics of 34 colorectal cancer cell lines - a resource for biomedical studies.

BACKGROUND: Colorectal cancer (CRC) cell lines are widely used pre-clinical model systems. Comprehensive insights into their molecular characteristics may improve model selection for biomedical studies. METHODS: We 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. RESULTS: The 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. CONCLUSIONS: This 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.

Transcriptomic subtyping of malignant peripheral nerve sheath tumours highlights immune signatures, genomic profiles, patient survival and therapeutic targets.

BACKGROUND: Malignant peripheral nerve sheath tumour (MPNST) is an aggressive orphan disease commonly affecting adolescents or young adults. Current knowledge of molecular tumour biology has been insufficient for development of rational treatment strategies. We aimed to discover molecular subtypes of potential clinical relevance. METHODS: Fresh frozen samples of MPNSTs (n = 94) and benign neurofibromas (n = 28) from 115 patients in a European multicentre study were analysed by DNA copy number and/or transcriptomic profiling. Unsupervised transcriptomic subtyping was performed and the subtypes characterized for genomic aberrations, clinicopathological associations and patient survival. FINDINGS: MPNSTs were classified into two transcriptomic subtypes defined primarily by immune signatures and proliferative processes. "Immune active" MPNSTs (44%) had sustained immune signals relative to neurofibromas, were more frequently low-grade (P = 0.01) and had favourable prognostic associations in a multivariable model of disease-specific survival with clinicopathological factors (hazard ratio 0.25, P = 0.003). "Immune deficient" MPNSTs were more aggressive and characterized by proliferative signatures, high genomic complexity, aberrant TP53 and PRC2 loss, as well as high relative expression of several potential actionable targets (EGFR, ERBB2, EZH2, KIF11, PLK1, RRM2). Integrated gene-wise analyses suggested a DNA copy number-basis for proliferative transcriptomic signatures in particular, and the tumour copy number burden further stratified the transcriptomic subtypes according to patient prognosis (P < 0.01). INTERPRETATION: Approximately half of MPNSTs belong to an "immune deficient" transcriptomic subtype associated with an aggressive disease course, PRC2 loss and expression of several potential therapeutic targets, providing a rationale for molecularly-guided intervention trials. FUNDING: Research grants from non-profit organizations, as stated in the Acknowledgements.

Metastatic heterogeneity of the consensus molecular subtypes of colorectal cancer.

Gene expression-based subtypes of colorectal cancer have clinical relevance, but the representativeness of primary tumors and the consensus molecular subtypes (CMS) for metastatic cancers is not well known. We investigated the metastatic heterogeneity of CMS. The best approach to subtype translation was delineated by comparisons of transcriptomic profiles from 317 primary tumors and 295 liver metastases, including multi-metastatic samples from 45 patients and 14 primary-metastasis sets. Associations were validated in an external data set (n = 618). Projection of metastases onto principal components of primary tumors showed that metastases were depleted of CMS1-immune/CMS3-metabolic signals, enriched for CMS4-mesenchymal/stromal signals, and heavily influenced by the microenvironment. The tailored CMS classifier (available in an updated version of the R package CMScaller) therefore implemented an approach to regress out the liver tissue background. The majority of classified metastases were either CMS2 or CMS4. Nonetheless, subtype switching and inter-metastatic CMS heterogeneity were frequent and increased with sampling intensity. Poor-prognostic value of CMS1/3 metastases was consistent in the context of intra-patient tumor heterogeneity.

Genomic and prognostic heterogeneity among RAS/BRAFV600E /TP53 co-mutated resectable colorectal liver metastases.

Hepatic resection is potentially curative for patients with colorectal liver metastases, but the treatment benefit varies. KRAS/NRAS (RAS)/TP53 co-mutations are associated with a poor prognosis after resection, but there is large variation in patient outcome within the mutation groups, and genetic testing is currently not used to evaluate benefit from surgery. We have investigated the potential for improved prognostic stratification by combined biomarker analysis with DNA copy number aberrations (CNAs), and taking tumor heterogeneity into account. We determined the mutation status of RAS, BRAFV600 , and TP53 in 441 liver lesions from 171 patients treated by partial hepatectomy for metastatic colorectal cancer. CNAs were profiled in 232 tumors from 67 of the patients. Mutations and high-level amplifications of cancer-critical genes, the latter including ERBB2 and EGFR, were predominantly homogeneous within patients. RAS/BRAFV600E and TP53 co-mutations were associated with a poor patient outcome (hazard ratio, HR, 3.9, 95% confidence interval, CI, 1.3-11.1, P = 0.012) in multivariable analyses with clinicopathological variables. The genome-wide CNA burden and intrapatient intermetastatic CNA heterogeneity varied within the mutation groups, and the CNA burden had prognostic associations in univariable analysis. Combined prognostic analyses of RAS/BRAFV600E /TP53 mutations and CNAs, either as a high CNA burden or high intermetastatic CNA heterogeneity, identified patients with a particularly poor outcome (co-mutation/high CNA burden: HR 2.7, 95% CI 1.2-5.9, P = 0.013; co-mutation/high CNA heterogeneity: HR 2.5, 95% CI 1.1-5.6, P = 0.022). In conclusion, DNA copy number profiling identified genomic and prognostic heterogeneity among patients with resectable colorectal liver metastases with co-mutated RAS/BRAFV600E /TP53.

De novo transcriptomic subtyping of colorectal cancer liver metastases in the context of tumor heterogeneity.

BACKGROUND: Gene expression-based subtyping has the potential to form a new paradigm for stratified treatment of colorectal cancer. However, current frameworks are based on the transcriptomic profiles of primary tumors, and metastatic heterogeneity is a challenge. Here we aimed to develop a de novo metastasis-oriented framework. METHODS: In total, 829 transcriptomic profiles from patients with colorectal cancer were analyzed, including primary tumors, liver metastases, and non-malignant liver samples. High-resolution microarray gene expression profiling was performed of 283 liver metastases from 171 patients treated by hepatic resection, including multiregional and/or multi-metastatic samples from each of 47 patients. A single randomly selected liver metastasis sample from each patient was used for unsupervised subtype discovery by nonnegative matrix factorization, and a random forest prediction model was trained to classify multi-metastatic samples, as well as liver metastases from two independent series of 308 additional patients. RESULTS: Initial comparisons with non-malignant liver samples and primary colorectal tumors showed a highly variable degree of influence from the liver microenvironment in metastases, which contributed to inter-metastatic transcriptomic heterogeneity, but did not define subtype distinctions. The de novo liver metastasis subtype (LMS) framework recapitulated the main distinction between epithelial-like and mesenchymal-like tumors, with a strong immune and stromal component only in the latter. We also identified biologically distinct epithelial-like subtypes originating from different progenitor cell types. LMS1 metastases had several transcriptomic features of cancer aggressiveness, including secretory progenitor cell origin, oncogenic addictions, and microsatellite instability in a microsatellite stable background, as well as frequent RAS/TP53 co-mutations. The poor-prognostic association of LMS1 metastases was independent of mutation status, clinicopathological variables, and current subtyping frameworks (consensus molecular subtypes and colorectal cancer intrinsic subtypes). LMS1 was also the least heterogeneous subtype in comparisons of multiple metastases per patient, and tumor heterogeneity did not confound the prognostic value of LMS1. CONCLUSIONS: We report the first large study of multi-metastatic gene expression profiling of colorectal cancer. The new metastasis-oriented subtyping framework showed potential for clinically relevant transcriptomic classification in the context of metastatic heterogeneity, and an LMS1 mini-classifier was constructed to facilitate prognostic stratification and further clinical testing.

Molecular correlates of sensitivity to PARP inhibition beyond homologous recombination deficiency in pre-clinical models of colorectal cancer point to wild-type TP53 activity.

BACKGROUND: PARP inhibitors are active in various tumour types beyond BRCA-mutant cancers, but their activity and molecular correlates in colorectal cancer (CRC) are not well studied. METHODS: Mutations and genome-wide mutational patterns associated with homologous recombination deficiency (HRD) were investigated in 255 primary CRCs with whole-exome sequencing and/or DNA copy number data. Efficacy of five PARP inhibitors and their molecular correlates were evaluated in 93 CRC cell lines partly annotated with mutational-, DNA copy number-, and/or gene expression profiles. Post-treatment gene expression profiling and specific protein expression analyses were performed in two pairs of PARP inhibitor sensitive and resistant cell lines. FINDINGS: A subset of microsatellite stable (MSS) CRCs had truncating mutations in homologous recombination-related genes, but these were not associated with genomic signatures of HRD. Eight CRC cell lines (9%) were sensitive to PARP inhibition, but sensitivity was not predicted by HRD-related genomic and transcriptomic signatures. In contrast, drug sensitivity in MSS cell lines was strongly associated with TP53 wild-type status (odds ratio 15.7, p = 0.023) and TP53-related expression signatures. Increased downstream TP53 activity was among the primary response mechanisms, and TP53 inhibition antagonized the effect of PARP inhibitors. Wild-type TP53-mediated suppression of RAD51 was identified as a possible mechanism of action for sensitivity to PARP inhibition. INTERPRETATION: PARP inhibitors are active in a subset of CRC cell lines and preserved TP53 function may increase the likelihood of response.

Frequent copy number gains of SLC2A3 and ETV1 in testicular embryonal carcinomas.

Testicular germ cell tumours (TGCTs) appear as different histological subtypes or mixtures of these. They show similar, multiple DNA copy number changes, where gain of 12p is pathognomonic. However, few high-resolution analyses have been performed and focal DNA copy number changes with corresponding candidate target genes remain poorly described for individual subtypes. We present the first high-resolution DNA copy number aberration (CNA) analysis on the subtype embryonal carcinomas (ECs), including 13 primary ECs and 5 EC cell lines. We identified recurrent gains and losses and allele-specific CNAs. Within these regions, we nominate 30 genes that may be of interest to the EC subtype. By in silico analysis of data from 150 TGCTs from The Cancer Genome Atlas (TCGA), we further investigated CNAs, RNA expression, somatic mutations and fusion transcripts of these genes. Among primary ECs, ploidy ranged between 2.3 and 5.0, and the most common aberrations were DNA copy number gains at chromosome (arm) 7, 8, 12p, and 17, losses at 4, 10, 11, and 18, replicating known TGCT genome characteristics. Gain of whole or parts of 12p was found in all samples, including a highly amplified 100 kbp segment at 12p13.31, containing SLC2A3. Gain at 7p21, encompassing ETV1, was the second most frequent aberration. In conclusion, we present novel CNAs and the genes located within these regions, where the copy number gain of SLC2A3 and ETV1 are of interest, and which copy number levels also correlate with expression in TGCTs.

High Concordance and Negative Prognostic Impact of RAS/BRAF/PIK3CA Mutations in Multiple Resected Colorectal Liver Metastases.

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: Intra-patient 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.

Gene expression profiles of CMS2-epithelial/canonical colorectal cancers are largely driven by DNA copy number gains.

About 80% of colorectal cancers (CRCs) have chromosomal instability, which is an integral part of aggressive malignancy development, but the importance of specific copy number aberrations (CNAs) in modulating gene expression, particularly within the framework of clinically relevant molecular subtypes, remains mostly elusive. We performed DNA copy number profiling of 257 stage I-IV primary CRCs and integrative gene expression analysis in 151 microsatellite stable (MSS) tumors, focusing on high-level amplifications and the effect of CNAs on the characteristics of the gene expression-based consensus molecular subtypes (CMS). The results were validated in 323 MSS tumors from TCGA. Novel recurrent high-level amplifications (≥15 additional copies) with a major impact on gene expression were found for TOX3 (16q) at 1.5% frequency, as well as for CCND2 (12p) and ANXA11 (10q) at 1% frequency, in addition to the well-known targets ERBB2 (17q) and MYC (8q). Focal amplifications with ≥15 or ≥5 additional copies of at least one of these regions were associated with a poor overall survival among patients with stage I-III MSS CRCs (multivariable hazard ratio ≥3.2, p ≤ 0.01). All high-level amplifications were focal and had a more consistent relationship with gene expression than lower amplitude and/or broad-range amplifications, suggesting specific targeting during carcinogenesis. Genome-wide, copy number driven gene expression was enriched for pathways characteristic of the CMS2-epithelial/canonical subtype, including DNA repair and cell cycle progression. Furthermore, 50% of upregulated genes in CMS2-epithelial/canonical MSS CRCs were driven by CNAs, an enrichment compared with the other CMS groups, and associated with the stronger correspondence between CNAs and gene expression in malignant epithelial cells than in the cells of the tumor microenvironment (fibroblasts, endothelial cells, leukocytes). In conclusion, we identify novel recurrent amplifications with impact on gene expression in CRC and provide the first evidence that CMS2 may have a stronger copy-number related genetic basis than subtypes more heavily influenced by gene expression signals from the tumor microenvironment.

A robust internal control for high-precision DNA methylation analyses by droplet digital PCR.

Background: Droplet digital PCR (ddPCR) allows absolute quantification of nucleic acids and has potential for improved non-invasive detection of DNA methylation. For increased precision of the methylation analysis, we aimed to develop a robust internal control for use in methylation-specific ddPCR. Methods: Two control design approaches were tested: (a) targeting a genomic region shared across members of a gene family and (b) combining multiple assays targeting different pericentromeric loci on different chromosomes. Through analyses of 34 colorectal cancer cell lines, the performance of the control assay candidates was optimized and evaluated, both individually and in various combinations, using the QX200™ droplet digital PCR platform (Bio-Rad). The best-performing control was tested in combination with assays targeting methylated CDO1, SEPT9, and VIM. Results: A 4Plex panel consisting of EPHA3, KBTBD4, PLEKHF1, and SYT10 was identified as the best-performing control. The use of the 4Plex for normalization reduced the variability in methylation values, corrected for differences in template amount, and diminished the effect of chromosomal aberrations. Positive Droplet Calling (PoDCall), an R-based algorithm for standardized threshold determination, was developed, ensuring consistency of the ddPCR results. Conclusion: Implementation of a robust internal control, i.e., the 4Plex, and an algorithm for automated threshold determination, PoDCall, in methylation-specific ddPCR increase the precision of DNA methylation analysis.

Multilevel genomics of colorectal cancers with microsatellite instability-clinical impact of JAK1 mutations and consensus molecular subtype 1.

BACKGROUND: Approximately 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. METHODS: A 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. RESULTS: Novel, 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). CONCLUSIONS: Multilevel 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.

Quantification of the concentration gradient of biomarkers between ovarian carcinoma interstitial fluid and blood.

BACKGROUND: Tumor interstitial fluid (TIF) rather than plasma should be used in cancer biomarker discovery because of the anticipated higher concentration of locally produced proteins in the tumor microenvironment. Nevertheless, the actual TIF-to-plasma gradient of tumor specific proteins has not been quantified. We present the proof-of-concept for the quantification of the postulated gradient between TIF and plasma. METHODS: TIF was collected by centrifugation from serous (n = 19), endometrioid (n = 9) and clear cell (n = 3) ovarian carcinomas with early (n = 15) and late stage (n = 16) disease in grades 1 (n = 2), 2 (n = 8) and 3 (n = 17), and ELISA was used for the determination of CA-125, osteopontin and VEGF-A. RESULTS: All three markers were significantly up-regulated in TIF compared with plasma (p < 0.0001). The TIF-to-plasma ratio of the ovarian cancer biomarker CA-125 ranged from 1.4 to 24,300 (median = 194) and was inversely correlated to stage (p = 0.0006). The cancer related osteopontin and VEGF-A had TIF-to-plasma ratios ranging from 1 to 62 (median = 15) and 2 to 1040 (median = 59), respectively. The ratios were not affected by tumor stage, indicative of more widespread protein expression. CONCLUSION: We present absolute quantitative data on the TIF-to-plasma gradient of selected proteins in the tumor microenvironment, and demonstrate a substantial and stage dependent gradient for CA-125 between TIF and plasma, suggesting a relation between total tumor burden and tissue-to-plasma gradient. GENERAL SIGNIFICANCE: We present novel quantitative data on biomarker concentration in the tumor microenvironment, and a new strategy for biomarker selection, applicable in future biomarker studies.