Ina A. Eilertsen

Epigenetic and genetic features of 24 colon cancer cell lines

Cell lines are invaluable biomedical research tools, and recent literature has emphasized the importance of genotype authentication and characterization. In the present study, 24 out of 27 cell line identities were confirmed by short tandem repeat profiling. The molecular phenotypes of the 24 colon cancer cell lines were examined, and microsatellite instability (MSI) and CpG island methylator phenotype (CIMP) were determined, using the Bethesda panel mononucleotide repeat loci and two epimarker panels, respectively. Furthermore, the BRAF, KRAS and PIK3CA oncogenes were analyzed for mutations in known hotspots, while the entire coding sequences of the PTEN and TP53 tumor suppressors were investigated. Nine cell lines showed MSI. Thirteen and nine cell lines were found to be CIMP positive, using the Issa panel and the Weisenberger et al. panel, respectively. The latter was found to be superior for CIMP classification of colon cancer cell lines. Seventeen cell lines harbored disrupting TP53 mutations. Altogether, 20/24 cell lines had the mitogen-activated protein kinase pathway activating mutually exclusive KRAS or BRAF mutations. PIK3CA and PTEN mutations leading to hyperactivation of the phosphoinositide 3-kinase/AKT pathway were observed in 13/24 cell lines. Interestingly, in four cell lines there were no mutations in neither BRAF, KRAS, PIK3CA nor in PTEN. In conclusion, this study presents molecular features of a large number of colon cancer cell lines to aid the selection of suitable in vitro models for descriptive and functional research.

NEDD4 is overexpressed in colorectal cancer and promotes colonic cell growth independently of the PI3K/PTEN/AKT pathway

Ubiquitination controls multiple cellular processes relevant to cancer pathogenesis. Using Gene Set Enrichment Analysis of an mRNA transcriptome dataset, we have identified genes encoding components of the ubiquitin system that are differentially expressed in colorectal cancers as compared to normal colonic mucosa. Among the significantly overexpressed genes was NEDD4 (neural precursor cell-expressed developmentally down-regulated 4), the prototype member of the HECT (homologous to E6AP C-terminus) E3 ubiquitin ligase family. Previous studies have shown that NEDD4 may act as an oncoprotein by inducing ubiquitination and degradation of the tumor suppressor protein PTEN (phosphatase and tensin homolog). To investigate its functional importance in colorectal cancer, HCT-15 and LoVo colon cancer cells were depleted of NEDD4 by small interfering RNA. The depletion resulted in reduced growth and altered cell morphology in both cell lines. However, NEDD4 depletion did not affect the PTEN protein level or PI3K/AKT signaling pathway activation. Moreover, ectopic expression of NEDD4 did not influence the PTEN subcellular localization or protein level. Collectively, these data demonstrate that NEDD4 is overexpressed in colorectal cancers, and suggest that NEDD4 promotes growth of colon cancer cells independently of PTEN and PI3K/AKT signaling.

the novel colorectal cancer biomarkers cdo1, Zscan18 and znf331 are frequently methylated across gastrointestinal cancers

We have previously shown that gastrointestinal cancers display similar epigenetic aberrations. In a recent study, we identified frequently methylated genes for cholangiocarcinoma (CDO1, DCLK1, SFRP1 and ZSCAN18), where one of these genes, DCLK1, was also confirmed to be highly methylated in colorectal cancer. The aim of the present study was to determine whether these four genes, in addition to one gene found to be methylated in colon cancer cell lines (ZNF331), are commonly methylated across gastrointestinal malignancies, as well as explore their role as potential biomarkers. Quantitative methylation specific PCR (qMSP) of colorectal cancer (n = 164) and normal colorectal mucosa (n = 106) samples showed that all genes were frequently methylated in colorectal cancer (71–92%) with little or no methylation in normal mucosa (0–3%). Methylation of minimum two of these five genes identified 95% of the tumors with a specificity of 98%, and an area under the receiver operating characteristics curve (AUC) of 0.98. For gastric (n = 25) and pancreatic (n = 20) cancer, the same panel detected 92% and 90% of the tumors, respectively. Seventy‐four cancer cell lines were further analyzed by qMSP and real time RT‐PCR. In addition to the previously reported DCLK1, a high negative correlation between promoter DNA methylation and gene expression was observed for CDO1, ZNF331 and ZSCAN18. In conclusion, the high methylation frequency of these genes in colorectal‐ as well as in gastric‐, pancreatic‐ and bile duct cancer confirmed an epigenetic similarity between gastrointestinal cancer types, and simultaneously demonstrated their potential as biomarkers, particularly for colorectal cancer detection.

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

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

Colorectal Cancer Consensus Molecular Subtypes Translated to Preclinical Models Uncover Potentially Targetable Cancer Cell Dependencies

Purpose: Response to standard oncologic treatment is limited in colorectal cancer. The gene expression–based consensus molecular subtypes (CMS) provide a new paradigm for stratified treatment and drug repurposing; however, drug discovery is currently limited by the lack of translation of CMS to preclinical models. Experimental Design: We analyzed CMS in primary colorectal cancers, cell lines, and patient-derived xenografts (PDX). For classification of preclinical models, we developed an optimized classifier enriched for cancer cell–intrinsic gene expression signals, and performed high-throughput in vitro drug screening (n = 459 drugs) to analyze subtype-specific drug sensitivities. Results: The distinct molecular and clinicopathologic characteristics of each CMS group were validated in a single-hospital series of 409 primary colorectal cancers. The new, cancer cell–adapted classifier was found to perform well in primary tumors, and applied to a panel of 148 cell lines and 32 PDXs, these colorectal cancer models were shown to recapitulate the biology of the CMS groups. Drug screening of 33 cell lines demonstrated subtype-dependent response profiles, confirming strong response to EGFR and HER2 inhibitors in the CMS2 epithelial/canonical group, and revealing strong sensitivity to HSP90 inhibitors in cells with the CMS1 microsatellite instability/immune and CMS4 mesenchymal phenotypes. This association was validated in vitro in additional CMS-predicted cell lines. Combination treatment with 5-fluorouracil and luminespib showed potential to alleviate chemoresistance in a CMS4 PDX model, an effect not seen in a chemosensitive CMS2 PDX model. Conclusions: We provide translation of CMS classification to preclinical models and uncover a potential for targeted treatment repurposing in the chemoresistant CMS4 group. Clin Cancer Res; 24(4); 794–806. ©2017 AACR.

CMS-dependent prognostic impact of KRAS and BRAFV600E mutations in primary colorectal cancer

Abstract Background The prognostic impact of KRAS and BRAFV600E mutations in primary colorectal cancer (CRC) varies with microsatellite instability (MSI) status. The gene expression–based consensus molecular subtypes (CMSs) of CRC define molecularly and clinically distinct subgroups, and represent a novel stratification framework in biomarker analysis. We investigated the prognostic value of these mutations within the CMS groups. Patients and methods Totally 1197 primary tumors from a Norwegian series of CRC stage I–IV were analyzed for MSI and mutation status in hotspots in KRAS (codons 12, 13 and 61) and BRAF (codon 600). A subset was analyzed for gene expression and confident CMS classification was obtained for 317 samples. This cohort was expanded with clinical and molecular data, including CMS classification, from 514 patients in the publically available dataset GSE39582. Gene expression signatures associated with KRAS and BRAFV600E mutations were used to evaluate differential impact of mutations on gene expression among the CMS groups. Results BRAF V600E and KRAS mutations were both associated with inferior 5-year overall survival (OS) exclusively in MSS tumors (BRAFV600E mutation versus KRAS/BRAF wild-type: Hazard ratio (HR) 2.85, P < 0.001; KRAS mutation versus KRAS/BRAF wild-type: HR 1.30, P = 0.013). BRAFV600E-mutated MSS tumors were strongly enriched and associated with metastatic disease in CMS1, leading to negative prognostic impact in this subtype (OS: BRAFV600E mutation versus wild-type: HR 7.73, P = 0.001). In contrast, the poor prognosis of KRAS mutations was limited to MSS tumors with CMS2/CMS3 epithelial-like gene expression profiles (OS: KRAS mutation versus wild-type: HR 1.51, P = 0.011). The subtype-specific prognostic associations were substantiated by differential effects of BRAFV600E and KRAS mutations on gene expression signatures according to the MSI status and CMS group. Conclusions BRAF V600E mutations are enriched and associated with metastatic disease in CMS1 MSS tumors, leading to poor prognosis in this subtype. KRAS mutations are associated with adverse outcome in epithelial (CMS2/CMS3) MSS tumors.

Drug sensitivity and resistance testing identifies PLK1 inhibitors and gemcitabine as potent drugs for malignant peripheral nerve sheath tumors

Patients with malignant peripheral nerve sheath tumor (MPNST), a rare soft tissue cancer associated with loss of the tumor suppressor neurofibromin (NF1), have poor prognosis and typically respond poorly to adjuvant therapy. We evaluated the effect of 299 clinical and investigational compounds on seven MPNST cell lines, two primary cultures of human Schwann cells, and five normal bone marrow aspirates, to identify potent drugs for MPNST treatment with few side effects. Top hits included Polo‐like kinase 1 (PLK1) inhibitors (volasertib and BI2536) and the fluoronucleoside gemcitabine, which were validated in orthogonal assays measuring viability, cytotoxicity, and apoptosis. DNA copy number, gene expression, and protein expression were determined for the cell lines to assess pharmacogenomic relationships. MPNST cells were more sensitive to BI2536 and gemcitabine compared to a reference set of 94 cancer cell lines. PLK1, RRM1, and RRM2 mRNA levels were increased in MPNST compared to benign neurofibroma tissue, and the protein level of PLK1 was increased in the MPNST cell lines compared to normal Schwann cells, indicating an increased dependence on these drug targets in malignant cells. Furthermore, we observed an association between increased mRNA expression of PLK1, RRM1, and RRM2 in patient samples and worse disease outcome, suggesting a selective benefit from inhibition of these genes in the most aggressive tumors.

Prognostic, predictive, and pharmacogenomic assessments of CDX2 refine stratification of colorectal cancer

We aimed to refine the value of CDX2 as an independent prognostic and predictive biomarker in colorectal cancer (CRC) according to disease stage and chemotherapy sensitivity in preclinical models. CDX2 expression was evaluated in 1045 stage I–IV primary CRCs by gene expression (n = 403) or immunohistochemistry (n = 642) and in relation to 5‐year relapse‐free survival (RFS), overall survival (OS), and chemotherapy. Pharmacogenomic associations between CDX2 expression and 69 chemotherapeutics were assessed by drug screening of 35 CRC cell lines. CDX2 expression was lost in 11.6% of cases and showed independent poor prognostic value in multivariable models. For individual stages, CDX2 was prognostic only in stage IV, independent of chemotherapy. Among stage I–III patients not treated in an adjuvant setting, CDX2 loss was associated with a particularly poor survival in the BRAF‐mutated subgroup, but prognostic value was independent of microsatellite instability status and the consensus molecular subtypes. In stage III, the 5‐year RFS rate was higher among patients with loss of CDX2 who received adjuvant chemotherapy than among patients who did not. The CDX2‐negative cell lines were significantly more sensitive to chemotherapeutics than CDX2‐positive cells, and the multidrug resistance genes MDR1 and CFTR were significantly downregulated both in CDX2‐negative cells and in patient tumors. Loss of CDX2 in CRC is an adverse prognostic biomarker only in stage IV disease and appears to be associated with benefit from adjuvant chemotherapy in stage III. Early‐stage patients not qualifying for chemotherapy might be reconsidered for such treatment if their tumor has loss of CDX2 and mutated BRAF.

Inferior survival for patients with malignant peripheral nerve sheath tumors defined by aberrant TP53

Malignant peripheral nerve sheath tumor is a rare and aggressive disease with poor treatment response, mainly affecting adolescents and young adults. Few molecular biomarkers are used in the management of this cancer type, and although TP53 is one of few recurrently mutated genes in malignant peripheral nerve sheath tumor, the mutation prevalence and the corresponding clinical value of the TP53 network remains unsettled. We present a multi-level molecular study focused on aberrations in the TP53 network in relation to patient outcome in a series of malignant peripheral nerve sheath tumors from 100 patients and 38 neurofibromas, including TP53 sequencing, high-resolution copy number analyses of TP53 and MDM2, and gene expression profiling. Point mutations in TP53 were accompanied by loss of heterozygosity, resulting in complete loss of protein function in 8.2% of the malignant peripheral nerve sheath tumors. Another 5.5% had MDM2 amplification. TP53 mutation and MDM2 amplification were mutually exclusive and patients with either type of aberration in their tumor had a worse prognosis, compared to those without (hazard ratio for 5-year disease-specific survival 3.5, 95% confidence interval 1.78–6.98). Both aberrations had similar consequences on the gene expression level, as analyzed by a TP53-associated gene signature, a property also shared with the copy number aberrations and/or loss of heterozygosity at the TP53 locus, suggesting a common “TP53-mutated phenotype” in as many as 60% of the tumors. This was a poor prognostic phenotype (hazard ratio = 4.1, confidence interval:1.7–9.8), thus revealing a TP53-non-aberrant patient subgroup with a favorable outcome. The frequency of the “TP53-mutated phenotype” warrants explorative studies of stratified treatment strategies in malignant peripheral nerve sheath tumor.

Alternative splicing expands the prognostic impact of KRAS in microsatellite stable primary colorectal cancer: KRAS splicing and prognosis in microsatellite stable colorectal cancer

KRAS mutation is a well‐known marker for poor response to targeted treatment and patient prognosis in microsatellite stable (MSS) colorectal cancer (CRC). However, variation in clinical outcomes among patients wild‐type for KRAS underlines that this is not a homogeneous population. Here, we evaluated the prognostic impact of KRAS alternative splicing in relation to mutation status in a single‐hospital series of primary MSS CRCs (N = 258). Using splicing‐sensitive microarrays and RNA sequencing, the relative expression of KRAS‐4A versus KRAS‐4B transcript variants was confirmed to be down‐regulated in CRC compared to normal colonic mucosa (N = 41; p ≤ 0.001). This was independent of mutation status, however, gene set enrichment analysis revealed that the effect of splicing on KRAS signaling was specific to the KRAS wild‐type subgroup, in which low relative KRAS‐4A expression was associated with a higher level of KRAS signaling (p = 0.005). In concordance, the prognostic value of KRAS splicing was also dependent on mutation status, and for patients with Stage I–III KRAS wild‐type MSS CRC, low relative KRAS‐4A expression was associated with inferior overall survival (HR: 2.36, 95% CI: 1.07–5.18, p = 0.033), a result not found in mutant cases (pinteraction = 0.026). The prognostic association in the wild‐type subgroup was independent of clinicopathological factors, including cancer stage in multivariable analysis (HR: 2.68, 95% CI: 1.18–6.09, p = 0.018). This suggests that KRAS has prognostic value beyond mutation status in MSS CRC, and highlights the importance of molecular heterogeneity in the clinically relevant KRAS wild‐type subgroup.