In colorectal cancer, oncogenic mutations transform a hierarchically organized and homeostatic epithelium into invasive cancer tissue lacking visible organization. We sought to define transcriptional states of colorectal cancer cells and signals controlling their development by performing single‐cell transcriptome analysis of tumors and matched non‐cancerous tissues of twelve colorectal cancer patients. We defined patient‐overarching colorectal cancer cell clusters characterized by differential activities of oncogenic signaling pathways such as mitogen‐activated protein kinase and oncogenic traits such as replication stress. RNA metabolic labeling and assessment of RNA velocity in patient‐derived organoids revealed developmental trajectories of colorectal cancer cells organized along a mitogen‐activated protein kinase activity gradient. This was in contrast to normal colon organoid cells developing along graded Wnt activity. Experimental targeting of EGFR‐BRAF‐MEK in cancer organoids affected signaling and gene expression contingent on predictive KRAS/BRAF mutations and induced cell plasticity overriding default developmental trajectories. Our results highlight directional cancer cell development as a driver of non‐genetic cancer cell heterogeneity and re‐routing of trajectories as a response to targeted therapy.

Background: Pancreatic neuroendocrine neoplasms (PanNENs) fall into two subclasses: the well-differentiated, low- to high-grade pancreatic neuroendocrine tumors (PanNETs), and the poorly-differentiated, high-grade pancreatic neuroendocrine carcinomas (PanNECs). While recent studies suggest an endocrine descent of PanNETs, the origin of PanNECs remains unknown.Methods: We performed DNA methylation analysis for 57 PanNEN samples and found that distinct methylation profiles separated PanNENs into two major groups, clearly distinguishing high-grade PanNECs from other PanNETs including high-grade NETG3. DNA alterations and immunohistochemistry of cell-type markers PDX1, ARX, and SOX9 were utilized to further characterize PanNECs and their cell of origin in the pancreas.Results: Phylo-epigenetic and cell-type signature features derived from alpha, beta, acinar, and ductal adult cells suggest an exocrine cell of origin for PanNECs, thus separating them in cell lineage from other PanNENs of endocrine origin.Conclusions: Our study provides a robust and clinically applicable method to clearly distinguish PanNECs from G3 PanNETs, improving patient stratification.

Most metastatic colorectal cancer (mCRC) patients succumb to refractory disease due to secondary chemotherapy resistance. To elucidate the molecular changes associated with secondary resistance, we recruited 64 patients with mCRC and hepatic metastases before standard first-line chemotherapy between 2014 and 2018. We subjected DNA from primary tumor specimens (P), hepatic metastasis specimens after treatment (M), and liquid biopsies (L) taken prior to (pre), during (intra), and after (post) treatment to next generation sequencing. We performed Nanostring expression analysis in P and M specimens. Comparative bioinformatics and statistical analysis revealed typical mutational patterns with frequent alterations in TP53, APC, and KRAS in P specimens (n = 48). P and pre-L (n = 42), as well as matched P and M (n = 30), displayed a similar mutation spectrum. In contrast, gene expression profiles classified P (n = 31) and M (n = 23), distinguishable by up-regulation of immune/cytokine receptor and autophagy programs. Switching of consensus molecular subtypes from P to M occurred in 58.3% of cases. M signature genes SFRP2 and SPP1 associated with inferior survival, as validated in an independent cohort. Molecular changes during first-line treatment were detectable by expression profiling rather than by mutational tumor and liquid biopsy analyses. SFRP2 and SPP1 may serve as biomarkers and/or actionable targets.

Background: Colorectal cancer (CRC) development is generally accepted as a sequential process, with genetic mutations determining phenotypic tumor progression. However, matching genetic profiles with histological transition requires the analyses of temporal samples from the same patient at key stages of progression. Results: Here, we compared the genetic profiles of 34 early carcinomas with their respective adenomatous precursors to assess timing and heterogeneity of driver alterations accompanying the switch from benign adenoma to malignant carcinoma. In almost half of the cases, driver mutations specific to the carcinoma stage were not observed. In samples where carcinoma-specific alterations were present, TP53 mutations and chromosome 20 copy gains commonly accompanied the switch from adenomatous tissue to carcinoma. Remarkably, 40% and 50% of high-grade adenomas shared TP53 mutations and chromosome 20 gains, respectively, with their matched carcinomas. In addition, multi-regional analyses revealed greater heterogeneity of driver mutations in adenomas compared to their matched carcinomas. Conclusion: Genetic alterations in TP53 and chromosome 20 occur at the earliest histological stage in colorectal carcinomas (pTis and pT1). However, high-grade adenomas can share these alterations despite their histological distinction. Based on the well-defined sequence of CRC development, we suggest that the timing of genetic changes during neoplastic progression is frequently uncoupled from histological progression.

Purpose: We evaluated additional mutations in RAS wild-type (WT) metastatic colorectal cancer (mCRC) as prognostic and predictive biomarkers for the efficacy of added panitumumab to a 5-fluorouracil plus folinic acid (FU/FA) maintenance as pre-specified analysis of the randomized PanaMa trial. Patients and Methods: Mutations (MUT) were identified using targeted next-generation sequencing (NGS; Illumina Cancer Hotspot Panel v2) and IHC. RAS/BRAF V600E/PIK3CA/AKT1/ALK1/ERBB2/PTEN MUT and HER2/neu overexpressions were negatively hyperselected and correlated with median progression-free survival (PFS) and overall survival (OS) since start of maintenance treatment, and objective response rates (ORR). Univariate/multivariate Cox regression estimated hazard ratios (HR) and 95% confidence intervals (CI). Results: 202 of 248 patients (81.5%) of the full analysis set (FAS) had available NGS data: hyperselection WT, 162 (80.2%); MUT, 40 (19.8%). From start of maintenance therapy, hyperselection WT tumors were associated with longer median PFS as compared with hyperselection MUT mCRC (7.5 vs. 5.4 months; HR, 0.75; 95% CI, 0.52–1.07; P = 0.11), OS (28.7 vs. 22.2 months; HR, 0.53; 95% CI, 0.36–0.77; P = 0.001), and higher ORR (35.8% vs. 25.0%, P = 0.26). The addition of panitumumab to maintenance was associated with significant benefit in hyperselection WT tumors for PFS (9.2 vs. 6.0 months; HR, 0.66; 95% CI, 0.47–0.93; P = 0.02) and numerically also for OS (36.9 vs. 24.9 months; HR, 0.91; 95% CI, 0.61–1.36; P = 0.50), but not in hyperselection MUT tumors. Hyperselection status interacted with maintenance treatment arms in terms of PFS (P = 0.06) and OS (P = 0.009). Conclusions: Extended molecular profiling beyond RAS may have the potential to improve the patient selection for anti-EGFR containing maintenance regimens.

Recurrence of tumor cells following local and systemic therapy is a significant hurdle in cancer. Most patients with metastatic colorectal cancer (mCRC) will relapse, despite resection of the metastatic lesions. A better understanding of the evolutionary history of recurrent lesions is required to identify the spatial and temporal patterns of metastatic progression and expose the genetic and evolutionary determinants of therapeutic resistance. With this goal in mind, here we leveraged a unique single-cell whole-genome sequencing dataset from recurrent hepatic lesions of an mCRC patient. Our phylogenetic analysis confirms that the treatment induced a severe demographic bottleneck in the liver metastasis but also that a previously diverged lineage survived this surgery, possibly after migration to a different site in the liver. This lineage evolved very slowly for two years under adjuvant drug therapy and diversified again in a very short period. We identified several non-silent mutations specific to this lineage and inferred a substantial contribution of chemotherapy to the overall, genome-wide mutational burden. All in all, our study suggests that mCRC subclones can migrate locally and evade resection, keep evolving despite rounds of chemotherapy, and re-expand explosively.