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

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.

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 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.

The immunogenic radiation and new players in immunotherapy and targeted therapy for head and neck cancer

Although treatment modalities for head and neck cancer have evolved considerably over the past decades, survival rates have plateaued. The treatment options remained limited to definitive surgery, surgery followed by fractionated radiotherapy with optional chemotherapy, and a definitive combination of fractionated radiotherapy and chemotherapy. Lately, immunotherapy has been introduced as the fourth modality of treatment, mainly administered as a single checkpoint inhibitor for recurrent or metastatic disease. While other regimens and combinations of immunotherapy and targeted therapy are being tested in clinical trials, adapting the appropriate regimens to patients and predicting their outcomes have yet to reach the clinical setting. Radiotherapy is mainly regarded as a means to target cancer cells while minimizing the unwanted peripheral effect. Radiotherapy regimens and fractionation are designed to serve this purpose, while the systemic effect of radiation on the immune response is rarely considered a factor while designing treatment. To bridge this gap, this review will highlight the effect of radiotherapy on the tumor microenvironment locally, and the immune response systemically. We will review the methodology to identify potential targets for therapy in the tumor microenvironment and the scientific basis for combining targeted therapy and radiotherapy. We will describe a current experience in preclinical models to test these combinations and propose how challenges in this realm may be faced. We will review new players in targeted therapy and their utilization to drive immunogenic response against head and neck cancer. We will outline the factors contributing to head and neck cancer heterogeneity and their effect on the response to radiotherapy. We will review in-silico methods to decipher intertumoral and intratumoral heterogeneity and how these algorithms can predict treatment outcomes. We propose that (a) the sequence of surgery, radiotherapy, chemotherapy, and targeted therapy should be designed not only to annul cancer directly, but to prime the immune response. (b) Fractionation of radiotherapy and the extent of the irradiated field should facilitate systemic immunity to develop. (c) New players in targeted therapy should be evaluated in translational studies toward clinical trials. (d) Head and neck cancer treatment should be personalized according to patients and tumor-specific factors.

Data driven model discovery and interpretation for CAR T-cell killing using sparse identification and latent variables

In the development of cell-based cancer therapies, quantitative mathematical models of cellular interactions are instrumental in understanding treatment efficacy. Efforts to validate and interpret mathematical models of cancer cell growth and death hinge first on proposing a precise mathematical model, then analyzing experimental data in the context of the chosen model. In this work, we present the first application of the sparse identification of non-linear dynamics (SINDy) algorithm to a real biological system in order discover cell-cell interaction dynamics in in vitro experimental data, using chimeric antigen receptor (CAR) T-cells and patient-derived glioblastoma cells. By combining the techniques of latent variable analysis and SINDy, we infer key aspects of the interaction dynamics of CAR T-cell populations and cancer. Importantly, we show how the model terms can be interpreted biologically in relation to different CAR T-cell functional responses, single or double CAR T-cell-cancer cell binding models, and density-dependent growth dynamics in either of the CAR T-cell or cancer cell populations. We show how this data-driven model-discovery based approach provides unique insight into CAR T-cell dynamics when compared to an established model-first approach. These results demonstrate the potential for SINDy to improve the implementation and efficacy of CAR T-cell therapy in the clinic through an improved understanding of CAR T-cell dynamics.