Functional interactions between neurofibromatosis tumor suppressors underlie Schwann cell tumor de-differentiation and treatment resistance.

Schwann cell tumors are the most common cancers of the peripheral nervous system and can arise in patients with neurofibromatosis type-1 (NF-1) or neurofibromatosis type-2 (NF-2). Functional interactions between NF1 and NF2 and broader mechanisms underlying malignant transformation of the Schwann lineage are unclear. Here we integrate bulk and single-cell genomics, biochemistry, and pharmacology across human samples, cell lines, and mouse allografts to identify cellular de-differentiation mechanisms driving malignant transformation and treatment resistance. We find DNA methylation groups of Schwann cell tumors can be distinguished by differentiation programs that correlate with response to the MEK inhibitor selumetinib. Functional genomic screening in NF1-mutant tumor cells reveals NF2 loss and PAK activation underlie selumetinib resistance, and we find that concurrent MEK and PAK inhibition is effective in vivo. These data support a de-differentiation paradigm underlying malignant transformation and treatment resistance of Schwann cell tumors and elucidate a functional link between NF1 and NF2.

Molecular Features of Resected Melanoma Brain Metastases, Clinical Outcomes, and Responses to Immunotherapy.

Importance: Central nervous system (CNS)-penetrant systemic therapies have significantly advanced care for patients with melanoma brain metastases. However, improved understanding of the molecular landscape and microenvironment of these lesions is needed to both optimize patient selection and advance treatment approaches. Objective: To evaluate how bulk and single-cell genomic features of melanoma brain metastases are associated with clinical outcome and treatment response. Design, Setting, and Participants: This cohort study analyzed bulk DNA sequencing and single nuclear RNA-sequencing data from resected melanoma brain metastases and included 94 consecutive patients with a histopathologically confirmed diagnosis of melanoma brain metastasis who underwent surgical resection at a single National Comprehensive Cancer Network cancer center in San Francisco, California, from January 1, 2009, to December 31, 2022. Exposure: A Clinical Laboratory Improvement Amendments-certified targeted sequencing assay was used to analyze tumor resection specimens, with a focus on BRAF V600E alteration. For frozen pathologic specimens from CNS treatment-naive patients undergoing surgical resection, commercial single nuclear RNA sequencing approaches were used. Main Outcomes and Measures: The primary outcome was overall survival (OS). Secondary outcomes included CNS progression-free survival (PFS), microenvironmental composition with decreased T-cell and macrophage populations, and responses to immunotherapy. Results: To correlate molecular status with clinical outcome, Kaplan-Meier survival analysis of 94 consecutive patients (median age, 64 years [range, 24-82 years]; 70 men [74%]) with targeted BRAF alteration testing showed worse median intracranial PFS (BRAF variant: 3.6 months [IQR, 0.1-30.6 months]; BRAF wildtype: 11.0 months [IQR, 0.8-81.5 months]; P < .001) and OS (BRAF variant: 9.8 months [IQR, 2.5-69.4 months]; BRAF wildtype: 23.2 months [IQR, 1.1-102.5 months]; P = .005; log-rank test) in BRAF V600E variant tumors. Multivariable Cox proportional hazards regression analysis revealed that BRAF V600E status was an independent variable significantly associated with both PFS (hazard ratio [HR], 2.65; 95% CI, 1.54-4.57; P < .001) and OS (HR, 1.96; 95% CI, 1.08-3.55; P = .03). For the 45 patients with resected melanoma brain metastases undergoing targeted DNA sequencing, molecular classification recapitulated The Cancer Genome Atlas groups (NRAS variant, BRAF variant, NF1 variant, and triple wildtype) with no subtype enrichment within the brain metastasis cohort. On a molecular level, BRAF V600E variant lesions were found to have a significantly decreased tumor mutation burden. Moreover, single nuclear RNA sequencing of treatment-naive BRAF V600E variant (n = 3) brain metastases compared with BRAF wildtype (n = 3) brain metastases revealed increased immune cell populations in BRAF wildtype tumors (mean [SD], 11% [4.1%] vs 3% [1.6%] CD45-positive cells; P = .04). Survival analysis of postoperative immunotherapy responses by BRAF status revealed that BRAF wildtype lesions were associated with a response to checkpoint inhibition (median OS: with immunotherapy, undefined; without immunotherapy, 13.0 months [range, 1.1-61.7 months]; P = .001; log-rank test) while BRAF variant lesions (median OS: with immunotherapy, 9.8 months [range, 2.9-39.8 months]; without immunotherapy, 9.5 months [range, 2.5-67.2 months]; P = .81; log-rank test) were not. Conclusions and Relevance: This molecular analysis of patients with resected melanoma brain metastases found that BRAF V600E alteration is an important translational biomarker associated with worse clinical outcomes, differential microenvironmental composition, and benefit from immunotherapy. Patients with BRAF V600E variant melanoma brain metastases may thus benefit from alternative CNS-penetrant systemic regimens.

Multiomic Single Cell Sequencing Identifies Stemlike Nature of Mixed Phenotype Acute Leukemia and Provides Novel Risk Stratification.

Mixed phenotype acute leukemia (MPAL) is a leukemia whose biologic drivers are poorly understood, therapeutic strategy remains unclear, and prognosis is poor. We performed multiomic single cell (SC) profiling of 14 newly diagnosed adult MPAL patients to characterize the immunophenotypic, genetic, and transcriptional landscapes of MPAL. We show that neither genetic profile nor transcriptome reliably correlate with specific MPAL immunophenotypes. However, progressive acquisition of mutations is associated with increased expression of immunophenotypic markers of immaturity. Using SC transcriptional profiling, we find that MPAL blasts express a stem cell-like transcriptional profile distinct from other acute leukemias and indicative of high differentiation potential. Further, patients with the highest differentiation potential demonstrated inferior survival in our dataset. A gene set score, MPAL95, derived from genes highly enriched in this cohort, is applicable to bulk RNA sequencing data and was predictive of survival in an independent patient cohort, suggesting utility for clinical risk stratification.

Microfluidics-free single-cell genomics with templated emulsification.

Current single-cell RNA-sequencing approaches have limitations that stem from the microfluidic devices or fluid handling steps required for sample processing. We develop a method that does not require specialized microfluidic devices, expertise or hardware. Our approach is based on particle-templated emulsification, which allows single-cell encapsulation and barcoding of cDNA in uniform droplet emulsions with only a vortexer. Particle-templated instant partition sequencing (PIP-seq) accommodates a wide range of emulsification formats, including microwell plates and large-volume conical tubes, enabling thousands of samples or millions of cells to be processed in minutes. We demonstrate that PIP-seq produces high-purity transcriptomes in mouse-human mixing studies, is compatible with multiomics measurements and can accurately characterize cell types in human breast tissue compared to a commercial microfluidic platform. Single-cell transcriptional profiling of mixed phenotype acute leukemia using PIP-seq reveals the emergence of heterogeneity within chemotherapy-resistant cell subsets that were hidden by standard immunophenotyping. PIP-seq is a simple, flexible and scalable next-generation workflow that extends single-cell sequencing to new applications.

Characterizing cell interactions at scale with made-to-order droplet ensembles (MODEs).

Cell-cell interactions are important to numerous biological systems, including tissue microenvironments, the immune system, and cancer. However, current methods for studying cell combinations and interactions are limited in scalability, allowing just hundreds to thousands of multicell assays per experiment; this limited throughput makes it difficult to characterize interactions at biologically relevant scales. Here, we describe a paradigm in cell interaction profiling that allows accurate grouping of cells and characterization of their interactions for tens to hundreds of thousands of combinations. Our approach leverages high-throughput droplet microfluidics to construct multicellular combinations in a deterministic process that allows inclusion of programmed reagent mixtures and beads. The combination droplets are compatible with common manipulation and measurement techniques, including imaging, barcode-based genomics, and sorting. We demonstrate the approach by using it to enrich for chimeric antigen receptor (CAR)-T cells that activate upon incubation with target cells, a bottleneck in the therapeutic T cell engineering pipeline. The speed and control of our approach should enable valuable cell interaction studies.

Modular barcode beads for microfluidic single cell genomics.

Barcode beads allow efficient nucleic acid tagging in single cell genomics. Current barcode designs, however, are fabricated with a particular application in mind. Repurposing to novel targets, or altering to add additional targets as information is obtained is possible but the result is suboptimal. Here, we describe a modular framework that simplifies generation of multifunctional beads and allows their easy extension to new targets.

Joint profiling of DNA and proteins in single cells to dissect genotype-phenotype associations in leukemia.

Studies of acute myeloid leukemia rely on DNA sequencing and immunophenotyping by flow cytometry as primary tools for disease characterization. However, leukemia tumor heterogeneity complicates integration of DNA variants and immunophenotypes from separate measurements. Here we introduce DAb-seq, a technology for simultaneous capture of DNA genotype and cell surface phenotype from single cells at high throughput, enabling direct profiling of proteogenomic states in tens of thousands of cells. To demonstrate the approach, we analyze the disease of three patients with leukemia over multiple treatment timepoints and disease recurrences. We observe complex genotype-phenotype dynamics that illustrate the subtlety of the disease process and the degree of incongruity between blast cell genotype and phenotype in different clinical scenarios. Our results highlight the importance of combined single-cell DNA and protein measurements to fully characterize the heterogeneity of leukemia.

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.

Targeted Single-Cell RNA and DNA Sequencing With Fluorescence-Activated Droplet Merger.

Analyzing every cell in a diverse sample provides insight into population-level heterogeneity, but abundant cell types dominate the analysis and rarer populations are scarcely represented in the data. To focus on specific cell types, the current paradigm is to physically isolate subsets of interest prior to analysis; however, it remains difficult to isolate and then single-cell sequence such populations because of compounding losses. Here, we describe an alternative approach that selectively merges cells with reagents to achieve enzymatic reactions without having to physically isolate cells. We apply this technique to perform single-cell transcriptome and genome sequencing of specific cell subsets. Our method for analyzing heterogeneous populations obviates the need for pre- or post-enrichment and simplifies single-cell workflows, making it useful for other applications in single-cell biology, combinatorial chemical synthesis, and drug screening.

Activity of the mycobacterial proteasomal ATPase Mpa is reversibly regulated by pupylation.

Pupylation is a bacterial post-translational modification of target proteins on lysine residues with prokaryotic ubiquitin-like protein Pup. Pup-tagged substrates are recognized by a proteasome-interacting ATPase termed Mpa in Mycobacterium tuberculosis. Mpa unfolds pupylated substrates and threads them into the proteasome core particle for degradation. Interestingly, Mpa itself is also a pupylation target. Here, we show that the Pup ligase PafA predominantly produces monopupylated Mpa modified homogeneously on a single lysine residue within its C-terminal region. We demonstrate that this modification renders Mpa functionally inactive. Pupylated Mpa can no longer support Pup-mediated proteasomal degradation due to its inability to associate with the proteasome core. Mpa is further inactivated by rapid Pup- and ATPase-driven deoligomerization of the hexameric Mpa ring. We show that pupylation of Mpa is chemically and functionally reversible. Mpa regains its enzymatic activity upon depupylation by the depupylase Dop, affording a rapid and reversible activity control over Mpa function.