Spatially confined sub-tumor microenvironments in pancreatic cancer.

Intratumoral heterogeneity is a critical frontier in understanding how the tumor microenvironment (TME) propels malignant progression. Here, we deconvolute the human pancreatic TME through large-scale integration of histology-guided regional multiOMICs with clinical data and patient-derived preclinical models. We discover "subTMEs," histologically definable tissue states anchored in fibroblast plasticity, with regional relationships to tumor immunity, subtypes, differentiation, and treatment response. "Reactive" subTMEs rich in complex but functionally coordinated fibroblast communities were immune hot and inhabited by aggressive tumor cell phenotypes. The matrix-rich "deserted" subTMEs harbored fewer activated fibroblasts and tumor-suppressive features yet were markedly chemoprotective and enriched upon chemotherapy. SubTMEs originated in fibroblast differentiation trajectories, and transitory states were notable both in single-cell transcriptomics and in situ. The intratumoral co-occurrence of subTMEs produced patient-specific phenotypic and computationally predictable heterogeneity tightly linked to malignant biology. Therefore, heterogeneity within the plentiful, notorious pancreatic TME is not random but marks fundamental tissue organizational units.

Ordered and deterministic cancer genome evolution after p53 loss.

Although p53 inactivation promotes genomic instability1 and presents a route to malignancy for more than half of all human cancers2,3, the patterns through which heterogenous TP53 (encoding human p53) mutant genomes emerge and influence tumorigenesis remain poorly understood. Here, in a mouse model of pancreatic ductal adenocarcinoma that reports sporadic p53 loss of heterozygosity before cancer onset, we find that malignant properties enabled by p53 inactivation are acquired through a predictable pattern of genome evolution. Single-cell sequencing and in situ genotyping of cells from the point of p53 inactivation through progression to frank cancer reveal that this deterministic behaviour involves four sequential phases-Trp53 (encoding mouse p53) loss of heterozygosity, accumulation of deletions, genome doubling, and the emergence of gains and amplifications-each associated with specific histological stages across the premalignant and malignant spectrum. Despite rampant heterogeneity, the deletion events that follow p53 inactivation target functionally relevant pathways that can shape genomic evolution and remain fixed as homogenous events in diverse malignant populations. Thus, loss of p53-the 'guardian of the genome'-is not merely a gateway to genetic chaos but, rather, can enable deterministic patterns of genome evolution that may point to new strategies for the treatment of TP53-mutant tumours.

Childhood cerebellar tumours mirror conserved fetal transcriptional programs.

Study of the origin and development of cerebellar tumours has been hampered by the complexity and heterogeneity of cerebellar cells that change over the course of development. Here we use single-cell transcriptomics to study more than 60,000 cells from the developing mouse cerebellum and show that different molecular subgroups of childhood cerebellar tumours mirror the transcription of cells from distinct, temporally restricted cerebellar lineages. The Sonic Hedgehog medulloblastoma subgroup transcriptionally mirrors the granule cell hierarchy as expected, while group 3 medulloblastoma resembles Nestin+ stem cells, group 4 medulloblastoma resembles unipolar brush cells, and PFA/PFB ependymoma and cerebellar pilocytic astrocytoma resemble the prenatal gliogenic progenitor cells. Furthermore, single-cell transcriptomics of human childhood cerebellar tumours demonstrates that many bulk tumours contain a mixed population of cells with divergent differentiation. Our data highlight cerebellar tumours as a disorder of early brain development and provide a proximate explanation for the peak incidence of cerebellar tumours in early childhood.

A TEAD2-Driven Endothelial-Like Program Shapes Basal-Like Differentiation and Metastasis of Pancreatic Cancer.

BACKGROUND & AIMS: Pancreatic ductal adenocarcinoma (PDA), with its highly metastatic propensity, is one of the most lethal subtypes of pancreatic cancer. Although recent large-scale transcriptomic studies have demonstrated that heterogeneous gene expressions play an essential role in determining molecular phenotypes of PDA, biological cues for and consequences of distinct transcriptional programs remain unclear. METHODS: We developed an experimental model that enforces the transition of PDA cells toward a basal-like subtype. We combined epigenome and transcriptome analyses with extensive in vitro and in vivo evaluations of tumorigenicity to demonstrate the validity of basal-like subtype differentiation in association with endothelial-like enhancer landscapes via TEA domain transcription factor 2 (TEAD2). Finally, we used loss-of-function experiments to investigate the importance of TEAD2 in regulating reprogrammed enhancer landscape and metastasis in basal-like PDA cells. RESULTS: Aggressive characteristics of the basal-like subtype are faithfully recapitulated in vitro and in vivo, demonstrating the physiological relevance of our model. Further, we showed that basal-like subtype PDA cells acquire a TEAD2-dependent proangiogenic enhancer landscape. Genetic and pharmacologic inhibitions of TEAD2 in basal-like subtype PDA cells impair their proangiogenic phenotypes in vitro and cancer progression in vivo. Last, we identify CD109 as a critical TEAD2 downstream mediator that maintains constitutively activated JAK-STAT signaling in basal-like PDA cells and tumors. CONCLUSIONS: Our findings implicate a TEAD2-CD109-JAK/STAT axis in the basal-like differentiated pancreatic cancer cells and as a potential therapeutic vulnerability.

GATA4 and GATA6 loss-of-expression is associated with extinction of the classical programme and poor outcome in pancreatic ductal adenocarcinoma.

OBJECTIVE: GATA6 is a key regulator of the classical phenotype in pancreatic ductal adenocarcinoma (PDAC). Low GATA6 expression associates with poor patient outcome. GATA4 is the second most expressed GATA factor in the pancreas. We assessed whether, and how, GATA4 contributes to PDAC phenotype and analysed the association of expression with outcome and response to chemotherapy. DESIGN: We analysed PDAC transcriptomic data, stratifying cases according to GATA4 and GATA6 expression and identified differentially expressed genes and pathways. The genome-wide distribution of GATA4 was assessed, as well as the effects of GATA4 knockdown. A multicentre tissue microarray study to assess GATA4 and GATA6 expression in samples (n=745) from patients with resectable was performed. GATA4 and GATA6 levels were dichotomised into high/low categorical variables; association with outcome was assessed using univariable and multivariable Cox regression models. RESULTS: GATA4 messenger RNA is enriched in classical, compared with basal-like tumours. We classified samples in 4 groups as high/low for GATA4 and GATA6. Reduced expression of GATA4 had a minor transcriptional impact but low expression of GATA4 enhanced the effects of GATA6 low expression. GATA4 and GATA6 display a partially overlapping genome-wide distribution, mainly at promoters. Reduced expression of both proteins in tumours was associated with the worst patient survival. GATA4 and GATA6 expression significantly decreased in metastases and negatively correlated with basal markers. CONCLUSIONS: GATA4 and GATA6 cooperate to maintain the classical phenotype. Our findings provide compelling rationale to assess their expression as biomarkers of poor prognosis and therapeutic response.

Systemic inflammatory prognostic scores in advanced pancreatic adenocarcinoma.

BACKGROUND: Systemic inflammatory scores may aid prognostication and patient selection for trials. We compared five scores in advanced pancreatic adenocarcinoma (PDAC). METHODS: Unresectable/metastatic PDAC patients enrolled in the Comprehensive Molecular Characterisation of Advanced Pancreatic Ductal Adenocarcinoma for Better Treatment Selection trial (NCT02750657) were included. Patients had pre-treatment biopsies for whole genome and RNA sequencing. CD8 immunohistochemistry was available in a subset. The neutrophil-to-lymphocyte ratio (NLR), platelet-to-lymphocyte ratio, Prognostic Nutritional Index, Gustave Roussy Immune Score (GRIm-S), and Memorial Sloan Kettering Prognostic Score (MPS) were calculated. Overall survival (OS) was estimated using Kaplan-Meier methods. Associations between inflammatory scores, clinical/genomic characteristics, and OS were analysed. RESULTS: We analysed 263 patients. High-risk NLR, GRIm-S and MPS were poorly prognostic. The GRIm-S had the highest predictive ability: median OS 6.4 vs. 10 months for high risk vs. low-risk (P < 0.001); HR 2.26 (P < 0.001). ECOG ≥ 1, the basal-like subtype, and low-HRDetect were additional poor prognostic factors (P < 0.01). Inflammatory scores did not associate with RNA-based classifiers or homologous recombination repair deficiency genotypes. High-risk MPS (P = 0.04) and GRIm-S (P = 0.02) patients had lower median CD8 + tumour-infiltrating lymphocytes. CONCLUSIONS: Inflammatory scores incorporating NLR have prognostic value in advanced PDAC. Understanding immunophenotypes of poor-risk patients and using these scores in trials will advance the field.

Machine Learning-Based Early Warning Systems for Acute Care Utilization During Systemic Therapy for Cancer.

BACKGROUND: Emergency department visits and hospitalizations frequently occur during systemic therapy for cancer. We developed and evaluated a longitudinal warning system for acute care use. METHODS: Using a retrospective population-based cohort of patients who started intravenous systemic therapy for nonhematologic cancers between July 1, 2014, and June 30, 2020, we randomly separated patients into cohorts for model training, hyperparameter tuning and model selection, and system testing. Predictive features included static features, such as demographics, cancer type, and treatment regimens, and dynamic features, such as patient-reported symptoms and laboratory values. The longitudinal warning system predicted the probability of acute care utilization within 30 days after each treatment session. Machine learning systems were developed in the training and tuning cohorts and evaluated in the testing cohort. Sensitivity analyses considered feature importance, other acute care endpoints, and performance within subgroups. RESULTS: The cohort included 105,129 patients who received 1,216,385 treatment sessions. Acute care followed 182,444 (15.0%) treatments within 30 days. The ensemble model achieved an area under the receiver operating characteristic curve of 0.742 (95% CI, 0.739-0.745) and was well calibrated in the test cohort. Important predictive features included prior acute care use, treatment regimen, and laboratory tests. If the system was set to alarm approximately once every 15 treatments, 25.5% of acute care events would be preceded by an alarm, and 47.4% of patients would experience acute care after an alarm. The system underestimated risk for some treatment regimens and potentially underserved populations such as females and non-English speakers. CONCLUSIONS: Machine learning warning systems can detect patients at risk for acute care utilization, which can aid in preventive intervention and facilitate tailored treatment. Future research should address potential biases and prospectively evaluate impact after system deployment.

Genomic Features and Classification of Homologous Recombination Deficient Pancreatic Ductal Adenocarcinoma.

BACKGROUND AND AIMS: Homologous recombination deficiency (HRD) in pancreatic ductal adenocarcinoma (PDAC), remains poorly defined beyond germline (g) alterations in BRCA1, BRCA2, and PALB2. METHODS: We interrogated whole genome sequencing (WGS) data on 391 patients, including 49 carriers of pathogenic variants (PVs) in gBRCA and PALB2. HRD classifiers were applied to the dataset and included (1) the genomic instability score (GIS) used by Myriad's MyChoice HRD assay; (2) substitution base signature 3 (SBS3); (3) HRDetect; and (4) structural variant (SV) burden. Clinical outcomes and responses to chemotherapy were correlated with HRD status. RESULTS: Biallelic tumor inactivation of gBRCA or PALB2 was evident in 43 of 49 germline carriers identifying HRD-PDAC. HRDetect (score ≥0.7) predicted gBRCA1/PALB2 deficiency with highest sensitivity (98%) and specificity (100%). HRD genomic tumor classifiers suggested that 7% to 10% of PDACs that do not harbor gBRCA/PALB2 have features of HRD. Of the somatic HRDetecthi cases, 69% were attributed to alterations in BRCA1/2, PALB2, RAD51C/D, and XRCC2, and a tandem duplicator phenotype. TP53 loss was more common in BRCA1- compared with BRCA2-associated HRD-PDAC. HRD status was not prognostic in resected PDAC; however in advanced disease the GIS (P = .02), SBS3 (P = .03), and HRDetect score (P = .005) were predictive of platinum response and superior survival. PVs in gATM (n = 6) or gCHEK2 (n = 2) did not result in HRD-PDAC by any of the classifiers. In 4 patients, BRCA2 reversion mutations associated with platinum resistance. CONCLUSIONS: Germline and parallel somatic profiling of PDAC outperforms germline testing alone in identifying HRD-PDAC. An additional 7% to 10% of patients without gBRCA/PALB2 mutations may benefit from DNA damage response agents.

Revised map of the human progenitor hierarchy shows the origin of macrophages and dendritic cells in early lymphoid development.

The classical model of hematopoiesis posits the segregation of lymphoid and myeloid lineages as the earliest fate decision. The validity of this model in the mouse has been questioned; however, little is known about the lineage potential of human progenitors. Here we provide a comprehensive analysis of the human hematopoietic hierarchy by clonally mapping the developmental potential of seven progenitor classes from neonatal cord blood and adult bone marrow. Human multilymphoid progenitors, identified as a distinct population of Thy-1(neg-lo)CD45RA(+) cells in the CD34(+)CD38(-) stem cell compartment, gave rise to all lymphoid cell types, as well as monocytes, macrophages and dendritic cells, which indicated that these myeloid lineages arise in early lymphoid lineage specification. Thus, as in the mouse, human hematopoiesis does not follow a rigid model of myeloid-lymphoid segregation.

A renewed model of pancreatic cancer evolution based on genomic rearrangement patterns.

Pancreatic cancer, a highly aggressive tumour type with uniformly poor prognosis, exemplifies the classically held view of stepwise cancer development. The current model of tumorigenesis, based on analyses of precursor lesions, termed pancreatic intraepithelial neoplasm (PanINs) lesions, makes two predictions: first, that pancreatic cancer develops through a particular sequence of genetic alterations (KRAS, followed by CDKN2A, then TP53 and SMAD4); and second, that the evolutionary trajectory of pancreatic cancer progression is gradual because each alteration is acquired independently. A shortcoming of this model is that clonally expanded precursor lesions do not always belong to the tumour lineage, indicating that the evolutionary trajectory of the tumour lineage and precursor lesions can be divergent. This prevailing model of tumorigenesis has contributed to the clinical notion that pancreatic cancer evolves slowly and presents at a late stage. However, the propensity for this disease to rapidly metastasize and the inability to improve patient outcomes, despite efforts aimed at early detection, suggest that pancreatic cancer progression is not gradual. Here, using newly developed informatics tools, we tracked changes in DNA copy number and their associated rearrangements in tumour-enriched genomes and found that pancreatic cancer tumorigenesis is neither gradual nor follows the accepted mutation order. Two-thirds of tumours harbour complex rearrangement patterns associated with mitotic errors, consistent with punctuated equilibrium as the principal evolutionary trajectory. In a subset of cases, the consequence of such errors is the simultaneous, rather than sequential, knockout of canonical preneoplastic genetic drivers that are likely to set-off invasive cancer growth. These findings challenge the current progression model of pancreatic cancer and provide insights into the mutational processes that give rise to these aggressive tumours.

Clinical and genomic characterisation of mismatch repair deficient pancreatic adenocarcinoma.

OBJECTIVE: To describe the clinical, pathological and genomic characteristics of pancreatic cancer with DNA mismatch repair deficiency (MMRD) and proficiency (MMRP). DESIGN: We identified patients with MMRD and MMRP pancreatic cancer in a clinical cohort (N=1213, 519 with genetic testing, 53 with immunohistochemistry (IHC)) and a genomic cohort (N=288 with whole-genome sequencing (WGS)). RESULTS: 12 out of 1213 (1.0%) in the clinical cohort were MMRD by IHC or WGS. Of the 14 patients with Lynch syndrome, 3 (21.4%) had an MMRP pancreatic cancer by IHC, and 4 (28.6%) were excluded because tissue was unavailable for testing. MMRD cancers had longer overall survival after surgery (weighted HR after coarsened exact matching 0.11, 95% CI 0.02 to 0.78, p=0.001). One patient with an unresectable MMRD cancer has an ongoing partial response 3 years after starting treatment with PD-L1/CTLA-4 inhibition. This tumour showed none of the classical histopathological features of MMRD. 9 out of 288 (3.1%) tumours with WGS were MMRD. Despite markedly higher tumour mutational burden and neoantigen loads, MMRD cancers were significantly less likely to have mutations in usual pancreatic cancer driver genes like KRAS and SMAD4, but more likely to have mutations in genes that drive cancers with microsatellite instability like ACV2RA and JAK1. MMRD tumours were significantly more likely to have a basal-like transcriptional programme and elevated transcriptional markers of immunogenicity. CONCLUSIONS: MMRD pancreatic cancers have distinct clinical, pathological and genomic profiles. Patients with MMRD pancreatic cancer should be considered for basket trials targeting enhanced immunogenicity or the unique genomic drivers in these malignancies.

Morphological classification of pancreatic ductal adenocarcinoma that predicts molecular subtypes and correlates with clinical outcome.

INTRODUCTION: Transcriptional analyses have identified several distinct molecular subtypes in pancreatic ductal adenocarcinoma (PDAC) that have prognostic and potential therapeutic significance. However, to date, an indepth, clinicomorphological correlation of these molecular subtypes has not been performed. We sought to identify specific morphological patterns to compare with known molecular subtypes, interrogate their biological significance, and furthermore reappraise the current grading system in PDAC. DESIGN: We first assessed 86 primary, chemotherapy-naive PDAC resection specimens with matched RNA-Seq data for specific, reproducible morphological patterns. Differential expression was applied to the gene expression data using the morphological features. We next compared the differentially expressed gene signatures with previously published molecular subtypes. Overall survival (OS) was correlated with the morphological and molecular subtypes. RESULTS: We identified four morphological patterns that segregated into two components ('gland forming' and 'non-gland forming') based on the presence/absence of well-formed glands. A morphological cut-off (≥40% 'non-gland forming') was established using RNA-Seq data, which identified two groups (A and B) with gene signatures that correlated with known molecular subtypes. There was a significant difference in OS between the groups. The morphological groups remained significantly prognostic within cancers that were moderately differentiated and classified as 'classical' using RNA-Seq. CONCLUSION: Our study has demonstrated that PDACs can be morphologically classified into distinct and biologically relevant categories which predict known molecular subtypes. These results provide the basis for an improved taxonomy of PDAC, which may lend itself to future treatment strategies and the development of deep learning models.

Mutations in Mitochondrial DNA From Pancreatic Ductal Adenocarcinomas Associate With Survival Times of Patients and Accumulate as Tumors Progress.

Somatic mutations have been found in the mitochondria in different types of cancer cells, but it is not clear whether these affect tumorigenesis or tumor progression. We analyzed mitochondrial genomes of 268 early-stage, resected pancreatic ductal adenocarcinoma tissues and paired non-tumor tissues. We defined a mitochondrial somatic mutation (mtSNV) as a position where the difference in heteroplasmy fraction between tumor and normal sample was ≥0.2. Our analysis identified 304 mtSNVs, with at least 1 mtSNV in 61% (164 of 268) of tumor samples. The noncoding control region had the greatest proportion of mtSNVs (60 of 304 mutations); this region contains sites that regulate mitochondrial DNA transcription and replication. Frequently mutated genes included ND5, RNR2, and CO1, plus 29 mutations in transfer RNA genes. mtSNVs in 2 separate mitochondrial genes (ND4 and ND6) were associated with shorter overall survival time. This association appeared to depend on the level of mtSNV heteroplasmy. Non-random co-occurrence between mtSNVs and mutations in nuclear genes indicates interactions between nuclear and mitochondrial DNA. In an analysis of primary tumors and metastases from 6 patients, we found tumors to accumulate mitochondrial mutational mutations as they progress.

The genetic basis of early T-cell precursor acute lymphoblastic leukaemia.

Early T-cell precursor acute lymphoblastic leukaemia (ETP ALL) is an aggressive malignancy of unknown genetic basis. We performed whole-genome sequencing of 12 ETP ALL cases and assessed the frequency of the identified somatic mutations in 94 T-cell acute lymphoblastic leukaemia cases. ETP ALL was characterized by activating mutations in genes regulating cytokine receptor and RAS signalling (67% of cases; NRAS, KRAS, FLT3, IL7R, JAK3, JAK1, SH2B3 and BRAF), inactivating lesions disrupting haematopoietic development (58%; GATA3, ETV6, RUNX1, IKZF1 and EP300) and histone-modifying genes (48%; EZH2, EED, SUZ12, SETD2 and EP300). We also identified new targets of recurrent mutation including DNM2, ECT2L and RELN. The mutational spectrum is similar to myeloid tumours, and moreover, the global transcriptional profile of ETP ALL was similar to that of normal and myeloid leukaemia haematopoietic stem cells. These findings suggest that addition of myeloid-directed therapies might improve the poor outcome of ETP ALL.

A genetic roadmap of pancreatic cancer: still evolving.

A diagnosis of pancreatic ductal adenocarcinoma (PDA) is often fatal. PDA is widely recognised as one of the 'incurable cancers' because therapies against this tumour type are generally ineffective. The fatal nature of this tumour is due to its aggressive clinical course. Pancreatic cancer commonly presents at the metastatic stage; even in cases where tumours are localised to the pancreas at diagnosis, metastatic seeds have often been invariably been spawned off, frustrating surgical attempts to cure the cancer. The key principles of pancreatic cancer mutational development were outlined nearly two decades ago using the genetics of precursor lesions to position the various stages of tumour progression. Since then, there has been a cavalcade of new data. How these recent studies impact the classical perceptions of pancreatic cancer development is a work in progress. Given that significant improvements in patient outcomes are not in sight for this disease, it is likely that broadening the current perspectives and acquiring deeper biological insights into the morphogenetic route of tumour development will be needed to foster new strategies for more effective cancer control.

Evolution of human BCR-ABL1 lymphoblastic leukaemia-initiating cells.

Many tumours are composed of genetically diverse cells; however, little is known about how diversity evolves or the impact that diversity has on functional properties. Here, using xenografting and DNA copy number alteration (CNA) profiling of human BCR-ABL1 lymphoblastic leukaemia, we demonstrate that genetic diversity occurs in functionally defined leukaemia-initiating cells and that many diagnostic patient samples contain multiple genetically distinct leukaemia-initiating cell subclones. Reconstructing the subclonal genetic ancestry of several samples by CNA profiling demonstrated a branching multi-clonal evolution model of leukaemogenesis, rather than linear succession. For some patient samples, the predominant diagnostic clone repopulated xenografts, whereas in others it was outcompeted by minor subclones. Reconstitution with the predominant diagnosis clone was associated with more aggressive growth properties in xenografts, deletion of CDKN2A and CDKN2B, and a trend towards poorer patient outcome. Our findings link clonal diversity with leukaemia-initiating-cell function and underscore the importance of developing therapies that eradicate all intratumoral subclones.

PLZF is a regulator of homeostatic and cytokine-induced myeloid development.

A major question in hematopoiesis is how the system maintains long-term homeostasis whereby the generation of large numbers of differentiated cells is balanced with the requirement for maintenance of progenitor pools, while remaining sufficiently flexible to respond to periods of perturbed cellular output during infection or stress. We focused on the development of the myeloid lineage and present evidence that promyelocytic leukemia zinc finger (PLZF) provides a novel function that is critical for both normal and stress-induced myelopoiesis. During homeostasis, PLZF restricts proliferation and differentiation of human cord blood-derived myeloid progenitors to maintain a balance between the progenitor and mature cell compartments. Analysis of PLZF promoter-binding sites revealed that it represses transcription factors involved in normal myeloid differentiation, including GFI-1, C/EBPalpha, and LEF-1, and induces negative regulators DUSP6 and ID2. Loss of ID2 relieves PLZF-mediated repression of differentiation identifying it as a functional target of PLZF in myelopoiesis. Furthermore, induction of ERK1/2 by myeloid cytokines, reflective of a stress response, leads to nuclear export and inactivation of PLZF, which augments mature cell production. Thus, negative regulators of differentiation can serve to maintain developmental systems in a primed state, so that their inactivation by extrinsic signals can induce proliferation and differentiation to rapidly satisfy increased demand for mature cells.

A small molecule screening strategy with validation on human leukemia stem cells uncovers the therapeutic efficacy of kinetin riboside.

Gene regulatory networks that govern hematopoietic stem cells (HSCs) and leukemia-initiating cells (L-ICs) are deeply entangled. Thus, the discovery of compounds that target L-ICs while sparing HSC is an attractive but difficult endeavor. Presently, most screening approaches fail to counter-screen compounds against normal hematopoietic stem/progenitor cells (HSPCs). Here, we present a multistep in vitro and in vivo approach to identify compounds that can target L-ICs in acute myeloid leukemia (AML). A high-throughput screen of 4000 compounds on novel leukemia cell lines derived from human experimental leukemogenesis models yielded 80 hits, of which 10 were less toxic to HSPC. We characterized a single compound, kinetin riboside (KR), on AML L-ICs and HSPCs. KR demonstrated comparable efficacy to standard therapies against blast cells in 63 primary leukemias. In vitro, KR targeted the L-IC-enriched CD34(+)CD38(-) AML fraction, while sparing HSPC-enriched fractions, although these effects were mitigated on HSC assayed in vivo. KR eliminated L-ICs in 2 of 4 primary AML samples when assayed in vivo and highlights the importance of in vivo L-IC and HSC assays to measure function. Overall, we provide a novel approach to screen large drug libraries for the discovery of anti-L-IC compounds for human leukemias.

Cytidine deaminases APOBEC3C and APOBEC3D promote DNA replication stress resistance in pancreatic cancer cells.

Gemcitabine is a potent inhibitor of DNA replication and is a mainstay therapeutic for diverse cancers, particularly pancreatic ductal adenocarcinoma (PDAC). However, most tumors remain refractory to gemcitabine therapies. Here, to define the cancer cell response to gemcitabine, we performed genome-scale CRISPR-Cas9 chemical-genetic screens in PDAC cells and found selective loss of cell fitness upon disruption of the cytidine deaminases APOBEC3C and APOBEC3D. Following gemcitabine treatment, APOBEC3C and APOBEC3D promote DNA replication stress resistance and cell survival by deaminating cytidines in the nuclear genome to ensure DNA replication fork restart and repair in PDAC cells. We provide evidence that the chemical-genetic interaction between APOBEC3C or APOBEC3D and gemcitabine is absent in nontransformed cells but is recapitulated across different PDAC cell lines, in PDAC organoids and in PDAC xenografts. Thus, we uncover roles for APOBEC3C and APOBEC3D in DNA replication stress resistance and offer plausible targets for improving gemcitabine-based therapies for PDAC.