AMULET: a novel read count-based method for effective multiplet detection from single nucleus ATAC-seq data

Mapping genotypes to chromatin accessibility profiles in single cells

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

Single cell regulatory architecture of human pancreatic islets suggests sex differences in β cell function and the pathogenesis of type 2 diabetes

Biological sex affects the pathogenesis of type 2 and type 1 diabetes (T2D, T1D) including the development of β cell failure observed more often in males. The mechanisms that drive sex differences in β cell failure is unknown. Studying sex differences in islet regulation and function represent a unique avenue to understand the sex-specific heterogeneity in β cell failure in diabetes. Here, we examined sex and race differences in human pancreatic islets from up to 52 donors with and without T2D (including 37 donors from the Human Pancreas Analysis Program [HPAP] dataset) using an orthogonal series of experiments including single cell RNA-seq (scRNA-seq), single nucleus assay for transposase-accessible chromatin sequencing (snATAC-seq), dynamic hormone secretion, and bioenergetics. In cultured islets from nondiabetic (ND) donors, in the absence of the in vivo hormonal environment, sex differences in islet cell type gene accessibility and expression predominantly involved sex chromosomes. Of particular interest were sex differences in the X-linked KDM6A and Y-linked KDM5D chromatin remodelers in female and male islet cells respectively. Islets from T2D donors exhibited similar sex differences in differentially expressed genes (DEGs) from sex chromosomes. However, in contrast to islets from ND donors, islets from T2D donors exhibited major sex differences in DEGs from autosomes. Comparing β cells from T2D and ND donors revealed that females had more DEGs from autosomes compared to male β cells. Gene set enrichment analysis of female β cell DEGs showed a suppression of oxidative phosphorylation and electron transport chain pathways, while male β cell had suppressed insulin secretion pathways. Thus, although sex-specific differences in gene accessibility and expression of cultured ND human islets predominantly affect sex chromosome genes, major differences in autosomal gene expression between sexes appear during the transition to T2D and which highlight mitochondrial failure in female β cells.

scifi-ATAC-seq: massive-scale single-cell chromatin accessibility sequencing using combinatorial fluidic indexing

Single-cell ATAC-seq has emerged as a powerful approach for revealing candidate cis-regulatory elements genome-wide at cell-type resolution. However, current single-cell methods suffer from limited throughput and high costs. Here, we present a novel technique called scifi-ATAC-seq, single-cell combinatorial fluidic indexing ATAC-sequencing, which combines a barcoded Tn5 pre-indexing step with droplet-based single-cell ATAC-seq using the 10X Genomics platform. With scifi-ATAC-seq, up to 200,000 nuclei across multiple samples can be indexed in a single emulsion reaction, representing an approximately 20-fold increase in throughput compared to the standard 10X Genomics workflow.

ISL1 and POU4F1 Directly Interact to Regulate the Differentiation and Survival of Inner Ear Sensory Neurons

The inner ear sensory neurons play a pivotal role in auditory processing and balance control. Though significant progresses have been made, the underlying mechanisms controlling the differentiation and survival of the inner ear sensory neurons remain largely unknown. During development, ISL1 and POU4F transcription factors are co-expressed and are required for terminal differentiation, pathfinding, axon outgrowth and the survival of neurons in the central and peripheral nervous systems. However, little is understood about their functional relationship and regulatory mechanism in neural development. Here, we have knocked out Isl1 or Pou4f1 or both in mice of both sexes. In the absence of Isl1, the differentiation of cochleovestibular ganglion (CVG) neurons is disturbed and with that Isl1-deficient CVG neurons display defects in migration and axon pathfinding. Compound deletion of Isl1 and Pou4f1 causes a delay in CVG differentiation and results in a more severe CVG defect with a loss of nearly all of spiral ganglion neurons (SGNs). Moreover, ISL1 and POU4F1 interact directly in developing CVG neurons and act cooperatively as well as independently in regulating the expression of unique sets of CVG-specific genes crucial for CVG development and survival by binding to the cis-regulatory elements including the promoters of Fgf10, Pou4f2, and Epha5 and enhancers of Eya1 and Ntng2 These findings demonstrate that Isl1 and Pou4f1 are indispensable for CVG development and maintenance by acting epistatically to regulate genes essential for CVG development.

Massive-scale single-cell chromatin accessibility sequencing using combinatorial fluidic indexing

Single-cell ATAC-seq has emerged as a powerful approach for revealing candidate cis-regulatory elements genome-wide at cell-type resolution. However, current single-cell methods suffer from limited throughput and high costs. Here, we present a novel technique called single-cell combinatorial fluidic indexing ATAC-sequencing ("scifi-ATAC-seq"), which combines a barcoded Tn5 pre-indexing step with droplet-based single-cell ATAC-seq using a widely commercialized microfluidics platform (10X Genomics). With scifi-ATAC-seq, up to 200,000 nuclei across multiple samples in a single emulsion reaction can be indexed, representing a ~20-fold increase in throughput compared to the standard 10X Genomics workflow.

Predicting proximal tubule failed repair drivers through regularized regression analysis of single cell multiomic sequencing

Renal proximal tubule epithelial cells have considerable intrinsic repair capacity following injury. However, a fraction of injured proximal tubule cells fails to undergo normal repair and assumes a proinflammatory and profibrotic phenotype that may promote fibrosis and chronic kidney disease. The healthy to failed repair change is marked by cell state-specific transcriptomic and epigenomic changes. Single nucleus joint RNA- and ATAC-seq sequencing offers an opportunity to study the gene regulatory networks underpinning these changes in order to identify key regulatory drivers. We develop a regularized regression approach to construct genome-wide parametric gene regulatory networks using multiomic datasets. We generate a single nucleus multiomic dataset from seven adult human kidney samples and apply our method to study drivers of a failed injury response associated with kidney disease. We demonstrate that our approach is a highly effective tool for predicting key cis- and trans-regulatory elements underpinning the healthy to failed repair transition and use it to identify NFAT5 as a driver of the maladaptive proximal tubule state.

Mosaic loss of Y chromosome is associated with aging and epithelial injury in chronic kidney disease

BACKGROUND

Mosaic loss of Y chromosome (LOY) is the most common chromosomal alteration in aging men. Here, we use single-cell RNA and ATAC sequencing to show that LOY is present in the kidney and increases with age and chronic kidney disease.

RESULTS

The likelihood of a cell having LOY varies depending on its location in the nephron. Cortical epithelial cell types have a greater proportion of LOY than medullary or glomerular cell types, which may reflect their proliferative history. Proximal tubule cells are the most abundant cell type in the cortex and are susceptible to hypoxic injury. A subset of these cells acquires a pro-inflammatory transcription and chromatin accessibility profile associated with expression of HAVCR1, VCAM1, and PROM1. These injured epithelial cells have the greatest proportion of LOY and their presence predicts future kidney function decline. Moreover, proximal tubule cells with LOY are more likely to harbor additional large chromosomal gains and express pro-survival pathways. Spatial transcriptomics localizes injured proximal tubule cells to a pro-fibrotic microenvironment where they adopt a secretory phenotype and likely communicate with infiltrating immune cells.

CONCLUSIONS

We hypothesize that LOY is an indicator of increased DNA damage and potential marker of cellular senescence that can be applied to single-cell datasets in other tissues.

Epigenetic reprogramming driving successful and failed repair in acute kidney injury

Acute kidney injury (AKI) causes epithelial damage followed by subsequent repair. While successful repair restores kidney function, this process is often incomplete and can lead to chronic kidney disease (CKD) in a process called failed repair. To better understand the epigenetic reprogramming driving this AKI-to-CKD transition we generated a single nucleus multiomic atlas for the full mouse AKI time course, consisting of ~280,000 single nucleus transcriptomes and epigenomes. We reveal cell-specific dynamic alterations in gene regulatory landscapes reflecting especially activation of proinflammatory pathways. We further generated single nucleus multiomic data from four human AKI samples including validation by genome-wide identification of NF-kB binding sites. A regularized regression analysis identifies key regulators involved in both successful and failed repair cell fate, identifying the transcription factor CREB5 as a regulator of both successful and failed tubular repair that also drives proximal tubule cell proliferation after injury. Our interspecies multiomic approach provides a foundation to comprehensively understand cell states in AKI.

ROR2 regulates cellular plasticity in pancreatic neoplasia and adenocarcinoma

Cellular plasticity is a hallmark of pancreatic ductal adenocarcinoma (PDAC) starting from the conversion of normal cells into precancerous lesions to the progression of carcinoma subtypes associated with aggressiveness and therapeutic response. We discovered that normal acinar cell differentiation, maintained by the transcription factor Pdx1, suppresses a broad gastric cell identity that is maintained in metaplasia, neoplasia, and the classical subtype of PDAC in mouse and human. We have identified the receptor tyrosine kinase Ror2 as marker of a gastric metaplasia (SPEM)-like identity in the pancreas. Ablation of Ror2 in a mouse model of pancreatic tumorigenesis promoted a switch to a gastric pit cell identity that largely persisted through progression to the classical subtype of PDAC. In both human and mouse pancreatic cancer, ROR2 activity continued to antagonize the gastric pit cell identity, strongly promoting an epithelial to mesenchymal transition, conferring resistance to KRAS inhibition, and vulnerability to AKT inhibition.

Uniform quantification of single-nucleus ATAC-seq data with Paired-Insertion Counting (PIC) and a model-based insertion rate estimator

Existing approaches to scoring single-nucleus assay for transposase-accessible chromatin with sequencing (snATAC-seq) feature matrices from sequencing reads are inconsistent, affecting downstream analyses and displaying artifacts. We show that, even with sparse single-cell data, quantitative counts are informative for estimating the regulatory state of a cell, which calls for a consistent treatment. We propose Paired-Insertion Counting as a uniform method for snATAC-seq feature characterization and provide a probability model for inferring latent insertion dynamics from snATAC-seq count matrices.

Single-nucleus multiomic atlas of frontal cortex in amyotrophic lateral sclerosis with a deep learning-based decoding of alternative polyadenylation mechanisms

The understanding of how different cell types contribute to amyotrophic lateral sclerosis (ALS) pathogenesis is limited. Here we generated a single-nucleus transcriptomic and epigenomic atlas of the frontal cortex of ALS cases with C9orf72 (C9) hexanucleotide repeat expansions and sporadic ALS (sALS). Our findings reveal shared pathways in C9-ALS and sALS, characterized by synaptic dysfunction in excitatory neurons and a disease-associated state in microglia. The disease subtypes diverge with loss of astrocyte homeostasis in C9-ALS, and a more substantial disturbance of inhibitory neurons in sALS. Leveraging high depth 3'-end sequencing, we found a widespread switch towards distal polyadenylation (PA) site usage across ALS subtypes relative to controls. To explore this differential alternative PA (APA), we developed APA-Net, a deep neural network model that uses transcript sequence and expression levels of RNA-binding proteins (RBPs) to predict cell-type specific APA usage and RBP interactions likely to regulate APA across disease subtypes.

Single-cell analysis of chromatin accessibility in the adult mouse brain

Recent advances in single-cell technologies have led to the discovery of thousands of brain cell types; however, our understanding of the gene regulatory programs in these cell types is far from complete1-4. Here we report a comprehensive atlas of candidate cis-regulatory DNA elements (cCREs) in the adult mouse brain, generated by analysing chromatin accessibility in 2.3 million individual brain cells from 117 anatomical dissections. The atlas includes approximately 1 million cCREs and their chromatin accessibility across 1,482 distinct brain cell populations, adding over 446,000 cCREs to the most recent such annotation in the mouse genome. The mouse brain cCREs are moderately conserved in the human brain. The mouse-specific cCREs-specifically, those identified from a subset of cortical excitatory neurons-are strongly enriched for transposable elements, suggesting a potential role for transposable elements in the emergence of new regulatory programs and neuronal diversity. Finally, we infer the gene regulatory networks in over 260 subclasses of mouse brain cells and develop deep-learning models to predict the activities of gene regulatory elements in different brain cell types from the DNA sequence alone. Our results provide a resource for the analysis of cell-type-specific gene regulation programs in both mouse and human brains.

Epicardioid single-cell genomics uncovers principles of human epicardium biology in heart development and disease

The epicardium, the mesothelial envelope of the vertebrate heart, is the source of multiple cardiac cell lineages during embryonic development and provides signals that are essential to myocardial growth and repair. Here we generate self-organizing human pluripotent stem cell-derived epicardioids that display retinoic acid-dependent morphological, molecular and functional patterning of the epicardium and myocardium typical of the left ventricular wall. By combining lineage tracing, single-cell transcriptomics and chromatin accessibility profiling, we describe the specification and differentiation process of different cell lineages in epicardioids and draw comparisons to human fetal development at the transcriptional and morphological levels. We then use epicardioids to investigate the functional cross-talk between cardiac cell types, gaining new insights into the role of IGF2/IGF1R and NRP2 signaling in human cardiogenesis. Finally, we show that epicardioids mimic the multicellular pathogenesis of congenital or stress-induced hypertrophy and fibrotic remodeling. As such, epicardioids offer a unique testing ground of epicardial activity in heart development, disease and regeneration.

Genetic risk converges on regulatory networks mediating early type 2 diabetes

Type 2 diabetes mellitus (T2D), a major cause of worldwide morbidity and mortality, is characterized by dysfunction of insulin-producing pancreatic islet β cells1,2. T2D genome-wide association studies (GWAS) have identified hundreds of signals in non-coding and β cell regulatory genomic regions, but deciphering their biological mechanisms remains challenging3-5. Here, to identify early disease-driving events, we performed traditional and multiplexed pancreatic tissue imaging, sorted-islet cell transcriptomics and islet functional analysis of early-stage T2D and control donors. By integrating diverse modalities, we show that early-stage T2D is characterized by β cell-intrinsic defects that can be proportioned into gene regulatory modules with enrichment in signals of genetic risk. After identifying the β cell hub gene and transcription factor RFX6 within one such module, we demonstrated multiple layers of genetic risk that converge on an RFX6-mediated network to reduce insulin secretion by β cells. RFX6 perturbation in primary human islet cells alters β cell chromatin architecture at regions enriched for T2D GWAS signals, and population-scale genetic analyses causally link genetically predicted reduced RFX6 expression with increased T2D risk. Understanding the molecular mechanisms of complex, systemic diseases necessitates integration of signals from multiple molecules, cells, organs and individuals, and thus we anticipate that this approach will be a useful template to identify and validate key regulatory networks and master hub genes for other diseases or traits using GWAS data.

Comprehensive single cell aging atlas of mammary tissues reveals shared epigenomic and transcriptomic signatures of aging and cancer

Aging is the greatest risk factor for breast cancer; however, how age-related cellular and molecular events impact cancer initiation is unknown. We investigate how aging rewires transcriptomic and epigenomic programs of mouse mammary glands at single cell resolution, yielding a comprehensive resource for aging and cancer biology. Aged epithelial cells exhibit epigenetic and transcriptional changes in metabolic, pro-inflammatory, or cancer-associated genes. Aged stromal cells downregulate fibroblast marker genes and upregulate markers of senescence and cancer-associated fibroblasts. Among immune cells, distinct T cell subsets (Gzmk+, memory CD4+, γδ) and M2-like macrophages expand with age. Spatial transcriptomics reveal co-localization of aged immune and epithelial cells in situ. Lastly, transcriptional signatures of aging mammary cells are found in human breast tumors, suggesting mechanistic links between aging and cancer. Together, these data uncover that epithelial, immune, and stromal cells shift in proportions and cell identity, potentially impacting cell plasticity, aged microenvironment, and neoplasia risk.

Transcription factor stoichiometry, motif affinity and syntax regulate single-cell chromatin dynamics during fibroblast reprogramming to pluripotency

Ectopic expression of OCT4, SOX2, KLF4 and MYC (OSKM) transforms differentiated cells into induced pluripotent stem cells. To refine our mechanistic understanding of reprogramming, especially during the early stages, we profiled chromatin accessibility and gene expression at single-cell resolution across a densely sampled time course of human fibroblast reprogramming. Using neural networks that map DNA sequence to ATAC-seq profiles at base-resolution, we annotated cell-state-specific predictive transcription factor (TF) motif syntax in regulatory elements, inferred affinity- and concentration-dependent dynamics of Tn5-bias corrected TF footprints, linked peaks to putative target genes, and elucidated rewiring of TF-to-gene cis-regulatory networks. Our models reveal that early in reprogramming, OSK, at supraphysiological concentrations, rapidly open transient regulatory elements by occupying non-canonical low-affinity binding sites. As OSK concentration falls, the accessibility of these transient elements decays as a function of motif affinity. We find that these OSK-dependent transient elements sequester the somatic TF AP-1. This redistribution is strongly associated with the silencing of fibroblast-specific genes within individual nuclei. Together, our integrated single-cell resource and models reveal insights into the cis-regulatory code of reprogramming at unprecedented resolution, connect TF stoichiometry and motif syntax to diversification of cell fate trajectories, and provide new perspectives on the dynamics and role of transient regulatory elements in somatic silencing.

scIBD: a self-supervised iterative-optimizing model for boosting the detection of heterotypic doublets in single-cell chromatin accessibility data

Application of the widely used droplet-based microfluidic technologies in single-cell sequencing often yields doublets, introducing bias to downstream analyses. Especially, doublet-detection methods for single-cell chromatin accessibility sequencing (scCAS) data have multiple assay-specific challenges. Therefore, we propose scIBD, a self-supervised iterative-optimizing model for boosting heterotypic doublet detection in scCAS data. scIBD introduces an adaptive strategy to simulate high-confident heterotypic doublets and self-supervise for doublet-detection in an iteratively optimizing manner. Comprehensive benchmarking on various simulated and real datasets demonstrates the outperformance and robustness of scIBD. Moreover, the downstream biological analyses suggest the efficacy of doublet-removal by scIBD.

Single Cell Multi-Omics of an iPSC Model of Human Sinoatrial Node Development Reveals Genetic Determinants of Heart Rate and Arrhythmia Susceptibility

Cellular heterogeneity within the sinoatrial node (SAN) is functionally important but has been difficult to model in vitro , presenting a major obstacle to studies of heart rate regulation and arrhythmias. Here we describe a scalable method to derive sinoatrial node pacemaker cardiomyocytes (PCs) from human induced pluripotent stem cells that recapitulates differentiation into distinct PC subtypes, including SAN Head, SAN Tail, transitional zone cells, and sinus venosus myocardium. Single cell (sc) RNA-sequencing, sc-ATAC-sequencing, and trajectory analyses were used to define epigenetic and transcriptomic signatures of each cell type, and to identify novel transcriptional pathways important for PC subtype differentiation. Integration of our multi-omics datasets with genome wide association studies uncovered cell type-specific regulatory elements that associated with heart rate regulation and susceptibility to atrial fibrillation. Taken together, these datasets validate a novel, robust, and realistic in vitro platform that will enable deeper mechanistic exploration of human cardiac automaticity and arrhythmia.

Single-cell multiomic understanding of HIV-1 reservoir at epigenetic, transcriptional, and protein levels

PURPOSE OF REVIEW

The success of HIV-1 eradication strategies relies on in-depth understanding of HIV-1-infected cells. However, HIV-1-infected cells are extremely heterogeneous and rare. Single-cell multiomic approaches resolve the heterogeneity and rarity of HIV-1-infected cells.

RECENT FINDINGS

Advancement in single-cell multiomic approaches enabled HIV-1 reservoir profiling across the epigenetic (ATAC-seq), transcriptional (RNA-seq), and protein levels (CITE-seq). Using HIV-1 RNA as a surrogate, ECCITE-seq identified enrichment of HIV-1-infected cells in clonally expanded cytotoxic CD4+ T cells. Using HIV-1 DNA PCR-activated microfluidic sorting, FIND-seq captured the bulk transcriptome of HIV-1 DNA+ cells. Using targeted HIV-1 DNA amplification, PheP-seq identified surface protein expression of intact versus defective HIV-1-infected cells. Using ATAC-seq to identify HIV-1 DNA, ASAP-seq captured transcription factor activity and surface protein expression of HIV-1 DNA+ cells. Combining HIV-1 mapping by ATAC-seq and HIV-1 RNA mapping by RNA-seq, DOGMA-seq captured the epigenetic, transcriptional, and surface protein expression of latent and transcriptionally active HIV-1-infected cells. To identify reproducible biological insights and authentic HIV-1-infected cells and avoid false-positive discovery of artifacts, we reviewed current practices of single-cell multiomic experimental design and bioinformatic analysis.

SUMMARY

Single-cell multiomic approaches may identify innovative mechanisms of HIV-1 persistence, nominate therapeutic strategies, and accelerate discoveries.

Epigenetic memory of coronavirus infection in innate immune cells and their progenitors

Inflammation can trigger lasting phenotypes in immune and non-immune cells. Whether and how human infections and associated inflammation can form innate immune memory in hematopoietic stem and progenitor cells (HSPC) has remained unclear. We found that circulating HSPC, enriched from peripheral blood, captured the diversity of bone marrow HSPC, enabling investigation of their epigenomic reprogramming following coronavirus disease 2019 (COVID-19). Alterations in innate immune phenotypes and epigenetic programs of HSPC persisted for months to 1 year following severe COVID-19 and were associated with distinct transcription factor (TF) activities, altered regulation of inflammatory programs, and durable increases in myelopoiesis. HSPC epigenomic alterations were conveyed, through differentiation, to progeny innate immune cells. Early activity of IL-6 contributed to these persistent phenotypes in human COVID-19 and a mouse coronavirus infection model. Epigenetic reprogramming of HSPC may underlie altered immune function following infection and be broadly relevant, especially for millions of COVID-19 survivors.

Best practices for single-cell analysis across modalities

Recent advances in single-cell technologies have enabled high-throughput molecular profiling of cells across modalities and locations. Single-cell transcriptomics data can now be complemented by chromatin accessibility, surface protein expression, adaptive immune receptor repertoire profiling and spatial information. The increasing availability of single-cell data across modalities has motivated the development of novel computational methods to help analysts derive biological insights. As the field grows, it becomes increasingly difficult to navigate the vast landscape of tools and analysis steps. Here, we summarize independent benchmarking studies of unimodal and multimodal single-cell analysis across modalities to suggest comprehensive best-practice workflows for the most common analysis steps. Where independent benchmarks are not available, we review and contrast popular methods. Our article serves as an entry point for novices in the field of single-cell (multi-)omic analysis and guides advanced users to the most recent best practices.

SWI/SNF Blockade Disrupts PU.1-Directed Enhancer Programs in Normal Hematopoietic Cells and Acute Myeloid Leukemia

In acute myeloid leukemia (AML), SWI/SNF chromatin remodeling complexes sustain leukemic identity by driving high levels of MYC. Previous studies have implicated the hematopoietic transcription factor PU.1 (SPI1) as an important target of SWI/SNF inhibition, but PU.1 is widely regarded to have pioneer-like activity. As a result, many questions have remained regarding the interplay between PU.1 and SWI/SNF in AML as well as normal hematopoiesis. Here we found that PU.1 binds to most of its targets in a SWI/SNF-independent manner and recruits SWI/SNF to promote accessibility for other AML core regulatory factors, including RUNX1, LMO2, and MEIS1. SWI/SNF inhibition in AML cells reduced DNA accessibility and binding of these factors at PU.1 sites and redistributed PU.1 to promoters. Analysis of nontumor hematopoietic cells revealed that similar effects also impair PU.1-dependent B-cell and monocyte populations. Nevertheless, SWI/SNF inhibition induced profound therapeutic response in an immunocompetent AML mouse model as well as in primary human AML samples. In vivo, SWI/SNF inhibition promoted leukemic differentiation and reduced the leukemic stem cell burden in bone marrow but also induced leukopenia. These results reveal a variable therapeutic window for SWI/SNF blockade in AML and highlight important off-tumor effects of such therapies in immunocompetent settings.

SIGNIFICANCE

Disruption of PU.1-directed enhancer programs upon SWI/SNF inhibition causes differentiation of AML cells and induces leukopenia of PU.1-dependent B cells and monocytes, revealing the on- and off-tumor effects of SWI/SNF blockade.

Profound phenotypic and epigenetic heterogeneity of the HIV-1-infected CD4+ T cell reservoir

Understanding the complexity of the long-lived HIV reservoir during antiretroviral therapy (ART) remains a considerable impediment in research towards a cure for HIV. To address this, we developed a single-cell strategy to precisely define the unperturbed peripheral blood HIV-infected memory CD4+ T cell reservoir from ART-treated people living with HIV (ART-PLWH) via the presence of integrated accessible proviral DNA in concert with epigenetic and cell surface protein profiling. We identified profound reservoir heterogeneity within and between ART-PLWH, characterized by new and known surface markers within total and individual memory CD4+ T cell subsets. We further uncovered new epigenetic profiles and transcription factor motifs enriched in HIV-infected cells that suggest infected cells with accessible provirus, irrespective of reservoir distribution, are poised for reactivation during ART treatment. Together, our findings reveal the extensive inter- and intrapersonal cellular heterogeneity of the HIV reservoir, and establish an initial multiomic atlas to develop targeted reservoir elimination strategies.

Characterizing cis-regulatory elements using single-cell epigenomics

Cell type-specific gene expression patterns and dynamics during development or in disease are controlled by cis-regulatory elements (CREs), such as promoters and enhancers. Distinct classes of CREs can be characterized by their epigenomic features, including DNA methylation, chromatin accessibility, combinations of histone modifications and conformation of local chromatin. Tremendous progress has been made in cataloguing CREs in the human genome using bulk transcriptomic and epigenomic methods. However, single-cell epigenomic and multi-omic technologies have the potential to provide deeper insight into cell type-specific gene regulatory programmes as well as into how they change during development, in response to environmental cues and through disease pathogenesis. Here, we highlight recent advances in single-cell epigenomic methods and analytical tools and discuss their readiness for human tissue profiling.

scATACpipe: A nextflow pipeline for comprehensive and reproducible analyses of single cell ATAC-seq data

Single cell ATAC-seq (scATAC-seq) has become the most widely used method for profiling open chromatin landscape of heterogeneous cell populations at a single-cell resolution. Although numerous software tools and pipelines have been developed, an easy-to-use, scalable, reproducible, and comprehensive pipeline for scATAC-seq data analyses is still lacking. To fill this gap, we developed scATACpipe, a Nextflow pipeline, for performing comprehensive analyses of scATAC-seq data including extensive quality assessment, preprocessing, dimension reduction, clustering, peak calling, differential accessibility inference, integration with scRNA-seq data, transcription factor activity and footprinting analysis, co-accessibility inference, and cell trajectory prediction. scATACpipe enables users to perform the end-to-end analysis of scATAC-seq data with three sub-workflow options for preprocessing that leverage 10x Genomics Cell Ranger ATAC software, the ultra-fast Chromap procedures, and a set of custom scripts implementing current best practices for scATAC-seq data preprocessing. The pipeline extends the R package ArchR for downstream analysis with added support to any eukaryotic species with an annotated reference genome. Importantly, scATACpipe generates an all-in-one HTML report for the entire analysis and outputs cluster-specific BAM, BED, and BigWig files for visualization in a genome browser. scATACpipe eliminates the need for users to chain different tools together and facilitates reproducible and comprehensive analyses of scATAC-seq data from raw reads to various biological insights with minimal changes of configuration settings for different computing environments or species. By applying it to public datasets, we illustrated the utility, flexibility, versatility, and reliability of our pipeline, and demonstrated that our scATACpipe outperforms other workflows.

Multimodal single cell sequencing implicates chromatin accessibility and genetic background in diabetic kidney disease progression

The proximal tubule is a key regulator of kidney function and glucose metabolism. Diabetic kidney disease leads to proximal tubule injury and changes in chromatin accessibility that modify the activity of transcription factors involved in glucose metabolism and inflammation. Here we use single nucleus RNA and ATAC sequencing to show that diabetic kidney disease leads to reduced accessibility of glucocorticoid receptor binding sites and an injury-associated expression signature in the proximal tubule. We hypothesize that chromatin accessibility is regulated by genetic background and closely-intertwined with metabolic memory, which pre-programs the proximal tubule to respond differently to external stimuli. Glucocorticoid excess has long been known to increase risk for type 2 diabetes, which raises the possibility that glucocorticoid receptor inhibition may mitigate the adverse metabolic effects of diabetic kidney disease.

Fundamental and practical approaches for single-cell ATAC-seq analysis

Assays for transposase-accessible chromatin through high-throughput sequencing (ATAC-seq) are effective tools in the study of genome-wide chromatin accessibility landscapes. With the rapid development of single-cell technology, open chromatin regions that play essential roles in epigenetic regulation have been measured at the single-cell level using single-cell ATAC-seq approaches. The application of scATAC-seq has become as popular as that of scRNA-seq. However, owing to the nature of scATAC-seq data, which are sparse and noisy, processing the data requires different methodologies and empirical experience. This review presents a practical guide for processing scATAC-seq data, from quality evaluation to downstream analysis, for various applications. In addition to the epigenomic profiling from scATAC-seq, we also discuss recent studies in which the function of non-coding variants has been investigated based on cell type-specific cis-regulatory elements and how to use the by-product genetic information obtained from scATAC-seq to infer single-cell copy number variants and trace cell lineage. We anticipate that this review will assist researchers in designing and implementing scATAC-seq assays to facilitate research in diverse fields.

Glucocorticoid Receptor-Regulated Enhancers Play a Central Role in the Gene Regulatory Networks Underlying Drug Addiction

Substance abuse and addiction represent a significant public health problem that impacts multiple dimensions of society, including healthcare, the economy, and the workforce. In 2021, over 100,000 drug overdose deaths were reported in the US, with an alarming increase in fatalities related to opioids and psychostimulants. Understanding the fundamental gene regulatory mechanisms underlying addiction and related behaviors could facilitate more effective treatments. To explore how repeated drug exposure alters gene regulatory networks in the brain, we combined capped small (cs)RNA-seq, which accurately captures nascent-like initiating transcripts from total RNA, with Hi-C and single nuclei (sn)ATAC-seq. We profiled initiating transcripts in two addiction-related brain regions, the prefrontal cortex (PFC) and the nucleus accumbens (NAc), from rats that were never exposed to drugs or were subjected to prolonged abstinence after oxycodone or cocaine intravenous self-administration (IVSA). Interrogating over 100,000 active transcription start regions (TSRs) revealed that most TSRs had hallmarks of bonafide enhancers and highlighted the KLF/SP1, RFX, and AP1 transcription factors families as central to establishing brain-specific gene regulatory programs. Analysis of rats with addiction-like behaviors versus controls identified addiction-associated repression of transcription at regulatory enhancers recognized by nuclear receptor subfamily 3 group C (NR3C) factors, including glucocorticoid receptors. Cell-type deconvolution analysis using snATAC-seq uncovered a potential role of glial cells in driving the gene regulatory programs associated with addiction-related phenotypes. These findings highlight the power of advanced transcriptomics methods to provide insight into how addiction perturbs gene regulatory programs in the brain.

Doublet identification in single-cell sequencing data using scDblFinder

Doublets are prevalent in single-cell sequencing data and can lead to artifactual findings. A number of strategies have therefore been proposed to detect them. Building on the strengths of existing

approaches, we developed scDblFinder, a fast, flexible and accurate Bioconductor-based doublet detection method. Here we present the method, justify its design choices, demonstrate its performance on both single-cell RNA and accessibility (ATAC) sequencing data, and provide some observations on doublet formation, detection, and enrichment analysis. Even in complex datasets, scDblFinder can accurately identify most heterotypic doublets, and was already found by an independent benchmark to outcompete alternatives.

Doublet identification in single-cell sequencing data using scDblFinder

Doublets are prevalent in single-cell sequencing data and can lead to artifactual findings. A number of strategies have therefore been proposed to detect them. Building on the strengths of existing approaches, we developed scDblFinder, a fast, flexible and accurate Bioconductor-based doublet detection method. Here we present the method, justify its design choices, demonstrate its performance on both single-cell RNA and accessibility (ATAC) sequencing data, and provide some observations on doublet formation, detection, and enrichment analysis. Even in complex datasets, scDblFinder can accurately identify most heterotypic doublets, and was already found by an independent benchmark to outcompete alternatives.