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?

SMaSH: Sample matching using SNPs in humans.

BACKGROUND: Inadvertent sample swaps are a real threat to data quality in any medium to large scale omics studies. While matches between samples from the same individual can in principle be identified from a few well characterized single nucleotide polymorphisms (SNPs), omics data types often only provide low to moderate coverage, thus requiring integration of evidence from a large number of SNPs to determine if two samples derive from the same individual or not. METHODS: We select about six thousand SNPs in the human genome and develop a Bayesian framework that is able to robustly identify sample matches between next generation sequencing data sets. RESULTS: We validate our approach on a variety of data sets. Most importantly, we show that our approach can establish identity between different omics data types such as Exome, RNA-Seq, and MethylCap-Seq. We demonstrate how identity detection degrades with sample quality and read coverage, but show that twenty million reads of a fairly low quality RNA-Seq sample are still sufficient for reliable sample identification. CONCLUSION: Our tool, SMASH, is able to identify sample mismatches in next generation sequencing data sets between different sequencing modalities and for low quality sequencing data.

FuSpot: a web-based tool for visual evaluation of fusion candidates.

BACKGROUND: Gene fusions often occur in cancer cells and in some cases are the main driver of oncogenesis. Correct identification of oncogenic gene fusions thus has implications for targeted cancer therapy. Recognition of this potential has led to the development of a myriad of sequencing-based fusion detection tools. However, given the same input, many of these detectors will find different fusion points or claim different sets of supporting data. Furthermore, the rate at which these tools falsely detect fusion events in data varies greatly. This discrepancy between tools underscores the fact that computation algorithms still cannot perfectly evaluate evidence; especially when provided with small amounts of supporting data as is typical in fusion detection. We assert that when evidence is provided in an easily digestible form, humans are more proficient in identifying true positives from false positives. RESULTS: We have developed a web tool that, given the genomic coordinates of a candidate fusion breakpoint, will extract fusion and non-fusion reads adjacent to the fusion point from partner transcripts, and color code reads by transcript origin and read orientation for ease of intuitive inspection by the user. Fusion partner transcript read alignments are performed using a novel variant of the Smith-Waterman algorithm. CONCLUSIONS: Combined with dynamic filtering parameters, the visualization provided by our tool introduces a powerful new investigative step that allows researchers to comprehensively evaluate fusion evidence. Additionally, this allows quick identification of false positives that may deceive most fusion detectors, thus eliminating unnecessary gene fusion validation. We apply our visualization tool to publicly available datasets and provide examples of true as well as false positives reported by open source fusion detection tools.

Statistical methods for detecting differentially methylated regions based on MethylCap-seq data.

DNA methylation is a well-established epigenetic mark, whose pattern throughout the genome, especially in the promoter or CpG islands, may be modified in a cell at a disease stage. Recently developed probabilistic approaches allow distributing methylation signals at nucleotide resolution from MethylCap-seq data. Standard statistical methods for detecting differential methylation suffer from 'curse of dimensionality' and sparsity in signals, resulting in high false-positive rates. Strong correlation of signals between CG sites also yields spurious results. In this article, we review applicability of high-dimensional mean vector tests for detection of differentially methylated regions (DMRs) and compare and contrast such tests with other methods for detecting DMRs. Comprehensive simulation studies are conducted to highlight the performance of these tests under different settings. Based on our observation, we make recommendations on the optimal test to use. We illustrate the superiority of mean vector tests in detecting cancer-related canonical gene pathways, which are significantly enriched for acute myeloid leukemia and ovarian cancer.

PrEMeR-CG: inferring nucleotide level DNA methylation values from MethylCap-seq data.

MOTIVATION: DNA methylation is an epigenetic change occurring in genomic CpG sequences that contribute to the regulation of gene transcription both in normal and malignant cells. Next-generation sequencing has been used to characterize DNA methylation status at the genome scale, but suffers from high sequencing cost in the case of whole-genome bisulfite sequencing, or from reduced resolution (inability to precisely define which of the CpGs are methylated) with capture-based techniques. RESULTS: Here we present a computational method that computes nucleotide-resolution methylation values from capture-based data by incorporating fragment length profiles into a model of methylation analysis. We demonstrate that it compares favorably with nucleotide-resolution bisulfite sequencing and has better predictive power with respect to a reference than window-based methods, often used for enrichment data. The described method was used to produce the methylation data used in tandem with gene expression to produce a novel and clinically significant gene signature in acute myeloid leukemia. In addition, we introduce a complementary statistical method that uses this nucleotide-resolution methylation data for detection of differentially methylated features.

Genome-wide methylation profiling in decitabine-treated patients with acute myeloid leukemia.

The outcome of older (≥ 60 years) acute myeloid leukemia (AML) patients is poor, and novel treatments are needed. In a phase 2 trial for older AML patients, low-dose (20 mg/m(2) per day for 10 days) decitabine, a DNA hypomethylating azanucleoside, produced 47% complete response rate with an excellent toxicity profile. To assess the genome-wide activity of decitabine, we profiled pretreatment and post treatment (day 25/course 1) methylomes of marrow samples from patients (n = 16) participating in the trial using deep-sequencing analysis of methylated DNA captured by methyl-binding protein (MBD2). Decitabine significantly reduced global methylation compared with pretreatment baseline (P = .001). Percent marrow blasts did not correlate with global methylation levels, suggesting that hypomethylation was related to the activity of decitabine rather than to a mere decrease in leukemia burden. Hypomethylation occurred predominantly in CpG islands and CpG island-associated regions (P ranged from .03 to .04) A significant concentration (P < .001) of the hypomehtylated CpG islands was found in chromosome subtelomeric regions, suggesting a differential activity of decitabine in distinct chromosome regions. Hypermethylation occurred much less frequently than hypomethylation and was associated with low CpG content regions. Decitabine-related methylation changes were concordant with those previously reported in distinct genes. In summary, our study supports the feasibility of methylome analyses as a pharmacodynamic endpoint for hypomethylating therapies.

Transcriptome free energy can serve as a dynamic patient-specific biomarker in acute myeloid leukemia.

Acute myeloid leukemia (AML) is prevalent in both adult and pediatric patients. Despite advances in patient categorization, the heterogeneity of AML remains a challenge. Recent studies have explored the use of gene expression data to enhance AML diagnosis and prognosis, however, alternative approaches rooted in physics and chemistry may provide another level of insight into AML transformation. Utilizing publicly available databases, we analyze 884 human and mouse blood and bone marrow samples. We employ a personalized medicine strategy, combining state-transition theory and surprisal analysis, to assess the RNA transcriptome of individual patients. The transcriptome is transformed into physical parameters that represent each sample's steady state and the free energy change (FEC) from that steady state, which is the state with the lowest free energy.We found the transcriptome steady state was invariant across normal and AML samples. FEC, representing active molecular processes, varied significantly between samples and was used to create patient-specific barcodes to characterize the biology of the disease. We discovered that AML samples that were in a transition state had the highest FEC. This disease state may be characterized as the most unstable and hence the most therapeutically targetable since a change in free energy is a thermodynamic requirement for disease progression. We also found that distinct sets of ongoing processes may be at the root of otherwise similar clinical phenotypes, implying that our integrated analysis of transcriptome profiles may facilitate a personalized medicine approach to cure AML and restore a steady state in each patient.

State-transition modeling of blood transcriptome predicts disease evolution and treatment response in chronic myeloid leukemia.

Chronic myeloid leukemia (CML) is initiated and maintained by BCR::ABL which is clinically targeted using tyrosine kinase inhibitors (TKIs). TKIs can induce long-term remission but are also not curative. Thus, CML is an ideal system to test our hypothesis that transcriptome-based state-transition models accurately predict cancer evolution and treatment response. We collected time-sequential blood samples from tetracycline-off (Tet-Off) BCR::ABL-inducible transgenic mice and wild-type controls. From the transcriptome, we constructed a CML state-space and a three-well leukemogenic potential landscape. The potential's stable critical points defined observable disease states. Early states were characterized by anti-CML genes opposing leukemia; late states were characterized by pro-CML genes. Genes with expression patterns shaped similarly to the potential landscape were identified as drivers of disease transition. Re-introduction of tetracycline to silence the BCR::ABL gene returned diseased mice transcriptomes to a near healthy state, without reaching it, suggesting parts of the transition are irreversible. TKI only reverted the transcriptome to an intermediate disease state, without approaching a state of health; disease relapse occurred soon after treatment. Using only the earliest time-point as initial conditions, our state-transition models accurately predicted both disease progression and treatment response, supporting this as a potentially valuable approach to time clinical intervention, before phenotypic changes become detectable.

Molecular Atlas of HER2+ Breast Cancer Cells Treated with Endogenous Ligands: Temporal Insights into Mechanisms of Trastuzumab Resistance.

Trastuzumab therapy in HER2+ breast cancer patients has mixed success owing to acquired resistance to therapy. A detailed understanding of downstream molecular cascades resulting from trastuzumab resistance is yet to emerge. In this study, we investigate the cellular mechanisms underlying acquired resistance using trastuzumab-sensitive and -resistant cancer cells (BT474 and BT474R) treated with endogenous ligands EGF and HRG across time. We probe early receptor organization through microscopy and signaling events through multiomics measurements and assess the bioenergetic state through mitochondrial measurements. Integrative analyses of our measurements reveal significant alterations in EGF-treated BT474 HER2 membrane dynamics and robust downstream activation of PI3K/AKT/mTORC1 signaling. EGF-treated BT474R shows a sustained interferon-independent activation of the IRF1/STAT1 cascade, potentially contributing to trastuzumab resistance. Both cell lines exhibit temporally divergent metabolic demands and HIF1A-mediated stress responses. BT474R demonstrates inherently increased mitochondrial activity. HRG treatment in BT474R leads to a pronounced reduction in AR expression, affecting downstream lipid metabolism with implications for treatment response. Our results provide novel insights into mechanistic changes underlying ligand treatment in BT474 and BT474R and emphasize the pivotal role of endogenous ligands. These results can serve as a framework for furthering the understanding of trastuzumab resistance, with therapeutic implications for women with acquired resistance.

Comparison of insertional RNA editing in Myxomycetes.

RNA editing describes the process in which individual or short stretches of nucleotides in a messenger or structural RNA are inserted, deleted, or substituted. A high level of RNA editing has been observed in the mitochondrial genome of Physarum polycephalum. The most frequent editing type in Physarum is the insertion of individual Cs. RNA editing is extremely accurate in Physarum; however, little is known about its mechanism. Here, we demonstrate how analyzing two organisms from the Myxomycetes, namely Physarum polycephalum and Didymium iridis, allows us to test hypotheses about the editing mechanism that can not be tested from a single organism alone. First, we show that using the recently determined full transcriptome information of Physarum dramatically improves the accuracy of computational editing site prediction in Didymium. We use this approach to predict genes in the mitochondrial genome of Didymium and identify six new edited genes as well as one new gene that appears unedited. Next we investigate sequence conservation in the vicinity of editing sites between the two organisms in order to identify sites that harbor the information for the location of editing sites based on increased conservation. Our results imply that the information contained within only nine or ten nucleotides on either side of the editing site (a distance previously suggested through experiments) is not enough to locate the editing sites. Finally, we show that the codon position bias in C insertional RNA editing of these two organisms is correlated with the selection pressure on the respective genes thereby directly testing an evolutionary theory on the origin of this codon bias. Beyond revealing interesting properties of insertional RNA editing in Myxomycetes, our work suggests possible approaches to be used when finding sequence motifs for any biological process fails.

State-transition Modeling of Blood Transcriptome Predicts Disease Evolution and Treatment Response in Chronic Myeloid Leukemia.

Chronic myeloid leukemia (CML) is initiated and maintained by BCR::ABL which is clinically targeted using tyrosine kinase inhibitors (TKIs). TKIs can induce long-term remission but are also not curative. Thus, CML is an ideal system to test our hypothesis that transcriptome-based state-transition models accurately predict cancer evolution and treatment response. We collected time-sequential blood samples from tetracycline-off (Tet-Off) BCR::ABL-inducible transgenic mice and wild-type controls. From the transcriptome, we constructed a CML state-space and a three-well leukemogenic potential landscape. The potential's stable critical points defined observable disease states. Early states were characterized by anti-CML genes opposing leukemia; late states were characterized by pro-CML genes. Genes with expression patterns shaped similarly to the potential landscape were identified as drivers of disease transition. Re-introduction of tetracycline to silence the BCR::ABL gene returned diseased mice transcriptomes to a near healthy state, without reaching it, suggesting parts of the transition are irreversible. TKI only reverted the transcriptome to an intermediate disease state, without approaching a state of health; disease relapse occurred soon after treatment. Using only the earliest time-point as initial conditions, our state-transition models accurately predicted both disease progression and treatment response, supporting this as a potentially valuable approach to time clinical intervention even before phenotypic changes become detectable.

Acquired miR-142 deficit in leukemic stem cells suffices to drive chronic myeloid leukemia into blast crisis.

The mechanisms underlying the transformation of chronic myeloid leukemia (CML) from chronic phase (CP) to blast crisis (BC) are not fully elucidated. Here, we show lower levels of miR-142 in CD34+CD38- blasts from BC CML patients than in those from CP CML patients, suggesting that miR-142 deficit is implicated in BC evolution. Thus, we create miR-142 knockout CML (i.e., miR-142-/-BCR-ABL) mice, which develop BC and die sooner than miR-142 wt CML (i.e., miR-142+/+BCR-ABL) mice, which instead remain in CP CML. Leukemic stem cells (LSCs) from miR-142-/-BCR-ABL mice recapitulate the BC phenotype in congenic recipients, supporting LSC transformation by miR-142 deficit. State-transition and mutual information analyses of "bulk" and single cell RNA-seq data, metabolomic profiling and functional metabolic assays identify enhanced fatty acid ß-oxidation, oxidative phosphorylation and mitochondrial fusion in LSCs as key steps in miR-142-driven BC evolution. A synthetic CpG-miR-142 mimic oligodeoxynucleotide rescues the BC phenotype in miR-142-/-BCR-ABL mice and patient-derived xenografts.

Enrichment-based DNA methylation analysis using next-generation sequencing: sample exclusion, estimating changes in global methylation, and the contribution of replicate lanes.

BACKGROUND: DNA methylation is an important epigenetic mark and dysregulation of DNA methylation is associated with many diseases including cancer. Advances in next-generation sequencing now allow unbiased methylome profiling of entire patient cohorts, greatly facilitating biomarker discovery and presenting new opportunities to understand the biological mechanisms by which changes in methylation contribute to disease. Enrichment-based sequencing assays such as MethylCap-seq are a cost effective solution for genome-wide determination of methylation status, but the technical reliability of methylation reconstruction from raw sequencing data has not been well characterized. METHODS: We analyze three MethylCap-seq data sets and perform two different analyses to assess data quality. First, we investigate how data quality is affected by excluding samples that do not meet quality control cutoff requirements. Second, we consider the effect of additional reads on enrichment score, saturation, and coverage. Lastly, we verify a method for the determination of the global amount of methylation from MethylCap-seq data by comparing to a spiked-in control DNA of known methylation status. RESULTS: We show that rejection of samples based on our quality control parameters leads to a significant improvement of methylation calling. Additional reads beyond ~13 million unique aligned reads improved coverage, modestly improved saturation, and did not impact enrichment score. Lastly, we find that a global methylation indicator calculated from MethylCap-seq data correlates well with the global methylation level of a sample as obtained from a spike-in DNA of known methylation level. CONCLUSIONS: We show that with appropriate quality control MethylCap-seq is a reliable tool, suitable for cohorts of hundreds of patients, that provides reproducible methylation information on a feature by feature basis as well as information about the global level of methylation.

Methods for high-throughput MethylCap-Seq data analysis.

BACKGROUND: Advances in whole genome profiling have revolutionized the cancer research field, but at the same time have raised new bioinformatics challenges. For next generation sequencing (NGS), these include data storage, computational costs, sequence processing and alignment, delineating appropriate statistical measures, and data visualization. Currently there is a lack of workflows for efficient analysis of large, MethylCap-seq datasets containing multiple sample groups. METHODS: The NGS application MethylCap-seq involves the in vitro capture of methylated DNA and subsequent analysis of enriched fragments by massively parallel sequencing. The workflow we describe performs MethylCap-seq experimental Quality Control (QC), sequence file processing and alignment, differential methylation analysis of multiple biological groups, hierarchical clustering, assessment of genome-wide methylation patterns, and preparation of files for data visualization. RESULTS: Here, we present a scalable, flexible workflow for MethylCap-seq QC, secondary data analysis, tertiary analysis of multiple experimental groups, and data visualization. We demonstrate the experimental QC procedure with results from a large ovarian cancer study dataset and propose parameters which can identify problematic experiments. Promoter methylation profiling and hierarchical clustering analyses are demonstrated for four groups of acute myeloid leukemia (AML) patients. We propose a Global Methylation Indicator (GMI) function to assess genome-wide changes in methylation patterns between experimental groups. We also show how the workflow facilitates data visualization in a web browser with the application Anno-J. CONCLUSIONS: This workflow and its suite of features will assist biologists in conducting methylation profiling projects and facilitate meaningful biological interpretation.

Dynamic patterns of microRNA expression during acute myeloid leukemia state-transition.

MicroRNAs (miRNAs) have been shown to hold prognostic value in acute myeloid leukemia (AML); however, the temporal dynamics of miRNA expression in AML are poorly understood. Using serial samples from a mouse model of AML to generate time-series miRNA sequencing data, we are the first to show that the miRNA transcriptome undergoes state-transition during AML initiation and progression. We modeled AML state-transition as a particle undergoing Brownian motion in a quasi-potential and validated the AML state-space and state-transition model to accurately predict time to AML in an independent cohort of mice. The critical points of the model provided a framework to align samples from mice that developed AML at different rates. Our mathematical approach allowed discovery of dynamic processes involved during AML development and, if translated to humans, has the potential to predict an individual's disease trajectory.

Dietary omega-3 fatty acid intake impacts peripheral blood DNA methylation -anti-inflammatory effects and individual variability in a pilot study.

Omega-3 or n-3 polyunsaturated fatty acids (PUFAs) are widely studied for health benefits that may relate to anti-inflammatory activity. However, mechanisms mediating an anti-inflammatory response to n-3 PUFA intake are not fully understood. Of interest is the emerging role of fatty acids to impact DNA methylation (DNAm) and thereby modulate mediating inflammatory processes. In this pilot study, we investigated the impact of n-3 PUFA intake on DNAm in inflammation-related signaling pathways in peripheral blood mononuclear cells (PBMCs) of women at high risk of breast cancer. PBMCs of women at high risk of breast cancer (n=10) were obtained at baseline and after 6 months of n-3 PUFA (5 g/d EPA+DHA dose arm) intake in a previously reported dose finding trial. DNA methylation of PBMCs was assayed by reduced representation bisulfite sequencing (RRBS) to obtain genome-wide methylation profiles at the single nucleotide level. We examined the impact of n-3 PUFA on genome-wide DNAm and focused upon a set of candidate genes associated with inflammation signaling pathways and breast cancer. We identified 24,842 differentially methylated CpGs (DMCs) in gene promoters of 5507 genes showing significant enrichment for hypermethylation in both the candidate gene and genome-wide analyses. Pathway analysis identified significantly hypermethylated signaling networks after n-3 PUFA treatment, such as the Toll-like Receptor inflammatory pathway. The DNAm pattern in individuals and the response to n-3 PUFA intake are heterogeneous. PBMC DNAm profiling suggests a mechanism whereby n-3 PUFAs may impact inflammatory cascades associated with disease processes including carcinogenesis.

Regulation of chromatin accessibility by the histone chaperone CAF-1 sustains lineage fidelity.

Cell fate commitment is driven by dynamic changes in chromatin architecture and activity of lineage-specific transcription factors (TFs). The chromatin assembly factor-1 (CAF-1) is a histone chaperone that regulates chromatin architecture by facilitating nucleosome assembly during DNA replication. Accumulating evidence supports a substantial role of CAF-1 in cell fate maintenance, but the mechanisms by which CAF-1 restricts lineage choice remain poorly understood. Here, we investigate how CAF-1 influences chromatin dynamics and TF activity during lineage differentiation. We show that CAF-1 suppression triggers rapid differentiation of myeloid stem and progenitor cells into a mixed lineage state. We find that CAF-1 sustains lineage fidelity by controlling chromatin accessibility at specific loci, and limiting the binding of ELF1 TF at newly-accessible diverging regulatory elements. Together, our findings decipher key traits of chromatin accessibility that sustain lineage integrity and point to a powerful strategy for dissecting transcriptional circuits central to cell fate commitment.

Vascularity and Dynamic Contrast-Enhanced Breast Magnetic Resonance Imaging.

Angiogenesis is a key step in the initiation and progression of an invasive breast cancer. High microvessel density by morphological characterization predicts metastasis and poor survival in women with invasive breast cancers. However, morphologic characterization is subject to variability and only can evaluate a limited portion of an invasive breast cancer. Consequently, breast Magnetic Resonance Imaging (MRI) is currently being evaluated to assess vascularity. Recently, through the new field of radiomics, dynamic contrast enhanced (DCE)-MRI is being used to evaluate vascular density, vascular morphology, and detection of aggressive breast cancer biology. While DCE-MRI is a highly sensitive tool, there are specific features that limit computational evaluation of blood vessels. These include (1) DCE-MRI evaluates gadolinium contrast and does not directly evaluate biology, (2) the resolution of DCE-MRI is insufficient for imaging small blood vessels, and (3) DCE-MRI images are very difficult to co-register. Here we review computational approaches for detection and analysis of blood vessels in DCE-MRI images and present some of the strategies we have developed for co-registry of DCE-MRI images and early detection of vascularization.

Treatment-induced arteriolar revascularization and miR-126 enhancement in bone marrow niche protect leukemic stem cells in AML.

BACKGROUND: During acute myeloid leukemia (AML) growth, the bone marrow (BM) niche acquires significant vascular changes that can be offset by therapeutic blast cytoreduction. The molecular mechanisms of this vascular plasticity remain to be fully elucidated. Herein, we report on the changes that occur in the vascular compartment of the FLT3-ITD+ AML BM niche pre and post treatment and their impact on leukemic stem cells (LSCs). METHODS: BM vasculature was evaluated in FLT3-ITD+ AML models (MllPTD/WT/Flt3ITD/ITD mouse and patient-derived xenograft) by 3D confocal imaging of long bones, calvarium vascular permeability assays, and flow cytometry analysis. Cytokine levels were measured by Luminex assay and miR-126 levels evaluated by Q-RT-PCR and miRNA staining. Wild-type (wt) and MllPTD/WT/Flt3ITD/ITD mice with endothelial cell (EC) miR-126 knockout or overexpression served as controls. The impact of treatment-induced BM vascular changes on LSC activity was evaluated by secondary transplantation of BM cells after administration of tyrosine kinase inhibitors (TKIs) to MllPTD/WT/Flt3ITD/ITD mice with/without either EC miR-126 KO or co-treatment with tumor necrosis factor alpha (TNFα) or anti-miR-126 miRisten. RESULTS: In the normal BM niche, CD31+Sca-1high ECs lining arterioles have miR-126 levels higher than CD31+Sca-1low ECs lining sinusoids. We noted that during FLT3-ITD+ AML growth, the BM niche lost arterioles and gained sinusoids. These changes were mediated by TNFα, a cytokine produced by AML blasts, which induced EC miR-126 downregulation and caused depletion of CD31+Sca-1high ECs and gain in CD31+Sca-1low ECs. Loss of miR-126high ECs led to a decreased EC miR-126 supply to LSCs, which then entered the cell cycle and promoted leukemia growth. Accordingly, antileukemic treatment with TKI decreased the BM blast-produced TNFα and increased miR-126high ECs and the EC miR-126 supply to LSCs. High miR-126 levels safeguarded LSCs, as shown by more severe disease in secondary transplanted mice. Conversely, EC miR-126 deprivation via genetic or pharmacological EC miR-126 knock-down prevented treatment-induced BM miR-126high EC expansion and in turn LSC protection. CONCLUSIONS: Treatment-induced CD31+Sca-1high EC re-vascularization of the leukemic BM niche may represent a LSC extrinsic mechanism of treatment resistance that can be overcome with therapeutic EC miR-126 deprivation.

Mathematical deconvolution of CAR T-cell proliferation and exhaustion from real-time killing assay data.

Chimeric antigen receptor (CAR) T-cell therapy has shown promise in the treatment of haematological cancers and is currently being investigated for solid tumours, including high-grade glioma brain tumours. There is a desperate need to quantitatively study the factors that contribute to the efficacy of CAR T-cell therapy in solid tumours. In this work, we use a mathematical model of predator-prey dynamics to explore the kinetics of CAR T-cell killing in glioma: the Chimeric Antigen Receptor T-cell treatment Response in GliOma (CARRGO) model. The model includes rates of cancer cell proliferation, CAR T-cell killing, proliferation, exhaustion, and persistence. We use patient-derived and engineered cancer cell lines with an in vitro real-time cell analyser to parametrize the CARRGO model. We observe that CAR T-cell dose correlates inversely with the killing rate and correlates directly with the net rate of proliferation and exhaustion. This suggests that at a lower dose of CAR T-cells, individual T-cells kill more cancer cells but become more exhausted when compared with higher doses. Furthermore, the exhaustion rate was observed to increase significantly with tumour growth rate and was dependent on level of antigen expression. The CARRGO model highlights nonlinear dynamics involved in CAR T-cell therapy and provides novel insights into the kinetics of CAR T-cell killing. The model suggests that CAR T-cell treatment may be tailored to individual tumour characteristics including tumour growth rate and antigen level to maximize therapeutic benefit.