Somatic mutations and cell identity linked by Genotyping of Transcriptomes

Human embryonic genetic mosaicism and its effects on development and disease

Nearly every mammalian cell division is accompanied by a mutational event that becomes fixed in a daughter cell. When carried forward to additional cell progeny, a clone of variant cells can emerge. As a result, mammals are complex mosaics of clones that are genetically distinct from one another. Recent high-throughput sequencing studies have revealed that mosaicism is common, clone sizes often increase with age and specific variants can affect tissue function and disease development. Variants that are acquired during early embryogenesis are shared by multiple cell types and can affect numerous tissues. Within tissues, variant clones compete, which can result in their expansion or elimination. Embryonic mosaicism has clinical implications for genetic disease severity and transmission but is likely an under-recognized phenomenon. To better understand its implications for mosaic individuals, it is essential to leverage research tools that can elucidate the mechanisms by which expanded embryonic variants influence development and disease.

Megakaryocyte maturation involves activation of the adaptive unfolded protein response

Endoplasmic reticulum stress triggers the unfolded protein response (UPR) to promote cell survival or apoptosis. Transient endoplasmic reticulum stress activation has been reported to trigger megakaryocyte production, and UPR activation has been reported as a feature of megakaryocytic cancers. However, the role of UPR signaling in megakaryocyte biology is not fully understood. We studied the involvement of UPR in human megakaryocytic differentiation using PMA (phorbol 12-myristate 13-acetate)-induced maturation of megakaryoblastic cell lines and thrombopoietin-induced differentiation of human peripheral blood-derived progenitors. Our results demonstrate that an adaptive UPR is a feature of megakaryocytic differentiation and that this response is not associated with ER stress-induced apoptosis. Differentiation did not alter the response to the canonical endoplasmic reticulum stressors DTT or thapsigargin. However, thapsigargin, but not DTT, inhibited differentiation, consistent with the involvement of Ca2+ signaling in megakaryocyte differentiation.

The evolutionary theory of cancer: challenges and potential solutions

The clonal evolution model of cancer was developed in the 1950s-1970s and became central to cancer biology in the twenty-first century, largely through studies of cancer genetics. Although it has proven its worth, its structure has been challenged by observations of phenotypic plasticity, non-genetic forms of inheritance, non-genetic determinants of clone fitness and non-tree-like transmission of genes. There is even confusion about the definition of a clone, which we aim to resolve. The performance and value of the clonal evolution model depends on the empirical extent to which evolutionary processes are involved in cancer, and on its theoretical ability to account for those evolutionary processes. Here, we identify limits in the theoretical performance of the clonal evolution model and provide solutions to overcome those limits. Although we do not claim that clonal evolution can explain everything about cancer, we show how many of the complexities that have been identified in the dynamics of cancer can be integrated into the model to improve our current understanding of cancer.

Unveiling aging dynamics in the hematopoietic system insights from single-cell technologies

As the demographic structure shifts towards an aging society, strategies aimed at slowing down or reversing the aging process become increasingly essential. Aging is a major predisposing factor for many chronic diseases in humans. The hematopoietic system, comprising blood cells and their associated bone marrow microenvironment, intricately participates in hematopoiesis, coagulation, immune regulation and other physiological phenomena. The aging process triggers various alterations within the hematopoietic system, serving as a spectrum of risk factors for hematopoietic disorders, including clonal hematopoiesis, immune senescence, myeloproliferative neoplasms and leukemia. The emerging single-cell technologies provide novel insights into age-related changes in the hematopoietic system. In this review, we summarize recent studies dissecting hematopoietic system aging using single-cell technologies. We discuss cellular changes occurring during aging in the hematopoietic system at the levels of the genomics, transcriptomics, epigenomics, proteomics, metabolomics and spatial multi-omics. Finally, we contemplate the future prospects of single-cell technologies, emphasizing the impact they may bring to the field of hematopoietic system aging research.

Differential modulation of mutant CALR and JAK2 V617F-driven oncogenesis by HLA genotype in myeloproliferative neoplasms

It has been demonstrated previously that human leukocyte antigen class I (HLA-I) and class II (HLA-II) alleles may modulate JAK2 V617F and CALR mutation (CALRmut)-associated oncogenesis in myeloproliferative neoplasms (MPNs). However, the role of immunogenetic factors in MPNs remains underexplored. We aimed to investigate the potential involvement of HLA genes in CALRmut+ MPNs. High-resolution genotyping of HLA-I and -II loci was conducted in 42 CALRmut+ and 158 JAK2 V617F+ MPN patients and 1,083 healthy controls. A global analysis of the diversity of HLA-I genotypes revealed no significant differences between CALRmut+ patients and controls. However, one HLA-I allele (C*06:02) showed an inverse correlation with presence of CALR mutation. A meta-analysis across independent cohorts and healthy individuals from the 1000 Genomes Project confirmed an inverse correlation between the presentation capabilities of the HLA-I loci for JAK2 V617F and CALRmut-derived peptides in both patients and healthy individuals. scRNA-Seq analysis revealed low expression of TAP1 and CIITA genes in CALRmut+ hematopoietic stem and progenitor cells. In conclusion, the HLA-I genotype differentially restricts JAK2 V617F and CALRmut-driven oncogenesis potentially explaining the mutual exclusivity of the two mutations and differences in their presentation latency. These findings have practical implications for the development of neoantigen-based vaccines in MPNs.

Transcripts of repetitive DNA elements signal to block phagocytosis of hematopoietic stem cells

Macrophages maintain hematopoietic stem cell (HSC) quality by assessing cell surface Calreticulin (Calr), an "eat-me" signal induced by reactive oxygen species (ROS). Using zebrafish genetics, we identified Beta-2-microglobulin (B2m) as a crucial "don't eat-me" signal on blood stem cells. A chemical screen revealed inducers of surface Calr that promoted HSC proliferation without triggering ROS or macrophage clearance. Whole-genome CRISPR-Cas9 screening showed that Toll-like receptor 3 (Tlr3) signaling regulated b2m expression. Targeting b2m or tlr3 reduced the HSC clonality. Elevated B2m levels correlated with high expression of repetitive element (RE) transcripts. Overall, our data suggest that RE-associated double-stranded RNA could interact with TLR3 to stimulate surface expression of B2m on hematopoietic stem and progenitor cells. These findings suggest that the balance of Calr and B2m regulates macrophage-HSC interactions and defines hematopoietic clonality.

Cancer cell states: Lessons from ten years of single-cell RNA-sequencing of human tumors

Human tumors are intricate ecosystems composed of diverse genetic clones and malignant cell states that evolve in a complex tumor micro-environment. Single-cell RNA-sequencing (scRNA-seq) provides a compelling strategy to dissect this intricate biology and has enabled a revolution in our ability to understand tumor biology over the last ten years. Here we reflect on this first decade of scRNA-seq in human tumors and highlight some of the powerful insights gleaned from these studies. We first focus on computational approaches for robustly defining cancer cell states and their diversity and highlight some of the most common patterns of gene expression intra-tumor heterogeneity (eITH) observed across cancer types. We then discuss ambiguities in the field in defining and naming such eITH programs. Finally, we highlight critical developments that will facilitate future research and the broader implementation of these technologies in clinical settings.

New frameworks for hematopoiesis derived from single-cell genomics

ABSTRACT

Recent advancements in single-cell genomics have enriched our understanding of hematopoiesis, providing intricate details about hematopoietic stem cell biology, differentiation, and lineage commitment. Technological advancements have highlighted extensive heterogeneity of cell populations and continuity of differentiation routes. Nevertheless, intermediate "attractor" states signify structure in stem and progenitor populations that link state transition dynamics to fate potential. We discuss how innovative model systems quantify lineage bias and how stress accelerates differentiation, thereby reducing fate plasticity compared with native hematopoiesis. We conclude by offering our perspective on the current model of hematopoiesis and discuss how a more precise understanding can translate to strategies that extend healthy hematopoiesis and prevent disease.

Identifying cancer cells from calling single-nucleotide variants in scRNA-seq data

MOTIVATION

Single-cell RNA sequencing (scRNA-seq) data are widely used to study cancer cell states and their heterogeneity. However, the tumour microenvironment is usually a mixture of healthy and cancerous cells and it can be difficult to fully separate these two populations based on transcriptomics alone. If available, somatic single-nucleotide variants (SNVs) observed in the scRNA-seq data could be used to identify the cancer population and match that information with the single cells' expression profile. However, calling somatic SNVs in scRNA-seq data is a challenging task, as most variants seen in the short-read data are not somatic, but can instead be germline variants, RNA edits or transcription, sequencing, or processing errors. In addition, only variants present in actively transcribed regions for each individual cell will be seen in the data.

RESULTS

To address these challenges, we develop CCLONE (Cancer Cell Labelling On Noisy Expression), an interpretable tool adapted to handle the uncertainty and sparsity of SNVs called from scRNA-seq data. CCLONE jointly identifies cancer clonal populations, and their associated variants. We apply CCLONE on two acute myeloid leukaemia datasets and one lung adenocarcinoma dataset and show that CCLONE captures both genetic clones and somatic events for multiple patients. These results show how CCLONE can be used to gather insight into the course of the disease and the origin of cancer cells in scRNA-seq data.

AVAILABILITY AND IMPLEMENTATION

Source code is available at github.com/HaghverdiLab/CCLONE.

The stress response regulator HSF1 modulates natural killer cell anti-tumour immunity

Diverse cellular insults converge on activation of the heat shock factor 1 (HSF1), which regulates the proteotoxic stress response to maintain protein homoeostasis. HSF1 regulates numerous gene programmes beyond the proteotoxic stress response in a cell-type- and context-specific manner to promote malignancy. However, the role(s) of HSF1 in immune populations of the tumour microenvironment remain elusive. Here, we leverage an in vivo model of HSF1 activation and single-cell transcriptomic tumour profiling to show that augmented HSF1 activity in natural killer (NK) cells impairs cytotoxicity, cytokine production and subsequent anti-tumour immunity. Mechanistically, HSF1 directly binds and regulates the expression of key mediators of NK cell effector function. This work demonstrates that HSF1 regulates the immune response under the stress conditions of the tumour microenvironment. These findings have important implications for enhancing the efficacy of adoptive NK cell therapies and for designing combinatorial strategies including modulators of NK cell-mediated tumour killing.

Unlocking the Complexity: Exploration of Acute Lymphoblastic Leukemia at the Single Cell Level

Acute lymphoblastic leukemia (ALL) is the most common cancer in children. ALL originates from precursor lymphocytes that acquire multiple genomic changes over time, including chromosomal rearrangements and point mutations. While a large variety of genomic defects was identified and characterized in ALL over the past 30 years, it was only in recent years that the clonal heterogeneity was recognized. Thanks to the latest advancements in single-cell sequencing techniques, which have evolved from the analysis of a few hundred cells to the analysis of thousands of cells simultaneously, the study of tumor heterogeneity now becomes possible. Different modalities can be explored at the single-cell level: DNA, RNA, epigenetic modifications, and intracellular and cell surface proteins. In this review, we describe these techniques and highlight their advantages and limitations in the study of ALL biology. Moreover, multiomics technologies and the incorporation of the spatial dimension can provide insight into intercellular communication. We describe how the different single-cell sequencing technologies help to unravel the molecular complexity of ALL, shedding light on its development, its heterogeneity, its interaction with the leukemia microenvironment and possible relapse mechanisms.

Transient Differentiation-State Plasticity Occurs during Acute Lymphoblastic Leukemia Initiation

Leukemia is characterized by oncogenic lesions that result in a block of differentiation, whereas phenotypic plasticity is retained. A better understanding of how these two phenomena arise during leukemogenesis in humans could help inform diagnosis and treatment strategies. Here, we leveraged the well-defined differentiation states during T-cell development to pinpoint the initiation of T-cell acute lymphoblastic leukemia (T-ALL), an aggressive form of childhood leukemia, and study the emergence of phenotypic plasticity. Single-cell whole genome sequencing of leukemic blasts was combined with multiparameter flow cytometry to couple cell identity and clonal lineages. Irrespective of genetic events, leukemia-initiating cells altered their phenotypes by differentiation and dedifferentiation. The construction of the phylogenies of individual leukemias using somatic mutations revealed that phenotypic diversity is reflected by the clonal structure of cancer. The analysis also indicated that the acquired phenotypes are heritable and stable. Together, these results demonstrate a transient period of plasticity during leukemia initiation, where phenotypic switches seem unidirectional. Significance: A method merging multicolor flow cytometry with single-cell whole genome sequencing to couple cell identity with clonal lineages uncovers differentiation-state plasticity in leukemia, reconciling blocked differentiation with phenotypic plasticity in cancer.

Data acquisition approaches for single cell proteomics

Single-cell proteomics (SCP) aims to characterize the proteome of individual cells, providing insights into complex biological systems. It reveals subtle differences in distinct cellular populations that bulk proteome analysis may overlook, which is essential for understanding disease mechanisms and developing targeted therapies. Mass spectrometry (MS) methods in SCP allow the identification and quantification of thousands of proteins from individual cells. Two major challenges in SCP are the limited material in single-cell samples necessitating highly sensitive analytical techniques and the efficient processing of samples, as each biological sample requires thousands of single cell measurements. This review discusses MS advancements to mitigate these challenges using data-dependent acquisition (DDA) and data-independent acquisition (DIA). Additionally, we examine the use of short liquid chromatography gradients and sample multiplexing methods that increase the sample throughput and scalability of SCP experiments. We believe these methods will pave the way for improving our understanding of cellular heterogeneity and its implications for systems biology.

Multiomic profiling of human clonal hematopoiesis reveals genotype and cell-specific inflammatory pathway activation

ABSTRACT

Clonal hematopoiesis (CH) is an age-associated phenomenon that increases the risk of hematologic malignancy and cardiovascular disease. CH is thought to enhance disease risk through inflammation in the peripheral blood.1 Here, we profile peripheral blood gene expression in 66 968 single cells from a cohort of 17 patients with CH and 7 controls. Using a novel mitochondrial DNA barcoding approach, we were able to identify and separately compare mutant Tet methylcytosine dioxygenase 2 (TET2) and DNA methyltransferase 3A (DNMT3A) cells with nonmutant counterparts. We discovered the vast majority of mutated cells were in the myeloid compartment. Additionally, patients harboring DNMT3A and TET2 CH mutations possessed a proinflammatory profile in CD14+ monocytes through previously unrecognized pathways such as galectin and macrophage inhibitory factor. We also found that T cells from patients with CH, although mostly unmutated, had decreased expression of GTPase of the immunity associated protein genes, which are critical to T-cell development, suggesting that CH impairs T-cell function.

Cell Type- and Age-Specific Expression of lncRNAs across Kidney Cell Types

Key Points

Background

Accumulated evidence demonstrates that long noncoding RNAs (lncRNAs) regulate cell differentiation and homeostasis, influencing kidney aging and disease. Despite their versatility, the function of lncRNA remains poorly understood because of the lack of a reference map of lncRNA transcriptome in various cell types.

Methods

In this study, we used a targeted single-cell RNA sequencing method to enrich and characterize lncRNAs in individual cells. We applied this method to various mouse tissues, including normal and aged kidneys.

Results

Through tissue-specific clustering analysis, we identified cell type–specific lncRNAs that showed a high correlation with known cell-type marker genes. Furthermore, we constructed gene regulatory networks to explore the functional roles of differentially expressed lncRNAs in each cell type. In the kidney, we observed dynamic expression changes of lncRNAs during aging, with specific changes in glomerular cells. These cell type– and age-specific expression patterns of lncRNAs suggest that lncRNAs may have a potential role in regulating cellular processes, such as immune response and energy metabolism, during kidney aging.

Conclusions

Our study sheds light on the comprehensive landscape of lncRNA expression and function and provides a valuable resource for future analysis of lncRNAs (https://gist-fgl.github.io/sc-lncrna-atlas/).

Discovering mechanisms of human genetic variation and controlling cell states at scale

Population-scale sequencing efforts have catalogued substantial genetic variation in humans such that variant discovery dramatically outpaces interpretation. We discuss how single-cell sequencing is poised to reveal genetic mechanisms at a rate that may soon approach that of variant discovery. The functional genomics toolkit is sufficiently modular to systematically profile almost any type of variation within increasingly diverse contexts and with molecularly comprehensive and unbiased readouts. As a result, we can construct deep phenotypic atlases of variant effects that span the entire regulatory cascade. The same conceptual approach to interpreting genetic variation should be applied to engineering therapeutic cell states. In this way, variant mechanism discovery and cell state engineering will become reciprocating and iterative processes towards genomic medicine.

GoT-Splice protocol for multi-omics profiling of gene expression, cell-surface proteins, mutational status, and RNA splicing in human cells

Studying RNA splicing factor mutations is challenging due to difficulties in distinguishing wild-type and mutant cells within complex human tissues and inaccuracies associated with reconstructing splicing signals from short-read sequencing data. Here, we present Genotyping of Transcriptomes (GoT)-Splice, a protocol that overcomes these limitations by combining GoT with enhanced long-read single-cell transcriptome and cell-surface proteomics profiling. We describe steps for long-read library preparation and analysis, followed by cDNA re-amplification, enrichment of mutation of interest, sample indexing, and GoT library preparation. For complete details on the use and execution of this protocol, please refer to Cortés-López et al.1.

Massively parallel in vivo Perturb-seq reveals cell-type-specific transcriptional networks in cortical development

Leveraging AAVs' versatile tropism and labeling capacity, we expanded the scale of in vivo CRISPR screening with single-cell transcriptomic phenotyping across embryonic to adult brains and peripheral nervous systems. Through extensive tests of 86 vectors across AAV serotypes combined with a transposon system, we substantially amplified labeling efficacy and accelerated in vivo gene delivery from weeks to days. Our proof-of-principle in utero screen identified the pleiotropic effects of Foxg1, highlighting its tight regulation of distinct networks essential for cell fate specification of Layer 6 corticothalamic neurons. Notably, our platform can label >6% of cerebral cells, surpassing the current state-of-the-art efficacy at <0.1% by lentivirus, to achieve analysis of over 30,000 cells in one experiment and enable massively parallel in vivo Perturb-seq. Compatible with various phenotypic measurements (single-cell or spatial multi-omics), it presents a flexible approach to interrogate gene function across cell types in vivo, translating gene variants to their causal function.

Single-cell genomics-based immune and disease monitoring in blood malignancies

Achieving long-term disease control using therapeutic immunomodulation is a long-standing concept with a strong tradition in blood malignancies. Besides allogeneic hematopoietic stem cell transplantation that continues to provide potentially curative treatment for otherwise challenging diagnoses, recent years have seen impressive progress in immunotherapies for leukemias and lymphomas with immune checkpoint blockade, bispecific monoclonal antibodies, and CAR T cell therapies. Despite their success, non-response, relapse, and immune toxicities remain frequent, thus prioritizing the elucidation of the underlying mechanisms and identifying predictive biomarkers. The increasing availability of single-cell genomic tools now provides a system's immunology view to resolve the molecular and cellular mechanisms of immunotherapies at unprecedented resolution. Here, we review recent studies that leverage these technological advancements for tracking immune responses, the emergence of immune resistance, and toxicities. As single-cell immune monitoring tools evolve and become more accessible, we expect their wide adoption for routine clinical applications to catalyze more precise therapeutic steering of personal immune responses.

Mutation and cell state compatibility is required and targetable in Ph+ acute lymphoblastic leukemia minimal residual disease

Efforts to cure BCR::ABL1 B cell acute lymphoblastic leukemia (Ph+ ALL) solely through inhibition of ABL1 kinase activity have thus far been insufficient despite the availability of tyrosine kinase inhibitors (TKIs) with broad activity against resistance mutants. The mechanisms that drive persistence within minimal residual disease (MRD) remain poorly understood and therefore untargeted. Utilizing 13 patient-derived xenograft (PDX) models and clinical trial specimens of Ph+ ALL, we examined how genetic and transcriptional features co-evolve to drive progression during prolonged TKI response. Our work reveals a landscape of cooperative mutational and transcriptional escape mechanisms that differ from those causing resistance to first generation TKIs. By analyzing MRD during remission, we show that the same resistance mutation can either increase or decrease cellular fitness depending on transcriptional state. We further demonstrate that directly targeting transcriptional state-associated vulnerabilities at MRD can overcome BCR::ABL1 independence, suggesting a new paradigm for rationally eradicating MRD prior to relapse. Finally, we illustrate how cell mass measurements of leukemia cells can be used to rapidly monitor dominant transcriptional features of Ph+ ALL to help rationally guide therapeutic selection from low-input samples.

Application of omics in the diagnosis, prognosis, and treatment of acute myeloid leukemia

Acute myeloid leukemia (AML) is the most frequent leukemia in adults with a high mortality rate. Current diagnostic criteria and selections of therapeutic strategies are generally based on gene mutations and cytogenetic abnormalities. Chemotherapy, targeted therapies, and hematopoietic stem cell transplantation (HSCT) are the major therapeutic strategies for AML. Two dilemmas in the clinical management of AML are related to its poor prognosis. One is the inaccurate risk stratification at diagnosis, leading to incorrect treatment selections. The other is the frequent resistance to chemotherapy and/or targeted therapies. Genomic features have been the focus of AML studies. However, the DNA-level aberrations do not always predict the expression levels of genes and proteins and the latter is more closely linked to disease phenotypes. With the development of high-throughput sequencing and mass spectrometry technologies, studying downstream effectors including RNA, proteins, and metabolites becomes possible. Transcriptomics can reveal gene expression and regulatory networks, proteomics can discover protein expression and signaling pathways intimately associated with the disease, and metabolomics can reflect precise changes in metabolites during disease progression. Moreover, omics profiling at the single-cell level enables studying cellular components and hierarchies of the AML microenvironment. The abundance of data from different omics layers enables the better risk stratification of AML by identifying prognosis-related biomarkers, and has the prospective application in identifying drug targets, therefore potentially discovering solutions to the two dilemmas. In this review, we summarize the existing AML studies using omics methods, both separately and combined, covering research fields of disease diagnosis, risk stratification, prognosis prediction, chemotherapy, as well as targeted therapy. Finally, we discuss the directions and challenges in the application of multi-omics in precision medicine of AML. Our review may inspire both omics researchers and clinical physicians to study AML from a different angle.

Genotype-to-phenotype mapping of somatic clonal mosaicism via single-cell co-capture of DNA mutations and mRNA transcripts

Somatic mosaicism is a hallmark of malignancy that is also pervasively observed in human physiological aging, with clonal expansions of cells harboring mutations in recurrently mutated driver genes. Bulk sequencing of tissue microdissection captures mutation frequencies, but cannot distinguish which mutations co-occur in the same clones to reconstruct clonal architectures, nor phenotypically profile clonal populations to delineate how driver mutations impact cellular behavior. To address these challenges, we developed single-cell Genotype-to-Phenotype sequencing (scG2P) for high-throughput, highly-multiplexed, single-cell joint capture of recurrently mutated genomic regions and mRNA phenotypic markers in cells or nuclei isolated from solid tissues. We applied scG2P to aged esophagus samples from five individuals with high alcohol and tobacco exposure and observed a clonal landscape dominated by a large number of clones with a single driver event, but only rare clones with two driver mutations. NOTCH1 mutants dominate the clonal landscape and are linked to stunted epithelial differentiation, while TP53 mutants and double-driver mutants promote clonal expansion through both differentiation biases and increased cell cycling. Thus, joint single-cell highly multiplexed capture of somatic mutations and mRNA transcripts enables high resolution reconstruction of clonal architecture and associated phenotypes in solid tissue somatic mosaicism.

Cell of origin epigenetic priming determines susceptibility to Tet2 mutation

Hematopoietic stem cell (HSC) mutations can result in clonal hematopoiesis (CH) with heterogeneous clinical outcomes. Here, we investigate how the cell state preceding Tet2 mutation impacts the pre-malignant phenotype. Using an inducible system for clonal analysis of myeloid progenitors, we find that the epigenetic features of clones at similar differentiation status are highly heterogeneous and functionally respond differently to Tet2 mutation. Cell differentiation stage also influences Tet2 mutation response indicating that the cell of origin's epigenome modulates clone-specific behaviors in CH. Molecular features associated with higher risk outcomes include Sox4 that sensitizes cells to Tet2 inactivation, inducing dedifferentiation, altered metabolism and increasing the in vivo clonal output of mutant cells, as confirmed in primary GMP and HSC models. Our findings validate the hypothesis that epigenetic features can predispose specific clones for dominance, explaining why identical genetic mutations can result in different phenotypes.

Analyzing the functional effects of DNA variants with gene editing

Continual advancements in genomics have led to an ever-widening disparity between the rate of discovery of genetic variants and our current understanding of their functions and potential roles in disease. Systematic methods for phenotyping DNA variants are required to effectively translate genomics data into improved outcomes for patients with genetic diseases. To make the biggest impact, these approaches must be scalable and accurate, faithfully reflect disease biology, and define complex disease mechanisms. We compare current methods to analyze the function of variants in their endogenous DNA context using genome editing strategies, such as saturation genome editing, base editing and prime editing. We discuss how these technologies can be linked to high-content readouts to gain deep mechanistic insights into variant effects. Finally, we highlight key challenges that need to be addressed to bridge the genotype to phenotype gap, and ultimately improve the diagnosis and treatment of genetic diseases.

De novo detection of somatic mutations in high-throughput single-cell profiling data sets

Characterization of somatic mutations at single-cell resolution is essential to study cancer evolution, clonal mosaicism and cell plasticity. Here, we describe SComatic, an algorithm designed for the detection of somatic mutations in single-cell transcriptomic and ATAC-seq (assay for transposase-accessible chromatin sequence) data sets directly without requiring matched bulk or single-cell DNA sequencing data. SComatic distinguishes somatic mutations from polymorphisms, RNA-editing events and artefacts using filters and statistical tests parameterized on non-neoplastic samples. Using >2.6 million single cells from 688 single-cell RNA-seq (scRNA-seq) and single-cell ATAC-seq (scATAC-seq) data sets spanning cancer and non-neoplastic samples, we show that SComatic detects mutations in single cells accurately, even in differentiated cells from polyclonal tissues that are not amenable to mutation detection using existing methods. Validated against matched genome sequencing and scRNA-seq data, SComatic achieves F1 scores between 0.6 and 0.7 across diverse data sets, in comparison to 0.2-0.4 for the second-best performing method. In summary, SComatic permits de novo mutational signature analysis, and the study of clonal heterogeneity and mutational burdens at single-cell resolution.

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.

The molecular map of CLL and Richter's syndrome

Clonal expansion of B-cells, from the early stages of monoclonal B-cell lymphocytosis through to chronic lymphocytic leukemia (CLL), and then in some cases to Richter's syndrome (RS) provides a comprehensive model of cancer evolution, notable for the marked morphological transformation and distinct clinical phenotypes. High-throughput sequencing of large cohorts of patients and single-cell studies have generated a molecular map of CLL and more recently, of RS, yielding fundamental insights into these diseases and of clonal evolution. A selection of CLL driver genes have been functionally interrogated to yield novel insights into the biology of CLL. Such findings have the potential to impact patient care through risk stratification, treatment selection and drug discovery. However, this molecular map remains incomplete, with extant questions concerning the origin of the B-cell clone, the role of the TME, inter- and intra-compartmental heterogeneity and of therapeutic resistance mechanisms. Through the application of multi-modal single-cell technologies across tissues, disease states and clinical contexts, these questions can now be addressed with the answers holding great promise of generating translatable knowledge to improve patient care.

Deciphering cell states and genealogies of human haematopoiesis

The human blood system is maintained through the differentiation and massive amplification of a limited number of long-lived haematopoietic stem cells (HSCs)1. Perturbations to this process underlie diverse diseases, but the clonal contributions to human haematopoiesis and how this changes with age remain incompletely understood. Although recent insights have emerged from barcoding studies in model systems2-5, simultaneous detection of cell states and phylogenies from natural barcodes in humans remains challenging. Here we introduce an improved, single-cell lineage-tracing system based on deep detection of naturally occurring mitochondrial DNA mutations with simultaneous readout of transcriptional states and chromatin accessibility. We use this system to define the clonal architecture of HSCs and map the physiological state and output of clones. We uncover functional heterogeneity in HSC clones, which is stable over months and manifests as both differences in total HSC output and biases towards the production of different mature cell types. We also find that the diversity of HSC clones decreases markedly with age, leading to an oligoclonal structure with multiple distinct clonal expansions. Our study thus provides a clonally resolved and cell-state-aware atlas of human haematopoiesis at single-cell resolution, showing an unappreciated functional diversity of human HSC clones and, more broadly, paving the way for refined studies of clonal dynamics across a range of tissues in human health and disease.

Computational and Systems Biology Advances to Enable Bioagent Agnostic Signatures

Enumerated threat agent lists have long driven biodefense priorities. The global SARS-CoV-2 pandemic demonstrated the limitations of searching for known threat agents as compared to a more agnostic approach. Recent technological advances are enabling agent-agnostic biodefense, especially through the integration of multi-modal observations of host-pathogen interactions directed by a human immunological model. Although well-developed technical assays exist for many aspects of human-pathogen interaction, the analytic methods and pipelines to combine and holistically interpret the results of such assays are immature and require further investments to exploit new technologies. In this manuscript, we discuss potential immunologically based bioagent-agnostic approaches and the computational tool gaps the community should prioritize filling.

Development and future of droplet microfluidics

Over the past two decades, advances in droplet-based microfluidics have facilitated new approaches to process and analyze samples with unprecedented levels of precision and throughput. A wide variety of applications has been inspired across multiple disciplines ranging from materials science to biology. Understanding the dynamics of droplets enables optimization of microfluidic operations and design of new techniques tailored to emerging demands. In this review, we discuss the underlying physics behind high-throughput generation and manipulation of droplets. We also summarize the applications in droplet-derived materials and droplet-based lab-on-a-chip biotechnology. In addition, we offer perspectives on future directions to realize wider use of droplet microfluidics in industrial production and biomedical analyses.

Molecular Genetic Profile of Myelofibrosis: Implications in the Diagnosis, Prognosis, and Treatment Advancements

Myelofibrosis (MF) is an essential element of primary myelofibrosis, whereas secondary MF may develop in the advanced stages of other myeloid neoplasms, especially polycythemia vera and essential thrombocythemia. Over the last two decades, advances in molecular diagnostic techniques, particularly the integration of next-generation sequencing in clinical laboratories, have revolutionized the diagnosis, classification, and clinical decision making of myelofibrosis. Driver mutations involving JAK2, CALR, and MPL induce hyperactivity in the JAK-STAT signaling pathway, which plays a central role in cell survival and proliferation. Approximately 80% of myelofibrosis cases harbor additional mutations, frequently in the genes responsible for epigenetic regulation and RNA splicing. Detecting these mutations is crucial for diagnosing myeloproliferative neoplasms (MPNs), especially in cases where no mutations are present in the three driver genes (triple-negative MPNs). While fibrosis in the bone marrow results from the disturbance of inflammatory cytokines, it is fundamentally associated with mutation-driven hematopoiesis. The mutation profile and order of acquiring diverse mutations influence the MPN phenotype. Mutation profiling reveals clonal diversity in MF, offering insights into the clonal evolution of neoplastic progression. Prognostic prediction plays a pivotal role in guiding the treatment of myelofibrosis. Mutation profiles and cytogenetic abnormalities have been integrated into advanced prognostic scoring systems and personalized risk stratification for MF. Presently, JAK inhibitors are part of the standard of care for MF, with newer generations developed for enhanced efficacy and reduced adverse effects. However, only a minority of patients have achieved a significant molecular-level response. Clinical trials exploring innovative approaches, such as combining hypomethylation agents that target epigenetic regulators, drugs proven effective in myelodysplastic syndrome, or immune and inflammatory modulators with JAK inhibitors, have demonstrated promising results. These combinations may be more effective in patients with high-risk mutations and complex mutation profiles. Expanding mutation profiling studies with more sensitive and specific molecular methods, as well as sequencing a broader spectrum of genes in clinical patients, may reveal molecular mechanisms in cases currently lacking detectable driver mutations, provide a better understanding of the association between genetic alterations and clinical phenotypes, and offer valuable information to advance personalized treatment protocols to improve long-term survival and eradicate mutant clones with the hope of curing MF.

scSNV-seq: high-throughput phenotyping of single nucleotide variants by coupled single-cell genotyping and transcriptomics

CRISPR screens with single-cell transcriptomic readouts are a valuable tool to understand the effect of genetic perturbations including single nucleotide variants (SNVs) associated with diseases. Interpretation of these data is currently limited as genotypes cannot be accurately inferred from guide RNA identity alone. scSNV-seq overcomes this limitation by coupling single-cell genotyping and transcriptomics of the same cells enabling accurate and high-throughput screening of SNVs. Analysis of variants across the JAK1 gene with scSNV-seq demonstrates the importance of determining the precise genetic perturbation and accurately classifies clinically observed missense variants into three functional categories: benign, loss of function, and separation of function.

Cancer Evolution: A Multifaceted Affair

Cancer cells adapt and survive through the acquisition and selection of molecular modifications. This process defines cancer evolution. Building on a theoretical framework based on heritable genetic changes has provided insights into the mechanisms supporting cancer evolution. However, cancer hallmarks also emerge via heritable nongenetic mechanisms, including epigenetic and chromatin topological changes, and interactions between tumor cells and the tumor microenvironment. Recent findings on tumor evolutionary mechanisms draw a multifaceted picture where heterogeneous forces interact and influence each other while shaping tumor progression. A comprehensive characterization of the cancer evolutionary toolkit is required to improve personalized medicine and biomarker discovery.

SIGNIFICANCE

Tumor evolution is fueled by multiple enabling mechanisms. Importantly, genetic instability, epigenetic reprogramming, and interactions with the tumor microenvironment are neither alternative nor independent evolutionary mechanisms. As demonstrated by findings highlighted in this perspective, experimental and theoretical approaches must account for multiple evolutionary mechanisms and their interactions to ultimately understand, predict, and steer tumor evolution.

Integrative genotyping of cancer and immune phenotypes by long-read sequencing

Single-cell transcriptomics has become the definitive method for classifying cell types and states, and can be augmented with genotype information to improve cell lineage identification. Due to constraints of short-read sequencing, current methods to detect natural genetic barcodes often require cumbersome primer panels and early commitment to targets. Here we devise a flexible long-read sequencing workflow and analysis pipeline, termed nanoranger, that starts from intermediate single-cell cDNA libraries to detect cell lineage-defining features, including single-nucleotide variants, fusion genes, isoforms, sequences of chimeric antigen and TCRs. Through systematic analysis of these classes of natural 'barcodes', we define the optimal targets for nanoranger, namely those loci close to the 5' end of highly expressed genes with transcript lengths shorter than 4 kB. As proof-of-concept, we apply nanoranger to longitudinal tracking of subclones of acute myeloid leukemia (AML) and describe the heterogeneous isoform landscape of thousands of marrow-infiltrating immune cells. We propose that enhanced cellular genotyping using nanoranger can improve the tracking of single-cell tumor and immune cell co-evolution.

De novo identification of expressed cancer somatic mutations from single-cell RNA sequencing data

Identifying expressed somatic mutations from single-cell RNA sequencing data de novo is challenging but highly valuable. We propose RESA - Recurrently Expressed SNV Analysis, a computational framework to identify expressed somatic mutations from scRNA-seq data. RESA achieves an average precision of 0.77 on three in silico spike-in datasets. In extensive benchmarking against existing methods using 19 datasets, RESA consistently outperforms them. Furthermore, we applied RESA to analyze intratumor mutational heterogeneity in a melanoma drug resistance dataset. By enabling high precision detection of expressed somatic mutations, RESA substantially enhances the reliability of mutational analysis in scRNA-seq. RESA is available at https://github.com/ShenLab-Genomics/RESA .

A protocol for simultaneous high-sensitivity genotyping and chromatin accessibility profiling in single cells

Single-cell assay for transposase-accessible chromatin with sequencing (scATAC-seq) resolves the heterogeneity of epigenetic states across cells but does not typically capture exonic mutations, which limits our knowledge of how somatic mutations alter chromatin landscapes. Here, we present a plate-based approach coupling high-sensitivity genotyping of genomic loci with high-content scATAC-seq libraries from the same single cells. We first describe steps for optimization of genotyping primers, followed by detailed guidance on the preparation of both scATAC-seq and single-cell genotyping libraries, fully automated on high-throughput liquid handling platforms. For complete details on the use and execution of this protocol, please refer to Turkalj, Jakobsen et al.1.

Preleukemic single-cell landscapes reveal mutation-specific mechanisms and gene programs predictive of AML patient outcomes

Acute myeloid leukemia (AML) and myeloid neoplasms develop through acquisition of somatic mutations that confer mutation-specific fitness advantages to hematopoietic stem and progenitor cells. However, our understanding of mutational effects remains limited to the resolution attainable within immunophenotypically and clinically accessible bulk cell populations. To decipher heterogeneous cellular fitness to preleukemic mutational perturbations, we performed single-cell RNA sequencing of eight different mouse models with driver mutations of myeloid malignancies, generating 269,048 single-cell profiles. Our analysis infers mutation-driven perturbations in cell abundance, cellular lineage fate, cellular metabolism, and gene expression at the continuous resolution, pinpointing cell populations with transcriptional alterations associated with differentiation bias. We further develop an 11-gene scoring system (Stem11) on the basis of preleukemic transcriptional signatures that predicts AML patient outcomes. Our results demonstrate that a single-cell-resolution deep characterization of preleukemic biology has the potential to enhance our understanding of AML heterogeneity and inform more effective risk stratification strategies.

Myeloproliferative Neoplasms: Contemporary Review and Molecular Landscape

Myelofibrosis (MF), Myeloproliferative neoplasms (MPNs), and MDS/MPN overlap syndromes have a broad range of clinical presentations and molecular abnormalities, making their diagnosis and classification complex. This paper reviews molecular aberration, epigenetic modifications, chromosomal anomalies, and their interactions with cellular and other immune mechanisms in the manifestations of these disease spectra, clinical features, classification, and treatment modalities. The advent of new-generation sequencing has broadened the understanding of the genetic factors involved. However, while great strides have been made in the pharmacological treatment of these diseases, treatment of advanced disease remains hematopoietic stem cell transplant.

A Bayesian method to infer copy number clones from single-cell RNA and ATAC sequencing

Single-cell RNA and ATAC sequencing technologies enable the examination of gene expression and chromatin accessibility in individual cells, providing insights into cellular phenotypes. In cancer research, it is important to consistently analyze these states within an evolutionary context on genetic clones. Here we present CONGAS+, a Bayesian model to map single-cell RNA and ATAC profiles onto the latent space of copy number clones. CONGAS+ clusters cells into tumour subclones with similar ploidy, rendering straightforward to compare their expression and chromatin profiles. The framework, implemented on GPU and tested on real and simulated data, scales to analyse seamlessly thousands of cells, demonstrating better performance than single-molecule models, and supporting new multi-omics assays. In prostate cancer, lymphoma and basal cell carcinoma, CONGAS+ successfully identifies complex subclonal architectures while providing a coherent mapping between ATAC and RNA, facilitating the study of genotype-phenotype maps and their connection to genomic instability.

A targeted sequencing extension for transcript genotyping in single-cell transcriptomics

As no existing methods within the single-cell RNA sequencing repertoire combine genotyping of specific genomic loci with high throughput, we evaluated a straightforward, targeted sequencing approach as an extension to high-throughput droplet-based single-cell RNA sequencing. Overlaying standard gene expression data with transcript level genotype information provides a strategy to study the impact of genetic variants. Here, we describe this targeted sequencing extension, explain how to process the data and evaluate how technical parameters such as amount of input cDNA, number of amplification rounds, and sequencing depth influence the number of transcripts detected. Finally, we demonstrate how targeted sequencing can be used in two contexts: (1) simultaneous investigation of the presence of a somatic variant and its potential impact on the transcriptome of affected cells and (2) evaluation of allele-specific expression of a germline variant in ad hoc cell subsets. Through these and other comparable applications, our targeted sequencing extension has the potential to improve our understanding of functional effects caused by genetic variation.

Proteogenetic drug response profiling elucidates targetable vulnerabilities of myelofibrosis

Myelofibrosis is a hematopoietic stem cell disorder belonging to the myeloproliferative neoplasms. Myelofibrosis patients frequently carry driver mutations in either JAK2 or Calreticulin (CALR) and have limited therapeutic options. Here, we integrate ex vivo drug response and proteotype analyses across myelofibrosis patient cohorts to discover targetable vulnerabilities and associated therapeutic strategies. Drug sensitivities of mutated and progenitor cells were measured in patient blood using high-content imaging and single-cell deep learning-based analyses. Integration with matched molecular profiling revealed three targetable vulnerabilities. First, CALR mutations drive BET and HDAC inhibitor sensitivity, particularly in the absence of high Ras pathway protein levels. Second, an MCM complex-high proliferative signature corresponds to advanced disease and sensitivity to drugs targeting pro-survival signaling and DNA replication. Third, homozygous CALR mutations result in high endoplasmic reticulum (ER) stress, responding to ER stressors and unfolded protein response inhibition. Overall, our integrated analyses provide a molecularly motivated roadmap for individualized myelofibrosis patient treatment.

Massively parallel in vivo Perturb-seq reveals cell type-specific transcriptional networks in cortical development

Systematic analysis of gene function across diverse cell types in vivo is hindered by two challenges: obtaining sufficient cells from live tissues and accurately identifying each cell's perturbation in high-throughput single-cell assays. Leveraging AAV's versatile cell type tropism and high labeling capacity, we expanded the resolution and scale of in vivo CRISPR screens: allowing phenotypic analysis at single-cell resolution across a multitude of cell types in the embryonic brain, adult brain, and peripheral nervous system. We undertook extensive tests of 86 AAV serotypes, combined with a transposon system, to substantially amplify labeling and accelerate in vivo gene delivery from weeks to days. Using this platform, we performed an in utero genetic screen as proof-of-principle and identified pleiotropic regulatory networks of Foxg1 in cortical development, including Layer 6 corticothalamic neurons where it tightly controls distinct networks essential for cell fate specification. Notably, our platform can label >6% of cerebral cells, surpassing the current state-of-the-art efficacy at <0.1% (mediated by lentivirus), and achieve analysis of over 30,000 cells in one experiment, thus enabling massively parallel in vivo Perturb-seq. Compatible with various perturbation techniques (CRISPRa/i) and phenotypic measurements (single-cell or spatial multi-omics), our platform presents a flexible, modular approach to interrogate gene function across diverse cell types in vivo, connecting gene variants to their causal functions.

Single-cell multi-omics defines the cell-type-specific impact of splicing aberrations in human hematopoietic clonal outgrowths

RNA splicing factors are recurrently mutated in clonal blood disorders, but the impact of dysregulated splicing in hematopoiesis remains unclear. To overcome technical limitations, we integrated genotyping of transcriptomes (GoT) with long-read single-cell transcriptomics and proteogenomics for single-cell profiling of transcriptomes, surface proteins, somatic mutations, and RNA splicing (GoT-Splice). We applied GoT-Splice to hematopoietic progenitors from myelodysplastic syndrome (MDS) patients with mutations in the core splicing factor SF3B1. SF3B1mut cells were enriched in the megakaryocytic-erythroid lineage, with expansion of SF3B1mut erythroid progenitor cells. We uncovered distinct cryptic 3' splice site usage in different progenitor populations and stage-specific aberrant splicing during erythroid differentiation. Profiling SF3B1-mutated clonal hematopoiesis samples revealed that erythroid bias and cell-type-specific cryptic 3' splice site usage in SF3B1mut cells precede overt MDS. Collectively, GoT-Splice defines the cell-type-specific impact of somatic mutations on RNA splicing, from early clonal outgrowths to overt neoplasia, directly in human samples.

Clonal evolution and hierarchy in myeloid malignancies

Myeloid malignancies, a group of hematopoietic disorders that includes acute myeloid leukemia (AML), myelodysplastic syndromes (MDS), and myeloproliferative neoplasms (MPNs), are caused by the accumulation of genetic and epigenetic changes in hematopoietic stem and progenitor cells (HSPCs) over time. Despite the relatively low number of genomic drivers compared with other forms of cancer, the process by which these changes shape the genomic architecture of myeloid malignancies remains elusive. Recent advancements in clonal hematopoiesis research and the use of cutting-edge single cell technologies have shed new light on the developmental process of myeloid malignancies. In this review, we delve into the intricacies of clonal evolution in myeloid malignancies and its implications for the development of new diagnostic and therapeutic approaches.

Loss of SYNCRIP unleashes APOBEC-driven mutagenesis, tumor heterogeneity, and AR-targeted therapy resistance in prostate cancer

Tumor mutational burden and heterogeneity has been suggested to fuel resistance to many targeted therapies. The cytosine deaminase APOBEC proteins have been implicated in the mutational signatures of more than 70% of human cancers. However, the mechanism underlying how cancer cells hijack the APOBEC mediated mutagenesis machinery to promote tumor heterogeneity, and thereby foster therapy resistance remains unclear. We identify SYNCRIP as an endogenous molecular brake which suppresses APOBEC-driven mutagenesis in prostate cancer (PCa). Overactivated APOBEC3B, in SYNCRIP-deficient PCa cells, is a key mutator, representing the molecular source of driver mutations in some frequently mutated genes in PCa, including FOXA1, EP300. Functional screening identifies eight crucial drivers for androgen receptor (AR)-targeted therapy resistance in PCa that are mutated by APOBEC3B: BRD7, CBX8, EP300, FOXA1, HDAC5, HSF4, STAT3, and AR. These results uncover a cell-intrinsic mechanism that unleashes APOBEC-driven mutagenesis, which plays a significant role in conferring AR-targeted therapy resistance in PCa.

Clonal hematopoiesis and inflammation - the perpetual cycle

Acquired genetic or cytogenetic alterations in a blood stem cell that confer clonal fitness promote its relative expansion leading to clonal hematopoiesis (CH). Despite a largely intact hematopoietic output, CH is associated with a heightened risk of progression to hematologic malignancies and with non-hematologic health manifestations, including cardiovascular disease and overall mortality. We focus on the evidence for the role of inflammation in establishing, maintaining and reciprocally being affected by CH. We describe the known pro-inflammatory signals associated with CH and preclinical studies that elucidated the cellular mechanisms involved. We review the evolving literature on early-onset CH in germline predisposition conditions and the possible role of immune dysregulation in this context.

Clonal hematopoiesis is associated with protection from Alzheimer's disease

Clonal hematopoiesis of indeterminate potential (CHIP) is a premalignant expansion of mutated hematopoietic stem cells. As CHIP-associated mutations are known to alter the development and function of myeloid cells, we hypothesized that CHIP may also be associated with the risk of Alzheimer's disease (AD), a disease in which brain-resident myeloid cells are thought to have a major role. To perform association tests between CHIP and AD dementia, we analyzed blood DNA sequencing data from 1,362 individuals with AD and 4,368 individuals without AD. Individuals with CHIP had a lower risk of AD dementia (meta-analysis odds ratio (OR) = 0.64, P = 3.8 × 10-5), and Mendelian randomization analyses supported a potential causal association. We observed that the same mutations found in blood were also detected in microglia-enriched fraction of the brain in seven of eight CHIP carriers. Single-nucleus chromatin accessibility profiling of brain-derived nuclei in six CHIP carriers revealed that the mutated cells comprised a large proportion of the microglial pool in the samples examined. While additional studies are required to validate the mechanistic findings, these results suggest that CHIP may have a role in attenuating the risk of AD.

Cell type-specific role of CBX2 and its disordered region in spermatogenesis

Polycomb group (PcG) proteins maintain the repressed state of lineage-inappropriate genes and are therefore essential for embryonic development and adult tissue homeostasis. One critical function of PcG complexes is modulating chromatin structure. Canonical Polycomb repressive complex 1 (cPRC1), particularly its component CBX2, can compact chromatin and phase-separate in vitro. These activities are hypothesized to be critical for forming a repressed physical environment in cells. While much has been learned by studying these PcG activities in cell culture models, it is largely unexplored how cPRC1 regulates adult stem cells and their subsequent differentiation in living animals. Here, we show in vivo evidence of a critical nonenzymatic repressive function of cPRC1 component CBX2 in the male germline. CBX2 is up-regulated as spermatogonial stem cells differentiate and is required to repress genes that were active in stem cells. CBX2 forms condensates (similar to previously described Polycomb bodies) that colocalize with target genes bound by CBX2 in differentiating spermatogonia. Single-cell analyses of mosaic Cbx2 mutant testes show that CBX2 is specifically required to produce differentiating A1 spermatogonia. Furthermore, the region of CBX2 responsible for compaction and phase separation is needed for the long-term maintenance of male germ cells in the animal. These results emphasize that the regulation of chromatin structure by CBX2 at a specific stage of spermatogenesis is critical, which distinguishes this from a mechanism that is reliant on histone modification.