A latent subset of human hematopoietic stem cells resists regenerative stress to preserve stemness.

Continuous supply of immune cells throughout life relies on the delicate balance in the hematopoietic stem cell (HSC) pool between long-term maintenance and meeting the demands of both normal blood production and unexpected stress conditions. Here we identified distinct subsets of human long-term (LT)-HSCs that responded differently to regeneration-mediated stress: an immune checkpoint ligand CD112lo subset that exhibited a transient engraftment restraint (termed latency) before contributing to hematopoietic reconstitution and a primed CD112hi subset that responded rapidly. This functional heterogeneity and CD112 expression are regulated by INKA1 through direct interaction with PAK4 and SIRT1, inducing epigenetic changes and defining an alternative state of LT-HSC quiescence that serves to preserve self-renewal and regenerative capacity upon regeneration-mediated stress. Collectively, our data uncovered the molecular intricacies underlying HSC heterogeneity and self-renewal regulation and point to latency as an orchestrated physiological response that balances blood cell demands with preserving a stem cell reservoir.

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.

BPG: Seamless, automated and interactive visualization of scientific data.

BACKGROUND: We introduce BPG, a framework for generating publication-quality, highly-customizable plots in the R statistical environment. RESULTS: This open-source package includes multiple methods of displaying high-dimensional datasets and facilitates generation of complex multi-panel figures, making it suitable for complex datasets. A web-based interactive tool allows online figure customization, from which R code can be downloaded for integration with computational pipelines. CONCLUSION: BPG provides a new approach for linking interactive and scripted data visualization and is available at http://labs.oicr.on.ca/boutros-lab/software/bpg or via CRAN at https://cran.r-project.org/web/packages/BoutrosLab.plotting.general.

Molecular heterogeneity of non-small cell lung carcinoma patient-derived xenografts closely reflect their primary tumors.

Availability of lung cancer models that closely mimic human tumors remains a significant gap in cancer research, as tumor cell lines and mouse models may not recapitulate the spectrum of lung cancer heterogeneity seen in patients. We aimed to establish a patient-derived tumor xenograft (PDX) resource from surgically resected non-small cell lung cancer (NSCLC). Fresh tumor tissue from surgical resection was implanted and grown in the subcutaneous pocket of non-obese severe combined immune deficient (NOD SCID) gamma mice. Subsequent passages were in NOD SCID mice. A subset of matched patient and PDX tumors and non-neoplastic lung tissues were profiled by whole exome sequencing, single nucleotide polymorphism (SNP) and methylation arrays, and phosphotyrosine (pY)-proteome by mass spectrometry. The data were compared to published NSCLC datasets of NSCLC primary and cell lines. 127 stable PDXs were established from 441 lung carcinomas representing all major histological subtypes: 52 adenocarcinomas, 62 squamous cell carcinomas, one adeno-squamous carcinoma, five sarcomatoid carcinomas, five large cell neuroendocrine carcinomas, and two small cell lung cancers. Somatic mutations, gene copy number and expression profiles, and pY-proteome landscape of 36 PDXs showed greater similarity with patient tumors than with established cell lines. Novel somatic mutations on cancer associated genes were identified but only in PDXs, likely due to selective clonal growth in the PDXs that allows detection of these low allelic frequency mutations. The results provide the strongest evidence yet that PDXs established from lung cancers closely mimic the characteristics of patient primary tumors.

SeqControl: process control for DNA sequencing.

As high-throughput sequencing continues to increase in speed and throughput, routine clinical and industrial application draws closer. These 'production' settings will require enhanced quality monitoring and quality control to optimize output and reduce costs. We developed SeqControl, a framework for predicting sequencing quality and coverage using a set of 15 metrics describing overall coverage, coverage distribution, basewise coverage and basewise quality. Using whole-genome sequences of 27 prostate cancers and 26 normal references, we derived multivariate models that predict sequencing quality and depth. SeqControl robustly predicted how much sequencing was required to reach a given coverage depth (area under the curve (AUC) = 0.993), accurately classified clinically relevant formalin-fixed, paraffin-embedded samples, and made predictions from as little as one-eighth of a sequencing lane (AUC = 0.967). These techniques can be immediately incorporated into existing sequencing pipelines to monitor data quality in real time. SeqControl is available at http://labs.oicr.on.ca/Boutros-lab/software/SeqControl/.

Identification of genes expressed by immune cells of the colon that are regulated by colorectal cancer-associated variants.

A locus on human chromosome 11q23 tagged by marker rs3802842 was associated with colorectal cancer (CRC) in a genome-wide association study; this finding has been replicated in case-control studies worldwide. In order to identify biologic factors at this locus that are related to the etiopathology of CRC, we used microarray-based target selection methods, coupled to next-generation sequencing, to study 103 kb at the 11q23 locus. We genotyped 369 putative variants from 1,030 patients with CRC (cases) and 1,061 individuals without CRC (controls) from the Ontario Familial Colorectal Cancer Registry. Two previously uncharacterized genes, COLCA1 and COLCA2, were found to be co-regulated genes that are transcribed from opposite strands. Expression levels of COLCA1 and COLCA2 transcripts correlate with rs3802842 genotypes. In colon tissues, COLCA1 co-localizes with crystalloid granules of eosinophils and granular organelles of mast cells, neutrophils, macrophages, dendritic cells and differentiated myeloid-derived cell lines. COLCA2 is present in the cytoplasm of normal epithelial, immune and other cell lineages, as well as tumor cells. Tissue microarray analysis demonstrates the association of rs3802842 with lymphocyte density in the lamina propria (p = 0.014) and levels of COLCA1 in the lamina propria (p = 0.00016) and COLCA2 (tumor cells, p = 0.0041 and lamina propria, p = 6 × 10(-5)). In conclusion, genetic, expression and immunohistochemical data implicate COLCA1 and COLCA2 in the pathogenesis of colon cancer. Histologic analyses indicate the involvement of immune pathways.

Identifying molecular features that distinguish fluvastatin-sensitive breast tumor cells.

Statins, routinely used to treat hypercholesterolemia, selectively induce apoptosis in some tumor cells by inhibiting the mevalonate pathway. Recent clinical studies suggest that a subset of breast tumors is particularly susceptible to lipophilic statins, such as fluvastatin. To quickly advance statins as effective anticancer agents for breast cancer treatment, it is critical to identify the molecular features defining this sensitive subset. We have therefore characterized fluvastatin sensitivity by MTT assay in a panel of 19 breast cell lines that reflect the molecular diversity of breast cancer, and have evaluated the association of sensitivity with several clinicopathological and molecular features. A wide range of fluvastatin sensitivity was observed across breast tumor cell lines, with fluvastatin triggering cell death in a subset of sensitive cell lines. Fluvastatin sensitivity was associated with an estrogen receptor alpha (ERα)-negative, basal-like tumor subtype, features that can be scored with routine and/or strong preclinical diagnostics. To ascertain additional candidate sensitivity-associated molecular features, we mined publicly available gene expression datasets, identifying genes encoding regulators of mevalonate production, non-sterol lipid homeostasis, and global cellular metabolism, including the oncogene MYC. Further exploration of this data allowed us to generate a 10-gene mRNA abundance signature predictive of fluvastatin sensitivity, which showed preliminary validation in an independent set of breast tumor cell lines. Here, we have therefore identified several candidate predictors of sensitivity to fluvastatin treatment in breast cancer, which warrant further preclinical and clinical evaluation.

The genomic and evolutionary landscapes of anaplastic thyroid carcinoma.

Anaplastic thyroid carcinoma is arguably the most lethal human malignancy. It often co-occurs with differentiated thyroid cancers, yet the molecular origins of its aggressivity are unknown. We sequenced tumor DNA from 329 regions of thyroid cancer, including 213 from patients with primary anaplastic thyroid carcinomas. We also whole genome sequenced 9 patients using multi-region sequencing of both differentiated and anaplastic thyroid cancer components. Using these data, we demonstrate thatanaplastic thyroid carcinomas have a higher burden of mutations than other thyroid cancers, with distinct mutational signatures and molecular subtypes. Further, different cancer driver genes are mutated in anaplastic and differentiated thyroid carcinomas, even those arising in a single patient. Finally, we unambiguously demonstrate that anaplastic thyroid carcinomas share a genomic origin with co-occurring differentiated carcinomas and emerge from a common malignant field through acquisition of characteristic clonal driver mutations.

Transcriptomic classes of BCR-ABL1 lymphoblastic leukemia.

In BCR-ABL1 lymphoblastic leukemia, treatment heterogeneity to tyrosine kinase inhibitors (TKIs), especially in the absence of kinase domain mutations in BCR-ABL1, is poorly understood. Through deep molecular profiling, we uncovered three transcriptomic subtypes of BCR-ABL1 lymphoblastic leukemia, each representing a maturation arrest at a stage of B-cell progenitor differentiation. An earlier arrest was associated with lineage promiscuity, treatment refractoriness and poor patient outcomes. A later arrest was associated with lineage fidelity, durable leukemia remissions and improved patient outcomes. Each maturation arrest was marked by specific genomic events that control different transition points in B-cell development. Interestingly, these events were absent in BCR-ABL1+ preleukemic stem cells isolated from patients regardless of subtype, which supports that transcriptomic phenotypes are determined downstream of the leukemia-initialing event. Overall, our data indicate that treatment response and TKI efficacy are unexpected outcomes of the differentiation stage at which this leukemia transforms.

SINGLE CELL DISSECTION OF DEVELOPMENTAL ORIGINS AND TRANSCRIPTIONAL HETEROGENEITY IN B-CELL ACUTE LYMPHOBLASTIC LEUKEMIA.

Sequencing of bulk tumor populations has improved genetic classification and risk assessment of B-ALL, but does not directly examine intratumor heterogeneity or infer leukemia cellular origins. We profiled 89 B-ALL samples by single-cell RNA-seq (scRNA-seq) and compared them to a reference map of normal human B-cell development established using both functional and molecular assays. Intra-sample heterogeneity was driven by cell cycle, metabolism, differentiation, and inflammation transcriptional programs. By inference of B lineage developmental state composition, nearly all samples possessed a high abundance of pro-B cells, with variation between samples mainly driven by sub-populations. However, ZNF384- r and DUX4- r B-ALL showed composition enrichment of hematopoietic stem cells, BCR::ABL1 and KMT2A -r ALL of Early Lymphoid progenitors, MEF2D -r and TCF3::PBX1 of Pre-B cells. Enrichment of Early Lymphoid progenitors correlated with high-risk clinical features. Understanding variation in transcriptional programs and developmental states of B-ALL by scRNA-seq refines existing clinical and genomic classifications and improves prediction of treatment outcome.

A primary hierarchically organized patient-derived model enables in depth interrogation of stemness driven by the coding and non-coding genome.

Many cancers are organized as cellular hierarchies sustained by cancer stem cells (CSC), whose eradication is crucial for achieving long-term remission. Difficulties to isolate and undertake in vitro and in vivo experimental studies of rare CSC under conditions that preserve their original properties currently constitute a bottleneck for identifying molecular mechanisms involving coding and non-coding genomic regions that govern stemness. We focussed on acute myeloid leukemia (AML) as a paradigm of the CSC model and developed a patient-derived system termed OCI-AML22 that recapitulates the cellular hierarchy driven by leukemia stem cells (LSC). Through classical flow sorting and functional analyses, we established that a single phenotypic population is highly enriched for LSC. The LSC fraction can be easily isolated and serially expanded in culture or in xenografts while faithfully recapitulating functional, transcriptional and epigenetic features of primary LSCs. A novel non-coding regulatory element was identified with a new computational approach using functionally validated primary AML LSC fractions and its role in LSC stemness validated through efficient CRISPR editing using methods optimized for OCI-AML22 LSC. Collectively, OCI-AML22 constitutes a valuable resource to uncover mechanisms governing CSC driven malignancies.

A cellular hierarchy framework for understanding heterogeneity and predicting drug response in acute myeloid leukemia.

The treatment landscape of acute myeloid leukemia (AML) is evolving, with promising therapies entering clinical translation, yet patient responses remain heterogeneous, and biomarkers for tailoring treatment are lacking. To understand how disease heterogeneity links with therapy response, we determined the leukemia cell hierarchy makeup from bulk transcriptomes of more than 1,000 patients through deconvolution using single-cell reference profiles of leukemia stem, progenitor and mature cell types. Leukemia hierarchy composition was associated with functional, genomic and clinical properties and converged into four overall classes, spanning Primitive, Mature, GMP and Intermediate. Critically, variation in hierarchy composition along the Primitive versus GMP or Primitive versus Mature axes were associated with response to chemotherapy or drug sensitivity profiles of targeted therapies, respectively. A seven-gene biomarker derived from the Primitive versus Mature axis was associated with response to 105 investigational drugs. Cellular hierarchy composition constitutes a novel framework for understanding disease biology and advancing precision medicine in AML.

Multiomic Profiling of Central Nervous System Leukemia Identifies mRNA Translation as a Therapeutic Target.

Central nervous system (CNS) dissemination of B-precursor acute lymphoblastic leukemia (B-ALL) has poor prognosis and remains a therapeutic challenge. Here we performed targeted DNA sequencing as well as transcriptional and proteomic profiling of paired leukemia-infiltrating cells in the bone marrow (BM) and CNS of xenografts. Genes governing mRNA translation were upregulated in CNS leukemia, and subclonal genetic profiling confirmed this in both BM-concordant and BM-discordant CNS mutational populations. CNS leukemia cells were exquisitely sensitive to the translation inhibitor omacetaxine mepesuccinate, which reduced xenograft leptomeningeal disease burden. Proteomics demonstrated greater abundance of secreted proteins in CNS-infiltrating cells, including complement component 3 (C3), and drug targeting of C3 influenced CNS disease in xenografts. CNS-infiltrating cells also exhibited selection for stemness traits and metabolic reprogramming. Overall, our study identifies targeting of mRNA translation as a potential therapeutic approach for B-ALL leptomeningeal disease. SIGNIFICANCE: Cancer metastases are often driven by distinct subclones with unique biological properties. Here we show that in B-ALL CNS disease, the leptomeningeal environment selects for cells with unique functional dependencies. Pharmacologic inhibition of mRNA translation signaling treats CNS disease and offers a new therapeutic approach for this condition.This article is highlighted in the In This Issue feature, p. 1.

The Transition from Quiescent to Activated States in Human Hematopoietic Stem Cells Is Governed by Dynamic 3D Genome Reorganization.

Lifelong blood production requires long-term hematopoietic stem cells (LT-HSCs), marked by stemness states involving quiescence and self-renewal, to transition into activated short-term HSCs (ST-HSCs) with reduced stemness. As few transcriptional changes underlie this transition, we used single-cell and bulk assay for transposase-accessible chromatin sequencing (ATAC-seq) on human HSCs and hematopoietic stem and progenitor cell (HSPC) subsets to uncover chromatin accessibility signatures, one including LT-HSCs (LT/HSPC signature) and another excluding LT-HSCs (activated HSPC [Act/HSPC] signature). These signatures inversely correlated during early hematopoietic commitment and differentiation. The Act/HSPC signature contains CCCTC-binding factor (CTCF) binding sites mediating 351 chromatin interactions engaged in ST-HSCs, but not LT-HSCs, enclosing multiple stemness pathway genes active in LT-HSCs and repressed in ST-HSCs. CTCF silencing derepressed stemness genes, restraining quiescent LT-HSCs from transitioning to activated ST-HSCs. Hence, 3D chromatin interactions centrally mediated by CTCF endow a gatekeeper function that governs the earliest fate transitions HSCs make by coordinating disparate stemness pathways linked to quiescence and self-renewal.

TFEB-mediated endolysosomal activity controls human hematopoietic stem cell fate.

It is critical to understand how human quiescent long-term hematopoietic stem cells (LT-HSCs) sense demand from daily and stress-mediated cues and then transition into bioenergetically active progeny to differentiate and meet these cellular needs. However, the demand-adapted regulatory circuits of these early steps of hematopoiesis are largely unknown. Here we show that lysosomes, sophisticated nutrient-sensing and signaling centers, are regulated dichotomously by transcription factor EB (TFEB) and MYC to balance catabolic and anabolic processes required for activating LT-HSCs and guiding their lineage fate. TFEB-mediated induction of the endolysosomal pathway causes membrane receptor degradation, limiting LT-HSC metabolic and mitogenic activation, promoting quiescence and self-renewal, and governing erythroid-myeloid commitment. In contrast, MYC engages biosynthetic processes while repressing lysosomal catabolism, driving LT-HSC activation. Our study identifies TFEB-mediated control of lysosomal activity as a central regulatory hub for proper and coordinated stem cell fate determination.

Sphingosine-1-phosphate receptor 3 potentiates inflammatory programs in normal and leukemia stem cells to promote differentiation.

Acute myeloid leukemia (AML) is a caricature of normal hematopoiesis, driven from leukemia stem cells (LSC) that share some hematopoietic stem cell (HSC) programs including responsiveness to inflammatory signaling. Although inflammation dysregulates mature myeloid cells and influences stemness programs and lineage determination in HSC by activating stress myelopoiesis, such roles in LSC are poorly understood. Here, we show that S1PR3, a receptor for the bioactive lipid sphingosine-1-phosphate, is a central regulator which drives myeloid differentiation and activates inflammatory programs in both HSC and LSC. S1PR3-mediated inflammatory signatures varied in a continuum from primitive to mature myeloid states across AML patient cohorts, each with distinct phenotypic and clinical properties. S1PR3 was high in LSC and blasts of mature myeloid samples with linkages to chemosensitivity, while S1PR3 activation in primitive samples promoted LSC differentiation leading to eradication. Our studies open new avenues for therapeutic target identification specific for each AML subset.

Mapping the cellular origin and early evolution of leukemia in Down syndrome.

Children with Down syndrome have a 150-fold increased risk of developing myeloid leukemia, but the mechanism of predisposition is unclear. Because Down syndrome leukemogenesis initiates during fetal development, we characterized the cellular and developmental context of preleukemic initiation and leukemic progression using gene editing in human disomic and trisomic fetal hematopoietic cells and xenotransplantation. GATA binding protein 1 (GATA1) mutations caused transient preleukemia when introduced into trisomy 21 long-term hematopoietic stem cells, where a subset of chromosome 21 microRNAs affected predisposition to preleukemia. By contrast, progression to leukemia was independent of trisomy 21 and originated in various stem and progenitor cells through additional mutations in cohesin genes. CD117+/KIT proto-oncogene (KIT) cells mediated the propagation of preleukemia and leukemia, and KIT inhibition targeted preleukemic stem cells.

Relapse-Fated Latent Diagnosis Subclones in Acute B Lineage Leukemia Are Drug Tolerant and Possess Distinct Metabolic Programs.

Disease recurrence causes significant mortality in B-progenitor acute lymphoblastic leukemia (B-ALL). Genomic analysis of matched diagnosis and relapse samples shows relapse often arising from minor diagnosis subclones. However, why therapy eradicates some subclones while others survive and progress to relapse remains obscure. Elucidation of mechanisms underlying these differing fates requires functional analysis of isolated subclones. Here, large-scale limiting dilution xenografting of diagnosis and relapse samples, combined with targeted sequencing, identified and isolated minor diagnosis subclones that initiate an evolutionary trajectory toward relapse [termed diagnosis Relapse Initiating clones (dRI)]. Compared with other diagnosis subclones, dRIs were drug-tolerant with distinct engraftment and metabolic properties. Transcriptionally, dRIs displayed enrichment for chromatin remodeling, mitochondrial metabolism, proteostasis programs, and an increase in stemness pathways. The isolation and characterization of dRI subclones reveals new avenues for eradicating dRI cells by targeting their distinct metabolic and transcriptional pathways before further evolution renders them fully therapy-resistant. SIGNIFICANCE: Isolation and characterization of subclones from diagnosis samples of patients with B-ALL who relapsed showed that relapse-fated subclones had increased drug tolerance and distinct metabolic and survival transcriptional programs compared with other diagnosis subclones. This study provides strategies to identify and target clinically relevant subclones before further evolution toward relapse.

Author Correction: Transcription phenotypes of pancreatic cancer are driven by genomic events during tumor evolution.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.