Reconstruction of single-cell lineage trajectories and identification of diversity in fates during the epithelial-to-mesenchymal transition.

Exploring the complexity of the epithelial-to-mesenchymal transition (EMT) unveils a diversity of potential cell fates; however, the exact timing and mechanisms by which early cell states diverge into distinct EMT trajectories remain unclear. Studying these EMT trajectories through single-cell RNA sequencing is challenging due to the necessity of sacrificing cells for each measurement. In this study, we employed optimal-transport analysis to reconstruct the past trajectories of different cell fates during TGF-beta-induced EMT in the MCF10A cell line.

Mathematical Modeling Identifies Optimum Palbociclib-fulvestrant Dose Administration Schedules for the Treatment of Patients with Estrogen Receptor-positive Breast Cancer.

Unlabelled: Cyclin-dependent kinases 4/6 (CDK4/6) inhibitors such as palbociclib are approved for the treatment of metastatic estrogen receptor-positive (ER+) breast cancer in combination with endocrine therapies and significantly improve outcomes in patients with this disease. However, given the large number of possible pairwise drug combinations and administration schedules, it remains unclear which clinical strategy would lead to best survival. Here, we developed a computational, cell cycle-explicit model to characterize the pharmacodynamic response to palbociclib-fulvestrant combination therapy.

Single-cell heterogeneity of EGFR and CDK4 co-amplification is linked to immune infiltration in glioblastoma.

Glioblastoma (GBM) is the most aggressive brain tumor, with a median survival of ∼15 months. Targeted approaches have not been successful in this tumor type due to the large extent of intratumor heterogeneity. Mosaic amplification of oncogenes suggests that multiple genetically distinct clones are present in each tumor.

Identification of optimal dosing schedules of dacomitinib and osimertinib for a phase I/II trial in advanced EGFR-mutant non-small cell lung cancer.

Despite the clinical success of the third-generation EGFR inhibitor osimertinib as a first-line treatment of EGFR-mutant non-small cell lung cancer (NSCLC), resistance arises due to the acquisition of EGFR second-site mutations and other mechanisms, which necessitates alternative therapies. Dacomitinib, a pan-HER inhibitor, is approved for first-line treatment and results in different acquired EGFR mutations than osimertinib that mediate on-target resistance. A combination of osimertinib and dacomitinib could therefore induce more durable responses by preventing the emergence of resistance.

Breast tumours maintain a reservoir of subclonal diversity during expansion.

Our knowledge of copy number evolution during the expansion of primary breast tumours is limited. Here, to investigate this process, we developed a single-cell, single-molecule DNA-sequencing method and performed copy number analysis of 16,178 single cells from 8 human triple-negative breast cancers and 4 cell lines. The results show that breast tumours and cell lines comprise a large milieu of subclones (7-22) that are organized into a few (3-5) major superclones.

Imaging dynamic mTORC1 pathway activity in vivo reveals marked shifts that support time-specific inhibitor therapy in AML.

Acute myeloid leukemia (AML) is a high remission, high relapse fatal blood cancer. Although mTORC1 is a master regulator of cell proliferation and survival, its inhibitors have not performed well as AML treatments. To uncover the dynamics of mTORC1 activity in vivo, fluorescent probes are developed to track single cell proliferation, apoptosis and mTORC1 activity of AML cells in the bone marrow of live animals and to quantify these activities in the context of microanatomical localization and intra-tumoral heterogeneity.

Beyond Deterministic Models in Drug Discovery and Development.

The model-informed drug discovery and development paradigm is now well established among the pharmaceutical industry and regulatory agencies. This success has been mainly due to the ability of pharmacometrics to bring together different modeling strategies, such as population pharmacokinetics/pharmacodynamics (PK/PD) and systems biology/pharmacology. However, there are promising quantitative approaches that are still seldom used by pharmacometricians and that deserve consideration.

Pharmacokinetic Profiles Determine Optimal Combination Treatment Schedules in Computational Models of Drug Resistance.

Identification of optimal schedules for combination drug administration relies on accurately estimating the correct pharmacokinetics, pharmacodynamics, and drug interaction effects. Misspecification of pharmacokinetics can lead to wrongly predicted timing or order of treatments, leading to schedules recommended on the basis of incorrect assumptions about absorption and elimination of a drug and its effect on tumor growth. Here, we developed a computational modeling platform and software package for combination treatment strategies with flexible pharmacokinetic profiles and multidrug interaction curves that are estimated from data.

ESTIpop: a computational tool to simulate and estimate parameters for continuous-time Markov branching processes.

Summary: ESTIpop is an R package designed to simulate and estimate parameters for continuous-time Markov branching processes with constant or time-dependent rates, a common model for asexually reproducing cell populations. Analytical approaches to parameter estimation quickly become intractable in complex branching processes. In ESTIpop, parameter estimation is based on a likelihood function with respect to a time series of cell counts, approximated by the Central Limit Theorem for multitype branching processes.

Synthetic Lethal and Resistance Interactions with BET Bromodomain Inhibitors in Triple-Negative Breast Cancer.

BET bromodomain inhibitors (BBDIs) are candidate therapeutic agents for triple-negative breast cancer (TNBC) and other cancer types, but inherent and acquired resistance to BBDIs limits their potential clinical use. Using CRISPR and small-molecule inhibitor screens combined with comprehensive molecular profiling of BBDI response and resistance, we identified synthetic lethal interactions with BBDIs and genes that, when deleted, confer resistance. We observed synergy with regulators of cell cycle progression, YAP, AXL, and SRC signaling, and chemotherapeutic agents.

Acquired resistance to combined BET and CDK4/6 inhibition in triple-negative breast cancer.

BET inhibitors are promising therapeutic agents for the treatment of triple-negative breast cancer (TNBC), but the rapid emergence of resistance necessitates investigation of combination therapies and their effects on tumor evolution. Here, we show that palbociclib, a CDK4/6 inhibitor, and paclitaxel, a microtubule inhibitor, synergize with the BET inhibitor JQ1 in TNBC lines. High-complexity DNA barcoding and mathematical modeling indicate a high rate of de novo acquired resistance to these drugs relative to pre-existing resistance.

DIFFpop: a stochastic computational approach to simulate differentiation hierarchies with single cell barcoding.

Summary: DIFFpop is an R package designed to simulate cellular differentiation hierarchies using either exponentially-expanding or fixed population sizes. The software includes functionalities to simulate clonal evolution due to the emergence of driver mutations under the infinite-allele assumption as well as options for simulation and analysis of single cell barcoding and labeling data. The software uses the Gillespie Stochastic Simulation Algorithm and a modification of expanding or fixed-size stochastic process models expanded to a large number of cell types and scenarios.

KDM5 Histone Demethylase Activity Links Cellular Transcriptomic Heterogeneity to Therapeutic Resistance.

Members of the KDM5 histone H3 lysine 4 demethylase family are associated with therapeutic resistance, including endocrine resistance in breast cancer, but the underlying mechanism is poorly defined. Here we show that genetic deletion of KDM5A/B or inhibition of KDM5 activity increases sensitivity to anti-estrogens by modulating estrogen receptor (ER) signaling and by decreasing cellular transcriptomic heterogeneity. Higher KDM5B expression levels are associated with higher transcriptomic heterogeneity and poor prognosis in ER breast tumors.

SIApopr: a computational method to simulate evolutionary branching trees for analysis of tumor clonal evolution.

Summary: SIApopr (Simulating Infinite-Allele populations) is an R package to simulate time-homogeneous and inhomogeneous stochastic branching processes under a very flexible set of assumptions using the speed of C ++. The software simulates clonal evolution with the emergence of driver and passenger mutations under the infinite-allele assumption. The software is an application of the Gillespie Stochastic Simulation Algorithm expanded to a large number of cell types and scenarios, with the intention of allowing users to easily modify existing models or create their own.

Punctuated copy number evolution and clonal stasis in triple-negative breast cancer.

Aneuploidy is a hallmark of breast cancer; however, knowledge of how these complex genomic rearrangements evolve during tumorigenesis is limited. In this study, we developed a highly multiplexed single-nucleus sequencing method to investigate copy number evolution in patients with triple-negative breast cancer. We sequenced 1,000 single cells from tumors in 12 patients and identified 1-3 major clonal subpopulations in each tumor that shared a common evolutionary lineage.

Evolutionary Action Score of TP53 Coding Variants Is Predictive of Platinum Response in Head and Neck Cancer Patients.

TP53 is the most frequently altered gene in head and neck squamous cell carcinoma (HNSCC), with mutations occurring in over two thirds of cases; however, the predictive response of these mutations to cisplatin-based therapy remains elusive. In the current study, we evaluate the ability of the Evolutionary Action score of TP53-coding variants (EAp53) to predict the impact of TP53 mutations on response to chemotherapy. The EAp53 approach clearly identifies a subset of high-risk TP53 mutations associated with decreased sensitivity to cisplatin both in vitro and in vivo in preclinical models of HNSCC.

Wee-1 kinase inhibition overcomes cisplatin resistance associated with high-risk TP53 mutations in head and neck cancer through mitotic arrest followed by senescence.

Although cisplatin has played a role in "standard-of-care" multimodality therapy for patients with advanced squamous cell carcinoma of the head and neck (HNSCC), the rate of treatment failure remains particularly high for patients receiving cisplatin whose tumors have mutations in the TP53 gene. We found that cisplatin treatment of HNSCC cells with mutant TP53 leads to arrest of cells in the G2 phase of the cell cycle, leading us to hypothesize that the wee-1 kinase inhibitor MK-1775 would abrogate the cisplatin-induced G2 block and thereby sensitize isogenic HNSCC cells with mutant TP53 or lacking p53 expression to cisplatin. We tested this hypothesis using clonogenic survival assays, flow cytometry, and in vivo tumor growth delay experiments with an orthotopic nude mouse model of oral tongue cancer.

The neuroimmune guidance cue netrin-1 promotes atherosclerosis by inhibiting the emigration of macrophages from plaques.

Atherosclerotic plaque formation is fueled by the persistence of lipid-laden macrophages in the artery wall. The mechanisms by which these cells become trapped, thereby establishing chronic inflammation, remain unknown. Here we found that netrin-1, a neuroimmune guidance cue, was secreted by macrophages in human and mouse atheroma, where it inactivated the migration of macrophages toward chemokines linked to their egress from plaques.

CREB downregulation in vascular disease: a common response to cardiovascular risk.

Objective: To examine the impact of low-density lipoprotein (LDL), an established mediator of atherosclerosis, on the transcription factor cAMP-response element-binding protein (CREB), which is a regulator of vascular smooth muscle cell (VSMC) quiescence.

Methods And Results: VSMC CREB content is diminished in rodent models of diabetes and pulmonary hypertension. We examined aortic CREB content in rodent models of aging, hypertension, and insulin resistance, and we determined nuclear CREB protein in the medial VSMC of high-fat-fed LDL receptor-null mice.

Mice lacking hepatic lipase are lean and protected against diet-induced obesity and hepatic steatosis.

Hepatic lipase (HL)-mediated lipoprotein hydrolysis provides free fatty acids for energy, storage, and nutrient signaling and may play a role in energy homeostasis. Because HL-activity increases with increased visceral fat, we hypothesized that increased HL-activity favors weight gain and obesity and consequently, that HL deficiency would reduce body fat stores and protect against diet-induced obesity. To test this hypothesis, we compared wild-type mice (with endogenous HL) and mice genetically deficient in HL with respect to daily body weight and food intake, body composition, and adipocyte size on both chow and high-fat (HF) diets.

A stochastic model of cell cycle desynchronization.

A general branching process model is suggested to describe cell cycle desynchronization. Cell cycle phase times are modeled as random variables and a formula for the expected fraction of cells in S phase as a function of time is established. The model is compared to data from the literature and is also compared to previously suggested deterministic and stochastic models.

Serum amyloid A3 does not contribute to circulating SAA levels.

Adipose tissue secretes proteins like serum amyloid A (SAA), which plays important roles in local and systemic inflammation. Circulating SAA levels increase in obese humans, but the roles of adipose-derived SAA and hyperlipidemia in this process are unclear. We took advantage of the difference in the inducible isoforms of SAA secreted by adipose tissue (SAA3) and liver (SAA1 and 2) of mice to evaluate whether adipose tissue contributes to the circulating pool of SAA in obesity and hyperlipidemia.