Single-cell multiomic analysis identifies regulatory programs in mixed-phenotype acute leukemia

Mapping the chromatin accessibility landscape of zebrafish embryogenesis at single-cell resolution by SPATAC-seq

Currently, the dynamic accessible elements that determine regulatory programs responsible for the unique identity and function of each cell type during zebrafish embryogenesis lack detailed study. Here we present SPATAC-seq: a split-pool ligation-based assay for transposase-accessible chromatin using sequencing. Using SPATAC-seq, we profiled chromatin accessibility in more than 800,000 individual nuclei across 20 developmental stages spanning the sphere stage to the early larval protruding mouth stage. Using this chromatin accessibility map, we identified 604 cell states and inferred their developmental relationships. We also identified 959,040 candidate cis-regulatory elements (cCREs) and delineated development-specific cCREs, as well as transcription factors defining diverse cell identities. Importantly, enhancer reporter assays confirmed that the majority of tested cCREs exhibited robust enhanced green fluorescent protein expression in restricted cell types or tissues. Finally, we explored gene regulatory programs that drive pigment and notochord cell differentiation. Our work provides a valuable open resource for exploring driver regulators of cell fate decisions in zebrafish embryogenesis.

Defining the Regulatory Logic of Breast Cancer Using Single-Cell Epigenetic and Transcriptome Profiling

Annotation of the cis-regulatory elements that drive transcriptional dysregulation in cancer cells is critical to improving our understanding of tumor biology. Herein, we present a compendium of matched chromatin accessibility (scATAC-seq) and transcriptome (scRNA-seq) profiles at single-cell resolution from human breast tumors and healthy mammary tissues processed immediately following surgical resection. We identify the most likely cell-of-origin for luminal breast tumors and basal breast tumors and then introduce a novel methodology that implements linear mixed-effects models to systematically quantify associations between regions of chromatin accessibility (i.e. regulatory elements) and gene expression in malignant cells versus normal mammary epithelial cells. These data unveil regulatory elements with that switch from silencers of gene expression in normal cells to enhancers of gene expression in cancer cells, leading to the upregulation of clinically relevant oncogenes. To translate the utility of this dataset into tractable models, we generated matched scATAC-seq and scRNA-seq profiles for breast cancer cell lines, revealing, for each subtype, a conserved oncogenic gene expression program between in vitro and in vivo cells. Together, this work highlights the importance of non-coding regulatory mechanisms that underlie oncogenic processes and the ability of single-cell multi-omics to define the regulatory logic of BC cells at single-cell resolution.

pTDP-43 levels correlate with cell type specific molecular alterations in the prefrontal cortex of C9orf72 ALS/FTD patients

Repeat expansions in the C9orf72 gene are the most common genetic cause of amyotrophic lateral sclerosis and familial frontotemporal dementia (ALS/FTD). To identify molecular defects that take place in the dorsolateral frontal cortex of patients with C9orf72 ALS/FTD, we compared healthy controls with C9orf72 ALS/FTD donor samples staged based on the levels of cortical phosphorylated TAR DNA binding protein (pTDP-43), a neuropathological hallmark of disease progression. We identified distinct molecular changes in different cell types that take place during FTD development. Loss of neurosurveillance microglia and activation of the complement cascade take place early, when pTDP-43 aggregates are absent or very low, and become more pronounced in late stages, suggesting an initial involvement of microglia in disease progression. Reduction of layer 2-3 cortical projection neurons with high expression of CUX2/LAMP5 also occurs early, and the reduction becomes more pronounced as pTDP-43 accumulates. Several unique features were observed only in samples with high levels of pTDP-43, including global alteration of chromatin accessibility in oligodendrocytes, microglia, and astrocytes; higher ratios of premature oligodendrocytes; increased levels of the noncoding RNA NEAT1 in astrocytes and neurons, and higher amount of phosphorylated ribosomal protein S6. Our findings reveal previously unknown progressive functional changes in major cell types found in the frontal cortex of C9orf72 ALS/FTD patients that shed light on the mechanisms underlying the pathology of this disease.

CAraCAl: CAMML with the integration of chromatin accessibility

BACKGROUND

A vital step in analyzing single-cell data is ascertaining which cell types are present in a dataset, and at what abundance. In many diseases, the proportions of varying cell types can have important implications for health and prognosis. Most approaches for cell type annotation have centered around cell typing for single-cell RNA-sequencing (scRNA-seq) and have had promising success. However, reliable methods are lacking for many other single-cell modalities such as single-cell sequencing assay for transposase-accessible chromatin (scATAC-seq), which quantifies the extent to which genes of interest in each cell are epigenetically "open" for expression.

RESULTS

To leverage the informative potential of scATAC-seq data, we developed CAMML with the integration of chromatin accessibility (CAraCAl), a bioinformatic method that performs cell typing on scATAC-seq data. CAraCAl performs cell typing by scoring each cell for its enrichment of cell type-specific gene sets. These gene sets are composed of the most upregulated or downregulated genes present in each cell type according to projected gene activity.

CONCLUSIONS

We found that CAraCAl does not improve performance beyond CAMML when scRNA-seq is present, but if only scATAC-seq is available, CAraCAl performs cell typing relatively successfully. As such, we also discuss best practices for cell typing and the strengths and weaknesses of various cell annotation options.

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.

Mapping the cellular biogeography of human bone marrow niches using single-cell transcriptomics and proteomic imaging

Non-hematopoietic cells are essential contributors to hematopoiesis. However, heterogeneity and spatial organization of these cells in human bone marrow remain largely uncharacterized. We used single-cell RNA sequencing (scRNA-seq) to profile 29,325 non-hematopoietic cells and discovered nine transcriptionally distinct subtypes. We simultaneously profiled 53,417 hematopoietic cells and predicted their interactions with non-hematopoietic subsets. We employed co-detection by indexing (CODEX) to spatially profile over 1.2 million cells. We integrated scRNA-seq and CODEX data to link predicted cellular signaling with spatial proximity. Our analysis revealed a hyperoxygenated arterio-endosteal neighborhood for early myelopoiesis, and an adipocytic localization for early hematopoietic stem and progenitor cells (HSPCs). We used our CODEX atlas to annotate new images and uncovered mesenchymal stromal cell (MSC) expansion and spatial neighborhoods co-enriched for leukemic blasts and MSCs in acute myeloid leukemia (AML) patient samples. This spatially resolved, multiomic atlas of human bone marrow provides a reference for investigation of cellular interactions that drive hematopoiesis.

Remodeling ceramide homeostasis promotes functional maturation of human pluripotent stem cell-derived β cells

Human pluripotent stem cell-derived β cells (hPSC-β cells) show the potential to restore euglycemia. However, the immature functionality of hPSC-β cells has limited their efficacy in application. Here, by deciphering the continuous maturation process of hPSC-β cells post transplantation via single-cell RNA sequencing (scRNA-seq) and single-cell assay for transposase-accessible chromatin sequencing (scATAC-seq), we show that functional maturation of hPSC-β cells is an orderly multistep process during which cells sequentially undergo metabolic adaption, removal of negative regulators of cell function, and establishment of a more specialized transcriptome and epigenome. Importantly, remodeling lipid metabolism, especially downregulating the metabolic activity of ceramides, the central hub of sphingolipid metabolism, is critical for β cell maturation. Limiting intracellular accumulation of ceramides in hPSC-β cells remarkably enhanced their function, as indicated by improvements in insulin processing and glucose-stimulated insulin secretion. In summary, our findings provide insights into the maturation of human pancreatic β cells and highlight the importance of ceramide homeostasis in function acquisition.

The non-canonical poly(A) polymerase FAM46C promotes erythropoiesis

The post-transcriptional regulation of mRNA is a crucial component of gene expression. The disruption of this process has detrimental effects on the normal development and gives rise to various diseases. Searching for novel post-transcriptional regulators and exploring their roles are essential for understanding development and disease. Through a multimodal analysis of red blood cell trait genome-wide association studies (GWAS) and transcriptomes of erythropoiesis, we identify FAM46C, a non-canonical RNA poly(A) polymerase, as a necessary factor for proper red blood cell development. FAM46C is highly expressed in the late stages of the erythroid lineage, and its developmental upregulation is controlled by an erythroid-specific enhancer. We demonstrate that FAM46C stabilizes mRNA and regulates erythroid differentiation in a polymerase activity-dependent manner. Furthermore, we identify transcripts of lysosome and mitochondria components as highly confident in vivo targets of FAM46C, which aligns with the need of maturing red blood cells for substantial clearance of organelles and maintenance of cellular redox homeostasis. In conclusion, our study unveils a unique role of FAM46C in positively regulating lysosome and mitochondria components, thereby promoting erythropoiesis.

Convergent epigenetic evolution drives relapse in acute myeloid leukemia

Relapse of acute myeloid leukemia (AML) is highly aggressive and often treatment refractory. We analyzed previously published AML relapse cohorts and found that 40% of relapses occur without changes in driver mutations, suggesting that non-genetic mechanisms drive relapse in a large proportion of cases. We therefore characterized epigenetic patterns of AML relapse using 26 matched diagnosis-relapse samples with ATAC-seq. This analysis identified a relapse-specific chromatin accessibility signature for mutationally stable AML, suggesting that AML undergoes epigenetic evolution at relapse independent of mutational changes. Analysis of leukemia stem cell (LSC) chromatin changes at relapse indicated that this leukemic compartment underwent significantly less epigenetic evolution than non-LSCs, while epigenetic changes in non-LSCs reflected overall evolution of the bulk leukemia. Finally, we used single-cell ATAC-seq paired with mitochondrial sequencing (mtscATAC) to map clones from diagnosis into relapse along with their epigenetic features. We found that distinct mitochondrially-defined clones exhibit more similar chromatin accessibility at relapse relative to diagnosis, demonstrating convergent epigenetic evolution in relapsed AML. These results demonstrate that epigenetic evolution is a feature of relapsed AML and that convergent epigenetic evolution can occur following treatment with induction chemotherapy.

Genomic and phenotypic correlates of mosaic loss of chromosome Y in blood

Mosaic loss of Y (mLOY) is the most common somatic chromosomal alteration detected in human blood. The presence of mLOY is associated with altered blood cell counts and increased risk of Alzheimer's disease, solid tumors, and other age-related diseases. We sought to gain a better understanding of genetic drivers and associated phenotypes of mLOY through analyses of whole genome sequencing of a large set of genetically diverse males from the Trans-Omics for Precision Medicine (TOPMed) program. This approach enabled us to identify differences in mLOY frequencies across populations defined by genetic similarity, revealing a higher frequency of mLOY in the European American (EA) ancestry group compared to those of Hispanic American (HA), African American (AA), and East Asian (EAS) ancestry. Further, we identified two genes ( CFHR1 and LRP6 ) that harbor multiple rare, putatively deleterious variants associated with mLOY susceptibility, show that subsets of human hematopoietic stem cells are enriched for activity of mLOY susceptibility variants, and that certain alleles on chromosome Y are more likely to be lost than others.

TIM-3+ CD8 T cells with a terminally exhausted phenotype retain functional capacity in hematological malignancies

Chronic antigen stimulation is thought to generate dysfunctional CD8 T cells. Here, we identify a CD8 T cell subset in the bone marrow tumor microenvironment that, despite an apparent terminally exhausted phenotype (TPHEX), expressed granzymes, perforin, and IFN-γ. Concurrent gene expression and DNA accessibility revealed that genes encoding these functional proteins correlated with BATF expression and motif accessibility. IFN-γ+ TPHEX effectively killed myeloma with comparable efficacy to transitory effectors, and disease progression correlated with numerical deficits in IFN-γ+ TPHEX. We also observed IFN-γ+ TPHEX within CD19-targeted chimeric antigen receptor T cells, which killed CD19+ leukemia cells. An IFN-γ+ TPHEX gene signature was recapitulated in TEX cells from human cancers, including myeloma and lymphoma. Here, we characterize a TEX subset in hematological malignancies that paradoxically retains function and is distinct from dysfunctional TEX found in chronic viral infections. Thus, IFN-γ+ TPHEX represent a potential target for immunotherapy of blood cancers.

A noncoding regulatory variant in IKZF1 increases acute lymphoblastic leukemia risk in Hispanic/Latino children

Hispanic/Latino children have the highest risk of acute lymphoblastic leukemia (ALL) in the US compared to other racial/ethnic groups, yet the basis of this remains incompletely understood. Through genetic fine-mapping analyses, we identified a new independent childhood ALL risk signal near IKZF1 in self-reported Hispanic/Latino individuals, but not in non-Hispanic White individuals, with an effect size of ∼1.44 (95% confidence interval = 1.33-1.55) and a risk allele frequency of ∼18% in Hispanic/Latino populations and <0.5% in European populations. This risk allele was positively associated with Indigenous American ancestry, showed evidence of selection in human history, and was associated with reduced IKZF1 expression. We identified a putative causal variant in a downstream enhancer that is most active in pro-B cells and interacts with the IKZF1 promoter. This variant disrupts IKZF1 autoregulation at this enhancer and results in reduced enhancer activity in B cell progenitors. Our study reveals a genetic basis for the increased ALL risk in Hispanic/Latino children.

Single-cell Landscape of Malignant Transition: Unraveling Cancer Cell-of-Origin and Heterogeneous Tissue Microenvironment

Understanding disease progression and sophisticated tumor ecosystems is imperative for investigating tumorigenesis mechanisms and developing novel prevention strategies. Here, we dissected heterogeneous microenvironments during malignant transitions by leveraging data from 1396 samples spanning 13 major tissues. Within transitional stem-like subpopulations highly enriched in precancers and cancers, we identified 30 recurring cellular states strongly linked to malignancy, including hypoxia and epithelial senescence, revealing a high degree of plasticity in epithelial stem cells. By characterizing dynamics in stem-cell crosstalk with the microenvironment along the pseudotime axis, we found differential roles of ANXA1 at different stages of tumor development. In precancerous stages, reduced ANXA1 levels promoted monocyte differentiation toward M1 macrophages and inflammatory responses, whereas during malignant progression, upregulated ANXA1 fostered M2 macrophage polarization and cancer-associated fibroblast transformation by increasing TGF-β production. Our spatiotemporal analysis further provided insights into mechanisms responsible for immunosuppression and a potential target to control evolution of precancer and mitigate the risk for cancer development.

Tissue-specific enhancer-gene maps from multimodal single-cell data identify causal disease alleles

Translating genome-wide association study (GWAS) loci into causal variants and genes requires accurate cell-type-specific enhancer-gene maps from disease-relevant tissues. Building enhancer-gene maps is essential but challenging with current experimental methods in primary human tissues. Here we developed a nonparametric statistical method, SCENT (single-cell enhancer target gene mapping), that models association between enhancer chromatin accessibility and gene expression in single-cell or nucleus multimodal RNA sequencing and ATAC sequencing data. We applied SCENT to 9 multimodal datasets including >120,000 single cells or nuclei and created 23 cell-type-specific enhancer-gene maps. These maps were highly enriched for causal variants in expression quantitative loci and GWAS for 1,143 diseases and traits. We identified likely causal genes for both common and rare diseases and linked somatic mutation hotspots to target genes. We demonstrate that application of SCENT to multimodal data from disease-relevant human tissue enables the scalable construction of accurate cell-type-specific enhancer-gene maps, essential for defining noncoding variant function.

An immunophenotype-coupled transcriptomic atlas of human hematopoietic progenitors

Analysis of the human hematopoietic progenitor compartment is being transformed by single-cell multimodal approaches. Cellular indexing of transcriptomes and epitopes by sequencing (CITE-seq) enables coupled surface protein and transcriptome profiling, thereby revealing genomic programs underlying progenitor states. To perform CITE-seq systematically on primary human bone marrow cells, we used titrations with 266 CITE-seq antibodies (antibody-derived tags) and machine learning to optimize a panel of 132 antibodies. Multimodal analysis resolved >80 stem, progenitor, immune, stromal and transitional cells defined by distinctive surface markers and transcriptomes. This dataset enables flow cytometry solutions for in silico-predicted cell states and identifies dozens of cell surface markers consistently detected across donors spanning race and sex. Finally, aligning annotations from this atlas, we nominate normal marrow equivalents for acute myeloid leukemia stem cell populations that differ in clinical response. This atlas serves as an advanced digital resource for hematopoietic progenitor analyses in human health and disease.

Beyond single cells: microfluidics empowering multiomics analysis

Single-cell multiomics technologies empower simultaneous measurement of multiple types of molecules within individual cells, providing a more profound comprehension compared with the analysis of discrete molecular layers from different cells. Microfluidic technology, on the other hand, has emerged as a pivotal facilitator for high-throughput single-cell analysis, offering precise control and manipulation of individual cells. The primary focus of this review encompasses an appraisal of cutting-edge microfluidic platforms employed in the realm of single-cell multiomics analysis. Furthermore, it discusses technological advancements in various single-cell omics such as genomics, transcriptomics, epigenomics, and proteomics, with their perspective applications. Finally, it provides future prospects of these integrated single-cell multiomics methodologies, shedding light on the possibilities for future biological research.

HyGAnno: hybrid graph neural network-based cell type annotation for single-cell ATAC sequencing data

Reliable cell type annotations are crucial for investigating cellular heterogeneity in single-cell omics data. Although various computational approaches have been proposed for single-cell RNA sequencing (scRNA-seq) annotation, high-quality cell labels are still lacking in single-cell sequencing assay for transposase-accessible chromatin (scATAC-seq) data, because of extreme sparsity and inconsistent chromatin accessibility between datasets. Here, we present a novel automated cell annotation method that transfers cell type information from a well-labeled scRNA-seq reference to an unlabeled scATAC-seq target, via a parallel graph neural network, in a semi-supervised manner. Unlike existing methods that utilize only gene expression or gene activity features, HyGAnno leverages genome-wide accessibility peak features to facilitate the training process. In addition, HyGAnno reconstructs a reference-target cell graph to detect cells with low prediction reliability, according to their specific graph connectivity patterns. HyGAnno was assessed across various datasets, showcasing its strengths in precise cell annotation, generating interpretable cell embeddings, robustness to noisy reference data and adaptability to tumor tissues.

Mapping the Cellular Biogeography of Human Bone Marrow Niches Using Single-Cell Transcriptomics and Proteomic Imaging

The bone marrow is the organ responsible for blood production. Diverse non-hematopoietic cells contribute essentially to hematopoiesis. However, these cells and their spatial organization remain largely uncharacterized as they have been technically challenging to study in humans. Here, we used fresh femoral head samples and performed single-cell RNA sequencing (scRNA-Seq) to profile 29,325 enriched non-hematopoietic bone marrow cells and discover nine transcriptionally distinct subtypes. We next employed CO-detection by inDEXing (CODEX) multiplexed imaging of 18 individuals, including both healthy and acute myeloid leukemia (AML) samples, to spatially profile over one million single cells with a novel 53-antibody panel. We discovered a relatively hyperoxygenated arterio-endosteal niche for early myelopoiesis, and an adipocytic, but not endosteal or perivascular, niche for early hematopoietic stem and progenitor cells. We used our atlas to predict cell type labels in new bone marrow images and used these predictions to uncover mesenchymal stromal cell (MSC) expansion and leukemic blast/MSC-enriched spatial neighborhoods in AML patient samples. Our work represents the first comprehensive, spatially-resolved multiomic atlas of human bone marrow and will serve as a reference for future investigation of cellular interactions that drive hematopoiesis.

Disruption of mitochondrial energy metabolism is a putative pathogenesis of Diamond-Blackfan anemia

Energy metabolism in the context of erythropoiesis and related diseases remains largely unexplored. Here, we developed a primary cell model by differentiating hematopoietic stem progenitor cells toward the erythroid lineage and suppressing the mitochondrial oxidative phosphorylation (OXPHOS) pathway. OXPHOS suppression led to differentiation failure of erythroid progenitors and defects in ribosome biogenesis. Ran GTPase-activating protein 1 (RanGAP1) was identified as a target of mitochondrial OXPHOS for ribosomal defects during erythropoiesis. Overexpression of RanGAP1 largely alleviated erythroid defects resulting from OXPHOS suppression. Coenzyme Q10, an activator of OXPHOS, largely rescued erythroid defects and increased RanGAP1 expression. Patients with Diamond-Blackfan anemia (DBA) exhibited OXPHOS suppression and a concomitant suppression of ribosome biogenesis. RNA-seq analysis implied that the substantial mutation (approximately 10%) in OXPHOS genes accounts for OXPHOS suppression in these patients. Conclusively, OXPHOS disruption and the associated disruptive mitochondrial energy metabolism are linked to the pathogenesis of DBA.

Temporally resolved early BMP-driven transcriptional cascade during human amnion specification

Amniogenesis, a process critical for continuation of healthy pregnancy, is triggered in a collection of pluripotent epiblast cells as the human embryo implants. Previous studies have established that BMP signaling is a major driver of this lineage specifying process, but the downstream BMP-dependent transcriptional networks that lead to successful amniogenesis remain to be identified. This is, in part, due to the current lack of a robust and reproducible model system that enables mechanistic investigations exclusively into amniogenesis. Here, we developed an improved model of early amnion specification, using a human pluripotent stem cell-based platform in which the activation of BMP signaling is controlled and synchronous. Uniform amniogenesis is seen within 48 hours after BMP activation, and the resulting cells share transcriptomic characteristics with amnion cells of a gastrulating human embryo. Using detailed time-course transcriptomic analyses, we established a previously uncharacterized BMP-dependent amniotic transcriptional cascade, and identified markers that represent five distinct stages of amnion fate specification; the expression of selected markers was validated in early post-implantation macaque embryos. Moreover, a cohort of factors that could potentially control specific stages of amniogenesis was identified, including the transcription factor TFAP2A. Functionally, we determined that, once amniogenesis is triggered by the BMP pathway, TFAP2A controls the progression of amniogenesis. This work presents a temporally resolved transcriptomic resource for several previously uncharacterized amniogenesis states and demonstrates a critical intermediate role for TFAP2A during amnion fate specification.

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.

Bacterial phenotypic heterogeneity through the lens of single-cell RNA sequencing

Bacterial transcription is not monolithic. Microbes exist in a wide variety of cell states that help them adapt to their environment, acquire and produce essential nutrients, and engage in both competition and cooperation with their neighbors. While we typically think of bacterial adaptation as a group behavior, where all cells respond in unison, there is often a mixture of phenotypic responses within a bacterial population, where distinct cell types arise. A primary phenomenon driving these distinct cell states is transcriptional heterogeneity. Given that bacterial mRNA transcripts are extremely short-lived compared to eukaryotes, their transcriptional state is closely associated with their physiology, and thus the transcriptome of a bacterial cell acts as a snapshot of the behavior of that bacterium. Therefore, the application of single-cell transcriptomics to microbial populations will provide novel insight into cellular differentiation and bacterial ecology. In this review, we provide an overview of transcriptional heterogeneity in microbial systems, discuss the findings already provided by single-cell approaches, and plot new avenues of inquiry in transcriptional regulation, cellular biology, and mechanisms of heterogeneity that are made possible when microbial communities are analyzed at single-cell resolution.

Dictionary learning for integrative, multimodal and scalable single-cell analysis

Mapping single-cell sequencing profiles to comprehensive reference datasets provides a powerful alternative to unsupervised analysis. However, most reference datasets are constructed from single-cell RNA-sequencing data and cannot be used to annotate datasets that do not measure gene expression. Here we introduce 'bridge integration', a method to integrate single-cell datasets across modalities using a multiomic dataset as a molecular bridge. Each cell in the multiomic dataset constitutes an element in a 'dictionary', which is used to reconstruct unimodal datasets and transform them into a shared space. Our procedure accurately integrates transcriptomic data with independent single-cell measurements of chromatin accessibility, histone modifications, DNA methylation and protein levels. Moreover, we demonstrate how dictionary learning can be combined with sketching techniques to improve computational scalability and harmonize 8.6 million human immune cell profiles from sequencing and mass cytometry experiments. Our approach, implemented in version 5 of our Seurat toolkit ( http://www.satijalab.org/seurat ), broadens the utility of single-cell reference datasets and facilitates comparisons across diverse molecular modalities.

Single-nucleus multi-omic profiling of human placental syncytiotrophoblasts identifies cellular trajectories during pregnancy

The human placenta has a vital role in ensuring a successful pregnancy. Despite the growing body of knowledge about its cellular compositions and functions, there has been limited research on the heterogeneity of the billions of nuclei within the syncytiotrophoblast (STB), a multinucleated entity primarily responsible for placental function. Here we conducted integrated single-nucleus RNA sequencing and single-nucleus ATAC sequencing analyses of human placentas from early and late pregnancy. Our findings demonstrate the dynamic heterogeneity and developmental trajectories of STB nuclei and their correspondence with human trophoblast stem cell (hTSC)-derived STB. Furthermore, we identified transcription factors associated with diverse STB nuclear lineages through their gene regulatory networks and experimentally confirmed their function in hTSC and trophoblast organoid-derived STBs. Together, our data provide insights into the heterogeneity of human STB and represent a valuable resource for interpreting associated pregnancy complications.

Decoding leukemia at the single-cell level: clonal architecture, classification, microenvironment, and drug resistance

Leukemias are refractory hematological malignancies, characterized by marked intrinsic heterogeneity which poses significant obstacles to effective treatment. However, traditional bulk sequencing techniques have not been able to effectively unravel the heterogeneity among individual tumor cells. With the emergence of single-cell sequencing technology, it has bestowed upon us an unprecedented resolution to comprehend the mechanisms underlying leukemogenesis and drug resistance across various levels, including the genome, epigenome, transcriptome and proteome. Here, we provide an overview of the currently prevalent single-cell sequencing technologies and a detailed summary of single-cell studies conducted on leukemia, with a specific focus on four key aspects: (1) leukemia's clonal architecture, (2) frameworks to determine leukemia subtypes, (3) tumor microenvironment (TME) and (4) the drug-resistant mechanisms of leukemia. This review provides a comprehensive summary of current single-cell studies on leukemia and highlights the markers and mechanisms that show promising clinical implications for the diagnosis and treatment of leukemia.

Scalable single-cell profiling of chromatin modifications with sciCUT&Tag

Cleavage under targets and tagmentation (CUT&Tag) is an antibody-directed in situ chromatin profiling strategy that is rapidly replacing immune precipitation-based methods, such as chromatin immunoprecipitation-sequencing. The efficiency of the method enables chromatin profiling in single cells but is limited by the numbers of cells that can be profiled. Here, we describe a combinatorial barcoding strategy for CUT&Tag that harnesses a nanowell dispenser for simple, high-resolution, high-throughput, single-cell chromatin profiling. In this single-cell combinatorial indexing CUT&Tag (sciCUT&Tag) protocol, lightly cross-linked nuclei are bound to magnetic beads and incubated with primary and secondary antibodies in bulk and then arrayed in a 96-well plate for a first round of cellular indexing by antibody-directed Tn5 tagmentation. The sample is then repooled, mixed and arrayed across 5,184 nanowells at a density of 12-24 nuclei per well for a second round of cellular indexing during PCR amplification of the sequencing-ready library. This protocol can be completed in 1.5 days by a research technician, and we illustrate the optimized protocol by profiling histone modifications associated with developmental gene repression (H3K27me3) as well as transcriptional activation (H3K4me1-2-3) in human peripheral blood mononuclear cells and use single-nucleotide polymorphisms to facilitate collision removal. We have also used sciCUT&Tag for simultaneous profiling of multiple chromatin epitopes in single cells. The reduced cost, improved resolution and scalability of sciCUT&Tag make it an attractive platform to profile chromatin features in single cells.

Chromatin structure and 3D architecture define differential functions of PU.1 cis regulatory elements in human blood cell lineages

The precise spatio-temporal expression of the hematopoietic ETS transcription factor PU.1 that determines the hematopoietic cell fates is tightly regulated at the chromatin level. However, it remains elusive as to how chromatin signatures are linked to this dynamic expression pattern of PU.1 across blood cell lineages. Here we performed an unbiased and in-depth analysis of the relationship between human PU.1 expression, the presence of trans-acting factors, and 3D architecture at various cis-regulatory elements (CRE) proximal to the PU.1 locus. We identified multiple novel CREs at the upstream region of the gene following an integrative inspection for conserved DNA elements at the chromatin-accessible regions in primary human blood lineages. We showed that a subset of CREs localize within a 10 kb-wide cluster that exhibits that exhibit molecular features of a myeloid-specific super-enhancer involved in mediating PU.1 autoregulation, including open chromatin, unmethylated DNA, histone enhancer marks, transcription of enhancer RNAs, and occupancy of the PU.1 protein itself. Importantly, we revealed the presence of common 35-kb-wide CTCF-bound insulated neighborhood that contains the CRE cluster, forming the chromatin territory for lineage-specific and CRE-mediated chromatin interactions. These include functional CRE-promoter interactions in myeloid and B cells but not in erythroid and T cells. Our findings also provide mechanistic insights into the interplay between dynamic chromatin structure and 3D architecture in defining certain CREs as enhancers or silencers in chromatin regulation of PU.1 expression. The study lays the groundwork for further examination of PU.1 CREs as well as epigenetic regulation in malignant hematopoiesis.

A cell type-aware framework for nominating non-coding variants in Mendelian regulatory disorders

Unsolved Mendelian cases often lack obvious pathogenic coding variants, suggesting potential non-coding etiologies. Here, we present a single cell multi-omic framework integrating embryonic mouse chromatin accessibility, histone modification, and gene expression assays to discover cranial motor neuron (cMN) cis-regulatory elements and subsequently nominate candidate non-coding variants in the congenital cranial dysinnervation disorders (CCDDs), a set of Mendelian disorders altering cMN development. We generated single cell epigenomic profiles for ~86,000 cMNs and related cell types, identifying ~250,000 accessible regulatory elements with cognate gene predictions for ~145,000 putative enhancers. Seventy-five percent of elements (44 of 59) validated in an in vivo transgenic reporter assay, demonstrating that single cell accessibility is a strong predictor of enhancer activity. Applying our cMN atlas to 899 whole genome sequences from 270 genetically unsolved CCDD pedigrees, we achieved significant reduction in our variant search space and nominated candidate variants predicted to regulate known CCDD disease genes MAFB, PHOX2A, CHN1, and EBF3 - as well as new candidates in recurrently mutated enhancers through peak- and gene-centric allelic aggregation. This work provides novel non-coding variant discoveries of relevance to CCDDs and a generalizable framework for nominating non-coding variants of potentially high functional impact in other Mendelian disorders.

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.

Therapeutic targeting of EP300/CBP by bromodomain inhibition in hematologic malignancies

CCS1477 (inobrodib) is a potent, selective EP300/CBP bromodomain inhibitor which induces cell-cycle arrest and differentiation in hematologic malignancy model systems. In myeloid leukemia cells, it promotes rapid eviction of EP300/CBP from an enhancer subset marked by strong MYB occupancy and high H3K27 acetylation, with downregulation of the subordinate oncogenic network and redistribution to sites close to differentiation genes. In myeloma cells, CCS1477 induces eviction of EP300/CBP from FGFR3, the target of the common (4; 14) translocation, with redistribution away from IRF4-occupied sites to TCF3/E2A-occupied sites. In a subset of patients with relapsed or refractory disease, CCS1477 monotherapy induces differentiation responses in AML and objective responses in heavily pre-treated multiple myeloma. In vivo preclinical combination studies reveal synergistic responses to treatment with standard-of-care agents. Thus, CCS1477 exhibits encouraging preclinical and early-phase clinical activity by disrupting recruitment of EP300/CBP to enhancer networks occupied by critical transcription factors.

Single cell multi-omics reveal intra-cell-line heterogeneity across human cancer cell lines

Human cancer cell lines have long served as tools for cancer research and drug discovery, but the presence and the source of intra-cell-line heterogeneity remain elusive. Here, we perform single-cell RNA-sequencing and ATAC-sequencing on 42 and 39 human cell lines, respectively, to illustrate both transcriptomic and epigenetic heterogeneity within individual cell lines. Our data reveal that transcriptomic heterogeneity is frequently observed in cancer cell lines of different tissue origins, often driven by multiple common transcriptional programs. Copy number variation, as well as epigenetic variation and extrachromosomal DNA distribution all contribute to the detected intra-cell-line heterogeneity. Using hypoxia treatment as an example, we demonstrate that transcriptomic heterogeneity could be reshaped by environmental stress. Overall, our study performs single-cell multi-omics of commonly used human cancer cell lines and offers mechanistic insights into the intra-cell-line heterogeneity and its dynamics, which would serve as an important resource for future cancer cell line-based studies.

SEACells infers transcriptional and epigenomic cellular states from single-cell genomics data

Metacells are cell groupings derived from single-cell sequencing data that represent highly granular, distinct cell states. Here we present single-cell aggregation of cell states (SEACells), an algorithm for identifying metacells that overcome the sparsity of single-cell data while retaining heterogeneity obscured by traditional cell clustering. SEACells outperforms existing algorithms in identifying comprehensive, compact and well-separated metacells in both RNA and assay for transposase-accessible chromatin (ATAC) modalities across datasets with discrete cell types and continuous trajectories. We demonstrate the use of SEACells to improve gene-peak associations, compute ATAC gene scores and infer the activities of critical regulators during differentiation. Metacell-level analysis scales to large datasets and is particularly well suited for patient cohorts, where per-patient aggregation provides more robust units for data integration. We use our metacells to reveal expression dynamics and gradual reconfiguration of the chromatin landscape during hematopoietic differentiation and to uniquely identify CD4 T cell differentiation and activation states associated with disease onset and severity in a Coronavirus Disease 2019 (COVID-19) patient cohort.

Resolving therapy resistance mechanisms in multiple myeloma by multiomics subclone analysis

Intratumor heterogeneity as a clinical challenge becomes most evident after several treatment lines, when multidrug-resistant subclones accumulate. To address this challenge, the characterization of resistance mechanisms at the subclonal level is key to identify common vulnerabilities. In this study, we integrate whole-genome sequencing, single-cell (sc) transcriptomics (scRNA sequencing), and chromatin accessibility (scATAC sequencing) together with mitochondrial DNA mutations to define subclonal architecture and evolution for longitudinal samples from 15 patients with relapsed or refractory multiple myeloma. We assess transcriptomic and epigenomic changes to resolve the multifactorial nature of therapy resistance and relate it to the parallel occurrence of different mechanisms: (1) preexisting epigenetic profiles of subclones associated with survival advantages, (2) converging phenotypic adaptation of genetically distinct subclones, and (3) subclone-specific interactions of myeloma and bone marrow microenvironment cells. Our study showcases how an integrative multiomics analysis can be applied to track and characterize distinct multidrug-resistant subclones over time for the identification of molecular targets against them.

Microfluidics-free single-cell genomics with templated emulsification

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

Identification and targeting of treatment resistant progenitor populations in T-cell Acute Lymphoblastic Leukemia

Refractoriness to initial chemotherapy and relapse after remission are the main obstacles to cure in T-cell Acute Lymphoblastic Leukemia (T-ALL). Biomarker guided risk stratification and targeted therapy have the potential to improve outcomes in high-risk T-ALL; however, cellular and genetic factors contributing to treatment resistance remain unknown. Previous bulk genomic studies in T-ALL have implicated tumor heterogeneity as an unexplored mechanism for treatment failure. To link tumor subpopulations with clinical outcome, we created an atlas of healthy pediatric hematopoiesis and applied single-cell multiomic (CITE-seq/snATAC-seq) analysis to a cohort of 40 cases of T-ALL treated on the Children's Oncology Group AALL0434 clinical trial. The cohort was carefully selected to capture the immunophenotypic diversity of T-ALL, with early T-cell precursor (ETP) and Near/Non-ETP subtypes represented, as well as enriched with both relapsed and treatment refractory cases. Integrated analyses of T-ALL blasts and normal T-cell precursors identified a bone-marrow progenitor-like (BMP-like) leukemia sub-population associated with treatment failure and poor overall survival. The single-cell-derived molecular signature of BMP-like blasts predicted poor outcome across multiple subtypes of T-ALL within two independent patient cohorts using bulk RNA-sequencing data from over 1300 patients. We defined the mutational landscape of BMP-like T-ALL, finding that NOTCH1 mutations additively drive T-ALL blasts away from the BMP-like state. We transcriptionally matched BMP-like blasts to early thymic seeding progenitors that have low NR3C1 expression and high stem cell gene expression, corresponding to a corticosteroid and conventional cytotoxic resistant phenotype we observed in ex vivo drug screening. To identify novel targets for BMP-like blasts, we performed in silico and in vitro drug screening against the BMP-like signature and prioritized BMP-like overexpressed cell-surface (CD44, ITGA4, LGALS1) and intracellular proteins (BCL-2, MCL-1, BTK, NF-κB) as candidates for precision targeted therapy. We established patient derived xenograft models of BMP-high and BMP-low leukemias, which revealed vulnerability of BMP-like blasts to apoptosis-inducing agents, TEC-kinase inhibitors, and proteasome inhibitors. Our study establishes the first multi-omic signatures for rapid risk-stratification and targeted treatment of high-risk T-ALL.

XMAP: Cross-population fine-mapping by leveraging genetic diversity and accounting for confounding bias

Fine-mapping prioritizes risk variants identified by genome-wide association studies (GWASs), serving as a critical step to uncover biological mechanisms underlying complex traits. However, several major challenges still remain for existing fine-mapping methods. First, the strong linkage disequilibrium among variants can limit the statistical power and resolution of fine-mapping. Second, it is computationally expensive to simultaneously search for multiple causal variants. Third, the confounding bias hidden in GWAS summary statistics can produce spurious signals. To address these challenges, we develop a statistical method for cross-population fine-mapping (XMAP) by leveraging genetic diversity and accounting for confounding bias. By using cross-population GWAS summary statistics from global biobanks and genomic consortia, we show that XMAP can achieve greater statistical power, better control of false positive rate, and substantially higher computational efficiency for identifying multiple causal signals, compared to existing methods. Importantly, we show that the output of XMAP can be integrated with single-cell datasets, which greatly improves the interpretation of putative causal variants in their cellular context at single-cell resolution.

Uptake of long-chain fatty acids from the bone marrow suppresses CD8+ T-cell metabolism and function in multiple myeloma

T cells demonstrate impaired function in multiple myeloma (MM) but suppressive mechanisms in the bone marrow microenvironment remain poorly defined. We observe that bone marrow CD8+ T-cell function is decreased in MM compared with controls, and is also consistently lower within bone marrow samples than in matched peripheral blood samples. These changes are accompanied by decreased mitochondrial mass and markedly elevated long-chain fatty acid uptake. In vitro modeling confirmed that uptake of bone marrow lipids suppresses CD8+ T function, which is impaired in autologous bone marrow plasma but rescued by lipid removal. Analysis of single-cell RNA-sequencing data identified expression of fatty acid transport protein 1 (FATP1) in bone marrow CD8+ T cells in MM, and FATP1 blockade also rescued CD8+ T-cell function, thereby identifying this as a novel target to augment T-cell activity in MM. Finally, analysis of samples from cohorts of patients who had received treatment identified that CD8+ T-cell metabolic dysfunction resolves in patients with MM who are responsive to treatment but not in patients with relapsed MM, and is associated with substantial T-cell functional restoration.

Single-cell chromatin accessibility and transcriptomic characterization of Behcet's disease

Behect's disease is a chronic vasculitis characterized by complex multi-organ immune aberrations. However, a comprehensive understanding of the gene-regulatory profile of peripheral autoimmunity and the diverse immune responses across distinct cell types in Behcet's disease (BD) is still lacking. Here, we present a multi-omic single-cell study of 424,817 cells in BD patients and non-BD individuals. This study maps chromatin accessibility and gene expression in the same biological samples, unraveling vast cellular heterogeneity. We identify widespread cell-type-specific, disease-associated active and pro-inflammatory immunity in both transcript and epigenomic aspects. Notably, integrative multi-omic analysis reveals putative TF regulators that might contribute to chromatin accessibility and gene expression in BD. Moreover, we predicted gene-regulatory networks within nominated TF activators, including AP-1, NF-kB, and ETS transcript factor families, which may regulate cellular interaction and govern inflammation. Our study illustrates the epigenetic and transcriptional landscape in BD peripheral blood and expands understanding of potential epigenomic immunopathology in this disease.

Using single-cell chromatin accessibility sequencing to characterize CD4+ T cells from murine tissues

The Assay for Transposase-Accessible Chromatin using sequencing (ATAC-seq) is a cutting-edge technology that enables researchers to assess genome-wide chromatin accessibility and to characterize cell type specific gene-regulatory programs. Recent technological progress allows for using this technology also on the single-cell level. In this article, we describe the whole value chain from the isolation of T cells from murine tissues to a complete bioinformatic analysis workflow. We start with methods for isolating scATAC-seq-ready CD4+ T cells from murine tissues such as visceral adipose tissue, skin, colon, and secondary lymphoid tissues such as the spleen. We describe the preparation of nuclei and quality control parameters during library preparation. Based on publicly available sequencing data that was generated using these protocols, we describe a step-by-step bioinformatic analysis pipeline for data pre-processing and downstream analysis. Our analysis workflow will follow the R-based bioinformatics framework ArchR, which is currently well established for scATAC-seq datasets. All in all, this work serves as a one-stop shop for generating and analyzing chromatin accessibility landscapes in T cells.

Single-cell RNA sequencing distinctly characterizes the wide heterogeneity in pediatric mixed phenotype acute leukemia

BACKGROUND

Mixed phenotype acute leukemia (MPAL), a rare subgroup of leukemia characterized by blast cells with myeloid and lymphoid lineage features, is difficult to diagnose and treat. A better characterization of MPAL is essential to understand the subtype heterogeneity and how it compares with acute myeloid leukemia (AML) and acute lymphoblastic leukemia (ALL). Therefore, we performed single-cell RNA sequencing (scRNAseq) on pediatric MPAL bone marrow (BM) samples to develop a granular map of the MPAL blasts and microenvironment landscape.

METHODS

We analyzed over 40,000 cells from nine pediatric MPAL BM samples to generate a single-cell transcriptomic landscape of B/myeloid (B/My) and T/myeloid (T/My) MPAL. Cells were clustered using unsupervised single-cell methods, and malignant blast and immune clusters were annotated. Differential expression analysis was performed to identify B/My and T/My MPAL blast-specific signatures by comparing transcriptome profiles of MPAL with normal BM, AML, and ALL. Gene set enrichment analysis (GSEA) was performed, and significantly enriched pathways were compared in MPAL subtypes.

RESULTS

B/My and T/My MPAL blasts displayed distinct blast signatures. Transcriptomic analysis revealed that B/My MPAL profile overlaps with B-ALL and AML samples. Similarly, T/My MPAL exhibited overlap with T-ALL and AML samples. Genes overexpressed in both MPAL subtypes' blast cells compared to AML, ALL, and healthy BM included MAP2K2 and CD81. Subtype-specific genes included HBEGF for B/My and PTEN for T/My. These marker sets segregated bulk RNA-seq AML, ALL, and MPAL samples based on expression profiles. Analysis comparing T/My MPAL to ETP, near-ETP, and non-ETP T-ALL, showed that T/My MPAL had greater overlap with ETP-ALL cases. Comparisons among MPAL subtypes between adult and pediatric samples showed analogous transcriptomic landscapes of corresponding subtypes. Transcriptomic differences were observed in the MPAL samples based on response to induction chemotherapy, including selective upregulation of the IL-16 pathway in relapsed samples.

CONCLUSIONS

We have for the first time described the single-cell transcriptomic landscape of pediatric MPAL and demonstrated that B/My and T/My MPAL have distinct scRNAseq profiles from each other, AML, and ALL. Differences in transcriptomic profiles were seen based on response to therapy, but larger studies will be needed to validate these findings.

Convergent Epigenetic Evolution Drives Relapse in Acute Myeloid Leukemia

Relapse of acute myeloid leukemia (AML) is highly aggressive and often treatment refractory. We analyzed previously published AML relapse cohorts and found that 40% of relapses occur without changes in driver mutations, suggesting that non-genetic mechanisms drive relapse in a large proportion of cases. We therefore characterized epigenetic patterns of AML relapse using 26 matched diagnosis-relapse samples with ATAC-seq. This analysis identified a relapse-specific chromatin accessibility signature for mutationally stable AML, suggesting that AML undergoes epigenetic evolution at relapse independent of mutational changes. Analysis of leukemia stem cell (LSC) chromatin changes at relapse indicated that this leukemic compartment underwent significantly less epigenetic evolution than non-LSCs, while epigenetic changes in non-LSCs reflected overall evolution of the bulk leukemia. Finally, we used single-cell ATAC-seq paired with mitochondrial sequencing (mtscATAC) to map clones from diagnosis into relapse along with their epigenetic features. We found that distinct mitochondrially-defined clones exhibit more similar chromatin accessibility at relapse relative to diagnosis, demonstrating convergent epigenetic evolution in relapsed AML. These results demonstrate that epigenetic evolution is a feature of relapsed AML and that convergent epigenetic evolution can occur following treatment with induction chemotherapy.

The technological landscape and applications of single-cell multi-omics

Single-cell multi-omics technologies and methods characterize cell states and activities by simultaneously integrating various single-modality omics methods that profile the transcriptome, genome, epigenome, epitranscriptome, proteome, metabolome and other (emerging) omics. Collectively, these methods are revolutionizing molecular cell biology research. In this comprehensive Review, we discuss established multi-omics technologies as well as cutting-edge and state-of-the-art methods in the field. We discuss how multi-omics technologies have been adapted and improved over the past decade using a framework characterized by optimization of throughput and resolution, modality integration, uniqueness and accuracy, and we also discuss multi-omics limitations. We highlight the impact that single-cell multi-omics technologies have had in cell lineage tracing, tissue-specific and cell-specific atlas production, tumour immunology and cancer genetics, and in mapping of cellular spatial information in fundamental and translational research. Finally, we discuss bioinformatics tools that have been developed to link different omics modalities and elucidate functionality through the use of better mathematical modelling and computational methods.

scBridge embraces cell heterogeneity in single-cell RNA-seq and ATAC-seq data integration

Single-cell multi-omics data integration aims to reduce the omics difference while keeping the cell type difference. However, it is daunting to model and distinguish the two differences due to cell heterogeneity. Namely, even cells of the same omics and type would have various features, making the two differences less significant. In this work, we reveal that instead of being an interference, cell heterogeneity could be exploited to improve data integration. Specifically, we observe that the omics difference varies in cells, and cells with smaller omics differences are easier to be integrated. Hence, unlike most existing works that homogeneously treat and integrate all cells, we propose a multi-omics data integration method (dubbed scBridge) that integrates cells in a heterogeneous manner. In brief, scBridge iterates between i) identifying reliable scATAC-seq cells that have smaller omics differences, and ii) integrating reliable scATAC-seq cells with scRNA-seq data to narrow the omics gap, thus benefiting the integration for the rest cells. Extensive experiments on seven multi-omics datasets demonstrate the superiority of scBridge compared with six representative baselines.

AtacAnnoR: a reference-based annotation tool for single cell ATAC-seq data

Here, we present AtacAnnoR, a two-round annotation method for scATAC-seq data using well-annotated scRNA-seq data as reference. We evaluate AtacAnnoR's performance against six competing methods on 11 benchmark datasets. Our results show that AtacAnnoR achieves the highest mean accuracy and the highest mean balanced accuracy and performs particularly well when unpaired scRNA-seq data are used as the reference. Furthermore, AtacAnnoR implements a 'Combine and Discard' strategy to further improve annotation accuracy when annotations of multiple references are available. AtacAnnoR has been implemented in an R package and can be directly integrated into currently popular scATAC-seq analysis pipelines.

Multimodal deep learning approaches for single-cell multi-omics data integration

Integrating single-cell multi-omics data is a challenging task that has led to new insights into complex cellular systems. Various computational methods have been proposed to effectively integrate these rapidly accumulating datasets, including deep learning. However, despite the proven success of deep learning in integrating multi-omics data and its better performance over classical computational methods, there has been no systematic study of its application to single-cell multi-omics data integration. To fill this gap, we conducted a literature review to explore the use of multimodal deep learning techniques in single-cell multi-omics data integration, taking into account recent studies from multiple perspectives. Specifically, we first summarized different modalities found in single-cell multi-omics data. We then reviewed current deep learning techniques for processing multimodal data and categorized deep learning-based integration methods for single-cell multi-omics data according to data modality, deep learning architecture, fusion strategy, key tasks and downstream analysis. Finally, we provided insights into using these deep learning models to integrate multi-omics data and better understand single-cell biological mechanisms.

Integrative analysis of single-cell RNA-seq and ATAC-seq reveals heterogeneity of induced pluripotent stem cell-derived hepatic organoids

To gain deeper insights into transcriptomes and epigenomes of organoids, liver organoids from two states (expandable and more differentiated) were subjected to single-cell RNA-seq (scRNA-seq) and single-cell ATAC-seq (scATAC-seq) analyses. Mitochondrial gene expression was higher in differentiated than in non-differentiated hepatocytes, with ATAC-seq peaks increasing near the mitochondrial control region. Differentiation of liver organoids resulted in the expression of transcription factors that act as enhancers and repressors. In addition, epigenetic mechanisms regulating the expression of alpha-fetoprotein (AFP) and albumin (ALB) differed in liver organoids and adult liver. Knockdown of PDX1, an essential transcription factor for pancreas development, led to the hepatic maturation of liver organoids through regulation of AFP and ALB expression. This integrative analysis of the transcriptomes and epigenomes of liver organoids at the single-cell level may contribute to a better understanding of the regulatory networks during liver development and the further development of mature in vitro human liver models.

Dynamic single-cell regulomes characterize human peripheral blood innate lymphoid cell subpopulations

Innate lymphoid cells (ILCs) are plastic immune cells divided into 3 main subsets, characterized by distinct phenotypic and functional profiles. Using single cell approaches, heightened heterogeneity of mouse ILCs has been appreciated, imprinted by tissue signals that shape their transcriptome and epigenome. Intra-subset diversity has also been observed in human ILCs. However, combined transcriptomic and epigenetic analyses of single ILCs in humans are lacking. Here, we show high transcriptional and epigenetic heterogeneity among human circulating ILCs in healthy individuals. We describe phenotypically distinct subclusters and diverse chromatin accessibility within main ILC populations, compatible with differentially poised states. We validate the use of this healthy donor-based analysis as resource dataset to help inferring ILC changes occurring in disease conditions. Overall, our work provides insights in the complex human ILC biology. We anticipate it to facilitate hypothesis-driven studies in patients, without the need to perform single cell OMICs using precious patients' material.

Polygenic regression uncovers trait-relevant cellular contexts through pathway activation transformation of single-cell RNA sequencing data

Advances in single-cell RNA sequencing (scRNA-seq) techniques have accelerated functional interpretation of disease-associated variants discovered from genome-wide association studies (GWASs). However, identification of trait-relevant cell populations is often impeded by inherent technical noise and high sparsity in scRNA-seq data. Here, we developed scPagwas, a computational approach that uncovers trait-relevant cellular context by integrating pathway activation transformation of scRNA-seq data and GWAS summary statistics. scPagwas effectively prioritizes trait-relevant genes, which facilitates identification of trait-relevant cell types/populations with high accuracy in extensive simulated and real datasets. Cellular-level association results identified a novel subpopulation of naive CD8+ T cells related to COVID-19 severity and oligodendrocyte progenitor cell and microglia subsets with critical pathways by which genetic variants influence Alzheimer's disease. Overall, our approach provides new insights for the discovery of trait-relevant cell types and improves the mechanistic understanding of disease variants from a pathway perspective.