Interrogation of the Microenvironmental Landscape in Brain Tumors Reveals Disease-Specific Alterations of Immune Cells.

Brain malignancies encompass a range of primary and metastatic cancers, including low-grade and high-grade gliomas and brain metastases (BrMs) originating from diverse extracranial tumors. Our understanding of the brain tumor microenvironment (TME) remains limited, and it is unknown whether it is sculpted differentially by primary versus metastatic disease. We therefore comprehensively analyzed the brain TME landscape via flow cytometry, RNA sequencing, protein arrays, culture assays, and spatial tissue characterization. This revealed disease-specific enrichment of immune cells with pronounced differences in proportional abundance of tissue-resident microglia, infiltrating monocyte-derived macrophages, neutrophils, and T cells. These integrated analyses also uncovered multifaceted immune cell activation within brain malignancies entailing converging transcriptional trajectories while maintaining disease- and cell-type-specific programs. Given the interest in developing TME-targeted therapies for brain malignancies, this comprehensive resource of the immune landscape offers insights into possible strategies to overcome tumor-supporting TME properties and instead harness the TME to fight cancer.

Mapping genotypes to chromatin accessibility profiles in single cells.

In somatic tissue differentiation, chromatin accessibility changes govern priming and precursor commitment towards cellular fates1-3. Therefore, somatic mutations are likely to alter chromatin accessibility patterns, as they disrupt differentiation topologies leading to abnormal clonal outgrowth. However, defining the impact of somatic mutations on the epigenome in human samples is challenging due to admixed mutated and wild-type cells. Here, to chart how somatic mutations disrupt epigenetic landscapes in human clonal outgrowths, we developed genotyping of targeted loci with single-cell chromatin accessibility (GoT-ChA). This high-throughput platform links genotypes to chromatin accessibility at single-cell resolution across thousands of cells within a single assay. We applied GoT-ChA to CD34+ cells from patients with myeloproliferative neoplasms with JAK2V617F-mutated haematopoiesis. Differential accessibility analysis between wild-type and JAK2V617F-mutant progenitors revealed both cell-intrinsic and cell-state-specific shifts within mutant haematopoietic precursors, including cell-intrinsic pro-inflammatory signatures in haematopoietic stem cells, and a distinct profibrotic inflammatory chromatin landscape in megakaryocytic progenitors. Integration of mitochondrial genome profiling and cell-surface protein expression measurement allowed expansion of genotyping onto DOGMA-seq through imputation, enabling single-cell capture of genotypes, chromatin accessibility, RNA expression and cell-surface protein expression. Collectively, we show that the JAK2V617F mutation leads to epigenetic rewiring in a cell-intrinsic and cell type-specific manner, influencing inflammation states and differentiation trajectories. We envision that GoT-ChA will empower broad future investigations of the critical link between somatic mutations and epigenetic alterations across clonal populations in malignant and non-malignant contexts.

Single-cell mutation analysis of clonal evolution in myeloid malignancies.

Myeloid malignancies, including acute myeloid leukaemia (AML), arise from the expansion of haematopoietic stem and progenitor cells that acquire somatic mutations. Bulk molecular profiling has suggested that mutations are acquired in a stepwise fashion: mutant genes with high variant allele frequencies appear early in leukaemogenesis, and mutations with lower variant allele frequencies are thought to be acquired later1-3. Although bulk sequencing can provide information about leukaemia biology and prognosis, it cannot distinguish which mutations occur in the same clone(s), accurately measure clonal complexity, or definitively elucidate the order of mutations. To delineate the clonal framework of myeloid malignancies, we performed single-cell mutational profiling on 146 samples from 123 patients. Here we show that AML is dominated by a small number of clones, which frequently harbour co-occurring mutations in epigenetic regulators. Conversely, mutations in signalling genes often occur more than once in distinct subclones, consistent with increasing clonal diversity. We mapped clonal trajectories for each sample and uncovered combinations of mutations that synergized to promote clonal expansion and dominance. Finally, we combined protein expression with mutational analysis to map somatic genotype and clonal architecture with immunophenotype. Our findings provide insights into the pathogenesis of myeloid transformation and how clonal complexity evolves with disease progression.

CXCL8/CXCR2 signaling mediates bone marrow fibrosis and is a therapeutic target in myelofibrosis.

Proinflammatory signaling is a hallmark feature of human cancer, including in myeloproliferative neoplasms (MPNs), most notably myelofibrosis (MF). Dysregulated inflammatory signaling contributes to fibrotic progression in MF; however, the individual cytokine mediators elicited by malignant MPN cells to promote collagen-producing fibrosis and disease evolution are yet to be fully elucidated. Previously, we identified a critical role for combined constitutive JAK/STAT and aberrant NF-κB proinflammatory signaling in MF development. Using single-cell transcriptional and cytokine-secretion studies of primary cells from patients with MF and the human MPLW515L (hMPLW515L) murine model of MF, we extend our previous work and delineate the role of CXCL8/CXCR2 signaling in MF pathogenesis and bone marrow fibrosis progression. Hematopoietic stem/progenitor cells from patients with MF are enriched for a CXCL8/CXCR2 gene signature and display enhanced proliferation and fitness in response to an exogenous CXCL8 ligand in vitro. Genetic deletion of Cxcr2 in the hMPLW515L-adoptive transfer model abrogates fibrosis and extends overall survival, and pharmacologic inhibition of the CXCR1/2 pathway improves hematologic parameters, attenuates bone marrow fibrosis, and synergizes with JAK inhibitor therapy. Our mechanistic insights provide a rationale for therapeutic targeting of the CXCL8/CXCR2 pathway among patients with MF.

Microglia promote glioblastoma via mTOR-mediated immunosuppression of the tumour microenvironment.

Tumour-associated microglia/macrophages (TAM) are the most numerous non-neoplastic populations in the tumour microenvironment in glioblastoma multiforme (GBM), the most common malignant brain tumour in adulthood. The mTOR pathway, an important regulator of cell survival/proliferation, is upregulated in GBM, but little is known about the potential role of this pathway in TAM. Here, we show that GBM-initiating cells induce mTOR signalling in the microglia but not bone marrow-derived macrophages in both in vitro and in vivo GBM mouse models. mTOR-dependent regulation of STAT3 and NF-κB activity promotes an immunosuppressive microglial phenotype. This hinders effector T-cell infiltration, proliferation and immune reactivity, thereby contributing to tumour immune evasion and promoting tumour growth in mouse models. The translational value of our results is demonstrated in whole transcriptome datasets of human GBM and in a novel in vitro model, whereby expanded-potential stem cells (EPSC)-derived microglia-like cells are conditioned by syngeneic patient-derived GBM-initiating cells. These results raise the possibility that microglia could be the primary target of mTOR inhibition, rather than the intrinsic tumour cells in GBM.

Combined deletion of cathepsin protease family members reveals compensatory mechanisms in cancer.

Proteases are important for regulating multiple tumorigenic processes, including angiogenesis, tumor growth, and invasion. Elevated protease expression is associated with poor patient prognosis across numerous tumor types. Several multigene protease families have been implicated in cancer, including cysteine cathepsins. However, whether individual family members have unique roles or are functionally redundant remains poorly understood. Here we demonstrate stage-dependent effects of simultaneously deleting cathepsin B (CtsB) and CtsS in a murine pancreatic neuroendocrine tumor model. Early in tumorigenesis, the double knockout results in an additive reduction in angiogenic switching, whereas at late stages, several tumorigenic phenotypes are unexpectedly restored to wild-type levels. We identified CtsZ, which is predominantly supplied by tumor-associated macrophages, as the compensatory protease that regulates the acquired tumor-promoting functions of lesions deficient in both CtsB and CtsS. Thus, deletion of multiple cathepsins can lead to stage-dependent, compensatory mechanisms in the tumor microenvironment, which has potential implications for the clinical consideration of selective versus pan-family cathepsin inhibitors in cancer.

Plasmacytoid dendritic cell expansion defines a distinct subset of RUNX1-mutated acute myeloid leukemia.

Plasmacytoid dendritic cells (pDCs) are the principal natural type I interferon-producing dendritic cells. Neoplastic expansion of pDCs and pDC precursors leads to blastic plasmacytoid dendritic cell neoplasm (BPDCN), and clonal expansion of mature pDCs has been described in chronic myelomonocytic leukemia. The role of pDC expansion in acute myeloid leukemia (AML) is poorly studied. Here, we characterize patients with AML with pDC expansion (pDC-AML), which we observe in ∼5% of AML cases. pDC-AMLs often possess cross-lineage antigen expression and have adverse risk stratification with poor outcome. RUNX1 mutations are the most common somatic alterations in pDC-AML (>70%) and are much more common than in AML without pDC expansion and BPDCN. We demonstrate that pDCs are clonally related to, as well as originate from, leukemic blasts in pDC-AML. We further demonstrate that leukemic blasts from RUNX1-mutated AML upregulate a pDC transcriptional program, poising the cells toward pDC differentiation and expansion. Finally, tagraxofusp, a targeted therapy directed to CD123, reduces leukemic burden and eliminates pDCs in a patient-derived xenograft model. In conclusion, pDC-AML is characterized by a high frequency of RUNX1 mutations and increased expression of a pDC transcriptional program. CD123 targeting represents a potential treatment approach for pDC-AML.

Molecular predictors of immunophenotypic measurable residual disease clearance in acute myeloid leukemia.

Measurable residual disease (MRD) is a powerful prognostic factor in acute myeloid leukemia (AML). However, pre-treatment molecular predictors of immunophenotypic MRD clearance remain unclear. We analyzed a dataset of 211 patients with pre-treatment next-generation sequencing who received induction chemotherapy and had MRD assessed by serial immunophenotypic monitoring after induction, subsequent therapy, and allogeneic stem cell transplant (allo-SCT). Induction chemotherapy led to MRD- remission, MRD+ remission, and persistent disease in 35%, 27%, and 38% of patients, respectively. With subsequent therapy, 34% of patients with MRD+ and 26% of patients with persistent disease converted to MRD-. Mutations in CEBPA, NRAS, KRAS, and NPM1 predicted high rates of MRD- remission, while mutations in TP53, SF3B1, ASXL1, and RUNX1 and karyotypic abnormalities including inv (3), monosomy 5 or 7 predicted low rates of MRD- remission. Patients with fewer individual clones were more likely to achieve MRD- remission. Among 132 patients who underwent allo-SCT, outcomes were favorable whether patients achieved early MRD- after induction or later MRD- after subsequent therapy prior to allo-SCT. As MRD conversion with chemotherapy prior to allo-SCT is rarely achieved in patients with specific baseline mutational patterns and high clone numbers, upfront inclusion of these patients into clinical trials should be considered.

Genetic and epigenetic evolution as a contributor to WT1-mutant leukemogenesis.

Genetic studies have identified recurrent somatic mutations in acute myeloid leukemia (AML) patients, including in the Wilms' tumor 1 (WT1) gene. The molecular mechanisms by which WT1 mutations contribute to leukemogenesis have not yet been fully elucidated. We investigated the role of Wt1 gene dosage in steady-state and pathologic hematopoiesis. Wt1 heterozygous loss enhanced stem cell self-renewal in an age-dependent manner, which increased stem cell function over time and resulted in age-dependent leukemic transformation. Wt1-haploinsufficient leukemias were characterized by progressive genetic and epigenetic alterations, including those in known leukemia-associated alleles, demonstrating a requirement for additional events to promote hematopoietic transformation. Consistent with this observation, we found that Wt1 depletion cooperates with Flt3-ITD mutation to induce fully penetrant AML. Our studies provide insight into mechanisms of Wt1-loss leukemogenesis and into the evolutionary events required to induce transformation of Wt1-haploinsufficient stem/progenitor cells.

PURE-seq identifies Egr1 as a Potential Master Regulator in Murine Aging by Sequencing Long-Term Hematopoietic Stem Cells.

Single-cell transcriptomics is valuable for uncovering individual cell properties, particularly in highly heterogeneous systems. However, this technique often results in the analysis of many well-characterized cells, increasing costs and diluting rare cell populations. To address this, we developed PURE-seq (PIP-seq for Rare-cell Enrichment and Sequencing) for scalable sequencing of rare cells. PURE-seq allows direct cell loading from FACS into PIP-seq reactions, minimizing handling and reducing cell loss. PURE-seq reliably captures rare cells, with 60 minutes of sorting capturing tens of cells at a rarity of 1 in 1,000,000. Using PURE-seq, we investigated murine long-term hematopoietic stem cells and their transcriptomes in the context of hematopoietic aging, identifying Egr1 as a potential master regulator of hematopoiesis in the aging context. PURE-seq offers an accessible and reliable method for isolating and sequencing cells that are currently too rare to capture successfully with existing methods.

Mutation order in acute myeloid leukemia identifies uncommon patterns of evolution and illuminates phenotypic heterogeneity.

Acute myeloid leukemia (AML) has a poor prognosis and a heterogeneous mutation landscape. Although common mutations are well-studied, little research has characterized how the sequence of mutations relates to clinical features. Using published, single-cell DNA sequencing data from three institutions, we compared clonal evolution patterns in AML to patient characteristics, disease phenotype, and outcomes. Mutation trees, which represent the order of select mutations, were created for 207 patients from targeted panel sequencing data using 1 639 162 cells, 823 mutations, and 275 samples. In 224 distinct orderings of mutated genes, mutations related to DNA methylation typically preceded those related to cell signaling, but signaling-first cases did occur, and had higher peripheral cell counts, increased signaling mutation homozygosity, and younger patient age. Serial sample analysis suggested that NPM1 and DNA methylation mutations provide an advantage to signaling mutations in AML. Interestingly, WT1 mutation evolution shared features with signaling mutations, such as WT1-early being proliferative and occurring in younger individuals, trends that remained in multivariable regression. Some mutation orderings had a worse prognosis, but this was mediated by unfavorable mutations, not mutation order. These findings add a dimension to the mutation landscape of AML, identifying uncommon patterns of leukemogenesis and shedding light on heterogeneous phenotypes.

FLT3 tyrosine kinase inhibition modulates PRC2 and promotes differentiation in acute myeloid leukemia.

Internal tandem duplication mutations in fms-like tyrosine kinase 3 (FLT3-ITD) are recurrent in acute myeloid leukemia (AML) and increase the risk of relapse. Clinical responses to FLT3 inhibitors (FLT3i) include myeloid differentiation of the FLT3-ITD clone in nearly half of patients through an unknown mechanism. We identified enhancer of zeste homolog 2 (EZH2), a component of polycomb repressive complex 2 (PRC2), as a mediator of this effect using a proteomic-based screen. FLT3i downregulated EZH2 protein expression and PRC2 activity on H3K27me3. FLT3-ITD and loss-of-function mutations in EZH2 are mutually exclusive in human AML. We demonstrated that FLT3i increase myeloid maturation with reduced stem/progenitor cell populations in murine Flt3-ITD AML. Combining EZH1/2 inhibitors with FLT3i increased terminal maturation of leukemic cells and reduced leukemic burden. Our data suggest that reduced EZH2 activity following FLT3 inhibition promotes myeloid differentiation of FLT3-ITD leukemic cells, providing a mechanistic explanation for the clinical observations. These results demonstrate that in addition to its known cell survival and proliferation signaling, FLT3-ITD has a second, previously undefined function to maintain a myeloid stem/progenitor cell state through modulation of PRC2 activity. Our findings support exploring EZH1/2 inhibitors as therapy for FLT3-ITD AML.

CRISPR Dependency Screens in Primary Hematopoietic Stem Cells Identify KDM3B as a Genotype Specific Vulnerability in IDH2- and TET2-Mutant Cells.

Clonal hematopoiesis (CH) is a common premalignant state in the blood and confers an increased risk of blood cancers and all-cause mortality. Identification of therapeutic targets in CH has been hindered by the lack of an ex vivo platform amenable for studying primary hematopoietic stem and progenitor cells (HSPCs). Here, we utilize an ex vivo co-culture system of HSPCs with bone marrow endothelial cells to perform CRISPR/Cas9 screens in mutant HSPCs. Our data reveal that loss of the histone demethylase family members Kdm3b and Jmjd1c specifically reduces the fitness of Idh2- and Tet2-mutant HSPCs. Kdm3b loss in mutant cells leads to decreased expression of critical cytokine receptors including Mpl, rendering mutant HSPCs preferentially susceptible to inhibition of downstream JAK2 signaling. Our study nominates an epigenetic regulator and an epigenetically regulated receptor signaling pathway as genotype-specific therapeutic targets and provides a scalable platform to identify genetic dependencies in mutant HSPCs.

Jak2V617F Reversible Activation Shows Its Essential Requirement in Myeloproliferative Neoplasms.

Gain-of-function mutations activating JAK/STAT signaling are seen in the majority of patients with myeloproliferative neoplasms (MPN), most commonly JAK2V617F. Although clinically approved JAK inhibitors improve symptoms and outcomes in MPNs, remissions are rare, and mutant allele burden does not substantively change with chronic therapy. We hypothesized this is due to limitations of current JAK inhibitors to potently and specifically abrogate mutant JAK2 signaling. We therefore developed a conditionally inducible mouse model allowing for sequential activation, and then inactivation, of Jak2V617F from its endogenous locus using a combined Dre-rox/Cre-lox dual-recombinase system. Jak2V617F deletion abrogates MPN features, induces depletion of mutant-specific hematopoietic stem/progenitor cells, and extends overall survival to an extent not observed with pharmacologic JAK inhibition, including when cooccurring with somatic Tet2 loss. Our data suggest JAK2V617F represents the best therapeutic target in MPNs and demonstrate the therapeutic relevance of a dual-recombinase system to assess mutant-specific oncogenic dependencies in vivo. SIGNIFICANCE: Current JAK inhibitors to treat myeloproliferative neoplasms are ineffective at eradicating mutant cells. We developed an endogenously expressed Jak2V617F dual-recombinase knock-in/knock-out model to investigate Jak2V617F oncogenic reversion in vivo. Jak2V617F deletion abrogates MPN features and depletes disease-sustaining MPN stem cells, suggesting improved Jak2V617F targeting offers the potential for greater therapeutic efficacy. See related commentary by Celik and Challen, p. 701. This article is featured in Selected Articles from This Issue, p. 695.

Mutation order in acute myeloid leukemia identifies uncommon patterns of evolution and illuminates phenotypic heterogeneity.

Acute myeloid leukemia (AML) has a poor prognosis and a heterogeneous mutation landscape. Although common mutations are well-studied, little research has characterized how the sequence of mutations relates to clinical features. Using published, single-cell DNA sequencing data from three institutions, we compared clonal evolution patterns in AML to patient characteristics, disease phenotype, and outcomes. Mutation trees, which represent the order of select mutations, were created for 207 patients from targeted panel sequencing data using 1 639 162 cells, 823 mutations, and 275 samples. In 224 distinct orderings of mutated genes, mutations related to DNA methylation typically preceded those related to cell signaling, but signaling-first cases did occur, and had higher peripheral cell counts, increased signaling mutation homozygosity, and younger patient age. Serial sample analysis suggested that NPM1 and DNA methylation mutations provide an advantage to signaling mutations in AML. Interestingly, WT1 mutation evolution shared features with signaling mutations, such as WT1-early being proliferative and occurring in younger individuals, trends that remained in multivariable regression. Some mutation orderings had a worse prognosis, but this was mediated by unfavorable mutations, not mutation order. These findings add a dimension to the mutation landscape of AML, identifying uncommon patterns of leukemogenesis and shedding light on heterogenous phenotypes.

Single-cell genotypic and phenotypic analysis of measurable residual disease in acute myeloid leukemia.

Measurable residual disease (MRD), defined as the population of cancer cells that persist following therapy, serves as the critical reservoir for disease relapse in acute myeloid leukemia and other malignancies. Understanding the biology enabling MRD clones to resist therapy is necessary to guide the development of more effective curative treatments. Discriminating between residual leukemic clones, preleukemic clones, and normal precursors remains a challenge with current MRD tools. Here, we developed a single-cell MRD (scMRD) assay by combining flow cytometric enrichment of the targeted precursor/blast population with integrated single-cell DNA sequencing and immunophenotyping. Our scMRD assay shows high sensitivity of approximately 0.01%, deconvolutes clonal architecture, and provides clone-specific immunophenotypic data. In summary, our scMRD assay enhances MRD detection and simultaneously illuminates the clonal architecture of clonal hematopoiesis/preleukemic and leukemic cells surviving acute myeloid leukemia therapy.

Deep generative model deciphers derailed trajectories in acute myeloid leukemia.

Single-cell genomics has the potential to map cell states and their dynamics in an unbiased way in response to perturbations like disease. However, elucidating the cell-state transitions from healthy to disease requires analyzing data from perturbed samples jointly with unperturbed reference samples. Existing methods for integrating and jointly visualizing single-cell datasets from distinct contexts tend to remove key biological differences or do not correctly harmonize shared mechanisms. We present Decipher, a model that combines variational autoencoders with deep exponential families to reconstruct derailed trajectories (https://github.com/azizilab/decipher). Decipher jointly represents normal and perturbed single-cell RNA-seq datasets, revealing shared and disrupted dynamics. It further introduces a novel approach to visualize data, without the need for methods such as UMAP or TSNE. We demonstrate Decipher on data from acute myeloid leukemia patient bone marrow specimens, showing that it successfully characterizes the divergence from normal hematopoiesis and identifies transcriptional programs that become disrupted in each patient when they acquire NPM1 driver mutations.

Hematopoietic-specific heterozygous loss of Dnmt3a exacerbates colitis-associated colon cancer.

Clonal hematopoiesis (CH) is defined as clonal expansion of mutant hematopoietic stem cells absent diagnosis of a hematologic malignancy. Presence of CH in solid tumor patients, including colon cancer, correlates with shorter survival. We hypothesized that bone marrow-derived cells with heterozygous loss-of-function mutations of DNMT3A, the most common genetic alteration in CH, contribute to the pathogenesis of colon cancer. In a mouse model that combines colitis-associated colon cancer (CAC) with experimental CH driven by Dnmt3a+/Δ, we found higher tumor penetrance and increased tumor burden compared with controls. Histopathological analysis revealed accentuated colonic epithelium injury, dysplasia, and adenocarcinoma formation. Transcriptome profiling of colon tumors identified enrichment of gene signatures associated with carcinogenesis, including angiogenesis. Treatment with the angiogenesis inhibitor axitinib eliminated the colon tumor-promoting effect of experimental CH driven by Dnmt3a haploinsufficiency and rebalanced hematopoiesis. This study provides conceptually novel insights into non-tumor-cell-autonomous effects of hematopoietic alterations on colon carcinogenesis and identifies potential therapeutic strategies.

eIF2α kinases control chalone production in Dictyostelium discoideum.

Growing Dictyostelium cells secrete CfaD and AprA, two proteins that have been characterized as chalones. They exist within a high-molecular-weight complex that reversibly inhibits cell proliferation, but not growth, via cell surface receptors and a signaling pathway that includes G proteins. How the production of these two proteins is regulated is unknown. Dictyostelium cells possess three GCN2-type eukaryotic initiation factor 2 α subunit (eIF2α) kinases, proteins that phosphorylate the translational initiation factor eIF2α and possess a tRNA binding domain involved in their regulation. The Dictyostelium kinases have been shown to function during development in regulating several processes. We show here that expression of an unregulated, activated kinase domain greatly inhibits cell proliferation. The inhibitory effect on proliferation is not due to a general inhibition of translation. Instead, it is due to enhanced production of a secreted factor(s). Indeed, extracellular CfaD and AprA proteins, but not their mRNAs, are overproduced in cells expressing the activated kinase domain. The inhibition of proliferation is not seen when the activated kinase domain is expressed in cells lacking CfaD or AprA or in cells that contain a nonphosphorylatable eIF2α. We conclude that production of the chalones CfaD and AprA is translationally regulated by eIF2α phosphorylation. Both proteins are upregulated at the culmination of development, and this enhanced production is lacking in a strain that possesses a nonphosphorylatable eIF2α.

An integrated pipeline for comprehensive analysis of immune cells in human brain tumor clinical samples.

Human tissue samples represent an invaluable source of information for the analysis of disease-specific cellular alterations and their variation between different pathologies. In cancer research, advancing a comprehensive understanding of the unique characteristics of individual tumor types and their microenvironment is of considerable importance for clinical translation. However, investigating human brain tumor tissue is challenging due to the often-limited availability of surgical specimens. Here we describe a multimodule integrated pipeline for the processing of freshly resected human brain tumor tissue and matched blood that enables analysis of the tumor microenvironment, with a particular focus on the tumor immune microenvironment (TIME). The protocol maximizes the information yield from limited tissue and includes both the preservation of bulk tissue, which can be performed within 1 h following surgical resection, as well as tissue dissociation for an in-depth characterization of individual TIME cell populations, which typically takes several hours depending on tissue quantity and further downstream processing. We also describe integrated modules for immunofluorescent staining of sectioned tissue, bulk tissue genomic analysis and fluorescence- or magnetic-activated cell sorting of digested tissue for subsequent culture or transcriptomic analysis by RNA sequencing. Applying this pipeline, we have previously described the overall TIME landscape across different human brain malignancies, and were able to delineate disease-specific alterations of tissue-resident versus recruited macrophage populations. This protocol will enable researchers to use this pipeline to address further research questions regarding the tumor microenvironment.