Integrated Single-Cell Analysis Maps the Continuous Regulatory Landscape of Human Hematopoietic Differentiation

A multi-omic single-cell landscape of cellular diversification in the developing human cerebral cortex

The vast neuronal diversity in the human neocortex is vital for high-order brain functions, necessitating elucidation of the regulatory mechanisms underlying such unparalleled diversity. However, recent studies have yet to comprehensively reveal the diversity of neurons and the molecular logic of neocortical origin in humans at single-cell resolution through profiling transcriptomic or epigenomic landscapes, owing to the application of unimodal data alone to depict exceedingly heterogeneous populations of neurons. In this study, we generated a comprehensive compendium of the developing human neocortex by simultaneously profiling gene expression and open chromatin from the same cell. We computationally reconstructed the differentiation trajectories of excitatory projection neurons of cortical origin and inferred the regulatory logic governing lineage bifurcation decisions for neuronal diversification. We demonstrated that neuronal diversity arises from progenitor cell lineage specificity and postmitotic differentiation at distinct stages. Our data paves the way for understanding the primarily coordinated regulatory logic for neuronal diversification in the neocortex.

Inferring cell differentiation maps from lineage tracing data

During development, mulitpotent cells differentiate through a hierarchy of increasingly restricted progenitor cell types until they realize specialized cell types. A cell differentiation map describes this hierarchy, and inferring these maps is an active area of research spanning traditional single marker lineage studies to data-driven trajectory inference methods on single-cell RNA-seq data. Recent high-throughput lineage tracing technologies profile lineages and cell types at scale, but current methods to infer cell differentiation maps from these data rely on simple models with restrictive assumptions about the developmental process. We introduce a mathematical framework for cell differentiation maps based on the concept of potency, and develop an algorithm, Carta , that infers an optimal cell differentiation map from single-cell lineage tracing data. The key insight in Carta is to balance the trade-off between the complexity of the cell differentiation map and the number of unobserved cell type transitions on the lineage tree. We show that Carta more accurately infers cell differentiation maps on both simulated and real data compared to existing methods. In models of mammalian trunk development and mouse hematopoiesis, Carta identifies important features of development that are not revealed by other methods including convergent differentiation of specialized cell types, progenitor differentiation dynamics, and the refinement of routes of differentiation via new intermediate progenitors.

Code availability

Carta software is available at https://github.com/raphael-group/CARTA.

Epigenetic regulation of hematopoietic stem cell fate

Hematopoietic stem cells (HSCs) sustain blood cell production throughout the mammalian life span. However, it has become clear that at the single cell level a subset of HSCs is stably biased in their lineage output, and that such heterogeneity may play a key role in physiological processes including aging and adaptive immunity. Analysis of chromatin accessibility, DNA methylation, and histone modifications has revealed that HSCs with different lineage bias exhibit distinct epigenetic traits inscribed at poised, lineage-specific enhancers. This allows for lineage priming without initiating lineage-specific gene expression in HSCs, controlling lineage bias while preserving self-renewal and multipotency. Here, we review our current understanding of epigenetic regulation in the establishment and maintenance of HSC fate decisions under different physiological conditions.

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.

Heterotypic interaction promotes asymmetric division of human hematopoietic progenitors

Hematopoietic stem and progenitor cells (HSPCs) give rise to all cell types of the hematopoietic system through various processes, including asymmetric divisions. However, the contribution of stromal cells of the hematopoietic niches in the control of HSPC asymmetric divisions remains unknown. Using polyacrylamide microwells as minimalist niches, we show that specific heterotypic interactions with osteoblast and endothelial cells promote asymmetric divisions of human HSPCs. Upon interaction, HSPCs polarize in interphase with the centrosome, the Golgi apparatus, and lysosomes positioned close to the site of contact. Subsequently, during mitosis, HSPCs orient their spindle perpendicular to the plane of contact. This division mode gives rise to siblings with unequal amounts of lysosomes and of the differentiation marker CD34. Such asymmetric inheritance generates heterogeneity in the progeny, which is likely to contribute to the plasticity of the early steps of hematopoiesis.

Single-nucleus transcriptomics and chromatin accessibility analysis of musk gland development in Chinese forest musk deer (Moschus berezovskii)

Musk secreted by male forest musk deer (Moschus berezovskii) musk glands is an invaluable component of medicine and perfume. Musk secretion depends on musk gland maturation; however, the mechanism of its development remains elusive. Herein, using single cell multiome ATAC + gene expression coupled with several bioinformatic analyses, a dynamic transcriptional cell atlas of musk gland development was revealed, and key genes and transcription factors affecting its development were determined. Twelve cell types, including two different types of acinar cells (Clusters 0 and 10) were identified. Single-nucleus RNA and single-nucleus ATAC sequencing analyses revealed that seven core target genes associated with musk secretion (Hsd17b2, Acacb, Lss, Vapa, Aldh16a1, Aldh7a1, and Sqle) were regulated by 12 core transcription factors (FOXO1, CUX2, RORA, RUNX1, KLF6, MGA, NFIC, FOXO3, ETV5, NR3C1, HSF4, and MITF) during the development of Cluster 0 acinar cells. Kyoto Encyclopedia of Genes and Genomes enrichment showed significant changes in the pathways associated with musk secretion during acinar cell development. Gene set variation analysis also revealed that certain pathways associated with musk secretion were enriched in 6-year-old acinar cells. A gene co-expression network was constructed during acinar cell development to provide a precise understanding of the connections between transcription factors, genes, and pathways. Finally, intercellular communication analysis showed that intercellular communication is involved in musk gland development. This study provides crucial insights into the changes and key factors underlying musk gland development, which serve as valuable resources for studying musk secretion mechanisms and promoting the protection of this endangered species.

Fast clustering and cell-type annotation of scATAC data using pre-trained embeddings

Data from the single-cell assay for transposase-accessible chromatin using sequencing (scATAC-seq) are now widely available. One major computational challenge is dealing with high dimensionality and inherent sparsity, which is typically addressed by producing lower dimensional representations of single cells for downstream clustering tasks. Current approaches produce such individual cell embeddings directly through a one-step learning process. Here, we propose an alternative approach by building embedding models pre-trained on reference data. We argue that this provides a more flexible analysis workflow that also has computational performance advantages through transfer learning. We implemented our approach in scEmbed, an unsupervised machine-learning framework that learns low-dimensional embeddings of genomic regulatory regions to represent and analyze scATAC-seq data. scEmbed performs well in terms of clustering ability and has the key advantage of learning patterns of region co-occurrence that can be transferred to other, unseen datasets. Moreover, models pre-trained on reference data can be exploited to build fast and accurate cell-type annotation systems without the need for other data modalities. scEmbed is implemented in Python and it is available to download from GitHub. We also make our pre-trained models available on huggingface for public use. scEmbed is open source and available at https://github.com/databio/geniml. Pre-trained models from this work can be obtained on huggingface: https://huggingface.co/databio.

RNA sequestration in P-bodies sustains myeloid leukaemia

Post-transcriptional mechanisms are fundamental safeguards of progenitor cell identity and are often dysregulated in cancer. Here, we identified regulators of P-bodies as crucial vulnerabilities in acute myeloid leukaemia (AML) through genome-wide CRISPR screens in normal and malignant haematopoietic progenitors. We found that leukaemia cells harbour aberrantly elevated numbers of P-bodies and show that P-body assembly is crucial for initiation and maintenance of AML. Notably, P-body loss had little effect upon homoeostatic haematopoiesis but impacted regenerative haematopoiesis. Molecular characterization of P-bodies purified from human AML cells unveiled their critical role in sequestering messenger RNAs encoding potent tumour suppressors from the translational machinery. P-body dissolution promoted translation of these mRNAs, which in turn rewired gene expression and chromatin architecture in leukaemia cells. Collectively, our findings highlight the contrasting and unique roles of RNA sequestration in P-bodies during tissue homoeostasis and oncogenesis. These insights open potential avenues for understanding myeloid leukaemia and future therapeutic interventions.

Profiling cell identity and tissue architecture with single-cell and spatial transcriptomics

Single-cell transcriptomics has broadened our understanding of cellular diversity and gene expression dynamics in healthy and diseased tissues. Recently, spatial transcriptomics has emerged as a tool to contextualize single cells in multicellular neighbourhoods and to identify spatially recurrent phenotypes, or ecotypes. These technologies have generated vast datasets with targeted-transcriptome and whole-transcriptome profiles of hundreds to millions of cells. Such data have provided new insights into developmental hierarchies, cellular plasticity and diverse tissue microenvironments, and spurred a burst of innovation in computational methods for single-cell analysis. In this Review, we discuss recent advancements, ongoing challenges and prospects in identifying and characterizing cell states and multicellular neighbourhoods. We discuss recent progress in sample processing, data integration, identification of subtle cell states, trajectory modelling, deconvolution and spatial analysis. Furthermore, we discuss the increasing application of deep learning, including foundation models, in analysing single-cell and spatial transcriptomics data. Finally, we discuss recent applications of these tools in the fields of stem cell biology, immunology, and tumour biology, and the future of single-cell and spatial transcriptomics in biological research and its translation to the clinic.

Benchmarking computational methods for single-cell chromatin data analysis

BACKGROUND

Single-cell chromatin accessibility assays, such as scATAC-seq, are increasingly employed in individual and joint multi-omic profiling of single cells. As the accumulation of scATAC-seq and multi-omics datasets continue, challenges in analyzing such sparse, noisy, and high-dimensional data become pressing. Specifically, one challenge relates to optimizing the processing of chromatin-level measurements and efficiently extracting information to discern cellular heterogeneity. This is of critical importance, since the identification of cell types is a fundamental step in current single-cell data analysis practices.

RESULTS

We benchmark 8 feature engineering pipelines derived from 5 recent methods to assess their ability to discover and discriminate cell types. By using 10 metrics calculated at the cell embedding, shared nearest neighbor graph, or partition levels, we evaluate the performance of each method at different data processing stages. This comprehensive approach allows us to thoroughly understand the strengths and weaknesses of each method and the influence of parameter selection.

CONCLUSIONS

Our analysis provides guidelines for choosing analysis methods for different datasets. Overall, feature aggregation, SnapATAC, and SnapATAC2 outperform latent semantic indexing-based methods. For datasets with complex cell-type structures, SnapATAC and SnapATAC2 are preferred. With large datasets, SnapATAC2 and ArchR are most scalable.

A Cell Cycle-Aware Network for Data Integration and Label Transferring of Single-Cell RNA-Seq and ATAC-Seq

In recent years, the integration of single-cell multi-omics data has provided a more comprehensive understanding of cell functions and internal regulatory mechanisms from a non-single omics perspective, but it still suffers many challenges, such as omics-variance, sparsity, cell heterogeneity, and confounding factors. As it is known, the cell cycle is regarded as a confounder when analyzing other factors in single-cell RNA-seq data, but it is not clear how it will work on the integrated single-cell multi-omics data. Here, a cell cycle-aware network (CCAN) is developed to remove cell cycle effects from the integrated single-cell multi-omics data while keeping the cell type-specific variations. This is the first computational model to study the cell-cycle effects in the integration of single-cell multi-omics data. Validations on several benchmark datasets show the outstanding performance of CCAN in a variety of downstream analyses and applications, including removing cell cycle effects and batch effects of scRNA-seq datasets from different protocols, integrating paired and unpaired scRNA-seq and scATAC-seq data, accurately transferring cell type labels from scRNA-seq to scATAC-seq data, and characterizing the differentiation process from hematopoietic stem cells to different lineages in the integration of differentiation data.

Terminal deoxynucleotidyl transferase and CD84 identify human multi-potent lymphoid progenitors

Lymphoid specification in human hematopoietic progenitors is not fully understood. To better associate lymphoid identity with protein-level cell features, we conduct a highly multiplexed single-cell proteomic screen on human bone marrow progenitors. This screen identifies terminal deoxynucleotidyl transferase (TdT), a specialized DNA polymerase intrinsic to VDJ recombination, broadly expressed within CD34+ progenitors prior to B/T cell emergence. While these TdT+ cells coincide with granulocyte-monocyte progenitor (GMP) immunophenotype, their accessible chromatin regions show enrichment for lymphoid-associated transcription factor (TF) motifs. TdT expression on GMPs is inversely related to the SLAM family member CD84. Prospective isolation of CD84lo GMPs demonstrates robust lymphoid potentials ex vivo, while still retaining significant myeloid differentiation capacity, akin to LMPPs. This multi-omic study identifies human bone marrow lymphoid-primed progenitors, further defining the lympho-myeloid axis in human hematopoiesis.

Epigenetic alterations affecting hematopoietic regulatory networks as drivers of mixed myeloid/lymphoid leukemia

Leukemias with ambiguous lineage comprise several loosely defined entities, often without a clear mechanistic basis. Here, we extensively profile the epigenome and transcriptome of a subgroup of such leukemias with CpG Island Methylator Phenotype. These leukemias exhibit comparable hybrid myeloid/lymphoid epigenetic landscapes, yet heterogeneous genetic alterations, suggesting they are defined by their shared epigenetic profile rather than common genetic lesions. Gene expression enrichment reveals similarity with early T-cell precursor acute lymphoblastic leukemia and a lymphoid progenitor cell of origin. In line with this, integration of differential DNA methylation and gene expression shows widespread silencing of myeloid transcription factors. Moreover, binding sites for hematopoietic transcription factors, including CEBPA, SPI1 and LEF1, are uniquely inaccessible in these leukemias. Hypermethylation also results in loss of CTCF binding, accompanied by changes in chromatin interactions involving key transcription factors. In conclusion, epigenetic dysregulation, and not genetic lesions, explains the mixed phenotype of this group of leukemias with ambiguous lineage. The data collected here constitute a useful and comprehensive epigenomic reference for subsequent studies of acute myeloid leukemias, T-cell acute lymphoblastic leukemias and mixed-phenotype leukemias.

Benchmarking Algorithms for Gene Set Scoring of Single-cell ATAC-seq Data

Gene set scoring (GSS) has been routinely conducted for gene expression analysis of bulk or single-cell RNA sequencing (RNA-seq) data, which helps to decipher single-cell heterogeneity and cell type-specific variability by incorporating prior knowledge from functional gene sets. Single-cell assay for transposase accessible chromatin using sequencing (scATAC-seq) is a powerful technique for interrogating single-cell chromatin-based gene regulation, and genes or gene sets with dynamic regulatory potentials can be regarded as cell type-specific markers as if in single-cell RNA-seq (scRNA-seq). However, there are few GSS tools specifically designed for scATAC-seq, and the applicability and performance of RNA-seq GSS tools on scATAC-seq data remain to be investigated. Here, we systematically benchmarked ten GSS tools, including four bulk RNA-seq tools, five scRNA-seq tools, and one scATAC-seq method. First, using matched scATAC-seq and scRNA-seq datasets, we found that the performance of GSS tools on scATAC-seq data was comparable to that on scRNA-seq, suggesting their applicability to scATAC-seq. Then, the performance of different GSS tools was extensively evaluated using up to ten scATAC-seq datasets. Moreover, we evaluated the impact of gene activity conversion, dropout imputation, and gene set collections on the results of GSS. Results show that dropout imputation can significantly promote the performance of almost all GSS tools, while the impact of gene activity conversion methods or gene set collections on GSS performance is more dependent on GSS tools or datasets. Finally, we provided practical guidelines for choosing appropriate preprocessing methods and GSS tools in different application scenarios.

A single-cell chromatin accessibility dataset of human primed and naïve pluripotent stem cell-derived teratoma

Teratoma, due to its remarkable ability to differentiate into multiple cell lineages, is a valuable model for studying human embryonic development. The similarity of the gene expression and chromatin accessibility patterns in these cells to those observed in vivo further underscores its potential as a research tool. Notably, teratomas derived from human naïve (pre-implantation epiblast-like) pluripotent stem cells (PSCs) have larger embryonic cell diversity and contain extraembryonic lineages, making them more suitable to study developmental processes. However, the cell type-specific epigenetic profiles of naïve PSC teratomas have not been yet characterized. Using single-cell assay for transposase-accessible chromatin sequencing (scATAC-seq), we analyzed 66,384 cell profiles from five teratomas derived from human naïve PSCs and their post-implantation epiblast-like (primed) counterparts. We observed 17 distinct cell types from both embryonic and extraembryonic lineages, resembling the corresponding cell types in human fetal tissues. Additionally, we identified key transcription factors specific to different cell types. Our dataset provides a resource for investigating gene regulatory programs in a relevant model of human embryonic development.

Decoding human in vitro terminal erythropoiesis originating from umbilical cord blood mononuclear cells and pluripotent stem cells

Ex vivo red blood cell (RBC) production generates unsatisfactory erythroid cells. A deep exploration into terminally differentiated cells is required to understand the impairments for RBC generation and the underlying mechanisms. Here, we mapped an atlas of terminally differentiated cells from umbilical cord blood mononuclear cells (UCBMN) and pluripotent stem cells (PSC) and observed their dynamic regulation of erythropoiesis at single-cell resolution. Interestingly, we detected a few progenitor cells and non-erythroid cells from both origins. In PSC-derived erythropoiesis (PSCE), the expression of haemoglobin switch regulators (BCL11A and ZBTB7A) were significantly absent, which could be the restraint for its adult globin expression. We also found that PSCE were less active in stress erythropoiesis than in UCBMN-derived erythropoiesis (UCBE), and explored an agonist of stress erythropoiesis gene, TRIB3, could enhance the expression of adult globin in PSCE. Compared with UCBE, there was a lower expression of epigenetic-related proteins (e.g., CASPASE 3 and UBE2O) and transcription factors (e.g., FOXO3 and TAL1) in PSCE, which might restrict PSCE's enucleation. Moreover, we characterized a subpopulation with high proliferation capacity marked by CD99high in colony-forming unit-erythroid cells. Inhibition of CD99 reduced the proliferation of PSC-derived cells and facilitated erythroid maturation. Furthermore, CD99-CD99 mediated the interaction between macrophages and erythroid cells, illustrating a mechanism by which macrophages participate in erythropoiesis. This study provided a reference for improving ex vivo RBC generation.

Trisomy 8 Defines a Distinct Subtype of Myeloproliferative Neoplasms Driven by the MYC-Alarmin Axis

Despite advances in understanding the genetic abnormalities in myeloproliferative neoplasms (MPN) and the development of JAK2 inhibitors, there is an urgent need to devise new treatment strategies, particularly for patients with triple-negative (TN) myelofibrosis (MF) who lack mutations in the JAK2 kinase pathway and have very poor clinical outcomes. Here we report that MYC copy number gain and increased MYC expression frequently occur in TN-MF and that MYC-directed activation of S100A9, an alarmin protein that plays pivotal roles in inflammation and innate immunity, is necessary and sufficient to drive development and progression of MF. Notably, the MYC-S100A9 circuit provokes a complex network of inflammatory signaling that involves numerous hematopoietic cell types in the bone marrow microenvironment. Accordingly, genetic ablation of S100A9 or treatment with small molecules targeting the MYC-S100A9 pathway effectively ameliorates MF phenotypes, highlighting the MYC-alarmin axis as a novel therapeutic vulnerability for this subgroup of MPNs. Significance: This study establishes that MYC expression is increased in TN-MPNs via trisomy 8, that a MYC-S100A9 circuit manifest in these cases is sufficient to provoke myelofibrosis and inflammation in diverse hematopoietic cell types in the BM niche, and that the MYC-S100A9 circuit is targetable in TN-MPNs.

Harnessing Single-Cell Technologies in the Search for New Therapies for Diamond-Blackfan Anemia Syndrome

The emergence of multiomic single-cell technologies over the last decade has led to improved insights into both normal hematopoiesis and its perturbation in a variety of hematological disorders. Diamond-Blackfan anemia (DBA) syndrome is one such disorder where single-cell assays have helped to delineate the cellular and molecular defects underlying the disease. DBA is caused by heterozygous loss-of-function germline variants in genes encoding ribosomal proteins (RPs). Despite the widespread role of ribosomes in hematopoiesis, the most frequent and severe cytopenia in DBA is anemia. In this review we discussed how single-cell studies, including clonogenic cell culture assays, fluorescence-activated cell sorting (FACS) and single-cell RNA sequencing (scRNA-seq), have led to insights into the pathogenesis of DBA. The main therapies are regular blood transfusions, glucocorticoids, or hematopoietic stem cell transplantation (HSCT) but all are associated with significant morbidity and mortality. We will therefore outline how single-cell studies can inform new therapies for DBA. Furthermore, we discussed how DBA serves as a useful model for understanding normal erythropoiesis in terms of its cellular hierarchy, molecular regulation during homeostasis, and response to "stress."

Quantifying cell-state densities in single-cell phenotypic landscapes using Mellon

Cell-state density characterizes the distribution of cells along phenotypic landscapes and is crucial for unraveling the mechanisms that drive diverse biological processes. Here, we present Mellon, an algorithm for estimation of cell-state densities from high-dimensional representations of single-cell data. We demonstrate Mellon's efficacy by dissecting the density landscape of differentiating systems, revealing a consistent pattern of high-density regions corresponding to major cell types intertwined with low-density, rare transitory states. We present evidence implicating enhancer priming and the activation of master regulators in emergence of these transitory states. Mellon offers the flexibility to perform temporal interpolation of time-series data, providing a detailed view of cell-state dynamics during developmental processes. Mellon facilitates density estimation across various single-cell data modalities, scaling linearly with the number of cells. Our work underscores the importance of cell-state density in understanding the differentiation processes, and the potential of Mellon to provide insights into mechanisms guiding biological trajectories.

Epigenetics of hematopoietic stem cell aging

PURPOSE OF REVIEW

The development of new antiaging medicines is of great interest to the current elderly and aging population. Aging of the hematopoietic system is attributed to the aging of hematopoietic stem cells (HSCs), and epigenetic alterations are the key effectors driving HSC aging. Understanding the epigenetics of HSC aging holds promise of providing new insights for combating HSC aging and age-related hematological malignancies.

RECENT FINDINGS

Aging is characterized by the progressive loss of physiological integrity, leading to impaired function and increased vulnerability to death. During aging, the HSCs undergo both quantitative and qualitative changes. These functional changes in HSCs cause dysregulated hematopoiesis, resulting in anemia, immune dysfunction, and an increased risk of hematological malignancies. Various cell-intrinsic and cell-extrinsic effectors influencing HSC aging have also been identified. Epigenetic alterations are one such mechanism.

SUMMARY

Cumulative epigenetic alterations in aged HSCs affect their fate, leading to aberrant self-renewal, differentiation, and function of aged HSCs. In turn, these factors provide an opportunity for aged HSCs to expand by modulating their self-renewal and differentiation balance, thereby contributing to the development of hematological malignancies.

Targeting MCL1-driven anti-apoptotic pathways overcomes blast progression after hypomethylating agent failure in chronic myelomonocytic leukemia

RAS pathway mutations, which are present in 30% of patients with chronic myelomonocytic leukemia (CMML) at diagnosis, confer a high risk of resistance to and progression after hypomethylating agent (HMA) therapy, the current standard of care for the disease. Here, using single-cell, multi-omics technologies, we seek to dissect the biological mechanisms underlying the initiation and progression of RAS pathway-mutated CMML. We identify that RAS pathway mutations induce transcriptional reprogramming of hematopoietic stem and progenitor cells (HSPCs) and downstream monocytic populations in response to cell-intrinsic and -extrinsic inflammatory signaling that also impair the functions of immune cells. HSPCs expand at disease progression after therapy with HMA or the BCL2 inhibitor venetoclax and rely on the NF-κB pathway effector MCL1 to maintain survival. Our study has implications for the development of therapies to improve the survival of patients with RAS pathway-mutated CMML.

Non-viral DNA delivery and TALEN editing correct the sickle cell mutation in hematopoietic stem cells

Sickle cell disease is a devastating blood disorder that originates from a single point mutation in the HBB gene coding for hemoglobin. Here, we develop a GMP-compatible TALEN-mediated gene editing process enabling efficient HBB correction via a DNA repair template while minimizing risks associated with HBB inactivation. Comparing viral versus non-viral DNA repair template delivery in hematopoietic stem and progenitor cells in vitro, both strategies achieve comparable HBB correction and result in over 50% expression of normal adult hemoglobin in red blood cells without inducing β-thalassemic phenotype. In an immunodeficient female mouse model, transplanted cells edited with the non-viral strategy exhibit higher engraftment and gene correction levels compared to those edited with the viral strategy. Transcriptomic analysis reveals that non-viral DNA repair template delivery mitigates P53-mediated toxicity and preserves high levels of long-term hematopoietic stem cells. This work paves the way for TALEN-based autologous gene therapy for sickle cell disease.

Single-cell lineage capture across genomic modalities with CellTag-multi reveals fate-specific gene regulatory changes

Complex gene regulatory mechanisms underlie differentiation and reprogramming. Contemporary single-cell lineage-tracing (scLT) methods use expressed, heritable DNA barcodes to combine cell lineage readout with single-cell transcriptomics. However, reliance on transcriptional profiling limits adaptation to other single-cell assays. With CellTag-multi, we present an approach that enables direct capture of heritable random barcodes expressed as polyadenylated transcripts, in both single-cell RNA sequencing and single-cell Assay for Transposase Accessible Chromatin using sequencing assays, allowing for independent clonal tracking of transcriptional and epigenomic cell states. We validate CellTag-multi to characterize progenitor cell lineage priming during mouse hematopoiesis. Additionally, in direct reprogramming of fibroblasts to endoderm progenitors, we identify core regulatory programs underlying on-target and off-target fates. Furthermore, we reveal the transcription factor Zfp281 as a regulator of reprogramming outcome, biasing cells toward an off-target mesenchymal fate. Our results establish CellTag-multi as a lineage-tracing method compatible with multiple single-cell modalities and demonstrate its utility in revealing fate-specifying gene regulatory changes across diverse paradigms of differentiation and reprogramming.

SIMBA: single-cell embedding along with features

Most current single-cell analysis pipelines are limited to cell embeddings and rely heavily on clustering, while lacking the ability to explicitly model interactions between different feature types. Furthermore, these methods are tailored to specific tasks, as distinct single-cell problems are formulated differently. To address these shortcomings, here we present SIMBA, a graph embedding method that jointly embeds single cells and their defining features, such as genes, chromatin-accessible regions and DNA sequences, into a common latent space. By leveraging the co-embedding of cells and features, SIMBA allows for the study of cellular heterogeneity, clustering-free marker discovery, gene regulation inference, batch effect removal and omics data integration. We show that SIMBA provides a single framework that allows diverse single-cell problems to be formulated in a unified way and thus simplifies the development of new analyses and extension to new single-cell modalities. SIMBA is implemented as a comprehensive Python library ( https://simba-bio.readthedocs.io ).

Scalable and unbiased sequence-informed embedding of single-cell ATAC-seq data with CellSpace

Standard scATAC sequencing (scATAC-seq) analysis pipelines represent cells as sparse numeric vectors relative to an atlas of peaks or genomic tiles and consequently ignore genomic sequence information at accessible loci. Here we present CellSpace, an efficient and scalable sequence-informed embedding algorithm for scATAC-seq that learns a mapping of DNA k-mers and cells to the same space, to address this limitation. We show that CellSpace captures meaningful latent structure in scATAC-seq datasets, including cell subpopulations and developmental hierarchies, and can score transcription factor activities in single cells based on proximity to binding motifs embedded in the same space. Importantly, CellSpace implicitly mitigates batch effects arising from multiple samples, donors or assays, even when individual datasets are processed relative to different peak atlases. Thus, CellSpace provides a powerful tool for integrating and interpreting large-scale scATAC-seq compendia.

Clock-dependent chromatin accessibility rhythms regulate circadian transcription

Chromatin organization plays a crucial role in gene regulation by controlling the accessibility of DNA to transcription machinery. While significant progress has been made in understanding the regulatory role of clock proteins in circadian rhythms, how chromatin organization affects circadian rhythms remains poorly understood. Here, we employed ATAC-seq (Assay for Transposase-Accessible Chromatin with Sequencing) on FAC-sorted Drosophila clock neurons to assess genome-wide chromatin accessibility at dawn and dusk over the circadian cycle. We observed significant oscillations in chromatin accessibility at promoter and enhancer regions of hundreds of genes, with enhanced accessibility either at dusk or dawn, which correlated with their peak transcriptional activity. Notably, genes with enhanced accessibility at dusk were enriched with E-box motifs, while those more accessible at dawn were enriched with VRI/PDP1-box motifs, indicating that they are regulated by the core circadian feedback loops, PER/CLK and VRI/PDP1, respectively. Further, we observed a complete loss of chromatin accessibility rhythms in per01 null mutants, with chromatin consistently accessible at both dawn and dusk, underscoring the critical role of Period protein in driving chromatin compaction during the repression phase at dawn. Together, this study demonstrates the significant role of chromatin organization in circadian regulation, revealing how the interplay between clock proteins and chromatin structure orchestrates the precise timing of biological processes throughout the day. This work further implies that variations in chromatin accessibility might play a central role in the generation of diverse circadian gene expression patterns in clock neurons.

Single-cell multi-omics analysis identifies context-specific gene regulatory gates and mechanisms

There is a growing interest in inferring context specific gene regulatory networks from single-cell RNA sequencing (scRNA-seq) data. This involves identifying the regulatory relationships between transcription factors (TFs) and genes in individual cells, and then characterizing these relationships at the level of specific cell types or cell states. In this study, we introduce scGATE (single-cell gene regulatory gate) as a novel computational tool for inferring TF-gene interaction networks and reconstructing Boolean logic gates involving regulatory TFs using scRNA-seq data. In contrast to current Boolean models, scGATE eliminates the need for individual formulations and likelihood calculations for each Boolean rule (e.g. AND, OR, XOR). By employing a Bayesian framework, scGATE infers the Boolean rule after fitting the model to the data, resulting in significant reductions in time-complexities for logic-based studies. We have applied assay for transposase-accessible chromatin with sequencing (scATAC-seq) data and TF DNA binding motifs to filter out non-relevant TFs in gene regulations. By integrating single-cell clustering with these external cues, scGATE is able to infer context specific networks. The performance of scGATE is evaluated using synthetic and real single-cell multi-omics data from mouse tissues and human blood, demonstrating its superiority over existing tools for reconstructing TF-gene networks. Additionally, scGATE provides a flexible framework for understanding the complex combinatorial and cooperative relationships among TFs regulating target genes by inferring Boolean logic gates among them.

Hematopoietic stem cells with granulo-monocytic differentiation state overcome venetoclax sensitivity in patients with myelodysplastic syndromes

The molecular mechanisms of venetoclax-based therapy failure in patients with acute myeloid leukemia were recently clarified, but the mechanisms by which patients with myelodysplastic syndromes (MDS) acquire secondary resistance to venetoclax after an initial response remain to be elucidated. Here, we show an expansion of MDS hematopoietic stem cells (HSCs) with a granulo-monocytic-biased transcriptional differentiation state in MDS patients who initially responded to venetoclax but eventually relapsed. While MDS HSCs in an undifferentiated cellular state are sensitive to venetoclax treatment, differentiation towards a granulo-monocytic-biased transcriptional state, through the acquisition or expansion of clones with STAG2 or RUNX1 mutations, affects HSCs' survival dependence from BCL2-mediated anti-apoptotic pathways to TNFα-induced pro-survival NF-κB signaling and drives resistance to venetoclax-mediated cytotoxicity. Our findings reveal how hematopoietic stem and progenitor cell (HSPC) can eventually overcome therapy-induced depletion and underscore the importance of using close molecular monitoring to prevent HSPC hierarchical change in MDS patients enrolled in clinical trials of venetoclax.

Machine Learning and Deep Learning in Synthetic Biology: Key Architectures, Applications, and Challenges

Machine learning (ML), particularly deep learning (DL), has made rapid and substantial progress in synthetic biology in recent years. Biotechnological applications of biosystems, including pathways, enzymes, and whole cells, are being probed frequently with time. The intricacy and interconnectedness of biosystems make it challenging to design them with the desired properties. ML and DL have a synergy with synthetic biology. Synthetic biology can be employed to produce large data sets for training models (for instance, by utilizing DNA synthesis), and ML/DL models can be employed to inform design (for example, by generating new parts or advising unrivaled experiments to perform). This potential has recently been brought to light by research at the intersection of engineering biology and ML/DL through achievements like the design of novel biological components, best experimental design, automated analysis of microscopy data, protein structure prediction, and biomolecular implementations of ANNs (Artificial Neural Networks). I have divided this review into three sections. In the first section, I describe predictive potential and basics of ML along with myriad applications in synthetic biology, especially in engineering cells, activity of proteins, and metabolic pathways. In the second section, I describe fundamental DL architectures and their applications in synthetic biology. Finally, I describe different challenges causing hurdles in the progress of ML/DL and synthetic biology along with their solutions.

scAuto as a comprehensive framework for single-cell chromatin accessibility data analysis

Interpreting single-cell chromatin accessibility data is crucial for understanding intercellular heterogeneity regulation. Despite the progress in computational methods for analyzing this data, there is still a lack of a comprehensive analytical framework and a user-friendly online analysis tool. To fill this gap, we developed a pre-trained deep learning-based framework, single-cell auto-correlation transformers (scAuto), to overcome the challenge. Following DNABERT's methodology of pre-training and fine-tuning, scAuto learns a general understanding of DNA sequence's grammar by being pre-trained on unlabeled human genome via self-supervision; it is then transferred to the single-cell chromatin accessibility analysis task of scATAC-seq data for supervised fine-tuning. We extensively validated scAuto on the Buenrostro2018 dataset, demonstrating its superior performance on chromatin accessibility prediction, single-cell clustering, and data denoising. Based on scAuto, we further developed an interactive web server for single-cell chromatin accessibility data analysis. It integrates tutorial-style interfaces for those with limited programming skills. The platform is accessible at http://zhanglab.icaup.cn. To our knowledge, this work is expected to help analyze single-cell chromatin accessibility data and facilitate the development of precision medicine.

scCASE: accurate and interpretable enhancement for single-cell chromatin accessibility sequencing data

Single-cell chromatin accessibility sequencing (scCAS) has emerged as a valuable tool for interrogating and elucidating epigenomic heterogeneity and gene regulation. However, scCAS data inherently suffers from limitations such as high sparsity and dimensionality, which pose significant challenges for downstream analyses. Although several methods are proposed to enhance scCAS data, there are still challenges and limitations that hinder the effectiveness of these methods. Here, we propose scCASE, a scCAS data enhancement method based on non-negative matrix factorization which incorporates an iteratively updating cell-to-cell similarity matrix. Through comprehensive experiments on multiple datasets, we demonstrate the advantages of scCASE over existing methods for scCAS data enhancement. The interpretable cell type-specific peaks identified by scCASE can provide valuable biological insights into cell subpopulations. Moreover, to leverage the large compendia of available omics data as a reference, we further expand scCASE to scCASER, which enables the incorporation of external reference data to improve enhancement performance.

Cross-Omic Transcription Factor Analysis: An Insight on Transcription Factor Accessibility and Expression Correlation

It is well known how sequencing technologies propelled cellular biology research in recent years, providing incredible insight into the basic mechanisms of cells. Single-cell RNA sequencing is at the front in this field, with single-cell ATAC sequencing supporting it and becoming more popular. In this regard, multi-modal technologies play a crucial role, allowing the possibility to simultaneously perform the mentioned sequencing modalities on the same cells. Yet, there still needs to be a clear and dedicated way to analyze these multi-modal data. One of the current methods is to calculate the Gene Activity Matrix (GAM), which summarizes the accessibility of the genes at the genomic level, to have a more direct link with the transcriptomic data. However, this concept is not well defined, and it is unclear how various accessible regions impact the expression of the genes. Moreover, the transcription process is highly regulated by the transcription factors that bind to the different DNA regions. Therefore, this work presents a continuation of the meta-analysis of Genomic-Annotated Gene Activity Matrix (GAGAM) contributions, aiming to investigate the correlation between the TF expression and motif information in the different functional genomic regions to understand the different Transcription Factors (TFs) dynamics involved in different cell types.

A synthetic differentiation circuit in Escherichia coli for suppressing mutant takeover

Differentiation is crucial for multicellularity. However, it is inherently susceptible to mutant cells that fail to differentiate. These mutants outcompete normal cells by excessive self-renewal. It remains unclear what mechanisms can resist such mutant expansion. Here, we demonstrate a solution by engineering a synthetic differentiation circuit in Escherichia coli that selects against these mutants via a biphasic fitness strategy. The circuit provides tunable production of synthetic analogs of stem, progenitor, and differentiated cells. It resists mutations by coupling differentiation to the production of an essential enzyme, thereby disadvantaging non-differentiating mutants. The circuit selected for and maintained a positive differentiation rate in long-term evolution. Surprisingly, this rate remained constant across vast changes in growth conditions. We found that transit-amplifying cells (fast-growing progenitors) underlie this environmental robustness. Our results provide insight into the stability of differentiation and demonstrate a powerful method for engineering evolutionarily stable multicellular consortia.

scMGCN: A Multi-View Graph Convolutional Network for Cell Type Identification in scRNA-seq Data

Single-cell RNA sequencing (scRNA-seq) data reveal the complexity and diversity of cellular ecosystems and molecular interactions in various biomedical research. Hence, identifying cell types from large-scale scRNA-seq data using existing annotations is challenging and requires stable and interpretable methods. However, the current cell type identification methods have limited performance, mainly due to the intrinsic heterogeneity among cell populations and extrinsic differences between datasets. Here, we present a robust graph artificial intelligence model, a multi-view graph convolutional network model (scMGCN) that integrates multiple graph structures from raw scRNA-seq data and applies graph convolutional networks with attention mechanisms to learn cell embeddings and predict cell labels. We evaluate our model on single-dataset, cross-species, and cross-platform experiments and compare it with other state-of-the-art methods. Our results show that scMGCN outperforms the other methods regarding stability, accuracy, and robustness to batch effects. Our main contributions are as follows: Firstly, we introduce multi-view learning and multiple graph construction methods to capture comprehensive cellular information from scRNA-seq data. Secondly, we construct a scMGCN that combines graph convolutional networks with attention mechanisms to extract shared, high-order information from cells. Finally, we demonstrate the effectiveness and superiority of the scMGCN on various datasets.

A Cell Cycle-aware Network for Data Integration and Label Transferring of Single-cell RNA-seq and ATAC-seq

In recent years, the integration of single-cell multi-omics data has provided a more comprehensive understanding of cell functions and internal regulatory mechanisms from a non-single omics perspective, but it still suffers many challenges, such as omics-variance, sparsity, cell heterogeneity and confounding factors. As we know, cell cycle is regarded as a confounder when analyzing other factors in single-cell RNA-seq data, but it's not clear how it will work on the integrated single-cell multi-omics data. Here, we developed a Cell Cycle-Aware Network (CCAN) to remove cell cycle effects from the integrated single-cell multi-omics data while keeping the cell type-specific variations. This is the first computational model to study the cell-cycle effects in the integration of single-cell multi-omics data. Validations on several benchmark datasets show the out-standing performance of CCAN in a variety of downstream analyses and applications, including removing cell cycle effects and batch effects of scRNA-seq datasets from different protocols, integrating paired and unpaired scRNA-seq and scATAC-seq data, accurately transferring cell type labels from scRNA-seq to scATAC-seq data, and characterizing the differentiation process from hematopoietic stem cells to different lineages in the integration of differentiation data.

A fast, scalable and versatile tool for analysis of single-cell omics data

Single-cell omics technologies have revolutionized the study of gene regulation in complex tissues. A major computational challenge in analyzing these datasets is to project the large-scale and high-dimensional data into low-dimensional space while retaining the relative relationships between cells. This low dimension embedding is necessary to decompose cellular heterogeneity and reconstruct cell-type-specific gene regulatory programs. Traditional dimensionality reduction techniques, however, face challenges in computational efficiency and in comprehensively addressing cellular diversity across varied molecular modalities. Here we introduce a nonlinear dimensionality reduction algorithm, embodied in the Python package SnapATAC2, which not only achieves a more precise capture of single-cell omics data heterogeneities but also ensures efficient runtime and memory usage, scaling linearly with the number of cells. Our algorithm demonstrates exceptional performance, scalability and versatility across diverse single-cell omics datasets, including single-cell assay for transposase-accessible chromatin using sequencing, single-cell RNA sequencing, single-cell Hi-C and single-cell multi-omics datasets, underscoring its utility in advancing single-cell analysis.

The cellular landscape of the endochondral bone during the transition to extrauterine life

The cellular complexity of the endochondral bone underlies its essential and pleiotropic roles during organismal life. While the adult bone has received significant attention, we still lack a deep understanding of the perinatal bone cellulome. Here, we have profiled the full composition of the murine endochondral bone at the single-cell level during the transition from fetal to newborn life and in comparison with the adult tissue, with particular emphasis on the mesenchymal compartment. The perinatal bone contains different fibroblastic clusters with blastema-like characteristics in organizing and supporting skeletogenesis, angiogenesis and hematopoiesis. Our data also suggest dynamic inter- and intra-compartment interactions, as well as a bone marrow milieu that seems prone to anti-inflammation, which we hypothesize is necessary to ensure the proper program of lymphopoiesis and the establishment of central and peripheral tolerance in early life. Our study provides an integrative roadmap for the future design of genetic and cellular functional assays to validate cellular interactions and lineage relationships within the perinatal bone.

Cellular carcinogenesis in preleukemic conditions:drivers and defenses

Acute myeloid leukemia (AML) arises from preleukemic conditions. We have investigated the pathogenesis of typical preleukemia, myeloproliferative neoplasms, and clonal hematopoiesis. Hematopoietic stem cells in both preleukemic conditions harbor recurrent driver mutations; additional mutation provokes further malignant transformation, leading to AML onset. Although genetic alterations are defined as the main cause of malignant transformation, non-genetic factors are also involved in disease progression. In this review, we focus on a non-histone chromatin protein, high mobility group AT-hook2 (HMGA2), and a physiological p53 inhibitor, murine double minute X (MDMX). HMGA2 is mainly overexpressed by dysregulation of microRNAs or mutations in polycomb components, and provokes expansion of preleukemic clones through stem cell signature disruption. MDMX is overexpressed by altered splicing balance in myeloid malignancies. MDMX induces leukemic transformation from preleukemia via suppression of p53 and p53-independent activation of WNT/β-catenin signaling. We also discuss how these non-genetic factors can be targeted for leukemia prevention therapy.

Chromatin accessibility and cell cycle progression are controlled by the HDAC-associated Sin3B protein in murine hematopoietic stem cells

BACKGROUND

Blood homeostasis requires the daily production of millions of terminally differentiated effector cells that all originate from hematopoietic stem cells (HSCs). HSCs are rare and exhibit unique self-renewal and multipotent properties, which depend on their ability to maintain quiescence through ill-defined processes. Defective control of cell cycle progression can eventually lead to bone marrow failure or malignancy. In particular, the molecular mechanisms tying cell cycle re-entry to cell fate commitment in HSCs remain elusive. Previous studies have identified chromatin coordination as a key regulator of differentiation in embryonic stem cells.

RESULTS

Here, we utilized genetic inactivation of the chromatin-associated Sin3B protein to manipulate cell cycle control and found dysregulated chromatin accessibility and cell cycle progression in HSCs. Single cell transcriptional profiling of hematopoietic stem and progenitor cells (HSPCs) inactivated for Sin3B reveals aberrant progression through the G1 phase of the cell cycle, which correlates with the engagement of specific signaling pathways, including aberrant expression of cell adhesion molecules and the interferon signaling program in LT-HSCs. In addition, we uncover the Sin3B-dependent accessibility of genomic elements controlling HSC differentiation, which points to cell cycle progression possibly dictating the priming of HSCs for differentiation.

CONCLUSIONS

Our findings provide new insights into controlled cell cycle progression as a potential regulator of HSC lineage commitment through the modulation of chromatin features.

Incorporating network diffusion and peak location information for better single-cell ATAC-seq data analysis

Single-cell assay for transposase-accessible chromatin using sequencing (scATAC-seq) data provided new insights into the understanding of epigenetic heterogeneity and transcriptional regulation. With the increasing abundance of dataset resources, there is an urgent need to extract more useful information through high-quality data analysis methods specifically designed for scATAC-seq. However, analyzing scATAC-seq data poses challenges due to its near binarization, high sparsity and ultra-high dimensionality properties. Here, we proposed a novel network diffusion-based computational method to comprehensively analyze scATAC-seq data, named Single-Cell ATAC-seq Analysis via Network Refinement with Peaks Location Information (SCARP). SCARP formulates the Network Refinement diffusion method under the graph theory framework to aggregate information from different network orders, effectively compensating for missing signals in the scATAC-seq data. By incorporating distance information between adjacent peaks on the genome, SCARP also contributes to depicting the co-accessibility of peaks. These two innovations empower SCARP to obtain lower-dimensional representations for both cells and peaks more effectively. We have demonstrated through sufficient experiments that SCARP facilitated superior analyses of scATAC-seq data. Specifically, SCARP exhibited outstanding cell clustering performance, enabling better elucidation of cell heterogeneity and the discovery of new biologically significant cell subpopulations. Additionally, SCARP was also instrumental in portraying co-accessibility relationships of accessible regions and providing new insight into transcriptional regulation. Consequently, SCARP identified genes that were involved in key Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways related to diseases and predicted reliable cis-regulatory interactions. To sum up, our studies suggested that SCARP is a promising tool to comprehensively analyze the scATAC-seq data.

Crosstalk between chromatin and Shavenbaby defines transcriptional output along the Drosophila intestinal stem cell lineage

The transcription factor Shavenbaby (Svb), the only member of the OvoL family in Drosophila, controls the fate of various epithelial embryonic cells and adult stem cells. Post-translational modification of Svb produces two protein isoforms, Svb-ACT and Svb-REP, which promote adult intestinal stem cell renewal or differentiation, respectively. To define Svb mode of action, we used engineered cell lines and develop an unbiased method to identify Svb target genes across different contexts. Within a given cell type, Svb-ACT and Svb-REP antagonistically regulate the expression of a set of target genes, binding specific enhancers whose accessibility is constrained by chromatin landscape. Reciprocally, Svb-REP can influence local chromatin marks of active enhancers to help repressing target genes. Along the intestinal lineage, the set of Svb target genes progressively changes, together with chromatin accessibility. We propose that Svb-ACT-to-REP transition promotes enterocyte differentiation of intestinal stem cells through direct gene regulation and chromatin remodeling.

Hallmarks of stemness in mammalian tissues

All adult tissues experience wear and tear. Most tissues can compensate for cell loss through the activity of resident stem cells. Although the cellular maintenance strategies vary greatly between different adult (read: postnatal) tissues, the function of stem cells is best defined by their capacity to replace lost tissue through division. We discuss a set of six complementary hallmarks that are key enabling features of this basic function. These include longevity and self-renewal, multipotency, transplantability, plasticity, dependence on niche signals, and maintenance of genome integrity. We discuss these hallmarks in the context of some of the best-understood adult stem cell niches.

Nlrc3 signaling is indispensable for hematopoietic stem cell emergence via Notch signaling in vertebrates

Hematopoietic stem and progenitor cells generate all the lineages of blood cells throughout the lifespan of vertebrates. The emergence of hematopoietic stem and progenitor cells is finely tuned by a variety of signaling pathways. Previous studies have revealed the roles of pattern-recognition receptors such as Toll-like receptors and RIG-I-like receptors in hematopoiesis. In this study, we find that Nlrc3, a nucleotide-binding domain leucine-rich repeat containing family gene, is highly expressed in hematopoietic differentiation stages in vivo and vitro and is required in hematopoiesis in zebrafish. Mechanistically, nlrc3 activates the Notch pathway and the downstream gene of Notch hey1. Furthermore, NF-kB signaling acts upstream of nlrc3 to enhance its transcriptional activity. Finally, we find that Nlrc3 signaling is conserved in the regulation of murine embryonic hematopoiesis. Taken together, our findings uncover an indispensable role of Nlrc3 signaling in hematopoietic stem and progenitor cell emergence and provide insights into inflammation-related hematopoietic ontogeny and the in vitro expansion of hematopoietic stem and progenitor cells.

The evolving views of hematopoiesis: from embryo to adulthood and from in vivo to in vitro

The hematopoietic system composed of hematopoietic stem and progenitor cells (HSPCs) and their differentiated lineages serves as an ideal model to uncover generic principles of cell fate transitions. From gastrulation onwards, there successively emerge primitive hematopoiesis (that produces specialized hematopoietic cells), pro-definitive hematopoiesis (that produces lineage-restricted progenitor cells), and definitive hematopoiesis (that produces multipotent HSPCs). These nascent lineages develop in several transient hematopoietic sites and finally colonize into lifelong hematopoietic sites. The development and maintenance of hematopoietic lineages are orchestrated by cell-intrinsic gene regulatory networks and cell-extrinsic microenvironmental cues. Owing to the progressive methodology (e.g., high-throughput lineage tracing and single-cell functional and omics analyses), our understanding of the developmental origin of hematopoietic lineages and functional properties of certain hematopoietic organs has been updated; meanwhile, new paradigms to characterize rare cell types, cell heterogeneity and its causes, and comprehensive regulatory landscapes have been provided. Here, we review the evolving views of HSPC biology during developmental and postnatal hematopoiesis. Moreover, we discuss recent advances in the in vitro induction and expansion of HSPCs, with a focus on the implications for clinical applications.

Deep generative modeling of transcriptional dynamics for RNA velocity analysis in single cells

RNA velocity has been rapidly adopted to guide interpretation of transcriptional dynamics in snapshot single-cell data; however, current approaches for estimating RNA velocity lack effective strategies for quantifying uncertainty and determining the overall applicability to the system of interest. Here, we present veloVI (velocity variational inference), a deep generative modeling framework for estimating RNA velocity. veloVI learns a gene-specific dynamical model of RNA metabolism and provides a transcriptome-wide quantification of velocity uncertainty. We show that veloVI compares favorably to previous approaches with respect to goodness of fit, consistency across transcriptionally similar cells and stability across preprocessing pipelines for quantifying RNA abundance. Further, we demonstrate that veloVI's posterior velocity uncertainty can be used to assess whether velocity analysis is appropriate for a given dataset. Finally, we highlight veloVI as a flexible framework for modeling transcriptional dynamics by adapting the underlying dynamical model to use time-dependent transcription rates.

A longitudinal single-cell atlas of treatment response in pediatric AML

Pediatric acute myeloid leukemia (pAML) is characterized by heterogeneous cellular composition, driver alterations and prognosis. Characterization of this heterogeneity and how it affects treatment response remains understudied in pediatric patients. We used single-cell RNA sequencing and single-cell ATAC sequencing to profile 28 patients representing different pAML subtypes at diagnosis, remission and relapse. At diagnosis, cellular composition differed between genetic subgroups. Upon relapse, cellular hierarchies transitioned toward a more primitive state regardless of subtype. Primitive cells in the relapsed tumor were distinct compared to cells at diagnosis, with under-representation of myeloid transcriptional programs and over-representation of other lineage programs. In some patients, this was accompanied by the appearance of a B-lymphoid-like hierarchy. Our data thus reveal the emergence of apparent subtype-specific plasticity upon treatment and inform on potentially targetable processes.

SingleScan: a comprehensive resource for single-cell sequencing data processing and mining

Single-cell sequencing has shed light on previously inaccessible biological questions from different fields of research, including organism development, immune function, and disease progression. The number of single-cell-based studies increased dramatically over the past decade. Several new methods and tools have been continuously developed, making it extremely tricky to navigate this research landscape and develop an up-to-date workflow to analyze single-cell sequencing data, particularly for researchers seeking to enter this field without computational experience. Moreover, choosing appropriate tools and optimal parameters to meet the demands of researchers represents a major challenge in processing single-cell sequencing data. However, a specific resource for easy access to detailed information on single-cell sequencing methods and data processing pipelines is still lacking. In the present study, an online resource called SingleScan was developed to curate all up-to-date single-cell transcriptome/genome analyzing tools and pipelines. All the available tools were categorized according to their main tasks, and several typical workflows for single-cell data analysis were summarized. In addition, spatial transcriptomics, which is a breakthrough molecular analysis method that enables researchers to measure all gene activity in tissue samples and map the site of activity, was included along with a portion of single-cell and spatial analysis solutions. For each processing step, the available tools and specific parameters used in published articles are provided and how these parameters affect the results is shown in the resource. All information used in the resource was manually extracted from related literature. An interactive website was designed for data retrieval, visualization, and download. By analyzing the included tools and literature, users can gain insights into the trends of single-cell studies and easily grasp the specific usage of a specific tool. SingleScan will facilitate the analysis of single-cell sequencing data and promote the development of new tools to meet the growing and diverse needs of the research community. The SingleScan database is publicly accessible via the website at http://cailab.labshare.cn/SingleScan .

Cofea: correlation-based feature selection for single-cell chromatin accessibility data

Single-cell chromatin accessibility sequencing (scCAS) technologies have enabled characterizing the epigenomic heterogeneity of individual cells. However, the identification of features of scCAS data that are relevant to underlying biological processes remains a significant gap. Here, we introduce a novel method Cofea, to fill this gap. Through comprehensive experiments on 5 simulated and 54 real datasets, Cofea demonstrates its superiority in capturing cellular heterogeneity and facilitating downstream analysis. Applying this method to identification of cell type-specific peaks and candidate enhancers, as well as pathway enrichment analysis and partitioned heritability analysis, we illustrate the potential of Cofea to uncover functional biological process.

scATAC-seq preprocessing and imputation evaluation system for visualization, clustering and digital footprinting

Single-cell ATAC-seq (scATAC-seq) is a recently developed approach that provides means to investigate open chromatin at single cell level, to assess epigenetic regulation and transcription factors binding landscapes. The sparsity of the scATAC-seq data calls for imputation. Similarly, preprocessing (filtering) may be required to reduce computational load due to the large number of open regions. However, optimal strategies for both imputation and preprocessing have not been yet evaluated together. We present SAPIEnS (scATAC-seq Preprocessing and Imputation Evaluation System), a benchmark for scATAC-seq imputation frameworks, a combination of state-of-the-art imputation methods with commonly used preprocessing techniques. We assess different types of scATAC-seq analysis, i.e. clustering, visualization and digital genomic footprinting, and attain optimal preprocessing-imputation strategies. We discuss the benefits of the imputation framework depending on the task and the number of the dataset features (peaks). We conclude that the preprocessing with the Boruta method is beneficial for the majority of tasks, while imputation is helpful mostly for small datasets. We also implement a SAPIEnS database with pre-computed transcription factor footprints based on imputed data with their activity scores in a specific cell type. SAPIEnS is published at: https://github.com/lab-medvedeva/SAPIEnS. SAPIEnS database is available at: https://sapiensdb.com.