Cell-level metadata are indispensable for documenting single-cell sequencing datasets.

Single-cell RNA sequencing (scRNA-seq) provides an unprecedented view of cellular diversity of biological systems. However, across the thousands of publications and datasets generated using this technology, we estimate that only a minority (<25%) of studies provide cell-level metadata information containing identified cell types and related findings of the published dataset. Metadata omission hinders reproduction, exploration, validation, and knowledge transfer and is a common problem across journals, data repositories, and publication dates. We encourage investigators, reviewers, journals, and data repositories to improve their standards and ensure proper documentation of these valuable datasets.

Fatal COVID-19 and Non-COVID-19 Acute Respiratory Distress Syndrome Is Associated with Incomplete Alveolar Type 1 Epithelial Cell Differentiation from the Transitional State without Fibrosis.

Acute respiratory distress syndrome (ARDS) due to coronavirus disease 2019 and other etiologies results from injury to the alveolar epithelial cell (AEC) barrier resulting in noncardiogenic pulmonary edema, which causes acute respiratory failure; recovery requires epithelial regeneration. During physiological regeneration in mice, type 2 AECs (AEC2s) proliferate, exit the cell cycle, transiently assume a transitional state, then differentiate into type 1 AECs (AEC1s); in humans, persistence of the transitional state is associated with pulmonary fibrosis. It is unknown whether transitional cells emerge and differentiate into AEC1s without fibrosis in human ARDS and why transitional cells differentiate into AEC1s during physiological regeneration but persist in fibrosis. We hypothesized that incomplete but ongoing AEC1 differentiation from transitional cells without fibrosis may underlie persistent barrier permeability and acute respiratory failure in ARDS. Immunostaining of postmortem ARDS lungs revealed abundant transitional cells without fibrosis. They were typically cuboidal or partially spread, sometimes flat, and occasionally expressed AEC1 markers. Immunostaining and/or single-cell RNA sequencing revealed that transitional cells in mouse models of physiological regeneration, ARDS, and fibrosis express markers of cell cycle exit but only in fibrosis express a specific senescence marker. Thus, in severe, fatal early ARDS, AEC1 differentiation from transitional cells is incomplete, underlying persistent barrier permeability and respiratory failure but ongoing without fibrosis; senescence of transitional cells may be associated with pulmonary fibrosis.

Simultaneous measurement of biochemical phenotypes and gene expression in single cells.

Methods to measure heterogeneity among cells are rapidly transforming our understanding of biology but are currently limited to molecular abundance measurements. We developed an approach to simultaneously measure biochemical activities and mRNA abundance in single cells to understand the heterogeneity of DNA repair activities across thousands of human lymphocytes, identifying known and novel cell-type-specific DNA repair phenotypes.

Recovery and analysis of transcriptome subsets from pooled single-cell RNA-seq libraries.

Single-cell RNA sequencing (scRNA-seq) methods generate sparse gene expression profiles for thousands of single cells in a single experiment. The information in these profiles is sufficient to classify cell types by distinct expression patterns but the high complexity of scRNA-seq libraries often prevents full characterization of transcriptomes from individual cells. To extract more focused gene expression information from scRNA-seq libraries, we developed a strategy to physically recover the DNA molecules comprising transcriptome subsets, enabling deeper interrogation of the isolated molecules by another round of DNA sequencing. We applied the method in cell-centric and gene-centric modes to isolate cDNA fragments from scRNA-seq libraries. First, we resampled the transcriptomes of rare, single megakaryocytes from a complex mixture of lymphocytes and analyzed them in a second round of DNA sequencing, yielding up to 20-fold greater sequencing depth per cell and increasing the number of genes detected per cell from a median of 1313 to 2002. We similarly isolated mRNAs from targeted T cells to improve the reconstruction of their VDJ-rearranged immune receptor mRNAs. Second, we isolated CD3D mRNA fragments expressed across cells in a scRNA-seq library prepared from a clonal T cell line, increasing the number of cells with detected CD3D expression from 59.7% to 100%. Transcriptome resampling is a general approach to recover targeted gene expression information from single-cell RNA sequencing libraries that enhances the utility of these costly experiments, and may be applicable to the targeted recovery of molecules from other single-cell assays.

Dynamic temperature-sensitive A-to-I RNA editing in the brain of a heterothermic mammal during hibernation.

RNA editing diversifies genomically encoded information to expand the complexity of the transcriptome. In ectothermic organisms, including Drosophila and Cephalopoda, where body temperature mirrors ambient temperature, decreases in environmental temperature lead to increases in A-to-I RNA editing and cause amino acid recoding events that are thought to be adaptive responses to temperature fluctuations. In contrast, endothermic mammals, including humans and mice, typically maintain a constant body temperature despite environmental changes. Here, A-to-I editing primarily targets repeat elements, rarely results in the recoding of amino acids, and plays a critical role in innate immune tolerance. Hibernating ground squirrels provide a unique opportunity to examine RNA editing in a heterothermic mammal whose body temperature varies over 30°C and can be maintained at 5°C for many days during torpor. We profiled the transcriptome in three brain regions at six physiological states to quantify RNA editing and determine whether cold-induced RNA editing modifies the transcriptome as a potential mechanism for neuroprotection at low temperature during hibernation. We identified 5165 A-to-I editing sites in 1205 genes with dynamically increased editing after prolonged cold exposure. The majority (99.6%) of the cold-increased editing sites are outside of previously annotated coding regions, 82.7% lie in SINE-derived repeats, and 12 sites are predicted to recode amino acids. Additionally, A-to-I editing frequencies increase with increasing cold-exposure, demonstrating that ADAR remains active during torpor. Our findings suggest that dynamic A-to-I editing at low body temperature may provide a neuroprotective mechanism to limit aberrant dsRNA accumulation during torpor in the mammalian hibernator.

Overexpression of pendrin in intercalated cells produces chloride-sensitive hypertension.

Inherited and acquired disorders that enhance the activity of transporters mediating renal tubular Na(+) reabsorption are well established causes of hypertension. It is unclear, however, whether primary activation of an Na(+)-independent chloride transporter in the kidney can also play a pathogenic role in this disease. Here, mice overexpressing the chloride transporter pendrin in intercalated cells of the distal nephron (Tg(B1-hPDS) mice) displayed increased renal absorption of chloride. Compared with normal mice, these transgenic mice exhibited a delayed increase in urinary NaCl and ultimately, developed hypertension when exposed to a high-salt diet. Administering the same sodium intake as NaHCO3 instead of NaCl did not significantly alter BP, indicating that the hypertension in the transgenic mice was chloride-sensitive. Moreover, excessive chloride absorption by pendrin drove parallel absorption of sodium through the epithelial sodium channel ENaC and the sodium-driven chloride/bicarbonate exchanger (Ndcbe), despite an appropriate downregulation of these sodium transporters in response to the expanded vascular volume and hypertension. In summary, chloride transport in the distal nephron can play a primary role in driving NaCl transport in this part of the kidney, and a primary abnormality in renal chloride transport can provoke arterial hypertension. Thus, we conclude that the chloride/bicarbonate exchanger pendrin plays a major role in controlling net NaCl absorption, thereby influencing BP under conditions of high salt intake.

Spatial transcriptomic analysis delineates epithelial and mesenchymal subpopulations and transition stages in childhood ependymoma.

BACKGROUND: The diverse cellular constituents of childhood brain tumor ependymoma, recently revealed by single cell RNA-sequencing, may underly therapeutic resistance. Here we use spatial transcriptomics to further advance our understanding of the tumor microenvironment, mapping cellular subpopulations to the tumor architecture of ependymoma posterior fossa subgroup A (PFA), the commonest and most deadly childhood ependymoma variant. METHODS: Spatial transcriptomics data from intact PFA sections was deconvoluted to resolve the histological arrangement of neoplastic and non-neoplastic cell types. Key findings were validated using immunohistochemistry, in vitro functional assays and outcome analysis in clinically-annotated PFA bulk transcriptomic data. RESULTS: PFA are comprised of epithelial and mesenchymal histological zones containing a diversity of cellular states, each zone including co-existing and spatially distinct undifferentiated progenitor-like cells; a quiescent mesenchymal zone population, and a second highly mitotic progenitor population that is restricted to hypercellular epithelial zones and that is more abundant in progressive tumors. We show that myeloid cell interaction is the leading cause of mesenchymal transition in PFA, occurring in zones spatially distinct from hypoxia-induced mesenchymal transition, and these distinct EMT-initiating processes were replicated using in vitro models of PFA. CONCLUSIONS: These insights demonstrate the utility of spatial transcriptomics to advance our understanding of ependymoma biology, revealing a clearer picture of the cellular constituents of PFA, their interactions and influence on tumor progression.

Integration of single-nuclei RNA-sequencing, spatial transcriptomics and histochemistry defines the complex microenvironment of NF1-associated plexiform neurofibromas.

Plexiform neurofibroma (PN) is a leading cause of morbidity in children with the genetic condition Neurofibromatosis Type 1 (NF1), often disfiguring or threatening vital structures. During formation of PN, a complex tumor microenvironment (TME) develops, with recruitment of neoplastic and non-neoplastic cell types being critical for growth and progression. Due to the cohesive cellularity of PN, single-cell RNA-sequencing is difficult and may result in a loss of detection of critical cellular subpopulations. To bypass this barrier, we performed single-nuclei RNA-sequencing (snRNA-seq) on 8 frozen PN samples, and integrated this with spatial transcriptomics (ST) in 4 PN samples and immunohistochemistry to provide morphological context to transcriptomic data. SnRNA-seq analysis definitively charted the heterogeneous cellular subpopulations in the PN TME, with the predominant fraction being fibroblast subtypes. PN showed a remarkable amount of inter-sample homogeneity regarding cellular subpopulation proportions despite being resected from a variety of anatomical locations. ST analysis identified distinct cellular subpopulations which were annotated using snRNA-seq data and correlated with histological features. Schwann cell/fibroblast interactions were identified by receptor/ligand interaction analysis demonstrating a high probability of Neurexin 1/Neuroligin 1 (NRXN1/NLGN1) receptor-ligand cross-talk predicted between fibroblasts and non-myelinated Schwann cells (NM-SC) and subtypes, respectively. We observed aberrant expression of NRXN1 and NLGN1 in our PN snRNA-seq data compared to a normal mouse sciatic nerve single-cell RNA-seq dataset. This pathway has never been described in PN and may indicate a clear and direct communication pathway between putative NM-SC cells of origin and surrounding fibroblasts, potentially driving disease progression. SnRNA-seq integrated with spatial transcriptomics advances our understanding of the complex cellular heterogeneity of PN TME and identify potential novel communication pathways that may drive disease progression, a finding that could provide translational therapy options for patients with these devastating tumors of childhood and early adulthood.

Significant increase of high-risk chromosome 1q gain and 6q loss at recurrence in posterior fossa group A ependymoma: A multicenter study.

BACKGROUND: Ependymoma (EPN) posterior fossa group A (PFA) has the highest rate of recurrence and the worst prognosis of all EPN molecular groups. At relapse, it is typically incurable even with re-resection and re-irradiation. The biology of recurrent PFA remains largely unknown; however, the increasing use of surgery at first recurrence has now provided access to clinical samples to facilitate a better understanding of this. METHODS: In this large longitudinal international multicenter study, we examined matched samples of primary and recurrent disease from PFA patients to investigate the biology of recurrence. RESULTS: DNA methylome derived copy number variants (CNVs) revealed large-scale chromosome gains and losses at recurrence in PFA. CNV changes were dominated by chromosome 1q gain and/or 6q loss, both previously identified as high-risk factors in PFA, which were present in 23% at presentation but increased to 61% at first recurrence. Multivariate survival analyses of this cohort showed that cases with 1q gain or 6q loss at first recurrence were significantly more likely to recur again. Predisposition to 1q+/6q- CNV changes at recurrence correlated with hypomethylation of heterochromatin-associated DNA at presentation. Cellular and molecular analyses revealed that 1q+/6q- PFA had significantly higher proportions of proliferative neuroepithelial undifferentiated progenitors and decreased differentiated neoplastic subpopulations. CONCLUSIONS: This study provides clinically and preclinically actionable insights into the biology of PFA recurrence. The hypomethylation predisposition signature in PFA is a potential risk-classifier for trial stratification. We show that the cellular heterogeneity of PFAs evolves largely because of genetic evolution of neoplastic cells.

Identification of HDAC4 as a potential therapeutic target and prognostic biomarker for ZFTA-fused ependymomas.

Members of the HDAC family are predictive biomarkers and regulate the tumorigenesis in several cancers. However, the role of these genes in the biology of intracranial ependymomas (EPNs) remains unexplored. Here, an analysis of eighteen HDACs genes in an EPN transcriptomic dataset, revealed significantly higher levels of HDAC4 in supratentorial ZFTA fusion (ST-ZFTA) compared with ST-YAP1 fusion and posterior fossa EPNs, while HDAC7 and SIRT2 were downregulated in ST-ZFTA. HDAC4 was also overexpressed in ST-ZFTA as measured by single-cell RNA-Seq, quantitative real time-polymerase chain reaction, and immunohistochemistry. Survival analyses showed a significantly worse outcome for EPNs with higher HDAC4 and SIRT1 mRNA levels. Ontology enrichment analysis showed an HDAC4-high signature consistent with viral processes while collagen-containing extracellular matrix and cell-cell junction were enriched in those with an HDAC4-low signature. Immune gene analysis demonstrated a correlation between HDAC4 expression and low levels of NK resting cells. Several small molecules compounds targeting HDAC4 and ABCG2, were predicted by in silico analysis to be effective against HDAC4-high ZFTA. Our results provide novel insights into the biology of the HDAC family in intracranial ependymomas and reveal HDAC4 as a prognostic marker and potential therapeutic target in ST-ZFTA.

Regulation of epithelial transitional states in murine and human pulmonary fibrosis.

Idiopathic pulmonary fibrosis (IPF) is a progressive scarring disease arising from impaired regeneration of the alveolar epithelium after injury. During regeneration, type 2 alveolar epithelial cells (AEC2s) assume a transitional state that upregulates multiple keratins and ultimately differentiate into AEC1s. In IPF, transitional AECs accumulate with ineffectual AEC1 differentiation. However, whether and how transitional cells cause fibrosis, whether keratins regulate transitional cell accumulation and fibrosis, and why transitional AECs and fibrosis resolve in mouse models but accumulate in IPF are unclear. Here, we show that human keratin 8 (KRT8) genetic variants were associated with IPF. Krt8-/- mice were protected from fibrosis and accumulation of the transitional state. Keratin 8 (K8) regulated the expression of macrophage chemokines and macrophage recruitment. Profibrotic macrophages and myofibroblasts promoted the accumulation of transitional AECs, establishing a K8-dependent positive feedback loop driving fibrogenesis. Finally, rare murine transitional AECs were highly senescent and basaloid and may not differentiate into AEC1s, recapitulating the aberrant basaloid state in human IPF. We conclude that transitional AECs induced and were maintained by fibrosis in a K8-dependent manner; in mice, most transitional cells and fibrosis resolved, whereas in human IPF, transitional AECs evolved into an aberrant basaloid state that persisted with progressive fibrosis.

Neoplastic and immune single-cell transcriptomics define subgroup-specific intra-tumoral heterogeneity of childhood medulloblastoma.

BACKGROUND: Medulloblastoma (MB) is a heterogeneous disease in which neoplastic cells and associated immune cells contribute to disease progression. We aimed to determine the influence of neoplastic and immune cell diversity on MB biology in patient samples and animal models. METHODS: To better characterize cellular heterogeneity in MB we used single-cell RNA sequencing, immunohistochemistry, and deconvolution of transcriptomic data to profile neoplastic and immune populations in patient samples and animal models across childhood MB subgroups. RESULTS: Neoplastic cells cluster primarily according to individual sample of origin which is influenced by chromosomal copy number variance. Harmony alignment reveals novel MB subgroup/subtype-associated subpopulations that recapitulate neurodevelopmental processes, including photoreceptor and glutamatergic neuron-like cells in molecular subgroups GP3 and GP4, and a specific nodule-associated neuronally differentiated subpopulation in the sonic hedgehog subgroup. We definitively chart the spectrum of MB immune cell infiltrates, which include subpopulations that recapitulate developmentally related neuron-pruning and antigen-presenting myeloid cells. MB cellular diversity matching human samples is mirrored in subgroup-specific mouse models of MB. CONCLUSIONS: These findings provide a clearer understanding of the diverse neoplastic and immune cell subpopulations that constitute the MB microenvironment.

Fatal COVID-19 and non-COVID-19 Acute Respiratory Distress Syndrome is Associated with Incomplete Alveolar Type 1 Epithelial Cell Differentiation from the Transitional State Without Fibrosis.

ARDS due to COVID-19 and other etiologies results from injury to the alveolar epithelial cell (AEC) barrier resulting in noncardiogenic pulmonary edema, which causes acute respiratory failure; clinical recovery requires epithelial regeneration. During physiologic regeneration in mice, AEC2s proliferate, exit the cell cycle, and transiently assume a transitional state before differentiating into AEC1s; persistence of the transitional state is associated with pulmonary fibrosis in humans. It is unknown whether transitional cells emerge and differentiate into AEC1s without fibrosis in human ARDS and why transitional cells differentiate into AEC1s during physiologic regeneration but persist in fibrosis. We hypothesized that incomplete but ongoing AEC1 differentiation from transitional cells without fibrosis may underlie persistent barrier permeability and fatal acute respiratory failure in ARDS. Immunostaining of postmortem ARDS lungs revealed abundant transitional cells in organized monolayers on alveolar septa without fibrosis. They were typically cuboidal or partially spread, sometimes flat, and occasionally expressed AEC1 markers. Immunostaining and/or interrogation of scRNAseq datasets revealed that transitional cells in mouse models of physiologic regeneration, ARDS, and fibrosis express markers of cell cycle exit but only in fibrosis express a specific senescence marker. Thus, in severe, fatal early ARDS, AEC1 differentiation from transitional cells is incomplete, underlying persistent barrier permeability and respiratory failure, but ongoing without fibrosis; senescence of transitional cells may be associated with pulmonary fibrosis.

Liver Transcriptome Dynamics During Hibernation Are Shaped by a Shifting Balance Between Transcription and RNA Stability.

Hibernators dramatically lower metabolism to save energy while fasting for months. Prolonged fasting challenges metabolic homeostasis, yet small-bodied hibernators emerge each spring ready to resume all aspects of active life, including immediate reproduction. The liver is the body's metabolic hub, processing and detoxifying macromolecules to provide essential fuels to brain, muscle and other organs throughout the body. Here we quantify changes in liver gene expression across several distinct physiological states of hibernation in 13-lined ground squirrels, using RNA-seq to measure the steady-state transcriptome and GRO-seq to measure transcription for the first time in a hibernator. Our data capture key timepoints in both the seasonal and torpor-arousal cycles of hibernation. Strong positive correlation between transcription and the transcriptome indicates that transcriptional control dominates the known seasonal reprogramming of metabolic gene expression in liver for hibernation. During the torpor-arousal cycle, however, discordance develops between transcription and the steady-state transcriptome by at least two mechanisms: 1) although not transcribed during torpor, some transcripts are unusually stable across the torpor bout; and 2) unexpectedly, on some genes, our data suggest continuing, slow elongation with a failure to terminate transcription across the torpor bout. While the steady-state RNAs corresponding to these read through transcripts did not increase during torpor, they did increase shortly after rewarming despite their simultaneously low transcription. Both of these mechanisms would assure the immediate availability of functional transcripts upon rewarming. Integration of transcriptional, post-transcriptional and RNA stability control mechanisms, all demonstrated in these data, likely initiate a serial gene expression program across the short euthermic period that restores the tissue and prepares the animal for the next bout of torpor.

ENTPD3 Marks Mature Stem Cell-Derived ß-Cells Formed by Self-Aggregation In Vitro.

Stem cell-derived ß-like cells (sBC) carry the promise of providing an abundant source of insulin-producing cells for use in cell replacement therapy for patients with diabetes, potentially allowing widespread implementation of a practical cure. To achieve their clinical promise, sBC need to function comparably with mature adult ß-cells, but as yet they display varying degrees of maturity. Indeed, detailed knowledge of the events resulting in human ß-cell maturation remains obscure. Here we show that sBC spontaneously self-enrich into discreet islet-like cap structures within in vitro cultures, independent of exogenous maturation conditions. Multiple complementary assays demonstrate that this process is accompanied by functional maturation of the self-enriched sBC (seBC); however, the seBC still contain distinct subpopulations displaying different maturation levels. Interestingly, the surface protein ENTPD3 (also known as nucleoside triphosphate diphosphohydrolase-3 [NDPTase3]) is a specific marker of the most mature seBC population and can be used for mature seBC identification and sorting. Our results illuminate critical aspects of in vitro sBC maturation and provide important insights toward developing functionally mature sBC for diabetes cell replacement therapy.

clustifyr: an R package for automated single-cell RNA sequencing cluster classification.

Assignment of cell types from single-cell RNA sequencing (scRNA-seq) data remains a time-consuming and error-prone process. Current packages for identity assignment use limited types of reference data and often have rigid data structure requirements. We developed the clustifyr R package to leverage several external data types, including gene expression profiles to assign likely cell types using data from scRNA-seq, bulk RNA-seq, microarray expression data, or signature gene lists. We benchmark various parameters of a correlation-based approach and implement gene list enrichment methods. clustifyr is a lightweight and effective cell-type assignment tool developed for compatibility with various scRNA-seq analysis workflows. clustifyr is publicly available at https://github.com/rnabioco/clustifyr.

Dynamic RNA Regulation in the Brain Underlies Physiological Plasticity in a Hibernating Mammal.

Hibernation is a physiological and behavioral phenotype that minimizes energy expenditure. Hibernators cycle between profound depression and rapid hyperactivation of multiple physiological processes, challenging our concept of mammalian homeostasis. How the hibernator orchestrates and survives these extremes while maintaining cell to organismal viability is unknown. Here, we enhance the genome integrity and annotation of a model hibernator, the 13-lined ground squirrel. Our new assembly brings this genome to near chromosome-level contiguity and adds thousands of previously unannotated genes. These new genomic resources were used to identify 6,505 hibernation-related, differentially-expressed and processed transcripts using RNA-seq data from three brain regions in animals whose physiological status was precisely defined using body temperature telemetry. A software tool, squirrelBox, was developed to foster further data analyses and visualization. SquirrelBox includes a comprehensive toolset for rapid visualization of gene level and cluster group dynamics, sequence scanning of k-mer and domains, and interactive exploration of gene lists. Using these new tools and data, we deconvolute seasonal from temperature-dependent effects on the brain transcriptome during hibernation for the first time, highlighting the importance of carefully timed samples for studies of differential gene expression in hibernation. The identified genes include a regulatory network of RNA binding proteins that are dynamic in hibernation along with the composition of the RNA pool. In addition to passive effects of temperature, we provide evidence for regulated transcription and RNA turnover during hibernation. Significant alternative splicing, largely temperature dependent, also occurs during hibernation. These findings form a crucial first step and provide a roadmap for future work toward defining novel mechanisms of tissue protection and metabolic depression that may 1 day be applied toward improving human health.

Single-Cell RNA Sequencing of Childhood Ependymoma Reveals Neoplastic Cell Subpopulations That Impact Molecular Classification and Etiology.

Ependymoma (EPN) is a brain tumor commonly presenting in childhood that remains fatal in most children. Intra-tumoral cellular heterogeneity in bulk-tumor samples significantly confounds our understanding of EPN biology, impeding development of effective therapy. We, therefore, use single-cell RNA sequencing, histology, and deconvolution to catalog cellular heterogeneity of the major childhood EPN subgroups. Analysis of PFA subgroup EPN reveals evidence of an undifferentiated progenitor subpopulation that either differentiates into subpopulations with ependymal cell characteristics or transitions into a mesenchymal subpopulation. Histological analysis reveals that progenitor and mesenchymal subpopulations co-localize in peri-necrotic zones. In conflict with current classification paradigms, relative PFA subpopulation proportions are shown to determine bulk-tumor-assigned subgroups. We provide an interactive online resource that facilitates exploration of the EPN single-cell dataset. This atlas of EPN cellular heterogeneity increases understanding of EPN biology.

Monocytic Subclones Confer Resistance to Venetoclax-Based Therapy in Patients with Acute Myeloid Leukemia.

Venetoclax-based therapy can induce responses in approximately 70% of older previously untreated patients with acute myeloid leukemia (AML). However, up-front resistance as well as relapse following initial response demonstrates the need for a deeper understanding of resistance mechanisms. In the present study, we report that responses to venetoclax +azacitidine in patients with AML correlate closely with developmental stage, where phenotypically primitive AML is sensitive, but monocytic AML is more resistant. Mechanistically, resistant monocytic AML has a distinct transcriptomic profile, loses expression of venetoclax target BCL2, and relies on MCL1 to mediate oxidative phosphorylation and survival. This differential sensitivity drives a selective process in patients which favors the outgrowth of monocytic subpopulations at relapse. Based on these findings, we conclude that resistance to venetoclax + azacitidine can arise due to biological properties intrinsic to monocytic differentiation. We propose that optimal AML therapies should be designed so as to independently target AML subclones that may arise at differing stages of pathogenesis. SIGNIFICANCE: Identifying characteristics of patients who respond poorly to venetoclax-based therapy and devising alternative therapeutic strategies for such patients are important topics in AML. We show that venetoclax resistance can arise due to intrinsic molecular/metabolic properties of monocytic AML cells and that such properties can potentially be targeted with alternative strategies.