chromVAR: inferring transcription-factor-associated accessibility from single-cell epigenomic data

A Novel Macrophage Subpopulation Conveys Increased Genetic Risk of Coronary Artery Disease

BACKGROUND

Coronary artery disease (CAD), the leading cause of death worldwide, is influenced by both environmental and genetic factors. Although over 250 genetic risk loci have been identified through genome-wide association studies, the specific causal variants and their regulatory mechanisms are still largely unknown, particularly in disease-relevant cell types such as macrophages.

METHODS

We utilized single-cell RNA-seq and single-cell multiomics approaches in primary human monocyte-derived macrophages to explore the transcriptional regulatory network involved in a critical pathogenic event of coronary atherosclerosis-the formation of lipid-laden foam cells. The relative genetic contribution to CAD was assessed by partitioning disease heritability across different macrophage subpopulations. Meta-analysis of single-cell RNA-seq data sets from 38 human atherosclerotic samples was conducted to provide high-resolution cross-referencing to macrophage subpopulations in vivo.

RESULTS

We identified 18 782 cis-regulatory elements by jointly profiling the gene expression and chromatin accessibility of >5000 macrophages. Integration with CAD genome-wide association study data prioritized 121 CAD-related genetic variants and 56 candidate causal genes. We showed that CAD heritability was not uniformly distributed and was particularly enriched in the gene programs of a novel CD52-hi lipid-handling macrophage subpopulation. These CD52-hi macrophages displayed significantly less lipoprotein accumulation and were also found in human atherosclerotic plaques. We investigated the cis-regulatory effect of a risk variant rs10488763 on FDX1, implicating the recruitment of AP-1 and C/EBP-β in the causal mechanisms at this locus.

CONCLUSIONS

Our results provide genetic evidence of the divergent roles of macrophage subsets in atherogenesis and highlight lipid-handling macrophages as a key subpopulation through which genetic variants operate to influence disease. These findings provide an unbiased framework for functional fine-mapping of genome-wide association study results using single-cell multiomics and offer new insights into the genotype-environment interactions underlying atherosclerotic disease.

Osr2 functions as a biomechanical checkpoint to aggravate CD8+ T cell exhaustion in tumor

Alterations in extracellular matrix (ECM) architecture and stiffness represent hallmarks of cancer. Whether the biomechanical property of ECM impacts the functionality of tumor-reactive CD8+ T cells remains largely unknown. Here, we reveal that the transcription factor (TF) Osr2 integrates biomechanical signaling and facilitates the terminal exhaustion of tumor-reactive CD8+ T cells. Osr2 expression is selectively induced in the terminally exhausted tumor-specific CD8+ T cell subset by coupled T cell receptor (TCR) signaling and biomechanical stress mediated by the Piezo1/calcium/CREB axis. Consistently, depletion of Osr2 alleviates the exhaustion of tumor-specific CD8+ T cells or CAR-T cells, whereas forced Osr2 expression aggravates their exhaustion in solid tumor models. Mechanistically, Osr2 recruits HDAC3 to rewire the epigenetic program for suppressing cytotoxic gene expression and promoting CD8+ T cell exhaustion. Thus, our results unravel Osr2 functions as a biomechanical checkpoint to exacerbate CD8+ T cell exhaustion and could be targeted to potentiate cancer immunotherapy.

Single-cell and spatial atlases of spinal cord injury in the Tabulae Paralytica

Here, we introduce the Tabulae Paralytica-a compilation of four atlases of spinal cord injury (SCI) comprising a single-nucleus transcriptome atlas of half a million cells, a multiome atlas pairing transcriptomic and epigenomic measurements within the same nuclei, and two spatial transcriptomic atlases of the injured spinal cord spanning four spatial and temporal dimensions. We integrated these atlases into a common framework to dissect the molecular logic that governs the responses to injury within the spinal cord1. The Tabulae Paralytica uncovered new biological principles that dictate the consequences of SCI, including conserved and divergent neuronal responses to injury; the priming of specific neuronal subpopulations to upregulate circuit-reorganizing programs after injury; an inverse relationship between neuronal stress responses and the activation of circuit reorganization programs; the necessity of re-establishing a tripartite neuroprotective barrier between immune-privileged and extra-neural environments after SCI and a failure to form this barrier in old mice. We leveraged the Tabulae Paralytica to develop a rejuvenative gene therapy that re-established this tripartite barrier, and restored the natural recovery of walking after paralysis in old mice. The Tabulae Paralytica provides a window into the pathobiology of SCI, while establishing a framework for integrating multimodal, genome-scale measurements in four dimensions to study biology and medicine.

Identification of endothelial and mesenchymal FOXF1 enhancers involved in alveolar capillary dysplasia

Mutations in the FOXF1 gene, a key transcriptional regulator of pulmonary vascular development, cause Alveolar Capillary Dysplasia with Misalignment of Pulmonary Veins, a lethal lung disease affecting newborns and infants. Identification of new FOXF1 upstream regulatory elements is critical to explain why frequent non-coding FOXF1 deletions are linked to the disease. Herein, we use multiome single-nuclei RNA and ATAC sequencing of mouse and human patient lungs to identify four conserved endothelial and mesenchymal FOXF1 enhancers. We demonstrate that endothelial FOXF1 enhancers are autoactivated, whereas mesenchymal FOXF1 enhancers are regulated by EBF1 and GLI1. The cell-specificity of FOXF1 enhancers is validated by disrupting these enhancers in mouse embryonic stem cells using CRISPR/Cpf1 genome editing followed by lineage-tracing of mutant embryonic stem cells in mouse embryos using blastocyst complementation. This study resolves an important clinical question why frequent non-coding FOXF1 deletions that interfere with endothelial and mesenchymal enhancers can lead to the disease.

Single-Cell Chromatin Accessibility Analysis Reveals Subgroup-Specific TF-NTR Regulatory Circuits in Medulloblastoma

Medulloblastoma (MB) stands as one of the prevalent malignant brain tumors among pediatric patients. Despite its prevalence, the intricate interplay between the regulatory program driving malignancy in MB cells and their interactions with the microenvironment remains insufficiently understood. Leveraging the capabilities of single-cell Assay for Transposase-Accessible Chromatin sequencing (scATAC-seq), the chromatin accessibility landscape is unveiled across 59,015 distinct MB cells. This expansive dataset encompasses cells belonging to discrete molecular subgroups, namely SHH, WNT, Group3, and Group4. Within these chromatin accessibility profiles, specific regulatory elements tied to individual subgroups are uncovered, shedding light on the distinct activities of transcription factors (TFs) that likely orchestrate the tumorigenesis process. Moreover, it is found that certain neurotransmitter receptors (NTRs) are subgroup-specific and can predict MB subgroup classification when combined with their associated transcription factors. Notably, targeting essential NTRs within tumors influences both the in vitro sphere-forming capability and the in vivo tumorigenic capacity of MB cells. These findings collectively provide fresh insights into comprehending the regulatory networks and cellular dynamics within MBs. Furthermore, the significance of the TF-NTR regulatory circuits is underscored as prospective biomarkers and viable therapeutic targets.

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.

Spatially resolved multiomics on the neuronal effects induced by spaceflight in mice

Impairment of the central nervous system (CNS) poses a significant health risk for astronauts during long-duration space missions. In this study, we employed an innovative approach by integrating single-cell multiomics (transcriptomics and chromatin accessibility) with spatial transcriptomics to elucidate the impact of spaceflight on the mouse brain in female mice. Our comparative analysis between ground control and spaceflight-exposed animals revealed significant alterations in essential brain processes including neurogenesis, synaptogenesis and synaptic transmission, particularly affecting the cortex, hippocampus, striatum and neuroendocrine structures. Additionally, we observed astrocyte activation and signs of immune dysfunction. At the pathway level, some spaceflight-induced changes in the brain exhibit similarities with neurodegenerative disorders, marked by oxidative stress and protein misfolding. Our integrated spatial multiomics approach serves as a stepping stone towards understanding spaceflight-induced CNS impairments at the level of individual brain regions and cell types, and provides a basis for comparison in future spaceflight studies. For broader scientific impact, all datasets from this study are available through an interactive data portal, as well as the National Aeronautics and Space Administration (NASA) Open Science Data Repository (OSDR).

Defining the cellular origin of seminoma by transcriptional and epigenetic mapping to the normal human germline

Aberrant male germline development can lead to the formation of seminoma, a testicular germ cell tumor. Seminomas are biologically similar to primordial germ cells (PGCs) and many bear an isochromosome 12p [i(12p)] with two additional copies of the short arm of chromosome 12. By mapping seminoma transcriptomes and open chromatin landscape onto a normal human male germline trajectory, we find that seminoma resembles premigratory/migratory PGCs; however, it exhibits enhanced germline and pluripotency programs and upregulation of genes involved in apoptosis, angiogenesis, and MAPK/ERK pathways. Using pluripotent stem cell-derived PGCs from Pallister-Killian syndrome patients mosaic for i(12p), we model seminoma and identify gene dosage effects that may contribute to transformation. As murine seminoma models do not exist, our analyses provide critical insights into genetic, cellular, and signaling programs driving seminoma transformation, and the in vitro platform developed herein permits evaluation of additional signals required for seminoma tumorigenesis.

Low Kit expression identifies Hematopoietic Stem Cell subsets with enhanced lymphoid potential in mice and humans

Hematopoietic stem cells (HSCs) with multilineage potential are critical for effective T cell reconstitution and restoration of the adaptive immune system after allogeneic Hematopoietic Cell Transplantation (allo-HCT). Our current understanding of lymphoid-poised HSC subsets is limited due to a lack of reliable markers. Through comprehensive molecular profiling of HSCs, we identify Kit lo HSCs-a subset with enhanced lymphoid potential, which declines with age. Using a preclinical allo-HCT model, we demonstrate that Kit lo HSCs exhibit superior lymphoid progenitor output, enhanced thymic recovery, and T-cell reconstitution resulting in improved T cell responses to infection post-HCT. Furthermore, Kit lo HSCs with augmented BM lymphopoiesis mitigate age-associated thymic alterations, thus enhancing T-cell recovery in middle-aged hosts. Chromatin profiling revealed that Kit lo HSCs exhibit higher activity of lymphoid-specifying transcription factors (TFs), including Zbtb1 . Deletion of Zbtb1 in Kit lo HSCs, diminished their T-cell potential, while reinstating Zbtb1 in megakaryocytic-biased Kit hi HSCs rescued T-cell potential, in vitro and in vivo . Finally, we discover an analogous Kit lo HSC subset with enhanced lymphoid potential in human bone marrow. Our results demonstrate that Kit lo HSCs identify HSC subsets with enhanced lymphoid potential mediated by an underlying epigenetic program.

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.

Single-Cell Chromatin Accessibility Analysis Reveals the Epigenetic Basis and Signature Transcription Factors for the Molecular Subtypes of Colorectal Cancers

Colorectal cancer is a highly heterogeneous disease, with well-characterized subtypes based on genome, DNA methylome, and transcriptome signatures. To chart the epigenetic landscape of colorectal cancers, we generated a high-quality single-cell chromatin accessibility atlas of epithelial cells for 29 patients. Abnormal chromatin states acquired in adenomas were largely retained in colorectal cancers, which were tightly accompanied by opposite changes of DNA methylation. Unsupervised analysis on malignant cells revealed two epigenetic subtypes, exactly matching the iCMS classification, and key iCMS-specific transcription factors (TFs) were identified, including HNF4A and PPARA for iCMS2 tumors and FOXA3 and MAFK for iCMS3 tumors. Notably, subtype-specific TFs bind to distinct target gene sets and contribute to both interpatient similarities and diversities for both chromatin accessibilities and RNA expressions. Moreover, we identified CpG-island methylator phenotypes and pinpointed chromatin state signatures and TF regulators for the CIMP-high subtype. Our work systematically revealed the epigenetic basis of the well-known iCMS and CIMP classifications of colorectal cancers.

SIGNIFICANCE

Our work revealed the epigenetic basis of the well-known iCMS and CIMP classifications of colorectal cancers. Moreover, interpatient minor similarities and major diversities of chromatin accessibility signatures of TF target genes can faithfully explain the corresponding interpatient minor similarities and major diversities of RNA expression signatures of colorectal cancers, respectively. This article is featured in Selected Articles from This Issue, p. 897.

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

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

Targeting pro-inflammatory T cells as a novel therapeutic approach to potentially resolve atherosclerosis in humans

Atherosclerosis (AS), a leading cause of cardio-cerebrovascular disease worldwide, is driven by the accumulation of lipid contents and chronic inflammation. Traditional strategies primarily focus on lipid reduction to control AS progression, leaving residual inflammatory risks for major adverse cardiovascular events (MACEs). While anti-inflammatory therapies targeting innate immunity have reduced MACEs, many patients continue to face significant risks. Another key component in AS progression is adaptive immunity, but its potential role in preventing AS remains unclear. To investigate this, we conducted a retrospective cohort study on tumor patients with AS plaques. We found that anti-programmed cell death protein 1 (PD-1) monoclonal antibody (mAb) significantly reduces AS plaque size. With multi-omics single-cell analyses, we comprehensively characterized AS plaque-specific PD-1+ T cells, which are activated and pro-inflammatory. We demonstrated that anti-PD-1 mAb, when captured by myeloid-expressed Fc gamma receptors (FcγRs), interacts with PD-1 expressed on T cells. This interaction turns the anti-PD-1 mAb into a substitute PD-1 ligand, suppressing T-cell functions in the PD-1 ligands-deficient context of AS plaques. Further, we conducted a prospective cohort study on tumor patients treated with anti-PD-1 mAb with or without Fc-binding capability. Our analysis shows that anti-PD-1 mAb with Fc-binding capability effectively reduces AS plaque size, while anti-PD-1 mAb without Fc-binding capability does not. Our work suggests that T cell-targeting immunotherapy can be an effective strategy to resolve AS in humans.

SIMBA: single-cell embedding along with features

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

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

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

Multi-omics profiling of mouse polycystic kidney disease progression at a single cell resolution

Autosomal dominant polycystic kidney disease (ADPKD) is the most common hereditary kidney disease and causes significant morbidity, ultimately leading to end-stage kidney disease. PKD pathogenesis is characterized by complex and dynamic alterations in multiple cell types during disease progression, hampering a deeper understanding of disease mechanism and the development of therapeutic approaches. Here, we generate a single nucleus multimodal atlas of an orthologous mouse PKD model at early, mid and late timepoints, consisting of 125,434 single-nucleus transcriptomic and epigenetic multiomes. We catalogue differentially expressed genes and activated epigenetic regions in each cell type during PKD progression, characterizing cell-type-specific responses to Pkd1 deletion. We describe heterogeneous, atypical collecting duct cells as well as proximal tubular cells that constitute cyst epithelia in PKD. The transcriptional regulation of the cyst lining cell marker GPRC5A is conserved between mouse and human PKD cystic epithelia, suggesting shared gene regulatory pathways. Our single nucleus multiomic analysis of mouse PKD provides a foundation to understand the earliest changes molecular deregulation in a mouse model of PKD at a single-cell resolution.

Sequential drug treatment targeting cell cycle and cell fate regulatory programs blocks non-genetic cancer evolution in acute lymphoblastic leukemia

BACKGROUND

Targeted therapies exploiting vulnerabilities of cancer cells hold promise for improving patient outcome and reducing side-effects of chemotherapy. However, efficacy of precision therapies is limited in part because of tumor cell heterogeneity. A better mechanistic understanding of how drug effect is linked to cancer cell state diversity is crucial for identifying effective combination therapies that can prevent disease recurrence.

RESULTS

Here, we characterize the effect of G2/M checkpoint inhibition in acute lymphoblastic leukemia (ALL) and demonstrate that WEE1 targeted therapy impinges on cell fate decision regulatory circuits. We find the highest inhibition of recovery of proliferation in ALL cells with KMT2A-rearrangements. Single-cell RNA-seq and ATAC-seq of RS4;11 cells harboring KMT2A::AFF1, treated with the WEE1 inhibitor AZD1775, reveal diversification of cell states, with a fraction of cells exhibiting strong activation of p53-driven processes linked to apoptosis and senescence, and disruption of a core KMT2A-RUNX1-MYC regulatory network. In this cell state diversification induced by WEE1 inhibition, a subpopulation transitions to a drug tolerant cell state characterized by activation of transcription factors regulating pre-B cell fate, lipid metabolism, and pre-BCR signaling in a reversible manner. Sequential treatment with BCR-signaling inhibitors dasatinib, ibrutinib, or perturbing metabolism by fatostatin or AZD2014 effectively counteracts drug tolerance by inducing cell death and repressing stemness markers.

CONCLUSIONS

Collectively, our findings provide new insights into the tight connectivity of gene regulatory programs associated with cell cycle and cell fate regulation, and a rationale for sequential administration of WEE1 inhibitors with low toxicity inhibitors of pre-BCR signaling or metabolism.

The chromatin landscape of pathogenic transcriptional cell states in rheumatoid arthritis

Synovial tissue inflammation is a hallmark of rheumatoid arthritis (RA). Recent work has identified prominent pathogenic cell states in inflamed RA synovial tissue, such as T peripheral helper cells; however, the epigenetic regulation of these states has yet to be defined. Here, we examine genome-wide open chromatin at single-cell resolution in 30 synovial tissue samples, including 12 samples with transcriptional data in multimodal experiments. We identify 24 chromatin classes and predict their associated transcription factors, including a CD8 + GZMK+ class associated with EOMES and a lining fibroblast class associated with AP-1. By integrating with an RA tissue transcriptional atlas, we propose that these chromatin classes represent 'superstates' corresponding to multiple transcriptional cell states. Finally, we demonstrate the utility of this RA tissue chromatin atlas through the associations between disease phenotypes and chromatin class abundance, as well as the nomination of classes mediating the effects of putatively causal RA genetic variants.

Nucleotide depletion promotes cell fate transitions by inducing DNA replication stress

Control of cellular identity requires coordination of developmental programs with environmental factors such as nutrient availability, suggesting that perturbing metabolism can alter cell state. Here, we find that nucleotide depletion and DNA replication stress drive differentiation in human and murine normal and transformed hematopoietic systems, including patient-derived acute myeloid leukemia (AML) xenografts. These cell state transitions begin during S phase and are independent of ATR/ATM checkpoint signaling, double-stranded DNA break formation, and changes in cell cycle length. In systems where differentiation is blocked by oncogenic transcription factor expression, replication stress activates primed regulatory loci and induces lineage-appropriate maturation genes despite the persistence of progenitor programs. Altering the baseline cell state by manipulating transcription factor expression causes replication stress to induce genes specific for alternative lineages. The ability of replication stress to selectively activate primed maturation programs across different contexts suggests a general mechanism by which changes in metabolism can promote lineage-appropriate cell state transitions.

Building and analyzing metacells in single-cell genomics data

The advent of high-throughput single-cell genomics technologies has fundamentally transformed biological sciences. Currently, millions of cells from complex biological tissues can be phenotypically profiled across multiple modalities. The scaling of computational methods to analyze and visualize such data is a constant challenge, and tools need to be regularly updated, if not redesigned, to cope with ever-growing numbers of cells. Over the last few years, metacells have been introduced to reduce the size and complexity of single-cell genomics data while preserving biologically relevant information and improving interpretability. Here, we review recent studies that capitalize on the concept of metacells-and the many variants in nomenclature that have been used. We further outline how and when metacells should (or should not) be used to analyze single-cell genomics data and what should be considered when analyzing such data at the metacell level. To facilitate the exploration of metacells, we provide a comprehensive tutorial on the construction and analysis of metacells from single-cell RNA-seq data ( https://github.com/GfellerLab/MetacellAnalysisTutorial ) as well as a fully integrated pipeline to rapidly build, visualize and evaluate metacells with different methods ( https://github.com/GfellerLab/MetacellAnalysisToolkit ).

Dependency-aware deep generative models for multitasking analysis of spatial omics data

Spatially resolved transcriptomics (SRT) technologies have significantly advanced biomedical research, but their data analysis remains challenging due to the discrete nature of the data and the high levels of noise, compounded by complex spatial dependencies. Here, we propose spaVAE, a dependency-aware, deep generative spatial variational autoencoder model that probabilistically characterizes count data while capturing spatial correlations. spaVAE introduces a hybrid embedding combining a Gaussian process prior with a Gaussian prior to explicitly capture spatial correlations among spots. It then optimizes the parameters of deep neural networks to approximate the distributions underlying the SRT data. With the approximated distributions, spaVAE can contribute to several analytical tasks that are essential for SRT data analysis, including dimensionality reduction, visualization, clustering, batch integration, denoising, differential expression, spatial interpolation, resolution enhancement and identification of spatially variable genes. Moreover, we have extended spaVAE to spaPeakVAE and spaMultiVAE to characterize spatial ATAC-seq (assay for transposase-accessible chromatin using sequencing) data and spatial multi-omics data, respectively.

Cell-type-specific effects of autism-associated chromosome 15q11.2-13.1 duplications in human brain

Recurrent copy number variation represents one of the most well-established genetic drivers in neurodevelopmental disorders, including autism spectrum disorder (ASD). Duplication of 15q11.2-13.1 (dup15q) is a well-described neurodevelopmental syndrome that increases the risk of ASD by over 40-fold. However, the effects of this duplication on gene expression and chromatin accessibility in specific cell types in the human brain remain unknown. To identify the cell-type-specific transcriptional and epigenetic effects of dup15q in the human frontal cortex we conducted single-nucleus RNA-sequencing and multi-omic sequencing on dup15q cases (n=6) as well as non-dup15q ASD (n=7) and neurotypical controls (n=7). Cell-type-specific differential expression analysis identified significantly regulated genes, critical biological pathways, and differentially accessible genomic regions. Although there was overall increased gene expression across the duplicated genomic region, cellular identity represented an important factor mediating gene expression changes. Neuronal subtypes, showed greater upregulation of gene expression across a critical region within the duplication as compared to other cell types. Genes within the duplicated region that had high baseline expression in control individuals showed only modest changes in dup15q, regardless of cell type. Of note, dup15q and ASD had largely distinct signatures of chromatin accessibility, but shared the majority of transcriptional regulatory motifs, suggesting convergent biological pathways. However, the transcriptional binding factor motifs implicated in each condition implicated distinct biological mechanisms; neuronal JUN/FOS networks in ASD vs. an inflammatory transcriptional network in dup15q microglia. This work provides a cell-type-specific analysis of how dup15q changes gene expression and chromatin accessibility in the human brain and finds evidence of marked cell-type-specific effects of this genetic driver. These findings have implications for guiding therapeutic development in dup15q syndrome, as well as understanding the functional effects CNVs more broadly in neurodevelopmental disorders.

Cell of origin epigenetic priming determines susceptibility to Tet2 mutation

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

Meiotic transcriptional reprogramming mediated by cell-cell communications in humans and mice revealed by scATAC-seq and scRNA-seq

Meiosis is a highly complex process significantly influenced by transcriptional regulation. However, studies on the mechanisms that govern transcriptomic changes during meiosis, especially in prophase I, are limited. Here, we performed single-cell ATAC-seq of human testis tissues and observed reprogramming during the transition from zygotene to pachytene spermatocytes. This event, conserved in mice, involved the deactivation of genes associated with meiosis after reprogramming and the activation of those related to spermatogenesis before their functional onset. Furthermore, we identified 282 transcriptional regulators (TRs) that underwent activation or deactivation subsequent to this process. Evidence suggested that physical contact signals from Sertoli cells may regulate these TRs in spermatocytes, while secreted ENHO signals may alter metabolic patterns in these cells. Our results further indicated that defective transcriptional reprogramming may be associated with non-obstructive azoospermia (NOA). This study revealed the importance of both physical contact and secreted signals between Sertoli cells and germ cells in meiotic progression.

Bidirectional epigenetic editing reveals hierarchies in gene regulation

CRISPR perturbation methods are limited in their ability to study non-coding elements and genetic interactions. In this study, we developed a system for bidirectional epigenetic editing, called CRISPRai, in which we apply activating (CRISPRa) and repressive (CRISPRi) perturbations to two loci simultaneously in the same cell. We developed CRISPRai Perturb-seq by coupling dual perturbation gRNA detection with single-cell RNA sequencing, enabling study of pooled perturbations in a mixed single-cell population. We applied this platform to study the genetic interaction between two hematopoietic lineage transcription factors, SPI1 and GATA1, and discovered novel characteristics of their co-regulation on downstream target genes, including differences in SPI1 and GATA1 occupancy at genes that are regulated through different modes. We also studied the regulatory landscape of IL2 (interleukin-2) in Jurkat T cells, primary T cells and chimeric antigen receptor (CAR) T cells and elucidated mechanisms of enhancer-mediated IL2 gene regulation. CRISPRai facilitates investigation of context-specific genetic interactions, provides new insights into gene regulation and will enable exploration of non-coding disease-associated variants.

CellWalker2: multi-omic discovery of hierarchical cell type relationships and their associations with genomic annotations

CellWalker2 is a graph diffusion-based method for single-cell genomics data integration. It extends the CellWalker model by incorporating hierarchical relationships between cell types, providing estimates of statistical significance, and adding data structures for analyzing multi-omics data so that gene expression and open chromatin can be jointly modeled. Our open-source software enables users to annotate cells using existing ontologies and to probabilistically match cell types between two or more contexts, including across species. CellWalker2 can also map genomic regions to cell ontologies, enabling precise annotation of elements derived from bulk data, such as enhancers, genetic variants, and sequence motifs. Through simulation studies, we show that CellWalker2 performs better than existing methods in cell type annotation and mapping. We then use data from the brain and immune system to demonstrate CellWalker2's ability to discover cell type-specific regulatory programs and both conserved and divergent cell type relationships in complex tissues.

Multimodal analyses reveal genes driving electrophysiological maturation of neurons in the primate prefrontal cortex

The prefrontal cortex (PFC) is critical for myriad high-cognitive functions and is associated with several neuropsychiatric disorders. Here, using Patch-seq and single-nucleus multiomic analyses, we identified genes and regulatory networks governing the maturation of distinct neuronal populations in the PFC of rhesus macaque. We discovered that specific electrophysiological properties exhibited distinct maturational kinetics and identified key genes underlying these properties. We unveiled that RAPGEF4 is important for the maturation of resting membrane potential and inward sodium current in both macaque and human. We demonstrated that knockdown of CHD8, a high-confidence autism risk gene, in human and macaque organotypic slices led to impaired maturation, via downregulation of key genes, including RAPGEF4. Restoring the expression of RAPGEF4 rescued the proper electrophysiological maturation of CHD8-deficient neurons. Our study revealed regulators of neuronal maturation during a critical period of PFC development in primates and implicated such regulators in molecular processes underlying autism.

Causal gene regulatory analysis with RNA velocity reveals an interplay between slow and fast transcription factors

Single-cell expression dynamics, from differentiation trajectories or RNA velocity, have the potential to reveal causal links between transcription factors (TFs) and their target genes in gene regulatory networks (GRNs). However, existing methods either overlook these expression dynamics or necessitate that cells be ordered along a linear pseudotemporal axis, which is incompatible with branching trajectories. We introduce Velorama, an approach to causal GRN inference that represents single-cell differentiation dynamics as a directed acyclic graph of cells, constructed from pseudotime or RNA velocity measurements. Additionally, Velorama enables the estimation of the speed at which TFs influence target genes. Applying Velorama, we uncover evidence that the speed of a TF's interactions is tied to its regulatory function. For human corticogenesis, we find that slow TFs are linked to gliomas, while fast TFs are associated with neuropsychiatric diseases. We expect Velorama to become a critical part of the RNA velocity toolkit for investigating the causal drivers of differentiation and disease.

Type I interferons induce an epigenetically distinct memory B cell subset in chronic viral infection

Memory B cells (MBCs) are key providers of long-lived immunity against infectious disease, yet in chronic viral infection, they do not produce effective protection. How chronic viral infection disrupts MBC development and whether such changes are reversible remain unknown. Through single-cell (sc)ATAC-seq and scRNA-seq during acute versus chronic lymphocytic choriomeningitis viral infection, we identified a memory subset enriched for interferon (IFN)-stimulated genes (ISGs) during chronic infection that was distinct from the T-bet+ subset normally associated with chronic infection. Blockade of IFNAR-1 early in infection transformed the chromatin landscape of chronic MBCs, decreasing accessibility at ISG-inducing transcription factor binding motifs and inducing phenotypic changes in the dominating MBC subset, with a decrease in the ISG subset and an increase in CD11c+CD80+ cells. However, timing was critical, with MBCs resistant to intervention at 4 weeks post-infection. Together, our research identifies a key mechanism to instruct MBC identity during viral infection.

Single-nucleus transcriptional and chromatin accessible profiles reveal critical cell types and molecular architecture underlying chicken sex determination

INTRODUCTION

Understanding the sex determination mechanisms in birds has great significance for the biological sciences and production in the poultry industry. Sex determination in chickens is a complex process that involves fate decisions of supporting cells such as granulosa or Sertoli cells. However, a systematic understanding of the genetic regulation and cell commitment process underlying sex determination in chickens is still lacking.

OBJECTIVES

We aimed to dissect the molecular characteristics associated with sex determination in the gonads of chicken embryos.

METHODS

Single-nucleus RNA-seq (snRNA-seq) and ATAC-seq (snATAC-seq) analysis were conducted on the gonads of female and male chickens at embryonic day 3.5 (E3.5), E4.5, and E5.5.

RESULTS

Here, we provided a time-course transcriptional and chromatin accessible profiling of gonads during chicken sex determination at single-cell resolution. We uncovered differences in cell composition and developmental trajectories between female and male gonads and found that the divergence of transcription and accessibility in gonadal cells first emerged at E5.5. Furthermore, we revealed key cell-type-specific transcription factors (TFs) and regulatory networks that drive lineage commitment. Sex determination signaling pathways, dominated by BMP signaling, are preferentially activated in males during gonadal development. Further pseudotime analysis of the supporting cells indicated that granulosa cells were regulated mainly by the TEAD gene family and that Sertoli cells were driven by the DMRT1 regulons. Cross-species analysis suggested high conservation of both cell types and cell-lineage-specific TFs across the six vertebrates.

CONCLUSIONS

Overall, our study will contribute to accelerating the development of sex manipulation technology in the poultry industry and the application of chickens as a unique model for studying cell fate decisions.

Single-Cell multiomics reveals ENL mutation perturbs kidney developmental trajectory by rewiring gene regulatory landscape

Cell differentiation during organogenesis relies on precise epigenetic and transcriptional control. Disruptions to this regulation can result in developmental abnormalities and malignancies, yet the underlying mechanisms are not well understood. Wilms tumors, a type of embryonal tumor closely linked to disrupted organogenesis, harbor mutations in epigenetic regulators in 30-50% of cases. However, the role of these regulators in kidney development and pathogenesis remains unexplored. By integrating mouse modeling, histological characterizations, and single-cell transcriptomics and chromatin accessibility profiling, we show that a Wilms tumor-associated mutation in the chromatin reader protein ENL disrupts kidney development trajectory by rewiring the gene regulatory landscape. Specifically, the mutant ENL promotes the commitment of nephron progenitors while simultaneously restricting their differentiation by dysregulating key transcription factor regulons, particularly the HOX clusters. It also induces the emergence of abnormal progenitor cells that lose their chromatin identity associated with kidney specification. Furthermore, the mutant ENL might modulate stroma-nephron interactions via paracrine Wnt signaling. These multifaceted effects caused by the mutation result in severe developmental defects in the kidney and early postnatal mortality in mice. Notably, transient inhibition of the histone acetylation binding activity of mutant ENL with a small molecule displaces transcriptional condensates formed by mutant ENL from target genes, abolishes its gene activation function, and restores developmental defects in mice. This work provides new insights into how mutations in epigenetic regulators can alter the gene regulatory landscape to disrupt kidney developmental programs at single-cell resolution in vivo . It also offers a proof-of-concept for the use of epigenetics-targeted agents to rectify developmental defects.

Integrated Single-cell Multiomic Analysis of HIV Latency Reversal Reveals Novel Regulators of Viral Reactivation

Despite the success of antiretroviral therapy, human immunodeficiency virus (HIV) cannot be cured because of a reservoir of latently infected cells that evades therapy. To understand the mechanisms of HIV latency, we employed an integrated single-cell RNA sequencing (scRNA-seq) and single-cell assay for transposase-accessible chromatin with sequencing (scATAC-seq) approach to simultaneously profile the transcriptomic and epigenomic characteristics of ∼ 125,000 latently infected primary CD4+ T cells after reactivation using three different latency reversing agents. Differentially expressed genes and differentially accessible motifs were used to examine transcriptional pathways and transcription factor (TF) activities across the cell population. We identified cellular transcripts and TFs whose expression/activity was correlated with viral reactivation and demonstrated that a machine learning model trained on these data was 75%-79% accurate at predicting viral reactivation. Finally, we validated the role of two candidate HIV-regulating factors, FOXP1 and GATA3, in viral transcription. These data demonstrate the power of integrated multimodal single-cell analysis to uncover novel relationships between host cell factors and HIV latency.

Disruption of maternal vascular remodeling by a fetal endoretrovirus-derived gene in preeclampsia

BACKGROUND

Preeclampsia, one of the most lethal pregnancy-related diseases, is associated with the disruption of uterine spiral artery remodeling during placentation. However, the early molecular events leading to preeclampsia remain unknown.

RESULTS

By analyzing placentas from preeclampsia, non-preeclampsia, and twin pregnancies with selective intrauterine growth restriction, we show that the pathogenesis of preeclampsia is attributed to immature trophoblast and maldeveloped endothelial cells. Delayed epigenetic reprogramming during early extraembryonic tissue development leads to generation of excessive immature trophoblast cells. We find reduction of de novo DNA methylation in these trophoblast cells results in selective overexpression of maternally imprinted genes, including the endoretrovirus-derived gene PEG10 (paternally expressed gene 10). PEG10 forms virus-like particles, which are transferred from the trophoblast to the closely proximate endothelial cells. In normal pregnancy, only a low amount of PEG10 is transferred to maternal cells; however, in preeclampsia, excessive PEG10 disrupts maternal vascular development by inhibiting TGF-beta signaling.

CONCLUSIONS

Our study reveals the intricate epigenetic mechanisms that regulate trans-generational genetic conflict and ultimately ensure proper maternal-fetal interface formation.

Transcriptomic, epigenomic, and spatial metabolomic cell profiling redefines regional human kidney anatomy

A large-scale multimodal atlas that includes major kidney regions is lacking. Here, we employed simultaneous high-throughput single-cell ATAC/RNA sequencing (SHARE-seq) and spatially resolved metabolomics to profile 54 human samples from distinct kidney anatomical regions. We generated transcriptomes of 446,267 cells and chromatin accessibility profiles of 401,875 cells and developed a package to analyze 408,218 spatially resolved metabolomes. We find that the same cell type, including thin limb, thick ascending limb loop of Henle and principal cells, display distinct transcriptomic, chromatin accessibility, and metabolomic signatures, depending on anatomic location. Surveying metabolism-associated gene profiles revealed non-overlapping metabolic signatures between nephron segments and dysregulated lipid metabolism in diseased proximal tubule (PT) cells. Integrating multimodal omics with clinical data identified PLEKHA1 as a disease marker, and its in vitro knockdown increased gene expression in PT differentiation, suggesting possible pathogenic roles. This study highlights previously underrepresented cellular heterogeneity underlying the human kidney anatomy.

A blueprint for tumor-infiltrating B cells across human cancers

B lymphocytes are essential mediators of humoral immunity and play multiple roles in human cancer. To decode the functions of tumor-infiltrating B cells, we generated a B cell blueprint encompassing single-cell transcriptome, B cell-receptor repertoire, and chromatin accessibility data across 20 different cancer types (477 samples, 269 patients). B cells harbored extraordinary heterogeneity and comprised 15 subsets, which could be grouped into two independent developmental paths (extrafollicular versus germinal center). Tumor types grouped into the extrafollicular pathway were linked with worse clinical outcomes and resistance to immunotherapy. The dysfunctional extrafollicular program was associated with glutamine-derived metabolites through epigenetic-metabolic cross-talk, which promoted a T cell-driven immunosuppressive program. These data suggest an intratumor B cell balance between extrafollicular and germinal-center responses and suggest that humoral immunity could possibly be harnessed for B cell-targeting immunotherapy.

Jak2V617F Reversible Activation Shows Its Essential Requirement in Myeloproliferative Neoplasms

Gain-of-function mutations activating JAK/STAT signaling are seen in the majority of patients with myeloproliferative neoplasms (MPN), most commonly JAK2V617F. Although clinically approved JAK inhibitors improve symptoms and outcomes in MPNs, remissions are rare, and mutant allele burden does not substantively change with chronic therapy. We hypothesized this is due to limitations of current JAK inhibitors to potently and specifically abrogate mutant JAK2 signaling. We therefore developed a conditionally inducible mouse model allowing for sequential activation, and then inactivation, of Jak2V617F from its endogenous locus using a combined Dre-rox/Cre-lox dual-recombinase system. Jak2V617F deletion abrogates MPN features, induces depletion of mutant-specific hematopoietic stem/progenitor cells, and extends overall survival to an extent not observed with pharmacologic JAK inhibition, including when cooccurring with somatic Tet2 loss. Our data suggest JAK2V617F represents the best therapeutic target in MPNs and demonstrate the therapeutic relevance of a dual-recombinase system to assess mutant-specific oncogenic dependencies in vivo.

SIGNIFICANCE

Current JAK inhibitors to treat myeloproliferative neoplasms are ineffective at eradicating mutant cells. We developed an endogenously expressed Jak2V617F dual-recombinase knock-in/knock-out model to investigate Jak2V617F oncogenic reversion in vivo. Jak2V617F deletion abrogates MPN features and depletes disease-sustaining MPN stem cells, suggesting improved Jak2V617F targeting offers the potential for greater therapeutic efficacy. See related commentary by Celik and Challen, p. 701. This article is featured in Selected Articles from This Issue, p. 695.

Single-nucleotide variant calling in single-cell sequencing data with Monopogen

Single-cell omics technologies enable molecular characterization of diverse cell types and states, but how the resulting transcriptional and epigenetic profiles depend on the cell's genetic background remains understudied. We describe Monopogen, a computational tool to detect single-nucleotide variants (SNVs) from single-cell sequencing data. Monopogen leverages linkage disequilibrium from external reference panels to identify germline SNVs and detects putative somatic SNVs using allele cosegregating patterns at the cell population level. It can identify 100 K to 3 M germline SNVs achieving a genotyping accuracy of 95%, together with hundreds of putative somatic SNVs. Monopogen-derived genotypes enable global and local ancestry inference and identification of admixed samples. It identifies variants associated with cardiomyocyte metabolic levels and epigenomic programs. It also improves putative somatic SNV detection that enables clonal lineage tracing in primary human clonal hematopoiesis. Monopogen brings together population genetics, cell lineage tracing and single-cell omics to uncover genetic determinants of cellular processes.

FOXO1 enhances CAR T cell stemness, metabolic fitness and efficacy

Chimeric antigen receptor (CAR) T cell therapy has transformed the treatment of haematological malignancies such as acute lymphoblastic leukaemia, B cell lymphoma and multiple myeloma1-4, but the efficacy of CAR T cell therapy in solid tumours has been limited5. This is owing to a number of factors, including the immunosuppressive tumour microenvironment that gives rise to poorly persisting and metabolically dysfunctional T cells. Analysis of anti-CD19 CAR T cells used clinically has shown that positive treatment outcomes are associated with a more 'stem-like' phenotype and increased mitochondrial mass6-8. We therefore sought to identify transcription factors that could enhance CAR T cell fitness and efficacy against solid tumours. Here we show that overexpression of FOXO1 promotes a stem-like phenotype in CAR T cells derived from either healthy human donors or patients, which correlates with improved mitochondrial fitness, persistence and therapeutic efficacy in vivo. This work thus reveals an engineering approach to genetically enforce a favourable metabolic phenotype that has high translational potential to improve the efficacy of CAR T cells against solid tumours.

Identification of lineage-specific epigenetic regulators FOXA1 and GRHL2 through chromatin accessibility profiling in breast cancer cell lines

Breast cancer is a heterogeneous disease, and breast cancer cell lines are invaluable for studying this heterogeneity. However, the epigenetic diversity across these cell lines remains poorly understood. In this study, we performed genome-wide chromatin accessibility analysis on 23 breast cancer cell lines, including 2 estrogen receptor (ER)-positive/human epidermal growth factor receptor 2 (HER2)-negative (ER+/HER2-), 3 ER+/HER2+, 3 HER2+, and 15 triple-negative breast cancer (TNBC) lines. These cell lines were classified into three groups based on their chromatin accessibility: the receptor-positive group (Group-P), TNBC basal group (Group-B), and TNBC mesenchymal group (Group-M). Motif enrichment analysis revealed that only Group-P exhibited coenrichment of forkhead box A1 (FOXA1) and grainyhead-like 2 (GRHL2) motifs, whereas Group-B was characterized by the presence of the GRHL2 motif without FOXA1. Notably, Group-M did not show enrichment of either FOXA1 or GRHL2 motifs. Furthermore, gene ontology analysis suggested that group-specific accessible regions were associated with their unique lineage characteristics. To investigate the epigenetic landscape regulatory roles of FOXA1 and GRHL2, we performed knockdown experiments targeting FOXA1 and GRHL2, followed by assay for transposase-accessible chromatin sequencing analysis. The findings revealed that FOXA1 maintains Group-P-specific regions while suppressing Group-B-specific regions in Group-P cells. In contrast, GRHL2 preserves commonly accessible regions shared between Group-P and Group-B in Group-B cells, suggesting that FOXA1 and GRHL2 play a pivotal role in preserving distinct chromatin accessibility patterns for each group. Specifically, FOXA1 distinguishes between receptor-positive and TNBC cell lines, whereas GRHL2 distinguishes between basal-like and mesenchymal subtypes in TNBC lines.

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.

RAG suppresses group 2 innate lymphoid cells

Antigen specificity is the central trait distinguishing adaptive from innate immune function. Assembly of antigen-specific T cell and B cell receptors occurs through V(D)J recombination mediated by the Recombinase Activating Gene endonucleases RAG1 and RAG2 (collectively called RAG). In the absence of RAG, mature T and B cells do not develop and thus RAG is critically associated with adaptive immune function. In addition to adaptive T helper 2 (Th2) cells, group 2 innate lymphoid cells (ILC2s) contribute to type 2 immune responses by producing cytokines like Interleukin-5 (IL-5) and IL-13. Although it has been reported that RAG expression modulates the function of innate natural killer (NK) cells, whether other innate immune cells such as ILC2s are affected by RAG remains unclear. We find that in RAG-deficient mice, ILC2 populations expand and produce increased IL-5 and IL-13 at steady state and contribute to increased inflammation in atopic dermatitis (AD)-like disease. Further, we show that RAG modulates ILC2 function in a cell-intrinsic manner independent of the absence or presence of adaptive T and B lymphocytes. Lastly, employing multiomic single cell analyses of RAG1 lineage-traced cells, we identify key transcriptional and epigenomic ILC2 functional programs that are suppressed by a history of RAG expression. Collectively, our data reveal a novel role for RAG in modulating innate type 2 immunity through suppression of ILC2s.

Single-cell insights: pioneering an integrated atlas of chromatin accessibility and transcriptomic landscapes in diabetic cardiomyopathy

BACKGROUND

Diabetic cardiomyopathy (DCM) poses a growing health threat, elevating heart failure risk in diabetic individuals. Understanding DCM is crucial, with fibroblasts and endothelial cells playing pivotal roles in driving myocardial fibrosis and contributing to cardiac dysfunction. Advances in Multimodal single-cell profiling, such as scRNA-seq and scATAC-seq, provide deeper insights into DCM's unique cell states and molecular landscape for targeted therapeutic interventions.

METHODS

Single-cell RNA and ATAC data from 10x Multiome libraries were processed using Cell Ranger ARC v2.0.1. Gene expression and ATAC data underwent Seurat and Signac filtration. Differential gene expression and accessible chromatin regions were identified. Transcription factor activity was estimated with chromVAR, and Cis-coaccessibility networks were calculated using Cicero. Coaccessibility connections were compared to the GeneHancer database. Gene Ontology analysis, biological process scoring, cell-cell communication analysis, and gene-motif correlation was performed to reveal intricate molecular changes. Immunofluorescent staining utilized various antibodies on paraffin-embedded tissues to verify the findings.

RESULTS

This study integrated scRNA-seq and scATAC-seq data obtained from hearts of WT and DCM mice, elucidating molecular changes at the single-cell level throughout the diabetic cardiomyopathy progression. Robust and accurate clustering analysis of the integrated data revealed altered cell proportions, showcasing decreased endothelial cells and macrophages, coupled with increased fibroblasts and myocardial cells in the DCM group, indicating enhanced fibrosis and endothelial damage. Chromatin accessibility analysis unveiled unique patterns in cell types, with heightened transcriptional activity in myocardial cells. Subpopulation analysis highlighted distinct changes in cardiomyocytes and fibroblasts, emphasizing pathways related to fatty acid metabolism and cardiac contraction. Fibroblast-centered communication analysis identified interactions with endothelial cells, implicating VEGF receptors. Endothelial cell subpopulations exhibited altered gene expressions, emphasizing contraction and growth-related pathways. Candidate regulators, including Tcf21, Arnt, Stat5a, and Stat5b, were identified, suggesting their pivotal roles in DCM development. Immunofluorescence staining validated marker genes of cell subpopulations, confirming PDK4, PPARγ and Tpm1 as markers for metabolic pattern-altered cardiomyocytes, activated fibroblasts and endothelial cells with compromised proliferation.

CONCLUSION

Our integrated scRNA-seq and scATAC-seq analysis unveils intricate cell states and molecular alterations in diabetic cardiomyopathy. Identified cell type-specific changes, transcription factors, and marker genes offer valuable insights. The study sheds light on potential therapeutic targets for DCM.

Decoding the universal human chromatin landscape through teratoma-based profiling

Teratoma formation is key for evaluating differentiation of human pluripotent stem cells into embryonic germ layers and serves as a model for understanding stem cell differentiation and developmental processes. Its potential for insights into epigenome and transcriptome profiling is significant. This study integrates the analysis of the epigenome and transcriptome of hESC-generated teratomas, comparing transcriptomes between hESCs and teratomas. It employs cell type-specific expression patterns from single-cell data to deconvolve RNA-Seq data and identify cell types within teratomas. Our results provide a catalog of activating and repressive histone modifications, while also elucidating distinctive features of chromatin states. Construction of an epigenetic signature matrix enabled the quantification of diverse cell populations in teratomas and enhanced the ability to unravel the epigenetic landscape in heterogeneous tissue contexts. This study also includes a single cell multiome atlas of expression (scRNA-Seq) and chromatin accessibility (scATAC-Seq) of human teratomas, further revealing the complexity of these tissues. A histology-based digital staining tool further complemented the annotation of cell types in teratomas, enhancing our understanding of their cellular composition. This research is a valuable resource for examining teratoma epigenomic and transcriptomic landscapes and serves as a model for epigenetic data comparison.

Temporal coordination of the transcription factor response to H2O2 stress

Oxidative stress from excess H2O2 activates transcription factors that restore redox balance and repair oxidative damage. Although many transcription factors are activated by H2O2, it is unclear whether they are activated at the same H2O2 concentration, or time. Dose-dependent activation is likely as oxidative stress is not a singular state and exhibits dose-dependent outcomes including cell-cycle arrest and cell death. Here, we show that transcription factor activation is both dose-dependent and coordinated over time. Low levels of H2O2 activate p53, NRF2 and JUN. Yet under high H2O2, these transcription factors are repressed, and FOXO1, NF-κB, and NFAT1 are activated. Time-lapse imaging revealed that the order in which these two groups of transcription factors are activated depends on whether H2O2 is administered acutely by bolus addition, or continuously through the glucose oxidase enzyme. Finally, we provide evidence that 2-Cys peroxiredoxins control which group of transcription factors are activated.

A germline point mutation in the MYC-FBW7 phosphodegron initiates hematopoietic malignancies

Oncogenic activation of MYC in cancers predominantly involves increased transcription rather than coding region mutations. However, MYC-dependent lymphomas frequently acquire point mutations in the MYC phosphodegron, including at threonine 58 (T58), where phosphorylation permits binding via the FBW7 ubiquitin ligase triggering MYC degradation. To understand how T58 phosphorylation functions in normal cell physiology, we introduced an alanine mutation at T58 (T58A) into the endogenous c-Myc locus in the mouse germline. While MYC-T58A mice develop normally, lymphomas and myeloid leukemias emerge in ∼60% of adult homozygous T58A mice. We found that primitive hematopoietic progenitor cells from MYC-T58A mice exhibit aberrant self-renewal normally associated with hematopoietic stem cells (HSCs) and up-regulate a subset of MYC target genes important in maintaining stem/progenitor cell balance. In lymphocytes, genomic occupancy by MYC-T58A was increased at all promoters compared with WT MYC, while genes differentially expressed in a T58A-dependent manner were significantly more proximal to MYC-bound enhancers. MYC-T58A lymphocyte progenitors exhibited metabolic alterations and decreased activation of inflammatory and apoptotic pathways. Our data demonstrate that a single point mutation stabilizing MYC is sufficient to skew target gene expression, producing a profound gain of function in multipotential hematopoietic progenitors associated with self-renewal and initiation of lymphomas and leukemias.

The chromatin accessibility dynamics during cell fate specifications in zebrafish early embryogenesis

Chromatin accessibility plays a critical role in the regulation of cell fate decisions. Although gene expression changes have been extensively profiled at the single-cell level during early embryogenesis, the dynamics of chromatin accessibility at cis-regulatory elements remain poorly studied. Here, we used a plate-based single-cell ATAC-seq method to profile the chromatin accessibility dynamics of over 10 000 nuclei from zebrafish embryos. We investigated several important time points immediately after zygotic genome activation (ZGA), covering key developmental stages up to dome. The results revealed key chromatin signatures in the first cell fate specifications when cells start to differentiate into enveloping layer (EVL) and yolk syncytial layer (YSL) cells. Finally, we uncovered many potential cell-type specific enhancers and transcription factor motifs that are important for the cell fate specifications.

SMARCA4 is a haploinsufficient B cell lymphoma tumor suppressor that fine-tunes centrocyte cell fate decisions

SMARCA4 encodes one of two mutually exclusive ATPase subunits in the BRG/BRM associated factor (BAF) complex that is recruited by transcription factors (TFs) to drive chromatin accessibility and transcriptional activation. SMARCA4 is among the most recurrently mutated genes in human cancer, including ∼30% of germinal center (GC)-derived Burkitt lymphomas. In mice, GC-specific Smarca4 haploinsufficiency cooperated with MYC over-expression to drive lymphomagenesis. Furthermore, monoallelic Smarca4 deletion drove GC hyperplasia with centroblast polarization via significantly increased rates of centrocyte recycling to the dark zone. Mechanistically, Smarca4 loss reduced the activity of TFs that are activated in centrocytes to drive GC-exit, including SPI1 (PU.1), IRF family, and NF-κB. Loss of activity for these factors phenocopied aberrant BCL6 activity within murine centrocytes and human Burkitt lymphoma cells. SMARCA4 therefore facilitates chromatin accessibility for TFs that shape centrocyte trajectories, and loss of fine-control of these programs biases toward centroblast cell-fate, GC hyperplasia and lymphoma.

Single-cell 'omic profiles of human aortic endothelial cells in vitro and human atherosclerotic lesions ex vivo reveal heterogeneity of endothelial subtype and response to activating perturbations

Heterogeneity in endothelial cell (EC) sub-phenotypes is becoming increasingly appreciated in atherosclerosis progression. Still, studies quantifying EC heterogeneity across whole transcriptomes and epigenomes in both in vitro and in vivo models are lacking. Multiomic profiling concurrently measuring transcriptomes and accessible chromatin in the same single cells was performed on six distinct primary cultures of human aortic ECs (HAECs) exposed to activating environments characteristic of the atherosclerotic microenvironment in vitro. Meta-analysis of single-cell transcriptomes across 17 human ex vivo arterial specimens was performed and two computational approaches quantitatively evaluated the similarity in molecular profiles between heterogeneous in vitro and ex vivo cell profiles. HAEC cultures were reproducibly populated by four major clusters with distinct pathway enrichment profiles and modest heterogeneous responses: EC1-angiogenic, EC2-proliferative, EC3-activated/mesenchymal-like, and EC4-mesenchymal. Quantitative comparisons between in vitro and ex vivo transcriptomes confirmed EC1 and EC2 as most canonically EC-like, and EC4 as most mesenchymal with minimal effects elicited by siERG and IL1B. Lastly, accessible chromatin regions unique to EC2 and EC4 were most enriched for coronary artery disease (CAD)-associated single-nucleotide polymorphisms from Genome Wide Association Studies (GWAS), suggesting that these cell phenotypes harbor CAD-modulating mechanisms. Primary EC cultures contain markedly heterogeneous cell subtypes defined by their molecular profiles. Surprisingly, the perturbations used here only modestly shifted cells between subpopulations, suggesting relatively stable molecular phenotypes in culture. Identifying consistently heterogeneous EC subpopulations between in vitro and ex vivo models should pave the way for improving in vitro systems while enabling the mechanisms governing heterogeneous cell state decisions.

Single-cell epigenomic dysregulation of Systemic Sclerosis fibroblasts via CREB1/EGR1 axis in self-assembled human skin equivalents

Systemic sclerosis (SSc) is an autoimmune disease characterized by skin fibrosis, internal organ involvement and vascular dropout. We previously developed and phenotypically characterized an in vitro 3D skin-like tissue model of SSc, and now analyze the transcriptomic (scRNA-seq) and epigenetic (scATAC-seq) characteristics of this model at single-cell resolution. SSc 3D skin-like tissues were fabricated using autologous fibroblasts, macrophages, and plasma from SSc patients or healthy control (HC) donors. SSc tissues displayed increased dermal thickness and contractility, as well as increased α-SMA staining. Single-cell transcriptomic and epigenomic analyses identified keratinocytes, macrophages, and five populations of fibroblasts (labeled FB1 - 5). Notably, FB1 APOE-expressing fibroblasts were 12-fold enriched in SSc tissues and were characterized by high EGR1 motif accessibility. Pseudotime analysis suggests that FB1 fibroblasts differentiate from a TGF-β1-responsive fibroblast population and ligand-receptor analysis indicates that the FB1 fibroblasts are active in macrophage crosstalk via soluble ligands including FGF2 and APP. These findings provide characterization of the 3D skin-like model at single cell resolution and establish that it recapitulates subsets of fibroblasts and macrophage phenotypes observed in skin biopsies.

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.