Abstract
In leukemia, a distinct subpopulation of cancer-initiating cells called leukemia stem cells (LSCs) is believed to drive population expansion and tumor growth. Failing to eliminate LSCs may result in disease relapse regardless of the amount of non-LSCs destroyed. The first step in targeting and eliminating LSCs is to identify and characterize them. Acute precursor B lymphoblastic leukemia (B-ALL) cells derived from patients were incubated with fluorescent glucose analog 2-(N-(7-Nitrobenz-2-oxa-1, 3-diazol-4-yl) Amino)-2-Deoxyglucose (NBDG) and sorted based on NBDG uptake. Cell subpopulations defined by glucose uptake were then serially transplanted into mice and evaluated for leukemia initiating capacity. Gene expression profiles of these cells were characterized using RNA-Sequencing (RNA-Seq). A distinct population of NBDG-low cells was identified in patient B-ALL samples. These cells are a small population (1.92% of the entire leukemia population), have lower HLA expression, and are smaller in size (4.0 to 7.0 μm) than the rest of the leukemia population. All mice transplanted with NBDG-low cells developed leukemia between 5 and 14 weeks, while those transplanted with NBDG-high cells did not develop leukemia (p ≤ 0.0001-0.002). Serial transplantation of the NBDG-low mouse model resulted in successful leukemia development. NBDG-medium (NBDG-med) populations also developed leukemia. Interestingly, comprehensive molecular characterization of NBDG-low and NBDG-med cells from patient-derived xenograft (PDX) models using RNA-Seq revealed a distinct profile of 2,162 differentially-expressed transcripts (DETs) (p<0.05) with 70.6% down-regulated in NBDG-low cells. Hierarchical clustering of DETs showed distinct segregation of NBDG-low from NBDG-med and NBDG-high groups with marked transcription expression alterations in the NBDG-low group consistent with cancer survival. In conclusion, A unique subpopulation of cells with low glucose uptake (NBDG-low) in B-ALL was discovered. These cells, despite their quiescence characteristics, once transplanted in mice, showed potent leukemia initiating capacity. Although NBDG-med cells also initiated leukemia, gene expression profiling revealed a distinct signature that clearly distinguishes NBDG-low cells from NBDG-med and the rest of the leukemia populations. These results suggest that NBDG-low cells may represent quiescent LSCs. These cells can be activated in the appropriate environment in vivo, showing leukemia initiating capacity. Our study provides insight into the biologic mechanisms of B-ALL initiation and survival.
Introduction
Leukemia is a classic example of a disease of aberrant differentiation. Specifically, leukemia results from disrupted differentiation of pluripotent bone marrow progenitors into mature and functional mononuclear cells. Accordingly, leukemias are classified based on the intermediate stage at which differentiation is halted. In addition to lacking a terminally-differentiated phenotype and associated functional attributes, leukemia cells retain proliferative potential to continuously expand the population (1, 2). In some leukemias, as well as in most other cancers, it is now recognized that a rare but distinct subpopulation of cancer-initiating cells is responsible for population expansion and maintenance of tumor growth (3, 4). These leukemia stem cells (LSCs) are believed to make up an extremely small percentage of tumor cells, but they possess robust proliferative and self-renewal properties to facilitate sustained survival.
Despite the significance of LSCs in the development of curative leukemia treatments, definitive identification and characterization of LSCs in acute lymphoblastic leukemia (ALL) remains controversial (5, 6). While multiple leukemic subpopulations have been shown to be capable of propagating leukemia, their hierarchical structure is unclear or inconsistent. In addition, demonstration of leukemia initiating capacity alone is insufficient to validate LSCs, as it has been suggested that multiple ALL subpopulations with leukemia initiating capacity switch between dormancy and active proliferation, indicating a dynamic architecture of ALL development (7). As a result, it is very possible that ALL LSCs are by nature heterogeneous.
There are also few reliable distinguishing markers to characterize potential ALL LSCs. For example, hematopoietic stem cells (HSC) and LSCs in acute myeloid leukemia (AML) are traditionally defined based on phenotypic markers, including CD34+/CD38− (8–11). Characterization of these markers allows stem cells to be distinguished from their more differentiated progeny, which is a particularly crucial undertaking in the context of the heterogeneity of LSCs. In contrast, the phenotypic markers of stem cells in ALL, either B-cell or T-cell type, remain unclear. Previous studies demonstrated that subpopulations of ALL cells with CD34, CD38, CD20 or CD19 positive and/or negative markers could give rise to a leukemia phenotype in immunodeficient mice (12–14). These conflicting results demonstrate that bona fide phenotypic markers to distinguish leukemia initiating cells (LICs) in ALL are not yet available.
The metabolic activity of ALL is also poorly characterized. Previous reports have shown that ALL cells rely on glycolysis as their energy resource (the so-called Warburg effect) and an increased glycolysis rate was shown to be directly related to resistance to glucocorticoid treatment, which caused negative clinical outcomes (15–20). Thus, glycolysis is thought to have essential roles in ALL survival and malignancy. However, HSCs in general are known to be metabolically quiescent, with different manifestations in normal versus malignant HSCs (21–23). While normal HSCs have decreased oxidative respiration and instead rely on glycolysis, malignant HSCs–for example stem cells in AML–have been shown to possess lower levels of both glycolysis and oxidative respiration compared to non-stem cells in AML (24). Moreover, cancer stem cells (CSCs) are known to demonstrate metabolic plasticity, for example oscillating between oxidative phosphorylation and glycolysis to best promote tumor growth depending on the environmental state (25, 26).
In summary, there is quite significant inconsistency in characterizing the hierarchical structure, phenotypic markers, and metabolic activity of LSCs or CSCs. In this paper, we identified a unique subpopulation of B-ALL cells with leukemia initiating capacity within human patient xenografts and cell lines. These cells possess very low glucose uptake and a distinctive molecular profile relative to the broader cell population while also displaying potent leukemia initiating capacity in vivo. We hypothesize that these cells could represent LICs.
Materials and methods
Reagents
2-(N-(7-Nitrobenz-2-oxa-1, 3-diazol-4-yl) Amino)-2-Deoxyglucose (NBDG) was purchased from Invitrogen.
Patient-derived leukemia cells
Primary patient leukemia samples were collected from patients with informed consent based on the UC Davis Institutional Review Board approved protocol and transplanted into 6–12-week-old female NOD/SCID/IL2Rg-/- (NSG) mice using our institutionally-approved animal care protocol. As mice developed leukemia, they were sacrificed and leukemia cells were harvested from the leukemia-infiltrated spleen, bone marrow or liver for experiments. Human leukemia cells were confirmed by flow cytometry using anti-HLA-ABC and anti-CD10, CD19, CD22, CD34, CD38, CD45 antibodies were used for phenotyping (Biolegends & BD Biosciences).
Cell isolation and sorting
Leukemia cells were resuspended in glucose-free DMEM and incubated with NBDG at 37°C in the dark for 30 minutes at 1ul of 5mg/ml NBDG per 3 million cells (27). After 30 minutes of incubation, leukemic cells were washed with PBS. The cells were then stained with anti-HLA antibody and 4’, 6-Diamidino-2-Phenylindole, Dilactate (DAPI) (Thermo Fisher Scientific). NBDG-low, NBDG-med, or NBDG-high cells were sorted by Cytopeia InFlux Cell Sorter (BD Biosciences) and BD Influx Cell Sorter (BD Biosciences) using 488nm excitation and 505-522nm emission wavelengths.
Leukemia cell transplantation
The human leukemia mouse model was created by intratibial cell injection. The cell numbers transplanted ranged from 100 to 2.5 million cells per mouse. Sorted cells were used for some in vivo studies at the indicated cell doses. To determine the differentiation capacity of NBDG-low cells in xenograft mice, phenotypes of engrafted leukemia cells were compared with those of the original patient sample and/or serially-transplanted xenografts. Serial transplantations were established using the same technique and leukemia cells harvested from bone marrow, leukemia-infiltrated spleen, and liver (up to fourth generation mice).
Mice were monitored daily and euthanized when they showed signs of sickness, such as unkempt fur, ataxia, or weight loss more than 20% of pre-treatment body weight, in accordance with Institutional Animal Care and Use Committee (IACUC) policy on Humane Endpoints. Healthy mice were observed up to several months when they were euthanized to check whether they developed leukemia.
Statistical significance for survival time (days) was determined by the log-rank test. Kaplan-Meier survival curves were plotted for the two groups. Mice that did not develop human leukemia were censored. Two mice that had signs of sickness, but negative expression of HLA, were excluded from the study. Analyses used Prism 8.3 software (GraphPad).
RNA-Sequencing library preparation and next-generation sequencing
Indexed, directional sequencing libraries were prepared directly from intact cells using the SMART-Seq Stranded Kit (Takara Bio, Inc.) according to the manufacturer’s standard protocol. Briefly, 300 cells were lysed and followed by random-primed first-strand cDNA synthesis, tailing, and template switching with reverse transcriptase (SMARTScribe RT). Illumina barcoded adapters were incorporated by PCR with SeqAmp DNA Polymerase and rRNA-derived cDNA was selectively depleted. The remaining mRNA-derived cDNA fragments were PCR-enriched and the resulting libraries purified with AMPure XP beads (Beckman Coulter, Indianapolis, IN). The RNA-Seq libraries were quantified using a Qubit dsDNA High Sensitivity Assay (Thermo Fisher Scientific) and sized by analysis with an Agilent 2100 Bioanalyzer (Agilent Technologies, Santa Clara, CA) using the Agilent High Sensitivity DNA Kit. Libraries were pooled and multiplex sequenced (2 x 150bp, paired-end) on an Illumina NovaSeq 6000 Sequencing System.
RNA-Seq data analysis
De-multiplexed raw sequencing data (FASTQ format) was evaluated using FastQC (Babraham Bioinformatics) (28) and passed standard quality metrics, such as per base quality scores with ≥ 85% of bases >Q30. Additional NGS data quality control, pre-processing, alignment, quantification, differential expression, and statistical analysis were performed with Partek Flow Genomic Analysis software (Partek Inc., Chesterfield, MO). Sequence reads were pre-processed to trim 3’ adapter sequences and SMART-Seq Stranded library-specific sequence (i.e., first three bases of read two). Reads were then mapped to the GRCh38/hg38 human reference genome assembly using the STAR (Spliced Transcripts Alignment to a Reference; version 2.7.8a) (29) aligner. Transcript expression quantification was performed using the quantify to annotation model (Partek E/M) with GENCODE Release 36 (GRCh38.p13) annotation. Differential expression analysis (e.g., group-wise comparison of NBDG-low vs. NBDG-med) was performed with DESeq2 (30) using the poscounts estimation method. Hierarchical clustering of differentially-expressed transcripts (DETs) was performed using average linkage cluster and Euclidean point distance metrics. Functional annotation enrichment analysis was performed on DETs using the ToppFun tool in the ToppGene Suite (31). Benjamini-Hochberg’s FDRs or Storey’s q-values were used to correct for multiple testing.
Results
A distinct NBDG-low population was identified in B-ALL cell lines and patient samples
To address the hypothesis that we could distinctly characterize B-ALL cells on the basis of metabolic features, we utilized the fluorescent glucose analog 2-(N-(7-Nitrobenz-2-oxa-1, 3-diazol-4-yl) Amino)-2-Deoxyglucose (NBDG) and distinguished low- and high-fluorescing cells with flow cytometry. These cells were referred to as NBDG-low (low glucose uptake) and NBDG-high (high glucose uptake) cells, respectively (32). A distinct population of NBDG-low cells was detected in the cell lines JM1 and Reh (Supplementary Figure 1) as well as in patient-derived B-ALL cells. The median rate of NBDG-low cell population was 1.92% (ranging from 0.36 to 5.75%) of the whole leukemia population, calculated by averaging the NBDG-low populations as a percentage of total leukemia cells, using 16 flow cytometry sorted samples (Supplementary Table 1). NBDG-low cells were smaller in size, ranging from 4.0 to 7.0 μm, compared to the counterpart NBDG-high cells, ranging from 5.0 to 14.0 μm (Figures 1A-C, and Supplementary Figure 2A). The NBDG-low cells also showed lower HLA expression than the NBDG-high cells (Figure 1A).
Figure 1
NBDG-low cells isolated from patient samples exhibit potent in vivo leukemia initiating capacity
Leukemia initiating capacity for LICs is defined based on a transplantation assay, which tests their capacity to initiate, propagate, and maintain bulk leukemia growth in vivo as a xenograft in NSG mice. To confirm our hypothesis that NBDG-low cells have leukemia initiating capacity, we isolated NBDG-low or NBDG-high cells from six patient samples (PS) and transplanted the cells into NSG mice. Sample information is summarized in Supplementary Table 2 (PS1-5 and PS10). All six patients were diagnosed with B-ALL; three patients were considered standard risk and three as high risk based on the currently used risk categorization (33). All the mice transplanted with NBDG-low cells developed leukemia between five and 14 weeks, whereas those transplanted with NBDG-high cells did not develop leukemia in the observed period, which was extended to more than three to four months after leukemia development in the NBDG-low group. The median survival time of the mice in the NBDG-low group was PS1: 36.0 days, PS2: 53.5 days, PS3: 98.5 days and PS4: 82.0 days. Mice in the NBDG-high groups of these four samples survived significantly longer than those of NBDG-low groups (p ≤ 0.0001-0.002) (Figure 1D).
The potency of NBDG-low cells for in vivo LIC was next examined by transplanting a very low number of cells (100 cells per mouse) using two samples, PS4 and PS5. In PS4, eight out of 14 mice transplanted with NBDG-low cells developed leukemia, whereas none of the 14 mice transplanted with NBDG-high cells developed leukemia by the end of the study. In PS5, six out of 11 mice transplanted with NBDG-low cells develop leukemia by the end of the study. However, none of the 12 mice transplanted with NBDG-high cells developed leukemia except that one mouse was found to have a partial leukemia phenotype at the end of the study. The median survival time of mice in the NBDG-low group is PS4: 263.5 days and PS5: 217.0 days. Mice in the NBDG-high groups of these two samples survived significantly longer than those in the NBDG-low groups (p=0.0009-0.02) (Supplementary Figure 2B).
Of note, cells with an intermediate level of NBDG uptake (NBDG-med), located between the NBDG-low and NBDG-high populations (Figure 2A), from three patient samples (PS1, PS2, PS6) were transplanted and these mice subsequently also developed leukemia. Survival curves compared to NBDG-low transplanted mice are shown in Figure 2B. Combined median survival for the NBDG-med mice (n=14) from the three samples was 52 days, compared to 51.5 days for the NBDG-low mice (n=14). One mouse in the NBDG-med group and two mice in the NBDG-low group showed no signs of leukemia at the end of the experiment (247 days). These results indicate that NBDG-med cells were also able to establish an in vivo leukemia phenotype, similar to NBDG-low cells.
Figure 2
NBDG-low cells show in vivo differentiation and self-renewal capacity
Overall, the leukemia developed from transplanted NBDG-low cells recapitulate the original B-ALL phenotypes (results of PS4 are shown in Supplementary Table 3). These results indicated that NBDG-low cells were able to reconstitute and reestablish the complete leukemic phenotype in vivo, similar to the original leukemia phenotype. We further confirmed the self-renewal capacity of the NBDG-low cells by serial transplantation of the PS1 xenograft model. All the serially-transplanted mice successfully developed leukemia (n=9, three independent experiments). Altogether, these results demonstrated the leukemia initiating capacity of NBDG-low cells in six different B-ALL xenografts (using five patient samples).
NBDG-low cells have a unique transcriptome profile
As shown above, we identified a subpopulation of B-ALL cells from patient-derived xenograft (PDX) models that is characterized by low glucose uptake (i.e., based on NBDG fluorescence) and potent leukemia initiating capacity. To better understand the mechanisms underlying this biologic property, comprehensive molecular characterization was performed using next-generation RNA-sequencing (RNA-Seq)-based transcriptome analysis. This was performed on leukemia samples from five different B-ALL PDX models that were flow sorted into subpopulations based on NBDG fluorescence intensity, and therefore glucose uptake: NBDG-low (Low), NBDG-med (Medium), and NBDG-high (High). The RNA-Seq data was processed with a STAR-Partek(E/M)-DESeq2 analysis pipeline for read alignment, annotation, transcript expression estimation, and differential expression (DE) analysis as described in Materials and Methods.
As shown in Figure 3A, unsupervised analysis of the normalized gene/transcript count data using principal component analysis (PCA) demonstrated that the data from the NBDG-low samples possessed sufficient variation from that of the NBDG-med and NBDG-high samples to distinguish them from the latter groups in the 3D plot and similarity in expression. Notably, the NBDG-med and NBDG-high uptake samples demonstrated tighter intra-group similarity of expression profiles while the NBDG-low sample group appeared to be more heterogeneous.
Figure 3
To identify genes/transcripts whose altered expression were associated with and/or mediated the leukemia initiating phenotype, differential expression analysis was conducted by group-wise comparison of the NBDG-low and NBDG-med samples. The results revealed an expression profile of 2,162 DETs (NBDG-low vs. NBDG-med, p<0.05) consisting of 895 up-regulated and 1,267 down-regulated transcripts (Supplementary Table 4). Higher stringency statistical filtering resulted in 588 DETs (NBDG-low vs. NBDG-med, false discovery rate [FDR] p<0.05) (Supplementary Table 5) with a stronger bias towards down-regulated genes (n=424, 72.1%). Subsequent hierarchical clustering of these DETs and heatmap visualization demonstrated robust clustering into two major groups of samples: 1) NBDG-low and 2) NBDG-med with NBDG-high. In addition, there were two major patterns, and sub-clusters, of gene expression highlighted by marked uniformity in the expression of down-regulated genes in the NBDG-low samples and heterogeneous patterns of relative increased expression of genes in both the NBDG-med and NBDG-high samples (Supplementary Figures 3A, B). Similar group comparison of NBDG-low and NBDG-high samples revealed 1,104 DETs, with 656 and 448 being up- or down-regulated (p<0.05), respectively. When combined (i.e., “NBDG-low vs. NBDG-med” and “NBDG-low vs. NBDG-high” DETs), there were a total of 2,512 DETs, with 754 of these transcripts having commonly altered expression compared to that of both the NBDG-med and NBDG-high groups. In this instance, hierarchical clustering of the combined set of DETs and subsequent heatmap visualization demonstrated distinct segregation of the samples into 1) NBDG-low and 2) NBDG-med plus NBDG-high groups under two major branches of the dendrogram (Figure 3B). Additionally, two general clusters of increased and decreased gene/transcript expression were apparent for each group of samples, as well as the heterogeneity in up-regulated gene expression in each sample. Consistent with the striking leukemia-initiating phenotype of the NBDG-low subpopulations, these cells also exhibited a number of marked transcript expression alterations in both directions (e.g., log2 ratios of -19.93 and 11.10) (Figure 4; Supplementary Tables 4, 5).
Figure 4
Insight into the potentially altered processes and pathways in the NBDG-low cells was gained by conducting functional annotation enrichment analysis of the DETs comprising the leukemia-initiating phenotype transcript profile using the ToppFun tool (31). The analysis yields over-represented genes in functional categories and gene sets from across multiple databases and collections, such as Gene Ontology Resource (34), KEGG (35), Reactome (36), and the Molecular Signatures Database (MSigDB) (37). The complete set of functional enrichment results for the NBDG-low transcript expression profile is provided as Supplementary Table 6 and representative functional categories are presented in Figure 5. Notably, multiple functional categories of RNA metabolic processes, including RNA processing and splicing, followed by chromatin organization, are highly enriched in NBDG-low cells.
Figure 5
Furthermore, we investigated the metabolism-associated transcriptomic profile of the NBDG-low cells as these cells do not uptake glucose. Pathway analysis indicated that glucose metabolic steps, such as the TCA cycle and oxidative phosphorylation, are negatively regulated in NBDG-low cells. (Supplementary Figure 6).
Discussion
LSCs are thought to be the root of cancer and are responsible for treatment resistance and disease relapse. However, LSCs have not yet been defined in B-ALL and no study that focuses on identifying B-ALL LSC populations by glucose metabolism has been reported yet. In this study, our group discovered a unique subpopulation of B-ALL cells by glucose uptake from primary B-ALL samples. We demonstrated that this cell population has in vivo leukemia initiating capacity and a distinct transcriptome profile. As detailed below, this population seems quiescent yet becomes activated in the appropriate in vivo environment, evidencing significant heterogeneity in leukemia initiation.
Our study showed that NBDG-low cells represent a small percentage of the entire leukemia population (median 1.92%) and are smaller in size (4.0 to 7.0 μm), which is consistent with phenotypes of quiescent stem cells (22, 38–43) (Figure 1B, Figure 1C, and Supplementary Figure 2A). The frequency of HSCs and AML LSCs has been reported to be 0.001-0.1% (44–48) and 0.00002-0.01% (10, 11), respectively. The frequency of CSCs, including those for colon cancer, neck tumors and breast cancer, has been reported to vary from 0.6 to 50%, depending on the type of cancer (49). Our finding of NBDG-low cells in the range 0.36 to 5.75%, with potentially only a fraction of these cells representing true cancer stem cells, seems to be more comparable to stem cells from hematopoietic origin.
There are two key functions that define stem cells: proliferation/differentiation and self-renewal (13, 14). The isolated NBDG-low cells consistently possess proliferative capacity in the six different series of PDX mouse models tested in this study (Figure 1D, Figure 2B, and Supplementary Figure 2B). The NBDG-low cells also retained the capacity to differentiate into leukemia phenotypically similar to that of the original patient leukemia as demonstrated by the expression of common cell surface markers, such as CD10 and CD19 (an example shown in Supplementary Table 3). Cell differentiation capacity, in that the NBDG-low cells form the same leukemia as the original type, was also demonstrated (an example shown in Supplementary Table 3). In addition, the self-renewal capacity of isolated NBDG-low cells was shown by successful serial transplantation of the isolated cells.
We additionally found that the NBDG-med population, similar to NBDG-low, has in vivo leukemia initiating capacity. As NBDG-low cells are quiescent yet possess leukemia initiating capacity, we speculate that these cells become activated into the NBDG-med population in vivo. This population may then terminally differentiate into the NBDG-high population and lose its leukemia initiating capacity, as well as potentially de-differentiate back into the NBDG-low population. This notion is supported by our transcriptome profiling showing that NBDG-low cells have a distinct profile, whereas the NBDG-med and NBDG-high populations cluster together (Figure 3B) (13). NBDG-low cells also cluster separately from NBDG-high and unsorted (Whole) cell populations using high throughput quantitative RT-PCR analysis (Supplementary Figure 4), establishing the NBDG-low profile as distinct from all other populations. Moreover, in vitro colony forming assays were performed using Reh and JM1 B-ALL cell lines, resulting in colony formation between 6-36 days in NBDG-low cells versus 5-25 days in NBDG-high cells (Supplementary Figure 5). The delayed colony formation observed in NBDG-low cells further supports the quiescent behavior of NBDG-low prior to activation and differentiation. Altogether, our findings suggest that leukemia initiation is a highly fluid process that involves requisite circumstances for activation and the capacity for continuous re- and de-differentiation.
Pathway analysis of glucose metabolism shows negatively regulated glucose metabolism process in NBDG-low cells (Supplementary Figure 6), which is consistent with the low glucose uptake phenotype of these cells. Of note, part of fat metabolism is upregulated in the NBDG-low cells, suggesting these cells may use fat for their energy source. Previous studies show that lipid metabolism is closely associated with maintenance of stemness in cancer stem cells (50, 51), indicating that LICs also utilize the lipidic metabolism pathway to maintain their properties, such as stemness and non-glucose dependency.
Interestingly, functional enrichment results show significant enrichment of genes associated with RNA processing, including mRNA splicing via spliceosome and the spliceosomal complex (Figure 5). Previous studies indicated that cancer cells utilize alternative splicing for their survival (52, 53). Interestingly, RNA processing alterations are common in the five NBDG-low samples, whereas PCA analysis of their transcriptome profiles are heterogeneous (Figure 3). This suggests that RNA processing changes commonly occur in NBDG-low cells affecting a variety of other pathways contributing to the heterogeneity of NBDG-low cells. Mutations in splicing have been shown to contribute to aberrant hematopoiesis and may be among the first mutations to occur in HSCs (54, 55). Recent studies have shown aberrant splicing to be highly prevalent in leukemia driver genes and to significantly contribute to leukemia pathogenesis in both AML and B-ALL (56, 57). In fact, splicing can also represent a critical target for future investigations on potential treatments and drug resistance. For example, spliceosome inhibitors including isoginkgetin and pladienolide B have been shown to be cytotoxic in both mouse and human tumor cell lines (58), supporting future in vivo studies on splicing inhibition in leukemia. Moreover, alteration of CD19 splicing and resultant impairment of recognition by CD19-specific T-cells has been shown to contribute to treatment resistance (59); our study provides a possible genetic explanation for this finding. These suggest that NBDG-low cells modify the transcribed RNAs through the RNA binding proteins as well as the gene expression by the formation of heterochromatin to maintain their quiescent and undifferentiated phenotype. Furthermore, there is enrichment of genes associated with protein translation in the NBDG-low group; this is consistent with the known effects of aberrant translation on the pathogenesis of cancer (60, 61). For example, assembly of the eukaryotic initiation factor 4F complex has been shown to be essential for translation of mRNAs that contribute to cell growth and inhibition of apoptosis (62). This represents an additional genetic factor potentially driving the leukemia initiating properties of the NBDG-low population, as well as a therapeutic target for future studies.
Genes regulating chromosome organization, such as SMC1, TP53, and PARP1 are also up-regulated in the NBDG-low, leukemia-initiating subpopulation (Supplementary Table 6). These genes have been reported to be associated with cell cycle regulation and maintenance of stem cell quiescence (63–65). Histone deacetylases (HDACs) play a major role in epigenetic gene regulation via histone modification and chromatin remodeling, and they maintain LSC survival in chronic myeloid leukemia (CML) and AML (66–69). HDAC5, up-regulated in the NBDG-low population, may lead to the maintenance of quiescence in leukemia initiating cells.
In summary, our findings demonstrate the robustness of utilizing NBDG to consistently and accurately identify a unique subpopulation of leukemia-initiating B-ALL cells based on glucose uptake. This suggests the future potential for using other indicators of metabolic activity to develop a more comprehensive characterization of LICs. Importantly, gene expression profiling of multiple B-ALL PDXs demonstrated that the NBDG-low population can be clearly distinguished from the NBDG-med, NBDG-high, and whole cell populations by the differential expression of a distinct set of NBDG-low signature genes. Therefore, the NBDG-low gene expression signature provides further insight into the biological mechanisms underlying their persistence and recurrence and represents a molecular tool for identifying a distinct subset of LICs in B-cell ALL PDXs.
Funding
This work was supported by research funding from the Keaton’s Child Cancer Alliance, The Hartwell Foundation, Hyundai Hope on Wheels, Academic Senate Research Grant, Children’s Miracle Network Grant, Cancer Research Coordinating Committee Grant, National Center for Advancing Translational Sciences, NIH through Grant #UL1 TR000002, CTSC-MCRTP, Mr. and Mrs. George Davis, Mr. and Mrs. Ray Diamondstone, and Salad Cosmo (Satake), the Walter and Loretta Tepper Memorial Cancer Research Fund (Tepper), and the St. Baldrick’s Foundation and American Society of Hematology (Lee). The UC Davis Comprehensive Cancer Center Genomics Shared Resource, Flow Cytometry Shared Resource and Biostatistics Shared Resource are supported by Cancer Center Support Grant (P30CA093373) from the National Cancer Institute.
Acknowledgments
The authors wish to thank Ms. Stephenie Y. Liu (Department of Pathology and Laboratory Medicine, Genomics Shared Resource, UC Davis Comprehensive Cancer Center) for her expert technical assistance for RNA sample processing.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Statements
Data availability statement
The datasets presented in this study can be found online in the National Center for Biotechnology Information Gene Expression Omnibus (NCBI GEO) repository with the accession number GSE206258.
Ethics statement
The studies involving human participants were reviewed and approved by UC Davis Institutional Review Board. Written informed consent to participate in this study was provided by the participants’ legal guardian/next of kin. The animal study was reviewed and approved by UC Davis Institutional Animal Care and Use Committee.
Author contributions
NN, CT and NS contributed to the conception and design of the study. AL, CD, SY, and MI conducted the mouse experiments. JD and BM assisted with cell sorting. HK, RD, JM, and CT performed the gene expression profiling and analyses. KK, YL and LB performed the statistical analysis. AL, HK, MI, and CT wrote sections of the manuscript. All authors contributed to the article and approved the submitted version
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fonc.2022.972323/full#supplementary-material
References
1
BuickRNMindenMDMcCullochEA. Self-renewal in culture of proliferative blast progenitor cells in acute myeloblastic leukemia. Blood. (1979) 54(1):95–104. doi: 10.1182/blood.V54.1.95.95
2
TsiftsoglouASPappasISVizirianakisIS. Mechanisms involved in the induced differentiation of leukemia cells. Pharmacol Ther (2003) 100(3):257–90. doi: 10.1016/j.pharmthera.2003.09.002
3
UdomsakdiCEavesCJSwolinBReidDSBarnettMJEavesAC. Rapid decline of chronic myeloid leukemic-cells in long-term culture due to a defect at the leukemic stem-cell level. P Natl Acad Sci USA (1992) 89(13):6192–6. doi: 10.1073/pnas.89.13.6192
4
AillesLEGerhardBHoggeDE. Detection and characterization of primitive malignant and normal progenitors in patients with acute myelogenous leukemia using long-term coculture with supportive feeder layers and cytokines. Blood. (1997) 90(7):2555–64. doi: 10.1182/blood.V90.7.2555
5
LangFWojcikBRiegerMA. Stem cell hierarchy and clonal evolution in acute lymphoblastic leukemia. Stem Cells Int (2015) 2015:137164. doi: 10.1155/2015/137164
6
CobaledaCSánchez-GarcíaI. B-cell acute lymphoblastic leukaemia: towards understanding its cellular origin. BioEssays. (2009) 31(6):600–9. doi: 10.1002/bies.200800234
7
BerntKMArmstrongSA. Leukemia stem cells and human acute lymphoblastic leukemia. Semin Hematology. (2009) 46(1):33–8. doi: 10.1053/j.seminhematol.2008.09.010
8
MillerJSMcCullarVPunzelMLemischkaIRMooreKA. Single adult human CD34(+)/Lin(-)/CD38(-) progenitors give rise to natural killer cells, b-lineage cells, dendritic cells, and myeloid cells. Blood. (1999) 93(1):96–106. doi: 10.1182/blood.V93.1.96
9
HenonPSovalatHBeckerMArkamYOjeda-UribeMRaidotJPet al. Primordial role of CD34(+)38(-) cells in early and late trilineage haemopoietic engraftment after autologous blood cell transplantation. Brit J Haematol (1998) 103(2):568–81. doi: 10.1046/j.1365-2141.1998.01066.x
10
LapidotTSirardCVormoorJMurdochBHoangTCacerescortesJet al. A cell initiating human acute myeloid-leukemia after transplantation into scid mice. Nature. (1994) 367(6464):645–8. doi: 10.1038/367645a0
11
BonnetDDickJE. Human acute myeloid leukemia is organized as a hierarchy that originates from a primitive hematopoietic cell. Nat Med (1997) 3(7):730–7. doi: 10.1038/nm0797-730
12
JiangZDengMWeiXYeWXiaoYLinSet al. Heterogeneity of CD34 and CD38 expression in acute b lymphoblastic leukemia cells is reversible and not hierarchically organized. J Hematol Oncol (2016) 9(1):94. doi: 10.1186/s13045-016-0310-1
13
le ViseurCHotfilderMBomkenSWilsonKRottgersSSchrauderAet al. In childhood acute lymphoblastic leukemia, blasts at different stages of immunophenotypic maturation have stem cell properties. Cancer Cell (2008) 14(1):47–58. doi: 10.1016/j.ccr.2008.05.015
14
Gutierrez-CiancaCFloresPLGonzalezGSSanchez-Garcia. A primitive hematopoietic cell is the target for the leukemic transformation in human Philadelphia-positive acute lymphoblastic leukemia (vol 95, pg 1007, 2000). Blood. (2000) 95(7):2211–. doi: 10.1182/blood.V95.3.1007.003k35_1007_1013
15
KimJWDangCV. Cancer's molecular sweet tooth and the warburg effect. Cancer Res (2006) 66(18):8927–30. doi: 10.1158/0008-5472.CAN-06-1501
16
BoagJMBeesleyAHFirthMJFreitasJRFordJHoffmannKet al. Altered glucose metabolism in childhood pre-b acute lymphoblastic leukaemia. Leukemia. (2006) 20(10):1731–7. doi: 10.1038/sj.leu.2404365
17
HullemanEKazemierKMHollemanAVanderWeeleDJRudinCMBroekhuisMJet al. Inhibition of glycolysis modulates prednisolone resistance in acute lymphoblastic leukemia cells. Blood. (2009) 113(9):2014–21. doi: 10.1182/blood-2008-05-157842
18
HongoTYajimaSSakuraiMHorikoshiYHanadaR. In vitro drug sensitivity testing can predict induction failure and early relapse of childhood acute lymphoblastic leukemia. Blood. (1997) 89(8):2959–65. doi: 10.1182/blood.V89.8.2959
19
KaspersGJPietersRVan ZantwijkCHVan WeringERvan der Does-Van Den BergAVeermanAJ. Prednisolone resistance in childhood acute lymphoblastic leukemia: vitro-vivo correlations and cross-resistance to other drugs. Blood. (1998) 92(1):259–66. doi: 10.1182/blood.V92.1.259.413k21_259_266
20
Den BoerMLHarmsDOPietersRKazemierKMGobelUKorholzDet al. Patient stratification based on prednisolone-vincristine-asparaginase resistance profiles in children with acute lymphoblastic leukemia. J Clin Oncol (2003) 21(17):3262–8. doi: 10.1200/JCO.2003.11.031
21
SignerRAMageeJASalicAMorrisonSJ. Haematopoietic stem cells require a highly regulated protein synthesis rate. Nature. (2014) 509(7498):49–54. doi: 10.1038/nature13035
22
van VelthovenCTJRandoTA. Stem cell quiescence: Dynamism, restraint, and cellular idling. Cell Stem Cell (2019) 24(2):213–25. doi: 10.1016/j.stem.2019.01.001
23
RietherCSchurchCMOchsenbeinAF. Regulation of hematopoietic and leukemic stem cells by the immune system. Cell Death Differ (2015) 22(2):187–98. doi: 10.1038/cdd.2014.89
24
LagadinouEDSachACallahanKRossiRMNeeringSJMinhajuddinMet al. BCL-2 inhibition targets oxidative phosphorylation and selectively eradicates quiescent human leukemia stem cells. Cell Stem Cell (2013) 12(3):329–41. doi: 10.1016/j.stem.2012.12.013
25
CipolleschiMGMarziISantiniRFredducciDVinciMCD'AmicoMet al. Hypoxia-resistant profile implies vulnerability of cancer stem cells to physiological agents, which suggests new therapeutic targets. Cell Cycle (2014) 13(2):268–78. doi: 10.4161/cc.27031
26
Peiris-PagesMMartinez-OutschoornUEPestellRGSotgiaFLisantiMP. Cancer stem cell metabolism. Breast Cancer Res (2016) 18(1):55. doi: 10.1186/s13058-016-0712-6
27
NitinNCarlsonALMuldoonTEl-NaggarAKGillenwaterARichards-KortumR. Molecular imaging of glucose uptake in oral neoplasia following topical application of fluorescently labeled deoxy-glucose. Int J Cancer. (2009) 124(11):2634–42. doi: 10.1002/ijc.24222
28
WingettSWAndrewsS. FastQ screen: A tool for multi-genome mapping and quality control. F1000Res. (2018) 7:1338. doi: 10.12688/f1000research.15931.1
29
DobinADavisCASchlesingerFDrenkowJZaleskiCJhaSet al. STAR: ultrafast universal RNA-seq aligner. Bioinformatics. (2013) 29(1):15–21. doi: 10.1093/bioinformatics/bts635
30
LoveMIHuberWAndersS. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol (2014) 15(12):1–21. doi: 10.1186/s13059-014-0550-8
31
ChenJBardesEEAronowBJJeggaAG. ToppGene suite for gene list enrichment analysis and candidate gene prioritization. Nucleic Acids Res (2009) 37(Web Server):W305–W11. doi: 10.1093/nar/gkp427
32
TeSlaaTTeitellMA. Techniques to monitor glycolysis. Method Enzymol (2014) 542:91–114. doi: 10.1016/B978-0-12-416618-9.00005-4
33
SchultzKRPullenDJSatherHNShusterJJDevidasMBorowitzMJet al. Risk- and response-based classification of childhood b-precursor acute lymphoblastic leukemia: a combined analysis of prognostic markers from the pediatric oncology group (POG) and children's cancer group (CCG). Blood. (2007) 109(3):926–35. doi: 10.1182/blood-2006-01-024729
34
CarbonSDouglassEGoodBMUnniDRHarrisNLMungallCJet al. The gene ontology resource: enriching a GOld mine. Nucleic Acids Res (2021) 49(D1):D325–D34. doi: 10.1093/nar/gkaa1113
35
KanehisaMSatoYKawashimaMFurumichiMTanabeM. KEGG as a reference resource for gene and protein annotation. Nucleic Acids Res (2016) 44(D1):D457–D62. doi: 10.1093/nar/gkv1070
36
GillespieMJassalBStephanRMilacicMRothfelsKSenff-RibeiroAet al. The reactome pathway knowledgebase 2022. Nucleic Acids Res (2022) 50(D1):D687–D92. doi: 10.1093/nar/gkab1028
37
SubramanianATamayoPMoothaVKMukherjeeSEbertBLGilletteMAet al. Gene set enrichment analysis: A knowledge-based approach for interpreting genome-wide expression profiles. Proc Natl Acad Sci (2005) 102(43):15545–50. doi: 10.1073/pnas.0506580102
38
ZismanovVChichkovVColangeloVJametSWangSSymeAet al. Phosphorylation of eIF2 alpha is a translational control mechanism regulating muscle stem cell quiescence and self-renewal. Cell Stem Cell (2016) 18(1):79–90. doi: 10.1016/j.stem.2015.09.020
39
Llorens-BobadillaEZhaoSBaserASaiz-CastroGZwadloKMartin-VillalbaA. Single-cell transcriptomics reveals a population of dormant neural stem cells that become activated upon brain injury. Cell Stem Cell (2015) 17(3):329–40. doi: 10.1016/j.stem.2015.07.002
40
BuszczakMSignerRAJMorrisonSJ. Cellular differences in protein synthesis regulate tissue homeostasis. Cell. (2014) 159(2):242–51. doi: 10.1016/j.cell.2014.09.016
41
LiuLInokiAFanKMaoFShiGJinXet al. ER associated degradation preserves hematopoietic stem cell quiescence and self-renewal by restricting mTOR activity. Blood. (2020) 136:2975–86. doi: 10.1182/blood.2020007975
42
ValcourtJRLemonsJMSHaleyEMKojimaMDemurenOOCollerHA. Staying alive metabolic adaptations to quiescence. Cell Cycle (2012) 11(9):1680–96. doi: 10.4161/cc.19879
43
GrayJVPetskoGAJohnstonGCRingeDSingerRAWerner-WashburneM. "Sleeping beauty": Quiescence in saccharomyces cerevisiae. Microbiol Mol Biol R (2004) 68(2):187–+. doi: 10.1128/MMBR.68.2.187-206.2004
44
HarrisonDE. Competitive repopulation: a new assay for long-term stem cell functional capacity. Blood. (1980) 55(1):77–81. doi: 10.1182/blood.V55.1.77.77
45
HarrisonDEJordanCTZhongRKAstleCM. Primitive hematopoietic stem-cells - direct assay of most productive populations by competitive repopulation with simple binomial, correlation and covariance calculations. Exp Hematol (1993) 21(2):206–19.
46
SzilvassySJHumphriesRKLansdorpPMEavesACEavesCJ. Quantitative assay for totipotent reconstituting hematopoietic stem-cells by a competitive repopulation strategy. P Natl Acad Sci USA (1990) 87(22):8736–40. doi: 10.1073/pnas.87.22.8736
47
HawleyRGRamezaniAHawleyTS. [7] - hematopoietic stem cells. In: KlimanskayaILanzaR, editors. Methods in enzymology, vol. 419 . Academic PressSan Diego, California, USA: Elsevier (2006). p. 149–79.
48
KurandaKVargaftigJde la RocherePDosquetCCharronDBardinFet al. Age-related changes in human hematopoietic stem/progenitor cells. Aging Cell (2011) 10(3):542–6. doi: 10.1111/j.1474-9726.2011.00675.x
49
BakerM. Cancer stem cells, becoming common. Nat Rep Stem Cells (2008). doi: 10.1038/stemcells.2008.153
50
ManciniRNotoAPisanuMEDe VitisCMaugeri-SaccàMCilibertoG. Metabolic features of cancer stem cells: the emerging role of lipid metabolism. Oncogene. (2018) 37(18):2367–78. doi: 10.1038/s41388-018-0141-3
51
LiuHZhangZSongLGaoJLiuY. Lipid metabolism of cancer stem cells. Oncol Lett (2022) 23(4):119. doi: 10.3892/ol.2022.13239
52
DavidCJChenMAssanahMCanollPManleyJL. HnRNP proteins controlled by c-myc deregulate pyruvate kinase mRNA splicing in cancer. Nature. (2010) 463(7279):364–8. doi: 10.1038/nature08697
53
ParonettoMPPassacantilliISetteC. Alternative splicing and cell survival: from tissue homeostasis to disease. Cell Death Differ (2016) 23(12):1919–29. doi: 10.1038/cdd.2016.91
54
TaylorJLeeSC. Mutations in spliceosome genes and therapeutic opportunities in myeloid malignancies. Genes Chromosomes Cancer. (2019) 58(12):889–902. doi: 10.1002/gcc.22784
55
YangHBeutlerBZhangD. Emerging roles of spliceosome in cancer and immunity. Protein Cell (2021) 13:559–579. doi: 10.1007/s13238-021-00856-5
56
OchiYOgawaS. Chromatin-spliceosome mutations in acute myeloid leukemia. Cancers (Basel). (2021) 13(6). doi: 10.3390/cancers13061232
57
BlackKLNaqviASAsnaniMHayerKEYangSYGillespieEet al. Aberrant splicing in b-cell acute lymphoblastic leukemia. Nucleic Acids Res (2018) 46(21):11357–69. doi: 10.1093/nar/gky946
58
VanzylEJSayedHBlackmoreABRickKRCFernandoPMcKayBC. The spliceosome inhibitors isoginkgetin and pladienolide b induce ATF3-dependent cell death. PLoS One (2020) 15(12):e0224953. doi: 10.1371/journal.pone.0224953
59
SotilloEBarrettDMBlackKLBagashevAOldridgeDWuGet al. Convergence of acquired mutations and alternative splicing of CD19 enables resistance to CART-19 immunotherapy. Cancer Discovery (2015) 5(12):1282–95. doi: 10.1158/2159-8290.CD-15-1020
60
SongPYangFJinHWangX. The regulation of protein translation and its implications for cancer. Signal Transduction Targeted Ther (2021) 6(1):1–68. doi: 10.1038/s41392-020-00444-9
61
SpilkaRErnstCMehtaAKHaybaeckJ. Eukaryotic translation initiation factors in cancer development and progression. Cancer Letters. (2013) 340(1):9–21. doi: 10.1016/j.canlet.2013.06.019
62
MartelliAMEvangelistiCChappellWAbramsSLBäseckeJStivalaFet al. Targeting the translational apparatus to improve leukemia therapy: roles of the PI3K/PTEN/Akt/mTOR pathway. Leukemia. (2011) 25(7):1064–79. doi: 10.1038/leu.2011.46
63
LiuYElfSEMiyataYSashidaGLiuYHHuangGet al. p53 regulates hematopoietic stem cell quiescence. Cell Stem Cell (2009) 4(1):37–48. doi: 10.1016/j.stem.2008.11.006
64
LiXSippleJPangQSDuW. Salidroside stimulates DNA repair enzyme parp-1 activity in mouse HSC maintenance. Blood. (2012) 119(18):4162–73. doi: 10.1182/blood-2011-10-387332
65
LiXErdenOLiLYeQDWilsonADuW. Binding to WGR domain by salidroside activates PARP1 and protects hematopoietic stem cells from oxidative stress. Antioxid Redox Sign. (2014) 20(12):1853–65. doi: 10.1089/ars.2013.5600
66
CeccacciEMinucciS. Inhibition of histone deacetylases in cancer therapy: lessons from leukaemia. Brit J Cancer. (2016) 114(6):605–11. doi: 10.1038/bjc.2016.36
67
BamoduOAKuoKTYuanLPChengWHLeeWHHoYSet al. HDAC inhibitor suppresses proliferation and tumorigenicity of drug resistant chronic myeloid leukemia stem cells through regulation of hsa-miR-196a targeting BCR/ABL1. Exp Cell Res (2018) 370(2):519–30. doi: 10.1016/j.yexcr.2018.07.017
68
Al BaghdadiTAbonourRBoswellHS. Novel combination treatments targeting chronic myeloid leukemia stem cells. Cl Lymph Myelom Leuk. (2012) 12(2):94–105. doi: 10.1016/j.clml.2011.10.003
69
LiYChenKZhouYXiaoYDengMJiangZet al. A new strategy to target acute myeloid leukemia stem and progenitor cells using chidamide, a histone 1z deacetylase inhibitor. Curr Cancer Drug Targets. (2015) 15(6):493–503. doi: 10.2174/156800961506150805153230
Summary
Keywords
B-ALL, glucose, metabolic activity, leukemia initiating capacity, leukemia initiating cells, transcriptome profiling
Citation
Lee AQ, Konishi H, Duong C, Yoshida S, Davis RR, Van Dyke JE, Ijiri M, McLaughlin B, Kim K, Li Y, Beckett L, Nitin N, McPherson JD, Tepper CG and Satake N (2022) A distinct subpopulation of leukemia initiating cells in acute precursor B lymphoblastic leukemia: quiescent phenotype and unique transcriptomic profile. Front. Oncol. 12:972323. doi: 10.3389/fonc.2022.972323
Received
18 June 2022
Accepted
24 August 2022
Published
21 September 2022
Volume
12 - 2022
Edited by
Natalia Baran, University of Texas MD Anderson Cancer Center, United States
Reviewed by
Elvin Wagenblast, Icahn School of Medicine at Mount Sinai, United States; Francesca Sacco, University of Rome Tor Vergata, Italy
Updates
Copyright
© 2022 Lee, Konishi, Duong, Yoshida, Davis, Van Dyke, Ijiri, McLaughlin, Kim, Li, Beckett, Nitin, McPherson, Tepper and Satake.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Noriko Satake, nsatake@ucdavis.edu; Clifford G. Tepper, cgtepper@ucdavis.edu
This article was submitted to Molecular and Cellular Oncology, a section of the journal Frontiers in Oncology
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.