Aberrant stem cell and developmental programs in pediatric leukemia.

Hematopoiesis is a finely orchestrated process, whereby hematopoietic stem cells give rise to all mature blood cells. Crucially, they maintain the ability to self-renew and/or differentiate to replenish downstream progeny. This process starts at an embryonic stage and continues throughout the human lifespan. Blood cancers such as leukemia occur when normal hematopoiesis is disrupted, leading to uncontrolled proliferation and a block in differentiation of progenitors of a particular lineage (myeloid or lymphoid). Although normal stem cell programs are crucial for tissue homeostasis, these can be co-opted in many cancers, including leukemia. Myeloid or lymphoid leukemias often display stem cell-like properties that not only allow proliferation and survival of leukemic blasts but also enable them to escape treatments currently employed to treat patients. In addition, some leukemias, especially in children, have a fetal stem cell profile, which may reflect the developmental origins of the disease. Aberrant fetal stem cell programs necessary for leukemia maintenance are particularly attractive therapeutic targets. Understanding how hijacked stem cell programs lead to aberrant gene expression in place and time, and drive the biology of leukemia, will help us develop the best treatment strategies for patients.

A human genome editing-based MLL::AF4 ALL model recapitulates key cellular and molecular leukemogenic features.

Cellular ontogeny and MLL breakpoint site influence the capacity of MLL-edited CD34+ hematopoietic cells to initiate and recapitulate infant patients' features in pro-B-cell acute lymphoblastic leukemia (B-ALL). We provide key insights into the leukemogenic determinants of MLL-AF4+ infant B-ALL.

MLL-AF4 cooperates with PAF1 and FACT to drive high-density enhancer interactions in leukemia.

Aberrant enhancer activation is a key mechanism driving oncogene expression in many cancers. While much is known about the regulation of larger chromosome domains in eukaryotes, the details of enhancer-promoter interactions remain poorly understood. Recent work suggests co-activators like BRD4 and Mediator have little impact on enhancer-promoter interactions. In leukemias controlled by the MLL-AF4 fusion protein, we use the ultra-high resolution technique Micro-Capture-C (MCC) to show that MLL-AF4 binding promotes broad, high-density regions of enhancer-promoter interactions at a subset of key targets. These enhancers are enriched for transcription elongation factors like PAF1C and FACT, and the loss of these factors abolishes enhancer-promoter contact. This work not only provides an additional model for how MLL-AF4 is able to drive high levels of transcription at key genes in leukemia but also suggests a more general model linking enhancer-promoter crosstalk and transcription elongation.

Human Bone Marrow Organoids for Disease Modeling, Discovery, and Validation of Therapeutic Targets in Hematologic Malignancies.

A lack of models that recapitulate the complexity of human bone marrow has hampered mechanistic studies of normal and malignant hematopoiesis and the validation of novel therapies. Here, we describe a step-wise, directed-differentiation protocol in which organoids are generated from induced pluripotent stem cells committed to mesenchymal, endothelial, and hematopoietic lineages. These 3D structures capture key features of human bone marrow-stroma, lumen-forming sinusoids, and myeloid cells including proplatelet-forming megakaryocytes. The organoids supported the engraftment and survival of cells from patients with blood malignancies, including cancer types notoriously difficult to maintain ex vivo. Fibrosis of the organoid occurred following TGFß stimulation and engraftment with myelofibrosis but not healthy donor-derived cells, validating this platform as a powerful tool for studies of malignant cells and their interactions within a human bone marrow-like milieu. This enabling technology is likely to accelerate the discovery and prioritization of novel targets for bone marrow disorders and blood cancers. SIGNIFICANCE: We present a human bone marrow organoid that supports the growth of primary cells from patients with myeloid and lymphoid blood cancers. This model allows for mechanistic studies of blood cancers in the context of their microenvironment and provides a much-needed ex vivo tool for the prioritization of new therapeutics. See related commentary by Derecka and Crispino, p. 263. This article is highlighted in the In This Issue feature, p. 247.

Scrutinizing Deleterious Nonsynonymous SNPs and Their Effect on Human POLD1 Gene.

POLD1 (DNA polymerase delta 1, catalytic subunit) is a protein-coding gene that encodes the large catalytic subunit (POLD1/p125) of the DNA polymerase delta (Polδ) complex. The consequence of missense or nonsynonymous SNPs (nsSNPs), which occur in the coding region of a specific gene, is the replacement of single amino acid. It may also change the structure, stability, and/or functions of the protein. Mutation in the POLD1 gene is associated with autosomal dominant predisposition to colonic adenomatous polyps, colon cancer, endometrial cancer (EDMC), breast cancer, and brain tumors. These de novo mutations in the POLD1 gene also result in autosomal dominant MDPL syndrome (mandibular hypoplasia, deafness, progeroid features, and lipodystrophy). In this study, genetic variations of POLD1 which may affect the structure and/or function were analyzed using different types of bioinformatics tools. A total of 17038 nsSNPs for POLD1 were collected from the NCBI database, among which 1317 were missense variants. Out of all missense nsSNPs, 28 were found to be deleterious functionally and structurally. Among these deleterious nsSNPs, 23 showed a conservation scale of >5, 2 were predicted to be associated with binding site formation, and one acted as a posttranslational modification site. All of them were involved in coil, extracellular structures, or helix formation, and some cause the change in size, charge, and hydrophobicity.

In utero origin of myelofibrosis presenting in adult monozygotic twins.

The latency between acquisition of an initiating somatic driver mutation by a single-cell and clinical presentation with cancer is largely unknown. We describe a remarkable case of monozygotic twins presenting with CALR mutation-positive myeloproliferative neoplasms (MPNs) (aged 37 and 38 years), with a clinical phenotype of primary myelofibrosis. The CALR mutation was absent in T cells and dermal fibroblasts, confirming somatic acquisition. Whole-genome sequencing lineage tracing revealed a common clonal origin of the CALR-mutant MPN clone, which occurred in utero followed by twin-to-twin transplacental transmission and subsequent similar disease latency. Index sorting and single-colony genotyping revealed phenotypic hematopoietic stem cells (HSCs) as the likely MPN-propagating cell. Furthermore, neonatal blood spot analysis confirmed in utero origin of the JAK2V617F mutation in a patient presenting with polycythemia vera (aged 34 years). These findings provide a unique window into the prolonged evolutionary dynamics of MPNs and fitness advantage exerted by MPN-associated driver mutations in HSCs.

Processing single-cell RNA-seq datasets using SingCellaR.

Single-cell RNA sequencing has led to unprecedented levels of data complexity. Although several computational platforms are available, performing data analyses for multiple datasets remains a significant challenge. Here, we provide a comprehensive analytical protocol to interrogate multiple datasets on SingCellaR, an analysis package in R. This tool can be applied to general single-cell transcriptome analyses. We demonstrate steps for data analyses and visualization using bespoke pipelines, in conjunction with existing analysis tools to study human hematopoietic stem and progenitor cells. For complete details on the use and execution of this protocol, please refer to Roy et al. (2021).

CD34+CD19-CD22+ B-cell progenitors may underlie phenotypic escape in patients treated with CD19-directed therapies.

CD19-directed immunotherapies have revolutionized the treatment of advanced B-cell acute lymphoblastic leukemia (B-ALL). Despite initial impressive rates of complete remission (CR) many patients ultimately relapse. Patients with B-ALL successfully treated with CD19-directed T cells eventually relapse, which, coupled with the early onset of CD22 expression during B-cell development, suggests that preexisting CD34+CD22+CD19- (pre)-leukemic cells represent an "early progenitor origin-related" mechanism underlying phenotypic escape to CD19-directed immunotherapies. We demonstrate that CD22 expression precedes CD19 expression during B-cell development. CD34+CD19-CD22+ cells are found in diagnostic and relapsed bone marrow samples of ∼70% of patients with B-ALL, and their frequency increases twofold in patients with B-ALL in CR after CD19 CAR T-cell therapy. The median of CD34+CD19-CD22+ cells before treatment was threefold higher in patients in whom B-ALL relapsed after CD19-directed immunotherapy (median follow-up, 24 months). Fluorescence in situ hybridization analysis in flow-sorted cell populations and xenograft modeling revealed that CD34+CD19-CD22+ cells harbor the genetic abnormalities present at diagnosis and initiate leukemogenesis in vivo. Our data suggest that preleukemic CD34+CD19-CD22+ progenitors underlie phenotypic escape after CD19-directed immunotherapies and reinforce ongoing clinical studies aimed at CD19/CD22 dual targeting as a strategy for reducing CD19- relapses. The implementation of CD34/CD19/CD22 immunophenotyping in clinical laboratories for initial diagnosis and subsequent monitoring of patients with B-ALL during CD19-targeted therapy is encouraged.

Efficient micropropagation of Thunbergia coccinea Wall. and genetic homogeneity assessment through RAPD and ISSR markers.

Thunbergia coccinea Wall. ex D. Don being a rare, ornamental and medicinal plant of India, is needed to propagate for conserving the germplasm and analyzing its phytochemical compounds in the future. A reliable protocol for direct in vitro propagation using nodal shoot meristem of T. coccinea as explant was standardized. The highest number of shoots per explant (22.17 ± 0.54) with maximum shoot length (2.36 ± 0.28) in cm was obtained in Murashige and Skoog (MS) medium supplemented with 9.70 µM of 6-furfurylaminopurine (Kinetin) and 0.053 µM of α-naphthaleneacetic acid (NAA) combination, among all the different plant growth regulators (PGR's) and concentrations tested. The aforesaid PGR's combination was optimum for axillary shoot bud induction and multiplication in T. coccinea. The best rooting was observed on the half-strength MS medium fortified with 2.68 µM NAA with the highest number of roots per shoot (3.75 ± 0.12) and maximum length (5.22 ± 0.32) in cm. All the in vitro raised plantlets were acclimatized in sterile sand and soil mixture (1:1) with a survival rate of 70% on earthen pots under greenhouse conditions. PCR-based RAPD (Random Amplified Polymorphic DNA) and ISSR (Inter-Simple Sequence Repeat) molecular markers were employed to determine the genetic homogeneity amongst the plantlets. Twelve (12) RAPD and nine (9) ISSR primers developed a total of 104 and 91 scorable bands, respectively. The band profiles of micropropagated plantlets were monomorphic to the mother, donor in vivo plant, and similarity values varied from 0.9542-1.000. The dendrogram generated through UPGMA (unweighted pair group method with arithmetic mean) showed 99% similarities amongst all tested plants confirming the genetic uniformity of in vitro raised plants.

Preparation of Nanosensors for Detecting the Activity of Glycosaminoglycan Cleaving Enzymes.

Glycosaminoglycans (GAGs) play crucial roles in several biological processes including cell division, angiogenesis, anticoagulation, neurogenesis, axon guidance and growth, and viral and bacterial infections among others. The GAG cleaving hydrolases/lyases play a major role in the control of GAG structures, functions, and turn over. Dysregulation of GAG cleaving enzymes in vivo are linked to a number of human diseases including cancer, diabetes, atherosclerosis, arthritis, inflammation, and cardiovascular diseases. Several GAG cleaving enzymes are widely used for studying GAG glycobiology: heparitinases, chondroitinases, heparanases, hyaluronidases, and keratanases. Herein, we describe a method to synthesize four distinct nanometal surface energy transfer (NSET)-based gold-GAG-dye conjugates (nanosensors). Heparin, chondroitin sulfate, heparan sulfate, and hyaluronic acid are covalently linked with distinct fluorescent dyes and then immobilized on gold nanoparticles (AuNPs) to build nanosensors that serve as excellent substrates for GAG cleaving enzymes. Upon treatment of nanosensors with their respective GAG cleaving enzymes, dye-labeled oligosaccharides/disaccharides are released from AuNPs resulting in enhanced fluorescence recovery. These nanosensors have a great promise as diagnostic tools in various human pathophysiological conditions for detecting dysregulated expression of GAG cleaving enzymes and also as a sensitive analytical tool for assessing the quality control of pharmaceutical grade heparin polysaccharides that are produced in millions of small- and medium-sized animal slaughter houses worldwide.

A human fetal liver-derived infant MLL-AF4 acute lymphoblastic leukemia model reveals a distinct fetal gene expression program.

Although 90% of children with acute lymphoblastic leukemia (ALL) are now cured, the prognosis for infant-ALL remains dismal. Infant-ALL is usually caused by a single genetic hit that arises in utero: an MLL/KMT2A gene rearrangement (MLL-r). This is sufficient to induce a uniquely aggressive and treatment-refractory leukemia compared to older children. The reasons for disparate outcomes in patients of different ages with identical driver mutations are unknown. Using the most common MLL-r in infant-ALL, MLL-AF4, as a disease model, we show that fetal-specific gene expression programs are maintained in MLL-AF4 infant-ALL but not in MLL-AF4 childhood-ALL. We use CRISPR-Cas9 gene editing of primary human fetal liver hematopoietic cells to produce a t(4;11)/MLL-AF4 translocation, which replicates the clinical features of infant-ALL and drives infant-ALL-specific and fetal-specific gene expression programs. These data support the hypothesis that fetal-specific gene expression programs cooperate with MLL-AF4 to initiate and maintain the distinct biology of infant-ALL.

Single-cell profiling of human bone marrow progenitors reveals mechanisms of failing erythropoiesis in Diamond-Blackfan anemia.

Ribosome dysfunction underlies the pathogenesis of many cancers and heritable ribosomopathies. Here, we investigate how mutations in either ribosomal protein large (RPL) or ribosomal protein small (RPS) subunit genes selectively affect erythroid progenitor development and clinical phenotypes in Diamond-Blackfan anemia (DBA), a rare ribosomopathy with limited therapeutic options. Using single-cell assays of patient-derived bone marrow, we delineated two distinct cellular trajectories segregating with ribosomal protein genotypes. Almost complete loss of erythroid specification was observed in RPS-DBA. In contrast, we observed relative preservation of qualitatively abnormal erythroid progenitors and precursors in RPL-DBA. Although both DBA genotypes exhibited a proinflammatory bone marrow milieu, RPS-DBA was characterized by erythroid differentiation arrest, whereas RPL-DBA was characterized by preserved GATA1 expression and activity. Compensatory stress erythropoiesis in RPL-DBA exhibited disordered differentiation underpinned by an altered glucocorticoid molecular signature, including reduced ZFP36L2 expression, leading to milder anemia and improved corticosteroid response. This integrative analysis approach identified distinct pathways of erythroid failure and defined genotype-phenotype correlations in DBA. These findings may help facilitate therapeutic target discovery.

Transitions in lineage specification and gene regulatory networks in hematopoietic stem/progenitor cells over human development.

Human hematopoiesis is a dynamic process that starts in utero 18-21 days post-conception. Understanding the site- and stage-specific variation in hematopoiesis is important if we are to understand the origin of hematological disorders, many of which occur at specific points in the human lifespan. To unravel how the hematopoietic stem/progenitor cell (HSPC) compartment changes during human ontogeny and the underlying gene regulatory mechanisms, we compare 57,489 HSPCs from 5 different tissues spanning 4 developmental stages through the human lifetime. Single-cell transcriptomic analysis identifies significant site- and developmental stage-specific transitions in cellular architecture and gene regulatory networks. Hematopoietic stem cells show progression from cycling to quiescence and increased inflammatory signaling during ontogeny. We demonstrate the utility of this dataset for understanding aberrant hematopoiesis through comparison to two cancers that present at distinct time points in postnatal life-juvenile myelomonocytic leukemia, a childhood cancer, and myelofibrosis, which classically presents in older adults.

A KMT2A-AFF1 gene regulatory network highlights the role of core transcription factors and reveals the regulatory logic of key downstream target genes.

Regulatory interactions mediated by transcription factors (TFs) make up complex networks that control cellular behavior. Fully understanding these gene regulatory networks (GRNs) offers greater insight into the consequences of disease-causing perturbations than can be achieved by studying single TF binding events in isolation. Chromosomal translocations of the lysine methyltransferase 2A (KMT2A) gene produce KMT2A fusion proteins such as KMT2A-AFF1 (previously MLL-AF4), causing poor prognosis acute lymphoblastic leukemias (ALLs) that sometimes relapse as acute myeloid leukemias (AMLs). KMT2A-AFF1 drives leukemogenesis through direct binding and inducing the aberrant overexpression of key genes, such as the anti-apoptotic factor BCL2 and the proto-oncogene MYC However, studying direct binding alone does not incorporate possible network-generated regulatory outputs, including the indirect induction of gene repression. To better understand the KMT2A-AFF1-driven regulatory landscape, we integrated ChIP-seq, patient RNA-seq, and CRISPR essentiality screens to generate a model GRN. This GRN identified several key transcription factors such as RUNX1 that regulate target genes downstream of KMT2A-AFF1 using feed-forward loop (FFL) and cascade motifs. A core set of nodes are present in both ALL and AML, and CRISPR screening revealed several factors that help mediate response to the drug venetoclax. Using our GRN, we then identified a KMT2A-AFF1:RUNX1 cascade that represses CASP9, as well as KMT2A-AFF1-driven FFLs that regulate BCL2 and MYC through combinatorial TF activity. This illustrates how our GRN can be used to better connect KMT2A-AFF1 behavior to downstream pathways that contribute to leukemogenesis, and potentially predict shifts in gene expression that mediate drug response.

The Origin of B-cells: Human Fetal B Cell Development and Implications for the Pathogenesis of Childhood Acute Lymphoblastic Leukemia.

Human B-lymphopoiesis is a dynamic life-long process that starts in utero by around six post-conception weeks. A detailed understanding of human fetal B-lymphopoiesis and how it changes in postnatal life is vital for building a complete picture of normal B-lymphoid development through ontogeny, and its relevance in disease. B-cell acute lymphoblastic leukemia (B-ALL) is one of the most common cancers in children, with many of the leukemia-initiating events originating in utero. It is likely that the biology of B-ALL, including leukemia initiation, maintenance and progression depends on the developmental stage and type of B-lymphoid cell in which it originates. This is particularly important for early life leukemias, where specific characteristics of fetal B-cells might be key to determining how the disease behaves, including response to treatment. These cellular, molecular and/or epigenetic features are likely to change with age in a cell intrinsic and/or microenvironment directed manner. Most of our understanding of fetal B-lymphopoiesis has been based on murine data, but many recent studies have focussed on characterizing human fetal B-cell development, including functional and molecular assays at a single cell level. In this mini-review we will give a short overview of the recent advances in the understanding of human fetal B-lymphopoiesis, including its relevance to infant/childhood leukemia, and highlight future questions in the field.

The genome-wide impact of trisomy 21 on DNA methylation and its implications for hematopoiesis.

Down syndrome is associated with genome-wide perturbation of gene expression, which may be mediated by epigenetic changes. We perform an epigenome-wide association study on neonatal bloodspots comparing 196 newborns with Down syndrome and 439 newborns without Down syndrome, adjusting for cell-type heterogeneity, which identifies 652 epigenome-wide significant CpGs (P < 7.67 × 10-8) and 1,052 differentially methylated regions. Differential methylation at promoter/enhancer regions correlates with gene expression changes in Down syndrome versus non-Down syndrome fetal liver hematopoietic stem/progenitor cells (P < 0.0001). The top two differentially methylated regions overlap RUNX1 and FLI1, both important regulators of megakaryopoiesis and hematopoietic development, with significant hypermethylation at promoter regions of these two genes. Excluding Down syndrome newborns harboring preleukemic GATA1 mutations (N = 30), identified by targeted sequencing, has minimal impact on the epigenome-wide association study results. Down syndrome has profound, genome-wide effects on DNA methylation in hematopoietic cells in early life, which may contribute to the high frequency of hematological problems, including leukemia, in children with Down syndrome.

Heterogeneous disease-propagating stem cells in juvenile myelomonocytic leukemia.

Juvenile myelomonocytic leukemia (JMML) is a poor-prognosis childhood leukemia usually caused by RAS-pathway mutations. The cellular hierarchy in JMML is poorly characterized, including the identity of leukemia stem cells (LSCs). FACS and single-cell RNA sequencing reveal marked heterogeneity of JMML hematopoietic stem/progenitor cells (HSPCs), including an aberrant Lin-CD34+CD38-CD90+CD45RA+ population. Single-cell HSPC index-sorting and clonogenic assays show that (1) all somatic mutations can be backtracked to the phenotypic HSC compartment, with RAS-pathway mutations as a "first hit," (2) mutations are acquired with both linear and branching patterns of clonal evolution, and (3) mutant HSPCs are present after allogeneic HSC transplant before molecular/clinical evidence of relapse. Stem cell assays reveal interpatient heterogeneity of JMML LSCs, which are present in, but not confined to, the phenotypic HSC compartment. RNA sequencing of JMML LSC reveals up-regulation of stem cell and fetal genes (HLF, MEIS1, CNN3, VNN2, and HMGA2) and candidate therapeutic targets/biomarkers (MTOR, SLC2A1, and CD96), paving the way for LSC-directed disease monitoring and therapy in this disease.

H3K79me2/3 controls enhancer-promoter interactions and activation of the pan-cancer stem cell marker PROM1/CD133 in MLL-AF4 leukemia cells.

MLL gene rearrangements (MLLr) are a common cause of aggressive, incurable acute lymphoblastic leukemias (ALL) in infants and children, most of which originate in utero. The most common MLLr produces an MLL-AF4 fusion protein. MLL-AF4 promotes leukemogenesis by activating key target genes, mainly through recruitment of DOT1L and increased histone H3 lysine-79 methylation (H3K79me2/3). One key MLL-AF4 target gene is PROM1, which encodes CD133 (Prominin-1). CD133 is a pentaspan transmembrane glycoprotein that represents a potential pan-cancer target as it is found on multiple cancer stem cells. Here we demonstrate that aberrant PROM1/CD133 expression is essential for leukemic cell growth, mediated by direct binding of MLL-AF4. Activation is controlled by an intragenic H3K79me2/3 enhancer element (KEE) leading to increased enhancer-promoter interactions between PROM1 and the nearby gene TAPT1. This dual locus regulation is reflected in a strong correlation of expression in leukemia. We find that in PROM1/CD133 non-expressing cells, the PROM1 locus is repressed by polycomb repressive complex 2 (PRC2) binding, associated with reduced expression of TAPT1, partially due to loss of interactions with the PROM1 locus. Together, these results provide the first detailed analysis of PROM1/CD133 regulation that explains CD133 expression in MLLr ALL.