Identification of a neural development gene expression signature in colon cancer stem cells reveals a role for EGR2 in tumorigenesis.

Recent evidence demonstrates that colon cancer stem cells (CSCs) can generate neurons that synapse with tumor innervating fibers required for tumorigenesis and disease progression. Greater understanding of the mechanisms that regulate CSC driven tumor neurogenesis may therefore lead to more effective treatments. RNA-sequencing analyses of ALDHPositive CSCs from colon cancer patient-derived organoids (PDOs) and xenografts (PDXs) showed CSCs to be enriched for neural development genes. Functional analyses of genes differentially expressed in CSCs from PDO and PDX models demonstrated the neural crest stem cell (NCSC) regulator EGR2 to be required for tumor growth and to control expression of homebox superfamily embryonic master transcriptional regulator HOX genes and the neural stem cell and master cell fate regulator SOX2. These data support CSCs as the source of tumor neurogenesis and suggest that targeting EGR2 may provide a therapeutic differentiation strategy to eliminate CSCs and block nervous system driven disease progression.

RNA sequencing of long-term label-retaining colon cancer stem cells identifies novel regulators of quiescence.

Recent data suggest that therapy-resistant quiescent cancer stem cells (qCSCs) are the source of relapse in colon cancer. Here, using colon cancer patient-derived organoids and xenografts, we identify rare long-term label-retaining qCSCs that can re-enter the cell cycle to generate new tumors. RNA sequencing analyses demonstrated that these cells display the molecular hallmarks of quiescent tissue stem cells, including expression of p53 signaling genes, and are enriched for transcripts common to damage-induced quiescent revival stem cells of the regenerating intestine. In addition, we identify negative regulators of cell cycle, downstream of p53, that we show are indicators of poor prognosis and may be targeted for qCSC abolition in both p53 wild-type and mutant tumors. These data support the temporal inhibition of downstream targets of p53 signaling, in combination with standard-of-care treatments, for the elimination of qCSCs and prevention of relapse in colon cancer.

Generation of Patient-Derived Colorectal Cancer Organoids for RAS Studies.

Human cell line models have been widely used for testing of novel anticancer compounds and for predicting clinical response to monotherapies and combinatorial therapies. For many years, standard monolayer culture conditions were used as gold standard, only surpassed by in vivo testing of mouse models. Recently, the incorporation of three-dimensional culture has been shown to further improve predictive compound testing. In view of the renewed interest in anti-RAS cancer therapy, we provide a protocol for establishing colorectal cancer organoids which are characterized by a high prevalence of KRAS mutations.

Heterogeneous pathway activation and drug response modelled in colorectal-tumor-derived 3D cultures.

Organoid cultures derived from colorectal cancer (CRC) samples are increasingly used as preclinical models for studying tumor biology and the effects of targeted therapies under conditions capturing in vitro the genetic make-up of heterogeneous and even individual neoplasms. While 3D cultures are initiated from surgical specimens comprising multiple cell populations, the impact of tumor heterogeneity on drug effects in organoid cultures has not been addressed systematically. Here we have used a cohort of well-characterized CRC organoids to study the influence of tumor heterogeneity on the activity of the KRAS/MAPK-signaling pathway and the consequences of treatment by inhibitors targeting EGFR and downstream effectors. MAPK signaling, analyzed by targeted proteomics, shows unexpected heterogeneity irrespective of RAS mutations and is associated with variable responses to EGFR inhibition. In addition, we obtained evidence for intratumoral heterogeneity in drug response among parallel "sibling" 3D cultures established from a single KRAS-mutant CRC. Our results imply that separate testing of drug effects in multiple subpopulations may help to elucidate molecular correlates of tumor heterogeneity and to improve therapy response prediction in patients.

Non-Canonical Hedgehog Signaling Is a Positive Regulator of the WNT Pathway and Is Required for the Survival of Colon Cancer Stem Cells.

Colon cancer is a heterogeneous tumor driven by a subpopulation of cancer stem cells (CSCs). To study CSCs in colon cancer, we used limiting dilution spheroid and serial xenotransplantation assays to functionally define the frequency of CSCs in a panel of patient-derived cancer organoids. These studies demonstrated cancer organoids to be enriched for CSCs, which varied in frequency between tumors. Whole-transcriptome analysis identified WNT and Hedgehog signaling components to be enhanced in CSC-enriched tumors and in aldehyde dehydrogenase (ALDH)-positive CSCs. Canonical GLI-dependent Hedgehog signaling is a negative regulator of WNT signaling in normal intestine and intestinal tumors. Here, we show that Hedgehog signaling in colon CSCs is autocrine SHH-dependent, non-canonical PTCH1 dependent, and GLI independent. In addition, using small-molecule inhibitors and RNAi against SHH-palmitoylating Hedgehog acyltransferase (HHAT), we demonstrate that non-canonical Hedgehog signaling is a positive regulator of WNT signaling and required for colon CSC survival.

Molecular dissection of colorectal cancer in pre-clinical models identifies biomarkers predicting sensitivity to EGFR inhibitors.

Colorectal carcinoma represents a heterogeneous entity, with only a fraction of the tumours responding to available therapies, requiring a better molecular understanding of the disease in precision oncology. To address this challenge, the OncoTrack consortium recruited 106 CRC patients (stages I-IV) and developed a pre-clinical platform generating a compendium of drug sensitivity data totalling >4,000 assays testing 16 clinical drugs on patient-derived in vivo and in vitro models. This large biobank of 106 tumours, 35 organoids and 59 xenografts, with extensive omics data comparing donor tumours and derived models provides a resource for advancing our understanding of CRC. Models recapitulate many of the genetic and transcriptomic features of the donors, but defined less complex molecular sub-groups because of the loss of human stroma. Linking molecular profiles with drug sensitivity patterns identifies novel biomarkers, including a signature outperforming RAS/RAF mutations in predicting sensitivity to the EGFR inhibitor cetuximab.

DNA copy number changes define spatial patterns of heterogeneity in colorectal cancer.

Genetic heterogeneity between and within tumours is a major factor determining cancer progression and therapy response. Here we examined DNA sequence and DNA copy-number heterogeneity in colorectal cancer (CRC) by targeted high-depth sequencing of 100 most frequently altered genes. In 97 samples, with primary tumours and matched metastases from 27 patients, we observe inter-tumour concordance for coding mutations; in contrast, gene copy numbers are highly discordant between primary tumours and metastases as validated by fluorescent in situ hybridization. To further investigate intra-tumour heterogeneity, we dissected a single tumour into 68 spatially defined samples and sequenced them separately. We identify evenly distributed coding mutations in APC and TP53 in all tumour areas, yet highly variable gene copy numbers in numerous genes. 3D morpho-molecular reconstruction reveals two clusters with divergent copy number aberrations along the proximal-distal axis indicating that DNA copy number variations are a major source of tumour heterogeneity in CRC.

Assay Establishment and Validation of a High-Throughput Screening Platform for Three-Dimensional Patient-Derived Colon Cancer Organoid Cultures.

The application of patient-derived three-dimensional culture systems as disease-specific drug sensitivity models has enormous potential to connect compound screening and clinical trials. However, the implementation of complex cell-based assay systems in drug discovery requires reliable and robust screening platforms. Here we describe the establishment of an automated platform in 384-well format for three-dimensional organoid cultures derived from colon cancer patients. Single cells were embedded in an extracellular matrix by an automated workflow and subsequently self-organized into organoid structures within 4 days of culture before being exposed to compound treatment. We performed validation of assay robustness and reproducibility via plate uniformity and replicate-experiment studies. After assay optimization, the patient-derived organoid platform passed all relevant validation criteria. In addition, we introduced a streamlined plate uniformity study to evaluate patient-derived colon cancer samples from different donors. Our results demonstrate the feasibility of using patient-derived tumor samples for high-throughput assays and their integration as disease-specific models in drug discovery.

Bi-PROF: bisulfite profiling of target regions using 454 GS FLX Titanium technology.

The use of next generation sequencing has expanded our view on whole mammalian methylome patterns. In particular, it provides a genome-wide insight of local DNA methylation diversity at single nucleotide level and enables the examination of single chromosome sequence sections at a sufficient statistical power. We describe a bisulfite-based sequence profiling pipeline, Bi-PROF, which is based on the 454 GS-FLX Titanium technology that allows to obtain up to one million sequence stretches at single base pair resolution without laborious subcloning. To illustrate the performance of the experimental workflow connected to a bioinformatics program pipeline (BiQ Analyzer HT) we present a test analysis set of 68 different epigenetic marker regions (amplicons) in five individual patient-derived xenograft tissue samples of colorectal cancer and one healthy colon epithelium sample as a control. After the 454 GS-FLX Titanium run, sequence read processing and sample decoding, the obtained alignments are quality controlled and statistically evaluated. Comprehensive methylation pattern interpretation (profiling) assessed by analyzing 10 (2)-10 (4) sequence reads per amplicon allows an unprecedented deep view on pattern formation and methylation marker heterogeneity in tissues concerned by complex diseases like cancer.

Network quantification of EGFR signaling unveils potential for targeted combination therapy.

The epidermal growth factor receptor (EGFR) signaling network is activated in most solid tumors, and small-molecule drugs targeting this network are increasingly available. However, often only specific combinations of inhibitors are effective. Therefore, the prediction of potent combinatorial treatments is a major challenge in targeted cancer therapy. In this study, we demonstrate how a model-based evaluation of signaling data can assist in finding the most suitable treatment combination. We generated a perturbation data set by monitoring the response of RAS/PI3K signaling to combined stimulations and inhibitions in a panel of colorectal cancer cell lines, which we analyzed using mathematical models. We detected that a negative feedback involving EGFR mediates strong cross talk from ERK to AKT. Consequently, when inhibiting MAPK, AKT activity is increased in an EGFR-dependent manner. Using the model, we predict that in contrast to single inhibition, combined inactivation of MEK and EGFR could inactivate both endpoints of RAS, ERK and AKT. We further could demonstrate that this combination blocked cell growth in BRAF- as well as KRAS-mutated tumor cells, which we confirmed using a xenograft model.