Understanding patient-derived tumor organoid growth through an integrated imaging and mathematical modeling framework.

Patient-derived tumor organoids (PDTOs) are novel cellular models that maintain the genetic, phenotypic and structural features of patient tumor tissue and are useful for studying tumorigenesis and drug response. When integrated with advanced 3D imaging and analysis techniques, PDTOs can be used to establish physiologically relevant high-throughput and high-content drug screening platforms that support the development of patient-specific treatment strategies. However, in order to effectively leverage high-throughput PDTO observations for clinical predictions, it is critical to establish a quantitative understanding of the basic properties and variability of organoid growth dynamics. In this work, we introduced an innovative workflow for analyzing and understanding PDTO growth dynamics, by integrating a high-throughput imaging deep learning platform with mathematical modeling, incorporating flexible growth laws and variable dormancy times. We applied the workflow to colon cancer organoids and demonstrated that organoid growth is well-described by the Gompertz model of growth. Our analysis showed significant intrapatient heterogeneity in PDTO growth dynamics, with the initial exponential growth rate of an organoid following a lognormal distribution within each dataset. The level of intrapatient heterogeneity varied between patients, as did organoid growth rates and dormancy times of single seeded cells. Our work contributes to an emerging understanding of the basic growth characteristics of PDTOs, and it highlights the heterogeneity in organoid growth both within and between patients. These results pave the way for further modeling efforts aimed at predicting treatment response dynamics and drug resistance timing.

Merging Metabolic Modeling and Imaging for Screening Therapeutic Targets in Colorectal Cancer.

Cancer-associated fibroblasts (CAFs) play a key role in metabolic reprogramming and are well-established contributors to drug resistance in colorectal cancer (CRC). To exploit this metabolic crosstalk, we integrated a systems biology approach that identified key metabolic targets in a data-driven method and validated them experimentally. This process involved high-throughput computational screening to investigate the effects of enzyme perturbations predicted by a computational model of CRC metabolism to understand system-wide effects efficiently. Our results highlighted hexokinase (HK) as one of the crucial targets, which subsequently became our focus for experimental validation using patient-derived tumor organoids (PDTOs). Through metabolic imaging and viability assays, we found that PDTOs cultured in CAF conditioned media exhibited increased sensitivity to HK inhibition. Our approach emphasizes the critical role of integrating computational and experimental techniques in exploring and exploiting CRC-CAF crosstalk.

Targeting circadian transcriptional programs in triple negative breast cancer through a cis-regulatory mechanism.

Circadian clock genes are emerging targets in many types of cancer, but their mechanistic contributions to tumor progression are still largely unknown. This makes it challenging to stratify patient populations and develop corresponding treatments. In this work, we show that in breast cancer, the disrupted expression of circadian genes has the potential to serve as biomarkers. We also show that the master circadian transcription factors (TFs) BMAL1 and CLOCK are required for the proliferation of metastatic mesenchymal stem-like (mMSL) triple-negative breast cancer (TNBC) cells. Using currently available small molecule modulators, we found that a stabilizer of cryptochrome 2 (CRY2), the direct repressor of BMAL1 and CLOCK transcriptional activity, synergizes with inhibitors of proteasome, which is required for BMAL1 and CLOCK function, to repress a transcriptional program comprising circadian cycling genes in mMSL TNBC cells. Omics analyses on drug-treated cells implied that this repression of transcription is mediated by the transcription factor binding sites (TFBSs) features in the cis-regulatory elements (CRE) of clock-controlled genes. Through a massive parallel reporter assay, we defined a set of CRE features that are potentially repressed by the specific drug combination. The identification of cis -element enrichment may serve as a new way of defining and targeting tumor types through the modulation of cis -regulatory programs, and ultimately provide a new paradigm of therapy design for cancer types with unclear drivers like TNBC.

A novel thin plate spline methodology to model tissue surfaces and quantify tumor cell invasion in organ-on-chip models.

Organ-on-chip (OOC) models can be useful tools for cancer drug discovery. Advances in OOC technology have led to the development of more complex assays, yet analysis of these systems does not always account for these advancements, resulting in technical challenges. A challenging task in the analysis of these two-channel microfluidic models is to define the boundary between the channels so objects moving within and between channels can be quantified. We propose a novel imaging-based application of a thin plate spline method - a generalized cubic spline that can be used to model coordinate transformations - to model a tissue boundary and define compartments for quantification of invaded objects, representing the early steps in cancer metastasis. To evaluate its performance, we applied our analytical approach to an adapted OOC developed by Emulate, Inc., utilizing a two-channel system with endothelial cells in the bottom channel and colorectal cancer (CRC) patient-derived organoids (PDOs) in the top channel. Initial application and visualization of this method revealed boundary variations due to microscope stage tilt and ridge and valley-like contours in the endothelial tissue surface. The method was functionalized into a reproducible analytical process and web tool - the Chip Invasion and Contour Analysis (ChICA) - to model the endothelial surface and quantify invading tumor cells across multiple chips. To illustrate applicability of the analytical method, we applied the tool to CRC organoid-chips seeded with two different endothelial cell types and measured distinct variations in endothelial surfaces and tumor cell invasion dynamics. Since ChICA utilizes only positional data output from imaging software, the method is applicable to and agnostic of the imaging tool and image analysis system used. The novel thin plate spline method developed in ChICA can account for variation introduced in OOC manufacturing or during the experimental workflow, can quickly and accurately measure tumor cell invasion, and can be used to explore biological mechanisms in drug discovery.

CCR5 and CCL5 gene expression in colorectal cancer: comprehensive profiling and clinical value.

BACKGROUND: The C-C motif chemokine receptor 5 (CCR5)/C-C motif chemokine ligand 5 (CCL5) axis plays a major role in colorectal cancer (CRC). We aimed to characterize the molecular features associated with CCR5/CCL5 expression in CRC and to determine whether CCR5/CCL5 levels could impact treatment outcomes. METHODS: 7604 CRCs tested with NextGen Sequencing on DNA and RNA were analyzed. Molecular features were evaluated according to CCR5 and CCL5 tumor gene expression quartiles. The impact on treatment outcomes was assessed in two cohorts, including 6341 real-world patients and 429 patients from the Cancer and Leukemia Group B (CALGB)/SWOG 80405 trial. RESULTS: CCR5/CCL5 expression was higher in right-sided versus left-sided tumors, and positively associated with consensus molecular subtypes 1 and 4. Higher CCR5/CCL5 expression was associated with higher tumor mutational burden, deficiency in mismatch repair and programmed cell death ligand 1 (PD-L1) levels. Additionally, high CCR5/CCL5 were associated with higher immune cell infiltration in the tumor microenvironment (TME) of MMR proficient tumors. Ingenuity pathway analysis revealed upregulation of the programmed cell death protein 1 (PD-1)/PD-L1 cancer immunotherapy pathway, phosphatase and tensin homolog (PTEN) and peroxisome proliferator-activated receptors (PPAR) signaling, and cytotoxic T-lymphocyte antigen 4 (CTLA-4) signaling in cytotoxic T lymphocytes, whereas several inflammation-related pathways were downregulated. Low CCR5/CCL5 expression was associated with increased benefit from cetuximab-FOLFOX treatment in the CALGB/SWOG 80405 trial, where significant treatment interaction was observed with biologic agents and chemotherapy backbone. CONCLUSIONS: Our data show a strong association between CCR5/CCL5 gene expression and distinct molecular features, gene expression profiles, TME cell infiltration, and treatment benefit in CRC. Targeting the CCR5/CCL5 axis may have clinical applications in selected CRC subgroups and may play a key role in developing and deploying strategies to modulate the immune TME for CRC treatment.

A novel thin plate spline methodology to model tissue surfaces and quantify tumor cell invasion in organ-on-chip models.

Organ-on-chip (OOC) models can be useful tools for cancer drug discovery. Advances in OOC technology have led to the development of more complex assays, yet analysis of these systems does not always account for these advancements, resulting in technical challenges. A challenging task in the analysis of these two-channel microfluidic models is to define the boundary between the channels so objects moving within and between channels can be quantified. We propose a novel imaging-based application of a thin plate spline method - a generalized cubic spline that can be used to model coordinate transformations - to model a tissue boundary and define compartments for quantification of invaded objects, representing the early steps in cancer metastasis. To evaluate its performance, we applied our analytical approach to an adapted OOC developed by Emulate, Inc., utilizing a two-channel system with endothelial cells in the bottom channel and colorectal cancer (CRC) patient-derived organoids (PDOs) in the top channel. Initial application and visualization of this method revealed boundary variations due to microscope stage tilt and ridge and valley-like contours in the endothelial tissue surface. The method was functionalized into a reproducible analytical process and web tool - the Chip Invasion and Contour Analysis (ChICA) - to model the endothelial surface and quantify invading tumor cells across multiple chips. To illustrate applicability of the analytical method, we applied the tool to CRC organoid-chips seeded with two different endothelial cell types and measured distinct variations in endothelial surfaces and tumor cell invasion dynamics. Since ChICA utilizes only positional data output from imaging software, the method is applicable to and agnostic of the imaging tool and image analysis system used. The novel thin plate spline method developed in ChICA can account for variation introduced in OOC manufacturing or during the experimental workflow, can quickly and accurately measure tumor cell invasion, and can be used to explore biological mechanisms in drug discovery.

Integration of Patient-Derived Organoids and Organ-on-Chip Systems: Investigating Colorectal Cancer Invasion within the Mechanical and GABAergic Tumor Microenvironment.

Three-dimensional (3D) in vitro models are essential in cancer research, but they often neglect physical forces. In our study, we combined patient-derived tumor organoids with a microfluidic organ-on-chip system to investigate colorectal cancer (CRC) invasion in the tumor microenvironment (TME). This allowed us to create patient-specific tumor models and assess the impact of physical forces on cancer biology. Our findings showed that the organoid-on-chip models more closely resembled patient tumors at the transcriptional level, surpassing organoids alone. Using 'omics' methods and live-cell imaging, we observed heightened responsiveness of KRAS mutant tumors to TME mechanical forces. These tumors also utilized the γ-aminobutyric acid (GABA) neurotransmitter as an energy source, increasing their invasiveness. This bioengineered model holds promise for advancing our understanding of cancer progression and improving CRC treatments.

Phosphorylation and stabilization of EZH2 by DCAF1/VprBP trigger aberrant gene silencing in colon cancer.

Our recent work has shown that DCAF1 (also known as VprBP) is overexpressed in colon cancer and phosphorylates histone H2AT120 to drive epigenetic gene inactivation and oncogenic transformation. We have extended these observations by investigating whether DCAF1 also phosphorylates non-histone proteins as an additional mechanism linking its kinase activity to colon cancer development. We now demonstrate that DCAF1 phosphorylates EZH2 at T367 to augment its nuclear stabilization and enzymatic activity in colon cancer cells. Consistent with this mechanistic role, DCAF1-mediated EZH2 phosphorylation leads to elevated levels of H3K27me3 and altered expression of growth regulatory genes in cancer cells. Furthermore, our preclinical studies using organoid and xenograft models revealed that EZH2 requires phosphorylation for its oncogenic function, which may have therapeutic implications for gene reactivation in colon cancer cells. Together, our data define a mechanism underlying DCAF1-driven colonic tumorigenesis by linking DCAF1-mediated EZH2 phosphorylation to EZH2 stability that is crucial for establishing H3K27me3 and gene silencing program.

Phenotypic deconvolution in heterogeneous cancer cell populations using drug-screening data.

Tumor heterogeneity is an important driver of treatment failure in cancer since therapies often select for drug-tolerant or drug-resistant cellular subpopulations that drive tumor growth and recurrence. Profiling the drug-response heterogeneity of tumor samples using traditional genomic deconvolution methods has yielded limited results, due in part to the imperfect mapping between genomic variation and functional characteristics. Here, we leverage mechanistic population modeling to develop a statistical framework for profiling phenotypic heterogeneity from standard drug-screen data on bulk tumor samples. This method, called PhenoPop, reliably identifies tumor subpopulations exhibiting differential drug responses and estimates their drug sensitivities and frequencies within the bulk population. We apply PhenoPop to synthetically generated cell populations, mixed cell-line experiments, and multiple myeloma patient samples and demonstrate how it can provide individualized predictions of tumor growth under candidate therapies. This methodology can also be applied to deconvolution problems in a variety of biological settings beyond cancer drug response.

Organs-on-chips: a decade of innovation.

Recreating 'living organs' with groundbreaking organ-on-a-chip (OOC) technologies is facilitating a new era of drug discovery. Studies by Huh et al. and Ronaldson-Bouchard et al. underscore advances made over a decade, spanning single organ functionality to interconnected organs, that enable examination of drug toxicities and disease pathogenesis in reconstituted tissues.

Breast cancer and neurotransmitters: emerging insights on mechanisms and therapeutic directions.

Exploring the relationship between various neurotransmitters and breast cancer cell growth has revealed their likely centrality to improving breast cancer treatment. Neurotransmitters play a key role in breast cancer biology through their effects on the cell cycle, epithelial mesenchymal transition, angiogenesis, inflammation, the tumor microenvironment and other pathways. Neurotransmitters and their receptors are vital to the initiation, progression and drug resistance of cancer and progress in our biological understanding may point the way to lower-cost and lower-risk antitumor therapeutic strategies. This review discusses multiple neurotransmitters in the context of breast cancer. It also discusses risk factors, repurposing of pharmaceuticals impacting neurotransmitter pathways, and the opportunity for better integrated models that encompass exercise, the intestinal microbiome, and other non-pharmacologic considerations. Neurotransmitters' role in breast cancer should no longer be ignored; it may appear to complicate the molecular picture but the ubiquity of neurotransmitters and their wide-ranging impacts provide an organizing framework upon which further understanding and progress against breast cancer can be based.

Neurotransmitter signaling: a new frontier in colorectal cancer biology and treatment.

The brain-gut axis, a bidirectional network between the central and enteric nervous system, plays a critical role in modulating the gastrointestinal tract function and homeostasis. Recently, increasing evidence suggests that neuronal signaling molecules can promote gastrointestinal cancers, however, the mechanisms remain unclear. Aberrant expression of neurotransmitter signaling genes in colorectal cancer supports the role of neurotransmitters to stimulate tumor growth and metastatic spread by promoting cell proliferation, migration, invasion, and angiogenesis. In addition, neurotransmitters can interact with immune and endothelial cells in the tumor microenvironment to promote inflammation and tumor progression. As such, pharmacological targeting of neurotransmitter signaling represent a promising novel anticancer approach. Here, we present an overview of the current evidence supporting the role of neurotransmitters in colorectal cancer biology and treatment.

Roadmap on plasticity and epigenetics in cancer.

The role of plasticity and epigenetics in shaping cancer evolution and response to therapy has taken center stage with recent technological advances including single cell sequencing. This roadmap article is focused on state-of-the-art mathematical and experimental approaches to interrogate plasticity in cancer, and addresses the following themes and questions: is there a formal overarching framework that encompasses both non-genetic plasticity and mutation-driven somatic evolution? How do we measure and model the role of the microenvironment in influencing/controlling non-genetic plasticity? How can we experimentally study non-genetic plasticity? Which mathematical techniques are required or best suited? What are the clinical and practical applications and implications of these concepts?

Elucidating tumor-stromal metabolic crosstalk in colorectal cancer through integration of constraint-based models and LC-MS metabolomics.

Colorectal cancer (CRC) is a major cause of morbidity and mortality in the United States. Tumor-stromal metabolic crosstalk in the tumor microenvironment promotes CRC development and progression, but exactly how stromal cells, in particular cancer-associated fibroblasts (CAFs), affect the metabolism of tumor cells remains unknown. Here we take a data-driven approach to investigate the metabolic interactions between CRC cells and CAFs, integrating constraint-based modeling and metabolomic profiling. Using metabolomics data, we perform unsteady-state parsimonious flux balance analysis to infer flux distributions for central carbon metabolism in CRC cells treated with or without CAF-conditioned media. We find that CAFs reprogram CRC metabolism through stimulation of glycolysis, the oxidative arm of the pentose phosphate pathway (PPP), and glutaminolysis, as well as inhibition of the tricarboxylic acid cycle. To identify potential therapeutic targets, we simulate enzyme knockouts and find that CAF-treated CRC cells are especially sensitive to inhibitions of hexokinase and glucose-6-phosphate, the rate limiting steps of glycolysis and oxidative PPP. Our work gives mechanistic insights into the metabolic interactions between CRC cells and CAFs and provides a framework for testing hypotheses towards CRC-targeted therapies.

Germline polymorphisms in genes maintaining the replication fork predict the efficacy of oxaliplatin and irinotecan in patients with metastatic colorectal cancer.

BACKGROUND: The TIMELESS-TIPIN complex protects the replication fork from replication stress induced by chemotherapeutic drugs. We hypothesised genetic polymorphisms of the TIMELESS-TIPIN complex may affect the response, progression-free survival (PFS), and overall survival (OS) of cytotoxic drugs in patients with metastatic colorectal cancer (mCRC). METHODS: We analysed data from the MAVERICC trial, which compared FOLFOX/bevacizumab and FOLFIRI/bevacizumab in untreated patients with mCRC. Genomic DNA extracted from blood samples was genotyped using an OncoArray. Eight functional single nucleotide polymorphisms (SNPs) in TIMELESS and TIPIN were tested for associations with clinical outcomes. RESULTS: In total, 324 patients were included (FOLFOX/bevacizumab arm, n = 161; FOLFIRI/bevacizumab arm, n = 163). In the FOLFOX/bevacizumab arm, no SNPs displayed confirmed associations with survival outcomes. In the FOLFIRI/bevacizumab arm, TIMELESS rs2291739 was significantly associated with OS in multivariate analysis (G/G vs. any A allele, hazard ratio = 3.06, 95% confidence interval = 1.49-6.25, p = 0.004). TIMELESS rs2291739 displayed significant interactions with treatment regarding both PFS and OS. CONCLUSIONS: TIMELESS rs2291739 might have different effects on therapeutic efficacy between oxaliplatin- and irinotecan-based chemotherapies. Upon further validation, our findings may be useful for personalised approaches in the first-line treatment of mCRC.

Reprogramming CBX8-PRC1 function with a positive allosteric modulator.

Canonical targeting of Polycomb repressive complex 1 (PRC1) to repress developmental genes is mediated by cell-type-specific, paralogous chromobox (CBX) proteins (CBX2, 4, 6, 7, and 8). Based on their central role in silencing and their dysregulation associated with human disease including cancer, CBX proteins are attractive targets for small-molecule chemical probe development. Here, we have used a quantitative and target-specific cellular assay to discover a potent positive allosteric modulator (PAM) of CBX8. The PAM activity of UNC7040 antagonizes H3K27me3 binding by CBX8 while increasing interactions with nucleic acids. We show that treatment with UNC7040 leads to efficient and selective eviction of CBX8-containing PRC1 from chromatin, loss of silencing, and reduced proliferation across different cancer cell lines. Our discovery and characterization of UNC7040 not only reveals the most cellularly potent CBX8-specific chemical probe to date, but also corroborates a mechanism of Polycomb regulation by non-specific CBX nucleotide binding activity.

Multiplexed imaging and effluent analysis to monitor cancer cell intravasation using a colorectal cancer-on-chip.

Despite colorectal cancer's (CRC) prevalence, its progression is not well understood. The microfluidic organ-on-chip (OOC) model described herein recreates the epithelial-endothelial tissue-tissue interface, fluid flow, and mechanical forces that exist in vivo , making it an attractive model to understand and ultimately disrupt CRC intravasation. This protocol provides step-by-step details for tumor cell seeding to create a CRC-on-chip model, chip effluent collection and analysis, and on-chip imaging to monitor tumor cell invasion within a more physiologically relevant microenvironment. For complete details on the use and execution of this protocol, please refer to Strelez et al. (2021).

VprBP directs epigenetic gene silencing through histone H2A phosphorylation in colon cancer.

Histone modification is aberrantly regulated in cancer and generates an unbalanced state of gene transcription. VprBP, a recently identified kinase, phosphorylates histone H2A on threonine 120 (T120) and is involved in oncogenic transcriptional dysregulation; however, its specific role in colon cancer is undefined. Here, we show that VprBP is overexpressed in colon cancer and directly contributes to epigenetic gene silencing and cancer pathogenesis. Mechanistically, the observed function of VprBP is mediated through H2AT120 phosphorylation (H2AT120p)-driven transcriptional repression of growth regulatory genes, resulting in a significantly higher proliferative capacity of colon cancer cells. Our preclinical studies using organoid and xenograft models demonstrate that treatment with the VprBP inhibitor B32B3 impairs colonic tumor growth by blocking H2AT120p and reactivating a transcriptional program resembling that of normal cells. Collectively, our work describes VprBP as a master kinase contributing to the development and progression of colon cancer, making it a new molecular target for novel therapeutic strategies.

Human colorectal cancer-on-chip model to study the microenvironmental influence on early metastatic spread.

Colorectal cancer (CRC) progression is a complex process that is not well understood. We describe an in vitro organ-on-chip model that emulates in vivo tissue structure and the tumor microenvironment (TME) to better understand intravasation, an early step in metastasis. The CRC-on-chip incorporates fluid flow and peristalsis-like cyclic stretching and consists of endothelial and epithelial compartments, separated by a porous membrane. On-chip imaging and effluent analyses are used to interrogate CRC progression and the resulting cellular heterogeneity. Mass spectrometry-based metabolite profiles are indicative of a CRC disease state. Tumor cells intravasate from the epithelial channel to the endothelial channel, revealing differences in invasion between aggressive and non-aggressive tumor cells. Tuning the TME by peristalsis-like mechanical forces, the epithelial:endothelial interface, and the addition of fibroblasts influences the invasive capabilities of tumor cells. The CRC-on-chip is a tunable human-relevant model system and a valuable tool to study early invasive events in cancer.

3D cell culture models: Drug pharmacokinetics, safety assessment, and regulatory consideration.

Nonclinical testing has served as a foundation for evaluating potential risks and effectiveness of investigational new drugs in humans. However, the current two-dimensional (2D) in vitro cell culture systems cannot accurately depict and simulate the rich environment and complex processes observed in vivo, whereas animal studies present significant drawbacks with inherited species-specific differences and low throughput for increased demands. To improve the nonclinical prediction of drug safety and efficacy, researchers continue to develop novel models to evaluate and promote the use of improved cell- and organ-based assays for more accurate representation of human susceptibility to drug response. Among others, the three-dimensional (3D) cell culture models present physiologically relevant cellular microenvironment and offer great promise for assessing drug disposition and pharmacokinetics (PKs) that influence drug safety and efficacy from an early stage of drug development. Currently, there are numerous different types of 3D culture systems, from simple spheroids to more complicated organoids and organs-on-chips, and from single-cell type static 3D models to cell co-culture 3D models equipped with microfluidic flow control as well as hybrid 3D systems that combine 2D culture with biomedical microelectromechanical systems. This article reviews the current application and challenges of 3D culture systems in drug PKs, safety, and efficacy assessment, and provides a focused discussion and regulatory perspectives on the liver-, intestine-, kidney-, and neuron-based 3D cellular models.